川北元坝地区长兴组白云石化作用机制及其对储层形成的影响*

2014-05-30孟万斌武恒志李国蓉张小青吕正祥

岩石学报 2014年3期

关键词：长兴白云石微晶

孟万斌武恒志李国蓉张小青吕正祥

1．油气藏地质及开发工程国家重点实验室，成都 610059

2．中国石化西南油气分公司，成都 610041

3．中国石化西南油气分公司勘探开发研究院，成都 610041

4．成都理工大学能源学院，成都 610059

碳酸盐岩储层是油气勘探领域极为重要的一类储层，世界上超过60%的油气产量都来自碳酸盐储层，其中多半是白云岩储层(马永生等，2007)。四川盆地东北部是中国近年来油气重点勘探地区，上二叠统长兴组－下三叠统飞仙关组发育的礁滩白云岩储层是该区油气的主力储层之一(Ma et al．，2007，2008)。元坝地区长兴组储层主要为台地边缘浅滩相及台地边缘礁滩相沉积，储层岩性以亮晶生屑灰岩及亮晶生屑白云岩为主，储层物性较好(郭彤楼，2011)，详细研究白云石化作用机制及其对储层发育的影响，对该区进一步勘探将起到很好的指导作用。

1 地质背景

元坝地区位于四川盆地北部，苍溪、旺苍、巴中市交汇区域。该区在构造上属于川中平缓褶皱带北缘的一部分，位于龙门山北段前缘、九龙山背斜构造带南翼、通南巴背斜构造带西南侧，地层产状平缓，构造变形弱。

关于长兴组时期的沉积格局，目前仍有分歧。王一刚等(2005，2006a，b，2009)认为，川东北地区自长兴组时期发育与鄂西海盆相连的“开江－梁平海槽”，其北部与深水陆棚相连，海槽两侧发育台地。马永生等(2006a，b，2007)则认为这一海槽并不存在，它只是水体相对较深的碳酸盐岩缓坡，长兴组只是水体相对台地较深的“台棚”环境产物，台棚西南向台地具有开阔海台地性质，东北向台地具有局限台地性质。显然，无论是“海槽”论者，还是“台棚”或“陆棚”论者，他们均认为旺苍－巴中－宣汉－开江－梁平这一相对较深水区域向东南方向仅延伸至梁平附近，未与更东南边的鄂西海槽(海盆)贯通，在东南方向，四川盆地西部碳酸盐台地与东北部的台地是连在一起的。通过对有关区调资料、钻井资料和地震资料的查阅和分析，开江－梁平这一区域长兴组－飞仙关组的沉积物反映出水体明显较鄂西海槽(海盆)区域浅。结合前人(马永生等，2006a，b，2007)的研究成果，本文认为使用“陆棚”来描述它较为适宜。元坝地区处于开江－梁平台棚体系的西北部，横跨了西部台地内部－台地边缘－斜坡－陆棚等沉积相带部位(图1)。

元坝地区二叠纪地层发育较为齐全，由下而上包括了下二叠统梁山组、中二叠统栖霞组和茅口组、上二叠统吴家坪组和长兴组。长兴组与下伏吴家坪组和上覆飞仙关组之间呈整合接触，是上二叠统上部的一套浅海台地相碳酸盐岩地层，厚度200～300m，以微晶灰岩、生物碎屑灰岩、礁灰岩、白云岩发育为特征。

表1 元坝地区长兴组不同类型白云石的X衍射有序度分析结果Table 1 XRD of different dolomites in the Changxing Formation in Yuanba area

2 白云石矿物学特征

元坝地区长兴组发育微晶白云石(岩)、粉－中晶他形白云石、粉－中晶自形白云石、异形白云石等四种类型。

2．1 微晶白云石(岩)特征

此类白云石在区内长兴组碳酸盐岩中较为少见，白云石晶体细小，岩石整体由微晶白云石组成构成微晶云岩(图2a)，也可以是微晶白云石在岩石中呈斑状富集形成微晶白云石条带或斑块(图2b)，但以斑块状分布为主。岩石可受到缝合线和各期破裂作用裂缝的切割以及各期溶蚀作用的改造。白云石有序度介于0．46～0．64，平均值为0．556(表1)，显示出有序度低的特征，表明白云石化作用发生于高盐度环境，是一个快速结晶的过程。

2．2 粉－中晶他形白云石特征

此类白云石在区内长兴组碳酸盐岩中常见，白云石晶体可由粉晶到中晶级，晶体自形程度普遍较差，多呈他形晶，晶体表面普遍较脏，个别保留了原来颗粒结构的阴影，晶体之间往往为镶嵌状接触(图2c，d)。白云石晶体仍可受到缝合线、裂缝以及溶蚀作用的改造。白云石有序度介于0．51～0．95，但主要在0．51 ～0．75 范围内变化，平均值为 0．678，有序度较低(表1)，表明此类白云石仍是高盐度环境下成岩作用的产物，白云石结晶速度仍然较快。

图1 四川盆地北部上二叠统长兴组时期沉积格局及研究区位置图Fig．1 Study area and sedimentary framework of the Upper Permian Changxing Formation in the norhthern Sichuan Basin

图2 元坝地区长兴组中常见白云石类型图谱(a)－微晶云岩，白云石化作用完全彻底YB204－13(2/65);(b)－砂屑微晶云岩，发育成岩缝，YB2－6(2/18);(c)－粉晶云岩，由粉晶他形白云石组成，发育缝状溶蚀孔洞，边缘有细晶自形白云石生长，YB9－7(3/27);(d)－生屑云岩，发育粉－中晶他形白云石，白云石化后海百合碎片保持原始形态，YB27－8(2/21);(e)－粉－细晶云岩，由粉－细晶自形白云石组成，溶蚀孔洞发育，YB29－2(5/18);(f)－藻粘结云岩，粘结部位为微晶白云石，窗状孔部位为粉晶自形白云石，YB101－3(30/30);(g)－藻粘结灰岩，具弱波状消光的异形白云石呈斑状分布，YB22－4(9/90);(h)－粉－中晶云岩，溶蚀孔洞边缘发育波状消光异形白云石，YB102－7(2/21)．除照片(g)为正交偏光外，其余均为单偏光Fig．2 Thin section photomicrographs for common dolomites in the Changxing Formation of the Yuanba area

表2 元坝地区长兴组白云石碳、氧稳定同位素分析结果统计Table 2 Carbon and oxygen isotope values of different dolomites in the Changxing Fromation of the Yuanba area

2．3 粉－中晶自形白云石特征

此类白云石在区内长兴组碳酸盐岩中最为常见。白云石以分散状、斑状为主要产出形式，也可呈整体云化、偶见分布等形式存在。白云石晶体多呈自形晶，部分为半自形晶，晶体较粗，多为粉－细晶级，部分达中晶级，表面较脏或较干净、或显雾心亮边结构(图2e，f)。白云石可交代改造胶结物方解石及缝洞充填的方解石，被缝合线及三期构造破裂作用及裂缝充填方解石切割，被第二期、第三期溶蚀作用改造。白云石有序度介于0．70 ～0．98，平均值为0．863，有序度较高(表1)，表明此类白云石可能形成于盐度较低环境条件下，白云石结晶速度较为缓慢。

2．4 异形白云石特征

此类白云石在区内长兴组部分井段见到，白云石在岩石中可呈斑状、分散状、偶见等三种形式存在，但以斑状分布最为常见。白云石晶体大小差别较大，可达粉晶－粗晶级，部分晶体晶面弧形特征清楚，表面较干净，正交光下具波状消光或弱波状消光特征(图2g，h)。

3 白云石(岩)地球化学特征

3．1 碳、氧稳定同位素组成特征

据45件样品碳、氧同位素分析结果(表2)，长兴组白云石的碳、氧同位素具有以下特点:

一是δ13C(PDB标准，下同)值均为正值且普遍较高。三种白云石(岩)的平均值分别为 2．65‰、4．23‰和 3．83‰，它们与同时期的正常海相灰岩背景值(3．156‰)(党录瑞等，2011)较为接近或略高，揭示了白云石的碳来源于长兴组时期的正常海水或者高盐度海水，或者来源于长兴组海相碳酸盐岩地层本身，缺乏有机质氧化产生的轻CO2和大气淡水的加入。

二是白云石的δ18O(PDB标准，下同)值总体较高，且由微晶白云石→粉－中晶他形白云石→粉－中晶自形白云石，δ18O值逐渐降低，但差距较小。长兴组白云石的δ18O值在－5．85‰ ～－3．32‰范围，主要在－6‰ ～－5‰区间，变化范围狭窄，平均值为－5．19‰;长兴组正常海相微晶灰岩的δ18O 值为－6．30‰ ～－4．79‰，平均值为－5．68‰。白云石的δ18O值与同时期的正常海相灰岩背景值仍非常接近，表明白云石形成与正常海水、高盐度海水，或者与正常海水、高盐度海水相似的地层水条件下，而且成岩环境温度不高。

Keith and Wbber(1964)(转引自郭金春等，2008)用δ13CPDB和δ18OPDB结合指示古盐度，即用公式Z(盐度指数)=2．048(δ13CPDB+50)+0．498(δ18OPDB+50)计算 Z 值，Z ＞120为海相，Z＜120为陆相。本次研究利用区内长兴组白云石的碳、氧同位素分析结果进行了Z值计算，结果表明Z值在128．57 ～136．02，均大于 120，也说明成岩流体为海水或与海水性质相似的地层水。

利用碳酸盐矿物氧同位素温度计算公式T1=13．85－4．54δ18OPDB+0．04(δ18OPDB)2(余志伟，1999)计算的长兴组白云石的形成温度在29．38～41．81℃范围内变化，也表明白云石形成温度不高，白云石化是发生在地表、近地表，或者非常浅的埋藏环境中。

上述白云石碳氧同位素特征反映出微晶白云石是在蒸发条件下盐度更高的环境中形成;粉－中晶他形白云石、粉－中晶自形白云石是在盐度正常或略高的环境中形成，而粉－中晶自形白云石即使是在埋藏成岩环境中形成，但埋藏深度也非常浅;由微晶白云石→粉－中晶他形白云石→粉－中晶自形白云石可能是一个沉积埋藏过程中多阶段白云石化作用的产物。

3．2 锶同位素组成特征

根据黄思静等(2011，2012)的研究，四川盆地上二叠统长兴组正常海相碳酸盐岩的锶同位素值非常低，在0．707048～0．707112范围内变化。

39件白云石和微晶灰岩样品的锶同位素地球化学特征表现出较为复杂的特征(表3):①2个正常海相微晶灰岩样品87Sr/86Sr比值为 0．706576 ～ 0．707713，平均值为0．707145，处于同期海水87Sr/86Sr比值范围，样品的岩石学观察发现样品缺乏后期改造特征，因此该值可以作为区内长兴组正常海相碳酸盐岩的锶同位素背景值。②2个微晶白云石样品的锶同位素值为0．707480和0．708314(平均值为0．707897)，与长兴组正常海相微晶灰岩锶同位素背景值较为接近，表明其继承或保存了同生期海水白云石化的锶同位素值特征。③25个粉－中晶他形白云石或以其为主的样品和9个粉－中晶自形白云石或以其为主的样品的87Sr/86Sr比值变化有三种情况。第一种是87Sr/86Sr比值小于0．706000，两类岩石样品的平均值分别为0．705337和0．705104，其锶同位素值较低，结合岩石学特征，这些层段或相邻层段受到热液作用的改造。第二种是87Sr/86Sr比值介于0．706000～0．710000，两类岩石样品的平均值分别为 0．707134 和0．707351，与长兴组正常海相微晶灰岩的锶同位素背景值非常接近，是同生期海水白云石化锶同位素特征的反映，它们形成于与长兴组海水相似的地层水环境中。第三种是87Sr/86Sr比值大于0．710000，锶同位素值较高，两类岩石样品的平均值分别为0．714747和0．713152，暗示还另有盆地内部碎屑岩地层热液流体对长兴组白云岩进行过改造。

表3 元坝地区长兴组不同类型白云石锶同位素分析结果统计Table 3 Strontium isotope values of different dolomites in the Changxing Formation of the Yuanba area

图3 元坝地区长兴组白云石稀土元素组成及配分模式(a)－微晶白云石类;(b)－粉－中晶他形白云石;(c)－溶蚀孔洞较发育的粉－中晶他形白云石;(d)－粉－中晶自形白云石Fig．3 REE patterns of different dolomites in the Changxing Formation of the Yuanba area

总之，锶同位素分析结果揭示出区内长兴组的微晶白云石、粉－中晶他形白云石、粉－中晶自形白云石都是在长兴组正常海水或者与长兴组正常海水相似的地层水环境中形成的，但形成后又遭受过深部热液或者盆地内部碎屑岩地层热液的改造作用。

3．3 稀土元素组成及配分模式

岩石或矿物中的稀土元素组成往往是成岩流体及物质来源、形成环境的指示剂，特别是识别热液作用的重要标志。稀土元素Eu正异常、Ce负异常通常指示了高温热液作用的发生(丁振举等，2000;冯胜斌等，2007;胡文塇等，2010)。

利用长兴组41件白云岩及白云石斑块样品的稀土元素分析结果，以球粒陨石为标准对各样品的稀土元素数据作图，并与正常海相微晶灰岩的结果进行比较，发现不同类型白云石的稀土元素组成具有明显的差异(图3)。

①微晶白云石类的稀土元素组成及配分模式特征与长兴组正常海相微晶灰岩相似，可相对于球粒陨石贫稀土元素或富稀土元素，具有相对富集轻稀土元素、贫重稀土元素、略显Ce负异常特征(图3a)，揭示了此类白云石是在正常海水环境或者海水浓缩环境中形成。在海相环境中沉积成岩的灰岩，除了重稀土有少量亏损外，较好地保存了原始海水的稀土元素配分型式(Webb and Kamber，2000;Kamber and Webb，2001)。②粉－中晶他形白云石类的稀土元素组成及配分模式也与正常海相微晶灰岩相似，但绝大多数样品相对于球粒陨石贫稀土元素，显示出相对富集轻稀土元素、贫重稀土元素、Ce负异常的特征(图3b，c)，揭示此类白云石应是在正常海水或者浓缩海水环境中形成的。而个别粉－中晶他形白云石样品，与前述特征有一定差异，突出表现在Eu、Sm正异常，Ce负异常，结合该层段存在热液成因天青石等矿物，认为它们是粉－中晶他形白云石后期遭受深部热液作用改造的反映。③粉－中晶自形白云石也具有两种稀土元素组成及配分模式(图3d)，有2个样品的稀土元素组成与正常海相微晶灰岩相似，相对于球粒陨石贫稀土元素，曲线较平滑，略显相对富集轻稀土元素、贫重稀土元素特征，揭示了其是正常海水环境或与正常海水环境相似的地层水环境成岩作用的产物。然而，大多数样品虽仍显示出相对球粒陨石贫稀土元素，但曲线特征及配分模式明显不同于正常海相微晶灰岩，突出表现在具有Eu、Sm正异常，Ce弱负异常特征，结合这些井段或相邻井段存在热液成因矿物的现象，认为它们是粉－中晶自形白云石形成后受到热液改造的响应特征。

综上所述，稀土元素分析揭示出区内长兴组白云石主要是在浓缩海水环境、正常海水环境、或者与正常海水相似的地层水环境中形成的，并在形成后遭受过热液地质作用的改造。

4 白云石化作用机制与模式

4．1 白云石化作用与层序、沉积相的关系

4．1．1 长兴组白云石化作用与层序的关系分析

通过分析白云岩分布与三级层序的关系以及不同类型白云石与高频层序的关系，发现长兴组白云石化作用与层序的关系具有以下特点:①白云岩及其白云石化作用主要在高位体系域发育(图4);②不同类型的白云石与高频层序的关系不一致(图5)。微晶白云石主要在五级层序的上部发育，其次是在五级层序的中部发育;粉－中晶他形白云石和粉－中晶自形白云石在五级层序的上部和中部都优势发育，但也可以延续到五级层序的下部;而异形白云石与五级层序不存在特定的关系。

4．1．2 长兴组白云石化与沉积相带的关系

图4 元坝地区长兴组白云岩及白云石化作用在三级层序体系域中发育情况统计Fig．4 Dolomites and dolomitization distribution in different system tracts

图5 元坝地区长兴组不同类型白云石与高频层序的关系Fig．5 The relationships of dolomite types with highfrequency sequence

认识白云石化与沉积相带的关系，有助于对白云岩成因机制的认识和对白云岩发育分布的预测。元坝地区白云岩发育分布与沉积相带具有以下关系:

①白云岩主要在台地相区发育，但与沉积相关系复杂。已有的钻井揭示长兴组白云岩主要在台地区发育，一般台地边缘白云岩厚度较大，台地内部白云岩厚度较小，但也存在着非均质性和复杂性。②不同类型白云石与沉积微相的关系存在明显差异。微晶白云石发育在生物礁、生屑滩、潮坪等沉积微相中，其中潮坪微相最有利于微晶白云石的发育。粉－中晶他形白云石发育在生物礁、生屑滩、潮坪、滩间、礁间微相中，生屑滩及滩间微相是最有利的沉积微相。粉－中晶自形白云石在生物礁、生屑滩、滩间、潮坪等沉积微相中发育，生屑滩及滩间微相最有利于粉－中晶自形白云石的发育。异形白云石在生物礁、生屑滩、滩间、礁间、潮坪微相中都有出现，但主要在生物礁及潮坪相沉积物中发育。

图6 长兴组白云石化模式Fig．6 Dolomitization models of the Changxing Formation

4．2 长兴组白云石化作用机制与模式

4．2．1 微晶白云石的形成机制与模式

微晶白云石的岩石学特征和地球化学特征揭示了其形成于盐度高的浓缩海水、地表或近地表低温条件下的早期成岩作用阶段，五级层序中上部的潮坪、生物礁、生屑滩部位是有利于微晶白云石化作用发生的部位。因此，可用同生期蒸发泵白云石化模式来解释该类白云石的形成(图6a)。进入五级层序中上部时期，随着海平面不断降低，礁滩的不断生长并部分出露海面，在局部形成潮坪环境，海水通过毛细管作用向生屑滩顶部、生物礁顶部、潮坪沉积物渗透运动，并发生蒸发浓缩，导致白云石化作用发生并形成微晶白云石。

4．2．2 粉－中晶他形白云石的形成机制与模式

岩石学特征表明，粉－中晶他形白云石形成于同生期的成岩作用阶段，其地球化学特征揭示，白云石作用发生于盐度较高或正常、地表或近地表低温条件。该类白云石在五级层序的上部和中部优势发育，在五级层序下部也较为发育，在生物礁、生屑滩、滩间、礁间、潮坪等部位均可发育。这些特征表明它们是同生期蒸发泵白云石化形成的微晶白云石经埋藏阶段重结晶作用的结果;而生物礁、生屑滩及潮坪下部及滩间、礁间部位的粉－中晶他形白云石则是渗透回流白云石化作用的结果(图6b)。在五级层序甚或六级层序的中晚期，在生物礁、生屑滩、潮坪沉积物上部发生蒸发泵白云石化作用的同时，由于环境受限，蒸发作用强，水体盐度增高、比重增大，水体在重力作用下向下运动，渗透进入到生物礁、生屑滩、潮坪沉积物的下部及滩间、礁间沉积物中，导致这些沉积物发生白云石化，由于成岩介质盐度较高，白云石化速度较快，因而形成的白云石往往晶形不好，表面较脏。

4．2．3 粉－中晶自形白云石的形成机制与模式

粉－中晶自形白云石的岩石学特征和地球化学特征揭示，该类白云石形成于成岩早期阶段流体性质与长兴组正常海水相近或盐度稍高、埋藏较浅的环境，白云石化作用可能是一个与渗透回流白云石化相连续的过程(图6c)。白云石化发生于沉积物刚刚沉积后的浅埋藏过程中，成岩流体为沉积物中封存于孔隙中的海水，可以是正常盐度海水或者与渗透回流白云石化流体相似的海水，镁离子来源于封存海水和由文石、高镁方解石稳定化过程中释放出来的镁离子，压实作用驱动孔隙水流体在沉积物中渗透运动，进而导致沉积物发生白云石化。由于镁离子的供给受到一定限制，白云石化作用过程缓慢且往往不彻底，从而形成晶形较好的白云石，并呈整体云化、斑状分布、分散状分布、偶见分布等多种产状。

图7 元坝地区长兴组白云岩孔隙度、渗透率分布直方图Fig．7 Porosity and permeability histogram of the dolomites in the Changxing Formation of the Yuanba area

4．2．4 异形白云石的形成机制与模式

异形白云石的产出状态和时空分布与其他几类白云石的分布、长兴组的层序和沉积相之间没有密切的关系，揭示了该类白云石的形成是成岩晚期的一次成岩事件，成岩流体来源于长兴组地层外部，流体来源及其性质与其他三类白云石明显不同。异形白云石通常被认为是热液地质作用的产物，结合区内长兴组地层溶蚀孔洞内见有石英、天青石等典型热液矿物的沉淀和交代现象，元素分析揭示部分样品Sr、Ba含量异常高，第二期构造裂缝方解石、溶蚀缝洞方解石、自形和异形白云石、受改造的他形白云石和微晶灰岩均显示了Eu正异常、Ce负异常特征，也指示了区内长兴组碳酸盐岩地层遭受过热液地质作用的改造，而有关白云石(白云岩)或灰岩的锶同位素值既有明显低于长兴组正常海相微晶灰岩的，也有明显高于长兴组正常海相微晶灰岩的，则暗示了这种外部来源热液既有地幔深部来源，又有盆地内部碎屑岩地层来源。因此，可以认为异形白云石由成岩晚期热液地质作用机制形成(图6d)。伴随着燕山运动的发生，断裂及裂缝系统及其减压带的形成，地下深部流体或盆地内部碎屑岩地层流体沿断裂和裂缝通道进入到长兴组碳酸盐岩地层，进而沿着先前的孔洞缝系统进行渗透流动，导致异形白云石的形成，并伴随石英、天青石、方解石的沉淀和交代作用，还可能产生热液溶蚀作用并形成溶蚀孔、洞、缝。

4．3 长兴组白云石化作用与储渗空间的形成

白云石化作用被认为有利于储层储渗空间的形成，是碳酸盐岩储层发育的建设性因素，白云岩往往构成一类重要的优质碳酸盐岩储层。或因白云石化作用是一个岩石体积减小或孔隙体积增加的作用过程(Murray，1960;Weyl，1960;Lucia and Major，1994;Sun，1995)，或因白云石化机制与成岩环境所控制的溶蚀作用(Runnells，1969;Badiozamani，1973;Hardie，1987;Warren，2000)，或者白云岩层抗压能力较强，有利于先期孔隙的保存，相同深度条件下，白云石化灰岩及白云岩相对于灰岩及泥质灰岩有更高的孔隙度(Purser et al．，1994;Brown 1997)。

已有的观察及研究表明，元坝地区长兴组储层及储渗空间主要在白云岩及白云石斑块部位优势发育，揭示了白云石化作用是区内长兴组优质储层形成的重要基础。据岩芯样品统计，区内长兴组白云岩孔隙度主要分布在5% ～10%区间，其次在2% ～5%区间;渗透率主要分布在0．1×10－3μm2～10×10－3μm2区间，其次在 0．001 ×10－3μm2～0．1 ×10－3μm2区间(图7)，反映了区内长兴组白云岩储层孔渗性较好。区内长兴组在同生期发育蒸发泵和渗透回流白云石化，并由此影响了沉积物孔隙水的性质，影响浅埋藏白云石化作用的发育和分布，导致不同成因类型白云石在一个统一的区域内共存分布，形成白云岩和白云石斑块，这些白云石化本身有利于晶间孔等储渗空间的形成，使其早期就有相对较好的储渗条件，有利于后期具溶蚀能力的流体发生溶蚀作用;另外，中深埋藏条件下，白云石的溶解速率大于方解石，强烈的重结晶作用和溶蚀作用使得岩石形成大量的次生孔隙(王瑞华等，2007)，从而形成很好的储层。

5 结论

(1)元坝地区长兴组发育微晶白云石(岩)、粉－中晶他形白云石、粉－中晶自形白云石、异形白云石等四种白云石类型。

(2)长兴组微晶白云石(岩)、粉－中晶他形白云石、粉－中晶自形白云石主要是在浓缩海水环境、正常海水环境、或者与正常海水相似的地层水环境中形成的，并在形成后遭受过热液地质作用的改造。由微晶白云石→粉－中晶他形白云石→粉－中晶自形白云石是一个沉积埋藏过程中多阶段白云石化作用的产物。

(3)各种不同类型的白云石由不同的机制形成，微晶白云石由同生期蒸发泵白云石化作用形成，粉－中晶他形白云石由同生期渗透回流白云石化作用形成，粉－中晶自形白云石由成岩早期浅埋藏状态下地层水白云石化作用形成，异形白云石是成岩晚期热液白云石化作用的产物。

(4)本区白云石化作用有利于早期储渗空间和后期次生孔隙的形成，因此是有利于长兴组储层发育的建设性成岩作用。

Badiozamani K．1973．The Dorag dolomitization model－application to the Middle Ordovician of Wisconsin．Journal of Sedimentary Petrology，43(4):965－984

Brown A．1997．Porosity variation in carbonates as a function of depth:Mississippian Madison Group，Williston Basin．In:Kupecz JA，Gluyas J and Bloch S(eds．)．Reservoir Quality Prediction in Sandstones and Carbonates:AAPG Memoir 69，29－46

Dang LR，Zheng RC，Zheng C，Wen QB，Chen SH and Liao J．2011．Origins and diagenetic system of dolomite reservoirs in the Upper Permian Changxing Fm．，eastern Sichuan Basin．Natural Gas Industry，31(11):17－52(in Chinese with English abstract)

Ding ZJ，Liu CQ，Yao SZ and Zhou ZG．2000．Rare earth elements compositions of high temperature hydrothermal fluids in sea floor and control factors．Advance in Earth Sciences，15(6):307 －312(in Chinese with English abstract)

Feng SB，Zhou HR，Yan CH，Peng Y，Yuan XQ and He J．2007．Geochemical characteristics of hydrothermal cherts of Erlangping Group in EastQinling and theirgeologic significance．Acta Sedimentologica Sinica，25(4):564－573(in Chinese with English abstract)

Guo JC，Ma HZ，Song EY and Chen N．2008．Applications of lacustrine carbonate in paleoenvironment research．JournalofSaltLake Research，16(2):66－72(in Chinese with English abstract)

Guo TL．2011．Reservoir characteristics and its controlling factors of the Changxing Formation reservoir in the Yuanba gas field，Sichuan Basin，China．Acta Petrologica Sinica，27(8):2381 －2391(in Chinese with English abstract)

Hardie LA．1987．Dolomitization:A critical view of some current views．Journal of Sedimentary Research，57(1):166－183

Hu WX，Chen Q，Wang XL and Cao J．2010．REE models for the discrimination of fluids in the formation and evolution of dolomite reservoirs．Oil＆ Gas Geology，31(6):810－818(in Chinese with English abstract)

Huang SJ，Huang Y，Lan YF and Huang KK．2011．A comparative study on strontium isotope composition of dolomites and their coeval seawater in the Late Permian－Early Triassic，NE Sichuan Basin．Acta Petrologica Sinica，27(12):3831－3842(in Chinese with English abstract)

Huang SJ，Huang KK，Lü J and Lan YF．2012．Carbon isotopic composition of Early Triassic marine carbonates，eastern Sichuan Basin，China．Scientia Sinica(Terrae)，42(10):1508－1522(in Chinese)

Kamber BS and Webb GE．2001．The geochemistry of Late Archaean microbial carbonate:Implications for ocean chemistry and continental erosion history．Geochim．Cosmochim．Acta，65(15):2509－2525

Lucia FJ and Major RP．1994．Porosity evolution through hypersaline reflux dolomitization．In:Purser B，Tucker M and Zenger D(eds．)．Dolomite:A Volume in Honor of Dolomieu．Oxford:Blackwell Publishing Ltd．，325 －341

Ma YS，Mu CL，Tan QY and Yu Q．2006a．A discussion on Kaijiang－Liangping ocean trough．Oil＆ Gas Geology，27(3):326－331(in Chinese with English abstract)

Ma YS，Mu CL，Guo XS，Tan QY and Yu Q．2006b．Characteristic and framework of the Changxingian sedientation in the northeastern Sichuan Basin．Geological Review，52(1):25－29(in Chinese with English abstract)

Ma YS，Guo XS，Guo TL，Huang R，Cai XY and Li GX．2007．The Puguang gas field:New giant discovery in the mature Sichuan Basin，Southwest China．AAPG Bulletin，91(5):627－643

Ma YS，Mou CL，Tan QY，Yu Q and Wang RH．2007．Reef－bank features of Permian Changxing Formation and Triassic Feixianguan Formation in the Daxian－Xuanhan area，Sichuan Province，South China and constraint for the reservoirs of natural gas．Earth Science Frontiers，14(1):182 －192(in Chinese with English abstract)

Ma YS，Zhang SC，Guo TL，Zhu GY，Cai XY and Li MW．2008．Petroleum geology of the Puguang sour gas field in the Sichuan Basin，SW China．Marine and Petroleum Geology，25(4－5):357－370

Murray RC．1960．Origin of porosity in carbonate rocks．Journal of Sedimentary Petrology，30(1):59－84

Purser BH，Brown A and Aissaoui DM．1994．Nature，origins and evolution of porosity in dolomites．In:Purser B，Tucker M and Zenger D(eds．)．Dolomites:A Volume in Honor of Dolomieu．Oxford，UK:Blackwell Publishing Ltd．，283 －308

Runnells DD．1969．Diagenesis，chemical sediments，and the mixing of natural waters．Journal of Sedimentary Petrology，39(3):1188－1201

Sun SQ．1995．Dolomite reservoirs:Porosity evolution and reservoir characteristics．AAPG Bulletin，79(2):186－204

Wang RH，Mou CL，Tan QY，Yu Q and Yan JF．2007．Porosity evolution during the diagenesis of the reef shoal dolostones from the Changxing Formation in the Daxian－Xuanhan region， Sichuan．Sedimentary Geology and Tethyan Geology，27(2):9 －13(in Chinese with English abstract)

Wang YG，Zhang J，Liu XG，Xu DH，Shi XR，Song SJ and Wen YC．2005．Sedimentary facies of evaporative carbonate platform of the Feixianguan Formation of Lower Triassic in northeastern Sichuan Basin．Journal of Palaeogeography，7(3):357 －370(in Chinese with English abstract)

Wang YG，Wen YC，Hong HT，Xia ML and Song SY．2006a．Dalong Formation found in Kaijiang－Liangping ocenic trough in the Sichuan Basin．Natural Gas Industry，26(9):32 －36(in Chinese with English abstract)

Wang YG，Wen YC，Hong HT，Xia ML，Zhang J，Song SY and Liu HY．2006b．Petroleum geological characteristics of deep water deposits in Upper Permian－Lower Triassic trough in Sichuan Basin and adjacent areas．Oil＆ Gas Geology，27(5):702－714(in Chinese with English abstract)

Wang YG，Wen YC，Hong HT，Xia ML，Fan Y，Wen L，Kong LX and Wu CH．2009．Carbonate slope facies sedimentary character istics of the Late Permian to Early Triassic in northern Sichuan Basin．Journal of Palaeogeography，11(2):143－156(in Chinese with English abstract)

Warren J．2000．Dolomite:Occurrence，evolution and economically important associations．Earth－Science Reviews，52(1 －3):1 －81

Webb GE and Kamber BS．2000．Rare earth elements in Holocene reefal microbialites:A new shallow seawater proxy．Geochim．Cosmochim．Acta，64(9):1557－1565

Weyl PK．1960．Porosity through dolomitization:Conservation of mass requirements．Journal of Sedimentary Research，30(1):85 －90

Yu ZW．1999．Application of oxygen and carbon isotope in petrogenesis of dolomite．Bulletin of Mineralogy，Petrology and Geochemistry，18(2):103－105(in Chinese with English abstract)

附中文参考文献

党录瑞，郑荣才，郑超，文其兵，陈守春，廖军．2011．川东地区长兴组白云岩储层成因与成岩系统．天然气工业，31(11):17－52

丁振举，刘丛强，姚书振，周宗桂．2000．海底热液系统高温流体的稀土元素组成及其控制因素．地球科学进展，15(6):307－312

冯胜斌，周洪瑞，燕长海，彭翼，袁效奇，贺静．2007．东秦岭二郎坪群硅质岩热水沉积地球化学特征及其地质意义．沉积学报，25(4):564－573

郭金春，马海州，宋恩玉，陈宁．2008．湖泊碳酸盐在过去环境变化研究中的应用．盐湖研究，16(2):66－72

郭彤楼．2011．元坝气田长兴组储层特征与形成主控因素研究．岩石学报，27(8):2381－2391

胡文塇，陈琪，王小林，曹剑．2010．白云岩储层形成演化过程中不同流体作用的稀土元素判别模式．石油与天然气地质，31(6):810－818

黄思静，黄喻，兰叶芳，黄可可．2011．四川盆地东北部晚二叠世－早三叠世白云岩与同期海水锶同位素组成的对比研究．岩石学报，27(12):3831－3842

黄思静，黄可可，吕杰，兰叶芳．2012．早三叠世海水的碳同位素组成与演化——来自四川盆地东部的研究．中国科学(地球科学)，42(10):1508－1522

马永生，牟传龙，谭钦银，余谦．2006a．关于开江－梁平海槽的认识．石油与天然气地质，27(3):326－331

马永生，牟传龙，郭旭升，谭钦银，余谦．2006b．四川盆地东北部长兴期沉积特征与沉积格局．地质论评，52(1):25－29

马永生，牟传龙，谭钦银，余谦，王瑞华．2007．达县－宣汉地区长兴组－飞仙关组礁滩相特征及其对储层的制约．地学前缘，14(1):182－192