煤层上覆地层含水不均匀性电法探测的可能性

煤田上覆岩层中含水直接对煤层的开采构成威胁，本文先从理论上阐述了电阻率勘探可以探测煤层上覆地层含水不均匀性的可能性，然后用简单的断层型薄层含水模型的CSAMT正演模拟结果展示出：地下900m深处，尽管只有20m厚的的含水岩层，但在电阻率剖面上该薄层的反映具有非常明显的特征，从理论上和数值模拟实例说明，一定深度内一定厚度的薄层含水性的电阻率法探测是可能的。     

应迅速开展对我国主要自然灾害的成因与规律的研究

一、引言 我国是一个多自然灾害的国家,在世界上也是一个气候脆弱区之一。每年国家由于自然灾害遭受的损失,少则几亿,多则几十亿,甚至几百亿。由于这些灾害的成因与规律目前还不是很清楚,因此,对这些灾害进行可靠的预测还比较困难。如果人类掌握了这些灾害的规律,并认识到它们的物理成因,这样就完全可以预测这些自然灾害的产生,从而最大限度地减轻它们造成的损失。     

各向异性弹性波有限元叠前逆时偏移

利用有限元法和最小走时射线追踪的界面点法，实现了各向异性弹性波的叠前逆时偏移.理论模型资料的偏移结果清晰准确，证实了该系统的可靠性.通过对各向异性模型资料的各向同性偏移处理，发现常规偏移剖面存在较大误差，地质体的垂向深度和横向位置与实际模型有偏离.     

二维大地电磁数据的整体反演

从麦克斯韦第一、二方程出发，经过变量替换，将二维电磁问题线性化，得到不同测点、不同浓度下关于导电率微小变化的线性方程组，用Gauss-Seidel迭代法求解。实现了整体对二维大地电磁数据的反演，由于对各个测点下的数据同时反演，克服了以往用一维方法来反演二维问题时固有的缺点，讨论了不同的初始模型对反演结果的影响，提出了一种简便实用的初始值的给定方法，数值模拟结果表明该方法有一定的抗噪能力。     

金沙江龙蟠右岸变形体的地震学研究

拟建中的虎跳峡水电站水库蓄水后距虎跳峡上游20 km处的龙蟠右岸变形体是否会突然失稳崩滑为许多人关注.本文在折射波勘探方法获得龙蟠变形体底界形态、埋深、规模的基础上，进一步利用折射波勘探已经采集的炮集记录识别的反射波、面波、折射波等来研究变形体的内部结构和弹性及非弹性参数，获得变形体内部存在纵向横向结构的依据以及面波和折射波衰减系数，进而估算了变形体介质横波和纵波的非弹性性质参数与弹性性质参数的比值分别为0000184和0000144.由此判断水库蓄水后因变形体介质大量孔隙充水，造成孔隙压增加和剪切强度降低，进而导致变形体可能发生稳定非弹性变形或滑动，但稳定变形或滑动是否会演化成突然的非稳定滑动尚缺乏资料的支持.     

南水北调西线千米深长隧洞围岩构造地球物理勘探

2004年对南水北调引水隧洞西线一期工程中玛柯河－贾曲段20 km深埋长隧洞进行了地质结构面、岩性特征、赋水性等地球物理勘探、工程地质勘查，并进行综合解释以查清勘探区影响引水隧洞工程建设的不良地质体和不良地质构造.本文介绍了CSAMT法隧洞勘探的特点、资料采集的过程及应特别注意的复杂地形资料处理技术，结合其他补充物探手段及地面工程地质勘察结果对20 km深埋长隧洞围岩介质的结构特征进行了解释，指出了可能影响工程的断层、破碎带及异常区     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.     

有限元法2.5维CSAMT数值模拟

目前国内大多数CSAMT资料处理都用一维反演方法,即假定地球介质是水平层状地电结构,这对于横向变化较大的介质结构的反演则是不真实的．为了合理的模拟可控信号源的三维特征应该使用三维有限元方法．受计算机内存和计算速度的限制,三维方法的实用化受到约束．在许多情况下,地电结构沿走向变化很小,只沿倾向发生变化．这种地电结构是二维的,而人工源是三维的,因此CSAMT资料的观测可用2．5维有限元方法进行数值模拟．本文从麦克斯韦方程组出发,建立了2．5维有限元CSAMT数值模拟方法,其核心是把地电参数变化小的走向方向转化成波数域,用一系列波数模拟三维源的特征．用一个横向均匀的三层地电结构模型展示了2．5维数值模拟的特征,并和一维模拟结果进行了比较,证实了有限元方法2．5维CSAMT数值模拟的可靠性．在此基础上对一个地电结构已知的实际模型进行了2．5维数值模拟并和该剖面的野外实测剖面数据进行了比较,进一步有效地说明了本文介绍的2．5维CSAMT数值模拟方法是仿真和可靠的,为在此基础上的2．5维反演打下了良好的基础．     

声波方程频率域有限元参数反演

推导出频率域有限元声波正演方程,为了消除边界反射,将Clayton-Engquist旁轴波动方程吸收边界条件引入频率域,并对有限元刚度矩阵和质量矩阵进行压缩存储,利用广义共轭梯度法求解有限元方程获得正演解.在此基础上,推导出在某一频率下波场数据残差δU与单元物性参数修改量δλ之间关系的Jacobi矩阵,反演方法允许利用地面二维炮集全波场资料与给出初始模型参数的正演值的差值δU,迭代求得δλ.由于计算机内存的限制,方法计算不允许有过多数目的未知数个数,因此还提出了对同一介质物性单元的Jacobi矩阵元素进行压缩组装的措施,从而使反演的未知量个数减少,结合采用共轭梯度迭代法,使得只需利用有效波频段的少数一些频率即可进行迭代反演.正演和反演理论模型的数值模拟结果表明方法是有效的.     

改进Moser法射线追踪

地震波场正演模拟、层析成像、偏移成像经常需要用到射线追踪,本文在基于图形理论的Ｍｏｓｅｒ法基础上,改进其网格节点的划分、追踪时路径的选取,并增加直射线追踪以消除在速度变化不大时射线受节点布置影响出现不合理的折曲现象．