CN110553609A - 一种确定深海压力仪海深测量精度的方法及系统 - Google Patents
一种确定深海压力仪海深测量精度的方法及系统 Download PDFInfo
- Publication number
- CN110553609A CN110553609A CN201910999426.5A CN201910999426A CN110553609A CN 110553609 A CN110553609 A CN 110553609A CN 201910999426 A CN201910999426 A CN 201910999426A CN 110553609 A CN110553609 A CN 110553609A
- Authority
- CN
- China
- Prior art keywords
- sea
- satellite
- value
- sea level
- deep sea
- 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
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B13/00—Measuring arrangements characterised by the use of fluids
- G01B13/14—Measuring arrangements characterised by the use of fluids for measuring depth
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F25/00—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
- G01F25/20—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of apparatus for measuring liquid level
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开一种确定深海压力仪海深测量精度的方法及系统。该方法包括:获取深海压力仪的布设位置;获取测高卫星与深海压力仪同步观测期间的卫星轨道数据和深海压力仪在测高卫星过顶时间的海深测量数据;根据卫星轨道数据和深海压力仪的布设位置,获取海面高差;根据卫星轨道数据,获得点PCA处卫星测量海面高值;根据深海压力仪在测高卫星过顶时间的海深测量数据,采用差值计算,获取卫星测量海面高值对应时刻的深海压力仪的海深测量值;获得深海压力仪的测量标准差,得到深海压力仪进行海深测量时的精度。本发明可以实现深海压力仪相对海深测量精度的评估,且评估结果可靠性高。
Description
技术领域
本发明涉及海洋测绘领域,特别是涉及一种确定深海压力仪海深测量精度的方法及系统。
背景技术
随着人类对深蓝海精确化认知需求的不断提升,高精度大量程海深测量设备应运而生。近年来,海深测量设备的性能进一步提高,对数千米的海深进行观测时其测量分辨率已达到毫米量级。然而,海深测量设备的精度评定长期以来缺乏有效的手段。这主要包含两方面原因:一方面,用于深海海深变化监测的压力仪其测量精度难以通过模拟实验环境进行评估;另一方面,对于大量程深海海深测量设备而言,难以找到高精度海深基准值用于其测量精度评估。
因而,寻找对布设于深海海底的压力仪(简称深海压力仪)进行高质量精度评定的方法具有十分重要的现实意义。
发明内容
本发明的目的是提供一种确定深海压力仪海深测量精度的方法及系统,利用海洋测高卫星数据确定深海压力仪相对海深测量的精度,以实现深海压力仪相对海深测量精度的评估。
为实现上述目的,本发明提供了如下方案:
一种确定深海压力仪海深测量精度的方法,包括:
获取深海压力仪的布设位置;
获取测高卫星与所述深海压力仪同步观测期间的卫星轨道数据和深海压力仪在所述测高卫星过顶时间的海深测量数据;
根据所述卫星轨道数据和所述深海压力仪的布设位置,获取海面高差;所述海面高差为所述深海压力仪布设位置对应海面点CP与卫星地面轨迹中距离点CP最近的点PCA之间的海面高差;
根据所述卫星轨道数据,采用差值计算,获得点PCA处卫星测量海面高值;
根据所述深海压力仪在所述测高卫星过顶时间的海深测量数据,获取所述卫星测量海面高值对应时刻的所述深海压力仪的海深测量值;
根据所述海面高差、所述点PCA处卫星测量海面高值和所述深海压力仪的海深测量值,获得所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度。
可选的,所述测高卫星与所述深海压力仪同步观测期间的卫星轨道数据包括卫星的轨道高度、卫星到海面的高度计的测距值、干对流层延迟改正数据、湿对流层延迟改正数据、电离层延迟改正数据、负荷潮汐改正数据、地球固体潮汐改正数据、极潮改正数据和海况偏差改正数据。
可选的,所述根据所述卫星轨道数据和所述深海压力仪的布设位置,获取海面高差,具体包括:
根据所述卫星轨道数据和所述深海压力仪的布设位置,计算海面点CP与点PCA之间的距离;
采用全球平均海面高模型计算点CP的平均海面高梯度值;
根据所述海面点CP与点PCA之间的距离和所述点CP的平均海面高梯度值,利用公式ΔMSS=GradMSS×l确定所述海面高差;其中,ΔMSS为海面高差,GradMSS为点CP的平均海面高梯度值,l为海面点CP与点PCA之间的距离。
可选的,所述根据所述卫星轨道数据,采用差值计算,获得点PCA处卫星测量海面高值,具体包括:
利用SSHsat2=h-(R+Rdry+Rwet+Riono+Rssb+Rlt+Rst+Rpt)计算得到点PCA处卫星测量海面高值;其中,SSHsat为点PCA处卫星测量海面高值,h为卫星的轨道高度,R为卫星到海面的高度计的测距值,Rdry表示干对流层延迟改正数据,Rwet表示湿对流层延迟改正数据,Riono表示电离层延迟改正数据,Rlt表示负荷潮汐改正数据,Rst表示地球固体潮汐改正数据,Rpt表示极潮改正数据,Rssb表示海况偏差改正数据。
可选的,所述根据所述海面高差、所述点PCA处卫星测量海面高值和所述深海压力仪的海深测量值,获得所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度,具体包括:
根据所述海面高差和所述深海压力仪的海深测量值,利用公式SSHsat1=D+dmooring+ΔMSS计算与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值;其中,SSHsat1为与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值,D为所述深海压力仪的零位高度,dmooring为所述深海压力仪的海深测量值,ΔMSS为所述海面高差;
根据所述深海压力仪测量海面高值和所述点PCA处卫星测量海面高值,利用公式计算单次同步观测期间的观测差值;其中,为单次同步观测期间的观测差值,SSHsat2为点PCA处卫星测量海面高值;
根据多个同步观测周期中所述观测差值的变化序列,得到序列标准差值;
获取所述测高卫星测量海面高的精度值;
根据所述海面高差、所述序列标准差值和所述测高卫星测量海面高的精度值,利用公式计算所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度;其中,为所述深海压力仪的测量标准差,为序列标准差值,为测高卫星测量海面高的精度值,σΔMSS为海面高差的标准差。
本发明还提供一种确定深海压力仪海深测量精度的系统,包括:
布设位置获取模块,用于获取深海压力仪的布设位置;
观测数据获取模块,用于获取测高卫星与所述深海压力仪同步观测期间的卫星轨道数据和深海压力仪在所述测高卫星过顶时间的海深测量数据;
海面高差获取模块,用于根据所述卫星轨道数据和所述深海压力仪的布设位置,获取海面高差;所述海面高差为所述深海压力仪布设位置对应海面点CP与卫星地面轨迹中距离点CP最近的点PCA之间的海面高差;
卫星测量海面高值获取模块,用于根据所述卫星轨道数据,采用差值计算,获得点PCA处卫星测量海面高值;
深海压力仪的海深测量值获取模块,用于根据所述深海压力仪在所述测高卫星过顶时间的海深测量数据,获取所述卫星测量海面高值对应时刻的所述深海压力仪的海深测量值;
测量标准差获取模块,用于根据所述海面高差、所述点PCA处卫星测量海面高值和所述深海压力仪的海深测量值,获得所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度。
可选的,所述观测数据获取模块获取的测高卫星与所述深海压力仪同步观测期间的卫星轨道数据,包括卫星的轨道高度、卫星到海面的高度计的测距值、干对流层延迟改正数据、湿对流层延迟改正数据、电离层延迟改正数据、负荷潮汐改正数据、地球固体潮汐改正数据、极潮改正数据和海况偏差改正数据。
可选的,所述海面高差获取模块具体包括:
距离计算单元,用于根据所述卫星轨道数据和所述深海压力仪的布设位置,计算海面点CP与点PCA之间的距离;
平均海面高梯度值计算单元,用于采用全球平均海面高模型计算点CP的平均海面高梯度值;
海面高差确定单元,用于根据所述海面点CP与点PCA之间的距离和所述点CP的平均海面高梯度值,利用公式ΔMSS=GradMSS×l确定所述海面高差;其中,ΔMSS为海面高差,GradMSS为点CP的平均海面高梯度值,l为海面点CP与点PCA之间的距离。
可选的,所述卫星测量海面高值获取模块利用SSHsat2=h-(R+Rdry+Rwet+Riono+Rssb+Rlt+Rst+Rpt)计算得到点PCA处卫星测量海面高值;其中,SSHsat为点PCA处卫星测量海面高值,h为卫星的轨道高度,R为卫星到海面的高度计的测距值,Rdry表示干对流层延迟改正数据,Rwet表示湿对流层延迟改正数据,Riono表示电离层延迟改正数据,Rlt表示负荷潮汐改正数据,Rst表示地球固体潮汐改正数据,Rpt表示极潮改正数据,Rssb表示海况偏差改正数据。
可选的,所述测量标准差获取模块具体包括:
深海压力仪测量海面高值计算单元,用于根据所述海面高差和所述深海压力仪的海深测量值,利用公式SSHsat1=D+dmooring+ΔMSS计算与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值;其中,SSHsat1为与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值,D为所述深海压力仪的零位高度,dmooring为所述深海压力仪的海深测量值,ΔMSS为所述海面高差;
观测差值计算单元,用于根据所述深海压力仪测量海面高值和所述点PCA处卫星测量海面高值,利用公式Δ=SSHsat2-SSHsat1计算单次同步观测期间的观测差值;其中,为单次同步观测期间的观测差值,SSHsat2为点PCA处卫星测量海面高值;
序列标准差值获取单元,用于根据多个同步观测周期中所述观测差值的变化序列,得到序列标准差值;
测高卫星测高精度值获取单元,用于获取所述测高卫星测量海面高的精度值;
深海压力仪测量标准差获取单元,用于根据所述海面高差、所述序列标准差值和所述测高卫星测量海面高的精度值,利用公式计算所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度;其中,为所述深海压力仪的测量标准差,为序列标准差值,为测高卫星测量海面高的精度值,σΔMSS为海面高差的标准差。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明使用高精度测高卫星数据来评估深海压力仪测量海深相对变化的精度,对评定压力仪在深海环境下的测量精度提供了可行的途径;充分利用了卫星测高技术在深海海域高精度海面高观测特性,精度评估结果具有高可靠性;方法原理简单,计算简便,可作为第三方对深海压力仪进行性能检验的工具。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明确定深海压力仪海深测量精度的方法的流程示意图;
图2为本发明确定深海压力仪海深测量精度的系统的结构示意图;
图3为本发明具体实施案例中的系统结构示意图;
图4为本发明具体实施案例中点CP与卫星星下点轨迹位置的关系示意图;
图5为本发明具体实施案例中测站与地面轨迹之间的位置关系示意图;
图6为本发明具体实施案例中观测差值的变化序列图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种确定深海压力仪海深测量精度的方法及系统,利用测高卫星对于海面高高精度测量的特点优势并以其观测值作为深海压力仪评估的基准参考值,从而评估深海压力仪对于相对海深的测量精度,给出深海压力仪在特定海深处的测量指标,作为深海压力仪测量性能的第三方评价手段。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明确定深海压力仪海深测量精度的方法的流程示意图。如图1所示,所述确定深海压力仪海深测量精度的方法包括以下步骤:
步骤100:获取深海压力仪的布设位置。
步骤200:获取测高卫星与深海压力仪同步观测期间的卫星轨道数据和深海压力仪在测高卫星过顶时间的海深测量数据。
根据确定的测高卫星,可以直接获得同步观测期间多个周期的卫星轨道数据,包括卫星的轨道高度、卫星到海面的高度计的测距值、干对流层延迟改正数据、湿对流层延迟改正数据、电离层延迟改正数据、负荷潮汐改正数据、地球固体潮汐改正数据、极潮改正数据和海况偏差改正数据。而且,还可以得到深海压力仪在卫星过顶时间的海深测量数据。
步骤300:根据卫星轨道数据和深海压力仪的布设位置,获取海面高差。所述海面高差为所述深海压力仪布设位置对应海面点CP与卫星地面轨迹中距离点CP最近的点PCA之间的海面高差。具体过程如下:
首先,根据卫星轨道数据以及压力仪布设位置的坐标,可以求得点PCA的位置,以及计算海面点CP与点PCA之间的距离l。
然后,根据压力仪布设位置的坐标,采用全球平均海面高模型计算点CP的平均海面高梯度值GradMSS。
最后,利用公式ΔMSS=GradMSS×l确定所述海面高差ΔMSS。
步骤400:根据卫星轨道数据,采用差值计算,获得点PCA处卫星测量海面高值。具体的,卫星轨道数据根据误差变化特性在计算过程中需要插值计算得到,插值计算方法如表1所示,采用插值计算方法可以得到卫星到达点PCA时的卫星轨道数据的具体取值。
表1不同参数项的拟合插值方法
参数项 | 拟合插值方法 |
电离层 | TCA-21s至+20s求均值 |
干对流层 | TCA-5s至+5s线性拟合;在TCA时刻插值 |
湿对流层 | TCA-15s至+15s之间线性拟合;在TCA时刻插值 |
海况偏差 | TCA-10s至+10s三次多项式拟合;在TCA时刻插值 |
Ku波段测距值 | TCA-20s至+19s之间五阶多项式拟合;在TCA时刻插值 |
卫星轨道高 | TCA-20s至+19s之间三次多项式拟合;在TCA时刻插值 |
深海压力仪测量值 | TCA-1100s至+1100s线性拟合;在TCA时刻插值 |
然后,利用公式SSHsat2=h-(R+Rdry+Rwet+Riono+Rssb+Rlt+Rst+Rpt)计算得到点PCA处卫星测量海面高值;其中,SSHsat为点PCA处卫星测量海面高值,h为卫星的轨道高度,R为卫星到海面的高度计的测距值,Rdry表示干对流层延迟改正数据,Rwet表示湿对流层延迟改正数据,Riono表示电离层延迟改正数据,Rlt表示负荷潮汐改正数据,Rst表示地球固体潮汐改正数据,Rpt表示极潮改正数据,Rssb表示海况偏差改正数据。
步骤500:根据深海压力仪在测高卫星过顶时间的海深测量数据,获取卫星测量海面高值对应时刻的深海压力仪的海深测量值。
将卫星到达点PCA的时刻定义为TCA,则根据深海压力仪在测高卫星过顶时间的海深测量数据,可以筛选得到TCA时刻对应的深海压力仪的海深测量值,即点PCA处卫星测量海面高值对应的深海压力仪的海深测量值。
步骤600:根据海面高差、点PCA处卫星测量海面高值和深海压力仪的海深测量值,获得深海压力仪的测量标准差,得到深海压力仪进行海深测量时的精度。具体过程如下:
Step1:根据所述海面高差和所述深海压力仪的海深测量值,利用公式SSHsat1=D+dmooring+ΔMSS计算与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值;其中,SSHsat1为与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值,D为所述深海压力仪的零位高度,dmooring为所述深海压力仪的海深测量值,ΔMSS为所述海面高差。
Step2:根据所述深海压力仪测量海面高值和所述点PCA处卫星测量海面高值,利用公式计算单次同步观测期间的观测差值;其中,为单次同步观测期间的观测差值,SSHsat2为点PCA处卫星测量海面高值。根据Step1的公式,可以得到
Step3:根据多个同步观测周期中所述观测差值的变化序列,得到序列标准差值
Step4:获取所述测高卫星测量海面高的精度值。此参数与所采用的测高卫星相关,根据所采用的测高卫星,经查阅便可得到测高卫星测量海面高精度
Step5:根据所述海面高差、所述序列标准差值和所述测高卫星测量海面高的精度值,利用公式计算所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度;其中,为所述深海压力仪的测量标准差,为序列标准差值,为测高卫星测量海面高的精度值,σΔMSS为海面高差的标准差。
公式对应的误差传播方程为由于D在零潮汐框架下为常值,因而有σD=0,故可得到因此,
对应于图1所示的确定深海压力仪海深测量精度的方法,本发明还提供一种确定深海压力仪海深测量精度的系统,图2为本发明确定深海压力仪海深测量精度的系统的结构示意图。如图2所示,所述深海压力仪海深测量精度的系统包括以下结构:
布设位置获取模块201,用于获取深海压力仪的布设位置。
观测数据获取模块202,用于获取测高卫星与所述深海压力仪同步观测期间的卫星轨道数据和深海压力仪在所述测高卫星过顶时间的海深测量数据。
海面高差获取模块203,用于根据所述卫星轨道数据和所述深海压力仪的布设位置,获取海面高差;所述海面高差为所述深海压力仪布设位置对应海面点CP与卫星地面轨迹中距离点CP最近的点PCA之间的海面高差。
卫星测量海面高值获取模块204,用于根据所述卫星轨道数据,采用差值计算,获得点PCA处卫星测量海面高值。
深海压力仪的海深测量值获取模块205,用于根据所述深海压力仪在所述测高卫星过顶时间的海深测量数据,获取所述卫星测量海面高值对应时刻的所述深海压力仪的海深测量值。
测量标准差获取模块206,用于根据所述海面高差、所述点PCA处卫星测量海面高值和所述深海压力仪的海深测量值,获得所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度。
作为另一实施例,本发明的所述观测数据获取模块202获取的测高卫星与所述深海压力仪同步观测期间的卫星轨道数据,包括卫星的轨道高度、卫星到海面的高度计的测距值、干对流层延迟改正数据、湿对流层延迟改正数据、电离层延迟改正数据、负荷潮汐改正数据、地球固体潮汐改正数据、极潮改正数据和海况偏差改正数据。
作为另一实施例,本发明的所述海面高差获取模块203具体包括:
距离计算单元,用于根据所述卫星轨道数据和所述深海压力仪的布设位置,计算海面点CP与点PCA之间的距离。
平均海面高梯度值计算单元,用于采用全球平均海面高模型计算点CP的平均海面高梯度值。
海面高差确定单元,用于根据所述海面点CP与点PCA之间的距离和所述点CP的平均海面高梯度值,利用公式ΔMSS=GradMSS×l确定所述海面高差;其中,ΔMSS为海面高差,GradMSS为点CP的平均海面高梯度值,l为海面点CP与点PCA之间的距离。
作为另一实施例,本发明的所述卫星测量海面高值获取模块204利用SSHsat2=h-(R+Rdry+Rwet+Riono+Rssb+Rlt+Rst+Rpt)计算得到点PCA处卫星测量海面高值;其中,SSHsat为点PCA处卫星测量海面高值,h为卫星的轨道高度,R为卫星到海面的高度计的测距值,Rdry表示干对流层延迟改正数据,Rwet表示湿对流层延迟改正数据,Riono表示电离层延迟改正数据,Rlt表示负荷潮汐改正数据,Rst表示地球固体潮汐改正数据,Rpt表示极潮改正数据,Rssb表示海况偏差改正数据。
作为另一实施例,本发明的所述测量标准差获取模块206具体包括:
深海压力仪测量海面高值计算单元,用于根据所述海面高差和所述深海压力仪的海深测量值,利用公式SSHsat1=D+dmooring+ΔMSS计算与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值;其中,SSHsat1为与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值,D为所述深海压力仪的零位高度,dmooring为所述深海压力仪的海深测量值,ΔMSS为所述海面高差。
观测差值计算单元,用于根据所述深海压力仪测量海面高值和所述点PCA处卫星测量海面高值,利用公式计算单次同步观测期间的观测差值;其中,为单次同步观测期间的观测差值,SSHsat2为点PCA处卫星测量海面高值。
序列标准差值获取单元,用于根据多个同步观测周期中所述观测差值的变化序列,得到序列标准差值。
测高卫星测高精度值获取单元,用于获取所述测高卫星测量海面高的精度值。
深海压力仪测量标准差获取单元,用于根据所述海面高差、所述序列标准差值和所述测高卫星测量海面高的精度值,利用公式计算所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度;其中,为所述深海压力仪的测量标准差,为序列标准差值,为测高卫星测量海面高的精度值,σΔMSS为海面高差的标准差。
下面提供一个具体实施案例,进一步说明本发明图1和图2所示的技术方案。
图3为本发明具体实施案例中的系统结构示意图,如图3所示,本实施案例采用卫星测高技术通过卫星搭载的高度计,辅以精密定轨技术,再进行高度计电磁信号的传播延迟改正、海况偏差效应改正等,可测定卫星星下点的瞬时海面高SSH。压力仪所测定的深度变化为瞬时海面高与压力仪零位高度之差。压力仪零位高度无法精确直接测定,在零潮汐框架下压力仪零位高度可认为是常值。
瞬时海面高与压力仪对应的海面位置(记为点CP)通常不完全重合,记卫星地面轨迹最接近压力仪布设位置最接近的点为PCA(Point of Closest Approach)点,它们的位置关系如图4所示,图4为本发明具体实施案例中点CP与卫星星下点轨迹位置的关系示意图,其中T为卫星地面轨迹运行方向。与点PCA相对应的时刻记为TCA(Time of ClosestApproach)。
点PCA与点CP之间的实际海面高差为瞬时海面高差ΔSSH,瞬时海面高差可通过使用平均海面高(Mean Sea Surface,MSS)之差来近似,平均海面高差记为ΔMSS,ΔMSS=MSSPCA-MSSCP。ΔMSS可以通过点CP的平均海面高梯度GradMSS以及PCA、CP两点之间的距离l之积求得,即:
ΔMSS=MSSPCA-MSSDART=GradMSS×l (1)
其中,GradMSS可通过全球平均海面高模型(例如CNES_CLS15)求得。
图4中,y轴与北向一致,轨迹法向与y轴正向之间的夹角为θ([-π,π]),以y轴绕原点顺时针方向为正。首先由全球平均海面高模型中点CP邻近网格点的平均海面高,求取点CP东西方向平均海面高梯度Gradx和南北方向平均海面高梯度Grady,再根据式(2)可求取GradMSS。
GradMSS=Gradxsinθ+Gradycosθ (2)
与点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值为:
SSHsat1=D+dmooring+ΔMSS (3)
零潮汐条件下的SSHsat2可以由式(4)计算得到,即:
SSHsat2=h-(R+Rdry+Rwet+Riono+Rssb+Rlt+Rst+Rpt) (4)
其中h为卫星的轨道高度,R为卫星到海面的高度计的测距值,Rdry表示干对流层延迟改正数据,Rwet表示湿对流层延迟改正数据,Riono表示电离层延迟改正数据,Rlt表示负荷潮汐改正数据,Rst表示地球固体潮汐改正数据,Rpt表示极潮改正数据,Rssb表示海况偏差改正数据。
本具体实施案例将深海压力仪布设于数千米的海底,这里以Jason-3测高卫星与DART测站对临近海域同步观测的数据为例,给出本发明的计算结果。21419号DART测站布设区域海深约5.3km,其布设位置与Jason-3卫星理论地面轨迹的距离小于1km,两者之间的位置关系如图5所示,图5为本发明具体实施案例中测站与地面轨迹之间的位置关系示意图,该DART测站选用Paroscientific Digiquartz 410K型压力传感器。具体过程如下:
步骤一:根据21419号DART测站的观测数据,该测站在2016年2月23日至2017年7月3日期间布设坐标位于(44.435°N,155.717°E),由该坐标,通过DTU13、CNES_CLS15、WHU2013三种全球平均海面高模型计算得到该区域布设点处卫星轨迹法向平均海面高梯度值GradMSS为西向东-1.43cm/km。
以Jason-3卫星第3周期的同步观测数据为例,步骤二至步骤五如下:
步骤二:根据Jason-3卫星的GDR数据产品,获取卫星过顶压力仪前后九项参数的变化数值(TCA前后30s);根据卫星过顶时刻,获取TCA前后1h压力仪测量海深变化数值。
步骤三:求得PCA点位置为(44.4396°N,155.7065°E),PCA点与压力仪对应海面点之间的距离为0.98km;根据步骤一计算的平均海面高梯度值GradMSS,求得两点之间ΔMSS为-1.4cm(PCA点高于压力仪点海面1.4cm)。
步骤四:根据表1,求得TCA时刻PCA点卫星测量海面高SSHsat为10.8234m,dmooring为10.9594m(设D为常数5336.586m);
步骤五:根据式单次同步观测期间观测差值为-12.2cm;
通过2016年2月23日至2017年7月3日期间,Jason-3卫星与21419号DART压力仪51次同步观测的结果,重复步骤二至步骤五,得到差值序列,如图6所示,图6为本发明具体实施案例中观测差值的变化序列图。
步骤六:计算得到差值序列标准差为10.78cm;
步骤七:查阅Jason-3卫星用户手册,该卫星海面高测量精度为3.53cm,σΔMSS设为0.3cm。通过式计算得到压力仪测量标准差为10.18cm,即21419号DART测站所使用的深海压力仪海深相对测量精度为10.18cm。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种确定深海压力仪海深测量精度的方法,其特征在于,包括:
获取深海压力仪的布设位置;
获取测高卫星与所述深海压力仪同步观测期间的卫星轨道数据和深海压力仪在所述测高卫星过顶时间的海深测量数据;
根据所述卫星轨道数据和所述深海压力仪的布设位置,获取海面高差;所述海面高差为所述深海压力仪布设位置对应海面点CP与卫星地面轨迹中距离点CP最近的点PCA之间的海面高差;
根据所述卫星轨道数据,采用差值计算,获得点PCA处卫星测量海面高值;
根据所述深海压力仪在所述测高卫星过顶时间的海深测量数据,获取所述卫星测量海面高值对应时刻的所述深海压力仪的海深测量值;
根据所述海面高差、所述点PCA处卫星测量海面高值和所述深海压力仪的海深测量值,获得所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度。
2.根据权利要求1所述的确定深海压力仪海深测量精度的方法,其特征在于,所述测高卫星与所述深海压力仪同步观测期间的卫星轨道数据包括卫星的轨道高度、卫星到海面的高度计的测距值、干对流层延迟改正数据、湿对流层延迟改正数据、电离层延迟改正数据、负荷潮汐改正数据、地球固体潮汐改正数据、极潮改正数据和海况偏差改正数据。
3.根据权利要求1所述的确定深海压力仪海深测量精度的方法,其特征在于,所述根据所述卫星轨道数据和所述深海压力仪的布设位置,获取海面高差,具体包括:
根据所述卫星轨道数据和所述深海压力仪的布设位置,计算海面点CP与点PCA之间的距离;
采用全球平均海面高模型计算点CP的平均海面高梯度值;
根据所述海面点CP与点PCA之间的距离和所述点CP的平均海面高梯度值,利用公式ΔMSS=GradMSS×l确定所述海面高差;其中,ΔMSS为海面高差,GradMSS为点CP的平均海面高梯度值,l为海面点CP与点PCA之间的距离。
4.根据权利要求2所述的确定深海压力仪海深测量精度的方法,其特征在于,所述根据所述卫星轨道数据,采用差值计算,获得点PCA处卫星测量海面高值,具体包括:
利用SSHsat2=h-(R+Rdry+Rwet+Riono+Rssb+Rlt+Rst+Rpt)计算得到点PCA处卫星测量海面高值;其中,SSHsat为点PCA处卫星测量海面高值,h为卫星的轨道高度,R为卫星到海面的高度计的测距值,Rdry表示干对流层延迟改正数据,Rwet表示湿对流层延迟改正数据,Riono表示电离层延迟改正数据,Rlt表示负荷潮汐改正数据,Rst表示地球固体潮汐改正数据,Rpt表示极潮改正数据,Rssb表示海况偏差改正数据。
5.根据权利要求1所述的确定深海压力仪海深测量精度的方法,其特征在于,所述根据所述海面高差、所述点PCA处卫星测量海面高值和所述深海压力仪的海深测量值,获得所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度,具体包括:
根据所述海面高差和所述深海压力仪的海深测量值,利用公式SSHsat1=D+dmooring+ΔMSS计算与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值;其中,SSHsat1为与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值,D为所述深海压力仪的零位高度,dmooring为所述深海压力仪的海深测量值,ΔMSS为所述海面高差;
根据所述深海压力仪测量海面高值和所述点PCA处卫星测量海面高值,利用公式计算单次同步观测期间的观测差值;其中,为单次同步观测期间的观测差值,SSHsat2为点PCA处卫星测量海面高值;
根据多个同步观测周期中所述观测差值的变化序列,得到序列标准差值;
获取所述测高卫星测量海面高的精度值;
根据所述海面高差、所述序列标准差值和所述测高卫星测量海面高的精度值,利用公式计算所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度;其中,为所述深海压力仪的测量标准差,为序列标准差值,为测高卫星测量海面高的精度值,σΔMSS为海面高差的标准差。
6.一种确定深海压力仪海深测量精度的系统,其特征在于,包括:
布设位置获取模块,用于获取深海压力仪的布设位置;
观测数据获取模块,用于获取测高卫星与所述深海压力仪同步观测期间的卫星轨道数据和深海压力仪在所述测高卫星过顶时间的海深测量数据;
海面高差获取模块,用于根据所述卫星轨道数据和所述深海压力仪的布设位置,获取海面高差;所述海面高差为所述深海压力仪布设位置对应海面点CP与卫星地面轨迹中距离点CP最近的点PCA之间的海面高差;
卫星测量海面高值获取模块,用于根据所述卫星轨道数据,采用差值计算,获得点PCA处卫星测量海面高值;
深海压力仪的海深测量值获取模块,用于根据所述深海压力仪在所述测高卫星过顶时间的海深测量数据,获取所述卫星测量海面高值对应时刻的所述深海压力仪的海深测量值;
测量标准差获取模块,用于根据所述海面高差、所述点PCA处卫星测量海面高值和所述深海压力仪的海深测量值,获得所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度。
7.根据权利要求6所述的确定深海压力仪海深测量精度的系统,其特征在于,所述观测数据获取模块获取的测高卫星与所述深海压力仪同步观测期间的卫星轨道数据,包括卫星的轨道高度、卫星到海面的高度计的测距值、干对流层延迟改正数据、湿对流层延迟改正数据、电离层延迟改正数据、负荷潮汐改正数据、地球固体潮汐改正数据、极潮改正数据和海况偏差改正数据。
8.根据权利要求6所述的确定深海压力仪海深测量精度的系统,其特征在于,所述海面高差获取模块具体包括:
距离计算单元,用于根据所述卫星轨道数据和所述深海压力仪的布设位置,计算海面点CP与点PCA之间的距离;
平均海面高梯度值计算单元,用于采用全球平均海面高模型计算点CP的平均海面高梯度值;
海面高差确定单元,用于根据所述海面点CP与点PCA之间的距离和所述点CP的平均海面高梯度值,利用公式ΔMSS=GradMSS×l确定所述海面高差;其中,ΔMSS为海面高差,GradMSS为点CP的平均海面高梯度值,l为海面点CP与点PCA之间的距离。
9.根据权利要求7所述的确定深海压力仪海深测量精度的系统,其特征在于,所述卫星测量海面高值获取模块利用SSHsat2=h-(R+Rdry+Rwet+Riono+Rssb+Rlt+Rst+Rpt)计算得到点PCA处卫星测量海面高值;其中,SSHsat为点PCA处卫星测量海面高值,h为卫星的轨道高度,R为卫星到海面的高度计的测距值,Rdry表示干对流层延迟改正数据,Rwet表示湿对流层延迟改正数据,Riono表示电离层延迟改正数据,Rlt表示负荷潮汐改正数据,Rst表示地球固体潮汐改正数据,Rpt表示极潮改正数据,Rssb表示海况偏差改正数据。
10.根据权利要求6所述的确定深海压力仪海深测量精度的系统,其特征在于,所述测量标准差获取模块具体包括:
深海压力仪测量海面高值计算单元,用于根据所述海面高差和所述深海压力仪的海深测量值,利用公式SSHsat1=D+dmooring+ΔMSS计算与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值;其中,SSHsat1为与所述点PCA处卫星测量海面高值相同时刻的深海压力仪测量海面高值,D为所述深海压力仪的零位高度,dmooring为所述深海压力仪的海深测量值,ΔMSS为所述海面高差;
观测差值计算单元,用于根据所述深海压力仪测量海面高值和所述点PCA处卫星测量海面高值,利用公式计算单次同步观测期间的观测差值;其中,为单次同步观测期间的观测差值,SSHsat2为点PCA处卫星测量海面高值;
序列标准差值获取单元,用于根据多个同步观测周期中所述观测差值的变化序列,得到序列标准差值;
测高卫星测高精度值获取单元,用于获取所述测高卫星测量海面高的精度值;
深海压力仪测量标准差获取单元,用于根据所述海面高差、所述序列标准差值和所述测高卫星测量海面高的精度值,利用公式计算所述深海压力仪的测量标准差,得到所述深海压力仪进行海深测量时的精度;其中,为所述深海压力仪的测量标准差,为序列标准差值,为测高卫星测量海面高的精度值,σΔMSS为海面高差的标准差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910999426.5A CN110553609B (zh) | 2019-10-21 | 2019-10-21 | 一种确定深海压力仪海深测量精度的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910999426.5A CN110553609B (zh) | 2019-10-21 | 2019-10-21 | 一种确定深海压力仪海深测量精度的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110553609A true CN110553609A (zh) | 2019-12-10 |
CN110553609B CN110553609B (zh) | 2021-07-23 |
Family
ID=68743140
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910999426.5A Active CN110553609B (zh) | 2019-10-21 | 2019-10-21 | 一种确定深海压力仪海深测量精度的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110553609B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111650608A (zh) * | 2020-06-15 | 2020-09-11 | 中国人民解放军61540部队 | 一种Ka/C双频高度计消除电离层误差方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5191719A (en) * | 1991-10-21 | 1993-03-09 | Kitt Steve B | Depth gauge for a fishing line |
CN2648433Y (zh) * | 2003-08-18 | 2004-10-13 | 浙江大学 | 一种深海热液气体的取样分析装置 |
CN101231167A (zh) * | 2008-02-20 | 2008-07-30 | 刘雁春 | 海洋测线网系统误差的检测与调整方法 |
CN102997856A (zh) * | 2012-12-12 | 2013-03-27 | 南京大学 | 一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法 |
CN103364766A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院电子学研究所 | 星载InSAR系统的外定标方法 |
CN104678354A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于gnss-r的海面目标双星定位方法 |
CN206459611U (zh) * | 2017-01-17 | 2017-09-01 | 中建一局集团建设发展有限公司 | 一种基于北斗卫星的超高层建筑高精度变形监测系统 |
CN208313246U (zh) * | 2018-05-18 | 2019-01-01 | 广州墨匠科技有限公司 | 一种液压式测量仪 |
CN109212524A (zh) * | 2018-10-24 | 2019-01-15 | 航天东方红卫星有限公司 | 一种考虑校正量误差相关性的卫星测高精度确定方法 |
-
2019
- 2019-10-21 CN CN201910999426.5A patent/CN110553609B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5191719A (en) * | 1991-10-21 | 1993-03-09 | Kitt Steve B | Depth gauge for a fishing line |
CN2648433Y (zh) * | 2003-08-18 | 2004-10-13 | 浙江大学 | 一种深海热液气体的取样分析装置 |
CN101231167A (zh) * | 2008-02-20 | 2008-07-30 | 刘雁春 | 海洋测线网系统误差的检测与调整方法 |
CN103364766A (zh) * | 2012-04-01 | 2013-10-23 | 中国科学院电子学研究所 | 星载InSAR系统的外定标方法 |
CN102997856A (zh) * | 2012-12-12 | 2013-03-27 | 南京大学 | 一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法 |
CN104678354A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种基于gnss-r的海面目标双星定位方法 |
CN206459611U (zh) * | 2017-01-17 | 2017-09-01 | 中建一局集团建设发展有限公司 | 一种基于北斗卫星的超高层建筑高精度变形监测系统 |
CN208313246U (zh) * | 2018-05-18 | 2019-01-01 | 广州墨匠科技有限公司 | 一种液压式测量仪 |
CN109212524A (zh) * | 2018-10-24 | 2019-01-15 | 航天东方红卫星有限公司 | 一种考虑校正量误差相关性的卫星测高精度确定方法 |
Non-Patent Citations (3)
Title |
---|
ANDY STOCK: "Satellite mapping of Baltic Sea Secchi depth with multiple regression models", 《INTERNATIONAL JOURNAL OF APPLIED EARTH OBSERVATION AND GEOINFORMATION》 * |
欧阳明达: "利用测高重力异常建立南中国海海底地形模型", 《大地测量与地球动力学》 * |
聂琳娟: "卫星激光测高ICESat 确定海面高精度分析", 《测绘信息与工程》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111650608A (zh) * | 2020-06-15 | 2020-09-11 | 中国人民解放军61540部队 | 一种Ka/C双频高度计消除电离层误差方法及系统 |
CN111650608B (zh) * | 2020-06-15 | 2023-01-10 | 中国人民解放军61540部队 | 一种Ka/C双频高度计消除电离层误差方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110553609B (zh) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106970404B (zh) | 一种基于Delaunay三角网的多冗余网络RTK大气误差内插方法 | |
Kaimal et al. | Some errors in the measurement of Reynolds stress | |
CN100535683C (zh) | 用于超短基线声学定位系统的校准方法 | |
CN109085655B (zh) | 一种水下平台重力测量方案与验证方法 | |
CN107479078A (zh) | 铁路测量中大地坐标转换为独立平面坐标的方法及系统 | |
CN106443722A (zh) | 一种天线相位中心偏差检测方法 | |
CN109856691B (zh) | 一种基于梯度法的航空重力矢量向下延拓方法及系统 | |
CN108254766B (zh) | 一种卫星抗多径误差方法 | |
CN108413983A (zh) | 一种sins/usbl一体化定位系统安装误差标定的机械方法 | |
CN110553609B (zh) | 一种确定深海压力仪海深测量精度的方法及系统 | |
CN114383578B (zh) | 一种基于海洋监测浮标的海面高测量系统及方法 | |
CN110082033B (zh) | 一种运动状态下的水上载体重心测量装置和方法 | |
CN116858290B (zh) | 基于大型无人机的深远海海面高度观测定标方法和系统 | |
CN105204047A (zh) | 一种卫星导航系统中观测量单个粗差的探测与修复方法 | |
CN109254270A (zh) | 一种星载x波段合成孔径雷达干涉定标方法 | |
CN110850382B (zh) | 一种评估干涉雷达高度计测量精度的方法及系统 | |
CN105352459A (zh) | 一种测边前方交会测点的方法 | |
CN111412930B (zh) | 一种联合测姿装置安装误差的标定运算方法 | |
JP2014048075A (ja) | Gnssによる位置計測装置 | |
CN113819882A (zh) | 一种跨海高程点间重力位差计算方法 | |
CN113819892A (zh) | 基于半参数估计和附加深度约束的深海基准网平差方法 | |
CN114063181B (zh) | 一种高精度海洋重力测量重力基点读数校正方法 | |
CN114125430B (zh) | 一种光电转塔上摄像机视轴偏离度的校准方法 | |
CN116753991B (zh) | 一种基于固定场区的卫星高度计定标方法及系统 | |
Alves et al. | Introduction of a geometry-based Network RTK quality indicator |
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 |