CN109783966A - Improved method and device of the Calpuff model to light gas diffusion simulations - Google Patents

Improved method and device of the Calpuff model to light gas diffusion simulations Download PDF

Info

Publication number
CN109783966A
CN109783966A CN201910071708.9A CN201910071708A CN109783966A CN 109783966 A CN109783966 A CN 109783966A CN 201910071708 A CN201910071708 A CN 201910071708A CN 109783966 A CN109783966 A CN 109783966A
Authority
CN
China
Prior art keywords
light gas
parameter
height
calpuff
calpuff model
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.)
Granted
Application number
CN201910071708.9A
Other languages
Chinese (zh)
Other versions
CN109783966B (en
Inventor
许嘉钰
周迅
蔡超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING SANYI SICHUAN TECHNOLOGY Co Ltd
Tsinghua University
Original Assignee
BEIJING SANYI SICHUAN TECHNOLOGY Co Ltd
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BEIJING SANYI SICHUAN TECHNOLOGY Co Ltd, Tsinghua University filed Critical BEIJING SANYI SICHUAN TECHNOLOGY Co Ltd
Priority to CN201910071708.9A priority Critical patent/CN109783966B/en
Publication of CN109783966A publication Critical patent/CN109783966A/en
Application granted granted Critical
Publication of CN109783966B publication Critical patent/CN109783966B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)

Abstract

Present disclose provides a kind of Calpuff models to the improved method of light gas diffusion simulations, obtains accident source information and environmental information, and be collected with environmental information to accident source information;Operation is optimized with sinking speed parameter to the height of release parameter during collection;Parameter after optimization is input in Calpuff model, light gas concentration field result is obtained.Calpuff model of this method after two parameter optimizations of height of release and sinking speed can be suitable for simulating the diffusion mobility of light gas, quickly carry out dynamic analog to the leakage accident of light gas so as to more convenient;Parameter inside the formula of the parameter optimization of proposition is easier to obtain, so as in leakage accident to collect necessary data in the shorter time, to realize the dynamic analog of leakage accident.The disclosure additionally provides Calpuff model to the improvement device of light gas diffusion simulations.

Description

Improved method and device of the Calpuff model to light gas diffusion simulations
Technical field
This disclosure relates to the present invention relates to the dynamic Simulation Techniques field of light gas Release and dispersion, in particular to The improved method and device of a kind of Calpuff model to light gas diffusion simulations.
Background technique
Nowadays there is various dangerous substances in production, storage, transport and use process for petroleum, chemical company.These Dangerous substance wherein just comprising light gas, due to discharge after it is smaller than atmospheric density, under the action of buoyancy, compared to common Gas has apparent lifting;It also can conveying and diffusion to light gas by the turbulent flow reflex from boundary layer Generate certain influence.Detailed process is as follows: when continuous release occurs for light gas, it is formed by plume due to depositing between air In density contrast, the phenomenon that plume has obvious rising under the action of buoyancy;And when the initial temperature for leaking light gas is higher than sky When gas, this plume effect of increasing be will become apparent from.Under the action of cardinal wind, the downward wind direction of the plume of rising is curved Song forms a typical bending plume.It needs to calculate effective source of leaks height of plume at this time;From ground surface and Obviously also the conveying to light gas and diffusion generate certain influence to the auspicious stream reflex in boundary layer.Therefore surface condition It is different with required calculating light gas concentration be it is closely bound up, such as the diffusion in high density settling pond region and wilderness area, finally The regional scope of influence is obviously different.
Calpuff model is the Gaussian cigarette group diffusion for being used to simulate the multilayer of unstable state, several species pollution Mode, it is suitable for simulating the migration, conversion and removing of space-time pollutant all under the meteorological condition of variation.It considers complexity The boundary of the influence of landform, seashore influences, and washes influence, dried wet deposition and simple chemical conversion, Ke Yiji under building Calculate the concentration and settling amount in preset.But there is centainly restricted for current Calpuff model.Specifically, For the lifting after light gas discharge: Calpuff model can consider to pollute when calculating gaseous pollutant height of release Object causes lifting to a certain extent since temperature is less than air themperature, but for light gas, not only its density of temperature with There are density contrasts for air, obvious rising are had under buoyancy, Calpuff model can not simulate this characteristic;For lightweight The diffusion phase of gas: Calpuff model is regarded pollutant air mass density as close with air after the lifting of gas has been calculated Quite, pollutant air mass is spread degree with aerial migration, but for light gas, is more vulnerable to from ground surface and boundary layer Turbulent flow reflex, cause light gas to have a degree of subsidence velocity, Calpuff model after the reflection in boundary layer This characteristic can not be simulated.
Summary of the invention
The technical issues of in order to solve in the prior art, the disclosure is directed to the light gas of different characteristic, to Calpuff mould The parameter of type is adjusted, so that Calpuff can be realized simulation light gas diffusion mobility.I.e. the embodiment of the present disclosure provides A kind of Calpuff model obtains accident source information and environmental information to the improved method and device of light gas diffusion simulations, and Accident source information is collected with environmental information;To during collection height of release parameter and sinking speed parameter carry out it is excellent Change operation;Parameter after optimization is input in Calpuff model, light gas concentration field result is obtained.By height of release With the Calpuff model after two parameter optimizations of sinking speed, it can be suitable for simulating the diffusion mobility of light gas, so as to It is enough more convenient quickly to the leakage accident progress dynamic analog of light gas;Parameter inside the formula of the parameter optimization of proposition It is easier to obtain, so as to, to collect necessary data in the shorter time, leak in leakage accident to realize The dynamic analog of accident.
In a first aspect, the embodiment of the present disclosure provides a kind of Calpuff model to the improvement side of light gas diffusion simulations Method, comprising the following steps: obtain accident source information and environmental information, and the accident source information and the environmental information are carried out It collects;Operation is optimized with sinking speed parameter to the height of release parameter during collection;By the parameter input after optimization Into Calpuff model, light gas concentration field result is obtained.
The accident source information includes: the row of the height of release of light gas, light gas in one of the embodiments, Put the discharge amount of temperature, the mass rate of emission of light gas and light gas.
In one of the embodiments, the environmental information include: ambient wind velocity, Mohs length, friction velocity, season/ Underlying surface.
It collects and the height of release parameter and sinking speed ginseng in preprocessing process for described pair in one of the embodiments, It includes: high for light gas calculating buoyance flux and the lifting for calculating rising for light gas that number, which optimizes operation, Degree.
It collects and the height of release parameter and sinking speed ginseng in preprocessing process for described pair in one of the embodiments, Number optimizes operation further include:
According to the source of leaks parameter of light gas, the rising height risen is calculated by the first preset formula, wherein described First preset formula are as follows:
Wherein, x is lower wind direction distance (m), and u is ambient wind velocity (m/s), and the additional rising height of light gas is with lower wind direction Distance increases and increases, and Δ H is the rising height of the rising as caused by light gas;
The input value of the height of release in Calpuff model is calculated by the second preset formula, wherein described second is default Formula are as follows: Hi=Hs+ΔH;
Wherein, HsFor actual discharge height (m), HiFor the height of release (m) inputted in Calpuff model;
Buoyance flux is calculated according to third preset formula, wherein the third preset formula are as follows:Its In, VsIt is acceleration of gravity (m/s for flue gas emission rate (m/s) at exhaust outlet, g2), TsFor flue gas exit temperature (K), Ta For ambient air temperature (K).
It collects and the height of release parameter and sinking speed ginseng in preprocessing process for described pair in one of the embodiments, Number optimizes operation further include:
According to the environmental information, the sinking speed of gas, the sinking speed V of calculating are calculated based on the 4th preset formulad As the input value of the gas USED TO ESTIMATE DRY DEPOSITION VELOCITY in Calpuff model, the 4th preset formula are as follows:
Wherein, Ra--- turbosphere aerodynamic drag, Rb--- quasi- laminar flow pair layer resistance;RaIt is to be drawn by turbulent motion The resistance being subject to when the pollutant risen is from Atmospheric Transportation to Near Ground, according to surface layer by the mass transport and heat of pollutant It conveys similar hypothesis and obtains the 5th preset formula, wherein the 5th preset formula are as follows:
Wherein, ZrIt is to calculate VdWhen the reference altitude (m) that selects, UaIt is friction velocity (m/s), k is that Von Karman is normal Number, takes 0.4, Z0Be underlying surface roughness it is long, (m) is empirically obtained by local features of terrain, ψcIt is related with mass flux Stability correction function,
Wherein,1/L is the inverse (1/m) of Monin-Obukhov length;
RbPollutant to when surficial deposit by resistance suffered by laminar film near the ground, and by connecting with earth's surface The transmission of class laminar film is associated;
Wherein, UaIt is friction velocity (m/s), k is Von Karman constant, takes 0.4, PrFor Prandtl number, value 0.72, ScFor Schmidt number, Sc0/Di, μ0For air movement viscosity, value 0.15cm2s-1, DiFor the 1st kind of dirt relevant to component Contaminate the molecule diffusance of object.
In one of the embodiments, further include: by judging that multiple underlying surface types of multiple seasons obtain receptor surface Resistance.
Second aspect, the embodiment of the present disclosure provide a kind of computer readable storage medium, are stored thereon with computer journey The step of sequence, which realizes above-mentioned method when being executed by processor.
The third aspect, the embodiment of the present disclosure provide a kind of computer equipment, including memory, processor and are stored in On reservoir and the computer program that can run on a processor, the processor realize above-mentioned method when executing described program Step.
Fourth aspect, the embodiment of the present disclosure provide a kind of Calpuff model and fill to the improvement of light gas diffusion simulations It sets, described device includes: acquisition and collection module, is believed for obtaining accident source information and environmental information, and to the accident source Breath is collected with the environmental information;Optimization module, for the height of release parameter and sinking speed ginseng during collection Number optimizes operation;As a result module is obtained, for the parameter after optimization to be input in Calpuff model, obtains lightweight gas Bulk concentration field result.
A kind of Calpuff model provided by the invention obtains thing to the improved method and device of light gas diffusion simulations Therefore source information and environmental information, and accident source information is collected with environmental information;To the height of release ginseng during collection It is several to optimize operation with sinking speed parameter;Parameter after optimization is input in Calpuff model, it is dense to obtain light gas Spend field result.Calpuff model of this method after two parameter optimizations of height of release and sinking speed, can be suitable for mould The diffusion mobility of quasi- light gas quickly carries out dynamic analog to the leakage accident of light gas so as to more convenient;It mentions Parameter inside the formula of parameter optimization out is easier to obtain, so as in leakage accident in the shorter time Necessary data were collected, to realize the dynamic analog of leakage accident.
Detailed description of the invention
In order to illustrate more clearly of the technical solution of the embodiment of the present disclosure, below to needed in embodiment description Attached drawing is briefly described:
Fig. 1 is improved method of one of the one embodiment of the invention Calpuff model to light gas diffusion simulations Step flow diagram;
Fig. 2 is improved method of one of another embodiment of the present invention Calpuff model to light gas diffusion simulations Step flow diagram;
Fig. 3 is improved method of one of the one embodiment of the invention Calpuff model to light gas diffusion simulations In lower wind direction center line on it is maximum when equal concentration change schematic diagram;
Fig. 4 is improved method of one of the one embodiment of the invention Calpuff model to light gas diffusion simulations In leeward to the ground when equal concentration field schematic diagram;
Fig. 5 is improvement device of one of the one embodiment of the invention Calpuff model to light gas diffusion simulations Structural schematic diagram.
Specific embodiment
The application is further discussed in detail with reference to the accompanying drawings and examples.
In following introductions, term " first ", " second " only for descriptive purposes, and should not be understood as instruction or dark Show relative importance.Following introductions provide multiple embodiments of the disclosure, can replace or merge between different embodiments Combination, therefore the application is it is also contemplated that all possible combinations comprising documented identical and/or different embodiments.Thus, such as Fruit one embodiment include feature A, B, C, another embodiment include feature B, D, then the application also should be regarded as include containing A, the every other possible combined embodiment of one or more of B, C, D, although the embodiment may be in the following contents In have specific literature record.
In order to make the objectives, technical solutions, and advantages of the present invention clearer, by the following examples, it and combines attached Figure, to a kind of Calpuff model of the present invention to the specific embodiment of the improved methods of light gas diffusion simulations and device into Row is further described.It should be appreciated that described herein, specific examples are only used to explain the present invention, is not used to limit The present invention.
As shown in Figure 1, being one of one embodiment Calpuff model to the improved method of light gas diffusion simulations Flow diagram, specifically includes the following steps:
Step 102, accident source information and environmental information are obtained, and accident source information is collected with environmental information.Tool Body, accident source information includes: the mass rate of emission of the height of release of light gas, the exhaust temperature of light gas, light gas And the discharge amount of light gas.The diversity and selectivity of the accident source information acquisition for light gas are improved as a result,. In addition, environmental information includes: ambient wind velocity, Mohs length, friction velocity, season/underlying surface.Environmental information is improved as a result, Multi-selection.
Step 104, operation is optimized with sinking speed parameter to the height of release parameter during collection.It needs to illustrate , include calculating buoyance flux and calculating the rising height risen in height of release parameter optimization;Sinking speed parameter is excellent Change includes calculating turbulent motion resistance, calculate resistance near the ground and calculating receptor surface resistance.Specifically, to collecting and pre- place It includes: to calculate buoyance flux for light gas that height of release parameter and sinking speed parameter during reason, which optimize operation, And the rising height risen is calculated for light gas.
Further, operation is optimized with sinking speed parameter to collection and the height of release parameter in preprocessing process Further include:
According to the source of leaks parameter of light gas, the rising height risen is calculated by the first preset formula, wherein first Preset formula are as follows:
Wherein, x is lower wind direction distance (m), and u is ambient wind velocity (m/s), and the additional rising height of light gas is with lower wind direction Distance increases and increases, and Δ H is the rising height of the rising as caused by light gas;
The input value of the height of release in Calpuff model is calculated by the second preset formula, wherein the second preset formula Are as follows: Hi=Hs+ΔH;
Wherein, HsFor actual discharge height (m), HiFor the height of release (m) inputted in Calpuff model;
Buoyance flux is calculated according to third preset formula, wherein the third preset formula are as follows:Its In, VsIt is acceleration of gravity (m/s for flue gas emission rate (m/s) at exhaust outlet, g2), TsFor flue gas exit temperature (K), Ta For ambient air temperature (K).
Further, behaviour is optimized with sinking speed parameter to collection and the height of release parameter in preprocessing process Make further include:
According to environmental information, the sinking speed of gas, the sinking speed V of calculating are calculated based on the 4th preset formuladAs The input value of gas USED TO ESTIMATE DRY DEPOSITION VELOCITY in Calpuff model, the 4th preset formula are as follows:
Wherein, Ra--- turbosphere aerodynamic drag, Rb--- quasi- laminar flow pair layer resistance;RaIt is to be drawn by turbulent motion The resistance being subject to when the pollutant risen is from Atmospheric Transportation to Near Ground, according to surface layer by the mass transport and heat of pollutant It conveys similar hypothesis and obtains the 5th preset formula, wherein the 5th preset formula are as follows:
Wherein, ZrIt is to calculate VdWhen the reference altitude (m) that selects, UaIt is friction velocity (m/s), k is that Von Karman is normal Number, takes 0.4, Z0Be underlying surface roughness it is long, (m) is empirically obtained by local features of terrain, ψcIt is related with mass flux Stability correction function,
Wherein,1/L is the inverse (1/m) of Monin-Obukhov length;
RbPollutant to when surficial deposit by resistance suffered by laminar film near the ground, and by connecting with earth's surface The transmission of class laminar film is associated;
Wherein, UaIt is friction velocity (m/s), k is Von Karman constant, takes 0.4, PrFor Prandtl number, value 0.72, ScFor Schmidt number, Sc0/Di, μ0For air movement viscosity, value 0.15cm2s-1, DiFor the 1st kind of dirt relevant to component Contaminate the molecule diffusance of object.
In addition it is also necessary to illustrate, improvement of the Calpuff model to light gas diffusion simulations involved in the disclosure Method further include: by judging that multiple underlying surface types of multiple seasons obtain receptor surface resistance.
Step 106, the parameter after optimization is input in Calpuff model, obtains light gas concentration field result.
A kind of Calpuff model provided by the invention obtains accident source letter to the improved method of light gas diffusion simulations Breath and environmental information, and accident source information is collected with environmental information;To the height of release parameter during collection and sink Drop speed parameter optimizes operation;Parameter after optimization is input in Calpuff model, light gas concentration field knot is obtained Fruit.Calpuff model of this method after two parameter optimizations of height of release and sinking speed can be suitable for simulation lightweight The diffusion mobility of gas quickly carries out dynamic analog to the leakage accident of light gas so as to more convenient;The ginseng of proposition Number optimization formula inside parameter be easier to obtain, so as in leakage accident to be collected in the shorter time Necessary data, to realize the dynamic analog of leakage accident.
Light gas is expanded with a kind of Calpuff model proposed with the disclosure in order to more precisely, neatly understand The improved method of simulation is dissipated, following example is carried out.It should be noted that the range that the disclosure is protected is not limited to following example.
It should be noted that Schmidt the and Prandtl number of usually gas is related to molecular weight, with specific reference to table 1, i.e., Table 1 is for Schmidt the and Prandtl number of several gases and relative to water diffusion coefficient ratio.
In addition it is also necessary to explanation, RcIt is receptor surface resistance, includes other a variety of resistances, usually not by judgement It is obtained with season different underlying surface type.Meaning as shown in table 2.The skin resistance suffered when depositing of the various underlying surfaces of table 2.
Specifically, i.e. Fig. 2 is one of another embodiment of the present invention Calpuff model to lightweight in conjunction with shown in Fig. 2-4 The step flow diagram of the improved method of gas diffusion simulation;Fig. 3 is one of one embodiment of the invention Calpuff mould Type is to equal concentration changes schematic diagram when maximum on the lower wind direction center line in the improved method of light gas diffusion simulations;Fig. 4 is One of one embodiment of the invention Calpuff model to the leeward in the improved method of light gas diffusion simulations to the ground When equal concentration field schematic diagram.For the test of ammonia toxic in document, toxic ammonia is on daytime because contingency generation is continuously let out Leakage, source strength 20kg/s, leakage initial rate are 1.5m/s, and practical leakage height is 3m, define ammonia and the atmospheric environment temperature difference The factor is 0.1, and wind speed is 4m/s, atmospheric stability C, surface roughness 0.5m when accident occurs.
Firstly, being calculated about height of release, the parameter arrived according to on-site collection, calculating buoyance flux F firstb:
Then it calculates and discharges height as caused by light gas Spend incrementss Δ H.
Δ H=1.6 × 1.471/3×4-1×(3×10)2/3=4.4;The then height of release H Jing Guo parameter optimizationi=Hs+ Δ H=3+4.4=7.4.
Further, it is calculated about sinking speed, the parameter arrived according to on-site collection, calculates caused by turbulent motion first Pollutant from Atmospheric Transportation to Near Ground when the resistance R that is subject toa:
Then calculate pollutant to when surficial deposit by resistance R suffered by laminar film near the groundb.I.e.Then receptor surface resistance R is judged according to live season and underlying surface typec;Rc=0; Finally calculate light gas USED TO ESTIMATE DRY DEPOSITION VELOCITY Vd
In conclusion Calpuff calculated result and with measured result compare, by the height of release of above-mentioned calculating and sedimentation speed Degree is inputted as Calpuff model parameter, finally obtains the concentration field of this leakage accident liquid chlorine.According to national health department's phase Regulation is closed, when ammonia concentration reaches 30mg/m3When, there is minor response in human body;Reach 350mg/m3When, it is anti-that moderate occurs in human body It answers;Reach 550mg/m3When, there is kickback in human body;Reach 1500mg/m3When, it may cause death after contacting a period of time; When concentration is greater than 3500mg/m3When, then it is dead immediately after may sucking.Simulation is calculated according to this five toxicity thresholds Ammonia concentration influence area is divided, away from source of leaks from closely to far being successively divided into: lethal zone, dead area, strong immediately Reaction zone, moderate reaction zone and minor response area.As shown in table 3, i.e. the simulation of table 3 calculates concentration and document experimental concentration pair Than.
As can be seen from the above results, in addition to minor response area, the length of the different threshold ranges of Calpuff estimation before improving Degree is all larger than document test value, and the length is coincide with document test value substantially after improving, and is illustrated into excessively improved Calpuff Model can be good at simulating light gas.
In summary, the invention discloses a kind of improved method based on Calpuff modeling light gas, propositions pair The input parameter height of release and USED TO ESTIMATE DRY DEPOSITION VELOCITY of Calpuff model carry out parameter optimization, finally utilize Calpuff model dynamic Simulate the diffusion mobility of light gas.Test result shows: the algorithm can simulate light gas diffusion mobility characteristic well. The disclosure is adjusted the parameter of Calpuff model for the light gas of different characteristic by a kind of algorithm, so that Calpuff can be realized simulation light gas diffusion mobility.To two parameters of height of release in Calpuff model and sinking speed It optimizes, to realize simulation of the Calpuff model for light gas.
Based on the same inventive concept, a kind of Calpuff model is additionally provided to the improvement device of light gas diffusion simulations. The principle solved the problems, such as due to this device is with a kind of aforementioned Calpuff model to the improved method phase of light gas diffusion simulations Seemingly, therefore, the implementation of the device can realize that overlaps will not be repeated according to the specific steps of preceding method.
As shown in figure 5, being one of one embodiment Calpuff model to the improvement device of light gas diffusion simulations Structural schematic diagram.The Calpuff model includes: acquisition and collection module to the improvement device 10 of light gas diffusion simulations 200, optimization module 400 and result obtain module 600.
Wherein, it obtains with collection module 200 for obtaining accident source information and environmental information, and to accident source information and ring Border information is collected;Optimization module 400 is used to carry out the height of release parameter during collection with sinking speed parameter excellent Change operation;As a result it obtains module 600 to be used to for the parameter after optimization being input in Calpuff model, obtains light gas concentration Field result.
A kind of Calpuff model provided by the invention passes through acquisition to the improvement device of light gas diffusion simulations first Accident source information and environmental information are obtained with collection module, and accident source information is collected with environmental information;Again by excellent Change module and operation is optimized with sinking speed parameter to the height of release parameter during collection;Mould is obtained eventually by result Parameter after optimization is input in Calpuff model by block, obtains light gas concentration field result.The device passes through height of release With the Calpuff model after two parameter optimizations of sinking speed, it can be suitable for simulating the diffusion mobility of light gas, so as to It is enough more convenient quickly to the leakage accident progress dynamic analog of light gas;Parameter inside the formula of the parameter optimization of proposition It is easier to obtain, so as to, to collect necessary data in the shorter time, leak in leakage accident to realize The dynamic analog of accident.
The embodiment of the invention also provides a kind of computer readable storage medium, stored on the computer readable storage medium There is computer program, which is executed by processor in Fig. 1 and Fig. 2.The embodiment of the invention also provides a kind of comprising instruction Computer program product.When the computer program product is run on computers, so that computer executes above-mentioned Fig. 1 and Fig. 2 Method.
Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with Relevant hardware is instructed to complete by computer program, the program can be stored in a computer-readable storage medium In, the program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein, the storage medium can be magnetic Dish, CD, read-only memory (Read-Only Memory, ROM) or random access memory (Random Access Memory, RAM) etc..The embodiments described above only express several embodiments of the present invention, and description is more specific and detailed Carefully, but it cannot be understood as limitations on the scope of the patent of the present invention.It should be pointed out that for the common skill of this field For art personnel, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to this hair Bright protection scope.Therefore, the scope of protection of the patent of the invention shall be subject to the appended claims.For the mesh of example and description Have been presented for above description.This description is not intended to for embodiment of the disclosure to be restricted to form disclosed herein.Although with On discussed multiple exemplary aspects and embodiment, but those skilled in the art will appreciate that its certain modification, modify, change Change, addition and sub-portfolio.

Claims (10)

1. a kind of Calpuff model is to the improved method of light gas diffusion simulations, which comprises the following steps:
Acquisition accident source information and environmental information, and the accident source information is collected with the environmental information;
Operation is optimized with sinking speed parameter to the height of release parameter during collection;
Parameter after optimization is input in Calpuff model, light gas concentration field result is obtained.
2. Calpuff model according to claim 1 is to the improved method of light gas diffusion simulations, which is characterized in that The accident source information include: the height of release of light gas, the exhaust temperature of light gas, light gas mass rate of emission with And the discharge amount of light gas.
3. Calpuff model according to claim 1 is to the improved method of light gas diffusion simulations, which is characterized in that The environmental information includes: ambient wind velocity, Mohs length, friction velocity, season/underlying surface.
4. Calpuff model according to claim 1 is to the improved method of light gas diffusion simulations, which is characterized in that It includes: for lightweight that described pair of collection, which optimizes operation with sinking speed parameter with the height of release parameter in preprocessing process, Gas calculates buoyance flux and calculates the rising height risen for light gas.
5. Calpuff model according to claim 1 is to the improved method of light gas diffusion simulations, which is characterized in that It collects for described pair and optimizes operation with sinking speed parameter with the height of release parameter in preprocessing process further include:
According to the source of leaks parameter of light gas, the rising height risen is calculated by the first preset formula, wherein described first Preset formula are as follows:
Wherein, x is lower wind direction distance (m), and u is ambient wind velocity (m/s), and the additional rising height of light gas is with lower wind direction distance Increase and increase, Δ H is the rising height of the rising as caused by light gas;
The input value of the height of release in Calpuff model is calculated by the second preset formula, wherein second preset formula Are as follows: Hi=Hs+ΔH;
Wherein, HsFor actual discharge height (m), HiFor the height of release (m) inputted in Calpuff model;
Buoyance flux is calculated according to third preset formula, wherein the third preset formula are as follows:Wherein, Vs It is acceleration of gravity (m/s for flue gas emission rate (m/s) at exhaust outlet, g2), TsFor flue gas exit temperature (K), TaFor environment Atmospheric temperature (K).
6. Calpuff model according to claim 1 is to the improved method of light gas diffusion simulations, which is characterized in that It collects for described pair and optimizes operation with sinking speed parameter with the height of release parameter in preprocessing process further include:
According to the environmental information, the sinking speed of gas, the sinking speed V of calculating are calculated based on the 4th preset formuladAs The input value of gas USED TO ESTIMATE DRY DEPOSITION VELOCITY in Calpuff model, the 4th preset formula are as follows:
Wherein, Ra--- turbosphere aerodynamic drag, Rb--- quasi- laminar flow pair layer resistance;RaIt is as caused by turbulent motion The resistance being subject to when pollutant is from Atmospheric Transportation to Near Ground is conveyed according to surface layer by the mass transport and heat of pollutant Similar hypothesis obtains the 5th preset formula, wherein the 5th preset formula are as follows:
Wherein, ZrIt is to calculate VdWhen the reference altitude (m) that selects, UaIt is friction velocity (m/s), k is Von Karman constant, is taken 0.4, Z0Be underlying surface roughness it is long, (m) is empirically obtained by local features of terrain, ψcFor stabilization related with mass flux Correction function is spent,
Wherein,1/L is the inverse (1/m) of Monin-Obukhov length;
RbPollutant to when surficial deposit by resistance suffered by laminar film near the ground, and the such sheets by connecting with earth's surface The transmission of fluid layer is associated;
Wherein, UaIt is friction velocity (m/s), k is Von Karman constant, takes 0.4, PrFor Prandtl number, value 0.72, ScFor Schmidt number, Sc0/Di, μ0For air movement viscosity, value 0.15cm2s-1, DiFor the 1st kind of pollutant relevant to component Molecule diffusance.
7. Calpuff model according to claim 6 is to the improved method of light gas diffusion simulations, which is characterized in that Further include: by judging that multiple underlying surface types of multiple seasons obtain receptor surface resistance.
8. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is held by processor The step of any one of claim 1-7 the method is realized when row.
9. a kind of computer equipment including memory, processor and stores the meter that can be run on a memory and on a processor Calculation machine program, which is characterized in that the processor realizes side described in any one of described claim 1-7 when executing described program The step of method.
10. a kind of Calpuff model is to the improvement device of light gas diffusion simulations accuracy, which is characterized in that described device Include:
Acquisition and collection module, for obtaining accident source information and environmental information, and to the accident source information and the environment Information is collected;
Optimization module, for optimizing operation with sinking speed parameter to the height of release parameter during collection;
As a result module is obtained, for the parameter after optimization to be input in Calpuff model, obtains light gas concentration field knot Fruit.
CN201910071708.9A 2019-01-25 2019-01-25 Method and device for improving simulation of light gas diffusion by Calpuff model Active CN109783966B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910071708.9A CN109783966B (en) 2019-01-25 2019-01-25 Method and device for improving simulation of light gas diffusion by Calpuff model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910071708.9A CN109783966B (en) 2019-01-25 2019-01-25 Method and device for improving simulation of light gas diffusion by Calpuff model

Publications (2)

Publication Number Publication Date
CN109783966A true CN109783966A (en) 2019-05-21
CN109783966B CN109783966B (en) 2020-12-18

Family

ID=66502500

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910071708.9A Active CN109783966B (en) 2019-01-25 2019-01-25 Method and device for improving simulation of light gas diffusion by Calpuff model

Country Status (1)

Country Link
CN (1) CN109783966B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111832933A (en) * 2020-07-10 2020-10-27 陕西金合信息科技股份有限公司 Emergency resource scheduling method and device
CN112330812A (en) * 2020-11-05 2021-02-05 桂林理工大学 Gas diffusion visualization method and system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130208262A1 (en) * 2010-10-06 2013-08-15 Tea Sistemi S.P.A. Method for monitoring fugitive gas emissions from the soil, via vertical concentration measurements
CN106021817A (en) * 2016-06-22 2016-10-12 西南石油大学 Rapid simulation method of leakage of high-sulfur natural gas gathering and transportation device for marine gas field
CN106407714A (en) * 2016-10-14 2017-02-15 珠海富鸿科技有限公司 Air pollution assessment method and device based on CALPUFF system
CN107229834A (en) * 2017-06-27 2017-10-03 国网江苏省电力公司电力科学研究院 A kind of complicated landform emergency response air pollution DIFFUSION PREDICTION method
CN107633148A (en) * 2017-10-16 2018-01-26 核工业北京地质研究院 A kind of method for numerical simulation for nuclear facilities exhaust gas diffusion
CN109117549A (en) * 2018-08-08 2019-01-01 浙江航天恒嘉数据科技有限公司 A kind of pollutant DIFFUSION PREDICTION method and system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130208262A1 (en) * 2010-10-06 2013-08-15 Tea Sistemi S.P.A. Method for monitoring fugitive gas emissions from the soil, via vertical concentration measurements
CN106021817A (en) * 2016-06-22 2016-10-12 西南石油大学 Rapid simulation method of leakage of high-sulfur natural gas gathering and transportation device for marine gas field
CN106407714A (en) * 2016-10-14 2017-02-15 珠海富鸿科技有限公司 Air pollution assessment method and device based on CALPUFF system
CN107229834A (en) * 2017-06-27 2017-10-03 国网江苏省电力公司电力科学研究院 A kind of complicated landform emergency response air pollution DIFFUSION PREDICTION method
CN107633148A (en) * 2017-10-16 2018-01-26 核工业北京地质研究院 A kind of method for numerical simulation for nuclear facilities exhaust gas diffusion
CN109117549A (en) * 2018-08-08 2019-01-01 浙江航天恒嘉数据科技有限公司 A kind of pollutant DIFFUSION PREDICTION method and system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
伯鑫 等: "CALPUFF模式的标准化应用技术研究", 《环境科学与技术》 *
张玉珍: "基于CALPUFF模型的核电站事故源项大气扩散研究", 《中国优秀硕士学位论文全文数据库》 *
方露瑶: "CALPUFF模型的应用与研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111832933A (en) * 2020-07-10 2020-10-27 陕西金合信息科技股份有限公司 Emergency resource scheduling method and device
CN112330812A (en) * 2020-11-05 2021-02-05 桂林理工大学 Gas diffusion visualization method and system
CN112330812B (en) * 2020-11-05 2023-01-10 桂林理工大学 Gas diffusion visualization method and system

Also Published As

Publication number Publication date
CN109783966B (en) 2020-12-18

Similar Documents

Publication Publication Date Title
O'Donnell Surface fronts in estuaries: a review
Sullivan et al. Large-eddy simulation of marine atmospheric boundary layers above a spectrum of moving waves
Lazure et al. 3D modelling of seasonal evolution of Loire and Gironde plumes on Biscay Bay continental shelf
Joffre Momentum and heat transfers in the surface layer over a frozen sea
Moraes et al. Surface layer turbulence parameters over a complex terrain
Fan et al. A new convective velocity scale for studying diurnal urban heat island circulation
CN107341341A (en) A kind of river mouth point source sudden water pollution event source tracing method
CN109783966A (en) Improved method and device of the Calpuff model to light gas diffusion simulations
Grachev et al. Turbulent fluxes and transfer of trace gases from ship‐based measurements during TexAQS 2006
CN110020448A (en) Boundary layer meteorology parameter improvement type Mesoscale photochemical pollution simulation and forecast algorithm
Kurose et al. Direct numerical simulation of turbulent heat transfer across a sheared wind-driven gas–liquid interface
Tedeschi et al. Development of a 2D marine aerosol transport model: Application to the influence of thermal stability in the marine atmospheric boundary layer
Udina et al. Pollutant vertical mixing in the nocturnal boundary layer enhanced by density currents and low-level jets: two representative case studies
Ozoe et al. Characteristics of air pollution in the presence of land and sea breeze—A numerical simulation
Gelderloos et al. Coastal trapped waves and other subinertial variability along the southeast Greenland coast in a realistic numerical simulation
CN109840374A (en) The improved method and device that Calpuff model simulates heavy gas diffusion
Tedeschi et al. Study of vertical transport of marine aerosol using an unsteady 2D model
Liu et al. A laboratory simulation of plume dispersion in stratified atmospheres over complex terrain
Turner et al. Application of wake shielding effects with a finite element net model in determining hydrodynamic loading on aquaculture net pens
CN108197422B (en) A kind of water age measuring method in semiclosed waters
Bristow et al. Topographic perturbation of turbulent boundary layers by low‐angle, early‐stage aeolian dunes
Barthelmie Evaluating the impact of wind induced roughness change and tidal range on extrapolation of offshore vertical wind speed profiles
Grønås et al. Numerical simulations of local winds over steep orography in the storm over north Norway on October 12, 1996
Zufall et al. Dry deposition of particles to wave surfaces: I. Mathematical modeling
Monti et al. Waves and turbulence in katabatic winds

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