A method and device for analyzing and calculating the sedimentation state of groundwater recharge particles

Abstract

The present invention discloses an analytical calculation method and device for the sedimentation state of groundwater recharge particles, and the analytical calculation method includes: S1. By combining the lattice Boltzmann method with the momentum equation, the motion variable equation of the particles is constructed to obtain the density distribution function of the particles moving with the water flow and the density distribution function of the particles after colliding with the pore wall; S2. Based on the two density distribution functions of the particles combined with the Euler equation model, an analytical model is constructed; S3. The particle density distribution function value in the pores inside the actual medium is detected and obtained, and the actual particle density distribution function value is input into the analytical model to calculate the sedimentation state of the suspended particles in the water body. Based on the particle movement process and force conditions inside the saturated porous medium, the present invention combines the lattice Boltzmann method and the momentum theorem, and can accurately simulate the movement trajectory of the particles in the pores, and then judge the sedimentation state of the suspended particles, so as to facilitate dynamic regulation during the artificial recharge of groundwater and avoid particle blockage.

Description

Technical Field

[0001] The invention relates to the technical field of groundwater recharge engineering, and in particular to an analytical calculation method and equipment for the sedimentation state of groundwater recharge particles. Background Art

[0002] The microscopic mechanism of the occurrence and development of particle blockage in artificial groundwater recharge is that the deposition of particles in the pores inside the medium reduces the pore space and blocks the pore connectivity, among which particle blockage is the macroscopic manifestation of particle deposition, and particle deposition is the microscopic mechanism of particle blockage. The particle deposition process in the pores inside the medium (microscopic scale) controls the particle blockage process of the medium as a whole (macroscopic scale). Therefore, studying the microscopic mechanism of particle deposition corresponding to the macroscopic phenomenon of particle blockage during groundwater recharge will help deepen the understanding of groundwater recharge blockage, and it is also beneficial to reasonably avoid the adverse effects of particle blockage on recharge projects during the design and operation of groundwater recharge projects.

[0003] Particle blockage is mainly caused by the deposition of suspended particles in the pores inside the medium when water flows through the porous medium. This process involves multiple factors such as particle movement, fluid mechanics properties of the medium, and interactions between particles. Particle deposition not only reduces the pore space, but may also lead to changes in the pore structure, thereby affecting the flow path of the fluid and the permeability of water.

[0004] Studies have found that the properties of particles (such as particle size, shape, density) and the characteristics of fluids (such as flow rate, viscosity) have a significant impact on the particle deposition state. In addition, environmental factors such as water quality, temperature and pH value also play an important role in the deposition behavior of particles. In the study of the microscopic mechanism of groundwater recharge, the process of particle deposition is considered to be a key factor affecting the seepage of the medium. The deposition of particles is not only related to the shear force, gravity and buoyancy of the fluid, but also closely related to the porosity and pore structure of the medium. Only by studying the deposition mechanism of particles in the medium can we effectively predict the particle clogging phenomenon and optimize the design of groundwater recharge projects to avoid the negative impact of particle clogging on groundwater recharge projects.

[0005] Although there are many numerical models and experimental studies in the prior art to explore the mechanism of particle deposition, the existing models are often not based on the movement process of the particles themselves. They often generalize the deposition process of a single particle as an overall movement process involving a huge number of particles. They cannot accurately describe the movement process of a single particle at the microscopic level, and thus cannot calculate and determine whether a single suspended particle in a porous medium is deposited or not.

[0006] To this end, the present application proposes a method and device for analyzing and calculating the deposition state of groundwater reinjection particles, which can solve the above-mentioned technical problems by calculating and judging the deposition state of particles in the pores inside the medium. Summary of the invention

[0007] The main purpose of the present invention is to provide a method and device for analytical calculation of the deposition state of groundwater reinjection particles. Based on the particle movement process and force conditions inside a saturated porous medium, the movement process of suspended particles in the pores can be calculated, and an analytical model for determining whether particles are deposited or not can be established to solve the technical problems raised in the background technology.

[0008] The present invention adopts the following technical solutions to solve the above technical problems:

[0009] An analytical calculation method for groundwater recharge particle deposition state, comprising:

[0010] S1. By combining the lattice Boltzmann method with the momentum equation, the equation of motion variables of the particles is constructed to obtain the density distribution function of the particles moving with the water flow and the density distribution function of the particles after colliding with the pore wall;

[0011] S2. An analytical model is constructed based on the density distribution function of particles moving with water flow and the density distribution function of particles after collision with pore walls combined with the Euler equation model;

[0012] S3. Detect and obtain the particle density distribution function value in the pores inside the actual medium, input the actual particle density distribution function value into the analytical model, and calculate the deposition state of suspended particles in the water body.

[0013] Preferably, the specific expression of the motion variable equation in step S1 is:

[0014]

[0015] in, For particles in a specified medium The density distribution function in is For particles in a specified medium Single pore inside The density distribution function in is For particles in a specified medium The forces acting on the inside is the spatial position of the particle, is the time node where the particle is located, is the time that the particle is pulled by the force The step length, For particles in a specified medium The speed of movement within a specified time step, is the particle size in time step Movement distance within Generally, it is set as the space step size, is the dimensionless relaxation time.

[0016] Preferably, the density distribution function of the particles moving with the water flow and the density distribution function of the particles after colliding with the pore wall in step S1 specifically include:

[0017] (1) Particles move in pores with water flow:

[0018]

[0019] (2) Particles collide with the pore wall along with the water flow:

[0020]

[0021] in, For the specified medium The space where the collision process occurs The step length, The time when the collision occurred The step length, is the density distribution function after the particles collide with the pore wall.

[0022] Preferably, the specific construction process of the analytical model in step S2 includes: converting the density distribution function and Substituting the solid-liquid two-phase flow continuity equation in the Euler-Euler model, the particle concentration at different heights of the sand column during the reinjection process can be calculated as the deposition state of suspended particles in the water body. The specific expression of the solid-liquid two-phase flow continuity equation is:

[0023]

[0024] in, represents the liquid or solid phase, for The volume concentration of the phase, for The velocity vector of the phase, for The density of the phase, For time, is the Hamiltonian differential operator.

[0025] Preferably, a device for analytical calculation of the sedimentation state of groundwater recharge particles has an analytical model stored therein, wherein the analytical model is constructed based on any of the above-mentioned analytical calculation methods for the sedimentation state of groundwater recharge particles, and the analytical model is also used to input the material properties of water bodies and suspended particles to calculate the movement process of suspended particles.

[0026] Preferably, the specific calculation process of the suspended particle movement process includes:

[0027] L1. Simulate the flow field of water in the pores inside the analytical medium, which includes the flow velocity and the interaction force on each spatial particle of the water flow;

[0028] L2. Calculate the total vector force on the particles in the pores of the sand column based on the water flow field. The total vector force includes the pore water flow drag force and the electric field force.

[0029] L3. Based on the total vector force on the particles, the particle acceleration of the entire time process is obtained as the movement process of the particles in the pores inside the sand column with the water flow.

[0030] Preferably, the specific operation process in the L1 step includes:

[0031] L11. Calculation of Reynolds number , the specific calculation formula is:

[0032]

[0033] in, is the flow velocity of the seepage inside the pores, is the pore diameter, is the kinematic viscosity of water;

[0034] L12. Determine the flow pattern of water inside the pores of the porous medium as Stokes flow with a Reynolds number less than 1, and set the microscopic pore water flow field in the pores inside the medium as Stokes flow with neglected inertial force, which is:

[0035]

[0036]

[0037]

[0038] in, Specify the velocity vector at the particle for the fluid, Specify the pressure vector at the particle for the fluid, is the Hamiltonian differential operator, , , are the three coordinate vectors for the velocity vector in the three-dimensional coordinate system;

[0039] L13. Based on the Stokes equation to characterize the water flow velocity at each spatial position in the pores inside the medium in the groundwater recharge project, there exists:

[0040]

[0041]

[0042] in, is the water body space step length, is the space step length The water flow vector velocity inside is the space step The vector force acting on the internal water flow is the sum of the external forces.

[0043] Preferably, the specific calculation expression of the water flow drag force in the L2 step is:

[0044]

[0045] in, is the drag force on the particle by the water flow, is the mass of the particle, is the density of the particle material, is the diameter of the particle, is the vector velocity of the specified particle.

[0046] Preferably, the electric field force in the L2 step is calculated based on Coulomb's law, and the differential formula of Coulomb's law used at this time is as follows:

[0047]

[0048]

[0049] in, is the electric displacement vector, is the charge density, is the vacuum dielectric constant, which is measured , is the dielectric constant of water, which is temperature-dependent. is the electric field strength.

[0050] Preferably, the particle acceleration during the entire time process in the L3 step is The specific calculation expression is:

[0051]

[0052] in, is the mass of the particle, is the vector velocity of the specified particle, For particles at a specified time step The total vector force on Including the water flow drag force of particles, gravity, water buoyancy and electric field force.

[0053] On the other hand, the present invention further discloses a computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the processor executes the steps of the above method.

[0054] On the other hand, the present invention further discloses a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the above method.

[0055] It can be seen from the above technical solution that the present invention provides a method and device for analyzing and calculating the sedimentation state of groundwater recharge particles. Compared with the prior art, the present invention has the following advantages:

[0056] 1. The present invention can accurately simulate the movement trajectory of particles in the pores, and by calculating the movement process of the particles, further determine the sedimentation state of the suspended particles, thereby facilitating dynamic regulation during the artificial recharge of groundwater and avoiding particle blockage.

[0057] 2. The present invention establishes an analytical model for judging whether particles are deposited based on the stress conditions and movement process of suspended particles inside saturated porous media. Compared with most current models based on medium particle blockage phenomena (indirect data), the present invention provides a particle deposition calculation method based on the particles themselves (direct data).

[0058] It should be understood that the content described in this section is not intended to identify the key or important features of the embodiments of the present invention, nor is it intended to limit the scope of the present invention. Other features of the present invention will become easy to understand through the following description. Of course, it is not necessary to achieve all of the advantages described above simultaneously for any product implementing the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] The drawings constituting a part of the present application are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0060] Figure 1 It is a schematic diagram of the overall process of the method of the present invention;

[0061] Figure 2 It is a schematic diagram of the overall process of the particle motion process calculation method of the present invention;

[0062] Figure 3 It is a schematic diagram of the correlation curve between water temperature and water body kinematic viscosity coefficient;

[0063] Figure 4 It is a schematic diagram of the trajectory of particle movement process and the three-dimensional path of the particle deposition state model of the present invention;

[0064] Figure 5 It is a schematic diagram of the particle movement process trajectory and particle deposition state model in the z-axis plane top view of the present invention;

[0065] Figure 6 It is a schematic diagram of the particle movement process trajectory and particle deposition state model of the present invention in the side view of the x-axis plane;

[0066] Figure 7 It is a schematic diagram of the particle movement process trajectory and particle deposition state model of the present invention in the y-axis plane side view;

[0067] Figure 8 A close-up view of the deposition state of some particles of the present invention;

[0068] Fig. 9 This is a close-up of the trajectory of some particles in the present invention. DETAILED DESCRIPTION

[0069] The technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all of the embodiments. In the absence of conflict, the embodiments in this application and the features in the embodiments can be combined with each other. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0070] In the embodiment, see Figures 1 to 9 .

[0071] In the prior art, during the artificial recharge of groundwater, the sedimentation behavior of particles under the action of hydrodynamics will affect the pore structure and flow characteristics of the medium, thereby causing particle clogging. Therefore, an embodiment of the present invention proposes an analytical calculation method for the sedimentation state of groundwater recharge particles. By constructing an analytical model based on the lattice Boltzmann method and the momentum theorem, the movement trajectory of particles in the pores can be accurately simulated, and the sedimentation state can be judged by calculation, thereby controlling and avoiding particle clogging.

[0072] Among them, Figure 1 As shown, the analytical calculation method is to calculate the deposition of suspended particles based on the individual resolution and overall number statistics of particles moving in the pores inside the medium. The specific steps include:

[0073] S1. Suspended particles in the pores of saturated porous media move in the pores of the medium with the water flow under the action of hydrodynamic force. Because their bulk density is greater than that of water, they will deviate from the pore water streamline under the action of inertia and collide with the pore wall to generate momentum loss. When the momentum of the particle is reduced to zero and the direction of the superimposed force points to the pore wall, the particle will be deposited on the pore wall. Therefore, the particle motion variable equation is constructed by combining the lattice Boltzmann method with the momentum equation to obtain the density distribution function of the particle moving with the water flow and the density distribution function of the particle after colliding with the pore wall.

[0074] The specific expression of the motion variable equation is:

[0075]

[0076] in, For particles in a specified medium The density distribution function in is For particles in a specified medium Single pore inside The density distribution function in is For particles in a specified medium The forces acting on the inside is the spatial position of the particle, is the time node where the particle is located, is the time that the particle is pulled by the force The step length, For particles in a specified medium The speed of movement within a specified time step, is the particle size in time step Movement distance within Generally, it is set as the space step size, is the dimensionless relaxation time;

[0077] At this time, the density distribution function of the particles moving with the water flow and the density distribution function of the particles after colliding with the pore wall specifically include:

[0078] (1) Particles move in pores with water flow:

[0079]

[0080] (2) Particles collide with the pore wall along with the water flow:

[0081]

[0082] in, For the specified medium The space where the collision process occurs The step length, The time when the collision occurred The step length, is the density distribution function after the particles collide with the pore wall.

[0083] S2. Construct an analytical model based on the density distribution function of two particle motions (density distribution function of particles moving with water flow and density distribution function of particles after collision with pore walls) combined with the Euler equation model;

[0084] The specific construction process of the analytical model includes: and By substituting the solid-liquid two-phase flow continuity equation in the Euler-Euler model, the particle concentration at different heights of the sand column during the reinjection process can be calculated as the deposition state of suspended particles in the water body. At this time, the Euler-Euler model treats substances of different phases as interpenetrating continua, which is suitable for calculating the research scenario of particles moving with water flow.

[0085] The specific expression of the solid-liquid two-phase flow continuity equation is:

[0086]

[0087] in, represents the liquid or solid phase, for The volume concentration of the phase, for The velocity vector of the phase, for The density of the phase, For time, is the Hamiltonian differential operator.

[0088] At this point, it can be further explained that since the particles in the water body are solid phase, the density is quite different from that of the water body as a liquid phase. In the process of describing and calculating the movement and collision of particles with the flow of water, it is necessary to define separate momentum equations for the solid phase and the liquid phase respectively. The Euler-Euler model can realize this independent momentum definition for each phase. The name of the Euler-Euler model comes from the fact that the water body (liquid phase) and the particles in the water body (solid phase) are generalized as continua and described by the Euler equation.

[0089] It is further explained here that water is a continuous phase in actual physical scenarios, but the particles in the water are not continuous in actual physical scenarios, but a dispersed phase in which each particle is distributed separately. The dispersed phase is generalized as a continuous phase in the Euler-Euler model only to make the computational model have better robustness and convergence, and it can also save computer computing resources.

[0090] For physical scenarios with a small number of particles and a relatively simple pore structure, if the computing power of the computer workstation used is sufficient, the Euler equation can be used to describe the water body as the continuous phase, and the Lagrangian equation can be used to describe the particles as the dispersed phase. Together, they form the Euler-Lagrangian model. This computational model is closer to the actual physical scenario than the Euler-Euler model, but the convergence conditions are more demanding.

[0091] Both the Euler equation and the Lagrange equation are classic equations in the finite element method. The finite element method, also known as the finite unit method, is a numerical solution method that uses mathematical approximate analysis methods to simulate the solution object, defines the basic elements of interaction in a complex system as units, and uses a finite number of known units to approximate the real system with infinite unknown quantities. Specifically, a finite element model is established based on the known quantities, and an approximate solution to the real system is obtained by solving the mathematical model. In theory, the more known quantities used in modeling, the closer the approximate solution is to the actual solution.

[0092] In summary, the Euler equation is to comprehensively analyze the overall motion of a certain item in the entire research space by defining the motion state of particles at different spatial nodes. The definition of the motion of particles is composed of the velocity component and acceleration component of each point at a certain moment.

[0093] Specifically, the particle in the Euler equation is Always , , The velocity components in the direction are as follows:

[0094]

[0095]

[0096]

[0097] In Euler's equation, the particle Always , , The components of acceleration in the direction are as follows:

[0098]

[0099]

[0100]

[0101] The Lagrange equation uses mathematical description to determine the motion trajectory of each particle that makes up a specified phase. By integrating the motion trajectory of a single particle, the motion law of the entire phase can be obtained. The definition of the particle motion trajectory is determined by the initial moment The spatial coordinates of , , ) and any moment The spatial coordinates of are composed together, where the Lagrange equation is as follows:

[0102]

[0103]

[0104]

[0105] In the above formula, , , , are collectively referred to as Lagrange variables, ( , , ) represents a particle of the phase at the initial moment The spatial coordinates of , , ) is the particle in The spatial coordinates at a certain moment in time, the Lagrange equation describes the position information of each particle of a specified phase throughout the entire study time interval, and the connection in series is the motion trajectory of all the particles that make up the phase.

[0106] At this time, high-precision non-destructive testing technology can be used to identify particles in the internal pores of the medium and thus obtain the number of particles in a certain spatial area, that is, the particle density distribution function value.

[0107] S3. Detect and obtain the particle density distribution function value in the pores inside the actual medium, input the actual particle density distribution function value into the analytical model, and calculate the deposition state of suspended particles in the water body.

[0108] In summary, by analyzing the stress conditions of suspended particles inside the pores of saturated porous media, based on the lattice Boltzmann method and the momentum theorem, we can accurately simulate the movement trajectory of particles in the pores, and by calculating the movement process of the particles, we can further determine the deposition state of the suspended particles, thereby facilitating dynamic regulation during artificial groundwater recharge and avoiding particle blockage.

[0109] On the other hand, the present invention also discloses an analytical calculation device for the sedimentation state of groundwater recharge particles, which stores an analytical model internally. The analytical model uses multi-physics field simulation software to establish a research area describing the pores in a saturated porous medium on a COMSOL Multiphysics platform, and is constructed based on the analytical calculation method for the sedimentation state of groundwater recharge particles. The analytical model is also used to input the material properties of water bodies and suspended particles to calculate the movement process of suspended particles. Figure 2 As shown in Figure 2, the specific calculation process of the suspended particle movement process includes:

[0110] L1. Simulate the water flow field in the pores inside the analytical medium. The water flow field includes the flow velocity and the interaction force on each spatial particle of the water flow.

[0111] The specific operation process at this time includes:

[0112] L11. Calculation of Reynolds number , the specific calculation formula is:

[0113]

[0114] in, is the flow velocity of the seepage inside the pores, is the pore diameter, is the kinematic viscosity of water;

[0115] Here, by substituting the temperature into Figure 3 The correlation curve between water temperature and water body kinematic viscosity coefficient can be used to obtain the kinematic viscosity coefficient. Specific value of

[0116] L12. Ensure that the flow pattern of water inside the pores of saturated porous media is Stokes flow with a Reynolds number less than 1, which is also the main flow pattern of saturated pore water flow inside geological bodies with particle blockage during artificial groundwater recharge;

[0117] As can be seen from the L11 step, the physical meaning of the Reynolds number is the ratio of the inertial force and the viscous force of the water flow. In the engineering scenario of artificial groundwater recharge, the viscous force of the water flow field in the medium is generally much greater than the inertial force. Therefore, the microscopic pore water flow field in the pores inside the medium is set to Stokes flow with neglected inertial force, which is:

[0118]

[0119]

[0120]

[0121] in, Specify the velocity vector at the particle for the fluid, Specify the pressure vector at the particle for the fluid, is the Hamiltonian differential operator, , , are the three coordinate vectors for the velocity vector in the three-dimensional coordinate system;

[0122] L13. Based on the Stokes equation to characterize the water flow velocity at each spatial position in the pores inside the medium in the groundwater recharge project, there exists:

[0123]

[0124]

[0125] in, is the water body space step length, is the space step length The water flow vector velocity inside is the space step The vector force acting on the internal water flow is the sum of the external forces.

[0126] L2. Calculate the total vector force on the particles in the pores of the sand column based on the water flow field. The total vector force includes the pore water flow drag force and the electric field force.

[0127] in:

[0128] (1) Calculate the drag force of pore water flow on particles in the pores of the sand column. The specific calculation expression is:

[0129]

[0130] in, is the drag force on the particle by the water flow, is the mass of the particle, is the density of the particle material, is the diameter of the particle, is the vector velocity of the specified particle. It is further defined as the vector velocity of the water body within the spatial range that interacts with the specified particle and is obtained based on the water flow field;

[0131] (2) The electric field force on the particle is calculated using Coulomb's law. The differential formula of Coulomb's law used is as follows:

[0132]

[0133]

[0134] in, is the electric displacement vector, is the charge density, is the vacuum dielectric constant, which is measured , is the dielectric constant of water, which is temperature-dependent. is the electric field strength.

[0135] L3. Based on the total vector force on the particles, the force on the particles is substituted into Newton's second law to obtain the particle acceleration over the entire time process, which is the movement process of the particles in the pores inside the sand column along with the water flow.

[0136] The particle acceleration during the whole time process is The specific calculation expression is:

[0137]

[0138] in, is the mass of the particle, is the vector velocity of the specified particle, For particles at a specified time step The total vector force on Including the water flow drag force of particles, gravity, water buoyancy and electric field force.

[0139] Through the above expression, the vector acceleration of the particles at each time step can be calculated. With the acceleration of all particles corresponding to each time step, the movement of particles in the whole process of groundwater recharge can be described.

[0140] In summary, the analytical model can calculate the particle deposition state under the condition of directly observing the movement of particles inside the medium during the test. An analytical model for calculating and judging whether particles are deposited is established based on the force conditions and movement process of suspended particles inside saturated porous media. Compared with most current models based on medium particle blockage phenomena (indirect data), this provides a particle deposition calculation method based on the particles themselves (direct data).

[0141] By accurately describing the movement process of individual particles at the microscopic level and calculating whether individual suspended particles in porous media are deposited, it is possible to establish a cross-scale mechanistic connection between the particle deposition process in the medium pores at the microscopic level and the particle blockage phenomenon in saturated porous media at the macroscopic level. This effectively proposes a microscopic mechanism for understanding the particle deposition of groundwater recharge, which facilitates the calculation of the deposition state of particles inside saturated porous media through direct data.

[0142] The above analytical algorithms can all be implemented in the corresponding modules of the multi-physics simulation software COMSOL Multiphysics and obtain converged solutions. Through the coupling between the modules in the COMSOL platform, the reproduction and simulation of the occurrence and development process of particle deposition in the pores inside the groundwater recharge test medium can be realized.

[0143] By setting the research area and boundary conditions in the COMSOL platform and inputting the material properties of the water body and suspended particles, the force and movement trajectory of the suspended particles within the time step can be calculated and simulated. The particle movement trajectory and particle deposition state in the research area are shown in the attached figure. Figure 4 , as attached Figure 4 The figure shows the trajectory of the particle movement process and the three-dimensional path schematic diagram of the particle deposition state model. The path simulation plane of the three-dimensional perspective (viewed from the top along the z-axis, the side along the x-axis, and the side along the y-axis) can be further referenced. Figure 5 , Figure 6 and Figure 7 , a close-up of the sedimentation state of some suspended particles is attached Figure 8 The trajectory of some suspended particles is shown in the attached Fig. 9 .

[0144] On the other hand, the present invention further discloses a computer-readable storage medium storing a computer program, wherein when the computer program is executed by a processor, the processor executes the steps of the above method.

[0145] On the other hand, the present invention further discloses a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the above method.

[0146] In another embodiment provided in the present application, a computer program product including instructions is also provided, which, when executed on a computer, enables the computer to execute any analytical calculation method for the deposition state of groundwater recharge particles in the above-mentioned embodiments.

[0147] It is understandable that the system provided by the embodiment of the present invention corresponds to the method provided by the embodiment of the present invention, and the explanation, examples and beneficial effects of the relevant contents can refer to the corresponding parts in the above method.

[0148] The embodiment of the present application also provides an electronic device, including 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.

[0149] Memory, used to store computer programs;

[0150] The processor is used to implement the analytical calculation method of the above-mentioned groundwater recharge particle deposition state when executing the program stored in the memory.

[0151] The communication bus mentioned in the above electronic device may be a Peripheral Component Interconnect (PCI) bus or an Extended Industrial Standard Architecture (EISA) bus, etc. The communication bus may be divided into an address bus, a data bus, a control bus, etc.

[0152] The communication interface is used for communication between the above electronic device and other devices.

[0153] The memory may include a random access memory (RAM) or a non-volatile memory (NVM), such as at least one disk memory. Optionally, the memory may also be at least one storage device located away from the aforementioned processor.

[0154] The above-mentioned processor can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

[0155] It should also be noted that electronic devices also include terminal devices, which can also be called terminals, user equipment, mobile stations, mobile terminals, etc. Terminal devices can be mobile phones, smart TVs, wearable devices, tablet computers, computers with wireless transceiver functions, virtual reality terminal devices, augmented reality terminal devices, wireless terminals in industrial control, wireless terminals in unmanned driving, wireless terminals in remote surgery, wireless terminals in smart grids, wireless terminals in transportation safety, wireless terminals in smart cities, wireless terminals in smart homes, etc. The embodiments of this application do not limit the specific technology and specific device form used by the terminal devices.

[0156] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented using software, it can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function described in the embodiment of the present application is generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions may be transmitted from one website site, computer, server or data center to another website site, computer, server or data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line) or wireless (e.g., infrared, wireless, microwave, etc.) means. The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more available media integrated. The available medium may be a magnetic medium (e.g., a floppy disk, a hard disk, a tape), an optical medium, or a semiconductor medium (e.g., a solid-state hard disk), etc.

[0157] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

[0158] In addition, it should be noted that if the embodiments of the present invention involve directional indications (such as up, down, left, right, front, back...), the directional indications are only used to explain the relative position relationship, movement status, etc. between the components in a certain specific posture. If the specific posture changes, the directional indication will also change accordingly.

[0159] In addition, if there are descriptions involving "first", "second", etc. in the embodiments of the present invention, the descriptions of "first", "second", etc. are only used for descriptive purposes and cannot be understood as indicating or implying their relative importance or implicitly indicating the number of technical features indicated. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In addition, the meaning of "and / or" appearing in the full text includes three parallel schemes. Taking "A and / or B" as an example, it includes scheme A, or scheme B, or schemes that A and B meet at the same time. In addition, in the embodiments of the present invention, "multiple" refers to more than two. In addition, the technical solutions between the various embodiments can be combined with each other, but it must be based on the ability of ordinary technicians in the field to implement. When the combination of technical solutions is contradictory or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

Claims

1. An analytical calculation method for the sedimentation state of groundwater recharge particles, characterized in that: include: S1. By combining the lattice Boltzmann method with the momentum equation, the equation of motion variables of the particles is constructed to obtain the density distribution function of the particles moving with the water flow and the density distribution function of the particles after colliding with the pore wall; S2. An analytical model is constructed based on the density distribution function of particles moving with water flow and the density distribution function of particles after collision with pore walls combined with the Euler equation model; S3. Detect and obtain the particle density distribution function value in the pores inside the actual medium, input the actual particle density distribution function value into the analytical model, and calculate the deposition state of suspended particles in the water body; in: The specific expression of the motion variable equation in step S1 is: in, For particles in a specified medium The density distribution function in is For particles in a specified medium Single pore inside The density distribution function in is For particles in a specified medium The forces acting on the inside is the spatial position of the particle, is the time node where the particle is located, is the time that the particle is pulled by the force The step length, For particles in a specified medium The speed of movement within a specified time step, is the dimensionless relaxation time; The density distribution function of the particles moving with the water flow in step S1 and the density distribution function of the particles after colliding with the pore wall specifically include: (1) Particles move in pores with water flow: (2) Particles collide with the pore wall along with the water flow: in, For the specified medium The space where the collision process occurs The step length, The time when the collision occurred The step length, is the density distribution function after the particles collide with the pore wall; The specific construction process of the analytical model in step S2 includes: and Substituting the solid-liquid two-phase flow continuity equation in the Euler-Euler model, the particle concentration at different heights of the sand column during the reinjection process can be calculated as the deposition state of suspended particles in the water body. The specific expression of the solid-liquid two-phase flow continuity equation is: in, represents the liquid or solid phase, for The volume concentration of the phase, for The velocity vector of the phase, for The density of the phase, For time, is the Hamiltonian differential operator.

2. An analytical calculation device for the sedimentation state of groundwater recharge particles, characterized in that: An analytical model is stored internally, and the analytical model is constructed based on the analytical calculation method for the sedimentation state of groundwater reinjection particles described in claim 1. The analytical model is also used to input material properties of water bodies and suspended particles to calculate the movement process of suspended particles.

3. The analytical calculation device for groundwater recharge particle deposition state according to claim 2, characterized in that: The specific calculation process of the suspended particle movement process includes: L1. Simulate the flow field of water in the pores inside the analytical medium, which includes the flow velocity and the interaction force on each spatial particle of the water flow; L2. Calculate the total vector force on the particles in the pores of the sand column based on the water flow field. The total vector force includes the pore water flow drag force and the electric field force. L3. Based on the total vector force on the particles, the particle acceleration of the entire time process is obtained as the movement process of the particles in the pores inside the sand column with the water flow.

4. The analytical calculation device for groundwater recharge particle deposition state according to claim 3, characterized in that: The specific operation process in the L1 step includes: L11. Calculation of Reynolds number , the specific calculation formula is: in, is the flow velocity of the seepage inside the pores, is the pore diameter, is the kinematic viscosity of water; L12. Determine the flow pattern of water inside the pores of the porous medium as Stokes flow with a Reynolds number less than 1, and set the microscopic pore water flow field in the pores inside the medium as Stokes flow with neglected inertial force, which is: in, Specify the velocity vector at the particle for the fluid, Specify the pressure vector at the particle for the fluid, is the Hamiltonian differential operator, , , are the three coordinate vectors for the velocity vector in the three-dimensional coordinate system; L13. Based on the Stokes equation to characterize the water flow velocity at each spatial position in the pores inside the medium in the groundwater recharge project, there exists: in, is the water body space step length, is the space step length The water flow vector velocity inside is the space step The vector force acting on the internal water flow is the sum of the external forces.

5. The analytical calculation device for groundwater recharge particle deposition state according to claim 4, characterized in that: The specific calculation expression of the water flow drag force in the L2 step is: in, is the drag force on the particle by the water flow, is the mass of the particle, is the density of the particle material, is the diameter of the particle, is the vector velocity of the specified particle.

6. The analytical calculation device for groundwater recharge particle deposition state according to claim 3, characterized in that: The electric field force in the L2 step is calculated based on Coulomb's law.

7. The analytical calculation device for groundwater recharge particle deposition state according to claim 3, characterized in that: The particle acceleration during the entire time course of the L3 step The specific calculation expression is: in, is the mass of the particle, is the vector velocity of the specified particle, For particles at a specified time step The total vector force on Including the water flow drag force of particles, gravity, water buoyancy and electric field force.

Images