Thermal evolution modeling and present geothermal field of the Lishui Sag of East China Sea Basin
-
摘要: 东海陆架盆地丽水凹陷烃源岩埋深大、钻井少、油气发现少,至今对其生油气能力具有较大的争议。为明确烃源岩热演化史,从少量具备实测地层温度和镜质体反射率(Ro)资料的钻井出发,利用正演模拟方法,通过地温与地层深度相关性分析,将钻井地温研究结果外推到无井区。在恢复丽水凹陷探井岩石圈结构及其演化史的基础上,首次系统建立了丽水凹陷的区域地温场并分析其对烃源岩演化的作用。研究表明,古新统下段月桂峰组底界(T100)烃源岩大部分都已成熟,次洼中心局部过成熟;月桂峰组顶界(T90)烃源岩大部分处于低熟—成熟阶段,次洼中心局部高—过成熟;灵峰组顶界(T85)烃源岩处于未熟—低熟状态,西次洼中心有成熟烃源岩。烃源岩特征及演化分析表明,古新统下段月桂峰组烃源岩是丽水凹陷的主力生油气烃源岩。Abstract: Hydrocarbon source rocks are deeply buried in the Lishui Sag of East China Sea Basin, where drilling wells are rare and encountered only a few oil and gas discoveries. There are hot debates over the hydrocarbon generating capacity of the basin. Careful study has been made for this time to clarify the thermal evolutionary history of source rocks. Starting from drilling wells with measured formation temperature and vitrinite reflectance (Ro) data, using the forward simulation method, and through the correlation between temperature and stratum depth, we extrapolated the results of drilling temperature research to non-well areas. On the basis of restoring lithospheric structure and its evolution history, the regional geothermal field of the Lishui Sag has been systematically established for the first time. The impacts of geothermal field onto source rock evolution are analyzed. Data shows, most of the source rocks at the bottom of Yueguifeng Formation in the lower part of Paleocene are mature and the sub-depression center is partly over-matured; most of the source rocks at the top of Yueguifeng Formation are in the stage of low-maturity and maturity, the sub-depression center is partly high-mature to over-mature; the source rocks at the top of the Lingfeng Formation are immature to low-maturity, and there are mature source rocks at the center of the western sub-depression. According to the analysis of source rock characteristics and evolution, the source rocks in the lower Paleocene Yueguifeng Formation should be the main source rocks for hydrocarbon generation in the Lishui Sag.
-
位于西北干旱区的新疆,现代气候受西风系统和东亚冬季风系统的共同调控[1-2]。随着全球变暖(尤其是最近几十年来),在新疆以夏季降水为主的区域,降水量普遍增加,甚至出现强夏季暴雨[3],这对当地的生态环境和居民生活产生了重大影响。因此,了解新疆长时间尺度降水变化,对理解现代西风影响降水的特征和预测未来降水的趋势具有重要的参考意义。
目前对新疆及其周边过去千年气候变化的研究主要基于树轮[4-5]、冰芯[6]、湖泊沉积物[7-8]和风沙沉积[9]。然而,不同指标得出的结论对气候变化与控制机制间关系的解释存在差异。Chen[10]等通过评估过去千年区域有效湿度变化发现,有效湿度与温度在百年时间尺度上存在负相关关系,即湿(干)气候与冷(暖)特征。这种“冷湿/暖干”模式得到了塔里木盆地[11]和玛纳斯湖[12]的数据支持。但咸海沉积记录[13]表明,在小冰期进入中亚的水汽少于中世纪暖期,即呈现“冷干/暖湿”模式,并得到了来自西伯利亚南部基于阿米巴虫重建的地下水位深度记录的支持[14]。
可见,关于新疆及其周边湿度特征仍存在明显的争论,这可能归咎于区域湿度受到气温和降水的双重调控[15]和缺乏对区域降水的定量重建。因此,本文利用已报道的玉什库勒泥炭孢粉资料[16]对过去2 000年降水变化进行定量重建,期望该工作可以加深对阿尔泰山过去2 000年降水量特征及其驱动因素的理解。
1. 研究区概况
玉什库勒泥炭(46°45′~46°57′N、90°46′~90°61′E,
2374 m a.s.l)位于新疆北部青河县东北44 km,发育于冰川作用形成的U形谷底部,是小青格里河的发源地(图1)。该泥炭的地下水位深约0.5 m,以降水和积雪融水为主要补给[17]。现代植物种类主要有帕米尔苔草、阿尔泰苔草、灰藓等[18]。玉什库勒泥炭表层10 cm有机质含量高达90%以上,整体呈酸性[18]。研究区现代气候为温带大陆性气候。根据青河气象站记录,年均气温为−0.2 ℃,最冷和最热的月份分别是1月(−23.5 ℃)和7月(18.5 ℃);年均降水量为177.7 mm,暖季(4—9月)降水量占年均降水量的61%,冷季(10月至次年3月)降水量占39%[15]。2. 材料与方法
2.1 采样和测年
在玉什库勒泥炭(46°49′N、90°52′E,
2636 m a.s.l)获得长度为88 cm的岩芯,按照1 cm间隔分样,共获得88个样品。以10 cm间隔挑选植物残体作为测年材料,送至亚利桑那大学AMS 14C实验室进行年代的测试,相关结果见表1[15]。在R软件中调用Bacon程序建立深度-年龄模型[22](图2)。2.2 孢粉分析
实验室孢粉提取采用常规的HCl-NaOH-HF处理法[23],在OlympusBX-53光学显微镜下开展孢粉的鉴定和统计工作,每个样品统计至少600粒陆生孢粉[16],记录所出现的石松孢子数。在计算花粉百分含量时,以陆生花粉分类群的总和作为分母,蕨类孢子的百分比则根据陆生孢粉与蕨类孢子的总和计算。得出孢粉百分比后,用Tilia1.7.16软件绘制孢粉图,并对其进行有序聚类分析(CONISS),为花粉分区提供参考[16]。需要注意的是,因本剖面莎草科(Cyperceae)是湿地局地生长的植被,在定量重建时去除莎草科重新计算了孢粉百分比。
2.3 基于孢粉的定量重建
基于花粉的气候定量重建的可靠性取决于现代花粉校准集的质量[24-27]。为确保重建结果的可信度,需要考虑三点:(1)现代花粉数据集包含足够和均匀分布的采样点;(2)从总体中选择现代花粉数据的适当“范围策略”;(3)现代花粉数据集生物气候变化具有良好的代表性。为保证重建的可信度,本研究从东亚孢粉数据集和中国现代孢粉数据集中选择用于古气候重建的现代孢粉数据集[28-29],采用基于半径距离的筛选策略充分选择现代孢粉数据点[27]。
本文收集了研究区800 km范围内的现代表土花粉作为气候-花粉校准集,该校准集由627个样本和115个花粉物种组成,涵盖了现代生物气候空间范围,这些空间范围涵盖了玉什库勒过去2 000年生物气候时空范围。通过设置气候变量的上限和下限,确保堆叠花粉数据集中生物气候变量的良好代表性[19]。本文筛选现代花粉点的气候变量的上限和下限设定为:年均降水量(Pann):36~596 mm;年均温(Tann):−9.96~9.92 ℃;最冷月平均气温(Tmin):−15.96~4.28 ℃;最热月平均气温(Tmax):−3.96~17.1 ℃。通过ArcGIS10.2软件,从WorldClim数据集(第二版)(http://www.worldclim.org)获取现代花粉点的气候信息。
排序技术通过计算气候因子对花粉种群变化的解释量来评估用于重建的气候变量的适用性[19]。首先,对校准集中的花粉组合做去趋势对应分析(Detrended Correspondence Analysis, DCA)。当DCA第一梯度轴的长度>4时,基于正态分布的典型对应分析(Canonical Correspondence Analysis, CCA)更适合评估气候因子与花粉组合的关系;当DCA第一梯度轴的长度<3时,基于线性分布模型的冗余分析(Redundancy Analysis, RDA)则更为适合;当DCA第一梯度轴的长度为3~4时,两者皆可[19]。对于每个气候变量,用约束特征值与第一个无约束特征值的比值(λ1/λ2)来评估气候变量是否适合重建。一般来说,比值越高,气候变量越适合重建[30]。
匹配度检验法被用来评估重建模型的可信度,是通过计算每个化石孢粉组合与现代孢粉组合之间的平方弦距离来评估二者之间的匹配度[31]。当平方弦<5%时,说明两者之间有“良好的类比型”;当平方弦距离>5%且小于10%时,表明该化石孢粉样品在现代校准集中“有类比型”,反之则认为是“无类比型”[32]。用预测值与实测值的决定系数R2和预测的均方根误(Root Mean Squared Error of Prediction, RMSEP)[31]评价模型性能。其中,R2值越大,RMSEP越小,模型的性能越好。
3. 结果
3.1 孢粉记录的重述
Zhang等[16]对玉什库勒泥炭孢粉数据进行了详细描述。在所分析的88个样品中,共鉴定出55种不同陆生化石孢粉,包括15种木本孢粉(主要为松属(Pinus)、落叶松属(Larix)和桦木属(Betula))和40种草本孢粉(主要为莎草科(Cyperaceae)、禾本科(Poaceae)、蒿属(Artemisia)、藜科(Chenopodiaceae)和石竹科(Caryophyllaceae))。根据孢粉组合变化和CONISS分析结果,玉什库勒泥炭孢粉记录可以分为5个孢粉带(图3)。
孢粉带A(83~62 cm;12—520 AD)以草本植物为主,包括莎草科(64.81%~84.99%)、蒿属(6.67%~15.44%)、禾本科(1.85%~9.81%)和藜科(0.72%~6.16%),木本植物孢粉含量较少(1.06%)。
孢粉带B(61~43 cm;520—900 AD)相对于孢粉带A,莎草科孢粉含量降低(51.46%)、禾本科(25.16%)和石竹科(2.01%)孢粉含量增加。
孢粉带C(42~29 cm;900—1230 AD)的特征是松属(0.28%)、落叶松属(0.26%)和云杉属(Picea, 0.23%)孢粉含量增加,但仍以莎草科植物为主。
孢粉带D(28~15 cm;
1230 —1580 AD)与上一个孢粉带相比,莎草科(61.59%)孢粉含量降低,蒿属(20.32%)、禾本科(7.65%)、藜科(4.97%)、松属(0.57%)、落叶松属(0.45%)以及云杉属(0.37%)孢粉含量增加。孢粉带E(14~1 cm;
1580 —2020 AD)的特征是木本孢粉达到最大值(1.93%),蒿属(30.79%)、藜科(6.03%)、禾本科(11.04%)、石竹科(1.38%)孢粉含量持续增加,而莎草科孢粉含量(46.04%)持续减少。3.2 定量重建的可行性
根据排序分析,DCA的第一轴长度为3 SD(标准差单位)。选择CCA对现代表土花粉组合与现代气象因子进行分析。使用方差膨胀因子(Variance Inflation Factors, VIF)检验环境变量之间的共线性,结果如表2所示:校准集中的Tann与Tmax、Tmin显著相关;在分别排除Tmax和Tmin后,Tmax和Tmin气候参数的VIF值均大于30;在排除Tann后,校准集中的Pann、Tmax和Tmin的VIF值均小于30,这表明共线性的影响已经大大降低。因此,Tann被排除在重建之外。CCA轴1(6.55%)和轴2(1.59%)解释了8.14%的3个气候变量花粉组合之间的变异。花粉分类群与气候变量主要分布在轴1上,Pann与轴1之间的夹角较小(λ1/λ2值较大),表明Pann是重建的最佳变量(图4a)。
表 2 典型对应分析(CCA)与现代孢粉数据的统计和4个气候变量的统计汇总Table 2. Statistics of canonical correspondence analysis (CCA) and modern pollen data, and four climatic variables for the Yushenkule Peat sequence气候变量 VIF(含Tann) VIF(不含Tmax) VIF(不含Tmin) VIF(不含Tann) Pann 1.83 1.66 1.66 1.82 Tmin 71037.7 101.55 19.14 Tmax 88268.9 103.79 22.29 Tann 312085.7 107.37 98.34 图 4 典型对应分析(CCA)轴显示孢粉分类群和3个气候变量(a)与YSKL序列气候重建的统计显著性检验(b)Figure 4. Axes of canonical correspondence analysis (CCA) showing the relationships between pollen taxa and four climatic variables (a) for the Yushenkule Peat and statistical significance test of climatic reconstructions of Yushenkule Peat sequence using the MAT method (b)显著性检验(96.5%)结果表明,重建的Pann超过了预期的95%显著性水平(图4b)。该模型具有较高的R2(0.52)和相对较低的RMSEP(90.29)(图5)。另外,用平方弦距离法对校准集做相似性检验(图6b),结果显示98.8%的化石孢粉组合的平方弦距离<5%,具有较好的现代相似性,1.2%的化石孢粉组合的现代相似性较差,其平方弦距离为5%~10%,较差的相似性发生在约150 AD;其中蓝色和红色竖线分别表示5%分位数和10%分位数。综上所述,该模型适用于降水量的定量重建。
3.3 孢粉定量重建结果
基于玉什库勒泥炭定量重建结果如图6a所示,降水量在132~300 mm之间波动,其均值为205 mm。玉什库勒降水量的变化趋势特征划分为两个阶段:第I阶段(0—1010 AD)平均降水量高,上下波动幅度大,呈不显著下降;第II阶段(1010—2020 AD)降水量显著下降。这两个阶段可细分5个亚阶段:第1阶段(0—500 AD)的变化特征是平均高降水量和高振幅波动且基本保持不变;第2阶段(500—1010 AD)的降水量最高(平均值228 mm),且波动频率小于上一阶段(0—500 AD),在约750 AD降水量出现最大值(约300 mm)。第3阶段(1010—1300 AD)降水量下降,平均值为189 mm;第4阶段(1300—1600 AD)降水量高于前一个阶段,平均值为197 mm;第5阶段(1600—2020 AD)是过去两千年降水量最少的阶段(约165 mm),近200年降水量有所回升。
4. 讨论
4.1 区域对比
基于玉什库勒孢粉定量重建的降水量表明,0—1010 AD降水量较高,年均降水量为约225 mm;1010—2000 AD降水量降低,平均值为约183 mm,记录了降水量逐步下降的过程(图7a)。Zhang[16]等对玉什库勒泥炭孢粉主成分分析(PCA)(图7b)和Yang[15]等对泥炭沉积中δ13C(图7c)的研究均表现出同样的变化趋势。需要注意的是,1900 AD之后δ13C表现出变湿的趋势,可能是由于大量的苔草植物代替了泥炭藓植物所导致的[15]。基于孢粉定量重建的乌伦古湖[19](图7d)、喀纳斯湖[20](图7e)和娜仁夏泥炭地[21](图7f)降水也呈现变干的趋势。值得注意的是,因各序列年代分辨率不同而引起降水多寡的时间略有不同。过去2 000年区域降水量逐渐下降的过程也得到了模拟结果的支持[33]。
图 7 阿尔泰山南部降水序列对比a:玉什库勒年均降水量(本研究), b. 玉什库勒δ13C含量[15], c. 玉什库勒PCA轴2[16], d. 乌伦古湖降水量[19],e. 喀纳斯湖降水量[20],f. 娜仁夏降水量[21], g. 北大西洋涛动[35], h. 总太阳辐射[37], i. 阿尔泰山温度[4]。Figure 7. Comparison of the precipitation sequence in the southern Altai Mountainsa: Precipitation in the Yushenkule Peat (this study),b: δ13C-indcated moisture in the Yushenkule Peat[15], c:PCA axis 2 in the Yushengkule Peat[16],d: Precipitation in the Wulungu Lake[19], e: Precipitation in the Kanas Lake[20], f: Precipitation in the Narenxia Peat[21], g: NAO[36], h. Total solar irradiance[37], i: Temperature data in the Altai Mountains[4].结合阿尔泰山树轮所记录的气温序列[4]可以发现,罗马暖期(RWP, 0—500 AD)降水量和温度均处于较高状态,该时期的气候特征为“暖湿”。进入黑暗时代冷期(DACP, 500—800 AD),气候转为“冷湿”。中世纪暖期(MWP, 800—1400 AD)降水量降低,气候是“暖干”的。小冰期(LIA, 1400—1850 AD)降水量较低,温度降低,气候特征可以描述为“冷干”。现代暖期(CWP, 1850年至今)温度和降水量均持续上升,气候由“冷干”转向“暖湿”。可见,过去2 000年水热模式呈现了暖湿-冷湿-暖干-冷干-暖湿的演变过程,呈现既不符合“暖干-冷湿”模式,又不符合“暖湿-冷干”的模式。
4.2 驱动机制
基于玉什库勒泥炭孢粉定量重建的降水量,不仅与区域记录的降水量重建(图7)显示出一致性,还与来自欧洲的树轮和海洋沉积记录显示出同步变化[34-35]。这种一致性促使我们研究水汽通过西风环流从欧洲穿过欧亚大陆达到新疆的路径以及其与北大西洋涛动(NAO)的关系。但当前关于NAO和区域降水量的关系有两种观点,第一种观点认为NAO正相(即亚速尔高压和冰岛低压之间的气压差增加)有利于新疆降水增加[8],另一种认为负相NAO(即亚速尔高压和冰岛低压之间的气压差降低)有利于新疆降水增加[36]。
通过与来自特隆赫姆斯峡湾记录的初级生产力所记录的高质量NAO指数[35]对比可以看出,在0—1010 AD期间,区域降水特征与NAO指数在百年尺度上具有正向关系(图7g),这表明NAO正相驱动了区域降水量的增加。从机理角度来说,当NAO正相时,亚速尔群岛高压和冰岛低压之间的海平面压强增加,使更多的水汽通过西风环流输送到阿尔泰山。加上较高的总太阳辐射[37]驱动西风环流强度和北大西洋蒸发水汽的增加,为中纬度地区提供更多的水汽,导致阿尔泰山降水量增加。
1010 AD之后,降水量变化与NAO的关系变得不明晰,甚至在多世纪尺度上显示出负相关关系(图7g),这可能与总太阳辐射减弱及其关联的区域水汽循环变化有关(图7h)。详细来说,过去千年总太阳辐射有所减弱,引起中纬度西风带强度下降,且引起西风带南移[38],导致更少的水汽被输送至阿尔泰山。尽管该时段厄尔尼诺-南方涛动(ENSO)强度增加,增强的ENSO导致中亚低压槽的发生,导致来自阿拉伯海的水汽供应增加,使得阿尔泰山降水量增加[21, 39]。但ENSO输送的水汽通量有限,并未很好的弥补西风带南移所造成的区域降水量下降。可见,过去两千年阿尔泰山降水量受到NAO、总太阳辐射和ENSO等因素共同作用的影响。5. 结论
(1)阿尔泰山过去2000年降水量最大值为300 mm,最小值为132 mm,其均值为205 mm。过去2000年可分为湿润的0—1010 AD(降水量约224 mm)和偏干的1010—2000 AD(降水量约182 mm)。最大降水量出现在600—800 AD,年均量为288 mm。最小降水量出现在1500—2000 AD,年均量为162 mm。
(2)过去2000年水热模式呈现了暖湿-冷湿-暖干-冷干-暖湿的演变过程,该变化过程与NAO、太阳辐射变化引起的西风带的南移密切相关。
-
表 2 LS5井岩石圈模型参数
Table 2 Lithospheric modeling parameters of LS5
岩石圈原始厚度ho=125 km[21] 岩石圈现今厚度ho=80 km 上地壳放射性生热率Qr0=2.8 μW/m3 原始地壳厚度hc=35 km 现今地壳厚度hc=21.5 km 上地壳热导率=3.5 W/mK 原始上地幔厚度hm=90 km 现今上地幔厚度hm=52 km 下地壳热导率=3.1 W/mK 研究区现今纬度=27° 基底深度db=6 499 m 上地幔热导率=2.5 W/mK 软流圈顶部温度Tb=1 330 ℃ 上地壳密度ρc上 =2.7 g/cm3 上地壳比热容=900 J·kg−1·K−1 拉张时间ts=27.1 Ma 下地壳密度ρc下 =2.9 g/cm3 下地壳比热容=1 050 J·kg−1·K−1 地壳拉张系数βc=1.63 上地幔密度ρm =3.4 g/cm3 上地幔比热容=1 200 J·kg−1·K−1 表 3 丽水凹陷岩石密度和热属性参数
Table 3 Density and thermal properties of rocks encountered in Lishui Sag
岩性 密度
/(g/cm3)热导率
/(W/mK)生热率
/(μW/m3)比热容
/(J·kg−1·K−1)砾岩 2.7 4.18 0.85 820 砂岩 2.65 3.85 0.7 855 粉砂岩 2.65 3.35 0.96 910 泥岩 2.72 1.62 1.5 860 煤 1.6 0.46 0.38 1 300 表 4 丽水凹陷热流值及相关参数
Table 4 Heat flow values and related parameters in Lishui Sag
井号 测温深度/m 平均地温梯度
/(℃/100 m)基底热流
/(mW/m2)LS1井 2 665.8 2.78 48.26 LS2井 2 000.5 3.17 50.16 LS3井 2 798 2.96 49.06 LS4井 2 801 3.14 53.95 LS5井 2 791 3.21 56.08 LS6井 2 902.7 2.59 48.53 LS7井 3 022.4 3.43 54.86 LS8井 2 574~2 587.5 2.60 54.30 LS9井 3 117 2.69 41.78 LS10井 − 3.24 54.88 LS11井 2 032.67 3.26 58.53 LS12井 2 467.26 3.27 56.81 LS13井 − 3.26 53.39 LS14井 2 698.1 3.33 53.48 LS15井 1 050 3.65 62.79 LS16井 2 325 3.21 61.46 LS17井 4 015 2.72 55.56 LS18井 3 803 2.78 47.55 LS19井 3 040 2.95 44.76 -
[1] 葛和平, 陈志勇, 方来富, 等. 丽水凹陷油气成藏期次探讨[J]. 中国海上油气(地质), 2003, 17(1):44-50. [GE Heping, CHEN Zhiyong, FANG Laifu, et al. A discussion on hydrocarbon accumulation periods in Lishui Sag, East China Sea Basin [J]. China Offshore Oil and Gas (Geology), 2003, 17(1): 44-50. [2] 孙玉梅, 席小应. 东海盆地丽水凹陷油气源对比与成藏史[J]. 石油勘探与开发, 2003, 30(6):24-28. [SUN Yumei, XI Xiaoying. Petroleum reservoir filling history and oil-source correlation in the Lishui Sag, East China Sea Basin [J]. Petroleum Exploration and Development, 2003, 30(6): 24-28. doi: 10.3321/j.issn:1000-0747.2003.06.007 [3] 陈志勇, 吴培康, 吴志轩. 丽水凹陷石油地质特征及勘探前景[J]. 中国海上油气: 地质, 2000, 14(6):384-391. [CHEN Zhiyong, WU Peikang, WU Zhixuan. Petroleum geology and exploration potential of Lishui Sag [J]. China Offshore Oil and Gas: Geology, 2000, 14(6): 384-391. [4] 陈斯忠. 东海盆地主要地质特点及找气方向[J]. 中国海上油气(地质), 2003, 17(1):6-13, 19. [CHEN Sizhong. Main geological characteristics and gas exploration directions in East China Sea Basin [J]. China Offshore Oil and Gas (Geology), 2003, 17(1): 6-13, 19. [5] 郭永华, 于水, 葛玲. 东海盆地丽水凹陷LS36-1构造成藏机理研究[J]. 石油勘探与开发, 2003, 30(6):29-31. [GUO Yonghua, YU Shui, GE Ling. Formation of the LS36-1 oil and gas structure in the Lishui Sag, East China Sea Basin [J]. Petroleum Exploration and Development, 2003, 30(6): 29-31. doi: 10.3321/j.issn:1000-0747.2003.06.008 [6] 李宁. 丽水凹陷西次凹古近系烃源岩热演化特征研究[J]. 长江大学学报: 自科版, 2013, 10(10):46-48. [LI Ning. Study on the thermal evolution characteristics of hydrocarbon source rocks in Paleogene system of Lishui Sag [J]. Journal of Yangtze University: Natural Science Edition, 2013, 10(10): 46-48. [7] 冯晓杰, 蔡东升, 王春修, 等. 东海陆架盆地中新生代构造演化特征[J]. 中国海上油气(地质), 2003, 17(1):33-37. [FENG Xiaojie, CAI Dongsheng, WANG Chunxiu, et al. The meso-cenozoic tectonic evolution in East China Sea Shelf Basin [J]. China Offshore Oil and Gas (Geology), 2003, 17(1): 33-37. [8] 张田, 张建培, 张绍亮, 等. 平衡剖面技术在东海丽水凹陷构造演化研究中的应用[J]. 上海国土资源, 2014, 35(1):92-96. [ZHANG Tian, ZHANG Jianpei, ZHANG Shaoliang, et al. An application of the balanced cross-section technique: the tectonic evolution of Lishui Sag, the East China Sea Shelf Basin [J]. Shanghai Land & Resources, 2014, 35(1): 92-96. doi: 10.3969/j.issn.2095-1329.2014.01.021 [9] 田兵, 李小燕, 庞国印, 等. 叠合断陷盆地沉积体系分析——以东海丽水-椒江凹陷为例[J]. 沉积学报, 2012, 30(4):696-705. [TIAN Bing, LI Xiaoyan, PANG Guoyin, et al. Sedimentary systems of the superimposed rift-subsidence basin: Taking Lishui-Jiaojiang Sag of the East China Sea as an example [J]. Acta Sedimentologica Sinica, 2012, 30(4): 696-705. [10] 陈国俊, 李超, 梁建设, 等. 东海陆架盆地瓯江凹陷明月峰组沉积相及沉积特征分析[J]. 天然气地球科学, 2011, 22(5):760-769. [CHEN Guojun, LI Chao, LIANG Jianshe, et al. Sedimentary facies of Mingyuefeng formation in Oujiang Sag, East China Sea Basin [J]. Natural Gas Geoscience, 2011, 22(5): 760-769. [11] Kauerauf A I, Hantschel T. Fundamentals of Basin and Petroleum Systems Modeling[M]. Berlin, Heidelberg: Springer-Verlag, 2009: 1-476.
[12] Lin C S, Zhang Y M, Li S T, et al. Quantitative modelling of multiphase lithospheric stretching and deep thermal history of some tertiary rift basins in eastern China [J]. Acta Geologica Sinica, 2002, 76(3): 324-330.
[13] He Z Y, Crews S G, Corrigan J. Rifting and heat flow: Why the McKenzie model is only part of the story[C]//Basin Modeling Perspectives: Innovative Developments and Novel Applications. Hague, Netherlands: AAPG, 2007.
[14] 石广仁. 油气盆地数值模拟方法[M]. 3版. 北京: 石油工业出版社, 2004. SHI Guangren. Approaches of Petroleum Basin Modeling[M]. 3rd ed. Beijing: Petroleum Industry Press, 2004.
[15] 郝天珧, 徐亚, 胥颐, 等. 对黄海-东海研究区深部结构的一些新认识[J]. 地球物理学报, 2006, 49(2):458-468. [HAO Tianyao, XU Ya, XU Yi, et al. Some new understandings on deep structure in Yellow Sea and East China Sea [J]. Chinese Journal of Geophysics, 2006, 49(2): 458-468. doi: 10.3321/j.issn:0001-5733.2006.02.019 [16] 蔡学林, 朱介寿, 曹家敏, 等. 中国大陆及邻区岩石圈地壳三维结构与动力学型式[J]. 中国地质, 2007, 34(4):543-557. [CAI Xuelin, ZHU Jieshou, CAO Jiamin, et al. 3D structure and dynamic types of the lithospheric crust in continental China and its adjacent regions [J]. Geology in China, 2007, 34(4): 543-557. doi: 10.3969/j.issn.1000-3657.2007.04.001 [17] 周海廷, 姜效典, 李德勇, 等. 东海陆架盆地西湖凹陷岩石圈热流变性质[J]. 大地构造与成矿学, 2017, 41(3):481-490. [ZHOU Haiting, JIANG Xiaodian, LI Deyong, et al. Thermal-rheological property of lithosphere beneath Xihu Sag, East China Sea Shelf Basin [J]. Geotectonica et Metallogenia, 2017, 41(3): 481-490. [18] 高德章, 赵金海, 薄玉玲, 等. 东海及邻近地区岩石圈三维结构研究[J]. 地质科学, 2006, 41(1):10-26. [GAO Dezhang, ZHAO Jinhai, BO Yuling, et al. A study on lithosphere 3D structures in the East China Sea and adjacent regions [J]. Chinese Journal of Geology, 2006, 41(1): 10-26. doi: 10.3321/j.issn:0563-5020.2006.01.002 [19] 江为为, 刘少华, 郝天珧, 等. 应用重力资料估算东海冲绳海槽地壳厚度[J]. 地球物理学进展, 2002, 17(1):35-41. [JIANG Weiwei, LIU Shaohua, HAO Tianyao, et al. Using gravity data to compute crustal thickness of East China Sea and Okinawa Trough [J]. Progress in Geophysics, 2002, 17(1): 35-41. doi: 10.3969/j.issn.1004-2903.2002.01.004 [20] 周志远, 高金耀, 吴招才, 等. 东海莫霍面起伏与地壳减薄特征初步分析[J]. 海洋学研究, 2013, 31(1):16-25. [ZHOU Zhiyuan, GAO Jinyao, WU Zhaocai, et al. Preliminary analyses of the characteristics of Moho undulation and crustal thinning in East China Sea [J]. Journal of Marine Sciences, 2013, 31(1): 16-25. doi: 10.3969/j.issn.1001-909X.2013.01.002 [21] 仝志刚, 赵志刚, 杨树春, 等. 低勘探程度盆地烃源岩热演化及排烃史研究——以东海椒江凹陷为例[J]. 石油实验地质, 2012, 34(3):319-324, 329. [TONG Zhigang, ZHAO Zhigang, YANG Shuchun, et al. Research on thermal evolution and hydrocarbon expulsion history of source rocks in low-exploration basins: a case study on Jiaojiang Sag, East China Sea Basin [J]. Petroleum Geology & Experiment, 2012, 34(3): 319-324, 329. doi: 10.3969/j.issn.1001-6112.2012.03.016 [22] 方银霞, 刘建华. 东海的地壳结构特征[J]. 东海海洋, 2004, 22(3):9-17. [FANG Yinxia, LIU Jianhua. The crustal structure character of East China Sea [J]. Donghai Marine Science, 2004, 22(3): 9-17. [23] 栾锡武, 高德章, 喻普之, 等. 我国东海陆架地区新生代地层的热导率[J]. 海洋与湖沼, 2002, 33(2):151-159. [LUAN Xiwu, GAO Dezhang, YU Puzhi, et al. Thermal conductivity of the Cenozoic layer of East China Sea Shelf [J]. Oceanologia et Limnologia Sinica, 2002, 33(2): 151-159. doi: 10.3321/j.issn:0029-814X.2002.02.006 [24] 栾锡武, 赵一阳, 秦蕴珊, 等. 我国东海陆架区新生代地层岩石生热率研究[J]. 沉积学报, 2003, 21(4):634-639. [LUAN Xiwu, ZHAO Yiyang, QIN Yunshan, et al. Heat generation of the cenozoic layer of East China Sea Shelf [J]. Acta Sedimentologica Sinica, 2003, 21(4): 634-639. doi: 10.3969/j.issn.1000-0550.2003.04.015 [25] 仝志刚, 杨树春, 郝建荣. 盆地模拟技术应用与发展(内部资料). 2011. TONG Zhigang, YANG Shuchun, HAO Jianrong. Application and development of basin simulation technology(Internal data), 2011.
[26] Wygrala B P. Integrated study of an oil field in the Southern Po Basin, Northern Italy[D]. Doctor Dissertation of University of Cologne, Germany, 1989.
[27] 王斌, 赵永强, 罗宇, 等. 塔里木盆地草湖凹陷热演化与生烃史——基于IES软件盆地模拟技术[J]. 石油实验地质, 2010, 32(6):605-609. [WANG Bin, ZHAO Yongqiang, LUO Yu, et al. Thermal evolution and hydrocarbon generation in the Caohu Sag of the Tarim Basin—based on IES basin simulation technology [J]. Petroleum Geology & Experiment, 2010, 32(6): 605-609. doi: 10.3969/j.issn.1001-6112.2010.06.017 [28] Waples D W, Ramly M. A statistical method for correcting log-derived temperatures [J]. Petroleum Geoscience, 2001, 7(3): 231-240. doi: 10.1144/petgeo.7.3.231
[29] Waples D W, Pacheco J, Vera A. A method for correcting log-derived temperatures in deep wells, calibrated in the Gulf of Mexico [J]. Petroleum Geoscience, 2004, 10(3): 239-245. doi: 10.1144/1354-079302-542
[30] 仝志刚, 贺清, 何仕斌, 等. 东海西湖凹陷地温场及其对烃源岩的作用[J]. 石油实验地质, 2009, 31(5):466-471, 484. [TONG Zhigang, HE Qing, HE Shibin, et al. Geothermal field and its effect on source rock in the Xihu Sag, the East China Sea Basin [J]. Petroleum Geology & Experiment, 2009, 31(5): 466-471, 484. doi: 10.3969/j.issn.1001-6112.2009.05.006