郯庐古岛弧构造系与金矿成矿的关系

本文从郯庐断裂带整体地质构造特征出发，应用板块构造理论和遥感技术、地质调查、物化探等综合研究方法，初步提出郯庐断裂带应归属于东亚大陆古岛孤－新裂谷式巨型构造系统，它与现代西太平洋岛弧具有“相对应性”及“等结性”控制金矿成矿集中区的新认识，并据此探索和总结了金矿（含其他矿产）的成矿规律，为寻找“招掖”式金矿进行成矿前景预测。     

用于速度图象重建的层析成象法

本文提出了一种用于速度图象重建的层析成象法。特别注意到地震学上的ST与医学上的CT的不尽相同之处,并处理了由此引起的困难。与速度图象重建有关的正问题的分析表明,不同频率范围的资料其分辨能力是不同的,因此,区域地震与远震资料的解释应予注意。本文有关速度图象重建的讨论不仅涉及到与成象方法有关的反演理论和方法,还给出了进行数值计算的内存量和运算量的分析。关于重建图象的质量评价,还给出了反演解的可靠性分析,为分辨特性的描述提供了一种既简便又形象直观的方案。 本文的方法与现行主要方法的比较表明:1.本文的方法允许存在速度间断面,有利于揭示不同构造区域的地壳厚度的差异;2.在常用的方法中,作为正问题假定的常速度块同反演结果解释时非常速度块之间不协调。由于我们把给定网格点的速度值的插值函数作为速度的空间函数,这一困难已被解决;3.相对于Aki等的一般方法,本文的内存量节省约一个量级、运算量减少约一半;4.相对于ART类方法本身不能给出成象结果的可靠性分析,本文的方法则具有显著的优点。 作者在分析了天然地震资料的状况之后相信,对充分利用我国地震台网的现有观测资料,本文的方法是有效的。     

电阻率测深在地电异常判定中的应用

塘沽台NS方向地电阻率观测值在观测极距缩短后出现大幅度上升变化，为了判断此变化的性质及其产生的可能原因，开展了多次恢复原极距的对比观测。并依据该台已有的电测深资料，对此变化进行了最佳拟合的理论反演计算。经比测及反、正演计算结果分析，认为该台NS方向大幅度的上升变化其主要原因与台站观测极距缩短有关，可能是随观测极距的缩短，受浅部电性结构层影响增大，而与深部电性结构无关，因而这一变化不具有震兆意义。     

Gravity/magnetic potential of uneven shell topography

 A fast spherical harmonic approach enables the computation of gravitational or magnetic potential created by a non-uniform shell of material bounded by uneven topographies. The resulting field can be evaluated outside or inside the sphere, assuming that density of the shell varies with latitude, longitude, and radial distance. To simplify, the density (or magnetization) source inside the sphere is assumed to be the product of a surface function and a power series expansion of the radial distance. This formalism is applied to compute the gravity signal of a steady, dry atmosphere. It provides geoid/gravity maps at sea level as well as satellite altitude. Results of this application agree closely with those of earlier studies, where the atmosphere contribution to the Earth's gravity field was determined using more time-consuming methods. Received: 14 August 2000 / Accepted: 19 March 2001     

变形监测中几个涉及测距问题的探讨

本文以平面位移监测为背景,阐述了边长前方交会法精度中的一个歧意,介绍了测距气象改正数的补偿计算方法,以及为便于比较各期距离观测值的变化而设计的仪器定高器。该文对高精度测距工作具有一定的参考价值。     

频率域线源大地电磁法有限元正演模拟

本文介绍了频率域线源大地电磁法有限元正演模拟的研究结果. 在外边界上统一应用适合于人工源的一阶吸收边界条件来形成边值问题,可减小基于平面波假设造成的人为截断边界的影响. 程序编辑中设计了两个二维数组分别存储总体系数矩阵的非零元素和在总体结点编号中的位置,使内存占用量减少,且物理意义明确,方便用高斯-赛德尔等迭代法解有限元方程时调用. 采用视δ函数模拟线源,提高了解方程组的稳定性. 最后通过对1个简单模型和1个复杂模型的模拟,证明所用的方法对异常体能够有明显的反映,说明了该方法的可靠性和有效性.     

一种由一维反演初步解释二维电阻率构造的方法

利用对非水平层状介质的地质结构首先进行二维模型的正演模拟,尔后对模拟曲线再进行一维反演来确定层状结构特征,实现由一维反演初步解释二维电阻率构造的解释过程。通过对三个不同结构特征的二维模型的解释结果表明,采用这种方法能够解释简单的二维电阻率构造,尤其是对于确定存在横向电阻率差异的界面效果更好;而且利用一个简单的经验公式在等值线断面图上通过求解截止电阻率值达到确定异常体的几何形态     

布格重力异常建模反演在二连盆地油气勘探中的应用

二连盆地由众多小凹陷组成,伊和凹陷是其中之一.重力资料在盆地"探边摸底"、"定凹选带"等方面发挥着重要作用.本文基于布格重力异常资料,利用2.5D建模反演技术对其中横跨主洼槽的两条测线进行可视化建模及正演拟合.其中,初始模型的几何参数取自二维地震解释结果;物性参数中的纵向密度取自测井声波波速与密度的相关公式.正演拟合为建模反演技术的关键环节,密度横向变化,是在拟合过程中通过人机交互试错调整实现.通过正演拟合,揭示了主洼槽沉积层横向和纵向均存在较大的密度不均匀性,主洼槽有利于油气生成.与地震解释结果对比,二者在赛汉组-第四纪、腾格尔组深度相差不大,但凹陷底部拟合结果与地震解释结果有明显差异.通过钻井结果验证,正演拟合各沉积层深度与钻井揭示的各层深度结果一致,二维地震解释结果偏浅,其可能原因是时深转换选取的速度值偏小造成.     

基于混合阶矢量基函数的海洋可控源电磁三维谱元法数值模拟

在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction(IDR(s))迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究.     

基于流体动力学数值模拟的海水层反射地震研究

由于多种海水运动同时存在以及地形的影响,海水层结构非常复杂,解释海水层地震相,分析海水运动过程是地震海洋学研究的新方向.本文提出结合流体动力学数值模拟与反射地震正演分析海水层地震相的方法.首先,对地形和流体建模,得到特定条件下流体运动状态;然后用反射地震正演将模拟获得的海水层温盐剖面转换为反射地震数据;进一步和实际测量得到的地震海洋学剖面进行对比,分析地形、海水运动对海水层地震相的影响.以内孤立波浅化过程为例,通过流体动力学数值模拟,获得其浅化过程中出现的下沉型、分裂、转换型三个阶段的海水层剖面;对温盐剖面进行反射地震正演,分析浅化不同阶段海水层反射几何形态、反射结构等特征.这种新方法有望解释复杂地震海洋学图像,深化海底地形对海水运动影响的认识.