吳星宇,李 雪,吳桂桔
(中國地震局地震研究所 地震大地測量重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430071)
地震是人類遭遇的最嚴(yán)重的自然災(zāi)害之一,其突發(fā)性和危害性位于其他自然災(zāi)害之首。地震災(zāi)害可導(dǎo)致人員傷亡,建筑物損毀,耕地、林地破壞等,對社會經(jīng)濟(jì)造成極其嚴(yán)重的影響。由于城市規(guī)模的擴(kuò)大和現(xiàn)代化水平的提高,人類社會在地震災(zāi)害下的暴露程度正在急劇增加,人類社會面臨的地震風(fēng)險(xiǎn)正逐步增大。
隨著我國城鎮(zhèn)化的不斷發(fā)展,人口高度集中且發(fā)展快速的城市正在逐漸取代農(nóng)村,成為更高層次的人居環(huán)境形式[1],以城市為承災(zāi)體的地震災(zāi)害所占比例也在不斷升高。破壞性地震的發(fā)震間隔往往長達(dá)幾十年甚至數(shù)百年,而城市化建設(shè)的發(fā)展速率遠(yuǎn)大于地震活動的孕育周期。若不能及時(shí)、有效地評估城市變化導(dǎo)致的地震災(zāi)害風(fēng)險(xiǎn)變化,將會對人民生命財(cái)產(chǎn)安全和社會可持續(xù)發(fā)展構(gòu)成重大威脅[2-4]。
當(dāng)前災(zāi)害風(fēng)險(xiǎn)研究領(lǐng)域?qū)τ陲L(fēng)險(xiǎn)的定義略有不同。Maskrey[5]將自然災(zāi)害風(fēng)險(xiǎn)定義為自然災(zāi)害的總損失;Morgan等[6]則認(rèn)為風(fēng)險(xiǎn)是災(zāi)害可能影響和損失的暴露性;Smith[7]則直接給出風(fēng)險(xiǎn)等于致災(zāi)因子發(fā)生概率和損失的乘積。Dela Cruz Reyan[8]定義風(fēng)險(xiǎn)=(致災(zāi)因子×暴露性×脆弱性)/災(zāi)害準(zhǔn)備。綜合各方觀點(diǎn),地震災(zāi)害風(fēng)險(xiǎn)可由地震災(zāi)害危險(xiǎn)性(hazard)、暴露性(exposure)和承災(zāi)體的脆弱性或易損性(vulnerability)三個(gè)因素相互綜合作用決定。
當(dāng)前地震災(zāi)害風(fēng)險(xiǎn)研究主要關(guān)注致災(zāi)因子和承災(zāi)體脆弱性等方面。李陶[9]采用地震烈度、地震動峰值加速度等地震參數(shù)模擬致災(zāi)因子,建立了地震風(fēng)險(xiǎn)模糊綜合評價(jià)體系,并用該方法評價(jià)了畢節(jié)市的地震風(fēng)險(xiǎn)。曾今艷等[10]采用層次分析法利用地震參數(shù)構(gòu)建地震危險(xiǎn)性指數(shù),研究了運(yùn)城市區(qū)域內(nèi)的地震風(fēng)險(xiǎn)。張蓓蕾等[11]通過結(jié)合德爾菲法與層次分析法利用地震動參數(shù)構(gòu)建致災(zāi)因子危險(xiǎn)性指標(biāo),并評估了寧波市的地震風(fēng)險(xiǎn)。Kijko等[12]利用對建筑物的易損性分析的結(jié)果以及建筑物的基本場地條件數(shù)據(jù)分析了Tulbagh地區(qū)的潛在地震風(fēng)險(xiǎn)。Preciado等[13]從建筑學(xué)角度對墨西哥城歷史城區(qū)(La Merced街區(qū))的166座歷史建筑進(jìn)行建筑特性分析,采用定性和定量相結(jié)合的簡化地震易損性評價(jià)方法,對所選區(qū)域進(jìn)行了地震易損性評價(jià)。Rajarathnam等[14]采用基于GIS平臺的航空影像輔助快速視覺篩選技術(shù)提取了Chennai地區(qū)的建筑物特征進(jìn)行城市地震風(fēng)險(xiǎn)評估。Pavié等[15]對Osijek市的建筑物創(chuàng)建了暴露模型,并分析了其潛在的地震風(fēng)險(xiǎn)。
傳統(tǒng)地震災(zāi)害風(fēng)險(xiǎn)研究對于暴露性的研究并不充分。然而,城市地區(qū)人口和社會經(jīng)濟(jì)活動的急劇增加導(dǎo)致了世界各地城市的快速擴(kuò)張[16,17],而城市的大規(guī)模擴(kuò)張同時(shí)也導(dǎo)致了承災(zāi)體暴露性的快速增長。
本文旨在利用遙感數(shù)據(jù)獲取中國西南川滇地區(qū)10年尺度的城市變化情況,結(jié)合地震危險(xiǎn)性數(shù)據(jù)分析川滇地區(qū)的地震風(fēng)險(xiǎn)變化,為當(dāng)?shù)胤罏?zāi)減災(zāi)工作提供決策支持。
川滇地區(qū)位于中國西南部地區(qū)地處青藏高原東南緣,受到印度板塊與亞歐板塊的碰撞及揚(yáng)子盆地的阻擋作用,是我國地震高發(fā)區(qū)域。川滇地區(qū)地貌豐富,山地、丘陵、平原等均有分布。主要包括青藏高原、云貴高原、四川盆地以及橫斷山區(qū)[18]。川滇地區(qū)常住人口13 233.3萬人,城鎮(zhèn)人口6 881.1萬人,城鎮(zhèn)化率為52.0 %。
該地區(qū)主要活動斷裂包括板塊邊界斷裂龍門山斷裂帶,構(gòu)成川滇菱形塊體東部邊界與北部邊界的鮮水河、則木河、安寧河、小江斷裂帶及南邊界的紅河斷裂帶等。
龍門山斷裂自東北向西南沿著四川盆地西北邊緣分布,是青藏高原和四川盆地的分界線。該斷裂帶整體為東北-西南走向,全長約 500 km,斷層傾向北西向,傾角 50~70°。
鮮水河—小江斷裂位于青藏高原東側(cè),是我國西南地區(qū)的一條極為特殊的地震構(gòu)造帶,斷裂帶全長 400 km,總體呈北東傾向,局部地區(qū)呈現(xiàn)南西向,傾角為 55~80°。安寧河斷裂帶全長約360 km,總體走向近南北向,斷層傾向?yàn)闁|,傾角約為 50~80°。則木河斷裂帶長約110 km,總體走向?yàn)楸北蔽?,斷層傾向?yàn)楸睎|向。小江斷裂帶長約 400 km,斷層傾向東方向,傾角約為 70°。
紅河斷裂位于華南地塊與印支地塊之間,是川滇菱形塊體的西南邊界,該斷裂帶整體呈北西走向,在中國境內(nèi)長約 600 km,斷層傾角為 60~80°。
1970~2020年川滇地區(qū)共發(fā)生6級以上地震66次,其中7級以上大地震共13次(圖1)。
圖1 研究區(qū)斷層與歷史地震分布Fig.1 Distribution map of faults and historical earthquakes in Sichuan-Yunnan region注:圖中統(tǒng)計(jì)了1950~2010年震級Ms5.0以上、震源深度小于30 km的地震。符號大小表示震級大小,符號越大,震級越大;顏色深淺表示時(shí)間變化,顏色越深,地震發(fā)生時(shí)間越近。
本文以川滇地區(qū)為例,通過土地利用數(shù)據(jù)獲取研究區(qū)內(nèi)城鎮(zhèn)建設(shè)用地,采用景觀指數(shù)法研究城市分布的時(shí)空變化,結(jié)合研究區(qū)地震危險(xiǎn)性數(shù)據(jù),構(gòu)建了地震風(fēng)險(xiǎn)變化指數(shù)來描述時(shí)間尺度內(nèi)研究區(qū)的地震風(fēng)險(xiǎn)變化。
本研究所采用的城市空間分布數(shù)據(jù)來源于2000版和2010版全球土地覆蓋數(shù)據(jù)庫(GlobeLand30),它是世界上第一個(gè)30 m分辨率的全球土地利用數(shù)據(jù)集[19,20]。該數(shù)據(jù)集由2000多幅Landsat衛(wèi)星和中國HJ-1衛(wèi)星影像提取,覆蓋南北緯80°的陸地范圍內(nèi)10種地表覆蓋類型?;贏rcGIS重分類、重采樣工具獲取人造地表覆蓋類型作為城市空間分布數(shù)據(jù)。為了避免研究區(qū)范圍過大造成數(shù)據(jù)量過大,重采樣的空間分辨率設(shè)為200 m。
地震對城市的威脅不僅來自于斷層附近的地表破裂,更多來自于地震動導(dǎo)致的地表位移和形變[21]。因此,本文使用的地震動峰值加速度 (Peak Ground Acceleration,PGA)數(shù)據(jù)來模擬地震致災(zāi)因子,數(shù)據(jù)源自《中國地震動參數(shù)區(qū)劃圖GB18306-2015》。通過ArcGIS將矢量PGA數(shù)據(jù)轉(zhuǎn)換為柵格數(shù)據(jù),并重采樣至城市空間分布數(shù)據(jù)相同的地面分辨率(圖2)。
圖2 川滇地區(qū)PGA數(shù)據(jù)Fig.2 PGA data of Sichuan-Yunnan region
城市發(fā)展的集聚效應(yīng)會影響城市空間分布模式的變化[22]。因此,相比于基于像素的分析方法,基于斑塊對象的分析方法更適合于開展城市空間分布變化分析。本文引入景觀生態(tài)學(xué)中的景觀指數(shù)斑塊總面積 (Class Area,CA)與斑塊數(shù)量 (Number of Patches,NP)開展城市空間分布變化分析。
CA代表城市建筑的面積:
CA=∑Pi×A
(1)
其中,Pi表示第i斑塊所包含的像元數(shù)量;A為像元大小。
NP代表城市圖斑數(shù)量。NP與CA可以反映城市的空間分布格局。根據(jù)兩個(gè)時(shí)期城市斑塊CA、NP的數(shù)值變化,可以確定城市空間分布的變化模式。
本文以區(qū)縣為統(tǒng)計(jì)單位,利用直方圖百分比方法[23]自動確定CA、NP指數(shù)的變化率區(qū)間。區(qū)間組距計(jì)算公式如下:
(2)
其中,H表示組距,IQR是數(shù)據(jù)的四分位距,N是樣本容量。
每個(gè)指標(biāo)定義了三個(gè)變化區(qū)間。CA變化介于[-H,H]之間的區(qū)縣視為城市面積穩(wěn)定;CA變化大于H的區(qū)縣視為城市面積增加;CA變化小于-H的區(qū)縣視為城市面積減少。同理對NP可進(jìn)行相同的分類。
根據(jù)同一區(qū)縣的CA和NP變化情況,可以確定城市空間分布的變化情況(圖3)。當(dāng)CA增加且NP增加時(shí),表示城市空間分布呈離散型增長模式。當(dāng)CA減少且NP減少時(shí),表示城市空間分布呈聚集型減少模式。
圖3 城市變化特征Fig.3 Characteristics of urban change
地震災(zāi)害風(fēng)險(xiǎn)是由地震災(zāi)害危險(xiǎn)性(Hazard)、暴露性(Exposure)和承災(zāi)體的脆弱性或易損性(Vulnerability)三個(gè)因素相互綜合作用而成,假設(shè)地震風(fēng)險(xiǎn)指數(shù)可以表示為:
EqR=∑HiViEi
(3)
式中,EqR為地震風(fēng)險(xiǎn)指數(shù)。Hi為像素i的地震危險(xiǎn)性,即地震動峰值加速度,Hi值越大表示該地區(qū)的地震危險(xiǎn)性越大。Vi為建筑的脆弱性指數(shù),代表建筑結(jié)構(gòu)在地震作用下發(fā)生某種程度破壞的可能性。Ei為每個(gè)像素的暴露性指數(shù)。
對于T1、T2時(shí)刻的地震風(fēng)險(xiǎn)變化指數(shù)可以表示為:
ΔEqR=EqRT2-EqRT1
(4)
其中,EqRT1和EqRT2分別表示T1,T2時(shí)刻的地震風(fēng)險(xiǎn)指數(shù)。
(5)
由式(5)可見,城市地震風(fēng)險(xiǎn)變化與城市建筑的空間分布變化及其所在位置的地震危險(xiǎn)性相關(guān)。
設(shè)地震風(fēng)險(xiǎn)變化率為η,則有:
(6)
由式(6)可見,城市地震風(fēng)險(xiǎn)變化率與城市的暴露性相關(guān)。對于地震風(fēng)險(xiǎn)變化指數(shù)和地震風(fēng)險(xiǎn)變化率,同樣以區(qū)縣為統(tǒng)計(jì)單位,利用直方圖百分比方法自動確定變化區(qū)間。
城市地震風(fēng)險(xiǎn)變化分析的技術(shù)路線圖見圖4。
圖4 技術(shù)路線Fig.4 Technical roadmap
圖5 研究區(qū)城市空間分布變化情況Fig.5 Changes of urban spatial distribution in Sichuan-Yunnan region
川滇地區(qū)2000~2010年間城市空間分布變化情況如圖6所示??梢?0.9 %的區(qū)縣城市面積有所增加,37.8 %保持穩(wěn)定,11.3 %CA減少,說明城市面積增加的區(qū)縣占大多數(shù)。而城市圖斑數(shù)量變化相對平均,其增加、穩(wěn)定、減少的區(qū)縣比例分別為30 %、37.5 %和32.5 %。城市空間分布未見明顯變化的區(qū)縣數(shù)量最多,其次為離散型增長的區(qū)縣和聚集型增長的區(qū)縣。其中離散型增長的區(qū)縣主要以黑水縣、美姑縣、大姚縣為中心集中分布,而聚集型增長的區(qū)縣分布較為均勻,而城市面積減少的區(qū)縣較少,主要有汶川縣、名山縣、富寧縣等。
川滇地區(qū)2000~2010年地震風(fēng)險(xiǎn)變化率以及地震風(fēng)險(xiǎn)變化指數(shù)分布情況及其統(tǒng)計(jì)規(guī)律如圖6所示。根據(jù)直方圖百分比方法確定的地震風(fēng)險(xiǎn)變化指數(shù)的統(tǒng)計(jì)結(jié)果為:變化指數(shù)小于-11.45視為減少,有25個(gè)區(qū)縣;介于[-11.45,11.45]之間的視為穩(wěn)定,有125個(gè)區(qū)縣;變化指數(shù)大于11.45視為增加,有133個(gè)區(qū)縣。根據(jù)圖6(a)可以看出,地震風(fēng)險(xiǎn)變化指數(shù)增長的區(qū)縣主要集中在川滇地區(qū)南部,如保山市與鶴慶縣。減少的區(qū)縣在名山縣周邊存在集中分布情況,而穩(wěn)定的區(qū)縣分布相對均勻。
地震風(fēng)險(xiǎn)變化率的統(tǒng)計(jì)結(jié)果分別為:變化率小于-10.80 %視為減少,有29個(gè)區(qū)縣;介于[-10.80 %, 10.80 %]視為穩(wěn)定,有115個(gè)區(qū)縣;變化率大于10.80 %視為增加,有139個(gè)區(qū)縣。其分布情況如圖6(b)所示,可見變化率增長明顯的區(qū)縣主要呈現(xiàn)出以黑水縣、美姑縣和大姚縣為中心的集中分布情況,而變化率穩(wěn)定與降低的區(qū)縣分布較為均勻。
圖6 2000~2010年研究區(qū)地震風(fēng)險(xiǎn)變化指數(shù)及變化率Fig.6 Change index and change rate of earthquake risk in Sichuan Yunnan region from 2000 to 2010
根據(jù)地震風(fēng)險(xiǎn)變化指數(shù)的分布情況可知,川滇地區(qū)2000~2010年地震風(fēng)險(xiǎn)變化指數(shù)增加的區(qū)縣共計(jì)133個(gè),約占總數(shù)的一半;穩(wěn)定的區(qū)縣有125個(gè)、減少的區(qū)縣有25個(gè)。在地震風(fēng)險(xiǎn)變化指數(shù)增加的區(qū)縣中,城市空間分布呈離散型增加的區(qū)縣有52個(gè),主要分布在川滇西南部地區(qū)和北部地區(qū);城市空間分布呈聚集型增加的區(qū)縣有33個(gè),主要分布在成都市、昆明市的周邊地區(qū)。在地震風(fēng)險(xiǎn)變化指數(shù)保持穩(wěn)定的區(qū)縣中,城市空間分布呈穩(wěn)定型的區(qū)縣有51個(gè),主要分布在川滇東南地區(qū)。城市空間分布呈現(xiàn)聚集型穩(wěn)定的區(qū)縣有20個(gè),在研究區(qū)內(nèi)分布較為均勻。而呈穩(wěn)定型減少的區(qū)縣有1個(gè),即元陽縣(元陽縣城市面積減少15.63 %,圖斑數(shù)量變化為0,屬穩(wěn)定型減少變化模式。而元陽縣的地震風(fēng)險(xiǎn)變化指數(shù)為-6.45,視為穩(wěn)定,地震風(fēng)險(xiǎn)變化率則為-17.09 %,視為減少)。在地震風(fēng)險(xiǎn)變化指數(shù)減少的區(qū)縣中,城市空間分布呈現(xiàn)聚集型減少的區(qū)縣有16個(gè),主要分布在汶川縣、名山縣周邊地區(qū)。
以上數(shù)據(jù)顯示,地震風(fēng)險(xiǎn)變化指數(shù)增長的地區(qū),其城市空間分布變化大多呈離散型增長模式。地震風(fēng)險(xiǎn)變化指數(shù)保持穩(wěn)定的地區(qū)多為穩(wěn)定型變化的地區(qū)。地震風(fēng)險(xiǎn)變化指數(shù)減少的地區(qū),多表現(xiàn)為聚集型減少的變化模式。
在地震風(fēng)險(xiǎn)變化指數(shù)視為增加的133個(gè)區(qū)縣中,有22個(gè)區(qū)縣的地震風(fēng)險(xiǎn)變化率屬于穩(wěn)定范疇。在地震風(fēng)險(xiǎn)變化指數(shù)視為穩(wěn)定的125個(gè)區(qū)縣中,有12個(gè)區(qū)縣的地震風(fēng)險(xiǎn)變化率視為減少,28個(gè)區(qū)縣的地震風(fēng)險(xiǎn)變化率視為增加。在地震風(fēng)險(xiǎn)變化指數(shù)視為減少的區(qū)縣中,有8個(gè)區(qū)縣的地震風(fēng)險(xiǎn)變化率視為穩(wěn)定。上述區(qū)縣中大多數(shù)都表現(xiàn)出地震風(fēng)險(xiǎn)變化率與城市聚集性變化的正相關(guān)性。
結(jié)合川滇地區(qū)風(fēng)險(xiǎn)變化指數(shù)散點(diǎn)圖(圖7),美姑縣、河口縣與色達(dá)縣的地震風(fēng)險(xiǎn)變化率均大于75.60 %,地震風(fēng)險(xiǎn)變化指數(shù)增加分別為39.35、4.85和32.85。且上述三者均屬于離散型增加城市。相反地,昆明市與成都市則是變化指數(shù)增加較多均大于57.25,地震風(fēng)險(xiǎn)變化率較低分別為38.45 %和26.41 %,且兩者均屬于聚集型增長城市。永勝縣和溫江縣的地震風(fēng)險(xiǎn)變化指數(shù)均大于57.25,地震風(fēng)險(xiǎn)變化率兩者也均大于75.6 %,且溫江縣的地震風(fēng)險(xiǎn)變化率要高于永勝縣,而永勝縣的城市空間分布變化屬于離散型增長,溫江縣則屬于聚集型增長。以上說明了城市城市的聚集性發(fā)展對于地震風(fēng)險(xiǎn)變化率的影響較大。
圖7 2000~2010年研究區(qū)地震風(fēng)險(xiǎn)變化指數(shù)散點(diǎn)圖Fig.7 Scatter diagram of seismic risk change index in Sichuan-Yunnan region from 2000 to 2010
結(jié)合研究區(qū)地貌分布圖(圖8),川滇地區(qū)地震高風(fēng)險(xiǎn)指數(shù)顯著升高地區(qū)的地形主要為中高海拔山地地區(qū),如美姑縣、大姚縣; 其次是高海拔山地地區(qū),如黑水縣;以及少數(shù)低海拔山地地區(qū),如蒼溪縣。四川盆地對應(yīng)的四川中東部的丘陵地區(qū)地震風(fēng)險(xiǎn)變化指數(shù)總體也呈現(xiàn)出增長趨勢,但漲幅相對中高海拔地區(qū)較低。川滇地區(qū)中部的風(fēng)險(xiǎn)上升地區(qū)的空間分布情況更為集中,總體呈現(xiàn)出以大姚縣為中心的,中心高邊緣低的空間分布情況。云南省南部及東部的低海拔山地地區(qū),以及四川西部大部分高海拔地區(qū),地震風(fēng)險(xiǎn)變化指數(shù)則保持一個(gè)較為穩(wěn)定的情況。
圖8 研究區(qū)地貌分布Fig.8 Landform distribution map of Sichuan-Yunnan area
壤塘縣、黑水縣以及汶川縣均屬于高海拔山地地形,但黑水縣的地震風(fēng)險(xiǎn)變化指數(shù)增長明顯高于壤塘縣和汶川縣,主要原因可能是因?yàn)楹谒h的CA變化率更高。同樣,蒼溪縣和名山縣同為低海拔山地,但蒼溪縣的CA變化率比名山縣高,地震風(fēng)險(xiǎn)變化指數(shù)也相應(yīng)更高。此外,還有同屬中海拔山地地區(qū)的青山縣和江城縣,以及同屬中高海拔地區(qū)的美姑縣與昆明市。說明相同地形條件的區(qū)縣,地震風(fēng)險(xiǎn)變化指數(shù)主要受城市空間分布變化影響。
2008年5月12日的汶川地震對川滇地區(qū)造成了嚴(yán)重影響,尤其是震中汶川縣一帶。本次實(shí)驗(yàn)結(jié)果表明,汶川縣城鎮(zhèn)建筑總面積減少43.48 %,城市圖斑數(shù)量減少58.33 %,屬于聚集型減少城市。城市暴露性的急劇減小使得汶川縣將來面臨的地震風(fēng)險(xiǎn)變化指數(shù)減少了44.44 %。而同為極重災(zāi)區(qū)的綿竹縣城市空間分布變化也表現(xiàn)為聚集型減少,對應(yīng)的地震風(fēng)險(xiǎn)變化指數(shù)減少了19.71 %。表明遭受過嚴(yán)重地震的地區(qū),相比于未發(fā)生過地震的地區(qū),地震風(fēng)險(xiǎn)變化指數(shù)降低更明顯。
通過使用GL30數(shù)據(jù)結(jié)合PGA數(shù)據(jù),計(jì)算了2000年至2010年川滇地區(qū)地震風(fēng)險(xiǎn)變化指數(shù),主要可以得出以下結(jié)論:
1)2000~2010年,地震風(fēng)險(xiǎn)變化率與城市聚集性變化呈正相關(guān)關(guān)系。地震風(fēng)險(xiǎn)變化指數(shù)增長明顯的區(qū)縣大多也表現(xiàn)出城市建筑面積增加以及圖斑數(shù)量增加的趨勢,即城市空間分布變化呈離散型增長的地區(qū),地震風(fēng)險(xiǎn)變化指數(shù)增長越明顯。
2)地震風(fēng)險(xiǎn)變化指數(shù)增長地區(qū)大多屬于中高海拔山地。山地相較于其他地形,對地震風(fēng)險(xiǎn)變化指數(shù)的影響更大。相同地形條件的區(qū)縣,地震風(fēng)險(xiǎn)變化指數(shù)主要受城市空間分布變化的影響。
3)遭受過嚴(yán)重地震的地區(qū)城市暴露性急劇減小,相比于未發(fā)生過地震的地區(qū),地震風(fēng)險(xiǎn)變化指數(shù)降低更明顯。