熱環(huán)境視角下的最小生態(tài)安全距離
——以大連北三市為例

楊 俊,白世豪,金 翠,*,喬瑩瑩,3

1 遼寧師范大學(xué) 人居環(huán)境研究中心, 大連 116029 2 東北大學(xué) 江河建筑學(xué)院,沈陽 110169 3 遼寧師范大學(xué) 海洋經(jīng)濟(jì)與可持續(xù)發(fā)展研究中心,大連 116029

2014年3月國家環(huán)保部在全國環(huán)境保護(hù)工作會議中首次提出“最小生態(tài)安全距離(Minimum Distance Forecological Security,MDES)”概念。最小生態(tài)安全距離,主要是為了保障城市內(nèi)部及城市間的生態(tài)安全,并滿足生態(tài)系統(tǒng)正常運(yùn)轉(zhuǎn)以及環(huán)境污染物稀釋擴(kuò)散的需要[1]。這是我國城鎮(zhèn)化進(jìn)程中治理城市間復(fù)合污染問題的全新思路,包含了生態(tài)保護(hù)、安全底線、資源承載力和環(huán)境容量的綜合性概念。城市熱環(huán)境作為城市生態(tài)系統(tǒng)的重要組成部分,隨著城市化進(jìn)程的加快,城市人口的快速增加,城市建筑密度的大幅提升,人們開始逐漸關(guān)注由城市熱島效應(yīng)帶來的一系列負(fù)面影響。目前城市熱環(huán)境不僅影響著城市氣候,還影響著城市化進(jìn)程和生態(tài)環(huán)境的良性發(fā)展。而基于熱環(huán)境的城市群最小生態(tài)安全距離是為防范城市群熱環(huán)境風(fēng)險(xiǎn)和保障人類宜居環(huán)境,要求城市間承載人類經(jīng)濟(jì)社會活動的城市用地間隔一定的距離作為生態(tài)緩沖帶,以最大限度地發(fā)揮生態(tài)用地降溫效益,從而消解城市間熱環(huán)境疊加污染。

國內(nèi)外學(xué)者已對城市規(guī)模和內(nèi)部結(jié)構(gòu)影響城市熱環(huán)境的機(jī)理、城市熱環(huán)境自身空間特征展開了大量理論和實(shí)證研究,對城市熱環(huán)境風(fēng)險(xiǎn)防范也展開了探索性研究。其中,城市規(guī)模顯著影響城市熱島強(qiáng)度已在世界多個(gè)城市被證實(shí)。Imhoff 等[2]分析美國38個(gè)城市熱島強(qiáng)度特征得出,熱島強(qiáng)度與城市面積呈正向線性相關(guān)關(guān)系ΔTUrban-Rural=3.48log(Area)+1.75;Zhang等[3]分析全球3000個(gè)城鎮(zhèn)居民點(diǎn)熱島強(qiáng)度特征指出,當(dāng)居民點(diǎn)面積大于500 km2時(shí),熱島強(qiáng)度達(dá)到4.7℃,居民點(diǎn)面積大于10 km2并小于50 km2時(shí),熱島強(qiáng)度為2.5℃;Peng等[4]分析了全球419個(gè)城市晝夜熱島強(qiáng)度特征,白天城市熱島強(qiáng)度與城市大小呈顯著的正相關(guān)關(guān)系(R2=0.16,P=0.003,56個(gè)歐洲城市),城市形態(tài)也在一定程度上影響城市熱島效應(yīng);Qiao等[5]通過研究1989—2010年北京市建成區(qū)與熱環(huán)境空間關(guān)系指出,城市熱島與城市建設(shè)用地重心轉(zhuǎn)移方向基本一致,當(dāng)城市規(guī)模相近,城市形態(tài)越緊湊,城市熱島效應(yīng)越明顯;Ke等[6]通過WRF(Weather Research and Forecast)模擬2020年武漢都市圈不同城市擴(kuò)展模式下的氣候效應(yīng),增溫效應(yīng)在城市用地集中擴(kuò)展模式下最強(qiáng),城市用地分散擴(kuò)展模式增溫最弱;王曉云[7]通過數(shù)值模式模擬四種500 km2理想條件下城市總體布局方案的氣候特征,組團(tuán)狀城市布局平均溫度最低,其次是十字狀和帶狀布局,點(diǎn)狀集中狀布局平均溫度最高。在城市熱環(huán)境空間結(jié)構(gòu)及形態(tài)研究方面;周紀(jì)等[8]提出并構(gòu)建基于高斯曲面的城市熱島容量模型充分反映了熱島信號的空間變化特征。針對單核心、空間分布基本對稱的城市熱島強(qiáng)度分布,使用二維高斯曲面來反映城市熱環(huán)境的空間形態(tài)具有一定的可操作性。查良松等[9]利用改進(jìn)的半徑法對合肥市的城市熱場狀況進(jìn)行研究,發(fā)現(xiàn)亮溫的水平分布在城市建成區(qū)與非城市建成區(qū)間有突變現(xiàn)象,并以此分析城市熱力場的空間分布,得到熱島強(qiáng)度計(jì)算公式I=TB-T1。此外,學(xué)者已注意到城市熱島足跡和城市熱島容量的晝夜性、季節(jié)性差異[10- 13],以及城市熱環(huán)境的區(qū)域性差異[14- 16],并通過建立熱環(huán)境分區(qū)方案,指導(dǎo)城市各區(qū)域針對引發(fā)城市熱環(huán)境的主因子進(jìn)行重點(diǎn)調(diào)控和防治[17- 19]。通過城市不同區(qū)域之間綠化帶、隔離帶等的設(shè)計(jì)緩解城市熱環(huán)境策略被廣泛提出,城市的建筑空間形態(tài)、通風(fēng)廊道、清潔空氣廊道等案例也被不斷運(yùn)用到緩解城市熱環(huán)境設(shè)計(jì)中[20- 23]。

但至今為止,國內(nèi)外學(xué)術(shù)界對于“最小生態(tài)安全距離”的實(shí)質(zhì)內(nèi)涵還未作出明確的界定,本文基于城市生態(tài)安全格局與緩沖帶等理論基礎(chǔ),定義熱環(huán)境視角下的最小生態(tài)安全距離。量測基于熱環(huán)境視角的城市間最小生態(tài)安全距離,不僅能夠?yàn)樵O(shè)置城市發(fā)展的“剛性界限”和“生態(tài)安全底線”提供理論依據(jù),還能夠?yàn)榻鉀Q城市集中發(fā)展所顯現(xiàn)出來的城市群熱環(huán)境沖突和疊加等新型環(huán)境問題提供科學(xué)參考。鑒此本文通過研究分析大連市普蘭店、瓦房店、莊河三市之間的熱環(huán)境疊加問題,量測城市間基于熱環(huán)境應(yīng)保持的最小生態(tài)安全距離,進(jìn)而判斷城市內(nèi)部用地是否受到了其他城市熱環(huán)境疊加污染。

1 研究區(qū)概況

本文以大連市的北三市,普蘭店市、瓦房店市、莊河市為研究樣區(qū),北三市位于大連市北部,總面積10370.5 km2。隨著大連市經(jīng)濟(jì)的快速發(fā)展,瓦房店與普蘭店兩市之間的建成區(qū)又相距較近,之間較易產(chǎn)生熱環(huán)境疊加問題;而普蘭店與莊河兩市建成區(qū)相距較遠(yuǎn),之間不易出現(xiàn)熱環(huán)境疊加問題,以上兩種情況都對研究“熱環(huán)境視角下最小生態(tài)安全距離”具有代表性意義。此外,瓦房店市西部及莊河市南部有部分島嶼,對城市內(nèi)部的熱島效應(yīng)不成構(gòu)影響,故不予考慮。

圖1 研究區(qū)概況Fig.1 The location of study area

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源

本文研究數(shù)據(jù)采用Landsat8 OLI_TIRS衛(wèi)星數(shù)據(jù)來反演城市的地表溫度,城市建成區(qū)邊界數(shù)據(jù)采用由中科院地理所提供的土地利用/覆被數(shù)據(jù)進(jìn)行提取,城市行政邊界數(shù)據(jù)由大連市規(guī)劃局提供,如表1所示。

表1 數(shù)據(jù)來源及說明

2.2 建成區(qū)邊界提取

建成區(qū)主要由城市不透水面構(gòu)成,是指城市市區(qū)連片集中及分散在近郊與城市有密切聯(lián)系、具有基本完善市政公用設(shè)施的城市建設(shè)用地。在我國土地資源遙感調(diào)查與監(jiān)測技術(shù)規(guī)程中,劉紀(jì)遠(yuǎn)等[24]提出了基于遙感TM影像的土地利用/覆蓋分類系統(tǒng)設(shè)計(jì)原則,并在此基礎(chǔ)上設(shè)計(jì)了完整的土地利用覆蓋分類系統(tǒng),將土地利用/覆蓋分為耕地、林地、草地、水域、建設(shè)用地和未利用土地6個(gè)一級類和有林地、灌木林、疏林地、其他林地以及高、中、低覆蓋度草地等25個(gè)二級類型。本文參照土地利用覆蓋分類系統(tǒng),采用中科院地理所提供的土地利用/覆被數(shù)據(jù),選擇各研究區(qū)的城市建設(shè)用地,并使用ArcGIS提取城市的建成區(qū)邊界。

2.3 地表溫度反演

由于受到大氣、地形、方向等因素的干擾,亮度溫度并不是地表物體輻射的真實(shí)溫度,覃志豪等[25]提出了根據(jù)地表熱輻射傳導(dǎo)方程而不需要大氣校正直接反演地表溫度的單窗算法來反演地表真實(shí)溫度,經(jīng)實(shí)地檢驗(yàn),其演算的精度為0.4—1.1℃,在目前各種反演模型中是一個(gè)比較簡單精度又較高的方法,故本文采用該方法進(jìn)行地表溫度反演(式1),公式如下：

Ts={a10(1-C10-D10)+[b10(1-C10-D10)+C10+D10]T10-D10Ta}/C10

(1)

其中,

C10=τ10ε10,D10=(1-τ10)[1+τ10(1-ε10)]

式中,a10=-67.355351,b10=0.458606,Ta為大氣平均作用溫度,τ10為大氣透射率,ε10為地表輻射率。

2.4 城市群熱環(huán)境空間形態(tài)曲面模擬

究其本質(zhì),影響地表溫度的因素主要有輻射總量、初始地溫和地物的熱慣量等[8]。初始地溫是整個(gè)地表熱力系統(tǒng)的背景溫度,地物熱性質(zhì)的差別是造成地物溫度差異的主要原因。因此,熱力系統(tǒng)的地表溫度可以表示為(式2)：

Ts(x,y)=Ts1(x,y)+Ts2(x,y)

(2)

式中,(x,y)分別為像元坐標(biāo);Ts(x,y)為地表溫度;Ts1(x,y)為背景溫度;Ts2(x,y)為熱力系統(tǒng)局部熱異常導(dǎo)致的溫度增量。

當(dāng)存在顯著的城市熱異常疊加在鄉(xiāng)村地表溫度場之上時(shí),即出現(xiàn)城市熱島效應(yīng)。由于鄉(xiāng)村地表的熱性質(zhì)相對穩(wěn)定,則其溫度R(x,y)可認(rèn)為熱力系統(tǒng)的背景溫度,將城市熱島導(dǎo)致的溫度增量Ts2(x,y)定義為熱島信號函數(shù)UHI(x,y)。因此,熱力系統(tǒng)的地表溫度可轉(zhuǎn)換為(式3)：

Ts(x,y)=R(x,y)+UHI(x,y)

(3)

其中,

UHI(x,y)=a0×e-U∕2

式中,a0為高斯曲面高度,代表熱島強(qiáng)度;ax,ay分別為高斯曲面底部橢圓的長半軸和短半軸,可以用UHI(x,y)測量內(nèi)的經(jīng)度范圍值和緯度范圍值代替;x0,y0分別為橢圓圓心,即熱島中心;φ為熱島方向角度[8]。

已有的模擬城市熱環(huán)境曲面的方法包括快速傅里葉變換的非參數(shù)化模型、高斯模型和支持向量機(jī)算法。高斯混合模型利用高斯概率密度函數(shù)(正態(tài)分布曲線)精確地量化事物,將一個(gè)事物分解為若干的基于高斯概率密度函數(shù)形成的模型。城市群空間擴(kuò)展及其對熱環(huán)境影響的動態(tài)性和背景場的不確定性會制約曲面模型精度,而高斯混合模型對于運(yùn)動目標(biāo)檢測提取是最為成功的方法之一。因此本研究采用高斯混合模型擬合城市群熱環(huán)境曲面[26],即通過高斯混合模型(GMM),分別抽象出鄉(xiāng)村地表溫度場R(x,y)和熱島信號函數(shù)UHI(x,y),即可擬合城市熱環(huán)境曲面。其數(shù)學(xué)定義如下(式4)：

(4)

(5)

EM(最大期望算法)算法主要用來計(jì)算基于不完全數(shù)據(jù)的極大似然估計(jì),它是一種迭代算法,每一次迭代都能保證似然函數(shù)值增加,并且收斂到一個(gè)局部極大值。故本文使用EM算法求解GMM參數(shù)(式5)[27]。

2.5 城市最小生態(tài)安全距離測算

根據(jù)城市生態(tài)安全格局與緩沖帶等理論基礎(chǔ)[28- 30],假定A和B分別為已建成和規(guī)劃在建的兩個(gè)相鄰大型城市,且A、B兩城市間相隔距離固定,城市功能定位、城市熱環(huán)境特征明確,熱島中心均處于城市中心位置且均勻向外輻射。城市A的熱輻射能力固定,城市B不定。在此條件下,根據(jù)A、B城市熱島信號函數(shù)UHIA(x,y)和UHIB(x,y)衰減距離判斷現(xiàn)有建設(shè)或規(guī)劃條件下城市群熱環(huán)境是否存在風(fēng)險(xiǎn)。構(gòu)建數(shù)學(xué)模型如下：

圖2 最小生態(tài)安全距離示意圖Fig.2 The sketch map of minimum distance forecological security

設(shè)城市A的中心坐標(biāo)(XA-YA),半徑為RA。圖2中黑色圓環(huán)為城市A的建成區(qū),綠色圓環(huán)表示城市A的熱島信號函數(shù)UHIA(x,y)衰減為0,即圓環(huán)外區(qū)域?qū)⒉辉偈艹鞘蠥的熱輻射影響。此時(shí),城市A的建成區(qū)邊緣到熱島衰減緩沖帶邊緣的距離為d1。設(shè)城市B的中心坐標(biāo)為(XB-YB),半徑為RB。假設(shè)城市B的熱環(huán)境特征存在三種情況,即城市B的熱島信號函數(shù)UHIA(x,y)衰減為0時(shí)距離B城市建成區(qū)邊緣的距離分別為d2、d3和d4。其中,當(dāng)城市A的熱島信號函數(shù)UHIA(x,y)與城市B的熱島信號函數(shù)UHIB(x,y)同時(shí)衰減到0時(shí),即UHIA(x,y)+UHIB(x,y)=0時(shí),此時(shí)兩城市熱島衰減緩沖帶圓環(huán)相切(綠色圓環(huán)),d1+d2即為城市A和B的最小生態(tài)安全距離,其中約束條件為,城市A、B之間的熱島信號函數(shù)疊加后的最小熱島強(qiáng)度不超過某一安全閾值I(式6)。

min{UHIA(x,y)+UHIB(x,y)}

(6)

式中,當(dāng)城市熱輻射等于背景場溫度時(shí),熱島信號函數(shù)為0,此時(shí)城市A、B的熱環(huán)境相互不干擾,即I= 0。

當(dāng)d3d2時(shí)(紅色圓環(huán)),即熱島信號函數(shù)UHIA(x,y)和UHIB(x,y)重疊,兩城市熱島衰減緩沖帶圓環(huán)相交,說明A、B兩城市熱環(huán)境污染疊加,已發(fā)生熱環(huán)境沖突,必須通過優(yōu)化生態(tài)用地空間布局改善城市熱環(huán)境問題。

3 結(jié)果與分析

3.1 地表溫度特征

本文采用單窗算法對Landsat8 OLI_TIRS衛(wèi)星數(shù)據(jù)進(jìn)行地表溫度反演,反演之后的結(jié)果數(shù)據(jù)通過ArcGIS建立格網(wǎng)來提取地表溫度值,格網(wǎng)精度為30m×30m。結(jié)果如圖3所示,城市建成區(qū)的地表溫度高于城市郊區(qū)地表溫度,呈現(xiàn)出明顯的城市熱島效應(yīng)現(xiàn)象,整體呈南部溫度高于城市北部的趨勢。2013年至2016年,瓦房店市與普蘭店市在2015年地表溫度值最高,莊河市則在2016年時(shí)地表溫度達(dá)到最高值。地表溫度的增長趨勢,逐漸由瓦房店市與普蘭店市建成區(qū)位置,逐漸向莊河市方向蔓延。

圖3 2013—2016年北三市地表溫度Fig.3 Land surface temperature of the northern three districts during 2013—2016

3.2 城市群熱環(huán)境空間形態(tài)曲面模擬

利用普蘭店市、瓦房店市、莊河市的地表溫度反演結(jié)果、高斯混合模型的聚類分析結(jié)果,同時(shí)結(jié)合熱島信號函數(shù),分別模擬了各研究區(qū)2013年至2016年間的背景場曲面以及基于熱島信號函數(shù)建立的熱環(huán)境曲面。由于比例縮放問題,圖上變化不大,故本文分別給出背景場曲面圖3幅(圖4),UHI曲面圖3幅(圖5)。從左至右,分別是普蘭店市、瓦房店市、莊河市。其中,通過擬合的背景場曲面,普蘭店市的結(jié)果精度參數(shù)如表2所示,當(dāng)誤差平方和(SSE)、均方根誤差(RMSE)越小,負(fù)相關(guān)系數(shù)(R-Square)越接近1時(shí)表明擬合的程度越好。

表2 背景場擬合精度參數(shù)

由圖4可明顯看出普蘭店市與莊河市的背景溫度場都在自北向南的方向上呈現(xiàn)溫度遞增的趨勢,而瓦房店市的背景溫度場在自東向西的方向上呈現(xiàn)溫度遞增的趨勢。同時(shí)結(jié)合圖5,可以明顯看出普蘭店市與莊河市所產(chǎn)生的熱島效應(yīng)都在城市的南部出現(xiàn)峰值,并向城市的四周進(jìn)行輻射的特征,而瓦房店市所產(chǎn)生的熱島效應(yīng)在城市的西部出現(xiàn)峰值,并向城市四周進(jìn)行輻射的特征,且擬合所得到模型與各年份熱島方向及分布趨勢相同。

圖4 北三市背景場曲面Fig.4 The Background surface of northern three districts

圖5 北三市UHI曲面Fig.5 The UHI surface of northern three districts

3.3 城市群間最小生態(tài)安全距離測算

城市建成區(qū)與UHI范圍線之間的部分即該城市基于熱環(huán)境條件下應(yīng)保持的最小生態(tài)安全距離,即通過對圖上量測相鄰城市間建成區(qū)邊界到歷年UHI范圍線間的距離之和。本文擬合得到的熱環(huán)境曲面與背景場溫度的差值,即為擬合后的熱島特征值,選取熱島特征值大于0時(shí)的平面位置坐標(biāo),并利用ArcGIS對其進(jìn)行矢量化,可判斷出各個(gè)研究區(qū)的UHI輻射范圍,即可求出相鄰城市間的最小生態(tài)安全距離。由圖6可知,在2013年至2014,普蘭店市與瓦房店市還未出現(xiàn)熱環(huán)境疊加的情況,直至2015年與2016年時(shí),普蘭店市與瓦房店市熱環(huán)境出現(xiàn)疊加的情況,則說明2015年之后普蘭店市與瓦房店市之間已經(jīng)超過了兩座城市基于熱環(huán)境應(yīng)該保持最小生態(tài)安全距離,兩城市之間生態(tài)環(huán)境受到了彼此城市的熱環(huán)境污染。而普蘭店市、莊河市在2013至2016年間,因建成區(qū)相隔較遠(yuǎn),所以沒有出現(xiàn)熱島效應(yīng)疊加的情況,所以兩座城市間的基于熱環(huán)境的最小生態(tài)安全距離也相隔較遠(yuǎn),兩城市間的生態(tài)環(huán)境相對安全。

圖6 2013—2016年最小生態(tài)安全距離空間分布Fig.6 The MDES spatial distribution from 2013 to 2016

圖7 2013—2016年最小生態(tài)安全距離隨時(shí)間變化Fig.7 The MDES varied with time between 2013 and 2016

在2013年至2016年,普蘭店市、瓦房店市、莊河市的城市熱島面積都有逐漸擴(kuò)大的趨勢(圖6)。這4年間,普蘭店市的熱島面積從264.20 km2擴(kuò)大到726.60 km2,其中,2013年至2014年擴(kuò)大最為明顯,擴(kuò)大了248.60 km2,平均每年擴(kuò)張115.50 km2;瓦房店市的熱島面積從217.40 km2擴(kuò)大到630.20 km2,平均每年擴(kuò)張103.40 km2;莊河市的熱島面積從212.30 km2擴(kuò)大到506.50 km2,平均每年擴(kuò)張73.50 km2。同時(shí),城市間的熱環(huán)境視角下的最小生態(tài)安全距離也出現(xiàn)逐漸衰減的趨勢。其中,普蘭店市與瓦房店市之間的最小生態(tài)安全距離由8.70 km衰減至-3.10 km,平均每年衰減2.93 km;而莊河市與普蘭店市之間的最小生態(tài)安全距離由73.47 km衰減至58.24 km,平均每年衰減3.80 km。

4 結(jié)論與討論

城市熱環(huán)境對局地氣候、人居環(huán)境質(zhì)量的影響日益加劇。定量計(jì)算熱環(huán)境視角下的最小生態(tài)安全距離能夠反映城市熱環(huán)境疊加的污染程度。通過多年城市建成區(qū)與UHI范圍線的比較,可以對城市間最小生態(tài)安全距離進(jìn)行歷史和區(qū)域比對,從而為城市管理、規(guī)劃決策部門提供客觀、連續(xù)的衡量指標(biāo),具有較高的科學(xué)意義和實(shí)用價(jià)值。本文在熱環(huán)境視角下的最小生態(tài)安全距離的實(shí)質(zhì)內(nèi)涵的基礎(chǔ)上,通過建立高斯混合模型和建立熱島信號函數(shù),對2013年至2016年間大連市北三市的城市熱島面積和最小生態(tài)安全距離的測算。通過分析得到以下結(jié)論：

(1)2013年至2016年,大連市北三市的城市熱島面積逐漸擴(kuò)張,總面積由694.10 km2擴(kuò)張到1864.30 km2,且普蘭店市與瓦房店市在2015年之后熱環(huán)境出現(xiàn)疊加的情況,表明兩城市間的生態(tài)環(huán)境受到了熱環(huán)境的污染。

(2)這4年間,大連市北三市熱環(huán)境視角下的最小生態(tài)安全距離逐漸衰減。其中,普蘭店市與瓦房店市之間的最小生態(tài)安全距離呈現(xiàn)逐漸惡化的現(xiàn)象,由8.70 km衰減至-3.10 km;而普蘭店市與莊河市之間的最小生態(tài)安全距離,由73.47 km衰減至58.24 km,仍很安全。

綜上所述,從2013年至2016年大連市北三市城市內(nèi)部熱島效應(yīng)逐漸明顯,熱島面積不斷擴(kuò)大,熱環(huán)境視角下的最小生態(tài)安全距離逐年縮減,甚至出現(xiàn)熱環(huán)境疊加部分,導(dǎo)致最小生態(tài)安全距離呈現(xiàn)出負(fù)值。為解決熱環(huán)境視角下得城市間最小生態(tài)安全距離,建議實(shí)施一些緩解城市內(nèi)部熱島效應(yīng)的措施,以確保城市間的最小生態(tài)安全距離。

本文選取的研究區(qū)較為分散,其熱環(huán)境疊加和未疊加現(xiàn)象較為明顯,通過分析即可判斷城市間的最小生態(tài)安全距離,還有很多大型城市存在多核現(xiàn)象以及城市相鄰式發(fā)展現(xiàn)象,所引發(fā)城市熱環(huán)境問題的機(jī)理更為復(fù)雜和繁瑣,這將在今后的工作中展開研究。此外,本研究存在很多因素限制實(shí)驗(yàn)精度,比如遙感影像的精度、通過ArcGIS對研究區(qū)重采樣格網(wǎng)的大小、以及通過UHI擬合熱環(huán)境曲面的精度等,這些問題會在今后的研究中改進(jìn)。