An adaptive pulse eddy current corrosion defect detection method and detection device

Abstract

The application discloses a kind of self-adapting pulse eddy current corrosion defect detection method and detection device, belong to electromagnetic nondestructive testing technical field.The method includes: respectively collecting air domain and the pulse eddy current signal of measured piece, and difference signal is constructed after filtering and logarithmic transformation;Based on the peak time of the difference signal, the slope ratio of rising and falling section and the initial change time difference, the best fitting interval of late slope is adaptively determined;Linear fitting is carried out to the measured piece signal in the interval, and the late slope characteristic value is extracted;Using the known thickness of calibration position and corresponding slope, the residual wall thickness of measured position is calculated by formula, and the quantitative evaluation of corrosion defect is realized.The application overcomes the problem that fitting interval depends on fixed threshold in prior art, and the adaptability is poor, can automatically adapt to different detection conditions, significantly improves the stability of feature extraction and the reliability of detection result, especially suitable for online corrosion detection of ferromagnetic component with cladding layer.

Description

Technical Field

[0001] This invention relates to the field of electromagnetic nondestructive testing technology, and in particular to an adaptive pulsed eddy current corrosion defect detection method and device. Background Technology

[0002] Pulsed eddy current (PEC) testing is a typical time-domain electromagnetic nondestructive testing method. By analyzing the attenuation process of eddy currents in the test piece after rapid changes in excitation current, it enables the detection of changes in component wall thickness and corrosion defects. Compared with traditional sinusoidal excitation eddy current testing methods, PEC testing has advantages such as a wide excitation spectrum, strong penetration capability, and insensitivity to surface conditions. It can perform online testing of magnetic pressure-bearing equipment without removing the cladding or performing surface treatment, thus showing broad application prospects in the fields of petrochemical, energy, and special equipment safety testing.

[0003] Pulsed eddy current testing signals typically exhibit a nonlinear voltage response curve that decays over time. Numerous studies have shown that information closely related to material wall thickness and corrosion state is concentrated in the later stages of signal decay, particularly the near-linear decay characteristics observed in logarithmic coordinates. The mechanism lies in the fact that as material thinning occurs, its equivalent wall thickness decreases, accelerating the eddy current decay process and leading to a larger slope in the later stages. Therefore, defects are manifested as differences in electromagnetic response caused by changes in remaining wall thickness, rather than directly presenting geometric morphology. Consequently, the slope of the later-stage signal has gradually become one of the commonly used characteristic quantities in pulsed eddy current testing, capable of characterizing changes in wall thickness or the degree of corrosion thinning in ferromagnetic components.

[0004] In existing technologies, some invention patents have attempted to extract the slope as a detection feature by linearly fitting the later attenuation segment of the pulsed eddy current detection signal. For example, invention patent CN104266579B proposes a method for least-squares fitting of the signal segment of the induced voltage measurement curve after 0.1 times the characteristic attenuation time in a semi-logarithmic coordinate system, and uses the slope of the fitted line as a detection feature for measuring the wall thickness of ferromagnetic components. This method, to some extent, utilizes the linear attenuation characteristics of the later signal, providing a new approach for feature extraction of pulsed eddy current signals.

[0005] However, in practical engineering applications, this method still has certain limitations: on the one hand, the “feature decay time” it uses lacks a unified and clear basis for determination, and the selection of the fitting start point depends on empirical judgment, making it difficult to maintain consistency under different detection objects and detection conditions; on the other hand, this method does not clearly define the fitting termination interval. When the detection signal further decays and tends to flatten or is affected by noise or system drift in the later stage, continuing to participate in the fitting is prone to introducing errors, thereby affecting the stability and reliability of the slope feature.

[0006] Furthermore, existing technologies have proposed determining late-signal intervals through time-window analysis. For example, invention patent CN112946065B discloses a pulse eddy current detection method based on the slope of late-signal signals. This method performs logarithmic transformation and normalization on the detected signal, selects a fixed number of late-signal time windows based on the difference between adjacent time windows, and performs linear fitting on the signal within this interval to obtain slope characteristic values. This method, to some extent, avoids directly using the characteristic decay time as the fitting starting point.

[0007] However, the determination of the late signal interval in this type of method still relies on a pre-set fixed time window rule, and the length and position of the fitting interval are not dynamically adjusted with the overall changes in signal characteristics. When the material, wall thickness range, coating thickness, or signal noise level of the detected object changes, the fixed time window rule cannot guarantee that the fitting interval is always in the effective linear decay range of the signal, which may still lead to unstable slope feature extraction results, and to some extent depends on the selection of manual empirical parameters.

[0008] In summary, existing feature extraction methods based on the slope of the later signal in pulsed eddy current detection generally suffer from problems such as fixed rules for selecting the fitting interval, insufficient adaptability to signal changes, and strong reliance on human experience. There is a lack of a pulsed eddy current detection method that can automatically and continuously determine the fitting interval of the later slope based on the signal's own variation characteristics, while balancing fitting accuracy and stability. Therefore, it is necessary to propose a new adaptive fitting method for the later signal slope to improve the reliability and engineering applicability of feature extraction in pulsed eddy current detection. Summary of the Invention

[0009] Therefore, the purpose of this invention is to provide an adaptive pulsed eddy current corrosion defect detection method. By analyzing the variation characteristics of the pulsed eddy current detection signal and its differential signal, the starting point and ending point of the later slope fitting interval are automatically determined. Without the need to manually set the endpoints of the fitting interval, the method can achieve stable and reliable extraction of the slope feature values ​​of the later signal, thereby improving the accuracy and engineering applicability of wall thickness and corrosion defect assessment in pulsed eddy current detection.

[0010] To achieve the above objectives, the present invention provides an adaptive pulsed eddy current corrosion defect detection method, comprising the following steps:

[0011] S1. Acquire data; including acquiring the detection signal of the pulse eddy current sensor when it is placed in air without load, as the reference signal in the air domain; and acquiring the detection signal of the test device when the pulse eddy current sensor is placed above the test device.

[0012] S2. After preprocessing the air domain reference signal and the test device detection signal, differential calculation is performed. The preprocessing includes filtering and logarithmic transformation.

[0013] S3. Plot the differential signal as a curve, and mark the coordinates of the plotted curve when it reaches the maximum peak value as the starting fitting point index1.

[0014] S4. Starting from the initial fitting point index1, adaptively determine the fitting endpoint index2 based on the variation characteristics of the differential signal, and use the initial fitting point index1 and the fitting endpoint index2 as the boundary of the subsequent slope fitting interval.

[0015] S5. In the preprocessed test component detection signal, extract the signal segment within the slope fitting interval and perform linear fitting to obtain the subsequent slope feature value.

[0016] S6. Calculate the corrosion thinning amount of the tested part based on the later slope characteristic value, and determine the defect location.

[0017] More preferably, in S2, the logarithmic transformation adopts the following formula:

[0018]

[0019] in, The original signal, This is the transformed signal.

[0020] More preferably, in S2, the data of the air domain reference signal after preprocessing is represented as follows: The data after preprocessing the detection signal of the test piece is represented as The difference calculation is expressed by the following formula:

[0021] ;

[0022] in, This represents the differential signal between the air domain reference signal and the detection signal of the device under test.

[0023] More preferably, in S4, the initial fitting point index1 and the fitting endpoint index2 are used as the boundaries of the later slope fitting interval, including the following steps:

[0024] S401. Using the initial fitting point index1 as the boundary, divide the differential signal into rising segments. and falling segment signal ;

[0025] S402. Using the sliding slope method to analyze the rising segment of the differential signal. and falling segment signal Perform local least squares linear fitting separately;

[0026] S403, Calculate the rising segment signal maximum slope With falling segment signal minimum slope ;

[0027] S404, Calculate the maximum slope With minimum slope The absolute value of the ratio M;

[0028] S405, for rising segment signals Find the second derivative to obtain the rising segment signal. The coordinates at which the maximum value is obtained are denoted as the starting point of change, index. 0；

[0029] S406. Calculate the coordinate difference ΔT between the initial change point index0 and the initial fitting point index1;

[0030] S407. Based on the coordinate difference ΔT, the absolute value of the ratio M, and the initial fitting point index... 1， Calculate the fitting endpoint The coordinates are obtained; the slope fitting interval is denoted as [ , ].

[0031] More preferably, in S407, the fitting endpoint The coordinates are calculated using the following formula:

[0032] ;

[0033] in, The coordinates of the fitting endpoint are given, where ΔT is the coordinate difference and M is the absolute value of the ratio.

[0034] More preferably, in S5, from the preprocessed test component detection signal, a signal segment within the slope fitting interval is extracted, and linear fitting is performed to obtain the subsequent slope feature value, including:

[0035] Preprocessed test signal of the test piece In the middle, the slope fitting interval is truncated. , The data within [ ] were linearly fitted using the least squares method, and the slope of the fitted line was used as the slope feature value for the later stage. ;

[0036] ;

[0037] Where slope is the slope function. This is the fitted curve of the detection signal of the test piece after preprocessing.

[0038] More preferably, in S6, calculating the corrosion thinning amount of the tested part based on the later slope characteristic value and determining the defect location includes:

[0039] A calibration position is set at an uncorroded area of ​​the test piece, and the slope characteristic value of the calibration position is obtained later. ; Obtain the measured slope feature value at the detection location in the later stage. ;

[0040] Based on the known thickness at the calibration location Calculate the thickness at the measured location using the following formula. :

[0041] Calculate the thickness at the measured location. .

[0042] The present invention also provides an adaptive pulsed eddy current corrosion defect detection device for performing the above-described adaptive pulsed eddy current corrosion defect detection method, comprising:

[0043] The signal excitation and acquisition module is used to generate a pulsed excitation magnetic field and acquire the induced voltage signal.

[0044] The signal preprocessing and feature extraction module is used to process the acquired signal and extract the slope feature value in the later stage;

[0045] The thickness calculation and output module is used to calculate and output the wall thickness or corrosion assessment results based on the characteristic values.

[0046] More preferably, the signal excitation and acquisition module includes:

[0047] A pulse signal generator is used to generate pulse currents with a set frequency and duty cycle.

[0048] An excitation coil is connected to the pulse signal generator to generate an electromagnetic field, which excites the device under test to generate a pulsed eddy current field.

[0049] The receiving coil is used to sense the change in the composite magnetic field formed by the electromagnetic field generated by the excitation coil and the secondary magnetic field generated by the test device in the pulsed eddy current field, and outputs the corresponding electrical signal.

[0050] The data acquisition unit is used to synchronously acquire the voltage signal output by the receiving coil to obtain the air domain reference signal. V air With the detection signal of the test piece V test .

[0051] The adaptive pulsed eddy current corrosion defect detection method and device disclosed in this application have the following advantages compared with the prior art:

[0052] 1. It can adaptively determine the slope fitting interval in the later stage, avoiding the uncertainty caused by fixed time thresholds or preset time window rules.

[0053] This invention automatically determines the start and end points of the slope fitting interval by analyzing the variation characteristics of the pulsed eddy current detection signal and its differential signal. It does not rely on the feature decay time or a fixed number of time windows, thus overcoming the problems of fuzzy selection of fitting interval and poor versatility in the prior art.

[0054] 2. Effectively improves the stability and repeatability of slope feature extraction in the later stage.

[0055] By determining continuous fitting intervals based on signal change characteristics, linear fitting always operates on the effective attenuation phase of the signal, avoiding the introduction of noise-enhanced sections or signal-flattening sections into the fitting process, thereby reducing the impact of noise, system drift, and changes in detection conditions on the slope feature value.

[0056] 3. Reduce human experience intervention and improve the objectivity and consistency of the testing process.

[0057] This invention eliminates the need for manual setting of fitting start and end points or adjustment of fitting interval parameters. The extraction process of slope feature values ​​can be completed automatically by the program, which is beneficial for obtaining consistent detection results under different operators, different detection objects, and different detection environments.

[0058] 4. Applicable to various testing conditions, with strong engineering applicability.

[0059] Since the fitting interval is determined based on the signal's own variation characteristics, rather than depending on the specific wall thickness range of the test piece, the material's electromagnetic parameters, or the sensor's lift-off height, this invention has good adaptability under different thicknesses and different coating layers.

[0060] 5. The method has a clear structure, low computational cost, and is easy to implement in engineering and integrate into systems.

[0061] The signal processing steps used in this invention are all based on conventional digital signal processing methods. The algorithm flow is clear, the computational complexity is low, and it is suitable for implementation in embedded devices or online detection systems. This is beneficial for the automation and intelligent application of pulse eddy current detection technology. Attached Figure Description

[0062] Figure 1 This is a schematic flowchart of the adaptive pulsed eddy current corrosion defect detection method provided by the present invention.

[0063] Figure 2 This is a schematic diagram of the system apparatus for the pulsed eddy current detection method of the present invention.

[0064] Figure 3 This is a schematic diagram of the time-domain curves of the air domain reference signal and the test device detection signal after preprocessing.

[0065] Figure 4 This is a schematic diagram of the differential signal curve constructed from the air domain reference signal and the test signal of the device under test.

[0066] Figure 5 This is a schematic diagram of key feature points for determining the slope fitting interval in the later stage based on the differential signal.

[0067] Figure 6 A comparative schematic diagram showing the adaptive fitting results of the later slope obtained by the method of this invention for test pieces of different thicknesses, 9mm and 12mm. Detailed Implementation

[0068] The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0069] The following detailed descriptions are illustrative and intended to further illustrate this application. Unless otherwise specified, the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the exemplary embodiments according to this application. Exemplary embodiments according to this application will now be described in more detail with reference to the accompanying drawings. However, these exemplary embodiments may be implemented in many different forms and should not be construed as being limited to the embodiments set forth herein.

[0070] like Figure 1 As shown, the adaptive pulsed eddy current corrosion defect detection method provided by one embodiment of the present invention includes the following steps:

[0071] S1. Acquire data; including acquiring the detection signal of the pulse eddy current sensor when it is placed in air without load, as the reference signal in the air domain; and acquiring the detection signal of the test device when the pulse eddy current sensor is placed above the test device.

[0072] During the testing process, a pulsed excitation current acts on the excitation coil, generating a pulsed eddy current field in the test piece. The receiving coil senses the combined magnetic field changes formed by the electromagnetic field generated by the excitation coil and the secondary magnetic field generated by the test piece in the pulsed eddy current field, and outputs a corresponding electrical signal to form a pulsed eddy current detection signal. The test piece is set as a Q235 steel plate with a certain thickness of cladding layer on top. The excitation signal is a square wave signal with a frequency of 2Hz and a duty cycle of 50%. Figure 2 As shown.

[0073] S2. After preprocessing the air domain reference signal and the measured object detection signal, differential calculation is performed. The preprocessing includes filtering and logarithmic transformation. The sensor is placed in the air and above the measured object, respectively, to acquire the air domain signal and the measured object signal. The acquired signals are processed by a filter, and a logarithmic transformation is performed on the filtered signals to convert the original exponentially decaying time-domain signal into an approximately linear form for subsequent feature extraction and linear fitting. The transformation formula is:

[0074]

[0075] in, The original signal, This is the transformed signal.

[0076] After the above preprocessing and analysis of the original signals, the detection signals of the tested component are obtained respectively. and air domain detection signal ,like Figure 3 As shown, the logarithmically transformed detection signal and the air domain reference signal are differentially processed to obtain the differential signal. ,like Figure 4 As shown.

[0077] S3. Perform curve fitting on the differential signal, and take the coordinates of the maximum peak value in the fitted curve as the starting fitting point index1;

[0078] S4. Starting from the initial fitting point index1, adaptively determine the fitting endpoint index2 based on the variation characteristics of the difference signal, and use the initial fitting point index1 and the fitting endpoint index2 as the boundaries of the subsequent slope fitting interval; for example... Figure 5 The specific steps shown are as follows:

[0079] S401. Using the initial fitting point index1 as the boundary, divide the differential signal into rising segments. and falling segment signal ;

[0080] S402. Using the sliding slope method to analyze the rising segment of the differential signal. and falling segment signal Perform local least squares linear fitting separately;

[0081] S403, Calculate the rising segment signal maximum slope With falling segment signal minimum slope ;in:

[0082] Maximum slope of the ascending segment:

[0083] Minimum slope of the descending segment:

[0084] Where t is the sampling time, used to characterize the moment when the induced voltage signal changes with time during pulse eddy current detection, specifically the time corresponding to each sampling point after pulse excitation. This indicates the solution for the rising segment signal. The first derivative with respect to time t; This indicates the solution for the falling segment signal. The first derivative with respect to time t;

[0085] S404, Calculate the maximum slope With minimum slope The absolute value of the ratio M; ;

[0086] S405, for rising segment signals Find the second derivative to obtain the rising segment signal. The coordinates at which the maximum value is obtained are denoted as the starting point of change, index0. ;

[0087] S406. Calculate the coordinate difference ΔT between the initial change point index0 and the initial fitting point index1; ;

[0088] S407. Based on the coordinate difference ΔT, the absolute value of the ratio M, and the initial fitting point index... 1， Calculate the coordinates of the fitting endpoint The slope fitting interval is denoted as [ , ].

[0089] The above slope fitting interval is as follows: Figure 5 As shown.

[0090] S5. In the preprocessed test component detection signal, extract the signal segment within the slope fitting interval and perform linear fitting to obtain the subsequent slope feature value.

[0091] Preprocessed test signal of the test piece In the middle, the slope fitting interval is truncated. , The data within [ ] were linearly fitted using the least squares method, and the slope of the fitted line was used as the slope feature value for the later stage. ;

[0092] ;

[0093] Where slope is the slope function. This is the fitted curve of the detection signal of the test piece after preprocessing.

[0094] S6. Calculate the corrosion thinning amount of the tested part based on the later slope characteristic value, and determine the defect location.

[0095] A calibration position is set at an uncorroded area of ​​the test piece, and the slope characteristic value of the calibration position is obtained later. ; Obtain the measured slope feature value at the detection location in the later stage. ;

[0096] Based on the known thickness at the calibration location Calculate the thickness at the measured location using the following formula. : .

[0097] To verify the effectiveness of the method of this invention, pulsed eddy current testing experiments were conducted on steel plate test pieces of different thicknesses. For example, the cladding layer height was set to 5mm for a 9mm thick test piece and 30mm for a 12mm thick test piece. Detection signals were acquired under different cladding layer thickness conditions, and the slope fitting interval for the later stage was automatically determined using the method of this invention. The characteristic slope was calculated based on the fitted straight line, such as... Figure 6 As shown. Under the above conditions, the slope of the 9mm test piece was measured to be -20.12, and the slope of the 12mm test piece was measured to be -11.28.

[0098] Using a 12mm thick test piece as the calibration test piece and a 9mm thick test piece as the test piece, the characteristic value slopes of -11.28 and -20.12 were obtained from the above measurements. Based on the published formula... Calculate the thickness at the measured location. The calculated thickness of the tested component was 8.98 mm, with a relative error of 0.22%. This allows for the assessment of corrosion thinning of the component.

[0099] Obviously, the above embodiments are merely illustrative examples for clear explanation and are not intended to limit the implementation. Those skilled in the art will recognize that other variations or modifications can be made based on the above description. It is neither necessary nor possible to exhaustively list all possible implementations here. However, obvious variations or modifications derived therefrom are still within the scope of protection of this invention.

Claims

1. An adaptive pulsed eddy current corrosion defect detection method, characterized in that, Includes the following steps: S1. Acquire the collected data; This includes acquiring the detection signal of the pulsed eddy current sensor when it is placed in air without load, as a reference signal for the air domain; and acquiring the detection signal of the test device when the pulsed eddy current sensor is placed above the test device. S2. After preprocessing the air domain reference signal and the test device detection signal, differential calculation is performed. The preprocessing includes filtering and logarithmic transformation. S3. Plot the differential signal as a curve, and mark the coordinates of the plotted curve when it reaches the maximum peak value as the starting fitting point index1. S4. Starting from the initial fitting point index1, adaptively determine the fitting endpoint index2 based on the variation characteristics of the differential signal, and use the initial fitting point index1 and the fitting endpoint index2 as the boundary of the subsequent slope fitting interval. S401. Using the initial fitting point index1 as the boundary, divide the differential signal into rising segments. and falling segment signal ; S402. Using the sliding slope method to analyze the rising segment of the differential signal. and falling segment signal Perform local least squares linear fitting separately; S403, Calculate the rising segment signal maximum slope With falling segment signal minimum slope ; S404, Calculate the maximum slope With minimum slope The absolute value of the ratio M; S405, for rising segment signals Find the second derivative to obtain the rising segment signal. The coordinates at which the maximum value is obtained are denoted as the starting point of change, index. 0； S406. Calculate the coordinate difference ΔT between the initial change point index0 and the initial fitting point index1; S407. Based on the coordinate difference ΔT, the absolute value of the ratio M, and the initial fitting point index... 1， Calculate the fitting endpoint The coordinates are obtained; the slope fitting interval is denoted as [ , ]; S5. In the preprocessed test component detection signal, extract the signal segment within the slope fitting interval and perform linear fitting to obtain the subsequent slope feature value. S6. Calculate the corrosion thinning amount of the tested part based on the later slope characteristic value, and determine the defect location.

2. The adaptive pulsed eddy current corrosion defect detection method according to claim 1, characterized in that, In S2, the logarithmic transformation is performed using the following formula: in, The original signal, This is the transformed signal.

3. The adaptive pulsed eddy current corrosion defect detection method according to claim 1, characterized in that, In S2, the preprocessed data of the air domain reference signal is represented as follows: The data after preprocessing the detection signal of the test piece is represented as The difference calculation is expressed by the following formula: ； in, This represents the differential signal between the air domain reference signal and the detection signal of the device under test.

4. The adaptive pulsed eddy current corrosion defect detection method according to claim 3, characterized in that, In S407, the fitting endpoint The coordinates are calculated using the following formula: ； in, The coordinates of the fitting endpoint are given, where ΔT is the coordinate difference and M is the absolute value of the ratio.

5. The adaptive pulsed eddy current corrosion defect detection method according to claim 4, characterized in that, In S5, from the preprocessed test component detection signal, a signal segment within the slope fitting interval is extracted, and linear fitting is performed to obtain the subsequent slope feature value, including: Preprocessed test signal of the test piece In the middle, the slope fitting interval is truncated. , The data within [ ] were linearly fitted using the least squares method, and the slope of the fitted line was used as the slope feature value for the later stage. ; ； Where slope is the slope function. This is the fitted curve of the detection signal of the test piece after preprocessing.

6. The adaptive pulsed eddy current corrosion defect detection method according to claim 4, characterized in that, In S6, the corrosion thinning amount of the tested part is calculated based on the later slope characteristic value, and the defect location is determined, including: A calibration position is set at an uncorroded area of ​​the test piece, and the slope characteristic value of the calibration position is obtained later. ; Obtain the measured slope feature value at the detection location in the later stage. ; Based on the known thickness at the calibration location Calculate the thickness at the measured location using the following formula. : Calculate the thickness at the measured location. .

7. An adaptive pulsed eddy current corrosion defect detection device, characterized in that, The adaptive pulsed eddy current corrosion defect detection method for performing any one of claims 1 to 6 includes: The signal excitation and acquisition module is used to generate a pulsed excitation magnetic field and acquire the induced voltage signal. The signal preprocessing and feature extraction module is used to process the acquired signal and extract the slope feature value in the later stage; The thickness calculation and output module is used to calculate and output the wall thickness or corrosion assessment results based on the characteristic values.

8. The adaptive pulsed eddy current corrosion defect detection device according to claim 7, characterized in that, The signal excitation and acquisition module includes: A pulse signal generator is used to generate pulse currents with a set frequency and duty cycle. An excitation coil is connected to the pulse signal generator to generate an electromagnetic field, which excites the device under test to generate a pulsed eddy current field. The receiving coil is used to sense the change in the composite magnetic field formed by the electromagnetic field generated by the excitation coil and the secondary magnetic field generated by the test device in the pulsed eddy current field, and outputs the corresponding electrical signal. The data acquisition unit is used to synchronously acquire the voltage signal output by the receiving coil to obtain the air domain reference signal. V air With the detection signal of the test piece V test .

Images