Method for determining groundwater level change amount based on gravity gradient, electronic device and medium

Abstract

The application provides a method for determining a change amount of a groundwater level based on a gravity gradient, an electronic device and a medium, and can be applied to the technical field of hydrographic geography. The method comprises the following steps: obtaining gravity gradient change amounts of a plurality of groundwater level measurement points; based on the correlation between the gravity gradient and the water level, obtaining the groundwater level change amounts of the plurality of groundwater level measurement points according to the three-dimensional underground network associated with the plurality of groundwater level measurement points and the gravity gradient change amounts of the plurality of groundwater level measurement points; the gravity gradient change amount of the gravity gradient measurement point and the correlation between the gravity gradient and the water level are obtained according to the attribute information of a plurality of grids included in the three-dimensional underground network associated with the gravity gradient measurement point; and the three-dimensional underground network associated with the surface measurement point is obtained by grid division on a three-dimensional space determined according to a preset gravity gradient resolution, a preset groundwater level resolution and a soil saturated water content rate associated with the surface measurement point.

Description

Technical Field

[0001] This application relates to the field of hydrogeography, and more specifically to a method, electronic device and medium for determining groundwater level changes based on gravity gradient. Background Technology

[0002] Urban building safety is closely related to groundwater changes, and the safety hazards of urban buildings caused by groundwater level fluctuations are becoming increasingly prominent. Therefore, groundwater level monitoring can capture dynamic trends in groundwater levels, providing a basis for early warning of urban building safety, foundation protection design, and risk prevention and control. Summary of the Invention

[0003] In view of the above problems, embodiments of this application provide a method, electronic device and medium for determining groundwater level changes based on gravity gradient.

[0004] According to a first aspect of the embodiments of this application, a method for determining groundwater level changes based on gravity gradient is provided, comprising: acquiring the gravity gradient changes of multiple groundwater level test points; based on the correlation between gravity gradient and water level, and according to a three-dimensional underground network associated with the multiple groundwater level test points and the gravity gradient changes of the multiple groundwater level test points, obtaining the groundwater level changes of each of the multiple groundwater level test points; wherein, the gravity gradient changes of the gravity gradient test points and the correlation between gravity gradient and water level are determined based on the correlation between gravity gradient and water level of the gravity gradient test points. The three-dimensional underground network is obtained by dividing the three-dimensional space into grids based on the attribute information of multiple grids. The three-dimensional underground network associated with the surface measurement points is obtained by dividing the three-dimensional space into grids according to the preset gravity gradient resolution, preset groundwater level resolution, and soil saturated moisture content associated with the surface measurement points. The surface measurement points include gravity gradient measurement points and groundwater level measurement points. The attribute information of the grids includes the grid position function and the preset moisture content change. The position function represents the change in gravity gradient at the gravity gradient measurement point caused by a unit density change in the grid.

[0005] According to an embodiment of this application, based on the correlation between gravity gradient and water level, and based on the three-dimensional underground network associated with multiple groundwater level measurement points and the gravity gradient change of the multiple groundwater level measurement points, the groundwater level change of each of the multiple groundwater level measurement points is obtained. This includes: for any groundwater level measurement point among the multiple groundwater level measurement points, constructing a three-dimensional underground network water content equation based on the gravity gradient change of the groundwater level measurement point and the three-dimensional underground network associated with the groundwater level measurement point; and solving the three-dimensional underground network water content equation set composed of the multiple three-dimensional underground network water content equations based on the correlation between gravity gradient and water level to obtain the groundwater level change of each of the multiple groundwater level measurement points.

[0006] According to an embodiment of this application, based on the correlation between gravity gradient and water level, the three-dimensional underground network water content equation system composed of multiple three-dimensional underground network water content equations is solved to obtain the groundwater level change of each of the multiple groundwater level test points. This includes: determining the spatial distribution range of groundwater change associated with each of the multiple groundwater level test points based on the correlation between gravity gradient and water level; and solving the three-dimensional underground grid water content equation system based on the spatial distribution range of groundwater change to obtain the groundwater level change of each of the multiple groundwater level test points.

[0007] According to an embodiment of this application, the method for determining the groundwater level change based on gravity gradient further includes: for any groundwater level test point among multiple groundwater level test points, obtaining the groundwater level of the groundwater level test point based on the groundwater level change of the groundwater level test point and the groundwater level of a known survey point, wherein the groundwater level represents the depth of the water level relative to a preset reference surface; and correcting the groundwater level of the groundwater level test point based on the preset reference surface to obtain the target groundwater level of the groundwater level test point, wherein the target groundwater level represents the depth of the water level relative to the ground surface.

[0008] According to an embodiment of this application, the three-dimensional underground network associated with the surface measurement points is obtained in the following manner: determining the vertical measurement range associated with the surface measurement points based on a preset gravity gradient resolution, a preset groundwater level resolution, and the soil saturation moisture content associated with the surface measurement points; determining the horizontal measurement range associated with the surface measurement points based on the vertical measurement range; and dividing the three-dimensional space formed by the vertical and horizontal measurement ranges into a grid to obtain the three-dimensional underground network associated with the surface measurement points.

[0009] According to an embodiment of this application, determining the vertical measurement range associated with a surface measurement point based on a preset gravity gradient resolution, a preset groundwater level resolution, and the soil saturated moisture content associated with the surface measurement point includes: determining the vertical distance associated with the surface measurement point based on the preset gravity gradient resolution and the soil saturated moisture content associated with the surface measurement point, wherein the vertical distance characterizes the upper limit of the vertical downward detection depth of the surface measurement point; and determining the vertical measurement range based on the multiple relationship between the vertical distance and the preset groundwater level resolution.

[0010] According to an embodiment of this application, a three-dimensional space composed of vertical and horizontal measurement ranges is divided into grids to obtain a three-dimensional underground network associated with surface measurement points. This includes: dividing the vertical measurement range into equal-interval layers based on a preset groundwater level resolution to obtain multiple horizontal layers; for any horizontal layer, starting from the center point of the horizontal layer, expanding horizontally according to preset reference grid units to obtain the actual gravity gradient generated by at least one reference grid within the horizontal layer on the surface measurement points, wherein the reference grid is constructed with side lengths that are integer multiples of the preset reference grid units and heights that are the distance between the horizontal layer and the surface measurement points; determining at least one grid for the horizontal layer when the actual gravity gradient and the gravity gradient resolution satisfy a preset multiple relationship; and stitching together at least one grid from each of the multiple horizontal layers to obtain a three-dimensional underground network associated with the surface measurement points.

[0011] According to an embodiment of this application, based on the attribute information of each of the multiple grids included in the three-dimensional underground network associated with the gravity gradient test point, the gravity gradient change of the gravity gradient test point and the correlation between the gravity gradient and the water level are obtained, including: for any grid among the multiple grids included in the three-dimensional underground network associated with the gravity gradient test point, determining the gravity gradient change component of the grid based on the product between the grid's position function and a preset water content change; accumulating the gravity gradient change components of each of the multiple grids to obtain the gravity gradient change of the gravity gradient test point; and determining the correlation between the gravity gradient and the water level based on the gravity gradient change of the gravity gradient test point and the preset water content change of each of the multiple grids included in the three-dimensional underground network associated with the gravity gradient test point.

[0012] A second aspect of this application provides an apparatus for determining groundwater level changes based on gravity gradient, comprising: an acquisition module for acquiring gravity gradient changes of multiple groundwater level test points; and an obtaining module for obtaining groundwater level changes of each of the multiple groundwater level test points based on the correlation between gravity gradient and water level, according to a three-dimensional underground network associated with the multiple groundwater level test points and the gravity gradient changes of the multiple groundwater level test points; wherein the gravity gradient changes of the gravity gradient test points and the correlation between gravity gradient and water level are determined based on the correlation between gravity gradient and water level. The three-dimensional underground network associated with the gradient measurement points is obtained by dividing the three-dimensional space into grids based on the preset gravity gradient resolution, preset groundwater level resolution, and soil saturated moisture content associated with the surface measurement points. The surface measurement points include gravity gradient measurement points and groundwater level measurement points. The grid attribute information includes the grid's position function and preset moisture content change. The position function characterizes the change in gravity gradient at the gravity gradient measurement points when the grid undergoes a unit density change.

[0013] A third aspect of this application provides an electronic device, including: one or more processors; and a memory for storing one or more computer programs, wherein the one or more processors execute the one or more computer programs to implement the steps of the method described above.

[0014] A fourth aspect of this application provides a computer-readable storage medium having a computer program or instructions stored thereon, which, when executed by a processor, implement the steps of the above-described method.

[0015] A fifth aspect of this application also provides a computer program product, including a computer program or instructions that, when executed by a processor, implement the steps of the above-described method.

[0016] By establishing a three-dimensional underground network associated with surface measurement points and combining the gravity gradient changes at the measurement points, the groundwater level changes at multiple groundwater level measurement points can be obtained. Based on this, by adaptively meshing the three-dimensional space determined by the preset gravity gradient resolution, preset groundwater level resolution, and soil saturation moisture content associated with the surface measurement points, a three-dimensional underground network associated with the surface measurement points is obtained. This improves the accuracy of the correlation between gravity gradient and water level, thereby enhancing data processing precision and efficiency. Consequently, the accuracy of the groundwater level changes at multiple groundwater level measurement points is improved, providing technical support for groundwater dynamic detection and inversion. Attached Figure Description

[0017] The above-mentioned contents, other objects, features and advantages of this application will become clearer from the following description of embodiments with reference to the accompanying drawings, in which:

[0018] Figure 1 The diagram illustrates an application scenario of the method for determining groundwater level changes based on gravity gradient according to an embodiment of this application.

[0019] Figure 2 A flowchart illustrating a method for determining groundwater level changes based on gravity gradient according to an embodiment of this application is shown schematically.

[0020] Figure 3 A flowchart illustrating the construction of a three-dimensional underground network associated with surface measurement points according to an embodiment of this application is shown.

[0021] Figure 4 A schematic diagram illustrating a three-dimensional underground network associated with surface measurement points according to an embodiment of this application is shown.

[0022] Figure 5 This schematically illustrates a flowchart of the forward modeling process for obtaining the change in gravity gradient at the point to be measured according to an embodiment of this application.

[0023] Figure 6 This illustration schematically shows a flowchart of the inversion process to obtain the groundwater level at the point to be measured according to an embodiment of this application;

[0024] Figure 7 A schematic diagram illustrates the structure of a device for determining groundwater level changes based on gravity gradient according to an embodiment of this application; and

[0025] Figure 8 A block diagram of an electronic device suitable for implementing a method for determining groundwater level changes based on gravity gradients, according to an embodiment of this application, is illustrated schematically. Detailed Implementation

[0026] The embodiments of this application will now be described with reference to the accompanying drawings. However, it should be understood that these descriptions are exemplary only and are not intended to limit the scope of this application. In the following detailed description, numerous specific details are set forth to provide a thorough understanding of the embodiments of this application for ease of explanation. However, it will be apparent that one or more embodiments may be implemented without these specific details. Furthermore, descriptions of well-known structures and technologies are omitted in the following description to avoid unnecessarily obscuring the concepts of this application.

[0027] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the scope of this application. The terms “comprising,” “including,” etc., as used herein indicate the presence of the stated features, steps, operations, and / or components, but do not exclude the presence or addition of one or more other features, steps, operations, or components.

[0028] All terms used herein (including technical and scientific terms) have the meanings commonly understood by those skilled in the art, unless otherwise defined. It should be noted that the terms used herein are to be interpreted in a manner consistent with the context of this specification, and not in an idealized or overly rigid way.

[0029] When using expressions such as "at least one of A, B and C", they should generally be interpreted in accordance with the meaning that is commonly understood by those skilled in the art (e.g., "a system having at least one of A, B and C" should include, but is not limited to, a system having A alone, a system having B alone, a system having C alone, a system having A and B, a system having A and C, a system having B and C, and / or a system having A, B and C, etc.).

[0030] Due to changes in groundwater levels, current urban development faces the safety issue of buildings "sinking" or "floating." According to statistics, nearly half of all cities are currently experiencing geological subsidence. Furthermore, especially in coastal and riverside areas, urban building safety is closely related to groundwater fluctuations. The safety hazards caused by groundwater level fluctuations are becoming increasingly prominent. Therefore, monitoring groundwater levels can capture dynamic trends and provide a basis for early warning systems, foundation protection design, and risk control in urban buildings.

[0031] However, in related technologies, the detection efficiency and accuracy of methods using gravity gradient technology to detect groundwater level changes are relatively low. In view of this, embodiments of this application provide a method for determining groundwater level changes based on gravity gradients, including: acquiring the gravity gradient changes of multiple groundwater level measurement points; based on the correlation between gravity gradients and water levels, and according to a three-dimensional underground network associated with the multiple groundwater level measurement points and the gravity gradient changes of the multiple groundwater level measurement points, obtaining the groundwater level changes of each of the multiple groundwater level measurement points; wherein, the gravity gradient changes of the gravity gradient measurement points and the correlation between gravity gradients and water levels are determined based on the correlation between gravity gradient measurement points and groundwater levels. The three-dimensional underground network is obtained by dividing the three-dimensional space into grids based on the attribute information of multiple grids. The three-dimensional underground network associated with the surface measurement points is obtained by dividing the three-dimensional space into grids according to the preset gravity gradient resolution, the preset groundwater level resolution, and the soil saturated moisture content associated with the surface measurement points. The surface measurement points include gravity gradient measurement points and groundwater level measurement points. The attribute information of the grids includes the grid position function and the preset moisture content change. The position function represents the change in gravity gradient at the gravity gradient measurement point caused by a unit density change in the grid.

[0032] Figure 1 The illustration shows an application scenario of the method for determining groundwater level changes based on gravity gradient according to an embodiment of this application.

[0033] like Figure 1 As shown, application scenario 100 according to this embodiment may include a terminal device 101, a second terminal device 102, a third terminal device 103, a network 104, a server 105, and a gravity gradiometer 106. The network 104 serves as a medium for providing communication links between the first terminal device 101, the second terminal device 102, the third terminal device 103, the server 105, and the gravity gradiometer 106. The network 104 may include various connection types, such as wired or wireless communication links, or fiber optic cables, etc.

[0034] Users can use the first terminal device 101, the second terminal device 102, and the third terminal device 103 to interact with the server 105 and the gravity gradiometer 106 via the network 104 to receive or send messages, etc. Various communication client applications can be installed on the first terminal device 101, the second terminal device 102, and the third terminal device 103, such as shopping applications, web browser applications, search applications, instant messaging tools, email clients, social media platform software, etc. (for example only).

[0035] The first terminal device 101, the second terminal device 102, and the third terminal device 103 can be various electronic devices with displays and support web browsing, including but not limited to smartphones, tablets, laptops, and desktop computers.

[0036] Server 105 can be a server that provides various services, such as a backend management server that supports websites browsed by users using the first terminal device 101, the second terminal device 102, and the third terminal device 103 (this is just an example). The backend management server can analyze and process data such as received user requests, and feed back the processing results (such as web pages, information, or data obtained or generated according to user requests) to the terminal devices.

[0037] The gravity gradient meter 106 can be used to collect the gravity gradient of the groundwater level at the point to be measured.

[0038] It should be noted that the determination of groundwater level change based on gravity gradient provided in this application embodiment can generally be performed by the first terminal device 101, the second terminal device 102, or the third terminal device 103. Correspondingly, the device for determining groundwater level change based on gravity gradient provided in this application embodiment can generally be installed in the first terminal device 101, the second terminal device 102, or the third terminal device 103. The method for determining groundwater level change based on gravity gradient provided in this application embodiment can also be performed by a terminal device different from the first terminal device 101, the second terminal device 102, or the third terminal device 103 and capable of communicating with the first terminal device 101, the second terminal device 102, the third terminal device 103, and / or the server 105. Accordingly, the device for determining groundwater level changes based on gravity gradient provided in this application embodiment can also be installed in a terminal device that is different from the first terminal device 101, the second terminal device 102, or the third terminal device 103 and is capable of communicating with the first terminal device 101, the second terminal device 102, the third terminal device 103 and / or the server 105.

[0039] It should be understood that Figure 1 The number of the first terminal device, second terminal device, third terminal device network, server, and gravity gradiometer in the diagram is merely illustrative. Depending on implementation needs, any number of terminal devices, networks, servers, and gravity gradiometers can be included.

[0040] Figure 2 A flowchart illustrating a method for determining groundwater level changes based on gravity gradient according to an embodiment of this application is shown.

[0041] like Figure 2As shown, the method 200 for determining groundwater level change based on gravity gradient in this embodiment includes operation S210 and operation S220.

[0042] In operation S210, the gravity gradient change of each of the multiple groundwater level test points is obtained.

[0043] In operation S220, based on the correlation between gravity gradient and water level, and according to the three-dimensional underground network associated with multiple groundwater level test points and the gravity gradient change of multiple groundwater level test points, the groundwater level change of each of the multiple groundwater level test points is obtained.

[0044] Among them, the change in gravity gradient at the gravity gradient test point and the correlation between gravity gradient and water level are obtained based on the attribute information of multiple grids in the three-dimensional underground network associated with the gravity gradient test point. The three-dimensional underground network associated with the surface measurement point is obtained by dividing the three-dimensional space into grids based on the preset gravity gradient resolution, the preset groundwater level resolution, and the soil saturated moisture content associated with the surface measurement point. The surface measurement point includes the gravity gradient test point and the groundwater level test point.

[0045] Groundwater level measurement points can represent the spatial locations on the Earth's surface where groundwater level data is to be determined.

[0046] In one implementation, a gravity gradiometer can be used to collect the gravity gradients of multiple groundwater level measurement points at different times, thereby obtaining the gravity gradient changes of each of the multiple groundwater level measurement points.

[0047] Before inverting the groundwater level changes at the designated groundwater level monitoring points, it is necessary to filter the gravity gradient changes obtained from each of the multiple monitoring points. The high-frequency, time-varying portion of the gravity gradient changes is identified as the gravity gradient change caused by groundwater variations, while the low-frequency, time-varying portion is identified as the gravity gradient change caused by geological or other sphere variations. The high-frequency, time-varying portion of the gravity gradient changes is then used for subsequent inversion.

[0048] The preset gravity gradient resolution can be set based on the observation accuracy of the gravity gradiometer or experience. The preset groundwater level resolution can be set based on engineering application requirements, hydrogeological parameters of surface measurement points, or experience. For example, hydrogeological parameters of surface measurement points may include aquifer thickness or medium density geological parameters.

[0049] A three-dimensional subsurface network can characterize an adaptive subsurface network corresponding to surface measurement points. The three-dimensional subsurface network includes multiple grids, for example, the grids can be cuboid grids, and the grid size can be non-uniform, for example, the grid size near the surface measurement points is smaller than the grid size far from the surface measurement points.

[0050] Surface measurement points can include gravity gradient measurement points, which can characterize the spatial locations on the surface where gravity gradient data is to be determined. In one implementation, there can be multiple gravity gradient measurement points, and different gravity gradient measurement points correspond to different three-dimensional subsurface networks.

[0051] Based on the preset gravity gradient resolution, preset groundwater level resolution, and soil saturation moisture content associated with the surface measurement points, the three-dimensional space associated with the surface measurement points can be determined. By meshing the three-dimensional space, a three-dimensional underground network associated with the surface measurement points can be obtained.

[0052] The attribute information of the mesh can include the position function of the mesh and the preset water content change of the mesh. The position function can characterize the change of gravity gradient at the point to be measured when the mesh undergoes a unit density change.

[0053] By accumulating the attribute information of multiple grids within a 3D underground network associated with the gravity gradient test point, the gravity gradient change at the test point can be obtained through forward modeling. The water content change at the test point can be determined by summing the preset water content changes of multiple grids. Furthermore, the relationship between the gravity gradient and water level can be established based on both the gravity gradient change and the water content change at the test point.

[0054] By establishing a three-dimensional underground network associated with surface measurement points and combining the gravity gradient changes at the measurement points, the groundwater level changes at multiple groundwater level measurement points can be obtained. Based on this, by adaptively meshing the three-dimensional space determined by the preset gravity gradient resolution, preset groundwater level resolution, and soil saturation moisture content associated with the surface measurement points, a three-dimensional underground network associated with the surface measurement points is obtained. This improves the accuracy of the correlation between gravity gradient and water level, thereby enhancing data processing precision and efficiency. Consequently, the accuracy of the groundwater level changes at multiple groundwater level measurement points is improved, providing technical support for groundwater dynamic detection and inversion.

[0055] Specifically, based on the correlation between gravity gradient and water level, and according to the three-dimensional underground network associated with multiple groundwater level measurement points and the gravity gradient changes of these points, the groundwater level changes of each of the multiple groundwater level measurement points can be obtained. This can include: for any groundwater level measurement point among the multiple groundwater level measurement points, constructing a three-dimensional underground network water content equation based on the gravity gradient changes of the groundwater level measurement point and the three-dimensional underground network associated with it; and solving the three-dimensional underground network water content equation system based on the gravity gradient and water level correlation to obtain the groundwater level changes of each of the multiple groundwater level measurement points.

[0056] There are multiple groundwater level monitoring points. For any one of these monitoring points, in one practical method, based on the gravity gradient change of the groundwater level monitoring points and the associated three-dimensional underground network, a three-dimensional underground network water content equation is constructed using the Eulerian deconvolution method. Then, the water content equations of the three-dimensional underground network constructed from multiple monitoring points are used to form a set of three-dimensional underground network water content equations.

[0057] Specifically, based on the correlation between gravity gradient and water level, the three-dimensional underground network water content equation system, which is composed of multiple three-dimensional underground network water content equations, is solved to obtain the groundwater level change of each of the multiple groundwater level test points. This includes: determining the spatial distribution range of groundwater changes associated with each of the multiple groundwater level test points based on the correlation between gravity gradient and water level; and solving the three-dimensional underground grid water content equation system based on the spatial distribution range of groundwater changes to obtain the groundwater level change of each of the multiple groundwater level test points.

[0058] Based on the correlation between gravity gradient and water level, the spatial distribution range of groundwater changes associated with each of the multiple groundwater level monitoring points can be determined.

[0059] Based on the spatial distribution range of groundwater changes associated with each of the multiple groundwater level monitoring points, the three-dimensional underground network water content equations are cross-solved, and the groundwater level changes of each of the multiple groundwater level monitoring points can be obtained by inversion.

[0060] Specifically, the inversion is a process of constructing a many-to-one mapping relationship between multiple groundwater level changes and a single gravity gradient change based on a three-dimensional underground network associated with the groundwater level measurement point, and obtaining the groundwater level change through fitting.

[0061] By constructing a set of three-dimensional underground network water content equations and performing cross-joint solutions for multiple groundwater level measurement points, the underdetermined or multiple-solution problems caused by many-to-one mapping can be reduced, the inversion constraints can be enhanced, and the accuracy of the groundwater level change inversion results can be improved.

[0062] The above-mentioned method for determining the change in groundwater level based on gravity gradient may further include: for any groundwater level to be measured among multiple groundwater level to be measured points, obtaining the groundwater level of the groundwater level to be measured point based on the change in groundwater level of the groundwater level to be measured point and the groundwater level of known survey points; and correcting the groundwater level of the groundwater level to be measured point based on a preset reference surface to obtain the target groundwater level of the groundwater level to be measured point.

[0063] Groundwater level can characterize the depth of water relative to a predetermined reference surface. The known groundwater level at a survey point can represent the groundwater level at that point obtained using survey wells or other methods. Based on the groundwater level at the survey point, the initial groundwater level at the point to be measured can be obtained using interpolation.

[0064] There can be multiple groundwater level monitoring points. For any one of these multiple groundwater level monitoring points, the groundwater level of the monitoring point can be obtained based on the groundwater level change at the monitoring point and the groundwater level at a known survey point. In other words, the groundwater level of the monitoring point can be obtained by accumulating the groundwater level change at the monitoring point and the initial groundwater level at the monitoring point.

[0065] The target groundwater level can characterize the depth of the water level relative to the Earth's surface. In one implementation, the target groundwater level at the point to be measured can be obtained by subtracting a preset reference surface from the groundwater level at the point to be measured. The preset reference surface can be a horizontal plane with a fixed elevation relative to the surface geoid selected during the forward modeling process.

[0066] By combining the water level change with the groundwater level at known survey points, a groundwater level with a preset reference surface is obtained, which improves the accuracy of the determined groundwater level. Then, based on the reference surface, the water level is corrected for elevation, converting the water level relative to the reference surface into the burial depth relative to the ground surface. This reduces the deviation caused by the reference surface and topographic undulations, making the final target groundwater level more consistent with the actual engineering application scenario and improving the practicality and accuracy of the target groundwater level.

[0067] Figure 3 The flowchart illustrating the inversion of groundwater level at the test point according to an embodiment of this application is shown in the schematic diagram.

[0068] like Figure 3As shown, the groundwater level at the point to be measured is obtained through inversion, including operations S310 to S340.

[0069] In operation S310, the gravity gradient changes of multiple groundwater level test points at different times are obtained and filtered.

[0070] In operation S320, based on the correlation between gravity gradient and water level, and according to the three-dimensional underground network associated with multiple groundwater level test points and the gravity gradient change of multiple groundwater level test points, the groundwater level change of each of the multiple groundwater level test points is obtained.

[0071] In operation S330, the groundwater level at the point to be measured is obtained based on the change in groundwater level at the point to be measured and the known groundwater level at the survey point.

[0072] In operation S340, based on the preset reference surface, the groundwater level at the groundwater level to be measured is corrected to obtain the target groundwater level at the groundwater level to be measured.

[0073] The three-dimensional underground network associated with the surface measurement points can be obtained as follows: Based on the preset gravity gradient resolution, preset groundwater level resolution, and soil saturation moisture content associated with the surface measurement points, determine the vertical measurement range associated with the surface measurement points; based on the vertical measurement range, determine the horizontal measurement range associated with the surface measurement points; and divide the three-dimensional space formed by the vertical and horizontal measurement ranges into a grid to obtain the three-dimensional underground network associated with the surface measurement points.

[0074] Based on the preset gravity gradient resolution, preset groundwater level resolution, and soil saturated moisture content associated with the surface measurement point, the vertical measurement range associated with the surface measurement point is determined, including: determining the vertical distance associated with the surface measurement point based on the preset gravity gradient resolution and soil saturated moisture content associated with the surface measurement point; and determining the vertical measurement range based on the multiple relationship between the vertical distance and the preset groundwater level resolution.

[0075] The soil saturated water content is the aquifer saturated water content corresponding to the surface measurement point, and can be used to characterize the maximum water storage capacity of the soil below the surface at that measurement point. In one implementation, the soil saturated water content can be set to 0.1.

[0076] The preset gravity gradient resolution characterizes the minimum change in gravity gradient that can be detected by the gravity gradiometer. The vertical distance characterizes the upper limit of the vertical depth that a surface measurement point can probe. The preset gravity gradient resolution and the soil saturated moisture content associated with the surface measurement point are positively correlated with the vertical distance; that is, the vertical distance associated with the surface measurement point can be calculated based on the preset gravity gradient resolution and the soil saturated moisture content associated with the surface measurement point. A higher preset gravity gradient resolution results in higher sensitivity of the gravity gradiometer. In one implementation, higher saturated moisture content or higher gravity gradiometer sensitivity corresponds to a larger vertical distance.

[0077] A preset groundwater level resolution characterizes the smallest identifiable change in groundwater level. Based on the ratio between the vertical distance and the preset groundwater level resolution, the vertical measurement range can be determined. For example, with a vertical distance of 40 meters and a preset groundwater level resolution of 0.3 meters, the vertical measurement range can be determined to be 39.9 meters. Similarly, with a vertical distance of 50 meters and a preset groundwater level resolution of 0.1 meters, the vertical measurement range can be determined to be 50 meters.

[0078] Based on the spatial distribution characteristics of the gravity gradient field, with the surface measurement point as the vertex and the vertical measurement range as the axis, there exists a conical surface dividing the positive and negative gravity gradients. The maximum horizontal measurement radius of this conical surface is equal to the horizontal measurement range associated with the surface measurement point. Specifically, the conical surface dividing the positive and negative gravity gradients satisfies the following condition: the maximum horizontal measurement radius of this conical surface is approximately equal to... The vertical measurement range is increased by a factor of 1.

[0079] Specifically, the process of dividing the three-dimensional space formed by the vertical and horizontal measurement ranges into a grid to obtain a three-dimensional underground network associated with surface measurement points can include: dividing the vertical measurement range into equal-interval layers according to a preset groundwater level resolution to obtain multiple horizontal layers; for any horizontal layer, starting from the center point of the horizontal layer, expanding along the horizontal direction according to preset reference grid units to obtain the actual gravity gradient generated by at least one reference grid within the horizontal layer for the surface measurement points, wherein the reference grid is composed of side lengths that are integer multiples of the preset reference grid units and heights that are the distance between the horizontal layer and the surface measurement points; determining at least one grid for the horizontal layer when the actual gravity gradient and the gravity gradient resolution satisfy a preset multiple relationship; and stitching together at least one grid from each of the multiple horizontal layers to obtain the three-dimensional underground network associated with the surface measurement points.

[0080] Based on the preset groundwater level resolution, the vertical measurement range is divided into equal-interval layers to obtain multiple horizontal layers. In other words, by dividing the vertical measurement range into equal-interval layers with the preset groundwater level resolution as the interval, multiple horizontal layers can be obtained.

[0081] In one implementation, the reference grid is formed by radiating four rays from the center point of the horizontal layer in the horizontal directions of east, west, south, and north. The side lengths are integer multiples of the preset reference grid cells, and the height is the distance between the horizontal layer and the surface measurement point. In another implementation, the preset reference grid cell can be a square. Using this square as the base, a cuboid grid is formed with the height equal to the distance between the horizontal layer and the surface measurement point.

[0082] For any one of multiple horizontal layers, starting from the center point of the horizontal layer, the grid is expanded horizontally according to a preset reference grid unit to obtain the actual gravity gradient generated by at least one reference grid at the surface measurement point within the horizontal layer. If a preset multiple relationship is satisfied between the actual gravity gradient and the gravity gradient resolution, at least one grid of the horizontal layer can be determined.

[0083] In one implementation, when the actual gravity gradient is an integer multiple of the gravity gradient resolution, the individual grids for each horizontal layer can be determined. For example, if the actual gravity gradient is equal to the gravity gradient resolution, one grid is determined; if the actual gravity gradient is twice the gravity gradient resolution, another grid is determined; if the actual gravity gradient is three times the gravity gradient resolution, another grid is determined, and so on, until the grid extends to the boundary of the horizontal measurement range, resulting in the complete grid for the horizontal layer. This process is repeated for multiple horizontal layers until the boundary of the vertical measurement range is reached, forming a three-dimensional subsurface network associated with the surface measurement points.

[0084] By pre-setting the gravity gradient resolution, the pre-set groundwater level resolution, and the soil saturated water content, the vertical and horizontal measurement ranges associated with surface measurement points are determined respectively. This ensures that the boundaries of the three-dimensional space match the required observation accuracy and hydrogeological characteristics, reducing redundant calculations and improving data processing efficiency. By first determining the upper limit of depth and then combining it with the multiple relationship of the pre-set groundwater level resolution, the vertical measurement range is determined so that it satisfies both the detectability condition of the gravity gradient and the accuracy of water level detection, improving the rationality of grid division. The vertical range is layered at equal intervals according to the pre-set groundwater level resolution. The reference grid is expanded outward from the center of the horizontal layer, and the relationship between the actual gravity gradient and the gravity gradient resolution is used as the criterion for grid validity. This achieves adaptive density of the grid in the horizontal direction, and finally, the grids are stitched together to form a three-dimensional underground network. This makes the grid distribution consistent with the spatial distribution of the gravity gradient. That is, in areas where the gravity gradient is less affected by groundwater changes, the calculation is simplified, and a rapid response to higher frequency environments can be achieved. In areas where the gravity gradient is sensitive to groundwater changes, or at the interface between positive and negative changes in the response, the calculation is focused to reduce errors, thereby improving the accuracy, stability and computational efficiency of the joint data processing of groundwater level and gravity gradient changes.

[0085] Figure 4 A flowchart illustrating the construction of a three-dimensional underground network associated with surface measurement points according to an embodiment of this application is shown.

[0086] like Figure 4 As shown, constructing a three-dimensional underground network associated with surface measurement points includes operations 410 to S460.

[0087] In operation S410, determine the vertical measurement range and the horizontal measurement range.

[0088] When operating the S420, the vertical measurement range is divided into equal-interval layers according to the preset groundwater level resolution, resulting in multiple horizontal layers. The center point of each horizontal layer is moved downwards with the horizontal layer as the step size.

[0089] In operation S430, starting from the center point of the horizontal layer, the system expands along the horizontal direction according to the preset reference grid cells to obtain the actual gravity gradient generated by at least one reference grid cell in the horizontal layer for the surface measurement point.

[0090] In operation S440, under the condition that a preset multiple relationship is satisfied between the actual gravity gradient and the gravity gradient resolution, at least one grid of the horizontal layer is determined until the boundary of the horizontal measurement range is reached.

[0091] In operation S450, check whether the center point of the horizontal layer reaches the boundary of the vertical measurement range. If yes, execute operation S460; otherwise, execute operation S420.

[0092] Using the S460, a three-dimensional underground network associated with surface measurement points is constructed.

[0093] Figure 5 A schematic diagram of a three-dimensional underground network associated with surface measurement points according to an embodiment of this application is shown.

[0094] like Figure 5 As shown, the horizontal axis represents the distance the center point of the horizontal layer moves northward, ranging from approximately -2 to +2. The vertical axis represents the gravity gradient generated by the grid at the center point. It can be seen that the three-dimensional subsurface network associated with surface measurement points can include multiple grids, and the grid size is not fixed.

[0095] The change in gravity gradient at the point to be measured and the relationship between gravity gradient and water level can be obtained as follows: For any grid among the multiple grids in the three-dimensional underground network associated with the point to be measured, the gravity gradient change component of the grid is determined according to the product between the grid's position function and the preset water content change; the gravity gradient change components of each of the multiple grids are accumulated to obtain the change in gravity gradient at the point to be measured; the relationship between gravity gradient and water level is determined based on the change in gravity gradient at the point to be measured and the preset water content change of each of the multiple grids in the three-dimensional underground network associated with the point to be measured.

[0096] Forward modeling can be a process of integrating the changes in gravity gradient generated by irregular groundwater bodies and accumulating the changes in gravity gradient generated within several cuboid grids to calculate the gravity gradient changes associated with the points to be measured that generate gravity gradients.

[0097] Before performing forward modeling on the gravity gradient changes at the points to be measured, a horizontal plane with a fixed elevation relative to the geoid needs to be selected as the preset reference surface for gravity gradient forward modeling. During the forward modeling calculation of gravity gradient changes, all gravity gradient points should be assumed to be on this preset reference surface to correct for the impact of topographic relief changes.

[0098] The three-dimensional subsurface network associated with the gravity gradient measurement point can include multiple grids. The product of the position function of each grid and the preset water content change can determine the gravity gradient change component of that grid. Then, by accumulating the gravity gradient change components of the multiple grids in the three-dimensional subsurface network associated with the gravity gradient measurement point, the gravity gradient change of the measurement point can be obtained through forward modeling.

[0099] In one implementation, if the initial water content change of each grid in the 3D underground network associated with the gravity gradient test point is incomplete, it is necessary to supplement the initial water content change to obtain the preset water content change of each grid in the 3D underground network associated with the gravity gradient test point. In another implementation, if the initial water content change of some grids in the 3D underground network is missing, a quadratic interpolation method can be used to interpolate the 3D underground network to obtain the preset water content change of each grid in the 3D underground network associated with the gravity gradient test point. In yet another implementation, if the initial water content change of many grids in the 3D underground network is missing, the existing initial water content change can be symmetrically reversed, and the boundaries smoothed to extrapolate and extend the 3D underground network to obtain the preset water content change of each grid in the 3D underground network associated with the gravity gradient test point.

[0100] The preset water content changes of multiple grids in the three-dimensional underground network associated with the gravity gradient test point can be accumulated to determine the water content change associated with the gravity gradient test point.

[0101] The relationship between gravity gradient and water level can be established based on the change in gravity gradient at the point to be measured and the change in water content associated with the point to be measured. This relationship can reflect the mapping law between water level change and gravity gradient response.

[0102] The gravity gradient change components of each grid are obtained by multiplying the position function with the preset water content change by the grid-by-grid calculation, and then summed to obtain the total gravity gradient change of the test point. Based on the total gravity gradient change and the preset water content change of each grid, the correlation between gravity gradient and water level is constructed, realizing the quantitative mapping between groundwater change and gravity gradient change, improving the reliability of the correlation and the accuracy of subsequent water level inversion.

[0103] Figure 6 The flowchart illustrating the forward modeling of the gravity gradient change at the test point according to an embodiment of this application is shown.

[0104] like Figure 6 As shown, the forward modeling yields the gravity gradient change at the point to be measured, including operations S610 to S640.

[0105] In operation S610, the initial water content change of each of the multiple grids in the three-dimensional underground network associated with the point to be measured by gravity gradient is obtained.

[0106] In operation S620, based on a preset reference surface, the initial water content change of each of the multiple grids in the three-dimensional underground network associated with the gravity gradient test point is corrected to obtain the intermediate water content change of each of the multiple grids in the three-dimensional underground network associated with the gravity gradient test point.

[0107] In operation S630, at least one of the following operations is performed on the intermediate water content change of each of the multiple grids in the three-dimensional underground network associated with the gravity gradient test point: interpolation, extrapolation, or boundary smoothing, to obtain the preset water content change of each of the multiple grids in the three-dimensional underground network associated with the gravity gradient test point.

[0108] In operation S640, the gravity gradient change component of the grid is determined based on the product between the grid position function and the preset moisture content change. The gravity gradient change components of multiple grids are then accumulated to obtain the gravity gradient change of the point to be measured.

[0109] Based on the above-mentioned method for determining groundwater level changes based on gravity gradient, this application also provides a device for determining groundwater level changes based on gravity gradient. The following will be combined with... Figure 7 The device is described in detail.

[0110] Figure 7 A schematic diagram of a device for determining groundwater level changes based on gravity gradient according to an embodiment of this application is shown.

[0111] like Figure 7 As shown, the groundwater level change determination device 700 based on gravity gradient in this embodiment includes an acquisition module 710 and a first acquisition module 720.

[0112] The acquisition module 710 is used to acquire the gravity gradient change of each of the multiple groundwater level measurement points. In one embodiment, the acquisition module 710 can be used to perform the operation S210 described above, which will not be repeated here.

[0113] The first acquisition module 720 is used to obtain the groundwater level change of each of the multiple groundwater level test points based on the correlation between gravity gradient and water level, according to the three-dimensional underground network associated with multiple groundwater level test points and the gravity gradient change of the multiple groundwater level test points. The gravity gradient change of the gravity gradient test points and the correlation between gravity gradient and water level are obtained based on the attribute information of the multiple grids included in the three-dimensional underground network associated with the gravity gradient test points. The three-dimensional underground network associated with the surface measurement points is obtained by dividing the three-dimensional space into grids according to the preset gravity gradient resolution, the preset groundwater level resolution and the soil saturated water content associated with the surface measurement points. The surface measurement points include gravity gradient test points and groundwater level test points. The attribute information of the grids includes the grid position function and the preset water content change. The position function represents the gravity gradient change caused at the gravity gradient test points when the grid undergoes a unit density change. In one embodiment, the first obtaining module 720 can be used to perform the operation S220 described above, which will not be repeated here.

[0114] According to an embodiment of this application, the first obtaining module 720 includes a construction unit and a solving unit.

[0115] The construction unit is used to construct a three-dimensional underground network water content equation for any groundwater level measurement point among multiple groundwater level measurement points, based on the change in gravity gradient of the groundwater level measurement point and the three-dimensional underground network associated with the groundwater level measurement point.

[0116] The solution unit is used to solve the three-dimensional underground network water content equation system based on the relationship between gravity gradient and water level, and to obtain the groundwater level change of each groundwater level test point.

[0117] According to an embodiment of this application, the solving unit includes a first determining subunit and a solving subunit.

[0118] The first determining sub-unit is used to determine the spatial distribution range of groundwater changes associated with each of the multiple groundwater level measurement points based on the correlation between gravity gradient and water level.

[0119] The solution sub-element is used to solve the three-dimensional underground grid water content equation system based on the spatial distribution range of groundwater changes, and obtain the groundwater level change of each groundwater level measurement point.

[0120] According to an embodiment of this application, the device 700 for determining groundwater level changes based on gravity gradient further includes a second acquisition module and a correction module.

[0121] The second acquisition module is used to obtain the groundwater level of any groundwater level test point among multiple groundwater level test points, based on the groundwater level change of the groundwater level test point and the groundwater level of known survey points. The groundwater level represents the depth of the water level relative to a preset reference surface.

[0122] The correction module is used to correct the groundwater level at the groundwater level test point based on a preset reference surface, so as to obtain the target groundwater level at the groundwater level test point, wherein the target groundwater level represents the depth of the water level relative to the ground surface.

[0123] According to an embodiment of this application, the first obtaining module 720 includes a first determining unit, a second determining unit, and a dividing unit. The first determining unit is used to determine the vertical measurement range associated with the surface measurement point based on a preset gravity gradient resolution, a preset groundwater level resolution, and the soil saturated moisture content associated with the surface measurement point.

[0124] The second determining unit is used to determine the horizontal measurement range associated with the ground measurement point based on the vertical measurement range.

[0125] The dividing unit is used to divide the three-dimensional space into grids based on the vertical and horizontal measurement ranges, thereby obtaining a three-dimensional underground network associated with the surface measurement points.

[0126] According to an embodiment of this application, the first determining unit includes a second determining subunit and a third determining subunit.

[0127] The second determining subunit is used to determine the vertical distance associated with the surface measurement point based on the preset gravity gradient resolution and the soil saturated water content associated with the surface measurement point, wherein the vertical distance characterizes the upper limit of the vertical downward probe depth of the surface measurement point.

[0128] The third determining sub-unit is used to determine the vertical measurement range based on the multiple relationship between the vertical distance and the preset groundwater level resolution.

[0129] According to embodiments of this application, the partitioning unit includes a hierarchical subunit, an extended subunit, a fourth determining subunit, and a splicing subunit.

[0130] The layered sub-unit is used to divide the vertical measurement range into equal-interval layers according to the preset groundwater level resolution, resulting in multiple horizontal layers.

[0131] An extended sub-unit is used to extend any one of multiple horizontal layers along the horizontal direction according to a preset reference grid unit, starting from the center point of the horizontal layer, to obtain the actual gravity gradient generated by at least one reference grid unit in the horizontal layer for the surface measurement point. The reference grid unit is composed of a side length that is an integer multiple of the preset reference grid unit unit and a height that is the distance between the horizontal layer and the surface measurement point.

[0132] The fourth determination sub-unit is used to determine at least one grid of the horizontal layer when a preset multiple relationship is satisfied between the actual gravity gradient and the gravity gradient resolution.

[0133] The splicing sub-unit is used to splice at least one grid from multiple horizontal layers to obtain a three-dimensional underground network associated with surface measurement points.

[0134] According to an embodiment of this application, the first obtaining module 720 further includes a third determining unit, a processing unit, and a fourth determining unit.

[0135] The third determining unit is used to determine the gravity gradient change component of any grid in the three-dimensional underground network associated with the point to be measured by gravity gradient, based on the product between the grid's position function and a preset water content change.

[0136] The processing unit is used to accumulate the gravity gradient change components of multiple grids to obtain the gravity gradient change of the point to be measured.

[0137] The fourth determining unit is used to determine the relationship between gravity gradient and water level based on the gravity gradient change of the gravity gradient test point and the preset water content change of each of the multiple grids in the three-dimensional underground network associated with the gravity gradient test point.

[0138] According to embodiments of this application, any plurality of modules in the acquisition module 710 and the first acquisition module 720 can be combined into one module, or any one of these modules can be split into multiple modules. Alternatively, at least a portion of the functionality of one or more of these modules can be combined with at least a portion of the functionality of other modules and implemented in one module. According to embodiments of this application, at least one of the acquisition module 710 and the first acquisition module 720 can be at least partially implemented as hardware circuitry, such as a field-programmable gate array (FPGA), a programmable logic array (PLA), a system-on-a-chip, a system-on-a-substrate, a system-on-package, an application-specific integrated circuit (ASIC), or any other reasonable means of integrating or packaging circuitry, or implemented in software, hardware, or firmware, or in any appropriate combination of any of these three implementation methods. Alternatively, at least one of the acquisition module 710 and the first acquisition module 720 can be at least partially implemented as a computer program module, which, when run, can perform corresponding functions.

[0139] Figure 8 A block diagram of an electronic device suitable for implementing a method for determining groundwater level changes based on gravity gradients, according to an embodiment of this application, is illustrated schematically.

[0140] like Figure 8 As shown, an electronic device 800 according to an embodiment of this application includes a processor 801, which can perform various appropriate actions and processes according to a program stored in ROM 802 or a program loaded from storage portion 808 into RAM 803. The processor 801 may include, for example, a general-purpose microprocessor (e.g., a CPU), an instruction set processor and / or an associated chipset and / or a special-purpose microprocessor (e.g., an application-specific integrated circuit (ASIC)), etc. The processor 801 may also include onboard memory for caching purposes. The processor 801 may include a single processing unit or multiple processing units for performing different actions of the method flow according to an embodiment of this application.

[0141] RAM 803 stores various programs and data required for the operation of electronic device 800. Processor 801, ROM 802, and RAM 803 are interconnected via bus 804. Processor 801 executes various operations of the method flow according to embodiments of this application by executing programs in ROM 802 and / or RAM 803. It should be noted that the programs may also be stored in one or more memories other than ROM 802 and RAM 803. Processor 801 may also execute various operations of the method flow according to embodiments of this application by executing programs stored in said one or more memories.

[0142] According to embodiments of this application, the electronic device 800 may further include an input / output (I / O) interface 805, which is also connected to a bus 804. The electronic device 800 may also include one or more of the following components connected to the input / output (I / O) interface 805: an input section 806 including a keyboard, mouse, etc.; an output section 807 including a cathode ray tube (CRT), liquid crystal display (LCD), etc., and a speaker, etc.; a storage section 808 including a hard disk, etc.; and a communication section 809 including a network interface card such as a LAN card, modem, etc. The communication section 809 performs communication processing via a network such as the Internet. A drive 810 is also connected to the input / output (I / O) interface 805 as needed. A removable medium 811, such as a disk, optical disk, magneto-optical disk, semiconductor memory, etc., is installed on the drive 810 as needed so that computer programs read from it can be installed into the storage section 808 as needed.

[0143] This application also provides a computer-readable storage medium, which may be included in the device / apparatus / system described in the above embodiments; or it may exist independently and not assembled into the device / apparatus / system. The computer-readable storage medium carries one or more programs, which, when executed, implement the method according to the embodiments of this application.

[0144] According to embodiments of this application, the computer-readable storage medium can be a non-volatile computer-readable storage medium, such as including but not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In this application, the computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device. For example, according to embodiments of this application, the computer-readable storage medium may include ROM 802 and / or RAM 803 and / or one or more memories other than ROM 802 and RAM 803 described above.

[0145] Embodiments of this application also include a computer program product comprising a computer program containing program code for performing the methods shown in the flowchart. When the computer program product is run on a computer system, the program code enables the computer system to implement the method for determining groundwater level changes based on gravity gradients provided in the embodiments of this application.

[0146] When the computer program is executed by the processor 801, it performs the functions defined in the system / apparatus of this application embodiment. According to the embodiments of this application, the systems, apparatuses, modules, units, etc., described above can be implemented by computer program modules.

[0147] In one embodiment, the computer program may rely on a tangible storage medium such as an optical storage device or a magnetic storage device. In another embodiment, the computer program may also be transmitted and distributed in the form of signals over a network medium, and may be downloaded and installed via the communication section 809, and / or installed from a removable medium 811. The program code contained in the computer program can be transmitted using any suitable network medium, including but not limited to: wireless, wired, etc., or any suitable combination thereof.

[0148] In such an embodiment, the computer program can be downloaded and installed from a network via the communication section 809, and / or installed from the removable medium 811. When the computer program is executed by the processor 801, it performs the functions defined in the system of this application embodiment. According to the embodiments of this application, the systems, devices, apparatuses, modules, units, etc., described above can be implemented by computer program modules.

[0149] According to embodiments of this application, program code for executing the computer programs provided in the embodiments of this application can be written in any combination of one or more programming languages. Specifically, these computational programs can be implemented using high-level procedural and / or object-oriented programming languages, and / or assembly / machine languages. Programming languages ​​include, but are not limited to, languages ​​such as Java, C++, Python, "C", or similar programming languages. The program code can be executed entirely on the user's computing device, partially on the user's device, partially on a remote computing device, or entirely on a remote computing device or server. In cases involving remote computing devices, the remote computing device can be connected to the user's computing device via any type of network, including a local area network (LAN) or a wide area network (WAN), or it can be connected to an external computing device (e.g., via the Internet using an Internet service provider).

[0150] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this application. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in a block diagram or flowchart, and combinations of blocks in a block diagram or flowchart, may be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.

[0151] Those skilled in the art will understand that the features described in the various embodiments of this application can be combined and / or combined in various ways, even if such combinations or combinations are not explicitly described in this application. In particular, the features described in the various embodiments of this application can be combined and / or combined in various ways without departing from the spirit and teachings of this application. All such combinations and / or combinations fall within the scope of this application.

[0152] The embodiments of this application have been described above. However, these embodiments are merely illustrative and not intended to limit the scope of this application. Although various embodiments have been described above, this does not mean that the measures in the various embodiments cannot be used advantageously in combination. Without departing from the scope of this application, those skilled in the art can make various substitutions and modifications, all of which should fall within the scope of this application.

Claims

1. A method for determining groundwater level changes based on gravity gradient, comprising: Obtain the gravity gradient change at multiple groundwater level measurement points; Based on the relationship between gravity gradient and water level, the groundwater level change of each of the multiple groundwater level test points is obtained according to the three-dimensional underground network associated with the multiple groundwater level test points and the gravity gradient change of the multiple groundwater level test points. The gravity gradient change at the gravity gradient test point and the correlation between gravity gradient and water level are obtained based on the attribute information of multiple grids in the three-dimensional underground network associated with the gravity gradient test point. The three-dimensional underground network associated with the surface measurement point is obtained by dividing a three-dimensional space into grids based on a preset gravity gradient resolution, a preset groundwater level resolution, and the soil saturated moisture content associated with the surface measurement point. The surface measurement point includes the gravity gradient test point and the groundwater level test point. The attribute information of the grid includes the grid's position function and a preset moisture content change. The position function characterizes the gravity gradient change caused at the gravity gradient test point when the grid undergoes a unit density change.

2. The method according to claim 1, wherein, The method, based on the correlation between gravity gradient and water level, obtains the groundwater level change of each of the multiple groundwater level measurement points according to the three-dimensional underground network associated with the multiple groundwater level measurement points and the gravity gradient change of the multiple groundwater level measurement points, including: For any groundwater level measurement point among multiple groundwater level measurement points, a three-dimensional underground network water content equation is constructed based on the gravity gradient change of the groundwater level measurement point and the three-dimensional underground network associated with the groundwater level measurement point. Based on the relationship between gravity gradient and water level, the three-dimensional underground network water content equation system, which is composed of multiple three-dimensional underground network water content equations, is solved to obtain the groundwater level change of each of the multiple groundwater level test points.

3. The method according to claim 2, wherein, Based on the correlation between the gravity gradient and the water level, the three-dimensional underground network water content equation system, composed of multiple three-dimensional underground network water content equations, is solved to obtain the groundwater level change at each of the multiple groundwater level measurement points, including: Based on the correlation between the gravity gradient and the water level, the spatial distribution range of groundwater changes associated with each of the multiple groundwater level measurement points is determined. Based on the spatial distribution range of groundwater changes, the water content equations of the three-dimensional underground grid are solved to obtain the groundwater level changes of each of the groundwater level measurement points.

4. The method according to any one of claims 1 to 3, further comprising: For any groundwater level to be measured among a plurality of groundwater level to be measured points, the groundwater level of the groundwater level to be measured point is obtained based on the groundwater level change of the groundwater level to be measured point and the groundwater level of the known survey points, wherein the groundwater level represents the depth of the water level relative to a preset reference surface. Based on a preset reference surface, the groundwater level at the groundwater level to be measured is corrected to obtain the target groundwater level at the groundwater level to be measured, wherein the target groundwater level represents the depth of the water level relative to the ground surface.

5. The method according to any one of claims 1 to 4, wherein, The three-dimensional subsurface network associated with the surface measurement points is obtained as follows: Based on the preset gravity gradient resolution, the preset groundwater level resolution, and the soil saturated moisture content associated with the surface measurement point, the vertical measurement range associated with the surface measurement point is determined. Based on the vertical measurement range, determine the horizontal measurement range associated with the surface measurement point; The three-dimensional space formed by the vertical and horizontal measurement ranges is divided into grids to obtain a three-dimensional underground network associated with the surface measurement points.

6. The method according to claim 5, wherein, The step of determining the vertical measurement range associated with the surface measurement point based on the preset gravity gradient resolution, the preset groundwater level resolution, and the soil saturated moisture content associated with the surface measurement point includes: Based on the preset gravity gradient resolution and the soil saturated water content associated with the surface measurement point, the vertical distance associated with the surface measurement point is determined, wherein the vertical distance characterizes the upper limit of the vertical downward probe depth of the surface measurement point. The vertical measurement range is determined based on the multiple relationship between the vertical distance and the preset groundwater level resolution.

7. The method according to claim 5 or 6, wherein, The step of dividing the three-dimensional space formed by the vertical and horizontal measurement ranges into a grid to obtain a three-dimensional underground network associated with the surface measurement points includes: Based on the preset groundwater level resolution, the vertical measurement range is divided into equal-interval layers to obtain multiple horizontal layers. For any of the plurality of horizontal layers, starting from the center point of the horizontal layer, the system expands along the horizontal direction according to a preset reference grid unit to obtain the actual gravity gradient generated by at least one reference grid in the horizontal layer on the surface measurement point. The reference grid is constructed with an integer multiple of the preset reference grid unit as the side length and the distance between the horizontal layer and the surface measurement point as the height. If a preset multiple relationship is satisfied between the actual gravity gradient and the gravity gradient resolution, at least one grid of the horizontal layer is determined; At least one grid from each of the multiple horizontal layers is stitched together to obtain a three-dimensional underground network associated with the surface measurement points.

8. The method according to any one of claims 5 to 7, wherein, The change in gravity gradient at the measured point and the correlation between gravity gradient and water level are obtained as follows: For any grid among multiple grids in the three-dimensional underground network associated with the gravity gradient test point, the gravity gradient change component of the grid is determined based on the product between the grid's position function and a preset water content change. The gravity gradient change components of each of the multiple grids are accumulated to obtain the gravity gradient change at the gravity gradient test point. The relationship between gravity gradient and water level is determined based on the change in gravity gradient at the gravity gradient test point and the preset water content change of each of the multiple grids in the three-dimensional underground network associated with the gravity gradient test point.

9. An electronic device, comprising: One or more processors; Memory, used to store one or more computer programs. The characteristic feature is that the one or more processors execute the one or more computer programs to implement the steps of the method according to any one of claims 1 to 8.

10. A computer-readable storage medium having stored thereon executable instructions that, when executed by a processor, cause the processor to perform the method of any one of claims 1 to 8.

Images