基于负荷形变的区域地面重力变化连续监测方法及装置与流程

1.本发明涉及地质测绘技术领域，尤其涉及一种基于负荷形变的区域地面重力变化连续监测方法及装置。


背景技术：

2.区域地面重力变化连续监测方法有地面重力台站观测、流动重力复测以及卫星重力观测等。
3.地面重力台站连续观测，是利用重力台站上的连续重力观测仪进行观测，获得准确、可靠、连续的重力随时间变化的观测资料及辅助观测资料，并利用平均滤波器将重力仪观测数据采样率降低，抽取整点值，开展间断、掉格、突跳等粗差探测和预处理，并完成零漂改正，形成重力变化时间序列。流动重力复测，是利用相对重力仪对重力测点进行定期复测，以获取重力测点之间的重力段差及复测间隔时间内的重力段差变化信息。重力变化由多期复测的差值计算。卫星重力观测，是以卫星为载体，利用卫星本身为重力传感器或卫星所携带的重力传感器，观测由地球重力场引起的卫星轨道摄动，以这些数据资料来反演和恢复地球重力场的方法和技术。
4.由于重力台站建站花费较大，且重力台站数目较少，虽然流动重力测量点相对重力台站要丰富，但是不能连续测量，一般每年复测1至2次。重力台站测量和流动重力测量主要是以测点形式开展，在空间分辨率上也存在不足。卫星重力测量能够提供全球重力场变化，并提供面上的重力变化，但是分辨率较低，一般空间分辨率300千米以上，时间分辨率1个月左右。
5.如何方便地对地面重力变化进行准确监测是一个亟需解决的问题。


技术实现要素：

6.本发明提供一种基于负荷形变的区域地面重力变化连续监测方法及装置，用以解决现有技术中监测成本高且监测分辨率低的缺陷，实现利用分布广泛的cors站的观测数据得到高分辨率的重力变化监测数据。
7.本发明提供一种基于负荷形变的区域地面重力变化连续监测方法，包括：确定目标区域内每个连续运行基准站cors的非线性大地高变化时间序列；基于每个所述cors的非线性大地高变化时间序列和地表负荷对应的大地高变化值，得到每个所述cors的残差大地高时间序列；基于所有所述cors的残差大地高时间序列，确定所述目标区域的残差等效水高变化数据；基于所述目标区域的残差等效水高变化数据，确定所述目标区域的残差重力变化数据；基于所述目标区域的残差重力变化数据和所述地表负荷对应的重力变化数据，确定所述目标区域的重力变化数据。
8.根据本发明提供的一种基于负荷形变的区域地面重力变化连续监测方法，所述基于所有所述cors的残差大地高时间序列，确定所述目标区域的残差等效水高变化数据，包括：将所有所述cors的残差大地高时间序列代入地面等效水高变化引起的大地高变化数据观测方程；以残差等效水高变化数据为未知数，对所述大地高变化数据观测方程进行求解；所述大地高变化数据观测方程表示为：其中，为地面计算点的地心纬度和经度，为地面积分流动点的地心纬度和经度，l为地面积分流动点到地面计算点的空间距离，为残差等效水高变化数据，为径向格林函数，为球面角距离，为构造垂直形变、地下水均衡垂直形变以及噪声影响因子，为水的密度，g为万有引力常数。
9.根据本发明提供的一种基于负荷形变的区域地面重力变化连续监测方法，所述对所述大地高变化数据观测方程进行求解，包括：以所述cors为中心按照高斯函数配置平滑系数；基于地面及外部重力位形变与地面其他形变量之间的动力学约束，地球内部密度调整与地面形变量之间的地球物理关系以及体元的瞬时应力与瞬时形变位移关系，设置所述大地高变化数据观测方程求解的边界条件；基于所述平滑系数、所述边界条件以及标准差最小原则约束，通过基于岭估计的正则化解算所述大地高变化数据观测方程。
10.根据本发明提供的一种基于负荷形变的区域地面重力变化连续监测方法，所述地表负荷对应的重力变化数据通过以下方式确定：获取所述目标区域的地表负荷变化格网数据；对所述地表负荷变化格网数据进行规格化负荷球谐化分析，确定每个所述cors的所述地表负荷对应的等效水高变化数据；基于每个所述cors的所述地表负荷对应的等效水高变化数据，确定所述地表负荷对应的重力变化数据。
11.根据本发明提供的一种基于负荷形变的区域地面重力变化连续监测方法，所述地表负荷对应的大地高变化值通过以下方式确定：获取所述目标区域的地表负荷变化格网数据；对所述地表负荷变化格网数据进行规格化负荷球谐化分析，确定每个所述cors的所述地表负荷对应的大地高变化值。
12.根据本发明提供的一种基于负荷形变的区域地面重力变化连续监测方法，所述确定目标区域内每个卫星导航定位基准站cors的非线性大地高变化时间序列，包括：获取每个所述cors在目标时间段内的观测数据；
解算每个所述cors的观测数据，得到每个所述cors的站点坐标变化时间序列；基于每个所述cors的站点坐标变化时间序列，确定每个所述cors的非线性大地高变化时序。
13.根据本发明提供的一种基于负荷形变的区域地面重力变化连续监测方法，所述地表负荷包括大气负荷、地表水负荷和非潮汐海潮负荷中的至少一种。
14.本发明还提供一种地面重力变化监测装置，包括：第一处理模块，用于确定目标区域内每个连续运行基准站cors的非线性大地高变化时间序列；第二处理模块，用于基于每个所述cors的非线性大地高变化时间序列和地表负荷对应的大地高变化值，得到每个所述cors的残差大地高时间序列；第三处理模块，用于基于所有所述cors的残差大地高时间序列，确定所述目标区域的残差等效水高变化数据；第四处理模块，用于基于所述目标区域的残差等效水高变化数据，确定所述目标区域的残差重力变化数据；第五处理模块，用于基于所述目标区域的残差重力变化数据和所述地表负荷对应的重力变化数据，确定所述目标区域的重力变化数据。
15.本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述基于负荷形变的区域地面重力变化连续监测方法。
16.本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述基于负荷形变的区域地面重力变化连续监测方法。
17.本发明还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如上述任一种所述基于负荷形变的区域地面重力变化连续监测方法。
18.本发明提供的基于负荷形变的区域地面重力变化连续监测方法及装置，通过获取离散分布的连续运行基准站cors的观测数据，将地表负荷对应的大地高变化值进行移去，可以得到残差大地高时间序列，进而反演得到准确的残差等效水高变化格网数据，降低已知的地表负荷对反演求解的干扰，进而得到更加准确的残差重力变化数据，最终和精确的地表负荷对应的重力变化数据叠加得到目标区域准确的重力变化数据，既能保证数据在时间和空间上的分辨率，又保证了数据的准确性，还有效降低了成本。
附图说明
19.为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
20.图1是本发明提供的基于负荷形变的区域地面重力变化连续监测方法的流程示意图之一；图2是本发明提供的基于负荷形变的区域地面重力变化连续监测方法的流程示意
图之二；图3是本发明提供的a地区的部分cors大地高变化时间序列示意图；图4是本发明提供的a地区的地面重力变化格网时序示意图；图5是本发明提供的地面重力变化监测装置的结构示意图；图6是本发明提供的电子设备的结构示意图。
具体实施方式
21.为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
22.下面结合图1-图6描述本发明的基于负荷形变的区域地面重力变化连续监测方法及装置。
23.本发明实施例的基于负荷形变的区域地面重力变化连续监测方法的执行主体可以是控制器，或者，还可以是服务器等，此处对执行主体的具体类型不作限制。下面以执行主体为地面重力变化监测平台的服务器为例来对本发明实施例的基于负荷形变的区域地面重力变化连续监测方法进行说明。
24.针对现有技术中地面重力变化监测中存在的不足，例如流动重力复测时间长、卫星重力观测空间分辨率低等，本发明实施例的地面重力监测方法提出基于负荷形变理论来建立地面形变和地面重力变化的关系模型，利用cors（continuously operating reference station，连续运行基准站）组网和gnss（global navigation satellite system，全球卫星导航定位系统）连续观测的技术优势，解算地面形变信息，基于形变动力学及移去恢复技术，同化水文、气象等观测资料，建立区域地面重力变化模型。本发明可以有效利用现有的gnss连续观测站点，即连续运行基准站cors，又称卫星导航定位基准站，来有效提升地面重力变化监测水平，降低监测成本。
25.如图1所示，本发明实施例的基于负荷形变的区域地面重力变化连续监测方法主要包括步骤110、步骤120、步骤130、步骤140和步骤150。
26.步骤110，确定目标区域内每个连续运行基准站cors的非线性大地高变化时间序列。
27.需要说明的是，目标区域为需要进行地面重力变化监测的区域。在本实施方式中，可以通过利用目标区域内cors网络中各个站点的观测数据来重构大地高时间序列。
28.在一些实施例中，确定目标区域内每个卫星导航定位基准站cors的非线性大地高变化时间序列包括以下步骤。
29.可以先获取每个cors在目标时间段内的观测数据。目标时间段为需要对地面重力变化进行监测的时间段。例如，目标时间段可以是3年-5年，目标时间段的长度可以根据实际需要进行确定。
30.在获取到观测数据后，需要对观测数据进行解算，得到各个cors的大地高时间序列。
31.cors观测数据的解算一般采用gamit/globk软件进行。gamit/globk软件对每站每
日的gnss观测数据进行解算，观测数据包括观测站位置、接收机钟差与卫星钟差等在内的多种参数。
32.在本实施方式中，可以获得测站和卫星轨道的单日区域松弛解，单日解给出了区域测站、极移和卫星等参数的松弛解和方差-协方差矩阵，再利用globk软件进行整体平差，得到国际框架下的站点坐标变化时间序列。
33.在此基础上，再对站点坐标变化时间序列进行粗差探测、线性项移去、非线性项多周期重构。对cors站gnss解大地高和连续重力时间序列进行预处理，再进行数据统计，统计最大最小值、标准差等信息，以3倍标准差（3σ）准则进行初步的粗差探测以及跳变修复。
34.再基于正交余弦、切比谢夫等不同基函数的低频参数重构算法，构筑时序线性项。通过分离时序线性项，对非线性部分再次利用3σ准则进行信噪分离，来实现二次粗差探测。
35.通过傅里叶分析方法，利用大地高变化时序的频谱特征，进行周期提取，并利用选定的周期进行信号重构，基于前后残差变化小于观测时序标准差的2%，获得非线性大地高变化时序。
36.换言之，可以解算每个cors的观测数据，得到每个cors的站点坐标变化时间序列，再对每个cors的站点坐标变化时间序列进一步解算，确定出每个cors的非线性大地高变化时序。
37.在本实施方式中，通过对cors的观测数据进行解算，可以得到目标区域内每个cors站点的非线性大地高变化时序，而大地高与地面重力均与地面等效水高之间存在关联关系，进而方便通过大地高数据来反演地面等效水高数据，再通过地面等效水高数据来确定地面重力变化数据。
38.步骤120，基于每个cors的非线性大地高变化时间序列和地表负荷对应的大地高变化值，得到每个cors的残差大地高时间序列。
39.受地表负荷变化影响，cors的大地高随时间变化。地表负荷变化引起的cors的大地高变化满足负荷格林函数积分式。
40.地表负荷包括大气负荷、地表水负荷和非潮汐海潮负荷中的至少一种。
41.大气中含有较多的大气水。地表水包括土壤水、江河湖库水以及冰川冰盖雪山水等。非潮汐海潮包括除了海洋潮汐以外的海平面变化等现象。
42.在地表负荷发生较大变化的情况下，当前区域容易遭遇极端天气、地震、潮汐等动力环境异常的情况，并引发地质灾害。地质灾害的发生与地面的垂直形变以及重力变化有较大的关联性，通过离散分布的cors能够实现对地面垂直形变的连续监测。
43.在一些实施例中，地表负荷对应的大地高变化值通过以下方式确定。
44.可以先获取目标区域的地表负荷变化格网数据，再对地表负荷变化格网数据进行规格化负荷球谐化分析，并确定每个cors的地表负荷对应的大地高变化值。
45.大气负荷、地表水负荷和非潮汐海潮负荷，可用地面等效水高表示。地面等效水高变化的规格化负荷球谐展开式为：根据地球负荷形变动力学理论，地面等效水高变化引起的大地高变
化为：其中，为地面计算点的地心纬度和经度，为n阶m次规格化负荷球谐系数，为n阶m次规格化缔合勒让德函数，为n阶位负荷（勒夫）数，为n阶径向负荷数，为水的密度，为固体地球平均密度，g为万有引力常数，m为地球总质量，r为地球平均半径，为地面平均重力。
46.在本实施方式中，通过对地表负荷变化格网数据的规格化负荷球谐化分析，能够确定出各个cors因地表负荷所产生的大地高变化，进而方便进一步确定出每个cors的残差大地高时间序列。
47.在一些实施例中，可以先后对大气负荷、地表水负荷和非潮汐海潮负荷所引起的大地高变化从非线性大地高变化时间序列中移去，最后得到每个cors的残差大地高时间序列。
48.对于大气负荷，可以将全球大气压月平均格网数据转换为等效水高格网，再将转换后的全球大气等效水高格网数据进行规格化球谐展开，获得全球大气等效水高球谐系数，并保存球谐系数到360阶。
49.在一些实施例中，可以由球谐系数计算地面大地高变化。从区域大气观测中移去全球大气等效水高值，计算区域残差大气等效水高。由残差大气等效水高，按格林函数法计算残差地面大地高变化的影响。将模型大气负荷影响与残差大气负荷影响相加，得到目标区域大气负荷影响，其包括对大地高变化的影响。
50.在此种情况下，再从非线性大地高变化时序中，移去大气负荷对大地高的影响。
51.对于地表水负荷，地表水负荷可利用全球陆面数据同化系统模型进行计算。地表水负荷可以根据地表径流情况来单独计算江河湖库水负荷影响。
52.在一些实施例中，主要考虑到土壤水负荷影响，可以考虑只对土壤水负荷进行确定。在此种情况下，可以将全球土壤水等效水高格网进行规格化球谐展开，得到360阶等效水高球谐系数，并由土壤水等效水高球谐系数计算对大地高的负荷影响。
53.在此种情况下，再从非线性大地高变化时序中，移去地表水负荷对大地高的影响。
54.对于非潮汐海洋负荷，主要考虑到对沿海的目标区域的影响。在一些实施例中，可以根据目标区域距离海洋的远近来确定是否需要考虑非潮汐海洋负荷的影响。
55.在本实施方式中，可以利用全球海平面变化数据，并基于负荷形变理论，采用球谐分析方法，计算非潮汐海洋负荷对大地高变化的影响。
56.在此种情况下，再从非线性大地高变化时序中，移去非潮汐海洋负荷对大地高的影响。
57.步骤130，基于所有cors的残差大地高时间序列，确定目标区域的残差等效水高变化数据。
58.可以理解的是，根据地球负荷形变理论，大地高与等效水高之间满足负荷格林函
数。
59.在本实施方式中，可以根据大地高与等效水高之间的关系式建立残差大地高与对应的残差等效水高之间的关系。
60.在一些实施例中，可以将所有cors的残差大地高时间序列代入地面等效水高变化引起的大地高变化数据观测方程，并以残差等效水高变化数据为未知数，对大地高变化数据观测方程进行求解。
61.大地高变化数据观测方程表示为：其中，为地面计算点的地心纬度和经度，为地面积分流动点的地心纬度和经度，l为地面积分流动点到地面计算点的空间距离，为残差等效水高变化数据，为径向格林函数，为球面角距离，为构造垂直形变、地下水均衡垂直形变以及噪声影响因子，为水的密度，g为万有引力常数。
62.在本实施方式中，可以根据以获取的cors的残差大地高时间序列反演推算出其所对应的残差等效水高变化数据。
63.在一些实施例中，对大地高变化数据观测方程进行求解，包括以下过程。
64.可以按照cors为中心按照高斯函数配置平滑系数。通过平滑系数可以控制解算收敛的速度，并根据实际收敛情况调整比例系数。以cors为中心能最大限度抑制格网化代表性误差。
65.在一些实施例中，平滑系数可以配置为：其中，为地面计算点的地心纬度和经度，r为地球平均半径，σ为各个地面计算点的离散程度。
66.在一些实施例中，可以根据弹性地球动力学方程组，设置边界条件。弹性地球动力学方程组，包括角动量守恒方程、连续性方程和本构方程。
67.角动量守恒方程是自转地球的线性弹性运动方程。可以通过利用角动量守恒方程，建立地面及外部重力位形变与地面其他形变量之间的动力学约束。
68.连续性方程中形变前后质量守恒。利用连续性方程可以建立地球内部密度调整与地面形变量之间的地球物理关系，位形变与密度形变满足泊松方程。
69.本构方程是体元的瞬时应力与瞬时形变位移关系。在静力平衡初始状态，假设各向同性，对自转地球微小弹性形变的完整动力学方程建立以下边界条件：解在地球质心处正则；地球表面上的应力为零；地球内部不连续面处位移形变与应力连续；地球内外位形变在地球表面及内部不连续面上相等。形变位移u各分量和位形变有相同调和级数形式。
70.在此种情况下，可以基于地面及外部重力位形变与地面其他形变量之间的动力学约束，地球内部密度调整与地面形变量之间的地球物理关系以及体元的瞬时应力与瞬时形变位移关系，设置大地高变化数据观测方程求解的边界条件。再根据平滑系数、边界条件以及标准差最小原则约束，通过基于岭估计的正则化解算大地高变化数据观测方程。
71.在解算过程中，可以根据积分半径，设置计算点远区影响为0，进行反演收敛约束。可以以标准差最小原则约束，调整迭代次数，并综合采用基于岭估计的正则化方法解算法方程。
72.岭估计是在解算中增加一个正常矩阵，使得解的奇异性可能性降低，提高解算稳定性。
73.在本实施方式中，以标准差最小和动力学规律等设置约束条件，整个解算过程一般需要进行多次迭代计算，直到满足限定条件，实现了目标区域高分辨负荷等效水高建模反演，进而可以获取更加准确的数据。
74.步骤140，基于目标区域的残差等效水高变化数据，确定目标区域的残差重力变化数据。
75.在此种情况下，可以根据等效水高数据与地面重力之间满足的格林函数积分式来计算目标区域的残差等效水高变化数据所对应的残差重力变化数据。
76.步骤150，基于目标区域的残差重力变化数据和地表负荷对应的重力变化数据，确定目标区域的重力变化数据。
77.在确定出残差重力变化数据后，再将地表负荷对应的重力变化数据进行恢复。具体地，可以将地表负荷对应的重力变化数据与残差重力变化数据进行叠加，最终得到目标区域的重力变化数据。
78.在一些实施例中，地表负荷对应的重力变化数据通过以下方式确定。
79.可以先获取目标区域的地表负荷变化格网数据，再对地表负荷变化格网数据进行规格化负荷球谐化分析，确定每个cors的地表负荷对应的等效水高变化数据。
80.在此基础上，再根据每个cors的地表负荷对应的等效水高变化数据，确定地表负荷对应的重力变化数据。
81.在一些实施例中，可以利用反演获得的残差等效水高变化数据，可以计算对地面重力变化，具体可以采用下面的公式计算：与上述实施例类似的，大气负荷、地表水负荷和非潮汐海潮负荷，可用地面等效水高表示。地面等效水高变化的规格化负荷球谐展开式为：地表负荷对应的重力变化数据可以采用地表负荷对应的地面等效水高变化引起的地面重力变化公式计算：
；其中，为重力格林函数，为地面计算点的地心纬度和经度，为n阶m次规格化负荷球谐系数，为n阶m次规格化缔合勒让德函数，为n阶位负荷（勒夫）数，为n阶径向负荷数，为水的密度，为固体地球平均密度，g为万有引力常数，m为地球总质量，r为地球平均半径。
82.在本实施方式中，地表负荷数据为已知的影响地面重力负荷的因素，根据地表负荷变化数据所求得的地表负荷对应的重力变化数据是地面重力变化中可以准确得到的一部分变化量。
83.在反演推算残差等效水高变化数据时，通过将非线性大地高变化时间序列中地表负荷对应的大地高变化值进行移去，可以将已知的地表负荷对应的等效水高变化数据进行去除，降低已知的地表负荷数据对未知数据在反演解算过程中的干扰，进而提高解算精度，从而提高数据的准确性。
84.需要说明的是，现有技术中，重力台站和流动重力复测是以点的形式对地面重力进行观测。重力台站由于站点稀少，空间分辨率低，而流动重力复测的时间分辨率较低。卫星重力的空间分辨率约在300km，时间分辨率1月左右，分辨率也偏低。本发明利用cors站网，综合水文气象数据，施加负荷动力学约束和移去恢复技术，构建了残差负荷等效水高模型，实现了地面重力变化的高时空分辨率监测，时间分辨率可达1周，空间分辨率5km左右。
85.传统区域地面重力变化主要利用流动重力复测手段，在固定点上开展多期观测，需要花费较大的人力物力。本发明实施例的地面重力变化检测方法可以利用现有的cors进行观测，通过解算三维坐标变化，进行数据综合反演，实现对地面重力变化监测，可以有效节约观测成本，提高监测效率。
86.根据本发明实施例提供的基于负荷形变的区域地面重力变化连续监测方法，通过获取离散分布的连续运行基准站cors的观测数据，将地表负荷对应的大地高变化值进行移去，可以得到残差大地高时间序列，进而反演得到准确的残差等效水高变化格网数据，降低已知的地表负荷对反演求解的干扰，进而得到更加准确的残差重力变化数据，最终和精确的地表负荷对应的重力变化数据叠加得到目标区域准确的重力变化数据，既能保证数据在时间和空间上的分辨率，又保证了数据的准确性，还有效降低了成本。
87.参照图2，下面以目标区域为内陆区域a地区为例来对本发明实施例的重力变化监测方法进行详细说明。
88.在本实施方式中，可以利用a地区的cors网进行观测并获取观测数据，再对cors网观测数据进行解算与数据预处理。可以利用gamit解算获得测站和卫星轨道的单日区域松弛解，这个单日解给出了区域测站、极移和卫星等参数的松弛解和方差-协方差矩阵，利用globk软件进行整体平差，得到国际框架下的站点坐标变化时间序列。
89.cors站天解大地高时序分析、线性项分离与非线性变化提取流程为：针对每座cors站大地高天解时间序列，以切比雪夫函数为基函数，对时序进行低频参数估计，分离线性项；按低频参数重构大地高非线性变化时序，作为粗差探测的基准，按3倍残差标准差，探
测并剔除粗差，生成大地高非线性变化时间序列。对时序散点进行周期估计，以60个周期，估计周期项参数，然后重构大地高非线性变化时序。
90.如图3所示，图3展示了部分cors的大地高变化时间序列，图3中，横坐标表示时间，纵坐标表示大地高。
91.在本实施方式中，由于a地区距离海洋较远，非潮汐海洋负荷影响很小，可忽略不计。但是在沿海区域，需要考虑非潮汐海洋负荷影响。因此，在a地区只需要移去大气负荷和地表水负荷对大地高和地面重力的影响，进而得到残差观测量，并通过残差观测方程反演解算得到残差等效水高，进而再求出残差地面重力变化。
92.在一些实施例中，可以利用全球大气模型和区域数字地形模型，采用球谐分析方法转换为等效水高球谐系数，通过负荷形变公式计算大气负荷对大地高变化和地面重力变化的影响。
93.在一些实施例中，可以利用地表水模型和数字地形模型数据，计算地表水负荷对大地高变化影响。全球地表水等效水高月变化采用全球大陆同化系统数据。将地表水等效水高采用球谐分析方法转换为球谐系数，分别计算负荷对大地高变化和地面重力变化的影响。
94.根据负荷动力学方程，将残差大地高变化周平均，利用每周点数据建立观测方程和法方程。以cors站为中心按高斯函数配置平滑系数，平滑半径选择1度。
95.在此基础上，可以利用反演获得的残差等效水高计算对地面重力变化影响，即得到残差地面重力变化，并恢复已知的大气负荷和地表水负荷对地面重力变化的影响，获得区域地面重力变化格网数值模型，实现区域地面重力变化连续监测，如图4所示，图中示出了a地区的地面重力变化格网时序示意图，其中横坐标表示时间，纵坐标表示重力，重力通过重力加速度的单位ugal来将进行表征，重力的最大值、平均值和最小值分别采用不同的曲线进行展示。
96.下面对本发明提供的地面重力变化监测装置进行描述，下文描述的地面重力变化监测装置与上文描述的基于负荷形变的区域地面重力变化连续监测方法可相互对应参照。
97.如图5所示，本发明实施例的地面重力变化监测装置包括第一处理模块510、第二处理模块520、第三处理模块530、第四处理模块540以及第五处理模块550。
98.第一处理模块510用于确定目标区域内每个连续运行基准站cors的非线性大地高变化时间序列；第二处理模块520用于基于每个cors的非线性大地高变化时间序列和地表负荷对应的大地高变化值，得到每个cors的残差大地高时间序列；第三处理模块530用于基于所有cors的残差大地高时间序列，确定目标区域的残差等效水高变化数据；第四处理模块540用于基于目标区域的残差等效水高变化数据，确定目标区域的残差重力变化数据；第五处理模块550用于基于目标区域的残差重力变化数据和地表负荷对应的重力变化数据，确定目标区域的重力变化数据。
99.根据本发明实施例提供的重力变化监测装置，通过获取离散分布的连续运行基准站cors的观测数据，将地表负荷对应的大地高变化值进行移去，可以得到残差大地高时间
序列，进而反演得到准确的残差等效水高变化格网数据，降低已知的地表负荷对反演求解的干扰，进而得到更加准确的残差重力变化数据，最终和精确的地表负荷对应的重力变化数据叠加得到目标区域准确的重力变化数据，既能保证数据在时间和空间上的分辨率，又保证了数据的准确性，还有效降低了成本。
100.在一些实施例中，第三处理模块530还用于将所有cors的残差大地高时间序列代入地面等效水高变化引起的大地高变化数据观测方程；以残差等效水高变化数据为未知数，对大地高变化数据观测方程进行求解；大地高变化数据观测方程表示为：其中，为地面计算点的地心纬度和经度，为地面积分流动点的地心纬度和经度，l为地面积分流动点到地面计算点的空间距离，为残差等效水高变化数据，为径向格林函数，为球面角距离，为构造垂直形变、地下水均衡垂直形变以及噪声影响因子，为水的密度，g为万有引力常数。
101.在一些实施例中，第三处理模块530还用于以cors为中心按照高斯函数配置平滑系数；基于地面及外部重力位形变与地面其他形变量之间的动力学约束，地球内部密度调整与地面形变量之间的地球物理关系以及体元的瞬时应力与瞬时形变位移关系，设置大地高变化数据观测方程求解的边界条件；基于平滑系数、边界条件以及标准差最小原则约束，通过基于岭估计的正则化解算大地高变化数据观测方程。
102.在一些实施例中，第五处理模块550还用于获取目标区域的地表负荷变化格网数据；对地表负荷变化格网数据进行规格化负荷球谐化分析，确定每个cors的地表负荷对应的等效水高变化数据；基于每个cors的地表负荷对应的等效水高变化数据，确定地表负荷对应的重力变化数据。
103.在一些实施例中，第二处理模块520还用于获取目标区域的地表负荷变化格网数据；对地表负荷变化格网数据进行规格化负荷球谐化分析，确定每个cors的地表负荷对应的大地高变化值。
104.在一些实施例中，第一处理模块510还用于获取每个cors在目标时间段内的观测数据；解算每个cors的观测数据，得到每个cors的站点坐标变化时间序列；基于每个cors的站点坐标变化时间序列，确定每个cors的非线性大地高变化时序。
105.在一些实施例中，地表负荷包括大气负荷、地表水负荷和非潮汐海潮负荷中的至少一种。
106.图6示例了一种电子设备的实体结构示意图，如图6所示，该电子设备可以包括：处理器（processor）610、通信接口（communications interface）620、存储器（memory）630和通信总线640，其中，处理器610，通信接口620，存储器630通过通信总线640完成相互间的通
信。处理器610可以调用存储器630中的逻辑指令，以执行基于负荷形变的区域地面重力变化连续监测方法，该方法包括：确定目标区域内每个连续运行基准站cors的非线性大地高变化时间序列；基于每个cors的非线性大地高变化时间序列和地表负荷对应的大地高变化值，得到每个cors的残差大地高时间序列；基于所有cors的残差大地高时间序列，确定目标区域的残差等效水高变化数据；基于目标区域的残差等效水高变化数据，确定目标区域的残差重力变化数据；基于目标区域的残差重力变化数据和地表负荷对应的重力变化数据，确定目标区域的重力变化数据。
107.此外，上述的存储器630中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器（rom，read-only memory）、随机存取存储器（ram，random access memory）、磁碟或者光盘等各种可以存储程序代码的介质。
108.另一方面，本发明还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述各方法所提供的基于负荷形变的区域地面重力变化连续监测方法，该方法包括：确定目标区域内每个连续运行基准站cors的非线性大地高变化时间序列；基于每个cors的非线性大地高变化时间序列和地表负荷对应的大地高变化值，得到每个cors的残差大地高时间序列；基于所有cors的残差大地高时间序列，确定目标区域的残差等效水高变化数据；基于目标区域的残差等效水高变化数据，确定目标区域的残差重力变化数据；基于目标区域的残差重力变化数据和地表负荷对应的重力变化数据，确定目标区域的重力变化数据。
109.又一方面，本发明还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各方法提供的基于负荷形变的区域地面重力变化连续监测方法，该方法包括：确定目标区域内每个连续运行基准站cors的非线性大地高变化时间序列；基于每个cors的非线性大地高变化时间序列和地表负荷对应的大地高变化值，得到每个cors的残差大地高时间序列；基于所有cors的残差大地高时间序列，确定目标区域的残差等效水高变化数据；基于目标区域的残差等效水高变化数据，确定目标区域的残差重力变化数据；基于目标区域的残差重力变化数据和地表负荷对应的重力变化数据，确定目标区域的重力变化数据。
110.以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
111.通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上
述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如rom/ram、磁碟、光盘等，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行各个实施例或者实施例的某些部分所述的方法。
112.最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。