Underground reservoir defluorination optimization method, device and equipment
By establishing a three-dimensional underground reservoir model and optimizing the dosing points and dosage of the defluorinating agent, the problem of poor defluorination effect in underground reservoirs of coal mines was solved, and a highly efficient defluorination effect in underground reservoirs was achieved.
Patent Information
- Application Number
- CN202410514750.4
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-04-26
- Publication Date
- 2025-10-28
AI Technical Summary
Existing technologies are ineffective in defluorination in underground coal mine reservoirs, mainly due to insufficient or inadequate dosing of defluorinating agents, resulting in insufficient residence time of the defluorinating agents and failing to meet the defluorination requirements of underground reservoirs.
By establishing a three-dimensional groundwater reservoir model and setting boundary conditions, including the defluoridation agent dosage, dosing point and dosing speed, simulation software is used to optimize the defluoridation agent addition. The model is adjusted according to the current removal rate, and the dosing point and dosing amount are optimized to improve the defluoridation effect.
It has achieved efficient removal of fluoride pollutants from underground reservoirs, optimized the dosing process of defluorinating agents, and improved the defluorination effect and utilization rate of defluorinating agents.
Smart Images

Figure CN120841602A_ABST
Abstract
Description
Technical Field
[0001] This disclosure relates to the field of underground reservoir wastewater treatment, and in particular to the field of underground reservoir defluoridation optimization technology. Background Technology
[0002] With the large-scale mining of coal, high concentrations of fluoride pollutants have been detected in mine water in many mining areas, even far exceeding the emission standards (F- < 1 mg / L) in Class III of the "GB3838-2002 Surface Water Environmental Quality Standard". However, field experiments in underground reservoirs face challenges such as large experimental spaces, limited open areas, difficulties in adding chemicals at the working face, and the lack of agitation or other power facilities in the reservoir. These issues can easily lead to inadequate dosing of the defluoridating agent, resulting in insufficient residence time and insufficient dosage, thus affecting the defluoridation effect of the underground reservoir. Therefore, optimizing the defluoridation effect of underground reservoirs has become an urgent problem to be solved. Summary of the Invention
[0003] This disclosure provides an optimized method, apparatus, equipment, and storage medium for defluoridation of underground reservoirs.
[0004] According to a first aspect of this disclosure, an optimized method for defluoridation in underground water reservoirs is provided. The method includes:
[0005] Obtain parameters of the underground water reservoir in the mining area;
[0006] Based on the parameters of the underground reservoir, a three-dimensional underground reservoir model is established;
[0007] Boundary conditions are set for the three-dimensional underground reservoir model, wherein the boundary conditions include: water flow parameter boundaries and defluoridator parameter boundaries, wherein the defluoridator parameter boundaries include: defluoridator dosage boundaries, defluoridator dosing points in the underground reservoir, and defluoridator dosing rate boundaries, and the defluoridator is used to remove fluoride pollutants in the underground reservoir;
[0008] Run the three-dimensional underground reservoir model to obtain the current removal rate of the fluoride pollutant;
[0009] Based on the current removal rate, determine whether to adjust the three-dimensional underground reservoir model.
[0010] In addition to the aspects and any possible implementations described above, a further implementation is provided in which, before running the three-dimensional underground reservoir model, the method further includes:
[0011] In the three-dimensional underground reservoir model, the concentration of fluoride pollutants at the inlet of the underground reservoir is set as a variable;
[0012] In the three-dimensional underground reservoir model, the concentration of fluoride pollutants at the outlet of the underground reservoir is set as a monitoring variable;
[0013] In the three-dimensional underground reservoir model, the first physical property parameter of the collapsed rock mass of the underground reservoir is assigned, the second physical property parameter of the grouting filling material of the collapsed rock mass is assigned, and the particle size of the defluorinating agent is set in the three-dimensional underground reservoir model;
[0014] In the three-dimensional underground reservoir model, the collapsed rock mass and sidewalls of the underground reservoir are set as solid, non-slip boundaries.
[0015] In addition to the aspects and any possible implementations described above, a further implementation is provided in which the water flow parameter boundaries include: the flow rate boundaries of the inlet and outlet of the underground reservoir or the pressure head boundaries of the inlet and outlet of the underground reservoir.
[0016] The first physical property parameters include: the volume of the collapsed rock mass and the porosity of the collapsed rock mass;
[0017] The second physical property parameters include: porosity, permeability coefficient, water yield, water storage rate, viscosity, flow rate, and strength of the grouting filling material.
[0018] In addition to the aspects and any possible implementations described above, a further implementation is provided, wherein determining whether to adjust the three-dimensional underground reservoir model based on the current removal rate includes:
[0019] Obtain the target removal rate of fluoride pollutants from the three-dimensional underground reservoir model;
[0020] Compare the current removal rate with the target removal rate;
[0021] If the current removal rate is less than the target removal rate, then the boundary conditions of the three-dimensional underground reservoir model are adjusted.
[0022] In addition to the aspects and any possible implementations described above, a further implementation is provided, wherein the method further includes:
[0023] Run the three-dimensional underground reservoir model to obtain the flow field of the underground reservoir and the dispersion of the defluoridating agent in the underground reservoir;
[0024] Based on the flow field conditions of the underground reservoir and / or the dispersion of the defluorinating agent in the underground reservoir, it is determined whether to adjust the defluorinating agent dosing point in the underground reservoir.
[0025] In addition to the aspects and any possible implementations described above, a further implementation is provided, wherein determining whether to adjust the defluoridator dosing point in the underground reservoir by utilizing the flow field conditions of the underground reservoir and / or the reagent dispersion of the defluoridator in the underground reservoir includes:
[0026] If the flow field in the first preset area of the underground reservoir is turbulent and the Reynolds number in the turbulent area is higher than a preset threshold, then the first preset area will be added as a defluoridation agent dosing point in the underground reservoir; or
[0027] If the flow field in the second preset area of the underground reservoir is a turbulent zone and the dispersion of the defluorinating agent in the preset area is higher than the preset agent dispersion, then the second preset area will be added as a defluorinating agent dosing point in the underground reservoir.
[0028] In addition to the aspects and any possible implementations described above, a further implementation is provided, wherein the method further includes:
[0029] The current removal rate is then visualized.
[0030] According to a second aspect of this disclosure, a defluoridation optimization device for underground water reservoirs is provided. The device includes:
[0031] The acquisition module is used to acquire parameters of the underground water reservoir in the mining area;
[0032] A module is established to create a three-dimensional underground reservoir model based on the underground reservoir parameters.
[0033] The setting module is used to set boundary conditions for the three-dimensional underground reservoir model. The boundary conditions include: water flow parameter boundaries and defluoridator parameter boundaries. The defluoridator parameter boundaries include: defluoridator dosage boundaries, defluoridator dosing points in the underground reservoir, and defluoridator dosing rate boundaries. The defluoridator is used to remove fluoride pollutants from the underground reservoir.
[0034] The operation module is used to run the three-dimensional underground reservoir model to obtain the current removal rate of the fluoride pollutants;
[0035] The judgment module is used to determine whether to adjust the three-dimensional underground reservoir model based on the current removal rate.
[0036] According to a third aspect of this disclosure, an electronic device is provided. The electronic device includes a memory and a processor, wherein the memory stores a computer program, and the processor executes the program to implement the method described above.
[0037] According to a fourth aspect of this disclosure, a computer-readable storage medium is provided having a computer program stored thereon that, when executed by a processor, implements the method according to a first aspect of this disclosure.
[0038] In this disclosure, after establishing a three-dimensional underground reservoir model based on the parameters of the underground reservoir, boundary conditions can be set for the three-dimensional underground reservoir model. Then, the three-dimensional underground reservoir model is run under the constraints of the boundary conditions to obtain the current removal rate of the fluoride pollutants. Based on the current removal rate, it is determined whether to adjust the three-dimensional underground reservoir model. In this way, by creating a simulated three-dimensional underground reservoir model for the underground reservoir and then continuously optimizing the three-dimensional underground reservoir model, the removal rate of fluoride pollutants can be continuously optimized, thereby continuously optimizing the fluoride removal effect of the underground reservoir.
[0039] It should be understood that the description in the Summary of the Invention section is not intended to limit the key or essential features of the embodiments of this disclosure, nor is it intended to restrict the scope of this disclosure. Other features of this disclosure will become readily apparent from the following description. Attached Figure Description
[0040] The above and other features, advantages, and aspects of the embodiments of this disclosure will become more apparent from the accompanying drawings and the following detailed description. The drawings are provided for a better understanding of the invention and are not intended to limit the scope of this disclosure. In the drawings, the same or similar reference numerals denote the same or similar elements, wherein:
[0041] Figure 1 A flowchart of an optimized method for defluoridation of an underground reservoir according to an embodiment of the present disclosure is shown;
[0042] Figure 2 A flowchart of another optimized method for defluoridation of underground reservoirs according to an embodiment of the present disclosure is shown;
[0043] Figure 3 A schematic diagram of a chemical dosing point in an underground reservoir according to an embodiment of the present disclosure is shown;
[0044] Figure 4 A block diagram of an underground reservoir defluoridation optimization device according to an embodiment of the present disclosure is shown;
[0045] Figure 5 A block diagram of an exemplary electronic device capable of implementing embodiments of the present disclosure is shown. Detailed Implementation
[0046] To make the objectives, technical solutions, and advantages of the embodiments of this disclosure clearer, the technical solutions of the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this disclosure, and not all embodiments. Based on the embodiments of this disclosure, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this disclosure.
[0047] Furthermore, the term "and / or" in this article is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, or B existing alone. Additionally, the character " / " in this article generally indicates that the preceding and following related objects have an "or" relationship.
[0048] Figure 1 A flowchart of an optimized method 100 for defluoridation of underground water reservoirs according to an embodiment of the present disclosure is shown. Method 100 may include:
[0049] Step 110: Obtain the parameters of the underground reservoir in the mining area;
[0050] The parameters of the underground reservoir include, but are not limited to: the size of the underground reservoir, the water area, the water level elevation, the lowest water level elevation, the average water depth, the water volume, the injection volume, and the drainage volume. They may also include: the first physical property parameters of the collapsed rock mass and the second physical property parameters of the grouting filling material of the collapsed rock mass. The first physical property parameters include, but are not limited to: the volume of the collapsed rock mass and the porosity of the collapsed rock mass. The second physical property parameters include, but are not limited to: the porosity, permeability coefficient, water yield, water storage rate, viscosity, flow velocity, and strength of the grouting filling material.
[0051] Step 120: Establish a three-dimensional underground reservoir model based on the underground reservoir parameters;
[0052] The 3D underground reservoir model is used to model the 3D form of an underground reservoir. Through this model, the simulated defluoridator dosing point, dosing amount, and dosing rate can be input, and the current removal rate of fluoride pollutants can be output. This allows us to measure whether the simulated parameters such as the dosing point, dosing amount, and dosing rate of the defluoridator are accurate. If accurate, the defluoridator can be added directly in the actual underground reservoir according to the simulated parameters.
[0053] Using groundwater and multiphase flow simulation software, including but not limited to MODFLOW, GMS, FEFLOW, and CFD, a three-dimensional groundwater reservoir model is established based on the aforementioned groundwater reservoir parameters. This model simulates three-dimensional units such as underground dams, pipelines, and filling materials within the groundwater reservoir. After adding defluoridating agents through well drilling and injection and inlet-side dosing (setting the dosing method and dosing point parameters for the three-dimensional groundwater reservoir model), the model is defined, including but not limited to groundwater seepage models, groundwater and agent solid-liquid multiphase flow models, adsorbent diffusion fluid dynamics models, adsorption kinetic models, and pollutant solute transport models.
[0054] Step 130: Set boundary conditions for the three-dimensional underground reservoir model. The boundary conditions include: water flow parameter boundaries and defluoridant parameter boundaries. The defluoridant parameter boundaries include: defluoridant dosage boundaries (i.e., the specific value of the defluoridant dosage or the range of the defluoridant dosage), defluoridant dosing points in the underground reservoir (i.e., the dosing locations where the defluoridant is placed in the underground reservoir), and defluoridant dosing rate boundaries (i.e., the defluoridant dosing rate value or the range of the defluoridant dosing rate). The defluoridant is used to remove fluoride pollutants from the underground reservoir.
[0055] The water flow parameter boundaries include, but are not limited to: the flow rate boundaries of the inlet and outlet of the underground reservoir or the pressure head boundaries of the inlet and outlet of the underground reservoir.
[0056] In the boundary of water flow parameters and the boundary of defluoridant parameters, the boundary can be understood as a range of values or a specific parameter value.
[0057] Step 140: Run the three-dimensional underground reservoir model to obtain the current removal rate of the fluoride pollutants;
[0058] The current removal rate can be equal to the ratio of the difference in (mass or) concentration of fluoride pollutants before and after adding the defluorinating agent to the ratio of the (mass or) concentration of fluoride pollutants before adding the defluorinating agent to the ratio of the (mass or) concentration of fluoride pollutants before adding the defluorinating agent, and is generally expressed as a percentage (%).
[0059] Step 150: Based on the current removal rate, determine whether to adjust the three-dimensional underground reservoir model.
[0060] After establishing a three-dimensional underground reservoir model based on the underground reservoir parameters, boundary conditions can be set for the three-dimensional underground reservoir model. Then, the three-dimensional underground reservoir model is run under the constraints of the boundary conditions to obtain the current removal rate of the fluoride pollutants. Based on the current removal rate, it is determined whether to adjust the three-dimensional underground reservoir model. In this way, by creating a simulated three-dimensional underground reservoir model for the underground reservoir and then continuously optimizing the three-dimensional underground reservoir model, the removal rate of fluoride pollutants can be continuously optimized, thereby continuously optimizing the fluoride removal effect of the underground reservoir.
[0061] In some embodiments, before running the three-dimensional underground reservoir model, the method further includes:
[0062] In the three-dimensional underground reservoir model, the concentration of fluoride pollutants at the inlet of the underground reservoir is set as a variable;
[0063] This variable can be assigned a value, and the value of this variable is adjustable.
[0064] In the three-dimensional underground reservoir model, the concentration of fluoride pollutants at the outlet of the underground reservoir is set as a monitoring variable;
[0065] Monitoring variables are quantities that need to be monitored regularly.
[0066] In the three-dimensional underground reservoir model, the first physical property parameter of the collapsed rock mass of the underground reservoir is assigned, the second physical property parameter of the grouting filling material of the collapsed rock mass is assigned, and the particle size of the defluorinating agent is set in the three-dimensional underground reservoir model;
[0067] Defluorinating agents are not limited to solid powders, granules, slurries, solutions, or other forms.
[0068] In the three-dimensional underground reservoir model, the collapsed rock mass and sidewalls of the underground reservoir are set as solid, non-slip boundaries.
[0069] The flow field of the underground reservoir can also be set to steady flow or unsteady flow in the 3D underground reservoir model, and other user-defined data can also be set in the 3D underground reservoir model.
[0070] In some embodiments, the water flow parameter boundaries include: the flow rate boundaries of the inlet and outlet of the underground reservoir (i.e., the specific values or ranges of the flow rate at the inlet and outlet) or the head boundaries of the inlet and outlet of the underground reservoir (i.e., the specific values or ranges of the head at the inlet and outlet). Head is a fluid mechanics concept that refers to the energy per unit weight of fluid.
[0071] The first physical property parameters include: the volume of the collapsed rock mass and the porosity of the collapsed rock mass;
[0072] The second physical property parameters include: porosity, permeability coefficient, water yield, water storage rate, viscosity, flow rate, and strength of the grouting filling material.
[0073] The collapsed rock mass refers to the rock that has collapsed from the underground reservoir itself, while the grouting filling material refers to the filling material injected into the rock mass.
[0074] In some embodiments, determining whether to adjust the three-dimensional underground reservoir model based on the current removal rate includes:
[0075] Obtain the target removal rate of fluoride pollutants from the three-dimensional underground reservoir model;
[0076] Compare the current removal rate with the target removal rate;
[0077] If the current removal rate is less than the target removal rate, then the boundary conditions of the three-dimensional underground reservoir model are adjusted.
[0078] If the current removal rate is less than the target removal rate, it indicates that the removal effect of fluoride pollutants is not good and the addition of defluorinating agent is insufficient. Therefore, the boundary conditions of the three-dimensional underground reservoir model can be adjusted to optimize the removal rate of fluoride pollutants by optimizing the three-dimensional underground reservoir model.
[0079] In some embodiments, the method further includes:
[0080] Run the three-dimensional underground reservoir model to obtain the flow field of the underground reservoir and the dispersion of the defluoridating agent in the underground reservoir;
[0081] The flow field can characterize which areas of the underground reservoir belong to the turbulent region and which belong to the laminar region, as well as the velocity and pressure changes at each location in the underground reservoir.
[0082] The dispersion of the defluorinating agent in the underground reservoir can be measured by the concentration of fluoride pollutants in the underground reservoir (assuming that the initial concentration of fluoride pollutants in the flow field of the underground reservoir is the same). After adding the defluorinating agent, the lower the concentration of fluoride pollutants, the better the defluorinating agent mixes with the water in the underground reservoir, and the better the dispersion of the defluorinating agent in the underground reservoir.
[0083] Based on the flow field conditions of the underground reservoir and / or the dispersion of the defluorinating agent in the underground reservoir, it is determined whether to adjust the defluorinating agent dosing point in the underground reservoir.
[0084] By utilizing the flow field conditions of the underground reservoir and / or the dispersion of the defluorinating agent in the underground reservoir, it can be accurately determined whether to adjust the defluorinating agent dosing point in the underground reservoir, such as whether it is necessary to add a new defluorinating agent dosing point in the underground reservoir.
[0085] In some embodiments, determining whether to adjust the defluoridator dosing point in the underground reservoir based on the flow field conditions of the underground reservoir and / or the dispersion of the defluoridating agent in the underground reservoir includes:
[0086] If the flow field in the first preset area of the underground reservoir is a turbulent region and the Reynolds number in the turbulent region is higher than a preset threshold, then the first preset area will be added as a defluoridation agent dosing point in the underground reservoir.
[0087] The Reynolds number is a similarity criterion number in fluid mechanics that characterizes the effects of viscosity.
[0088] If the flow field in the first preset area of the underground reservoir is a turbulent region and the Reynolds number in the turbulent region is higher than a preset threshold, it indicates that the defluorinating agent and water are highly mixed in the turbulent region and the defluorination effect is good. Therefore, the first preset area can be added as a defluorinating agent dosing point in the underground reservoir, such as adding a defluorinating agent dosing point in the first preset area.
[0089] or
[0090] In some embodiments, determining whether to adjust the defluoridator dosing point in the underground reservoir based on the flow field conditions of the underground reservoir and / or the dispersion of the defluoridating agent in the underground reservoir includes:
[0091] If the flow field in the second preset area of the underground reservoir is turbulent and the dispersion of the defluorinating agent in that preset area is higher than the preset agent dispersion, then the second preset area will be added as a defluorinating agent dosing point in the underground reservoir. The first preset area and the second preset area can be different areas of the underground reservoir.
[0092] If the flow field in the second preset area of the underground reservoir is a turbulent zone and the dispersion of the defluorinating agent in the preset area is higher than the preset agent dispersion, it indicates that the defluorinating agent and water are highly mixed in the turbulent zone and the defluorination effect is good. Therefore, the second preset area can be added as a defluorinating agent dosing point in the underground reservoir, such as adding a defluorinating agent dosing point in the second preset area.
[0093] In some embodiments, the method further includes:
[0094] The current removal rate is then visualized.
[0095] By visualizing the current removal rate, users can clearly see the current removal effect of the defluorinating agent on fluoride pollutants in the underground reservoir. It also makes it easier to statistically analyze the change curve of the current removal rate, which is beneficial for statistically analyzing the correspondence between various parameters such as the boundary conditions of the three-dimensional underground reservoir model and the removal rate of fluoride pollutants in the underground reservoir.
[0096] In summary, this disclosure provides a simulation solution for the dosing process and chemical reaction in the seepage field, and proposes that the pollutant change law of the underground reservoir at the dosing point of the defluorination agent is equivalent to the reverse migration law of pollutants migrating from the dosing point out of the flow field in the three-dimensional underground reservoir model.
[0097] The following will combine Figure 2 Further details of the technical solution disclosed herein:
[0098] 1. Data collection:
[0099] Based on the actual underground reservoir conditions in the mining area and the distribution of mined-out areas after mining activities, and combined with the underground reservoir water circulation and scheduling system, data on the underground reservoirs were obtained, including parameters such as reservoir size, water accumulation area, water level elevation, minimum water level elevation, average water accumulation depth, water volume, injection volume, and drainage volume. The physical properties of the underground reservoir's collapsed rock mass and grouting filling materials were also defined, including porosity, permeability coefficient, specific yield, and storage capacity.
[0100] 2. Mathematical model construction of underground reservoir:
[0101] Groundwater and multiphase flow simulation software, including but not limited to MODFLOW, GMS, FEFLOW, and CFD, were used to establish governing equations. A three-dimensional underground reservoir model was defined to simulate three-dimensional units such as underground dams, pipelines, and filling materials. After adding adsorbent via well drilling and inlet-side dosing, the following models were defined: including but not limited to groundwater seepage models, groundwater-adsorbent solid-liquid multiphase flow models, adsorbent diffusion fluid dynamics models, adsorption kinetic models, and solute transport models of pollutants.
[0102] Geometric modeling and mesh generation: Unstructured meshes are used for the flow field edges within the underground reservoir and adsorbent, with finer meshing applied to key areas such as dosing wells and inlet / outlet points. Discretization methods for the governing equations are not limited to the finite volume method, finite element method, or finite difference method.
[0103] 3. Set boundary conditions and assign parameter values in the 3D underground reservoir model:
[0104] The boundary conditions mainly consist of flow boundary and concentration boundary, determined based on the data collected in step 1. The flow boundary is set as the adjustable flow rate or head at the inlet and outlet of the groundwater reservoir; the pollutant concentration at the groundwater inlet is a variable; the well (hole) boundary is the chemical concentration boundary, with adjustable dosage, dosing rate, and dosing point; the fluoride concentration at the outlet is set as a monitored value; parameters such as reagent concentration, particle size, and porosity are set; the internal fractured rock mass and sidewalls of the reservoir are all set as solid, non-slip boundaries; based on chemical transport models, including diffusion, convection, and dispersion, the flow field is determined to be either steady or unsteady; and other user-defined data are also included.
[0105] 4. Simulation calculation & identification analysis:
[0106] By adjusting boundary conditions such as influent and effluent flow velocities, the groundwater flow field, pollutant removal rate (with local concentration monitoring), reagent dispersion, and dosage were initially determined. Specific elevations and locations of laminar and turbulent flow zones were identified. Based on the characteristic regions of laminar and turbulent flow in the flow field, and considering the defluorination efficiency per unit residence time of the defluorination agent (measured by fluoride concentration; the faster the fluoride concentration decreases, the higher the defluorination efficiency), the locations of newly added wells were designed. The model was continuously optimized through boundary parameter callbacks until a solution meeting the requirements was obtained.
[0107] 5. Post-processing data analysis and visualization of simulation results.
[0108] The technical solution of this disclosure will be further illustrated below through embodiments:
[0109] Example 1: A CFD three-dimensional single-phase flow model was established to simulate the single-phase flow of fluid in an underground reservoir containing collapsed rock mass structures. The flow field of the groundwater was obtained, and the locations of typical laminar and turbulent regions in the underground reservoir were characterized. The dosing point was set at the fluid inlet; the boundary conditions were the fluid inlet and outlet velocities and the dosing rate. The outlet boundary adopted a free outflow boundary; the internal fractured rock mass and sidewalls of the reservoir were all set as solid walls with no slippage; the finite volume method was used to discretize the governing equations.
[0110] In this embodiment, the adsorbent is a powdered agent. The settling rate of the adsorbent powder in the laminar flow zone and the mixing rate in the turbulent flow zone of the underground reservoir can be reasonably estimated. The residence time of the agent can be calculated by combining the fluid velocity in the flow field, and the defluorination efficiency of a certain dosage can be further calculated.
[0111] Example 2: A CFD three-dimensional multiphase flow model was established to simulate the solid-liquid two-phase flow model of fluid and particulate adsorbent in an underground reservoir containing collapsed rock mass. The water phase of the underground reservoir is the main continuous phase, and the adsorbent particles are the secondary dispersed phase. The material properties are isotropic. The flow field of the groundwater was obtained, and the mixing of the solid and liquid phases and the collision distribution of particles in the underground reservoir were characterized. The outlet boundary adopts a free outflow boundary. The internal fractured rock mass and sidewalls of the reservoir are all set as solid walls without slippage. The finite volume method is used as the discretization method for the governing equations.
[0112] In this embodiment, the adsorbent is a granular agent, and the dosing method is pulsed gas-solid delivery dosing. The adjustable parameters are the dosing amount, dosing frequency, and pulse airflow velocity. The boundary conditions are the fluid inlet and outlet flow rates and the dosing amount. Combining the defluorination performance of the adsorbent itself, the contact residence time between a certain flow rate of groundwater and the agent in the underground reservoir is determined, and the defluorination efficiency of a certain dosing amount is further calculated.
[0113] Example 3: A three-dimensional FEFLOW groundwater reservoir seepage model was established to simulate the fluid flow in a groundwater reservoir filled with gangue, obtaining a visualized groundwater flow field. Based on the visualized flow field, combined with a solute transport model, the transport of pollutants from the flow field to the dosing point was simulated, the pollutant concentration in the flow field was calculated, and the dosing point and dosage were determined accordingly. The boundary conditions in this example are the fluid inlet and outlet velocities and the pollutant concentration. Combining the defluoridation performance of the adsorbent itself, the contact residence time between groundwater at a certain flow rate and the reagent in the groundwater reservoir was determined, and the defluoridation efficiency at a certain dosage was further calculated. The model uses an unstructured mesh, with finer meshes in key areas such as dosing wells and inlet / outlet points. The finite element method was used to discretize the governing equations.
[0114] Figure 3 The diagram simply illustrates the dosing point for defluoridating agents in the underground water reservoir of a coal mine.
[0115] Figure 3Points 1-8 are designated as defluoridation agent dosing points for the underground reservoir, which can be used or closed. Point 1 is the reservoir inlet, and point 9 is the reservoir outlet. The chemicals used at the dosing points come from the dosing system. The defluoridation agent is not limited to solid powder, granules, slurry, or solution forms, and the dosing method is not limited to pneumatic conveying, screw pumps, or water pumps.
[0116] It should be noted that, for the sake of simplicity, the foregoing method embodiments are all described as a series of actions. However, those skilled in the art should understand that this disclosure is not limited to the described order of actions, because according to this disclosure, some steps can be performed in other orders or simultaneously. Furthermore, those skilled in the art should also understand that the embodiments described in the specification are all optional embodiments, and the actions and modules involved are not necessarily essential to this disclosure.
[0117] The above is an introduction to the method embodiments. The following describes the solution described in this disclosure further through device embodiments.
[0118] Figure 4 A block diagram of an underground reservoir defluoridation optimization device 400 according to an embodiment of the present disclosure is shown. Figure 4 As shown, the device 400 includes:
[0119] Module 410 is used to acquire parameters of the underground water reservoir in the mining area;
[0120] Module 420 is used to establish a three-dimensional underground reservoir model based on the underground reservoir parameters;
[0121] Setting module 430 is used to set boundary conditions for the three-dimensional underground reservoir model. The boundary conditions include: water flow parameter boundaries and defluoridator parameter boundaries. The defluoridator parameter boundaries include: defluoridator dosage boundaries, defluoridator dosing points in the underground reservoir, and defluoridator dosing rate boundaries. The defluoridator is used to remove fluoride pollutants from the underground reservoir.
[0122] The running module 440 is used to run the three-dimensional underground reservoir model to obtain the current removal rate of the fluoride pollutants;
[0123] The judgment module 450 is used to determine whether to adjust the three-dimensional underground reservoir model based on the current removal rate.
[0124] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working process of the described module can be referred to the corresponding process in the foregoing method embodiments, and will not be repeated here.
[0125] According to embodiments of the present disclosure, the present disclosure also provides an electronic device and a non-transitory computer-readable storage medium storing computer instructions.
[0126] Figure 5 A schematic block diagram of an electronic device 800 that can be used to implement embodiments of the present disclosure is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device may also represent various forms of mobile devices, such as personal digital processors, cellular phones, smartphones, wearable devices, and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the present disclosure described and / or claimed herein.
[0127] Device 800 includes a computing unit 801, which can perform various appropriate actions and processes based on a computer program stored in read-only memory (ROM) 802 or a computer program loaded from storage unit 808 into random access memory (RAM) 803. RAM 803 may also store various programs and data required for the operation of device 800. The computing unit 801, ROM 802, and RAM 803 are interconnected via bus 804. Input / output (I / O) interface 805 is also connected to bus 804.
[0128] Multiple components in device 800 are connected to I / O interface 805, including: input unit 806, such as keyboard, mouse, etc.; output unit 807, such as various types of monitors, speakers, etc.; storage unit 808, such as disk, optical disk, etc.; and communication unit 809, such as network card, modem, wireless transceiver, etc. Communication unit 809 allows device 800 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.
[0129] The computing unit 801 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of the computing unit 801 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various computing units running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. The computing unit 801 performs the various methods and processes described above, such as method 100. For example, in some embodiments, method 100 may be implemented as a computer software program tangibly contained in a machine-readable medium, such as storage unit 808. In some embodiments, part or all of the computer program may be loaded and / or installed on device 800 via ROM 802 and / or communication unit 809. When the computer program is loaded into RAM 803 and executed by the computing unit 801, one or more steps of method 100 described above may be performed. Alternatively, in other embodiments, the computing unit 801 may be configured to perform method 100 by any other suitable means (e.g., by means of firmware).
[0130] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems-on-a-chip (SoCs), payload-programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.
[0131] The program code used to implement the methods of this disclosure may be written in any combination of one or more programming languages. This program code may be provided to a processor or controller of a general-purpose computer, special-purpose computer, or other programmable data processing apparatus, such that when executed by the processor or controller, the program code causes the functions / operations specified in the flowcharts and / or block diagrams to be implemented. The program code may be executed entirely on a machine, partially on a machine, as a standalone software package partially on a machine and partially on a remote machine, or entirely on a remote machine or server.
[0132] In the context of this disclosure, a machine-readable medium can be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, apparatus, or device. A machine-readable medium can be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium can be, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination of the foregoing. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the foregoing.
[0133] To provide interaction with a user, the systems and techniques described herein can be implemented on a computer having: a display device for displaying information to the user (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor); and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the computer. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).
[0134] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as a data server), or computing systems that include middleware components (e.g., an application server), or computing systems that include frontend components (e.g., a user computer with a graphical user interface or web browser through which a user can interact with embodiments of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., a communication network). Examples of communication networks include local area networks (LANs), wide area networks (WANs), and the Internet.
[0135] Computing systems can include clients and servers. Clients and servers are generally located far apart and typically interact via communication networks. Client-server relationships are created by computer programs running on the respective computers and having a client-server relationship with each other. Servers can be cloud servers, servers in distributed systems, or servers incorporating blockchain technology.
[0136] It should be understood that the various forms of processes shown above can be used to rearrange, add, or delete steps. For example, the steps described in this disclosure can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution disclosed in this disclosure can be achieved, and this is not limited herein.
[0137] The specific embodiments described above do not constitute a limitation on the scope of protection of this disclosure. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this disclosure should be included within the scope of protection of this disclosure.
Claims
1. An optimized method for defluoridation in underground water reservoirs, characterized in that, include: Obtain parameters of the underground water reservoir in the mining area; Based on the parameters of the underground reservoir, a three-dimensional underground reservoir model is established; Boundary conditions are set for the three-dimensional underground reservoir model, wherein the boundary conditions include: water flow parameter boundaries and defluoridator parameter boundaries, wherein the defluoridator parameter boundaries include: defluoridator dosage boundaries, defluoridator dosing points in the underground reservoir, and defluoridator dosing rate boundaries, and the defluoridator is used to remove fluoride pollutants in the underground reservoir; Run the three-dimensional underground reservoir model to obtain the current removal rate of the fluoride pollutant; Based on the current removal rate, determine whether to adjust the three-dimensional underground reservoir model.
2. The method according to claim 1, characterized in that, Before running the three-dimensional underground reservoir model, the method further includes: In the three-dimensional underground reservoir model, the concentration of fluoride pollutants at the inlet of the underground reservoir is set as a variable; In the three-dimensional underground reservoir model, the concentration of fluoride pollutants at the outlet of the underground reservoir is set as a monitoring variable; In the three-dimensional underground reservoir model, the first physical property parameter of the collapsed rock mass of the underground reservoir is assigned, the second physical property parameter of the grouting filling material of the collapsed rock mass is assigned, and the particle size of the defluorinating agent is set in the three-dimensional underground reservoir model; In the three-dimensional underground reservoir model, the collapsed rock mass and sidewalls of the underground reservoir are set as solid, non-slip boundaries.
3. The method according to claim 2, characterized in that, The water flow parameter boundaries include: the flow rate boundaries of the inlet and outlet of the underground reservoir or the pressure head boundaries of the inlet and outlet of the underground reservoir. The first physical property parameters include: the volume of the collapsed rock mass and the porosity of the collapsed rock mass; The second physical property parameters include: porosity, permeability coefficient, water yield, water storage rate, viscosity, flow rate, and strength of the grouting filling material.
4. The method according to claim 1, characterized in that, The step of determining whether to adjust the three-dimensional underground reservoir model based on the current removal rate includes: Obtain the target removal rate of fluoride pollutants from the three-dimensional underground reservoir model; Compare the current removal rate with the target removal rate; If the current removal rate is less than the target removal rate, then the boundary conditions of the three-dimensional underground reservoir model are adjusted.
5. The method according to claim 1, characterized in that, The method further includes: Run the three-dimensional underground reservoir model to obtain the flow field of the underground reservoir and the dispersion of the defluoridating agent in the underground reservoir; Based on the flow field conditions of the underground reservoir and / or the dispersion of the defluorinating agent in the underground reservoir, it is determined whether to adjust the defluorinating agent dosing point in the underground reservoir.
6. The method according to claim 1, characterized in that, The step of determining whether to adjust the defluoridator dosing point in the underground reservoir based on the flow field conditions of the underground reservoir and / or the dispersion of the defluoridating agent in the underground reservoir includes: If the flow field in the first preset area of the underground reservoir is turbulent and the Reynolds number in the turbulent area is higher than a preset threshold, then the first preset area will be added as a defluoridation agent dosing point in the underground reservoir; or If the flow field in the second preset area of the underground reservoir is a turbulent zone and the dispersion of the defluorinating agent in the preset area is higher than the preset agent dispersion, then the second preset area will be added as a defluorinating agent dosing point in the underground reservoir.
7. The method according to any one of claims 1 to 6, characterized in that, The method further includes: The current removal rate is then visualized.
8. A defluoridation optimization device for underground reservoirs, characterized in that, include: The acquisition module is used to acquire parameters of the underground water reservoir in the mining area; A module is established to create a three-dimensional underground reservoir model based on the underground reservoir parameters. The setting module is used to set boundary conditions for the three-dimensional underground reservoir model. The boundary conditions include: water flow parameter boundaries and defluoridator parameter boundaries. The defluoridator parameter boundaries include: defluoridator dosage boundaries, defluoridator dosing points in the underground reservoir, and defluoridator dosing rate boundaries. The defluoridator is used to remove fluoride pollutants from the underground reservoir. The operation module is used to run the three-dimensional underground reservoir model to obtain the current removal rate of the fluoride pollutants; The judgment module is used to determine whether to adjust the three-dimensional underground reservoir model based on the current removal rate.
9. An electronic device, characterized in that, include: At least one processor; as well as A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the method of any one of claims 1-7.
10. A non-transitory computer-readable storage medium storing computer instructions, characterized in that, The computer instructions are used to cause the computer to perform the method according to any one of claims 1-7.
Citation Information
Patent Citations
Treatment method of water body in coal mine underground reservoir
CN104291399A
Method for determining water storage coefficient of coal mine underground water reservoir
CN108169093A
Construction method for three-dimensional numerical model of underground water in mining area
CN114357839A
Fluorine ion concentration prediction method, fluorine removal method and fluorine removal device
CN115712990A
System and method for dissolving gases in fluids and for delivery of dissolved gases
WO2006127979A2