牛小軍
河北省煤田地質(zhì)局物測地質(zhì)隊,河北 邢臺 054000
目前煤礦突水事故大多數(shù)是由于直接與附近老窯采空積水區(qū)打通或直接揭露隱伏的導(dǎo)含水構(gòu)造而引起瞬間大量突水所致,造成了較大的經(jīng)濟損失,同時對煤礦工人的人身安全也帶來了較大的威脅。因而,在巷道迎頭掘進或采煤工作面開始回采之間查明其掘進或開采影響范圍內(nèi)是否存在隱伏的導(dǎo)含水構(gòu)造或老窯采空積水區(qū),提前采取治理措施成為礦井防治水工作的關(guān)鍵所在。目前煤礦防治水領(lǐng)域的探測主要包括直接法和間接法兩類[1][2],直接法常采用鉆探實現(xiàn),該技術(shù)具有結(jié)果直觀、可靠等明顯優(yōu)點。間接法常采用礦井地球物理方法實現(xiàn),目前礦井地球物理探測方法主要有:礦井直流電法、瑞雷波法、礦井地質(zhì)雷達和礦井瞬變電磁法等。
本文研究礦井瞬變電磁法的基本理論與相關(guān)技術(shù),首先在理論上對礦井瞬變電磁場的全空間分布特征和水平層狀介質(zhì)中低阻異常體的瞬變電磁響應(yīng)特征進行了研究,在此基礎(chǔ)上設(shè)計了井下全空間條件下多方向超前探測技術(shù),該技術(shù)可快速有效的探測出掘進巷道正前方及其采動影響范圍內(nèi)是否存在賦含水構(gòu)造異?;蚶细G采空積水區(qū)的空間方位。
在導(dǎo)磁率為μ、導(dǎo)電率為σ的巷道迎頭上安放一個矩形天線,其面積為S,在天線中供以階躍脈沖電流[3][4]
在t<0時,也就是電流斷開之前,發(fā)射電流在回線空間中形成一個比較穩(wěn)定的磁場;如果將電流斷開(既時),回線空間中的一次場會隨即消失。磁場的變化會在圍巖中產(chǎn)生二次場(如圖1所示),在此過程存在的電磁場就是礦井瞬變電磁場。
圖1 TEM探測方向示意圖
礦井瞬變電磁超前探測一般在煤礦的采掘巷道內(nèi)進行,其瞬變電磁探測裝置接收回線中的感應(yīng)電位是采掘巷道周圍一定范圍內(nèi)所有介質(zhì)的綜合電性反映。所以,所采集的視電阻率值為迎頭范圍內(nèi)全部巖層介質(zhì)的綜合電性反映,視電阻率按照以下公式[6-7]來計算:
式中:C為響應(yīng)系數(shù);S是接收線圈的面積;N指線圈的匝數(shù);t是指二次場的衰減時間;V/I指的是接歸一化后二次場的電位值。在巷道的迎頭部位,通過線圈采集二次場隨時間變化的感應(yīng)電動勢,就可以推測迎頭附近介質(zhì)的電磁性質(zhì)在空間和時間上的變化規(guī)律,從而確定構(gòu)造異?;蚋凰惓5漠a(chǎn)狀和發(fā)育的位置等。
圖2 平行層理激發(fā)三層模型示意圖
由于煤系地層成層分布,各巖層電阻率不同,礦井瞬變電磁法進行定向超前探測時,發(fā)射線圈平面垂直于層理面(地層界面),線圈法線方向(即探測方向)平行層理(圖2所示),激發(fā)方向為平行層理。設(shè)計平行層理激發(fā)三層模型如圖2所示,由上至下第一層頂板為高阻層,電阻率為500Ω·m;第二層電阻率100Ω·m,厚度80m;第三層底板為低阻層,電阻率10Ω·m。激發(fā)源距第一層與第二層分界面、第二層與第三層分界面均為40m。圖3為利用2.5D時間域有限差分模擬算法得到的激發(fā)源位于第二層中間平行層理激發(fā)時三層地電斷面各時刻瞬變電場等值線的分布圖。由圖可知,在整個擴散衰減過程中,最強感應(yīng)電流密度中心始終位于激發(fā)源所在平面上,由于電性斷面的非對稱性,上下兩個感應(yīng)電流密度中心向外擴散速度不一致,高阻層一側(cè)的電流密度中心擴散速度較快,很快擴散至高阻層中;而低阻層一側(cè)的電流密度中心擴散速度較慢,衰減過程中,逐漸擴散至整個空間中,并最終占據(jù)主導(dǎo)地位。
圖3 平行層理激發(fā)三層地電斷面瞬變電場擴散圖
通過對圖2模型中各參數(shù)不同時的模擬計算分析可知,均勻全空間瞬變電磁場的擴散規(guī)律與均勻半空間不同,對于磁偶源激發(fā)的均勻全空間瞬變電磁場而言,最大電流密度區(qū)始終位于發(fā)射線圈所在平面上,瞬變電磁場以球形方式勻速向空間各個方向擴散,因而激發(fā)源附近接收線圈接收到的信號中不僅包含發(fā)射線圈前后介質(zhì)的電性信息,還包含發(fā)射線圈所在面方向上介質(zhì)的信息。平行層理激發(fā)條件下,最初時期瞬變電磁場僅分布于激發(fā)源所在層,低阻層對瞬變電磁場的分布尚未產(chǎn)生影響,隨著時間的推移,瞬變電磁場開始擴散到低阻層,低阻層中感應(yīng)的渦旋電流逐漸升高而后衰減,低阻層的衰減信息開始占據(jù)主導(dǎo)地位,源點處瞬變磁場響應(yīng)曲線出現(xiàn)符號翻轉(zhuǎn)現(xiàn)象,這種表現(xiàn)跟激發(fā)源所在層厚度以及所在層與低阻層的電性差異有關(guān),所在層厚度越小,電性差異越大,表現(xiàn)越明顯;早期瞬變電磁場的分布與低阻層厚度無關(guān),當(dāng)瞬變電磁場擴散至低阻層后,晚期階段源點處瞬變磁場時間變化率的響應(yīng)強度隨低阻層厚度的增加而略微減小;源點處瞬變磁場的響應(yīng)強度受地電斷面成層分布的影響隨激發(fā)源所在層厚度的增大而減小。瞬變電磁超前探測結(jié)果與激發(fā)源所在層的厚度、電阻率、源所在層外側(cè)電性層厚度、電阻率有關(guān)系,激發(fā)源所在層厚度與超前探測距離的比值越大,其探測效果越接近于傳統(tǒng)的垂直層理激發(fā)方式的探測效果,進而采用垂直層理激發(fā)方式下數(shù)據(jù)處理與解釋方法,結(jié)果與真實情況越貼切。
平行層理激發(fā)過程中,感應(yīng)電流密度中心分布傾向于低阻層,當(dāng)感應(yīng)電流密度中心移動到低阻層內(nèi)部時,瞬變電場開始擴散到低阻層外側(cè)電性層,并很快擴散至整個外層空間,與此同時,各電性層內(nèi)瞬變電場分別沿層理方向上擴散衰減,電流密度中心拉伸逐漸變長,激發(fā)源前后介質(zhì)逐漸影響瞬變電場的擴散,進而實現(xiàn)超前探測。
由于井下空間的限制,礦井瞬變電磁法的發(fā)射和接收線圈都不能使用太大的,只能使用匝束較多回線邊長較小的發(fā)射和接收線圈(邊長為2~3m)。在巷道迎頭布置測點(圖4),即在巷道迎頭從左側(cè)開始,首先使發(fā)射和接收線圈的法線方向與巷道的左側(cè)面垂直進行探測,然后按一定角度旋轉(zhuǎn)線圈,使線圈的法線方向和巷道的左側(cè)面形成15°、30°、45°、60°,75°和90°的角度進行測量,當(dāng)線圈法線的方向垂直迎頭面時,據(jù)迎頭面的實際寬度設(shè)置2~3個點,到迎頭右側(cè)時再使線圈旋轉(zhuǎn),使線圈的法線和巷道右面形成90°,75°、60°、45°、30°和15°的夾角實施測量。采用上述主法實現(xiàn)多角度數(shù)據(jù)采集,以獲得盡可能完整的前方介質(zhì)電性分布信息。
圖4 測點水平布置示意圖
在實際工作過程中對于圖4中的每個發(fā)射點,可調(diào)整天線的法線與巷道底板的夾角大小,以探測巷道頂板、順層和底板方向的圍巖變化情況,探測方法如圖5所示。這樣可得到位于巷道迎頭前方一個錐體范圍內(nèi)地層介質(zhì)的電性變化情況如圖6所示。
圖5 定向超前探測剖面方向示意圖
圖6 定向超前探測體積置示意圖
瞬變電磁超前探的資料解釋方法和過程首先是對所采集數(shù)據(jù)中的噪音去除處理,依據(jù)公式(2)將瞬變電磁儀所測的電流值歸一化后轉(zhuǎn)換成不同時間值所對應(yīng)的視電阻率的值,隨后再實施時間和深度轉(zhuǎn)換計算處理,得到的結(jié)果是每一條測線的視電阻率斷面圖。一般采用2種顯示成果的方法:矩形成圖方法是其中的一種,在圖5中的每一個測點都可以看作是等間距分布在一條直線上的點,在視電阻率等值線圖中測點標(biāo)號作為橫坐標(biāo)(或者是相對于00度測點的距離),把探測深度作為縱坐標(biāo);扇形成圖法是第2種顯示方法,在圖4中的每一個測點對應(yīng)的視電阻率值展布在實際平面上,以扇形模式繪制視電阻率等值線圖,把在迎頭位置向正前方測量的距離(m)作為縱坐標(biāo)、把圖5中兩側(cè)00測點的測量距離加上迎頭的寬度(m)作為橫坐標(biāo),其中迎頭的寬度可依據(jù)正前方測點的間距按照一定的比例放大。在定向超前探測結(jié)果中使用扇形方式顯示可提高超前探測資料的解釋精度,能更準(zhǔn)確反映出異常體的空間特征。
山西長治王莊煤業(yè)有限責(zé)任公司附近小煤窯比較多,現(xiàn)在經(jīng)過煤炭企業(yè)整合重組后小煤礦都已經(jīng)停產(chǎn)關(guān)閉。由于缺少附近小煤礦的生產(chǎn)相關(guān)圖件,特別是防治水方面的相關(guān)地質(zhì)資料,為整合煤礦的后繼生產(chǎn)帶來了潛在的水災(zāi)安全隱患。為了查清老窯采空積水區(qū)的具體位置及其影響范圍,山西長治王莊煤業(yè)有限責(zé)任公司采用礦井瞬變電磁定向超前探測技術(shù)在巷道掘進過程中進行跟蹤超前探測。
圖7為巷道掘進迎頭按圖4中方法所示沿煤層方向采用礦井瞬變電磁法定向超前探測的視電阻率等值線圖。從圖中可知巷道迎頭前方55~90米范圍內(nèi)存在相對低阻異常區(qū),巷道迎頭左前方和右前方也存在低阻異常區(qū)。從已掌握的相關(guān)水文地質(zhì)資料分析可知,該低阻異常區(qū)為老窯采空積水區(qū)的反應(yīng),其具體影響范圍見圖8中的陰影區(qū)域所示。
圖7 順層方向視電阻率等值線圖
圖8 超前探測成果示意圖
圖9 探測方向剖面示意圖
該異常區(qū)域經(jīng)礦方鉆探驗證,當(dāng)鉆探到巷道迎頭正前方約48米時鉆孔出水,水壓、水量較大,經(jīng)過多也鉆孔和擴大孔徑的方法歷時一個多月時間放水量約12萬m3。隨后巷道繼續(xù)向前掘進,揭露老窯廢棄巷道位置如圖8中黑色條帶區(qū)域所示。
冀中能源股份有限公司某掘進巷道前方存在一落差約30m的斷層。為了給巷道的安全掘進提供相關(guān)技術(shù)資料,礦方要求采用礦井瞬變電磁法定向超前探測技術(shù)進一步確定F7斷層的具體影響范圍。在實際工作過程中對于圖4中的每個發(fā)射點,調(diào)整天線的法線與巷道底板的夾角大?。?5°,0°,-15°和-30°四個方向依次探測(如圖9所示),以獲得巷道迎頭前方錐體范圍內(nèi)地層介質(zhì)的電性情況,并由此進一步研究相關(guān)構(gòu)造及含水性情況。
圖10至13分別為在巷道掘進迎頭處15°,0°,-15°和-30°四個方向探測的視電阻率等值線圖。從圖11中對應(yīng)的視電阻率等值線分析可知,在沿探測方向(順煤煤層)巷道迎頭正前方約80m處出現(xiàn)有相對明顯電性差異界面,根據(jù)已知地質(zhì)資料可解釋為探巷前方F7斷層的邊界反應(yīng)。從其對應(yīng)視電阻率值相對大小分析,可知F7斷層賦含水性較弱或不含水。從圖11、圖12和圖13中可知,視電阻率等值線的變化相對較為均勻,無明顯的低阻異常擾動,說明在探測范圍內(nèi)沒有明顯的賦含水地質(zhì)構(gòu)造異常。
根據(jù)本次探測結(jié)果,巷道繼續(xù)掘進過程中及時調(diào)整坡度方向,實現(xiàn)了斷層上下兩盤煤層的對接。圖14為實際探巷素描剖面示意圖,圖中顯示在探測迎頭前方約78m即為F17斷層邊界,與探測結(jié)果基本吻合。
圖10 迎頭15度方向視電阻率等值線圖
圖11 迎頭0度方向視電阻率等值線圖
圖12 迎頭-15度方向視電阻率等值線圖
圖13 迎頭-30度方向視電阻率等值線圖
礦井瞬變電磁法定向超前探測技術(shù)可快速查明巷道前方及其采動影響范圍內(nèi)不同方位的賦含水構(gòu)造、老窯采空積水區(qū)。該方法也可精確查明采煤工作面頂、底板采動影響范圍內(nèi)不同空間位置隱伏導(dǎo)含水通道或局部富含水區(qū)域的具體影響范圍,并對其富含水性進行定性評價,為煤礦的防治水工作提供相關(guān)技術(shù)資料,防止突水災(zāi)害事故的發(fā)生。
圖14 實際探巷素描剖面示意圖
[1]劉志新.礦井瞬變電磁場分布規(guī)律與應(yīng)用研究[D].徐州:中國礦業(yè)大學(xué),2008
[2]姜志海.巷道掘進工作面瞬變電磁超前探測機理與技術(shù)研究[D].徐州:中國礦業(yè)大學(xué),2008
[3]李貅.瞬變電磁測深的理論與應(yīng)用[M].西安:陜西科學(xué)技術(shù)出版社,2002
[4]蔣邦遠主編.實用近區(qū)磁源瞬變電磁法勘探[M].北京:地質(zhì)出版社,1998
[5]LIU Zhixin,YUE Jianhua,LIU Shucai.《Application of Synthetic Mine Geophysical Prospecting Technique in Detecting Collapse Columns》[J],International Conference“Waste Management,Environmental Geotechnology and Global Sustainable Development”,Ljubljana,SLOVENIA,28-30 August 2007
[6]劉志新,岳建華,劉仰光.扇形探測技術(shù)在超前探測中的應(yīng)用研究[J].中國礦業(yè)大學(xué)學(xué)報,2007,36(6):822-825.LIU Zhi-xin,YUE Jian-hua,LIU Yang-guang.Research on Application of Sector Scanning Technology on Advance Detecting[J].Journal of China University of Mining &Technology,2007,36(6):822-825.
[7]劉志新,于景邨,張華.小波變換在礦井瞬變電磁法中的應(yīng)用[J].煤田地質(zhì)與勘探,2007,35(4):70~72.LIU Zhi-xin,Yu Jing-cun,Zhang hua.The application of wavelet transfer technology in coal mine transient electromagnetic data[J].coal geology &exploration,2007,35(4):70~72.
[8]劉樹才,劉志新,姜志海.瞬變電磁法在煤礦采區(qū)水文勘探中的應(yīng)用[J].中國礦業(yè)大學(xué)學(xué)報,2005,34(4):414~417.LIU Shu-cai,LIU Zhi-xin,JIANG Zhi-hai.Application of TEM in hudrogeological prospecting of mining district .Journal of China University of Mining &Technology,2005,34(4):414-417
[9]劉志新,劉樹才,劉仰光.礦井富水體的瞬變電磁場物理模型實驗研究。巖石力學(xué)與工程學(xué)報[J].2009,28(2):259-266 LIU Zhi-xin,LIU Shu-cai,LIU Yang-guang.Research on transient electromagnetic field of mine water-bearing structure by physics model esperiment[J].chinese journal of rock mechanics and engineering,2009,28(2):259-266.
[10]于景邨,劉志新,劉樹才.深部采場突水構(gòu)造礦井瞬變電磁法探查理論及應(yīng)用[J].煤炭學(xué)報,2007,32(8)818-821.YU Jing-cun,LIU Zhi-xin,LIU Shu-cai.Theoretical analysis of mine transient electromagnetic method and its application in detecting water burst structures in deep coal stope.Journal of China Coal Society,2007,32(8)818-821
[11]YU Jing-cun,LIU Zhi-xin,TANG Jin-yun.Research on Full Space Transient Electromagnetism Technique for Detecting Aqueous Structures in Coal Mines[J].Journal of china university of mining and technology,2007,17(1):58-62
[12]于景邨,劉志新,岳建華.煤礦深部開采中的地球物理技術(shù)現(xiàn)狀及展望[J].地球物理學(xué)進展,2007,22(2)586~592 YU Jing-cun,LIU Zhi-xin,YUE Jian-hua.Present situation and prespect of geophysical technology in deep coal mining[J].Progress in Geophysics,2007,22(2)586~592
[13]Jiang Zhi-hai,Yue Jian-hua,Liu Shu-cai.Prediction technology of buried water-bearing structures in coal mines using t ransient electromagnetic method[J].Journal of china university of mining and technology,2007,vol.17(2):164-167