李彥星

（山西省煤炭地質(zhì)物探測繪院，山西晉中030600）

瞬變電磁及氡值測量法探測煤層積水采空區(qū)應(yīng)用與研究

李彥星*

（山西省煤炭地質(zhì)物探測繪院，山西晉中030600）

基于瞬變電磁法和氡值測量法的基礎(chǔ)理論，通過分析2種方法的特點及局限性，提出綜合2種探測方法應(yīng)用于煤礦采空區(qū)探測。結(jié)合煤礦積水采空區(qū)探測實例，通過采集、處理、成果推斷解釋等各流程的對比探討得出結(jié)論。結(jié)果表明利用綜合地球物理探測方法，可以有效彌補單一方法受干擾等問題出現(xiàn)，造成假異常導(dǎo)致解釋成果精度降低的缺陷，提高了勘探精度，為煤礦安全生產(chǎn)提供保障。

瞬變電磁法；氡值測量法；積水采空區(qū)；干擾；安全生產(chǎn)

瞬變電磁法具有穿透能力強、干擾小、施工效率高、探測深度大等很多優(yōu)點，廣泛應(yīng)用于金屬礦產(chǎn)探測、水文地質(zhì)勘查及煤層采空區(qū)探測。氡值測量法也應(yīng)用于煤層采空區(qū)探測和煤層自燃區(qū)探測。以上2種方法在實際工作中均取得了很好的探測效果。但2種方法在實際應(yīng)用中均存在一定局限性。如今城鎮(zhèn)化、工業(yè)化進程不斷加快，地面瞬變電磁數(shù)據(jù)采集時受工頻及諧波干擾嚴(yán)重［1］。目前，各種去噪方法消除工頻干擾的效果并不理想［2］。瞬變電磁法探測高阻體反映不靈敏，且會受“邊框效應(yīng)”影響導(dǎo)致數(shù)據(jù)質(zhì)量變差［3］。淺層積水采空區(qū)附加效應(yīng)對瞬變電磁信號影響也很大［4］。而氡值測量法在一個測點只有一個測定值，沒有垂向的分辨能力［5］?；谝陨戏治?，考慮綜合應(yīng)用2種方法對采空區(qū)進行探測。這樣對照分析，相互驗證可以提高資料解釋精度。

1　瞬變電磁與氡值測量法基本原理

1.1瞬變電磁法基本原理

使用不接地回線通以脈沖電流作為場源，以激勵探測目的物感生二次電流，在脈沖間隙測量二次場隨時間變化的響應(yīng)。由于二次場從產(chǎn)生到結(jié)束的時間是短暫的，即是“瞬變的”，所以稱之為瞬變電磁法，其屬于時間域電磁法。

瞬變電磁法的物理基礎(chǔ)為電磁感應(yīng)原理，即導(dǎo)電介質(zhì)在階躍變化的激勵磁場激發(fā)下產(chǎn)生渦流場的問題。對于定源回線裝置，瞬變電磁場表達式為：

式中：q——接收線圈有效面積；

R——發(fā)送回線半徑，當(dāng)發(fā)送回線為矩形時R=L/π；

a——地下導(dǎo)體的半徑；

h——導(dǎo)體的中心埋深；

τ——時間常數(shù)，k=2πR/τ。

1.2氡值測量法基本原理

煤層采空區(qū)形成后，促使地質(zhì)體發(fā)生變形，改變了地下地質(zhì)體的應(yīng)力分布狀態(tài)，從而使地下液體氣體的運移與集聚環(huán)境發(fā)生改變，對氡氣的運移與富集具有一定的控制作用。首先煤層采空在地下形成較大空洞，成為地下水與氣體的儲存場所；其次采空區(qū)煤層頂板塌陷冒落，地下應(yīng)力場改變，采空區(qū)及周邊應(yīng)力減小，地下氣體及液體自然由高應(yīng)力區(qū)域向低應(yīng)力區(qū)域運移，造成氡氣在采空區(qū)的聚集；最后在采空區(qū)冒落帶及其派生裂隙形成后，氣體自下向上運移，在地溫與地壓作用下，氡氣必然與其他氣體（CO、CO2、CH4、H2S等）一起自地下深處向地表遷移，在地表形成與采空區(qū)域平面形態(tài)基本相同的氡異常區(qū)。氡為非極性單原子分子，而活性炭為非極性吸附劑。當(dāng)氡運移到活性炭表面時，這2種物質(zhì)的分子或原子相互接近時，在色散力的主要作用下，氡則很快被活性炭吸附。在電離室測量氡及其子體衰變輻射出的γ射線強度，就可測量提取氡值數(shù)據(jù)。

2　工作方法

2.1瞬變電磁工作方法

瞬變電磁法施工裝置類型較多［6］，根據(jù)探測目標(biāo)體的不同，在施工時視具體情況選擇相適應(yīng)的裝置。裝置確定后通過現(xiàn)場試驗確定發(fā)射和采集參數(shù)。國內(nèi)常用的工作裝置主要有重疊回線、定源回線和偶極裝置。本次擬選定源回線裝置施工。工程測線布設(shè)時垂直于主要構(gòu)造走向，按工程要求比例尺度設(shè)計測網(wǎng)密度［7］。施工時接收機按測網(wǎng)設(shè)計測量點位依次順點采集。

2.2氡值測量工作方法

在每個測量點的位置挖40cm深的倒錐形坑，將裝有活性炭的瓶子編號后埋入坑中，在地表留記號便于取杯時找到埋置地點，工作筆記上記錄該點埋杯時間等相關(guān)信息。埋置的活性炭瓶5d之后統(tǒng)一取出，記錄取杯時間。送回實驗室，放入活性炭測氡儀電離室中測量氡及其子體輻射出的γ射線強度，然后對測氡數(shù)據(jù)進行提取，對原始數(shù)據(jù)進行預(yù)處理后成圖解釋。

3　工程應(yīng)用

3.1勘查地區(qū)地質(zhì)概況

施工地區(qū)位于山西省西北部，地層由老至新分別是奧陶系、石炭系、二疊系、第四系，區(qū)內(nèi)地表基本被第四系地層所覆蓋，基巖僅在溝谷中出露。主要可采6、8號煤層，為保證煤礦安全生產(chǎn)，開展了地面綜合地球物理探測研究。

3.2地球物理特征

瞬變電磁勘探是從測量地表電場的分布狀況來分析地下電性情況的［8］，地下異常體的存在和分布必須使地表電場有明顯變化才能測出。本次探查區(qū)域地層由老至新在縱向上各層電性差異較大，而橫向上未開采煤層與積水采空區(qū)的電性差異同樣較大［9］。由此可見，積水采空區(qū)的視電阻率明顯區(qū)別于完整煤層的視電阻率，是瞬變電磁法探測積水采空區(qū)并分析積水情況的物理前提。

測氡是基于地層中含有豐富的天然放射性元素，其中又以鈾的同位素所占比例最大，鈾經(jīng)一系列衰變后形成氡，氡及氡的衰變母體鐳沿著構(gòu)造帶、裂隙和地下水的垂向搬運在地表富集，形成氡異常。放射性元素隨水中的SiO2含量增加而增加，且地溫的升高加快了氡向地表的遷移。本區(qū)勘查積水采空區(qū)，應(yīng)用氡值測量法探測是可行的。

3.3數(shù)據(jù)采集

本次勘查設(shè)計測網(wǎng)為測線距40m，測點距20m，根據(jù)探測目的層布設(shè)320m×320m發(fā)射線框。瞬變電磁使用目前較為主流的加拿大鳳凰V8電法儀。通過試驗確定發(fā)射頻率10Hz，發(fā)射電流6A，使用等效100m2專用探頭，600次疊加接收。氡值測量儀器采用活性炭測氡儀，此方法測網(wǎng)與瞬變電磁測網(wǎng)一致，這樣便于成果的對比分析。

3.4資料處理與解釋

3.4.1瞬變電磁正反演數(shù)值模擬

針對本地區(qū)正常地層地電條件及采空區(qū)地電條件，建立相應(yīng)地電模型進行瞬變電磁正反演數(shù)值模擬。旨在經(jīng)過數(shù)值模擬得到地下采空區(qū)的瞬變電磁響應(yīng)，結(jié)合理論模型分析瞬變電磁法探測煤層采空區(qū)的可行性。

正演：通過解頻率域亥姆霍茲方程求得諧變場，通過域的轉(zhuǎn)換技術(shù)將頻率域轉(zhuǎn)換到時間域，求得晚期視電阻率。并將計算結(jié)果通過與解析解的對照，驗證其精確性。

假設(shè)3層層狀地電模型，如圖1-a所示：第一層為蓋層厚度為100m，電阻率100Ω·m；第二層為煤層采空區(qū)厚度為20m，電阻率分別為5（Model 1）、20（Model 2）、80（Model 3）、400Ω·m（Model 4），分別代表采空強積水、采空積水、未采空及采空未積水4種情況；第三層為底層，電阻率為150Ω·m。模型參數(shù)：發(fā)射線框400m×400m，發(fā)射電流1A，采樣時間為10-6～10-1s，指數(shù)等分50個時刻。

通過正演計算結(jié)果對比，圖1-b所示：瞬變電磁探測低阻靈敏度高，但電阻率越低對深部影響越大，有屏蔽現(xiàn)象。如圖1中Model 1米字間隔曲線所示。對高阻體即采空未積水，如Model 4三角間隔曲線與黑色均勻半空間曲線對比后，發(fā)現(xiàn)此方法對高阻體反映不靈敏。

反演：利用上文地電模型正演數(shù)據(jù)，使用一維Occam反演法計算，進行反演結(jié)果分析，得出瞬變電磁法探測煤層采空區(qū)的可行性結(jié)論。

從圖2中看出2-a反演結(jié)果清晰，反演異常位置與給定模型基本對應(yīng)，但是低阻層位下部200m以內(nèi)數(shù)據(jù)均受上部低阻體影響，導(dǎo)致低阻體下部數(shù)據(jù)無法真實反演；圖2-b看出反演異常結(jié)果位置準(zhǔn)確，與圖2-a比較圖2-b采空煤層電阻率較大，所以對下部屏蔽不如圖2-a嚴(yán)重；圖2-c與圖2-d均可看出瞬變電磁對高阻體反演不靈敏，其反演結(jié)果幾乎沒有高阻反映，無法準(zhǔn)確識別。

圖1　瞬變電磁采空區(qū)正演結(jié)果對比圖

圖2　瞬變電磁采空區(qū)反演結(jié)果對比圖

通過針對本地區(qū)的數(shù)值模擬，得出瞬變電磁在本地區(qū)探測的理論結(jié)論：①可以有效探測呈低阻的煤層積水采空區(qū)；②無法準(zhǔn)確探測采空未積水的高阻區(qū)域；③若上部有低阻積水采空區(qū)，則與其間距較近的下部地質(zhì)體無法真實反演（即受條件所限瞬變電磁分層探測能力有限）。

3.4.2成果解釋

瞬變電磁視電阻率擬斷面分析。從已知資料得知在，鉆孔ZK-1302（對應(yīng)本次勘探工程布置為730線1030樁號）發(fā)現(xiàn)6號煤缺失，推斷鉆孔周圍是煤層采空區(qū)域。從瞬變電磁730測線視電阻率擬斷面圖3看出此測線1000附近，剖面反映為低阻區(qū)域，根據(jù)電阻率變化，可準(zhǔn)確圈定低阻范圍。從剖面看出瞬變數(shù)據(jù)成層性較好，且反演結(jié)果與鉆孔資料從縱向到橫向基本吻合。瞬變電磁測線700～900、1200～1300附近推斷為正常區(qū)域。

圖3　瞬變電磁視電阻率擬斷面圖與鉆孔對比圖

氡值剖面曲線分析。如圖4，490線由西向東布置，共76個點。由上圖分析在350點處有個高值，而附近變化沒規(guī)律，所以10～1090、1150～1390點推斷為正常區(qū)，1090～1150、1390～1570點氡值起伏比較大，有規(guī)律性，氡值普遍偏高，推斷為采空區(qū)。

圖4　490線氡值剖面曲線圖

瞬變電磁及氡值平面疊合分析。通過測線斷面分析校正，進入平面分析階段。本次勘查使用的2種方法數(shù)據(jù)處理，成果均按照數(shù)理統(tǒng)計法結(jié)合已知資料，通過人機交互來確定異常閾值。從統(tǒng)計概率角度，在具體計算算術(shù)平均值和標(biāo)準(zhǔn)偏差時，應(yīng)該注意剔除極大、極小的值的影響，保證數(shù)據(jù)呈正態(tài)分布。使用處理后數(shù)據(jù)體生成斷面及平面圖，以視電阻率低于閾值200Ω·m劃為瞬變電磁異常區(qū)；氡值測定值則將高于閾值620個計數(shù)/3min區(qū)域圈為氡值異常區(qū)。

從圖5-a分析氡值異常，顏色由淺至深代表氡值由高至低。整體上北高南低，按氡值異常分析煤層采空區(qū)，則采空主要處于勘探區(qū)域北部，南部未見異常。從5-b看出瞬變電磁順煤層切片平面異常基本與氡值異常相符，只是在勘探區(qū)南部瞬變電磁存在低阻異常。對照已知資料并結(jié)合實際采集情況分析，在實地采集數(shù)據(jù)時，勘探區(qū)南部有高壓線經(jīng)過，因此瞬變電磁數(shù)據(jù)受干擾，從而導(dǎo)致反演結(jié)果顯示南部低阻異常，推斷此處為假異常。又由于前文針對本測區(qū)所做數(shù)值模擬得出結(jié)論，本次瞬變電磁探測6號、8號煤層間隔較近，反演結(jié)果無法準(zhǔn)確分層，因此無法分層推斷各煤層低阻異常，只能將2組煤層合并為一層，按合并后順層切片推斷分析低阻異常。

3.5探測成果

圖5　氡值異常及視電阻率平面圖

結(jié)合已知水文地質(zhì)資料及其他相關(guān)資料，并通過上文分析，得出本次綜合物探勘探成果。本次勘探積水采空區(qū)位于勘探區(qū)北部，2種探測方法可相互印證；勘探區(qū)南部推斷不存在積水采空區(qū)，因為氡值未顯示異常，同時瞬變電磁在此處異常可信度較差。本區(qū)共劃出積水采空區(qū)3處，礦方針對本次物探成果提前做出施工部署，且北中部圈定積水異常已得到驗證。

4　結(jié)論

綜合瞬變電磁法與氡值測量法應(yīng)用于煤層積水采空區(qū)探測，可以彌補瞬變電磁在施工時遇到各種無法避免干擾時，造成數(shù)據(jù)質(zhì)量下降解釋成果精度降低的缺陷。又可以2種方法相互驗證，排除假異常。這樣可以從根本上提高地球物理勘探精度，達到優(yōu)勢互補。應(yīng)用結(jié)果表明：此技術(shù)可有效圈定富水異常區(qū)域，推斷積水采空區(qū)，指導(dǎo)礦方安全生產(chǎn)。同時可推廣應(yīng)用至其他勘探領(lǐng)域。

［1］牛之璉.時間域電磁法原理［M］.長沙：中南大學(xué)出版社,2007：37-103.

［2］曹建章，陳民新，宋建平，等.瞬變電磁測量中抗工頻干擾問題的研究［J］.工程勘查，2000（6）：35-37.

［3］解海軍.煤礦積水采空區(qū)瞬變電磁法探測技術(shù)研究（博士論文）［D］.中國地質(zhì)大學(xué),2009.

［4］解海軍，孟小紅，王信文，等.煤礦積水采空區(qū)瞬變電磁法探測的附加效應(yīng)［J］.煤田地質(zhì)與勘探，2009,37（2）：71-74.

［5］楊建軍，申燕，劉鴻福.測氡法和瞬變電磁法在探測煤礦采空區(qū)的應(yīng)用［J］.物探與化探，2008,32（6）：661-664.

［6］劉國興.電法勘探原理與方法［M］.北京：地質(zhì)出版社,2005：192-211.

［7］齊文秀，牛之鏈.DZ/T0187-1997地面瞬變電磁法技術(shù)規(guī)程［S］.中華人民共和國地質(zhì)礦產(chǎn)部，1997.

［8］朱若軍，程久龍，張要田，等.礦井工作面開采水文地質(zhì)條件綜合地球物理探測研究［J］.中國礦業(yè)，2010,19（7）：98-101.

［9］張曉峰.瞬變電磁法探測煤田采空區(qū)的應(yīng)用研究（碩士論文）［D］.長安大學(xué),2007.

［10］李世峰，梁盛軍，周建雄.強干擾背景條件下TEM探測鐵礦資料解釋方法研究［J］.地質(zhì)與勘探，2006,42（3）：72-75.

［11］王信文，馮宏，解海軍，等.TEM技術(shù)地形校正初探［C］//陜西省地球物理文集（八）——深部礦井災(zāi)害源探測實踐，2008：127-131.

P631.3

A

1004-5716（2016）01-0153-05

2015-01-13

2015-01-14

李彥星（1980-），男（漢族），山西長治人，工程師，現(xiàn)從事地球物理勘探、地質(zhì)災(zāi)害調(diào)查工作。