CN111833208A - A groundwater storage monitoring method and system based on vertical line deviation disturbance - Google Patents

Abstract

The invention discloses a method and a system for monitoring underground water reserves based on vertical deviation disturbance. The underground water reserve monitoring method based on vertical deviation disturbance comprises the following steps: acquiring the running speed of a measuring carrier; constructing a vertical deviation disturbance two-order random process optimization model based on the driving speed; the vertical deviation disturbance two-order random process optimization model comprises a two-order random process model and a differentiator; calculating a vertical deviation disturbance component based on a vertical deviation disturbance two-step random process optimization model; the vertical deviation disturbance component comprises a vertical deviation disturbance north-south component and a vertical deviation disturbance east-west component; and obtaining the underground water reserve variation of the area where the measurement carrier is located by inversion of the vertical deviation disturbance component. The invention solves the problem that large-scale and large-scale underground water data cannot be acquired, and improves the monitoring efficiency of underground water reserves.

Description

A groundwater storage monitoring method and system based on vertical line deviation disturbance

technical field

The invention relates to the technical field of gravity measurement, in particular to a groundwater storage monitoring method and system based on vertical line deviation disturbance.

Background technique

With the development of social economy and population growth, industrial and agricultural water consumption is increasing day by day, and there is a contradiction between the supply and demand of water resources used by human beings. Groundwater has become the main water supply source for large and medium-sized cities in my country. Groundwater storage has a great impact on natural ecological balance and sustainable development, and it is of great significance to carry out monitoring of groundwater.

Due to the irregular shape of the Earth's surface and the uneven density of its interior, there is a difference between true and normal gravity. The magnitude of the difference appears as gravity anomaly, and the direction appears as deflection of the vertical (DOV). The vertical deviation signal can be decomposed into two parts, high frequency and medium and low frequency. The vertical line deviation is the basic observation of the earth's gravitational field, and it contains rich high-frequency information of the gravitational field. The vertical line deviation can better reflect the real information of the gravity field, and has very important applications in the fields of resource exploration, geophysical inversion, satellite precise orbit, volcano observation, earthquake and auxiliary navigation. High-frequency gravity components are important information urgently needed in many research fields. Therefore, it is of great significance to obtain high-resolution and high-precision high-frequency vertical line deviation information, and to study its measurement and approximation methods.

Most of the traditional groundwater storage monitoring methods are realized by obtaining data from fixed monitoring network sites of water resources on the earth's surface. This method is affected by the number and distribution of monitoring sites, and cannot obtain large-scale and large-scale groundwater data, resulting in low monitoring efficiency.

SUMMARY OF THE INVENTION

Based on this, it is necessary to provide a groundwater storage monitoring method and system based on vertical deviation disturbance, which solves the problem of inability to obtain large-scale and large-scale groundwater data and improves the monitoring efficiency of groundwater storage.

For achieving the above object, the present invention provides the following scheme:

A groundwater storage monitoring method based on vertical line deviation disturbance, comprising:

Obtain the traveling speed of the measurement carrier;

A vertical line deviation disturbance two-order stochastic process optimization model is constructed based on the travel speed; the vertical line deviation disturbance two-order stochastic process optimization model includes a two-order stochastic process model and a differentiator; the vertical line deviation disturbance two-order stochastic process optimization model The output of the model is the output result after the output of the two-order stochastic process model passes through the differentiator;

Calculate the vertical deviation disturbance component based on the vertical deviation disturbance two-order stochastic process optimization model; the vertical deviation disturbance component includes the vertical deviation disturbance north-south component and the vertical deviation disturbance east-west component;

The variation of groundwater storage in the area where the measurement carrier is located is obtained by inversion from the vertical deviation disturbance component.

Optionally, the vertical line deviation disturbance two-order stochastic process optimization model is:

为x₁(t)的一阶导数，

为x₁(t)的二阶导数，ω₀为中心频率，

为阻尼系数，

为q(t)的一阶导数，q(t)为高斯白噪声，ω₀＝2πV/λ₀，V为行驶速度，λ₀为中心波长。x ₁ (t) is the output of the two-order stochastic process optimization model perturbed by vertical deviation,

is the first derivative of x ₁ (t),

is the second derivative of x ₁ (t), ω ₀ is the center frequency,

is the damping coefficient,

is the first derivative of q(t), q(t) is white Gaussian noise, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength.

Optionally, the calculation formula of the vertical line deviation disturbance component is:

为阻尼系数，ω₀＝2πV/λ₀，V为行驶速度，λ₀为中心波长。δξ(t) is the north-south component of the vertical deviation disturbance, δη(t) is the east-west component of the vertical deviation disturbance, x _ξ (t) is the intermediate variable in the north-south direction, x _η (t) is the intermediate variable in the east-west direction, x _ξ The derivative of (t) is δξ(t), the derivative of x _η (t) is δη(t), q _η (t) is the process noise in the east-west direction, q _ξ (t) is the process noise in the north-south direction, ω ₀ is the center frequency,

is the damping coefficient, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength.

Optionally, the inversion of the variation of groundwater storage in the area where the measurement carrier is located from the vertical line deviation disturbance component specifically includes:

Inverting the water storage variation corresponding to the vertical deviation disturbance component from the vertical deviation disturbance component;

The variation of groundwater storage in the area where the measurement carrier is located is calculated from the variation of water storage, the variation of snow water equivalent and the variation of water content in the soil corresponding to the perturbation component of the vertical line deviation.

Optionally, the calculation formula of the groundwater storage change is:

ΔGN=ΔTNS-ΔSN-ΔSNE;

ΔGN is the change in groundwater storage, ΔTNS is the change in water storage corresponding to the vertical deviation disturbance component, ΔSN is the change in snow water equivalent, and ΔSNE is the change in soil water content.

The present invention also provides a groundwater storage monitoring system based on vertical deviation disturbance, comprising:

The data acquisition module is used to acquire the traveling speed of the measurement carrier;

a model building module for constructing a vertical deviation perturbation two-order stochastic process optimization model based on the travel speed; the vertical deviation perturbation two-order stochastic process optimization model includes a two-order stochastic process model and a differentiator; the vertical deviation The output of the perturbed two-order stochastic process optimization model is the output result after the output of the two-order stochastic process model passes through the differentiator;

The vertical deviation disturbance calculation module is used to calculate the vertical deviation disturbance component based on the vertical deviation disturbance two-order stochastic process optimization model; the vertical deviation disturbance component includes the vertical deviation disturbance north-south component and the vertical deviation disturbance east-west component ;

The inversion module is configured to invert the variation of groundwater storage in the area where the measurement carrier is located from the vertical deviation disturbance component.

Optionally, the vertical line deviation disturbance two-order stochastic process optimization model in the model building module is:

为x₁(t)的一阶导数，

为x₁(t)的二阶导数，ω₀为中心频率，

为阻尼系数，

为q(t)的一阶导数，q(t)为高斯白噪声，ω₀＝2πV/λ₀，V为行驶速度，λ₀为中心波长。x ₁ (t) is the output of the two-order stochastic process optimization model perturbed by vertical deviation,

is the first derivative of x ₁ (t),

is the second derivative of x ₁ (t), ω ₀ is the center frequency,

is the damping coefficient,

is the first derivative of q(t), q(t) is white Gaussian noise, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength.

Optionally, the calculation formula of the vertical deviation disturbance component in the vertical deviation disturbance calculation module is:

为阻尼系数，ω₀＝2πV/λ₀，V为行驶速度，λ₀为中心波长。δξ(t) is the north-south component of the vertical deviation disturbance, δη(t) is the east-west component of the vertical deviation disturbance, x _ξ (t) is the intermediate variable in the north-south direction, x _η (t) is the intermediate variable in the east-west direction, x _ξ The derivative of (t) is δξ(t), the derivative of x _η (t) is δη(t), q _η (t) is the process noise in the east-west direction, q _ξ (t) is the process noise in the north-south direction, ω ₀ is the center frequency,

is the damping coefficient, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength.

Optionally, the inversion module specifically includes:

a first inversion unit, configured to invert the variation of water storage corresponding to the vertical deviation disturbance component from the vertical deviation disturbance component;

A calculation unit, configured to calculate the variation of groundwater storage in the area where the measurement carrier is located from the variation of water storage, the variation of snow water equivalent, and the variation of water content in the soil corresponding to the vertical line deviation disturbance component.

Optionally, the calculation formula of the groundwater storage variation in the calculation unit is:

ΔGN=ΔTNS-ΔSN-ΔSNE;

ΔGN is the change in groundwater storage, ΔTNS is the change in water storage corresponding to the vertical deviation disturbance component, ΔSN is the change in snow water equivalent, and ΔSNE is the change in soil water content.

Compared with the prior art, the beneficial effects of the present invention are:

The invention proposes a groundwater storage monitoring method and system based on vertical deviation disturbance, constructs a vertical deviation disturbance two-order stochastic process optimization model including a two-order stochastic process model and a differentiator; The optimization model calculates the perturbation component of vertical line deviation, so as to obtain the variation of groundwater storage in the area where the measurement carrier is located. The present invention passes the output of the random process model through a differentiator, which can reduce the gain of the high-frequency disturbance part of the vertical deviation in the low-frequency region, so that under the action of the differentiator, the high-frequency vertical deviation perturbation has a stronger effect in the low-frequency region. It improves the approximation accuracy of vertical line deviation disturbance, solves the problem of inability to obtain large-scale and large-scale groundwater data, and improves the monitoring efficiency of groundwater storage.

Description of drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the accompanying drawings required in the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some of the present invention. In the embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without creative labor.

1 is a flowchart of a method for monitoring groundwater storage based on vertical line deviation disturbance provided by an embodiment of the present invention;

2 is a schematic diagram of a vertical line deviation perturbing a north-south component provided by an embodiment of the present invention;

3 is a schematic diagram of a vertical line deviation perturbing an east-west component provided by an embodiment of the present invention;

FIG. 4 is a schematic structural diagram of a groundwater storage monitoring system based on vertical line deviation disturbance provided by an embodiment of the present invention.

Detailed ways

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

In order to make the above objects, features and advantages of the present invention more clearly understood, the present invention will be described in further detail below with reference to the accompanying drawings and specific embodiments.

Excessive exploitation of groundwater can cause problems such as water depletion, land subsidence, and seawater intrusion. Most of the traditional groundwater storage monitoring methods are realized by acquiring the data of fixed monitoring points of water resources on the earth's surface, which will be affected by the number and distribution of monitoring stations. Changes in groundwater storage will cause changes in local area mass, and changes in Earth's mass distribution will lead to changes in the Earth's gravitational field. The vertical deviation disturbance is the high frequency component of the gravity field, and the vertical deviation disturbance has a strong correlation with groundwater storage. Therefore, groundwater storage can be inverted by perturbation of vertical line deviation. This method can efficiently monitor the groundwater level and solve the problem that large-scale and large-scale groundwater data cannot be obtained.

其功率谱密度在低频区域的增益较大致使高频垂线偏差扰动的模型与相关误差存在严重的耦合现象。为了抑制该耦合效应，使短波垂线偏差扰动模型的功率谱密度在低频区域内具有强的衰减特性，减小低频增益，本实施例将两阶随机过程模型的输出x(t)通过一个微分器，得到随机过程x₁(t)。这样减小了垂线偏差高频扰动部分在低频区域的增益。在微分器的作用下，高频垂线偏差扰动在低频区域内具有更强烈的衰减特性，提高了垂线偏差扰动的逼近精度。The two-order stochastic process model can be written as:

The gain of its power spectral density in the low-frequency region is large, which leads to a serious coupling phenomenon between the model perturbed by the high-frequency vertical line deviation and the related errors. In order to suppress the coupling effect, make the power spectral density of the short-wave vertical line deviation disturbance model have strong attenuation characteristics in the low-frequency region, and reduce the low-frequency gain, in this embodiment, the output x(t) of the two-order stochastic process model is passed through a differential , to obtain a random process x ₁ (t). This reduces the gain of the high frequency disturbance part of the vertical deviation in the low frequency region. Under the action of the differentiator, the high frequency vertical deviation disturbance has stronger attenuation characteristics in the low frequency region, which improves the approximation accuracy of the vertical deviation disturbance.

FIG. 1 is a flowchart of a groundwater storage monitoring method based on vertical line deviation disturbance provided by an embodiment of the present invention. Referring to Fig. 1, the groundwater storage monitoring method based on vertical deviation disturbance in this embodiment includes:

Step 101: Obtain the traveling speed of the measurement carrier. The measurement carrier is a measurement vehicle or a measurement ship.

Step 102 : constructing a vertical deviation perturbation two-order stochastic process optimization model based on the traveling speed; the vertical deviation perturbation two-order stochastic process optimization model includes a two-order stochastic process model and a differentiator.

The output of the vertical line deviation disturbance two-order stochastic process optimization model is the output result after the output of the two-order stochastic process model passes through the differentiator.

Step 103: Calculate the vertical deviation disturbance component based on the vertical deviation disturbance two-order stochastic process optimization model; the vertical deviation disturbance component includes the vertical deviation disturbance north-south component and the vertical deviation disturbance east-west component.

Step 104: Inverting the variation of groundwater storage in the area where the measurement carrier is located from the vertical deviation disturbance component.

Wherein, in step 102, the two-order stochastic process model is:

The differential equation of the vertical deviation perturbation two-order stochastic process is:

Simplifying the differential equation of the above-mentioned vertical line deviation perturbation two-order stochastic process can be obtained:

为x(t)的一阶导数，

为x(t)的二阶导数，x₁(t)为垂线偏差扰动两阶随机过程优化模型的输出，

为x₁(t)的一阶导数，

为x₁(t)的二阶导数，ω₀为中心频率，

为阻尼系数，

为q(t)的一阶导数，q(t)为高斯白噪声，ω₀＝2πV/λ₀，V为行驶速度，λ₀为中心波长。

和q(t)可以由行驶速度和测量环境条件确定。The simplified differential equation of the two-order stochastic process perturbed by vertical deviation is the optimization model of the two-order stochastic process perturbed by vertical deviation. where x(t) is a two-order stochastic process model,

is the first derivative of x(t),

is the second derivative of x(t), x ₁ (t) is the output of the optimization model of the vertical deviation perturbation two-order stochastic process,

is the first derivative of x ₁ (t),

is the second derivative of x ₁ (t), ω ₀ is the center frequency,

is the damping coefficient,

is the first derivative of q(t), q(t) is white Gaussian noise, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength.

and q(t) can be determined from the driving speed and the measured environmental conditions.

In order to facilitate subsequent calculations, the vertical line deviation disturbance two-order stochastic process optimization model can be written in the form of state equation:

Wherein, in step 103, the calculation formula of the vertical deviation disturbance component is:

为阻尼系数，ω₀＝2πV/λ₀，V为行驶速度，λ₀为中心波长。δξ(t) is the north-south component of the vertical deviation disturbance (the component of the vertical deviation disturbance on the meridian plane, also known as the vertical deviation disturbance meridian component), and δη(t) is the east-west component of the vertical deviation disturbance (the vertical deviation disturbance The component on the 卯寉 plane, also known as the vertical deviation disturbance 卯寉 component), x _ξ (t) is the intermediate variable in the north-south direction, x _η (t) is the intermediate variable in the east-west direction, and x _ξ (t) is the intermediate variable. The derivative is δξ(t), the derivative of x _η (t) is δη(t), q _η (t) is the process noise in the east-west direction, q _ξ (t) is the process noise in the north-south direction, ω ₀ is the center frequency,

is the damping coefficient, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength.

Among them, the vertical line deviation disturbance component in step 104 contains rich high-frequency information of the spatial variation of the gravitational field, and is very sensitive to changes in the mass of the earth's surface. The change of groundwater storage will inevitably lead to the change of the quality of the local area, and the change of the earth's mass distribution will lead to the change of the earth's gravitational field. Therefore, groundwater storage can be inverted from the perturbation component of vertical deviation. The local mass change will lead to the change of the gravity field of the point to be measured. The local surface area belongs to the high frequency category of the gravity field. The vertical deviation disturbance component is the high frequency component of the gravity field, which is related to the regional mass change, and can be determined by the vertical deviation disturbance component. The change in groundwater storage is obtained by inversion. Step 104 specifically includes:

1) Inverting the water storage variation corresponding to the vertical deviation disturbance component from the vertical deviation disturbance component.

2) Calculate the variation of groundwater storage in the area where the measurement carrier is located from the variation of water storage, the variation of snow water equivalent, and the variation of water content in the soil corresponding to the vertical line deviation disturbance component. The calculation formula of the groundwater storage change is:

ΔGN=ΔTNS-ΔSN-ΔSNE;

ΔGN is the change in groundwater storage, ΔTNS is the change in water storage corresponding to the vertical deviation disturbance component, ΔSN is the change in snow water equivalent, and ΔSNE is the change in soil water content.

After step 102, this embodiment further includes: performing Laplace transformation on the vertical deviation disturbance two-order stochastic process optimization model to obtain the transfer function of the stochastic process x ₁ (t) and the corresponding power spectral density.

The transfer function of x ₁ (t) is:

where s represents the complex frequency.

The power spectral density is obtained by the transfer function, and the corresponding power spectral density is:

为q(t)的方差。where ω is the frequency, j is the imaginary unit,

is the variance of q(t).

By comparing and analyzing the output x(t) of the two-order stochastic process model and the power spectral density distribution of the two-order stochastic process optimization model x ₁ (t) perturbed by vertical line deviation, it is concluded that x ₁ (t) is lower than x(t) at low frequencies. The attenuation is stronger in the region. The vertical deviation perturbation two-order stochastic process optimization model x ₁ (t) reduces the gain of the high frequency perturbation part of the vertical deviation in the low frequency region. Under the action of the differentiator, the high-frequency vertical deviation disturbance has stronger attenuation characteristics in the low-frequency region, which improves the approximation accuracy of the vertical deviation, and solves the problem that large-scale and large-scale groundwater data cannot be obtained, thereby improving the Monitoring efficiency of groundwater storage.

In order to verify the effectiveness of the above groundwater storage monitoring method based on vertical deviation disturbance, the above method is used to approximate the high frequency disturbance components of vertical deviation δη(t) and δξ(t). It can be seen from Fig. 2 and Fig. 3 that δη(t) and δξ(t) fluctuate randomly around 0, the approximation accuracy is within 2", and the accuracy is high. Therefore, the reliability of the above method in this embodiment is relatively high. It shows that the above method can effectively measure the high-frequency disturbance of vertical deviation, and the accurate high-frequency disturbance of vertical deviation can efficiently invert groundwater storage.

The present invention also provides a groundwater storage monitoring system based on vertical deviation disturbance. FIG. 4 is a schematic structural diagram of the groundwater storage monitoring system based on vertical deviation disturbance provided by an embodiment of the present invention. Referring to Figure 4, the groundwater storage monitoring system based on vertical deviation disturbance includes:

The data acquisition module 201 is used for acquiring the traveling speed of the measurement carrier.

A model building module 202, configured to build a vertical deviation disturbance two-order stochastic process optimization model based on the traveling speed; the vertical deviation disturbance two-order stochastic process optimization model includes a two-order stochastic process model and a differentiator; the vertical line The output of the deviation disturbance two-order stochastic process optimization model is the output result after the output of the two-order stochastic process model passes through the differentiator.

The vertical deviation disturbance calculation module 203 is configured to calculate the vertical deviation disturbance component based on the vertical deviation disturbance two-order stochastic process optimization model; the vertical deviation disturbance component includes the vertical deviation disturbance north-south component and the vertical deviation disturbance east-west component weight.

The inversion module 204 is configured to invert the variation of groundwater storage in the area where the measurement carrier is located from the vertical deviation disturbance component.

As an optional implementation manner, the vertical line deviation disturbance two-order stochastic process optimization model in the model building module 202 is:

为x₁(t)的一阶导数，

为x₁(t)的二阶导数，ω₀为中心频率，

为阻尼系数，

为q(t)的一阶导数，q(t)为高斯白噪声，ω₀＝2πV/λ₀，V为行驶速度，λ₀为中心波长。x ₁ (t) is the output of the two-order stochastic process optimization model perturbed by vertical deviation,

is the first derivative of x ₁ (t),

is the second derivative of x ₁ (t), ω ₀ is the center frequency,

is the damping coefficient,

is the first derivative of q(t), q(t) is white Gaussian noise, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength.

As an optional implementation manner, the calculation formula of the vertical deviation disturbance component in the vertical deviation disturbance calculation module 203 is:

为阻尼系数，ω₀＝2πV/λ₀，V为行驶速度，λ₀为中心波长。δξ(t) is the north-south component of the vertical deviation disturbance, δη(t) is the east-west component of the vertical deviation disturbance, x _ξ (t) is the intermediate variable in the north-south direction, x _η (t) is the intermediate variable in the east-west direction, x _ξ The derivative of (t) is δξ(t), the derivative of x _η (t) is δη(t), q _η (t) is the process noise in the east-west direction, q _ξ (t) is the process noise in the north-south direction, ω ₀ is the center frequency,

is the damping coefficient, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength.

As an optional implementation manner, the inversion module 204 specifically includes:

The first inversion unit is configured to invert the variation of water storage corresponding to the vertical deviation disturbance component from the vertical deviation disturbance component.

A calculation unit, configured to calculate the variation of groundwater storage in the area where the measurement carrier is located from the variation of water storage, the variation of snow water equivalent, and the variation of water content in the soil corresponding to the vertical line deviation disturbance component.

As an optional implementation manner, the calculation formula of the groundwater storage variation in the calculation unit is:

ΔGN=ΔTNS-ΔSN-ΔSNE;

ΔGN is the change in groundwater storage, ΔTNS is the change in water storage corresponding to the vertical deviation disturbance component, ΔSN is the change in snow water equivalent, and ΔSNE is the change in soil water content.

The groundwater storage monitoring system based on the vertical deviation disturbance provided by this embodiment can reduce the gain of the high frequency disturbance part of the vertical deviation in the low frequency region, so that the high frequency vertical deviation disturbance has a stronger attenuation characteristic in the low frequency region, Improve the approximation accuracy of vertical line deviation disturbance and efficiently invert groundwater storage.

The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments, and the same and similar parts between the various embodiments can be referred to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant part can be referred to the description of the method.

In this paper, specific examples are used to illustrate the principles and implementations of the present invention. The descriptions of the above embodiments are only used to help understand the methods and core ideas of the present invention; meanwhile, for those skilled in the art, according to the present invention There will be changes in the specific implementation and application scope. In conclusion, the contents of this specification should not be construed as limiting the present invention.

Claims (10)

1. a groundwater storage monitoring method based on vertical deviation disturbance, is characterized in that, comprises: Obtain the traveling speed of the measurement carrier; A vertical line deviation disturbance two-order stochastic process optimization model is constructed based on the travel speed; the vertical line deviation disturbance two-order stochastic process optimization model includes a two-order stochastic process model and a differentiator; the vertical line deviation disturbance two-order stochastic process optimization model The output of the model is the output result after the output of the two-order stochastic process model passes through the differentiator; Calculate the vertical deviation disturbance component based on the vertical deviation disturbance two-order stochastic process optimization model; the vertical deviation disturbance component includes the vertical deviation disturbance north-south component and the vertical deviation disturbance east-west component; The variation of groundwater storage in the area where the measurement carrier is located is obtained by inversion from the vertical deviation disturbance component. 2. a kind of groundwater storage monitoring method based on vertical line deviation disturbance according to claim 1, is characterized in that, described vertical line deviation disturbance two-order stochastic process optimization model is:

x ₁ (t) is the output of the two-order stochastic process optimization model perturbed by vertical deviation,

is the first derivative of x ₁ (t),

is the second derivative of x ₁ (t), ω ₀ is the center frequency,

is the damping coefficient,

is the first derivative of q(t), q(t) is white Gaussian noise, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength. 3. a kind of groundwater storage monitoring method based on vertical deviation disturbance according to claim 1, is characterized in that, the calculation formula of described vertical deviation disturbance component is:

δξ(t) is the north-south component of the vertical deviation disturbance, δη(t) is the east-west component of the vertical deviation disturbance, x _ξ (t) is the intermediate variable in the north-south direction, x _η (t) is the intermediate variable in the east-west direction, x _ξ The derivative of (t) is δξ(t), the derivative of x _η (t) is δη(t), q _η (t) is the process noise in the east-west direction, q _ξ (t) is the process noise in the north-south direction, ω ₀ is the center frequency,

is the damping coefficient, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength. 4. The method for monitoring groundwater storage based on vertical deviation disturbance according to claim 1, wherein the groundwater storage variation in the area where the measurement carrier is located is obtained by inversion from the vertical deviation disturbance component amount, including: Inverting the water storage variation corresponding to the vertical deviation disturbance component from the vertical deviation disturbance component; The variation of groundwater storage in the area where the measurement carrier is located is calculated from the variation of water storage, the variation of snow water equivalent and the variation of water content in the soil corresponding to the perturbation component of the vertical line deviation. 5. a kind of groundwater storage monitoring method based on vertical line deviation disturbance according to claim 4, is characterized in that, the calculation formula of described groundwater storage variation is: ΔGN=ΔTNS-ΔSN-ΔSNE; ΔGN is the change in groundwater storage, ΔTNS is the change in water storage corresponding to the vertical deviation disturbance component, ΔSN is the change in snow water equivalent, and ΔSNE is the change in soil water content. 6. A groundwater storage monitoring system based on vertical deviation disturbance, characterized in that, comprising: The data acquisition module is used to acquire the traveling speed of the measurement carrier; a model building module for constructing a vertical deviation perturbation two-order stochastic process optimization model based on the travel speed; the vertical deviation perturbation two-order stochastic process optimization model includes a two-order stochastic process model and a differentiator; the vertical deviation The output of the perturbed two-order stochastic process optimization model is the output result after the output of the two-order stochastic process model passes through the differentiator; The vertical deviation disturbance calculation module is used to calculate the vertical deviation disturbance component based on the vertical deviation disturbance two-order stochastic process optimization model; the vertical deviation disturbance component includes the vertical deviation disturbance north-south component and the vertical deviation disturbance east-west component ; The inversion module is configured to invert the variation of groundwater storage in the area where the measurement carrier is located from the vertical deviation disturbance component. 7. a kind of groundwater storage monitoring system based on vertical deviation disturbance according to claim 6, is characterized in that, described vertical deviation disturbance two-order stochastic process optimization model in described model building module is:

x ₁ (t) is the output of the two-order stochastic process optimization model perturbed by vertical deviation,

is the first derivative of x ₁ (t),

is the second derivative of x ₁ (t), ω ₀ is the center frequency,

is the damping coefficient,

is the first derivative of q(t), q(t) is white Gaussian noise, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength. 8 . The groundwater storage monitoring system based on vertical deviation disturbance according to claim 6 , wherein the calculation formula of the vertical deviation disturbance component in the vertical deviation disturbance calculation module is: 8 .

δξ(t) is the north-south component of the vertical deviation disturbance, δη(t) is the east-west component of the vertical deviation disturbance, x _ξ (t) is the intermediate variable in the north-south direction, x _η (t) is the intermediate variable in the east-west direction, x _ξ The derivative of (t) is δξ(t), the derivative of x _η (t) is δη(t), q _η (t) is the process noise in the east-west direction, q _ξ (t) is the process noise in the north-south direction, ω ₀ is the center frequency,

is the damping coefficient, ω ₀ =2πV/λ ₀ , V is the traveling speed, and λ ₀ is the center wavelength. 9. A groundwater storage monitoring system based on vertical line deviation disturbance according to claim 6, wherein the inversion module specifically comprises: a first inversion unit, configured to invert the variation of water storage corresponding to the vertical deviation disturbance component from the vertical deviation disturbance component; A calculation unit, configured to calculate the variation of groundwater storage in the area where the measurement carrier is located from the variation of water storage, the variation of snow water equivalent, and the variation of water content in the soil corresponding to the vertical line deviation disturbance component. 10. A groundwater storage monitoring system based on vertical line deviation disturbance according to claim 9, wherein the calculation formula of the groundwater storage variation in the calculation unit is: ΔGN=ΔTNS-ΔSN-ΔSNE; ΔGN is the change in groundwater storage, ΔTNS is the change in water storage corresponding to the vertical deviation disturbance component, ΔSN is the change in snow water equivalent, and ΔSNE is the change in soil water content.

Images