[發明專利]地震數據速度模型建立方法和裝置

在審

• 申請號: 201510547221.5
• 申請日: 2015-08-31
• 公開（公告）號: CN105204063A
• 公開（公告）日: 2015-12-30
• 發明（設計）人: 章威
• 申請（專利權）人: 中國石油天然氣集團公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 黨曉林

1.一種地震數據速度模型建立方法，其特征在于，包括：

確定地震數據速度模型的目標函數；

基于層位約束計算得到所述目標函數的梯度；

根據計算得到的梯度對所述目標函數進行調整，以得到調整后的速度模型。

2.如權利要求1所述的方法，其特征在于，所述基于層位約束計算得到所述速度模型的梯度，包括：

計算所述目標函數的曲率；

將所述層位轉化為空間中的掩碼，其中，所述層位表示N維空間中的一組N-1維的流形；

將所述目標函數的曲率與所述掩碼相乘得到掩碼曲率，；

利用單程波動方程，將震源子波向下傳播，得到下行震源波場；

將所述下行震源波場與所述掩碼相乘得到第一相乘結果，將所述第一相乘結果用單程波動方程向上傳播，得到模擬數據；

將所述下行震源波場與所述掩碼曲率相乘得到第二相乘結果，將所述第二相乘結果用所述單程波動方程向上傳播，得到第一上行震源波場；

利用所述單程波動方程，將所述模擬數據向下傳播，得到下行數據波場；

將所述下行數據波場與所述掩碼曲率相乘得到第三相乘結果，將所述第三相乘結果用所述單程波動方程向上傳播，得到第二上行震源波場；

通過所述下行震源波場、所述第二上行震源波場、所述下行數據波場、和所述第一上行震源波場計算得到所述目標函數的梯度。

3.如權利要求2所述的方法，其特征在于，通過所述下行震源波場、所述第二上行震源波場、所述下行數據波場、和所述第一上行震源波場計算得到所述目標函數的梯度，包括：

按照以下公式計算所述目標函數的梯度：

G＝Re(U_sd*U_ru+U_rd*U_su)

其中，G表示所述目標函數的梯度，Re表示求復數的實部，U_sd、U_ru、U_rd、U_su都是復數，其中，*表示復數乘法。

4.如權利要求2所述的方法，其特征在于，將所述層位轉化為空間中的掩碼，包括：

對空間中的每個點，判斷該點是否在層位上，如果在，則將該點的掩碼置為1，如果不在，則將該點的掩碼置為0。

5.如權利要求4所述的方法，其特征在于：

在將所述第一相乘結果用單程波動方程向上傳播的過程中，當遇到掩碼為1時，停止傳播，記錄當前波場值，并將所述下行震源波場與掩碼的相乘結果作為震源繼續傳播；

在將所述第二相乘結果用所述單程波動方程向上傳播的過程中，當遇到掩碼曲率不為0時，停止傳播當前波場值，用當前的掩碼曲率替代波場值繼續傳播；

在將所述模擬數據向下傳播的過程中，當遇到模擬數據不為0時，停止傳播當前的波場值，用當前的模擬數據替代波場值繼續傳播；

在將所述第三相乘結果用所述單程波動方程向上傳播的過程中，當遇到掩碼曲率不為0時，停止傳播當前的波場值，用下行數據波場與掩碼曲率相乘的結果替代波場值繼續傳播。

6.如權利要求1至5中任一項所述的方法，其特征在于，確定地震數據速度模型的目標函數，包括：

在地震數據的偏移圖像域計算所述目標函數；

或者，在地震數據的數據域計算所述目標函數。

7.一種地震數據速度模型建立裝置，其特征在于，包括：

目標函數確定單元，用于確定地震數據速度模型的目標函數；

梯度計算單元，用于基于層位約束計算得到所述目標函數的梯度；

模型調整單元，用于根據計算得到的梯度對所述目標函數進行調整，以得到調整后的速度模型。