川中地区下二叠统栖霞组白云岩储层特征和成因

2022-01-07段军茂郑剑锋沈安江姚倩颖

海相油气地质 2021年4期

段军茂，郑剑锋，沈安江，朱茂，姚倩颖，郝毅

1中国石油勘探开发研究院；2中国石油杭州地质研究院；3中国石油集团碳酸盐岩储层重点实验室

0前言

近年来，四川盆地下二叠统栖霞组不断获得勘探发现，被认为是四川盆地下一步勘探的重要接替领域［1］。四川盆地中部（简称川中）地区也相继在MX42井、MX31X1井、HS2井、HS4井等井钻遇白云岩储层并连获高产工业气流，揭示了栖霞组巨大的勘探潜力。前人主要对四川盆地栖霞组白云岩储层的成因开展了大量研究，主要观点包括准同生白云石化［2］、埋藏成因［3］、浅埋藏期的构造热液白云石化［4-5］等，以及多期白云石化叠合作用，如回流渗透白云石化与断控热液白云石化的叠合［6-7］，准同生白云石化、浅埋藏白云石化和热液白云石化的叠合［8-9］。也有不少学者开展了白云岩储层主控因素的相关研究，认为白云岩储层的发育受控于沉积相带［10-11］、构造-热液断裂［12］、层序界面控制的早期岩溶［13］、快速埋藏［3］等等，其中以受控于沉积相带为主流观点。然而这些观点却难以进一步解释白云岩储层单层厚度薄、横向非均质性强等现象，而这种现象在川中地区尤为明显。另一方面，前人针对栖霞组白云岩储层的研究多为盆地尺度或集中于川西地区，川中地区的相关研究相对较少。基于此，本文将川中地区作为研究对象，选取MX42井、MX108井等9口井，共取岩心样品402块，结合钻遇栖霞组的钻井测井资料，进行了岩心、薄片观察及物性分析，指出了川中栖霞组白云岩储层特征和分布规律；在此基础上，结合样品白云石有序度、稳定碳氧锶同位素、主量元素、微量元素、稀土元素和阴极发光等地球化学分析和碳酸盐岩激光U-Pb定年测试，对川中栖霞组白云岩成岩作用和储层成因展开了系统分析，并指出了川中地区栖霞组白云岩储层主控因素。所获得的认识不仅对川中地区，也对盆地内其他区域的栖霞组白云岩储层的勘探有所裨益。

1 区域地质背景

四川盆地在大地构造上位于扬子板块西北部，二叠纪处于赤道附近，构造稳定，为被动大陆边缘背景［14］，满足生物大量繁殖的湿热环境。志留纪末的加里东运动促使了在晚震旦世就已初具雏形的乐山—龙女寺古隆起的最终形成，并持续影响到二叠纪沉积［15］，从而控制了栖霞组西南高、东北低的沉积格局，也为盆地中部地区的台内滩和岩溶作用发育奠定了基础。

研究区位于四川盆地川中低缓褶皱带中东部（图1a）。研究区内下二叠统由下往上依次包括梁山组、栖霞组，往上为中二叠统茅口组［16］（图1b）。梁山组与下伏地层呈不整合接触，为一套滨岸—沼泽相砂泥岩夹煤线沉积；茅口组底部为一套富含泥质的生物碎屑灰岩和泥晶生物碎屑灰岩，与栖霞组呈Ⅱ型不整合接触。研究区内栖霞组厚度介于100～150 m，分为栖霞组一段（简称栖一段）和栖霞组二段（简称栖二段）。四川盆地从晚泥盆世到二叠纪空谷阶（即栖霞阶）为持续性的晚古生代大冰期［17］，栖霞阶处于冰室期向温室期过渡阶段，受冰川消融、构造活动和海底扩张等因素的影响，发生了2期三级海平面升降［18-19］，因此栖霞组可划分为两个三级层序SQ1和SQ2，并分别对应栖一段和栖二段，每一段又包括海侵体系域（TST）和高位体系域（HST）（图1b）。进入栖霞期后，盆地范围内海平面上升引发大规模海侵，使区内由滨岸砂泥岩沉积演变为碳酸盐台地沉积［20］，TST时期为开阔台地沉积，发育深灰色—黑灰色含泥泥晶灰岩与生物碎屑泥晶灰岩，其中SQ1的TST沉积整合于梁山组之上；HST时期海平面下降，川中地区以台内颗粒滩和滩间海沉积为主［21-22］，在海平面下降到一定程度后，易形成局限环境，为早期暴露溶蚀及白云石化作用奠定了基础。另一方面，二叠纪也是泛大陆形成的时期［17］，与之相伴生的中二叠世末期的峨眉山玄武岩喷发［23-24］，其形成的断裂系统和热液流体为二叠纪地层提供了热成岩背景。

2 储层地质特征

2.1 白云岩岩石学特征

川中地区下二叠统栖霞组岩性以生物碎屑灰岩和泥质生物碎屑灰岩为主，部分井段于HST时期发育储集能力较好的白云岩，按晶体大小可分为泥粉晶白云岩、粉细晶白云岩和中粗晶白云岩，另外在部分溶蚀孔洞和裂缝中可见鞍状白云石。研究区栖霞组储层主要发育于粉细晶白云岩和中粗晶白云岩（粉—粗晶白云岩）中。

晶粒白云岩宏观上呈现浅灰色—灰色，以中—薄层状为主。显微镜下白云石晶体按晶形可分为他形晶和半自形—自形晶。他形晶以非平直晶面为主，晶面较脏，晶体形态和轮廓多不清晰（图2a），部分需在正交光下方可识别颗粒形态,多具有波状消光的特点；偶见残余生物碎屑（图2b），部分他形晶白云岩经原岩恢复可以发现较为明显的生物碎屑结构（图2c，2d），包括有孔虫、类、棘屑、珊瑚等；阴极发光呈昏暗发光或不发光（图2e，2f）。半自形—自形晶以平直晶面为主，多见于中粗晶白云岩；晶粒核心较脏而环边颜色常较亮，晶形较好者可见较为明显的亮边雾心结构（图2g，2h）；阴极发光具有核心昏暗发光和环边明亮发光的特征（图2i）。

本刊讯 12月5日下午，省人大常委会党组理论学习中心组集体学习研讨会在济南举行，集体学习《习近平谈治国理政》第二卷和《习近平新时代中国特色社会主义思想三十讲》部分章节，并结合学习体会和工作实际进行研讨交流。省人大常委会党组书记、副主任于晓明主持会议并讲话；党组副书记、副主任王良，党组成员、副主任齐涛，党组成员、副主任王云鹏出席；副主任王随莲列席。

图2 川中地区栖霞组二段白云岩微观特征Fig.2 Microscopic characteristics of dolomite of of the Qixia Member 2 in central Sichuan Basin

研究区栖霞组除粉—粗晶白云岩外，还可见泥粉晶白云岩和鞍状白云石。泥粉晶白云岩在研究区较少见，目前仅见于HS2井中上部，泥粉晶颗粒呈砂屑状，并多为亮晶白云石环边包围（图2j），紧邻残余生物碎屑的粉细晶白云岩，孔隙被沥青充填或半充填。另外川中地区栖霞组于岩心和显微镜下均可见非平直晶面的鞍状白云石充填于溶洞或溶缝中，白云石晶体粗大，晶面弯曲程度较大，正交偏光下具波状消光等特点（图2k）。

此外，笔者建议将高校法学院法律诊所的学生吸纳到社区矫正工作中。相较于其他的志愿者，法律诊所的学生具备专业的法律知识背景，对社区矫正工作也有一定的了解，具有明显的优势。如果吸纳他们参与社区矫正工作，他们定会快速地掌握工作流程和方法，负责相应的工作。

2.2 储集空间特征

研究区栖霞组晶粒白云岩储层的储集空间以晶间（溶）孔、溶蚀孔洞为主，同时发育少量残余粒间/内孔、微裂缝、溶蚀扩大缝等。

晶间（溶）孔存在于各类晶粒白云岩中，是晶粒白云岩最重要的储集空间类型。孔隙直径为0.02～0.5 mm，一般孔隙周缘晶粒越大，晶间孔趋向于越大；晶粒自形程度越高，孔隙分布均质性越强。因此晶间孔在半自形—自形中粗晶白云岩中最为发育（图2g），而一些镶嵌状白云岩则基本不发育孔隙（图2b，2e），或仅发育少量残余粒间/内孔（图2a）。宏观上晶间孔呈针孔状，镜下面孔率为2%～8%,局部可大于10%。粉细晶白云岩中的晶间孔较小，并多被沥青充填或半充填（图2b）。中粗晶白云岩中的晶间孔或者被颜色明亮的白云石颗粒部分占据而呈平直状（图2g），推测这些白云石为孔隙形成后在孔隙空间内自由生长而成，或者内边缘具明显溶蚀加大痕迹（这种又可称为晶间溶孔，图2l）。

此篇文章主要是讲述了马铃薯脱毒种薯的情况，分析存在的不足，从而提出增强种薯生产经营管理，规范脱毒种薯生产情况，健全种薯繁育、推广和质量保证，从而增强建设种薯检测体系，改良检测方法，确保种薯质量的建议等。

3.2.1 白云石有序度

图3 川中地区栖霞组白云岩岩石学宏观特征Fig.3 Macroscopic petrological characteristics of dolomite of Qixia Formation in central Sichuan Basin

2.3 物性特征

对川中地区MX42井、MX117井等9口取心井的387个栖霞组样品进行物性分析（其中部分白云岩数据引自文献［13］），结果显示：栖霞组石灰岩孔隙度最大值为13.29%，最小值为0.04%，小于2%的样品占据石灰岩的85.7%；渗透率最大值为3.85×10-3μm2，最小值趋近于0，小于0.1×10-3μm2的样品占石灰岩的92.8%。少量邻近白云岩的石灰岩显示出较好的孔渗特征，但石灰岩整体致密，几乎没有储渗能力。白云岩孔隙度最大值为10.05%，最小值为0.42%，平均值为4.07%；渗透率最大值为24.37×10-3μm2，最小值趋近于0，平均值为3.12×10-3μm2。白云岩整体呈现为中等孔隙度、中高渗透率的特点。孔隙度和渗透率变化范围较大，非均质性较强，孔渗相关性较弱（图4），反映除了晶间（溶）孔外，溶蚀孔洞和裂缝也是非常重要的储集和渗滤空间。

图4 川中栖霞组岩石孔隙度与渗透率交会图(部分白云岩数据引自文献［13］)Fig.4 Crossplot of rock porosity and permeability of Qixia Formation in central Sichuan Basin（part of dolomite data from reference［13］）

2.4 白云岩分布特征

测井曲线特征及岩心特征显示栖霞组白云岩主要发育于2个三级层序的HST内（如MX42井，图5a），每个HST又划分为若干个高频层序，单一高频层序在岩心上表现为由上往下能量逐渐减弱。以MX42井4 651.85～4 653.90 m段为例（图5b），该段包括了顶部1.55 m的晶粒白云岩（图5b的照片①和②）和底部0.5 m的生物碎屑泥晶灰岩（图5b的照片③），越靠近顶部，溶蚀孔洞越发育，往下则逐渐减弱，至底部生物碎屑泥晶灰岩基本不发育溶蚀孔洞，整体构成了一个典型的高频旋回。需要注意的是，并不是每一次的高频旋回顶部都会发育白云岩，但白云岩均发育于高频旋回的顶部。这使得川中地区栖霞组白云岩在纵向宏观上呈现为单层厚度薄（单层厚度为2～9 m），与石灰岩叠置发育（总厚度3～20 m）的特点。如MX42井栖霞组自下往上发育了6层白云岩，但累计厚度却仅约17 m（图5）。

然而挑战同时存在，这一年Wine-Searcher经历了变革：网站创始人马丁·布朗（Martin Brown）不再参与公司运营。“我们拥有不少优秀的员工，经历了这件事情后，还要确保大家能像往常一样高标准地工作，对我而言是很大挑战。”David说道。现在，他们也在期盼着新的团队可以对网站进行一系列的完善与推进。

图5 川中MX42井栖霞组综合柱状图及取心段特征Fig.5 Comprehensive stratigraphic column of Qixia Formation of Well MX42 in central Sichuan Basin

前人研究普遍认为四川盆地栖霞组白云岩横向分布连续性较差［25］。但是，在川中地区，如GS112井、GS113井、HS2井、HS4井等井所在区域，以及GS001-X3井、GS1井、GS2井、GS18等井所在区域，栖霞组均发育数米到十几米的白云岩（图1a）。连井对比剖面（图6）显示，GS18井、GS16井、HS2井等井于SQ1的HST发育白云岩，其一般位于栖霞组中部；HS4井、MX31X1井等井则于SQ2的HST发育白云岩，一般位于距离栖霞组顶部10～20 m的位置。相邻井位白云岩发育的位置相对稳定，可对比性较强，表明川中地区栖霞组白云岩横向上在一定区域内仍具较好的连续性。

习近平指出：我们要建设的现代化是人与自然和谐共生的现代化，既要创造更多物质财富和精神财富以满足人民日益增长的美好生活需要，也要提供更多优质生态产品以满足人民日益增长的优美生态环境需要。

图6 川中地区栖霞组白云岩分布连井对比剖面图(连井剖面位置见图1a)Fig.6 Correlation section of the dolomite distribution of Qixia Formation in central Sichuan Basin(section position is shown in Fig.1a)

3 储层地球化学特征

3.1 样品与实验方法

U-Pb同位素测年被认为是目前唯一能用于测试碳酸盐岩矿物绝对地质年龄的方法［35］，尤其是近年来随着各种国际标样的建立而兴起的激光原位U-Pb同位素测年技术，以其高分辨率、高精度、快速等优势在解决多个地区的成岩演化问题中发挥了重要作用［26,36］。本次研究挑选了MX108-2样品（图9a），分别对其白云岩围岩和孔隙充填鞍状白云石胶结物进行了年龄测试，得到的白云岩围岩年龄为274.5±9.9 Ma（图9b），对应前人研究的栖霞阶年龄［16］；孔隙充填鞍状白云石胶结物年龄为262.7±7.5 Ma（图9c），与前人研究的茅口末期的峨眉山玄武岩喷发时期［37］为同一时期。

3.2 测试结果

钻井岩心中可见川中栖霞组中上部白云岩层段发育不同尺度的溶蚀孔洞，这是栖霞组另一种主要的储集空间。溶蚀孔洞孔径差异较大，主要孔径在0.2～5 cm之间不等。毫米级针状溶蚀孔洞局部均匀分布（图3a）；厘米级溶蚀孔洞则主要呈不规则状非均质分布，孔洞中多充填有巨晶方解石、中粗晶白云石、鞍状白云石、石英及沥青等（图3b）。成像测井揭示大溶蚀孔洞沿裂缝发育,以似水平状为主（图3c）。

白云石有序度与白云石形成环境的温度、盐度、Mg2+/Ca2+及成岩流体的其他性质等都密切相关，不同环境下白云石结晶的速度不同。白云石结晶速度越慢，有序度越高，反之则越低，故准同生白云岩具有较低的白云石有序度，而埋藏和热液环境下形成的白云石具有较高的有序度［27］。研究区栖霞组白云石有序度介于0.58～0.67，平均为0.62，白云石有序度整体较低。

3.2.2 微量元素、稀土元素

稳定碳氧同位素分析显示：栖二段石灰岩的δ13C值为2.85‰～5.33‰，均值为3.97‰；δ18O值为-7.05‰～-5.57‰，均值为-6.30‰。栖二段白云岩的δ13C值为3.84‰～5.35‰，均值为4.79‰；δ18O值为-8.34‰～-7.83‰，均值为-8.08‰。整体而言，白云岩和石灰岩的δ13C值均位于早二叠世海水范围［32］内，而δ18O值大致位于或略低于同期海水（图8）。

一个国家的制度安排与制度设计，只有深深植根于该国的文化历史与社会现实中，才能够具有旺盛的生命力。同样，一个国家的审计体制也必须与该国的政治体制、历史文化传统以及现实社会实践等相适应。中国特色国家审计制度是国家审计与中国特色社会主义制度相结合的产物，不仅要合乎审计的本质要求，而且要符合中国国情，体现中国特色社会主义政治制度的本质要求。组建中央审计委员会，为党统一指挥审计工作提供了制度依据和体制保障，是加强党的全面领导这一政治原则在审计工作的体现，是符合新时代中国特色社会主义政治制度本质要求的审计制度创新。〔本文受到南京审计大学首批政府审计学院中标课题(GASA161001)资助〕

表1 川中栖二段白云岩和石灰岩部分微量元素、稀土元素分析结果Table 1 Analysis results of trace and rare earth elements in dolomite and limestone of the Qixia Member 2 in central Sichuan Basin

图7 川中地区栖霞组石灰岩、白云岩稀土元素配分模式图Fig.7 REE distribution patterns of limestone and dolomite of Qixia Formation in central Sichuan Basin

δCe和δEu能较为准确地反映流体来源和成岩环境。因为La元素会明显影响δCe值，因此可结合δPr值判断Ce值异常的真实性［30］。研究区白云岩和石灰岩均具有δPr＞1、δCe＜1的特征，代表La元素正异常,Ce元素为轻微的负异常［30］。研究区白云岩δEu值为0.97～1.98，平均值为1.22，因此认为栖霞组白云岩存在δEu轻微正异常。

岩石中Y/Ho值与δGd值可以用来反映岩石成岩流体来源。海水中由于Y元素和Ho元素络合行为不同导致Y/Ho值较高，处于44～74之间，而受热液影响或近岸大气淡水影响的碳酸盐岩会显著低于该值［31］。研究区栖霞组白云岩Y/Ho值为43.41～67.65，平均值为56.66；石灰岩Y/Ho值为45.86～80.70，平均值为64.56。白云岩的δGd值为0.97～1.29，平均值为1.05；石灰岩的δGd值为0.93～1.33，平均值为1.14。二者基本符合海水成岩流体的特征，白云岩的Y/Ho值、δGd值略低于石灰岩。

3.2.3 稳定碳氧同位素

不同的成岩流体具有不同的微量元素、稀土元素特征，后期成岩作用的影响会导致海相碳酸盐岩微量元素、稀土元素的变化［28］。研究区栖霞组白云岩稀土元素值（∑REE）为0.86～3.75μg/g，平均值为1.85μg/g，轻、重稀土元素比值∑LREESN/∑HREESN（SN表示经过PAAS标准化后的稀土值，下同）为0.49～0.91，平均值为0.63（表1）；石灰岩的∑REE为0.58～2.20μg/g，平均值为1.37μg/g，∑LREESN/∑HREESN为0.32～0.69，平均值为0.51。整体而言，白云岩和石灰岩∑REE均偏低，白云岩呈现和石灰岩大体一致的稀土元素（REE）配分模式（图7），石灰岩∑REE和∑LREESN/∑HREESN值略低于白云岩（石灰岩表现为更明显的HREE富集）。

4.2.1 微古地貌高部位的颗粒滩是储层发育的物质基础

在天然湿度状态下的土层中开挖沟槽，且地下水位低于槽底时可开直槽，不设支撑，但对槽深有限制要求：砂土和砂砾土土层的深度不大于1.0 m；亚砂土和亚黏土土层的深度不大于1.25 m；黏土土层的深度不大于1.5 m。

图8 川中地区栖霞组石灰岩、白云岩碳氧同位素特征Fig.8 Carbon and oxygen isotopic characteristics of limestone and dolomite of Qixia Formation in central Sichuan Basin

需要注意的是，白云岩的Mn/Sr平均值为1.45（表1），平均值小于2，指示其更多保存的是原始海水的信息［40］，但远大于石灰岩（均值为0.12）的Mn/Sr值，指示其轻微受到后期成岩流体的改造。白云岩较石灰岩更低的δ18O值和更高的87Sr/86Sr值同样证明了这一点。研究区白云岩的δEu值呈现仅在极端还原高温环境下才会出现的正异常［41］，指示后期改造流体为热液流体。而白云岩阴极发光特征（图2i）和定年年龄表明，热液来自于茅口末期的峨眉山火山岩事件，且并没有对白云岩围岩造成重结晶程度的影响，仅于先存孔隙中形成亮色白云石环边（图2g，2h）和鞍状白云石（图2k）。

锶同位素87Sr/86Sr值不受温度、压力和微生物作用影响，其长期不变且容易辨别，因而在推断碳酸盐岩成岩流体来源中得以广泛应用［33］。川中地区栖二段石灰岩的87Sr/86Sr值在0.707 1～0.708 0之间，平均值为0.707 5；白云岩的87Sr/86Sr值在0.707 7～0.708 7之间，平均值为0.708 0（部分白云岩数据引自文献［5］）。通过对比全球二叠纪海水87Sr/86Sr值［34］发现，石灰岩87Sr/86Sr值落于同期海水范围内，白云岩则基本落于或略高于同期海水范围。

“裂缝”夫妻中有很多都是这样的情况：妻子早年跟随丈夫吃过很多苦头，到晚年丈夫仍没有任何回报的表示，让妻子继续为自己忙上忙下，因此，妻子对丈夫的不满有点外力就会爆发。所以，男士要对糟糠之妻常怀感恩之心，时不时地在生日时买点礼物，在家庭聚餐时当众表达感谢之情等。双方不要都暴躁，这样的局面不好控制。

3.2.5 碳酸盐岩U-Pb同位素年龄

本次地球化学分析测试样品取自于研究区MX42井、MX108井、NJ井等9口取心井。为方便对比，对白云岩和石灰岩都进行了地球化学测试。所有分析测试均完成于中国石油集团碳酸盐岩储层重点实验室。其中，阴极发光：样品数23，仪器型号为CL8200 MK5(配以Leica偏光显微镜)，工作条件为9～11 kV和400～450 mA；白云石有序度：样品数13，全岩样品，测试仪器为X射线衍射仪X’pert Pro，在温度25℃、湿度40%RH的环境下测试完成；微量元素和稀土元素：石灰岩样品数8，白云岩样品数9，全岩样品，分析仪器为PANalytical Axios XRFX射线荧光光谱仪，所有样品稀土值采用后太古宙澳大利亚页岩（PASS）为标准进行标准化；碳氧稳定同位素：石灰岩样品数9，白云岩样品数8，全岩样品，仪器为同位素比质谱仪Delta V Advantage，应用标准100%磷酸法测试，并通过白垩纪Pee Dee组箭石（VPDB）标准化，测试精度为δ13C＜0.01‰、δ18O＜0.02‰；锶同位素：石灰岩样品8，白云岩样品数2，全岩样品，仪器为热电离同位素比质谱仪TRITON PLUS，检测环境温度为22℃、相对湿度为66%，误差范围为0.000 014～0.000 084。上述地球化学实验样品均不包含裂缝和溶洞矿物。激光U-Pb定年分析测试采用激光微区选点，单一组构内进行70次单点剥蚀，具体仪器、方法和流程见文献［26］。

图9 川中地区栖霞组白云岩围岩和孔隙充填鞍状白云石U-Pb年龄Fig.9 U-Pb age of dolomite of surrounding rock and pore-filling saddle dolomite of Qixia Formation in central Sichuan Basin

4 讨论

4.1 白云岩成因

综合白云岩宏观、微观特征及地球化学与定年分析，认为川中栖霞组白云岩主要形成于准同生期，浅埋藏期受到热液作用改造。

首先，从岩石学特征来讲，可见白云岩和石灰岩截然接触的现象，如MX42井取心段4 653.46 m处（图5b），这说明白云石化作用至少发生在大规模压溶之前。在HS2井4 537～4 538 m井段，发育有代表准同生环境的泥粉晶白云岩［38］（图2j），虽然泥粉晶白云岩在川中栖霞组中少见，但它指示川中地区存在准同生白云岩发育环境。从地球化学角度来看，川中地区栖霞组白云岩有序度整体较低，符合准同生白云石化的特点。白云岩和石灰岩的稀土含量整体较低，且均为HREE富集，整体表现为与海水大体一致的稀土元素配分模式（图7），另外白云岩、石灰岩均与海水的La富集、Ce亏损、Gd轻微富集、Y显著富集等典型稀土元素特征相吻合［39］，反映石灰岩的成岩流体和白云石化的成岩流体来源相同，都来源于海水。基本落于二叠纪海水范围的δ13C值和87Sr/86Sr值特征也能很好地证明这一点。另外海源流体因为贫Fe和Mn（表1），形成的岩石以昏暗发光或不发光为主［28］，研究区白云岩的阴极发光特征与其相吻合（图2f）。在年代学方面，白云岩围岩定年年龄为（274.5±9.9）Ma，栖霞阶年龄为（283.5±0.6）Ma～（272.95±0.11）Ma［16］，白云岩年龄落于栖霞阶内，这是白云岩围岩形成于准同生期的最直接证据。

3.2.4 锶同位素

4.2 白云岩储层发育的主控因素

XRF法测定硅钙合金中Si、Ca、Fe、Al已有报道，本文结合本实验室设备条件和Si、Ca、Fe、Al的XRF分析报道[10,12]的实验条件，对各元素的分析条件进行了选择，所确定的分析条件见表1。

由于白云石化作用以准同生为主，白云石化持续时间较短，部分生物碎屑难以白云石化，使得其在白云岩中得以残留（图2b），而完全白云石化的白云石在一定的透射光处理后也能恢复其生物碎屑和砂屑结构（图2c，2d），指示白云岩原岩为颗粒滩生物灰岩。颗粒滩生物灰岩具有较高原始孔隙度和大量不稳定矿物溶蚀孔隙，这些孔隙是成岩流体进入的通道［42］，白云岩晶间（溶）孔也继承自这些孔隙（图2a，2c），因此颗粒滩是白云岩储层发育的物质基础。另一方面，沉积时期微古地貌高部位因其水体能量较高，因而是颗粒滩发育的优势地形。

4.2.2 高频层序界面控制了准同生岩溶作用和早期白云石化

白云岩单层厚度薄、层数多，与石灰岩交替叠置发育，整体构成多个旋回，单一旋回对应四级乃至更高级次的海平面升降（图5）。高位域海平面下降期间，微古地貌高部位颗粒滩容易暴露，接受大气淡水溶蚀，形成孔隙发育层［43］和似水平状发育的溶蚀孔洞（图3c），并呈现层序界面顶部溶蚀孔洞较大，往下逐渐变小的特点（图5b）。白云岩较石灰岩偏负的δ18O值较好地证实了这一点。同时，海平面下降必然导致古地貌高部位水体逐渐受限，因蒸发作用海水中的Mg2+逐渐富集，富Mg2+流体沿着先存溶蚀孔洞进入颗粒滩灰岩，引发颗粒滩灰岩的渗透回流白云石化（图10a），这也是研究区溶蚀孔洞主要存在于白云岩中的原因。因此是先存孔隙的存在导致白云石化更易发生，而不是白云石化作用形成孔隙［44］。

涉众型经济犯罪指的是在经济案件中涉及的受害人比较多，而且受害群体并非特定性，其中包括传销、集资诈骗等犯罪行为。这种犯罪行为影响范围广，造成的损失大，所以必须要加强对犯罪的防控措施进行分析。

4.2.3 准同生白云石化作用有利于孔隙保存

石灰岩发生白云石化作用后，其成分和结构组成均被改变，并变得更加抗压实、抗压溶，故准同生白云石化作用能有效保存原生孔隙和次生溶孔。另外前人研究表明，深埋藏时期白云岩较石灰岩更易溶从而更容易进一步改善孔渗［45］。前文叙及，研究区白云岩形成于准同生时期，这也是栖霞组自沉积后经历了复杂的成岩作用和构造作用，白云岩晶间（溶）孔和溶蚀孔洞仍能得以保存的原因。这些孔隙也为茅口末期的热液进入栖霞组提供了通道，形成了以破坏孔隙为主的白云石环边和鞍状白云石（图10b）。

图10 川中栖霞组白云岩储层发育模式Fig.10 Development model of dolomite reservoir of Qixia Formation in central Sichuan Basin

4.3 白云岩储层的发育规律

基于前文论述，沉积时期微古地貌高地是颗粒滩发育的优势地形，高位域海平面下降期间，高部位颗粒滩最容易暴露被溶蚀形成早期孔隙，这些孔隙为早期白云石化提供了通道，而早期白云石化又有利于早期孔隙的保存。另一方面，沉积时期微古地貌高部位水体能量高，其沉积厚度会显著高于周围，因此栖霞组高位域沉积厚值区是白云岩储层有利勘探区。

以川中地区的合川三维地震工区为例，分别采用栖霞组内部最大波峰振幅平面属性图和顶界最大波峰振幅平面属性图来预测栖霞组SQ1和SQ2两个三级层序高位域时期的白云岩储层。基于最大波峰振幅越强，地层厚度越小的原则［46］，栖霞组内部最大波峰振幅平面属性图（图11a）指示栖霞组沉积厚值区位于GS18井、MX31X1井、GS16井附近以及HS4井以东区域，即这些区域在沉积时期为微古地貌高部位，据此认为这些区域是栖霞组SQ1高位域白云岩储层发育的有利区域。同理，SQ2高位域白云岩储层主要发育于HS4井和MX31X1井附近（图11b）。实际钻井结果也证实，这些井发育多层数米厚、累厚达十多米的白云岩（图6），其中，GS18井栖霞组日产气41.74×104m3，MX31X1井栖霞组日产气70×104m3，GS16井栖霞组日产气5.46×104m3，HS4井栖霞组日产气45.65×104m3。钻探结果证实了栖霞组通过寻找沉积厚值区从而预测白云岩储层的勘探思路的可靠性，该思路也从另一个角度解决了栖霞组白云岩储层单层厚度薄，地震相上难以识别的局限性。

图11 川中栖霞组最大波峰振幅属性平面图（三维工区位置见图1a）Fig.11 Seismic attribute maps of maximum peak amplitude of Qixia Formation in central Sichuan Basin(the 3D seismic survey is shown in Fig.1a)