李 彬 金 銘 白 明 李志平

(1.中國科學(xué)院國家空間科學(xué)中心中國科學(xué)院微波遙感技術(shù)重點(diǎn)實(shí)驗(yàn)室，北京 100190；2.中國科學(xué)院遙感與數(shù)字地球研究所遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，北京 100101；3.北京航空航天大學(xué)電子信息工程學(xué)院，北京 100191)

1 引 言

微波黑體定標(biāo)源是微波輻射計(jì)系統(tǒng)中的關(guān)鍵部件，為微波輻射定量測量提供標(biāo)準(zhǔn)的亮溫輻射參考[1～4]。特別對于長期定量遙感觀測應(yīng)用的星載微波輻射計(jì)，微波黑體定標(biāo)源是必不可少的在軌定標(biāo)器件[1,2]。為了同時(shí)實(shí)現(xiàn)接近于1的發(fā)射率和均勻的表面溫度分布，微波黑體定標(biāo)源常被設(shè)計(jì)為周期性涂覆尖錐陣列的結(jié)構(gòu)形式[5～8]。根據(jù)基爾霍夫熱平衡定律，在熱平衡狀態(tài)下，物體的電磁輻射與電磁吸收相同。因此，設(shè)計(jì)高發(fā)射率的微波黑體定標(biāo)源，也即設(shè)計(jì)高電磁吸收率低電磁反射率的結(jié)構(gòu)[8, 9]。同時(shí)，也這是微波黑體定標(biāo)源的發(fā)射率測定時(shí)所采取的基本思路，即通過測量微波黑體定標(biāo)源的電磁散射推知其電磁反射率[3,4,10～13]，進(jìn)而根據(jù)e=1-r得到其發(fā)射率測值。

微波輻射計(jì)定標(biāo)源在應(yīng)用之前，需要準(zhǔn)確的獲知其真實(shí)發(fā)射率，已為高精度的輻射計(jì)定標(biāo)提供依據(jù)。然而，在目前最可行的基于基爾霍夫熱平衡定律的發(fā)射率測定方法體系中(基于電磁散射測量)，存在兩個(gè)方面的難點(diǎn)問題。首先，是定標(biāo)源本身的發(fā)射率高，即其電磁反射率很低，電磁散射十分微弱，對散射測量系統(tǒng)和測量方式提出了很高的要求[10,11]。然后，無論采用單站還是雙站的測量方式，都難以按照理論概念獲得上半空間完備的散射分布信息，進(jìn)而直接積分獲取其電磁反射率，也即需要基于有限的測得散射信息推知積分反射率(r)，進(jìn)而得到發(fā)射率(e)[13]。

國內(nèi)外在微波黑體定標(biāo)源發(fā)射率/反射率測定方面開展了一系列研究工作[3,4,6,7,9,10～15]。其中文獻(xiàn)[10,13,14]，對基于單站和雙站散射測量的發(fā)射率測定方法都有所討論，而文獻(xiàn)[4,11,12]的研究主要集中在單站方式。雙站散射測量方式，雖然其測量系統(tǒng)構(gòu)成和實(shí)測過程相比單站更為復(fù)雜，但其對小散射信號的測量能力更強(qiáng)，獲取的散射信息也更多。目前，受科技部重點(diǎn)研發(fā)計(jì)劃的資助，中科院國家空間科學(xué)中心微波系統(tǒng)部正研制和構(gòu)建具備單站(鏡像)和雙站散射測量能力的發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置，頻率覆蓋6GHz～500GHz，已滿足未來長期的微波輻射計(jì)定標(biāo)源計(jì)量需求。

本文中，將依據(jù)定標(biāo)源發(fā)射率-電磁散射的理論關(guān)系，針對微波發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置的應(yīng)用需求，特別是應(yīng)對發(fā)射率測定方法體系中的難點(diǎn)問題，在毫米波和亞毫米波段(40GHz以上)探討對其照射場控制的原理和要求。并基于準(zhǔn)光(Quas-Optical)技術(shù)的原形天線設(shè)計(jì)[16]，通過仿真手段驗(yàn)證和總結(jié)準(zhǔn)光天線設(shè)計(jì)的應(yīng)用效能，同時(shí)也可對其它近場測試系統(tǒng)的研制和應(yīng)用提供參考。

2 微波黑體定標(biāo)源的發(fā)射率和電磁散射特性

2.1 發(fā)射率、反射率、電磁散射

根據(jù)基爾霍夫熱平衡定律，考察微波黑體定標(biāo)源沿正上方的發(fā)射率e，可以通過考察其對來自正上方的電磁來波的電磁反射率來實(shí)現(xiàn)(e=1-r)，而電磁反射率原理上由定標(biāo)源在上半空間散射系數(shù)的積分得到，如公式(1)所示[5～7,13～15]。

圖1 涂覆錐體陣列型定標(biāo)源發(fā)射率、電磁散射關(guān)系和結(jié)構(gòu)示意圖Fig.1 Configurations on the emissivity, electro-magnetic scattering and geometry of calibration target in shape of coated pyramid array

(1)

式中：S表示散射功率角分布(遠(yuǎn)場)，其根據(jù)入射功率歸一化即得到微分散射系數(shù)γ，進(jìn)一步進(jìn)行空間積分即可得到積分反射率r。

由此可知，理論上，要獲得定標(biāo)源結(jié)構(gòu)的發(fā)射率，需要獲知其在上半空間的全部散射分布信息進(jìn)行空間積分。在實(shí)際操作中，這幾乎是不可能完成的，更常見的情況是只能獲得部分的散射分布信息。

圖2高頻段(40GHz以上)發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置結(jié)構(gòu)示意圖Fig.2 Configuration of the Standard Measurement Facility on the Emissivity for the 40GHz+ frequency range

圖2是正在研制的發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置結(jié)構(gòu)示意圖。這是一個(gè)典型的雙站散射測量系統(tǒng)，其中發(fā)射天線相對待測黑體定標(biāo)源保持固定，而接收天線進(jìn)行俯仰角度(θ)的旋轉(zhuǎn)。發(fā)射天線呈傾斜照射(10°以內(nèi))，這樣接收天線在旋轉(zhuǎn)掃描時(shí)可以捕捉到鏡像散射瓣。發(fā)射/接收天線到目標(biāo)區(qū)的距離在1m左右。在該標(biāo)準(zhǔn)裝置中，典型的掃描測量模式是獲得俯仰掃描面內(nèi)的掃角散射結(jié)果(線掃)。

2.2 定標(biāo)源的散射特性

目前，微波黑體定標(biāo)源的典型結(jié)構(gòu)形式是周期性的涂覆尖錐陣列。前期的研究表明，由于定標(biāo)源結(jié)構(gòu)的周期性，其在上半空間的散射分布特性符合Floquet定律，即散射能量集中在各個(gè)散射瓣之中[5,6,8,15]。例如，正入射條件下各個(gè)散射瓣的角度位置，可以由公式(2)得到。傾斜入射的情況與正入射類似，只是需要計(jì)入入射波方向本身引起的波矢量投影。

(2)

具體而言，當(dāng)頻率足夠低，定標(biāo)源錐體周期p小于波長λ時(shí)，定標(biāo)源結(jié)構(gòu)只會產(chǎn)生鏡像瓣(后向瓣)一個(gè)散射主瓣，對應(yīng)上式中m，n=0。當(dāng)頻率逐漸升高，定標(biāo)源錐體周期p大于λ時(shí)，定標(biāo)源結(jié)構(gòu)會產(chǎn)生不止一個(gè)散射主瓣，即對應(yīng)上式中m，n≠0。此時(shí)，各個(gè)散射瓣皆攜帶散射功率，評估定標(biāo)源結(jié)構(gòu)的電磁反射率，需要計(jì)入所有散射瓣攜帶的總散射功率。

對于周期性結(jié)構(gòu)的定標(biāo)源而言，其在平面波的照射下，當(dāng)頻率足夠高時(shí)，將在遠(yuǎn)場散射區(qū)域呈現(xiàn)顯著的Floquet散射瓣，散射功率都集中在這些散射瓣之中。換言之，對于定標(biāo)源的積分反射率，并不需要完整地獲取上半空間的散射角分布；如能通過實(shí)測確定散射功率在各個(gè)散射瓣中的分配情況，足以準(zhǔn)確地獲取積分反射率進(jìn)而確定發(fā)射率真值[13～15]。

2.3 期望的天線照射效果

在發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置中，對收發(fā)天線的照射效果有以下兩個(gè)方面的要求：

(1) 聚束性。定標(biāo)源結(jié)構(gòu)是典型的低散射結(jié)構(gòu)。以往的研究表明，在毫米波和亞毫米波段(40GHz以上)定標(biāo)源的整體反射率都處于很低的水平。因此從保證小散射信號測量的角度，希望天線能將電磁波束聚束到待測目標(biāo)區(qū)上，這樣也同時(shí)避免了目標(biāo)下部支撐結(jié)構(gòu)的散射干擾。

(2)相位平坦性。從散射特征分析和反射率獲取的角度，期望得到定標(biāo)源結(jié)構(gòu)的遠(yuǎn)場散射特性，獲取其Floquet散射特征。然而，特別是在高頻段，發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置是典型的近場雙站測量系統(tǒng)。此時(shí)，期望天線能在目標(biāo)區(qū)形成相位平坦的照射效果，從后面的仿真結(jié)果中可看出，這樣的設(shè)置可近似測得周期結(jié)構(gòu)的遠(yuǎn)場Floquet散射特征。

上面兩個(gè)方面的要求在圖3中進(jìn)行了直觀的展示：即期望電磁波束能聚束在待測定標(biāo)源結(jié)構(gòu)之上，并且在照射范圍內(nèi)相位平坦。

圖3在待測微波黑體定標(biāo)源處期待達(dá)到的天線照射效果示意Fig.3 Preferred illumination performance at the microwave blackbody calibration target under test

在40GHz以上的頻率，為實(shí)現(xiàn)上面分析的照射效果，基于準(zhǔn)光學(xué)技術(shù)設(shè)計(jì)的聚焦式反射面式天線系統(tǒng)是可行的選擇。下面，將基于準(zhǔn)光天線的原形設(shè)計(jì)，驗(yàn)證和分析其在此近場雙站散射測量系統(tǒng)中的應(yīng)用效能。

3 準(zhǔn)光照射天線的適用性驗(yàn)證與討論

本小節(jié)中，將基于準(zhǔn)光天線原形設(shè)計(jì)的場照射效果，建立發(fā)射天線-假定周期目標(biāo)-接收天線的模擬仿真鏈路，模擬分析基于準(zhǔn)光照射天線的掃角測值曲線特性，如圖4所示。

圖4 準(zhǔn)光照射天線原形設(shè)計(jì)在目標(biāo)區(qū)處的照射效果，包括照射場的幅值(上)和相位(下)，50GHz(40～60GHz設(shè)計(jì))和60GHz(50～75GHz設(shè)計(jì))，線極化，主極化Y方向Fig.4 Illumination fields in the target zone of the quasi-optical antenna prototype, in amplitude (upside) and phase (downside), at 50GHz and 60GHz, linear polarization in Y

3.1 準(zhǔn)光天線原形設(shè)計(jì)的照射效果

首先，在本小節(jié)中，給出原形天線設(shè)計(jì)的照射場分布結(jié)果，如圖5所示。其中，天線反射鏡口徑為200mm(50GHz設(shè)計(jì)，40GHz～60GHz)和180mm(60GHz設(shè)計(jì)，50GHz～75GHz)，目標(biāo)區(qū)觀察口面的橫向尺寸為400mm× 400mm?？梢钥闯觯ㄟ^準(zhǔn)光反射面設(shè)計(jì)，可以在目標(biāo)區(qū)形成較為理想的相位平坦的聚束照射。本文中對天線照射場和后續(xù)目標(biāo)散射場的仿真，基于常用的商業(yè)軟件FEKO[17]，主要采用了其中的MFLMA(多層快速多極子算法)和Large element PO(大網(wǎng)格快速物理光學(xué)積分)算法。

3.2 天線-目標(biāo)-天線近場測試鏈路仿真

為分析準(zhǔn)光照射天線的應(yīng)用效能，建立了天線-目標(biāo)-天線的近場測試仿真分析鏈路，如圖6所示。其中，照射和接收天線效果通過準(zhǔn)光反射面出射的口面場分布引入。在目標(biāo)區(qū)設(shè)定了假定散射目標(biāo)，即周期性金屬小球陣列，用于模擬周期性錐體定標(biāo)源結(jié)構(gòu)的電磁散射特性。為了模擬錐體定標(biāo)源的立體結(jié)構(gòu)特性，設(shè)置了兩層小球陣列，分別對應(yīng)錐尖和錐底的位置，高寬比4∶1。

圖5 假定散射目標(biāo)示意，周期性小球陣列(周期p=15mm, 高度h=60mm, 小球半徑r=λ/2@50GHz)Fig.5 Assumed scattering target, 2 layer periodic metal sphere array, parameters include: period p=15mm, heighth = 60mm, r =λ/2@50GHz

圖6 照射天線-散射目標(biāo)-接收天線的近場散射測量仿真鏈路示意。其中照射和接收天線的效果由其口面場分布引入Fig.6 Configuration of simulations on the near-field scattering test link, starting from the transmitting antenna, to the assumed scatters, and ends at the receiving antenna.The transmitting and receiving antennas are included by using their radiated aperture field

鏈路分析的方法算法方面，首先，發(fā)射天線到目標(biāo)區(qū)的散射計(jì)算，基于FEKO軟件的MFLMA算法完成，得到目標(biāo)在準(zhǔn)光天線照射條件下在接收天線角度掃描接收位置處的近場散射場分布。然后，基于口徑耦合積分公式[13]見公式(3)，基于接收天線的口面場分布，得到接收天線在不同角度掃描位置處的等效接收功率。

(3)

式中：Pr——口徑天線的接收功率，即散射測值;Pr——口徑天線作為發(fā)射天線時(shí)的發(fā)射功率;Escat和Hscat——在耦合積分口面S上的散射場分布;Et和Ht——在耦合積分口面S上口徑天線等效為發(fā)射狀態(tài)時(shí)的口面場分布。

3.3 近場雙站散射測量鏈路的仿真結(jié)果

首先，在目標(biāo)區(qū)內(nèi)設(shè)置金屬平板(100mm×100mm)作為假定目標(biāo)，進(jìn)行上一小節(jié)中所描述的測試鏈路仿真。發(fā)射天線傾斜7.5°入射(-7.5°為后向接收角度位置，7.5°為鏡像接收角度位置)，考慮50GHz，鏈路仿真的結(jié)果如圖7所示。其中口徑接收測值曲線的結(jié)果形式是接收功率系數(shù)，即接收天線接收到的功率值除以發(fā)射天線的發(fā)射功率值。由于金屬平板的散射機(jī)理為鏡面反射，從掃角測值曲線的最大值(大于-5dB，約-2dB～-3dB)可以看出，發(fā)射天線-接收天線的鏈路衰減較小，這從測試靈敏度的角度對于低散射目標(biāo)的測量是十分有益的。這說明準(zhǔn)光照射天線可將主要能量集中在目標(biāo)之上，同時(shí)接收天線的口徑效應(yīng)也起到了能量收集的效果。

圖7 金屬平板目標(biāo)的雙站掃角測值曲線(鏈路仿真得到：在50GHz，測試距離1m)，發(fā)射天線7.5°傾斜照射Fig.7 Simulated bi-static scanning curves from metal plate, parameters include: at 50GHz, the test range is 1000mm, the incident angle is 7.5°for the transmitter

然后，考慮金屬小球陣列目標(biāo)，周期15mm，在圖8中給出了雙站掃角測值曲線結(jié)果，同時(shí)近場散射分布和遠(yuǎn)場散射方向圖也作為參考。由于天線照射場在目標(biāo)區(qū)的相位平坦性，金屬小球陣列的遠(yuǎn)場散射場分布體現(xiàn)了顯著的Floquet散射特性，即散射能量集中在有限的散射瓣之中。同時(shí)經(jīng)準(zhǔn)光天線口徑接收的測值曲線，可以很好地符合陣列的遠(yuǎn)場散射分布特性，如角度位置、測值動(dòng)態(tài)范圍等。因此，借助相位平坦的準(zhǔn)光照射天線，可以在近場條件下的微波發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置中近似測得周期性錐體定標(biāo)源結(jié)構(gòu)的遠(yuǎn)場散射特征。

圖8 金屬小球陣列目標(biāo)(p=15mm, 兩層小球)的雙站掃角測值曲線(鏈路仿真得到：在50GHz，測試距離1m)，發(fā)射天線7.5°傾斜照射Fig.8 Simulated bi-static scanning curves from two-layer metal spheres array (p = 15mm), parameters include: at 50GHz, the test range is 1000mm, the incident angle is 7.5°for the transmitter

以上兩個(gè)算例，一方面驗(yàn)證了準(zhǔn)光天線原形設(shè)計(jì)的照射性能，另一方面直觀地驗(yàn)證了發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置天線照射設(shè)計(jì)思路的有效性和準(zhǔn)確性。

下面，進(jìn)一步擴(kuò)展頻率，考慮500GHz處的情況，此時(shí)假定目標(biāo)金屬小球陣列的參數(shù)為周期10mm(依舊為兩層小球)。其在角度掃描軌道附近的散射場分布如圖9所示。對比50GHz時(shí)的結(jié)果，可以看出在500GHz，由于波長(0.6mm)遠(yuǎn)小于陣列周期(10mm)，散射場分布中體現(xiàn)出了密集的散射瓣，并出現(xiàn)了散射瓣和散射瓣間干涉現(xiàn)象。同樣進(jìn)行測試鏈路的仿真，如圖10所示，得益于準(zhǔn)光天線的口徑接收效果，口徑接收曲線中可以很好地體現(xiàn)各個(gè)散射瓣，得到互不干涉，動(dòng)態(tài)范圍大的Floquet散射分布特征。

圖9 金屬小球陣列目標(biāo)在近場雙站掃角軌道附近的近場散射場分布(軌道半徑1m)，50GHz和500GHz對比Fig.9 Simulated near-field scattering fields around the bi-static angular scanning track (radius=1000mm), at 50GHz and 500GHz

圖10 金屬小球陣列目標(biāo)(p=10mm, 兩層小球)的雙站掃角測值曲線(鏈路仿真得到：在500GHz，測試距離1m)，發(fā)射天線7.5°傾斜照射Fig.10 Simulated bi-static scanning curves from two-layer metal spheres array (p=10mm), parameters include: at 500GHz, the test range is 1000mm, the incident angle is 7.5°for the transmitter

4 結(jié)束語

針對微波發(fā)射率計(jì)量標(biāo)準(zhǔn)裝置的研制和發(fā)射率計(jì)量測試方法的研究需要，本文探討了其近場測試場景中照射天線的設(shè)計(jì)思路和方法。并借助商業(yè)仿真軟件FEKO，建立了天線-目標(biāo)-天線的測值仿真鏈路，模擬雙站的掃角測值量值特性。仿真結(jié)果驗(yàn)證了本文提出的照射天線設(shè)計(jì)思路的準(zhǔn)確性和有效性：通過聚束和相位平坦設(shè)計(jì)，可以降低測試鏈路中的能量衰減，同時(shí)可以在近場測試環(huán)境中近似獲得周期性定標(biāo)源結(jié)構(gòu)的遠(yuǎn)場散射特性，準(zhǔn)確地獲取其Floquet散射瓣特征以利于后續(xù)的反射率測值獲取。同時(shí)，本文中介紹的照射天線設(shè)計(jì)思路，對于一般性的近場散射測量也有一定的參考意義。