Optical fiber coupling medium quality factor inversion method for hole sealing quality detection

Abstract

The invention relates to an optical fiber coupling medium quality factor inversion method for hole sealing quality detection, and the method comprises the steps: carrying out the vibration signal direct wave first arrival pickup of optical fiber monitoring data, and carrying out the interception of the optical fiber monitoring data based on the direct wave first arrival pickup, and obtaining a vibration signal; generating a joint loss function based on an amplitude Q value constructor function, a phase constructor function and an envelope constructor function; and carrying out inversion on the joint loss function by adopting an L-BFGS-B optimization inversion algorithm to obtain a Q value of amplitude, phase and envelope joint constraint. By fusing the distributed optical fiber sensing technology and the multi-parameter Q value inversion algorithm, important breakthroughs of technical performance and engineering applicability are realized in the field of hole sealing quality detection. Compared with a traditional method, the method has the advantages that full-well-section, real-time and dynamic quantitative diagnosis can be carried out on the coupling state of the shaft hole sealing stratum through the continuous space sensing capacity of the distributed optical fiber and the high sensitivity characteristic of the Q value to the medium attenuation characteristic.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of hole sealing quality detection, in particular to a fiber-coupled medium quality factor inversion method for hole sealing quality detection. BACKGROUND

[0002] With the rapid development of underground engineering construction, drilling and grouting technology is widely used in tunnel surrounding rock reinforcement, mine curtain water plugging and foundation treatment fields. The grouting quality directly affects the safety of the project. Traditional monitoring methods such as resistivity method, acoustic detection or drilling core have defects such as low spatial resolution, poor real-time performance and high cost. SUMMARY

[0003] In order to overcome at least one of the deficiencies in the prior art, the present application provides a fiber-coupled medium quality factor inversion method for hole sealing quality detection.

[0004] In a first aspect, a fiber-coupled medium quality factor inversion method for hole sealing quality detection is provided, comprising:

[0005] The fiber monitoring data is subjected to direct wave first arrival picking of the vibration signal, and the fiber monitoring data is subjected to cutting based on the direct wave first arrival picking, to obtain the vibration signal;

[0006] The vibration signal is subjected to amplitude compensation based on the spherical diffusion compensation principle, to obtain the compensated vibration signal;

[0007] The compensated vibration signal is subjected to band-pass filtering, to obtain the filtered vibration signal;

[0008] The filtered vibration signal is cut according to the calculation window length, to obtain the cut signal; the cut signal is subjected to Hilbert transform, to generate the analytic signal; and the instantaneous amplitude, instantaneous phase and spectrum are calculated according to the analytic signal;

[0009] A Q value construction function based on amplitude is generated according to the instantaneous amplitude and amplitude attenuation model; a construction function based on phase is generated according to the instantaneous phase; an envelope-based construction function is generated according to the relative delay broadening effect of the envelope peak position; a joint loss function is generated according to the Q value construction function based on amplitude, the construction function based on phase and the envelope-based construction function; and the Q value is the fiber-coupled medium quality factor;

[0010] The joint loss function is inversed by using L-BFGS-B optimization inversion algorithm, to obtain the Q value jointly constrained by amplitude, phase and envelope.

[0011] In one embodiment, the fiber monitoring data is subjected to direct wave first arrival picking of the vibration signal, and the fiber monitoring data is subjected to cutting based on the direct wave first arrival picking, to obtain the vibration signal, comprising:

[0012] The fiber monitoring data comprises data of multiple channels;

[0013] The data of each channel is subjected to a direct wave first arrival picking of a vibration signal to obtain a direct wave first arrival of each channel.

[0014] An average value of the direct wave first arrivals of all channels is taken as a signal intercepting marker time point.

[0015] A time period before and after the signal intercepting marker time point is set as an intercepting time window, and the data of each channel is subjected to intercepting to obtain a vibration signal of each channel.

[0016] In one embodiment, the vibration signal is subjected to amplitude compensation based on a spherical divergence compensation principle to obtain a compensated vibration signal, comprising:

[0017] The vibration signal comprises vibration signals of multiple channels.

[0018] An amplitude compensation coefficient of each channel is calculated:

[0019]

[0020] wherein G(i) is the amplitude compensation coefficient of the i th channel, r i is the offset distance of the i th channel.

[0021] The vibration signal of each channel is subjected to amplitude compensation using the corresponding amplitude compensation coefficient to obtain a compensated vibration signal of each channel.

[0022] In one embodiment, a Q value construction function based on amplitude is generated according to an instantaneous amplitude and an amplitude decay model, comprising:

[0023] The instantaneous amplitude comprises instantaneous amplitudes of multiple channels, and for each channel, the instantaneous amplitudes at all time points within a time window are calculated to form an amplitude spectrum of the channel.

[0024] Any channel is taken as an observation channel; according to the amplitude decay model, a theoretical decayed reference amplitude spectrum of the observation channel is determined, and the following formula is used:

[0025]

[0026] wherein A(f) is the theoretical decayed reference amplitude spectrum of the observation channel, A0(f) is the amplitude spectrum of the reference channel, f is the frequency, Δt is the time difference of wave propagation from the reference channel to the observation channel, and Q is the quality factor of the fiber-coupled medium.

[0027] The mean square error of the amplitude spectrum of the observation channel and the theoretical decayed reference amplitude spectrum A(f) of the observation channel is calculated as the Q value construction function L amp based on amplitude of the observation channel.

[0028] In an embodiment, a phase-based constructor is generated according to the instantaneous phase, comprising:

[0029] The instantaneous phase comprises the instantaneous phase of each trace;

[0030] Any trace is regarded as an observation trace; the mean square error between the instantaneous phase of the observation trace and the instantaneous phase of the reference trace is calculated as the phase-based constructor L phase .

[0031] In an embodiment, an envelope-based constructor is generated according to the relative delay broadening effect of the spectrum and the envelope peak position, comprising:

[0032] The spectrum comprises the spectrum of each trace;

[0033] For each trace, the signal envelope is determined according to the spectrum, and the peak time of the signal envelope is determined;

[0034] Any trace is regarded as an observation trace; the ratio of the peak time of the signal envelope of the observation trace to the peak time of the signal envelope of the reference trace is regarded as the observation trace ratio;

[0035] When Q→∞, the theoretical ratio is 1, and Q is the fiber coupling medium quality factor;

[0036] The mean square error between the observation trace ratio and the theoretical ratio is calculated as the envelope-based constructor L env .

[0037] In an embodiment, the joint loss function is:

[0038] L Q =ω1L amp +ω2L phase +ω3L env

[0039] Wherein, L Q is the joint loss function, L amp is the amplitude-based Q value constructor, L phase is the phase-based constructor, L env is the envelope-based constructor, ω1+ω2+ω3=1, ω1 is the weight of L amp , ω2 is the weight of L phase , and ω3 is the weight of L env .

[0040] In an embodiment, the method further comprises:

[0041] Based on the fiber monitoring data at different monitoring times, the amplitude, phase, and envelope joint-constrained Q value at different monitoring times are obtained.

[0042] In a second aspect, a fiber-coupled medium quality factor inversion device for borehole sealing quality detection is provided, comprising:

[0043] A first arrival picking and intercepting module picks up a direct wave first arrival of a vibration signal from fiber monitoring data, and intercepts the fiber monitoring data based on the direct wave first arrival picking to obtain the vibration signal;

[0044] An amplitude compensation module compensates the amplitude of the vibration signal based on a spherical diffusion compensation principle to obtain a compensated vibration signal;

[0045] A band-pass filtering module is configured to perform band-pass filtering on the compensated vibration signal to obtain a filtered vibration signal;

[0046] A Hilbert transform module is configured to intercept the filtered vibration signal according to a calculation time window length to obtain an intercepted signal, perform Hilbert transform on the intercepted signal to generate an analytic signal, and calculate an instantaneous amplitude, an instantaneous phase, and a spectrum from the analytic signal;

[0047] A joint loss function generation module is configured to generate an amplitude-based Q value construction function from the instantaneous amplitude and an amplitude attenuation model, generate a phase-based construction function from the instantaneous phase, generate an envelope-based construction function from a relative delay broadening effect of the spectrum and an envelope peak position, and generate a joint loss function from the amplitude-based Q value construction function, the phase-based construction function, and the envelope-based construction function. The Q value is a fiber-coupled medium quality factor;

[0048] An inversion module is configured to perform inversion on the joint loss function using an L-BFGS-B optimization inversion algorithm to obtain a Q value jointly constrained by the amplitude, the phase, and the envelope.

[0049] In a third aspect, a computer-readable storage medium is provided, which stores a computer program. When the computer program is executed by a processor, the fiber-coupled medium quality factor inversion method for borehole sealing quality detection described above is implemented.

[0050] Compared with the prior art, the present application has the following beneficial effects: by fusing distributed optical fiber sensing technology and multi-parameter Q value inversion algorithm, the present application realizes important breakthroughs in technical performance and engineering applicability in the field of borehole sealing quality detection. Compared with traditional methods, by using the continuous spatial sensing capability of distributed optical fiber and the high sensitivity characteristics of Q value to medium attenuation characteristics, the present application can quantitatively diagnose the coupling state of the wellbore sealing formation in the whole well section in real time and dynamically.

[0051] The core advantages are as follows: firstly, the technical performance is significantly improved, the Q value inversion is realized through the joint modeling of multiple physical fields such as amplitude attenuation, phase delay and envelope expansion, the anti-noise capability is stronger, and the inversion precision is higher; secondly, the distributed optical fiber has the characteristics of corrosion resistance, anti-electromagnetic interference and long service life, supports the whole life cycle monitoring of grouting, well cementing and wellbore service period, does not need to stop drilling or segmented operation, and can realize real-time early warning and accurate positioning of the sealing defect; thirdly, the detection cost is reduced by replacing multiple traditional logging with one laying. Compared with the prior art, the application breaks through the bottleneck of traditional methods in spatial continuity, detection accuracy and real-time performance, provides a high-precision, low-cost and full-cycle technical solution for sealing quality evaluation, and has a significant industry application prospect. BRIEF DESCRIPTION OF DRAWINGS

[0052] The application can be better understood by referring to the following description in conjunction with the accompanying drawings, which are incorporated in and form a part of the specification, and together with the detailed description, serve to explain the principles of the application. In the drawings:

[0053] Figure 1 A flow chart of a fiber-coupled medium quality factor inversion method for sealing quality detection is shown;

[0054] Figure 2 A structural block diagram of a fiber-coupled medium quality factor inversion device for sealing quality detection is shown. DETAILED DESCRIPTION

[0055] In the following, exemplary embodiments of the application will be described with reference to the accompanying drawings. In the specification, not all the features of the actual embodiments are described for the sake of clarity and conciseness. However, it should be understood that many embodiment-specific decisions can be made in the process of developing any such actual embodiment in order to achieve the specific goals of the developer, and these decisions can vary from embodiment to embodiment.

[0056] It should also be noted here that, in order to avoid obscuring the application due to unnecessary details, only the device structure closely related to the scheme according to the application is shown in the drawings, and other details not closely related to the application are omitted.

[0057] It should be understood that the application is not limited to the described embodiments only by the following description with reference to the drawings. In this text, the embodiments can be combined with each other, the features between different embodiments can be replaced or borrowed, and one or more features can be omitted in one embodiment, if possible.

[0058] In recent years, distributed fiber optic sensing technology has provided a new approach for monitoring grouting effects due to its advantages such as high precision, long distance, and resistance to electromagnetic interference. By pre-embedding the sensing fiber into the grouting hole, vibration signals can be acquired in real time during the grouting process and the solidification process of the slurry, and then the fiber-coupled medium quality factor can be inverted. The coupling medium quality factor reflects the wave velocity, density, porosity, and fluid saturation of the coupling medium. As the slurry solidifies over time, the density and wave velocity of the optical cable coupling medium increase, but the porosity and fluid saturation decrease significantly. Therefore, the sealing quality of the hole is evaluated by inverting the Q value. Traditional Q value estimation methods usually rely on a single physical field characteristic and are prone to large deviations under low signal-to-noise ratio or complex wave field interference.

[0059] This application achieves significant breakthroughs in both technical performance and engineering applicability in the field of wellbore sealing quality detection by integrating distributed fiber optic sensing technology with a multi-parameter Q-value inversion algorithm. Compared to traditional methods, this application leverages the continuous spatial sensing capabilities of distributed fiber optics and the high sensitivity of Q-values ​​to dielectric attenuation characteristics to enable full-well, dynamic, quantitative diagnosis of the wellbore sealing-stratum coupling state in real time.

[0060] The present invention provides a method for inverting the quality factor of a fiber-coupled medium for detecting the sealing quality. Figure 1 The flowchart of the fiber-coupled medium quality factor inversion method for sealing quality detection is shown in FIG. Figure 1 , the method mainly includes the following steps:

[0061] Step S1 , performing first arrival picking of a direct wave of a vibration signal on optical fiber monitoring data, and intercepting the optical fiber monitoring data based on the first arrival picking of the direct wave to obtain a vibration signal.

[0062] Here, the fiber optic monitoring data is first subjected to conventional amplitude normalization, and the calculation process adopts conventional seismic data processing algorithms.

[0063] Then, the Morlet wavelet is used as the wavelet basis function for wavelet transform. The maximum point of the wavelet coefficient modulus at the smallest scale is calculated. The maximum value is checked to see if it appears synchronously at adjacent scales. The first maximum point that satisfies multi-scale consistency is selected as the first arrival position of the direct wave. The first arrival of the direct wave of the vibration signal is picked up for each channel data to obtain the first arrival of the direct wave of each channel.

[0064] Find the average of the first arrival of the direct wave of all channels as the signal interception mark time;

[0065] The data of each channel is intercepted with the set time period before and after the signal interception mark moment (500 milliseconds before and after the interception mark moment) as the interception time window to obtain the vibration signal of each channel.

[0066] Step S2, amplitude compensation is performed on the vibration signal based on the spherical divergence compensation principle to obtain a compensated vibration signal.

[0067] Specifically, the vibration signal includes a plurality of channel vibration signals; an amplitude compensation coefficient of each channel is calculated:

[0068]

[0069] Wherein, G(i) is the amplitude compensation coefficient of the i-th channel, r i is the offset distance of the i-th channel;

[0070] The amplitude compensation is performed on the vibration signal of each channel using the corresponding amplitude compensation coefficient to obtain the compensated vibration signal of each channel.

[0071] Step S3, band-pass filtering is performed on the compensated vibration signal to obtain a filtered vibration signal.

[0072] Here, according to the main frequency of the vibration signal, a suitable effective frequency band is selected for band-pass filtering to obtain a vibration signal with high-frequency noise and low-frequency drift suppressed. It can be realized by using mature seismic data processing module.

[0073] Step S4, the filtered vibration signal is intercepted according to the calculation window length to obtain an intercepted signal; the Hilbert transform is performed on the intercepted signal to generate an analytic signal; and the instantaneous amplitude, instantaneous phase and spectrum are calculated according to the analytic signal.

[0074] Specifically, the analytic signal S a (t) is:

[0075] S a (t) = S win (t) + jS h (t)

[0076] Wherein, S win (t) is the intercepted signal, S h (t) is the Hilbert transform result, and t is the time within the calculation window.

[0077]

[0078] Wherein, τ is the integration time.

[0079] According to the analytic signal, the instantaneous amplitude A(t), the instantaneous phase and the spectrum spectrum are calculated, and the following formula is used:

[0080]

[0081] spectrum = |FFT(S win(t))

[0082] Step S5, generating an amplitude-based Q-value constructor according to the instantaneous amplitude and the amplitude attenuation model; generating a phase-based constructor according to the instantaneous phase; generating an envelope-based constructor according to the relative delay broadening effect of the spectrum and the envelope peak position; generating a joint loss function according to the amplitude-based Q-value constructor, the phase-based constructor and the envelope-based constructor; the Q-value is a fiber-coupled medium quality factor.

[0083] Step S6, using an L-BFGS-B optimization inversion algorithm to invert the joint loss function to obtain a Q-value jointly constrained by amplitude, phase and envelope.

[0084] Here, the Q-value jointly constrained by amplitude, phase and envelope includes a Q-value corresponding to each trace, and a Q-value variation curve with well depth can be obtained according to the relationship between the offset distance of each trace and the well depth, so that the sealing quality at different well depths can be evaluated.

[0085] In addition, based on the fiber monitoring data at different monitoring moments after grouting, the Q-value at different moments can be obtained according to steps S1-S6, and the sealing quality effect with the progress of the slurry solidification time can be evaluated.

[0086] The embodiment takes the grouting hole orifice knocking signal detected by the distributed optical fiber as the starting point, performs denoising, filtering and other pretreatments on the vibration signal, and inverses the coupling medium quality factor according to the amplitude, phase and envelope of the seismic wave on the basis of the pretreatment. The solidification state of the slurry at different depths in the hole before and after grouting is reflected by quantifying the energy absorption characteristics of the coupling medium, and the sealing quality is quantified. The traditional sealing quality evaluation method (such as acoustic logging and pressure testing) has the defects of low spatial resolution, high cost and inability to monitor in real time. Compared with the traditional technology, the application breaks through the bottleneck of the traditional method in resolution, real-time and quantification, and realizes high-precision, low-cost and full-life-cycle monitoring of the sealing quality.

[0087] In one embodiment, step S5, generating an amplitude-based Q-value constructor according to the instantaneous amplitude and the amplitude attenuation model, includes:

[0088] The instantaneous amplitude includes the instantaneous amplitudes of multiple traces, and for each trace, the instantaneous amplitudes at all moments in the time window are calculated to form the amplitude spectrum of the trace;

[0089] Any trace is regarded as an observation trace; according to the amplitude attenuation model, the theoretical reference amplitude spectrum after attenuation of the observation trace is determined, and the following formula is used:

[0090]

[0091] Wherein, A(f) is the reference amplitude spectrum of the theoretical attenuation of the observation channel, A0(f) is the amplitude spectrum of the reference channel, f is the frequency, Δt is the time difference of wave propagation from the reference channel to the observation channel, and Q is the quality factor of the optical fiber coupling medium; the reference channel is one of all channels, for example, it can be the first channel.

[0092] The mean square error between the amplitude spectrum of the observation channel and the reference amplitude spectrum A(f) of the theoretical attenuation of the observation channel is calculated as the amplitude-based Q value constructor L amp .

[0093] In one embodiment, step S5, the phase-based constructor is generated according to the instantaneous phase, including:

[0094] The instantaneous phase includes the instantaneous phase of each channel;

[0095] Any channel is taken as the observation channel; the mean square error between the instantaneous phase of the observation channel and the instantaneous phase of the reference channel is calculated as the phase-based constructor L phase .

[0096] In one embodiment, step S5, the envelope-based constructor is generated according to the relative delay broadening effect of the spectrum and the envelope peak position, including:

[0097] The spectrum includes the spectrum of each channel;

[0098] For each channel, the signal envelope is determined according to the spectrum, and the peak time of the signal envelope is determined; any channel is taken as the observation channel; the ratio of the peak time of the signal envelope of the observation channel to the peak time of the signal envelope of the reference channel is determined, which is denoted as the observation channel ratio;

[0099] When Q→∞, the theoretical ratio is 1, and Q is the quality factor of the optical fiber coupling medium;

[0100] The mean square error between the observation channel ratio and the theoretical ratio is calculated as the envelope-based constructor L env of the observation channel as a penalty term.

[0101] In one embodiment, step S5, the joint loss function is:

[0102] L Q =ω1L amp +ω2L phase +ω3L env

[0103] Wherein, L Q is the joint loss function, L amp is the amplitude-based Q value constructor, L phase is the phase-based constructor, and L envFor the envelope-based constructor, ω1+ω2+ω3=1, ω1 is the weight of L amp , ω2 is the weight of L phase , ω3 is the weight of L env , and the specific weight distribution is adjusted according to the data quality, with ω1 increased for high signal-to-noise ratio data and ω3 increased for low signal-to-noise ratio data.

[0104] In the above embodiment, the multi-physical field joint modeling of amplitude, phase, and envelope spread is stronger in anti-noise ability and higher in inversion accuracy.

[0105] Based on the same inventive concept as the optical fiber coupling medium quality factor inversion method for sealing quality detection, the embodiment also provides a corresponding optical fiber coupling medium quality factor inversion device for sealing quality detection, Figure 2 The structural block diagram of the optical fiber coupling medium quality factor inversion device for sealing quality detection is shown in Figure 2 , which comprises:

[0106] The first arrival picking and intercepting module picks up the direct wave first arrival of the vibration signal from the optical fiber monitoring data, and intercepts the optical fiber monitoring data based on the direct wave first arrival picking to obtain the vibration signal;

[0107] The amplitude compensation module compensates the amplitude of the vibration signal based on the spherical diffusion compensation principle to obtain the compensated vibration signal;

[0108] The band-pass filtering module is used for band-pass filtering the compensated vibration signal to obtain the filtered vibration signal;

[0109] The Hilbert transform module is used for intercepting the filtered vibration signal according to the calculation window length to obtain the intercepted signal; performing Hilbert transform on the intercepted signal to generate the analytic signal; and calculating the instantaneous amplitude, instantaneous phase, and spectrum according to the analytic signal;

[0110] The joint loss function generation module is used for generating a Q value construction function based on amplitude according to the instantaneous amplitude and amplitude decay model; generating a construction function based on phase according to the instantaneous phase; generating an envelope-based construction function according to the relative delay spread effect of the envelope peak position; and generating a joint loss function according to the Q value construction function based on amplitude, the construction function based on phase, and the envelope-based construction function; and the Q value is the quality factor of the optical fiber coupling medium;

[0111] The inversion module is used for performing inversion on the joint loss function by using the L-BFGS-B optimization inversion algorithm to obtain the Q value jointly constrained by amplitude, phase, and envelope.

[0112] The device for inversed calculation of the quality factor of the fiber-coupled medium for the sealing quality detection of the embodiment has the same inventive concept as the method for inversed calculation of the quality factor of the fiber-coupled medium for the sealing quality detection described above, and thus the specific implementation of the device can be seen from the embodiment part of the method for inversed calculation of the quality factor of the fiber-coupled medium for the sealing quality detection described above, and the technical effects thereof correspond to those of the method described above, which will not be repeated here.

[0113] The embodiment of the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the method for inversed calculation of the quality factor of the fiber-coupled medium for the sealing quality detection.

[0114] Compared with the prior art, the present application has the following technical effects:

[0115] 1. The distributed optical fiber has the characteristics of corrosion resistance, electromagnetic interference resistance and long service life, supports grouting, well cementing and whole life cycle monitoring of wellbore service period, uses the distributed optical fiber to evaluate the sealing quality, does not need to stop drilling or segmented operation, can realize real-time early warning and accurate positioning of sealing defects, and the engineering applicability is enhanced.

[0116] 2. The traditional Q value estimation method usually depends on a single physical field characteristic, and large deviation is easy to occur under low signal-to-noise ratio or complex wave field interference. By constructing a weighted loss function coupled with multiple physical fields, the amplitude attenuation model (revealing the high-frequency energy absorption law), the phase delay model (reflecting the dispersion effect) and the envelope broadening model (characterizing the wave shape time domain distortion) are optimized collaboratively. Compared with the traditional single field inversion method, the sensitive directions of the multiple physical field characteristics to the noise are orthogonal, the amplitude characteristics are significantly disturbed by random noise, while the phase delay and the envelope broadening are more sensitive to systematic errors, and by weight distribution, the noise interference in the DAS data can be effectively suppressed.

[0117] 3. Compared with the prior art, the present application breaks through the bottleneck of the traditional method in spatial continuity, detection accuracy and real-time performance, replaces multiple traditional logging by one-time laying, reduces the detection cost, provides a high-precision, low-cost and whole-cycle technical solution for sealing quality evaluation, and has a significant industry application prospect.

[0118] The above describes only various embodiments of the present application, but the protection scope of the present application is not limited thereto, and any person skilled in the art can easily think of changes or replacements within the technical scope disclosed in the present application, which should be covered in the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A fiber-coupled medium quality factor inversion method for sealing quality detection, characterized in that: include: Performing vibration signal direct wave first arrival picking on the optical fiber monitoring data, and intercepting the optical fiber monitoring data based on the direct wave first arrival picking to obtain a vibration signal; performing amplitude compensation on the vibration signal based on the spherical diffusion compensation principle to obtain a compensated vibration signal; performing bandpass filtering on the compensated vibration signal to obtain a filtered vibration signal; intercepting the filtered vibration signal according to the calculated time window length to obtain an intercepted signal; performing a Hilbert transform on the intercepted signal to generate an analytical signal; and calculating the instantaneous amplitude, instantaneous phase, and frequency spectrum based on the analytical signal; According to the instantaneous amplitude and amplitude attenuation model, an amplitude-based Q value constructor is generated; according to the instantaneous phase, a phase-based constructor is generated; according to the relative delay broadening effect of the spectrum and envelope peak position, an envelope-based constructor is generated; according to the amplitude-based Q value constructor, the phase-based constructor, and the envelope-based constructor, a joint loss function is generated; the Q value is the fiber coupling medium quality factor; The L-BFGS-B optimization inversion algorithm is used to invert the joint loss function to obtain the Q value of the amplitude, phase and envelope joint constraints.

2. The method according to claim 1, wherein in, Picking up the first arrival of a direct wave of a vibration signal on the optical fiber monitoring data, and intercepting the optical fiber monitoring data based on the first arrival of the direct wave to obtain a vibration signal, including: The optical fiber monitoring data includes data of multiple channels; Pick up the first arrival of the direct wave of the vibration signal for each channel data to obtain the first arrival of the direct wave of each channel; Find the average of the first arrival of the direct wave of all channels as the signal interception mark time; The time period before and after the signal interception mark is set as the interception time window, and the data of each channel is intercepted to obtain the vibration signal of each channel.

3. The method according to claim 1, wherein in, Performing amplitude compensation on the vibration signal based on the spherical diffusion compensation principle to obtain a compensated vibration signal includes: The vibration signal includes vibration signals of multiple channels; Calculate the amplitude compensation coefficient for each trace: Where G(i) is the amplitude compensation coefficient of the i-th track, r i is the offset of the ith track; The vibration signal of each track is amplitude compensated using the corresponding amplitude compensation coefficient to obtain the compensated vibration signal of each track.

4. The method according to claim 1, wherein in, According to the instantaneous amplitude and amplitude attenuation model, an amplitude-based Q value constructor is generated, including: The instantaneous amplitude includes the instantaneous amplitudes of multiple channels. For each channel, the instantaneous amplitudes at all moments within the calculation time window constitute the amplitude spectrum of the channel. Any channel is recorded as the observation channel; according to the amplitude attenuation model, the reference amplitude spectrum of the observation channel after theoretical attenuation is determined using the following formula: Where A(f) is the reference amplitude spectrum of the theoretically attenuated observation channel, A0(f) is the amplitude spectrum of the reference channel, f is the frequency, Δt is the time difference of the wave propagating from the reference channel to the observation channel, and Q is the quality factor of the fiber coupling medium. Calculate the mean square error between the amplitude spectrum of the observation channel and the reference amplitude spectrum A(f) after theoretical attenuation of the observation channel as the amplitude-based Q value constructor L of the observation channel amp .

5. The method according to claim 1, wherein in, According to the instantaneous phase, a phase-based constructor is generated; comprising: The instantaneous phase includes the instantaneous phase of each channel; Denote any channel as the observation channel; calculate the mean square error between the instantaneous phase of the observation channel and the instantaneous phase of the reference channel as the phase-based constructor L of the observation channel phase .

6. The method according to claim 1, wherein in, According to the relative delay broadening effect of the spectrum and envelope peak position, an envelope-based constructor is generated, including: The spectrum includes a spectrum of each channel; For each channel, determining a signal envelope according to the frequency spectrum, and determining a peak time of the signal envelope; Record any channel as an observation channel; determine the ratio of the peak time of the signal envelope of the observation channel to the peak time of the signal envelope of the reference channel, and record it as the observation channel ratio; When Q→∞, the theoretical ratio is 1, and Q is the quality factor of the fiber coupling medium; Calculate the mean square error between the observed trace ratio and the theoretical ratio as the envelope-based constructor L of the observed trace env .

7. The method according to claim 1, wherein The joint loss function is: L Q =ω1L amp +ω2L phase +ω3L env Among them, L Q is the joint loss function, L amp is the Q value constructor based on amplitude, L phase is a phase-based constructor, L env is the envelope-based constructor, ω1+ω2+ω3=1, ω1 is L amp The weight of L phase The weight of L env The weight of .

8. The method according to claim 1, wherein The method further comprises: Based on the optical fiber monitoring data at different monitoring times, the Q value of the amplitude, phase and envelope joint constraints at different monitoring times is obtained.

9. A fiber-coupled medium quality factor inversion device for sealing quality detection, characterized in that: include: A first arrival picking and intercepting module is used to pick up the first arrival of the direct wave of the vibration signal from the optical fiber monitoring data, and intercept the optical fiber monitoring data based on the first arrival picking of the direct wave to obtain the vibration signal; an amplitude compensation module, performing amplitude compensation on the vibration signal based on the spherical diffusion compensation principle to obtain a compensated vibration signal; a bandpass filtering module, configured to perform bandpass filtering on the compensated vibration signal to obtain a filtered vibration signal; A Hilbert transform module is configured to intercept the filtered vibration signal according to a calculated time window length to obtain an intercepted signal; perform a Hilbert transform on the intercepted signal to generate an analytical signal; and calculate the instantaneous amplitude, instantaneous phase, and frequency spectrum based on the analytical signal; A joint loss function generation module is configured to generate an amplitude-based Q value constructor according to the instantaneous amplitude and amplitude attenuation model; generate a phase-based constructor according to the instantaneous phase; generate an envelope-based constructor according to the relative delay broadening effect of the spectrum and envelope peak position; and generate a joint loss function according to the amplitude-based Q value constructor, the phase-based constructor, and the envelope-based constructor; the Q value is the fiber coupling medium quality factor; The inversion module is used to invert the joint loss function using the L-BFGS-B optimization inversion algorithm to obtain the Q value of the amplitude, phase and envelope joint constraints.

10. A computer-readable storage medium, characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method for inverting the quality factor of a fiber-coupled medium for sealing quality detection according to any one of claims 1 to 8 is implemented.

Images