Method for determining the number of gas leak sources

Abstract

The application provides a method for determining the number of gas leakage sources, comprising the following steps: (A1) monitoring the first movement of a vehicle on the same path, recording the peak-to-peak value of the gas concentration, the position and the wind direction angle, (A2) performing Gaussian fitting on the concentration peaks to obtain the full width at half maximum corresponding to each concentration peak; obtaining the effective concentration peak area of each concentration peak; (A3) counting the distance between different concentration peaks, establishing a concentration peak distance set, and establishing an effective concentration peak area set; (A4) measuring multiple times in the manner of steps (A1)-(A3), respectively obtaining the concentration peak distance set; establishing the effective concentration peak area set; (A5) extracting one element from each set {S 1i} and {S ji}, each element corresponding to a wind direction angle, forming K1 elements; fitting the relationship between the extracted elements and the wind direction angle to determine whether there are two leakage sources. The application has the advantages of accurate results.

Description

Technical Field

[0001] This invention relates to gas analysis, and more particularly to a method for determining the number of gas leak sources. Background Technology

[0002] Natural gas is an essential energy source in our daily lives; however, frequent gas leaks pose significant safety hazards. my country has 110,000 kilometers of natural gas pipelines, most of which are buried underground. Daily monitoring of gas pipeline leaks presents a large workload and high level of difficulty. The deployment of mobile leak gas monitoring equipment can significantly reduce the workload of inspection personnel, and through rapid leak location and accurate concentration measurement, efficiently eliminate safety hazards and protect the lives and property of the public.

[0003] Mobile monitoring of leaked gas can combine information from multiple measurement sessions within a region using wind direction, wind speed, atmospheric stability, gas concentration, and vehicle location. This allows for high-speed, large-scale collection and monitoring of plumes of leaked gas in the air. Ultimately, through statistical analysis of the measurement data, the source of the leaked gas can be traced, and the emission rate of the leaking source can be estimated. In practice, mobile monitoring of leaked gas often measures multiple consecutive concentration peaks. Current technologies often treat these as being caused by a single leaking source in turbulent gas, neglecting the possibility of multiple leaking sources causing a multi-peak distribution. Therefore, traditional algorithms cannot estimate the number of leaking sources in a multi-peak distribution of leaked gas concentration. Summary of the Invention

[0004] To address the shortcomings of the existing technical solutions, the present invention provides a method for determining the number of gas leakage sources.

[0005] The objective of this invention is achieved through the following technical solution:

[0006] A method for determining the number of gas leak sources, comprising the following steps:

[0007] (A1) Monitor the vehicle's first movement along the same path and record the corresponding peak-to-peak values ​​of the leaked gas concentration C. 1i Position x 1i And the wind direction angle θ1, i = 1, 2, ..., N1, where N1 is a positive integer;

[0008] (A2) Perform Gaussian fitting on the concentration peaks to obtain the full width at half maximum (FWHM) W corresponding to each concentration peak. 1i i = 1, 2, ..., N1, where N1 is a positive integer;

[0009] The effective concentration peak area of ​​each concentration peak was obtained. N1 is a positive integer;

[0010] (A3) Statistical analysis of the spacing L between peaks of different concentrations 1mn =x 1m -x 1n If m < n ≤ N1, establish a set of concentration peak spacings {L}. 1mn}, this set consists of Composed of the spacing between concentration peaks;

[0011] Establish the effective concentration peak area set {S} 1i The set consists of N1 valid areas, i = 1, 2, ..., N1, where N1 is a positive integer;

[0012] (A4) Perform multiple measurements following steps (A1)-(A3) to obtain the concentration peak spacing set {L}. jmn}, this set consists of The concentration peaks are composed of inter-peak spacing, j = 2, 3, ..., K1, where K1 is a positive integer;

[0013] Establish the effective concentration peak area set {S} ji}, the set consists of N j The effective areas are composed of j = 2, 3, ..., K1, where K1 is a positive integer;

[0014] (A5) From each set {S 1i}、{S ji Each element is extracted from the group, and each element corresponds to a wind direction angle, forming K1 elements;

[0015] Fit the relationship between extracted elements and wind direction angle y is the extracted element, θ is the wind direction angle, and d 12 α and α are the fitting parameters, respectively;

[0016] If the fitted correlation coefficient is higher than the first threshold, calculate the relative standard deviation of the effective concentration peak area of ​​the two possible leakage sources. If the relative standard deviation is lower than the second threshold, then it is determined that there are two leakage sources.

[0017] If the fitted correlation coefficient is lower than the first threshold or the relative standard deviation is higher than the second threshold, it indicates that there are no two leakage sources.

[0018] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0019] 1. Traditional algorithms cannot estimate the number of leakage sources in a multi-peak distribution of leaked gas concentration. This invention utilizes the relationship between the concentration peak spacing and the wind direction angle θ to estimate the number of leakage sources in a multi-peak distribution of leaked gas concentration.

[0020] 2. This invention proposes the concept of effective concentration peak area, distinguishing between multiple peaks caused by turbulence and multiple peaks caused by multiple leakage sources. Attached Figure Description

[0021] The disclosure of this invention will become more readily understood with reference to the accompanying drawings. It will be readily understood by those skilled in the art that these drawings are merely illustrative of the technical solutions of this invention and are not intended to limit the scope of protection of this invention. In the drawings:

[0022] Figure 1 This is a schematic diagram showing the correspondence between concentration peaks and positions according to an embodiment of the present invention;

[0023] Figure 2 This is a schematic diagram of the parameter relationship between two leakage sources according to an embodiment of the present invention. Detailed Implementation

[0024] Figure 1-2 The following description illustrates optional embodiments of the invention to teach those skilled in the art how to implement and reproduce the invention. Some conventional aspects have been simplified or omitted to explain the technical solutions of the invention. Those skilled in the art should understand that variations or substitutions derived from these embodiments will be within the scope of the invention. Those skilled in the art should understand that the following features can be combined in various ways to form multiple variations of the invention. Therefore, the invention is not limited to the following optional embodiments, but is defined only by the claims and their equivalents.

[0025] Example 1:

[0026] The method for determining the number of gas leakage sources according to an embodiment of the present invention includes the following steps:

[0027] (A1) Monitor the vehicle's first movement along the same path and record the corresponding peak-to-peak values ​​of the leaked gas concentration C. 1i Position x 1i And the wind direction angle θ1, i = 1, 2...N1, where N1 is a positive integer; such as Figure 1 As shown,

[0028] (A2) Perform Gaussian fitting on the concentration peaks to obtain the full width at half maximum (FWHM) W corresponding to each concentration peak. 1i i = 1, 2, ..., N1, where N1 is a positive integer;

[0029] The effective concentration peak area of ​​each concentration peak was obtained. N1 is a positive integer;

[0030] (A3) Statistical analysis of the spacing L between peaks of different concentrations 1mn =x 1m -x 1n If m < n ≤ N1, establish a set of concentration peak spacings {L}. 1mn}, this set consists of Composed of the spacing between concentration peaks;

[0031] Establish the effective concentration peak area set {S} 1i The set consists of N1 valid areas, i = 1, 2, ..., N1, where N1 is a positive integer;

[0032] (A4) Perform multiple measurements following steps (A1)-(A3) to obtain the concentration peak spacing set {L}. jmn}, this set consists of The concentration peaks are composed of inter-peak spacing, j = 2, 3, ..., K1, where K1 is a positive integer;

[0033] Establish the effective concentration peak area set {S} ji}, the set consists of N j The effective areas are composed of j = 2, 3, ..., K1, where K1 is a positive integer;

[0034] (A5) From each set {S 1i}、{S ji Each element is extracted from the group, and each element corresponds to a wind direction angle, forming K1 elements;

[0035] Fit the relationship between extracted elements and wind direction angle y is the extracted element, θ is the wind direction angle, and d 12 α and α are the fitting parameters, such as Figure 2 As shown;

[0036] If the fitted correlation coefficient is higher than the first threshold of 0.9, calculate the relative standard deviation of the effective concentration peak area of ​​the two possible leakage sources. If the relative standard deviation is lower than the second threshold of 0.2, then it is determined that there are two leakage sources.

[0037] If the fitted correlation coefficient is lower than the first threshold or the relative standard deviation is higher than the second threshold, it indicates that there are no two leakage sources.

[0038] (A6) Repeat step (A5), each time extracting a different element from the set, to obtain the number of leakage sources on the path.

[0039] Example 2:

[0040] An application example of the method for determining the number of gas leakage sources according to Embodiment 1 of the present invention.

[0041] In this application example, the method for determining the number of gas leak sources includes the following steps:

[0042] (A1) Monitor the vehicle's first movement along the same path and record the corresponding peak-to-peak values ​​of the leaked gas concentration C. 1i Position x 1i And the wind direction angle θ1, i = 1, 2, ..., N1, where N1 is a positive integer;

[0043] (A2) Perform Gaussian fitting on the concentration peaks to obtain the full width at half maximum (FWHM) W corresponding to each concentration peak. 1i i = 1, 2, ..., N1, where N1 is a positive integer;

[0044] The effective concentration peak area of ​​each concentration peak was obtained. N1 is a positive integer;

[0045] (A3) Statistical analysis of the spacing L between peaks of different concentrations 1mn =x 1m -x 1n If m < n ≤ N1, establish a set of concentration peak spacings {L}. 1mn}, this set consists of Composed of the spacing between concentration peaks;

[0046] Establish the effective concentration peak area set {S} 1i The set consists of N1 valid areas, i = 1, 2, ..., N1, where N1 is a positive integer;

[0047] (A4) Perform 10 measurements following steps (A1)-(A3) to obtain the concentration peak spacing set {L}. jmn}, this set consists of It consists of concentration peak spacing, j = 2, 3, ..., 10;

[0048] Establish the effective concentration peak area set {S} ji}, the set consists of N j Composed of several effective areas, j = 2, 3, ..., 10;

[0049] (A5) From each set {S 1i}、{S ji One element is extracted from each of the following groups, and each element corresponds to a wind direction angle, forming 10 elements;

[0050] Fit the relationship between extracted elements and wind direction angle y is the extracted element, θ is the wind direction angle, and d 12 α and α are the fitting parameters, respectively;

[0051] If the fitted correlation coefficient is higher than the first threshold, calculate the relative standard deviation of the effective concentration peak area of ​​the two possible leakage sources. If the relative standard deviation is lower than the second threshold, then it is determined that there are two leakage sources.

[0052] If the fitted correlation coefficient is lower than the first threshold or the relative standard deviation is higher than the second threshold, it indicates that there are no two leakage sources.

[0053] (A6) Repeat step (A5), each time extracting a different element from the set, to obtain the number of leakage sources on the path.

Claims

1. A method for determining the number of gas leak sources, the method comprising the following steps: (A1) When the monitoring vehicle moves along the same path for the first time, record the corresponding peak-to-peak values ​​of the leaked gas concentration C. 1i Position x 1i And the wind direction angle θ1, i=1,2···N1, where N1 is a positive integer; (A2) Perform Gaussian fitting on the concentration peaks to obtain the full width at half maximum (WW) corresponding to each concentration peak. 1i i = 1, 2, ..., N1, where N1 is a positive integer; The effective concentration peak area S of each concentration peak was obtained. 1i =C 1i ·W 1i / sinθ1, i=1,2···N1, N1 is a positive integer; (A3) Statistically analyze the spacing L between peaks at different concentrations 1mn = x 1m - x 1n , where m < n ≤ N1, establish a set of peak spacing of concentrations {L 1mn}, and this set consists of N1(N1 - 1) / 2 peak spacings of concentrations; Establish the effective concentration peak area set {S} 1i The set consists of N1 valid areas, i = 1, 2, ..., N1, where N1 is a positive integer; (A4) Perform multiple measurements following steps (A1)-(A3) to obtain the concentration peak spacing set {L}. jmn The set consists of N j (N j The concentration peaks consist of -1) / 2 inter-peak spacings, j=2,3···K1, where K1 is a positive integer; Establish the effective concentration peak area set {S} ji }, the set consists of N j The effective areas are composed of j = 2, 3, ..., K1, where K1 is a positive integer; (A5) From each set {S 1i }、{S ji Each element is extracted from the group, and each element corresponds to a wind direction angle, forming K1 elements; The relationship between the extracted element S and the wind direction angle θ is S=d 12 ·sin[abs(θ-α)] / sinθ,d 12 It is the straight-line distance between the two leakage sources fitted, and α is the angle between the straight line connecting the two leakage sources and the direction of vehicle movement; If the fitted correlation coefficient is higher than the first threshold, calculate the relative standard deviation of the effective concentration peak area of ​​the two possible leakage sources. If the relative standard deviation is lower than the second threshold, then it is determined that there are two leakage sources. If the fitted correlation coefficient is lower than the first threshold or the relative standard deviation is higher than the second threshold, it indicates that there are no two leakage sources.

2. The method for determining the number of gas leakage sources according to claim 1, characterized in that, The method for determining the number of gas leakage sources also includes the following steps: (A6) Repeat step (A5), each time extracting a different element from the set, to obtain the number of leakage sources on the path.

Images