Two-dimensional discrete element simulation of plate subduction deformation process: An insight into the genesis of East China Sea Shelf Basin
-
摘要: 俯冲变形作用是板块汇聚过程中存在的构造地质现象，是当前构造地质研究的热点。目前对板块俯冲变形的研究尚不完善，俯冲角度变化对变形过程造成的影响还需进一步研究。位于欧亚板块东南部的中国东海大陆架盆地的构造演化特征及动力学机制与菲律宾板块向欧亚板块的俯冲作用有关。利用构造地质领域中新兴的离散元模拟方法，通过构建离散元模型模拟研究板块俯冲变形演化过程，并将实验结果与菲律宾板块向中国东海俯冲部位地层相比较，结果表明：① 板块俯冲变形特征与俯冲角度有关，俯冲角度不同，其最终形成的变形样式也不同；② 断层数目随着俯冲角度的减缓而增加，断层所扩展的水平距离随俯冲角度的减缓而增大，且不同俯冲角度下相同位置所形成的断距不同；③ 俯冲楔高度随着俯冲角度的减缓而增大，地壳变形幅度越大，且最终形成俯冲楔形态类型不同；④ 反冲断层形成时间随着俯冲角度的减缓而越来越晚；⑤ 实验模拟结果与实例具有相似的构造特征。研究结果解析了不同俯冲角度下板块俯冲变形的演化过程，有助于对板块汇聚过程中俯冲变形作用的进一步认识。Abstract: Subduction deformation is a tectonic phenomenon in the process of plates convergence, which remains a hot topic. The East China Sea Shelf Basin is located in the south-eastern border of the Eurasian plate, whose tectonic evolution and dynamic mechanism are related to the subduction of the Philippine plate and the Pacific plate beneath the Eurasian plate. At present, the research on the subduction deformation is still not so perfect, and the influence of subduction angle changes on the deformation process needs further study. The authors, in this paper, tried to use the method of discrete element simulation, a kind of emerging method in the field of tectonics, to simulate the evolution process of plate subduction deformation. By constructing discrete element models and comparing the experimental results with the strata in the subduction zone between the Philippine plate and the East China Sea, the research shows that (1) The subduction deformation characteristics are related to the subduction angle. Different subduction angles result in different deformation patterns. (2) The number of faults and the horizontal offsets of faults increase as the subduction angle decreases. Fault displacement at the same position under different subduction angles is different. (3) The height of the subduction wedge and deformation degree increases as the subduction angle decreases. And the shape of subduction wedge is different. (4) Back-thrusts form later as the subduction angle decreases. (5) Experimental simulation results have similar structural characteristics with instance. The experimental simulation results may explain the evolution process of plate subduction deformation under different subduction angles, which is helpful to further understand the subduction deformation during plates convergence process.
-
地震绕射波是地下构造细节和异常的反映，地下构造的不连续点会产生绕射波，如裂缝、断层、溶洞和地层尖灭等。在原始地震勘探数据中，既存在绕射波也存在反射波，反射波是地质背景的综合反映，绕射波则是地质细节的小构造的反映。相对反射波而言，绕射波的能量往往很弱，甚至难以识别，导致绕射信息常常被高能量的反射信息所覆盖，而绕射波信息是提高地震勘探分辨率的重要信息载体。在地震资料处理解释中，正确识别断层、尖灭和小尺寸散射体等地质不连续性是一个重要问题，绕射波包含有关产生绕射波的介质结构和组成的有价值信息。然而，在标准地震资料处理中，绕射波常被视为噪声而被抑制[1-3]。把绕射波从原始地震勘探记录中分离出来进行单独成像，可以加强小型地震异常体的成像，从而提高复杂断块的地震成像精度。
在过去的几十年，许多地球物理学者对地震绕射波分离方法做了相关研究。1988年Kanasewich和Phadke[4]提出了一种适用于共偏移距道集的绕射波走时曲线直接拾取绕射波的方法来分离绕射波。1992年Claerbout[5]首次提出了平面波预测滤波器的概念，并在信噪分离方面得到了很好的应用。Fomel[6]进一步完善了平面波预测滤波器，并提出了平面预测滤波器的残差可以作为绕射波的近似。Nowak等[7]根据反射波和绕射波同相轴的几何差异，采用加权双曲Radon滤波器分离地震绕射波。Berkovitch等[8]利用绕射波多次聚焦叠加剖面将绕射波聚焦到绕射点位置，同时将反射波反聚焦到整个道集记录，从而分离绕射波。Koren和Ravve[9]提出了方位相关反射波属性，通过构建加权叠加滤波器将绕射波和反射波分离。Klokov等[10]根据角度域共成像点道集中反射波和绕射波的差异，首先压制去除反射顶点附近反射能量，然后利用混合Radon变换去除剩余反射波能量，进而分离出绕射波。Moser和Howard[11]提出基于深度域的反射聚焦，从叠前数据中过滤反射进而对绕射波分离成像的方法。Decker等[12]使用人工模拟数据和从皮森盆地获得的野外实际数据进行绕射波分离，评估解决绕射特征的多种方法，证明平面波预测滤波器方法能使图像保存更多绕射能量和较少的迁移伪影，但不足是无法消除斜率不连续的反射变量。上述绕射波分离算法主要是利用了绕射波和反射波在不同地震数据道集中的差异，通过数学变换或者信号预测的方法将绕射波识别和分离出来。在国内研究方面，赵娟娟等[13]利用F-K滤波方法分离成像绕射波，通过F-K域斜率滤波去除反射波从而分离出绕射波信息。马永生[14]从油藏中具有代表性的地质模型出发，分析VSP中绕射波的基本特点，理论上总结了绕射波运动学的部分认识规律。黄建平等[15]着重介绍了平面波预测滤波器及其在叠前域、叠后域分离并成像绕射波的方法原理，探讨了绕射波分离成像方法的应用前景和改进方向。朱生旺[16]提出一种局部倾角滤波和预测反演联合起来分离绕射波的方法。蒋波等[17]提出一种基于反射波层拉平的绕射波分离与成像方法用于碳酸盐岩缝洞储集体地质目标，提高对溶洞、裂缝等特殊异常体的刻画精度。李正伟等[18]根据倾角域共成像点道集中菲涅耳带变化特点，通过倾角-偏移距道集中精确切除菲涅耳带的方法来分离绕射波。本文在平面波分解绕射波分离[6]基础上，分析了绕射波与反射波信号特征的差异，深入研究了平面波预测的绕射波分离方法的影响因素，并探讨了噪声和平面波预测滤波器平滑参数对地震绕射波分离结果的影响，为在实际地震数据处理中应用平面波分解方法进行绕射波分离提供了理论指导。
1. 平面波预测绕射波分离
地震记录既含有反射波同相轴也含有绕射波同相轴，在点震源激发时，反射波和绕射波波形在炮记录中都是双曲线，不易于分离反射波和绕射波。在平面波震源激发时，来自反射界面的反射波是平面波，而来自绕射体的绕射波是双曲线。因此，在平面波震源条件下反射波和绕射波具有明显的几何差异，更有利于绕射波分离。在实际地震勘探中大都是点源激发，反射波和绕射波的同相轴不易区分，但我们根据爆炸反射面理论将共偏移距剖面近似为平面波震源获得的地震记录，这样绕射波和反射波有较好的几何差异，根据其几何差异进行绕射波分离。
平面波预测滤波器利用未知的局部倾角构造最优非稳态预测误差滤波器，根据相邻地震道的数据，在局部倾角平滑变化约束下，使预测误差最小化以此预测目标道数据。通过最优化方法不断迭代可获得局部倾角场，倾角场能够反映反射波场几何信息。
局部平面波的微分方程表示为[6]：
$$\frac{{\partial P}}{{\partial {{x}}}} + \sigma \frac{{\partial P}}{{\partial t}} = 0$$ (1) 其中，
$P\left( {t,x} \right)$ 是波场，$x$ 、$t$ 分别是炮检距和旅行时，$\sigma $ 是局部倾角，是$x$ 和$t$ 的函数，局部倾角是旅行时对炮检距的导数。Fomel给出了Z变换域平面波预测的形式[6]：
$${\hat P_{x + 1}}\left( {{Z_t}} \right) = {\hat P_x}\left( {{Z_t}} \right)\frac{{B\left( {{Z_t}} \right)}}{{B\left( {{\rm{1/}}{Z_t}} \right)}}$$ (2) 式中，
${{\hat P}_x}\left( {{Z_t}} \right)$ 是相应记录道的Z变换形式，$B\left( {{Z_t}} \right){\rm{/}}B\left( {{\rm{1/}}{Z_t}} \right)$ 等价于全通数字滤波器。用Taylor级数展开拟合得到${B\left( {{Z_t}} \right)}$ 五阶中心滤波器：$$ \begin{split} {B_5}\left( {{Z_t}} \right) =& \frac{{(1 - \sigma )(2 - \sigma )(3 - \sigma )(4 - \sigma )}}{{1\;680}}Z_t^{ - 2} + \\ &{\frac{{(4 - \sigma )(2 - \sigma )(3 - \sigma )(4 + \sigma )}}{{420}}Z_t^{ - 1} + }\\ &{\frac{{(4 - \sigma )(3 - \sigma )(3 + \sigma )(4 + \sigma )}}{{280}} + }\\ &{\frac{{(4 - \sigma )(2 + \sigma )(3 + \sigma )(4 + \sigma )}}{{420}}{Z_t} + }\\ &{\frac{{(1 + \sigma )(2 + \sigma )(3 + \sigma )(4 + \sigma )}}{{1\;680}}Z_t^2} \end{split} $$ (3) 利用平面波预测滤波器预测分离绕射波的过程中，关键一步是反射波同相轴局部倾角的估计。如果局部倾角已知的情况下，可以用
$C\left( \sigma \right)$ 表示地震数据与二维滤波器的褶积算子。因此，在最小平方意义下，可通过求解如下最小二乘目标函数估计局部倾角：$$C\left( \sigma \right)d \approx 0$$ (4) 其中，
$d$ 为已知的地震数据。优化问题（3）可以利用高斯—牛顿线性化迭代的方法求解：$$C'\left( {{\sigma _0}} \right)\Delta \sigma d + C\left( {{\sigma _0}} \right)d \approx 0$$ (5) 式中，
$\Delta \sigma $ 为倾角增量，${{\sigma _0}}$ 为初始倾角，$C\left( {{\sigma _0}} \right)$ 为对应初始倾角${{\sigma _0}}$ 估计的褶积算子。$C'\left( {{\sigma _0}} \right)$ 为$C\left( {{\sigma _0}} \right)$ 的导数值。为了保证倾角场稳定可靠，在反演求解过程中引入正则化约束项，控制倾角场沿着空间和时间两个方向平滑。由于平面预测是针对平面波形的预测，而绕射波难于满足其预测方程，因此一旦估计出了倾角场，将倾角参数代入预测方程（4），预测的残差就是绕射波的近似。图1给出了平面波预测绕射波分离算法的流程图。在平面波预测绕射波分离方法中，倾角场估计的正则化平滑参数是一个重要的参数，它控制着倾角的空间平滑性，从而约束着预测的范围。同时，实际地震勘探中噪声的存在也会影响倾角的估计和绕射波的分离，因此需要对噪声和正则化参数进行深入的分析讨论。
2. 噪音水平和光滑半径影响因素
以地震勘探中被广泛采用的Sigsbee 2A模型为例对绕射波分离算法的影响因素进行分析。在绕射波分离的过程中，倾角场预测结果对绕射波分解效果影响较大，而噪声和平面波预测平滑参数均对倾角场的估计产生影响。图2-4中，在目标数据中增添了不同方差的均值为0的正态分布随机噪音，给出了不同噪声水平下的倾角场估计结果和绕射波分离结果。随着噪音水平的增大，原始数据共偏移距地震数据变得复杂，估计的倾角场准确性降低，尤其在低信噪比情况下倾角场估计的误差变大，但是高信噪情况下倾角场估计的结果较好，尽管在原始数据中可以观察到噪声，但估计的倾角几乎和没有噪声数据估计的倾角一致，其主要原因是正则化平滑技术约束了倾角场的平滑性，提高估计的稳定性和可靠性。在高信噪比情况下，尽管倾角场估计比较准确，但是绕射波的分离仍然受到了噪声影响，其原因为在平面波预测中我们只把可预测的平面波提取出来，而噪声也被留在了残差中，而残差正是绕射波的估计，因此在分离的绕射波中含有噪声。噪声水平越强，分离的绕射波效果越差，在低信噪比条件下，倾角场的估计也变差，绕射波分离结果几乎被淹没在噪声和假象下，无法得到较好的绕射波信息。由于噪声会降低倾角场估计的精度，同时噪声也会影响分离算法的效果，因此，在利用平面波预测法进行绕射波分离时，必须要先压制随机噪音，从而提高分离结果的精度。
![]()
图 2 不含噪声地震数据的倾角场估计和绕射波分离a. 原始含有绕射波的数据，b. 估计的倾角场，c. 分离的绕射波。Figure 2. The dip field estimation and diffraction separation of seismic data without noisea. The original data containing diffraction wave, b. estimated dip field, c. separated diffraction wave.![]()
图 3 高信噪比地震数据的倾角场估计和绕射波分离（噪声方差为0.001）a.原始含有绕射波的数据，b.估计的倾角场，c.分离的绕射波。Figure 3. The dip field estimation and diffraction separation of high SNR seismic data（noise variance is 0.001）a. The original data containing diffraction wave, b. estimated dip field, c. separated diffraction wave.![]()
图 4 低信噪比地震数据的倾角估计和绕射波分离（噪声方差为0.025）a.原始含有绕射波的数据，b.估计的倾角场，c.分离的绕射波。Figure 4. The dip estimation and diffraction separation of low SNR seismic data（noise variance is 0.025）a. The original data containing diffraction wave, b. estimated dip field, c. separated diffraction wave.倾角估计问题属于反问题范畴，存在多解性和不适定性，为了解决这一问题，需要采用正则化技术对估计的局部倾角场进行约束。在反问题正则化中光滑半径参数控制着正则化程度，下面讨论在倾角估计过程中，沿着水平方向平滑和沿垂直方向平滑半径大小对倾角预测和绕射波分离结果的影响。图5显示了不同平滑半径下倾角估计和绕射波分离的结果。随着水平或垂直方向光滑半径的不断增大，倾角估计值的范围越来越小，倾角估计的结果在水平或垂直方向上越来越平滑，随着光滑半径的进一步增大，理论上当光滑半径趋于无穷大时，倾角估计结果也会趋于与坐标轴平行。从对应的绕射波分离结果来看，随着光滑半径的增大，倾角估计值变小，绕射波分离结果显得越来越干净，信噪比明显提高。具体原因有两方面：一是大的光滑半径将部分噪音平滑掉了，减少了噪音对分离结果的干扰；二是大的光滑半径将小的绕射波忽略了，这样大尺度的绕射波就显得很干净。光滑半径太小会使倾角估计的范围大于理论值，使绕射波分离结果中出现许多全波场中没有的小尺度绕射波，使得绕射波分离结果信噪比降低，利用这些小尺度绕射波偏移成像后可能出现许多小尺度地质构造体假象，即造成偏移噪音和偏移假象[19]。
![]()
图 5 不同平滑半径对倾角估计及绕射分离结果的影响a-b是光滑半径为3的结果，c-d是光滑半径为10的结果，e-f是光滑半径为50的结果。Figure 5. Influence of different smoothing radius on dip angle estimation and diffraction separation results.a-b is the result of smooth radius of 3, c-d is the result of smooth radius of 10, e-f is the result of smooth radius of 50.因此，在应用平面波预测进行倾角估计和绕射波分离时应选择合适的光滑半径。光滑半径较小时，算法对小型地质体的绕射波识别更敏感。光滑半径较大时，小尺度绕射波没有很好地分离出来，不利于小尺度地质体的识别与成像。但较大的光滑半径对大型的绕射波聚焦顶点的识别非常清晰。在处理实际数据时，应根据具体的研究目的对平滑参数做适当的调整。例如，识别的目标是地下小尺度构造体、小型缝洞时，可以选择较小的光滑半径；若研究的目标是大型裂缝、断层、不整合面以及地层缺失时，应当选择较大的光滑半径。
3. 地震数据绕射波分离算例
为了说明平面波预测地震绕射波分离的效果，我们将其应用于一个模拟数据和两个实际数据绕射波分离中。图6a是一个模拟的共偏移距道集，可以清楚地看到绕射波和反射波互相干扰。图6b是预测出的倾角场，这里的光滑半径为5。图6c是经过分离后的结果，可以看出分离的绕射波中反射波同相轴几乎不存在，绕射波信息明显，该算例表明平面波预测滤波器可以很好地预测出反射波信息并对其进行压制。
![]()
图 6 模拟的具有绕射波的共偏移距道集a.原始数据，b.预测倾角场，c.分离的绕射波。Figure 6. The synthetic common offset gather with diffracted wavesa. original data, b. prediction of dip field, c. separated diffracted waves.然后，我们将其应用于实际海洋地震勘探的共偏移距叠前道集（图7a），可以发现该数据绕射波非常发育，从叠前共偏移距道集上可以清楚看到一些绕射波，尤其是对于3.5 s处复杂断块引起的绕射波非常清楚，我们首先估计其局部倾角场（图7b），光滑半径确定为5。在反射波发育的地方，倾角场比较平滑且数值较小，而在绕射波发育的地方倾角较大，从分离的绕射波结果（图7c）可以看出平面波预测方法可以很好地分离绕射波，这些被分离出的绕射波可以很好地进行断块构造的成像[20]及储层表征[21]。
![]()
图 7 实际海洋地震数据的共偏移距道集a.原始数据，b.估计的倾角场，c.分离的绕射波。Figure 7. The common offset gather of real marine seismic dataa. original data, b. estimated dip field, c. separated diffracted waves.最后，我们对叠后的实际地震数据进行绕射波分离（图8a），由于叠加会破坏一些绕射波的信息，因此在叠后的剖面上绕射波能量较弱，但是叠后的剖面信噪比较高，对估计倾角场比较有利。图8b为估计的局部倾角场信息，由于叠加对绕射波的压制，同时因为该数据剖面上的绕射波分布散乱，为了得到正确的倾角，此处的横向光滑半径为50，纵向光滑半径为20，图8c显示了叠后地震数据绕射波分离后的结果，可以发现绕射波很好地被分离出来。模拟和实际数据表明平面波预测方法可以有效分离绕射波。
![]()
图 8 实际地震数据叠加剖面a.原始数据，b.估计的局部倾角，c.分离的绕射波。Figure 8. The field poststack seismic sectiona. original data, b. estimated local dip field, c. separated diffraction wave.4. 结论
平面波预测滤波器在反射波同相轴光滑连续的条件下，可以有效地估计地震数据的局部倾角场，进而分离出绕射波，该方法对隐藏在强反射能量下的小型绕射波信息起到很好的分离效果，为后期绕射波单独成像提供绕射波信息。通过对噪声和平滑参数对地震绕射波分离结果的影响分析，发现噪声水平严重影响绕射波分离效果，噪音过大分离结果会不准确；平面波预测滤波器光滑半径偏小会引入噪音和假象，光滑半径偏大会导致绕射波分离不彻底，小断块地震绕射波响应难以分离。在绕射波分离算法中，振幅保真性也是一个重要问题，因此提高算法的抗噪性和保真性是今后研究的重要方向。
-
![]()
图 1 离散元数值模拟俯冲变形实验装置示意图
Figure 1. Schematic diagram of subduction deformation experimental device for discrete element numerical simulation
![]()
图 2 30°俯冲角离散元数值模拟实验变形过程解释图
Figure 2. Interpretation of the deformation process in the discrete element simulation experiment at 30° subduction angle
![]()
图 3 15°俯冲角数值模拟实验变形过程解释图
Figure 3. Interpretation of the deformation process in the discrete element simulation experiment at 15° subduction angle
![]()
图 4 5°俯冲角离散元模拟实验变形过程解释图
Figure 4. Interpretation of the deformation process in the discrete element simulation experiment at 5° subduction angle
![]()
图 5 离散元模拟实验累计位移量图
Figure 5. Cumulative displacement from the discrete element simulation experiments
![]()
图 6 3组实验中F1断层断距随俯冲量变化图
Figure 6. Fault displacement variation of F1 with subduction in three experiments
![]()
图 7 3组实验中俯冲楔高度随俯冲量变化图
Figure 7. Variation of wedge height with subduction in three experiments
![]()
![]()
![]()
图 10 实验结果与马克兰增生楔对比
a. 马克兰增生楔，b. 5°俯冲角离散元数值模拟实验结果构造图（据Grando等改, 2007）
Figure 10. Comparison of the experimental results with Makran's accretion wedge
a. Maklan accretion wedge, b. Result of 30° subduction angle simulation experiment (from Grando et al., 2007)
![]()
表 1 离散元数值模型参数
Table 1 Parameters of the discrete element numerical simulation
实验数据 弹性模量/Pa 摩擦系数 法向与切向黏结强度/Pa 速度/（m·s−1） 俯冲角/（°） 实验1 陆壳地层2×107
洋壳地层4×1080.6 3.6×105 0.5 30 实验2 陆壳地层2×107
洋壳地层4×1080.6 3.6×105 0.5 15 实验3 陆壳地层2×107
洋壳地层4×1080.6 3.6×105 0.5 5 
下载: 导出CSV
-
[1] Cundall P A. A computer model for simulating progressive large scale movements in blocky rock systems[C]//Proceedings of International Proceedings Symposium on Rock Fracture. Nancy, France, 1971: 128-132.
[2] 蔡申阳, 尹宏伟, 李长圣, 等. 基于离散元数值模拟的应变分析和裂缝预测技术[J]. 高校地质学报, 2016, 22(1):183-193. [CAI Shenyang, YIN Hongwei, LI Changsheng, et al. Technology of strain analysis and fracture prediction based on DEM numerical simulation [J]. Geological Journal of China Universities, 2016, 22(1): 183-193. [3] Iwashita K, Oda M. Micro-deformation mechanism of shear banding process based on modified distinct element method [J]. Powder Technology, 2000, 109(1-3): 192-205. doi: 10.1016/S0032-5910(99)00236-3
[4] Strayer L M, Suppe J. Out-of-plane motion of a thrust sheet during along-strike propagation of a thrust ramp: a distinct-element approach [J]. Journal of Structural Geology, 2002, 24(4): 637-650. doi: 10.1016/S0191-8141(01)00115-8
[5] Finch E, Hardy S, Gawthorpe R. Discrete element modelling of contractional fault-propagation folding above rigid basement fault blocks [J]. Journal of Structural Geology, 2003, 25(4): 515-528. doi: 10.1016/S0191-8141(02)00053-6
[6] Strayer L M, Erickson S G, Suppe J. Influence of Growth Strata on the Evolution of Fault-Related Folds—Distinct-Element Models[M]//McClay K R. Thrust Tectonics and Hydrocarbon Systems. American Association of Petroleum Geologists, 2004: 413-437.
[7] Hardy S, Finch E. Discrete-element modelling of detachment folding [J]. Basin Research, 2005, 17(4): 507-520. doi: 10.1111/j.1365-2117.2005.00280.x
[8] González G, Gerbault M, Martinod J, et al. Crack formation on top of propagating reverse faults of the Chuculay Fault System, northern Chile: insights from field data and numerical modelling [J]. Journal of Structural Geology, 2008, 30(6): 791-808. doi: 10.1016/j.jsg.2008.02.008
[9] 张洁, 尹宏伟, 徐士进. 用离散元方法讨论岩石强度对主动底辟盐构造断层分布模式的影响[J]. 南京大学学报: 自然科学, 2008, 44(6):642-652. [ZHANG Jie, YIN Hongwei, XU Shijin. Influence of rock strength on fault patterns above active salt domes: insights from 2D discrete element simulations [J]. Journal of Nanjing University: Natural Sciences, 2008, 44(6): 642-652. [10] Abe S, Van Gent H, Urai J L. DEM simulation of normal faults in cohesive materials [J]. Tectonophysics, 2011, 512(1-4): 12-21. doi: 10.1016/j.tecto.2011.09.008
[11] Egholm D L, Sandiford M, Clausen O R, et al. A new strategy for discrete element numerical models: 2. Sandbox applications [J]. Journal of Geophysical Research: Solid Earth, 2007, 112(B5): B05204.
[12] Hardy S. The Devil truly is in the detail. A cautionary note on computational determinism: Implications for structural geology numerical codes and interpretation of their results [J]. Interpretation, 2015, 3(4): SAA29-SAA35. doi: 10.1190/INT-2015-0052.1
[13] Botter C, Cardozo N, Hardy S, et al. Seismic characterisation of fault damage in 3D using mechanical and seismic modelling [J]. Marine and Petroleum Geology, 2016, 77: 973-990. doi: 10.1016/j.marpetgeo.2016.08.002
[14] De Franco R, Govers R, Wortel R. Nature of the plate contact and subduction zones diversity [J]. Earth and Planetary Science Letters, 2008, 271(1-4): 245-253. doi: 10.1016/j.jpgl.2008.04.019
[15] Lallemand S, Heuret A, Boutelier D. On the relationships between slab dip, back-arc stress, upper plate absolute motion, and crustal nature in subduction zones [J]. Geochemistry, Geophysics, Geosystems, 2005, 6(9).
[16] Billen M I, Hirth G. Rheologic controls on slab dynamics [J]. Geochemistry, Geophysics, Geosystems, 2007, 8(8): Q08012.
[17] Hardy S, McClay K, Muñoz J A. Deformation and fault activity in space and time in high-resolution numerical models of doubly vergent thrust wedges [J]. Marine and Petroleum Geology, 2009, 26(2): 232-248. doi: 10.1016/j.marpetgeo.2007.12.003
[18] Burbidge D R, Braun J. Numerical models of the evolution of accretionary wedges and fold-and-thrust belts using the distinct-element method [J]. Geophysical Journal International, 2002, 148(3): 542-561. doi: 10.1046/j.1365-246x.2002.01579.x
[19] Naylor M, Sinclair H D, Willett S, et al. A discrete element model for orogenesis and accretionary wedge growth [J]. Journal of Geophysical Research: Solid Earth, 2005, 110(B12): B12403. doi: 10.1029/2003JB002940
[20] Egholm D L. A new strategy for discrete element numerical models: 1. Theory [J]. Journal of Geophysical Research: Solid Earth, 2007, 112(B5): B05203.
[21] Liu Z N, Koyi H A. The impact of a weak horizon on kinematics and internal deformation of a failure mass using discrete element method [J]. Tectonophysics, 2013, 586: 95-111. doi: 10.1016/j.tecto.2012.11.009
[22] Okamura Y, Nishizawa A, Oikawa M, et al. Differential subsidence of the forearc wedge of the Ryukyu (Nansei-Shoto) Arc caused by subduction of ridges on the Philippine Sea Plate [J]. Tectonophysics, 2017, 717: 399-412. doi: 10.1016/j.tecto.2017.08.025
[23] Nishizawa A, Kaneda K, Oikawa M, et al. Variations in seismic velocity distribution along the Ryukyu (Nansei-Shoto) Trench subduction zone at the northwestern end of the Philippine Sea plate [J]. Earth, Planets and Space, 2017, 69(1): 86. doi: 10.1186/s40623-017-0674-7
[24] Lallemand S, Heuret A, Boutelier D. On the relationships between slab dip, back‐arc stress, upper plate absolute motion, and crustal nature in subduction zones [J]. Geochemistry, Geophysics, Geosystems, 2015, 6(9): Q09006.
[25] Skinner S M, Clayton R W. An evaluation of proposed mechanisms of slab flattening in central Mexico [J]. Pure and Applied Geophysics, 2011, 168(8-9): 1461-1474. doi: 10.1007/s00024-010-0200-3
[26] 皇甫鹏鹏, 王岳军, 范蔚茗, 等. 大洋平板俯冲的数值模拟再现: 洋–陆汇聚速率影响[J]. 大地构造与成矿学, 2016, 40(3):429-445. [HUANGFU Pengpeng, WANG Yuejun, FAN Weiming, et al. Numerical modeling of flat subduction: constraints from the ocean-continent convergence velocity [J]. Geotectonica et Metallogenia, 2016, 40(3): 429-445. [27] Yang C Q, Yang Y Q, Li G, et al. The Mesozoic basin-mountain coupling process of the southern East China Sea shelf basin and its adjacent land area [J]. Acta Geologica Sinica - English Edition, 2016, 90(3): 1051-1052. doi: 10.1111/1755-6724.12748
[28] 侯方辉, 张训华, 李刚, 等. 从被动陆缘到主动陆缘——东海陆架盆地中生代构造体制转换的盆地记录[J]. 石油地球物理勘探, 2015, 50(5):980-990. [HOU Fanghui, ZHANG Xunhua, LI Gang, et al. From passive continental margin to active continental margin: basin recordings of Mesozoic tectonic regime transition of the East China Sea Shelf Basin [J]. Oil Geophysical Prospecting, 2015, 50(5): 980-990. [29] Shang L N, Zhang X H, Jia Y G, et al. Late Cenozoic evolution of the East China continental margin: Insights from seismic, gravity, and magnetic analyses [J]. Tectonophysics, 2017, 698: 1-15. doi: 10.1016/j.tecto.2017.01.003
[30] 杨长清, 杨传胜, 李刚, 等. 东海陆架盆地南部中生代构造演化与原型盆地性质[J]. 海洋地质与第四纪地质, 2012, 32(3):105-111. [YANG Changqing, YANG Chuansheng, LI Gang, et al. Mesozoic tectonic evolution and prototype basin characters in the southern East China Sea Shelf Basin [J]. Marine Geology & Quaternary Geology, 2012, 32(3): 105-111. [31] 崔幸, 王亮亮, 罗洪明, 等. 东海陆架盆地南部中生代盆地性质与演化: 砂箱物理模拟检验[J]. 海洋地质与第四纪地质, 2017, 37(4):181-192. [CUI Xing, WANG Liangliang, LUO Hongming, et al. Sandbox modeling test for Mesozoic basins in southern East China Sea Shelf Basin [J]. Marine Geology & Quaternary Geology, 2017, 37(4): 181-192. [32] Cui X, Dai L M, Li S D, et al. Control of strike‐slip and pull‐apart processes to tectonic transition of the southern East China Sea Shelf Basin [J]. Geological Journal, 2019, 54(2): 850-561. doi: 10.1002/gj.3363
[33] Liu Z, Dai L M, Li S Z, et al. Mesozoic magmatic activity and tectonic evolution in the southern East China Sea Continental Shelf Basin: thermo‐mechanical modelling [J]. Geological Journal, 2018, 53(S1): 240-251.
[34] Hu M Y, Li S Z, Dai L M, et al. Dynamic mechanism of tectonic inversion and implications for oil–gas accumulation in the Xihu Sag, East China Sea Shelf basin: insights from numerical modelling [J]. Geological Journal, 2018, 53(S1): 225-239.
[35] Abedi M, Bahroudi A. A geophysical potential field study to image the Makran subduction zone in SE of Iran [J]. Tectonophysics, 2016, 688: 119-134. doi: 10.1016/j.tecto.2016.09.025
[36] Kopp C, Fruehn J, Flueh E R, et al. Structure of the makran subduction zone from wide-angle and reflection seismic data [J]. Tectonophysics, 2000, 329(1-4): 171-191. doi: 10.1016/S0040-1951(00)00195-5
[37] Grando G, McClay K. Morphotectonics domains and structural styles in the Makran accretionary prism, offshore Iran [J]. Sedimentary Geology, 2007, 196(1-4): 157-179. doi: 10.1016/j.sedgeo.2006.05.030
-
期刊类型引用(2)
1. 许胜利，张雨，孙玉红，许靖康，滑彩虹，李灿灿. 射线原理下中-深层裂缝地震可预测性分析及应用. 石油工业技术监督. 2024(11): 1-5 . 
百度学术
2. 栾锡武，杨佳佳. 地震绕射波波场分离与成像方法综述. 石油物探. 2022(05): 761-770 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 3166
- HTML全文浏览量: 594
- PDF下载量: 37
- 被引次数: 2
