A method, system, computer equipment and medium for calculating the impact range of groundwater circulation well pollution remediation

Abstract

The present invention provides a method, system, computer equipment and medium for calculating the impact range of groundwater circulation well pollution remediation, which belongs to the field of groundwater pollution remediation technology. The method includes the following steps: setting a circulation well; measuring the groundwater variables of all position point sets T1 before and after the circulation well is operated; setting a threshold range for the difference in groundwater variables; calculating the difference in groundwater variables at all positions before and after the circulation well is operated; constructing a three-dimensional point set T2 of the groundwater circulation well impact range; performing Delaunay triangulation on the three-dimensional point set T2; and calculating the impact volume of the Delaunay triangulated circulation well impact range. The present invention adopts the above-mentioned method, system, computer equipment and medium for calculating the impact range of groundwater circulation well pollution remediation to solve the defect that the impact range of the circulation well in heterogeneous aquifers cannot be accurately calculated due to the complex hydrogeological conditions.

Description

Technical Field

[0001] The present invention relates to the technical field of groundwater pollution remediation, and in particular to a method, system, computer equipment and medium for calculating the impact range of groundwater circulation well pollution remediation. Background Art

[0002] Carrying out effective groundwater pollution remediation work is a hot and difficult issue in current research. Among groundwater pollution remediation technologies, groundwater circulation well (GCW) technology has a high efficiency in groundwater pollution remediation, but this technology currently faces the problem of being unable to accurately calculate the impact range. At present, most of the calculations of the impact range of groundwater circulation wells are based on the perspective of impact radius, and most of the research only uses analytical solutions to calculate the impact radius in homogeneous aquifers, while aquifers in nature are often heterogeneous and have hydraulic gradients. Some hydrogeological parameters such as permeability coefficient and porosity will have a huge impact on the impact range of groundwater circulation wells. The currently commonly used calculation method of the impact radius of circulation wells cannot obtain the actual impact range distribution, resulting in poor accuracy and large errors in the calculation of the impact range, which affects the remediation of groundwater pollution. Summary of the Invention

[0003] The purpose of the present invention is to provide a method, system, computer equipment and medium for calculating the influence range of groundwater circulation well pollution remediation, which combines the calculation of the influence range of the circulation well with the Delaunay triangulation method, solves the defect of large errors in the calculation of the influence range of the circulation well in heterogeneous aquifers, and realizes accurate calculation of the influence range of the circulation well under complex (heterogeneous and natural gradient) conditions, laying a solid foundation for research in the field of groundwater pollution remediation technology.

[0004] To achieve the above objectives, the present invention provides a method for calculating the impact range of groundwater circulation well pollution remediation, comprising the following steps:

[0005] Step S1: setting a circulation well in the contaminated area;

[0006] Step S2, measuring the groundwater variables of all position point sets T1 before the circulation well is operated;

[0007] Step S3, measuring the groundwater variables of all location point sets T1 after the circulation well is operated;

[0008] Step S4, setting a groundwater variable difference threshold range;

[0009] Step S5, calculating the groundwater variable differences at all locations before and after the circulation well is operated;

[0010] Step S6: extracting points that meet the variable difference threshold range from all point sets T1 after calculating the groundwater variable difference to form a three-dimensional point set T2 of the groundwater circulation well influence range;

[0011] Step S7, dividing the three-dimensional point set T2 of the groundwater circulation well influence range into two parts containing positive values ​​and negative values ​​according to the positive and negative values ​​of the groundwater variable difference;

[0012] Step S8, performing Delaunay triangulation on the point set containing positive values ​​and the point set containing negative values ​​in the three-dimensional point set T2 of the influence range of the groundwater circulation well respectively;

[0013] Step S9: Calculate the influence volume of the circulation well influence range after completing the Delaunay triangulation.

[0014] Preferably, in step S2 and step S3, the groundwater variables include groundwater head and pollutant concentration.

[0015] Preferably, in step S6, the three-dimensional point set T2 of the influence range of the groundwater circulation well is as follows:

[0016] T2={p1,p2,p3,...,p n};

[0017] Among them, p i =(x i ,y i ,z i ) is a point in the three-dimensional point set T2 within the influence range of the groundwater circulation well, x i For point p i x coordinate, y i For point p i The y coordinate, z i For point p i The z coordinate of

[0018] Preferably, in step S8, Delaunay triangulation is performed on the point set containing positive values ​​and the point set containing negative values ​​in the three-dimensional point set T2 of the influence range of the groundwater circulation well, respectively. The specific operation is:

[0019] Step S81: Initially create a tetrahedron set L and an empty point set P;

[0020] Step S82: add the first point p1 in the three-dimensional point set T2 of the groundwater circulation well influence range to the set P, initialize a tetrahedron containing the point p1, and add the tetrahedron containing the point p1 to the tetrahedron set L;

[0021] Step S83: For point p in the three-dimensional point set T2 within the influence range of the groundwater circulation well,i , find the file containing p i Tetrahedron

[0022] Step S84: Determine the point p in the three-dimensional point set T2 within the influence range of the groundwater circulation well. i With tetrahedron The relationship between the outside ball:

[0023] Step S841: Determine the tetrahedron The center of the circumscribed ball and the circumscribed sphere radius

[0024] Step S842: Check the point p adjacent to p1, p2, p3, and p4 in the three-dimensional point set T2 within the influence range of the groundwater circulation well. i Is it in the tetrahedron Inside the outside ball:

[0025] When point p i In the tetrahedron The outside of the circumscribed ball, that is, When , the tetrahedron satisfies the conditions of Delaunay triangulation;

[0026] When point p i In the tetrahedron The inside of the circumscribed ball satisfies When tetrahedron Does not meet the conditions of Delaunay triangulation, for tetrahedron Refactoring

[0027] Step S85: for each tetrahedron Adjacent tetrahedrons When the tetrahedron With p i When the circumscribed spheres intersect, the tetrahedron Remove, forming a p i and Each edge of is a new tetrahedron formed by the edges, that is, a tetrahedron that satisfies the Delaunay triangulation condition;

[0028] Step S86: All points p in the three-dimensional point set T2 affected by the groundwater circulation well i The newly formed tetrahedron that meets the Delaunay condition is added to the tetrahedron set L;

[0029] Step S87: When all points p in the three-dimensional point set T2 within the influence range of the groundwater circulation well are iAfter all the additions are completed, the Delaunay triangulation of all points in the influence range of the circulation well is completed to form a three-dimensional circulation well influence range.

[0030] Preferably, in step S9, the influence volume of the circulation well influence range of the Delaunay triangulation is calculated, and the specific operation is:

[0031] Step S91: After completing the Delaunay triangulation of all points in the influence range of the circulation well, a tetrahedron set is obtained. Each tetrahedron The vertex p is the vertex of the tetrahedron composed of every four points that conform to the Delaunay triangulation rule in the three-dimensional point set T2 within the influence range of the groundwater circulation well. i1 ,p i2 ,p i3 ,p i4 ;

[0032] Step S92: Calculate the three edge vectors v associated with the vertices of the tetrahedron:

[0033] v i1 =p i2 -p i1 , v i2 =p i3 -p i1 , v i3 =p i4 -p i1 ;

[0034] Step S93, calculate the p i1 ,p i2 ,p i3 ,p i4 The directed volume v of the parallelepiped i1 ·(v i2 ×v i3 ):

[0035]

[0036] Where e1, e2, and e3 are unit basis vectors; · represents vector dot product;

[0037] Step S94: Calculate the volume V of the influence range of the three-dimensional circulation well:

[0038]

[0039] Where m is the number of tetrahedrons that satisfy the Delaunay condition; i is the index of the tetrahedron that satisfies the Delaunay condition.

[0040] The present invention also provides a groundwater circulation well pollution remediation impact range calculation system, comprising:

[0041] The data processing module is used to set the groundwater variable difference threshold range, calculate the groundwater variable difference at all locations before and after the circulation well is operated, and construct the three-dimensional point set T2 of the groundwater circulation well influence range;

[0042] Delaunay triangulation module, used to perform Delaunay triangulation on the three-dimensional point set T2 within the influence range of the groundwater circulation well;

[0043] The volume calculation module is used to calculate the influence volume of the circulation well influence range that completes the Delaunay triangulation.

[0044] The present invention also provides a computer device, comprising: a memory and a processor; the memory stores a computer program, and when the processor executes the computer program, it implements the steps of the above-mentioned method for calculating the impact range of groundwater circulation well pollution remediation.

[0045] The present invention also provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps of the above-mentioned method for calculating the impact range of groundwater circulation well pollution remediation are implemented.

[0046] Therefore, the present invention adopts the above-mentioned method, system, computer equipment and medium for calculating the impact range of groundwater circulation well pollution remediation, and the beneficial technical effects are as follows:

[0047] (1) The present invention is not limited to the evaluation of the influence range of groundwater circulation wells in homogeneous aquifers, but is applicable to all aquifer conditions. It does not only calculate the maximum influence radius of the circulation well, but also collects data points by measuring groundwater variables (groundwater head, pollutant concentration) at different positions before and after the circulation well is operated in a natural (heterogeneous) aquifer, and performs Delaunay triangulation on these point sets. The actual shape of the influence range of the groundwater circulation well is constructed through Delaunay triangulation.

[0048] (2) The present invention can accurately calculate the volume of the actual shape of the influence range of the groundwater circulation well through Delaunay triangulation. The calculated volume is the influence volume of the groundwater circulation well. The influence range calculated using this method will not be affected by the limitations of the analytical solution calculation. Moreover, the present invention can predict and calculate the influence range of the groundwater circulation well in the actual aquifer. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 This is a schematic diagram of the circulation well structure;

[0050] Figure 2Schematic diagram of the aquifer model grid division;

[0051] Figure 3 is the distribution diagram of the hydraulic conductivity of the heterogeneous aquifer a;

[0052] Figure 4 The pollution concentration distribution diagram of the polluted heterogeneous aquifer a after the circulation well is operated;

[0053] Figure 5 Schematic diagram of the three-dimensional point set T2 in the influence range of the groundwater circulation well in the heterogeneous aquifer a;

[0054] Figure 6 3D schematic diagram of the influence range of the circulation well in heterogeneous aquifer a;

[0055] Figure 7 is the distribution diagram of the hydraulic conductivity of the heterogeneous aquifer b;

[0056] Figure 8 The groundwater head distribution diagram of the contaminated heterogeneous aquifer b after the circulation well is operated;

[0057] Figure 9 Schematic diagram of the three-dimensional point set T2 in the influence range of the groundwater circulation well in the heterogeneous aquifer b;

[0058] Figure 10 3D schematic diagram of the influence range of the circulation well in the heterogeneous aquifer b.

[0059] Reference numerals

[0060] 1. Pumping screen section; 2. Injection screen section; 3. Actual impact range; 4. Predicted impact range. DETAILED DESCRIPTION

[0061] The technical solution of the present invention is further described below with reference to the accompanying drawings and embodiments.

[0062] Unless otherwise defined, technical or scientific terms used in the present invention shall have the same meaning as commonly understood by one of ordinary skill in the art to which the present invention belongs.

[0063] Example 1

[0064] Circulation well configuration.

[0065] A method for calculating the impact range of groundwater circulation well pollution remediation is to arrange the circulation wells as follows: Figure 1 As shown, the circulation well includes a water injection screen section and a water pumping screen section, wherein the water injection screen section 2 is set at the lower end of the circulation well, and the height of the water injection screen section 2 is 1m, and the water pumping screen section 1 is set at the upper end of the circulation well, and the height of the water pumping screen section 1 is 1m. The water injection volume of the water injection screen section 2 is 50m 3 / d, the pumping capacity of the pumping screen section 1 is 50m3 The actual impact range 3 of the pollution remediation after the circulation well operation can be converted into the predicted impact range 4 using the Delaunay triangulation method.

[0066] Heterogeneous aquifer model configuration and influence range calculation method.

[0067] This example primarily uses the Monte Carlo method to set the permeability coefficient (K) of the study area's aquifers to simulate different heterogeneous aquifers. VSAFT3 (variably saturated aquifer flow and transport in 3D) is also used to simulate groundwater flow and solute transport in these heterogeneous aquifers. The parameters for generating the permeability coefficient (K) of random heterogeneous aquifers using the Monte Carlo method are shown in Table 1 below.

[0068] Table 1 Monte Carlo parameters for the permeability coefficient K of random heterogeneous aquifers

[0069]

[0070] The heterogeneous aquifer model is set as a 50m×50m×20m rectangular aquifer, and the model is divided into 20 layers, each layer is divided into 50×50 cells, and the length, width and height of each cell are 1m. The plane range of the aquifer is 50m×50m, the top elevation of the aquifer is 20m, and the bottom elevation is 0m. It is evenly divided into 50,000 1m×1m×1m cubes. The boundary at x=0m in the x-direction of the aquifer is a constant head boundary of 20.25m, and the boundary at x=50m in the x-direction is a constant head boundary of 20m, forming a natural water flow from the boundary at x=0m of the aquifer to the boundary at x=50m of the aquifer. The pollutant in the groundwater is oxytetracycline. The initial pollutant concentration in the model is set to 1mg / L, and the surrounding concentration boundaries are all constant concentration boundaries of 1mg / L, with a diffusion coefficient of 0.001m 2 / d. Aquifer model grid division diagram is shown in Figure 2 .

[0071] After completing the above model configuration, the steps to calculate the influence range of the circulation well are as follows:

[0072] ①Construct a three-dimensional point set T2 within the influence range of the groundwater circulation well:

[0073] T2={p1,p2,p3,...,p n};

[0074] Among them, p i =(x i ,y i ,z i) is a point in the three-dimensional point set T2 within the influence range of the groundwater circulation well; x i For point p i x coordinate; y i For point p i The y coordinate of i For point p i The z coordinate of

[0075] ② According to the positive and negative values ​​of the groundwater variable difference, the three-dimensional point set T2 of the groundwater circulation well influence range is divided into two parts containing positive values ​​and negative values.

[0076] ③ Initially create a tetrahedron set L and an empty point set P.

[0077] ④ Add the first point p1 in the three-dimensional point set T2 of the groundwater circulation well influence range to the set P, initialize a tetrahedron containing point p1, and add the tetrahedron containing point p1 to the tetrahedron set L;

[0078] ⑤ For point p in the three-dimensional point set T2 within the influence range of the groundwater circulation well i , find the file containing p i Tetrahedron

[0079] ⑥ Determine the point p in the three-dimensional point set T2 within the influence range of the groundwater circulation well i With tetrahedron The relationship between the outside ball:

[0080] 1) Determine the tetrahedron The center of the circumscribed ball and the circumscribed sphere radius

[0081] 2) Check the point p adjacent to p1, p2, p3, and p4 in the three-dimensional point set T2 within the influence range of the groundwater circulation well. i Is it in the tetrahedron Inside the outside ball:

[0082] When point p i In the tetrahedron The outside of the circumscribed ball, that is, When , the tetrahedron satisfies the Delaunay condition;

[0083] When point p i In the tetrahedron The inside of the circumscribed ball satisfies When tetrahedron Does not satisfy the Delaunay condition for tetrahedrons Refactoring

[0084] ⑦Reconstruct tetrahedron: For each tetrahedron Adjacent tetrahedrons When the tetrahedron With p i When the circumscribed spheres intersect, the tetrahedron Remove, forming a p i and Each edge of is a new tetrahedron formed by the edges, which is a tetrahedron that satisfies the Delaunay condition;

[0085] ⑧ Update the tetrahedron set L: All points p in the three-dimensional point set T2 affected by the groundwater circulation well i The newly formed tetrahedron that meets the Delaunay condition is added to the tetrahedron set L;

[0086] ⑨ When all points p in the three-dimensional point set T2 affected by the groundwater circulation well i After all the additions are completed, the Delaunay triangulation of all points in the influence range of the circulation well is completed to form a three-dimensional circulation well influence range.

[0087] ⑩ After completing the Delaunay triangulation of all points within the influence range of the circulation well, a set of tetrahedrons is obtained. Each tetrahedron The vertex p is the vertex of the tetrahedron composed of every four points that conform to the Delaunay triangulation rule in the three-dimensional point set T2 within the influence range of the groundwater circulation well. i1 ,p i2 ,p i3 ,p i4 .

[0088] Compute the three edge vectors v associated with the vertices of the tetrahedron:

[0089] v i1 =p i2 -p i1 , v i2 =p i3 -p i1 , v i3 =p i4 -p i1 ;

[0090] Calculated by p i1 ,p i2 ,p i3 ,p i4 The directed volume v of the parallelepiped i1 ·(v i2 ×v i3 ):

[0091]

[0092] Where e1, e2, and e3 are unit basis vectors; · represents vector dot product;

[0093] Calculate the volume V of the influence range of the three-dimensional circulation well:

[0094]

[0095] Where m is the number of tetrahedrons that satisfy the Delaunay condition; i is the index of the tetrahedron that satisfies the Delaunay condition.

[0096] Example 2

[0097] Concentration-based method for calculating the remediation scope of groundwater circulation well contamination.

[0098] The influence range of the circulation well is calculated by using the pollution concentration variable for the polluted heterogeneous aquifer a. The distribution of the permeability coefficient of the heterogeneous aquifer a is shown in Figure 3 .

[0099] By measuring the pollution concentration before and after the operation of the circulation well in the heterogeneous aquifer a, the pollution remediation results of this heterogeneous aquifer under the operation of the groundwater circulation well can be obtained. The pollution concentration distribution of the polluted heterogeneous aquifer a after the operation of the circulation well is shown in Figure 2. Figure 4 .

[0100] The three-dimensional point set T2 of the circulation well influence range is formed by the position points where the pollution concentration changes within the threshold range of 1%. The schematic diagram of the three-dimensional point set T2 is shown in Figure 5 (A 1% pollution concentration change threshold range is used in this example. In other practical situations, other threshold ranges may be selected based on the actual aquifer).

[0101] Perform Delaunay triangulation on the three-dimensional point set T2 of the circulation well influence range that meets the threshold, and obtain a three-dimensional schematic diagram of the circulation well influence range of the heterogeneous aquifer a. Figure 6 .

[0102] After Delaunay triangulation of the influence range of the circulation well in the heterogeneous aquifer a, the influence volume of the circulation well can be calculated. It can be found that the total influence volume of the circulation well in the heterogeneous aquifer a is 5382m 3 .

[0103] Example 3

[0104] A method for calculating the remediation scope of groundwater circulation well contamination based on groundwater head.

[0105] The influence range of the circulation well is calculated by using the groundwater head variable for the contaminated heterogeneous aquifer b. The distribution of the permeability coefficient of the heterogeneous aquifer b is shown in Figure 7 .

[0106] The heterogeneous aquifer b has a natural hydraulic gradient in the positive x direction throughout the aquifer area. By measuring the groundwater head before and after the circulation well in the heterogeneous aquifer b, the water head distribution of the heterogeneous aquifer b after the circulation well is operated can be obtained. Figure 8 .

[0107] The three-dimensional point set T2 of the circulation well influence range is formed by the position points where the groundwater head change is within the threshold range of 1%. The schematic diagram of the three-dimensional point set T2 of this embodiment is shown in FIG. Figure 9 (In this embodiment, a groundwater head change threshold range of 1% is used. In other practical situations, other threshold ranges may be selected according to the actual aquifer).

[0108] Perform Delaunay triangulation on the three-dimensional point set T2 of the circulation well influence range that meets the threshold, and obtain a three-dimensional schematic diagram of the circulation well influence range of the heterogeneous aquifer b as shown in Figure 10 .

[0109] After Delaunay triangulation of the influence range of the circulation well in the heterogeneous aquifer b, the influence volume of the circulation well can be calculated. It can be obtained that the total influence volume of the circulation well in the heterogeneous aquifer b until it stabilizes is 3843m 3 .

[0110] Example 4

[0111] A groundwater circulation well pollution remediation impact range calculation system, comprising:

[0112] The data processing module is used to set the groundwater variable difference threshold range, calculate the groundwater variable difference at all locations before and after the circulation well is operated, and construct the three-dimensional point set T2 of the groundwater circulation well influence range;

[0113] Delaunay triangulation module, used to perform Delaunay triangulation on the three-dimensional point set T2 within the influence range of the groundwater circulation well;

[0114] The volume calculation module is used to calculate the influence volume of the circulation well influence range that completes the Delaunay triangulation.

[0115] If the above functions are implemented in the form of 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 the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0116] The logic and / or steps represented in the flowcharts or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing the logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (e.g., a computer-based system, a system including a processor, or other system that can fetch and execute instructions from an instruction execution system, apparatus, or device). For purposes of this specification, a "computer-readable medium" can be any device that can contain, store, communicate, propagate, or transport a program for use by, or in conjunction with, an instruction execution system, apparatus, or device.

[0117] More specific examples (a non-exhaustive list) of computer-readable media include the following: a portable computer disk cartridge (magnetic device), random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), fiber optic devices, and portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program can be obtained electronically, for example, by optically scanning the paper or other medium, and then editing, interpreting, or processing in another suitable manner as necessary, and then storing it in a computer memory.

[0118] It is worth noting that the contents not elaborated in detail in the present invention are all prior art and are well known to those skilled in the art.

[0119] Therefore, the present invention adopts the above-mentioned method, system, computer equipment and medium for calculating the influence range of groundwater circulation well pollution remediation, combines the calculation of the influence range of the circulation well with the Delaunay triangulation method, and uses variables such as the pollution concentration and groundwater head of the aquifer where the circulation well is located to calculate the influence range of the circulation well, thereby solving the current technical bottleneck of large errors in the prediction calculation of the influence range of the circulation well in heterogeneous aquifers, and realizing the accurate calculation of the influence range of the circulation well under complex conditions, laying the foundation for research in the field of groundwater pollution remediation technology, especially in the field of groundwater circulation wells.

[0120] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit the same. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that they can still modify or replace the technical solutions of the present invention with equivalents, and these modifications or equivalent replacements cannot cause the modified technical solutions to deviate from the spirit and scope of the technical solutions of the present invention.

Claims

1. A method for calculating the impact range of groundwater circulation well pollution remediation, characterized in that: The following steps are involved: Step S1: setting a circulation well in the contaminated area; Step S2, measuring the groundwater variables of all position point sets T1 before the circulation well is operated; Step S3, measuring the groundwater variables of all location point sets T1 after the circulation well is operated; Step S4, setting a groundwater variable difference threshold range; Step S5, calculating the groundwater variable differences at all locations before and after the circulation well is operated; Step S6: extracting points that meet the variable difference threshold range from all point sets T1 after calculating the groundwater variable difference to form a three-dimensional point set T2 of the groundwater circulation well influence range; Step S7, dividing the three-dimensional point set T2 of the groundwater circulation well influence range into two parts containing positive values ​​and negative values ​​according to the positive and negative values ​​of the groundwater variable difference; Step S8, performing Delaunay triangulation on the point set containing positive values ​​and the point set containing negative values ​​in the three-dimensional point set T2 of the influence range of the groundwater circulation well respectively; Step S9, calculating the influence volume of the circulation well influence range after completing the Delaunay triangulation; In step S6, the three-dimensional point set T2 of the influence range of the groundwater circulation well is as follows: T2={p1,p2,p3,...,p n }; Among them, p i =(x i ,y i ,z i ) is a point in the three-dimensional point set T2 within the influence range of the groundwater circulation well, x i For point p i x coordinate, y i For point p i The y coordinate, z i For point p i The z coordinate of In step S9, the influence volume of the circulation well influence range after Delaunay triangulation is calculated. The specific operation is: Step S91: After completing the Delaunay triangulation of all points in the influence range of the circulation well, a tetrahedron set L = (l1, l2, l3, ..., l m ), where each tetrahedron l i The vertex p is the vertex of the tetrahedron composed of every four points that conform to the Delaunay triangulation rule in the three-dimensional point set T2 within the influence range of the groundwater circulation well. i1 ,p i2 ,p i3 ,p i4 ; Step S92: Calculate the three edge vectors v associated with the vertices of the tetrahedron: v i1 =p i2 -p i1 ,v i2 =p i3 -p i1 ,v i3 =p i4 -p i1 ; Step S93, calculate the p i1 ,p i2 ,p i3 ,p i4 The directed volume v of the parallelepiped i1 ·(v i2 ×v i3 ): Where e1, e2, and e3 are unit basis vectors; · represents vector dot product; Step S94: Calculate the volume V of the influence range of the three-dimensional circulation well: Where m is the number of tetrahedrons that satisfy the Delaunay condition; i is the index of the tetrahedron that satisfies the Delaunay condition.

2. The method for calculating the impact range of groundwater circulation well pollution remediation according to claim 1 is characterized in that: In step S2 and step S3, groundwater variables include groundwater head and pollutant concentration.

3. The method for calculating the impact range of groundwater circulation well pollution remediation according to claim 2 is characterized in that: In step S8, Delaunay triangulation is performed on the point set containing positive values ​​and the point set containing negative values ​​in the three-dimensional point set T2 of the influence range of the groundwater circulation well. The specific operation is as follows: Step S81: Initially create a tetrahedron set L and an empty point set P; Step S82: add the first point p1 in the three-dimensional point set T2 of the groundwater circulation well influence range to the set P, initialize a tetrahedron containing the point p1, and add the tetrahedron containing the point p1 to the tetrahedron set L; Step S83: For point p in the three-dimensional point set T2 within the influence range of the groundwater circulation well, i , find the file containing p i The tetrahedron l∈L; Step S84: Determine the point p in the three-dimensional point set T2 within the influence range of the groundwater circulation well. i Relationship with the circumscribed sphere of tetrahedron l: Step S841, determine the center C(l) and radius r(l) of the circumscribed sphere of tetrahedron l; Step S842: Check the point p adjacent to p1, p2, p3, and p4 in the three-dimensional point set T2 within the influence range of the groundwater circulation well. i Is it inside the circumsphere of tetrahedron l? When point p i Outside the circumscribed sphere of tetrahedron l, that is, satisfying ||p i When -C(l)||>r(l), the tetrahedron satisfies the conditions for Delaunay triangulation; When point p i Inside the circumscribed sphere of tetrahedron l, that is, satisfying ||p i When -C(l)||<r(l), tetrahedron l does not meet the conditions of Delaunay triangulation, and tetrahedron l is reconstructed; Step S85: For each tetrahedron l' adjacent to tetrahedron l, when tetrahedron l' is adjacent to p i When the circumscribed spheres intersect, remove the tetrahedron l' to form a i The new tetrahedron formed by each edge of l′ is a tetrahedron that satisfies the Delaunay triangulation condition; Step S86: All points p in the three-dimensional point set T2 affected by the groundwater circulation well i The newly formed tetrahedron that meets the Delaunay condition is added to the tetrahedron set L; Step S87: When all points p in the three-dimensional point set T2 within the influence range of the groundwater circulation well are i After all the additions are completed, the Delaunay triangulation of all points in the influence range of the circulation well is completed to form a three-dimensional circulation well influence range.

4. A groundwater circulation well pollution remediation impact range calculation system, characterized in that: The method for calculating the impact range of groundwater circulation well pollution remediation according to any one of claims 1 to 3 comprises: The data processing module is used to set the groundwater variable difference threshold range, calculate the groundwater variable difference at all locations before and after the circulation well is operated, and construct the three-dimensional point set T2 of the groundwater circulation well influence range; Delaunay triangulation module, used to perform Delaunay triangulation on the three-dimensional point set T2 within the influence range of the groundwater circulation well; The volume calculation module is used to calculate the influence volume of the circulation well influence range that completes the Delaunay triangulation.

5. A computer device comprising: memory and processor; The memory stores a computer program, characterized in that when the processor executes the computer program, the steps of the method for calculating the impact range of groundwater circulation well pollution remediation described in any one of claims 1 to 3 are implemented.

6. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for calculating the impact range of groundwater circulation well pollution remediation described in any one of claims 1 to 3 are implemented.

Images