[发明专利]泥页岩储层矿物成分含量解析与岩相识别表征方法

权利要求书

1.一种泥页岩储层矿物成分含量解析与岩相识别表征方法，其特征在于，包括：

建立泥页岩储层矿物成分含量的反演解析模型，并对所述反演解析模型求解，获得各矿物的解析计算结果；

基于各矿物的所述解析计算结果，通过矿物组分决策树获得岩相识别结果；

其中，所述基于各矿物的所述解析计算结果之前，还包括：

基于各矿物的体积含量，统计黏土矿物体积相对含量、长英质体积相对含量、灰质体积相对含量的最大值、最小值和平均值，同时统计出各矿物含量分布区间的频数，以构建所述矿物组分决策树；

其中，通过矿物组分决策树获得岩相识别结果，包括：

根据黏土矿物体积相对含量最大值、最小值、平均值及分布区间频数图，划分确定出两个区间，两个区别分别对应黏土含量相对体积为14％-48％和48％-72％；

对于黏土相对体积含量为48％-72％的区间，识别为泥岩相；

对于黏土相对体积含量为14％-48％的区间，结合平均值和频数分布特征，黏土相对体积含量为40％作为细分的界限，则根据黏土相对体积含量细分为14％-40％和40％-48％两个区间；

识别黏土含量为14％-40％的区间，根据灰质相对体积含量的最大值、最小值和平均值，划分出灰质相对体积含量为0％-18％和18％-31％两个区间，其中灰质含量为18％-31％时识别为泥质灰岩相；而灰质含量为0％-18％时，则根据长英质含量的最大值、最小值和平均值，当黏土含量为14％-40％、灰质含量为0％-18％、长英质含量为6％-38％时识别为泥质粉砂岩相，当黏土含量为14％-40％、灰质含量为0％-18％、长英质含量为38％-60％时识别为粉细砂岩相；

识别黏土含量为40％-48％的区间，根据长英质含量的最大值、最小值和平均值细分为三个区间6％-30％、30％-38％和38％-60％，当长英质含量为30-38％时识别为粉砂质泥岩相，当长英质含量为38％-60％时识别为粉细砂岩相；当长英质含量为6％-30％时需结合灰质含量进行细分，当黏土含量为40％-48％、长英质含量为6％-30％、灰质含量为0％-18％时识别为灰质泥岩相，当黏土含量为40％-48％、长英质含量为6％-30％、灰质含量为18％-31％时识别泥质灰岩相。

2.根据权利要求1所述的方法，其特征在于，所述建立泥页岩储层矿物成分含量的反演解析模型，包括：

选取敏感测井曲线建立泥页岩储层矿物成分含量的反演解析模型；其中，所述敏感测井曲线包括声波时差AC、补偿密度DEN和补偿中子CNL。

3.根据权利要求1所述的方法，其特征在于，对所述反演解析模型求解，包括：

根据构建的目标函数，采用“盲源”模拟退火方法进行匹配迭代求解。

4.根据权利要求3所述的方法，其特征在于，所述目标函数为：

式中，F(v)为误差指数曲线；P_j为权重因子系数；K为范数；A_ji为各矿物骨架及孔隙流体的测井响应值；Y_j为测得实际地层的测井值；V_i为各矿物及孔隙的相对体积百分比，即要解析出的各矿物体积含量及孔隙度。

5.根据权利要求3所述的方法，其特征在于，所述采用“盲源”模拟退火方法进行匹配迭代求解，还包括：

针对每种矿物的测井响应值，在所述迭代求解时增加最大值与最小值的限制条件。

6.根据权利要求2所述的方法，其特征在于，所述反演解析模型包括：

式中，DEN为地层的补偿密度测井值；AC为地层的声波时差测井值；CNL为地层的补偿中子测井值；V_石英、V_长石、V_方解石、V_黏土、V_其余矿物、V_孔隙分别为石英、长石、方解石、黏土、其余矿物及孔隙所占体积百分比；AC_石英、DEN_石英、CNL_石英为石英的声波时差、密度、中子测井响应值；AC_长石、DEN_长石、CNL_长石为长石的声波时差、密度、中子测井响应值；AC_方解石、DEN_方解石、CNL_方解石为方解石的声波时差、密度、中子测井响应值；AC_黏土、DEN_黏土、CNL_黏土为黏土的声波时差、密度、中子测井响应值，AC_其余矿物、DEN_其余矿物、CNL_其余矿物为其余矿物的声波时差、密度、中子测井响应值，AC_孔隙、DEN_孔隙、CNL_孔隙为孔隙流体的声波时差、密度、中子测井响应值。