GB2601639A - High speed sampling of sensors - Google Patents

Abstract

The method involves receiving a set of bursts of samples of an optical signal, where each burst of samples comprises a set of samples of each pulse in a group of discrete optical pulses (A1-E1, A2-E2), groups are regularly spaced in time with a group repetition period larger than group width, and the optical pulses comprise a phase modulated signal. The bursts of samples are processed to obtain a series of regularly spaced samples. The regularly spaced samples are demodulated to retrieve phase information from the detected pulses. The processing step comprises the steps of down-sampling the received samples, up-sampling the down-sampled samples at a rate of the group width multiplied by the sample rate in each burst, then filtering the up-sampled samples and combining the filtered samples. The method increases the effective sampling rate to improve sensor bandwidth and enables locking cycles within a pulse for phase-locked loop to improve locking performance.

Description

HIGH SPEED SAMPLING OF SENSORS

Background

[0001] Optical sensors provide a convenient means for remote sensing systems, and enable the monitoring of locations large distances from an interrogation location.

[0002] Figure 1 shows a schematic diagram of an optical sensing system comprising four sensing regions X, Y, Z, P. The optical system comprises a main optical fibre 100 through which light propagates. Each sensing region produces an optical phase change in light propagating through the sensor region in response to a change in the natural phenomena being monitored. For example, a hydrophone may detect changes in pressure or accelerometers may detect a change in velocity. Typically such sensors produce a change in physical length of the optical fibre and hence an optical phase change at the output of the sensor, but other mechanisms may be utilised.

[0003] In the system of Figure 1, optical couplers 101a, 101b, 101c, 101d are provided between each sensor region to couple a portion of the light propagating in the main fibre 100 into a branch fibre 102a, 102b, 102c, 102d. The branch fibres are each terminated by a mirror 103a, 103b, 103c, 103d. A mirror 103e terminates the end of the main optical fibre 100.

[0004] An interrogation system 104 comprises at least one light source coupled to transmit light into the main optical fibre 100, and at least one optical receiver coupled to receive light returning from that fibre. The interrogation system may be provided in a distributed manner such that elements are located in different locations. For example, it may be desirable position the optical receiver and sampling systems in proximity to the optical sensors (for example on the sea bed) and to transmit the data to a processing location in a convenient place for analysis.

[0005] Light reflected by pairs of mirrors either side of each sensor region interferes at the receiver, with the power of the resulting signal being dependent on the optical phase difference between the returning signals. That phase difference is dependent on the optical phase length of the relevant sensor region, and hence will change depending on that phase length. Changes in the detected power can thus be used to detect changes in sensor phase length and hence the physical phenomena to which the sensor responds.

[0006] The output of such optical sensors is cyclical with the optical phase length and thus the change between subsequent samples that can be unambiguously resolved is limited to +/-7 radians.

[0007] Figure 2 shows an example of a series of 8 samples from an optical sensing system. Figures 3a and 3b show these on a phase circle. At sample #0 the phase is assumed to be 0 as the starting point for the measurement. At samples #1 -#3 the phase increases in relatively small increments and can be tracked unambiguously. Between sample #3 and #4 the measured phase "wraps" around from +0.957 to -0.37 (shown by the hollow circle in Figure 2). However, by making the assumption that the phase never changes by more than 7 between samples, we can unambiguously resolve the actual total phase change and infer that the phase has continued to increase to +1.77 in relation to the starting phase. Sample #5 shows a measured phase of +0.657, and by making the above assumption the total phase change can be unambiguously resolved to 2.657.

[0008] Sample #6 shows a sampled phase of -0.187. Applying the assumption that the phase change is less than 7 would result in the calculated overall phase change decreasing to +1.827 (point 201). However, in fact, the phase may have continued to increase and has actually increased to +3.827 (point 202). There is therefore an error of 27 in the overall phase change, and the direction of change has been calculated incorrectly. The error between samples 5 and 6 can never be recovered because it is only by continuous tracking that the absolute change can be resolved.

[0009] Optical sensors thus have a limited range between samples which limits the rate of change of phase in the output signal that can be detected unambiguously to less than 7/T5, where T, = sample period. The limit that this applies to the rate of change of the sensed phenomena will depend on the characteristics of the sensor, but the rate is limited.

[0010] There is therefore a need for an optical sensing system capable of resolving larger rates of change.

Summary

[0011] This Summary is provided to introduce a selection of concepts in a simplified form that are further described below in the Detailed Description. This Summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used as an aid in determining the scope of the claimed subject matter.

[0012] Systems and methods for interrogating sensing systems utilising bursts of samples. Bursts of samples correspond to optical pulses returning from optical sensors, where pulses are spaced at a period significantly longer than the pulse width, giving irregular sample spacing. The interrogation system and method processes the irregular busts of samples to recover phase information from received signals.

[0013] There is provided a method of interrogating an optical sensing system performed at an interrogation system, comprising the steps of receiving a plurality of bursts of samples of an optical signal, each burst of samples comprising a plurality of samples of each pulse in a group of discrete optical pulses, wherein the groups are regularly spaced in time with a group repetition period larger than the group width, and the optical pulses comprise a phase modulated signal; processing the bursts of samples to obtain a series of regularly spaced samples; and demodulating the regularly spaced samples to retrieve the phase information from the detected pulses. [0014] The samples may be processed utilising a filter bank.

[0015] The filter bank may comprise a number of parallel paths equal to the number of samples in each burst.

[0016] The step of processing may comprise the steps of downsampling the received samples, upsampling the downsampled samples, filter the upsampled samples, and combining the filtered samples.

[0017] The downsampling ratio may be equal to the sample rate in each burst divided by the group repetition frequency.

[0018] The upsampling may be performed at a rate of the group width multiplied by the sample rate in each burst.

[0019] There is also provided a method of interrogating an optical sensing system performed at an interrogation system, comprising the steps of receiving a plurality of bursts of samples of an optical signal, each burst of samples comprising a plurality of samples of each pulse in a group of discrete optical pulses, wherein the groups are regularly spaced in time with a group repetition period larger than the group width, and the optical pulses comprise a phase modulated signal; calculating the instantaneous frequency of the signal represented by each burst of samples; and recovering phase information of the detected pulses.

[0020] The step of calculating may comprise filtering the samples utilising an FIR filter having a length less than the number of samples in each burst of samples. [0021] There is also provided a method of interrogating an optical sensing system performed at an interrogation system, comprising the steps of receiving a plurality of bursts of samples of an optical signal, each burst of samples comprising a plurality of samples of each pulse in a group of discrete optical pulses, wherein the groups are regularly spaced in time with a group repetition period larger than the group width, and the optical pulses comprise a phase modulated signal; phase locking a signal to the samples of each burst of samples; and utilise the phase locked signal to recover phase information of the detected pulses.

[0022] Each group of pulses may consist of a single pulse.

[0023] Each burst of samples may comprise 8 to 80 samples.

[0024] The samples within each burst may be at a sampling frequency of 40 to 400 MHz.

[0025] The group width may be 200ns.

[0026] The group repetition period may be in the range of 0.5 jts to 10 ps.

[0027] The regularly spaced samples may have a sample rate in the range of 10MHz to 20MHz.

[0028] The method may further comprise the steps of detecting the optical pulses, and sampling the pulses to obtain the bursts of samples, and transmitting the bursts of samples to the interrogation system.

[0029] The steps of detecting, sampling, and transmitting may be performed at a different location to the other steps of the method.

[0030] The optical pulses may be returning from an optical sensor.

[0031] The method may further comprise the steps of detecting other optical pulses returning from other optical sensors, wherein the optical pulses returning from different sensors are time multiplexed.

[0032] The phase change between adjacent samples of the regularly spaced samples may be less than 7E.

[0033] There is also provided an interrogation system configured to perform the methods described herein (e.g. an interrogation system configured to perform any of the methods described above).

[0034] The preferred features may be combined as appropriate, as would be apparent to a skilled person, and may be combined with any of the aspects of the invention.

Brief Description of the Drawings

[0035] Embodiments of the invention will be described, by way of example, with reference to the following drawings, in which:-Figure 1 shows a schematic diagram of an optical sensing system; Figure 2 shows an example output of an optical sensor; Figures 3a and 3b show phase plots relating to Figure 2; Figure 4a shows an example of pulses transmitted into a system to interrogate the sensors, Figure 4b shows an example of pulses returning from a sensor of the system of Figure 1; Figure 5 shows an example of pulses returning from the system of Figure 1; Figure 6 shows a frequency chart; Figure 7 shows a series of returning pulses from a sensor system; Figure 8 shows the pulses of Figure 7 sampled at a high rate; Figure 9 shows a single sensor's pulses from Figure 8; Figure 10 shows a reconstructed signal; Figure 11 a shows a filter bank for processing irregular samples; Figure 11 b shows an example phase demodulation apparatus; Figure 12 shows a filter bank for designing the filters utilised in the filter bank of Figure 11a; Figure 13 shows a system of processing data utilising instantaneous frequency; Figure 14 shows a chart of instantaneous frequency; Figure 15 shows a further system utilising instantaneous frequency; Figure 16 shows a chart displaying extrapolation techniques; Figure 17 shows a chart displaying PLL behaviour; and Figure 18 shows a PLL system.

Common reference numerals are used throughout the figures to indicate similar features

Detailed Description

[0036] Embodiments of the present invention are described below by way of example only. These examples represent the best ways of putting the invention into practice that are currently known to the Applicant although they are not the only ways in which this could be achieved. The description sets forth the functions of the example and the sequence of steps for constructing and operating the example. However, the same or equivalent functions and sequences may be accomplished by different examples.

[0037] Figure 4a shows a series of optical pulses which may be used to interrogate the optical sensor system of Figure 1. Each pair of pulses is timed such that a reflection of the first pulse from the second of a pair of mirrors interferes with a reflection of the second pulse of the pair from the first of that pair of mirrors at the interrogation system. That is, a reflection of pulse Al from mirror 103b returns to the interrogation system at the same time as a reflection of pulse A2 from mirror 103a. The delay between pulses Al and A2 is thus equal to the difference in round-trip time from the interrogator to mirrors 103a and 103b, 103b and 103c, 103c and 103d, 103d and 103e. A single overlapping pulse is thus received at the receiver for each sensor region X, Y, Z, and P. As the optical phase length of the sensor regions changes the result of the interference between Al and A2 changes, which can be used as explained above to determine changes in the physical phenomena detected by the sensors.

[0038] Subsequent pairs of pulses 131, 82 etc can be transmitted into the fibre after an appropriate delay such that the final useful reflection from the first pair of pulses has returned to the interrogation system before the first useful reflection of the second pair of pulses. The delay between pairs of pulses defines the optical sampling period Ts (T=1/F, where Fs is the optical sampling frequency) for each sensor region.

[0039] Sensors X, Y, Z, and P are located sequentially along the main fibre 100 and accordingly can be interrogated using a single pair of pulses, with reflected pulses being time-multiplexed in the returning signals. Figure 4b shows the return signal from a single pair of pulses (e.g. Al and A2) being launched into the sensor system. The first return pulse is the reflection of pulse Al from mirror 103a, and the last pulse is a reflection of pulse A2 from mirror 103e. The central pulses are pairs of overlapping pulses reflected from the pairs of mirrors 103a/103b, 103b/103c, 103c/103d, and 103d/103e. Each of these central pulses contains information on the optical phase length of the sensors X, Y, Z, and P. [0040] In a typical system the pulse width, Wp, may be 200ns, with a spacing of 200ns between pulses of a pair. The spacing between pairs of pulses varies strongly depending on the length (and hence round-trip time) of each system. In Figure 4a this is lgs which is sufficient for the 4 sensors.

[0041] The pairs of pulses may be at different optical frequencies in a heterodyne system, or the same frequency in a homodyne system. In general a heterodyne system is preferred to avoid difficulties with zero signals in a homodyne arrangement. The following description will concentrate on a heterodyne system, but the principles are also applicable to homodyne systems.

[0042] Figure 5 shows an example of the signals received from transmitting the signals shown in Figure 3 into the system of Figure 1. Each bar represents a returning pulse from the system. The hatched bars represent non-useful single-pulse reflections (e.g. A1/B1/C1/D1/E1 from mirror 103a and A2/B2/C2/D2/E2 from mirror 103e), as described above with reference to Figure 4. In Figure 5, the timing is such that the final non-useful pulse from one pulse pair coincides with the first non-useful pulse from the subsequent pulse pair (e.g. A2/103e coincides with B1/103a etc.). This minimises the percentage of unusable time.

[0043] The dashed line shows the signal interpolated between the signals returning from sensor Z. The signal has the form sin(2ft(f2-fi)t + IMO), where f2 is the optical frequency of the second pulse in each pair (A2, B2, 02 etc.) and fl is the optical frequency of the first pulse in each pair (Al, B1, Cl etc.), and (1),(t) is the optical phase shift of sensor z with time, t. The resulting signal is thus a phase-modulated sine wave with a carrier frequency Fe = (f241). The detected pulses represent samples along that signal, with an optical sampling period of 1/Fs defined by the spacing between pairs of pulses, as described with reference to Figure 3.

[0044] It is attractive to increase the number of sensors which are time-multiplexed. A dense system may therefore instead have a pulse-pair spacing of 5,,is (giving an optical sample period of T5=51is and an optical sample frequency of Fs=200kHz). Such a dense system would allow 4 sensor units similar to that of Figure 1 and an additional single reference sensor to be interrogated by appropriately multiplexing them with delays to give 17 sensors in total.

[0045] In accordance with the Nyquist sampling theorem, samples at a sampling frequency of Fs can only unambiguously recover signals in a bandwidth of Fs/2, if there are no other constraints on the signal. The rate of change of (t)must thus be limited such that the total bandwidth of sin(2ft(f2-fi)t + IMO) is less than Fs/2. That is:-F0 + (-1/27)4/dt < F5/2.

[0046] This is equivalent to the explanation above of a limit of a ft radian phase change between samples. The instantaneous frequency is defined F, = F0+(-1/2rOddildt.

[0047] Figure 6 shows the frequency domain with the carrier signal Fs positioned at F314, which provides the maximum use of the 0 -F5/2 bandwidth having equal bandwidths of Fs/4 for positive and negative frequency shifts. Beyond +/-Fs/4 (0 to Fs/2) aliasing occurs and the original signal cannot be unambiguously recreated. In the example of Figure 5, Fs is 1MHz, but for a dense system a typical Fs may be 200kHz, giving a bandwidth of +/-50kHz.

[0048] The impact of this is that the carrier frequency Fe=f2-f, and rate of change of clpz is limited by the optical sampling rate. In turn the upper limit of the optical sampling rate is limited by the period required between pairs of pulses which is set by the optical length of the system to avoid overlap between signal returning from the first and last sensors. For a given system there is therefore a maximum achievable pulse repetition rate, which gives a maximum sample rate, which in turns defines the maximum achievable bandwidth of each sensor.

[0049] Figure 7 shows optical power received at the interrogation unit in a system of Figure 1 in which the instantaneous frequency, FJ, is much larger than discussed above. Due to the higher instantaneous frequency the optical power varies significantly within the period of the returning pulse. Figure 7 shows an optical signal with an instantaneous frequency of 4MHz. Hence 1/F; is only slightly longer than the pulse width (a 250ns cycle against the 200ns pulse). In an example, Fp may be defined as 3MHz, which has a cycle of 333ns. Under static conditions, FFFo, so a 200ns pulse would contain around 60% of one cycle. Negative phase ramps will result in an increase in F; (e.g. up to close to 6MHz) and positive phase ramps will result in a decrease in F, (e.g. down to close to 0Hz). When a phase change is applied to that signal the instantaneous frequency may increase to 4MHz as shown in Figure 7. In example systems the unmodulated carrier frequency, Fd, may be such that each optical pulse may contain > 50% of a cycle or > a whole cycle.

[0050] In the conventional sampling system described above there would be a single digital sample for each optical pulse such that the digital sample rate matches the optical sample rate; the digital sampling rate would not be high enough to recover the signal and information would be lost.

[0051] Figure 8 shows the signal of Figure 7 digitally sampled with an analogue to digital converter at a higher rate than in previous systems.

[0052] Figure 9 shows the sampled signals of Figure 8 for only a single sensor. As shown in this figure the digitally sampled data comprises discontinuous bursts of samples, with significant periods with no information on the signals. In this example system the pulse width Wp=200ns with a pulse repetition period of 11.ts. In the example shown, a digital sample rate Fd=60MHz provides 12 samples within each pulse. There are thus M=12 samples unevenly distributed over a period of 1tis. Conventional sampling and processing systems utilise regular sample spacing to recover the input signal. However, special processing techniques may be utilised to recover the signal from this form of data. Processing the discontinuous samples of Figure 9, as described below, yields a set of regularly spaced samples at a rate equivalent to the average sample rate of Figure 9, as shown in Figure 10. These samples can then be utilised to unambiguously recover the signal, provided the criteria discussed below are met. The effective sampling rate of this technique is much higher than in the conventional system described above and accordingly the sensor's bandwidth is improved.

[0053] The signal of Figure 10 is sin(2-rr(f2-fi)t + 07(t)) and thus the sensor's phase change can be extracted by demodulating the signal. The average sampling frequency is given by Fa, = M.F, = Wp. Fd. Fs [0054] In the example discussed above 12 samples are obtained every 1p.s, giving an average sample rate in Figure 10 of Fa=12MHz. The Nyquist theorem applies to the average sampling rate, and therefore the bandwidth for perfect reconstruction is 0-6M Hz.

[0055] As above, selecting the centre of the bandwidth as the carrier frequency, Fc, maximises the bandwidth available for phase changes, and thus Fe may be selected as 3MHz, giving a bandwidth of +1-3MHz for phase change measurements.

[0056] As previously described it is desirable to increase the multiplexing to allow a dense system of 17 sensors by increasing the optical sample period to 5,is. This may be achieved by increasing the digital sample rate to Fd=300MHz so that there are 60 pulses within each returned signal. With a 51.1s optical sample period and 200ns pulse width, this corresponds to an average sampling rate of F5v=12MHz as before. Thus the same carrier frequency of Fe=3MHz can be used allowing a range of +1-3MHz as before. This compares favourably with a bandwidth of +1-50kHz in previous dense systems.

[0057] In certain circumstances it may be possible to reduce the sample rate by making assumptions about the input signal. For example, it may be possible to assume that the input signal has the form sin(2rr(f2-fl)t + (PAO), with the amplitude and frequency of (k assumed to be limited to known ranges.

[0058] Various methods are known for converting a discontinuous set of samples into regularly spaced samples. Recovery of the original signal then requires conventional processing. Figure 11 shows a schematic diagram of a reconstruction system utilising a filter bank to convert the set of discontinuous samples w[n] into a regular set of samples x^[n]. Delay elements 'T' provide a delay equal to the input digital sample period, 1/Fd, such that each path of the filter bank starts processing from a different sample. The number of paths in the bank is equal to the number of samples in each burst. The signal of each path is first downsampled (by N=Fd/Fs) to select one sample from each burst. Due to the delay elements T each path processes a different sample of the bursts, such that all samples are processed by one path. The signals are then upsampled (by M=Wo.Fd) to the average sampling rate, Fav, to produce a regular set of samples. The regular samples are then filtered by a respective filter Fo(z) (z) and the paths are combined to provide a regular set of samples representing the signal of one sensor, which can be processed to obtain the phase data for that sensor. As explained above the bandwidth is significantly increased compared to previous sensors and hence the sensor system bandwidth is also increased.

[0059] In summary, a series of bursts of samples are obtained, each burst corresponding to one of a series of detected, regularly spaced, optical pulses. The bursts of samples are converted to a regular series of samples which may then be demodulated to retrieve a signal with which the detected pulses were modulated.

[0060] There are a number of conventional techniques which can be applied to demodulate the phase once it has been transformed to a uniform sampling pattern. One technique, shown in Figure 11b, is to split the signal into two and generate inphase I and quadrature Q signals. The I signal may be generated by multiplying (110) by sin(2rrFct) and low pass filtering (111) to remove components of 2Fc and above. Similarly the Q signal may be generated by multiplying (112) by cos(27Fct) and low pass filtering (113). In an example Fc is locked to the oscillator that generates Fi and F2 for example, the same oscillator is utilised. The phase is then demodulated using a four quadrant arc tangent function (114) as criz(t)=arctan(Q/1). The demodulated phase is then unwrapped (115). However other phase demodulation schemes (e.g. Hilbert Transformer, phase-locked loop) are equally applicable. Once the phase has been demodulated and unwrapped it can be decimated to the required output sample rate. It can be converted to the correct units by multiplying by an appropriate calibration coefficient. The paper MuItirate Filter Bank Reconstruction of Signals from Bunched samples, Conference record of the Thirty-Sixth Asilomar Conference on Signals, Systems and Computers, 2002, Volume 1, Pages 781-785, provides a method of calculating the filter coefficients Fo(z) Fm_1(z) in Figure 11 for near-perfect reconstruction of signals and is incorporated herein in its entirety. Briefly, the input w[n] is modelled by an interpolator which consists of the cascade of an upsample by P operation followed by interpolation filter Hiat(z), where P=M/N. The interpolation filter is chosen to have a normalised cut-off of irr/2 rad. It is chosen as a compromise, such that it has sufficient order to give a reasonable approximation to x[n] where x[n] is the true input signal at sample rate Fav. However, the order should be minimised to reduce the complexity of the resulting synthesis filters. The interpolation filter, H(z), is then transformed to be in polyphase form in terms of zP.. The Noble Identity is then used, to swap the cascade order, such that the polyphase interpolation filter (now in terms of z) precedes the downsample by P operation. Once this has been done, the downsample by P is combined with downsample by N, to yield a single downsample by M. The transformed interpolation filter coefficients are combined with the delay terms to generate filters Fo(z) FR/1_1(z). These transformations result in the maximally decimated filter bank shown in Figure 12.This filter bank takes the actual input signal x[n], hypothetically sampled at rate Fav, and generates the approximated signal xA[n] sampled at Fa,. Near perfect reconstruction of the whole system is achieved if xA[n] = K.x[n-d] [0061] where K is a constant scaling factor and d is a delay. A matrix approach is then used to solve this equation for the maximally decimated filter bank of Figure 12, to determine the filters Fo(z) Fm_i(z). These filters can then be used for the implementation of the system of Figure 11. A more general technique can be found in Reconstruction of Bandlimited Periodic Nonuniformly Sampled Signals through Multirate Filter Banks, IEEE Transactions on Circuits and Systems I: Regular Papers, Volume 51, Issue 8, Prendergast, RS., and is incorporated herein in its entirety; further techniques are disclosed in references [7], [8] and [9] of this work. Using this more general technique it is possible to cover situations whereby the samples are more complicated than a periodic bunched pattern. For example, rather than have a simple pair of pulses, each pulse within the pair could themselves be divided up. One configuration would be to have the pulse at fi separated into a pulse 100ns long, followed by a 100ns gap, followed by a 100ns long pulse. Similarly the pulse at f2 is separated into a pulse 100ns long, followed by a 100ns gap, followed by a 100ns long pulse. Now every sample period, Ts there are two bunches of samples rather than a single bunch. E.g. for fd=300MHz, there will be 30 samples, followed by a 100ns gap followed by another 30 samples, then a gap of 4700ns and then repeating. Within the gaps another string of sensors could be inserted, by appropriately multiplexing them in with delays.

[0062] The pulse signals transmitted into the system may be generated in any conventional manner. For example, a known approach is to use an Acousto-Optic Modulator (AOM) which modulates both amplitude and frequency of light. When the AOM is "on" to allow light to pass the output optical frequency is shifted by the SF frequency used to drive the modulator. By varying that frequency different frequency shifts can be applied. For example, the two pulses used in the examples may be generated using an AOM driven at 198.5MHz for the first pulse and 201.5MHz for the second pulse. As will be appreciated by the reader, the above description has been given in relation to a limited number of sensors, but the concepts described can be expanded to apply to any number of multiplexed sensors and light sources. For example, it is known to utilise an interrogator system to generate multiple wavelength signals, with different wavelengths being directed to different sensors by wavelength selective couplers in the sensor system. Similarly the receiver system at the interrogator utilises known techniques and configurations, but with a higher sampling frequency to allow the processing described above.

[0063] In alternative systems the discontinuous signal may be processed in other ways to recover the signal. Figure 13 shows a system for determining the instantaneous frequency, and calculating a phase value. Figure 14 shows the amplitude of the signal and its instantaneous frequency. The signal at the interrogation system comprises bursts of high frequency carrier signal, phase modulated by the optical sensor regions. The instantaneous frequency of each signal burst contains information on the optical phase change during that burst and can be utilised to recover the phase signal. The signal is sampled at block 140 at a high digital sampling frequency, Fd (e.g. Fd=300MHz) and a set of M samples for each pulse gathered in buffer 141. The samples are passed through a Matched FIR filter 142 to give the real part and a Hilbert FIR filter 143 to give the imaginary part, each filter having a length less than M. The filters 142 and 143 are designed so that they have substantially the same spectral amplitude response, but with filter 143 having a -Tr/2 phase difference across the passband. The outputs of these filters are used to calculate the instantaneous frequency, Fi, at block 144:-M-1

F = Arg * ;Zn+1 -7 n n=1 [0064] Where* denotes the complex conjugate, Zd and Zni are consecutive complex samples arising out of the Hilbert FIR filter/matched filter combination, and Arg is the angle of the resulting complex number. The signal is low-pass filtered at block 145. It is then integrated at block 146 to recover the phase using 0 = -27 Fidt [0065] Other processing operations may be performed such as decimation and multiplication by sensor calibration coefficients.

[0066] Figure 15 shows a further system for measuring the instantaneous frequency and recovering the phase information. The system generates a set of complex samples for each pulse as before. For each pulse, the instantaneous frequency is estimated in 144 as before. In 145 the instantaneous phase of each sample is calculated and the average of all samples within the pulse taken. The phase unwrap block 146 utilises an estimate of the actual phase from the previous phase and frequency measurement to unwrap the actual total phase change, as shown in Figure 16. From sample 1 the instantaneous frequency (160) and phase (161) is calculated and used to extrapolate the phase (162) to the nearest 7 radians, using =-2-rrFiTs [0067] where A(1) is the phase change between samples, F, is the instantaneous frequency and Ts is the optical sample period. The estimated phase change is used to determine the number of 2rrs that have been traversed between samples and thus unwrap the phase obtained in sample 2 to calculate the actual phase 163.

[0068] There are other ways in which the measure of instantaneous frequency may be utilised in order to extend the allowable dynamic range. For example in W02008/110780 incorporated herein in its entirety, there is a method disclosed in which measurements of the instantaneous frequency are combined with the measured phase by stitching together segments in the time domain. A determination is made to see over what time periods the phase measurement is overscaled' (i.e. aliased). During periods of time in which the phase measurement is overscaled', an integrated version of the instantaneous frequency is substituted. In this way a signal is produced which is not overscaled, but has good signal to noise ratio.

[0069] A Phase-Locked Loop (PLL) system is an alternative method which may be utilised to recover the signal. A PLL may be used which can lock to the sinusoidal signal which is detected within each pulse. The PLL must be capable of acquiring and locking to the received signal within the pulse width. In between pulses the PLL 'free-runs' at the same frequency at which it locked to. VVhen it reaches the next pulse, it re-locks to the sinusoidal signal within that pulse. During the free-run' period, the PLL is effectively tracking through multiple 2rr phase changes. Simulation has shown that a PLL is capable of locking to a signal within 2 cycles. This requires that Fs is higher than in previous examples such that there are sufficient cycles of the carrier within each pulse for the PLL to acquire and lock to the signal. Simulation, as shown in Figure 17 has demonstrated that a PLL can acquire lock within 2 cycles. A further 2 cycles within a pulse for the PLL to remain locked improves performance. Based on the previously discussed pulse width of 200ns, this suggests an Fs of 20MHz. Utilising a high speed sampling rate of 200MHz, and a pulse repetition frequency of 200kHz simulations show that a PLL with a natural frequency fn of 8MHz and damping factor c=0.9 will settle within 20 samples. The phase output of the PLL can be used as a direct measurement of phase, or the frequency can be integrated to provide the measurement.

[0070] Figure 18 shows a schematic diagram of a PLL that may be utilised in this system. The exemplary PLL utilises a Hilbert phase detector 180 which can provide an instantaneous phase error without subsequent filtering (which would cause time delay) to provide a very fast lock time. The PLL of Figure 18 is given for example only and any PLL with appropriate characteristics can be utilised.

[0071] The above description has been given with respect to a particular sensor system and interrogation method. However, as will be appreciated the processing and analysis techniques are applicable to any sensor system or interrogation method in which a phase modulated signal must be demodulated to recover sensor data. In particular, the systems and methods described herein are not restricted to optical sensing systems.

[0072] Any range or device value given herein may be extended or altered without losing the effect sought, as will be apparent to the skilled person.

[0073] It will be understood that the benefits and advantages described above may relate to one embodiment or may relate to several embodiments. The embodiments are not limited to those that solve any or all of the stated problems or those that have any or all of the stated benefits and advantages.

[0074] Any reference to 'an' item refers to one or more of those items. The term 'comprising' is used herein to mean including the method blocks or elements identified, but that such blocks or elements do not comprise an exclusive list and a method or apparatus may contain additional blocks or elements.

[0075] The steps of the methods described herein may be carried out in any suitable order, or simultaneously where appropriate. Additionally, individual blocks may be deleted from any of the methods without departing from the spirit and scope of the subject matter described herein. Aspects of any of the examples described above may be combined with aspects of any of the other examples described to form further examples without losing the effect sought.

[0076] It will be understood that the above description of a preferred embodiment is given by way of example only and that various modifications may be made by those skilled in the art. Although various embodiments have been described above with a certain degree of particularity, or with reference to one or more individual embodiments, those skilled in the art could make numerous alterations to the disclosed embodiments without departing from the spirit or scope of this invention.

Claims (19)

1. Claims 1. A method of interrogating an optical sensing system performed at an interrogation system, comprising the steps of receiving a plurality of bursts of samples of an optical signal, each burst of samples comprising a plurality of samples of each pulse in a group of discrete optical pulses, wherein the groups are regularly spaced in time with a group repetition period larger than the group width, and the optical pulses comprise a phase modulated signal; processing the bursts of samples to obtain a series of regularly spaced samples; and demodulating the regularly spaced samples to retrieve the phase information from the detected pulses, wherein the step of processing comprises the steps of downsampling the received samples, upsampling the downsampled samples at a rate of the group width multiplied by the sample rate in each burst, filter the upsampled samples, and combining the filtered samples.
2. 2. A method according claim 1, wherein processing the bursts of samples comprises processing the samples utilising a filter bank.
3. 3. A method according to claim 2, wherein the filter bank comprises a number of parallel paths equal to the number of samples in each burst.
4. 4. A method according to any preceding claim, wherein the downsampling ratio is equal to the sample rate in each burst divided by the group repetition frequency.
5. 5. A method of interrogating an optical sensing system performed at an interrogation system, comprising the steps of receiving a plurality of bursts of samples of an optical signal, each burst of samples comprising a plurality of samples of each pulse in a group of discrete optical pulses, wherein the groups are regularly spaced in time with a group repetition period larger than the group width, and the optical pulses comprise a phase modulated signal; calculating the instantaneous frequency of the signal represented by each burst of samples; and recovering phase information of the detected pulses.
6. 6. A method according to claim 5, wherein the step of calculating comprises filtering the samples utilising an FIR filter having a length less than the number of samples in each burst of samples.
7. 7. A method of interrogating an optical sensing system performed at an interrogation system, comprising the steps of receiving a plurality of bursts of samples of an optical signal, each burst of samples comprising a plurality of samples of each pulse in a group of discrete optical pulses, wherein the groups are regularly spaced in time with a group repetition period larger than the group width, and the optical pulses comprise a phase modulated signal; phase locking a signal to the samples of each burst of samples; and utilise the phase locked signal to recover phase information of the detected pulses.
8. 8. A method according to any preceding claim wherein each group of pulses consists of a single pulse.
9. 9. A method according to any preceding claim, wherein each burst of samples comprises 8 to 80 samples.
10. 10. A method according to any preceding claim wherein the samples within each burst are at a sampling frequency of 40 to 400 MHz.
11. 11. A method according to any preceding claim, wherein the group width is 200ns.
12. 12. A method according to any preceding claim, wherein the group repetition period is in the range of 0.5 Rs to 10 is.
13. 13. A method according to any preceding claim, wherein the regularly spaced samples have a sample rate in the range of 10MHz to 20MHz.
14. 14. A method according to any preceding claim, further comprising the steps of detecting the optical pulses, and sampling the pulses to obtain the bursts of samples, and transmitting the bursts of samples to the interrogation system.
15. 15. A method according to claim 14, wherein the steps of detecting, sampling, and transmitting are performed at a different location to the other steps of the method.
16. 16. A method according to claim 14, wherein the optical pulses are returning from an optical sensor.
17. 17. A method according to claim 14 to 16, further comprising the steps of detecting other optical pulses returning from other optical sensors, wherein the optical pulses returning from different sensors are time multiplexed.
18. 18. A method according to any preceding claim, wherein the phase change between adjacent samples of the regularly spaced samples is less than 7C.
19. 19. An interrogation system configured to perform the methods of any of claims 1 to 18.