摘 要：針對(duì)中國西部地區(qū)輻射資源充沛但觀測(cè)資料匱乏的特點(diǎn)，提出一種基于輻照度觀測(cè)數(shù)據(jù)、遙感數(shù)據(jù)、 McClear和隨機(jī)森林算法的太陽輻照度超短期預(yù)測(cè)方法，并重點(diǎn)分析遙感數(shù)據(jù)對(duì)輻照度預(yù)測(cè)效果的影響。結(jié)果表明：添加遙感數(shù)據(jù)能夠優(yōu)化不同時(shí)間步長的輻照度預(yù)測(cè)效果，并能顯著降低平均絕對(duì)百分比誤差（MAPE）值高于40%的預(yù)測(cè)大誤差出現(xiàn)概率。同時(shí)，添加遙感數(shù)據(jù)對(duì)預(yù)測(cè)效果的提升隨時(shí)間步長呈線性增加關(guān)系，nRMSE的差值變化范圍從2.08%變?yōu)?3.81%；nMAE的差值從1.64%變化為14.52%；R2的差值隨時(shí)間步長的變化最為明顯，從-0.03變?yōu)?0.43。但值得注意的是，添加衛(wèi)星數(shù)據(jù)會(huì)顯著增加模型的建立和超參尋優(yōu)時(shí)間。

關(guān)鍵詞：太陽輻照度；預(yù)測(cè)；機(jī)器學(xué)習(xí)；FY-4A；晴空模型

中圖分類號(hào)：P40 文獻(xiàn)標(biāo)志碼：A

0 引 言

在中國，超過75%的電力仍由煤炭供應(yīng)，這意味著現(xiàn)在仍然向環(huán)境中排放著大量的NOx、CO2和SO2。太陽能是為實(shí)現(xiàn)到2030年中國化石能源使用率為20%這一目標(biāo)的重要戰(zhàn)略技術(shù)［1-3］。近年來，作為應(yīng)用光伏發(fā)電系統(tǒng)最重要的光伏組件和相關(guān)設(shè)備的成本均大幅下降，這為未來大規(guī)模應(yīng)用光伏發(fā)電并將太陽能作為世界最主要能源之一創(chuàng)造了條件［4-7］。

太陽輻照度預(yù)測(cè)能夠?yàn)楣夥到y(tǒng)的順利運(yùn)行提供技術(shù)支持。研究表明［1］，根據(jù)數(shù)據(jù)類型和方法，太陽輻照度預(yù)測(cè)方法可簡(jiǎn)單分為4種：1）數(shù)據(jù)驅(qū)動(dòng)方法［8-9］；2）基于全天空成像儀的方法［10-11］；3）基于衛(wèi)星的方法［12-13］；4）數(shù)值天氣預(yù)報(bào)（numerical weather prediction，NWP）［14］。數(shù)據(jù)驅(qū)動(dòng)方法適用于所有預(yù)測(cè)時(shí)間步長，但通常適用于時(shí)間步長大于1 h的預(yù)測(cè)。雖然基于全天空成像儀和衛(wèi)星的方法通常可利用連續(xù)圖像進(jìn)行預(yù)測(cè)，但基于全天空成像儀的預(yù)測(cè)時(shí)間步長通常在1 h以內(nèi)，而基于衛(wèi)星的預(yù)測(cè)時(shí)間步長可從0.5 h延長到6 h。衛(wèi)星數(shù)據(jù)或數(shù)值模式由于其網(wǎng)格大小和輸出時(shí)間的原因，在時(shí)間和空間分辨率上受到一定限制。同時(shí)，數(shù)據(jù)驅(qū)動(dòng)預(yù)測(cè)方法需要來自一級(jí)（氣象）站的多數(shù)據(jù)輸入，有潛在的校準(zhǔn)和數(shù)據(jù)質(zhì)量控制問題 ［15］。因此，如何以低成本提供高精度空間和時(shí)間分辨率的太陽輻照度預(yù)報(bào)成為一個(gè)重要研究課題。對(duì)于輻照度預(yù)測(cè)，研究方法的選擇也至關(guān)重要。現(xiàn)有大多數(shù)方法不能解決噪聲環(huán)境中的復(fù)雜非線性關(guān)系［16］。因此，隨著人工智能的發(fā)展，越來越多的研究人員選擇將機(jī)器學(xué)習(xí)方法應(yīng)用于太陽輻照度預(yù)測(cè)。Deo等［16］開發(fā)了支持向量機(jī)小波耦合模型（WSVM），采用日照時(shí)數(shù)、溫度、風(fēng)速、蒸發(fā)量和降水量作為模型輸入?yún)?shù)，該模型可輸出可靠的日預(yù)測(cè)數(shù)值；Yagli等［17］結(jié)合機(jī)器學(xué)習(xí)方法，通過一種稱為核條件密度估計(jì)的新技術(shù)糾正了偏差，結(jié)果表明，與使用地面數(shù)據(jù)生成的輻照度預(yù)測(cè)相比，使用偏差校正衛(wèi)星衍生數(shù)據(jù)生成的預(yù)報(bào)具有更高的精度。

上述研究證明了機(jī)器學(xué)習(xí)在輻照度預(yù)測(cè)中的可行性。然而，如何利用有限的觀測(cè)數(shù)據(jù)建立西北地區(qū)高時(shí)空分辨率的輻照度預(yù)測(cè)模型仍需進(jìn)一步研究。因此，本研究基于FY-4A衛(wèi)星的高分辨率圖像信息，結(jié)合輻照度觀測(cè)、機(jī)器學(xué)習(xí)方法和晴空模型，構(gòu)建超短時(shí)（0～3 h）太陽輻照度預(yù)測(cè)模型。并分別利用敏感性實(shí)驗(yàn)，測(cè)試在輸入數(shù)據(jù)中包含（或去除）遙感數(shù)據(jù)對(duì)不同時(shí)間步長輻照度預(yù)測(cè)效果的影響，并將二者的預(yù)測(cè)效果與持久性模型進(jìn)行比較。最終通過相關(guān)指標(biāo)，分析遙感數(shù)據(jù)在高精度太陽輻照度超短期預(yù)測(cè)中的作用。

1 數(shù)據(jù)和方法

1.1 觀測(cè)數(shù)據(jù)

地面太陽輻照度數(shù)據(jù)來自中國西北的3個(gè)氣象觀測(cè)站，即榆中、民勤和敦煌。西北地區(qū)植被類型簡(jiǎn)單、人口稀少、氣候干旱，占中國總裝機(jī)容量和并網(wǎng)電力的60%。這些臺(tái)站提供了使用FS-S6太陽輻射傳感器測(cè)量的短波向下全球輻照度（W/m2）數(shù)據(jù)。儀器每兩年校準(zhǔn)一次，符合ISO 9060：1990和WMO標(biāo)準(zhǔn)的一流技術(shù)指標(biāo)。所使用數(shù)據(jù)集為2019年7月—2020年6月，時(shí)間分辨率為1 min。表1列出了3個(gè)觀測(cè)點(diǎn)的地理信息（包括經(jīng)度、緯度和海拔高度）、平均和最大總輻照度（GHI）以及數(shù)據(jù)量。

1.2 風(fēng)云四號(hào)衛(wèi)星數(shù)據(jù)（FY-4A）

2016年，中國的新一代靜止氣象衛(wèi)星風(fēng)云四號(hào)A成功發(fā)射，其采用SAST5000平臺(tái)和六角圓柱結(jié)構(gòu)，具有地面面積大、質(zhì)心低的優(yōu)點(diǎn)。風(fēng)云四號(hào)A衛(wèi)星裝載了先進(jìn)的觀測(cè)儀器，包括先進(jìn)的干涉式大氣垂直探測(cè)器、空間環(huán)境監(jiān)測(cè)儀器、地球靜止軌道輻射成像儀和閃電成像儀。風(fēng)云衛(wèi)星數(shù)據(jù)已證明了其在中國區(qū)域進(jìn)行輻照度預(yù)報(bào)的可行性［18］。

風(fēng)云四號(hào)A衛(wèi)星云圖數(shù)據(jù)來自國家氣象科學(xué)數(shù)據(jù)中心（http：//satellite.nsmc.org.cn/），本研究使用的衛(wèi)星圖像的空間分辨率為2 km×2 km，并選擇7個(gè)通道（0.45～4.00 μm）。首先，刪除太陽高度角小于10°的衛(wèi)星圖像（由于光線較弱，觀測(cè)誤差較大）。其次，對(duì)圖像進(jìn)行幾何定標(biāo)和輻射定標(biāo)。圖像的裁剪面積為32 km×32 km，反演區(qū)域的像素位于中間。最后，通過線性插值對(duì)缺失的圖像進(jìn)行插值，獲得7個(gè)通道的10 min時(shí)間分辨率（與觀測(cè)數(shù)據(jù)的分辨率相匹配）的區(qū)域反照率平均值。

1.3 McClear數(shù)據(jù)

McClear是一種新的晴空模型，它可利用MACC數(shù)據(jù)源來估計(jì)地面接收的下行短波直接輻射和全球輻射。近年來，許多研究證明了McCle-ar的可靠性（包括在中國）［19-20］。其中，通過輸入經(jīng)度、緯度、海拔和輸出格式，可直接獲得1 min、15 min、1 h、日和月時(shí)間尺度的全球總輻射、直接輻射和散射輻射 （https：//www.soda-pro.com/web-services/radiation/cams-mcclear）。本研究下載了2019年6月—2020年7月，榆中、民勤和敦煌時(shí)間分辨率為1 min的太陽總輻照度，作為機(jī)器學(xué)習(xí)輸入?yún)?shù)之一。

1.4 隨機(jī)森林（RF）

隨機(jī)森林（random forest，RF）方法是一種回歸樹（regression tree，RT），它是一種改進(jìn)的套袋回歸樹模型。在20世紀(jì)90年代，Breiman［21］建議通過樹形結(jié)構(gòu)并對(duì)其預(yù)測(cè)值進(jìn)行平均，以產(chǎn)生更穩(wěn)定的最終預(yù)測(cè)。由于隨機(jī)森林的適用范圍廣、精度高、難以過度擬合且能夠處理非線性數(shù)據(jù)，因而其廣泛應(yīng)用于不同領(lǐng)域的預(yù)測(cè)和分類問題中。文獻(xiàn)［22］介紹了RF的詳細(xì)理論。

在機(jī)器學(xué)習(xí)過程中進(jìn)行超參數(shù)尋優(yōu)，能夠快速找到特定的超參數(shù)數(shù)值，進(jìn)而實(shí)現(xiàn)模式在驗(yàn)證數(shù)據(jù)集上的最優(yōu)表現(xiàn)。在本研究中，選擇網(wǎng)格搜索法來進(jìn)行超參尋優(yōu)。對(duì)于RF，超參數(shù)優(yōu)化主要針對(duì)以下6個(gè)超參數(shù)：決策樹的數(shù)量（n_估計(jì)器）、決策樹的最大深度（max_depth）、分離樣本的最小數(shù)量（min_samples_split）、分離特征的最大數(shù)目（max_features）、葉節(jié)點(diǎn)樣本的最小數(shù)目（mi-n_ssamples_leaf）和是否進(jìn)行隨機(jī)抽樣（bootstrap）。

1.5 參考模型

為確保預(yù)測(cè)模型的可比性，引入楊大智［24］提出的參考模型。根據(jù)研究，單值氣候?qū)W和持久性是“數(shù)據(jù)聚合譜”的兩端。其中，持久性算法只利用最近的觀測(cè)值。通過計(jì)算滯后自相關(guān)，楊大智［24］的方法可將經(jīng)典的氣候?qū)W和持久性方法最佳地結(jié)合起來（在最小二乘意義上）。參考模型的核心參數(shù)如式（1）～式（4）所示。

[RRMSEc=1ni=1n（μ-xi）2=σ] （1）

[α=cov（xi-h，xi）σ2=γ（h）] （2）

[RRMSEcp=1ni=1n[αxi-h+（1-α）μ-xi]2=α2σ2+σ2-2αcov（xi-h，xi）] （3）

[yi=αxi-h+（1-α）μ] （4）

式中：[xi-h]——滯后自相關(guān)（假設(shè)末端效應(yīng)可忽略不計(jì)）；[μ]和[σ]——樣本均值和標(biāo)準(zhǔn)差（假設(shè)[n]較大）；[yi]——預(yù)測(cè)值。

1.6 評(píng)價(jià)方法

本研究選取4個(gè)統(tǒng)計(jì)指標(biāo)，包括歸一化均方根誤差（nRMSE）、歸一化平均絕對(duì)誤差（nMAE）、歸一化均偏誤差（nMBE）和相關(guān)系數(shù)（R2）來評(píng)估模型的性能［24］。

首先，對(duì)數(shù)據(jù)進(jìn)行預(yù)處理。數(shù)據(jù)被隨機(jī)分為訓(xùn)練數(shù)據(jù)集（70%）和測(cè)試數(shù)據(jù)集（30%）。接下來，將RF算法引入模型構(gòu)建，并引入隨機(jī)網(wǎng)格搜索方法來進(jìn)行超參數(shù)尋優(yōu)。此外，根據(jù)相關(guān)指標(biāo)（nRMSE、nMAE、nMBE、Skillscore、MAPE和R2），進(jìn)行預(yù)測(cè)效果評(píng)價(jià)。最后，分析得出遙感數(shù)據(jù)對(duì)不同時(shí)間步長輻射預(yù)測(cè)效果的影響，具體流程如圖1所示。

2 結(jié)果分析

為分析遙感數(shù)據(jù)添加與否對(duì)不同時(shí)間步長輻照度預(yù)測(cè)效果的影響，本研究中分別對(duì)添加遙感數(shù)據(jù)和不添加遙感數(shù)據(jù)的太陽輻照度超短期（10～180 min）預(yù)測(cè)效果進(jìn)行比較，預(yù)測(cè)結(jié)果如圖2所示。云的移動(dòng)和形成消散過程會(huì)影響地表輻射。遙感數(shù)據(jù)中的反照率信息能反映云的形成和消散過程。研究結(jié)果表明，在未添加風(fēng)云衛(wèi)星遙感數(shù)據(jù)時(shí)，機(jī)器學(xué)習(xí)方法和參考模型（CP）互有優(yōu)劣，但當(dāng)添加風(fēng)云衛(wèi)星遙感數(shù)據(jù)后，機(jī)器學(xué)習(xí)具備顯著優(yōu)勢(shì)。因此，遙感數(shù)據(jù)對(duì)輻照度預(yù)測(cè)效果的影響值得關(guān)注。

總體來說，添加衛(wèi)星遙感數(shù)據(jù)后不同時(shí)間步長的輻照度預(yù)測(cè)效果都比參考模型和未添加衛(wèi)星遙感數(shù)據(jù)時(shí)好。同時(shí)結(jié)合表2可知，添加FY-4A衛(wèi)星遙感數(shù)據(jù)對(duì)輻照度預(yù)測(cè)效果的影響隨著時(shí)間步長逐漸增加。其中，nRMSE的差值（無遙感數(shù)據(jù)減去添加遙感數(shù)據(jù)）變化范圍從10 min的2.08%變化為180 min的13.81%；nMAE的差值從10 min的1.64%變化為180 min的14.52%；[R2]的差值隨時(shí)間步長的變化最為明顯，從10 min的-0.03變化為180 min的-0.43；nMBE的變化與其他變量不同，呈現(xiàn)波動(dòng)的變化趨勢(shì)，其變化范圍為[-0.43～]0.10個(gè)百分點(diǎn)。值得注意的是，從60 min時(shí)間步長開始，輻照度預(yù)測(cè)效果的差別逐漸增大。

為更準(zhǔn)確、深入評(píng)價(jià)遙感數(shù)據(jù)對(duì)不同步長輻照度預(yù)測(cè)效果的影響，本研究參考文獻(xiàn)［2］方法，計(jì)算MAPE（平均絕對(duì)百分比誤差）數(shù)值，并統(tǒng)計(jì)MAPE數(shù)值大于40%在預(yù)測(cè)結(jié)果中出現(xiàn)的概率，直觀的評(píng)估預(yù)測(cè)模型的有效性，其計(jì)算式為：

[SMAPE=i=1nxi-yixi×100%n] （5）

式中：[n]——樣本數(shù)量；[xi]——實(shí)際值；[yi]——預(yù)測(cè)值。

結(jié)合圖2和圖3可知，在增加衛(wèi)星遙感數(shù)據(jù)后，不僅能夠提升總體的預(yù)測(cè)效果（參考nRMSE，nMAE和R2的數(shù)值），而且在不同時(shí)間步長的輻照度預(yù)測(cè)中能降低MAPE高于40%的大誤差出現(xiàn)概率。并且，遙感數(shù)據(jù)對(duì)大誤差出現(xiàn)概率的降低幅度與時(shí)間步長呈現(xiàn)線性增加關(guān)系。其中，3個(gè)站點(diǎn)10 min時(shí)間步長輻照度預(yù)測(cè)大誤差出現(xiàn)概率的降低幅度變化范圍為1.56%～1.95%，而180 min時(shí)間步長的大誤差出現(xiàn)概率的降低幅度迅速增加至12.59%～16.39%。同時(shí)，綜合圖4可知，在添加遙感數(shù)據(jù)后，不同觀測(cè)站點(diǎn)輻照度預(yù)測(cè)的總體效果

都得到明顯提升，具體表現(xiàn)為nMAE和nRMSE均下降7%左右，而R2提高約0.1。然而，值得注意的是，增加輸入變量（遙感數(shù)據(jù)）會(huì)顯著增加RF算法進(jìn)行模式構(gòu)建和超參尋優(yōu)的時(shí)間，進(jìn)而影響預(yù)測(cè)結(jié)果的輸出時(shí)間。其中，在i7-10700 CPU@2.90 GHz、16 GB RAM的計(jì)算機(jī)上運(yùn)行模式（10～180 min）的時(shí)間差值（未輸入遙感數(shù)據(jù)減去輸入遙感數(shù)據(jù)）變化范圍為[-1958.14～-1245.72 s]，運(yùn)行時(shí)間的差值并未隨預(yù)測(cè)時(shí)間步長呈現(xiàn)線性變化。

3 結(jié) 論

本研究針對(duì)高分辨率太陽輻照度超短期預(yù)測(cè)中氣象數(shù)據(jù)獲取困難，預(yù)測(cè)精度偏低等特點(diǎn)，提出一種基于輻照度觀測(cè)、McClear，遙感數(shù)據(jù)和隨機(jī)森林算法的太陽輻照度超短期預(yù)測(cè)方法。該預(yù)測(cè)方法通過輸入輻照度觀測(cè)和遙感數(shù)據(jù)，利用隨機(jī)森林和超參數(shù)尋優(yōu)構(gòu)建預(yù)測(cè)模型。結(jié)果表明在所有輻照度監(jiān)測(cè)站點(diǎn)，添加遙感數(shù)據(jù)的預(yù)測(cè)效果與未添加遙感數(shù)據(jù)相比，預(yù)測(cè)效果均得到顯著改善。其中，在180 min時(shí)間步長的預(yù)測(cè)中，3個(gè)站點(diǎn)添加遙感數(shù)據(jù)比未添加遙感數(shù)據(jù)的nRMSE平均下降13.81個(gè)百分點(diǎn)，nMAE平均下降14.51個(gè)百分點(diǎn)，R2則平均提高了43.46%。同時(shí)，添加遙感數(shù)據(jù)也能降低在輻照度預(yù)測(cè)中MAPE高于40%的大誤差出現(xiàn)概率，且降低幅度與時(shí)間步長呈現(xiàn)線性增加趨勢(shì)。同時(shí)，不同觀測(cè)站點(diǎn)的輻照度預(yù)測(cè)總體效果在添加遙感數(shù)據(jù)后均得到顯著提升。但值得注意的是，添加衛(wèi)星數(shù)據(jù)會(huì)顯著增加模型的建立和超參尋優(yōu)時(shí)間。因此，高性能計(jì)算機(jī)是實(shí)現(xiàn)利用衛(wèi)星遙感數(shù)據(jù)進(jìn)行高精度輻照度預(yù)測(cè)的前提。本研究能夠利用有限的輻照度觀測(cè)提供可靠的太陽輻照度超短期預(yù)測(cè)結(jié)果，可為太陽能的利用和新能源發(fā)展提供新思路。

［參考文獻(xiàn)］

［1］ YANG D E， KLEISSL J， GUEYMARD C A， et al. History and trends in solar irradiance and PV power forecasting： a preliminary assessment and review using text mining［J］.Solar energy， 2018， 168： 60-101.

［2］ JIA D Y，YANG L N， LV T， et al. Evaluation of machine learning models for predicting daily global and diffuse solar radiation under different weather/pollution conditions［J］. Renewable energy， 2022， 187： 896-906.

［3］ 胡斯勒?qǐng)D， 施建成， 李明， 等. 基于衛(wèi)星數(shù)據(jù)的地表下行短波輻射估算： 方法、 進(jìn)展及問題［J］. 中國科學(xué)： 地球科學(xué)， 2020， 50（7）： 887-902.

HU S L T， SHI J C，LI M， et al. A review of the estimation of downward surface shortwave radiation based on satellite data：methods， progress and problems［J］. Scientia sinica eerrae， 2020， 50（7）： 887-902.

［4］ 蔣俊霞， 高曉清. 光伏系統(tǒng)氣候效應(yīng)及影響機(jī)理研究進(jìn)展［J］. 高原氣象， 2022， 41（4）： 953-962.

JIANG J X， GAO X Q. Research progress on climate effect and" influence" mechanism" of" photovoltaic" systems［J］. Plateau meteorology， 2022， 41（4）： 953-962.

［5］ MANUEL á J， GABRIEL R J， ANTONIO O S， et al. Hybrid techniques to predict solar radiation using support vector machine and search optimization algorithms：a review［J］. Applied sciences， 2021， 11（3）： 1044.

［6］ ZHANG L， MA W， ZHANG D. Stacked sparse autoencoder in PolSAR data classification using local spatial" "information［J］." IEEE" geoscience" and" remote sensing letters， 2016， 13（9）： 1359-1363.

［7］ HINTON G E， SALAKHUTDINOV R R. Reducing the dimensionality of data with neural networks［J］. Science，2006， 313（5786）： 504-507.

［8］ YAGLI G M， YANG D E， SRINIVASAN D. Automatic hourly solar forecasting using machine learning models［J］. Renewable and sustainable energy reviews， 2019， 105： 487-498.

［9］ MEGIA F A， KURTZ B， LEVIS A， et al. Cloud tomography applied to sky images： a virtual testbed［J］.Solar energy， 2018， 176： 287-300.

［10］ KUHN P， WILBERT S， PRAHL C， et al. Shadow camera system for the generation of solar irradiance maps［J］. Solar energy， 2017， 157： 157-170.

［11］ ANDRé M， PEREZ R， SOUBDHAN T， et al. Preliminary assessment of two spatio-temporal forecasting technics for hourly" " satellite-derived" " irradiance" " in" " a" " "complex meteorological context［J］. Solar energy， 2019， 177： 703-712.

［12］ HARTY T M， HOLMGREN W F， LORENZO A T， et al. Intra-hour cloud index forecasting with data assimilation［J］. Solar energy， 2019， 185： 270-282.

［13］ WU E， CLEMESHA R E S， KLEISSL J. Coastal Stratocumulus cloud edge forecasts［J］. Solar energy， 2018， 164： 355-369.

［14］ MARION B， SIMITH B. Photovoltaic system derived data for determining the solar resource and for modeling the performance of other photovoltaic systems［J］. Solar energy， 2017， 147： 349-357.

［15］ KISI O， PARMAR K S. Application of least square support vector machine and multivariate adaptive regression spline models in long term prediction of river water pollution［J］.Journal of hydrology， 2016， 534： 104-112.

［16］ DEO R C，WEN X H， Q F I. A wavelet-coupled support vector machine model for forecasting global incident solar radiation using limited meteorological dataset［J］. Applied energy， 2016， 168： 568-593.

［17］ YAGLI G M， YANG D Z， GANDHI O， et al. Can we justify producing univariate machine-learning forecasts with satellite-derived solar irradiance？［J］. Applied energy，2020， 259： 114122.

［18］ HUANG C L，SHI H R， GAO L， et al. Fengyun-4 geostationary" nbsp;satellite-based" "solar" "energy" "nowcasting system and its application in North China［J］. Advances in atmospheric sciences， 2022， 39（8）： 1316-1328.

［19］ YANG D Z. Choice of clear-sky model in solar forecasting［J］. Journal of renewable and sustainable energy， 2020， 12（2）： 026101.

［20］ JIA D Y， HUA J J， WANG L W， et al. Estimations of global horizontal irradiance and direct normal irradiance by using fengyun-4A satellite data in Northern China［J］.Remote sensing， 2021， 13（4）： 790.

［21］ BREIMAN L. Bagging predictors［J］. Machine learning，1996， 24（2）： 123-140.

［22］ FENG Y，GONG D Z， ZHANG Q W， et al. Evaluation of temperature-based machine learning and empirical models for" predicting" daily" global" solar" radiation［J］." Energy conversion and management， 2019， 198： 111780.

［23］ YANG D Z. Making reference solar forecasts with climatology， persistence， and their optimal convex combination［J］. Solar energy， 2019， 193： 981-985.

［24］ YANG D Z. A universal benchmarking method for probabilistic solar irradiance forecasting［J］. Solar energy，2019， 184： 410-416.

NOWCASTING PREDICTION OF SOLAR IRRADIANCE BASED ON

FY-4A AND MACHINE LEARNING

Jia Dongyu1，Li Kaiming1，Gao Xiaoqing2，Gao Yumeng3

（1. College of Urban Environment， Lanzhou City University， Lanzhou 730070， China；

2. Northwest Institute of Eco-Environment and Resources， CAS/Key Laboratory of Land Process and Climate Change in Cold and Arid Regions，

Chinese Academy of Sciences， Lanzhou 730000， China； 3. Unit 94754 of the People’s Liberation Army， Jiaxing 314000， China）

Abstract：In view of the characteristics of abundant radiation resources but lack of observation data in China， this study proposes a short-term solar irradiance forecasting method based on radiation observation data， remote sensing data， McClear， and random forest algorithm， and focuses on analyzing the impact of remote sensing data on radiation forecasting effectiveness. The results show that adding remote sensing data can optimize the forecasting effectiveness at different time horizons and significantly reduce the probability of large prediction errors with a mean absolute percentage error （MAPE） value exceeding 40%. Additionally， the improvement of the forecasting effectiveness with the addition of remote sensing data increases linearly with the time horizon. The difference range of normalized root mean square error （nRMSE） changes from 2.08% to 13.81%， the difference of normalized mean absolute error （nMAE） changes from 1.64% to 14.52%， the difference of R2 shows the most significant change with the time step， changing from -0.03 to -0.43. However， it is worth noting that adding satellite data will significantly increase the time required for model establishment and hyperparameter optimization.

Keywords：solar irradiance； forecasting； machine learning； FY-4A； clear sky model