激波在磁盔－电流片位场中的传播：（Ⅰ）激波与磁盔间的相互作用

利用新近获得的子午面磁盔－电流片背景太阳风稳态解，对激波从盔底沿电流片方向往外传播时与磁盔间的相互作用进行了数值模拟研究，重要新结果是：１．磁盔的存在使受扰介质速度跃变中央出现下凹，随着激波传出磁盔区并沿电流片方向传播，速度下凹逐渐减弱以致消失；２．激波将磁盔拉长并把盔顶的环形（垂直赤道面）磁场带到行星际空间，成为行星际磁场南向分量的来源之一；３．５个太阳半径（Ｒ⊙）内的磁盔部分将出现精细结构，沿盔外边界形成两条高速带，以及马蹄形密度（亮）环形结构等．这些结果表明，太阳附近高速等离子体与磁盔间存在重要的动力学相互作用过程，对行星际空间的太阳风三维结构有重要影响．     

层状各向同性多孔介质中井间地震的数值计算

依据Ｂｉｏｔ理论，提出计算理论地震图的相应传播矩阵算法，并针对井间地震震源和检波器的不同配置进行数值计算，结果表明，本文算法高频可达１０００Ｈ２以上．文中算例说明利用震源和检波器相对于储油层的不同位置，可从反射震相记录波形上的差异，判定储油层存在与否．     

裂隙各向异性介质波场VSP多分量记录的数值模拟

从各向异性介质中最为普遍的裂隙各向异性介质出发，推导了同零井源距ＶＳＰ三分量观测方式相一致的一维三分量波动方程．为了更符合实际观测，引入了坐标旋转，可以模拟井中任意方位的多层各向异性介质的波场响应．用有限元方法进行了波场数值模拟，重点突出了获得凡Ｒ_３×３矩阵剖面的方法，同时对Ｒ_３×３，矩阵剖面进行了一般性的分析解释，事实证明，Ｒ_３×３矩阵剖面含有介质丰富的各向异性信息．     

成层介质中交流电测井响应

本文给出了纵向成层、径向不均匀介质中轴对称电磁波的反射阵和透射阵的递推公式．用这些阵可以计算交流电测并响应，包括感应测并和深探测电磁波传播测井（ＤＰＴ）．同时还给出了幅度基函数和斜度基函数，它们分别＂支撑＂起结点上函数的幅度和斜度．     

粘弹性场地地形对地震动影响分析的显式有限元-有限差分方法

推导出了分析二维粘及弹性场地地形对地震动影响的显式有限元-有限差分方法．这一方法中，首先利用人工边界及有限元离散方法，给出问题分析的有限元离散网格计算力学模型，并利用一种类似于差分方法的有限元方法，建立局部网格节点的动力方程，而后利用笔者提出的有阻尼体系动力方程求解的显式差分格式，及推广的多次透射边界公式，给出网格节点运动量计算的时域显式逐步积分公式．利用计算机程序实现这一方法的计算具有所需计算机内存量小及计算时间量小的优势，而且，这一方法适用于任意地形情况，具有较高的计算精度及较好的计算稳定性．      

利用地震层析成象技术确定横向非均匀介质中的速度分布和界面位置

本文提出一种横向非均匀介质中的地震层析成象方法，由地震走时资料确定介质的速度分布和界面位置。计算中采用了参数分离、变阻尼、奇异值分解等技术，方程的解较稳定，且可得到解的分办。除反射波外，该方法还可同时利用分析波的资料。因此，适宜于折射及宽角度反射地震测深剖面的资料处理。数字试验结果表明，在模型内部，该方法能得到可靠的反演结果。在边界附近效果较差，受初始值影响较大。     

横向非均匀介质中弹性地震波传播的数值模拟

应用混合变量弹性动力学方程和线性常微分方程组的矩阵指数解法，将层状介质中广泛应用的弹性波传播矩阵解法推广至横向非均匀介质，给出了一种可计算复杂地质体中弹性波传播的广义传播矩阵数值解法。该方法可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度。     

粘弹性组合介质中应力波的数值模型

对分块均匀和各向同性的粘弹性介质,本文提出一种适于应力波数值计算的守恒型有限元法,旨在用以研究各类地质体（包括油、气藏）对于地震波的响应。     

地震波在湿颗粒介质中的传播

本文利用拉格朗日运动方程建立了湿颗粒介质中的波场方程组,由此研究了平面简谐波在其中的传播。文中还导出了湿颗粒介质中的支承拉梅系数和耦合弹性系数的表示式。     

考虑土体固液二相性质的地下管线地震反应研究

本文将管线周围的土体视为固-液二相介质，采用考虑地震过程中土体内孔隙水压力的增长和消散的有效应力分析方法，并以非线性的土体本构模型模拟土体动力特性的改变以及在孔隙水压作用下的变形，以求接近土体的真实地震反应。研究表明随着土体软化程度的不同，受到土体约束的管线应力亦发生相应变化。文中详细研究了地震动峰值、管线直径、行波波速和入射角度等对管线应力的影响。