碎屑岩储层成岩相测井识别方法综述及研究进展

赖锦，王贵文,2，王书南，信毅，吴庆宽，郑懿琼，李鉴伦，苍丹

(1. 中国石油大学(北京) 地球科学学院，北京，102249；2. 中国石油大学(北京) 油气资源与探测国家重点实验室，北京，102249；3. 中国石油塔里木油田公司勘探开发研究院，新疆 库尔勒，841000；4. 大庆油田有限责任公司第八采油厂，新疆 克拉玛依，834000)

成岩相是成岩环境的物质表现，是沉积物在特定沉积和物理化学环境中，在成岩与构造等作用下，经历一定成岩作用和演化阶段的产物，包括岩石颗粒、胶结物、组构和孔洞缝等综合特征，通常成岩相包括2个方面的内容，即成岩环境及在该环境下的成岩产物[1−4]。现阶段我国大多数陆上油田已进入以岩性—地层等隐蔽油气藏为主的勘探阶段[5−6]，在隐蔽油气藏勘探的实践中人们逐渐认识到储集体的质量即物性受构造相、沉积相和成岩相的综合影响和控制[7−8]。其中，构造相、沉积相是决定储层物性的地质基础，而成岩相是决定优质储层分布的关键[9]。在一定的构造和沉积背景下，成岩相是决定储集层性能和油气富集的核心要素之一，通过成岩相的研究有助于储层的区域评价和预测[1]。国内外专家学者们有鉴于此，针对自己研究工区的地质特点，充分利用手头所拥有的资料从不同的角度出发对储层成岩相展开了卓有成效的研究工作，在成岩相研究的理论和方法上取得重要成果、认识和进展，在油气勘探的实践中也得到了广泛应用，概括起来主要体现在 4个方面：(1) 研究孔隙演化特征，为恢复构造演化史提供一定依据[10−11]；(2) 寻找预测优质储层发育带及含油气有利区[7−8,12−17]；(3) 识别层序边界，指导层序地层学研究[18−22]；(4) 应用于岩相古地理学，辅助以古生物研究，指示特殊的沉积环境和岩相类型[23−24]。

勘探目标的转变推动了成岩相的研究进展，而成岩相的研究对成岩相的测井识别方法提出了严峻挑战。近年来众多专家学者致力于成岩相的研究也逐渐达成了这样一个共识，即只有搞清楚成岩相的展布规律才能实现优质储集体预测的目标[25]。但是限于取心技术和出于成本的考虑，一个地区的取心井和取心井段总是有限的，要将成岩相的研究由点向线和面上拓展，就必须依赖作为地下地质信息载体的测井、地震资料。因此成岩相研究的关键就是通过测井相和地震相分析取心井段以外的成岩相类型，通过测井、地震资料得到成岩相剖面与平面展布规律。地震资料由于垂向分辨率差，直接用地震资料来划分储层成岩相目前还难以做到，有根据地震波反射结构特征和反射单元的几何形态与沉积相之间相关性建立地震相模式来判别成岩特征的，但其精度一般不高，尚无法达到储层的精细评价和预测[2,25]。相比较而言，测井资料在成岩相的研究中具有不可取代的作用，一方面是由于测井资料纵向分辨率较高、连续性好、获取时间短和成本低，另一方面成岩作用差异造成的岩石成分和结构的变化在常规及成像等测井资料上具有明显不同的响应特征[26−27]。总体而言，目前国内外关于储层成岩相的测井识别研究完全处于初始探索阶段。因此，总结成岩相的测井响应特征以及建立不同成岩相的测井识别模式与准则，利用测井资料将成岩相推到面上是现阶段成岩相研究亟待解决的关键问题。

针对以上研究现状，本文作者总结了前人对成岩相的测井识别方法的研究，认为成岩相的测井识别方法最关键的是在储层岩性识别的基础上成岩作用类型和强度的定性判定以及成岩矿物含量的定量精确计算，并提出了利用蜘蛛网图识别成岩相的方法。以期理论认识与技术的进步能够抵消勘探复杂性增加所带来的不利影响。

1 成岩相与构造相、沉积相关系探讨

研究表明，成岩相在一定的构造、沉积背景下形成的，要进行成岩相的测井识别，首先要搞清楚成岩相的形成及分布的控制因素，确定沉积物沉积、成岩条件及其岩性岩相特征。在现今油气勘探阶段，主导的“相控论”油气勘探理论强调“选凹定相”，认为油气藏都处于有利相带中，是若干有利“相”耦合作用的产物，寻找富油气聚集区带的核心是寻找有利储集相带[28]。也就是说，在烃源岩和构造沉积相等成藏地质条件基本明确的情况下，寻找受成岩相控制的优质储集体是现阶段油气勘探工作的主要目标[29]。

构造作用和沉积作用直接决定了成岩作用的形式和内容[4]。构造演化控制了沉积作用及各种沉积体的发育[30]，而沉积相是影响成岩相分布以及储集性能的地质基础[31]，不同的沉积体系形成于不同构造背景和气候环境[4]，往往发育不同的成岩相，而在同一沉积体系中，沉积物原始碎屑组分和结构以及其他沉积岩石学参数(包括粒度、分选性、磨圆度等)的差异都会对成岩相产生影响[32]。研究表明，优质储层一般分布于水动力较强的高能沉积环境中，因此成岩相研究应与沉积环境分析相结合才能够更好地实现优质储集体的预测，应将有利沉积相带与有利成岩相带结合来划分有利的沉积—成岩相带，这样可以将控制优质储层形成的主要因素均考虑进去，从而更加客观地评价和预测优质储层[33]。构造作用对成岩相的直接控制作用，体现在2个方面：一是通过控制沉积物埋藏史、古地温和酸碱度来影响成岩演化史；二是构造应力作用形成大型断裂或者裂缝从而为孔隙流体流动提供运移通道[30]。

综上所述，构造相和沉积相控制着宏观上油气藏的分布规律，是宏观尺度的地质相概念，而成岩相是相对微观尺度的地质相概念，控制着微观上油气藏内部的非均质性和含油气性，通常构造相控制着沉积相、沉积相控制着岩相、岩相控制成岩相，而四者又共同控制着岩石物理相，并最终控制油气的富集[5−6]。

2 成岩相研究内容及其分类命名

2.1 成岩相研究内容

成岩相的研究一般围绕其内涵和外延的展开，内涵即成岩作用、成岩矿物和成岩环境等，外延则包括成岩阶段、成岩事件和成岩演化序列等[1]。一般要考虑的是沉积物所经历的成岩作用，所处的成岩阶段、成岩环境、成岩过程中具有指示意义的矿物标志、主要成岩事件和成岩演化序列，核心是成岩矿物和成岩作用。

成岩作用可分为建设性成岩作用(溶蚀作用和破裂作用等)和破坏性成岩作用(压实作用和胶结作用等)，但也有些胶结作用，如早期的绿泥石环边胶结一般能增加岩石抗压实能力并且抑制石英的次生加大，能有效地保护原生孔隙，对储层储集物性起建设性作用。成岩矿物既包括成岩作用以前的矿物，也包括成岩过程中新生成的，且能够对成岩相起到指示作用的矿物，如石英、长石、方解石、白云石、沸石以及自生黏土矿物。成岩环境则主要指沉积物埋深和温度条件下，也包括地层水酸碱度、氧化还原性质等。成岩阶段一般根据石油天然气行业成岩阶段划分标准可分为同生期、早成岩、中成岩和晚成岩阶段等。可结合岩石结构、孔隙类型、自生矿物分布及形成顺序、黏土矿物类型及混层黏土矿物的演化结果、镜质体反射率、古地温特征等标志，确定沉积物所处的成岩阶段[9]。成岩演化序列是碎屑组分的成岩变化以及自生矿物的类型、产出形式、先后次序以及各种自生矿物相互的交代关系。成岩事件指成岩的过程、顺序和强度等，一般具共存性和继承性[33]。

2.2 成岩相划分及命名

对于储层成岩相研究的一般思路和基本步骤是：首先分析沉积物构造、沉积和成岩条件，确定成岩相形成背景→微观成岩特征的镜下观察→成岩作用、成岩矿物、成岩阶段、成岩环境和成岩演化序列等的确定→成岩相的分类及命名→通过岩心薄片资料刻度测井确定各成岩相测井、地震响应特征→利用测井、地震资料对成岩相进行剖面展开和平面成图→根据有利成岩相的展布预测有利储集体分布。在前期工作如镜下微观成岩特征的确定的基础上对成岩相进行分类命名是成岩相测井识别的前提条件，只有按照一定的法则将储层成岩相划分为不同的类型，并且划分好的不同的成岩相要具有不同的成岩作用、成岩矿物及成岩演化特征，表现在其岩石学、矿物学特征和孔隙微观几何特征的差异[9]，才可以找到不同成岩相的测井响应特征，建立成岩相的识测井别标准，利用测井资料来识别不同的成岩相。目前虽尚未形成完全统一的成岩相分类命名方案，但大多涉及成岩作用、成岩矿物和成岩环境等，分歧不是很大[1]。总的看来，前人划分成岩相时主要考虑的是成岩作用类型及强度、成岩矿物及其对储层储集物性的影响，一般按照控制成岩特征的主要成岩作用，成岩环境等以优势相或特殊相对成岩相进行分类命名，如表1所示[1,2,4,26,34−40]。

表1 国内外成岩相分类研究现状Table 1 Division and research status of diagenetic facies at home and abroad

3 成岩相测井识别原理

沉积岩是由沉积物被埋藏以后经过一系列物理的、化学的、生物的作用即成岩作用形成的，主要的成岩作用类型有压实、胶结、重结晶、交代、水化等，成岩作用的差异导致了沉积物在粒度、形状、表面结构、取向、矿物成分、孔隙度和渗透率等诸方面发生变化[35]。成岩作用不仅影响了储集空间的演化过程和孔隙结构特征，还最终决定了现今储层的面貌和分布规律[41]。

目前通过测井资料直接识别储层成岩相尚存在一定局限性，但通过综合精细分析不同的测井曲线信息，可以总结出测井资料对不同储层成岩相的响应特征从而达到测井识别成岩相的目的。这是由于测井技术获取的地层信息主要包括地层岩石各种物理性质，如密度、电阻率、含氢指数、声波传播速度等[26]。在岩石骨架矿物基本确定的情况下，地层的密度、电阻率、声速和放射性等物理特性是岩石成岩强度的表现，即岩石胶结程度、压实程度和次生孔缝发育程度的表现[25]。根据成岩作用、成岩矿物等划分出的不同类型成岩相在结构、矿物成分、物性上的差异，导致了它们在测井曲线上具有不同响应特征[42]。因此，可根据同类岩石的测井特性的差异，可以定性和定量判别岩石的成岩强度，找到成岩相的测井响应[25]。通过测井相分析取心井段以外的成岩相类型，利用测井资料实现成岩相的连续划分，得到成岩相的剖面与平面展布规律，从而利用有利成岩相的展布优选优质储层分布区，最终达到储层的区域评价与预测的目标，以便更好地指导油气勘探。与传统的方法相比，用测井资料识别成岩相的方法具有快速经济精确的优势[43]。

4 典型成岩相的识别标志和测井响应特征

前已述及，成岩作用是沉积环境的继续和发展[47]，是储集层发育和形成的必经过程，决定了储层的最终面貌。不同成岩相的本质区别在于其成岩序列、成岩环境、成岩模式不同[44]，是沉积物经历不同的成岩作用、成岩环境和成岩演化序列而形成的[26]，最终表现为其岩石学特征及成岩特征(包括压实和溶蚀组构和强度、胶结物成分与胶结类型、孔隙类型及分布等)不同，物性以及孔隙结构特征具明显差异[45]。因此，利用测井资料识别成岩相的关键一方面就是找到沉积物由于经历不同成岩作用过程以及成岩演化序列之后现今的地质特征的差异，另一方面就是通过沉积物现今地质特征的差异反演找到导致这些差异的根本，即成岩作用、成岩矿物、成岩环境和成岩演化序列等在测井上的响应特征。

成岩作用类型和强度以及成岩矿物含量的差异通常在三孔隙度、自然伽马、自然电位、电阻率曲线上具有不同的响应特征，可通过岩心资料刻度测井定性确定储层成岩作用类型和强度以及计算成岩矿物的含量。

由于成岩矿物对成岩环境具有一定的指示意义，如伊蒙混层的出现一般指示埋藏较深的碱性环境，高岭石指示酸性成岩环境，钙质胶结一般存在于碱性环境，在酸性条件下将发生部分溶解，而长石、岩屑的溶解也存在于酸性环境，黄铁矿的存在一般指示还原环境，因此可通过成岩矿物及其组合交代关系的指示性作用确定成岩相成岩环境变化。

由于一个地区某个层段沉积物所处成岩阶段一般差异不太大，一般可根据镜下观察到的成岩特征，包括自生矿物分布形成顺序、黏土矿物组合、岩石的结构构造特点、有机质成熟度和古温度特征确定不同成岩相的成岩阶段(SY/T5477—2003)。与成岩作用和成岩阶段相关的成岩演化序列也可以通过薄片镜下观察，根据自生矿物生成顺序以及矿物穿插关系来确定[46]。

因此，成岩相测井识别关键在于测井资料中所包含的成岩相有关的信息的提取，其核心的技术就是通过测井资料确定储层经历成岩作用和强度以及精确计算指示性成岩矿物含量。本文即以碎屑岩主要的成岩作用(压实、压溶、胶结和溶解作用)为主线并兼顾成岩矿物等对不同成岩相的识别标志以及测井响应特征进行深入探讨与详细总结。

4.1 压实压溶相

识别标志：颗粒之间以点−线、线接触为主，塑性颗粒如云母的变形，脆性颗粒如石英等的破裂。压溶作用鉴别最典型的特征是颗粒之间呈缝合线接触。测井响应特征：典型压实压溶相的测井响应特征是体积密度大、声波时差小。强压实相发育层段往往泥质以及塑性矿物含量高、或者岩性以泥质粉砂岩、粉砂质泥岩等为主，因此，压实压溶相另一重要的鉴定特征就是补偿中子以及自然伽马值均较高。

4.2 胶结相

识别标志：胶结作用类型主要有钙质胶结、硅质胶结、泥质(黏土矿物)胶结、铁质胶结(黄铁矿等)、沸石以及硫酸盐胶结，其中以前三者为主。镜下各胶结相的识别主要就是鉴定胶结物的类型以及含量。由于长石类胶结物、硫酸盐胶结物(主要包括石膏和硬石膏)和铁质胶结物在各类碎屑岩中含量一般都非常少，对储集物性总体上影响不大，且剖面上和平面上的分布规模也不广，本文主要考虑钙质、硅质以及泥质胶结相。

钙质胶结物类型有方解石、白云石、文石等，在成岩早晚期均可以形成，早期的多以粒状、镶嵌状或栉壳状的结构出现，晚期一般形成较自形的铁方解石和铁白云石。碳酸盐胶结物强胶结时还可以使碎屑岩形成“悬浮砂”构造，即连晶胶结的方解石或白云石导致碎屑颗粒呈点接触、漂浮状，掩盖了强压实作用下颗粒线接触特征。测井上钙质胶结相具有高电阻的特征，即“钙尖峰”，一般声波时差小，体积密度大，自然伽马值很低。

硅质胶结物最常见的是石英的次生加大和孔隙中充填的自生石英，也有蛋白石、玉髓和燧石。测井曲线上表现为中−高密度，中−低声波时差，中−低电阻率，自然伽马值低，深浅侧向曲线为正幅度差[42]。

黏土矿物胶结物类型主要有高岭石、绿泥石、蒙脱石、伊利石以及伊蒙混层、绿蒙混层等，这里指的主要就是成岩过程中孔隙中沉淀生成或者再生的类型。黏土矿物胶结相最典型的特征的自然伽马值相对高值，中子孔隙度大、密度低、声波时差中-高，电阻率低。胶结作用总体是对储层储集物性起破坏性的成岩作用，然而，需要注意的是早期形成的环边状绿泥石胶结物由于一方面能够增加颗粒的抗压实能力抵御压实作用的影响，另一方面还可以阻止石英次生加大的进行，因此，绿泥石环边胶结是对储层储集物性起建设性的成岩作用[48]。绿泥石环边胶结相一般中子孔隙度相对较低，且中子−密度孔隙度值差异不大，根据这一点可通过三孔隙度曲线的差异将绿泥石环边胶结相与其它胶结成岩相区分开来。

4.3 溶蚀相

识别标志：溶蚀作用的结果是形成了砂岩中的次生孔隙，因此溶蚀相的识别主要是依据粒内溶孔、粒间溶孔、粒缘溶孔和铸模孔等。酸性条件下易溶蚀的组分主要有长石、岩屑、沸石和碳酸盐胶结物，在碱性成岩条件下，石英也可以发生一定的溶蚀。溶蚀的机理一是有机质生烃脱羧产生的有机酸和CO2，二是同生成岩期或者由于构造抬升沉积物至地表，从而接受的大气淡水淋滤。

测井曲线上溶蚀相具中等自然伽马、低密度，中等中子孔隙度和低声波时差的特征[25−26]。且一般情况下中子和密度测井反映储层总孔隙度，而声波时差则反映基质孔隙度，不反应次生孔隙，因此溶蚀相较发育的层段其中子−密度孔隙度往往比声波孔隙度高[49]，且深浅侧向组合正幅度差大[50]。

4.4 破碎裂缝相

识别标志：镜下可见粒缘缝等微裂缝，岩心上可见宏观裂缝。

测井响应特征：只利用常规测井资料对破碎裂缝(特别是微裂缝)成岩相储层的识别比较难以准确把握。声波时差曲线的周波跳跃现象，双侧向曲线呈正幅度差(指示高角度裂缝)和负幅度差(指示低角度裂缝)能一定程度上识别和探测裂缝[51]。在一些非常规测井如裂缝识别测井(FIL)上裂缝发育层段明显低阻异常，岩性密度测井(LDT)在光电截面曲线出现光电吸收截面指数高值，电成像测井上表现为暗色正弦曲线[25]。

基于上述构造相、沉积相对成岩相的控制作用的探讨，精细分析成岩相与构造相、沉积相和测井相之间的关系[25]。并根据不同成岩相储层的构造相、沉积环境和岩性特征的差异，以及成岩作用过程导致储层物性和孔隙结构特征的不同，找到不同成岩相的测井响应特征和识别标志，即可以达到利用常规测井资料识别成岩相的目的[26]。

5 测井识别成岩相关键问题探讨

测井识别成岩相首先应该重视基础地质研究，即在前期工作上把好关，包括：(1) 构造特征、沉积地层背景的研究和岩性的识别，只有在岩性正确识别的基础上，才能对岩石的胶结程度、压实程度和次生孔缝发育程度作出正确的解释和定量识别[2]。而成岩相与构造和沉积相是息息相关的，一定的成岩相总是与沉积相相联系的，因此，成岩相测井识别过程中应把握成岩相与构造-沉积相的相关关系，如山前褶皱带由于压实作用强，易发育强压实相，形成于水动力较弱沉积环境的层段，泥质含量高，抗压实能力差，也易发育强压实相。在地质历史时期沉积物被抬升至近地表接受大气淡水淋滤层段易发育溶蚀相，接近生烃洼陷中心层段由于易于有机酸接触也易发育溶蚀相等。(2) 镜下成岩特征的精细判别与鉴定，只有通过取心薄片的镜下观察，才能确定沉积物经历的成岩作用及其演化序列、成岩矿物类型、所处成岩环境和成岩阶段，在此基础上，才可以划分出不同的成岩相类型。研究时应充分利用有限的取心薄片资料，除较常用的铸体薄片、扫描电镜之外，更应充分利用地质分析测试新技术，如全岩矿物分析可用于分析储集岩矿物组合特征，而流体包裹体通过分析流体性质的变化为成岩环境变化分析提供有力的证据，拉曼光谱更能准确识别矿物的组合特征，反映成岩环境变化[18]。(3)测井资料的多解性要求我们建立不同成岩相测井识别模式时应对多口井测井资料进行综合对比分析。

在研究区构造、沉积背确定以及镜下成岩特征观察鉴定的前提下，可根据深度归位后的岩心薄片资料刻度测井原则，精细厘定不同成岩相在不同测井曲线及其组合上的响应特征，并通过岩心观察结果和岩电等资料的验证、修正和标定，从而建立测井信息和地质信息之间的正演、反演模型，提高测井识别成岩相解释精度。如此便可以建立适合特定地区特定层段的成岩相测井识别标志和测井识别方法序列，实现成岩相的测井识别[2]。

前已述及，成岩作用类型和强度以及成岩矿物的差异通常在测井曲线上具有不同的响应特征，当然复杂岩性储层由于其矿物成分复杂且随深度变化大，造成经历不同成岩作用类型和强度的沉积物测井响应特征相同或者近似，或者经历相同成岩作用的层段由于成岩矿物成分及含量不同造成其测井响应特征的明显差异，这二者将增加测井识别成岩相的难度。因此，在利用测井资料识别成岩相时，应在岩心资料刻度测井基础上确定成岩作用类型、强度的和计算成岩矿物的含量，并通过岩电资料的修正和标定来确定不同层段成岩相类型。

6 成岩相测井识别方法进展

6.1 测井识别成岩相主要方法

一般来说，对成岩相较敏感的常规测井曲线主要是密度、中子、声波时差、自然伽马、电阻率和自然电位等，在实际的操作过程中，应结合地区的地质特征，优选对成岩相灵敏性较好的曲线组合。不同的学者提出了不同成岩相的测井识别方法，有交会图法[26,34,42,52−54]，也有将成岩相跟岩电、压汞实验参数相对应建立判别公式的方法[43]。段新国等[55]通过成岩作用类型确定、归纳统计不同成岩相特征矿物含量、孔隙类型及孔隙度，并参考了各井电测解释的储层孔渗情况，绘制了四川盆地须二段储层成岩相平面分布图。

本文作者认为成岩相应是沉积物经历种种成岩作用后特征的综合，应通过对应实际取芯段岩电、孔渗资料标定。为了使成岩相直观地表达出来，可用蜘蛛网图(或梯形图)来表示成岩相，即以能够反映相特征的各种参数值作为辐射轴，以不同成岩相之间的差别为依据，以图形区分不同成岩相，图1所示为姬塬地区长8成岩相测井响应蜘蛛图，是根据石玉江等[26]总结的5种不同成岩相在自然伽马、自然电位、密度、补偿中子、声波时差和中子−密度孔隙度差 6条曲线上的响应特征投点形成的识别成岩相的蜘蛛网图，具有简洁明了的特点。

对于成岩矿物含量的计算，则可利用测井资料中的自然电位、自然伽玛等参数，或通过孔隙度系列中子-密度交会图技术确定碎屑岩中泥质含量。在确定泥质含量基础上，充分应用组合测井资料及岩性−密度测井、自然伽玛能谱测井等资料与黏土矿物关系，应用交会图技术进行黏土矿物分析，或进行黏土矿物含量与自然伽马能谱测井(U，Th，K和TH/K)的多元线性回归，得到回归方程，确定出黏土矿物类型及其含量[55−56]。

图1 姬塬地区长8成岩相测井响应蜘蛛图Fig. 1 Spider web diagram of logging response characteristics of diagenetic facies of Chang 8 reservoirs in Jiyuan regiom

综合利用储层中不同矿物成份在不同曲线上的响应特征求取各矿物的含量，一般要采用多参数优化技术而非单一方程求解地层组分，即控制不同输入曲线的权值等参数，这样对测井信息的利用率高，计算精度也最大[57]。即要使得输入的测井曲线条数大于等于要求解的矿物组分数，即保证方程组为超定线性方程组，如下式：

式中：m为测井曲线条数；n为待计算储集岩矿物组分的个数；Xi为第i种矿物组分的相对含量；Aij为第i种矿物组分对第j种曲线的响应值；Bj为第j种测井曲线的响应值。

以四川盆地中部须家河组须二段低孔低渗复杂碎屑岩储层为例，须二段储层经历压实压溶、胶结和溶蚀等成岩作用，胶结物类型主要为碳酸盐和黏土矿物，黏土矿物尤以伊利石最为发育，少量为伊蒙混层。在对须二段储层成岩作用、成岩矿物、成岩阶段和成岩演化序列等研究的基础上，按成岩作用类型和强度、成岩矿物及其对储层储集物性的影响，将须二段储层成岩相划分为强压实胶结相、强压实次生溶蚀相和破碎裂缝相3种成岩相。

通过输入AC，RT，GR，DEN和CNL等曲线，借助FORWARD平台的多矿物模型分析，计算得到了砂质、泥质和方解石等不同岩石、矿物成分类型和其相应含量如图2所示。再通过自然伽马能谱测井计算出了不同层段的伊利石和伊蒙混层黏土矿物的相对含量，由此识别出了强压实胶结成岩相。根据溶蚀成岩相的测井响应机理，利用声波孔隙度和密度−中子孔隙度二者之差，定性判别出了溶蚀作用发育层段，再剔除方解石胶结层段以及黏土矿物胶结层段，可确定出强压实次生溶蚀相发育层段(图2)。破碎裂缝相则主要通过 EMI成像测井图上的裂缝发育层段拾取来确定。该区各单井纵向上的成岩相的测井识别为后期工作如优质储集体的预测提供了先决条件。

6.2 成岩相测井识别未来发展方向

成岩相测井识别是测井地质学与沉积成岩学的有机渗透与融合，成岩相测井识别方法的进展一是体现在测井技术的发展，如除常规测井之外的成像测井、元素俘获测井等非常规测井技术的发展，二是测井资料解释处理新方法的进展，挖掘隐藏在丰富的测井资料中的关于成岩相相关的地质信息也是成岩相测井识别方法未来发展方向。

由于成岩环境代表了一定的沉积环境，成岩相与沉积相是紧密关联、相互渗透的。目前通过测井资料识别沉积相的方法、理论与技术已经比较成熟，并已形成一门新的边缘学科—测井沉积学[58]。利用测井资料识别成岩相时也可以借鉴测井沉积学中的理论和方法。如测井沉积相分析方法中的多元统计、神经网络模式识别技术等与岩石组合(成分)粒度相关的岩性相分析程序[27]在成岩相的测井识别中也具有较好的应用前景。

图2 蓬莱地区须二段储层成岩相测井识别(PORA为声波孔隙度，PORT为密度和中子孔隙度平均值，1为气水同层)Fig. 2 Logging recognition of reservoir diagenetic facies of Xujiahe Formation Member 2 in Penglai Area

自然伽玛能谱、岩性密度测井、激发伽玛能谱测井等地球化学测井技术可直接测量岩石中Al，Ca，Si，Fe，S，Cl和H等12种元素成分，并经过计算和处理可用其合成的6种矿物，几乎能测出岩石骨架的全部成分，并利用其与岩性之间的相关关系可推算成岩作用并精确计算部分成岩矿物含量。而光电吸收指数、多元回归分析及非线性理论等方法在研究矿物成分中的广泛应用，使得岩石成分的识别能力得到提高[59]，再通过实验室岩心薄片分析结果的验证和标定，结合反应胶结程度、次生孔缝发育程度的孔隙度测井序列，即可以提高成岩相测井识别的精度。

总之，虽然成岩相测井识别这一领域相对起步较晚、研究基础薄弱，要实现这一方法的广泛应用尚需长时间探索，但可以预料随着时间的推移，应用需求的日趋广泛和测井新技术的不断发展，测井资料在储层成岩相研究中的应用将更加广阔，而其解释的精度也将不断增强。

[1]邹才能, 陶士振, 周慧, 等．成岩相的形成、分类与定量评价方法[J]. 石油勘探与开发, 2008, 35(5): 526−540.ZOU Caineng, TAO Shizhen, ZHOU Hui, et al. Genesis,classification and evaluation method of diagenetic facies[J].Petroleum Exploration and Debelopment, 2008, 35(5): 526−540.

[2]陈彦华, 刘莺. 成岩相: 储集体预测的新途径[J]. 石油实验地质, 1994, 16(3): 274−280.CHEN Yanhua, LIU Ying. Diagenetic facies: A new approach to the prediction of reservoir rocks[J]. Experimental Petroleum Geology, 1994, 16(3): 274−280.

[3]Railsback L B. Carbonate diagenetic facies in the upper pennsylvanian dennis formation in iowa, missouri, and kansas[J].Journal of Sedimentary Petrology, 1984, 54(3): 986−999.

[4]杜业波, 季汉成, 吴因业, 等. 前陆层序致密储层的单因素成岩相分析[J]. 石油学报, 2006, 27(2): 48−52.DU Yebo, JI Hancheng, WU Yinye, et al. Single factor diagenetic facies analysis of tight reservoir in western Sichuan foreland basin[J]. Acta Petroleum Sinica, 2006, 27(2): 48−52.

[5]庞雄奇, 李丕龙, 陈冬霞, 等. 陆相断陷盆地相控油气特征及其基本模式[J]. 古地理学报, 2011, 13(2): 55−73.PANG Xiongqi, LI Peilong, CHEN Dongxia, et al.Characteristics and basicmodel of facies controlling oil and gas in continental fault basin[J]. Journal of Palaeogeography, 2011,13(2): 55−73.

[6]庞雄奇, 李丕龙, 张善文, 等. 陆相断陷盆地相-势耦合控藏作用及其基本模式[J]. 石油与天然气地质, 2007, 28(5):641−651.PANG Xiongqi, LI Pilong, ZHANG Shanwen, et al. Control of facies potential coupling on hydrocarbon accumulation in continental faulted basins and its basic geologicalmodels[J]. Oil and Gas Geology, 2007, 28(5): 641−651.

[7]孟元林, 高建军, 刘德来, 等. 渤海湾盆地西部凹陷南段成岩相分析与优质储层预测[J]. 沉积学报, 2006, 24(2): 185−189.MENG Yuanlin, GAO Jianjun, LIU Delai, et al. Diagenetic fac ies analysis and high quality reservoir prediction in the southern Xibu depression of the boha iwan basin[J]. Acta Sedimentologica Sinica, 2006, 24(2): 185−189.

[8]孟元林, 李娜, 黄文彪, 等. 辽河坳陷西部斜坡带南段新生界成岩相分析与优质储集层预测[J]. 古地理学报, 2008, 10(1):33−40.MENG Yuanlin, LI Na, HUANG Wenbiao, et al. Diagenetic facies analysis and high quality reservoir prediction of the Cenozoic in southern segment ofwestern slope of Liaohe Depression[J]. Journal of Palaeogeography, 2008, 10(1): 33−40.

[9]刘清俊, 于炳松, 周芳芳, 等. 阿克库勒凸起东河砂岩成岩作用与成岩相[J]. 西南石油大学学报: 自然科学版, 2010, 33(5):54−63.LIU Qingjun, YU Bingsong, ZHOU Fenfang, et al. Diagenesis and diagenetic facies of donghe sandstone of akekule uplift in tarim basin[J]. Journal of Southwest Petroleum University:Science ＆ Technology Edition, 2011, 33(5): 54−63.

[10]Rudolf B, Chagnon Ré, Malo M. Sedimentologic, diagenetic and tectonic evolution of the saint-flavien gas reservoir at the structural front of the quebec appalachians[J]. Bulletin of Canadian Petroleum Geology, 2003, 51(2): 126−154.

[11]李晓光, 陈振岩. 辽河坳陷古近系碎屑岩储层孔隙演化特征探讨[J]. 古地理学报, 2006, 8(2): 252−259.LI Xiaoguang, CHEN Zhenyan. Discussion on pore evolution features of clastic rock reservoir of the paleogene in liaohe depression[J]. Journal of Palaeogeography, 2006, 8(2): 252−259.

[12]王琪, 禚喜准, 陈国俊, 等. 鄂尔多斯西部长6砂岩成岩演化与优质储层[J]. 石油学报, 2005, 26(5): 17−25.WANG Qi, ZHUO Xizhun, CHEN Guojun, et al. Diagenetic evolution and high-quality reservoir in Chang6 Sandstone in the western Ordos basin[J]. Acta petroleum Sinica, 2005, 26(5):17−25.

[13]陈振岩, 孟元林, 高建军, 等. 鸳鸯沟洼陷西斜坡成岩作用定量表征及有利区带预测[J]. 地学前缘, 2008, 15(1): 71−76.CHEN Zhenyan, MENG Yuanlin, GAO Jianjun, et al.Quantitative characterization of diagenesis and prediction of favorable areas in the west slope of Yuanyanggou subsag[J].Earth Science Frontiers, 2008, 15(1): 71−79.

[14]Sylvia M C, Anjos D, de Ros L F. Depositional and diagenetic controls on the reservoir quality of lower cretaceous pendeˆncia sandstones, potiguar rift basin, brazil[J]. AAPG Bulletin, 2000,84(11): 1719−1742.

[15]Maas T E, Jens J, Knut B. Diagenetic controls on reservoir quality in middle to upper jurassic sandstones in the south viking graben, north sea[J]. AAPG Bulletin, 2011, 95(11): 1937−1958.

[16]Roger J B, Oetting G C, GAO Guoqiu. Strontium isotopic signatures of oil-field waters: Applications for reservoir characterization[J]. AAPG Bulletin, 2004, 88(12): 1677−1704.

[17]Aysen O, Cumella S P, Milliken K L, et al. Prediction of lithofacies and reservoir quality using well logs, late cretaceous williams fork formation, mamm creek field, piceance basin,colorado[J]. AAPG Bulletin, 2011, 95(10): 1699−1723.

[18]周瑶琪, 周振柱, 陈勇, 等. 东营凹陷民丰地区深部储层成岩环境变化研究[J]. 地学前缘, 2011, 18(2): 268−276.ZHOU Yaoqi, ZHOU Zhenzhu, CHEN Yong, et al. Research on diagenetic environmental changes of deep reservoir in Minfeng Area, Dongying Sag[J]. Earth Science Frontiers, 2011, 18(2):268−276.

[19]Doherty P D, Soreghan G S, Castagna J P. Outcrop-based reservoir characterization: A composite phylloid-algal mound,western orogrande basin(New Mexico)[J]. AAPG Bulletin, 2002,86(5): 779−795.

[20]Nakazawa T, Ueno K. Sequence boundary and related sedimentary and diagenetic facies formed on middle permian mid-oceanic carbonate platform: Core observation of Akiyoshi Limestone, Southwest Japan[J]. Facies, 2004, 50: 301−311.

[21]Kalefa M A, Ghali E, Howri M, et al. Distribution of diagenetic alterations in glaciogenic sandstones within a depositional facies and sequence stratigraphic framework: Evidence from the upper ordovician of the Murzuq Basin, SW Libya[J]. Sedimentary Geology, 2006, 190: 323−351.

[22]Kalefa M A, Ghali El, Sadoon M, et al. Distribution of diagenetic alterations within depositional facies and sequence stratigraphic framework of fluvial sandstones: Evidence from the petrohan terrigenous group, lower triassic, NW Bulgaria[J].Marine and Petroleum Geology, 2009, 26: 1212−1227.

[23]Guidry S A, Chafetz H S. Depositional facies and diagenetic alteration in a relict siliceous hot-spring accumulation: Examples from Yellowstone National Park, USA[J]. Journal of Sedimentary Research, 2003, 73(5): 806−823.

[24]赖锦, 王贵文, 郑懿琼, 等. 金秋区块须四段储层成岩相及测井识别[J]. 西南石油大学学报: 自然科学版, 2013, 35(5):41−49.LAI Jin, WANG Guiwen, ZHENG Yiqiong, et al. Diagenetic facies of Xujiahe Formation Member 4 reservoir in Jinqiu area and the logging recognition[J]. Journal of Southwest Petroleum University: Science ＆ Technology Edition, 2013, 35(5): 41−49.

[25]杨威, 魏国齐, 金惠, 等. 碳酸盐岩成岩相研究方法及其应用—以扬子地块北缘飞仙关组鲕滩储层为例[J]. 岩石学报,2011, 27(3): 749−756.YANG Wei, WEI Guoqi, JIN Hui, et al. Methods for diagenetic facies research on carbonate rocks and its application: Example from oolitic beach reservoir of the Feixianguan Formation in The Northern Margin of Yangtze Block[J]. Acta Petrologica Sinica,2011, 27(3): 749−756.

[26]石玉江, 肖亮, 毛志强, 等. 低渗透砂岩储层成岩相测井识别方法及其地质意义: 以鄂尔多斯盆地姬塬地区长 8段储层为例[J]. 石油学报, 2011, 32(5): 820−827.SHI Yujiang, XIAO Liang, MAO Zhiqiang, et al. An identification method for diagenetic facies with well logs and its geological significance in low-permeability sandstones: A case study on Chang8 reservoirs in the Jiyuan region, Ordos basin[J].Acta Petroleum Sinica, 2011, 32(5): 820−827.

[27]王贵文, 郭荣坤. 测井地质学[M]. 北京: 石油工业出版社,2000: 1−9.WANG Guiwen, GUO Rongkun. Well logging geology[M].Beijing: Petroleum Industry Press, 2000: 1−9.

[28]邹才能, 陶士振, 薛叔浩. “相控论”的内涵及其勘探意义[J].石油勘探与开发, 2005, 32(6): 7−12.ZOU Caineng, TAO Shizhen, XUE Shuhao. Connotation of“facies control theory” and its significance for exploration[J].Petroleum Exploration and Development, 2005, 32(6): 7−12.

[29]谭开俊, 许多年, 尹路, 等. 准噶尔盆地乌夏地区三叠系成岩相定量研究[J]. 岩性油气藏, 2011, 23(6): 24−28.TAN Kaijun, XU Duonian, YIN Lu, et al. Quantitative study on Triassic diagenetic facies in Wuxia area, Junggar Basin[J].Lithologic Resercoirs, 2011, 23(6): 24−28.

[30]王留奇, 赵微林, 刘孟慧, 等. 油区构造沉积成岩相带研究[J].石油勘探与开发, 1992, 19(5): 26−30.WANG Liuqi, ZHAO Weilin, LIU Menghui, et al. A study on the structural sedimentary diagenetic zone in an oil field[J].Petroleum Exploration and Development, 1992, 19(5): 26−31.

[31]王岚, 邹才能, 林潼, 等. 鄂尔多斯盆地白豹−华池地区延长组长 6油层组沉积环境及成岩相分析[J]. 天然气地球科学,2011, 22(5): 796−805.WANG Lan, ZOU Caineng, LIN Tong, et al. Sedimentary environment and diagenesis facies of Chang6 Yanchang Formation in Baibao-Huachi, Ordos Basin[J]. Natural Gas Geosciences, 2011, 22(5): 796−805.

[32]张明松, 雷卞军, 黄有根, 等. 鄂尔多斯盆地余兴庄−子洲地区上古生界山 32储层砂岩成岩作用与成岩相[J]. 沉积学报,2011, 29(6): 1031−1040.ZHANG Mingsong, LEI Bianjun, HUANG Yougen, et al.Diagenesis and diagenetic facies of sandstones reservoir in the upper palaeozoic Shan32of Yuxingzhuang-Zizhou Area, Ordos Basin[J]. Acta Sedimentologica Sinica, 2011, 29(6): 1031−1040.

[33]赖锦, 王贵文, 王书南, 等. 碎屑岩储层成岩相研究现状及进展[J]. 地球科学进展, 2013, 28(1): 39−50.LAI Jin, WANG Guiwen, WANG Shu’nan, et al. Research status and advances in the diagenetic facies of clastic reservoirs[J]. Advances in Earth Science, 2013, 28(1): 39−50.

[34]张响响, 邹才能, 陶士振, 等. 四川盆地广安地区上三叠统须家河组四段低孔渗砂岩成岩相类型划分及半定量评价[J]. 沉积学报, 2010, 28(1): 50−57.ZHANG Xiangxiang, ZOU Caineng, TAO Shizhen, et al.Diagenetic facies types and semiquantitative evaluation of low porosity and permeability sandstones of the Fourth Member Xujiahe Formation Guangan Area, Sichuan Basin[J]. Acta Sedimentologica Sinica, 2010, 28(1): 50−57.

[35]钟广法, 邬宁芬. 成岩岩相分析及其在泌阳凹陷核三下亚段砂岩储层中的初步应用[J]. 江汉石油学院学报, 1995, 17(1):37−40.ZHONG Guangfa, WU Ningfen. Diagenetic lithofacies analysis and its application in sandstone reservoirs of Biyang depression[J]. Journal of Jianghan Petroleum Institute,1995,17(1): 37−40.

[36]Aleta A. Mineralogical descriptions of the bentonite in balamban,cebu province, philippines[J]. Clay Science, 2000, 11(3):299−316.

[37]Grigsby J D, Langsford R P. Effects of diagenesis on enhanced-resolution bulk density logs in tertiary gulf coast sandstones: An example from the Lower Vicksburg Formation,Mcalien. Ranch Field, South Texas[J]. AAPG Bulletin, 1996,80(11):1801−1819.

[38]Elfigih O B. Regional diagenesis and its relation to facies Change in the Upper Silurian, Lower Acacus Formation,Hamada (Ghadames)Basin, Northwestern Libya[J]. PHD Thesis,1999, 62(10): 4426−4427.

[39]Abercrombie H J, Hutcheon I E, Bloch J D, et a1. Silica activity and the smectite-illite reaction[J]. Geology, 1994, 22(6):539−542.

[40]Turner J.R. Recognition of low resistivity, high permeability reservoir beds in the travis peak and cotton valley of east texas[J].Gulf Coast Association of Geological Societies Transactions, 1997, 4: 585−593.

[41]张胜斌, 王琪, 李小燕, 等.川中南河包场须家河组砂岩沉积−成岩作用[J].石油学报, 2009, 30(2): 225−232.ZHANG Shengbin, WANG Qi, LI Xiaoyan, et al. Depositional diagenetic coupling complex of Xujiahe sandstone in Hebaochang Block in the south part of the Central Sichuan Basin[J]. Acta Petroleum Sinica, 2009, 30(2): 225−230.

[42]张响响, 邹才能, 朱如凯, 等. 川中地区上三叠统须家河组储层成岩相[J]. 石油学报, 2011, 32(2): 257−264.ZHANG Xiangxiang, ZOU Caineng, ZHU Rukai, et al.Reservoir diagenetic facies of the upper Triassic Xujiahe formation in the central Sichuan basin[J]. Acta Petroleum Sinica,2011, 32(2): 257−264.

[43]时卓, 张海涛, 王永莉. 苏里格气田储层成岩相类型及测井识别[J]. 国外测井技术, 2011, 22(4): 13−18.SHI Zhuo, ZHANG Haitao, WANG Yongli. Diagenetic facies type of the reservoirs in sulige gas field and logging identification[J]. World Well Logging Technology, 2011, 22(4):13−18.

[44]张顺存, 刘振宇, 刘巍, 等. 准噶尔盆地西北缘克－百断裂下盘二叠系砂砾岩储层成岩相研究[J]. 岩性油气藏, 2010, 22(4):43−51.ZHANG Shuncun, LIU Zhenyu, LIU Wei, et al. Diagenesis facies of Permian sandy conglomerate reservoir in footwall of Kebai Fault in northwestern margin of Junggar Basin[J].Lithologic Reservoirs, 2010, 22(4): 43−51.

[45]张琴, 朱筱敏, 陈祥, 等. 南华北盆地谭庄凹陷下白垩统成岩相分布及优质储层预测[J]. 石油与天然气地质, 2010, 31(4):472−481.ZHANG Qin, ZHU Xiaomin, CHEN Xiang, et al. Distribution of diagenetic facies and prediction of high-quality reservoirs in the Lower Cretaceous of the Tanzhuang Sag, the Southern North China Basin[J]. Oil and Gas Geology, 2010, 31(4): 472−481.

[46]郑荣才, 耿威, 周刚, 等. 鄂尔多斯盆地白豹地区长6砂岩成岩作用与成岩相研究[J]. 岩性油气藏, 2007, 19(2): 1−7.ZHENG Rongcai, GENG Wei, ZHOU Gang, et al. Diagenesis and diagenetic facies of Chang6 sandstone of Yanchang Formation in Baibao Area, Ordos Basin[J]. Lithologic Reservoirs, 2007, 19(2): 1−7.

[47]何周, 史基安, 唐勇, 等. 准噶尔盆地西北缘二叠系碎屑岩储层成岩相与成岩演化研究[J]. 沉积学报, 2011, 29(6):1069−1077.HE Zhou, SHI Jian, TANG Yong, et al. Characteristics of diagenesis and diagenetic facies of permian clastic reservoir in Northwest Margin of Junggar Basin[J]. Acta Sedimentologica Sinica, 2011, 29(6): 1069−1077.

[48]刘金库, 彭军, 刘建军, 等. 绿泥石环边胶结物对致密砂岩孔隙的保存机制: 以川中-川南过渡带包界地区须家河组储层为例[J]. 石油与天然气地质, 2009, 30(1): 53−58.LIU Jinku, PENG Jun, LIU Jianjun, et al. Pore preserving mechanism of chlorite rims in tight sandstone: An example from the T3x Formation of Baojie Area in the transitional zone from the central to Southern Sichuan Basin[J]. Oil and Gas Geology,2009, 30(1): 53−58.

[49]唐文生, 赵德勇, 刘兴礼. 塔中地区生物碎屑灰岩储层分析与评价[J]. 测井技术, 1996, 20(4): 250−255.TANG Wensheng, ZHAO Deyong, LIU Xingli. Analysis and evaluation of bioclastic limestone reservoir in the central region of Tarim Basin[J]. Well Logging Technology, 1996, 20(4):250−255.

[50]林文姬, 汤达祯, 徐凤银, 等. 苏里格气田盒8段成岩相类型及其测井标志[J]. 石油天然气学报, 2010, 32(2): 271−274.LIN Wenji, TANG Dazhen, XU Fengyin, et al. Diagenesis facies type of H8 reservoir in Sulige Gas Field and its logging identification marker[J]. Journal of Oil and Gas Technology,2010, 32(2): 271−274.

[51]陈莹, 谭茂金. 利用测井技术识别和探测裂缝[J]. 测井技术,2003, 27(增刊): 11−14.CHEN Ying, TAN Maojin. Fracture detection and identification using logging techniques[J]. Well Logging Technology, 2003,27(supplement): 11−14.

[52]孙玉善, 申银民, 徐迅, 等. 应用成岩岩相分析法评价和预测非均质性储层及其含油性: 以塔里木盆地哈得逊地区为例[J].沉积学报, 2002, 20(1): 55−60.SUN Yushan, SHEN Yinmin, XU Xun, et al. Evaluating and predicting heterogeneous reservoirs and its oil-bearing properties by the analysis technique of the diagenetic lithofacies: Taking Hadexun Area in Tarim Basin as an example[J]. Acta Sedimentologica Sinica, 2002, 20(1): 55−60.

[53]孙玉善, 申银民, 徐迅, 等. 应用成岩岩相分析法评价和预测非均质性储层及其含油性: 以塔里木盆地哈得逊地区为例[J].沉积学报, 2002, 20(1): 55−60.SUN Yushan, SHEN Yinmin, XU Xun, et al. Evaluating and predicting heterogeneous reservoirs and its oil-bearing properties by the analysis technique of the diagenetic lithofacies: Taking Hadexun Area in Tarim Basin as an example[J]. Acta Sedimentologica Sinica, 2002, 20(1): 55−60.

[54]吴少波, 李亮, 余小雷, 等. 靖安油田大路沟区长6油层组储层成岩作用与成岩相[J]. 西安石油大学学报: 自然科学版,2009, 24(1): 8−13.WU Shaobo, LI Liang, YU Xiaolei, et al. Diagenesis and diagenetic facies of the Upper Triassic Chang26 reservoir in Dalugou Area, Jing’an Oilfield[J]. Journal of Xi’an Shiyou University: Natural Science Edition, 2009, 24(1): 8−13.

[55]段新国, 宋荣彩, 李国辉, 等. 四川盆地须二段综合成岩相特征研究[J]. 西南石油大学学报: 自然科学版, 2011, 33(1):7−14.DUAN Xinguo, SONG Rongcai, LI Guohui, et al. Research of integrated diagenetic facies characteristics of t3x2 reservoir in Sichuan Basin[J]. Journal of Southwest Petroleum University:Science ＆ Technology Edition, 2011, 33(1): 7−14.

[56]张小莉, 沈英, 陈文学, 等. 利用测井资料分析成岩作用对储集层的影响[J]. 沉积学报, 2000, 18(1): 127−131.ZHANG Xiaoli, SHEN Ying, CHEN Wenxue, et al. Application of well logging information to analysis on the effect of diagenesis in reservoirs[J]. Acta Sedimentologica Sinica, 2000,18(1): 127−131.

[57]刘菁华, 王祝文, 易清平, 等. 利用自然γ能谱测井资料确定黏土矿物的含量及其应用[J]. 吉林大学学报: 地球科学版,2010, 40(1): 215−222.LIU Jinghua, WANG Zhuwen, YI Qingping, et al. Determinating clay mineral content by natural Gamma-ray spectral logging data[J]. Journal of Jilin University: Earth Science Edition, 2010,40(1): 215−222.

[58]李高仁, 郭清娅, 石玉江, 等. 鄂尔多斯盆地高自然伽马储层识别研究[J]. 测井技术, 2006, 30(6): 511−515.LI Gaoren, GUO Qingya, SHI Yujiang, et al. Identification of high Gamma ray reservoir in Ordos Basin[J]. Well Logging Technology, 2006, 30(6): 511−515.

[59]尹寿鹏, 王贵文. 测井沉积学研究综述[J]. 地球科学进展,1999, 14(5): 440−445.YIN Shoupeng, WANG Guiwen. A summary of sedimentology from well logs[J]. Advance in Earth Sciences, 1999, 14(5):440−445.

[60]张兆辉, 高楚桥, 刘娟娟, 等. 基于地层组分分析的储层孔隙度计算方法研究[J]. 岩性油气藏, 2012, 24(1): 97−99.ZHANG Zhaohui, GAO Chuqiao, LIU Juanjuan, et al.Calculation method of porosity based on formation component analysis[J]. Lithologic Reservoirs, 2012, 24(1): 97−99.