Improved normalization method for fiber optic sensor systems

Abstract: Quadrature signals Q and in-phase signals I output from a pathlength mismatched interferometric sensor array (10) are normalized by sampling the signals to obtain a first set of values Q k = Q 0 B k sin  k corresponding to the sine of the phase shift and a second set of values I k = I 0 B k cos  k corresponding to the cosine of the phase shift where k is a running index of the set of consecutive samples of Q and I. A quantity DQ k =|Q u(k) ) 2 -(Q s(k) ) 2 | is calculated for the quadrature signals; and a quantity DI k = |(I u(k) ) 2 -(I s(k) ) 2 | is calculated for the in-phase signals, where u(k) and s(k) are parameters used to determine indices that select specific sampled signals and k is a running index of the set of consecutive samples of Q and I. The normalization factor is then calculated as the ratio R = Q N / I N where and

Description: Background of the Invention This invention relates generally to signal processing techniques for fiber optic sensor systems. This invention relates particularly to processing signals indicative of a phase shift in an optical signal in both frequency-division-multiplexed (FDM) and time-division-multiplexed (TDM) sensor systems. Still more particularly, this invention relates to processing signals proportional to the sine and cosine of the phase shift in an optical sensor system to determine the magnitude of the phase shift. Fiber optic sensor systems as disclosed e.g. in US 5227857 of both the FDM and TDM varieties acquire in the demodulation process one term proportional to the sine of the sensor phase shift and another term proportional to the cosine of the phase shift. The sine of the sensor phase shift is referred to as the quadrature term, Q; and the cosine of the sensor phase shift is referred to as the in-phase term, I. The angle of the phase shift is determined by calculating ratio Q/I, which is the tangent of the sensor phase shift. The amplitudes of the sine and cosine terms must be set equal by a normalization procedure to ensure the successful implementation of an arctangent routine to find the sensor phase shift. Summary of the Invention An improved normalization method is proposed that does not use peak detection. It can be used to update the normalization ratio for each sensor channel on an ongoing basis. The acoustic environment can be relatively quiet. Acoustic tones with peak-to-peak amplitudes of about 0.5 radian are sufficient. In the absence of these tones phase spinning at a few hertz due to static pressure or temperature variations is also sufficient. A method for processing signals output from a pathlength mismatched interferometric sensor array that produces phase shifts in sensors included in the array to indicate changes in a parameter being monitored, comprises the steps of sampling the signal output from each sensor in the array to obtain for each sensor a first set of signal terms Q k = Q 0 B k sin  k corresponding to the sine of the phase shift and a second set of signal terms I k = I 0 B k cos  k corresponding to the cosine of the phase shift. For each sensor in the array the method includes the step of calculating a first quantity DQ k =|(Q u(k) ) 2 -(Q s(k) ) 2 | and a second quantity DI k =|(I u(k) ) 2 -(I s(k) ) 2 | where u(k) and s(k) are functions used to determine indices that select specific sampled signals and k is a running index of the set of consecutive samples of Q and L The method further includes the steps of calculating a quantity and calculating a quantity The normalization factor is then determined by the step of calculating the ratio R = QN /IN . The method according to the present invention preferably further includes the step of modulating the phase of optical signals input to at least one of the sensors in the array to normalize the signal terms Q N and I N so that the ratio R=QN /IN is unity for the selected sensor. An appreciation of the objectives of the present invention and a more complete understanding of its structure and method of operation may be had by studying the following description of the preferred embodiment and by referring to the accompanying drawings. Brief Description of the Drawings FIG. 1 illustrates an embodiment of a fiber optic sensor array with which the normalization method according to the present invention may be used; and FIG. 2 is a flow chart of the normalization method according to the invention. Detailed Description of the Invention FIG. 1 illustrates an exemplary sensor array architecture with which the signal processing method of the present invention may be used. The invention is not limited to use with this architecture. The invention may be used with any sensor system that provides two output signals that are in quadrature. FIG. 1 illustrates an FDM architecture 10 that includes a pair of lasers 12 and 14 and three polarization diversity detectors 16-18. The lasers 12 and 14 have corresponding phase modulators 20 and 22. The phase modulators 20 and 22 have phase generated carrier frequencies f 1 and f 2 . The array 10 includes a plurality of Mach-Zehnder interferometer sensors 25-30. The Mach-Zehnder interferometer sensors 25-30 comprise an example of a sensor array with which the present invention can be used. The Mach-Zehnder interferometers 25-27 receive optical signals from the laser 14 via an optical fiber 32. Optical couplers 34 and 36 couple the optical signals from the optical fiber 32 to the Mach-Zehnder interferometers 25 and 26. The Mach-Zehnder interferometer sensors 28-30 receive optical signals from the laser 12 via an optical fiber 40. Optical couplers 42 and 44 couple the optical signals from the optical fiber 40 into the Mach-Zehnder interferometer sensors 28 and 29. The Mach-Zehnder interferometer sensors 28-30 provide output signals to optical fibers 50-52, respectively. The optical fibers 50-52 are arranged to provide optical signals to the detectors 16-18, respectively. The detectors 16-18 are connected to a demodulator 56. The Mach-Zehnder interferometer sensors 25-27 provide output signals to optical couplers 60-62, respectively. The optical couplers 60-62 couple the output signals from the Mach-Zehnder interferometer sensors 25-27 into the optical fibers 50-52, respectively, so that these output signals are also input to the detectors 16-18, respectively. The Mach-Zehnder interferometer sensors 25-30 may be substantially identical. Therefore, only the Mach-Zehnder interferometer sensor 25 is described in detail. The reference numerals for components of the Mach-Zehnder interferometer sensor 25 end with the letter "a." The remaining Mach-Zehnder interferometer sensors 26-30 have reference numerals that end with the letters b-f, respectively. The Mach-Zehnder interferometer sensor 25 includes a pair of optical couplers 70a and 72a each arranged to couple light between a pair of optical fibers 74a and 76a. The optical fibers 74a and 76a include portions wound around corresponding mandrels 78a and 80a, respectively. The coupler 34 couples light from the optical fiber 32 into the optical fiber 74a. The coupler 70a divides the light between the optical fibers 74a and 76a. The coupler 72a combines the light that has propagated in the optical fibers 74a and 76a. The combined light beams are then input to the fiber 50 via the coupler 60. The combined light beams produce an interference pattern that indicates the phase difference between light guided by the optical fibers 74a and 76a between the two couplers 70a and 72a. The phase difference can be used to monitor acoustic pressure, for example, where each sensor is located. The demodulator 56 receives electrical signals from the photodetectors 16-18 and then produces electrical signals that indicate the two terms Q and I that are processed to determine the phase shift for each of the sensors 25-30. There are several techniques for obtaining electrical signals that indicate the terms Q and I that are used to determine the phase shift for each sensor in a fiber optic sensor array. One suitable technique involves sinusoidally modulating or dithering the optical signal output by the optical signal source. The terms Q and I may then be determined by synchronous detection of an odd harmonic and an even harmonic of the modulation or dither frequency. Modulation of the optical signal may be accomplished by modulating the current drive to a semiconductor laser. Another technique involves external phase modulation of the laser output. External phase modulation is preferred when low noise is a significant consideration. In a mismatched path fiber interferometer, the optical frequency modulation of the source light becomes a "phase-generated carrier" modulation that the optical receiver converts into an electrical signal. A pair of voltage sources 84 and 86 are connected to the phase modulators 22 and 20, respectively. Voltage sources 84 and 86 are connected to demodulator 56 via corresponding feedback loops. Using the feedback, the voltage applied to the phase modulator 22 by the voltage source 84 is adjusted using the algorithm described below and illustrated in the flow chart of FIG. 2 to control the modulation depth of the phase modulator 22 such that the ratio Q/I is 1.0 for one selected sensor, for example, the sensor 25. The voltage applied to the phase modulator 20 by the voltage source 86 is similarly adjusted to control the modulation depth of the phase modulator 20 such that the ratio Q/I is 1.0 for another selected sensor, for example, the sensor 28. If the other sensors 26, 27, 29 and 30 are similar to the sensors 25 and 28, then adjusting the phase modulator voltages to normalize the signal terms Q and I for the sensors 25 and 28 will ordinarily cause the signal terms Q and I from all sensors in the array to be approximately normalized. The signal terms Q and I for the remaining sensors in the array are then completely normalized using the algorithm illustrated in the flow chart of FIG. 2. The quadrature and in phase terms for signal processing are considered as discrete samples in the digital domain. The quadrature and in-phase terms for two samples designated 1 and 2 are given by Q 1 =Q 0 B 1 sinϕ 1 ; I 1 = I 0 B 1 cosϕ 1 ; Q 2 = Q 0 B 2 sinϕ 2 ; and I 2 = I 0 B 2 cosϕ 2 . In Eqs. 1-4, the phase shift to be measured is ϕ, and the signal amplitudes that change over time due to polarization fading are Q 0 B and I 0 B. The coefficients Q 0 and I 0 are multiplicative constants whose ratio is needed for signal normalization. Over longer times the ratio of Q 0 /I 0 varies slowly due to changes in the path mismatch of the interferometer sensor. The fading term B can be considered a constant during a sufficiently short time span between samples 1 and 2. In this limit Eqs. 1-4 can be solved to obtain the ratio R = Q 0 /I 0 , which may be written as R= Q 0 I 0 = |Q 2 1 -Q 2 2 | |I 2 1 -I 2 2 | In addition to the assumption of stationary B, there is the assumption of some movement in ϕ to obtain the ratio R. FIG. 2 illustrates a more general algorithm that is described below. First, a quantity DQ k is defined as: DQ k =|(Q u(k) ) 2 -(Q s(k) ) 2 |; and a quantity DI k is defined as: DI k = |(I u(k) ) 2 - (I s(k) ) 2 |. The index k is a running index of the set of consecutive samples of Q and I. Expressions u(k) and s(k) are used to determine indices that select specific sampled signals. Any Q and I in the sampled data may be used to determine the difference terms DQ k and DI k . Some examples of u(k) and s(k) with the constraints N 1 ≤ k ≤ N 2 and N 1 ≥ 0 where N 1 and N 2 are integers used to label the sampled data are listed below. (a) u(k)=(k+s) 2 and s(k)=k2 , with s being an integer constant. (b) u(k) =2 k +(-1) k and s(k) = N 1 . (c) u(k)=u·k+t and s(k)=s·k, with u, s and t being integer constants. (d) u(k)=s·k+t and s(k)=s·k with t being the spacing between the signals determined by s·k. A value is then obtained for Q N as follows: Similarly, a value is obtained for I N : The normalization factor relating the amplitude terms or the expressions for the quadrature and in-phase terms may then be calculated as: R = Q N I N Usually only a small percentage of all samples of Q and I need be used. Using the parameters u(k) and s(k) as defined in example (d) above, differences are taken of squared terms t samples apart. Sums of these differences s samples apart are generated in estimating Q 0 /I 0 . Example (d), which is a special case of example (c) has been used in a simulation of the normalization method according to the present invention. Using u(k) and s(k) as defined in example (c), the general expressions of Equations (6) and (7) become: DQ k =|(Q u·k+t ) 2 -(Q s·k ) 2 | and DI k = |(I u·k+t ) 2 -(I s·k ) 2 |. For example (d), Equations (11) and (12) become DQ k =|(Q s·k+t ) 2 -(Q s·k ) 2 | and. DI k = |(I s·k+t ) 2 -(I s·k ) 2 |, Simulations of the normalization method have been performed with Math Cad with the index u = s where s is an integer. The model uses 32768 samples at a sample rate of 512,000/s. Polarization fading is incorporated into the model with drifts in the signal amplitude of about 15 percent over a run of 32768 samples lasting 64 ms. In the first trial, a small acoustic tone with an amplitude of 0.4 radians peak-to-peak at 500 Hz is introduced with no spinning phase due to static pressure or temperature variations. The algorithm outlined in for the special case of example (d) was used with s = 100 and t = 200 and was repeated 301 times to form the appropriate sums. In all, 303 samples containing 303 sine and 303 cosine terms were used out of the total of 32768 samples. This result was an estimate for Q 0 /I 0 that was in error by only 0.5 percent. In the second trial, a large acoustic tone of 100 radians peak-to-peak at 500 Hz was introduced. Integers s and t have the same values and only 53 samples were used. In this case the estimate for Q 0 /I 0 was in error by only 0.16 percent. The simulations demonstrate the ability of the algorithm to estimate the normalization ratio to within a few tenths of one percent for small and large acoustic tones. Additional simulations were done at the lower sampling rate of 50 thousand samples per second with a total run time of 655 milliseconds and a drift in signal amplitude of 15 percent due to polarization fading. Acoustic tones of only 0.6 radians peak-to-peak were introduced at a range of different frequencies from 1 kHz to 10 kHz. The worst case estimate for Q 0 /I 0 was in error by 3 percent. In the vast majority of cases the error was less than 0.5 percent. Integer s is set to 10 and t is set to 3 and a total of 202 samples was used in all cases. The same acoustic tones are introduced again in the presence of phase spinning at one Hz. With s set at 300 and t set at 150, a total of 202 samples was used. The worst case estimate for Q 0 /I 0 was in error by only 1.5 percent when the tone amplitude was set to zero. These examples show that the algorithm parameters s, t, and N 1 - N 2 can be varied to suit the operating conditions encountered.

Claims: A method for processing signals output from a pathlength mismatched interferometric sensor array (10) that produces phase shifts in sensors (25-30) included in the array to indicate changes in a parameter being monitored, comprising the steps of: sampling the signal output from each sensor (25-30) in the array to obtain for each sensor (25-30) a first set of signal terms Q k = Q 0 B k sin  k corresponding to the sine of the phase shift and a second set of signal terms I k =I 0 B k cos  k corresponding to the cosine of the phase shift; calculating a first quantity DQ k =|(Q u(k) ) 2 -(Q s(k) ) 2 |and a second quantity DI k =|(I u(k) ) 2 -(I s(k) ) 2 |where u(k) and s(k) are parameters used to determine indices that select specific sampled signals and k is a running index of the set of consecutive samples of Q and I; calculating a quantity calculating a quantity and calculating the ratio R = Q N /I N . The method of claim 1 further including the step of modulating the phase of optical signals input to a selected one of the sensors (25-30) in the array (10) to normalize the signal terms Q N and I N so that the ratio R = Q N / I N is unity for the selected sensor (25-30).

Inventors: HALL DAVID B, HALL, DAVID B.
Applicants: LITTON SYSTEMS INC, LITTON SYSTEMS, INC.
Application year: 2000
Main classification: G01D5
IPC codes: G01D5