Method and System for Correcting Imbalance in In-phase/Quadrature Demodulation of Optical fiber DAS Data

20240337525 ยท 2024-10-10

Assignee

Inventors

Cpc classification

International classification

Abstract

A method for correcting imbalance in IQ demodulation of optical fiber DAS data is disclosed. The method includes obtaining I/Q signals I0 and Q0 of each sampling on an optical fiber from a DAS acquisition instrument; performing DC bias correction on the I/Q signals I0 and Q0 to obtain I/Q signals I1 and Q1 after DC bias correction; performing amplitude imbalance correction on the I/Q signals I1 and Q1 after DC bias correction by using Hilbert transform to obtain I/Q signals I2 and Q2 after amplitude imbalance correction; and performing phase imbalance correction on the I/Q signals I2 and Q2 after amplitude imbalance correction by using Hilbert transform to obtain I/Q signals I3 and Q3 after phase imbalance correction. It is possible to perform accurate correction on I/Q signals, accurately and efficiently suppress optical fiber demodulation noise, and improve optical fiber DAS data acquisition quality.

Claims

1. A method for correcting imbalance in IQ demodulation of optical fiber DAS data, comprising: obtaining I/Q signals I.sub.0 and Q.sub.0 of each sampling; performing DC bias correction on the I/Q signals Io and Q.sub.0 to obtain I/Q signals I.sub.1 and Q.sub.1 after DC bias correction; performing amplitude imbalance correction on the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction by using Hilbert transform to obtain I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction; and performing phase imbalance correction on the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction by using Hilbert transform to obtain I/Q signals I.sub.3 and Q.sub.3 after phase imbalance correction.

2. The method according to claim 1, wherein the step of performing DC bias correction on the I/Q signals I.sub.0 and Q.sub.0 to obtain I/Q signals I.sub.1 and Q.sub.1 after DC bias correction comprises: calculating corresponding DC bias factors D.sub.i and D.sub.q according to the I/Q signals I.sub.0 and Q.sub.0; and performing DC bias correction on the I/Q signals I.sub.0 and Q.sub.0 by using the DC bias factors D.sub.i and D.sub.q to obtain the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction.

3. The method according to claim 2, wherein the step of calculating corresponding DC bias factors D.sub.i and D.sub.q according to the I/Q signals I.sub.0 and Q.sub.0 comprises: calculating the corresponding DC bias factors D.sub.i and D.sub.q according to a mean value of the I/Q signals I.sub.0 and Q.sub.0 at different time.

4. The method according to claim 1, wherein the step of performing amplitude imbalance correction on the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction by using Hilbert transform to obtain I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction comprises: calculating corresponding amplitude imbalance factors A.sub.i and A.sub.q by using Hilbert transform and the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction; and performing amplitude imbalance correction on the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction according to the amplitude imbalance factors A.sub.i and A.sub.q to obtain the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction.

5. The method according to claim 4, wherein the step of calculating corresponding amplitude imbalance factors A.sub.i and A.sub.q by using Hilbert transform and the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction comprises: calculating corresponding instantaneous amplitudes of the I/Q signals I.sub.1 and Q.sub.1 by Hilbert transform; and calculating the corresponding amplitude imbalance factors A.sub.i and A.sub.q according to a mean value of corresponding instantaneous amplitudes at different time.

6. The method according to claim 1, wherein the step of performing phase imbalance correction on the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction by using Hilbert transform to obtain I/Q signals I.sub.3 and Q.sub.3 after imbalance phase correction comprises: calculating phases ?.sub.i and ?.sub.q corresponding to the I/Q signals I.sub.2 and Q.sub.2 by using Hilbert transform and the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction; calculating a phase imbalance factor ? according to the phases ?.sub.i and ?.sub.q corresponding to the I/Q signals I.sub.2 and Q.sub.2; and performing phase imbalance correction on the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction according to the phase imbalance factor ? to obtain the I/Q signals I.sub.3 and Q.sub.3 after phase imbalance correction.

7. The method according to claim 6, wherein the step of calculating phases ?.sub.i and ?.sub.q corresponding to the I/Q signals I.sub.2 and Q.sub.2 by using Hilbert transform and the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction comprises: obtaining the corresponding phases ?.sub.i and ?.sub.q by an arctangent function of a ratio of the I/Q signals I.sub.2 and Q.sub.2 to the corresponding I/Q signals after Hilbert transform.

8. The method according to claim 6, wherein the step of calculating a phase imbalance factor ? according to the phases ?.sub.i and ?.sub.q corresponding to the I/Q signals comprises: determining the phase imbalance factor ? by using a mean value of differences of the phases ?.sub.i and ?.sub.q at different time.

9. The method according to claim 6, wherein performing phase imbalance correction on the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction according to the phase imbalance factor ? by the following manner to obtain the I/Q signals I.sub.3 and Q.sub.3 after phase imbalance correction: I 3 = I 2 ; Q 3 = Q 2 - I 2 sin ? cos ? .

10. A system for correcting imbalance in IQ demodulation of optical fiber DAS data, comprising: an obtaining module, configured to obtain I/Q signals I.sub.0 and Q.sub.0 of each sampling on an optical fiber; a DC bias correction module, configured to perform DC bias correction on the I/Q signals I.sub.0 and Q.sub.0 to obtain I/Q signals I.sub.1 and Q.sub.1 after DC bias correction; an amplitude imbalance correction module, configured to perform amplitude imbalance correction on the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction by using Hilbert transform to obtain I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction; and a phase imbalance correction module, configured to perform phase imbalance correction on the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction by using Hilbert transform to obtain I/Q signals I.sub.3 and Q.sub.3 after phase imbalance correction.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0018] The drawings are used to provide a further understanding of embodiments of the present invention and constitute a part of this specification, and to explain the embodiments of the present invention together with specific implementation modes, but do not limit the embodiments of the present invention. In the drawings:

[0019] FIG. 1 is a flowchart illustrating a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0020] FIG. 2 is a schematic diagram illustrating I/Q signals I.sub.0(t.sub.i) and Q.sub.0(t.sub.i) at sample points when there is no external acoustic signal in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0021] FIG. 3 is a schematic diagram illustrating I/Q signals I.sub.1(t.sub.i) and Q.sub.1(t.sub.i) after DC bias correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0022] FIG. 4 is a schematic diagram illustrating I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) after amplitude imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0023] FIG. 5 is a schematic diagram illustrating a phase curve and a phase difference curve of I/Q signals obtained by Hilbert transform in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0024] FIG. 6 is a schematic diagram illustrating I/Q signals I.sub.3(t.sub.i) and Q.sub.3(t.sub.i) after imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0025] FIG. 7 is a zoomed view illustrating I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) after amplitude imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0026] FIG. 8 is a partially enlarged view illustrating I/Q signals I.sub.3(t.sub.i) and Q.sub.3(t.sub.i) after phase imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0027] FIG. 9 is a schematic diagram illustrating a demodulation phase curve and a demodulation phase difference curve before and after imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0028] FIG. 10 is a schematic diagram illustrating amplitude curves before and after imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0029] FIG. 11 is a schematic diagram illustrating I/Q cross plots before and after imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0030] FIG. 12 is a schematic diagram illustrating IQ demodulation background noise before and after imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention;

[0031] FIG. 13 is a schematic diagram illustrating IQ demodulation seismic wavefields before and after imbalance correction in a method for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention; and

[0032] FIG. 14 is a schematic diagram illustrating a system for correcting imbalance in IQ demodulation of optical fiber DAS data according to an embodiment of the present invention.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0033] The specific implementation modes of the embodiments of the present invention are described in detail below in combination with the drawings. It should be understood that the specific implementation modes described herein are only used for describing and explaining the embodiments of the present invention and are not used for limiting the embodiments of the present invention.

[0034] The obtainment, storage, use, processing, or the like of data in the technical solutions of the present application comply with the relevant provisions of national laws and regulations.

[0035] Referring to FIG. 1, the present invention, in a first aspect, provides a method for correcting imbalance in IQ demodulation of optical fiber DAS data, including the following steps.

[0036] S101: I/Q signals I.sub.0 and Q.sub.0 of each sampling on an optical fiber are obtained.

[0037] Specifically, in an embodiment of the present invention, an optical fiber DAS acquisition instrument is connected to an optical fiber coupled to a to-be-measured object, and I/Q signals I.sub.0 and Q.sub.0 are output at each sampling on the optical fiber in the absence of external acoustic signals around the optical fiber in the case of emitting optical pulses at a certain frequency. The optical fiber DAS acquisition instrument is an instrument with a phase-sensitive optical time domain reflectometer (?-OTDR) as a core. The I/Q signals are in-phase and quadrature signals, and it is assumed that there is no external acoustic signals near each sampling on the optical fiber.

[0038] S102: DC bias correction is performed on the I/Q signals I.sub.0 and Q.sub.0 to obtain I/Q signals I.sub.1 and Q.sub.1 after DC bias correction.

[0039] Further, the step that DC bias correction is performed on the I/Q signals I.sub.0 and Q.sub.0 to obtain I/Q signals I.sub.1 and Q.sub.1 after DC bias correction includes the following step: Corresponding DC bias factors D.sub.i and D.sub.q are calculated according to a mean value of the I/Q signals I.sub.0 and Q.sub.0 at different time:

[00002] D i = 1 n .Math. i = 1 n I 0 ( t i ) ; D q = 1 n .Math. i = 1 n Q 0 ( t i ) ;

wherein I.sub.0(t.sub.i) and Q.sub.0(t.sub.i) are I/Q signals I.sub.0 and Q.sub.0 at time t.sub.i, n is the number of samples for calculation.

[0040] DC bias correction is performed on the I/Q signals I.sub.0 and Q.sub.0 by using the DC bias factors D.sub.i and D.sub.q to obtain the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction:

[00003] I 1 ( t i ) = I 0 ( t i ) - D i ; Q 1 ( t i ) = Q 0 ( t i ) - D q ;

wherein I.sub.1 (t.sub.i) and Q.sub.1 (t.sub.i) are I/Q signals I.sub.1 and Q.sub.1 at time t.sub.i.

[0041] S103: Amplitude imbalance correction is performed on the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction by using Hilbert transform to obtain I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction.

[0042] Further, the step that amplitude imbalance correction is performed on the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction by using Hilbert transform to obtain I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction includes the following steps: Corresponding instantaneous amplitudes of the I/Q signals I.sub.1 and Q.sub.1 are calculated by Hilbert transform; and Corresponding amplitude imbalance factors A.sub.i and A.sub.q are calculated according to a mean value of corresponding instantaneous amplitudes at different time:

[00004] A i = 1 n .Math. i = 1 n I 1 2 ( t i ) + I ? 1 2 ( t i ) ; A q = 1 n .Math. i = 1 n Q 1 2 ( t i ) + Q ~ 1 2 ( t i ) ;

wherein custom-character(t.sub.i) and custom-character(t.sub.i) are I/Q signals obtained by Hilbert transform on the I/Q signals I.sub.1(t.sub.i) and Q.sub.1(t.sub.i) after DC bias correction, n is the number of samples for calculation.

[0043] Amplitude imbalance correction is performed on the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction according to the amplitude imbalance factors A.sub.i and A.sub.q to obtain the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction:

[00005] I 2 { t i ) = I 1 ( t i ) / A i ; Q 2 { t i ) = Q 1 ( t i ) / A q ;

wherein I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) are I/Q signals I.sub.2 and Q.sub.2 at time t.sub.i.

[0044] S104: Phase imbalance correction is performed on the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction by using Hilbert transform to obtain I/Q signals I.sub.3 and Q.sub.3 after phase imbalance correction.

[0045] Further, the step that phase imbalance correction is performed on the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction by using Hilbert transform to obtain I/Q signals I.sub.3 and Q.sub.3 after imbalance correction includes the following steps: Corresponding phases ?.sub.i and ?.sub.q are obtained by an arctangent function of a ratio of the I/Q signals I.sub.2 and Q.sub.2 to the corresponding I/Q signals after Hilbert transform:

[00006] ? i ( t i ) = arctan ( I 2 ~ ( t i ) I 2 ( t i ) ) ; ? q ( t i ) = arctan ( Q 2 ~ ( t i ) Q 2 ( t i ) ) ;

wherein custom-character(t.sub.i) and custom-character(t.sub.i) are I/Q signals obtained by Hilbert transform on the I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) after amplitude imbalance correction.

[0046] A phase imbalance factor ? is determined by using a mean value of differences of the phases ?.sub.i and ?.sub.q at different time:

[00007] ? = 1 n .Math. i = 1 n ( ? q ( t i ) - ? i ( t i ) - ? 2 ) .

[0047] Phase imbalance correction is performed on the I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) after amplitude imbalance correction according to the phase imbalance factor ? to obtain I/Q signals I.sub.3(t.sub.i) and Q.sub.3(t.sub.i) after imbalance correction:

[00008] I 3 ( t i ) = I 2 ( t i ) ; Q 3 ( t i ) = Q 2 ( t i ) - I 2 ( t i ) sin ? cos ? .

[0048] Then, subsequent phase demodulation and phase unwrapping may be performed by using the above data to obtain a demodulation phase ? and an amplitude Amp, thereby obtaining acoustic wavefield data. The demodulation phase ? and the amplitude Amp are calculated as follows:

[00009] ? ( t i ) = arctan ( Q 3 ( t i ) I 3 ( t i ) ) ; Amp ( t i ) = I 3 2 ( t i ) + Q 3 2 ( t i ) .

[0049] In an embodiment of the present invention, the DC bias factors D.sub.i and D.sub.q, the amplitude imbalance factors A.sub.i and A.sub.q, and the phase imbalance factor ? may be obtained by selecting positions without external acoustic signals uniformly after the I/Q signals I.sub.0 and Q.sub.0 at each sampling on the optical fiber are obtained, and then are applied to all acquired data at all samplings.

Embodiment

[0050] A optical fiber DAS acquisition instrument is connected to an optical fiber hung in a cased well, in the case of emitting optical pulses at a repetition rate of 10 K, I/Q signals I.sub.0(t.sub.i) and Q.sub.0(t.sub.i) output at time t.sub.i(i=1, 2, 3, . . . ) at each sampling on the optical fiber without external acoustic signals around the optical fiber are acquired when a vibroseis is not excited, as shown in FIG. 2.

[0051] Corresponding DC bias factors D.sub.i and D.sub.q are calculated according to a mean value of the I/Q signals I.sub.0(t.sub.i) and Q.sub.0(t.sub.i) at multiple points:

[00010] D i = 1 n .Math. i = 1 n I 0 ( t i ) ; D q = 1 n .Math. i = 1 n Q 0 ( t i ) ;

[0052] DC bias correction is performed on the I/Q signals I.sub.0(t.sub.i) and Q.sub.0(t.sub.i) by using the DC bias factors D.sub.i and D.sub.q to obtain I/Q signals I.sub.1(t.sub.i) and Q.sub.1(t.sub.i) after DC bias correction, as shown in FIG. 3:

[00011] I 1 ( t i ) = I 0 ( t i ) - D i ; Q 1 ( t i ) = Q 0 ( t i ) - D q .

[0053] Then, corresponding amplitude imbalance factors A.sub.i and A.sub.q are calculated by using Hilbert transform and the I/Q signals I.sub.1(t) and Q.sub.1(t) after DC bias correction:

[00012] A i = 1 n .Math. i = 1 n I 1 2 ( t i ) + I ? 1 2 ( t i ) ; A q = 1 n .Math. i = 1 n Q 1 2 ( t i ) + Q ~ 1 2 ( t i ) ;

wherein custom-character(t.sub.i) and custom-character(t.sub.i) are I/Q signals obtained by Hilbert transform on the I/Q signals I.sub.1(t.sub.i) and Q.sub.1(t.sub.i) after DC bias correction.

[0054] Amplitude imbalance correction is performed on the I/Q signals I.sub.1(t.sub.i) and Q.sub.1(t.sub.i) after DC bias correction according to the amplitude imbalance factors A.sub.i and A.sub.q to obtain I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) after amplitude imbalance correction, as shown in FIG. 4:

[00013] I 2 ( t i ) = I 1 ( t i ) / A i ; Q 2 ( t i ) = Q 1 ( t i ) / A q .

[0055] Then, phases ?.sub.i(t.sub.i) and ?.sub.q(t.sub.i) corresponding to the the I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) are calculated by using Hilbert transform and the I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) after amplitude imbalance correction, as shown in FIG. 5:

[00014] ? i ( t i ) = arctan ( I 2 ~ ( t i ) I 2 ( t i ) ) ; ? q ( t i ) = arctan ( Q 2 ~ ( t i ) Q 2 ( t i ) ) ;

wherein, custom-character(t.sub.i) and custom-character(t.sub.i) are I/Q signals obtained by Hilbert transform on I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) after amplitude imbalance correction.

[0056] Referring to FIG. 5, a phase imbalance factor is calculated according to the phases ?.sub.i(t.sub.i) and ?.sub.q(t.sub.i) corresponding to the I/Q signals, and the I/Q signals themselves are quadrature signals, and have a phase difference

[00015] ? 2 ,

so that the phase imbalance factor (phase difference) ? is:

[00016] ? = 1 n .Math. i = 1 n ( ? q ( t i ) - ? i ( t i ) - ? 2 ) .

[0057] In this embodiment, the mean value of the phase difference is ?14 degrees.

[0058] Phase imbalance correction is performed on the I/Q signals I.sub.2(t.sub.i) and Q.sub.2(t.sub.i) after amplitude imbalance correction according to the phase imbalance factor to obtain I/Q signals I.sub.3(t.sub.i) and Q.sub.3(t.sub.i) after imbalance correction, as shown in FIGS. 6-8, wherein FIG. 7 and FIG. 8 are zoomed views of the I/Q signals before and after phase imbalance correction:

[00017] I 3 ( t i ) = I 2 ( t i ) ; Q 3 ( t i ) = Q 2 ( t i ) - I 2 ( t i ) sin ? cos ? .

[0059] Subsequent phase demodulation and phase unwrapping may be performed by using the above data after the I/Q signals imbalance correction are obtained so as to obtain a demodulation phase ? and an amplitude Amp, thereby obtaining acoustic wavefield data. The demodulation phase ? and the amplitude Amp are calculated as follows:

[00018] ? ( t i ) = arctan ( Q 3 ( t i ) I 3 ( t i ) ) ; Amp ( t i ) = I 3 2 ( t i ) + Q 3 2 ( t i ) .

[0060] Referring to FIG. 9-FIG. 12, a demodulation phase curve and a demodulation phase difference curve before and after phase imbalance correction are as shown in FIG. 9. Amplitude curves before and after phase imbalance correction are as shown in FIG. 10. I/Q crossplots before and after imbalance correction are as shown in FIG. 11, and it can be seen that the crossplot is elliptical before correction and is circular after correction. IQ demodulation background noise before and after imbalance correction is as shown in FIG. 12, and the noise level is reduced by 60% or above by the method of the present invention. IQ demodulation seismic wavefields before and after imbalance correction are as shown in FIG. 13, and a signal-to-noise ratio is improved by 3-4 times by the method of the present invention.

[0061] Referring to FIG. 14, the present invention, in a second aspect, provides a system for correcting imbalance in IQ demodulation of optical fiber DAS data. The system includes: an obtaining module, configured to obtain I/Q signals I.sub.0 and Q.sub.0 of each sampling on an optical fiber; a DC bias correction module, configured to perform DC bias correction on the I/Q signals I.sub.0 and Q.sub.0 to obtain I/Q signals I.sub.1 and Q.sub.1 after DC bias correction; an amplitude imbalance correction module, configured to perform amplitude imbalance correction on the I/Q signals I.sub.1 and Q.sub.1 after DC bias correction by using Hilbert transform to obtain I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction; and a phase imbalance correction module, configured to perform phase imbalance correction on the I/Q signals I.sub.2 and Q.sub.2 after amplitude imbalance correction by using Hilbert transform to obtain I/Q signals I.sub.3 and Q.sub.3 after phase imbalance correction.

[0062] The preferred embodiments of the present invention are described above in detail in combination with the drawings, but the present invention is not limited to the specific details in the above embodiments. Various simple variations can be made to the technical solutions of the present invention in the scope of the technical concept of the present invention, and these simple variations all fall within the protection scope of the present invention.

[0063] In addition, it should be noted that the specific technical features described in the above specific embodiments can be combined in any appropriate mode under the non-contradictory condition, and various possible combination modes will not be described separately in order to avoid unnecessary repetition.

[0064] Furthermore, various different embodiments of the present invention can also be combined optionally, and the combinations also should be regarded as the content disclosed by the present invention as long as the combinations do not violate the concept of the present invention.