CN105372388B - For the gas concentration monitoring method and device of finite time steady state leakage - Google Patents

Abstract

The invention discloses gas concentration monitoring method and device for finite time steady state leakage.Method includes：Gather the leakage intensity of coordinate and source of leaks of the monitoring point on source of leaks and geographic area in geographic coordinate system；The mean wind direction angle of coordinate and geographic area according to source of leaks and monitoring point in geographic coordinate system, determines coordinate of the monitoring point in target-based coordinate system；The floor type of coordinate and geographic area of the meteorologic parameter, leakage time of origin, monitoring point of coordinate, geographic area according to source of leaks in geographic coordinate system in target-based coordinate system, determines the actual dispersion coefficient at monitoring point；The gas of continuous release is separated into gas cigarette group, and determines the concentration that each gas cigarette group is formed at monitoring point；According to the concentration that each gas cigarette group is formed at monitoring point, the gas leakage concentration that monitoring point is in total monitoring duration is determined.Thus, amount of calculation can be reduced, the distribution of gas leakage concentration is obtained in time.

Description

For the gas concentration monitoring method and device of finite time steady state leakage

Technical field

The present invention relates to gas leakage concentration monitor field, in particular it relates to a kind of for finite time steady state leakage Gas concentration monitoring method and device.

Background technology

With the continuous progress of industrial technology, the mankind will face increasing Safety production questions, disposal of pollutants problem Deng gas (particularly toxic gas) leakage problem is then one of them.The toxic gas of leakage can be fast under the action of the forces of the wind Speed is propagated, and causes large area personnel poisoning, injures and deaths, even results in great, serious accident.

The leakage of toxic gas can be divided into continuous steady state leakage, finite time steady state leakage, transient state leakage three types. Wherein, finite time steady state leakage refers to：The quantity of the toxic gas of leakage is limited, after leakage there occurs certain hour Leakage stops, or during leaking, the situation of leakage etc. is stopped by manual control.For finite time steady state leakage problem, Currently used mainly integration cigarette group model.It is mainly using integration cigarette group model and will leak out gas by between the regular hour Every independent cigarette group is divided into, the change of monitoring gas concentration in the space that each cigarette group causes is investigated, then take the time The mode of integration obtains the distribution of SS gas concentration after finite time leakage.However, for integration cigarette group model, More long when the monitoring time, monitoring time interval is smaller, when monitoring geographical coverage area is larger, can be brought using integration cigarette group model The very big problem of amount of calculation is (for example, the amount of calculation in integration method forms quadratic power and increases with the increase for solving time step number It is long), so as to cause computational efficiency low, it is impossible to obtain the situation of gas leakage concentration distribution in time.

The content of the invention

Amount of calculation can be saved, improve computational efficiency and let out for finite time stable state it is an object of the invention to provide a kind of The gas concentration monitoring method and device of leakage.

To achieve these goals, the present invention provides a kind of gas concentration monitoring side for finite time steady state leakage Method, the method includes：Gather coordinate of the monitoring point on source of leaks and geographic area in geographic coordinate system, and the leakage The leakage intensity in source；Coordinate and the geography according to the source of leaks and the monitoring point in the geographic coordinate system The mean wind direction angle in region, determines coordinate of the monitoring point in target-based coordinate system, wherein, the target-based coordinate system is with institute Source of leaks vertical subpoint earthward is stated for origin, is the positive coordinate system of x-axis with the mean wind direction of the geographic area；Root It is coordinate, the meteorologic parameter of the geographic area according to the source of leaks in the geographic coordinate system, leakage time of origin, described The floor type of coordinate and the geographic area of the monitoring point in the target-based coordinate system, determines at the monitoring point Actual dispersion coefficient；The gas of continuous release is separated into gas cigarette group, and during leakage intensity, monitoring according to the source of leaks Between interval, always monitoring duration, the seat of actual dispersion coefficient, the monitoring point at the monitoring point in the target-based coordinate system The effective depth of mark, the wind speed of the geographic area and the source of leaks, determines each gas cigarette group shape at the monitoring point Into concentration；According to the concentration that described each gas cigarette group is formed at the monitoring point, determine that the monitoring point is in chief inspector The gas leakage concentration surveyed in duration.

The present invention also provides a kind of gas concentration monitoring device for finite time steady state leakage, and the device includes：Adopt Collection module, for gathering coordinate of the monitoring point on source of leaks and geographic area in geographic coordinate system, and the source of leaks Leakage intensity；Coordinate determining module, for the seat according to the source of leaks and the monitoring point in the geographic coordinate system Mark and the mean wind direction angle of the geographic area, determine coordinate of the monitoring point in target-based coordinate system, wherein, it is described Target-based coordinate system is that the mean wind direction with the geographic area is as x with the source of leaks vertical subpoint earthward as origin The positive coordinate system of axle；Diffusion coefficient determining module, for the coordinate according to the source of leaks in the geographic coordinate system, institute State the meteorologic parameter of geographic area, coordinate in the target-based coordinate system of leakage time of origin, the monitoring point and described The floor type of geographic area, determines the actual dispersion coefficient at the monitoring point；First concentration determination module, for by continuously The gas of leakage is separated into gas cigarette group, and leakage intensity according to the source of leaks, monitoring time interval, always monitoring duration, Actual dispersion coefficient, the monitoring point at monitoring point coordinate in the target-based coordinate system, the geographic area The effective depth of wind speed and the source of leaks, determines the concentration that each gas cigarette group is formed at the monitoring point；Second concentration Determining module, for the concentration formed at the monitoring point according to each gas cigarette group, determines that the monitoring point is in Gas leakage concentration in total monitoring duration.

By above-mentioned technical proposal, compared to integration cigarette group model, amount of calculation can be greatly reduced.By by time domain Integral Problem is converted into the superposition problem of spatial concentration, can greatly save the calculating time, accelerates calculating speed, improves and calculates effect Rate, so as to ensure after generation is leaked, can in time obtain the situation of gas leakage concentration distribution, and then can be toxic gas The work such as rescue, evacuation after leakage try to gain time precious to one.

Other features and advantages of the present invention will be described in detail in subsequent specific embodiment part.

Brief description of the drawings

Accompanying drawing is, for providing a further understanding of the present invention, and to constitute the part of specification, with following tool Body implementation method is used to explain the present invention together, but is not construed as limiting the invention.In the accompanying drawings：

Fig. 1 shows the gas concentration monitoring method for finite time steady state leakage according to the embodiment of the present invention Flow chart.

Fig. 2 shows the gas concentration monitoring for finite time steady state leakage of another implementation method of the invention The flow chart of method.

Fig. 3 shows the gas concentration monitoring for finite time steady state leakage of another implementation method of the invention The flow chart of method.

Fig. 4 shows the gas concentration monitoring for finite time steady state leakage of another implementation method of the invention The flow chart of method.

Fig. 5 shows the gas concentration monitoring device for finite time steady state leakage according to the embodiment of the present invention Structured flowchart.

Specific embodiment

Specific embodiment of the invention is described in detail below in conjunction with accompanying drawing.It should be appreciated that this place is retouched The specific embodiment stated is merely to illustrate and explain the present invention, and is not intended to limit the invention.

Fig. 1 shows the gas concentration monitoring method for finite time steady state leakage according to the embodiment of the present invention Flow chart.As shown in figure 1, the method can include：

In step S101, coordinate of the monitoring point on source of leaks and geographic area in geographic coordinate system is gathered, and The leakage intensity of source of leaks.Wherein, the scope of geographic area can be preset, for example, can be according to actual monitoring need Ask and set.

Coordinate of the source of leaks in geographic coordinate system can be expressed as (Lo₀, La₀), wherein, Lo₀Represent the longitude of source of leaks (unit is degree, °), La₀Represent the latitude (unit for degree, °) of source of leaks.The effective depth of source of leaks can use H₀Represent (unit For m).Lo₀、La₀Can be measured by GPS (global positioning system) receiver, or can be positioned by mobile phone and map Carry out numerical value extraction.H₀Can actual height based on source of leaks calculate, computational methods can with for example, see：《Chemical process Safety theory and application (second edition)》, Jiang Juncheng etc. translates.

Coordinate of the monitoring point in geographic coordinate system can be expressed as (Lo_ij, La_ij), wherein, Lo_ijRepresent longitude (unit It is degree, °), La_ijExpression latitude (unit is degree, °), Lo_ij, La_ijCan be measured by GPS, or can led to Crossing mobile phone positioning and map carries out numerical value extraction.

At least one monitoring point can be included in geographic area.Can be measured by GPS (or can pass through Mobile phone is positioned and map carries out numerical value and extracts) coordinate of each monitoring point in geographic coordinate system in geographic area, i.e. can To obtain the longitude and latitude of each monitoring point.The longitude of each monitoring point can constitute a second order longitude square of M × N Battle array Lo, i.e.

Wherein, i=1,2 ... ..., M；J=1,2 ... ..., N, M × N can represent the total of the monitoring point on geographic area Number.Similarly, the latitude of each monitoring point can constitute a second order latitude matrix L a of M × N, i.e.

In step s 102, according to source of leaks and monitoring point in geographic coordinate system coordinate and geographic area it is flat Equal wind angle, determines coordinate of the monitoring point in target-based coordinate system, wherein, the target-based coordinate system is the lead with source of leaks earthward Straight subpoint is origin, is the positive coordinate system of x-axis with the mean wind direction of geographic area.

By the step for, coordinate (Lo that can be by monitoring point in geographic coordinate system_ij, La_ij) be converted to and sat in target Coordinate (X in mark system_ij, Y_ij).Also, x-axis coordinate and y-axis coordinate of the source of leaks in target-based coordinate system are 0.

In step s 103, according to source of leaks in geographic coordinate system coordinate, the meteorologic parameter of geographic area, leakage hair The floor type of the coordinate and geographic area of raw time, monitoring point in target-based coordinate system, determines the actual expansion at monitoring point Dissipate coefficient.Wherein, the meteorologic parameter can for example including total amount of cloud, low cloud cover, wind speed etc. parameter.The leakage time of origin can be with For example including information such as the leakage moment for occurring, the dates for leaking generation.

In step S104, the gas of continuous release is separated into gas cigarette group, and leakage intensity, prison according to source of leaks Survey time interval, always monitoring duration, coordinate, the geography of actual dispersion coefficient, monitoring point at monitoring point in target-based coordinate system The wind speed in region and the effective depth of source of leaks, determine concentration (that is, each gas that each gas cigarette group is formed at monitoring point Contribution of the cigarette group to the concentration of monitoring point).

In step S105, according to the concentration that each gas cigarette group is formed at monitoring point, determine that monitoring point is in chief inspector The gas leakage concentration surveyed in duration.

By above-mentioned technical proposal, compared to integration cigarette group model, amount of calculation can be greatly reduced.By by time domain Integral Problem is converted into the superposition problem of spatial concentration, can greatly save the calculating time, accelerates calculating speed, improves and calculates effect Rate, so as to ensure after generation is leaked, can in time obtain the situation of gas leakage concentration distribution, and then can be toxic gas The work such as rescue, evacuation after leakage try to gain time precious to one.

Fig. 2 shows the gas concentration monitoring for finite time steady state leakage of another implementation method of the invention The flow chart of method.As shown in Fig. 2 the step S102 can include：

In step s 201, according to monitoring point in geographic coordinate system coordinate, the mean sea level of geographic area, Yi Ji One geographical coordinate and the transfer function of global coordinates, determine coordinate of the monitoring point in global coordinates system.

In the present invention, global coordinates system is with mean sea level H that height above sea level in geographic area is the geographic area This almost plane of earth ellipsoid surface is xo₁The plane coordinate system of y plane.

In one embodiment, the first geographical coordinate can include with the transfer function of global coordinates：

Wherein, LX_ij、LY_ijRepresent coordinate of the monitoring point in global coordinates system；H represents the mean sea level of geographic area； Lo_ij、La_ijRepresent coordinate of the monitoring point in geographic coordinate system；N_glThe radius of curvature in prime vertical of the earth is represented, anda_glRepresent earth ellipsoid major semiaxis (for example, a_gl=6378.137 × 10³M), W_glIt is intermediate variable, andE represents the first eccentricity of earth ellipsoid,b_glRepresent earth ellipsoid Semi-minor axis (for example, b_gl=6356.7523141 × 10³m)。

In this manner it is possible to the Coordinate Conversion by monitoring point in geographic coordinate system is the coordinate in global coordinates system.

In step S202, mean sea level, the Yi Ji of coordinate, geographic area according to source of leaks in geographic coordinate system Two geographical coordinates and the transfer function of global coordinates, determine coordinate of the source of leaks in global coordinates system.

In one embodiment, the second geographical coordinate can include with the transfer function of global coordinates：

Wherein, LX₀、LY₀Represent coordinate of the source of leaks in global coordinates system；Lo₀、La₀Represent source of leaks in geographical coordinate Coordinate in system.

In this manner it is possible to will leak out Coordinate Conversion of the source in geographic coordinate system is the coordinate in global coordinates system.

In step S203, the seat of coordinate, source of leaks according to monitoring point in global coordinates system in global coordinates system Mark, the mean wind direction angle of geographic area and global coordinates and coordinates of targets transfer function, determine monitoring point in target-based coordinate system In coordinate.

In one embodiment, global coordinates can include with coordinates of targets transfer function：

Wherein, X_ij、Y_ijRepresent coordinate of the monitoring point in target-based coordinate system；θ₀Represent global coordinates systems x-axis it is positive with The angle in due east direction；α represents the mean wind direction angle (that is, the angle of wind direction and due east direction) of geographic area, wherein, α ∈ [0, 2 π), unit is radian (rad).

In one embodiment, the mean wind direction angle of geographic area can be measured by wind direction analyzer.Or, another In one implementation method, the wind direction of the geographic area that can be obtained by meteorological data or weather forecast determines geographic area Mean wind direction angle, wherein, the mean wind direction angle α corresponding to wind direction is as follows：

In one embodiment, θ can be determined in the following manner₀：

First, coordinate of any two monitoring point in collection geographic area in geographic coordinate system, for example, the first monitoring Point (Lo₁₁, La₁₁) and the second monitoring point (Lo₁₂, La₁₂)。

Next, according to equation (1), it can be deduced that coordinate (LX of the two monitoring points in global coordinates system₁₁, LY₁₁) (LX₁₂, LY₁₂)。

Next, building vector v using below equation (4) and equation (5)₁And v₂, wherein, v₁Represent in geographic coordinate system In the second monitoring point to the first monitoring point vector, v₂Represent in global coordinates system the second monitoring point to the first monitoring point to Amount：

v₁=(Lo₁₂-Lo₁₁,La₁₂-La₁₁)(|v₁|≠0) (4)

v₂=(LX₁₂-LX₁₁,LY₁₂-LY₁₁)(|v₂|≠0) (5)

Afterwards, can be according to vector v₁And v₂, using equation (6), determine global coordinates system xo₁The x-axis of y it is positive with just East to angle theta₀(unit is degree, °)：

By xo₁Y coordinate system bypasses origin and perpendicular to xo₁The axle of y is rotated to the positive mean wind direction with monitored area of x-axis Consistent direction, and at origin translation to the vertical projection in ground of source of leaks, obtain new coordinate system xo₂Y, i.e. the coordinates of targets System.

Afterwards, equation (3) is recycled, can be by monitoring point in global coordinates system xo₁Coordinate Conversion in y is to be sat in target Mark system xo₂Coordinate in y.

In this manner it is possible to determine that coordinate of the monitoring point in target-based coordinate system is (X_ij, Y_ij)。

Fig. 3 shows the gas concentration monitoring for finite time steady state leakage of another implementation method of the invention The flow chart of method.As shown in figure 3, the step S103 can include：

In step S301, coordinate, the moment of leakage generation according to source of leaks in geographic coordinate system, leakage occur Date ordinal number and sun altitude calculate function, determine sun altitude.

In one embodiment, sun altitude calculates function and can include：

θ=360d_n/365

h_sun=arcsin [sinLo₀sinσ+cosLo₀cosσcos(15T₀)+La₀-300]

Wherein, h_sunExpression sun altitude (unit is degree, °)；d_nRepresent the date ordinal number (d that leakage occurs_n=0~ 364, wherein, for the d of January 1_n=0, by that analogy)；θ is intermediate variable；σ represents the sun declination in source of leaks location； Lo₀、La₀Represent coordinate of the source of leaks in geographic coordinate system, i.e. the longitude and latitude of source of leaks；T₀Represent leakage occur when Carve (unit is hour), for example, 30 points, T when being directed to moment 15₀=15.5 hours.

In step s 302, at the moment for occurring according to sun altitude, total amount of cloud, low cloud cover and leakage, sun spoke is determined Penetrate grade.

It is possible, firstly, to meteorological condition and table 1 according to geographic area, determine the sky of geographic area total amount of cloud and Low cloud cover.

Table 1

Afterwards, can be according to the sun altitude h determined in step S301_sun, according to table 1 determine total amount of cloud, low clouds At the moment that amount and leakage occur, according to table 2, determine solar radiation grade SL.

Table 2

" night " in table 2 in a few days falls behind the time range of 1 hour to 1 hour before sunrise.

For example, in the case of total amount of cloud≤4 and low cloud cover≤4, if the moment that leakage occurs belongs to night, then Determine solar radiation grade SL=-2.If leaking the moment for occurring is not belonging to night, then if sun altitude h_sun≤ 15 °, it is determined that solar radiation grade SL=-1.

Next, in step S303, wind speed and solar radiation grade according to geographic area determine atmospheric stability.

The wind velocity U (unit is m/s) of geographic area can be obtained for example by anemometry instrument.Obtaining geographic area After wind velocity U, solar radiation grade SL can be combined, atmospheric stability is determined according to table 3, wherein, A, A in table 3~B, C, C ~D, D, E, F represent atmospheric stability.

Table 3

In step s 304, according to atmospheric stability, the floor type of geographic area and monitoring point in target-based coordinate system In coordinate, determine the actual dispersion coefficient at monitoring point.

In one embodiment, the step S304 may comprise steps of：

First, according to atmospheric stability, monitoring point in target-based coordinate system coordinate and basic diffusion coefficient calculates letter Number, determines the basic diffusion coefficient at monitoring point(unit is m, and x/y/z represents x/y/z directions respectively).

In one embodiment, basic diffusion coefficient calculates function and can include：

Wherein,Represent the basic diffusion coefficient at monitoring point；A, A~B, B, B~C, C, C~D, D, D ~E, E, E~F, F represent atmospheric stability；X_ijRepresent x-axis coordinate of the monitoring point in target-based coordinate system.

So, after atmospheric stability is determined, can be determined at monitoring point according to equation (8), (9) and (10) first Basic diffusion coefficient

Next, according to the floor type of geographic area, determining the roughness Z of geographic area_r(unit：m).

In one embodiment, the roughness Z of geographic area can be for example determined according to table 4_r：

Table 4

Floor type
		Grassland, flat open ground	≤0.1
Crop plant area	(0.1,0.3]
		Village, the scattered woods	(0.3,1]
Scattered height building (city)	(1,4)
		Intensive height building (big city)	4

Next, the basic diffusion coefficient at roughness and monitoring point according to geographic area, determines the reality at monitoring point Border diffusion coefficient.

Specifically, in one embodiment, in the case where the roughness satisfaction first of geographic area is pre-conditioned, prison Actual dispersion coefficient at measuring point is the basic diffusion coefficient at monitoring point.For example, described first pre-conditioned is：z_r≤0.1。 That is, meet this first it is pre-conditioned in the case of, actual dispersion coefficient at monitoring point is the basis at monitoring point Shown in diffusion coefficient, such as below equation (11)：

Wherein,The actual dispersion coefficient at monitoring point can be represented.

In another embodiment, in the case where the roughness satisfaction second of geographic area is pre-conditioned, according to big Gas stability, the roughness of geographic area and correction factor calculate function, determine correction factor.For example, described second is default Condition is：z_r＞ 0.1.

In one embodiment, correction factor calculates function and can include：

f_y=1+a₀Z_r (12)

Wherein, f_yRepresent correction factor；Z_rRepresent the roughness of geographic area；a₀Represent pre-determined factor, and a₀With air Stability is corresponding.For example, a can be determined according to atmospheric stability, according to below equation (13)₀Value：

Afterwards, can be according to correction factor f_y, basic diffusion coefficient at monitoring pointAnd basis Diffusion coefficient correction function, draws the actual dispersion coefficient at monitoring point

In one embodiment, basic diffusion coefficient correction function can include：

Thus, it is possible to determine the actual dispersion coefficient at monitoring point(unit is m, x/y/z difference Represent x/y/z directions).

Fig. 4 shows the gas concentration monitoring for finite time steady state leakage of another implementation method of the invention The flow chart of method.As shown in figure 4, the step S104 can include：

In step S401, according to total monitoring duration T_mt(unit：Second or s) and monitoring time interval Δ T (units：Second or S), monitoring cycle sum n is determined.Wherein, T_mt=n Δs T.

Total monitoring duration T_mtCan be pre-arranged with monitoring time interval Δ T.For example, in one embodiment, always Monitoring duration T_mtIt it is 12 hours, Δ T=600s then has T_mt=72 Δ T, n=72.

In step S402, height hm of the monitoring point far from ground is determined_ij.In one embodiment, this highly can be by The spatial dimension for having human living determines, for example, the height of 30 layers of residential building is about 84 meters (by the equal 2.8 meters of calculating of layer).Generally In the case of, hm_ij∈(0,100].It will be understood however, that taking different value also dependent on being actually needed.

In step S403, for each monitoring cycle, according to actual dispersion coefficient, the monitoring point at monitoring point in target The leakage intensity of coordinate, source of leaks in coordinate system, monitoring time interval, the counting of each monitoring cycle, the wind speed of geographic area, Height and first concentration of the effective depth, monitoring point of source of leaks far from ground calculate function, and each monitoring cycle is determined respectively The gas cigarette of internal leakage rolls into a ball the concentration formed at monitoring point, i.e. the gas cigarette of each monitoring cycle internal leakage is rolled into a ball to monitoring point Contribution of concentration.

In one embodiment, the first concentration calculates function and can include：

Wherein,Represent the concentration that the gas cigarette group of k-th monitoring cycle internal leakage is formed at monitoring point, i.e., k-th Contribution of concentration of the gas cigarette group of monitoring cycle internal leakage to monitoring point；K represents that monitoring cycle is counted, and wherein, k=1, 2 ..., n, wherein, n represents monitoring cycle sum；X_ij、Y_ijRepresent coordinate of the monitoring point in target-based coordinate system；U represents geographical The wind speed in region；Represent the actual dispersion coefficient at monitoring point；H₀Represent the effective depth of source of leaks；Q tables Show the leakage intensity of source of leaks；Δ T represents monitoring time interval；hm_ijRepresent height of the monitoring point far from ground.

So, the gas of continuous release can be divided into some discrete cigarette groups, and is calculated at monitoring point by each cigarette The gas leakage concentration of group's contribution.

Additionally, as shown in figure 4, the step S105 can include：In step s 404, according to being let out in each monitoring cycle The gas cigarette of leakage rolls into a ball the concentration and the second concentration calculating function formed at monitoring point, when determining that monitoring point is in total monitoring Gas leakage concentration in length.

In one embodiment, the second concentration calculates function and can include：

Wherein, c_ijRepresent the gas leakage concentration that monitoring point is in total monitoring duration.

The basic thought of this step is, to avoid the amount of calculation in traditional quadrature method with the increase for solving time step number Into the problem that quadratic power increases, the diffusion process of each cigarette group is considered as identical.In this case, can be by k-th Δ T Moment monitoring point (X_ij,Y_ij, hm_ij) concentration at place is considered the concentration that preceding k identical cigarette group, identical diffusion process cause Superposition.

In this manner it is possible to detect gas leakage concentration of the monitoring point in total monitoring duration.Can be by the above method Gas leakage concentration of the monitoring point of any desired on the geographic area of any desired in total monitoring duration is detected, this Sample, can carry out the monitoring of gas leakage diffusion concentration on a large scale.

Additionally, by above-mentioned technical proposal, compared to integration cigarette group model, amount of calculation can be greatly reduced.By inciting somebody to action The Integral Problem of time domain is converted into the superposition problem of spatial concentration, can greatly save the calculating time, accelerates calculating speed, improves Computational efficiency, so as to ensure after generation is leaked, can in time obtain the situation of gas leakage concentration distribution, and then can be have The work such as rescue, evacuation after poisonous gas leakage try to gain time precious to one.

Fig. 5 shows the gas concentration monitoring device for finite time steady state leakage according to the embodiment of the present invention Structured flowchart.As shown in figure 5, the device can include：Acquisition module 501, for gathering on source of leaks and geographic area Coordinate of the monitoring point in geographic coordinate system, and the source of leaks leakage intensity；Coordinate determining module 502, for basis The mean wind direction angle of the coordinate and the geographic area of the source of leaks and the monitoring point in the geographic coordinate system, Determine coordinate of the monitoring point in target-based coordinate system, wherein, the target-based coordinate system be with the source of leaks earthward Vertical subpoint is origin, is the positive coordinate system of x-axis with the mean wind direction of the geographic area；Diffusion coefficient determining module 503, for the coordinate according to the source of leaks in the geographic coordinate system, the meteorologic parameter of the geographic area, leakage hair The floor type of the coordinate and the geographic area of raw time, the monitoring point in the target-based coordinate system, it is determined that described Actual dispersion coefficient at monitoring point；First concentration determination module 504, for the gas of continuous release to be separated into gas cigarette Group, and leakage intensity according to the source of leaks, monitoring time interval, always monitor duration, the actual dispersion at the monitoring point Coordinate in the target-based coordinate system of coefficient, the monitoring point, the wind speed of the geographic area and the source of leaks it is effective Highly, the concentration that each gas cigarette group is formed at the monitoring point is determined；Second concentration determination module 505, for according to institute The concentration that each gas cigarette group is formed at the monitoring point is stated, the leakage gas that the monitoring point is in total monitoring duration is determined Bulk concentration.

Gas concentration monitoring device for finite time steady state leakage of the invention corresponds to is used for finite time stable state The gas concentration monitoring method of leakage, therefore identical content repeats no more.

The preferred embodiment of the present invention is described in detail above in association with accompanying drawing, but, the present invention is not limited to above-mentioned reality The detail in mode is applied, in range of the technology design of the invention, various letters can be carried out to technical scheme Monotropic type, these simple variants belong to protection scope of the present invention.

Additionally, although the operation of the inventive method is described with particular order in the accompanying drawings, this do not require that or Hint must be performed according to the particular order these operation, or the operation having to carry out shown in whole could realize it is desired As a result.Additionally or alternatively, it is convenient to omit some steps, multiple steps are merged into a step to perform, and/or by one Step is decomposed into execution of multiple steps.

One skilled in the art will appreciate that embodiments of the present invention can be implemented as a kind of system, device, equipment, method Or computer program product.Therefore, embodiments of the present invention can be implemented as following form, i.e.,：It is complete hardware, complete Full software (including firmware, resident software, microcode etc.), or the form that hardware and software is combined.

It is further to note that each particular technique feature described in above-mentioned specific embodiment, in not lance In the case of shield, can be combined by any suitable means.In order to avoid unnecessary repetition, the present invention to it is various can The combination of energy is no longer separately illustrated.

Additionally, can also be combined between a variety of implementation methods of the invention, as long as it is without prejudice to originally The thought of invention, it should equally be considered as content disclosed in this invention.

Claims (11)

1. a kind of gas concentration monitoring method for finite time steady state leakage, it is characterised in that the method includes：

Coordinate of the monitoring point in geographic coordinate system on collection source of leaks and geographic area, and the leakage of the source of leaks is strong Degree；

Coordinate and the geographic area according to the source of leaks and the monitoring point in the geographic coordinate system it is average Wind angle, determines coordinate of the monitoring point in target-based coordinate system, wherein, the target-based coordinate system be with the source of leaks to The vertical subpoint on ground is origin, is the positive coordinate system of x-axis with the mean wind direction of the geographic area；

When coordinate, the meteorologic parameter of the geographic area, leakage according to the source of leaks in the geographic coordinate system occur Between, the floor type of coordinate and the geographic area of the monitoring point in the target-based coordinate system, determine the monitoring Actual dispersion coefficient at point；

The gas of continuous release is separated into gas cigarette group, and it is leakage intensity according to the source of leaks, monitoring time interval, total It is coordinate in the target-based coordinate system of monitoring duration, actual dispersion coefficient, the monitoring point at the monitoring point, described The wind speed in region and the effective depth of the source of leaks are managed, determines the gas cigarette group of each monitoring cycle internal leakage in the monitoring The concentration formed at point；

Gas cigarette according to each monitoring cycle internal leakage rolls into a ball the concentration and the second concentration formed at the monitoring point Function is calculated, the gas leakage concentration that the monitoring point is in total monitoring duration is determined, second concentration calculates function bag Include：Wherein, c_ijRepresent the gas leakage concentration that the monitoring point is in total monitoring duration；Table Show the concentration that the gas cigarette group of k-th monitoring cycle internal leakage is formed at the monitoring point；K represents that monitoring cycle is counted, and And, k=1,2 ..., n, wherein, n represents the monitoring cycle sum.

2. method according to claim 1, it is characterised in that it is described according to the source of leaks and the monitoring point described The mean wind direction angle of coordinate and the geographic area in geographic coordinate system, determines the monitoring point in target-based coordinate system Coordinate include：

Coordinate, the mean sea level of the geographic area and the first ground according to the monitoring point in the geographic coordinate system Reason coordinate and the transfer function of global coordinates, determine coordinate of the monitoring point in global coordinates system；

Coordinate, the mean sea level of the geographic area and the second ground according to the source of leaks in the geographic coordinate system Reason coordinate and the transfer function of global coordinates, determine coordinate of the source of leaks in global coordinates system；And

The seat of coordinate, the source of leaks according to the monitoring point in the global coordinates system in the global coordinates system Mark, the mean wind direction angle of the geographic area and global coordinates and coordinates of targets transfer function, determine the monitoring point in institute State the coordinate in target-based coordinate system.

3. method according to claim 2, it is characterised in that the transfer function of first geographical coordinate and global coordinates Including：

$\left\{\begin{array}{c}L{X}_{ij}=(N_{gl}+H)cos(\mathrm{\π}L{o}_{ij}/180)cos(\mathrm{\π}L{a}_{ij}/180)\\ L{Y}_{ij}=(N_{gl}+H)cos(\mathrm{\π}L{o}_{ij}/180)sin(\mathrm{\π}L{a}_{ij}/180)\end{array}\right.$

Wherein, LX_ij、LY_ijRepresent coordinate of the monitoring point in the global coordinates system；H represents the mean sea level； Lo_ij、La_ijRepresent coordinate of the monitoring point in the geographic coordinate system；N_glThe radius of curvature in prime vertical of the earth is represented, and Anda_glRepresent earth ellipsoid major semiaxis, W_glIt is intermediate variable, andE tables Show the first eccentricity of earth ellipsoid,b_glRepresent the semi-minor axis of earth ellipsoid；

Second geographical coordinate includes with the transfer function of global coordinates：

$\left\{\begin{array}{c}L{X}_0=(N_{gl}+H)cos(\mathrm{\π}L{o}_0/180)cos(\mathrm{\π}L{a}_0/180)\\ L{Y}_0=(N_{gl}+H)cos(\mathrm{\π}L{o}_0/180)sin(\mathrm{\π}L{a}_0/180)\end{array}\right.$

Wherein, LX₀、LY₀Represent coordinate of the source of leaks in the global coordinates system；Lo₀、La₀Represent that the source of leaks exists Coordinate in the geographic coordinate system；

The global coordinates include with coordinates of targets transfer function：

$\left\{\begin{array}{c}{X}_{ij}={\mathrm{LX}}_{ij}\cos (\mathrm{\α}-{\mathrm{\θ}}_0)+{\mathrm{LY}}_{ij}\sin (\mathrm{\α}-{\mathrm{\θ}}_0)-\left({\mathrm{LX}}_0\cos \right(\mathrm{\α}-{\mathrm{\θ}}_0)+{\mathrm{LY}}_0\sin (\mathrm{\α}-{\mathrm{\θ}}_0\left)\right)-{\mathrm{LX}}_0\\ Y_{ij}={\mathrm{LY}}_{ij}\cos (\mathrm{\α}-{\mathrm{\θ}}_0)-{\mathrm{LX}}_{ij}\sin (\mathrm{\α}-{\mathrm{\θ}}_0)-\left({\mathrm{LY}}_0\cos \right(-{\mathrm{\θ}}_0)-{\mathrm{LX}}_0\sin (\mathrm{\α}-{\mathrm{\θ}}_0\left)\right)-{\mathrm{LY}}_0\end{array}\right.$

Wherein, X_ij、Y_ijRepresent coordinate of the monitoring point in the target-based coordinate system；θ₀Represent the x of the global coordinates system The positive angle with due east direction of axle；α represents the mean wind direction angle of the geographic area.

4. method according to claim 1, it is characterised in that the meteorologic parameter of the geographic area includes total amount of cloud, low Cloud amount and wind speed；And

Coordinate, the meteorologic parameter of the geographic area, the leakage hair according to the source of leaks in the geographic coordinate system The floor type of the coordinate and the geographic area of raw time, the monitoring point in the target-based coordinate system, it is determined that described Actual dispersion coefficient at monitoring point includes：

Date ordinal number that moment that coordinate according to the source of leaks in the geographic coordinate system, leakage occur, leakage occur, And sun altitude calculates function, determines sun altitude；

According to the moment that the sun altitude, the total amount of cloud, the low cloud cover and the leakage occur, solar radiation is determined Grade；

According to the wind speed and the solar radiation grade, atmospheric stability is determined；

According to the atmospheric stability, the floor type of the geographic area and the monitoring point in the target-based coordinate system In coordinate, determine the actual dispersion coefficient at the monitoring point.

5. method according to claim 4, it is characterised in that the sun altitude calculates function to be included：

θ=360d_n/365

$\mathrm{\σ}=\left[\begin{array}{c}0.006918-0.39912\cos \mathrm{\θ}+0.070257\sin \mathrm{\θ}-0.006758\cos 2\mathrm{\θ}\\ +0.000907\sin 2\mathrm{\θ}-0.002697\cos 3\mathrm{\θ}+0.00148\sin 3\mathrm{\θ}\end{array}\right]180/\mathrm{\π}$

h_sun=arcsin [sinLo₀ sinσ+cosLo₀ cosσcos(15T₀)+La₀-300]

Wherein, h_sunRepresent the sun altitude；d_nRepresent the date ordinal number that the leakage occurs；θ is intermediate variable；σ is represented The sun declination in source of leaks location；Lo₀、La₀Represent coordinate of the source of leaks in the geographic coordinate system；T₀Represent institute State the moment that leakage occurs.

6. method according to claim 4, it is characterised in that described according to the atmospheric stability, the geographic area Coordinate in the target-based coordinate system of floor type and the monitoring point, determine the actual dispersion at the monitoring point Coefficient includes：

Coordinate and basic diffusion coefficient meter according to the atmospheric stability, the monitoring point in the target-based coordinate system Function is calculated, the basic diffusion coefficient at the monitoring point is determined；

According to the floor type of the geographic area, the roughness of the geographic area is determined；

The basic diffusion coefficient at roughness and the monitoring point according to the geographic area, determines the reality at the monitoring point Border diffusion coefficient.

7. method according to claim 6, it is characterised in that the roughness and the prison according to the geographic area Basic diffusion coefficient at measuring point, determines that the actual dispersion coefficient at the monitoring point includes：

In the case where the roughness satisfaction first of the geographic area is pre-conditioned, the actual dispersion coefficient at the monitoring point It is the basic diffusion coefficient at the monitoring point；

The geographic area roughness meet second it is pre-conditioned in the case of, according to the atmospheric stability, described The roughness and correction factor for managing region calculate function, determine correction factor；

According to the basic diffusion coefficient at the correction factor, the monitoring point and basic diffusion coefficient correction function, draw Actual dispersion coefficient at the monitoring point.

8. method according to claim 7, it is characterised in that the basic diffusion coefficient calculates function to be included：

${\mathrm{\σ}}_{ij}^{bx}={\mathrm{\σ}}_{ij}^{by}$

${\mathrm{\σ}}_{ij}^{by}=\left\{\begin{array}{cc}0.22X_{ij}{(1+0.001X_{ij})}^{-1/2}& A\\ \frac{0.22X_{ij}{(1+0.001X_{ij})}^{-1/2}+0.16X_{ij}{(1+0.001X_{ij})}^{-1/2}}{2}& A~B\\ 0.16X_{ij}{(1+0.001X_{ij})}^{-1/2}& B\\ \frac{0.16X_{ij}{(1+0.001X_{ij})}^{-1/2}+0.11X_{ij}{(1+0.001X_{ij})}^{-1/2}}{2}& B~C\\ 0.11X_{ij}{(1+0.001X_{ij})}^{-1/2}& C\\ \frac{0.11X_{ij}{(1+0.001X_{ij})}^{-1/2}+0.08X_{ij}{(1+0.001X_{ij})}^{-1/2}}{2}& C~D\\ 0.08X_{ij}{(1+0.001X_{ij})}^{-1/2},& D\\ \frac{0.08X_{ij}{(1+0.001X_{ij})}^{-1/2}+0.06X_{ij}{(1+0.001X_{ij})}^{-1/2}}{2},& D~E\\ 0.06X_{ij}{(1+0.001X_{ij})}^{-1/2},& E\\ \frac{0.06X_{ij}{(1+0.001X_{ij})}^{-1/2}+0.04X_{ij}{(1+0.001X_{ij})}^{-1/2}}{2},& E~F\\ 0.04X_{ij}{(1+0.001X_{ij})}^{-1/2},& F\end{array}\right.$

${\mathrm{\σ}}_{ij}^{bz}=\left\{\begin{array}{cc}0.20X_{ij},& A\\ \frac{0.20X_{ij}+0.12X_{ij}}{2},& A~B\\ 0.12X_{ij},& B\\ \frac{0.12X_{ij}+0.08X_{ij}{(1+0.002X_{ij})}^{-1/2}}{2},& B~C\\ 0.08X_{ij}{(1+0.002X_{ij})}^{-1/2},& C\\ \frac{0.08X_{ij}{(1+0.002X_{ij})}^{-1/2}+0.06X_{ij}{(1+0.0015X_{ij})}^{-1/2}}{2},& C~D\\ 0.06X_{ij}{(1+0.0015X_{ij})}^{-1/2},& D\\ \frac{0.06X_{ij}{(1+0.0015X_{ij})}^{-1/2}+0.03X_{ij}{(1+0.0003X_{ij})}^{-1/2}}{2},& D~E\\ 0.03X_{ij}{(1+0.0003X_{ij})}^{-1/2},& E\\ \frac{0.03X_{ij}{(1+0.0003X_{ij})}^{-1/2}+0.016X_{ij}{(1+0.0003X_{ij})}^{-1/2}}{2},& E~F\\ 0.016X_{ij}{(1+0.0003X_{ij})}^{-1/2},& F\end{array}\right.$

Wherein,Represent the basic diffusion coefficient at the monitoring point；A, A~B, B, B~C, C, C~D, D, D ~E, E, E~F, F represent the atmospheric stability；X_ijRepresent x-axis coordinate of the monitoring point in the target-based coordinate system；

The correction factor calculates function to be included：

f_y=1+a₀Z_r

Wherein, f_yRepresent the correction factor；Z_rRepresent the roughness of the geographic area；a₀Represent pre-determined factor, and a₀With The atmospheric stability is corresponding；

The basic diffusion coefficient correction function includes：

${\mathrm{\σ}}_{ij}^x={\mathrm{\σ}}_{ij}^y$

$z_j^y={\mathrm{\σ}}_{ij}^{by}*f_y$

${\mathrm{\σ}}_{ij}^z={\mathrm{\σ}}_{ij}^{zy}*f_y$

Wherein,Represent the actual dispersion coefficient at the monitoring point.

9. the method according to any claim in claim 1-8, it is characterised in that the gas by continuous release It is separated into gas cigarette group, and leakage intensity, always monitoring time interval, monitoring duration, the monitoring point according to the source of leaks Coordinate of the actual dispersion coefficient, the monitoring point at place in the target-based coordinate system, the wind speed of the geographic area and described The effective depth of source of leaks, determines the concentration that the gas cigarette group of each monitoring cycle internal leakage is formed at the monitoring point：

According to total monitoring duration and the monitoring time interval, monitoring cycle sum is determined；

Determine height of the monitoring point far from ground；

For each monitoring cycle, according to the actual dispersion coefficient at the monitoring point, the monitoring point in the coordinates of targets The leakage intensity of coordinate, the source of leaks in system, the monitoring time interval, each monitoring cycle are counted, the geographic region The height and the first concentration of the wind speed in domain, the effective depth of the source of leaks, the monitoring point far from ground calculate function, point Do not determine that the gas cigarette of each monitoring cycle internal leakage rolls into a ball the concentration formed at the monitoring point.

10. method according to claim 9, it is characterised in that first concentration calculates function to be included：

$\begin{array}{c}{c}_{ij}^k=\frac{Q\·\mathrm{\Δ}T}{{\left(2\mathrm{\π}\right)}^{3/2}}\exp \left(-\frac{{(X_{ij}-U\·k\·\mathrm{\Δ}T)}^2}{2{\left({\mathrm{\σ}}_{ij}^x\right)}^2}\right)\exp \left(-\frac{{Y_{ij}}^2}{2{\left({\mathrm{\σ}}_{ij}^y\right)}^2}\right)\\ \·\left(\exp \left(-\frac{{({\mathrm{hm}}_{ij}-H_0)}^2}{2{\left({\mathrm{\σ}}_{ij}^z\right)}^2}\right)+\exp \left(-\frac{{\left({\mathrm{hm}}_{ij}+H_0\right)}^2}{2{\left({\mathrm{\σ}}_{ij}^z\right)}^2}\right)\right)\end{array}$

Wherein, X_ij、Y_ijRepresent coordinate of the monitoring point in target-based coordinate system；U represents the wind speed of the geographic area；Represent the actual dispersion coefficient at the monitoring point；H₀Represent the effective depth of the source of leaks；Q represents institute State the leakage intensity of source of leaks；Δ T represents the monitoring time interval；hm_ijRepresent height of the monitoring point far from ground.

11. a kind of gas concentration monitoring devices for finite time steady state leakage, it is characterised in that the device includes：

Acquisition module is for gathering coordinate of the monitoring point on source of leaks and geographic area in geographic coordinate system and described The leakage intensity of source of leaks；

Coordinate determining module, for the coordinate according to the source of leaks and the monitoring point in the geographic coordinate system and The mean wind direction angle of the geographic area, determines coordinate of the monitoring point in target-based coordinate system, wherein, the coordinates of targets It, with the source of leaks vertical subpoint earthward as origin, is positive x-axis with the mean wind direction of the geographic area that system is Coordinate system；

Diffusion coefficient determining module, for the coordinate according to the source of leaks in the geographic coordinate system, the geographic area Meteorologic parameter, leakage coordinate and the geographic area in the target-based coordinate system of time of origin, the monitoring point Floor type, determines the actual dispersion coefficient at the monitoring point；

First concentration determination module, for the gas of continuous release to be separated into gas cigarette group, and letting out according to the source of leaks Leakage intensity, monitoring time interval, always monitoring duration, actual dispersion coefficient, the monitoring point at the monitoring point be in the mesh The effective depth of coordinate, the wind speed of the geographic area and the source of leaks in mark coordinate system, determines in each monitoring cycle The gas cigarette of leakage rolls into a ball the concentration formed at the monitoring point；

Second concentration determination module, for rolling into a ball the shape at the monitoring point according to the gas cigarette of each monitoring cycle internal leakage Into concentration and the second concentration calculate function, determine the gas leakage concentration that the monitoring point is in total monitoring duration, institute Stating the second concentration calculating function includes：Wherein, c_ijRepresent that the monitoring point is in total monitoring duration Interior gas leakage concentration；Represent the concentration that the gas cigarette group of k-th monitoring cycle internal leakage is formed at the monitoring point； K represents that monitoring cycle is counted, also, k=1,2 ..., n, wherein, n represents the monitoring cycle sum.