CN111256043A - Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array - Google Patents

Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array Download PDF

Info

Publication number
CN111256043A
CN111256043A CN202010193729.0A CN202010193729A CN111256043A CN 111256043 A CN111256043 A CN 111256043A CN 202010193729 A CN202010193729 A CN 202010193729A CN 111256043 A CN111256043 A CN 111256043A
Authority
CN
China
Prior art keywords
leakage source
source
leakage
potential
concentration
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
Application number
CN202010193729.0A
Other languages
Chinese (zh)
Inventor
周康
李斐
庄俊亿
蔡浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Tech University
Original Assignee
Nanjing Tech University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Tech University filed Critical Nanjing Tech University
Priority to CN202010193729.0A priority Critical patent/CN111256043A/en
Publication of CN111256043A publication Critical patent/CN111256043A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17D—PIPE-LINE SYSTEMS; PIPE-LINES
    • F17D5/00—Protection or supervision of installations
    • F17D5/02—Preventing, monitoring, or locating loss
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004—Gaseous mixtures, e.g. polluted air
    • G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
    • G01N33/0062—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method, e.g. intermittent, or the display, e.g. digital
    • G01N33/0063—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method, e.g. intermittent, or the display, e.g. digital using a threshold to release an alarm or displaying means
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/0004—Gaseous mixtures, e.g. polluted air
    • G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
    • G01N33/0062—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method, e.g. intermittent, or the display, e.g. digital
    • G01N33/0067—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method, e.g. intermittent, or the display, e.g. digital by measuring the rate of variation of the concentration

Abstract

The invention discloses a method for rapidly identifying the position and release amount of a leakage source when a gas pipeline of a comprehensive pipe rack leaks by using an impulse response and concentration sensor array. According to the method, firstly, pulse source release is carried out at the position of a potential leakage source to obtain a corresponding response matrix, then, the strength and the position of the leakage source can be calculated according to the real-time monitoring feedback numerical value of a concentration sensor in a comprehensive pipe gallery cabin by combining a regularization method and a Bayesian probability model. The method can be used for quickly arranging and accurately performing inversion estimation on the position of the leakage source of the gas pipeline of the comprehensive pipe gallery, has great significance to safety engineering, and can be expanded to the identification of pollution sources in indoor space.

Description

Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array
(I) technical field
The invention belongs to the field of urban public safety, and particularly relates to a method for identifying a gas leakage source of a gas pipeline in a comprehensive pipe gallery cabin.
(II) background of the invention
City utility tunnel is various pipelines such as catchment, electricity, communication, gas centralized pipeline laying mode in an organic whole, compares with the direct buried pipe gallery of tradition, and the underground space who occupies is few, has avoided the construction of digging repeatedly on ground, plays fine guarantee effect to urban traffic and city development. But the utility tunnel belongs to underground narrow and small space, and the gas pipeline is gone into the corridor and has also caused great hidden danger to utility tunnel's security simultaneously, and the danger that the gas leakage explosion caused is also huge, has also constituted the threat to the safe operation in city and resident's life safety. At present, the emergency measure when the gas pipeline of the comprehensive pipe gallery leaks is mainly to perform gas alarm, and when the gas concentration reaches the explosion limit, a gas pipeline sectional valve is cut off in an emergency or a ventilation device is started in time. Therefore, if the intensity and the position of a leakage source are identified when gas leaks, the source is maintained and treated, and the method plays an important role in reducing secondary accidents and accelerating maintenance and recovery.
Disclosure of the invention
Solves the technical problem
Aiming at the defects of the existing pipe gallery gas leakage monitoring and preventing method, the invention provides a comprehensive pipe gallery gas leakage source identification method based on an impulse response method and a sensor array, and the intensity and the position of the leakage source can be quickly and accurately identified
Technical scheme
In order to achieve the purpose, the invention is realized by the following technical scheme:
(1) firstly, determining the position of each concentration sensor in the pipe gallery and the position of a potential leakage source;
(2) a quantity of gas is released in pulses at the location of the potential leakage source, respectively, a response matrix is obtained for the different concentration sensors, which represents the relationship between the leakage source intensity q and the sensor concentration C, with a linear relationship C-a · q,
Figure BSA0000204120330000011
in the formula
Figure BSA0000204120330000014
Is shown at tnMeasuring the value of a concentration sensor at a certain position at a certain moment;
Figure BSA0000204120330000012
is shown at tnA response factor at a time;
Figure BSA0000204120330000013
is shown at tnLeakage intensity of leakage source at a time;
(3) after a response matrix of a potential leakage source is obtained, combining measured values of part of concentration sensors, and calculating the strength of the leakage source reversely by using a regularization method;
(4) obtaining posterior probability of each potential leakage source based on Bayesian criterion by using monitoring data of other concentration sensors
Figure BSA0000204120330000021
Wherein N is the number of potential leakage sources, k represents the number of potential leakage sources, and p (Y)k) To a priori probability, L (O | Y)k) In order to be a function of the likelihood,
Figure BSA0000204120330000022
wherein C isMWhich is indicative of the measured value of the concentration sensor,
Figure BSA0000204120330000023
representing the concentration sensor value calculated by the response matrix for the kth potential leakage source, wherein sigma is the data standard deviation;
(5) finally using the posterior probability p (Y)kAnd | O) calculating the probability of each potential leakage source, wherein the potential leakage source position with the maximum probability value is the real leakage source position. In the step (2), the pulse release time is set as the calculated time step, and the amount of the gas released by the pulse reaches the lower measurement limit of the concentration sensor.
In the step (4), the rest concentration sensors are different from the sensors used for back-calculating the intensity q of the leakage source.
In the step (4), the prior probability p (Y)k) The probability of leakage sources at different potential locations can be specified to be equal.
Advantageous effects
According to the invention, the pulse response method is used for identifying the intensity of the leakage source and the position of the leakage source in the pipe gallery, the defect that the intensity of the leakage source and the accurate position of the leakage source cannot be identified by cutting off the fireproof subarea when the concentration sensor reaches the alarm concentration in the traditional gas leakage monitoring and controlling method is solved, and the pulse response method avoids the step of CFD numerical simulation in the traditional pollution source identification method, so that the method is more rapid and convenient. The method can effectively, conveniently and quickly identify the intensity and the position of the leakage source in the underground comprehensive pipe gallery, and realize quick identification and early warning.
(IV) description of the drawings
FIG. 1 is a schematic diagram of the arrangement of the release position and concentration sensors of the experimental model constructed by the invention
FIG. 2 is a frame diagram of a technical solution for identifying the location and strength of a leakage source in a utility tunnel by using an impulse response method according to the present invention
FIG. 3 is a graph showing the result of inverse calculation of leakage source intensity in a utility tunnel by using an impulse response method
(V) detailed description of the preferred embodiments
The following detailed description of specific embodiments of the present invention is provided in connection with the accompanying drawings and examples, which are set forth to illustrate, but are not to be construed to limit the scope of the invention. After reading the teaching of the present invention, those skilled in the art can make various changes or modifications to the invention, and these equivalents also fall within the scope of the claims appended to the present application.
As shown in fig. 1, an underground comprehensive pipe gallery model with the length of 9 meters, the width of 1.7 meters and the height of 2.1 meters is established, and a gas pipeline with the length of 9 meters traverses the model. The position of a real leakage source (potential source 1) is arranged below a gas pipeline and at a position 1.63 meters away from an air inlet, and two potential sources 2 and 3 are arranged. The concentration sensors are arranged on the section of the 'S' 9 meters away from the air inlet, and the concentration sensors are four concentration sensors S1, S2, S3 and S4.
In an experiment, at the position of a potential leakage source 1, tracer gas (carbon dioxide) is released for 10 seconds (in a pulse form) at a flow rate of 45L/min to obtain data of the four concentration sensors changing along with time, and the data is led into MATLAB for data processing to obtain a response matrix A corresponding to different concentration sensors at the position of the potential leakage source 1; and then constantly releasing the tracer gas at the real leakage source 1 with the leakage intensity of 30L/min, and finally determining the intensity q of the leakage source corresponding to different sensors by combining the concentration value C fed back by the position of each sensor and the calculated concentration response matrix A.
Since the density response matrix a is a pathologic matrix, that is, q cannot be obtained by directly inverting C ═ a · q, the stability of matrix operation is enhanced by the Tikhonov regularization method. The Tikhonov regularization method is a problem of converting the relation between the intensity of a leakage source and the concentration, C ═ a · q, into the minimum objective function of the following formula:
Figure BSA0000204120330000031
in the formula (I), the compound is shown in the specification,
Figure BSA0000204120330000032
is the matrix two norm: l is a regularization matrix and λ is a regularization parameter. Derivative of z (q) determines the leakage source intensity q that causes the above equation to take a local minimum. From the above equation, if the concentration time series C and the response matrix a of the concentration sensor are known, the magnitude q of the leakage source intensity of the corresponding sensor can be determined, and the result of the leakage source intensity is calculated back, as shown in fig. 3. Wherein the leakage intensity of the real source, i.e. 30L/min; s1, S2, S3 and S4 are respectively the back calculation results of different sensor positions. The results of the calculations show that each sensor position is able to substantially back-calculate the intensity of the leak source, but is generally less than the actual leak intensity. Meanwhile, the feasibility of reverse identification under the condition of actual leakage by using the method for acquiring the response matrix by using the impulse response method is verified.
A Bayesian probability model is adopted for determining the position of the leakage source, a concentration sensor S1 is taken as an example, an S1 concentration sensor is selected for identifying the intensity of the leakage source, and an S2 concentration sensor is used for identifying the position of the leakage source.
The leakage intensities of the three potential sources are first back-calculated using the data from the S1 concentration sensor and the response matrices obtained for the leakage at the three potential source locations, respectively.
Then, using three potential sourcesThe leakage intensity and the corresponding response matrix of the S2 concentration sensor calculate the likelihood probability L (O | Y) of the potential sources 1, 2 and 31)、L(O|Y2)、L(O|Y3)。
Finally, the posterior probabilities p (Y) of the three potential leakage sources are calculated1|O)，p(Y2|O)，p(Y3And | O), comparing the three posterior probability values, and finding that the posterior probability value of the potential leakage source 1 is the maximum value, namely the predicted real source position. In conclusion, the identification of the release strength and the position of the leakage source is completed.

Claims (4)

1. Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array, its characterized in that includes the following steps:
(1) firstly, determining the position of each concentration sensor in the pipe gallery and the position of a potential leakage source;
(2) a quantity of gas is released in pulses at the location of the potential leakage source, respectively, a response matrix is obtained for the different concentration sensors, which represents the relationship between the leakage source intensity q and the sensor concentration C, with a linear relationship C-a · q,
Figure FSA0000204120320000011
in the formula
Figure FSA0000204120320000012
Is shown at tnMeasuring the value of a concentration sensor at a certain position at a certain moment;
Figure FSA0000204120320000013
is shown at tnA response factor at a time;
Figure FSA0000204120320000014
is shown at tnLeakage strength of leakage source at time instant.
(3) And after the response matrix of the potential leakage source is obtained, combining the measured values of part of the concentration sensors, and calculating the strength of the leakage source reversely by using a regularization method.
(4) Obtaining posterior probability p (Y) of each potential leakage source based on Bayesian criterion by using monitoring data of other concentration sensorsk|O)，
Figure FSA0000204120320000015
Wherein N is the number of potential leakage sources, k represents the number of potential leakage sources, and p (Y)k) To a priori probability, L (O | Y)k) In order to be a function of the likelihood,
Figure FSA0000204120320000016
wherein C isMWhich is indicative of the measured value of the concentration sensor,
Figure FSA0000204120320000017
and (4) representing the sensor value calculated by the response matrix of the kth potential leakage source, wherein the sigma is the standard deviation of the data.
(5) Finally using the posterior probability p (Y)kAnd | O) calculating the probability of each potential leakage source, wherein the potential leakage source position with the maximum probability value is the real leakage source position.
2. The method of claim 1, wherein in step (2), the pulse release time is set to a calculated time step, and the amount of gas released by the pulse is to reach a lower measurement limit of the concentration sensor.
3. The method of claim 1, wherein in step (4), the remaining concentration sensors are sensors other than the sensor used to back-calculate the intensity q of the leak source.
4. The method as claimed in claim 1, wherein in step (4), the prior probability p (Y) isk) The probability of leakage sources at different potential locations can be specified to be equal.
CN202010193729.0A 2020-03-18 2020-03-18 Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array Pending CN111256043A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010193729.0A CN111256043A (en) 2020-03-18 2020-03-18 Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010193729.0A CN111256043A (en) 2020-03-18 2020-03-18 Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array

Publications (1)

Publication Number Publication Date
CN111256043A true CN111256043A (en) 2020-06-09

Family

ID=70951769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010193729.0A Pending CN111256043A (en) 2020-03-18 2020-03-18 Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array

Country Status (1)

Country Link
CN (1) CN111256043A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113586968A (en) * 2021-06-01 2021-11-02 北京市燃气集团有限责任公司 Natural gas leakage source positioning method and device
CN113808372A (en) * 2021-09-26 2021-12-17 深圳市智联云控科技有限公司 Combustible gas safety monitoring method and device and electronic equipment

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551064A (en) * 2009-05-22 2009-10-07 重庆大学 Water supply pipe leakage detection locating signal processing method
CN107145737A (en) * 2017-04-28 2017-09-08 南京工业大学 One kind is based on the markovian unsteady state flow reverse recognizer of pollution sources off field
EP3327206A1 (en) * 2016-11-25 2018-05-30 Tata Consultancy Services Limited Ranking pipes for maintenance in pipe networks using approximate hydraulic metrics
CN108416478A (en) * 2018-03-19 2018-08-17 国通广达（北京）技术有限公司 Gas concentration LSTM prediction techniques and system inside a kind of pipe gallery
CN111257519A (en) * 2020-03-18 2020-06-09 南京工业大学 Leakage source identification device and method in building space based on sensor array and impulse response principle

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101551064A (en) * 2009-05-22 2009-10-07 重庆大学 Water supply pipe leakage detection locating signal processing method
EP3327206A1 (en) * 2016-11-25 2018-05-30 Tata Consultancy Services Limited Ranking pipes for maintenance in pipe networks using approximate hydraulic metrics
CN107145737A (en) * 2017-04-28 2017-09-08 南京工业大学 One kind is based on the markovian unsteady state flow reverse recognizer of pollution sources off field
CN108416478A (en) * 2018-03-19 2018-08-17 国通广达（北京）技术有限公司 Gas concentration LSTM prediction techniques and system inside a kind of pipe gallery
CN111257519A (en) * 2020-03-18 2020-06-09 南京工业大学 Leakage source identification device and method in building space based on sensor array and impulse response principle

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XIAORAN LIU等: "Dynamical source term estimation in a multi-compartment building under", 《BUILDING AND ENVIRONMENT》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113586968A (en) * 2021-06-01 2021-11-02 北京市燃气集团有限责任公司 Natural gas leakage source positioning method and device
CN113808372A (en) * 2021-09-26 2021-12-17 深圳市智联云控科技有限公司 Combustible gas safety monitoring method and device and electronic equipment

Similar Documents

Publication Publication Date Title
Pospelov et al. Short-term fire forecast based on air state gain recurrence and zero-order brown model
CN111256043A (en) Comprehensive pipe gallery gas leakage source identification method based on pulse response method and sensor array
CN110851981B (en) Method for realizing rapid tracing of sudden water pollution
CN102072409B (en) Pipe network leakage monitoring method combining leakage probability calculation and recorder monitoring
Wang et al. Numerical study on the diffusion process of pinhole leakage of natural gas from underground pipelines to the soil
CN205748814U (en) Valves leakage monitoring system
CN103147799B (en) Gas discharge characteristic-based outburst danger early warning method
CN103062628A (en) Urban gas buried pipeline leakage detection and locating method and applications thereof
CN106065773A (en) A kind of gas well liquid loading tubing string leakage detection method based on gas tracer
CN111257519A (en) Leakage source identification device and method in building space based on sensor array and impulse response principle
Wang et al. Study and application of a new gas pressure inversion model in coal seam while drilling based on directional drilling technology
CN108071941B (en) A kind of booster identification and localization method based on time series vector similarity analysis
Lei et al. A method for identifying the fire status through ventilation systems using tracer gas for improved rescue effectiveness in roadway drivage of coal mines
CA3062555C (en) Porous medium permeability prediction method based on machine image intelligent learning
Hu et al. Inverse tracking of an airborne pollutant source location in a residential apartment by joint simulation of CFD and a multizone model
Liu et al. City pipe network intelligent service based on GIS and internet of things
Piazza et al. Optimal positioning of water quality sensors in water distribution networks: Comparison of numerical and experimental results
Zhou et al. Risk index assessment for urban natural gas pipeline leakage based on artificial neural network
Neumann et al. Tomographic reconstruction of soil gas distribution from multiple gas sources based on sparse sampling
CN112664839B (en) Method and system for predicting and tracing combustible gas diffusion of communication pipeline
Wu et al. BI-IEnKF coupling model for effective source term estimation of natural gas leakage in urban utility tunnels
CN115127036A (en) Municipal gas pipe network leakage positioning method and system
CN111442194B (en) Method for measuring pipeline leakage by using unmanned aerial vehicle carrying hygrometer
Wang et al. Leakage risk quantitative calculation model and its application for anaerobic reactor
CN106250578A (en) Coal mine gas detection method

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200609