用视模型空间对比进行地球物理反演

针对以数据拟合为目标的非线性反演问题逐次线性化方法存在不适定性和计算量大等问题，提出了从模型的角度来研究非线性反演问题求解方法的新思路，建立了视模型空间对比方法．该方法避开了非线性问题直接线性化的处理，在求取模型参数修改量时有较好的适定性，并且计算量小．文中还研究了方法的收敛性条件和收敛速度，并用大地电磁测深（ＭＴ）合成数据进行了２９１２个参数的二维成像，对一条实测ＭＴ剖面进行了二维反演．结果证实了方法的有效性和解决大规模反演问题的能力．

青岛即墨东北部地区基岩富水区新探

通过对青岛即墨东北部大官庄乡这一剥蚀残丘垅岗区进行地下水资源探测，发现了延伸规模较大的东西向张性断裂破碎带，并在大旱之年的枯水季节（１９９７年１２月）于有利部位打出日产４００ｍ３的富水井，确认该断裂带富水。由此，为在青岛地区寻找基岩裂隙水开辟了新的断裂构造领域。

划分磁异常场的插值切割法在研究火成岩体分布中的应用

插值切割法通常被用于分离磁异常的区域场和局部场。本文进一步将插值切割法用于研究火成岩体的分布，因为大多数火成岩具有较强的磁性。其方法是：首先选取一较大的切割半径，对经过化权并消除干扰处理的实测磁异常进行分离，划分出由火成者引起的局部磁异常；然后用从小到大的一系列切割半径对局部磁异常进行抵值切割处理，得到相应的一系列异常幅值；当异常幅值出现稳定极值时，所对应的切割半径与该局部磁异常对应的火成岩体的中心埋深相当。这一规律是对球状火成岩体而言的。至于其它形状火成岩体，经过特殊处理，也能得到类似的结果。

磁(重力)异常局部场与区域场分离处理中的扩边方法研究

本文针对局部场与区域场分离处理的特点，提出了区域场扩边的思想，并给出了一种区域场扩边方法。与其它两种数据处理中常用的扩边方法进行了对比研究。模型和实测资料的处理结果表明，区域场扩边方法效果最好，对折扩边方法次之，余弦衰减至零的方法效果最差。

利用插值切割法研究磁性基底局部起伏特征

首先提出了明显改善插值切割法异常分离效果的改进切割方案，然后着重研究了该方法在研究盆地磁性基底局部起伏特征中的应用。模型与实际资料的研究表明，分离结果可以清晰直观地反映基底起伏特征，并在有利条件下可用来比较基底起伏幅度的相对大小，因而可以取得明显的地质效果，是盆地基底起伏特征研究的一种有效方法。

多次插值切割法分场的一个实例

根据地矿部航空物探遥感中心、胜利石油管理局的研究结果，胜利油田某区的磁性层可分为浅、中、深三层。本文采用了多次插位切割法对该区的航磁异常进行了划分。具体作法是：①首先采用4次插值切割法，计算出区域异常，然后用磁异常减去区域异常得到局都异常，这就是浅部异常。经过化极垂向一次导数和垂向二次导数的验证，吻合很好，说明浅部异常的可信度高。②将新的区域异常继续进行6次插值切割法，再次计算区域异常和局部异常，将此局部异常作为中部异常，区域异常作为深都异常。在浅、中、深部异常划分的基础上，作者进行了细致的解释工作，包括：①确定了第四系地层中的9处玄武岩、凝灰岩或干扰磁性源的分布范围、顶深和磁化率；②确定了沙三段地层中的玄武岩体的合理分布范围、顶深和磁化率，并得到区内两口钻井资料的证实；③推断了一条基底断层。

各向异性介质中波动方程有限元法模拟及其稳定性

利用六面体单元和三线性插值函数，推导出在任意弹性各向异性介质中三维三分量波动方程所满足的有限元方程．给出用时域有限差分求解该有限元方程时，时间、空间采样率和弹性参数应满足的稳定性条件．提出避免每一时间步都求解一大型线性代数方程组的方法．在开始求波场前，先求质量矩阵的逆的因子形式，以后每时间步只需作刚度矩阵的逆的因子形式与向量的乘积运算，可大大提高有限元素法在微机上进行波动方程数值模拟的运算速度．结合两个数值模拟算例，说明该方法的可行性，并与伪谱法作了计算速度的比较，有限元法的计算速度高于伪谱法的计算速度．

电导率分块连续变化的二维MT有限元模拟（Ⅰ）

用有限元法进行电导率分块连续变化的二维大地电磁的模拟。本文第一部分给出有限元法系数矩阵的计算公式，可用于编制有限元法的计算程序。模型上的计算结果和计算精度的分析将在本文的第二部分介绍。

电导率分块连续变化的二维MT有限元模拟（Ⅱ）

用有限单元法进行电导率分块连续变化的二维大地电磁场的模拟，本文的第一部分给出有限单元法系数矩阵的计算公式，这部分给出四个模型的计算结果。在有限元中，分别采用双线性插值和双二次插值。计算结果与解析法或有限差分法进行对比，吻合得很好。算例表明，本方法具有计算精度高、网格剖分简单。

大地电磁H_x型波二维地形改正的方法与效果

用边界无法计算大地电磁Hx型波二维地形影响，并用有限元法计算了起伏地形下水平层状介质的大地电磁响应，用边界元法的结果对畸变后的大地电磁响应进行地形改正，基本恢复了水平地形下层状介质的大地电磁响应的特点．本文给出了物理意义明确、形式简单的视电阻车和阻抗相位地形改正公式．

0-1编码遗传算法

本文分析了常规二进制编码遗传算法中二进制编码方法的特点,总结出二进制编码方法存在占用内存多、实现不灵活和译码运算量相对大的缺点,使较大规模的多参数优化问题难于用二进制编码遗传算法在较小内存的计算机上实现。为了克服二进制编码方法的这一缺点,我们提出一种0-1编码方法。文中介绍了0-1编码的方法和特点,并从定义的图式概念出发,证明了0-1编码遗传算法的收敛性。实际算例也表明,0-1编码遗传算法是可行的。

二维直流电阻率测深曲线的快速反演

本文采用作者提出的一维直流电阻率测深的直接反演公式，结合二维正演和渐近迭代反演，推出了一种二维直流电阻率测深的快速反演方法，经实际应用，效果很好。

点源二维电场中纵剖面上电位的计算

将总电位ν分解为正常电位．ｕｏ和异常电位ｕ·本文绘出点源二维电场中，异常电位ｕ的傅氏变换Ｕ的变分方程，以及点源二维电场地形影响的异常电位ｕ的傅氏变换Ｕ的积分方程，前者可用有限单元法求解，后者可用边界单元法求解．根据电源点的Ｕ的傅氏反交换，解决了前人尚未解决的直接计算点源二维电场中纵剖面（通过电源点、平行走向的剖面）的电位问题．

灾害损失等级划分的模糊灾度判别法

作者应用模糊模式识别理论于灾害损失等级划分的研究，提出了模糊灾度概念，并建立了模糊灾度等级的隶属函数，从而给出了一种可用于灾害损失等级划分的模糊灾度判别法。

垂直断层附近三维不均匀体点源电场的边界单元法

本文将给出垂直断层附近三维不均匀体电阻率法的数值解，所用方法是边界单元法，但基本解的选择与以往的边界单元法不同。

灾害损失定量评估的模糊综合评判方法

讨论了灾害损失的构成要素以及灾害损失定量评估的指标体系；依据灾害损失评估指标体系的一般形式，给出了由灾害损失单一因子计算灾级的公式；基于笔者提出的模糊灾度概念，给出了可用于各类灾害损失的等级划分的模湖综合评判方法。

岩浆岩储层的形成及油气藏类型

对岩浆岩这一特殊油气藏的储层形成及其孔隙结构特征进行了较系统的阐述 ,认为造成该类储层的决定性因素为构造运动、风化及热液作用等 ,形成了以次生的裂缝为主的复杂孔隙结构类型 ,其储渗空间的分布具有局限性和一定的分散性。由此 ,也造成该类油气藏在分类上的特殊性 ,即按其储层成因的主导因素可分为 :构造裂缝油气藏、古风化壳油气藏及火山碎屑岩油气藏等三大类、若干种 ;提出岩浆岩油气藏在地壳中分布的规律 ,为勘探和开发该类油气藏提供一定的理论基础和评价的依据。

区域均衡补偿模式在青藏高原地区的应用

采用区域补偿模式中的粘弹性板模型，利用相关函数技术，对青藏高原地区的均衡补偿间题进行研究，得到结果：一等效弹性厚度约为３０ｋｍ的粘弹性板的挠曲完成了较长波长上的区域补偿，未被补偿的短波长地形及地下负载由该板所承载．

用样条函数法求重力异常二阶垂向导数和向上延拓计算

对于二元三次样条函数的插值问题，本文给出了插值系数的求法及利用样条函数求导数、不定积分和定积分的统一计算公式，并用二元三次样条函数对地质模型进行了重力异常二阶垂向导数和向上延拓的计算。模型试其值和理论值的比较结果证明了本文提出的方法精度高且方便实用。

用边界元法模拟三维地形对MT场的影响

二维地形对大地电磁场(MT)的影响的模拟,已经解决。三维地形影响的模拟,至今未见文献报道。在三维问题中,H_x和E_x波的解法相同,本文只讨论H_x波的边界无法解。 1 边值问题 在均匀介质中,定态波(时间因子e^(-iwt))