計量級雷達(dá)液位計高精度測量技術(shù)
摘要：計量級高精度測量對雷達(dá)液位計的系統(tǒng)設(shè)計和處理方法等提出了諸多挑戰(zhàn), 調(diào)頻連續(xù)波雷達(dá)通過測量發(fā)射信號和回波信號的頻率差來測量距離, 能夠獲得較高的測量精度, 同時由于其易于實現(xiàn), 是雷達(dá)液位計的主要體制之一。在介紹FMCW雷達(dá)液位計的系統(tǒng)原理基礎(chǔ)上, 對諸如晶振頻率、調(diào)頻斜率、工作環(huán)境等影響測量精度的因素進(jìn)行了分析, 提出了相應(yīng)的解決措施, 并開展了高精度頻譜估計技術(shù)研究。該項技術(shù)在實際的雷達(dá)液位計上得到了應(yīng)用, 其測量精度滿足系統(tǒng)要求。
High Accuracy Measurement Technology of Metrological Radar Level Meter
XIE Jinhua GUO Mingming
AVIC Leihua Electronic Technology Research Institute
Abstract：
Metrological high accuracy measurement presents many challenges to the system design and processing of radar level meter.LFCW radar measures distance by estimating the difference frequency between the transmitted wave and echo, high precision measurement can be obtained.And also it is easy to implement, so it becomes one of the main systems of radar level meters.On the basis of introducing the system principle of radar level meter, the influences of the crystal-controlled frequency, the slope of frequency modulation, the speed of light, and the work environment on the measurement accuracy are analyzed, and the corresponding solutions are put forward.Finally, a method of high accuracy range measurement technology is presented.This technology has been applied to the practical radar level meter and the measurement accuracy meets the system requirement.
0、引言：
油庫儲油罐感知層由溫度、壓力、流量、液位等信息傳感器組成, 雷達(dá)液位計是儲油罐油量的非接觸式測量裝置, 采用發(fā)射—反射—接收的工作模式探測液位高度, 具有安全可靠、精度高、壽命長、適應(yīng)各種工作環(huán)境的特點, 廣泛應(yīng)用于石化、倉儲行業(yè), 是儲油罐的核心傳感器[1]。
雷達(dá)液位計的發(fā)展是從20世紀(jì)80年代開始的, 雷達(dá)液位計主要由調(diào)頻連續(xù)波 (FMCW) 雷達(dá)、脈沖雷達(dá)兩類[2], 一般工作在C, X, K波段。眾所周知, 雷達(dá)的測距精度主要取決于信號的帶寬和信噪比, 脈沖雷達(dá)液位計由于受脈沖寬度和高速計時器件的限制, 測量精度一般在3m～10mm, 主要用于過程控制, 為中低端雷達(dá)液位計。與脈沖雷達(dá)相比, FMCW雷達(dá)液位計通過測量發(fā)射和目標(biāo)反射的FMCW信號的頻率差來測量距離, 可以形成大帶寬的FMCW信號且不受脈寬的限制, 因此能達(dá)到更高的測量精度和分辨率, 目前FMCW雷達(dá)液位計測量精度可達(dá)0.4m～3mm, 計量級精度則在±1mm以內(nèi), 主要用于計量級和商品交易, 為中高端應(yīng)用場合使用[3]。
對于中遠(yuǎn)程雷達(dá)來說, 其測距精度一般在米級水平, 諸如頻率精度、波形線性度、大氣溫度和壓力等因素對精度的影響微乎其微, 在設(shè)計時這些因素一般不予考慮。而在超近程高精度的雷達(dá)液位計中, 這些因素則是影響測距精度的主要方面。例如, 在雷達(dá)液位計中的100MHz晶振具有10-5的頻率精度, 經(jīng)過上行鏈路 10倍頻到X或K波段, 則頻率精度變?yōu)?0-4, 在測量20m的目標(biāo)時就有2mm的誤差。因此, 計量級高精度測量對雷達(dá)液位計的系統(tǒng)設(shè)計和處理方法上提出了諸多挑戰(zhàn)。本文從FMCW雷達(dá)液位計系統(tǒng)角度出發(fā), 針對這些因素對測距精度的影響進(jìn)行了深入研究, 并提出了可行的技術(shù)手段。
1、系統(tǒng)原理和組成：
1.1、FMCW雷達(dá)測距原理：
線性調(diào)頻連續(xù)波 (FMCW) 雷達(dá)是一種通過對發(fā)射FMCW信號及相應(yīng)的接收處理獲得距離與速度信息的雷達(dá)體制, 一方面發(fā)射信號通過天線向外輻射, 另一方面作為混頻器的基準(zhǔn)信號解調(diào)回波信號, 其得到的差頻信號頻率與目標(biāo)的距離成正比, 式 (1) 表示差頻信號頻率與目標(biāo)的距離的關(guān)系。雷達(dá)波形及測距原理如圖1所示。
圖1 FMCW雷達(dá)波形及測距原理圖
式中, fb為差頻信號頻率, R為目標(biāo)距離, B為調(diào)頻帶寬, T為調(diào)頻周期, c為光速。
1.2、FMCW雷達(dá)液位計組成：
FMCW雷達(dá)液位計由喇叭天線、收發(fā)模塊、電源模塊、信號處理模塊、顯示模塊等單元組成。收發(fā)模塊產(chǎn)生FMCW信號, 經(jīng)環(huán)行器至喇叭天線向外輻射, 回波信號由喇叭天線接收經(jīng)環(huán)行器送到接收通道, 接收通道對回波信號進(jìn)行諧波混頻、放大、低通濾波后, 輸出包含距離信息的差頻信號, 經(jīng)A/D采樣形成數(shù)字信號, 信號處理模塊接收收發(fā)模塊的差頻數(shù)字信號進(jìn)行FFT、雜波抑制、目標(biāo)檢測、頻譜估計等處理后, 得到高精度的液位數(shù)據(jù), 通過現(xiàn)場總線送到顯示單元和上位機(jī), 其原理框圖如圖2所示。
圖2 FMCW雷達(dá)液位計原理框圖
2、測距精度影響因素分析及相應(yīng)的解決措施：
根據(jù)FMCW信號測距精度公式[4]:
測距精度主要取決于信號帶寬和信噪比, 假定信號帶寬B為2GHz, 要達(dá)到1mm的精度, 則信噪比 (SNR) J=29dB, 式 (2) 為雷達(dá)測距精度的Cramer-Rao界, 因此在實際系統(tǒng)中信噪比應(yīng)大于29dB。由此看出, 要滿足1mm的測距精度, 雷達(dá)液位計應(yīng)具有大帶寬、高信噪比的特征。
在雷達(dá)液位計系統(tǒng)設(shè)計滿足測距精度情況下, 最終得到的距離精度還受其他因素影響, 由式 (1) 變換得
式中, K=B/T為調(diào)頻斜率。從式 (3) 可以看出, 差頻頻率、調(diào)頻斜率將會影響測距精度, 誤差表現(xiàn)形式為比例誤差, 下面對這些因素對測距精度的影響進(jìn)行分析, 并提出相應(yīng)的技術(shù)措施。
2.1、晶振頻率：
差頻頻率fb和帶寬B的誤差的主要來源于晶振頻率的誤差, 該誤差主要由中心頻率偏差和溫度變化引起的頻率偏差兩部分構(gòu)成, 前者用頻率準(zhǔn)確度表示, 后者用頻率溫度穩(wěn)定度表示。對于晶振中心頻率偏差, 可在雷達(dá)系統(tǒng)校正時進(jìn)行比例修正, 盡量選擇頻率準(zhǔn)確度較高的器件, 目前晶振的頻率準(zhǔn)確度 (與晶振輸出頻率的比值) 一般在1×10-5以內(nèi);對于溫度變化引起的頻率偏差, 為了減小溫度變化的影響, 需選用恒溫晶振, 其溫度穩(wěn)定度指標(biāo)可達(dá)1×10-9;考慮到價格因素, 也可采用價格便宜的溫補(bǔ)晶振, 其溫度穩(wěn)定度指標(biāo)可達(dá)5×10-7, 可以滿足精度要求。
2.2、FMCW調(diào)頻斜率：
FMCW信號的調(diào)頻斜率的線性度指標(biāo)直接關(guān)系到距離測量精度, 早期的FMCW雷達(dá)液位計大多采用VCO產(chǎn)生FMCW信號, 由于VCO是模擬電壓控制頻率產(chǎn)生, 易受模擬電壓非線性和溫度變化的影響, FMCW信號的線性度不高, 要采用復(fù)雜的校正控制電路和溫度補(bǔ)償電路。而數(shù)字直接頻率合成技術(shù) (DDS) 采用數(shù)字方法產(chǎn)生頻率, 線性度極高且受溫度影響較小, 在雷達(dá)和通信系統(tǒng)中得到了廣泛的應(yīng)用, 且芯片價格已下降到合理的水平。因此, 在設(shè)計中需采用DDS來產(chǎn)生FMCW信號, 以保證FMCW信號的線性度。
DDS具有良好的線性度, 其等效線性度在10-6以下, 根據(jù)線性度誤差的經(jīng)驗公式[3]:
式中, L為線性度, R為目標(biāo)距離, 設(shè)L=10-6, R=20m, 則δL=0.08mm。
采用DDS后, 其調(diào)頻線性度對測距精度基本沒有影響, 但DDS的量化誤差對測距精度有影響, 該影響也可通過系統(tǒng)校正進(jìn)行比例修正。
2.3、 系統(tǒng)校正：
晶振頻率、調(diào)頻斜率等的誤差對測距精度有較大的影響, 其表現(xiàn)形式為比例誤差, 在系統(tǒng)設(shè)計時要進(jìn)行定量分析, 其定量分析數(shù)據(jù)對系統(tǒng)調(diào)試有指導(dǎo)作用, 但最終需要通過測量系統(tǒng)來對雷達(dá)液位計進(jìn)行精確的比例校正, 考慮到比例誤差的影響, 式 (3) 修正為
式中, g (fr, K) 為校正系數(shù), 與晶振頻率、調(diào)頻斜率等的誤差有關(guān)。
雷達(dá)液位計測量系統(tǒng)主要由激光跟蹤儀、儲油罐模擬板、液位計支架 (液位計安裝車) 、監(jiān)控計算機(jī)等組成, 放置在滿足一定長度要求的暗室內(nèi), 其原理框圖如圖3所示。在該系統(tǒng)中, 激光跟蹤儀測量雷達(dá)液位計與目標(biāo)板之間的距離, 并以此距離值為基準(zhǔn), 與雷達(dá)液位計測量的距離值進(jìn)行比對來計量雷達(dá)液位計的精度, 通過設(shè)置多個測量點來進(jìn)行比例校正, 具體方法如下:
1) 在量程范圍內(nèi)均勻設(shè)置N個測量點;
2) 記錄激光跟蹤儀的N點測量值Rl, n, n=1, 2, …, N, 計算N-1點距離差分值:
3) 記錄雷達(dá)液位計的N點測量值Rr, n, n=1, 2, …, N, 計算N-1點距離差分值:
4) 計算比例校正系數(shù):
5) 計算比例校正系數(shù)的均值, 代入式 (5) , 完成比例校正。
由于雷達(dá)液位計內(nèi)部的射頻電路和電纜在傳輸信號時有一定的延時, 因此雷達(dá)液位計基準(zhǔn)面的距離值不等于零, 存在固定誤差, 該固定誤差可通過雷達(dá)液位計測量系統(tǒng)進(jìn)行測量, 在系統(tǒng)計算時予以消除。
圖3 雷達(dá)液位計測量系統(tǒng)
3、高精度頻譜估計技術(shù)：
從前面的分析可知, 要達(dá)到1mm的測量精度, 雷達(dá)液位計需要具備大帶寬、高信噪比的特點。傳統(tǒng)測量方法通過數(shù)字化和FFT計算得到距離量化值, 但是FFT的分辨率是顯然不能滿足測量精度要求的。例如, 假定信號帶寬為2GHz, 則FFT的距離分辨率為75mm。這是因為檢測頻譜上的高點不一定是回波信號頻譜的 實際峰值, 直接采用FFT的離散頻譜點來計算距離會有較大的誤差。因此, 需要通過高精度頻譜估計來找到兩離散頻譜高點中間的峰值位置, 以得到高精度的距離測量值。
為了滿足高精度測量需求, 學(xué)者們提出了多種提高FMCW雷達(dá)距離測量精度的方法。例如:
1) 利用DFT頻譜相位信息實現(xiàn)高精度測距[5], 該方法運算量大, 對硬件要求高, 工程化實現(xiàn)難度較大;
2) 結(jié)合Chirp-Z變換, “由粗到細(xì)”逐步找到頻譜實際峰值的位置[6], 該方法的本質(zhì)是通過增加采樣點的方式提高精度, 增加了運算量;
3) 找到頻譜主瓣內(nèi)最大和次大譜線, 使用插值法求出中心頻率[7], 該方法原理簡單, 計算量小, 易于工程實現(xiàn), 但其提高仿真精度的能力有限。
針對上述問題, 本文提出一種基于迭代插值的高精度譜估計算法, 該算法首先采用FFT粗測差頻信號的頻率, 然后在FFT峰值附近取兩條譜線進(jìn)行頻率細(xì)化。其核心思想為:找到回波序列x (n) 的頻譜X (k) 幅度最大值A(chǔ)和次大值B, 及其所在的位置k1, k2;若A=B, 則目標(biāo)真實位置所對應(yīng)的頻率就在k1, k2的中心。其原理如圖4所示。
圖4 高精度譜估計算法示意圖
實現(xiàn)方法為:找到X (k) 幅度最大值所對應(yīng)的位置k1, 選出k1左右兩點中幅度較大的那個點位置記為k2, 即k2=k1+1或k2=k1-1;根據(jù)這兩個分量, 利用式 (8) 計算出第一次迭代結(jié)果k1 (1) :
式中, 用上標(biāo)來表示迭代序號。根據(jù)式 (8) 計算出的位置, 利用式 (9) 可以進(jìn)一步確定一組新的
通過式 (10) 計算k1 (1) , k2 (1) 這兩點對應(yīng)的頻率分量的幅度值, 分別為
得到X (k1 (1) ) 和X (k2 (1) ) 之后, 繼續(xù)利用式 (8) 計算出k (2) , 按照同樣的步驟, 迭代3次, 最后得到k1 (4) , k2 (4) , 而最優(yōu)頻率分量對應(yīng)的位置k*由式 (11) 計算得到:
這種方法測量精度高, 而且最大的優(yōu)點是運算量小, 單次測量只需進(jìn)行6次DFT。
設(shè)定一組仿真參數(shù):帶寬B=2GHz, 頻率調(diào)制時間T=5.12ms, 采樣率fs=200kHz, 信噪比SNR=29dB。使用上述算法進(jìn)行仿真, 對20.8～21m以1mm采樣間隔得到的距離進(jìn)行測試, 得到誤差統(tǒng)計結(jié)果, 如圖5所示。
圖5 本文算法測距誤差統(tǒng)計結(jié)果
從圖5可知, 利用本文提出的基于迭代插值的高精度譜估計算法進(jìn)行測距, 不同距離的測量誤差都控制在1mm以內(nèi)。經(jīng)統(tǒng)計, 誤差絕對值的最大值為0.153 2 mm, 誤差絕對值的最小值為2.987 8×10-5 mm;誤差均值為0.001 2mm, 標(biāo)準(zhǔn)差為0.045 8 mm;誤差絕對值的均值為0.036 4mm, 誤差絕對值的標(biāo)準(zhǔn)差為0.027 7mm;證明了該算法可以滿足高精度測距的要求。
4、結(jié)束語：
本文從FMCW雷達(dá)液位計系統(tǒng)設(shè)計角度出發(fā), 根據(jù)計量級測量精度的要求, 分析得出雷達(dá)液位計應(yīng)具有大帶寬高信噪比的特點;在此基礎(chǔ)上, 對諸如晶振頻率、調(diào)頻斜率、工作環(huán)境等影響測量精度的因素進(jìn)行了分析, 提出了相應(yīng)的技術(shù)措施, 得到了其影響為比例誤差的結(jié)論;通過設(shè)計雷達(dá)液位計測量系統(tǒng), 對比例誤差進(jìn)行校正, 給出了具體的測試方法;最后研究了一種高精度頻譜估計算法。本文研究的FMCW雷達(dá)液位計高精度測量技術(shù)在實際的雷達(dá)液位計得到了應(yīng)用, 并通過了1mm精度的國家計量認(rèn)證, 驗證了該技術(shù)的有效性。
相關(guān)產(chǎn)品推薦：
