重力地改数据处理系统

地质部为了编制1∶100万全国区域重力异常图,要求各省、市、自治区的物探队在85年以前完成本地区的区域重力测量工作。该项任务包括对旧有重力资料的收集、整理和数据处理,以及对广大空白区逐块填测,在处理新旧重力资料时都要进行地形改正值的计算,而这是一项十分繁重的工作。通常,为了计算一个重力测点的地改值需要用到其周围几百个,甚至上千个点的高程数据进行数值积分,而区域重力测量在一个省区中的重力测点是数以千计的。     

航空重力地改最大半径的研究及地改快速计算方法比较

航空重力地形改正是获得航空布格重力异常的重要环节,是航空重力勘探数据处理中的重点和难点问题。本文针对航空重力特点,分析了地改最大半径的选择与地形特点及计算精度的关系。为满足大数据量网格数据的计算要求,对全分辨率地形剖分方法、远近分区地形剖分方法、自适应四叉树地形剖分方法对航空重力地形改正计算的产生影响进行了对比。其中自适应四叉树地形剖分法可以对地形网格距离和高程进行综合考虑,达到最佳分辨率的地形剖分,既保证计算精度,又提高运算速度。     

重力勘探不规则测网的一种地改电算方法

一、概述消除地形起伏的影响是山区重力勘探的资料整理中最繁重的工作。随着电子计算机的广泛应用,传统的手算方法将很快被电算所代替。目前,已有多种大比例尺、规则测网的重力地改电算方案问世,对小比例尺不规则网的重力地改电算方法也正在纷纷的试验中。我队结合1/(10万)的重力找盐工作,边生产边探索不规则测网的地改电算方案,七七年底初步试验成功。对43个重力点进行了地改电算,地改范围算至30公里,与手算对比的均方差为±0.15mgl,七八年,地改电算方法投入生产,共计算965个点,随机抽29     

结点高程误差对重力地改的影响

本文目的在讨论结点高程误差对重力地形改正的影响;分析按查图法、数据图法确定结点高程可能达到的精度,以便采用适当方法测量结点高程,满足重力地改的需要.一结点高程误差对重力地改影响的计算公式1.按圆形域地改一般公式进行分析     

浅水域及其边缘重力测点的地改问题

随着国民经济的发展,区域重力工作已由纯陆域扩展到了江、河、湖、滨海等浅水域及其边缘陆地地区,这对区域重力工作在方法技术上提出了新的要求。这类工作目前主要用于石油资源勘查,其布格异常总精度要求较高,一般为100～200×10~(-8)m.s~(-2),有的甚至高于100×10~(-8)m.s~(-2)。虽然其工作比例尺属区域重力范围,但精度要求之高已达到大比例尺重力工作的精度,一般称之为重力细测工作,具有两个特点: 作为区域重力工作,根据“五统一”的     

一个计算重力地改值的 TI-59程序

近中区重力地改值计算,是一项繁杂的工作。本文在充分利用TI-59型程序计算器内存的基础上,设计了由六十八个直立扇形柱块组成六个环带,进行双重嵌套、多重循环的计算程序。程序采用一次输入、系统运算、中间显示与自动打印的方案,大幅度的简化了近中区圆域地改的计算工作。为了便于掌握,文中除绘出了计算公式、程序框图及使用方法外,并附有算例。     

方形域重力地改方法技术改进的探讨

在重力地改工作中,取数任务极为繁重,长期满足不了野外生产需要.因此,改进地改工作方法,提高精度和效率,是当前重力勘探急需研究的一个重要课题.本文首先简要讨论方形域重力地改方法,然后根据我们几年的实践经验,提出以下几种改进方法:     

空间相关阵的一种快速算法

空间因子分析和空间聚类分析是近五年发展起来的新方法，对于区域化变量的数据处理质量的提高有着重要意义。它们的标准算法运算速度极慢，影响了其实际应用。为此本文介绍一种快速计算空间相关阵以提高空间因子分析和空间聚类分析运算速度的算法。     

航测解析法在山区重力地改中的应用

山区重力勘探中,最繁重的工作是地形改正,它直接影响重力工作的进展和成果的质量。自1966年以来,我们与一些兄弟单位对山区重力工作方法进行了试验研究,在测地工作中曾先后引进过地面和航空立体摄影测量技术。     

一种高频面波频散函数的快速算法—改进的Abo-Zena法

作者在本文中首先介绍了改进的Anas Abo-Zena传递矩阵法面波频散函 数的计算问题。这是一种高频时稳定的算法，适用于一般地基勘查和无损检测的面波频散函数计算。其次，还讨论了用Monte Carlo法求解面波频散函数的问题。并用这种算法和Newton迭代法进行了对比，显示了对于Newton迭代法不能计算的速度逆转剖面（即层速度自上而下逐层减递），Monte Carlo法也适用。     

不完全伽玛函数的快速算法

不完全伽玛函数的快速算法吴明官，李彦兴（黑龙江省水利勘测设计院）一、前言在概率统计中皮尔逊Ⅲ型曲线是比较常见的线型（以下简称Ｐ─Ⅲ型），尤其在工程水文的设计暴雨、洪水、年径流等频率计算中的应用更为广泛。应用Ｐ─Ⅲ型曲线的关键就是如何推求离均系数Φ＿ｐ...     

基于并行预处理算法的三维重力快速反演

随着地球物理设备和探测技术的不断发展,快速处理大规模地球物理数据的需求也随之增长。为了解决三维重力数据密度反演的耗时问题,提出一种并行的预处理共轭梯度算法来提高计算效率。本文分别采用两种不同的预处理算子通过组合模型数据反演进行测试比较,并利用迭代残差和计算用时共同评价其加速效果。结果表明：对称逐次超松弛预处理方法比对角预处理方法反演计算速度快,密度结果更贴近实际模型;与传统串行的共轭梯度算法相比,本文并行预处理快速算法可以获得近19倍的加速比。将该算法应用于美国Vinton盐丘的实测重力数据中,反演结果能够很好地圈定出岩体的位置,验证了本文并行预处理共轭梯度法在三维重力数据快速反演中的高效性和可行性。     

三角形网格自动生成的一种快速算法

提出了一种三角形网格自动生成的快速算法:从边界点出发逐步向内部联网,不仅联网速度快,而且很容易处理凹边界。     

估计重力地改最佳地层密度值的面积相关法

为了较准确地选取重力地改中的地层密度值,针对“剖面法”中选取剖面条数的局限性,笔者采用“面积法”,考虑到了全区范围内地层密度对重力值的影响,使得选用的地层密度值更具代表性.该方法在四川西部某重力工区重力地改中的应用中取得了较好效果,极大地降低了地改后的布格重力异常与地形的相关性.该方法通过选用一定数值区间内的不同密度值,对实测重力值重新进行改算和地形改正,得到不同地层密度下的布格重力异常,借用重磁相关分析的方法原理,对布格重力异常与地形高程进行相关性分析,选取相关性最小的地层密度作为该区地形改正密度.     

区域重力地改“测点高程平移法”精度探讨

在讨论测点高程平移法区域重力地改的精度问题基础上，提出了系统差概念，并推荐一种新的计算方法，即把测点高程平移法和节点高程插值法有机地结合起来实现地改，通过理论模型和实例计算，表明该方法计算精度优于节点高程插值法和测点高程平移法。     

一种简易快速的时域地震波衰减补偿算法

针对地震波的能量随着传播时间的增加而递减,提出了一种时间补偿方法。通过调试与传播时间有关的两个可调参数进行地震波振幅补偿。模型试验证实了该方法可以有效地补偿能量的损失,而且在地层中地震波振幅衰减与品质因子Q值有关时,其中的一个可调参数可以用另一个可调参数和Q值来表示。也就是说,本文方法同样是一种反Q滤波方法。     

一种新的基数为2的DFT快速算法

本文介绍一种基于γ－循环阵快速分解方法的、新的基数为２的ＤＦＴ快速算法；该算法与具有最小乘法数的Ｐｒｅｕｓｓ算法比较，新算法的乘法数接近Ｐｒｅｕｓｓ算法，但其加法数仅为它的７５％。在Ｓｕｎ－ＳＰＡＲＣＳｔａｔｉｏｎ２工作站上，采用Ｆｏｒｔｒａｎ７７语言编程实现，结果表明：新算法在阶数（Ｎ）小于１０２４时，速度比分裂基（Ｓｐｌｉｔ—Ｒａｄｉｘ）算法快。如果进一步考虑算法中的固定常数乘法，那么新算法将会更加有效。