A groundwater storage inversion method and system considering deformation of ground soil layer

Abstract

The present application relates to the technical field of electronic digital data processing, in particular to a groundwater reserve inversion method and system considering ground soil layer deformation. The method comprises: determining a plurality of groundwater areas and total groundwater content of each groundwater area according to groundwater content of each point in the underground; determining instability degree of each adjacent soil area according to gas content of each point in each adjacent soil area of the groundwater area for each groundwater area; determining environmental stability index of the groundwater area according to the total groundwater content of the groundwater area and the instability degree of each adjacent soil area; and inverting future groundwater content of the groundwater area according to change data of the groundwater content of the groundwater area and the environmental stability index. The method improves the accuracy of groundwater reserve inversion.

Description

Technical Field

[0001] This invention relates to the field of electronic digital data processing technology, and specifically to a method and system for groundwater storage inversion that takes into account ground soil deformation. Background Technology

[0002] Water resources are vital natural resources and an indispensable foundation for human production and daily life. Groundwater reserves are a crucial component of the water cycle. In recent years, the enormous depletion of global groundwater resources has posed a significant threat to global water security and may lead to problems such as declining agricultural productivity and land subsidence. Therefore, it is necessary to conduct inversion analysis of groundwater reserves and develop rational planning for their use.

[0003] Currently, the inversion of groundwater reserves typically relies on GRACE (Gravity Recovery and Climate Experiment) satellite data and hydro-meteorological data to observe changes in groundwater reserves and predict future changes. However, human activities involving groundwater extraction significantly impact changes in groundwater reserves, making the inversion data less accurate. Summary of the Invention

[0004] To address the technical problem of inaccurate groundwater storage inversion, the present invention aims to provide a groundwater storage inversion method and system that considers surface soil deformation. The specific technical solution adopted is as follows:

[0005] This invention provides a method for groundwater storage inversion considering surface soil deformation, the method comprising:

[0006] Based on the groundwater content at various underground points, determine several underground water bodies and the total groundwater content of each of the underground water bodies;

[0007] For each of the underground water bodies, the instability of each adjacent soil region is determined based on the gas content at each point in each adjacent soil region of the underground water body.

[0008] The environmental stability index of the underground water body is determined based on the total groundwater content of the underground water body and the instability of each of the adjacent soil regions.

[0009] Based on the data on changes in groundwater content in the underground water area and the environmental stability index, the future groundwater content of the underground water area can be inverted.

[0010] According to the groundwater storage inversion method considering surface soil deformation provided by the present invention, the step of determining a number of underground water bodies and the total groundwater content of each underground water body based on the groundwater content at various underground points includes:

[0011] Based on the groundwater content at various underground points, several underground water bodies were identified;

[0012] For each of the underground water bodies, the location coordinates of a representative point of the underground water body are determined based on the groundwater content at each point in the underground water body.

[0013] The total groundwater content of the underground water body is determined based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, as well as the groundwater content at each point.

[0014] According to the groundwater storage inversion method considering surface soil deformation provided by the present invention, the determination of several underground water bodies based on the groundwater content at various underground points includes:

[0015] Based on the groundwater content at each underground point, a set of target points is selected from these points;

[0016] Each time, a point is randomly selected from the unmarked points in the target point set as a seed point. Starting from the seed point, regional growth is performed based on the groundwater content of each point in the neighborhood of the seed point to obtain the underground water area corresponding to the seed point, and each point in the underground water area is marked.

[0017] This continues until all points in the target point set are marked to underground water bodies, resulting in several underground water bodies.

[0018] According to the groundwater storage inversion method considering surface soil deformation provided by the present invention, the step of selecting a target point set from each point based on the groundwater content at each point includes:

[0019] The average groundwater content at each underground point is calculated and used as the groundwater content threshold.

[0020] Points with groundwater content greater than or equal to the groundwater content threshold are identified as target points, forming a set of target points.

[0021] According to the groundwater storage inversion method considering surface soil deformation provided by the present invention, determining the location coordinates of representative points of the groundwater area based on the groundwater content at each point in the groundwater area includes:

[0022] The extreme points are determined based on the groundwater content at each point in the underground water area;

[0023] Based on the groundwater content at each extreme point, a first weight is determined for each extreme point.

[0024] Based on the first weight of each extreme point, the position coordinates of each extreme point are weighted and summed to obtain the position coordinates of the representative point of the underground water area.

[0025] According to the groundwater storage inversion method considering surface soil deformation provided by the present invention, determining the total groundwater content of the groundwater area based on the distance between the location coordinates of each point in the groundwater area and the location coordinates of the representative point, and the groundwater content of each point, includes:

[0026] The second weight of each point in the underground water area is determined based on the distance between the position coordinates of each point in the underground water area and the position coordinates of the representative point.

[0027] Based on the second weight of each point in the underground water body, the groundwater content of each point in the underground water body is weighted and summed to obtain the total groundwater content of the underground water body.

[0028] According to the groundwater storage inversion method considering surface soil deformation provided by the present invention, determining the instability degree of each adjacent soil region based on the gas content at each point in each adjacent soil region of the groundwater includes:

[0029] For each adjacent soil region of the underground water body, the average gas content at each point in the adjacent soil region is determined;

[0030] Determine the number of gas phase target points in the adjacent soil regions whose gas content is greater than or equal to the average gas content;

[0031] The instability of the adjacent soil regions is determined based on the average gas content, the number of gas phase target points, and the distance between each gas phase target point.

[0032] According to the groundwater storage inversion method considering surface soil deformation provided by the present invention, the step of determining the environmental stability index of the groundwater area based on the total groundwater content of the groundwater area and the instability of each of the adjacent soil regions includes:

[0033] The degree of influence of each adjacent soil region on the environmental stability of the groundwater is determined based on the distance between the centroid of each adjacent soil region and the representative point of the groundwater, as well as the number of points in each adjacent soil region.

[0034] Based on the degree of influence of each of the adjacent soil regions on the environmental stability of the groundwater, the instability of each of the adjacent soil regions is weighted and summed to obtain the degree of soil instability surrounding the groundwater.

[0035] The environmental stability index of the underground water body is determined based on the total groundwater content and the instability of the surrounding soil.

[0036] According to the groundwater storage inversion method considering surface soil deformation provided by the present invention, the step of inverting the future groundwater content of the underground water body based on the groundwater content change data of the underground water body and the environmental stability index includes:

[0037] Based on the data on the change in groundwater content in the underground water area, determine the initial rate of change in groundwater content in the underground water area at the current moment;

[0038] The initial groundwater content change rate is weighted according to the environmental stability index of the underground water body to obtain the groundwater content change rate of the underground water body at the current moment.

[0039] Based on the rate of change of groundwater content in the underground water body at the current moment, the future groundwater content of the underground water body can be inverted.

[0040] This invention provides a groundwater storage inversion system that considers surface soil deformation. The system includes a memory and a processor. The memory is used to store executable program code. The processor is used to call and run the executable program code from the memory to implement the groundwater storage inversion method considering surface soil deformation provided by this invention.

[0041] This invention has the following beneficial effects: Based on the groundwater content at various underground points, several underground water bodies and the total groundwater content of each underground water body are determined. For each underground water body, the instability of each adjacent soil region is determined based on the gas content at each point in each adjacent soil region. Based on the total groundwater content of the underground water body and the instability of each adjacent soil region, the environmental stability index of the underground water body is determined. This takes into account the soil conditions around the groundwater and their impact on the underground water body. Based on the changes in groundwater content and the environmental stability index of the underground water body, the future groundwater content of the underground water body is inverted. This can offset the impact of human mining activities on the accuracy of the inversion results and improve the accuracy of groundwater storage inversion. Attached Figure Description

[0042] To more clearly illustrate the technical solutions and advantages in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0043] Figure 1 This is a flowchart illustrating a groundwater storage inversion method considering surface soil deformation, provided in one embodiment of the present invention.

[0044] Figure 2 A schematic diagram of an adjacent soil region of an underground water body provided in one embodiment of the present invention;

[0045] Figure 3 This is a schematic diagram of a process for determining the total groundwater content of an underground water body according to an embodiment of the present invention;

[0046] Figure 4 This is a schematic diagram of a process for determining the instability of adjacent soil regions according to an embodiment of the present invention;

[0047] Figure 5 This is a schematic diagram of a groundwater storage inversion system that takes into account the deformation of the ground soil layer, provided as an embodiment of the present invention. Detailed Implementation

[0048] To further illustrate the technical means and effects adopted by the present invention to achieve its intended purpose, the following, in conjunction with the accompanying drawings and preferred embodiments, details the specific implementation, structure, features, and effects of a groundwater storage inversion method and system considering surface soil deformation proposed according to the present invention. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. Furthermore, specific features, structures, or characteristics in one or more embodiments can be combined in any suitable form.

[0049] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains.

[0050] The following description, in conjunction with the accompanying drawings, details the specific scheme of the groundwater storage inversion method and system that takes into account surface soil deformation provided by the present invention.

[0051] Please see Figure 1 The diagram illustrates a flowchart of a groundwater storage inversion method considering surface soil deformation according to an embodiment of the present invention, comprising the following steps:

[0052] Step 101: Determine the groundwater content of several underground water bodies and the total groundwater content of all underground water bodies based on the groundwater content at various underground points.

[0053] The groundwater content at each underground point refers to the groundwater content within the local area where that point is located; specifically, it represents the percentage of groundwater volume within that local area relative to the total volume of that area. The values ​​for the groundwater content at each underground point are determined by... The total groundwater content of an underground water body represents the overall groundwater content of that body, i.e., the proportion of groundwater volume to the total volume of the underground water body. A higher total groundwater content value indicates a higher overall water content in the underground water body. The total groundwater content of an underground water body is measured within... Inside.

[0054] In one embodiment, the area requiring groundwater storage inversion is first determined, and then ground-penetrating radar (GPR) technology is used to survey the underground. GPR technology transmits high-frequency, polarized radio waves underground. When these electromagnetic waves strike various objects buried underground, they generate reflected waves. The antenna receives these reflected waves, and the content of various substances in the local area where each point is located is determined based on the reflected waves.

[0055] In one embodiment, an arbitrary point at a horizontal level in the surveyed area is taken as the origin O of the coordinate system. The horizontal plane is designated as the xOy plane, and the vertically downward direction is defined as the positive z-axis. A three-dimensional spatial coordinate system is established, discretizing the underground space. Three-dimensional point cloud data is obtained through sampling using ground-penetrating radar technology. The three-dimensional point cloud data includes the location of each point within the underground three-dimensional space and the content of various substances within the local area where each point is located.

[0056] In one embodiment, the content of each substance within the local area where the point is located includes the content of groundwater (liquid phase), non-groundwater substances (solid phase), and gas (gas phase) within the local area where the point is located. The groundwater content, non-groundwater substance content, and gas content respectively represent the proportion of the volume of groundwater, the volume of non-groundwater substances, and the volume of gas within the local area where the point is located to the total volume of the local area. The value ranges of groundwater content, non-groundwater substance content, and gas content are all within... Within this range, the sum of groundwater content, non-groundwater substance content, and gas content equals 1. Let the first... The groundwater content, non-groundwater substance content, and gas content within the local area of ​​each three-dimensional spatial point are respectively... , , ,but .

[0057] In one embodiment, several underground water bodies can be determined based on the groundwater content at various underground points.

[0058] In one embodiment, a set of target points can be selected from each underground point based on the groundwater content of each point, and the points in the target point set can be used as seed points for regional growth to obtain several underground water bodies.

[0059] In one embodiment, points with groundwater content greater than or equal to a groundwater content threshold can be identified as target points, forming a set of target points. In one embodiment, the groundwater content threshold can be preset, determined based on the average groundwater content at various underground points, or determined in other ways, without limitation.

[0060] In one embodiment, for each underground water body, the location coordinates of a representative point of the underground water body are determined based on the groundwater content at each point in the underground water body, and the total groundwater content of the underground water body is determined based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, as well as the groundwater content at each point.

[0061] Step 102: For each underground water body, determine the instability of each adjacent soil region based on the gas content at each point in each adjacent soil region.

[0062] In one embodiment, the soil region surrounding the edge of the underground water body is taken as the initial surrounding soil region. A seed point is randomly selected from the initial surrounding soil region. Starting from the seed point, region growth is performed based on the similarity between the seed point and the surrounding points in terms of groundwater content, non-groundwater substance content, and gas content, to obtain several adjacent soil regions of the underground water body.

[0063] It is understandable that the extraction of groundwater resources from underground water bodies cannot solely rely on their water content; the soil structure of the surrounding area must also be considered. Some groundwater areas are surrounded by soil with a fragile structure, a low solid content, a high gas content, and large internal voids. Extracting groundwater from such soils reduces the liquid content, making it difficult to support the soil's weight, easily leading to ground subsidence and damaging the ecological environment. Therefore, it is necessary to further assess the environmental stability of the groundwater area in conjunction with the surrounding soil conditions. The higher the gas content within the adjacent soil area of ​​a groundwater area, and the more concentrated the high gas content in localized areas, the larger the internal voids in that adjacent soil area and the more prone it is to collapse, making the surrounding groundwater environment more unstable.

[0064] like Figure 2 As shown, the adjacent soil regions around the underground water body can be identified based on the groundwater (liquid phase), non-groundwater substance (solid phase), and gas (gas phase) content at various points around the underground water body.

[0065] Step 103: Determine the environmental stability index of the groundwater area based on the total groundwater content and the instability of each adjacent soil region.

[0066] Among them, the environmental stability index is used to characterize the environmental stability of underground water bodies.

[0067] In one embodiment, the degree of influence of each adjacent soil region on the environmental stability of the groundwater is determined. Based on the degree of influence of each adjacent soil region on the environmental stability of the groundwater, the instability of each adjacent soil region is weighted and summed to obtain the instability of the surrounding soil of the groundwater. Based on the total groundwater content of the groundwater and the instability of the surrounding soil, the environmental stability index of the groundwater is determined. The environmental stability index of the groundwater is positively correlated with the total groundwater content of the groundwater and negatively correlated with the instability of the surrounding soil.

[0068] Step 104: Based on the data on changes in groundwater content and environmental stability index of the underground water area, invert the future groundwater content of the underground water area.

[0069] The data on changes in groundwater content in underground water bodies includes the groundwater content at each specific moment, which can be obtained by satellites that observe changes in the Earth's gravitational field. For example, it can be obtained through the GRACE satellite.

[0070] The aforementioned groundwater storage inversion method, which considers surface soil deformation, determines several underground water bodies and the total groundwater content of each underground water body based on the groundwater content at various underground points. For each underground water body, the instability of each adjacent soil region is determined based on the gas content at each point in the adjacent soil region. Based on the total groundwater content of the underground water body and the instability of each adjacent soil region, the environmental stability index of the underground water body is determined. This method takes into account the surrounding soil conditions and their impact on the underground water body. Based on the changes in groundwater content and the environmental stability index, the future groundwater content of the underground water body is inverted. This method can offset the impact of human mining activities on the accuracy of the inversion results and improve the accuracy of groundwater storage inversion.

[0071] In one embodiment, see Figure 3 Step 101 determines the groundwater content of several underground water bodies and the total groundwater content of each underground water body based on the groundwater content at various underground points, including the following steps:

[0072] Step 1011: Determine several underground water bodies based on the groundwater content at various underground points.

[0073] Step 1012: For each underground water body, determine the location coordinates of the representative point of the underground water body based on the groundwater content at each point in the underground water body.

[0074] In one embodiment, extreme points are determined based on the groundwater content at each point in the underground water body, and then the location coordinates of representative points in the underground water body are determined based on the groundwater content and location coordinates of each extreme point.

[0075] In one embodiment, the extreme point can be the point with the highest groundwater content within a neighborhood of a preset size. For example, the extreme point can be the point with the highest groundwater content within a 27-neighborhood.

[0076] Step 1013: Determine the total groundwater content of the underground water body based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, as well as the groundwater content of each point.

[0077] In the above embodiments, several underground water bodies are determined based on the groundwater content at various underground points. Then, based on the groundwater content at each point in the underground water bodies, the location coordinates of the representative point of the underground water body are determined. Based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, as well as the groundwater content at each point, the total groundwater content of the underground water bodies can be accurately determined.

[0078] In one embodiment, determining several underground water bodies based on the groundwater content at various underground points includes: selecting a set of target points from each point based on the groundwater content; randomly selecting a point from the unmarked points in the target point set as a seed point each time; starting from the seed point, performing region growth based on the groundwater content of each point in the neighboring area of ​​the seed point to obtain the underground water body corresponding to the seed point, and marking each point in the underground water body; until all points in the target point set are marked as underground water bodies, resulting in several underground water bodies.

[0079] In the above embodiments, since each underground point contains a certain amount of liquid phase in its local area, the groundwater content in the underground water area is relatively large compared to other areas, and the change in groundwater content within the same underground water area is progressive and non-abrupt. Therefore, several underground water areas can be accurately obtained through regional growth.

[0080] In one embodiment, selecting a set of target points from the groundwater content at each underground point includes: calculating the average groundwater content at each underground point as a groundwater content threshold; and identifying points with groundwater content greater than or equal to the groundwater content threshold as target points to form a set of target points.

[0081] The average groundwater content at various underground points can be calculated using the following formula:

[0082]

[0083] in, This represents the average groundwater content at various underground points, i.e., the groundwater content threshold. Indicates the underground The groundwater content at each point. This represents the total number of points in the underground 3D point cloud data.

[0084] In the above embodiments, the average groundwater content at each underground point is calculated and used as a groundwater content threshold to segment the point cloud to obtain a target point set. This allows for the accurate selection of target points with relatively higher groundwater content, thereby enabling the accurate determination of underground water bodies based on the target points.

[0085] In one embodiment, determining the location coordinates of a representative point in the underground water body based on the groundwater content at each point in the underground water body includes: determining extreme points based on the groundwater content at each point in the underground water body; determining a first weight for each extreme point based on the groundwater content at each extreme point; and weighting and summing the location coordinates of each extreme point based on the first weight of each extreme point to obtain the location coordinates of the representative point in the underground water body.

[0086] In one embodiment, the sum of groundwater content at all extreme points in the underground water body is calculated, and the first weight corresponding to each extreme point is determined based on the ratio between the groundwater content at each extreme point and the sum of groundwater content.

[0087] In one embodiment, based on the first weight of each extreme point, the x-coordinate, y-coordinate, and z-coordinate values ​​in the location coordinates of each extreme point are weighted and summed to obtain the x-coordinate, y-coordinate, and z-coordinate values ​​in the location coordinates of the representative point of the underground water area.

[0088] For example, the x-coordinate value of the location coordinates of a representative point in an underground water body can be determined using the following formula:

[0089]

[0090] in, Indicates the first The x-coordinate value of the location coordinates of the representative point of each underground water area. Indicates the first The first underground water area The groundwater content at each extreme point. Indicates the first The sum of groundwater content at all extreme points within a single underground water body, i.e. . Indicates the first The first underground water area The first weight of each extreme point. Indicates the first The total number of extreme points within a single underground water body. Indicates the first The first underground water area The x-coordinate value of the position coordinates of each extreme point.

[0091] Similarly, the y-coordinate and z-coordinate values ​​of the representative point in the underground water area can be calculated, and denoted as follows: and .

[0092] In the above embodiments, since the location coordinates of extreme points with higher groundwater content have a greater impact on the location coordinates of representative points in the underground water area, extreme points are determined based on the groundwater content of each point in the underground water area. The first weight of each extreme point is determined based on the groundwater content of each extreme point. The location coordinates of each extreme point are weighted and summed based on the first weight of each extreme point to obtain the location coordinates of representative points in the underground water area. This method can accurately determine the location coordinates of representative points in the underground water area.

[0093] In one embodiment, determining the total groundwater content of the underground water body based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, and the groundwater content of each point, includes: determining a second weight for each point in the underground water body based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point; and weighting and summing the groundwater content of each point in the underground water body based on the second weight to obtain the total groundwater content of the underground water body.

[0094] In one embodiment, the second weight of a point in the groundwater area is negatively correlated with the distance between the point's location coordinates and the location coordinates of a representative point. That is, the closer the distance, the greater the second weight. A larger second weight indicates a greater correlation between the groundwater content at that point and the total groundwater content of the groundwater area.

[0095] In one embodiment, the sum of the distances between the position coordinates of each point in the underground water area and the position coordinates of the representative point can be calculated. Then, the difference between the sum of the distances and the corresponding distances of each point in the underground water area can be calculated. Then, the differences corresponding to each point in the underground water area can be normalized to obtain the second weight corresponding to each point (the sum of the second weights corresponding to each point is equal to 1).

[0096] In one embodiment, the total groundwater content of an underground water body can be determined according to the following formula:

[0097]

[0098] in, Indicates the first The total groundwater content of each underground water body. Indicates the first The distance between the location coordinates of the b-th point in the underground water area and the location coordinates of the representative point. Indicates the first The sum of the distances from the coordinates of each point in the underground water area to the coordinates of the representative point, i.e. . Indicates the first The total number of points in the underground water area. Indicates the first The second weight of the b-th point in the underground water area. Indicates the first The groundwater content at point b in the underground water area.

[0099] In the above embodiments, since the total groundwater content of an underground water body is directly related to the groundwater content of each point within it, the closer to the representative point of the underground water body, the higher the correlation between the groundwater content of that point and the total groundwater content of the underground water body. Therefore, based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, the second weight of each point in the underground water body is determined. Based on the second weight of each point in the underground water body, the groundwater content of each point in the underground water body is weighted and summed to accurately determine the total groundwater content of the underground water body.

[0100] In one embodiment, see Figure 4 Step 102 determines the instability of each adjacent soil region based on the gas content at each point in each adjacent soil region of the groundwater, including the following steps:

[0101] Step 1021: For each adjacent soil region of the underground water body, determine the average gas content at each point in the adjacent soil region.

[0102] Step 1022: Determine the number of gas phase target points in adjacent soil regions whose gas content is greater than or equal to the average gas content.

[0103] Among them, the gas phase target point is the point in the adjacent soil region where the gas content is greater than or equal to the average gas content.

[0104] Step 1023: Determine the instability of adjacent soil areas based on the average gas content, the number of gas phase target points, and the distance between each gas phase target point.

[0105] In one embodiment, the two closest gas phase target points can be selected sequentially from the gas phase target points, and the distance between these two gas phase target points can be calculated. This process continues until all gas phase target points have been selected at least once. The sum of the distances between all selected pairs of gas phase target points is then calculated and denoted as the total distance. The instability of adjacent soil regions is determined based on the mean gas content, the number of gas phase target points, and the total distance. The instability of adjacent soil regions is positively correlated with the mean gas content and the number of gas phase target points, and negatively correlated with the total distance.

[0106] In one embodiment, the product of the mean gas content at each point in adjacent soil regions and the number of gas phase target points is calculated, and then the product is divided by the sum of distances to obtain the instability level of adjacent soil regions. The formula is as follows:

[0107]

[0108] in, Indicates the first The first underground water area The degree of instability of adjacent soil regions. Indicates the first The first underground water area The average gas content at each point in adjacent soil regions. Indicates the first The first underground water area The number of gas phase target points in adjacent soil regions whose gas content is greater than or equal to the average gas content. Indicates the first The first underground water area The sum of the distances between adjacent soil regions.

[0109] Mean gas content at points in adjacent soil regions The number of gas phase target points in adjacent soil regions with gas content greater than or equal to the average gas content. The larger the value, the higher the gas content and the larger the soil gaps within the adjacent soil region, thus indicating a higher degree of instability in the adjacent soil region. The larger; the sum of the distances between adjacent soil regions The smaller the value, the more concentrated the local areas with higher gas content are within adjacent soil regions, thus indicating a higher degree of instability in adjacent soil regions. The larger the value, the greater the instability of the adjacent soil region. The larger the value, the higher the instability of the adjacent soil area, and the more unstable the environment of the groundwater.

[0110] In the above embodiments, for each adjacent soil region of the underground water body, the average gas content, the number of gas phase target points, and the distance between each gas phase target point are determined for each point in the adjacent soil region, which can accurately determine the degree of instability of the adjacent soil region.

[0111] In one embodiment, the environmental stability index of the groundwater is determined based on the total groundwater content and the instability of each adjacent soil region. This includes: determining the degree of influence of each adjacent soil region on the environmental stability of the groundwater based on the distance between the centroid of each adjacent soil region and the representative point of the groundwater, and the number of points in each adjacent soil region; weighting and summing the instability of each adjacent soil region based on the degree of influence of each adjacent soil region on the environmental stability of the groundwater to obtain the instability of the surrounding soil of the groundwater; and determining the environmental stability index of the groundwater based on the total groundwater content and the instability of the surrounding soil.

[0112] In one embodiment, the centroid coordinates of each adjacent soil region can be determined using the point mass method. Based on the centroid coordinates and the coordinates of a representative point of the groundwater, the distance between the centroid of each adjacent soil region and the representative point of the groundwater can be determined.

[0113] In one embodiment, the degree of influence of adjacent soil regions on the environmental stability of groundwater is negatively correlated with the distance between the centroid of the adjacent soil region and the representative point of the groundwater, and positively correlated with the number of points in the adjacent soil region.

[0114] In one embodiment, the degree of influence of an adjacent soil region on the environmental stability of the groundwater can be determined based on the ratio of the number of points in an adjacent soil region to the distance between the centroid of that adjacent soil region and a representative point in the groundwater. The formula is as follows:

[0115]

[0116] in, Indicates the first The first underground water area The degree of influence of adjacent soil regions on the environmental stability of groundwater. Indicates the first The first underground water area The number of points in the three-dimensional point cloud of adjacent soil regions. Indicates the first The first underground water area The centroid of the adjacent soil regions and the first The distance between representative points of the underground water bodies. This represents the minimum value hyperparameter, used to avoid the case where the denominator is 0.

[0117] Number of points in the 3D point cloud of adjacent soil regions The more [specific soil area] there is, the larger the adjacent soil region, and therefore the greater the impact on the environmental stability of the groundwater. The larger the distance between the centroid of adjacent soil regions and the representative point of the groundwater. The smaller the value, the closer the adjacent soil area is to the groundwater, and therefore the less impact it has on the environmental stability of the groundwater. The larger.

[0118] In one embodiment, the sum of the impacts of each adjacent soil region on the environmental stability of the groundwater is calculated. Based on the ratio between the impact of each adjacent soil region on the environmental stability of the groundwater and the sum of the impacts, a third weight is determined for each adjacent soil region. The instability of each adjacent soil region is weighted and summed based on the third weights for each adjacent soil region to obtain the instability of the soil surrounding the groundwater.

[0119] In one embodiment, the environmental stability index of the groundwater is positively correlated with the total groundwater content of the groundwater and negatively correlated with the instability of the surrounding soil.

[0120] In one embodiment, the environmental stability index of a groundwater body can be determined based on the ratio of its total groundwater content to the instability of the surrounding soil. The formula is as follows:

[0121]

[0122] in, Indicates the first An environmental stability index for underground water bodies. Indicates the first The total groundwater content of each underground water body. Indicates the first The first underground water area The degree of instability of adjacent soil regions. Indicates the first The first underground water area The degree of influence of adjacent soil regions on the environmental stability of groundwater. Indicates the first The sum of the impacts of all adjacent soil regions on the environmental stability of a groundwater body, i.e. . Indicates the first The total number of adjacent soil areas of each underground water body. This indicates the degree of instability of the soil surrounding underground water bodies.

[0123] Total groundwater content in underground water bodies The larger the value, the richer the water content in the underground water body, and the more suitable it is for extraction. The instability of the adjacent soil area of ​​the underground water body. The larger the value, the higher the instability of the adjacent soil area. This means that mining the groundwater is more likely to cause ground subsidence in the surrounding area, and the more unstable the environment of the groundwater is, the less suitable it is for mining. The overall soil instability around the groundwater is obtained by weighted summing the instability values ​​of each adjacent soil area based on their impact on the environmental stability of the groundwater. .therefore Characterizing the first The degree of influence of adjacent soil regions on the environmental stability of a groundwater body is used to determine the total groundwater content of the groundwater body. After adjustments, the final environmental stability index of the underground water body is obtained. Environmental stability index of underground water The larger the value, the higher the stability of the underground water environment, and the more likely it is to be mined.

[0124] In the above embodiments, based on the distance between the centroid of each adjacent soil region and the representative point of the groundwater, and the number of points in each adjacent soil region, the degree of influence of each adjacent soil region on the environmental stability of the groundwater can be accurately determined. Then, based on the total groundwater content of the groundwater and the instability of the surrounding soil, the environmental stability index of the groundwater can be accurately determined.

[0125] In one embodiment, the method of retrieving the future groundwater content of a groundwater body based on groundwater content change data and an environmental stability index includes: determining the initial groundwater content change rate at the current moment based on the groundwater content change data; weighting the initial groundwater content change rate according to the environmental stability index of the groundwater body to obtain the groundwater content change rate at the current moment; and retrieving the future groundwater content based on the groundwater content change rate at the current moment.

[0126] The data on changes in groundwater content includes the groundwater content at various points in time.

[0127] In one embodiment, data on changes in groundwater content in underground water bodies can be obtained using the GRACE satellite.

[0128] In one embodiment, the curve of the change in groundwater content in the underground water body can be fitted using the least squares method to obtain the slope of the curve at the current moment, which can be used as the initial rate of change of groundwater content in the underground water body at the current moment.

[0129] Assume the rate of change of the initial groundwater content in the underground water body at the current moment is denoted as . Based on the environmental stability index of the groundwater area, the initial groundwater content change rate is weighted to obtain the groundwater content change rate at the current moment, which can be expressed as: ,in, This indicates the environmental stability index of underground water bodies.

[0130] In the above embodiments, by weighting the initial groundwater content change rate according to the environmental stability index of the underground water body, a more accurate groundwater content change rate at the current moment can be obtained. Based on the groundwater content change rate at the current moment, the future groundwater content of the underground water body can be accurately inverted.

[0131] See Figure 5 This invention provides a groundwater storage inversion system that considers surface soil deformation. The system includes a memory and a processor. The memory stores executable program code. The processor calls and runs the executable program code from the memory to perform the following steps: determining several underground water bodies and the total groundwater content of each underground water body based on the groundwater content at various points underground; determining the instability of each adjacent soil region for each underground water body based on the gas content at each point in each adjacent soil region; determining the environmental stability index of the underground water body based on the total groundwater content and the instability of each adjacent soil region; and inverting the future groundwater content of the underground water body based on the groundwater content change data and the environmental stability index.

[0132] In one embodiment, determining a number of underground water bodies and the total groundwater content of each underground water body based on the groundwater content at various underground points includes: determining a number of underground water bodies based on the groundwater content at various underground points; determining the location coordinates of a representative point of each underground water body based on the groundwater content at each point in the underground water body; and determining the total groundwater content of the underground water body based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, as well as the groundwater content at each point.

[0133] In one embodiment, determining several underground water bodies based on the groundwater content at various underground points includes: selecting a set of target points from each point based on the groundwater content; randomly selecting a point from the unmarked points in the target point set as a seed point each time; starting from the seed point, performing region growth based on the groundwater content of each point in the neighboring area of ​​the seed point to obtain the underground water body corresponding to the seed point, and marking each point in the underground water body; until all points in the target point set are marked as underground water bodies, resulting in several underground water bodies.

[0134] In one embodiment, selecting a set of target points from the groundwater content at each underground point includes: calculating the average groundwater content at each underground point as a groundwater content threshold; and identifying points with groundwater content greater than or equal to the groundwater content threshold as target points to form a set of target points.

[0135] In one embodiment, determining the location coordinates of a representative point in the underground water body based on the groundwater content at each point in the underground water body includes: determining extreme points based on the groundwater content at each point in the underground water body; determining a first weight for each extreme point based on the groundwater content at each extreme point; and weighting and summing the location coordinates of each extreme point based on the first weight of each extreme point to obtain the location coordinates of the representative point in the underground water body.

[0136] In one embodiment, determining the total groundwater content of the underground water body based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, and the groundwater content of each point, includes: determining a second weight for each point in the underground water body based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point; and weighting and summing the groundwater content of each point in the underground water body based on the second weight to obtain the total groundwater content of the underground water body.

[0137] In one embodiment, determining the instability of each adjacent soil region based on the gas content at each point in each adjacent soil region of the groundwater includes: determining the mean gas content at each point in each adjacent soil region of the groundwater; determining the number of gas phase target points in the adjacent soil regions whose gas content is greater than or equal to the mean gas content; and determining the instability of the adjacent soil regions based on the mean gas content, the number of gas phase target points, and the distance between each gas phase target point.

[0138] In one embodiment, the environmental stability index of the groundwater is determined based on the total groundwater content and the instability of each adjacent soil region. This includes: determining the degree of influence of each adjacent soil region on the environmental stability of the groundwater based on the distance between the centroid of each adjacent soil region and the representative point of the groundwater, and the number of points in each adjacent soil region; weighting and summing the instability of each adjacent soil region based on the degree of influence of each adjacent soil region on the environmental stability of the groundwater to obtain the instability of the surrounding soil of the groundwater; and determining the environmental stability index of the groundwater based on the total groundwater content and the instability of the surrounding soil.

[0139] In one embodiment, the method of retrieving the future groundwater content of a groundwater body based on groundwater content change data and an environmental stability index includes: determining the initial groundwater content change rate at the current moment based on the groundwater content change data; weighting the initial groundwater content change rate according to the environmental stability index of the groundwater body to obtain the groundwater content change rate at the current moment; and retrieving the future groundwater content based on the groundwater content change rate at the current moment.

[0140] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0141] The embodiments described above are merely examples of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these modifications and improvements all fall within the protection scope of this application.

[0142] It should be noted that the order of the above embodiments of the present invention is merely for descriptive purposes and does not represent the superiority or inferiority of the embodiments. The processes depicted in the accompanying drawings do not necessarily require a specific or sequential order to achieve the desired result. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0143] The various embodiments in this specification are described in a progressive manner. The same or similar parts between the various embodiments can be referred to each other. Each embodiment focuses on describing the differences from other embodiments.

Claims

1. A method for groundwater storage inversion considering surface soil deformation, characterized in that, The method includes: Based on the groundwater content at various underground points, determine several underground water bodies and the total groundwater content of each of the underground water bodies; For each of the underground water bodies, the instability of each adjacent soil region is determined based on the gas content at each point in each adjacent soil region of the underground water body. The environmental stability index of the underground water body is determined based on the total groundwater content of the underground water body and the instability of each of the adjacent soil regions. Based on the data on changes in groundwater content in the underground water area and the environmental stability index, the future groundwater content of the underground water area can be inverted. The determination of the instability level of each adjacent soil region based on the gas content at each point in each adjacent soil region of the underground water area includes: For each adjacent soil region of the underground water body, the average gas content at each point in the adjacent soil region is determined; Determine the number of gas phase target points in the adjacent soil regions whose gas content is greater than or equal to the average gas content; The instability of the adjacent soil regions is determined based on the average gas content, the number of gas phase target points, and the distance between each gas phase target point. The determination of the environmental stability index of the underground water body based on the total groundwater content and the instability of each of the adjacent soil regions includes: The degree of influence of each adjacent soil region on the environmental stability of the groundwater is determined based on the distance between the centroid of each adjacent soil region and the representative point of the groundwater, as well as the number of points in each adjacent soil region. Based on the degree of influence of each of the adjacent soil regions on the environmental stability of the groundwater, the instability of each of the adjacent soil regions is weighted and summed to obtain the degree of soil instability surrounding the groundwater. The environmental stability index of the underground water body is determined based on the total groundwater content and the instability of the surrounding soil.

2. The groundwater storage inversion method considering surface soil deformation according to claim 1, characterized in that, The determination of several underground water bodies and the total groundwater content of each underground water body based on the groundwater content at various underground points includes: Based on the groundwater content at various underground points, several underground water bodies were identified; For each of the underground water bodies, the location coordinates of a representative point of the underground water body are determined based on the groundwater content at each point in the underground water body. The total groundwater content of the underground water body is determined based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, as well as the groundwater content at each point.

3. The groundwater storage inversion method considering surface soil deformation according to claim 2, characterized in that, Based on the groundwater content at various underground points, several underground water bodies are identified, including: Based on the groundwater content at each underground point, a set of target points is selected from these points; Each time, a point is randomly selected from the unmarked points in the target point set as a seed point. Starting from the seed point, regional growth is performed based on the groundwater content of each point in the neighborhood of the seed point to obtain the underground water area corresponding to the seed point, and each point in the underground water area is marked. This continues until all points in the target point set are marked to underground water bodies, resulting in several underground water bodies.

4. The groundwater storage inversion method considering surface soil deformation according to claim 3, characterized in that, The selection of a target point set based on the groundwater content at each underground point includes: The average groundwater content at each underground point is calculated and used as the groundwater content threshold. Points with groundwater content greater than or equal to the groundwater content threshold are identified as target points, forming a set of target points.

5. The groundwater storage inversion method considering surface soil deformation according to claim 2, characterized in that, The step of determining the location coordinates of representative points in the underground water body based on the groundwater content at each point in the underground water body includes: The extreme points are determined based on the groundwater content at each point in the underground water area; Based on the groundwater content at each extreme point, a first weight is determined for each extreme point. Based on the first weight of each extreme point, the position coordinates of each extreme point are weighted and summed to obtain the position coordinates of the representative point of the underground water area.

6. The groundwater storage inversion method considering surface soil deformation according to claim 2, characterized in that, The step of determining the total groundwater content of the underground water body based on the distance between the location coordinates of each point in the underground water body and the location coordinates of the representative point, and the groundwater content at each point, includes: The second weight of each point in the underground water area is determined based on the distance between the position coordinates of each point in the underground water area and the position coordinates of the representative point. Based on the second weight of each point in the underground water body, the groundwater content of each point in the underground water body is weighted and summed to obtain the total groundwater content of the underground water body.

7. The groundwater storage inversion method considering surface soil deformation according to claim 1, characterized in that, The step of reversing the future groundwater content of the underground water body based on the groundwater content change data of the underground water body and the environmental stability index includes: Based on the data on the change in groundwater content in the underground water area, determine the initial rate of change in groundwater content in the underground water area at the current moment; The initial groundwater content change rate is weighted according to the environmental stability index of the underground water body to obtain the groundwater content change rate of the underground water body at the current moment. Based on the rate of change of groundwater content in the underground water body at the current moment, the future groundwater content of the underground water body can be inverted.

8. A groundwater storage inversion system considering surface soil deformation, characterized in that, The system includes a memory and a processor; the memory is used to store executable program code; the processor is used to call and run the executable program code from the memory to implement the groundwater storage inversion method considering ground soil deformation as described in any one of claims 1 to 7.

Images