2. 中国矿业大学, 深部岩土力学与地下工程国家重点实验室, 徐州 221116
2. State Key Laboratory for Geomechanics & Deep Underground Engineering, Jiangsu Xuzhou 221116, China
随着国民经济的快速发展,瞬变电磁法(TEM)亦发展迅速,在矿井开采、资源勘察、水文地质调查、工程与环境等方面应用广泛(刘树才等,2005;刘盛东等,2014),随着探测精度要求的不断提高,对瞬变电磁的理论研究与资料解释也提出了更高的要求(牛之琏,2007;薛国强等,2007).目前,瞬变电磁法主要采用视电阻率值进行探测成果的推断解释,所采用的早期、晚期计算公式经过了简化、近似(苏朱刘和胡文宝,2002;严良俊等,2003),导致部分时间段视电阻率计算结果的畸变,为了使计算结果更加准确可靠,很多学者从不同方向进行了大量研究,包括全区视电阻率计算方法(严良俊等,1999;白登海等,2003;翁爱华等,2003;杨海燕等,2010)、拟地震波成像方法(李貅等,2007;张军等,2011;程久龙等,2013;戚志鹏等,2013)和微分电导解释方法(薛国强等,2003;苏茂鑫等,2010;戚志鹏等,2015;张军等,2015)等.
本文基于瞬变电磁场在不同介质中衰减扩散速度不同的物理现象,提出对探测到的感应电动势进行微分计算,通过分析计算结果来划分电性层以及判断低阻异常体.瞬变电磁勘探在每一个时间节点所获得的感应电动势是最原始、最直接、最基础的数据,因此在不做任何近似的情况下,对其衰减特征的研究可以保证数据的真实性与完整性.首先研究不同地电断面的感应电动势经微分计算所得值的特征,力求总结和归纳出瞬变电磁响应的衰减规律,然后将这一规律应用于瞬变电磁探测结果的电性层划分与异常解释,以期避开视电阻率计算,提高瞬变电磁解释的精度和准确性,丰富瞬变电磁解释方法.
1 瞬变电磁响应微分解释基本方法一般而言,瞬变电磁感应电动势ε总是随时间t的增大以指数规律衰减的,本文在此特征的基础上,引入TEM微分衰减速度K,公式为
(1) |
其中K为无量纲量.
在微分计算过程中,首先在双对数坐标系中对ε做关于t的插值,弥补采样间隔较大时带来的影响,继而将插值后的ε、t值代入公式(2):
(2) |
式(2) 中,εn代表第n个时间窗口接收到的感应电动势值,εn-1代表第n-1个时间窗口接收到的感应电动势值,tn代表第n个时间窗口,tn-1代表第n-1个时间窗口.
2 瞬变电磁响应微分衰减 2.1 均匀半空间瞬变电磁感应电动势经微分计算所得衰减速度与介质的电阻率有关,不同的电阻率在同一采样时间对应的衰减速度不同.为此,建立了5种不同电阻率的均匀半空间地电模型(表 1),数值计算(文中所有数值计算均基于重叠回线装置,回线边长100 m×100 m)其相应的瞬变电磁感应电动势响应,而后基于式(2) 对其微分计算,得到感应电动势的一阶微分,即微分衰减速度,研究衰减速度随电导率的变化关系有助于分析不同电阻率条件下瞬变电磁响应信号微分衰减速度的基本特征.
图 1为上述5种模型计算得到的微分衰减曲线.综合对比5条微分衰减曲线,不同电性介质的微分衰减趋势基本一致,衰减值在1.5~2.5之间,早期衰减值较小但变化幅度较大,随着时间的推移,衰减值增大,变化幅度逐渐减小,最终曲线平稳,趋于定值2.5.同一时间窗口,曲线上升部分,电阻率值越大则微分衰减值越大,且随时间的推移,这种趋势慢慢减小.对于同一种电性介质,瞬变电磁二次场在早期、晚期微分衰减速度不同,但是每个时间窗口的衰减值总是小于等于下一个时间窗口的衰减值,这种现象可看作瞬变电磁二次场衰减的固有规律.
均匀半空间介质中瞬变电磁二次场微分衰减速度存在固有的规律:随时间的推移,总是先上升后平稳,且在上升阶段,同一时间窗口下介质电阻率越低其微分衰减值就越低.基于此,以三层地电模型为基本模型,研究瞬变电磁二次场在非均匀介质中扩散的微分衰减响应规律.
图 2为同一观测系统下,改变覆盖层层厚时得到的三层模型(模型参数:H:电阻率ρ1=100 Ω·m,ρ2=10 Ω·m,ρ3=100 Ω·m;层厚:h1=50、60、70、80、90、100 m,h2=5 m, h3=∞)TEM微分衰减曲线图.图 3为同一观测系统下,改变中间层层厚时得到的三层模型(模型参数:H:电阻率ρ1=100 Ω·m,ρ2=10 Ω·m,ρ3=100 Ω·m;层厚:h1=100 m,h2=5、10、15、20、25、30 m, h3=∞)TEM微分衰减曲线图.
由图 2和图 3可知,整体上,曲线首支上升、中部下降、尾支上升并趋于定值,存在两个极值点.覆盖层厚度增加时,两个极值点均延时间轴正方向推移;中间层厚度增加时,极大值点出现时间基本一致,极小值点延时间轴正方向推移.根据瞬变电磁理论可知,当涡流环在覆盖层中向下及向外扩散时,仅受覆盖层电阻率影响,微分衰减速度逐渐增大;当涡流环扩散到覆盖层和中间层的交界面时,因中间层电阻率远低于覆盖层,微分衰减速度开始减小,表现在微分衰减曲线上为极大值点;随着涡流环继续向下及向外扩散,微分衰减速度受中间层影响程度加大,其值逐渐减小;当涡流环扩散到中间层与第三层交界面时,因第三层电阻率远高于中间层,微分衰减速度开始增大,表现在微分衰减曲线上为极小值点;随着涡流环继续向下及向外扩散,微分衰减速度受第三层影响程度加大,其值逐渐增加,因第三层厚度趋于无穷大,当涡流环扩散到一定深度,可以看作仅受第三层电阻率影响,此时涡流环的扩散可以看成是在均匀半空间中进行,微分衰减速度趋于定值.
由图 1、图 2及图 3可知,瞬变电磁二次场在同一种介质中扩散时,其微分衰减速度存在固有的规律:先上升后平稳,即每一个时间窗口的微分衰减值皆不小于前一个时间窗口.然而,当二次场扩散受到相对低阻异常的影响时,微分衰减速度发生突变,其值将会减小;同理,当二次场扩散受到相对高阻异常的影响时,微分衰减值将增加.总之,当瞬变电磁二次场由一种介质扩散到另一种介质时,其微分衰减速度将改变固有的规律,发生突变(对应图 2及图 3中微分衰减曲线的极值点),可依据该突变点进行电性划分.
3 瞬变电磁响应二阶微分衰减通过以上数值计算及分析结果可知,可依据瞬变电磁微分衰减曲线突变点进行电性划分,然而突变点并非等值点,不利于进行大批量数据的快速处理.从数学的角度发现微分衰减曲线的突变点的导数为0, 故,在瞬变电磁微分衰减速度的基础上,引入TEM二阶微分衰减速度K′,公式为
(3) |
其中K′为无量纲量.
在二阶微分计算时,将微分衰减速度K代入公式(4):
(4) |
式(4) 中,Kn代表第n个时间窗口对应的微分衰减速度值,Kn-1代表第n-1个时间窗口对应的微分衰减速度值,tn代表第n个时间窗口,tn-1代表第n-1个时间窗口.
3.1 水平层状地电模型基于上述分析,建立不同类型的水平层状地电模型(见表 2),通过数值计算得到其二阶微分衰减速度,并绘制二阶微分衰减曲线图进行对比分析.
图 4即为在表 2中不同地电模型的基础上得到的计算结果,由图可知:
(1) 均匀半空间模型下所得微分衰减曲线先上升后平稳,不存在突变点,对应的二阶微分衰减值均不小于零;
(2) 二层地电模型下所得微分衰减曲线先上升,后下降,再上升,存在两个突变点,即极值点,对应的二阶微分衰减曲线存在两个零值点,其中,极大值点可作为电性划分的依据,极小值点并非由二次场扩散到新的介质中产生,而是二次场为恢复其固有衰减规律的表现;
(3) 三层地电模型下所得微分衰减曲线先上升,后下降,再上升(上升幅度大于二层模型),存在两个突变点,即极值点,对应的二阶微分衰减曲线存在两个零值点,均可作为电性划分的依据;
(4) 四层地电模型下所得微分衰减曲线存在两个突变点(一为拐点,一为极小值点),对应的二阶微分衰减曲线存在两个零值点,然而拐点为前两层高阻电性层与第三层低阻电性层的综合反映,不能作为电性划分的依据,极小值点及其对应的零值点可作为电性划分的依据.
3.2 倾斜层状地电模型基于上述分析,三层倾斜层状地电模型(如图 5,模型参数:H:ρ1=100 Ω·m,ρ2=10 Ω·m,ρ3=100 Ω·m),其中覆盖层与中间层交界面为倾斜界面,中间层与第三层交界面为水平界面,覆盖层与中间层深度总和为150 m(1号测点:h1=60 m,h2=90 m;9号测点:h1=60 m,h2=90 m).
图 6为在此模型的基础上得到的研究结果,图 6b为微分衰减速度关于时间的曲线图,由图可见,曲线极大值点出现的时间与覆盖层厚度呈正相关,极小值点出现的时间基本保持不变,与中间层底面相吻合;图 6c为二阶微分衰减速度关于时间的曲线图,由图可见,图 6b中曲线极值点与图 6c中曲线的“0”值点相对应,代表不同电性层的分界面,为电性分层提供了一个固定的参数.
为了更直观的划分各电性层,参考目前常用的视电阻率拟断面图,绘制TEM二阶微分衰减拟断面图(图 7).由图 7可知,以“0”等值线为界面划分出3层电性层,与图 5中的电性模型相吻合.其中蓝色部分值为负,代表低阻层,红色部分值为正,代表高阻层,高、低阻电性层交界面为“0”.因此,以TEM二阶微分衰减拟断面图为基础进行数据解释是可行的.
为了验证瞬变电磁微分解释方法的效果,设计物理模拟实验进行验证,本次物理实验模型采用铜棒(参数:直径=15 cm,长度=30 cm)为目标体,水平放置于水槽(参数:200 cm×130 cm×90 cm)中,铜棒中心线距水面50 cm(图 7),使用terraTEM瞬变电磁仪沿水槽长边中心线探测,测点间距10 cm,实验装置为重叠回线装置(参数:20 cm×20 cm,50匝).
图 9即为在图 8中的模型基础上得到的TEM视电阻率拟断面图和二阶微分衰减拟断面图,由图 9a可知,在8~12号测点对应位置存在横向上的相对低阻异常区,但却需要依据个人经验来圈定异常体.由图 9b可以清晰地看到由“0”等值线划分的水槽底面和水面,在8~13号测点对应位置可以看到受铜棒影响横向变化明显,可较为清晰的圈定异常体,红色线框即为圈定的铜棒位置,由此可见,瞬变电磁微分解释方法在电性界面划分方面有显著效果,可以应用于瞬变电磁资料解释.
为了验证瞬变电磁微分解释方法的实际应用效果,以山西盂县兴峪煤矿地面瞬变电磁探测试验为例进行验证,本次探测的目的是查明15105工作面15#煤层中因小煤窑无规则开采产生的积水采空区的影响范围.据已知资料显示,探测区域内大部分为第四系松散覆盖层,覆盖层以下为二叠系泥、砂岩与石炭系石灰岩,之下为煤层,部分区域为采空充水区,煤层以下为石灰岩基底.煤层采空充水区相对顶底板为较明显低阻异常.基于此,以大定源回线装置进行物探,以其中一条测线进行对比分析,发射线圈边长为600 m×600 m,接收线圈面积为20000 m2,图 10为TEM视电阻率拟断面图(a)与TEM二阶微分衰减速度拟断面图(b).由图可见,在图 10a中测线7-11号测点下方存在低阻异常区域,然而,在异常范围圈定时,一方面依靠相应地质资料,另一方面则依赖个人经验,没有固定的标准;在图 10b中可以看到蓝色区域所代表的低阻异常区域,两者对比可以看出TEM二阶微分衰减拟断面图,可为瞬变电磁电性异常具体划分提供较为清晰的边界,能够作为瞬变电磁资料解释的参考图件.
瞬变电磁法二次场在不同电性介质中具有不同的扩散衰减速度,对于同一种电性介质,瞬变电磁场在每个时间窗口的衰减值总是小于等于下一个时间窗口的衰减值.当感应的涡旋电流扩散到电性交界面时,其微分衰减速度改变固有的规律,发生突变,产生突变点,可作为电性划分的依据.
6.2瞬变电磁微分解释方法对瞬变电磁采集到的大地响应信号直接进行计算,避免了视电阻率公式计算中早、晚期的划分,从新的角度对全区信号进行分析,简单方便,可嵌入仪器中,便于对采样数据实时监测.
6.3模拟结果与应用实例表明,以TEM微分衰减速度和TEM二阶微分衰减速度两个参数为基础进行瞬变电磁资料解释,能够为电性异常具体划分提供较为清晰的边界,可作为瞬变电磁法的资料解释参照图件.
6.4目前瞬变电磁微分解释方法尚不能进行时深转换,即不能进行异常体深度上的定位,后续研究将逐渐完善这一缺陷.
致谢 感谢审稿专家提供的宝贵意见,感谢吕阿谈硕士提供的应用数据,同时感谢课题组的师兄弟对于物理实验的帮助.[] | BAI Deng-Hai, Meju M A, LU Jian, et al. 2003. Numerical calculation of all-time apparent resistivity for the central loop transient electromagnetic method[J]. Chinese Journal of Geophysics, 46(5): 697–704. DOI:10.3321/j.issn:0001-5733.2003.05.018 |
[] | CHENG Jiu-Long, QIU Hao, YE Yun-Tao, et al. 2013. Research on wave-field transformation and data processing of the mine transient electromagnetic method[J]. Journal of China Coal Society, 38(9): 1646–1650. |
[] | LI Xiu, QUAN Hong-Juan, XU A-Xiang, et al. 2003. Differential coefficient imaging of the longitudinal conductance in the transient electromagnetic sounding[J]. Coal Geology & Exploration, 31(3): 59–61. |
[] | LI Xiu, XUE Guo-Qiang, GUO Wen-Bo. 2007. Research progress in TEM pseudo-seismic imaging[J]. Progress in Geophysics, 22(3): 811–816. DOI:10.3969/j.issn.1004-2903.2007.03.023 |
[] | LIU Sheng-Dong, LIU Jing, YUE Jian-Hua. 2014. Development status and key problems of Chinese mining geophysicaltechnology[J]. Journal of China Coal Society, 39(1): 19–25. |
[] | LIU Shu-Cai, LIU Zhi-Xin, JIANG Zhi-Hai. 2005. Application of TEM in hydrogeological prospecting of mining district[J]. Journal of China University of Mining & Technology, 34(4): 414–417. |
[] | NIU Zhi-Lian. 2007. Principle of Time Domain Electromagnetic Method[M]. Changsha: Central South University Press. |
[] | QI Zhi-Peng, LI Xiu, GUO Jian-Lei, et al. 2015. Fast imaging of correlation stack for airborne TEM based on differential conductance[J]. Progress in Geophysics, 30(4): 1903–1911. DOI:10.6038/pg20150451 |
[] | QI Zhi-Peng, LI Xiu, WU Qiong, et al. 2013. A new algorithm for full-time-domain wave-field transformation based on transient electromagnetic method[J]. Chinese Journal of Geophysics, 56(10): 3581–3595. DOI:10.6038/cjg20131033 |
[] | SU Mao-Xin, LI Shu-Cai, XUE Yi-Guo, et al. 2010. TEM apparent longitudinal conductance interpretation in tunnel geological forecast[J]. Chinese Journal of Geotechnical Engineering, 32(11): 1722–1726. |
[] | SU Zhu-Liu, HU Wen-Bao. 2002. Pseudo-full-region apparent resistivity and its one-dimensional inversion for center-loop-line configuration TEM data[J]. Geophysical Prospecting for Petroleum, 41(2): 216–221. |
[] | WENG Ai-Hua, LU Dong-Hua, LIU Guo-Xing. 2003. Definition of whole zone apparent resistivity for transient electromagnetic method of current dipole source[J]. Coal Geology & Exploration, 31(3): 56–58, 59. |
[] | XUE Guo-Qiang, LI Xiu. 2008. The technology of TEM tunnel prediction imaging[J]. Chinese Journal of Geophysics, 51(3): 894–900. DOI:10.3321/j.issn:0001-5733.2008.03.031 |
[] | XUE Guo-Qiang, Li Xiu, Di Qing-Yun. 2007. The progress of TEM in theory and application[J]. Progress in Geophysics, 22(4): 1195–1200. DOI:10.3969/j.issn.1004-2903.2007.04.026 |
[] | YAN Liang-Jun, HU Wen-Bao, CHEN Qing-Li, et al. 1999. The estimation and fast inversion of all-time apparent resistivities in long-offset transient electromagnetic sounding[J]. Oil Geophysical Prospecting, 34(5): 532–538. |
[] | YAN Liang-Jun, XU Shi-Zhe, HU Wen-Bao, et al. 2003. All-time apparent vertical conductance interpretation method for central loop transient electromagnetic sounding[J]. Journal of Zhejiang University (Science Edition), 30(2): 236–240. |
[] | YANG Hai-Yan, DENG Ju-Zhi, ZHANG Hua, et al. 2010. Research on full-space apparent resistivity interpretation technique in mine transient electromagneticmethod[J]. Chinese Journal of Geophysics, 53(3): 651–656. DOI:10.3969/j.issn.0001-5733.2010.03.020 |
[] | ZHANG Jun, LI Xiu, ZHAO Ying, et al. 2011. A technology research of high-resolation imaging for the transient electromagnetic pseudo wave field[J]. Progress in Geophysics, 26(3): 1077–1084. DOI:10.3969/j.issn.1004-2903.2011.03.038 |
[] | ZHANG Jun, MA Suo-Ni, LI Xiu. 2013. A study on differential conductance depth correction based on TEM[J]. Chinese Journal of Engineering Geophysics, 12(4): 438–444. |
[] | 白登海, MejuMA, 卢健, 等. 2003. 时间域瞬变电磁法中心方式全程视电阻率的数值计算[J]. 地球物理学报, 46(5): 697–704. DOI:10.3321/j.issn:0001-5733.2003.05.018 |
[] | 程久龙, 邱浩, 叶云涛, 等. 2013. 矿井瞬变电磁法波场变换与数据处理方法研究[J]. 煤炭学报, 38(9): 1646–1650. |
[] | 李貅, 全红娟, 许阿祥, 等. 2003. 瞬变电磁测深的微分电导成像[J]. 煤田地质与勘探, 31(3): 59–61. |
[] | 李貅, 薛国强, 郭文波. 2007. 瞬变电磁法拟地震成像研究进展[J]. 地球物理学进展, 22(3): 811–816. DOI:10.3969/j.issn.1004-2903.2007.03.023 |
[] | 刘树才, 刘志新, 姜志海. 2005. 瞬变电磁法在煤矿采区水文勘探中的应用[J]. 中国矿业大学学报, 34(4): 414–417. |
[] | 刘盛东, 刘静, 岳建华. 2014. 中国矿井物探技术发展现状和关键问题[J]. 煤炭学报, 39(1): 19–25. |
[] | 牛之琏. 2007. 时间域电磁法原理[M]. 长沙: 中南大学出版社. |
[] | 戚志鹏, 李貅, 郭建磊, 等. 2015. 基于微分电导的航空瞬变电磁合成孔径快速成像方法研究[J]. 地球物理学进展, 30(4): 1903–1911. DOI:10.6038/pg20150451 |
[] | 戚志鹏, 李貅, 吴琼, 等. 2013. 从瞬变电磁扩散场到拟地震波场的全时域反变换算法[J]. 地球物理学报, 56(10): 3581–3595. DOI:10.6038/cjg20131033 |
[] | 苏茂鑫, 李术才, 薛翊国, 等. 2010. 隧道地质预报中的瞬变电磁视纵向电导解释方法研究[J]. 岩土工程学报, 32(11): 1722–1726. |
[] | 苏朱刘, 胡文宝. 2002. 中心回线方式瞬变电磁测深虚拟全区视电阻率和一维反演方法[J]. 石油物探, 41(2): 216–221. |
[] | 翁爱华, 陆冬华, 刘国兴. 2003. 利用连分式定义瞬变电磁法全区视电阻率研究[J]. 煤田地质与勘探, 31(3): 56–58, 59. |
[] | 薛国强, 李貅. 2008. 瞬变电磁隧道超前预报成像技术[J]. 地球物理学报, 51(3): 894–900. DOI:10.3321/j.issn:0001-5733.2008.03.031 |
[] | 薛国强, 李貅, 底青云. 2007. 瞬变电磁法理论与应用研究进展[J]. 地球物理学进展, 22(4): 1195–1200. DOI:10.3969/j.issn.1004-2903.2007.04.026 |
[] | 严良俊, 胡文宝, 陈清礼, 等. 1999. 长偏移距瞬变电磁测深的全区视电阻率求取及快速反演方法[J]. 石油地球物理勘探, 34(5): 532–538. |
[] | 严良俊, 徐世浙, 胡文宝, 等. 2003. 中心回线瞬变电磁测深全区视纵向电导解释方法[J]. 浙江大学学报(理学版), 30(2): 236–240. |
[] | 杨海燕, 邓居智, 张华, 等. 2010. 矿井瞬变电磁法全空间视电阻率解释方法研究[J]. 地球物理学报, 53(3): 651–656. DOI:10.3969/j.issn.0001-5733.2010.03.020 |
[] | 张军, 李貅, 赵莹, 等. 2011. 瞬变电磁虚拟波场高分辨成像技术研究[J]. 地球物理学进展, 26(3): 1077–1084. DOI:10.3969/j.issn.1004-2903.2011.03.038 |
[] | 张军, 马索尼, 李貅. 2015. 瞬变电磁微分电导深度校正技术研究[J]. 工程地球物理学报, 12(4): 438–444. |