重力地改数据处理系统

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

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

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

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

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

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

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

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

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

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

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

高程数据网格间距对重力中区地形改正精度的影响

为了解高程数据网格间距对表面积分法、直立长方体法和平均高程直立长方体法计算的中区地形改正值精度的影响,笔者选择某地区450个测点,并使用不同网格间距高程数据计算中区地改值,通过对比发现表面积分法计算精度受高程数据网格间距影响较小,而直立长方体法反之。然后将中区地改50~2 000 m分为10个区间段进行计算,通过统计得出误差的45%和30%左右都分布在50~200 m和200~500 m段,因此提出提高中区地形改正精度必须提高50~200 m和200~500 m内高程数据网格密度。     

竖盘固定误差对三角高程的影响及对策

本文研究了竖盘固定误差对三角高程的影响，着重介绍测定方法及消除竖盘固定误差的措施。     

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

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

重力地改被积函数的一种快速算法

对于重力地形改正的计算,当积分公式与数字地形(即高程网)确定之后,其计算速度和精度主要取决于被积函数的处理方法。一、问题的提出众所周知,方域地形改正的计算公式为:     

重力异常与高程关系的统计分析

本文采用数理统计方法,对我国东南部约300万km²范围内的布伽重力异常和地形高程的关系进行了统计分析。通过对布伽异常和高程关系的讨论,提出了一种分离重力异常的方法。利用这种方法对某区与高程有关的重力异常和因地壳地质构造引起的剩余重力异常进行分解,编制了相应的图件。这为不同构造单元的密度结构特征及其可能的深部构造轮廓的推断提供了信息。     

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

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

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

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

大比例尺重力勘探GPS高程转换

在大比例尺重力勘探工作中,当布格重力异常设计总精度提高到0.025×10-5m/s2时,要求测点高程均方误差为0.05m.为使GPS高程能够达到该精度要求,提出利用EGM2008模型和GPS水准数据拟合局部似大地水准面的方法以实现测点GPS高程的转换,并利用某地区GPS水准资料,对转换方法的可行性进行验证分析.结果表明,这种方法完全适用于大比例尺重力勘探工作.     

重力地形改正误差的讨论

利用重力寻找金属矿(如铬铁矿),要求精度高,而且必须消除地形影响.本文简要讨论重力地改误差的主要来源及其大小. 众所周知,在直角坐标系中,重力地形改正公式为: ■■■■■为便于计算,将上式改用柱坐标表示: ■■■■■式中,f-万有引力常数;σ-围岩密度;r_m-内环半径;r_(m+1)-外环半径;h-扇形柱高度;R-近区半径为方便计算起见, 在近区,采用扇形锥地改公式:     

DEM网格间距对重力远区地改精度的影响及效果

1:5万重力远区地改一般由测区的1:5万DEM高程模型改正获得。1:5万DEM可按不同网格间距拼接而成,不同网格间距对应了不同的地改精度;RGIS是由自带高程库参与计算,完成测区1:5万重力远区地改。本文以山东省栖霞市臧家庄幅1:5万重力远Ⅰ区(2～20km)地改为例,通过采用25、50、100、200m四种网格节点距进行了1:5万DEM数据拼接和改进的双线性插值法重力远Ⅰ区地改、均方误差计算,并与RGIS自带高程库重力远Ⅰ区地改均方误差比较分析,证实基于1:5万DEM高程模型重力远Ⅰ区地改精度优于RGIS自带高程库重力远Ⅰ区地改精度。基于DEM改进的双线性插值法远区地改布格重力异常和高程相关度更高,对线性构造和地层、岩体的边界的识别精度高。     

青藏高原重力异常与地形高程的负相关校正

针对消除青藏高原重力异常与地形高程呈负相关关系所引起的假异常,提出了均衡校正和残余负相关校正相结合的一套具有可操作性的校正方法,并对青藏高原可可西里地区重力异常数据进行了处理,效果明显,验证了方法的有效性.