一種二維網(wǎng)格數(shù)據(jù)的補償方向濾波方法
【專利摘要】本發(fā)明涉及地球物理勘探的重力、磁力勘探數(shù)據(jù)處理技術(shù)，是一種二維網(wǎng)格數(shù)據(jù)的補償方向濾波方法，將原始數(shù)據(jù)進(jìn)行網(wǎng)格化處理，在窗口內(nèi)劃分方位角和每個方位對應(yīng)的窗口內(nèi)的節(jié)點分布，設(shè)置每個方位上背景場濾波時節(jié)點系數(shù)和剩余場補償濾波求取時的節(jié)點系數(shù)，計算窗口內(nèi)濾波方位的背景場值，循環(huán)迭代完成補償方向濾波。本發(fā)明濾波的方向性強(qiáng)，滿足了重磁處理提取方向異常信息的需要，同時也能滿足壓制某方向異常信息的需要。本發(fā)明推動了重磁勘探數(shù)據(jù)處理解釋技術(shù)的進(jìn)步，同時也為其它學(xué)科的二維網(wǎng)格數(shù)據(jù)的方向濾波處理提供了新方法。
【專利說明】一種二維網(wǎng)格數(shù)據(jù)的補償方向濾波方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地球物理勘探的重力、磁力勘探數(shù)據(jù)處理領(lǐng)域，屬于數(shù)字濾波技術(shù)，是一種二維網(wǎng)格數(shù)據(jù)的補償方向濾波方法。
【背景技術(shù)】
[0002]在磁力(包括航磁)資料采集中，由于測線間的一致性差異等原因引起的沿線的線性異常是具有方向的，這些方向性異常是地質(zhì)研究的干擾(即數(shù)據(jù)噪音)，因而，在重磁勘探資料處理解釋中需要提取具有方向性的異常信息，以滿足解釋斷層、構(gòu)造的需要，需要在資料處理解釋前進(jìn)行處理，即利用方向濾波方法壓制重磁數(shù)據(jù)中方向性的重磁噪音數(shù)據(jù)，使得濾波處理后的重磁數(shù)據(jù)更好地反映地質(zhì)信息是需要解決的問題。
[0003]目前地球物理勘探的重力、磁力勘探數(shù)據(jù)處理中常用的頻率域橢圓方向濾波是在頻率域內(nèi)實現(xiàn)，它的作用是提取出沿選擇方向的方向性異常。這種方法常見的步驟是:
1.將重力(或磁力)數(shù)據(jù)網(wǎng)格化，2.對重力網(wǎng)格數(shù)據(jù)進(jìn)行傅立葉變換，3.選擇濾波方向，選取橢圓方向濾波器，4.對傅立葉變換后的數(shù)據(jù)進(jìn)行橢圓方向濾波，獲得濾波后的頻譜數(shù)據(jù)，
5.對濾波后的頻譜數(shù)據(jù)進(jìn)行反傅立葉變換，即獲得橢圓方向濾波后的重力(或磁力)異常網(wǎng)格數(shù)據(jù)。這種方法不具有迭代性和局部場補償作用，而且方向性差。
[0004]地質(zhì)出版社，2005.6，P194_204公開了一種《重力場與重力勘探》方法，利用空間域滑動窗口加權(quán)平均法建立一種圓滑濾波器，對數(shù)據(jù)具有圓滑濾波作用，可獲得區(qū)域場，但是這種方法的方向濾波方法性能差，方向性異常往往被圓滑，無法提取重磁勘探資料處理解釋中具有方向性的異常信息，不能利用方向濾波方法壓制重磁數(shù)據(jù)中方向性的重磁噪音數(shù)據(jù)。

【發(fā)明內(nèi)容】

[0005]本發(fā)明目的在于為重磁數(shù)據(jù)處理提供一種具有較好方向濾波功能，可提取方向性重磁異常數(shù)據(jù)或壓制具有方向性的噪音數(shù)據(jù)的二維網(wǎng)格數(shù)據(jù)的補償方向濾波方法。
[0006]本發(fā)明通過以下步驟實現(xiàn):
[0007]I)進(jìn)行野外重力或磁力勘探，采集得到原始數(shù)據(jù)并進(jìn)行網(wǎng)格化處理；
[0008]2)運算窗口采用7行7列窗口，迭代次數(shù)取I?10 ;
[0009]步驟2)所述的運算窗口大小為7行7列。
[0010]3)在窗口內(nèi)劃分方位角和每個方位對應(yīng)的窗口內(nèi)的節(jié)點分布；
[0011]步驟3)所述的窗口內(nèi)劃分方位角是劃分為南和北方位、東和西方位、北東和西南方位、北西和東南方位、北北東和南南西方位、北東東和南西西方位、北北西和南南東方位、北西西和南東東方位16個。
[0012]步驟3)所述的窗口內(nèi)每對方位的節(jié)點分布相同，按以下分布:
[0013]南和北這對方位的節(jié)點分布是一組，窗口節(jié)點的行序為從下到上，窗口節(jié)點的列的序為從左到右，坐標(biāo)組行序，列序為:南和北方位的窗口內(nèi)節(jié)點分布為(I，4)、(2,4)、(3,4)、(4，4)、(5，4)、(6，4)、(7，4)；
[0014]東和西方位的節(jié)點分布是一組，東和西方位的窗口內(nèi)節(jié)點分布為(4，I)、(4,2),(4，3)、(4，4)、(4，5)、(4，6)、(4，7)；
[0015]北東和南西方位的節(jié)點分布是一組，北東和南西方位的窗口內(nèi)節(jié)點分布為(1，I)、(2，2)、(3，3)、(4，4)、(5，5)、(6，6)、(7，7)；
[0016]北西和南東方位的節(jié)點分布是一組，北西和南東方位的窗口內(nèi)節(jié)點分布為(7，I )、(6，2)、(5，3)、(4，4)、(3，5)、(2，6)、(1，7)；
[0017]北北東和南南西方位的節(jié)點分布是一組(圖1的5)，北北東和南南西方位的窗口內(nèi)節(jié)點分布為(1，3)、(2，3)、(3，3)、(3，4)、(4，4)、(5，4)、(5，5)、(6，5)、(7，5)；
[0018]北東東和南西西方位的節(jié)點分布是一組(圖1的6)，北東東和南西西方位的窗口內(nèi)節(jié)點分布為(3，I)、(3，2)、(3，3)、(4，3)、(4，4)、(4，5)、(5，5)、(5，6)、(5，7)；
[0019]北北西和南南東方位的節(jié)點分布是一組，北北西和南南東方位的窗口內(nèi)節(jié)點分布為(1，5)、(2，5)、(3，5)、(3，4)、(4，4)、(5，4)、(5，3)、(6，3)、(7，3)；
[0020]北西西和南東東方位的節(jié)點分布是一組，北西西和南東東方位的窗口內(nèi)節(jié)點分布為(5，I)、(5，2)、(5，3)、(4，3)、(4，4)、(4，5)、(3，5)、(3，6)、(3，7)。
[0021]4)設(shè)置每個方位上背景場濾波時節(jié)點系數(shù)和剩余場補償濾波求取時的節(jié)點系數(shù)；
[0022]步驟4)所述的節(jié)點系數(shù)為:
[0023]東西向濾波時，背景場濾波時節(jié)點系數(shù)和剩余場補償濾波求取時的節(jié)點系數(shù)為:(4，1)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)節(jié)點的背景場濾波時節(jié)點系數(shù)自左向右分別為 3、5、7、12、7、5、3 ;剩余場償濾波求取時的(1，4)、(2，4)、(3，4)、(4，4)、(5，4)、(6，4)、(7，4)節(jié)點系數(shù)分別為1、1、1、1、1、1、1 ；
[0024]南北向濾波時，背景場濾波(1，4)、(2，4)、(3，4)、(4，4)、(5，4)、(6，4)、(7，4)節(jié)點系數(shù)分別為 3、5、7、12、7、5、3，剩余場償濾波(4，1)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)節(jié)點的系數(shù)分別為1、1、1、1、1、1、1 ；
[0025]其它各個方位的背景場濾波時節(jié)點系數(shù)和剩余場補償濾波求取時的節(jié)點系數(shù)設(shè)置采用東西向濾波和南北向濾波設(shè)置方法類推。
[0026]5)按照以下公式計算窗口內(nèi)濾波方位的背景場值gw:
[0027]gbJ = (I/ Σ (k (i, j))) X (Σ (k (i, j) Xg(i, j)))
[0028]式中:
[0029]k (i, j)為節(jié)點(i，j)的系數(shù)值，g(l, j)為節(jié)點(i，j)的重力或磁力值；
[0030]步驟5)所述的濾波方位的背景場值gw，當(dāng)東西向濾波時，先計算計算點(4，4)的背景場值gw (4，4):
[0031]gbJ (4，4)= (1/42)X (3Xg(4, I)+5Xg(4，2)+7Xg(4，3)+12Xg(4，4)+7Xg(4，5) +5Xg(4, 6)+3Xg(4, 7))；
[0032]式中g(shù)(4，I)、g(4，2)、g(4，3)、g(4，4)、g(4，5)、g(4，6)、g(4，7)分別表示計算窗 P內(nèi)節(jié)點(4，I)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)的重力值或者是磁力值；
[0033]依次移動運算窗口，計算工區(qū)平面網(wǎng)格的每一個計算點的背景場值；
[0034]其它各方位濾波，替換所選坐標(biāo)節(jié)點和對應(yīng)節(jié)點的系數(shù)、按各方位的系數(shù)加權(quán)平均計算。
[0035]6)按照以下公式計算窗口內(nèi)濾波方位的剩余場補償濾波值gbc;h:
[0036]gbch= (1/ Σ (kbch (i, j))) X (Σ (kbch (i, j) Xgbch(i, j)))
[0037]式中:
[0038]kbch (i, j)為節(jié)點(i，j)的系數(shù)值kbc;h，g(i, j)為節(jié)點(i，j)的重力或磁力值；
[0039]步驟6)所述的剩余場補償濾波值gbc;h,東西向濾波時:
[0040]gbch (4,4) =gbJ (4,4)-( (1/7) X (gbJ(1, 4) +gbJ (2，4) +gbJ (3，4) +gbJ (4，4) +gbJ (5，4) +gbJ(6,4)+gbJ(7,4)))；
[0041]然后,計算計算點的補償方向濾波輸出值gf(4,4),
[0042]gf (4,4) =g (4, 4) -gbJ (4,4) +gbch (4,4)；
[0043]再計算該點的剩余異常值gs(4，4)，
[0044]gs (4，4) =g (4，4) -gf (4,4)；
[0045]依次移動運算窗口，計算每一個計算點的補償方向濾波值gfl(i，j)、剩余異常值
Ss1 (i，j)。
[0046]7)循環(huán)迭代直到完成迭代運算次數(shù)，完成補償方向濾波。
[0047]所述的循環(huán)迭代是將第一次計算得到的各個節(jié)點的補償方向濾波值gf計為gfl (i，j)，將第一次計算得到的各個節(jié)點的剩余異常值gs計為gs1 (i，j)，對gs1 (i，j)進(jìn)行下一次的補償方向濾波運算，即把gs1(i，j)視為步驟5)中的原始數(shù)據(jù)g(i，j)，進(jìn)行迭代，計算步驟5)、步驟6)，獲得新的g’ f (i，j)、g’ s(i，j)，則第二次迭代的結(jié)果為:
[0048]gf2(i, j)=gfl(i，j)+g，f(i，j).[0049]gs2(i, j)=gsl(i，j)_g，f(i，j).[0050]式中(i，j)表示計算點的節(jié)點在工區(qū)坐標(biāo)，g’ f、g’ s分別表示本次迭代中計算得到的計算點的補償方向濾波、剩余異常值。
[0051]本發(fā)明對實測數(shù)據(jù)處理，獲得了多次迭代后的補償方向濾波異常和補償方向濾波剩余異常，濾波的方向性強(qiáng)，滿足了重磁處理提取方向異常信息的需要，同時也能滿足壓制某方向異常信息的需要。本發(fā)明推動了重磁勘探數(shù)據(jù)處理解釋技術(shù)的進(jìn)步，同時也為其它學(xué)科的二維網(wǎng)格數(shù)據(jù)的方向濾波處理提供了一種新方法。
【專利附圖】

【附圖說明】
[0052]圖1是本發(fā)明二維網(wǎng)格的窗口選取和16方位節(jié)點分布圖；
[0053]圖2是本發(fā)明背景場濾波節(jié)點系數(shù)和補償場濾波節(jié)點系數(shù)分布圖；
[0054]圖3是原始磁異常平面圖；
[0055]圖4是本發(fā)明南北向濾波磁異常平面圖；
[0056]圖5是本發(fā)明南北向濾波剩余磁異常平面圖。
【具體實施方式】
[0057]以下結(jié)合實例詳細(xì)說明本發(fā)明。
[0058]本實例通過以下步驟實現(xiàn):
[0059]1)進(jìn)行野外航空磁力勘探，采集得到航空磁力原始數(shù)據(jù)并進(jìn)行網(wǎng)格化處理，見圖3 ；
[0060]2)運算窗口采用7行7列窗口(見圖1)，迭代次數(shù)取2;
[0061]步驟2)所述的運算窗口大小為7行7列，見圖1。
[0062]3)在窗口內(nèi)劃分方位角和每個方位對應(yīng)的窗口內(nèi)的節(jié)點分布，見圖1 ;
[0063]步驟3)所述的窗口內(nèi)劃分方位角是劃分為南和北方位、東和西方位、北東和西南方位、北西和東南方位、北北東和南南西方位、北東東和南西西方位、北北西和南南東方位、北西西和南東東方位16個，見圖1。
[0064]步驟3)所述的窗口內(nèi)每對方位的節(jié)點分布相同，按以下分布:
[0065]南和北這對方位的節(jié)點分布是一組，見圖1的1，窗口節(jié)點的行序為從下到上，窗口節(jié)點的列的序為從左到右，坐標(biāo)組行序，列序為:南和北方位的窗口內(nèi)節(jié)點分布為(I，
4)、(2，4)、(3，4)、(4，4)、(5，4)、(6，4)、(7，4)；
[0066]東和西方位的 節(jié)點分布是一組，見圖1的2，東和西方位的窗口內(nèi)節(jié)點分布為(4，I)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)；
[0067]北東和南西方位的節(jié)點分布是一組，見圖1的3，北東和南西方位的窗口內(nèi)節(jié)點分布為(1，I)、(2，2)、(3，3)、(4，4)、(5，5)、(6，6)、(7，7)；
[0068]北西和南東方位的 節(jié)點分布是一組，見圖1的4，北西和南東方位的窗口內(nèi)節(jié)點分布為(7，I)、(6，2)、(5，3)、(4，4)、(3，5)、(2，6)、(1，7)；
[0069]北北東和南南西方位的節(jié)點分布是一組，見圖1的5，北北東和南南西方位的窗口內(nèi)節(jié)點分布為(1，3)、(2，3)、(3，3)、(3，4)、(4，4)、(5，4)、(5，5)、(6，5)、(7，5)；
[0070]北東東和南西西方位的節(jié)點分布是一組，見圖1的6，北東東和南西西方位的窗口內(nèi)節(jié)點分布為(3，I)、(3，2)、(3，3)、(4，3)、(4，4)、(4，5)、(5，5)、(5，6)、(5，7)；
[0071]北北西和南南東方位的節(jié)點分布是一組，見圖1的7，北北西和南南東方位的窗口內(nèi)節(jié)點分布為(1，5)、(2，5)、(3，5)、(3，4)、(4，4)、(5，4)、(5，3)、(6，3)、(7，3)；
[0072]北西西和南東東方位的節(jié)點分布是一組，見圖1的8，北西西和南東東方位的窗口內(nèi)節(jié)點分布為(5，I)、(5，2)、(5，3)、(4，3)、(4，4)、(4，5)、(3，5)、(3，6)、(3，7)。
[0073]4)設(shè)置南和北方位上背景場濾波時節(jié)點系數(shù)和剩余場補償濾波求取時的節(jié)點系數(shù)；
[0074]圖2是背景場濾波節(jié)點系數(shù)和補償場濾波節(jié)點系數(shù)分布圖。圖中I為東-西向，2為南-北向，3為北西-東南向,4為東北-西南向，5為北西西-南南東向,6為北北西-南南東向，7為北東東-西西南向，8為北北東-南南西向。
[0075]步驟4)所述的節(jié)點系數(shù)為:
[0076]南北向濾波時，背景場濾波(1，4)、(2，4)、(3，4)、(4，4)、(5，4)、(6，4)、(7，4)節(jié)點系數(shù)分別為 3、5、7、12、7、5、3，剩余場償濾波(4，1)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，
7)節(jié)點的系數(shù)分別為1、1、1、1、1、1、1，節(jié)點系數(shù)分布見圖2的I ;
[0077]5)按照以下公式計算窗口內(nèi)濾波方位的背景場值gw:
[0078]gbJ = (I/ Σ (k (i, j))) X (Σ (k (i, j) Xg(i, j)))
[0079]式中:k (i，j)為節(jié)點(i，j)的系數(shù)值，g(i，j)為節(jié)點(i，j)的重力或磁力值；
[0080]步驟5)所述的濾波方位的背景場值gw，南北方向濾波時，先計算計算點(4，4)的背景場值gw (4，4):[0081]
【權(quán)利要求】
1.一種二維網(wǎng)格數(shù)據(jù)的補償方向濾波方法，特點是通過以下步驟實現(xiàn): 1)進(jìn)行野外重力或磁力勘探，采集得到原始數(shù)據(jù)并進(jìn)行網(wǎng)格化處理； 2)運算窗口采用7行7列窗口，迭代次數(shù)取I~10; 3)在窗口內(nèi)劃分方位角和每個方位對應(yīng)的窗口內(nèi)的節(jié)點分布； 4)設(shè)置每個方位上背景場濾波時節(jié)點系數(shù)和剩余場補償濾波求取時的節(jié)點系數(shù)； 5)按照以下公式計算窗口內(nèi)濾波方位的背景場值gw: gbJ =(1/ Σ (k (i, j))) X (Σ (k (i, j) Xg(i, j))) 式中: k (i, j)為節(jié)點(i，j)的系數(shù)值，g(i, j)為節(jié)點(i，j)的重力或磁力值； 6)按照以下公式計算窗口內(nèi)濾波方位的剩余場補償濾波值gbc;h:
gbch= (I/ Σ (kbch (i，j)))X (Σ (kbch (i，j) Xgbch(i，j))) 式中: kbch (i, j)為節(jié)點(i，j)的系數(shù)值kbc;h，g(i, j)為節(jié)點(i，j)的重力或磁力值； 7 )循環(huán)迭代直到完成迭代運算次數(shù)，完成補償方向濾波。
2.根據(jù)權(quán)利要求1所述的方法，特點是步驟2)所述的運算窗口大小為7行7列。`
3.根據(jù)權(quán)利要求1所述的方法，特點是步驟3)所述的窗口內(nèi)劃分方位角是劃分為南和北方位、東和西方位、北東和西南方位、北西和東南方位、北北東和南南西方位、北東東和南西西方位、北北西和南南東方位、北西西和南東東方位16個。
4.根據(jù)權(quán)利要求1所述的方法，特點是步驟3)所述的窗口內(nèi)每對方位的節(jié)點分布相同，按以下分布: 南和北這對方位的節(jié)點分布是一組，窗口節(jié)點的行序為從下到上，窗口節(jié)點的列的序為從左到右，坐標(biāo)組行序，列序為:南和北方位的窗口內(nèi)節(jié)點分布為(1，4)、(2,4)、(3,4)、(4，4)、(5，4)、(6，4)、(7，4);東和西方位的節(jié)點分布是一組，東和西方位的窗口內(nèi)節(jié)點分布為(4，I)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)； 北東和南西方位的節(jié)點分布是一組，北東和南西方位的窗口內(nèi)節(jié)點分布為(1，I)、(2，2)、(3，3)、(4，4)、(5，5)、(6，6)、(7，7)； 北西和南東方位的節(jié)點分布是一組，北西和南東方位的窗口內(nèi)節(jié)點分布為(7，I)、(6，2)、(5，3)、(4，4)、(3，5)、(2，6)、(1，7)； 北北東和南南西方位的節(jié)點分布是一組(圖1的5)，北北東和南南西方位的窗口內(nèi)節(jié)點分布為(1，3)、(2，3)、(3，3)、(3，4)、(4，4)、(5，4)、(5，5)、(6，5)、(7，5)； 北東東和南西西方位的節(jié)點分布是一組(圖1的6)，北東東和南西西方位的窗口內(nèi)節(jié)點分布為(3，I)、(3，2)、(3，3)、(4，3)、(4，4)、(4，5)、(5，5)、(5，6)、(5，7)； 北北西和南南東方位的節(jié)點分布是一組，北北西和南南東方位的窗口內(nèi)節(jié)點分布為(1，5)、(2，5)、(3，5)、(3，4)、(4，4)、(5，4)、(5，3)、(6，3)、(7，3)； 北西西和南東東方位的節(jié)點分布是一組，北西西和南東東方位的窗口內(nèi)節(jié)點分布為(5，I)、(5，2)、(5，3)、(4，3)、(4，4)、(4，5)、(3，5)、(3，6)、(3，7)。
5.根據(jù)權(quán)利要求1所述的方法，特點是步驟4)所述的節(jié)點系數(shù)為: 東西向濾波時，背景場濾波時節(jié)點系數(shù)和剩余場補償濾波求取時的節(jié)點系數(shù)為:(4，I)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)節(jié)點的背景場濾波時節(jié)點系數(shù)自左向右分別為 3、5、7、12、7、5、3 ;剩余場償濾波求取時的(1，4)、(2，4)、(3，4)、(4，4)、(5，4)、(6，4)、(7，4)節(jié)點系數(shù)分別為1、1、1、1、1、1、1 ； 南北向濾波時，背景場濾波(1，4)、(2，4)、(3，4)、(4，4)、(5，4)、(6，4)、(7，4)節(jié)點系數(shù)分別為 3、5、7、12、7、5、3，剩余場償濾波(4，1)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)節(jié)點的系數(shù)分別為1、1、1、1、1、1、I ； 其它各個方位的背景場濾波時節(jié)點系數(shù)和剩余場補償濾波求取時的節(jié)點系數(shù)設(shè)置采用東西向濾波和南北向濾波設(shè)置方法類推。
6.根據(jù)權(quán)利要求1所述的方法，特點是步驟5)所述的濾波方位的背景場值gw，當(dāng)東西向濾波時，先計算計算點(4，4)的背景場值gw (4，4):
gbJ(4，4)= (1/42)X (3Xg(4, I)+5Xg(4, 2)+7Xg(4, 3)+12Xg(4, 4) +7Xg(4, 5)+5Xg(4，6)+3 Xg (4, 7))； 式中 g(4, I)、g(4, 2)、g(4, 3)、g(4, 4)、g(4, 5)、g(4, 6)、g(4, 7)分別表示計算窗口內(nèi)節(jié)點(4，I)、(4，2)、(4，3)、(4，4)、(4，5)、(4，6)、(4，7)的重力值或者是磁力值； 依次移動運算窗口，計算工區(qū)平面網(wǎng)格的每一個計算點的背景場值； 其它各方位濾波，替換所選坐標(biāo)節(jié)點和對應(yīng)節(jié)點的系數(shù)、按各方位的系數(shù)加權(quán)平均計笪
7.根據(jù)權(quán)利要求1所述的方法，特點是步驟6)所述的剩余場補償濾波值gbc;h，東西向濾波時:`
gbch (4，4) =gbJ (4,4)-( (1/7) X (gbj (I, 4) +gbJ (2，4) +gbJ (3，4) +gbJ (4，4) +gbJ (5，4) +gbJ (6，4)+gbJ(7,4)))； 然后，計算計算點的補償方向濾波輸出值gf (4，4)， gf (4，4) =g (4，4) -gbJ (4，4) +gbch (4,4)； 再計算該點的剩余異常值gs(4，4)， gs (4,4) =g (4, 4) -gf (4,4)； 依次移動運算窗口，計算每一個計算點的補償方向濾波值gfl(i，j)、剩余異常值gsl(i，j)。
8.根據(jù)權(quán)利要求1所述的方法，特點是步驟7)所述的循環(huán)迭代是將第一次計算得到的各個節(jié)點的補償方向濾波值gf計為gfl(i，j)，將第一次計算得到的各個節(jié)點的剩余異常值gs計為gsi(i，j)，對gsl(i，j)進(jìn)行下一次的補償方向濾波運算，即把gsl(i，j)視為步驟5)中的原始數(shù)據(jù)g(i，j)，進(jìn)行迭代，計算步驟5)、步驟6)，獲得新的g’ f(i，j)、g’ s(i，j)，則第二次迭代的結(jié)果為:
gf2(i，j)=gn(i，j)+g，f(i，j);
gs2(i，j)=gsi(i，j)-g，f(i，j); 式中(i，j)表示計算點的節(jié)點在工區(qū)坐標(biāo)，g’ f、g’s分別表示本次迭代中計算得到的計算點的補償方向濾波、剩余異常值。
【文檔編號】G01V7/02GK103513288SQ201210211223
【公開日】2014年1月15日 申請日期:2012年6月21日 優(yōu)先權(quán)日:2012年6月21日
【發(fā)明者】劉云祥 申請人:中國石油天然氣集團(tuán)公司, 中國石油集團(tuán)東方地球物理勘探有限責(zé)任公司