[发明专利]一种针对无井地震反演的地震速度扰动建模方法

说明书

技术领域

本发明涉及一种建模方法，特别是涉及一种针对无井地震反演的地震速度扰动建模方法，属于石油勘探类地震反演建模技术。

背景技术

地震叠前同步反演技术常用于识别储层及烃类检测，所谓叠前同步反演即利用不同炮检距道集数据以及纵波、横波、密度等测井资料联合反演，得到与岩性、含气性相关的多种弹性参数，综合判别储层物性及含油气性的一种新技术。目前主流的叠前同步反演方法分为两种：基于模型的叠前同步反演和基于地震道的叠前同步反演。

在反演前，这两类方法都要建立一个相对准确的反演模型，初始模型的构建必须以地质理论为依据，以地震层位解释为约束，结合测井资料，才能得到准确的初始模型。在建立模型的过程中，就地质条件而言，应当考虑工区内界面产状的起伏变化，地层厚度的变化及砂层的尖灭、断层、不整合，砂泥岩薄互层中的精细旋回性及韵律变化，沉积相、沉积模式等地质现象的空间展布规律，还应当考虑具体制作的工艺。就内插、外推的数学方法而言，也不能一成不变，应当针对不同的构造与沉积模式，采用不同的数学方法。

在地震剖面上，应当精细地解释好层位，使初始模型顺解释层位横向递推。这是因为以声波测井曲线构造低频信息过程，实际上是利用己知井点的声波测井资料构造任意CDP位置声波测井曲线的过程，属于同一层内的声波时差值，在横向上应当具有相似性和稳定性。层位解释连线时，还要仔细顺波峰或波谷走，不能窜相位。窜相位对层位解释影响不大，但对反演就意味着人工合成地震记录与地震道不匹配，会导致反演后误差增大，对非线性反演，甚至可能导致迭代发散。同一断块上的两井间低频信息内插时，首先，要由解释人员选择反射较好的几个标准反射层进行解释，给出准确的地震解释结果，解释结果越详细，初始模型构造的准确度就越高；任意CDP位置的初始模型由两井的初始模型进行拉伸或压缩后，再由权函数控制其各自在这点的贡献率，然后相加求和，这个权函数一般与井点到任意CDP位置的距离成反比。

上述建模方法普遍应用于有井区域，然而对于很多勘探新区，如南海深水区，水深在2000ms左右，钻井费用一般在亿元以上，整个工区没有一口钻井，这给建模带来了难度，而建模对于后续的反演至关重要。

发明内容

本发明的主要目的在于，克服现有技术中的不足，提供一种针对无井地震反演的地震速度扰动建模方法，特别适用于南海深水区。

本发明所要解决的技术问题是提供操作简单、建模方便、准确可靠、实用性强的针对无井地震反演的地震速度扰动建模方法，依赖速度场的扰动法进行建模，即从纵波模型出发，考虑背景泥岩vp/vs在一定范围内变化，根据vp/vs的取值为常数，转换获得横波模型vs，再让横波模型vs有规律的递变，通过garder公式获得密度模型，以及计算获得纵波阻抗模型和横波阻抗模型，将获得的诸多模型作为反演模型的输入进行反演，从而获得稳定的反演结果，实现有效快速建模，极具有产业上的利用价值。

为了达到上述目的，本发明所采用的技术方案是：

一种针对无井地震反演的地震速度扰动建模方法，包括以下步骤：

1)获得深水无井工区的深度偏移速度场，即纵波模型，用vp表示；根据背景泥岩vp/vs的取值为常数，计算出深水无井工区的横波模型vs；

2)根据深水工区的经验信息，选取vp/vs为[1.90～2.10]之间的常数，选定起点值从vp/vs＝1.90开始，依据纵波模型vp，按照设定步长按序递增方式进行扰动变化，形成纵波模型vp不变、以设定步长为间隔的vp/vs模型阵列，进而获得每一递增序号对应扰动的横波模型vs及其组成的vs横波模型阵列；

3)根据gardner公式求取密度模型den，den＝c1*vp**c2，其中c1、c2为参数；

4)根据公式zp＝vp*den计算纵波阻抗模型zp；

5)根据公式zs＝vs*den计算横波阻抗模型zs；

6)将密度模型、纵波阻抗模型、横波阻抗模型均作为反演模型的输入进行反演。