多种物探方法在煤矿回采工作面防治中的应用.pdf

第3 2 卷增刊12023年6 月文章编号：10 0 4-4 0 5 1(2 0 2 3)S1-0355-06中国矿业CHINA MINING MAGAZINED0I:10.12075/j.issn.1004-4051.20230166Vol.32,Suppl 1June2023多种物探方法在煤矿回采工作面防治中的应用杨增林，杜平（国能神东煤炭地测公司，内蒙古鄂尔多斯0 17 0 0 0）摘要：煤矿回采工作面在回采过程中，顶底板的水害、回采煤层及回采煤层顶底板岩层的地质构造等容易造成瓦斯突出、水害等严重事故。为了煤矿的安全高效生产，必须探测回采工作面顶板的含水性和顶板以及回采工作面内部的地质构造情况。以神东某矿12 4 0 7 工作面为例，利用音频电穿透法、坑道无线电波透视法和瞬变电磁法三种工作面物探方法对顶板和工作面内部进行探测。通过音频电穿透法和瞬变电磁法相结合探查回采工作面顶板富水性以及地质构造情况；通过坑道无线电波透视法和瞬变电磁法相结合探查回采工作面内部富水性以及地质构造情况。三种物探方法探测成果进行相互对比分析、验证和总结，有效地探查了12 4 0 7 工作面顶板和工作面内部的富水性以及地质构造情况。三种物探方法施工简单、施工效率高，探测成果相互补充，有效地圈定了地质异常体的位置和范围，为煤矿安全高效生产、灾害防治提供可靠的依据。关键词：音频电穿透法；瞬变电磁法；坑道无线电波透视法；富水性；地质构造；灾害防治中图分类号：TD166Application of various geophysical methods in the prevention and(Geological Survey Company,Shendong Coal,CHN Energy,Ordos 017000,China)Abstract:During the mining process of coal mining working face,the water damage of roof and floor,thegeological structure of the coal seam and the rock formation of the roof and floor of the mining coal seam arelikely to cause serious accidents such as gas outburst and water damage.For the safe and efficient productionof coal mines,it is necessary to detect the water content of the roof of the working face and the geologicalstructure of the roof and the interior of the working face.Taking the 12407 working face of a mine inShendong as an example,the roof and the interior of the working face are detected by three working facegeophysical exploration methods:audio frequency electric penetration method,tunnel radio wave perspectivemethod and transient electromagnetic method.The water-rich roof and geological structure of the miningface are explored by the combination of audio frequency electric penetration method and transientelectromagnetic method;geological structure.The detection results of the three geophysical methods arecompared,analyzed,verified and summarized,and the water richness and geological structure of the roof andinterior of the 12407 working face are effectively probed.The three geophysical exploration methods havesimple construction,high construction efficiency,and complementary detection results,effectively delineatingthe location and scope of geological anomalies,providing a reliable basis for safe and efficient production anddisaster prevention in coal mines.Keywords:audio frequency electric penetration method;transient electromagnetic method;tunnel radiowave perspective method;water-rich;geological structure;disaster prevention and control文献标识码：Acontrol of coal mining faceYANG Zenglin,DU Ping收稿日期：2 0 2 3-0 3-2 0引用格式：杨增林，杜平.多种物探方法在煤矿回采工作面防治中的应用 J.中国矿业，2 0 2 3，3 2（S1）：3 5 5-3 6 0.YANG Zenglin,DU Ping.Application of various geophysical methods in the prevention and control of coal mining faceJJ.China MiningMagazine,2023,32(S1):355-360.责任编辑：刘硕3560 引 言在煤矿回采过程中，矿井水害和复杂的地质构造都会影响煤矿安全高效生产，因此，煤矿的防治水工作和探查复杂的地质构造情况的工作相当重要 1-3 。煤矿防治水细则规定：严格执行井下探放水“两探”要求，采掘工作面超前探放水应当同时采用钻探、物探两种方法，做到相互验证，查清采掘工作面及采空区区域、周边含水层富水性、以及积水等情况4。目前，井下回采工作面常用的物探方法有音频电穿透法和坑道无线电波透视法 5 。音频电穿透法在探查工作面顶底板富水性有一定的效果，但是由于井下的干扰因素较多，探测效果会有所下降 6-7 ；坑道无线电波透视法在探查工作面内部的地质构造情况有一定的效果，但是由于回采工作面宽度逐渐增大，无线电波要穿透工作面必须降低频率，但是降低频率以后其分辨率会下降，探测效果随之减弱 8 10 1。瞬变电磁法在矿井物探中的应用越来越广泛，理论研究也飞速发展。牛之琏等 11-14 对瞬变电磁的理论、数值模拟计算等进行了大量研究，并且在煤田地质勘探、采空区探查中进行针对性应用，取得了很好的效果；李松峰等 15 利用瞬变电磁法和无线电波透视法探查了煤矿水害情况；刘平喜等 16 1运用瞬变电磁法在井下富水采空区进行了探查应用，达到了很好的效果。根据常用的两种工作面物探方法音频电穿透和坑道无线电波透视法进行探测，再利用瞬变电磁法进行探测，以音频电穿透和坑道无线电波透视法为主，瞬变电磁法为辅，三种物探方法探测成果彼此相互对比、分析、借鉴，得出更精确的探测结果。以神东某矿12 4 0 7 工作面为例，运用三种物探方法探测工作顶板上方基岩含水层的富水性及工作面内部地质构造分布情况，为下一步防治水与工作面回采工作提供参考依据。1方法原理1.1音频电穿透法由于井下各种岩（矿)石导电性不同，进而影响人工电场的分布形态。煤系地层在正常组合条件下,在纵向与横向上的物理特性都有固定的变化规律 17 。当构造发育时，一般会出现层位错动和裂隙充水。裂隙充水类似于局部出现有明显的含水构造，井下裂隙水的导电性良好，会造成含水构造与围岩在电性上表现出较大差异，从而在横向与纵向上都打破了原有地层电性的变化规律。电法探测技术是研究以导电性差异变化特征为基础的，矿井音频中国矿业电穿透法就是利用专门的仪器在井下观测井下人工场源的分布规律和变化特征来达到解决地质问题的目的，矿井音频电穿透法仍属矿井直流电法 18-19 ，但因其施工方法、资料处理技术的差别及主要针对性（探测回采工作面内部含水性异常）等和直流电法的不同，从而形成矿井音频电穿透法分支。1.2坑道无线电波透视法坑道无线电波透视法又称坑透法，电磁波在地下岩层中传播时，由于各岩矿石的介电常数和电阻率不同，使得其对电磁波的能量吸收不同，低阻岩层对电磁波吸收较强，高阻岩层对电磁波吸收较弱2 0 1。当电磁波在地下岩层传播时遇见构造（如断层、陷落柱等)会在其断面产生反射和折射作用，造成电磁波能量的损耗，故通过探测电磁波的能量损耗可以推断对应的高、低阻岩层以及复杂的地质构造情况。1.3瞬变电磁法瞬变电磁法以岩石的导电性差异为基础，利用接地或不接地回线向探测目标体发送脉冲电流作为场源，激励探测目的物感生二次电流，在脉冲间隙测量二次场随时间的响应，从而了解探测目标体介质的电性变化情况 2 1。在电流断开之前，发射电流在回线周围的空间中建立起一个稳定的磁场。将电流突然断开，由该电流产生的磁场也立即消失，一次磁场的这一剧烈变化激发出感应二次场。由于感应二次场的衰减规律与地下地质体的导电性有关，导电性越好，二次场的衰减越慢；导电性越差，二次场衰减越快。因此，通过分析研究二次场的衰减规律即可达到探测地下地质异常体的目的 2 2 。2工程量2.1音频电穿透法工程量12407工作面，工作面长2 5 5 4.2 m，工作面宽300.2m（图1）。音频电穿透法物探工程起点为工作面切眼，终点为工作面回撤通道，总施工长度2554.2m。运输顺槽、回风顺槽均进行施工，运输顺槽布置测量点2 5 6 个，回风顺槽布置测量点2 5 6 个。采用f=120Hz和f=15Hz双频点施工，探查工作面顶板0 6 0 m、6 0 12 0 m 深度层段含水层的富水异常区的平面位置、形态及含水性的相对强弱。2.2坑道无线电波透视法工程量坑道无线电波透视法施工是从切眼到回撤通道方向进行，探测长度为2 5 5 4.2 m。发射点距5 0 m，接收点距10 m，采集供电点对应前后10 0 m范围内的数据，对应巷道的一定区段进行扇形扫描接收，探查工作面内断层、陷落柱及煤层变薄带的延展范围。第3 2 卷增刊12.3眼瞬变电磁法工程量瞬变电磁法施工共有四条测线：第一条测线是在回风顺槽施工，从切眼到回撤通道，一共有2 5 6 个测点，探测方向为正邦顶板（发射线圈平面和顶板平面夹角3 0）；第二条测线是在回风顺槽施工，从切眼到回撤通道，一共有2 5 6 个测点，探测方向为正邦顶板（发射线圈平面和顶板平面夹角6 0）；第三条测线是在回风顺槽施工，从切眼到回撤通道，一共有256个测点，探测方向为正邦（发射线圈平面和顶板平面夹角9 0）；第四条测线是在主运顺槽施工，从切眼到回撤通道，一共有2 5 6 个测点，探测方向为正邦（发射线圈平面和顶板平面夹角9 0）。3物探成果分析3.1音频电穿透资料解释根据层析成像原理，利用供电点到测量点的电压降，重建回采面的电性变化，通过多次迭代求出电导率近似解。通过计算成图得出12 4 0 7 工作面顶板300-号异常200-10000Fig.2 Results of 0-60 m audio frequency electric penetration of the roof of 12407 working face in a coal mine300-200常10000Fig.3 Results of 60-120 m audio frequency electric penetration of the roof of 12407 working face in a coal mine电磁波场强衰减系数300-PDC5200DC1100QS600杨增林，等：多种物探方法在煤矿回采工作面防治中的应用某煤矿12 4 0 7 工作面图1某煤矿12 4 0 7 工作面示意图Fig.1 Schematic diagram of 12407 working face of a coal mine12408回风顺槽?号异常?号味常5001000图2 某煤矿12 4 0 7 工作面顶板0 6 0 m音频电穿透成果图号异常号常号异常5001000图3 某煤矿12 4 0 7 工作面顶板6 0 12 0 m音频电穿透成果图DC25001000图4 某煤矿12 4 0 7 工作面无线电波透视成果图Fig.41Radio wave perspective results of 12407 working face in a coal mine35706 0 m 音频电穿透成果图（图2），12 4 0 7 工作面顶板6 0 12 0 m音频电穿透成果图（图3）。通过图2可以圈定6 个异常区，位置相对集中,其中，号异常区、号异常区异常幅值中等，号异常区、号异常区、号异常区、号异常区异常范围较大且幅值相对较强。通过图3 可以圈定10 个异常，位置较为分散，其中，号异常区号异常区异常幅值中等，号异常区异常幅值相对较强。3.2坑道无线电波透视资料解释无线电波透视相对衰减层析成像技术对数据进行处理，相对衰减层析成像是用相对衰减算法进行层析处理，以电磁波场强衰减系数在工作面中的相对分布进行异常图像的重建，得出12 4 0 7 工作面无线电波透视成果图（图4）。通过图4 可以看出7 处明显的相对衰减较强异常区，分别编号为DC1号异常、DC2号异常、DC3号异常、DC4号异常、DC5号异常、QS6号异常、QS7号异常。12407主运顺槽12407回风顺槽5号异常?号常150012406运输顺槽12408回风顺槽?号异常号昇常号常150012406运输顺槽2 0 0 012408回风顺槽DC3DC4150012406运输顺槽2 0 0 05.2号常3.73.22.72.0002.500岩石电导率色标5.20号异常2.725000.0.008.800.004QS7.0038.0020.001250000.0013583.3瞬变电磁资料解释对采集的数据进行去噪处理，根据晚期场公式将电流归一化值转换成不同的时间窗口所对应的视电阻率值，然后进行时深转换处理，得到各个测线视电阻率断面图。第一条测线的瞬变成果图为12 4 0 7100十500100+500300-292001000Fig.7 Transient electromagnetic results of 12407 return air and main transport trough roof at 90 in a coal mine通过图5 可以看出有5 个低阻异常区，分别为I异常区、异常区、异常区、IV异常区、V异常区,其中，I异常区、异常区、V异常区、V异常范围大，幅值强；号异常范围相对较小，幅值中等。通过图6 可以看出有5 个低阻异常，分别为I异常区、异常区、异常区、IV异常区、V异常区,其中，I异常区、异常区、IV异常范围大，幅值强；异常区、V号异常范围相对较小，幅值中等。通过图7 可以看出有11个低阻异常区,分别为1号异常区11号异常，其中，1号异常区、2 号异常区、4 号异常区、5号异常区、6 号异常区、11号异常区异常范围较大，幅值较强；3 号异常区、7 号异常区、8 号异常区、9号异常区、10 号异常区异常范围较大，但是幅值相对较小。3.4综合异常分析12407工作面末采段位于母花海子古冲沟中上游，回采范围内上覆基岩厚7 3 13 1m，第四系松散层厚6 4 2 m，松散含水层厚0 3 2 m。上覆基岩厚度最小处位于12 4 0 7 回顺15 联巷左右，对应松散含水层厚度最大。12 4 0 7 回顺靠近12 4 0 6 采空区，中国矿业回风顺槽顶板3 0 瞬变电磁成果图（图5）；第二条测线的瞬变成果图为12 4 0 7 回风顺槽顶板6 0 瞬变电磁成果图（图6）；第三条测线和第四条测线的瞬变成果图为12 4 0 7 回风顺槽、主运顺槽顶板9 0 瞬变电磁成果图（图7)。200号异常区号异常区号异常区V号异常区5001000图5 某煤矿12 4 0 7 回风顺槽顶板3 0 瞬变电磁成果图Fig.5 Transient electromagnetic results at 30 on the roof of 12407 return air in a coal mine1号异常区号异常区500Fig.6 Transient electromagnetic results at 60 of 12407 return air along the roof of a coal mine500图7 某煤矿12 4 0 7 回风顺槽、主运顺槽顶板9 0 瞬变电磁成果图第3 2 卷岩层视电阻率V昌异常区501500200012406运顺亚号异常区IV号异常区10001500图6 某煤矿12 4 0 7 回风顺槽顶板6 0 瞬变电磁成果图12408回风顺槽某矿12 4 0 7 工作面1000150012406运输顺槽2 0 0 0局部低洼点可能有积水。因此，图2 的6 个异常区可能与顶板上方基岩含水层含水性有关，图3 的10个异常区可能与顶板上方上覆松散层富水性有关。图5 和图6 中的低阻异常位置大概重合，说明两次瞬变电磁探测的异常位置基本在同一方向，图5 瞬变电磁成果图中的I号异常基本和图2 工作面顶板060m音频电穿透成果图中的号异常、图3 工作面顶板6 0 12 0 m音频电穿透成果图中的号异常、号异常、号异常、号异常基本吻合。图5瞬变电磁成果图中的号异常基本和图2 工作面顶板0 6 0 m音频电穿透成果图中的号异常、图3工作面顶板6 0 12 0 m音频电穿透成果图中的号异常、号异常基本吻合。图5 瞬变电磁成果图中的IV号异常基本和图2 工作面顶板0 6 0 m音频电穿透成果图中的号异常、图3 工作面顶板6 0 120m音频电穿透成果图中的号异常、号异常基本吻合。图5 瞬变电磁成果图中的V号异常基本和图2 工作面顶板0 6 0 m音频电穿透成果图中的号异常、号异常、图3 工作面顶板6 0 12 0 m音频电穿透成果图中的号异常基本吻合。2500V号异常区岩层视电阻率200150100502.00012406运顺2.500三10岩层视电南200250010增刊1根据图4 坑透成果图并结合已知两顺槽地质面图和地质资料，判断DC1号异常区可能是已揭露断层引起的异常;DC2号异常区可能为已揭露断层向工作面内部延伸引起的异常;DC3号异常区可能是已揭露断层引起的异常;DC4号异常区可能是已揭露断层向工作面内部延伸引起的异常；DC5号异常区可能为已揭露断层向工作面内部延伸引起的异常；QS6号异常区范围较大，幅值相对较强，致使异常区范围内电磁波数值整体变化明显，初步推测异常区可能为冲刷构造延伸引起的异常；QS7号异常区范围较大，幅值相对较强，异常区在靠近回风顺槽侧平行分布，异常区范围内电磁波数值整体变化明显，初步推测异常区可能为巷道起伏石厚度变化引起的异常。在图7 瞬变电磁成果图中1号低阻异常、5 号低阻异常、6 号低阻异常整体视电阻率较低基本和QS6号异常吻合，3 号低阻异常、8 号低阻异常基本和DC3号异常吻合，4 号低阻异常基本和DC5号异常吻合，7 号低阻异常基本和DC2号异常吻合，9 号低阻异常、10 号低阻异常、11号异常基本和QS7号异常吻合。4验证情况通过物探划分的顶板上方低阻富水异常区，进行了打钻疏放，图2 中号异常、号异常、号异常、号异常、号异常所在位置都打钻出水，涌水量较大。在回采过程中DC2号异常、DC3号异常、DC5号异常、QS6号异常、QS7号异常基本和推断的地质情况相吻合。5结语矿井回采工作面仅靠一种物探方法探测，数据不会太准确，如果利用多种物探方法对探测的数据进行互相对比、分析和验证，可以的得出更精确的探测结果。对回采工作面的顶板探测以音频电穿透法为主，瞬变电磁为辅可以有效、精确的探查出富水异常区；对回采工作面内部探测以坑道无线电波透视法为主，瞬变电磁法为辅可以有效探查工作面内的地质构造。参考文献1刘清宝，丁湘，冯洁，等.“八位一体”闭环式煤矿顶板水害管控模式的构建与应用 J.煤田地质与勘探，2 0 2 1，4 9（4）：17 0-17 7.LIU Qingbao,DING Xiang,FENG Jie,et al.Construction andapplication of eight in one closed-loop management and con-trol mode of roof water disaster in coal minesJJ.Coal Geolo-gy&Exploration,2021,49(4):170-177.2 马良，郭瑞.准格尔煤田不连沟矿井构造特征及其对煤矿生产的影响 J.煤田地质与勘探，2 0 2 0 4 8（3）：3 5-4 4,5 0.杨增林，等：多种物探方法在煤矿回采工作面防治中的应用gy,2020,48(S2);24-28.4 刘守强，武强，曾一凡.煤矿防治水细则修订要点解析 .煤炭工程，2 0 19,5 1（3)：1-4LIU Shouqiang,WU Qiang,ZENG Yifan.Analysis of the re-vision points of Detailed Rules for Water Control in CoalMineJJ.Coal Engineering,2019,51(3):1-4.5张国恩，解振华，史洪恺.音频电穿透技术在探查工作面顶板富水性的应用 J.煤矿安全，2 0 2 15 2（5）：8 2-8 6.ZHANG Guoen,XIE Zhenhua,SHI Hongkai.Application ofaudio frequency electric penetrating technology in exploringwater-rich roof of working faceJl.Safety in Coal Mines,2021,52(5):82-86.6 王鹏飞，代凤强，王根盛，等.音频电穿透技术在呼吉尔特矿区的应用 J.煤炭技术，2 0 15，3 4（7）：12 4-12 6.WANG Pengfei,DAI Fengqiang,WANG Gensheng,et al.Application of voice frequency electric perspective techniquein Hjert Coal Mine AreaJJ.Coal Technology,2015,34(7):124-126.7 高靖学，高，闵鹏.音频电透视技术在顶板富水性探测中的应用 J.陕西煤炭，2 0 2 1（3）：9 2-9 6.GAO Jingxue,GAO Tuo,MIN Peng.Application of audio fre-quency electric perspective technology in roof water abun-dance detectionJ.Shaanxi Coal,2021(3):92-96.8 温亨聪，杨海涛，张辉辉.多频无线电坑透技术在煤矿地质构造探测中的应用 J.煤炭技术，2 0 2 1，4 0（7）：7 1-7 4.WEN Hengcong,YANG Haitao,ZHANG Huihui.Applica-tion of multi-frequency radio wave perspective technology ingeological structure exploration of coal mineJJ.Coal Tech-nology,2021,40(7):71-74.9 温亨聪，李学臣，杨海涛，等.无线电坑透技术探查矿井老旧空巷效果研究 J.煤炭技术，2 0 2 1，4 0（4）：6 3-6 5.WEN Hengcong,LI Xuechen,YANG Haitao,et al.Study oneffect of radio penetrating technology in exploration of oldempty roadway in mineJJ.Coal Technology,2021,40(4):63-65.10 朱立.无线电坑透在潞宁矿2 2 116 工作面应用分析JJ.山东煤炭科技，2 0 2 0 3 8（12）：16 7-16 9.ZHU Li.Application analysis of radio pit penetration in the22116 working face of Luning MineJJ.Shandong Coal Sci-ence and Technology,2020,38(12):167-169.11牛之琏.瞬变电磁测深法对导电层的探测能力 J.地质与勘探,19 9 2(1):4 2-4 7.NIU Zhilian.Detectability of an intermediate conductive layerby transient electromagnetic soundingJJ.Geology and Explo-ration,1992(1):42-47.12 牛之琏.瞬变电磁测深曲线微分参数解释方法 J.中南矿冶359MA Liang,GUO Rui.The structural characteristics and theirinfluence on production of Buliangou coal mineJ.Coal Geol-ogy&Exploration,2020,48(3):35-44,50.3 靳玉强，杨建飞，王迎，等.煤层火烧区边界“两探”措施的探索及应用J.煤炭科学技术，2 0 2 0，4 8（S2）：2 4-2 8.JIN Yuqiang,YANG Jianfei,WANG Ying,et al.Explorationand application of“two explorations measures for the bound-ary of coal seam burning areaJl.Coal Science and Technolo-360学院学报，19 9 2（2）：13 0-13 5.NIU Zhilian.The interpretive method of transint electromag-netic sounding by curve differential parameterJJ.Journal ofCentral South University of Mining and Metallurgy,1992(2):130-135.13牛之琏，李爱华.瞬变电磁法在铅锌金矿床上的找矿效果 J.地质与勘探，19 9 1,19(11)：4 1-4 3.NIU Zhilian,LI Aihua.Exploration results of TEM over a Pb-Zn-Au deposit in Hunan ProvinceJ.Coal Geology&.Explo-ration,1991,19(11):41-43.14 牛之琏.近区瞬变电磁测深曲线畸变段的消除方法 J.煤田地质与勘探，19 9 2,2 0（1)：6 0-6 3.NIU Zhilian.Eliminate method for distorted stage of transientelectromagnetic sounding curve in near-zoneJJ.Coal Geology&.Exploration,1992,20(1):60-63.15李松峰，王栋，魏军军.综合物探在煤矿水害防治中的应用J.华北自然资源，2 0 2 1（6）：15-18.LI Songfeng,WANG Dong,WEI Junjun.The application ofcomprehensive geophysical prospecting in the prevention andcontrol of water hazards in coal minesJ.Huabei Natural Re-sources,2021(6):15-18.16刘平喜，徐贵阳，曹伟康，等.瞬变电磁法在井下富水采空区探测中的应用 J.煤炭与化工，2 0 2 0，4 3（10）：4 9-5 1.LIU Pingxi,XU Guiyang,CAO Weikang,et al.Application oftransient electromagnetic detection in underground water-richair goafJJ.Coal and Chemical Industry,2020,43(10):49-51.中国矿业17张志巍.含隐伏断层煤层底板高承压水导升及突水机制研究D.北京：煤炭科学研究总院，2 0 2 1.18李静蕊，解海军，孙展雄.富水构造地质异常体矿井直流电法超前探查响应特征 J.能源与环保，2 0 2 1，4 3（5）：5 1-5 6.LI Jingrui,XIE Haijun,SUN Zhanxiong.Response character-istics of direct current method advanced exploration in minewith water rich structural geological anomalyLJJ.China Ener-gy and Environmental Protection,2021,43(5):51-56.19张平松，欧元超，李圣林，我国矿井物探技术及装备的发展现状与思考 J.煤炭科学技术，2 0 2 1，4 9（7)：10-15.ZHANG Pingsong,OU Yuanchao,LI Shenglin.Developmentquo-status and thinking of mine geophysical prospecting tech-nology and equipment in ChinaJJ.Coal Science and Technol-0gy,2021,49(7):10-15.20屈小兵，鲍俊睿.基于电磁波层析成像技术的隐伏构造区边界判定 J.煤矿安全，2 0 2 0,5 1（4）：14 7-15 1.QU Xiaobing,BAO Junrui.Boundary determination of hiddenstructural areas based on electromagnetic wave tomographyJ.Safety in Coal Mines,2020,51(4):147-151.21李貅.瞬变电磁测深的理论与应用 M.西安：陕西科学技术出版社，2 0 0 2:10-2 5.22薛国强，李貅，底青云.瞬变电磁法正反演问题研究进展 .地球物理学进展，2 0 0 8,2 3（4）：116 5-117 2.XUE Guoqiang,LI Xiu,DI Qingyun.Research progress inTEM forward modeling and inversion calculationJJ.Progressin Geophysics,2008,23(4):1165-1172.第3 2 卷