KR102018281B1 - Method and system for analysing behavior of hazardous chemical - Google Patents

Abstract

According to an aspect of the present invention, a hazardous chemical behavior analysis system driven in an electronic device comprises a memory for storing a hazardous chemical behavior analysis program and a processor for executing the program. When the program is executed, the processor provides a user interface for setting leak information on leak locations, leak chemicals, and leak timing, and container information on the size of a container in which the leak chemicals are stored, the container pressure, and the leak height of the container, generates wind fields for a surrounding area of the leak location based on the leak information, and terrain information and weather observation information obtained from a server, calculates a leak rate of the leak chemical based on the container information, and performs modeling for a diffusion state of the leaked chemical based on the wind field and the calculated leak rate, wherein the leak rate is calculated repeatedly until a liquid level is lowered below the leak height of the container.

Description

Hazardous Chemical Behavior Analysis and System {METHOD AND SYSTEM FOR ANALYSING BEHAVIOR OF HAZARDOUS CHEMICAL}

A method and system for analyzing the behavior of hazardous chemicals.

More than 15 million chemicals are used commercially around the world, with about 200 to 1,000 chemicals being produced in more than 1 tonne per year. About 43,000 chemicals are distributed in Korea. Chemicals can be exposed to humans and the environment through various routes throughout the manufacturing, use, and disposal process.In particular, chemical accidents have the potential to cause catastrophic disasters, requiring thorough management and response measures. Do. However, since the ministries that manage these hazardous chemicals are operated independently, it is difficult to operate an efficient system through the organic relationship of disaster management.

In this regard, the wind direction and wind speed of the area where a chemical leakage accident occurs by driving a meteorological diagnosis model is performed, and then the air diffusion model is driven to apply the diagnosed wind direction and wind speed and chemical information to the driven air diffusion module. A method for predicting the spread of is proposed. However, since the air diffusion model takes a considerable time to drive the weather diagnosis model, there is a problem in that an immediate chemical leakage accident does not occur immediately.

Therefore, it is necessary to introduce a unified system that can unify related information of related ministries and pursue effectiveness of corresponding information for quick initial response.

Republic of Korea Patent Application Publication No. 10-2016-0116675 (Name of the invention: Leak incident response system and method)

The present invention is to solve the above-mentioned problems of the prior art, to determine the topographic information and weather observation information based on the information on the harmful chemicals and the location of the leak to predict the wind field for the surrounding area of the leak location The aim is to provide a diffusion result of the leaking chemical based on the leak rate from the storage container in which the leaking chemical is stored.

However, the technical problem to be achieved by the present embodiment is not limited to the technical problem as described above, and other technical problems may exist.

As a technical means for achieving the above technical problem, the hazardous chemical behavior analysis system driven in the electronic device according to an aspect of the present invention includes a memory for storing the hazardous chemical behavior analysis program and a processor for executing the program; The processor provides a user interface that, as the program runs, sets leak information for leak locations, leak chemicals, and timing of leaks, as well as container information for the container size, container size, container pressure, and leak height of the container where the leak chemicals are stored. Generate a wind field for the surrounding area of the leak location based on the leak information and the terrain information and weather observation information obtained from the server, calculate the leak rate of the leaked chemical based on the container information, Model the diffusion of leaking chemicals based on the leak rate At the same time, the leak rate is repeated until the liquid level drops below the leak height of the vessel.

In addition, the operating method of the hazardous chemical behavior analysis system according to another aspect of the present invention is the leakage information on the leak location, the leak chemical and the timing of the leak and the size of the container in which the leak chemical is stored, the container pressure and the leak height of the container Providing a user interface for setting container information; Generating a wind field for the surrounding area of the leak location based on the leak information and the terrain information and the weather observation information obtained from the server; Calculating a leak rate of the leaked chemical based on the container information; And modeling the diffusion state of the leaking chemical based on the wind field and the calculated leak rate, wherein calculating the leak rate repeatedly calculates the leak rate until the liquid level drops below the leak height of the vessel.

According to the above-described problem solving means of the present invention, the hazardous chemical behavior analysis system according to an embodiment of the present invention and the terrain information obtained from the server as the leak location and the leaked chemicals (that is, leak information) and Intuitive and unifying of leak location and leaked chemicals for users accessing the system by predicting wind fields for the surrounding area of the leak location based on weather observation information and providing a real-time visualization of the diffusion of the leaked material Information can be provided.

In addition, by calculating and visually providing not only the leak rate of the container in which the leaked chemical is stored but also the full evaporation rate of the leaked chemical introduced into the ground, it is possible to intuitively grasp the leakage result of the harmful chemical over time.

1 is a block diagram showing the configuration of a hazardous chemical substance behavior analysis system according to an embodiment of the present invention.
2A is a view for explaining that the wind field generation according to an embodiment of the present invention is divided into an observation region and an unobserved region.
Figure 2b is a diagram showing each parameter affecting the wind (wind direction / wind speed) according to the terrain conditions according to an embodiment of the present invention.
3 is a view for explaining the diffusion state of the leak chemical in accordance with an embodiment of the present invention.
4A illustrates an example of a screen on which a leak rate calculation result of a container according to an embodiment of the present invention is output.
Figure 4b is an example showing a screen on which the evaporation rate calculation result of the pool according to an embodiment of the present invention is output.
5 is an example showing a screen output from the display unit as the hazardous chemical behavior analysis process is performed according to an embodiment of the present invention.
FIG. 6 is a flowchart illustrating an operation method of performing a hazardous chemical behavior analysis process by the processor of FIG. 1 according to an exemplary embodiment of the present invention.

DETAILED DESCRIPTION Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those skilled in the art may easily implement the present invention. As those skilled in the art would realize, the described embodiments may be modified in various different ways, all without departing from the spirit or scope of the present invention. In the drawings, parts irrelevant to the description are omitted in order to clearly describe the present invention, and like reference numerals designate like parts throughout the specification.

Throughout the specification, when a part is "connected" to another part, this includes not only "directly connected" but also "electrically connected" with another element in between. . In addition, when a part is said to "include" a certain component, which means that it may further include other components, except to exclude other components, unless specifically stated otherwise, one or more other features It is to be understood that the present disclosure does not exclude the possibility of the presence or the addition of numbers, steps, operations, components, parts or combinations thereof.

In the present specification, the term 'unit' includes a unit realized by hardware, a unit realized by software, and a unit realized by both. In addition, one unit may be realized using two or more pieces of hardware, and two or more units may be realized by one piece of hardware. Meanwhile, '~' is not limited to software or hardware, and '~' may be configured to be in an addressable storage medium or may be configured to reproduce one or more processors. Thus, as an example, '~' means components such as software components, object-oriented software components, class components, and task components, and processes, functions, properties, procedures, and the like. Subroutines, segments of program code, drivers, firmware, microcode, circuits, data, databases, data structures, tables, arrays, and variables. The functionality provided within the components and the 'parts' may be combined into a smaller number of components and the 'parts' or further separated into additional components and the 'parts'. In addition, the components and '~' may be implemented to play one or more CPUs in the device or secure multimedia card.

Hereinafter, an embodiment of the present invention will be described in detail with reference to the accompanying drawings.

1 is a block diagram showing the configuration of a hazardous chemical substance behavior analysis system according to an embodiment of the present invention.

Referring to FIG. 1, the hazardous chemical behavior analysis system 100 of the present invention includes a memory 110 for storing a hazardous chemical behavior analysis program and a processor 130 for executing the program. As the program runs, it provides a user interface for setting leak information on leak locations, leak chemicals, and timing of leaks, as well as container information on the container size, container pressure, and the leak height of the container where the leak chemicals are stored. Generate wind fields for the surrounding area of the leak location based on the information and the topographical information and weather observation information obtained from the server 10, calculate the leak rate of the leaked chemicals based on the container information, Model the diffusion of the leaking chemicals based on the leak rate, but calculate the leak rate until the liquid level drops below the leak level of the vessel. It can be calculated repeatedly.

Here, the map information on the leak location of the leak information may be obtained from an external geographic information providing server (not shown) as GIS (Geographic Information System) information. For example, the leak location may be composed of a leak point (x, y) and a leak height (m) which will be described later with reference to FIG. 5, where the leak height means a height on a terrain or a building where a container is located. The leak height of one container and the other. In addition, the terrain information and weather observation information may be obtained from an external server 10 such as a geographic information providing server and a weather observation server, and the terrain information includes land use topology information including land use and altitude.

Specifically, the hazardous chemical behavior analysis system 100 includes a memory 110, a communication unit 120, and a processor (or control unit) 130, such as a desktop computer, a laptop computer, a server, a portable terminal, and the like. It may be driven in an electronic device having a communication function.

The memory 110 stores a hazardous chemical behavior analysis program (or at least one instruction). In this case, the memory 110 refers to a nonvolatile storage device that maintains stored information even when power is not supplied, and a volatile storage device that requires power to maintain stored information.

The communicator 120 may include one or more components for communicating with the outside, such as the server 10 providing terrain information and weather observation information. For example, the communication unit 120 performs a local area network (LAN) or a Bluetooth (Bluetooth), BLE (Bluetooth Low Energy), infrared (IrDA), or Zigbee communication that performs wired or wireless data communication. It may include a component capable of performing at least one of. Alternatively, the communication unit 120 may transmit and receive a wireless signal through a mobile communication network or a broadcast communication network.

The processor 130 controls the overall operation of the hazardous chemical behavior analysis system 100. To this end, the processor 130 may include at least one of a random access memory (RAM), a read-only memory (ROM), a CPU, a graphic processing unit (GPU), and a bus.

In addition, the processor 130 executes a program stored in the memory 110 to perform a process of analyzing harmful chemical behavior.

2A is a view for explaining that the wind field generation according to an embodiment of the present invention is divided into an observation region and an unobserved region. Figure 2b is a diagram showing each parameter affecting the wind (wind direction / wind speed) according to the terrain conditions according to an embodiment of the present invention. 3 is a view for explaining the diffusion state of the leak chemical in accordance with an embodiment of the present invention.

The processor 130 provides a user interface for setting leak information on the leak location, the leak chemical and the timing of the leak, and the container information on the container size, the container pressure and the leak height of the container in which the leak chemical is stored, the leak information, Based on the information obtained from the server 10 providing vessel information, terrain information and weather observation information, a wind field is generated for the surrounding area of the leak location, and the liquid level is lowered below the leak height of the container based on the container information. Iterate over to calculate the leak rate of the leaking chemical, and model the diffusional state of the leaking chemical based on the wind field and the calculated leak rate. Here, the leak rate of the leaked chemical includes not only the leak rate of flammable liquid leaked from the container but also the evaporation rate when a pool is formed on the surface, and a detailed description of the user interface for outputting container information is shown in FIG. 4A. And it will be described later with reference to Figure 4b. In addition, the wind field means that the spatial distribution of the wind speed / wind direction is represented by the position function, and a detailed description of the user interface for outputting leak information will be described later with reference to FIG. 5.

Referring to FIG. 2A, the processor 130 divides the surrounding area of the leak location into a grid having a predetermined size when generating the wind field, and divides the grid divided area into the observation area 322 and the unobserved area 321. The atmospheric information of the unobserved region 321 may be predicted by using the terrain information and the weather observation information of the observation region 322 obtained from the server 10 and the terrain information of the non-observation region 321. Here, the grid may be set in meters (m) to have an area of square meters (m ² ).

The processor 130 may predict the wind direction and the wind speed of the unobserved region based on at least one or more of kinematic effects, slope flow, and terrain conditions when predicting atmospheric information of the unobserved region. Here, the kinematic effect refers to the amount of wind direction / wind speed change generated according to the ground height (altitude). Slop flow refers to the amount of wind direction / wind speed change caused by the change of the boundary layer structure in addition to the simple change of direction as the terrain conditions change. Blocking Effect refers to the amount of wind direction / wind change that occurs according to the use of land, such as lakes, buildings, grasslands, for example.

For example, the processor 130 may use the terrain information (land use and altitude) of the observation area 322 and the weather observation information (atmosphere information) and the terrain information of the non-observation area 321. The atmospheric information of 321 is predicted. In this case, the processor 130 may predict the atmospheric information including the wind direction / wind speed of the unobserved region 321 based on the kinematic influence of the wind, the slope flow, the terrain condition (that is, the lake, the building, the meadow, etc.). have.

Hereinafter, the equations according to kinematic influences, slopes and topographic conditions will be described in detail.

Kinematic influence is calculated by the following equation (1).

<Equation 1>

V stands for domain-mean wind, h _t stands for terrain height, k stands for stability-dependent coefficient of exponential decay, z stands for vertical Represents a vertical coordinate.

In addition, the inclined flow is calculated by the following expression (2).

<Equation 2>

S _e represents the equilibrium speed of the slope flow, L _e represents the equilibrium length scale, and x represents the distance to the crest of the hill. Δθ represents the potential temperature deficit with respect to the environment, θ represents the potential temperature of the environment, C _D represents the surface drag coefficient, h Denotes the depth of the slope flow, α denotes the angle of the terrain relative to the horizontal, and k denotes the entrainment coefficient at the top of the slope flow layer. And g represents gravitational acceleration (9.8 m / s ² ).

Figure 2b shows each parameter of the impact on the wind (wind direction / wind speed) according to the topographic conditions (land use), for example, the surface roughness according to the land use of the city, farmland, lakes, wetlands, wasteland, etc. Respective parameter values for Surface Roughness, Albedo, Bowen Raito, Soil Heat Flux Parameter, Anthropogenic Heat Flux, and Leaf Area Index Indicates.

Next, the processor 130 may calculate the concentration of the leaking chemical in the wind field when modeling the diffusion state of the leaking chemical.

In exemplary embodiments, the processor 130 calculates a leak rate (kg / sec) of leaked chemicals, that is, a leak amount per second (kg), and generates a particle diffusion model based on the calculated leak rate and wind field. It can be applied to model the diffusion state of leaking chemicals (ie, chemical concentrations in the wind field).

Equation 3 below represents a particle diffusion model.

<Equation 3>

C _p represents the concentration of the leaking chemical (kg / m ³ ) at time t, position (x, y, z), H _E represents the effective height, and Q represents the leakage. Further, σ _x is the diffusion coefficient in the wind direction (ie, wind direction), σ _y is the diffusion coefficient in the wind cross section direction, and σ _z is the diffusion coefficient in the wind cross section direction. In addition, u represents the wind speed at the actual leak height (H _s ), H _m is the mixing height, t is the time measured after the chemical leaks.

In this case, the particle diffusion model does not consider the change in the wind field over time. Through this, the processor 130 may model the diffusion state of the leaked chemical in real time. Here, the real time is a time within several minutes after the incident scenario is input, and may be a threshold time for initial response to the leaked chemical.

FIG. 3 is an example showing the diffusion states 410, 420 of the leaking chemical obtained through the particle diffusion model on the wind field 302 for the area around the leak location. As shown in FIG. 3, the processor 130 may predict the diffusion state over time after the chemical leakage through the particle diffusion model. Diffusion concentrations over time can be calculated by summing the concentrations of leaking chemicals for each lattice of the wind field 302. For example, the leaked particles are diffused 410 according to the first leak, and the leaked particles are spaced at a time (preferably 1 second) from the leaked particles according to the second leak, and the first and second leaks are diffused. Has independent concentration distribution by the particle diffusion model described above. Because of this, the total concentration of the leaking chemical may consist of the sum 420 of the primary and secondary leak concentrations.

In addition, the toxic chemical behavior analysis system 100 may optionally include a display unit 140. The display unit 140 receives an image signal processed by the processor 130 and outputs the image signal to the display area, and may be implemented as a PDP, an LCD, an OLED, a flexible display, or a 3D display.

On the other hand, when the toxic chemical behavior analysis system 100 does not include the display unit 140, the image signal processed by the processor 130 is provided to the external display device through the communication unit 120, the external display device It can also be output.

4A illustrates an example of a screen on which a leak rate calculation result of a container according to an embodiment of the present invention is output.

Figure 4b is an example showing a screen on which the evaporation rate calculation result of the pool according to an embodiment of the present invention is output.

The processor 130 may repeatedly calculate the leak rate of the leaked chemical material until the liquid level is lowered below the leak height of the container based on the container information.

Specifically, the leak rate (kg / sec) of the container is calculated based on the previously set container information. Then, the density of the changed container is calculated as the liquid level inside the container is lowered. The vessel temperature and vessel pressure can then be recalculated from the changed density and the leak rate can be calculated repeatedly until the liquid level drops below the vessel's leak height based on the calculated temperature and pressure again. In addition, referring to FIG. 4A, the container information includes a container size (volume of a container), a container pressure (pressure inside a container), and a leak height of the container, and may be set by a user, but is not limited thereto. For example, the leak chemical, the liquid level, the leak area (leak diameter on the container) and the storage temperature (temperature) may be automatically input or set by a user according to the set value.

The processor 130 may calculate the leak rate of the container by the following Equation 4 based on the density, temperature, and pressure change amount of the container according to the change of the liquid level.

<Equation 4>

, V는 용기의 부피,

은 누출률, T_o는 용기의 초기온도, T₁은 용기의 변화된 온도, P_o는 용기의 초기 압력, P₁은 용기의 변화된 압력, ρ_o는 용기의 초기밀도, ρ₁은 용기의 변화된 밀도를 나타낸다. Where C _d is the leakage coefficient, A _h is the leakage area, P ₀ is the internal pressure of the container, T is the storage temperature, C _p is the static specific heat of the material, C _v is the static specific heat of the material, and Γ is the specific heat ratio C _v / C _p , R is gas constant, K is

, V is the volume of the vessel,

Is the leak rate, T _o is the vessel's initial temperature, T ₁ is the vessel's changed temperature, P _o is the vessel's initial pressure, P ₁ is the vessel's changed pressure, ρ _o is the vessel's initial density, ρ ₁ is the vessel's varied Indicates density.

When the leak rate is calculated, the processor 130 calculates the surface area of the pool and calculates the evaporation rate of the pool based on the surface area of the pool when the leaked chemical leaks to the ground to form a pool. It may be repeated until the diameter does not change.

Specifically, referring first to FIG. 4B, the surface area of the pool formed on the ground surface is calculated based on the diameter of the pool, and the evaporation rate of the pool is calculated. Subsequently, the evaporation rate of the pool can be calculated again from the surface area of the changed pool, and the evaporation rate of the pool can be calculated by repeating until the diameter of the pool formed on the ground surface does not change (when all materials evaporate).

For example, the pool formation model is a model that can be calculated for the expansion of the pool and the evaporation of the liquid substance on the ground or the water surface, and in case of leakage on the ground, heat conduction from the ground, convection from the atmosphere, and radiant heat And reflects on the diffusion of steam. These factors are the main mechanisms for boiling and evaporation.

It also takes into account the possible reaction and dissolution of liquid substances in case of spillage on the water surface, which is important for some chemicals. These effects can be modeled numerically by maintaining the heat and mass balance equations for the boiling or evaporating pool. For example, when boiling on the ground, the heat conduction from the ground can be modeled assuming a uniform semi-infinite medium in which the grass expands, and has a fixed radius (the grass touches the wall of the effluent or has a minimum thickness). Evaporation rate can be calculated by considering parameters (coarse factor, thermal conductivity, thermal diffusivity) and parameters according to the type of surface and surface type (wet soil, dry soil, concrete, insulating concrete, etc.). When evaporated, it can be modeled by taking into account the factors that affect the pool as a result of convection on the ground or surface, solar radiation, boiling on the surface, and evaporation from the surface.

The temperature of the pool is calculated by considering both the heat absorbed by the leaked liquid or the heat lost during evaporation. The principle of pool evaporation can be explained by the law of energy conservation, and the maximum energy that a liquid pool can have is its energy at boiling point. In other words, energy from outside (surface, sun, atmosphere, etc.) cannot accumulate continuously in the liquid phase, and thus, to conserve energy, it releases energy coming back through a mechanism called evaporation.

The evaporation rate of the pool is calculated by the following expression (5).

<Equation 5>

Specifically, at each time step of the numerical calculation, the heat balance of the pool is evaluated and the evaporation mechanism, evaporation rate and pool temperature are determined. The initial pool temperature, T _pool, is assumed to be equal to the T _spill temperature, the temperature of the leaked liquid. The total heat flow of the pool is calculated as shown in equation (1).

Where Q _cond (W) is heat flux from conduction, Q _conv (W) is heat flux from convection, Q _rad (W) is heat flux from solar radiation (W), Q _sol (W) is heat flux from melt, Q _spill (W ) Is the heat flux from the leaked liquid and Q _evap (W) is the heat flux from the evaporation.

If there is a leak on the ground, the heat flux Q _sol due to the melting furnace is zero, and for continuous and time-varying leaks, the Q _spill is calculated as shown in Eq. (2).

C _pL (T _spill ) is the heat capacity of the liquid at the time of leakage (J / kgK), and Q _spill is zero if the leak ends or is a transient leak. The change in pool temperature is calculated from the total heat flux entering the pool.

In equation (3), M _pool is the mass of the liquid pool at the time of calculation. This equation gives the pool temperature as a function of time and cannot exceed the boiling point of the liquid. If the temperature of the pool reaches the boiling point, the liquid boils, and the heat flux entering the pool at this time is given by equation (4), and the evaporation rate is given by equation (5).

ΔHv is the heat of evaporation (J / kg). If the temperature of the pool is lower than the boiling point, it is assumed that evaporation occurs at that temperature, and the evaporation rate is calculated as shown in Eq. (6).

4A and 4B, the user interface 300 includes an image output area 210 for outputting a simulation image and a first chart area 221 and a second chart area 222 for providing a simulation result. Can be configured. Here, the first and second chart areas 221 and 222 may be configured to be switched in the same area by the container liquid leakage and the pool evaporation identification information.

For example, the image output area 210 may be located at an upper area of the user interface 300 to display container information and result values according to a user setting, and may include a process of leaking hazardous chemicals (evaporation of a leak / pool of a container). ) Can be used to output an image that simulates. In addition, the user interface 300 may be configured to include a button for setting play, stop, section movement, and output speed of the simulation image. For example, when a leak simulation image is reproduced, a material (leak chemical), a tank size (a container size), a leak location (a leak height of a container), and an air temperature (a container temperature) are displayed on a screen of the image output area 210. The container information and the resulting value including the elapsed time of accident occurrence, cumulative evaporation amount, leakage amount and the pool diameter can be displayed in conjunction with the elapsed time of accident occurrence.

In addition, as shown in Figure 4a, the first chart area 221 is located in the lower area of the user interface 300 and the liquid level (m), leakage of the container over time (s) with respect to the liquid leakage of the container Results may include the rate (kg / s), temperature (° C.) and pressure (atm). Similarly, as shown in FIG. 4B, the second chart area 222 is located in the same area as the first chart area 221 and the pool diameter m and the pool over time (s) with respect to the evaporation of the pool. The results including the evaporation rate (kg / s) can be displayed. For example, the first and second chart areas 221 and 222 may display the result in the form of a chart or graph.

Thus, by calculating and visually providing not only the leak rate of the leaked chemicals but also the full evaporation rate of the leaked chemicals introduced into the ground, it is possible to intuitively grasp the leakage of the harmful chemicals over time.

5 is an example showing a screen output from the display unit as the hazardous chemical behavior analysis process is performed according to an embodiment of the present invention.

Referring to FIG. 5, the output user interface 300 includes a first region 310 provided to input leakage information, each wind field (spatial distribution of wind velocity / wind direction as a position function) and concentration distribution in a grid form. It may be configured to include a second area 320 provided with information.

For example, as shown in FIG. 5, through the configuration of the user interface 300, the hazardous chemical behavior analysis system viewed through the center of longitude and latitude of the region of interest, the target range (Km), and the observation time (min) Through this, you can set the region of interest information that you want to predict the diffusion state of the modeled leaked material.

Then enter the leak location, such as leak point (x, y), leak height (m), and specify the leak chemical (e.g. chlorine) and the leak point (e.g. year, month, day, time). You can enter Meanwhile, the size of the grid in the user interface 300 may be changed according to a user input through the diffusion model driving grid (m) item, and specific weather information (wind speed and wind direction) may be applied.

Thus, the system according to the disclosed embodiment sets the wind field for the surrounding area of the leak location based on topographical information and weather observation information obtained from the server as it sets the leak location, leak chemicals and the time of the leak (ie, leak information). By predicting and providing real-time visualization of the diffusion of leaked materials, users can access the system with intuitive, unified information about leaked locations and leaked chemicals.

In the case of the configuration performing the same function among the components shown in FIGS. 1 to 5 described above, description thereof will be omitted.

FIG. 6 is a flowchart illustrating an operation method of performing a hazardous chemical behavior analysis process by the processor of FIG. 1 according to an exemplary embodiment of the present invention.

The operation method of the hazardous chemical behavior analysis system of the present invention is a user who sets the leak information on the leak location, the leak chemical and the time of the leak and the container information on the container size, the container pressure and the leak height of the container in which the leak chemical is stored. Providing an interface (S110), generating a wind field for the surrounding area of the leak location based on the information obtained from the server 10 providing the leak information and terrain information and weather observation information (S120), the container Calculating a leak rate of the leaked chemical based on the information (S130) and modeling a diffusion state of the leaked chemical based on the wind field and the calculated leak rate (S140), and calculating the leak rate. Step S130 may repeat the calculation of the leak rate until the liquid level is lowered below the leak height of the container.

Generating the wind field (S120) is a step of dividing the surrounding area of the leak location into a grid of a predetermined size, dividing the grid divided area into an observation area and an unobserved area, and obtained from the server 10 The method may include predicting atmospheric information of the unobserved region by using the topographic information of the observation region and the weather observation information and the topographic information of the non-observation region.

Predicting atmospheric information of the unobserved region may include predicting wind direction and wind speed of the unobserved region based on at least one of kinematic effects, slope flow, and terrain conditions.

Calculating the leak rate of the leaked chemical (S130) is based on the container information, the unit time of the leaked chemical based on at least one of the storage amount of the leaked chemical of the container, the leak diameter (area) of the container and the storage container pressure The amount of sugar leak can be calculated. In addition, step S130 includes calculating the surface area of the pool, and calculating the evaporation rate of the pool based on the surface area of the pool, when the leaked chemicals leak to the ground to form a pool. The calculating step may be repeated until the diameter of the pool does not change.

Modeling the diffusion state of the leaked chemical (S140) may include calculating the concentration of the leaked chemical in the wind field.

Meanwhile, an embodiment of the present invention may also be implemented in the form of a recording medium including instructions executable by a computer, such as a program module executed by the computer. Computer readable media can be any available media that can be accessed by a computer and includes both volatile and nonvolatile media, removable and non-removable media. In addition, the computer readable medium may include a computer storage medium. Computer storage media includes both volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data.

Although the methods and systems of the present invention have been described in connection with specific embodiments, some or all of their components or operations may be implemented using a computer system having a general purpose hardware architecture.

The foregoing description of the present invention is intended for illustration, and it will be understood by those skilled in the art that the present invention may be easily modified in other specific forms without changing the technical spirit or essential features of the present invention. will be. Therefore, it should be understood that the embodiments described above are exemplary in all respects and not restrictive. For example, each component described as a single type may be implemented in a distributed manner, and similarly, components described as distributed may be implemented in a combined form.

The scope of the present invention is shown by the following claims rather than the above description, and all changes or modifications derived from the meaning and scope of the claims and their equivalents should be construed as being included in the scope of the present invention. do.

100: Hazardous Chemical Behavior Analysis System
110: memory
120: communication unit
130: processor
140: display unit

Claims (12)

In a hazardous chemical behavior analysis system driven by an electronic device,
A memory for storing the hazardous chemical behavior analysis program; And
Including; processor for executing the program;
The processor, as the program is executed,
Providing a user interface for setting leak information on leak locations, leak chemicals and timing of leaks and container information on the container size, container pressure and leak height of the container where the leak chemicals are stored,
Generate a wind field for the surrounding area of the leak location based on the leak information and terrain information and weather observation information obtained from the server, calculate a leak rate of the leak chemical based on the container information, Model the diffusion state of the leaking chemical based on the field and the calculated leak rate,
The leakage rate of the container is calculated by the following Equation 1 based on the density, temperature, and pressure change of the container according to the change of the liquid level.
When calculating the leak rate is to repeat the calculation until the liquid level is lower than the leak height of the container,
When the leaking chemical leaks to the ground surface to form a pool when calculating the leak rate,
Calculating the surface area of the pool based on the diameter of the pool and calculating the evaporation rate of the pool based on the surface area of the pool,
It is to calculate the evaporation rate of the pool based on the heat balance, evaporation rate and temperature of the pool over time, iterate repeatedly until the diameter of the pool does not change,
The user interface is
An image output area for outputting a simulation image and a first chart area and a second chart area for providing a simulation result,
The image output area is located in the upper area of the user interface and displays the container information, the diameter and the resultant value of the pool according to a user setting, and the leakage of the hazardous chemicals including the leak of the container and the evaporation of the pool. Outputs a simulation image of the process, and is configured to include buttons for setting the play, stop, section movement and output speed of the simulated image,
The first and second chart areas are positioned in the lower area of the user interface but are configured to be switched in the same area.
The first chart area is to display a result including the liquid level, leak rate, temperature and pressure of the container over time with respect to the leak of the container,
Wherein the second chart area displays a result including the diameter and evaporation rate of the pool over time for the evaporation of the pool.
[Equation 1]

Where C _d is the leak coefficient, A _h is the leak area, P ₀ is the pressure inside the vessel, T is the storage temperature, C _p is the static specific heat of the material, C _v is the static specific heat of the material, and Γ is the specific heat ratio C _v / C _p , R is gas constant, K is

, V is the volume of the vessel,

Is the leak rate, T _o is the vessel's initial temperature, T ₁ is the vessel's changed temperature, P _o is the vessel's initial pressure, P ₁ is the vessel's changed pressure, ρ _o is the vessel's initial density, ρ ₁ is the vessel's varied density)
The method of claim 1,
The processor is
When generating the wind field, the surrounding area of the leak location is divided into a grid having a predetermined size, and the area divided by the grid is divided into an observation area and an unobserved area,
To predict the atmospheric information of the unobserved region using the terrain information and weather observation information of the observation area and the observation information obtained from the server, hazardous chemicals behavior analysis system.
The method of claim 2,
The processor is
And when predicting atmospheric information of the unobserved region, predicting the wind direction and the wind speed of the unobserved region based on at least one of kinematic effects, slope flow and terrain conditions.
delete delete The method of claim 1,
The processor is
To calculate the concentration of leaking chemicals in the wind field when modeling the diffusion state of the leaking chemicals, hazardous chemical behavior analysis system.
In the operation method of the toxic chemical behavior analysis system,
Providing a user interface for setting leak information on a leak location, leak chemical and a time of leak and container information on the container size, container pressure and leak height of the container in which the leak chemical is stored;
Generating a wind field for a surrounding area of the leak location based on the leak information, terrain information obtained from the server, and weather observation information;
Calculating a leak rate of the leaked chemical based on the container information; And
Modeling a diffusion state of the leaking chemical based on the wind field and the calculated leak rate,
The step of calculating the leak rate
The leakage rate of the container is calculated by the following Equation 1 based on the density, temperature, and pressure change of the container according to the change of the liquid level.
Repeat the calculation until the liquid level drops below the leak height of the vessel,
If the leaking chemical leaks to the ground and forms a pool,
Calculating a surface area of the pool, and calculating an evaporation rate of the pool based on the surface area of the pool,
Calculating the evaporation rate of the pool
Calculating the evaporation rate of the pool based on the heat balance, evaporation rate and temperature of the pool over time,
Repeat until the diameter of the pool does not change,
The user interface is
An image output area for outputting a simulation image and a first chart area and a second chart area for providing a simulation result,
The image output area is located in the upper area of the user interface and displays the container information, the diameter and the resultant value of the pool according to a user setting, and the leakage of the hazardous chemicals including the leak of the container and the evaporation of the pool. Outputs a simulation image of the process, and is configured to include buttons for setting the play, stop, section movement and output speed of the simulated image,
The first and second chart areas are positioned in the lower area of the user interface but are configured to be switched in the same area.
The first chart area is to display a result including the liquid level, leak rate, temperature and pressure of the container over time with respect to the leak of the container,
And the second chart area displays a result including the diameter and evaporation rate of the pool over time for the evaporation of the pool.
[Equation 1]

Where C _d is the leak coefficient, A _h is the leak area, P ₀ is the pressure inside the vessel, T is the storage temperature, C _p is the static specific heat of the material, C _v is the static specific heat of the material, and Γ is the specific heat ratio C _v / C _p , R is gas constant, K is

, V is the volume of the vessel,

Is the leak rate, T _o is the vessel's initial temperature, T ₁ is the vessel's changed temperature, P _o is the vessel's initial pressure, P ₁ is the vessel's changed pressure, ρ _o is the vessel's initial density, ρ ₁ is the vessel's varied density)
The method of claim 7, wherein
Generating the wind field
Dividing the area around the leak location into a grid of a predetermined size;
Dividing the grid-divided area into an observation area and an unobserved area; And
Predicting atmospheric information of the unobserved region by using the terrain information and weather observation information of the observation region and the terrain information of the non-observation region obtained from the server.
The method of claim 8,
Predicting atmospheric information of the unobserved region
Predicting wind direction and wind speed of the unobserved region based on at least one or more of kinematic effects, slope flow and terrain conditions.
delete delete The method of claim 7, wherein
Modeling the diffusion state of the leaking chemicals
Calculating a concentration of leaking chemical in the wind field.

Images