The present application is related to Netherlands Patent Application No. 1031209 filed Feb. 22, 2006, entitled “METHOD AND DEVICE FOR ACCURATELY DETERMINING THE LEVEL L OF A LIQUID BY MEANS OF RADAR SIGNALS EMITTED TO THE LIQUID SURFACE AND RADAR SIGNALS REFLECTED FROM THE LIQUID SURFACE.” Netherlands Patent Application No. 1031209 is hereby incorporated by reference into the present disclosure as if fully set forth herein. The present application hereby claims priority under 35 U.S.C. §119(a) to Netherlands Patent Application No. 1031209.
The invention relates to a phase-based method for accurately determining the level L of a liquid by means of radar signals emitted to the liquid surface and radar signals reflected from the liquid surface.
The invention further relates to a device for accurately determining the level of a liquid by means of the method according to the invention, which device comprises at least a radar antenna disposed above the liquid for emitting radar signals to the liquid and receiving radar signals reflected from the liquid surface, as well as means for determining the liquid level on the basis of the emitted radar signals and the reflected radar signals.
Radar (Radio Detection And Ranging) is widely used for contactless distance measurement. A very well-known principle is the time difference method. According to this method, a radar antenna emits a radar signal which strikes an object, for example a liquid surface. The object reflects part of the emitted radar signal/wave back in the direction of the radar antenna, which receives the reflected radar signal/wave.
Although it is possible to use separate radar antennas for the emitted radar signal and for the reflected radar signal, it is common practice to use the same radar antenna for emitting as well as for receiving. The radar system measures the time difference Δt between the emitted radar signal and the received radar signal. If the velocity of the emitted radar signal is known, the distance to the surface of the liquid can be determined in a simple manner, using suitable measuring means.
A device as described in the introduction which employs the above measuring principle is frequently used for accurately determining the level of a liquid, for example water or oil, in a storage tank in the process industry or in an oil refinery by means of radar signals. The radar signals that are used are generally pulsed radar signals.
The methods that are currently known are in particular based on the phase difference between the emitted radar signal and the reflected radar signal. Unfortunately, the actual phase difference cannot be measured directly by the radar receivers that are currently used. Due to the unreliability and uncertainty as regards the measured phase difference, it is not possible to determine an unequivocal and reliable value of the actual level of the liquid.
The object of the method according to the invention is to obviate the above drawback and to introduce a more accurate measuring principle, in which the aforesaid inaccuracies in the measuring methods that are currently used are taken into account. In order to achieve that object, the method according to the invention comprises the steps of
More specifically, the method according to the invention is characterised in that step iv) comprises the steps of
According to the invention, a less accurate level value L′ is then derived in step vii on the basis of the frequency difference Δf1-2 and the phase shift Δψ1-2 as determined in steps v and vi.
Subsequently the actual phase difference can be determined by deriving the number k between the emitted radar signals and the reflected radar signals in step viii according to the invention on the basis of the less accurate level value L′ as determined in step vii, on the basis of which actual phase difference the level L of the liquid is determined with a high degree of accuracy on the basis of the number k as determined in step viii and the phase difference as determined in step iii.
This provides a more accurate picture of the actual liquid level than is possible with the usual measuring methods. The phase-based signal processing method makes it possible to realise a level measurement with a margin of error of less than 1 mm. The invention makes use of the phase differences in radar signals at different frequencies.
In accordance with the measuring method according to the invention, the level L is determined by:
L=Δφ1tv/(4πf1), wherein
According to the method, an inaccurate level measurement is first carried out, wherein the inaccurate level L′ is determined by:
L′=Δψ12v/(4πΔf12), wherein
The actual phase difference between the emitted radar signal and the reflected radar signal having frequency f1 is determined by:
Δφ1t=Δφ1+2kπ, wherein
The number k may be determined by means of the measuring method according to the invention on the basis of k=|2fL′/v|, after the determination of the number k the actual phase difference between the emitted radar signal and the reflected radar signal having frequency f1 can be determined, on the basis of which the actual level measurement L can be realised.
The device according to the invention is characterised in that the radar antenna is arranged for emitting radar signals having frequencies f1, f2, . . . and phases φ1i, φ2i, . . . to the liquid surface in time sequence; as well as for receiving radar signals having frequencies f1, f2, . . . and phases φ1j being reflected from the liquid surface in time sequence; wherein said means are arranged for determining the phase difference Δφ1, Δφ2, . . . between the emitted radar signal and the reflected radar signals and deriving the level L partially on the basis of the determined phase difference.
More specifically said means are arranged for determining the frequency difference Δf1-2, . . . between the successive emitted radar signals having frequencies f1, f2, . . . and determining the phase shift Δψ1-2, . . . between the successively determined phase differences Δφ1, Δφ2, . . . , wherein said means are furthermore arranged for deriving a less accurate level value L′ on the basis of the determined frequency difference Δf1-2 and the phase shift Δψ1-2, and for furthermore deriving the number k between the emitted radar signals and the reflected radar signals on the basis of the less accurate level value L′.
Analogously to the above-described measuring method, said means are according to the invention arranged for determining the level L on the basis of the number k as determined and the phase difference as determined.
The method and invention according to the invention will now be explained in more detail with reference to drawings.
As already set forth above, the invention relates to a phase-based method for determining the level of a liquid in a reliable and very accurate manner. The method uses radar signals for determining the level of a product that is stored in, for example, a storage tank, wherein the presence of all kinds of obstacles or parts in the tank does not interfere with the measuring method.
a-8b-8c illustrate obstacle test results obtained with different measuring methods;
a-9b illustrate test results involving tank wall effects, obstacles and near effects obtained with different measuring methods; and
A well-known principle for level measurement is to use pulsed radar signals.
Present in the tank 1 is an amount of liquid 2, and the height of the liquid level 3 is indicated by the letter L.
The device 10 has at least one radar antenna 12, which is provided with an emitting surface 11 for emitting a radar signal 4a in the direction of the liquid surface 3. The radar signal 4a (shown in
The device 10 according to the prior art is also provided with means for determining the liquid level 3 (L) on the basis of the emitted radar signal 4a and the reflected radar signal 4b, wherein the measuring system is based on the determination of the time difference Δt between the emitted signal pulse and the received signal pulse. Since the velocity of the radar signal is known, the distance to the measuring object or, in this case to the liquid surface, can be determined by
L=H−h=½.v.Δt (1)
wherein
A drawback of level measurement based on pulsed radar signals is that the time measurement between the emitted radar signal and the received radar signal must be very high. An inaccurate time measurement will inevitably result in an inaccurate level measurement. The measurement on the basis of pulsed radar signals depends inter alia on the pulse shape of the radar signal and possibly also on the pulse amplitude. In addition to that, reflections of the radar signal from objects other than the liquid surface interfere with this method.
A more extensive method uses frequency modulation (Frequency Modulation Continuous-Wave radar), wherein the frequency of the radar signal changes over time. The signal may for example be triangular in shape, as shown in
wherein
The FMCW technology does not have the drawbacks of the pulsed radar measurement. The time delay measurement has been substituted for a more accurate frequency count. It is possible to use a Fast Fourier Transformation (FFT) in the digital signal processing so as to obtain a power spectrum density (PSD) in which the power distribution, seen over the frequency band of the beat signal, is shown as a single peak. If reflections occur, as caused by the presence of objects (obstructions) in the path of the signals, the reflections will appear as multiple peaks in the PSD diagram, as also shown in
The above-described measuring method of the FMCW radar has a number of significant drawbacks, however. In the first place, heavy demands are made on the stability of the slope of the frequency “sweep”, which must be constant to a high degree. In the second place it is difficult to maintain the high linearity of the “sweep” shape, and thus the medium frequency of the frequency band is indistinct. In the third-place, the calculation by means of Fourier Transformation of the exact centre of the beat frequency is very sensitive to interference by reflections (indicated at 4′ in
For this purpose a measuring technique referred to as Stepped Frequency-Continuous Wave (SF-CW) has been developed. The SF-CW radar method emits and receives a series of sinusoidal signals at discrete frequencies that fill the frequency band, which is required in the measurement and the control (
wherein
This actual phase difference Δφt cannot be measured directly by means of a radar instrument due to the phase indistinctness 2π and indistinct disturbances in the phase characteristic. The measured phase difference is represented as Δφw. Usually the phase signal is represented as a sinus signal, for example as cos(Δφw) or sin(Δφw). As sin(Δφ)=sin(Δφ+2kπ), with k being an integer, accurate information is lost. Since the exact value k is not known, the measurement will be inaccurate as well. The factor k is the so-called ‘wrapping’ factor and results in a phase indistinctness ‘wrapped phase’ Δφw. See also
Δφw=Δφt−2kπ (4)
In (4) k is an integer, such that 0≦Δφw<2π. The inaccuracy caused by the value of the k-number is best illustrated by expressing the distance L as a function of the k-number:
For example, RADAR-based level measurement devices operate predominantly in the frequency bandwidth of 8-12.5 Ghz (the X-band). A frequency of 10 GHz corresponds to a wavelength λ that equals 30 mm in vacuum. If the number k changes by 1, this corresponds to a distance change of 15 mm. Consequently, the number k is of major importance for an accurate measurement of the level distance L.
In digitised data the phase indistinctness can be calculated at the central frequency or medium frequency of the frequency bandwidth. A usual method (the so-called PSD method, as discussed above) for determining the distance L from a cos(Δφw) signal is to determine the number of periods of the sinus waveform during a frequency measurement as shown in
wherein m is the number of periods of the signal. By using the amplitude-based (PSD) starting distance LPSD, the indistinctness factor of the number k can be determined on the basis of:
The factor int[..] represents a rounding-off factor, which rounds off k to the nearest integer. By means of the indistinctness factor ('wrapping factor') k thus obtained, the normal phase based distance Lconv can be obtained analogously to (5):
However, if the error in the PSD-based distance LPSD is greater than a quarter of the wavelength, i.e. 7.5 mm at 10 GHz, this will lead to an incorrect indistinctness factor kPSD and accordingly also to an incorrect phase distance LCONV with an error of a magnitude several times that of a half wavelength. This means that the accuracy of the normal phase-based level depends strongly on the PSD-based level. It is generally known that the PSD-based method is sensitive to various interferences. If an interfering object is present in the path of the radar beam, an error of a magnitude of tenths of millimeters in the PSD-based level measurement is likely to occur. The interfering objects or obstacles may be the wall 1a, the bottom 1c etc of the storage vessel 1, for example, as well as mixers (5-5a), heat coils or ladders. See
It has been found that an error of a few % near a peak in the PSD spectrum already suffices for generating an incorrect indistinctness factor kPS. In addition to that the interference from the nearby antenna area may also lead to significant errors if the level of the medium in the storage tank rises to a height near the antenna. Generally the usual PSD method is very sensitive, and this leads to an unstable and inaccurate calculation method. As a result, the error comprises a “level leap” of half the wavelength in the normal “phase-based” method, which is 15 mm in the X-bandwidth. Such level leaps are highly undesirable in some applications.
Moreover, some known FMCW and/or SFCW phase-based measuring methods use the relative phase measurement for the purpose of correcting for the distance change between two successive measurements.
LCONV=L0+ΔL1+ΔL2+ . . . +ΔLj (9)
wherein
It will be understood that accumulating errors may become very large errors, even if only the one-off incorrect indistinctness factor is represented by LPSD. Thus the vulnerability of the normal phase-based distance calculation is affected by an incorrect PSD-based method, resulting in a poor accuracy performance.
The method according to the invention is aimed at neutralising this measuring error. From the above equation (3) it has already been determined that if the actual phase at one or more frequencies is known, the absolute distance to the target (L) can be determined on the basis of:
A known characteristic of the digitally controlled SF-CW radar technology is that every step frequency that is generated is known. According to the invention, the method uses the phase variations at different frequencies for solving the phase indistinctness of a half wavelength and thus accurately determining the absolute distance L. The method according to the invention relates to the use of the measured or formed phase at two different frequencies for determining the rough distance on the basis of:
wherein Δf=f1−f2 and Δψw=Δφw(f1)−Δφw(f2). The integer k can subsequently be determined by means of the rough level measurement Lcoarse, in accordance with:
In this case, too, the factor int[..] indicates a rounding-off factor, which rounds off k to the nearest integer. Thus the phase and the number can be utilised for determining whether the absolute phase between the emitted radar signal and the reflected signal and consequently also for determining the absolute, very accurate level distance LINV:
The measuring method according to the invention is characterised as an independent, phase-based signal processing method. Using the above-described method, very reliable and repeatable accuracies of ±1 mm or better can be realised, even in complex measuring conditions.
Consequently, the method according to the invention does not use the known, less accurate PSD as a reference for the phase indistinctness. In contrast to the known methods, the method according to the invention does not make use of the relative phase distance change with respect to the preceding measurement that influences the current distance. It calculates the absolute phase and distance with every measurement, which gives an absolute measure of the current target distance. Thus, the phase error accumulation from preceding measurements is entirely prevented by using this method.
According to the invention the level determining means 13 of the device 10 shown in
A few test results or measurement data are shown in
This so-called obstacle test is carried out in the level measurement to determine the existence of undesirable objects that will appear in the radar beam during the target measurement. Said undesirable objects may be so-called tank mixers 5-5a, ladders, heat coils, the tank bottom 1c, the tank wall 1a, etcetera (see
To illustrate this,
The results of the test carried out to determine the effect of the bottom on the measurement are shown in
Number | Date | Country | Kind |
---|---|---|---|
1031209 | Feb 2006 | NL | national |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/NL2007/000048 | 2/19/2007 | WO | 00 | 8/19/2008 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2007/111498 | 10/4/2007 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
3217267 | Loposer | Nov 1965 | A |
3337814 | Brase et al. | Aug 1967 | A |
3353104 | Loposer | Nov 1967 | A |
3579281 | Kam et al. | May 1971 | A |
3789302 | Rearwin et al. | Jan 1974 | A |
3903482 | Pausini et al. | Sep 1975 | A |
4000476 | Walker et al. | Dec 1976 | A |
4024464 | Underhill et al. | May 1977 | A |
4027274 | Fukui et al. | May 1977 | A |
4068199 | Madoff | Jan 1978 | A |
4072947 | Johnson | Feb 1978 | A |
4114110 | Nossen | Sep 1978 | A |
4361801 | Meyer et al. | Nov 1982 | A |
4451930 | Chapman et al. | May 1984 | A |
4510461 | Dickes et al. | Apr 1985 | A |
4516084 | Crowley | May 1985 | A |
4567448 | Ikeda | Jan 1986 | A |
4675617 | Martin | Jun 1987 | A |
4691176 | Hsiung et al. | Sep 1987 | A |
4800341 | Johnson | Jan 1989 | A |
4823399 | George | Apr 1989 | A |
4868494 | Ryder et al. | Sep 1989 | A |
4928525 | Aderholt et al. | May 1990 | A |
4958228 | Kutsuki | Sep 1990 | A |
4972160 | Sylvain | Nov 1990 | A |
5027526 | Crane | Jul 1991 | A |
5034703 | Schumacher | Jul 1991 | A |
5036291 | Marz | Jul 1991 | A |
5052028 | Zwack | Sep 1991 | A |
5210539 | Voyce | May 1993 | A |
5270669 | Jokura | Dec 1993 | A |
5400253 | O'Connor | Mar 1995 | A |
5406842 | Locke | Apr 1995 | A |
5428361 | Hightower et al. | Jun 1995 | A |
5442359 | Rubin | Aug 1995 | A |
5446416 | Lin et al. | Aug 1995 | A |
5572167 | Alder et al. | Nov 1996 | A |
5708424 | Orlando et al. | Jan 1998 | A |
5734302 | Teng et al. | Mar 1998 | A |
5773913 | Casselden | Jun 1998 | A |
5774089 | Bamler et al. | Jun 1998 | A |
5923284 | Artis et al. | Jul 1999 | A |
5994905 | Franchville | Nov 1999 | A |
6107957 | Cramer et al. | Aug 2000 | A |
6114987 | Bjornholt | Sep 2000 | A |
6374187 | Knight et al. | Apr 2002 | B1 |
6404288 | Bletz et al. | Jun 2002 | B1 |
6621449 | Kunert | Sep 2003 | B1 |
6629458 | Westerling et al. | Oct 2003 | B1 |
6636575 | Ott | Oct 2003 | B1 |
6662649 | Knight et al. | Dec 2003 | B1 |
6671500 | Damgaard et al. | Dec 2003 | B2 |
6762634 | Hattori | Jul 2004 | B1 |
6774732 | Harnishfeger et al. | Aug 2004 | B1 |
6876621 | Ohuchida et al. | Apr 2005 | B2 |
7135870 | Mohajer et al. | Nov 2006 | B2 |
7891229 | Sai | Feb 2011 | B2 |
20020101373 | Arndt et al. | Aug 2002 | A1 |
20020183030 | Damgaard et al. | Dec 2002 | A1 |
20030167839 | Burger et al. | Sep 2003 | A1 |
20040196177 | Billington et al. | Oct 2004 | A1 |
20040207477 | Gumm | Oct 2004 | A1 |
20050052314 | Spanke et al. | Mar 2005 | A1 |
20050190098 | Bridgelall et al. | Sep 2005 | A1 |
20050206555 | Bridgelall et al. | Sep 2005 | A1 |
20060044145 | Akerstrom et al. | Mar 2006 | A1 |
20060044147 | Knox et al. | Mar 2006 | A1 |
20060143000 | Setoguchi | Jun 2006 | A1 |
20070046528 | Larsson et al. | Mar 2007 | A1 |
20070109177 | Baath et al. | May 2007 | A1 |
20100002912 | Solinsky | Jan 2010 | A1 |
20100037673 | Wicht et al. | Feb 2010 | A1 |
20100066589 | Sai et al. | Mar 2010 | A1 |
20100070207 | Sai | Mar 2010 | A1 |
20100070208 | Sai | Mar 2010 | A1 |
20100070209 | Sai | Mar 2010 | A1 |
20100175470 | Schrier et al. | Jul 2010 | A1 |
20100241369 | Wicht et al. | Sep 2010 | A1 |
Number | Date | Country |
---|---|---|
296 08 551 | Jul 1996 | DE |
100 40 180 | Apr 2002 | DE |
10 2004 06144 | Jun 2006 | DE |
10 2005 019 0 | Oct 2006 | DE |
0 573 034 | Dec 1993 | EP |
0 887 658 | Dec 1998 | EP |
2 064 188 | Jun 1981 | GB |
2 342 995 | Apr 2000 | GB |
1032192 | Jul 2006 | NL |
WO 9812514 | Mar 1998 | WO |
WO 2004018978 | Mar 2004 | WO |
WO 2004053521 | Jun 2004 | WO |
WO 2007053007 | May 2007 | WO |
WO 2007111498 | Oct 2007 | WO |
WO 2008010702 | Jan 2008 | WO |
WO 2008104967 | Sep 2008 | WO |
WO 2010019427 | Feb 2010 | WO |
Number | Date | Country | |
---|---|---|---|
20110163910 A1 | Jul 2011 | US |