CN115453587B - 顾及潮族响应关系的gnss海潮负荷位移反演方法及系统 - Google Patents
顾及潮族响应关系的gnss海潮负荷位移反演方法及系统 Download PDFInfo
- Publication number
- CN115453587B CN115453587B CN202210926138.9A CN202210926138A CN115453587B CN 115453587 B CN115453587 B CN 115453587B CN 202210926138 A CN202210926138 A CN 202210926138A CN 115453587 B CN115453587 B CN 115453587B
- Authority
- CN
- China
- Prior art keywords
- tide
- gnss
- family
- load displacement
- main
- 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
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000004044 response Effects 0.000 title claims abstract description 51
- 238000004458 analytical method Methods 0.000 claims abstract description 41
- 238000012937 correction Methods 0.000 claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000003089 Pariser Parr Pople method Methods 0.000 claims description 9
- 238000012935 Averaging Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 230000007774 longterm Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000018199 S phase Effects 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000007787 solid Substances 0.000 description 4
- 230000009897 systematic effect Effects 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 208000025174 PANDAS Diseases 0.000 description 1
- 208000021155 Paediatric autoimmune neuropsychiatric disorders associated with streptococcal infection Diseases 0.000 description 1
- 240000000220 Panda oleosa Species 0.000 description 1
- 235000016496 Panda oleosa Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000003797 telogen phase Effects 0.000 description 1
- 238000004613 tight binding model Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/40—Correcting position, velocity or attitude
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/393—Trajectory determination or predictive tracking, e.g. Kalman filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Probability & Statistics with Applications (AREA)
- Operations Research (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种顾及潮族响应关系的GNSS海潮负荷位移反演方法及系统，包括GNSS动态PPP数据采集与处理，得到潮汐调和分析的输入数据；明确周日半周日潮族响应关系，直接对周日半周日主分潮导纳在频率域进行线性插值而获取小分潮的导纳；构建顾及潮族响应关系的调和分析模型；求解主分潮的待估参数，进而获取调和常数；利用所求周日半周日主分潮的调和常数联合海潮模型估计长周期分潮的调和常数，在进行PPP动态解算时，进行海潮负荷位移修正，进一步提高GNSS定位精度。本发明利用潮族响应关系进行潮汐调和分析，不仅考虑同潮族其余小分潮影响，还避免了现有交点改正方法所用平衡潮理论值与实际潮汐成分不一致而导致的误差，从而提高GNSS海潮负荷位移反演精度。
Description
技术领域
本发明属于GNSS高精度定位及其应用领域，在利用高精度GNSS观测值提取海潮负荷信号领域中应用，具体涉及一种顾及潮族响应关系的GNSS海潮负荷位移反演方案。
背景技术
在日、月引潮力的作用下，固定地球会产生周期性形变，这一现象称为固体潮，实际海平面也会发生周期性的涨落，即海潮。海潮引起的海水质量重新分布同样会导致固体地球发生弹性形变，称之为海潮负荷效应。低纬度地区测站因固体潮引起的位移变化可高达几十厘米，可根据IERS协议2010对其进行有效改正，精度已达到1mm量级。沿海区域测站海潮负荷效应可达十几厘米，可由海潮模型和格林函数通过离散褶积积分计算改正，但海潮模型在浅海海域与高纬度海域精度欠佳，导致海潮负荷效应改正无法达到mm级精度的需求。因此，如何提高海潮负荷位移计算精度已成为高精度大地测量亟待解决的关键问题。
随着GNSS技术的不断发展和完善，利用GNSS观测值可以准确反演海潮负荷位移参数已被证实。在静态PPP求解单天解测站坐标等参数时，将8个周日半周日主分潮三维方向的调和常数也作为待估参数一同估计，存在待估参数过多降低解算强度，采样间隔24h容易混叠高频海潮信号等弊端。动态PPP可以逐历元(一般为30s)获取海潮负荷引起的测站位移变化，再通过调和分析估算主要分潮的调和常数。相比于静态估计法，动态估计法能够更加准确地提取海潮负荷信号。
采用动态PPP获取的测站坐标时间序列的时间跨度一般小于18.61年，在潮汐调和分析中，为了削弱或消除临近频率同一亚群小分潮的影响，通常采用交点改正获取更加准确的8个周日半周日主分潮的调和常数。实际上，中国、欧洲南部沿海等地区观测到的潮汐交点改正与其相应的平衡潮理论值存在一定的差异，可能受浅水区域海底摩擦、辐射潮等因素影响。由此可见，交点改正进行调和分析时在某些特定区域会产生系统性误差，从而直接导致估算的分潮调和常数准确度不足。同时，现有的调和分析方法因GNSS观测时段仅数年甚至更短，也无法直接估算同潮族其余小分潮的调和常数。IERS协议2010推荐的海潮负荷改正方法已考虑周日与半周日潮族共255个、长周期潮族76个小分潮的影响，相关研究表明小分潮会引起沿海或岛屿测站垂向位移变化可达8mm，是否考虑小分潮可能会影响调和分析建模精度。
发明内容
为了能够充分考虑小分潮的影响，并且避免交点改正方法在部分区域存在系统性误差的问题，进一步提高GNSS海潮负荷位移反演精度，本发明提出了一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
为了实现本发明的目的，本发明采用的技术方案是一种顾及潮族响应关系的GNSS海潮负荷位移反演方法，包括以下步骤：
步骤1，GNSS动态PPP数据采集与处理，包括采集GNSS连续观测站的长期观测数据，利用精密轨道和钟差产品进行动态PPP解算，解算时未采用海潮负荷位移改正，获取逐历元地心三维坐标时间序列，取平均得到时段解后转化为站心地平坐标时间序列，进行粗差探测与剔除后作为潮汐调和分析的输入数据；
步骤2，明确周日半周日潮族响应关系，包括直接对周日半周日主分潮导纳在频率域进行线性插值而获取小分潮的导纳；
步骤3，构建顾及潮族响应关系的调和分析模型，所述顾及潮族响应关系的调和分析模型针对周日半周日主分潮进行调和分析时，在考虑交点改正的影响同时，考虑同潮族其余小分潮的影响；计算出小分潮对主分潮待估参数的累积影响系数，优化观测方程的系数矩阵；
步骤4，利用最小二乘方法求解主分潮的待估参数，进而获取振幅和格林尼治相位，作为调和常数；
步骤5，利用所求周日半周日主分潮的调和常数联合海潮模型估计长周期分潮的调和常数，在进行PPP动态解算时，进行海潮负荷位移修正，进一步提高GNSS定位精度。
而且，步骤1中，将逐历元解取平均得到时段解，取平均得到时段解时，时段为5、10、15、20或30分钟。
而且，步骤2中，忽略插值方法异同对海潮负荷位移预测结果影响，将三次样条插值方法简化为以下线性关系，
其中，满足f2>fx>f1，表示小分潮频率fx位于同潮族相邻两个主分潮f1、f2之间，u1和u2为对应主分潮导纳，ux为小分潮导纳(实部或虚部)。
而且，步骤3中，建立求解主分潮调和常数的观测方程如下，
其中，Z0为常数项，a为线性趋势项，n表示主分潮个数，参数为第i个分潮对应的Cartwright-Edden振幅/>的绝对值，Vi'(t)为第i个分潮t时刻的原始天文相角，分别对周日潮族各个分潮的原始天文相角作加90°处理，αj,i、βj,i为小分潮对第i个主分潮在j方向待估的参数/>的累积影响系数，R(t)为残差；
参数uc、us分别为主要周日半周日分潮导纳u的实部和虚部。
另一方面，本发明提供一种顾及潮族响应关系的GNSS海潮负荷位移反演系统，用于实现如上所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
而且，包括以下模块，
第一模块，用于GNSS动态PPP数据采集与处理，包括采集GNSS连续观测站的长期观测数据，利用精密轨道和钟差产品进行动态PPP解算，解算时未采用海潮负荷位移改正，获取逐历元地心三维坐标时间序列，取平均得到时段解后转化为站心地平坐标时间序列，进行粗差探测与剔除，为潮汐调和分析的输入数据；
第二模块，用于明确周日半周日潮族响应关系，包括直接对周日半周日主分潮导纳在频率域进行线性插值而获取小分潮的导纳；
第三模块，用于构建顾及潮族响应关系的调和分析模型，所述顾及潮族响应关系的调和分析模型针对周日半周日主分潮进行调和分析时，在考虑交点改正的影响同时，考虑同潮族其余小分潮的影响；计算出小分潮对主分潮待估参数的累积影响系数，优化观测方程的系数矩阵；
第四模块，用于利用最小二乘方法求解主分潮的待估参数，进而获取振幅和格林尼治相位，作为调和常数；
第五模块，用于利用所求周日半周日主分潮的调和常数联合海潮模型估计长周期分潮的调和常数，在进行PPP动态解算时，进行海潮负荷位移修正，进一步提高GNSS定位精度。
或者，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
或者，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
与现有的技术相比，本发明具有如下有益效果：
本发明所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方案，构建调和分析模型时，创新性发现同潮族主分潮与小分潮之间可通过简化线性插值替代三次样条插值，不仅考虑了交点改正的影响，还充分考虑了同潮族其余小分潮的影响，使观测方程建模更加准确，潮汐残余信号大大减少。本发明考虑潮族响应关系而非交点改正方法进行调和分析，有效避免了传统方法因实际观测到的交点改正分量与其对应的平衡潮理论值不一致而引起的系统性误差问题，提高了GNSS海潮负荷位移的反演精度。
附图说明
图1为本发明实施例的方法流程图；
图2为本发明实施例ZJWZ台站预测的三维方向海潮负荷位移因采用三次样条插值与线性插值方法导致的差值序列。
具体实施方式
为了更加清楚地说明本发明的目的、技术方案与有益效果，下面结合附图与实施例对本发明具体实施方式做进一步的详细说明。
参见图1，本发明实施例提供了一种顾及潮族响应关系的GNSS海潮负荷位移反演方法，包括以下步骤：
步骤1：GNSS动态PPP数据采集与处理
具体实施时，取平均得到时段解时，时段可为5、10、15、20或30分钟。实施例中，优选采集GNSS连续观测站的长期观测数据，利用IGS分析中心提供的精密轨道、钟差等产品进行动态PPP解算(未采用海潮负荷位移改正)，获取逐历元地心三维坐标时间序列，取平均为半小时解后转化为站心地平坐标时间序列，进行粗差探测与剔除，作为潮汐调和分析的输入数据。
本实施例针对陆态网ZJWZ台站(位于浙江温州)2019年(001-365)的观测数据，采用PRIDE PPP-AR II软件动态PPP模式求解GPS/GLONASS组合系统模糊度固定解，精密轨道、卫星钟差、地球自转参数、相位偏差产品、天线相位中心改正等采用CODE产品。在解算过程中，固体潮、极潮已按照IERS协议2010进行改正，尚未考虑海潮负荷位移改正，方便获取海潮负荷位移信号。将解算的逐历元地心三维坐标时间序列，取平均为半小时解后转化为站心地平坐标时间序列，进行粗差探测与剔除后作为潮汐调和分析的输入数据。
其中步骤1中采用不同的GNSS数据处理软件(如BERNESE、GIPSY、PANDA等)与不同机构提供的精密产品(如GFZ、JPL、WHU等)均会不可避免地引入不同的解算误差。在现有软件逐步完善，解算策略愈发合理的情况下，GNSS观测值能够较为准确地反映海潮负荷位移信号。此外，单GPS观测值容易受轨道周期与重访周期引起的卫星轨道误差影响，导致K1和K2分潮反演精度较差，本发明实施例采用的GPS/GLONASS组合观测值能够有效削弱这一误差，同时模糊度固定解相比于浮点解也可进一步提高解算精度。
步骤2：明确周日半周日潮族响应关系
本发明进一步提出，可以直接对周日半周日主分潮导纳在频率域进行线性插值而获取小分潮的导纳。
IERS协议2010推荐的海潮负荷位移改正方法，精度可达1％，分别对8个主要周日半周日分潮(M2、S2、N2、K2、K1、O1、P1和Q1)构建导纳u，分为实部和虚部Aj,i和φj,i为第i个分潮在j方向的调和常数，即振幅和格林尼治相位，为第i个分潮对应的Cartwright-Edden振幅/>的绝对值，通过三次样条插值可获取同潮族其余255个小分潮导纳的实部或虚部：
式中，满足f2>fx>f1，表示小分潮频率fx位于同潮族相邻两个主分潮f1、f2之间，s1和s2则为对应主分潮导纳u1和u2(实部或虚部)的二阶导数，ux为小分潮导纳的实部或虚部，当小分潮频率fx小于同潮族最小频率的主分潮fmin时，即fx<fmin，ux＝umin，当小分潮频率大于同潮族最大频率的主分潮fmax时，即fx>fmax，ux＝umax，其中，fmin和fmax对应的分潮导纳的实部或虚部分别为umin和umax。实际上，只有108个小分潮的导纳是通过主分潮三次样条插值获取的，其余直接等于频率最接近的主分潮导纳。同时，式(1)中二阶导数s1和s2对海潮负荷位移预测结果影响很小，若忽略其影响则为简化线性关系：
其中步骤2实施例根据NAO.99b模型估算的ZJWZ台站的海潮负荷位移参数，分别采用三次样条插值与线性插值方法对该台站进行海潮负荷位移预测，图2显示两方法结果在U、W、S方向的海潮负荷位移差值序列，水平方向差值控制在0.06mm以内，垂直方向在0.2mm以内，说明插值方法不同引起的海潮负荷位移差异很小，在后续调和分析建模中可以被残差吸收，验证了步骤2算法利用简化线性插值替代三次样条插值的可行性。由公式(2)可以直接对8个周日半周日主分潮导纳在频率域进行线性插值而获取255个小分潮的导纳，为实现顾及潮族响应关系的调和分析模型提供重要基础。
步骤3：构建顾及潮族响应关系的调和分析模型
本发明提出的顾及潮族响应关系的调和分析模型，针对8个周日半周日主分潮进行调和分析时，考虑了同潮族其余255个小分潮的影响，该模型不仅考虑了交点改正的影响，而且考虑了同潮族其余小分潮的影响，与常用的调和分析模型仅考虑交点改正相比，进一步提高了建模精度。
由步骤2可知，255个小分潮的导纳可由8个主分潮的导纳线性插值获取，而各个分潮的导纳是由振幅和格林尼治相位组合构成(除以Cartwright-Edden振幅的绝对值进行标准化)。因此，可以计算出255个小分潮对8个主分潮待估参数的累积影响系数，优化了观测方程的系数矩阵，待估参数仍不变。
针对t时刻GNSS在j方向的观测值hj(t)进行调和分析：
式中，参数即当/>为正数时，δi＝1，为负数时，δi＝-1，Z0为常数项，a为线性趋势项(GNSS观测值水平方向较为明显)，n表示主分潮个数(n＝8)，m则为小分潮个数(m＝255)，Vi'(t)为第i个分潮t时刻的原始天文相角(未对引潮力位展开式中余弦与正弦函数进行统一)，需要注意的是，分别对周日潮族各个分潮的原始天文相角作加90°处理，R(t)为残差(包含非海潮信号和噪声)，将式(3)线性化可转化为式(4)：
将式(2)代入式(5)，即可将周日半周日潮族255个小分潮的参数或/>由同潮族频率相近的一个或两个主分潮的该参数表示，计算出255个小分潮对第i个主分潮在j方向待估参数/>的累积影响系数αj,i，βj,i，将公式(4)转化为求解8个主分潮调和常数的观测方程：
步骤4：利用最小二乘方法求解8个主分潮的待估参数然后通过下式获取振幅和格林尼治相位，即获取调和常数Aj,i和φj,i：
通过上述步骤，ZJWZ台站的三维海潮负荷位移参数解算结果列于表1，注：下文表格中振幅单位均为m，相位单位为度。
表1ZJWZ台站海潮负荷位移参数的解算结果
| 分潮名称 | M2 | S2 | N2 | K2 | K1 | O1 | P1 | Q1 |
| U振幅 | 0.01740 | 0.00600 | 0.00394 | 0.00296 | 0.00917 | 0.00674 | 0.00353 | 0.00133 |
| W振幅 | 0.00469 | 0.00204 | 0.00105 | 0.00051 | 0.00311 | 0.00250 | 0.00087 | 0.00024 |
| S振幅 | 0.00285 | 0.00158 | 0.00059 | 0.00058 | 0.00238 | 0.00119 | 0.00125 | 0.00036 |
| U相位 | -170.7 | -123.8 | 176.8 | 93.8 | -74.9 | -105.2 | -71.2 | -121.7 |
| W相位 | 164.2 | 170.3 | 153.6 | -138.4 | -147.8 | -134.8 | -106.9 | -149.4 |
| S相位 | 16.5 | 26.0 | -1.3 | -134.2 | 131.7 | 121.3 | 150.7 | 94.5 |
将NAO.99b模型估算的海潮负荷位移参数作为参考，在中国沿海区域，NAO.99b模型的精度通常高于其他全球海潮模型，主要是因为NAO.99b模型同化了日本、韩国沿岸验潮站和5年沿轨海面高数据。与NAO.99b模型估算结果的差值统计列于表2。
表2ZJWZ台站海潮负荷位移参数解算结果与NAO.99b模型估值的差值
| 分潮名称 | M2 | S2 | N2 | K2 | K1 | O1 | P1 | Q1 |
| U振幅 | -0.00097 | 0.00019 | 0.00020 | 0.00141 | 0.00017 | -0.00061 | 0.00063 | -0.00022 |
| W振幅 | -0.00072 | 0.00025 | -0.00007 | 0.00003 | 0.00026 | 0.00046 | -0.00004 | -0.00016 |
| S振幅 | -0.00036 | 0.00064 | -0.00004 | 0.00033 | 0.00038 | -0.00050 | 0.00060 | -0.00001 |
| U相位 | -5.2 | 14.3 | 3.3 | -122.9 | 10.5 | 0.1 | 16.3 | -8.1 |
| W相位 | 0.3 | -16.0 | 5.6 | 41.0 | -36.2 | -0.6 | 6.2 | -3.8 |
| S相位 | -6.8 | -10.8 | -1.5 | -164.8 | -7.6 | 8.9 | 12.2 | -6.8 |
由表2可知本发明所公开的方法可以有效提取海潮负荷位移参数，除K2分潮U方向振幅差值达1.4mm，其余振幅差值均控制在1mm以内，除K2分潮三方向、K1分潮W方向外，其余相位亦符合较好。K2与K1分潮分别与GPS卫星轨道周期和重访周期相近，容易受卫星轨道误差和多路径误差影响，因此K2、K1分潮精度不及其他分潮。同时利用传统方法即顾及交点改正进行调和分析，其结果与模型估值的差值列于表3。
表3传统方法ZJWZ台站海潮负荷位移参数解算结果与NAO.99b模型估值的差值
对比表2和表3可知，与NAO.99b模型估值相比，两方法解算结果大部分精度相当，即振幅差异0.1mm以内，相位差异5°以内，但相比于传统方法解算结果，本发明所提方法解算结果精度普遍略有提升，特别在N2分潮U方向振幅精度优于0.28mm，Q1分潮W、S方向相位精度分别优于13.1°、7.4°。
步骤5：利用所求8个周日半周日主分潮的调和常数联合NAO.99b海潮模型估计三个长周期主分潮(MF、MM和SSA)的调和常数，对ZJWZ台站进行PPP动态解算时采用海潮负荷位移修正(其它解算策略与步骤1相同)。修正后，ZJWZ台站ENU方向RMS值由8.4、8.6、27.2mm变化为6.5、7.7、20.6mm，分别降低了22.6％、10.5％、24.3％，可见该站三维方向精度均不同程度提高。
以上流程利用潮族响应关系进行潮汐调和分析，不仅考虑了同潮族其余小分潮的影响，还避免了现有的交点改正方法所用平衡潮理论值与实际潮汐成分不一致而导致的误差，从而提高了GNSS海潮负荷位移反演精度。
具体实施时，本发明技术方案提出的方法可由本领域技术人员采用计算机软件技术实现自动运行流程，实现方法的系统装置例如存储本发明技术方案相应计算机程序的计算机可读存储介质以及包括运行相应计算机程序的计算机设备，也应当在本发明的保护范围内。
在一些可能的实施例中，提供一种顾及潮族响应关系的GNSS海潮负荷位移反演系统，包括以下模块，
第一模块，用于GNSS动态PPP数据采集与处理，包括采集GNSS连续观测站的长期观测数据，利用精密轨道和钟差产品进行动态PPP解算，解算时未采用海潮负荷位移改正，获取逐历元地心三维坐标时间序列，取平均得到时段解后转化为站心地平坐标时间序列，进行粗差探测与剔除，为潮汐调和分析的输入数据；
第二模块，用于明确周日半周日潮族响应关系，包括直接对周日半周日主分潮导纳在频率域进行线性插值而获取小分潮的导纳；
第三模块，用于构建顾及潮族响应关系的调和分析模型，所述顾及潮族响应关系的调和分析模型针对周日半周日主分潮进行调和分析时，在考虑交点改正的影响同时，考虑同潮族其余小分潮的影响；计算出小分潮对主分潮待估参数的累积影响系数，优化观测方程的系数矩阵；
第四模块，用于利用最小二乘方法求解主分潮的待估参数，进而获取振幅和格林尼治相位，作为调和常数；
第五模块，用于利用所求周日半周日主分潮的调和常数联合海潮模型估计长周期分潮的调和常数，在进行PPP动态解算时，进行海潮负荷位移修正，进一步提高GNSS定位精度。
在一些可能的实施例中，提供一种顾及潮族响应关系的GNSS海潮负荷位移反演系统，包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如上所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
在一些可能的实施例中，提供一种顾及潮族响应关系的GNSS海潮负荷位移反演系统，包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如上所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
以上所述的具体实施方式仅是对本发明所公开的方法作举例说明。本发明所属技术领域的技术人员可以对具体实施例进行合理修改、补充或替换，但并不会脱离本发明所附权利要求书所定义的范围。
Claims (7)
1.一种顾及潮族响应关系的GNSS海潮负荷位移反演方法，其特征在于，包括以下步骤：
步骤1，GNSS动态PPP数据采集与处理，包括采集GNSS连续观测站的长期观测数据，利用精密轨道和钟差产品进行动态PPP解算，解算时未采用海潮负荷位移改正，获取逐历元地心三维坐标时间序列，取平均得到时段解后转化为站心地平坐标时间序列，进行粗差探测与剔除后作为潮汐调和分析的输入数据；
步骤2，明确周日半周日潮族响应关系，包括直接对周日半周日主分潮导纳在频率域进行线性插值而获取小分潮的导纳；
步骤3，构建顾及潮族响应关系的调和分析模型，所述顾及潮族响应关系的调和分析模型针对周日半周日主分潮进行调和分析时，在考虑交点改正的影响同时，考虑同潮族其余小分潮的影响；实现方式为，建立求解主分潮调和常数的观测方程如下，
其中，Z0为常数项，a为线性趋势项，n表示主分潮个数，参数 为第i个分潮对应的Cartwright-Edden振幅/>的绝对值，Vi'(t)为第i个分潮t时刻的原始天文相角，分别对周日潮族各个分潮的原始天文相角作加90°处理，αj,i、βj,i为小分潮对第i个主分潮在j方向待估参数/>的累积影响系数，R(t)为残差；hj(t)表示t时刻GNSS在j方向的观测值；
参数uc、us分别为主要周日半周日分潮导纳u的实部和虚部；
计算出小分潮对主分潮待估参数的累积影响系数，优化观测方程的系数矩阵；
步骤4，利用最小二乘方法求解主分潮的待估参数，进而获取振幅和格林尼治相位，作为调和常数；
步骤5，利用所求周日半周日主分潮的调和常数联合海潮模型估计长周期分潮的调和常数，在进行PPP动态解算时，进行海潮负荷位移修正，进一步提高GNSS定位精度。
2.根据权利要求1所述顾及潮族响应关系的GNSS海潮负荷位移反演方法，其特征在于：步骤1中，将逐历元解取平均得到时段解，取平均得到时段解时，时段为5、10、15、20或30分钟。
3.根据权利要求1所述顾及潮族响应关系的GNSS海潮负荷位移反演方法，其特征在于：步骤2中，忽略插值方法异同对海潮负荷位移预测结果影响，将三次样条插值方法简化为以下线性关系，
其中，满足f2>fx>f1，表示小分潮频率fx位于同潮族相邻两个主分潮f1、f2之间，u1和u2为对应主分潮导纳，ux为小分潮导纳实部或虚部。
4.一种顾及潮族响应关系的GNSS海潮负荷位移反演系统，其特征在于：用于实现如权利要求1-3任一项所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
5.根据权利要求4所述顾及潮族响应关系的GNSS海潮负荷位移反演系统，其特征在于：包括以下模块，
第一模块，用于GNSS动态PPP数据采集与处理，包括采集GNSS连续观测站的长期观测数据，利用精密轨道和钟差产品进行动态PPP解算，解算时未采用海潮负荷位移改正，获取逐历元地心三维坐标时间序列，取平均得到时段解后转化为站心地平坐标时间序列，进行粗差探测与剔除，为潮汐调和分析的输入数据；
第二模块，用于明确周日半周日潮族响应关系，包括直接对周日半周日主分潮导纳在频率域进行线性插值而获取小分潮的导纳；
第三模块，用于构建顾及潮族响应关系的调和分析模型，所述顾及潮族响应关系的调和分析模型针对周日半周日主分潮进行调和分析时，在考虑交点改正的影响同时，考虑同潮族其余小分潮的影响；计算出小分潮对主分潮待估参数的累积影响系数，优化观测方程的系数矩阵；
第四模块，用于利用最小二乘方法求解主分潮的待估参数，进而获取振幅和格林尼治相位，作为调和常数；
第五模块，用于利用所求周日半周日主分潮的调和常数联合海潮模型估计长周期分潮的调和常数，在进行PPP动态解算时，进行海潮负荷位移修正，进一步提高GNSS定位精度。
6.根据权利要求4所述顾及潮族响应关系的GNSS海潮负荷位移反演系统，其特征在于：包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如权利要求1-3任一项所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
7.根据权利要求4所述顾及潮族响应关系的GNSS海潮负荷位移反演系统，其特征在于：包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如权利要求1-3任一项所述的一种顾及潮族响应关系的GNSS海潮负荷位移反演方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210926138.9A CN115453587B (zh) | 2022-08-03 | 2022-08-03 | 顾及潮族响应关系的gnss海潮负荷位移反演方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210926138.9A CN115453587B (zh) | 2022-08-03 | 2022-08-03 | 顾及潮族响应关系的gnss海潮负荷位移反演方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115453587A CN115453587A (zh) | 2022-12-09 |
| CN115453587B true CN115453587B (zh) | 2024-02-20 |
Family
ID=84296299
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210926138.9A Active CN115453587B (zh) | 2022-08-03 | 2022-08-03 | 顾及潮族响应关系的gnss海潮负荷位移反演方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115453587B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006330884A (ja) * | 2005-05-24 | 2006-12-07 | Chugoku Electric Power Co Inc:The | 潮流予測システム |
| JP2011145139A (ja) * | 2010-01-14 | 2011-07-28 | Hitachi Zosen Corp | 海面変位計測システム |
| CN106772498A (zh) * | 2016-11-21 | 2017-05-31 | 华东交通大学 | 一种gps位置时间序列噪声模型建立方法 |
| CN107544075A (zh) * | 2017-08-22 | 2018-01-05 | 中国科学院国家授时中心 | 基于精密单点定位与调和分析估计海潮负荷位移参数方法 |
| CN110082787A (zh) * | 2019-04-11 | 2019-08-02 | 华东师范大学 | 一种从gnss时间序列中提取周日半周日海潮信号的方法 |
| CN113761808A (zh) * | 2021-11-08 | 2021-12-07 | 长沙理工大学 | 一种基于gps和经验潮汐模型的地表潮汐位移获取方法及应用方法和系统 |
| CN114444330A (zh) * | 2022-04-02 | 2022-05-06 | 自然资源部第一海洋研究所 | 一种顾及平均效应的gnss海潮负荷位移的预测方法 |
| CN114595946A (zh) * | 2022-02-17 | 2022-06-07 | 广州市城市规划勘测设计研究院 | 一种海域理论最低潮面计算方法、系统、设备及介质 |
| CN114662348A (zh) * | 2022-05-25 | 2022-06-24 | 山东省科学院海洋仪器仪表研究所 | 一种基于gnss技术建立时变海潮负载位移模型的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030167124A1 (en) * | 2001-06-16 | 2003-09-04 | Bowden Dohn A. | Electronic tide tables |
-
2022
- 2022-08-03 CN CN202210926138.9A patent/CN115453587B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006330884A (ja) * | 2005-05-24 | 2006-12-07 | Chugoku Electric Power Co Inc:The | 潮流予測システム |
| JP2011145139A (ja) * | 2010-01-14 | 2011-07-28 | Hitachi Zosen Corp | 海面変位計測システム |
| CN106772498A (zh) * | 2016-11-21 | 2017-05-31 | 华东交通大学 | 一种gps位置时间序列噪声模型建立方法 |
| CN107544075A (zh) * | 2017-08-22 | 2018-01-05 | 中国科学院国家授时中心 | 基于精密单点定位与调和分析估计海潮负荷位移参数方法 |
| CN110082787A (zh) * | 2019-04-11 | 2019-08-02 | 华东师范大学 | 一种从gnss时间序列中提取周日半周日海潮信号的方法 |
| CN113761808A (zh) * | 2021-11-08 | 2021-12-07 | 长沙理工大学 | 一种基于gps和经验潮汐模型的地表潮汐位移获取方法及应用方法和系统 |
| CN114595946A (zh) * | 2022-02-17 | 2022-06-07 | 广州市城市规划勘测设计研究院 | 一种海域理论最低潮面计算方法、系统、设备及介质 |
| CN114444330A (zh) * | 2022-04-02 | 2022-05-06 | 自然资源部第一海洋研究所 | 一种顾及平均效应的gnss海潮负荷位移的预测方法 |
| CN114662348A (zh) * | 2022-05-25 | 2022-06-24 | 山东省科学院海洋仪器仪表研究所 | 一种基于gnss技术建立时变海潮负载位移模型的方法 |
Non-Patent Citations (10)
| Title |
|---|
| Improved Method for Estimating the Ocean Tide Loading Displacement Parameters by GNSS Precise Point Positioning and Harmonic Analysis;Rui Tu et al.;Journal of Surveying Engineering;全文 * |
| The ocean pole tide loading and its effect on GPS position time-series;Yujiao Niu et al.;Geophysical Journal International;全文 * |
| 一种建立南海浅海海域高精度潮汐模型方法的研究;雷宁;付延光;杨龙;周兴华;阳凡林;王朝阳;孙维康;;海洋科学进展(第03期);全文 * |
| 中国沿岸验潮站潮汐调和常数的精度评估;暴景阳;许军;;海洋测绘;20130125(第01期);全文 * |
| 利用GPS动态PPP技术求解海潮负荷位移;赵红;涂锐;刘智;蒋光伟;;测绘学报(第08期);全文 * |
| 利用短期静态PPP结果反演海潮负荷位移;魏国光等;大地测量与地球动力学;第39卷(第6期);全文 * |
| 南海TOPEX/POSEIDON高度计资料的正交响应法潮汐分析;李培良, 左军成, 李磊, 黄娟;海洋与湖沼(第03期);全文 * |
| 短期潮汐潮流数据的正交潮响应分析研究;唐岩 等;武汉大学学报信息科学版;第35卷(第10期);全文 * |
| 面向城市复杂环境的3种多频多系统GNSS单点高精度定位方法及性能分析;耿江辉等;测绘学报;第49卷(第1期);全文 * |
| 魏国光 ; 王琪洁 ; 彭葳 ; 郭良琦 ; .利用短期静态PPP结果反演海潮负荷位移.大地测量与地球动力学.2019,(第06期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115453587A (zh) | 2022-12-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111596321B (zh) | 利用非差改正数的多gnss多路径误差恒星日滤波方法及系统 | |
| CN108387912B (zh) | 一种Multi-GNSS精密单点定位的解算方法 | |
| CN111125609A (zh) | 一种基于双指数驱动的电离层三维电子密度重构方法 | |
| CN107544075A (zh) | 基于精密单点定位与调和分析估计海潮负荷位移参数方法 | |
| CN116540303B (zh) | 一种历元间差分观测方程建立方法和地震位移计算方法 | |
| CN103869344A (zh) | 一种抗差估计方法 | |
| CN110146904B (zh) | 一种适用于区域电离层tec的精确建模方法 | |
| CN115453587B (zh) | 顾及潮族响应关系的gnss海潮负荷位移反演方法及系统 | |
| CN112230254B (zh) | 一种gps载波相位多径误差的校正方法及装置 | |
| Choy et al. | An evaluation of various ionospheric error mitigation methods used in single frequency PPP | |
| CN114779301B (zh) | 一种基于广播星历的卫星导航实时精密单点定位方法 | |
| CN115343734B (zh) | 一种基于双线性内插半球模型的gnss变形监测方法 | |
| CN115980317B (zh) | 基于修正相位的地基gnss-r数据土壤水分估算方法 | |
| CN115308781B (zh) | 基于bdgim辅助的相位平滑伪距高精度时间传递方法 | |
| CN115421172B (zh) | 一种基于实时与准实时结合的北斗变形监测方法 | |
| CN114488227B (zh) | 一种基于空间相关性的多路径误差改正方法 | |
| CN116147622A (zh) | 一种基于图优化的组合导航系统融合定位方法 | |
| CN111766614B (zh) | 一种守时实验室联合组网的亚纳秒实时时间比对方法 | |
| CN114048585A (zh) | 一种电离层模型事后分析方法及装置 | |
| Zhou et al. | Ocean Tide Loading Displacement Parameters Estimated From GNSS-Derived Coordinate Time Series Considering the Effect of Mass Loading in Hong Kong | |
| CN110764122A (zh) | 一种单频gps接收机的精密单点定位方法 | |
| Kass et al. | Global GPS data analysis at the National Geodetic Survey | |
| Kuang et al. | Assessment of L5 frequency on real-time GPS orbit and clock estimation | |
| CN113899957B (zh) | 基于双差相位观测值的室外天线绝对标定方法及系统 | |
| CN115712130B (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 |