WO2016053978A1 - Level finding using multiple search steps - Google Patents
Level finding using multiple search steps Download PDFInfo
- Publication number
- WO2016053978A1 WO2016053978A1 PCT/US2015/052846 US2015052846W WO2016053978A1 WO 2016053978 A1 WO2016053978 A1 WO 2016053978A1 US 2015052846 W US2015052846 W US 2015052846W WO 2016053978 A1 WO2016053978 A1 WO 2016053978A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- level
- pulse
- waveguide
- minimum
- echo
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
- G01F23/22—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
- G01F23/28—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
- G01F23/284—Electromagnetic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
- G01F23/22—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
- G01F23/28—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
- G01F23/296—Acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/523—Details of pulse systems
- G01S7/526—Receivers
- G01S7/527—Extracting wanted echo signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/483—Details of pulse systems
- G01S7/486—Receivers
- G01S7/4865—Time delay measurement, e.g. time-of-flight measurement, time of arrival measurement or determining the exact position of a peak
- G01S7/4866—Time delay measurement, e.g. time-of-flight measurement, time of arrival measurement or determining the exact position of a peak by fitting a model or function to the received signal
Definitions
- Disclosed embodiments relate to level finding using time domain reflectometry (TDR) and level measurement gauges and systems therefrom.
- TDR time domain reflectometry
- Level measurement technologies include ultrasound and radar typically being
- GWR Guided Wave Radar
- GWR works by generating a stream of pulses of electromagnetic energy and propagating transmitted pulses down a transmission line formed into a level sensing probe.
- the probe is generally placed vertically in a tank or other container and the electromagnetic pulse is launched downward from the top of the probe.
- the probe is open to both the air and the material to be sensed in such a way that the electromagnetic fields of the propagating pulse penetrate the air until, they reach the level of the material.
- the electromagnetic fields see the higher dielectric constant of the material.
- This higher dielectric constant causes a reduction in the impedance of the transmission line, resulting in a pulsed reflected (echo) signal back to the top of the probe.
- the pulse travels through the air dielectric portion of the probe at a known velocity. This allows the material level and an interface thickness (if present) to be determined by a peak finding algorithms that identifies peak(s) in the pulse echo curve.
- Disclosed embodiments recognize known peak finding algorithms used to find peaks in echo curves for a time domain reflectometry (TDR)-based level measurement system for determining the distance to a material level and interface thickness (if present) operate without a-priori knowledge of the transmitted pulse being used, the physical properties of the waveguide, and medium in the tank being into account. Such methods then leave it to the user to then define a threshold array to pick out peaks of interest in the echo curve to determine the material level and interface(s) if present.
- TDR time domain reflectometry
- Disclosed multiple search step level finding algorithms automatically determine the level of a material (material level) in a tank and interface(s) if present from a reflection signal (echo curve) and a sampled pulse model echo, which generally include two search steps, an initial coarse search, then and a fine search.
- Level readings provided by disclosed level finding algorithms are more accurate and require less expert knowledge of the user as compared to known peak finding algorithms, in one embodiment being for relative!' difficult guided wave radar (GWR) applications such as for measuring an interface or for determining the material level in turbulent conditions.
- GWR difficult guided wave radar
- FIG. 1 is a flow chart that shows steps in an example method of level finding, according to an example embodiment
- FIG. 2 shows a measured versus a modelled radar reflection pulse, according to an example embodiment.
- FIG. 3 depicts an example pulsed radar system shown as a GWR system that a disclosed level finding algorithm which is shown implemented in the firmware of a memory associated with a processor, according to an example embodiment.
- FIG. 4 is an example illustration of a measured echo curve and a parametric model of the surface reflection, according to an example embodiment.
- FIG. 5A is an overlay of an example surface reflection model and measured echo curve
- FIG. 5B shows an example objective function and its threshold, according to example embodiments.
- Coupled to or “couples with” (and the like) as used herein without further qualification are intended to describe either an indirect or direct electrical connection.
- a first device “couples” to a second device, that connection can be through a direct electrical connection where there are only parasitics in the pathway, or through an indirect electrical connection via intervening items including other devices and connections.
- the intervening item generally does not modif the information of a signal but may adjust its current level, voltage level, and/or power level.
- FIG. 1 is a flow chart that shows steps in an example method 100 of level finding, according to an example embodiment.
- Disclosed embodiments determine the material level in a tank including two search steps from a reflection signal (echo curve) and a sampled pulse model echo, using an initial coarse search, then and a finer search. Although generally not needed, third or more successively finer searches may also be added.
- level finding is generally described for GWR applications, disclosed embodiments can also be applied to ultrasound and non-contacting radar.
- Step 101 comprises providing parameters including characteristics of a shape of a transmitted pulse (modeled or actual) in the time domain (e.g., pulse width and amplitude) launched by a transmitter onto a waveguide (or probe) into a tank having at least one material therem, physical properties of the waveguide, and dielectric characteristics of the material.
- a shape of a transmitted pulse modeled or actual
- the time domain e.g., pulse width and amplitude
- Step 102 comprises using a processor having these parameters supplied, implementing a level finding algorithm having a coarse search and then a fine search that is stored in an associated memory.
- the processor can comprise a microprocessor, microcontroller, digital signal processor (DSP), field programmable gate array (FPGA), application specific integrated circuit (ASIC), or discrete logic devices.
- the memory can be on the processor chip, or can be a separate memory.
- k is the discrete sample distance
- k € f3 ⁇ 4J ⁇ ' — I] is an integer
- ⁇ is the sample distance
- p is the sampledinstalle model echo.
- Step 104 comprises the fine search calculating at least one minimum or maximum using the objective function K3 ⁇ 43 in the vicinity of .3 ⁇ 4' , wherein the minimum or maximum %* corresponds to a level of the material or an interface involving the material in sub-sample resolution where the sample resolution is given by 3 ⁇ 4.,unless where the center of the reflections, not peaks, are identified.
- step 105 can be added comprising using a selection criteria (such as the minimum value of the objective function) to select the best candidate in the set of candidates.
- Step 101 to 104 and optionally step 105 are evaluated for all echo reflections of interest.
- the distance to the upper product surface ⁇ g ⁇ - eS is calculated as the difference between the best reflection candidate for the upper product and the reference plane: simi larly for lower product and end of
- a parameterized model or a measured nearly-perfectly- reflected ultrasound or radar pulse (e.g., transmitted pulse measured after passing through only the receive circuitry of the gauge) is obtained.
- g is a gain or amplitude parameter
- w is a width parameter
- m&3 ⁇ 4s is the discrete sample distance
- >&3 ⁇ 4 is the sample distance.
- An example of a parametrized model radar pulse reflection versus the actual measured radar pulse reflection is shown in Error! Reference source not found..
- the gain or amplitude ⁇ (shown in FIG. 2 as the gain) is governed by the amplitude/gain of the transmitted pulse and the dielectric characteristics of the materia! in the tank that the reflected pulse travels through.
- a is a function of the waveguide and the material in the tank in which the pulse is travelling.
- F&* smm models the radar or ultrasound pulse energy that is lost to its surroundings or is lost due to being absorbed by the waveguide.
- Reflection candidates are found in the measured echo curve by evaluating a performance index / ⁇ 3 ⁇ 4. ⁇ that is a norm of the prediction error vector For an example with the Euclidean norm, see in eq.
- ⁇ ⁇ is the sample distance
- M is the half window width of the reflection model
- 3 ⁇ 4 is the set of discrete distances for the reflection candidates
- & is the performance index threshold
- N is number of samples in the echo curve. It is noted that norms other than the Euclidean norm can be used in this step.
- are optionally down sampled or decimated for efficiency in the coarse search.
- the objective function 3 ⁇ 43 at is then subject to a fine search.
- An example fine search can include:
- an additional candidate selection criterion can be applied as described above relative to step 105.
- the selection criterion may be different depending on the reflection in question, such as whether reflection in question is from the reference plane, upper product surface, lower product surface (interface), or the end of probe.
- the selection criterion corresponds to the minimum value of the objective function:
- selection criteria can be used, such as using prior knowledge of where the level was in previous time intervals to make the best decision. Instead of the lowest minimum, one can select the first minimum below threshold as it is likely a level as opposed to a second reflection. One can also reject reflections at distance/time integer multiples of the first reflection as they are likely false.
- FIG. 3 depicts an example pulsed radar or ultrasound level finding system shown as a GW system 300 that includes a pulsed radar or ultrasound level gauge shown as a pulse radar level gauge (PRO) 360 including a disclosed level finding algorithm 310a is shown implemented in the firmware of memory 310 associated with a processor 315, according to an example embodiment. Also shown is a transceiver 320 and coaxial connector 325 that is on the top of the tank 305. The transmitter and the receiver provided by transceiver 320 may be implemented as separate blocks. Accordingly, a transceiver as used herein includes both of these arrangements.
- a transceiver as used herein includes both of these arrangements.
- the processor 315 provides digital signal levels to a digital to analog converter (DAC) 322 which is connected to an input of a pulse generator (PGen) block 321 that is coupled to the transmitter of the transceiver 320.
- DAC digital to analog converter
- Pen pulse generator
- the receiver of the transceiver 220 receives reflected echo signals that are transduced by the sensor 326, where the output signal from the sensor 326 is coupled to an analog-to-digital converter (ADC) 327 which converts analog signals from the sensor 326 into digital signals for the processor 315 to process which functions as a signal analyzer. Many pulses are generated and the ADC 327 generates the waveform comprising many discrete voltage values that are synchronized by a clock circuit.
- a flange may also be present.
- the waveguide (or probe) is shown as 340. As noted above, although generally described for GWR applications, disclosed level finding can also be applied to ultrasound and non-contacting radar.
- FIG, 4 illustrates results from an example parametric model and an actual measured echo curve. This example also shows the found reference reflection and the surface (upper product) reflection for a single liquid application, the linear attenuation ⁇ # ' ⁇ ⁇ 8 3 ⁇ 4
- model gain g has a step change at the surface boundary due to the difference in dielectric constant between the vapor and the measured material .
- FIG. 5A shows the upper product surface reflection model overlaid on the measured echo curve at the location of the upper product surface.
- FIG. 5B shows the corresponding objective function and its threshold.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Fluid Mechanics (AREA)
- Thermal Sciences (AREA)
- Acoustics & Sound (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Levels Of Liquids Or Fluent Solid Materials (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
A method (100) of level finding includes providing (101) characteristics of a shape of a transmitted pulse in time domain launched onto a waveguide (340) into a tank (305) having at least one material therein, physical properties of the waveguide and real and imaginary dielectric characteristics of the material at a frequency of the pulse. A level finding algorithm (310a) having a coarse search and a fine search is implemented, where the coarse search minimizes (103) a prediction error between an echo signal (echo curve y(k)) and a sampled pulse model echo p(k) to obtain an objective function J(k) in a vicinity of a minimum prediction error (k*). The fine search calculates (104) at least one minimum or maximum using J(k) in the vicinity of k*. The minimum or the maximum corresponds to a level of the material or an interface involving the material.
Description
LEVEL FINDING USING MULTIPLE SEARCH STEPS
CROSS REFERENCE TO RELATED APPLICATIONS
[0001] This application claims the benefit of Provisional Application Serial No.
62/058,358 entitled "TIME DOMAIN REFLECTOMETRY (TDR) SIGNAL PROCESSING", filed October 1, 2014, which is herein incorporated by reference in its entirety,
FIELD
[0002] Disclosed embodiments relate to level finding using time domain reflectometry (TDR) and level measurement gauges and systems therefrom.
BACKGROUND
[0Θ03] Level measurement technologies include ultrasound and radar typically being
Guided Wave Radar (GWR). Newer level measurement devices are TDR-based that utilize sound or echo-based measurement whether ultrasound or GWR. GWR is a contact radar technique to measure the level of liquids or solids in a tank using time TDR principles.
[0004] GWR works by generating a stream of pulses of electromagnetic energy and propagating transmitted pulses down a transmission line formed into a level sensing probe. The probe is generally placed vertically in a tank or other container and the electromagnetic pulse is launched downward from the top of the probe. The probe is open to both the air and the material to be sensed in such a way that the electromagnetic fields of the propagating pulse penetrate the air until, they reach the level of the material. At that point, the electromagnetic fields see the higher dielectric constant of the material. This higher dielectric constant causes a reduction in the impedance of the transmission line, resulting in a pulsed reflected (echo) signal back to the top of the probe. The pulse travels through the air dielectric portion of the probe at a known velocity. This allows the material level and an
interface thickness (if present) to be determined by a peak finding algorithms that identifies peak(s) in the pulse echo curve.
SUMMARY
[0005] This Summar is provided to introduce a brief selection of disclosed concepts in a simplified form that are further described below in the Detailed Description including the drawings provided. This Summary is not intended to limit the claimed subject matter's scope.
[0006J Disclosed embodiments recognize known peak finding algorithms used to find peaks in echo curves for a time domain reflectometry (TDR)-based level measurement system for determining the distance to a material level and interface thickness (if present) operate without a-priori knowledge of the transmitted pulse being used, the physical properties of the waveguide, and medium in the tank being into account. Such methods then leave it to the user to then define a threshold array to pick out peaks of interest in the echo curve to determine the material level and interface(s) if present.
[0007] Disclosed multiple search step level finding algorithms automatically determine the level of a material (material level) in a tank and interface(s) if present from a reflection signal (echo curve) and a sampled pulse model echo, which generally include two search steps, an initial coarse search, then and a fine search. Level readings provided by disclosed level finding algorithms are more accurate and require less expert knowledge of the user as compared to known peak finding algorithms, in one embodiment being for relative!)' difficult guided wave radar (GWR) applications such as for measuring an interface or for determining the material level in turbulent conditions. Although generally described for GWR applications (contact pulsed radar), disclosed level finding can also be applied to ultrasound and non-contacting radar.
BRIEF DESCRIPTION OF THE DRAWINGS
[0008] FIG. 1 is a flow chart that shows steps in an example method of level finding, according to an example embodiment,
[0009] FIG. 2 shows a measured versus a modelled radar reflection pulse, according to an example embodiment.
[0010] FIG. 3 depicts an example pulsed radar system shown as a GWR system that a disclosed level finding algorithm which is shown implemented in the firmware of a memory associated with a processor, according to an example embodiment.
[0011] FIG. 4 is an example illustration of a measured echo curve and a parametric model of the surface reflection, according to an example embodiment.
[0012] FIG. 5A is an overlay of an example surface reflection model and measured echo curve, and FIG. 5B shows an example objective function and its threshold, according to example embodiments.
DETAILED DESCRIPTION
[0013] Disclosed embodiments are described with reference to the attached figures, wherein like reference numerals are used throughout the figures to designate similar or equivalent elements. The figures are not drawn to scale and they are provided merely to illustrate certain disclosed aspects. Several disclosed aspects are described below with reference to example applications for illustration. It should be understood that numerous specific details, relationships, and methods are set forth to provide a full understanding of the disclosed embodiments.
[0014] One having ordinar skill in the relevant art, however, will read ly recognize that the subject matter disclosed herein can be practiced without one or more of the specific details or with other methods. In other instances, well-known structures or operations are not shown in detail to avoid obscuring certain aspects. This Disclosure is not limited by the illustrated ordering of acts or events, as some acts may occur in different orders and/or
concurrently with other acts or events. Furthermore, not all illustrated acts or events are required to implement a methodology in accordance with the embodiments disclosed herein.
[0015] Also, the terms "coupled to" or "couples with" (and the like) as used herein without further qualification are intended to describe either an indirect or direct electrical connection. Thus, if a first device "couples" to a second device, that connection can be through a direct electrical connection where there are only parasitics in the pathway, or through an indirect electrical connection via intervening items including other devices and connections. For indirect coupling, the intervening item generally does not modif the information of a signal but may adjust its current level, voltage level, and/or power level.
[0Θ16] FIG. 1 is a flow chart that shows steps in an example method 100 of level finding, according to an example embodiment. Disclosed embodiments determine the material level in a tank including two search steps from a reflection signal (echo curve) and a sampled pulse model echo, using an initial coarse search, then and a finer search. Although generally not needed, third or more successively finer searches may also be added. As noted above, although disclosed level finding is generally described for GWR applications, disclosed embodiments can also be applied to ultrasound and non-contacting radar.
[0017] Step 101 comprises providing parameters including characteristics of a shape of a transmitted pulse (modeled or actual) in the time domain (e.g., pulse width and amplitude) launched by a transmitter onto a waveguide (or probe) into a tank having at least one material therem, physical properties of the waveguide, and dielectric characteristics of the material. To obtain the dielectric characteristics of the material one can simulate or more typically measure the dielectric characteristics of the material in the tank at the desired operating pulse frequency or frequencies. Step 102 comprises using a processor having these parameters supplied, implementing a level finding algorithm having a coarse search and then a fine search that is stored in an associated memory. The processor can comprise a
microprocessor, microcontroller, digital signal processor (DSP), field programmable gate array (FPGA), application specific integrated circuit (ASIC), or discrete logic devices. The memory can be on the processor chip, or can be a separate memory.
[0018] Step 103 comprises the coarse search where reflection candidates ar are found in the measured echo curve } (¾) by evaluating a performance index /£%} that can be a norm of the prediction error vector — where ¾ = . k is the discrete sample distance, k€ f¾J¥'— I] is an integer, is the number of samples in the echo curve, Δχ is the sample distance, and p is the sampled puise model echo. Step 104 comprises the fine search calculating at least one minimum or maximum using the objective function K¾3 in the vicinity of .¾' , wherein the minimum or maximum %* corresponds to a level of the material or an interface involving the material in sub-sample resolution where the sample resolution is given by ¾.,„ where the center of the reflections, not peaks, are identified. In case there are a set (2 or more) reflection candidates, step 105 can be added comprising using a selection criteria (such as the minimum value of the objective function) to select the best candidate in the set of candidates.
[0019] Step 101 to 104 and optionally step 105 are evaluated for all echo reflections of interest. This includes but is not limited to a reference plane (flange), upper product surface, lower product surface, and the end of probe (or waveguide). The distance to the upper product surface ^^g^^-eS is calculated as the difference between the best reflection candidate for the upper product and the reference plane: simi larly for lower product and end of
probe.
For example, for the coarse search:
A parameterized model = or a measured nearly-perfectly- reflected ultrasound or radar pulse (e.g., transmitted pulse measured after passing through only the receive circuitry of the gauge) is obtained. Here g is a gain or amplitude parameter, w is a width parameter, = m&¾s is the discrete sample distance, m€ f— ^ ] is an integer, and >&¾ is the sample distance. An example of a parametrized model radar pulse reflection versus the actual measured radar pulse reflection is shown in Error! Reference source not found.. The gain or amplitude § (shown in FIG. 2 as the gain) is governed by the amplitude/gain of the transmitted pulse and the dielectric characteristics of the materia! in the tank that the reflected pulse travels through. A time/distance pulse model " S¾
is formed that adjusts an amplitude of the modeled pulse echo signal using a linear attenuation parameter a, where a is a function of the waveguide and the material in the tank in which the pulse is travelling. F&* smm models the radar or ultrasound pulse energy that is lost to its surroundings or is lost due to being absorbed by the waveguide. Reflection candidates (for the center of the reflections) are found in the measured echo curve by evaluating a performance index /{¾.} that is a norm of the prediction error vector For an example with the Euclidean norm, see in eq.
Δχ is the sample distance, M is the half window width of the reflection model, ¾ is the set of discrete distances for the reflection candidates, & is the performance index threshold, and N is number of samples in the echo curve. It is noted that norms other than the Euclidean norm can be used in this step.
5. J¾¾-3 and %| are optionally down sampled or decimated for efficiency in the coarse search.
6. The objective function ¾3 at is then subject to a fine search.
[0Θ20] In the fine search an interpolation technique is used to find a continuous representation of ¾¾ which can then be minimized analytically. An example fine search can include:
1. Forming a 2nd order polynomial representation (curve fit) of ii^ -h Jci $ = ¾«£ % in the vicinity of
3. Calculating the fine search reflection distance candidates —
€ f , which can be solved analytically and = ^ The fine search candidates are the reflection candidates of interest at sub-sample accuracy.
It is noted that other interpolation and curve fitting techniques can be used, such as higher order polynomial or non-polynomial functions.
[0021] In case the set of final reflection distance candidates x* contains more than one candidate an additional candidate selection criterion can be applied as described above relative to step 105. The selection criterion may be different depending on the reflection in question, such as whether reflection in question is from the reference plane, upper product
surface, lower product surface (interface), or the end of probe. In its simplest form, the selection criterion corresponds to the minimum value of the objective function:
Other possible selection criteria can be used, such as using prior knowledge of where the level was in previous time intervals to make the best decision. Instead of the lowest minimum, one can select the first minimum below threshold as it is likely a level as opposed to a second reflection. One can also reject reflections at distance/time integer multiples of the first reflection as they are likely false.
[0022] FIG. 3 depicts an example pulsed radar or ultrasound level finding system shown as a GW system 300 that includes a pulsed radar or ultrasound level gauge shown as a pulse radar level gauge (PRO) 360 including a disclosed level finding algorithm 310a is shown implemented in the firmware of memory 310 associated with a processor 315, according to an example embodiment. Also shown is a transceiver 320 and coaxial connector 325 that is on the top of the tank 305. The transmitter and the receiver provided by transceiver 320 may be implemented as separate blocks. Accordingly, a transceiver as used herein includes both of these arrangements.
[0023] In the transmit direction the processor 315 provides digital signal levels to a digital to analog converter (DAC) 322 which is connected to an input of a pulse generator (PGen) block 321 that is coupled to the transmitter of the transceiver 320. In the receive direction, the receiver of the transceiver 220 receives reflected echo signals that are transduced by the sensor 326, where the output signal from the sensor 326 is coupled to an analog-to-digital converter (ADC) 327 which converts analog signals from the sensor 326 into digital signals for the processor 315 to process which functions as a signal analyzer. Many pulses are generated and the ADC 327 generates the waveform comprising many discrete voltage values that are synchronized by a clock circuit.
[0024] A flange (not shown) may also be present. The waveguide (or probe) is shown as 340. As noted above, although generally described for GWR applications, disclosed level finding can also be applied to ultrasound and non-contacting radar.
EXAMPLES
[0025] Disclosed embodiments are further illustrated by the following specific
Examples, which should not be construed as limiting the scope or content of this Disclosure in any way.
[0026] FIG, 4 illustrates results from an example parametric model and an actual measured echo curve. This example also shows the found reference reflection and the surface (upper product) reflection for a single liquid application, the linear attenuation ·#'~Κ8¾
, and how the model gain g has a step change at the surface boundary due to the difference in dielectric constant between the vapor and the measured material .
[0027] FIG. 5A shows the upper product surface reflection model overlaid on the measured echo curve at the location of the upper product surface. FIG. 5B shows the corresponding objective function and its threshold.
[0028] While various disclosed embodiments have been described above, it should be understood that they have been presented by way of example only, and not limitation. Numerous changes to the subject matter disclosed herein can be made in accordance with this Disclosure without departing from the spirit or scope of this Disclosure. In addition, while a particular feature may have been disclosed with respect to only one of several implementations, such feature may be combined with one or more other features of the other implementations as may be desired and advantageous for any given or particular application.
Claims
1. A method (100) of level finding, comprising: providing (101) parameters including characteristics of a shape of a transmitted pulse in time domain launched by a transmitter onto a waveguide (340) into a tank (305) having at least one material therein, physical properties of said waveguide (340) and real (phase) and imaginar (attenuation) dielectric characterist cs of said material at a frequency of said transmitted pulse; using (102) a processor (315) having said parameters supplied and implementing a level finding algorithm (310a) having a coarse search and a fine search stored in an associated memory (310):
said coarse search minimizing (103) a prediction error between a reflection signal (echo curve y(k)) and a sampled pulse model echo p(k) to obtain an objective function J(k) in a vicinity of a minimum prediction error (k*), and
said fine search calculating (104) at least one minimum or maximum using said J(k) in said vicinity of said k*, wherein said minimum or said maximum corresponds to a level of said material or an interface involving said material.
2. The method of claim 1, wherein said transmitted pulse is a modeled pulse and said y(k) is a modeled echo signal.
3. The method of claim 1, wherein said waveguide (340) contacts said material, and wherein said method comprise guided wave radar (GWR).
4, The method of claim 1 , wherein said fine search calculates a plurality of said minimums or maximums, further comprising selecting a best one from said plurality of said
minimums or maximums using a selection criteria and using said best one for said determining said level of said material or said interface involving said material.
5. A computer program product, comprising:
a non-transitory data storage medium that includes program instructions executable by a processor (315) to enable said processor to execute a method of level finding having a coarse search and a fine search, said computer program product including: code for accessing parameters including characteristics of a shape of a transmitted pulse in time domain launched by a transmitter onto a waveguide (340) into a tank (305 ) having at least one material therein, physical properties of said waveguide and real (phase) and imaginary (attenuation) dielectric characteristics of said material at a frequency of said transmitted pulse; code for said coarse search including minimizing a prediction error between a reflection signal (echo curve y(k)) and a sampled pulse model echo p(k) to obtain an objective function J(k) in a vicinity of a minimum prediction error (k*), and
code for said fine search after said coarse search calculating at least one minimum or maximum using said J(k) in said vicinity of said k*, wherein said minimum or said maximum corresponds to a level of said material or an interface involving said material.
6. The computer program product of claim 5, wherein said transmitted pulse is a modeled pulse and said is a modeled echo signal.
7. The computer program product of claim 5, wherein said waveguide (340) contacts said material, and wherein said method comprises guided wave radar (GWR).
8. The computer program product of claim 5, wherein said fine search calculates a plurality of said minimums or maximums, further comprising selecting a best one from said plurality of said minimums or maximums using a selection criteria and using said best one for said determining said level of said material or said interface involving said material.
9. A pulsed radar or ultrasound level gauge (level gauge) (360), comprising: a processor (315) having an associated memory (310) storing a level finding algorithm (310a) having a coarse search and a fine search; a pulse generator (PGen) block (321) that is coupled to receive control signals originating from said processor; a transceiver (320) for coupling to a waveguide (340) in a tank (305) having at least one material therein having an input coupled to an output of said PGen block (321 ) for transmitting transmitted pulses and an output coupled through a sensor (326) to an input of said processor for processing echo signals received responsive to said transmitted pulses; said level finding algorithm when implemented by said processor having provided parameters including characteristics of a shape of said transmitted pulses in time domain, physical properties of said waveguide and real (phase) and imaginary (attenuation) dielectric characteristics of said material at a frequency of said transmitted pulses causing: said coarse search minimizing a prediction error between said echo signal (echo curve y(k)) and a sampled pulse model echo p(k) to obtain an objective function J(k) in a vicinity of a minimum prediction error (k*), and
said fine search calculating at least one minimum or maximum using said J(k) in said vicinity of said k*, wherein said minimum or said maximum corresponds to a level of said material or an interface involving said material.
10. The level gauge of claim 9, wherein said transmitted pulse is a modeled pulse and said y(k) is a modeled echo signal,
1 1 . The level gauge of claim 9, wherein said transmitted pulse comprises a radar pulse, wherein said waveguide (340) is for contacting said material, and wherein said level gauge comprises guided wave radar (GWR) level gauge.
12. The level gauge of claim 9, wherein said fine search calculates a plurality of said minimums or maximums, further comprising selecting a best one from said plurality of said minimums or maximums using a selection criteria and using said best one for said determining said level of said material or said interface involving said material.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201580053609.8A CN106716081B (en) | 2014-10-01 | 2015-09-29 | It is found using the level of multiple search steps |
EP15847708.3A EP3201580A4 (en) | 2014-10-01 | 2015-09-29 | Level finding using multiple search steps |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201462058358P | 2014-10-01 | 2014-10-01 | |
US62/058,358 | 2014-10-01 | ||
US14/859,694 | 2015-09-21 | ||
US14/859,694 US10180342B2 (en) | 2014-10-01 | 2015-09-21 | Level finding using multiple search steps |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2016053978A1 true WO2016053978A1 (en) | 2016-04-07 |
Family
ID=55631335
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2015/052846 WO2016053978A1 (en) | 2014-10-01 | 2015-09-29 | Level finding using multiple search steps |
Country Status (4)
Country | Link |
---|---|
US (1) | US10180342B2 (en) |
EP (1) | EP3201580A4 (en) |
CN (1) | CN106716081B (en) |
WO (1) | WO2016053978A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10309821B2 (en) | 2016-12-07 | 2019-06-04 | Honeywell International Inc. | Sensor for inventory management applications with remote mounting and asymmetric reflection modeling |
US10386180B2 (en) | 2016-11-28 | 2019-08-20 | Honeywell International Inc. | Apparatus and method for measuring thin material thicknesses in inventory management applications |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10802111B2 (en) * | 2017-02-06 | 2020-10-13 | Magnetrol International, Inc. | Through air radar level transmitter with measurement of first moving echo |
US11221406B2 (en) | 2017-05-09 | 2022-01-11 | Honeywell International Inc. | Guided wave radar for consumable particle monitoring |
RU2653578C1 (en) * | 2017-06-08 | 2018-05-11 | Акционерное общество "ЛИМАКО" | Radar-location level gage for measuring volume of bulk product in tanks |
CN107063392B (en) * | 2017-06-21 | 2019-11-29 | 电子科技大学 | A kind of intelligent liquid level measuring appliance based on microwave transmission line and DSP technology |
US11193809B2 (en) * | 2019-04-01 | 2021-12-07 | Abb Schweiz Ag | Expert control systems and methods for level measurement |
CN110609293B (en) * | 2019-09-19 | 2022-05-27 | 深圳奥锐达科技有限公司 | Distance detection system and method based on flight time |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040257269A1 (en) * | 2002-07-19 | 2004-12-23 | Robert Laun | Method and device for determining an expectancy range for a filling level echo and a false echo |
US20050285777A1 (en) * | 2004-06-24 | 2005-12-29 | Hakan Karlsson | Near zone detection in radar level gauge system |
US20100070209A1 (en) * | 2008-09-18 | 2010-03-18 | Enraf B.V. | Method and apparatus for adaptively handling level measurements under unstable conditions |
US8159386B2 (en) * | 2008-12-30 | 2012-04-17 | Endress + Hauser Gmbh + Co. Kg | Method for ascertaining and monitoring fill level of a medium in a container by a travel time measuring method |
US20140104098A1 (en) * | 2012-10-16 | 2014-04-17 | Magnetrol International, Incorporated | Guided wave radar interface measurement medium identification |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB9211086D0 (en) * | 1992-05-23 | 1992-07-15 | Cambridge Consultants | Short range electromagnetic sensing signal processing |
US5315562A (en) * | 1992-10-23 | 1994-05-24 | Rowe, Deines Instruments Inc. | Correlation sonar system |
US5430445A (en) * | 1992-12-31 | 1995-07-04 | Raytheon Company | Synthetic aperture radar guidance system and method of operating same |
US5424742A (en) * | 1992-12-31 | 1995-06-13 | Raytheon Company | Synthetic aperture radar guidance system and method of operating same |
US5486833A (en) * | 1993-04-02 | 1996-01-23 | Barrett; Terence W. | Active signalling systems |
US6504793B2 (en) * | 2000-09-11 | 2003-01-07 | Vega Grieshaber Kg | Method and device for range measurement |
CA2350816A1 (en) * | 2001-06-14 | 2002-12-14 | Michael L. Cowan | Dynamic sound scattering |
JP4095960B2 (en) * | 2001-07-27 | 2008-06-04 | エンドレス ウント ハウザー ゲーエムベーハー ウント コンパニー コマンディートゲゼルシャフト | Method for measuring filling level of medium in container using filling level measuring device operating based on propagation time principle |
JP3988571B2 (en) * | 2001-09-17 | 2007-10-10 | 株式会社デンソー | Radar equipment |
DE102006019191A1 (en) * | 2006-04-21 | 2007-10-25 | Endress + Hauser Gmbh + Co. Kg | Method for detecting and monitoring the level of a medium in a container |
NL1034327C2 (en) * | 2007-09-04 | 2009-03-05 | Enraf Bv | Method and device for determining the level L of a liquid within a certain measuring range with the aid of radar signals radiated to the liquid level and radar signals reflected by the liquid level. |
DE102007060579B4 (en) * | 2007-12-13 | 2019-04-25 | Endress+Hauser SE+Co. KG | Method for determining and / or assessing the filling state of a container filled with at least one medium |
FR2925169B1 (en) * | 2007-12-18 | 2012-03-09 | Commissariat Energie Atomique | METHOD AND SYSTEM FOR ASSISTING ENVIRONMENTAL CHARACTERIZATION WITH ULTRA BROADBAND RADIO FREQUENCY SIGNALS. |
RU2410650C2 (en) * | 2008-11-01 | 2011-01-27 | ООО "Предприятие "Контакт-1" | Method to measure level of material in reservoir |
EP2442129B1 (en) * | 2010-10-18 | 2016-03-23 | Siemens Aktiengesellschaft | A method for processing an echo amplitude profile generated by a pulse-echo ranging system |
BR112014001545A2 (en) * | 2011-08-04 | 2017-06-13 | Vega Grieshaber Kg | linear relationship between tracks |
-
2015
- 2015-09-21 US US14/859,694 patent/US10180342B2/en active Active
- 2015-09-29 CN CN201580053609.8A patent/CN106716081B/en active Active
- 2015-09-29 WO PCT/US2015/052846 patent/WO2016053978A1/en active Application Filing
- 2015-09-29 EP EP15847708.3A patent/EP3201580A4/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040257269A1 (en) * | 2002-07-19 | 2004-12-23 | Robert Laun | Method and device for determining an expectancy range for a filling level echo and a false echo |
US20050285777A1 (en) * | 2004-06-24 | 2005-12-29 | Hakan Karlsson | Near zone detection in radar level gauge system |
US20100070209A1 (en) * | 2008-09-18 | 2010-03-18 | Enraf B.V. | Method and apparatus for adaptively handling level measurements under unstable conditions |
US8159386B2 (en) * | 2008-12-30 | 2012-04-17 | Endress + Hauser Gmbh + Co. Kg | Method for ascertaining and monitoring fill level of a medium in a container by a travel time measuring method |
US20140104098A1 (en) * | 2012-10-16 | 2014-04-17 | Magnetrol International, Incorporated | Guided wave radar interface measurement medium identification |
Non-Patent Citations (1)
Title |
---|
See also references of EP3201580A4 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10386180B2 (en) | 2016-11-28 | 2019-08-20 | Honeywell International Inc. | Apparatus and method for measuring thin material thicknesses in inventory management applications |
US10309821B2 (en) | 2016-12-07 | 2019-06-04 | Honeywell International Inc. | Sensor for inventory management applications with remote mounting and asymmetric reflection modeling |
Also Published As
Publication number | Publication date |
---|---|
CN106716081A (en) | 2017-05-24 |
US20160097669A1 (en) | 2016-04-07 |
US10180342B2 (en) | 2019-01-15 |
CN106716081B (en) | 2019-11-19 |
EP3201580A1 (en) | 2017-08-09 |
EP3201580A4 (en) | 2018-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10180342B2 (en) | Level finding using multiple search steps | |
RU2678933C2 (en) | Self-diagnosing frequency-modulated continuous-wave radar level gauge | |
US9534944B2 (en) | Method for determining and/or monitoring fill level of a medium in a container | |
US10055519B2 (en) | Pulse shape change for interface determination | |
US9612147B2 (en) | Radar level gauge system with multiple receiver branches | |
US20120265486A1 (en) | Method for ascertaining and monitoring fill level of a medium in a container with a travel time measuring method | |
US9069056B2 (en) | Guided wave radar probe reference target | |
US20120158321A1 (en) | Wireless Liquid Quantity Measurement System | |
US9442000B2 (en) | Method for ascertaining and monitoring fill level of a medium in a container by means of a fill-level measuring device using a travel time measuring method | |
US20090146867A1 (en) | Radar level gauge system and method providing a signal indicative of process reliability | |
JP2008089583A (en) | Radar/level measurement | |
US9618617B2 (en) | Level measurement using correlation between a pair of secondary reference signals | |
US9671488B2 (en) | Radar level gauge with signal division | |
CN103261854B (en) | The determination of the dielectric property when level gauging | |
CN109357729B (en) | Method for measuring liquid interface level based on guided wave radar liquid level meter | |
WO2016054000A1 (en) | Resolution mode switching for pulsed radar | |
CA2286439A1 (en) | Improvements in time domain reflectometry | |
CN117452399A (en) | Marking box obstacles using electronic level | |
US10113901B2 (en) | Method for evaluating a TDR limit level switch | |
WO2019086196A1 (en) | Radar level gauge system and method for interface measurement | |
US10145719B2 (en) | Echo curve determination at a resolution that differs on area-by-area basis | |
US20230228612A1 (en) | Level gauge and method for sensing material levels in tanks | |
CN106989794B (en) | Method for determining the filling level of a product in a tank and radar level gauge system | |
KR20150119623A (en) | Method for oil level changing tracking and detailed dectection using radar sensor | |
RU2413180C1 (en) | Apparatus for measuring thickness of dielectric coating |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 15847708 Country of ref document: EP Kind code of ref document: A1 |
|
REEP | Request for entry into the european phase |
Ref document number: 2015847708 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2015847708 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |