论文部分内容阅读
大陆动力学是当今固体地球科学研究的主导方向之一,其核心问题是大陆构造变形及其动力学机制.造山带是大陆变形最强烈、地表形态最明显的构造,它所涉及的地表地质结构、地球物理场特征、深部作用过程以及地表作用,一直是人们关注的热点.造山动力学研究对推动大陆构造变形及其动力学机制研究具有不可替代的作用.天山造山带远离板块边界,是现今世界上最为活跃的陆内造山带,被公认为研究陆内造山的天然实验场.刘启元(2004)提出了天山造山带动力学的10个关键科学问题.其中,塔里木板块的水平挤压与地幔对流作用对天山造山带的影响作用,以及二者之间相互关系是最核心的问题.现今天山造山带的活动不属于原始的碰撞造山过程,而是在被剥蚀夷平的古老造山带上的“复活”过程.同时,深部地球物理探测显示,地幔对流可能对天山造山带的复活起到了相当重要的作用.因此,确定中国境内天山造山带下方地幔中是否存在对流,是天山陆内造山动力学研究需要首先解答的关键问题.随之,如果存在地幔对流,那么对流的形态如何?地幔对流对造山过程产生了什么影响?塔里木板块的水平挤压与地幔对流构成了怎样的相互关系?塔里木板块的水平挤压作用是否足以造就现今天山的变形?古老的天山造山带内部结构与物性在造山复活过程中起到了怎样的作用?这些问题都有待于在资料解释和定量分析的基础上进行深入探讨.针对上述问题,本文首先确立了大陆构造变形及动力学定量研究的数学模型.采用热-流体连续介质的动量守恒方程、能量守恒方程和质量守恒方程,并引入了流体热力学的状态方程描述地幔对流和造山动力学过程.实现了求解这些方程的“ALE FEM+MIC”方法,即“任意拉格朗日-欧拉(ALE)描述的有限单元方法(FEM)”和“网格-粒子(MIC)”技术相结合的数值求解新技术.在热-蠕变耦合的方程中,基本未知量是速度场和温度场.未知量通过ALE描述的有限单元方法求解.同时,方法中引入了“网格-粒子”技术,即在有限单元内部设置若干代表物质的“粒子”,通过粒子的移动跟踪物质的运动和变形.有限单元方法求解动量方程和连续性方程时采用了速度场和压力场等阶插值的压力场稳定化Petrov-Galerkin(PSPG)方法,求解能量方程时采用了流线迎风Petrov-Galerkin(SUPG)方法.ALE描述通过定义合理的网格移动方式跟踪计算区域边界的变化.网格-粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格-粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.ALE FEM+MIC方法使得包含边界变化的极其强烈的变形流动问题和强烈不连续边界条件问题均可以正常描述.算例表明计算程序具有较好的计算精度和数值稳定性.