Simulation method for calculating underground water balance based on three-dimensional geological data

Abstract

The application discloses a simulation method for calculating underground water balance based on three-dimensional geological data, and relates to the technical field of groundwater calculation, and comprises the following steps: constructing an independent analysis space based on an aquifer three-dimensional structure model; dividing the aquifer three-dimensional structure model into a first type space and a second type space based on the independent analysis space; constructing an in-space reference point in the second type space; obtaining a third type space and a fourth type space in the second type space based on the in-space reference point; obtaining a first proportion based on the fourth type space; calculating the collectable water amount of the whole aquifer based on the number of the first type space, the number of the third type space, the number of the fourth type space, the collectable water amount in the first type space and the first proportion; and the application is used to solve the problem that the collectable amount of groundwater cannot be accurately calculated in the existing groundwater calculation technology due to the complex aquifer structure of the superposition of the fourth system and the Cretaceous system.

Description

Technical Field

[0001] This invention relates to the field of groundwater calculation technology, specifically a groundwater balance calculation and simulation method based on three-dimensional geological data. Background Technology

[0002] Accurate assessment, rational development, and effective protection of groundwater resources are crucial for ensuring water supply security, maintaining the ecological environment, and supporting sustainable development. Therefore, it is necessary to estimate the amount of groundwater that can be extracted to prevent excessive extraction from causing damage to the ecological environment.

[0003] Traditional groundwater resource assessment and simulation methods typically rely on two-dimensional planar maps and simplified hydrogeological profiles, simplifying complex underground media into homogeneous, isotropic layered models. This method struggles to accurately depict the heterogeneity, anisotropy, and complex internal geological structures of aquifers. For example, in hilly and gully areas, aquifers exhibit fragmented structures, specifically a superposition of Quaternary and Cretaceous strata. This leads to inaccurate calculations of recoverable groundwater. Therefore, methods based on three-dimensional models have been proposed for calculating recoverable groundwater. However, due to the fragmented and complex aquifer structures in hilly and gully areas, obtaining the aquifer volume is difficult and straightforward, further complicating the calculation of recoverable groundwater. Existing groundwater calculation techniques, due to the complex aquifer structures involving the superposition of Quaternary and Cretaceous strata, cannot accurately calculate recoverable groundwater. Summary of the Invention

[0004] The present invention aims to at least partially solve one of the technical problems in the prior art. By calculating the total collectable water volume of the aquifer based on the number of first-type spaces, the number of third-type spaces, the number of fourth-type spaces, the collectable water volume in the first-type spaces, and a first ratio, the present invention addresses the problem that existing groundwater calculation techniques cannot accurately calculate the collectable groundwater volume due to the complex aquifer structure with Quaternary and Cretaceous superposition.

[0005] To achieve the above objectives, this application provides a groundwater balance calculation and simulation method based on three-dimensional geological data, comprising the following steps:

[0006] Obtain the three-dimensional structural model of the aquifer to be detected and mark it as the three-dimensional structural model of the aquifer;

[0007] An independent analysis space was constructed based on the three-dimensional structural model of the aquifer.

[0008] Based on the independent analysis space, the three-dimensional structural model of the aquifer is divided into a first type of space and a second type of space;

[0009] Construct spatial reference points within the second type of space;

[0010] Based on reference points within the space, obtain the third and fourth types of spaces in the second type of space;

[0011] The spatial average water content is obtained based on the water content in the first type of space of the third number of aquifers to be detected;

[0012] The amount of collectable water in the first type of space is obtained based on the spatial average water content.

[0013] The first proportion is obtained based on the fourth number of fourth-class spaces;

[0014] The total collectable water volume of the aquifer is calculated based on the number of first-type spaces, the number of third-type spaces, the number of fourth-type spaces, the collectable water volume in the first-type spaces, and the first ratio.

[0015] Furthermore, constructing an independent analysis space based on the three-dimensional structural model of the aquifer includes the following sub-steps:

[0016] Establish a three-dimensional coordinate system and label it as the model's three-dimensional coordinate system; place the aquifer three-dimensional structure model in the model's three-dimensional coordinate system;

[0017] The plane containing both the X-axis and Y-axis is designated as the XY plane; the plane parallel to the XY plane is designated as the first plane.

[0018] Starting with the XY plane, the interval between each first plane is the first length, and the direction is the positive and negative directions of the Z axis, and a first number of first planes are drawn;

[0019] The plane containing both the X-axis and Z-axis is designated as the XZ plane; the plane parallel to the XZ plane is designated as the second plane.

[0020] Starting with the XZ plane, the interval between each second plane is the first length, and the direction is the positive and negative directions of the Y axis, and a first number of second planes are drawn;

[0021] The plane containing both the Y-axis and Z-axis is designated as the YZ plane; the plane parallel to the YZ plane is designated as the third plane.

[0022] Starting with the YZ plane, the interval between each third plane is the first length, and the direction is the positive and negative directions of the Z axis, and the first number of third planes are drawn;

[0023] The model's three-dimensional coordinate system is divided into cubes with a side length of the first length by all the first, second, and third planes drawn, and these cubes are marked as independent analysis spaces.

[0024] Furthermore, based on the independent analysis space, the three-dimensional structural model of the aquifer is divided into a first type of space and a second type of space, including the following sub-steps:

[0025] Obtain an independent analysis space containing only three-dimensional structural models of aquifers, and label it as the first type of space;

[0026] The independent analysis space is divided into three-dimensional structural models of the aquifer, and is marked as the second type of space.

[0027] Furthermore, constructing an intra-space reference point within the second type of space includes the following sub-steps:

[0028] The top surface of the second type of space is marked as the second top surface, and the surface parallel to the second top surface is marked as the parallel top surface;

[0029] Starting with the parallel top surface, with the interval between each surface being the second length and the direction being the bottom surface of the second type of space, draw the second number of parallel top surfaces;

[0030] The side of the second type of space is labeled as the second side, and the surface parallel to the second side is labeled as the parallel side.

[0031] Starting with the second side, with the interval between each parallel side being the second length and the direction being another side of the second type of space, draw the second number of parallel sides;

[0032] The front of the second type of space is marked as the second front, and the face parallel to the second front is marked as the parallel front;

[0033] Starting with the second front, the interval between each parallel front is the second length, and the direction is the back of the second type of space, and a second number of parallel fronts are drawn;

[0034] Obtain the intersection point of the three surfaces parallel to the top surface, parallel to the side surface, and parallel to the front surface in the second type of space, and mark it as the reference point in the space.

[0035] Furthermore, obtaining the third and fourth types of spaces within the second type of space based on reference points within the space includes the following sub-steps:

[0036] With a reference point in space as the center and a radius of the third length, draw a sphere and mark it as the sphere reference space;

[0037] Obtain a spherical reference space containing only three-dimensional structural models of aquifers, and label it as the third type of space;

[0038] The sphere reference space is obtained by dividing the interior of the aquifer into three-dimensional structural models, and is marked as the fourth type of space.

[0039] Furthermore, obtaining the spatial average water content based on the water content within the first type of space of the third number of aquifers to be detected includes the following sub-steps:

[0040] Obtain the water content in the first type of space of the third number of aquifers to be tested, and mark it as the experimental water content;

[0041] Obtain the range of experimental water content, divide the range of experimental water content into k equally spaced intervals, and mark them as water content division intervals;

[0042] Count the frequency of each moisture content division interval and label it as the moisture content division frequency;

[0043] A histogram was plotted with the experimental moisture content as the X-axis, the frequency of moisture content division as the Y-axis, and the intervals of moisture content division as the histogram intervals. This histogram was then labeled as the moisture content histogram.

[0044] Furthermore, obtaining the spatial average water content based on the water content within the first type of space of a third number of aquifers to be detected also includes the following sub-steps:

[0045] The water content frequency threshold is calculated as: Psh = u × g3 ÷ k; where Psh is the water content frequency threshold, u is a value between 0 and 1, and g3 is the third quantity;

[0046] Water content less than or equal to the water content frequency threshold is classified as an abnormal frequency.

[0047] Delete the water content division intervals with the leftmost and rightmost abnormal division frequencies in the water content histogram, and obtain the minimum and maximum values ​​of the water content division intervals in the water content histogram after deletion, and mark them as the first water content threshold and the second water content threshold.

[0048] The mean of experimental water content that is greater than or equal to the first water content threshold and less than or equal to the second water content threshold is calculated and marked as the spatial average water content.

[0049] Furthermore, obtaining the collectable water volume within the first type of space based on the spatial average water content includes the following sub-steps:

[0050] Obtain the lowest safe water content within the first type of space, denoted as Ha;

[0051] The collectable water volume in each first-type space is calculated as: Hc = Hp - Ha; where Hc is the collectable water volume in each first-type space, and Hp is the average water content of the space.

[0052] Furthermore, obtaining the first proportion based on the fourth number of fourth-class spaces includes the following sub-steps:

[0053] Obtain the fourth number of fourth-class spaces and mark them as historical spaces;

[0054] Obtain the volume of the three-dimensional structural model of the aquifer in historical space and mark it as the historical volume;

[0055] Obtain the total volume of historical space and the total volume of all historical spaces;

[0056] The ratio of the total historical volume to the total historical spatial volume is calculated and marked as the first ratio.

[0057] Furthermore, calculating the total collectable water volume of the aquifer based on the number of first-type spaces, the number of third-type spaces, the number of fourth-type spaces, the collectable water volume within the first-type spaces, and the first ratio includes the following sub-steps:

[0058] Get the number of first-class spaces, labeled D1; get the number of all third-class spaces, labeled D3; get the number of all fourth-class spaces, labeled D4; get the number of reference points within each second-class space, labeled D5;

[0059] The total collectable water volume of the aquifer is calculated as: Hz = D1 × Hc + (D3 + D4 × b1) / D5 × Hc; where Hz is the total collectable water volume of the aquifer and b1 is the first proportion.

[0060] The beneficial effects of the present invention are as follows: The present invention calculates the total collectable water volume of the aquifer based on the number of first-type spaces, the number of third-type spaces, the number of fourth-type spaces, the collectable water volume in the first-type spaces, and the first ratio. The advantage is that it can cope with the complex aquifer structure of Quaternary and Cretaceous superposition and improve the accuracy of the calculated collectable groundwater volume.

[0061] This invention obtains the third and fourth types of space in the second type of space based on reference points within the space. Its advantage lies in its ability to cope with complex aquifer structures, optimize the acquisition of aquifer volume, and improve the accuracy of aquifer volume acquisition. Attached Figure Description

[0062] Figure 1 This is a flowchart illustrating the steps of the method of the present invention;

[0063] Figure 2 This is a schematic diagram of the independent analysis space of the present invention;

[0064] Figure 3 This is a schematic diagram of the spatial reference point of the present invention;

[0065] Figure 4 This is a schematic diagram of the water content histogram of the present invention. Detailed Implementation

[0066] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0067] Example 1, please refer to Figure 1 As shown, this application provides a groundwater balance calculation and simulation method based on three-dimensional geological data, including the following steps:

[0068] Step S1: Obtain the three-dimensional structural model of the aquifer to be detected and mark it as the three-dimensional structural model of the aquifer; obtain the three-dimensional structural model of the aquifer based on existing instruments.

[0069] Step S2 involves constructing an independent analysis space based on the three-dimensional structural model of the aquifer. Step S2 includes the following sub-steps:

[0070] Step S201: Establish a three-dimensional coordinate system and label it as the model's three-dimensional coordinate system; place the aquifer three-dimensional structure model in the model's three-dimensional coordinate system;

[0071] Step S202: The plane containing both the X-axis and Y-axis is marked as the XY plane; the plane parallel to the XY plane is marked as the first plane;

[0072] Step S203: Starting from the XY plane, the interval between each first plane is the first length, and the direction is the positive and negative directions of the Z-axis, drawing a first number of first planes; due to the complex aquifer structure of the superimposed Quaternary and Cretaceous, the volume of the aquifer cannot be directly obtained based on the three-dimensional structure model of the aquifer, so an independent analysis space is defined. Based on the volume of the three-dimensional structure model of the aquifer obtained by the independent analysis, the first length is the side length of the independent analysis space, so the first length should not be too large, for example, the first length is 2 meters; the defined independent analysis space should include all the three-dimensional structure models of the aquifer, so the first number is set based on the first length and the size of the three-dimensional structure model of the aquifer, for example, the first number is 100.

[0073] Step S204: The plane in which both the X-axis and Z-axis lie is marked as the XZ plane; the plane parallel to the XZ plane is marked as the second plane;

[0074] Step S205: Starting with the XZ plane, the interval between each second plane is the first length, and the direction is the positive and negative directions of the Y axis, draw the first number of second planes;

[0075] Step S206: The plane containing both the Y-axis and Z-axis is marked as the YZ plane; the plane parallel to the YZ plane is marked as the third plane;

[0076] Step S207: Starting with the YZ plane, the interval between each third plane is the first length, and the direction is the positive and negative directions of the Z axis, draw the first number of third planes;

[0077] Step S208: The three-dimensional coordinate system of the model is divided into cubes with a side length of the first length by drawing all the first, second and third planes, and marked as independent analysis space;

[0078] For practical applications, please refer to Figure 2 As shown, an independent analysis space is plotted.

[0079] Step S3: Based on the independent analysis space, the three-dimensional structural model of the aquifer is divided into a first type of space and a second type of space; Step S3 includes the following sub-steps:

[0080] Step S301: Obtain the independent analysis space containing only three-dimensional structural models of aquifers, and mark it as the first type of space;

[0081] Step S302: Obtain the independent analysis space, which is divided into three-dimensional structural models of the aquifer, and mark it as the second type of space; because the complex aquifer structure of the Quaternary and Cretaceous superposition will result in the second type of space.

[0082] The first type of space has a cubic structure, so the volume of the three-dimensional structure model of the water layer in the first type of space is easy to obtain. The second type of space is an incomplete three-dimensional structure model of the aquifer, so the three-dimensional structure model of the water layer in the second type of space needs further analysis. Therefore, it is divided into the first type of space and the second type of space.

[0083] Step S4: Construct a spatial reference point within the second type of space; Step S4 includes the following sub-steps:

[0084] Step S401: Mark the top surface of the second type of space as the second top surface, and mark the surface parallel to the second top surface as the parallel top surface;

[0085] Step S402: Starting with the parallel top surface, the interval between each surface is the second length, and the direction is the bottom surface of the second type of space, draw a second number of parallel top surfaces; where the second length is the interval between the spherical reference spaces in subsequent steps, and all spherical reference spaces must be guaranteed to be within the second type of space, so the second length must be much smaller than the first length, for example, the second length is 0.4m; the second number is set based on the first length and the second length. Since the first length is 2m and the second length is 0.4m, in order to ensure that the reference points in the space are within the second type of space, the second number is 4;

[0086] Step S403: Mark the side of the second type of space as the second side and mark the surface parallel to the second side as the parallel side.

[0087] Step S404: Starting with the second side, with the interval between each parallel side being the second length and the direction being another side of the second type of space, draw the second number of parallel sides.

[0088] Step S405: Mark the front of the second type of space as the second front, and mark the surface parallel to the second front as the parallel front;

[0089] Step S406: Starting with the second front, the interval between each parallel front is the second length, and the direction is the back of the second type of space, draw the second number of parallel fronts;

[0090] Step S407: Obtain the intersection point of the three surfaces parallel to the top surface, parallel to the side surface, and parallel to the front surface in the second type of space, and mark them as reference points in the space; constructing reference points in the space is to facilitate the subsequent establishment of the spherical reference space;

[0091] For practical applications, please refer to Figure 3 As shown, a reference point within a second-class space.

[0092] Step S5: Obtain the third and fourth type spaces in the second type of space based on reference points within the space; Step S5 includes the following sub-steps:

[0093] Step S501: Using the reference point in space as the center and the radius as the third length, draw a sphere and mark it as the sphere reference space. The third length must be set to ensure that the sphere reference space does not exceed the second type of space and that the sphere reference spaces do not overlap. Therefore, the third length should be less than half of the second length, for example, the third length is 0.06m. At the same time, the sphere reference space established based on the reference point in space can make the sphere reference space evenly distributed in the second type of space. Establishing the sphere reference space can not only analyze the volume of the aquifer in the second type of space relatively accurately, but also reduce the amount of calculation by dividing the second type of space completely.

[0094] Step S502: Obtain the sphere reference space, which is entirely composed of three-dimensional aquifer structure models, and mark it as the third type of space;

[0095] Step S503: Obtain the spherical reference space, which is divided into a three-dimensional structural model of the aquifer, and mark it as the fourth type of space; here, the third type of space and the fourth type of space are divided in order to further analyze the volume of the aquifer in the second type of space.

[0096] Step S6: Obtain the spatial average water content based on the water content within the first type of space of the third number of aquifers to be detected; Step S6 includes the following sub-steps:

[0097] Step S601: Obtain the water content in the first type of space of the third number of aquifers to be tested, and mark it as the experimental water content; the experimental water content can be obtained by sampling the aquifers to be tested; the third number is obtained in order to reduce the error of the average water content obtained, for example, the third number is 297.

[0098] Step S602: Obtain the range of experimental moisture content, divide the range of experimental moisture content into k equally spaced intervals, and mark them as moisture content division intervals; the setting of k is to obtain the distribution of experimental moisture content, so the setting of k should not be too large or too small, for example, k is 9;

[0099] Step S603: Count the frequency of each moisture content division interval and mark it as the moisture content division frequency;

[0100] Step S604: Plot a histogram with the experimental moisture content as the X-axis, the frequency of moisture content division as the Y-axis, and the moisture content division intervals as the histogram interval intervals, and mark it as the moisture content histogram.

[0101] For practical applications, please refer to Figure 4 As shown, the water content histogram was drawn.

[0102] Step S605, calculate the water content frequency threshold as: Psh=u×g3÷k; where Psh is the water content frequency threshold, u is a value between 0 and 1, and g3 is the third quantity; the water content frequency threshold is set in order to filter out the smaller water content and divide the frequency, so the proportion of u is set to be small, for example, u is 0.1;

[0103] Step S606: Classify water content frequencies that are less than or equal to the water content frequency threshold as abnormal classification frequencies;

[0104] Step S607: Delete the water content division intervals with abnormal frequency divisions on the leftmost and rightmost sides of the water content histogram, and obtain the minimum and maximum values ​​of the water content division intervals in the water content histogram after deletion, and mark them as the first water content threshold and the second water content threshold.

[0105] Step S608: Calculate the average experimental water content that is greater than or equal to the first water content threshold and less than or equal to the second water content threshold, and mark it as the spatial average water content; the spatial average water content is the average water content within each first type of space.

[0106] For practical applications, please refer to Figure 4As shown, the calculated water content frequency threshold is: Psh = u × g³ ÷ k = 0.1 × 297 ÷ 9 = 3.3. Water content intervals less than or equal to 3.3 are marked as abnormal intervals. The leftmost water content interval with a frequency of 3 in the water content histogram is also an abnormal interval, and the corresponding interval is deleted (1.2m). 3 up to 1.3m 3 If the rightmost water content interval in the water content histogram has a frequency of 3, it is considered an abnormal interval, and the corresponding water content interval 2.1m should be deleted. 3 up to 2.2m 3 The minimum and maximum values ​​of the water content intervals in the water content histogram after deletion are 1.3m. 3 With 2.2m 3 Therefore, the first and second water cut thresholds are 1.3m. 3 With 2.2m 3 Find the value greater than or equal to 1.3m. 3 At the same time, less than or equal to 2.2m 3 The average water content in the experiment was 1.8 m. 3 Therefore, the average water content in the space is 1.8m. 3 .

[0107] Step S7: Obtain the collectable water volume within the first type of space based on the spatial average water content; Step S7 includes the following sub-steps:

[0108] Step S701: Obtain the lowest safe water content within the first type of space, denoted as Ha; the safe water content is the groundwater that can maintain normal ecology and geological structure, for example, Ha is 1.6m. 3 ;

[0109] Step S702, calculate the collectable water volume in each first-type space as: Hc = Hp - Ha; where Hc is the collectable water volume in each first-type space, Hp is the average water content of the space; Hc is the average collectable water volume in each first-type space;

[0110] In practical applications, the collectable water volume within each type I space is calculated as: Hc = Hp - Ha = 1.8 - 1.6 = 0.2 m³ 3 This means that, on average, 0.2 m³ of water can be collected per first-category space. 3 .

[0111] Step S8: Obtain the first proportion based on the fourth number of fourth-class spaces; Step S8 includes the following sub-steps:

[0112] Step S801: Obtain the fourth number of fourth-class spaces and mark them as historical spaces; here, the fourth-class spaces are obtained based on the three-dimensional structure model of the aquifer. Obtaining the fourth number of fourth-class spaces is to analyze the average volume of the three-dimensional structure model of the aquifer within the fourth-class spaces, so as to facilitate the subsequent calculation of the overall volume of the three-dimensional structure model of the aquifer.

[0113] Step S802: Obtain the volume of the three-dimensional aquifer structure model in the historical space and mark it as the historical volume. The historical volume can be obtained by any method, such as constructing the solid volume of the three-dimensional aquifer structure model, using the liquid filling method, and calculating the difference between the overall volume of the fourth type of space and the volume of liquid filling to obtain the historical volume.

[0114] Step S803: Obtain the total volume of the historical space and the total volume of the historical space;

[0115] Step S804: Calculate the ratio of the total historical volume to the total historical spatial volume, and mark it as the first ratio; the first ratio is the average proportion of the three-dimensional structural model of the aquifer in the fourth type of space;

[0116] In practical applications, for example, the first ratio is calculated to be 0.32.

[0117] Step S9 involves calculating the total collectable water volume of the aquifer based on the number of spaces of the first type, the number of spaces of the third type, the number of spaces of the fourth type, the collectable water volume within the spaces of the first type, and the first ratio. Step S9 includes the following sub-steps:

[0118] Step S901: Obtain the number of first-type spaces, marked as D1; ​​obtain the number of all third-type spaces, marked as D3; obtain the number of all fourth-type spaces, marked as D4; obtain the number of reference points within each second-type space, marked as D5;

[0119] Step S902, calculate the total collectable water volume of the aquifer as: Hz = D1 × Hc + (D3 + D4 × b1) / D5 × Hc; where Hz is the total collectable water volume of the aquifer, and b1 is the first ratio; here (D3 + D4 × b1) / D5 is equivalent to the number of spaces that convert the second type of space into the first type of space.

[0120] In practical applications, for example, if D1 is 29212, D3 is 62421, D4 is 82431, and D5 is 329600, the total collectable water volume of the aquifer can be calculated as: Hz = 29212 × 0.2 + (62421 + 82431 × 0.32) / 329600 × 0.2 = 5842 m³ 3 If the calculation result is rounded to the nearest integer, then the total collectable water volume of the aquifer in this area is 5842 m³. 3The water content in the three-dimensional structural model of the aquifer to be tested cannot exceed 5842m. 3 .

[0121] Example 2: This application also provides an electronic device, which may include: a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus. The memory stores computer-readable instructions, and the processor can call the instructions in the memory. When the computer-readable instructions are executed by the processor, the steps in the groundwater balance calculation simulation method based on three-dimensional geological data are performed to achieve the following functions: acquiring a three-dimensional structural model of the aquifer to be detected and marking it as the aquifer three-dimensional structural model; constructing an independent analysis space based on the aquifer three-dimensional structural model; dividing the aquifer three-dimensional structural model into a first type of space and a second type of space based on the independent analysis space; constructing reference points in the second type of space; acquiring the third type of space and the fourth type of space in the second type of space based on the reference points; obtaining the spatial average water content based on the water content in the first type of space of a third number of aquifers to be detected; obtaining the collectable water content in the first type of space based on the spatial average water content; obtaining a first proportion based on a fourth number of fourth type of spaces; calculating the overall collectable water content of the aquifer based on the number of first type of spaces, the number of third type of spaces, the number of fourth type of spaces, the collectable water content in the first type of space, and the first proportion.

[0122] Furthermore, when the logical instructions in the aforementioned memory can be implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0123] Example 3: This application also provides a computer program product, which includes a computer program stored on a computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the groundwater balance calculation and simulation method based on three-dimensional geological data provided by the above methods. The method includes: acquiring a three-dimensional structural model of the aquifer to be detected and marking it as an aquifer three-dimensional structural model; constructing an independent analysis space based on the aquifer three-dimensional structural model; dividing the aquifer three-dimensional structural model into a first type of space and a second type of space based on the independent analysis space; constructing reference points in the second type of space; acquiring a third type of space and a fourth type of space in the second type of space based on the reference points; acquiring the spatial average water content based on the water content in the first type of space of a third number of aquifers to be detected; acquiring the collectable water content in the first type of space based on the spatial average water content; acquiring a first proportion based on a fourth number of fourth type of spaces; and calculating the overall collectable water content of the aquifer based on the number of first type of spaces, the number of third type of spaces, the number of fourth type of spaces, the collectable water content in the first type of space, and the first proportion.

[0124] Example 4: This application also provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, it performs the steps of the above-described groundwater balance calculation and simulation method based on three-dimensional geological data to achieve the following functions: acquiring a three-dimensional structural model of the aquifer to be detected and labeling it as the aquifer three-dimensional structural model; constructing an independent analysis space based on the aquifer three-dimensional structural model; dividing the aquifer three-dimensional structural model into a first type of space and a second type of space based on the independent analysis space; constructing reference points within the second type of space; acquiring a third type of space and a fourth type of space within the second type of space based on the reference points; obtaining the spatial average water content based on the water content in the first type of space of a third number of aquifers to be detected; obtaining the collectable water content in the first type of space based on the spatial average water content; obtaining a first proportion based on a fourth number of fourth type of spaces; and calculating the overall collectable water content of the aquifer based on the number of first type of spaces, the number of third type of spaces, the number of fourth type of spaces, the collectable water content in the first type of space, and the first proportion.

[0125] Based on the above description of the embodiments, the embodiments of the present invention can be provided as methods, systems, or computer program products. Based on this understanding, the above technical solutions, in essence or in terms of their contribution to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or certain parts of the embodiments.

[0126] In the embodiments provided in this application, it should be understood that the disclosed system or method can be implemented in other ways. The embodiments described above are merely illustrative. For example, the division of modules or units is only a logical functional division, and there may be other division methods in actual implementation. Furthermore, multiple modules or units may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the coupling or direct coupling or communication connection shown or discussed may be through some communication interfaces. The indirect coupling or communication connection between systems, modules, and units may be electrical, mechanical, or other forms.

[0127] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application.

Claims

1. A groundwater balance calculation and simulation method based on three-dimensional geological data, characterized in that, Includes the following steps: Obtain the three-dimensional structural model of the aquifer to be detected and mark it as the three-dimensional structural model of the aquifer; An independent analysis space was constructed based on the three-dimensional structural model of the aquifer. Based on the independent analysis space, the three-dimensional structural model of the aquifer is divided into a first type of space and a second type of space; Construct spatial reference points within the second type of space; Based on reference points within the space, obtain the third and fourth types of spaces in the second type of space; The spatial average water content is obtained based on the water content of the first type of space in the third number of aquifers to be detected; The amount of collectable water in the first type of space is obtained based on the spatial average water content. The first proportion is obtained based on the fourth number of fourth-type spaces; The total collectable water volume of the aquifer is calculated based on the number of first-type spaces, the number of third-type spaces, the number of fourth-type spaces, the collectable water volume in the first-type spaces, and the first ratio. Constructing an independent analysis space based on a three-dimensional aquifer structural model includes the following sub-steps: establishing a three-dimensional coordinate system, labeled as the model's three-dimensional coordinate system; placing the aquifer three-dimensional structural model within the model's three-dimensional coordinate system; labeling the plane containing both the X and Y axes as the XY plane; labeling the plane parallel to the XY plane as the first plane; starting from the XY plane, drawing a first number of first planes with a first length interval and directions corresponding to the positive and negative directions of the Z-axis; labeling the plane containing both the X and Z axes as the XZ plane; labeling the plane parallel to the XZ plane as the second plane. Starting with the XZ plane, draw a first number of second planes with a first length interval between each second plane and the positive and negative directions of the Y-axis. Mark the plane containing both the Y-axis and Z-axis as the YZ plane. Mark the plane parallel to the YZ plane as the third plane. Starting with the YZ plane, draw a first number of third planes with a first length interval between each third plane and the positive and negative directions of the Z-axis. Divide the model's three-dimensional coordinate system into cubes with a first length side using all the drawn first, second, and third planes, marking them as independent analysis spaces. Dividing the aquifer 3D structure model into a first type of space and a second type of space based on the independent analysis space includes the following sub-steps: obtaining the independent analysis space that consists entirely of aquifer 3D structure models, and marking it as the first type of space; obtaining the independent analysis space that is divided into aquifer 3D structure models within the independent analysis space, and marking it as the second type of space; Obtaining the third and fourth types of spaces within the second type of space based on reference points within the space includes the following sub-steps: drawing a sphere with the reference point within the space as the center and the radius as the third length, and marking it as the sphere reference space; obtaining the sphere reference space within which all elements are aquifer 3D structural models, and marking it as the third type of space; obtaining the sphere reference space within which the interior is divided into aquifer 3D structural models, and marking it as the fourth type of space; Obtaining the first proportion based on the fourth number of fourth-type spaces includes the following sub-steps: obtaining the fourth number of fourth-type spaces and marking them as historical spaces; obtaining the volume of the three-dimensional structural model of the aquifer in the historical space and marking it as the historical volume; obtaining the sum of the volumes of the historical spaces and the sum of the historical volumes. Calculate the ratio of the total historical volume to the total historical spatial volume, and mark it as the first ratio; The overall collectable water volume of the aquifer is calculated based on the number of first-type spaces, the number of third-type spaces, the number of fourth-type spaces, the collectable water volume in the first-type spaces, and the first ratio. This includes the following sub-steps: obtaining the number of first-type spaces, marked as D1; ​​obtaining the total number of third-type spaces, marked as D3; obtaining the total number of fourth-type spaces, marked as D4. Obtain the number of reference points within each second-class space and label it as D5; The total collectable water volume of the aquifer is calculated as: Hz = D1 × Hc + (D3 + D4 × b1) / D5 × Hc; where Hz is the total collectable water volume of the aquifer and b1 is the first proportion.

2. The groundwater balance calculation and simulation method based on three-dimensional geological data according to claim 1, characterized in that, Constructing an intra-space reference point within the second type of space includes the following sub-steps: The top surface of the second type of space is marked as the second top surface, and the surface parallel to the second top surface is marked as the parallel top surface; Starting with the parallel top surface, with the interval between each surface being the second length and the direction being the bottom surface of the second type of space, draw the second number of parallel top surfaces; The side of the second type of space is labeled as the second side, and the surface parallel to the second side is labeled as the parallel side. Starting with the second side, with the interval between each parallel side being the second length and the direction being another side of the second type of space, draw the second number of parallel sides; The front of the second type of space is marked as the second front, and the face parallel to the second front is marked as the parallel front; Starting with the second front, the interval between each parallel front is the second length, and the direction is the back of the second type of space, and a second number of parallel fronts are drawn; Obtain the intersection point of the three surfaces parallel to the top surface, parallel to the side surface, and parallel to the front surface in the second type of space, and mark it as the reference point in the space.

3. The groundwater balance calculation and simulation method based on three-dimensional geological data according to claim 2, characterized in that, Obtaining the spatial average water content based on the water content of the first type of space in the third number of aquifers to be detected includes the following sub-steps: Obtain the water content of the first type of space in the third number of aquifers to be tested, and mark it as the experimental water content; Obtain the range of experimental water content, divide the range of experimental water content into k equally spaced intervals, and mark them as water content division intervals; Count the frequency of each moisture content division interval and label it as the moisture content division frequency; A histogram was plotted with the experimental moisture content as the X-axis, the frequency of moisture content division as the Y-axis, and the intervals of moisture content division as the histogram intervals. This histogram was then labeled as the moisture content histogram.

4. The groundwater balance calculation and simulation method based on three-dimensional geological data according to claim 3, characterized in that, Obtaining the spatial average water content based on the water content of the first type of space in the third number of aquifers to be detected also includes the following sub-steps: The water content frequency threshold is calculated as: Psh = u × g3 ÷ k; where Psh is the water content frequency threshold, u is a value between 0 and 1, and g3 is the third quantity; Water content less than or equal to the water content frequency threshold is classified as an abnormal frequency. Delete the water content division intervals with the leftmost and rightmost abnormal division frequencies in the water content histogram, and obtain the minimum and maximum values ​​of the water content division intervals in the water content histogram after deletion, and mark them as the first water content threshold and the second water content threshold. The mean of experimental water content that is greater than or equal to the first water content threshold and less than or equal to the second water content threshold is calculated and marked as the spatial average water content.

5. The groundwater balance calculation and simulation method based on three-dimensional geological data according to claim 4, characterized in that, Obtaining the collectable water volume of the first type of space based on the spatial average water content includes the following sub-steps: Obtain the lowest safe water content within the first type of space, denoted as Ha; The collectable water volume in each first-type space is calculated as: Hc = Hp - Ha; where Hc is the collectable water volume in each first-type space, and Hp is the average water content of the space.

Images