[發明專利]一種二維網格數據的補償方向濾波方法

有效

• 申請號: 201210211223.3
• 申請日: 2012-06-21
• 公開（公告）號: CN103513288A
• 公開（公告）日: 2014-01-15
• 發明（設計）人: 劉云祥
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V7/02
• 分類號: G01V7/02;G01V3/38
• 代理公司: 北京市中實友知識產權代理有限責任公司 11013
• 代理人: 劉天語;武尚修

技術領域

本發明涉及地球物理勘探的重力、磁力勘探數據處理領域，屬于數字濾波技術，是一種二維網格數據的補償方向濾波方法。

背景技術

在磁力（包括航磁）資料采集中，由于測線間的一致性差異等原因引起的沿線的線性異常是具有方向的，這些方向性異常是地質研究的干擾（即數據噪音），因而，在重磁勘探資料處理解釋中需要提取具有方向性的異常信息，以滿足解釋斷層、構造的需要，需要在資料處理解釋前進行處理，即利用方向濾波方法壓制重磁數據中方向性的重磁噪音數據，使得濾波處理后的重磁數據更好地反映地質信息是需要解決的問題。

目前地球物理勘探的重力、磁力勘探數據處理中常用的頻率域橢圓方向濾波是在頻率域內實現，它的作用是提取出沿選擇方向的方向性異常。這種方法常見的步驟是：1.將重力（或磁力）數據網格化，2.對重力網格數據進行傅立葉變換，3.選擇濾波方向，選取橢圓方向濾波器，4.對傅立葉變換后的數據進行橢圓方向濾波，獲得濾波后的頻譜數據，5.對濾波后的頻譜數據進行反傅立葉變換，即獲得橢圓方向濾波后的重力（或磁力）異常網格數據。這種方法不具有迭代性和局部場補償作用，而且方向性差。

地質出版社，2005.6，P194-204公開了一種《重力場與重力勘探》方法，利用空間域滑動窗口加權平均法建立一種圓滑濾波器，對數據具有圓滑濾波作用，可獲得區域場，但是這種方法的方向濾波方法性能差，方向性異常往往被圓滑，無法提取重磁勘探資料處理解釋中具有方向性的異常信息，不能利用方向濾波方法壓制重磁數據中方向性的重磁噪音數據。

發明內容

本發明目的在于為重磁數據處理提供一種具有較好方向濾波功能，可提取方向性重磁異常數據或壓制具有方向性的噪音數據的二維網格數據的補償方向濾波方法。

本發明通過以下步驟實現：

1）進行野外重力或磁力勘探，采集得到原始數據并進行網格化處理；

2）運算窗口采用7行7列窗口，迭代次數取1～10；

步驟2）所述的運算窗口大小為7行7列。

3）在窗口內劃分方位角和每個方位對應的窗口內的節點分布；

步驟3）所述的窗口內劃分方位角是劃分為南和北方位、東和西方位、北東和西南方位、北西和東南方位、北北東和南南西方位、北東東和南西西方位、北北西和南南東方位、北西西和南東東方位16個。

步驟3）所述的窗口內每對方位的節點分布相同，按以下分布：

南和北這對方位的節點分布是一組，窗口節點的行序為從下到上，窗口節點的列的序為從左到右，坐標組行序，列序為：南和北方位的窗口內節點分布為（1，4）、（2，4）、（3，4）、（4，4）、（5，4）、（6，4）、（7，4）；

東和西方位的節點分布是一組，東和西方位的窗口內節點分布為（4，1）、（4，2）、（4，3）、（4，4）、（4，5）、（4，6）、（4，7）；

北東和南西方位的節點分布是一組，北東和南西方位的窗口內節點分布為（1，1）、（2，2）、（3，3）、（4，4）、（5，5）、（6，6）、（7，7）；

北西和南東方位的節點分布是一組，北西和南東方位的窗口內節點分布為（7，1）、（6，2）、（5，3）、（4，4）、（3，5）、（2，6）、（1，7）；

北北東和南南西方位的節點分布是一組(圖1的5)，北北東和南南西方位的窗口內節點分布為（1，3）、（2，3）、（3，3）、（3，4）、（4，4）、（5，4）、（5，5）、（6，5）、（7，5）；

北東東和南西西方位的節點分布是一組(圖1的6)，北東東和南西西方位的窗口內節點分布為（3，1）、（3，2）、（3，3）、（4，3）、（4，4）、（4，5）、（5，5）、（5，6）、（5，7）；

北北西和南南東方位的節點分布是一組，北北西和南南東方位的窗口內節點分布為（1，5）、（2，5）、（3，5）、（3，4）、（4，4）、（5，4）、（5，3）、（6，3）、（7，3）；

北西西和南東東方位的節點分布是一組，北西西和南東東方位的窗口內節點分布為（5，1）、（5，2）、（5，3）、（4，3）、（4，4）、（4，5）、（3，5）、（3，6）、（3，7）。

4）設置每個方位上背景場濾波時節點系數和剩余場補償濾波求取時的節點系數；

步驟4）所述的節點系數為：