-
“反转构造”概念的首次提出是在20世纪80年代,被定义为“盆地逆转成构造隆起”[1, 2]。随后几十年中,随着对反转构造研究的不断加强和深入以及地震技术的发展,这一概念逐渐得到重视并不断完善发展:是指同一地质体在不同地质历史时期,由于应力场的变化,造成挤压构造与拉伸构造垂向上叠加的一种复合构造形式[1, 3]。反转作用过程伴随盆地的发育演化过程,研究反转构造有助于深刻认识盆地的发育演化过程;同时,构造反转强度和反转期与油气生成和运聚的时间匹配关系等对油气聚集成藏有重要影响[1],更使其在油气勘探构造分析中愈来愈为人们所关注[3-10]。
盆地构造反转的性质可以通过大量的多阶段反转断裂体现出来,尤其是先存的伸展正断层和伴生的背斜或断褶等构造行迹。为了探索盆地构造反转特征及其最主要影响因素等问题,前人也已展开了很多数值模拟和物理模拟实验[11-20]:Buiter等[13]通过数值模拟实验探讨铲形断层对盆地构造反转的影响,得出先存软弱断裂带是盆地构造反转大量发生的必要条件;Ventisette等[15]实验室模拟了逆断层与早期先存正断层的耦合性的影响,认为挤压构造的形成强烈受控于先存正断层及其倾角大小;Marques等[16]通过沙箱实验探讨了盆地收缩量与收缩方向对盆地反转作用的影响;Gomes等[18]通过物理模拟实验,探讨铲形断层在盆地构造反转中的作用,得到半地堑盆地边界铲式断层的形态对于构造反转中正断层的再活动有着强烈的影响。上述数值模拟和物理模拟的实验,探讨盆地反转构造样式的一些影响因素,但还存在一些问题。如物理模拟实验中缺乏考虑先存控盆断层的形态变化;数值模拟实验多为二维平面模拟,缺乏对盆地构造反转过程的重现以及在此基础上对其动力学机制的分析等。
位于中国东部陆缘的东海陆架盆地,自晚中生代以来经历了多期构造运动的改造,形成了复杂多样的反转构造,其中以西湖凹陷最为典型。对东海陆架盆地西湖凹陷反转构造体系特征及其动力学机制的研究,有助于进一步了解盆地的形成发育历史及构造分布规律,同时,盆地反转构造与油气聚集成藏之间关系的研究,对提高海域盆地油气的勘探效益具有重要的指导意义。西湖凹陷自晚中生代以来长期处于中国东部陆缘伸展构造带内,区域应力场存在数次由伸展到挤压的转换,使得凹陷内普遍发育显著的反转构造。自20世纪90年代东海海域第一个油气田——西湖凹陷平湖油气田投入开发以来,众多学者和油气勘探工作者对西湖凹陷区域反转构造特征、沉积特征以及油气聚集规律都进行了进一步的探索和研究工作。根据大量的二维地震反射剖面以及钻孔资料等,系统定性地探讨了反转构造的分布与样式[21-23],断裂系统的几何学特征及其成因机制[19, 24-27],以及反转构造与油气聚集之间关系等[28-30]。虽然前人对西湖凹陷的研究已相当充分,但还存在一些问题,如对于西湖凹陷反转构造特征的精细刻画不足、反转过程和成因机制不明确等。本文的研究正是针对其构造反转演化过程及其成因机制不明确,根据前人对于西湖凹陷反转成因机制的不同观点,利用数值模拟的方法,建立了三种边界条件的模型,对西湖凹陷构造反转的过程及其成因机制进行系统的模拟,以此重现盆地构造演化过程,从而探讨了盆地演化的迁移规律及其动力学机制。本文的研究模型考虑了盆地岩石圈横向上的不均匀性,断层的设置也是依据现今详细地震数据构建,并且充分考虑了盆地三维空间上的不连续性。
-
东海陆架盆地位于西太平洋活动大陆边缘,是印度-澳大利亚板块和太平洋板块与欧亚板块汇聚的地带,占据了全球板块汇聚的独特位置,同时受到了印度板块的挤入、太平洋板块的后退式俯冲、台湾造山带的楔入的联合作用[31-33],使得该区域具有一个复杂的构造背景。因此,东海陆架盆地记录了因太平洋板块和菲律宾海板块俯冲速率和方向的变化、欧亚板块和印度板块碰撞形成的盆地构造反转的大量信息[19, 31-38]。盆地地理上位于闽浙隆起区以东,钓鱼岛隆褶带以西(图 1),基底为华南陆块在东海陆架的延伸,是晚白垩世以来发育起来的中、新生代盆地,盆地自西向东可划分为3个主要的构造单元,分别为西部坳陷带、中央隆起带以及东部坳陷带[38],其东部坳陷带具有沉降幅度大,新生代沉积巨厚的特点,是东海陆架盆地的主要油气聚集区[39]。
西湖凹陷位于东海陆架盆地东部坳陷带,东侧以钓鱼岛隆褶带为界,西连海礁隆起,南北分别与钓北凹陷和福江凹陷相接(图 1)[22],总体呈NNE—SSW向长条状展布,长约440km,宽约110km,总面积约5.9万km2,是以新生代碎屑沉积为主的沉积凹陷,沉积物最大厚度15 000m,具有自西向东逐渐变厚的趋势。西湖凹陷自西向东可划分为西部斜坡带、中央反转构造带及东部断阶带3个一级构造单元,以及多个二级构造单元(图 2)。其构造演化主要受太平洋板块与菲律宾海板块俯冲方向和速率的影响,大致可以分为3个阶段[26, 29](表 1):盆地形成—始新世裂陷期、渐新世—中新世拗陷反转期、中新世至今沉降期。
表 1 西湖凹陷地层柱及构造演化
Table 1. Stratigraphic column and tectonic evolution of the Xihu Sag
晚白垩世到始新世期间,由于太平洋板块俯冲速率由高到低、角度由低到高的变化,导致欧亚东南陆缘处于右旋张扭应力场下,自晚白垩世开始大规模伸展裂陷,东海陆架盆地西部开始张裂[31, 32, 34, 38]。始新世末(约32 Ma),西湖凹陷经历了第一期构造反转,即玉泉运动造成的构造反转,西部斜坡带隆升并遭受严重剥蚀,东部坳陷带局部地区则发生了并不强烈的构造反转,形成了T30构造不整合界面(表 1)。
渐新世至中新世,太平洋板块俯冲方向由斜向俯冲转变为垂向俯冲,俯冲速率逐渐上升[34, 35, 41],钓鱼岛隆褶带发生强烈隆升,西湖凹陷整体进入拗陷-反转阶段。菲律宾岛弧发生向西仰冲,使得东海陆架盆地开始处于挤压或左旋压扭的构造背景,渐新世末(约23 Ma),西湖凹陷发生第二期构造反转,即花港运动造成的构造反转,主要沿西湖-基隆大断裂发生反转,并使得断裂带附近地层褶皱、抬升并遭受剥蚀,形成了T20构造不整合界面(表 1)。
中新世晚期,菲律宾海板块以北西西向朝着欧亚板块俯冲,且速率增加,太平洋板块俯冲速率继续加大,使得东亚陆缘继续处于左行压扭应力状态下,随着冲绳海槽的持续扩张,钓鱼岛隆褶带遭受挤压变形和岩浆侵入发生隆升作用并向西挤压[42-47]。中新世末(约10 Ma)西湖凹陷发生第三期构造反转,即龙井运动造成的构造反转。该期反转十分强烈,范围广、幅度大,地层剥蚀量高,在凹陷中东部形成了大规模的地层褶皱、抬升和剥蚀等,以及下部断层的逆向再活动[19, 22, 38, 48],形成了T12构造不整合面(表 1)。中新世末以来,冲绳海槽继续扩张裂陷,向西挤压开始不明显,东海陆架盆地发生整体沉降[49-52]。
前人研究结果显示,西湖凹陷构造反转在时间上具有从玉泉运动(T30)、花港运动(T20)至龙井运动(T12)由老到新相对增强的趋势,空间上具有由南向北逐渐增强的趋势[26]。因此本文数值模拟研究主要集中在构造反转强度较大的中北部地区。
-
为了计算出西湖凹陷的应力和地形特点的时空演化过程,从而讨论其构造反转动力学机制,本文采用了三维有限元模拟的方法,使用有限元软件ANSYS来耦合弹性应变、粘性应变和塑性应变。三维模型的建立主要是基于西湖凹陷中北部已有地震剖面、区域地质成果以及钻井等资料,在横向上,为了减小边界效应对模拟结果的影响,模型设计的范围较大,涵盖了西湖凹陷内部主要的一级和二级构造单元及其周缘地区,包括西部的海礁隆起带和东部的钓鱼岛隆褶带;在纵向上,利用已有地震剖面拟合出了各套地层,从上至下分为6层,其中基底深度最大处约15km,充分考虑了盆地基底横向上的不均匀性,更加真实地拟合出了盆地内各套地层的空间形态(图 3)。
断层的滑移速率以及周边应力场的分布特征同样受到断层的空间几何形态影响[53-57],因此,为了反映东海陆架盆地西湖凹陷内基底和地层的空间不连续性,模型中加入了9条主要断裂,同样采用实际地震剖面数据拟合出其空间形态并加入模型中,包括西湖-基隆大断裂及平湖斜坡大断裂两条主断裂带(图 3)。对于这9条断层接触行为的判定,本文使用了库仑摩擦模型[56, 58-64],断层两盘间的抗剪切应力(τ)可表达为:
$ \tau=\mu p $
(1) 其中,p为断层接触面间的正应力(Pa),μ表示断层接触面间的摩擦系数。根据前人数值模拟结果发现[65-67],大型走滑断层的摩擦系数一般较小,约为0.08,因此,本文模型将其摩擦系数统一设置为0.08。西湖坳陷基底断层上部的沉积盖层多以砂岩或泥岩类为主,断层接触面相对粗糙,因此,摩擦系数设置为0.1。本文的有限元接触算法主要采用了惩罚函数[68, 69],这是为了防止模拟过程中断层两盘的穿透量过大而影响计算结果,其表达式为:
$ f_{n}=k_{n} h_{n} $
(2) 其中,n表示接触面法线方向,fn表示接触力,kn表示惩罚函数,hn表示穿透距离。因为该模型为准静态模型,所以忽略了地震等其他可能因素对断层及其周边地区位移场造成的影响。
前人的岩石和矿物的变形实验研究和大量的野外观察证明,中-深度地壳和上地幔的大多数岩石变形,都是通过组成岩石的矿物塑性变形和流动变形完成,表现在微观尺度上主要为位错蠕变和扩散蠕变两种过程。为了能够正确描述这种塑性流变学行为,本文采用了稳态幂指数流变速率公式[70, 71],其表达式如下:
$ \dot{\varepsilon}=A \sigma^{n} \exp \left(\frac{-Q}{R T}\right) $
(3) 其中,ε是应变率,A是物质结构常数,R是气体常数,Q是活化能,σ是差应力,n是应力指数,上述参数的选取主要是依据前人的研究成果设定(表 2)。
-
关于模型边界条件的设定,主要是根据前人的研究成果[26, 32, 73, 74],西湖凹陷的构造演化主要经历了三期构造反转,考虑了盆地形成现今空间结构的三种观点,因而设定了三种边界条件,各开展了三期构造反转的模拟(图 4)。需要强调的是,每种边界条件施加过程中,前一期模拟结果均作为下一期的初始条件。
第一种模型边界条件强调西湖凹陷构造应力场的分布特征主要受北东向边界断层拉张作用的影响,走滑只起到调节性作用。因此在模型的东南边界设置了早期垂直盆地长轴的伸展,即SE向位移;中晚期垂直盆地走向挤压,即NW向位移的边界条件,其中模型北西边界为固定边界。第二种边界条件认为西湖凹陷为一个张扭性质的盆地,模型边界条件强调西湖凹陷构造应力场的分布特征除主要受北东向边界断层拉张作用的影响外,也受到了走滑作用对凹陷的影响,因此对模型施加了早期右行张扭剪切作用、中晚期左行压扭剪切作用,相对第一种条件来说,增大了走滑分量的影响,即在第一期模型的西边界和东边界分别设置了NNE向和SSW向位移,第二、三期在模型的西边界和东边界分别设置了SSW向和NNE向位移。第三种模型边界条件强调西湖凹陷构造应力场的分布特征明显受北东向边界断层走滑作用的影响,进一步增大了走滑分量对模型中凹陷的影响,因此对模型施加了早期右行剪切作用、中晚期左行剪切作用,即第一期在模型西边界和东边界分别设置了NE向和SW向位移,第二、三期在模型西边界和东边界分别设置了SW向和NE向位移。
-
盆内构造反转的最直接现象就是形成褶皱[75],这种构造现象在模型中反应在垂向位移上。在本文模拟结果中,首先提取了不同模型不同构造期次垂直位移场的演化过程,这是因为垂直位移场能够直观反映盆地的反转样式及过程[19]。实验过程中,本文分别从西湖凹陷3个模型模拟结果中选取了贯穿盆地的3个曲面和2个测线,从而实现测量地壳不均匀性和断裂的发育对地表变形的影响。本文用4个参数来描述褶皱的几何形态(图 5A, 图 6A, 图 7A):波长、波幅、面积和均衡抬升。波长代表褶皱两翼的最大宽度,波幅代表褶皱谷底到波峰的最大距离,面积代表褶皱反转发育的最大范围。
-
图 5为第一种边界条件作用下研究区垂向位移的分布特征,垂向位移高值为红色,低值为蓝色(以下类似图件相同)。第1期SE向伸展过后,西湖凹陷北部整体上为一宽缓隆起,地形变化表现为整体上的沉降,沉降幅度较小;凹陷南部的局部地区表现为受先存断层控制的断阶带或堑-垒构造,基底断陷形成两个隆起带和一个坳陷区(图 5C)。
第2期以第1期的模拟结果作为初始条件,应力环境由伸展变为挤压,此时凹陷表现出整体反转的特征。凹陷北部地形隆升,高度约60m(图 5C),主要发生在凹陷中央隆起带及西部斜坡带,东部断阶带靠近嘉兴构造带区域却表现为地形上的向下挠曲,这可能受到先存断层的控制;凹陷南部垂向上则表现为受中央隆起带先存断层控制的地形隆起,抬升量较大,达110m左右,形成了两个隆起带和一个坳陷区,波长较第1期反转后大。反转隆升中心位于南部中央反转构造带,推测处于宁波构造带附近。
第3期继承第2期的模拟结果,仍为挤压应力,结果与第2期基本一致(图 5C),依然表现为逐渐反转隆升的特征。坳陷北部嘉兴构造带附近由向下挠曲开始逐渐转变为隆升,发生地形的反转隆升;南部的宁波构造带附近反转强度则持续增强。
基于以上模拟结果可见,在该种边界模型下,西湖凹陷中部地区构造反转过程表现为伸展期的整体沉降和挤压期的整体抬升(图 5)。
-
第1期在右行张扭应力环境作用下,西湖凹陷中部的中央隆起带,包括宁波构造带及嘉兴构造带,首先开始抬升,成为局部的抬升中心,抬升量约70m(图 6C)。同时,在凹陷的西部斜坡带和东部断阶带则表现为整体的沉降。由此可见,在该种边界条件作用下凹陷应整体表现为地堑-地垒-地堑的垂向分布特征。而产生这种现象的原因应该是受到了先存东倾中央断裂带和凹陷西部一系列西倾次级断裂带的联合影响。
第2期以第1期的模拟结果作为初始条件,边界条件则由右行张扭环境逐渐转变为与盆地走向呈一定夹角的左行压扭环境。可以看到,在该种边界条件下,凹陷的西部斜坡带和中央隆起带整体发生沉降,沉降幅度约60m(图 6C),并分别形成了一个隆起带和一个坳陷区;而在靠近宁波构造带的东部断阶带则表现为强烈的隆升反转,北部反转抬升量约150m,南部则达到近250 m,成为反转隆升的中心。
第3期继承第2期的模拟结果,在持续性的左行压扭作用下,凹陷的反转过程发生了明显的变化。凹陷的西部斜坡带南北部均形成隆起带,抬升至约200~250m高,而中央隆起带和东部断阶带发生了明显的挠曲沉降过程(图 6C)。反转中心由东部的断阶带逐渐迁移到了西部斜坡带。
在大角度张/压扭转换边界条件下,凹陷的演化过程应表现为:在右行张扭环境下,凹陷由于差异沉降导致了中央隆起带的整体抬升;到张/压扭转换期,在两种不同剪切环境的联合作用下导致了凹陷东部断阶带的反转隆升和西部的沉降;最后,在左行压扭控制下,凹陷发生了西部的反转抬升和东部的整体沉降(图 6)。
-
在第三种模型中,本文继续加大了走滑分量对模型中凹陷的影响。图 7为小角度右行张扭和左行压扭条件作用下研究区垂向位移的分布特征。
第1期小角度右行张扭环境下,介于中央断裂带和西部斜坡带之间的整个西部凹陷带出现了垂向上的明显抬升,抬升量约150m(图 7C),隆升的中心位于两条北东向次级断层之间,反转强度由北至南明显减小;与此同时,介于中央断裂带和西湖基隆大断裂之间的东部凹陷带则出现了明显的沉降过程,沉降量局部可达80m左右,局部区域还出现先存断层控制的堑-垒构造。
第2期凹陷边界条件则由小角度右行张扭作用逐渐变为小角度左行压扭作用。在这种应力转换环境下,介于中央断裂带和西部斜坡断裂带之间的西部凹陷带开始逐渐沉降,垂向位移表现为负值;而东部凹陷带的南部靠近钓鱼岛隆褶带的地区则发生了明显的隆起抬升,抬升量不大,约30 m,呈现由北至南抬升量增大的特征(图 7C)。
第3期凹陷的构造应力场主要受控于左行小角度压扭环境控制之下。在该种边界条件下,凹陷的西部斜坡带进入了整体沉降期,沉降量达到110m左右(图 7C);而东部断阶带,尤其是在中央隆起带,由于受中央断裂带的控制发生了明显的构造抬升,抬升量约195m,即发生了构造反转。隆升中心主要位于中央断裂带的嘉兴构造带及宁波构造带,反转强度由北至南逐渐减小。
在小角度张/压扭应力环境下,凹陷的构造演化过程表现为:在小角度右行张扭环境下,凹陷中央隆起带西部为隆升中心,东部断阶带则发生了区域沉降;到张/压扭转换期,在两种不同剪切环境的联合作用下导致了凹陷东部断阶带的反转隆升和西部的沉降;第3期在左行压扭应力环境控制下,凹陷西部斜坡带进一步沉降,而东部断阶带隆升则进一步加强(图 7C)。
-
不同的构造反转过程会引起盆地内部复杂的构造地貌抬升与沉降,抬升导致地层相应的剥蚀,进而形成与反转构造成因相关的区域角度不整合。构造反转剥蚀量的恢复不仅有利于直接了解盆地构造反转隆升的强度、抬升剥露历史,而且也为研究地层中油气的运移和保存提供了重要依据。因此,本文将对研究区反转剥蚀量进行恢复,与模型反转隆升结果进行对比,从而得到符合实际地质演化的较为合理的边界条件。目前剥蚀量计算的方法很多,根据西湖凹陷的具体情况,本文主要采用了地震剖面显示的地层结构外延计算地层剥蚀厚度,并结合前人的研究成果互相校验,以期达到半定量的剥蚀量估算结果。
图 8为西湖凹陷不整合面剥蚀量估算结果,不同界面剥蚀量大致呈NE向条带状分布,与盆地的挤压反转构造带走向分布一致,反映了反转构造作用对盆地剥蚀界面及剥蚀量分布的控制作用。反转带的隆升时间越早,隆升幅度越大,遭受剥蚀的时间越长,导致剥蚀量越大。T30时期地层剥蚀作用主要发育在凹陷西部斜坡带北端,剥蚀量南部较小,向北增大,而东缘和中央构造带剥蚀量较小。这一不整合界面是玉泉运动的产物,它的形成与裂后的构造格局差异调整作用和热隆起有关,代表了盆地的初次反转。这一时期的剥蚀量分布表明这一期的构造反转主要发生在西部斜坡带,且自北向南反转强度减小。T20时期地层剥蚀作用主要发生在盆地东缘,分布不均一,由南向北变小。这一不整合界面是花港运动的结果,是坳陷期盆地遭受强烈挤压的产物。剥蚀量分布表明这一时期构造反转主要发生在凹陷东部边缘带,且由北向南反转强度增大。T12时期地层剥蚀主要分布在中央反转构造带和东部边缘带,北部剥蚀量大于南部。这一不整合界面是龙井运动的结果,是凹陷发生的第三次大规模构造挤压反转。剥蚀量的分布表明,这一时期发生强烈的挤压构造反转作用,且主要分布在中央反转构造带和东部边缘带,且北部反转强度明显大于南部。
由模型中西湖凹陷模拟结果与西湖凹陷不同界面剥蚀厚度图(图 8)对比明显可以看到,小角度张/压扭剪切模型在第1期构造反转后,中央隆起带西部为隆升区,隆升并遭受剥蚀,东部则为沉降区,沉降并接受沉积,这与西湖凹陷始新世晚期(T30)时期剥蚀区域显示的结果较为一致;同样,第2期构造反转作用后,中央隆起带西部整体为沉降区,东部则发生构造反转,地形隆升,接受剥蚀,这与西湖凹陷在早中新世(T20—T16)时期剥蚀区域显示的结果一致;而在第3期构造反转后,东部隆升进一步加强,剥蚀作用明显增强,与中中新世晚期(T12)时期剥蚀区域显示的结果较为一致。因此,在凹陷中以剪切为主并辅以拉伸或挤压的凹陷演化模式应该是控制凹陷构造反转过程自西向东逐渐发生迁移的主要影响因素。
-
空间中任意3个相互垂直的应力组成了区域的应力状态,其大小和方向决定了区域上的构造变形特征。因此,不同应力及其组合环境下不同类型的构造对区域的变形特征起到重要的控制作用,但单就对盆地中具有流体性质的油气资源来说,压应力的大小可能决定了油气运移的基本方向,即由高势区(高压区)向低势区(低压区)沿主要断裂发生运移[76]。即对于在构造应力作用下沉积盆地中先存断裂对油气运移的控制作用来说,压应力起着主要的控制作用。
图 9为第三种边界条件小角度张/压扭剪切作用下盆地的压应力分布特征与区域已知油气分布对比图。西湖凹陷压应力低值区,主要发育于中央隆起带,其分布范围明显受控于中央断裂带空间形态的变化。在中南部宁波构造带附近,压应力低值区主要位于宁波主断裂的两盘,在其周缘高压应力作用下,宁波构造带是一个很好的油气聚集区;在北部嘉兴构造带附近,压应力低值区位于嘉兴主断裂及其西部次级断裂以及嘉兴主断裂与其东部断裂之间,其周缘地区均具有较高的压应力值,周缘油气很可能会在高流体势的作用下沿着嘉兴主断裂及其次级断裂运移至嘉兴构造带内。对比区域已知油气分布(图 9)可见,在宁波构造带和嘉兴构造带内具有很好的油气显示,这进一步证明了第三种边界条件下西湖凹陷的构造演化可能更加符合该地区的真实情况。但相比于宁波构造带,嘉兴构造带的勘探程度明显较低。因此,本文认为在具有相同的生、储、盖和应力条件下,嘉兴构造带值得重点关注。
-
西湖凹陷构造反转模式主要考虑了3条主要的大型断层(图 10),分别为西部海礁隆起带边界断层、东部钓鱼岛隆褶带的西湖-基隆大断层以及中央隆起带断层。3条断层之间的空间组合关系分别为:西部海礁隆起带边界断层与中央隆起带断层之间为左阶断裂组合关系,而东部边界断层与中央隆起带断层之间则呈现为右阶断裂组合关系,在这种断裂组合特征之下,盆地在T30时期发生右行走滑作用,西部凹陷带的动力学演化特征表现为在右行左阶环境下的地形隆升,遭受剥蚀,并很可能在西部断阶带发生构造反转,成为其反转中心。与此相反,东部凹陷带在右行右阶环境下则表现为局部上的整体沉降;T20—T12时期,边界条件由小角度的右行张扭变为小角度的左行压扭,在这种边界环境转换模式下,反转中心同时也发生迁移,由西部断阶带迁移至东部带,西部由反转变为沉降,东部则由沉降变为反转,地形隆升,遭受剥蚀。由此可见,盆地内部先存断裂格局对盆地反转构造时空迁移具有显著影响,是反转构造动力学必须考虑的重要因素。
-
模型中第1期反转对应于盆地T30时期反转阶段。在太平洋板块NWW向俯冲作用下,西湖凹陷处于右行张扭应力场作用下,中央隆起带西部发生构造隆升,垂向位移表现为正值;东部则发生沉降,局部地区出现先存断层控制的堑-垒构造,对未来反转区的进一步扩展抬升起到了一定的限制性作用。
第2期反转对应于盆地T20时期的反转阶段。太平洋板块俯冲速率持续增强,菲律宾岛弧发生向西仰冲,西湖凹陷处于左旋压扭的应力场作用下,中央隆起带西部区域开始逐渐沉降,垂向位移表现为负值;构造反转向东部区域传播,发生了明显的隆起抬升,由北至南抬升量不断增大。
第3期反转对应于盆地T12时期的反转阶段。菲律宾海板块以NWW向朝着欧亚板块俯冲,导致冲绳海槽持续扩张,钓鱼岛隆褶带隆升并向西挤压,西湖凹陷处于左旋压扭应力场作用下,产生三期反转中的最大反转强度。中央隆起带西部区域进入整体沉降期,沉降幅度较大;东部区域,尤其是中央隆起带北部,由于受到先存中央断裂带的控制发生了明显的构造抬升,产生了一套新的断裂-逆冲断裂系统,并伴随一批大型挤压背斜的形成(图 10)。反转由北至南传播,强度逐渐减小,这一点在地震剖面上也很明显。
-
基于有限元数值模拟技术,结合西湖凹陷区域地质特征,通过对西湖凹陷反转过程的三维数值模拟,本文得到了晚白垩世以来西湖凹陷反转构造的动态演化机制。尽管模型还存在一些问题,但模拟结果为深刻理解西湖凹陷变形机制提供了基础。本文得到了如下结论:
(1) 模拟结果显示,小角度张/压扭模型模拟结果与西湖凹陷不同界面剥蚀厚度揭示的反转过程一致,因此,凹陷中以剪切为主并辅以拉伸或挤压的凹陷演化模式是控制凹陷构造反转过程自西向东逐渐发生迁移的主要影响因素。
(2) 西湖凹陷构造反转模式主要考虑了3条主要的大型断层,表现为T30时期的小角度的右行张扭到T20—T12时期的小角度的左行压扭边界条件的转换,导致反转中心由西部断阶带迁移至东部边缘带。
(3) 西湖凹陷的中北部反转构造动力学机制主要受太平洋板块、菲律宾海板块向欧亚板块的俯冲作用控制。由于俯冲速率及方向的改变,导致区域构造应力场由右行张扭转变为左行压扭,控制了西湖凹陷构造反转的形成。
NUMERICAL DYNAMIC MODELING OF TECTONIC INVERSION IN THE NORTHEASTERN XIHU SAG
More Information-
摘要: 东海陆架盆地位于欧亚板块东缘,受太平洋板块俯冲、印度板块与欧亚板块碰撞影响,经历了3期构造反转作用,尤其是位于其东部拗陷带的西湖凹陷,记录了新生代以来的构造反转过程。选择反转作用最为强烈的西湖凹陷中北部区域,进行有限元数值模拟,以此为基础来探讨反转构造的动力学机制。在大量地震剖面的综合分析的基础上,基于西湖凹陷的构造几何学特征,建立了具有黏弹性性质的6套地层以及9条主要断裂组成的地质模型。模拟结果表明,在西湖凹陷中以剪切为主并辅以早期拉伸、晚期挤压的凹陷演化模式是控制凹陷构造反转过程自西向东逐渐发生迁移的主要因素。西湖凹陷垂向位移场演化特征的定量化分析表明,反转构造过程表现为第一期反转导致西部隆升、东部沉降,第二、三期反转导致西部沉降、中东部隆升的特征,模拟结果与西湖凹陷不同界面剥蚀厚度显示的区域隆升及沉降演化过程一致,且模型中压性应力场低值区的分布与已知油气藏的分布也具有较好的对应关系。因此,西湖凹陷构造反转过程主要受区域应力场由早期的右行张扭逐渐调整为晚期的左行压扭影响。造成这种调整的原因主要是太平洋板块、菲律宾海板块向欧亚板块的俯冲速率及方向的改变。Abstract: The East China Sea Continental Shelf Basin is a large basin located in the eastern border of the Eurasian Plate. Under the effects of the India-Eurasia collision and the subduction of the Pacific Plate, this region has experienced three-stages of tectonic inversions since Cenozoic, especially in the Xuhu Sag. This paper presents a numerical simulation results by finite element method to the north-central Xihu Sag in an attempt to investigate the mechanism of tectonic inversion. Based on the comprehensive structural analysis and the structural geometry of the sag established with a large number of seismic profiles, a viscoelastic geologic model of 6 layers and 9 major faults was constructed for the Xihu Sag. Simulation results show that the subsidence caused by the shearing movement characterized by transtension in early stage and transpression in later stage is the major force controling the inversion process and its propagation from west to east. Based on quantitative analysis of vertical displacement field of the Sag, this paper identified a tectonic inversion process with uplifting in the west and subsiding in the east after the first-stage of inversion; and then subsiding in the west and uplifting in the middle and east after the second and third-stages. The simulation results are well consistent with the process of uplifting or subsidence obtained from the denudation thickness and the distribution of low compression stress, displaying certain correlation with the distribution of oil reservoirs. The results indicate the plausibility of similar deformation controls. In conclusion, the formation of the tectonic inversion in the Xihu Sag is controlled by the adjustment of the stress field from dextral transtension to sinistral transpression caused by the change in subduction rates and direction of the Pacific Plate and the Philippine Sea Plate.
-
表 1 西湖凹陷地层柱及构造演化
Table 1. Stratigraphic column and tectonic evolution of the Xihu Sag
-
[1] Cooper M A, Williams G D. Inversion Tectonics[M]. London: Geological Society by Blackwell Scientific Publications, 1989. [2] Buchanan P G, Mcclay K R. Sandbox experiments of inverted listric and planar fault systems[J]. Tectonophysics, 1991, 188(1-2): 97-115. doi: 10.1016/0040-1951(91)90317-L [3] 姚超, 焦贵浩, 王同和, 等.中国含油气构造样式[M].北京:石油工业出版社, 2004. YAO Chao, JIAO Guihao, WANG Tonghe, et al. Structural Styles of Oil-Bearing Basins in China[M]. Beijing: Petroleum Industry Press, 2004. [4] 陈昭年, 陈发景.反转构造与油气圈闭[J].地学前缘, 1995, 2(3-4): 96-102.http://www.earthsciencefrontiers.net.cn/CN/article/searchArticle.do CHEN Zhaonian, CHEN Fajing. Inversion structures and their relationship to traps of oil and gas[J]. Earth Science Frontiers, 1995, 2(3-4): 96-102.http://www.earthsciencefrontiers.net.cn/CN/article/searchArticle.do [5] 胡望水, 刘学锋, 吕新华, 等.论正反转构造的分类[J].新疆石油地质, 2000, 21(1): 5-8. doi: 10.3969/j.issn.1001-3873.2000.01.002 HU Wangshui, LIU Xuefeng, LV Xinhua, et al. On the classifications of positive reversed structures[J]. Xinjiang Petroleum Geology, 2000, 21(1): 5-8. doi: 10.3969/j.issn.1001-3873.2000.01.002 [6] 胡望水, 吕炳全, 毛治国, 等.中国东部中新生代含油气盆地的反转构造[J].同济大学学报, 2004, 32(2): 182-186. doi: 10.3321/j.issn:0253-374X.2004.02.009 HU Wangshui, LV Bingquan, MAO Zhiguo, et al. Inversion structure characteristic of petroleum basin in Mesozoic and Cenozoic in middle and East China[J]. Journal of Tongji University, 2004, 32(2): 182-186. doi: 10.3321/j.issn:0253-374X.2004.02.009 [7] 李军生, 林春明.反转背斜构造自生自储油藏成藏模式[J].石油学报, 2006, 27(2): 34-37. doi: 10.3321/j.issn:0253-2697.2006.02.007 LI Junsheng, LIN Chunming. Forming pattern of self-generating and self-preserving reservoirs in reversal anticline structure[J]. Acta Petrolei Sinica, 2006, 27(2): 34-37. doi: 10.3321/j.issn:0253-2697.2006.02.007 [8] 胡望水, 李瑞升, 李涛, 等.正反转构造动力学成因探讨[J].石油天然气学报, 2007, 29(4): 23-27. doi: 10.3969/j.issn.1000-9752.2007.04.004 HU Wangshui, LI Ruisheng, LI Tao, et al. Discussion on dynamic genesis of positive inversion structure[J]. Journal of Oil and Gas Technology, 2007, 29(4): 23-27. doi: 10.3969/j.issn.1000-9752.2007.04.004 [9] 杨风丽, 王敏雪, 庄建建, 等.西湖凹陷反转构造定量运动学过程及对油气的控制作用[J].石油学报, 2010, 31(4): 596-601. doi: 10.7623/syxb201004013 YANG Fengli, WANG Minxue, ZHUANG Jianjian, et al. Kinematical processes of inversion structure and its contribution to hydrocarbon accumulation in Xihu Depression of East China Sea Basin[J]. Acta Petrolei Sinica, 2010, 31(4): 596-601. doi: 10.7623/syxb201004013 [10] 徐子英, 孙珍, 张云帆, 等.南海北部陆缘盆地反转构造研究:以莺歌海盆地临高隆起和珠江口盆地琼海凹陷为例[J].地学前缘, 2010, 17(4): 90-98.http://www.earthsciencefrontiers.net.cn/CN/article/searchArticle.do XU Ziying, SUN Zhen, ZHANG Yunfan, et al. Inversion structures in the northern continental margin of the South China Sea: taking Lingao Uplift in the Yinggehai Basin and Qionghai Sag in the Pearl River Mouth Basin as examples[J]. Earth Science Frontiers, 2010, 17(4): 90-98.http://www.earthsciencefrontiers.net.cn/CN/article/searchArticle.do [11] Sibson R H. A note on fault reactivation[J]. Journal of Structural Geology, 1985, 7(6): 751-754. doi: 10.1016/0191-8141(85)90150-6 [12] Sandiford M. Mechanics of basin inversion[J]. Tectonophysics, 1999, 305(1-3): 109-120. doi: 10.1016/S0040-1951(99)00023-2 [13] Buiter S J H, Pfiffner O A. Numerical models of the inversion of half-graben basins[J]. Tectonics, 2003, 22(5), doi: 10.1029/2002TC001417. [14] Panien M, Schreurs G, Pfiffner A. Sandbox experiments on basin inversion: testing the influence of basin orientation and basin fill[J]. Journal of Structural Geology, 2005, 27(3): 433-445. doi: 10.1016/j.jsg.2004.11.001 [15] Ventisette D C, Montanari D, Sani F, et al. Basin inversion and fault reactivation in laboratory experiments[J]. Journal of Structural Geology, 2006, 28(11): 2067-2083. doi: 10.1016/j.jsg.2006.07.012 [16] Marques F O, Nogueira C R. Normal fault inversion by orthogonal compression: Sandbox experiments with weak faults[J]. Journal of Structural Geology, 2008, 30(6): 761-766. doi: 10.1016/j.jsg.2008.02.015 [17] Buiter S J H, Pfiffner O A, Beaumont C. Inversion of extensional sedimentary basins: a numerical evaluation of the localisation of shortening[J]. Earth and Planetary Science Letters, 2009, 288(3-4): 492-504. doi: 10.1016/j.epsl.2009.10.011 [18] Gomes C J S, Danderfer Filho A, Posada A M A, et al. The role of backstop shape during inversion tectonics physical models[J]. Anais da Academia Brasileira de Ciências, 2010, 82(4): 997-1012. doi: 10.1590/S0001-37652010000400021 [19] Dai L M, Li S Z, Lou D, et al. Numerical modeling of Late Miocene tectonic inversion in the Xihu Sag, East China Sea Shelf Basin, China[J]. Journal of Asian Earth Sciences, 2014, 86: 25-37. doi: 10.1016/j.jseaes.2013.09.033 [20] Dai L M, Li Q W, Li S Z, et al. Numerical modelling of stress fields and earthquakes jointly controlled by NE- and NW-trending fault zones in the Central North China Block[J]. Journal of Asian Earth Sciences, 2015, 114: 28-40. doi: 10.1016/j.jseaes.2015.05.021 [21] 张建培, 张涛, 刘景彦, 等.西湖凹陷反转构造分布与样式[J].海洋石油, 2008, 28(4): 14-20. doi: 10.3969/j.issn.1008-2336.2008.04.003 ZHANG Jianpei, ZHANG Tao, LIU Jingyan, et al. Distribution and style of inversed structures in Xihu Depression[J]. Offshore Oil, 2008, 28(4): 14-20. doi: 10.3969/j.issn.1008-2336.2008.04.003 [22] Yang F L, Yu H X, Zhang Q L, et al. Correlations between shortening rate, uplift rate, and inversion rate in central inversion zone of Xihu depression, East China Sea Basin[J]. Journal of Earth Science, 2009, 20(4): 699-708. doi: 10.1007/s12583-009-0056-6 [23] Zhang G H, Li S Z, Suo Y H, et al. Cenozoic positive inversion tectonics and its migration in the East China Sea Shelf Basin[J]. Geological Journal, 2016, 51(S1): 176-187. doi: 10.1002/gj.2809 [24] 张田, 张建培, 张绍亮, 等.有限元数值模拟技术在西湖凹陷中央反转构造带形成机制研究中的应用[J].海洋石油, 2012, 32(4): 11-16. doi: 10.3969/j.issn.1008-2336.2012.04.011 ZHANG Tian, ZHANG Jianpei, ZHANG Shaoliang, et al. Application of the finite element numerical simulation method to the study on formation mechanism of the central inverted structure zone, Xihu Depression[J]. Offshore Oil, 2012, 32(4): 11-16. doi: 10.3969/j.issn.1008-2336.2012.04.011 [25] 张绍亮, 张建培, 唐贤君, 等.东海西湖凹陷断裂系统几何学特征及其成因机制[J].海洋地质与第四纪地质, 2014, 34(1): 87-94.http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hydzydsjdz201401010 ZHANG Shaoliang, ZHANG Jianpei, TANG Xianjun, et al. Geometry characteristic of the fault system in Xihu Sag in East China Sea and its formation mechanism[J]. Marine Geology and Quaternary Geology, 2014, 34(1): 87-94.http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hydzydsjdz201401010 [26] 张国华, 张建培.东海陆架盆地构造反转特征及成因机制探讨[J].地学前缘, 2015, 22(1): 260-270. doi: 10.13745/j.esf.2015.01.022 ZHANG Guohua, ZHANG Jianpei. A discussion on the tectonic inversion and its genetic mechanism in the East China Sea Shelf Basin[J]. Earth Science Frontiers, 2015, 22(1): 260-270. doi: 10.13745/j.esf.2015.01.022 [27] Zhang J P, Li S Z, Suo Y H. Formation, tectonic evolution and dynamics of the East China Sea Shelf Basin[J]. Geological Journal, 2016, 51(S1): 162-175. doi: 10.1002/gj.2808 [28] 刘景彦, 林畅松, 姜亮, 等.东海西湖凹陷第三系反转构造及其对油气聚集的影响[J].地球学报, 2000, 21(4): 350-355. doi: 10.3321/j.issn:1006-3021.2000.04.003 LIU Jingyan, LIN Changsong, JIANG Liang, et al. Characteristics of tertiary inversion structures and their influence on oil-gas accumulation in Xihu Trough, East China Sea[J]. Acta Geoscientia Sinica, 2000, 21(4): 350-355. doi: 10.3321/j.issn:1006-3021.2000.04.003 [29] 陈志勇, 葛和平.西湖凹陷反转构造与油气聚集[J].中国海上油气(地质), 2003, 17(1): 20-24.http://d.old.wanfangdata.com.cn/Periodical/zghsyq-dz200301004 CHEN Zhiyong, GE Heping. Inversion structures and hydrocarbon accumulation in Xihu Sag, East China Sea Basin[J]. China Offshore Oil and Gas (Geology), 2003, 17(1): 20-24.http://d.old.wanfangdata.com.cn/Periodical/zghsyq-dz200301004 [30] 张敏强, 钟志洪, 夏斌, 等.东海西湖凹陷中南部晚中新世构造反转与油气运聚[J].中国海上油气, 2005, 17(2): 73-79. doi: 10.3969/j.issn.1673-1506.2005.02.001 ZHANG Minqiang, ZHONG Zhihong, XIA Bin, et al. Late Miocene tectonic inversion and hydrocarbon migration and accumulation in central and southern Xihu sag, East China Sea[J]. China Offshore Oil and Gas, 2005, 17(2): 73-79. doi: 10.3969/j.issn.1673-1506.2005.02.001 [31] Ren J Y, Tamaki K, Li S T, et al. Late Mesozoic and Cenozoic rifting and its dynamic setting in Eastern China and adjacent areas[J]. Tectonophysics, 2002, 344(3-4): 175-205. doi: 10.1016/S0040-1951(01)00271-2 [32] 索艳慧, 李三忠, 戴黎明, 等.东亚及其大陆边缘新生代构造迁移与盆地演化[J].岩石学报, 2012, 28(8): 2602-2618.http://www.ysxb.ac.cn/ysxb/ch/reader/key_query.aspx SUO Yanhui, LI Sanzhong, DAI Liming, et al. Cenozoic tectonic migration and basin evolution in East Asia and its continental margins[J]. Acta Petrologica Sinica, 2012, 28(8): 2602-2618.http://www.ysxb.ac.cn/ysxb/ch/reader/key_query.aspx [33] Li S Z, Zhao G C, Dai L M, et al. Cenozoic faulting of the Bohai Bay Basin and its bearing on the destruction of the eastern North China Craton[J]. Journal of Asian Earth Sciences, 2012, 47: 80-93. doi: 10.1016/j.jseaes.2011.06.011 [34] Northrup C J, Royden L H, Burchfiel B C. Motion of the Pacific plate relative to Eurasia and its potential relation to Cenozoic extension along the eastern margin of Eurasia[J]. Geology, 1995, 23(8): 719-722. doi: 10.1130/0091-7613(1995)023<0719:MOTPPR>2.3.CO;2 [35] 郑求根, 周祖翼, 蔡立国, 等.东海陆架盆地中新生代构造背景及演化[J].石油与天然气地质, 2005, 26(2): 197-201. doi: 10.3321/j.issn:0253-9985.2005.02.012 ZHENG Qiugen, ZHOU Zuyi, CAI Liguo, et al. Meso-Cenozoic tectonic setting and evolution of East China Sea shelf basin[J]. Oil & Gas Geology, 2005, 26(2): 197-201. doi: 10.3321/j.issn:0253-9985.2005.02.012 [36] Li C F, Zhou Z Y, Ge H P, et al. Correlations between erosions and relative uplifts from the central inversion zone of the Xihu depression, East China Sea Basin[J]. Terrestrial, Atmospheric and Oceanic Sciences, 2007, 18(4): 757-776. doi: 10.3319/TAO.2007.18.4.757(TT) [37] Li S Z, Zhao G C, Dai L M, et al. Mesozoic basins in eastern China and their bearing on the deconstruction of the North China Craton[J]. Journal of Asian Earth Sciences, 2012, 47: 64-79. doi: 10.1016/j.jseaes.2011.06.008 [38] Suo Y H, Li S Z, Yu S, et al. Cenozoic tectonic jumping and implications for hydrocarbon accumulation in basins in the East Asia Continental Margin[J]. Journal of Asian Earth Sciences, 2014, 88: 28-40. doi: 10.1016/j.jseaes.2014.02.019 [39] 杨文达, 崔征科, 张异彪.东海地质与矿产[M].北京:海洋出版社, 2010. YANG Wenda, CUI Zhengke, ZHANG Yibiao, et al. Geology and Mineral Resources in East China Sea[M]. Beijing: Ocean Press, 2010. [40] 徐发.东海陆架盆地新生界结构特征及迁移规律[J].石油天然气学报, 2012, 34(6): 1-7. doi: 10.3969/j.issn.1000-9752.2012.06.001 XU Fa. Characteristics of Caenozoic structure and tectonic migration of the East China Sea Shelf Basin[J]. Journal of Oil and Gas Technology, 2012, 34(6): 1-7. doi: 10.3969/j.issn.1000-9752.2012.06.001 [41] 郭令智, 钟志洪, 王良书, 等.莺歌海盆地周边区域构造演化[J].高校地质学报, 2001, 7(1): 1-12. doi: 10.3969/j.issn.1006-7493.2001.01.001 GUO Lingzhi, ZHONG Zhihong, WANG Liangshu, et al. Regional tectonic evolution around Yinggehai Basin of South China Sea[J]. Geological Journal of China Universities, 2001, 7(1): 1-12. doi: 10.3969/j.issn.1006-7493.2001.01.001 [42] 臧绍先, 宁杰远.菲律宾海板块与欧亚板块的相互作用及其对东亚构造运动的影响[J].地球物理学报, 2002, 45(2): 188-197. doi: 10.3321/j.issn:0001-5733.2002.02.005 ZANG Shaoxian, NING Jieyuan. Interaction between Philippine Sea Plate (PH) and Eurasia (EU) Plate and its influence on the Movement Eastern Asia[J]. Chinese Journal of Geophysics, 2002, 45(2): 188-197. doi: 10.3321/j.issn:0001-5733.2002.02.005 [43] 赵金海, 唐建, 薄玉玲.关于东海新生代盆地油气勘探的若干见解[J].中国海上油气, 2003, 17(1): 14-19.http://d.old.wanfangdata.com.cn/Periodical/zghsyq-dz200301003 ZHAO Jinhai, TANG Jian, BO Yuling. Some opinions on petroleum exploration in Cenozoic East China Sea Basin[J]. China Offshore Oil and Gas, 2003, 17(1): 14-19.http://d.old.wanfangdata.com.cn/Periodical/zghsyq-dz200301003 [44] 丁巍伟, 杨树锋, 陈汉林, 等.台湾岛以南海域新近纪的弧-陆碰撞造山作用[J].地质科学, 2006, 41(2): 195-201. doi: 10.3321/j.issn:0563-5020.2006.02.002 DING Weiwei, YANG Shufeng, CHEN Hanlin, et al. Arc-continent collision orogeny in offshore Taiwan during Neogene[J]. Chinese Journal of Geology, 2006, 41(2): 195-201. doi: 10.3321/j.issn:0563-5020.2006.02.002 [45] 王鹏, 赵志刚, 张功成, 等.东海盆地钓鱼岛隆褶带构造演化分析及对西湖凹陷油气勘探的意义[J].地质科技情报, 2011, 30(4): 65-72. doi: 10.3969/j.issn.1000-7849.2011.04.009 WANG Peng, ZHAO Zhigang, ZHANG Gongcheng, et al. Analysis on structural evolution in Diaoyu Islands Folded-Uplift Belt, East China Sea Basin and its impact on the hydrocarbon exploration in Xihu Sag[J]. Geological Science and Technology Information, 2011, 30(4): 65-72. doi: 10.3969/j.issn.1000-7849.2011.04.009 [46] 徐杰, 周本刚, 计凤桔, 等.中国东部海域及其邻区现代构造应力场研究[J].地学前缘, 2012, 19(4): 1-7.http://www.earthsciencefrontiers.net.cn/CN/article/searchArticle.do XU Jie, ZHOU Bengang, JI Fengju, et al. The recent tectonic stress field of offshore of China mainland and adjacent areas[J]. Earth Science Frontiers, 2012, 19(4): 1-7.http://www.earthsciencefrontiers.net.cn/CN/article/searchArticle.do [47] 郭真, 刘池洋, 田建锋.东海陆架盆地龙井运动构造影响及其发育背景[J].西北大学学报:自然科学版, 2015, 45(5): 801-810. doi: 10.16152/j.cnki.xdxbzr.2015-05-022 GUO Zhen, LIU Chiyang, TIAN Jianfeng. Longjing Movement structural effect and developmental background in East China Sea Basin[J]. Journal of Northwest University: Natural Science Edition, 2015, 45(5): 801-810. doi: 10.16152/j.cnki.xdxbzr.2015-05-022 [48] Zhou Z Y, Jiang J Y, Liao Z T, et al. Basin inversion in Xihu Depression, East China Sea[J]. Gondwana Research, 2001, 4(4): 844-845. doi: 10.1016/S1342-937X(05)70627-4 [49] 许浚远, 张凌云.西北太平洋边缘新生代盆地成因(下):后裂谷期构造演化[J].石油与天然气地质, 2000, 21(4): 287-292. doi: 10.3321/j.issn:0253-9985.2000.04.001 XU Junyuan, ZHANG Lingyun. Genesis of Cenozoic basins in northwest Pacific Ocean Margin (3): tectonic evolution of post rifting period[J]. Oil & Gas Geology, 2000, 21(4): 287-292. doi: 10.3321/j.issn:0253-9985.2000.04.001 [50] 吴根耀, 矢野孝雄.东亚大陆边缘的构造格架及其中-新生代演化[J].地质通报, 2007, 26(7): 787-800. doi: 10.3969/j.issn.1671-2552.2007.07.002 WU Genyao, YANO T. Tectonic framework and Meso-Cenozoic evolution of the East Asian continental margin[J]. Geological Bulletin of China, 2007, 26(7): 787-800. doi: 10.3969/j.issn.1671-2552.2007.07.002 [51] 李三忠, 余珊, 赵淑娟, 等.东亚大陆边缘的板块重建与构造转换[J].海洋地质与第四纪地质, 2013, 33(3): 65-94.http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hydzydsjdz201303008 LI Sanzhong, YU Shan, ZHAO Shujuan, et al. Tectonic transition and plate reconstructions of the East Asian Continental Margin[J]. Marine Geology and Quaternary Geology, 2013, 33(3): 65-94.http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hydzydsjdz201303008 [52] 温志新, 童晓光, 张光亚, 等.全球板块构造演化过程中五大成盆期原型盆地的形成、改造及叠加过程[J].地学前缘, 2014, 21(3): 26-37. doi: 10.13745/j.esf.2014.03.004 WEN Zhixin, TONG Xiaoguang, ZHANG Guangya, et al. The transformation and stacking process of prototype basin in five global plate tectonic evolution stages[J]. Earth Science Frontiers, 2014, 21(3): 26-37. doi: 10.13745/j.esf.2014.03.004 [53] Chester F M, Chester J S. Stress and deformation along wavy frictional faults[J]. Journal of Geophysical Research, 2000, 105(B10): 23421-23430. doi: 10.1029/2000JB900241 [54] Li Q S, Liu M. Geometrical impact of the San Andreas Fault on stress and seismicity in California[J]. Geophysical Research Letters, 2006, 33(8): L08302. doi: 10.1029/2005gl025661 [55] 朱守彪, 张培震. 2008年汶川M S8.0地震发生过程的动力学机制研究[J].地球物理学报, 2009, 52(2): 418-427. ZHU Shoubiao, ZHANG Peizhen. A study on the dynamical mechanisms of the Wenchuan M S8.0 earthquake, 2008[J]. Chinese Journal of Geophysics, 2009, 52(2): 418-427. [56] 戴黎明, 李三忠, 陶春辉, 等.俯冲带耦合作用对苏门答腊地区应变场影响的三维数值模拟[J].地球物理学报, 2010, 53(8): 1837-1851. doi: 10.3969/j.issn.0001-5733.2010.08.009 DAI Liming, LI Sanzhong, TAO Chunhui, et al. 3D numerical modeling of strain field in Sumatra area influenced by the coupling effect of subduction zone[J]. Chinese Journal of Geophysics, 2010, 53(8): 1837-1851. doi: 10.3969/j.issn.0001-5733.2010.08.009 [57] Marshall S T, Cooke M L, Owen S E. Effects of nonplanar fault topology and mechanical interaction on fault-slip distributions in the Ventura Basin, California[J]. Bulletin of the Seismological Society of America, 2008, 98(3): 1113-1127. doi: 10.1785/0120070159 [58] He J K, Chéry J. Slip rates of the Altyn Tagh, Kunlun and Karakorum faults (Tibet) from 3D mechanical modeling[J]. Earth and Planetary Science Letters, 2008, 274(1-2): 50-58. doi: 10.1016/j.epsl.2008.06.049 [59] Wang J, Ye Z R, He J K. Three-dimensional mechanical modeling of large-scale crustal deformation in China constrained by the GPS velocity field[J]. Tectonophysics, 2008, 446(1-4): 51-60. doi: 10.1016/j.tecto.2007.11.006 [60] 朱守彪, 邢会林, 谢富仁, 等.地震发生过程的有限单元法模拟-以苏门答腊俯冲带上的大地震为例[J].地球物理学报, 2008, 51(2): 460-468. doi: 10.3321/j.issn:0001-5733.2008.02.018 ZHU Shoubiao, XING Huilin, XIE Furen, et al. Simulation of earthquake processes by finite element method: the case of megathrust earthquakes on the Sumatra subduction zone[J]. Chinese Journal of Geophysics, 2008, 51(2): 460-468. doi: 10.3321/j.issn:0001-5733.2008.02.018 [61] Zhu S B, Zhang P Z. Numeric modeling of the strain accumulation and release of the 2008 Wenchuan, Sichuan, China, Earthquake[J]. Bulletin of the Seismological Society of America, 2010, 100(5B): 2825-2839. doi: 10.1785/0120090351 [62] Zhu S B, Zhang P Z. FEM simulation of interseismic and coseismic deformation associated with the 2008 Wenchuan Earthquake[J]. Tectonophysics, 2013, 584: 64-80. doi: 10.1016/j.tecto.2012.06.024 [63] 戴黎明, 李三忠, 陶春辉, 等.印度板块挤压驱动龙门山断裂带活动的三维数值模型[J].地球物理学进展, 2011, 26(1): 41-51. doi: 10.3969/j.issn.1004-2903.2011.01.004 DAI Liming, LI Sanzhong, TAO Chunhui, et al. Three-dimensional numberical modeling of activity of the Longmenshan fault zone driven by the India Plate[J]. Progress in Geophysics, 2011, 26(1): 41-51. doi: 10.3969/j.issn.1004-2903.2011.01.004 [64] 戴黎明, 李三忠, 楼达, 等.渤海湾盆地黄骅坳陷应力场的三维数值模拟分析[J].地球物理学报, 2013, 56(3): 929-942. doi: 10.6038/cjg20130321 DAI Liming, LI Sanzhong, LOU Da, et al. Numerical modeling on the stress field in the Huanghua depression, Bohai Bay basin[J]. Chinese Journal of Geophysics, 2013, 56(3): 929-942. doi: 10.6038/cjg20130321 [65] Mount V S, Suppe J. State of stress near the San Andreas fault: implications for wrench tectonics[J]. Geology, 1987, 15(12): 1143-1146. doi: 10.1130/0091-7613(1987)15<1143:SOSNTS>2.0.CO;2 [66] Kong X, Bird P. Neotectonics of Asia: Thin-shell finite-element models with faults[M]//Yin A, Harrison T M. The Tectonic Evolution of Asia. Cambridge: Cambridge University Press, 1996: 18-34. [67] Wang J, Ye Z R. Dynamic modeling for crustal deformation in China: comparisons between the theoretical prediction and the recent GPS data[J]. Physics of the Earth and Planetary Interiors, 2006, 155(3-4): 201-207. doi: 10.1016/j.pepi.2005.11.003 [68] 朱守彪, 张培震, 石耀霖.华北盆地强震孕育的动力学机制研究[J].地球物理学报, 2010, 53(6): 1409-1417. doi: 10.3969/j.issn.0001-5733.2010.06.019 ZHU Shoubiao, ZHANG Peizhen, SHI Yaolin. A study on the mechanisms of strong earthquake occurrence in the North China Basin[J]. Chinese Journal of Geophysics, 2010, 53(6): 1409-1417. doi: 10.3969/j.issn.0001-5733.2010.06.019 [69] 戴黎明, 李三忠, 刘泽, 等.车镇凹陷大1断层对区域油气成藏控制作用的三维数值模拟[J].大地构造与成矿学, 2016, 40(1): 47-57. doi: 10.16539/j.ddgzyckx.2016.01.005 DAI Liming, LI Sanzhong, LIU Ze, et al. 3D numerical modeling of oil and gas accumulation process controlled by the Da 1 Fault, Chezhen Depression[J]. Geotectonica et Metallogenia, 2016, 40(1): 47-57. doi: 10.16539/j.ddgzyckx.2016.01.005 [70] Karato S I, Wu P. Rheology of the upper mantle: a synthesis[J]. Science, 1993, 260(5109): 771-778. doi: 10.1126/science.260.5109.771 [71] Freed A M, Bürgmann R. Evidence of power-law flow in the Mojave desert mantle[J]. Nature, 2004, 430(6999): 548-551. doi: 10.1038/nature02784 [72] Franco D, Govers R, Wortel R. Numerical comparison of different convergent plate contacts: subduction channel and subduction fault[J]. Geophysical Journal International, 2007, 171(1): 435-450. doi: 10.1111/j.1365-246X.2006.03498.x [73] 张建培, 张田, 唐贤君.东海陆架盆地类型及其形成的动力学环境[J].地质学报, 2014, 88(11): 2033-2043. doi: 10.3969/j.issn.0001-5717.2014.11.002 ZHANG Jianpei, ZHANG Tian, TANG Xianjun. Basin type and dynamic environment in the East China Sea Shelf Basin[J]. Acta Geologica Sinica, 2014, 88(11): 2033-2043. doi: 10.3969/j.issn.0001-5717.2014.11.002 [74] Suo Y H, Li S Z, Zhao S J, et al. Continental margin basins in East Asia: tectonic implications of the Meso-Cenozoic East China Sea pull-apart basins[J]. Geological Journal, 2015, 50(2): 139-156. doi: 10.1002/gj.2535 [75] Jarosinski M, Beekman F, Matenco L, et al. Mechanics of basin inversion: finite element modelling of the Pannonian Basin System[J]. Tectonophysics, 2011, 502(1-2): 121-145. doi: 10.1016/j.tecto.2009.09.015 [76] 于来刚.东营凹陷油气运移数值模拟研究[D].中国石油大学硕士学位论文, 2007. YU Laigang. The FEM simulating of hydrocarbon migration in Dongying Depression[D]. Master's Thesis of China University of Petroleum, 2007. [77] 宋小勇, 储呈林, 芮志峰.东海盆地西湖凹陷构造样式及其演化[J].高校地质学报, 2010, 16(1): 86-93. doi: 10.3969/j.issn.1006-7493.2010.01.013 SONG Xiaoyong, CHU Chenglin, RUI Zhifeng. Structural framework and evolution of Xihu Sag in East China sea basin[J]. Geological Journal of China Universities, 2010, 16(1): 86-93. doi: 10.3969/j.issn.1006-7493.2010.01.013 -