CN104913830A - Radar level instrument signal processing method based on curve fitting - Google Patents

Radar level instrument signal processing method based on curve fitting Download PDF

Info

Publication number
CN104913830A
CN104913830A CN201510252502.8A CN201510252502A CN104913830A CN 104913830 A CN104913830 A CN 104913830A CN 201510252502 A CN201510252502 A CN 201510252502A CN 104913830 A CN104913830 A CN 104913830A
Authority
CN
China
Prior art keywords
rising edge
echo
time
slope
crest
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201510252502.8A
Other languages
Chinese (zh)
Inventor
毕欣
高洁
杜劲松
王伟
赵越南
田星
仝盼盼
张青石
李想
徐洪庆
丛日刚
高扬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang Institute of Automation of CAS
Original Assignee
Shenyang Institute of Automation of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenyang Institute of Automation of CAS filed Critical Shenyang Institute of Automation of CAS
Priority to CN201510252502.8A priority Critical patent/CN104913830A/en
Publication of CN104913830A publication Critical patent/CN104913830A/en
Pending legal-status Critical Current

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to a radar level instrument signal processing method based on curve fitting. The method includes the steps of calculating the slope of each point of an echo rising edge, judging whether the duration time of the rising edge meets the minimum keeping time or not, determining a certain rising edge as a target estimation position through the maximum slope, fitting a parabola through the least square method according to data of a 60% to 90% section of a crest of the rising edge at the target estimation position, calculating the peak value and corresponding time, and finally accurately positioning a level instrument. By means of the method, the purpose of positioning the multi-diameter multi-echo crest is effectively achieved, the influences of the measurement environment and a target object itself on an echo rear edge are eliminated, and the positioning accuracy is improved.

Description

A kind of radar levelmeter signal processing method based on curve
Technical field
The present invention relates to field of industrial measurement, specifically a kind of radar levelmeter signal processing method based on curve.
Background technology
Thing position is industrial important parameter.The level measurement instruments such as existing static pressure method, float-ball type, ultrasound wave, have respective features and application scope.Radar levelmeter, by advanced treatment technology, has unique advantage, demonstrates the performance of its brilliance especially, play more and more important effect in the industrial production under the harsh measurement conditions such as high temperature, high pressure, steam, corrosive medium are strong.
Industry thing position, liquid level gauge operating mode complexity, cause the problem of measuring accuracy under the measuring condition as chaff interference, guided wave bar coating, foam, steam and turbulent flow.Typical algorithm is detection peak, and the prerequisite of algorithm application is that echo waveform is enough good.In actual application, waveform multipath effect often has multiple echo, and such as, wall echo causes the peak value of more than, is exactly in addition badly to determine peak value.
Existing placement technology is by carrying out linear interpolation realization to echo, along the impact being subject to measurement environment and destination object itself after peak value and echo, linear interpolation is difficult to compound actual waveform characteristic, and this makes to be difficult to determine correct echo time of arrival, and precision improvement is limited in one's ability.
Summary of the invention
For the deficiencies in the prior art, curve is fused in radar levelmeter location algorithm by the present invention, proposes a kind of radar levelmeter signal processing method based on curve.
The technical scheme that the present invention is adopted for achieving the above object is: a kind of radar levelmeter signal processing method based on curve, comprises the following steps:
Step 1: using the reflection echo under aimless original state as Noise Background N;
Step 2: compared with Noise Background N by the reflection echo X of actual measurement and disappear mutually, disappeared echo S mutually;
Step 3: to the smoothing filtering of echo S that disappears mutually, with the small amplitude fluctuations of filtering, obtains filtered echo S w, calculate echo S win the slope k often put of each rising edge, determine rising edge trough S wtnto crest S wcninitial time t snwith end time t en, i.e. the duration T of rising edge hn=t en-t snif, rising edge duration T hnbe less than minimum hold time T hmin, then this rising edge is invalid;
Step 4: in the effective situation of rising edge, using the rising edge at the maximum some place of slope as target state estimator position;
Step 5: according to the crest S of target state estimator position wcnthe data in 60% to 90% interval, utilize least square fitting para-curve;
Step 6: obtain peak value by the Fitting Calculation, and the time t of correspondence;
Step 7: R=(t-t accurately located by level meter c) c/2, wherein, R be the thing distance of positions from, c is the light velocity, t cfor the clock initial correction time.
The present invention has following beneficial effect and advantage:
1. the present invention is directed to multipath effect in actual application and produce multiple echo, by the feature determination target echo of rising edge;
2. along the impact being subject to measurement environment and destination object itself after the present invention utilizes least square fitting para-curve to solve peak value and echo;
3. the present invention utilizes least square fitting para-curve to promote the positioning precision of level meter.
Accompanying drawing explanation
Fig. 1 is method flow diagram of the present invention;
Fig. 2 (a) shows radar levelmeter echo situation 1;
Fig. 2 (b) shows radar levelmeter echo situation 2;
Fig. 2 (c) shows radar levelmeter echo situation 3;
Fig. 3 is curve schematic diagram of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail.
As shown in Figure 1, be schematic flow sheet of the present invention.First calculate the slope that echo rising edge is often put, and judge whether the rising edge duration meets minimum hold time; Then determine that certain rising edge is as target state estimator position by maximum slope; Again by the data separate least square fitting para-curve in 60% to 90% interval of the rising edge crest of target state estimator position, calculate the time of peak value and correspondence, finally obtain level meter and accurately locate.Concrete steps are as follows:
Step 1: using the reflection echo under aimless original state as Noise Background N;
Step 2: compared with Noise Background N by the reflection echo X of actual measurement and disappear mutually, disappeared echo S mutually;
Step 3: to the smoothing filtering of echo S that disappears mutually, the small amplitude fluctuations of filtering, obtains filtered echo S w, calculate echo S win the slope k often put of the n-th rising edge, determine rising edge trough S wtnto crest S wcninitial time t snwith end time t en, i.e. the duration T of rising edge hn=t en-t snif, rising edge duration T hnbe less than minimum hold time T hmin, by invalid for this rising edge;
Step 4: maximum slope k maxbe the slope of the n-th rising edge point, then using the n-th rising edge as target state estimator position;
As shown in Fig. 2 (a), there are 3 rising edges in radar levelmeter echo, the duration of 3 rising edges is all greater than minimum hold time, maximum slope k maxbe the slope of the 2nd rising edge point, then using the 2nd rising edge as target state estimator position;
As shown in Fig. 2 (b), there are 3 rising edges in radar levelmeter echo, and the duration of the 3rd rising edge is all less than minimum hold time, then the 3rd rising edge is invalid, maximum slope k maxbe the slope of the 2nd rising edge point, then using the 2nd rising edge as target state estimator position;
As shown in Fig. 2 (c), there are 3 rising edges in radar levelmeter echo, and the duration of the 3rd rising edge is all less than minimum hold time, then the 3rd rising edge is invalid, and the slope of the 3rd rising edge point is maximum slope k max, because the 3rd rising edge is invalid, after excluding invalid rising edge, maximum slope k maxbe the slope of the 2nd rising edge point, then using the 2nd rising edge as target state estimator position;
Step 5: the n-th rising edge crest S wcnthe data separate least square fitting para-curve in 60% to 90% interval;
Step 6: obtain peak value by the Fitting Calculation, and the time t of correspondence;
As shown in Figure 3, radar n-th rising edge point A is rising edge crest S with the value of some B wcn60% and 90%, by the data separate least square fitting para-curve in this interval, obtain dotted line, dotted line fixed point C is the peak value that curve calculates, and the corresponding time is t;
Step 7: R=(t-t accurately located by level meter c) c/2, wherein, R be the thing distance of positions from, c is the light velocity, t cfor the clock initial correction time.

Claims (1)

1., based on a radar levelmeter signal processing method for curve, it is characterized in that: comprise the following steps:
Step 1: using the reflection echo under aimless original state as Noise Background N;
Step 2: compared with Noise Background N by the reflection echo X of actual measurement and disappear mutually, disappeared echo S mutually;
Step 3: to the smoothing filtering of echo S that disappears mutually, with the small amplitude fluctuations of filtering, obtains filtered echo S w, calculate echo S win the slope k often put of each rising edge, determine rising edge trough S wtnto crest S wcninitial time t snwith end time t en, i.e. the duration T of rising edge hn=t en-t snif, rising edge duration T hnbe less than minimum hold time T hmin, then this rising edge is invalid;
Step 4: in the effective situation of rising edge, using the rising edge at the maximum some place of slope as target state estimator position;
Step 5: according to the crest S of target state estimator position wcnthe data in 60% to 90% interval, utilize least square fitting para-curve;
Step 6: obtain peak value by the Fitting Calculation, and the time t of correspondence;
Step 7: R=(t-t accurately located by level meter c) c/2, wherein, R be the thing distance of positions from, c is the light velocity, t cfor the clock initial correction time.
CN201510252502.8A 2015-05-15 2015-05-15 Radar level instrument signal processing method based on curve fitting Pending CN104913830A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510252502.8A CN104913830A (en) 2015-05-15 2015-05-15 Radar level instrument signal processing method based on curve fitting

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510252502.8A CN104913830A (en) 2015-05-15 2015-05-15 Radar level instrument signal processing method based on curve fitting

Publications (1)

Publication Number Publication Date
CN104913830A true CN104913830A (en) 2015-09-16

Family

ID=54083091

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510252502.8A Pending CN104913830A (en) 2015-05-15 2015-05-15 Radar level instrument signal processing method based on curve fitting

Country Status (1)

Country Link
CN (1) CN104913830A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106443646A (en) * 2016-11-21 2017-02-22 重庆兆洲科技发展有限公司 Ultrasonic ranging system and echo processing method and device
CN107505027A (en) * 2016-06-14 2017-12-22 Vega格里沙贝两合公司 Microwave reflection barrier
CN108426520A (en) * 2018-05-15 2018-08-21 天津大学 The detection device in pipeline stolen hole and the detection method at stolen hole center, diameter
CN110596662A (en) * 2019-10-21 2019-12-20 富临精工先进传感器科技(成都)有限责任公司 Distance deviation correction method of MIMO radar
CN111812606A (en) * 2020-06-03 2020-10-23 西安电子科技大学 Guided wave radar-based level extraction method
CN113296066A (en) * 2021-04-27 2021-08-24 北京锐达仪表有限公司 Echo correction method, device, readable storage medium, electronic device and system
CN113639816A (en) * 2021-07-06 2021-11-12 燕山大学 Non-invasive method for monitoring the level of a liquid in a container and system for implementing said method
CN114637021A (en) * 2022-05-18 2022-06-17 四川吉埃智能科技有限公司 Sub-centimeter-level full-waveform laser radar ranging method and device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5969666A (en) * 1997-06-06 1999-10-19 Endress + Hauser GmbH +Co. Radar-based method of measuring the level of a material in a container
US20070201309A1 (en) * 2006-02-21 2007-08-30 Robert Brown Method for processing an echo profile and a pulse-echo ranging system
CN102661773A (en) * 2012-05-15 2012-09-12 合肥工业大学 System and method for processing echo signal of guided wave type radar level meter
CN203177928U (en) * 2012-09-21 2013-09-04 清华大学 Echo sampling unit and pulse type radar material level meter with the same
CN103542911A (en) * 2013-11-01 2014-01-29 合肥工业大学 Guided wave radar level meter echo signal processing system and method based on first-order derivative

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5969666A (en) * 1997-06-06 1999-10-19 Endress + Hauser GmbH +Co. Radar-based method of measuring the level of a material in a container
US20070201309A1 (en) * 2006-02-21 2007-08-30 Robert Brown Method for processing an echo profile and a pulse-echo ranging system
CN102661773A (en) * 2012-05-15 2012-09-12 合肥工业大学 System and method for processing echo signal of guided wave type radar level meter
CN203177928U (en) * 2012-09-21 2013-09-04 清华大学 Echo sampling unit and pulse type radar material level meter with the same
CN103542911A (en) * 2013-11-01 2014-01-29 合肥工业大学 Guided wave radar level meter echo signal processing system and method based on first-order derivative

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107505027A (en) * 2016-06-14 2017-12-22 Vega格里沙贝两合公司 Microwave reflection barrier
US10704948B2 (en) 2016-06-14 2020-07-07 Vega Grieshaber Kg Reflection microwave barrier
CN106443646A (en) * 2016-11-21 2017-02-22 重庆兆洲科技发展有限公司 Ultrasonic ranging system and echo processing method and device
CN106443646B (en) * 2016-11-21 2019-03-22 重庆兆洲科技发展有限公司 A kind of ultrasonic ranging system, echo processing techniques and device
CN108426520A (en) * 2018-05-15 2018-08-21 天津大学 The detection device in pipeline stolen hole and the detection method at stolen hole center, diameter
CN110596662A (en) * 2019-10-21 2019-12-20 富临精工先进传感器科技(成都)有限责任公司 Distance deviation correction method of MIMO radar
CN111812606A (en) * 2020-06-03 2020-10-23 西安电子科技大学 Guided wave radar-based level extraction method
CN111812606B (en) * 2020-06-03 2023-12-22 西安电子科技大学 Material level extraction method based on guided wave radar
CN113296066A (en) * 2021-04-27 2021-08-24 北京锐达仪表有限公司 Echo correction method, device, readable storage medium, electronic device and system
CN113296066B (en) * 2021-04-27 2023-06-23 北京锐达仪表有限公司 Echo correction method, device, readable storage medium, electronic equipment and system
CN113639816A (en) * 2021-07-06 2021-11-12 燕山大学 Non-invasive method for monitoring the level of a liquid in a container and system for implementing said method
CN114637021A (en) * 2022-05-18 2022-06-17 四川吉埃智能科技有限公司 Sub-centimeter-level full-waveform laser radar ranging method and device

Similar Documents

Publication Publication Date Title
CN104913830A (en) Radar level instrument signal processing method based on curve fitting
CN106643939B (en) The method for calculating ultrasonic propagation time for ultrasonic flowmeter
CN100565241C (en) Ultrasonic echo frontier inspection based on modulation domain measurement
CN107860430B (en) Time difference measurement method of ultrasonic gas flowmeter based on time difference method
CN104501889B (en) Ultrasonic flow detection method based on cross-correlation time difference method
CN102297712A (en) Method for measuring propagation time of ultrasonic echo
CN110297250B (en) Initial grazing angle solving method based on Taylor expansion, sound ray bending correction method and equipment
CN106706056B (en) A kind of compensation method of heavy caliber ultrasonic water meter flow measurement
RU2012133315A (en) VARIABILITY TRACKING
CN105180853A (en) Electromagnetic ultrasonic metal material thickness measurement method
CN104515560A (en) Ultrasonic flow meter and dirt detection and error compensation method thereof
CN104748704A (en) Thin-walled structure ultrasonic resonance thickness measurement frequency spectrum analysis interpolation correction method
US11346694B2 (en) Method for measuring the velocity of a fluid using ultrasound
KR101764870B1 (en) Signal processing system for ultrasonic floemeter
CN109470320B (en) Step filtering method for ultrasonic water meter
CN112304376A (en) Ultrasonic flowmeter flow measuring method based on data fusion
CN105785123A (en) Radar signal frequency calculation method based on apFFT phase difference
KR101386641B1 (en) Distance estimation device and method using the difference of wave speed between waves
CN114355360B (en) Echo frequency estimation method of Doppler log
CN104502627B (en) Method is solved based on transmitting Design of Signal and the phase ambiguity of processing in ADCP
CN110346783B (en) Sound line correction method based on power series equivalent ranging model
JP4392463B2 (en) Film thickness measurement method
JP4904099B2 (en) Pulse signal propagation time measurement device and ultrasonic flow measurement device
JP7174574B2 (en) Ultrasonic flowmeter and determination method of zero-crossing time in ultrasonic flowmeter
CN109031260B (en) LFM signal time delay measurement method based on fractional Fourier modulation rate analysis

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20150916