Device and Method for Determining the Composition of a Mixture of Fluids

Abstract

A device for determining the composition of a mixture of fluids that flow along a pipe includes: a radiation source for illuminating the mixture with radiation; a detector for detecting radiation that has been attenuated by the mixture; and a device for monitoring the flow rate of fluid along the pipe and outputting a signal indicative of the flow rate. The device for determining further includes a device for adjusting the intensity of radiation emitting by the radiation source in response to the signal indicative of the flow rate so that the intensity of the radiation source is reduced if the flow rate reduces.

Claims

1-15. (canceled)

16. A device for determining the composition of a mixture of fluids that flow along a pipe, which comprises: a radiation source for illuminating the mixture with radiation; a detector for detecting radiation that has been attenuated by the mixture; a device for monitoring the flow rate of fluid along the pipe and outputting a signal indicative of the flow rate; wherein it includes a device for adjusting the intensity of radiation emitted by the radiation source in response to the signal indicative of the flow rate so that the intensity of the radiation source is reduced if the flow rate reduces.

17. A device according to claim 16 wherein said radiation source is a broadband NIR source synthesised from solid state devices.

18. A device according to claim 16, wherein the radiation source is a broadband radiation source that is operable to emit radiation over the entire range of absorption bands of the fluids.

19. A device according to claim 16, wherein it includes a digital mirror device for modulating radiation from the radiation source.

20. A device according to claim 16, wherein said radiation source has a range of wavelengths that extends over a predefined absorption band of each one of said fluids and over at least one predefined reference band outside, but adjacent to, each and every one of said absorption bands, and wherein it comprises a separating device configured for separating the radiation that has been attenuated by the mixture before said radiation reaches said detector, and wherein said separating device comprises: a first radiation divider configured for separating said radiation into different wavelength bands, each one of which extending over the absorption band of a corresponding one of said fluids and over at least one of said reference bands adjacent to said absorption band; a second radiation divider configured for receiving from said first radiation divider the radiation corresponding to each one of said different wavelength bands and for further separating it into a wavelength band corresponding to the respective said absorption band and a wavelength band corresponding to the respective at least one of said reference bands.

21. A device as claimed in claim 20, wherein the separating device comprises a dichroic beam splitter and/or a rugate notch filter.

22. A device as claimed in claim 21, wherein said first radiation divider includes a plurality of dichroic beam splitters to split the radiation into a plurality of wavelength bands, and said second radiation divider includes a rugate notch filter configured for splitting the radiation in each wavelength band into an absorption band and a reference band having a wavelength range adjacent to the absorption band.

23. A device as claimed in claim 22, wherein it is arranged so that radiation impinges on the dichroic beam splitters and/or rugate notch filters at an angle of incidence of not more than 20°.

24. A device as claimed in claim 20, wherein it is operable to multiplex radiation from a plurality of narrowband radiation sources to generate the radiation, and the detector is operable to demultiplex the detected radiation.

25. A method of determining the composition of a mixture of fluids that flow along a pipe, comprising the steps of: illuminating the mixture with radiation from a radiation source; detecting radiation that has been attenuated by the mixture; monitoring the flow rate of fluid along the pipe; wherein it further comprises the step of adjusting the intensity of radiation emitted by the radiation source in response to the flow rate determined so that the intensity of the radiation source is reduced if the flow rate reduces.

26. A method as claimed in claim 25, wherein said radiation is emitted from a radiation source having a range of wavelengths that extends over an absorption band of each one of said fluids and over at least one reference band outside, but adjacent to, each one of said absorption bands, and wherein it further comprises the steps of: separating the radiation that has been attenuated by the mixture into different wavelength bands, each one of which extending over the absorption band of a corresponding one of said fluids and over at least one of said reference bands adjacent to said absorption band; further separating the radiation corresponding to each one of said different wavelength bands into a wavelength band corresponding to the respective said absorption band and a wavelength band corresponding to the respective at least one of said reference bands; said separation steps being performed before said step of detecting the radiation.

27. A method as claimed in claim 26, wherein said step of separating the radiation comprises separating the radiation into two of said reference bands located adjacent and on opposite sides of the absorption band of one of the fluids.

28. A method as claimed in claim 27, wherein said two reference bands are located substantially symmetrically on either side of said absorption band.

29. A method as claimed in claim 27, wherein it includes the step of interpolating the level of absorption of the radiation in said two reference bands in order to generate a baseline for the corresponding absorption band.

30. A method as claimed in claim 25, wherein the mixture comprises water and hydrocarbons.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Various forms of device and method according to the invention will now be described by way of example, with reference to the accompanying drawings in which:

[0060] FIG. 1 is a schematic view of a measurement probe for determining the composition of a composition flowing in a pipe;

[0061] FIG. 2 shows a typical range of Near Infra Red (NIR) absorption spectra for oil, water methane;

[0062] FIG. 3 is a schematic view showing a beam splitter employed in the device according to the invention;

[0063] FIG. 4a is a schematic view showing a generic beam splitter employed in the device according to the invention;

[0064] FIG. 4b shows the characteristic of a broadband radiation source passed to a first dichroic beam splitter DBn that transmits radiation of wavelength grater than a specified value λn;

[0065] FIG. 4c shows the characteristic of a broadband radiation source passed to a dichroic beam splitter DBn+1 that transmits radiation of wavelength greater than a second wavelength λn+1;

[0066] FIG. 4d shows the characteristic of the radiation reflected by the second beam splitter DBn+1;

[0067] FIG. 5a shows the bands A, B and C reflected by a Rugate beam splitter employed in the device according to the invention;

[0068] FIG. 5b shows filter characteristics of a Rugate beam splitter employed in the device according to the invention;

[0069] FIG. 6 is a schematic view showing the principle of an alternative form of beam splitter according to the invention;

[0070] FIGS. 7a and 7b are schematic figures showing the operation of a digital mirror device employed in an embodiment of the invention;

[0071] FIG. 8 is a schematic view of an alternative embodiment in which the angle of the light beam at the filters is reduced.

[0072] FIG. 9 is a graphical representation of the transmissivities of the system at the detectors of the embodiment of FIG. 8;

[0073] FIG. 10a shows the theoretical values of the errors in the water cut determination for low and high background scatter theoretically obtained by the embodiment according to the invention;

[0074] FIG. 10b shows the values of the errors in the water cut determination for low and high background scatter obtained by the manufactured embodiment according to the invention;

[0075] FIG. 10c shows the values for low and high background scatter obtained using a system that employs a remote reference;

[0076] FIG. 11 is a graphical representation of the absorption of the mixtures used in generating the errors shown in FIG. 10b;

[0077] FIG. 12 is a schematic view of an alternative form of measurement probe in which a number of narrowband radiation sources are time multiplexed; and

[0078] FIG. 13 is a schematic of the signal processing circuit for determining the absorption of radiation by a fluid using a device shown in FIGS. 1 to 9.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0079] FIG. 1 is a schematic representation of a spectroscopic optical fluid analyser (SOFA) probe in which a spectrally broadband input source light distribution, Pi(λ), is generated by transceiver T and sent along an optical fibre link F1 to a Single Side Transmission (SST) retro-reflective measurement probe RP that is immersed in the specimen S. The SST probe may be replaced by a double sided transmission probe depending on the specific installation. In the measurement probe, the specimen is constrained to a depth d in the flow channel of the probe. After absorption of the radiation by the specimen, the radiation is reflected by the head of the probe and returned to the detector of the transceiver T by a second optical fibre F2. The output power Po(λ) is spectrally attenuated as the result of broadband and molecular resonant absorption, and the output current of i.sub.o(λ) of the detected signal is given by,

l.sub.o=P.sub.i(λ)Γ(λ)[1−σ.sub.m(λ)]  (7)

where, [0080] Γ(λ) is the detector responsivity (Amp/Watt) as function of wavelength λ, [0081] t.sub.s(λ) is the Broadband spectral transmission of the specimen dependent on non-resonant effects (e.g., scatter); and [0082] σ.sub.m(λ) is the resonant molecular adsorption of the specimen.

[0083] A conventional spectrometer may be used for general detection and measurement of the absorption spectra as, for example, in a utility test system. However, a general spectrometer is preferably replaced, for reasons discussed above, by a component specific DOAS/rugate detection unit shown in FIGS. 3 and 4a described below.

[0084] The probe shown in FIG. 1 may employ a broadband radiation source derived, for example, from a conventional spectrally broadband incandescent source such as a quartz halogen lamp, or alternatively the broadband source may be synthesized by combining radiation from different sources.

[0085] For example, radiation from a plurality of narrowband sources such as LEDs, super luminescent light emitting diodes (SLEDs), or thermorestive sources may be superimposed using a dichroic beam combiner, or any other suitable means for combining beams and coupled into the illumination optical fibre. Although such sources will couple less power into the fibre than a conventional lamp, and therefore may require the system to measure more slowly, they do have intrinsically longer lifetimes than which makes them suitable for applications such as sub-sea or subterranean measurement where the light source cannot be changed routinely. They may also be employed in other applications where the light source cannot easily be changed, for instance where the device is enclosed in a sealed unit that may be provided to prevent the risk of fire or explosion due to the presence of the hydrocarbon gases. They are also particularly suitable for life time extension using the modulation of the output power in response to variations in flow rate as discussed earlier. Yet another form of radiation source is a rugged incandescent ceramic radiation source, for example using a silicon nitride heating element. Such sources are compact, mechanically robust and low cost and emit radiation at wavelengths from 1 to 2 μm. They may typically be employed as gas igniters in domestic cookers and hobs.

[0086] FIG. 2 is an absorption spectrum that is typical for a mixture of hydrocarbons, optionally including gaseous hydrocarbons, and water that may be obtained in the bore of an oil well. The absorption spectrum may be made up of spectra from a number of components of the mixture. For example, spectrum 2 that is due to water exhibits an absorption peak 4 at approximately 1450 nm and extending from about 1400 to 1530 nm, and a second, larger, absorption peak at approximately 1930 nm. Three absorption spectra are shown for oil, curve 10 is the spectrum for filtered and shaken oil which may contain a number of small bubbles, for example containing air or gaseous hydrocarbons, that contribute to scattering of the radiation. Since the various different mixtures will be recorded by the spectrometer in rapid succession as the mixture flows along the pipe, the spectrum therefore will have a relatively large baseline caused by noise. Curve 12 is the spectrum for filtered and settled oil which has a significantly lower baseline due to the lack of air bubbles, while curve 14 is the spectrum of unfiltered and settled oil (type 3) with a relatively low baseline. Each of the spectra exhibits a hydrocarbon absorption peak 16 at approximately 1730 nm, extending from about 1700 to 1750 nm. In addition, it can be seen that the baseline of the curves is not flat, but is significantly larger at shorter wavelengths. Finally, curve 18 is the spectrum for methane or other low molecular weight gaseous hydrocarbons.

[0087] This curve exhibits an absorption peak at approximately 1670 nm extending from about 1650 to 1690 nm.

[0088] FIG. 3 is a schematic view showing the principle of the present invention in which typical filter and Rugate wavelength bands (as indicated as a guideline below) enable referenced measurements of the NIR absorption bands of water, methane and oil defined above. This device comprises a series of dichroic beam splitters DB1 to DB4 that split a beam 32 of radiation from a spectrally broadband source (not shown) into a number of beams 34 to 40 of different wavelengths. In operation, the light from this source is first delivered to the measurement fluid by the optical fibre F1 and then delivered at 32 in FIG. 3 by the fibre F2 as shown in FIG. 1 after being spectrally modulated as a result of spectral absorption by the by the measurement fluid.

[0089] The principle of beam splitting shown in FIG. 4a is shown in FIG. 4b where a broadband radiation source is passed to a first dichroic beam splitter DBn that transmits radiation of wavelength greater that a specified value λn, in this case 1345 nm, and reflects radiation of a wavelength below that value into beam 33.

[0090] The transmitted beam is then sent to dichroic beam splitter DBn+1 that transmits radiation of wavelength greater than a second wavelength λn+1, in this case greater than 1590 nm, and reflects radiation of wavelength less than λn+1 as shown in FIG. 4c so that the radiation reflected by the second beam splitter will have a wavelength from 1345 to 1590 nm as shown in FIG. 4d. In this way, the original beam is split into four beams having wavelengths less than 1345 nm, 1345 to 1590 nm, 1590 to 1690 nm and 1690 to 1825 nm respectively. The first beam 34 of wavelength less than 1345 nm is an overall reference channel which may be bandwidth limited by an additional notch filter Fr, and provides an auxiliary, non-resonant adsorption reference at nominally 1300 nm as required.

[0091] The remaining beams 36, 38 and 40 each have wavelengths corresponding to an absorption band of interest and also wavelengths that are outside, but adjacent to, the absorption bands. These beams are passed to a Rugate beam splitter 42 having a characteristic shown schematically in FIGS. 5 (a) and (b). In principle, any optical notch filter having the appropriate characteristic may be employed, for example a dielectric quarter-wave stack, but Rugate notch filters have the advantage that they exhibit significantly lower amplitude higher-order harmonic structure in the Rugate notch. The Rugate notch filter 42 will reflect that part of beam 36 having wavelengths of 1400 to 1530 nm, shown as region A in FIG. 5a, corresponding to the water absorption band, into detector D2a and transmit the remaining parts of the beam, i.e., of wavelength 1345 to 1400 nm and 1530 to 1590 nm into detector D2r to provide a reference.

[0092] Similarly, the Rugate notch filter will reflect that part of beam 38 having wavelengths of 1650 to 1690 nm, shown as region B in FIG. 5a corresponding to the gaseous hydrocarbon (methane) absorption band, into detector D3a while transmitting the remaining parts of the beam into the reference detector D3r. Finally, the Rugate notch filter will reflect that part of beam 40 having wavelengths of 1700 to 1750 nm, shown as region C in FIG. 5a, corresponding to the oil absorption band, into detector D4a, while transmitting 15 the remaining parts of the beam, i.e., of wavelengths 1690 to 1700 nm and 1750 to 1800 nm into reference detector Dar. In this way, the device according to the invention allows the intensity of the radiation attenuated by the water, methane and oil absorption bands respectively to be detected, and to detect the intensity of non-resonant reference signals on each side of the absorption bands.

[0093] The above wavelengths give a general indication of those that would be used specifically for DOAS water, methane and oil measurement. They may in practice be modified for optimum operation and in particular may be modified to enable the measurement of water using the absorption band centred at nominally 1950 nm. The wavelengths will be changed entirely when the same general principle is used for the measurement of different molecular constituents.

[0094] The values for attenuation of the radiation detected by the reference detectors on either side of the absorption bands may simply be averaged in order to 30 provide a baseline for the absorption bands. In view of the wavelength dependency of the background absorption of the radiation shown in FIG. 2, the baseline for the absorption bands may not be the same for each band.

[0095] Further, where the absorption band extends over a significant wavelength range, the baseline for the absorption may not be horizontal but may be 35 formed as a straight-line interpolation of the intensity measured by the reference detector on either side of the absorption band.

[0096] It is not necessary to employ a number of dichroic beam splitters to divide an original beam into a number of separate beams of smaller wavelength range, and to pass each divided beam through a Rugate filter as described above.

[0097] Other arrangements may be employed instead. For example, one or more Rugate filters may be used to divide the original beam into beams of different wavelengths and the separate beams may be passed to one or more dichroic beam splitters or to a further Rugate filter that is different from the first Rugate filter in order to form separate absorption and reference beams. Similarly, it is possible to split the original beam into the appropriate wavelength ranges to using only dichroic beam splitters. On the other hand, it is possible to use dielectric quarter-wave stacks instead of rugate notch filters as mentioned above.

[0098] It is not essential to the invention that the reference detectors detect the radiation intensity on both sides of each absorption peak. It is possible, for 15 example for the reference detectors to detect the radiation on each side of a pair of absorption bands or on each side of all three absorption bands and to average or interpolate the baseline if necessary using the reference signals.

[0099] FIG. 6 shows an alternative arrangement in which reference absorption is detected on each side of a pair of absorption bands. The purpose of this arrangement is to simplify the design of the Rugate by reducing the number of spectral notches. In this arrangement, the original broadband radiation beam 32 is split by three dichroic beam splitters DB1 to DB3 into beams 60, 62 and 64, the first of which has a wavelength less than Ai (1345 nm) and is filtered by a notch filter Fr as described above to provide an auxiliary, non-resonant 25 absorption reference at nominally 1300 nm as required. The transmitted beam is split by beam splitter DB2 to reflect radiation of wavelength less than λ2 (1590 nm) and is passed to the rugate beam splitter 62 in the manner described above in order to split the radiation so that radiation having a wavelength of 1400 to 1530 nm corresponding to the water absorption band can be detected by detector D2a and reference bands having wavelengths of 1345 to 1400 nm and 1530 to 1590 nm are detected by detector D2r as described above. The radiation that is transmitted by dichroic beam splitter DB2 having a wavelength greater than λ2 (1590 nm) is passed to dichroic beam splitter DB4 which reflects radiation of wavelength below λ4 (1825 nm) corresponding to both the oil and gaseous hydrocarbon component absorption bands onto the rugate beam splitter 62. Radiation of wavelength below 1650 nm and above 1750 nm is transmitted by the Rugate beam splitter 62, into detector D3r to detect the reference background radiation intensity on each side of the pair of absorption bands for oil and gaseous hydrocarbons.

[0100] Beam 64 that is reflected by the rugate beam splitter 62 and has wavelengths of from 1650 to 1750 nm is passed to dichroic beam splitter DB3 that reflects 5 radiation of wavelength below λ3 (1690 nm), corresponding to the methane absorption band, into detector D3a1 and transmits radiation of wavelength above λ3 corresponding to the oil absorption band, into detector D3a2.

[0101] This simplification of the Rugate component is expected to reduce the cost of the overall system.

[0102] It is noted that an identical source and detection system may be used in combination with other designs of probe such as double pass transmission and ATR (Attenuated Total-internal Reflection).

[0103] One limitation with employing broadband light sources such as incandescent sources is that it is difficult to modulate them electronically at the frequencies required for effective compensation of the dark current noise of the near infrared detector (for example an InGaAs detector). To overcome this [imitation an external device may be used to modulate the light field. For example, the invention as described may, for this purpose, use a compact spatial light modulator (SLM) in which a digital mirror device (DMD) provides an intrinsically zo high speed and high contrast on/off ratio. This is shown in principle in FIG. 7a and FIG. 7b.

[0104] In FIG. 7a, lens L1 forms an image of the output aperture of the source S in the plane of the aperture L2 via reflection of the beam by the DMD, and L2 images the output aperture of L1 in the plane of the output optical fibre F that couples light to the sensor probe immersed in the fluid. The apertures and focal lengths of L1 and L2 are chosen to match optimally the phase volume of the source to that of the fibre.

[0105] The DMD consists of an array of multiple square mirrors (in this case approximately 7 μm×7 μm) each of which can be individually actuated as 30 shown in FIG. 7b to switch between angular orientations +θ (position a) and −θ (position b) via the neutral position b (θ=θ). In operation the light field coupled into the output fibre is modulated by simultaneously switching all of the mirrors in the array between the +θ and −θ0 position corresponding to the on state in which the radiation is coupled into the output fibre F and the off state in which the light is reflected into a beam dump D. The modulation rate may be as high as 6.5 kHz or more by virtue of the DMD. The DMD may be employed with any form of incandescent light source, but it is particularly convenient to use it with the ceramic radiation source referred to above.

[0106] The design of device shown in FIGS. 3 and 6 in which the radiation is reflected by the rugate filters at a 45° degree angle of incidence can have the disadvantage that the s and p polarisation components of randomly polarised light delivered from the sensing probe may have different filter transmission responses that result in spectral blurring and ringing of the light transmitted or reflected by a given filter. This could be eliminated by the use of a polarising filter, but this would entail a minimum 50% light loss and efficacy over a limited to wavelength range. These effects can, however, be eliminated or substantially eliminated over the full operating wavelength range without significant light loss by limiting the angle of incidence to a maximum value of not more than 20° and preferably not more than 15° (to the normal to the plane of the filter).

[0107] In one preferred design of device an image of the light field 32 at the aperture of the input fibre connected to the probe output is formed at the detector apertures S1 to S5 by selected lenses in the lens chain L1 to L8 and hence via transmission and reflection at selected filters and Rugates within this chain.

[0108] In this arrangement the light first passes through NRF1 and then NRF2 via the broadband, high efficiency front reflecting mirrors IM to be incident on the long pass filter LP1. The purpose of the filters NRF1 and NRF2 is to reject low and high out of band spectral noise at the sensor input. The long wavelength component of the light transmitted by LP1 passes to the Rugate filter R1 where it is divided into the transmitted and reflected components. The former forms the reference bands for the oil/methane measurement after transmission through the trimming short pass filter TSP1 which suppresses the longer wavelengths present. The latter is reflected by the double sided mirror DM via lens L4. Light reflected by the mirror DM is again reflected by the Rugate R1 and mirror DM before being passed to low pass filter LP2 which transmits light of wavelength of the oil absorption band to the trimming band pass filter TPB3 which spectrally limits light of the oil absorption band, and finally to sensor S3 for detecting the oil concentration. Light reflected by the low pass filter LP2 has a wavelength corresponding to the methane absorption band and is passed to the trimming band pass filter TBP2 and hence to sensor S2 for determining the amplitude of the methane peak.

[0109] Light reflected by the long pass filter LP1 is incident on Rugate R2 which transmits light of the water reference wavelength to the water reference frequency detector S4 via trimming short pass filter TSP4 which suppresses long wavelength of the water reference band. Light of wavelength of the water absorption band is reflected by the Rugate R2 back to the double sided front reflection mirror DM whereupon it is reflected back to the trimming band pass filter TBP5 for the water measurement and detector SS for detecting the water absorption peak.

[0110] In this detector, the sensor may be configured to measure either methane and oil combined at detector S2/3 by an optical flat, broadband anti-reflection element coated on both faces which has a high transmission for both methane to and oil absorption bands, or to measure oil at detector S2 and methane at detector S2 by the use of the low pass filter LP2.

[0111] In this design of detector, the Rugate R1 and the double sided front mirror DM may be separated by a significant distance, for example in the range of 30 to 50 mm so that the light beams will be spaced apart laterally by a sufficient distance to enable detectors S1 and S2 or S3 to be used, and to enable different detectors S4 and SS to be used while keeping the angle of incidence of the light rays at the various filters and mirrors at a low value of about 15°.

[0112] FIG. 9 shows the concatenated spectral transmissivity of each of the channels S1 to S5 of the arrangement shown in FIG. 8 obtained by multiplying the spectral reflectances/transmissivities of each component within a given channel. As can be seen, the transmissivity peak for the detector S5 (1400 to 1510 nm) corresponding to water absorption is bounded on each side by transmissivity peaks at 1300 to 1380 nm and at 1510 to 1590 nm for the detector S4 which peaks are used as references for the water detection peak observed by detector S5. Similarly, detector S2 has a transmissivity peak at 1650 to 1690 nm, which detects the presence of methane, and detector S3 has a transmissivity peak at 1700 to 1750 nm which detects the presence of oil. This pair of transmissivities is bounded by a pair of transmissivity peaks occurring at 1590 to 1650 nm and at about 1770 to 1850 nm which are used as references for determining the value of the methane and oil absorption peaks.

[0113] As can be seen, the arrangement according to the invention enables the reference bands employed to determine the intensity of the water and the hydrocarbon absorption bands to abut the bands, and the peak value for the transmissivity values for the references will normally be not more than 150 nm and especially not more than 100 nm from the peak value for the absorption bands. A general aim is to make the reference beams symmetrical about the absorption bands in order to minimise errors due to the spectral gradient of the background signal.

[0114] The improvement in accuracy obtained by local referencing, i.e., by determining the absorption due to the presence of water and hydrocarbons with reference to the absorption at a wavelength that is adjacent to the absorption peaks is shown in FIGS. 10a, 10b and 11. FIG. 11 is an absorption spectrograph for six mixtures of hydrocarbon and water at three different values for the water cut (WC) and at two different values for the background scatter due to turbulence, as follows:

[0115] Mix 1: W.sub.c=50%; High background

[0116] Mix 2: W.sub.c=50%; Low background

[0117] Mix 3: W.sub.c=99%; High background

[0118] Mix 4: W.sub.c=99%; Low background

[0119] Mix 5: W.sub.c=1%; High background

[0120] Mix 6: W.sub.c=1%; Low background.

[0121] These fluids were used for determining the filter/reference performance simulations shown in FIGS. 10a, 10b and 10c. The graph shown in FIG. 11 is similar to that shown in FIG. 2 but somewhat more accurate. The curve for mix 6 has the lowest value because the hydrocarbons have lower absorption than that of water, and there is relatively little scatter due to turbulence. The curve for mix 5 shows a strong wavelength dependency due to Rayleigh and Mie scatter. The remaining curves show water absorption peaks depending on the water cut (Wc) superposed on a baseline that decreases with wavelength due to scatter.

[0122] The results for the measured water cut values given against reference water cut values at low and high values of background scatter are shown graphically in FIGS. 10a, 10b and 10c. FIGS. 10a and 10b show predicted noise-free values of the water cut against reference values for both low and high background scatter (caused by small particles or bubbles in the mixture) employing an algorithm based on equation 1 incorporating calibration terms for 100% water and 100% oil in the probe.

[0123] FIG. 10c is a similar graph showing the values for low and high background scatter obtained using a system that employs a remote reference as in U.S. Pat. No. 6,292,756. The dotted straight line in each of the graphs extending linearly from 0% to 100% measured and reference water cut shows the actual value of the water cut. The upper two graphs, demonstrating a system according to the present invention with local referencing, have an error in the water cut value due to the presence of scatter background which is less than 4% even when obtained with high background scatter, whereas the equivalent error in the water cut value using a system with remote referencing is in excess of 40%. For results obtained using non-local referencing, even in the case of low background scatter the output departs from the “truth” line (where the measured water cut is equal to the reference water cut). The close correlation between the outputs shown in FIGS. 10a and 10b obtained respectively for the theoretical and manufactured spectral transmissivity of the filters (c.f. FIG. 9) demonstrates the practical feasibility of the overall system.

[0124] While it is preferred for the measurement probe to employ a broadband radiation source, either formed from a single device or from a number of narrowband devices, as shown in FIG. 1, it is possible to employ a number of time-multiplexed narrowband or quasi narrowband sources as shown in FIG. 12. Radiation from a number of discrete sources centred at wavelengths λ1 to λn with corresponding bandwidth Δλn, is time and wavelength multiplexed to propagate sequentially at times t1 to tn through the system. The power transmitted by the specimen at each wavelength is measured sequentially at times t1 to tn using time and wavelength de-multiplexed detection.

[0125] The detection electronics are synchronised with the time multiplexed input signal. This enables the source and hence wavelength at which the absorption is measured to be identified. The spectral distribution of the input sources (i.e., either a Light Emitting Diode: LED, Laser Diode: LD, or Super Luminescent Light Emitting Diode: SLED), may be selected to match the spectral ranges employed in the measurement. Under these conditions a Rugate beam splitter system may be used in combination with time de-multiplexing of the signals detected by the photodetectors for absorption in the relevant bands of the liquid, the reference bands and subsidiary reference detection (˜1300 nm). The ratio of the absorption to reference signals recorded in corresponding time slots generate the DOAS signal required for the measurement of the concentration of water, methane and oil in the presence of spectrally broad band variations in absorption.

[0126] FIG. 13 is a schematic showing an electronic circuit that may be employed to process the radiation intensity signals that are generated by the detectors Dna and Dnr of the device. One problem with the output of the detectors of the device is that the background signal on which the absorption band signals are superposed can have a very large dynamic range that may be significantly greater than the differences between the absorption bands, and which varies rapidly with time. This has the result that, if a conventional circuit were used, an extremely high resolution analogue-to-digital converter would be required to cover the full dynamic range of the signals while being sensitive enough to measure small differences between the signals. The circuit that is employed 5 according to this aspect of the invention has the advantage that lower cost components may be used while having the required bandwidth do deal with changes in the signal intensities.

[0127] The circuit comprises a pair of a pair of photodetectors 71 and 72 forming absorption detectors Dna and reference detectors Dnr respectively for detecting to the absorption band attenuated radiation intensity and the reference radiation intensity for one absorption band. Additional circuits corresponding to the circuit shown in FIG. 4a are required to detect the intensity of the other absorption bands. The detectors may be for example a pin photodiode or an avalanche photodiode as is conventional. The outputs from the photodetectors 70 and 72 are fed into preamplifiers 74 and 76 whose output is fed into time division mulitiplexer 78. The outputs I1 and I2 of the preamplifiers 74 and 76 respectively are chopped alternately with the output of a high frequency oscillator 80 to give a square wave alternately proportional to I2 and I1, and the output is fed into the input of a variable gain differential amplifier of a closed-loop gain block 84. Time-division multiplexing of the outputs I1 and I2 is performed so that the same gain can be applied both outputs by the differential amplifier 82. The output from gain block 84 is fed into a second gain block 86.

[0128] As shown, three gain blocks are employed, but more or fewer may be used as is necessary in order that the circuit is capable of handling the dynamic range of the variations in intensity of the detected radiation. The purpose of the gain blocks is to reduce the dynamic range of the signals to a manageable level while preserving the ratio information of I1 and I2.

[0129] In the initial gain stage 84, where the input is the chopped square wave of I1 and I2, an amplified version of the input formed by the amplifier 82 is demodulated with the oscillator signal by demodulator 90 to give just I2.

[0130] The signal I2 from the demultiplexer is filtered with respect to a reference signal 92 by means of an integrating amplifier 94. The output of the integrating amplifier, which will be proportional to I2, is used to set the gain of the amplifier 82, for example by means of a servo, so that:

G.Math.I.sub.1=ref

where G is the gain of amplifier 82 and ref is the reference voltage.

[0131] The gain G is thus inversely proportional to the value of I.sub.2, or

[00006]$\begin{array}{cc}G=\frac{\mathrm{ref}}{I_2}& \left(8\right)\end{array}$

[0132] In the case of an amplifier circuit formed from a number of gain stages 84, 86 and 88, in which the gain of each of the gain stages is set with reference to the reference voltage ref, the total gain, G.sub.r is given by equation 2.

[0133] The signal is then demodulated by demodulator 96 using the anti-phase of the oscillator 80 formed by inverter 98 to remove I.sub.2 and leave I.sub.1 as the input to the to the amplifier. The signal output by the overall amplifier is thus G.sub.T.Math.I.sub.1 which is given by

[00007]$\begin{array}{cc}{G}_T{I}_1=\frac{I_{1,\mathrm{ref}}}{I_2}& \left(9\right)\end{array}$

[0134] Since ref is known, and the same gain has been applied to both I.sub.1 and I.sub.2, the is required ratio I.sub.1/I.sub.2 can be calculated in the digits from the output signal G.sub.TI.sub.1.

[0135] The output of the circuit is thus not dependent on the rapidly time-varying background signal, but only on the relatively slowly varying ratio I.sub.1/I.sub.2 which has a much smaller dynamic range than the absolute value of the individual absorption bands.