浮島光伏電場對地表溫度空間分布特征的影響

伯英杰， 曾業(yè)隆， 李國慶， 曹興文， 姚清秀

(1.魯東大學(xué)資源與環(huán)境工程學(xué)院，煙臺 264025； 2.中國地質(zhì)大學(xué)(北京)海洋學(xué)院，北京 100083；3.中國科學(xué)院空天信息創(chuàng)新研究院遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，北京 100101； 4.中國科學(xué)院新疆生態(tài)與地理研究所荒漠與綠洲國家重點(diǎn)實(shí)驗(yàn)室,烏魯木齊 830011)

0 引言

與煤炭、石油、天然氣等化石能源不同，太陽能光伏發(fā)電技術(shù)具有清潔安全、永續(xù)可再生和取用方便的優(yōu)勢[1]，其產(chǎn)業(yè)在全球范圍的發(fā)展，有望改變自工業(yè)革命以來以化石能源為基礎(chǔ)的能源生產(chǎn)和消費(fèi)方式，能有效緩解能源危機(jī)、全球碳排放和環(huán)境污染[2-5]。近年來，中國光伏產(chǎn)業(yè)在多樣化政策激勵下發(fā)展迅速，在2013年和2015年我國相繼超越意大利和德國，成為全球光伏發(fā)電裝機(jī)量最多的國家[6-7]。據(jù)國家能源局統(tǒng)計(jì)，截至2020年上半年，中國大陸的光伏發(fā)電裝機(jī)容量達(dá)到2.16億kW，其中大型集中并網(wǎng)式電站裝機(jī)占69%以上。

然而，大型光伏電場工程的建設(shè)不可避免地對所在區(qū)域的土地、生態(tài)以及能量循環(huán)等產(chǎn)生特定的環(huán)境效應(yīng)[8-11]。大型光伏電場工程在建設(shè)時(shí)對環(huán)境的影響與其他工業(yè)制造過程類似，例如，平整土地、構(gòu)筑新結(jié)構(gòu)等工程活動直接作用于土壤、植被、地貌和土地利用，使其原有的形態(tài)發(fā)生變化，甚至可能會破壞原有環(huán)境要素的穩(wěn)定結(jié)構(gòu)，造成土地?cái)_動、植被破壞和水土流失等問題； 電磁干擾、固體廢棄物污染等也是光伏電場工程環(huán)境影響評估的重要內(nèi)容[9-12]。同時(shí)，作為光伏電場核心的太陽能光伏陣列在建成投產(chǎn)后，一方面吸收短波輻射將太陽能轉(zhuǎn)換為電能，另一方面向外輻射長波輻射改變了局部地表的能量流動方式，形成熱環(huán)境效應(yīng)[13]。即使熱環(huán)境效應(yīng)尚未納入常規(guī)的光伏電場建設(shè)的環(huán)境影響評價(jià)體系，但是該熱環(huán)境效應(yīng)會影響光伏系統(tǒng)的工作溫度，進(jìn)而影響光伏的發(fā)電效率，研究表明較低的光伏電池或模塊的工作溫度可提高發(fā)電效率[14]。另外，光伏電場產(chǎn)生的熱環(huán)境效應(yīng)改變了原有的地表輻射、反照率和蒸散發(fā)狀態(tài)，對不同下墊面的局部氣候和生態(tài)系統(tǒng)功能產(chǎn)生影響[13,15]。地表溫度(land surface temperature, LST)是研究地表熱環(huán)境的重要因素，在大規(guī)模、連續(xù)性的長時(shí)間序列研究具有獨(dú)特優(yōu)勢的遙感觀測法已成為主要的熱環(huán)境評價(jià)方法之一[16-17]。因此，采用遙感技術(shù)研究大型光伏電場工程建設(shè)的熱環(huán)境對光伏發(fā)電產(chǎn)業(yè)的環(huán)境影響評估和當(dāng)?shù)丨h(huán)境保護(hù)具有重要意義。

目前，光伏電場對局地?zé)岘h(huán)境的影響主要集中在干旱荒漠區(qū)。已有研究表明，光伏電場對以干旱荒漠地帶為典型的大規(guī)模建設(shè)區(qū)的地表熱環(huán)境產(chǎn)生影響： 野外觀測與遙感觀測均發(fā)現(xiàn)干旱荒漠區(qū)的光伏電場具有降溫效應(yīng)，主要發(fā)生在白天[18-19]，其影響范圍約為100～600 m[20-21]。近年來，隨著光伏產(chǎn)業(yè)規(guī)模的不斷擴(kuò)大，繼荒漠、草原和居民屋頂之后，閑置水面成為建設(shè)光伏電場的可利用環(huán)境之一，為我國新能源和清潔生產(chǎn)的發(fā)展探索了新的開發(fā)模式和領(lǐng)域[22]。然而，由于下墊面性質(zhì)與以往干旱區(qū)光伏電場不同，浮島光伏電場對建設(shè)區(qū)及其周圍環(huán)境的LST空間分布特征的影響尚不明晰，判斷光伏電場對局地溫度的影響模式、空間范圍和變化幅度均沒有進(jìn)行有效的分析。鑒于此，本研究以目前全球最大的浮島光伏電場之一的安徽省淮南市典型浮島光伏電場為研究區(qū)，為了明晰浮島光伏電場對周圍LST的影響,開展以下研究： ①基于Google Earth Engine(GEE)在線云平臺和Landsat8遙感影像，采用單通道算法提取典型浮島光伏電場及其鄰近區(qū)域的LST，建立LST數(shù)據(jù)集； ②利用LST數(shù)據(jù)集，判斷光伏電場對LST空間分布特征的影響模式； ③在此基礎(chǔ)上，明確光伏電場對LST空間分布特征影響的空間尺度和增/降溫能力。

1 研究區(qū)概況與數(shù)據(jù)源

1.1 研究區(qū)概況

研究區(qū)位于安徽省淮南市謝家集區(qū)，地處亞熱帶與暖溫帶的過渡帶，屬暖溫帶半濕潤季風(fēng)氣候區(qū)，夏季高溫多雨、冬季寒冷干燥； 全年日照時(shí)數(shù)為2 000 h左右，光照充足； 年太陽輻射量在4 800 MJ/m2附近，太陽能資源豐富。近年來，當(dāng)?shù)匾劳匈Y源和政策優(yōu)勢，在采煤沉陷區(qū)的閑置水面規(guī)劃和建設(shè)大規(guī)模的浮島電場。本研究的浮島光伏電場的位置如圖1所示，選址于采煤塌陷區(qū)的閑置水面，采取水上發(fā)電、水下養(yǎng)殖的模式，總裝機(jī)容量為20 MWp，是目前全球最大的浮島光伏電場之一[22]。該浮島光伏電場于2015年11月開工建設(shè)，2016年3月正式并網(wǎng)發(fā)電； 2015年11月—2016年10月為光伏電場的穩(wěn)定建設(shè)時(shí)段。一般來說，光伏電場對LST的影響范圍在100～600 m以內(nèi)[20-21]，因此，本研究大致以浮島光伏電場為中心，面積為25 km2的區(qū)域?yàn)檠芯糠秶瑒澏ü夥妶鲋車? km的緩沖區(qū)作為光伏電場的潛在影響區(qū)。根據(jù)浮島光伏電場工程的施工時(shí)間段，對照Landsat8和谷歌地球的歷史影像，對光伏電場1 km范圍內(nèi)未發(fā)生類型變化的典型地物進(jìn)行目視解譯，最終得到總面積約為0.258 km2的浮島光伏電場區(qū)以及水體、林地、耕地、裸地和建設(shè)用地。研究區(qū)如圖 1 所示。

圖1 研究區(qū)位置示意圖

1.2 數(shù)據(jù)源

1.2.1 遙感數(shù)據(jù)

為了盡可能多地獲取無云遙感影像，提高數(shù)據(jù)的連續(xù)性，本研究使用GEE在線云平臺獲取研究區(qū)2013年3月—2019年2月“行-列號”為“121-38”和“121-37”的Landsat8遙感影像數(shù)據(jù)集(USGS Landsat8 Surface Reflectance Tier 1,ImageCollection ID,Landsat/LC08/C01/T1_SR)，共292景。該數(shù)據(jù)集的空間分辨率為30 m，包含5個(gè)可見光與近紅外(VNIR)波段和2個(gè)短波紅外(SWIR)波段，這些數(shù)據(jù)經(jīng)過大氣校正處理后得到地表反射率數(shù)據(jù)； 還有2個(gè)已被處理為亮度溫度的熱紅外(TIR)波段和質(zhì)量信息波段(pixel_qa)。

1.2.2 氣候數(shù)據(jù)

地面臺站逐月氣溫來源于中國氣象數(shù)據(jù)網(wǎng)的地面氣候資料月值數(shù)據(jù)集(http: //data.cma.cn)。選取以浮島光伏建設(shè)區(qū)為中心60 km緩沖區(qū)范圍內(nèi)2個(gè)氣象臺站的逐月平均氣溫，用于后續(xù)的統(tǒng)計(jì)分析，即壽縣(距離光伏電場11 km)和蚌埠市(距離光伏電場58 km)。

2 技術(shù)路線與研究方法

2.1 技術(shù)路線

本文的技術(shù)路線如圖2所示。首先，利用GEE在線云平臺檢索2013年3月—2019年2月全時(shí)段的Landsat8遙感影像，采用單通道算法提取研究區(qū)無云時(shí)間序列LST和逐月LST，通過計(jì)算逐月LST與氣溫之差(LSTs-a)以減小氣溫波動對LST的影響； 然后，選取夏季和冬季的LSTs-a，計(jì)算它們在光伏電場建成前后的變化幅度以判斷影響LST空間分布的模式，再利用局部Moran’s I方法進(jìn)行空間自相關(guān)分析以明確光伏電場影響LST空間分布的范圍； 最后，在解譯提取浮島光伏電場周圍典型地物的基礎(chǔ)上，對LSTs-a進(jìn)行變化幅度提取，分析光伏電場的增/降溫能力，從而明確浮島光伏電場在不同建設(shè)階段對LST空間分布的影響程度。

圖2 技術(shù)路線

2.2 研究方法

2.2.1 基于單通道算法的無云LST提取

由于Landsat8熱紅外波段B11(11.50～12.51 μm)的雜散光噪聲使該波段LST反演的精度較低[23]，因此本研究選擇單通道算法(single-channel algorithms,SC)對B10(10.60～11.19 μm)進(jìn)行LST反演。在本研究中，使用Jimenez-Munoz等提出的經(jīng)典單通道方法估算LST[24-25]。該方法在保證反演精度的基礎(chǔ)上簡單易行，適用于Landsat數(shù)據(jù)長時(shí)間序列的LST反演，廣泛應(yīng)用于城市地表環(huán)境評估和人類活動影響評估[26-28]，其計(jì)算公式為：

(1)

(2)

式中:ε為地表發(fā)射率；Lsen為傳感器端輻射亮度(W·m-2·sr-1·μm-1)；Tsen為傳感器亮度溫度； B10和B11的參數(shù)bγ分別為： 1 324 K和1 199 K；ψ1，ψ2和ψ3為大氣函數(shù)?？衫孟率将@得，即

(3)

式中τ,Lu和Ld分別為大氣透過率、上行輻射(W·m-2·sr-1·μm-1)和下行輻射(W·m-2·sr-1·μm-1)，這些大氣剖面參數(shù)可在網(wǎng)站(http: //atmcorr.gsfc.nasa.gov/)上輸入每幅遙感影像中心經(jīng)緯度與成像時(shí)間信息并由MODTRAN 模型估算獲取[29-30]。雖然地表比輻射率因不同地表結(jié)構(gòu)的材料而不同，但仍可利用多光譜數(shù)據(jù)對不同土地覆蓋條件建立查找表估算像元尺度的地表比輻射率[25,31]，公式為：

(4)

式中: 水體的地表發(fā)射率為0.99;εs為裸土(0.05≤NDVI<0.15)的地表發(fā)射率0.97;εv為均勻植被(NDVI≥0.45)的地表比輻射率0.985；Pv為地表植被覆蓋度。Pv計(jì)算公式為：

(5)

式中：NDVIi為像元的歸一化植被指數(shù)；NDVIs為裸土的NDVI值，取值0.15；NDVIv為植被的NDVI值，取值0.45。NDVI的計(jì)算公式為：

(6)

式中ρNIR和ρR分別為近紅外波段和紅光波段的地表反射率。

為了在統(tǒng)計(jì)分析中減少云對LST的影響，本研究將pixel_qa波段用于過濾云像素，得到逐景影像的無云LST，再計(jì)算2013年3月—2019年2月逐月平均LST用于后續(xù)的統(tǒng)計(jì)分析(由于云量過多，缺少2014年4月、2014年7月和2015年5月的數(shù)據(jù))。另外，本文研究浮島光伏電場不同建設(shè)狀態(tài)(建成后與建設(shè)前)同時(shí)期的無云LST差異，并不嚴(yán)格要求LST計(jì)算的絕對精度，因此該算法的誤差對本研究的影響基本可以忽略。

2.2.2 浮島光伏電場對LSTs-a的影響范圍與季節(jié)差異

為了判斷浮島光伏電場對LST的影響范圍和季節(jié)差異，以浮島光伏建成后與建設(shè)前相同時(shí)期的無云LST差異作為LST的變化幅度是識別該模式和程度的有效手段。然而，建設(shè)區(qū)LST的變化除了受地表覆蓋差異的影響之外，主要受到地表能量平衡(如太陽輻射和云量)的控制和人類活動因素(土地覆蓋變化、地形的人為改變等)的影響[32]，直接將建成后與建設(shè)前的LST相減，其LST差異不能準(zhǔn)確反映浮島光伏電場對LST場的變化情況。所以，在分析的過程中，需要盡量剔除光伏電場之外的其他影響因素。為了盡量剔除人類活動對LST的影響，筆者選擇了與光伏電場建設(shè)前土地覆蓋類型一致且在研究時(shí)段內(nèi)沒有發(fā)生變化的水體作為參照。圖3(a)顯示，伴隨著氣溫波動，建設(shè)前某些時(shí)段的LST普遍高于其他年份相同時(shí)期的LST(如2013年)，浮島光伏電場建成前后的LST變化幅度存在不確定性，不能準(zhǔn)確反映浮島光伏建設(shè)區(qū)LST的變化情況。為了進(jìn)一步剔除自然要素變化對LST的影響，考慮到太陽輻射、地形、地表覆蓋等要素的差異所導(dǎo)致的LST的變化可以在氣溫上體現(xiàn)出來[32-34]，本研究構(gòu)建LSTs-a(逐月LST與對應(yīng)月份的氣溫之差)在一定程度上消除自然要素變化對LST變化幅度的干擾。LSTs-a的計(jì)算公式為：

LSTs-a=LSTi-Ti，

(7)

式中：LSTi為第i月的LST；Ti為第i月的氣溫； 氣溫取以浮島光伏建設(shè)區(qū)為中心60 km緩沖區(qū)范圍內(nèi)2個(gè)氣象臺站的平均值。通過構(gòu)建LSTs-a，一定程度上消除氣溫變化對LST變化幅度的干擾(圖3(b))。

(a) 氣溫與LST隨時(shí)間的變化(b) LSTs-a隨時(shí)間的變化

為了判斷浮島光伏電場對LST的影響模式，選取“夏季”和“冬季”2個(gè)溫度變化最為明顯的季節(jié)進(jìn)行對比分析。按照式(7)提取浮島光伏電場建設(shè)區(qū)及其鄰近區(qū)域夏季(6—8月)、冬季(12月—次年2月)的LSTs-a，計(jì)算浮島光伏電場建成后與建設(shè)前的LSTs-a差異。在此基礎(chǔ)上，利用ArcGIS軟件中的局部Moran’s I來分析浮島光伏電場對鄰域的影響程度，進(jìn)而判斷光伏電場對LST的具體影響距離。

局部Moran’s I是研究空間依賴關(guān)系的重要空間統(tǒng)計(jì)分析方法，它通過特定地點(diǎn)或區(qū)域的空間位置與每個(gè)觀測值的貢獻(xiàn)來識別“集群”位置[35-36]。同時(shí)，該指數(shù)在識別極端(高或低)值的空間格局時(shí)特別有效，從而使其對于評估LST的空間變化格局非常敏感[37-39]。局部Moran’s I取值在-1～1之間，>0表示正自相關(guān)，即表示要素具有包含同樣高或同樣低的屬性值的鄰近要素； <0表示負(fù)自相關(guān)，即要素具有包含不同值的鄰近要素； 等于0則表示屬性值是隨機(jī)分布的。通過統(tǒng)計(jì)學(xué)上的顯著性檢驗(yàn)(P≤0.05)的局部Moran’s I有4種輸出模式： 高值(HH)聚類、低值(LL)聚類、高值主要由低值圍繞的異常值(HL)以及低值主要由高值圍繞的異常值(LH)[37-39]。在本研究中， HH模式代表LSTs-a變化幅度以高值為中心的空間集聚，LL模式代表LSTs-a變化幅度以低值為中心的空間集聚，HL模式代表LSTs-a變化幅度高值主要由低值圍繞的異常值，LH模式代表LSTs-a變化幅度低值主要由高值圍繞的異常值，通過獲取LSTs-a的高值(HH)聚類分布，進(jìn)而得到光伏電場對LST的影響距離。需要說明的是，為了去除單一年份對分析結(jié)果造成的偏差，建成前、后的夏季LST分別取2013—2014年和2017—2018年夏季的平均值，建成前、后階段的冬季LST分別取2013—2014年和2017—2018年冬季的平均值。

2.2.3 浮島光伏電場不同建設(shè)階段對LSTs-a的影響程度

為了明確浮島光伏電場的增/降溫能力，本研究利用ArcGIS軟件的分區(qū)統(tǒng)計(jì)功能提取典型地物類型逐月LST與LSTs-a的平均值和標(biāo)準(zhǔn)差，分析浮島光伏電場不同建設(shè)階段對LST的影響。需要注意的是，典型土地利用類型是結(jié)合Google Earth的歷史影像與無云Landsat8影像無變化區(qū)域進(jìn)行目視解譯得到的。根據(jù)浮島光伏工程的施工時(shí)間段，選取2013年9月和2018年9月無云覆蓋的Landsat8地表反射率影像，利用ENVI 軟件提供的Spectral Angle Difference變化檢測方法得到未變化區(qū)域； 再對照Google Earth的歷史影像，對光伏電場1 km范圍內(nèi)未變化區(qū)域的典型地物進(jìn)行目視解譯。

3 結(jié)果與分析

3.1 浮島光伏電場對LSTs-a影響的范圍與季節(jié)差異

夏季和冬季，浮島光伏電場對LST的影響存在增溫效應(yīng)。如圖4(a)和(b)所示，夏季，原來的水體對照區(qū)和未建浮島光伏電場前的水體區(qū)域形成的“冷島”區(qū)域被切斷，在浮島光伏電場建設(shè)區(qū)形成溫暖區(qū)域； 浮島光伏電場的建設(shè)明顯增加了其占地范圍內(nèi)的LST(圖4(c))，LSTs-a最大增溫幅度超過7 ℃，LSTs-a平均增溫幅度也能達(dá)到4.9 ℃。林地范圍的LST有所降低，LSTs-a平均降幅達(dá)1.1 ℃，形成了一個(gè)局部低溫中心； 耕地靠近浮島光伏電場的邊緣部分的區(qū)域LST有所升高，LSTs-a增溫幅度在3～5 ℃，遠(yuǎn)離浮島光伏電場溫暖區(qū)的耕地LST無明顯變化； 裸地受林地阻擋，其LST無明顯變化。由局部Moran’s I(圖4(d))得知，夏季浮島光伏電場的LST符合HH模式，對鄰近區(qū)域的增溫作用集中在其200 m以內(nèi)區(qū)域，主要為與浮島光伏電場毗鄰的北部的建設(shè)區(qū)水域、西部的耕地以及東南部的水陸邊緣。

(a) 建設(shè)前階段的LSTs-a空間分布(b) 建成階段的LSTs-a空間分布

(c) LSTs-a變化幅度(d) LSTs-a變化幅度的局部Moran’s I模式

冬季，浮島光伏電場的建立增加了其占地范圍內(nèi)的LST(圖5)，LSTs-a平均增溫幅度為3.5 ℃，較夏季增溫幅度??； 原來的水體對照區(qū)和未建浮島光伏電場前的水體區(qū)域形成的“冷島”區(qū)域被切斷，在浮島光伏電場建設(shè)區(qū)形成溫暖區(qū)域； 浮島光伏電場建設(shè)區(qū)與林地、耕地的LST連成一片； 裸地受林地阻擋，其LST無明顯變化。由局部Moran’s I(圖5(d))可知，冬季浮島光伏電場增溫作用較夏季弱，僅停留在浮島光伏電場建設(shè)區(qū)及其200 m范圍內(nèi)的水域。

(a) 建設(shè)前階段的LSTs-a空間分布(b) 建成階段的LSTs-a空間分布

(c) LSTs-a變化幅度(d) LSTs-a變化幅度的局部Moran’s I模式

3.2 浮島光伏電場不同建設(shè)階段對LSTs-a的影響程度

基于單通道算法提取了研究區(qū)典型地物的逐月平均LST，再根據(jù)浮島光伏電場工程的建設(shè)周期，得到建設(shè)前階段(2013年3月—2015年10月)、建設(shè)階段(2015年11月—2016年10月)和建成階段(2016年11月—2019年2月)的LST變化曲線(圖6)。浮島光伏電場建設(shè)前階段，擬建設(shè)區(qū)與水體區(qū)域同為采煤塌陷區(qū)的低洼積水塘，LST的變化曲線一致，且月平均LST為最低(分別為20.2 ℃和19.6 ℃)； 其余由低到高依次為林地、耕地、建設(shè)用地和裸地。浮島光伏電場建設(shè)階段和建成階段，建設(shè)區(qū)的月均LST普遍高于水體，接近于林地的LST，甚至在某些時(shí)段與建設(shè)用地的LST相當(dāng)。根據(jù)局部Moran’s I的結(jié)果，浮島光伏電場對LST的影響范圍主要集中在建設(shè)區(qū)200 m范圍內(nèi)，對其周圍典型地類的增溫效應(yīng)非常微弱。因此，可以采用建設(shè)區(qū)的LSTs-a變化幅度來明確浮島光伏電場不同建設(shè)階段的增溫能力。圖7為浮島光伏電場不同建設(shè)階段的LSTs-a變化，其中建設(shè)前階段的LSTs-a取2013年3月—2015年2月的平均值，建設(shè)階段的LSTs-a取2015年11月—2016年10月的平均值，建成階段的LSTs-a取2017年3月—2019年2月的平均值，陰影部分為LST的標(biāo)準(zhǔn)差。結(jié)果顯示，浮島光伏電場建設(shè)前，擬建設(shè)區(qū)與水體的年均LSTs-a差別很微弱，僅為0.06 ℃； 建設(shè)階段差異逐漸增大(2.77 ℃)； 而浮島光伏電場建成后，二者的年均LSTs-a之差達(dá)到了3.96 ℃。在浮島光伏電場建設(shè)區(qū)，建設(shè)階段的年均LSTs-a之差為3.26 ℃，建成后與建設(shè)前的年均LSTs-a之差為4.50 ℃，而水體區(qū)域僅為0.55 ℃，這說明了浮島光伏電場的建設(shè)使建設(shè)區(qū)低洼積水塘的年均LST升高了4.50 ℃。

(a) 浮島光伏建設(shè)區(qū)LST(b) 水體與浮島光伏建設(shè)區(qū)LST

(c) 林地與浮島光伏建設(shè)區(qū)LST(d) 耕地與浮島光伏建設(shè)區(qū)LST

(e) 裸地與浮島光伏建設(shè)區(qū)LST(f) 建設(shè)用地與浮島光伏建設(shè)區(qū)LST

圖7 浮島光伏電場不同建設(shè)階段的LSTs-a變化

4 結(jié)論與討論

4.1 結(jié)論

本文基于Landsat8時(shí)間序列遙感數(shù)據(jù)，利用單通道算法提取了淮南市典型浮島光伏電場及其鄰近區(qū)域的LST數(shù)據(jù)集，對光伏電場影響LST空間分布的模式和影響的范圍進(jìn)行了分析，并明確了建設(shè)區(qū)不同建設(shè)階段對LST空間分布的影響程度。主要結(jié)論為：

1)浮島光伏電場的建設(shè)明顯改變了建設(shè)區(qū)的熱環(huán)境，在溫度變化最明顯的夏季和冬季都存在增溫效應(yīng)，其中夏季的平均增溫幅度達(dá)到4.9 ℃，冬季平均增溫幅度則為3.5℃。增溫效應(yīng)主要集中在建設(shè)區(qū)200 m范圍內(nèi)，對其周圍典型地類的增溫效應(yīng)非常微弱。

2)浮島光伏電場建設(shè)階段和建成階段，建設(shè)區(qū)的月均LST普遍高于水體，接近于林地的LST； 2個(gè)階段的年均增溫幅度分別為3.26 ℃和4.50 ℃。

4.2 討論

本文對浮島光伏電場不同建設(shè)階段及建成前后不同季節(jié)的熱環(huán)境變化進(jìn)行分析，結(jié)果表明浮島光伏電場對建設(shè)區(qū)及其鄰近200 m以內(nèi)具有增溫效應(yīng)，其增溫幅度夏季大于冬季，同時(shí)為光伏電場對局地環(huán)境影響評價(jià)的相關(guān)研究提供了借鑒。然而，LST受諸多因素影響，其變化是復(fù)雜的，存在眾多不確定性的問題，這對光伏電場的熱環(huán)境效應(yīng)評估提出了新的挑戰(zhàn)。

1)無云時(shí)間序列LST構(gòu)建。大范圍云蓋的存在對Landsat8影像LST反演帶來挑戰(zhàn)，一方面降低了LST的準(zhǔn)確度，另一方面減少了LST時(shí)間序列的連續(xù)性，這些都為后續(xù)光伏電場LST的統(tǒng)計(jì)分析帶來偏差。為了減少這一偏差，以往的研究通常在LST反演前進(jìn)行數(shù)據(jù)篩選，即整幅影像的云量在一定閾值以內(nèi)[18,20]。但是這樣的篩選容易造成漏選研究區(qū)的無云遙感影像，也容易減少LST時(shí)間序列的連續(xù)性。實(shí)際上，筆者也曾考慮過利用時(shí)間分辨率較優(yōu)的MODIS LST 產(chǎn)品用于分析光伏電場的熱環(huán)境變化，但受限于MODIS較低的空間分辨率(1 km)，光伏電場的面積尚不足MODIS產(chǎn)品的一個(gè)像元，不能像Landsat8數(shù)據(jù)一樣展現(xiàn)更多的LST空間細(xì)節(jié)，更不能很好地分析光伏電場對LST的影響范圍。綜合時(shí)間連續(xù)性和空間分辨率，本研究利用GEE在線云平臺獲取研究區(qū)全時(shí)段的Landsat8數(shù)據(jù)集，將pixel_qa波段用于過濾云像素，得到逐景影像的無云LST，再計(jì)算逐月平均LST用于后續(xù)的統(tǒng)計(jì)分析； 在確保LST精度的前提下，可以有效增加LST時(shí)間序列的連續(xù)性。然而，由于云量過多，月尺度下的LST仍然存在缺失。在后續(xù)的研究中，可以考慮引入Savitzky-Golay濾波、HANTS(harmonic analysis of time series)[40-41]等時(shí)間序列重建方法構(gòu)建Landsat無云時(shí)間序列LST。

2)光伏電場增/降溫幅度的分離。由于LST受到地表能量平衡等自然要素的控制和人類活動因素的影響[32-34]，光伏電場建成前后的LST差值不能準(zhǔn)確反映增/降溫幅度?？紤]到太陽輻射、地形、地表覆蓋等要素的差異所導(dǎo)致的LST的變化可以在氣溫上體現(xiàn)出來，本研究計(jì)算逐月LST與對應(yīng)月份的氣溫之差(LSTs-a)以減少光伏電場之外的其他因素對LST的干擾，氣溫取以浮島光伏建設(shè)區(qū)為中心60 km緩沖區(qū)范圍內(nèi)2個(gè)氣象臺站的平均值。當(dāng)然，本文的計(jì)算方法并不能完全消除自然要素和人類活動對LST的影響。LSTs-a中，氣溫取以浮島光伏建設(shè)區(qū)為中心60 km緩沖區(qū)范圍內(nèi)2個(gè)氣象臺站的平均值，與每個(gè)Landsat像元相比較，空間分辨率確實(shí)很低，也會存在誤差。該方法對不同下墊面，特別是異質(zhì)性較強(qiáng)或土地利用類型易于改變的地表，適用性尚不清晰，仍需深入研究。在后續(xù)的研究中，將進(jìn)一步考慮采用氣候?qū)W模型和實(shí)地觀測的方法，使得結(jié)果更加準(zhǔn)確。

3)不同下墊面光伏電場熱環(huán)境效應(yīng)對鄰域的影響范圍、程度與歸因分析。本文的光伏電場位于漂浮水面，它對建設(shè)區(qū)及其鄰近200 m以內(nèi)具有增溫效應(yīng)，但不同下墊面光伏電場影響LST空間分布的模式及影響范圍、程度存在差異，不同地形條件的荒漠區(qū)表現(xiàn)不同的降溫、增溫效應(yīng)[19-20,42-43]。例如，青海共和盆地荒漠區(qū)大型光伏電站的白天LST年均下降幅度為4 ℃[40]； 甘肅省金昌市西坡戈壁光伏電場白天LST年均下降幅度為0.65 ℃[21]。此外，光伏電場影響LST空間分布的驅(qū)動機(jī)制尚不明晰。在后續(xù)研究中，應(yīng)發(fā)展優(yōu)化適用于不同下墊面條件的光伏電場自動化遙感提取方法，并考慮結(jié)合地面觀測、遙感觀測與能量平衡模型對光伏電場影響LST空間分布的模式進(jìn)行綜合分析。