KR20140125056A

KR20140125056A - Method for contaminant motion standard construct

Info

Publication number: KR20140125056A
Application number: KR20130042680A
Authority: KR
Inventors: 이용주; 이성규; 정자은; 김윤관
Original assignee: 그린에코스 주식회사
Priority date: 2013-04-18
Filing date: 2013-04-18
Publication date: 2014-10-28

Abstract

The present invention discloses a method for constructing a standard model for the operation of a pollutant. The present invention establishes a standard modeling model of pollutant behavior that estimates the behavior of pollutants in the environment, The pollutant concentration in the environment composed of the atmosphere, soil, vegetation and water by the pollutant diffusion can be estimated and the response according to the pollutant diffusion can be made quickly.

Description

A method for contaminant motion standard construction

The present invention relates to a technique for constructing a standard model of pollutant behavior, and more particularly, to construct a model of a pollutant behavior standard model for estimating the behavior of pollutants in the environment, , Vegetation, and water to estimate the pollutant concentration in the environment.

Generally, an insecticide is composed of a composition containing one insecticidal active ingredient and a composition containing two or more insecticidal active ingredients. In case of containing one insecticidal component, it exhibits a single efficacy, and in case of insecticide prepared by using two or more components, it has been widely used because of its superior efficacy over a single preparation.

However, such a mixed preparation generally causes degradation of the active ingredient when a long period of time elapses after the preparation, and thus there is a problem in its stability. In order to solve this problem, JP-A-10-2003-37989 ) Has developed a pesticide composition capable of rescuing larvae as well as larvae of insect pests by mixing seaweed margins and deltamethine as effective ingredients and mixing them with suitable excipients.

However, there has not been proposed a modeling technique for estimating the behavior of insecticides applied in the past when the insecticide is applied to the insecticide as described above. If the concentration of the insecticide can not be estimated based on the behavior of the insecticide, There were many.

SUMMARY OF THE INVENTION Accordingly, the present invention has been made to solve the above-mentioned conventional problems, and it is an object of the present invention to provide a method for predicting the behavior of pollutants (for example, insecticides) (Pollution), which enables to precisely estimate the pollution concentration in the environment composed of air, soil, vegetation and water caused by pollutants, And to provide a method for constructing a standard model of material behavior.

In order to achieve the above object, the method for constructing a standard model for the behavior of pollutants according to the present invention can be classified into intra-media processes and inter-media processes considering air, soil, vegetation, A first step of deriving an input parameter of an environmental behavior mechanism classified into a first step; A second step in which the environmental behavior mechanism derived from the first step measures the amount of change in the behavior of contaminants with time from input variables; A third step of constructing a mass balance equation for the amount of change in behavior measured from the second step; And a fourth step of constructing modeling of a pollutant behavior standard through a mass balance equation composed of the third step; .

In addition, the pollutant is one of an insecticide, a gas, and a chemical substance that is diffused into the atmosphere and pollutes the atmospheric, soil, vegetation, and water of the target medium.

In addition, the input parameters of the media mechanism include advection and wind drift of the applied pollutants, degradation of the pollutants, leaching of the pollutants.

Further, the input parameters of the inter-media mechanism may include volatilization of the contaminants, dry deposition of the contaminants, wet deposition of the contaminants, run-off of the contaminants, Uptake of the substance to the vegetation roots, wash-off by rainfall at the surface of the vegetation for contaminants, and diffusion of the gaseous matter from the atmosphere into the soil and water.

In the third step, the mass balance equation for the media kinetics model is constructed so that the amount of contaminant change in the entire system is represented by the sum of the amount of material change over time in each medium (i-the compartment) as shown in the following equation .

Also, the change amount of the pollutant in the medium is configured to be expressed as a value obtained by subtracting the output rate from the input rate of the pollutant as shown in the following equation.

Further, the mass balance equation for the air is constructed as shown in the following equation.

In addition, the mass balance equation of the soil (Soil) to the bare soil is expressed by the following equation.

In addition, the mass balance equation for the forest soil in the soil (Soil) is as shown in the following equation.

In addition, the mass balance equation for the vegetation cuts in the above vegetation is expressed by the following equation.

In addition, the mass balance equation for the vegetation surfaces in the above vegetation is composed as shown in the following equation.

In addition, the mass balance equation for the water is constructed as shown in the following equation.

In addition, the model constructed by the mass balance equation in the fourth step is a model program that solves the Euler method.

In addition, the drift (%) of the airborne pollutant flow is calculated by a proportional approach based on the empirical formula or a method of solving the amount of movement of the windward wind by applying the meteorological theory. .

In addition, the advection of the pollutants is to check the inflow and outflow of pollutants diffused into the atmosphere to apply the advection and diffusion by the wind at predetermined time intervals.

In the above model, the distance between the position of the soil where the rainfall runoff occurs and the exit point of the watershed is calculated to calculate the reach distance, and the rate of reaching the runoff rate according to the reach distance is calculated according to the simulation time of the model .

In addition, the environmental behavior mechanism, which is the basic data of the pollutant behavior, is configured to utilize the geographic information system (GIS).

As described above, the present invention establishes a standard modeling model of the pollutant behavior that estimates the behavior of pollutants in the environment, and thus, it is possible to estimate the environment of the air, soil, vegetation and water caused by pollutants diffused into the air It is possible to estimate the concentration of the pollutants accurately so that the response to the pollutant diffusion can be promptly performed.

1 is a flow chart illustrating a method for establishing a standard model for contaminant behavior according to an embodiment of the present invention.
2 is a conceptual diagram of a method for constructing a standard model of pollutant behavior according to an embodiment of the present invention.
FIG. 3 is a diagram illustrating a model of a pollutant behavior standard model using wind-induced lean ratios according to an embodiment of the present invention.
FIG. 4 is a diagram illustrating a model of a pollutant behavior standard model applying wind-induced weathering characteristics as an embodiment of the present invention.
FIG. 5 is a diagram of a standard model for pollutant behavior considering atmospheric diffusion, according to an embodiment of the present invention. FIG.
FIG. 6 is an experimental chart showing the spatial composition of a model construction in a river according to an embodiment of the present invention. FIG.
FIG. 7 is an experimental chart showing the spatial structure of the model building in the trench as an embodiment of the present invention. FIG.
8 is an experimental chart showing the spatial composition of a model construction in a pond according to an embodiment of the present invention.

Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings.

2 is a conceptual diagram of a method for constructing a standard model of pollutant behavior according to an embodiment of the present invention, and FIG. 3 is a schematic view of an embodiment of the present invention FIG. 4 is a diagram showing a construction of a standard model for contaminant behavior using the wind-induced microfossilism according to an embodiment of the present invention, and FIG. For example, it shows the construction of a standard model of pollutant behavior considering atmospheric diffusion.

As shown in FIGS. 1 to 5, the method for constructing a standard model for contaminant behavior according to an embodiment of the present invention can estimate the behavior of pollutants (eg, pesticides, gases, chemicals, etc.) In order to develop a multi-media model composed of the atmosphere, water, vegetation, and soil, the first to fourth steps are carried out.

In the first stage, the environmental behavior mechanism, which is divided into intra-media processes and inter-media processes considering the target medium of the model, such as air, soil, vegetation, and water, It will be done.

In other words, the atmosphere, water, and soil that constitute the model area are divided into four large categories such as air, water, soil, and vegetation. The soil is divided into bare soil and vegetation covered soil.

Here, the input parameters of the environmental behavior mechanisms considered in the model are classified into intra-media processes and inter-media processes, and the input parameters of the media are the input / advection and wind drift, degradation of pollutants and leaching of pollutants, and the input parameters of the inter-media mechanism include volatilization of pollutants, pollutants Dry deposition of pollutants, wet deposition of pollutants, run-off of pollutants, uptake of pollutants to vegetation roots, Wash-off by rainfall, and molecular diffusion of gaseous matter from the atmosphere into the soil and water.

Wherein the second step is to measure the amount of change in the behavior of the pollutant according to the time of the environmental behavior mechanism input parameter input from the first step and the third step is to measure the change amount of the pollutant with respect to the mass balance of the behavior change amount model measured from the second step And the fourth step is to construct the modeling of the pollutant behavior standard through the material balance equation composed of the third step.

That is, the present invention establishes a model of an unsteady state that explains the change of the pollutant over time in order to evaluate the reliability of prediction of the inter-medium concentration. To this end, After the measurement of the change in the behavior of the material, the material balance equation including the input / output parameters of the in-medium / medium distribution, movement and deformation mechanism is established as in the third step, Euler method) to create a model program.

Here, although the above model program uses Microsoft Visual C # .Net 2008, it is not necessarily limited to this.

In this case, in the atmospheric lattice of the input variables, the wind caused by wind after the application of the pollutant is stochastically or microscopically expressed, and spreading is carried out independently according to the wind direction and velocity. .

That is, the atmospheric lattice supports the adjustment of the lattice size according to the changed target area or environmental characteristic, while the contaminant moves to the other lattice by the advection and diffusion due to wind by each lattice, Airborne pollutants are to be taken into account by moving to soil, vegetation, and water by dry and wet deposition.

On the other hand, the pollutants that are sprayed into the atmosphere are changed from the atmosphere to the soil in the form of rainfall, and then, when they reach the water, their amount changes depending on the characteristics of the watershed or the physicochemical properties of the chemical. In order to take into account the geographical characteristics, altitude and direction of the water flow in the area, the geographical information system (GIS) is used to take into account the directionality of the dynamic terrain. And the rainfall runoff coefficient for each grid was calculated and applied to the model.

For example, in the modeling of wind, when the wind is applied, the pollutants are not sprayed all over the spraying area after the spraying of the pollutants, and there is a phenomenon that the winds are blown toward the adjacent region and the stream. Is the ratio approach based on the empirical formula based on the existing monitoring data as shown in Fig.

This proportional approach provides the degree of pollutant blowing as a ratio of distance and wind speed, which is usually applied in the European EPPO and FOCUS model, where the terrain and weather of Europe and Korea There is a limit to apply immediately because of the characteristics, but it has the merit of convenient application of model.

That is, the rate of drift by the wind can be expressed by an empirical equation as shown in Equation 1 below,

[Equation 1]

drift (%): percentage of drift

drift_dist: Distance to drift

V: wind speed (m / sec)

T: Temperature (° C)

X _H2O : Relative humidity (%)

When the above formula is applied, it is calculated that the concentration of the spraying point is about 20% at 1 m from the spraying point and diluted 5 times at the spraying point, 1% is drafted at the spraying point of 10 m, .

However, the actual monitoring data were 1% drafted at 1m point and diluted to 100 times, and the actual monitoring data tended to increase as distance increased.

However, this does not measure the maximum atmospheric concentration at the spraying point and the atmospheric concentration for each point distance after spraying, so that it is difficult to compare with the actual data and calculation results, and it is assumed that the total amount of the spraying points are all in the atmosphere in the spray area And its concentration was calculated. The plate test after 1 m of the spray point. Assuming that the detected concentration was the concentration of the atmospheric point after 1 m of the spray point, it was confirmed that 25% was drafted and diluted 4 times.

Secondly, as shown in the attached FIG. 4 and the following equation (2), there is a method of solving the movement amount of the windward wind by applying the meteorological theory, which is usually applied to the US AgDrift model It is based on many factors that represent the meteorological state.

However, although the second method may reflect geographical and meteorological characteristics, many factors must be constituted to apply to the model according to the embodiment of the present invention, It is also true that there is also a force in the following.

&Quot; (2) "

On the other hand, looking against diffusion and convection of air, air diffusion means to be set in conventional Lagrangian based on the model, a particle size distribution of the contaminants uncertainty in the Lagrange model, and the distribution by σ _x, σ _y .

To represent the diffusion of the emitted pollutant distribution, the particles are considered as cylinders with the standard hogs associated with them, and the diffusion is assumed to have a Gaussian distribution.

Since the emission, transport and diffusion of air pollutants in the model are made in lattice units, diffusion needs to be calculated in units of lattice, and since it is assumed that the concentrations of pollutants are constant in one lattice, By increasing the size of the lattice over time, diffusion was expressed.

In the model according to the embodiment of the present invention, the atmospheric advection system considering the geographical information about the air flow in the medium among the input parameters of the medium is taken into consideration. The wind direction and the wind velocity field are implemented on the map, The GIS is used for the advection and diffusion of pollutants moving by the wind, which is based on the data for 30 years, and the average wind speed is calculated for each day.

At the same time, the atmospheric movement considering the turbulent diffusion was realized, and the inflow and outflow of pollutants in the air were constructed at the interval of every minute by applying the advection and diffusion by the wind. After the particles are moved and diffused, they are divided into portions overlapping with other grids, and the ratio of the contaminants in a certain lattice is determined as a ratio to finally calculate the concentration of contaminants flowing into the lattice.

The vertical distribution of pollutants in a grid is constant in the designated atmospheric altitude and only the horizontal diffusion is considered. Vertical movement of pollutants takes into consideration only the phenomenon of pollutants in the air moving to soil, water and vegetation through deposition and adsorption. .

In the model according to the embodiment of the present invention, since the discharge, migration and diffusion of air pollutants are made in a lattice unit, diffusion needs to be calculated in units of lattices, and the concentration of contaminants in one lattice is constant As a result, diffusion is considered by increasing the size of the lattice over time.

As a result, the length L of one side of the diffused grating grows to (L + 1.54? _Y ) as shown in FIG.

The size of the lattice after diffusion is changed to (L + 1.54σ _y ) X (L + 1.54σ _x ), and the horizontal diffusion coefficients σ _x and σ _y depend on the level of Horizon of Briggs Diffusion coefficient calculation method was used.

In the case of the prediction of soil loss using the rainfall runoff MUSLE equation, in the present invention, in order to estimate the soil loss in the watershed caused by the rainfall, in the non-urban area where the soil is exposed to rainfall, MUSLE was applied (Williams, 1995).

The MUSLE was developed to estimate the amount of soil loss by a single rainfall event. Unlike the USLE, the rainfall energy factor is changed to a runoff energy parameter and the soil loss To be calculated.

&Quot; (3) "

The unit of the soil loss (Y) is 1000 kg / day (or metric ton / day), Q _surf is the direct surface runoff (mmH ₂ O / ha), Area _t is the area of the watershed (ha) _P is a flood amount (m ³ / s), K is a soil erodibility factor, LS is a slope length and steepness factor, C is a crop factor, P is a soil conservation factor is a conservation practice factor.

The soil erosion factor, K, is the relative erosion rate (Morgan, 1986), with a slope of 22.13 m and a slope of 5 ° with a standard slope of 5 °, indicating that the sand and organic matter content, soil structure and permeability (Wischmeier, Johnson & Cross, 1971; Morgan, 1986), depending on the soil characteristics.

The slope and slope factor, LS, is used as a single indicator by a combination of slope length and slope, as in Equation (4) below, which is the expected soil loss per unit area from a given slope and length (Wischmeier & Smith, 1978), and is found by the following equation (4).

At this time, when the slope is large, the erosion increases because the water is superimposed and the speed increases, and when the slope is doubled, the erosion increase rate is about 20-40%.

&Quot; (4) "

Here, the unit of L is m, and the unit of S is%.

The vegetation index, C, is expressed as the ratio of soil loss under uncultivated conditions to the conditions under which certain crops grow (Morgan, 1986). Vegetation protects the surface of the soil from the impact of rainfall, reduces soil runoff by rainfall, improves permeability, and reduces runoff velocity to prevent erosion of the soil. Vegetation impacts include crop species, (Morgan, 1986), as well as the time, as well as the time, of change, depending on the state, cultivation type and management factors.

The value of C is about 1 in highland before the crop grows, and 0.1 in lowland or high density regions. C values for each vegetation condition were used with reference to the values given by Wischmeier and Smith (1978) and the US Soil Conservation Bureau.

Finally, P, the soil conservation factor, is the ratio of soil loss to various arable land on the surface. Cultivated land type elements are divided into contour cultivation, strip cropping, and terrace systems, which vary according to the degree of slope. In the absence of preservation behavior, P is 1.0 (Morgan, 1986) and grassland is 0.17 (Ministry of Environment, 1995).

In addition, the amount of pollutants entering the water system is different from the amount held before the rainfall runoff in the watershed, which is expressed as a ratio (Iyo Sang, 1996; Song Dongha, 1999). The definition of the delivery rate is shown in Equation 5 below (Novotny and Olem, 1994).

&Quot; (5) "

In this case, DR is the delivery ratio (%), Y is the amount of pollutant (kg / day) at the watershed outlet point, and A is the amount of pollutant (kg / day) .

In the model according to the embodiment of the present invention, the distance is calculated by calculating the distance between the position of the soil where the rainfall runoff occurs and the exit point of the watershed, and the rate of reaching the runoff rate according to the reach distance is calculated Respectively.

The GIS was used to easily convert the environmental data, which is the basic data of hazardous chemical behavior, into the input data of the water quality model, and the water runoff and the solid runoff were calculated Combined with a multi-media dynamics model. Through this combination, a topographical model similar to the actual terrain can be created in the model, and a model that allows the toxic chemicals to flow through the terrain can be sorted through this terrain.

As a third step, the mass balance equation of the multicomponent dynamics model can be expressed by the following equation (6): i-the compartment ) Of the amount of material change with time.

&Quot; (6) "

Where Vsys represents the volume of the entire system (m ³ ), and Csys and Msys represent the molar concentration (ng / m ³ ) and mass (ng) of contaminants in the system, respectively. Subscript i represents each medium such as atmospheric, water body, soil, and vegetation.

The amount of contaminant change in each medium can be expressed as the input rate of pollutant minus the output rate as shown in Equation (7) below. The influent rates of pollutants include direct emissions (Ei) of pollutants to the medium, reaction production (Si) of pollutants in the medium, inflow (Finp) through atmospheric and water advection from adjacent areas, (Tj, i), and the outflow rate includes outflow (Fout) by advection and removal (Ri) by decomposition and chemical reaction.

&Quot; (7) "

The mass balance equation for each media is summarized.

The subscripts A, W, S, FS, V, X and GW represent the atmospheric, water body, soil, vegetation soils, vegetation, particulate phase and vapor phase, respectively.

First,

Second, soil (Soil)

- Bare soil

- Forest soil

Third,

- Leaves (Vegetation cuticles)

- Vegetation surfaces

Fourth,

The following describes the terms for the above equations.

winddrift = wind drift flux (unit: ng / sec)

E = emission flux (unit: ng / sec) in which the insecticide is sprayed into the atmosphere in the target area

SO = sorption flux (unit: ng / sec) in which the gas phase insecticidal component is absorbed into vegetation or soil

DD = dry deposition flux (unit: ng / sec) in sediments of soil, vegetation surface,

WD = wet deposition flux (unit: ng / sec), which is dissolved or swept away by atmospheric insecticide in the soil, vegetation surface and river.

VO = volatilization flux (unit: ng / sec) due to the difference of concentration of insecticides in soil,

SRO = the amount of solid run-off flux (unit: ng / sec) that is attached to the particle when the particles in the soil are introduced into the soil by rainfall runoff.

WRO = the amount of water run-off flux (unit: ng / sec)

R = amount of reaction (unit: ng / sec) decomposed by chemical reaction

LC = Water in soil is introduced into groundwater by rainfall. The leaching flux (unit: ng / sec)

LT = litterfall flux (unit: ng / sec) of the insecticidal component sprayed on the leaf of a plant,

WO = wash off flux (unit: ng / sec), which is the amount of insecticide sprayed on the leaf of vegetation,

UP = uptake by root (unit: ng / sec) in which the insecticides dissolved in the soil moisture are absorbed by the vegetation roots.

[Experimental Example]

1. Scope of Model

The media of the model considers the four media of air, soil, vegetation, and water, and the target substance was selected as an insecticide such as Deltamethrin. The spatial range includes streams, ditches, and ponds.

a. Stream

The target stream (stream) is first exemplified in the Daejeon Metropolitan City Rock Stone. The rock was used as a standard target area including river (200m X 1m) and river side (200m X 14m). The spatial extent of this region prediction is divided into a spray area where the spraying is performed, an adjacent area between the spray area and the river, and a river area, as shown in FIG. 6 attached hereto.

Because the river side of a thousand rocks is larger than 7m, which is the actual spraying distance of the insecticide, the adjacent area between the spray zone and the stream was formed.

b. Ditch

Daejeon Metropolitan City (600m X 1m) ditch was used as the target ditch. 7, the total width of the trench including the trench is 4 m, and the contaminant is directly sprayed on the trench water within 7 m of the spraying.

The spatial composition of the model consisted of soil, vegetation, and water in the spray area, and the process of escape by winds, advection and diffusion into the rest of the area. Since bi-directional spraying is in progress, the bi-directional spraying method is considered together.

c. Pond

It is a pond of a pond of 1.5m X 1m in the Mudong area of Daejeon metropolitan city. As shown in FIG. 8, the pond is sprayed directly with water within 7 m of the total width of the pond including the pond side of 1.5 m. The spatial composition of the model is the same as the ditch and varies in size. Only one direction spraying is carried out.

2. Scenario composition of the model

The scenarios were created as shown in the following Table 1 based on the spraying method actually applied to the target area. We predicted the behavior of Deltamethrin in the environment.

In the actual field, 0.125L / min was sprayed at a speed of 10-15km / hour, and the dilution ratio of the weakener was 100: 1. Based on the information obtained from the field, the application amount was estimated with the dilution ratio of undiluted drug of 15g / L.

The number of spraying times was 4 times in one week interval in June. In addition, the meteorological data input to the model is based on 30 - day average data of Daejeon meteorological office.

In the case of rivers, bi-directional spraying is underway, so this is considered, and the actual model is assumed to be sprayed simultaneously in both directions. For the remaining trenches and ponds, only one-way spraying is considered. Since the actual air spray method was applied, the aerosol-type insecticide component was once sprayed into the atmosphere, without considering the direct introduction of the insecticide component into the soil and vegetation, and the insecticide component was immersed in the soil and vegetation . This is to reflect more of the effects of the movement and diffusion of the atmosphere.

division Item Input value

Volume of application Target substance Deltamethrin Insecticide application rate 60L / 8min Insecticide application amount 50 l / ha Deltamethrin effective dose 8.76 g / ha Dilution ratio of undiluted solution 100: 1 Ingredient content ratio 15g / L Injection particle size 100 to 250 Spray width 7m Spraying time Spray distance (200m) / Vehicle travel speed
Spraying scenario Number of Sprays 4 times / year Spray season About a week in June Simulation period 1 year Weather data Monthly / daily average data for 30 years of Daejeon

3. Input variable composition of the model

Input variables that should be integrated into the model are classified into physical property data and parameters of the target material, topographic data of the target area, and weather data as shown in Table 2 below. The main input data are shown in the following table.

The meteorological data are based on the 30 - year average data of Daejeon metropolitan city meteorological office to take into account monthly change and day change. Month and day of the week.

division Item Input value division Item Input value

Target substance
Property data Molecular Weight 505.21
Topographic data LAI, RAI 2.0, 5.4 Vapor pressure (Pa) 0.000002 Flow rate Monitoring data Solubility (g / m ³ ) 0.002

Weather data Wind direction Month change Henry constant 0.0002 Wind velocity Day change Log Kow 6.2 Precipitation Day change oc (L / Kg) 12038 Temperature Day change Half-life (days) 40 Absolute humidity Day change

Topographic data Cloth material GIS data

parameter Deposition rate Diameter change dependence slope GIS data Volatility rate Temperature change dependence Organic carbon content ratio 0.02 Decomposition rate Half-life dependence Porosity 0.3 Diffusion coefficient Atmospheric stability dependence Soil density (kg / m ³ ) 2600 Adsorption coefficient Property dependence

4. Construction of models

When the insecticide delthamethrin was sprayed on the target area, a multi - media model was constructed to predict the multi - media behavior of diffusion from the atmosphere to water. The model includes the process of diffusion by wind and the process of movement between media. We have developed a model that solves the material problem including Euler numerical analysis method (Euler method) including the distribution in the medium and the medium, and the moving device. The model utilized Visual Basic 2009.

It was also constructed as an unsteady state model that describes the change of pollutants over time after spraying. The model is calculated in units of 1 second, and the output unit of the result of the model is output in 1 minute for the atmosphere with large change in concentration over time, and 1 hour for the soil with relatively low concentration change. We also constructed the input variables described above that should be reflected in the mass balance equation.

In order to reflect some of the physical properties and parameters of the input variables as the inherent data depending on the material, other insecticidal components can be predicted in connection with the database constructed in the general task. Thus, in the present invention, Since the database is constructed and integrated with the topographic data, it is possible to predict the contamination of other insecticides in the database.

In addition, since the model according to the embodiment of the present invention can input other information about the terrain data, information about the application amount, and the main information about the river, other values than the default value applied in the actual standard area in the basic screen, It can be confirmed that generalized model can be applied by inputting environmental conditions.

While the present invention has been described with reference to the accompanying drawings and accompanying drawings, it is to be understood that the invention is not limited to the disclosed exemplary embodiments.

It will be apparent to those skilled in the art that various modifications and variations can be made in the present invention without departing from the spirit and scope of the invention as defined by the appended claims. It is to be understood that such changes and modifications are within the scope of the claims.

Claims

The first step is to derive the input parameters of the environmental behavior mechanism, which is divided into intra-media processes and inter-media processes considering the target media of the model, air, soil, vegetation and water.
A second step in which the environmental behavior mechanism derived from the first step measures the amount of change in the behavior of contaminants with time from input variables;
A third step of constructing a mass balance equation for the amount of change in behavior measured from the second step; And
A fourth step of establishing modeling of the pollutant behavior standard through the material balance equation composed of the third step; Wherein the method comprises constructing a model of a pollutant behavior standard model.

The method according to claim 1, wherein the pollutant is one of an insecticide, a gas, and a chemical substance that is diffused into the atmosphere and pollutes the atmospheric, atmospheric, soil, vegetation, and water of a target medium.

The method of claim 1, wherein the input parameters of the media include advection and wind drift of applied contaminants, degradation of contaminants, leaching of contaminants, A method for constructing a standard model of pollutant behavior.

The method of claim 1, wherein the input parameters of the inter-media mechanism are selected from the group consisting of volatilization of pollutants, dry deposition of contaminants, wet deposition of contaminants, run- off uptake of contaminants to vegetation roots, wash-off by rainfall on vegetation surfaces to contaminants, and diffusion of gaseous matter from the atmosphere into the soil and water. A method for constructing a standard model of pollutant behavior.

The method according to claim 3, wherein the amount of drift (%) of the wind-induced pollutant flow is calculated by a numerical method using empirical equations or a meteorological theory, Wherein the pollutant concentration is calculated by the following formula.

4. The method according to claim 3, wherein the advection of the pollutant is performed by checking the inflow and outflow of the pollutant diffused into the atmosphere so as to apply advection and diffusion by the wind at predetermined time intervals. How to build.

The method of claim 1, wherein the model calculates the distance between the location of the soil where the rainfall runoff occurs and the exit point of the watershed, calculates the reach distance, and calculates the rate of reaching the runoff rate according to the reach distance according to the model's simulation time Wherein the method comprises the steps of:

The method according to claim 1, wherein the input parameters of the environmental behavior mechanism, which is the basic data of the pollutant behavior, are configured to utilize a geographic information system (GIS).

The method according to claim 1, wherein in the third step, the mass balance equation for the medium kinetics model is calculated by multiplying the amount of change in the pollutant of the entire system by the sum of the mass change amounts with time in each medium (i-the compartment) Wherein the method comprises the steps of:

10. The method according to claim 9, wherein the change amount of the pollutant in the medium is expressed by a value obtained by subtracting the output rate from the input rate of the pollutant, How to build.

The method according to claim 1, wherein the mass balance equation for the air is expressed by the following equation.

The method according to claim 1, wherein the mass balance equation for the soil in the soil is constructed as shown in the following equation.

2. The method according to claim 1, wherein the mass balance equation for the soil in the soil is constructed as shown in the following equation.

2. The method according to claim 1, wherein the mass balance equation for the vegetation cuts in the vegetation is expressed by the following equation.

The method according to claim 1, wherein the material balance equation for the vegetation surfaces in the vegetation is expressed by the following equation.

The method according to claim 1, wherein the mass balance equation for the water is expressed by the following equation.

The method according to any one of claims 11 to 16, wherein the model constructed by the mass balance equation in the fourth step is a model program that solves the Euler method. Model building method.