垂直及傾斜接觸面電測深曲線特征

湯井田， 譚 潔, 潘克家

(中南大學 地球科學與信息物理學院， 長沙 410083)

0 引言

直流電阻率法是以地殼中巖石和礦石的導電性差異為物理基礎，通過觀測與研究人工建立的地中穩定電流場的分布規律,進行找礦和解決地質問題的一組電法勘探分支方法，無論在普查金屬、非金屬礦產和研究地質構造方面，還是在水文、工程、環境地質調查以及勘查能源等方面，均取得了良好的地質效果，發揮著重要作用。

卡連諾夫[1]早在上世紀三十年代提出，當地下存在傾斜接觸面時，研究點源電流場的分布具有重要的理論與實踐意義。此后R.F.Aldredge[2], ，B.H.達連諾夫[3]等人也對該問題進行了研究，K.Maeda[4]利用點源多次鏡像進行求解，但僅適用于π是α的整數倍，或者(ρ2-ρ1)/(ρ2+ρ1)=-1(如ρ2=0)且π/2是α的整數倍，不及解析法的應用廣泛。

作者詳細推導了傾斜接觸面存在時電測深視電阻率的實用計算公式[5-11,14]，對進一步研究提供了一定的基礎，分析了對稱四極測深裝置[12-13]下不同的布極方位時的測深視電阻率曲線特征，對實際工作仍有一定的指導作用。

1 公式推導

利用解析法求解傾斜接觸面存在時的電位，是通過邊界條件求解與實際問題相應的拉普拉斯方程。

圖1為點源場中傾斜接觸面模型結構示意圖，A為點源，M為測點。采用柱坐標系，以斜面與地表的交線為Z軸，于是A點的柱坐標為(r0,0,0)。設任意一點M的柱坐標為(r,θ,z)，當0<θ<α時，M位于二面角α內；當α<θ<π時，M位于二面角α外。

圖1 點源場中斜接觸面模型結構示意圖Fig．1 Field of a point source situated on the earth above a dipping bed

問題在于求大地任一點M(r,θ,z)處的電位

U(r,θ,z)：

顯然電位滿足拉普拉斯方程:

(1)

U1、U2滿足以下邊界條件：

(2)U1=U2(θ=α)。

令U(r,θ,z)=P(θ)Q(r,z)，用分離變量法求解該定解問題可以得到：

(2)

(3)

為了便于編程計算，將解化為級數形式。

(1)當r>r0，且Res>-1/2時，交換積分次序，可改寫成式(4)。

Is(r0t)Ks(rt)ds

(4)

式中

g(s)=cos[s(π-θ)]+kcos[s(a-θ)]；h(s)=sin(sπ)+ksin(sa)。

(2)當r-1/2時，同理可得式(5)。

Is(rt)Ks(r0t)ds

(5)

U2的改寫同理可得。

設F(s)表示上面兩式中內層積分的被積函數。為了應用留數定理來計算積分，先考慮F(s)的奇點。由于g(s)與第一、第二類Bessel函數Is、Ks是s的解析函數，故F(s)的奇點就是h(s)的零點，即方程h(s)=0的根，該方程sin(sπ)+ksin(sa)=0(式中a=2α-π)依具體的k值及s的值而定。

數學上可以證明，在0<α<π的情況下，方程的根必為實數，故方程寫為h(x)=0。且|k|<1的情形下，h(x)=0的根必為單根。此時h(x)=sin(xπ)+ksin(xa)=0的根x滿足sin2(xπ)=k2sin2(xa)。

當0<|x|<1/2時，因0≤|α|<π，故0≤|xα|<|x|π<1/2。又在0≤φ<π/2時，(sinφ)′=cosφ>0，故sinφ是遞增函數。

于是sin|xα|

由此可見，在以原點為中心、以δ(<1/2)為半徑的圓內，僅有s=0是方程的單根，從而原點是F(s)的單極點，故從iδ到-iδ沿這圓的左半部分c積分：

(6)

(7)

同樣，

(8)

于是

(9)

設|k|<1的情形，將[-iym,iδ]與[iδ,iym]二線段以及半圓c與折線Lm接成的圍線記為Cm，據留數定理，依正向積分：

(10)

式中N是在Cm內h(x)=0的正根xi的個數，而

(11)

故：

(12)

令m→∞(從而ym→∞，N→∞)得到Cauchy主值：

(13)

將求和下標j改寫成m，可得：

(1)當r>r0時，

(14)

(2)當r

(15)

代入式(2)與式(3)得到：

(1)當r>r0時，

(16)

(2)當r

由于

(17)

則

(18)

(19)

方程h(x)=0的單根xm就是F(s)的單極點，也就是g(s)/h(s)及cos[s(π-θ)]/h(s)的單極點，在這些點處：

(20)

(21)

由此可見，當|k|<1時，計算公式為：

(22)

(23)

式中a=2α-π。

2 不同裝置下傾斜接觸面電測深曲線

在公式推導的基礎上，繪制了二極裝置下點源場中傾斜接觸面電測深曲線[19-21]，對垂直接觸面上不同布極方位電測深曲線特征進行了詳細討論。

2.1 二極裝置下傾斜接觸面電測深曲線

圖2為二極測深裝置的模型結構示意圖。模型參數：裝置軸與接觸面垂直，傾斜接觸面的角度分別取為30°、45°、60°以及90°，測點M橫坐標設為d，取為10 m，源點A橫坐標分別取10.2 m、10.3 m、10.4 m、10.5 m、10.6 m、10.7 m、10.8 m、10.9 m、11.0 m、12.0 m、13.0 m、14.0 m、15.0 m、16.0 m、17.0 m、18.0 m、19.0 m、20.0 m、25.0 m、30.0 m、40.0 m、50.0 m、60.0 m、70.0 m、80.0 m、100.0 m、150.0 m、200.0 m、300.0 m、400.0 m、500.0 m、600.0 m、800.0 m、1 000.0 m。

圖2 二極測深裝置模型結構示意圖Fig．2 The sketch of geoelectricity model with pole-pole array

二極裝置下，視電阻率定義為：

視電阻率曲線見圖3-圖6。

圖3 傾斜角度φ為30°，μ12分別取10、5、2、0.5、0.2、0.1Fig．3 Apparent resistivity curves for a dipping angle of 30° at constant values of μ12 as would be obtained by pole-pole array, where search electrodes are laid out in the down-dip direction

圖4 傾斜角度φ為45°，μ12分別取10、5、2、0.5、0.2、0.1Fig．4 Same as Fig.3 but for a dipping angle of 45°

圖5 傾斜角度φ為60°，μ12分別取10、5、2、0.5、0.2、0.1Fig．5 Same as Fig.3 but for a dipping angle of 60°

圖6 傾斜角度φ為90°，μ12分別取為10、5、2、0.5、0.2、0.1Fig．6 Same as but for a dipping angle of 90°

當點源跨過接觸面時，測點M橫坐標為10 m，源點A橫坐標為9.5 m、9.4 m、9.3 m、9.2 m、9.1 m、9.0 m、8.0 m、7.0 m、6.0 m、5.0 m、4.0 m、3.0 m、2.5 m、2.0 m、1.5 m、1.0 m、0.5 m、-0.5 m、-1.0 m、-2.0 m、-3.0 m、-4.0 m、-5.0 m、-6.0 m、-7.0 m、-8.0 m、-9.0 m、-10.0 m、-11.0 m、-12.0 m、-15.0 m、-20.0 m、-30.0 m、-50.0 m、-100.0 m、-200.0 m、-500.0 m，相應的曲線見圖7-圖10。

圖7 傾斜角度φ為30°，μ12分別取10、5、2、0.5、0.2、0.1Fig．7 Apparent resistivity curves for a dipping angle of 30° at constant values of μ12 as would be obtained by pole-pole array, where search electrodes are laid out perpendicular to the strike

圖8 傾斜角度φ為45°，μ12分別取為10、5、2、0.5、0.2、0.1Fig．8 Same as Fig.7 but for a dipping angle of 45°

圖9 傾斜角度φ為60°，μ12分別取10、5、2、0.5、0.2、0.1Fig．9 Same as Fig.7 but for a dipping angle of 60°

圖10 傾斜角度φ為90°，μ12分別取10、5、2、0.5、0.2、0.1Fig．10 Same as Fig.7 but for a dipping angle of 90°

由圖3-圖10可以看出，測深軸跨過接觸面時，在接觸面附近視電阻率曲線突變，而未跨過接觸面的曲線變化平緩，對比兩組曲線，可確定接觸面的位置；同時，接觸面兩側介質的電阻率相差越大，接觸面對視電阻率曲線的影響也越大。

2.2 不同布極方位下垂直接觸面的電測深曲線

本節重點討論了對稱四極裝置下垂直接觸面存在時的電測深曲線(圖11、圖12)，通過求解不同模型下的視電阻率，得到對稱四極裝置下視電阻率曲線隨各個參數的變化規律。

圖11為對稱四極測深裝置的模型示意圖,對稱四級裝置下視電阻率定義為：

圖11 對稱四極測深裝置模型結構示意圖Fig．11 The sketch of geoelectricity model with schlumberger array

圖12 對稱四極測深裝置在垂直接觸面上布極示意圖Fig．12 The sketch of geoelectricity model with schlumberger array over a vertical sheet

圖12的模型參數：垂直接觸面兩側電阻率之比ρ2/ρ1分別取為0.01、0.05、0.5、5、50、100；裝置軸與接觸面所成的角度分別取為30°、45°、60°以及90°，測點與接觸面的距離設為d，取為10 m。視電阻率曲見圖13-圖17。

圖13 φ=0°時視電阻率曲線，μ12分別取0.01、0.05、0.5、5、50、100Fig．13 Apparent resistivity curves for a intersection angle of 0° at constant values of μ12 as would be obtained by schlumberger array, where search electrodes are laid out perpendicular to the strike

圖14 φ=30°時視電阻率曲線，μ12分別取0.01、0.05、0.5、5、50、100Fig．14 Same as Fig.13 but for a intersection angle of 30°

圖15 φ=45°時視電阻率曲線，μ12分別取0.01、0.05、0.5、5、50、100Fig．15 Same as Fig.13 but for a intersection angle of 45°

圖16 φ=60°時視電阻率曲線，μ12分別取0.01、0.05、0.5、5、50、100Fig．16 Same as Fig.13 but for a intersection angle of 60°

圖17 φ=90°時視電阻率曲線，μ12分別取0.01、0.05、0.5、5、50、100Fig．17 Same as Fig.13 but for a intersection angle of 90°

從圖13-圖17可以看出，對稱四極裝置在垂直接觸面上的測深曲線，總體呈現先增大后趨于定值的趨勢。φ=0°即裝置軸AMNB平行接觸面時，曲線的左側漸近值為ρ1，右側漸近值為ρs=2ρ1ρ2/(ρ2+ρ1)，在AB/d>>20時，界面對測深視電阻率的影響已可忽略；垂直接觸面時，曲線在AB/d=2處畸變達到極大值，垂直界面位置對應于ρs極大(當ρ2/ρ1>1)或極小(ρ2/ρ1<1)；當ρ2/ρ1→∞時，ρs極大值等于1.45ρ1；當ρ2/ρ1=0時，ρs極大值等于0.55ρ1；右支漸近線ρs→ρ1，與ρ2/ρ1值無關；當布極方位與處置接觸面走向為不同夾角時，測深曲線形狀將介于上述兩種情況之間。

在同一垂直接觸面上，對比對稱四極裝置的裝置軸平行接觸面走向和垂直接觸面走向時的測深視電阻率曲線可見，因垂直接觸面的存在,給測深曲線帶來了一定的畸變，給實際異常曲線的解釋帶來了干擾。因此對其電測深理論正演曲線的研究，有助于對實測電測深曲線的分析認識和推斷解釋。顯然，垂直接觸面走向布極時,測深曲線的畸變比平行接觸面走向布極時復雜，在這種情形下，實際勘探中一般要求平行構造走向布極。

界面附加進行十字測深，由兩個方向ρs曲線的差異可以肯定垂直分界面的存在，而且可以大致確定其位置。

3 結論

通過作者對視電阻率曲線的分析及誤差對比可得到如下結論：

(1)平行接觸面時，曲線的左側漸近值為ρ1，右側漸近值為ρ∞=2ρ1ρ2/(ρ2+ρ1)；垂直接觸面時，曲線在AB/d=2處畸變達到極大值，垂直界面位置對應于ρs極大(ρ2/ρ1>1)或極小(ρ2/ρ1<1)；當ρ2/ρ1=∞時，ρs極大值等于1.45ρ1；當ρ2/ρ1=0時，ρs極大值等于0.55ρ1,右支漸近線ρ∞→ρ1，與ρ2/ρ1值無關。

(2)當布極方位與處置接觸面走向為不同夾角時，測深曲線形狀將介于上述兩種情況之間。

(3) 界面附加進行十字測深，由兩個方向ρs曲線的差異可以肯定垂直分界面的存在，而且可以大致確定其位置。

(4)實際勘探中一般要求平行構造走向布極。

參考文獻：

[1] 卡連諾夫.E.H.垂向電測深曲線解釋[M].北京：石油出版社，1959.

[2] ALDREDGE R F .The Effect of Dipping Strata on Earth Resistivity Determination[J].Colorado School of Mines Quart,1937，32(1)：169-186.

[3] B.H.達連諾夫.石油與天然氣電法勘探[M]. 林振民,譯,北京：地質出版社，1953.

[4] MAEDA K . Apparent resistivity curves for dipping bnd[J] .Geophysics,1955,120(1):123-139.

[5] 史磷華.三層傾斜層的K型地電斷面直流電測深的正反演問題[J].礦產及地質，1990，15( 4):63-68.

[6] 葛為中，阮百堯.直立多層介質地面點源電場的遞推公式[J].桂林冶金地質學院報，1990,10(2):173-181.

[7] 陳福萍，于安良.直立接觸面視電阻率量板[J].地質與勘探,1978(7):73-78.

[8] 劉象勤,李純志.直立巖層ρsηs測深理論曲線正演及程序[J].物探化探計算技術，1990, 12(2)：159-167.

[9] 劉象勤,李純志.二維直立巖層電阻率測深異常曲線的某些特征[J].礦產及地質，1992(3)：207-212.

[10] 屈超群.二度體地形對點源場影響的數學分析[J].應用數學學報,1977(2)：48-54.

[11] 屈超群，張靜.點源場計算方法[M].昆明：云南科學技術出版社，1999.

[12] 李金銘.地電場與電法勘探[M].北京：地質出版社，2005.

[13] 傅良魁.電法勘探教程[M].北京：地質出版社，1983.

[14] 王磊. 點電源場中垂直接觸帶地表節點電位求解分析[J].電腦開發與應用,2008，21(3)：50-52.

[15] 張秋光. 場論(中冊)[M]. 北京：地質出版社，1985.

[16] 王元明. 數學物理方程及特殊函數[M]. 北京：高等教育出版社，2004.

[17] 任懷宗，師先進.特殊函數概論[M].長沙：中南工業大學出版社，1986.

[18] 李壽貴，尹水仿.復變函數與積分變換[M].北京：科學出版社，2009.

[19] 劉衛國，易剛.C++程序設計實踐教程[M].北京：中國鐵道出版社，2008.

[20] IVOR HORTON. Visual C++2008入門經典[M].北京：清華大學出版社,2009.

[21] 羅延中 張桂青．電子計算機在電法勘探中的應用[M]．武漢：武漢地質學院出版社，1987．