[发明专利]一种变粘土骨架的孔隙度计算方法

权利要求书

1.一种变粘土骨架的孔隙度计算方法，其特征在于，包括以下步骤：

步骤一：对岩心进行物性实验、X衍射全岩矿物实验，获取岩心孔隙度、矿物组分、矿物含量数据；

步骤二：用步骤一得到的X衍射全岩矿物实验数据刻度岩性扫描测井，准确计算矿物剖面；

步骤三：根据步骤一和步骤二得到的岩心孔隙度、矿物剖面与密度曲线、混合骨架密度计算模型、密度孔隙度体积模型相结合反演得到粘土骨架密度值；

步骤四：将步骤三得到的粘土骨架密度值与能谱曲线、密度曲线结合反演得到粘土骨架密度计算模型；

步骤五：将步骤四得到的粘土骨架密度计算模型、步骤二得到的矿物剖面与混合骨架密度计算模型相结合反演得到混合骨架密度值；

步骤六：将步骤五得到的混合骨架密度值结合密度孔隙度体积模型反演得到孔隙度。

2.根据权利要求1所述的一种变粘土骨架的孔隙度计算方法，其特征在于，用步骤一得到的X衍射全岩矿物实验数据刻度岩性扫描测井，得到随深度变化的连续的矿物剖面。

3.根据权利要求1所述的一种变粘土骨架的孔隙度计算方法，其特征在于，根据步骤一和步骤二得到的岩心孔隙度、矿物剖面与密度曲线、混合骨架密度计算模型、密度孔隙度体积模型相结合反演得到粘土骨架密度值；具体将按照如下公式进行：

其中，ρ_ma为混合骨架密度，ρ_sh表示粘土骨架密度，v_sh表示粘土矿物含量，ρ_mai表示第i种矿物的密度骨架且不包括粘土矿物，v_mai表示第i种矿物含量且不包括粘土矿物，Φ表示岩心孔隙度，ρ_b表示测井测量密度值，ρ_f表示地层水的密度，n＞1；

式中已知参数ρ_mai、ρ_f由常用骨架密度参数表获得，v_mai、v_sh由岩性扫描测井确定的矿物剖面获得，Φ由岩心物性实验数据获得，ρ_b由密度曲线获得，求解方程组得到ρ_sh的值。

4.根据权利要求1所述的一种变粘土骨架的孔隙度计算方法，其特征在于，由步骤三得到的粘土骨架密度值，结合其对应深度的能谱曲线的钍元素和钾元素的含量及对应深度的密度曲线的密度值，通过统计分析建立ρ_sh与THOR、POTA、ρ_b的关系式，即粘土骨架密度计算模型

ρ_sh＝a*ρ_b+b*THOR+c*POTA+M

式中为粘土骨架密度，THOR为元素钍的含量、POTA为元素钾的含量、为测井测量密度值,a、b、c为系数,M为常数项；

系数a、b、c的求取根据方程组；如下：

方程组中已知参数ρ_sh1、ρ_sh2，…，ρ_shn为粘土骨架密度值，ρ_b1、ρ_b2，…，ρ_bn为粘土骨架密度值对应的测井测量密度值，(THOR)₁，(THOR)₂，…,(THOR)_n为粘土骨架密度值对应的钍含量，(POTA)₁、(POTA)₂，…，(POTA)_n为粘土骨架密度值对应的钾含量，求解方程组得到a、b、c、M值。

5.根据权利要求1所述的一种变粘土骨架的孔隙度计算方法，其特征在于，将步骤四得到的粘土骨架密度计算模型、步骤二得到的矿物剖面与混合骨架密度计算模型相结合反演得到混合骨架密度值；

式中ρ_ma为混合骨架密度值。

6.根据权利要求1所述的一种变粘土骨架的孔隙度计算方法，其特征在于，将步骤(5)得到的混合骨架密度值结合密度孔隙度体积模型计算孔隙度：

式中Φ为计算孔隙度。