李 元，張 勇，胡麗琴，陸其峰，盧乃錳

（中國氣象局 國家衛(wèi)星氣象中心 中國遙感衛(wèi)星輻射測量和定標(biāo)重點(diǎn)開放實(shí)驗室，北京 100081）

1 引 言

自20世紀(jì)80年代以來，衛(wèi)星遙感逐漸向定量化方向發(fā)展。衛(wèi)星圖像不只用于目視解譯，還用來定量反演大氣、陸地和海洋等地球物理參數(shù)。為了確保在軌衛(wèi)星觀測的準(zhǔn)確性，各航天大國著手選擇輻射特性均勻穩(wěn)定的地面目標(biāo)，將其作為“標(biāo)準(zhǔn)”，通過衛(wèi)星過境時的地面同步觀測和輻射傳輸規(guī)律，對在軌衛(wèi)星遙感器進(jìn)行標(biāo)定，從而實(shí)現(xiàn)輻射校正。輻射校正已經(jīng)成為遙感領(lǐng)域的重大基礎(chǔ)研究課題，得到了國際遙感界的高度重視。我國在20世紀(jì)90年代啟動了“中國遙感衛(wèi)星輻射校正場”項目的科研和建設(shè)工作，開啟了我國遙感衛(wèi)星定量化技術(shù)系統(tǒng)性攻關(guān)。中國遙感衛(wèi)星輻射校正技術(shù)以輻射校正場為基礎(chǔ)，圍繞我國全系列國產(chǎn)遙感衛(wèi)星發(fā)展和定量化應(yīng)用需求，在輻射校正技術(shù)指標(biāo)體系、理論與方法、標(biāo)準(zhǔn)與測量技術(shù)、觀測與系統(tǒng)集成技術(shù)等方面開展研究，建立了由敦煌戈壁、青海湖水體、思茅雨林組成的全球第一個多尺度、全譜段輻射校正場。發(fā)展了遙感衛(wèi)星輻射校正理論與方法，建立了高絕對精度、高空間分辨率和全譜段在軌輻射校正模型。突破了高精度輻射標(biāo)準(zhǔn)傳遞技術(shù)，實(shí)現(xiàn)了具有國際先進(jìn)水平的場地和大氣光譜輻射觀測。建立了國際先進(jìn)的輻射校正業(yè)務(wù)運(yùn)行系統(tǒng)，形成了輻射校正質(zhì)量保障體系[1-4]。

經(jīng)過近20年的運(yùn)行，中國遙感衛(wèi)星輻射校正系統(tǒng)針對我國各系列遙感衛(wèi)星在軌輻射校正需求，開展了近50多次大規(guī)模在軌遙感衛(wèi)星輻射校正試驗，完成了30多顆國產(chǎn)衛(wèi)星和10余顆國外衛(wèi)星的輻射校正觀測，解決了我國遙感衛(wèi)星因缺少穩(wěn)定的星上定標(biāo)系統(tǒng)而產(chǎn)生的定量應(yīng)用問題，在我國現(xiàn)代衛(wèi)星遙感系統(tǒng)發(fā)展和應(yīng)用方面發(fā)揮了不可替代的重大作用。中國遙感衛(wèi)星輻射校正系統(tǒng)極大地推動了我國遙感衛(wèi)星定量化應(yīng)用水平的提高，確保了遙感衛(wèi)星效益的發(fā)揮，使我國在這一領(lǐng)域進(jìn)入國際先進(jìn)行列[5-10]。

中國遙感衛(wèi)星輻射校正場敦煌戈壁場區(qū)位于敦煌市以西30 km的戈壁灘上（緯度N40.00°~40.33°，經(jīng)度E94.00°~94.50°），自2002年通過國家驗收并轉(zhuǎn)入業(yè)務(wù)運(yùn)行以來，先后為我國所有對地觀測遙感衛(wèi)星提供了在軌輻射定標(biāo)支撐。作為國際對地觀測組織輻射定標(biāo)與真實(shí)性檢驗的有地面儀器支撐的重要場地之一，敦煌場同時為國外對地觀測衛(wèi)星的輻射定標(biāo)和真實(shí)性檢驗提供地面場地同步測量數(shù)據(jù)[11-15]，在國內(nèi)外遙感定量化應(yīng)用領(lǐng)域具有不可替代的作用。

2015年12月國家衛(wèi)星氣象中心完成了敦煌場自動化觀測基地建設(shè)，實(shí)現(xiàn)了敦煌場地表反射波段輻射特性參數(shù)（包括表面雙向反射分布函數(shù)（BRDF）、輻亮度）與大氣狀態(tài)參數(shù)（包括大氣輻照度、透過率、氣溶膠光學(xué)厚度、水汽總量、全天空云量、臭氧總量、溫度、濕度、壓強(qiáng)等）的全天候自動觀測。所有過境遙感衛(wèi)星反射波段均可利用自動觀測數(shù)據(jù)實(shí)現(xiàn)連續(xù)自動輻射定標(biāo)[16-18]。

敦煌場因其得天獨(dú)厚的表面均勻、氣候干燥等特性而被選作遙感輻射校正場，同時也具有豐富的太陽能、風(fēng)能資源。2016年12月，位于敦煌場外圍南側(cè)的10 MW熔鹽塔式光熱發(fā)電項目（小電站）建成并投入使用[19]。2018年12月100 MW熔鹽塔式光熱發(fā)電項目（大電站）投入使用[20]。

然而，光熱電站集熱塔的散射輻射對敦煌場區(qū)光環(huán)境變化帶來的影響需要開展深入細(xì)致的研究。本文著重研究集熱塔頂吸熱器的散射輻射，使用模型評估與實(shí)測分析相結(jié)合的方案對其影響進(jìn)行定量分析。

2 散射輻射引起的光環(huán)境變化分析

2.1 模型模擬

利用Monte Carlo三維輻射傳輸模擬方法[21]，模擬計算集熱塔頂吸熱器反射光經(jīng)大氣多次散射后被衛(wèi)星觀測到的輻射能量分布情況。

由于集熱塔頂吸熱器的反射光不同于太陽等自然平行光源，而是非平行光源（圖1）；同時，探求反射光對鄰近像元的影響本質(zhì)上是一個三維輻射傳輸問題，因此，大氣科學(xué)領(lǐng)域常用的許多針對太陽和地球輻射的平面平行輻射傳輸求解方法，并不適用于集熱塔反射光的輻射傳輸模擬研究。在非平行穩(wěn)態(tài)光源問題和三維輻射傳輸問題方面，Monte Carlo方法是一種常用的理論模擬工具。

圖 1 在自動觀測基地附近觀察到的工作狀態(tài)的集熱塔Fig. 1 Heat collection tower in working condition observed from the automatic observation base

圖 2 等高線測量角度示意圖Fig. 2 Schematic diagram of the almucantar measurement angle

Monte Carlo方法是一種隨機(jī)模擬方法，它并不求解輻射傳輸方程，而是對輻射傳輸?shù)膶?shí)際過程進(jìn)行直接模擬。在應(yīng)用Monte Carlo方法模擬集熱塔反射光的輻射傳輸過程時，模擬程序從光源處連續(xù)釋放具有初始能量權(quán)重的光子，并且在大氣介質(zhì)中逐個跟蹤這些光子的運(yùn)動路徑。由光源發(fā)出的光子在大氣介質(zhì)中將被散射和吸收。結(jié)合計算機(jī)產(chǎn)生的在適當(dāng)區(qū)間上均勻分布的隨機(jī)數(shù)以及構(gòu)建的與散射過程有關(guān)的概率分布函數(shù)，可以確定光子每一次散射后的行進(jìn)路程和方向。光子若被介質(zhì)吸收，Monte Carlo方法根據(jù)介質(zhì)的單次散射吸收比的數(shù)值來調(diào)整該光子的能量權(quán)重值，然后光子以新的能量權(quán)重在新的方向上繼續(xù)前進(jìn)。追蹤過程一直重復(fù)下去，直到所用的光子達(dá)到一定數(shù)目，符合預(yù)期的精度，則結(jié)束此次隨機(jī)模擬過程，并統(tǒng)計最終的大氣頂出射能量分布，其用出射的光子數(shù)表示。根據(jù)上述過程，可以模擬得到衛(wèi)星觀測到的由集熱塔光源發(fā)出并經(jīng)大氣多次散射后到達(dá)衛(wèi)星入瞳處的輻射能量分布情況，所需要的相關(guān)參數(shù)如表1定義。

表 1 模擬的相關(guān)參數(shù)定義Tab. 1 Definitions of parameters in the simulation

2.2 實(shí)測分析

集熱塔頂吸熱器的表面材料為白色，經(jīng)陽光匯聚后目視明亮。但與太陽相比，吸熱器所造成的漫射光較為微弱，如何將其與原有的天空漫射輻射相剝離是測量的關(guān)鍵。如果與早期集熱塔未建成，不存在光源的時候進(jìn)行對比，即使在相同的觀測幾何下，由于大氣狀態(tài)的差異，也很難保證天空漫射輻射不發(fā)生變化。

敦煌場區(qū)自身的漫射輻射近似滿足以太陽到觀測點(diǎn)連線為中心的旋轉(zhuǎn)對稱分布[22]。所以當(dāng)大氣足夠穩(wěn)定、晴空無云、氣溶膠與下墊面分布均勻時，可以用靠近塔頂吸熱器的漫射輻射減去遠(yuǎn)離吸熱器、相對于太陽對稱位置的漫射輻射，以二者的差值代表單純由吸熱器漫射輻射產(chǎn)生的影響。

上述方法可利用CE318等高線觀測模式獲取的數(shù)據(jù)進(jìn)行驗證。CE318可實(shí)現(xiàn)在不同太陽高度位置，360°旋轉(zhuǎn)測量等高線上的天空漫射輻射。這樣的觀測模式為實(shí)際測量吸熱器漫射輻射提供了可能。為了排除云的影響，需要使用云量自動觀測儀輻射選取晴空數(shù)據(jù)參與計算驗證。

2.2.1 觀測儀器

法國Cimel公司生產(chǎn)的CE318太陽光度計[23]相關(guān)參數(shù)如表2、表3所示。表2 中S代表對太陽觀測，A代表對太陽周邊光暈觀測，K代表對天空光觀測。

表 2 CE318通道設(shè)置Tab. 2 CE318 channel specifications

表 3 CE318性能參數(shù)Tab. 3 CE318 performance parameters

表3 說明CE318的動態(tài)范圍增益可調(diào)，量化等級達(dá)15位，觀測天空漫射輻射的輸出計數(shù)值在103量級，可以滿足定量分析吸熱器散射輻射的需求。

為了確保天空漫射光觀測數(shù)據(jù)中沒有云底輻射的干擾，使用安徽云能天智能科技有限責(zé)任公司生產(chǎn)的ASC200云量自動觀測儀對觀測時刻的云量進(jìn)行監(jiān)控。儀器相關(guān)參數(shù)如表4所示[24]。

表 4 ASC200性能參數(shù)Tab. 4 ASC200 performance parameters

CE318是目前國際上非常通用的一款自動太陽光度計，它采用濾光片分光的通道式設(shè)計。與光柵或棱鏡分光的高光譜儀器相比可獲得更可靠的觀測數(shù)據(jù)，但無法實(shí)現(xiàn)全波長范圍的觀測。由于現(xiàn)階段缺少整個太陽反射波段高光譜等高線觀測設(shè)備，使用半球天空漫射輻射高光譜數(shù)據(jù)描述天空漫射光的譜形特征。所使用的自動觀測儀器為由安徽光學(xué)精密機(jī)械研究所開發(fā)的漫射/總輻射比測量設(shè)備HIM。HIM的相關(guān)參數(shù)如表5所示。

表 5 HIM性能參數(shù)Tab. 5 HIM performance parameters

2.2.2 等高線測量數(shù)據(jù)篩選依據(jù)

利用CE318自帶的等高線測量模式[23]，在每天的固定時刻與大氣質(zhì)量數(shù)下，收集等高線測量數(shù)據(jù)。測量角度間隔如圖2(彩圖見期刊電子版)所示。圖2中Φ表示角度間隔，CE318等高線測量模式首先從太陽位置開始，保持高度角不變，方位角從0°掃描至+180°（almucantar right，ALR）。然后回到太陽位置，保持高度角不變，方位角從0°掃描至?180°（almucantar left，ALL）。其中方位角是太陽方位角與觀測方位角的差值。ALR方位角的變化間隔為0°S、3.0°A、3.5°A、4.0°A、5.0°A、6.0°A、6.0°K、7.0°K、8.0°K、10.0°K、12.0°K、14.0°K、16.0°K、18.0°K、20.0°K、25.0°K、30.0°K、35.0°K、40.0°K、45.0°K、50.0°K、60.0°K、70.0°K、80.0°K、90.0°K、100.0°K、120.0°K、140.0°K、160.0°K、180.0°K。ALL方位角變化間隔與其互為相反數(shù)。

觀測數(shù)據(jù)的質(zhì)量控制閾值如表6所示。為保證天空漫射輻射符合對稱分布，僅選擇云量為零的觀測數(shù)據(jù)參與計算。利用CE318的觀測數(shù)據(jù)計算氣溶膠光學(xué)厚度[25]，刪除光學(xué)厚度大于0.2的點(diǎn)。同時要求ALL與ALR的觀測時間差異在1 min以內(nèi)。當(dāng)距離集熱塔的投影大于5 km時，要求ALL與ALR觀測數(shù)據(jù)相對偏差的標(biāo)準(zhǔn)差要小于1%。這一限制不適用于距離小于等于5 km的數(shù)據(jù)，這是為了保留塔頂吸熱器輻射帶來的偏差。但仍需要通過人工判讀的方式剔除5 km內(nèi)與電站無關(guān)的明顯抖動數(shù)據(jù)。當(dāng)太陽與吸熱器的夾角過小時，與吸熱器對稱的太陽位置也可能會受到影響，所以要求太陽與吸熱器的地面投影距離大于2.5 km。同時要求大小電站位于太陽的同側(cè)，以避免大小電站彼此出現(xiàn)在對稱位置或離對方的對稱位置過近的情況。

表 6 觀測數(shù)據(jù)的質(zhì)量控制閾值Tab. 6 Quality control threshold of observation data

為了盡可能地利用有效數(shù)據(jù)，未設(shè)置電站在等高線上的投影角度與觀測節(jié)點(diǎn)角度的差異閾值。但是需要明確的是投影角度越接近節(jié)點(diǎn)角度，所得結(jié)果越接近真實(shí)情況。

2.2.3 四象限定位誤差的修正

由于存在四象限定位誤差，晴空時ALL與ALR測量得到的光暈輻射差異較大。如表3所示CE318四象限定位精度為0.1°，對應(yīng)光暈觀測值的相對差異達(dá)10%，這將嚴(yán)重影響對天空漫射輻射對稱性的判斷。

本文提出了對四象限定位誤差進(jìn)行修正的算法。在敦煌天空漫射輻射符合旋轉(zhuǎn)對稱分布的前提下，設(shè)計如下定位修正公式。

式（1）中x=[?6 ?5 ?4 ?3.5 ?3 3 3.5 4 5 6]，其為ALL與ALR觀測太陽光暈時的相對方位角，單位是“°”。之所以不使用更大的角度參與計算是因為電站相對方位角一般大于50°，如果將電站附近的數(shù)據(jù)納入回歸將減小電站可能存在的影響。v(x)是光暈?zāi)Ｊ较屡cx對應(yīng)的輸出計數(shù)值，v0a是定標(biāo)系數(shù)，單位是W·m?2·sr?1·μm?1。a1、a2、a3是最小二乘回歸系數(shù)，其中a2的數(shù)值代表了四象限定位誤差。

2.3 模型模擬與實(shí)測分析相結(jié)合的定量分析方法

模型模擬是以中心點(diǎn)（衛(wèi)星觀測天頂角為0°，距離塔頂為0 km）為基準(zhǔn)值1求取相對衰減分布，無法計算得到絕對的輻射量值。實(shí)測方法可以得到真實(shí)的漫射輻射量值，但僅在某幾個位置有實(shí)測數(shù)據(jù)的觀測幾何，無法遍及每個觀測天頂角與距離。特別是極端靠近或遠(yuǎn)離塔頂時無法獲取有效的實(shí)測數(shù)據(jù)。

為了實(shí)現(xiàn)對集熱塔散射輻射的定量分析，采用模型模擬與實(shí)測分析相結(jié)合的方法，即先借助模型模擬得到散射輻射相對于中心點(diǎn)的相對分布，再利用實(shí)測數(shù)據(jù)得到這個相對分布上某個點(diǎn)的實(shí)際散射輻射量值；繼而推算出其他位置上的實(shí)際量值。

3 結(jié)果與分析

3.1 集 熱 塔 頂 吸 熱 器 散 射 輻 射 模 型 模 擬 結(jié) 果

使用Monte Carlo方法開展輻射傳輸模擬計算，對集熱塔周圍大氣的影響開展定量分析。

圖3 ( 彩圖見期刊電子版)為衛(wèi)星觀測到的集熱塔反射光的輻射能量水平分布示意圖（衛(wèi)星觀測天頂角為0°，觀測方位為正北方向）。模擬結(jié)果如表7所示。

圖 3 衛(wèi)星觀測到的集熱塔反射光的輻射能量水平分布示例Fig. 3 An example of the radiant energy level distribution of the light reflected by the heat collection tower observed by the satellite

表 7 衛(wèi)星觀測到的集熱塔反射光的輻射能量相對分布情況Tab. 7 The relative distribution of the radiant energy of the light reflected by the heat collection tower observed by the satellite

結(jié)果顯示集熱塔對大氣輻亮度的貢獻(xiàn)與距離及衛(wèi)星觀測天頂角相關(guān)。設(shè)集熱塔所在像元在衛(wèi)星觀測天頂角為0°時的觀測輻射為基準(zhǔn)值1，距離集熱塔3 km時，垂直觀測時對大氣的影響為1.97%，觀測天頂角為50°時對大氣的影響可達(dá)2.69%。

3.2 對集熱塔頂吸熱器散射輻射的實(shí)測結(jié)果

圖4 為2020年3月2日協(xié)調(diào)世界時（UTC）7:00:05時刻ASC200云量自動觀測儀拍攝到的一幅云量為零時的紅外云圖。

圖 4 2020年3月2日UTC 7:00:05時刻ASC200云量自動觀測儀紅外云圖Fig. 4 Infrared cloud image obtained from ASC200 cloud cover automatic observer at 7:00:05 UTC on March 2, 2020

與傳統(tǒng)的可見光成像儀相比，紅外成像儀云圖太陽周圍沒有光暈，邊界清晰，可避免光暈被誤判為云的情況。圖4左下角邊緣處的兩個亮點(diǎn)中，稍亮的為工作中的大電站，稍暗的為待機(jī)中的小電站。

為了保證天空漫射輻射的絕對測量精度，于2019年12月11日至20日在中國氣象科學(xué)研究院頂層平臺對所使用的CE318進(jìn)行了標(biāo)定。采用基準(zhǔn)儀器MASTER傳遞溯源方式[26]，使用標(biāo)準(zhǔn)機(jī)#746對儀器進(jìn)行了定標(biāo)。定標(biāo)后，反演的光學(xué)厚度與MASTER的相對差異在2%以內(nèi)。

完成定標(biāo)后，將CE318安裝于敦煌自動觀測基地，開展了為期3個月的自動觀測。對2020年1月8日—3月31日收集到的CE318等高線觀測數(shù)據(jù)進(jìn)行質(zhì)量控制。共采集等高線觀測數(shù)據(jù)4 011組（觀測間隔小于1 min的ALL與ALR組合）。按照表6對觀測數(shù)據(jù)進(jìn)行篩選后，獲得有效數(shù)據(jù)26組，如表8所示。

表 8 各級篩選后有效觀測數(shù)據(jù)量值統(tǒng)計（組）Tab. 8 Statistics of effective observation data after being selected at different levels (group)

其中篩選項目分別為：

(a)觀測間隔小于1 min的ALL與ALR組合；(b)云 量<0；(c)光 學(xué) 厚 度<0.2；(d)ALL與ALR相對偏差的標(biāo)準(zhǔn)差<1%（與大電站的投影距離>5 km)；(e)ALL與ALR相 對 偏 差 的 標(biāo) 準(zhǔn) 差<1%（與小電站的投影距離>5 km)；(f)太陽與大電站的投影距離>2.5 km；(g)太陽與小電站的投影距離>2.5 km；(h)大小電站在太陽的同側(cè)。

CE318在1 020、1 640、870、675、440 nm和500 nm這6個通道開展等高線觀測。比較CE318測量的6個通道天空漫射輻射數(shù)據(jù)，發(fā)現(xiàn)波長越長，通過（d）和（e)篩選的有效數(shù)據(jù)組數(shù)越少。這由該波段漫射輻射抖動的劇烈程度與儀器的靈敏程度共同決定。1 640 nm有效數(shù)據(jù)多于1 020 nm的原因是1 640 nm采用的是銦鎵砷探測器，符合所設(shè)定的質(zhì)控標(biāo)準(zhǔn)且有效數(shù)據(jù)大于3組的通道只有440 nm和500 nm。因此針對這兩個通道展開具體的分析，統(tǒng)計結(jié)果如表9和表10所示。其中小電站相對于場地自動觀測基地距離2.86 km，方位角為110.48°。大電站相對于場地自動觀測基地距離3.83 km，方位角為149.98°。

表 9 吸熱器散射輻射帶來的相對偏差（440 nm）Tab. 9 Relative deviation due to scattered radiation from the heat absorber (440 nm)

表 10 吸熱器散射輻射帶來的相對偏差（500 nm）Tab. 10 Relative deviation due to scattered radiation from the heat absorber (500 nm)

在表9和表10中，AOD表示氣溶膠光學(xué)厚度，Sz與Sa分別是太陽天頂角與方位角，Beta=a_2是四象限定位修正角。HB與HS分別代表等高線處于大小電站上方的距離，on表示電站的工作情況，1為正常工作，0為待機(jī)。DegB與DegS表示太陽方位角與大小電站方位角的差值。DdegB與DdegS表示大小電站與等高線觀測節(jié)點(diǎn)角度差。RdevB與RdevS表示四象限定位修正前大小電站與太陽軸對稱位置處的輻亮度的相對偏差，RdevBr與RdevSr則表示四象限定位修正后大小電站相對偏差。

圖5 給出表9與表10中有效數(shù)據(jù)對應(yīng)時刻的全天空漫射輻射相對光譜分布圖，并與CE318不同通道光譜響應(yīng)函數(shù)進(jìn)行了對比。由于HIM在2020年2月27日出現(xiàn)故障，所以缺少這天的數(shù)據(jù)。從圖中可以發(fā)現(xiàn)在440 nm和500 nm這兩個通道天空散射輻射能量的分布是最強(qiáng)的，相對而言，信噪比高，受大氣波動影響小。這樣可以從譜形分布角度說明僅440 nm和500 nm通道滿足篩選要求是合理的。

圖 5 歸一化的全天空漫射輻射光譜與CE318不同通道光譜響應(yīng)函數(shù)Fig. 5 Normalized full sky diffuse irradiance and CE318 spectral response function

太陽的天頂角越大，等高線距離塔頂?shù)木嚯x也就越近，但很難獲得等高線離吸熱器距離非常近的點(diǎn)。已知大電站集熱塔高度為260 m，太陽在這個高度時已經(jīng)接近傍晚，吸熱器的亮度有限，所造成的散射輻射與天空漫射輻射都淹沒在觀測噪聲中無法測量。在表9和表10中，觀測距離范圍是0.65~3.55 km?？梢园l(fā)現(xiàn)滿足質(zhì)量控制的數(shù)據(jù)中，小電站正常工作的數(shù)據(jù)較少，僅有4條。下面僅以大電站正常工作時的數(shù)據(jù)進(jìn)行分析。

疊加吸熱器輻射的天空漫射輻射理論上應(yīng)該大于沒有疊加的那一邊。但是440 nm和500 nm通道的RdevBr有為負(fù)的情況，最小值分別為?2.07%（03?09 09:59）和?3.65%（03?03 08:03）。負(fù)值的存在說明雖然制定了嚴(yán)格的質(zhì)量控制標(biāo)準(zhǔn)，但所確定的有效觀測中，等高線上的天空漫射輻射還是存在小范圍的波動。具體比較上述兩個通道相對偏差RdevBr最小的情況如圖6與圖7(彩圖見期刊電子版)所示。

圖 6 440 nm通道RdevBr為最小值（?2.07%）時對應(yīng)時刻（03?09 09:59）的輻亮度(a)與相對偏差(b)Fig. 6 Radiance (a) and relative deviation (b) at the corresponding time (03?09 09:59) when the 440nm channel RdevBr is the minimum (?2.07%)

圖6 與圖7中圖例RdB、RdBr、RdS、RdSr、Rd、Rdr分別代表RdevB、RdevBr、RdevS、RdevSr和修正前后電站所在半等高線與對稱半等高線天空漫射輻射的相對偏差Rdev和Rdevr。Std是1 min內(nèi)ALL與ALR觀測天空漫射輻射相對偏差的標(biāo)準(zhǔn)差。StdB與StdS分別是距離大、小電站投影距離大于5 km的ALL與ALR的標(biāo)準(zhǔn)差。圖6的Rdr明顯存在一個獨(dú)立于電站影響的波動，圖7的Rdr波動則更為明顯。電站對漫射輻射的增益無法抵消波動帶來的減小，因此對應(yīng)的RdevBr為負(fù)值。為了保留吸熱器輻射帶來的偏差，未對方差過大的距電站投影距離≤5 km的數(shù)據(jù)進(jìn)行篩除，因此需要進(jìn)一步通過人工判讀的方式剔除5 km內(nèi)與電站無關(guān)的明顯抖動數(shù)據(jù)。據(jù)此人工剔 除440 nm 03?09 09:59時 刻，500 nm 03?03 08:03、02?27 07:03時刻的數(shù)據(jù)。

圖 7 500 nm通道RdevBr為最小值（?3.65%）時對應(yīng)時刻（03?03 08:03）的輻亮度(a)與相對偏差(b)Fig. 7 Radiance (a) and relative deviation (b) at the corresponding time (03?03 08:03) when the 500 nm channel RdevBr is the minimum (?3.65%)

對篩選后的數(shù)據(jù)進(jìn)行統(tǒng)計對比，440 nm大電站工作時的RdevBr均值為?0.57%，方差為0.51%，小電站待機(jī)時的RdevSr均值為?0.59%，方差為0.82%。說明大電站附近的天空漫射輻射平均值略小于太陽對稱位置的輻射值，疊加吸熱器的輻射也未能改變。待機(jī)狀態(tài)小電站的RdevSr反映的是自然情況下天空漫射輻射的差異。小電站附近天空漫射輻射與對稱位置輻射的差異以0.82%的方差略有波動。在440 nm工作狀態(tài)大電站未對天空漫射輻射帶來明顯影響。

500 nm大電站工作時的RdevBr均值為0.27%，方差為0.67%。小電站待機(jī)時的RdevSr均值為?0.15%，方差為0.66%。說明500 nm工作狀態(tài)大電站會對天空漫射輻射帶來微小影響。對應(yīng)的RdevBr最大值為0.93%（01?25 09:04），如圖8（彩圖見期刊電子版）所示。說明在有效數(shù)據(jù)對應(yīng)的觀測幾何下（距離HB 0.87~3.07 km，觀測天頂角Sz 77.30°~51.32°）吸熱器對天空漫射輻射的影響不超過0.93%。

圖 8 500 nm通道RdevBr為最大值（0.93%）時對應(yīng)時刻（01?25 09:04）的輻亮度(a)與相對偏差(b)Fig. 8 Radiance (a) and relative deviation (b) at the corresponding time (01?25 09:04) when the 500 nm channel RdevBr is the maximum (0.93%)

3.3 模型模擬與實(shí)測分析相結(jié)合分析

表10 說明RdevBr隨距離的變化規(guī)律不明顯。為了進(jìn)一步分析大電站吸熱器散射輻射的影響，下面與模型模擬的結(jié)果結(jié)合起來分析。

以500 nm為例，實(shí)測有效數(shù)據(jù)中觀測天頂角最小為51.32°，對應(yīng)與大電站的距離為3.07 km，最接近于表7中3 km、50°時的情況（0.027）。觀測天頂角最大為77.30°，與大電站的距離0.87 km，對表7數(shù)據(jù)開展二維插值外推，可計算得到相應(yīng)值為0.10。根據(jù)實(shí)測結(jié)果吸熱器對天空漫射輻射的影響不超過0.93%，與0.027等比例換算后可得到大電站漫射輻射與太陽對稱位置的相對偏差隨距離與觀測角度變化情況，如表11所示。與0.10等比例換算時表11的數(shù)值還要再減小3.7倍，保守起見以對漫射輻射影響最大的數(shù)據(jù)作為分析依據(jù)。

表 11 大電站散射輻射帶來的天空漫射輻射相對變化隨距離與觀測角度變化(%)Tab. 11 The relative change of the sky diffuse radiation caused by the scattered radiation from the large power station changes with distance and the observation angle (%)

4 結(jié) 論

文中定位修正對均勻性判斷結(jié)果的影響是明顯的。有效數(shù)據(jù)中四象限定位修正角平均值為0.11°，接近儀器官方標(biāo)稱精度（0.1°）。修正前后數(shù)據(jù)對的相對偏差在方位角小于50°時較為明顯。以3°光暈處電站所在半等高線與無干擾半等高線相對偏差為例，修正前的平均相對偏差達(dá)?9.18%，修正后僅為?1.34%。如果不使用修正，有效數(shù)據(jù)個數(shù)將從26組降低為5組。此方法對CE318天空漫射測量數(shù)據(jù)的準(zhǔn)確使用具有借鑒意義。

收集2020年1至3月滿足設(shè)定閾值的CE318數(shù)據(jù)顯示對于440 nm通道，工作狀態(tài)大電站未對天空漫射輻射帶來明顯影響。除了500 nm通道，在有效數(shù)據(jù)對應(yīng)的觀測幾何下（距離HB 0.87~3.07 km，觀測天頂角Sz 77.30°~51.32°）吸熱器對天空漫射輻射的影響不超過0.93%。與模型模擬的結(jié)果結(jié)合起來分析后，可以得出當(dāng)距離電站2 km時大電站散射輻射帶來的天空漫射輻射相對變化<2%，大于等于3 km時相對變化<1%的結(jié)論。

綜上，中國遙感衛(wèi)星輻射場敦煌場區(qū)周邊的光電項目建設(shè)和運(yùn)行對場區(qū)光環(huán)境變化影響極為有限，不會對已有的場區(qū)地表、大氣觀測設(shè)施和觀測數(shù)據(jù)造成影響。

致謝：感謝敦煌市氣象局和敦煌光電產(chǎn)業(yè)園的各位同仁在收集相關(guān)信息資料時給予的大力幫助。