李帝銓,肖教育,張繼峰,胡艷芳,劉最亮,張新
(1.中南大學(xué) 有色金屬成礦預(yù)測(cè)與地質(zhì)環(huán)境監(jiān)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室,湖南 長(zhǎng)沙 410083; 2.有色資源與地質(zhì)災(zāi)害探測(cè)湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 長(zhǎng)沙 410083; 3.中南大學(xué) 地球科學(xué)與信息物理學(xué)院,湖南 長(zhǎng)沙 410083; 4.長(zhǎng)安大學(xué) 地球科學(xué)與資源學(xué)院,陜西 西安 710000; 5.陽(yáng)泉煤業(yè)(集團(tuán))有限責(zé)任公司 地質(zhì)測(cè)量部,山西 陽(yáng)泉 045000)
煤礦開采時(shí),煤礦隱蔽致災(zāi)地質(zhì)因素是影響煤炭安全開采的首要問題,這些地質(zhì)因素包括地質(zhì)構(gòu)造、礦井水、煤層瓦斯、頂?shù)装鍡l件等,煤層中富水區(qū)是煤礦開采中重大隱蔽地質(zhì)隱患之一。地球物理探測(cè)是查明這些地質(zhì)因素的重要手段,常用的電磁勘探方法有可控源音頻大地電磁法(CSAMT)、瞬變電磁法、巷道無(wú)線電波透視等。廣域電磁法是相對(duì)于傳統(tǒng)的CSAMT法和MELOS方法提出來(lái)的,使用人工場(chǎng)源發(fā)射電流,采用單分量或多分量測(cè)量和適合于全域的公式來(lái)計(jì)算廣域視電阻率,供電頻率的帶寬選擇較大,適合大深度勘探[1-4],在油氣勘探、固體礦產(chǎn)勘探、地?zé)峥碧絒5-10]等領(lǐng)域應(yīng)用較為廣泛。廣域電磁法有可能在煤炭水害富水區(qū)探測(cè)領(lǐng)域做出貢獻(xiàn),但是之前還沒有具體實(shí)施應(yīng)用。為了比較廣域電磁法和可控源音頻大地電磁法的探測(cè)效果,證實(shí)廣域電磁法在煤田中探查水的可行性,在山西沁水煤田新元煤礦進(jìn)行了對(duì)比實(shí)驗(yàn)。這是廣域電磁法在煤礦系統(tǒng)富水探測(cè)中的首次應(yīng)用。
本次廣域電磁法采取了多維多分量的方式,即采用旁側(cè)裝置和軸向裝置分別測(cè)量電場(chǎng)的Ex、Ey分量,從而獲得同一測(cè)點(diǎn)4個(gè)不同的電場(chǎng)分量,每個(gè)電場(chǎng)分量均求出相應(yīng)的廣域視電阻率。理論分析與試驗(yàn)數(shù)據(jù)表明,在野外采用旁側(cè)裝置測(cè)量電場(chǎng)Ex分量,測(cè)線敷設(shè)和收攏方便,移動(dòng)靈活,信號(hào)強(qiáng)度大,更適用于實(shí)際勘探[11]。本文主要對(duì)廣域電磁法的旁側(cè)裝置Ex、軸向裝置Ex的頻率視電阻率擬斷面圖和CSAMT的頻率視電阻率擬斷面進(jìn)行對(duì)比;同時(shí)對(duì)信號(hào)強(qiáng)度較大、分辨效果好的旁側(cè)Ex裝置的反演結(jié)果與CSAMT的反演結(jié)果進(jìn)行對(duì)比,從而說明廣域電磁法在煤炭領(lǐng)域富水區(qū)探測(cè)的可行性。
廣域電磁法是一種采用2n序列偽隨機(jī)信號(hào)作為激勵(lì)信號(hào)的人工源頻率域電磁法,該方法采用水平電流源或者垂直磁偶源作為場(chǎng)源,通過測(cè)量電磁場(chǎng)的某個(gè)分量以獲得廣域視電阻率。
根據(jù)均勻大地表面水平電流源產(chǎn)生的電場(chǎng)Ex的精確表達(dá)式如下:
(1)
由式(1)得到廣域視電阻率表達(dá)式:
(2)
式中:I為發(fā)送電流;dL為場(chǎng)源尺寸;r為收發(fā)距;φ為方位角;k為波數(shù);σ為介質(zhì)電導(dǎo)率;ΔVMN為觀測(cè)電位差,MN為電極距;KE-Ex是一個(gè)只與觀測(cè)裝置的幾何尺寸有關(guān)的系數(shù);FE-Ex(ikr)是由發(fā)送電流頻率ω、地下電阻率ρ以及收發(fā)距r構(gòu)成的復(fù)函數(shù)。采用式(2)定義的廣域視電阻率可不受觀察點(diǎn)到場(chǎng)源距離的限制,因此,利用計(jì)算機(jī)迭代方法求解的視電阻率最佳值適用于廣大區(qū)域[10,12]。
可控源音頻大地電磁法是利用人工源激發(fā)地下巖石產(chǎn)生的電導(dǎo)率差異觀測(cè)一次場(chǎng)電位和磁場(chǎng)強(qiáng)度變化的一種電磁勘探方法,其頻率范圍一般為0.25~8 192 Hz。主要針對(duì)大地電磁測(cè)深法場(chǎng)源的隨機(jī)性和信號(hào)微弱,以致觀測(cè)十分困難的狀況提出的一種改進(jìn)方案[12-13]。
CSAMT法采用水平電偶極子或垂直磁偶極子作為場(chǎng)源,在遠(yuǎn)區(qū)測(cè)量一對(duì)相互正交的電場(chǎng)和磁場(chǎng)分量,沿用MT法采用波阻抗提取視電阻率的公式:
(3)
式中:ρa(bǔ)為視電阻率,ω為電磁波的頻率,μ為介質(zhì)的磁導(dǎo)率。
CSAMT和WFEM在野外施工場(chǎng)源布置和測(cè)線布置方面有所差異,相同的目標(biāo)體探測(cè)深度下,WFEM所需的收發(fā)距較小。CSAMT測(cè)量相互正交的電場(chǎng)和磁場(chǎng),WFEM只需測(cè)量單一的電場(chǎng)或磁場(chǎng)即可進(jìn)行廣域視電阻率的換算以及反演。CSAMT和WFEM在數(shù)據(jù)處理方面由于各方法視電阻率的定義不同,處理方法也有所差異。
新元井田31004工作面屬華北石炭—二疊系煤田,井田基本形態(tài)為一單斜,井田位于沁水煤田的西北部,陽(yáng)煤集團(tuán)壽陽(yáng)區(qū)中南部。新元井田范圍內(nèi)自上而下的含水層組主要包括:第四系砂礫石層孔隙含水層組,二疊系石盒子組、石千峰組、三疊系劉家溝組砂巖裂隙含水層組,二疊系下統(tǒng)山西組砂巖裂隙含水層組,石炭系上統(tǒng)太原組石灰?guī)r溶及砂巖裂隙含水層和中奧陶統(tǒng)峰峰組與上馬家溝組石灰?guī)r巖溶水含水層組(圖1)。
研究區(qū)主要地層從上到下依次有第四系(Q)、二疊系上石盒子組(P2s)、二疊系下石盒子組(P1x)、二疊系山西組(P1s)、石炭系太原組(C2t)、石炭系本溪組(C2b)、奧陶系峰峰組(O2f)、奧陶系上馬家溝組(O2s)、奧陶系下馬家溝組(O2x)。根據(jù)已有的測(cè)井電阻率信息統(tǒng)計(jì)分析,測(cè)區(qū)地層電阻率總體呈現(xiàn)高阻—低阻—高阻—低阻—高阻的趨勢(shì)(表1)。研究區(qū)主要可采煤層為3號(hào)煤、9號(hào)煤、15號(hào)煤。富水區(qū)的目的層位于山西組和太原組中,均表現(xiàn)為低阻特征。
表1 地層及圍巖物性參數(shù)統(tǒng)計(jì)
此次在新元31004工作面所屬地面共布置2條廣域電磁法測(cè)線——280線、285線,其中280線為廣域電磁法和可控源音頻大地電磁法在相同地質(zhì)條件下的共同測(cè)線。280線總長(zhǎng)3 km,點(diǎn)距20 m,共150個(gè)測(cè)點(diǎn)(圖2)。
圖2 測(cè)線布置示意Fig.2 Survey line layout
廣域電磁法野外采集使用中南大學(xué)自主研發(fā)的廣域電磁發(fā)射機(jī)和廣域電磁接收儀,采用旁側(cè)裝置(E-Ex、E-Ey)和軸向裝置(E-Ex、E-Ey),發(fā)射電流為120 A、發(fā)射頻率為0.75~8 192 Hz,共28個(gè)頻點(diǎn)。此次采用的旁側(cè)裝置和軸向裝置(圖3)與傳統(tǒng)的廣域電磁法區(qū)別在于利用2個(gè)相互正交的源進(jìn)行供電,例如:當(dāng)采用軸向裝置進(jìn)行觀測(cè)時(shí),分別沿著平行于源的方向采集Ex分量,沿著垂直于源的方向采集Ey分量;旁側(cè)裝置亦是如此。CSAMT野外工作采用美國(guó)ZONG公司研發(fā)的GDP32- II電法工作站,信號(hào)發(fā)射和接收采用與之配套的裝置,發(fā)射電流20 A,測(cè)量頻率0.125~8 192 Hz,共30個(gè)頻點(diǎn)。廣域電磁法和可控源音頻大地電磁法收發(fā)距同為8.99 km,場(chǎng)源位置相同,供電電極距1.0 km。
圖3 廣域電磁法觀測(cè)裝置示意Fig.3 Schematic diagram of WFEM's observation device
采用中南大學(xué)自主設(shè)計(jì)的廣域電磁法數(shù)據(jù)處理和反演一體化解釋平臺(tái),經(jīng)過對(duì)數(shù)據(jù)進(jìn)行預(yù)處理、“飛點(diǎn)”剔除、靜態(tài)校正、數(shù)據(jù)濾波、定性分析、參數(shù)點(diǎn)分析以及曲線類型分析,再建立模型作數(shù)據(jù)反演,最后依據(jù)其他物探成果進(jìn)行綜合推斷解釋??煽卦匆纛l大地電磁法采用Mtsoft2D軟件進(jìn)行處理,包括近場(chǎng)校正、去除“飛點(diǎn)”、靜態(tài)校正以及反演電阻率斷面成圖,最后進(jìn)行結(jié)果解釋。
依據(jù)勘探區(qū)地球物理特征,將上石盒子組(P2s)電性特征為高阻,下石盒子組電性特征為低阻(P1x),以及山西組(P1s)、太原組(C2t)整體電性特征呈現(xiàn)高阻,本溪組(C2b)的電性特征為低阻,作為反演電阻率斷面圖的地質(zhì)解釋原則。
圖4、圖5分別為新元礦31004工作面280線廣域電磁法和可控源音頻大地電磁法所測(cè)得山頂處4個(gè)號(hào)位的視電阻率曲線??梢钥闯觯篊SAMT法卡尼亞視電阻率曲線由于受到場(chǎng)源效應(yīng)的影響在低頻段進(jìn)入過渡區(qū)和近區(qū),導(dǎo)致CSAMT法卡尼亞視電阻率曲線在低頻段出現(xiàn)明顯的上翹趨勢(shì),在10~60 Hz處會(huì)出現(xiàn)一個(gè)電阻率極小值,此時(shí)處于過渡區(qū),隨后CSAMT法卡尼亞視電阻率曲線由極小值呈45°上升的趨勢(shì)。CSAMT法卡尼亞視電阻率在235號(hào)點(diǎn)和237號(hào)點(diǎn)高頻處出現(xiàn)“飛點(diǎn)” 而廣域電磁法整體的電阻率曲線較為平滑,整體呈現(xiàn)高阻—低阻—高阻—低阻—高阻的趨勢(shì),這與實(shí)際的地質(zhì)資料較為吻合。
圖4 235、237號(hào)點(diǎn)CSAMT和WFEM視電阻曲線對(duì)比Fig.4 Comparison of apparent resistance curve between CSAMT and WFEM at point 235 and 237
圖5 381、383號(hào)點(diǎn)CSAMT和WFEM視電阻率曲線對(duì)比Fig.5 Comparison of apparent resistance curve between CSAMT and WFEM at point 381 and 383
圖6為新元礦31004工作面280線廣域電磁法和CSAMT的視電阻率擬斷面對(duì)比??梢钥闯鰪V域電磁法的旁側(cè)裝置Ex(圖6a)和軸向裝置Ex(圖6b)的原始數(shù)據(jù)頻率視電阻率擬斷面圖在低頻段依然有響應(yīng),說明受到場(chǎng)源效應(yīng)的影響很小,不受過渡區(qū)和近區(qū)的影響;而CSAMT法卡尼亞視電阻率曲線(圖6c)由于受到場(chǎng)源效應(yīng)的影響,在低頻段進(jìn)入過渡區(qū)和近區(qū)導(dǎo)致CSAMT法卡尼亞視電阻率曲線在低頻段出現(xiàn)明顯的上翹趨勢(shì),所以導(dǎo)致在低頻段無(wú)法對(duì)地層進(jìn)行劃分。同時(shí),可以認(rèn)為廣域電磁法旁側(cè)裝置Ex和軸向裝置Ex對(duì)富水體分辨范圍更大;綜合比較圖6a和圖6b,可以認(rèn)為旁側(cè)裝置Ex對(duì)富水體的分層能力優(yōu)于軸向裝置Ex。
圖6 280線WFEM和CSAMT視電阻率擬斷面對(duì)比Fig.6 Comparison of apparent resistivity of original data of WFEM and CSAMT
最新淋水點(diǎn)靠近輔助進(jìn)風(fēng)巷,285線與之較為接近。最新的31004工作面淋水點(diǎn)最大出水量達(dá)到25 m3/h,約靠近輔助進(jìn)風(fēng)巷60 m,與285測(cè)線距離為1 620 m。在285線電阻率反演斷面(圖7)中進(jìn)行了標(biāo)注,這與實(shí)際的淋水點(diǎn)較為吻合,證明了電阻率反演斷面的準(zhǔn)確性;淋水點(diǎn)處電阻率約為35 Ω·m,這與我們統(tǒng)計(jì)的圍巖電阻率和出水點(diǎn)電阻率(ρs≤35 Ω·m)也是吻合的,推斷水來(lái)自3號(hào)煤層上部K8含水層。
圖7 285線反演結(jié)果與最大淋水點(diǎn)驗(yàn)證Fig.7 285 line inversion results and water outlet verification
在對(duì)CSAMT數(shù)據(jù)反演之前需要對(duì)原始數(shù)據(jù)進(jìn)行整理,剔除個(gè)別畸變點(diǎn),通過可視化MT-Pioneer5.1數(shù)據(jù)處理軟件進(jìn)行下一步操作,不進(jìn)行人為數(shù)據(jù)的平滑,對(duì)因干擾產(chǎn)生的“飛點(diǎn)”直接剔除以提高反演的準(zhǔn)確性。選擇二維非線性共軛梯度法進(jìn)行反演,模型光滑度設(shè)置為80,約束視電阻率的誤差為9%。
廣域電磁法資料處理解釋以巖石物性為基礎(chǔ),通過對(duì)原始數(shù)據(jù)飛點(diǎn)剔除、靜態(tài)校正等預(yù)處理并結(jié)合原始數(shù)據(jù)“頻率—視電阻率”曲線完成定性分析,然后建立合理的地質(zhì)模型,進(jìn)行一維連續(xù)介質(zhì)反演,并在一維反演的基礎(chǔ)上進(jìn)行二維反演,結(jié)合電阻率測(cè)井資料解釋,提高廣域電磁法解釋成果的可靠性和解釋精度。利用一維反演的結(jié)果可進(jìn)行斷層劃分,利用二維連續(xù)介質(zhì)反演結(jié)果可進(jìn)行地層劃分[10],最后綜合一維連續(xù)介質(zhì)和二維連續(xù)介質(zhì)反演結(jié)果,對(duì)富水區(qū)進(jìn)行精確的識(shí)別。
圖8為WFEM和CSAMT的反演結(jié)果對(duì)比。通過對(duì)圍巖電阻率和巷道出水點(diǎn)的電阻率分析,可將ρ≤35 Ω·m劃為低阻異常區(qū)。從反演電阻率斷面看,K10含水層在WFEM反演斷面圖中距離0~500 m、1 300~3 200 m富水性較強(qiáng),且該富水區(qū)連續(xù)性較強(qiáng)。從WFEM反演電阻率斷面看,K8含水層在WFEM反演斷面圖中距離0~1 200 m、1 300~1 700m、2 000~2 100 m處富水性較強(qiáng)。
圖8 280線WFEM(a)和CSAMT(b)的反演結(jié)果對(duì)比Fig.8 Comparison of inversion results of WFEM (a) and CSAMT (b) of line 280
通過對(duì)WFEM的反演斷面圖和CSAMT的反演斷面圖對(duì)比分析可看出:WFEM反演結(jié)果呈現(xiàn)低阻—高阻—低阻—高阻—低阻—高阻的關(guān)系,對(duì)第四系和上石盒子組、下石盒子組、山西組、太原組—奧陶系有很好的分辨能力,而CSAMT的反演結(jié)果呈現(xiàn)低阻—高阻—低阻—高阻的關(guān)系;WFEM反演結(jié)果對(duì)K10含水層和K8含水層的分辨較為清晰,而CSAMT的反演結(jié)果對(duì)K10含水層和K8含水層的分辨較為模糊,分辨不清。
1)廣域電磁法是基于全區(qū)視電阻率來(lái)定義的,對(duì)視電阻率的求解公式?jīng)]有近似取舍,因此測(cè)量范圍大、測(cè)量精度高??煽卦匆纛l大地電磁法受到場(chǎng)源效應(yīng)和近區(qū)、過渡區(qū)的影響導(dǎo)致卡尼亞視電阻率曲線在低頻段產(chǎn)生45°上揚(yáng)的趨勢(shì)。
2)廣域電磁法采用2n偽隨機(jī)信號(hào)多頻波發(fā)射,相對(duì)于可控源音頻大地電磁法的單頻波發(fā)射,其深部的分辨率較高。同時(shí)只需觀測(cè)單一電場(chǎng)Ex分量即可求得廣域視電阻率。
3)廣域電磁法首次在煤礦中進(jìn)行富水區(qū)探查應(yīng)用,經(jīng)過此次對(duì)比可發(fā)現(xiàn)廣域電磁法在煤礦中對(duì)富水體的探查高效、準(zhǔn)確度高。也證實(shí)了廣域電磁法在煤炭領(lǐng)域富水區(qū)探測(cè)具有推廣應(yīng)用的巨大潛力,可為煤炭領(lǐng)域富水區(qū)探測(cè)提供新技術(shù)。