復雜地形下瞬變電磁法探測多層采空區積水

李章波，陳　露

復雜地形下瞬變電磁法探測多層采空區積水

李章波1，陳露2

（1.中煤科工集團重慶研究院，重慶 400039， 2.西南石油大學 地球科學與技術學院，成都610500）

根據井田構造形態勘探線呈NE-SW向布置，瞬變電磁法測量采用40m×20m的網度400×400m發射線框，15A發射電流，8.333Hz發射頻率，420μ s采樣延遲時間等參數能夠取得最佳數據效果。資料處理時采用直接剔除畸變數據和非線性濾波兩種方法消除噪聲干擾，處理結果為5號煤層含水異常區有14個，其中2個富水性強，5個富水性中等；8號煤層含水異常區有10個，其中1個富水性強；1個富水性中等；富水性中等以上的區域面積為113 281m2。經礦方實際掘進證實了研究成果的準確性。

瞬變電磁；多目標層；老空水；技術參數

瞬變電磁法（TEM）是利用不接地回線向地下發射階躍波形電磁脈沖，激發地下導體產生感應渦旋電流，在渦流電流激發的二次磁場沒有立即隨一次場的消失而消失時，及時接收二次磁場并研究其與時間的變化規律，從而確定地下目標體的電性分布結構及空間形態的地球物理勘探方法[1]-[5]。它具有探測深度大、體積效應小、異常形態簡單、分層能力強等優點，它對低阻充水斷層、充泥充水溶洞、破碎帶等不良地質體反應靈敏，易于探測出低弱的電阻率異常，可有效的區分采空區及含水區，因而在煤炭行業應用廣泛。

探測研究區地理位置為東經110°57′51″~110°59′56″，北緯37°41′03″~37°42′43″。為呂梁山系梁峁狀中低山黃土地貌，侵蝕切割強烈、沖溝發育，地形復雜、溝谷縱橫，不利于探測方法的實施和解譯。同時本次工作的目標層為5號、8號煤層采空區，埋深跨度大（180~520m），噪聲干擾強。因而需要對探測技術參數進行對比試驗和優化選取。

1　地質特征及其地球物理特性

1.1地形地貌

探測區屬呂梁山系，井田內大部被黃土所覆蓋，黃土層受強烈的侵蝕切割，形成為典型的梁峁狀中低山黃土地貌，地形復雜、溝谷縱橫，地勢總體呈現為北高南低，最高點位于井田西南部，標高+1 114.6m，最低點位于井田中部的溝谷，標高+885.0m，最大相對高差229.60m。

1.2地質特征

探測區地表大部區域為黃土覆蓋，基巖零星出露。井田內發育的地層由老至新有：奧陶系中統峰峰組灰巖、泥灰巖、泥巖、角礫狀灰巖；石炭系中統本溪組鐵鋁巖、灰巖、泥巖、砂巖；上統太原組石英砂巖、砂質泥巖、泥巖及煤層；二疊系下統山西組砂巖、砂質泥巖、泥巖及煤層；石盒子組中細粒砂巖、灰綠色泥巖、砂質泥巖；上第三系上新統粘土、亞粘土（含鈣質結核層）；第四系中上更新統黃土；全新統砂礫石。

5號煤層位于山西組，埋深約180～200m，平均厚度7.21m；8號、9號煤層位于太原組，埋深約280～320m，平均厚度分別為3.21m、3.78m，8號煤層上距5號煤層約47m。

1.3地層地球物理特性

一般而言，灰巖、煤層、礫石層、中粗砂巖、粉砂巖、泥巖電阻率依次降低，煤層的電阻率略高泥灰巖[2-5]。在橫向上，煤系地層導電性具有均一性，當煤層中出現采空區時，則表現為異常高阻，當采空區充水時，則表現為異常低阻；通過對比探測到的電阻率的變化、異常，可得到異常分布形態 ，進從而達到探測采空區和老空水的目的[2]-[5]。

探測區地層的電性特征表、圖1所示。

勘探區地層巖性、厚度及電阻率一覽表

2　工作參數研究

圖1　鉆孔柱狀圖與旁側測點視電阻率曲線對比圖

不同的探測區、不同的埋深、不同的地質背景、不同的目的層與探測目標層數需選擇不同的工作參數，合理的工作參數是完成探測目的的關鍵因素。本次選取45個試驗點進行參數研究，檢查物理點61個。

2.1發射線框規格

一般而言發射線框面積越大，探測深度越深，但縱向分辨率隨之降低，同時降低施工效率。本次工作目的層埋深180m～320m，故選擇邊長為200×200m和400 ×400m的發射線框進行優選。

采用同頻率25Hz、同發射電流15A、同采樣時間3min對邊長200m×200m和400m× 400m的發射邊框進行試驗，結果（圖2）表明兩種發射線框的起始有效數據點時間相同，400×400m線框的最后一個有效數據點時間比200×200m線框的晚十幾毫秒，說明前者對于晚期信號的干擾壓制效果好，能取得時間窗口更寬的原始數據，提高了勘探深度。因而應采用400×400m規格的發射線框。

2.2發射電流

圖2　發射線框對比試驗圖

圖3　發射電流對比試驗結果

一般而言發射電流越大，一次場的影響越深，激發的二次場的強度也越大、持續時間也越久、勘測深度也更深。在發射線框為400×400m，同發射頻率、同采樣時間為3min時，對10A和15A的發射電流進行比選，結果如圖3所示。當發射電流為10A時（圖3a），原始衰減曲線末端數據出現蹦跳，數據離差較大，5Hz晚期道甚至出現負值（空缺，曲線不顯示）；而15A時（圖3b）晚期信號光滑、連續?？梢?，發射電流為10A時激發的二次場強度不夠，不能有效的壓制晚期信號的干擾，因而選擇發射電流為15A。

2.3發射頻率

發射頻率影響充電時間和量測時間，發射頻率低，充電時間和量測時間長，能接收到更深的信號（晚期信號），但淺部分辨率低；發射頻率高，充電時間和量測時間短，晚期信號效果差，但淺部分辨率高。在同發射線框（400×400m）、同發射電流（15A）、同采樣時間（3min）的情況下，對25Hz、8.333Hz和5Hz三種發射頻率進行比選，結果如圖4所示。

由圖4可知，發射頻率為25Hz時在100m深度處信號較密集，對淺部探測分辨率較高；發射頻率為8.333Hz和5Hz時，200m～400m深度處信號較密集，能很好的探測中深部目的層，但發射頻率為5Hz晚期信號強度低于背景噪聲，抗干擾能力差。根據本次工作目的層的埋深情況（180m-320m），最終確定發射頻率為8.333Hz。

2.4采樣延遲研究

采樣延遲影響不同深度信號的接收，采樣延遲越大采集的深部信息越可靠，在同發射邊框（400m× 400m）、同發射電流（15A）、同發射頻率（8.333Hz）情況下對采樣延遲時間為380μ s、400μ s、420μ s、440μ s進行比選試驗，試驗結果如圖5所示，可見380μ s、400μ s的采樣延遲信號很不穩定，440μ s的采樣延遲信號衰減較快，420μ s的采樣延遲信號穩定，故選擇采樣延遲時間為420μ s。

圖4　發射頻率對比試驗結果

圖5　采樣延遲對比試驗曲線圖

3　現場工作

3.1儀器設備

現場儀器采用加拿大鳳凰地球物理公司研制的第八代多功能電法系統（V8），由發射系統、采集（接收）系統、定位系統、數據記錄處理系統組成；具有3個磁道、3個電道， 頻率范圍10 000Hz-0.000 05Hz （20000秒），工作溫度-20℃-+50℃；采用GPS同步，當GPS信號不好時，系統內晶振時鐘會自動啟動同步；采用24位A/D轉換器，可保持最高動態范圍和分辨率；可接收多頻點信號，提高了測量垂向分辨率和勘探精度；可高精度同步疊加、掃頻，任意增加疊加次數和掃頻時間；采用大功率發射機，頻率高，可有效的避免干擾信號。

3.2測線布置

探測區地質構造形態為一走向北東、向北西緩傾的單斜，褶曲不發育，地層走向變化不大，地層傾角為5°~8°。因此本次工作勘探線沿NE-SW向展布，采用40m×20m網度進行一次高密度面性量測，以提高對中小異常的分辨能力。

本次工作布置測線58條，剖面總長23 260m，測線物理點1 212個，另有檢查物理點61個，實驗點45個。

4　資料處理

圖6　資料處理流程圖

4.1技術措施與流程

探測區為黃土覆蓋的丘陵地貌，沖溝縱橫，交通繁忙，干擾噪聲大，因而采用直接剔除畸變數據和非線性濾波兩種方法進行數據預處理。最終采用TEMPRO及IX1D瞬變電磁法處理軟件進行數據處理，處理主要流程如圖6所示。

4.2資料解釋

資料解釋采取以下步驟：①解析多測道斷面圖和對應的視電阻率斷面圖；②將視電阻率、視深度做成一個瞬變電磁測深的三維數據體，根據從不同角度切取的視電阻率等值線平面圖進行富水異常的綜合解釋；③圈出含水異常帶。

1）多測道曲線：根據二次感應電壓表現為中間高兩側低的倒“U”形或倒“V”形特征解譯含水異常體，富水性越強，中間高與兩側低的變化幅度越大（圖7）。

2）等視電阻率斷面圖：當采空區或裂隙帶充水時，導電性加強，在電法剖面圖上表現為局部高電導、低電阻異常，視電阻率等值線上表現為“高-低-高”或圈閉狀低阻異常。

3）視電阻率順層水平切面圖：多測道曲線和等視電阻率斷面均是從測線剖面上進行的解釋，未考慮含水體的平面聯系，一般情況下，同一分布穩定、完整的層位其電阻率值變化不大，反映在視電阻率圖上即等值線分布均勻、平緩；反之，若局部地層為富水和含導水構造區域，電性的均勻分布規律被打破，在順層切片圖上異常區就會呈現視電阻率值減小、等值線扭曲、變形圈閉或呈密集條帶狀等。

4.3地質成果

5號、8號煤層探測結果如圖9所示（藍色為富水區），由圖9可知，5號含水異常區域14個，其中強含水異常區2個（深藍色）、中等含水異常區5個（淺綠色區域）、弱含水異常區7個（黑色圈定區域），強含水區面積22 563m2，中等含水區面積83 898m2，弱含水區面積69 535m2；8號含水異常區域10個，其中強含水異常區個（深藍色）、中等含水異常區個（淺綠色區域）、弱含水異常區8個（黑色圈定區域），強含水區面積3 587m2，中等含水區面積3 233m2，弱含水區面積37 005m2。

圖7　多測道曲線含水異常示意圖

圖8　含水異常的瞬變電磁等視電阻率斷面圖

5　結論

1）由于場地地形起伏大，溝壑縱橫，對比研究發現，發射線框邊長為400×400m即可有效的壓制晚期信號的干擾，再放大線框邊長已無必要且影響施工效率。

2）考慮探測深度的要求和信號干擾的因素，瞬變電磁供電電流為15A時探測效果較好。

3）本次勘探目的層深度180~320m，對比研究發現發射頻率為8.333Hz時對目標層探測的分辨率較高，抗干擾能力強。

4）為了采集到可靠的深部信息，采樣延遲時間選擇為420μ s。

5）工區溝壑縱橫，交通繁忙，干擾信號多且較強，因而在數據處理時采用直接剔除畸變數據和非線性濾波兩種方法。

6）通過實驗選擇的以上瞬變電磁法參數，在研究礦井取得了良好的地質效果。勘探結果表明，工區內5號煤層含水異常區有14個，其中2個富水性強，5個富水性中等；8號煤層含水異常區10個，其中1個富水性強，1個富水性中等。

7）成果驗證：在勘探研究報告提交半年后，據礦方反饋，在掘進8101運輸巷時，在8號煤層II號含水異常區外圍警戒線發現巖壁潮濕，打超前鉆，20m見老空（鉆桿推進速度明顯加快），22m見煤，25m見木屑卡鉆終孔，涌水量1.35m3每小時。表明本次勘探在研究礦井取得了良好的地質效果，有效的指導了礦井巷道開拓，為煤礦安全生產提供了參考資料。

圖9　煤順層切片立體圖（上圖5號煤層，下圖8號煤層）

[1] 陳載林，黃臨平，陳玉梁. 我國瞬變電磁法應用綜述[J]. 鈾礦地質，2010.26（1）：51～54

[2] 范濤，李文剛，王鵬 等. 瞬變電磁擬M T深度反演方法精細解釋煤礦巖層富水性研究[J]. 煤炭學報，2013.S1：129～135.

[3] 張占斌，孫世國， 董培鑫. 大定源瞬變電磁法探測撫順煤礦采空區的研究[J]. 北方工業大學學報，2013.25(1)：90～94.

[4] 張運霞，牛向東. 瞬變電磁勘探技術在深部礦井水害防治工作中的應用[J]. 礦業安全與環保，2013.40(3)：100～106.

[5] 劉云. 起伏地形大地電磁、時間域瞬變電磁二維數值模擬及直接反演法[D]. [博士學位論文]. 成都理工大學，2012

[6] 黃啟春. 瞬變電磁法在棗莊市采空區工程勘察中的應用[J]. 勘察科學與技術，2013.181（1）： 61～64

[7] 王凱, 劉鴻福. 瞬變電磁法運用 GDP-32Ⅱ探測煤層采空區的效果[J]. 工程勘察, 2012.2：92～95.

The Application of TEM to the Detection of Multilayer Goaf Water under Condition of Complex Landform

LI Zhang-bo1CHEN Lu2
（Chongqing Research Institute， China Coal Technology & Engineering Group， Chongqing400039；2-Institute of Earth Science and Technology， Southwest Petroleum University， Chengdu610500）

the feasibility tests and the optimal selection of emission frame size， emission current， emission frequency， sampling delay time， exploration line arrangement， and data processing techniques are taken in order to detect the deeply buried and multi-layer goaf water under condition of complex landform by the TEM. The results show that it is efficient that the emission frame size is 400×400m， the ?emission current is 15A， the emission frequency is 8.333Hz， the ?sampling delay time 420μs. The Exploration line was arranged NE-SW direction according to the structural form in the mine field. Detection results show that there are 14 abnormal aqueous areas in the fifth coal seam area， 2 among which are strong rich water areas and 5—medium rich water areas. There are 10 abnormal aqueous areas in the eighth coal seam area， one among which is strong rich water areas and one medium rich water areas. The total size of medium to strong rich water areas is 113281m2. The mine actual excavation confirms the accuracy of the research results.

complex landform； multi-target layer； TEM； goaf water； technical parameter

P631,3+25

A

1006-0995（2015）03-0443-05

10.3969/j.issn.1006-0995.2015.03.030

2014-05-29

李章波（1983-）男，四川內江人，工程師，主要從事煤礦防治水方面物探工作