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

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

地壳对海洋潮汐的响应

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

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

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

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

中国大陆科学钻探先导孔零偏VSP资料解释

用六级三分量检波器在中国大陆科学钻探先导孔中实施了零偏VSP测量.数据处理结果表明,中国大陆科学钻探孔区超高压变质岩石的地震波速度主要介于4500～7000m/s之间,显著高于一般的沉积岩地区,而且随深度变化不明显.声波测井速度系统地稍低于VSP层速度,可能是由于井壁处岩石的完整性受到破坏而造成的.地震波速度与岩石密度和岩性存在明确的对应关系,榴辉岩的密度和地震波速度均显著高于片麻岩类岩石;由榴辉岩退变生成的斜长角闪岩类岩石,其密度和地震波速度均呈现出较大的变化,主要与其退变质程度有关;超基性岩中的裂隙系统导致其密度和地震波速度大幅度下降.由于榴辉岩与其他岩石类型之间存在较大的波阻抗差异,因此用零偏VSP资料标定该区地震波的地质层位是有效的.关于地震波反射的原因,通过综合研究地震波(包括反射纵波、上行转换横波、井筒波)的特征、岩石速度和密度分布以及井径变化,认为主要是岩性分界面、韧性剪切带和断裂(带),但还有一些因素尚待进一步研究.     

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

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

青藏块体东北缘地壳水平运动状态

应用青藏块体东北缘1999～2003年多期GPS观测资料，计算了不同时段GPS点水平运动速率。通过分析发现：甘青块体可分为东部块体和西部块体，东、西部块体的运动状态存在明显的差异；受2001年11月14日昆仑Ms8．1地震的影响，震后地壳运动状态发生了明显的改变。     

高准直双极分子外流速度的“哈勃”律

在高准许直双极分子外流中，气体速度和离中心恒星距离的“哈勃”关系，可能是分子外流物理性质中最关键的特征。高速喷流驱动弓激波的方式能够解释这种速度结构。本建立起一组喷射角可变的微分方程组。用数值方法得出分子外流的速度结构，这些结构非常类似于观测结果。