Meteorology-driven pollutant multi-medium migration simulation method and system

Abstract

The application provides a meteorology-driven pollutant multi-medium migration simulation method and system, and relates to the field of pollutant migration simulation. The method comprises the following steps: determining a plurality of environmental media of a simulation area; determining a plurality of key input parameters of the simulation area based on the plurality of environmental media of the simulation area; obtaining meteorological data of the simulation area; correcting the plurality of key input parameters of the simulation area based on the meteorological data of the simulation area; establishing a multi-medium migration simulation model of the simulation area based on the corrected plurality of key input parameters; and solving the multi-medium migration simulation model of the simulation area to determine the pollutant simulation concentration distribution of each environmental medium, which has the advantage of realizing dynamic migration simulation of pollutants in various environmental media.

Description

Technical Field

[0001] This invention relates to the field of pollutant migration simulation, and in particular to a meteorological-driven multi-media migration and migration simulation method and system for pollutants. Background Technology

[0002] With the continuous advancement of industrialization and urbanization, various pollutants accumulate and interact in environmental media such as the atmosphere, water bodies, and soil, seriously threatening ecosystem security and human health. To study the migration and fate of pollutants in multi-media environments, fugacity models, as an important tool for environmental pre-assessment and risk analysis, have been widely used in simulation studies of pollutant distribution and their transfer pathways. Existing multi-media fugacity models mainly include Level III steady-state models and Level IV dynamic models. Level III models assume the system is in a steady state and are suitable for analyzing the final distribution and fate of pollutants in various media; Level IV models consider time variables and can simulate the change of pollutant concentration over time, making them more suitable for dynamic prediction under complex environmental scenarios.

[0003] Existing multi-media pollutant simulation technologies suffer from insufficient consideration of meteorological factors. In real-world environments, the migration and transformation of pollutants are highly dependent on the dynamic changes in meteorological conditions such as temperature, humidity, wind speed, and precipitation. However, traditional fugacity models mostly use fixed parameters or simplified treatments, failing to accurately reflect the dynamic behavior of pollutants under variable meteorological conditions. These issues limit the effectiveness of existing models in regional environmental management and pollution control strategy development.

[0004] Therefore, there is a need to provide a meteorological-driven multi-media migration and homing simulation method and system for pollutants to achieve dynamic homing simulation of pollutants in various environmental media. Summary of the Invention

[0005] This invention provides a meteorological-driven multi-media migration simulation method for pollutants, comprising: identifying multiple environmental media in a simulation area; determining multiple key input parameters of the simulation area based on the multiple environmental media; acquiring meteorological data of the simulation area; correcting the multiple key input parameters of the simulation area based on the meteorological data; establishing a multi-media migration simulation model of the simulation area based on the corrected multiple key input parameters; and solving the multi-media migration simulation model of the simulation area to determine the simulated concentration distribution of pollutants in each environmental media.

[0006] Furthermore, the simulated area includes at least air, water, soil, and sediments.

[0007] Furthermore, the key input parameters of the simulation region include at least the basic physical parameters, environmental migration parameters, fugacity capacity, and interphase migration parameters of each environmental medium.

[0008] Furthermore, the meteorological data of the simulated area includes at least temperature, humidity, and solar radiation intensity; based on the meteorological data of the simulated area, several key input parameters of the simulated area are corrected, including: based on temperature, the pollutant diffusion coefficient is corrected; based on temperature, the pollutant chemical reaction half-life is corrected; based on temperature and solar radiation intensity, the saturated vapor pressure is corrected; based on temperature, the solubility equilibrium constant and Henry's constant are corrected.

[0009] Furthermore, a multi-media migration simulation model for the simulated area is established, including: establishing a mass transport model characterizing the pollutant migration process; constructing migration rate and flow path models for multiple environmental media; simulating the diffusion process of pollutants within the media based on the adjusted pollutant diffusion coefficient and the diffusion path length of pollutants in the media, determining the pollutant removal mechanism, and establishing a multi-media migration simulation model for the simulated area.

[0010] Furthermore, the multi-media convergence migration simulation model is as follows:

[0011]

[0012] in, For the volume of the atmosphere, The fugacity capacity of the atmosphere. This represents the atmospheric fugacity value. This represents the migration coefficient of water bodies to the atmosphere. The fugacity value of the water body. The soil phase degradation reaction rate coefficient, The migration coefficient from the atmosphere to water bodies. The migration coefficient from the atmosphere to the soil. The advection mobility coefficient in the atmospheric phase. This represents the degradation rate coefficient in the atmospheric phase. The volume of the water body Fugacity capacity of water bodies The fugacity value of the water body. The rate at which pollutants are discharged into water bodies. The pollutant flux is a horizontal flow input. for, The migration coefficient of sediments into water bodies. The migration coefficient from water to sediment. The advection migration coefficient in the water phase. For the volume of soil, For soil fugacity capacity, This represents the fugacity value of the soil. The migration coefficient from the atmosphere to the soil. This represents the soil migration coefficient to the atmosphere. The migration coefficient from soil to water. The volume of the sediment. Fugacity capacity of sediments, The fugacity value of the sediment. The migration coefficient from water to sediment. This is the migration coefficient of sediments into the water body.

[0013] Furthermore, the multi-media migration simulation model of the simulated area is solved to determine the simulated concentration distribution of pollutants in each environmental medium, including: solving the fugacity of each environmental medium; and calculating the simulated concentration distribution of pollutants in each environmental medium through the concentration-fugacity relationship.

[0014] Furthermore, the concentration-fugacity relationship is as follows:

[0015]

[0016] in, This refers to the concentration of pollutants in the environmental medium. The fugacity capacity of the environmental medium. It is called "future degree".

[0017] Furthermore, the method also includes visualizing the simulated concentration distribution of pollutants for each environmental medium.

[0018] This invention provides a meteorological-driven multi-media pollutant migration and facilitation simulation system for use in the aforementioned meteorological-driven multi-media pollutant migration and facilitation simulation method. The system includes: a multi-media environmental modeling module for determining multiple environmental media in the simulation area; the multi-media environmental modeling module is also used to determine multiple key input parameters of the simulation area based on the multiple environmental media; a meteorological data input module for acquiring meteorological data of the simulation area; the multi-media environmental modeling module is also used to correct the multiple key input parameters of the simulation area based on the meteorological data; the multi-media environmental modeling module is also used to establish a multi-media migration and facilitation simulation model of the simulation area based on the corrected multiple key input parameters; and a pollutant migration and transformation simulation module for solving the multi-media migration and facilitation simulation model of the simulation area to determine the simulated pollutant concentration distribution of each environmental media.

[0019] Compared with existing technologies, the meteorological-driven multi-media migration simulation method and system for pollutants provided by this invention have at least the following beneficial effects:

[0020] By introducing high spatiotemporal resolution meteorological data as a driving factor, the migration and transformation processes of pollutants in various environmental media such as the atmosphere, water bodies, and soil are coupled and simulated, effectively improving simulation accuracy and dynamic response capabilities. Compared with existing methods that only model based on static parameters or single media, this method and system can comprehensively reflect the time-varying evolution process of pollutants in complex environmental systems, maintaining high simulation stability and accuracy even under extreme meteorological conditions. Furthermore, the modular design of this method and system facilitates integration with different types of pollution source inventories and geographic information systems, significantly reducing model deployment and operation costs and enhancing the model's versatility and practicality. The implementation of the above technical solutions not only improves the scientific rigor and predictive ability of pollutant fate simulation but also provides reliable technical support for environmental risk assessment and pollution prevention and control. Attached Figure Description

[0021] This specification will be further described by way of exemplary embodiments, which will be described in detail with reference to the accompanying drawings. These embodiments are not limiting; in these embodiments, the same reference numerals denote the same structures, wherein:

[0022] Figure 1 This is a flowchart illustrating a meteorological-driven multi-media migration simulation method for pollutants, as shown in some embodiments of this specification.

[0023] Figure 2 This is a schematic flowchart illustrating the process of establishing a multi-media convergence migration simulation model of a simulation region according to some embodiments of this specification;

[0024] Figure 3 This is a schematic diagram of a meteorology-driven multi-media migration and homing simulation system for pollutants, as shown in some embodiments of this specification.

[0025] Figure 4 This is a schematic diagram illustrating the transport of contaminants at the air / water interface according to some embodiments of this specification;

[0026] Figure 5 This is a schematic diagram illustrating the transport of contaminants at the water / sediment interface according to some embodiments of this specification;

[0027] Figure 6 This is a schematic diagram illustrating the transport of pollutants at the air / soil interface according to some embodiments of this specification;

[0028] Figure 7 This is a flowchart illustrating a hierarchical path filtering strategy according to some embodiments of this specification;

[0029] Figure 8This is a schematic diagram illustrating the migration process and flux of nitrate contaminants according to some embodiments of this specification. Detailed Implementation

[0030] To more clearly illustrate the technical solutions of the embodiments in this specification, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are merely some examples or embodiments of this specification. For those skilled in the art, these drawings can be applied to other similar scenarios without creative effort. Unless obvious from the context or otherwise specified, the same reference numerals in the drawings represent the same structures or operations.

[0031] Figure 1 This is a flowchart illustrating a meteorology-driven multi-media migration simulation method for pollutants, as shown in some embodiments of this specification. Figure 1 As shown, the meteorological-driven multi-media migration simulation method for pollutants may include the following steps.

[0032] Step 110: Determine multiple environmental media in the simulation area.

[0033] The simulated area includes at least several environmental media, including air, water, soil, and sediment. Air refers to the atmosphere of the simulated area, including the troposphere and other near-surface atmospheres, which is the primary site for pollutant volatilization, diffusion, and deposition. Water includes surface water bodies such as rivers, lakes, and oceans, and is an important medium for pollutant dissolution, migration, and transformation. Soil is a loose layer of material covering the earth's surface, including minerals, organic matter, air, and moisture, and is an important medium for pollutant adsorption, degradation, and migration. Sediment refers to solid particulate matter accumulated at the bottom of water bodies (such as rivers, lakes, and oceans), mainly composed of inorganic matter (such as silt and clay) and organic matter (such as animal and plant remains). Through physical, chemical, and biological processes, suspended particulate matter in water gradually settles and accumulates at the bottom, forming a sediment layer.

[0034] Multiple environmental media in the simulated area can be identified based on the migration and transformation processes of pollutants in the environment. For example, pollutants may volatilize from the soil into the atmosphere, or seep into the soil and sediments from water bodies, and then accumulate in organisms through the food chain. Therefore, the multiple environmental media in the simulated area can include not only air, water, soil, and sediments, but also organisms.

[0035] Step 120: Determine several key input parameters of the simulation area based on multiple environmental media within the simulation area.

[0036] The key input parameters for the simulation region include at least the basic physical parameters, environmental migration parameters, fugacity capacity, and interphase migration parameters for each environmental medium.

[0037] Specifically, the basic physical parameters of environmental media can include surface area, depth, volume, density, volume fraction, and organic carbon content, as shown in Table 1. Spatial analysis and processing are performed using the ArcGIS platform to construct a multi-media environmental system encompassing the atmosphere, water, and soil within the region. Geographic information data such as topography and land use are incorporated to achieve three-dimensional spatial modeling. Finally, regional area parameters are obtained through landscape pattern statistical calculations using Fragstats software.

[0038] Table 1

[0039] parameter unit Atmospheric area <![CDATA[m 2 ]]> water area <![CDATA[m 2 ]]> Soil area <![CDATA[m 2 ]]> Sediment area <![CDATA[m 2 ]]> Mean Atmospheric Height m Average water depth m Average soil depth m Mean sediment depth m Gas volume fraction in water Gas volume fraction in soil Gas volume fraction in sediments Volume fraction of particulate matter in the atmosphere Volume fraction of particulate matter in urban atmosphere Volume fraction of particulate matter in water Solid volume fraction in soil Solid volume fraction in sediments Water phase volume fraction in the atmosphere Water phase volume fraction in soil Water phase volume fraction in sediments Atmospheric density <![CDATA[kg / m 3 ]]> Water density <![CDATA[kg / m 3 ]]> Density of particulate matter in the atmosphere <![CDATA[kg / m 3 ]]> Density of particulate matter in water <![CDATA[kg / m 3 ]]> Density of solids in soil <![CDATA[kg / m 3 ]]> Density of solids in sediments <![CDATA[kg / m 3 <!-- 4 -->]]> Atmospheric organic carbon content Organic carbon content in water Soil organic carbon content Organic carbon content in sediments

[0040] The fugacity capacity can be calculated as shown in Table 2.

[0041] Table 2

[0042]

[0043] The Z-value represents the proportional relationship between the concentration (mol / m³) of a pollutant in a specific phase and its partial pressure (Pa) in the atmosphere under equilibrium conditions. A Z represents the partition coefficient of pollutants in the gas phase. PA Z is the distribution coefficient of pollutants in the atmospheric particulate phase. BA v is the distribution coefficient of pollutants in the entire atmospheric phase. PA Z represents the volume fraction of particulate matter in the atmosphere (m³ particles / m³ air), a dimensionless quantity. It indicates the proportion of particulate matter volume in a unit volume of air. W Z is the partition coefficient of pollutants in the aqueous phase, H is the Henry's constant (Pa·m³ / mol), and Z is the distribution coefficient of pollutants in the aqueous phase. PW ρ is the distribution coefficient of pollutants in the suspended particulate phase of water. PW K represents the density of suspended particles (kg / m³). OC For organic carbon-water distribution, f OC . PW Z represents the mass fraction of organic carbon in suspended particulate matter (kg OC / kg particles). BW v is the distribution coefficient of pollutants in the entire water phase. PW Z represents the volume fraction of suspended particles. S ρ is the distribution coefficient of pollutants in the soil phase. S Soil particle density (kg / m³), K OW f is the octanol-water partition coefficient. OC.SZ represents the mass fraction of soil organic carbon, 0.41 is an empirical coefficient representing the ratio of organic carbon adsorption capacity to octanol adsorption capacity. BS v is the distribution coefficient of pollutants in the entire soil phase. A v represents the volume fraction of the gas phase in the soil. W v represents the volume fraction of water in the soil. S Z represents the volume fraction of the solid phase in the soil. Sed ρ is the partition coefficient of pollutants in the sedimentary phase. Sed f is the density of sediment particles. OC.Sed Z represents the mass fraction of organic carbon in sediments. BT v is the distribution coefficient of pollutants in the entire sediment phase. Sed This represents the volume fraction of the solid phase in the sediment.

[0044] Environmental migration parameters can characterize the migration rate of pollutants between any two environmental media (such as air and water, water and sediment, air and soil, etc.), as shown in Table 3.

[0045] Table 3

[0046] parameter symbol unit Gas-side mass transfer coefficient at the gas / water interface <![CDATA[k VA ]]> m / h Water-side mass transfer coefficient at the air / water interface <![CDATA[k VW ]]> m / h Precipitation rate <![CDATA[U R ]]> m / h Clearance Q Dry sedimentation rate <![CDATA[U Q ]]> m / h Gas-side mass transport coefficient at the gas / soil interface <![CDATA[k EA ]]> m / h Molecules have long diffusion paths in soil <![CDATA[Y3]]> m Molecular diffusion coefficient in air <![CDATA[B MA ]]> <![CDATA[m 2 / h]]> Molecular diffusion coefficient in water <![CDATA[B MW ]]> <![CDATA[m 2 / h]]> Effective diffusion coefficient of molecules in sediment pore water <![CDATA[m 2 / h]]> Water runoff rate from soil <![CDATA[U WW ]]> m / h Loss rate of solids from soil <![CDATA[U EW ]]> m / h Water-side mass transport coefficient at the water / sediment interface <![CDATA[k SW ]]> m / h Length of diffusion path of molecules in sediments <![CDATA[Y4]]> m sediment settling rate <![CDATA[U DP ]]> m / h sediment resuspension rate <![CDATA[U RS ]]> m / h Sediment burial rate <![CDATA[U BS ]]> m / h Gas-side mass transfer coefficient at the gas / impermeable layer interface Kaf m / h Falling rate constant of litter Rlf <![CDATA[h -1 ]]> Leaf wax erosion mass transfer coefficient Kwe m / h Impermeable layer erosion rate constant W Impermeable layer - water mass transfer coefficient Krw m / h Air / pine needle interface boundary layer diffusion mass transport coefficient kav m / h Leaf epidermal penetration rate Pcv Leaf cutoff fraction FrUF Particle dry deposition rate to leaf surface Uaf m / h Impermeable layer permeability rate Pcf m / s Migration rate to higher altitudes <![CDATA[U S ]]> m / h Volume fraction of aerosols <![CDATA[Vf a ]]> Filtration rate of water from bottom sediment to groundwater <![CDATA[U L ]]> <![CDATA[m 3 / m 2 ×h]]>

[0047] The calculation methods for environmental migration parameters obtained through literature review, empirical value calibration, or field measurement are shown in Table 4.

[0048] Table 4

[0049]

[0050] Among them, D VW K is the atmospheric-water diffusion and migration parameter. VA The atmospheric gas-phase mass transfer coefficient (m / h) reflects the diffusion efficiency of pollutants from the bulk atmosphere to the interface. W K represents the surface area of ​​the water body. VW D is the water phase mass transfer coefficient (m / h) on the water body side, reflecting the diffusion efficiency of pollutants from the interface to the main body of the water body. RW U represents the atmospheric-water vapor phase wet deposition migration parameter. R For the rainfall rate, D PW Z represents the atmospheric-water particle wet deposition migration parameter. PA Z is the atmospheric particulate phase value, and Q is the scour ratio, an empirical parameter reflecting the efficiency of rainwater in removing particulate matter (usually 1). ), v PA D represents the volume fraction of particulate phase in the atmosphere. DW U represents the atmospheric-water particle dry deposition migration parameter. P The dry settling rate of particulate matter (m / h) is related to particle size, density, and meteorological conditions.A-W D represents the environmental migration parameters for the overall atmosphere-water process. VS k is the atmospheric-soil diffusion and migration parameter. AS The atmospheric gas-phase mass transfer coefficient (m / h) is related to wind speed and surface roughness. A S Y is the soil surface area (m²), and Y3 is the effective thickness of the soil air layer (m), representing the soil pore depth. A3 B W3 To separate the gas / water phase mass transfer coefficients (m / h) in soil pores, reflecting the migration efficiency of pollutants within the soil pores, D RS D represents the atmospheric-soil vapor phase wet deposition migration parameter. PS D represents the atmospheric-soil particle wet deposition and migration parameter. DS D represents the atmospheric-soil particle dry deposition migration parameter. A-S D represents the diffusion and migration parameters of the overall atmosphere-soil process. VS D is an environmental migration parameter between soil and atmosphere. SRW U is a soil loss and migration parameter for soil-water exchange. SW D represents the soil loss rate (m / h). WRW U is a parameter for rainwater runoff migration in soil-water systems. WW D is the runoff rate (m / h). SW D is the environmental migration parameter for the overall soil-water process. Y k is a parameter for sediment-water diffusion and migration. Sed Let A be the mass transfer coefficient of the water phase on the water body side (m / h), describing the resistance of the overlying water boundary layer. Sed A is the surface area of ​​the sediment (m²). S Y is the soil surface area (m²), Y4 is the effective thickness of the sediment layer (m), and B is the depth to which pollutants diffuse in the sediment pores. W4 D is the mass transfer coefficient of water in sediment pores (m / h), reflecting the diffusion efficiency of pore water. RSed Z represents the sediment-water resuspension migration parameter. PSed U represents the Z-value of sediment grain phase (mol / (m³·Pa)). RSed A represents the sediment resuspension rate (m / h), the volume of sediment resuspended per unit area per unit time, and A's resuspension rate (m / h). Sed D is the surface area of ​​the sediment (m²). Sed-W D is an environmental migration parameter for the overall sediment-water process. DSed D is the sedimentation and migration parameter of water and sediment. W-Sed Here, i represents the environmental migration parameter for the overall water-sediment process, and D represents the medium number (e.g., water, atmosphere, soil). A(i)To represent the amount of mass transported in medium i due to advection (unit: mol / s), G (i) Z is the volumetric flow rate (m³ / h) of the phase entering medium i. (i) D is the total phase Z value of medium i (mol / (m³·Pa)). R(i) K represents the degradation rate of a substance in medium i (unit: mol / s). It indicates the amount of substance that degrades in medium i per unit time, reflecting the rate of degradation of the substance in the medium. (i) V is the advection output rate constant of medium i. (i Z is the volume (m³) of medium i. (i) Let Z be the total phase Z value of medium i (mol / (m³·Pa)).

[0051] Step 130: Obtain meteorological data for the simulated area.

[0052] The meteorological data for the simulated area includes at least temperature, humidity, and solar radiation intensity.

[0053] Meteorological data for the simulated area may also include precipitation, wind speed, etc.

[0054] Step 140: Based on the meteorological data of the simulated area, correct several key input parameters of the simulated area.

[0055] Specifically, several key input parameters are highly sensitive to changes in meteorological conditions, and their values ​​may fluctuate significantly with changes in meteorological conditions. Therefore, it is necessary to correct several key input parameters of the simulated area based on meteorological data of the simulated area.

[0056] In some embodiments, step 140 specifically includes:

[0057] The pollutant diffusion coefficient is corrected based on temperature;

[0058] The half-life of pollutant chemical reactions is corrected based on temperature;

[0059] The saturated vapor pressure is corrected based on temperature and solar radiation intensity.

[0060] The solubility equilibrium constant and Henry's constant are corrected based on temperature.

[0061] For example, the pollutant diffusion coefficient can be corrected based on the following formula:

[0062]

[0063] in, The diffusion coefficient is given by the temperature T. Reference temperature The diffusion coefficient, E, is the activation energy or empirical energy term, reflecting the sensitivity of diffusion behavior to temperature. It is a core parameter for quantifying the temperature dependence of organic matter diffusion processes. E is an empirical value; the range of E for air is 5-15, for water it is 15-25, and for soil it is 30-80. R is the ideal gas constant. Let be the absolute temperature, and exp be the natural exponential function.

[0064] This method uses a unified formula to correct the diffusion coefficient. For different environmental media (air, water, soil, sediment), a specific molecular diffusion coefficient representing the characteristics of each medium is set in this formula. This enables targeted modifications to the diffusion behavior of different media.

[0065] The half-life of pollutant chemical reactions is corrected based on the following formula:

[0066]

[0067] in, The half-life at temperature T. Reference temperature The half-life of the next.

[0068] The saturated vapor pressure is corrected based on the following formula:

[0069]

[0070] in, Let T be the saturated vapor pressure at temperature T. Reference temperature The saturated vapor pressure below It is the enthalpy of vaporization (unit: J / mol).

[0071] The solubility equilibrium constant and Henry's constant are corrected based on the following formula:

[0072]

[0073] in, Here is Henry's constant at temperature T. Reference temperature Henry's constant under the following conditions This is the enthalpy of dissolution (unit: J / mol).

[0074] A quantitative response relationship between meteorological factors and migration flux is established by solving multiple equations simultaneously.

[0075] Step 150: Based on the corrected key input parameters, establish a multi-media migration simulation model for the simulation area.

[0076] Figure 2 This is a flowchart illustrating the process of establishing a multi-media convergence migration simulation model of a simulation region according to some embodiments of this specification, such as... Figure 2 As shown, step 150 specifically includes:

[0077] Establish a mass transfer model to characterize the migration process of pollutants. Specifically, define and calculate the mass transfer coefficient for different media interfaces (such as air / water, water / sediment, air / soil, etc.). The mass transfer coefficient is dynamically corrected according to climatic conditions such as temperature, wind speed, and humidity to reflect the changes in the pollutant transfer capacity of the interface.

[0078] Construct migration rate and flow path models for multiple environmental media. Specifically, calculate the migration rate between multiple media, as well as various migration paths such as rainfall washout, evaporation, sedimentation, runoff, resuspension, and burial.

[0079] Based on the adjusted pollutant diffusion coefficient and the diffusion path length of pollutants in the medium, the diffusion process of pollutants inside the medium is simulated, the pollutant removal mechanism is determined, and a multi-media migration simulation model of the simulation area is established. Specifically, the diffusion coefficient is adjusted according to temperature and humidity to simulate the diffusion process of pollutants inside the medium, and processes such as removal rate and leaf surface deposition are set to reflect the effect of pollutant adsorption, degradation and other mechanisms.

[0080] The mass transport coefficient represents the flux rate (mol / h) of a contaminant from medium A across the interface to medium B.

[0081] Figure 4 This is a schematic diagram illustrating the transport of pollutants at the air / water interface according to some embodiments of this specification, such as... Figure 4 As shown, the total mass transfer coefficient at the air / water interface (Air-Water) It can be calculated based on the following formula:

[0082]

[0083] in, is the gas-phase mass transfer coefficient (m / h). The gas constant is Temperature (K) is the mass transfer coefficient on the aqueous side (m / h), and H is the Henry's constant (Pa·m³ / mol).

[0084] Figure 5 This is a schematic diagram illustrating the transport of contaminants at the water / sediment interface according to some embodiments of this specification, such as... Figure 5 As shown, the total mass transport coefficient at the water-sediment interface. It can be calculated based on the following formula:

[0085]

[0086] in, The effective diffusion coefficient of sediment (m² / h) is given. It is the thickness of the bottom water layer (m, usually 0.1~1 mm). The diffusion coefficient of molecules in water (m² / h) is given. The porosity of the sediment is 0.6~0.9. The tortuosity factor is (1.5~3.0).

[0087] Figure 6 This is a schematic diagram illustrating the transport of pollutants at the air / soil interface according to some embodiments shown in this specification, such as... Figure 6 As shown, the total mass transport coefficient at the air-soil interface (Air-Soil) It can be calculated based on the following formula:

[0088]

[0089] in, The effective diffusion coefficient of air in the soil (m² / h) is given. The thickness of the static air layer on the soil surface (m, approximately 0.5~2 cm). It is the molecular diffusion coefficient in the atmosphere (m² / h). This refers to the air porosity of the soil.

[0090] A four-level media system (air-water-soil-sediment) is constructed, with migration paths based on physical connectivity and pollutant characteristics. The air / water interface achieves bidirectional mass exchange through molecular diffusion and turbulent mixing, while the water / sediment interface achieves bidirectional exchange through pore diffusion and resuspension. From air to soil, dry and wet deposition dominate. From soil to water, rainfall-driven runoff occurs; from water to sediment, gravity deposition occurs. The path network satisfies the law of conservation of mass.

[0091] The transport flux of each migration path is characterized by the mass transport coefficient (D value), the calculation of which integrates physical mechanisms and chemical partitioning principles:

[0092] 1. Dominant diffusion pathways (e.g., gas / water exchange):

[0093]

[0094] in, The mass transfer coefficient of the gas / water exchange migration path. For air fugacity capacity, Let A be the fugacity capacity of water, and let A be the area of ​​the air-water interface (in m²), which is the common contact interface between air and water. For gas mass transfer coefficient, Let H be the mass transfer coefficient on the water side, H be Henry's constant, R be the gas constant, and T be the absolute temperature.

[0095] 2. Dominant advection pathways (e.g., soil runoff):

[0096]

[0097] in, The mass transport coefficient of soil runoff. For water phase fugacity capacity, The soil-water contact area. The runoff rate (m / h) is dynamically adjusted by rainfall intensity.

[0098] 3. Particle settling path:

[0099]

[0100] in, The mass transfer coefficient is the particle settling path. For sediment fugacity capacity, The area of ​​the water-sediment interface. The value represents the settling rate.

[0101] The migration path intensity responds in real time to changes in environmental conditions, with key regulating factors including temperature, wind speed, precipitation, humidity, and sunshine.

[0102] The effectiveness of migration pathways is constrained by the physicochemical properties of pollutants. For example, a Henry's constant > 10-3 atm·m³ / mol activates significant gas / water exchange pathways. Another example is log KOW (i.e., n-octanol-water partition coefficient) > 4, which enhances sediment adsorption and bioaccumulation pathways. Yet another example is that substances with long half-lives activate long-distance migration pathways.

[0103] Figure 7 This is a flowchart illustrating a hierarchical path filtering strategy according to some embodiments of this specification, such as... Figure 7 As shown, a hierarchical path filtering strategy is used to filter paths. As an example only, the migration process and flux of nitrate pollutants are shown... Figure 8 As shown.

[0104] In some embodiments, the multi-media convergence migration simulation model is as follows:

[0105]

[0106] in, For the volume of the atmosphere, The fugacity capacity of the atmosphere. This represents the atmospheric fugacity value. This represents the migration coefficient of water bodies to the atmosphere. The fugacity value of the water body. The soil phase degradation reaction rate coefficient, The migration coefficient from the atmosphere to water bodies. The migration coefficient from the atmosphere to the soil. The advection mobility coefficient in the atmospheric phase. This represents the degradation rate coefficient in the atmospheric phase. The volume of the water body Fugacity capacity of water bodies The fugacity value of the water body. The rate at which pollutants are discharged into water bodies. The pollutant flux is a horizontal flow input. for, The migration coefficient of sediments into water bodies. The migration coefficient from water to sediment. The advection migration coefficient in the water phase. For the volume of soil, For soil fugacity capacity, This represents the fugacity value of the soil. The migration coefficient from the atmosphere to the soil. This represents the soil migration coefficient to the atmosphere. The migration coefficient from soil to water. The volume of the sediment. Fugacity capacity of sediments, The fugacity value of the sediment. The migration coefficient from water to sediment. This is the migration coefficient of sediments into the water body.

[0107] The initial condition is set to t=0. , , and All are f0, representing the initial fugacity value of pollutants in each environmental medium, i.e., the initial background value of the multi-media migration simulation model.

[0108] The mass balance equations describing the equilibrium state of pollutants in a multi-media environment are solved using the matrix operation capabilities of Matlab software. By transforming the equilibrium equations into matrix form, and leveraging Matlab's efficient matrix computation capabilities, the solutions to the equation system can be obtained quickly, thereby determining the equilibrium concentration distribution of pollutants in each medium.

[0109] For the system of differential equations, such as the dynamic mass balance equation, which describes the time-varying changes in the migration and transformation of pollutants, the Runge-Kutta method in Matlab is used for solution. The Runge-Kutta method is a commonly used numerical solution method for differential equations. It can calculate the numerical solutions of the system of differential equations at different time points with relatively high accuracy, thereby simulating the dynamic migration and transformation process of pollutants in a multi-media environment.

[0110] Step 160: Solve the multi-media migration simulation model of the simulated area to determine the simulated concentration distribution of pollutants in each environmental medium.

[0111] Specifically, it includes:

[0112] Solve for the fugacity of each environmental medium;

[0113] The simulated concentration distribution of pollutants in each environmental medium is calculated using the concentration-fugacity relationship.

[0114] In some embodiments, the concentration-fugacity relationship is as follows:

[0115]

[0116] in, This refers to the concentration of pollutants in the environmental medium. Fugacity capacity refers to the fugacity capacity of an environmental medium, reflecting the medium's "storage capacity" for pollutants. It is called "future degree".

[0117] In some embodiments, the method further includes:

[0118] The simulated concentration distribution of pollutants in each environmental medium is visualized, specifically presented in the form of pollutant concentration distribution maps, providing an intuitive basis for environmental management. For example, the concentration distribution of pollutants in environmental media can be displayed using two-dimensional or three-dimensional graphics. Two-dimensional maps use contour maps, heat maps, or filled maps, while three-dimensional maps use surface maps or volume rendering maps. The maps should clearly label coordinate axes (such as spatial location, depth, etc.), concentration scales (such as color or height to indicate concentration magnitude), and the type of environmental medium (such as air, water, soil, etc.). Example: On a two-dimensional map, different colors are used to represent the concentration distribution of pollutants in soil, with red areas representing high concentrations and blue areas representing low concentrations.

[0119] For time-series data, dynamic videos or animations can be created to show how pollutant concentrations change over time. Dynamic visualization helps users understand the migration trends and fate of pollutants.

[0120] Visualized results provide environmental management departments with intuitive decision-making support. For example, identifying high-concentration areas allows for the determination of priority remediation zones. Observing the migration paths of pollutants enables the prediction of future distribution trends, allowing for proactive prevention and control measures.

[0121] As an example, suppose a chemical spill occurs in a certain area, requiring migration simulation and risk assessment of pollutants in the soil. This method can simulate the concentration distribution of pollutants in the soil and present it as a concentration distribution map. Environmental management departments can use the visualization results to identify high-concentration areas and take corresponding remediation measures. Furthermore, dynamic visualization can show the migration trends of pollutants, providing a basis for long-term environmental monitoring and control.

[0122] Understandably, in practical applications, such as simulating the multi-media fate processes of typical industrial pollutants (e.g., polycyclic aromatic hydrocarbons) in a region, by integrating real-time meteorological driving data, the migration trends and concentration change trajectories of pollutants in the atmosphere-water-soil-sediment system can be reflected in real time. Compared with traditional static or simplified simulation methods, this method can more accurately capture the pollution response behavior of short-term abrupt events (e.g., runoff scouring after heavy rainfall), exhibiting higher simulation accuracy and response capabilities. Furthermore, users can flexibly change or add pollutant types and meteorological sources, enabling the reuse of a single model across multiple scenarios, significantly reducing development and deployment costs.

[0123] This method not only theoretically achieves refined and dynamic simulation of pollutant multi-media migration and homing, but also has good scalability and adaptability in practical operation, and can be widely applied to multiple environmental science and engineering management fields such as pollution risk assessment, pollution source tracing, and environmental emergency response.

[0124] This method innovatively incorporates real-time or historical high spatiotemporal resolution meteorological data (including precipitation, temperature, wind speed, wind direction, humidity, etc.) into the pollutant migration simulation process as the main external driving factor. This enables dynamic driving of pollutant migration and transformation processes among media such as the atmosphere, water, and soil. By integrating multi-source data including regional topography, land use, and hydrology, a multi-media environmental model encompassing the atmosphere, water, soil, and sediments is constructed, providing high-precision spatial support for pollutant migration and transformation across different media. It possesses powerful data integration and boundary condition configuration capabilities, making it suitable for multi-field coupled simulation and dynamic change analysis of complex environmental systems. The multi-media fate migration simulation model integrates the migration paths and process mechanisms of pollutants among multiple media, such as diffusion, deposition, volatilization, adsorption, and biodegradation. Quantitative calculations are performed based on the coupling of pollutant physicochemical properties, environmental media characteristics, and meteorological conditions, enabling dynamic simulation of fate under complex conditions, particularly suitable for predicting the fate paths of organic pollutants and heavy metals.

[0125] Figure 3 This is a schematic diagram of a module of a meteorology-driven multi-media pollutant fate and migration simulation system, as shown in some embodiments of this specification. Figure 3 As shown, the meteorological-driven multi-media migration and migration simulation system for pollutants can include a multi-media environmental modeling module, a meteorological data input module, a pollutant migration and transformation simulation module, and a visualization module.

[0126] The multi-media environment modeling module is used to determine multiple environmental media in the simulation area;

[0127] The multi-media environment modeling module is also used to determine multiple key input parameters of the simulation area based on multiple environmental media within the simulation area;

[0128] The meteorological data input module is used to acquire meteorological data for the simulated area;

[0129] The multi-media environment modeling module is also used to correct multiple key input parameters of the simulated area based on meteorological data of the simulated area;

[0130] The multi-media environment modeling module is also used to establish a multi-media migration simulation model of the simulation area based on several corrected key input parameters;

[0131] The pollutant migration and transformation simulation module is used to solve the multi-media migration simulation model of the simulation area and determine the simulated concentration distribution of pollutants in each environmental medium.

[0132] The visualization module can be used to visualize the simulated concentration distribution of pollutants in each environmental medium.

[0133] The meteorological-driven multi-media migration simulation system for pollutants can be used to perform meteorological-driven multi-media migration simulation methods for pollutants, which will not be elaborated here.

[0134] Finally, it should be understood that the embodiments described in this specification are merely illustrative of the principles of the embodiments described herein. Other variations may also fall within the scope of this specification. Therefore, alternative configurations of the embodiments described herein are intended to be illustrative rather than limiting, and should be considered consistent with the teachings of this specification. Accordingly, the embodiments described herein are not limited to those explicitly introduced and described herein.

Claims

1. A meteorological-driven multi-media migration simulation method for pollutants, characterized in that, include: Determine multiple environmental media within the simulation area; Based on multiple environmental media in the simulation area, several key input parameters of the simulation area are determined. Obtain meteorological data for the simulated area; Based on meteorological data from the simulated region, several key input parameters of the simulated region were corrected. Based on the corrected key input parameters, a multi-media migration simulation model for the simulation area is established. The multi-media migration simulation model of the simulated area was solved to determine the simulated concentration distribution of pollutants in each environmental medium; Based on several corrected key input parameters, a multi-media migration simulation model for the simulation area is established, including: Establish a mass transfer model to characterize the pollutant migration process. Specifically, define and calculate the mass transfer coefficient for different media interfaces. The mass transfer coefficient is dynamically corrected according to temperature, wind speed, and humidity. The media interfaces include air / water interface, water / sediment interface, and air / soil interface. Construct migration rate and flow path models for multiple environmental media. Specifically, calculate the migration rate between multiple media, as well as the migration paths of rainfall washout, evaporation, sedimentation, runoff, resuspension, and burial. Based on the adjusted pollutant diffusion coefficient and the diffusion path length of pollutants in the medium, the diffusion process of pollutants inside the medium is simulated, the pollutant removal mechanism is determined, and a multi-media migration simulation model of the simulation area is established. Specifically, the diffusion coefficient is adjusted according to temperature and humidity to simulate the diffusion process of pollutants inside the medium, and the removal rate and leaf surface deposition process are set. The mass transport coefficient represents the flux rate at which a contaminant travels from one medium across an interface to another. Total mass transfer coefficient at the air / water interface Calculated based on the following formula: in, The mass transfer coefficient on the gas phase side is... The gas constant is... For temperature, Here, H is the mass transfer coefficient on the aqueous side, and H is Henry's constant. Total mass transport coefficient at the water / sediment interface Calculated based on the following formula: in, The effective diffusion coefficient of sediments, The thickness of the bottom water layer, The diffusion coefficient of molecules in water, The value represents the sediment porosity, ranging from 0.6 to 0.

9. This is the tortuosity factor, with a value range of 1.5 to 3.0; Total mass transport coefficient at the air / soil interface Calculated based on the following formula: in, The effective diffusion coefficient of air in soil. The thickness of the static air layer on the soil surface. The diffusion coefficient of molecules in the atmosphere. The air porosity of the soil; A four-level media system corresponding to air-water-soil-sediment is constructed. The migration path is based on physical connectivity and pollutant characteristics. The air / water interface realizes bidirectional mass exchange through molecular diffusion and turbulent mixing. The water / sediment interface realizes bidirectional exchange through pore diffusion and resuspension. From air to soil, dry and wet deposition dominate. From soil to water, the runoff process is driven by rainfall. From water to sediment, the gravity deposition process is used. The path network satisfies the law of conservation of mass. The transport flux of each migration path is characterized by the mass transport coefficient, the calculation of which integrates physical mechanisms and chemical partitioning principles: (1) Dominant diffusion pathway: in, The mass transfer coefficient of the gas / water exchange migration path. For air fugacity capacity, Let A be the fugacity capacity of water, and A be the area of ​​the air-water interface, i.e., the common contact interface between air and water. For gas mass transfer coefficient, Here, H is the water-side mass transfer coefficient, H is Henry's constant, R is the gas constant, and T is the absolute temperature. (2) Dominant advection path: in, The mass transport coefficient of soil runoff. For water phase fugacity capacity, The soil-water contact area. The runoff rate is dynamically adjusted by rainfall intensity. (3) Particle settling path: in, The mass transfer coefficient is the particle settling path. For sediment fugacity capacity, The area of ​​the water-sediment interface. The value represents the settling rate.

2. The meteorological-driven multi-media migration simulation method for pollutants according to claim 1, characterized in that, The simulated area includes at least air, water, soil, and sediments.

3. The meteorological-driven multi-media migration simulation method for pollutants according to claim 2, characterized in that, The key input parameters of the simulation region include at least the basic physical parameters, environmental migration parameters, fugacity capacity, and interphase migration parameters of each environmental medium.

4. The meteorological-driven multi-media migration simulation method for pollutants according to any one of claims 1-3, characterized in that, The multi-media migration simulation model of the simulated region is solved to determine the simulated pollutant concentration distribution of each environmental medium, including: Solve for the fugacity of each environmental medium; The simulated concentration distribution of pollutants in each environmental medium is calculated using the concentration-fugacity relationship.

5. The meteorological-driven multi-media migration simulation method for pollutants according to any one of claims 1-3, characterized in that, Also includes: The simulated concentration distribution of pollutants for each environmental medium is visualized.

6. A meteorological-driven multi-media migration and homing simulation system for pollutants, characterized in that, The method for performing the meteorology-driven multi-media migration and homing simulation of pollutants according to any one of claims 1-5 includes: The multi-media environment modeling module is used to determine multiple environmental media in the simulation area; The multi-media environment modeling module is also used to determine multiple key input parameters of the simulation area based on multiple environmental media within the simulation area; The meteorological data input module is used to acquire meteorological data for the simulated area; The multi-media environment modeling module is also used to correct multiple key input parameters of the simulated area based on meteorological data of the simulated area; The multi-media environment modeling module is also used to establish a multi-media trend migration simulation model of the simulation area based on several corrected key input parameters; The pollutant migration and transformation simulation module is used to solve the multi-media migration simulation model of the simulation area and determine the simulated concentration distribution of pollutants in each environmental medium.

Images