论文部分内容阅读
该文把大陆岩石层视为一形变服从牛顿或者幂指数流变律的连续介质,用有限差分的数值方法了其大时间、大空间尺度的形变.模型假设大陆岩石层的运动为上伏在一非粘滞性或者一粘滞性相对大陆岩石层很小的基底之上流动,并假设流场的水平速度的垂直分一可忽略,在这种假设下数值模拟首先在一简单的正方形界框架下给出了一组一般性的结果,通过对比这组结果与Engiand and McMkenzieⅰ1982,1983ⅱ的结果,验证了该文数值方法及程序的正确性.在上述简单模型的基础上,进一步给出了在一与青藏高原及其外围地区的实际地质框架比较相信的边界模型中的数值结果,该边界模型把现代板块构造中比较明显的欧亚古陆块与东亚三解地带的分界线作为一强边界限制条件,此时边界形态为一梯形区域.在该数值模拟中,研究人员用一边界简单正变换和一网格的密反变换把大陆形变在上梯形区域的求解变换到对一正方形区域上求解.首先为了验证该变换的有效性和正确性,文中对上简单正方形边界形态下的结果进行了恢复,然后在上梯形边界形态的基础上,给出了四种边界模型,并给出了他数值模拟结果.在上述边界形态和边界模型条件下,大陆岩石的流动受应力指数n及Argand数Ar控制.Ar数表示地壳的厚度差所引起的应力与维持形变以周围环境应变率进行的应力的比.由数值结果可见,在一挤压嵌入边界条件下,大陆岩石层视为一粘性薄层的变形会产生地壳厚度,速度场及应变率场的分解.与England的结果不同之处也是更符合实际观测结果之处在于在梯形边界形态下,地壳厚度及速度场会出现更大程度的向东的扩展而且应变革场也会出现东西向的扩展.同时数值结果还显示一旦高原形成,由地壳厚度差所引起的浮力将限制地地壳的进一步加厚,岩石层的形变则主要表现为向物质流入边界近前方区域的外围的扩展.该文致力于提供一定量方法及其数值模拟结果,给出把岩石层作为连续介质薄层时其形变的速度场、应变率场和地壳厚度分布及其一般性的特征,以后将可以用这一定量方法对青藏高原地区的岩石形变的其它特征应力场的分布、地表的隆升、断层的展布等进行模拟.