青藏高原隆升的非线性动态有限元仿真研究

根据青藏高原的地质特征建立分析模型，采用3维动态有限元方法，在计算仿真板块速度场的基础上，计算在青藏高原的隆升过程中该地区地壳岩石的等效应力和位移随时间的变化，计算仿真得到的速度场与1998年GPS观测的速度场吻合良好；与过去一贯的假设相反，计算结果反映出地壳应力场不是静态的，而是此起彼伏，不断变化的，应力值最大且变化最剧烈的地区在克什米尔地区、鄂尔多斯地区和鲜水河－小江断裂带，与地震多发区域吻合。     

对青藏东北缘现今块体划分、运动及变形的初步研究

利用2维非连续变形分析方法（DDA），以位移代替围压作为边界约束力，研究青藏东北缘现今块体划分及其运动变形。根据该地区地质构造及地震活动，以GPS点测量位移作为模拟结果约束点，得出了较合理的块体划分模型和随时间演化的主应变分布图，并把应变高值区与近几年来发生的5级以上地震作对比，得出了研究区内地震危险性可能较大的区域。另外，对模拟的甘青块体与阿拉善块体的边缘带断裂左旋运动做了大概计算。     

地壳对海洋潮汐的响应

应用三维动态有限元方法研究了中国北部地区的地壳对邻近的渤海与黄海海平面变化的响应。虽然此应力场过于微弱不足以引发地震，但发现应力集中的位置及应力场变化较大的位置恰好与某些现代地震的震中一致。这一结果表明研究地壳对广泛分布的载荷的响应对研究区域地震构造是有帮助的。     

昆仑山8.1级地震前中国大陆的构造应变背景

利用“网络工程”1998～2001年累积的1181个测站的GPS重复观测资料，采用双三次样条函数模型建立中国大陆水平运动模型速度场，用大地坐标在椭球面上计算各类应变场，详细分析了2001年昆仑山8．1级地震前中国大陆水平构造应变场空间分布特征。各类构造应变场的最高值都出现在喜马拉雅构造带与昆仑山地块内(地震断裂带南侧)，鲜水河—安宁河断裂带次之。分析表明，昆仑山8．1级地震正好发生在张性面膨胀应变率的高值区，第一、第二和最大剪应变率高值区边缘的突变区和最大、最小主应变率的高值区。     

地震前后垂直形变场动态演化的量化指标

提出一种表达垂直形变场动态演化过程的量化指标——区域应变率、应变集中度。在此基础上，对南北地震带各水准监测区近30年的垂直形变资料进行了实际计算．并结合具体震例进行了对比研究。结果表明：量化指标在一定程度上能够反映地震前后形变场的动态演化过程和地壳运动状态，对地震的中长期预报有一定的积极意义。     

描述地壳垂直运动过程的一种函数形式

本文在已有数据处理方法的基础上，利用近代数值逼近理论，给出了从时空域角度描述地壳垂直运动过程的一种具体的函数解析形式。最后给出了一个实际算例。     

定水头注水引起的含水层水平运动和应变

基于含水层固体颗粒与孔隙水不可压缩的假设 ,本文导出了单井注水情况下泰斯承压含水层水平运动速度与水头之间的基本关系式。然后利用注水井壁处的应力、应变边界条件 ,进一步导出了单井定水头注水引起的泰斯承压含水层水平运动速度、位移和应变解析表达式。该水平位移与应变由两部分组成 :一部分为由注水压力本身引起的经典弹性力学解项 ,它仅随半径而变化 ,与注水时间无关 ;另一部分为由地下水头变化引起的水动力学位移和应变解项。其中 ,含水层水动力学水平位移随时间加长呈指数增长特征 ,水动力学径向应变则表现为近井处拉张、远井处挤压的分区特征 ,且近井拉张区随时间加长逐渐向外扩展。单井注水含水层水动力学水平位移、应变解的导出 ,完善和发展了单孔内压经典弹性平面力学问题解     

富氏谱分析法模拟地壳垂直运动的地学依据及其策应力模型

应用地壳波浪与镶嵌构造学说对富氏谱分析法提取地壳垂直形变信息的科学性做了地质学意义上的阐释 ,并提出了根据多期形变资料提取特定波段上构造策应力的数学模型