纵波和转换波叠前联合反演方法及系统技术方案

公开了一种纵波和转换波叠前联合反演方法和系统。该方法包括：输入纵波角度道集和转换波角度道集；提取各角度纵波子波和转换波子波；建立纵波时间域的初始弹性参数模型；计算纵波时间域的纵波合成记录误差；计算时移算子矩阵和误差反传播算子矩阵；计算转换波时间域的转换波合成记录误差；计算纵波时间域的纵波和转换波的总误差，如果总误差小于或等于设定误差门槛值，输出弹性参数反演结果；如果总误差大于设定误差门槛值，更新纵波时间域的对数弹性参数模型。本发明专利技术实现了纵横波精确自动匹配的纵波和转换波叠前联合反演，获取精确稳定的弹性参数反演结果，提高储层预测精度。

Prestack joint inversion method and system of P-wave and converted wave

【技术实现步骤摘要】
纵波和转换波叠前联合反演方法及系统
本专利技术涉及油气地球物理
，更具体地，涉及一种纵波和转换波叠前联合反演方法及系统。
技术介绍
多波地震勘探是进行岩性油气藏和隐蔽油气藏勘探的一种非常有潜力的手段。虽然全球多波多分量地震勘探技术研究和实际应用取得了长足的进步，但一直没有取得显著进展，并且已经成为制约多波地震勘探技术进一步发展的“瓶颈”。多波资料的匹配、多波资料的联合反演、多波综合解释和全波属性的地质应用等问题一直是多波地震解释和储层预测、含气性预测的重点和难点。纵波井震标定主要是用纵波自激自收合成记录匹配纵波叠加记录，但是当纵波垂直入射时，没有转换波，因此不能用转换波自激自收合成记录匹配转换波叠加记录。Stewart(1992)利用某一角度转换波记录与横波自激自收振幅关系把转换波叠加剖面转换为横波叠加剖面，再利用横波自激自收合成记录匹配，实现井震标定。Hampson&Russell软件可以实现某一角度转换波合成记录制作和对全偏移距合成记录叠加两种方法进行转换波井震标定。JamesEG(1996)详细介绍了最大相关法求取纵本文档来自技高网...

【技术保护点】
1.一种纵波和转换波叠前联合反演方法，其特征在于，包括：/n1)输入纵波时间域的纵波角度道集

【技术特征摘要】
1.一种纵波和转换波叠前联合反演方法，其特征在于，包括：
1)输入纵波时间域的纵波角度道集和转换波时间域的转换波角度道集
2)提取各角度纵波子波Wpp和转换波子波Wps；
3)建立纵波时间域的初始弹性参数模型其中，k为迭代次数；
4)计算纵波时间域的纵波合成记录误差
5)计算时移算子矩阵Q和误差反传播算子矩阵QT；
6)计算转换波时间域的转换波合成记录误差
7)计算纵波时间域的纵波和转换波的总误差rk，如果总误差小于或等于设定误差门槛值，转至步骤9)，如果总误差大于设定误差门槛值，执行步骤8)；
8)更新纵波时间域的对数弹性参数模型，重复执行步骤4)～步骤7)；
9)输出弹性参数反演结果
其中，为弹性参数模型，分别为在纵波时间域的对数纵波速度向量、对数横波速度向量、对数密度向量。


2.根据权利要求1所述的纵波和转换波叠前联合反演方法，其特征在于，所述纵波时间域的纵波角度道集和所述转换波波时间域的转换波角度道集分别通过公式(1)和公式(2)表示：






其中，Gpp＝[WppADWppBDWppCD]，Gps＝[0WpsQEDWpsQFD]，

A(θi)、B(θi)、C(θi)为在纵波入射角为θi时纵波反射系数中的纵波速度变化率、横波速度变化率和密度变化率的系数，

E(θi)、F(θi)为在纵波入射角为θi时转换波反射系数中的纵波速度变化率、横波速度变化率和密度变化率的系数，K为接收地震道数目；D为差分算子矩阵，


3.根据权利要求2所述的纵波和转换波叠前联合反演方法，其特征在于，在步骤8)中，采用误差反传播算子变权重共轭梯度迭代法更新纵波时间域的对数弹性参数模型，包括：
8.1)令其中，p(k)为弹性参数更新方向参数，p(0)＝r(0)；和是k次迭代的时移算子矩阵Q和误差反传播算子矩阵QT的函数，即
8.2)计算与p(k)正交的更新方向参数p(k+1)＝r(k+1)+αkp(k)，
其中，

和是k+1次迭代的时移算子矩阵Q和误反传播算子矩阵QT的函数，即



...

【专利技术属性】
技术研发人员：陈天胜，张珺，季玉新，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京;11