[發明專利]一種構造增強的梯度預處理全波形反演方法及系統

在審

• 申請號: 202011058205.7
• 申請日: 2020-09-30
• 公開（公告）號: CN114428338A
• 公開（公告）日: 2022-05-03
• 發明（設計）人: 胡光輝;孫思宇;何兵紅;杜澤源
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 北京知舟專利事務所(普通合伙) 11550
• 代理人: 郭韞

1.一種構造增強的梯度預處理全波形反演方法，其特征在于：所述方法根據噪音分布特點統計噪音分布，同時記錄噪音的位置，并對噪音的梯度進行重新賦值，實現對噪音的去除。

2.根據權利要求1所述的構造增強的梯度預處理全波形反演方法，其特征在于：所述方法包括：

(1)采集炮集記錄，獲得梯度G，并計算獲得梯度平均值；

(2)計算梯度G中每一個元素與梯度平均值的平方差，記為A(i，j)；

(3)將A(i,j)中的各個元素重新排序并記錄，得到一維數組A’(K)；

(4)設置要去除的噪音個數百分比p，計算異常點的個數c，并標記各個異常點；

(5)依次對梯度G中的每個異常點的梯度進行重新賦值，得到每個異常點的新的梯度。

3.根據權利要求2所述的構造增強的梯度預處理全波形反演方法，其特征在于：所述步驟(1)的操作包括：

利用下式計算獲得梯度平均值S：

其中，S為梯度平均值，M，N為梯度在橫向和豎向的采樣點個數，i，j為計數符號，表示梯度G中的某一元素。

4.根據權利要求3所述的構造增強的梯度預處理全波形反演方法，其特征在于：所述步驟(2)的操作包括：

利用下式計算獲得梯度G中每一個元素與梯度平均值的平方差：

A(i,j)＝(G(i,j)-S)²,i＝1,...,M；j＝1,...,N (2)。

5.根據權利要求4所述的構造增強的梯度預處理全波形反演方法，其特征在于：所述步驟(3)的操作包括：

利用下式對數組A(i，j)進行重新排序，得到一維數組A’(K)：

A′(k)＝arr[A(i,j)] (3)

其中，arr表示將A(i，j)按照從小到大的順序進行排列；

K為一維數組的計數符號，取值區間為[1M*N]。

6.根據權利要求5所述的構造增強的梯度預處理全波形反演方法，其特征在于：所述步驟(4)的操作包括：

利用下式計算異常點的個數c：

c＝INT(M*N*p/2) (4)

其中，INT表示取整數；

將A’(K)中的前c個、后c個元素均標記為異常點，并記錄各個異常點在梯度G中的坐標(i,j)。

7.根據權利要求6所述的構造增強的梯度預處理全波形反演方法，其特征在于：所述步驟(4)中的p的取值為1.0％到5.0％之間。

8.根據權利要求7所述的構造增強的梯度預處理全波形反演方法，其特征在于：所述步驟(5)的操作包括：

找到該異常點周圍的點中所有沒有被標記為異常點的點，并對這些點的梯度取算術平均值，然后將該算術平均值賦值給該異常點，得到該異常點的新的梯度。

9.一種構造增強的梯度預處理全波形反演系統，其特征在于：所述系統包括：

數據采集單元，用于采集炮集記錄，獲得梯度G，并計算獲得梯度平均值；

平方差計算單元，與所述數據采集單元連接，用于計算梯度G中每一個元素與梯度平均值的平方差，記為A(i，j)；

排序單元，與所述平方差計算單元連接，用于將A(i,j)中的各個元素重新排序并記錄，得到一維數組A’(K)；

異常點標記單元，與所述排序單元連接，用于設置要去除的噪音個數百分比p，計算異常點的個數c，并標記各個異常點；

梯度賦值單元，與所述異常點標記單元連接，用于依次對梯度G中的每個異常點的梯度進行重新賦值，得到每個異常點的新的梯度。

10.一種計算機可讀存儲介質，其特征在于：所述計算機可讀存儲介質存儲有計算機可執行的至少一個程序，所述至少一個程序被所述計算機執行時使所述計算機執行權利要求1-8之任一項所述的構造增強的梯度預處理全波形反演方法中的步驟。