散乱离散点数据的三角形网格化快速成图

提出了一种基于三角网格的等值线成图线性插值方法。常规等值线成图都是基于矩形网格进行网格化插值的,三角网格能够更好地逼近地球物理场的形态和散乱离散点数据的边界,得到的等值线图更加光滑。通过搜索边界、三角网格化剖分、线性插值、搜索等值线、Bezier曲线光滑等值线等5个步骤,可以对任意散乱离散点数据进行快速成图。实际数据的成图结果表明:该方法插值效果好,不进行数据外推,得到的等值线图能直接反映散乱离散点数据的空间位置且成图速度快,可大大提高实际工作效率。     

岩土体渗流自由面问题的重心插值无网格方法

岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。     

基于单位分解的扩展径向点插值无网格法

针对线弹性断裂力学问题,提出扩展径向点插值无网格法（X-RPIM）。该方法基于单位分解思想,在传统径向点插值无网格法的位移模式中加入扩展项来描述裂纹两侧的不连续位移场和裂尖奇异场。由于其形函数具有Kronecker ? 函数性质,易于施加本质边界条件。详细描述了X-RPIM不连续位移模式的建立,支配方程的离散形式以及J积分计算混合模式裂纹的应力强度因子的实现过程,讨论了不同积分区域对应力强度因子的影响。数值算例分析证明了该方法在求解断裂问题时的可行性和有效性,同时说明扩展径向点插值无网格法在模拟裂纹扩展问题时具有良好的前景。     

波动方程垂向插值斜缆采集正演模拟

为了解决海洋斜缆地震数据正演模拟时大多数检波器不在网格点上,不能直接记录波场值的问题,在使用有限差分法求解波动方程时,增加了垂向插值计算,得到了实际检波点的波场值。通过建立精细模型对正演模拟结果进行了误差分析,验证了该垂向插值方法的可靠性。分析对比了正演模拟数据和实际生产数据的波场特征,进一步证明了该波动方程与垂向插值正演模拟方法的正确性。该方法在倾斜电缆地震采集参数设计方面,具有较高的实用价值。     

基于紧支径向基函数的点插值无网格方法

紧支径向基函数能使支配方程中的刚度矩阵具有稀疏性,很适合应用于无网格方法中,其缺点是在插值计算时精度不高.点插值方法的插值函数具有Delta函数性质,可以很方便的施加本质边界条件,但在计算插值函数时矩阵易出现奇异.为了提高计算精度并避免点插值法的局限性,首先对紧支径向基函数进行完备性修正,然后用完备性修正的紧支径向基函数代替多项式来形成插值函数,建立了紧支径向基函数点插值方法.由于该方法中的形函数满足Delta函数性质,因此本质边界条件可以像传统的有限元方法一样很容易施加.然后将该方法用于二维弹性静力问题的求解,导出了其相应的离散方程.最后将该方法应用于一个悬臂梁的分析中,初步验证了该方法的有效性与合理性.     

二维起伏地表直流电场插值法数值模拟

本文提出了一种采用插值法实现边界条件的有限差分数值模拟新方法, 地表函数通过三次样条插值能很好的拟合真实地表形状。其地表边界条件的实现思想是在地表法线方向上寻找一个与地表点位于同一等位线上的点, 该点通过周围点插值来建立离散方程。采用不规则四边形距离加权插值法使插值形式统一化。在水平地表情况下做的数值解与解析解的对比分析表明: 该方法的计算精度适合要求, 对山谷模型的实例计算也得到了视电阻率与地形成镜像对称的结果。     

地层三维建模及面模型插值方法对比研究

通过面模型进行实体建模是构建三维地层模型的有效方法之一,为了更好的选取面模型插值方法,本文对距离幂次反比方法、克里格插值方法以及离散光滑插值方法进行了对比研究。实例结果显示,离散光滑插值方法既可以保持插值区面的光滑,又能与原始数据最大程度的拟合,较其他两种插值方法对于表现局部的非连续空间特征具有更大优势。同时本文通过将离散光滑插值方法与地质剖面信息相结合,通过面模型与三维实体模型之间的空间运算来获得倾入、尖灭等地层之间的空间关系,并结合地层的年代信息,建立了一套可以减少人为干预,提高建模精度和合理性的三维地层建模流程。实践证明,该流程能较好应用于安徽省钟姑矿田白象山矿区的地层三维建模工作。     

非规则网格下二维浅水流动的数值模拟

建立了一种在非规则结构化网格上求解平面二维浅水流动的有限体积方法。通过采用地形在离散网格内双线性变化及离散网格界面间地形连续的地形逼近方法和应用可以有效处理间断问题的Roe格式来离散浅水方程中的对流项,并通过VanLeer提出的状态插值法提高格式精度。在计算原始变量在网格内的插值梯度时,采用最小二乘方法求变量的最优梯度代替差分计算梯度,从而可采用任意形状的不规则四边形网格离散计算域。计算实例表明,该方法能够计算间断问题并能够处理各种复杂流态的过渡,具有较好适应性和计算精度,能够满足不同实际问题的计算要求。     

弹性波场和稳恒电流场的Frechet导数计算方法

在地球物理参数反演中，观测的物理场量，例如波场中的位移或电流场中的电位，对于介质参数的编数的计算是非常重要的。本文采用自伴随算子和格林函数的方法给出了弹性波场和稳恒电流场问题的Ｆｒｅｃｈｅｔ导数解析式。该式表明，Ｆｒｅｃｈｅｔ导数可由格林函数Ｇ及其梯度ΔＧ求得。把任意场源正演问题的解和Ｆｒｅｃｈｅｔ导数的计算统一起来，并提出用有限元方法计算任意背景场下的格林函数，从而得到Ｆｒｅｃｈｅｔ导数的数值解     

有理逼近及其在数值优化中的应用

文章利用倒差商研究了数值有理逼近的理论和方法，证明了有理插值涵数的唯一性，并应用该方法求解数值优化问题，建立了直线搜索的计算方案。计算结果表明，用有理插值方法求解数值优化问题是速度快、精度高的算法。     

快速样条函数插值网格化法

网格化是物化探数据处理的重要环节。本文采用首次提出的快速样条函数插值方法对二维离散数据进行曲面插值，进而作网格化计算。算例表明，本文网格化方法精度高、计算速度快。     

大型对称变带宽方程组的Cholesky分解法

作者针对地球物理数值模拟中常碰到的大型稀疏变带宽方程组的求解问题,介绍了 大型稀疏变带宽矩阵的存储方法及Ｃｈｏｌｅｓｋｙ分解法,该方法的特点是用二个一维数组,其中一个输入时存储对称稀疏矩阵变带宽内的元素,输出时存储Ｃｈｌｏｅｓｋｙ下三角短阵带宽内的元素;另一个存储对称变带宽矩阵对角线元素在前一维数组中的位置,大大节约了所占的计算内存空间。     

超高密度电阻率法在岩溶及破碎带探测中的应用

结合工程实例,介绍了超高密度电阻率法在勘探岩溶及破碎带中的应用,总结了数据采集过程中的注意事项,并将成果与钻孔资料作了对比分析。结果表明,该方法在岩溶及破碎带的分布位置、大小、埋藏深度及发育情况等方面的勘察上都取得了较好的效果,是一种能有效提高工程勘察精度的物探方法。     

X射线荧光分析技术在四川龙塘铅锌矿成矿规律研究与资源评价中的应用

本文以四川龙塘铅锌矿为例,对X射线荧光分析的技术条件、资料处理、方法应用等方面进行较系统地研究,并把现场测量成果与地质研究密切结合,总结成矿规律,指导找矿预测和资源评价,取得良好效果。     

均匀电流场中极化体上方放电曲线形态及其与充电持续时间的关系

作者采用二阶线性系统来描述均匀电流场中形状简单的极化体的特征。这种地电断面对阶跃讯号的响应是两个指数函数的线性组合,研究表明,系统的两个时间常数只与地球物理参数有关,对整个观测剖面来说它们的值是不变的,两指数函数前之系数,除与地球物理参数有关外还与几何因子有关,在剖面上不同测点,两系数的数值与符号是可变化的。因此,可能有8种不同类型的放电曲线。作者还研究了充电持续时间对放电曲线的影响。理论计算结果可以很好地解释J．伯廷、J．洛布和傅良魁的野外实测资料和模型实验结果,并指出可能存在而尚未被发现的现象。作者认为本文所采用的理论研究方法可应用于更复杂的情况．     

X射线荧光分析技术在四川龙塘铅锌矿成矿规律研究与资源评价中的应用

本文以四川龙塘铅锌矿为例 ,对X射线荧光分析的技术条件、资料处理、方法应用等方面进行较系统地研究 ,并把现场测量成果与地质研究密切结合 ,总结成矿规律 ,指导找矿预测和资源评价 ,取得良好效果     

美国大陆地区震中图的构造解释

按照地震线的判定原则,利用美国大陆地区历史地震和现代地震的震中分布图,并结合地质和地球物理资料,画出了美国大陆地区的地震线和地震密集带,分析了美国大陆总体的和分区的地震构造格局特征,并与中国大陆的地震构造格局作了比较。结果说明,地震震中的空间分布特征反映了地壳内部活动构造的信息,大陆的地震构造格局不仅受原有构造格架和现今边界动力环境的控制,而且还可能受更高层次的全球构造动力因素的影响。     

用交流电阻率法探测煤矿导(含)水构造

分析了煤矿井下地球物理场的特殊性,指出电阻率法的频率域测量方法具有压制工业干扰能力强、克服电极接地条件差的优点,适用于煤矿井下作业环境的需求。煤矿井下实际探测结果说明,交流电阻率法在三维空间内异常定位准确,适用于煤层是否存有导(含)水构造的采前探测。     

基于信号重构的地震道插值

在简要回顾几种插值方法的基础上 ,提出了基于信号重构的地震道插值方法。这种方法无需任何地质和地球物理假设 ,主要用于恢复缺失的地震道。本方法的基础是傅立叶变换理论和最小平方反演技术。文中给出了原理、算法、实现步骤及简要的流程框图。如果对变换域的参数进行规则采样 ,则该方法可以用高效、快速的算法来实现。本插值方法适用于规则和非规则采样的数据。文中给出了理论合成数据和实际数据的试验结果 ,并对结果进行了讨论。同时 ,本方法对于 AVO分析、DMO处理、三维叠前偏移成像等研究也将是有益的     

塔河油田下奥陶统碳酸盐岩古岩溶储层研究--以塔河油田6区为例

通过对研究区FMI成像测井和XMAC，DSI测井等资料的综合，并结合该区的油气地质、地球物理资料，对下奥陶统缝洞型碳酸盐岩古岩溶储层的剖面结构及展布特征进行了较为深入地分析结果表明，研究区古岩溶储层主要分布于下奥陶统顶部不整合面之下约250m的范围内，储层垂向发育具有明显分带性，岩溶剖面的纵向分带与该区多期次岩淀旋回作用密切相关。古岩溶储层的发育与分布受多种因素控制。