[发明专利]一种基于孔隙替换的储层孔隙类型评价的方法

权利要求书

1.一种基于孔隙替换的储层孔隙类型评价的方法，其特征在于，包括以下步骤：

步骤一，根据已知储层的纵波速度、横波速度和密度，计算储层含有原始孔隙类型1的饱和岩石弹性模量；选取储层岩石基质矿物组分的弹性模量和孔隙流体组成成分的弹性模量，结合孔隙度、饱和度和矿物组分体积含量，计算储层岩石基质的弹性模量；计算储层含有原始孔隙类型1的岩石骨架弹性模量；反演含有孔隙类型1时储层的孔隙纵横比；

步骤二，通过给定新的孔隙纵横比，将储层原始孔隙类型1替换成另外一种新的孔隙类型2；计算孔隙类型替换后的储层新的纵波速度、横波速度；

步骤三，对比孔隙类型替换前后储层纵波速度、横波速度的变化情况，进行储层孔隙类型评价；

其中，所述步骤三包括如下步骤：

步骤3-1，根据公式(16)和(17)计算孔隙类型替换前后储层的纵波速度差异和横波速度差异：

ΔV_P＝(V_P2-V_P1)/V_P1(16)

ΔV_S＝(V_S2-V_S1)/V_S1(17)

公式(16)和(17)中，ΔV_P、ΔV_S分别是孔隙类型替换前后储层的纵波速度差异和横波速度差异；

步骤3-2，利用孔隙类型替换前后储层的纵波速度差异ΔV_P和横波速度差异ΔV_S来评价储层孔隙类型，如果纵波速度差异ΔV_P和横波速度差异ΔV_S小于20％，说明原来储层中含有较多孔隙类型2，否则说明原来储层中含有较少孔隙类型2。

2.根据权利要求1所述的方法，其特征在于，步骤一包括如下步骤：

步骤1-1，收集储层的纵波速度、横波速度、密度、孔隙度、矿物组分体积含量和流体饱和度参数；

步骤1-2，根据公式(1)和(2)计算储层含有原始孔隙类型1的饱和岩石体积模量和剪切模量：

Ksat1=ρsat1(VP12-43VS12)---(1)]]>

μsat1=ρsat1VS12---(2)]]>

公式(1)和(2)中，V_P1是纵波速度，V_S1是横波速度，ρ_sat1是密度，K_sat1、μ_sat1分别是含有原始孔隙类型1的饱和岩石体积模量和剪切模量；

步骤1-3，选取储层岩石基质矿物组分的弹性模量和孔隙流体组成成分的弹性模量，结合矿物组分体积含量和饱和度，根据公式(3)、(4)和(5)计算储层岩石基质的体积模量和剪切模量以及孔隙流体的体积模量：

Km=(Σi=1NfiKi+1/Σi=1NfiKi)/2---(3)]]>

μm=(Σi=1Nfiμi+1/Σi=1Nfiμi)/2---(4)]]>

1Kfl=SoilKoil+SgasKgas+SwaterKwater---(5)]]>

式中，K_m是储层岩石基质的体积模量，μ_m是储层岩石基质的剪切模量，K_fl是孔隙流体的体积模量，K_i是第i种矿物组分的体积模量，μ_i是第i种矿物组分的剪切模量，f_i是第i种矿物组分的体积含量，满足N是岩石所有矿物组分的总数，K_oil，K_gas，K_water分别为油、气、水的体积模量，S_oil，S_gas，S_water分别为油、气、水的饱和度，且满足S_oil+S_gas+S_water＝1；

步骤1-4，利用公式(6)和(7)计算含有原始孔隙类型1的储层岩石骨架的体积模量K_dry1和剪切模量μ_dry1：

Kdry1=Ksat(φKmKfl+1-φ)-KmφKmKfl+KsatKm-1-φ---(6)]]>

μ_dry1＝μ_sat1(7)

公式(6)和(7)中，K_dry1、μ_dry1分别是含有原始孔隙类型1的储层岩石骨架的体积模量和剪切模量，φ是储层的孔隙度；

步骤1-5，利用公式(8)反演含有原始孔隙类型1时储层的孔隙纵横比α₁：

Kdry1-Km(1-φφcK(α1))=0---(8)]]>

公式(8)中，P^mi(α₁)为储层原始孔隙类型1与储层岩石基质之间的极化因子，是孔隙纵横比α₁的函数。