REDUCED COMPLEXITY POLARIZATION COMBINING METHOD FOR COHERENT DAS

Abstract

Aspects of the present disclosure are directed to improved systems, methods, and structures providing coherent detection of DAS. In sharp contrast to the prior art, systems, methods, and structures according to aspects of the present disclosure advantageously reduce the beating diversity terms such that required memory and bandwidth are reduced over the art.

Claims

1. A reduced complexity polarization combining method for a coherent distributed acoustic sensing (DAS) system including a length of optical sensing fiber, an optical interrogator that generates optical pulses, introduces them into the optical fiber and receives Rayleigh reflected signals from the fiber, and a coherent receiver configured to extract information from the Rayleigh reflected signals, the method comprising: operating the DAS system to obtain a set of in-phase and quadrature time-domain data; x.sub.i(t), x.sub.q(t), y.sub.i(t), y.sub.q(t); digitizing the set of in-phase and quadrature time-domain data to generate a set of in-phase and quadrature data corresponding to the in-phase and quadrature time-domain data, x.sub.i(t), x.sub.q(t), y.sub.i(t), y.sub.q(t); optionally pre-processing the digitized set of in-phase and quadrature data; to obtain a pre-processed set of in-phase and quadrature data x.sub.i′(t), x.sub.q′(t), y.sub.i′(t), y.sub.q′(t); x-y-combining the polarizations of optionally pre-processed digitized set of in-phase and quadrature data to produce an x-y-combined output S(n); performing a differential beating on the xy-combined output to generate a single polarization diversity (n); and post-processing the generated single polarization diversity to generate an output indicative of acoustic impulse(s) interacting with the optical sensing fiber.

2. The method of claim 1 further comprising: the x-y-combining includes a two-step rotation having an x-y alignment step followed by a phase continuity step.

3. The method of claim 2 wherein the x-y alignment step rotates a polarization of a lower averaged power data to a higher averaged power data.

4. The method of claim 3 wherein the phase continuity step rotates an updated angle such that phase continuity is maintained.

5. The method of claim 4 wherein an averaged phase difference P.sub.diff(n) used in the rotations is the normalized to P.sub.diff(n)/|P.sub.diff(n)|, where P.sub.diff(n) is the phase difference at data n.

Description

BRIEF DESCRIPTION OF THE DRAWING

[0011] A more complete understanding of the present disclosure may be realized by reference to the accompanying drawing in which:

[0012] FIG. 1(A) is a schematic diagram of an illustrative prior art DAS arrangement;

[0013] FIG. 1(B) is a schematic diagram of an illustrative coherent DAS arrangement according to aspects of the present disclosure;

[0014] FIG. 2(A) is a schematic diagram of an illustrative prior art coherent DAS operational/processing arrangement;

[0015] FIG. 2(B) is a schematic diagram of an illustrative coherent DAS operational/processing arrangement according to aspects of the present disclosure;

[0016] FIG. 3 is a flow diagram outlining operation of an illustrative coherent DAS according to aspects of the present disclosure; and

[0017] FIG. 4 is a flow diagram describing an illustrative method according to aspects of the present disclosure;

[0018] The illustrative embodiments are described more fully by the Figures and detailed description. Embodiments according to this disclosure may, however, be embodied in various forms and are not limited to specific or illustrative embodiments described in the drawing and detailed description.

DESCRIPTION

[0019] The following merely illustrates the principles of the disclosure. It will thus be appreciated that those skilled in the art will be able to devise various arrangements which, although not explicitly described or shown herein, embody the principles of the disclosure and are included within its spirit and scope.

[0020] Furthermore, all examples and conditional language recited herein are intended to be only for pedagogical purposes to aid the reader in understanding the principles of the disclosure and the concepts contributed by the inventor(s) to furthering the art and are to be construed as being without limitation to such specifically recited examples and conditions.

[0021] Moreover, all statements herein reciting principles, aspects, and embodiments of the disclosure, as well as specific examples thereof, are intended to encompass both structural and functional equivalents thereof. Additionally, it is intended that such equivalents include both currently known equivalents as well as equivalents developed in the future, i.e., any elements developed that perform the same function, regardless of structure.

[0022] Thus, for example, it will be appreciated by those skilled in the art that any block diagrams herein represent conceptual views of illustrative circuitry embodying the principles of the disclosure.

[0023] Unless otherwise explicitly specified herein, the FIGs comprising the drawing are not drawn to scale.

[0024] By way of some additional background—we again note that in recent years, distributed fiber optic sensing (DFOS) systems including distributed vibration sensing (DVS) and distributed acoustic sensing (DAS) have found widespread acceptance in numerous applications including—but not limited to—infrastructure monitoring, intrusion detection, and earthquake detection. For DAS and DVS, backward Rayleigh scattering effects are used to detect changes in the fiber strain, while the fiber itself acts as the transmission medium for conveying the optical sensing signal back to an interrogator for subsequent analysis.

[0025] FIG. 1(A) shows a simplified schematic diagram of a prior art DFOS/DVS/DAS system employing an interrogator/coherent receiver/detection/analysis system. Operationally, such a system generates/applies an optical Tx signal into an optical sensing fiber which results in reflected/scattered optical signals being directed back to the receiver/analysis system which receives/detects/analyzes the reflected/scattered and subsequently received signals. The signals are analyzed, and an output is generated which is indicative of the environmental conditions encountered along the length of the fiber.

[0026] FIG. 1(B) is a schematic block diagram of an illustrative coherent DAS system according to aspects of the present disclosure. Since a coherent DAS system uses Rayleigh scattering to detect acoustic signal along a sensing fiber, it generally employs an interrogator that includes a transmitter (Tx) that generates an optical signal (pulse or code) periodically. The optical signal is directed into a distributed sensing fiber. Each location along the fiber reflects a small portion of the optical signal back to the interrogator. The backscattered signal is processed in the Rx section to recover the acoustic signal or detect the vibration.

[0027] As previously noted, coherent DAS uses differential beating for every two selected locations along an optical fiber to detect fiber stress at location(s) in between the two selected locations. Coherent optical detection has X and Y polarization diversities, which changes randomly due to fiber movement or other factors. For this reason, the beating may use X-X, X-Y, Y X, and Y-Y to fully utilize all the power, which results in 4 polarization diversities ζ.sub.xx, ζ.sub.xy, ζ.sub.yx, and ζ.sub.yy. Subsequent processing is required to combine the 4 diversity terms into a single term.

[0028] DAS received signal samples are received in sequence of location-by-location within each frame, while the polarization diversity combining process requires a frame-by-frame processing for each location. The sequence conversion requires large amount of memory and bandwidth. Doubling the diversity terms from beating process further doubles the memory and bandwidth needed.

[0029] Systems, methods and structures according to aspects of the present disclosure generally operate within or in conjunction with the receiver, and advantageously reduces the memory and bandwidth required by reducing beating diversity terms.

[0030] According to aspects of the present disclosure, X and Y polarizations are merged before beating, since polarization switching is a slow process as compared to location sampling rate (i.e., DAS pulse or frame repetition rate). Operationally, the two polarizations are first aligned to the same direction before merging, by rotating one of the polarizations (X or Y) to the other (Y or X), then rotated to maintain phase continuity.

[0031] The two polarizations first align to the one having higher averaged power (say pol-P). The X-Y combined signal is then passed to the beating module for differential beating, followed by phase extraction or other additional processing.

[0032] Advantageously, systems, methods, and structures according to aspects of the present disclosure combines the two polarizations into one output before beating, such that there is only a single input to a beating module and only one output from beating. This overall inventive operation advantageously reduces the processing complexity and memory size.

[0033] According to additional aspects of the present disclosure, a first rotation is performed and used to align the two polarizations; uses a second rotation to maintain phase continuity when the higher power polarization changes from one to the other. This second rotation guarantees that the phase difference between two samples at the beating output completely reflects the signal plus noise, not other effect added by the processing.

[0034] For this discussion, we use x for temporal average (or low-pass filtered value), x* for complex conjugate. As previously noted, A coherent DAS system uses Rayleigh scattering to detect acoustic signal along the connected fiber. It includes a transmitter (Tx) section that periodically generates an optical pulse signal or code signal periodically. The signal is directed into a distributed sensing fiber. Each location along the fiber reflects/scatters a small portion of the optical signal back to the interrogator. The backscattered signal is processed in a Rx section to recover acoustic signals or detect vibrations. Our focus of this disclosure is the receiver section, in particular the combination of signals exhibiting different polarization diversities into a single signal.

[0035] FIGS. 2(A) and 2(B) illustrate coherent DAS receiver processing according to the prior art (FIG. 2(A) and the present disclosure (FIG. 2(B)). As shown in the figure(s), an A/D is an analog-to-digital converter/process that samples a continuous time domain signal (analog) and converts same to a discrete digitized signal (digital). Next, the resulting signals are pre-processed which may include frequency conversion, filtering, down-sampling, and/or decoding. Differential beating is performed on the pre-processed signals and generates an inner product of two signals which are representative of represent the phase change. Polarization combining merges multiple polarization diversity terms into a single signal. The resulting single signal is post processed which may include phase unwrapping from complex signal(s), signal filtering, and subsequent vibration detection provided as output.

[0036] As shown in FIG. 2(B), DAS Rx processing according to aspects of the present disclosure advantageously performs X-Y combining prior to the differential beating processing. Importantly—rather than processing four polarization diversity terms ζ.sub.xx, ζ.sub.xy, ζ.sub.yx, and ζ.sub.yy following differential beating output—systems, methods and structures according to aspects of the present disclosure performs polarization combining before differential beating such that it (differential beating) receives a single input S(n) and outputs only a a single ζ(n) output.

[0037] In coherent DAS, a signal R.sub.s at time n has an associated location z, which is the relative time to the start of the interrogating pulse or code. Polarization combining such as that illustrated in FIG. 2(A) treats each location independently, which effectively processes R.sub.s(z,m) for each location z over time m.

[0038] To simplify the discussion, the following sections focus on a single location z, unless otherwise specified.

[0039] The procedure to achieve this is performed through a two steps' rotation. The first step is an X-Y alignment, that rotates the polarization of the one exhibiting a lower averaged power to the polarization of the higher one. The second step performed involves rotating by a dynamically updated angle to maintain phase continuity. The detailed flow chart of the procedure is given in FIG. 3, with simplified Matlab code shown in the listing that follows. Note that the averaged phase difference P.sub.diff(n) used in the rotations is the normalized to P.sub.diff(n)/|P.sub.diff(n)|, as described in the Matlab pseudocode listing.

[0040] A simplified MATLAB code listing showing this first rotation and second rotation according to aspects of the present disclosure is as follows:

TABLE-US-00001 MATLAB Program Listing      % calculate x-y phase difference and averaged difference diff = x.* conj(y); % phase difference Pdiff(1) = a*diff(1); % averaged phase difference. ‘a’ is low pass filter coefficient, using % y(n)=(1−a)*y(n−1)+a*x(n) for i=2:siglen, % siglen: length of signal x and y  Pdiff(i) = (1−a)*Pdiff(i−1) + a*diff(i); end      % first rotation, to align x and y xr = zeros(1, siglen); % aligned and combined x-y signals for i=2:siglen,  if xa(i)>ya(i), % xa / ya: X / Y polarization average power % rotate y to x   xr(i) = x(i) + y(i)*Pdiff(i)/abs(Pdiff(i));  else % rotate x to y   xr(i) = y(i) + x(i)*conj(Pdiff(i))/abs(Pdiff(i));  end end      % second rotation, to keep phase continuity R = zeros(1, siglen); % complex signal for 2.sup.nd rotation R(1) = 1+0*j; for i=2:siglen, % angle to rotate  if xa(i)>ya(i) && xa(i−1)<=ya(i−1), % higher power polarization changes from y to x   R(i) = R(i−1) * conj(Pdiff(i))/abs(Pdiff(i));  elseif xa(i)<=ya(i) && xa(i−1)>ya(i−1), % higher power polarization changes from x to y   R(i) = R(i−1) * Pdiff(i)/abs(Pdiff(i));  else   R(i) = R(i−1);  end      % rotate  xr(i) = xr(i)*R(i); end

[0041] The polarization combining illustratively shown in the flow chart of FIG. 3 uses the original amplitude of x and y rotated by normalized P.sub.diff. In one embodiment, the combining may apply a weight, using the averaged amplitude as:

xr(i)=x(i)*xa(i)/(xa(i)+ya(i))+y(i)*Pdiff(i)/abs(Pdiff(i))*ya(i)/(xa(i)+ya(i));

or:

xr(i)=y(i)*ya(i)/(xa(i)+ya(i))+x(i)*conj(Pdiff(i))/abs(Doff(i))*xa(i)/(xa(i)+ya(i)).

[0042] To simplify the implementation, the first rotation and combining can be reduced by using the higher power signal only, namely:

[00001]$xr\left(i\right)=\{\begin{array}{c}x\left(i\right),\mathrm{if}\mathrm{xa}\left(i\right)\ge \mathrm{ya}\left(i\right)\\ y\left(i\right),\mathrm{otherwise}\end{array}.$

[0043] To further reduce the processing complexity, the second rotation can be updated by the instant phase difference, using

[00002]$R\left(n\right)=\{\begin{array}{c}R\left(n\right).Math.\left(x\left(n\right).Math.{y\left(n\right)}^{*}\right),\mathrm{if}\mathrm{xa}\left(n-1\right)\ge ya\left(n-1\right)\mathrm{and}\mathrm{xa}\left(n\right)<ya\left(n\right)\\ R\left(n\right).Math.\left(y\left(n\right).Math.{x\left(n\right)}^{*}\right),\mathrm{if}\mathrm{xa}\left(n-1\right)<\mathrm{ya}\left(n-1\right)\mathrm{and}\mathrm{xa}\left(n\right)\ge \mathrm{ya}\left(n\right),\\ R\left(n\right),\mathrm{otherwise}\end{array}$

as shown in the the corresponding flow chart shown in FIG. 4.

[0044] At this point, while we have presented this disclosure using some specific examples, those skilled in the art will recognize that our teachings are not so limited. Accordingly, this disclosure should be only limited by the scope of the claims attached hereto.