CN112526625B - 航空重力测量点的布格重力异常值的计算装置 - Google Patents
航空重力测量点的布格重力异常值的计算装置 Download PDFInfo
- Publication number
- CN112526625B CN112526625B CN202011329110.4A CN202011329110A CN112526625B CN 112526625 B CN112526625 B CN 112526625B CN 202011329110 A CN202011329110 A CN 202011329110A CN 112526625 B CN112526625 B CN 112526625B
- Authority
- CN
- China
- Prior art keywords
- measurement point
- cuboid
- elevation
- gravity
- aviation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000005484 gravity Effects 0.000 title claims abstract description 203
- 238000005259 measurement Methods 0.000 title claims abstract description 121
- 230000002159 abnormal effect Effects 0.000 title abstract description 9
- 238000012937 correction Methods 0.000 claims abstract description 55
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 238000010276 construction Methods 0.000 claims description 6
- 238000000034 method Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 238000012876 topography Methods 0.000 description 11
- 238000009826 distribution Methods 0.000 description 8
- 238000001914 filtration Methods 0.000 description 8
- 239000011435 rock Substances 0.000 description 8
- 238000003860 storage Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/16—Measuring gravitational fields or waves; Gravimetric prospecting or detecting specially adapted for use on moving platforms, e.g. ship, aircraft
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Ocean & Marine Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种航空重力测量点的布格重力异常值的计算装置,包括:获取模块,用于获取航空重力测量点的地理位置信息、所述航空重力测量点的GNSS高程以及重力异常值;投影模块,用于将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,得到高程投影点;构建模块;地形改正值计算模块;布格重力异常计算模块,用于根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常。本发明提供的航空重力测量点的布格重力异常值的计算装置,可以获得航空布格重力异常。
Description
技术领域
本发明是涉及航空重力测量技术领域,特别是涉及一种航空重力测量点的布格重力异常值的计算装置。
背景技术
随着科技的飞速发展,在重力勘探领域重力测量的精度要求越来越高。
重力测量所得到的数据为重力观测数据。重力观测数据包括因地下密度不均匀的地质体引起的异常,也包含因各测点周围地形不同、所处纬度不同等因素的影响,因此需要对重力观测数据进行校正。而航空重力测量是重力测量的一种,具体是把航空重力测量***装在飞机上进行连续测量的一种重力测量方法。
目前,布格重力异常(Bouguer gravity anomaly)对重力数据的处理具有很重要的意义。地面布格重力异常通常通过地面重力测量获得。地面重力仪的观测结果经过纬度改正、高度改正,中间层改正和地形改正以后,再减去正常重力值后所得到的重力称为地面布格重力异常。
基于此,本申请的发明人发现,现有技术中通过地面重力测量值获取布格重力异常的计算方法在航空重力异常获取航空布格重力异常的计算中存在一定的不适用性。
公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
为解决上述问题,本发明实施例的目的在于提供一种航空重力测量点的布格重力异常值的计算装置。
第一方面,本发明实施例提供了一种航空重力测量点的布格重力异常值的计算方法,包括:获取航空重力测量点的地理位置信息、所述航空重力测量点的GNSS高程以及重力异常值;将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,得到高程投影点;在地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体,所述第一长方体的高度为四个节点高程值的平均值,所述地形高程数据网格中包括多个节点的坐标值构成的栅格数据,所述地形高程数据网格中包括多个节点的高程值;分别以所述四个节点作为起点,以地形高程数据网格中的网格步长作为边长,向背离投影点的方向延伸,构建第二长方体,所述第二长方体的高度为底面上四个节点对应的高程值的平均值,重复构建第二长方体的步骤,直至节点的坐标值超出预设范围为止;根据所述航空重力测量点到每一个长方体的高度差、所述航空重力测量点到每一个长方体的距离差、每一个长方体的密度以及万有引力常数计算每一个长方体的地形改正值,其中,所述每一个长方体为第一长方体或第二长方体;对所有地形改正值进行求和,获得该航空重力测量点的地形改正量;根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常。
在一种可能的实现方式中,所述地形高程数据网格包括局部地形高程数据网格以及区域地形高程数据网格;局部地形高程数据网格中节点坐标值的预设范围为0-20km,区域地形数据网格中节点坐标值的预设范围为20km-166.7km。
在一种可能的实现方式中,所述节点的坐标值超出预设范围为:节点的坐标值超出166.7km。
在一种可能的实现方式中,所述在地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体之前,还包括:将测量点的地理位置信息投影至岩石密度分布数据网格中,得到密度投影点,所述岩石密度分布数据网格中包括多个节点的坐标值构成的栅格数据,所述岩石密度分布数据网格中包括多个节点的密度值;所述第一长方体的密度为密度投影点对应的密度值。
在一种可能的实现方式中,所述第二长方体的密度值为所述第二长方体底面上四个节点分别投影至岩石密度分布数据网格中的平均值。
在一种可能的实现方式中,局部地形高程数据网格的划分规则为高精度和大比例尺;区域地形高程数据网格的划分规则为低等级精度和小比例尺。
在一种可能的实现方式中,所述计算方法还包括:获取下一航空重力测量点的地理位置信息、所述下一航空重力测量点的GNSS高程以及重力测量值,并重复上述对下一航空重力测量点对应的地形改正量进行计算的步骤,直至完成所有航空重力测量点的地形改正量的计算;对所有航空重力测量点的地形改正量进行逐点低通滤波,所述低通滤波的滤波尺度与航空重力测量点重力异常的滤波尺度相对应;所述根据该航空重力测量点的重力异常值以及地形改正量获得该航空重力测量点的布格重力异常值包括:根据航空重力测量点的重力异常值以及滤波后的地形改正量获得该航空重力测量点的布格重力异常值。
第二方面,本发明实施例还提供了一种航空重力测量点的布格重力异常值的计算装置,包括:获取模块,用于获取航空重力测量点的地理位置信息、所述航空重力测量点的GNSS高程以及重力异常值;投影模块,用于将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,得到高程投影点;构建模块,用于在地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体,所述第一长方体的高度为四个节点高程值的平均值,所述地形高程数据网格中包括多个节点的坐标值构成的栅格数据,所述地形高程数据网格中包括多个节点的高程值;所述构建模块还用于分别以所述四个节点作为起点,以地形高程数据网格中的网格步长作为边长,向背离投影点的方向延伸,构建第二长方体,所述第二长方体的高度为底面上四个节点对应的高程值的平均值,重复构建第二长方体的步骤,直至节点的坐标值超出预设范围为止;地形改正值计算模块,用于根据所述航空重力测量点到每一个长方体的高度差、所述航空重力测量点到每一个长方体的距离差、每一个长方体的密度以及万有引力常数计算每一个长方体的地形改正值,其中,所述每一个长方体为第一长方体或第二长方体;地形改正量计算模块,用于对所有地形改正值进行求和,获得该航空重力测量点的地形改正量;布格重力异常计算模块,用于根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常。
第三方面,本发明实施例还提供了一种存储介质,所述存储介质存储有计算机可执行指令,所述计算机可执行指令用于执行上述的航空重力测量点的布格重力异常值的计算。
本发明实施例提供的航空重力测量点的布格重力异常值的计算装置,通过将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体,分别以所述四个节点作为起点,重复构建长方体;计算每一个长方体的地形改正值,对所有地形改正值进行求和,获得该航空重力测量点的地形改正量;根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常,可以实现航空布格重力异常的计算。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出了本发明实施例所提供的航空重力测量点的布格重力异常值的计算方法的流程图;
图2示出了本发明实施例所提供的地形高程数据网格的示意图;
图3示出了本发明实施例所提供的长方体计算方法的示意图;
图4示出了本发明实施例所提供的航空重力测量点的布格重力异常值的计算装置的结构示意图。
具体实施方式
下面结合附图,对本发明的具体实施方式进行详细描述,但应当理解本发明的保护范围并不受具体实施方式的限制。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
除非另有其它明确表示,否则在整个说明书和权利要求书中,术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分,而并未排除其它元件或其它组成部分。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。
由于航空重力测量的实际地形起伏变化复杂,因此本实施例将航空重力测量点周围的地形剖分成许多长方体,计算出每个长方体地形质量对航空重力测量点的重力影响值,最后累加求和便得到该点的地形影响值。由于测量点周围的地形是极不规则的曲面,无法用一个空间坐标的函数式来表达,采用地形高程数据实现对地面地形的数字化模拟,并建立有规则的地形高程的网格数据。在进行航空重力地形校正计算时,在一定的范围内分别在局部地形高程网格数据和区域地形高程网格数据中按照网格节点来读取高程数据。
本发明实施例提供了一种航空重力测量点的布格重力异常值的计算方法的流程图,参见图1所示,包括:步骤1-步骤7。
步骤1,获取航空重力测量点的地理位置信息、所述航空重力测量点的GNSS高程以及重力异常值;
步骤2,将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,得到高程投影点。
在一种实现方式中,所述地形高程数据网格包括局部地形高程数据网格以及区域地形高程数据网格;局部地形高程数据网格中节点坐标值的预设范围为0-20km,区域地形数据网格中节点坐标值的预设范围为20km-166.7km。参见图2所示,其为地形高程数据网格的示意图,EFGH为局部地形区,IJKL为区域地形,中间黑色区域为航空重力测量点的范围。
局部地形高程数据网格的划分规则为高精度和大比例尺;区域地形高程数据网格的划分规则为低等级精度和小比例尺。
步骤3,在地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体,所述第一长方体的高度为四个节点高程值的平均值,所述地形高程数据网格中包括多个节点的坐标值构成的栅格数据,所述地形高程数据网格中包括多个节点的高程值;
步骤4,分别以所述四个节点作为起点,以地形高程数据网格中的网格长度作为边长,向背离投影点的方向延伸,构建第二长方体,所述第二长方体的高度为底面上四个节点对应的高程值的平均值,重复构建第二长方体的步骤,直至节点的坐标值超出预设范围为止;
其中,所述节点的坐标值超出预设范围为可以是节点的坐标值超出166.7km。
步骤5,根据所述航空重力测量点到每一个长方体的高度差、所述航空重力测量点到每一个长方体的距离差、每一个长方体的密度以及万有引力常数计算每一个长方体的地形改正值,其中,所述每一个长方体为第一长方体或第二长方体;
步骤6,对所有地形改正值进行求和,获得该航空重力测量点的地形改正量;
步骤7,根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常。
本实施例提供的航空重力测量点的布格重力异常值的计算方法,通过将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体,分别以所述四个节点作为起点,重复构建长方体;计算每一个长方体的地形改正值,对所有地形改正值进行求和,获得该航空重力测量点的地形改正量;根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常,可以实现航空布格重力异常的计算。
在一种实现方式中,步骤3之前还包括:将测量点的地理位置信息投影至岩石密度分布数据网格中,得到密度投影点,所述岩石密度分布数据网格中包括多个节点的坐标值构成的栅格数据,所述岩石密度分布数据网格中包括多个节点的密度值;所述第一长方体的密度为密度投影点对应的密度值。
相应的,所述第二长方体的密度值为所述第二长方体底面上四个节点分别投影至岩石密度分布数据网格中的平均值。
在一种实现方式中,步骤5可以具体包括:
针对每一个长方体,根据公式一计算该长方体的航空重力地形改正值Δg,所述公式一包括:
其中,G为万有引力常数,ρ为每个长方体的密度,x,y,z为观测点空间位置,(X1,X2)、(Y1,Y2)、(Z1,Z2)为长方体空间位置。R为观测点到长方体各个顶点的距离。
进一步的,步骤6以及步骤7的实现方式如下。参阅图3所示,其为本实施例提供的长方体计算方法的示意图,航空重力观测点投影到地形高程网格数据的位置为P(X0,Y0)点,其中A、B、C、D点为地形网格数据中距离P点最近的节点。以A、B、C、D点为长方体的边,其高度为节点A、B、C、D的高程值的平均值,计算长方体到航空重力测量点的地形影响值Δg1。向网格的两个方向扩展,以此类推计算每四个节点对应的拟合长方体到航空重力测量点A的地形影响值Δgi。
航空重力测量点的地形校正值之和为:
航空重力测量点的布格重力异常值为:
G布格=G自空-G地 (2)
其中G自空为航空重力测量点A的空间重力异常值,即通过航空重力测量获得。
在一种实现方式中,步骤7之后还可以包括:获取下一航空重力测量点的地理位置信息、所述下一航空重力测量点的GNSS高程以及重力测量值,并重复上述对下一航空重力测量点对应的地形改正量进行计算的步骤,直至完成所有航空重力测量点的地形改正量的计算。
具体的,针对每一航空重力测量点,重复上述步骤1-步骤6的计算步骤,直至所有航空重力测量点的地形改正量的计算。
在一种实现方式中,还可以包括:对所有航空重力测量点的地形改正量进行逐点低通滤波,所述低通滤波的滤波尺度与航空空间重力异常的滤波尺度相对应。
相应的,步骤7可以为:根据航空重力测量点的重力异常值以及滤波后的地形改正量获得该航空重力测量点的布格重力异常值。
本发明实施例提供还了一种航空重力测量点的布格重力异常值的计算装置的结构示意图,参见图4所示,包括:获取模块1、投影模块2、构建模块3、地形改正值计算模块4、地形改正量计算模块5以及布格重力异常计算模块6。
获取模块1用于获取航空重力测量点的地理位置信息、所述航空重力测量点的GNSS高程以及重力异常值;
投影模块2用于将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,得到高程投影点;
构建模块3用于在地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体,所述第一长方体的高度为四个节点高程值的平均值,所述地形高程数据网格中包括多个节点的坐标值构成的栅格数据,所述地形高程数据网格中包括多个节点的高程值;
所述构建模块3还用于分别以所述四个节点作为起点,以地形高程数据网格中的网格长度作为边长,向背离投影点的方向延伸,构建第二长方体,所述第二长方体的高度为底面上四个节点对应的高程值的平均值,重复构建第二长方体的步骤,直至节点的坐标值超出预设范围为止;
地形改正值计算模块4用于根据所述航空重力测量点到每一个长方体的高度差、所述航空重力测量点到每一个长方体的距离差、每一个长方体的密度以及万有引力常数计算每一个长方体的地形改正值,其中,所述每一个长方体为第一长方体或第二长方体;
地形改正量计算模块5,用于对所有地形改正值进行求和,获得该航空重力测量点的地形改正量;
布格重力异常计算模块6,用于根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常。
本实施例提供的航空重力测量点的布格重力异常值的计算装置,通过将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体,分别以所述四个节点作为起点,重复构建长方体;计算每一个长方体的地形改正值,对所有地形改正值进行求和,获得该航空重力测量点的地形改正量;根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常,可以实现航空布格重力异常的计算。
本发明实施例还提供了一种存储介质,所述存储介质存储有计算机可执行指令,其包含用于执行上述航空重力测量点的布格重力异常值的计算方法的程序,该计算机可执行指令可执行上述任意方法实施例中的方法。
其中,所述存储介质可以是计算机能够存取的任何可用介质或数据存储设备,包括但不限于磁性存储器(例如软盘、硬盘、磁带、磁光盘(MO)等)、光学存储器(例如CD、DVD、BD、HVD等)、以及半导体存储器(例如ROM、EPROM、EEPROM、非易失性存储器(NAND FLASH)、固态硬盘(SSD))等。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
前述对本发明的具体示例性实施方案的描述是为了说明和例证的目的。这些描述并非想将本发明限定为所公开的精确形式,并且很显然,根据上述教导,可以进行很多改变和变化。对示例性实施例进行选择和描述的目的在于解释本发明的特定原理及其实际应用,从而使得本领域的技术人员能够实现并利用本发明的各种不同的示例性实施方案以及各种不同的选择和改变。本发明的范围意在由权利要求书及其等同形式所限定。
Claims (1)
1.一种航空重力测量点的布格重力异常值的计算装置,其特征在于,包括:
获取模块,用于获取航空重力测量点的地理位置信息、所述航空重力测量点的GNSS高程以及重力异常值;
投影模块,用于将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,得到高程投影点;
构建模块,用于在地形高程数据网格中,根据距离所述高程投影点最近的四个节点作为底面构建第一长方体,所述第一长方体的高度为四个节点高程值的平均值,所述地形高程数据网格中包括多个节点的坐标值构成的栅格数据,所述地形高程数据网格中包括多个节点的高程值;
所述构建模块还用于分别以所述四个节点作为起点,以地形高程数据网格中的网格步长作为边长,向背离高程投影点的方向延伸,构建第二长方体,所述第二长方体的高度为底面上四个节点对应的高程值的平均值;重复构建第二长方体的步骤,直至节点的坐标值超出预设范围为止;
地形改正值计算模块,用于根据所述航空重力测量点到每一个长方体的高度差、所述航空重力测量点到每一个长方体的距离差、每一个长方体的密度以及万有引力常数计算每一个长方体的地形改正值,其中,所述每一个长方体为第一长方体或第二长方体;
地形改正量计算模块,用于对所有地形改正值进行求和,获得航空重力测量点的地形改正量;
布格重力异常计算模块,用于根据航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011329110.4A CN112526625B (zh) | 2020-06-08 | 2020-06-08 | 航空重力测量点的布格重力异常值的计算装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010511402.3A CN111679336B (zh) | 2020-06-08 | 2020-06-08 | 航空重力测量点的布格重力异常值的计算方法、装置 |
CN202011329110.4A CN112526625B (zh) | 2020-06-08 | 2020-06-08 | 航空重力测量点的布格重力异常值的计算装置 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010511402.3A Division CN111679336B (zh) | 2020-06-08 | 2020-06-08 | 航空重力测量点的布格重力异常值的计算方法、装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112526625A CN112526625A (zh) | 2021-03-19 |
CN112526625B true CN112526625B (zh) | 2024-01-30 |
Family
ID=72454001
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011329110.4A Active CN112526625B (zh) | 2020-06-08 | 2020-06-08 | 航空重力测量点的布格重力异常值的计算装置 |
CN202010511402.3A Expired - Fee Related CN111679336B (zh) | 2020-06-08 | 2020-06-08 | 航空重力测量点的布格重力异常值的计算方法、装置 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010511402.3A Expired - Fee Related CN111679336B (zh) | 2020-06-08 | 2020-06-08 | 航空重力测量点的布格重力异常值的计算方法、装置 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN112526625B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112462443B (zh) * | 2020-11-13 | 2021-07-20 | 中国自然资源航空物探遥感中心 | 用于航空重力测量中的同步地形改正方法及装置 |
RU2764108C1 (ru) * | 2020-12-03 | 2022-01-13 | Иван Владимирович Бригадин | Способ маскировки подземных объектов |
CN115629426A (zh) * | 2022-10-24 | 2023-01-20 | 中国自然资源航空物探遥感中心 | 一种提取局部重力异常的方法、***、装置及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104166170A (zh) * | 2014-08-13 | 2014-11-26 | 昆明理工大学 | 一种坑道重力全空间域定位探测低密度隐伏矿体的方法 |
WO2015104633A2 (en) * | 2014-01-13 | 2015-07-16 | Cgg Services Sa | System and method for two dimensional gravity modeling with variable densities |
CN110646858A (zh) * | 2019-10-12 | 2020-01-03 | 山东省物化探勘查院 | 一种海底重力测量中区—远一区地形改正计算方法 |
CN110706275A (zh) * | 2019-10-16 | 2020-01-17 | 中国石油大学(华东) | 一种基于卫星测高重力数据的局部重力异常提取方法 |
CN110967778A (zh) * | 2019-10-24 | 2020-04-07 | 西北大学 | 一种动态坐标系多面体剖分重力布格校正方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170108612A1 (en) * | 2015-10-15 | 2017-04-20 | King Saud University | Inertial system for gravity difference measurement |
CN108008459B (zh) * | 2017-11-28 | 2019-11-01 | 北京中科地物能源技术有限公司 | 一种获得剩余重力异常的方法及装置 |
-
2020
- 2020-06-08 CN CN202011329110.4A patent/CN112526625B/zh active Active
- 2020-06-08 CN CN202010511402.3A patent/CN111679336B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015104633A2 (en) * | 2014-01-13 | 2015-07-16 | Cgg Services Sa | System and method for two dimensional gravity modeling with variable densities |
CN104166170A (zh) * | 2014-08-13 | 2014-11-26 | 昆明理工大学 | 一种坑道重力全空间域定位探测低密度隐伏矿体的方法 |
CN110646858A (zh) * | 2019-10-12 | 2020-01-03 | 山东省物化探勘查院 | 一种海底重力测量中区—远一区地形改正计算方法 |
CN110706275A (zh) * | 2019-10-16 | 2020-01-17 | 中国石油大学(华东) | 一种基于卫星测高重力数据的局部重力异常提取方法 |
CN110967778A (zh) * | 2019-10-24 | 2020-04-07 | 西北大学 | 一种动态坐标系多面体剖分重力布格校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111679336B (zh) | 2020-12-15 |
CN112526625A (zh) | 2021-03-19 |
CN111679336A (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112526625B (zh) | 航空重力测量点的布格重力异常值的计算装置 | |
AU2016305571B2 (en) | System and method for gravity and/or gravity gradient terrain corrections | |
CN101839710B (zh) | 一种似大地水准面的确定方法 | |
CN102645209B (zh) | 机载LiDAR点云和高分辨率影像进行空间点的联合定位方法 | |
CN112363236B (zh) | 一种基于pde的重力场数据等效源延拓与数据类型转换方法 | |
CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
CN109934911B (zh) | 基于OpenGL的移动端高精度倾斜摄影三维建模方法 | |
CN103256914B (zh) | 一种基于dem计算淤地坝淹没面积的方法及*** | |
CN103278115A (zh) | 一种基于dem计算淤地坝淤积量的方法及*** | |
CN111257956A (zh) | 一种基于Matlab的区域似大地水准面精化方法 | |
CN113311494A (zh) | 一种卫星重力场反演方法 | |
CN112578450A (zh) | 基于频散介质标量波方程有限差分的正演模拟方法及装置 | |
Yildiz | A study of regional gravity field recovery from GOCE vertical gravity gradient data in the Auvergne test area using collocation | |
Kwon et al. | Gravity requirements for compensation of ultra-precise inertial navigation | |
CN113671570B (zh) | 一种地震面波走时和重力异常联合反演方法与*** | |
CN113109883B (zh) | 基于等参变换全球离散网格球坐标下卫星重力场正演方法 | |
CN102650703B (zh) | Cgcs2000板块模型cpm-cgcs2000构建方法 | |
Jekeli et al. | Geoid determination in South Korea from a combination of terrestrial and airborne gravity anomaly data | |
CN115112090A (zh) | 一种基于少量gnss水准点的山区数字高程基准建立方法 | |
CN113552644B (zh) | 密度确定方法、装置及电子设备 | |
Jekeli | Deflections of the vertical from full-tensor and single-instrument gravity gradiometry | |
US20240085182A1 (en) | Methods for lifting and lowering magnetic anomalies for magnetic navigation | |
Grzyb et al. | The effect of topography and quality of a digital terrain model on the accuracy of terrain corrections for centimetre quasigeoid modelling | |
CN114964169B (zh) | 像方物方协同改正的遥感影像平差方法 | |
CN104614730B (zh) | 一种利用月面激光反射阵坐标标校月球高程模型的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |