CN108761565A - 航空重力测量的质量评估方法、装置、设备及存储介质 - Google Patents
航空重力测量的质量评估方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN108761565A CN108761565A CN201810295970.7A CN201810295970A CN108761565A CN 108761565 A CN108761565 A CN 108761565A CN 201810295970 A CN201810295970 A CN 201810295970A CN 108761565 A CN108761565 A CN 108761565A
- Authority
- CN
- China
- Prior art keywords
- value
- current
- airborne gravity
- measurement
- line
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明实施例提供了一种航空重力测量的质量评估方法、装置、设备及存储介质,其中,该方法包括:测量测区内各切割线的航空重力异常值;测量当前架次内当前测线的航空重力异常值;计算当前架次内当前测线与各切割线的交叉点处的航空重力异常不符值,进而计算当前架次内当前测线的原始测量精度值;对各切割线的航空重力异常值进行水平偏差校正,根据各切割线的校正后的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与各切割线的交叉点处的校正后的航空重力异常不符值;进而计算当前架次内当前测线的校正后的测量精度值;根据原始测量精度值和校正后的测量精度值,评估当前架次内当前测线的测量质量。
Description
技术领域
本发明涉及重力测量技术领域,特别涉及一种航空重力测量的质量评估方法、装置、设备及存储介质。
背景技术
航空重力测量是航空地球物理勘探方法之一,其过程是利用小型固定翼飞机、直升机或其他飞行平台搭载航空重力仪(airborne gravimeter)系统、GPS导航定位系统沿设计航线和高度进行地球重力场(地球重力场是地球周围空间任何一点存在的一种重力作用或重力效应,或为地球表面或附近一点处单位质量所受的重力,数值上等于重力加速度。重力场是空间中的一种力或力场,分布于地球表面及其邻近的空间;空间中任一质点都受到重力的作用。重力场是引力场和惯性离心力场的合成场)测量。通常航空重力测量即测量地球重力场的垂向分量。航空重力测量数据经过解算处理,可以得到沿测线的航空重力异常(重力异常(gravity anomaly),即由于地下物质密度分布不均匀引起的重力随空间位置的变化。在重力勘探中,将由于地下岩石、矿物质密度分布不均匀所引起的重力变化,或地质体与围岩密度的差异引起的重力变化,称为重力异常)值和沿切割线的航空重力异常值。
航空重力测量是一种涉及惯性导航、地球物理、GPS等多学科的综合性航空地球物理测量方法,航空重力测量数据通常用于大地水准面测量、地质构造分析、油气与矿产资源勘查等领域。与地面重力测量相比,航空重力测量不仅快速经济,而且能够在一些难以开展地面重力测量的特殊区域(如沙漠、冰川、沼泽、高山、森林、海陆交互地带等区域)进行作业。
由于航空重力测量技术复杂,测量过程有较多的不确定性,数据测量质量的量化评价与控制较困难。如果航空重力测区内地面重力测量工作充分,有较高质量的地面重力测量数据,利用航空重力测量数据与地面重力测量数据进行比对,可以达到检验、控制航空重力测量数据的质量的目的。但在实际工作中,航空重力测量主要用于在地面工作程度不高的海陆交互区或者沙漠戈壁等区域进行普查,这些测区内高质量的地面重力测量数据资料缺失或不易搜集。目前,在地面重力测量数据缺失的情况下,通常采用根据飞行高度、姿态、GPS数据质量、重复线(在航空地球物理测量中,重复线是指在同一条测线位置重复测量,所获取的多次重复测量的数据为重复线数据。重复线一般用来评价航空地球物理勘探仪器(如航空重力仪)动态测量的重复性和一致性,检测动态精度指标,或用于对某条测线进行重复检查)质量等辅助数据结合经验,人为判断航空重力测量数据的测量质量。缺少一种经济、可行的技术方法来量化考核航空重力测量数据的精度或质量。
目前,在航空重力测量过程中,通常用重复线来评价航空重力仪动态测量的重复性和一致性。通过统计重复线的航空重力异常内符合精度(该精度通过以下公式(1)计算得到),检测航空重力仪的动态精度,或用于对某条测线进行重复检查。
公式(1)中:j为重复线序号,j=1,2,3…m;i为重复线观测点序号,i=1,2,3…n;δij为第j条重复线上第i个观测点测量值与所有重复线在该位置平均值之差;εj为第j条重复线内符合精度。
如果使用重复线方法评估航空重力测量的精度或质量,要求对每条测线都要进行重复测量,即每条测线都需要重复飞行不少于2次,即每条测线都采用重复线检查的方法作业,这样成本将会翻倍,使得成本剧增,而且造成测量进度和效率急剧下降。因此,虽然该重复线方法可以用于评估航空重力测量的精度或质量,但是,由于成本、测量进度和效率等因素,重复线方法从未被采用,而且以后也不会被采用。
目前,在航空磁法测量过程中,利用测量数据的四阶差分均方根值的(如公式(2)所示)作为评价磁力仪噪声的指标,进而分辨测量数据中的错误数据,以便可以对错误数据进行校正或消除。
公式(2)中:n为参与计算的四阶差分值的总个数;Bi为第i点的四阶差分值,i=1,2,3…n;为n个四阶差分值的平均值;S为计算得到的噪声水平。
该四阶差分方法可以识别测量数据中的尖峰型、阶跃型等高频假值,如图1所示,使得可以评价测量数据的质量,进而实现对每条航磁测线的质量控制。四阶差分方法也是目前航空磁法测量领域中常用的对测线进行质量评价和控制的技术手段。
但是,四阶差分法只能鉴别高频假值,无法识别低频假值扰动。而在航空重力测量中存在量级较大的高频加速度扰动,在数据处理过程中必须采用大周期滤波处理,最终的航空重力测量数据是波长达数千米的低频信号(低于0.01Hz),也就是说,航空重力测量数据中已经不存在尖峰型、阶跃型等高频假值。因此,无法借鉴四阶差分法来对航空重力测量的质量进行评估。
目前,在航空磁法及航空重力测量等领域,通常使用如下公式(3)计算测量总精度:
式中:δ为测量总精度;δi为第i个切割线与测线的交叉点上的场值之差(即不符值);n为参加计算的切割线与测线的交叉点的个数。
但是,该测量总精度的计算方法只能是在对全测区的测量工作全部结束之后,整体评估整个测区的测量精度,无法实现在测量过程中及时对每个架次的每条测线的测量数据进行质量评估。
发明内容
本发明实施例提供了一种航空重力测量的质量评估方法,以解决现有技术中无法对每个架次的每条测线的测量数据进行质量评估的技术问题。该方法包括:
测量测区内各切割线的航空重力异常值;
测量当前架次内当前测线的航空重力异常值;
根据所述测区内各切割线的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的航空重力异常不符值;
根据计算得到的所有所述航空重力异常不符值,计算当前架次内当前测线的原始测量精度值;
对所述测区内各切割线的航空重力异常值进行水平偏差校正,根据所述测区内各切割线的校正后的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的校正后的航空重力异常不符值;
根据计算得到的所有所述校正后的航空重力异常不符值,计算当前架次内当前测线的校正后的测量精度值;
根据所述原始测量精度值和所述校正后的测量精度值,评估当前架次内当前测线的测量质量。
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述任意的航空重力测量的质量评估方法。以解决现有技术中无法对每个架次的每条测线的测量数据进行质量评估的技术问题。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行权上述任意的航空重力测量的质量评估方法的计算机程序。以解决现有技术中无法对每个架次的每条测线的测量数据进行质量评估的技术问题。
本发明实施例还提供了一种航空重力测量的质量评估装置,以解决现有技术中无法对每个架次的每条测线的测量数据进行质量评估的技术问题。该装置包括:
第一数据测量模块,用于测量测区内各切割线的航空重力异常值;
第二数据测量模块,用于测量当前架次内当前测线的航空重力异常值;
第一不符值计算模块,用于根据所述测区内各切割线的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的航空重力异常不符值;
第一精度值计算模块,用于根据计算得到的所有所述航空重力异常不符值,计算当前架次内当前测线的原始测量精度值;
第二不符值计算模块,用于对所述测区内各切割线的航空重力异常值进行水平偏差校正,根据所述测区内各切割线的校正后的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的校正后的航空重力异常不符值;
第二精度值计算模块,用于根据计算得到的所有所述校正后的航空重力异常不符值,计算当前架次内当前测线的校正后的测量精度值;
评估模块,用于根据所述原始测量精度值和所述校正后的测量精度值,评估当前架次内当前测线的测量质量。
在本发明实施例中,通过分别计算当前架次内当前测线的原始测量精度值和校正后的测量精度值,进而根据原始测量精度值和校正后的测量精度值,评估当前架次内当前测线的测量质量,即提出了一种在航空重力测量中对每条测线的测量质量进行评估的方法。本申请提出了低成本、易实施的评估架次内每条测线的测量质量的方法,有利于保障航空重力测量数据的质量;此外,由于结合了原始测量精度值和校正后的测量精度值来评估当前架次内当前测线的测量质量,且校正后的测量精度值相对而言更准确,使得有利于更准确地评估当前架次内当前测线的测量质量。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,并不构成对本发明的限定。在附图中:
图1是现有技术中的一种采用四阶差分法识别高频噪声的效果图;
图2是本发明实施例提供的一种航空重力测量的质量评估方法的流程图;
图3是本发明实施例提供的一种航空重力测量中计算线测的测量精度的流程图;
图4是本发明实施例提供的一种架次测线及全区切割线实测航迹的示意图;
图5是本发明实施例提供的一种航空重力仪测得的三个测区内测线的原始测量精度值的散点图;
图6是本发明实施例提供的一种航空重力仪测得的三个测区内测线的修正后的测量精度值的散点图;
图7是本发明实施例提供的一种航空重力仪测得的三个测区内测线的修正后的测量精度值的分布直方图;
图8是本发明实施例提供的一种上述航空重力测量的质量评估方法的具体实施流程图;
图9是本发明实施例提供的一种航空重力测量的质量评估装置的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施方式和附图,对本发明做进一步详细说明。在此,本发明的示意性实施方式及其说明用于解释本发明,但并不作为对本发明的限定。
在本发明实施例中,提供了一种航空重力测量的质量评估方法,如图2所示,该方法包括:
步骤201:测量测区内各切割线的航空重力异常值;
步骤202:测量当前架次内当前测线的航空重力异常值;即在测完测区内所有切割线的航空重力异常值之后,再逐架次测量每条测线的航空重力异常值。
步骤203:根据所述测区内各切割线的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的航空重力异常不符值;
步骤204:根据计算得到的所有所述航空重力异常不符值,计算当前架次内当前测线的原始测量精度值;
步骤205:对所述测区内各切割线的航空重力异常值进行水平偏差校正,根据所述测区内各切割线的校正后的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的校正后的航空重力异常不符值;
步骤206:根据计算得到的所有所述校正后的航空重力异常不符值,计算当前架次内当前测线的校正后的测量精度值;
步骤207:根据所述原始测量精度值和所述校正后的测量精度值,评估当前架次内当前测线的测量质量。
由图2所示的流程可知,在本发明实施例中,通过分别计算当前架次内当前测线的原始测量精度值和校正后的测量精度值,进而根据原始测量精度值和校正后的测量精度值,评估当前架次内当前测线的测量质量,即提出了一种在航空重力测量中对每条测线的测量质量进行评估的方法。本申请提出了低成本、易实施的评估架次内每条测线的测量质量的方法,有利于保障航空重力测量数据的质量;此外,由于结合了原始测量精度值和校正后的测量精度值来评估当前架次内当前测线的测量质量,且校正后的测量精度值相对而言更准确,使得有利于更准确地评估当前架次内当前测线的测量质量,测线的测量质量评估结果可以作为控制航空重力测量的依据。
具体实施时,如图3所示,可以通过以下步骤来计算当前架次内当前测线的原始测量精度值:
S1:计算当前架次内当前测线与横切该测线的各切割线的交叉点位置;
具体的,测区内横纵垂直相交的测线和切割线组成测网,其中,需要根据地质构造规划测点分布,通常垂直于地质构造方向安排测线(survey line),测线间距少则几十米,多则几公里;沿垂直于测线方向安排切割线(tie line,或称为控制线),切割线间距通常为5公里或10公里。如图4所示,L1、L2、L3、L4为测线,T1、T2...TN为全测区的切割线,p11为L1与T1的交叉点位置,p12为L1与T2的交叉点位置,以此类推;p21为L2与T1的交叉点位置,p31为L3与T1的交叉点位置,以此类推。可以根据投影坐标计算测线与切割线的交叉点位置。
S2:计算各切割线在交叉点位置上的航空重力异常值,并建立切割线上交叉点位置的航空重力异常值数据库;
具体的,航空重力测量的测区通常需要飞行多个架次才能完成测量,测线的测量精度按架次进行计算,即选择针对测线的一个飞行架次,计算该当前架次内各条测线的测量精度值,逐架次完成全部测线的测量精度值的计算。为了提高测量数据的准确性,在本实施例中,首先完成所有切割线的飞行测量,即按照图4中的T1、T2...TN等切割线进行飞行测量,可通过多个飞行架次完成所有切割线的飞行测量。交叉点通常不与测点重合,可以采取线性插值的方式计算交叉点位置上测线及切割线对应的航空重力异常值。
S3:计算当前测线在交叉点位置上的航空重力异常值,并建立测线上交叉点位置的航空重力异常值数据库;
具体的,完成所有切割线的飞行测量之后,再针对测线逐架次进行飞行测量。
S4:计算当前架次内当前测线与所述测区内各切割线的交叉点处的航空重力异常不符值,该航空重力异常不符值即交叉点处当前测线的航空重力异常值与切割线的航空重力异常值的差值;
S5:计算得到的所有航空重力异常不符值的标准差,该标准差即为当前架次内当前测线的原始测量精度值。具体的,还可以选择其他算法实现根据所有航空重力异常不符值计算当前架次内当前测线的原始测量精度值,本实施例优选的,选择计算所有航空重力异常不符值的标准差,并将该标准差作为当前架次内当前测线的原始测量精度值。
例如,可以通过以下公式(4)来计算所有航空重力异常不符值的标准差:
其中,i=1,2,3……N,N为该测线与切割线的交叉点总数;σ为标准差;δi为某测线上与切割线的第i个交叉点处的航空重力异常不符值。以图4中L1测线为例,δi计算公式为:
其中:为L1测线上p1i交叉点位置处的航空重力异常值;为Ti切割线上p1i交叉点位置处的航空重力异常值,i=1,2,3……N。
具体实施时,不同架次飞行时,由于航空重力仪器状态及飞行条件不同,如飞行高度、测量点位等存在差异,会造成航空重力仪器测量值存在水平差异/偏差。在一个架次飞行期间,由于飞行条件、航空重力仪器状态基本一致,测线间水平偏差量可以忽略。全测区切割线数据采集通常分多个架次完成,且架次时间间隔可能较长,所以全测区切割线与每个架次测线间通常会存在水平偏差。为了准确计算测线的测量精度,本申请提出利用每个架次内的测线数据对全测区的切割线数据进行水平偏差校正,以便根据测区各切割线的校正后的航空重力异常值和当前架次内当前测线的航空重力异常值,来计算得到更准确的当前测线的校正后的测量精度值。具体可以通过以下步骤来对测区内各切割线的航空重力异常值进行水平偏差校正:
针对每条切割线,计算该切割线与当前架次内各测线的交叉点处的航空重力异常不符值;
计算所有所述航空重力异常不符值的平均值,将所述平均值确定为该切割线的水平偏差校正量;
将该切割线的航空重力异常值减去所述水平偏差校正量,得到该切割线的校正后的航空重力异常值。
具体实施时,可以通过以下公式(6)来计算上述水平偏差校正量:
式中:δji T为切割线Tj上第i个交叉点位置处的航空重力异常不符值;i=1,2,3……N,N为与切割线Tj相交的测线数(也为切割线Tj上交叉点的个数);为水平偏差校正量。例如,以图4中切割线T1为例,δ1i T计算公式为:
其中:为切割线T1上pi1交叉点位置处的切割线T1的航空重力异常值,为Li测线上pi1交叉点位置处的测线Li的航空重力异常值。
通过以下公式(8)来对第j条切割线的数据进行水平偏差校正:
式中:为第j条切割线的原始航空重力异常值,为公式(6)计算得到的第j条切割线的水平偏差校正量,为第j条切割线的水平偏差校正后的航空重力异常值。
具体实施时,得到测区内各切割线的校正后的航空重力异常值之后,在交叉点处当前架次内当前测线的航空重力异常值与各切割线的校正后的航空重力异常值之差即为校正后的航空重力异常不符值,然后,可以通过公式(4)来计算当前架次内当前测线的校正后的测量精度值,此时,将δi替换为当前测线上与切割线的第i个交叉点处的校正后的航空重力异常不符值代入公式(4)即可。
具体实施时,得到原始测量精度值和校正后的测量精度值后,可以人为依据经验来根据原始测量精度值和校正后的测量精度值评估当前架次内当前测线的测量质量,为了进一步提高质量评估的准确性,还可以将原始测量精度值和校正后的测量精度值与预存值比较大小,以评估当前架次内当前测线的测量质量。
具体实施时,为了进一步提高质量评估的准确性,本实施例提出了通过以下步骤来实现根据所述原始测量精度值和所述校正后的测量精度值,评估当前架次内当前测线的测量质量:
获取预设等级评估规则,其中,所述预设等级评估规则包括原始测量精度值和校正后的测量精度值与质量等级的对应关系;
将所述原始测量精度值和所述校正后的测量精度值与所述预设等级评估规则进行匹配,得出当前架次内当前测线的测量质量的等级。
具体的,预设等级评估规则中原始测量精度值和校正后的测量精度值可以设置为单独的阈值,来直接与当前测线的所述原始测量精度值和所述校正后的测量精度值比较大小;也可以这样设置,例如,测线的原始测量精度值小于预设值,对应于校正后的测量精度值设置有多个数值区间,不同的数值区间对应不同的质量等级,其中,质量等级的级别与校正后的测量精度值的大小成反比。该预设等级评估规则可以通过表格的形式存在,例如,如下表1所示。
测线的测量精度(单位:mGal,毫伽) | 质量等级 |
校正后的测量精度≤0.8,且原始测量精度≤3.5 | 一级 |
0.8<校正后的测量精度≤1.0,且原始测量精度≤3.5 | 二级 |
1.0<校正后的测量精度≤1.5,且原始测量精度≤3.5 | 三级 |
校正后的测量精度>1.5,或原始测量精度>3.5 | 作废 |
表1
具体的,预设值可以根据不同实际需求来确定,主要用于根据原始测量精度值判断该测线的测量是否为正常测量,例如,表1中预设值为3.5。可以根据不同实际需求针对校正后的测量精度值设置不同个数值区间(例如,三个数值区间、四个数值区间等),数值区间的上限值和下限值也可以根据不同实际需求来确定。
具体实施时,为了提高预设等级评估规则的准确性,在本实施例中,可以通过以下步骤来预先得到上述预设等级评估规则:
针对某种型号的航空重力仪(获取预设等级评估规则与测量当前测线的过程使用相同型号的航空重力仪),搜集不少于3个测区(或工作量不少于6万公里)的实测数据,建立测线、切割线的大样本数据库。如果没有实测数据或实测数据样本不充分,则需要待完成几个测区或实测样本充分后再建立数据库。
例如,俄罗斯生产的GT-1/2A型航空重力仪在国内已经应用超过10年,在国内多个地区完成了大面积测量。可以搜集国内三个GT-1/2A航空重力测区的实测数据建立大样本数据库,三个测区编号为A1、A2、A3。其中,A1为近海海域测区,A2、A3为以平缓地形为主的内陆测区,使用飞机均为小型固定翼飞机,飞机飞行速度、高度均一致。按架次统计每条测线的原始测量精度和水平偏差校正后的测量精度,每条测线的原始测量精度和水平偏差校正后的测量精度的散点图如图5、6(图5、6中横坐标的测线的序号,纵坐标为测量精度值)所示。
从图5、6不难看出,原始测量精度由于水平偏差的干扰,精度值较离散,而且有许多严重偏离平均值的样本。由于水平偏差校正基本消除了测线与切割线间由于测量时间不同存在的整体水平差异(个别样本消减幅度超过2mGal),使得水平偏差校正后的测量精度值更能反映真实的测量精度。因此,使用水平偏差校正后的测量精度值数据库作为建立质量等级表的大样本数据库。水平偏差校正后的测量精度值分布如图7以及下表2所示。
表2
根据测线的测量精度值的分布情况,本申请提出以水平偏差校正后的测量精度值为主、以原始测量精度值为辅建立GT-1/2A航空重力仪的测线质量等级表(即上述预设等级评估规则,如表1所示)。
以下详细介绍上述航空重力测量的质量评估方法的应用流程,如图8所示,该流程包括以下步骤:
首先,通过以下步骤获取预设等级评估规则。
针对某一型号航空重力仪,获取多个测区合格的实测数据,作为建立质量等级标准(即预设等级评估规则)的航空重力测线、切割线的大样本数据库;
利用航空重力测线、切割线的大样本数据库,按架次统计每条测线与横切测线的所有切割线的交叉点位置处的航空重力异常不符值的均方差,得到每条测线的测量精度(该测量精度包括原始测量精度和校正后的测量精度),形成航空重力测线的测量精度值的大样本数据库;
对航空重力测线的测量精度值的大样本数据库进行统计,根据统计规律对测量精度值进行分档,建立航空重力测线的测量质量等级表(即预设等级评估规则)。
其次,应用上述航空重力测量的质量评估方法对每条测线进行质量评估。
在测区飞行开始后,首先安排完成所有切割线的飞行测量。为了获取高质量的切割线数据,可以采用重复飞行的方式,即每条切割线至少飞行两次,利用重复测量求平均的方式提高切割线数据质量,建立全区高质量切割线数据库。
完成所有切割线的飞行测量后,进行测线飞行架次。每完成一个测线飞行架次,利用单架次的测线的测量数据库及切割线的测量数据库,计算每条测线的原始精度值和水平偏差校正后的测量精度值,将测线的原始精度值和水平偏差校正后的测量精度值与质量等级表比对,以进行质量评估,确定出每条测线的质量等级。例如,某架次飞行4条测线,L1010、L1020、L1030、L1040,计算的测量精度如表3所示。
表3
与表1对比后,4条测线的质量等级划分如下:
L1010:一级;
L1020:一级;
L1030:三级;
L1040:作废,需要重飞。
每飞行一个测线架次,按如上流程进行质量评估,将质量合格数据分级并存入数据库存档,作废测线记录,安排日后飞行,直至所有测线飞行测量完成。
在本实施例中,提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述任意的航空重力测量的质量评估方法。
在本实施例中,提供了一种计算机可读存储介质,所述计算机可读存储介质存储有执行权上述任意的航空重力测量的质量评估方法的计算机程序。
基于同一发明构思,本发明实施例中还提供了一种航空重力测量的质量评估装置,如下面的实施例所述。由于航空重力测量的质量评估装置解决问题的原理与航空重力测量的质量评估方法相似,因此航空重力测量的质量评估装置的实施可以参见航空重力测量的质量评估方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图9是本发明实施例的航空重力测量的质量评估装置的一种结构框图,如图9所示,该装置包括:
第一数据测量模块901,用于测量测区内各切割线的航空重力异常值;
第二数据测量模块902,用于测量当前架次内当前测线的航空重力异常值;
第一不符值计算模块903,用于根据所述测区内各切割线的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的航空重力异常不符值;
第一精度值计算模块904,用于根据计算得到的所有所述航空重力异常不符值,计算当前架次内当前测线的原始测量精度值;
第二不符值计算模块905,用于对所述测区内各切割线的航空重力异常值进行水平偏差校正,根据所述测区内各切割线的校正后的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的校正后的航空重力异常不符值;
第二精度值计算模块906,用于根据计算得到的所有所述校正后的航空重力异常不符值,计算当前架次内当前测线的校正后的测量精度值;
评估模块907,用于根据所述原始测量精度值和所述校正后的测量精度值,评估当前架次内当前测线的测量质量。
在一个实施例中,所述第二精度值计算模块,包括:
不符值计算单元,用于针对每条切割线,计算该切割线与当前架次内各测线的交叉点处的航空重力异常不符值;
校正量计算单元,用于计算所有所述航空重力异常不符值的平均值,将所述平均值确定为该切割线的水平偏差校正量;
校正单元,用于将该切割线的航空重力异常值减去所述水平偏差校正量,得到该切割线的校正后的航空重力异常值。
在一个实施例中,所述评估模块,包括:
获取单元,用于获取预设等级评估规则,其中,所述预设等级评估规则包括原始测量精度值和校正后的测量精度值与质量等级的对应关系;
评估单元,用于将所述原始测量精度值和所述校正后的测量精度值与所述预设等级评估规则进行匹配,得出当前架次内当前测线的测量质量的等级。
在一个实施例中,在所述预设等级评估规则中,测线的原始测量精度值小于预设值,对应于校正后的测量精度值设置有多个数值区间,不同的数值区间对应不同的质量等级,其中,质量等级的级别与校正后的测量精度值的大小成反比。
在另外一个实施例中,还提供了一种软件,该软件用于执行上述实施例及优选实施方式中描述的技术方案。
在另外一个实施例中,还提供了一种存储介质,该存储介质中存储有上述软件,该存储介质包括但不限于:光盘、软盘、硬盘、可擦写存储器等。
本发明实施例实现了如下技术效果:通过分别计算当前架次内当前测线的原始测量精度值和校正后的测量精度值,进而根据原始测量精度值和校正后的测量精度值,评估当前架次内当前测线的测量质量,即提出了一种在航空重力测量中对每条测线的测量质量进行评估的方法。本申请提出了低成本、易实施的评估架次内每条测线的测量质量的方法,有利于保障航空重力测量数据的质量;此外,由于结合了原始测量精度值和校正后的测量精度值来评估当前架次内当前测线的测量质量,且校正后的测量精度值相对而言更准确,使得有利于更准确地评估当前架次内当前测线的测量质量。
显然,本领域的技术人员应该明白,上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种航空重力测量的质量评估方法,其特征在于,包括:
测量测区内各切割线的航空重力异常值;
测量当前架次内当前测线的航空重力异常值;
根据所述测区内各切割线的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的航空重力异常不符值;
根据计算得到的所有所述航空重力异常不符值,计算当前架次内当前测线的原始测量精度值;
对所述测区内各切割线的航空重力异常值进行水平偏差校正,根据所述测区内各切割线的校正后的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的校正后的航空重力异常不符值;
根据计算得到的所有所述校正后的航空重力异常不符值,计算当前架次内当前测线的校正后的测量精度值;
根据所述原始测量精度值和所述校正后的测量精度值,评估当前架次内当前测线的测量质量。
2.如权利要求1所述的航空重力测量的质量评估方法,其特征在于,对所述测区内各切割线的航空重力异常值进行水平偏差校正,包括:
针对每条切割线,计算该切割线与当前架次内各测线的交叉点处的航空重力异常不符值;
计算所有所述航空重力异常不符值的平均值,将所述平均值确定为该切割线的水平偏差校正量;
将该切割线的航空重力异常值减去所述水平偏差校正量,得到该切割线的校正后的航空重力异常值。
3.如权利要求1或2所述的航空重力测量的质量评估方法,其特征在于,根据所述原始测量精度值和所述校正后的测量精度值,评估当前架次内当前测线的测量质量,包括:
获取预设等级评估规则,其中,所述预设等级评估规则包括原始测量精度值和校正后的测量精度值与质量等级的对应关系;
将所述原始测量精度值和所述校正后的测量精度值与所述预设等级评估规则进行匹配,得出当前架次内当前测线的测量质量的等级。
4.如权利要求3所述的航空重力测量的质量评估方法,其特征在于,在所述预设等级评估规则中,测线的原始测量精度值小于预设值,对应于校正后的测量精度值设置有多个数值区间,不同的数值区间对应不同的质量等级,其中,质量等级的级别与校正后的测量精度值的大小成反比。
5.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至4中任一项所述的航空重力测量的质量评估方法。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至4中任一项所述的航空重力测量的质量评估方法的计算机程序。
7.一种航空重力测量的质量评估装置,其特征在于,包括:
第一数据测量模块,用于测量测区内各切割线的航空重力异常值;
第二数据测量模块,用于测量当前架次内当前测线的航空重力异常值;
第一不符值计算模块,用于根据所述测区内各切割线的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的航空重力异常不符值;
第一精度值计算模块,用于根据计算得到的所有所述航空重力异常不符值,计算当前架次内当前测线的原始测量精度值;
第二不符值计算模块,用于对所述测区内各切割线的航空重力异常值进行水平偏差校正,根据所述测区内各切割线的校正后的航空重力异常值和当前架次内当前测线的航空重力异常值,计算当前架次内当前测线与所述测区内各切割线的交叉点处的校正后的航空重力异常不符值;
第二精度值计算模块,用于根据计算得到的所有所述校正后的航空重力异常不符值,计算当前架次内当前测线的校正后的测量精度值;
评估模块,用于根据所述原始测量精度值和所述校正后的测量精度值,评估当前架次内当前测线的测量质量。
8.如权利要求7所述的航空重力测量的质量评估装置,其特征在于,所述第二精度值计算模块,包括:
不符值计算单元,用于针对每条切割线,计算该切割线与当前架次内各测线的交叉点处的航空重力异常不符值;
校正量计算单元,用于计算所有所述航空重力异常不符值的平均值,将所述平均值确定为该切割线的水平偏差校正量;
校正单元,用于将该切割线的航空重力异常值减去所述水平偏差校正量,得到该切割线的校正后的航空重力异常值。
9.如权利要求7或8所述的航空重力测量的质量评估装置,其特征在于,所述评估模块,包括:
获取单元,用于获取预设等级评估规则,其中,所述预设等级评估规则包括原始测量精度值和校正后的测量精度值与质量等级的对应关系;
评估单元,用于将所述原始测量精度值和所述校正后的测量精度值与所述预设等级评估规则进行匹配,得出当前架次内当前测线的测量质量的等级。
10.如权利要求9所述的航空重力测量的质量评估装置,其特征在于,在所述预设等级评估规则中,测线的原始测量精度值小于预设值,对应于校正后的测量精度值设置有多个数值区间,不同的数值区间对应不同的质量等级,其中,质量等级的级别与校正后的测量精度值的大小成反比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810295970.7A CN108761565B (zh) | 2018-03-30 | 2018-03-30 | 航空重力测量的质量评估方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810295970.7A CN108761565B (zh) | 2018-03-30 | 2018-03-30 | 航空重力测量的质量评估方法、装置、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108761565A true CN108761565A (zh) | 2018-11-06 |
CN108761565B CN108761565B (zh) | 2019-11-12 |
Family
ID=63981258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810295970.7A Expired - Fee Related CN108761565B (zh) | 2018-03-30 | 2018-03-30 | 航空重力测量的质量评估方法、装置、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108761565B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109443386A (zh) * | 2018-11-21 | 2019-03-08 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种实时高精度重力补偿方法 |
CN110334174A (zh) * | 2019-07-15 | 2019-10-15 | 武汉中海庭数据技术有限公司 | 一种pos数据纠偏方法及系统 |
CN113311505A (zh) * | 2021-04-13 | 2021-08-27 | 中国自然资源航空物探遥感中心 | 航空重力仪器飞行高度质量评估方法及装置 |
CN114019578A (zh) * | 2021-09-22 | 2022-02-08 | 中国科学院空天信息创新研究院 | 一种超导航磁梯度张量重复测线一致性检验方法和装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060117848A1 (en) * | 2004-12-07 | 2006-06-08 | Felix Raffalt | Method and apparatus for material identification by means of gravitational field analysis |
CN102608668A (zh) * | 2011-12-19 | 2012-07-25 | 华中科技大学 | 重力梯度测量系统及其测量方法 |
CN102901989A (zh) * | 2011-07-29 | 2013-01-30 | 中国国土资源航空物探遥感中心 | 一种基于重力场或磁场数据的地质体三维可视化建模与解释方法 |
CN103364841A (zh) * | 2013-07-29 | 2013-10-23 | 中国人民解放军国防科学技术大学 | 一种用于航空重力测量中星座跳变误差的平滑消除方法 |
CN103852799A (zh) * | 2014-02-25 | 2014-06-11 | 中国人民解放军92859部队 | 一种基于ppp技术的船载重力测量方法 |
CN104618033A (zh) * | 2015-01-08 | 2015-05-13 | 东南大学 | 一种多层自适应形态滤波重力信号噪声抑制方法 |
-
2018
- 2018-03-30 CN CN201810295970.7A patent/CN108761565B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060117848A1 (en) * | 2004-12-07 | 2006-06-08 | Felix Raffalt | Method and apparatus for material identification by means of gravitational field analysis |
CN102901989A (zh) * | 2011-07-29 | 2013-01-30 | 中国国土资源航空物探遥感中心 | 一种基于重力场或磁场数据的地质体三维可视化建模与解释方法 |
CN102608668A (zh) * | 2011-12-19 | 2012-07-25 | 华中科技大学 | 重力梯度测量系统及其测量方法 |
CN103364841A (zh) * | 2013-07-29 | 2013-10-23 | 中国人民解放军国防科学技术大学 | 一种用于航空重力测量中星座跳变误差的平滑消除方法 |
CN103852799A (zh) * | 2014-02-25 | 2014-06-11 | 中国人民解放军92859部队 | 一种基于ppp技术的船载重力测量方法 |
CN104618033A (zh) * | 2015-01-08 | 2015-05-13 | 东南大学 | 一种多层自适应形态滤波重力信号噪声抑制方法 |
Non-Patent Citations (1)
Title |
---|
郭志宏 等: "航空重力重复线测试数据质量评价方法研究", 《地球物理学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109443386A (zh) * | 2018-11-21 | 2019-03-08 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种实时高精度重力补偿方法 |
CN109443386B (zh) * | 2018-11-21 | 2022-06-28 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种实时高精度重力补偿方法 |
CN110334174A (zh) * | 2019-07-15 | 2019-10-15 | 武汉中海庭数据技术有限公司 | 一种pos数据纠偏方法及系统 |
CN110334174B (zh) * | 2019-07-15 | 2021-04-09 | 武汉中海庭数据技术有限公司 | 一种pos数据纠偏方法及系统 |
CN113311505A (zh) * | 2021-04-13 | 2021-08-27 | 中国自然资源航空物探遥感中心 | 航空重力仪器飞行高度质量评估方法及装置 |
CN113311505B (zh) * | 2021-04-13 | 2021-10-22 | 中国自然资源航空物探遥感中心 | 航空重力仪器飞行高度质量评估方法及装置 |
CN114019578A (zh) * | 2021-09-22 | 2022-02-08 | 中国科学院空天信息创新研究院 | 一种超导航磁梯度张量重复测线一致性检验方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108761565B (zh) | 2019-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108761565B (zh) | 航空重力测量的质量评估方法、装置、设备及存储介质 | |
Font et al. | Hypocentre determination offshore of eastern Taiwan using the Maximum Intersection method | |
CN103824133B (zh) | 一种花岗岩型铀矿田远景区综合预测方法 | |
McEvilly et al. | The Parkfield, California, earthquakes of 1966 | |
CN101371165B (zh) | 用于地球物理勘测的地形校正方法 | |
US20130223187A1 (en) | Geological Structure Contour Modeling and Imaging | |
NO335931B1 (no) | System og fremgangsmåte for automatisert reservoarutpeking | |
CN108333627B (zh) | 火成岩地区断裂真伪的识别方法及装置 | |
EA019807B1 (ru) | Системы обработки геофизических данных | |
CN105158796B (zh) | 确定toc含量的方法和装置 | |
SA520410962B1 (ar) | تقدير متعدد داخلي للظواهر الزلزالية المنعكسة | |
CN103562753A (zh) | 重力梯度计勘测技术 | |
CN107436452A (zh) | 基于概率神经网络算法的烃源岩预测方法及装置 | |
CN106443776A (zh) | 一种基于时间切片法的海底地震仪重定位方法 | |
CN109375283A (zh) | 一种砂岩储层3d渗透率演化史的分析方法 | |
Ayala | A new compilation of gravity data over the Iberian Peninsula and surrounding areas | |
CN103954996A (zh) | 一种基于旅行时法确定地层裂隙裂缝走向的方法及装置 | |
CN110244363A (zh) | 一种用于预测缝洞型储集层资源量的方法 | |
Ray et al. | Characterizing and modeling natural fracture networks in a tight carbonate reservoir in the Middle East: A methodology | |
CN100349014C (zh) | 重力勘探数据处理变密度地形校正方法 | |
Gumert | Airborne gravity measurements | |
CN106054252A (zh) | 一种叠前时间偏移的方法及装置 | |
Çıvgın et al. | Investigating the consistency of online earthquake catalogs of Turkey and surroundings | |
Miroslav et al. | Tutorial: The gravity-stripping process as applied to gravity interpretation in the eastern Mediterranean | |
CN107015289B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191112 Termination date: 20200330 |