欢迎来到知嘟嘟! 联系电话:13095918853 卖家免费入驻,海量在线求购! 卖家免费入驻,海量在线求购!
知嘟嘟
我要发布
联系电话:13095918853
知嘟嘟经纪人
收藏
专利号: 2020102039675
申请人: 西华师范大学
专利类型:发明专利
专利状态:已下证
专利领域: 测量;测试
更新日期:2024-01-05
缴费截止日期: 暂无
价格&联系人
年费信息
委托购买

摘要:

权利要求书:

1.基于交叉梯度正则化约束的叠前地震AVA反演方法,其特征在于,包括以下步骤:(1)输入叠前地震数据体;

(2)提取角度地震子波;分入射角进行地震子波提取,获取不同入射角叠前地震数据上的地震子波;

(3)构建初始模型与先验低频约束模型;

具体方法为,基于预处理后的测井数据,由克里金插值法建立初始模型,再进行低通滤波,得到先验低频约束模型;

(4)建立反演目标函数;

具体方法为:

在叠前地震AVA反演目标函数中加入待反演模型参数间的交叉梯度正则化约束项,则有

其中F(R)为目标函数,G=WA,为叠前地震AVA反演的正演矩阵,W为角度地震子波褶积矩阵,A为由系数因子组成的正演系数矩阵,D为叠前地震数据矩阵,R为由Rp、Rs与Rρ组成的模型参数,Rp、Rs与Rρ分别为纵波速度、横波速度与密度的相对变化率,即弹性参数相对变化率, 为矩阵GR‑D的Frobenius范数,Cauchy(R)为模型参数的柯西正则化约束,C为积分算子矩阵,ξprior为对数先验低频约束模型,CG(R)为模型参数的交叉梯度正则化约束,λ、μ与β分别为相应的正则化约束系数;交叉梯度的表达式为:

2 2 2

CG(R)=∫(|▽Rp×▽Rs|+|▽Rp×▽Rρ|+|▽Rs×▽Rρ|)    (3)2

其中▽为梯度算子,×为叉积运算,|▽Rp×▽Rρ|为▽Rp×▽Rρ的绝对值平方;

(5)将叠前地震数据按Inline线号进行矩阵排列;不同入射角的时间采样点按行排列,CMP号按列排列;

(6)采用梯度下降算法对反演目标函数进行迭代求解;

具体方法为,输入叠前地震数据矩阵,角度地震子波褶积矩阵,正演系数矩阵,积分算子矩阵,对数先验低频约束模型,初始模型参数,以及相应的正则化约束系数,采用梯度下降算法对反演目标函数进行迭代求解;

(7)判断迭代是否停止,若停止迭代则输出反演得到的模型参数;

(8)通过弹性参数相对变化率与弹性参数之间的积分关系,将步骤(7)得到的模型参数换算成最终的反演结果,即弹性参数;

(9)对每条Inline线叠前地震数据进行AVA反演,得到整个叠前地震数据体的反演结果。

2.根据权利要求1所述的基于交叉梯度正则化约束的叠前地震AVA反演方法,其特征在于,步骤(1)所述的输入叠前地震数据体,可以是非叠加叠前地震数据,也可以是部分角度叠加地震数据。