CN114595553B - Construction method and application of dissolution-volatilization type hazardous chemical substance leakage diffusion model - Google Patents
Construction method and application of dissolution-volatilization type hazardous chemical substance leakage diffusion model Download PDFInfo
- Publication number
- CN114595553B CN114595553B CN202111653914.4A CN202111653914A CN114595553B CN 114595553 B CN114595553 B CN 114595553B CN 202111653914 A CN202111653914 A CN 202111653914A CN 114595553 B CN114595553 B CN 114595553B
- Authority
- CN
- China
- Prior art keywords
- hazardous chemical
- ocean
- atmosphere
- surface layer
- mode
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000383 hazardous chemical Substances 0.000 title claims abstract description 153
- 238000009792 diffusion process Methods 0.000 title claims abstract description 61
- 239000000126 substance Substances 0.000 title claims description 101
- 238000010276 construction Methods 0.000 title abstract description 3
- 239000002344 surface layer Substances 0.000 claims abstract description 59
- 239000010410 layer Substances 0.000 claims abstract description 33
- 238000004088 simulation Methods 0.000 claims abstract description 32
- 230000008878 coupling Effects 0.000 claims abstract description 14
- 238000010168 coupling process Methods 0.000 claims abstract description 14
- 238000005859 coupling reaction Methods 0.000 claims abstract description 14
- 239000013535 sea water Substances 0.000 claims abstract description 11
- 230000002238 attenuated effect Effects 0.000 claims abstract description 8
- 230000008859 change Effects 0.000 claims abstract description 6
- 239000000700 radioactive tracer Substances 0.000 claims description 35
- 238000000034 method Methods 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 20
- 239000003344 environmental pollutant Substances 0.000 claims description 14
- 231100000719 pollutant Toxicity 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 6
- 239000000356 contaminant Substances 0.000 claims description 6
- 238000001704 evaporation Methods 0.000 claims description 5
- 230000008020 evaporation Effects 0.000 claims description 5
- 238000002054 transplantation Methods 0.000 claims 1
- 238000004090 dissolution Methods 0.000 description 5
- 239000000725 suspension Substances 0.000 description 5
- PPBRXRYQALVLMV-UHFFFAOYSA-N Styrene Chemical compound C=CC1=CC=CC=C1 PPBRXRYQALVLMV-UHFFFAOYSA-N 0.000 description 4
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- MVPPADPHJFYWMZ-UHFFFAOYSA-N chlorobenzene Chemical compound ClC1=CC=CC=C1 MVPPADPHJFYWMZ-UHFFFAOYSA-N 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 239000009759 San-Chi Substances 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 150000007522 mineralic acids Chemical class 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/08—Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
- G06Q10/083—Shipping
- G06Q10/0832—Special goods or special handling procedures, e.g. handling of hazardous or fragile goods
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Evolutionary Computation (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- Entrepreneurship & Innovation (AREA)
- Human Resources & Organizations (AREA)
- Geometry (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention belongs to the technical field of warehouse logistics and discloses a construction method and application of a dissolution-volatilization type hazardous chemical leakage diffusion model, after a dissolution-volatilization type hazardous chemical leaks into the sea, the dissolution-volatilization type hazardous chemical is quickly dissolved in seawater and is transported and diffused in the sea in a three-dimensional motion mode, the hazardous chemical on the surface layer of the sea volatilizes into the atmosphere at a certain speed, so that the concentration of the hazardous chemical on the surface layer is attenuated, and the volatilized hazardous chemical enters into the atmosphere bottom layer and is further transported and diffused in the atmosphere in a three-dimensional mode; the model is embedded into an atmosphere-ocean coupling mode system, numerical simulation of transportation and diffusion of the dissolved-volatile hazardous chemical in the ocean and the atmosphere after leakage is achieved, and the simulation result comprises three-dimensional concentration distribution and time change of the hazardous chemical in the ocean and the atmosphere after leakage.
Description
Technical Field
The invention belongs to the technical field of storage logistics, and particularly relates to a method for constructing a soluble and volatile hazardous chemical substance leakage diffusion model, a storage medium for receiving a user input program, a computer program product stored on a computer readable medium, and application of the model in marine transportation of hazardous chemical substances.
Background
In general, hazardous chemicals are transported by sea in liquid form, but after the hazardous chemicals leak into the sea, the hazardous chemicals have different movement forms due to different properties such as density, solubility, volatility and the like. Generally, the dangers with similar properties are mainly in a certain main movement form, and the current leakage diffusion types of the dangers at sea mainly comprise 4 types: (1) Water-insoluble, transporting and diffusing in two dimensions on the water surface, such as oil substances; (2) Easily soluble in water, uniformly distributed in water body, transported and diffused in three-dimensional form, including various inorganic acids and alkalis; (3) Insoluble in water and having a specific gravity greater than or close to that of water, and is deposited on the bottom of water or suspended in a water body for transportation, such as chlorobenzene; (4) High volatility, and can be directly transferred into atmosphere for diffusion, such as liquid ammonia. Then, according to the four types of movement forms, the marine hazardous chemical substances are correspondingly of four types, namely a drift type, a dissolution type, a suspension type and a volatile type (see table 1).
TABLE 1 summary of hazardous chemicals model
The first three models are basically coupled in a module form in an ocean mode (such as a regional ocean mode of a ROMS and a FVOM), wherein a particle (floats) module is adopted in a sea surface drifting type and a suspension transport type, and an inert tracer (tracer) module is mostly adopted in a dissolving diffusion type. And the volatile type is basically coupled in an atmosphere mode in the form of a tracer module, such as an inert tracer (passive tracer) module in a WRF mode and various pollutant modules in a WRF-chem mode.
The previous researches are mainly based on the four models, numerical simulation work of leakage of drifting type, suspension type and dissolving type hazardous chemicals is carried out by independently adopting an ocean mode, or numerical simulation work of leakage of volatile hazardous chemicals is carried out by independently adopting an atmospheric mode, and the simulation condition is single. In the event of a real hazardous chemical leakage, there are types of hazardous chemicals that may have volatility at the surface of the sea. For example, styrene belongs to sea surface drift type hazardous chemicals, but is volatile at the same time. Hazardous chemicals such as liquid ammonia and condensate oil belong to dissolved hazardous chemicals, but are volatile on the sea surface. Therefore, consideration is given to the co-transport and diffusion of hazardous chemicals in the sea and in the atmosphere after they leak into the sea. For example, in the event of leakage of the heroic wheel at the estuary in 4 months in 2001, a large amount of styrene is leaked by the heroic wheel, floats on the sea surface, and is continuously diffused on the sea surface along with the flow of seawater to pollute the estuary sea area, and is continuously volatilized to enter the atmosphere to pollute the atmosphere around the leakage area. In the event of the Sangji wheel in 2018, 1 month, the Sangji wheel sinks into the sea and carries a large amount of condensate oil, and the condensate oil is light oil, can be dissolved in water, can be quickly volatilized on the sea surface after entering the water and is easily naturally weathered. Domestic and foreign researches show that the volatile amount of condensate oil reaches 99% approximately after about 5 hours; within 24 hours, the volatile was almost complete. Therefore, the condensate oil leaked after the Sanchi wheel is sunk is transported and diffused in a three-dimensional mode in seawater, and when the condensate oil vertically moves to the surface layer of the ocean, the condensate oil can be quickly volatilized into the atmosphere, so that the condensate oil can drift and diffuse in the atmosphere.
The dangerous chemical substance leakage process is simply summarized as follows: after leaking into the sea, the hazardous chemical substances are quickly dissolved in the sea water, transported and diffused in the sea, and the hazardous chemical substances reaching the surface layer of the sea volatilize at a certain speed and enter the bottom layer of the atmosphere, so that the hazardous chemical substances are transported and diffused in the atmosphere in a three-dimensional manner. At present, a single hazardous chemical substance leakage model is difficult to depict the hazardous chemical substance leakage transportation diffusion process of the sea-air interface.
Disclosure of Invention
In order to overcome the problems in the related art, the disclosed embodiments of the present invention provide a method for constructing a diffusion model of easily soluble and volatile hazardous chemical substance leakage, a storage medium for receiving a user input program, a computer program product stored on a computer readable medium, and an application in marine transportation of hazardous chemical substances. The technical scheme is as follows:
the method for constructing the easy-to-dissolve and volatile hazardous chemical substance leakage diffusion model comprises the following steps:
after leaking into the sea, the hazardous chemical substances which are easy to dissolve and volatilize are quickly dissolved in the sea water and transported and diffused in the sea in a three-dimensional motion mode, the hazardous chemical substances on the surface layer of the sea volatilize into the atmosphere at a certain speed to cause the concentration attenuation of the hazardous chemical substances on the surface layer, and the volatilized hazardous chemical substances enter the bottom layer of the atmosphere to be transported and diffused in the atmosphere in a three-dimensional mode;
the model is embedded into an atmosphere-ocean coupling mode system, numerical simulation of transportation and diffusion of the dissolved-volatile hazardous chemical in the ocean and the atmosphere after leakage is achieved, and the simulation result comprises three-dimensional concentration distribution and time change of the hazardous chemical in the ocean and the atmosphere after leakage.
In one embodiment, the diffusion equation of the ocean surface for three-dimensional diffusion of dissolution-volatilization type hazardous chemicals is:
wherein t is time; x, y, z spatial position coordinates; u, v and w are flow velocity components (m/s) in the directions of x, y and z; k x Swirl diffusion coefficient (m) in x-direction 2 /s)；K y Vortex diffusion coefficient (m) in y-direction 2 /s)；K z Vortex diffusion coefficient (m) in z direction 2 /s)；v θ Is the molecular diffusion coefficient (m) 2 S); c is the concentration (g/m) of pollutants in the water body 3 )；F C Is a strong source of contaminants (g/m) 3 ·s)，D C Is the dissipative term (g/m) of the contaminant 3 S), reduced V C1 Concentration reduction caused by volatilization of dangerous chemicals on the surface layer of the ocean (vollatilization effect, unit is g/m) 3 ·s)。
In one embodiment, the diffusion equation of the atmospheric bottom layer in the three-dimensional diffusion of the dissolution-volatilization type hazardous chemical is as follows:
increased V C2 The volatile hazardous chemical substances in the sea enter the atmosphere to cause the concentration of the hazardous chemical substances to increase (in g/m) 3 ·s)。
In one embodiment, the three-dimensional diffusion equation for the solution-volatile hazardous chemical is as follows:
in one embodiment, the diffusion equation for the upper atmosphere for the three-dimensional diffusion of the dissolution-volatilization type hazardous chemical is:
the calculation of the volatilization amount of the ocean surface layer hazardous chemical substance refers to a half-life period concentration calculation formula, and the concentration of the ocean surface layer hazardous chemical substance is C at the time t t The volatilization half-life period of the hazardous chemical is T, dt is the calculation time step length, and the concentration of the hazardous chemical on the ocean surface layer is T + dt at the moment of T + dt
The amount of the hazardous chemical substances volatilized in the single grid is as follows:
Wherein h is the thickness of the ocean surface layer;
the volatile danger article in ocean top layer all enters into the atmosphere bottom because atmosphere bottom thickness is generally inconsistent with ocean top layer thickness, therefore atmosphere bottom danger article concentration increases to be:
wherein H is the thickness of the ocean surface layer and H is the thickness of the atmospheric bottom layer.
In one embodiment, after the dissolution-volatilization type hazardous chemical substance model is constructed, the model is embedded into an atmosphere-ocean coupling mode system, and the specific flow is as follows:
step one, in a WRF mode, a passive tracker module option is enabled;
step two, in a WRF mode, adding a pollutant variable defined as tracer _ atm;
step three, in a WRF mode, increasing evaporation variables of the marine surface pollutants, which are defined as ctracers;
step four, in the ROMS mode, a passive tracker module option is started;
step five, in the ROMS mode, adding a pollutant variable defined as tracer _ oc;
step six, in a calculation program of the Tracer _ oc in the ROMS mode, a volatilization attenuation item is added to the Tracer _ oc of the surface layer, and the Tracer _ oc of the ocean surface layer is continuously attenuated along with the increase of the time step;
seventhly, in the MCT coupler, transmitting the amount of tracer _ oc volatilization of the surface layer in the ROMS mode to a ctracer variable;
step eight, reading in a ctracer variable from a coupler when calculating the bottom layer tracer _ atm in a tracer _ atm calculation module in a WRF mode, and enabling the tracer _ atm = tracer _ atm + ctracer to realize calculation of hazardous chemicals volatilized from the surface layer of the ocean entering the atmosphere bottom layer;
step nine, adding an output item in the WRF mode: tracer atm and ctracer.
In one embodiment, after the code is transplanted, the mode coupling calculation process realizes that the dissolution-volatilization type hazardous chemical substances are transported and diffused in the ocean in a three-dimensional motion mode after entering the ocean, the hazardous chemical substances on the surface layer of the ocean volatilize at a certain speed, the concentration of the hazardous chemical substances on the surface layer is attenuated, and the volatilized hazardous chemical substances enter the bottom layer of the atmosphere and are transported and diffused in the atmosphere in a three-dimensional motion mode.
Another object of the present invention is to provide a program storage medium for receiving user input, wherein a stored computer program enables an electronic device to execute a method for constructing a diffusion model of leakage of easily soluble and volatile hazardous chemical substances, comprising:
after leaking into the sea, the hazardous chemical substances which are easy to dissolve and volatilize are quickly dissolved in the sea water and transported and diffused in the sea in a three-dimensional motion mode, the hazardous chemical substances on the surface layer of the sea volatilize into the atmosphere at a certain speed to cause the concentration attenuation of the hazardous chemical substances on the surface layer, and the volatilized hazardous chemical substances enter the bottom layer of the atmosphere to be transported and diffused in the atmosphere in a three-dimensional mode;
the model is embedded into an atmosphere-ocean coupling mode system, numerical simulation of transportation and diffusion of the dissolved-volatile hazardous chemical substances in the ocean and the atmosphere after leakage is achieved, and simulation results comprise three-dimensional concentration distribution and time change of the hazardous chemical substances in the ocean and the atmosphere after leakage.
Another object of the present invention is to provide a computer program product stored on a computer readable medium, which includes a computer readable program for providing a user input interface to implement the method for constructing the easily dissolved and volatilized hazardous chemical leakage diffusion model when the computer program product is executed on an electronic device.
The invention also aims to provide application of the method for constructing the easy-to-dissolve and volatile hazardous chemical substance leakage diffusion model in marine transportation of hazardous chemical substances.
By combining all the technical schemes, the invention has the advantages and positive effects that:
at present, simulation prediction for marine hazardous chemical substance leakage is limited in simulation situations, and simulation prediction can be only carried out for single-type (drift type, suspension type, dissolution type and volatilization type) hazardous chemical substance leakage events, wherein numerical simulation work for drift type, suspension type and dissolution type hazardous chemical substance leakage is carried out by independently adopting an ocean mode, or numerical simulation work for volatile hazardous chemical substances is carried out by independently adopting an atmospheric mode. In a real dangerous chemical leakage event, dangerous chemicals have various attributes, such as liquid ammonia, condensate oil and the like, belong to dissolved dangerous chemicals, are dissolved in seawater, are transported and diffused under the action of ocean current, and are simultaneously volatilized on a sea surface to enter the atmosphere and diffuse in the atmosphere. The single atmospheric hazardous chemical substance leakage simulation and the marine hazardous chemical substance simulation can only simulate the transportation and diffusion of the hazardous chemical substance in the atmosphere or the sea after leakage, and cannot consider the process that the hazardous chemical substance continuously volatilizes to the atmosphere in the transportation and transportation process in the sea, namely, for the atmospheric hazardous chemical substance simulation, a leakage surface source which continuously and dynamically moves exists. Therefore, the co-transport and diffusion of the compound hazardous chemical in the sea and the atmosphere after the compound hazardous chemical leaks into the sea need to be considered. The invention constructs a leakage diffusion model of a dissolving-volatile hazardous chemical, after the hazardous chemical leaks into the sea, the hazardous chemical is quickly dissolved in seawater and is transported and diffused in the sea in a three-dimensional motion mode, and the hazardous chemical on the surface layer of the sea volatilizes at a certain speed and enters the atmosphere to be transported and diffused in the atmosphere. Furthermore, the model is embedded into an atmosphere-ocean coupling mode system, ideal case simulation is carried out, and numerical simulation of transportation and diffusion of the dissolved-volatile hazardous chemical in the ocean and the atmosphere after leakage is achieved. The model mainly aims at the working condition that hazardous chemicals are drifted and diffused after entering the sea and are volatilized into the atmosphere at the same time, has better innovation and practical application value, and has more accurate process portrayal for evaluating the influence of the dissolved-volatilized hazardous chemicals on the marine environment and the atmospheric environment after leaking into the sea.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present disclosure and together with the description, serve to explain the principles of the disclosure.
Fig. 1 is a schematic view of a dissolution-volatilization type hazardous chemical leakage diffusion model provided by an embodiment of the invention.
Wherein, S,
The sub-tables represent ocean vertical stratification and atmospheric vertical stratification, H represents the thickness of the lowest atmosphere layer, and H represents the thickness of the highest ocean layer.
FIG. 2 is a graph showing the time variation of evaporation amount and residual amount in the sea of the dissolution-volatilization type hazardous chemical substances provided by the embodiment of the present invention.
Fig. 3 is a flowchart of a simulated dissolution-volatilization type hazardous chemical model embedded in an atmosphere-ocean coupling mode system according to an embodiment of the present invention.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in detail below. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein, but rather should be construed as broadly as the present invention is capable of modification in various respects, all without departing from the spirit and scope of the present invention.
It will be understood that when an element is referred to as being "secured to" another element, it can be directly on the other element or intervening elements may also be present. When an element is referred to as being "connected" to another element, it can be directly connected to the other element or intervening elements may also be present. As used herein, the terms "vertical," "horizontal," "left," "right," and the like are for purposes of illustration only and are not intended to represent the only embodiments.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
The invention constructs a method for constructing a model for leakage and diffusion of easily soluble and volatile hazardous chemicals, after the easily soluble and volatile hazardous chemicals leak into the sea, the easily soluble and volatile hazardous chemicals are quickly dissolved in seawater and transported and diffused in the sea in a three-dimensional motion mode, and the hazardous chemicals on the surface layer of the sea volatilize into the atmosphere at a certain speed and are further transported and diffused in the atmosphere. The model is embedded into an atmosphere-ocean coupling mode system, numerical simulation of transportation and diffusion of dissolved-volatile hazardous chemicals in the ocean and the atmosphere after leakage is achieved, and the simulation result mainly comprises three-dimensional concentration distribution and time change of the hazardous chemicals in the ocean and the atmosphere after leakage.
As shown in fig. 1, after the dangerous chemical substance leaks from the ocean, the dangerous chemical substance is rapidly dissolved in seawater and is transported and diffused in the ocean in a three-dimensional motion mode, the dangerous chemical substance on the surface layer of the ocean volatilizes at a certain speed, so that the concentration of the dangerous chemical substance on the surface layer is attenuated, and the volatilized dangerous chemical substance enters the atmospheric bottom layer and is further transported and diffused in the atmosphere in a three-dimensional mode. The three-dimensional diffusion equation of the dissolution-volatilization type hazardous chemical is as follows:
1. ocean surface layer:
wherein t is time; x, y, z spatial position coordinates; u, v and w are flow velocity components (m/s) in the directions of x, y and z; k x Swirl diffusion coefficient (m) in x-direction 2 /s)；K y Vortex diffusion coefficient (m) in y-direction 2 /s)；K z The vortex diffusion coefficient (m) in the z direction 2 /s)；v θ Is the molecular diffusion coefficient (m) 2 S); c is the concentration (g/m) of pollutants in the water body 3 )；F C Is a strong source of contaminants (g/m) 3 ·s)，D C Is the dissipative term (g/m) of the contaminant 3 S), reduced V C1 Concentration reduction caused by volatilization of dangerous chemicals on the surface layer of the ocean (vollatilization effect, unit is g/m) 3 ·s)
2. Atmosphere bottom layer:
increased V C2 The volatile hazardous chemical substance in the ocean enters the atmosphere to cause the concentration of the hazardous chemical substance to increase (the unit is g/m) 3 ·s)
3. Ocean lower layer:
4. an upper atmosphere layer:
wherein the calculation of the volatilization amount of the ocean surface layer hazardous chemical substance refers to a half-life period concentration calculation formula, and the concentration of the ocean surface layer hazardous chemical substance is C at the moment t t The half-life of the volatile hazardous chemical is T, the T + dt moment (dt is the calculation time step) is the T, and the concentration of the hazardous chemical on the ocean surface layer is
The amount of the hazardous chemical substances volatilized in the single grid is as follows:
Wherein h is the thickness of the ocean surface layer. The volatile danger article in ocean top layer all enters into the atmosphere bottom because atmosphere bottom thickness is generally inconsistent with ocean top layer thickness, therefore atmosphere bottom danger article concentration increases to be:
Wherein H is the thickness of the ocean surface layer and H is the thickness of the atmospheric bottom layer.
FIG. 2 shows the ideal time-varying graphs of the evaporation amount and the residual amount of the dissolution-volatilization type hazardous chemical, wherein the half-life period is set to 12 hours, the time calculation step is 1 minute, the calculation time is 120 hours, and the initial concentration is set to 100.
After the dissolution-volatilization type hazardous chemical substance model is constructed, the model is embedded into an atmosphere-ocean coupling mode system (a COAWST mode which is composed of a WRF atmosphere mode and an ROMS ocean mode), and the specific flow is as follows:
s101, enabling a passive tracker module option in a WRF mode;
s102, adding a pollutant variable (defined as tracer-atm) in a WRF mode;
s103, increasing the evaporation variable (defined as ctracer) of the marine surface pollutant in a WRF mode;
s104, enabling a passive tracker module option in the ROMS mode;
s105, adding a pollutant variable (defined as tracerOc) in the ROMS mode;
s106, in a program for calculating the tracer _ oc in the ROMS mode, a volatilization attenuation item is added to the tracer _ oc of the surface layer, and the tracer _ oc of the ocean surface layer is continuously attenuated along with the increase of the time step;
s107, in the MCT coupler, transmitting the amount of concentrator _ oc volatilization of the surface layer in the ROMS mode to a ctracer variable;
and S108, in a Tracer _ atm calculation module in the WRF mode, reading in a ctracer variable from a coupler when calculating the lowest Tracer _ atm in each time step, and enabling the Tracer _ atm to be = Tracer _ atm + ctracer. And the calculation of dangerous chemicals volatilized from the surface layer of the ocean entering the bottom layer of the atmosphere is realized.
S109, adding an output item in the WRF mode: tracer _ atm and ctracer.
After the code is transplanted, the mode coupling calculation process realizes that after the dissolution-volatilization type hazardous chemical substances enter the sea, the dissolution is rapidly transported and diffused in the sea in a three-dimensional motion mode, the hazardous chemical substances on the surface layer of the sea volatilize at a certain speed, so that the concentration of the hazardous chemical substances on the surface layer is attenuated, and the volatilized hazardous chemical substances enter the atmosphere bottom layer and are transported and diffused in the atmosphere in a three-dimensional motion mode.
Case simulation:
after the dissolution-volatilization type hazardous chemical substances are implanted into the atmosphere-ocean coupling mode, ideal case simulation is carried out for testing the simulation effect of the model, assuming that the dissolution-volatilization type hazardous chemical substance leakage event occurs in the east black tide region in Taiwan province, the release position of the hazardous chemical substances is selected, the hazardous chemical substances are uniformly released from the surface to the bottom, the initial concentration is 100, the simulation time is 2018, 1 month and 1 day, and the simulation time is 9 days.
According to the concentration distribution of the simulated dissolved-volatile hazardous chemical substances on the ocean mode surface layer and the atmosphere mode bottom layer at different moments after leakage, after the hazardous chemical substances are leaked, the hazardous chemical substances given on the ocean surface layer expand in the northeast direction along with the flow of black tide, the concentration is continuously reduced due to the dissolution and diffusion in the ocean and the volatilization on the surface layer, and the black tide flow axis is drawn by the distribution of the ocean surface layer hazardous chemical substances in the 192 th hour. Simulation shows that the hazardous chemical substances on the surface layer of the ocean continuously volatilize to enter the bottom layer of the atmosphere in the process of drifting along with ocean currents, and then drift and diffuse in the atmosphere, and the hazardous chemical substances drift and diffuse in different directions along with the change of a wind field, for example, in 24 hours of simulation, pollutants volatilized into the atmosphere first drift to the northeast of Taiwan province in China to the west, and then turn to the southwest; in 48 hours, the dangerous chemicals drifting to the south of Taiwan in China drift to the northwest and enter the upper part of Guangdong province; and at the 96 th hour, the wind direction is changed to the northwest, so that the dangerous chemicals on the ocean surface layer mainly drift to the northwest, the concentration of the dangerous chemicals volatilized before is reduced to be lower, at the 192 th hour, the dangerous chemicals on the ocean surface layer are distributed along the main axis of the black tide, and the volatilized dangerous chemicals mainly drift and diffuse to the southwest along the main axis of the black tide due to the northwest wind.
The invention discloses a method for constructing an easily soluble and volatile hazardous chemical substance leakage diffusion model. Further, the model is embedded into an atmosphere-ocean coupling mode system, ideal case simulation is carried out, and numerical simulation of transportation and diffusion of the dissolved-volatile hazardous chemical in the ocean and the atmosphere after leakage is achieved. The model mainly aims at the working condition that hazardous chemicals are drifted and diffused after entering the sea and are volatilized into the atmosphere at the same time, has better innovation and practical application value, and has more accurate process portrayal for evaluating the influence of the dissolved-volatilized hazardous chemicals on the marine environment and the atmospheric environment after leaking into the sea.
Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure disclosed herein. This application is intended to cover any variations, uses, or adaptations of the disclosure following, in general, the principles of the disclosure and including such departures from the present disclosure as come within known or customary practice within the art to which the disclosure pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the disclosure being indicated by the following claims.
It will be understood that the present disclosure is not limited to the precise arrangements described above and shown in the drawings and that various modifications and changes may be made without departing from the scope thereof. The scope of the present disclosure should be limited only by the attached claims.
Claims (7)
1. A method for constructing a soluble and volatile hazardous chemical substance leakage diffusion model is characterized by comprising the following steps:
after leaking into the sea, the hazardous chemical substances which are easy to dissolve and volatilize are quickly dissolved in the sea water and transported and diffused in the sea in a three-dimensional motion mode, the hazardous chemical substances on the surface layer of the sea volatilize into the atmosphere to cause the concentration attenuation of the hazardous chemical substances on the surface layer, and the volatilized hazardous chemical substances enter into the atmosphere bottom layer to be transported and diffused in the atmosphere in a three-dimensional mode;
after the dissolution-volatilization type hazardous chemical substance model is constructed, the model is embedded into an atmosphere-ocean coupling mode system, and the specific flow is as follows:
step one, in a WRF mode, a passive tracker module option is enabled;
step two, in a WRF mode, adding a pollutant variable defined as tracer _ atm;
step three, in a WRF mode, increasing evaporation variables of the marine surface pollutants, which are defined as ctracers;
step four, in the ROMS mode, a passive tracker module option is started;
step five, in the ROMS mode, adding a pollutant variable defined as tracer _ oc;
step six, in a calculation program of the Tracer _ oc in the ROMS mode, a volatilization attenuation item is added to the Tracer _ oc of the surface layer, and the Tracer _ oc of the ocean surface layer is continuously attenuated along with the increase of the time step;
step seven, in the MCT coupler, transferring the volume of tracer _ oc volatilization of the surface layer in the ROMS mode to a ctracer variable;
step eight, reading a ctracer variable from the coupler when calculating the tracer _ atm of the bottommost layer in a tracer _ atm calculation module in a WRF mode and enabling the tracer _ atm = tracer _ atm + ctracer when calculating each time step, so that calculation of hazardous chemical substances volatilized from the ocean surface layer entering the atmosphere bottom layer is realized;
step nine, adding an output item in the WRF mode: tracer _ atm and ctracer;
numerical simulation of transportation and diffusion of the dissolved-volatile hazardous chemical substances in the ocean and the atmosphere after leakage is realized, and simulation results comprise three-dimensional concentration distribution and time change of the hazardous chemical substances in the ocean and the atmosphere after leakage;
the diffusion equation of the ocean surface layer during the three-dimensional diffusion of the dissolution-volatilization type hazardous chemical substances is as follows:
wherein t is time; x, y, z spatial position coordinates; u, v and w are flow velocity components in the directions of x, y and z; k x Swirl diffusion coefficient in x direction; k y Swirl diffusion coefficient in y direction; k is z The swirl diffusion coefficient in the z direction; v. of θ Is the molecular diffusion coefficient; c is the concentration of pollutants in the water body; f C Being a strong source of contaminants, D C Reduced V for dissipative term of contaminant C1 The concentration of hazardous chemicals on the surface layer of the ocean is reduced due to volatilization.
2. The method for constructing the easy-to-dissolve and volatile hazardous chemical leakage diffusion model according to claim 1, wherein the diffusion equation of the atmosphere bottom layer during the three-dimensional diffusion of the dissolution-volatile hazardous chemical is as follows:
increased V C2 The volatile hazardous chemical substance in the ocean enters the atmosphere to cause the concentration of the hazardous chemical substance to increase (the unit is g/m) 3 ·s)。
3. The method for constructing the easy-to-dissolve and volatile hazardous chemical leakage diffusion model according to claim 1, wherein when the dissolution-volatile hazardous chemical is diffused in three dimensions, the diffusion equation of the lower ocean layer is as follows:
4. the method for constructing the easy-to-dissolve and volatile hazardous chemical leakage diffusion model according to claim 1, wherein when the dissolution-volatile hazardous chemical is diffused in three dimensions, the diffusion equation of the upper atmosphere is as follows:
the calculation of the volatilization amount of the ocean surface layer hazardous chemical substance refers to a half-life period concentration calculation formula, and the concentration of the ocean surface layer hazardous chemical substance is C at the time t t The volatilization half-life period of the hazardous chemical is T, dt is the calculation time step length, and the concentration of the hazardous chemical on the ocean surface layer is T + dt at the moment of T + dt
The amount of the hazardous chemical substances volatilized in the single grid is as follows:
Wherein h is the thickness of the ocean surface layer;
the volatile danger article in ocean top layer all enters into the atmosphere bottom because atmosphere bottom thickness is generally inconsistent with ocean top layer thickness, therefore atmosphere bottom danger article concentration increases to be:
wherein H is the thickness of the ocean surface layer and H is the thickness of the atmospheric bottom layer.
5. The method for constructing the easy-to-dissolve and volatile hazardous chemical leakage diffusion model according to claim 1, further comprising after the code transplantation is completed, the mode coupling calculation process realizes that the dissolution-volatile hazardous chemical is transported and diffused in the ocean in a three-dimensional motion form after entering the ocean, while the hazardous chemical on the surface layer of the ocean volatilizes at a certain rate to cause the concentration of the hazardous chemical on the surface layer to be attenuated, and the volatilized hazardous chemical enters the atmospheric bottom layer to be transported and diffused in the atmosphere in a three-dimensional motion form.
6. A program storage medium for receiving user input, wherein the stored computer program causes an electronic device to execute the method for constructing the model of diffusion of leakage of easily soluble and volatile hazardous chemical substances according to any one of claims 1 to 5.
7. The application of the method for constructing the easy-to-dissolve and volatile hazardous chemical leakage diffusion model according to any one of claims 1 to 5 in marine transportation of hazardous chemicals.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111653914.4A CN114595553B (en) | 2021-12-30 | 2021-12-30 | Construction method and application of dissolution-volatilization type hazardous chemical substance leakage diffusion model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111653914.4A CN114595553B (en) | 2021-12-30 | 2021-12-30 | Construction method and application of dissolution-volatilization type hazardous chemical substance leakage diffusion model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114595553A CN114595553A (en) | 2022-06-07 |
| CN114595553B true CN114595553B (en) | 2023-04-07 |
Family
ID=81803871
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111653914.4A Active CN114595553B (en) | 2021-12-30 | 2021-12-30 | Construction method and application of dissolution-volatilization type hazardous chemical substance leakage diffusion model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114595553B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118430670B (en) * | 2024-07-05 | 2024-10-11 | 东海实验室 | Rapid dangerous chemical leakage simulation method and system applied to marine pollution prediction deduction |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR102018281B1 (en) * | 2019-01-29 | 2019-09-05 | 주식회사 세이프티아 | Method and system for analysing behavior of hazardous chemical |
| CN110569514A (en) * | 2018-06-06 | 2019-12-13 | 中国石油化工股份有限公司 | Method for determining leakage pollution range of hazardous chemical substances |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112016226B (en) * | 2020-08-13 | 2021-05-18 | 中国地质大学深圳研究院 | Three-dimensional simulation method and device for offshore oil spill |
| CN113468815A (en) * | 2021-07-12 | 2021-10-01 | 南京航空航天大学 | Liquid hazardous chemical substance volatilization concentration estimation method based on Gaussian diffusion model and deep neural network |
-
2021
- 2021-12-30 CN CN202111653914.4A patent/CN114595553B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110569514A (en) * | 2018-06-06 | 2019-12-13 | 中国石油化工股份有限公司 | Method for determining leakage pollution range of hazardous chemical substances |
| KR102018281B1 (en) * | 2019-01-29 | 2019-09-05 | 주식회사 세이프티아 | Method and system for analysing behavior of hazardous chemical |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114595553A (en) | 2022-06-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114595553B (en) | Construction method and application of dissolution-volatilization type hazardous chemical substance leakage diffusion model | |
| Forthofer et al. | 4.4 simulating diurnally driven slope winds with windninja | |
| Garratt | The stably stratified internal boundary layer for steady and diurnally varying offshore flow | |
| Li et al. | Ecological stability evaluation of tidal flat in coastal estuary: A case study of Liaohe estuary wetland, China | |
| Li et al. | The Impacts of Different PBL Schemes on the Simulation of PM2. 5 during Severe Haze Episodes in the Jing‐Jin‐Ji Region and Its Surroundings in China | |
| Liu et al. | A coupled model of hydrodynamics and water quality for Yuqiao reservoir in Haihe river basin | |
| Wei et al. | Effects of flood on thermal structure of a stratified reservoir | |
| Jin et al. | Application and validation of three-dimensional model in a shallow lake | |
| Koračin et al. | Numerical simulations of the response of the marine atmosphere to ocean forcing | |
| Häfele | A Systems Approach to Energy: Handling large amounts of energy in a way that is safe, clean, cheap, and efficient is a more serious long-range problem than producing an adequate fuel supply | |
| Reese et al. | Local mixing determines spatial structure of diahaline exchange flow in a mesotidal estuary: A study of extreme runoff conditions | |
| Stoner et al. | Application of stereolithographic custom models for studying the impact of biofilms and mineral precipitation on fluid flow | |
| CN113705131B (en) | Water body pollution diffusion simulation prediction method and device based on particle motion | |
| Killworth | A continuously stratified nonlinear ventilated thermocline | |
| Tang et al. | Two-dimensional water environment numerical simulation research based on EFDC in Mudan River, Northeast China | |
| Li | Chinese Science Fiction and Environmental Criticism: From the Anthropocentric to the Cosmocentric | |
| Pyanova et al. | Modeling the propagation of impurities from point sources in the winter atmosphere of the Baikal region | |
| Achiari | Study of Temporal and Spatial changes in the Estuarine Area of the Mayangan River, West Java–Indonesia | |
| Paul et al. | A shallow water numerical method for assessing impacts of hydrodynamics and nutrient transport processes on water quality values of Lake Victoria | |
| Zou et al. | Hydrodynamic and eutrophication modeling for a tidal marsh impacted estuarine system using EFDC | |
| Perianez | An operative Lagrangian model for simulating radioactivity dispersion in the Strait of Gibraltar | |
| Saito | Numerical experiments of the land and sea breeze circulation | |
| Yu | Comments on" Virtual Simulation of Pollution Diffusion in Oil Pollution Accidents". | |
| Ahmad et al. | A solution-adaptive grid generation scheme for atmospheric flow simulations | |
| Lai | Research and development of Poyang Lake hydrodynamic modeling system based on MapWinGIS |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: No. 6 Xianxialing Road, Laoshan District, Qingdao City, Shandong Province Patentee after: FIRST INSTITUTE OF OCEANOGRAPHY, MNR Country or region after: China Patentee after: Qingdao Marine Science and Technology Center Address before: Laoshan District xianxialing road 266061 Shandong city of Qingdao province No. 6 Patentee before: FIRST INSTITUTE OF OCEANOGRAPHY, MNR Country or region before: China Patentee before: QINGDAO NATIONAL LABORATORY FOR MARINE SCIENCE AND TECHNOLOGY DEVELOPMENT CENTER |