[发明专利]一种不依赖于测井资料的绝对波阻抗反演方法

摘要

本发明公开了一种不依赖于测井资料的绝对波阻抗反演方法。首先，从叠前炮集地震数据出发，利用地震数据全部信息通过全波形反演的手段得到长波长P波速度模型；其次，基于长波长P波速度模型通过深时转换以及Gardner经验公式得到低频波阻抗模型；最后，使用低频波阻抗模型约束基于多组变异差分进化算法的波阻抗反演，得到绝对波阻抗模型。本发明直接从地震数据出发，通过利用全波形反演与全局优化算法，实现了在不依赖于测井资料的前提下，对横向强非均质介质绝对波阻抗参数的反演。因为使用了叠前地震数据全部信息，所以可以更好的构建横向连续性好的低频波阻抗模型；并且因为使用了全局优化算法求解波阻抗反演问题，能够反演得到全局最优波阻抗参数。

主权项

1.一种不依赖于测井资料的绝对波阻抗反演方法，其特征在于，首先，从叠前炮集地震数据出发，利用地震数据全部信息通过全波形反演的手段得到长波长P波速度模型；其次，基于长波长P波速度模型通过深时转换以及Gardner经验公式得到低频波阻抗模型；最后，使用低频波阻抗模型约束基于多组变异差分进化算法的波阻抗反演，得到绝对波阻抗模型，具体包括以下步骤：1)采集原始地震数据，然后对采集到的地震数据进行预处理，处理后得到叠前炮集地震数据，记为d_obs(t,x_r；x_s)，其中t表示时间变量；x_r,x_s分别为检波器和震源的位置；处理得到叠后反射地震数据，记为，其中t表示时间变量；n为共深度点道集序号；2)构建长方形网格地质模型，设定离散网格的空间离散间隔、网格规模以及正演模型时间采样间隔；3)给定地质模型P波速度v(x)的初始模型，并指定待优化的目标函数J(v(x))；4)通过最速下降法迭代得到目标函数的全局极小值，进而得到地质模型P波速度的长波长分量vl(x)；5)通过“深时转换”得到低频P波速度模型vlf(t,n)，进而通过Gardner公式得到低频波阻抗模型Ilf(t,n)；6)构建水平层状地质模型，给定地质模型的层数和层的厚度，每层的介质模型参数为波阻抗I；7)通过使用步骤5)的结果确定参数I的搜索空间，并给定波阻抗反演待优化的目标函数F(I(t))；8)使用全局优化算法对步骤7)的目标函数进行优化，求取其全局极小值点；9)选出步骤8)迭代结束后全局优化算法群体中目标函数值最小的个体，将目标函数值最小的个体解码后即得到最终搜索到的最优绝对波阻抗模型。