邱衛(wèi)忠，張 輝，唐新功，嚴(yán)良俊，張雷潔

(1.山西省煤田地質(zhì)115勘察院，山西 大同 037039;2.長江大學(xué) 地球物理與石油資源學(xué)院，湖北 武漢 430100)

?

煤田積水采空區(qū)井地電位法三維有限差分正演模擬

邱衛(wèi)忠1，張 輝2，唐新功2，嚴(yán)良俊2，張雷潔1

(1.山西省煤田地質(zhì)115勘察院，山西 大同 037039;2.長江大學(xué) 地球物理與石油資源學(xué)院，湖北 武漢 430100)

煤田積水采空區(qū)是威脅礦井安全的主要因素之一。針對井地電位法特定的邊值問題，推導(dǎo)出了相應(yīng)的有限差分格式，結(jié)合煤田積水采空區(qū)的電性特征及分布特點，設(shè)計了多個地電模型，分別利用有限差分法對其進(jìn)行了正演模擬。正演結(jié)果表明，井地電位法對單個積水采空區(qū)及橫向上分布的多個積水采空區(qū)識別能力較好，而對垂向分布的多個積水采空區(qū)識別能力較弱。研究結(jié)果對于煤田積水采空區(qū)的探測研究具有一定的參考價值。

煤田積水采空區(qū)；井地電位法；有限差分法；三維正演

1 引 言

礦井的水害事故是煤田生產(chǎn)中重要的災(zāi)害之一，是僅次于瓦斯事故的第二大煤礦事故，其中采空區(qū)積水在水害事故占有相當(dāng)大的比重[1]。同時，采空區(qū)積水具有突發(fā)性強(qiáng)，突水量大，破壞性強(qiáng)等特點，是威脅礦井安全的主要因素之一[2]。因此，做好煤田積水采空區(qū)的探測工作，對礦井水害事故的防治具有積極的意義。

作為無損探測的電磁探測方法，由于其對含水體敏感、分辨能力較高，因此在煤田水文地質(zhì)勘探中獲得了廣泛的應(yīng)用。但在地下結(jié)構(gòu)復(fù)雜，目的層埋深又較深的情況下，地面物探方法會由于探測深度不夠、信號不強(qiáng)、干擾太大等原因，無法采集到高質(zhì)量的物探數(shù)據(jù)，給后期的資料處理與解釋造成很大的困難，不能得出精確的結(jié)果。

為克服常規(guī)地面物探方法的缺陷，地—井物探方法開始逐漸興起，井地電位法便是其中的代表之一。它是近幾年發(fā)展起來的一種電磁勘探新方法和新技術(shù)。井地電位法是通過套管井向地下供入大功率電流，在地表觀測由于地下介質(zhì)的電性變化引起的電位異常，據(jù)此來分析地下導(dǎo)電體的電性分布情況。相比于傳統(tǒng)的電阻率法，井地電位法具有勘探深度大、探測范圍廣、分辨率高等特點[3]。

從20世紀(jì)70年代開始，一直到本世紀(jì)初，井地電位法在圈定煤層、油氣邊界以及確定深部礦藏等領(lǐng)域不斷取得成果[4]。這一階段對井地電位法的研究，還集中在定性解釋和物理模擬[11，12]上，數(shù)值模擬研究較少。此后，徐凱軍等[13]、劉地淵等[14]、蘇朱劉等[15]、王志剛等[16]、屈有恒等[17]、柯敢攀等[18]先后分別使用多種數(shù)值方法，對油田井地電位法的三維正反演問題展開了研究。近年來，戴前偉等[19]、李蕓蕓[3]還開展了復(fù)雜地形條件下的井地電位法三維數(shù)值模擬研究，取得了一定的效果。

隨著井地電位法在煤田積水采空區(qū)探測中的進(jìn)一步應(yīng)用，結(jié)合煤田積水采空區(qū)的特點，繼續(xù)深入開展煤田井地電位法的三維數(shù)值模擬是非常必要的。有限差分法作為應(yīng)用最早、最經(jīng)典的數(shù)值模擬方法，依靠其相對簡單的理論基礎(chǔ)及簡便的計算過程，在地球物理工作中應(yīng)用非常廣泛[20-22]。本文利用有限差分算法，結(jié)合井地電位方法在煤田積水采空區(qū)中的應(yīng)用，進(jìn)行三維正演模擬研究，探討地井電位法對不同位置積水采空區(qū)的探測效果。

2 井地電位法的邊值問題

井地電位法的工作原理如圖1所示，在地表向被測井A施加一個大功率電流，同時在無窮遠(yuǎn)處選擇另一口井B作回流電極。兩個供電電極的金屬套管，全部用水泥封固。由于水泥膠結(jié)致密，因此電流主要通過地表及地下目的層流向地層，造成地層中的電流密度分布不均勻，但是電流幾乎垂直于套管壁流出。同時，由于套管直徑與套管長度相比可以忽略，因此在井地電位法的正演模擬、反演成像以及后期的數(shù)據(jù)處理過程中，一般都將套管當(dāng)作線電流源來處理。

圖1 井地電位法工作原理示意圖Fig.1 Schematic map of borehole-surface potential method

在非均勻各向同性的三維介質(zhì)中，點電流源產(chǎn)生的電位U可以用如下的邊值問題來描述[3]：

(1-1)

(1-2)

(1-3)

式(1)中，I為供電電流大小，σ為電阻率ρ的倒數(shù)，δ(A)為δ函數(shù)，Γs和Γ∞分別是研究區(qū)域Ω的地表邊界和地下邊界，n為地下邊界的外法線方向，r為點電流源A到邊界點的距離，cos(r,n)為n與r之間夾角的余弦值。

線電流源的本質(zhì)，其實可以看成是無數(shù)個點電流源的疊加。而根據(jù)疊加原理，幾個點電源在空間中某處電位的大小，等于各個點電源單獨存在時在該點產(chǎn)生的電位大小的代數(shù)和。因此，在井地電位法對線電流源的正演模擬中，可以將線電流源剖分成多個點電流源，分別對點電流源產(chǎn)生的電位U做正演模擬，再將得到的結(jié)果進(jìn)行疊加處理，便得到線電流源的響應(yīng)。

3 有限差分法的格式推導(dǎo)

對三維求解區(qū)域進(jìn)行網(wǎng)格離散化處理，采用直角坐標(biāo)系，如圖2所示。

圖2 求解區(qū)域網(wǎng)格離散示意圖Fig.2 Schematic map of mesh grid discretization

利用有限差分法的原理：用差商來代替導(dǎo)數(shù)，可以寫出任意網(wǎng)格節(jié)點(i,j,k)處電位U的一階差商和二階差商分別為：

(2-1)

(2-2)

(2-3)

(3-1)

(3-2)

(3-3)

因此，對于線電流源L的任一線元dL，電流大小設(shè)為dI，式(1)可以寫為：

=-dIδ(i,j,k),∈Ω

(4-1)

(4-2)

(4-3)

利用邊界條件，在每一個網(wǎng)格節(jié)點處得到一個方程，聯(lián)立方程組，數(shù)值求解即可得到區(qū)域內(nèi)任意節(jié)點處的電位U。疊加所有線元dL正演得到的U，便可得到線電流源L在空間各處產(chǎn)生的電位大小。

4 三維模型的正演計算

利用上述的三維有限差分算法，計算均勻半空間條件下線電流源在地表產(chǎn)生的電位，與解析解之間的誤差大小如圖3所示。此檢驗?zāi)Ｐ偷木W(wǎng)格步長為1 m，網(wǎng)格區(qū)域大小為15 m×15 m×30 m，線源L位于地表中心處，長度為10 m。

從圖3中可以看出，除源點附近電位值有所畸變外，區(qū)域內(nèi)的其他測點誤差大小都在可以接受的范圍內(nèi)。因此可以認(rèn)為，本文采用的三維有限差分算法是正確的。

由于煤田的積水采空區(qū)往往不止一層，所以開展三維空間中多個異常體的異常響應(yīng)的研究是很有必要的?？紤]到積水區(qū)范圍有限，結(jié)合采空區(qū)積水的電性特征，設(shè)計了如下3個地質(zhì)模型：①單個低阻異常體，代表單個積水采空區(qū)；②垂向的兩個低阻異常體，代表不同深度的兩個積水采空區(qū)；③橫向的兩個低阻異常體，代表埋深相同、但橫向不同位置的2個積水采空區(qū)模型。分別研究在井中激發(fā)、地表接收時的視電阻率值，以此分析井地電位法對地下不同位置煤田積水采空區(qū)的識別能力。

圖3 數(shù)值解與解析解的誤差Fig.3 The rms error of numerical solution and analytical solution

在本文所有的模型中，背景電阻率均取為100 Ω·m，異常體電阻率為1 Ω·m，線源長度L為10 m，有限差分網(wǎng)格步長為1 m。模型①、②的異常體大小為5 m×5 m×1 m，網(wǎng)格區(qū)域大小為60 m×15 m×30 m，測距MN=1.414 m，指向為平行于y=x的方向。模型③的異常體大小為5 m×5 m×5 m，網(wǎng)格區(qū)域大小為60 m×15 m×30 m，測距MN=1 m，指向為平行于y=0的方向。

4.1 單個低阻異常體

圖4為均勻半空間中單個異常體模型的示意圖，圖中實線小長方體為異常體，垂直實線L代表放置于井中的垂直線源。圖5(a)、(b)分別為圖4中異常體埋深為3 m、5 m時的視電阻率等值線圖，圖5中黑色十字與方框分別代表垂直線源和異常體在地面上的投影位置。從圖5中可以看到對應(yīng)于視電阻率明顯偏低的區(qū)域與低阻異常體位置吻合較好，而隨著埋深的增加，異常響應(yīng)會有所減小。在低阻異常區(qū)域的兩端，出現(xiàn)了一個不太明顯的高阻假象，這可以利用視電阻率的微分形式來解釋[23]。同時由于源的奇異性，在源附近視電阻率值不準(zhǔn)確。

圖4 單個異常體模型示意圖Fig.4 Schematic map of one anomaly body model

圖5 模型①視電阻率的等值線Fig.5 The apparent resistivity contour map of the first model

4.2 垂向的兩個低阻異常體

圖6為垂向上兩個異常體模型的分布示意圖，來模擬地下三維空間中不同深度上的多個積水采空區(qū)情況。兩個異常體的埋深分別為5 m和10 m，水平方向頂點坐標(biāo)均為(10,10)。相應(yīng)的視電阻率等值線圖如圖7所示。從圖7可以看到，與圖5相同，在低阻體相應(yīng)位置有一個低阻異常區(qū)域，同時與只有一個埋深為5 m的異常體(圖5b)相比，低阻異常數(shù)值上略微增加，這是由于疊加了10 m的低阻異常體的響應(yīng)的原因。但僅從地表觀測的電位上，很難直接分辨埋深為5 m與10 m的兩個異常體。這說明，對于煤田中深度不同的多個積水采空區(qū)，井地電位法的分辨能力并不理想。

圖7 模型②視電阻率等值線Fig.7 The apparent resistivity contour map of the second model

圖6 垂向兩個異常體模型示意圖Fig.6 Schematic model of two anomaly bodies distributed vertically

圖8 橫向兩個異常體模型示意圖Fig.8 Schematic model of two anomaly bodies distributed transversely

圖9 模型③視電阻率等值線Fig.9 The apparent resistivity contour map of the third model

4.3 橫向的兩個低阻異常體

圖8為橫向兩個異常體模型的示意圖，以此來模擬地下空間中水平位置不同的積水采空區(qū)。兩個異常體的埋深均為5 m，水平方向頂點坐標(biāo)分別為(20,6)和(45,6)。其相應(yīng)的視電阻率等值線圖如圖9所示。從圖9可以看到，在兩個低阻體的對應(yīng)位置都有一個明顯的低阻異常區(qū)域，并且數(shù)值大小也非常接近。這一模型的正演結(jié)果表明，井地電位法對水平位置不同的多個積水采空區(qū)有較好的識別效果。

5 結(jié) 論

針對廣泛應(yīng)用于煤田積水采空區(qū)探測的井地電位法，本文利用有限差分法做了三維正演模擬研究，研究表明：井地電位法對一定深度范圍內(nèi)的積水采空區(qū)有良好的識別能力，但對垂向分布的多個積水采空區(qū)的分辨能力較弱，而對水平位置不同的多個積水采空區(qū)分辨能力較好。

[1]董書寧,虎維岳.中國煤礦水害基本特征及其主要影響因素[J].煤田地質(zhì)與勘探,2007,35(5):34-38.

[2]宋東日,譚湘龍,李杰.煤礦礦井水害的分析和防治方法[J].科技和產(chǎn)業(yè),2009,9(4):65-66.

[3]李蕓蕓.起伏地形井地電位三維有限元數(shù)值模擬[D].長沙:中南大學(xué),2014.

[4]Daniels J. Three-dimensional resistivity and induced-polarization modeling using buried electrodes[J]. Geophysics, 1977, 42(5):1 006-1 019.

[5]Sill W R, Ward S H. Electrical Energizing of Well Casings[M]. Utah: University of Utah, Department of Geology and Geophysics, 1978.

[6]Rocroi J P, Koulikov A V. The use of vertical line sources in electrical prospecting for hydrocarbon[J]. Geophysical Prospecting, 1985, 33(1): 138-152.

[7]張金成.電位法井間監(jiān)測技術(shù)[J].地震地質(zhì),2001,23(2):292-300.

[8]凌明友,郝新武,程同軍，等.井—地電位成像技術(shù)研究剩余油分布[J].斷塊油氣田,2003,10(4):55-58.

[9]譚河清,沈金松,周超，等.井—地電位成像技術(shù)及其在孤東八區(qū)剩余油分布研究中的應(yīng)用[J].石油大學(xué)學(xué)報(自然科學(xué)版),2004,28(2):31-37.

[10]何展翔,劉雪軍,裘尉庭，等.大功率井—地電法油藏邊界預(yù)測技術(shù)及效果[J].石油勘探與開發(fā),2004,31(5):74-76.

[11]王志剛,何展翔,魏文博，等.井地電法三維物理模型試驗[J].石油地球物理勘探,2005,40(5):594-597.

[12]劉昱,王志剛,何展翔.井地電法供電電場分布模擬研究[J].工程地球物理學(xué)報,2006,3(5):331-336.

[13]徐凱軍,李桐林.垂直有限線源三維地電場有限差分正演研究[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2006,36(1):137-141.

[14]劉地淵,徐凱軍,趙廣茂，等.任意形狀線電流源三維地電場研究[J].地球物理學(xué)進(jìn)展,2006,21(2):395-399.

[15]蘇朱劉,楊志冬,李宏，等.電位法在油儲動態(tài)監(jiān)測中的應(yīng)用[J].石油天然氣學(xué)報(江漢石油學(xué)院學(xué)報),2006,28(2):56-59.

[16]王志剛,何展翔,魏文博，等.井地電法的準(zhǔn)解析近似三維反演研究[J].石油地球物理勘探,2007,42(2):220-225.

[17]屈有恒,張貴賓,晉風(fēng)明.傾斜線源的三維電場數(shù)值模擬研究[J].物探化探計算技術(shù),2007,29(5):431-435.

[18]柯敢攀,黃清華.井地電法的三維正反演研究[J].北京大學(xué)學(xué)報(自然科學(xué)版),2009,45(2):264-272.

[19]戴前偉,陳德鵬,劉海飛，等.復(fù)雜管網(wǎng)線源井地電位三維有限元模擬[J].石油地球物理勘探,2012,47(6):998-1 005.

[20]李雙喜,曾昭發(fā),李靜，等.探地雷達(dá)極化探測時域有限差分法模擬效果分析[J].工程地球物理學(xué)報,2014,11(4):513-521.

[21]鄭冰,李柳德.高密度電法不同裝置的探測效果對比[J].工程地球物理學(xué)報,2015,12(1):33-39.

[22]徐正玉,楊海燕,鄧居智，等.回線源三維地—井瞬變電磁法FDTD數(shù)值模擬[J].工程地球物理學(xué)報,2015,12(3):327-332.

[23]李金銘.地電場與電法勘探[M].北京:地質(zhì)出版社,2005.

Three-dimensional Finite Difference Forward Modeling of Borehole-Surface Potential Method in Water Accumulating Gob in Coalfield

Qiu Weizhong1，Zhang Hui2,Tang Xingong2,Yan Liangjun2,Zhang Leijie1

(1.115CoalfieldGeologicalSurveyInstituteofShanxiProvince,DatongShanxi037039,China;2.SchoolofGeophysicsandOilResources,YangtzeUniversity,WuhanHubei430100,China)

Water accumulating gob in coalfield is one of the factors threatening mining safety. According to the special boundary value problem of borehole-surface potential method, the corresponding finite difference scheme is derived in this paper. Combining with the electrical characteristics and the distribution features of water-accumulating gob, a series of geo-electric models are designed to calculate their electric responses by using 3D finite difference approach. The forward modeling results show that the detection ability of borehole-surface potential method for a single water accumulating gob and multiple gobs which distributed transversely is good, while the detection ability for multiple gobs distributed vertically is relatively weak. The result has reference value for the detection of water accumulating gobs in coalfield.

Water accumulating gob; Borehole-surface potential method; Finite difference method; Three-dimensional forward modeling

1672—7940(2016)02—0149—06

10.3969/j.issn.1672-7940.2016.02.002

國家自然科學(xué)基金項目(編號：41274115，41274082，41404087，41504096)；973計劃項目(編號：2013CB228605)

邱衛(wèi)忠(1963-)，男，碩士，高級工程師，主要從事電磁法勘探與研究。E-mail：115qwz@sina.com

唐新功(1968-)，男，博士，教授，博士生導(dǎo)師，主要從事電磁法勘探，重磁勘探和巖石物理學(xué)的教學(xué)與研究。

E-mail：tangxingong@163.com

P631.3

A

2015-11-10