CN106842145A - 一种抑制fmcw雷达液位计回波冲击响应干扰的方法 - Google Patents

一种抑制fmcw雷达液位计回波冲击响应干扰的方法 Download PDF

Info

Publication number
CN106842145A
CN106842145A CN201510880669.9A CN201510880669A CN106842145A CN 106842145 A CN106842145 A CN 106842145A CN 201510880669 A CN201510880669 A CN 201510880669A CN 106842145 A CN106842145 A CN 106842145A
Authority
CN
China
Prior art keywords
detection
point
shock response
formula
level gauge
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
CN201510880669.9A
Other languages
English (en)
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.)
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Original Assignee
Leihua Electronic Technology Research Institute Aviation Industry Corp of China
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 Leihua Electronic Technology Research Institute Aviation Industry Corp of China filed Critical Leihua Electronic Technology Research Institute Aviation Industry Corp of China
Priority to CN201510880669.9A priority Critical patent/CN106842145A/zh
Publication of CN106842145A publication Critical patent/CN106842145A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01FMEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00Indicating 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/22Indicating 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/28Indicating 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/284Electromagnetic waves

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Thermal Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种抑制FMCW雷达液位计回波冲击响应干扰的方法,其特征在于,包括以下步骤:冲击响应干扰区的检测:将实测的检测采样点x0,x1,…,xN-1进行解析;在时域内,检测满足式(1)的检测采样点其中, 式中,xi,(i=1,2,…,N-1)为第i个采样点;将在检测采样点中的幅度最大值点左右各取15个点作为冲击响应干扰区,若检测采样点左侧不足15点,则左侧取从第一个点算起至幅度最大值点之间所有采样点;在时域原始信号中,挖除冲击干扰区的检测采样点,补零;基于AR模型预测内插。

Description

一种抑制FMCW雷达液位计回波冲击响应干扰的方法
技术领域
雷达液位计测距技术
背景技术
在国外,FMCW雷达液位计被广泛应用于石化、仓储行业,是储油罐的核心传感器,FMCW采用的是连续的信号发射及处理而非脉冲雷达的信号发射—等待—接收的间歇式工作方式,也就更适合于过程罐、搅拌罐等液面剧烈的场合,为中高端应用场合使用(罗斯蒙特雷达液位股份公司.雷达液位测量.发明专利,专利申请号:CN200710153641.0)。
但在雷达液位计研发过程中,由于雷达接收机接收通道的瞬态冲击响应,导致对雷达回波带来了冲击响应干扰,这种干扰主要为冲击噪声,这种干扰远强于油罐液位回波,这种干扰使得目标的SNR急剧下降从而对目标检测带来误判,进而影响雷达液位计的测量精度,在cm级别,不满足测位精度mm级别的量级。
因此为了满足以上需求,本发明提出了一种抑制FMCW雷达液位计回波冲击响应干扰的算法,并经过实测数据验证时有效的,并已应用到中航工业雷华电子技术研究所自主研发高精度雷达液位计上(中国航空工业集团公司雷华电子技术研究所.一种雷达液位计.发明专利,专利申请号:201430463897.2),得出的SNR提高了很多,得出精度为1mm量级。
发明内容
本发明的目的是在国产高精度雷达液位计研发需求,在雷达液位计锯齿波回波存在冲击响应干扰影响液位计测量精度和误测情况下,且冲击响应干扰区与检测目标区为不同区域,提出了一种在时域检测出冲击响应干扰的位置后,在原始信号中挖除存在冲击响应干扰的回波点,基于AR模型预测内插(AR模型系数采用Burg算法估计)算法,消除冲击响应对目标掩盖作用,从而保证液位检测出来足够的SNR。
本发明在研发的FMCW雷达液位计实测数据上进行了定量验证,算法是有效的,SNR提高25dB以上,从而保证雷达液位计高精度测量1mm奠定基础。
实现本发明目的的技术解决方案为:本发明提供了1.一种抑制FMCW雷达液位计回波冲击响应干扰的方法,其特征在于,包括以下步骤:
1.冲击响应干扰区的检测
a)将实测的检测采样点x0,x1,…,xN-1进行解析;
b)在时域内,检测满足式(1)的检测采样点
其中,
式中,xi,(i=1,2,…,N-1)为第i个采样点;
c)将在检测采样点中的幅度最大值点左右各取15个点作为冲击响应干扰区,若检测采样点左侧不足15点,则左侧取从第一个点算起至幅度最大值点之间所有采样点;
2.在时域原始信号中,挖除冲击干扰区的检测采样点,补零。
3.基于AR模型预测内插
采用AR模型预测内插对挖除冲击干扰区的检测采样点的原始回波信号进行恢复,AR模型系数采用Burg算法估计。
进一步的,
Burg算法是通过前N点的检测采样点序列x0,x1,…,xN-1,估计AR模型参数ap,m,(m=1,2,…,p),预测第N个检测采样点xN,有下列关系:
式中p是AR模型的阶数。由莱文生-杜平递推公式计算AR参数
ap,m=ap-1,m+ap,pap-1,p-m (4)
式中ap,p由下式给出
式中,kp称为反射系数,是预测误差的均方误差:
x(n-1)为采样点序列,通过式(3)~(7)可以确定AR模型的一组参数{ap,m},其中AR模型的阶数预先给定。
发明创造的优点和用途
(1)抑制冲击响应干扰效果好。对于FMCW雷达液位计采用锯齿波模式下,存在冲击响应干扰会将目标的噪底大大提升,目标的SNR急剧减少,难以满足储油罐液位1mm精度,本发明提出的方法可以有效抑制冲击响应干扰,恢复目标的SNR,基于本发明方法可使得SNR提升25dB以上,使得最终SNR为40dB及以上。
(2)条件约束少,工程易于实现,因为本发明提出的是在时域检测干扰,时域挖空干扰补零,并基于AR模型预测内插,不需要消耗过多的资源,且较容易实现;
(3)测量精度高,基于本发明提出抑制冲击响应干扰后,FMCW雷达液位计测量精度为1mm,这与去除干扰前测量精度为cm量级,提升了一个数量级;
(4)算法稳定性好。通过雷达液位计在某油罐外场试验基地连续30测试,雷达液位计测试结果良好,雷达液位计的稳定性<=1mm(依据JJG971-2002 7.2.3.5),雷达液位计的重复性指标为0.32mm<0.4mm(依据GB/T 13983-1992 3.35、GB/T 21117-2007 6.1.4的重复性规定:液位计的重复性应不超过最大允许示值误差绝对值的五分之二,最大允许误差绝对值为1mm),反映了雷达液位计具有较好的稳定性和重复性指标。
(5)代价小,运算量小。
附图说明
图1是本发明给出的雷达液位计测量储油罐液面场景示意图;
图2是本发明给出的抑制FMCW雷达液位计冲击响应干扰的算法流程图;
图3是本发明给出的是含有冲击响应干扰的雷达液位计时域回波图;
图4是本发明给出的是含有冲击响应干扰的雷达液位计回波(图3)的频域图;
图5是基于本发明提出的算法对图3所示的含有冲击响应干扰进行抑制效果的时域图;
图6是基于本发明提出的算法对图3所示的含有冲击响应干扰进行抑制效果的频域图;
图7是基于本发明抑制冲击响应干扰后的雷达液位计连续31小时实测油罐液位图。
具体实施方式
图1是本发明给出的是雷达液位计工作场景示意图,图2是本本发明给出的抑制FMCW雷达液位计冲击响应干扰的算法流程图。结合图2,于AR模型预测内插的算法来抑制FMCW雷达液位计的冲击响应干扰包括击响应干扰区的检测、在原始信号中挖除存在冲击响应干扰的回波点、基于AR模型预测内插三部分构成。
如图3所示,冲击响应干扰在雷达液位计回波前段非常明显,图4是本发明给出的是含有冲击响应干扰的雷达液位计回波(图3)的频域图,可以看出干扰的存在使得目标的SNR急剧下降,噪底抬升,需要抑制干扰。
本发明提供了一种抑制FMCW雷达液位计回波冲击响应干扰的方法,其特征在于,包括以下步骤:
1.冲击响应干扰区的检测
(1)将实测的检测采样点x0,x1,…,xN-1进行解析;
(2)在时域内,检测满足式(1)的检测采样点
其中,
式中,xi,(i=1,2,…,N-1)为第i个采样点;
(3)将在检测采样点中的幅度最大值点左右各取15个点作为冲击响应干扰区,若检测采样点左侧不足15点,则左侧取从第一个点算起至幅度最大值点之间所有采样点;
2.在时域原始信号中,挖除冲击干扰区的检测采样点,补零。
3.基于AR模型预测内插
采用AR模型预测内插对挖除冲击干扰区的检测采样点的原始回波信号进行恢复,AR模型系数采用Burg算法估计。
图5给出的是基于本发明提出的算法对图3所示的含有冲击响应干扰进行抑制效果的时域图;由图5可知,雷达液位计的时域回波图已经恢复接近正常的三角波形。
图6是基于本发明提出的算法对图3所示的含有冲击响应干扰进行抑制效果的频域图;由图6可知,抑制干扰后目标的SNR为60dB,而由图4给出的抑制干扰前的目标的SNR为30dB,可见此算法大大提高了抑制干扰效果。
抑制干扰效果评判
图7给出基于本发明提出的算法抑制干扰后,雷达液位计在某油罐外场试验基地的油罐(高8m*直径7m)连续测量31小时,稳定性和重复性能较好,图7给出了统计的稳定性和重复性试验结果,统计10000个数据点,雷达液位计的稳定性<=1mm(依据JJG971-20027.2.3.5),雷达液位计的重复性指标为0.32mm<0.4mm(依据GB/T 13983-1992 3.35、GB/T21117-2007 6.1.4的重复性规定:液位计的重复性应不超过最大允许示值误差绝对值的五分之二,最大允许误差绝对值为1mm),反映了雷达液位计具有较好的稳定性和重复性指标。
同时基于本发明抑制冲击响应干扰后的雷达液位计在总后油研所1号油罐测量数据与基于Honeywell的标尺HERMetic UTImeter(具有国际计量认证,精度为0.5mm)进行比对,如图7所示,基于雷达液位计实测数据进行数理统计得出雷达液位计精度为1mm(依据JJG971-2002 5.2、GB/T 25964-2010),满足设计要求,这进一步验证了本发明在工程应用的有效性。

Claims (2)

1.一种抑制FMCW雷达液位计回波冲击响应干扰的方法,其特征在于,包括以下步骤:
1.冲击响应干扰区的检测
a)将实测的检测采样点x0,x1,…,xN-1进行解析;
b)在时域内,检测满足式(1)的检测采样点
| 20 * log 10 ( x i ) - 20 * log 10 ( x &OverBar; ) | > 3 - - - ( 1 )
其中, x &OverBar; = 1 N &Sigma; i = 0 N - 1 x i - - - ( 2 )
式中,xi,(i=1,2,…,N-1)为第i个采样点;
c)将在检测采样点中的幅度最大值点左右各取15个点作为冲击响应干扰区,若检测采样点左侧不足15点,则左侧取从第一个点算起至幅度最大值点之间所有采样点;
2.在时域原始信号中,挖除冲击干扰区的检测采样点,补零,
3.基于AR模型预测内插
采用AR模型预测内插对挖除冲击干扰区的检测采样点的原始回波信号进行恢复,AR模型系数采用Burg算法估计。
2.根据权利要求1所述的抑制FMCW雷达液位计回波冲击响应干扰的方法,其特征在于,
Burg算法是通过前N点的检测采样点序列x0,x1,…,xN-1,估计AR模型参数ap,m,(m=1,2,…,p),预测第N个检测采样点xN,有下列关系:
x ^ N = - &Sigma; m = 1 p a p , m x N - m - - - ( 3 )
式中p是AR模型的阶数。由莱文生-杜平递推公式计算AR参数
ap,m=ap-1,m+ap,pap-1,p-m (4)
式中ap,p由下式给出
a p , p = k p = - r x x ( p ) + &Sigma; k = 1 p - 1 a p - 1 r x x ( p - k ) &sigma; p - 1 2 - - - ( 5 )
式中,kp称为反射系数, &sigma; p - 1 2 是预测误差的均方误差:
r x x ( p ) = &sigma; p 2 = ( 1 - k p 2 ) &sigma; p - 1 2 - - - ( 6 )
&sigma; 0 2 = r x x ( 0 ) = E &lsqb; ( x ( n - 1 ) ) 2 &rsqb; - - - ( 7 )
x(n-1)为采样点序列,通过式(3)~(7)可以确定AR模型的一组参数{ap,m},其中AR模型的阶数预先给定。
CN201510880669.9A 2015-12-03 2015-12-03 一种抑制fmcw雷达液位计回波冲击响应干扰的方法 Pending CN106842145A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510880669.9A CN106842145A (zh) 2015-12-03 2015-12-03 一种抑制fmcw雷达液位计回波冲击响应干扰的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510880669.9A CN106842145A (zh) 2015-12-03 2015-12-03 一种抑制fmcw雷达液位计回波冲击响应干扰的方法

Publications (1)

Publication Number Publication Date
CN106842145A true CN106842145A (zh) 2017-06-13

Family

ID=59149761

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510880669.9A Pending CN106842145A (zh) 2015-12-03 2015-12-03 一种抑制fmcw雷达液位计回波冲击响应干扰的方法

Country Status (1)

Country Link
CN (1) CN106842145A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113049061A (zh) * 2021-03-04 2021-06-29 国家石油天然气管网集团有限公司华南分公司 一种油罐液位计高精度距离估计方法
CN116500553A (zh) * 2023-06-07 2023-07-28 北京交通大学 雷达互相干扰抑制方法、装置、计算机设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005078475A1 (en) * 2004-02-13 2005-08-25 Saab Rosemount Tank Radar Ab A method and an arrangement in a radar level gauging system
CN101140180A (zh) * 2006-09-07 2008-03-12 罗斯蒙特雷达液位股份公司 雷达液位测量
CN101907703A (zh) * 2010-07-20 2010-12-08 中国人民解放军空军雷达学院 天波超视距雷达空时联合自适应抗干扰方法
US20150035695A1 (en) * 2013-08-05 2015-02-05 Finetek Co., Ltd. Frequency modulation continuous wave radar level meter and signal-tracking and phase-locking method for the same

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005078475A1 (en) * 2004-02-13 2005-08-25 Saab Rosemount Tank Radar Ab A method and an arrangement in a radar level gauging system
CN101140180A (zh) * 2006-09-07 2008-03-12 罗斯蒙特雷达液位股份公司 雷达液位测量
CN101907703A (zh) * 2010-07-20 2010-12-08 中国人民解放军空军雷达学院 天波超视距雷达空时联合自适应抗干扰方法
US20150035695A1 (en) * 2013-08-05 2015-02-05 Finetek Co., Ltd. Frequency modulation continuous wave radar level meter and signal-tracking and phase-locking method for the same

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周忠根: "天波超视距雷达瞬态干扰抑制和目标检测方法", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113049061A (zh) * 2021-03-04 2021-06-29 国家石油天然气管网集团有限公司华南分公司 一种油罐液位计高精度距离估计方法
CN113049061B (zh) * 2021-03-04 2024-02-13 国家石油天然气管网集团有限公司华南分公司 一种油罐液位计高精度距离估计方法
CN116500553A (zh) * 2023-06-07 2023-07-28 北京交通大学 雷达互相干扰抑制方法、装置、计算机设备及存储介质
CN116500553B (zh) * 2023-06-07 2023-10-03 北京交通大学 雷达互相干扰抑制方法、装置、计算机设备及存储介质

Similar Documents

Publication Publication Date Title
US10310056B2 (en) Apparatus and method for adjusting guided wave radar pulse width to optimize measurements
US8510065B2 (en) Wireless liquid quantity measurement system
US9557208B2 (en) Liquid level measuring apparatus, method, and program
CN102798434B (zh) 确定容器中边界表面的位置的特征变量的评估设备和方法
WO2014192326A1 (ja) 津波監視システム
CN102192772B (zh) 填充水平测量设备中的移动性检测
CN103675922A (zh) 基于探地雷达的运营期地下管道管径测定方法
US20140109677A1 (en) Ultrasonic measurement apparatus and method
CN102621529A (zh) 一种脉冲回波测距系统中的回波处理方法
WO2016053978A1 (en) Level finding using multiple search steps
CN102928828B (zh) 一种基于正交波形的分布式雷达的相位差估计方法
CN104913830A (zh) 一种基于曲线拟合的雷达物位计信号处理方法
CN106842145A (zh) 一种抑制fmcw雷达液位计回波冲击响应干扰的方法
Holland Seabed reflection measurement uncertainty
SA518391069B1 (ar) فحص صدى النبض المدمج بأنظمة خطوط الأنابيب
CN102546052B (zh) 一种自适应门限控制信号参数统计测量方法
CN103328939A (zh) 基于相位的跟踪
Angrisani et al. Estimating ultrasonic time-of-flight through quadrature demodulation
Moon et al. Assessment on shock pressure acquisition from underwater explosion using uncertainty of measurement
CN103733033A (zh) 轨迹之间的线性关系
AU2012325067B2 (en) Filling level measuring device and method for determining a functional relationship between different tracks
US20170268921A1 (en) Method and apparatus for detecting the level of a medium
Jain et al. Evaluation of SAMOSA3 adapted retracker using Cryosat-2 SAR altimetry data over the Arctic ocean
Angrisani et al. Ultrasonic-based distance measurement through discrete extended kalman filter
Gerasimov et al. Increasing the Efficiency of the Probing Signals Correlation Processing of an Ultrasonic Liquid Flowmeter

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for 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: 20170613