沙區光伏陣列對近地層風沙輸移的干擾效應

唐國棟，蒙仲舉，高 永，黨曉宏，史芮嘉

沙區光伏陣列對近地層風沙輸移的干擾效應

唐國棟1,2，蒙仲舉3※，高 永3，黨曉宏3，史芮嘉4

（1. 中國水利水電科學研究院內蒙古陰山北麓荒漠草原生態水文野外科學觀測研究站，北京 100038； 2. 水利部牧區水利科學研究所，呼和浩特 010020； 3. 內蒙古農業大學沙漠治理學院/內蒙古自治區風沙物理與防沙治沙工程重點實驗室，呼和浩特 010011； 4. 鄂爾多斯市水土保持工作站，鄂爾多斯 010300）

沙漠地區建設光伏陣列后，地表吹蝕和堆積過程引起的地貌變化不僅嚴重威脅到了光伏組件固定結構的穩定，而且間接加速損耗了電板發電功率。為探究沙區光伏陣列擾動下近地層風沙輸移特征，在庫布齊沙漠中段的200 MW光伏電站腹地區域，通過同步測定光伏陣列腹地電板不同部位（板間、板前和板后）和上風向無光伏設施覆蓋的流動沙地近地層輸沙率，同時利用HOBO小型移動氣象站記錄觀測期風速和風向信息，分析不同風速風向條件下光伏陣列整體阻沙率、局部不同部位風沙流結構及通量模型。結果顯示：光伏陣列與風向夾角在–12.30°～82.19°范圍內，光伏陣列阻沙率為35.34%～93.02%，當夾角超過45°時，光伏陣列平均阻沙率可達84.63%；隨光伏陣列與風向之間夾角增大，板間和板后位置風沙輸移高度有向上層移動的趨勢，而板前位置則更加貼近地表；雙參數指數函數可以較好地模擬光伏陣列內不同部位近地層30 cm高度范圍內輸沙率隨高度的變化規律。研究結果有助于認識沙漠地區建設光伏陣列后近地層風沙輸移規律，可為科學制定次生風沙危害防治技術方案提供依據與參考。

太陽能；沙粒；光伏陣列；輸沙率；風沙流結構；擬合模型；庫布齊沙漠

0 引 言

太陽能光伏發電是解決傳統化石能源日益減少與經濟高速發展對能源需求日益增長之間矛盾的可靠途徑之一[1-3]。據國家能源局報道，2020年中國太陽能發電量達到2.53億kW·h，同比增長24.10%。依據安裝方式，太陽能光伏電板一般可分為地面式光伏和屋頂式光伏[4]。屋頂式光伏裝機量十分有限，目前中國主要以地面式光伏為主，占到總裝機量的70%[5]。中國西北干旱半干旱荒漠地區由于太陽能資源豐富、土地占用成本低，成為了規?；孛媸焦夥娬窘ㄔO的理想場所[6]。

在沙漠地區建設太陽能光伏電站雖然有著諸多好處，但同時也并存著一些挑戰。沙漠地區原本氣候干旱、風大沙多、水資源短缺、植被稀少，是生態環境最為脆弱地區[7]，隨著沙漠地區大規模太陽能光伏電站興建，施工對地表造成擾動，土壤活化為風沙活動提供了豐富的沙源[8]，而光伏設施擾動下流場格局發生變異，打破了原有維持地表平衡狀態的動力分布，地表風沙運動規律隨之改變[9]。研究表明，庫布齊沙漠110 MW光伏陣列風沙活動導致北側邊緣區域光伏電板板下出現掏蝕現象，形成以光伏電板下沿為軸線的風蝕坑（溝）[10]，在電板背風側板間區域形成堆積沙壟。地表形態變化一方面導致按照平坦沙表面風荷載設計的結構強度不足，嚴重風蝕區光伏組件有倒塌風險，另一方面地表形態變化勢必會加劇地表粉塵釋放速率，增加沙塵在光伏電板上的沉積，從而加速降低光伏電板發電效率[11]。有效防治沙區光伏電站地表次生風沙危害成為亟待解決的關鍵問題。

現階段關于沙區光伏陣列擾動下近地層輸沙通量相關研究較少，而且僅關注了單風向環境下的風沙運動規律[12-13]。Jubayer等[14-15]數值模擬和野外觀測研究結果證明，環境風向條件對光伏陣列近地層風速流場分布規律有重要影響，風是沙物質運移的動力因素，可見不同風向下光伏陣列內近地層風沙輸移特征也將存在差異。風向改變的實質是光伏陣列排布方向與風向之間的夾角（下文簡稱“夾角”）發生改變，不同夾角條件的氣流共同作用塑造了沙區光伏陣列地表形態。因此，研究不同夾角條件下沙區光伏陣列近地層風沙流結構，可以進一步深入探究風沙地貌形態的形成發育、演變發展等規律，進而采取客觀有效的措施控制或促進風沙運動，達到轉害為利的目的。

鑒于此，本研究以庫布齊沙漠200 MW光伏電站為研究對象，野外觀測不同夾角下光伏電板不同部位風輸沙通量，明確光伏陣列整體對近地層風沙輸移的干擾效應，闡明光伏陣列局部近地層風沙流固體流量結構、濃度分布和數學模型及其與夾角之間的關系，旨在揭示沙區建設太陽能光伏電站對近地層沙物質輸移的影響，為科學制定沙區光伏電站地表次生風沙危害防治技術方案提供理論依據。

1　材料與方法

1.1 研究區概況

研究區位于庫布齊沙漠中段，地屬內蒙古自治區鄂爾多斯市杭錦旗獨貴塔拉鎮，地理坐標為37°20′～39°50′N，107°10′～111°45′E。該區域屬于溫帶大陸性氣候，年平均氣溫為5～8 ℃，年降水量在258.30 mm左右，年潛在蒸發量2 400 mm左右，年太陽總輻射量597.90 kJ/cm，無霜期149 d。風沙活動主要集中時間在3－5月，全年大風日數為25～35 d，全年盛行西北風和西風，夏季盛行東南風。平均風速≥5 m/s的次數為323.40次，最多年份達418次。該地區的沙丘主要呈西北-東南走向的新月形沙丘、新月形沙丘鏈和格狀沙丘鏈等，沙丘高度10～60 m，其中較為平緩的沙地占10%左右，半固定沙地占30%，其余60%為流動性沙地，且流動性較強。

1.2　試驗設計與數據處理

1.2.1 試驗布設

研究光伏電站于2018年底完成安裝，整體占地面積約為5.20 km2，峰值發電量為200 MW。光伏電板面向正南，沿東西方向排布，南北相鄰兩排光伏電板間距800 cm，光伏電板傾角為36°，面板上沿距地面垂直高度270 cm，下沿距地面垂直高度35 cm。單組光伏電板由2排18列99 cm×195 cm基本光伏電板單元組成，單組光伏電板整體規格為400 cm×1 800 cm，地面投影寬度約為320 cm。

研究區域3－5月份平均風速較大，且降雨量稀少，相對濕度較低，導致該時期風沙活動最為活躍，因此野外試驗觀測選擇在2019年3月25日—4月27日期間開展，測試期間試驗光伏電板周圍地表植被蓋度為0。在野外試驗觀測前，相較于建設初期，試驗光伏電板周圍已經發生嚴重的吹蝕和堆積現象，光伏電板下掏蝕形成以光伏電板下沿為軸線的風蝕溝槽，板間形成堆積沙壟地貌。在觀測儀器布設前，首先將試驗光伏電板周圍地表進行平整，整平后的光伏電板下沿高度距離地表約35 cm。具體的儀器布設方法如圖1所示，在光伏電板板間、板前和板后3個位置同時放置集沙儀，同時在上風向無光伏電板覆蓋區域設置對照觀測。集沙儀采集高度為30 cm，共15個進沙口，進沙口規格為2 cm×2 cm。輸沙觀測依據風況每次觀測時間為20～60 min不等，風速越大觀測時間越短。所收集到的沙物質用自封袋分層取樣，帶回實驗室用0.01 g電子天平稱質量后，計算單位時間單位面積沙物質傳輸速率，即輸沙率，g/(cm2·min)。進而得出光伏電板不同位置各層（0～30 cm，每2 cm一層，共15層）輸沙率，累加得出總輸沙率。

與風沙輸移觀測同步進行，利用HOBO小型氣象站記錄對照觀測點2 m高度處風速和風向數據，風速風向的數據記錄間隔設定為1 s，數據采集間隔設置為3 s。為了探究風向改變對光伏陣列內近地層風沙活動規律的影響，將風沙輸移觀測期風向情況用光伏陣列與風向之間的夾角值表示。與光伏陣列平行風向W或E，夾角記為0°；與光伏陣列垂直風向為N，夾角記為+90°（夾角值的“±”沒有物理意義，僅用于區別方向，下同）；風向為S，夾角記為–90°。本研究共收集16組不同夾角條件下輸沙同步觀測數據（圖2），覆蓋了研究區域主風向與光伏陣列間的夾角范圍。

1.2.2 數據處理

1）阻沙率

阻沙率用于反映沙區光伏陣列內近地層水平風沙通量相較于流動沙丘對照的降低幅度，具體計算公式[16]如下：

式中q為光伏陣列阻沙率；0為流動沙丘上方輸沙率，g/(cm2·min)；PV為光伏陣列內輸沙率，g/(cm2·min)；TP、FP和HP分別為光伏陣列內板間、板前和板后同步觀測輸沙率，g/(cm2·min)。

2）輸沙通量擬合函數模型

目前對于水平輸沙通量分布函數仍然沒有定論，但大量的風洞模擬和野外實地觀測試驗結果顯示，指數函數[17-20]、冪函數[21-23]及其修正函數[24]能夠很好地模擬輸沙率隨著高度增加的衰減規律。因此，本研究采用以下4種模型對沙區光伏陣列近地層風沙流通量進行擬合分析：

式中為輸沙率，g/(cm2·min)；為集沙儀進沙口中心距離地表高度，cm；1、2、3、4、1、2、3、4和2均為模型擬合系數。

3）赤池信息量準則

風沙流通量最佳擬合模型采用赤池信息量準則（Akaike Information Criterion，AIC）確定，該準則是衡量統計模型擬合效果的一種標準，AIC值越小，表明模型擬合效果越好。計算公式為

注：為光伏陣列與環境風向之間的夾角。下同。圖中刻度線表示風速段對應的頻數，同心圓之間的間隔代表頻數差值（圖g、h為200，其余小圖均為100）。

Note:represents the intersection angle between solar photovoltaic array and wind direction. Same as below. The scale line in figures represents the frequency corresponding to the wind speed section, and the interval between concentric circles represents the frequency difference (200 in figures g and h, and 100 in other small figures).

圖2 無光伏覆蓋區的風向玫瑰圖

Fig.2 Wind direction rose diagrams of no photovoltaic coverage area

2 結果與分析

2.1 光伏陣列近地層輸沙率

不同風況條件光伏陣列局部典型部位和流動沙丘輸沙率如圖3所示，所有風況條件下流動沙丘近地層輸沙率總是高于光伏陣列內，在光伏陣列內的3個典型部位輸沙率大小關系隨夾角變化存在差異。如圖3a所示，當夾角為–12.30°時，不同典型位置輸沙率從大到小為：板前、板間、板后，板間和板后分別為板前輸沙的80.35%和27.77%。如圖3b～圖3d所示，當夾角為從9.13°～10.71°時，板間輸沙率增大，板前和板后位置輸沙率大致相當，平均為板間輸沙率的49.85%。當夾角為18.07°時（圖3e），板前輸沙率增大，輸沙率從大到小為：板間、板前、板后。隨著夾角的繼續增大，夾角在29.67°～61.61°范圍內時（圖3f～圖3p），板前輸沙率上升至光伏陣列內最大，且光伏陣列內不同部位輸沙率差異較小。

進一步探究光伏陣列對近地層沙物質輸移的擾動效應及其與環境風向之間的關系。運用光伏陣列內3個部位輸沙率均值來表征光伏陣列內近地層輸沙率，與上風向無光伏電板覆蓋的流動沙丘下墊面同步觀測的輸沙率數據進行對比分析，計算得到光伏陣列阻沙率如表1所示，可以看出不同風向條件下光伏陣列阻沙率有明顯差別，夾角由9.13°向18.07°變化時，阻沙率呈下降趨勢，由59.93%下降至35.34%。當夾角為29.67°時，阻沙率迅速上升至62.73%。隨著夾角的進一步增大，夾角在56.97°～82.19°范圍時，光伏陣列阻沙率為84.63%±6.11%（均值±標準差）。換言之，該風向條件下光伏陣列可降低近地層過境風沙流84.63%的輸沙量。此外，夾角為10.71°時阻沙率為43.84%，夾角為18.07°時阻沙率為35.34%，可知夾角為12.30°時的阻沙率值在35.34%～43.84%之間，而夾角為–12.30°時，阻沙率高達56.00%。由此推斷在夾角絕對值相等條件下，夾角為負值時的光伏陣列阻沙率較高。

表1 不同風況下光伏陣列阻沙率

2.2 光伏電板典型部位風沙流結構

為更加直觀地認識光伏電板不同位置沙粒濃度與高度之間的關系，將累積輸沙率隨高度的變化做成柱狀圖。由圖4可以看出，流動沙丘下墊面風沙流運移基本貼近地表，不同風況條件下90%以上的輸沙量均分布在0～10 cm高度范圍內，且隨著風速的增大，沙物質輸移的高度有所上升。

光伏陣列內不同部位風沙流垂直結構隨風況變化差異明顯。比較相似風速不同風向條件下的風沙流結構特征（圖4e、圖4f、圖4h、圖4l、圖4n風況下，風速在7～8 m/s范圍內，對應夾角差異較大，分別為18.07°、29.67°、61.61°、77.48°和79.39°），結果顯示夾角較小時（圖4e和圖4f），光伏陣列不同部位沙物質輸移高度較低且差異較小，75%以上的輸沙量在0～6 cm高度范圍內傳輸，90%以上的輸沙量在0～10 cm高度范圍內傳輸。夾角較大時（圖4h、圖4l和圖4n風況），板間和板后位置沙物質輸移高度增加，75%以上的輸沙量在0～10 cm高度范圍內傳輸，90%以上的輸沙量在0～20 cm高度范圍內傳輸。而板前位置沙物質輸移高度則略有降低，75%以上的輸沙量在0～4 cm高度范圍內傳輸，90%以上的輸沙量在0～6 cm高度范圍內傳輸。

比較不同風速和不同風向條件下的風沙流結構特征，在夾角較小且風速較大風況條件下（圖4a～圖4d，風速在10.35～11.57 m/s范圍內，夾角絕對值在9.13°～12.30°范圍內），75%以上的輸沙量在0～8 cm高度范圍內傳輸，90%以上的輸沙量在0～12 cm高度范圍內傳輸。與之相比，夾角較大且風速較小風況條件下（圖4m～圖4p，風速在7.21～9.75 m/s范圍內，夾角在78.45°～82.19°范圍內），板間和板后位置沙物質輸移的高度明顯上升，75%以上的輸沙量在0～11 cm高度范圍內傳輸，90%以上的輸沙量在0～20 cm高度范圍內傳輸，板前位置沙物質輸移高度則表現出小幅下降趨勢，75%以上的輸沙量在0～4 cm高度范圍內傳輸，90%以上的輸沙量在0～5 cm高度范圍內傳輸。一般來說，風速越大，沙物質輸移高度越高，然而通過上述比較發現，即使風速相對較小，仍然呈現出隨著夾角的增大，板間和板后位置沙物質輸移高度表現出較強的向上層移動的趨勢。

綜上所述，隨著光伏陣列與風向夾角的增大，板間和板后位置沙物質輸移高度表現出較強的向上層移動趨勢，而板前位置沙物質輸移則相對更加貼近地表。

2.3 光伏電板典型部位輸沙通量模型

為研究輸沙率隨高度的變化規律，利用Origin軟件做出輸沙率隨高度變化的散點圖，結果顯示輸沙率隨著高度增加均呈現下降趨勢，不同風況下大致呈現3種變化規律，如圖5所示，隨著夾角的增大，光伏陣列內不同部位各高度層輸沙率大小關系發生改變，且相較于流動沙丘，光伏陣列內近地層10 cm高度范圍內輸沙率隨夾角增大而下降明顯。

為探究沙區光伏陣列擾動下近地層最佳風沙通量模型，運用指數函數、冪函數及其修正函數4種模型，對4個觀測點在16種不同風況下的64組數據擬合共256次，其中14次擬合失敗。擬合的2、AIC值和RMSE統計結果如表2所示，指數函數及其修正函數均方根誤差較小、2較大和AIC值更小。從擬合相關度2來看，有60.94%的情況下雙參數指數模型（模型1）與三參數指數模型（模型2）擬合結果一致，有39.06%的情況下雙參數指數模型略低于三參數指數模型擬合結果，而且這種差異性非常小，一般表現在0.001量級上；從AIC值來看，有64.06%情況下雙參數指數函數擬合結果更佳，有32.81%情況下三參數指數函數擬合結果更佳；綜合得出，雙參數指數模型可以較好地模擬光伏陣列內3個典型部位近地層30 cm高度范圍內輸沙率隨高度的變化規律。

表2 4種函數模型輸沙率擬合結果對比

注：AIC為赤池信息量準則；(z)為高度處輸沙率，g·cm-2·min-1；為集沙儀進沙口中心距離地表高度，cm；1、2、3、4、1、2、3、4和2均為擬合系數。下同。

Note: AIC is Akaike Information Criterion;(z)is the sand transport rate at height of, g·cm-2·min-1;represents the arithmetic mean of the top and bottom of each chamber in the sampler, cm;1,2,3,4,1,2,3,4, and2are regression coefficients. Same as below.

對于指數函數而言，系數表征風沙流中沙粒濃度的最大值或蠕移輸沙量[26]，系數一般認為是反映沙粒濃度隨高度的衰減程度[27]，即遞減率，值越大衰減速度越慢，值越小衰減速度越快。如表3所示，流動沙丘、板間、板前和板后位置系數值范圍為0.278～1.633、0.009～0.945、0.026～0.729和0.004～0.687。流動沙丘處的值高于光伏陣列內不同部位，表明光伏陣列的存在使得近地表風沙流的塵粒濃度或輸沙率降低。光伏陣列內不同部位的值隨風況條件變化存在差異，當夾角為9.13°～10.17°時，值表現為光伏電板板間較高，板前和板后相當且較低。當夾角增大至18.07°時，板前位置值略高于板間位置增至光伏陣列內最大。隨夾角繼續增大，夾角在29.67°～82.19°范圍內時，板前位置值最大，板間和板后位置值相當，平均僅為板前位置的0.33倍。

不同風況條件下流動沙丘、板間、板前和板后位置系數值范圍為2.103～3.826、3.820～6.120、2.163～3.821和2.525～3.953。隨著夾角的增大，板間位置值有一定的增大趨勢，表明輸沙率隨高度衰減速度越慢，即沙物質傳輸向高層移動，板前位置則隨夾角的增大值有減小的趨勢，表明輸沙率隨高度衰減速度變快，即沙物質傳輸向貼地層方向移動。板后位置值隨著夾角變化無明顯的規律。

表3 風沙流通量系數a，b與夾角之間的關系

注：和為函數模型擬合系數。當=56.97°時，由于風速低且輸沙率太小而擬合失敗。

Note:andare the fitting coefficients of function models. When the angle is 56.97°, the fitting fails because the wind speed is low and the sand transport rate is too small.

3 討 論

目前，關于沙區光伏陣列擾動下近地層風沙通量模型的研究相對較少，陳曦等[12]在烏蘭布和沙漠東南緣光伏電站的研究結果顯示，電板下沿、上沿及板間位置輸沙率隨高度增加均符合多項式函數形式。楊世榮[13]在庫布齊沙漠中部光伏電站的研究結果顯示，不同風速條件下電板不同部位最佳模型存在差異，6.5 m/s風速條件板下和板前輸沙通量指數函數模型表現最佳，而板間則冪函數模型最佳；8.4 m/s時板下和板前表現為冪函數模型最佳，而板間則為多項式函數模型最佳。不同的研究者觀測時環境地形地貌和風況不同，研究光伏電板規格、安裝角度、高度、間距等條件均存在差異，這將導致得出不同的甚至是截然相反的結論。本研究基于野外觀測數據，采用指數函數和冪函數及其修正函數共4種模型進行擬合分析，運用擬合相關度2和赤池信息量準則綜合對擬合結果進行評價，結果顯示雙參數指數函數模擬效果較好。本研究結果與上述他人研究結果產生差異的原因，主要是由于采用集沙儀的高度差異所致。陳曦等人和楊世榮觀測時均采用25層50 cm高度集沙儀，本研究采用的是15層30 cm高度集沙儀。而本研究光伏電板下沿野外觀測時距離地表約35～40 cm，因此在板前位置布設較高的集沙儀，當來流方向與光伏電板布置方向呈一定角度時，高出電板下沿以上的集沙盒子必定收集不到沙物質，即光伏電板對風沙流垂直結構有阻斷影響，由此可見，光伏陣列地表不同高度范圍內不同部位沙塵通量模型可能存在差別。綜上所述，本研究結果可以證實，不同風向光伏陣列擾動下地表30 cm高度范圍內雙參數指數能夠很好模擬輸沙率隨高度的變化趨勢。

風沙運動是一種貼近地表的沙粒搬運現象。張正偲等[25]通過451次風沙流觀測數據研究得出，騰格里沙漠東南緣平坦沙地上0～20 cm高度范圍內輸沙量占總輸沙量的比例大于95%。陳新闖等[28]對烏蘭布和沙漠流動沙丘野外觀測得出，80%以上輸沙量集中于0～10 cm高度，90%以上的輸沙量集中于0～20 cm高度。本研究發現，對照平坦裸沙地90%以上的輸沙量均分布在0～10 cm高度范圍內，與上述研究結果基本一致。在光伏陣列的擾動下，風向變化對沙粒傳輸高度有影響，夾角較大時，板間和板后位置沙物質輸移高度增加，90%以上的輸沙量在0～20 cm高度范圍內傳輸。而板前位置沙物質輸移高度則有下降趨勢，90%以上的輸沙量在0～6 cm高度范圍內傳輸。這是由于過境氣流運動方向與光伏陣列排布方向呈一定角度時，傾斜的光伏電板發揮導流作用使得下沿出風口氣流運動方向更加貼近于地表，因此板前沙物質輸移高度降低。而板間和板后區域此時容易產生渦旋[14]，從而增加了垂直向上的升力，導致風沙流表現出較強的向上層移動趨勢。

在流動沙丘地區，任何植被或機械沙障的覆蓋都會減少地表裸露，從而降低土壤可蝕性，達到保護土壤不被風蝕的作用[29-31]。通過與無光伏電板覆蓋區域對比研究發現，夾角在–12.30°～82.19°范圍內，近地層水平輸沙通量降低幅度在35.34%～93.02%之間，隨夾角增大而增大。當夾角超過45°時，光伏陣列平均可降低近地層84.63%的水平輸沙通量。光伏陣列與風向夾角相同而來流方向不同時，即夾角絕對值相等符號不同，光伏陣列阻沙率存在差異。本研究結果顯示，夾角為負值時光伏陣列的阻沙率更強，這是由于光伏陣列本身結構特性所致，光伏電板面朝正南，北高南低呈36°角傾斜。對于過境氣流的作用機理與風沙危害防治中的下導風工程相似。已有下導風工程領域相關模擬研究[32-33]表明，當夾角為+90°時，傾角為60°的導板迎面來向氣流在距離下沿約7/10的位置產生分離點，也就是說導板層流的70%導向地面端出風口，而30%導向大氣端出風口。夾角為–90°時則正好相反，導板層流的30%導向地面端出風口。以此類比，可以很好地解釋夾角為負值時光伏陣列的阻沙率強于夾角為正值時的情況。

綜合來看，光伏陣列與風向夾角是影響陣列內近地層風沙輸移強度的重要因素，夾角較大時，陣列本身可以具有較強的截流阻沙效應，此時迎風側邊緣區是防護重點；夾角較小時，陣列截流阻沙作用有限，挾沙氣流更加容易進入陣列內，進而引起地表的形態變化。此時陣列邊緣區域不僅要做好地表固沙措施，同時應設置具有防風作用的高立式機械沙障或建設防護林等，以削弱氣流進入光伏陣列時的初始動能，從而降低光伏陣列內風沙活動強度。

4 結 論

本研究主要分析不同風況條件下，庫布齊沙漠200 MW光伏電站整體對近地層水平輸沙通量的攔截效能和腹地區域電板不同部位（板間、板前和板后）輸沙率、風沙流結構及通量模型。得到如下結論：

1）沙漠地區建設光伏陣列后對近地層風沙輸移有一定攔截作用，尤其當光伏陣列與風向夾角超過45°時，光伏陣列平均可降低近地層84.63%的水平輸沙通量。此外，夾角絕對值相等條件下，夾角為負值比正值時光伏陣列對風沙輸移的攔截效能更加明顯。

2）光伏陣列局部不同部位沙粒濃度與高度之間的關系表現為：隨夾角增大，板間和板后位置風沙輸移高度表現出較強的向上層移動趨勢，而板前位置則相對更加貼近地表。雙參數指數模型可以較好地模擬光伏陣列內3個典型部位近地層30 cm高度范圍內輸沙率隨高度的變化規律，擬合相關度2達0.923～1.000，RMSE僅為0.740×10-3～567.262×10-3g/(cm2·min)。

[1] 崔榮國，郭娟，程立海，等. 全球清潔能源發展現狀與趨勢分析[J]. 地球學報，2021，42(2)：179-186.

Cui Rongguo, Guo Juan, Cheng Lihai, et al.Status and trends analysis of global clean energies[J].Acta Geoscientica Sinica, 2021, 42(2): 179-186.(in Chinese with English abstract)

[2] 鄭艷婷，徐利剛. 發達國家推動綠色能源發展的歷程及啟示[J]. 資源科學，2012，34(10)：1855-1863.

Zheng Yanting, Xu Ligang. The process of developing green energy in developed countriesand its lessons[J]. Resources Science, 2012, 34(10): 1855-1863. (in Chinese with English abstract)

[3] 林智欽. 中國能源環境中長期發展戰略[J]. 中國軟科學，2013(12)：45-57.

Lin Zhiqin. Medium and long-term development strategy of china’s energy & environment[J]. China Soft Science, 2013(12): 45-57. (in Chinese with English abstract)

[4] Huang B, Li Z N, Zhao Z F, et al. Near-ground impurity-free wind and wind-driven sand of photovoltaic power stations in a desert area[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2018, 179: 483-502.

[5] 曾澤榮，李進，羅多，等.中國建筑分布式光伏利用現狀及未來趨勢分析[J]. 建設科技，2020(20)：12-16，20.

Zeng Zerong, Li Jin, Luo Duo, et al. Analysis of Current Situation and Future Trend of Distributed Photovoltaic Utilization in Chinese Buildings[J]. Construction Scirncr and Trchnology, 2020(20): 12-16, 20. (in Chinese with English abstract)

[6] 文小航，尚可政，王式功，等. 1961－2000年中國太陽輻射區域特征的初步研究[J]. 中國沙漠，2008，28(3)：554-561.

Wen Xiaohang, Shang Kezheng, Wang Shigong, et al. Primary study on regional characteristics of solar radiation in China during 1961-2000[J]. Journal of Desert Research, 2008, 28(3): 554-561. (in Chinese with English abstract)

[7] 王濤，趙哈林. 中國沙漠科學的五十年[J]. 中國沙漠，2005，25(2)：3-23.

Wang Tao, Zhao Halin. Fifty-yearhistoryofChinadesertscience[J]. Journal of Desert Research, 2005, 25(2): 3-23. (in Chinese with English abstract)

[8] 袁方，張振師，卜崇峰，等. 毛烏素沙地光伏電站項目區風速流場及風蝕防治措施[J]. 中國沙漠，2016，36(2)：287-294.

Yuan Fang, Zhang Zhenshi, Bu Chongfeng, et al. Wind speed flow field and wind erosion control measures at photovoltaic power plant project area in Mu Us Sandy Land[J]. Journal of Desert Research, 2016, 36(2): 287-294. (in Chinese with English abstract)

[9] Macpherson T, Nickling W G, Gillies J A, et al. Dust emissions from undisturbed and disturbed supply-limited desert surfaces[J]. Journal of Geophysical Research, 2008, 113(2): 1-16.

[10] 郭彩贇，韓致文，李愛敏，等. 庫布齊沙漠110MW光伏基地次生風沙危害的動力學機制[J]. 中國沙漠，2018，38(2)：225-232.

Guo Caiyun, Han Zhiwen, Li Aimin, et al. Dynamic mechanism research on the secondary blown sand disaster in the 110MW photovoltaic arrays of the Hobq Desert[J]. Journal of Desert Research, 2018, 38(2): 225-232. (in Chinese with English abstract)

[11] 孫家歡，王濤，薛嫻，等. 降塵對中國北方主要城市光伏電站發電量的影響[J]. 中國沙漠，2016，36(4)：932-942.

Sun Jiahuan, Wang Tao, Xue Xian, et al. The effect of dustfall on the power output of photovoltaic plants located in major cities in northearn China[J]. Journal of Desert Research, 2016, 36(4): 932-942. (in Chinese with English abstract)

[12] 陳曦，高永，翟波，等. 沙區光伏電場的風沙流輸移特征[J]. 干旱區研究，2019，36(3)：684-690.

Chen Xi, Gao Yong, Zhai Bo, et al. Transport of wind-drift sand in photovoltaic electric field in sand area[J]. Arid Zone Research, 2019, 36(3): 684-690. (in Chinese with English abstract)

[13] 楊世榮. 庫布齊沙漠光伏電站風沙運動及蝕積特征研究[D].呼和浩特：內蒙古農業大學，2019.

Ynag Shirong. Study on Wind-blown Movement and Erosion Characteristics of the Hobq Desert Photovoltaic Power Station[D]. Hohhot: Inner Mongolia Agricultural University, 2009. (in Chinese with English abstract)

[14] Jubayer C M, Hangan H. A numerical approach to the investigation of wind loading on an array of ground mounted Solar Photovoltaic(PV) panels[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2016, 153: 60-70

[15] Etyemezian V, Nikolich G, Gillies J A. Mean flow through utility scale solar facilities and preliminary insights on dust impacts[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2017, 162: 45-56.

[16] 張登山，吳汪洋，田麗慧，等. 青海湖沙地麥草方格沙障的蝕積效應與規格選取[J]. 地理科學，2014，34(5)：627-634.

Zhang Dengshan, Wu Wangyang, Tian Lihui, et al. Effects of Erosion and Deposition and Dimensions Selection of Straw-checkerboard Barriers in the Desert of Qinghai Lake[J]. Scientia Geographica Sinica, 2014, 34(5): 627-634. (in Chinese with English abstract)

[17] 劉芳，郝玉光，辛智鳴，等. 烏蘭布和沙漠東北緣地表風沙流結構特征[J]. 中國沙漠，2014，34(5)：1200-1207.

Liu Fang, Hao Yuguang, Xin Zhiming, et al. The surface aeolian-sand flow structure in the northeastern margin of the Ulanbuh Desert[J]. Journal of Desert Research, 2014, 34(5): 1200-1207. (in Chinese with English abstract)

[18] Ellis J T, Li B, Farrell E J, et al. Protocols for characterizing aeolian mass-flux profiles[J]. Aeolian Research, 2009, 1(1/2): 19-26.

[19] Zhang Z C, Dong Z B, Wu G X. Field observations of sand transport over the crest of a transverse dune in northwestern China Tengger Desert[J]. Soil and Tillage Research, 2017, 166: 67-75.

[20] Zhang Z C, Dong Z B. Characteristics of aeolian sediment transport over different land surfaces in northern China[J]. Soil and Tillage Research, 2014, 143: 106-115.

[21] Mertia R S, Santra P, Kandpal B K, et al. Mass–height profile and total mass transport of wind eroded aeolian sediments from rangelands of the Indian Thar Desert[J]. Aeolian Research, 2010, 2(2/3): 135-142.

[22] Dong Z B, Lyu P, Zhang Z C, et al. Aeolian transport in the field: A comparison of the effects of different surface treatments[J]. Journal of Geophysical Research: Atmospheres, 2012, 117: 9210.

[23] Lyu P, Dong Z B, Ma X M. Aeolian sand transport above three desert surfaces in northern China with different characteristics (shifting sand, straw checkerboard, and gravel): field observations[J]. Environmental Earth Sciences, 2016, 75(7): 577.

[24] Dong Z B, Lu J F, Man D Q, et al. Equations for the near-surface mass flux density profile of wind-blown sediments[J]. Earth Surface Processes and Landforms, 2011, 36(10): 1292-1299.

[25] 張正偲，董治寶. 騰格里沙漠東南部野外風沙流觀測[J]. 中國沙漠，2013，33(4)：973-980.

Zhang Zhengcai, Dong Zhibao. Field observation of aeolian sediment flux in the southeast tengger desert[J]. Journal of Desert Research, 2013, 33(4): 973-980. (in Chinese with English abstract)

[26] Dong Z B, Liu X P, Wang H T, et al. The flux profile of a blowing sand cloud: A wind tunnelinvestigation[J]. Geomorphology, 2002, 49: 219-230.

[27] 王洪濤，董治寶，張曉航. 風沙流中沙粒濃度分布的實驗研究[J]. 地球科學進展，2004，19(5)：732-735.

Wang Hongtao, Dong Zhibao, Zhang Xiaohang. A study of the particle concentration in the sand flux a-wind tunnel investigation[J]. Advances in Earth Science, 2004, 19(5): 732-735. (in Chinese with English abstract)

[28] 陳新闖，董智，李錦榮，等. 烏蘭布和沙漠不同下墊面冬季沙塵通量[J]. 中國沙漠，2016，36(6)：1527-1527.

Chen Xinchuang, Dong Zhi, Li Jinrong, et al. Winter sand-dust flux over different underlying surfaces in the Ulan Buh Desert[J]. Journal of Desert Research, 2016, 36(6): 1527-1532. (in Chinese with English abstract)

[29] Fang H, Wu X, Zou X Y, et al. An integrated simulation- assessment study for optimizing wind barrier design[J]. Agricultural and Forest Meteorology, 2018,263: 198-206.

[30] Wang T, Qu J J, Niu Q H. Comparative study of the shelter efficacy of straw checkerboard barriers and rocky checkerboard barriers in a wind tunnel[J]. Aeolian Research, 2020, 43: 100575.

[31] Sankey J B, Law D J, Breshears D D, et al. Employing lidar to detail vegetation canopy architecture for prediction of aeolian transport[J]. Geophysical Research Letters, 2013, 40(9): 1724-1728.

[32] 辛林桂，程建軍，王連，等. 基于Design Exploration方法對鐵路下導風工程關鍵設計參數的優化[J]. 中國沙漠，2018，38(6)：1193-1199.

Xin Lingui, Cheng Jianjun, Wang Lian, et al. Optimization of key design parameters of the aviation baffle engineering of railway based on design exploration method[J]. Journal of Desert Research, 2018, 38(6): 1193-1199. (in Chinese with English abstract)

[33] 陳柏羽，程建軍，辛林桂，等. 基于離散伴隨求解器的鐵路下導風工程外形優化研究[J]. 鐵道科學與工程學報，2019，16(8)：1923-1930.

Chen Boyu, Cheng Jianjun, Xin Lingui, et al. Shape optimization of blower fences of railway based on discrete adjoint solver[J]. Journal of Railway Science and Engineering, 2019, 16(8): 1923-1930. (in Chinese with English abstract)

Interference effect of solar photovoltaic array on near surface aeolian sand transport in sandy areas

Tang Guodong1,2, Meng Zhongju3※, Gao Yong3, Dang Xiaohong3, Shi Ruijia4

(1.,,100038,; 2.,010020,; 3./,,010011,; 4.,010300,)

Deserts are ideal places to develop ground-mounted large-scale solar photovoltaic (PV) power stations. However, it is evitable surface erosion that may occur after the construction of a solar PV power station, where solar energy production, operation, and maintenance depend mainly on geomorphological changes in sandy areas. This study aims to investigate the characteristics of wind-sand movement under the interference of solar PV array, thereby reducing the damage to solar energy. The study area was located in the middle part of the Hobq Desert in China. The observation field was 300 m from the west edge of the test solar PV power station that was built at the end of 2018. There were no any protective measures on the surface of the solar PV power station during the test, such as sand-binding plants or sand-barriers. Field observations were conducted from 20 March to 13 April 2019. The reason was that the aeolian sand activity was the strongest in the study area during from March to May, due to the frequent occurrence of strong wind, extended drought, and limited rain. Thus, wind erosion led to the formation of trenches in the immediate vicinity of the downwind of panels, and sand ripples between adjacent north-south panels. The surface surrounding the test PV panels was smoothened, while the underlying surface was flattened before the experimental instruments were arranged. The flatting operation allowed for the comparison of experimental data. Sediment transport was measured in different wind directions above shifting dunes at three observation sites around the PV panels, such as between, in front of, and behind the panels. Meanwhile, the wind speed and direction were recorded using a HOBO sensor at the observation sites of shifting dunes. The sediment transport data was also collected at sixteen wind regimes. Sand-fixation of solar photovoltaic array, aeolian-sand flow structure and fitted model around the PV panels were then analyzed under the different wind regimes. The results showed that the near-surface sand transport rate above shifting dunes was always larger than that in the solar PV array in all cases, where the intersection angle between the solar PV array and wind direction was a key parameter to dominate the sand inhibition rate of solar PV array. Specifically, the sand inhibition rate ranged from 35.34% to 93.02% at the angle range from -12.30° to 82.19°. The mean value of sand transport rate above the solar PV array reduced to 84.63%, compared with the shifting dunes, especially with the angle exceeding 45°. There was also no change in the sand transport rate model when applying the solar PV array, similar to the shifting dunes. A two-parameter exponential function was better fit for the measured profiles of flux density on the near-surface of solar PV array. Wind-sand flow between and behind the panels tended to evidently move towards a high layer with the angle increased, where the rising range was 8-10 cm, whereas, the saltation height at the observation site before the panels tended to move towards a low layer, where the decrease range was 4-5 cm. The finding can contribute to the understanding of the wind-sand movement characteristics under the interference of solar PV array, providing insightful ideas to plan better technical schemes against wind-sand hazards at solar PV power stations.

solar energy; sand; photovoltaic array; sand transport rate; aeolian-sand flow structure; fitting model; the Hobq Desert

唐國棟，蒙仲舉，高永，等. 沙區光伏陣列對近地層風沙輸移的干擾效應[J]. 農業工程學報，2021，37(13)：101-110.

10.11975/j.issn.1002-6819.2021.13.012 http://www.tcsae.org

Tang Guodong, Meng Zhongju, Gao Yong, et al. Interference effect of solar photovoltaic array on near surface aeolian sand transport in sandy areas[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(13): 101-110. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2021.13.012 http://www.tcsae.org

2021-04-12

2021-06-08

國家重點研發計劃項目（2018YFC0507101）；內蒙古自治區科技重大專項（zdzx2018058-3）

唐國棟，博士，研究方向為荒漠化防治。Email：18247158690@163.com

蒙仲舉，博士，教授，研究方向為荒漠化防治。Email：mengzhongju@126.com

10.11975/j.issn.1002-6819.2021.13.012

TM615

A

1002-6819(2021)-13-0101-10