[發明專利]一種時移重力增量的去噪方法和裝置

有效

• 申請號: 201610517205.6
• 申請日: 2016-07-04
• 公開（公告）號: CN105974488B
• 公開（公告）日: 2018-02-02
• 發明（設計）人: 劉云祥;趙文舉;鄭勝賢
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V7/06
• 分類號: G01V7/06
• 代理公司: 北京三友知識產權代理有限公司11127
• 代理人: 李輝,徐煥

1.一種時移重力增量的去噪方法，其特征在于，包括：

獲取工區內各個測點的時移重力增量，其中，所述時移重力增量為測點在第一時間段內多次測量的重力平均值和在第二時間段內多次測量的重力平均值之間的差值；

獲取所述各個測點在第一時間段內多次測量的重力方差以及在第二時間段內多次測量的重力方差；

對所述各個測點中每個測點采用以下操作得到所述各個測點中每個測點去噪后的時移重力增量：選擇當前測點作為中心點，將以所述中心點為中心，以預設去噪區間作為邊界所限定的范圍內的所有測點在第一時間段內多次測量的重力方差、在第二時間段內多次測量的重力方差和時移重力增量進行加權處理，得到去噪后的當前測點的時移重力增量。

2.如權利要求1所述的方法，其特征在于，在得到所述各個測點中每個測點去噪后的時移重力增量之后，所述方法還包括：

將所述各個測點中每個測點去噪后的時移重力增量作為未去噪的時移重力增量，對所述各個測點中每個測點重復執行以下去噪操作直至達到預設的重復次數：選擇當前測點作為中心點，將以所述中心點為中心，以預設去噪區間作為邊界所限定的范圍內的所有測點在第一時間段內多次測量的重力方差、在第二時間段內多次測量的重力方差和未去噪的時移重力增量進行加權處理，得到去噪后的當前測點的時移重力增量。

3.如權利要求2所述的方法，其特征在于，所述預設的重復次數不小于10次。

4.如權利要求1所述的方法，其特征在于，按照以下公式計算得到去噪后的所述當前測點的時移重力增量：

當A<(D/2)和B<(D/2)時，M＝N；

當(D/2)<A<(3D/4)和(D/2)<B<(3D/4)時，M＝N×a+P×b；

當(D/2)<A<(3D/4)和B>(3D/4)、(D/2)<B<(3D/4)和A>(3D/4)滿足二者之一時，M＝N×c+P×d；

當A>(3D/4)且B>(3D/4)時，M＝N×e+P×f；

其中，A表示以所述預設去噪區間作為邊界所限定的范圍內的所有測點在第一時間段內多次測量的重力方差，B表示在第二時間段內多次測量的重力方差，D表示以所述預設去噪區間作為邊界所限定的范圍內的所有測點在第一時間段內多次測量的重力方差以及在第二時間段內多次測量的重力方差中的最大值，M表示去噪后的所述當前測點的時移重力增量，N表示所述時移重力增量，P表示以所述預設去噪區間作為邊界所限定的范圍內的所有測點的時移重力增量的平均值，a、b、c、d、e、f均表示在0％至100％之間的加權系數。

5.如權利要求1所述的方法，其特征在于，獲取工區內各個測點的時移重力增量，包括：對所述時移動重力增量進行網格化處理；

獲取所述各個測點在第一時間段內多次測量的重力方差以及在第二時間段內多次測量的重力方差，包括：對所述各個測點在第一時間段內多次測量的重力方差以及在第二時間段內多次測量的重力方差進行網格化處理。

6.如權利要求5所述的方法，其特征在于，所述網格化處理方法包括：克里金插值法。

7.如權利要求1所述的方法，其特征在于，所述預設去噪區間包括：n×n的矩陣，其中，所述n表示所述預設去噪區間中的測點個數。

8.一種時移重力增量的去噪裝置，其特征在于，包括：

平均值獲取模塊，用于獲取工區內各個測點的時移重力增量，其中，所述時移重力增量為測點在第一時間段內多次測量的重力平均值和在第二時間段內多次測量的重力平均值之間的差值；

方差獲取模塊，用于獲取所述各個測點在第一時間段內多次測量的重力方差以及在第二時間段內多次測量的重力方差；

加權處理模塊，用于對所述各個測點中每個測點采用以下操作得到所述各個測點中每個測點去噪后的時移重力增量：選擇當前測點作為中心點，將以所述中心點為中心，以預設去噪區間作為邊界所限定的范圍內的所有測點在第一時間段內多次測量的重力方差、在第二時間段內多次測量的重力方差和時移重力增量進行加權處理，得到去噪后的所述當前測點的時移重力增量。