裂谷盆地构造－热演化模拟中几个问题的讨论

刘琼颖，何丽娟

1中国科学院地质与地球物理研究所岩石圈演化国家重点实验室，北京 100029

2中国科学院大学，北京 100049

1 引言

沉积盆地构造－热演化模拟是盆地模拟的主要内容之一.构造沉降史与热流演化史构成其两大核心研究内容.拉张模型在描述拉张盆地沉降和热流演化方面取得了极大的成功，实现了构造和热的完美结合（何丽娟和汪集暘，2007）.

McKenzie（1978）提出了一维瞬时均匀拉张模型.该模型认为，岩石圈的拉张在瞬时完成，由于岩石圈减薄，引起热的软流圈物质被动上涌，导致岩石圈被加热，在Airy均衡作用下形成初始沉降.拉张结束后，岩石圈进入热冷却阶段，地表在初始沉降的基础上继续发生热沉降.该模型是运动学模型的代表模型，能很好地预测盆地主要观测数据，为以后研究拉张盆地构造热演化奠定了基础.考虑到岩石圈拉张不可能瞬时完成，学者们进一步提出有限匀速拉张模型（Jarvis and McKenzie，1980）.何丽娟（1999）提出了匀减速拉张速率模式，更接近实际情况.一些学者在均匀拉张模型的基础上提出了纵向非均匀拉张模型，如双层非均匀拉张模型（Royden and Keen，1980；Hellinger and Sclater，1983；Dehler et al.，1997）、与深度相关的拉张模型（Rowley and Sahagian，1986；Davis and Kusznir，2004；Kusznir and Karner，2007；Rosenbaum et al.，2008）.上述拉张模型多假设岩石圈拉张变形为纯剪变形，但实际上变形并非总是对称的.Wernicke（1985）提出了简单剪切拉张模型，但迄今研究者并未观测到切穿整个岩石圈的断裂.为此，学者们用耦合的简单剪切／纯剪切模型来描述不对称拉张的情形，具有一定的实用性（Kusznir and Ziegler，1992；Clift et al.，2002；Baur et al.，2010；Chen et al.，2013）.一些盆地在地质历史时期曾经历过多期拉张，且具有继承性，多期拉张模型由此被提出（Wang et al.，1996；Dehler et al.，1997；He et al.，2002；Lin et al.，2002）.White（1994）以及 White和Bellingham（2002）提出了利用沉降数据反演岩石圈拉张期间一维和二维应变速率的方法，该方法不需要知道岩石圈拉张持续时间、拉张期次等先验信息，在一些地区得到了较好的应用（Xie et al.，2006；Shillington et al.，2008；Song et al.2010；Chen，2014）.考虑到盆地拉张过程中可能的构造反转，研究者们进一步建立了适用于多期伸展和挤压交替的复杂过程的数值计算模型（Negredo et al.，1999；He and Wang，2004；赵长煜等，2010）.以上模型均为运动学模型，岩石圈变形通过定义的速度场来实现，它们在预测盆地热流演化方面具有优势，但岩石圈流变性隐含在其均衡补偿模式中，无法反映岩石圈流变性及其变化对盆地构造热演化的影响，在预测盆地几何形态方面存在一定的缺陷（何丽娟，2000）.

而在拉张模型的另一类模型，动力学模型中，岩石圈流变性则是通过本构方程明确体现，能够实现热学与流变学的耦合，深入理解盆地动力学的基本过程（Fernàndez and Ranalli，1997）.根据其本构方程，动力学模型可分为平面应变模型（Braun and Beaumont，1989；Bassi，1991；Chéry et al.，1992；Huismans and Beaumont，2008；Van Avendonk et al.，2009；Huismans and Beaumont，2011）和平面应力模型（Sonder and England，1989；Bassi and Sabadini，1994；Polyansky，2002）.由于目前对盆地成因的动力学根本机制还不甚了解，以及模型初始扰动的给定、岩石圈流变学参数的不确定性等都在一定程度上阻碍了动力学模型的发展.

现有的拉张模型多采用固定的初始岩石圈、地壳厚度，软流圈多被赋予和岩石圈地幔相同的热参数，模型上边界的给定存在争议，且模型中较少地考虑沉积盖层的热效应.本文拟重点探讨拉张模型中初始地壳、岩石圈厚度、软流圈对流、模型上边界对构造热演化的影响，以及沉积盖层存在的情况下拉张盆地侧翼抬升的演化.文中的模拟均使用运动学模型，假定岩石圈为纯剪变形，利用有限元方法，在拉格朗日坐标系下，通过求解二维瞬态热传导方程来模拟盆地构造位移发生变化的同时，温度场和热流在时间、空间上的演化历史.

2 初始岩石圈、地壳厚度对构造热演化的影响

多数拉张模型都假定拉张开始前模型的初始岩石圈厚度为125km （McKenzie，1978；Jarvis and McKenzie，1980；He et al.，2002；Lin et al.，2002；He and Wang，2004；陈林等，2008；Song et al.，2010；Chen，2014），这也是大洋岩石圈厚度的渐进值（Parsons and Sclater，1977），该值在拉张开始前岩石圈未经历拉张或破坏的情况下往往是适用的.然而，在一些地区拉张开始前岩石圈厚度已经发生了较大的变化，比如中生代的华北克拉通.这时，需要根据实际需要选取合适的初始岩石圈厚度.

对于初始地壳厚度的选取，往往依赖于现今的Moho面深度.Dehler等（1997）在计算加拿大西缘夏洛特皇后盆地的构造热演化时使用的初始地壳厚度为34km，He和Wang（2004）在模拟济阳坳陷的热史时给定的初始地壳厚度是36km，Chen（2014）模拟南海白云凹陷时采用的初始地壳厚度为30km，这些模型预测的剖面现今Moho面深度都和实际观测值接近.

图1 初始地壳、岩石圈厚度对构造－热演化的影响Fig.1 Effects of initial crustal and lithospheric thicknesses on tectono－thermal modeling

表1 模型中使用的参数Table 1 Parameters and values used in modeling

为了了解初始地壳、岩石圈厚度对构造热演化的影响，本文建立了两组模型（图1）.图1a和b固定初始地壳厚度为35km，模型拉张10Ma，拉张系数为1.5，热沉降20Ma，计算不同初始岩石圈厚度下的构造沉降史及热流演化史.图1c和d是固定初始岩石圈厚度为110km，改变初始地壳厚度，计算构造沉降量及热流值.模拟中使用的参数见表1.结果显示，当初始地壳厚度恒定，初始岩石圈厚度从70km增加到120km时，计算得到的构造沉降量会随之减小，最大差异达660m以上；热流也随着初始岩石圈厚度的增加而降低，最大差异近25mW·m－2.当初始岩石圈厚度恒定，初始地壳厚度从31km增加到41km时，构造沉降量随之增大，最大差异超过900m；热流也与初始地壳厚度成正比，但最大差异不到6mW·m－2.由此看出，初始岩石圈、地壳厚度的选择，对构造热演化的影响非常大，尤其是对构造沉降量的影响尤为显著.

实际计算中，当拉张开始前的初始地壳、岩石圈厚度未知时，可给出多组初始厚度值.每给出一组初始地壳、岩石圈厚度值，通过模拟就可以计算得到一组现今热流、Moho面深度及热岩石圈厚度值，将计算值与观测值相对比，通过不断调整初始地壳、岩石圈厚度值，即可反推得到模拟剖面的最佳初始地壳、岩石圈厚度值.

3 软流圈对流对构造热演化的影响

大部分拉张模型是将软流圈和岩石圈地幔赋予相同的热导率，在计算中仅考虑热传导作用（McKenzie，1978；Keen and Dehler，1993；He et al.，2002；Lin et al.，2009；Chen，2014）.而 Lachenbruch（1978）首次指出了拉张背景下对流在热传递中的重要性.Steckler（1985）、Buck（1986）、van Wijk等（2008）都进一步指出被动裂谷作用会引起软流圈的小范围对流，其作为额外的热源会对沉降、抬升史及温度场产生重要影响.

为了讨论软流圈对流对构造热演化的影响，本文采用等效热传导率综合考虑软流圈热传导和热对流因素（范桃园和安美建，2009）.事实上，通过在热传导方程中提高热导率以等效模拟对流的热效应的方法已广泛被学者们采用（Morgan and Chen，1993；Chen，2001；Spinelli and Harris，2011；杨少华和石耀霖，2013；Schmeling and Marquart，2014）.在纯对流的情况下，相对于正常的热导率，提高的倍数往往为表征对流传热效率与单纯的热传导效率之比的努塞尔特数，Nu（Morgan and Chen，1993；Chen，2001；Cochran and Buck，2001；杨少华和石耀霖，2013）.在研究Nu随系统总体热状态变化特征中，实验和数值模拟都给出如下规律（Schubert et al.，2001；Turcotte and Schubert，2002）：

式中，c，β为常数，β值一般在0.3左右，Ra为瑞利数.研究表明，当Ra接近5×105～106时，c＝0.27，β＝0.3185（Schubert et al.，2001）.而在上地幔对流的情况下，Ra通常为106左右（Schubert et al.，2001；Jaupart and Mareschal，2011）.利用上式可以求得此时Nu接近22，也就是说，当认为软流圈地幔以纯对流的形式热传递的话，其等效热导率应近似为22倍的正常地幔热导率，约为70W·m－1K－1.本文计算了软流圈等效热导率从3.2W·m－1K－1（纯热传导，和岩石圈地幔热导率相同）到70W·m－1K－1（等效于纯热对流）变化时的构造沉降量、热流及温度场.结果显示，软流圈等效热导率对热沉降期的构造沉降量及热流影响很大，而在拉张期其影响相对微弱.考虑了软流圈的对流效应之后，相比于单纯的热传导模型，构造沉降量变小，而热流则增大.对流越强，这种影响越明显（图2a和b）.在温度场方面，整个过程中考虑软流圈对流的模型其下地壳至模型底部的温度明显比仅考虑热传导的模型温度高.在热沉降阶段，软流圈纯对流模型中温度场变化甚微，而在软流圈纯热传导模型中则能明显的观察到温度场的冷却（图2c）.

图2 软流圈等效热导率对构造热演化的影响（拉张期为10Ma，拉张系数为1.5，热沉降期为20Ma）（a）构造沉降史；（b）热流随时间的变化；（c）温度－深度曲线.Fig.2 Effects of asthenosphere equivalent thermal conductivity on tectono－thermal modeling（The extension lasts 10Ma with a stretching factor of 1.5，and then enters thermal subsidence until 30Ma）（a）Tectonic subsidence history；（b）Heat flow history；（c）Temperature－depth curve.

软流圈纯传导模型的热流随时间的降低速度要比实际情况快得多，而软流圈纯对流模型计算的沉降量要比实际小得多，二者都与实际偏差较大.事实上，纯传导的固体岩石圈与纯对流的流体软流圈之间存在一过渡层，即流变边界层（Sleep，2006；Artemieva，2011；图3），其间传导与对流共同作用来传递热量.在拉张模型中，同时考虑软流圈的热传导与对流更接近真实的情况，纯传导与纯对流应作为其下限与上限.在模拟中，得知地幔对流的努塞尔特数Nu或瑞利数Ra时，可相应的判断软流圈等效热导率的上限.流变边界层的厚度对岩石圈本身的结构特征并不敏感，而主要受软流圈黏性系数η控制，其厚度与lg（η）成正比（He，2014）.考虑到软流圈对流的瑞利数Ra与黏性系数η成反比，而由式（1）可知努塞尔特数Nu正比于Raβ，由此可推断地幔对流的努塞尔特数Nu越大，则流变边界层的厚度越小.这也意味着在拉张模型中，流变边界层厚度小的地方，流变边界层中等效热传导率的上限变大，而其下纯对流软流圈的等效热导率为Nu倍的正常地幔热导率，也相应大于流变边界层厚的地方的软流圈等效热导率.其对构造热演化的影响由上述模拟结果可以得知，流变边界层厚度小的地方，构造沉降量速率变慢，总构造沉降量变小，而热沉降过程中热衰减变缓，热流降低的总幅度变小.

4 模型上边界对构造热演化的影响

图3 温度（T）、地温梯度（G）随深度的变化（纯传导的岩石圈与纯对流的软流圈之间为流变边界层RBL，改自He，2014）Fig.3 Variations in temperature （T）and temperature gradient（G）with depth（The transition layer between the conductive lithosphere and the convective asthenosphere is the rheological boundary layer RBL，modified from He，2014）.

现有的拉张模型对其上边界采用两种处理方法.一种是整个过程一直固定上边界，其并不随计算出的沉降量而作调整（McKenzie，1978；Rowley and Sahagian，1986；Keen and Dehler，1993），是对模型的一种简化；另一种则是根据均衡补偿得到的构造沉降量对模型上边界在垂向上不断进行调整，更真实的还原盆地的演化过程（Wang et al.，1996；He and Wang，2004；Wangen et al.，2008）.为了探讨固定与移动边界对构造热演化的影响，我们建立了以下两组模型.一组模型的拉张时间为10Ma，之后裂后热松弛持续100Ma（图4a和b）.另一组模型的拉张时间为40Ma，裂后热沉降时间同为100Ma（图4c和d）.两组模型都在拉张系数分别为2、4和6时，对比固定与移动上边界模型的计算结果.结果显示，当拉张时间为10Ma时，固定与移动上边界对拉张期的热流及沉降量影响甚微；但当拉张时间增加为40Ma时，模型上边界的处理方式对拉张期的热流及初始沉降产生了明显的影响，尤其是在拉张系数较大的情况下，但这种影响小于热沉降期两种模型的构造－热演化差异.当拉张系数为6时，拉张结束时两种模型的热流差异为0.3mW·m－2，初始沉降量差异为30m（图4c和d）.在热沉降阶段，无论拉张时间的长短，固定与移动上边界模型的热流及构造沉降的差异随热沉降的持续均逐渐变大，移动边界模型的热流大于固定边界的，而构造沉降量小于固定边界模型求得的值.但上述两组拉张时间不同、热沉降时间相同的模型在热沉降结束时，由于上边界处理方式不同造成的热流、总沉降量及温度场的差异相差并不大，以下结果以拉张时间为40Ma的一组模拟为例.拉张系数为2时，固定与移动上边界模型的热流及构造沉降的最大差值分别为0.7mW·m－2及85m；当拉张系数为4时，二者的热流及构造沉降最大差异增大为1.7mW·m－2和162m；当拉张系数进一步增大到6时，两种模型的热流最大差异变为2.4mW·m－2，构造沉降量的最大差异为198m（图4c和d）.在温度场方面，从图4e—g可以看出，固定与移动上边界模型的温度结构有着明显的差别，尤其是在岩石圈部分，相同深度上，固定边界模型的温度比移动边界模型的温度高几十到上百摄氏度.当拉张系数分别为2、4和6时，两种模型在同一深度上温度的最大差异分别为95、161及185℃.

总体上看，模型上边界对构造热演化的影响随着拉张系数的增大而增强.相比于模型上边界的给定对热流及构造沉降量的影响，其对模型温度场有着更为显著的影响.从上述讨论也可以看到，当计算主要针对热流及构造沉降量时，在拉张系数小于2或是拉张时间小于10Ma且只关注于拉张期的热流及构造沉降变化时，两种模型的计算结果差异并不大，使用固定上边界模型的计算结果也是可接受的.但在拉张强度较大的地区，如中国南海的莺歌海盆地等（拉张系数大于4，He et al.，2002），使用简化的固定上边界模型将对计算带来较大的误差.

5 拉张盆地的侧翼隆起

裂谷的侧翼隆起常作为拉张盆地及张裂大陆边缘的边界，其是拉张盆地模拟的重要组成.学者们提出了多种机制来解释拉张盆地的侧翼隆起：（1）由拉张岩石圈向未拉张区域的侧向热传导（Alvarez et al.，1984；Buck et al.，1988；Leroy et al.，2008）；（2）由拉张引起的小规模地幔对流造成的裂谷之下向侧翼的热传递（Keen，1985；Buck，1986）；（3）裂谷侧翼之下相比于地壳减薄更强的岩石圈地幔减薄（Hellinger and Sclater，1983；White and McKenzie，1988；Dehler et al.，1997）；（4）拉张过程中通过向上挠曲对裂谷侧翼的力学支持（Braun and Beaumont，1989；Weissel and Karner，1989；Chéry et al.，1992；Kooi et al.，1992；van der Beek et al.，1994）.一般认为裂谷的侧翼隆起会在拉张结束后的几百万年内随着热衰减而消失（Weissel and Karner，1989；ten Brink and Stern，1992）.而另有观测数据显示一些裂谷侧翼的隆升可以长期维持，如非洲南部的大陡崖，有学者用具较大强度的岩石圈的挠曲均衡回弹来解释这一现象（Weissel and Karner，1989；ten Brink and Stern，1992）.但这些模型是在空盆载水的前提下讨论的，忽略了沉积物的热效应.本文重点讨论在载水及载沉积物两种情况下由侧向热传导造成的盆地边缘抬升的差异.

图4 固定与移动上边界模型的构造热演化结果对比图a和b是拉张时间为10Ma、热沉降时间为100Ma时的计算结果，图c—g是拉张时间为40Ma、热沉降时间为100Ma的结果，图a和c为两种模型计算的热流差值的绝对值随时间的变化，图b和d为构造沉降量差值的绝对值，图e—g为不同拉张系数下热沉降结束时的温度－深度曲线.Fig.4 Comparison of tectono－thermal evolution between fixed and moving upper boundary models（a）and（b）The extension lasts 10Ma.（c）—（g）The stretching period is 40Ma.In all figures，the thermal subsidence period lasts 100 Ma.（a）and（c）are absolute values of differences of the calculated heat flow between the two models.（b）and（d）are absolute values of differences of the calculated tectonic subsidence.（e），（f）and（g）are temperature－depth curves at the end of post－rift thermal subsidence stage with different stretching factors.

图5 载水和载沉积物模型的盆地侧翼隆起效应拉张时间为10Ma，拉张系数为a图所示，图b—d为不同模型的构造沉降演化图，图e和f分别为拉张前x＝－40km及x＝－30km处两点的抬升量随时间的变化.Fig.5 Rift flank uplift of water and sediment loading modelsThe extension lasts 10Ma.（a）Distribution of stretching factor.（b），（c）and（d）Graphs of tectonic subsidence history for different models.（e）and（f）Variations of uplift with time at two points x＝－40km and x＝－30km，respectively.

模拟采用二维非均匀拉张模型，拉张系数分布如图5a所示，拉张10Ma，之后为热沉降期.本文计算了载水、载热导率分别为2.2W·m－1K－1及1.6W·m－1K－1的沉积物三组模型的构造沉降／抬升史.结果显示，三组模型的盆地侧翼随着拉张逐渐抬升，但最大抬升量都并非出现在拉张刚结束时，而是在拉张结束后的6—10Ma内侧翼抬升到最高，之后慢慢消减（图5b—f）.上述计算中，载水模型的盆地侧翼最大抬升量为273m，小于载沉积物模型的侧翼最大抬升量；而热导率为1.6W·m－1K－1的沉积物模型比热导率为2.2W·m－1K－1的沉积物模型的侧翼最大抬升量略大，二者分别为497m及489m（图5e）.模拟结果同时显示，离拉张区域越远，盆地侧翼的抬升量越小（图5e和f）.值得注意的是，载水模型的盆地侧翼隆升会在200Ma左右消失，而考虑了沉积物的模型其侧翼隆升则会在200Ma左右稳定于某个数值，也就是说侧翼隆升并没有随时间消失.沉积物热导率为2.2W·m－1K－1的模型的盆地侧翼隆升最终维持在180m左右，而沉积物热导率为1.6W·m－1K－1的模型的侧翼隆升量则会稳定于约245m（图5e）.这可能是由于沉积物与盆地基底热导率差异产生的热折射效应造成的.由于沉积物的热导率较小，盆地基底的热导率大于沉积盖层的，隆起处的温度低于凹陷区的温度，凹陷处的热会向相邻隆起处转移，使得盆地凹陷处的热流降低，而隆起处的热流增加（熊亮萍和高维安，1982；熊亮萍和张菊明，1984；Wang et al.，1985）.在沉积盖层与基底的交界处，具有最大的水平热流值.经过长时间的热松弛后，温度场达到稳定.载水模型的温度场可以在一定时间的热冷却后恢复到初始值.而由于沉积层的热效应，载沉积物的模型盆地边缘的温度明显比载水模型的温度高，在不考虑剥蚀的情况下，前者的盆地侧翼抬升能维持在某一水平.沉积盖层与基底热导率差异越大，这种效应越明显，也因此产生更大的盆地侧翼抬升.

6 结论

（1）拉张模型中改变初始地壳、岩石圈厚度会对计算的构造沉降量及热流值产生很大影响.模拟中可根据需要选取多组初始地壳、岩石圈厚度值，将模拟结果与观测值进行对比以确定最佳初始值.

（2）软流圈对流对构造热演化的影响很大.软流圈纯对流及纯传导模型是拉张模型的两种极限，在流变边界层中应同时考虑热传导与对流.得知地幔对流的努塞尔特数Nu时，可相应的判断软流圈等效热导率的上限.

（3）模型上边界对构造热演化的影响随着拉张系数的增大而增强，固定与移动上边界对热沉降期的热流及构造沉降量的影响比拉张期的大.两种模型得到的温度场有着较为明显的差别.

（4）沉积盖层的热效应对裂谷盆地侧翼隆起影响显著.载沉积物模型的盆地边缘抬升可以长时间维持，而载水模型的侧翼抬升则会在拉张结束后逐渐消失.

Alvarez F，Virieux J，Le Pichon X.1984.Thermal consequences of lithosphere extension over continental margins：the initial stretching phase.Geophys.J.Int.，78（2）：389－411，doi：10.1111／j.1365－246X.1984.tb01956.x.

Artemieva I.2011.The Lithosphere：an Interdisciplinary Approach.New York：Cambridge Univ.Press.

Bassi G.1991.Factors controlling the style of continental rifting：insights from numerical modelling.EarthPlanet.Sci.Lett.，105（4）：430－452，doi：10.1016／0012－821X（91）90183－I.

Bassi G，Sabadini R.1994.The importance of subduction for the modern stress field in the Tyrrhenian area.Geophys.Res.Lett.，21（5）：329－332，doi：10.1029／93GL03543.

Baur F，Littke R，Wielens H，et al.2010.Basin modeling meets rift analysis—A numerical modeling study from the Jeanne d′Arc basin，offshore Newfoundland，Canada.Mar.Pet.Geol.，27（3）：585－599，doi：10.1016／j.marpetgeo.2009.06.003.

Braun J，Beaumont C.1989.A physical explanation of the relation between flank uplifts and the breakup unconformity at rifted continental margins.Geology，17（8）：760－764，doi：10.1130／0091－7613（1989）017〈0760：APEOTR〉2.3.CO；2.

Buck W R.1986.Small－scale convection induced by passive rifting：the cause for uplift of rift shoulders.EarthPlanet.Sci.Lett.，77（3－4）：362－372，doi：10.1016／0012－821X（86）90146－9.

Buck W R，Martinez F，Steckler M S，et al.1988.Thermal consequences of lithospheric extension：pure and simple.Tectonics，7（2）：213－234，doi：10.1029／TC007i002p00213.

Chen L.2014.Stretching factor estimation for the long－duration and multi－stage continental extensional tectonics：Application to the Baiyun Sag in the northern margin of the South China Sea.Tectonophysics，611（2014）：167－180，doi：10.1016／j.tecto.2013.11.026.

Chen L，Song H B，Dong C Z，et al.2008.2－D strain rate forward modeling.ChineseJ.Geophys.（in Chinese），51（6）：1819－1825，doi：10.3321／j.issn：0001－5733.2008.06.023.

Chen L，Zhang Z J，Song H B，et al.2013.Numerical modeling of extensional sedimentary basin formation with MATLAB：Application to the northern margin of the South China Sea.Comput.Geosci.，51：153－165，doi：10.1016／j.cageo.2012.07.014.

Chen Y J.2001.Thermal effects of gabbro accretion from a deeper second melt lens at the fast spreading East Pacific Rise.J.Geophys.Res.：SolidEarth，106（B5）：8581－8588，doi：10.1029／2000JB900420.

Chéry J，Lucazeau F，Daignières M，et al.1992.Large uplift of rift flanks：A genetic link with lithospheric rigidity？EarthPlanet.Sci.Lett.，112（1－4）：195－211，doi：10.1016／0012－821X（92）90016－O.

Clift P，Lin J，Barckhausen U.2002.Evidence of low flexural rigidity and low viscosity lower continental crust during continental break－up in the South China Sea.Mar.Pet.Geol.，19（8）：951－970，doi：10.1016／S0264－8172（02）00108－3.

Cochran J R，Buck W R.2001.Near－axis subsidence rates，hydrothermal circulation，and thermal structure of mid－ocean ridge crests.J.Geophys.Res.：SolidEarth，106（B9）：19233－19258，doi：10.1029／2001JB000379.

Davis M，Kusznir N.2004.Depth－dependent lithospheric stretching at rifted continental margins.／／ Proc.NSF Rifted Margins Theoretical Institute，1：192－137.

Dehler S A，Keen C E，Rohr K M M.1997.Tectonic and thermal evolution of Queen Charlotte Basin：lithospheric deformation and subsidence models.BasinRes.，9（3）：243－261，doi：10.1046／j.1365－2117.1997.00043.x.

Fan T Y，An M J.2009.Discussion on thickness of the Tarim Craton lithosphere and its mechanism.GeologicalReview（in Chinese），55（3）：329－334.

Fernàndez M，Ranalli G.1997.The role of rheology in extensional basin formation modelling.Tectonophysics，282（1－4）：129－145，doi：10.1016／S0040－1951（97）00216－3.

He L J.1999.Mutiple tectono－thermal modeling of Liaohe Bain in the Cenozoic.ChineseJ.Geophys.（in Chinese），42（1）：62－68.

He L J.2000.Advance in tectono－thermal modelling of sedimentary basins.Adv.EarthSci.（in Chinese），15（6）：661－665.

He L J.2014.Numerical modeling of convective erosion and peridotite－melt interaction in big mantle wedge：Implications for the destruction of the North China Craton.J.Geophys.Res.：SolidEarth，119（4）：3662－3677，doi：10.1002／2013JB010657.

He L J，Wang J Y.2004.Tectono－thermal modelling of sedimentary basins with episodic extension and inversion，a case history of the Jiyang Basin，North China.BasinRes.，16（4）：587－599，doi：10.1111／j.1365－2117.2004.00245.x.

He L J，Wang J Y.2007.Tectono－thermal modeling of sedimentary basins：review and outlook.ProgressinGeophys.（in Chinese），22（4）：1215－1219.

He L J，Xiong L P，Wang J Y.2002.Heat flow and thermal modeling of the Yinggehai Basin，South China Sea.Tectonophysics，351（3）：245－253，doi：10.1016／s0040－1951（02）00160－9.

Hellinger S J，Sclater J G.1983.Some comments on two－layer extensional models for the evolution of sedimentary basins.J.Geophys.Res.：SolidEarth，88（B10）：8251－8269，doi：10.1029／JB088iB10p08251.

Huismans R，Beaumont C.2011.Depth－dependent extension，twostage breakup and cratonic underplating at rifted margins.Nature，473（7345）：74－78，doi：10.1038／nature09988.

Huismans R S，Beaumont C.2008.Complex rifted continental margins explained by dynamical models of depth－dependent lithospheric extension.Geology，36（2）：163－166，doi：10.1130／G24231A.1.

Jarvis G T，McKenzie D P.1980.Sedimentary basin formation with finite extension rates.EarthPlanet.Sci.Lett.，48（1）：42－52，doi：10.1016／0012－821X（80）90168－5.

Jaupart C，Mareschal J C.2011.Heat Generation and Transport in the Earth.New York：Cambridge Univ.Press.

Keen C.1985.The dynamics of rifting：deformation of the lithosphere by active and passive driving forces.Geophys.J.Int.，80（1）：95－120，doi：10.1111／j.1365－246X.1985.tb05080.x.

Keen C E，Dehler S A.1993.Stretching and subsidence：Rifting of conjugate margins in the North Atlantic Region.Tectonics，12（5）：1209－1229，doi：10.1029／93tc00915.

Kooi H，Cloetingh S，Burrus J.1992.Lithospheric necking and regional isostasy at extensional basins 1.Subsidence and gravity modeling with an application to the Gulf of Lions Margin（SE France）.J.Geophys.Res.：SolidEarth，97（B12）：17553－17571，doi：10.1029／92jb01377.

Kusznir N，Ziegler P.1992.The mechanics of continental extension and sedimentary basin formation：a simple－shear／pure－shear flexural cantilever model.Tectonophysics，215（1）：117－131，doi：10.1016／0040－1951（92）90077－J.

Kusznir N，Karner G.2007.Continental lithospheric thinning and breakup in response to upwelling divergent mantle flow：application to the Woodlark，Newfoundland and Iberia margins.Geol.Soc.Spec.Publ.，282（1）：389－419，doi：10.1144／SP282.16.Lachenbruch A H.1978.Heat flow in the Basin and Range province and thermal effects of tectonic extension.PureAppl.Geophys.，117（1－2）：34－50，doi：10.1007／BF00879732.

Leroy M，Gueydan F，Dauteuil O.2008.Uplift and strength evolution of passive margins inferred from 2－D conductive modelling.Geophys.J.Int.，172（1）：464－476，doi：10.1111／j.1365－246X.2007.03566.x.

Lin C S，Zhang Y M，Li S T，et al.2002.Quantitative modelling of multiphase lithospheric stretching and deep thermal history of some Tertiary rift basins in eastern China.ActaGeol.Sin.－Engl.，76（3）：324－330，doi：10.1111／j.1755－6724.2002.tb00548.x.

Lin G，Peng M L，Zhao C B，et al.2009.Numerical analysis and simulation experiment of lithospheric thermal structures in the South China Sea and the Western Pacific.J.EarthSci.，20（1）：85－94，doi：10.1007／s12583－009－0009－0.

McKenzie D.1978.Some remarks on the development of sedimentary basins.EarthPlanet.Sci.Lett.，40（1）：25－32，doi：10.1016／0012－821X（78）90071－7.

Morgan J P，Chen Y J.1993.The genesis of oceanic crust：Magma injection，hydrothermal circulation，and crustal flow.J.Geophys.Res.：SolidEarth，98（B4）：6283－6297，doi：10.1029／92JB02650.

Negredo A，Fernandez M，Torné M，et al.1999.Numerical modeling of simultaneous extension and compression：The Valencia Trough（western Mediterranean）.Tectonics，18（2）：361－374，doi：10.1029／1998TC900026.

Parsons B，Sclater J G.1977.An analysis of the variation of ocean floor bathymetry and heat flow with age.J.Geophys.Res.，82（5）：803－827，doi：10.1029／JB082i005p00803.

Polyansky O P.2002.Dynamic causes for the opening of the Baikal Rift Zone：a numerical modelling approach.Tectonophysics，351（1－2）：91－117，doi：10.1016／S0040－1951（02）00127－0.

Rosenbaum G，Weinberg R F，Regenauer－Lieb K.2008.The geodynamics of lithospheric extension.Tectonophysics，458（1－4）：1－8，doi：10.1016／j.tecto.2008.07.016.

Rowley D B，Sahagian D.1986.Depth－dependent stretching：a different approach.Geology，14（1）：32－35，doi：10.1130／0091－7613（1986）14〈32：DSADA〉2.0.CO；2.

Royden L，Keen C.1980.Rifting process and thermal evolution of the continental margin of eastern Canada determined from subsidence curves.EarthPlanet.Sci.Lett.，51（2）：343－361，doi：10.1016／0012－821X（80）90216－2.

Schmeling H，Marquart G.2014.A scaling law for approximating porous hydrothermal convection by an equivalent thermal conductivity：theory and application to the cooling oceanic lithosphere.Geophys.J.Int.，197（2）：645－664，doi：10.1093／gji／ggu022.

Schubert G，Turcotte D L，Olson P.2001.Mantle Convection in the Earth and Planets.Cambridge：Cambridge Univ.Press.

Shillington D J，White N，Minshull T A，et al.2008.Cenozoic evolution of the eastern Black Sea：A test of depth－dependent stretching models.EarthPlanet.Sci.Lett.，265（3－4）：360－378，doi：10.1016／j.epsl.2007.10.033.

Sleep N H.2006.Mantle plumes from top to bottom.EarthSci.Rev.，77（4）：231－271，doi：10.1016／j.earscirev.2006.03.007.

Sonder L J，England P C.1989.Effects of a temperature－dependent rheology on large－scale continental extension.J.Geophys.Res.：SolidEarth，94（B6）：7603－7619，doi：10.1029／JB094iB06p07603.

Song H B，Chen L，Zhang J，et al.2010.A MATLAB program for 1Dstrain rate inversion.Comput.Geosci.，36（1）：16－23，doi：10.1016／j.cageo.2009.08.002.

Spinelli G A，Harris R N.2011.Effects of the legacy of axial cooling on partitioning of hydrothermal heat extraction from oceanic lithosphere.J.Geophys.Res.：SolidEarth，116（B9）：B09102，doi：10.1029／2011JB008248.

Steckler M S.1985.Uplift and extension at the Gulf of Suez：indications of induced mantle convection.Nature，317（6033）：135－139，doi：10.1038／317135a0.

ten Brink U，Stern T.1992.Rift flank uplifts and hinterland basins：comparison of the Transantarctic Mountains with the Great Escarpment of southern Africa.J.Geophys.Res.：SolidEarth，97（B1）：569－585，doi：10.1029／91JB02231.

Turcotte D L，Schubert G.2002.Geodynamics.New York：Cambridge Univ.Press.

Van Avendonk H J A，Lavier L L，Shillington D J，et al.2009.Extension of continental crust at the margin of the eastern Grand Banks，Newfoundland.Tectonophysics，468（1－4）：131－148，doi：10.1016／j.tecto.2008.05.030.

van der Beek P，Cloetingh S，Andriessen P.1994.Mechanisms of extensional basin formation and vertical motions at rift flanks：Constraints from tectonic modelling and fission－track thermochronology.EarthPlanet.Sci.Lett.，121（3－4）：417－433，doi：10.1016／0012－821X（94）90081－7.

van Wijk J，van Hunen J，Goes S.2008.Small－scale convection during continental rifting：Evidence from the Rio Grande rift.Geology，36（7）：575－578，doi：10.1130／G24691A.1.

Wang J Y，Wang J A，Xiong L P，et al.1985.Analysis of factors affecting heat flow density determination in the Liaohe Basin，North China.Tectonophysics，121（1）：63－78，doi：10.1016／0040－1951（85）90268－9.

Wang K L，He L J，Xiong L P，et al.1996.Thermal history model of a passive continental margin：Northern South China Sea.／／Wang J Y，Li N S，eds.Geothermics in China.Beijing：Seismological Press，275－292.

Wangen M，Fjeldskaar W，Faleide J I，et al.2008.Forward modeling of stretching episodes and paleo heat flow of the Voring margin，NE Atlantic.J.Geodyn.，45（2－3）：83－98，doi：10.1016／j.jog.2007.07.002.

Weissel J K，Karner G D.1989.Flexural uplift of rift flanks due to mechanical unloading of the lithosphere during extension.J.Geophys.Res.：SolidEarth，94（B10）：13919－13950，doi：10.1029／JB094iB10p13919.

Wernicke B.1985.Uniform－sense normal simple shear of the continental lithosphere.Can.J.EarthSci.，22（1）：108－125，doi：10.1139／e85－009.

White N，McKenzie D.1988.Formation of the“steer′s head”geometry of sedimentary basins by differential stretching of the crust and mantle.Geology，16（3）：250－253，doi：10.1130／0091－7613（1988）016〈0250：FOTSSH〉2.3.CO；2.

White N.1994.An inverse method for determining lithospheric strain rate variation on geological timescales.EarthPlanet.Sci.Lett.，122（3）：351－371，doi：10.1016／0012－821X（94）90008－6.

White N，Bellingham P.2002.A two－dimensional inverse model for extensional sedimentary basins 1.Theory.J.Geophys.Res.：SolidEarth，107（B10）：ETG 17－1－ETG 17－20，doi：10.1029／2001JB000173.

Xie X N，Müller R D，Li S T，et al.2006.Origin of anomalous subsidence along the Northern South China Sea margin and its relationship to dynamic topography.Mar.Pet.Geol.，23（7）：745－765，doi：10.1016／j.marpetgeo.2006.03.004.

Xiong L P，Gao W A.1982.Characteristics of geotherm in uplift and depression.ChineseJ.Geophys.（in Chinese），25（5）：448－456.

Xiong L P，Zhang J M.1984.Mathematical simulation of refraction and redistribution of heat flow.ChineseJ.Geol.（in Chinese），19（4）：445－454.

Yang S H，Shi Y L.2013.Numerical simulation of Ningwu Ice Cave preserved for thousands of years in Shanxi Province.Sci.Technol.Rev.（in Chinese），31（13）：23－28，doi：10.3981／j.issn.1000－7857.2013.13.002.

Zhao C Y，Song H B，Qian R Y，et al.2010.Research on tectonothermal evolution modeling method for superimoposed basin with the Jianghan Basin as an example.ChineseJ.Geophys.（in Chinese），53（1）：128－137，doi：10.3969／j.issn.0001－5733.2010.01.014.

附中文参考文献

陈林，宋海斌，董崇志等.2008.二维应变速率正演模拟.地球物理学报，51（6）：1819－1825，doi：10.3321／j.issn：0001－5733.2008.06.023.

范桃园，安美建.2009.现今塔里木克拉通岩石圈厚度分析及机制探讨.地质论评，55（3）：329－334.

何丽娟.1999.辽河盆地新生代多期构造热演化模拟.地球物理学报，42（1）：62－68.

何丽娟.2000.沉积盆地构造热演化模拟的研究进展.地球科学进展，15（6）：661－665.

何丽娟，汪集暘.2007.沉积盆地构造热演化研究进展：回顾与展望.地球物理学进展，22（4）：1215－1219.

熊亮萍，高维安.1982.隆起与拗陷地区地温场的特点.地球物理学报，25（5）：448－456.

熊亮萍，张菊明.1984.热流的折射和再分配的数学模拟.地质科学，19（4）：445－454.

杨少华，石耀霖.2013.山西省宁武县万年冰洞持续存在机制数值模拟.科技导报，31（13）：23－28，doi：10.3981／j.issn.1000－7857.2013.13.002.

赵长煜，宋海斌，钱荣毅等.2010.叠合盆地构造热演化模拟方法研究——以江汉盆地为例.地球物理学报，53（1）：128－137，doi：10.3969／j.issn.0001－5733.2010.01.014.