[发明专利]一种叠前地震数据弹性参数反演问题的改进粒子群算法

权利要求书

1.一种叠前地震数据弹性参数反演问题的改进粒子群算法，其特征在于，包括如下步骤：

步骤S1：对油田地下储层测井中的若干点进行采样，获取若干组弹性参数，测量出对应采样点的若干组地震记录值；

步骤S2：对步骤S1得到的所述若干组弹性参数进行初始化操作以获取取值范围，在所述取值范围内随机选取若干组弹性参数；

步骤S3：利用步骤S2中得到的所述随机选取的若干组弹性参数，测量其对应的若干组地震记录值；

步骤S4：基于反演目标函数，将步骤S3得到的地震记录值与所述步骤S1中的地震记录值比较，获得反演目标函数值，若反演目标函数值小于目标函数预设值，则停止计算，输出此时的弹性参数以及对应的地震记录值，否则进入步骤S5；

步骤S5：对步骤S2中得到的所述弹性参数进行迭代运算，更新所述弹性参数的数值，迭代次数增加一次，若迭代次数超过最大迭代次数，则停止计算，输出此时的弹性参数以及对应的地震记录值，否则回到步骤S3。

2.根据权利要求1所述的叠前地震数据弹性参数反演问题的改进粒子群算法，其特征在于：每一组所述弹性参数为纵波速度V_p、横波速度V_s和密度ρ。

3.根据权利要求1所述的叠前地震数据弹性参数反演问题的改进粒子群算法，其特征在于：所述步骤S1中对油田地下储层测井中采样n+1个点，获取n+1组弹性参数：

[V_pi，V_si，ρ_i]，(其中i＝1,2…,n+1)；

测量获得m组地震记录值：

[s(θ_ij)](其中i＝1,2…,n；j＝1,2…m)

其中，V_pi，V_si，ρ_i依次为第i组弹性参数中的纵波速度、横波速度和密度，[s(θ_ij)]为第i层采样点中第j个角度对应的地震记录值，θ为角度。

4.根据权利要求1所述的叠前地震数据弹性参数反演问题的改进粒子群算法，其特征在于：所述改进的粒子群优化算法中的每个粒子是长度为3n+3的实数型一维数组，所述步骤S2中初始化操作获取的取值范围是根据所述步骤S1中测井获取的n+1组共3n+3个弹性参数所定的。