CN104570127B - 一种利用地震波速度同时反演孔隙度和泥质含量的方法 - Google Patents

Abstract

本发明提供了一种利用地震波速度同时反演孔隙度和泥质含量的方法，属于石油与天然气勘探开发领域。本方法包括:(1)输入基础数据;(2)在孔隙度和泥质含量取值范围内生成孔隙度φ和泥质含量β的可行解;(3)计算岩石骨架的等效体积模量K_p00(4)计算岩石骨架的纵波速度V_p0;(5)计算岩石的纵波速度V_p;(6)由地震振幅数据反演得到的纵波速度V_pi和岩石的纵波速度V_p建立同时反演泥质含量及孔隙度的目标函数f(φ,β);(7)利用非线性全局寻优算法对所述目标函数进行迭代反演，直到产生的孔隙度和泥质含量的可行解使得目标函数值满足终止条件为止;(8)输出上述满足终止条件的孔隙度和泥质含量的可行解，即为最终求得的孔隙度和泥质含量。

Description

一种利用地震波速度同时反演孔隙度和泥质含量的方法

技术领域

本发明属于石油与天然气勘探开发领域，具体涉及一种利用地震波速度同时反演孔隙度和泥质含量的方法。

背景技术

地震勘探方法是利用地震波在地下不同介质中传播的速度、振幅、频率、相位、波形等参数的变化来分析、预测油气储层分布范围及储层物性特征的一种地球物理学方法，是目前地球物理勘探中应用最为广泛的方法。随油气勘探开发的不断推进，常规油气勘探已由原来的寻找构造油气藏转向寻找岩性油气藏和隐蔽油气藏，非常规(煤层气、页岩气和致密砂岩气等)油气勘探开发逐渐受到人们的重视，而这类油气储层通常是波阻抗差异小且非均质性强烈，利用传统的波阻抗反演方法不能有效地实现储层预测。孔隙度和泥质含量作为储层参数中的重要参数，是定性、定量评价储层属性的重要参数。它们与地震波速度之间的关系对于揭示地震信息与岩石物性参数的关系意义重大，对储层评价、预测含油气远景区、储量估计及生产开发具有重要的意义。在确定岩石速度时，孔隙度可能是最重要的因素，而孔隙度与埋深和压力之间的关系使速度也随着这两个参数的变化而改变。泥质含量是影响速度的另外一个重要因素，有学者认为：如果岩石孔隙内充填的是泥质，其纵波速度要比充填流体的速度慢30％，会使横波的速度降低40％。在岩石速度与孔隙度等岩性参数之间关系的研究方面，史謌等利用回归分析得出纵、横波速度与孔隙度、泥质含量线性相关。李忠、贺振华等基于地层的孔隙度与速度、波阻抗等地震属性之间存在内在的联系，建立孔隙度和相应地震属性之间的相关关系，将地震属性数据转换成孔隙度数据，从而实现孔隙度预测。吴国平等提出了测井自然伽马求取泥质含量的新原理，利用测井自然伽马信号和提取的自然伽马真值子波，采用维纳滤波法求取测井***逆因子函数，根据逆因子函数，用自然伽马测井信号做卷积直接计算泥质含量。

然而地震波在岩石中的传播速度与岩石的孔隙度、孔隙的几何形状、岩石的矿物组分、岩石颗粒的胶结程度等有关，也与孔隙流体的饱和状态、环境的温度、压力有关，还与波本身的频率有关。在考察泥质含量与地震波速度复杂关系的同时，不可避免地要分析孔隙度对其的影响。现在主要的技术手段还以实验室标本分析为主，得到不同孔隙度、不同饱和度、不同泥质含量的实验室岩石波速值，再经过统计平均得到经验公式应用于实践。泥质含量的计算通常跟孔隙度的计算是不可分割的，目前求泥质含量主要采用测井方法。然而测井数据虽然具有较高分辨率，但是探区内往往钻井数量较少，插值得到的泥质含量剖面横向分辨率无法保证。

发明内容

本发明的目的在于解决上述现有技术中存在的难题，提供一种利用地震波速度同时反演孔隙度和泥质含量的方法,测井数据虽然具有较高分辨率，但是探区内往往钻井数量较少，为此本发明将具有较高纵向分辨率的测井数据和具有较高横向分辨率的地震数据结合，利用岩石物理学方法原理，充分挖掘地震数据中所包含的岩性参数信息，反演得到孔隙度和泥质含量剖面，以辅助岩性(和地层)圈闭油气藏和非常规油气藏(煤层气、页岩气等)的精细储层描述和钻前预测，为油气勘探开发优选目标提供可靠的基础性资料。

本发明是通过以下技术方案实现的：

一种利用地震波速度同时反演孔隙度和泥质含量的方法，包括：

(1)输入基础数据，包括砂岩体积模量K_s、泥岩体积模量K_sh、密度ρ、孔隙流体纵波速度V_pf和从地震振幅数据反演得到的纵波速度V_pi；输入终止条件的值ε，ε由用户给定，是一个较小的数，比如ε＝0.001。

(2)在孔隙度和泥质含量取值范围内生成孔隙度φ和泥质含量β的可行解；

(3)利用步骤(1)输入的砂岩体积模量K_s、泥岩体积模量K_sh、密度ρ及步骤(2)生成的孔隙度φ和泥质含量β的可行解，计算岩石骨架的等效体积模量K_p0；

(4)利用所述岩石骨架的等效体积模量K_p0和密度ρ，计算岩石骨架的纵波速度V_p0；

(5)利用所述岩石骨架的纵波速度V_p0和孔隙流体纵波速度V_pf及步骤(2)生成的孔隙度φ计算岩石的纵波速度V_p；

(6)由地震振幅数据反演得到的纵波速度V_pi和上述步骤(5)计算得到的岩石的纵波速度V_p建立同时反演泥质含量及孔隙度的目标函数f(φ,β)；

(7)利用非线性全局寻优算法对所述目标函数进行迭代反演，直到产生的孔隙度和泥质含量的可行解使得目标函数值满足终止条件为止；

(8)输出上述满足终止条件的孔隙度和泥质含量的可行解，即为最终求得的孔隙度和泥质含量。

所述步骤(2)是这样实现的：

孔隙度φ和泥质含量β是通过式(7)和式(8)生成的：

φ_pos∈[φ_min,φ_max] (7)

β_pos∈[β_min,β_max] (8)

其中，[φ_min,φ_max]是孔隙度的取值范围，[β_min,β_max]是泥质含量的取值范围，φ_pos和β_pos分别为孔隙度和泥质含量的可行解，φ_pos和β_pos在各自的取值范围[φ_min,φ_max]和[β_min,β_max]内随机产生。

所述步骤(3)是利用下式计算岩石骨架的等效体积模量K_p0：

K_p0＝[V_shK_sh+V_sK_s+1/(V_sh/K_sh+V_s/K_s)]/2 (1)

其中，V_sh和V_s分别为泥岩和砂岩的体积百分含量，由公式(2)和(3)计算得到：

V_sh＝(1-φ)β (2)

V_s＝(1-φ)(1-β) (3)。

所述步骤(4)是采用公式(4)计算岩石骨架的纵波速度V_p0：

所述步骤(5)是采用公式(5)计算岩石的纵波速度V_p：

V_p＝1/[φ/V_pf+(1-φ)/V_p0] (5)。

所述步骤(6)的目标函数f(φ,β)为：

f(φ,β)＝||V_p-V_pi||_L (6)，

式中，||·||_L表示范数，L通常取2。

所述步骤(7)中的终止条件为：

f(φ,β)≤ε (9)

所述步骤(7)是这样实现的：

判断f(φ,β)≤ε是否成立，如果否，则返回步骤(2)，如果是，则停止迭代反演。

与现有技术相比，本发明的有益效果是：本发明基于岩石物理学方法原理，建立了速度和孔隙度及泥质含量之间的泛函关系，构建了实测速度与理论计算速度之间的目标函数，利用具有全局寻优能力的优化算法对目标函数进行反演，最后得到孔隙度和泥质含量。

附图说明

图1-1为孔隙度理论模型数据。

图1-2为泥质含量理论模型数据。

图2为纵波速度模型数据。

图3-1为利用本发明方法得到的孔隙度反演结果。

图3-2为利用本发明方法得到的泥质含量反演结果。

图4为本发明方法的步骤框图。

图5为某地区的实际资料反演得到的纵波速度剖面。

图6为本发明反演得到的孔隙度剖面。

图7为本发明反演得到的泥质含量剖面。

具体实施方式

下面结合附图对本发明作进一步详细描述：

本发明基于岩石物理学方法原理，建立速度和孔隙度及泥质含量之间的泛函关系，构建实测速度与理论计算速度之间的目标函数，利用具有全局寻优能力的优化算法对目标函数进行反演，最后得到孔隙度和泥质含量。

如图4所示，本发明具体的步骤如下：

(1)准备基础数据，包括砂岩体积模量K_s、泥岩体积模量K_sh、密度ρ、孔隙流体纵波速度V_pf和从地震振幅数据反演得到的纵波速度V_pi。砂岩及泥岩的体积模量及孔隙流体纵波速度由地球物理实验室做实验观测得到，密度从测井数据得到。

(2)在孔隙度和泥质含量取值范围内生成孔隙度φ和泥质含量β。

(3)计算由砂岩和泥岩组成的岩石骨架的等效体积模量，由公式(1)计算：

K_p0＝[V_shK_sh+V_sK_s+1/(V_sh/K_sh+V_s/K_s)]/2 (1)

式中，K_sh和K_s分别为泥岩和砂岩的体积模量，V_sh和V_s分别为泥岩和砂岩的体积百分含量，K_p0为岩石骨架的等效体积模量。

所述公式(1)中的泥岩和砂岩的体积模量为已知的，体积百分含量V_sh和V_s分别由公式(2)和(3)计算：

V_sh＝(1-φ)β (2)

V_s＝(1-φ)(1-β) (3)。

所述公式(2)和(3)中，φ为孔隙度，β为泥质含量，V_sh和V_s分别为泥岩和砂岩的体积百分含量。

(4)由步骤(3)中计算的岩石骨架等效体积模量K_p0计算岩石骨架的纵波速度V_p0，由公式(4)计算：

所述公式(4)中，ρ是平均密度，由测井数据得到。

(5)根据Wyllie时间平均方程建立岩石速度与孔隙度的函数关系，由公式(5)计算：

V_p＝1/[φ/V_pf+(1-φ)/V_p0] (5)

式中，V_p，V_p0和V_pf(由地球物理实验室做实验观测得到)分别为饱含流体岩石、岩石骨架和岩石孔隙流体的纵波速度，φ为孔隙度。

(6)构建从地震资料同时反演泥质含量及孔隙度的目标函数f(φ,β)，由公式(6)计算：

f(φ,β)＝||V_p-V_pi||_L (6)

式中，||·||_L表示范数，L通常取2。(二范数表示欧几里德空间的距离度量，该误差表示所有离差平方和然后再开根号。本发明取二范数，即L取2。)。V_pi为地震反演得到的纵波速度，此处是作为输入数据给定的。

(7)利用非线性全局寻优算法对上述步骤(6)中的目标函数进行迭代反演。通过与地震反演得到的纵波速度V_pi的比较，来调整孔隙度和泥质含量。具体方式为在孔隙度和泥质含量的取值范围内以随机的方式产生孔隙度和泥质含量的可行解，利用步骤(3)-(4)计算岩石骨架的体积模量和纵波速度，利用步骤(5)计算岩石的纵波速度，重复上述过程，直到产生的孔隙度和泥质含量可行解使得目标函数值满足终止条件为止。孔隙度和泥质含量可行解生成方式和终止条件分别由下式(7)、(8)及(9)给出：

φ_pos∈[φ_min,φ_max] (7)

β_pos∈[β_min,β_max] (8)

f(φ,β)≤ε (9)

所述公式(7)和(8)中，[φ_min,φ_max]是孔隙度的取值范围(储层孔隙度大于0，小于1，[φ_min,φ_max]在0和1之间由用户任意取值)，[β_min,β_max]是泥质含量的取值范围(储层泥质含量大于0，小于1，[β_min,β_max]在0和1之间由用户任意取值)，φ_pos和β_pos分别为孔隙度和泥质含量的可行解，φ_pos和β_pos在各自的取值范围[φ_min,φ_max]和[β_min,β_max]内随机产生。

所述公式(9)中，ε为较小的数，如ε＝0.001，当孔隙度和泥质含量的可行解使得终止条件公式(9)成立，那么算法终止，此时的孔隙度和泥质含量的可行解就是反演得到的解。

本发明方法利用岩石物理学方法原理，充分挖掘地震数据中所包含的岩性参数信息，以地震波速度数据为基础，反演得到孔隙度和泥质含量剖面，以辅助岩性(和地层)圈闭油气藏和非常规油气藏的精细储层描述和钻前预测，为油气勘探开发优选目标提供可靠的基础性资料。

下面以某地区的实际地震数据为例，进行孔隙度和泥质含量的反演，说明本发明的效果。

图1-1至图3-2为理论模型数据的反演，其中，图1-1和图1-2分别为孔隙度剖面和泥质含量模型剖面，图2为纵波速度，图3-1和图3-2分别为反演得到的孔隙度和泥质含量剖面。通过将图3-1、图3-2的反演结果与图1-1、图1-2中的模型对比，可以看出，本发明方法能够在岩性指示及岩性分界面方面有良好反演效果。

图5-图7为实际地震数据的反演。其中，图5是经过地震反演的得到的速度剖面，本发明方法以该速度剖面为基础同时反演得到孔隙度和泥质含量剖面。图6是反演得到的孔隙度剖面，图7是反演得到的泥质含量剖面。

实施例：

本实施例是本发明用于某地区的实际地震数据的孔隙度和泥质含量反演，具体步骤如下：

(1)准备基础数据，包括砂岩体积模量K_s、泥岩体积模量K_sh、密度ρ、孔隙流体纵波速度V_pf和从地震振幅数据反演得到的纵波速度V_pi。砂岩及泥岩的体积模量及孔隙流体纵波速度由地球物理实验室做实验观测得到，密度从测井数据得到。

(2)在孔隙度和泥质含量取值范围内生成孔隙度φ和泥质含量β。

(3)利用上述步骤得到的砂、泥岩体积模量、密度及生成的孔隙度、泥质含量数据，计算岩石骨架的等效体积模量K_p0。

(4)利用岩石骨架等效体积模量K_p0和密度ρ，计算岩石骨架的纵波速度V_p0。

(5)由上述步骤计算得到的岩石骨架的纵波速度V_p0和岩石孔隙流体的纵波速度及孔隙度计算岩石的纵波速度V_p。

(6)由从地震资料反演得到的纵波速度V_pi和上述步骤(5)计算得到的岩石的纵波速度建立同时反演泥质含量及孔隙度的目标函数f(φ,β)。

(7)利用非线性全局寻优算法对上述步骤(6)中的目标函数进行迭代反演。通过与地震反演得到的纵波速度V_pi的比较，来调整孔隙度和泥质含量。重复上述步骤(2)-(6)，直到产生的孔隙度和泥质含量可行解使得目标函数值满足终止条件为止。

(8)输出上述满足条件的孔隙度和泥质含量。

本发明可应用于石油地球物理勘探中的储层预测领域。本发明方法基于岩石物理学方法原理，建立了地震波速度与孔隙度和泥质含量之间的泛函关系，并结合实测速度构建了同时反演孔隙度和泥质含量的目标函数，采用非线性全局寻优算法对目标函数进行迭代优化计算。该方法首先通过岩石骨架等效体积模量建立了岩石骨架纵波速度与泥质含量的关系，然后利用Wyllie时间平均方程建立了岩石速度与孔隙度的关系，最后结合速度剖面构建了同时反演孔隙度和泥质含量的目标函数。该方法能够从叠后地震剖面反演得到孔隙度和泥质含量剖面，对于岩性油气藏及隐蔽油气藏的勘探具有重大的应用价值。

上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发明上述具体实施方式所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的意义。

Claims (8)

1.一种利用地震波速度同时反演孔隙度和泥质含量的方法，其特征在于:所述方法包括:

(1)输入基础数据，包括砂岩体积模量K_s、泥岩体积模量K_sh、密度ρ、孔隙流体纵波速度V_pf和从地震振幅数据反演得到的纵波速度V_pi;输入终止条件的值ε;

(2)在孔隙度和泥质含量取值范围内生成孔隙度φ和泥质含量β的可行解;

(3)利用步骤(1)输入的砂岩体积模量K_s、泥岩体积模量K_sh、密度ρ及步骤(2)生成的孔隙度φ和泥质含量β的可行解，计算岩石骨架的等效体积模量K_p0;

(4)利用所述岩石骨架的等效体积模量K_p0和密度ρ，计算岩石骨架的纵波速度V_p0;

(5)利用所述岩石骨架的纵波速度V_p0和孔隙流体纵波速度V_pf及步骤(2)生成的孔隙度φ计算岩石的纵波速度V_p;

(6)由地震振幅数据反演得到的纵波速度V_pi和上述步骤(5)计算得到的岩石的纵波速度V_p建立同时反演泥质含量及孔隙度的目标函数f(φ,β);

(7)利用非线性全局寻优算法对所述目标函数进行迭代反演，直到产生的孔隙度和泥质含量的可行解使得目标函数值满足终止条件为止;

(8)输出上述满足终止条件的孔隙度和泥质含量的可行解，即为最终求得的孔隙度和泥质含量。

2.根据权利要求1所述的利用地震波速度同时反演孔隙度和泥质含量的方法，其特征在于:所述步骤(2)是这样实现的:

孔隙度φ和泥质含量β是通过式(7)和式(8)生成的:

φ_pos∈[φ_min,φ_max](7)

β_pos∈[β_min,β_max](8)

其中，[φ_min,φm_ax]是孔隙度的取值范围，[β_min,β_max]是泥质含量的取值范围，φ_pos和β_pos分别为孔隙度和泥质含量的可行解，φ_pos和β_pos在各自的取值范围[φ_min,φ_max]和[β_min,β_max]内随机产生。

3.根据权利要求2所述的利用地震波速度同时反演孔隙度和泥质含量的方法，其特征在于:所述步骤(3)是利用下式计算岩石骨架的等效体积模量K_p0:

K_p0＝[V_shK_sh+V_sK_s+1/(V_sh/K_sh+V_s/K_s)]/2(1)

其中，V_sh和V_s分别为泥岩和砂岩的体积百分含量，由公式(2)和(3)计算得到:

V_sh＝(1-φ)β(2)

V_s＝(1-φ)(1-β)(3)。

4.根据权利要求3所述的利用地震波速度同时反演孔隙度和泥质含量的方法，其特征在于:所述步骤(4)是采用公式(4)计算岩石骨架的纵波速度V_p0:

$V_{p0}=\sqrt{K_{p0}/\mathrm{\ρ}}\·\·\·\left(4\right)\·$

5.根据权利要求4所述的利用地震波速度同时反演孔隙度和泥质含量的方法，其特征在于:所述步骤(5)是采用公式(5)计算岩石的纵波速度V_p:

V_p＝1/[φ/V_pf+(1-φ)/V_p0](5)。

6.根据权利要求5所述的利用地震波速度同时反演孔隙度和泥质含量的方法，其特征在于:所述步骤(6)的目标函数f(φ,β)为:

f(φ,β)=||V_p-V_pi||_L

其中，L取2。

7.根据权利要求6所述的利用地震波速度同时反演孔隙度和泥质含量的方法，其特征在于:所述步骤(7)中的终止条件为:

f(φ,β)≤ε(9)。

8.根据权利要求7所述的利用地震波速度同时反演孔隙度和泥质含量的方法，其特征在于:所述步骤(7)是这样实现的:

判断f(φ,β)≤ε是否成立，如果否，则返回步骤(2)，如果是，则停止迭代反演。