CN-121994427-A - Infrasonic wave signal positioning analysis method and system

CN121994427ACN 121994427 ACN121994427 ACN 121994427ACN-121994427-A

Abstract

The invention relates to the technical field of safety monitoring and signal processing, in particular to a method and a system for positioning and analyzing an infrasonic wave signal, wherein the method comprises the steps of acquiring multichannel infrasonic wave signals, extracting multipath propagation characteristics, generating propagation compensation parameters, and completing signal correction to obtain an infrasonic wave signal set; the method comprises the steps of calculating an arrival time difference set by an infrasonic wave signal set, estimating an initial arrival azimuth, updating array geometric parameters based on the initial arrival azimuth and the arrival time difference set to estimate an arrival azimuth, carrying out weighted fusion according to consistency evaluation to obtain a fusion azimuth, inputting the fusion azimuth into an azimuth filter, updating noise parameters by combining a historical arrival azimuth reference sequence to output a refined arrival azimuth, forming an arrival time difference residual error based on the refined arrival azimuth, correcting to obtain a final arrival azimuth, fusing with pipeline geographic information to determine a leakage position, and improving positioning stability and high-precision positioning capability.

Inventors

HUANG YU
LIU JINLIANG
FU RUWEN
YUAN XINAN
LI WEI
LI XINHONG

Assignees

广东省特种设备检测研究院茂名检测院

Dates

Publication Date: 20260508
Application Date: 20260211

Claims (10)

1. The infrasonic wave signal positioning analysis method is characterized by comprising the following steps of: Acquiring multichannel infrasonic wave signals distributed in an oil and gas pipeline, extracting multipath propagation characteristics, determining propagation compensation parameters, and correcting the multichannel infrasonic wave signals based on the propagation compensation parameters to obtain an infrasonic wave signal set; calculating an arrival time difference set based on the infrasonic wave signal set, estimating an initial arrival azimuth based on the arrival time difference set, updating array geometric parameters based on the initial arrival azimuth and the arrival time difference set, estimating an arrival azimuth, and obtaining a fusion azimuth through weighted fusion; inputting the fusion azimuth into an azimuth filter, updating noise parameters of the azimuth filter based on a historical arrival azimuth reference sequence, and outputting a refined arrival azimuth; And calculating an arrival time difference residual error based on the refined arrival azimuth, correcting the refined arrival azimuth according to the arrival time difference residual error to obtain a final arrival azimuth, determining a leakage position based on the final arrival azimuth and pipeline geographic information, and outputting a leakage positioning result.
2. The method of claim 1, wherein extracting the multipath propagation characteristics comprises: Calculating a cross-correlation main peak width, a cross-correlation main peak to cross-correlation sub-peak energy ratio, and cross-channel coherence based on the multi-channel infrasonic signal to form the multipath propagation characteristics, determining the propagation compensation parameters includes determining an arrival time difference offset correction amount and a channel weighting coefficient based on the multipath propagation characteristics, correcting the multi-channel infrasonic signal based on the propagation compensation parameters includes performing an arrival time difference offset correction based on the arrival time difference offset correction amount and performing channel weighting based on the channel weighting coefficient to obtain the infrasonic signal set.
3. The infrasound signal localization analysis method of claim 1, wherein calculating the set of arrival time differences based on the set of infrasound signals comprises: and executing cross-correlation operation on any two channels of infrasonic signals in the infrasonic signal set, determining arrival time differences corresponding to cross-correlation main peaks, and collecting the arrival time differences to form the arrival time difference set.
4. The infrasonic wave signal localization analysis method of claim 1, wherein estimating the initial arrival azimuth based on the set of arrival time differences comprises: and constructing a arrival azimuth candidate set based on the array geometric parameters, calculating a theoretical arrival time difference set aiming at the arrival azimuth in the arrival azimuth candidate set, and taking the arrival azimuth with the smallest difference between the arrival time difference set and the theoretical arrival time difference set as the initial arrival azimuth.
5. The infrasonic wave signal localization analysis method of claim 1, wherein updating the array geometry parameters based on the initial arrival azimuth and the arrival time difference set includes: The method comprises the steps of constructing an array configuration set, wherein the array configuration set is a plurality of array geometric parameter candidate combinations generated based on sensor installation position constraint, calculating theoretical arrival time difference sets aiming at array configurations in the array configuration set, calculating differences between the arrival time difference sets and corresponding theoretical arrival time difference sets to obtain consistency evaluation values, and selecting array geometric parameters corresponding to the array configuration with the minimum consistency evaluation values as updated array geometric parameters.
6. The method of infrasonic wave signal localization analysis of claim 5, wherein obtaining the fused orientation by weighted fusion includes: and estimating the arrival azimuth corresponding to each array configuration in the array configuration set based on the array configuration set, determining configuration weight based on the consistency evaluation value corresponding to each array configuration, and carrying out weighted fusion on each arrival azimuth based on the configuration weight to obtain the fusion azimuth.
7. The method of claim 1, wherein the azimuth filter is a kalman filter, and the azimuth filter noise parameters include a process noise parameter and a measurement noise parameter.
8. The infrasonic wave signal localization analysis method of claim 7, wherein updating the azimuth filter noise parameters based on the historical arrival azimuth reference sequence includes: and forming an observation residual sequence by using a difference value between the fusion azimuth and the azimuth predicted by the azimuth filter, and updating the process noise parameter and the measurement noise parameter based on the observation residual sequence so as to output the refined arrival azimuth.
9. The infrasonic wave signal localization analysis method of claim 1, wherein calculating the arrival time difference residual based on the refined arrival direction includes: Calculating a theoretical arrival time difference set based on the refined arrival azimuth and the updated array geometric parameters, forming an arrival time difference residual by the difference between the arrival time difference set and the theoretical arrival time difference set, and correcting the refined arrival azimuth according to the arrival time difference residual, wherein the method comprises the steps of adopting a least square method to determine an arrival azimuth correction amount based on the arrival time difference residual; determining the leakage position based on the final arrival azimuth and the pipeline geographic information comprises taking the sensor installation position corresponding to the multichannel infrasonic wave signal as a monitoring point position, generating an azimuth line along the final arrival azimuth, and determining an intersection point of the azimuth line and a pipeline central line represented by the pipeline geographic information as the leakage position.
10. A system for locating and analyzing an infrasonic wave signal, for implementing the method for locating and analyzing an infrasonic wave signal according to any one of claims 1 to 9, characterized in that the system comprises: the signal correction module is used for acquiring multichannel infrasonic signals distributed in the oil and gas pipeline, extracting multipath propagation characteristics, determining propagation compensation parameters, and correcting the multichannel infrasonic signals based on the propagation compensation parameters to obtain an infrasonic signal set; The azimuth fusion module is used for calculating an arrival time difference set based on the infrasonic wave signal set, estimating an initial arrival azimuth based on the arrival time difference set, updating array geometric parameters based on the initial arrival azimuth and the arrival time difference set, estimating an arrival azimuth, and obtaining a fusion azimuth through weighted fusion; The filtering updating module is used for inputting the fusion azimuth into an azimuth filter, updating noise parameters of the azimuth filter based on a historical arrival azimuth reference sequence and outputting a refined arrival azimuth; and the positioning output module is used for calculating an arrival time difference residual error based on the refined arrival azimuth, correcting the refined arrival azimuth according to the arrival time difference residual error to obtain a final arrival azimuth, determining a leakage position based on the final arrival azimuth and pipeline geographic information and outputting a leakage positioning result.

Description

Infrasonic wave signal positioning analysis method and system Technical Field The invention relates to the technical field of safety monitoring and signal processing, in particular to an infrasonic wave signal positioning analysis method and system. Background The oil gas pipeline is laid in a long-distance cross-region mode, once a leakage event occurs, environmental pollution, fire explosion and transportation stoppage loss are easy to cause, the industry puts higher requirements on early discovery and quick positioning, and leakage monitoring is gradually brought into a fault prediction and health management system to realize full life cycle risk management and control. The existing leakage positioning means often face engineering problems of strong site noise, remarkable multipath propagation, limited site deployment density and energy supply communication conditions of an intelligent sensor, deviation of the installation position of the intelligent sensor and the like, so that the arrival time difference is unstable, the azimuth estimation is easy to shake or jump, the positioning result is sensitive to single observation and geometric errors, meanwhile, the fusion of the azimuth result and the geographic information of a pipeline lacks residual constraint and self-correction mechanisms, a consistency link capable of rechecking is difficult to form, high-precision positioning and reliable warning under a complex pipe section scene are influenced, and further closed-loop landing of fault prediction and health management is restricted. Disclosure of Invention The invention provides a method and a system for positioning and analyzing an infrasonic wave signal, which are used for at least solving the problems that under the conditions of multipath noise and array geometric errors, the leakage azimuth is unstable and high-precision positioning can not be realized easily. In a first aspect, the present invention provides a method for positioning and analyzing an infrasonic wave signal, including the following steps: Acquiring multichannel infrasonic wave signals distributed in an oil and gas pipeline, extracting multipath propagation characteristics, determining propagation compensation parameters, and correcting the multichannel infrasonic wave signals based on the propagation compensation parameters to obtain an infrasonic wave signal set; Calculating an arrival time difference set based on the infrasonic wave signal set, estimating an initial arrival azimuth based on the arrival time difference set, updating array geometric parameters based on the initial arrival azimuth and the arrival time difference set, estimating the arrival azimuth, and obtaining a fusion azimuth through weighted fusion; Inputting the fusion azimuth into an azimuth filter, updating noise parameters of the azimuth filter based on a historical arrival azimuth reference sequence, and outputting a refined arrival azimuth; And calculating an arrival time difference residual error based on the refined arrival azimuth, correcting the refined arrival azimuth according to the arrival time difference residual error to obtain a final arrival azimuth, determining a leakage position based on the final arrival azimuth and pipeline geographic information, and outputting a leakage positioning result. According to some embodiments of the invention, extracting multipath propagation characteristics comprises: Calculating a cross-correlation main peak width, a cross-correlation main peak to cross-correlation sub-peak energy ratio, and cross-channel coherence based on the multi-channel infrasonic signal to form a multipath propagation characteristic, determining a propagation compensation parameter includes determining an arrival time difference offset correction amount and a channel weighting coefficient based on the multipath propagation characteristic, and correcting the multi-channel infrasonic signal based on the propagation compensation parameter includes performing the arrival time difference offset correction amount and the channel weighting based on the channel weighting coefficient to obtain an infrasonic signal set. According to some embodiments of the invention, calculating the set of arrival time differences based on the set of infrasound signals comprises: And executing cross-correlation operation on any two channels of infrasonic signals in the infrasonic signal set, determining arrival time differences corresponding to the main peaks of the cross-correlation, and collecting the arrival time differences to form an arrival time difference set. According to some embodiments of the invention, estimating an initial arrival azimuth based on a set of arrival time differences comprises: and constructing a arrival azimuth candidate set based on the array geometric parameters, calculating a theoretical arrival time difference set aiming at the arrival azimuth in the arrival azimuth candidate set, and taking the arrival azimuth with the smallest di