测井原理毕业论文

路 菁1，2 李 军1

（1.中国石化石油勘探开发研究院，北京 100083；

2.中国石油大学（北京）博士后流动站，北京 102249）

摘 要 鄂西渝东地区下侏罗统为四川盆地典型的陆相页岩气藏，储层矿物组分及有机碳含量是确定该类气藏工程开采难度与有效性的重要指标。为突破常规储层测井评价方法在复杂矿物储层评价中存在的多解性问题，本研究充分挖掘常规测井资料中蕴含的地质信息，以非线性反演与最优化算法为核心思想，综合评价包含有机碳在内的页岩岩石组分与含量，取得了较好的测井评价结果。研究结果完善了页岩气储层测井评价手段，为推进页岩气勘探开发相关技术发展起到了积极的作用。

关键词 常规测井响应 矿物组分 有机碳含量 非线性反演 最优化方法 测井评价

Logging Evaluation of Mineralogical Constituent and

Total Organic Contents for Gas Shale

LU Jing1，2，LI Jun1

（ and Production Research Institute，SINOPEC，Beijing 100083，China；

Center，China University of Petroleum，Beijing 102249，China）

Abstract Formation of lower Jurassic in the western Hubei and East Chongqing is an typical continental facies shale gas reservoir in Sichuan mineral constituents and total organic contents（TOC）are important indicator of the engineering difficulty and its effectiveness for such gas breakthrough the problem of multi －solutions，that always occur when the conventional reservoir logging evaluation methods are used to solve the gas shale reservoir evaluation，this study fully tap the geological information concealed in conventional logging response，use the nonlinear joint inversion and optimization as the core ideas，to evaluate both the mineralogical constituent contents and TOC for gas shale，and achieve a desirable result .This research supplements the logging evaluation methods for gas shale and play a positive role in related technology developments for gas shale exploration and development.

Key words conventional logging；mineralogical constituent；TOC；nonlinear joint inversion；optimization； logging evaluation

鄂西渝东地区是四川盆地周缘页岩气藏有利目标区之一。建南构造位于四川盆地川东褶皱带石柱复向斜中北部，下侏罗统自流井组发育的深湖－半深湖页岩属于典型的陆相页岩气藏。该套页岩区域分布稳定、厚度大、埋藏浅，但相较于海相页岩具有更加频繁的相变特征，储层矿物组分复杂多变。准确把握页岩气储层矿物组分与有机碳含量是后续储层关键参数——脆性与含气性评价的重要基础，也是页岩气测井评价亟待攻克的重点及难点问题。依靠固定的解释模型，采用少部分测井曲线确定储层矿物含量的评价方法，在岩性及矿物较为单一的常规储层中评价效果较好，却无法妥善解决页岩气储层复杂矿物组分与含量的多解性问题。笔者通过深入挖掘各项常规测井资料中蕴含的丰富地质信息，分析建立储层矿物组分模型，以非线性反演与最优化算法评价包含有机碳在内的复杂岩石组分含量，突破了常规测井储层评价的思想，拓展了非常规页岩气储层矿物组分与有机碳含量的测井评价方法，通过实验室岩心全岩组分数据验证，该方法已取得了较好的评价效果。

1 东岳庙段含气页岩岩性及岩石矿物学特征

目标层下侏罗统自流井组东岳庙段泥页岩，区域横向分布稳定，厚度较大，暗色泥页岩厚约60～100m。储层岩性以含灰泥页岩为主，多见灰色粉砂质泥页岩、介壳泥页岩与介壳灰岩夹层（图1）；储层矿物成分以黏土矿物、石英及方解石为主（平均含量分别为％、％、％），同时含有少量长石与黄铁矿。自生矿物的存在，表明东岳庙段所处的沉积环境为有利于有机质富集与保存的还原环境，实验室分析结果显示，储层有机碳以Ⅱ型干酪根为主，平均含量2％～3％；储层孔隙结构以矿物粒间孔为主，同时发育少量粒内孔及溶蚀裂隙，大量因有机质热解产生的纳米孔隙，使储层具有较好的天然气吸附与储集性能。

图1 研究区东岳庙段泥页岩典型岩性

2 常规测井响应评价储层岩石组分

测井响应是被测地层物理特性的宏观表现［1］，在排除井眼与泥浆侵入等影响的情况下，测井响应本质是测井仪器探测范围内所有岩石微观组分物理特性的综合表现，故各类测井响应实际上涵盖了被测地层所有组分的岩石物理信息。充分挖掘、利用常规测井响应中蕴藏的储层信息评价页岩岩石组分，提供了一条除实验室分析和元素俘获能谱（ECS）测井之外的储层评价思路，同时，弥补了岩心实验室分析无法全井段连续、ECS测井数据采集与解释评价成本高昂等问题［2，3］。

常规曲线非线性联合最优化反演算法

目标函数

区别于利用单一或少数测井曲线与储层某一矿物含量建立函数关系、用以评价其含量的方法，利用常规测井信息开展非线性联合最优化反演评价储层矿物组分的方法与步骤，可简要概括如下：首先，需要对实测响应进行预处理，以期得到接近原始储层真实物理特性的校正测井响应；其次，依据岩心观察与常规评价结果得到的初步认识，圈定解释评价井段内存在的岩石组分类型，并确定其初始含量，形成完整的基于原始假设的储层岩石物理体积模型；再次，依据地区经验或理论参数合理选取各组分的测井响应骨架值，以非线性测井响应方程正演各个常规测井响应，并计算关于校正曲线与正演模拟曲线如式（1）所示的目标函数T（Xj）；最后，通过反复迭代调整各矿物组分含量，使目标函数T（X）达到最小值，并将此时的岩石组分与含量模型作为反演的最终结果，即通过解决图2所示的最优化问题，达到求解复杂矿物储层岩石组分与含量问题的目的［4］。

油气成藏理论与勘探开发技术（五）

图2 非线性联合最优化反演算法简图

式中：loggings为第j次迭代后产生的正演曲线组；loggingc为实测曲线经校正产生的校正曲线组；Xj为第j次迭代确定的各个岩石组分含量；W为各测井曲线在目标函数中的权重；α为迭代稳定性控制参数；T（Xj）为反映正演曲线与校正曲线相似程度的目标函数，当该函数达到最小值时，表明正演曲线已逼近校正曲线，此时，即可认为模型求解得到的岩石组分与含量与地层真实情况最为接近。需要说明的是，采用更丰富的测井响应信息，以及岩心分析、常规储层评价取得的地层初步认识等，能够在更大的程度上降低反演算法的多解性。

共轭梯度最优化算法

从上述分析可知，求解页岩复杂岩石组分的测井评价问题，已被转化为求解目标函数T（Xj）最小值的最优化问题。本研究综合考虑目标函数属于多元函数，且测井响应的非线性关系决定了目标函数的非线性特性，故采用共轭梯度法解决目标函数的最优化问题［5］。

对目标函数T（Xj），在极值点X*处作Taylor展开，忽略高效项时，有

油气成藏理论与勘探开发技术（五）

式中：H＝▽2T（X*）为T（X）在X*处的二阶偏导数矩阵。因为X*为极值点，故▽T（X*）＝0，因而

油气成藏理论与勘探开发技术（五）

可见，任何次的函数T（X）在其极值点附近具有二次函数的特征。设T（X）可以表示为如下所示二次函数

油气成藏理论与勘探开发技术（五）

可以证明具有N阶正定矩阵A的n元二次函数，最多可在n维空间中找到n个彼此关于A的共轭方向（向量），且从任意的初始点出发，依次沿这n个共轭方向作不超过n次的一维搜索，就可以求得目标函数T（X）在n维空间内的极小点。采用上述共轭梯度算法回避了因牛顿法及其改进算法需要计算二阶偏导数矩阵的逆矩阵而带来的巨大运算量，且克服了最速下降法在接近极小点时收敛速度很慢的缺陷，妥善地解决了研究建立的非线性反演算法的求解问题。

东岳庙段页岩岩石组分反演

初始模型假设

图3为研究区某井东岳庙段泥页岩常规测井响应，该井含气页岩岩石组分评价的目的在于，明确包括有机碳在内的岩石重要组分的具体含量。初始模型假设的建立，需要分别确定待求解的储层岩石组分及其初始含量，以及参与岩石组分评价的测井曲线。

依据上节所述实验室全岩分析结果（图3），初始模型假设页岩中不存在除干酪根之外的其他固体有机碳；脆性矿物包括石英、方解石、长石，塑性矿物即为黏土；另外，由于相关研究表明，页岩成岩过程中自生的黄铁矿常结晶于储层层理界面之间，在一定程度上有利于水力压裂形成网状缝，且黄铁矿物具有极好的导电特性、极高的光电俘获截面指数以及较高的密度，即使含量较小，对电阻率、光电截面指数与体积密度等测井响应的影响也十分明显，因此，作为影响页岩力学性质与岩石物理特性的重要矿物，黄铁矿在岩石组分模型中不可忽略；最后，由于该段泥页岩黏土矿物含量较高、有效孔隙度较低，且地层水矿化度不高，自由水对测井响应影响不大，故模型仅考虑黏土束缚水存在且假设页岩储层有效孔隙全部被游离气占据的情况。综合上述考虑，最终确定该井东岳庙段泥页岩需要反演计算的岩石组分如图4所示，依次包含黏土（含黏土束缚水）、石英、方解石、长石、黄铁矿、孔隙（游离气）与有机碳（干酪根）。

图3 建南地区某井东岳庙段泥页岩常规测井响应特征与岩心分析结果

图4 页岩岩石体积模型

综合考查本井可参考的测井曲线条数，以及上述页岩岩石体积模型需要涵盖的组分种类，确定利用光电截面指数（PEF）、自然伽马（GR）、中子孔隙度（NPHI ）、体积密度（DEN）、声波时差（DT）、浅侧向电阻率（LLS）、深侧向电阻率（LLD）、铀（URAN）、钍（TH）共9条曲线（图3），反演8种岩石（图4中）组分的含量。可以注意到，如不考虑欠定求解，参加非线性反演的测井曲线条数理论上最多可处理10种岩石组分含量的求解问题，此数大于本模型求解的岩石组分数量，故模型求解结果属于非线性超定解，能够有效降低评价结果的多解性，确保评价结果更加接近页岩气储层的真实情况。

依靠常规储层评价方法，如自然伽马泥质含量Vsb评价方法［6］、密度中子孔隙度Phi评价方法［7］、Pessay有机碳TOC含量评价方法等［8］，可以取得黏土、孔隙度、干酪根含量的初步评价结果，对本井岩石组分初始含量 进行赋值， 剩余组分的初始含量——石英含量 根据实验室岩心分析确定的平均含量（石英％、方解石％、长石％、黄铁矿％）按比例分配，结果如图5第2～8道内实线所示。可以注意到，各岩石组分初始含量（棕色实线）与岩心分析结果（黑色圆点）相比，均存在不同程度的偏差。其中方解石、长石两种矿物含量的偏差最为明显；利用中子－密度孔隙度评价的孔隙度结果也明显偏高；此外，利用自然伽马泥质含量评价方法计算的黏土矿物含量，以及电阻率－声波重叠Passey法计算的有机碳含量，在局部深度上还存在一定误差。本研究将通过随后的反演计算逐步降低这些误差，以得到最接近真实地层岩石组分的评价结果。

图5 建南地区某井东岳庙段泥页岩非线性反演初始模型

模型反演结果

经过非线性反演计算，最终确定该井岩石组分的含量如图6所示，图中第2～8道依次为黏土矿物（含黏土束缚水）、石英、方解石、长石、黄铁矿、孔隙及有机碳含量的评价结果（实线）与对应组分实验室分析结果（黑色圆点），图中第9与第10道分别为页岩岩石组分非线性反演结果与岩心实验室分析结果。

图6 建南地区某井东岳庙段泥页岩岩石组分非线性反演成果图

通过图7各组分初始评价结果（黑色方块）与非线性反演计算结果（三角）的对比分析可以发现，非线性反演结果与实验室分析结果具有更好的线性相关性，与初始评价结果相比更集中于45°对角线附近。图6与图7均显示，非线性反演算法显著提高了石英与方解石含量的评价精度；使孔隙度评价结果更加接近实验室分析结果；此外，黏土矿物与有机碳含量各自在局部位置上的误差也得到了较好的修正；在初始模型中，以平均含量为依据粗略估算的长石与黄铁矿含量，这里也得到了进一步细化，评价结果与实验室分析结果在整体趋势上更为吻合。至此，本研究利用建立的非线性反演方法，同时完成了研究区东岳庙段页岩气储层复杂矿物组分与有机碳含量测井评价两个问题，且取得了较高的评价精度，本研究将进一步定量分析测井评价结果，以验证该方法的可靠性与有效性。

非线性反演结果分析

考虑到各项实测测井响应其本质是被测储层岩石组分反映在各类物理场中的宏观物理特性，因此，为验证非线性反演算法及其反演结果的可靠性与有效性，本研究同时分析了非线性反演结果并在反演结果下模拟了测井响应的误差。

图8展示了非线性反演结果下的模拟测井响应（虚线）与环境校正后的测井响应（黑色实线），涉及的测井项目依次为自然伽马GR、铀Uran、钍Th、中子孔隙度Nphi、体积密度DEN、宏观截面指数U、声波时差DT、冲洗带电导率CXO与原状地层电导率CT。从两组测井响应的对比看，非线性反演结果下的模拟测井响应与实测测井响应具有良好的一致性。表1中定量评分析了两组测井响应间的相关系数，各项测井响应的相关系数在～之间，相关系数均值达到，充分反映了反演结果下的岩石组分宏观物理特性与真实储层物理特性的相似性，即说明通过非线性反演得到的岩石组分及其含量已十分接近页岩气储层的实际情况。此外，以实验室分析结果为标准，表2分别统计分析了图8中初始评价结果与非线性反演结果对实验室结果的相关系数，两组相关系数的对比可以说明，本研究建立的非线性反演算法明显提高了页岩各岩石组分评价的精确度。因此，上述两方面分析充分证明，本研究建立的非线性反演算法在解决页岩储层复杂岩石组分与含量评价问题方面的可靠性与有效性。

该方法能够同时解决页岩气储层岩石矿物组分与有机碳含量评价的两大问题，这两项问题的顺利解决对于后续储层脆性、吸附气含量等重要储层参数评价提供了科学的依据与技术保障。

图7 页岩岩石组分初始评价结果与非线性反演计算结果对比

表1 模拟测井响应与实测响应相关系数

图8 建南地区某井东岳庙段泥页岩复杂岩石组分反演质量控制

表2 初始评价及非线性反演评价较岩心分析结果的相关性对比

3 结论

本研究以非线性反演与最优化算法为核心思想建立的页岩气储层岩石组分测井评价方法，在鄂西渝东建南构造东岳庙段的页岩气储层评价中取得了较好的评价效果。该方法充分挖掘了常规测井资料中蕴含的丰富地质信息，同时解决了页岩储层重要矿物与有机碳含量评价两大问题，弥补了岩心分析深度不连续、ECS测井代价高昂的弊端，且极大地提高了测井评价结果的精度，为后续储层脆性与含气性的综合评价提供了科学的依据与重要的技术保障。

参考文献

［1］洪有密.测井原理与综合解释.东营：石油大学出版社，1993：1 ～365.

［2］Pemper R R，Sommer A，Guo P，et new pulsed neutron sonde for derivation of formation lithology and Annual Technical Conference and Exhibition，San Antonio，Texas，USA，2006：1～13.

［3］Sondergeld C H，Newsham K E，Comisky J T，et considerations in evaluating and producing shale gas Unconventional Gas Conference，Pittsburgh，Pennsylvania，USA，2010：1～34.

［4］Heidari Z，Torres－Verd C，Preeg W method for estimating total organic carbon and porosity，and for diagnosing mineral constituents from well logs in shale－gas ：1 ～15.

［5］雍世和.最优化测井解释.东营：石油大学出版社，1995：87～119.

［6］潘仁芳，伍媛，宋争.页岩气勘探的地球化学指标及测井分析方法初探.中国石油勘探.2009，3：6～28.

［7］朱华，姜文利，边瑞康，等.页岩气资源评价方法体系及其应用——以川西坳陷为例.天然气工业.2009，29（12）：130～134.

［8］Passey Q R，Creaney S，Kulla J B，et practical model for organic richness from porosity and resistivity ，74：1777～1794.