CN115638938A - Leakage protection distance measuring system and method based on positioning monitoring - Google Patents
Leakage protection distance measuring system and method based on positioning monitoring Download PDFInfo
- Publication number
- CN115638938A CN115638938A CN202211331979.1A CN202211331979A CN115638938A CN 115638938 A CN115638938 A CN 115638938A CN 202211331979 A CN202211331979 A CN 202211331979A CN 115638938 A CN115638938 A CN 115638938A
- Authority
- CN
- China
- Prior art keywords
- diffusion
- leakage
- ground
- parameter
- reference value
- 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.)
- Pending
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
The invention discloses a leakage protection distance measuring system and method based on positioning monitoring, which comprises an input module, a diffusion parameter reference value determining module, an atmospheric diffusion parameter determining module related to the effective roughness length of the ground, an atmospheric diffusion parameter determining module related to leakage time, a gas leakage diffusion model module, a leakage diffusion hazard radius determining module and an output module, wherein the diffusion parameter reference value is obtained by adopting a Tennel method according to meteorological conditions and geographic conditions, and a gas leakage diffusion model is established according to wind speed, wind direction, leakage sources and atmospheric diffusion parameters; and determining the leakage diffusion hazard radius according to the gas leakage diffusion model. According to the invention, the concentration of the leaked gas and the leakage diffusion damage radius can be obtained according to the meteorological condition, the geographic condition and the diffusion time, so that the problem of difficulty in acquiring field monitoring data is avoided.
Description
Technical Field
The invention relates to a leakage protection distance measuring system and method based on positioning monitoring, and belongs to the technical field of leakage monitoring and prediction.
Background
When a dangerous gas leakage accident occurs, the determination of the leakage position and the leakage source intensity is one of the bases and bases for establishing an emergency scheme. In general, emergency personnel cannot directly determine the leakage position and measure the leakage source strength due to the limiting factors such as safety, terrain and time. At the moment, the accident site and the periphery are monitored through the inherent monitoring device or the portable movable monitoring equipment, the position and the intensity of the accident source, the influence range, the degree and the like are determined by combining the acquired monitoring data with a back calculation method, reasonable guidance can be provided for formulating disaster escape routes and rescue measures, and the method has certain guiding significance on occasions such as gas leakage emergency monitoring, air pollution source tracing back check and the like.
(1) Method of probability statistics
And (3) carrying out statistical analysis by using different algorithms according to the numerical distribution rule of atmospheric diffusion, thereby reflecting the position or source intensity of the leakage source. The theoretical basis is bayes' theorem. And obtaining posterior probability distribution of the parameters by using the prior information of the parameters to be solved and actual observation data, and sampling the posterior probability distribution to obtain an estimated value of the parameters to be solved.
(2) Optimization theory method
And constructing a proper objective function and calculating an optimal solution of the objective function. In the optimized leakage source parameter back calculation, the accident is simulated through a diffusion model, and the simulation result of the model is compared with the measurement result to establish an objective function. The objective function is optimized through various optimization algorithms, and the source parameters are continuously adjusted, so that the calculation result and the measurement result can be better matched (the matching degree is quantified by the objective function).
The optimization model method is used for researching the leakage source inverse calculation problem and is summarized as the matching degree of the calculated concentration obtained by the diffusion mode (Gaussian model) and the actual observed concentration, namely the optimization solution of the objective function. The optimization algorithm involved is as follows: genetic Algorithm (GA), simulated Annealing (SA), particle swarm algorithm (PSO), pattern search, and the like.
(3) Machine learning algorithm
For example, an artificial neural network is established, adaptive dimensions and weights are set, and the correlation between the source strength and the concentration is obtained through a training network (a large amount of reliable data is needed to train the network) according to a large amount of actual measurement data (the measured concentration and the corresponding source strength value). And after the training is finished, calculating the source intensity reversely. The model and meteorological data such as diffusion parameters, stability and the like are not required to be known in advance.
The input vector is monitoring data of the sensor array, and the output vector is a description variable of the leakage condition. Each leakage source is provided with corresponding output variables which respectively represent leakage position, leakage area, leakage time and leakage quantity (leakage rate when continuous leakage occurs). In the process of research, a group of output vectors are used for representing the leakage probability of a plurality of storage tanks, and in the case of single leakage, if only one output quantity close to 1 appears in the output vectors, and the other output quantities are all close to 0, the position of the leakage tank can be judged.
Disclosure of Invention
The purpose of the invention is as follows: in order to solve the problem of difficulty in acquiring field monitoring data, the invention provides a leakage protection distance measuring system and method based on positioning monitoring.
The technical scheme is as follows: in order to achieve the purpose, the invention adopts the technical scheme that:
a leakage protection distance measuring method based on positioning monitoring comprises the following steps:
step 1, obtaining a diffusion parameter reference value by adopting a Ternal method according to meteorological conditions and geographic conditions.
And 2, determining atmospheric diffusion parameters related to the effective roughness length of the ground according to the effective roughness length of the ground and the diffusion parameter reference value.
And 3, determining the atmospheric diffusion parameter related to the leakage time according to the leakage time and the diffusion parameter reference value.
And 4, establishing a gas leakage diffusion model according to the wind speed, the wind direction, the leakage source, the atmospheric diffusion parameter related to the effective roughness length of the ground and the atmospheric diffusion parameter related to the leakage time.
And 5, determining the leakage diffusion damage radius according to the gas leakage diffusion model.
Preferably: in the step 1, a method for obtaining a diffusion parameter reference value by adopting a Ternal method according to meteorological conditions and geographic conditions comprises the following steps:
wherein the content of the first and second substances,indicating a reference value of a lateral wind direction dispersion parameter、Representing a lateral diffusion parameter power function coefficient;、representing the coefficients of a power function of a vertical diffusion parameter,the distance in the down-wind direction is indicated,indicating a vertical diffusion parameter reference value.
Preferably: in step 2, determining atmospheric diffusion parameters related to the effective roughness length of the ground according to the effective roughness length of the ground and the reference value of the diffusion parameters:
(1) When the ground is effectiveLength of roughness Z 0 When the grain size is less than or equal to 0.1 m:
wherein, the first and the second end of the pipe are connected with each other,represents a downwind atmospheric diffusion parameter related to the effective roughness length of the ground,represents a lateral wind direction diffusion parameter related to the effective roughness length of the ground,represents a vertical diffusion parameter related to the effective roughness length of the ground,indicating a reference value of a lateral wind direction dispersion parameterA reference value of a vertical diffusion parameter is indicated,representing the effective roughness length of the ground.
(2) As the effective roughness length Z of the ground 0 When the particle size is more than or equal to 0.1 m:
wherein, the first and the second end of the pipe are connected with each other,a feature fit function representing the extent of diffusion and the magnitude of the velocity in the x-direction,the distance in the down-wind direction is indicated,the effective roughness length of the ground is expressed,、、、、、、which represents the coefficient of atmospheric stability of the air,a feature quantity fitting function representing the extent of diffusion in the y-direction and the magnitude of the velocity,a feature quantity fitting function representing the extent of diffusion and the magnitude of velocity in the z direction,represents a downwind atmospheric diffusion parameter related to the effective roughness length of the ground,represents a lateral wind direction diffusion parameter related to the effective roughness length of the ground,represents a vertical diffusion parameter related to the effective roughness length of the ground,a reference value of a lateral wind direction dispersion parameter is shown,a reference value of a vertical diffusion parameter is indicated,representing the ground effective roughness length.
Preferably: and 3, determining the atmospheric diffusion parameter related to the leakage time according to the leakage time and the diffusion parameter reference value:
wherein, the first and the second end of the pipe are connected with each other,indicating a downwind atmospheric diffusion parameter associated with the leak time,represents a lateral wind direction dispersion parameter related to the leakage time,a vertical diffusion parameter related to the leakage time is indicated,the time of the leak is indicated,a reference value of a lateral wind direction dispersion parameter is shown,indicating a vertical diffusion parameter reference value.
Preferably: and 4, the gas leakage diffusion model in the step 4 comprises a continuous point source gas leakage diffusion model and an instant gas leakage diffusion model.
Continuous point source gas leakage diffusion model:
in the formula:representing a continuous point source gas leak diffusion concentration,the coordinates of the spatial position are represented,xindicating the distance from the discharge point of the leakage source to any point in the downwind direction;yThe distance from the central axis of the smoke to any point in the right-angle horizontal direction is represented;zindicating the height from the surface of the earth to any point,representing the mass flow of material being discharged continuously,represents a lateral wind direction dispersion parameter related to the leakage time,a vertical diffusion parameter related to the leakage time is indicated,which is indicative of the average wind speed of the environment,indicating the effective height of the source of the leak.
Instantaneous gas leakage diffusion model:
wherein, the first and the second end of the pipe are connected with each other,indicating the instantaneous gas leak diffusion concentration,which is indicative of the mass of material being discharged instantaneously,which is indicative of the time of day,
the coordinates of the spatial position are represented,xthe distance from the discharge point of the leakage source to any point in the downwind direction is shown;ythe distance from the central axis of the smoke to any point in the right-angle horizontal direction is represented;zindicating the height from the surface of the earth to any point,indicating a downwind atmospheric diffusion parameter associated with the leak time,represents a lateral wind direction dispersion parameter related to the leakage time,a vertical diffusion parameter related to the leakage time is indicated,which is indicative of the average wind speed of the environment,indicating the effective height of the source of the leak.
Preferably, the following components: the method for determining the leakage diffusion hazard radius according to the gas leakage diffusion model in the step 5 comprises the following steps:
downwind distance:
wherein, the first and the second end of the pipe are connected with each other,the distance in the down-wind direction is indicated,which is indicative of the average wind speed of the environment,indicating the diffusion time.
Distance from downwindSubstituting the model into the gas leakage diffusion model, and making y =0, the gas leakage diffusion concentration C is obtained.
Leak diffusion hazard radius squared:
wherein the content of the first and second substances,represents the leak diffusion hazard radius squared (m) 2 ),Is the diffusion coefficient in the y-direction,the concentration of the hazardous substances is shown,indicating gas leak diffusion concentration
Leak diffusion hazard radius:
wherein, the first and the second end of the pipe are connected with each other,indicating the leak diffusion hazard radius.
The utility model provides a leakage protection distance survey system based on positioning monitoring, includes input module, diffusion parameter benchmark confirm module, atmospheric diffusion parameter confirm module relevant with the effective roughness length in ground, atmospheric diffusion parameter confirm module relevant with the leakage time, gas leakage diffusion model module, leakage diffusion harm radius confirm module, output module, wherein:
the input module is used for inputting meteorological conditions, geographic conditions and diffusion time.
The diffusion parameter reference value determining module is used for obtaining a diffusion parameter reference value by adopting a Ternal method according to meteorological conditions and geographic conditions.
The atmospheric diffusion parameter determination module related to the effective roughness length of the ground is used for determining atmospheric diffusion parameters related to the effective roughness length of the ground according to the effective roughness length of the ground and a diffusion parameter reference value.
And the atmospheric diffusion parameter determining module related to the leakage time is used for determining the atmospheric diffusion parameter related to the leakage time according to the leakage time and the diffusion parameter reference value.
The gas leakage diffusion model module is used for obtaining the gas leakage diffusion concentration by utilizing a gas leakage diffusion model according to the input meteorological condition, the diffusion time, the atmospheric diffusion parameter related to the effective roughness length of the ground and the atmospheric diffusion parameter related to the leakage time.
The leakage diffusion damage radius determining module is used for obtaining the leakage diffusion damage radius according to the gas leakage diffusion concentration.
Preferably, the following components: the output module is used for outputting the downwind damage influence areaDownwind hazard distanceCross wind direction hazard distanceAnd forming time of the hazard zone.
Compared with the prior art, the invention has the following beneficial effects:
according to the invention, the concentration of the leaked gas and the leakage diffusion damage radius can be obtained according to the meteorological condition, the geographic condition and the diffusion time, so that the problem of difficulty in acquiring field monitoring data is avoided.
Drawings
FIG. 1 is a schematic diagram of warning distances for generating leakage protection on a map.
Detailed Description
The present invention is further illustrated by the following description in conjunction with the accompanying drawings and the specific embodiments, it is to be understood that these examples are given solely for the purpose of illustration and are not intended as a definition of the limits of the invention, since various equivalent modifications will occur to those skilled in the art upon reading the present invention and fall within the limits of the appended claims.
A leakage protection distance measuring method based on positioning monitoring comprises the following steps:
step 1, obtaining a diffusion parameter reference value by adopting a Ternal method according to meteorological conditions and geographic conditions.
And dividing the atmospheric stability by adopting a Terner method and combining factors such as meteorological conditions, geographic conditions and the like. The atmospheric stability refers to the stability of a certain air mass in the atmosphere in the vertical direction, and affects the diffusion of pollutants in the atmosphere, and the more unstable the atmosphere, the faster the diffusion rate of pollutants. The atmospheric stability is closely related to weather phenomena, space-time scales and geographical conditions, and is divided by adopting a Tennel method and combining factors such as meteorological conditions, geographical conditions and the like. The Tenner method applies to the place, time and height of the sunθ h And cloud cover (the whole sky is 10 minutes), determining the solar radiation level, and determining the atmospheric stability level according to the solar radiation level and the average wind speed 10m away from the ground height.
Solar altitude angleθ h Calculated according to the following formula:
in the formula:φ-local geographical latitude.
λ-local geographical longitude.
t-Beijing time at observation.
δSun tilt angle, calculated using the following formula.
In the formula:θ 0 =2πd n /365;
d n date ordinal number in one year, 0,1,2, \8230;, 364.
According to Table 1, from the height of the sunθ h And cloud cover (the whole sky is a tenth system, the total cloud cover and the low cloud cover are determined by local meteorological observation data, and the solar radiation grade is determined by referring to chapter 4.3 of ground meteorological observation Specification of Central weather Bureau of China).
TABLE 1 solar radiation rating (Chinese)
According to Table 2, the atmospheric stability rating is determined from the solar irradiance level and the average wind speed at 10m above the ground (average wind speed at 10min above 10m above ground, using weather station data, see Chapter 9 of the Central weather agency "ground Meteorological Observation Specification").
TABLE 2 atmospheric stability rating
Diffusion parameter σ y 、σ z Is a standard deviation of a normal distribution function, which is a characteristic quantity indicating the diffusion range and the rate size. Considering the technical method for establishing local atmospheric pollutant emission standard (GB/T-3840-91) according to national standardThe effect of the ground conditions on the diffusion was corrected for the atmospheric stability rating determined by the tanner method as follows:
(1) rural and urban suburban areas in plain areas: A. directly checking the table for determining the B-level stability and the C-level stability; D. e and F levels of stability need to be determined by half-level lookup to the unstable direction;
(2) in rural areas or cities in urban areas of industrial areas and hilly mountain areas, A and B grades are not mentioned, C grade is mentioned, D, E and F are mentioned in the unstable direction by one and a half grades, and then the table is looked up for determination; in cities in non-industrial areas, A and B grades are not mentioned, C grade is mentioned from B to C grade, D, E and F are mentioned in unstable directions, and then the urban cities are determined by looking up a table. The details are shown in Table 3.
TABLE 3 different area diffusion parameter stability level correction
The diffusion parameter calculation formula is as follows:
wherein the content of the first and second substances,represents a reference value of a lateral wind direction dispersion parameter,、representing lateral diffusion parameter power function coefficients;、representing the coefficients of a power function of a vertical diffusion parameter,the distance in the down-wind direction is indicated,indicating a vertical diffusion parameter reference value.
TABLE 4 coefficient value of power function expression of transverse diffusion parameter (sampling time 0.5 h)
TABLE 5 vertical diffusion parameter Power function expression coefficient values (sampling time 0.5 h)
According to the diffusion parameter power function expression coefficient values of tables 4 and 5, according to the downwind distancexDetermining the diffusion coefficient gamma 1 、a 1 、γ 2 、a 2 Determining diffusion parameters based on the established empirical function、。
And 2, determining atmospheric diffusion parameters related to the effective roughness length of the ground according to the effective roughness length of the ground and the diffusion parameter reference value.
(1) As the effective roughness length Z of the ground 0 When the grain size is less than or equal to 0.1 m:
wherein the content of the first and second substances,represents a downwind atmospheric diffusion parameter related to the effective roughness length of the ground,represents a lateral wind direction diffusion parameter related to the effective roughness length of the ground,represents a vertical diffusion parameter related to the effective roughness length of the ground,indicating a reference value of a lateral wind direction dispersion parameterA reference value of a vertical diffusion parameter is indicated,representing the ground effective roughness length.
(2) As the effective roughness length Z of the ground 0 When the particle size is more than or equal to 0.1 m:
wherein, the first and the second end of the pipe are connected with each other,a feature quantity fitting function representing the extent of diffusion and the magnitude of the velocity in the x-direction,the distance in the downwind direction is indicated,the effective roughness length of the ground is represented,、、、、、、which represents the coefficient of stability of the atmosphere,a feature quantity fitting function representing the extent of diffusion in the y-direction and the magnitude of the velocity,features indicating extent of diffusion in the z-direction and magnitude of velocityThe characteristic quantity fitting function is used for fitting the characteristic quantity,represents a downwind atmospheric diffusion parameter related to the effective roughness length of the ground,represents a lateral wind direction dispersion parameter related to the effective roughness length of the ground,represents a vertical diffusion parameter related to the effective roughness length of the ground,represents a reference value of a lateral wind direction dispersion parameter,a reference value of a vertical diffusion parameter is indicated,representing the effective roughness length of the ground.
And 3, determining the atmospheric diffusion parameter related to the leakage time according to the leakage time and the diffusion parameter reference value.
Wherein the content of the first and second substances,indicating a downwind atmospheric diffusion parameter related to the leak time,presentation versus leakage timeThe off-side wind direction spread parameter,a vertical diffusion parameter related to the leakage time is indicated,the time to leak is indicated by the time to leak,represents a reference value of a lateral wind direction dispersion parameter,indicating a vertical diffusion parameter reference value.
And 4, establishing a gas leakage diffusion model according to the wind speed, the wind direction, the leakage source, the atmospheric diffusion parameter related to the effective roughness length of the ground and the atmospheric diffusion parameter related to the leakage time.
The gas leakage diffusion model comprises a continuous point source gas leakage diffusion model and an instant gas leakage diffusion model, and the gas leakage diffusion model calculates leakage concentration data of different leakage sources under different conditions (mainly considering wind speed, wind direction, leakage source strength and the like) in a forward direction.
And establishing a continuous point source gas leakage diffusion model according to the wind speed, the wind direction, the leakage source and the atmospheric diffusion parameter related to the effective roughness length of the ground.
Continuous point source gas leakage diffusion model:
in the formula:representing points of continuityThe source gas leaks out of the diffusion concentration,the coordinates of the spatial position are represented,xthe distance (m) from the discharge point of the leakage source to any point in the downwind direction is shown;yrepresents the distance (m) from the central axis of the flue gas to any point in the right-angle horizontal direction;zrepresents the height (m) from the surface to any point;representing the mass flow of material being discharged continuously,represents a lateral wind direction dispersion parameter related to the leakage time,a vertical diffusion parameter related to the leakage time is indicated,、effective roughness length Z with the ground in continuous leakage 0 In connection with this, the present invention is,、when calculating, the atmospheric diffusion coefficient (calculated according to Z0 ≦ 0.1、) Reuse of the effective roughness Z of the ground 0 Whether it is greater than 0.1 or not,representing the ambient mean wind speed (m/s),indicating the effective height of the source of the leak.
Calculate concentration at z = Hr:
in the formula:indicating the leakage source lift height (m) (out of the "significant hazard classification criteria. Doc"),representing the gas cloud exit velocity (m/s),the diameter (m) of the outlet is shown,representing the ambient wind speed (m/s),represents the effective height (m) of the leakage source,indicating the leakage source height (m).
The calculation process for determining y in the formula is as follows:
a) Step size selection (solving the iterative performance problem, setting target concentration n) d (mg/m 3 ) Selecting a step size according to the target concentration,
when n is d <=0.01, the step size may be 1000;
when n is d <If =0.1, the step size may be 700;
when n is d <If =1, the step size may be 300; when n is d <=50, the step size may be taken as 100;
when n is d <=100, then the step size may be 00; when n is d <=300, the step size may be 1;
b) The downwind distance x =0 is iteratively increased in steps, and for a fixed x, the above equation may be usedThe reverse-push formula of (1) calculates the y value corresponding to X (i.e. corresponding to the cross wind direction distance, at which the substituted concentration is the poisoning concentration of the critical value), records the maximum cross wind direction distance under X which has been traversed at present by using one value, and records the X corresponding to the critical value at which a certain position y begins to become smaller as X increases as X is X mid In this case, y is the farthest distance affected by the transverse wind direction. The corresponding cross-wind distance is then calculated as x increases until y =0, calculated at the corresponding x, when x is the downwind poisoning hazard farthest distance. Calculating the area of poisoning hazard in downwind direction: dividing the hazard region into n regions according to step length, calculating the area of each region, and superposing, wherein when the step length is greater than 1, calculating the area of the current small region, the step length should be specifically calculated in a thinning manner to obtain
and (3) introducing time factors for analysis, and being suitable for instantaneous gas leakage diffusion, namely establishing an instantaneous gas leakage diffusion model according to wind speed, wind direction, leakage sources and atmospheric diffusion parameters related to leakage time.
Instantaneous gas leakage diffusion model:
wherein, the first and the second end of the pipe are connected with each other,indicating the instantaneous gas leak diffusion concentration,which is indicative of the mass of material being discharged instantaneously,which is indicative of the time of day,
the coordinates of the spatial position are represented,xthe distance (m) from the discharge point of the leakage source to any point in the downwind direction is shown;yrepresents the distance (m) from the central axis of the flue gas to any point in the right-angle horizontal direction;zrepresenting the height (m) from the surface to any point,representing the mass flow of material being discharged continuously,indicating a downwind atmospheric diffusion parameter related to the leak time,represents a lateral wind direction dispersion parameter related to the leakage time,a vertical diffusion parameter related to the leakage time is indicated,which represents the average wind speed of the environment,indicating the effective height of the source of the leak.
Length Z of effective roughness of ground surface during instantaneous leakage of atmospheric diffusion parameter 0 Independently, the calculated atmospheric diffusion parameter (calculated as Z0 ≦ 0.1) is reused for the leakage time T l It is corrected as shown in step 3.
the calculation process for determining the farthest distance of influence L of downwind is as follows:
from the above, the atmospheric diffusion parameter is related to x (x = ut), and the iteratively calculated concentration using time is compared with the target concentration until the calculated concentration<The time obtained by the target concentration (kg/m 3) is the farthest downwind poisoning influence forming time(s), and then the farthest downwind influence distance (m) L is calculated max =T max * u; where u represents the ambient mean wind speed (m/s).
And 5, determining the leakage diffusion damage radius according to the gas leakage diffusion model.
The method for determining the leak diffusion hazard radius according to the gas leak diffusion model is as follows:
downwind distance:
wherein the content of the first and second substances,the distance in the down-wind direction is indicated,which represents the average wind speed of the environment,indicating the diffusion time.
Distance from downwindSubstituting the model into the gas leakage diffusion model, and making y =0, the gas leakage diffusion concentration C is obtained.
Leak diffusion hazard radius squared:
wherein, the first and the second end of the pipe are connected with each other,represents the leakage spread hazard radius squared (m) 2 ),Is the diffusion coefficient in the y-direction,the concentration of the hazardous substances is shown,indicating the gas leak diffusion concentration.
Leak diffusion hazard radius:
When R >0, the output result is:
1. diffusion T k After second, the hazard radius at X meters downwind is r meters.
3. Diffusion hazard conditions (shown at the simulation interface).
3.4 Hazard zone formation time (seconds): t is k 。
When R < =0, then the resulting parameters are all 0.
The utility model provides a leakage protection distance survey system based on positioning monitoring, includes input module, diffusion parameter benchmark confirm module, atmospheric diffusion parameter confirm module relevant with the effective roughness length in ground, atmospheric diffusion parameter confirm module relevant with the leakage time, gas leakage diffusion model module, leakage diffusion harm radius confirm module, output module, wherein:
the input module is used for inputting meteorological conditions, geographic conditions and diffusion time.
The diffusion parameter reference value determining module is used for obtaining a diffusion parameter reference value by adopting a Ternal method according to meteorological conditions and geographic conditions.
The atmospheric diffusion parameter determining module related to the ground effective roughness length is used for determining atmospheric diffusion parameters related to the ground effective roughness length according to the ground effective roughness length and a diffusion parameter reference value.
And the atmospheric diffusion parameter determining module related to the leakage time is used for determining the atmospheric diffusion parameter related to the leakage time according to the leakage time and the diffusion parameter reference value.
The gas leakage diffusion model module is used for obtaining gas leakage diffusion concentration by utilizing a gas leakage diffusion model according to input meteorological conditions, diffusion time, atmospheric diffusion parameters related to the effective roughness length of the ground and atmospheric diffusion parameters related to the leakage time.
The leakage diffusion hazard radius determination module is used for obtaining a leakage diffusion hazard radius according to the gas leakage diffusion concentration.
The output module is used for outputting the leakage diffusion damage radiusArea affected by downwind damageDownwind hazard distanceCross wind direction hazard distanceAnd forming time of the hazard zone.
Selecting a dangerous chemical leakage scene, combining the monitoring results of the fixed sensor and the mobile control monitoring results, comparing a group of monitoring values with the concentration value calculated by the model, and when the standard deviation of the concentration is metLess than 10 -5 The position of a leakage source (release source) is located by back calculation, the relevant hazard distance calculated by the model is combined, and the ERG definition data is combined to generate a warning distance graph of leakage protection on a map, as shown in figure 1, so that the measurement of the leakage protection distance and the division of a warning isolation area are realized.
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and these are intended to be within the scope of the invention.
Claims (8)
1. A leakage protection distance measuring method based on positioning monitoring is characterized by comprising the following steps:
step 1, obtaining a diffusion parameter reference value by adopting a Ternal method according to meteorological conditions and geographic conditions;
step 2, determining atmospheric diffusion parameters related to the effective roughness length of the ground according to the effective roughness length of the ground and the reference value of the diffusion parameters;
step 3, determining atmospheric diffusion parameters related to the leakage time according to the leakage time and the diffusion parameter reference value;
step 4, establishing a gas leakage diffusion model according to the wind speed, the wind direction, the leakage source, the atmospheric diffusion parameter related to the effective roughness length of the ground and the atmospheric diffusion parameter related to the leakage time;
and 5, determining the leakage diffusion damage radius according to the gas leakage diffusion model.
2. The method for determining the leakage protection distance based on the positioning monitoring as claimed in claim 1, wherein: in the step 1, a method for obtaining a diffusion parameter reference value by adopting a Ternal method according to meteorological conditions and geographic conditions comprises the following steps:
wherein, the first and the second end of the pipe are connected with each other,represents a reference value of a lateral wind direction dispersion parameter,、representing lateral diffusion parameter power function coefficients;、representing the coefficients of a power function of a vertical diffusion parameter,the distance in the down-wind direction is indicated,indicating a vertical diffusion parameter reference value.
3. The method for determining the distance to prevent leakage based on location monitoring as claimed in claim 2, wherein: in step 2, determining an atmospheric diffusion parameter related to the effective roughness length of the ground according to the effective roughness length of the ground and a diffusion parameter reference value:
(1) As the effective roughness length Z of the ground 0 When the grain size is less than or equal to 0.1 m:
wherein the content of the first and second substances,represents a downwind atmospheric diffusion parameter related to the effective roughness length of the ground,represents a lateral wind direction diffusion parameter related to the effective roughness length of the ground,represents a vertical diffusion parameter related to the effective roughness length of the ground,indicating a reference value of a lateral wind direction dispersion parameterRepresents a reference value of a vertical diffusion parameter,representing the effective roughness length of the ground;
(2) As the effective roughness length Z of the ground 0 When the particle size is more than or equal to 0.1 m:
wherein the content of the first and second substances,a feature fit function representing the extent of diffusion and the magnitude of the velocity in the x-direction,the distance in the downwind direction is indicated,the effective roughness length of the ground is represented,、、、、、、which represents the coefficient of atmospheric stability of the air,a feature fit function representing the extent of diffusion in the y-direction and the magnitude of the velocity,a feature quantity fitting function representing the extent of diffusion and the magnitude of velocity in the z direction,represents a downwind atmospheric diffusion parameter related to the effective roughness length of the ground,represents a lateral wind direction dispersion parameter related to the effective roughness length of the ground,represents a vertical diffusion parameter related to the effective roughness length of the ground,represents a reference value of a lateral wind direction dispersion parameter,a reference value of a vertical diffusion parameter is indicated,representing the effective roughness length of the ground.
4. The method for determining the leakage protection distance based on the positioning monitoring as claimed in claim 3, wherein: and 3, determining an atmospheric diffusion parameter related to the leakage time according to the leakage time and the diffusion parameter reference value:
wherein, the first and the second end of the pipe are connected with each other,indicating a downwind atmospheric diffusion parameter associated with the leak time,represents a lateral wind direction dispersion parameter related to the leakage time,a vertical diffusion parameter related to the leakage time is indicated,the time of the leak is indicated,represents a reference value of a lateral wind direction dispersion parameter,indicating a vertical diffusion parameter reference value.
5. The method for determining the leakage protection distance based on the positioning monitoring as claimed in claim 4, wherein: the gas leakage diffusion model in the step 4 comprises a continuous point source gas leakage diffusion model and an instant gas leakage diffusion model;
continuous point source gas leakage diffusion model:
in the formula:representing a continuous point source gas leak diffusion concentration,the coordinates of the spatial position are represented,xindicating the distance from the leakage source discharge point to any point up the downwind,yindicating the distance of the central axis of the flue gas to any point in the right-angle horizontal direction,zindicating the height from the surface of the earth to any point,representing the mass flow of material being discharged continuously,represents a lateral wind direction dispersion parameter related to the leakage time,a vertical diffusion parameter related to the leakage time is indicated,which is indicative of the average wind speed of the environment,represents the effective height of the leakage source;
instantaneous gas leakage diffusion model:
wherein the content of the first and second substances,indicating the instantaneous gas leak diffusion concentration,which is indicative of the mass of material being discharged instantaneously,which is indicative of the time of day,the coordinates of the spatial position are represented,xindicating the distance from the leakage source discharge point to any point up the downwind,yindicating the distance of the central axis of the flue gas to any point in the right-angle horizontal direction,zindicating the height from the surface of the earth to any point,indicating a downwind atmospheric diffusion parameter related to the leak time,represents a lateral wind direction dispersion parameter related to the leakage time,a vertical diffusion parameter related to the leakage time is indicated,which represents the average wind speed of the environment,indicating the effective height of the source of the leak.
6. The method for determining the distance between a leak protection device and a base station based on positioning monitoring according to claim 5, wherein: the method for determining the leakage diffusion hazard radius according to the gas leakage diffusion model in the step 5 comprises the following steps:
downwind distance:
wherein, the first and the second end of the pipe are connected with each other,the distance in the down-wind direction is indicated,which represents the average wind speed of the environment,represents the diffusion time;
distance from downwindSubstituting the model into a gas leakage diffusion model, and enabling y =0 to obtain a gas leakage diffusion concentration C;
leakage spread hazard radius squared:
wherein, the first and the second end of the pipe are connected with each other,indicating the leak spread compromise radius squared,is a y directionThe upward diffusion coefficient of the light is,the concentration of the hazardous substances is shown,indicating a gas leak diffusion concentration;
leak diffusion hazard radius:
7. A leakage protection distance measuring system based on the leakage protection distance measuring method based on positioning monitoring of claim 1, characterized in that: the device comprises an input module, a diffusion parameter reference value determining module, an atmospheric diffusion parameter determining module related to the effective roughness length of the ground, an atmospheric diffusion parameter determining module related to the leakage time, a gas leakage diffusion model module, a leakage diffusion hazard radius determining module and an output module, wherein:
the input module is used for inputting meteorological conditions, geographic conditions and diffusion time;
the diffusion parameter reference value determining module is used for obtaining a diffusion parameter reference value by adopting a Ternal method according to meteorological conditions and geographic conditions;
the atmospheric diffusion parameter determining module related to the ground effective roughness length is used for determining atmospheric diffusion parameters related to the ground effective roughness length according to the ground effective roughness length and a diffusion parameter reference value;
the atmospheric diffusion parameter determining module related to the leakage time is used for determining atmospheric diffusion parameters related to the leakage time according to the leakage time and the diffusion parameter reference value;
the gas leakage diffusion model module is used for obtaining gas leakage diffusion concentration by utilizing a gas leakage diffusion model according to input meteorological conditions, diffusion time, atmospheric diffusion parameters related to the effective roughness length of the ground and atmospheric diffusion parameters related to the leakage time;
the leakage diffusion damage radius determining module is used for obtaining a leakage diffusion damage radius according to the gas leakage diffusion concentration;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211331979.1A CN115638938A (en) | 2022-10-28 | 2022-10-28 | Leakage protection distance measuring system and method based on positioning monitoring |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211331979.1A CN115638938A (en) | 2022-10-28 | 2022-10-28 | Leakage protection distance measuring system and method based on positioning monitoring |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115638938A true CN115638938A (en) | 2023-01-24 |
Family
ID=84945881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211331979.1A Pending CN115638938A (en) | 2022-10-28 | 2022-10-28 | Leakage protection distance measuring system and method based on positioning monitoring |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115638938A (en) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02190733A (en) * | 1989-01-20 | 1990-07-26 | Jgc Corp | Method for predicting diffusion area of leakage gas |
JPH06117600A (en) * | 1992-09-30 | 1994-04-26 | Mitsui Toatsu Chem Inc | System for estimating leakage point and leakage volume of gas, steam or the like |
JPH07198524A (en) * | 1993-12-28 | 1995-08-01 | Chiyoda Corp | Estimation method for gas leakage point and amount from gas concentration data |
KR20000002127A (en) * | 1998-06-17 | 2000-01-15 | 최인영 | Method of modeling diffusion of toxic gas |
CN102411669A (en) * | 2011-07-27 | 2012-04-11 | 南京林业大学 | Prediction method of leakage rate of contact mechanical seal |
CN104834017A (en) * | 2015-04-30 | 2015-08-12 | 华南理工大学 | Positioning method for hazardous gas leakage diffusion accident source |
CN105372388A (en) * | 2015-10-14 | 2016-03-02 | 北京千安哲信息技术有限公司 | Gas concentration monitoring method and device for finite time steady leakage |
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 |
US20170176182A1 (en) * | 2012-04-20 | 2017-06-22 | Picarro, Inc. | Leak position estimation method using wind direction and correlation measurements at isolated spatial points in a plume |
CN108508147A (en) * | 2018-03-13 | 2018-09-07 | 中国特种设备检测研究院 | A kind of measuring system and measurement method of liquid gas steady state leakage rate |
CN108932394A (en) * | 2018-07-19 | 2018-12-04 | 成都鼎安华智慧物联网股份有限公司 | Toxic gas leakage hazard method of determining range |
CN109029847A (en) * | 2018-09-29 | 2018-12-18 | 中国石油大学(华东) | Gas leakage monitoring system, FLNG system and its gas leakage monitoring method |
CN112345201A (en) * | 2020-11-10 | 2021-02-09 | 中国石油大学(北京) | Gas leakage detection method and system |
CN112529240A (en) * | 2020-09-14 | 2021-03-19 | 桂林电子科技大学 | Prediction method, system, device and storage medium of atmospheric environment data |
CN114896901A (en) * | 2022-04-25 | 2022-08-12 | 杭州立方控股股份有限公司 | Hazardous chemical atmospheric diffusion assessment method |
-
2022
- 2022-10-28 CN CN202211331979.1A patent/CN115638938A/en active Pending
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH02190733A (en) * | 1989-01-20 | 1990-07-26 | Jgc Corp | Method for predicting diffusion area of leakage gas |
JPH06117600A (en) * | 1992-09-30 | 1994-04-26 | Mitsui Toatsu Chem Inc | System for estimating leakage point and leakage volume of gas, steam or the like |
JPH07198524A (en) * | 1993-12-28 | 1995-08-01 | Chiyoda Corp | Estimation method for gas leakage point and amount from gas concentration data |
KR20000002127A (en) * | 1998-06-17 | 2000-01-15 | 최인영 | Method of modeling diffusion of toxic gas |
CN102411669A (en) * | 2011-07-27 | 2012-04-11 | 南京林业大学 | Prediction method of leakage rate of contact mechanical seal |
US20170176182A1 (en) * | 2012-04-20 | 2017-06-22 | Picarro, Inc. | Leak position estimation method using wind direction and correlation measurements at isolated spatial points in a plume |
CN104834017A (en) * | 2015-04-30 | 2015-08-12 | 华南理工大学 | Positioning method for hazardous gas leakage diffusion accident source |
CN105372388A (en) * | 2015-10-14 | 2016-03-02 | 北京千安哲信息技术有限公司 | Gas concentration monitoring method and device for finite time steady leakage |
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 |
CN108508147A (en) * | 2018-03-13 | 2018-09-07 | 中国特种设备检测研究院 | A kind of measuring system and measurement method of liquid gas steady state leakage rate |
CN108932394A (en) * | 2018-07-19 | 2018-12-04 | 成都鼎安华智慧物联网股份有限公司 | Toxic gas leakage hazard method of determining range |
CN109029847A (en) * | 2018-09-29 | 2018-12-18 | 中国石油大学(华东) | Gas leakage monitoring system, FLNG system and its gas leakage monitoring method |
CN112529240A (en) * | 2020-09-14 | 2021-03-19 | 桂林电子科技大学 | Prediction method, system, device and storage medium of atmospheric environment data |
CN112345201A (en) * | 2020-11-10 | 2021-02-09 | 中国石油大学(北京) | Gas leakage detection method and system |
CN114896901A (en) * | 2022-04-25 | 2022-08-12 | 杭州立方控股股份有限公司 | Hazardous chemical atmospheric diffusion assessment method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ansmann et al. | Extreme levels of Canadian wildfire smoke in the stratosphere over central Europe on 21-22 August 2017 | |
CN106023530B (en) | A kind of heavy rain type diluted debris flow monitoring and prediction method for early warning | |
Liu et al. | Aircraft study of aerosol vertical distributions over Beijing and their optical properties | |
Eckhardt et al. | Estimation of the vertical profile of sulfur dioxide injection into the atmosphere by a volcanic eruption using satellite column measurements and inverse transport modeling | |
CN104834017A (en) | Positioning method for hazardous gas leakage diffusion accident source | |
CN103914622A (en) | Quick chemical leakage predicating and warning emergency response decision-making method | |
CN108932394A (en) | Toxic gas leakage hazard method of determining range | |
KR102443718B1 (en) | System for producing 3-dimensional diffusion prediction map based on cyber physical for monitoring air pollution gas diffusion of environmental facilities, and method for the same | |
CN108133514B (en) | Atmospheric influence evaluation method for airborne emission of liquid radioactive effluent of nuclear power plant | |
CN111339665A (en) | Troposphere ozone profile calculation method | |
CN107367773A (en) | A kind of single-point meteorological element based on Big Dipper positioning determines method | |
Ivanov et al. | On the possibility of estimating the volume of NO2 emissions in cities using zenith spectral observations of diffuse solar radiation near 450 nm | |
CN105372388B (en) | For the gas concentration monitoring method and device of finite time steady state leakage | |
CN101739790B (en) | Forecasting, early warning and emergency controlling method for a plurality of fixed chemical risk sources | |
Zhang et al. | Study on optimization of shelter locations and evacuation routes of gas leakage accidents in chemical industrial park | |
KR20160078920A (en) | Method and apparatus for calculating irradiance of target point | |
CN113313291B (en) | Power grid wind disaster ensemble forecasting and verifying method based on initial disturbance | |
Al-Rashidi et al. | Investigation of the efficiency of existing air pollution monitoring sites in the state of Kuwait | |
CN115638938A (en) | Leakage protection distance measuring system and method based on positioning monitoring | |
CN116413838A (en) | Urban pollutant diffusion rapid early warning, prevention and control platform | |
US20100042332A1 (en) | System and method for source identification for a chemical release | |
CN114777030A (en) | Hazardous gas monitoring method based on NB-IOT technology | |
Kang et al. | Wind resource assessment and potential development of wind farms along the entire coast of South Korea using public data from the Korea meteorological administration | |
CN112016738B (en) | Harm prediction method for inhalable toxic substances | |
Zhou et al. | Source characterization of airborne emissions using a sensor network: Examining the impact of sensor quality, quantity, and wind climatology |
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 |