Guided wave signal path optimization selection method based on signal difference coefficient and energy

Abstract

The invention relates to the technical field of ultrasonic guided wave path optimization, and particularly provides a guided wave signal path optimization selection method based on a signal difference coefficient and energy. The method comprises the following steps: constructing an ultrasonic guided wave structure defect detection system; according to the built system, sensing path grouping is carried out; the health signals and the defect signals are collected and grouped based on sensing path grouping; according to acquisition and grouping of health signals and defect signals, effective sensing path identification is carried out, invalid paths are eliminated, effective defect information is highlighted, ultrasonic guided wave defect imaging accuracy is improved, and optimization of guided wave signal paths is effectively improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of ultrasonic guided wave path, and particularly relates to a guided wave signal path optimization selection method based on signal difference coefficient and energy. BACKGROUND

[0002] In the field of structural health monitoring, pipelines, as a key component of energy transportation, chemical production and urban infrastructure, their structural integrity is directly related to the operation safety and environmental protection, so it is of great significance to effectively detect pipeline defects. At present, the common detection methods of pipeline defects include internal detector (PIG), ray detection, penetration detection and conventional ultrasonic detection. However, these methods have certain limitations: the internal detector is limited in the application of pipe diameter and has high requirements for the cleanliness of the pipeline; the ray detection has safety protection problems and is difficult to identify small defects; the penetration detection is only suitable for surface defects; and the conventional ultrasonic detection is low in efficiency and difficult to achieve wide coverage.

[0003] The ultrasonic guided wave detection technology uses the reflection, scattering or mode conversion of high-frequency mechanical waves in the pipe wall when encountering defects to identify defects, and has the advantages of long propagation distance, small attenuation, wide detection range and high sensitivity, and is especially suitable for detecting hidden or inaccessible pipe sections. At present, the pipeline defect positioning method based on ultrasonic guided wave usually arranges a sensor array on the pipe wall, collects guided wave signals under all excitation-receiving paths, and constructs a defect probability image based on the full path data. However, this method has obvious problems in practical application: among all the signal paths, a large number of paths are far away from the defect area, and the signal changes of these paths are mainly caused by environmental noise, temperature fluctuations or measurement errors, rather than real defects. Including these invalid paths (or useless paths) in the imaging calculation will introduce significant noise interference and dilute the real defect information, thereby reducing the positioning accuracy and reliability. SUMMARY

[0004] Therefore, the present application provides a guided wave signal path optimization selection method based on signal difference coefficient and energy, which can eliminate invalid paths and highlight effective defect information, and improve the accuracy of ultrasonic guided wave defect imaging, so as to effectively improve the guided wave signal path optimization.

[0005] In the first aspect, the present application provides a guided wave signal path optimization selection method based on signal difference coefficient and energy, which comprises the following steps: Step 1, building an ultrasonic guided wave structure defect detection system; Step 2, grouping the sensing paths according to the built system; Step 3, collecting and grouping the healthy signals and defect signals based on the sensing path grouping; Step 4, effective sensing path recognition is performed according to the collection and grouping of the health signals and the defect signals.

[0006] Optionally, the step 1 comprises: An ultrasonic guided wave structure defect detection system is built by using a ring-shaped piezoelectric sensing array, an ultrasonic guided wave detector, a host computer and a pipeline structure to be detected. a. The ring-shaped piezoelectric sensing array is composed of 2N piezoelectric sensors: N piezoelectric sensors are used to excite ultrasonic guided wave signals, which are defined as excitation sensors and marked as A1, A2, …, AN respectively; the other N piezoelectric sensors are used to receive ultrasonic guided wave signals, which are defined as receiving sensors and marked as R1, R2, …, RN respectively; The excitation sensors A1, A2, …, AN are pasted on the outer surface of the pipeline structure at equal intervals around the circumferential direction of the pipeline structure, that is, if the outer radius of the pipeline structure is r, the arc length between the excitation sensor Ai1 and the excitation sensor Ai1+1 is 2πr / N, i1=1, 2, …, N-1; and the plane formed by A1, A2, …, AN is perpendicular to the central axis direction of the pipeline structure, and the sensing array formed by the excitation sensors A1, A2, …, AN is called the excitation sensing array; The receiving sensors R1, R2, …, RN are pasted on the outer surface of the pipeline structure at equal intervals around the circumferential direction of the pipeline structure, that is, if the outer radius of the pipeline structure is r, the arc length between the receiving sensor Ri1 and the receiving sensor Ri1+1 is 2πr / N, i1=1, 2, …, N-1; and the plane formed by R1, R2, …, RN is perpendicular to the central axis direction of the pipeline structure, and the sensing array formed by the receiving sensors R1, R2, …, RN is called the receiving sensing array; The line between the excitation sensor Ai' and the receiving sensor Ri' of the corresponding number is parallel to the central axis of the pipeline structure, and the axial distance between the excitation sensing array and the receiving sensing array is d, i'=1, 2, …, N; b. One end of the ultrasonic guided wave detector is connected with the piezoelectric sensor, and the other end is connected with the host computer, which is used for transmitting the excitation and collected guided wave signals; c. The host computer is used for setting the excitation signal parameters and displaying and storing the collected signals; d. The pipeline structure to be detected is a circular tubular structure made of steel or aluminum alloy material.

[0007] Optionally, the step 2 comprises: The signals used are the direct signal segments collected by the receiving sensors, and when the structure is expanded, only half of the expansion plane needs to be expanded on each side; Assuming that the number of unilateral sensors is even, when the number of unilateral sensors is 8, the original plane of the pipeline expansion structure is denoted as Ω0, the 1st order plane expanded to the left is denoted as Ω-1, and the 1st order plane expanded to the right is denoted as Ω1; according to the different directions, there are totally 9 sensor path groups, and since the path of the same type has the same path length of the left-handed wave and the right-handed wave in the pipeline, only one direction is considered, i.e. there are totally 8 sensor path groups; similarly, assuming that the number of unilateral sensors is odd, when the number of unilateral sensors is 7, there are totally 7 sensor path groups; According to the different propagation directions of the paths, all the sensor paths are divided into N groups, i.e. path group 1, path group 2, …, path group N; different path groups represent different propagation directions, i.e. direction 1, direction 2, …, direction N; the number of sensor paths in the kth group is denoted as Nk, and the total number of sensor paths in the sensor array is denoted as Ntotal.

[0008] Optionally, the step 3 comprises: Step 31, health signal acquisition and grouping; Under the health state of the pipeline structure to be measured, three groups of guided wave signals are acquired by using an ultrasonic guided wave detector, denoted as health signals HS, which are HS1, HS2 and HS3 respectively, and combined with the grouping result of the sensor paths, HS1, HS2 and HS3 are divided into N signal sets, i.e. {HS11}, {HS12}, …, {HS1N}, {HS21}, {HS22}, …, {HS2N}, and {HS31}, {HS32}, …, {HS3N} respectively; the N signal sets correspond to the N path groups respectively; if the number of sampling points is denoted as M, the dimension of each signal matrix in the kth signal set is MxNk. N N N The number of sensor paths in the kth group is denoted as Nk, and the total number of sensor paths in the sensor array is denoted as Ntotal. Step 32, defect signal acquisition and grouping; Under the defect state of the pipeline structure to be measured, guided wave signals in the structure are acquired by using an ultrasonic guided wave detector, denoted as original defect signals DS; combined with the grouping result of the sensor paths, the defect signals DS are divided into N signal sets, i.e. {DS1}, {DS2}, …, {DSN}; the N signal sets correspond to the N path groups respectively; the number of signals and the data dimension in the signal set {DSk} are the same as those of the signal sets {HS1k}, {HS2k} and {HS3k}, k = 1, 2, …, N. N k k k k ​​​​​​​​​​​​

[0009] Optionally, the step 4 comprises: Step 41, effective sensing path identification based on signal difference coefficient SDC; Based on the threshold calculation of the health signal, the baseline SDC value is calculated using the health signal HS1 and the health signal HS2, and the expression is: ; Wherein, , The correlation coefficient of the two health signals is represented by The transmitter is represented by The receiver is represented by The covariance is represented by The standard deviation is represented by And The mean of the two health signals is represented by The length of the signal is represented by The sampling time is represented by And The health signals collected at different times between the transmitter and the receiver are called baseline signals. After obtaining the baseline SDC, the mean And the standard deviation The expressions are respectively: ; ; Wherein, M represents all the baseline SDCs calculated; After obtaining the mean And the standard deviation The first threshold The expression is: ; Based on the threshold calculation of the defect signal, the defect SDC value is calculated using the original defect signal DS and the health signal HS1, and the expression is: ; Wherein, , The correlation coefficient of the health signal and the defect signal is represented by The transmitter is represented by The receiver is represented by The defect signal is represented by The mean of the defect signal is represented by Compare all path defect SDC values with the first threshold When SDC k > When the signal obtained from the path differs significantly from the baseline signal obtained when the structure is healthy, it indicates that there is a high probability of a defect near the current sensing path, so as to select the path that is most affected by the defect, k=1,2,…,N; Step 42: Effective sensing path identification based on windowed scattering signal energy (WSE); Defect energy calculation uses the defect signal DS and the healthy signal HS1 to calculate the scattered signal. Based on the sensing path grouping results, the scattered signal SS1 is divided into N signal sets {SS11}, {SS12}, ..., {SS1...} according to the path direction. N The above N signal sets correspond to N path groups; then, the envelope of the scattered signal SS1 is obtained through Hilbert transform, and its expression is: ; in, Indicates signal Hilbert transform; The symbol for Hilbert transform is represented; PV represents Cauchy principal value integral. The original signal Its Hilbert transform is the real part. As the imaginary part, it constitutes a complex signal. ,Will Called The analytic signal is expressed as follows: ; Analyzing signals The modulus is the original signal envelope Its expression is: After obtaining the envelope of the scattered signal SS1, the envelope is integrated, and an integration time window is selected. In order to avoid the influence of the spiral wave, the integration start time is the arrival time of the first direct wave of the healthy signal, and the end time is the arrival time of the spiral wave. , The time window is ; Where a is the axial distance between the pipe excitation sensor array and the receiving sensor array, b is the spacing between the sensors, d1 is the distance of the straight path Am-Rm+2, and d2 is the distance of the spiral path Am-Rm+2. -1 The distance, v is the velocity of the waveform; Since the path directions and distances differ, time windows are selected for each of the N path groups of the health signal HS1, and these are applied to the N path groups corresponding to the scattered signal SS1 for integration. The energy En of each path is calculated through integration, and its expression is as follows: ; Based on the threshold value calculation of the health signal, the signals S12, S13 and S23 are obtained by making difference between each two of the three groups of health signals HS1, HS2 and HS3; then the envelopes of the signals S12, S13 and S23 after the difference are calculated, and then the integrals of the envelopes are calculated according to the steps in the above defect energy calculation, to obtain three groups of reference energies E1, E2 and E3 containing each path signal; the mean value of all the reference energies is calculated to obtain the mean value energy Eh; According to the mean value energy Eh, the second threshold value , and its expression is: ; The values of all path energies En are compared with the second threshold value , when En k > , it indicates that the path is greatly affected by the defect, and then the required path is selected, k = 1, 2, …, N; Step 43, obtaining an effective sensing path set; The effective sensing path selected based on the signal difference coefficient SDC is fused with the effective sensing path selected based on the windowed scattered signal energy WES to obtain the final optimized total selected path.

[0010] In the second aspect, the embodiments of the present application provide a computer readable storage medium, the computer readable storage medium includes a stored program, wherein when the program runs, the device where the computer readable storage medium is located is controlled to execute the guided wave signal path optimization selection method based on the signal difference coefficient and the energy in the first aspect or any possible implementation manner of the first aspect.

[0011] In the third aspect, the embodiments of the present application provide an electronic device, including: one or more processors; a memory; and one or more computer programs, wherein the one or more computer programs are stored in the memory, and the one or more computer programs include instructions, when the instructions are executed by the device, the device is caused to execute the guided wave signal path optimization selection method based on the signal difference coefficient and the energy in the first aspect or any possible implementation manner of the first aspect.

[0012] In the technical scheme provided by the present application, the method includes building an ultrasonic guided wave structure defect detection system; according to the built system, the sensing path grouping is performed; based on the sensing path grouping, the health signal and the defect signal are collected and grouped; according to the collection and grouping of the health signal and the defect signal, the effective sensing path is identified, the invalid path is eliminated, the effective defect information is highlighted, the accuracy of the ultrasonic guided wave defect imaging is improved, and the optimization of the guided wave signal path is effectively improved. BRIEF DESCRIPTION OF DRAWINGS

[0013] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the embodiments will be briefly introduced as follows. Obviously, the drawings in the following description only constitute some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0014] Figure 1 The flow chart of the guided wave signal path optimization selection method based on signal difference coefficient and energy provided by the embodiments of the present application is shown in the figure. Figure 2 The pipeline unfolding schematic diagram of 16 sensors provided by the embodiments of the present application is shown in the figure. Figure 3 The pipeline unfolding schematic diagram of 14 sensors provided by the embodiments of the present application is shown in the figure. Figure 4 The Am-Rm+2 path schematic diagram provided by the embodiments of the present application is shown in the figure. Figure 5 The schematic diagram of an electronic device provided by the embodiments of the present application is shown in the figure. DETAILED DESCRIPTION

[0015] In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely in combination with the drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, but not all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without any creative effort are within the scope of protection of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without any creative effort are within the scope of protection of the present application.

[0016] The terms used in the embodiments of the present application are only for the purpose of describing specific embodiments, and are not intended to limit the present application. The singular forms "a", "an" and "the" used in the embodiments of the present application are also intended to include the plural forms, unless the context clearly indicates otherwise.

[0017] It should be understood that the term "and / or" used herein is only to describe the association relationship of the associated objects, which means that there can be three relationships, for example, A and / or B, which can represent the three cases of A alone, A and B together, and B alone. In addition, the character " / " in this paper generally represents that the front and rear associated objects are a "or" relationship.

[0018] Depending on context, the word "if" as used herein can be interpreted to mean "when" or "while" or "in response to determining" or "in response to detecting." Similarly, the phrase "if it is determined" or "if [a stated condition or event] is detected" can be interpreted to mean "when it is determined" or "in response to determining" or "when [the stated condition or event] is detected" or "in response to detecting [the stated condition or event]."

[0019] Figure 1 A flowchart of a method for guided wave signal path optimization selection based on signal difference coefficient and energy according to an embodiment of the present application is shown in FIG. 1, which comprises the following steps: Figure 1 Step 1, building an ultrasonic guided wave structure defect detection system.

[0020] In the embodiment of the present application, step 1 comprises: building an ultrasonic guided wave structure defect detection system by using a ring-shaped piezoelectric sensing array, an ultrasonic guided wave detector, a host computer and a pipeline structure to be detected; a. the ring-shaped piezoelectric sensing array is composed of 2N piezoelectric sensors: N piezoelectric sensors are used to excite ultrasonic guided wave signals, which are defined as excitation sensors and marked as A1, A2, …, AN respectively; the other N piezoelectric sensors are used to receive ultrasonic guided wave signals, which are defined as receiving sensors and marked as R1, R2, …, RN respectively; the excitation sensors A1, A2, …, AN are pasted on the outer surface of the pipeline structure at equal intervals around the circumferential direction, i.e. if the outer radius of the pipeline structure is r, the arc length between the excitation sensor Ai1 and the excitation sensor Ai1+1 is 2πr / N, i1=1, 2, …, N-1; and the plane composed of A1, A2, …, AN is perpendicular to the central axis direction of the pipeline structure, which is called the excitation sensor array composed of A1, A2, …, AN; the receiving sensors R1, R2, …, RN are pasted on the outer surface of the pipeline structure at equal intervals around the circumferential direction, i.e. if the outer radius of the pipeline structure is r, the arc length between the receiving sensor Ri1 and the receiving sensor Ri1+1 is 2πr / N, i1=1, 2, …, N-1; and the plane composed of R1, R2, …, RN is perpendicular to the central axis direction of the pipeline structure, which is called the receiving sensor array composed of R1, R2, …, RN; the line between the excitation sensor Ai' and the receiving sensor Ri' of the corresponding number is parallel to the central axis of the pipeline structure, and the axial distance between the excitation sensor array and the receiving sensor array is d, i'=1, 2, …, N; ​b, one end of the ultrasonic guided wave detector is connected with the piezoelectric sensor, and the other end is connected with the upper computer, and is used for transmitting the excited and collected guided wave signals, and has the functions of filtering and amplification; c, the upper computer is used for setting the excitation signal parameters and displaying and storing the collected signals; d, the pipeline structure to be detected is a circular tubular structure made of steel or aluminum alloy material.

[0021] Step 2, according to the system built, the sensing path grouping is carried out.

[0022] In the research of spiral guided wave propagation path, the pipeline structure is usually unfolded into a plane. The analysis of high-order spiral guided wave propagation path is realized by copying an infinite number of extension planes on the left and right sides of the pipeline unfolding plane. This method converts the three-dimensional pipeline problem into a two-dimensional plate problem.

[0023] In the embodiment of the application, step 2 comprises: The signal adopted is the direct signal segment collected by the receiving sensor. When the structure is expanded, only half of the unfolding plane needs to be expanded on each side; Suppose that the unilateral sensor is even, such as Figure 2 As shown in the figure, when the unilateral sensor is 8, as shown by the green area, the original plane of the pipeline unfolding structure is marked as Ω0, wherein the excitation sensors are marked as A1, A2, …, A8; the receiving sensors are marked as R1, R2, …, R8; the 1st plane extended to the left is marked as Ω-1, wherein the excitation sensors are marked as A1 -1 , A2 -1 , …, A8 -1 ; the receiving sensors are marked as R1 -1 , R2 -1 , …, R8 -1 , and the 1st plane extended to the right is marked as Ω1, wherein the excitation sensors are marked as A11, A21, …, A81; the receiving sensors are marked as R11, R21, …, R81, as shown by the gray area; the green dashed area represents the complete area adopted, which can be seen from Figure 2 According to the different directions, there are a total of 9 sensor path groups. Since the type path has the same length of left-handed wave and right-handed wave in the pipeline, only one direction is considered, that is, a total of 8 sensor path groups; similarly, suppose that the unilateral sensor is odd, such as Figure 3 As shown in the figure, when the unilateral sensor is 7, as shown by the green area, the original plane of the pipeline unfolding structure is marked as Ω0, wherein the excitation sensors are marked as A1, A2, …, A7; the receiving sensors are marked as R1, R2, …, R7; the 1st plane extended to the left is marked as Ω-1, wherein the excitation sensors are marked as A1 -1 , A2 -1 , …, A7-1 ; receiving sensors are denoted as R1 -1 , R2 -1 , …, R7 -1 , the first-order plane expanding to the right is denoted as Ω1, wherein the excitation sensors are denoted as A11, A21, …, A71; the receiving sensors are denoted as R11, R21, …, R71, as shown in the gray area; the green dashed line area represents the complete area adopted, which has a total of 7 sensor path groups; According to the different propagation directions of the paths, all the sensing paths are divided into N groups, i.e., path group 1, path group 2, …, path group N; different path groups represent different propagation directions, i.e., direction 1, direction 2, …, direction N; the number of sensing paths in the kth group is denoted as , and the total number of sensing paths in the sensing array is .

[0024] Step 3, based on the grouping of the sensing paths, the healthy signals and the defect signals are collected and grouped.

[0025] In the embodiment of the present application, step 3 comprises: Step 31, healthy signal collection and grouping; Under the health state of the to-be-tested pipeline structure, three groups of guided wave signals are collected by using an ultrasonic guided wave detector, denoted as healthy signals (HS), which are HS1, HS2 and HS3 respectively, and combined with the grouping results of the sensing paths, HS1, HS2 and HS3 are divided into N signal sets, i.e., {HS11}, {HS12}, …, {HS1 N}; {HS21}, {HS22}, …, {HS2 N}; {HS31}, {HS32}, …, {HS3 N}; the N signal sets correspond to N path groups respectively; if the number of sampling points is denoted as , the number of sensing paths in the kth group is denoted as , then the dimension of each signal matrix in the kth signal set is . Step 32, defect signal collection and grouping; Under the defect state of the to-be-tested pipeline structure, guided wave signals in the structure are collected by using an ultrasonic guided wave detector, denoted as original defect signals (DS); combined with the grouping results of the sensing paths, the defect signals DS are divided into N signal sets, i.e., {DS1}, {DS2}, …, {DS N}, the above N signal sets correspond to N path groups respectively; the number of signals and the data dimension in the signal set {DS k} are the same as those in the signal set {HS1 k}.k}、{HS3 k} the number of signals and data dimensions are the same, k = 1, 2, …, N.

[0026] Step 4, according to the collection and grouping of health signals and defect signals, effective sensing path recognition is performed.

[0027] In the embodiment of the application, step 4 comprises: Step 41, effective sensing path recognition based on signal difference coefficient (SDC); Based on the threshold calculation of the health signal, the baseline SDC value is calculated by using the health signal HS1 and the health signal HS2, and the expression is: ; Wherein, , denotes the correlation coefficient of two health signals, denotes the transmitter, denotes the receiver; denotes the covariance; denotes the standard deviation; and denotes the mean of two health signals; denotes the length of the signal; denotes the sampling time; and denotes the health signal collected at different times between the transmitter and the receiver, which is called the baseline signal; After obtaining the baseline SDC, the mean and the standard deviation are calculated, and the expressions are respectively: ; ; Wherein, M denotes all the calculated baseline SDCs; After obtaining the mean and the standard deviation , the first threshold is calculated, and the expression is: ; Based on the threshold calculation of the defect signal, the defect SDC value is calculated by using the original defect signal DS and the health signal HS1, and the expression is: ; Wherein, , denotes the correlation coefficient of the health signal and the defect signal, denotes the transmitter, Indicates the receiver; Indicates a defect signal; This represents the mean of the defect signal; Compare all path defect SDC values ​​with the first threshold When comparing, when SDC k > When the signal obtained from the path differs significantly from the baseline signal obtained when the structure is healthy, it indicates that there is a high probability of a defect near the current sensing path, so as to select the path that is most affected by the defect, k=1,2,…,N; Step 42: Effective sensing path identification based on windowed scattered signal energy (WSE); Defect energy calculation uses the defect signal DS and the healthy signal HS1 to calculate the scattered signal. Based on the sensing path grouping results, the scattering signal (SS1) is divided into N signal sets {SS11}, {SS12}, ..., {SS1...} according to the path direction. N The above N signal sets correspond to N path groups; then, the envelope of the scattered signal SS1 is obtained through Hilbert transform, and its expression is: ; in, Indicates signal Hilbert transform; The symbol for Hilbert transform is represented; PV represents Cauchy principal value integral. The original signal Its Hilbert transform is the real part. As the imaginary part, it constitutes a complex signal. ,Will Called The analytic signal is expressed as follows: ; Analyzing signals The magnitude (amplitude) is the original signal. envelope Its expression is: After obtaining the envelope of the scattered signal SS1, the envelope is integrated, and an integration time window is selected. In order to avoid the influence of the spiral wave, the integration start time is the arrival time of the first direct wave of the healthy signal, and the end time is the arrival time of the spiral wave. In embodiments of the present invention, such as Figure 4As shown, the green area represents the original plane of the pipe development structure. The 0th-order plane is denoted as Ω0, and Am-Rm+2 in the 0th-order plane is a straight path. The excitation sensors are denoted as A1, ..., Am, Am+1, Am+2, ..., AN; the receiving sensors are denoted as R1, ..., Rm, Rm+1, Rm+2, ..., RN. The 1st-order plane extending to the left is denoted as Ω-1, and the 1st-order plane Am-Rm+2... -1 The path is a spiral, where the excitation sensor is denoted as A1. -1 ..., Am+1 -1 Am+2 -1 , ..., AN -1 The receiving sensor is denoted as R1. -1 ..., Rm+1 -1 Rm+2 -1 ..., RN -1 The first-order plane extending to the right is denoted as Ω1, where the excitation sensors are denoted as A11, ..., Am+11, Am+21, ..., AN1; and the receiving sensors are denoted as R11, ..., Rm+11, Rm+21, ..., RN1, as shown in the gray area. , The time window is ; Where a is the axial distance between the pipe excitation sensor array and the receiving sensor array, b is the spacing between the sensors, d1 is the distance of the straight path Am-Rm+2, and d2 is the distance of the spiral path Am-Rm+2. -1 The distance, v is the velocity of the waveform; Since the path directions and distances differ, time windows are selected for each of the N path groups of the health signal HS1, and these are applied to the N path groups corresponding to the scattered signal SS1 for integration. The energy En of each path is calculated through integration, and its expression is as follows: ; Based on the threshold calculation of the health signals, signals S12, S13, and S23 are obtained by subtracting each pair of the three sets of health signals HS1, HS2, and HS3. Then, the envelopes of the subtracted signals S12, S13, and S23 are calculated respectively. Then, the envelopes are integrated according to the steps in the above defect energy calculation to obtain three sets of reference energies E1, E2, and E3 containing the signals of each path. The mean energy Eh is obtained by averaging all the reference energies. Calculate the second threshold based on the mean energy Eh. Its expression is: ; Compare all path energy En values ​​with the second threshold Comparison, when En k > When the path is affected by the defect, the required path is selected, k=1, 2, …, N; Step 43, obtaining an effective sensing path set; The effective sensing path selected based on the signal difference coefficient SDC is fused with the effective sensing path selected based on the windowed scattering signal energy WES to obtain a final optimized total selected path.

[0028] Compared with the prior art, the present application has the following beneficial effects: (1) The method greatly improves the analysis efficiency: the core advantage is to significantly reduce the computational complexity and time cost. It quickly filters out a few key paths that are most affected by the defect from all sensor paths through simple signal difference coefficient and envelope area calculation, thereby avoiding complex full waveform analysis on all paths. This realizes the change from comprehensive net casting to key fishing, making the defect positioning process more efficient and fast, especially suitable for real-time or near real-time monitoring engineering applications.

[0029] (2) The method effectively enhances the positioning reliability: it can actively focus on the sensitive signal of the defect, thereby improving the accuracy and anti-interference ability of positioning. Since the defect has the greatest impact on the path signal passing through its vicinity, the high-value path selected itself clearly indicates the possible area of the defect. This not only narrows the search range of the fine algorithm and improves the positioning accuracy, but also performs a signal filtering, effectively suppressing the interference of environmental noise and other irrelevant factors, making the final positioning result clearer and more reliable.

[0030] The technical scheme provided by the present application, the method comprises building an ultrasonic guided wave structure defect detection system; according to the built system, the sensing path grouping is carried out; based on the sensing path grouping, the health signal and the defect signal are collected and grouped; according to the collection and grouping of the health signal and the defect signal, the effective sensing path identification is carried out, which eliminates the invalid path, highlights the effective defect information, improves the ultrasonic guided wave defect imaging accuracy, and effectively improves the optimization of the guided wave signal path.

[0031] The various steps of the embodiments of the present application can be executed by an electronic device. The electronic device includes but is not limited to a tablet computer, a portable PC, a desktop computer, etc.

[0032] The embodiment of the present application provides a computer readable storage medium, which comprises a stored program, wherein when the program runs, the electronic device where the computer readable storage medium is located executes the embodiment of the guided wave signal path optimization selection method based on the signal difference coefficient and energy.

[0033] Figure 5A schematic diagram of an electronic device is provided for an embodiment of the present application, as shown in Figure 5 The electronic device 21 includes a processor 211, a memory 212, and a computer program 213 stored in the memory 212 and executable on the processor 211, which, when executed by the processor 211, implements the method for optimizing and selecting a guided wave signal path based on a signal difference coefficient and energy in the embodiment. To avoid repetition, details are not described here.

[0034] The electronic device 21 includes, but is not limited to, the processor 211 and the memory 212. Those skilled in the art can understand that Figure 5 The electronic device 21 is only an example and does not constitute a limitation on the electronic device 21, and can include more or fewer components than shown, or combine certain components, or different components, for example, the electronic device can also include an input / output device, a network access device, a bus, etc.

[0035] The processor 211 can be a central processing unit (CPU), and can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor.

[0036] The memory 212 can be an internal storage unit of the electronic device 21, such as a hard disk or a memory of the electronic device 21. The memory 212 can also be an external storage device of the electronic device 21, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. Further, the memory 212 can include both the internal storage unit and the external storage device of the electronic device 21. The memory 212 is used to store computer programs and other programs and data required by the network device. The memory 212 can also be used to temporarily store data that has been output or will be output.

[0037] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the system, device and unit described above can refer to the corresponding process in the foregoing method embodiments, which will not be described here.

[0038] The above description is only the preferred embodiment of the present application, and is not used to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the scope of protection of the present application.

Claims

1. A method for optimizing the selection of guided signal paths based on signal difference coefficient and energy, characterized in that, The method includes: Step 1: Build an ultrasonic guided wave structure defect detection system; Step 2: Group the sensor paths according to the system being built; Step 3: Based on sensor path grouping, collect and group health signals and defect signals; Step 4: Based on the collection and grouping of health signals and defect signals, identify effective sensing paths.

2. The method according to claim 1, characterized in that, Step 1 includes: An ultrasonic guided wave structural defect detection system was built using a ring piezoelectric sensor array, an ultrasonic guided wave detector, a host computer, and the pipe structure under test. a. The ring piezoelectric sensing array consists of 2N piezoelectric sensors: N piezoelectric sensors are used to excite ultrasonic guided wave signals, and these N piezoelectric sensors are defined as excitation sensors, labeled A1, A2, ..., AN respectively; the other N piezoelectric sensors are used to receive ultrasonic guided wave signals, and these N piezoelectric sensors are defined as receiving sensors, labeled R1, R2, ..., RN respectively. Excitation sensors A1, A2, ..., AN are attached at equal intervals to the outer surface of the pipe structure along the circumference. That is, if the outer radius of the pipe structure is r, the arc length between excitation sensors Ai1 and Ai1+1 is 2πr / N, i1=1,2,...,N-1; and the plane formed by A1, A2, ..., AN is perpendicular to the central axis of the pipe structure. The sensor array formed by excitation sensors A1, A2, ..., AN is called the excitation sensor array. Receiving sensors R1, R2, ..., RN are attached at equal intervals around the circumference of the pipe structure on the outer surface of the pipe. That is, if the outer radius of the pipe structure is r, the arc length between receiving sensors Ri1 and Ri1+1 is 2πr / N, i1=1,2,...,N-1; and the plane formed by R1, R2, ..., RN is perpendicular to the central axis of the pipe structure. The sensing array formed by receiving sensors R1, R2, ..., RN is called the receiving sensor array. The line connecting the excitation sensor Ai' and the corresponding numbered receiving sensor Ri' is parallel to the central axis of the pipe structure. The axial distance between the excitation sensor array and the receiving sensor array is d, and i' = 1, 2, ..., N. b. One end of the ultrasonic guided wave detector is connected to the piezoelectric sensor, and the other end is connected to the host computer to transmit the guided wave signals for excitation and acquisition. c. The host computer is used for setting excitation signal parameters and displaying and storing acquired signals; d. The pipe structure to be tested is a circular tube made of steel or aluminum alloy.

3. The method according to claim 2, characterized in that, Step 2 includes: The signal used is the direct signal segment collected by the receiving sensor. When expanding the structure, only half of the unfolded plane needs to be expanded on each side. Assuming there are an even number of sensors on one side, when there are 8 sensors on one side, the original plane of the pipe development structure, the 0th order plane, is denoted as Ω0; the 1st order plane extending to the left is denoted as Ω-1, and the 1st order plane extending to the right is denoted as Ω1; depending on the direction, there are a total of 9 sensor path groups. For a path of this type, the left-handed and right-handed wave paths in the pipe have the same length. Only one direction is considered, meaning there are a total of 8 sensor path groups. Similarly, assuming an odd number of sensors on one side, when there are 7 sensors on one side, there are a total of 7 sensor path groups. Based on the different propagation directions of the paths, all sensing paths are divided into N groups, namely path group 1, path group 2, ..., path group N; different path groups represent different propagation directions, namely direction 1, direction 2, ..., direction N; using Let represent the number of sensor paths within the k-th group. This represents the number of all sensing paths within the sensor array.

4. The method according to claim 3, characterized in that, Step 3 includes: Step 31: Health signal acquisition and grouping; Under the condition of a healthy pipeline structure, three sets of guided wave signals, denoted as health signals HS1, HS2, and HS3, are collected using an ultrasonic guided wave detector. Based on the sensor path grouping results, HS1, HS2, and HS3 are further divided into N signal sets, namely {HS11}, {HS12}, ..., {HS1...}. N };{HS21},{HS22},…,{HS2 N }; {HS31}, {HS32},…, {HS3 N }; N signal sets correspond to N path groups; if the number of sampling points is used express, Let represent the number of sensing paths within the k-th group. Then, the dimension of each signal matrix in the k-th signal set is... ; Step 32: Defect signal acquisition and grouping; Under the condition of defects in the pipeline structure, guided wave signals in the structure are collected using an ultrasonic guided wave detector, denoted as the original defect signal DS. Based on the sensor path grouping results, the defect signal DS is divided into N signal sets, namely {DS1}, {DS2}, ..., {DS...}. N The above N signal sets correspond to N path groups; signal set {DS} k The number of signals and data dimensions in} are related to the signal set {HS1}. k }、{HS2 k }、{HS3 k The number of signals and the data dimension are the same, k=1,2,…,N.

5. The method according to claim 4, characterized in that, Step 4 includes: Step 41: Effective sensing path identification based on signal difference coefficient SDC; Threshold calculation based on health signals: The baseline SDC value is calculated using health signal HS1 and health signal HS2, and its expression is as follows: ； in, , The correlation coefficient represents the correlation between two health signals. Indicates the transmitter. Indicates the receiver; Represents covariance; Indicates standard deviation; and This represents the mean of two health signals; Indicates the length of the signal; Indicates the sampling time; and The health signal collected at different times between the transmitter and receiver is called the baseline signal; After obtaining the baseline SDC, calculate its mean. and standard deviation Their expressions are as follows: ； ； Where M represents all calculated baseline SDCs; Get the mean and standard deviation Then, calculate the first threshold. Its expression is: ； Threshold calculation based on defect signals uses the original defect signal DS and the healthy signal HS1 to calculate the defect SDC value, and its expression is as follows: ； in, , The correlation coefficient represents the correlation between health signals and defect signals. Indicates the transmitter. Indicates the receiver; Indicates a defect signal; This represents the mean of the defect signal; Compare all path defect SDC values ​​with the first threshold When comparing, when SDC k > When the signal obtained from the path differs significantly from the baseline signal obtained when the structure is healthy, it indicates that there is a high probability of a defect near the current sensing path, so as to select the path that is most affected by the defect, k=1,2,…,N; Step 42: Effective sensing path identification based on windowed scattering signal energy (WSE); Defect energy calculation uses the defect signal DS and the healthy signal HS1 to calculate the scattered signal. Based on the sensing path grouping results, the scattered signal SS1 is divided into N signal sets {SS11}, {SS12}, ..., {SS1...} according to the path direction. N The above N signal sets correspond to N path groups; then, the envelope of the scattered signal SS1 is obtained through Hilbert transform, and its expression is: ； in, Indicates signal Hilbert transform; The symbol for Hilbert transform is represented; PV represents Cauchy principal value integral. The original signal Its Hilbert transform is the real part. As the imaginary part, it constitutes a complex signal. ,Will Called The analytic signal is expressed as follows: ； Analyzing signals The modulus is the original signal envelope Its expression is: After obtaining the envelope of the scattered signal SS1, the envelope is integrated, and an integration time window is selected. In order to avoid the influence of the spiral wave, the integration start time is the arrival time of the first direct wave of the healthy signal, and the end time is the arrival time of the spiral wave. , The time window is ; Where a is the axial distance between the pipe excitation sensor array and the receiving sensor array, b is the spacing between the sensors, d1 is the distance of the straight path Am-Rm+2, and d2 is the distance of the spiral path Am-Rm+2. -1 The distance, v is the velocity of the waveform; Since the path directions and distances differ, time windows are selected for each of the N path groups of the health signal HS1, and these are applied to the N path groups corresponding to the scattered signal SS1 for integration. The energy En of each path is calculated through integration, and its expression is as follows: ； Based on the threshold calculation of the health signals, signals S12, S13, and S23 are obtained by subtracting each pair of the three sets of health signals HS1, HS2, and HS3. Then, the envelopes of the subtracted signals S12, S13, and S23 are calculated respectively. Then, the envelopes are integrated according to the steps in the above defect energy calculation to obtain three sets of reference energies E1, E2, and E3 containing the signals of each path. The mean energy Eh is obtained by averaging all the reference energies. Calculate the second threshold based on the mean energy Eh. Its expression is: ； Compare all path energy En values ​​with the second threshold Comparison, when En k > If the path is significantly affected by the defect, then the required path is selected, k=1, 2, ..., N; Step 43: Obtain the set of valid sensor paths; The effective sensing path selected based on the signal difference coefficient (SDC) is fused with the effective sensing path selected based on the windowed scattering signal energy (WES) to obtain the final optimized overall selection path.

6. A computer-readable storage medium, characterized in that, The computer-readable storage medium includes a stored program, wherein, when the program is executed, it controls the device containing the computer-readable storage medium to perform the guided wave signal path optimization selection method based on signal difference coefficient and energy as described in any one of claims 1 to 5.

7. An electronic device, characterized in that, include: One or more processors; Memory; And one or more computer programs, wherein the one or more computer programs are stored in the memory, the one or more computer programs including instructions that, when executed by the device, cause the device to perform the guided wave signal path optimization selection method based on signal difference coefficient and energy as described in any one of claims 1 to 5.

Images