Patent Application
20240011824
The present invention relates to optical sensors and, in particular, distributed optical fibre sensors and applications thereof.
The benefits of optical fibres have been demonstrated in a number of sensing applications. The two major areas are: (i) distributed optical fibre sensors, and (ii) multiplexed point sensor arrays.
Distributed sensors utilise the intensity of backscatter light, with Raman and/or Brillouin peaks in the light signal utilised to measure temperature, strain or pressure. Distributed sensors offer a number of advantages including continuous sensing along the entire length of fibre, and flexibility and simplicity of the sensor, which may be standard telecoms optical fibre. For example, a distributed sensor may provide 10,000 measurement points along 10 km of optical fibre with a 1 m spatial resolution. Distributed sensor systems therefore offer low installation and ownership costs.
However, due to their slow response, distributed sensors are usually only used in applications where measurements taking in order of several seconds to hours are acceptable. The most common sensors of this type are the distributed temperature sensors (DTS), which are made by a number of companies. A typical performance of a DTS is 1 m spatial resolution and 1° C. temperature resolution in 60 seconds over a 10 km range.
Distributed sensors have also been used to measure strain by utilising Brillouin shifts in reflected or backscattered light, as described in U.S. Pat. No. 6,555,807 [1] or WO 98/27406 [2]. The frequency of the Brillouin shift is about 1 MHz/10 με and its linewidth is about 30 MHz. The strain in an order of 10 με can be determined along an optical fibre using the narrow frequency scanning methods described. However, using these approaches, the scanning rate is much slower than the pulse repetition rate and measurement times are typically in the order of few seconds to few minutes.
More recently, a technique for faster measurement of Brillouin frequency shift has been proposed in U.S. Pat. No. 7,355,163 [3]. This technique uses a frequency to amplitude convertor which may be in a form of an optical fibre Mach-Zehnder interferometer with a 3×3 coupler at its output. However, the strain resolution is limited by the linewidth of the Brillouin light and therefore the optical path length difference in the interferometer should be kept within the coherence length of the Brillouin light. Also, the polarisation fading between the two paths of the interferometer, the offset and gain variations of the photodetector receivers would significantly limit the strain measurement. Measurement times of around 0.1 seconds (10 Hz) with strain resolution of 50 με have been recently reported using this technique.
For many applications, such as acoustic sensing, much higher sensitivities and faster a measurement time in the order of 1 millisecond (1 kHz), 0.1 millisecond (10kHz) or 0.01 millisecond (100 kHz) is required.
Multiplexed point sensors offer fast measurements with high sensitivity and are used, for example, in hydrophone arrays. The main application for these in the energy market is for towed and seafloor seismic arrays. However, unlike with distributed sensors, multiplexed point sensors cannot be used where full coverage is required. The size and the position of the sensing elements are fixed and the number of sensors multiplexed on a single fibre is typically limited to 50 to 100 elements. Furthermore, the sensor design relies on additional optical fibre components leading to bulky and expensive array architectures. There is also considerable effort to increase the number of sensors that can be efficiently multiplexed on a single length of fibre.
Optical-time-domain reflectometry (OTDR) is a well known technique that has been used to test optical fibre communications cables. In order to reduce the effect of coherent backscatter interference, which is sometime is referred to as Coherent Rayleigh Noise, a broadband light source is normally used. However, proposals have also been made in U.S. Pat. No. 5,194,847 [4] to use coherent OTDR for sensing intrusion by detecting the fast changes in a coherent backscatter Rayleigh signal. In addition, Shatalin et al. [5] describes using coherent Rayleigh as a distributed optical fibre alarm sensor.
WO 2008/056143 [6] describes a disturbance sensor similar to that of U.S. Pat. No. 5,194,847 [4] using a semiconductor distributed feedback laser source. A fibre Bragg grating filter of preferably 7.5 GHz is used to reject out-of-band chirped light and, thereby, improve the coherence of the laser pulse sent into the fibre. However, this requires matching of the laser wavelength with the narrow band optical filter, which results in the signal visibility variation being reduced compared to a system which uses a very high coherent source as proposed by U.S. Pat. No. 5,194,847.
Similar techniques have also been proposed for the detection of buried optical fibre telecommunication cables (for example in WO 2004/102840 [7]), in perimeter security (GB 2445364 [8] and US2009/0114386 [9]) and downhole vibration monitoring (WO 2009/056855 [10]). However, the response of these coherent Rayleigh backscatter systems has been limited by a number of parameters such as polarisation and signal fading phenomena; the random variation of the backscatter light; and non-linear coherent Rayleigh response. Therefore these techniques are mainly used for event detection and do not provide quantitative measurements, such as the measurement of acoustic amplitude, frequency and phase over a wide range of frequency and dynamic range.
The present invention provides novel apparatus and methods for fast quantitative measurement of perturbation of optical fields transmitted, reflected and or scattered along a length of an optical fibre.
The present invention can be used for distributed sensors, point sensors, or the combination of both.
In particular this technique can be applied to distributed sensors while extending dramatically the speed and sensitivity to allow the detection of acoustic perturbations anywhere along a length of an optical fibre while achieving fine spatial resolution. The present invention offers unique advantages in a broad range of acoustic sensing and imaging applications. Typical uses are for monitoring oil and gas wells, for applications such as for distributed flow metering and/or imaging; seismic imaging, monitoring long cables and pipelines; acoustic imaging inside large vessels as well as security applications.
It is an object of the present invention to provide apparatus for highly sensitive and fast quantitative measurement of the phase, frequency and amplitude of the light transmitted, reflected or scattered along a length of an optical fibre.
In the prior art, optical couplers have been used in Michelson or Mach-Zehnder interferometer configurations where the polarisation between the two arms of the interferometer has to be carefully controlled. The novel interferometer in the present invention allows an m×m coupler to be utilised using non-reciprocal devices, such as Faraday rotator mirrors and an optical circulator, to provide compensated light interference with a given phase shift that can be measured at all ports of the optical coupler and analysed very quickly, such as at several tens of kilohertz.
The embodiments of the invention can be used for multiplexed acoustic point sensors, distributed sensors or a combination of both. In the case of distributed sensors, light pulses are injected into the fibre and the phase modulation of the backscattered light is measured along the fibre at several tens of kilohertz. The fibre can be standard telecommunication fibre and/or cable. Using the techniques described herein, the sensing system can thereby detect the acoustic field along the fibre to provide a distributed acoustic sensor whereby the lengths of the sensing elements can be selected by a combination of adjusting the modulation of the light pulse, the path length in the interferometer as well as the sensing fibre configuration.
The data collected along the fibre are automatically synchronised and they may be combined to provide coherent field images.
Embodiments of the invention and how to put it into practice are described by way of example with reference to the accompanying drawings in which:—
First and second interference signal components are directed by the optical coupler 104 to photodetectors 113 and 114, which measure the intensity of the respective interference signal components.
The circulator 103 provides an efficient path for the input light and the returning (third) interference signal component through the same port of the coupler 104. The interference signal component incident on the optical circulator 103 is directed towards photodetector 115 to measure the intensity of the interference signal component.
The outputs of the photodetectors 113, 114 and 115 are combined to measure the relative phase of the incoming light, as described in more detail below with reference to
Optionally, frequency shifters 110 and 111 and/or optical modulator 109 may be used along the paths 105 and 106 for heterodyne signal processing. In addition, the frequency shift of 110 and 111 may be alternated from f1, f2 to f2, f1 respectively to reduce any frequency-dependent effect between the two portions of the light propagating through optical paths 105 and 106.
The above-described embodiment provides a novel apparatus suitable for fast quantitative measurement of perturbation of optical fields, and in particular can be used for distributed and multiplexed sensors with high sensitivity and fast response times to meet requirements of applications such as acoustic sensing.
In this embodiment 700, light emitted by a laser 701 is modulated by a pulse signal 702. An optical amplifier 705 is used to boost the pulsed laser light, and this is followed by a band-pass filter 706 to filter out the ASE noise of the amplifier. The optical signal is then sent to an optical circulator 707. An additional optical filter 708 may be used at one port of the circulator 707. The light is sent to sensing fibre 712, which is for example a single mode fibre or a multimode fibre deployed in an environment in which acoustic perturbations are desired to be monitored. A length of the fibre may be isolated and used as a reference section 710, for example in a “quiet” location. The reference section 710 may be formed between reflectors or a combination of beam splitters and reflectors 709 and 711.
The reflected and the backscattered light generated along the sensing fibre 712 is directed through the circulator 707 and into the interferometer 713. The detailed operation of the interferometer 713 is described earlier with reference to
It may be desirable to change the optical frequency of the light slightly to improve the sensitivity of the backscattered or reflected signals. The optical modulator 703 may be driven by a microwave frequency of around 10-40 GHz to generate optical carrier modulation sidebands. The optical filter 708 can be used to select the modulation sidebands which are shifted relative to the carrier. By changing the modulation frequency it is possible to rapidly modulate the selected optical frequency.
The three intensity measurements I1, I2, I3, from the photodetectors 113, 114, 115 are combined at step 1102 to calculate the relative phase and amplitude of the reflected or backscattered light from the sensing fibre. The relative phase is calculated (step 1104) at each sampling point, and the method employs oversampling such that more data points are available than are needed for the required spatial resolution of the system.
Methods for calculating the relative phase and amplitude from three phase shifted components of an interference signal are known from the literature. For example, Zhiqiang Zhao et al. and U.S. Pat. No. 5,946,429 describe techniques for demodulating the outputs of 3×3 couplers in continuous wave multiplexing applications. The described techniques can be applied to the time series data of the present embodiment.
For each sampling point, a visibility factor V is calculated at step 1106 from the three intensity measurements I1, I2, I3, from the photodetectors 113, 114, 115, according to equation (1), for each pulse.
V=(I1−I2)2+(I2−I3)2+(I3−I1)2 Equation (1)
At a point of low visibility, the intensity values at respective phase shifts are similar, and therefore the value of V is low. Characterising the sampling point according the V allows a weighted average of the phase angle to be determined (step 1108), weighted towards the sampling points with good visibility. This methodology improves the quality of the phase angle data 1110.
Optionally, the visibility factor V may also be used to adjust (step 1112) the timing of the digital sampling of the light for the maximum signal sensitivity positions. Such embodiments include a digitiser with dynamically varying clock cycles, (which may be referred to herein as “iclock”). The dynamically varying clock may be used to adjust the timing of the digitised samples at the photodetector outputs for the position of maximum signal sensitivity and or shifted away from positions where light signal fading occurs.
The phase angle data is sensitive to acoustic perturbations experienced by the sensing fibre. As the acoustic wave passes through the optical fibre, it causes the glass structure to contract and expand. This varies the optical path length between the backscattered light reflected from two locations in the fibre (i.e. the light propagating down the two paths in the interferometer), which is measured in the interferometer as a relative phase change. In this way, the optical phase angle data can be processed at 1114 to measure the acoustic signal at the point at which the light is generated.
In preferred embodiments of the invention, the data processing method 1100 is performed utilising a dedicated processor such as a Field Programmable Gate Array.
For accurate phase measurement, it is important to measure the offset signals and the relative gains of the photo-detectors 113,114 and 115. These can be measured and corrected for by method 1200, described with reference to
Each photodetector has electrical offset of the photodetectors, i.e. the voltage output of the photodetector when no light is incident on the photodetector (which may be referred to as a “zero-light level” offset. As a first step (at 1202) switching off the incoming light from the optical fibre and the optical amplifier 101. When switched off, the optical amplifier 101 acts as an efficient attenuator, allowing no significant light to reach the photodetectors. The outputs of the photodetectors are measured (step 1204) in this condition to determine the electrical offset, which forms a base level for the calibration.
The relative gains of the photodetectors can be measured, at step 1208, after switching on the optical amplifier 101 while the input light is switched off (step 1206). The in-band spontaneous emission (i.e. the Amplified Spontaneous Emission which falls within the band of the bandpass filter 102), which behaves as an incoherent light source, can then be used to determine normalisation and offset corrections (step 1210) to calibrate the combination of the coupling efficiency between the interferometer arms and the trans-impedance gains of the photodetectors 113, 114 and 115. This signal can also be used to measure the signal offset, which is caused by the in-band spontaneous emission.
Conveniently, the optical amplifier, which is a component of the interferometer, is used as in incoherent light source without a requirement for an auxiliary source. The incoherence of the source is necessary to avoid interference effects at the photodetectors, i.e. the coherence length of the light should be shorter than the optical path lengh of the interferometer. However, for accurate calibration it is preferable for the frequency band of the source to be close to, or centred around, the frequency of light from the light source. The bandpass filter 102 is therefore selected to filter out light with frequencies outside of the desired bandwidth from the Amplified Spontaneous Emission.
When used in a pulsed system, such as may be used in a distributed sensor, the above-described method can be used between optical pulses from the light source, to effectively calibrate the system during use, before each (or selected) pulses from the light source with substantively no interruption to the measurement process.
Variations to the above-described embodiments are within the scope of the invention, and some alternative embodiments are described below.
The 4×4 optical coupler of this arrangement generates four interference signal components at relative phase shifts of −90 degrees, 0 degrees, 90 degrees, 180 degrees.
In all of the above-described embodiments, optical switches may be used to change and/or select different combinations of optical path lengths through the interferometer. This facilitates switching between different spatial resolution measurements (corresponding to the selected path length differences in the optical path lengths).
The embodiments described above relate to apparatus and methods for fast quantitative measurement of acoustic perturbations of optical fields transmitted, reflected and or scattered along a length of an optical fibre. The invention in its various aspects can be applied or implemented in other ways, for example to monitor an optical signal generated by a laser, and/or to monitor the performance of a heterodyne signal generator, and to generate optical pulses for transmission into an optical signal. An example is described with reference to
In this embodiment, the reflected and backscattered light is not detected by an interferometer according to the invention. Rather, the reflected and backscattered light is passed through an optical amplifier 804 and an optical filter 806 and are then sent to a fast, low-noise photodetector 807. The electrical signal is split and then down-converted to baseband signals by mixing the RF signal 805 at different phase angles, in a manner known in the art. The electrical signals are digitised and the relative optical phase modulation at each section of the fibre is computed by combining the digitised signals using a fast processor 808.
The optical filter 708 selects two modulation frequency sidebands 1202/1203 and 1204/1205 generated by the optical modulator as indicated in
These arrangements are particularly useful to increase the sensing sensitivity, frequency response and the spatial resolution of the sensing system, while simplifying installation techniques and minimising bending losses.
Also, since we do not require any discrete components, a higher operating temperature can be achieved with proper coating protection applied on to the fibre. The fibre sensitivity can also be enhanced or reduced using different coatings or jackets. Also, the fibre can be made into a continuous cable with an enhanced sensing sensitivity while proving a protection for the fibre in harsh environments.
In one aspect, the invention provides an optical interferometer apparatus which can provide multiple path differences between the optical signals and provide interference signals between different optical paths with fixed and/or variable phase shifts. The interferometer utilises beam splitting components, circulating devices and Faraday rotator mirrors in a novel configuration. The optical signals at the output of the interferometer are converted to electrical signals which digitised for fast processing. The offset levels of the electrical signals are removed and their amplitude are normalised. The relative phase shifts of optical signals are accurately determined by combining the normalised electrical signals.
In another aspect, the invention relates to an interferometer apparatus that utilises beam splitters and non-reciprocal devices to provide light interference with given phase shifts and path length differences that can be measured at all ports of the beam splitters whereby the relative phase modulation of the light can be computed very accurately and quickly, such as at every few nanoseconds. The interferometer may use optical fibre components such as an m×m fused optical fibre coupler that is connected to an optical fibre circulator at one of its ports; Faraday-rotator mirrors that reflect and, at the same time, provide polarisation compensation for the light propagating through the different paths of the interferometer and photodetectors that are used to measure the interference light signals. The incoming optical light may be amplified using an optical fibre amplifier, and preferably the interferometer has a pass band optical filter to filter out the out of band Amplified Spontaneous Emission noise (ASE). The interferometer may provide birefringence compensation for light propagating along different optical paths through the interferometer. This provides sufficiently high visibility at the outputs of the interferometer.
In another of its aspects, the invention provides a method for compensating the offset and the gain of the photo-detectors, and the coupling ratio of the interferometer arms, to normalise the resultant interference signals used to measure the relative phase of the modulated input light in any of preceding claims where the detector offset is measured by switching off the optical amplifier in the backscatter path; the resultant photo-detector offset and gain then being determined by switching on the amplifier while the input light is switched off; the ASE of the optical amplifier then acts as an independent incoherent light source and thereby the offsets and relative gains of the photo-detectors can be determined and the detected light signals normalised. The method may therefore use incoherent light that enters the input of the interferometer to normalise the relative signal amplitudes at the output of the photo-detectors. For example, when an optical preamplifier is used at the input of the interferometer, the spontaneous light emission can be used to measure the combination of the splitting ratio of the interferometer arms and the relative gains of the photo-detectors and thereby normalise the relative signal amplitudes accordingly.
Another additional feature of the present invention is to use phase modulators and/or frequency shifters to shift the relative frequency and or vary the phase between the optical paths of the interferometer. Frequency shifters and/or phase modulators may be used to provide heterodyne signals and/or to separate the resultant interference light signal from different paths through the interferometer.
An additional feature of an embodiment of the invention is selecting the frequency of the frequency shifter sufficiently high so that at least one cycle of the beat frequency results within one light pulse resolution. Different frequency shifts may be used between different optical paths of the interferometer for the separation and/or heterodyne detection of the phase between different optical paths. The frequency shifts between different optical paths may be alternated to correct for any frequency dependency of the interferometer output signals.
An additional feature of an embodiment of the invention is the selection of different optical paths through the interferometer such as by using optical switches. The optical switches may be used to select different optical paths through the interferometer and thereby select a different spatial resolution measurement. Another aspect of the invention relates to a system comprising a number of interferometers cascaded in a series or in a star configuration or a combination of both.
The invention also provides a system that utilises a light pulse for multiplexed and/or distributed sensors by measuring the phase modulation of the reflected and/or the backscattered light along a length of fibre with high sensitivity, high dynamic range and a high speed of over tens of kilohertz. In this way, the invention can provide a multiplexed and/or distributed acoustic sensing system.
An additional feature of an embodiment of the invention is digitising the outputs of the interferometer, or the photodetectors of the interferometer, at least twice over a spatial resolution interval. An additional feature of an embodiment of the invention is combining the outputs of the interferometer to determine the insensitive measurement sample points resulting from any signal fading of the light in order to reject and/or provide a weighted signal average of the multiple samples of the light over a given spatial resolution measurement or interval. Embodiments of the invention use a digitiser with dynamically varying clock cycles, (which may be referred to herein as “iclock”), to adjust the timing of the digital sampling of the light for the maximum signal sensitivity positions. The dynamically varying clock may be used to adjust the timing of the digitised samples at the photo-detector outputs for the position of maximum signal sensitivity and or shifted away where light signal fading occurs.
A further aspect of the invention provides frequency shifted light, using a fast optical modulator to generate sidebands, preferably with a suppressed carrier spectrum, and a band-pass optical filter to select the modulation sidebands whereby the modulation frequency is varied rapidly between two portions of light pulse propagating through the optical modulator. The optical modulator may also chop off a portion of light pulse in the middle so as to generate two pulses with different frequencies. In this case the reflected and/or the backscattered light generated by the two pulses are combined to result in a heterodyne signal whose phase is determined to measure the relative optical phase modulation along the sensing fibre.
Providing multiple heterodyne signals can improve the dynamic range and reduce the effect of signal fading. When the scattered and/or the reflected light from the two pulses are combined, the modulation sidebands generate different beat frequencies which are proportional to the modulation frequency difference and to the order of the sidebands. The frequency of the light may be changed to optimise the signal sensitivity over a given section of the fibre. The frequency of the light passing through the optical modulator may be changed rapidly so that at least two portions of light pulse have different modulation sideband frequencies and, in addition, part of the light pulse may be chopped to generate two distinct portions of light pulses with different modulation sideband frequencies. The modulation sidebands between the two portions of the light pulse scattered or reflected from a sensing fibre may beat together to generate multiple heterodyne signals at multiples of the frequency difference between the two pulses that are proportional to the order of the modulation sidebands.
Embodiments of the invention may use a laser light or a broadband light source. Coherent matching of the light with the same delay results in an interference signal that can be used to measure the relative phase modulation of the scattered or reflected light along the fibre. The invention may use wavelength division multiplexed components to utilise multiple laser light pulses with different wavelengths and, preferably, varying time shift with respect to each to control the cross-phase modulation between the light pulses and to allow the processing of multiple pulses in the sensing fibre without and cross-sensitivity to allow the system to achieve a higher measure and frequency response. This may be the acoustic frequency response of the system to provide a different spatial sampling resolutions and/or positions, and/or to allow the efficient rejection of any points with low sensitivity.
An additional feature of an embodiment of the invention is the selection of different spatial resolutions whereby the sensitivity and the frequency response along the sensing fibre can be adjusted, and the dynamic range can be widened.
The sensing fibre may be standard single mode fibre, polarisation maintaining fibre, a single polarisation fibre, and or a ribbon fibre, and it can be coated and or cabled to enhance or to suppress its sensitivity.
An additional feature of an embodiment of the invention is the selection of different configurations of the fibre to optimise the sensitivity, the frequency and the directionality of the sensing fibre at different locations. The fibre may be deployed as linear sensors, direction sensors or multidimensional array sensors. The fibre may be placed on a surface area in a continuous path without crossing over another part of the fibre to increase the sensitivity, the frequency response and or the spatial resolution of the sensor system such as in a folded three-Omegas (ΩΩΩ) and or double eights (88) configurations. This is particularly useful to increase the sensing sensitivity, frequency response and the spatial resolution of the sensing system, while simplifying installation techniques and minimising bending losses.
The fibre may be attached on a surface of a vessel to listen to the noise generated within the vessel to monitor the changes in the process, acoustically image the process, as well to detect any leaks.
A further aspect provides an apparatus using acoustic sensors for distributed flow measurement and imaging, in-well perforated zones monitoring and sand production monitoring. For example, for in-well applications, the acoustic noise profile can be used to measure the flow by noise logging at every location along the well. In addition, the noise spectrum can be used to identify the phase of the fluid. Further noise spectrum correlation techniques can be used over a long section of the well to determine the speed of sound as well as tracking eddies generated within the flow to accurately determine the flow rates.
The sensor systems may be used as a distributed acoustic sensor, enabling the determination of distributed flow measurement and imaging, perforated zones monitoring and sand production monitoring in oil and gas wells and flowlines. The distributed temperature and strain measurements may be combined to enhance the data interpretation of the distributed acoustic sensor.
A further application is listening along previously installed optical fibres for surveillance applications. This includes measurements along fibres installed along boreholes, pipelines, perimeters, ports and borders.
An additional aspect provides a dynamic positioning apparatus using acoustic fibre optic sensors and acoustic reference sources whereby the optical fibre sensor measures the time of flight of acoustic signals received at different locations along the structure and thereby determines its position.
A further aspect provides pipeline structure monitoring apparatus using an acoustic fibre sensor and a pig that emits a sound (known as a “whistling pig”). The optical fibre sensor measures the acoustic transmission through the wall of the pipe for diagnostics as well as for tracking the position of the pig.
Another aspect provides pipeline monitoring apparatus where the sensing fibre is deployed inside the pipeline and carried along the pipeline by the fluid drag to provide a measurement of the noise flow for diagnostics of the pipeline as well as for flow characterisation and/or imaging.
Another aspect provides an apparatus using a fibre sensor used for acoustic sensing and an energy harvesting self-powered acoustic source to generate sufficient acoustic emission that can be picked up by a nearby sensing fibre for data communication, measurement, diagnostics and surveillance applications including along long pipelines, in-well and in other remote applications.
Another aspect of the invention provides an apparatus using acoustic fibre sensors to measure seepage rates along dams and dykes by generating an acoustic noise source in the upstream reservoir or in the core of the dam and measuring the acoustic signal strength detected along the fibre whereby areas of seepage act as low acoustic impedance paths for acoustic wave transmission and thereby exhibiting louder signal levels.
Other advantages and applications of the invention will be apparent to those skilled in the art. Any of the additional or optional features can be combined together and combined with any of the aspects, as would be apparent to those skilled in the art.
As has been described above, apparatus and methods for fast quantitative measurement of perturbations of optical fields transmitted, reflected and/or scattered along a length of an optical fibre. In particular, the invention can be used for distributed sensing while extending dramatically the speed and sensitivity to allow the detection of acoustic perturbations anywhere along a length of an optical fibre while achieving fine spatial resolution. The present invention offers unique advantages in a broad range of acoustic sensing and imaging applications. Typical uses are for monitoring oil and gas wells such as for distributed flow metering and/or imaging, monitoring long cables and pipelines, imaging of large vessels as well as security applications.
There follows a set of numbered features describing particular embodiments of the invention. Where a feature refers to another numbered feature then those features may be considered in combination.
| Number | Date | Country | Kind |
|---|---|---|---|
| 0908990.5 | May 2009 | GB | national |
| 0912051.0 | Jul 2009 | GB | national |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 17364595 | Jun 2021 | US |
| Child | 18370974 | US | |
| Parent | 16508522 | Jul 2019 | US |
| Child | 17364595 | US | |
| Parent | 15368945 | Dec 2016 | US |
| Child | 16508522 | US | |
| Parent | 13322449 | Nov 2011 | US |
| Child | 15368945 | US |
