王 宇,于靜娜,陳先中,侯慶文
(1.北京科技大學(xué) 自動(dòng)化學(xué)院,北京 100083;2.北京科技大學(xué) 工業(yè)過程知識(shí)自動(dòng)化教育部重點(diǎn)實(shí)驗(yàn)室,北京 100083)
進(jìn)入21世紀(jì),鋼鐵工業(yè)高能耗生產(chǎn)的現(xiàn)狀影響著中國經(jīng)濟(jì)的可持續(xù)發(fā)展。高爐煉鐵是高耗能、高排放的產(chǎn)業(yè),生產(chǎn)過程中料面高度等信息參數(shù)是節(jié)能減排的參考指標(biāo)和安全生產(chǎn)的有力保證[1]。
傳統(tǒng)料面監(jiān)測(cè)方法有紅外成像[2]、激光技術(shù)[3]等,但由于受高爐內(nèi)部高粉塵、高溫、高壓的復(fù)雜惡劣環(huán)境影響,致使紅外線測(cè)量精度低、成像質(zhì)量差;激光穿透粉塵能力較低且成本較高。傳統(tǒng)監(jiān)測(cè)方法不足以獲取高爐料面的有效信息。微波雷達(dá)具有傳輸距離長、穿透力強(qiáng)、可集成、適用范圍廣等優(yōu)點(diǎn),受溫度、粉塵影響相對(duì)較小,技術(shù)相對(duì)成熟,被廣泛應(yīng)用在精確成像中?,F(xiàn)有的高爐成像雷達(dá)中,ZANKL等[4]設(shè)計(jì)的高爐雷達(dá)使用貼片單元構(gòu)成,天線增益較低,在高爐內(nèi)的惡劣環(huán)境中效果不好。使用喇叭天線的機(jī)械擺動(dòng)雷達(dá)[5]尺寸較大,而高爐開孔的空間有限,無法組成陣列,在料面測(cè)量中僅能測(cè)得距離向信息,且1.6 GHz的工作帶寬較窄,距離向分辨率為9.375 cm。使用介質(zhì)天線的多輸入多輸出(Multi-Input Multi-Output,MIMO)雷達(dá)[6]造價(jià)昂貴,在低頻下的高爐雷達(dá)測(cè)量效果不理想[7],且組成天線陣列后整體尺寸較大,復(fù)雜的機(jī)械結(jié)構(gòu)極大地增大了雷達(dá)維護(hù)的難度。天線帶寬直接影響高爐料面成像的距離向分辨率,低帶寬不能滿足料面成像的高分辨率要求,且傳統(tǒng)高爐雷達(dá)大多體積較大,而高爐爐頂開口空間受限,所以傳統(tǒng)高爐雷達(dá)不能滿足對(duì)高爐料面監(jiān)測(cè)高集成度、小型化、高分辨率等需求。
微帶天線具有重量輕、成本低、體積小的特點(diǎn),能夠有效減小傳感器尺寸,易與電路集成,可通過陣列設(shè)計(jì)提高增益。YIN等[8]針對(duì)5G毫米波應(yīng)用,設(shè)計(jì)了一種對(duì)稱E形貼片天線,外型小巧,復(fù)雜度低,可用于無線通信系統(tǒng)。劉凡等[9]設(shè)計(jì)了一款新型雙圓極化平面微帶縫隙天線,通過切槽和寄生元件設(shè)計(jì)改進(jìn)了天線的性能。但是微帶天線的帶寬一般較窄,不適用于高爐料面檢測(cè)。為擴(kuò)寬帶寬,陳青青等[10]設(shè)計(jì)了一種寬帶圓極化微帶天線,展寬了天線軸比帶寬,且天線增益性能較好。鄭庚琪等[11]設(shè)計(jì)了一種新型寬帶小型化微帶天線,通過刻蝕縫隙和缺陷的結(jié)構(gòu)展寬了天線帶寬。CAO等[12]通過平面寄生貼片設(shè)計(jì),設(shè)計(jì)了具有寬帶寬和高增益的新型微帶貼片天線。但是現(xiàn)有的寬帶微帶天線不適合高爐料面檢測(cè)所需要的組陣需求,針對(duì)多輸入多輸出陣列的高集成度寬帶設(shè)計(jì)研究較少,且微帶貼片單元增益多為5~7 dBi,增益較低。串聯(lián)饋電的陣列形式可以提高天線增益,便于實(shí)現(xiàn)高集成度的線性多輸入多輸出陣列設(shè)計(jì)需求。多輸入多輸出合成孔徑雷達(dá)(Multi-Input Multi-Output Synthetic Aperture Radar,MIMO-SAR)是一種新體制雷達(dá),兼顧了多輸入多輸出雷達(dá)和合成孔徑雷達(dá)的技術(shù)優(yōu)點(diǎn)[13-14],提升了料面數(shù)據(jù)采集的效率,波數(shù)域算法[15]在低波段下通用性強(qiáng),成像精度高,成像質(zhì)量好,能夠有效地監(jiān)測(cè)高爐料面的實(shí)時(shí)信息。
考慮到高爐開口空間受限的問題,筆者選用微帶天線實(shí)現(xiàn)高集成度、小型化雷達(dá)的天線設(shè)計(jì)。針對(duì)微帶天線增益較低的問題,設(shè)計(jì)了串饋微帶陣列天線,提高增益至約14.05 dBi。從高爐料面檢測(cè)的高分辨率成像需求出發(fā),創(chuàng)新性地通過寄生貼片和空氣層設(shè)計(jì)并結(jié)合串饋陣列展寬了帶寬,提高了距離向分辨率。通過波數(shù)域成像算法對(duì)MIMO-SAR系統(tǒng)進(jìn)行了仿真分析,對(duì)比了不同帶寬下對(duì)模擬高爐料線的成像結(jié)果,驗(yàn)證了天線模型和成像監(jiān)測(cè)的有效性,為高爐料面監(jiān)測(cè)提供了一條新的解決思路。
串聯(lián)饋電天線的結(jié)構(gòu)簡單,空間占用少,易調(diào)節(jié)阻抗匹配和電流分布,是常用的微帶陣列天線設(shè)計(jì)方法。通過矩形微帶貼片單元組成串饋微帶陣列天線,在提高天線增益的同時(shí),串饋天線的物理形狀便于組成線性多輸入多輸出陣列;通過寄生貼片和空氣層設(shè)計(jì)可展寬帶寬,進(jìn)而提高高爐料面監(jiān)測(cè)的距離向分辨率。
圖1 微帶陣列天線示意圖
微帶陣列天線示意圖如圖1所示。選用Rogers 4350B作為介質(zhì)基板材料,介質(zhì)基板厚度h為0.508 mm,相對(duì)介電常數(shù)εr為3.66,損耗角正切為0.004?;迳戏綖橘N片單元,單元長寬分別為L和W,基板下方為地層,饋電方式選擇微帶線饋電,饋電點(diǎn)選擇輻射貼片中點(diǎn),通過饋電線連接各陣元。
采用傳輸線模型理論,可以用下式[16]計(jì)算出矩形貼片單元的初始尺寸。假設(shè)介質(zhì)的介電常數(shù)為εr,光速為c,對(duì)于工作頻率為f的矩形微帶天線,設(shè)計(jì)出輻射貼片的寬度W為
(1)
輻射單元長度L一般取λe/2,λe是介質(zhì)內(nèi)的導(dǎo)波波長??紤]到邊緣縮短效應(yīng)后,實(shí)際上的輻射單元長度L應(yīng)為
(2)
其中,εe是有效介電常數(shù),ΔL是等效輻射縫隙長度。它們可分別用下式計(jì)算:
(3)
(4)
在25 GHz的工作頻率下,計(jì)算得到初始值W=3.93 mm,L=2.91 mm。為保證陣列天線陣元各相位相同,陣元間饋線長度L2=λe/2,寬度初始值W2=0.3 mm。圖1中天線輸入端為四分之一波長阻抗變換器,長度L1=λe/4,寬度初始值W1=1 mm。四分之一波長阻抗變換器的寬度會(huì)影響陣列阻抗匹配程度,對(duì)回波損耗造成巨大影響;天線陣元長度會(huì)影響天線諧振頻率,陣元間饋線寬度會(huì)影響回波損耗,不利于阻抗匹配。
通過HFSS軟件對(duì)陣列天線進(jìn)行仿真和優(yōu)化設(shè)計(jì)。陣列天線回波損耗圖如圖2所示,陣列天線在23.8~25.4 GHz的工作頻率內(nèi)反射系數(shù)S11<-10 dB,絕對(duì)帶寬為1.6 GHz。陣列天線方向圖如圖3所示,天線增益為14.5 dBi,歸一化副瓣電平為-12.6 dB,E面半功率波束寬度為11.8°,H面半功率波束寬度為82°。
圖2 陣列天線反射系數(shù)
圖3 陣列天線方向圖
天線帶寬直接影響雷達(dá)成像的距離向分辨率。天線帶寬越寬,雷達(dá)距離向分辨率越高。筆者通過在陣列天線正上方覆蓋寄生貼片和增加空氣層展寬帶寬來改進(jìn)天線結(jié)構(gòu),如圖4所示。改進(jìn)后的天線為雙層結(jié)構(gòu),上下兩層均為Rogers 4350B介質(zhì)基板,下層介質(zhì)基板的上表面是1.1節(jié)所設(shè)計(jì)的微帶陣列天線(如圖1 所示),下表面是金屬接地板。上層介質(zhì)基板的上表面是8個(gè)寄生貼片,寄生貼片分別在陣列天線陣元的正上方,寄生貼片長寬為L3×W3,中間層為空氣層。上下層介質(zhì)基板厚度h、h1均為0.508 mm,空氣層高度h2為2 mm。
(a) 側(cè)視圖
(b) 俯視圖
對(duì)天線影響參數(shù)進(jìn)行分析。通過圖5可以看到,下層微帶陣列天線陣元長度L直接影響陣列天線的工作頻率,當(dāng)增大L時(shí),天線工作頻率向低頻偏移。通過圖6可以看到,天線上層寄生貼片長度會(huì)影響陣列天線的高頻諧振點(diǎn),當(dāng)增大L3時(shí),天線高頻諧振點(diǎn)向低頻偏移。
圖5 陣元長度對(duì)反射系數(shù)的影響
圖6 寄生貼片長度對(duì)反射系數(shù)的影響
對(duì)改進(jìn)后的微帶陣列天線進(jìn)行參數(shù)優(yōu)化,最后得到各參數(shù):L=2.96 mm,W=3.6 mm,L1=3.3 mm,W1=0.2 mm,L2=1.687 mm,W2=1.2 mm,L3=2.5 mm,W3=3.2 mm。
加入寄生貼片和空氣層的微帶陣列反射系數(shù)如圖7所示,反射系數(shù)S11低于-10 dB的工作頻率為20.67~25.92 GHz,絕對(duì)帶寬為5.25 GHz,較原天線帶寬展寬了3.65 GHz。寬帶陣列天線方向圖如圖8所示,天線在主瓣方向增益約為14.05 dBi,較原陣列天線略有降低。E面半功率波束寬度約為11.2°,歸一化副瓣電平約為-10.7 dB,H面半功率波束寬度約為84.6°。
圖7 寬帶陣列天線反射系數(shù)圖
圖8 寬帶陣列天線方向圖
圖9 天線主瓣增益隨頻率變化圖
對(duì)天線增益與工作頻率關(guān)系進(jìn)行統(tǒng)計(jì)。如圖9所示,天線在20.67~25.92 GHz工作頻率內(nèi),主瓣方向增益最大為14.05 dBi,最小為12.05 dBi,增益較高,可滿足測(cè)量需求。
圖10 等效相位中心示意圖
為驗(yàn)證寬帶微帶陣列天線的成像效果,根據(jù)等效相位中心原理[17],設(shè)計(jì)了線性多輸入多輸出陣列。通過波數(shù)域算法,對(duì)模擬高爐料線進(jìn)行了成像仿真實(shí)驗(yàn),為MIMO-SAR系統(tǒng)的成像測(cè)量提供理論支撐。
根據(jù)等效相位中心原理,如圖10所示,可將多輸入多輸出陣列由收發(fā)分置等效成收發(fā)同置的虛擬陣列,采用合成孔徑雷達(dá)成像算法進(jìn)行成像分析。一對(duì)發(fā)射、接收分置的陣元,可以等效為一個(gè)位于中心位置、收發(fā)同置的相位中心。等效的距離差會(huì)導(dǎo)致相位誤差,因此需要對(duì)等效相位中心進(jìn)行相位校正。一個(gè)1發(fā)N收間距為d的多輸入多輸出陣列,可等效為N個(gè)間距為d/2的收發(fā)同置的相位中心。
天線布陣選擇8發(fā)5收方式。定義間距d等于導(dǎo)波波長,即d=λe=6.75 mm。如圖11所示,圖(a)為多輸入多輸出雷達(dá)上層寄生貼片的俯視圖,圖(b)為下層微帶陣列天線的俯視圖,自左到右分別為發(fā)射天線T1、T2、T3、T4;接收天線R1、R2、R3、R4、R5;發(fā)射天線T5、T6、T7、T8。5個(gè)接收天線R1~R5在中心等間距均勻分布,間距為d;相鄰收發(fā)天線T4與R1、T5與R5分別間距為d。發(fā)射天線對(duì)稱排列在接收天線兩側(cè),T1~T4在接收陣列左側(cè),間距為5d稀疏分布;T5~T8在接收陣列右側(cè),間距為5d稀疏分布。8個(gè)發(fā)射天線與5個(gè)接收天線形成包含40個(gè)虛擬陣元的均勻陣列[18],虛擬陣元間距為0.5d。通過軟件仿真和優(yōu)化,當(dāng)天線模型長寬高尺寸約為26 cm×5.5 cm×0.4 cm時(shí),能夠滿足高爐開口的有限空間限制。
(a) 上層寄生貼片俯視圖
(b) 下層微帶陣列天線俯視圖
可通過式(5)計(jì)算出寬帶微帶陣列天線距離向分辨率:
(5)
其中,c為光速,B為天線帶寬。筆者選用20.9~25.9 GHz的頻率作為成像工作頻率,帶寬為5 GHz,可得δy=3 cm。
如表1所示,相比較于傳統(tǒng)高爐雷達(dá)中使用喇叭天線的機(jī)械擺動(dòng)雷達(dá)[5]及使用介質(zhì)天線的T型多輸入多輸出雷達(dá)[6],筆者使用寬帶微帶陣列天線設(shè)計(jì)的多輸入多輸出雷達(dá)拓寬了帶寬,提高了距離向分辨率,且雷達(dá)物理尺寸較小,集成度較高,能夠更好地滿足高爐爐頂開口空間有限的實(shí)際要求。
表1 高爐雷達(dá)對(duì)比
對(duì)于一個(gè)收發(fā)分置的M發(fā)N收陣列天線,天線方向圖受天線陣元激勵(lì)、位置影響。若不考慮天線單元輻射特性,則陣列中發(fā)射與接收天線的陣因子方向圖為
(6)
其中,(xti,yti)和(xri,yri)分別是發(fā)射陣列和接收陣列第i個(gè)天線陣元的位置坐標(biāo);ati和ari分別是發(fā)射陣列和接收陣列中第i個(gè)天線陣元的激勵(lì)電流;u=sinθcosφ,v=sinθsinφ,θ和φ分別是陣列天線遠(yuǎn)場輻射平面波與陣列法線方向的夾角。多輸入多輸出雷達(dá)方向圖等效于發(fā)射陣列與接收陣列方向圖的克羅內(nèi)克積[19],即
F(θ,φ)=Ft(θ,φ)?Fr(θ,φ) 。
(7)
(b) 接收陣列方向圖
(c) 收發(fā)陣列方向圖
對(duì)2.1節(jié)中的多輸入多輸出陣列進(jìn)行仿真,得到發(fā)射陣列、接收陣列以及收發(fā)陣列方向圖,如圖12所示。從圖12(c)中可以看到,方向圖主瓣寬度較窄,-3 dB波束寬度約為2°,旁瓣電平低于-20 dB,具有較強(qiáng)的指向性,能量集中,可滿足多輸入多輸出陣列設(shè)計(jì)需求。
圖13 MIMO-SAR成像模型
為驗(yàn)證多輸入多輸出陣列布陣的有效性,采用波數(shù)域算法[20]對(duì)目標(biāo)進(jìn)行成像仿真實(shí)驗(yàn)。如圖13所示,被測(cè)目標(biāo)在x-y平面上。天線陣列在方位向上沿x′軸對(duì)被測(cè)目標(biāo)進(jìn)行采樣,通過等效相位中心原理,在x′軸上形成間距為0.5d的等效虛擬采樣陣列;在距離向上,天線發(fā)射步進(jìn)頻率信號(hào)h(t)=exp(-j2πft),其中,f=f0+(n-1)Δf,為瞬時(shí)頻率,f0為初始頻率,Δf為步進(jìn)頻率。結(jié)合筆者設(shè)計(jì)的寬帶微帶陣列天線,選擇20.9~25.9 GHz作為成像工作頻率。通過天線陣列接收到的回波數(shù)據(jù),計(jì)算出被測(cè)目標(biāo)的散射強(qiáng)弱分布,進(jìn)行方位向和距離向的二維成像仿真實(shí)驗(yàn)。
波數(shù)域算法計(jì)算流程如下式[7]所示:天線在采樣位置測(cè)量的目標(biāo)回波信號(hào)S(xt,xr,y0)為
(8)
其中,k=2πf/c。xt、xr分別是發(fā)射天線、接收天線在x′軸位置坐標(biāo),y0為天線陣列對(duì)目標(biāo)中心的垂直距離。待測(cè)目標(biāo)位于x軸上,忽略目標(biāo)散射點(diǎn)相互作用后,設(shè)目標(biāo)散射強(qiáng)度為σ(x,y=0)。
對(duì)回波信號(hào)進(jìn)行相位補(bǔ)償校正,可得
(9)
其中,xm=(xt+xr)/2,(xm,y0)是收發(fā)陣元的等效相位中心位置坐標(biāo)。
對(duì)回波信號(hào)進(jìn)行降采樣后做方位向快速傅里葉變換(FFT),得到
S(kx,k)=TFFT{S(xm,k)} 。
(10)
計(jì)算等效傳播的空間濾波器函數(shù):
(11)
由于信號(hào)經(jīng)收發(fā)天線為雙程作用,變量滿足色散關(guān)系式:
(12)
對(duì)測(cè)量信號(hào)的空間譜信號(hào)進(jìn)行空間濾波和插值處理,可得兩個(gè)正交空間頻率方向的均勻譜域分布數(shù)據(jù)。根據(jù)頻域譜圖,截取出目標(biāo)特征最集中的區(qū)域,通過頻率補(bǔ)零對(duì)圖像進(jìn)行細(xì)化:
(13)
對(duì)S(kx,ky)做二維逆快速傅里葉變換(IFFT),求得對(duì)應(yīng)空間坐標(biāo)位置的被測(cè)目標(biāo)二維散射強(qiáng)度Sr(x,y),畫圖可得被測(cè)目標(biāo)的二維成像圖:
(14)
為模擬高爐料線形狀,筆者對(duì)V-平臺(tái)形狀[7]進(jìn)行模擬料線成像仿真實(shí)驗(yàn)。被測(cè)目標(biāo)中心作為空間坐標(biāo)原點(diǎn),陣列在距離向上距離坐標(biāo)原點(diǎn)的距離為y0。模擬料線由7個(gè)點(diǎn)目標(biāo)構(gòu)成,點(diǎn)目標(biāo)1號(hào)到點(diǎn)目標(biāo)7號(hào)的坐標(biāo)分別為(-0.2,-0.05),(-0.125,-0.05),(-0.0625,0),(0,0.05),(0.0625,0),(0.125,-0.05),(0.2,-0.05)。多輸入多輸出雷達(dá)的模擬料線成像仿真圖如圖14(a)所示,對(duì)被測(cè)目標(biāo)成功聚焦,得到了清晰的點(diǎn)目標(biāo)成像,能夠分辨出不同點(diǎn)目標(biāo),并得到模擬料線形狀信息,驗(yàn)證了天線陣列成像的可行性。帶寬直接影響高爐料面監(jiān)測(cè)的距離向分辨率。為直觀展示帶寬對(duì)高爐料面監(jiān)測(cè)的重要意義,以1.6 GHz帶寬為例,圖14(b)為多輸入多輸出雷達(dá)工作在窄帶寬下的成像效果,當(dāng)帶寬較窄時(shí),3號(hào)、4號(hào)、5號(hào)點(diǎn)目標(biāo)的成像有重疊部分,難以區(qū)分不同的點(diǎn)目標(biāo),不能對(duì)模擬料線進(jìn)行有效成像。對(duì)比圖(a)和圖(b)可知,筆者設(shè)計(jì)的寬帶微帶陣列天線拓寬了帶寬,能夠面向高爐料面進(jìn)行精度更高的成像監(jiān)測(cè)。
(a) 多輸入多輸出雷達(dá)成像結(jié)果
(b) 1.6 GHz帶寬雷達(dá)成像結(jié)果
圖14(c)對(duì)理想模擬料線和仿真模擬料線進(jìn)行了直觀對(duì)比。筆者設(shè)計(jì)的多輸入多輸出雷達(dá)對(duì)模擬料線的仿真點(diǎn)目標(biāo)成像坐標(biāo)如表2所示,可以看到,與理想模擬料線相比,仿真模擬料線的方位向平均誤差為0.008 0 m,距離向平均誤差為0.001 1 m,誤差較小,精度較高,可以準(zhǔn)確地獲得料面形狀信息。
表2 模擬料線成像結(jié)果對(duì)比 m
針對(duì)高爐料面監(jiān)測(cè)的高分辨率和高集成度需求,通過寄生貼片和空氣層展寬了微帶陣列天線的帶寬。通過仿真結(jié)果分析,天線在20.67~25.92 GHz頻率內(nèi)反射系數(shù)S11<-10 dB,絕對(duì)帶寬為5.25 GHz,主瓣增益可達(dá)14.05 dBi。雷達(dá)距離向分辨率提高至3 cm,通過多輸入多輸出陣列設(shè)計(jì)和波數(shù)域算法實(shí)現(xiàn)了MIMO-SAR成像仿真驗(yàn)證。從模擬料線的仿真成像對(duì)比結(jié)果可知,對(duì)比傳統(tǒng)高爐雷達(dá),筆者設(shè)計(jì)的天線可以準(zhǔn)確地反映料面形狀信息,準(zhǔn)確地區(qū)分不同點(diǎn)目標(biāo),成像方位向平均誤差為0.008 0 m,距離向平均誤差約為0.001 1 m,成像誤差小,且天線尺寸小,集成度高,具有一定的工程意義,為今后高爐雷達(dá)設(shè)計(jì)和料面監(jiān)測(cè)提供了一種新思路。