城市地区充水薄层煤炭采空区瞬变电磁法探测效能分析.pdf

第5 3卷 第1 1期 2 0 2 3年1 1月中 国 海 洋 大 学 学 报P E R I O D I C A L O F O C E A N U N I V E R S I T Y O F C H I N A5 3(1 1):1 3 1 1 3 9N o v.,2 0 2 3城市地区充水薄层煤炭采空区瞬变电磁法探测效能分析孙兴昊1,罗平凡2,姜巧巧3 ,钱振嵩2,郭秀军1,4(1.中国海洋大学环境科学与工程学院,山东 青岛 2 6 6 1 0 0;2.山东轨道交通勘察设计院有限公司,山东 济南 2 5 0 0 1 4;3.山东省地矿工程勘察院,山东 济南 2 5 0 0 1 4;4.山东省海洋环境地质工程重点实验室,山东 青岛 2 6 6 1 0 0)摘 要:本文通过瞬变电磁法对薄层充水采空区进行实地探测,并使用天然源面波法进行二次探测,通过天然源面波法的探测结果和钻探的结果验证了瞬变电磁法对薄层充水采空区探测的可行性。通过数值模拟,构建不同采空区含水率、不同深度和不同厚度的薄层充水采空区的电阻率模型,研究了瞬变电磁法对薄层充水采空区的探测效果。研究结果表明,瞬变电磁法可以对薄层充水采空区进行有效的探测;随着薄层充水采空区含水率的降低,瞬变电磁法探测出的异常区范围和电阻率异常幅值变小,当采空区未充水时,瞬变电磁法无法对采空区进行有效的探测;随着薄层充水采空区深度的增加,瞬变电磁法探测出的异常范围和电阻率异常幅值变小,在本文模拟的工况条件下,1 0 0 m为探测的极限深度;随着薄层充水采空区厚度的降低,探测出的异常区范围和电阻率异常幅值变小。关键词:薄层充水采空区;瞬变电磁法;数值模拟;探测效果中图法分类号:P 6 3 1.2 文献标志码:A 文章编号:1 6 7 2-5 1 7 4(2 0 2 3)1 1-1 3 1-0 9D O I:1 0.1 6 4 4 1/j.c n k i.h d x b.2 0 2 2 0 2 9 5引用格式:孙兴昊,罗平凡,姜巧巧,等.城市地区充水薄层煤炭采空区瞬变电磁法探测效能分析J.中国海洋大学学报(自然科学版),2 0 2 3,5 3(1 1):1 3 1-1 3 9.S u n X i n g h a o,L u o P i n g f a n,J i a n g Q i a o q i a o,e t a l.A n a l y s i s o n d e t e c t i o n e f f i c i e n c y o f t r a n s i e n t e l e c t r o m a g n e t i c m e t h o d i n g o a f o f w a t e r f i l l e d t h i n c o a l i n u r b a n a r e aJ.P e r i o d i c a l o f O c e a n U n i v e r s i t y o f C h i n a,2 0 2 3,5 3(1 1):1 3 1-1 3 9.基金项目:济南城市轨道交通8号线一期工程沿线专项勘察与地面变形监测项目技术总结与创新研究项目(2 0 2 1 0 4 3 6)资助S u p p o r t e d b y t h e T e c h n i c a l S u mm a r y a n d I n n o v a t i o n R e s e a r c h o f S p e c i a l I n v e s t i g a t i o n a n d G r o u n d D e f o r m a t i o n M o n i t o r i n g P r o j e c t A l o n g t h e F i r s t P h a s e o f J i n a n U r b a n R a i l T r a n s i t L i n e 8(2 0 2 1 0 4 3 6)收稿日期:2 0 2 2-0 5-2 6;修订日期:2 0 2 2-0 6-3 0作者简介:孙兴昊(1 9 9 9),男,硕士生,从事地球物理方向研究。E-m a i l:3 1 2 6 0 4 2 8 9 q q.c o m 通信作者:E-m a i l:8 2 6 9 6 5 8 1 9 q q.c o m 山东省济南市圣井煤矿矿区位于济南市章丘区境内,属于冲洪积平原地貌,地面较为平整,地形起伏较小。区域地层由老到新依次为:奥陶系中统马家沟组,石炭系中统本溪组、上统太原组,二叠系下统山西组及第四系1。本文研究煤层为3煤,该煤层位于二叠系下统山西组中,煤层埋深为5 0 7 0 m,煤层厚度为0.7 1.5 m,开采厚度平均为1 m,顶板为砂岩、粉砂岩,底板为砂岩。根据前期勘探资料,该煤层已被完全开采完,形成采空区。根据规划,济南地铁8号线拟从该区域穿过。采空区的存在会严重制约地下工程建设,甚至造成地表开裂、沉陷、变形,结构物开裂、损毁。在地下水资源丰富的地区,采空区充水现象普遍存在。由于充水采空区与周围地质环境存在明显的物性差异,因此,可以利用地球物理勘探技术对充水采空区进行探测2。近年来,对充水采空区的地球物理探测技术较为成熟,但已有研究中采空区厚度多在5 2 0 m之间3-4,对于厚度小于2 m的薄层煤矿采空区的探测和研究较少,因此本文针对薄层煤层采空区探测的可行性展开研究。目前,常用的充水采空区地球物理探测方法主要包括瞬变电磁法5-8、高密度电法9、大地电磁法1 0、天然源面波法1 1等,各种物探方法均在探测充水采空区方面取得了一定的效果4。当煤层被采空后,若采空区充水则采空区表现为低阻体异常特征。瞬变电磁法对低阻异常体反应极其灵敏,横向探测分辨率较高,施工方便,工作效率高5。天然源面波法分辨率较高,抗干扰能力强,探测深度较大,不受地形限制,野外施工便捷,适用于交通繁忙、人口密集、电磁干扰严重和震动复杂的城市环境1 1。本文以济南圣井矿区作为研究区,选择瞬变电磁法为主要探测方法对该矿区进行探测,并用天然源面波法进行二次探测,通过天然源面波法的探测结果对瞬变电磁法的探测结果进行验证,同时,结合钻探结果进行验证,分析瞬变电磁法的探测能力。然后根据已知资料建立简化模型,采用有限元法针对充水采空区进行正演模拟,通过二分法计算全区视电阻率,使用烟圈法对全区视电阻率进行反演,分析瞬变电磁法对薄层煤炭充水采空区探测的剖面的异常特征。然后分别模拟不同深度、不中 国 海 洋 大 学 学 报2 0 2 3 年同充水率和不同厚度下瞬变电磁法对薄层采空区的探测,确定瞬变电磁法对薄层采空区的探测能力。1 探测原理和探测设备1.1 瞬变电磁法瞬变电磁法是通过不接地回线或者接地线源向地下发射脉冲电磁波产生一次电磁场,在一次电磁场的作用下,地下导体内部感应产生感应涡流,感应涡流产生的二次电磁场随时间以等效涡流环的形式向下传播、向外扩展,通过线圈或者接地电极观测二次电磁场,得到二次电磁场随时间衰减的特性曲线,根据二次电磁场随着时间的衰减变化特征曲线判断地下地质体的电性特征,从而探明不良地质体的分布情况1 2-1 3。通过对探测区进行现场勘查,本文选取中心回线装置进行探测。中心回线装置是指接收线框位于发射线框中心,在测量过程中两个线框沿探测剖面同步移动进行数据采集的一种瞬变电磁法装置。本次瞬变电磁探测选用加拿大P h o e n i x公司生产的V 8多功能电法仪,设定发射电流1 0 A;关断时间为1 0 m s;5 0 H z频率干扰压制;选用2 0 m 2 0 m的矩形发射线框,共3匝;选用直径为1 m的圆形接收线框,共1 0 0匝。1.2 天然源面波法天然源面波法的信号场源为地球内部系统产生的震动,根据面波频散原理,波速会随着频率和波长的变化而变化,通过从天然场的微动信号中提取出瑞利面波的频散曲线,反演频散曲线得到地下横波波速剖面,从而推测地下不良地质体的分布1 1,1 4。通过对探测区进行现场勘查,采用“L”型台阵布设检波器进行探测,在水平方向和垂直方向共布设8个检波器进行探测,道间距为1 0 m。本次天然源面波探测选用北京市水电物探研究所生产的WD-2无线智能微动勘探系统,检波器布置L型台阵,道间距为1 0 m,长边5 0 m,短边2 0 m;接收频率段为23 1.5 H z;接收道数为8道;采样间隔为4 m s;采样时间不少于1 0 m i n。2 瞬变电磁法实地探测效果2.1 探测测线布设为探明圣井矿区采空区的分布情况,根据现场探测条件,共布设2条测线进行探测,其中,S J 1测线沿济南市经十路方向横向布设,长为2 4 2 0 m,进行瞬变电磁法和天然源面波法探测;S J 3测线垂直经十路方向纵向布设,长为3 8 0 m,进行瞬变电磁法探测,S J 1和S J 3测线布设如图1所示。图1 探测工作部署图 F i g.1 D e p l o y m e n t d i a g r a m o f d e t e c t i o n w o r k2.2 探测成果解译将通过瞬变电磁法采集的各个测点数据进行视电阻率和视深度转换,并对数据进行滤波、一维反演,绘制视电阻率等值线图。探测区域S J 1线1 6 0 0 2 2 0 0 m视电阻率剖面图如图2所示。探测区域S J 3线4 03 2 0 m视电阻率剖面图如图3所示。图2 S J 1测线实测视电阻率等值线图F i g.2 C o n t o u r m a p o f m e a s u r e d a p p a r e n t r e s i s t i v i t y o f S J 1 s u r v e y l i n e2311 1期孙兴昊,等:城市地区充水薄层煤炭采空区瞬变电磁法探测效能分析图3 S J 3测线实测视电阻率等值线图F i g.3 C o n t o u r m a p o f m e a s u r e d a p p a r e n t r e s i s t i v i t y o f S J 3 s u r v e y l i n e由图2中可以看出,在S J 1线视电阻率等值线图中,从地表到高程7 5 m之间的空白区域属于瞬变电磁法探测的盲区,无法获取有效信息。结合地质资料,综合分析认为,在高程7 5 m以下,等值线图横向上视电阻率呈近似平行分布,垂向上视电阻率值呈现“高-低-高”的变化,上部的高阻层对应3煤顶部砂岩、泥岩,中部的低阻层对应3煤采空区,下部的高阻层对应3煤采空区下部的砂岩、泥岩。3煤被开采完,当存在充水采空区时,在视电阻率等值线图上将反映为低阻的电性标志层。水平位置1 6 0 01 9 3 0 m、1 9 6 02 2 0 0 m,高程3 5 5 0 m,出现低阻异常,结合地质资料,分析认为这些异常区存在充水采空区。结合地质资料解译出1条断层,为F 1 1断层。由图3可以看出,在S J 3线水平位置4 03 2 0 m,高程2 5 4 0 m,出现低阻异常区,结合地质资料,分析认为该异常区存在充水采空区。从图2中可看出,在S J 1和S J 3两条测线的交汇处位置,高程3 0 5 5 m出现低阻异常区。从图3中可看出,在两测线交汇处位置,高程3 54 0 m出现低阻异常区。结合图2和图3可看出,S J 1和S J 3两条测线在交汇处探测出的低阻异常区相互对应。从天然源面波法采集的原始微动信号中提取出频散曲线,对频散曲线反演,绘制横波波速等值线图。探测区域S J 1线横波波速等值线图如图4。图4 S J 1测线横波波速等值线图F i g.4 T r a n s v e r s e w a v e v e l o c i t y c o n t o u r m a p o f S J 1 s u r v e y l i n e 从图4中可以看出,在水平位置1 6 5 0 1 9 0 0 m、1 9 3 0 2 2 0 0 m,高程3 5 5 5 m,出现“U”型异常区,横波波速值对比周围出现明显降低,该异常与图2中瞬变电磁法探测出的异常相互对应,结合地质资料,分析认为这些异常区存在充水采空区。结合地质资料解译出1条断层,为F 1 1断层。通过天然源面波法解译出的异常区对通过瞬变电磁法解译出的异常区进行了较好的验证,证实了瞬变电磁法对薄层充水采空区探测的可行性,说明瞬变电磁法可以对薄层充水采空区进行探测。对比瞬变电磁法和天然源面波法的探测结果,可以看出,瞬变电磁法探测结果中的异常区更为明显,这说明瞬变电磁法对薄层充水采空区的探测更为灵敏。为验证S J 1和S J 3线通过瞬变电磁法和天然源面波法探测出的采空区的准确性,布设地质钻探点进行钻孔验证,布设S J K 1、S J K 2钻孔作为本次探测验证钻孔,如图所示。其中,S J K 1孔在高程3 5.43 7.5 m,S J K 2孔在高程2 52 6.2 m出现漏水、掉钻、进尺变331中 国 海 洋 大 学 学 报2 0 2 3 年快、岩心破碎、塌孔、漏浆,认为这两个钻孔均揭露3煤采空区,且采空区充水。瞬变电磁法和天然源面波法解译出的充水采空区和钻孔验证的结果基本吻合,因此认为瞬变电磁法和天然源面波法对薄层充水采空区的探测结果准确、可靠。3 瞬变电磁法探测效果数值模拟分析3.1 全区视电阻率计算和烟圈法反演3.1.1全区视电阻率计算 本文通过感应电动势对全区视电阻率进行定义,利用二分搜索法对全区视电阻率进行计算。二分搜索法的基本思想,是在探测区域电阻率可能的取值范围m i n,m a x 内,二等分取值范围m i n,m a x,将其中点处的电阻率值m i d=(m i n+m a x)/2代入均匀半空间解析式,计算中点的感应电动势Vm i d(t)。如果实测t时刻的感应电动势V0(t)大于中点处的理论感应电动势Vm i d(t),根据感应电动势随电阻率增加而单调递减的特性,V0(t)对应的视电阻率值应在m i n,m i d 中,反之,则在m i d,m a x 中。以此类推,直到V0(t)与Vm i d(t)之差小于某一个阈值,停止计算,就求出了t时刻感应电动势V0(t)对应的全区视电阻率值1 5-1 6。3.1.2 烟圈法反演 M.N.N a b i g h i a n指出,当一次场的电流关断后,任一时刻,地下半空间介质中形成的感应涡流在地表产生的二次磁场可以等效为一个水平环状电流产生的磁场,地下涡流向外、向下扩散的现象,即为“烟圈效应”1 6-1 8。根据烟圈理论,t时刻烟圈的垂向传播深度d和垂向传播速度v为:d=4 t0。(1)v=2 0t。(2)式中:表示均匀半空间的电阻率;t表示采样时间;0表示真空磁导率。通过差分公式近似计算传播速度v:v=dt=dj-ditj-ti=4 0 tjj-tiitj-ti。(3)由(2)可得=v2t04。(4)将(3)代入(4),可得似电阻率r=4 tjj-tiitj-ti2tj i。(5)似电阻率对应的深度为:Hr=0.4 4 1dj+di2。(6)式中:ti、tj表示两个相邻的采样时间;tjti,i、j为相邻两个采样时间对应的全区视电阻率;tj i是ti和tj的算术平均值。3.2 数值模型建立根据探测研究区的地层特性和电性特征,通过C OM S O L软件系统建立简化的含薄层充水采空区的三维地电模型,进行正演计算。模型自上至下依次为空气层、第四系覆盖层和砂页岩层,采空区在砂页岩层中,位于第四系覆盖层正下方,模型具体情况如图5所示。参考该地区以往电测井及电测深资料对模型进行电阻率赋值。其中,空气层大小为3 0 0 m3 0 0 m2 0 0 m,电阻率设定为1 06 m。土层大小为3 0 0 m3 0 0 m 5 m,电阻率设定为1 0 0 m。砂岩层分为两层,采空区上方砂岩层大小为3 0 0 m3 0 0 m5 5 m,采空区下方砂岩层大小为3 0 0 m3 0 0 m1 4 0 m,电阻率设定为2 0 0 m。采空区大小为2 0 0 m 3 0 0 m1 m,位于砂岩层中部,采空区充水后因浸润效应使得电阻率下降1 3,电阻率设定为5 m,采空区左右两侧各有5 0 m 3 0 0 m1 m的煤层,煤层电阻率设定为2 5 0 m,采空区顶界埋深6 0 m,采空区X Z方向平面视图如图5所示。图5 采空区模型X Z方向平面视图F i g.5 P l a n e v i e w o f g o a f i n X Z d i r e c t i o n3.3 计算参数设置选用装置为中心回线装置,发射线框大小为2 0 m2 0 m,匝数为3匝,供电电流为1 0 A,关断时间为1 0 m s,接收线框半径为0.5 m,匝数为1 0 0匝。测线沿4311 1期孙兴昊,等:城市地区充水薄层煤炭采空区瞬变电磁法探测效能分析采空区正上方布设,测线长为2 0 0 m。测点位于接收线框中心,第1个测点坐标为(-1 0 0,0,0),测点沿x轴正方向移动,点距为2 0 m,共计1 1个测点。正演计算采用有限元法,计算得到各个测点的感应电动势衰减曲线,通过二分搜索法计算全区视电阻率,并对全区视电阻率进行烟圈反演。3.4 数值模拟结果及分析烟圈反演结果如图6所示。从图6中可以看出,在水平位置-1 0 01 0 0 m,埋深5 08 0 m,出现低阻异常,该低阻异常即为模型正演设定的充水采空区。瞬变电磁法探测的横向边界与模型预设较为吻合,纵向边界较模型预设边界扩大3 0倍,横向分辨率较高,纵向分辨率较低,这是瞬变电磁法固有的缺点。当1 m厚度的充水采空区埋深6 0 m时,瞬变电磁法能对充水采空区进行有效的探测。图6 厚度为1 m、顶界埋深6 0 m的充水采空区模型视电阻率剖面图F i g.6 A p p a r e n t r e s i s t i v i t y p r o f i l e o f w a t e r f i l l e d g o a f m o d e l w i t h t h i c k n e s s o f 1 m a n d t o p b o u n d a r y b u r i e d d e p t h o f 6 0 m3.5 采空区不同充水率对探测效果的影响为探究薄层充水采空区充水率不同时,瞬变电磁法对充水采空区的探测效果的差异,通过改变采空区的充水率,模拟了采空区充水率为5 0%和0的瞬变电磁探测。本次模拟将采空区分为上下两层,上层为空气层,电阻率为1 06 m,下层为水层,电阻率为5 m,当采空区充水率为5 0%时,上层空气层厚度和下层水层厚度均为0.5 m,当采空区充水率为0时,空气层厚度为1 m,不含水层。反演结果如图7和图8所示。从图7可以看出,在水平位置-4 04 0 m,埋深5 0 6 5 m,出现低阻异常,该低阻异常与模型设定的6 0 m埋深的充水率为5 0%的采空区相对应。从图中可以看出,当采空区充水率降到5 0%时,异常区仍为低阻异常区,异常区相较充水率为1 0 0%时变小,瞬变电磁法探测的左右侧边界相较模型预设边界偏差6 0 m,纵向边界较模型预设边界扩大1 5倍。从图8可以看出,图中无明显异常区。这说明当采空区充水率降到0时,瞬变电磁法无法对薄层采空区进行探测。图7 5 0%充水率的采空区模型视电阻率剖面图F i g.7 A p p a r e n t r e s i s t i v i t y p r o f i l e o f g o a f m o d e l w i t h 5 0%w a t e r f i l l i n g r a t e图8 不充水的采空区模型视电阻率剖面图F i g.8 A p p a r e n t r e s i s t i v i t y p r o f i l e o f g o a f m o d e l w i t h o u t w a t e r f i l l i n g3.6 采空区深度对探测效果的影响为探究薄层充水采空区不同埋深时,瞬变电磁法对充水采空区的探测效果的差异,本文进一步增加薄531中 国 海 洋 大 学 学 报2 0 2 3 年层充水采空区的埋深,模拟了充水采空区顶部埋深在8 0、1 0 0和1 2 0 m的瞬变电磁探测,反演结果如图9、图1 0和图1 1所示。图9 厚度为1 m、顶界埋深8 0 m的充水采空区模型视电阻率剖面图F i g.9 A p p a r e n t r e s i s t i v i t y p r o f i l e o f w a t e r f i l l e d g o a f m o d e l w i t h t h i c k n e s s o f 1 m a n d t o p b o u n d a r y b u r i e d d e p t h o f 8 0 m图1 0 厚度为1 m、顶界埋深1 0 0 m的充水采空区模型视电阻率剖面图F i g.1 0 A p p a r e n t r e s i s t i v i t y p r o f i l e o f w a t e r f i l l e d g o a f m o d e l w i t h t h i c k n e s s o f 1 m a n d t o p b o u n d a r y b u r i e d d e p t h o f 1 0 0 m从图9可以看出,在水平位置-8 08 0 m,埋深6 5 8 5 m,出现低阻异常,该低阻异常与模型设定的8 0 m埋深的充水采空区相对应。从图中可以看出,当模型的采空区埋深增加到8 0 m时,异常区相较采空区埋深在6 0 m时变小,瞬变电磁法探测的左右侧边界较模型预设边界偏差2 0 m,纵向边界较模型预设边界扩大2 0倍。从图1 0可以看出,在水平位置-6 0 6 0 m,埋深8 0 1 0 0 m,出现低阻异常,该低阻异常与模型设定的1 0 0 m埋深的充水采空区相对应。从图1 0中可以看出,当模型的采空区埋深增大到1 0 0 m时,异常区相较采空区埋深在6 0 m时进一步变小,瞬变电磁法探测的左右侧边界较模型预设边界偏差4 0 m,纵向边界较模型预设边界扩大2 0倍。从图1 1可以看出,当模型的采空区埋深增大到1 2 0 m时,图中无明显异常区,通过瞬变电磁法已经难以对采空区进行探测。综上,在本文设置的模拟条件下,可以看出瞬变电磁法对1 m厚的薄层充水采空区的探测极限深度为1 0 0 m。图1 1 厚度为1 m、顶界埋深1 2 0 m的充水采空区模型视电阻率剖面图F i g.1 1 A p p a r e n t r e s i s t i v i t y p r o f i l e o f w a t e r f i l l e d g o a f m o d e l w i t h t h i c k n e s s o f 1 m a n d t o p b o u n d a r y b u r i e d d e p t h o f 1 2 0 m3.7 采空区厚度对探测结果的影响为探究薄层充水采空区不同厚度时,瞬变电磁法对充水采空区的探测效果的差异,通过改变采空区的厚度,模拟了采空区厚度为0.5和2 m的瞬变电磁探测。反演结果如图1 2和图1 3所示。从图1 2可以看出,在水平位置-4 04 0 m,埋深5 0 6 5 m,出现低阻异常,该低阻异常与模型设定的0.5 m厚度的充水采空区相对应。从图中可以看出,当模型采空区的厚度为0.5 m时,异常区相较采空区厚度为1 m时变小,瞬变电磁法探测的左右侧边界较模型预设边界偏差6 0 m,纵向边界较模型预设边界扩大3 0倍。从下图1 3可以看出,在水平位置-1 0 0 1 0 0 m,6311 1期孙兴昊,等:城市地区充水薄层煤炭采空区瞬变电磁法探测效能分析埋深5 0 1 0 0 m,出现低阻异常,该低阻异常与模型设定的2 m厚度的充水采空区相对应。从图中可以看出,当模型采空区的厚度为2 m时,异常区相较采空区厚度为1 m时变大,瞬变电磁法探测的横向边界与模型预设较为吻合,纵向边界较模型预设边界扩大2 5倍。在本文的模拟条件下,当采空区厚度降为0.5 m时,瞬变电磁法仍能对其进行较为有效的探测。图1 2 厚度为0.5 m、顶界埋深6 0 m的充水采空区模型视电阻率F i g.1 2 A p p a r e n t r e s i s t i v i t y p r o f i l e o f w a t e r f i l l e d g o a f m o d e l w i t h t h i c k n e s s o f 0.5 m a n d t o p b o u n d a r y b u r i e d d e p t h o f 6 0 m图1 3 厚度为2 m、顶界埋深6 0 m的充水采空区模型视电阻率F i g.1 3 A p p a r e n t r e s i s t i v i t y p r o f i l e o f w a t e r f i l l e d g o a f m o d e l w i t h t h i c k n e s s o f 2 m a n d t o p b o u n d a r y b u r i e d d e p t h o f 6 0 m4 结论和建议(1)通过瞬变电磁法和天然源面波法两种方法相结合的探测方法能够实现对薄层充水采空区的探测,天然源面波法探测的结果很好地验证了瞬变电磁法探测的结果,探测结果经过钻孔验证,取得了较好的效果。在瞬变电磁法和天然源面波法的数据解译过程中,需要结合已有的地质资料来对异常区进行判定,通过综合分析的思路来解译成果,才能确保探测的准确性。(2)通过已知资料建立简化模型,模拟实际探测的计算结果表明,当厚度为1 m的充水采空区顶界埋深为6 0 m时,瞬变电磁法能对充水采空区进行有效的探测。(3)模拟计算结果表明,当薄层采空区埋深6 0 m时,随着采空区充水率的降低,瞬变电磁法的探测效果越来越差,当充水率降到5 0%时,反演结果异常区的范围变小,当充水率降到0时,无法通过反演结果对采空区的位置进行划分。这说明瞬变电磁法对低阻异常区的探测极其灵敏,对高阻区的探测极不灵敏。(4)模拟计算结果表明,当1 m厚度的薄层充水采空区埋深6 0 m时,瞬变电磁法的探测效果较好,反演的结果能较为准确地分辨出横向边界,反演出的异常区的厚度比比正演设定的厚度大。随着采空区深度增加,异常区范围逐渐变小,横向边界的分辨率逐渐降低,当采空区深度达到1 2 0 m时,瞬变电磁法无法探测出薄层采空区的存在。(5)模拟计算结果表明,当薄层充水采空区埋深为6 0 m时,随着采空区厚度的减小,瞬变电磁法的探测效果越来越差,异常区范围越来越小。当厚度为2 m时,能较为准确地分辨出横向边界,随采空区厚度的降低,异常区的范围逐渐减小,横向边界的分辨率逐渐降低。参考文献:1 徐小连.瞬变电磁法在山东章丘某产业园采空区勘探中的应用J.低碳世界,2 0 1 8(2):2 8-2 9.X u X i a o l i a n.A p p l i c a t i o n o f t r a n s i e n t e l e c t r o m a g n e t i c m e t h o d i n g o a f e x p l o r a t i o n o f a n i n d u s t r i a l p a r k i n Z h a n g q i u,S h a n d o n g P r o v-i n c e J.L o w C a r b o n W o r l d,2 0 1 8(2):2 8-2 9.2 平守国,王永增,张忠海,等.齐大山铁矿南帮含水构造带探测研究J.金属矿山,2 0 2 0(1):5 6-6 2.P i n g S h o u g u o,W a n g Z e n g p i n g,Z h a n g Z h o n g h a i,e t a l.S t u d y o n t h e d e t e c t i o n o f w a t e r-b e a r i n g s t r u c t u r a l z o n e s i n t h e s o u t h s i d e o f Q i d a s h a n i r o n m i n eJ.M e t a l M i n e,2 0 2 0(1):5 6-6 2.3 覃庆炎.瞬变电磁法在充水采空区探测中的应用J.煤炭科学技术,2 0 1 4,4 2(8):1 0 9-1 1 2.Q i n Q i n g y a n.A p p l i c a t i o n o f t r a n s i e n t e l e c t r o m a g n e t i c m e t h o d f o r w a t e r a c c u m u l a t e d g o a f i n c o a l m i n e sJ.C o a l S c i e n c e a n d T e c h-731中 国 海 洋 大 学 学 报2 0 2 3 年n o l o g y,2 0 1 4,4 2(8):1 0 9-1 1 2.4 付天光.综合物探方法探测煤矿采空区及充水区技术研究J.煤炭科学技术,2 0 1 4,4 2(8):9 0-9 4.F u T i a n g u a n g.S t u d y o n t e c h n o l o g y o f c o m p r e h e n s i v e g e o p h y s i c a l m e t h o d e x p l o r a t i o n o f m i n e g o a f a n d w a t e r a c c u m u l a t e d a r e aJ.C o a l S c i e n c e a n d T e c h n o l o g y,2 0 1 4,4 2(8):9 0-9 4.5 李宏杰.浅层地震和瞬变电磁法在采空区探测中的应用研究J.煤矿开采,2 0 1 3(1):1 7-1 9+2 7.L i H o n g j i e.A p p l i c a t i o n o f s h a l l o w e a r t h q u a k e a n d t r a n s i e n t e l e c-t r o m a g n e t i c m e t h o d i n g o b e x p l o r a t i o nJ.C o a l M i n i n g T e c h n o l o-g y,2 0 1 3(1):1 7-1 9+2 7.6 郄卫东.瞬变电磁法在探测煤矿采空区及塌陷区中的应用J.物探与化探,2 0 1 2(3 6):9 2-9 5.X i W e i d o n g.T h e a p p l i c a t i o n o f t r a n s i e n t e l e c t r o m a g n e t i c m e t h o d t o t h e d e t e c t i o n o f c o a l m i n e g o a f a n d c o l l a p s e a r e a sJ.G e o p h y s i-c a l a n d G e o c h e m i c a l E x p l o r a t i o n,2 0 1 2(3 6):9 2-9 5.7 梁爽,李志民.瞬变电磁法在阳泉二矿探测积水采空区效果分析J.煤田地质与勘探,2 0 0 3,3 1(4):4 9-5 1.L i a n g S h u a n g,L i Z h i m i n.A n a l y s i s o f t h e e f f e c t b y u s i n g T E M d e-t e c t i n g w a t e r f i l l e d g o b a r e a i n N o.2 M i n e,Y a n g q u a n M i n e G r o u pJ.C o a l G e o l o g y&E x p l o r a t i o n,2 0 0 3,3 1(4):4 9-5 1.8 陈江军,刘波,王鹏来,等.综合物探方法在石膏矿采空区探测中的应用研究J.地球物理学进展,2 0 1 8,3 3(2):7 8 3-7 8 9.C h e n J i a n g j u n,L i u B o,W a n g P e n g l a i,e t a l.A p p l i c a t i o n r e s e a r c h o f i n t e g r a t e d g e o p h y s i c a l e x p l o r a t i o n t e c h n i q u e s t o g y p s u m m i n e g o bJ.P r o g r e s s i n G e o p h y s i c s,2 0 1 8,3 3(2):7 8 3-7 8 9.9 林井祥,张立明,张继忠.基于瞬变电磁法和大地电磁法的煤矿积水采空区探测J.煤炭工程,2 0 2 1,5 3(6):1 5 2-1 5 6.L i n J i n g x i a n g,Z h a n g L i m i n g,Z h a n g J i z h o n g.D e t e c t i o n o f g o a f h y d r o p s i n c o a l m i n e b a s e d o n t r a n s i e n t e l e c t r o m a g n e t i c m e t h o d a n d m a g n e t o t e l l u r i c m e t h o dJ.C o a l E n g i n e e r i n g,2 0 2 1,5 3(6):1 5 2-1 5 6.1 0 郭伟立,薛国强,周楠楠,等.利用瞬变电磁法监测煤矿含水采空区J.物探与化探,2 0 1 2(3 6):1 1 4-1 1 8.G u o W e i l i,X u o G u o q i a n g,Z h o u N a n n a n,e t a l.M o n i t o r i n g w a-t e r-f i l l e d g o a f b y u e i n g T E M t e c h n o l o g y.w i t h t r a n s m i s s i o n e l e c-t r o n m i c r o s c o p eJ.G e o p h y s i c a l a n d G e o c h e m i c a l E x p l o r a t i o n,2 0 1 2(3 6):1 1 4-1 1 8.1 1 陈实,李延清,李同贺,等.天然源面波技术在乌鲁木齐城市地质调查中的应用J.物探与化探,2 0 1 9,4 3(6):1 3 8 9-1 3 9 8.C h e n S h i,L i Y a n q i n g,L i T o n g h e,e t a l.T h e a p p l i c a t i o n o f n a t u-r a l s o u r c e s u r f a c e w a v e t e c h n o l o g y t o u r b a n