顾及下地壳拖曳力的川滇地区地表形变有限元数值模拟

摘要

在川滇地区,GPS观测资料揭示了其运动方向为东向—东南向,呈现顺时针旋转的趋势。已有研究表明该区域的下地壳是柔性的,且运动速度快于脆性的上地壳并对其具有拖曳作用,这是造成该区域独特运动特征的主要因素。因此在数值模拟区域形变场时,仅考虑弹性作用是不够的,还需要考虑柔性下地壳的拖曳力作用。拖曳力的分布通常采用人工分区和试错法得到,主观性较强。根据川滇地区的构造活动特征建立二维有限元弹性平面模型,以GPS观测数据为约束,利用未施加拖曳力时的计算结果与实际观测结果之间的位移残差ΔUx,ΔUy与各点位坐标x,y的拟合关系式反映研究区域柔性下地壳对上地壳的拖曳力作用。计算结果表明:拖曳力主要集中在川滇菱形块体内,且块体北部的值最大,高达1.0e8 N,并以南向为主;块体南部拖曳力逐渐减小,约为4.0e7 N,方向上逐渐向东西两侧分流。拖曳力以等效体力的形式加载至模型中,可得到与实测形变场吻合度较高的结果。该方法可以为川滇地区长期地表形变动力学机制的研究提供一种物理上更加合理、方法上更加简洁的新思路。