岩石学报  2020, Vol. 36 Issue (11): 3477-3492, doi: 10.18654/1000-0569/2020.11.14   PDF    
塔里木盆地古城地区奥陶系鹰山组白云岩特征及孔隙成因
王珊1, 曹颖辉1, 杜德道1, 张亚金2, 刘策1     
1. 中国石油勘探开发研究院, 北京 100083;
2. 大庆油田有限责任公司勘探开发研究院, 大庆 163712
摘要: 奥陶系鹰山组白云岩是塔里木盆地的重点勘探领域,优质白云岩储层是制约其勘探的关键难题。本文在前人研究基础之上,结合多口井岩心薄片观察、同位素、微量元素等地球化学分析,对古城地区奥陶系鹰山组白云岩成因、孔隙成因以及储层主控因素进行了深入研究,认为古城地区鹰山组发育粉晶、粉细晶、细晶、细中晶、中晶、粗晶6种白云岩。粉晶和粉细晶白云岩碳、锶同位素分布范围均与同时期灰岩一致;氧同位素值与同时期白云石一致或偏正;锶元素含量相对较高;稀土元素与同时期灰岩分配模式相同。以上特征表明该类白云岩形成于蒸发海水环境,为准同生海源成因。细晶、细中晶、中晶及粗晶白云岩碳同位素值与同期灰岩总体一致;氧同位素值与正常海相白云石一致或偏负;锶同位素值与灰岩背景值一致或偏高;锶元素含量相对较低;大部分样品Eu正异常。以上特征揭示该类白云岩主要为中浅埋藏成因,并且叠加了早期淡水溶蚀及晚期埋藏-热液改造。孔隙成因分析表明高能缓坡滩是孔隙形成的物质基础、准同生暴露溶蚀是孔隙形成的关键、早期白云石化有利于孔隙的继承和保存,构造破裂和埋藏-热液溶蚀对孔隙起到有利的改造作用。白云岩规模优质储层受控于云化滩和断裂热液溶蚀改造,二者叠合之处为规模优质储层发育区。这一认识为研究区下一步勘探部署提供了重要依据。
关键词: 白云岩    孔隙成因    奥陶系鹰山组    塔里木盆地    古城地区    
Characteristics and pore genesis of dolomite in Ordovician Yingshan Formation in Gucheng area, Tarim Basin
WANG Shan1, CAO YingHui1, DU DeDao1, ZHANG YaJin2, LIU Ce1     
1. PetroChina Research Institute of Petroleum Exploration & Development, Beijing 100083, China;
2. PetroChina DaQing Oilfield Company, Daqing, Heilongjiang 163712, China
Abstract: The dolomite of Ordovician Yingshan Formation is a significant exploration field in Tarim Basin. High-quality dolomite reservoir is the key factor restricting exploration. On the basis of previous studies, combined with core, thin section observation, isotope analysis and trace element analysis, this paper makes a deep study on the dolomitization genesis, pore genesis and controlling factors of dolomite reservoir in Yingshan Formation. The study shows that micrite dolomite, micrite-to fine-grained dolomite, fine-grained dolomite, fine-to medium-grained dolomite, medium-grained dolomite and coarse-grained dolomite are the 6 dolomite types in research area. Geochemical analysis shows that the distribution range of carbon and strontium isotopes of the micrite dolomite and micrite-to fine-grained dolomite are consistent with that of the limestone in the same period; the oxygen isotopic value is consistent with that of dolomite in the same period; the content of strontium element is relatively high; and the distribution pattern of rare earth elements is the same with that of the limestone in the same period, which indicates the dolomitization took place in penecontemporaneous phase. The carbon isotope values of fine-grained dolomite, fine-to medium-grained dolomite, medium-grained dolomite and coarse-grained dolomite are consistent with limestone in the same period, and the oxygen isotope are consistent or negative compared with normal marine dolomites; strontium isotope values are consistent with or higher than those of limestone background values; the content of strontium element is relatively low; most of the samples have positive Eu anomalies, The above characteristics reveal that the dolomitization mainly took place in shallow-medium buried phase, and partially influenced by paracontemporaneous exposure dissolution and hydrothermal. The study of pore genesis shows that the high-energy grainstone shoal is the material basis of pore formation, paracontemporaneous exposure dissolution is the key factor of pore formation, early dolomitization is conducive to the inheritance and preservation of pores, while structural fracture and buried hydrothermal dissolution play a favorable role in pore transformation. Large-scale and high-quality reservoir is controlled by shoal facies and fault hydrothermal solution, the overlap area develops high-quality reservoirs. The above understanding provides a new basis for the next exploration deployment of the study area.
Key words: Dolomite    Pore genesis    Ordovician Yingshan Formation    Tarim Basin    Gucheng area    

塔里木盆地奥陶系是油气勘探的主力层段,目前已发现了多个深层碳酸盐岩大油气田(杨海军等, 2011; 朱光有等, 2012)。古城地区奥陶系鹰山组储层发育,是塔东探区重点勘探领域,但是与塔北和塔中不同(Zhu et al., 2019a, b, c; 杨海军等, 2011; 张水昌等, 2011; 武芳芳等, 2010; 张保涛等, 2015; 孙崇浩等, 2016),该区域鹰山组储层以白云岩为主,前期钻探的GC6、GC8、GC9井相继在该层段获得高产气流,揭示了良好的勘探潜力。但是随后部署的3口探井勘探成果较差,显示了该区的复杂性。基于成功及失利井分析并结合前期研究成果(沈安江等, 2020; Cao et al., 2019; 刘伟等, 2016; 刘卫红等, 2017; 张月巧等, 2007),笔者认为规模优质储层是古城地区奥陶系白云岩勘探成功的关键,深入研究该区白云岩特征及储层成因机理对油气勘探具有重要意义。

前人对古城地区奥陶系鹰山组白云岩储层进行了研究并取得了较为丰富的认识,但尚存在以下问题:(1)对于白云岩成因的探讨相对较少。前人多从单一地球化学指标出发,分析其地球化学特征及地质意义,如郭春涛等(2018)刘策等(2017)通过稀土元素分析,探讨了白云石化流体的性质及成因机制;张振伟等(2016)通过碳、氧同位素分析,证实了早奥陶世晚期存在一期海退。目前尚缺乏利用多种手段对白云岩成因的综合研究。(2)对孔隙的成因认识存在分歧。周波等(2018)认为中等程度的白云石化为孔隙形成的主要原因;张哨楠等(2019)齐井顺等(2016)认为晚期热液溶蚀是孔隙形成的关键;邵红梅等(2014)则认为热液白云石化及热液溶蚀是有效孔隙形成的主要因素。(3)在储层主控因素方面,前人观点相对统一,总体为相控论,认为高能相带为储层形成的首要因素(沈安江等, 2018; 刘伟等, 2016; 张静和罗平, 2010; 胡九珍等, 2009)。然而新钻井情况揭示优质储层并非单纯相控,断裂溶蚀起了极为重要的作用。研究区储层究竟是相控、断控、还是两者综合控制,两种因素各发挥何种作用,需要从更宏观的角度进行综合论述。

鉴于以上原因,本文在前人研究基础之上,通过岩心薄片观察、物性分析、同位素、微量元素等地球化学综合分析,对白云岩成因、孔隙成因以及储层主控因素做了深入探讨,以期对古城地区鹰山组白云岩的勘探部署提供科学依据。

1 地质背景

古城低凸起西南部以塔中Ⅰ号断裂与塔中隆起相邻,东部以上寒武统-中下奥陶统坡折带与东南隆起相接,北部与满西低凸起毗邻(图 1),总体为北西倾的下古生界大型宽缓鼻状构造。古城低凸起形成于加里东中期,定型于加里东末期,海西期受南部挤压作用影响而进一步隆升,印支-喜山期构造活动相对较弱(王招明等, 2014)。

图 1 古城地区构造区划图及奥陶系综合地层柱状图 Fig. 1 Tectonic units and Ordovician stratigraphic histogram of Gucheng area

根据区域地质资料及前人研究成果,塔里木盆地奥陶系经历了碳酸盐岩弱镶边台地-远端变陡的缓坡-混积陆棚的演化过程(Cao et al., 2019)。奥陶系自下而上发育蓬莱坝组、鹰山组、一间房组、吐木休克组和却尔却克组5套地层。鹰山组沉积时期,古城地区为远端变陡的碳酸盐岩缓坡,广泛发育颗粒滩沉积。鹰山组自下而上分为4段,总体厚度为600~800m,鹰四段和鹰三段为鹰下段,厚度为350~550m,为一套颗粒滩相白云岩沉积;鹰二段和鹰一段为鹰上段,厚度为180~280m,为一套开阔台地滩间海相灰岩沉积。

2 白云岩成因 2.1 岩石类型及特征

岩心及薄片观察显示,古城地区鹰山组白云岩主要发育在鹰下段,以结晶白云岩为主,按照晶粒大小,可划分为粉晶、粉细晶、细晶、细中晶、中晶、粗晶白云岩6种类型。

粉晶白云岩直径为0.03~0.1mm,自形-半自形为主,晶体紧密镶嵌,多见纹层构造(图 2a),晶间常被泥质充填,总体致密,孔隙不发育。

图 2 古城地区鹰山组白云岩类型及储集空间特征 (a)纹层状粉晶白云岩,GC17,O1+2y3,6283.30m;(b)纹层状粉细晶白云岩,GC601,O1+2y3,6085.00m;(c)残余颗粒结构细晶白云岩,晶间孔、晶间溶孔发育,GC601,O1+2y3,6075.86m;(d)细晶白云岩,晶间沥青充填,去白云石化发育,GC14,O1+2y3,3-5/15;(e)细中晶白云岩,晶间溶孔发育,局部沥青充填孔隙,GC13,O1+2y3,6268.00m;(f)残余颗粒结构细中晶白云岩,晶间孔顺层发育,局部孔隙被方解石、沥青充填,GC601,O1+2y3,6061.54m;(g)残余颗粒结构中晶白云岩,晶间孔发育,部分晶间孔方解石充填,微裂缝发育,GC601,O1+2y3,6148.10m;(h)中晶白云岩,晶间孔、晶间溶孔发育,部分晶间孔方解石充填,微裂缝沟通孔隙,GC601,O1+2y3,6128.67m;(i)粗晶白云岩,晶间孔、晶间溶孔发育,石英充填,GC9,O1+2y3,6114.20m;(j)粗晶白云岩,渗流粉砂,GC601,O1+2y3,6048.95m;(k)粗晶白云岩中夹泥质条带,GC601,O1+2y3,6083.85m;(l)溶蚀孔洞沿裂缝发育,洞中见白云石、方解石充填,GC601,O1+2y3,10-46/61 Fig. 2 Petrology and reservoir space of dolomite in Yingshan Formation of Gucheng area

粉细晶白云岩由粉晶和细晶白云石组成,可见纹层构造和残余颗粒结构,孔隙不发育(图 2b)。

细晶白云岩晶粒直径多为0.01~0.25mm,自形-半自形为主,常见雾心亮边结构和去白云石化现象。残余颗粒结构明显,孔隙较发育(图 2c),局部见沥青充填(图 2d)。

细中晶白云岩晶形以自形-半自形为主,晶粒紧密镶嵌,局部晶间孔极发育(图 2e),可见孔隙顺层分布(图 2f)。

中晶白云岩晶粒直径主要为0.25~0.5mm,自形-半自形为主,普遍观察到雾心亮边和残余颗粒结构。局部晶间孔极发育,孔中常被方解石、石英等充填(图 2g, h)。

粗晶白云岩晶粒直径大于0.5mm,晶形多样,半自形-他形均有发育,呈紧密镶嵌分布(图 2i, j)。局部层段孔隙极为发育,常见方解石、石英等充填物(图 2i)。

2.2 地球化学特征及成因

本文研究的样品来自古城地区GC601、GC7、GC8、GC9、GC10、GC12、GC13、GC14、GC16共9口井的鹰山组岩心。锶同位素、锶元素及稀土元素分析在核工业北京地质研究院完成,其中锶元素及稀土元素分析仪器为PerkinElmer Elan DCR-e等离子体质谱分析仪,分析标样为GBW07104,分析误差为5%;锶同位素分析仪器为ISOPROBE-T热表面电离质谱仪,分析标样为NBS987,分析误差为±0.00007。碳、氧同位素由核工业北京地质研究院MAT-253气体同位素质谱计和大庆油田有机地球化学实验室ISOPRIM完成,分析标样为GBW04406,分析误差分别为0.1‰和0.2‰。实验结果见表 1表 2

表 1 古城地区奥陶系鹰山组碳酸盐岩碳、氧、锶同位素值 Table 1 The values of carbon, oxygen and strontium isotope for carbonate rocks in Yingshan Formation of Gucheng area

表 2 古城地区奥陶系鹰山组碳酸盐岩稀土元素分析结果(×10-6) Table 2 The REE concentrations and related parameters of carbonate rocks in Yingshan Formation of Gucheng area (×10-6)
2.2.1 碳、氧同位素

碳、氧同位素是判识白云岩成因的重要地球化学标志,其组成主要受控于原岩的碳、氧同位素组成以及成岩流体的盐度、温度及组分等(强子同, 1998)。

白云岩的碳同位素主要继承自其所交代的原岩(Land, 1980),沉积期流体盐度和有机质埋藏率对其影响较为显著,高盐度蒸发环境会导致碳同位素偏正,白云石化过程中若有大量富12C的有机质参与会导致碳同位素明显偏负(黄思静, 2010)。

相比之下,氧同位素受白云石化流体的组分和温度共同控制,更能反映白云岩的成因。高盐度蒸发环境会导致白云岩的氧同位素偏正,较高的形成温度(如深埋环境、热液发育环境)和稀释的白云石化流体(如大气淡水)会导致氧同位素偏负(黄思静, 2010)。

研究区灰岩和各类白云岩的δ13CPDBδ18OPDB值分布情况见表 3。碳同位素分析结果显示(表 3图 3),研究区灰岩δ13CPDB值为-2.40‰~0.68‰,平均值为-1.32‰,与Veizer et al. (1999)所建立的奥陶系海相方解石的δ13CPDB值分布范围(-1.5‰~0.5‰)一致,因此能代表研究区鹰山组沉积时期的海水信息,可将其作为各类白云岩成岩流体性质探讨的背景值。研究区白云岩的δ13CPDB值基本全部处于灰岩碳同位素分布范围之内,仅有3个细中晶和中晶白云岩样品略微偏负,说明白云岩的大部分碳继承自其原岩灰岩,白云石化流体与同期海水有关,局部可能有大气淡水的混入。

表 3 古城地区鹰山组碳酸盐岩碳、氧同位素值 Table 3 The values of carbon and oxygen isotope for carbonate rocks in Yingshan Formation of Gucheng area

图 3 古城地区奥陶系鹰山组碳酸盐岩碳、氧同位素交会图 Fig. 3 The relationship of δ18O vs. δ13C of carbonate rocks in Yingshan Formation of Gucheng area

氧同位素方面,根据Veizer et al. (1999)的研究,寒武-奥陶系海相方解石的δ18OPDB分布范围为-9.50‰~-7.50‰。Land (1985)Mckenzie (1981)研究认为,从相同流体中沉淀出的白云石氧同位素值比方解石高约2‰~3‰;Major et al. (1992)研究则认为该值为1.5‰~3.5‰,本文取二者结论的平均值2.5‰。根据Veizer et al. (1999)的研究数据计算出寒武-奥陶系与方解石同源流体的白云石氧同位素值分布范围约为-7.00‰~-5.00‰。

氧同位素分析结果揭示(表 3图 3),粉晶和粉细晶白云岩δ18OPDB值分布相对集中,分布范围为-8.00‰~3.60‰,平均约为-4.95‰;细晶、细中晶和中晶白云岩δ18OPDB值也较为接近,主要分布在-9.00‰~-6.00‰,平均约为-7.42‰;粗晶白云岩δ18OPDB值分布在-11.2‰~-8.31‰,平均为-9.77‰。可见,随着白云石晶体的增大,δ18OPDB值呈逐渐偏负的趋势,反映了淡水或深埋环境高温的影响。

与奥陶系正常海相白云石氧同位素分布范围相比,31个粉晶和粉细晶白云岩样品中,10个样品与之分布范围一致,20个样品偏正,仅有1个样品略偏负,表明该类白云岩总体形成于蒸发海水或同期正常海水环境,受后期成岩作用的影响较弱。

45个细晶、细中晶和中晶(后文简称细晶-中晶)白云岩样品中,20个样品与同期海相白云石氧同位素分布范围一致,23个样品略偏负(-9.50‰~-7.03‰),δ18OPDB值均>-10.0‰,仅有1个样品偏正(-4.69‰),1个样品明显偏负(-10.08‰)。以上特征说明,一部分细晶-中晶白云岩保留了同沉积海水的性质,一部分受到了其他成岩流体的改造导致氧同位素偏负。由于热液环境下形成的白云岩δ18OPDB值通常 < -10.0‰(Machel and Mountjoy, 1986; Emery and Robinson, 1993),故基本排除热液流体的改造,或者受热液流体改造程度微弱。δ18OPDB值偏负可能是大气淡水改造或深埋环境高温作用所致。前人研究表明(Meyers et al., 1997; 任影等, 2016),大气淡水和海水的混合流体通常具较低δ13CPDB值,并且δ13CPDBδ18OPDB存在较好的相关性。虽然研究区细晶-中晶白云岩δ13CPDB值偏负特征不明显,但与粉晶和粉细晶白云岩相比,δ13CPDB值偏低,且大部分样品的δ13CPDBδ18OPDB相关性较好(图 3)。同时岩石学及地震、测井资料(丁寒生等, 2017)上均有准同生期大气淡水淋滤溶蚀的证据,因此其白云石化过程中应有大气淡水的混入。同时,该类白云岩经历了深埋,故深埋增温效应必定对其有影响。综上,偏轻的氧同位素特征应是大气淡水和埋藏升温叠加作用的结果。

13件粗晶白云岩样品中δ18OPDB值与同期海相白云石相比均偏负,其中4件样品δ18OPDB值< -10.0‰(-11.20‰~-10.05‰),反映了热液流体的改造。9件样品δ18OPDB值>-10.0‰(-9.96‰~-8.31‰),δ13CPDB不具偏负特征,且δ13CPDBδ18OPDB不具明显相关性(图 3),因此理论上该类岩石未经受过大气淡水溶蚀。但是存在另一种可能,即热液或埋藏升温作用导致原始地球化学特征发生蚀变,淡水改造的地球化学信息被掩盖,因此不能完全排除淡水改造的可能性。

方解石脉的δ18OPDB值分布范围为-13.33‰~-10.91‰,平均为-11.94‰,严重偏负。岩心上的产状表明其为最后一期充填,进一步证实研究区白云岩成岩晚期受到了热液的改造。

综上,碳、氧同位素分析表明粉晶和粉细晶白云岩形成于蒸发海水或同期正常海水环境,为准同生海源成因;细晶-中晶白云岩总体为埋藏成因,且早期受到了大气淡水的改造;粗晶白云岩主要反映了埋藏-热液改造的特征。

2.2.2 锶同位素

锶同位素对研究白云岩成岩环境具有重要意义(黄思静等, 2001; 黄思静, 2010; 史忠生等, 2003)。一般认为,87Sr/86Sr值主要受锶来源的控制而不因物理、化学和生物作用发生同位素分馏,因此在一定地质时期内,海水的87Sr/86Sr值是不变的(向芳和王成善, 2001; 黄思静, 2010)。这一特性使锶同位素广泛用于流体来源分析(Roger et al., 2004)。海相碳酸盐岩和海水的锶同位素主要有3个来源:(Veizer et al., 1999; 黄思静, 2010),即壳源、幔源以及古老海相碳酸盐岩重溶。上述三种物质的87Sr/86Sr值分别为0.720±0.005、0.704±0.002及0.708±0.001,三者共同控制了碳酸盐岩锶同位素组成。

奥陶世海水87Sr/86Sr值约为0.7079~0.7092(Veizer et al., 1999; Denison RE et al., 1998),研究区灰岩87Sr/86Sr值为0.708568~0.709197,平均为0.708954,处于同期海水的87Sr/86Sr值范围内,岩心薄片资料上未见明显的后期改造特征,因此该值可作为研究区正常海相碳酸盐岩的锶同位素背景值。

研究区灰岩及各类白云岩87Sr/86Sr值分布情况见表 4。从87Sr/86Sr分布图(图 4)可以直观的得出,粉晶、粉细晶、细晶及细中晶(后文简称粉晶-细中晶)白云岩87Sr/86Sr值十分接近,反映了同一种流体来源,且87Sr/86Sr值均处在灰岩背景值内,说明其继承或保存了同生期海水的锶同位素特征,白云石化介质与同期海源流体有关。

表 4 古城地区鹰山组碳酸盐岩87Sr/86Sr值与Sr元素值分布特征 Table 4 The values of strontium element and strontium isotope for carbonate rocks in Yingshan Formation of Gucheng area

图 4 古城地区奥陶系鹰山组碳酸盐岩87Sr/86Sr分布图 Fig. 4 Diagram showing 87Sr/86Sr values for carbonate rocks in Yingshan Formation of Gucheng area

粗晶白云岩和方解石脉的87Sr/86Sr平均值分别为0.709309和0.709533,总体略高于灰岩背景值,说明盆地内部另有富含放射性87Sr的地层水或热液流体对白云岩进行过改造。根据王坤等(2016)对古城地区流体介质的研究,这种流体可能为二叠纪裂谷拉张背景下沿走滑断裂上涌的中酸性岩浆水。其流经深部碎屑岩疏导层时吸附了大量放射性成因87Sr而导致87Sr/86Sr值较高。

层间含泥粗晶白云岩和层间泥质的87Sr/86Sr平均值分别为0.712366和0.815304,远高于灰岩背景值。结合薄片上观察到的渗滤砂等暴露证据(图 2j),推测受到了壳源锶的注入,应为大气淡水溶蚀淋滤所致。

总体来看,粉晶-细中晶白云岩显示出继承原始灰岩锶同位素的特征,反映了白云石化流体主要与同期海水有关;大部分粗晶白云岩和方解石脉略高于灰岩背景值,说明晚期成岩过程中遭受了深部热液的改造;层间含泥粗晶白云岩和层间泥质具有远高于灰岩背景值的特征,说明其为岩溶面沉积物,研究区成岩早期受到了大气淡水的淋滤溶蚀。

2.2.3 锶元素

碳酸盐岩矿物的Sr主要来源于海水,因而Sr元素含量越高,样品对海水的代表性越好(Yang et al., 1999)。白云石化是Sr消耗的过程(黄思静, 2010),因此Sr含量与白云石化程度有关,一般情况下,早期形成的白云石Sr丰度较晚期形成的白云石要高(Warren, 2000)。蒸发环境中形成的白云石Sr含量可达500×10-6~700×10-6(Bein and Land, 1983);正常海水中形成的白云石的Sr含量理论上为470×10-6~550×10-6(Qing and Mountjoy, 1989),晚期埋藏成因的白云岩Sr含量平均值多在100×10-6以内(Morrow, 1982)。

研究区各类岩石的Sr含量及分布特征见表 4表 5。灰岩的Sr含量值为131.9×10-6~205.1×10-6,平均为160.2×10-6。低于粉晶和粉细晶白云岩的Sr含量,亦显著低于正常海水沉淀方解石的Sr含量(600×10-6)(Baker and Burns, 1985),说明灰岩在新生变形过程中发生了强烈的脱锶作用。数据显示随着白云石晶体的增大,Sr含量逐渐减少,说明从粉晶到中晶,白云石经历的成岩蚀变增强,Sr逐渐流失。同时数据显示研究区各类白云岩和灰岩的Sr含量普遍较低,其中中晶白云岩Sr含量最低,平均为118.3×10-6,但仍大于晚期埋藏成因的白云岩值,说明研究区白云岩并非深埋阶段形成,应主要为中浅埋藏期的产物,晚期埋藏对其有一定的改造作用。层间泥质沉积物Sr含量均值仅为35.20×10-6,较白云岩与灰岩明显偏低,进一步证实了大气淡水的淋滤溶蚀。

表 5 古城地区鹰山组不同岩石类型Sr元素值 Table 5 The values of strontium element for carbonate rocks in Yingshan Formation of Gucheng area
2.2.4 稀土元素

白云岩中稀土元素(REE)组成及地球化学性质主要取决于原岩及白云石化流体的稀土元素(Lottermoser, 1992),常用作成岩流体的来源及形成环境的指示剂。

研究表明,PAAS均一化的海水稀土元素具有轻稀土元素(LREE)亏损、重稀土元素(HREE)富集、Ce亏损、La富集、Gd轻微富集等主要特征(zhang and Nozaki, 1996)。研究区灰岩稀土元素表现为LREE富集、HREE亏损(表 6图 5),这一特征与海水不同。前人研究揭示碳酸盐岩在成岩过程LREE优先富集,且HREE迁移能力强于LREE,造成了岩石中LREE高于HREE的特征(胡忠贵等, 2009)。灰岩δCe平均值为0.78,与全新世古海水的δCe值(0.75)接近(Webb and Kamber, 2000),La具有显著正异常,Gd也表现出轻微的正异常,均与古海水特征相似。因此灰岩总体上代表了沉积时期海水的稀土元素特征。

表 6 古城地区鹰山组碳酸盐岩稀土元素特征 Table 6 The characteristics of rare earth elements of carbonate rocks in Yingshan Formation of Gucheng area

图 5 古城地区奥陶系鹰山组碳酸盐岩PASS标准化(标准化值据Taylor and McLennan, 1985)稀土元素配分图 (a)灰岩、粉晶、粉细晶白云岩稀土元素配分图;(b)灰岩、细晶-粗晶白云岩稀土元素配分图 Fig. 5 PASS-normalized REE distribution of carbonate rocks in Yingshan Formation of Gucheng area

粉晶白云岩和粉细晶白云岩REE配分模式相似,且与灰岩REE配分模式基本相同,均显示轻稀土元素(LREE)富集,重稀土元素(HREE)亏损及Eu、Ce的负异常(表 6图 5a),说明二者的白云石化流体具有同源性,其稀土元素主要继承自原岩灰岩,白云石化流体源于海水。

细晶、细中晶、中晶及粗晶白云岩(后文简称细晶-粗晶白云岩)同样显示LREE富集、HREE亏损的特征(表 6图 5b),说明成岩过程中未受到强烈的陆源物质的影响。这4类白云岩的Eu和Ce具有一定的差异。研究表明稀土元素Eu正异常通常指示高温热液的作用(Frimmel, 2009; 丁振举等, 2000; 胡文瑄等, 2010)。Ce由于其敏感的氧化还原性,常用来指示水体的氧化程度,在强氧化环境中,Ce3+易氧化为Ce4+呈难溶的CeO2迁出而出现负异常(Frimmel, 2009)。研究区22个细晶-粗晶白云岩样品中,5个样品表现为Eu弱负异常(δEu为0.80~0.96,平均0.92)和Ce负异常(δCe值为0.74~0.96,平均为0.83),与灰岩特征相似,反映了对原岩的继承性;17个样品Eu为明显的正异常(δEu值为1.03~3.62,平均为1.62),Ce则有两种情况,11个样品表现为Ce负异常(δCe为0.09~0.94,平均为0.68),6个样品为Ce正异常(δCe为1.21~1.28,平均为1.25),这说明可能有两期性质不同的热液或者同一期热液在不同演化阶段对该类岩石进行了改造。目的层段岩心薄片中可见石英等热液矿物,且闫博等(2018)研究也揭示塔东地区存在多期热液活动。因此,从稀土元素的特征来看,研究区大部分细晶-粗晶白云岩均不同程度受到了热液的改造。

至此,根据上述岩石学及碳、氧、锶同位素、锶元素及稀土元素的分析可得,研究区粉晶和粉细晶白云岩总体形成于蒸发环境,为准同生海源成因;细晶-粗晶白云岩主要为中浅埋藏成因,并且受到了准同生期大气淡水的淋滤及晚期热液流体的改造。

3 孔隙成因及储层主控因素 3.1 孔隙类型及特征

岩心、薄片及成像测井资料表明,研究区鹰山组白云岩储集空间主要有晶间孔、晶间溶孔、溶蚀孔洞及裂缝。

晶间孔  晶间孔为研究区储层段最重要的孔隙类型,孔隙直径较小,约0.05~0.5mm,孔隙边界平直,呈多边形,在细晶-粗晶白云岩中最为发育(图 2c, e-i)。晶间溶孔为晶间孔溶蚀扩大形成,直径相对晶间孔较大,一般约为0.5~1mm之间,通常形态为不规则港湾状(图 2c, h),在岩心上呈针孔状。晶间孔和晶间溶孔通常被白云石、方解石等充填(图 2f-h),充填物中亦见沥青(图 2d)及石英(图 2i)。大部分晶间孔和晶间溶孔具有一定的连通性,渗透率较好。需要说明的是,部分晶间孔或晶间溶孔发育在残余颗粒结构白云岩中,其分布具有规律性,应为粒间孔或者粒内孔演变而来(图 2f, g)。

溶蚀孔洞  直径在2mm以上,形状不规则,具有明显的非均质性,在研究区中晶及粗晶白云岩中最为发育。部分孔洞沿层面发育,推测与原始沉积环境相关,部分孔洞沿裂缝分布,应为裂缝扩溶孔洞,洞壁上多发育多期充填物,一般为白云石、方解石等(图 2l)。

裂缝  在研究区各类岩石中均有发育,是重要的储渗空间,既可作为储集空间,又可作为流体运移的通道。研究区裂缝宽度通常从微米级到毫米级大小不一(图 2g-i, l),部分被白云石、方解石、硅质等充填。

3.2 物性特征

本次研究对部分井的鹰山组白云岩岩心物性进行统计分析,分析结果见表 7。物性数据分析表明,古城地区鹰山组白云岩孔隙度普遍不高,主要分布范围在0.60%~3.0%,大部分样品渗透率也较低,通常小于0.10×10-3μm2。物性数据图表显示(图 6表 7),粗晶、细中晶、中晶白云岩孔隙度相对较高,细晶白云岩次之,粉晶及粉细晶白云岩孔隙度最低,这也同岩心薄片的观察结果一致。由此可得,研究区储层主要发育在细中晶、中晶和粗晶(后文简称细中晶-粗晶)白云岩中,整体为一套中低孔中低渗储层。

表 7 古城地区鹰山组不同类型白云岩孔隙度及渗透率分析结果 Table 7 The values of porosity and permeability of different dolomites in Yingshan Formation of Gucheng area

图 6 古城地区奥陶系鹰山组白云岩孔隙度直方图 Fig. 6 Histogram showing the porosity ranges of dolomite in Yingshan Formation of Gucheng area

从孔渗相关性上来看(图 7),粉细晶和细晶白云岩孔渗相关性较好,说明储层以基质孔为主,相对均质。粉晶、细中晶-粗晶白云岩孔渗相关性较差,部分较低孔隙度的样品具有较高的渗透率,说明裂缝对岩石储集性改造明显。

图 7 古城地区奥陶系鹰山组白云岩孔隙度-渗透率交会图 Fig. 7 The relationship of porosity and permeability of dolomite in Yingshan Formation of Gucheng area
3.3 孔隙成因

基于前文分析并结合沉积构造演化背景总结了古城地区鹰山组储层的白云石化及孔隙演化过程(图 8):早-中奥陶世,研究区为碳酸盐岩缓坡,海水潜流环境下沉积灰岩颗粒滩,由于海水对方解石过饱和,发生胶结作用使原生孔隙被封堵;缓坡环境下,海平面周期性升降,中等盐度的海水在密度差以及海平面变化的驱动下渗透回流,沉积物发生白云石化,形成粉晶-粉细晶白云岩。同时颗粒滩不断暴露出水面,受到准同生期大气淡水淋滤溶蚀,产生一定数量的次生溶孔;浅-中埋藏时期,孔隙水为主要成岩流体,残余海水提供Mg2+,沉积物发生大规模白云石化作用,形成自形-半自形细晶-中晶白云岩,晶间孔发育,连通性较好;深埋藏时期,地层流体为主要成岩流体,随着温度持续增高,发生过度白云石化作用,部分白云石晶体转变为他形中粗晶白云石,早期形成的晶间孔隙被部分封堵;晚奥陶-三叠纪构造运动形成断裂及大量伴生裂缝(冯子辉, 2019),岩石的储集性能提高,多期热液沿裂缝进入成岩体系发生溶蚀,储集空间再次增多;随着溶蚀作用的进行,热液流体的Mg2+对白云石逐渐饱和,基质白云岩再次被改造形成中粗晶他形白云岩,同时孔洞缝中也沉淀一定量的中粗晶鞍形白云石、石英等热液矿物,孔隙减少;之后发生一期或多期方解石充填,孔隙再次减少,形成现今的孔隙发育状态。

图 8 古城地区鹰山组储层白云石化及孔隙演化模式 Fig. 8 Dolomitization and pore evolution model of dolomite reservoir in Yingshan Formation of Gucheng area

基于上述分析,笔者认为古城地区奥陶系鹰山组白云石孔隙形成受以下因素控制。

(1) 高能缓坡滩沉积是孔隙形成的物质基础

根据前文所述,鹰下段白云岩中物性最好的岩石类型为细中晶-粗晶白云岩。薄片观察显示,大部分细中晶-粗晶白云岩中具有明显的残余颗粒结构(图 2f, g),表明其原岩为颗粒灰岩。奥陶系碳酸盐岩缓坡沉积环境下,颗粒灰岩为高能缓坡滩沉积。与之相对,发育在低能环境的纹层状粉晶白云岩、粉细晶白云岩则普遍致密,即使同样经历了早期云化,仍不能发育好储层(图 2a, b)。由此可见,孔隙的发育受控于原始沉积相带,缓坡颗粒滩为鹰山组白云岩孔隙发育的物质基础。

(2) 准同生期暴露溶蚀是孔隙形成的关键

岩心薄片资料显示并非所有白云岩均能发育成好储层。古城地区多口井具有如下特征:致密的滩相白云岩与多孔的滩相白云岩间互发育形成多个向上变浅的旋回,旋回顶部孔隙发育。这种现象应为沉积物受到准同生期周期性暴露溶蚀所致。古城地区鹰三段内幕发育多期与海平面下降相关的暴露,丁寒生等(2017)在鹰山组识别出4个四级层序界面。岩心薄片上顺层理发育的较均匀的基质孔隙(图 2f)、较多渗流粉砂(图 2j)、示底构造及渗流沉积现象(图 2k)均反映了淡水淋滤的特征,前文相关的地化分析亦进一步证实了准同生期的暴露溶蚀。准同生期大气淡水溶蚀可形成良好的晶间孔隙,这些孔隙既是储层孔隙的基础,也为成岩流体提供了运移通道,有利于后期进一步改造。

(3) 早期白云石化作用有利于孔隙的继承和保存

白云石化对研究区孔隙发育究竟起何种作用需要结合岩石类型、孔隙特征等详细分析。古城地区多口钻井揭示,鹰下段灰岩储层总体较差,无论是低能的泥晶灰岩还是高能的亮晶颗粒灰岩,均胶结致密,孔隙基本不发育,优质储层普遍发育在白云岩地层中,这揭示了白云石化作用对孔隙形成的重要性。前文述及研究区储层主要发育在细中晶-粗晶白云岩中,该类白云岩总体为早期中浅埋藏成因,局部受晚期埋藏-热液改造。前人研究表明,早期适度的白云石化有利于原生孔隙的继承,且白云石化使岩石具有更高的强度和脆性,能够更好的抵抗埋藏期的压实作用,这既有利于已有孔隙的保存,也有利于裂缝的形成,同时为后期的热液等成岩流体提供了运移通道和储集空间(Glover, 1968; Hugman and Friedman, 1979)。而晚期埋藏-热液白云石化则主要表现为重结晶及过度白云石化,可部分或者全部封堵孔隙,对储层起破坏作用。综上,白云石化是研究区孔隙大规模存在的前提条件,早期白云石化作用对孔隙的继承和保存最为有利。

(4) 构造破裂及埋藏-热液溶蚀对孔隙起到有利的改造作用

岩心、薄片资料显示,古城地区非组构选择性溶蚀孔洞发育,并且多沿组构选择性孔隙及裂缝分布。这一特征揭示了构造破裂和埋藏-热液溶蚀的叠加改造作用。研究区经历了多期构造运动,产生了大量裂缝,裂缝可沟通原有孔隙,埋藏-热液流体后期可将其溶蚀扩大。可见构造破裂和埋藏-热液溶蚀能极大的改善储层物性,有助于形成裂缝-孔洞型优质储层。

3.4 储层主控因素

研究区成功井及失利井测井资料、地震资料、酸压以及系统试井资料综合分析表明,规模优质储层是古城奥陶系白云岩勘探成功的关键,储层既要质量好,又要成规模,勘探才有可能获得成功。本次研究在明确白云岩孔隙成因的基础之上,重点梳理了研究区3口高产井和2口失利井的钻井及地震资料,从更宏观更综合的角度,对规模优质储层的主控因素进行了初步分析。

GC6、GC8、GC9井3口高产井均发育裂缝-孔洞型规模优质储层(Cao et al., 2019),钻井资料表明3口井均为云化滩相沉积(沈安江等, 2018),储层规模大,地震资料显示断裂发育,呈典型的断溶体特征(图 9),储层质量好。GC601井孔隙发育较普遍,钻井资料表明为云化滩沉积(图 2f),储层规模大,测井解释Ⅱ类储层厚度49.4m,孔隙度为0.10%~3.30%,但地震资料显示该井缺乏断裂溶蚀改造,断溶体特征不明显(图 10a, b),渗透率低,导致储层质量较差从而失利。GC18井地震资料显示断溶体发育(图 10c, d),测井解释Ⅰ类储层厚约110m,孔隙度为3.60%~9.50%,为裂缝-孔洞型优质储层,但是该井断溶体规模较小,井周能量较低物性较差,无大的滩体、缝洞体,储层规模有限,导致该井产水量有限。

图 9 过GC9-GC6-GC8井地震剖面反射特征(a)与结构张量属性剖面反射特征(b) Fig. 9 The reflection characteristics of seismic profile (a) and structural tensor attribute profile (b) through well GC9-GC6-GC8

图 10 过GC601、GC18井地震剖面反射特征(a、c)与结构张量属性剖面反射特征(b、d) Fig. 10 The reflection characteristics of seismic profile (a, c) and structural tensor attribute profile (b, d) through well GC601, GC18

综上可得,古城地区鹰山组规模优质储层受云化滩和断裂热液溶蚀(断溶体)双重控制,其中云化滩控制了储层规模,断裂热液溶蚀控制了储层质量,二者叠合区域应为规模优质储层发育区。这一结论初步回答了沉积相带和断裂溶蚀在储层发育中的作用,也对该区域下一步勘探部署提供了重要依据。

4 结论

(1) 古城地区鹰山组发育粉晶、粉细晶、细晶、细中晶、中晶、粗晶6种白云岩,孔隙主要发育在细中晶-粗晶白云岩中。碳、氧、锶同位素、锶元素及稀土元素分析表明粉晶、粉细晶白云岩总体形成于蒸发海水环境,为准同生海源成因。细晶-粗晶白云岩主要为中浅埋藏成因,并且叠加了早期淡水溶蚀及晚期埋藏-热液改造作用。

(2) 高能缓坡滩是研究区白云岩孔隙形成的物质基础、准同生暴露溶蚀是孔隙形成的关键、早期白云石化有利于孔隙的继承和保存,构造破裂和埋藏-热液溶蚀作用对孔隙起到有利的改造作用。

(3) 鹰山组白云岩规模优质储层受控于云化滩和断裂热液溶蚀改造,云化滩控制储层规模,断裂热液溶蚀控制储层质量,云化滩和断溶体叠合之处为规模优质储层发育区。

参考文献
Baker PA and Burns SJ. 1985. Occurrence and formation of dolomite in organic-rich continental margin sediments. AAPG Bulletin, 69(11): 1917-1930
Bein A and Land LS. 1983. Carbonate sedimentation and diagenesis associated with Mg-Ca-chloride brines:The Permian San Andres Formation in the Texas Panhandle. Journal of Sedimentary Petrology, 53(1): 243-260
Cao YH, Wang S, Zhang YJ, Yang M, Yan L, Zhao YM, Zhang JL, Wang XD, Zhou XX and Wang HJ. 2019. Petroleum geological conditions and exploration potential of Lower Paleozoic carbonate rocks in Gucheng area, Tarim Basin, China. Petroleum Exploration and Development, 46(6): 1165-1181 DOI:10.1016/S1876-3804(19)60271-5
Denison RE, Koepnick RB, Burke WH and Hetherington EA. 1998. Construction of the Cambrian and Ordovician seawater 87Sr/86Sr curve. Chemical Geology, 152(3-4): 325-340 DOI:10.1016/S0009-2541(98)00119-3
Ding HS. 2017. Identification of the dolomite exposure marks in eagle mountain formation and its petroleum geological significance in Gucheng area of East Tarim Basin. Petroleum Geology and Oilfield Development in Daqing, 36(4): 15-19 (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(3): 307-312 (in Chinese with English abstract)
Emery D and Robinson A. 1993. Inorganic Geochemistry:Applications to Petroleum Geology. London: Blackwell Scientific Publications
Feng ZH, Li Q, Zhang YJ, Wang XD and Zhang ZW. 2019. Accumulating conditions and distribution laws of Ordovician hydrocarbon in Gucheng low uplift. Petroleum Geology & Oilfield Development in Daqing, 38(5): 87-93 (in Chinese with English abstract)
Frimmel HE. 2009. Trace element distribution in Neoproterozoic carbonates as palaeoenvironmental indicator. Chemical Geology, 258(3/4): 338-353
Glover JE. 1968. Significance of stylolites in dolomitic limestones. Nature, 217(5131): 835-836 DOI:10.1038/217835a0
Guo CT, Li RY and Chen SM. 2018. Rare earth element geochemistry and genetic model of dolomite of Yingshan Formation in Gucheng area, Tarim Basin. Journal of Jilin University (Earth Science Edition), 48(4): 1121-1134 (in Chinese with English abstract)
Hu JZ, Liu SG, Ran QG, Sun W and Zhang CJ. 2009. Diagenetic characteristics and their effect on the formation of good-quality reservoirs of the Cambrian system to Lower Ordovician in the east of Tarim Basin, Xinjiang, China. Journal of Chengdu University of Technology (Science & Technology Edition), 36(2): 138-146 (in Chinese with English abstract)
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)
Hu ZG, Zheng RC, Hu JZ, Wen HG, Li Y, Wen QB and Xu FB. 2009. Geochemical characteristics of rare earth elements of Huanglong Formation dolomites reservoirs in eastern Sichuan-northern Chongqing area. Acta Geologica Sinica, 83(6): 782-790 (in Chinese with English abstract)
Huang SJ, Shi H, Liu J and Shen LC. 2001. Progress in strontium isotope stratigraphy. Advance in Earth Sciences, 16(2): 194-200 (in Chinese with English abstract)
Huang SJ. 2010. Carbonate Diagenesis. Beijing: Geological Publishing House, 121-122 (in Chinese)
Hugman Ⅲ RHH and Friedman M. 1979. Effects of texture and composition on mechanical behavior of experimentally deformed carbonate rocks. AAPG Bulletin, 63(9): 1478-1489
Land LS. 1980. The isotopic and trace element geochemistry of dolomite:The state of the art. SEPM Special Publication, 28: 87-110
Land LS. 1985. The origin of massive dolomite. Journal of Geological Education, 33(2): 112-125 DOI:10.5408/0022-1368-33.2.112
Liu C, Zhang YJ, Li HH, Cao YH, Zhou B and Yang M. 2017. Nature of dolomitizing fluids of Middle-Low Ordovician dolomites in the Gucheng area, Tarim Basin:Evidence from rare earth elements geochemistry. Bulletin of Mineralogy, Petrology and Geochemistry, 36(4): 602-610 (in Chinese with English abstract)
Liu W, Shen AJ, Liu GD, Zheng XP, Chen YN and Zhang Y. 2016. Characteristics and exploration domains of Lower Paleozoic carbonate reservoirs in eastern Tarim Basin. Marine Origin Petroleum Geology, 21(2): 1-12 (in Chinese with English abstract)
Liu WH, Gao XZ, Ye XL, Ran QG and Cheng HG. 2017. Conditions for gas accumulation and potential areas in Lower Paleozoic eastern Tarim Basin. Marine Origin Petroleum Geology, 22(4): 61-68 (in Chinese with English abstract)
Lottermoser BG. 1992. Rare Earth elements and hydrothermal ore formation processes. Ore Geology Reviews, 7(1): 25-41
Machel HG and Mountjoy EW. 1986. Chemistry and environments of dolomitization:A reappraisal. Earth Science Reviews, 23(3): 175-222 DOI:10.1016/0012-8252(86)90017-6
Major RP, Lloyd RM and Lucia FJ. 1992. Oxygen isotope composition of Holocene dolomite formed in a humid hypersaline setting. Geology, 20(7): 586-588 DOI:10.1130/0091-7613(1992)020<0586:OICOHD>2.3.CO;2
McKenzie JA. 1981. Holocene dolomitization of calcium carbonate sediments from the coastal sabkhas of Abu Dhabi, U.A.E.:A stable isotope study. The Journal of Geology, 89(2): 185-198 DOI:10.1086/628579
Meyers WJ, Lu FH and Achariah JK. 1997. Dolomitization by mixed evaporative brines and freshwater, Upper Miocene carbonates, Nijar, Spain. Journal of Sedimentary Research, 67(5): 898-912
Morrow DW. 1982. Diagenesis 1. Dolomite-Part 1:The chemistry of dolomitization and dolomite precipitation. Geoscience Canada, 9(1): 5-13
Qi JS, Zheng XP, Huang SW, Zhang Y and Shao GM. 2016. Characteristics of Yingshan dolomite reservoir in Gucheng low swell of Tarim Basin with favorable area prediction. Xinjiang Petroleum Geology, 37(1): 7-12 (in Chinese with English abstract)
Qiang ZT. 1998. Carbonate Reservoir Geology. Dongying: China University of Petroleum Press (in Chinese)
Qing HR and Mountjoy EW. 1989. Multistage dolomitization in Rainbow buildups, Middle Devonian Keg River Formation, Alberta, Canada. Journal of Sedimentary Research, 59(1): 114-126
Ren Y, Zhong DK, Gao CL, Yang XQ, Xie R, Li ZP, Deng MX and Zhou YC. 2016. Geochemical characteristics, genesis and hydrocarbon significance of dolomite in the Cambrian Longwangmiao Formation, eastern Sichuan Basin. Acta Petrolei Sinica, 37(9): 1102-1115 (in Chinese with English abstract)
Roger JB, Gregg CO and Guo QG. 2004. Strontium isotopic signatures of oil-field waters:Applications for reservoir characterization. AAPG Bulletin, 88(2): 1677-1704
Shao HM, Feng ZH, Wang C, Lu Xi, Wang YK and Gao B. 2014. Types and geneses of the dolomite reservoirs in Lower Paleozoic of Gucheng area of Tarim Basin. Petroleum Geology and Oilfield Development in Daqing, 33(5): 111-116 (in Chinese with English abstract)
Shen AJ, Fu XT, Zhang Y, Zheng XP, Liu W, Shao GM and Cao YQ. 2018. A study of source rocks & carbonate reservoirs and its implication on exploration plays from Sinian to Lower Paleozoic in the east of Tarim Basin, Northwest China. Natural Gas Geoscience, 29(1): 1-16 (in Chinese with English abstract)
Shen AJ, Zhang Y, Feng ZH, Zheng XP, Zhu M, Shao GM, Song X, Zhu KD and She M. 2020. Geological understandings and exploration prospects of carbonate reservoirs in Gucheng area, Tadong, Tarim Basin. China Petroleum Exploration, 25(3): 96-106 (in Chinese with English abstract)
Shi ZS, Chen KY, He HJ, Shi J, Liu BJ and Liu G. 2003. Application of Sr isotope to sedimentary environment analysis. Oil & Gas Geology, 24(2): 187-190 (in Chinese with English abstract)
Sun CH, Zhu GY, Zheng DM, Zhao KZ and Hu FJ. 2016. Characteristics and controlling factors of fracture-cavity carbonate reservoirs in the Halahatang Area, Tarim Basin. Bulletin of Mineralogy, Petrology and Geochemistry, 35(5): 1028-1036 (in Chinese with English abstract)
Taylor SR and McLennan SM. 1985. The Continental Crust:Its Composition and Evolution. Oxford: Blackwell, 1-312
Veizer J, Ala D, Azmy K, Bruckschen P, Buhl D, Bruhn F, Carden GAF, Diener A, Ebneth S, Godderis Y, Jasper T, Korte C, Pawellek F, Podlaha OG and Strauss H. 1999. 87Sr/86Sr, δ13C and δ18O evolution of Phanerozoic seawater. Chemical Geology, 161(1-3): 59-88 DOI:10.1016/S0009-2541(99)00081-9
Wang K, Hu SY, Hu ZY, Liu W, Huang QY, Shi SY, Ma K and Li M. 2016. Cambrian hydrothermal action in Gucheng area, Tarim Basin and its influences on reservoir development. Acta Petrolei Sinica, 37(4): 439-453 (in Chinese with English abstract)
Wang ZM, Yang HJ, Qi YM, Cheng YQ and Xu YL. 2014. Ordovician gas exploration breakthrough in the Gucheng lower uplift of the Tarim Basin and its enlightenment. Natural Gas Industry, 34(1): 1-9 (in Chinese with English abstract)
Warren J. 2000. Dolomite:Occurrence, evolution and economically important associations. Earth-Science Reviews, 52(1-3): 1-81 DOI:10.1016/S0012-8252(00)00022-2
Webb CE and Kamber BS. 2000. Rare earth elements in Holocene reefal microbialitites:A new shallow seawater proxy. Geochimica et Cosmochimica Acta, 64: 1557-1565 DOI:10.1016/S0016-7037(99)00400-7
Wu FF, Zhu GY, Zhang SC, Jin Q, Han JF, Chen ZW and Meng SC. 2010. Ordovician hydrocarbon accumulation and main controlling factors in the western margin of Tazhong No. 1 slope-break zone, Tarim Basin. Geological Review, 56(3): 339-348 (in Chinese with English abstract)
Xiang F and Wang CS. 2001. New progress in the application of strontium isotopes in sedimentology. Geology-Geochemistry, 29(1): 79-82 (in Chinese with English abstract)
Yan B. 2018. Mechanism of geological activities of eastern Tarim basin hydrothermal fluids and its significance in reservoir transformation. Journal of Southwest Petroleum University (Science & Technology Edition), 40(4): 17-28 (in Chinese with English abstract)
Yang HJ, Zhu GY, Han JF, Wu FF, Ji YG, Su J, Zhang HZ and Wang Y. 2011. Conditions and mechanism of hydrocarbon accumulation in large reef-bank karst oil/gas fields of Tazhong area, Tarim Basin. Acta Petrologica Sinica, 27(6): 1865-1885 (in Chinese with English abstract)
Yang JD, Sun WG, Wang ZZ, Xue YS and Tao XC. 1999. Variations in Sr and C isotopes and Ce anomalies in successions from China:Evidence for the oxygenation of Neoproterozoic seawater?. Precambrian Research, 93(2-3): 215-233 DOI:10.1016/S0301-9268(98)00092-8
Zhang BT, Yu BS, Zhu GY, Wang Y, Su J, Wang K and Liu XW. 2015. Research on hydrocarbon distribution regularity and main accumulation controlling factors of Tazhong northern slope hydrocarbon-rich region. Earth Science Frontiers, 22(1): 271-280 (in Chinese with English abstract)
Zhang J and Luo P. 2010. Genesis of Ordovician matrix-porosity dolomite reservoir in the Tarim Basin. Petroleum Geology and Experiment, 32(5): 470-474, 503 (in Chinese with English abstract)
Zhang J and Nozaki Y. 1996. Rare earth elements and yttrium in seawater:ICP-MS determinations in the East Caroline, Coral Sea, and South Fiji basins of the western South Pacific Ocean. Geochimica et Cosmochimica Acta, 60: 4631-4644 DOI:10.1016/S0016-7037(96)00276-1
Zhang SC, Zhu GY, Yang HJ, Su J, Yang DB, Zhu YF, Zhang B and Cui J. 2011. The phases of Ordovician hydrocarbon and their origin in the Tabei uplift, Tarim Basin. Acta Petrologica Sinica, 27(8): 2447-2460 (in Chinese with English abstract)
Zhang SN, Huang BW, Sui H, Ye N and Li YT. 2019. The reservoir characteristics and the pore evolution of Yingshan Formation in Gucheng region, Tarim Basin. Journal of Southwest Petroleum University (Science & Technology Edition), 41(1): 33-46 (in Chinese with English abstract)
Zhang YQ, Jia JD, Jin JQ and Liu QY. 2007. Characteristics of Cambrian-Ordovician sedimentary facies in Tadong region and its sedimentary model. Natural Gas Geoscience, 18(2): 229-234 (in Chinese with English abstract)
Zhang ZW, Ding HS and Zhang YJ. 2016. Carbon and oxygen stable isotope features of the Lower Ordovician in Gucheng region of Tarim basin. Marine Geology & Quaternary Geology, 36(2): 59-64 (in Chinese with English abstract)
Zhou B, Cao YH, Qi JS, Huang SW, Liu C, Jia JH and Chen XY. 2018. Development mechanism of Lower Ordovician reservoir in the Gucheng area, Tarim Basin, China. Natural Gas Geoscience, 29(6): 773-783 (in Chinese with English abstract)
Zhu GY, Yang HJ, Su J, He K, Han JF, Gu LJ, Zhu YF and Zhang BT. 2012. True exploration potential of marine oils in the Tarim Basin. Acta Petrologica Sinica, 28(4): 1333-1347 (in Chinese with English abstract)
Zhu GY, Milkov AV, Zhang ZY, Sun CH, Zhou XX, Chen FR, Han JF and Zhu YF. 2019a. Formation and preservation of a giant petroleum accumulation in superdeep carbonate reservoirs in the southern Halahatang oil field area, Tarim Basin, China. AAPG Bulletin, 103(7): 1703-1743 DOI:10.1306/11211817132
Zhu GY, Zhang Y, Zhou XX, Zhang ZY, Du DD, Shi SB, Li TT, Chen WY and Han JF. 2019b. TSR, deep oil cracking and exploration potential in the Hetianhe gas field, Tarim Basin, China. Fuel, 236: 1078-1092 DOI:10.1016/j.fuel.2018.08.119
Zhu GY, Zhang ZY, Milkov AV, Zhou XX, Yang HJ and Han JF. 2019c. Diamondoids as tracers of late gas charge in oil reservoirs:Example from the Tazhong area, Tarim Basin, China. Fuel, 253: 998-1017 DOI:10.1016/j.fuel.2019.05.030
丁寒生. 2017. 塔东古城地区鹰山组白云岩暴露标志识别及其油气地质意义. 大庆石油地质与开发, 36(4): 15-19.
丁振举, 刘丛强, 姚书振, 周宗桂. 2000. 海底热液系统高温流体的稀土元素组成及其控制因素. 地球科学进展, 15(3): 307-312.
冯子辉, 李强, 张亚金, 王显东, 张振伟. 2019. 古城低凸起奥陶系油气成藏条件与分布规律. 大庆石油地质与开发, 38(5): 87-93.
郭春涛, 李如一, 陈树民. 2018. 塔里木盆地古城地区鹰山组白云岩稀土元素地球化学特征及成因. 吉林大学学报(地球科学版), 48(4): 1121-1134.
胡九珍, 刘树根, 冉启贵, 孙玮, 张长俊. 2009. 塔东地区寒武系-下奥陶统成岩作用特征及对优质储层形成的影响. 成都理工大学学报(自然科学版), 36(2): 138-146.
胡文瑄, 陈琪, 王小林, 曹剑. 2010. 白云岩储层形成演化过程中不同流体作用的稀土元素判别模式. 石油与天然气地质, 31(6): 810-818.
胡忠贵, 郑荣才, 胡九珍, 文华国, 李瑜, 文其兵, 徐发波. 2009. 川东-渝北地区黄龙组白云岩储层稀土元素地球化学特征. 地质学报, 83(6): 782-790.
黄思静, 石和, 刘洁, 沈立成. 2001. 锶同位素地层学研究进展. 地球科学进展, 16(2): 194-200.
黄思静. 2010. 碳酸盐岩的成岩作用. 北京: 地质出版社, 121-122.
刘策, 张义杰, 李洪辉, 曹颖辉, 周波, 杨敏. 2017. 塔里木盆地古城地区中下奥陶统白云化流体性质厘定——来自稀土元素的证据. 矿物岩石地球化学通报, 36(4): 602-610.
刘伟, 沈安江, 柳广弟, 郑兴平, 陈娅娜, 张友. 2016. 塔里木盆地塔东地区下古生界碳酸盐岩储层特征与勘探领域. 海相油气地质, 21(2): 1-12.
刘卫红, 高先志, 叶信林, 冉启贵, 程宏岗. 2017. 塔东地区下古生界天然气成藏主控因素及勘探领域. 海相油气地质, 22(4): 61-68.
齐井顺, 郑兴平, 黄世伟, 张友, 邵冠铭. 2016. 塔里木盆地古城低凸起鹰山组白云岩储集层特征及有利区预测. 新疆石油地质, 37(1): 7-12.
强子同. 1998. 碳酸盐岩储层地质学. 东营: 中国石油大学出版社.
任影, 钟大康, 高崇龙, 杨雪琪, 谢瑞, 李卓沛, 邓闵心, 周勇成. 2016. 川东寒武系龙王庙组白云岩地球化学特征、成因及油气意义. 石油学报, 37(9): 1102-1115.
邵红梅, 冯子辉, 王成, 卢曦, 王彦凯, 高波. 2014. 塔里木盆地古城地区下古生界白云岩储层类型及成因. 大庆石油地质与开发, 33(5): 111-116.
沈安江, 付小东, 张友, 郑兴平, 刘伟, 邵冠铭, 曹彦清. 2018. 塔里木盆地塔东地区震旦系-下古生界碳酸盐岩油气生储条件与勘探领域. 天然气地球科学, 29(1): 1-16.
沈安江, 张友, 冯子辉, 郑兴平, 朱茂, 邵冠铭, 宋叙, 朱可丹, 佘敏. 2020. 塔东古城地区碳酸盐岩储层地质认识与勘探领域. 中国石油勘探, 25(3): 96-106.
史忠生, 陈开运, 何胡军, 史军, 刘保军, 刘刚. 2003. 锶同位素在沉积环境分析方面的应用. 石油与天然气地质, 24(2): 187-190.
孙崇浩, 朱光有, 郑多明, 赵宽志, 胡方杰. 2016. 塔里木盆地哈拉哈塘地区超深碳酸盐岩缝洞型储集层特征与控制因素. 矿物岩石地球化学通报, 35(5): 1028-1036.
王坤, 胡素云, 胡再元, 刘伟, 黄擎宇, 石书缘, 马奎, 李梅. 2016. 塔里木盆地古城地区寒武系热液作用及其对储层发育的影响. 石油学报, 37(4): 439-453.
王招明, 杨海军, 齐英敏, 陈永权, 徐彦龙. 2014. 塔里木盆地古城地区奥陶系天然气勘探重大突破及其启示. 天然气工业, 34(1): 1-9.
武芳芳, 朱光有, 张水昌, 金强, 韩剑发, 陈志伟, 孟书翠. 2010. 塔里木盆地塔中Ⅰ号断裂带西缘奥陶系油气成藏与主控因素研究. 地质论评, 56(3): 339-348.
向芳, 王成善. 2001. 锶同位素在沉积学中的应用新进展. 地质地球化学, 29(1): 79-82.
闫博. 2018. 塔东热液地质作用机制及对储层的改造意义. 西南石油大学学报(自然科学版), 40(4): 17-28.
杨海军, 朱光有, 韩剑发, 武芳芳, 吉云刚, 苏劲, 张海祖, 王宇. 2011. 塔里木盆地塔中礁滩体大油气田成藏条件与成藏机制研究. 岩石学报, 27(6): 1865-1885.
张保涛, 于炳松, 朱光有, 王宇, 苏劲, 王凯, 刘星旺. 2015. 塔中北斜坡富油气区油气分布规律与富集主控因素研究. 地学前缘, 22(1): 271-280.
张静, 罗平. 2010. 塔里木盆地奥陶系孔隙型白云岩储层成因. 石油实验地质, 32(5): 470-474, 503.
张水昌, 朱光有, 杨海军, 苏劲, 杨德彬, 朱永峰, 张斌, 崔杰. 2011. 塔里木盆地北部奥陶系油气相态及其成因分析. 岩石学报, 27(8): 2447-2460.
张哨楠, 黄柏文, 隋欢, 叶宁, 李映涛. 2019. 古城地区鹰山组储层特征及成岩孔隙演化. 西南石油大学学报(自然科学版), 41(1): 33-46.
张月巧, 贾进斗, 靳久强, 刘全有. 2007. 塔东地区寒武-奥陶系沉积相与沉积演化模式. 天然气地球科学, 18(2): 229-234.
张振伟, 丁寒生, 张亚金. 2016. 塔里木盆地古城地区下奥陶统碳酸盐岩碳氧同位素特征. 海洋地质与第四纪地质, 36(2): 59-64.
周波, 曹颖辉, 齐井顺, 黄世伟, 刘策, 贾进华, 陈秀艳. 2018. 塔里木盆地古城地区下奥陶统储层发育机制. 天然气地球科学, 29(6): 773-783.
朱光有, 杨海军, 苏劲, 何坤, 韩剑发, 顾礼敬, 朱永峰, 张保涛. 2012. 塔里木盆地海相石油的真实勘探潜力. 岩石学报, 28(4): 1333-1347.