[發明專利]一種陸地地震數據的全波形反演方法

有效

• 申請號: 201911327139.6
• 申請日: 2019-12-20
• 公開（公告）號: CN111045077B
• 公開（公告）日: 2022-03-18
• 發明（設計）人: 喬寶平;賀偉光;潘自強;曹成寅;黃昱丞;胡英才
• 申請（專利權）人: 核工業北京地質研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 核工業專利中心 11007
• 代理人: 閆兆梅

1.一種陸地地震數據的全波形反演方法，其特征在于，該方法具體包括以下步驟：

步驟一，建立陸地地震的光滑的初始化各向異性參數模型；

步驟二，從單炮地震記錄中獲取分頻帶的地震觀測記錄；

步驟三，根據參數模型和實際的地震觀測系統，獲取地震正演模擬合成數據集；

步驟四，計算合成地震數據與實際觀測地震數據的殘差δd_s,r；

步驟五，對于每一地震道，以體波第一個波包的中心作為高斯函數的中心點，計算高斯加權函數W_s,r(t)；其中，t代表時間變量；

步驟六，基于上述步驟四求取的數據殘差δd_s,r，計算基于KR-OT目標函數的優化函數其中，x_r代表檢波器的空間位置坐標；

步驟七，利用上述步驟五求取的高斯加權函數W_s,r(t)，計算優化后的KR-OT伴隨場震源D_d；

步驟八，設定基于優化的KR-OT伴隨場震源的目標函數計算目標函數參數模型更新的方向和大小，并利用目標函數參數模型代替步驟三中使用的參數模型；所述的步驟八中的基于優化的KR-OT伴隨場震源的目標函數其中，max代表求取其中最大值，Δx_r代表檢波器間距，代表地震觀測數據集，代表地震合成數據集，下標s代表震源角標，下標r代表檢波器角標；

步驟九，重復上述步驟三至步驟八，直至最終的數據殘差滿足精度要求，此時的參數模型即為全波形反演在當前頻帶所獲得的陸地地震數據的各項異性參數模型。

2.根據權利要求1所述的一種陸地地震數據的全波形反演方法，其特征在于：所述的步驟二的具體步驟如下：從單炮地震記錄中提取每炮記錄的震源子波，并從低頻到高頻，對觀測數據進行帶通濾波，獲取分頻帶的地震觀測記錄。

3.根據權利要求2所述的一種陸地地震數據的全波形反演方法，其特征在于：所述的步驟三的具體步驟如下：從最低頻帶開始，采用交錯網格，開展復雜地質情形下的地震波正演模擬計算，空間精度8階，時間精度2階，獲取地震正演模擬合成數據集。

4.根據權利要求3所述的一種陸地地震數據的全波形反演方法，其特征在于：所述的步驟四中的合成地震數據與實際觀測地震數據的殘差

5.根據權利要求4所述的一種陸地地震數據的全波形反演方法，其特征在于：所述的步驟五中的高斯加權函數其中，代表體波第一個波包的中心位置，σ代表方差。

6.根據權利要求5所述的一種陸地地震數據的全波形反演方法，其特征在于：所述的步驟七中的伴隨場震源其中，為Blip1，1-Lipschitz函數。