激光原位U-Pb同位素定年技术及其在碳酸盐岩成岩-孔隙演化中的应用
沈安江1,2, 胡安平1,2, 程婷3,4, 梁峰1,2, 潘文庆5, 俸月星3, 赵建新3
1. 中国石油天然气集团有限公司碳酸盐岩储层重点实验室,杭州 310023
2. 中国石油杭州地质研究院,杭州 310023
3. Radiogenic Isotope Facility, School of Earth and Environmental Sciences, The University of Queensland, Brisbane, QLD 4072, Australia
4. 中国地质科学院北京离子探针中心,北京 100037
5. 中国石油塔里木油田公司勘探开发研究院,新疆库尔勒 841000
联系作者简介:胡安平(1982-),女,浙江天台人,博士,中国石油杭州地质研究院高级工程师,主要从事碳酸盐岩储集层研究与地球化学实验技术研发工作。地址:浙江省杭州市西湖区西溪路920号,中国石油杭州地质研究院,邮政编码:310023。E-mail: huap_hz@petrochina.com.cn

第一作者简介:沈安江(1965-),男,浙江兰溪人,博士,中国石油杭州地质研究院教授级高级工程师,主要从事碳酸盐岩沉积储集层研究。地址:浙江省杭州市西湖区西溪路920号,中国石油杭州地质研究院,邮政编码:310023。E-mail: shenaj_hz@petrochina.com.cn

摘要

为解决溶液法难以实现的古老海相碳酸盐岩取样和超低U、Pb含量样品测年的难题,通过激光剥蚀进样系统与多接收电感耦合等离子体质谱仪的连用和对年龄为209.8 Ma实验室工作标样的开发和标定,建立了适用于古老海相碳酸盐岩的激光原位U-Pb同位素定年技术,并应用于四川盆地震旦系灯影组成岩-孔隙演化研究。通过充填孔洞、孔隙和裂缝中不同期次白云石胶结物的测年,指出灯影组白云岩储集层的埋藏成岩过程主要是原生孔隙和表生溶蚀孔洞逐渐被充填的过程。孔洞的充填作用发生在早加里东、晚海西—印支、燕山—喜马拉雅期3个阶段,孔隙的充填作用主要发生于早加里东期,未被胶结物充填的残留孔洞、孔隙和裂缝构成了主要储集空间,明确了四川盆地灯影组白云岩储集层的成岩-孔隙演化史。这些认识与该地区的构造-埋藏史、盆地热史具有很高的吻合度,说明测年数据的可靠性和激光原位U-Pb同位素定年技术的有效性,为古老海相碳酸盐岩成岩-孔隙演化研究和油气运移前有效孔隙评价提供了新的方法。图5表4参46

关键词: 激光原位剥蚀; U-Pb同位素定年; 四川盆地; 震旦系灯影组; 成岩-孔隙演化; 碳酸盐岩
中图分类号:TE122.1 文献标志码:A 文章编号:1000-0747(2019)06-1062-13
Laser ablation in situ U-Pb dating and its application to diagenesis-porosity evolution of carbonate reservoirs
SHEN Anjiang1,2, HU Anping1,2, CHENG Ting3,4, LIANG Feng1,2, PAN Wenqing5, FENG Yuexing3, ZHAO Jianxin3
1. Key Laboratory of Carbonate reservoirs, CNPC, Xihu District, Hangzhou 310023, China
2. PetroChina Hangzhou Research Institute of Geology (HIPG), Xihu District, Hangzhou 310023, China
3. Radiogenic Isotope Facility, School of Earth and Environmental Sciences, The University of Queensland, Brisbane, QLD 4072, Australia
4. Beijing SHRIMP Center, Institute of Geology, Chinese Academy of Geological Sciences (CAGS), Haidian District, Beijing 100037, China
5. Research Institute of Exploration and Development, PetroChina Tarim Oilfield, Korla 841000, China
Abstract

Through the development and calibration of a standard carbonate sample with an age of 209.8 Ma using a newly-developed Laser Ablation (LA) Multi-Collector Inductively Coupled Plasma Mass Spectrometry (MC-ICP-MS) technique, we successfully overcome the difficulty in sampling and dating ultra-low U-Pb ancient marine carbonates, which was previously untenable by isotope dilution (ID) methods. We developed the LA-MC-ICP-MS in situ U-Pb dating technique for ancient marine carbonates for the study of diagenesis-porosity evolution history in Sinian Dengying Formation, Sichuan Basin. By systematically dating of dolomitic cements from vugs, matrix pores and fractures, we found that the burial and diagenetic process of dolomite reservoirs in Sinian Dengying Formation was characterized by progressive filling-up of primary pores and dissolution vugs formed by epigenic karstification. The filling of vugs happened in three stages, early Caledonian, late Hercynian-Indosinian and Yanshanian-Himalayan, while the filling of matrix pores mainly took place in early Caledonian. The unfilled residual vugs, pores and fractures constitute the main storage space in the reservoir. Based on the above knowledge, we established the diagenesis-porosity evolution history of the dolomite reservoir in Sinian Dengying Formation, Sichuan Basin. These findings are highly consistent with the tectonic-burial and basin thermal histories of the study area. Our study demonstrates the reliability of in situ U-Pb dating for providing an effective way for the investigation of diagenesis-porosity evolution history and evaluation of porosity in ancient marine carbonate reservoirs before hydrocarbon migration.

Keyword: laser ablation; in-situ U-Pb dating; Sichuan Basin; Sinian Dengying Formation; diagenesis-porosity evolution; carbonates
0 引言

碳酸盐岩在全球油气勘探中占有重要地位, 近50%的油气资源分布在碳酸盐岩中。储集层成因和分布预测是碳酸盐岩油气勘探面临的关键问题之一, 前人在这方面做了大量研究工作, 取得了很多认识。碳酸盐岩孔隙成因问题, Kerans等[1]和Moore等[2]认为早成岩和高化学活动性导致碳酸盐岩以次生孔隙为主, 而且主要形成于埋藏溶蚀作用, James和Choquette[3]、Lucia等[4]则提出碳酸盐岩以沉积原生孔为主。白云石化对孔隙的贡献问题, 由于白云岩储集层的物性普遍好于灰岩, 储集层主要发育于白云岩中, 故大多数学者认为白云石化对孔隙的发育有重要贡献[5], 而Lucia等[6]、Purser等[7]则认为只有CO32-来源很受局限的成岩环境, 白云石化作用才能导致孔隙增加。热液作用对孔隙贡献问题, Davis等[8]认为受构造控制的热液活动导致白云岩储集层的发育。层序格架中储集层分布问题, Moore[2]指出台缘带和蒸发台地是碳酸盐岩储集层最有利的发育相带。

中国海相碳酸盐岩具年代古老、埋藏深和经历多期成岩叠加改造的特点[9], 成储和成藏历程均非常复杂; 勘探实践证实, 优质储集层发育段并不总是油气层段, 也有可能是水层或干层, 除缺烃源外, 还和孔隙发育时间与油气运移时间不匹配有关。这就需要开展储集层成岩-孔隙演化研究, 评价油气运移前的有效孔隙, 碳酸盐岩成岩矿物绝对年龄的确定是储集层成岩-孔隙演化史恢复的关键。碳酸盐U-Th溶液法可以进行0~500 000年碳酸盐矿物的绝对年龄测定, 精度可达1~2年[10, 11], 碳酸盐U-Pb溶液法同位素定年在中新生代年轻的孔洞和洞穴充填物定年研究中也有不少报道[12, 13, 14, 15, 16], 但由于U含量普遍偏低、缺乏合适的标样、难以钻取足够量的粉末样品(6~8个平行样品, 每份200 mg)等问题, 导致古老海相碳酸盐岩溶液法定年费时, 适合定年的样品不多, 测试成功率低, 无法推广。

本文通过对碳酸盐矿物的激光MC-ICPMS原位U-Pb同位素定年技术的建立和新的实验室工作标样的开发和标定, 解决了同位素稀释溶液法在采样和测试方面难以解决的一系列技术难题, 成功建立了适用于古老海相碳酸盐岩的同位素定年新技术。该技术应用于四川盆地震旦系灯影组白云岩成岩-孔隙演化研究, 取得的认识与工区的构造-埋藏史、盆地热演化史非常吻合, 说明测年数据的可靠性和激光原位U-Pb同位素定年技术的有效性, 为古老海相碳酸盐岩成岩-孔隙演化研究和孔隙评价提供了技术支撑。

1 地质背景和样品描述
1.1 区域地质背景

四川盆地震旦系灯影组自下而上可划分为灯一、灯二、灯三和灯四段[17], 以台地相沉积为主[18, 19]。灯一段沉积是晚震旦世早期海侵的产物, 主要为浅灰— 深灰色层状泥粉晶白云岩, 夹砂屑和藻屑白云岩, 与下震旦统陡山陀组呈整合或假整合接触, 厚300~450 m。灯二段沉积由早期至晚期, 由浅水台地藻纹层和藻砂屑白云岩(重结晶后呈粉— 细晶白云岩)转变为膏云岩及膏盐岩沉积, 海水盐度的增加有利于微生物的繁殖, 发育葡萄花边状构造, 残留孔洞发育, 受桐湾运动Ⅰ 幕影响, 使灯二段抬升遭受风化剥蚀, 形成近南北向展布的侵蚀谷[20], 与上覆地层呈假整合接触, 厚400~800 m。灯三段沉积早期发育海侵相的泥岩, 向南西方向泥岩逐渐减薄消失, 晚期发育浅水台地泥粉晶白云岩和颗粒滩沉积。灯四段沉积期是台内裂陷发育的鼎盛期[21], 台缘和台内微生物丘滩复合体发育, 岩性主要为藻纹层或藻叠层白云岩, 基质孔和孔洞发育, 受桐湾运动Ⅱ 幕影响, 使灯四段遭受不同程度的淋滤和剥蚀, 与上覆地层呈假整合接触, 残留厚度30~400 m。灯影组构造-岩相古地理特征对成储有重要的控制作用。

四川盆地加里东古隆起自灯影组沉积以来, 经历了5期构造演化阶段[22]。①加里东旋回早期构造演化阶段:发生桐湾Ⅰ 幕和桐湾Ⅱ 幕两期构造运动, 分别导致灯二段和灯四段的抬升和剥蚀。②加里东旋回中晚期构造演化阶段:寒武纪— 奥陶纪发生了3次超覆沉积与3次隆升剥蚀, 分别为兴凯运动、郁南运动和都匀运动; 志留纪末期的广西运动导致川中加里东古隆起整体抬升剥蚀, 并与二叠系呈平行不整合接触。③海西期构造演化阶段:上扬子区泥盆系— 石炭系整体隆升剥蚀, 石炭纪末受云南运动影响, 川中进一步遭受剥蚀; 二叠纪四川盆地主体处于沉降沉积期, 东吴运动导致二叠系茅口组遭受剥蚀。④印支— 燕山期构造演化阶段:中晚三叠世之交的印支运动完成了四川盆地由海相向陆相沉积的转换, 中下三叠统遭受不同程度的剥蚀。⑤喜马拉雅期构造演化阶段:古隆起东段的高石梯— 龙女寺相对稳定, 埋深大, 而古隆起西段的乐山— 资阳强烈褶皱, 埋深小。构造演化对四川盆地灯影组白云岩储集层改造、油气成藏及演化具重要的控制作用[19, 23]

1.2 样品特征与产状

灯影组白云岩是四川盆地重要的油气储集层, 主要发育于灯二段和灯四段, 累计厚20~100 m, 台缘带厚度明显大于台内, 孔隙(藻格架孔、粒间孔)、孔洞和裂缝构成主要的储集空间, 裂缝-孔隙-孔洞型储集层, 并有沥青充填[19]。选取充填孔洞的各期胶结物(葡萄花边状白云石胶结物)及围岩样品、充填藻纹层格架孔和藻砂屑粒间孔的各期胶结物及围岩样品、充填裂缝的各期胶结物样品, 开展定年研究, 为成岩-孔隙演化史的建立和油气运移前有效孔隙的确定提供重要信息。定年检测样品来源、层位、产状和检测目的信息见表1和图1。

表1 定年检测样品来源、层位、产状和检测目的信息表

图1 定年检测样品产状和特征
(a)(b)充填溶蚀孔洞的各期胶结物, 由①至⑤分别代表同心环边状白云石胶结物、放射状白云石胶结物、纹层状浅色和暗色白云石胶结物、晚期中— 粗晶白云石胶结物、最晚一期粗晶白云石胶结物, 先锋剖面灯二段; (c)充填溶蚀孔洞的各期胶结物, 由①至⑤分别代表围岩、叶片状白云石胶结物、粗晶白云石胶结物、石英、沥青, 磨溪22井灯二段5 418.70 m; (d)充填藻丘格架孔的等轴粒状白云石胶结物, 高石6井灯二段5 363.04 m; (e)泥晶白云岩, 由①至②分别代表围岩和裂缝中充填的白云石胶结物, 高石1井灯四段4 985.00 m; (f)裂缝及充填的白云石胶结物, 磨溪9井灯二段5 422.10 m

2 分析方法
2.1 碳酸盐矿物铀铅同位素定年方法和标样

使用传统的同位素稀释法对碳酸盐矿物进行U-Pb溶液法定年, 要求待测样品具有足够高的U、Pb含量, 能够从一块手标本上获得足够量的一组小样, 并且这组小样的U/Pb值有足够的变化范围, 能够拟合出一条207Pb/206Pb与238U/206Pb等时线, 使用这组数据拟合出的等时线和Tera-Wassenburg谐和线的下交点计算下交点年龄, 代表碳酸盐矿物的结晶年龄。然而, 碳酸盐矿物通常成因复杂, 具有多期多阶段性, 且易遭受后期改造, 故选择同源、同时、封闭体系且U/Pb值有一定变化范围能拟合出等时线的理想定年样品非常困难。而传统的同位素稀释法不适用于低铀碳酸盐矿物, 一是所需样品量大, 而大样品量全溶后U/Pb值均一化, 不同小样的U/Pb值变化范围小, 在Tera-Wassenburg谐和图上难以构筑理想的等时线以获取精确的下交点年龄; 二是取样和化学分离过程中非常容易造成污染, 对超净实验室本底要求非常高。因此, 过去20多年来, 这个方法虽然少数研究领域有几项成功的例子[13, 15], 但由于U含量、样品量、U/Pb值和实验室本底等诸多因素的限制, 难以应用于古老海相碳酸盐岩, 至今未见这方面任何报道。

近20年来, 随着激光剥蚀技术的日益兴盛, 激光原位U-Pb同位素定年技术已经广泛应用于测定高U矿物比如锆石、独居石、磷钇矿、榍石、金红石、磷灰石、石榴石等矿物的高精度年龄上, 成为地质年代学研究领域中最常用的测年方法。近年来, 一些低U矿物的激光原位U-Pb同位素定年工作也逐渐受到关注, 尤其碳酸盐矿物[24, 25, 26, 27, 28, 29, 30]。相比较于同位素稀释溶液法U-Pb测年, LA-(MC)-ICP-MS碳酸盐矿物微区U-Pb同位素定年技术具有原位、制样流程简单、样品消耗量小、低本底、空间分辨率高、分析速度快(单点分析仅需3~5 min)等优点。一般碳酸盐矿物具有明显低U特征(通常比锆石低2~4个数量级, 即其测量信号只有锆石的1/10 000~1/100), 故检测难度非常大。但碳酸盐本身铀含量有好几个数量级的变化, 分布范围从(1~2)× 10-6 mg/g到(1~10)× 10-3 mg/g, 一般在(0.05~0.50)× 10-3 mg/g。

前人使用扇形场(Sector-field)单接收ICP-MS(如Elements 2, Elements XR, Attom等), 或多接收MC-ICPMS(如Nu Plasma或Neptune)进行碳酸盐矿物U-Pb同位素测定, 初见成效[24, 25, 26, 27, 28, 29, 30]。但是, 扇形场单接收ICP-MS的灵敏度没有MC-ICPMS高, 而且因为采用跳峰法测量效率低。大多数MC-ICPMS虽然实现238U-208Pb-207Pb-206Pb的高效率静态测量, 但只有208Pb-207Pb-206Pb在离子计数器上测量, 238U往往在1011 Ω 的通用法拉第杯上测量, 通常当U含量小于0.2× 10-3 mg/g时, 238U是很难精确测量的。为了解决这个问题, 笔者在昆士兰大学的Nu Plasma II MC-ICPMS上最高质量端H10法拉第杯安装了1012 Ω 的高灵敏度前置放大器(使灵敏度高10倍)和一个离散打拿极倍增器(discrete dynode multiplier), 用于静态测量238U同位素, 前者使测试238U离子流的灵敏性比普通法拉第杯提高10倍, 后者提高100倍。当U含量足够高时(如大于0.1× 10-3 mg/g), 使用高灵敏度法拉第杯测量; 当U含量低时(如小于0.1× 10-3 mg/g), 使用倍增器测量。由于倍增器的背景极低, 一般U含量高于1× 10-6 mg/g时就能精确测量。除此以外, 笔者的Nu Plasma II MC-ICPMS在低质量范围还配有5个倍增器, 分别用于静态测量208Pb、207Pb、206Pb、204Pb(和204Hg干扰峰)和202Hg值。目前, 如果待测样品的U/Pb值足够高, 即使其U含量低至1× 10-6 mg/g, 笔者也能够对其有效定年。

笔者使用的激光剥蚀系统ASI RESOlution SE由澳大利亚科学仪器公司(Australian Scientific Instruments, 简称ASI)生产, 包括193 nm的ArF准分子激光器和Laurin Technic双室样品室。笔者首先使用RESOlution SE激光器和Thermo iCap-RQ四极杆ICPMS联动对样品靶的微量元素含量进行预扫描, 在微区尺度下能够快速清晰辨别样品中U、Pb、U/Pb值的变化范围。然后选择U/Pb值高普通Pb含量低的样品, 并根据U、Pb含量, 选择激光束斑直径, 在Nu Plasma II MC-ICPMS上进行年代测定。

使用LA-MC-ICP-MS法进行碳酸盐矿物U-Pb同位素年龄测定的一个重要难题是寻找合适的天然矿物标样, 并且为了避免不同矿物间的基体效应, 使用的矿物标样在成分和结构上应尽可能与待测样品相近。目前被应用于碳酸盐矿物测年的两个标样分别是ASH15E(采自以色列Negev沙漠的石笋)[28, 31, 32]和WC-1(美国新墨西哥城Whites City以西0.5 km的Walnut峡谷的方解石脉)[24, 25, 27, 30], 用它们来校正古老海相碳酸盐岩样品, 均存在一定局限性。ASH15E标样的同位素稀释溶液法标定年龄仅为3.001 Ma, 与古老海相碳酸盐岩年龄相比明显偏年轻得多, 不是古老海相碳酸盐岩的理想标样。WC-1标样的推荐年龄为254.4 Ma, 虽然与古老海相碳酸盐岩年龄相近, 但该标样的不均一性导致其自身可能存在3%~5%的不确定性[26, 28]。另外, 这两个国际标样极其有限, 不适合作为消耗量极大的笔者实验室的工作标样。

但是, 笔者从塔里木盆地阿克苏地区下寒武统肖尔布拉克组中发现了更加适合古老海相碳酸盐岩定年的标样AHX-1, 经标定后可作为实验室内部标准。AHX-1标样为充填沿断裂分布之孔洞的纯净方解石晶体, 样品纯净、均一、广泛分布, 是同源同期的成岩产物。过去一年以来, 笔者用LA-(MC)ICP-MS对该样品的微量元素和U-Pb同位素比值进行了反复检测, 获取数千组数据(见表2), 发现该样品非常适合激光原位U-Pb同位素测年, 绝大多数部位208Pb的含量几乎为零, 这说明样品没有普通Pb的影响, 206Pb和207Pb均为U衰变的产物, 238U的平均含量大约0.15× 10-3 mg/g, 虽然锆石标样91500的U含量比其大约70多倍, 但与笔者的MC-ICPMS仪器的U含量检测极限1× 10-6 mg/g相比, 完全满足检测需求。

表2 AHX-1标样微量和稀土元素分析

笔者在不同时段, 应用激光法在Nu Plasma II上将标样ASH15E与AHX-1反复对测20多次, 使用ASH15E的推荐年龄3.001 Ma[24]来标定AHX-1的年龄, 所获得的20多组年龄的加权平均值为(209.8± 1.3) Ma。图2a为以ASH15E为标样, 检测AHX-1的其中一次年龄为(209.1± 1.2)Ma。同时, 笔者假定AHX-1为标样, 检测ASH15E和WC-1的年龄, 测得的年龄分别为(2.960± 0.025)Ma(见图2b)和(259.0± 4.1)Ma(见图2c)。笔者的多次测量结果表明, AHX-1是一种更加理想的碳酸盐矿物标样, 每次测量不但具有很高的稳定性, 绝大多数数据点分布在谐和线的下交点附近, 有很低的普通Pb, 少数数据的分布在等时线上, 其结果是给出更小的年龄误差和上交点误差(上交点代表样品的普通铅组成), 更适合做标样。但是, 目前该标样目前只作为实验室工作标样, 不同测年方法(溶液法和激光法)间的标定工作正在开展中。笔者以AHX-1为假定的标样, 测量四川盆地五花洞北剖面泥盆系观雾山组方解石胶结物, U平均含量0.01× 10-3 mg/g, 测得的年龄为(244.3± 2.1)Ma(见图2d), 与以ASH15E和WC-1为标样测得的年龄相当(分别为(246.1± 2.3)Ma和(245.4± 1.7)Ma), 进一步验证了AHX-1标样的有效性。

图2 以AHX-1为工作标样、ASH15E和WC-1作为未知样品的一次LA-MC-ICPMS U-Pb定年结果

2.2 未知样品制备、测试和数据处理

符合激光原位U-Pb同位素定年的样品需同源、同期和均一的方解石或白云石晶体, 可以是围岩, 也可以充填于孔隙、孔洞或裂缝中。首先挑选符合定年要求的样品, 然后切割、清洁后进行制靶。本研究激光靶的制备在杭州地质研究院碳酸盐岩储层重点实验室完成。激光剥蚀(LA)碳酸盐靶的制备方法与SHRIMP锆石靶的制备相似[33]。为了消除制样中的Pb污染, 样品测试前还需要在超净室对待测样靶做进一步超净处理。

检测前先用激光剥蚀电感耦合等离子体质谱仪(LA-ICP-MS)对样品进行微量和稀土元素原位分析, 尤其是238U、206Pb、207Pb、208Pb的含量, 帮助笔者判断样品是否适合激光原位U-Pb同位素定年。激光系统工作条件是输出能量为3 J/cm2, 激光束斑直径视结构组分大小和U的含量可做不同选择, 一般使用100 µ m束斑, 剥蚀频率10 Hz。装入样品靶后, 气体连续冲洗样品池约2 h, 除去样品池和气路中可能存在的普通Pb。采样方式为单点剥蚀, 单点分析时间一般为3 min, 单点剥蚀之前对样品点预剥蚀2 s以消除表面Pb污染。信号较高的238U使用法拉第杯(Faraday cup)接收, 信号较低的207Pb、206Pb、208Pb、232Th、204Pb(+204Hg)、202Hg用离散打拿极倍增器接收, 超低238U含量的信号用IC5接收。在联接激光之前必须用Pb、Th、U的混合测试液对质谱仪Nu Plasma II进行质量标定和杯结构、透镜参数进行优化。连接激光之后采用对NIST612线扫描的方式调节仪器参数, 使Th/U值接近于1, UO/U值小于0.3%, 并且最优化仪器灵敏度。

数据处理可以在线或离线完成, 先使用Iolite 3.6软件[34]处理原始数据, 以获得相应的同位素比值, 然后在Isoplot 3.0软件[35]上完成谐和图绘制及年龄计算。

3 应用实例

以四川盆地灯影组为例, 开展围岩和不同期次白云石胶结物的激光原位U-Pb同位素定年分析, 同时对平行样开展二元同位素(D47)、碳氧稳定同位素组成、微量和稀土元素、锶同位素组成、阴极发光、包裹体均一温度分析。定年分析结果见表3、图3和图4, D47或包裹体均一温度、碳氧稳定同位素组成、阴极发光分析结果见表3, 微量和稀土元素、锶同位素组成分析结果见表4。微量和稀土元素、锶同位素组成数据来自昆士兰大学地球科学学院放射性同位素实验室, 二元同位素(D47)数据来自加利福尼亚大学洛杉矶分校地球与空间科学系同位素实验室, 阴极发光、碳氧稳定同位素、包裹体均一温度数据来自中国石油集团碳酸盐岩储层重点实验室。激光原位U-Pb同位素定年数据来自中国石油集团碳酸盐岩储层重点实验室和昆士兰大学地球科学学院放射性同位素实验室, 两家共同开发激光原位U-Pb同位素定年技术, 标样测定在澳方MC-ICP-MS(Nu Plasma Ⅱ )上完成, 灯影组样品年龄测定在中方LA-ICP-MS(Element XR)上完成。

表3 围岩和白云石胶结物U-Pb同位素年龄和地球化学特征[19]

图3 围岩和充填孔洞的白云石胶结物激光原位U-Pb同位素年龄
(a)围岩年龄(592± 24)Ma, 样品号GC-Z2dn2-B3; (b)同心环边状白云石胶结物年龄(546± 7.6)Ma, 样品号XF-Z2dn2-S5; (c)放射状白云石胶结物年龄(516± 10)Ma, 样品号XF-Z2dn2-S4; (d)纹层状浅灰— 暗色白云石胶结物年龄(482± 14)Ma, 样品号XF-Z2dn2-S4; (e)充填孔洞的晚期中— 粗晶白云石胶结物年龄(248± 27)Ma, 样品号GC-Z2dn2-B5; (f)充填孔洞的最晚一期粗晶白云石胶结物年龄(115± 69)Ma, 样品号GC-Z2dn2-B3

图4 充填孔隙和裂缝的白云石胶结物激光原位U-Pb同位素年龄
(a)充填藻丘格架孔的等轴粒状白云石胶结物年龄(545.7± 8.5)Ma, 样品号GS6-Z2d2-1; (b)充填孔隙的叶片状白云石胶结物年龄(499± 25)Ma, 样品号MX22-Z2d2-2; (c)充填孔隙的中— 粗晶白云石胶结物年龄(468± 12)Ma, 样品号MX22-Z2d2-3; (d)泥晶白云岩中充填裂缝的白云石胶结物年龄(472± 21)Ma, 样品号MX22-Z2d2-4; (e)藻纹层白云岩中充填裂缝的白云石胶结物年龄(41± 10)Ma, 样品号GS1-Z2d4-2; (f)藻纹层白云岩中裂缝的粗晶白云石胶结物年龄(41± 35)Ma, 样品号MX9-Z2d2-1

表4 围岩及白云石胶结物微量和稀土元素、锶同位素组成地球化学特征

表4虽然只罗列了6个测年样品的微量和稀土元素、锶同位素的地球化学特征, 但基本代表了本次用于测年的20个样品的地球化学特征。从表4的地球化学特征可知, 238U的含量较高, 6类结构组分的均值分别达到3.50× 10-3, 3.67× 10-3, 1.97× 10-3, 1.57× 10-3, 5.88× 10-3, 2.69× 10-3 mg/g, 远高于检测的极限值, 来自于238U衰变的206Pb和207Pb的含量也较高, 是理想的测年样品, 只是普通Pb含量偏高, 不如AHX-1标样理想, 在计算年龄时需考虑普通Pb的干扰。

4 讨论

四川盆地灯影组白云岩储集空间主要发育于藻纹层、藻叠层和藻砂屑白云岩中, 具有相控性, 类型主要有孔洞(2~100 mm)、孔隙(0.01~2.00 mm)和裂缝(见图3)。孔洞和孔隙均有原生、表生溶蚀和埋藏-热液溶蚀3种成因观点[19, 36, 37], 裂缝的期次及对孔洞、孔隙发育与充填的影响更是认识的空白。本文根据成岩产物的激光原位U-Pb同位素定年数据及二元同位素、碳氧同位素、锶同位素、微量和稀土元素地球化学特征, 结合构造-埋藏史、盆地热史和烃源岩生烃史, 分析四川盆地灯影组白云岩储集层的孔隙成因和成岩-孔隙演化史, 为油气运移前有效孔隙评价提供依据。

4.1 围岩和孔洞缝充填物绝对年龄

白云岩围岩:测得2个围岩(藻纹层或藻叠层白云岩)年龄的数据分别为(584± 32)Ma和(592± 24)Ma, 与Ediacaran系年龄(542~635 Ma)相当, 代表了地层年龄, 还可能反映了早期白云石化的年龄, 白云石化被认为与同沉积期蒸发的气候背景有关, 保留藻纹层和藻叠层的原岩结构[38]。阴极发光下不发光、碳同位素低正值(1‰ ~3‰ )、氧同位素低负值(-4‰ ~-1‰ )、锶同位素均值0.708 718(与同期海水相当), 是早期干旱氧化的萨布哈背景典型的地球化学特征[2]

孔洞充填物:孔洞中发育5期白云石胶结物, 构成葡萄花边状构造, 代表了完整的胶结充填序列, 大孔洞的胶结充填序列完整, 甚至有残留孔洞, 小孔洞的胶结充填序列不完整, 无残留孔洞。大多数学者认为葡萄花边状构造是沉积成因的, 孔洞为沉积原生孔、表生溶扩孔或埋藏溶蚀孔[39, 40, 41]

同心环边状白云石胶结物测得3个年龄数据, 分别为(546± 7.6)Ma、(545± 6)Ma和(545± 12)Ma。放射状白云石胶结物测得2个年龄数据, 分别为(516± 10)Ma和(514± 14)Ma。纹层状浅灰— 暗色白云石胶结物测得2个年龄数据, 分别为(482± 14)Ma和(487± 21)Ma。这3期胶结物均形成于早加里东期, 推测与工区的兴凯运动、郁南运动和都匀运动有关。中— 粗晶白云石胶结物测得的年龄为(248± 27)Ma, 可能与晚海西期东吴运动导致的热液活动有关。最晚一期粗晶白云石胶结物测得2个年龄数据, 分别为(20± 130)Ma和(115± 69)Ma, 误差虽然很大(误差大的原因是铀含量低但普通铅很高所致), 但均反映很年轻的胶结物, 可能与喜马拉雅期乐山— 资阳强烈褶皱导致的热液活动和充填有关。工区尚未发现晚加里东— 早海西期与广西运动相关的成岩产物。这些认识说明孔洞形成时间早于同心环边状白云石胶结物年龄, 应该形成于埋藏前, 与多幕次的桐湾运动导致的地层暴露和大气淡水淋溶有关, 非埋藏或热液溶蚀的产物。五期白云石胶结物由早至晚, 阴极发光强度逐渐增强、D47温度或包裹体均一温度逐渐升高、氧同位素逐渐偏负(见表3)的趋势佐证了孔洞充填过程中深埋逐渐加大、温度逐渐升高的特征[2], 并遭受晚海西期和喜马拉雅期热液活动的改造。

需要指出的是, 这5期胶结物的锶同素变化反映了中— 粗晶和粗晶白云石胶结物的热液活动痕迹。同心环边状白云石胶结物、放射状白云石胶结物、纹层状浅灰— 暗色白云石胶结物的锶同位素均值分别为0.708 732、0.709 351、0.709 186, 基本发映同期海水的锶同位素特征, 而中— 粗晶白云石胶结物、粗晶白云石胶结物的锶同位素均值分别为0.711 473和0.712 029, 与同期海水相比明显偏高, 是受到来自深部热液影响。深部热液在上升过程中会受到通道中围岩的放射性成因Sr的影响而偏高。

5期孔洞充填物的绝对年龄数据为葡萄花边状白云岩和孔洞的成因新解提供了证据, 葡萄花边状构造显然不是沉积成因的, 是埋藏后多期胶结作用的产物, 孔洞形成于沉积期(原生孔)或经历早表生的大气淡水溶蚀作用(溶扩孔)。

孔隙充填物:孔隙中主要发育3期白云石胶结物, 部分或完全充填孔隙。等轴粒状白云石胶结物测得的年龄为(545.7± 8.5)Ma, 与充填孔洞的同心环边状白云石胶结物年龄相当。叶片状白云石胶结物测得的年龄为(499± 25)Ma, 与充填孔洞的放射状白云石胶结物年龄相当。中— 粗晶白云石胶结物测得2个年龄数据, 分别为(457± 17)Ma和(468± 12)Ma, 略早于充填孔洞的纹层状浅灰— 暗色白云石胶结物年龄。这3期胶结物均形成于早加里东期, 但反映了连续充填的过程, 未见晚加里东、海西和喜马拉雅期热液矿物(鞍状白云石、闪锌矿等)对孔隙的充填, 与同期孔洞充填物具相似的阴极发光、D47温度或包裹体均一温度、氧同位素特征。晚加里东、海西和喜马拉雅期胶结充填物的缺失是因为孔隙不如孔洞大的缘故, 早加里东期的3期胶结物就足以把孔隙充填满, 并非同期没有埋藏-热液活动。

裂缝充填物:裂缝中白云石胶结物测得4个年龄数据, 分别为(472± 21)Ma、(20± 130)Ma、(41± 10) Ma和(41± 35)Ma, 第1个年龄数据代表早加里东期断裂活动和白云石胶结物的年龄, 充填孔隙的中— 粗晶白云石胶结物和充填孔洞的纹层状浅灰— 暗色白云石胶结物可能与这期断裂活动有关。后3个年龄数据差别虽然很大, 但均反映很年轻胶结物的年龄, 代表喜马拉雅期断裂活动和白云石胶结物的年龄, 充填孔洞的最晚一期粗晶白云石胶结物可能与这期断裂活动有关。虽然未测得与晚海西期年龄相当的充填裂缝的白云石胶结物年龄, 但并不意味着工区没有晚海西期的断裂活动, 也可能没取到对应的样品, 充填孔洞的中— 粗晶白云石胶结物应该与这期断裂活动有关。阴极发光下呈橙黄色中等发光、高D47温度及氧同位素高负值均佐证了充填裂缝白云石胶结物的热液, 受断裂活动控制。

孔隙和孔洞中早加里东3期白云石胶结物在产状和特征上有很大差异, 同样, 裂缝中充填的白云石胶结物与孔隙、孔洞中充填的白云石胶结物在产状和特征上也有很大差异, 常规技术手段很难将它们的成因期次相关联, 碳酸盐矿物的激光原位U-Pb同位素定年技术为不同产状和特征胶结物成因期次对比研究提供了技术手段。

孔洞缝中缺晚加里东— 早海西期白云石胶结物的绝对年龄, 这可能与志留纪— 石炭纪的广西运动、云南运动以整体抬升和剥蚀有关, 以平行不整合接触为特征, 构造挤压不强烈。孔洞缝中同样缺印支— 燕山期白云石胶结物的绝对年龄, 这可能与中晚三叠世之交的印支运动完成了四川盆地由海相向陆相沉积的转换, 随后进入持续和稳定地沉降, 并被陆相沉积充填有关, 构造总体稳定, 构造挤压不强烈。晚海西期是四川盆地乃至全球火山岩的发育期, 喜马拉雅期是四川盆地构造调整期, 构造活动强烈, 白云石胶结物的出现具有必然性。

4.2 成岩-孔隙演化史重建及应用

白云石胶结物绝对年龄的获得为四川盆地灯影组白云岩储集层的成因提供了新证据, 更为储集层成岩-孔隙演化史和有效孔隙评价提供了手段。灯影组白云岩储集层的储集空间形成于埋藏前的沉积和表生环境, 既有沉积原生孔隙(藻格架孔、粒间孔), 又有溶蚀扩大的孔洞, 孔隙和孔洞中充填的第1期白云石胶结物年龄与地层年龄相当, 足以说明这些孔隙和孔洞不是埋藏溶蚀作用的产物, 灯影组白云岩储集层的埋藏过程实际上是孔隙和孔洞逐渐被充填的过程。根据白云石胶结物的绝对年龄, 孔洞的充填作用发生在早加里东、晚海西、喜马拉雅期3个阶段, 孔隙的充填作用主要发生于早加里东期, 裂缝作为成岩介质的运移通道为充填孔洞和孔隙的胶结物提供了物源, 未被胶结物充填的残留孔洞、孔隙和裂缝构成了主要储集空间。

参照Clyde H Moore[2]的初始孔隙度值和镜下残留孔隙、胶结物分布面积累加, 初始孔隙度选30%, 表生溶蚀作用增孔10%, 埋藏前的总孔隙度达40%。早加里东期3期白云石胶结物使平均孔隙度由40%下降到15%, 晚海西期白云石胶结物使平均孔隙度由15%下降到10%, 并一直保持到燕山末期, 喜马拉雅期白云石胶结物使平均孔隙度由10%下降到8%, 裂缝对孔隙的贡献不大, 主要是作为埋藏-热液的通道, 且大多被充填。据此建立了灯影组白云岩储集层的成岩-孔隙演化史(见图5), 综合灯影组的构造-埋藏史[42]、盆地热史[43]和寒武系筇竹寺组烃源岩的生烃史[44], 就可对油气运移时间、油气运移前孔隙和成藏期次作出评价。

图5 四川盆地灯影组构造-埋藏史、成烃史、成岩-孔隙演化史和油气成藏史图

四川盆地高石梯— 磨溪构造灯影组气藏的烃源被认为来自寒武系筇竹寺组[44], 随着加里东期的持续埋藏, 志留纪末筇竹寺组烃源岩开始生烃并发生初次运移和成藏, 此时的有效孔隙度可达15%, 以残留孔洞为主。志留纪末油藏随着泥盆纪— 石炭纪的构造抬升发生裂解, 原油裂解气逸散或聚集成藏, 残留薄膜状沥青主要分布于孔洞中, 此时的孔隙度可以达到12%~15%。石炭纪末的持续深埋, 筇竹寺组烃源岩进入生烃高峰期, 在海西末期— 印支期(晚二叠世末— 三叠纪)聚集成藏, 是主成藏期, 此时的孔隙度可达12%。随着埋深的持续加大, 原油发生裂解, 形成斑块状沥青充填于孔洞中, 原油裂解气逸散或聚集形成燕山期气藏, 安岳气田就属于燕山期定型的气藏, 此时的孔隙度可达8%~10%。喜马拉雅期为气藏的调整期, 燕山期气藏经喜马拉雅构造运动改造发生调整, 部分被调整到喜马拉雅期构造圈闭中聚集成藏, 威远气田就属于喜马拉雅期定型的气藏, 此时的孔隙度可达8%。

综上所述, 碳酸盐矿物的激光原位U-Pb同位素定年技术的开发不仅解决了碳酸盐矿物的绝对年龄问题, 而且还可以通过碳酸盐胶结物的绝对年龄和含量, 恢复储集层的成岩-孔隙演化史, 结合构造-埋藏史和生烃史, 评价油气运移前的有效孔隙度和成藏有效性问题。前人[45, 46]主要通过油气包裹体均一温度确定成藏期次, 油气运移前有效孔隙度和成藏有效性评价为油气成藏期次的确定开辟了更为有效的途径。

5 结论

通过标定年龄为209.8Ma实验室工作标样AHX-1的开发、激光剥蚀进样系统与多接收电感耦合等离子体质谱仪的连用, 解决了溶液法难以实现的古老海相碳酸盐岩标样、取样和超低U、Pb含量(极限值大于1× 10-6 mg/g)样品的测年难题, 建立了适用于古老海相碳酸盐岩的激光原位U-Pb同位素定年技术。

通过充填孔洞、孔隙和裂缝中不同期次白云石胶结物的测年, 指出灯影组白云岩储集层的埋藏成岩过程主要是原生孔隙(藻格架孔、粒间孔)和表生溶蚀孔洞逐渐被充填的过程。孔洞的充填作用发生在早加里东、晚海西、喜马拉雅期3个阶段, 孔隙的充填作用主要发生于早加里东期, 3期裂缝白云石充填物与孔洞、孔隙中3期白云石胶结物年龄有很高的吻合度, 明确了充填孔洞、孔隙和裂缝各期白云石胶结物的成因对应关系。未被胶结物充填的残留孔洞、孔隙和裂缝构成了主要储集空间, 据此, 建立了川中地区灯影组白云岩储集层的成岩-孔隙演化史。

四川盆地灯影组白云岩储集层的成岩-孔隙演化史认识与该地区的构造-埋藏史、盆地热史具有很高的吻合度, 说明测年数据的可靠性和激光原位U-Pb同位素定年技术的有效性, 结合寒武系筇竹寺组烃源岩生烃史, 为古老海相碳酸盐岩储集层胶结物形成时间确定、成岩-孔隙演化研究和油气运移前有效孔隙评价提供了新的方法。

参考文献
[1] KERANS C. Karst-controlled reservoir heterogeneity in Ellenburger Group carbonates of west Texas[J]. AAPG Bulletin, 1988, 72(10): 1160-1183. [本文引用:1]
[2] MOORE C H. Carbonate reservoirs-porosity evolution and diagenesis in a sequence stratigraphic framework[M]. Amsterdam: Elsevier, 2001: 293-298. [本文引用:5]
[3] JAMES N P, CHOQUETTE P W. Paleokars[M]. Berlin: Springer- Verlag, 1988. [本文引用:1]
[4] LUCIA F J. Carbonate reservoir characterization[M]. Berlin: Springer-Verlag, 1999. [本文引用:1]
[5] CHOQUETTE P W, PRAY L C. Geologic nomenclature and classification of porosity in sedimentary carbonates[J]. AAPG Bulletin, 1970, 54(2): 207-250. [本文引用:1]
[6] LUCIA F J, MAJOR R P. Porosity evolution through hypersaline reflux dolomitization[J]. Sedimentary Geology, 1994, 21: 325-341. [本文引用:1]
[7] PURSER B H, AISSAOUI B A. Nature, origins and evolution of porosity in dolomites[J]. Sedimentary Geology, 1994, 21: 283-308. [本文引用:1]
[8] DAVIS G R, SMITH J L B. Structurally controlled hydrothermal dolomite reservoir facies: An overview[J]. AAPG Bulletin, 2007, 91(9): 1342-1344. [本文引用:1]
[9] 赵文智, 沈安江, 乔占峰, . 白云岩成因类型、识别特征及储集空间成因[J]. 石油勘探与开发, 2018, 45(6): 1-13.
ZHAO Wenzhi, SHEN Anjiang, QIAO Zhanfeng, et al. A research on genetic types and distinguished characteristics of dolostone, and the origin of dolostone reservoirs[J]. Petroleum Exploration and Development, 2018, 45(6): 1-13. [本文引用:1]
[10] 王兆荣, 彭子成, 陈文寄, . 腾冲地区年轻火山岩高精度热电离质谱(HP-TIMS)铀系法年龄研究[J]. 科学通报, 1999, 44(17): 1878-1881.
WANG Zhaorong, PENG Zicheng, CHEN Wenji, et al. Dating by high-precision Thermal Ionization Mass Spectrometry (TIMS) to the Paleotemperature stalagmites from cave meaning[J]. Chinese Science Bulletin, 1999, 44(17): 1878-1881. [本文引用:1]
[11] ZHAO Jianxin, NEIL D T, FENG Yuexing, et al. High-precision U-series dating of very young cyclone-transported coral reef blocks from Heron and Wistari reefs, southern Great Barrier Reef, Australia[J]. Quaternary International, 2009, 195(1/2): 122-127. [本文引用:1]
[12] WOODHEAD J, HELLSTROM J, MAAS R, et al. Taylor "U-Pb geochronology of speleothems by MC-ICP-MS[J]. Quaternary Geochronology, 2009, 1(3): 208-221. [本文引用:1]
[13] RASBURY E T, COLE J M. Directly dating geologic events: U-Pb dating of carbonates[J]. Reviews of Geophysics, 2009, 47(3): 4288-4309. [本文引用:2]
[14] PICKERING R, KRAMERS J D. A re-appraisal of the stratigraphy and new U-Pb dates at the Sterkfontein hominin site, South Africa[J]. Journal of Human Evolution, 2010, 56: 70-86. [本文引用:1]
[15] WOODHEAD J, PICKERING R. Beyond 500ka: Progress and prospects in the UPb chronology of speleothems, and their application to studies in palaeoclimate, human evolution, biodiversity and tectonics[J]. Chemical Geology, 2012, 322/323: 290-299. [本文引用:2]
[16] HILL C A, POLYAK V J, ASMEROM Y. Constraints on a Late Cretaceous uplift, denudation, and incision of the Grand Canyon region, southwestern Colorado Plateau, USA, from U-Pb dating of lacustrine limestone[J]. Tectonics, 2016, 35(4): 896-906. [本文引用:1]
[17] 邓胜徽, 樊茹, 李鑫, . 四川盆地及周缘地区震旦(埃迪卡拉)系划分与对比[J]. 地层学杂志, 2015, 39(3): 239-254.
DENG Shenghui, FAN Ru, LI Xin, et al. Subdivision and correlation of the Sinian (Ediacaran) system in the Sichuan basin and its adjacent area[J]. Journal of Stratigraphy, 2015, 39(3): 239-254. [本文引用:1]
[18] 李英强, 何登发, 文竹. 四川盆地及邻区晚震旦世古地理与构造-沉积环境演化[J]. 古地理学报, 2013, 15(2): 231-245.
LI Yingqiang, HE Dengfa, WEN Zhu, et al. Palaeogeography and tectonic-depositional environment evolution of the Late Sinian in Sichuan Basin and adjacent areas[J]. Journal of Palaeogeography, 2013, 15(2): 231-245. [本文引用:1]
[19] 陈娅娜, 沈安江, 潘立银, . 微生物白云岩储集层特征、成因和分布: 以四川盆地震旦系灯影组四段为例[J]. 石油勘探与开发, 2017, 44(5): 704-715.
CHEN Ya’na, SHEN Anjiang, PAN Liyin, et al. Characteristics, origin and distribution of microbial dolomite reservoirs: A case study of Sinian Dengying Formation in the Sichuan Basin[J]. Petroleum Exploration and Development, 2017, 44(5): 704-715. [本文引用:4]
[20] 汪泽成, 姜华, 王铜山, . 四川盆地桐湾期古地貌特征及成藏意义[J]. 石油勘探与开发, 2014, 41(3): 305-312.
WANG Zecheng, JIANG Hua, WANG Tongshan, et al. Paleo-geomorphology formed during Tongwan tectonization in Sichuan Basinand its significance for hydrocarbon accumulation[J]. Petroleum Exploration and Development, 2014, 41(3): 305-312. [本文引用:1]
[21] 刘树根, 王一刚, 孙玮, . 拉张槽对四川盆地海相油气分布的控制作用[J]. 成都理工大学学报(自然科学版), 2016, 43(1): 1-23.
LIU Shugen, WANG Yigang, SUN Wei, et al. Control of intracratonicsags on the hydrocarbon accumulations inthe marine strata acrossthe Sichuan Basin, China[J]. Journal of Chengdu University of Technology(Science & Technology Edition), 2016, 43(1): 1-23. [本文引用:1]
[22] 李伟, 易海永, 胡望水, . 四川盆地加里东古隆起构造演化与油气聚集的关系[J]. 天然气工业, 2014, 34(3): 8-15.
LI Wei, YI Haiyong, HU Wangshui, et al. Tectonic evolution of Caledonian paleohigh in the Sichuan Basin and its relationship with hydrocarbon accumulation[J]. Natural Gas Industy, 2014, 34(3): 8-15. [本文引用:1]
[23] 魏国齐, 杨威, 杜金虎, . 四川盆地震旦纪—早寒武世克拉通内裂陷地质特征[J]. 天然气工业, 2015, 35(1): 24-35.
WEI Guoqi, YANG Wei, DU Jinhu, et al. Geological characteristics of the Sinian-Early Cambrian intracratonic rift, Sichuan Basin[J]. Natural Gas Industry, 2015, 35(1): 24-35. [本文引用:1]
[24] LI Q, PARRISH R R, HORSTWOOD M S, et al. U-Pb dating of cements in Mesozoic ammonites[J]. Chemical Geology, 2014, 376: 76-83. [本文引用:4]
[25] COOGAN L A, PARRISH R R, ROBERTS N M. Early hydrothermal carbon uptake by the upper oceanic crust: Insight from in situ U-Pb dating[J]. Geology, 2016, 44(2): 147-150. [本文引用:3]
[26] ROBERTS N M, WALKER R J. U-Pb geochronology of calcite-mineralized faults: Absolute timing of rift-related fault events on the northeast Atlantic margin[J]. Geology, 2016, 44(7): 531-534. [本文引用:3]
[27] ROBERTS N M, RASBURY E T, PARRISH R R, et al. A calcite reference material for LA-ICP-MS U-Pb geochronology[J]. Geochemistry, Geophysics, Geosystems, 2017, 18(7): 2807-2814. [本文引用:3]
[28] NURIEL P, WEINBERGER R, KYLANDER-CLARK A R C, et al. The onset of the Dead Sea transform based on calcite age-strain analyses[J]. Geology, 2017, 45(7): 587-590. [本文引用:4]
[29] HANSMAN R J, ALBERT R, GERDEX A, et al. Absolutely ages of multiple generations of brittle structures by U-Pb dating of calcite[J]. Geology, 2018, 46(3): 207-210. [本文引用:2]
[30] GODEAU N, DESCHAMPS P, GUIHOU A, et al. U-Pb dating of calcite cement and diagenetic history in microporous carbonate reservoirs: Case of the Urgonian Limestone, France[J]. Geology, 2018, 46(3): 247-250. [本文引用:3]
[31] VAKS A, WOODHEAD J, BAR-MATTHEWS M, et al. Pliocene-Pleistocene climate of the northern margin of Saharan-Arabian Desert recorded in speleothems from the Negev Desert, Israel[J]. Earth and Planetary Science Letters, 2013, 368: 88-100. [本文引用:1]
[32] MASON A J, HENDERSON G M, VAKS A. An acetic acid-based extraction protocol for the recovery of U, Th and Pb from calcium carbonates for U-(Th)-Pb geochronology[J]. Geostand ards and Geoanalytical Research, 2013, 37(3): 261-275. [本文引用:1]
[33] 宋彪, 张玉海, 刘敦一. 微量原位分析仪器SHRIMP的产生与锆石同位素地质年代学[J]. 质谱学报, 2002, 23(1): 58-62.
SONG Biao, ZHANG Yuhai, LIU Dunyi. Introduction to the naissance of SHRIMP and its contribution to isotop geology[J]. Journal of Chinese Mass Spectrometry Society, 2002, 23(1): 58-62. [本文引用:1]
[34] PATON C, HELLSTROM J, PAUL B, et al. Iolite: Freeware for the visualisation and processing of mass spectrometric data[J]. Journal of Analytical Atomic Spectrometry, 2011, 26(12): 2508-2518. [本文引用:1]
[35] LUDWIG K R. User’s manual for ISOPLOT 3. 00: A geochronological toolkit for Microsoft excel[R]. Berkeley, California: Berkeley Geochronology Center, 2003. [本文引用:1]
[36] 施泽进, 梁平, 王勇, . 川东南地区灯影组葡萄石地球化学特征及成因分析[J]. 岩石学报, 2011, 27(8): 2263-2271.
SHI Zejin, LIANG Ping, WANG Yong, et al. Geochemical characteristics and genesis of grapestone in Sinian Dengying Formation in south-eastern Sichuan Basin[J]. Acta Petrologica Sinica, 2011, 27(8): 2263-2271. [本文引用:1]
[37] 冯明友, 强子同, 沈平, . 四川盆地高石梯—磨溪地区震旦系灯影组热液白云岩证据[J]. 石油学报, 2016, 37(5): 587-598.
FENG Mingyou, QIANG Zitong, SHEN Ping, et al. Evidence for hydrothermal dolomite of Sinian Dengying Formation in Gaoshiti-Moxiarea, Sichuan Basin[J]. Acta Petrolei Sinica, 2016, 37(5): 587-598. [本文引用:1]
[38] 沈安江, 赵文智, 胡安平, . 海相碳酸盐岩储集层发育主控因素[J]. 石油勘探与开发, 2015, 42(5): 545-554.
SHEN Anjiang, ZHAO Wenzhi, HU Anping, et al. Major factors controlling the development of marine carbonate reservoirs[J]. Petroleum Exploration and Development, 2015, 42(5): 545-554. [本文引用:1]
[39] 郝毅, 周进高, 陈旭, . 四川盆地灯影组“葡萄花边”状白云岩成因及地质意义[J]. 海相油气地质, 2015, 20(4): 57-64.
HAO Yi, ZHOU Jingao, CHEN Xu, et al. Genesis and geological significance of Upper Sinian Dengying dolostone with grape-lace shaped cement, Sichuan Basin[J]. Marine Petroluem Geology, 2015, 20(4): 57-64. [本文引用:1]
[40] 郝毅, 杨迅, 王宇峰, . 四川盆地震旦系灯影组表生岩溶作用研究[J]. 沉积与特提斯地质, 2017, 37(1): 48-54.
HAO Yi, YANG Xun, WANG Yufeng, et al. Supergene karstification in the Sinian Dengying Formation, Sichuan Basin[J]. Sedimentary Geology and Tethyan Geology, 2017, 37(1): 48-54. [本文引用:1]
[41] 邓韦克, 刘翔, 李翼杉. 川中震旦系灯影组储集层形成及演化研究[J]. 天然气勘探与开发, 2015, 38(3): 12-16.
DENG Weike, LIU Xiang, LI Yishan. Reservoir formation and evelution of Sinsian Dengying formation, Central Sichuan Basin[J]. Natural Gas Exploration & Development. 2015, 38(3): 12-16. [本文引用:1]
[42] 李伟, 刘静江, 邓胜徽, . 四川盆地及邻区震旦纪末—寒武纪早期构造运动性质与作用[J]. 石油学报, 2015, 36(5): 546-556.
LI Wei, LIU Jingjiang, DENG Shenghui, et al. The nature and role of Late Sinian-Early Cambrian tectonic movement in Sichuan Basin and its adjacent areas[J]. Acta Petrolei Sinica, 2015, 36(5): 546-556. [本文引用:1]
[43] 袁海锋, 刘勇, 徐昉昊, . 川中安平店—高石梯构造震旦系灯影组流体充注特征及油气成藏过程[J]. 岩石学报, 2014, 30(3): 727-736.
YUAN Haifeng, LIU Yong, XU Fanghao, et al. The fluid charge and hydrocarbon accumulation, Sinian reservoir, Anpingdian-Gaoshiti structure, central Sichuan Basin[J]. Acta Petrologica Sinica, 2014, 30(3): 727-736. [本文引用:1]
[44] 魏国齐, 王志宏, 李剑, . 四川盆地震旦系、寒武系烃源岩特征、资源潜力与勘探方向[J]. 天然气地球科学, 2017, 28(1): 1-13.
WEI Guoqi, WANG Zhihong, LI Jian, et al. Characteristics of source rocks, resource potential and exploration direction of Sinian and Cambrian in Sichuan Basin[J]. Natural Gas Geoscience, 2017, 28(1): 1-13. [本文引用:2]
[45] 李宏卫, 曹建劲, 李红中, . 油气包裹体在确定油气成藏年代及期次中的应用[J]. 中山大学研究生学刊(自然科学、医学版), 2008, 29(4): 29-35.
LI Hongwei, CAO Jianjin, LI Hongzhong, et al. Application of oil-gas inclusion to the reservoir formation time and accumulation periods of the oil-gas reservoirs[J]. Journal of the Graduates Sun Tai-sen University (Natural Science, Medicine), 2008, 29(4): 29-35. [本文引用:1]
[46] 刘文汇, 王杰, 陶成, . 中国海相层系油气成藏年代学[J]. 天然气地球科学, 2013, 24(2): 199-209.
LIU Wenhui, WANG Jie, TAO Cheng, et al. The geochronology of petroleum accumulation of China marine sequence[J]. Natural Gas Geosciences, 2013, 24(2): 199-209. [本文引用:1]