青藏高原及其邻区岩石层三维密度结构

搜集了青藏高原及其邻近区域的S波速度三维层析成像结果和2万多个实测重力点资料，将重力资料进行各种改正并网格化为30′×30′的布格重力异常．首先采用密度差与S波速度差之间的经验关系式，建立青藏地区岩石层密度的初始模型，再利用布格重力异常进行阻尼最小二乘法反演，得到青藏地区岩石层三维密度分布结果．反演结果表明：（1）青藏高原岩石层密度分布不仅在纵向上不均匀，而且在横向存在明显的不均匀．在深度10－70km范围内，高原整体呈低密度特性，在50－70km深度范围内低密度特征更加突出，与周缘地区存在0，1g／cm³的密度差．而在90－110km深度范围内，高原岩石层地幔显示密度高．（2）岩层密度分布与大地构造有明显相关的分区性，显示出青藏块体、巴颜喀拉块体、塔里木块体和印度块体．     

提高断裂构造研究水平之我见

根据近些年来国内外断裂构造及相关领域研究进展，阐述了提高断裂构造研究水平及扩展其应用前景的途径。主要讨论了利用多源地学信息，多侧面地识别断裂构造；借助多层次地学资料，从三维和动态上研究断裂构造；加强流变学和应力场研究，深化断裂动力学分析；提高断裂研究定量化水平，扩展应用前景。     

地下水稳定运动基本微分方程的三维解法──地下水向承压完整井的运动

以映射和叠加原理为基础，探讨了均质、等厚承压含水层中地下水向完整井作三维运动时，地下水稳定运动基本微分方程的求解方法，并与裘布衣公式作了对比。由于求解时考虑了地下水流速的水平、垂直分量，故所得结果较裘布衣公式更为准确与合理。     

郑州北郊水源地地下水流三维有限元模拟

本文从傍河水源地地下水流动的三维特征分析入手，结合三维有限元数值模型，模拟了两组大型抽水试验，识别和反求含水层参数，进一步对地下水开采方案进行了评价。     

地层模型的三维实体显示与剖切

本文介绍了地层模型的三维实体显示的原理与方法，同时又提出了沿任意方向垂直剖切地层模型的一种实用有效算法，它已用于交互式三维地层矿床显示系统ＧＤ－３ＤＷＳ，是三维地质解释的很好工具     

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

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

三维分块倾斜界面的反演及其应用

介绍了利用反射波走时反演介质层速度和包含直立断层的三维界面的计算方法．对各层界面利用多个倾斜或平界面方程来描述．给出了数值计算实例．结果表明，反演解与真模型很接近，说明计算方法是有效的．利用该方法处理了唐山地震区的ＰｍＰ反射波资料，获得了该区域莫霍界面的三维分块形态．唐山附近莫霍界面错断与该区域１９７６年唐山地震及一系列余震的发生有密切的联系．     

电法勘探的发展和展望

电法勘探在经历了近一个世纪的发展后，其方法理论、仪器设备、野外数据采集、处理和解释等方面都经历了一系列重大变化．本文以方法理论的进展为主线，回顾、展望了目前电法勘探中几个重要而令人关注的研究焦点．这些问题的研究进展将会对21世纪的电法勘探产生深远的影响．     

三维溶质运移问题的分步广义迎风解法

对对流占优的三维溶质运移问题提出了分步广义迎风解法，首先利用Ｎ．Ｎ．Ｙａ－ｎｅｎｋｏ对水动力弥散方程分步求解的思想，将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题，对对流定解问题采用广义迎风对偶单元均衡法求解，对扩散定解问题采用一般的Ｇａｌｅｒｋｉｎ有限元法求解，不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题，而且避免了求节点速度这一步，节省运算步骤，对井点的浓度变化给出了更合适的求解方法。     

中国中西部地区地震的重新定位和三维地壳速度结构

采用中国中西部地区(21~36N， 98~112E)193个地震台在1992～1999年间记录到的9 988次地震的Pg和Sg震相走时的读数资料，用Roecker的SPHYPIT90程序，反演了该地区三维地壳P波速度结构，并用SPHREL3D90程序进行了地震的重新定位．反演结果揭示了中国中西部地区地震P波速度结构明显的横向不均匀性，这些不同深度上波速的横向变化多以该地区的活动断裂为分界线．可以看出活动断裂两侧存在明显的速度反差．通过重新定位，得到了6 459次地震的震源参数，这些精确定位的地震震中明显沿该区活动断裂呈现条带状分布，其范围和尺度清晰地表示了这一地区地震活动与活动断裂的紧密关系．其中，82％重新精确定位的事件的震源深度在20 km以内. 这一结果与笔者用双差地震定位法得到的重新定位的震源深度分布相一致．