[发明专利]一种用于叠前地震数据参数反演的改进遗传算法

权利要求书

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近似方程计算对应的地震记录计算值。