一种快速求取地层横波速度的方法技术

本发明专利技术提供了一种快速求取地层横波速度的方法，属于岩石物理领域。本发明专利技术方法利用常规测井数据和实验室岩石物理数据，将表征岩石孔隙形状的孔隙纵横比作为岩石孔隙度和泥质含量的函数，并基于Gassmann方程和Xu-White岩石物理速度模型求取地层横波速度；所述常规测井数据包括纵波速度Vp、密度ρ、孔隙度φ和泥质含量Vsh；所述实验室岩石物理数据包括弹性模量和密度。利用本发明专利技术方法提高了横波速度估算的精度，便于测井资料的解释和精细储层描述，同时，本发明专利技术方法具有计算速度快、稳定性好的优点，可直接用于测井数据解释工作。

【技术实现步骤摘要】
一种快速求取地层横波速度的方法
本专利技术属于岩石物理领域，具体涉及一种快速求取地层横波速度的方法。
技术介绍
在进行叠前地震反演和叠前地震属性分析中，需要得到准确的纵、横波速度测井曲线，然而由于测量仪器的限制和横波测井费用高昂，在常规测井系列中往往缺乏横波测井，即使进行了横波测井，其数据质量也很差，这给储层预测工作带来了一定的困难。为了从地震资料中获取纵波和横波速度，很多地球物理学者提出了理论模型公式和经验公式。如Han(1986)、Klimentors(1991)、Greenberg(1992)、Castagna(1993)、Goldberg(1998)给出的岩石速度、孔隙度及岩石其它参数之间的经验公式，其缺陷是受地域性影响大，精度较低。Kuster和Toksoz(1974)、Xu和White(1996)给出了储层岩石弹性参数和纵、横波速度之间的理论模型公式，其中Kuster-Toksoz模型在假设孔隙形状为椭球体的前提下通过引入可以任意调整的二维的孔隙表面比将各种尺寸的孔隙考虑到了模型计算中，但是它要求岩石内孔隙是稀疏而孤立的，这就限制了孔隙和孔隙内流体之间的相互作用；其中Xu本文档来自技高网...

【技术保护点】
一种快速求取地层横波速度的方法，其特征在于：所述方法利用常规测井数据和实验室岩石物理数据，将表征岩石孔隙形状的孔隙纵横比作为岩石孔隙度和泥质含量的函数，并基于Gassmann方程和Xu‑White岩石物理速度模型求取地层横波速度；所述常规测井数据包括纵波速度Vp、密度ρ、孔隙度φ和泥质含量Vsh；所述实验室岩石物理数据包括弹性模量和密度。

【技术特征摘要】
1.一种快速求取地层横波速度的方法，其特征在于：所述方法利用常规测井数据和实验室岩石物理数据，将表征岩石孔隙形状的孔隙纵横比作为岩石孔隙度和泥质含量的函数，并基于Gassmann方程和Xu-White岩石物理速度模型求取地层横波速度；所述常规测井数据包括纵波速度Vp、密度ρ、孔隙度φ和泥质含量Vsh；所述实验室岩石物理数据包括弹性模量和密度，所述方法包括以下步骤：(1)输入纵波速度Vp、密度ρ、孔隙度φ和泥质含量Vsh数据，并将数据中的异常值剔除；(2)由固体矿物的弹性模量计算给定孔隙度条件下的干岩石的弹性模量；(3)建立砂岩孔隙纵横比与孔隙度及泥质含量的函数关系，如公式(3)所示：αs＝f(φ,Vsh)＝a+b·φ+c·Vsh(3)其中a，b，c为系数，a的初始值取为0.2，b和c的初始值取为0；(4)计算流体饱和岩石的弹性模量；(5)计算地层纵波速度和地层横波速度(6)调整系数a,b,c的取值，重复步骤(2)至步骤(6)，直到理论计算得到的地层纵波速度与测井数据中实测的纵波速度Vp误差满足终止条件时为止，tol是算法终止条件；(7)使用步骤(6)得到的系数a,b,c，重复步骤(2)至(5)得到地层横波速度所述步骤(6)是一个函数优化的过程，目标函数取为以下形式：其中Vp为测井数据中的纵波速度，为理论计算的纵波速度。2.根据权利要求1所述的快速求取地层横波速度的方法，其特征在于：所述步骤(2)具体如下：使用Xu-White岩石物理速度模型的近似公式(1)计算：K(φ)＝K0(1...

【专利技术属性】
技术研发人员：白俊雨，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11