Pipeline leakage positioning method and device, chip and terminal
By installing sensors on chemical plant pipelines to acquire acoustic signals, and using cross-correlation technology and pipeline attribute parameters for frequency band segmentation and weighted averaging, the problem of large positioning errors in chemical plant pipeline leak location was solved, achieving high-precision leak location.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- SEVNCE ROBOTICS CO LTD
- Filing Date
- 2024-01-11
- Publication Date
- 2026-04-14
AI Technical Summary
Existing acoustic detection and location methods cannot effectively improve the location accuracy in pipeline leak location in chemical plants, especially in complex and variable non-stationary noise environments, where the location error is large and it is difficult to meet engineering requirements.
The first and second sensors are used to acquire the sound signal of pipeline leakage. The signal is divided into equal-width frequency bands by cross-correlation technology. The leakage location is located by combining the target pipeline's attribute parameters such as average sound velocity and weighting coefficient. The leakage signal is adaptively separated by using linear relationship and weighting coefficient for weighted averaging.
It effectively solves the error problem of pipeline leak location under high background noise, improves the location accuracy, and adapts to the leak location needs in the complex environment of chemical plants.
Smart Images

Figure CN117869809B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of fault diagnosis technology, and in particular to a method, device, chip, and terminal for locating pipeline leaks. Background Technology
[0002] Pipeline networks in chemical plants are prone to creep cracks and corrosion aging due to long-term exposure to high temperature, high pressure and highly corrosive environments, which can lead to pipeline leaks and subsequently cause safety accidents such as fires, explosions and poisoning.
[0003] Acoustic detection and localization methods are widely used in pipeline leak location due to their advantages such as high detection sensitivity, small location error, fast response speed, and low cost. However, they are limited by the signal-to-noise ratio of the leaking acoustic signal, especially in chemical production environments where noise is complex and variable, leading to increased location errors. To address these issues, common approaches include using wavelet transform and blind source separation to eliminate interference and separate audio sources, thereby improving location accuracy; or using methods based on empirical mode decomposition and cross-time spectrum to improve feature extraction accuracy, thus enhancing location precision. However, these methods are only effective for specific or a few noise types, and are insufficient to achieve the leak location results required by engineering standards for the complex and variable non-stationary noise in chemical plants. Summary of the Invention
[0004] Based on this, the present invention provides a pipeline leak location method, device, chip and terminal, which can solve the problem that existing acoustic detection and location methods cannot effectively improve the leak location effect in pipeline leak location.
[0005] In a first aspect, a pipeline leak location method is provided, which is applied to a leak location system. The leak location system includes a first sensor, a second sensor, and a data acquisition unit connected to the first sensor and the second sensor. The first sensor is set at a first preset position on the outer surface of the target pipeline, and the second sensor is set at a second preset position on the outer surface of the target pipeline. The straight line where the first sensor and the second sensor are located coincides with the vertical projection of the centerline of the target pipeline. The distance between the first preset position and the second preset position is greater than 0.
[0006] The pipeline leak location method includes:
[0007] S11: Obtain the attribute parameters of the target pipe, including the average sound velocity in the target pipe, the sound velocity values based on I frequency bands, and the weight coefficient of each frequency band. The weight coefficient of the i-th frequency band is used to adjust the weight of the leak location result based on the i-th frequency band in the leak location result based on I frequency bands. The average sound velocity in the target pipe is solved according to the linear relationship between the leak location and the time delay.
[0008] S12: Read the first sound signal collected by the first sensor and the second sound signal collected by the second sensor from the data acquisition unit; the time delay in S11 is the difference between the time of receiving the first sound signal and the time of receiving the second sound signal;
[0009] S13: Divide the first sound signal and the second sound signal into I frequency bands of equal width; in the i-th frequency band, restore the signal to the time domain signal and divide it into different time periods to obtain the segmented first sound signal and the segmented second sound signal, and generate J first time delay intervals from different time periods by cross-correlation technology.
[0010] S14: The J first time delay intervals of the i-th frequency band are initially filtered by the average sound speed. After the initial filtering, a second filtering is performed. The principle of the second filtering is to retain the first time delay interval with the highest frequency as the actual sound signal delay interval.
[0011] S15: For the i-th frequency band, the leakage location is determined based on the sound velocity value of the i-th frequency band, the linear relationship, and the actual sound signal delay interval of the i-th frequency band, to obtain the leakage location result based on I frequency bands;
[0012] S16: Based on the weighting coefficients, perform a weighted average processing on the leakage location results based on I frequency bands to obtain the leakage location of the target pipeline;
[0013] Where I and J are positive integers, and i is a positive integer less than or equal to I.
[0014] Optionally, prior to S11, the following are included:
[0015] Check whether the target pipeline is installing a leak location system for the first time;
[0016] If so, proceed to step S11;
[0017] If not, measure the property parameters of the target pipeline.
[0018] Optionally, the property parameters of the target pipeline are measured, including:
[0019] S21: Create a simulated sound wave at the marked position of the target pipe, and read the first simulated sound wave signal collected by the first sensor and the second simulated sound wave signal collected by the second sensor from the data acquisition unit;
[0020] S22: Repeat S1 M times to obtain M sets of simulation data, the simulation data including the marker position, the first simulated sound wave signal, and the second simulated sound wave signal, where M is a positive integer;
[0021] S23: Fit a linear relationship between the marker position and the time delay based on the M sets of simulated data, wherein the time delay is the difference between the time of receiving the first simulated acoustic signal and the time of receiving the second simulated acoustic signal;
[0022] S24: Divide the first analog sound wave signal and the second analog sound wave signal into I frequency bands of equal width. In the i-th frequency band, restore the signal to the time domain signal and divide it into different time periods to obtain the divided first analog sound wave signal and the divided second analog sound wave signal. Use cross-correlation technology to associate and generate J second time delay intervals from different time periods.
[0023] S25: Calculate the average sound velocity in the target pipe and the sound velocity value based on I frequency bands according to the linear relationship;
[0024] S26: The J second time delay intervals of the i-th frequency band are initially filtered by the average sound speed. After the initial filtering, a second filtering is performed. The filtering principle of the second filtering is to retain the second time delay interval with the highest frequency as the simulated weight calculation index.
[0025] S27: For the i-th frequency band, the sound velocity value of the i-th frequency band, the linear relationship, and the simulated weight calculation index of the i-th frequency band are used to locate the leakage location, and the simulated positioning error rate of the i-th frequency band is recorded through the marked location. The simulated positioning error rate of the i-th frequency band is used to calculate the weight coefficient of the i-th frequency band.
[0026] Optionally, the linear relationship between the leak location and the time delay is calculated using the following formula:
[0027]
[0028] Where L is the distance between the first preset position and the second preset position, c is the average sound speed in the target pipe, Δt is the time delay, and x is the leak location.
[0029] Optionally, in step S13, the formula for calculating any first time delay interval is:
[0030]
[0031] Where R(n) is the CC function, s1(m+n) is the first sound signal, s2(m) is the second sound signal in the i-th frequency band, m represents the discrete time sequence index of the signal, n is the number of the first time delay intervals of the two sound sequences in any time period within the i-th frequency band, N is the total length of the signal, and the time delay between the segmented first sound signal and the segmented second sound signal is equal to the value of n when the CC function R(n) reaches its maximum value.
[0032] Optionally, in step S14, the preliminary filtering of the J first time delay intervals of the i-th frequency band using the average sound velocity includes:
[0033]
[0034] Where L is the distance between the first preset position and the second preset position, and F s The sampling frequency is represented by c, and the average speed of sound in the target pipe is c.
[0035] Optionally, S15 includes:
[0036] The time delay between the segmented first and second audio signals, represented by the actual audio signal delay interval, is as follows:
[0037]
[0038] Among them, F s Indicates the sampling frequency;
[0039] For the i-th frequency band, when locating the leakage position based on the sound velocity value of the i-th frequency band, the aforementioned linear relationship, and the actual sound signal delay interval number of the i-th frequency band, n is used. i This represents the actual audio signal delay interval in the i-th frequency band.
[0040] Where, n ij R represents the number of the first time delay intervals in the i-th frequency band and the j-th time period. ij The number of time delay intervals when the CC function in the i-th frequency band and j-th time period reaches its maximum value is n. ij MF represents selecting the most frequent of the J first time delay intervals from different time periods as the actual audio signal delay interval number for the i-th frequency band, where j is a positive integer less than or equal to J;
[0041] Based on the linear relationship between the leak location and the time delay, the leak location is located using the i-th frequency band.
[0042] The second aspect provides a pipeline leak location device for use in a leak location system. The leak location system includes a first sensor, a second sensor, and a data acquisition unit connected to the first sensor and the second sensor. The first sensor is set at a first preset position on the outer surface of the target pipeline, and the second sensor is set at a second preset position on the outer surface of the target pipeline. The straight line where the first sensor and the second sensor are located coincides with the vertical projection of the center line of the target pipeline. The distance between the first preset position and the second preset position is greater than 0.
[0043] The pipeline leak location device includes:
[0044] The attribute parameter acquisition module is used to acquire the attribute parameters of the target pipe. The attribute parameters include the average sound velocity in the target pipe, the sound velocity values based on I frequency bands, and the weight coefficient of each frequency band. The weight coefficient of the i-th frequency band is used to adjust the weight of the leak location result based on the i-th frequency band in the leak location result based on I frequency bands. The average sound velocity in the target pipe is calculated based on the linear relationship between the leak location and the time delay.
[0045] The sound signal reading module is used to read the first sound signal collected by the first sensor and the second sound signal collected by the second sensor from the data acquisition unit; the time delay in the attribute parameter acquisition module is the difference between the time of receiving the first sound signal and the time of receiving the second sound signal.
[0046] The signal segmentation module is used to segment both the first sound signal and the second sound signal into I frequency bands of equal width; after restoring the signal to the time domain signal in the i-th frequency band, it is segmented into different time periods to obtain the segmented first sound signal and the segmented second sound signal, and the J first time delay intervals from different time periods are generated by cross-correlation technology.
[0047] The filtering module is used to perform preliminary filtering on the J first time delay intervals of the i-th frequency band based on the average sound speed. After the preliminary filtering, a second filtering is performed. The principle of the second filtering is to retain the first time delay interval with the highest frequency as the actual sound signal delay interval.
[0048] The i-th frequency band leakage location module is used to locate the leakage location for the i-th frequency band based on the sound velocity value of the i-th frequency band, the linear relationship, and the actual sound signal delay interval of the i-th frequency band, and obtain leakage location results based on I frequency bands.
[0049] The target pipeline leakage location calculation module is used to process the leakage location results based on I frequency bands by weighted average according to the weight coefficient, and obtain the leakage location of the target pipeline.
[0050] Thirdly, a chip is provided, including a first processor for calling and running a computer program from a first memory, causing a device equipped with the chip to perform the steps of the pipeline leak location method described above.
[0051] Fourthly, a terminal is provided, including a second memory, a second processor, and a computer program stored in the second memory and executable on the second processor. When the second processor executes the computer program, it implements the various steps of the pipeline leak location method described above.
[0052] The aforementioned pipeline leak location method, device, chip, and terminal first acquire the attribute parameters of the target pipeline. Then, they read the current leak acoustic data (i.e., the first and second sound signals) collected by the first and second sensors. These signals are then divided into I equal-width frequency bands according to their spectrum. Cross-correlation techniques are used to correlate these bands and generate J first time delay intervals to calculate the difference between the time of receiving the first and second sound signals. This, combined with the attribute parameters, allows for leak location location based on the data in the i-th frequency band. After I iterations of leak location location, a leak location result based on the I frequency bands is obtained. Finally, a weighted average of the leak location results based on the I frequency bands is applied to obtain the leak location of the target pipeline. The first and second sound signals in different frequency bands have different statistical characteristics such as correlation, time delay, sound velocity, and frequency components. This embodiment of the invention uses weight coefficients for each frequency band to represent these statistical characteristics, thereby adaptively separating the leak signal from high background noise and effectively solving the problem of large location errors caused by the dispersion of the leak sound signal. Attached Figure Description
[0053] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0054] Figure 1 This is a schematic diagram of the basic process of the pipeline leak location method according to an embodiment of the present invention;
[0055] Figure 2 This is a schematic diagram of the structure of the leak location system according to an embodiment of the present invention;
[0056] Figure 3 This is a schematic diagram of the basic process of the pipeline leak location method according to an embodiment of the present invention;
[0057] Figure 4 This is a schematic diagram of the complete process of the pipeline leak location method according to an embodiment of the present invention;
[0058] Figure 5 This is a basic structural block diagram of the pipeline leakage locating device according to an embodiment of the present invention;
[0059] Figure 6 This is a basic structural block diagram of a terminal provided in an embodiment of the present invention. Detailed Implementation
[0060] To enable those skilled in the art to better understand the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings.
[0061] In some of the processes described in the specification, claims, and accompanying drawings of this invention, multiple operations appearing in a specific order are included. However, it should be clearly understood that these operations may not be executed in the order they appear herein, or may be executed in parallel. The operation numbers, such as 101, 102, etc., are merely used to distinguish different operations and do not represent any execution order. Furthermore, these processes may include more or fewer operations, and these operations may be executed sequentially or in parallel. It should be noted that the descriptions such as "first," "second," etc., in this document are used to distinguish different messages, devices, modules, etc., and do not represent a sequential order, nor do they limit "first" and "second" to different types.
[0062] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.
[0063] The embodiments of this application can acquire and process relevant data based on artificial intelligence technology. Artificial intelligence (AI) is the theory, method, technology, and application system that uses digital computers or machines controlled by digital computers to simulate, extend, and expand human intelligence, perceive the environment, acquire knowledge, and use that knowledge to obtain optimal results.
[0064] Foundational technologies for artificial intelligence generally include sensors, dedicated AI chips, cloud computing, distributed storage, big data processing, operating / interactive systems, and mechatronics. AI software technologies mainly encompass computer vision, robotics, biometrics, speech processing, natural language processing, and machine learning / deep learning.
[0065] Please refer to the details. Figure 1 , Figure 1 This is a basic flowchart illustrating the pipeline leak location method in this embodiment, which is applied to, for example... Figure 2 The leak location system 100 shown includes at least a first sensor 101, a second sensor 102, and a data acquisition unit 103 connected to the first sensor 101 and the second sensor 102. Figure 2In this configuration, the first sensor 101 is set at a first preset position D1 on the outer surface of the target pipe A, and the second sensor 102 is set at a second preset position D2 on the outer surface of the target pipe A. The straight line containing the first sensor 101 and the second sensor 102, such as L1, coincides with the vertical projection of the center line of the target pipe A, such as L2. It should be noted that when the first sensor 101 and the second sensor 102 are deployed, the distance L between the first preset position D1 and the second preset position D2 must be greater than 0.
[0066] In one embodiment, the first sensor 101 and the second sensor 102 are acceleration sensors.
[0067] like Figure 1 As shown, a method for locating pipeline leaks includes:
[0068] S11: Obtain the attribute parameters of the target pipe, including the average sound velocity in the target pipe, the sound velocity values based on I frequency bands, and the weighting coefficient of each frequency band.
[0069] The weighting coefficient of the i-th frequency band is used to adjust the weight of the leak location result based on the i-th frequency band in the leak location result based on the I frequency bands; the average sound velocity in the target pipe is solved according to the linear relationship between the leak location and the time delay.
[0070] It should be noted that when calculating the average sound velocity inside the target pipe, the leak location is known, such as through simulating the leak location. Furthermore, the aforementioned time delay is the difference between the time it takes to receive the first sound signal and the time it takes to receive the second sound signal.
[0071] In this embodiment of the invention, if a leak occurs at a distance x from the first sensor, the sound waves generated at the leak point will propagate along both directions of the pipe towards the two sensors at a certain speed c. The propagation of the sound waves can be described by the following two equations:
[0072] x = c·t1;
[0073] For an acoustic wave directed toward the first sensor, where t1 is the time required for the wave to cover the distance between the leak point and the location of the first sensor, we have:
[0074] Lx = c·t2;
[0075] For the sound wave propagating towards the second sensor, where t2 is the time required for the sound wave (from the leak point) to reach the second sensor. Therefore, by subtracting (1) and (2) to solve for x, the linear relationship between the leak location and the time delay is obtained, and the calculation formula is:
[0076]
[0077] Where L is the distance between the first preset position and the second preset position, c is the average sound speed in the target pipe, Δt is the time delay, and x is the leak location.
[0078] Based on the above linear relationship, for the sound velocity values of I frequency bands, the sound velocity value of the i-th frequency band is c. i The calculation method is as follows:
[0079]
[0080] Where, Δt i Time delay of the i-th frequency band.
[0081] In this embodiment of the invention, the average sound velocity, the sound velocity value based on I frequency bands, and the weighting coefficient of each frequency band are attribute parameters unique to the target pipeline. If the leak location system 100 is installed on other pipelines, then those other pipelines are the target pipelines, and the attribute parameters of different pipelines are all different. Therefore, if the leak location system 100 is installed on a certain pipeline for the first time, the attribute parameters of this pipeline need to be measured and stored so that the attribute parameters can be directly called when performing pipeline leak location on this pipeline.
[0082] S12: Read the first sound signal collected by the first sensor and the second sound signal collected by the second sensor from the data acquisition unit.
[0083] It is understandable that before step S12 above, the sampling frequency of the data acquisition unit is also set.
[0084] S13: Divide the first sound signal and the second sound signal into I frequency bands of equal width; in the i-th frequency band, restore the signal to the time domain signal and divide it into different time periods to obtain the segmented first sound signal and the segmented second sound signal, and generate J first time delay intervals from different time periods by cross-correlation technology.
[0085] In step S13 above, the signal is first divided into different frequency bands. Within each frequency band, the signal is restored to a time-domain signal and then divided into different time periods. Signals from different time periods within the same frequency band are cross-correlated to obtain the first time delay interval for each time period. In this embodiment of the invention, the sound signal is divided into multiple frequency bands, allowing for individual processing of each frequency band, which is adaptable to the complex and variable environment in chemical pipelines.
[0086] It should be noted that, based on the linear relationship between the leak location and the time delay, the time delay is a crucial parameter for determining the leak location and thus locating the pipeline leak. Therefore, this embodiment of the invention treats pipeline leak location as a problem of finding the leak location, and the answer comes from the time delay between the signals received by the two sensors. Since the first and second sound signals are two similar but time-shifted discrete-time signals, this embodiment of the invention uses the CC (CrossCorrelation) function to calculate the aforementioned time delay.
[0087] In the i-th frequency band, the first segmented audio signal and the second segmented audio signal are first obtained. Then, the first segmented audio signal and the second segmented audio signal are correlated using cross-correlation technology to generate multiple first time delay intervals from different time periods. For any first time delay interval, the calculation formula is as follows:
[0088]
[0089] Where R(n) is the CC function, s1(m+n) is the first sound signal, s2(m) is the second sound signal in the i-th frequency band, m represents the discrete time sequence index of the signal, n is the number of the first time delay intervals of the two sound sequences in any time period within the i-th frequency band, N is the total length of the signal, and the time delay between the segmented first sound signal and the segmented second sound signal is equal to the value of n when the CC function R(n) reaches its maximum value.
[0090] S14: The J first time delay intervals of the i-th frequency band are initially filtered by the average sound speed. After the initial filtering, a second filtering is performed. The principle of the second filtering is to retain the first time delay interval with the highest frequency as the actual sound signal delay interval.
[0091] In practical applications, leakage location is determined for each frequency band, and the obtained leakage location (x) value should be between 0 and L. Therefore, in step S14 above, the initial filtering principle is to reject outliers of n, which can be expressed by the formula:
[0092]
[0093] Where L is the distance between the first preset position and the second preset position, and F s The sampling frequency is represented by c, and the average speed of sound in the target pipe is c.
[0094] In step S14 above, the first time delay interval with the highest frequency is taken as the actual sound signal delay interval, which means that the time period in which the most frequently occurring value is located is selected as the representative of the i-th frequency band.
[0095] S15: For the i-th frequency band, the leakage location is determined based on the sound velocity value of the i-th frequency band, the linear relationship, and the actual sound signal delay interval of the i-th frequency band, thus obtaining the leakage location result based on I frequency bands.
[0096] In step S15 above, after I times of leak location localization, a leak location localization result based on I frequency bands is obtained. Furthermore, the leak location is determined based on the sound velocity value of the i-th frequency band, the linear relationship, and the actual sound signal delay interval of the i-th frequency band. This fully considers the time delay estimation and corresponding sound wave velocity under different frequency bands, and can effectively solve the problem of leak location under high background noise and sound wave dispersion.
[0097] It should be noted that the actual audio signal delay interval number in the i-th frequency band, or the value of the first time delay interval n when the CC function R(n) reaches its maximum value, cannot be directly equated to the time delay. Therefore, S15 includes:
[0098] The time delay between the segmented first and second audio signals, represented by the actual audio signal delay interval, is as follows:
[0099]
[0100] Among them, F s Indicates the sampling frequency;
[0101] For the i-th frequency band, when locating the leakage position based on the sound velocity value of the i-th frequency band, the aforementioned linear relationship, and the actual sound signal delay interval number of the i-th frequency band, n is used. i This represents the actual audio signal delay interval in the i-th frequency band.
[0102] Where, n ij R represents the number of the first time delay intervals in the i-th frequency band and the j-th time period. ij The number of time delay intervals when the CC function in the i-th frequency band and j-th time period reaches its maximum value is n. ij MF represents selecting the most frequent of the J first time delay intervals from different time periods as the actual audio signal delay interval number for the i-th frequency band, where j is a positive integer less than or equal to J;
[0103] Based on the linear relationship between the leak location and the time delay, the leak location is located based on the i-th frequency band, and finally the leak location result based on I frequency bands is obtained.
[0104] For example, x is located based on the leakage location of the i-th frequency band. i as follows:
[0105]
[0106] Among them, c i Let be the sound velocity value for the i-th frequency band.
[0107] S16: Based on the weighting coefficients, perform a weighted average processing on the leakage location results based on I frequency bands to obtain the leakage location of the target pipeline;
[0108] Where I and J are positive integers, and i is a positive integer less than or equal to I.
[0109] According to step S11 above, the weighting coefficient of the i-th frequency band is used to adjust the weight of the leakage location location result based on the i-th frequency band in the leakage location location result based on I frequency bands. The weighting coefficient represents the reliability of each band.
[0110] Combining step S15 above, the final value for locating the leak location in step S16 above, i.e., the leak location x of the target pipeline, is:
[0111]
[0112]
[0113] Through the above steps S11 to S16, including frequency band segmentation and time segmentation of the sound signal, the frequency band segmentation is actually a spectrum-based segmentation, which helps to reduce the negative impact of sound dispersion on positioning accuracy. The time segmentation is actually a segmentation of the time-domain leakage signal, which helps to improve the randomness of processing such signals. The weighting coefficient increases the adaptability of the leakage location method, making it applicable to pipeline systems in different environments.
[0114] This embodiment of the invention also describes the method for measuring the attribute parameters of a pipeline when the leak location system 100 is first installed on a certain pipeline.
[0115] Combining steps S11 to S16 above, the pipeline leak location method provided in this embodiment of the invention includes, before step S11:
[0116] Check whether the target pipeline is installing a leak location system for the first time;
[0117] If so, proceed to step S11;
[0118] If not, measure the property parameters of the target pipeline.
[0119] like Figure 3 As shown, the property parameters of the target pipeline are measured, including:
[0120] S21: Create a simulated sound wave at the marked position of the target pipe, and read the first simulated sound wave signal collected by the first sensor and the second simulated sound wave signal collected by the second sensor from the data acquisition unit.
[0121] In step S21 above, the marked location simulates the location of a pipe leak. The simulated sound wave can be generated by producing certain impacts at the assumed pipe leak location, i.e., the marked location, and then collected by the first and second sensors.
[0122] S22: Repeat S1 M times to obtain M sets of simulation data, the simulation data including the marker position, the first simulated sound wave signal, and the second simulated sound wave signal, where M is a positive integer.
[0123] S23: Fit a linear relationship between the marker position and the time delay based on the M sets of simulated data, wherein the time delay is the difference between the time of receiving the first simulated sound wave signal and the time of receiving the second simulated sound wave signal.
[0124] Based on the linear relationship between the leak location and the time delay, in step S23 above, the calculation formula for the linear relationship between the marked location and the time delay fitted from the M sets of simulation data is as follows:
[0125]
[0126] Where L is the distance between the first preset position and the second preset position, c is the average sound speed in the target pipe, Δt is the time delay, and x is the marker position.
[0127] It is understandable that x represents the location of the leak, which means that the location of the leak is considered a known number. Therefore, the average sound speed in the target pipe is an unknown number. Here, a linear relationship is used to solve for the average sound speed in the target pipe.
[0128] S24: Divide the first analog sound wave signal and the second analog sound wave signal into I frequency bands of equal width. In the i-th frequency band, after restoring the signal to the time domain signal, divide it into different time periods to obtain the divided first analog sound wave signal and the divided second analog sound wave signal. Use cross-correlation technology to associate and generate J second time delay intervals from different time periods.
[0129] The second time delay interval number obtained in step S24 is the same as the first time delay interval number obtained in step S13, and is used to calculate the time delay between the signals received by the two sensors, such as the first sound signal and the second sound signal, the first analog sound wave signal and the second analog sound wave signal.
[0130] S25: Calculate the average sound velocity within the target pipe and the sound velocity values based on I frequency bands according to the linear relationship.
[0131] Based on the above linear relationship, for the sound velocity values of I frequency bands, the sound velocity value of the i-th frequency band is c. i The calculation method is as follows:
[0132]
[0133] Where, Δt i Time delay of the i-th frequency band.
[0134] S26: The J second time delay intervals of the i-th frequency band are initially filtered using the average sound speed. After the initial filtering, a second filtering is performed. The principle of the second filtering is to retain the second time delay interval with the highest frequency as the simulated weight calculation index.
[0135] In step S26 above, the second time delay interval number that appears most frequently is used as the simulation weight calculation index to represent the i-th frequency band and calculate the time delay. The time delay is an important parameter for solving the leak location and completing the pipeline leak location. In step S21 above, the pipeline leak location is simulated by marking the location. Therefore, as in step S27 below, the embodiment of the present invention evaluates the simulation weight calculation index to obtain the weight coefficients used in steps S11 to S16 above.
[0136] S27: For the i-th frequency band, the sound velocity value of the i-th frequency band, the linear relationship, and the simulated weight calculation index of the i-th frequency band are used to locate the leakage location, and the simulated positioning error rate of the i-th frequency band is recorded through the marked location. The simulated positioning error rate of the i-th frequency band is used to calculate the weight coefficient of the i-th frequency band.
[0137] Step S27 above records the simulated positioning error rate of the i-th frequency band by marking the position, which fully reflects the influence factors of each frequency band after frequency band segmentation. That is, the influence of high background noise under different environments and different pipeline conditions on the statistical characteristics such as correlation, time delay, sound velocity value and frequency components in the i-th frequency band. In this embodiment of the invention, the signals of different time periods are weighted according to the calculated weight coefficients, so as to adapt to the leakage positioning of different pipelines in complex and ever-changing chemical environments.
[0138] For example, the formula for calculating the weighting coefficient of the i-th frequency band is:
[0139]
[0140] Where ω(i) is the weighting coefficient of the i-th frequency band, and err(i) is the error rate of the i-th frequency band when estimating the reference leakage location. The inverse square root in the above formula can prevent the difference between the values of the weighting coefficients from becoming too high, thus avoiding the disappearance of some low-coefficient frequency bands.
[0141] Furthermore, it should be noted that in S24, the segmented first and second analog sound wave signals are correlated using cross-correlation technology to generate J second time delay intervals, which is the same as in S13 where the segmented first and second sound signals are correlated using cross-correlation technology to generate J first time delay intervals; in S26, the J second time delay intervals of the i-th frequency band are initially filtered using the average sound velocity, which is the same as in S14 where the J first time delay intervals of the i-th frequency band are initially filtered using the average sound velocity; in S27, the leakage location is determined by calculating the exponent using the sound velocity value of the i-th frequency band, the linear relationship, and the analog weight of the i-th frequency band, which is the same as in S15 where the leakage location is determined by the sound velocity value of the i-th frequency band, the linear relationship, and the actual sound signal delay intervals of the i-th frequency band. Therefore, the detailed implementation of the above steps will not be elaborated here.
[0142] Combination Figure 1 and Figure 3 The pipeline leak location method shown is as follows: Figure 4 As shown in the figure, the embodiment of the present invention also shows a flowchart of the implementation of steps S11 to S16 and steps S21 to S27. Figure 4 Steps S21 to S27 are referred to as the initial stage, and steps S11 to S16 are referred to as the main stage. In the initial stage, the main purpose is to measure the attribute parameters of the target pipe, such as the calculation of the average sound velocity of the target pipe, the calculation of the weighting coefficient, and the sound velocity value of each frequency band. In the main stage, the primary objective is to locate pipeline leaks. First, the attribute parameters measured in the initial stage are retrieved. Then, the sound signal is divided into I frequency bands of equal width, identical to those in the initial stage. Within the i-th frequency band, the signal is restored to a time-domain signal and then divided into different time periods, yielding the first and second segmented sound signals. These are correlated using cross-correlation techniques to generate J first time delay intervals from different time periods. These J first time delay intervals are filtered twice to obtain crucial parameters for pipeline leak location: the actual sound signal delay intervals used to calculate time delays. Therefore, the leak location can be determined based on the sound velocity value, linear relationship, and actual sound signal delay intervals of the i-th frequency band. After I location operations, leak location results based on I frequency bands are obtained. Finally, the leak location results based on I frequency bands are weighted and averaged according to weighting coefficients to obtain the leak location of the target pipeline.
[0143] To address the aforementioned technical problems, embodiments of the present invention also provide a pipeline leak location device. Please refer to [link / reference needed] for details. Figure 5 , Figure 5 This is a basic structural block diagram of the pipeline leak locating device implemented in this embodiment, which is applied to, for example... Figure 2 The leak location system shown includes a pipeline leak location device 50 comprising:
[0144] The attribute parameter acquisition module 51 is used to acquire the attribute parameters of the target pipe. The attribute parameters include the average sound velocity in the target pipe, the sound velocity values based on I frequency bands, and the weight coefficient of each frequency band. The weight coefficient of the i-th frequency band is used to adjust the weight of the leak location result based on the i-th frequency band in the leak location result based on I frequency bands. The average sound velocity in the target pipe is calculated based on the linear relationship between the leak location and the time delay.
[0145] The sound signal reading module 52 is used to read the first sound signal collected by the first sensor and the second sound signal collected by the second sensor from the data acquisition unit; the time delay in the attribute parameter acquisition module is the difference between the time of receiving the first sound signal and the time of receiving the second sound signal.
[0146] The signal segmentation module 53 is used to segment the first sound signal and the second sound signal into I frequency bands of equal width; after restoring the signal to the time domain signal in the i-th frequency band, it is segmented into different time periods to obtain the segmented first sound signal and the segmented second sound signal, and the J first time delay intervals from different time periods are generated by cross-correlation technology.
[0147] The filtering module 54 is used to perform preliminary filtering on the J first time delay intervals of the i-th frequency band based on the average sound speed. After the preliminary filtering, a second filtering is performed. The principle of the second filtering is to retain the first time delay interval with the highest frequency as the actual sound signal delay interval.
[0148] The i-th frequency band leakage location module 55 is used to locate the leakage location for the i-th frequency band based on the sound velocity value of the i-th frequency band, the linear relationship and the actual sound signal delay interval of the i-th frequency band, and obtain leakage location results based on I frequency bands.
[0149] The target pipeline leakage location calculation module 56 is used to process the leakage location results based on I frequency bands by weighted average according to the weight coefficient, and obtain the leakage location of the target pipeline.
[0150] To address the aforementioned technical problems, embodiments of the present invention also provide a chip, which can be a general-purpose processor or a dedicated processor. The chip includes a processor that supports a terminal in executing the aforementioned related steps, such as retrieving and running a computer program from memory, causing a device equipped with the chip to execute the pipeline leak location method described in the various embodiments above.
[0151] Optionally, in some examples, the chip may also include a transceiver for receiving control from the processor to support the terminal in performing the aforementioned steps to implement the pipeline leak location method in the various embodiments described above.
[0152] Optionally, the chip may also include a storage medium.
[0153] It should be noted that the chip can be implemented using the following circuits or devices: one or more field programmable gate arrays (FPGAs), programmable logic devices (PLDs), controllers, state machines, gate logic, discrete hardware components, any other suitable circuits, or any combination of circuits capable of performing the various functions described throughout this application.
[0154] The present invention also provides a terminal, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the pipeline leak location method as described in any one of claims 1 to 7.
[0155] Please refer to the details. Figure 6 , Figure 6 The diagram illustrates a basic structural block diagram of a terminal, which includes a processor, a non-volatile storage medium, a memory, and a network interface connected via a system bus. The non-volatile storage medium stores an operating system, a database, and computer-readable instructions. The database may store a sequence of control information. When executed by the processor, the computer-readable instructions enable the processor to implement a pipeline leak location method. The processor provides computational and control capabilities to support the operation of the entire terminal. The memory stores computer-readable instructions, which, when executed by the processor, enable the processor to perform a pipeline leak location method. The network interface is used for communication with the terminal. Those skilled in the art will understand that the structure shown in the diagram is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the terminal to which the present application is applied. A specific terminal may include more or fewer components than shown in the diagram, or combine certain components, or have different component arrangements.
[0156] Those skilled in the art will understand that the terms "terminal" and "terminal device" as used herein include both devices that receive wireless signals, devices that only possess wireless signal receiver capabilities without transmission capabilities, and devices with receiving and transmitting hardware, electronic devices having receiving and transmitting hardware capable of performing bidirectional communication on a bidirectional communication link. Such electronic devices may include: cellular or other communication devices having a single-line display, a multi-line display, or a cellular or other communication device without a multi-line display; PCS (Personal Communications Service) that can combine voice, data processing, fax, and / or data communication capabilities; PDA (Personal Digital Assistant) that may include a radio frequency receiver, pager, Internet / intranet access, web browser, notepad, calendar, and / or GPS (Global Positioning System) receiver; and conventional laptop and / or handheld computers or other devices that have and / or include radio frequency receivers. As used herein, "terminal" or "terminal device" can be portable, transportable, installed in a means of transportation (air, sea, and / or land), or suitable and / or configured to operate locally and / or in a distributed manner, operating in any other location on Earth and / or in space. "Terminal" or "terminal device" as used herein can also be a communication terminal, an internet access terminal, or a music / video playback terminal, such as a PDA, a MID (Mobile Internet Device), and / or a mobile phone with music / video playback capabilities, or a smart TV, set-top box, etc.
[0157] The present invention also provides a storage medium storing computer-readable instructions, which, when executed by one or more processors, cause the one or more processors to perform the steps of the pipeline leak location method described in any of the above embodiments.
[0158] This embodiment also provides a computer program that can be distributed on a computer-readable medium and executed by a computing device to implement at least one step of the pipeline leak location method described above; and in some cases, at least one step shown or described may be executed in a different order than that described in the above embodiments.
[0159] This embodiment also provides a computer program product, including a computer-readable device on which the computer program as shown above is stored. In this embodiment, the computer-readable device may include the computer-readable storage medium as shown above.
[0160] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. This computer program can be stored in a computer-readable storage medium, and when executed, it can include the processes of the embodiments of the methods described above. The aforementioned storage medium can be a non-volatile storage medium such as a magnetic disk, optical disk, or read-only memory (ROM), or random access memory (RAM).
[0161] It should be understood that although the steps in the flowcharts of the accompanying figures are shown sequentially as indicated by the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the accompanying figures may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times, and their execution order is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the sub-steps or stages of other steps.
[0162] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.
[0163] The embodiments described above are merely illustrative of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention. Therefore, the scope of protection of this patent should be determined by the appended claims.
Claims
1. A method for locating pipeline leaks, characterized in that, The system is used in a leak location system, which includes a first sensor, a second sensor, and a data acquisition unit connected to the first sensor and the second sensor. The first sensor is set at a first preset position on the outer surface of the target pipe, and the second sensor is set at a second preset position on the outer surface of the target pipe. The straight line where the first sensor and the second sensor are located coincides with the vertical projection of the center line of the target pipe. The distance between the first preset position and the second preset position is greater than 0. The pipeline leak location method includes: S11: Obtain the attribute parameters of the target pipe, including the average sound velocity in the target pipe, the sound velocity values based on I frequency bands, and the weight coefficient of each frequency band. The weight coefficient of the i-th frequency band is used to adjust the weight of the leak location result based on the i-th frequency band in the leak location result based on I frequency bands. The average sound velocity in the target pipe is solved according to the linear relationship between the leak location and the time delay. S12: Read the first sound signal collected by the first sensor and the second sound signal collected by the second sensor from the data acquisition unit; the time delay in S11 is the difference between the time of receiving the first sound signal and the time of receiving the second sound signal; S13: Divide the first sound signal and the second sound signal into I frequency bands of equal width; in the i-th frequency band, restore the signal to the time domain signal and divide it into different time periods to obtain the segmented first sound signal and the segmented second sound signal, and generate J first time delay intervals from different time periods by cross-correlation technology. S14: The J first time delay intervals of the i-th frequency band are initially filtered by the average sound speed. After the initial filtering, a second filtering is performed. The principle of the second filtering is to retain the first time delay interval with the highest frequency as the actual sound signal delay interval. S15: For the i-th frequency band, the leakage location is determined based on the sound velocity value of the i-th frequency band, the linear relationship, and the actual sound signal delay interval of the i-th frequency band, to obtain the leakage location result based on I frequency bands; S16: Based on the weighting coefficients, perform a weighted average processing on the leakage location results based on I frequency bands to obtain the leakage location of the target pipeline; Where I and J are positive integers, and i is a positive integer less than or equal to I.
2. The pipeline leak location method as described in claim 1, characterized in that, Prior to S11, including: Check whether the target pipeline is installing a leak location system for the first time; If so, proceed to step S11; If not, measure the property parameters of the target pipeline.
3. The pipeline leak location method as described in claim 2, characterized in that, Measuring the property parameters of the target pipeline, including: S21: Create a simulated sound wave at the marked position of the target pipe, and read the first simulated sound wave signal collected by the first sensor and the second simulated sound wave signal collected by the second sensor from the data acquisition unit; S22: Repeat S1 M times to obtain M sets of simulation data, the simulation data including the marker position, the first simulated sound wave signal, and the second simulated sound wave signal, where M is a positive integer; S23: Fit a linear relationship between the marker position and the time delay based on the M sets of simulated data, wherein the time delay is the difference between the time of receiving the first simulated acoustic signal and the time of receiving the second simulated acoustic signal; S24: Divide the first analog sound wave signal and the second analog sound wave signal into I frequency bands of equal width. In the i-th frequency band, restore the signal to the time domain signal and divide it into different time periods to obtain the divided first analog sound wave signal and the divided second analog sound wave signal. Use cross-correlation technology to associate and generate J second time delay intervals from different time periods. S25: Calculate the average sound velocity in the target pipe and the sound velocity value based on I frequency bands according to the linear relationship; S26: The J second time delay intervals of the i-th frequency band are initially filtered by the average sound speed. After the initial filtering, a second filtering is performed. The filtering principle of the second filtering is to retain the second time delay interval with the highest frequency as the simulated weight calculation index. S27: For the i-th frequency band, the sound velocity value of the i-th frequency band, the linear relationship, and the simulated weight calculation index of the i-th frequency band are used to locate the leakage location, and the simulated positioning error rate of the i-th frequency band is recorded through the marked location. The simulated positioning error rate of the i-th frequency band is used to calculate the weight coefficient of the i-th frequency band.
4. The pipeline leak location method as described in claim 1, characterized in that, The linear relationship between the leak location and the time delay is calculated using the following formula: Where L is the distance between the first preset position and the second preset position, c is the average sound speed in the target pipe, Δt is the time delay, and x is the leak location.
5. The pipeline leak location method as described in claim 1, characterized in that, In step S13, the formula for calculating any first time delay interval is: Where R(n) is the CC function, s1(m+n) is the first sound signal, s2(m) is the second sound signal in the i-th frequency band, m represents the discrete time sequence index of the signal, n is the number of the first time delay intervals of the two sound sequences in any time period within the i-th frequency band, N is the total length of the signal, and the time delay between the segmented first sound signal and the segmented second sound signal is equal to the value of n when the CC function R(n) reaches its maximum value.
6. The pipeline leak location method as described in claim 1, characterized in that, In step S14, the initial filtering of the J first time delay intervals of the i-th frequency band using the average sound velocity includes: Where L is the distance between the first preset position and the second preset position, and F s The sampling frequency is represented by c, and the average speed of sound in the target pipe is c.
7. The pipeline leak location method as described in claim 1, characterized in that, S15 includes: The time delay between the segmented first and second audio signals, represented by the actual audio signal delay interval, is as follows: Among them, F s Indicates the sampling frequency; For the i-th frequency band, when locating the leakage position based on the sound velocity value of the i-th frequency band, the aforementioned linear relationship, and the actual sound signal delay interval number of the i-th frequency band, n is used. i This represents the actual audio signal delay interval in the i-th frequency band. Where, n ij R represents the number of the first time delay intervals in the i-th frequency band and the j-th time period. ij The number of time delay intervals when the CC function in the i-th frequency band and j-th time period reaches its maximum value is n. ij MF represents selecting the most frequent of the J first time delay intervals from different time periods as the actual audio signal delay interval number for the i-th frequency band, where j is a positive integer less than or equal to J; Based on the linear relationship between the leak location and the time delay, the leak location is located using the i-th frequency band.
8. A pipeline leak location device, characterized in that, The system is used in a leak location system, which includes a first sensor, a second sensor, and a data acquisition unit connected to the first sensor and the second sensor. The first sensor is set at a first preset position on the outer surface of the target pipe, and the second sensor is set at a second preset position on the outer surface of the target pipe. The straight line where the first sensor and the second sensor are located coincides with the vertical projection of the center line of the target pipe. The distance between the first preset position and the second preset position is greater than 0. The pipeline leak location device includes: The attribute parameter acquisition module is used to acquire the attribute parameters of the target pipe. The attribute parameters include the average sound velocity in the target pipe, the sound velocity values based on I frequency bands, and the weight coefficient of each frequency band. The weight coefficient of the i-th frequency band is used to adjust the weight of the leak location result based on the i-th frequency band in the leak location result based on I frequency bands. The average sound velocity in the target pipe is calculated based on the linear relationship between the leak location and the time delay. The sound signal reading module is used to read the first sound signal collected by the first sensor and the second sound signal collected by the second sensor from the data acquisition unit; the time delay in the attribute parameter acquisition module is the difference between the time of receiving the first sound signal and the time of receiving the second sound signal. The signal segmentation module is used to segment both the first sound signal and the second sound signal into I frequency bands of equal width; after restoring the signal to the time domain signal in the i-th frequency band, it is segmented into different time periods to obtain the segmented first sound signal and the segmented second sound signal, and the J first time delay intervals from different time periods are generated by cross-correlation technology. The filtering module is used to perform preliminary filtering on the J first time delay intervals of the i-th frequency band based on the average sound speed. After the preliminary filtering, a second filtering is performed. The principle of the second filtering is to retain the first time delay interval with the highest frequency as the actual sound signal delay interval. The i-th frequency band leakage location module is used to locate the leakage location for the i-th frequency band based on the sound velocity value of the i-th frequency band, the linear relationship, and the actual sound signal delay interval of the i-th frequency band, and obtain leakage location results based on I frequency bands. The target pipeline leakage location calculation module is used to process the leakage location results based on I frequency bands by weighted average according to the weight coefficient, and obtain the leakage location of the target pipeline.
9. A chip, characterized in that, include: A first processor is configured to retrieve and run a computer program from a first memory, causing a device equipped with the chip to perform the steps of the pipeline leak location method as described in any one of claims 1 to 7.
10. A terminal, characterized in that, The method includes a second memory, a second processor, and a computer program stored in the second memory and executable on the second processor, characterized in that the second processor, when executing the computer program, implements the steps of the pipeline leak location method as described in any one of claims 1 to 7.
Citation Information
Patent Citations
Method and system for remote sound speed measurement
US20100302907A1
Leak detector
US20120007744A1