[發明專利]一種地球物理勘探中的地震反演方法

有效

• 申請號: 201810903194.4
• 申請日: 2018-08-09
• 公開（公告）號: CN108957554B
• 公開（公告）日: 2021-03-12
• 發明（設計）人: 王揚州;邢濤;李國法;朱廣輝;許開卿
• 申請（專利權）人: 北京探創資源科技有限公司
• 主分類號: G01V1/40
• 分類號: G01V1/40
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種地球物理勘探中地震反演方法，其特征在于，包括：

步驟a，獲取地震工區內的測井數據，對每個地震層相應點的力和力矩數據進行采集和編輯；

步驟b，將每層地震數據進行層層解釋，并將解釋結果數據網格化，獲得基于深度的地質框架模型；

步驟c，根據確定的地質框架模型，確定測井曲線，并對測井曲線的分形參數體進行分形隨機插值，獲得初始測井曲線；

步驟d，根據地震資料確定實際測井曲線，并根據比較結果，對初始測井曲線進行修正，根據修正后的曲線確定地震反演模型；

在上述步驟a中，每個地震層的若干采樣點按照下述模型建立：其中一采樣點建立模型為：

其中，為該地震層采樣點的檢測力和力矩零偏值，為采樣點處自身的重力和重力矩，和為采樣點處所受到的外力和外力矩，其中，力矩為采樣點相對于地震層空洞的中心線的力矩；

設該地震層空洞的中心坐標系原點到采樣點質心的矢量為則

式(1)在基于地震層空洞的中心的測量坐標系中的詳細表達式為：

其中，為從世界坐標系到地震層空洞的中心坐標系的轉換矩陣，由字節初始位置決定，x_c,y_c,z_c；

當采樣點位于該地震層空洞的上表層時：

檢測力與外力存在關系：

F₁表示上表層位置時的力檢測值，m_load表示地質自身的實時重力，如沙礫、沙土自身重量，f_ext外力，F_0x,F_0y,F_0z表示該位置各個方向上的合力，T_0x,T_0y,T_0z表示該位置各個方向上的合力矩，x_c,y_c,z_c表示該位置處的坐標，能夠唯一確定采樣點的位置；

當采樣點位于該地震層空洞的下表層時：

檢測力與外力存在關系：

當采樣點位于該地震層空洞的中間層時：

檢測力與外力存在關系：

通過采集運動過程力檢測測量數據，可同步得到要辨識的力檢測零偏和自身負載重量參數：

通過上述公式(8)最終確定x_c,y_c,z_c位置處的坐標，能夠唯一確定采樣點所受力、力矩及對應采樣點的位置。

2.根據權利要求1所述的地球物理勘探中地震反演方法，其特征在于，在上述步驟d中，按照對選擇的每個點進行修正，按下述公式(10)進行修正；

i_m＝ρ×i₀ (10)

其中，i_m表示修正后的采樣點的瞬時值，ρ表示修正系數，i₀表示采樣點的瞬時值；修正系數ρ按下述公式(11)計算，

式中，ρ表示修正系數，i₀₁和i₀₂表示出現偏差時，偏差位置處的兩個點的瞬時值，N表示采樣次數，k表示采樣序列。

3.根據權利要求1所述的地球物理勘探中地震反演方法，其特征在于，在上述步驟b中，對地震資料進行多套反射層位的解釋和追蹤，解釋層位的質量檢查和修改，并對層位進行全三維追蹤解釋，使層位解釋成果比較精細，做到每道地震都有精細的層位解釋成果。

4.根據權利要求1所述的地球物理勘探中地震反演方法，其特征在于，在上述步驟c中，計算單井的測井曲線模型分形參數，將每一個地震道的分形參數轉換為測井曲線的分形參數體；分別計算單井的測井曲線深度域模型的分形參數的分形參數，通過歸一化和比例加權等方法，將每一地震道地震分形參數標定為測井曲線的分形參數體，這樣就得到了測井曲線每一地震道的分形參數，然后依據下列公式，計算出每一地震道的分形隨機插值的隨機分量funine：

其中，K是標定系數，H是反應層級域疊后地震數據細節和粗糙度的分維參數，σ是層級域疊后地震數據正態分布的方差，G是服從標準正態分布的Gauss隨機變量，ΔX是已知點和插值點之間的距離，a是插值點的一維坐標。