姚陳飛　王國(guó)富　李紅立

摘 要：直流電法超前探測(cè)技術(shù)是隧道掘進(jìn)工作面超前預(yù)報(bào)的重要手段之一，為研究異常體變化對(duì)直流電阻率法超前探測(cè)效果的影響，建立幾種含水低阻地質(zhì)異常體三維模型，采用三極法探討了在三維空間中不同距離、不同大小、不同方位、不同測(cè)線布置下的低阻地質(zhì)異常體直流電法超前探測(cè)響應(yīng)特征，以及地質(zhì)異常體與圍巖之間電阻率差異變化對(duì)于直流電阻率法超前探測(cè)的影響。研究結(jié)果表明：當(dāng)異常體體積越大，距離掘進(jìn)面越近，異常體與圍巖之間電阻率差異越大；異常體位于掘進(jìn)面正前方不發(fā)生偏移時(shí)，直流電法對(duì)異常體的響應(yīng)越明顯，三極法得到的探測(cè)效果越好；當(dāng)異常體位于不同方位，相應(yīng)的測(cè)線也位于該方位時(shí)，得到的探測(cè)效果最好，可以有效避免對(duì)隧道的影響。

關(guān)鍵詞：直流電法；超前探測(cè)；低阻地質(zhì)異常體；數(shù)值模擬

中圖分類號(hào)：U456.32；U452.11 DOI：10.16375/j.cnki.cn45-1395/t.2023.04.007

0 引言

電阻率法是根據(jù)巖石、礦石之間導(dǎo)電性的差別，研究巖石與礦石的電阻率變化，進(jìn)行地質(zhì)勘探的一種有效方法[1-2]。電阻率法對(duì)含水?dāng)鄬?、富水溶洞等低阻地質(zhì)體反應(yīng)敏感，分辨率高，在邊坡治理和超前探測(cè)方面具有獨(dú)特的優(yōu)勢(shì)[3-9]。

對(duì)于地質(zhì)異常體超前探測(cè)問(wèn)題，不少學(xué)者利用直流電阻率法對(duì)該問(wèn)題進(jìn)行了大量的研究。黃俊革等[10]通過(guò)構(gòu)建全空間板狀體模型，計(jì)算得出了點(diǎn)電源電位解析解，根據(jù)有限單元法得出坑道空腔對(duì)于極距較小的視電阻率幅值影響較大，而不會(huì)對(duì)迎頭前方地質(zhì)體的異常響應(yīng)造成較大的影響 。魯晶津等[11]應(yīng)用代數(shù)多重網(wǎng)格快速算法對(duì)二次場(chǎng)的有限差分問(wèn)題進(jìn)行求解得出巷道空腔對(duì)測(cè)量結(jié)果存在一定影響，采用比值曲線消除該影響，利用定點(diǎn)源對(duì)不同巷道面進(jìn)行測(cè)量的方式來(lái)識(shí)別旁側(cè)影響，達(dá)到準(zhǔn)確探測(cè)前方異常的目的。劉斌等[12]針對(duì)直流電阻率法超前探測(cè)研究中理論研究落后于工程實(shí)踐的現(xiàn)狀和探查效果較差的問(wèn)題，利用有限單元法進(jìn)行了正演，并提出了一套干擾識(shí)別與去除技術(shù)，形成了掌子面前方有用信息提取的有效途徑。馬炳鎮(zhèn)等[13]通過(guò)巷道影響因子研究了巷道對(duì)全空間穩(wěn)定電流場(chǎng)分布的影響，得出巷道影響與巷道幾何尺寸大小、供電電極布置及巷道圍巖導(dǎo)電性有關(guān)。王小龍等[14]利用COMSOL Multiphysics進(jìn)行正演模擬，根據(jù)巷道前方存在不同異常體時(shí)礦井超前探測(cè)視電阻率曲線的形態(tài)及極值點(diǎn)位置的不同，得出相對(duì)于板狀異常體，超前探測(cè)對(duì)球狀異常體探測(cè)效果更好。劉路[15]從數(shù)值模擬、物理模擬與工程實(shí)驗(yàn)等3個(gè)方面，對(duì)不同影響因素進(jìn)行了研究，得出離異常體較近測(cè)線的視電阻率曲線具有極值更明顯、幅度變化更大、尾支收斂速度更快的特點(diǎn)；通過(guò)多點(diǎn)源偏移處理，可有效排除旁側(cè)異常。張淼淼等 [16]進(jìn)行了不同空間方位異常體的探測(cè)，得到異常體存在傾斜角度非直立時(shí)，探測(cè)的結(jié)果與真實(shí)的位置之間存在偏差，導(dǎo)致異常的范圍有所擴(kuò)大，應(yīng)采取后期校正的方法來(lái)實(shí)現(xiàn)對(duì)異常體的準(zhǔn)確定位。

上述研究主要分析了巷道的影響以及針對(duì)消除巷道影響提出的對(duì)策，但對(duì)異常體的自身屬性未進(jìn)行有效分析，因此，有必要針對(duì)地質(zhì)異常體電阻率值的大小、地質(zhì)異常體方位、異常體的尺寸以及測(cè)線的布置等影響因素對(duì)異常體的響應(yīng)特征進(jìn)行研究，而且存在一定的可行性，可以有效地為直流電阻率法超前探測(cè)在實(shí)際生產(chǎn)中的應(yīng)用提供一定的參考。

1 直流電法的基本性質(zhì)

1.1 穩(wěn)定電流場(chǎng)基本定律

根據(jù)場(chǎng)論的知識(shí)可知，點(diǎn)電源產(chǎn)生的電場(chǎng)遵循歐姆定律可以看成是穩(wěn)定場(chǎng)，穩(wěn)定電流場(chǎng)與靜電場(chǎng)一樣只存在數(shù)值上的改變，都是勢(shì)場(chǎng)，[E]為電場(chǎng)強(qiáng)度，[σ]為電導(dǎo)率，[ρ]為電阻率，[j]為電流密度，[U]為電位，用公式表示為

為了說(shuō)明數(shù)值解的誤差能滿足精度要求，表1通過(guò)COMSOL Multiphysics軟件仿真的數(shù)值解與根據(jù)公式理論計(jì)算的解析解對(duì)比，電位的絕對(duì)誤差約為0.139 4 V，存在該誤差的原因與仿真軟件各個(gè)物理量耦合計(jì)算有關(guān)，在直流電法超前探測(cè)方法中觀測(cè)的是觀測(cè)點(diǎn)的電位差，因此，該誤差對(duì)數(shù)據(jù)處理的結(jié)果沒(méi)有影響，通過(guò)對(duì)數(shù)值解及解析解的分析，應(yīng)用COMSOL Multiphysics軟件對(duì)點(diǎn)源場(chǎng)模型的電位進(jìn)行模擬是可靠的。

3 直流電法數(shù)值模擬

由于含水溶洞具有低阻性質(zhì)，為了模擬這種地質(zhì)情況，用一個(gè)具有較低電阻率的塊狀體來(lái)代替含水溶洞進(jìn)行模擬，另外假設(shè)圍巖是由均勻介質(zhì)組成的地質(zhì)體，具有均勻電阻率，而且除了該異常體之外不含有其他地質(zhì)異常體情況，塊狀異常體與隧道軸線垂直且是直立的。

3.1 不同電阻率值異常體數(shù)值模擬

為了研究圍巖電阻率與地質(zhì)異常體電阻率差異對(duì)直流電阻率法探測(cè)視電阻率響應(yīng)的影響，含水構(gòu)造異常體電阻率分別取10、50、100 [Ω·m]，對(duì)其進(jìn)行數(shù)值模擬。

建立模型一圍巖大小為500 m×500 m×500 m，圍巖電阻率為500 [Ω·m]，截面積為20 m×20 m，厚度20 m，異常體位于掘進(jìn)面正前方10 m遠(yuǎn)處，與掘進(jìn)面中心在一條直線上，隧道空腔截面積為4 m×4 m，長(zhǎng)度100 m，電阻率為10 [Ω·m]，采用三極裝置，以掘進(jìn)面中心O點(diǎn)為原點(diǎn)，在地下A點(diǎn)（0，0，-2）供電，A點(diǎn)位于掘進(jìn)面與底板相交的中軸線上，測(cè)量點(diǎn)M點(diǎn)和N點(diǎn)也位于中軸線上，電流為1 A，AM極距6 m，MN極距為4 m，異常體不同電阻率模型示意圖如圖1所示。

圖2為不同電阻率異常體情況下的視電阻率曲線對(duì)比圖。從圖2視電阻率響應(yīng)曲線分析來(lái)看，圍巖電阻率與地質(zhì)異常體電阻率差異越大，曲線異常幅值越大，對(duì)應(yīng)的極值點(diǎn)越低，對(duì)異常體的響應(yīng)越明顯；隨著圍巖電阻率與地質(zhì)異常體電阻率之間的差異減小，曲線異常幅值減小，對(duì)應(yīng)的極值點(diǎn)上移，且尾支位置數(shù)值趨于圍巖電阻率值。說(shuō)明圍巖電阻率與地質(zhì)異常體電阻率差異越大，探測(cè)效果越好，當(dāng)異常體電阻率為10 [Ω·m]時(shí)，探測(cè)效果最佳。

3.2 異常體距掘進(jìn)面不同距離數(shù)值模擬

為了研究異常體距掘進(jìn)面不同距離對(duì)超前探測(cè)視電阻率響應(yīng)的影響，當(dāng)異常體電阻率為10 [Ω·m]，大小為20 m×20 m×20 m時(shí)，改變異常體正對(duì)掘進(jìn)面距離為10、30、50、70 m位置，對(duì)其進(jìn)行數(shù)值模擬，其他物性條件與模型一相同，異常體距掘進(jìn)面不同距離模型示意圖如圖3所示。

圖4為異常體距掘進(jìn)面不同距離的視電阻率曲線對(duì)比圖。由圖4視電阻率響應(yīng)曲線可知，地質(zhì)異常體距點(diǎn)電源越近，極值點(diǎn)越低，異常幅值越大。當(dāng)?shù)刭|(zhì)異常體在距點(diǎn)電源50 m和70 m處，異常幅值變化一致，極值點(diǎn)重合，說(shuō)明大于50 m之后，對(duì)于異常體的探測(cè)距離遠(yuǎn)近已經(jīng)分辨不出，通過(guò)比較極值點(diǎn)的位置以及異常幅值的大小可以得出，當(dāng)異常體電阻率為10 [Ω·m]時(shí)，在塊狀異常體位于掘進(jìn)面正前方10 m處，探測(cè)效果最好。

3.3 不同異常體大小數(shù)值模擬

為了研究異常體大小對(duì)超前探測(cè)視電阻率響應(yīng)曲線的影響，當(dāng)異常體電阻率為10 [Ω·m]時(shí)，在塊狀異常體位于掘進(jìn)面正前方10 m處，建立不同尺寸大小的異常體模型進(jìn)行對(duì)比分析，改變正方體邊長(zhǎng)a 分別為2、5 、10 、20 m，進(jìn)行數(shù)值模擬，其他物性條件與模型一相同。掘進(jìn)面前方不同尺寸異常體示意圖如圖5所示。

圖6為采用三極裝置點(diǎn)電源正前方不同尺寸大小的異常體正演結(jié)果（視電阻率）曲線對(duì)比圖。從圖6的5條曲線可以看出，4條曲線都出現(xiàn)異常幅值不同程度的變化，除了前方無(wú)異常體的情況外，其他4條曲線都呈現(xiàn)先下降后上升的趨勢(shì)，隨著異常體的體積逐漸減小，極值點(diǎn)不斷上移，異常幅值逐漸減小。根據(jù)異常幅值和極值點(diǎn)的變化，尺寸為2 m×2 m×2 m、5 m×5 m×5 m、10 m×10 m×10 m的異常體最大幅值分別只有25 、28 和40 [Ω·m]，尺寸為20 m×20 m×20 m的異常體幅值最大達(dá)到64 [Ω·m]，是其他尺寸大小異常體的2倍左右，由此可以得出，當(dāng)異常體電阻率為10 [Ω·m]時(shí)，在塊狀異常體位于掘進(jìn)面正前方10 m處，直流電法對(duì)于20 m×20 m×20 m尺寸的異常體探測(cè)效果明顯好于其他尺寸的異常體探測(cè)效果。

3.4 不同偏移方位數(shù)值模擬

為了研究掘進(jìn)面前方異常體不同方位對(duì)超前探測(cè)視電阻率響應(yīng)的影響，當(dāng)異常體電阻率為10 [Ω·m]，大小為20 m×20 m×20 m時(shí)，在塊狀異常體位于掘進(jìn)面正前方10 m處，將塊狀異常體向上偏移15 m、不偏移以及向下偏移15 m進(jìn)行對(duì)比分析，其他物性條件與模型一相同，掘進(jìn)面前方不同方位異常體模型示意圖如圖7所示。

圖8為異常體分別發(fā)生向上、不偏移、向下偏移15 m視電阻率曲線對(duì)比圖。從圖8可以看出，當(dāng)塊狀異常體未發(fā)生偏移時(shí)，視電阻率曲線極值點(diǎn)最小，異常幅值最大，隨著異常體發(fā)生偏移，視電阻率曲線極值點(diǎn)上移，異常幅值逐漸減小，根據(jù)正演結(jié)果可以得出，當(dāng)異常體大小為20 m×20 m×20 m、電阻率為10 [Ω·m]時(shí)，在塊狀異常體位于掘進(jìn)面正前方10 m，不發(fā)生上下偏移時(shí)，直流電法對(duì)于異常體的識(shí)別效果最好。

3.5 不同測(cè)線布置數(shù)值模擬

為了研究測(cè)線布置位置對(duì)超前探測(cè)視電阻率響應(yīng)的影響，假設(shè)平行隧道空腔左幫旁側(cè)15 m處，有一個(gè)邊長(zhǎng)為10 m的正方形異常體，位于掘進(jìn)面后方20 m處，且掘進(jìn)面前方無(wú)異常體，測(cè)線1布置于隧道底板中軸線，測(cè)線2布置于隧道底板與左幫相交位置，測(cè)線3布置于左幫中軸線位置，其他物性條件與模型一相同，不同測(cè)線布置模型示意圖如圖9所示。

圖10為不同測(cè)線布置視電阻率曲線圖。從圖10可以看出，測(cè)線3異常幅值大于測(cè)線2異常幅值，而測(cè)線2異常幅值大于測(cè)線1異常幅值，測(cè)線3的極值點(diǎn)相較于測(cè)線1和測(cè)線2的極值點(diǎn)更低，說(shuō)明測(cè)線3的布置對(duì)于異常體探測(cè)效果更好。從上面的測(cè)線布置結(jié)果可以得出，為了避免隧道空腔對(duì)探測(cè)結(jié)果產(chǎn)生較大的干擾，針對(duì)不同方位的異常體，需要在該方位布置測(cè)線進(jìn)行觀測(cè)，可以得出有效的異常響應(yīng)結(jié)果。

4 結(jié)論

通過(guò)數(shù)值模擬，分析研究了異常體電阻率大小、距離遠(yuǎn)近、方位偏移、體積大小以及測(cè)線布置等影響因素對(duì)直流電法超前探測(cè)效果的影響，得出以下結(jié)論：

1）采用COMSOL Multiphysics軟件進(jìn)行三維數(shù)值模擬，對(duì)比數(shù)值解與解析解之間的誤差，直流電阻率法數(shù)值模擬精度較高，能夠滿足精度要求。

2）從不同電阻率地質(zhì)異常體的視電阻率曲線來(lái)看，當(dāng)?shù)刭|(zhì)異常體與圍巖之間的電阻率差異較大時(shí)，異常幅值大，極小值點(diǎn)小，探測(cè)效果比較好。

3）通過(guò)探測(cè)距點(diǎn)電源不同距離的地質(zhì)異常體發(fā)現(xiàn)，距離點(diǎn)電源50 m和70 m所得到的視電阻率曲線幾乎重合，無(wú)法分辨出它們的差異，由此可以得出三極裝置探測(cè)距離大于50 m，分辨率大大降低。

4）對(duì)比不同尺寸大小的地質(zhì)異常體，當(dāng)異常體的邊長(zhǎng)為20 m時(shí)，所得到的異常幅值差異最大，探測(cè)效果明顯更好。

5）對(duì)于正前方地質(zhì)異常體是否發(fā)生偏移的情況，當(dāng)異常體未發(fā)生偏移時(shí)，異常幅值最大，極小值點(diǎn)最小，探測(cè)效果最佳。

6）對(duì)于異常體位于不同方位時(shí)，測(cè)線需布置于該方位，以減少隧道空腔對(duì)于異常體探測(cè)效果的影響，以得到更好的探測(cè)效果。

參考文獻(xiàn)

[1] 岳建華，劉樹才. 礦井直流電法勘探[M].徐州：中國(guó)礦業(yè)大學(xué)出版社，1999.

[2] 李金銘. 地電場(chǎng)與電法勘探[M].北京：地質(zhì)出版社，2005.

[3] 張平松，劉盛東，曹煜. 坑道掘進(jìn)立體電法超前預(yù)報(bào)技術(shù)研究[J].中國(guó)煤炭地質(zhì)，2009，21（2）：50-53.

[4] 石學(xué)鋒. 礦井直流電法超前探測(cè)影響因素?cái)?shù)值模擬[J].煤炭技術(shù)，2016，35（11）：122-124.

[5] 占文鋒，武玉梁，李文. 礦井直流電法全空間電場(chǎng)分布數(shù)值模擬及影響因素[J].煤田地質(zhì)與勘探，2018，46（1）：139-147.

[6] 周圓兀，韓曉偉，陳鮑. 基于不平衡推力法的邊坡穩(wěn)定性系數(shù)優(yōu)化計(jì)算方法[J].廣西科技大學(xué)學(xué)報(bào)，2016，27（1）：71-74.

[7] 趙冠宇，王敬，翟培合，等. 基于COMSOL Multiphysics數(shù)值模擬的礦井直流電法超前探測(cè)研究[J].煤炭技術(shù)，2018，37（8）：177-179.

[8] 胡雄武，張平松. 坑道隱伏陷落柱直流電阻率法超前探測(cè)分析[J].地球物理學(xué)進(jìn)展，2019，34（3）：1176-1183.

[9] 周圓兀，陳鮑，ANNIL M. 基于地震作用下傳遞系數(shù)法的邊坡穩(wěn)定性分析[J].廣西科技大學(xué)學(xué)報(bào)，2017，28（1）：25-28.

[10] 黃俊革，王家林，阮百堯. 坑道直流電阻率法超前探測(cè)研究[J].地球物理學(xué)報(bào)，2006，49（5）：1529-1538.

[11] 魯晶津，吳小平. 巷道直流電阻率法超前探測(cè)三維數(shù)值模擬[J].煤田地質(zhì)與勘探，2013，41（6）：83-86.

[12] 劉斌，李術(shù)才，李樹忱，等. 隧道含水構(gòu)造電阻率法超前探測(cè)正演模擬與應(yīng)用[J].吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2012，42（1）：246-253.

[13] 馬炳鎮(zhèn)，李貅. 礦井直流電法超前探中巷道影響的數(shù)值模擬[J].煤田地質(zhì)與勘探，2013，41（1）：78-81.

[14] 王小龍，馮宏，李萍，等. 礦井電阻率超前探測(cè)正演模擬[J].煤炭科學(xué)技術(shù)，2011，39（11）：112-117.

[15] 劉路. 礦井直流電阻率法三維超前探測(cè)技術(shù)研究[D].徐州：中國(guó)礦業(yè)大學(xué)，2014.

[16] 張淼淼，石顯新. 不同空間方位異常體的直流電阻率法分析[J].科學(xué)技術(shù)與工程，2022，22（2）：481-487.

Advanced detection analysis of water-rich structural anomalies

by DC resistivity method

YAO Chenfeia， WANG Guofub， LI Hongli*b

（a. School of Automation， Guangxi University of Science and Technology， Liuzhou 545616， China; b. School of Electronic Engineering， Guangxi University of Science and Technology， Liuzhou 545006， China）

Abstract： Direct current （DC） advanced detection technology is one of the important means of advanced prediction of tunnel facing. In order to study the influence of anomalous changes on the advanced detection effect of DC resistivity method， several three-dimensional models of water-bearing low-resistance geological anomalies were established， and the tripolar method was used to explore the response characteristics of the advanced detection of low-resistance geological anomalies under different distances， sizes， orientations and line arrangements in three-dimensional space， and the influence of the change of resistivity difference between geological anomalies and surrounding rock on the advanced detection of DC resistivity method. The results show that the larger the volume of the anomaly and the closer to the boring surface， the greater the resistivity difference between the anomaly and the surrounding rock. When the anomaly is located directly in front of the boring surface and does not shift， the more obvious the response of the direct current method to the anomaly， the better the detection effect obtained by the tripolar method. When the anomalies are located in different orientations and the corresponding measurement lines are also located in this orientation， the detection effect is the best， which can effectively avoid the impact on the tunnel.

Key words： direct current method; advanced detection; low-resistance geological anomalies; numerical simulation

（責(zé)任編輯：羅小芬）

收稿日期：2022-12-07

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目（61761009）資助

第一作者：姚陳飛，在讀碩士研究生

*通信作者：李紅立，博士，助理工程師，研究方向：綜合地球物理的教學(xué)及研究，E-mail：22723071@qq.com