预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共53页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN114076978A(43)申请公布日2022.02.22(21)申请号202010849459.4(22)申请日2020.08.21(71)申请人中国石油天然气股份有限公司地址100007北京市东城区东直门北大街9号(72)发明人徐旺林李建忠白斌高建荣孙远实付玲(74)专利代理机构北京三友知识产权代理有限公司11127代理人姚亮张德斌(51)Int.Cl.G01V1/28(2006.01)G01V1/30(2006.01)权利要求书8页说明书31页附图13页(54)发明名称一种反演方法、系统、装置及存储介质(57)摘要本发明提供了一种反演方法、系统、装置及存储介质。该反演方法包括:利用基于Fatti近似的叠前稀疏层反演获取目标区尖脉冲化的纵波阻抗反射率反演结果;基于所述纵波阻抗反射率反演结果,确定非零反射率的垂向位置向量;基于所述非零反射率的垂向位置向量,构建分界面约束矩阵;以构建的分界面约束矩阵作为约束进行基于Russell近似的最小二乘反演获取流体因子反演结果、剪切模量反演结果和密度反演结果中的一种或两种的组合。该方法有效解决了常规的基于Russell近似的叠前稀疏层反演稳定性较差的问题,能够获取高精度、高鲁棒性的流体因子和/或剪切模量和/或密度反演结果。CN114076978ACN114076978A权利要求书1/8页1.一种反演方法,其中,该方法包括:利用基于Fatti近似的叠前稀疏层反演获取目标区尖脉冲化的纵波阻抗反射率反演结果;基于所述纵波阻抗反射率反演结果,确定非零反射率的垂向位置向量;基于所述非零反射率的垂向位置向量,构建分界面约束矩阵;以构建的分界面约束矩阵作为约束进行基于Russell近似的最小二乘反演获取流体因子反演结果、剪切模量反演结果和密度反演结果中的一种或两种的组合。2.根据权利要求1所述的反演方法,其中,所述利用基于Fatti近似的叠前稀疏层反演获取目标区尖脉冲化的纵波阻抗反射率反演结果包括:基于Fatti近似方程和角度子波褶积矩阵确定第一AVO正演矩阵;基于所述第一AVO正演矩阵并引入反射系数奇偶分解理论,构建关于弹性参数反射率对应的奇偶分量系数的第一地震数据拟合差项;基于反射系数奇偶分解理论,构建关于弹性参数反射率对应的奇偶分量系数的稀疏层约束项;基于所述第一地震数据拟合差项、所述稀疏层约束项构建基于Fatti近似的叠前稀疏层反演目标函数即基于Fatti近似的关于弹性参数反射率对应的奇偶分量系数的叠前反演目标函数;基于构建的所述目标函数进行反演得到纵波阻抗反射率对应的奇偶分量系数;基于所述纵波阻抗反射率对应的奇偶分量系数,得到所述纵波阻抗反射率数据体即为所述目标区尖脉冲化的纵波阻抗反射率反演结果。3.根据权利要求2所述的反演方法,其中,所述第一AVO正演矩阵为:其中,GFatti为基于Fatti近似方程和角度子波褶积矩阵确定的第一AVO正演矩阵;W(θ1),W(θ2),…W(θK)为角度子波褶积矩阵,K指K个所述角度叠加地震数据;A(θj)、B(θj)和C(θj)均为对角矩阵,其中j=1,2,…K;θj为第j个入射角,其中j=1,2,…K;其中,2CN114076978A权利要求书2/8页4.根据权利要求2或3所述的反演方法,其中,所述第一地震数据拟合差项为:||d-GFattiΓx||2;其中,GFatti为由Fatti近似和子波褶积矩阵确定的第一AVO正演矩阵;d为角度叠加地震数据;Γ为由奇偶分解矩阵确定的块化矩阵;x为弹性参数反射率对应的奇偶分量系数;其中,K指K个所述角度叠加地震数据;θj为第j个入射角,其中j=1,2,…K;其中D为奇偶分解矩阵;其中,D=[De,Do];其中,其中,xI为纵波阻抗反射率对应的奇偶分量系数,xJ为横波阻抗反射率对应的奇偶分量系数,xρ为密度反射率对应的奇偶分量系数;其中,其中,xI_e为纵波阻抗反射率对应的偶分量系数,xI_o为纵波阻抗反射率对应的奇分量系数;3CN114076978A权利要求书3/8页其中,xJ_e为横波阻抗反射率对应的偶分量系数,xJ_o为横波阻抗反射率对应的奇分量系数;其中,xρ_e为密度反射率对应的偶分量系数,xρ_o为密度反射率对应的奇分量系数。5.根据权利要求2-4任一项所述的反演方法,其中,所述稀疏层约束项为:λ|x|;其中,λ为稀疏层约束权重;x为弹性参数反射率对应的奇偶分量系数;其中,其中,xI为纵波阻抗反射率对应的奇偶分量系数,xJ为横波阻抗反射率对应的奇偶分量系数,xρ为密度反射率对应的奇偶分量系数;其中,其中,xI_e为纵波阻抗反射率对应的偶分量系数,xI_o为纵波阻抗反射率对应的奇分量系数;其中,xJ_e为横波阻抗反射率对应的偶分量