地震势能剖面

地质体在应力场的作用下发生弹性形变,从而具有应变位能（势能）。应用应变位能公式,可以将位或位移转化为势能,用势能剖面代替常规地震剖面。由于势能转换能拓宽地震信号的有效频谱,具有高通滤波性质,同时又不损害低频信号,所以,同常规地震剖面相比,势能剖面在保持信噪比的前提下,有更高的纵、横向分辨率。通过势能转换的理论分析及合成数据和实际数据的势能转换处理,本文系统地研究了势能转化问题和实现方法,从而证明势能转换确实是行之有效的高分辨率地震处理方法,可用于实际地震资料的处理。     

二维线性流变体波的有限元模拟

用加里津方法推导了二维线性流变体内波传播的有限元方程，方程内包含了边界项，使得处理人工边界影响有较好的结果．同时考虑了其它一些因素，例如，吸收介质、爆炸源、自激自收等．实例表明，作者在仿真性方面的努力是成功的，计算的垂直地震剖面理论地震图质量比较好，结果也显示了地表低速层对地表水平记录剖面的干扰．     

弹性波克希霍夫积分偏移法

P、S波同时偏移法是为了适应复杂构造地区勘探情况而开始研究的一种偏移方法。在前人工作的基础上,应用弹性波克希霍夫积分偏移的基本原理,改进了纵波（PP）及转换波（PS）的同时克希霍夫积分偏移方法,给出了一种近于实用的计算机流程和程序;并用其分别对射线理论合成反射波资料、有限元模拟地震波传播资料进行了检验,证明此种方法是有效的。在解决复杂构造地区的地震波成象方面,该方法具有很大的潜力。     

一种最短路径射线追踪的快速算法

为提高最短路径射线追踪的精度,需要增加模型的剖分网格和离散节点,并增加子波传播方向,或者采用其他方法改善计算结果,这些处理会带来大量的额外计算.本文的快速算法改进了波前点的管理和子波传播的计算这两项耗时的工作,较大幅度地提高了传统算法的效率.在波前点的管理上,采用按时间步划分区间的方法,实现了波前点的桶排序管理,其效率高于传统方法中常用的堆排序算法. 在子波传播的计算上,利用斯奈尔定律,同时参考来自邻近节点的波的走时,来限定当前子波传播的有效区域,排除大量不需要计算的子波传播方向. 模型实算表明,本文快速算法的计算速度是传统方法的几倍至十多倍.     

伪随机编码源激发下的时域电磁信号合成

将伪随机编码技术引入到人工源电磁法后,可以通过加大发射功率以及应用后续的相关处理技术来达到压制噪声、加大探测深度及提高分辨率的目的,因此引起了越来越多学者的关注及研究,但大多数研究集中在资料处理的相关技术上,对模拟电磁信号关注较少.然而,资料处理工作大多是从电磁信号出发的,模拟伪随机编码源激发下的电磁信号不但可以为资料处理环节提供理论数据,而且可以为检测资料处理的效果提供中间结果,因此,模拟电磁信号工作必不可少.本文根据获得接收信号的物理过程来实现伪随机编码源激发下的电磁信号合成.首先用解析公式获得特定地电结构的大地频率域响应,然后通过余弦变换得到时间域阶跃响应,接下来用阶跃响应的时间导数得到大地脉冲响应,通过将大地脉冲响应与伪随机编码源的褶积得到理想接收信号,最后,用低通滤波器来模拟发射设备和接收设备的频带限制,将之和噪声一起加到理想接收信号上,最终模拟出仿真的合成信号.通过和野外实际接收信号对比发现本文合成信号仿真度较高,可以服务于后续的数据处理环节.     

雷达波有限元仿真模拟

为了突出雷达波自身的动力学特点,本文首先用Ｇａｌｅｒｋｉｎ’ｓ方法推导了含衰减项的雷达波有限元方程,用有限元方法实现了管状体、弯曲界面等复杂形体的雷达波场的仿真模拟．通过对比同一模型合衰减和不含衰减两种情况的波场,明显地展示了考虑雷达波自身特点后的波场仿真性较好．     

复杂介质有限元法2.5维可控源音频大地电磁法数值模拟

利用有限元2.5维可控源音频大地电磁法（简称CSAMT）数值模拟方法,对100Ωm均匀半空间介质中有限长度的电偶极源产生的电场、磁场及视电阻率、相位特征进行了数值模拟,研究了场的空间变化规律. 在一个象限中,场的特征存在双叶现象,当收发距大于4个趋肤深度时,电阻率较接近介质真实的电阻率,这些结果为观测系统和收发距的选择提供了依据. 波数域场的特征表明,低波数对源的贡献占较大的比例. 有限元法2.5维CSAMT数值模拟的优势在于能较准确地获得复杂介质结构的波场特征. 本文结合直立异常体、倾斜异常体及断陷模型对CSAMT电阻率、相位剖面特征及频率曲线特征的可靠性进行了研究. 数值模拟结果直观地给出了异常体的剖面异常形态. 通过对比研究异常体的剖面异常形态和半空间场的特征进一步说明本文方法和软件在模拟复杂介质结构场特征时是可靠的. 这为认识观测数据,指导反演解释提供了较好的依据.     

利用地质规则块体建模方法的频率域有限元弹性波速度反演

在频率域弹性波有限元正演方程的基础上,依据匹配函数（也就是观测数据和正演数据残差的二次范数）最小的准则,用矩阵压缩存储与LU分解技术来存储和求解频率域正演方程中的大型稀疏复系数矩阵、用可调阻尼因子的Levenberg Marquard方法求解反演方程组,直接求取地下介质的弹性波速度,导出了频率域弹性波有限元最小二乘反演算法. 为了利用地下地质体的分布规律,减少反演所求的未知数个数,本文又提出了规则地质块体建模方法引入到反演中来. 经数值模型验证,在噪声干扰很大（噪声达到50髎）或初始模型与真实模型相差很大的情况下,反演也能取得很满意的效果,证明本方法具有很好的抗噪性与“强壮性”.     

衰减雷达波有限元偏移

高频雷达波在地球介质中有较强的衰减,反演中不可忽略．为此文中首先给出了含衰减项的雷达波的有限元方程及其偏移理论．用有限差分法或有限元法可正演合成雷达波资料,加入一定的扰动后用含衰减项的雷达波有限元方程做偏移,实例结果表明,考虑衰减项的偏移结果能使界面更好地归位,这为提高探地雷达地质解释的分辨率提供了可能性,为逐渐地实现符合雷达波自身动力学特点的处理系统奠定了基础．     

CSAMT法和高密度电法探测地下水资源

我国乃至世界目前都面临着严重的水资源的短缺，利用新技术、新方法合理勘探开发并利用水资源，是经济腾飞的关键。本文概述了CSAMT方法和近几年来发展起来的高密度电阻率找水方法，列举了两种方法应用的的实例。实践表明这两种方法是水资源勘查的有效手段。