方 勇,李光達,李紹震
(深圳中廣核工程設計有限公司,深圳 518100)
隨著中國碳達峰、碳中和目標的確定,到2030 年中國風電、太陽能發電的總裝機容量將達到12 億kW 以上[1],光伏行業將迎來巨大發展機遇。光伏行業發展日趨成熟,以光伏發電系統精細化設計推動光伏發電成本下降成為市場的必然選擇,而選擇最優容配比是光伏發電系統精細化設計的重要手段之一,可為推動光伏電站全生命周期收益提升和平準化度電成本(LCOE)下降做出重要貢獻。
光伏發電系統的容配比,即光伏發電系統直流側光伏組件安裝容量與交流側逆變器額定輸出功率的比值。早期,光伏發電系統通常采用1:1 的容配比進行設計,在太陽輻照度較低的區域,受溫度及功率損耗等因素的影響,逆變器長期不能滿載運行,一定程度上造成了其輸出功率的浪費[2]。因此,在光伏發電系統優化設計中,計算與分析最優容配比成為國內外眾多學者的研究熱點之一。在國外,Velasco 等[3]提出了采用MATLAB 軟件仿真計算最優容配比的方法,Mounetou 等[4]提出使用PVsyst 仿真軟件深入研究影響容配比的因素。在國內,很多學者基于軟件或經驗計算最優容配比,并分析超配后的技術經濟特性,比如:牟娟等[5]根據Meteonorm 軟件數據庫,以系統投入產出比最低為尋優目標,計算分析最優容配比;馬慶虎等[6]通過分析不同太陽能資源區采用不同容配比對光伏電站發電量的影響,得到了不同地區的建議容配比;劉家鼎[7]通過實時監測數據,分析了在不同擴容工程的工程造價影響下,容配比與LCOE 的敏感性關系;謝磊等[8]分析了日照資源及LCOE 計算中各參數變化對LCOE 和最優容配比的影響,并給出了平單軸跟蹤光伏發電系統的容配比;陸煒等[9]通過分析各種情況下雙面組件光伏電站的內部收益率,得到了雙面組件光伏電站在不同背景反射率下典型的最優容配比;陳建國等[10]運用PVsyst仿真軟件,分析了不同容配比下的限光損失、光伏組件逐年衰減因素及限電因素對容配比設計的影響等,通過累積成本與發電量比值的最小化來確定最優容配比;王淑娟等[11]通過計算不同容配比下光伏組件輸出功率經折減后到達逆變器后的輸出功率,分析了不同容配比下的功率限制比,并進一步對增加部分投資的收益率進行了分析;王磊等[12]通過建立光伏逆變器失效概率的計算方法來確定最優容配比。
雖然很多學者對光伏發電系統的容配比進行了分析,但尚無通過曲線擬合方法計算得到最優容配比的相關論述。本文首先以光資源實測數據的分析及處理為基礎,通過解析計算及函數擬合曲線獲得光伏組件的理論最佳傾角,并將實測的水平面太陽輻照數據修正到傾斜面,從而計算光伏發電系統的發電量及棄光率;然后通過分析計算,獲得光伏發電系統的棄光率KQ、年均發電量棄光電量yQ等與容配比r的函數擬合關系,并以LCOE 為尋優目標,綜合考慮棄光率、光伏發電系統衰減率、直流側損耗、項目增值稅抵扣、運營成本等因素,定量分析計算出LCOE 與容配比之間的函數關系,進而確定最優容配比。
以某光伏發電項目為例進行分析,通過搜集該項目短期內的實測光資源數據及某數據庫中該項目地多年的光資源數據,統計分析并修正數據,將短期實測數據修正為水平年數據,以修正后的水平年實測數據作為光伏發電系統發電量計算的基礎。
某數據庫中項目所在地近15 年的水平面太陽輻照量數據如圖1 所示。
從圖1 可以看出: 1、2、11、12 月的太陽輻照量值較小,3~8 月的太陽輻照量值較大,9~10 月太陽輻照量值逐漸降低,基本符合項目所在地的太陽輻照量實際分布規律。

圖1 某數據庫中項目所在地近15 年的水平面太陽輻照量數據Fig. 1 Horizontal solar irradiation data of project location in a certain database over the past 15 years
通過統計分析某數據庫及短期實測的水平面太陽輻照量數據,并對短期實測的數據進行修正,修正后的項目所在地水平年實測水平面的太陽輻照量數據分布曲線如圖2 所示。

圖2 修正后的項目所在地水平年實測水平面的太陽輻照量數據分布曲線Fig. 2 Revised distribution curve of solar irradiation data measured on the horizontal plane in the level year of project location
從圖2 可以發現:某數據庫中2021 年的水平面年太陽輻照量為1460 kWh/m2,多年平均水平面年太陽輻照量為1498 kWh/m2;2021 年實測的水平面年太陽輻照量為1518 kWh/m2。實測值大于相應時段的數據庫值,說明數據庫值相對保守,通過計算分析,修正實測得到的1518 kWh/m2到水平年的水平面年太陽輻照量,修正后的水平年實測水平面年太陽輻照量為1563.6 kWh/m2。
結合該項目實際情況,計算得到不同光伏組件傾角(35°~50°)下光伏組件接收的年太陽輻照量數據[13],具體如圖3 所示。圖中:Ht為年太陽輻照量;S為光伏組件傾角;R2為決定系數。

圖3 不同光伏組件傾角下光伏組件接收的年太陽輻照量數據Fig. 3 Annual solar irradiation data received by PV modules at different tilt angles
由圖3 可知:R2=0.997,這說明該擬合曲線具有非常好的相關性。利用擬合函數進一步計算,可得到本項目的光伏組件理論最佳傾角為40.8°,此時傾斜面接收的年太陽輻照量為1923.8 kWh/m2。
本項目中光伏發電系統的基礎方案邊界條件如表1 所示。

表1 本項目中光伏發電系統的基礎方案邊界條件Table 1 Boundary conditions for basic scheme of PV power generation system in this project
光伏發電系統的實際發電量ym和棄光率可通過式(1)~式(7)計算得到。
式中:t為計算的起始時間;m為計算的截止時間,折算為數據最小采集間隔(1 min)的時間,如1 天、1 年、全生命周期等;P2為交流側的實時發電功率,MW;P0為棄光電量折算成的電功率或限發的功率,MW;ht為傾斜面的實時太陽輻照度,W/m2;Pdc為直流側的裝機容量,MW;Pac為交流側的額定容量,MW;KG為逆變器長期過載能力倍數;η為光伏發電系統的系統效率;P1為直流側實時發電功率,MW;ηdc為光伏發電系統直流側效率,不考慮逆變器之后的線路及設備損耗。
以本項目所在地某天傾斜面接收的太陽輻照度數據(如圖4 所示)為例,計算不同容配比下光伏發電系統直流側、交流側的實時發電功率及限發功率情況。

圖4 項目所在地某天傾斜面接收的太陽輻照度數據Fig. 4 Date of solar irradiance received by inclined plane of project location on a certain day
從圖4 可以看出:該項目中光伏發電系統單日傾斜面接收的太陽輻照度在10:00~14:00 之間的大部分時間均超過1000 W/m2,部分時間由于環境條件的變化,會低于1000 W/m2。
通過計算得到容配比為1:1 時,某天光伏發電系統直流側及交流側的發電功率與限發功率情況,如圖5 所示。

圖5 某天光伏發電系統直流側及交流側的發電功率與限發功率情況Fig. 5 Power and limiting power of DC and AC sides of PV power generation system on a certain day
結合圖4、圖5 中的數據,經過計算得到光伏發電系統單日限發總電量為0.5 MWh,單日實際發電量為665.86 MWh。
假設光伏發電系統交流側的額定功率保持100 MW 不變,分別按容配比1.00:1~1.40:1(以0.05 為間隔)計算某天光伏發電系統交流側發電功率及限發功率,結果如圖6~圖7 所示。

圖6 不同容配比時,光伏發電系統的交流側發電功率Fig. 6 AC side power of PV power generation system with different PV power to inverter power ratios

圖7 不同容配比時,光伏發電系統的交流側限發功率Fig. 7 AC side limited power of PV power generation system with different PV power to inverter power ratios
由圖6、圖7 可看出:當容配比較小(1.00:1~1.10:1)時,光伏發電系統交流側發電功率基本在限制值以下,限發功率很小,幾乎為零。當太陽輻照條件較好時,限發功率較大,且隨著容配比的增大,限發功率也隨之增大,當容配比為1.40:1時,限發功率達到峰值,為53 MW。
進一步計算實際年發電量及棄光電量,再根據光伏發電系統的衰減特性計算出不同容配比情況下的棄光率、全生命周期(25 年)發電量、年均發電量、棄光電量等參數。不同容配比時,光伏發電系統25 年內的年發電量如圖8 所示。

圖8 不同容配比時,光伏發電系統25 年內的年發電量Fig. 8 Annual power generation capacity of PV power generation system within 25 years with different PV power to inverter power ratios
從圖8 可以看出:隨著運營年限的增多,光伏發電系統的年發電量逐年遞減,且容配比越大,年發電量越高,但由于棄光電量隨著容配比的增大而增大,導致年發電量的增幅越來越小。
棄光率、年均發電量、棄光電量等可分別與容配比通過二次函數進行擬合,擬合結果如圖9~圖11 所示。
在圖9~圖11 擬合結果的基礎上可進一步求出其他容配比時擬合公式的相關參數,公式如式(8)~式(10)所示。不同容配比時得到的棄光率、年均發電量、棄光電量擬合公式的系數a、b、c及決定系數如表2 所示。

圖9 不同容配比下棄光率的擬合曲線Fig. 9 Fitting curves of ungrid-connected rate under different PV power to inverter power ratios

圖11 不同容配比下棄光電量的擬合曲線Fig. 11 Fitting curves of ungrid-connected power generation capacity under different PV power to inverter power ratios

圖10 不同容配比下年均發電量的擬合曲線Fig. 10 Fitting curves of average annual power generation capacity under different PV power to inverter power ratios
從表2 可看出:棄光率、年均發電量、棄光電量與容配比的二次函數擬合相關性非常高,說明該方法同樣適用于其他光伏發電項目不同容配比下的棄光率等因素擬合公式的參數計算,根據不同項目的實際情況,確定不同的擬合公式系數a、b、c即可。

表2 不同容配比時得到的棄光率、年均發電量、棄光電量擬合公式的參數a、b、c 及決定系數Table 2 Parameters a,b,c and coefficients of determination of fitting formula for ungrid-connected rate,average annual power generation capacity,and ungrid-connected power generation capacity under different PV power to inverter power ratios
LCOE 的計算式[14]如式(11)所示。
式中:i為折現率;n為光伏發電系統運行年數(n=1,2,…,N);N為光伏發電系統的全生命周期,年;I0為光伏發電系統的靜態初始投資,元;It為項目增值稅抵扣,元;VR為光伏發電系統殘值,元;Mn為第n年的運營成本(含維修、保險、材料、人工、輔助服務費等,不含利息),元;Yn為第n年的年上網電量,kWh。
經過計算,得到光伏發電系統LCOE 隨容配比的變化情況,如圖12 所示。

圖12 光伏發電系統LCOE 隨容配比的變化情況Fig. 12 Changes in LCOE of PV power generation system with PV power to inverter power ratios
從圖12 可以看出:在容配比為1.00:1 時,LCOE 最高;隨著容配比增大,LCOE 逐漸減小,在容配比約為1.20:1 時,LCOE 最低;隨著容配比進一步增大,LCOE 有所回升。由于本文的容配比間隔為0.05,從圖中趨勢可判斷最優容配比在1.15:1~1.25:1 之間,但無法通過圖獲得具體數值,為進一步確定最優容配比,可通過曲線擬合進行求解。擬合曲線公式如式(12)所示。
通過二次函數擬合,可計算得到最優容配比為1.228:1,對應的LCOE 為0.2574 元/kWh。從擬合結果來看,二次函數擬合決定系數為0.9946,非常接近1,且擬合的偏差(如表3 所示)非常小,基本在0.10%以下,說明二次函數擬合相關性很好,擬合結果可信。

表3 二次函數擬合值的偏差Table 3 Deviation of fitting value of quadratic function
本文通過對光資源數據分析及處理,采用解析方法計算光伏發電系統的棄光率及發電量,創新性地通過函數曲線擬合方法得到棄光率、年均發電量、棄光電量等與容配比的關系,并準確計算出任意容配比下的棄光率、年均發電量、棄光電量等參數。在此基礎上,分析計算全生命周期內的LCOE,得到LCOE 與容配比的函數關系曲線。結果顯示:通過二次函數擬合計算得到最優容配比的偏差基本在0.1%以下,決定系數大于0.99,擬合結果可信。由此可知,光伏發電系統的LCOE 與容配比之間一般存在強相關聯系,可通過曲線擬合方法求解得到最小LCOE 及最優容配比,該方法不僅適用于本項目,同樣可適用于不同地區、不同光資源的光伏發電系統。