Soluble pollutant diffusion and deposition fluid simulation prediction method and device and server

Abstract

The invention provides a soluble pollutant diffusion deposition fluid simulation prediction method and device and a server, and relates to the technical field of environmental protection simulation, and the method comprises the steps: obtaining topographic parameter information, meteorological parameter information and pollutant emission parameter information of a to-be-detected region, and determining a discrete domain according to the topographic parameter information; determining boundary conditions of a discrete domain according to the wind regime information, and determining diffusion distribution information of pollutants according to the pollutant emission parameter information and the boundary conditions; calculating pollutant dissolution information and pollutant deposition information in a discrete domain based on rainfall information and diffusion distribution information of pollutants, updating a preset diffusion equation so as to perform iterative processing in a next discrete time step, and when iteration is performed to preset maximum simulation time, calculating the pollutant dissolution information and the pollutant deposition information in the discrete domain; and determining target pollutant concentration distribution information and target pollutant surface deposition information. According to the method, the prediction accuracy of the diffusion and deposition conditions of the pollutants can be remarkably improved.

Description

Technical Field

[0001] The present invention relates to the technical field of environmental simulation, and in particular to a method, device and server for simulating and predicting the diffusion and deposition of soluble pollutants in a fluid. Background Art

[0002] With the advancement of urbanization and industrialization, human industrial production activities have generated a large amount of waste pollutants entering the atmospheric environment, and the pressure on the ecosystem has become increasingly heavy. At present, relevant technologies have proposed that the Gaussian plume diffusion theory model can be used to simulate the emission process and diffusion law of point, line or volume source pollutants, so as to predict the distribution of pollutants. However, the existing Gaussian diffusion model uses many assumptions, such as wind flow stability, wind speed uniformity, and uniform and continuous pollution sources, which often leads to large errors with the actual prediction effect and requires a lot of corrections. In addition, in rainy weather, for some highly water-soluble pollutants, their diffusion distribution will be seriously affected, and the existing commonly used Gaussian diffusion model is difficult to analyze this. At present, there are few studies and methods to analyze the diffusion and deposition of pollutants in rainy climates. Summary of the invention

[0003] In view of this, an object of the present invention is to provide a method, device and server for simulating and predicting the diffusion and deposition of soluble pollutants in fluids, which can significantly improve the accuracy of predicting the diffusion and deposition of pollutants.

[0004] In the first aspect, an embodiment of the present invention provides a simulation prediction method for diffusion and deposition of soluble pollutants, the method comprising: obtaining terrain parameter information, meteorological parameter information and pollutant emission parameter information of a test area, and determining a region to be analyzed based on the terrain parameter information, so as to discretize the region to be analyzed and determine a discrete domain, wherein the meteorological parameter information comprises wind condition information and rainfall information; determining boundary conditions of the discrete domain based on the wind condition information, and simulating and solving a preset flow equation and a preset diffusion equation based on the pollutant emission parameter information and the boundary conditions to determine the diffusion distribution information of the pollutants, wherein the boundary conditions comprise velocity magnitude, direction, height distribution and surface roughness of the discrete domain boundary, and the preset flow equation comprises a continuity equation and a Navier-Stokes equation; calculating and processing pollutant dissolution information and pollutant deposition information in the discrete domain based on rainfall information and the diffusion distribution information of pollutants, and updating the preset diffusion equation using the pollutant dissolution information, so as to perform iterative processing within the next discrete time step, and determining the target pollutant concentration distribution information corresponding to each time step and the target pollutant surface deposition information at each position when iterating to the preset maximum simulation time.

[0005] In one embodiment, the step of discretizing the area to be analyzed and determining the discrete domain includes: obtaining the total length of the horizontal plane, the total width of the horizontal plane and the terrain elevation parameter information in the area to be analyzed, and using the total length of the horizontal plane and the total width of the horizontal plane to discretize the area to be analyzed to determine a set of two-dimensional discrete Cartesian coordinate points; using a preset linear difference model, based on the terrain elevation parameter information, linear difference calculation processing is performed on the two-dimensional discrete Cartesian coordinate points to determine the elevation information corresponding to each discrete point in the two-dimensional discrete Cartesian coordinate points; using the elevation information, mapping the two-dimensional discrete Cartesian coordinate points in the height direction to determine a set of three-dimensional discrete Cartesian coordinate points, and determining a three-dimensional geometric discrete domain with actual terrain information based on the three-dimensional discrete Cartesian coordinate points.

[0006] In one embodiment, the step of determining the boundary conditions of the discrete domain according to the wind condition information includes: determining the change function of the wind speed with the height above the ground at the boundary of the discrete domain based on the wind condition information and the wind profile theory, so as to determine the boundary conditions by using the change function; wherein the change function is expressed as:

[0007] in, To solve the wind speed on the boundary of the domain, the wind speed and time and height Related, is the reference height, is the initial wind speed at the reference height, and a is the power exponent that affects wind distribution.

[0008] In one embodiment, based on rainfall information and diffusion distribution information of pollutants, the pollutant dissolution information and pollutant deposition information in the discrete domain are calculated and processed, and the preset diffusion equation is updated using the pollutant dissolution information, including: based on rainfall information, for each discrete unit at the position of each discrete unit in the horizontal plane of the discrete domain, the amount of soluble pollutants absorbed by raindrops in a single discrete time step is calculated and processed to determine the pollutant dissolution information, and the amount of pollutants absorbed by raindrops is accumulated and calculated to determine the pollutant deposition information; the absorption amount is used as the absorption source term to update the preset diffusion equation for iterative processing within the next discrete time step.

[0009] In one embodiment, based on rainfall information, the absorption amount of soluble pollutants by raindrops at a single discrete time step is calculated and processed for each discrete unit at the location of each discrete unit in the horizontal plane of the discrete domain, and the step of determining the pollutant dissolution information includes: converting the rainfall information into precipitation at each discrete time step, and using the precipitation to determine the number of raindrops and the volume swept by the raindrops in the space of the single unit passed by; calculating the total mass of pollutants absorbed by the raindrops at each discrete time based on the volume and pollutant concentration to determine the absorption amount of pollutants by raindrops in each unit, and determining the absorption amount as the pollutant dissolution information.

[0010] In one embodiment, before the step of converting the rainfall information into the precipitation at each discrete time step, the method includes: using a preset time step number calculation model to calculate the time required for the raindrops generated at each discrete time step to fall to the ground, determining the number of retention time steps of the raindrops generated at each time step in the discrete domain, and converting the rainfall information into the precipitation at each discrete time step according to the retention time steps, wherein the preset time step number calculation model is expressed as:

[0011] in, is the number of time steps to be retained, is the ceiling function, is the total elevation of the discrete unit on the horizontal plane, is a discrete time step, is the average radius of raindrops, is the resistance coefficient.

[0012] In one embodiment, the step of accumulating and calculating the amount of pollutants absorbed by raindrops to determine the pollutant deposition information includes: accumulating and calculating the amount of pollutants absorbed by raindrops based on the number of retention time steps through a pollutant deposition mass calculation model to determine the pollutant deposition information; wherein the pollutant deposition mass calculation model is expressed as:

[0013] in, is the pollutant deposition mass, is the total mass of pollutants absorbed by raindrops at each discrete time.

[0014] In a second aspect, an embodiment of the present invention further provides a simulation prediction device for diffusion and deposition of soluble pollutants, the device comprising: a region division module, which obtains terrain parameter information, meteorological parameter information and pollutant emission parameter information of the area to be measured, and determines the area to be analyzed according to the terrain parameter information, so as to discretize the area to be analyzed and determine the discrete domain, wherein the meteorological parameter information includes: wind condition information and rainfall information; a diffusion analysis module, which determines the boundary conditions of the discrete domain according to the wind condition information, and simulates and solves the preset flow equation and the preset diffusion equation according to the pollutant emission parameter information and the boundary conditions to determine the diffusion distribution of the pollutants. The preset flow equations include: the continuity equation and the Navier-Stokes equation; the deposition prediction module calculates and processes the pollutant dissolution information and the pollutant deposition information in the discrete domain based on the rainfall information and the diffusion distribution information of the pollutants, and uses the pollutant dissolution information to update the preset diffusion equation, so as to iterate within the next discrete time step, and determine the target pollutant concentration distribution information corresponding to each time step and the target pollutant surface deposition information at each position when iterating to the preset maximum simulation time.

[0015] In a third aspect, an embodiment of the present invention further provides a server, comprising a processor and a memory, wherein the memory stores computer executable instructions that can be executed by the processor, and the processor executes the computer executable instructions to implement any one of the methods provided in the first aspect.

[0016] In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium, which stores computer-executable instructions. When the computer-executable instructions are called and executed by a processor, the computer-executable instructions prompt the processor to implement any one of the methods provided in the first aspect.

[0017] The embodiments of the present invention bring the following beneficial effects: The embodiment of the present invention provides a method, device and server for simulating and predicting diffusion and deposition fluid of soluble pollutants. The method obtains terrain parameter information, meteorological parameter information and pollutant emission parameter information of the area to be measured, and determines the area to be analyzed according to the terrain parameter information, so as to discretize the area to be analyzed. After determining the discrete domain, the boundary conditions of the discrete domain are determined according to the wind condition information, and the preset flow equation and the preset diffusion equation are simulated and solved according to the pollutant emission parameter information and the boundary conditions to determine the diffusion distribution information of the pollutants. Finally, based on the rainfall information and the diffusion distribution information of the pollutants, the pollutant dissolution information and the pollutant deposition information in the discrete domain are calculated and processed, and the preset diffusion equation is updated using the pollutant dissolution information to perform iterative processing within the next discrete time step, and when iterating to the preset maximum simulation time, the target pollutant concentration distribution information corresponding to each time step and the target pollutant surface deposition information at each position are determined. The embodiment of the present invention can build a three-dimensional atmospheric flow field based on the terrain of the plant area to be analyzed, use meteorological data, perform large-scale calculations on complex atmospheric flow conditions and the diffusion of pollutants, and consider the influence of rainwater on the diffusion of soluble pollutants to accurately analyze the pollutant deposition.

[0018] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention are realized and obtained by the structures particularly pointed out in the description, claims and drawings.

[0019] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0021] Figure 1 A schematic flow chart of a method for simulating and predicting the diffusion and deposition of soluble pollutants in a fluid according to an embodiment of the present invention; Figure 2 A schematic diagram of a specific process of a method for simulating and predicting the diffusion and deposition of soluble pollutants in a fluid according to an embodiment of the present invention; Figure 3 A schematic diagram of the structure of a device for simulating and predicting the diffusion and deposition of soluble pollutants in a fluid provided by an embodiment of the present invention; Figure 4 A schematic diagram of the structure of a server provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0022] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution of the present invention will be clearly and completely described in combination with the embodiments below. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0023] At present, with the advancement of urbanization and industrialization, human industrial production activities have generated a large amount of waste pollutants entering the atmospheric environment, and the pressure on the ecosystem has become increasingly heavy. Common atmospheric pollutants emitted from industrial plants are VOCs, DMF, benzene, toluene, formaldehyde, etc., which have an unpleasant odor and are toxic, irritating and carcinogenic, causing great harm to the human body. Among them, some soluble pollutants, such as formaldehyde and DMF, may enter the urban water system with rainfall, affecting the safety of residents' lives. For this reason, a method for predicting the diffusion and deposition of atmospheric pollutants emitted by enterprises is proposed, which has important practical significance for the detection and control of pollutant emissions from polluting enterprises in plant areas.

[0024] In the existing pollutant diffusion prediction technology, the Gaussian plume diffusion theory model is often used to simulate the emission process and diffusion law of point, line or volume source pollutants, so as to predict the distribution of pollutants. However, the existing Gaussian diffusion model adopts many assumptions, such as wind flow stability, wind speed uniformity, and uniform and continuous pollution sources, which often leads to large errors with the actual prediction effect and requires a lot of corrections. In addition, in rainy weather, for some water-soluble pollutants, their diffusion distribution will be seriously affected, and the existing commonly used Gaussian diffusion model is difficult to analyze this. There are also few studies and methods for analyzing the diffusion and sedimentation of pollutants in rainy climates. Based on this, the simulation prediction method, device and server of the diffusion and deposition fluid of soluble pollutants provided by the present invention can be based on the terrain of the plant area to be analyzed, build a three-dimensional atmospheric flow field, use meteorological data, and perform large-scale calculations on complex atmospheric flow conditions and the diffusion of pollutants, and consider the influence of rainwater on the diffusion of soluble pollutants to accurately analyze the sedimentation of pollutants.

[0025] See also Figure 1 The flowchart of a method for simulating and predicting the diffusion and deposition of soluble pollutants in fluid is shown in FIG. 1 , and the method mainly includes the following steps S102 to S106: Step S102, obtaining terrain parameter information, meteorological parameter information and pollutant emission parameter information of the area to be measured, and determining the area to be analyzed based on the terrain parameter information, so as to discretize the area to be analyzed and determine the discrete domain, wherein the meteorological parameter information includes: wind condition information and rainfall information, the terrain parameters include: the total length, total width and total height of the area to be analyzed, as well as the surface elevation information and surface roughness in the area; the wind condition information includes: the incoming wind direction, incoming wind speed and wind measurement reference height that change with time in the time period to be analyzed; the rainfall information includes: the precipitation that changes with time in the time period to be analyzed; the pollutant emission parameters include: the pollutant emission position, pollutant emission concentration and diffusion coefficient of the pollutant in the time period to be analyzed.

[0026] In a real-time manner, the total length X, total width Y and terrain elevation parameter information of the horizontal plane in the area to be analyzed can be obtained, and the total length and total width of the horizontal plane are used to discretize the area to be analyzed to determine a set of two-dimensional discrete Cartesian coordinate points. Then, through the preset linear difference model, based on the terrain elevation parameter information, the linear difference calculation processing is performed on the two-dimensional discrete Cartesian coordinate points to determine the elevation information corresponding to each discrete point in the two-dimensional discrete Cartesian coordinate points Finally, using the elevation information, the two-dimensional discrete Cartesian coordinate points are mapped in the height direction to determine a set of three-dimensional discrete Cartesian coordinate points , and determine the three-dimensional geometric discrete domain with actual terrain information based on the three-dimensional discrete Cartesian coordinate points.

[0027] Step S104, determining the boundary conditions of the discrete domain according to the wind condition information, and performing simulation solution processing on the preset flow equation and the preset diffusion equation according to the pollutant emission parameter information and the boundary conditions to determine the diffusion distribution information of the pollutants, wherein the boundary conditions include: the velocity magnitude, direction, height distribution and surface roughness of the discrete domain boundary, the preset flow equations include: the continuity equation and the Navier-Stokes equations, and the simulation solution refers to replacing the original integral form or partial differential form of the flow equation and the diffusion equation by the algebraic form of time and space discretization and performing iterative solution.

[0028] In one implementation, the wind condition information can be used to determine the variation function of wind speed with height above the ground at the boundary of the discrete domain based on the wind profile theory, so as to determine the boundary condition using the variation function; wherein the variation function is expressed as:

[0029] in, To solve the wind speed on the boundary of the domain, the wind speed and time and height Related, is the reference height, is the initial wind speed at the reference height, and a is the power exponent that affects wind distribution.

[0030] In another embodiment, the diffusion equation calculates each atmospheric pollutant as a custom scalar in order to calculate the process of pollutant diffusion; the diffusion equation makes the following assumptions about atmospheric pollutants: 1. The composition of the pollutant is stable and will not chemically react with other components in the air; 2. The pollutant has no effect on the flow of the entire atmospheric flow field, where: Continuity equation:

[0031] Navier-Stokes equations:

[0032] in, For time, For pressure, is the fluid velocity, is the fluid density, is the Reynolds number, is the source term of the momentum conservation equation.

[0033] Furthermore, since the basic mode of mass transfer of the custom scalar field is convective mass transfer, which is similar to convective heat transfer, the solution to its mass transfer problem is also related to factors such as the flow pattern and velocity distribution of the fluid. At the same time, the heat, momentum and mass transfer processes are affected by the turbulent transfer caused by the macroscopic movement of fluid micro-groups. The diffusion equation is expressed as follows:

[0034]

[0035] in, is the mass average density, is the subscript of the component to be analyzed, and the component is the atmospheric pollutant to be analyzed. For the The density of the component, For the The mass content of the components, is the turbulent viscosity, For the The diffusion coefficient of the component, is the Schmidt number of the fluid, It is The source term for the creation or absorption of a component.

[0036] Among them, The source term for the production or absorption of a component can be further expanded:

[0037] in, The discretized equation The source term of a component is directly determined by the emission location and emission concentration of the component in the pollutant emission parameters; The discretized equation Absorption source term of the component.

[0038] Step S106, based on the rainfall information and the diffusion distribution information of the pollutants, the pollutant dissolution information and the pollutant deposition information in the discrete domain are calculated and processed, and the preset diffusion equation is updated using the pollutant dissolution information, so as to perform iterative processing within the next discrete time step, and when iterating to the preset maximum simulation time, the target pollutant concentration distribution information corresponding to each time step and the target pollutant surface deposition information at each position are determined. In one embodiment, the dissolution and deposition conditions in the discrete domain are statistically analyzed, which means that for each discrete body unit at the position of each discrete unit in the horizontal plane of the calculation domain, the absorption of easily soluble pollutants by raindrops is calculated based on the rainfall information in a single discrete time step, and the absorption of pollutants is used as an absorption source term and applied to the discretized equation.

[0039] The above-mentioned soluble pollutant diffusion and deposition fluid simulation prediction method provided by the embodiment of the present invention can significantly improve the accuracy of the prediction of the diffusion and deposition of pollutants.

[0040] The embodiment of the present invention also provides an implementation method for determining the dissolution and deposition of pollutants, as shown in (A) to (D) below: (A) Furthermore, the calculation of the amount of soluble pollutants absorbed by raindrops requires calculating the time required for the raindrops generated at each discrete time step to fall to the ground, and obtaining the number of time steps of the raindrops generated at each time step in the discrete domain. That is, using a preset time step calculation model, the time required for the raindrops generated at each discrete time step to fall to the ground is calculated and processed, and the number of time steps of the raindrops generated at each time step in the discrete domain is determined, so that the rainfall information is converted into the precipitation at each discrete time step according to the number of time steps. The preset time step calculation model is expressed as:

[0041] in, is the number of time steps to be retained, is the ceiling function, is the total elevation of the discrete unit on the horizontal plane, is a discrete time step, is the average radius of raindrops, is the resistance coefficient.

[0042] (B) Based on rainfall information, the amount of soluble pollutants absorbed by raindrops at a single discrete time step is calculated for each discrete unit at the location of each discrete unit in the horizontal plane of the discrete domain to determine pollutant dissolution information. In one embodiment, the rainfall information is converted into precipitation at each discrete time step, and the precipitation is used to determine the number of raindrops and the volume swept by the raindrops in the space of the single unit. Then, based on the volume and pollutant concentration, the total mass of pollutants absorbed by the raindrops at each discrete time is calculated to determine the amount of pollutants absorbed by the raindrops in each unit, and the absorption amount is determined as the pollutant dissolution information. Specifically, it includes the following: (1) to (4): (1) Convert the precipitation into the precipitation at each discrete time step, that is:

[0043] in, is the rainfall at each discrete time step, It is the amount of rainfall in a certain period of time within the meteorological information. is the rainfall frequency represented by this time period, is the discrete time step.

[0044] (2) At each discrete time step, all raindrops in the discrete region are in space, and the amount of pollutants absorbed by them by sweeping the discrete space is calculated. To do this, it is necessary to first calculate the number of raindrops. Calculate, the relationship is as follows:

[0045] in, is the area of ​​each discrete surface element in the horizontal plane at the top of the computational domain.

[0046] (3) Calculate all raindrops in a discrete time step The volume swept through the space of a single volume unit is , its main relationships are as follows:

[0047] in, is the maximum cross-sectional area of ​​a raindrop, It is the elevation swept by the raindrop in the volume unit.

[0048] (4) Calculate and accumulate the mass of pollutants absorbed by all volume units passed by the raindrops at the discrete time step to obtain the total mass of pollutants absorbed by the rainfall at the time step. The relationship is as follows:

[0049] in, To obtain the pollutant concentration within the volume unit at this discrete time step during the solution process.

[0050] (C) Accumulate and calculate the amount of pollutants absorbed by raindrops to determine the pollutant deposition information, use the absorbed amount as the absorption source term, and update the preset diffusion equation to iteratively process in the next discrete time step. Specifically, the pollutant deposition mass calculation model can be used to accumulate the amount of pollutants absorbed by raindrops based on the number of retention time steps to determine the pollutant deposition information. In other words, for the calculation of pollutant deposition in each area, it is necessary to accumulate all the pollutants in this time step and the previous time step. The mass of pollutants absorbed by raindrops generated in a time step and not yet falling to the ground, that is, the mass of pollutant deposition in the discrete area at this time step , where the pollutant deposition mass calculation model is expressed as:

[0051] in, is the pollutant deposition mass, is the total mass of pollutants absorbed by raindrops at each discrete time.

[0052] (D) Finally, the absorption amount is used as the absorption source term to update the preset diffusion equation so that it can be iterated in the next discrete time step, and the mass of pollutants absorbed by rainwater in each volume unit is reapplied as the absorption source term to the discrete equation:

[0053] Through the above simulation solution, the pollutant deposition concentration of each discrete domain in the discretized calculation domain at each time step can be obtained, and the impact of rainfall on pollutant diffusion can be updated as a source term in the discrete equation for the calculation of the next time step.

[0054] See also Figure 2The specific flow chart of a simulation prediction method for diffusion and deposition of soluble pollutants in fluid is shown. The present invention can accurately calculate the diffusion of pollutants and the deposition of soluble pollutants with rainfall in the complex atmospheric flow field in the region under natural wind conditions and after considering the effect of rainfall. The present invention includes the following four steps: 1. Based on the terrain parameter information, define the area to be analyzed and discretize it into discrete domains; 2. Based on the wind condition information in the meteorological parameter information, define the velocity, direction and height distribution of the discrete domain boundary; 3. Based on the pollutant emission parameter information, simulate and solve the flow equation and the diffusion equation to calculate the diffusion distribution of pollutants; 4. Based on the rainfall information in the meteorological parameter information, calculate the absorption of soluble pollutants in the discrete domain and the deposition of pollutants at various locations on the surface step by time; starting from the initial time of the simulation, in each discrete time step, repeat the above process 2-4 until the maximum simulation time is reached, so as to obtain the transient distribution of pollutant concentrations at various locations in the area to be analyzed, and obtain the surface deposition of pollutants at various locations.

[0055] In actual application, based on the production conditions and meteorological conditions of a certain economic development zone, three-day emission calculations of production enterprises were carried out, and the diffusion and precipitation deposition of a certain pollutant were analyzed. The range to be analyzed in the economic development zone is 10km×10km, and the height range is 300m; the production conditions to be analyzed include the location information, height information, outlet inner diameter, exhaust volume and emission gas concentration of the exhaust pipes of 27 emission enterprises and 142 production lines in the economic development zone; the meteorological conditions used include the local meteorological monitoring information every 3 hours within 3 days, including wind direction, wind speed and rainfall information: Step 1: Discretize the economic development area to be analyzed based on its total length of 10km and total width of 10km into a set of two-dimensional discrete Cartesian coordinate points Based on the elevation parameter information of the area, the elevation of each two-dimensional discrete point is obtained through linear interpolation. It is best to map the total height of 300m based on the elevation parameter information of each discrete point in the height direction, and finally obtain a set of three-dimensional discrete Cartesian coordinate points. , thereby obtaining the three-dimensional geometric discrete domain of the area.

[0056] Step 2: Based on the wind speed and wind direction information that changes with time in the meteorological parameters, the wind profile theory is used to solve the wind speed at the boundary of the domain with the height above the ground. The change function is defined as follows:

[0057] in, To solve the wind speed on the boundary of the domain, its relationship with time and height Related, is the reference height, is the initial wind speed at the reference height, and a is the power exponent that affects wind distribution.

[0058] Based on the relationship between the wind speed and the height z above the ground, the wind speed and direction at a specific height are directly defined on the boundary of the three-dimensional discrete domain: Step 3: This specific embodiment only involves the diffusion analysis of one type of pollutant gas, and it is necessary to define the source term of the pollutant gas in each discrete unit in the solution domain:

[0059] in, The coordinates of a center are The source term of pollutant gas generation in discrete units, It is the subscript of the production line number of a certain production area. is the pollutant gas emission concentration of the exhaust pipe of the production line, is the location of the exhaust pipe of the production line. is the radius of the exhaust chimney of the production line.

[0060] Step 4: Based on the source term defined in step 3, the pollutant absorption source term of each discrete unit in the solution domain can be considered to be 0, and the first time step is solved. Based on the pollutant concentration of the discrete unit in the solution result, for each discrete unit in the horizontal plane of the calculation domain, the absorption of pollutants by raindrops is calculated based on rainfall information in a single discrete time step, and the absorption source term is applied to the discretized equation for subsequent iterative solution.

[0061] Specifically, by calculating the time required for the raindrops generated at each time step to fall to the ground, the number of time steps that the raindrops generated at each time step remain in the discrete domain is obtained. ,Right now:

[0062] in, is the ceiling function, is the total elevation of the discrete unit on the horizontal plane, is a discrete time step, is the average radius of raindrops, is the resistance coefficient.

[0063] The rainfall information in the above single discrete time step refers to the need to convert the precipitation into the precipitation in each discrete time step, that is:

[0064] in, is the rainfall at each discrete time step, It is the amount of rainfall in a certain period of time within the meteorological information. is the rainfall frequency represented by this time period, is the discrete time step.

[0065] The calculation of the amount of pollutants absorbed by raindrops mentioned above refers to calculating the amount of pollutants absorbed by all raindrops in the discrete area by sweeping the discrete space at each discrete time step. To this end, it is necessary to first calculate the number of raindrops. Calculate, the relationship is as follows:

[0066] in, is the area of ​​each discrete surface element in the horizontal plane at the top of the computational domain.

[0067] Second, we need to calculate all raindrops in a discrete time step The volume swept through the space of a single volume unit is , its main relationships are as follows:

[0068] in, is the maximum cross-sectional area of ​​a raindrop, It is the elevation swept by the raindrop within the volume unit.

[0069] Then, the mass of pollutants absorbed by sweeping in all volume units passed by the raindrops at the discrete time step is calculated and accumulated to obtain the total mass of pollutants absorbed by the rainfall at the time step. The relationship is as follows:

[0070] in, To obtain the pollutant concentration within the volume unit at this discrete time step during the solution process.

[0071] Finally, accumulate the time step and the previous The mass of pollutants absorbed by raindrops that have not yet fallen to the ground is generated in a time step, that is, the mass of pollutant deposition in the discrete area at this time step , the relationship is as follows:

[0072] The aforementioned use of the pollutant absorption as an absorption source term to act on the discretized equation means that the mass of pollutants absorbed by rainwater in each volume unit is reapplied as an absorption source term to the discrete equation:

[0073] Through the above simulation solution, the pollutant deposition concentration of each discrete domain in the discretized calculation domain at each time step can be obtained, and the impact of rainfall on pollutant diffusion can be updated as a source term in the discrete equation for the calculation of the next time step until the maximum simulation time is reached, that is, 3 days.

[0074] In summary, compared with the commonly used theoretical models, the analysis method disclosed in the present invention has the ability to accurately calculate the diffusion of pollutants at any position in three dimensions and on a large scale, and can take into account rainfall conditions and the surface deposition of pollutants, and combined with the local plant topography and meteorological information, it can provide theoretical support for environmental governance-related departments to monitor, control and deal with the production activities of pollutant-emitting plants and their impacts. In addition, the present invention has a wide range of applicable working conditions and strong scalability. It can be combined with the local meteorological information system to conduct scientific analysis and prediction of pollutant diffusion and deposition, and assist in tracing and managing pollution sources.

[0075] For the simulation prediction method of the soluble pollutant diffusion deposition fluid provided in the above embodiment, the embodiment of the present invention provides a simulation prediction device for the soluble pollutant diffusion deposition fluid, see Figure 3 The schematic diagram of the structure of a soluble pollutant diffusion deposition fluid simulation prediction device is shown, and the device includes the following parts: The region division module 302 obtains terrain parameter information, meteorological parameter information and pollutant emission parameter information of the region to be measured, and determines the region to be analyzed according to the terrain parameter information, so as to discretize the region to be analyzed and determine a discrete domain, wherein the meteorological parameter information includes: wind condition information and rainfall information; The diffusion analysis module 304 determines the boundary conditions of the discrete domain according to the wind condition information, and performs simulation and solution processing on the preset flow equation and the preset diffusion equation according to the pollutant emission parameter information and the boundary conditions to determine the diffusion distribution information of the pollutants, wherein the boundary conditions include: the velocity magnitude, direction, height distribution and surface roughness of the discrete domain boundary, and the preset flow equation includes: the continuity equation and the Navier-Stokes equation; The deposition prediction module 306 calculates and processes the pollutant dissolution information and pollutant deposition information in the discrete domain based on rainfall information and the diffusion distribution information of pollutants, and uses the pollutant dissolution information to update the preset diffusion equation so as to perform iterative processing within the next discrete time step, and when iterating to the preset maximum simulation time, determine the target pollutant concentration distribution information corresponding to each time step and the target pollutant surface deposition information at each position.

[0076] The above-mentioned soluble pollutant diffusion and deposition fluid simulation prediction device provided in the embodiment of the present application can significantly improve the accuracy of the prediction of the diffusion and deposition of pollutants.

[0077] In one embodiment, when discretizing the area to be analyzed and determining the discrete domain, the above-mentioned area division module 302 is also used to: obtain the total length of the horizontal plane, the total width of the horizontal plane and the terrain elevation parameter information in the area to be analyzed, and use the total length of the horizontal plane and the total width of the horizontal plane to discretize the area to be analyzed to determine a set of two-dimensional discrete Cartesian coordinate points; through a preset linear difference model, based on the terrain elevation parameter information, linear difference calculation processing is performed on the two-dimensional discrete Cartesian coordinate points to determine the elevation information corresponding to each discrete point in the two-dimensional discrete Cartesian coordinate points; using the elevation information, the two-dimensional discrete Cartesian coordinate points are mapped in the height direction to determine a set of three-dimensional discrete Cartesian coordinate points, and a three-dimensional geometric discrete domain with actual terrain information is determined based on the three-dimensional discrete Cartesian coordinate points.

[0078] In one implementation, when performing the step of determining the boundary conditions of the discrete domain according to the wind condition information, the diffusion analysis module 304 is further used to: determine the change function of the wind speed with the height above the ground at the boundary of the discrete domain through the wind condition information based on the wind profile theory, so as to determine the boundary conditions using the change function; wherein the change function is expressed as:

[0079] in, To solve the wind speed on the boundary of the domain, the wind speed and time and height Related, is the reference height, is the initial wind speed at the reference height, and a is the power exponent that affects wind distribution.

[0080] In one embodiment, when calculating and processing pollutant dissolution information and pollutant deposition information in a discrete domain based on rainfall information and diffusion distribution information of pollutants, and updating a preset diffusion equation using the pollutant dissolution information, the above-mentioned deposition prediction module 306 is also used to: calculate and process the absorption amount of soluble pollutants by raindrops in a single discrete time step for each discrete unit at the location of each discrete unit in the horizontal plane of the discrete domain based on rainfall information, determine pollutant dissolution information, and accumulate and calculate the amount of pollutants absorbed by raindrops to determine pollutant deposition information; use the absorption amount as the absorption source term, and update the preset diffusion equation for iterative processing within the next discrete time step.

[0081] In one embodiment, when calculating and processing the absorption of soluble pollutants by raindrops at a single discrete time step for each discrete unit at the location of each discrete unit in the horizontal plane of the discrete domain based on rainfall information to determine the pollutant dissolution information, the above-mentioned deposition prediction module 306 is also used to: convert the rainfall information into precipitation at each discrete time step, and use the precipitation to determine the number of raindrops and the volume swept by the raindrops in the space of the single unit passed by; calculate the total mass of pollutants absorbed by the raindrops at each discrete time based on the volume and pollutant concentration to determine the absorption of pollutants by raindrops in each unit, and determine the absorption as the pollutant dissolution information.

[0082] In one embodiment, before converting the rainfall information into the precipitation at each discrete time step, the deposition prediction module 306 is further used to: use a preset time step number calculation model to calculate the time required for the raindrops generated at each discrete time step to fall to the ground, determine the number of retention time steps of the raindrops generated at each time step in the discrete domain, and convert the rainfall information into the precipitation at each discrete time step according to the retention time steps, wherein the preset time step number calculation model is expressed as:

[0083] in, is the number of time steps to be retained, is the ceiling function, is the total elevation of the discrete unit on the horizontal plane, is a discrete time step, is the average radius of raindrops, is the resistance coefficient.

[0084] In one embodiment, when performing the step of accumulating and calculating the amount of pollutants absorbed by raindrops to determine the pollutant deposition information, the above-mentioned deposition prediction module 306 is further used to: perform accumulative calculation of the amount of pollutants absorbed by raindrops based on the number of retention time steps through a pollutant deposition mass calculation model to determine the pollutant deposition information; wherein the pollutant deposition mass calculation model is expressed as:

[0085] in, is the pollutant deposition mass, is the total mass of pollutants absorbed by raindrops at each discrete time.

[0086] The device provided in the embodiment of the present invention has the same implementation principle and technical effects as those of the aforementioned method embodiment. For the sake of brief description, for matters not mentioned in the device embodiment, reference may be made to the corresponding contents in the aforementioned method embodiment.

[0087] An embodiment of the present invention provides a server. Specifically, the server includes a processor and a storage device. The storage device stores a computer program, and when the computer program is executed by the processor, it executes the method described in any one of the above-mentioned embodiments.

[0088] Figure 4 A structural diagram of a server provided in an embodiment of the present invention, the server 100 includes: a processor 40, a memory 41, a bus 42 and a communication interface 43, wherein the processor 40, the communication interface 43 and the memory 41 are connected via the bus 42; the processor 40 is used to execute an executable module stored in the memory 41, such as a computer program.

[0089] The memory 41 may include a high-speed random access memory (RAM), and may also include a non-volatile memory, such as at least one disk memory. The communication connection between the system network element and at least one other network element is realized through at least one communication interface 43 (which may be wired or wireless), and the Internet, wide area network, local area network, metropolitan area network, etc. may be used.

[0090] The bus 42 may be an ISA bus, a PCI bus, or an EISA bus, etc. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 4 Only one bidirectional arrow is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0091] Among them, the memory 41 is used to store programs, and the processor 40 executes the program after receiving the execution instruction. The method executed by the device for flow process definition disclosed in any embodiment of the above-mentioned embodiment of the present invention can be applied to the processor 40 or implemented by the processor 40.

[0092] The processor 40 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by the hardware integrated logic circuit or software instructions in the processor 40. The above processor 40 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 gates or transistor logic devices, discrete hardware components. The methods, steps and logic block diagrams disclosed in the embodiments of the present invention can be implemented or executed. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor. The steps of the method disclosed in the embodiment of the present invention can be directly embodied as a hardware decoding processor to be executed, or the hardware and software modules in the decoding processor can be executed. The software module can be located in a mature storage medium in the field such as a random access memory, a flash memory, a read-only memory, a programmable read-only memory, or an electrically erasable programmable memory, a register, etc. The storage medium is located in the memory 41, and the processor 40 reads the information in the memory 41 and completes the steps of the above method in combination with its hardware.

[0093] The computer program product of the readable storage medium provided in the embodiment of the present invention includes a computer-readable storage medium storing program code, and the instructions included in the program code can be used to execute the methods described in the previous method embodiments. The specific implementation can be referred to the previous method embodiments, which will not be repeated here.

[0094] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc., which can store program codes.

[0095] Finally, it should be noted that the above-described embodiments are only specific implementations of the present invention, which are used to illustrate the technical solutions of the present invention, rather than to limit them. The protection scope of the present invention is not limited thereto. Although the present invention is described in detail with reference to the above-described embodiments, ordinary technicians in the field should understand that any technician familiar with the technical field can still modify the technical solutions recorded in the above-described embodiments within the technical scope disclosed by the present invention, or can easily think of changes, or make equivalent replacements for some of the technical features therein; and these modifications, changes or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

1. A simulation prediction method for diffusion and deposition of soluble pollutants in fluids, characterized in that: The method comprises: Acquire terrain parameter information, meteorological parameter information and pollutant emission parameter information of the area to be measured, and determine the area to be analyzed according to the terrain parameter information, so as to discretize the area to be analyzed and determine a discrete domain, wherein the meteorological parameter information includes: wind condition information and rainfall information; Determine the boundary conditions of the discrete domain according to the wind condition information, and perform simulation and solution processing on the preset flow equation and the preset diffusion equation according to the pollutant emission parameter information and the boundary conditions to determine the diffusion distribution information of the pollutants, wherein the boundary conditions include: the velocity magnitude, direction, height distribution and surface roughness of the discrete domain boundary, and the preset flow equation includes: the continuity equation and the Navier-Stokes equation; Based on the rainfall information and the diffusion distribution information of the pollutants, the pollutant dissolution information and the pollutant deposition information in the discrete domain are calculated and processed, and the preset diffusion equation is updated using the pollutant dissolution information, so as to perform iterative processing within the next discrete time step, and when iterating to the preset maximum simulation time, the target pollutant concentration distribution information corresponding to each time step and the target pollutant surface deposition information at each position are determined.

2. The method for simulating and predicting the diffusion and deposition of soluble pollutants according to claim 1, characterized in that: The step of discretizing the area to be analyzed to determine a discrete domain includes: Obtaining the total length of the horizontal plane, the total width of the horizontal plane and the terrain elevation parameter information in the area to be analyzed, and using the total length of the horizontal plane and the total width of the horizontal plane to discretize the area to be analyzed to determine a set of two-dimensional discrete Cartesian coordinate points; By using a preset linear difference model, based on the terrain elevation parameter information, a linear difference calculation process is performed on the two-dimensional discrete Cartesian coordinate points to determine the elevation information corresponding to each discrete point in the two-dimensional discrete Cartesian coordinate points; The elevation information is used to map the two-dimensional discrete Cartesian coordinate points in the height direction to determine a set of three-dimensional discrete Cartesian coordinate points, and a three-dimensional geometric discrete domain with actual terrain information is determined based on the three-dimensional discrete Cartesian coordinate points.

3. The method for simulating and predicting the diffusion and deposition of soluble pollutants according to claim 1, characterized in that: The step of determining the boundary conditions of the discrete domain according to the wind condition information comprises: By using the wind condition information and based on the wind profile theory, a variation function of wind speed with height above the ground at the boundary of the discrete domain is determined, so as to determine the boundary condition by using the variation function; wherein the variation function is expressed as: in, To solve the wind speed on the boundary of the domain, the wind speed and time and height Related, is the reference height, is the initial wind speed at the reference height, and a is the power exponent that affects wind distribution.

4. The method for simulating and predicting the diffusion and deposition of soluble pollutants according to claim 1, characterized in that: The step of calculating and processing the pollutant dissolution information and the pollutant deposition information in the discrete domain based on the rainfall information and the diffusion distribution information of the pollutants, and updating the preset diffusion equation using the pollutant dissolution information includes: Based on the rainfall information, for each discrete unit at the location of each discrete unit in the horizontal plane of the discrete domain, the amount of soluble pollutants absorbed by raindrops in a single discrete time step is calculated to determine the pollutant dissolution information, and the amount of pollutants absorbed by raindrops is accumulated to determine the pollutant deposition information; The absorption amount is used as an absorption source term to update the preset diffusion equation so as to perform iterative processing within the next discrete time step.

5. The method for simulating and predicting the diffusion and deposition of soluble pollutants according to claim 4, characterized in that: The step of calculating and processing the absorption amount of easily soluble pollutants by raindrops at a single discrete time step for each discrete unit at the location of each discrete unit in the horizontal plane of the discrete domain based on the rainfall information to determine the pollutant dissolution information includes: Converting the rainfall information into precipitation at each discrete time step, and using the precipitation to determine the number of raindrops and the volume swept by the raindrops in the space of a single body unit; The total mass of pollutants absorbed by raindrops in each discrete time is calculated according to the volume and pollutant concentration to determine the amount of pollutants absorbed by raindrops in each unit cell, and the absorption amount is determined as the pollutant dissolution information.

6. The method for simulating and predicting the diffusion and deposition of soluble pollutants according to claim 5, characterized in that: Before the step of converting the rainfall information into the precipitation amount at each discrete time step, the method includes: The preset time step calculation model is used to calculate the time required for the raindrops generated in each discrete time step to fall to the ground, and the retention time steps of the raindrops generated in each time step in the discrete domain are determined, so as to convert the rainfall information into the precipitation in each discrete time step according to the retention time steps, wherein the preset time step calculation model is expressed as: in, is the number of time steps to be retained, is the ceiling function, is the total elevation of the discrete unit on the horizontal plane, is a discrete time step, is the average radius of raindrops, is the resistance coefficient.

7. The method for simulating and predicting the diffusion and deposition of soluble pollutants according to claim 4, characterized in that: The step of accumulating and calculating the amount of pollutants absorbed by raindrops to determine the pollutant deposition information includes: The pollutant deposition mass calculation model is used to calculate the amount of pollutants absorbed by raindrops based on the number of retention time steps to determine the pollutant deposition information; wherein the pollutant deposition mass calculation model is expressed as: in, is the pollutant deposition mass, is the total mass of pollutants absorbed by raindrops at each discrete time.

8. A simulation prediction device for diffusion and deposition of soluble pollutants, characterized in that: The device comprises: A region division module is used to obtain terrain parameter information, meteorological parameter information and pollutant emission parameter information of the area to be measured, and determine the area to be analyzed according to the terrain parameter information, so as to discretize the area to be analyzed and determine a discrete domain, wherein the meteorological parameter information includes: wind condition information and rainfall information; A diffusion analysis module, which determines the boundary conditions of the discrete domain according to the wind condition information, and performs simulation and solution processing on the preset flow equation and the preset diffusion equation according to the pollutant emission parameter information and the boundary conditions to determine the diffusion distribution information of the pollutants, wherein the boundary conditions include: the velocity magnitude, direction, height distribution and surface roughness of the discrete domain boundary, and the preset flow equation includes: the continuity equation and the Navier-Stokes equation; The deposition prediction module calculates and processes the pollutant dissolution information and the pollutant deposition information in the discrete domain based on the rainfall information and the diffusion distribution information of the pollutants, and uses the pollutant dissolution information to update the preset diffusion equation, so as to perform iterative processing within the next discrete time step, and when iterating to the preset maximum simulation time, determine the target pollutant concentration distribution information corresponding to each time step and the target pollutant surface deposition information at each position.

9. A server, characterized in that: The method comprises a processor and a memory, wherein the memory stores computer executable instructions that can be executed by the processor, and the processor executes the computer executable instructions to implement the method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions, and when the computer-executable instructions are called and executed by a processor, the computer-executable instructions prompt the processor to implement the method according to any one of claims 1 to 7.

Images