专利名称：利用纵波、转换波数据联合反演储层纵横波速度比的方法
技术领域：
本发明涉及地球物理勘探技术，更具体地讲，涉及ー种利用纵波、转换波数据联合反演储层纵横波速度比參数并进行储层含气性或含水性识别的方法。
背景技术：
在油气检测和油藏描述中，介质的弹性參数十分重要，这些參数与岩性和流体成分有夫。描述弹性介质的最主要的3个參数是纵波速度、横波速度和密度。目前地球物理勘探的常规技术是利用纵波AVCKAmplitude Versus O ffset,振幅随炮检距变化)叠前道集反演上述3个參数，反演时考虑了纵波反射振幅随入射角变化的因素。但是大量理论研究和实例应用表明，在有限偏移距的纵波资料中，这3个參数的准确求解存在一定的多解性。随着多波多分量地震勘探技术的不断发展，高质量转换波资料的获�。沟米莶ā⒆徊戏囱菁际醭晌嬷⒄蛊鹄吹男录寄�。纵波、转换波联合反演方法的基本原理是基于Smith和Gidlow( 1987)提出的加权叠加方法，通过PP数据的加权叠加来计算纵波和横波波阻抗以及其它的相关參数。Fatti等(1994)做了进ー步的研究，取消了对于Gardner公式的依赖并且包含了密度项。Stewart等(2001)、Landro等(2006)将加权叠加的思想用于联合反演，给出了实用的纵横波联合反演方法。但利用上述反演方法，反演结果得到的是纵波速度、横波速度和密度的变化率AVp/Vp、AVS/VS、Δ P/P，不能直接获取储层的纵横波速度比值，对储层含不同流体成分的弹性參数特征描述不够直观精确，解释复杂程度较大。
发明内容
本发明的示例性实施例的目的在于克服现有技术中的问题，例如，在于克服现有纵波、转换波储层岩石物理弾性參数联合反演方法存在的上述不足。因此，本发明的示例性实施例提出了ー种利用纵波、转换波数据直接反演储层纵横波速度比的方法，从而能够对储层不同含气性或含水性条件下地球物理响应特征进行精细刻画，并对纵波、转换波数据联合反演的计算方法提出了改进方案。此外，本发明的示例性实施例采用了基于阻尼最小ニ乘解的广义线性反演方法，提高了反演算法的稳定性，提高了计算效率和反演精度。本发明的示例性实施例提供ー种利用纵波、转换波数据联合反演储层纵横波速度比的方法，所述利用纵波、转换波数据联合反演储层纵横波速度比的方法可以包括下述步骤。(I)输入地震数据，从而基于相同成像点位置的纵波成像道集和转换波成像道集，根据目的储层所在时间范围来选取纵波反演时窗[n]和转换波反演时窗[m(2)根据选取的时窗的每ー个纵波成像道集和转换波成像道集，建立纵波反射系数表达式和转换波反射系数表达式
r n D ⑷ 1 avP Δρ f I ΔΓ 4 AVs 2 Δρ) . 2ηRpp(O) = - ^Tii + - +Sin26>
2V Vp P j 2 vp r Fs r P
权利要求
1.ー种利用纵波、转换波数据联合反演储层纵横波速度比的方法，其特征在于，所述方法包括下述步骤 输入地震数据，从而基于相同成像点位置的纵波成像道集和转换波成像道集，根据目的储层所在时间范围来选取纵波反演时窗[Ta，Tb]和转换波反演时窗[Ud]； 根据选取的时窗的每ー个纵波成像道集和转换波成像道集，建立纵波反射系数表达式和转换波反射系数表达式
2.根据权利要求I所述的方法，其特征在于，输入的地震数据包括与ニ维纵波和转换波相关的数据以及与三维纵波和转换波相关的数据中的至少ー种。
3.根据权利要求2所述的方法，其特征在干， 所述成像点位置是输入的地震数据中的与纵波和转换波相关的数据体中的任意ー个相同的共深度点或共转换点位置，或者 所述成像点位置是输入的地震数据中的与纵波和转换波相关的数据体中的一段共深度点或共转换点成像范围，或者 所述成像点位置是输入的地震数据中的全部数据体的所有成像点范围。
4.根据权利要求I所述的方法，其特征在干， 根据地震エ区目标储层的地质层位来进行选取时窗的步骤。
5.根据权利要求I所述的方法，其特征在于，根据输入的地震数据中的与纵波成像道集和转换波成像道集相关的振幅随炮检距变化信息来计算纵波成像道集在选定时窗内给定样点的截距、纵波成像道集在选定时窗内给定样点的梯度以及转换波成像道集在选定时窗内给定样点的梯度。
6.根据权利要求I至权利要求5中的任意一项权利要求所述的方法，其特征在于，所述方法还包括下述步骤 采集包括ニ维多分量和三维多分量中的至少ー种的地震数据资料； 通过ニ维或三维纵波、转换波地震资料处理，以得到包括能够进行储层纵横波速度比反演的纵波、转换波叠前时间偏移道集数据的将被输入的地震数据。
7.根据权利要求I所述的方法，其特征在于，所述方法还包括下述步骤 分别对于其他的成像点位置进行反演以得到相应的纵横波速度比，从而得到整个地震数据的纵横波速度比反演结果； 根据得到的整个地震数据的纵横波速度比反演结果，进行储层预测和气水识别。
全文摘要
本发明提供一种利用纵波、转换波数据联合反演储层纵横波速度比的方法，所述方法采用道集地震反射信息获取的γ变化率与旅行时信息获取的γ变化率无限逼近这一思想，构建反演目标函数，目标函数的关键参数包括纵波截距、纵波梯度、转换波梯度及旅行时信息计算得到的纵横波速度比，这些参数均可由实际输入的纵波、转换波道集的AVO信息和旅行时信息获�。虼四勘旰墓菇ǜ雍侠恚钪辗囱莸淖莺岵ㄋ俣缺鹊圆问哺泳�。
文档编号G01V1/28GK102692645SQ20121017791
公开日2012年9月26日 申请日期2012年6月1日 优先权日2012年6月1日
发明者何光明, 周阿波, 李忠, 杨海涛, 王栋, 谢万学 申请人:中国石油集团川庆钻探工程有限公司地球物理勘探公司
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
