Выделить слова: 


Патент США №

8319680

Автор(ы)

Sai

Дата выдачи

27 ноября 2012 г.


Radar liquid level detection using stepped frequency pulses



РЕФЕРАТ

The invention relates to a 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.


Авторы:

Bin Sai (Den Haag, NL)

Заявитель:

Enraf B.V. (Delft, NL)

ID семейства патентов

37067612

Номер заявки:

12/279,976

Дата регистрации:

19 февраля 2007 г.
PCT Filed: February 19, 2007
PCT No.: PCT/NL2007/000048
371(c)(1),(2),(4) Date: August 19, 2008
PCT Pub. No.: WO2007/111498
PCT Pub. Date: October 04, 2007

Приоритет зарубежной заявки

Feb 22, 2006 [NL]1031209


Класс патентной классификации США:

342/124; 73/290R; 73/291

Класс международной патентной классификации (МПК):

G01S 13/08

Класс совместной патентной классификации:

G01F 23/284 (20130101); G01S 13/34 (20130101); G01S 13/36 (20130101)

Область поиска:

342/124

Использованные источники

[Referenced By]

Патентные документы США

3217267November 1965Loposer
3337814August 1967Brase et al.
3353104November 1967Loposer
3579281May 1971Kam et al.
3789302January 1974Rearwin et al.
3903482September 1975Pausini et al.
4000476December 1976Walker et al.
4024464May 1977Underhill et al.
4027274May 1977Fukui et al.
4068199January 1978Madoff
4072947February 1978Johnson
4114110September 1978Nossen
4361801November 1982Meyer et al.
4451930May 1984Chapman et al.
4510461April 1985Dickes et al.
4516084May 1985Crowley
4567448January 1986Ikeda
4675617June 1987Martin
4691176September 1987Hsiung et al.
4800341January 1989Johnson
4823399April 1989George
4868494September 1989Ryder et al.
4928525May 1990Aderholt et al.
4958228September 1990Kutsuki
4972160November 1990Sylvain
5027526July 1991Crane
5034703July 1991Schumacher
5036291July 1991Marz
5052028September 1991Zwack
5210539May 1993Voyce
5270669December 1993Jokura
5400253March 1995O'Connor
5406842April 1995Locke
5428361June 1995Hightower et al.
5442359August 1995Rubin
5446416August 1995Lin et al.
5572167November 1996Alder et al.
5708424January 1998Orlando et al.
5734302March 1998Teng et al.
5773913June 1998Casselden
5774089June 1998Bamler et al.
5923284July 1999Artis et al.
5994905November 1999Franchville
6107957August 2000Cramer et al.
6114987September 2000Bjornholt
6374187April 2002Knight et al.
6404288June 2002Bletz et al.
6621449September 2003Kunert
6629458October 2003Westerling et al.
6636575October 2003Ott
6662649December 2003Knight et al.
6671500December 2003Damgaard et al.
6762634July 2004Hattori
6774732August 2004Harnishfeger et al.
6876621April 2005Ohuchida et al.
7135870November 2006Mohajer et al.
7891229February 2011Sai
2002/0101373August 2002Arndt et al.
2002/0183030December 2002Damgaard et al.
2003/0167839September 2003Burger et al.
2004/0196177October 2004Billington et al.
2004/0207477October 2004Gumm
2005/0052314March 2005Spanke et al.
2005/0190098September 2005Bridgelall et al.
2005/0206555September 2005Bridgelall et al.
2006/0044145March 2006Akerstrom et al.
2006/0044147March 2006Knox et al.
2006/0143000June 2006Setoguchi
2007/0046528March 2007Larsson et al.
2007/0109177May 2007Baath et al.
2010/0002912January 2010Solinsky
2010/0037673February 2010Wicht et al.
2010/0066589March 2010Sai et al.
2010/0070207March 2010Sai
2010/0070208March 2010Sai
2010/0070209March 2010Sai
2010/0175470July 2010Schrier et al.
2010/0241369September 2010Wicht et al.

Зарубежные патентные документы

296 08 551Jul., 1996DE
100 40 180Apr., 2002DE
10 2004 06144Jun., 2006DE
10 2005 019 0Oct., 2006DE
0 573 034Dec., 1993EP
0 887 658Dec., 1998EP
2 064 188Jun., 1981GB
2 342 995Apr., 2000GB
1032192Jul., 2006NL
WO 98/12514Mar., 1998WO
WO 2004/018978Mar., 2004WO
WO 2004/053521Jun., 2004WO
WO 2007/053007May., 2007WO
WO 2007/111498Oct., 2007WO
WO 2008/010702Jan., 2008WO
WO 2008/104967Sep., 2008WO
WO 2010/019427Feb., 2010WO

Другие источники


"HERMetic Sensor", Honeywell Enraf, (no month) 207, 5 pages. cited by other .
"Servo Gauge 854ATG" Product Sheet, Enraf B.V., May 2008, 4 pages. cited by other .
"Starrett Digitape 25", 1999 and possibly earlier, 4 pages. cited by other .
European Search Report dated Aug. 17, 2010 in connection with European Patent Application No. EP 10 15 5559. cited by other .
International Standard, Petroleum and liquid petroleum products--Measurement of level and temperature in storage tanks by automatic methods, Part 3: Measurement of level in pressurized storage tanks (non-refrigerated), 1st Edition, Nov. 15, 2002, 21 pages. cited by other .
"Level sensor", wikipedia.org, Aug. 2, 2008, 8 pages. cited by other .
Communication pursuant to Article 94(3) EPC dated Jan. 11, 2010 in Communication with European Patent Application No. EP 09169534. cited by other .
European Search Report dated Dec. 14, 2009 in connection with European Patent Application No. EP 09169534. cited by other .
Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration dated Oct. 28, 2008 in connection with PCT Patent Application No. PCT/NL2008/000196. cited by other .
Bai Sin, et al., "Advanced High Precision Radar Gauge for Industrial Applications", Proceedings of the 2006 CIE International Conference on Radar, vol. 1, Oct. 16, 2006, pp. 463-466. cited by other.

Главный эксперт: Brainard; Timothy A


ФОРМУЛА ИЗОБРЕТЕНИЯ



The invention claimed is:

1. A method for accurately determining a level of a liquid, the method comprising: emitting radar signals towards the liquid, the emitted radar signals having multiple frequencies and multiple first phases; receiving radar signals reflected from the liquid, the reflected radar signals having the multiple frequencies and multiple second phases; determining multiple phase differences between the emitted radar signals and the reflected radar signals; and deriving the level of the liquid based on an actual phase difference between the emitted radar signal having a first frequency and the reflected radar signal having the first frequency; wherein the actual phase difference is determined based on a measured phase difference between the emitted radar signal having the first frequency and the reflected radar signal having the first frequency; and wherein the measured phase difference is adjusted using a wrapping factor determined by: deriving a second level value based on a phase shift between phase differences associated with successively emitted and received radar signals having different frequencies; and deriving a value of the wrapping factor based on the second level value.

2. The method according to claim 1, wherein deriving the second level value comprises: determining a frequency difference between the successively emitted and received radar signals having the different frequencies; and using the phase shift and the frequency difference to determine the second level value.

3. The method according to claim 2, wherein the second level value is determined using a formula of: L'=.DELTA..psi..sub.1-2v/(4.pi..DELTA.f.sub.12) wherein L' represents the second level value; wherein v represents a velocity of the radar signals through a medium, wherein .DELTA.f.sub.12 represents the frequency difference between the first frequency of the radar signals and a second frequency of the radar signals; and wherein .DELTA..psi..sub.1-2 represents the phase shift between the phase difference associated with the first frequency of the radar signals and the phase difference associated with the second frequency of the radar signals.

4. The method according to claim 1, wherein the value of the wrapping factor is rounded to a nearest integer.

5. The method according to claim 1, wherein the level of the liquid is determined using a formula of: L=v(.DELTA..phi..sub.w+2k.pi.)/(4.pi.f) wherein L represents the level of the liquid; wherein v represents a velocity of the radar signals through a medium; wherein .DELTA..phi..sub.w represents the measured phase difference; wherein k represents the value of the wrapping factor; and wherein f represents the first frequency.

6. The method according to claim 1, wherein the level of the liquid is determined using a formula of: L=.DELTA..phi..sub.1tv/(4.pi.f.sub.1) wherein L represents the level of the liquid; wherein v represents a velocity of the radar signals through a medium; wherein f.sub.1 represents the first frequency; and wherein .DELTA..phi..sub.1t represents the actual phase difference.

7. The method according to claim 6, wherein the actual phase difference is determined using a formula of: .DELTA..phi..sub.1t=.DELTA..phi..sub.1+2k.pi. wherein .DELTA..phi..sub.1 represents the measured phase difference between the emitted radar signal having the first frequency and the reflected radar signal having the first frequency; and wherein k represents the value of the wrapping factor.

8. The method according to claim 1, wherein the value of the wrapping factor is determined using a formula of: k=|2fL'/v| wherein k represents the value of the wrapping factor; wherein f represents one of the frequencies of the radar signals; wherein L' represents the second level value; and wherein v represents a velocity of the radar signals through a medium.

9. A system for accurately determining a level of a liquid, the system comprising: a radar antenna configured to emit radar signals to the liquid and to receive radar signals reflected from the liquid, the emitted radar signals having multiple frequencies and multiple first phases, the reflected radar signals having the multiple frequencies and multiple second phases; and a processing unit configured to determine the level of the liquid based on the emitted radar signals and the reflected radar signals, wherein the processing unit is configured to determine the level of the liquid by: determining multiple phase differences between the emitted radar signals and the reflected radar signals; and deriving the level of the liquid based on an actual phase difference between the emitted radar signal having a first frequency and the reflected radar signal having the first frequency; wherein the processing unit is configured to determine the actual phase difference based on a measured phase difference between the emitted radar signal having the first frequency and the reflected radar signal having the first frequency; wherein the processing unit is configured to adjust the measured phase difference using a wrapping factor to generate the actual phase difference; and wherein the processing unit is configured to determine the wrapping factor by: deriving a second level value based on a phase shift between phase differences associated with successively emitted and received radar signals having different frequencies; and deriving a value of the wrapping factor based on the second level value.

10. The system according to claim 9, wherein the processing unit is configured to derive the second level value by: determining a frequency difference between the successively emitted and received radar signals having the different frequencies; and using the phase shift and the frequency difference to determine the second level value.

11. The system according to claim 9, wherein the processing unit is configured to determine the level of the liquid using an equation: L=v(.DELTA..phi..sub.w+2k.pi.)/(4.pi.f) wherein L represents the level of the liquid; wherein v represents a velocity of the radar signals through a medium; wherein .DELTA..phi..sub.w represents the measured phase difference; wherein k represents the value of the wrapping factor; and wherein f represents the first frequency.

12. The system according to claim 9, wherein the processing unit is configured to determine the level of the liquid using an equation: L=.DELTA..phi..sub.1tv/(4.pi.f.sub.1) wherein L represents the level of the liquid; wherein v represents a velocity of the radar signals through a medium; wherein f.sub.1 represents the first frequency; and wherein .DELTA..phi..sub.1t represents the actual phase difference.

13. The system according to claim 12, wherein the processing unit is configured to determine the actual phase difference using an equation: .DELTA..phi..sub.1t=.DELTA..phi..sub.1+2k.pi. wherein .DELTA..phi..sub.1 represents the measured phase difference between the emitted radar signal having the first frequency and the reflected radar signal having the first frequency; and wherein k represents the value of the wrapping factor.

14. The system according to claim 9, wherein the processing unit is configured to determine the second level value using an equation: L'=.DELTA..psi..sub.1-2v/(4.pi..DELTA.f.sub.12) wherein L' represents the second level value; wherein v represents a velocity of the radar signals through a medium; wherein .DELTA.f.sub.12 represents the frequency difference between the first frequency of the radar signals and a second frequency of the radar signals; and wherein .DELTA..psi..sub.1-2 represents the phase shift between the phase difference associated with the first frequency of the radar signals and the phase difference associated with the second frequency of the radar signals.

15. The system according to claim 9, wherein the processing unit is configured to determine the value of the wrapping factor using an equation: k=|2fL'/v| wherein k represents the value of the wrapping factor; wherein f represents one of the frequencies of the radar signals; wherein L' represents the second level value; and wherein v represents a velocity of the radar signals through a medium.

16. A device comprising: a processing unit configured to determine a level of a liquid based on emitted radar signals and reflected radar signals, the emitted radar signals having multiple frequencies and multiple first phases, the received radar signals having the multiple frequencies and multiple second phases; wherein the processing unit is configured to determine the level of the liquid by: determining multiple phase differences between the emitted radar signals and the reflected radar signals; and deriving the level of the liquid based on an actual phase difference between the emitted radar signal having a first frequency and the reflected radar signal having the first frequency; wherein the processing unit is configured to determine the actual phase difference based on a measured phase difference between the emitted radar signal having the first frequency and the reflected radar signal having the first frequency; wherein the processing unit is configured to adjust the measured phase difference using a wrapping factor to generate the actual phase difference; and wherein the processing unit is configured to determine the wrapping factor by: deriving a second level value based on a phase shift between phase differences associated with successively emitted and received radar signals having different frequencies; and deriving a value of the wrapping factor based on the second level value.

17. The device of claim 16, wherein the processing unit is configured to derive the second level value by: determining a frequency difference between the successively emitted and received radar signals having the different frequencies; and deriving the second level value based on the frequency difference and the phase shift.

18. The device of claim 16, wherein the processing unit is configured to determine the value of the wrapping factor using a formula of: k=|2fL'/v| wherein k represents the value of the wrapping factor; wherein f represents one of the frequencies of the radar signals; wherein L' represents the second level value; and wherein v represents a velocity of the radar signals through a medium.

19. A non-transitory machine readable medium containing instructions that when executed cause a processing unit to: determine multiple phase differences between emitted radar signals emitted towards a liquid and reflected radar signals reflected from the liquid, the emitted radar signals having multiple frequencies and multiple first phases, the reflected radar signals having the multiple frequencies and multiple second phases; and derive a level of the liquid based on an actual phase difference between the emitted radar signal having a first frequency and the reflected radar signal having the first frequency, the actual phase difference based on a measured phase difference between the emitted radar signal having the first frequency and the reflected radar signal having the first frequency; wherein the instructions that when executed cause the processing unit to derive the liquid level comprise instructions that when executed cause the processing unit to adjust the measured phase difference using a wrapping factor; and wherein instructions that when executed cause the processing unit to determine the wrapping factor comprise instructions that when executed cause the processing unit to: derive a second level value based on a phase shift between phase differences associated with successively emitted and received radar signals having different frequencies; and derive a value of the wrapping factor based on the second level value.

20. The machine readable medium of claim 19, wherein the instructions that when executed cause the processing unit to derive the value of the wrapping factor comprise instructions that when executed cause the processing unit to derive the value of the wrapping factor using a formula of: k=|2fL'/v| wherein k represents the value of the wrapping factor; wherein f represents one of the frequencies of the radar signals; wherein L' represents the second level value; and wherein v represents a velocity of the radar signals through a medium.

21. The machine readable medium of claim 19, wherein the instructions that when executed cause the processing unit to derive the second level value comprise instructions that when executed cause the processing unit to: determine a frequency difference between the successively emitted and received radar signals having the different frequencies; and use the phase shift and the frequency difference to determine the second level value.


ОПИСАНИЕ



CROSS-REFERENCE TO RELATED APPLICATION(S)

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. .sctn.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 .DELTA.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 i) emitting radar signals having frequencies f.sub.1, f.sub.2, . . . and phases .phi..sub.1i, .phi..sub.2i, . . . in time sequence to the liquid surface; ii) receiving the radar signals having frequencies f.sub.1, f.sub.2, . . . and phases .phi..sub.1i, .phi..sub.2i, . . . being reflected from the liquid surface in time sequence; iii) determining the phase difference .DELTA..phi..sub.1, .DELTA..phi..sub.2, . . . between the emitted radar signals and the reflected radar signals; iv) deriving the level L partially on the basis of the phase difference as determined.

More specifically, the method according to the invention is characterised in that step iv) comprises the steps of v) determining the frequency difference .DELTA.f.sub.1-2, . . . between the successive emitted radar signals having frequencies f.sub.1, f.sub.2, . . . ; vi) determining phase shift .DELTA..psi..sub.1-2, . . . between the successively determined phase differences .DELTA..phi..sub.1, .DELTA..phi..sub.2, . . . .

According to the invention, a less accurate level value L' is then derived in step vii on the basis of the frequency difference .DELTA.f.sub.1-2 and the phase shift .DELTA..psi..sub.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=.DELTA..phi..sub.1tv/(4.pi.f.sub.1), wherein v the velocity of the radar signal through the medium; f.sub.1 the frequency of the radar signal; .DELTA..phi..sub.1t the actual phase difference between the emitted radar signal and the reflected radar signal having frequency f.sub.1.

According to the method, an inaccurate level measurement is first carried out, wherein the inaccurate level L' is determined by: L'=.DELTA..psi..sub.12v/(4.pi..DELTA.f.sub.12), wherein v the velocity of the radar signal through the medium; .DELTA.f.sub.12 the frequency difference between the radar signals f.sub.1 and f.sub.2; .DELTA..psi..sub.12 the phase shift between the phase differences .DELTA..phi..sub.1 and .DELTA..phi..sub.2.

The actual phase difference between the emitted radar signal and the reflected radar signal having frequency f.sub.1 is determined by: .DELTA..phi..sub.1t=.DELTA..phi..sub.1+2k.pi., wherein .DELTA..phi..sub.1 the measured phase difference between the emitted radar signal and the reflected radar signal having frequency f.sub.1; k the number.

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 f.sub.1 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 f.sub.1, f.sub.2, . . . and phases .phi..sub.1i, .phi..sub.2i, . . . to the liquid surface in time sequence; as well as for receiving radar signals having frequencies f.sub.1, f.sub.2, . . . and phases .phi..sub.1j being reflected from the liquid surface in time sequence; wherein said means are arranged for determining the phase difference .DELTA..phi..sub.1, .DELTA..phi..sub.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 .DELTA.f.sub.1-2, . . . between the successive emitted radar signals having frequencies f.sub.1, f.sub.2, . . . and determining the phase shift .DELTA..psi..sub.1-2, . . . between the successively determined phase differences .DELTA..phi..sub.1, .DELTA..phi..sub.2, . . ., wherein said means are furthermore arranged for deriving a less accurate level value L' on the basis of the determined frequency difference .DELTA.f.sub.1-2 and the phase shift .DELTA..psi..sub.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.


КРАТКОЕ ОПИСАНИЕ РИСУНКОВ



FIG. 1 schematically shows a device for determining a level of a liquid in a tank;

FIGS. 2 and 3 illustrate radar signals used to determine a level of a liquid in a tank;

FIG. 4 illustrates a power spectral density (PSD) diagram of radar signals;

FIG. 5 illustrates operation of a Stepped Frequency-Continuous Wave (SF-CW) radar;

FIG. 6 illustrates phase indistinctness in radar signals;

FIG. 7 illustrates a number of periods of a sinus waveform during a frequency measurement;

FIGS. 8a-8b-8c illustrate obstacle test results obtained with different measuring methods;

FIGS. 9a-9b illustrate test results involving tank wall effects, obstacles and near effects obtained with different measuring methods; and

FIG. 10 illustrates test results involving tank bottom effects.

A well-known principle for level measurement is to use pulsed radar signals. FIG. 1 schematically shows a device for determining the level of the liquid in a tank in accordance with this known measuring principle. The device 10 is disposed in the upper part of a tank 1, which is built up of walls 1a, a roof 1b and a bottom 1c. The height of the tank 1 is indicated by the letter H.

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 FIG. 2) is partially reflected from the liquid surface 3 and the reflected radar signal 4b (also shown in FIG. 2) is in turn intercepted by the radar antenna 12. It is also possible, of course, to use a separate antenna for emitting the radar signals to the liquid surface and a separate receiver for intercepting the reflected radar signals.

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 .DELTA.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=1/2.v..DELTA.t (1)

wherein H=the height of the tank [m] L=the distance between the radar antenna and the liquid surface [m] h=the height of the liquid surface [m] v=the propagation velocity of the radar waves through the medium [m/sec] .DELTA.t=the time difference between the emitted radar signal and the reflected radar signal [sec]

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 FIG. 3. As a result of the time delay caused by the distance to be covered between the antenna and the target surface, a frequency difference f.sub.beat exists between the emitted radar signal 4a and the reflected radar signal 4b. Said frequency f.sub.beat can be calculated by Fourier transformation (see FIG. 4). The distance L can thus be determined for the triangular signal:

.times..times..DELTA..times..times..times..times. ##EQU00001##

wherein f.sub.m=the modulation frequency [Hz] .DELTA.F=`sweep` frequency band [Hz] f.sub.beat=the `beat` frequency between the emitted radar signal and the received radar signal [Hz] v=the propagation velocity of the radar waves through the medium [m/sec]

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 FIG. 4. An ordinary software algorithm may be used for determining the medium frequence f.sub.beat of the correctly selected peak, which corresponds to the reflection from the liquid surface.

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 FIG. 1) caused by obstacles, such as tank mixers (indicated at the 5 in FIG. 1 and provided with stirring elements 5a), the tank bottom (1c), the tank wall (1a), ladders, heat exchangers, said drive, and this may therefore lead to errors in the range of a few millimeters.

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 (FIG. 5). Using an SF-CW radar installation, the target distance can be established by determining the phase-distance relationship, or:

.DELTA..phi..times..pi..times..times..times..times. ##EQU00002##

wherein .DELTA..phi..sub.t=the phase difference between the reference signal and the reflected signal L=the distance between the radar antenna and the object to be measured [m] f=one of the `discrete` frequencies [Hz] v=the propagation velocity of the radar waves through the medium [m/sec]

This actual phase difference .DELTA..phi..sub.t cannot be measured directly by means of a radar instrument due to the phase indistinctness 2.pi. and indistinct disturbances in the phase characteristic. The measured phase difference is represented as .DELTA..phi..sub.w. Usually the phase signal is represented as a sinus signal, for example as cos(.DELTA..phi..sub.w) or sin(.DELTA..phi..sub.w). As sin(.DELTA..phi.)=sin(.DELTA..phi.+2k.pi.), 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` .DELTA..phi..sub.w. See also FIG. 6. In the case of a discrete signal this phase is represented as: .DELTA..phi..sub.w=.DELTA..phi..sub.t-2k.pi. (4)

In (4) k is an integer, such that 0.ltoreq..DELTA..phi..sub.w<2.pi.. 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:

.times..pi..times..times..times..DELTA..phi..times..times..times..pi. ##EQU00003##

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 .lamda. 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(.DELTA..phi..sub.w) signal is to determine the number of periods of the sinus waveform during a frequency measurement as shown in FIG. 7. This is the same as the Fourier transformation in the PSD method. Thus it can be derived that the following applies with regard to the distance:

.times. ##EQU00004##

wherein m is the number of periods of the signal. By using the amplitude-based (PSD) starting distance L.sub.PSD, the indistinctness factor of the number k can be determined on the basis of:

.times..times..times. ##EQU00005## 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 L.sub.conv can be obtained analogously to (5):

.times..pi..times..times..times..DELTA..phi..times..times..pi. ##EQU00006##

However, if the error in the PSD-based distance L.sub.PSD is greater than a quarter of the wavelength, i.e. 7.5 mm at 10 GHz, this will lead to an incorrect indistinctness factor k.sub.PSD and accordingly also to an incorrect phase distance L.sub.CONV 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 FIG. 1.

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 k.sub.PS. 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. L.sub.CONV=L.sub.0+.DELTA.L.sub.1+.DELTA.L.sub.2+ . . . +.DELTA.L.sub.j (9)

wherein L.sub.0=the starting distance [m] .DELTA.L.sub.1, .DELTA.L.sub.2=the difference distances between two successive measurements [m]

It will be understood that accumulating errors may become very large errors, even if only the one-off incorrect indistinctness factor is represented by L.sub.PSD. 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:

.times..pi..times..times..times..DELTA..phi. ##EQU00007##

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:

.times..pi..DELTA..times..times..times..DELTA..psi. ##EQU00008##

wherein .DELTA.f=f.sub.1-f.sub.2 and .DELTA..psi..sub.w=.DELTA..phi..sub.w(f.sub.1)-.DELTA..phi..sub.w(f.sub.2- ). The integer k can subsequently be determined by means of the rough level measurement L.sub.coarse, in accordance with:

.psi..times..times..times. ##EQU00009##

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 L.sub.INV:

.times..pi..times..times..times..DELTA..phi..times..psi..times..pi. ##EQU00010##

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 FIG. 1 further comprise an information processing unit 13a, which is arranged for determining whether the phase difference .DELTA..phi..sub.1, .DELTA..phi..sub.2, . . . between the emitted signals 4a and the reflected radar signals 4b and deriving the level L, partially on the basis of the determined phase difference, in accordance with the steps of the method according to the invention.

A few test results or measurement data are shown in FIGS. 8, 9 and 10. The obstacle test is shown in FIGS. 8a-8b-8c.

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 FIG. 1). These undesirable objects or obstacles may interfere electromagnetically with the normal target detection and target measurements.

To illustrate this, FIGS. 8a-8c show test results obtained with different measuring methods for determining the level in a storage tank: the known, amplitude-based method (FIG. 8a entitled "Level error when using the PSD-amplitude method L.sub.PSD"), the phase-based method (FIG. 8b entitled "Level error when using the normal phase method L.sub.CONV") and the method according to the invention (FIG. 8c entitled "Level error when using the method according to the invention L.sub.INV"). It is apparent that the accuracy and the repeatability the new phase-independent method as described in the present patent application are much higher than those obtained with the known phase-based and/or PSD-based methods. Its accuracy is about 50 times better than that achieved with the known method. The tests aimed at the tank wall effects, obstacles and near effects are shown in FIG. 9a L "INV--new method", in which the accuracy of the method according to the invention is compared with that of normal PSD-based methods (FIG. 9b entitled "L.sub.PSD--known Fourier FMCW method"). Its accuracy is about 55 times better than that achieved with the usual method.

The results of the test carried out to determine the effect of the bottom on the measurement are shown in FIG. 10 (entitled "Obstacle (heating element) near the bottom"). The empty tank is filled with a chemical liquid. The bottom reflection has a major effect on the liquid level measurement as carried out by means of the known method, whereas the effect on the level measurement as carried out by means of the method according to the invention is only minor.

* * * * *