[發明專利]一種用于疊前地震數據參數反演的改進遺傳算法

有效

• 申請號: 201710392537.0
• 申請日: 2017-05-27
• 公開（公告）號: CN107024717B
• 公開（公告）日: 2019-03-01
• 發明（設計）人: 伍慶華;王麗萍;朱志新;顏雪松
• 申請（專利權）人: 伍慶華
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京輕創知識產權代理有限公司 11212
• 代理人: 楊立;李蕾

1.一種用于疊前地震數據參數反演的改進遺傳算法，其特征在于，包括如下步驟：

S1、對油田地下儲層測井區域進行采樣，獲取至少兩個測井采樣點數據，獲取至少兩組彈性參數；對所述油田地下儲層測井區域與所述至少兩個測井采樣點對應的地震記錄采樣點進行地震勘探，獲取至少一組實際地震記錄；

S2、運用所述彈性參數進行迭代運算，更新彈性參數；

S3、通過S2計算得到的彈性參數計算對應的地震記錄計算值，并與所述S1中的實際地震記錄比較，得到驗證計算結果；

S4、當驗證計算結果小于閾值時，驗證成功得到所述彈性參數的反演結果；當驗證計算結果大于或者等于閾值時，驗證失敗則返回S2；

改進的遺傳算法的種群初始化選擇改進策略，初始種群的個體中的前三個元素的值根據所述S1中測井獲取的第一組彈性參數上下浮動隨機選取；

所述初始種群的個體中的第四個元素至最后一個元素的值按如下約束隨機選取：

其中，i是測井層數，n是采樣點的數目，Vp_iwell為步驟S1中的測井第i層的縱波波速，Vs_iwell為步驟S1中的測井第i層的橫波波速，ρ_iwell為步驟S1中的測井第i層的巖石密度，Vp_i為初始化生成的第i層的縱波波速，ΔVp_i為初始化生成的第i層和第i+1層的縱波波速差值，Vs_i為初始化生成的第i層的橫波波速，ΔVs_i為初始化生成的第i層和第i+1層的橫波波速差值，ρ_i為初始化生成的第i層的巖石密度，Δρ_i為初始化生成的第i層和第i+1層的巖石密度差值。

2.根據權利要求1所述的用于疊前地震數據參數反演的改進遺傳算法，其特征在于：

所述S2中通過實數型編碼的遺傳算法對所述彈性參數進行迭代運算，所述遺傳算法中初始種群中的個體的基因位的數值根據所述S1中獲取的若干組彈性參數在預設范圍內上下浮動生成；所述遺傳算法中個體適應度通過所述驗證計算結果進行計算。

3.根據權利要求2所述的用于疊前地震數據參數反演的改進遺傳算法，其特征在于：每一組所述彈性參數均包括地震波縱波速度V_p、地震波橫波速度V_s和巖石密度ρ。

4.根據權利要求3所述的用于疊前地震數據參數反演的改進遺傳算法，其特征在于,

所述S1中對所述油田地下儲層測井區域進行采樣，獲取n+1個測井采樣點的共n+1組彈性參數：

[V_piwell，V_siwell，ρ_iwell]，其中i＝1,2…,n+1，n為正整數；

對所述油田地下儲層測井區域中進行地震勘探獲取與所述n+1個測井采樣點對應的n個地震記錄采樣點的m種不同角度的實際地震記錄如下：

[s(θ_i,j)]，其中i＝1,2…,n；j＝1,2…m，m為正整數。

5.根據權利要求4所述的用于疊前地震數據參數反演的改進遺傳算法，其特征在于：

所述遺傳算法中的每個個體為長度為3n+3的實數型一維數組；初始種群中的個體中每個元素的值根據所述S1中獲取的n+1組共3n+3個彈性參數在預設范圍內上下浮動隨機選取。

6.根據權利要求5所述的用于疊前地震數據參數反演的改進遺傳算法，其特征在于：所述S3中，

通過S1獲取的角度以及S2計算得到的彈性參數運用Zoeppritz方程式計算對應的地震記錄計算值：

[s'(θ_i,j)]，其中i＝1,2…,n；j＝1,2…m。

7.根據權利要求6所述的用于疊前地震數據參數反演的改進遺傳算法，其特征在于：所述S3中，

所述驗證計算結果的計算方式如下：

8.根據權利要求6或7所述的用于疊前地震數據參數反演的改進遺傳算法，其特征在于：所述S3中，

通過S1獲取的角度以及S2計算得到的彈性參數運用Zoeppritz方程式中的Aki&Richard近似方程計算對應的地震記錄計算值。