第一作者简介:彭靖淞(1984-),男,重庆人,硕士,中海石油(中国)有限公司天津分公司高级工程师,主要从事构造与石油地质综合研究。地址:天津市滨海新区海川路2121号,中海石油(中国)有限公司天津分公司,邮政编码:300459。E-mail: pengjs@cnooc.com.cn
以三维地震和钻井资料为基础,通过地震属性分析、构造几何学分析、断层活动性分析、构造演化史分析和油气运移模拟等方法,对张家口—蓬莱断裂带渤海沙垒田凸起东北段的构造样式、断裂活动性、走滑位移量、中央构造带的形成机制及其油气运聚进行分析。研究区是早期(古新世—始新世)伸展和晚期(渐新世、上新世—第四纪)走滑-拉分的叠合构造区;张家口—蓬莱断裂渤海沙垒田凸起东北段左旋走滑活动有两个阶段,即渐新世和上新世—第四纪,渤海海域活动较陆地段更早;新生代以来左旋走滑位移量共4 km,其中渐新世活动较弱(1 km),上新世较强(3 km);大规模的走滑活动导致了断裂带围区断裂活动的加剧和盆岭再造,促进了中央构造带的形成和北东向构造带的协同演化;张家口—蓬莱断裂渤海沙垒田凸起东北段及其共轭走滑断裂带剧烈的新构造运动控制了油气从深层往浅层运移的充注点/段,从而控制了其浅层油气运移及成藏。
The structural style, fault activity, strike-slip displacement, and the formation mechanism and hydrocarbon migration and accumulation in the center tectonic zone in the northeast Shaleitian Bulge of Zhangjiakou-Penglai Fault Zone were studied by seismic attribute analysis, structural geometric analysis, fault activity analysis, structure evolution history and simulation of hydrocarbon migration, based on 3-D seismic and drilling data. The main results are as follows: (1) The study area is a superimposed tectonic zone, which experienced early (Paleocene and Eocene) extension and late (Oligocene and Pliocene-Quaternary) strike-slip and pull-apart. (2) The sinistral strike slip of the northeast Shaleitian Bulge of Zhangjiakou-Penglai Fault Zone went through two periods, Oligocene and Pliocene-Quaternary, and the Bohai section was active earlier than the inland section. (3) The sinistral strike slip displacement of Zhangjiakou-Penglai Fault is 4 km since Cenozoic, including 1 km in the Oligocene, and 3 km in the Pliocene-Quaternary. (4) The strike-slip movements have resulted in the increase of fault activity and basin-mountain restructure in the fault zone, also contributed to the formation of the central tectonic belt and the conjugate evolution in north-east structural belt. (5) The conjugate strike slip of the Zhangjiakou-Penglai Fault Zone dominated the migration and accumulation of hydrocarbon in shallow formations by controlling the injection points and segments of hydrocarbon from the deep layers to shallow layers.
走滑断裂对造山和成盆的作用一直是大地构造学研究的热点[1, 2, 3, 4, 5, 6]。渤海湾盆地被中国东部的两条巨型走滑断裂:郯城— 庐江断裂(下文简称郯庐断裂)和张家口— 蓬莱断裂带(下文简称张蓬断裂)所切割, 走滑样式丰富, 盆岭关系复杂[3, 4, 5, 6, 7, 8, 9]。近年来, 渤海海域油气勘探积累了大量钻井和三维地震资料, 为走滑断裂的大地构造学研究提供了良好的资料基础。
张蓬断裂西起张家口附近, 向东南经北京、天津、渤海, 东南段沿山东半岛北缘和威海及沿海海域延伸, 是一条走向近北西向的活动断裂带, 总长度超过700 km, 与北北东向的郯庐断裂带构成渤海海域新生代两个主要的活动构造带[6, 7, 8, 9]。与郯庐断裂相似, 张蓬断裂渤海段亦是渤海地区主要的油气聚集区带[10, 11, 12], 累计发现石油地质储量超过12× 108 t, 包括蓬莱9-1、渤中19-4、渤中8-4、曹妃甸12-6和曹妃甸6-4等断块油田。
张蓬断裂具有左旋走滑特征, 现今仍非常活跃, 地震活动频繁, 其主干断裂被多条北东— 南西向断裂所切割, 分段差异性明显[6, 7, 8, 9]。前人研究成果主要集中在张蓬断裂的陆上地区(张北— 南口段、南口— 宁河段和蓬莱— 烟台段)[6, 7, 8, 9, 10, 11, 12], 对渤海海域段的研究相对缺乏, 其走滑构造特征和走滑尺度不明。
渤海海域左旋张蓬断裂是右旋郯庐断裂的共轭走滑断裂[10, 11, 12], 共轭断裂具有同时、同步、协调活动的特征[10, 13, 14]。索艳慧等学者认为张蓬断裂陆地段第四纪(距今约2.6 Ma)开始左旋[6, 7, 8, 9], 其渤海段的发育时间是否与陆地段相同?如果相同, 这比郯庐断裂渐新世[15, 16, 17]右旋走滑开始时间(距今约30 Ma)要滞后近27 Ma, 有悖于“ 共轭” 断层同生、共生的认识, 因此张蓬断裂渤海段的形成时间还需进一步落实。
受太平洋板块俯冲和喜马拉雅造山运动远程效应的共同影响[3, 5, 6, 18, 19, 20, 21], 渤海海域新生代以来具有伸展和走滑的双动力源, 表现出早期(古新世— 始新世)伸展、晚期(渐新世至今)走滑-拉分的叠合构造样
式[18, 19, 20, 21, 22, 23, 24]。伸展构造会形成箕状断陷[18, 19, 20, 21, 22, 23, 24], 而走滑-拉分构造往往形成盆岭构造[1, 2, 3, 4, 5, 6]。在早期先存的伸展构造基础上发生走滑, 渤海地区的盆地和凸起的格局必然发生变化, 笔者称之为“ 盆岭再造” 。前人认为走滑弯曲、伸缩变形及其差异沉降/隆升是走滑带造山、造盆的主要机制[1, 2, 3, 4, 5, 6, 19, 20, 21]。笔者研究发现:除了以上机制, 对于平移断距较大的走滑断层, 两盘长距离错动形成的“ 海豚效应” [25, 26, 27]对盆岭再造的影响也值得重视, 这将在下文中详细论述。
晚期成藏是渤海地区油气成藏的重要特征, 新构造运动对活动断裂带的浅层成藏至关重要[10, 11, 12, 28, 29, 30, 31, 32, 33, 34, 35]。张蓬断裂是渤海湾盆地现今最为活跃的断裂带和油气富集区[5, 6], 前人对其复杂断块油气藏的石油地质条件进行过一定的研究, 认为张蓬断裂控制了圈闭的形成, 晚期活跃的断层把油气从深层输导到浅层在新近系成藏[10, 11, 12]。实际钻井和三维地震资料揭示, 张蓬断裂带深浅断裂叠合程度不高, 多期次构造发育, 盆地叠合/再造特征明显, 各个断块圈闭之间成藏差异较大, 张蓬断裂带对油气成藏的控制作用有待进一步研究。
本文研究以张蓬断裂渤海段的沙垒田凸起东北段为例, 以三维地震资料和钻井资料为基础, 通过构造解释、地震属性分析、构造几何分析、断层活动性分析、构造演化史分析和油气的运移模拟等方法, 分析张蓬断裂带的构造样式、走滑位移量、中央构造带的形成机制、断裂走向上活动性差异及其对油气运移的影响, 以期理清张蓬断裂的构造特征、活动期次和走滑规模, 探索张蓬断裂带盆岭再造的机制, 总结走滑带构造运动及其断裂活动对油气成藏的控制作用。
研究区位于渤海湾盆地渤中坳陷沙垒田凸起东北部, 张蓬断裂位于研究区西侧, 北西— 南东向贯通全区。该区域构造复杂, 发育依附于张蓬断裂的中央构造带、与张蓬断裂斜交的渤中8-4和渤中8-3构造带(见图1)。地层可划分为深浅两部分, 深层:古近系由浅到深依次发育东营组、沙河街组和孔店组, 主要为多期扇三角洲— 辫状河三角洲砂体和湖相泥岩的旋回沉积; 浅层:新近系由浅到深依次发育明化镇组和馆陶组, 主要为河流相和浅水三角洲相的砂砾岩和砂(泥)岩沉积[35, 36]。
2.1.1 典型剖面特征
张蓬断裂渤海沙垒田凸起东北段(下文简称张蓬断裂沙东北段)具有典型的走滑断裂带特征:剖面上, 其主走滑断裂及其派生断裂组合成负花状形态, 表现出伸展-走滑的应力背景[19, 22, 25, 26, 27]; 主走滑断面近乎直立, 断裂倾向在研究区南段朝东、北段朝西, 表现出走滑断层典型的丝带效应[25, 26, 27](见图2)。同时, 研究区的东侧为正断层、铲式正断层和“ Y” 形断层组合等伸展构造控制的箕状断陷(见图2)。可见研究区是伸展构造和走滑-拉分构造的叠合发育区, 体现了渤海地区的双动力源特征[18, 19, 20, 21, 22, 23, 24]。
2.1.2 走滑派生断层与走滑方向
张蓬断裂沙东北段在东营组以下的深层为一条北西— 南东向贯穿全区的断层(见图3a、图3b), 但在东营组及其上覆的浅层, 张蓬断裂不再是一条断层, 而是逐渐变成一系列具有走滑-伸展特征的派生断层, 这些派生断层剖面上表现为花状构造的花瓣, 平面上呈北东东— 南西西雁列展布, 形成了一个宽度为10~20 km的断裂带(见图3c、图3d)。同时, 张蓬断裂沙东北段主断裂与派生断裂之间的夹角指示其新生代以来为左旋走滑构造[17, 37]。
2.2.1 区域断裂演化
长期活动断裂的活动期次反映了构造的演化期次[38, 39]。通过对研究区长期活动断裂的研究表明, 研究区主要发育3期构造运动[17, 18, 19, 20, 21, 22, 28, 35](见图1、图4):①裂陷早幕, 古新统孔店组— 始新统沙河街组沉积时期, 该阶段主要发育了一些控凹的正断层, 断层的数量较少, 但单条断层的活动速率较大, 断层走向主要为北北西— 南南东向; ②裂陷晚幕, 渐新统东营组沉积时期, 该时期断裂的活动速率最高, 走滑-伸展断层的数量有一定程度的增加, 断层的主要走向变成北东东— 南西西向; ③新构造运动, 上新统明化镇组— 第四系平原组沉积时期, 经过中新世的整体热沉降后, 断裂又重新活化, 张蓬断裂带走滑-伸展特征派生断层的数量开始大量增加, 其派生断层数量占到了全区断层的72%, 且断层活动速率较大, 断层的走向与裂陷晚幕基本一致, 主要为北东东— 南西西向(见图5)。
2.2.2 走滑构造的活动时间
在了解区域断裂活动史的基础上, 对张蓬断裂沙东北段的走滑活动史进行了分析。由于几乎所有的走滑断裂活动, 都会形成一系列张扭性的派生断裂, 所以派生断裂活动史可以间接反映主走滑断裂的活动史[17]。通过对张蓬断裂沙东北段的走滑派生断层的分析表明, 其走滑活动开始于渐新世裂陷晚幕, 主要有两期走滑活动, 即渐新世裂陷晚幕和上新世— 第四纪新构造运动期, 后者的断裂活动强度几乎是前者的3倍, 形成了大量的派生断层(见图6)。
张蓬断裂的走滑活动发生在裂陷晚幕和新构造运动期, 但“ 缺席” 了伸展构造为主的裂陷早幕。笔者认为, 研究区是早期(裂陷早幕)伸展和晚期(裂陷晚幕, 新构造运动)走滑-拉分的叠合构造发育区。
对比张蓬断裂陆地段的活动时间[6, 7, 8, 9], 渤海沙东北段活动较陆地段活动更早, 张蓬断裂可能是由渤海海域向陆地发育。
2.2.3 与郯庐断裂带的演化关系
结合郯庐断裂带[15, 16, 17], 笔者认为张蓬断裂的左旋走滑与郯庐断裂新生代的右旋走滑都没有发生在渤海海域的裂陷早幕, 但几乎同时发育于渐新世裂陷晚幕和新构造运动时期, 它们具有共轭走滑断裂同时、协同发育的特征, 其共轭走滑活动促进了渤海海域在古— 始新世伸展背景下, 渐新世走滑-拉分再裂陷和上新世构造继承性活化(见表1)。
张蓬断裂现今仍处于左旋走滑的构造活动中[6, 7, 8, 9], 但对其走滑位移量尚缺乏研究。走滑断距的估算一直是走滑构造研究的难点[40, 41, 42, 43]。笔者综合采用了基底参考点走滑断距分析法、基底地貌相似性走滑断距分析法、沉积体系对比走滑断距分析和走滑活动断距累加原理, 对张蓬断裂沙东北段渐新世、上新世— 第四纪两个主要走滑时期的走滑量和新生代总走滑量进行了估算。
2.3.1 新生代总走滑量
新生界基底在张蓬断裂沙东北段左旋走滑过程中发生错动, 其错动的距离是张蓬断裂沙东北段新生代的总断距。笔者采用两种方法对其走滑量进行了计算。
①基底参考点走滑断距分析法。走滑活动前, 东西两盘的基底地貌应该是相似度最高的[25]。根据此原理, 在走滑断裂的两盘各选取两组以上地貌特征比较相似的参考点, 量取各组断点在走滑断裂上垂直投影位置之间的距离, 并求取平均值。通过这个方法, 在张蓬断裂沙东北段主走滑断层两盘基底上选取了3组参考点(A和A° 、B和B° 、C和C° ), 量取3组断距, 计算平均值, 得到新生代平均断距为4 km。同时, 通过反向走滑(右旋)4 km, 得到张蓬断裂沙东北段活动之前的古地貌, 古地貌过渡平顺, 也证明了走滑量估计合理(见图7)。
②基底地貌相似性走滑断距分析法。同样根据走滑前两盘地貌相似原理[25], 通过基底地貌相似性断距分析法计算断距:选取走滑断裂附近两盘沿走滑断裂走向的4 km古地貌线为样本, 对样本线进行200等分, 并将两盘匀速按着每步200 m反左旋, 即右旋错动恢复古地貌(见图8), 依次计算每一次错动对应两盘海拔差平方和:
\[{{P}_{k}}=\sum\limits_{i=1}
{200}{{{\left( {{H}_{i}}-{{H}_{i+k}} \right)}
{2}}} \ \ (1)\]
当Pk为最小值时, 两盘的地貌最为接近, 则走滑断距为:
$L=0.2{{k}_{\text{min}}}$ (2)
经过计算机编程计算发现, 当右旋步数为20时Pk最小, 进而计算出新生代走滑断距为4 km。这与参考点分析法估算的断距相等, 证明了新生代走滑断距应该为4 km。通过反向走滑(右旋)4 km, 得到张蓬断裂沙东北段活动之前的古地貌, 古地貌过渡平顺, 证明新生代走滑量估计合理。
2.3.2 上新世— 第四纪走滑量
如果走滑断层切割沉积体, 那么沉积体会发生错动, 因此沉积体记录了沉积之后的走滑位移量[43]。据此原理, 笔者采用沉积体系对比走滑断距分析法估算了馆陶组上新世— 第四纪的走滑量:馆陶组砂砾岩体的地震相为中— 弱振幅弱连续性反射, 根据地震属性刻画了该砂砾岩体的范围, 发现张蓬断裂沙东北段主走滑断裂东西两侧的砂砾岩体错动明显(见图9)。通过东西两盘砂砾岩体及其沉积体系的参考点分析与测量, 馆陶组沉积之后上新世— 第四纪的走滑断距是3 km(见图9)。
2.3.3 渐新世走滑位移量
根据走滑活动断距的累加原理, 新生代走滑位移量等于渐新世走滑位移量与上新世以来走滑位移量之和, 因此, 渐新世走滑位移量等于新生代走滑位移量(4 km)减去上新世以来走滑位移量(3 km), 由此可得渐新世走滑位移量为1 km。
通过以上研究表明:新生代张蓬断裂沙东北段总的走滑位移量为4 km, 其中渐新世走滑位移量为1 km, 上新世— 第四纪走滑错动3 km, 上新世— 第四纪走滑活动较渐新世更为强烈。这与走滑派生断层的断层活动变化速率趋势一致, 也符合张蓬断裂第四纪以来活动较强、地震活动频繁的构造背景[6, 7, 8, 9]。
与其共轭走滑断裂— — 郯庐断裂新生代的右旋走滑特征相比[15, 16, 17](见表1), 张蓬断裂沙东北段总的左旋走滑位移量为4 km, 远小于郯庐断裂10~20 km的右旋走滑位移量[3, 5, 19, 40], 这说明张蓬断裂活动尺度弱于郯庐断裂。郯庐断裂渐新世的活动强于上新世— 第四纪, 张蓬断裂上新世— 第四纪活动强于渐新世, 即前者早强晚弱, 而后者早弱晚强。渐新世以来, 张蓬断裂与郯庐断裂活动强度互为消长。笔者推测, 从古近纪到新近纪, 喜马拉雅运动的远端效应和太平洋板块俯冲方向变化的共同作用使华北板块由向北东向挤出变为向东挤出[3, 5, 6, 19], 由此导致共轭断裂之间活动强度的相互转化, 同时也造成现今张蓬断裂带的地震活动较郯庐断裂更为活跃。
综合研究表明, 研究区是早期(裂陷早幕)伸展、晚期(裂陷晚幕— 新构造运动)走滑-拉分的叠合区域。晚期的走滑-拉分运动是在早期伸展构造背景下的盆地格局的重要变革, 发生盆岭再造, 形成了中央构造带。
张蓬断裂沙东北段的左旋活动控制了中央构造带的演化(见图10):①走滑错动形成的“ 海豚效应” 奠定了中央构造带的基础。张蓬断裂左旋斜切了沙垒田凸起高部位的“ 一角” , 这“ 一角” 相对于西盘沙垒田凸起的主体向北错动了4 km, 平移到了凹陷区, 发生“ 海豚效应” [25, 26, 27], 形成了凹陷区相对较高的中央构造带的雏形(见图7、图10)。可见左旋走滑的“ 海豚效应” 是形成中央构造带的基础。②走滑拉分、差异升降促使中央构造带进一步强化。在张蓬断裂沙东北段主走滑断裂斜切凸起“ 一角” 的同时, 由于研究区张蓬断裂的断面为“ S” 形弯曲, 因此在凸起边缘产生了一个局部的“ 左旋左阶” 的走滑拉分区[44, 45](见图10), 地层下陷, 盆地边缘向西再造, 促使中央构造带西侧剧烈沉降, 也导致了中央构造带的相对隆升。
由此可见, 除了差异沉降, 走滑错动的“ 海豚效应” 也是断裂带内伸展、走滑-拉分叠合和盆岭再造的重要原因。
通过对北西向中央构造带与北东向渤中8-4和渤中8-3构造带演化剖面的对比分析(见图11)发现, 渤中8-4、渤中8-3构造带与中央构造带早期独立演化, 晚期协同演化。①裂陷早幕, 古新世— 始新世时期, 发生伸展背景下的早期裂陷, 形成渤中8-4和渤中8-3构造带, 此时中央构造带还是沙垒田凸起的一部分; ②裂陷晚幕, 渐新世时期, 伴随着张蓬断裂左旋走滑的剧烈活动, 在早期伸展构造背景下叠合发育走滑构造, 中央构造带走滑错动形成, 并与渤中8-4和渤中8-3构造带凹中隆协同演化; ③新构造运动, 上新世— 第四纪时期, 伴随着张蓬断裂走滑-拉分作用的加强, 正断层大量发育, 研究区构造进一步破碎化并最终定形, 上新世是张蓬断裂带圈闭群形成的主要时期, 这为油气的“ 晚期成藏” [10, 11, 12, 28, 29, 30, 31, 32, 33, 34, 35]奠定了基础。
由此可见, 研究区的盆岭再造是早期伸展, 晚期走滑-拉分这个盆地格局剧烈改变过程中构造叠加的结果, 它既包含了新构造带— — 中央构造带的形成, 也影响了先存构造带— — 渤中8-4、渤中8-3构造带的继承性发育和定形。
综合研究表明, 张蓬断裂沙东北段渐新世以来的走滑运动及其断裂活动, 导致了围区的盆岭再造, 这一过程对研究区的晚期成藏[10, 11, 12]有重要影响。
①盆岭再造形成低势区, 从而控制了深层油气运移。张蓬断裂沙东北段的走滑活动发生盆岭再造, 导致了北西向中央构造带在渐新世的形成和上新世流体低势的强化, 同时控制了北东向的渤中8-4和渤中8-3这两个先存构造带在渐新世— 上新世以来的继承性发育(见图11), 由此形成了研究区主要的几个流体低势区, 从而控制了深层油气运移。实际钻探资料揭示:研究区整体上为深、浅复合成藏, 浅层富集区与深层聚集区叠合较好, 说明油气深层聚集区应该为浅层富集区的油气中转站[13](见图5、图12), 可见张蓬断裂沙东北段的剧烈活动导致的盆岭再造控制了研究区油气运移主要方向。
②重塑了纵向上的油气运移通道, 从而控制了“ 浅层油气运移” 。研究区主要的成藏时间为上新世以来[10, 11, 12, 28, 29, 30, 31, 32, 33, 34, 35], 这恰好是张蓬断裂左旋最为强烈的新构造运动时期, 形成了大量的走滑-伸展断层。但并不是所有的断层都是油气运移通道, 而且断层在走向上其活动速率会发生变化, 达到一定活动速率的断裂点/段才是油气充注点/段[28, 46, 47, 48]。这些油源断裂充注点/段是浅层油气横向运移的“ 起始点” , 对于浅层的油气运移非常关键(见图5、图13)。
对研究区及其周围油田的断层活动速率和浅层储量占比的关系统计表明, 当斜坡带断层活动速率低于20 m/Ma, 凹陷区断层活动速率低于30 m/Ma时, 油气主要在深层古近系成藏; 当斜坡带断层活动速率高于20 m/Ma, 凹陷区断层活动速率低于30 m/Ma时, 油气主要在浅层成藏。所以断层活动速率20 m/Ma和30 m/Ma分别是斜坡带和凹陷区在浅层形成油气有效充注的门限。
根据这两个门限, 对研究区的油源断裂的油气充注点/段进行了识别(见图5), 分析表明:①张蓬断裂及其派生断裂控制了研究区85%的充注点/段分布; ②浅层断裂虽然有137条, 但只有31条断层是油源断层, 且只有17%的断面成为了充注点/段, 研究区断层虽多, 但能作为油气运移通道的却很少。
通过这些充注点/段的约束进行了油气运移模拟(见图13), 预测了浅层的油气运移路径与成藏部位。模拟结果与现今浅层油气藏的吻合度达到85%(见图5、图13)。这说明, 一方面本文分析方法、模拟结果可靠, 能很好地预测复杂断块区的浅层油气运移路径, 具有很好的实用性; 另一方面证明了张蓬断裂带沙东北段控制浅层油气运移的起始点、充注点/段, 从而控制了研究区浅层的油气成藏。
研究区是早期(古— 始新世)伸展和晚期(渐新世、上新世— 第四纪)走滑-拉分的叠合构造发育区, 张蓬断裂渤海沙东北段左旋走滑的主要时期为渐新世和上新世— 第四纪。基底地貌和沉积体系对比表明:张蓬断裂沙东北段新生代以来左旋走滑位移量共计4 km, 其中渐新世活动较弱为1 km, 上新世— 第四纪较强为3 km。上新世— 第四纪的新构造运动是张蓬断裂左旋活动最为剧烈的时期。
张蓬断裂与郯庐断裂活动强度互为消长:张蓬断裂沙东北段上新世— 第四纪活动强于渐新世, 郯庐断裂渐新世活动强于上新世— 第四纪, 前者早弱晚强, 而后者早强晚弱。研究区的盆岭再造是早期伸展、晚期走滑-拉分盆地格局剧烈改变过程中构造叠加的结果, 它既包含了新构造带— — 中央构造带的形成, 也影
响了先存构造带— — 渤中8-4、渤中8-3构造带的继承性发育和定形。走滑带除了差异沉降, 走滑错动的海豚效应也是断裂带内盆岭再造的重要机制。
张蓬断裂沙东北段的走滑活动及其新构造运动导致了盆岭再造和低势区的形成, 从而控制了深层油气运移, 其大量发育的派生断裂重塑了纵向上的油气运移通道, 从而控制了浅层油气运移。这一认识及其分析方法对于复杂断裂带油气的深浅复合勘探有重要借鉴意义。
致谢:感谢中海油天津分公司徐长贵、牛成民、王德英、杜晓峰和张新涛等专家在张蓬断裂带构造演化研究工作中给予的支持和帮助。
符号注释:
H— — 海拔, m; i— — 样本点序号, i=1, 2, …, 200, 无因次; k— — 断层错动次数, 无因次; kmin— — Pk为最小值时对应的右旋步数, 无因次; L— — 走滑断距, km; Pk— — 第k次错动对应断层两盘海拔差平方和, m2。
The authors have declared that no competing interests exist.