大定源回線(xiàn)TEM法在煤礦采空區(qū)的探測(cè)效果分析

摘 要：大定源回線(xiàn)TEM法在煤礦采空區(qū)起至關(guān)重要的作用，當(dāng)利用偶極子合成回線(xiàn)源方式，計(jì)算大定源回線(xiàn)TEM全區(qū)視電阻率處理實(shí)測(cè)采空區(qū)TEM數(shù)據(jù)時(shí)，可以有效消除晚期視電阻率存在的“邊界效應(yīng)”，得到的晚期視電阻率不存在視電阻率假高值，接近地表真實(shí)電阻率，這對(duì)探測(cè)地下煤礦采空區(qū)有重要意義。本文利用大定源回線(xiàn)TEM法對(duì)山西某煤礦探測(cè)成果進(jìn)行分析，根據(jù)視電阻率的相對(duì)高低分析采空異常區(qū)的分布規(guī)律，圈定低阻異常區(qū)域，推斷采空區(qū)范圍，并對(duì)視電阻率斷面圖和順層切片圖進(jìn)行共同解譯。

關(guān)鍵詞：瞬變電磁法；大定源回線(xiàn)；采空區(qū)；視電阻率

中圖分類(lèi)號(hào)：TD 82" " 文獻(xiàn)標(biāo)志碼：A

煤礦采空區(qū)是指在煤礦的作業(yè)過(guò)程中，將地下的煤炭或煤矸石等開(kāi)采完后留下的空洞或空腔。在煤礦開(kāi)采或者水利工程、公路工程、鐵路工程等的施工過(guò)程中，采空區(qū)常給正常施工帶來(lái)危害。因此，對(duì)地下采空區(qū)的分布范圍和空間形態(tài)特征進(jìn)行評(píng)判，對(duì)潛在危害性進(jìn)行評(píng)價(jià)以及合理確定采空區(qū)治理對(duì)策一直是困擾工程技術(shù)人員的問(wèn)題。目前國(guó)內(nèi)用于煤田采空區(qū)勘探的地球物理方法較多，主要有瞬變電磁法、高密度電法、人工地震法、地質(zhì)雷達(dá)法等[1-2]。在煤礦采空區(qū)的勘查中，不同地層具有不同的巖性及電性特征。在采空區(qū)地段，由于采空使周邊巖層結(jié)構(gòu)發(fā)生改變，因此該部位的視電阻率與圍巖會(huì)有一定的差異。采空區(qū)空洞對(duì)電流有排斥作用，當(dāng)?shù)孛嫣綔y(cè)時(shí)會(huì)表現(xiàn)相對(duì)高阻的特性，但當(dāng)采空區(qū)塌陷充水后，會(huì)表現(xiàn)相對(duì)低阻。煤礦采空區(qū)在塌陷前，通常很難預(yù)知，而以電性差異為基礎(chǔ)的大定源回線(xiàn)瞬變電磁探測(cè)技術(shù)具有勘探深度深、效應(yīng)高、可以穿透高阻屏蔽層、橫向分辨率高等特點(diǎn)，對(duì)煤礦采空區(qū)反應(yīng)靈敏，在煤礦采空區(qū)的探測(cè)中效果良好[3]。

1 TEM原理

瞬變電磁法又叫脈沖瞬變法，簡(jiǎn)稱(chēng)TEM，是一種時(shí)間域的電磁探測(cè)方法。利用不接地回線(xiàn)或接地電偶源向地下發(fā)送一次場(chǎng)，在一次場(chǎng)的間歇期觀(guān)測(cè)二次漩渦場(chǎng)的方法。其基本工作方法是在地面布設(shè)通以一定波形電流的發(fā)射線(xiàn)圈，從而在其周?chē)臻g產(chǎn)生一次磁場(chǎng)，并在地下導(dǎo)電巖礦體中產(chǎn)生感應(yīng)電流，斷電后感應(yīng)電流由于熱損耗而隨時(shí)間衰減。衰減過(guò)程分為早、中和晚期。早期的電磁場(chǎng)相當(dāng)于頻率域中的高頻成分，衰減快且趨膚深度淺。而晚期成分則相當(dāng)于頻率域中的低頻成分，衰減慢，趨附深度較深。通過(guò)測(cè)量斷電后各時(shí)間段的二次場(chǎng)隨時(shí)間變化規(guī)律，可以得到不同的地電特征。由于該方法是在一次場(chǎng)斷電的瞬間測(cè)純二次場(chǎng)，具有體積效應(yīng)小、橫向分辨率高、對(duì)低阻反應(yīng)靈敏、與探測(cè)地質(zhì)體具有最佳耦合和受旁側(cè)地質(zhì)體影響小的特點(diǎn)，因此可根據(jù)二次場(chǎng)衰減曲線(xiàn)特征判斷地下地質(zhì)體的電性、規(guī)模和產(chǎn)狀等[3-4]。

2 全區(qū)視電阻率定義

視電阻率是電法勘探中的重要參數(shù)之一，而全區(qū)視電阻率能夠在早期趨于第一層視電阻率真值，晚期趨于最后一層視電阻率真值，過(guò)渡期也能夠很客觀(guān)地反映電性的變化規(guī)律，因此研究全區(qū)視電阻率具有重要意義。

本文采用的是框內(nèi)全區(qū)視電阻率計(jì)算方法。首先，沿邊界將線(xiàn)框剖分為多個(gè)能視為水平電偶極子的電流元，其次，對(duì)均勻半空間電偶極子瞬變場(chǎng)及視電阻率進(jìn)行定義。最后，用偶極子疊加法即可計(jì)算大定源回線(xiàn)在地面一點(diǎn)的時(shí)域垂直磁場(chǎng)。

對(duì)水平電偶極子發(fā)射，均勻半空間地面任意一點(diǎn)（xc，yc）的時(shí)域垂直磁場(chǎng)如公式（1）所示。

（1）

對(duì)回線(xiàn)的某一條邊，均勻半空間地面任意一點(diǎn)（xc，yc）的時(shí)域垂直磁場(chǎng)如公式（2）所示。

（2）

式中：I為供電電流；dL為水平電偶極子的長(zhǎng)度；erf（u）為誤差函數(shù)。對(duì)整個(gè)回線(xiàn)激發(fā)，均勻半空間地面任意一點(diǎn)（xc，yc）的時(shí)域垂直磁場(chǎng)如公式（3）所示。

（3）

式中：HZLi為某個(gè)電偶極子產(chǎn)生的垂直磁場(chǎng)。在公式（1）～公式（3）中，垂直磁場(chǎng)與均勻半空間的視電阻率是隱函數(shù)關(guān)系，為不失一般性，將回線(xiàn)源的瞬變電磁場(chǎng)公式表示為H=φ（ρ，t，r），給定初值ρ0，在ρ0的鄰域內(nèi)對(duì)?（ρ，t，r）進(jìn)行泰勒展開(kāi)，如公式（4）所示。

（4）

省略后面的高階無(wú)窮小量，只保留前兩項(xiàng)，經(jīng)過(guò)整理可得視電阻率的近似表達(dá)式如公式（5）所示。

（5）

ρτ為視電阻率，并將公式（5）寫(xiě)成迭代格式如公式（6）、公式（7）所示。

ρτ（i+1）=?ρτi+ρτi，i=0，1，2，...n （6）

（7）

反復(fù)迭代下去，直到滿(mǎn)足公式（8）。

，其中，ξ[10-6，106]" （8）

由此可以得到任意一點(diǎn)的全區(qū)視電阻率。

式中：，r為收發(fā)距；x為測(cè)點(diǎn)橫坐標(biāo)；y為測(cè)點(diǎn)縱坐標(biāo)；?為垂直磁場(chǎng)H的相關(guān)函數(shù)；ρ為地層平均電阻率；t為時(shí)間。

3 應(yīng)用實(shí)例解釋分析

3.1 勘探區(qū)地層特征

勘探區(qū)煤系地層基底為奧陶系中統(tǒng)灰?guī)r，出露及鉆孔揭露地層由老到新依次為奧陶系中統(tǒng)、石炭系、二疊系以及第四系，將地層由老到新敘述如下。奧陶系中統(tǒng)峰峰組（O2f）為灰黃色-灰白色-深灰色，巨厚層狀，以石灰?guī)r為主，中夾泥質(zhì)灰?guī)r、角礫狀灰?guī)r、白云質(zhì)灰?guī)r，電性表現(xiàn)為高阻，晶粒結(jié)構(gòu)、泥晶結(jié)構(gòu)，溶蝕裂隙發(fā)育。與上覆煤系地層呈平行不整合接觸，厚度為139m。下伏奧陶系中統(tǒng)上馬家溝組，鉆孔未揭穿本組。石炭系中統(tǒng)本溪組（C2b）為底部黃綠色粉砂巖，黑灰色泥巖，中部為灰白色鋁質(zhì)泥巖，上部為灰色巨厚層狀石灰?guī)r，電性表現(xiàn)為中高阻，本溪組厚10～15m。石炭系上統(tǒng)太原組（C3t）為井田主含煤地層之一，由深灰色泥巖、石灰?guī)r、灰黑色砂質(zhì)泥巖、細(xì)砂巖、粉砂巖和煤組成，電性表現(xiàn)為中低阻。含煤6層，其中15#煤層穩(wěn)定，與下伏地層呈整合接觸，厚度110～115m。二疊系下統(tǒng)山西組（P1s）為井田主要含煤地層之一，由細(xì)粒砂巖、含泥粉砂巖、砂質(zhì)泥巖、泥巖及煤組成，電性表現(xiàn)為中低阻。含煤6層，其中3#煤層穩(wěn)定可采，3#煤厚4.5～5.7m，山西組厚52～59.5m。與下伏地層呈整合接觸。二疊系下統(tǒng)下石盒子組（P1x）的下部由黃綠色、杏黃色、灰黃色的頁(yè)巖及砂質(zhì)頁(yè)巖，夾中－細(xì)粒杏黃色、黃綠色長(zhǎng)石石英砂巖組成，上部由黃綠色、杏黃、黃色、雜色、紫色泥巖、頁(yè)巖夾黃色—杏黃色中細(xì)粒長(zhǎng)石石英砂巖組成，電性表現(xiàn)為中低阻。本組75～78m，與下伏地層呈整合接觸。二疊系上統(tǒng)上石盒子組（P2s）區(qū)內(nèi)出露廣泛，主要由砂巖、砂質(zhì)泥巖、泥巖、粉砂巖，含鋁泥巖互層組成，電性表現(xiàn)為中低阻。根據(jù)鉆孔揭露本組厚110～170m，最大厚度約295m。第四系（Q）的棕紅色亞黏土（含鈣質(zhì)結(jié)核），黃土，砂質(zhì)黏土，電性表現(xiàn)為低阻。不整合于各時(shí)代的地層上，主要分布于山梁及溝谷中，厚度一般為2～10m。

3.2 采集瞬變電磁數(shù)據(jù)

區(qū)內(nèi)共布置37條測(cè)線(xiàn)，測(cè)線(xiàn)方向均按NE28.7°布置，線(xiàn)距40m，點(diǎn)距20m。線(xiàn)號(hào)為147～ 267 ，幾條線(xiàn)的點(diǎn)號(hào)范圍不同。野外工作點(diǎn)位示意如圖1所示。采用瞬變電磁儀，工程測(cè)量采用GPS系統(tǒng)。經(jīng)過(guò)實(shí)地試驗(yàn)后，確定工作參數(shù)。大回線(xiàn)框?yàn)?00m×400m，發(fā)射頻率為25Hz，疊加次數(shù)為256，采樣窗口時(shí)間為0.036ms～27.915ms，延時(shí)為130μs，測(cè)道數(shù)為40道，發(fā)射電流約為3.5A，接收線(xiàn)圈等效面積為2200㎡。各項(xiàng)質(zhì)量檢查結(jié)果符合相關(guān)規(guī)范及設(shè)計(jì)要求，數(shù)據(jù)采集質(zhì)量可靠。

3.3 資料處理

首先對(duì)原始實(shí)測(cè)的數(shù)據(jù)進(jìn)行整理和預(yù)處理，其次計(jì)算全區(qū)視電阻率，經(jīng)過(guò)時(shí)深轉(zhuǎn)換、地形校正、濾波、白化和反演后，繪制全區(qū)視電阻率斷面圖及平面圖切片圖，最后進(jìn)行綜合解釋分析。

3.4 瞬變電磁視電阻率斷面圖解釋

為更好地分析TEM斷面異常以及測(cè)線(xiàn)間的相關(guān)性，本文共選取相鄰的兩條測(cè)線(xiàn)，對(duì)兩條測(cè)線(xiàn)的相關(guān)性以及每條測(cè)線(xiàn)的斷面異常進(jìn)行分析。圖2和圖3分別為勘探區(qū)199測(cè)線(xiàn)、203測(cè)線(xiàn)全區(qū)視電阻率斷面圖。圖中橫向長(zhǎng)虛線(xiàn)分別為3#、9#和15#煤層底板高程，下部的實(shí)線(xiàn)矩形框?yàn)橐阎?#煤層采空區(qū)范圍。虛線(xiàn)矩形框?yàn)?#和15#煤層推斷采空區(qū)，豎直粗實(shí)線(xiàn)為鉆孔。

根據(jù)甲方提供的資料，9#煤已基本全部采空，從兩幅斷面圖上可以看出3# 采空區(qū)的解釋結(jié)果與甲方提供的資料吻合，煤層在視電阻率斷面圖中反映的位置在中阻向高阻過(guò)渡的區(qū)域。3#煤層比9#和15#煤層的高程要高，并且基本全部采空，在斷面圖上的異常特征為低阻向中阻的過(guò)渡帶，形態(tài)上表現(xiàn)為等值線(xiàn)寬范圍、緩慢向下彎曲，原因是采空區(qū)塌陷所形成的電阻率變化帶較煤層的影響范圍大。由于工區(qū)內(nèi)煤層埋深在160m～360m，因此采空區(qū)解譯應(yīng)以低阻為主要特征，在形態(tài)上表現(xiàn)為向下彎曲的局部異常，范圍的大小與采空區(qū)的大小成正變關(guān)系。與周邊未采空的區(qū)域有較明顯差別。

圖2和圖3視電阻率斷面圖中提供了9#煤的已知采空區(qū)的范圍，兩幅圖都為200點(diǎn)～360點(diǎn)、400點(diǎn)～560點(diǎn)和600點(diǎn)～760點(diǎn)。局部有向下凹陷的低阻異常，如199線(xiàn)的200點(diǎn)～280點(diǎn)，330點(diǎn)～540點(diǎn)，720點(diǎn)～800點(diǎn)，900點(diǎn)～1010點(diǎn)，1070點(diǎn)～1270點(diǎn)，例如203線(xiàn)下部200點(diǎn)～280點(diǎn)形成一個(gè)低阻范圍，與旁邊的背景形成明顯的低阻局部異常，推斷該異常為采空區(qū)所引起。同樣在480點(diǎn)～700點(diǎn)及1140點(diǎn)～1300點(diǎn)間也形成了一個(gè)低阻范圍，在480點(diǎn)附近形態(tài)有向上彎曲的趨勢(shì)，應(yīng)該為未采段，也與甲方提供的采空區(qū)范圍相對(duì)應(yīng)（460點(diǎn)～500點(diǎn)段沒(méi)有采空區(qū)），而低阻異常范圍比甲方提供的采空區(qū)范圍小。該試驗(yàn)線(xiàn)煤層底板低，因此水滲流到采空區(qū)后，應(yīng)向該位置集中，其應(yīng)以低阻異常為特征，實(shí)測(cè)結(jié)果證實(shí)了該推斷。根據(jù)9#、15#煤采空區(qū)異常特征，推斷上述低阻異常為采空區(qū)所引起。從已知采空區(qū)和解釋采空區(qū)的位置來(lái)看，雖然二者有些小差別，但是總體的對(duì)應(yīng)性較好。9#煤在斷面圖中的異常特征為中阻向中高阻的過(guò)渡帶，形態(tài)上表現(xiàn)為等值線(xiàn)較3#煤小范圍、向下彎曲較明顯。與周邊未采空的區(qū)域有明顯差別。15#煤在斷面圖中的異常特征表現(xiàn)為中阻向中高阻的過(guò)渡帶，形態(tài)上表現(xiàn)為等值線(xiàn)較9#煤更小范圍、向下彎曲較明顯。與周邊未采空的區(qū)域有明顯差別。

從縱向上看，視電阻率從淺到深表現(xiàn)為中高-低-中高-高的電性特征（HA型），第一層反映了近地表黃土、二疊系上下石盒子地層。第二層反映了二疊系的下石盒子及山西組地層。第三層反映了石炭系的太原組及本溪組地層。第四層主要反映的是奧陶系的地層。從橫向上看，在調(diào)查采空區(qū)上方，由于采煤引起的上覆巖層的物性變化，因此視電阻率表現(xiàn)橫向的分帶性，具有明顯的視電阻率等值線(xiàn)變化的痕跡，圖中表現(xiàn)為以低阻為異常特征。從整體上看，這兩條測(cè)線(xiàn)的相關(guān)性比較好，推斷采空區(qū)和已知采空區(qū)對(duì)應(yīng)關(guān)系較好，異常比較可靠。

3.5 瞬變電磁斷面及平面綜合解釋

地下的煤層被采出后，在巖層內(nèi)形成一個(gè)有一定規(guī)模的空間，隨著時(shí)間的推移，上覆巖層就會(huì)在地球重力的作用下逐漸斷裂、塌陷，地下水就會(huì)侵入，采空區(qū)的視電阻率隨著發(fā)生變化，會(huì)出現(xiàn)相對(duì)低阻異常區(qū)。平面等值線(xiàn)圖是選取斷面圖上底板等高線(xiàn)上下各10m內(nèi)視電阻率值取平均，把這些斷面圖上選取的視電阻率做成平面等值線(xiàn)圖，其反映了煤系地層的視電阻率變化情況，通過(guò)視電阻率的相對(duì)變化，可以對(duì)煤系地層的采空區(qū)進(jìn)行推斷。為解煤層（平均厚10m）橫向范圍內(nèi)的電性情況，順層對(duì)典型的9#煤的底板平面視電阻率平面等值線(xiàn)圖（圖4）進(jìn)行分析，并根據(jù)前面所述原則劃分采空區(qū)范圍（橫軸為線(xiàn)號(hào)，縱軸為點(diǎn)號(hào)），平面等值線(xiàn)圖與視電阻率斷面圖較好，異常較可靠。

低阻異常區(qū)是采空區(qū)的主要特征，結(jié)合所收集的地質(zhì)資料結(jié)果，在此勘探區(qū)中總共圈定了7個(gè)采空區(qū)，采空區(qū)分布在工區(qū)南部、東南部和東北部和西北部，由南向北，由東向西依次編號(hào)為9c-1、9c-2、9c-3、9c-4、9c-5、9c-6和9c-7。因此推測(cè)勘探區(qū)的南部、西北部、東北部均有大片的采空區(qū)，在勘探區(qū)的中部存在一處較小的采空區(qū)域，應(yīng)引起重視，并及時(shí)打鉆驗(yàn)證。

4 結(jié)論

基于低阻異常區(qū)是探測(cè)采空區(qū)主要特征的地球物理前提，199測(cè)線(xiàn)和203測(cè)線(xiàn)視電阻率斷面圖的處理結(jié)果表明，斷面圖與實(shí)際的地電斷面對(duì)應(yīng)很好，可以做進(jìn)一步精細(xì)解釋。順層抽取9#煤底板的平面視電阻率平面等值線(xiàn)圖，根據(jù)視電阻率的相對(duì)高低分析采空異常區(qū)的分布規(guī)律，對(duì)斷面圖和平面圖進(jìn)行共同解釋。

文中根據(jù)低阻異常區(qū)是采空區(qū)的主要特征對(duì)采空區(qū)進(jìn)行劃分，但視電阻率的影響因素眾多，再加上地質(zhì)因素的復(fù)雜性和地球物理解釋的多解性，推測(cè)結(jié)果可能與實(shí)際情況有偏差，需要用鉆探及其他物探方法進(jìn)行驗(yàn)證。

參考文獻(xiàn)

[1]張俊英，王翰鋒，張彬，等.煤礦采空區(qū)勘查與安全隱患綜合治理技術(shù)[J].煤炭科學(xué)技術(shù)，2013，41（10）：76-80.

[2]董杰，孟慶生，尹明泉，等.綜合物探方法在采空塌陷區(qū)調(diào)查中的應(yīng)用[J].物探與化探，2013，37（3）：557-560.

[3]楊建軍，吳漢寧，馮兵，等.煤礦采空區(qū)探測(cè)效果研究[J].煤田地質(zhì)與勘探，2006，34（1）：67-70．

[4]張開(kāi)元，韓自豪，周韜.瞬變電磁法在探測(cè)煤礦采空區(qū)中的應(yīng)用[J].工程地球物理報(bào)，2007，4（4）：341-344.