[发明专利]基于椭球体双重孔隙模型的岩石纵波速度预测方法

说明书

技术领域

本发明涉及一种地震岩石物理领域的速度频散和衰减的预测方法，特别涉及一种包含椭球形状非均匀掺杂体/流体饱和区域的孔隙介质的斑块模型的纵波速度预测方法。

背景技术

近年来工业高速发展，油气资源需求日益增长，寻找新的资源，勘探油区成为我们的首要目标。随着油气资源的大幅度开采，地下油呈现“整体高度分散，局部相对富集”的状态，给开采带来了很大的困难。这就对勘探技术提出了更高的要求，需要提高声波勘探的分辨率，为最佳井位部署提供依据。利用声波的方法来圈闭油气藏和描述储层特性，是目前我们广泛使用的方法之一，其中纵横波速度是我们的重点描述对象。

在石油勘探过程中，常常需要研究声波在复杂地质构造中的传播问题。岩石中孔隙流体的性质发生变化时，会改变纵、横波速度，地震波诱导的孔隙流体局部流动也会产生影响，并造成弹性波的速度频散和衰减。弹性波穿过饱含流体的多孔介质时，在孔隙流体中会诱发压力梯度，导致孔隙流体流动，直到孔隙压力平衡。同时孔隙流体的流动可能发生在不同的尺度下，因此利用岩石物理模型预测储层纵、横波速度，关键在于明确不同尺度下的流体速度场。大部分油气储层都可看为具有双重孔隙特征的介质。在水-油、水-气、气-油饱和模型中，都可以通过其弹性模量、粘性、密度等参数来计算动能，势能和耗散函数，得到岩石纵、横波速度。

然而，现有技术中岩石纵波速度的预测方法的所有模型都将孔隙介质中的斑块饱和单元体等同为一致的同心球体。但现实中的孔隙介质饱和斑块区域往往不会是理想球体，这导致我们建立的动能与耗散方程与实际的有所出入，不可避免地引起地震波速度预测误差，造成预测的不准确。

发明内容

综上所述，确有必要提供一种能够提高预测准确度的岩石纵波速度的预测方法。

一种基于椭球体双重孔隙模型的岩石纵波速度预测方法，包括以下步骤：步骤S10，通过获取测井数据，实验观测数据，得到岩石物理参数，并生成岩石的干骨架模型和孔隙流体模型；步骤S20，以岩石的干骨架模型和孔隙流体模型为基础，建立双重孔隙介质椭球斑块饱和模型，所述双重孔隙介质椭球斑块饱和模型包括三层具有不同孔隙度和不同流体饱和的椭球壳区域模型，每层球壳代表具有不同孔隙特征和不同流体特征的区域，并计算孔隙介质的势能/动能，以及内嵌入体中流体的动能和耗散方程，推导出拉格朗日方程组，并求取纵横波速度；步骤S30，根据平面波分析方法，得到波动方程的频散关系，并得到纵波速度频散和衰减计算公式。

与现有技术相比较，本发明通过引入椭球形斑块饱和模型，并且包含了固体骨架和流体的空间非均匀性，因此其速度预测更具灵活性和贴近实际的特点，可以解决较为复杂的多流体、多孔隙类型岩石纵波预测问题。

附图说明

图1为本发明第一实施例提供的椭球体斑块饱和模型的结构示意图。

图2为椭球体三层斑块孔隙度和流体分布示意图；三层椭球壳具有不同的孔隙度φ₁,φ₂,φ₃和不同流体f₁,f₂,f₃，不同区域之间交界面为S₁,S₂。

图3为椭球体流体饱和斑块主轴半径和角度示意图。a,b,c为椭球的主半径，为椭球面上任意一点对应的空间角度

图4为三层椭球体斑块饱和模型的孔隙介质纵波速度预测实施方式流程图。

图5为球体斑块饱和模型与其他三种方法对Fort Union砂岩纵波速度预测结果对比。纵波频率为5千赫兹，其中BR代表Biot-Rayleigh方法；Johnson代表Johnson2001年提出的方法，White代表修正后的Whtie方法；TLP代表球体斑块饱和方法；Sw表示含水饱和度，Vp代表纵波速度，单位是米/秒(m/s)。

图6为椭球形斑块部分饱和模型与其他三种Fort Union砂岩5千赫兹纵波速度预测结果对比；其中BR代表Biot-Rayleigh方法；Johnson代表Johnson的方法，White代表修正后的Whtie方法；代表椭球体斑块饱和方法，椭球主轴比例r_ba采用统一等效值0.005;代表椭球体斑块饱和方法，椭球主轴比例r_ba随含水饱和度变化；Sw表示含水饱和度，Vp代表纵波速度，单位是米/秒(m/s)。