A study on the hydrocarbon generation potential of the coal-bearing source rocks in the Pinghu Formation of Pingbei area, the Xihu depression
-
摘要: 西湖凹陷平北地区平湖组发育海陆过渡相的煤系烃源岩，以生烃潜力为主要依据对煤岩以及暗色泥岩的有机质丰度进行了评价，并建立了相应的评价标准，评价结果显示其有机质丰度均为中等―好，煤岩表现为S1+S2＞90 mg/g，HI＞100 mg/g，氯仿沥青“A”＞1.0%；暗色泥岩表现为TOC＞1%, S1+S2＞2 mg/g，氯仿沥青“A”＞0.05%，总烃“HC”＞150×10-6。干酪根元素分析显示煤的干酪根类型主要为Ⅱ1-Ⅱ2型，富氢的显微组分含量较高，壳质组中的树脂体、孢子体及角质体含量达到10%～20%，且含有约20%的富氢镜质体，表明煤具有较高的生油潜力；暗色泥岩的富氢显微组分低于煤岩，有机质类型主要为Ⅱ2-Ⅲ型，以生气为主。根据镜质体反射率Ro以及S1+S2/TOC等指标共同确定了平北地区烃源岩的生烃门限（3 100 m）和排烃门限（3 500 m）。热模拟实验显示平湖组煤具有油气共生的特点，具有较高的产油率，可达130 mg/g TOC，产气率可达125 mg/g TOC；暗色泥岩以生气为主，产油率最高仅为50 mg/g TOC，产气率可达185 mg/g TOC，且生气周期长。基于平北地区较重的原油碳同位素和生物标志物特征对比，认为煤成油在研究区广泛存在。Abstract: The source rocks of the Pinghu Formation in the Pingbei area of the Xihu depression consist of a set of coal-bearing deposits of marine-continental transitional facies. Organic matter abundance of the coal and dark mudstone were evaluated mainly based on potential hydrocarbon generation amount according to the evaluation criteria established. Results show that both the coal rocks and dark mudstones have high hydrocarbon generation potential and medium to high abundance of organic matter. The coal rocks are featured by S1+S2＞90 mg/g, HI＞100 mg/g,“A”＞1.0%, while the mudstones characterized by TOC＞1%, S1+S2＞2 mg/g, “A”＞0.05%, “HC”＞150×10−6. It is identified that Kerogen of the coal rocks belongs to typeⅡ1-Ⅱ2 according to the elemental data, and the content of hydrogen-rich macerals is very high. The total content of resinite, sporophyte and cuticle in exinite is up to 10%～20%, and the content of hydrogen-rich vitrinite is as high as 20%, indicating a high oil-generating potential. The content of hydrogen-rich macerals in mudstone is lower than that in coal rocks. The kerogen belongs to typeⅡ2 with higher gas-generating potential. Threshold depth of oil generation and expulsion in the Pingbei area were determined by Ro values and (S1+S2)/TOC values, and they are 3 100 m and 3 500 m respectively. The results of thermal simulation show that the coal rocks are prone to generating both oil and gas, with higher proportion of oil. Oil generation rate can reach 130 mg/g TOC and Gas generation rate about 125 mg/g TOC. The mudstone is gas-prone for a long history of gas generation. Data suggests that the oil generation rate is 50 mg/g TOC only, while the gas generation rate can reach 185 mg/g TOC. Coal-derived oil is widely existed in the Xihu depression according to the heavier crude oil carbon isotope and correlation of biomarkers between oils and source rocks.
-
大洋洋底广泛分布着海山等粗糙地形隆起，如以马尼亚纳俯冲带等为典型代表的西太平洋俯冲边缘分布着大量海山[1-2]。这些洋壳基底粗糙地形的俯冲对上覆板块具有重要的影响，包括导致上盘发生复杂的构造变形、乱序型逆冲活动、地貌侵蚀、影响板块间的震间耦合、流变特性、岩浆活动等[3-7]。前人通过构造物理模拟、数值模拟、大洋钻探和地震剖面解析等方法研究了海山或海底高原俯冲对增生楔内逆冲断层结构和活动性的影响[3,8-10]。海山等粗糙海底地形的俯冲可以导致反序型断裂作用和断层的再活化，促进增生楔底部的构造侵蚀和沉积物的俯冲，并可能导致富含流体的沉积物进入深部孕震带[11-12]。因此，研究海山俯冲可以更好地理解汇聚板块边缘弧前增生楔构造侵蚀过程，对理解上覆板块断裂体系的形成、演化以及对现今俯冲带地震活动性和海啸孕育的影响具有重要意义。
希库朗伊俯冲边缘（Hikurangi subduction margin）位于新西兰北岛外海以东，是南太平洋克马德克-汤加海沟俯冲带（Kermadec-Tonga subduction zone）的南延部分[13]（图1）。希库朗伊高原是目前全球规模最大的洋底大火成岩省之一，发育了数量众多的大型海山[14-16]。这些大型海山的俯冲如何影响希库朗伊俯冲带北缘的构造侵蚀过程和局部的板块应力体制是目前的前沿热点问题[17-23]。海山的位置和特征对于理解现今与海山俯冲相关的板间地震活动至关重要，并且对俯冲带的孕震行为、水文地质过程和长期的构造演化也有重要的控制作用[1,17,25-29]。目前希库朗伊边缘开展了大洋钻探和地球物理观测，包括国际大洋发现计划（IODP 372、375航次）、主动源地震探测（二维和NZ3D OBS实验）、大地电磁探测（MT），但是对于俯冲带增生楔内部断裂和深部海山的位置、结构及其演化还存在一定的争议[30-33]。希库朗伊俯冲带现今存在两种截然不同的海山构造模型：①低速海山模型和②高速海山模型（图2），两种模型中海山的位置存在显著差异，前者位于现今增生楔底部且海山严重破碎化，而后者位于更下倾处并表现为一个地震高速体[17,24,30,34]。
![]()
![]()
离散元方法（discrete element method, DEM）是一种基于颗粒性质的数值模拟方法，可以用来研究构造变形并定量观测和分析变形系统的应力和应变演化，能够帮助揭示断层活动性和力学性质的演化[35-36]。本文拟以希库朗伊俯冲带北缘典型构造剖面为例，利用2-D离散元数值模拟方法，研究海山俯冲与增生楔、板间巨型逆冲断层（megathrust）及巨型分支断层（mega-splay fault）之间的相互作用，探讨断裂体系的结构及其活动性的时空演化规律。本项研究将为更好地理解海山俯冲对增生楔的构造演化和断层活动性的影响机制，对探讨板间地震活动和慢滑移事件的形成提供重要启示。
1. 区域地质背景
希库朗伊俯冲带延伸长约500 km，位于新西兰北岛外海以东，在此区域太平洋板块向西俯冲至澳大利亚板块之下[14]。本文选定的研究区位于希库朗伊俯冲带北缘，在此区域板块汇聚速率在52 mm/a左右[17,28]。即将俯冲的洋壳在地形上表现为10～15 km厚的一个希库朗伊高原，由海山、熔岩流及白垩纪沉积物盆地组成，为一个白垩纪的大火成岩省[37]。希库朗伊边缘是全球目前规模较大的大火成岩省正在发生俯冲的独特活动大陆边缘。希库朗伊高原的上覆地层主要由晚白垩纪的火山碎屑岩、古新世—上新世半深海碳酸盐和更新世的硅质碎屑岩海沟沉积三部分构成[17,31,33]。希库朗伊俯冲带的构造特征沿走向上存在显著的差异，其增生楔宽度从北向南逐渐变宽，并且南缘具有较宽（30～70 km）的增生楔，而北缘则为较窄（＜20 km）的增生楔[38]。此外，俯冲的希库朗伊高原洋壳厚度从北部的10 km向南逐渐增厚，在查塔姆隆起（Chatham Rise）厚度达15 km；而上覆的沉积物厚度从北部的约1 km向南部增加到＞5 km[18,34]。
希库朗伊俯冲带地震监测表明北缘存在有深度较浅（＜15 km）、每1～2年重复出现、持续时间较短（1～3周）的慢滑移事件[17,22]。海底测地学监测表明传播至海底以下2 km深度范围内的慢滑移事件也可能会出现在海沟附近[29]。历史记录表明北缘的板块界面逆冲断层可能发生了1947年的两个大地震事件，其矩震级（MW）范围为7.0～7.2，并造成了北岛沿岸发生高度8～10 m的海啸[17,26]。国际大洋发现计划（IODP）372和375航次在希库朗伊俯冲带北缘开展了一系列钻探工作，包括在U1518站位对增生楔内部发育的一个主要的、高活动性的分支断层Pāpaku断层的科学钻探工作。地震剖面显示Pāpaku断层为发源于板块边界的巨型分支断层之一，断层滑移量约为6 km，显著大于其他逆冲断层（图3）。钻探显示其包含了一个18 m厚的主断层带和一个约10 m厚的次级断层带。断层带内部包含断层角砾，小褶皱和非对称的岩屑，被肢解的、相对均匀的、淤泥质的半深海沉积物[31,39-40]。
![]()
Hikurangi北缘05CM-04地震反射剖面显示，俯冲海山上覆板块发生显著的变形，楔体后缘的缩短主要由内部的叠瓦逆冲断层及起源于板块边界的分支断层所吸收。俯冲海山上方发育的巨型分支断层（Pāpaku断层）则与下方的板块边界逆冲断层组成了一组坡-坪式“双重构造”（thrust duplex）（图3）。此外，地震反射剖面显示在深部俯冲海山的前缘还发育一个约 3 km 厚的高振幅反射带（HRZ）[17,24]。
2. 实验方法和设置
2.1 离散元方法原理
离散元方法（discrete element method，DEM）的基本原理是通过给构建的弹性颗粒系统施加外力后，应用时间-位移有限差分方法，计算弹性-摩擦接触颗粒在牛顿第二运动定律下的移动。颗粒间的接触关系遵循Hertz-Mindlin接触准则[41]。不同于基于相似性原理的砂箱物理模拟方法，数值模拟方法对系统内单一变量的精确控制和应用，往往具备可操作性强和重复性高等特点，并且可以获取更多的系统内部信息，如应力与应变等。因此可以通过调整材料的力学参数和设定边界条件等方式模拟与自然界相近的断裂构造变形特征[42-44]。
2.2 初始模型与参数设置
本次开展的模拟实验方案设计主要参考希库朗伊北缘05CM-04地震反射剖面（图3）的构造解释方案及地层序列。初始模型设置为一个长100 km、高15 km、倾角为2.8°的楔形结构模型。楔体底部倾角的设置依据主要源于05CM-04剖面中增生楔底部的平均倾角测量。模型内设置了一套最厚5 km的沉积地层单元，1.5 km厚的洋壳基底，其间设置一套厚约250 m的韧性滑脱层。板块边界岩石随着与俯冲洋壳的相互作用，碳酸盐岩及火山碎屑沉积岩较多蚀变为蒙脱石黏土，形成韧性的滑脱层[17,45]。尽管目前没有钻穿增生楔底部的滑脱层，但对于在即将俯冲板块中钻探的原型滑脱层（proto-detachment）的钻探揭示表明其对应碳酸盐和火山碎屑岩的界面。本文以希库朗伊北缘05CM-04剖面存在的“平顶山型”海山为例，将海山的宽度设置为15 km，高度设置为1 km[24,46]。同时参考希库朗伊北缘观测到的俯冲海山规模及现今剖面的平衡剖面恢复结果，将俯冲海山设置于剖面44 km处（图4B）[17-18,47]。
![]()
图 4 本文采用的离散元数值模拟实验模型设置A: 光滑海底模型，B: 海山俯冲模型；箭头表示俯冲基底的相对运动方向。Figure 4. The model setups of discrete-element numerical simulationA: Smooth seafloor model, B: seamount subduction model; the arrow indicates the relative direction of motion.模型中颗粒的力学参数设置参考Morgan的离散元实验设置[36]，其中上覆地层的颗粒密度设置为2500 kg/m3, 滑脱层内颗粒密度为2300 kg/m3。颗粒的泊松比为0.2，剪切模量为2.9 GPa，局部阻尼系数为0.4，时间步长为50 ms，重力加速度为9.8 m/s2。模型中颗粒黏结参数设置参考表1，地层的黏聚力设置为24.4 MPa，内摩擦角φ为22.5°，颗粒摩擦系数选取0.3。本文参考Morgan 和 Bang[44]中的模型设置，将左侧挡板以2 m/s的速度向右移动挤压颗粒，而海山相对于洋壳基底固定，用以模拟海山相对向左俯冲于增生楔之下。对于冲断增生楔而言，消减模型和推挤模型尽管在剖面上具有相反的速度场特征，但在平面和整体上具有一致的速度场和应力状态[48]。其次考虑到板块边界摩擦系数随深度的变化，将深部下倾的沿板块边界20 km长度的滑脱层的摩擦系数设置为0.3，其余上倾部分滑脱层的摩擦系数设置为0.1，模型的总缩短量设置为50 km。
单元 宏观参数 岩层颗粒间的黏结参数 内聚力
/MPa内摩擦角
/(°)摩擦系数μ 弹性模量/MPa 剪切模量/MPa 抗拉强度/MPa 剪切强度/MPa μ1 24.4 22.5 0.3 2.0×108 2.0×108 3.0×107 6.0×107 μ2 24.4 22.5 0.1 2.0×108 2.0×108 3.0×107 6.0×107 上覆地层 24.4 22.5 0.3 2.0×108 2.0×108 3.0×107 6.0×107 光滑海底模型中颗粒的力学性质与海山俯冲模型一致，作为参考模型对比使用（图4A）。通过对两组实验中得到的构造变形、应变及应力结果进行对比，进一步讨论有/无海山模型对增生楔形态及内部断层演化的影响。将颗粒变形与应变的结果相叠加，得到实验模拟的构造变形结果如图5所示。变形应变的定义参考Morgan（2015），其中色标中的红色和蓝色分别代表顺时针剪切和逆时针剪切，颜色深浅表示剪切应变的大小。
![]()
图 5 楔体构造变形与应变的演化左：光滑海底模型，右：海山俯冲模型。Figure 5. The experimental evolution of structural deformation and strain in the wedgesLeft: smooth seafloor model, right: seamount subduction model.3. 实验结果
3.1 光滑海底模型
光滑海底俯冲模型实验过程中增生楔的构造变形与应变如图5所示，左侧颗粒不断向右运动，当模型缩短量为2、5、11 km时，上覆地层发生缩短变形并依次发育F1、F2、F3断层。在挤压的初始阶段（＜10 km），楔体表面坡角（α）快速增大（＞20°），当缩短量为11 km时，楔体前缘F3断层的出现使表面坡角逐渐减小。在缩短量为18、23、27、32 km时，楔体内部依次形成了F4—F7断层，此阶段楔体的表面坡角无明显变化并趋于稳定（13°±1°）。当缩短量≥36 km时，楔体表面坡角稳定在12°（图6A）。F1—F5断层具有相似的演化特征，挤压初期断层位移量快速增加，后期趋于稳定（图7A）。通过对楔体内发育的断层位移量进行统计发现，楔体内部断层间距相近，且位移量随着缩短量增加而逐渐增大，总体符合前展式逆冲扩展顺序（图5、6）。由图6B可见，楔体宽度表现出显著的周期性变化特征，即在俯冲过程中楔体宽度先是快速增加而后表现出缓慢的缩短过程，楔体内部的缩短由新发育的断层所吸收。
![]()
图 6 楔体表面坡角（A）与楔体宽度（B）随缩短量变化统计图Figure 6. Statistical graphs of surface slope (A) and wedge width (B) versus shortening in seamount versus smooth seafloor models![]()
图 7 楔体内部断层位移量与缩短量统计A: 光滑海底模型，B: 海山俯冲模型。Figure 7. Statistical graphs of fault displacements versus shortening in wedgesA: Smooth seafloor model, B: seamount subduction model.3.2 海山俯冲模型
根据增生楔的形态和断层位移量特征，本文将DEM模拟的海山俯冲划分为3个阶段：I海山未接触阶段，对应缩短量0～17 km；II海山与增生楔相互作用阶段，对应缩短量17～36 km；III海山俯冲效应减弱，增生楔重新向前扩展阶段，对应缩短量36～50 km。
在第I阶段，在海山俯冲模型中的增生楔的总体演化特征与光滑海底模型相似，楔体中同样发育两条前展式的逆冲断层（F1和F2）并对应较大的楔体表面坡角（20°～40°）（图6A）。在此阶段，两个对比模型的楔体宽度变化趋势也基本一致，并随着缩短量的增加而逐渐变宽至12 km左右。当缩短量为12～16 km时，F3断层的发育导致海山模型的楔体宽度大于光滑模型的楔体宽度（图6B），而表面坡角相近并呈现出相似的变化趋势。
在第II阶段，当缩短量增至17 km时，模拟显示俯冲的海山顶部开始发育一条逆冲断层（F4），其作为新的板块间滑脱层沿着海山向前扩展，并最终突破至海底作大位移量（最大约10 km）、低倾角（近水平）的逆冲推覆。受海山俯冲的影响，滑脱层上移导致楔体底部倾角增加，而楔体的表面坡角则从17°骤减至9°左右。在第II阶段初期，楔体宽度从25 km快速增加至35 km，而随后则不断发生内部缩短；而同阶段的光滑模型中，楔体的宽度则经历周期性的增长（最大约37 km）。模拟显示在第II阶段，除F4之外，楔体内不再发育新的断层，表面坡角则基本稳定（8°～9°），但总体上低于同阶段的光滑海底俯冲模型（12°～13°）（图6A）。断层位移量统计显示，海山模型中的F4断层的位移量约为光滑模型中F4断层的3倍，且活动时间也远超过楔体中的其他断层（图7）。
在第III阶段，当缩短量大于36 km时，滑脱层沿着海山后缘下移并快速向前扩展，发育前缘逆冲断层系。增生楔前缘依次发育3条逆冲断层（F5、F6、F7），形态上组成了一个叠瓦状逆冲构造（图5、7）。在此阶段，海山模型中楔体的表面坡角和宽度显著增大并反超光滑模型。断层位移量统计显示，新形成的逆冲断层（F5、F6、F7）相较于光滑模型具有较小的位移量，并且其形成时间相对延迟（图7B）。平顶山型海山的俯冲导致板间滑脱层的上移，形成的逆冲断层（F4）沿着海底面向前发生低角度、大位移量的逆冲推覆，保护了海山后缘的局部地层遭受构造变形，进而抑制了海山后缘逆冲断层的发育。
本文还对比研究了光滑海底与海山俯冲两种模型在变形过程中应力状态的变化特征。图8展示了模型缩短量为28 km时，楔体内部的最大剪切应力的分布状况，其中左侧为光滑海底模型，右侧为海山俯冲模型。模型结果对比显示：相较于光滑海底俯冲模型，海山俯冲模型内部表现出最大剪切应力分布较强的非均质性特征。海山俯冲导致在其前缘区域形成较强的最大剪切应力分布，而海山后缘则形成应力较小的阴影区。
![]()
图 8 光滑海底（A）与海山俯冲（B）模型中楔体最大剪切应力（τmax）分布对比Figure 8. Results of maximum shear stress in the smooth seafloor (A) and seamount (B) subduction models4. 讨论
4.1 海山俯冲与巨型分支断层
巨型分支断层起源于板块边界巨型逆冲断层（megathrust），一直连续延伸至海底面，具有较长的断层长度，是一种具有较大位移量的 “反序型”逆冲断层（out-of-sequence thrust, OOST）[49]。巨型分支断层广泛发育在世界主要俯冲带，如日本南海（Nankai）、阿拉斯加（Alaska）和希库朗伊（Hikurangi）等[12,26,50-53]。由于其切穿了增生楔内部并且具有较大的断层面积，可以在同震期间造成海床大范围的快速垂直隆升，因此被认为对海啸的发生扮演重要的角色。但是，关于巨型分支断层的起源、断层活动性演化及其与海山俯冲之间的关系，目前还尚不完全明确。
IODP大洋钻探结合地震反射剖面和磁性地层定年的研究，揭示了日本南海海槽巨型分支断层与板间滑脱层之间的活动性交替演化规律：巨型分支断层在约1.95 Ma启动，在高活动性阶段之后，断层滑移向下沿着滑脱层向前，巨型分支断层的活动性下降；在约1.55 Ma以来，该分支断层重新活动和抬升[12,51]。希库朗伊俯冲带的IODP钻孔U1518钻穿了巨型分支断层Pāpaku断层，岩芯构造分析表明断层上盘主要为脆性结构特征，而下盘表现出韧性和脆性变形随深度降低的特征。同时，孔隙水的地球化学剖面也显示断层之下存在近海底面成岩序列的重复，表明下盘地层未经历明显埋藏过程[30-31]。
本文的研究表明，上述讨论的巨型分支断层的断层岩结构特征、断层活动性、孔隙水地球化学特征，可以用海山俯冲的模型来建立合理的解释方案。正如DEM模拟结果所展示的，巨型分支断层起源于板间滑脱层与海山的作用阶段，由滑脱层快速沿海山上移演变而成，断层最终延伸至海底，并作低角度、大位移量的逆冲推覆（图9B、C）。这就合理解释了希库朗伊俯冲带巨型分支断层的上盘脆性变形与下盘韧性变形的上下配置关系。结合下盘保存的成岩剖面，最新的研究推测巨型分支断层在50～100 ka沿着海底面发生大规模的逆冲推覆，对下盘地层造成“泥岩涂抹”（clay smear）的封堵效应，限制了下盘的成岩-压实作用[54]。而关于巨型分支断层的活动性问题，DEM模拟揭示了其活动性的快速增加与海山俯冲直接相关。随着海山效应的减弱，滑脱层下移沿海山后缘向前扩展形成叠瓦状逆冲构造，而巨型分支断层则最终停止活动。在实际的俯冲带中，多个海山的俯冲可以用来解释巨型分支断层活动性的周期性演化规律。因此，本文开展的DEM模拟实验支持在新西兰希库朗伊俯冲带和日本南海海槽中发育的巨型分支断层的形成及其演化主要受海山俯冲控制。
![]()
图 9 海山俯冲过程中的增生楔内部断裂演化的三阶段模式图A: 海山与增生楔作用初始阶段，滑脱层沿海山顶面上移；B: 海山与增生楔持续作用，发生沿海底的远距离推覆作用；C: 海山俯冲后期，滑脱层下移并向前传播形成前缘逆冲断层系。Figure 9. Three-stage model of faults evolution related to seamount subductionA: The initial stage of seamount-accretionary wedge interaction, the detachment moves up along the top of seamount; B: Seamounts and accretionary wedges continue to interact, and long-distance, low-angle thrusting along the seafloor occurs; C: In the late stage of seamount subduction, the detachment moves down and propagates forward to form the frontal thrusts.4.2 希库朗伊俯冲带“低速”与“高速”海山模型的争议
基于地震反射和磁法勘探数据，前人推测了希库朗伊俯冲带深部的海山位置和形态，其位于现今剖面增生楔底部滑脱层变缓部位之下（图2A）[24,46]。对于这一海山的位置，三维地震速度模型仅获取了相对低速的（VP=4.0～5.0 km/s）结果，并没有探测到任何随深度增加的速度现象。对于这一差异，他们提出两种可能的海山模型：① “低速海山”俯冲模型；② “高速海山”俯冲模型。对于第一种低速海山模型，如果Barker 等[24]模型中海山的位置是正确的话，那么正在现今俯冲的海山内部必定具有很强的速度非均质性并且大部分区域是低速性质的。而第二种“高速海山”模型中，海山位于“低速海山”模型沿滑脱层下倾方向约10 km，内部速度超过5.0 km/s（图2B）。这与之前地震反射剖面所揭示的形成于海山前缘高振幅反射带（high-amplitude reflectivity zone, HRZ）所解释的低地震波速、富含流体的沉积物的结论截然相反[24,34]。
Sun等[27]利用有限元数值模拟方法，模拟了汇聚边界海山俯冲对板间逆冲断层应力及滑移的力学行为及水文地质效应。在俯冲海山对上覆板块的碰撞破坏下，海山前缘广泛发育的断层及断裂作用致使其下倾端沉积层严重脱水且过度固结。沉积物孔隙度及合成剖面VP的经验关系表明海山前缘出现了异常于光滑模型的VP急剧增加的现象，且较低的孔隙度及硬度增强的围岩对板块边界逆冲断层的剪切强度具有增强作用。Clavert在对卡斯卡迪亚北部俯冲带深部地震资料解释也发现深部沉积层的过度固结脱水可呈现高振幅反射特征，且深部高反射带上的剪切作用形成了一个双重构造[55]。
Hikurangi北缘记录到的2002—2008年期间的慢滑移事件在位置上与高振幅反射区密切相关，一般认为高振幅反射区的流体主要来源于层间变质岩及俯冲沉积物的带入，或俯冲到深部的压实脱水作用[46,56]。显然，低速流体作用并不能很好地解释该区域的高VP波速度模型。因此，结合数值模拟实验的结果，本次研究显示在海山持续向深部俯冲过程中，海山前缘与增生楔底部存在着长期剪切滑动，最终形成了一个具有双重构造（duplex）特征的剪切带（图10B、C）。DEM数值模拟展示了海山前缘双重构造剪切带不断加宽的过程，最终的结果（图10B、C）与实际地震反射剖面的解释方案具有较好的一致性。在该双重构造带区域，持续的剪切作用可能导致海山前缘的沉积物发生严重脱水和过度固结，因此导致其VP速度显著增加。而海山顶部发育的低地震波速区域，可能是海山在俯冲过程中与增生楔底部碰撞作用使其发生破碎化，流体沿断裂进入内部所导致的[30]。因此，结合电磁反演的孔隙度分布特征，本文的DEM数值模拟研究支持希库朗伊俯冲带“低速”海山的模型。
![]()
图 10 DEM数值模拟结果与希库朗伊俯冲带构造剖面对比A: 05CM-04构造剖面图（修改自文献[17]），深度放大2倍， B: 海山俯冲DEM模拟的应变分布，C: 与海山俯冲相关构造解析模型。Figure 10. Comparison between DEM simulation results and structural profile in Hikurangi subduction zoneA: Structural interpretation of seismic profile 05CM-04 (modified from reference[17]), B: strain derived from DEM simulation related to seamount subduction; C: structural model related to seamount subduction.5. 结论
（1）海山俯冲改变了增生楔内部逆冲断层的应变分配，形成的巨型分支断层沿海底面发生低角度、远距离逆冲推覆，这一现象可以用来解释岩芯揭示的上盘脆性变形与下盘韧性变形的二元结构特征。
（2）海山俯冲导致在增生楔内部应力分布的非均质性，在其前缘形成较高的应力聚集区并形成双重构造特征的剪切带，而在其后缘则形成应力阴影区。
（3）海山俯冲前缘形成的双重构造剪切带有利于沉积物发生脱水和固结，因此该区域并非富含流体的低地震波速区域。本文的模拟结果结合电磁探测的证据支持希库朗伊俯冲带的“低速海山”俯冲模型。
致谢：本文的数值计算是在南京大学高性能计算中心的计算集群上完成的, 数值模拟实验使用东华理工大学李长圣博士研发的离散元数值模拟软件完成。文中采用的应变计算代码修改自莱斯大学Julia K Morgan和Thomas Fournier的脚本, 在此表示感谢。
-
![]()
图 2 平湖组烃源岩有机质丰度评价
Figure 2. Assessment on organic matter abundance of Pinghu Formation
![]()
图 4 平湖组某井煤岩壳质组有机显微组分
a.树脂体，油浸荧光，500×；b.孢子体，油浸荧光，500×；c. 层状藻类体，油浸荧光，500×；d. 角质体，油浸荧光，500×。
Figure 4. Organic micro-compositions of exinite of coal rocks in Pinghu Formation
a.resinite，immersion fluorescence，500×；b. sporophyte，immersion fluorescence，500×；c. lamellar algae，immersion fluorescence，500×；d. cutinite，immersion fluorescence，500×.
![]()
图 5 平北地区平湖组烃源岩热演化剖面
Figure 5. Thermal evolution profile of source rock of Pinghu Formation in Pingbei area
![]()
图 6 平湖组烃源岩热模拟试验产率与热演化剖面
Figure 6. Product yields and thermal evolution profile of the thermal simulation experiment in Pinghu Formation
![]()
图 7 平湖组煤岩、暗色泥岩以及原油的碳同位素特征
Figure 7. Character of carbon isotope of coal rocks, mudstone and oils in Pinghu Formation
![]()
图 8 B-4井烃源岩及原油生物标志物特征
Figure 8. Character of source rocks and oil biomarkers of B-4 Well
表 1 研究区煤岩有机质丰度评价标准
Table 1 Evaluation criteria of organic matter abundance of coal rocks in study area
评价标准 S1+S2/（mg/g） HI/（mg/g） A/% 好 ＞150 ＞300 ＞2.0 中 90～150 150～300 1.0～2.0 差 60～90 100～150 0.5～1.0 非 ＜60 ＜100 ＜0.5 
下载: 导出CSV
表 2 研究区泥岩有机质丰度评价标准
Table 2 Evaluation critera of organic matter abundance of mudstone in study area.
评价标准 TOC/% S1+S2/（mg/g） A/% HC/10-6 好 ＞2.5 ＞6.0 ＞0.15 ＞400 中 1.0～2.5 2.0～6.0 0.05～0.15 150～400 差 0.6～1.0 0.6～2.0 0.02～0.05 50～150 非 ＜0.6 ＜0.6 ＜0.02 ＜50
下载: 导出CSV
表 3 平湖组煤及泥岩热模拟试验样品地化参数
Table 3 Geochemical index of the thermal simulation experiment of coal rock and mudstone in Pinghu Formation
井号 岩性 Ro/% TOC/% HI/（mg/g） S1/（mg/g） S2/（mg/g） A-X井 煤 0.4 61 320 2.34 127.44 B-X井 泥岩 0.6 2.5 236 0.54 5.86
下载: 导出CSV
-
[1] 史训知, 戴金星, 王则民, 等. 联邦德国煤成气的甲烷碳同位素研究和对我们的启示[J]. 天然气工业, 1985, 5(2):1-9. [SHI Xunzhi, DAI Jinxing, WANG Zemin, et al. A study of methane carbon isotope of coal-formed gas in FRG and its inspiration to us [J]. Natural Gas Industry, 1985, 5(2): 1-9. [2] Brooks J D, Smith J W. The diagenesis of plant lipids during the formation of coal, petroleum and natural gas-I. Changes in the n-paraffin hydrocarbons [J]. Geochimica et Cosmochimica Acta, 1967, 31(12): 2389-2397.
[3] 戴金星. 煤成气及鉴别理论研究进展[J]. 科学通报, 2018, 63(14):1290-1305. [DAI Jinxing. Coal-derived gas theory and its discrimination [J]. Chinese Science Bulletin, 2018, 63(14): 1290-1305. [4] 黄第藩. 成烃理论的发展-(Ⅱ)煤成油及其初次运移模式[J]. 地球科学进展, 1996, 11(5):432-438. [HUANG Difan. Advances in hydrocarbon generation theory-(Ⅱ) oils from coal and its primary migration model [J]. Advance in Earth Sciences, 1996, 11(5): 432-438. doi: 10.3321/j.issn:1001-8166.1996.05.002 [5] Huang D F, Li J C, Zhang D J. Immature oil in Chinese nonmarine sediments and significance[M]//Organic Geochemistry Advances and Application in Energy and the Natural Environment. Manchester: Manchester University Press, 1991: 39-47.
[6] 黄第藩, 卢双舫. 煤成油地球化学研究现状与展望[J]. 地学前缘, 1999, 6(S1):183-194. [HUANG Difan, LU Shuangfang. Current situation and prospect in geochemistry of oil from coals [J]. Earth Science Frontiers, 1999, 6(S1): 183-194. [7] 周洁. 东海盆地西湖凹陷始新统平湖组煤系烃源岩地球化学特征[D]. 浙江大学硕士学位论文, 2012: 1-51. ZHOU Jie. Geochemical characteristics of Eocene Pinghu formation coal-bearing source rocks in the Xihu depression, East China Sea Basin[D]. Master Dissertation of Zhejiang University, 2012: 1-51.
[8] 钱门辉. 西湖凹陷煤系烃源岩生烃特征研究[D]. 中国地质大学(北京)硕士学位论文, 2010: 1-99. QIAN Menhui. A model research of generation of hydrocarbon coal source rock in Xihu depression[D]. Master Dissertation of China University of Geosciences (Beijing), 2010: 1-99.
[9] 孙伯强, 陈剑新. 东海西湖凹陷保俶斜坡下第三系煤成油的研究[J]. 海洋石油, 2000(1):1-6. [SUN Boqiang, CHEN Jianxin. Preliminary study on coal-generating oil of Eogene system in the Baochu slope of Xihu sag, East China Sea [J]. Offshore Oil, 2000(1): 1-6. [10] 魏恒飞, 陈践发, 陈晓东, 等. 西湖凹陷平湖组滨海型煤系烃源岩发育环境及其控制因素[J]. 中国地质, 2013, 40(2):487-497. [WEI Hengfei, CHEN Jianfa, CHEN Xiaodong, et al. The controlling factors and sedimentary environment for developing coastal coal-bearing source rock of Pinghu Formation in Xihu depression [J]. Geology in China, 2013, 40(2): 487-497. doi: 10.3969/j.issn.1000-3657.2013.02.013 [11] 陈建平, 赵长毅, 何忠华. 煤系有机质生烃潜力评价标准探讨[J]. 石油勘探与开发, 1997, 24(1):1-5. [CHEN Jianping, ZHAO Changyi, HE Zhonghua. Criteria for evaluating the hydrocarbon generating potential of organic matter in coal measures [J]. Petroleum Exploration and Development, 1997, 24(1): 1-5. [12] 黄第藩. 中国陆相油气生成[M]. 北京: 石油工业出版社, 1982: 1-355. HUANG Difan. Oil and Gas Generation in China Continental Deposit[M]. Beijing: China Petroleum Industry Press, 1982: 1-355.
[13] 徐春华, 樊庆真, 孔凡顺, 等. 煤系烃源岩评价中的误区[J]. 油气地质与采收率, 2001, 8(1):78-80. [XU Chunhua, FAN Qingzhen, KONG Fanshun, et al. Misunderstanding on evaluation of coal-measure hydrocarbon source rock [J]. Petroleum Geology and Recovery Efficiency, 2001, 8(1): 78-80. doi: 10.3969/j.issn.1009-9603.2001.01.023 [14] 秦建中. 中国烃源岩[M]. 北京: 科学出版社, 2005: 1-618. QIN Jianzhong. Source Rocks in China[M]. Beijing: Science Press, 2005: 1-618.
[15] 王飞宇, 傅家谟, 刘德汉. 煤和陆源有机质生油岩有机岩石学特点及评价[J]. 石油勘探与开发, 1994, 21(4):30-35. [WANG Feiyu, FU Jiamo, LIU Dehan. Organic petrological characteristics of coal and terrestrical organic matter and their assessment as a oil source rock [J]. Petroleum Exploration and Development, 1994, 21(4): 30-35. [16] 刘德汉, 傅家谟, 肖贤明. 煤成烃的成因与评价[J]. 石油勘探与开发, 2005, 32(4):137-141. [LIU Dehan, FU Jiamo, XIAO Xianming. Origin and appraisal of coal derived gas and oil [J]. Petroleum Exploration and Development, 2005, 32(4): 137-141. [17] 郭黔杰. 煤作为源岩研究时对显微组分的重新认识[J]. 石油勘探与开发, 1993, 20(2):37-40. [GUO Qianjie. Macerals reconsideration in the studies of “oil from coal” [J]. Petroleum Exploration and Development, 1993, 20(2): 37-40. [18] 杨鹏程, 余学兵, 张传运. 涠西南凹陷D次凹流二段优质烃源岩特征[J]. 海洋石油, 2017, 37(1):1-6, 79. [YANG Pengcheng, YU Xuebing, ZHANG Chuanyun. Characteristics of high-quality source rock in Member Liu-2 of D sub-sag in Wexinan sag [J]. Offshore Oil, 2017, 37(1): 1-6, 79. doi: 10.3969/j.issn.1008-2336.2017.01.001 [19] 张景廉. 中国侏罗系煤成油质疑[J]. 新疆石油地质, 2001, 22(1):1-8. [ZHANG Jinglian. On the query of coal-formed oil of Jurassic in China [J]. Xinjiang Petroleum Geology, 2001, 22(1): 1-8. doi: 10.3969/j.issn.1001-3873.2001.01.001 [20] 苏传国, 朱建国, 孟旺才, 等. 吐哈盆地“煤成油”形成机制探讨[J]. 吐哈油气, 2005, 10(1):14-20. [SU Chuanguo, ZHU Jianguo, MENG Wangcai, et al. Discussion on generation mechanism of coal-formed oil in Turpan-Hami Basin [J]. Tuha Oil & Gas, 2005, 10(1): 14-20. [21] 傅宁. 东海盆地西湖凹陷煤系烃源岩及凝析油中的二萜化合物[J]. 中国海上油气(地质), 1994, 8(1):21-28. [FU Ning. Diterpenoid compounds in coal measures and condensates in Xihu sag of East China Sea Basin [J]. China Offshore Oil and Gas (Geology), 1994, 8(1): 21-28. [22] 戴卿林, 郝石生, 卢双舫, 等. 煤成油油源对比问题讨论[J]. 地球化学, 1996, 25(4):324-330. [DAI Qinglin, HAO Shisheng, LU Shuangfang, et al. Discussions about the correlation of oils from coals and their sources [J]. Geochimica, 1996, 25(4): 324-330. -
期刊类型引用(4)
1. 李朗，黄晓燕，何伟，秦甜甜. 水文地质单元交界区潜水水化学特征及演化分析：以泰州市区为例. 科学技术与工程. 2024(06): 2265-2275 . 
百度学术
2. 甘恢元. 江淮东部史前文化格局与相关问题探讨. 南方文物. 2023(01): 119-131 . 百度学术
3. 李保华，张楷，徐烨，范代读. 大丰野鹿荡中全新世以来有孔虫群落及古环境变化. 微体古生物学报. 2022(04): 339-347 . 百度学术
4. 李宗盟，朱文敏，高红山，刘芬良，邢伟. 黄淮平原区晚新生代气候变迁. 海洋地质与第四纪地质. 2021(04): 179-191 . 本站查看
其他类型引用(2)
计量
- 文章访问数: 2653
- HTML全文浏览量: 431
- PDF下载量: 38
- 被引次数: 6
