王君恒， 耿煜

1 地下信息探測技術(shù)與儀器教育部重點(diǎn)實(shí)驗(yàn)室(中國地質(zhì)大學(xué)，北京)， 北京 100083 2 中國地質(zhì)大學(xué) 地質(zhì)過程與礦產(chǎn)資源國家重點(diǎn)實(shí)驗(yàn)室， 北京 100083 3 中國地質(zhì)大學(xué)(北京)地球物理與信息技術(shù)學(xué)院， 北京 100083 4 Center for Earthquake Research and Information, University of Memphis Memphis 38152, USA

?

自然電位梯度測量中的地形改正

王君恒1,2,3， 耿煜4*

1 地下信息探測技術(shù)與儀器教育部重點(diǎn)實(shí)驗(yàn)室(中國地質(zhì)大學(xué)，北京)， 北京 100083 2 中國地質(zhì)大學(xué) 地質(zhì)過程與礦產(chǎn)資源國家重點(diǎn)實(shí)驗(yàn)室， 北京 100083 3 中國地質(zhì)大學(xué)(北京)地球物理與信息技術(shù)學(xué)院， 北京 100083 4 Center for Earthquake Research and Information, University of Memphis Memphis 38152, USA

地形改正是重力和電法勘探中用于消除測點(diǎn)周圍地形起伏對觀測結(jié)果影響的技術(shù).由于地形形成自然電位異常的微觀機(jī)理較為復(fù)雜，解析公式難以得到，可以采用最為直觀的方法，尋找兩者之間的擬合公式.根據(jù)研究，在自然電位梯度測量中，地形引起的自然電位異常分為兩種：一是異常與地形起伏是鏡像關(guān)系，這種情況占多數(shù)；另一種情況是異常與地形起伏為反鏡像關(guān)系.多數(shù)地形起伏與異常存在線性關(guān)系、二次關(guān)系和指數(shù)關(guān)系，以線性關(guān)系為主.具體采用哪種擬合關(guān)系，可選取典型的地形和自然電位曲線，通過其變化規(guī)律進(jìn)行判斷；也可以采取試算的辦法，通過改正效果的比較選擇合適的擬合公式.本文選擇鄂爾多斯盆地作為研究地點(diǎn)，提出了適用于自然電位梯度測量的地形改正方法.通過三種擬合公式改正效果的比較，得出線性擬合公式能更好地描述該測區(qū)內(nèi)相對高差與自然電位改正量之間的關(guān)系.地形改正有效地消除了相對高差與自然電位之間的相關(guān)性，改正后的自然電位曲線能夠較好地反映地下地質(zhì)信息，自然電位等值線圖異常更加突出.自然電位梯度測量中，地形起伏可影響異常形態(tài)，特別是起伏地形的面積與目標(biāo)地質(zhì)體的面積大小相當(dāng)時(shí)，兩者異常有可能難以區(qū)分，因此地形改正的研究具有重要的實(shí)用性.

自然電位； 梯度測量； 最小二乘法； 地形改正

The Ordos Basin was selected to be the research site in this paper and the terrain correction method which is applicable to spontaneous potential gradient calculation was put forward. Due to many ridges and valleys, the relative elevation in this region is up to 250～300 m. Slopes are steep and regular, so that it is one of the best places within China to study terrain correction. Five measuring lines were deployed in east-west direction and seven measuring lines were deployed in north-south direction, respectively. Through the comparison of the correction results of three kinds of fitting formulas, it was concluded that the linear fitting formula can better describe the relationship between relative height and spontaneous potential correction amount in the surveyed area.

Because of the great impact of the terrain of this region, the spontaneous potential flat contour map before terrain correction shows significant positive anomalies that are associated with concave terrains and negative anomalies that are associated with raised terrains. Terrain correction can effectively eliminate the correlation between relative height and spontaneous potential. The corrected spontaneous potential curves can reflect subsurface geological structure better, and the anomalies shown by the spontaneous potential contour map can be more prominent.

The quality of terrain correction is evaluated by whether the correlation between terrain fluctuation and spontaneous potential could be eliminated. Terrain correction is effective when there is no significant mirror or anti-mirror relationship between terrain fluctuation and the spontaneous potential after terrain correction. Since terrain relief can influence the appearance of anomalies in spontaneous potential gradient calculation, especially when the area of rugged terrain is comparable to the area of oil and gas reservoirs, which may make it difficult to distinguish these two anomalies. Thus the study on terrain correction has great importance in practice.Keywords Spontaneous potential; Gradient calculation; Least square method; Terrain correction

1 引言

1.1 重力勘探的地形改正

地形改正一直被應(yīng)用地球物理工作者所關(guān)注，它是重力測量資料數(shù)據(jù)處理中最繁重、最重要的工作之一.物理大地測量學(xué)的理論與相關(guān)實(shí)踐表明：地球表層及內(nèi)部物質(zhì)的空間分布、運(yùn)動(dòng)和變化決定著地球重力場及其隨時(shí)間的變化.因此，精確確定地形起伏至關(guān)重要.

1912年，Hayford和Bowie、Hammer都在地形對重力測量的影響方面進(jìn)行過研究，提出了用同心圓和放射線將地形扇形分區(qū)進(jìn)行手算的重力地改方法，被稱為圓域法地改.隨著計(jì)算機(jī)的應(yīng)用，Kane (1962)等人不斷發(fā)展了適用于計(jì)算機(jī)計(jì)算的一系列地改方法，利用方形網(wǎng)格數(shù)字化地形高程數(shù)據(jù)、用方形楞柱體單元離散地形，有的稱其為方域法.大地水準(zhǔn)面外部不存在質(zhì)量是應(yīng)用Stokes公式計(jì)算的前提，這就要求大地水準(zhǔn)面以外地形質(zhì)量的影響必須被移去，因此局部地形改正的研究具有重要意義.周熙襄和李雄(1987)利用高斯公式提出了重力地改的表面積分法，得到了一個(gè)高精度計(jì)算重力地改的簡單公式.劉立言和辛舒暢(1992)將地球視為一表面起伏不規(guī)則的旋轉(zhuǎn)橢球體，對整個(gè)地球用三梭錐體元剖分，給出一種精確地形改正計(jì)算方法.用多面體精確公式計(jì)算地改，改正范圍任意直至全球.

赫爾默特凝聚校正、平衡校正和殘余地形模型(RTM)是地形校正和計(jì)算的實(shí)際方法.其中，赫爾默特凝聚校正修正大地水準(zhǔn)面，平衡校正和殘余地形模型(RTM)修正高度異常.寧津生等(1996)在某山區(qū)的試驗(yàn)表明：地面重力異常受地殼密度的橫向不均勻變化的影響在幾個(gè)毫伽，大地水準(zhǔn)面差距受到的影響可以達(dá)到厘米級.李?yuàn)檴櫤蛣⒀阌?2003)在對經(jīng)典赫爾默特校正、RTM和赫爾默特/RTM組合方法的研究中表明：赫爾默特/RTM組合方法和GPS/水準(zhǔn)數(shù)據(jù)擬合效果是最好的.王增利和文琳(2011)分別對地球表面和大地水準(zhǔn)面進(jìn)行物理受力分析，研究發(fā)現(xiàn)傳統(tǒng)的地形改正計(jì)算方法存在誤差，并給出了能有效改善地形校正計(jì)算精度的算法.黃謨濤、郭春喜、李建成等都就平面近似地形改正方法進(jìn)行過探討(郭春喜等, 2002; 李建成等, 2003; 龐振興, 2008; 章傳銀等, 2009, 2012; 郭東美, 2010).章傳銀認(rèn)為：傳統(tǒng)的平面布格改正破壞了重力場的調(diào)和性質(zhì)，只能用于地面重力異常，不宜推廣到地球表面以外的其他高度(章傳銀等, 2009, 2012).榮敏和周巍(2015)引入這種劃分方式，將大地水準(zhǔn)面視為球面，研究大地水準(zhǔn)面以上地形，并與傳統(tǒng)平面近似下的地形改正結(jié)果進(jìn)行了比較.

1.2 電法勘探的地形改正

在重力勘探中，地形起伏是由物質(zhì)的虧空和盈余所造成的.在明確微觀機(jī)理的前提下，地形異?？梢杂媒馕龉竭M(jìn)行表達(dá).與重力地改不同的是，電法勘探中的地形異常沒有這種解析公式，具體滿足何種關(guān)系要根據(jù)實(shí)際變化規(guī)律分析總結(jié).

電剖面法的觀測值受到地形起伏的不同影響.地形所引起的電場干擾存在三種表現(xiàn)形式：研究對象的異常全部消失，研究對象的某些異常特征消失和研究對象的異常位置移動(dòng).這些表現(xiàn)往往使人不能識別研究對象所引起的異常，或者使驗(yàn)證工作選錯(cuò)位置；此外，并不存在地質(zhì)對象的地點(diǎn)也會(huì)產(chǎn)生假異常.對于地形對電剖面法的影響，國內(nèi)外曾經(jīng)進(jìn)行了理論計(jì)算及模型實(shí)驗(yàn)工作.林崇遠(yuǎn)(1966)從1963年起先后分別進(jìn)行了實(shí)驗(yàn)工作，用土槽或薄水槽模擬野外的實(shí)際地形，用模擬地形上的觀測值對野外相應(yīng)測點(diǎn)的觀測值進(jìn)行改正，并取得了一些效果.

為了搞清視電阻率法的地形影響，國內(nèi)外一直在應(yīng)用模型實(shí)驗(yàn)和理論計(jì)算進(jìn)行研究.徐世浙(1966)推導(dǎo)并計(jì)算出有限長山谷和山脊的平面平行電場曲線，以及某些簡單光滑地形的線源場電流密度曲線.Maeda (1955)對無限長山脊和山谷的點(diǎn)源場電位公式有一些研究.葛為中(1977)提出用理論計(jì)算來消除復(fù)雜地形對視電阻率影響的方法.黃蘭珍等(1986)用解析法驗(yàn)證了角域地形上視電阻率數(shù)值解的正確性.實(shí)踐表明：用邊界元法獲得的地形改正的視電阻率曲線突出了目標(biāo)物異常；點(diǎn)源場電阻率法2D地形改正的邊界元解法是一種經(jīng)濟(jì)有效的方法.

理論研究和大量野外實(shí)踐現(xiàn)已充分證明：用比較法進(jìn)行地形改正可削弱地形干擾，突出有用異常，從而有效地提高電法勘探的地質(zhì)效果(何繼善和曾憲明, 1984).獲得純地形異常的方法目前大致分兩類：用保角映射或薄水層模型實(shí)驗(yàn)和二度電阻網(wǎng)絡(luò)來獲得線電源二度地形異常曲線；以及用純地形電位的解析表達(dá)式進(jìn)行計(jì)算，或者用物理模擬及數(shù)字模擬(如有限元法)獲得點(diǎn)源二度或三度地形異常曲線.

對于大地電磁的二維地形影響問題，徐世浙和趙生凱(1985)、Wannamaker等(1986)、Chouteau和Bouchard (1988)已用有限單元法進(jìn)行過研究.研究表明：Hx波受二維地形影響顯著，而Ex波受其影響較小.徐世浙(1995)用邊界單元法研究同一問題，得出了相同的結(jié)論.在此基礎(chǔ)上，徐世浙等(1997)簡化了Hx波二維地形影響的計(jì)算方法，實(shí)驗(yàn)表明邊界元法對Hx型波的二維地形改正是有效的.

充電法的地形影響來自兩個(gè)方面的因素(何裕盛和夏萬芳, 1978)：一是地形起伏改變了地下電流場的分布規(guī)律；二是在起伏地形上觀測改變了測點(diǎn)至場源的距離.何裕盛和夏萬芳(1978)、傅良魁(1983)討論了多種地形改正的方法來消除地形影響.戴光明和羅延鐘(1997)討論了充電法地形改正的已有方法存在的問題，基于比值法思想提出了一種地改方法.

2 自然電位梯度測量中的地形改正

在自然電位梯度測量的地形改正方法方面，目前尚無其他作者進(jìn)行研究.根據(jù)本文作者的研究，地形起伏可影響異常形態(tài)，特別是起伏地形的面積與地質(zhì)體的面積大小相當(dāng)時(shí)，兩者異常有可能難以區(qū)分，需做地形改正(王君恒等, 2000, 2003; 厲玉樂等, 2005).

2.1 梯度測量方法流程

(1)

其中V代表兩點(diǎn)間電位差，單位為mV.上一個(gè)點(diǎn)的V23與下一個(gè)點(diǎn)的V12近似相等，允許存在1 mV左右的誤差.

由于觀測得到是兩點(diǎn)之間的電位差，它不能反映地下的情況，需要把電位差轉(zhuǎn)化為電位.具體的處理方法為

(2)

其中Ui為疊加后的自然電位，單位為mV；a為測線起點(diǎn)處的自然電位，可設(shè)為0 mV或某一特定數(shù)值.

經(jīng)過(2)式，實(shí)測數(shù)據(jù)所對應(yīng)的V23和V12的平均值作為測點(diǎn)間的電位差，該電位差再進(jìn)行疊加得到每點(diǎn)的自然電位.處理時(shí)，需按照極差校正→數(shù)據(jù)疊加→日變改正→地形改正的工作流程.首先消除極差影響，然后對測量數(shù)據(jù)進(jìn)行疊加從而將電位差轉(zhuǎn)化為電位，再消除日變影響，最后在極差校正和日變改正的基礎(chǔ)上進(jìn)行地形改正.

2.2 地形改正原理

如圖1所示，在自然電位梯度測量中，地形異常分為兩種：一是異常與地形影響是鏡像關(guān)系，這種情況占多數(shù)；另一種情況是異常與地形為反鏡像關(guān)系，在油氣勘探中少見(Wang et al., 2007; 王君恒等, 2007).多數(shù)地形起伏與異常存在三種關(guān)系：線性關(guān)系、二次關(guān)系或指數(shù)關(guān)系，可利用其進(jìn)行改正，關(guān)系式為

ΔUi=a0+a1ΔHi，

(3)

(4)

ΔUi=AeBΔHi+C，

(5)

其中ΔUi為電位改正值，單位為mV；ΔHi為地形相對高差，單位為m；a0、a1、a2、A、B、C為無量綱常系數(shù).具體滿足何種關(guān)系，要根據(jù)觀察的實(shí)際變化規(guī)律決定.有少數(shù)地形起伏與異常關(guān)系復(fù)雜，可用相關(guān)分析的方法進(jìn)行改正.

圖1 地形起伏電位示意圖Fig.1 The potential caused by topography

數(shù)據(jù)處理時(shí)，可得到記錄高程的點(diǎn)所對應(yīng)的自然電位.改正時(shí)，首先繪制出每條測線的相對高差曲線和自然電位曲線.從相對高差曲線中選取一段變化幅度較大的地形，通過測點(diǎn)號找到其對應(yīng)的自然電位.這樣，就得到了一組相對高差{ΔH1, ΔH2, ΔH3, …, ΔHn}及其對應(yīng)的自然電位{ΔU1, ΔU2, ΔU3, …, ΔUn}，n為所取擬合數(shù)據(jù)的個(gè)數(shù).

如果設(shè)自然電位是關(guān)于相對高差的線性函數(shù)，即ΔU=f(ΔH)=a0+a1ΔH.應(yīng)用最小二乘擬合公式為

(6)

求出系數(shù)a0和a1的值，便得到自然電位改正量關(guān)于相對高差的線性擬合公式(2.3).

如果設(shè)自然電位是關(guān)于相對高差的二次函數(shù)，即ΔU=f(ΔH)=a0+a1ΔH+a2ΔH2.應(yīng)用最小二乘擬合公式為

(7)

求出系數(shù)a0、a1和a2的值，便得到自然電位改正量關(guān)于相對高差的二次擬合公式(4).

如果設(shè)自然電位是關(guān)于相對高差的指數(shù)函數(shù)，即ΔU=f(ΔH)=AeBΔH+C，首先將等號兩邊同時(shí)取對數(shù)，可得

lnΔU=ln(AeBΔH+C)=lnA+BΔH+C，

(8)

又因?yàn)閘nA和C可合并為同一常數(shù)，令lnA=a0、B=a1、C=0，于是得

lnΔU=a0+a1ΔH，

(9)

(10)

通過該式計(jì)算出常數(shù)a0和a1的值，并利用A=ea0、B=a1、C=0，便得到自然電位改正量關(guān)于相對高差的指數(shù)擬合公式(5).應(yīng)用式(3)、(4)和(5)可計(jì)算出測區(qū)內(nèi)每一測點(diǎn)的自然電位改正量，再通過該改正量對自然電位曲線進(jìn)行地形改正.

3 鄂爾多斯盆地的地形改正

3.1 工區(qū)概況

工區(qū)仍屬黃土塬區(qū)，區(qū)內(nèi)有數(shù)塊面積較大的黃土殘塬，其余地區(qū)被雨水沖刷切割成樹枝狀水系.大溝中出露白堊系，支溝中為第三系紅土或第四系黃土層，黃土層及紅土總厚約100～150余米.溝峁梁塬縱橫交錯(cuò)，相對高差可達(dá)250～300 m.斜坡坡度陡峭且形態(tài)規(guī)整，是我國境內(nèi)研究地形改正的最佳地點(diǎn)之一.在東西向部署了5條測線，南北向部署了7條測線.

3.2 擬合數(shù)據(jù)的選取

可用于地形改正的擬合數(shù)據(jù)需遵照以下原則進(jìn)行選?。?1)地形起伏較大；(2)地形引起的自然電位異常非常明顯；(3)起伏地形與異常形態(tài)之間存在較好的對應(yīng)關(guān)系.N1號測線是一條地形起伏較大的測線，測線近似南北向分布，中途經(jīng)過了兩道凹陷的地帶.設(shè)1300 m為該地區(qū)的平均高程，選取高程變化最大的一段數(shù)據(jù)，繪制其相對高差曲線和自然電位曲線，如圖2所示.

圖2 N1號測線地形起伏較大段相對高差變化和自然電位異常鏡像關(guān)系Fig.2 The mirror image relationship between the relative height curve and the spontaneous potential curve of the large undulating terrain section of N1

可以看到：起伏的地形引起了顯著的自然電位變化，對地下地質(zhì)信息的反映造成了較大干擾.兩者之間大致呈現(xiàn)鏡像關(guān)系，地形每下降200 m，引起的自然電位正異常約為15 mV.其中，測點(diǎn)號310到326、測點(diǎn)號340到352、測點(diǎn)號374到392、測點(diǎn)號400到418之間的地形坡度最大、引起的自然電位異常最為明顯，可作為最小二乘擬合所用數(shù)據(jù).應(yīng)用同樣的方法從該測區(qū)內(nèi)其他幾條測線中選取擬合數(shù)據(jù)后，得到了一系列相對高差和與之對應(yīng)的自然電位，如表1所示.

表1 用于擬合的相對高差和自然電位原始數(shù)據(jù)Table 1 The relative height and their corresponding spontaneous potential that are used for fitting

3.3 三種擬合公式的改正效果

若假設(shè)相對高差與自然電位之間滿足線性函數(shù)關(guān)系，應(yīng)用最小二乘擬合公式(6)，可計(jì)算出參數(shù)的值a0= -2.729、a1= -0.06497，于是得自然電位改正量與相對高差之間的線性擬合公式為

ΔUi=-2.729-0.06497ΔHi.

(11)

若假設(shè)相對高差與自然電位之間滿足二次函數(shù)關(guān)系，應(yīng)用最小二乘擬合公式(7)，可計(jì)算出參數(shù)的值a0=-3.901、a1=-0.1099和a2=-0.0002387，于是得自然電位改正量與相對高差之間的二次擬合公式為

(12)

若假設(shè)相對高差與自然電位之間滿足指數(shù)函數(shù)關(guān)系，應(yīng)用最小二乘擬合公式(10)，可計(jì)算出參數(shù)的值a0=-0.08927、a1=-0.01193，亦即A=0.9146、B=-0.01193、C=0，于是得自然電位改正量與相對高差之間的指數(shù)擬合公式為

ΔUi=0.9146e-0.01193ΔHi.

(13)

擬合曲線如圖3所示.

圖3 應(yīng)用最小二乘公式(6)、(7)和(10)所實(shí)現(xiàn)的線性擬合、二次擬合及指數(shù)擬合曲線Fig.3 The linear fitting, quadratic fitting and exponential fitting curves implemented by applying the least squares formula (6), (7) and (10)

從圖中可以看出：由于在二次擬合公式(12)中二次項(xiàng)的系數(shù)a2很小，因此二次擬合與線性擬合兩條曲線非常接近.指數(shù)擬合曲線與線性擬合曲線在相對高差 (-200 m, -50 m)的區(qū)間內(nèi)較為接近；當(dāng)相對高差在-50 m以上時(shí)，指數(shù)擬合公式的地形改正量接近于0.以N5號測線地形起伏較大的一段曲線為例，分析并比較三種不同擬合公式的改正效果.測點(diǎn)號250之后的部分相對高差約為-250 m.分別應(yīng)用線性、二次和指數(shù)擬合公式(11)、(12)和(13)，對測線進(jìn)行改正，改正效果如圖4所示.

圖4 分別應(yīng)用線性擬合公式、二次擬合公式和指數(shù)擬合公式對N5號測線地形起伏較大段所做地形改正前、后自然電位曲線(虛線：改正前自然電位曲線；實(shí)線：改正后自然電位曲線)Fig.4 The spontaneous potential curves of the large undulating terrain section of N5 before and after topography correction by applying the linear fitting formula, the quadratic fitting formula and the exponential fitting formula respectively(dashed line is the spontaneous potential curve before topography correction; solid line is the spontaneous potential curve after topography correction)

可以看到：受下降地形的影響，地形改正前N5號測線自然電位曲線測點(diǎn)號250之后的部分顯著向上抬升.對于測點(diǎn)號250之后的部分，三種擬合公式的改正效果無太大差別；對于測點(diǎn)號250之前的部分，線性擬合公式的改正效果最為顯著，指數(shù)擬合公式幾乎沒有改正效果.這是由于當(dāng)相對高差約為-250 m時(shí)，三種公式的擬合曲線非常接近(圖3)；而當(dāng)相對高差大于-50 m時(shí)，指數(shù)擬合公式的地形改正量接近于0(圖3).綜上所述：能更好地描述該段測線地形起伏與自然電位之間關(guān)系的是線性擬合公式.

3.4 測區(qū)的地形改正

對測區(qū)內(nèi)全部12條測線(東西向5條+南北向7條)進(jìn)行了地形改正，繪制出該測區(qū)的相對高差等值線圖、改正前自然電位等值線圖和改正后自然電位等值線圖，如圖5—圖7所示.

圖5 該測區(qū)的相對高差平面等值線(測量點(diǎn)位以黑色十字標(biāo)出)Fig.5 The flat contour map of relative height in the surveyed area (point positions have been marked with black crosses)

圖6 僅經(jīng)過極差校正和日變改正的自然電位平面等值線(測量點(diǎn)位以黑色十字標(biāo)出)Fig.6 The flat contour map of spontaneous potential after contrast correction and diurnal correction (point positions have been marked with black crosses)

圖7 經(jīng)過極差校正、日變改正和地形改正的自然電位平面等值線(測量點(diǎn)位以黑色十字標(biāo)出)Fig.7 The flat contour map of spontaneous potential after contrast correction, diurnal correction and topography correction (point positions have been marked with black crosses)

由圖可見：由于該地區(qū)的地形起伏影響較大，地形改正前的自然電位平面等值線圖中呈現(xiàn)出顯著的與凹陷地形相關(guān)的正異常(圖6中紅色線框標(biāo)注)及與突起地形相關(guān)的負(fù)異常(圖6中藍(lán)色線框標(biāo)注).地形改正有效地消除了相對高差與自然電位之間的相關(guān)性，從而排除了地形因素給地下地質(zhì)信息反映所帶來的干擾.

4 結(jié)論

(1) 在自然電位梯度測量中，地形異常分為兩種：一是異常與地形影響是鏡像關(guān)系；另一種情況是異常與地形為反鏡像關(guān)系.多數(shù)地形起伏與異常存在一種近似線性、二次或指數(shù)關(guān)系，可利用該關(guān)系進(jìn)行擬合和改正.

(2) 由于地形產(chǎn)生自然電位異常的機(jī)理尚不清楚，通過解析方法推導(dǎo)出的地改公式難以得到，可以通過更加直觀的方法尋找地形與自然電位之間的擬合公式.本文給出了相對高差與自然電位改正量之間的線性擬合公式、二次擬合公式和指數(shù)擬合公式.

(3) 自然電位梯度測量中的地形改正，是通過先假定公式、再計(jì)算公式中系數(shù)的方法求得的.具體采用哪種擬合關(guān)系，可以先選取典型的地形及其對應(yīng)的自然電位曲線，通過其變化規(guī)律進(jìn)行推斷；也可以采取試算的辦法，通過改正效果的比較選取最為合適的擬合公式.本文對比了三種擬合公式的改正效果，得出線性擬合公式能更好地描述該工區(qū)的相對高差與自然電位改正量之間的關(guān)系.

(4) 地形改正的優(yōu)劣評價(jià)取決于能否消除地形起伏與自然電位之間的相關(guān)性.若改正后的自然電位與地形起伏明顯不存在鏡像關(guān)系(或反鏡像關(guān)系)，則說明地形改正是有效的.

Bian S F, Zhang K F. 1991. On the role of terrain in the refinement of the geoid in China.ProgressinGeophysics(in Chinese), 6(4): 64-74. Chouteau M, Bouchard K. 1988. Two-dimensional terrain correction in magnetotelluric surveys.Geophysics, 53(6): 854-862.Dai G M, Luo T Z. 1997. A way of terrain correction for mise—a-la-masse method in 2 dimention.ComputingTechniquesforGeophysicalandGeochemicalExploration(in Chinese), 19(1): 36-40.

Fu L K. 1983. Electrical Prospecting Tutorials (in Chinese). Beijing: Geological Publishing House.

Ge W Z. 1977. Topographic correction for the DC apparent resistivity method of electrical prospecting.ActaGeophysicaSinica(in Chinese), 20(4): 299-311.

Gou C X, Wang H M, Wang B. 2002. Determination of high resolution grid terrain and isostatic corrections in all China area.ActaGeodaeticaetCartographicaSinica(in Chinese), 31(3): 201-205.

Guo D M. 2010. Refined theories, methods and applications of local quasi geoid[Ph. D. thesis] (in Chinese). Beijing: Graduate School of Chinese Academy of Sciences.

Guo D M, Xu H Z. 2011. Research on the singular integral of local terrain correction computation.ChineseJ.Geophys. (in Chinese), 54(4): 977-983, doi: 10.3969/j.issn.0001-5733.2011.04.012.He J S, Zeng X M. 1984. Fast calculation using polynomes for topographic correction in resistivity method.GeophysicalandGeochemicalExploration(in Chinese), 8(1): 27-33.

He Y S, Xia W F. 1978. The Excitation-at-the-Mass Method (in Chinese). Beijing: Geological Publishing House.

Hofmann-Wellenhof B, Moritz H. 2006. Physical Geodesy. Wien: Springer Science & Business Media.

Huang L Z, Tian X M, Cun S C. 1986. The boundary element method of the 2-D topographic correction for the resistivity method in the point source field.ComputingTechniquesforGeophysicalandGeochemicalExploration(in Chinese), 8(3): 201-208.

Huang M T, Zhai G J, Guan Z, et al. 2000. On computations of terrain corrections and indirect effect using FFT.JournalofInstituteofSurveyingandMapping(in Chinese), 17(4): 242-246.

Kane M F. 1962. A comprehensive system of terrain corrections using a digital computer.Geophysics, 27(4): 455-462.

Li J C, Chen J Y, Ning J S, et al. 2003. The Approximation Theory of the Earth′s Gravitational Field and the Determination of China 2000 Quasi-Geoid (in Chinese). Wuhan: Wuhan University Press.

Li S S, Liu Y Y. 2003. Calculation of the reference orbit of CHAMP by gravitational potential model.JournalofInstituteofSurveyingandMapping(in Chinese), 20(3): 165-167.

Li Y C. 1993. Optimized spectral geoid determination. UCSE Report No 20050. Calgary, Alberta, Canada: Geomatics Engineering, University of Calgary.Li Y L, Wang J H, Pan Z P, et al. 2005. Wavelet multiscale analysis and geophysical significance of the spontaneous potential of oil and gas reservoir.PetroleumExplorationandDevelopment(in Chinese), 32(1): 80-83.

Lin C Y. 1966. The model experimental method of topographic correction in electric profiling method.ActaGeophysicaSinica(in Chinese), 15(1): 66-72.

Lin Z M, Shi Z S. 1984. A discussion on several terrain correction methods in regional gravity survey.GeophysicalandGeochemicalExploration(in Chinese), 8(4): 193-198.Liu L Y, Xin S C. 1992. A topographical correction method for gravity measurement based on a triangular polyhedron.ComputingTechniquesforGeophysicalandGeochemicalExploration(in Chinese), 14(4): 293-300.Lu Z L. 1989. Improving the Stokes solution of the disturbing potential with the help of topographic data.ActaGeodaeticaetCartographicaSinica(in Chinese), 18(1): 60-71.

Lu Z L. 1996. Theories and Methods of the Earth′s Gravitational Field (in Chinese). Beijing: The People′s Liberation Army Press.

Maeda K. 1955. Apparent resistivity for dipping beds.Geophysics, 20(1): 123-139.

Ning J S, Luo Z C, Chao D B. 1996. The present situation on satellite gravity gradiometry and its vistas in the application of physical geodesy.JournalofWuhanTechnicalUniversityofSurveyingandMapping(in Chinese), 21(4): 309-314.

Pang Z X. 2008. Research on the refining methods of gravity anomaly[Master′s thesis] (in Chinese). Zhengzhou: The PLA Information Engineering University.

Rong M, Zhou W. 2015. Study on topography correction based on spherical approximation.JournalofGeodesyandGeodynamics(in Chinese), 35(1): 58-61.

Schwarz K P, Sideris M G, Forsberg R. 1990. The use of FFT techniques in physical geodesy.GeophysicalJournalInternational, 100(3): 485-514.Sideris M G. 1985. A fast Fourier transform method for computing terrain corrections.ManuscriptaGeodaetica, 10(1): 66-73.

Sideris M G. 1990. Rigorous gravimetric terrain modelling using Molodensky′s operator.ManuscriptaGeodaetica, 15(2): 97-106.

Sideris M G, Li Y C. 1993. Gravity field convolutions without windowing and edge effects.BulletinGéodésique, 67(2): 107-118.

Wang J H, Pan Z P, Guan Z N, et al. 2000. Self-potential detection principle of oil/gas reservoir and its application to field development.PetroleumExplorationandDevelopment(in Chinese), 27(3): 96-99, 102.Wang J H, Pan Z P, Deng M S, et al. 2003. Relationship between spontaneous potential (SP) and apparent effective thickness (AET) of reservoirs.PetroleumExplorationandDevelopment(in Chinese), 30(5): 65-67.

Wang J H, Pan Z P, Sun S W, et al. 2007. Apparent effective thickness prevision through spontaneous potential method and its application in oil development.EarthScience—JournalofChinaUniversityofGeosciences(in Chinese), 32(4): 461-468.

Wang J H, Pan Z P, Sun S W, et al. 2007. Prediction of apparent equivalent thickness using the spontaneous potential method and its application to oilfield development.JournalofChinaUniversityofGeosciences, 18(3): 269-279.

Wang Z L, Wen L. 2011. A new terrain correction method.JournalofGeodesyandGeodynamics(in Chinese), 31(3): 115-119.

Wannamaker P E, Stodt J A, Rijo L. 1986. Two-dimensional topographic responses in magnetotellurics modeled using finite elements.Geophysics, 51(11): 2131-2144.

Wichiencharoen C. 1982. The indirect effects on the computation of geoid undulations. Report No. 336. Columbus, OH ,United States: Dept. of Geodetic Science and Surveying, The Ohio State University.

Xu H Z. 2006. Some problems on the precise determination of regional geoid in China.GeospatialInformation(in Chinese), 4(5): 1-3.

Xu S Z. 1966. The impact of two-dimensional terrain on resistivity method.ActaGeophysicaSinica(in Chinese), 15(1): 73-82.

Xu S Z, Zhao S K. 1985. The topographic effects on magnetotelluric response.NorthwesternSeismologicalJournal(in Chinese), 7(4): 69-78.

Xu S Z. 1995. The Boundary Element Method in Geophysics (in Chinese). Beijing: Science Press.

Xu S Z, Li Y G, Liu B. 1997. The method and efficiency of 2-D terrain correction forHX-polarization of MT.ActaGeophysicaSinica(in Chinese), 40(6): 842-846.

Zhang C Y, Chao D B, Ding J, et al. 2006. Arithmetic and characters analysis of terrain effects for CM-order precision height anomaly.ActaGeodaeticaetCartographicaSinica(in Chinese), 35(4): 308-314.

Zhang C Y, Chao D B, Ding J, et al. 2009. Precision topographical effects for any kind of field quantities for any altitude.ActaGeodaeticaetCartographicaSinica(in Chinese), 38(1): 28-34.Zhang C Y, Dang Y M, Ke B G, et al. 2012. Discussion of some lssues in the coastal precise quasi-geoid determination.GeomaticsandInformationScienceofWuhanUniversity(in Chinese), 37(10): 1150-1154.

Zhou D Y. 2005. An innovation of 3-D non-singular boundary integral equations (in Chinese). Taiwan: National Taiwan University.

Zhou X X, Li X. 1987. Gravimetric terrain corrections by surface integration method.ComputingTechniquesforGeophysicalandGeochemicalExploration(in Chinese), 9(4): 273-279.

附中文參考文獻(xiàn)

邊少鋒, 張克非. 1991. 試論地形在我國大地水準(zhǔn)面精化中的作用. 地球物理學(xué)進(jìn)展, 6(4): 64-74.

戴光明, 羅延鐘. 1997. 二維地電條件下充電法地形改正的一種方法. 物探化探計(jì)算技術(shù), 19(1): 36-40.

傅良魁. 1983. 電法勘探教程. 北京: 地質(zhì)出版社.

葛為中. 1977. 視電阻率法的地形改正. 地球物理學(xué)報(bào), 20(4): 299-311.

郭春喜, 王惠民, 王斌. 2002. 全國高分辨率格網(wǎng)地形和均衡改正的確定. 測繪學(xué)報(bào), 31(3): 201-205.

郭東美. 2010. 局部擬大地水準(zhǔn)面精化理論、方法及應(yīng)用[博士論文]. 北京: 中國科學(xué)院研究生院.

郭東美, 許厚澤. 2011. 局部地形改正的奇異積分研究. 地球物理學(xué)報(bào), 54(4): 977-983, doi: 10.3969/j.issn.0001-5733.2011.04.012.

何繼善, 曾憲明. 1984. 電阻率法地形改正的多項(xiàng)式快速計(jì)算方法. 物探與化探, 8(1): 27-33.

何裕盛, 夏萬芳. 1978. 充電法. 北京: 地質(zhì)出版社.

黃蘭珍, 田憲謨, 寸樹蒼. 1986. 點(diǎn)源場電阻率法二維地形改正的邊界元法. 物探化探計(jì)算技術(shù), 8(3): 201-208.

黃謨濤, 翟國君, 管錚等. 2000. 利用FFT技術(shù)計(jì)算地形改正和間接效應(yīng). 測繪學(xué)院學(xué)報(bào), 17(4): 242-246.

李建成, 陳俊勇, 寧津生等. 2003. 地球重力場逼近理論與中國2000似大地水準(zhǔn)面的確定. 武漢: 武漢大學(xué)出版社.

李?yuàn)檴? 劉雁雨. 2003. 利用地球重力場模型計(jì)算CHAMP衛(wèi)星參考軌道. 測繪學(xué)院學(xué)報(bào), 20(3): 165-167.

厲玉樂, 王君恒, 潘竹平等. 2005. 油氣藏自然電位小波多尺度分解及其地球物理意義. 石油勘探與開發(fā), 32(1): 80-83.

林崇遠(yuǎn). 1966. 電剖面法地形改正的模型實(shí)驗(yàn)方法. 地球物理學(xué)報(bào),15(1): 66-72.

劉立言, 辛舒暢. 1992. 三角面多面體重力地形改正方法. 物探化探計(jì)算技術(shù), 14(4): 293-300.

陸仲連. 1989. 借助地形數(shù)據(jù)改進(jìn)擾動(dòng)位的司托克斯解. 測繪學(xué)報(bào), 18(1): 60-71.

陸仲連. 1996. 地球重力場理論與方法. 北京: 解放軍出版社.

寧津生, 羅志才, 晁定波. 1996. 衛(wèi)星重力梯度測量的研究現(xiàn)狀及其在物理大地測量中的應(yīng)用前景. 武漢測繪科技大學(xué)學(xué)報(bào), 21(4): 309-314.

龐振興. 2008. 平均空間重力異常精細(xì)構(gòu)制方法研究[碩士論文]. 鄭州: 解放軍信息工程大學(xué).

榮敏, 周巍. 2015. 球近似地形改正的研究分析. 大地測量與地球動(dòng)力學(xué), 35(1): 58-61.

王君恒, 鄧明勝, 潘竹平等. 2000. 油氣藏自然電位勘探原理及在開發(fā)中的應(yīng)用. 石油勘探與開發(fā), 27(3): 96-99, 102.

王君恒, 潘竹平, 鄧明勝等. 2003. 油氣藏自然電位與油層視有效厚度的關(guān)系研究. 石油勘探與開發(fā), 30(5): 65-67.

王君恒, 潘竹平, 孫樹文等. 2007. 自然電位視有效厚度預(yù)測及在油田開發(fā)中的應(yīng)用. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào), 32(4): 461-468.

王增利, 文琳. 2011. 一種地形改正新算法. 大地測量與地球動(dòng)力學(xué), 31(3): 115-119.

許厚澤. 2006. 我國精化大地水準(zhǔn)面工作中若干問題的討論. 地理空間信息, 4(5): 1-3.

徐世浙. 1966. 二度地形對電阻率法的影響. 地球物理學(xué)報(bào), 15(1): 73-82.

徐世浙, 趙生凱. 1985. 地形對大地電磁勘探的影響. 西北地震學(xué)報(bào), 7(4): 69-78.

徐世浙. 1995. 地球物理中的邊界單元法. 北京: 科學(xué)出版社.

徐世浙, 李予國, 劉斌. 1997. 大地電磁HX型波二維地形改正的方法與效果. 地球物理學(xué)報(bào), 40(6): 842-846.

章傳銀, 晁定波, 丁劍等. 2006. 厘米級高程異常地形影響的算法及特征分析. 測繪學(xué)報(bào), 35(4): 308-314.

章傳銀, 晁定波, 丁劍等. 2009. 球近似下地球外空間任意類型場元的地形影響. 測繪學(xué)報(bào), 38(1): 28-34.

章傳銀, 黨亞民, 柯寶貴等. 2012. 高精度海岸帶重力似大地水準(zhǔn)面的若干問題討論. 武漢大學(xué)學(xué)報(bào)·信息科學(xué)版, 37(10): 1150-1154.

Chou T Y. 以非奇異性暇積分改善三維邊界積分法. 臺灣大學(xué)工程科學(xué)及海洋工程學(xué)研究所學(xué)位論文, 2005: 1-179.

周熙襄, 李雄. 1987. 重力地形改正的表面積分法. 物探化探計(jì)算技術(shù), 9(4): 273-279.

(本文編輯 張正峰)

Terrain correction in the gradient calculation of spontaneous potential data

WANG Jun-Heng1,2,3， GENG Yu4*

1KeyLaboratoryofGeo-detection(ChinaUniversityofGeosciences,Beijing),MinistryofEducation,Beijing100083,China2ChinaUniversityofGeosciences,StateKeyLaboratoryofGeologicalProcessesandMineralResources,Beijing100083,China3SchoolofGeophysicsandInformationTechnology,ChinaUniversityofGeosciences,Beijing100083,China4CenterforEarthquakeResearchandInformation,UniversityofMemphisMemphis38152,USA

Terrain correction is to eliminate the effect of terrain relief around measurement points on observation results in gravitational and geoelectrical prospecting. Because the microscopic mechanism of spontaneous potential anomalies formed by terrains is rather complex, it is difficult to obtain an analytic formula, whereas the most visualized method can be adopted instead to find out the fitting formula between these two variables. According to the study, in spontaneous potential gradient calculation, spontaneous potential anomalies caused by topography can be divided into two kinds: one is that anomalies have a mirror-image relation with terrain relief, which accounts for the majority of cases; the other is that anomalies have an inverted mirror-image relation with terrain relief. The relation between a majority of anomalies and terrain relief could be linear, quadratic or exponential, among which linear relation is the primary. As for which fitting relation should be adopted specifically, it can be decided by selecting a typical terrain and spontaneous potential curve first and judging through its change rule; and trial calculation method can be adopted as well to choose the appropriate fitting formula by comparing the correction effects.

10.6038/cjg20151031.

Wang J H, Geng Y. 2015. Terrain correction in the gradient calculation of spontaneous potential data.ChineseJ.Geophys. (in Chinese),58(10):3815-3824,doi:10.6038/cjg20151031.

王君恒，男，1962年生，山東青島人，博士，中國地質(zhì)大學(xué)(北京)地球物理與信息技術(shù)學(xué)院副教授，主要從事應(yīng)用地球物理和理論地球物理研究.E-mail: w1128@cugb.edu.cn

*通訊作者 耿煜，男，1989年生，天津人，碩士，孟菲斯大學(xué)地震研究信息中心研究生，主要從事天然地震研究.E-mail: ygeng1@memphis.edu

10.6038/cjg20151031

P631

2015-08-04，2015-10-08收修定稿

王君恒，耿煜. 2015. 自然電位梯度測量中的地形改正.地球物理學(xué)報(bào)，58(10):3815-3824,