李孝波,常晁瑜,蔡曉光,吳義文,周興浩
(1.防災(zāi)科技學(xué)院,河北 三河 065201;2.河北省地震災(zāi)害防御與風(fēng)險(xiǎn)評(píng)價(jià)重點(diǎn)實(shí)驗(yàn)室,河北 三河 065201)
據(jù)中國(guó)地震臺(tái)網(wǎng)資料,2022年9月5日12時(shí)52分在四川省甘孜州瀘定縣發(fā)生MS6.8地震,震中位于瀘定縣磨西鎮(zhèn)海螺溝冰川森林公園內(nèi)(102.08°E,29.59°N),震源深度16km。地震震中距瀘定縣39km,距康定市47km,距石棉縣48km,距漢源縣62km,距滎經(jīng)縣78km,距成都市226km。地震誘發(fā)了大量的崩塌、滑坡以及滾石等地質(zhì)災(zāi)害,造成了嚴(yán)重的財(cái)產(chǎn)損失和人員傷亡[1-2]。從圖1可知:瀘定6.8級(jí)地震最高烈度為Ⅸ度(9度),等震線長(zhǎng)軸總體呈北西走向,長(zhǎng)軸195km,短軸112km,Ⅵ度(6度)區(qū)及以上面積19 089km2,主要涉及四川省3個(gè)市州12個(gè)縣(市),82個(gè)鄉(xiāng)鎮(zhèn)(街道)。其中,Ⅸ度(9度)區(qū)面積280km2,主要包含甘孜州瀘定縣磨西鎮(zhèn)、得妥鎮(zhèn)、燕子溝鎮(zhèn)、德威鎮(zhèn);雅安市石棉縣王崗坪彝族藏族鄉(xiāng)、草科藏族鄉(xiāng)、新民藏族彝族鄉(xiāng),共計(jì)7個(gè)鄉(xiāng)鎮(zhèn)。
圖1 四川瀘定6.8級(jí)地震烈度圖[4]Fig.1 Earthquake intensity map of Luding,Sichuan M S6.8 earthquake in 2022[4]
為更加全面地了解震害現(xiàn)象,深入探究震害機(jī)理,防災(zāi)科技學(xué)院瀘定6.8級(jí)地震科考隊(duì)于9月24日抵達(dá)震害現(xiàn)場(chǎng),深入磨西鎮(zhèn)、得妥鎮(zhèn)、王崗坪鄉(xiāng)以及草科鄉(xiāng)等震害較重鄉(xiāng)鎮(zhèn)開展了12天的震害調(diào)查與現(xiàn)場(chǎng)測(cè)試工作,取得了大量的第一手震害資料。位于磨西臺(tái)地之上的磨西鎮(zhèn),距震中7km 左右,地震動(dòng)強(qiáng)度大,房屋破壞嚴(yán)重(圖2),臺(tái)地兩側(cè)邊緣崩塌、滑坡等地質(zhì)災(zāi)害發(fā)育。圖3給出了瀘定6.8級(jí)地震作用下磨西臺(tái)地邊緣發(fā)育的滑坡、崩塌災(zāi)害,整體呈現(xiàn)了分布面積較大、順坡向發(fā)育和群發(fā)性明顯的特征。然而,由于滑動(dòng)面埋深較淺,滑體厚度不大,滑坡規(guī)模整體以中小型為主。結(jié)合范宣梅等[1]、鐵永波等[3]的研究成果,我們認(rèn)為磨西臺(tái)地震害之所以嚴(yán)重,除與地形地貌、地震動(dòng)強(qiáng)度相關(guān)以外,特殊的地層巖性也是一個(gè)不可忽視的影響因素。
圖2 磨西鎮(zhèn)的房屋震害Fig.2 Photos of earthquake damage of buildings in Moxi town
圖3 磨西臺(tái)地邊緣的滑坡災(zāi)害Fig.3 Photos of landslides on the edge of the Moxi Platform
磨西臺(tái)地長(zhǎng)約10km、寬1~2km,主要由燕子溝、磨子溝和海螺溝的冰川、冰水和洪積泥石流沉積而成,形成演化歷史較為復(fù)雜[5]。臺(tái)地屬于典型的冰川地貌,主體地層為倒數(shù)第二次冰期的冰川堆積物,其次為沖洪積物和老泥石流堆積物,厚120m左右,臺(tái)地前緣有較為明顯的層理現(xiàn)象(圖4)[1,5,6]。冰川堆積物的礫石成分主要為花崗巖、花崗閃長(zhǎng)巖等,磨圓度和分選性均較差,粒徑一般20~80cm,隨高程變化而有所差異[6,7]。長(zhǎng)期的河流地質(zhì)作用,使臺(tái)地邊緣逐漸被切割成高陡斜坡(圖4),在磨西臺(tái)地以東的雅家梗流域形成4級(jí)階地,沉積物包含礫石層、含礫石層、淤泥層等13層,新興鄉(xiāng)東側(cè)的堆積物分別高出河床5~10m、35m、65m和85m;磨西臺(tái)地以西燕子溝的多級(jí)河流階地主要保存在藥王廟和蔡陽坪附近,藥王廟附近發(fā)育的4級(jí)階地分別高出河床5~10m、30m、50m和100m[7,8]?,F(xiàn)有研究成果表明:磨西臺(tái)地巖土體的膠結(jié)性良好,但在強(qiáng)烈的地震作用下仍然容易產(chǎn)生破壞,特別是臺(tái)地兩側(cè)的臨空邊坡極易發(fā)生垮塌,產(chǎn)生滑坡、崩塌等地質(zhì)災(zāi)害[3]。
圖4 磨西臺(tái)地前緣高陡斜坡與地層層理現(xiàn)象Fig.4 High steep slope and str atigraphic bedding in the front edge of the Moxi Platform
為進(jìn)一步探究磨西臺(tái)地巖土體動(dòng)力響應(yīng)規(guī)律,填補(bǔ)磨西臺(tái)地場(chǎng)地地震效應(yīng)研究空白,本文基于瀘定6.8級(jí)地震科考獲取的地脈動(dòng)數(shù)據(jù),采用水平與垂直譜比法(Horizontal to Vertical Spectral Ratio,簡(jiǎn)稱HVSR法)[9],分析磨西臺(tái)地巖土體的頻譜、方向性特征,探究磨西臺(tái)地場(chǎng)地地震效應(yīng)規(guī)律,以期為磨西臺(tái)地所在鄉(xiāng)鎮(zhèn)災(zāi)后的規(guī)劃與重建提供參考。
地脈動(dòng)是地面的一種連續(xù)的非重復(fù)性的隨機(jī)波動(dòng),主要由氣象、海洋、地殼構(gòu)造活動(dòng)的自然力和交通等人為因素所引起的地球表面固有的微 弱 振 動(dòng)[10]。 依 據(jù) SESAME(Site EffectS assessment using AMbient Excitations)使 用 指南[11]和 地 基 動(dòng) 力 特 性 測(cè) 試 規(guī) 范(GB/T 50269-2015)[10],我們于9月29日至10月1日在磨西臺(tái)地先后開展了3次地脈動(dòng)測(cè)試工作,測(cè)試時(shí)間主要集中于晚上10點(diǎn)至次日早上8點(diǎn)。測(cè)試儀器采用Tromino數(shù)字式地脈動(dòng)儀(具體參數(shù)詳 見https://moho.world/en/tromino/),采 樣頻率128Hz,采集時(shí)間20min。
圖5給出了21個(gè)地脈動(dòng)測(cè)點(diǎn)的空間分布位置,可以看出測(cè)點(diǎn)的展布方向近似平行于磨西臺(tái)地的延伸方向,測(cè)點(diǎn)間距500m左右,測(cè)試總長(zhǎng)度近10km。表1列出了各點(diǎn)的具體地理位置,從M01至M21各點(diǎn)的高程逐漸增加,前后高差累積640m。M01位于磨西臺(tái)地最前端(龔家河電站后面),高程1310m,上覆土層很薄,其后各點(diǎn)均有一定厚度的土層覆蓋。為保證地脈動(dòng)數(shù)據(jù)質(zhì)量的可靠性,數(shù)據(jù)采集前先平整場(chǎng)地,保證測(cè)試儀器與地面接觸良好,以同時(shí)獲取東西(EW)、南北(NS)和豎直(UD)3個(gè)方向的地脈動(dòng)數(shù)據(jù)。
表1 地脈動(dòng)測(cè)點(diǎn)信息Tab.1 Information of microtremor measuring points
圖5 地脈動(dòng)測(cè)點(diǎn)的空間分布Fig.5 Spatial distribution of microtr emor measuring points
Kanai等[12]最早將地脈動(dòng)用于場(chǎng)地地震效應(yīng)研究,Nakamura[9]提出HVSR法后則更加廣泛應(yīng)用于場(chǎng)地地震反應(yīng)分析領(lǐng)域。HVSR法是一種非參考場(chǎng)地法,其最大的優(yōu)點(diǎn)是在使用過程中不受參考基巖臺(tái)站的限制。HVSR法的提出主要基于兩點(diǎn)假設(shè)[9,13]:(1)基巖或堅(jiān)硬土層上,震動(dòng)水平分量與垂直分量傅氏譜的比值近似為1;(2)震動(dòng)垂直分量受局部場(chǎng)地條件的影響不大。一般認(rèn)為,HVSR法能較準(zhǔn)確地估計(jì)場(chǎng)地卓越頻率,但對(duì)場(chǎng)地放大效應(yīng)有所低估[11,13]。
地脈動(dòng)數(shù)據(jù)獲取后,基于HVSR法,采用開源軟件Geopsy(http://www.geopsy.org)對(duì)其進(jìn)行處理與分析,主要包括數(shù)據(jù)選取、帶通濾波、傅里葉譜計(jì)算、平滑處理以及H/V值計(jì)算等5個(gè)步驟[11]。
(1)數(shù)據(jù)選取:設(shè)定數(shù)據(jù)窗時(shí)長(zhǎng)20s,同時(shí)截取EW、NS、UD 3個(gè)方向的地脈動(dòng)數(shù)據(jù),為保證數(shù)據(jù)窗內(nèi)信號(hào)平穩(wěn),每個(gè)時(shí)窗內(nèi)3個(gè)方向的數(shù)據(jù)均需滿足 0.5<STA/LTA<2.0的判別條件(STA/LTA,Short-Term Average and Long-Term Average ratio,數(shù)據(jù)長(zhǎng)短時(shí)比值),本文取STA=1.00s,LTA=20.00s,各20s時(shí)窗內(nèi)的數(shù)據(jù)允許5%的重疊。
(2)帶通濾波:為保證分析結(jié)果的可靠性,選用4階Butterworth帶通濾波函數(shù)對(duì)每個(gè)數(shù)據(jù)時(shí)窗內(nèi)的原始地脈動(dòng)數(shù)據(jù)進(jìn)行濾波處理,濾波帶寬0.05~20.00Hz,濾波后的數(shù)據(jù)再次進(jìn)行0.5<STA/LTA<2.0條件的檢驗(yàn)。
(3)傅里葉譜計(jì)算:濾波處理完成后,對(duì)每個(gè)數(shù)據(jù)時(shí)窗內(nèi)的數(shù)據(jù)進(jìn)行快速傅里葉變換(FFT),得出EW、NS、UD 3個(gè)方向的傅里葉譜值,完成地脈動(dòng)數(shù)據(jù)從時(shí)域到頻域的轉(zhuǎn)換。
(4)平滑處理:采用Konno和Ohmachi函數(shù)[14,15],對(duì)每個(gè)傅里葉譜窗內(nèi)的數(shù)據(jù)進(jìn)行平滑處理,帶寬系數(shù)b取為40。
(5)H/V值計(jì)算:基于HVSR法分別計(jì)算每個(gè)傅里葉譜窗的H/V值,得出HVSR曲線,其中水平方向(H向)的數(shù)值取EW、NS向數(shù)值的均方根值,每個(gè)測(cè)點(diǎn)的HVSR曲線取各傅里葉譜窗HVSR曲線的平均值。
HVSR曲線的形狀與場(chǎng)地條件相關(guān),常有單峰、雙峰以及多峰等多種形式。不同形狀的HVSR曲線,擁有不同的峰值頻率f0和峰值放大系數(shù)A0。各測(cè)點(diǎn)峰值頻率f0能否作為該點(diǎn)土層的卓越頻率fd,常需進(jìn)行如下5個(gè)條件的驗(yàn)證[11]:(1)A0>2;(2)在[f0/4,f0]頻率范圍內(nèi),至少存在一個(gè)H/V值(AH/V)小于A0/2;(3)在[f0,4f0]頻率范圍內(nèi),同樣至少存在一個(gè)H/V值(AH/V)小于A0/2;(4)在峰值頻率f0所處的頻率范圍內(nèi),峰值頻率f0的標(biāo)準(zhǔn)差σf應(yīng)小于表2中規(guī)定的閾值ε(f0);(5)在峰值頻率f0所處的頻率范圍內(nèi),峰值放大系數(shù)A0的標(biāo)準(zhǔn)差σA同樣應(yīng)小于表2中規(guī)定的閾值θ(f0)。
表2 穩(wěn)定性條件的閾值[11]Tab.2 Threshold values for stability conditions[11]
一般而言,如果一個(gè)測(cè)點(diǎn)的HVSR曲線能滿足上述5個(gè)條件中的4個(gè),就可將該測(cè)點(diǎn)的f0值作為該處土層卓越頻率fd值的可靠估計(jì)。
HVSR法除能較準(zhǔn)確地估計(jì)場(chǎng)地卓越頻率外,還能進(jìn)行場(chǎng)地地震效應(yīng)優(yōu)勢(shì)方向的研究。意大利學(xué)者Del Gaudio教授最早基于HVSR法開展地脈動(dòng)方向性特征研究,并通過與實(shí)際強(qiáng)震觀測(cè)數(shù)據(jù)結(jié)果的對(duì)比,認(rèn)為根據(jù)場(chǎng)地地脈動(dòng)的HVSR分析可以很好地辨識(shí)出場(chǎng)地地震效應(yīng)的優(yōu)勢(shì)方向[16]。此后,眾多學(xué)者的研究成果也都證實(shí)了這個(gè)結(jié)論[17-22]。
基于地脈動(dòng)開展場(chǎng)地地震效應(yīng)方向性研究的流程,除在水平方向(H向)數(shù)值的計(jì)算上有所差別外,其余均與前述步驟相同。依據(jù)SESAME使用指南[11],方向性研究中水平方向(H向)的數(shù)值取EW、NS向在不同方位角上的分量之和,方位角范圍常取0°~180°(與180°~360°呈對(duì)稱關(guān)系),計(jì)算時(shí)以10°間隔逐級(jí)取值。每個(gè)測(cè)點(diǎn)的方向性由H/V值與頻率、方位角的相互關(guān)系來體現(xiàn)。顯著的方向性特征常需滿足以下兩條標(biāo)準(zhǔn)[16,23]:(1)H/V值的最大值大于2;(2)同一頻率下,H/V值最大值與最小值(通常在近似正交的方向上?。┑谋戎荡笥?.5。
M01~M21的HVSR曲線如圖6所示??梢钥闯?,在0.2~15Hz的頻率范圍內(nèi),各測(cè)點(diǎn)的HVSR曲線以單峰、雙峰型為主,少量為多峰或近似平坦型。各測(cè)點(diǎn)都具有較為明顯的峰值頻帶,絕大部分測(cè)點(diǎn)的峰值放大系數(shù)大于2,土層或地形放大效應(yīng)較為顯著。
圖6 M 01~M 21的HVSR曲線(黑色粗實(shí)線為平均的HVSR曲線,黑色細(xì)虛線為HVSR曲線乘以或除以1倍標(biāo)準(zhǔn)差)Fig.6 HVSR curves for microtremor measuring points M 01~M 21(The thick black solid line is the average HVSR curve,the thin black dashed line is the HVSR curve multiplied or divided by 1-time standard deviation)
表3統(tǒng)計(jì)出了M01~M21的峰值頻率f0,除M01、M11和M13的峰值頻率較高外,其余各測(cè)點(diǎn)的峰值頻率都在1Hz左右,集中于中頻段(0.80~1.34Hz),受地層巖性和覆蓋層厚度的共同控制作用明顯。M01測(cè)點(diǎn)由于位于基巖表面,上覆土層很薄,峰值頻率最高達(dá)到了7.21Hz;M11測(cè)點(diǎn)的HVSR曲線較為平緩,峰值頻率不明顯,顯示下伏地層的結(jié)構(gòu)較為均勻;M13測(cè)點(diǎn)的HVSR曲線為多峰型,峰值頻率3.11Hz在第二個(gè)峰值點(diǎn)取得,第一個(gè)峰值點(diǎn)對(duì)應(yīng)的頻率為1.19Hz,介于0.80~1.34Hz。依據(jù)SESAME使用指南[11]判定土層卓越頻率,結(jié)果表明:M01~M05、M10、M14、M16~M18以及M21等11個(gè)測(cè)點(diǎn)的峰值頻率可作為測(cè)點(diǎn)處地層卓越頻率的可靠估計(jì),若不考慮基巖場(chǎng)地測(cè)點(diǎn)M01,則可得出磨西臺(tái)地場(chǎng)地土層卓越頻率的范圍為0.80~1.34Hz。
放大系數(shù)方面,除M11的峰值放大系數(shù)小于2以外,其余各測(cè)點(diǎn)的峰值放大系數(shù)都大于2。放大系數(shù)與測(cè)點(diǎn)分布位置的相關(guān)性不強(qiáng),再次體現(xiàn)了各測(cè)點(diǎn)場(chǎng)地地震效應(yīng)受地層巖性和覆蓋層厚度共同控制的特點(diǎn)。表3的統(tǒng)計(jì)結(jié)果表明:M03、M05~M10、M13、M15、M17以及M19~M21等13個(gè)測(cè)點(diǎn)的峰值放大系數(shù)位于2.25~4.00;M01、M04、M12、M14、M16以及M18等6個(gè)測(cè)點(diǎn)的峰值放大系數(shù)位于4.30~7.00;M02測(cè)點(diǎn)的峰值放大系數(shù)最大,達(dá)到了14.49。需要說明的是,M02測(cè)點(diǎn)的峰值放大系數(shù)之所以超過10,主要是由距測(cè)點(diǎn)50m左右的電機(jī)振動(dòng)造成的。為保障災(zāi)區(qū)應(yīng)急信號(hào)塔的正常運(yùn)行,地脈動(dòng)測(cè)試時(shí)該電機(jī)一直處于工作狀態(tài),從而導(dǎo)致M02測(cè)點(diǎn)的峰值放大系數(shù)較其他測(cè)點(diǎn)都明顯偏大。此外,M01測(cè)點(diǎn)的峰值放大系數(shù)5.76,體現(xiàn)出了明顯的地形放大效應(yīng);M11測(cè)點(diǎn)的HVSR曲線近似平坦,0.2~15Hz頻率范圍內(nèi)的H/V值都小于2,峰值放大系數(shù)只有1.34,推測(cè)與該測(cè)點(diǎn)處的地層結(jié)構(gòu)相關(guān),但具體原因還需進(jìn)一步深入研究??偟脕碚f,磨西臺(tái)地場(chǎng)地放大效應(yīng)較為明顯,與場(chǎng)地土層卓越頻率范圍相對(duì)應(yīng)的峰值放大系數(shù)為2.69~6.97。
表3 各測(cè)點(diǎn)峰值頻率f 0、峰值放大系數(shù)A 0以及地層卓越頻率f d 的估計(jì)Tab.3 The peak frequency f 0,the peak amplification factor A 0 and the predominant frequency of formation f d for each microtr emor measuring point
圖7給出了M01~M21的H/V值與頻率、方位角的變化關(guān)系。為能更好地體現(xiàn)方向性特征,圖中M11測(cè)點(diǎn)的頻率范圍為0.3~15Hz,M19測(cè)點(diǎn)的頻率范圍為0.4~15Hz,其余各測(cè)點(diǎn)的頻率范圍均為0.2~15Hz??梢钥闯?,各圖顏色的變化規(guī)律與圖6各點(diǎn)HVSR曲線的變化趨勢(shì)大體一致,顏色越紅H/V的數(shù)值越大,對(duì)應(yīng)的頻率越接近峰值頻率,對(duì)應(yīng)的方位角越接近場(chǎng)地地震效應(yīng)的優(yōu)勢(shì)方向。整體來看,M01~M05、M07~M10、M12~M16以及M18~M21等18個(gè)測(cè)點(diǎn)都表現(xiàn)出了良好的方向性特征,優(yōu)勢(shì)方向范圍較易識(shí)別,而M06、M11以及M17等3個(gè)測(cè)點(diǎn)由于在0°~180°范圍內(nèi)全部或多處放大,方向性特征表現(xiàn)不突出,難以確定場(chǎng)地地震效應(yīng)的優(yōu)勢(shì)方向。
圖7 M 01~M 21的H/V值與頻率、方位角的變化關(guān)系Fig.7 The relationship between H/V val ues wit h f r equency and azim uth for m icrotrem or m easuring points M 01~M 21
依據(jù)Del Gaudio等[16,23]提出的判定標(biāo)準(zhǔn),18個(gè)測(cè)點(diǎn)優(yōu)勢(shì)方向的判定結(jié)果列于表4。判定結(jié)果表明:除M07、M08、M12、M14、M16以及M18等6個(gè)測(cè)點(diǎn)以外,其余12個(gè)測(cè)點(diǎn)都具有顯著的方向性特征,優(yōu)勢(shì)方位角主要集中于100°~150°和140°~210°兩個(gè)區(qū)間,整體具有隨高程增加而逐漸增大的趨勢(shì),體現(xiàn)出了地層巖性差異對(duì)場(chǎng)地地震效應(yīng)優(yōu)勢(shì)方向的影響。結(jié)合表3的結(jié)論,表4中M01~M05、M10以及M21等7個(gè)測(cè)點(diǎn)都是地層卓越頻率的可靠估計(jì)點(diǎn),M01~M05、M10的優(yōu)勢(shì)方位角為100°~150°,M21的優(yōu)勢(shì)方位角為150°~190°。因此,可以得出與磨西臺(tái)地場(chǎng)地土層卓越頻率相對(duì)應(yīng)的優(yōu)勢(shì)方向范圍為100°~150°和150°~190°。
表4 12個(gè)測(cè)點(diǎn)的優(yōu)勢(shì)方向Tab.4 The dom inant direction of 12 microt rem or m easuring points
本文基于瀘定磨西臺(tái)地的地脈動(dòng)數(shù)據(jù),采用水平與垂直譜比法,從峰值頻率、放大系數(shù)以及優(yōu)勢(shì)方向等3個(gè)方面,探究了磨西臺(tái)地場(chǎng)地地震效應(yīng)特征,得出了一些初步的規(guī)律性認(rèn)識(shí)。總得來說,磨西臺(tái)地場(chǎng)地地震效應(yīng)特征較為明顯,土層的卓越頻率集中于0.80~1.34Hz,與卓越頻率相對(duì)應(yīng)的峰值放大系數(shù)為2.69~6.97,100°~150°和150°~190°是磨西臺(tái)地場(chǎng)地地震效應(yīng)優(yōu)勢(shì)方向的兩個(gè)區(qū)間。磨西臺(tái)地場(chǎng)地地震效應(yīng)特征參數(shù)的確定,填補(bǔ)了磨西臺(tái)地場(chǎng)地地震效應(yīng)研究的空白,為深入開展冰川堆積物地層動(dòng)力響應(yīng)研究奠定了基礎(chǔ),也可為磨西臺(tái)地所在鄉(xiāng)鎮(zhèn)災(zāi)后的規(guī)劃與重建提供技術(shù)支撐。
然而,囿于基礎(chǔ)數(shù)據(jù)資料的缺乏,本文的研究工作還有待進(jìn)一步完善。例如,現(xiàn)有工作只沿著磨西臺(tái)地的中軸線開展了地脈動(dòng)測(cè)試,且各點(diǎn)測(cè)試間距較大,接下來應(yīng)開展更大范圍、更小間距的地脈動(dòng)測(cè)試工作,進(jìn)行整個(gè)磨西臺(tái)地的地震效應(yīng)分析,深入探究冰川堆積物地層的動(dòng)力響應(yīng)規(guī)律,更加科學(xué)合理的指導(dǎo)同類型場(chǎng)地城鎮(zhèn)的防災(zāi)減災(zāi)規(guī)劃與工程建設(shè)。此外,磨西臺(tái)地作為一個(gè)高120m左右的臺(tái)地,臺(tái)地邊緣高陡斜坡的地形效應(yīng)也不能忽視,后續(xù)還應(yīng)在臺(tái)地兩側(cè)沖溝場(chǎng)地開展地脈動(dòng)測(cè)試,進(jìn)一步探究臺(tái)地的地形效應(yīng),加深對(duì)臺(tái)地地震響應(yīng)規(guī)律的認(rèn)識(shí),更好地服務(wù)防災(zāi)減災(zāi)事業(yè)的發(fā)展。