2. 国土资源部煤炭资源勘查与综合利用重点实验室, 西安 710054
2. Key Laboratory of Coal Resources Exploration and Comprehensive Utilization, Ministry of Land and Resources, Xi'an 710054, China
随着我国浅部资源的消耗殆尽,国家提出加强“第二深度空间”探矿以及“攻深探盲”的资源战略(滕吉文,2006).为实现地下500~2000 m深度范围内矿产资源的有效探测,需在已有基础之上完善或创新更多行之有效的地球物理勘探方法.瞬变电磁法(transient electromagnetic method,TEM)是一种建立在电磁感应原理基础上的时间域人工源电磁探测方法.它是利用不接地回线(磁性源)或接地线源(电性源)向地下发送一次脉冲场,在一次脉冲场的间歇期间,利用线圈或接地电极观测二次涡流场的方法.TEM可以在距离发射源很近的区域实现大深度的探测,并可以实现频率域电磁法无法实 现的同点观测(Nabighian,1979;Kaufman and Keller,1983; 朴化荣,1990;牛之琏,1992).时间域电磁法是近年来发展很快的电法勘探分支方法,在国际上有人称作是电法的“二次革命”.目前TEM已广泛应用于金属矿勘探、构造填图、油气田、煤田、地下水、冻土带、海洋地质、水文工程地质及工程检测等领域(李貅,2002;嵇艳菊等,2005;薛国强等,2007;陈卫营,2013; 韩自强等,2015).
传统的TEM工作多采用磁性源装置形式或电性源长偏移距形式(LOTEM).磁性源TEM体积效应小、探测精度高,但是探测深度较小(一般不超过500 m)并且对高阻目标体不敏感;LOTEM探测深度大(可达10 km),但是精度一般不高.为此,最近人们对采用电性源激发并在小偏移距范围内观测纯二次场信号的短偏移TEM装置越来越感兴趣,薛国强等(2013)将该种装置命名为SOTEM.大量的理论与实践证明,SOTEM装置相较于磁性源瞬变 电磁法,大大提升了探测深度,相较于LOTEM提高了探测精度(Ziolkowski,2010;Cuevas and Alumbaugh,2011;薛国强等,2013;Xue et al.,2014; Chen et al.,2015;Zhou et al.,2015;王显祥等,2015).
但是,目前对该方法的研究和应用尚不普遍也不完善,主要原因是相较于磁性源TEM,SOTEM的电磁场理论更为复杂,并且缺乏配套的数据处理解释软件,涉及二、三维的研究和应用几乎没有进行.而在一维理论下,详细研究电性源电磁场的各种特性是深入理解方法内涵、推动方法发展的有效途径.本文基于电性源瞬变电磁一维正演理论,研究了SOTEM地下感应电流扩散、电磁场响应平面分布、多偏移距衰减等特性,以及SOTEM的最佳观测区 域问题.研究结果对进一步研究和推广应用SOTEM具有指导意义.
2 SOTEM方法介绍电性源短偏移距瞬变电磁法(Short-offset transient electromagnetic method,简称SOTEM)是在传统LOTEM法基础上提出的一种新型瞬变电磁工作装置.它利用长约500~2000 m的接地长导线为发射源,供以强度一般为10~40 A的双极性矩形阶跃电流,并在小于2倍探测深度的偏移距范围内观测瞬变电磁场(图 1).与LOTEM采用连续波形激励、在大偏移距处(一般3~8倍探测深度)观测总场响应不同,SOTEM在小偏移距范围内观测纯二次场响应.这种工作方式一方面提高了观测信号的信噪比,另一方面减小了体积效应的影响,从而大大降低了数据处理的难度并提高了处理结果的准确度(薛国强等,2013;陈卫营和薛国强,2013).实际工作中,一般观测垂直磁场分量随时间的导数(感应电压)和水平电场分量.目前SOTEM的主要数据处理流程包括预处理、全期视电阻率计算和一维反演三个步骤.笔者将SOTEM应用于河南、山东、安徽等地的 盐腔、煤田、金属矿等领域的实际勘探证明,SOTEM 能够很好地实现地下2000 m深度以内的精细探测(薛国强等,2013; Xue et al.,2014; Chen et al.,2015; Li,2015).
目前,SOTEM的研究与应用主要基于1D正反演理论.纳比吉安(1992)给出了电偶极子源在层状1D大地模型表面产生的瞬变响应的计算方法.鉴于SOTEM偏移距较小,偶极子假设条件(一般认为当偏移距大于5倍的发射源长度时)不成立,因此在计算时需先对发射源进行偶极子分割,分段计算响应再进行积分求和得到整个源的响应(周楠楠等,2012).下文中对SOTEM响应特性与观测区域特性的研究便是基于上述1D正演理论.
3 SOTEM地下感应电流扩散特性地面观测到瞬变电磁场的响应是由大地感应涡 流产生的,根据Lewis和Lee(1978)、Nabighian(1979)、Reid和Macnae(1998)、Wang(2002)、闫述等(2002)的研究,我们知道大地中的感应电流可近似地用圆形电流环来等效,即等效电流可以形象地表示为由发射源吹出的“烟圈”.“烟圈”随时间的扩散过程可以用地下电场等值线来描述.在电性源瞬变电磁中,由于接地电极的存在,大地成为发射源的一部分,与长导线共同组成回路.相较于水平或垂直不接地回线源中仅存在沿发射源流动的电流,长导线-大地系统在地下电流方向要复杂的多,不仅存在水平感应电流还存在垂直感应电流,这也就导致了对应电磁场扩散及分布的复杂性(Key,2009).
为了研究这两个方向的感应电流的扩散特性,计算了长度为200 m的电性源(发射电流10 A)在均匀半空间(电阻率为100 Ωm)产生的水平感应电流和垂直感应电流.图 2为两个不同时刻(1 ms和3 ms)的水平感应电流的地下扩散图.电性源的水平感应电流由两部分组成:上部水平感应电流和下部水平感应电流(图 1a中红色虚线圈定区域)(Um et al.,2012).上部水平电流最大值在全期范围内都集中在场源附近.下部感应水平电流又称作“返回电流”(Gunderson et al.,1986),向下扩散速度较快但是振幅要比上部电流小.由于水平电流的 最大值一直离场源较近,所以当在场源附近存在不均匀体时电性源的电磁响应会产生静态偏移效应(Newman,1989).
图 3为电性源在两个不同时刻(1 ms和3 ms)激发的垂直感应电流的地下扩散图.垂直感应电流表现出了许多重要的特性.首先,垂直感应电流极大值集中于两个接地电极的下方,并逐渐向外、向下扩散;而在发射源中点的正下方不存在垂直感应电流,这是由于垂直电场分量仅由相互对称的两个接地项激励,因此在对称点上正负抵消为零(纳比吉安,1992).其次,垂直感应电流的强度随时间的衰减速度较水平感应电流更快,3 ms时的垂直电流强度比1 ms时衰减了近10倍,而水平感应电流仅衰减了7倍左右.另外,通过对比两种感应电流的振幅强度,发现在给定时间,垂直感应电流极大值的振幅差不多仅为水平感应电流的20%.也就是说垂直感应电流较水平感应电流具有振幅小、衰减快的特性.
图 4为回线源产生的水平感应电流和电性源产生的水平、垂直感应电流在不同时刻的极大值扩散路径图.Nabighian(1979)和牛之琏(1992)对于回线源的感应电流都做了比较详细的研究,指出感应电流的极大值大致沿与地面呈30°角的方向扩散.而图 4所示的电性源垂直感应电流极大值大致沿与地面呈45°角的方向扩散,水平感应电流极大值则是垂直向下扩散.图 4更加形象地显示出在三种感应电流中,电性源的垂直感应电流扩散速度最快,其次为回线源水平感应电流,扩散最慢的为电性源水平感应电流.
目前在陆地环境的SOTEM工作中,直接测量垂直电场分量具有一定困难,但是,Kaufman和Keller(1983)指出垂直感应电流在高阻薄层界面产生的积累电荷对于高阻薄层的探测极为有利.正是基于电性源这种在高阻体探测方面的优越性,本方法在油气田、采空区等高阻目标体勘探领域发挥着不可替代的作用.
3 SOTEM瞬变电磁场平面分布特性地下感应电流的扩散特性决定了地面电磁场的扩散与分布,研究电性源瞬变电磁场随时间在地面的扩散、分布规律对分析电磁场特性、选择合适的电磁分量进行观测具有重要的意义.置于地表的电偶极源在直角坐标系下可以产生全部六个电磁场分量,它们分别是Ex、Ey、Ez和Hx、Hy、Hz,虽然这六个场量都包含了地下介质的电性结构信息,对地层都具有探测能力,但在实际的地面工作中,Ez分量不易观测,因此在下述计算、分析中,暂不考虑Ez分量.计算了电阻率为100 Ωm的均匀半空间(由于Ey在均匀半空间表面产生的响应为零,因此计算Ey时所用模型为含一低阻薄层(埋深500 m,厚度50 m,电阻率10 Ωm)的均匀半空间)情况下Ex、Ey、Hx、Hy、Hz五个分量在不同时刻、不同位置的响应值,并绘制了两个时刻的响应平面分布图,上述计算中发射源参数为长度1000 m,电流强度10 A.
平面分布图(图 5)更清晰地显示出Ex的扩散与分布特性:在平面各个方向分布均匀,全期探测范围内电场极大值集中于发射源附近,离发射源越远信号强度就越低,这与前面分析的上部水平感应电流集中于发射源附近相吻合.
图 6所示的平面分布图表现出,Ey场值的分布极不均匀,极值集中于四个象限内.根据纳比吉安(1992)可知水平电场Ey的表达式中仅包含接地项,因此根据电磁场矢量叠加原理,两个极性相反的电极导致场在赤道向Φ=90°以及轴向Φ=0°区域等于零,在Φ=45°时场值最大.图 7所示的水平磁场Hx分布特征与Ey的分布特征基本上类似.
从平面分布图 8也可以看出,正向的水平磁场集中于发射源附近,这是因为Hy是由发射源中与电流方向一致的正向感应电流激发而来;而负向水平磁场则集中于两侧离发射源一定距离的范围,这是因为此处的Hy是由负向的“返回电流”感应得到的.Hy场随时间扩散,但是极值一直集中在线源处,这是因为产生水平磁场的垂向感应电流的极值是沿发射源所在平面垂直向下扩散的.
图 9说明Hz的极值分布于发射源两侧,并随时间推移沿发射源中垂线逐渐向远处移动.分析对比上述五个分量,可以发现,Ex和Hz这两个分量在平面内分布较为均匀、简单,特别是在沿发射源导 线和其中垂线的方向上,这对保证同一测线上各测点所测响应幅值趋于均匀非常重要.而Ey和Hx由于响应仅由接地项激发,致使场分布呈明显的四象限分布,在与发射源平行的方向上极不均匀.Hy则由于返回电流的作用,在全期内出现方向相反的两种场,使测量响应存在变号现象.所以,为了使观测得到的电磁响应简单、易于处理,我们要尽量选择场分布均匀、简单的分量.
上面对电磁场各分量在固定时刻的平面分布特性进行了分析,而指定位置处在整个观测时间范围内的衰减曲线是实际观测时主要分析的对象.为了说明不同偏移距处各分量的响应衰减特性,计算了与上节中所用模型和发射源相同时(对于Ey同样采用含中间低阻薄层的半空间模型)四个偏移距处(分别为1000 m,2000m,3000 m和5000 m)的响应衰减曲线,如图 10所示.
从图 10a可以看出,水平电场Ex随偏移距增大,信号强度在早期时间段内急剧减小,而在晚期时间段内不同偏移距处的响应强度趋于一致.Ey的计算结果如图 10b所示,可以看出,随着偏移距的增大,Ey早期响应的强度明显降低,而晚期道的信号强度则相差不大,这说明Ey场值的极值在全期内都集中于发射源附近,并不随时间推移向外扩散,这与Ex表现出类似的性质.图 10b中曲线最明显的特征是在1~10 ms之间每条衰减曲线都出现一个突然下降的极小值点,而其他地质模型(含高阻薄层的半空间)则不会出现该现象,因此可以推断,模型中的低阻薄层是造成上述信号急速衰减的原因.Ey这种对低阻薄层反应灵敏、表象剧烈的特性可被利用.水平磁场Hx也仅由接地项激发,因此场的衰减与Ey非常类似,但是从图 10c可以看出,随着偏移距增大,Hx早期信号强度减小而晚期信号强度增强,这说明Hx的极值是随时间推移逐渐向外扩散的,这点与Ey是不同的.水平磁场Hy的计算结果表明,在计算时间范围内Hy存在正、负两个方向的场量,早期道场值为负,晚期道场值为正(图 10d).垂直磁场Hz的信号强度在早期随偏移距增大而减小,在晚期随偏移距增大而增大(图 10e).
需要说明的是,虽然目前瞬变电磁工作中,还是以观测磁场随时间的导数(感应电压)为主,但是随着高温超导技术的发展,直接观测磁场强度已成为现实,并显示出许多优点(戴远东和王世光,1996;陈晓东等,2002).作者通过对Hz与∂B/∂t之间的对比发现,两种物理量在分布与衰减方面表现出相似的特性,因此本文对磁场分量的讨论仅以磁场强度为例.
5 SOTEM最佳观测区域前面已经分析了SOTEM地下感应电流的扩散特性以及目前主要观测的Ex和Hz(∂B/∂t)分量的信号衰减和分布特性,指出需在合适的区域观测相应的电磁场分量,以确保信号的强度和分布适于观测和处理.而在实际工作中,针对不同的探测目标体,不仅要考虑观测合适的电磁场分量,还需选择合适的观测区域,以达到最佳的探测效果.为了定量分析Ex和Hz对不同地层敏感区域的分布特性,我们以如图 11所示的中间薄层模型为例,计算了薄层相对围岩分别为低阻和高阻时的电磁场响应,并利用(1)式定义的相对异常进行定量对比.计算中发射源参数为发射源长1000 m,发送电流10 A.
式中Fia代表存在薄层时模型在第i时刻产生电磁场响应,Fi0代表不存在薄层的均匀半空间在第i时刻产生的电磁场响应.Pi值越大表示薄层产生的相对异常越大,代表信号对薄层越敏感.
图 12为以每个观测点处的最大Pi值绘制而成的平面分布图.分析图 12可以得出如下的结论:Hz对低阻薄层的敏感程度远大于对高阻薄层的,其敏感区域集中于发射源附近;Ex对低阻、高阻薄层的 敏感程度相当,但是对低阻薄层的敏感区域主要集 中于赤道区域,而对高阻薄层的敏感区域主要集中于轴向区域,并且这两个敏感区域都处于离发射源一定距离的位置.该距离与目标体的深度与围岩的电阻率有关,因此在利用Ex分量进行探测时,应先根据已知信息估算围岩电阻率和目标体埋深,从而选择合适的观测区域.
接地导线源在地下既能激发水平向感应电流也能激发垂直向感应电流,水平感应电流又分为上部 水平感应电流和下部水平感应电流(又称作返回电 流),水平感应电流的极大值主要集中于发射源附近
并垂直向下扩散;垂直感应电流极大值沿与地面呈45°角的方向向下、向外扩散,并且具有较低的振幅和较快的扩散速度.电性源激发的六个方向的电磁场分量都具有对地的探测能力,但是考虑到地面观测的方便性和各分量的传播、分布特点,大多数情况仅利用垂直磁场分量Hz(∂B/∂t)和水平电场分量Ex.其中,Hz仅对低阻目标体敏感,且敏感区域集 中于赤道向区域;Ex既对低阻体敏感也对高阻体敏感,对低阻体的敏感区域集中于赤道向区域,而对高阻体的敏感区域集中于轴向区域,实际工作中需要根据目标体深度及围岩电性合理选择观测区域.
目前,随着深部探测需求的增加,与“电性源”及“短偏移距”相关的研究工作也越来越多,表明该工作装置具有很强的探测能力有待开发.本文仅对电性源瞬变场响应的基本特性做出了研究和分析,旨在实现对SOTEM方法的深入理解,为后续工作建立理论基础.今后的研究应着重于SOTEM响应精确计算、二、三维正演、反演成像技术、多分量综合应用、附加效应消除等领域.
[1] | Chen W Y, Xue G Q. 2013. Application on coal-mine voids detection with multidevice TEM technology. Progress in Geophys. (in Chinese), 28(5):2709-2717, doi:10.6038/pg20130554. |
[2] | Chen W Y, Xue G Q. 2014. Effective skin depth of whole EM field due to a grounded wire source. Chinese Journal of Geophysics (in Chinese), 57(7):2314-2320, doi:10.6038/cjg20140725. |
[3] | Chen W Y, Xue G Q, Muhammad Y K, et al. 2015. Application of short-offset TEM (SOTEM) technique in mapping water-enriched zones of coal stratum, an example from East China. Pure and Applied Geophysics, 172(6):1643-1651, doi:10.1007/s00024-014-1028-z. |
[4] | Chen X D, Zhao Y, Wang C J, et al. 2002. The development of HTc RF SQUID magnetometer and its field test in TEM. Acta Geoscientia Sinica (in Chinese), 23(2):179-182. |
[5] | Cuevas N H, Alumbaugh D. 2011. Near-source response of a resistive layer to a vertical or horizontal electric dipole excitation. Geophysics, 76(6):F353-F371. |
[6] | Dai Y D, Wang S G. 1996. Comparison between high Tc SQUID sensor and conventional inductive sensor in transient magnetoelectric method for magnetotelluric sounding. Chinese Journal of Low Temperature Physics (in Chinese), 18(1):11-19. |
[7] | Gunderson B M, Newman G A, Hohmann G W. 1986. Three-dimensional transient electromagnetic responses for a grounded source. Geophysics, 51(11):2117-2130. |
[8] | Han Z Q, Liu T, Ouyang J, et al. 2015. Application of large rectangular loop source TEM all-time apparent resistivity on colliery goaf exploration. Progress in Geophysics (in Chinese), 30(1):343-349, doi:10.6038/pg20150150. |
[9] | Ji Y J, Lin J, Zhu K G, et al. 2005. Underground water prospecting by transient electromagnetic method. Progress in Geophysics (in Chinese), 20(3):828-833. |
[10] | Kaufman A A, Keller G V. 1983. Frequency and Transient Soundings. New York:Elsevier. |
[11] | Key K. 2009. 1D inversion of multicomponent, multifrequency marine CSEM data:Methodology and synthetic studies for resolving thin resistive layers. Geophysics, 74(2):F9-F20. |
[12] | Lewis R, Lee T. 1978. The transient electric fields about a loop on a halfspace. Bull. Austral. Soc. Expl. Geophys., 9(4):173-177. |
[13] | Li X. 2002. Theory and Application of Transient Electromagnetic Method (in Chinese). Xi'an:Shanxi Science and Technology Press. |
[14] | Li H, Xue G Q, Zhou N N, et al. 2015. Appraisal of an array TEM method in detecting a mined-out area beneath a conductive layer. Pure and Applied Geophysics, 172(10):2917-2929, doi:10.1007/s00024-015-1075-0. |
[15] | Nabighian M N. 1979. Quasi-static transient response of a conducting half-space-An approximate representation. Geophysics, 44(10):1700-1705. |
[16] | Nabighian M N, translated by Zhao J X, Wang Y J. 1992. Electromagnetic Methods in Applied Geophysics Volume 1 Theory (in Chinese). Beijing:Geological Press. |
[17] | Newman G A. 1989. Deep transient electromagnetic soundings with a grounded source over near-surface conductors. Geophysical Journal International, 98(3):587-601. |
[18] | Niu Z L. 1992. Principle of Time-domain Electromagnetic Method (in Chinese). Changsha:Central South University of Technology Press. |
[19] | Piao H R. 1990. Theory of Electromagnetic Sounding Method (in Chinese). Beijing:Geological Publishing House. |
[20] | Reid J E, Macnae J C. 1998. Comments on the electromagnetic "smoke ring" concept. Geophysics, 63(6):1908-1913. |
[21] | Streich R, Becken M. 2011. Electromagnetic fields generated by finite-length wire sources:comparison with point dipole solutions. Geophysical Prospecting, 59(2):361-374, doi:10.1111/j.1365-2478.2010.00926.x. |
[22] | Teng J W. 2006. Strengthening geophysical exploration and exploitation of metallic minerals in the second deep space of the crustal interior. Geological Bulletin of China (in Chinese), 25(7):767-771. |
[23] | Um E S, Alumbaugh D A, Harris J M, et al. 2012. Numerical modeling analysis of short-offset electric-field measurements with a vertical electric dipole source in complex offshore environments. Geophysics, 77(5):E329-E341. |
[24] | Wang T. 2002. The electromagnetic smoke ring in a transversely isotropic medium. Geophysics, 67(6):1779-1789. |
[25] | Wang X X, Di Q Y, Tang J. 2015. The preliminary study of the ground source transient. Progress in Geophysics (in Chinese), 30(2):872-877, doi:10.6038/pg20150253. |
[26] | Xue G Q, Li X, Di Q Y. 2007. The progress of TEM in theory and application. Progress in Geophysics (in Chinese), 22(4):1195-1200. |
[27] | Xue G Q, Chen W Y, Zhou N N, et al. 2013. Short-offset TEM technique with a grounded wire source for deep sounding. Chinese Journal of Geophysics (in Chinese), 56(1):255-261, doi:10.6038/cjg20130126. |
[28] | Xue G Q, Gelius L J, Sakyi P A, et al. 2014. Discovery of a hidden BIF deposit in Anhui province, China by integrated geological and geophysical investigations. Ore Geology Reviews, 63:470-477. |
[29] | Yan S, Chen M S, Fu J M. 2002. Direct time-domain numerical analysis of transient electromagnetic fields. Chinese J. Geophys. (in Chinese), 45(2):275-284. |
[30] | Zhou N N, Xue G Q, Kong X R. 2012. Response erro of electric-source transient electromagnetic method. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese), 34(6):640-645. |
[31] | Zhou N N, Xue G Q, Gelius L J, et al. 2015. Analysis of the near-source error in TEM due to the dipole hypothesis. Journal of Applied Geophysics, 116:75-83. |
[32] | Ziolkowski A. 2010. Short-offset transient electromagnetic geophysical survey. US Patent Application Publication, 2010/0201367. |
[33] | 陈卫营, 薛国强. 2013. 瞬变电磁法多装置探测技术在煤矿采空区调查中的应用. 地球物理学进展, 28(5):2709-2717, doi:10.6038/pg20130554. |
[34] | 陈卫营, 薛国强. 2014. 接地导线源电磁场全域有效趋肤深度. 地球物理学报, 57(7):2314-2320, doi:10.6038/cjg20140725. |
[35] | 陈晓东, 赵毅, 王赤军等. 2002. 高温超导磁强计的研制及在TEM上的野外试验. 地球学报, 23(2):179-182. |
[36] | 戴远东, 王世光. 1996. 在瞬变电磁法大地电磁测量中高Tc SQUID探头与常规探头的比较. 低温物理学报, 18(1):11-19. |
[37] | 韩自强, 刘涛, 欧阳进等. 2015. 矩形大定源回线TEM法全区视电阻率在煤田采空区勘探中的应用. 地球物理学进展, 30(1):343-349, doi:10.6038/pg20150150. |
[38] | 嵇艳菊, 林君, 朱凯光等. 2005. 利用瞬变电磁技术进行地下水资源勘察. 地球物理学进展, 20(3):828-833. |
[39] | 李貅. 2002. 瞬变电磁测深的理论与应用. 西安:陕西科学技术出版社. |
[40] | 米萨克 N 纳比吉安. 赵经祥, 王艳君译. 1992. 勘查地球物理 电磁法 第一卷 理论. 北京:地质出版社. |
[41] | 牛之琏. 1992. 时间域电磁法原理. 长沙:中南工业大学出版社. |
[42] | 朴化荣. 1990. 电磁测深法原理. 北京:地质出版社. |
[43] | 滕吉文. 2006. 强化开展地壳内部第二深度空间金属矿产资源地球物理找矿、勘探和开发. 地质通报, 25(7):767-771. |
[44] | 王显祥, 底青云, 唐静. 2015. 电性源瞬变电磁初探. 地球物理学进展, 30(2):872-877, doi:10.6038/pg20150253. |
[45] | 薛国强, 李貅, 底青云. 2007. 瞬变电磁法理论与应用研究进展. 地球物理学进展, 22(4):1195-1200. |
[46] | 薛国强, 陈卫营, 周楠楠等. 2013. 接地源瞬变电磁短偏移深部探测技术. 地球物理学报, 56(1):255-261, doi:10.6038/cjg20130126. |
[47] | 闫述, 陈明生, 傅君眉. 2002. 瞬变电磁场的直接时域数值分析. 地球物理学报, 45(2):275-284. |
[48] | 周楠楠, 薛国强, 孔祥儒. 2012. 电性源瞬变电磁响应计算中的误差研究. 物探化探计算技术, 34(6):640-645. |