[发明专利]一种同时提高解析式稳定性和精度的射线弹性参数反演方法

说明书

技术领域

本发明属于油气田勘探、开发过程中的储层预测和油气检测技术，是一种同时提高解析式稳定性和精度的射线弹性参数反演方法

背景技术

目前叠前弹性反演较多采用的是Connolly(1999）提出的基于Zoeppritz方程近似的、与角度有关的弹性阻抗反演方法(Elastic Inversion，简称EI反演)，该方法假设整个地震剖面不同深度上纵、横波速度比为常数，不符合实际；假设地震波在地层中按照常入射角度传播；随着纵波入射角度的增大，EI反演的误差逐渐增大；在实际应用中弹性阻抗（EI）的幅值随入射角度变化很大，用的归一化方法难以控制（Whitcombe，2002；马劲风，2003），幅值的变化易导致岩性和流体识别的错误。王仰华(2003)和马劲风(2004)分别公开了一种射线路径弹性阻抗反演方法，将弹性阻抗表示为纵波阻抗和纵横波速度比的函数关系式，射线弹性阻抗不需要假设纵横波速度比为常数，也无需进行入射角归一化处理，但在求解弹性参数时，解析式不稳定，抗噪能力差，因此限制了射线弹性阻抗的实际应用。刘力辉等(2011)对射线弹性阻抗方法进行二项式展开，将射线弹性阻抗表达成纵波阻抗和横波阻抗的函数，通过误差推导和试验分析，在小入射角度（小于25°）的情况下，与射线弹性阻抗方法相比，刘力辉的近似射线弹性阻抗方法在稳定性和抗噪能力方面有所改善，但是在入射角大于25°时，与射线弹性阻抗误差大于5%，而且随着入射角的增大误差越来越大。而对于入射角大于25°的弹性反演，这种方法不能有效区分储层物性特征。上述射线弹性阻抗反演具有诸多优势，但都存在解析式不稳定的问题。

发明内容

本发明目的是提供一种提高解析式稳定性和抗噪能力，同时提高入射角大于25°的弹性反演的精度，满足有些地区需要入射角大于25°弹性反演的同时提高解析式稳定性和精度的射线弹性参数反演方法。

本发明通过以下步骤实现：

1）在探区采集地震资料，处理得到地震剖面；

2）在探区井中激发地震纵波和横波，得到纵波速度、横波速度数据；

3）在探区测井得到地层密度；

4）选择标准地震子波或者从井旁地震剖面提取子波，用纵波速度、地层密度合成地震记录；

所述的标准地震子波是雷克子波或克劳德子波.

5）用合成地震记录对步骤1）地震剖面进行层位标定，确定目的层段，并拾取地震反射层位；

6）采用下式计算随入射角变化的射线弹性阻抗合成记录模型，确定射线弹性阻抗随入射角变化的类型；

REI(θ)=αρcosθ(1-4β2α2sin2θ+mβ4α4sin4θ);]]>

公式中：θ为入射角;α为纵波速度；β为横波速度；ρ为密度；m为调节系数；REI（θ）为入射角为θ的射线弹性阻抗；

所述的调节系数取2～6。

以上公式是用标准砂岩模型正演测试改进的射线弹性阻抗公式，式中略去了射线弹性阻抗公式的高次项，提高了解析式稳定性；与刘力辉等人提出的方法比较，提高了入射角大于25°的弹性反演的精度。

7）由步骤5）层位标定得到的时深关系，用均匀采样或者非均匀采样把测井数据由深度域转换为时间域；

8）在时间域采用下式由测井数据计算三个入射角对应的射线弹性阻抗；