1409號(hào)超強(qiáng)臺(tái)風(fēng)“威馬遜”對(duì)海南橡膠園的風(fēng)害影響評(píng)估

張京紅 張明潔 張亞杰 車(chē)秀芬

(1 海南省氣候中心，?？?570203；2 海南省南海氣象防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室，海口 570203）

0 引言

中國(guó)天然橡膠樹(shù)種植主要分布在海南、云南、廣東、廣西、福建等地，其中海南是中國(guó)主要的天然橡膠產(chǎn)地之一。截至2016年末，海南省橡膠種植面積達(dá)54.09萬(wàn) hm2，收獲面積37.94萬(wàn) hm2，年干膠產(chǎn)量35.14萬(wàn)噸[1]，占全國(guó)產(chǎn)量的43.07%；海南島的橡膠林覆蓋了全島本土面積約16.0%，其生態(tài)服務(wù)功能所產(chǎn)生的效益及其影響巨大。橡膠樹(shù)起源于高溫、靜風(fēng)的巴西亞馬遜河流域，抗風(fēng)能力差[2-3]。海南地區(qū)天然橡膠樹(shù)一年中的高產(chǎn)期在8—10月，基本處于臺(tái)風(fēng)影響活躍期[4]。因此，臺(tái)風(fēng)是海南省橡膠種植的首要?dú)庀鬄?zāi)害，嚴(yán)重影響橡膠樹(shù)的生長(zhǎng)和產(chǎn)量[4-6]。海南島橡膠林遭受熱帶氣旋災(zāi)害風(fēng)險(xiǎn)巨大，且成災(zāi)迅速、損失慘重。以往的橡膠風(fēng)害評(píng)估通常以災(zāi)情實(shí)地調(diào)查、逐級(jí)上報(bào)的方式[7-8]為主，費(fèi)時(shí)費(fèi)力，近年來(lái)遙感和GIS技術(shù)被逐漸應(yīng)用到橡膠林臺(tái)風(fēng)災(zāi)害監(jiān)測(cè)中[9-12]。

為了有效防范和減少橡膠臺(tái)風(fēng)災(zāi)害損失，十分有必要運(yùn)用科學(xué)的技術(shù)方法提高海南橡膠園臺(tái)風(fēng)災(zāi)害影響評(píng)估的及時(shí)性和準(zhǔn)確性。橡膠園臺(tái)風(fēng)災(zāi)害的成災(zāi)因素具有復(fù)雜性和不確定性，而可拓理論具有極大的處理空間數(shù)據(jù)的優(yōu)勢(shì)。因此，本文依據(jù)可拓理論，探索海南橡膠風(fēng)害評(píng)估技術(shù)，形成了一套科學(xué)性、可操作性強(qiáng)的快速及時(shí)評(píng)估方法，可為合理規(guī)劃和科學(xué)栽培橡膠，防災(zāi)減災(zāi)高產(chǎn)穩(wěn)產(chǎn)，增強(qiáng)海南橡膠園應(yīng)對(duì)臺(tái)風(fēng)災(zāi)害的能力提供科學(xué)支撐。

1 數(shù)據(jù)來(lái)源與研究方法

1.1 數(shù)據(jù)來(lái)源

所用資料來(lái)源于海南省氣象臺(tái)的1950—2018年的臺(tái)風(fēng)影響海南期間的國(guó)家基本站和區(qū)域自動(dòng)氣象站的觀測(cè)數(shù)據(jù)；臺(tái)風(fēng)災(zāi)情數(shù)據(jù)來(lái)源于海南省民政廳的1987—2018年臺(tái)風(fēng)災(zāi)情統(tǒng)計(jì)數(shù)據(jù)；天然橡膠資料來(lái)源于《海南省統(tǒng)計(jì)年鑒》[1]以及部分農(nóng)場(chǎng)記錄數(shù)據(jù)；地理信息數(shù)據(jù)來(lái)源于海南省測(cè)繪局的地理信息資料和土地利用類(lèi)型數(shù)據(jù)；橡膠園災(zāi)情資料來(lái)自海南省天然性橡膠產(chǎn)業(yè)集團(tuán)的災(zāi)情調(diào)查數(shù)據(jù)。

1.2 研究方法

橡膠園臺(tái)風(fēng)災(zāi)害成災(zāi)的控制因素多，復(fù)雜而且不確定，而可拓理論是空間數(shù)據(jù)挖掘技術(shù)中的一種新型智能工具[13]，具有處理空間數(shù)據(jù)的優(yōu)勢(shì)，所以選擇可拓理論作為臺(tái)風(fēng)災(zāi)害評(píng)估的數(shù)學(xué)方法?？赏啬Ｐ驮u(píng)估的基本思路是：1）根據(jù)歷史災(zāi)害數(shù)據(jù)集確定各評(píng)估因子特征的經(jīng)典域和節(jié)域；2）確定各特征的關(guān)聯(lián)函數(shù)；3）將待評(píng)估評(píng)價(jià)因子的具體信息數(shù)據(jù)帶入各數(shù)據(jù)類(lèi)中，計(jì)算關(guān)于各特征的關(guān)聯(lián)度，根據(jù)特征的權(quán)重系數(shù)計(jì)算綜合關(guān)聯(lián)度；4）最后與各類(lèi)別集合的關(guān)聯(lián)度大小進(jìn)行比較，關(guān)聯(lián)度越大，與某類(lèi)別的符合程度就越高[14]。

橡膠園臺(tái)風(fēng)災(zāi)害是孕災(zāi)環(huán)境、致災(zāi)因子、承災(zāi)體共同作用的結(jié)果，研究中從這三方面綜合考慮來(lái)確定橡膠園臺(tái)風(fēng)災(zāi)害評(píng)估的因子[15]。研究過(guò)程中，主要以高程、坡度、坡向、坡位等地形因子來(lái)衡量孕災(zāi)環(huán)境對(duì)橡膠臺(tái)風(fēng)災(zāi)害造成的影響。用臺(tái)風(fēng)影響過(guò)程的降水量、降水強(qiáng)度和風(fēng)向、風(fēng)速等氣象因子來(lái)衡量致災(zāi)因子對(duì)橡膠臺(tái)風(fēng)災(zāi)害造成的影響。品種、林齡、種植形式和種植密度、整形修枝和防護(hù)林等的不同，會(huì)造成橡膠園在同等強(qiáng)度的氣象災(zāi)害強(qiáng)度影響下，損失程度的不同。另外，種植面積越大，可能遭受的災(zāi)害損失就越大。因此，用橡膠品種、林齡、種植形式和種植密度、整形修枝和防護(hù)林、種植面積等橡膠因子來(lái)衡量承載體對(duì)橡膠臺(tái)風(fēng)災(zāi)害造成的影響。其中，橡膠因子和地形因子屬于靜態(tài)因子，降水量、降水強(qiáng)度和風(fēng)向、風(fēng)速等氣象因子屬于動(dòng)態(tài)因子。

依據(jù)可拓理論，建立海南橡膠風(fēng)害評(píng)估模型的步驟為：?jiǎn)卧翊_定→影響因素選取→評(píng)估結(jié)果的分級(jí)界定→各因子域值確定→標(biāo)準(zhǔn)物元模型→因子權(quán)重系數(shù)確定→評(píng)估模型。研究中根據(jù)在一定條件下某一段時(shí)間該區(qū)域內(nèi)橡膠風(fēng)害損失程度的大小，將橡膠風(fēng)害評(píng)估結(jié)果表述為 5個(gè)等級(jí)：橡膠風(fēng)害率指標(biāo)＞0.3（5級(jí)，嚴(yán)重）、0.2～0.3（4級(jí)，較嚴(yán)重）、0.1～0.2（3級(jí)，中等）、0.05～0.1（2級(jí)，較輕）、0～0.05（1級(jí)，輕），分別描述橡膠林成災(zāi)的大小。因子權(quán)重系數(shù)是反映各因子影響橡膠林成災(zāi)大小的重要系數(shù)，根據(jù)前人研究成果[16-17]和研究區(qū)歷史災(zāi)害資料，利用熵權(quán)法確定各因子的權(quán)重系數(shù)，取α=（0.297，0.1634，0.5396）表述橡膠綜合指數(shù)、地形綜合指數(shù)、氣象綜合指數(shù)等3個(gè)因子在橡膠林風(fēng)害評(píng)估中的權(quán)重。

2 結(jié)果與分析

2.1 1409號(hào)超強(qiáng)臺(tái)風(fēng)“威馬遜”風(fēng)雨概況

2014年第9號(hào)臺(tái)風(fēng)“威馬遜”于7月12日在西北太平洋洋面上生成，18日15：30（北京時(shí)，下同）在文昌翁田鎮(zhèn)沿海登陸，登陸時(shí)中心附近最大風(fēng)力17級(jí)（60 m/s），中心最低氣壓910 hPa。登陸后穿過(guò)文昌東北角進(jìn)入瓊州海峽，分別于18日19：30和19日07：10在廣東徐聞縣沿海和廣西防城港市沿海再次登陸?！巴R遜”是1949年以來(lái)登陸我國(guó)最強(qiáng)的臺(tái)風(fēng)，也是華南三?。▍^(qū)）有氣象記錄以來(lái)登陸最強(qiáng)的臺(tái)風(fēng)。

受“威馬遜”影響，海南普降暴雨和大暴雨，部分地區(qū)出現(xiàn)特大暴雨。18日08時(shí)—19日07時(shí)，全島共有39個(gè)鄉(xiāng)鎮(zhèn)雨量超過(guò)300 mm，19個(gè)鄉(xiāng)鎮(zhèn)雨量超過(guò)400 mm，7個(gè)鄉(xiāng)鎮(zhèn)雨量超過(guò)500 mm，昌江有2個(gè)鄉(xiāng)鎮(zhèn)雨量達(dá)到600 mm以上（昌化鎮(zhèn)678.2 mm和七叉鎮(zhèn)630.6 mm，圖1）。海南島東北部陸地普遍出現(xiàn)平均風(fēng)10～12級(jí)，陣風(fēng)13～16級(jí)，最大風(fēng)力出現(xiàn)在文昌翁田鎮(zhèn)，測(cè)得陣風(fēng)17級(jí)（58.8 m/s），平均風(fēng)12級(jí)（36.2 m/s）。海南島東部海面浮標(biāo)站和文昌七洲列島均測(cè)得最大陣風(fēng)17級(jí)以上（分別為74.1和72.4 m/s），另外文昌七洲列島測(cè)得最大平均風(fēng)17級(jí)（58.7 m/s）。

2.2 評(píng)估模型中降水和風(fēng)數(shù)據(jù)的處理分析

評(píng)估中用到的氣象數(shù)據(jù)來(lái)自海南省氣象局布設(shè)的自動(dòng)氣象站網(wǎng)的觀測(cè)數(shù)據(jù)，海南島目前有400多套自動(dòng)氣象站，站點(diǎn)密集且分布均勻，運(yùn)行穩(wěn)定，能夠提供穩(wěn)定且精確度較高的觀測(cè)數(shù)據(jù)，可以滿(mǎn)足橡膠風(fēng)害實(shí)時(shí)評(píng)估的需求。研究中需將數(shù)據(jù)處理為GIS可讀的格式，即在GIS中將降水、風(fēng)速等點(diǎn)數(shù)據(jù)文件采用通用克里金法進(jìn)行插值，獲取與靜態(tài)數(shù)據(jù)具有相同投影信息、相同分辨率的柵格文件，輸入模型進(jìn)行相應(yīng)的運(yùn)算。

圖1 2014年7月18日08時(shí)—19日07時(shí)海南各地降雨量（單位：mm）（圖中紅色數(shù)字為城區(qū)測(cè)站雨量，黑色數(shù)字為鄉(xiāng)鎮(zhèn)自動(dòng)站雨量）Fig. 1 Rainfall in Hainan from 08:00 to 07:00 on 18 July 2014 (unit: mm)(The red figures represent the rainfall of urban stations,and the black ones represent the rainfall of automatic stations)

2.3 “威馬遜”對(duì)橡膠園影響的評(píng)估

臺(tái)風(fēng)災(zāi)害過(guò)程中的大風(fēng)和降水是主要致災(zāi)因子，根據(jù)建立的海南橡膠園臺(tái)風(fēng)災(zāi)害評(píng)估模型計(jì)算結(jié)果分析，1409號(hào)強(qiáng)臺(tái)風(fēng)“威馬遜”影響過(guò)程（2014年7月17日08時(shí)—19日08時(shí)）對(duì)海南全島橡膠園都造成了影響（圖2），五指山以北地區(qū)受到的災(zāi)害影響重于五指山以南地區(qū)，災(zāi)害等級(jí)達(dá)到3級(jí)及以上，其中?？?、澄邁、臨高、儋州、白沙、昌江的大部分地區(qū)災(zāi)害等級(jí)達(dá)到五級(jí)，五指山以南和東南部沿海地區(qū)災(zāi)害損失等級(jí)主要為1～2級(jí)。由表1可知，評(píng)估結(jié)果1～5級(jí)單元格數(shù)占海南橡膠區(qū)域總單元格數(shù)的比例分別為19.81%、26.41%、14.35%、19.10%、20.33%。

在強(qiáng)臺(tái)風(fēng)“威馬遜”登陸后，海南天然橡膠產(chǎn)業(yè)集團(tuán)股份有限公司對(duì)其各分公司的橡膠園橡膠損失情況進(jìn)行了實(shí)際調(diào)查。根據(jù)其調(diào)查數(shù)據(jù)，海南島橡膠北部地區(qū)受災(zāi)嚴(yán)重，實(shí)際風(fēng)害率在20%以上（圖3），其中東北部的臨高、澄邁、海口、定安、文昌、瓊海的部分地區(qū)≥30%，五指山以北的其他地區(qū)風(fēng)害率在10%～20%，五指山以南地區(qū)風(fēng)害率在10%以下。表明模型評(píng)估結(jié)果與災(zāi)情實(shí)際調(diào)查結(jié)果總體趨勢(shì)一致，即東部和北部地區(qū)遭受的風(fēng)害最為嚴(yán)重，五指山以南地區(qū)災(zāi)情相對(duì)較輕，且從圖中可以看出，在災(zāi)害損失相對(duì)較輕的區(qū)域，零星局部地區(qū)損失較嚴(yán)重，顯示了模型評(píng)估結(jié)果的局地差異性，說(shuō)明建立的評(píng)估模型能夠客觀精細(xì)的對(duì)橡膠園臺(tái)風(fēng)災(zāi)害風(fēng)災(zāi)分布及受災(zāi)程度進(jìn)行評(píng)估?；谝陨戏治觯?409號(hào)強(qiáng)臺(tái)風(fēng)“威馬遜”對(duì)海南橡膠園的災(zāi)害影響經(jīng)綜合評(píng)估屬于重度災(zāi)害。

圖2 海南島橡膠園風(fēng)害評(píng)估結(jié)果Fig. 2 Wind damage assessment results of rubber plantation in Hainan Islands

表1 評(píng)估結(jié)果1～5級(jí)單元格數(shù)占區(qū)域總單元格數(shù)的比例

圖3 海南島橡膠園實(shí)際風(fēng)害率分布Fig. 3 Distribution of actual wind damage rate in rubber plantation in Hainan Islands

3 結(jié)論和建議

本研究根據(jù)海南橡膠園臺(tái)風(fēng)災(zāi)害形成的機(jī)理，從孕災(zāi)環(huán)境、致災(zāi)因子、承載體三個(gè)方面，確定橡膠園臺(tái)風(fēng)災(zāi)害評(píng)估的靜態(tài)和動(dòng)態(tài)因子，利用熵權(quán)法確定各評(píng)估因子的權(quán)重，依據(jù)可拓理論建立了海南橡膠園臺(tái)風(fēng)災(zāi)害評(píng)估模型。探索出了一套以可拓模型理論為基礎(chǔ)，從“單元格確定—影響因素選取—評(píng)估結(jié)果的分級(jí)界定—各因子域值確定—標(biāo)準(zhǔn)物元模型—因子權(quán)重系數(shù)確定—評(píng)估模型”的評(píng)估方法，用可拓集合的關(guān)聯(lián)函數(shù)值—關(guān)聯(lián)度大小描述各種特征參數(shù)與所研究對(duì)象的從屬關(guān)系，把屬于或不屬于的定性描述擴(kuò)展為定量描述，使評(píng)估從定性走向定量化。并對(duì)1409號(hào)強(qiáng)臺(tái)風(fēng)“威馬遜”造成的災(zāi)害損失進(jìn)行了評(píng)估，評(píng)估結(jié)果與實(shí)際災(zāi)情較一致，說(shuō)明研究依據(jù)可拓理論建立的橡膠園臺(tái)風(fēng)災(zāi)害評(píng)估模型是可行的。

進(jìn)一步分析發(fā)現(xiàn)，中部地區(qū)的災(zāi)害評(píng)估等級(jí)重于橡膠園風(fēng)害實(shí)際調(diào)查結(jié)果，分析原因，可能是與橡膠園的地形、品種和種植形式等因素有關(guān)，后續(xù)還將在更多的橡膠園臺(tái)風(fēng)災(zāi)害實(shí)例評(píng)估中，在獲得更加翔實(shí)的橡膠信息（包括海南橡膠園的種植面積、品種抗風(fēng)性、種植形式和種植密度、整形修枝、防護(hù)林等）并及時(shí)更新的基礎(chǔ)上不斷改進(jìn)評(píng)估模型，提高評(píng)估結(jié)果的精確度。