一种高分辨率的叠后地震波阻抗反演的方法技术

本发明专利技术公开了一种高分辨率的叠后地震波阻抗反演的方法，包括：根据一维非平稳地震信号的时域表达式和反射系数与波阻抗的近似表示式，构建正演模型；利用纵波阻抗与反射系数关系式和所述反射系数与波阻抗的近似表示式，建立纵波阻抗低频模型；根据拉格朗日乘子法、所述正演模型以及所述纵波阻抗低频模型，构建带低频约束的L

【技术实现步骤摘要】
一种高分辨率的叠后地震波阻抗反演的方法
本专利技术涉及叠后地震反演
，特别涉及一种高分辨率的叠后地震波阻抗反演的方法。
技术介绍
叠后地震反演，是利用经过叠加后的地震数据，经反演得出反射系数、波阻抗信息，侧重于对储层特征进行描述。近30年来，叠后地震反演取得了巨大的进展，并被应用于储层预测、油气藏特征描述等方面。总的来说，叠后地震反演分为基于反射系数逆公式的直接反演和基于正演模型的迭代反演，其中，直接反演法的应用最为广泛，尤其在Theune等人(2010)指出“块状”的反演结果有利于刻画地层的边界后，基于L1正则化的方法及其改进方法在叠后地震反演中得到了成功的应用。Wang等人采用L1正则化和梯度下降法实现了地震波阻抗的反演；Zhang等人发展了基于全变差分的地震波阻抗反演方法，该方法分别获得了相对波阻抗和绝对波阻抗；汪玲玲等人采用Toeplitz矩阵稀疏分解的方法实现了相对波阻抗的反演，并取得了较高的精度，但上述方法采用的都是平稳变化的单极子子波库。平稳变化的单极子子波库会导致反演结果垂向分辨率低。为了改善这一问题，Zhang等本文档来自技高网...

【技术保护点】
1.一种高分辨率的叠后地震波阻抗反演的方法，其特征在于，包括：/n根据一维非平稳地震信号的时域表达式和反射系数与波阻抗的近似表示式，构建正演模型；/n利用纵波阻抗与反射系数关系式和所述反射系数与波阻抗的近似表示式，建立纵波阻抗低频模型；/n根据拉格朗日乘子法、所述正演模型以及所述纵波阻抗低频模型，构建带低频约束的L

【技术特征摘要】
1.一种高分辨率的叠后地震波阻抗反演的方法，其特征在于，包括：
根据一维非平稳地震信号的时域表达式和反射系数与波阻抗的近似表示式，构建正演模型；
利用纵波阻抗与反射系数关系式和所述反射系数与波阻抗的近似表示式，建立纵波阻抗低频模型；
根据拉格朗日乘子法、所述正演模型以及所述纵波阻抗低频模型，构建带低频约束的L1-2最小化反演目标函数：
利用凸函数差异算法DCA和交替方法乘子法ADMM对所构建的带低频约束的L1-2最小化反演目标函数进行求解，得到纵波阻抗。


2.根据权利要求1所述的方法，其特征在于，所述根据一维非平稳地震信号的时域表达式和反射系数与波阻抗的近似表示式，构建正演模型包括：
根据一维非平稳地震信号的频域表达式，得到一维非平稳地震信号的时域表达式：
根据所述一维非平稳地震信号的时域表达式和反射系数与波阻抗的近似表示式，构建正演模型，该正演模型的表达式为：s＝WA(Q)Dm+n；
式中，所述ω代表角频率，τ代表时移，代表w(t)的傅里叶变换，r(t)是反射系数，i是虚数单位，t是采样点，a(τ,ω)代表地震信号的振幅衰减，Q表示层间Q值，W是由子波生成的Toepliz矩阵，A(Q)代表黏弹性衰减矩阵，m是波阻抗的自然对数，D为一阶差分矩阵，s表示地震记录，n表示噪声序列。


3.根据权利要求2所述的方法，其特征在于，所述利用纵波阻抗与反射系数关系式和所述反射系数与波阻抗的近似表示式，建立纵波阻抗低频模型包括：
根据所述纵波阻抗与反射系数关系式和所述反射系数与波阻抗的近似表示式，构建纵波阻抗的自然对数序列与纵波阻抗序列的关系式；
根据低通滤波矩阵和所述纵波阻抗的自然对数序列与纵波阻抗序列的关系式，建立纵波阻抗低频模型。


4.根据权利要求3所述的方法，其特征在于，所述低通滤波矩阵包括：



式中：n、m为低通滤波矩阵L中对应元素的行和列；N为待滤波信...

【专利技术属性】
技术研发人员：聂文亮，冯丽源，相菲，曹湘，刘毓，文晓涛，聂祥飞，
申请(专利权)人：重庆三峡学院，
类型：发明
国别省市：重庆;50