• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    熱帶氣旋災(zāi)害風(fēng)險評估現(xiàn)狀綜述*

    2017-10-11 09:11:09陳文方端義宏方佳毅石先武任福民
    災(zāi)害學(xué) 2017年4期
    關(guān)鍵詞:風(fēng)暴潮氣旋脆弱性

    陳文方,端義宏,陸 逸,方佳毅,石先武,任福民

    (1.中國氣象科學(xué)研究院災(zāi)害天氣國家重點實驗室,北京 100081;2.中國氣象局上海臺風(fēng)研究所,上海 200030;3.北京師范大學(xué) 減災(zāi)與應(yīng)急管理研究院,北京 100875;4.英國南安普頓大學(xué) 環(huán)境工程學(xué)院,南安普頓 英國 SO17 1BJ;5.國家海洋局海洋減災(zāi)中心,北京 100194)

    熱帶氣旋災(zāi)害風(fēng)險評估現(xiàn)狀綜述*

    陳文方1,端義宏1,陸 逸2,方佳毅3,4,石先武5,任福民1

    (1.中國氣象科學(xué)研究院災(zāi)害天氣國家重點實驗室,北京 100081;2.中國氣象局上海臺風(fēng)研究所,上海 200030;3.北京師范大學(xué) 減災(zāi)與應(yīng)急管理研究院,北京 100875;4.英國南安普頓大學(xué) 環(huán)境工程學(xué)院,南安普頓 英國 SO17 1BJ;5.國家海洋局海洋減災(zāi)中心,北京 100194)

    全球熱帶氣旋災(zāi)害風(fēng)險隨著社會經(jīng)濟的發(fā)展、沿海地區(qū)人口和財富的加速增長與聚集而持續(xù)上升,針對沿海地區(qū)熱帶氣旋災(zāi)害風(fēng)險開展精細化的定量評估是實現(xiàn)前瞻性災(zāi)害風(fēng)險管理的前提。該文從面向災(zāi)害管理部門開展熱帶氣旋災(zāi)害風(fēng)險管理的角度出發(fā),分別對國際上熱帶氣旋風(fēng)險評估中致災(zāi)因子危險性評估和承災(zāi)體脆弱性評估這兩個核心內(nèi)容的研究進展進行綜述,分析了其各自的現(xiàn)狀、存在的問題及其挑戰(zhàn)。然后對我國熱帶氣旋(臺風(fēng))災(zāi)害風(fēng)險管理體系中風(fēng)險定量評估的研究現(xiàn)狀進行了討論,針對我國熱帶氣旋風(fēng)險模型的研究現(xiàn)狀提出建議,為未來進一步開展沿海地區(qū)臺風(fēng)災(zāi)害風(fēng)險的定量評估研究提供了參考。

    熱帶氣旋;風(fēng)險評估;災(zāi)害;致災(zāi)因子;脆弱性

    全球范圍內(nèi),熱帶氣旋是給人類社會帶來損失最大的自然災(zāi)害之一。據(jù)統(tǒng)計,最近10年全球熱帶氣旋災(zāi)害平均每年直接經(jīng)濟損失約為550億美元(2006-2015,2015年價格),超過洪水、地震等其他自然災(zāi)害[1],而中國平均每年因熱帶氣旋直接經(jīng)濟損失約為90億美元,約為全球損失的15%。并且,沿海地區(qū)人口和財富的快速增長促使熱帶氣旋造成的經(jīng)濟損失持續(xù)增長,這些地區(qū)的人口和社會經(jīng)濟都將面臨著越來越高的災(zāi)害風(fēng)險[2-3],為此,熱帶氣旋災(zāi)害風(fēng)險的評估與管理受到了科學(xué)界越來越廣泛的關(guān)注[4-5]。

    如圖1所示,熱帶氣旋風(fēng)險的大小由熱帶氣旋致災(zāi)因子(大風(fēng)、降雨、風(fēng)暴潮等)本身的強度(即危險性)、承災(zāi)體(人口、房屋、財產(chǎn)等)的數(shù)量或價值(即暴露度)以及承災(zāi)體對致災(zāi)因子強度的敏感度(即脆弱性)等三個因素共同決定[6]?,F(xiàn)代災(zāi)害風(fēng)險管理迫切需要系統(tǒng)化、精細化和定量化的災(zāi)害風(fēng)險模型,即通過對熱帶氣旋致災(zāi)因子危險性、特定承災(zāi)體的暴露度和脆弱性分別開展定量的模擬和評估,獲得精細化的定量的災(zāi)害風(fēng)險預(yù)估值,從而為有針對性地開展災(zāi)害風(fēng)險管理和風(fēng)險分散提供科學(xué)參考(圖1)。早在19世紀(jì)末期,災(zāi)害觀測技術(shù)的提高以及風(fēng)險制圖技術(shù)的發(fā)展促使基于計算機技術(shù)進行巨災(zāi)損失概率模擬和評估的巨災(zāi)模型(catastrophe models)應(yīng)運而生[7-8]。受益于熱帶氣旋風(fēng)險研究的進步,國際上的巨災(zāi)模型公司(如AIR、RMS等)、保險或再保險公司(如Swiss Re、Munich Re等)以及美國政府(如FEMA開發(fā)了HAZUS-MH颶風(fēng)模型)、澳大利亞政府等紛紛開發(fā)了巨災(zāi)風(fēng)險模型[9]。

    針對熱帶氣旋巨災(zāi)風(fēng)險的評估和模擬問題,一些學(xué)者從致災(zāi)因子危險性[10-12]、建筑物脆弱性[13]等不同的角度進行了綜述。本文在已有研究的基礎(chǔ)上,從面向災(zāi)害管理部門開展熱帶氣旋災(zāi)害風(fēng)險管理的角度出發(fā),分別對國際上熱帶氣旋致災(zāi)因子危險性評估和承災(zāi)體脆弱性評估這兩項核心內(nèi)容的研究進展進行綜述,總結(jié)熱帶氣旋災(zāi)害風(fēng)險定量評估研究中存在的問題及其挑戰(zhàn)。此外,對我國熱帶氣旋(臺風(fēng))災(zāi)害風(fēng)險管理體系中風(fēng)險定量評估的研究現(xiàn)狀進行了分析和討論,針對其不足給出了具體的建議,為我國沿海地區(qū)進一步加強臺風(fēng)災(zāi)害風(fēng)險管理提供了科學(xué)參考。

    圖1 熱帶氣旋巨災(zāi)風(fēng)險模型的結(jié)構(gòu)

    1 熱帶氣旋致災(zāi)因子危險性評估

    相較于其他自然災(zāi)害,熱帶氣旋的致災(zāi)因子具有多樣性和復(fù)雜性。多樣性體現(xiàn)在,既包含大風(fēng)和暴雨等兩大直接致災(zāi)因子,又包含洪水、風(fēng)暴潮(及其引發(fā)的沿海洪水)等次生致災(zāi)因子。復(fù)雜性體現(xiàn)在,直接致災(zāi)因子和次生致災(zāi)因子可能引發(fā)“災(zāi)害鏈”現(xiàn)象,而風(fēng)、雨、潮以及天文大潮等多種致災(zāi)因子“并發(fā)”形成“復(fù)合災(zāi)害”,可能導(dǎo)致生命、財產(chǎn)、社會經(jīng)濟系統(tǒng)的巨大損失,較單一致災(zāi)因子的情況更為復(fù)雜。

    1.1 熱帶氣旋大風(fēng)危險性評估

    熱帶氣旋大風(fēng)危險性評估主要采用數(shù)值模擬方法,根據(jù)歷史上發(fā)生的熱帶氣旋災(zāi)害資料,確定特定地點遭遇不同強度(通常用風(fēng)速來表征)熱帶氣旋災(zāi)害的超越概率(exceedance probabilities)。自基于概率的熱帶氣旋危險性模擬方法于1960年代末被Russell[14]提出以來,其基本框架延續(xù)至今,其大致流程如圖2所示。

    在路徑模型方面,早期的研究主要采用“樣本圓或海岸段”法,以研究點(海岸段)為中心,選取附近一定范圍內(nèi)的歷史熱帶氣旋作為分析樣本,構(gòu)建熱帶氣旋關(guān)鍵參數(shù)的概率分布模型,從中抽樣生成大量模擬熱帶氣旋樣本[15-17]。這種方法的不足之處在于各關(guān)鍵參數(shù)的分布僅適用于特定地點或小范圍區(qū)域,由于局部區(qū)域歷史熱帶氣旋樣本數(shù)量、尤其是極端熱帶氣旋災(zāi)害樣本較少,可能低估熱帶氣旋致災(zāi)因子的危險性。Vickery等[18]提出了“全路徑”模擬法(full-track modelling),以整個海域的歷史熱帶氣旋記錄為樣本,分析各關(guān)鍵參數(shù)的時空分布特征,使用隨機抽樣法模擬大量熱帶氣旋從生成到消亡的整個生命過程中的路徑和強度的變化。這種方法的優(yōu)勢在于其模擬樣本可用于分析任意感興趣點的熱帶氣旋大風(fēng)危險性,但計算較為復(fù)雜,計算量較大,對數(shù)據(jù)的要求較高。對于如何模擬各路徑點的位置,不同學(xué)者提出了不同的方法,包括參數(shù)化方法[18]、半?yún)?shù)化方法[19]、統(tǒng)計與動力相結(jié)合的方法等[20]。

    對于強度的模擬主要包括純統(tǒng)計方法和動力統(tǒng)計法,后者通過建立絕對強度(中心氣壓)與海表溫度、對流層溫度等環(huán)境參數(shù)之間的轉(zhuǎn)換公式來考慮風(fēng)切變和海洋混合等因素對熱帶氣旋強度變化的影響[21,18]。

    風(fēng)場模型描述上層梯度風(fēng)速的徑向分布,一般分為動力數(shù)值模型(dynamical model)和參數(shù)模型(parametric model)兩種類型。由于參數(shù)模型則能夠使用少量關(guān)鍵參數(shù)理想化地表達熱帶氣旋的風(fēng)場結(jié)構(gòu),計算量遠遠小于動力數(shù)值模型,因而在熱帶氣旋風(fēng)險評估模型中得到了廣泛的應(yīng)用。簡單的參數(shù)風(fēng)場模型有Rankine渦旋模型、修正的Rankine渦旋模型[22]、SLOSH風(fēng)場模型[23]等圓對稱風(fēng)場模型。Holland模型[24]通過引入形狀參數(shù)B來描述不同熱帶氣旋氣壓剖面的形狀差異,并由氣壓分布根據(jù)梯度平衡方程推導(dǎo)最大風(fēng)速及梯度風(fēng)速的徑向分布,在熱帶氣旋風(fēng)險研究中得到了廣泛的應(yīng)用[16,25-26]。近期的一些研究探討使用B值與其他臺風(fēng)參數(shù)之間的函數(shù)關(guān)系式來計算B值[27-29]。除了參數(shù)風(fēng)場模型以外,一些學(xué)者通過將復(fù)雜的三維數(shù)值風(fēng)場模型(如文獻[30])參數(shù)化以達到既能夠考慮地表摩擦力、海陸相互作用等因素的影響,又能夠?qū)崿F(xiàn)大量隨機模擬的目的[31-32]。

    圖2 熱帶氣旋大風(fēng)危險性概率模擬過程

    路徑模型、強度模型和風(fēng)場模型結(jié)合其他模型(圖2)生成感興趣研究點(區(qū))的風(fēng)速樣本,最后使用參數(shù)或非參數(shù)概率曲線擬合方法得到大風(fēng)風(fēng)速的超越概率曲線,可以計算任意重現(xiàn)期下的大風(fēng)風(fēng)速。總體而言,局部范圍的路徑模擬方法和全海域的全路徑模擬方法各有其優(yōu)缺點,前者適用于局部地區(qū)的熱帶氣旋危險性模擬,計算量小,后者計算量大,對數(shù)據(jù)的長度和質(zhì)量要求更高,但能夠更全面而精確地考慮整個海岸線上的熱帶氣旋危險性。

    1.2 熱帶氣旋暴雨危險性評估

    熱帶氣旋暴雨形成機制較為復(fù)雜,既包含熱帶氣旋本身環(huán)流形成的降水,又包括熱帶氣旋外圍環(huán)流與中緯度系統(tǒng)相互作用產(chǎn)生的遠距離暴雨,同時受到地形、下墊面、中小尺度系統(tǒng)等因素的影響,因此也是熱帶氣旋系統(tǒng)中最難模擬的變量之一[33-34]。一些學(xué)者利用衛(wèi)星遙感降水產(chǎn)品或是地表降水觀測數(shù)據(jù)探討了熱帶氣旋降水的氣候特征,從中發(fā)現(xiàn)熱帶氣旋降水的年代際變化規(guī)律、長期變化趨勢,以及氣候變化對熱帶氣旋降水的影響等[35-37]。熱帶氣旋暴雨危險性評估則需要建立定量的熱帶氣旋暴雨預(yù)報和預(yù)測模型,盡管這對于復(fù)雜的熱帶氣旋暴雨系統(tǒng)來說非常具有挑戰(zhàn)性。

    HAZUS-MH颶風(fēng)模型使用從專題微波輻射成像儀(SSM/I)中提取的颶風(fēng)降水率數(shù)據(jù)構(gòu)建了颶風(fēng)降水率與離颶風(fēng)中心的距離之間的簡單經(jīng)驗?zāi)P停⒖紤]颶風(fēng)中心氣壓、颶風(fēng)移動速度、臺風(fēng)登陸位置等因素對模型進行修正[38]。Marks等利用NASA熱帶降雨測量計劃(TRMM)衛(wèi)星提供的熱帶氣旋降雨數(shù)據(jù)改進了原有的基于雨量計數(shù)據(jù)的降水氣候與持續(xù)性模型(R-CLIPER),該模型將降水率表示為徑向距離和登陸時間的函數(shù)[39]。Lonfat等在R-CLIPER的基礎(chǔ)上考慮垂直切變和地形對降水場不對稱性的影響提出了參數(shù)化颶風(fēng)降水模型(PHRaM),改善了颶風(fēng)降水的模擬結(jié)果[40,34]。

    相較于簡單的經(jīng)驗?zāi)P?,?fù)雜的三維數(shù)值預(yù)報模型能夠更精確地模擬熱帶氣旋降水場在時間和空間上的變化,但這些模型在啟動方案、分辨率、物理參數(shù)以及最終的預(yù)報結(jié)果上均存在較大差異[41,34]。對比研究發(fā)現(xiàn),全球預(yù)報系統(tǒng)(GFS)、GFDL颶風(fēng)預(yù)報模型及北美中尺度(NAM)模型這三種主要的數(shù)值模型在颶風(fēng)降水預(yù)報上的效果都相當(dāng)于或好于R-CLIPER模型,但各數(shù)值模型又有其各自的局限性[41]。

    1.3 熱帶氣旋風(fēng)暴潮危險性評估

    對熱帶氣旋風(fēng)暴潮危險性的評估主要有兩種方式:基于典型重現(xiàn)期的風(fēng)暴潮危險性評估和基于熱帶氣旋風(fēng)場模擬和風(fēng)暴潮數(shù)值模型的危險性評估。前者是對風(fēng)暴潮的歷史實測資料進行概率分析,獲得不同重現(xiàn)期下的風(fēng)暴潮增水,如百年一遇風(fēng)暴潮作為風(fēng)暴潮危險性的標(biāo)準(zhǔn),該方法一方面受限于歷史資料長度和精讀的限制,另一方面無法考慮單個風(fēng)暴潮事件的持續(xù)時間、影響范圍等因素。而后者則是基于歷史或假想的熱帶氣旋風(fēng)場使用風(fēng)暴潮數(shù)值模型模擬風(fēng)暴潮增水,以預(yù)測風(fēng)暴潮個例的實時變化、估計風(fēng)暴潮的重現(xiàn)期或最大可能風(fēng)暴潮(Probable Maximum Storm Surge, PMSS)[42-43]。

    在風(fēng)暴潮的實時預(yù)報中,通常使用高精度的數(shù)值天氣預(yù)報模式對熱帶氣旋風(fēng)場及其引發(fā)的風(fēng)暴潮增水進行估計,在長期的風(fēng)暴潮風(fēng)險評估中則主要使用簡單高效的參數(shù)化風(fēng)場模型進行風(fēng)場模擬[42]。美國國家氣象局(NWS)使用SLOSH模式模擬颶風(fēng)風(fēng)場驅(qū)動下的風(fēng)暴潮增水,其生成的產(chǎn)品不僅包括基于歷史颶風(fēng)風(fēng)場的概率性風(fēng)暴潮增水、還包括基于假想颶風(fēng)風(fēng)場的風(fēng)暴潮增水最大包絡(luò)圖和最大風(fēng)暴潮增水最大包絡(luò)圖等。我國也建立了覆蓋整個中國海的高分辨率風(fēng)暴潮集合數(shù)值預(yù)報系統(tǒng),用于沿海風(fēng)暴潮危險性評估的業(yè)務(wù)中[44]。

    2 承災(zāi)體脆弱性評估

    承災(zāi)體脆弱性評估是連接致災(zāi)因子危險性和承災(zāi)體的損失程度之間的紐帶,也是風(fēng)險評估中最為復(fù)雜的一部分。反映損失率和致災(zāi)因子危險性之間的確定性函數(shù)被稱為脆弱性函數(shù)(vulnerability function)或災(zāi)損函數(shù)(damage function),一般表現(xiàn)為曲線的形式,被稱為脆弱性曲線(vulnerability curve)或易損性曲線(fragility curve)。國內(nèi)外熱帶氣旋脆弱性評估研究中的承災(zāi)體對象主要為建筑物、基礎(chǔ)設(shè)施等造成直接經(jīng)濟損失的有形資產(chǎn)。

    2.1 建筑物脆弱性評估

    目前針對建筑結(jié)構(gòu)的熱帶氣旋脆弱性研究主要以大風(fēng)對建筑的破壞為研究對象(熱帶氣旋引發(fā)的洪水災(zāi)害單獨模擬),采用歷史損失數(shù)據(jù)法、專家經(jīng)驗法、歷史損失數(shù)據(jù)和專家經(jīng)驗結(jié)合法、建筑構(gòu)件可靠性法以及工程實驗和數(shù)值模擬法等方法建立建筑的脆弱性模型[13]。其中前三種方法依賴于歷史數(shù)據(jù)和主觀經(jīng)驗,能夠在一定程度上描述一組特征類似的建筑物損失程度與熱帶氣旋強度(通常是風(fēng)速)之間的關(guān)系[45-48];后兩種方法則對由建筑結(jié)構(gòu)特征、熱帶氣旋風(fēng)壓與建筑各構(gòu)件之間相互作用的過程進行更為細致的模擬,脆弱性模型更為精確,但所需的計算和模擬資源更多[49-53]。

    建筑物的熱帶氣旋脆弱性研究集中在美國和澳大利亞這兩個國家,此外日本和其他國家有少量的研究[13]。美國早在1970年代就針對居民住宅和財產(chǎn)推出了包括颶風(fēng)、洪水、地震等自然災(zāi)害在內(nèi)的房屋和財產(chǎn)保險。HAZUS-MH颶風(fēng)模型以建筑物(包括居民住宅、活動房屋、商業(yè)建筑等)承災(zāi)體引起的直接經(jīng)濟損失為模擬對象,針對各類建筑的結(jié)構(gòu)特點從建筑外部結(jié)構(gòu)破壞、建筑內(nèi)部破壞、建筑內(nèi)部財產(chǎn)損失及使用損失等各方面分別構(gòu)建損失函數(shù)[38]。澳大利亞的建筑脆弱性曲線最早是根據(jù)澳大利亞歷史上最嚴(yán)重的熱帶氣旋災(zāi)害之一Tracy的災(zāi)后損失調(diào)查數(shù)據(jù)建立的[54],Walker的模型在當(dāng)時被認為是最好的模型,被其他學(xué)者廣泛使用[55]。近些年來,一些學(xué)者則從建筑構(gòu)件的可靠性出發(fā)構(gòu)建了澳大利亞典型建筑的脆弱性曲線[56-57]。

    此外,一些研究認為社會脆弱性因素對建筑的損失也具有重要影響,例如Burton[58]將社會脆弱性指標(biāo)與颶風(fēng)大風(fēng)和風(fēng)暴潮一起與卡特里娜颶風(fēng)損失數(shù)據(jù)做多元回歸分析,發(fā)現(xiàn)社會脆弱性因子與颶風(fēng)損失之間存在顯著的關(guān)系。

    2.2 其他承災(zāi)體的脆弱性評估

    隨著預(yù)報預(yù)警、應(yīng)急疏散機制的完善,熱帶氣旋災(zāi)害對人口的影響越來越小,死亡和失蹤人口較早期有明顯的下降。熱帶氣旋災(zāi)害中的人員傷亡情況大多與建筑倒塌、洪澇災(zāi)害或飛濺物等次生災(zāi)害有關(guān),與熱帶氣旋災(zāi)害強度之間不存在明顯的相關(guān)關(guān)系,因此對人口的物理脆弱性研究較少,主要集中在社會脆弱性層面,例如Wang等[59]研究了老年人口對颶風(fēng)引起的風(fēng)暴潮和洪澇災(zāi)害的社會脆弱性在程度和結(jié)構(gòu)上的空間差異。盡管人口脆弱性研究無法給出定量的結(jié)果,探究造成熱帶氣旋災(zāi)害中人員傷亡差異的影響因素對減輕人口災(zāi)害風(fēng)險仍然有著重要的意義。為此,近年來學(xué)術(shù)界呼吁加強社會經(jīng)濟要素、人口主觀風(fēng)險認知、災(zāi)害應(yīng)急疏散等災(zāi)害社會科學(xué)的研究[60]。

    3 我國的熱帶氣旋風(fēng)險評估現(xiàn)狀

    從我國目前的臺風(fēng)災(zāi)害研究來看,研究內(nèi)容主要集中在區(qū)域臺風(fēng)災(zāi)害風(fēng)險的區(qū)劃研究[61-64]、臺風(fēng)綜合災(zāi)情指數(shù)的預(yù)測研究[65-67]和臺風(fēng)致災(zāi)因子危險性的預(yù)測和評估研究[68-71]等三個方面,在系統(tǒng)性的、針對具體承災(zāi)體的臺風(fēng)災(zāi)害風(fēng)險評估研究上的工作非常有限。

    在臺風(fēng)概率危險性模擬方面,我國學(xué)者嘗試?yán)梦鞅碧窖笈_風(fēng)歷史記錄對國際上已有的風(fēng)場模型和陸地衰減模型等進行參數(shù)校正,進而針對我國東南沿海地區(qū)的臺風(fēng)危險性開展一系列研究[68-71]。然而,在路徑隨機模擬、臺風(fēng)風(fēng)場構(gòu)建、臺風(fēng)概率危險性評估等研究工作數(shù)量偏少,減災(zāi)工作大多依賴于氣象部門臺風(fēng)模式的短期預(yù)報結(jié)果作為風(fēng)、雨等致災(zāi)因子危險性的輸入。

    承災(zāi)體脆弱性方面,我國學(xué)者在土木工程領(lǐng)域相繼開展了一些有關(guān)建筑大風(fēng)災(zāi)害脆弱性的研究[72-73]、建筑單體和區(qū)域建筑群的脆弱性等級研究[74-75],盡管尚未系統(tǒng)地應(yīng)用于臺風(fēng)災(zāi)害的風(fēng)險評估中,這些研究為我國的臺風(fēng)巨災(zāi)風(fēng)險的脆弱性評估研究提供了參考。但總體來看,我國的熱帶氣旋災(zāi)害風(fēng)險研究還有待加強和深入,針對特定承災(zāi)體的臺風(fēng)災(zāi)害風(fēng)險定量研究幾乎處于空白狀態(tài),尤其是特定承災(zāi)體的臺風(fēng)災(zāi)害脆弱性研究工作非常稀少。

    4 總結(jié)與討論

    本文從災(zāi)害風(fēng)險定量評估的角度分別評述了熱帶氣旋致災(zāi)因子危險性模型和脆弱性模型的研究現(xiàn)狀與發(fā)展方向,為未來進一步開展沿海地區(qū)熱帶氣旋災(zāi)害風(fēng)險的定量評估研究提供了參考。熱帶氣旋致災(zāi)因子包括大風(fēng)、暴雨、風(fēng)暴潮等三種主要類型,其中熱帶氣旋大風(fēng)危險性評估方法相對成熟,暴雨和風(fēng)暴潮通常與洪水災(zāi)害風(fēng)險聯(lián)系起來開展風(fēng)險評估。主要難點在于如何更準(zhǔn)確地模擬熱帶氣旋的路徑、強度、風(fēng)場和降水等過程或特征,而這有賴于熱帶氣旋結(jié)構(gòu)和風(fēng)雨過程的精細化觀測的進一步加強,以及對熱帶氣旋形成和發(fā)展的物理機制的深入理解。當(dāng)前熱帶氣旋風(fēng)險評估中以建筑物(直接經(jīng)濟損失)為主要承災(zāi)體(風(fēng)險指標(biāo)),對其他承災(zāi)體的脆弱性和風(fēng)險的定量評估研究還有待加強。

    盡管針對熱帶氣旋災(zāi)害風(fēng)險模型的研究已經(jīng)取得了顯著的進展,還有眾多相關(guān)科學(xué)問題需要進一步地研究,例如氣候變化背景下熱帶氣旋危險性的變化、承災(zāi)體脆弱性隨時間的變化、熱帶氣旋引發(fā)的風(fēng)暴潮及暴雨洪澇等多災(zāi)種的綜合風(fēng)險評估等。在眾多相關(guān)問題中,全球氣候變化對熱帶氣旋災(zāi)害風(fēng)險的影響正受到廣泛的關(guān)注[76]。研究表明,熱帶氣旋的長期活動特征(如強度、頻次和持續(xù)時間等)并沒有顯著的增加,沒有證據(jù)表明人類活動(氣候變化)對熱帶氣旋活動變化產(chǎn)生了顯著影響,但熱帶氣旋的最大風(fēng)速平均值有可能(66-100%的概率)上升,極端熱帶氣旋的熱帶氣旋引起的降水有可能增加;此外,人口社會暴露度的增加將導(dǎo)致更高的熱帶氣旋災(zāi)害損失的產(chǎn)生[5]。一些學(xué)者通過對歷史熱帶氣旋損失進行歸一化研究發(fā)現(xiàn),社會經(jīng)濟暴露度的增加是導(dǎo)致熱帶氣旋損失增加的主要原因[77-80]。最近,Estrada等使用基于回歸的歸一化方法代替?zhèn)鹘y(tǒng)的比例損失歸一化方法分析美國歷史颶風(fēng)損失,發(fā)現(xiàn)歸一化后的損失呈現(xiàn)與氣候變化信號相一致的增加趨勢[81],但該方法的可靠性仍然存在爭議[82]。

    我國正處于社會經(jīng)濟快速發(fā)展的階段,在全球變暖背景下,臺風(fēng)強度和強臺風(fēng)的頻次可能增加,沿海地區(qū)臺風(fēng)災(zāi)害風(fēng)險將日益增加,人口、社會、經(jīng)濟面臨嚴(yán)峻挑戰(zhàn)。臺風(fēng)災(zāi)害風(fēng)險評估對于政府防災(zāi)減災(zāi)、巨災(zāi)保險的開展以及建筑工程的規(guī)劃和設(shè)計都具有重要的意義,因此我國急需探索開展熱帶氣旋巨災(zāi)風(fēng)險評估的綜合研究,尤其是針對不同區(qū)域建筑物、基礎(chǔ)設(shè)施等特定承災(zāi)體的脆弱性研究,使學(xué)術(shù)研究、巨災(zāi)保險和政府災(zāi)害管理三者之間形成合力,為沿海地區(qū)減輕臺風(fēng)災(zāi)害風(fēng)險提供科學(xué)工具和有效手段。同時,加強對臺風(fēng)氣候變化的觀測研究,深入開展氣候變化對臺風(fēng)活動影響預(yù)測和模擬研究,以便為制定長期的臺風(fēng)風(fēng)險管理決策提供科學(xué)依據(jù)。

    [1] Aon Benfield. 2015 Annual Global Climate and Catastrophe Report[R]. Chicago, Illinois, 2016.

    [2] Pielke Jr R A and Landsea C W. Normalized hurricane damages in the United States: 1925-95[J]. Weather and Forecasting, 1998, 13(3): 621-631.

    [3] Raghavan S and Rajesh S. Trends in tropical cyclone impact: A study in Andhra Pradesh, India[J]. Bulletin of the American Meteorological Society, 2003, 84(5): 635-644.

    [4] WMO. Seventh International Workshop on Tropical Cyclones (IWTC-VII)[EB/OL]. (2009-12-31) [2017-02-10]. http://www.wmo.int/pages/prog/arep/wwrp/tmr/IWTC-VII.html.

    [5] IPCC. Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation. A Special Report of the Intergovernmental Panel on Climate Change[R]. Cambridge, UK, and New York, NY, USA, Cambridge University Press, 2012: 582.

    [6] UNISDR. Living with risk: A Global Review of Disaster Reduction Initiatives[R]. Geneva, Switzerland, United Nations, 2004.

    [7] Grossi P and TeHennepe C. Catastrophe Modelling Fundamentals[M]//The Review-A Guide to Catastrophe Modelling, 2008: 6-11.

    [8] Walker G R. Modelling the vulnerability of buildings to wind-a review[J]. Canadian Journal of Civil Engineering, 2011, 38(9): 1031-1039.

    [9] 陳克平. 災(zāi)難模型化及其國外主要開發(fā)商[J]. 自然災(zāi)害學(xué)報, 2004, 13(2): 1-8.

    [10] Vickery P J, Masters F J, Powell M D, et al. Hurricane hazard modeling: The past, present, and future[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2009, 97(7): 392-405.

    [11] 方偉華,石先武. 面向災(zāi)害風(fēng)險評估的熱帶氣旋路徑及強度隨機模擬綜述[J]. 地球科學(xué)進展, 2012, 27(8): 866-875.

    [12] 方偉華,林偉. 面向災(zāi)害風(fēng)險評估的臺風(fēng)風(fēng)場模型研究綜述[J]. 地理科學(xué)進展, 2013, 32(6): 852-867.

    [13] Pita G, Pinelli J P, Gurley K, et al. State of the art of hurricane vulnerability estimation methods: A review[J]. Natural Hazards Review, 2014, 16(2): 04014022.

    [14] Russell L R. Probability Distributions for Texas Gulf Coast Hurricane Effects of Engineering Interest[D]. Palo Alto: Stanford University,1968.

    [15] Batts M, Simiu E, and Russell L. Hurricane wind speeds in the United States[J]. Journal of the Structural Division, 1980, 106(10): 2001-2016.

    [16] Georgiou P, Davenport A G, and Vickery B. Design wind speeds in regions dominated by tropical cyclones[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1983, 13(1): 139-152.

    [17] Neumann C J. The National Hurricane Center Risk Analysis Program (HURISK)[C]// NOAA Technical Memorandum NWS NHC 38. Washington, DC, National Oceanic and Atmospheric Administration (NOAA).1991.

    [18] Vickery P, Skerlj P and Twisdale L. Simulation of hurricane risk in the U.S. using empirical track model[J]. Journal of Structural Engineering, 2000, 126(10): 1222-1237.

    [19] Hall T and Jewson S. Statistical modelling of tropical cyclone tracks: a semi-parametric model for the mean trajectory[J]. Physics, 0503231, 2005.

    [20] Emanuel K, Ravela S, Vivant E, et al. A statistical deterministic approach to hurricane risk assessment[J]. Bulletin of the American Meteorological Society, 2006, 87(3): 299-314.

    [21] Darling R W R. Estimating probabilities of hurricane wind speeds using a large-scale empirical model[J]. Journal of Climate, 1991, 4(10): 1035-1046.

    [22] Phadke A C, Martino C D, Cheung K F, et al. Modeling of tropical cyclone winds and waves for emergency management[J]. Ocean Engineering, 2003, 30(4): 553-578.

    [23] Jelesnianski C P, Chen J, and Shaffer W A. SLOSH: Sea, Lake, and Overland Surges from Hurricanes[C]// NOAA Technical Report NWS 48. Silver Spring, MD, National Weather Service.1992.

    [24] Holland G. An analytic model of the wind and pressure profiles in hurricanes[J]. Monthly weather review, 1980, 108(8): 1212-1218.

    [25] Harper B A and Holland G J. An updated parametric model of tropical cyclone[C]// Proceedings of 23rd Conference on Hurricanes and Tropical Meteorology. Dallas, Texas, American Meteorological Society.1999.

    [26] Lee K H and Rosowsky D V. Synthetic hurricane wind speed records: development of a database for hazard analyses and risk studies[J]. Natural Hazards Review, 2007, 8(2): 23-34.

    [27] Powell M, Soukup G, Cocke S, et al. State of Florida hurricane loss projection model: Atmospheric science component[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2005, 93(8): 651-674.

    [28] Holland G. A revised hurricane pressure-wind model[J]. Monthly Weather Review, 2008, 136(9): 3432-3445.

    [29] Vickery P J and Wadhera D. Statistical models of Holland pressure profile parameter and radius to maximum winds of hurricanes from flight-level pressure and H*Wind data[J]. Journal of Applied Meteorology and Climatology, 2008, 47(10): 2497-2517.

    [30] Shapiro L J. The asymmetric boundary layer flow under a translating hurricane[J]. Journal of the Atmospheric Sciences, 1983, 40(8): 1984-1998.

    [31] Georgiou P. Design Windspeeds in Tropical Cyclone-Prone Regions[D]. London, Ontario, Canada: University of Western Ontario,1985.

    [32] Vickery P, Skerlj P, Steckley A, et al. Hurricane wind field model for use in hurricane simulations[J]. Journal of Structural Engineering, 2000, 126(10): 1203-1221.

    [33] 程正泉, 陳聯(lián)壽, 徐祥德, 等. 近10年中國臺風(fēng)暴雨研究進展[J]. 氣象, 2005, 31(12): 3-9.

    [34] Rogers R, Marks F, and Marchok T. Tropical cyclone rainfall[J]. Encyclopedia of hydrological sciences, 2009.

    [35] Ren F, Wu G, Dong W, et al. Changes in tropical cyclone precipitation over China[J]. Geophysical Research Letters, 2006, 332(20): 131-145.

    [36] David K B and Davis R E. Climatology of Tropical Cyclone Rainfall in the Southeastern United States[J]. Physical Geography, 2007, 28(2): 126-147.

    [37] Shepherd J M, Grundstein A, and Mote T L. Quantifying the contribution of tropical cyclones to extreme rainfall along the coastal southeastern United States[J]. Geophysical Research Letters, 2007, 34(23): L23810.

    [38] FEMA. Hazus-MH 2.1 Hurricane Model Technical Manual[R]. Washington, DC, Department of Homeland Security, Federal Emergency Management Agency, Mitigation Division.2016.

    [39] Marks F, Kappler G and Demaria M. Development of a tropical cyclone rainfall climatology and persistence (R-CLIPER) model[C]// Proceedings of Conference on Hurricanes & Tropical Meteorology, 2002.

    [40] Lonfat M, Rogers R, Marchok T, et al. A parametric model for predicting hurricane rainfall[J]. Monthly Weather Review, 2007, 135(9): 3086-3097.

    [41] Marchok T, Rogers R and Tuleya R. Validation schemes for tropical cyclone quantitative precipitation forecasts: evaluation of operational models for U.S. landfalling cases[J]. Weather & Forecasting, 2007, 22(4): 726-746.

    [42] Lin N and Chavas D. On hurricane parametric wind and applications in storm surge modeling[J]. Journal of Geophysical Research, 2012, 117(D9): D09120.

    [43] 石先武, 譚駿, 國志興, 等. 風(fēng)暴潮災(zāi)害風(fēng)險評估研究綜述[J]. 地球科學(xué)進展, 2013, 28(8): 866-874.

    [44] 張月霞. 基于ADCIRC模式的臺風(fēng)風(fēng)暴潮災(zāi)害風(fēng)險評估方法研究與應(yīng)用——以浙江臺州為例[D]. 青島: 中國海洋大學(xué), 2014.

    [45] Hendrick R L and Friedman D G. Potential impacts of storm modification on the insurance industry. Human Dimensions of Weather Modification[C]//W. Sewell, University of Chicago Press, Chicago, IL.1966: 227-248.

    [46] Friedman D G. Natural hazard risk assessment for an insurance program[C]//Geneva Papers on Risk and Insurance, 1984: 57-128.

    [47] Jain V K, Guin J and He H. Statistical Analysis of 2004 and 2005 Hurricane Claims Data[C]//Proceedings of the 11th Americas Conference on Wind Engineering. San Juan, Puerto Rico.2009.

    [48] Kafali C and Jain V. A Methodology for Estimating Typhoon Losses to Building Environment in Japan[C]//Proceedings of the 11th Americas Conference on Wind Engineering. San Juan, Puerto Rico.2009.

    [49] Stubbs N and Perry D. A damage simulation model for buildings and contents in a hurricane environment[C]// Building an International Community of Structural Engineers. ASCE, 1996.

    [50] Applied Research Associates (ARA). Development of loss relativities for wind resistive features of residential structures[C]//Version 2.2. Raleigh, NC, Florida Department of Community Affairs, 2002.

    [51] Li Y and Ellingwood B R. Hurricane damage to residential construction in the US: Importance of uncertainty modeling in risk assessment[J]. Engineering Structures, 2006, 28(7): 1009-1018.

    [52] Vickery P J, Skerlj P F, Lin J, et al. HAZUS-MH hurricane model methodology. II: Damage and loss estimation[J]. Natural Hazards Review, 2006, 7(2): 94-103.

    [53] Hamid S, Kibria B G, Gulati S, et al. Predicting losses of residential structures in the state of Florida by the public hurricane loss evaluation model[J]. Statistical Methodology, 2010, 7(5): 552-573.

    [54] Walker G. CSIRO Division of Building, Construction and Engineering, Dr[R]. BA Harper, 1994.

    [55] Stewart M. Risk and economic analysis of residential housing climate adaptation strategies for wind hazards in southeast Australia[C]//CAEx Report 1/2013, Climate Adaptation Engineering for Extreme Events Cluster, CSIRO Climate Adaptation Flagship, Canberra.2013.

    [56] Henderson D J and Ginger J D. Vulnerability model of an Australian high-set house subjected to cyclonic wind loading[J]. Wind & Structures An International Journal, 2007, 10(3): 269-285.

    [57] Wehner M, Ginger J, Holmes J, et al. Development of methods for assessing the vulnerability of Australian residential building stock to severe wind[C]//IOP Conference Series: Earth and Environmental Science, 2010, 11(2010) 012017.

    [58] Burton C G. Social Vulnerability and Hurricane Impact Modeling[J]. Natural Hazards Review, 2010, 11(2): 58-68.

    [59] Wang C and Yarnal B. The vulnerability of the elderly to hurricane hazards in Sarasota, Florida[J]. Natural Hazards, 2012, 2012(63): 349-373.

    [60] Jung K, Shavitt S, Viswanathan M, et al. Female hurricanes are deadlier than male hurricanes[J]. Proceedings of the National Academy of Sciences of the United States of America, 2014, 111(24): 8782-8787.

    [61] 史培軍. 三論災(zāi)害研究的理論與實踐[J]. 自然災(zāi)害學(xué)報, 2002, 11(3): 1-9.

    [62] 陳香. 福建省臺風(fēng)災(zāi)害風(fēng)險評估與區(qū)劃[J]. 生態(tài)學(xué)雜志, 2007, 26(6): 961-966.

    [63] 牛海燕, 劉敏, 陸敏, 等. 中國沿海地區(qū)臺風(fēng)災(zāi)害損失評估研究[J]. 災(zāi)害學(xué), 2011, 26(3): 61-64.

    [64] 殷潔, 戴爾阜, 吳紹洪. 中國臺風(fēng)災(zāi)害綜合風(fēng)險評估與區(qū)劃[J]. 地理科學(xué), 2013, 33(11): 1370-1376.

    [65] 樊琦, 梁必騏. 熱帶氣旋災(zāi)情的預(yù)測及評估[J]. 地理學(xué)報, 2000, 55(S1): 52-56.

    [66] 魏章進, 隋廣軍, 唐丹玲. 臺風(fēng)災(zāi)情評估及方法綜述[J]. 災(zāi)害學(xué), 2012, 27(4): 107-113.

    [67] 葉小嶺, 施珮, 匡亮. 基于粒子群優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的臺風(fēng)災(zāi)損預(yù)測模型研究[J]. 災(zāi)害學(xué), 2013, 28(4): 11-15.

    [68] 歐進萍, 段忠東, 常亮. 中國東南沿海重點城市臺風(fēng)危險性分析[J]. 自然災(zāi)害學(xué)報, 2002, 11(4): 9-17.

    [69] 李瑞龍. 基于改進的臺風(fēng)關(guān)鍵參數(shù)的臺風(fēng)極值風(fēng)速預(yù)測[D]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2007.

    [70] 程鴻. 基于歷史觀測數(shù)據(jù)的中國臺風(fēng)致災(zāi)因子過程重建與蒙特拉洛模擬[D]. 北京: 北京師范大學(xué),2009.

    [71] 林偉, 方偉華. 西北太平洋臺風(fēng)風(fēng)場模型中 Holland B 系數(shù)區(qū)域特征研究[J]. 熱帶地理, 2013, 33(2): 124-132.

    [72] 王韶文. 廣東沿海典型低矮建筑調(diào)研與風(fēng)荷載特性的試驗研究[D]. 汕頭: 汕頭大學(xué), 2009.

    [73] 何春保, 嚴(yán)英杰, 譚文怡, 等. 華南地區(qū)農(nóng)村建筑抗風(fēng)減災(zāi)情況調(diào)查與對策分析[J]. 南方建筑, 2011(5): 53-55.

    [74] 辜智慧, 王娟, 葛怡, 等. 基于PTVA修正模型的建筑單體臺風(fēng)脆弱性評估方法[J]. 中國安全科學(xué)學(xué)報, 2015, 25(11): 99-105.

    [75] 辜智慧, 王娟, 葛怡, 等. 區(qū)域建筑群臺風(fēng)災(zāi)害綜合脆弱性評估方法研究:以廣東省為例[J]. 中國安全科學(xué)學(xué)報, 2016, 26(2): 108-113.

    [76] 雷小途, 徐明, 任福民. 全球變暖對臺風(fēng)活動影響的研究進展[J]. 氣象學(xué)報, 2009, 67(5): 679-688.

    [77] Miller S, Muir-Wood R, and Boissonnade A. An exploration of trends in normalized weather-related catastrophe losses [M]//H. F. Diaz and R. J. Murnane (Editor). Climate Extremes and Society. New York: Cambridge University Press, 2008: 225-247.

    [78] Pielke Jr R A, Gratz J, Landsea C W, et al. Normalized hurricane damage in the United States: 1900-2005[J]. Natural Hazards Review, 2008, 9(1): 29-42.

    [79] Schmidt S, Kemfert C, and H?ppe P. Tropical cyclone losses in the USA and the impact of climate change—A trend analysis based on data from a new approach to adjusting storm losses[J]. Environmental Impact Assessment Review, 2009, 29(6): 359-369.

    [80] 張嬌艷, 吳立廣, 張強. 全球變暖背景下我國熱帶氣旋災(zāi)害趨勢分析[J]. 熱帶氣象學(xué)報, 2011, 27(4): 442-454.

    [81] Estrada F, Botzen W J W and Tol R S J. Economic losses from US hurricanes consistent with an influence from climate change[J]. Nature Geosci, 2015, 8(11): 880-884.

    [82] Hallegatte S. Climate change: Unattributed hurricane damage[J]. Nature Geosci, 2015, 8(11): 819-820.

    Abstract:Global tropical cyclone disaster risks keep rising with the accelerated growth and aggregation of populations and wealth in the coastal areas. In recent decades, refined quantitative assessment of coastal tropical cyclone disaster risk has been regarded as a prerequisite to achieve prospective disaster risk management; therefore, developing the quantitative risk models, which incorporates multi-disciplinary theories and technologies including meteorology, geography, engineering, stochastic mathematics, insurance and so on, is the future direction of tropical cyclone risk management. From the perspective of integrated risk management for tropical cyclone disasters, this paper reviews the status quo of the tropical cyclone risk assessment, including the current development, problems of the two essential composing elements of risk, as the hazards and vulnerability respectively. The paper finally gives a brief overview of research on typhoon risk modelling in China, and provides responding suggestions for further improving quantitative typhoon risk research in coastal areas in China.

    Key words:tropical cyclone; risk assessment; disaster; hazard; vulnerability

    Review on Tropical Cyclone Risk Assessment

    CHEN Wenfang1, DUAN Yihong1, LU Yi2, FANG Jiayi3, 4, SHI Xianwu5and REN Fumin1

    (1.StateKeyLaboratoryofSevereWeather,ChineseAcademyofMeteorologicalScience,Beijing100081,China; 2.ShanghaiTyphoonInstituteofChinaMeteorologicalAdministration,Shanghai200030,China; 3.AcademyofDisasterReductionandEmergencyManagement,MinistryofCivilAffairs&MinistryofEducation,BeijingNormalUniversity,Beijing100875,China; 4.FacultyofEnvironmentandEngineering,UniversityofSouthampton,Southampton,SO17 1BJ; 5.NationalMarineHazardMitigationService,Beijing100194,China)

    陳文方,端義宏,陸逸,等. 熱帶氣旋災(zāi)害風(fēng)險評估現(xiàn)狀綜述[J]. 災(zāi)害學(xué),2017,32(4):146-152. [CHEN Wenfang, DUAN Yihong, LU Yi,et al. Review on Tropical Cyclone Risk Assessment[J]. Journal of Catastrophology,2017,32(4):146-152.

    10.3969/j.issn.1000-811X.2017.04.025.]

    X43;P42

    A

    1000-811X(2017)04-0146-07

    2017-04-14

    2017-06-01

    國家重點基礎(chǔ)研究發(fā)展計劃項目(2015CB452806)

    陳文方(1987-),女,江西南昌人,博士后,主要從事熱帶氣旋災(zāi)害風(fēng)險評估與管理研究. E-mail: wenfang@camscma.cn

    端義宏(1963-),男,江蘇南京人,研究員,主要從事熱帶氣旋機理和海氣相互作用研究. E-mail: duanyh@cma.gov.cn

    10.3969/j.issn.1000-811X.2017.04.025

    猜你喜歡
    風(fēng)暴潮氣旋脆弱性
    2003年3月北大西洋上兩個爆發(fā)性氣旋的“吞并”過程及發(fā)展機制分析?
    2012年“蘇拉”和“達維”雙臺風(fēng)影響的近海風(fēng)暴潮過程
    海洋通報(2021年2期)2021-07-22 07:55:24
    氨還原劑噴槍中保護氣旋流噴射作用的研究
    能源工程(2021年1期)2021-04-13 02:05:56
    防范未來風(fēng)暴潮災(zāi)害的綠色海堤藍圖
    科學(xué)(2020年4期)2020-11-26 08:27:00
    基于多變量LSTM神經(jīng)網(wǎng)絡(luò)模型的風(fēng)暴潮臨近預(yù)報
    海洋通報(2020年6期)2020-03-19 02:10:18
    煤礦電網(wǎng)脆弱性評估
    電子制作(2017年10期)2017-04-18 07:23:09
    北太平洋上一個爆發(fā)性氣旋族的結(jié)構(gòu)分析?
    殺毒軟件中指令虛擬機的脆弱性分析
    基于攻擊圖的工控系統(tǒng)脆弱性量化方法
    2014年3月大西洋上一個爆發(fā)性氣旋的研究
    少妇猛男粗大的猛烈进出视频| 久久综合国产亚洲精品| 亚洲欧美一区二区三区黑人 | 免费久久久久久久精品成人欧美视频 | 日韩一区二区视频免费看| a 毛片基地| 欧美国产精品一级二级三级 | 人妻少妇偷人精品九色| 亚洲国产精品一区三区| 少妇丰满av| 制服丝袜香蕉在线| 蜜桃在线观看..| 麻豆乱淫一区二区| 中文字幕av电影在线播放| 色吧在线观看| 男女免费视频国产| 日韩电影二区| 日韩一区二区三区影片| 卡戴珊不雅视频在线播放| 成人国产av品久久久| 在线天堂最新版资源| 成年女人在线观看亚洲视频| 天堂俺去俺来也www色官网| 国产一级毛片在线| 韩国av在线不卡| 久久久久久人妻| 精品久久久噜噜| 国产在线男女| 欧美变态另类bdsm刘玥| av有码第一页| 欧美 日韩 精品 国产| 国产日韩欧美在线精品| 又黄又爽又刺激的免费视频.| 成年人午夜在线观看视频| 免费大片18禁| 国产精品国产三级专区第一集| 深夜a级毛片| 亚洲电影在线观看av| 少妇人妻精品综合一区二区| 18禁在线无遮挡免费观看视频| 中国美白少妇内射xxxbb| 男人爽女人下面视频在线观看| 久久综合国产亚洲精品| 日韩中文字幕视频在线看片| 少妇人妻一区二区三区视频| 日本黄大片高清| www.av在线官网国产| 亚洲av二区三区四区| 国产毛片在线视频| 国产色婷婷99| 免费黄网站久久成人精品| av黄色大香蕉| av专区在线播放| 欧美人与善性xxx| 久久狼人影院| 蜜桃久久精品国产亚洲av| 人妻少妇偷人精品九色| 免费观看的影片在线观看| 少妇裸体淫交视频免费看高清| 久久人人爽人人爽人人片va| 国产免费一区二区三区四区乱码| 亚洲国产成人一精品久久久| 极品少妇高潮喷水抽搐| h视频一区二区三区| 久久久久久久久久人人人人人人| 午夜福利在线观看免费完整高清在| 成人黄色视频免费在线看| 亚洲美女视频黄频| 97在线视频观看| 久久6这里有精品| 久久久久精品性色| 久久精品久久久久久久性| 久久精品国产a三级三级三级| 日韩一区二区三区影片| 亚洲精品乱码久久久v下载方式| 又爽又黄a免费视频| 如何舔出高潮| 久久久久网色| 国产精品一区二区在线不卡| 七月丁香在线播放| 91aial.com中文字幕在线观看| 亚洲欧美精品自产自拍| 国产探花极品一区二区| 男女无遮挡免费网站观看| 久久久久久久国产电影| 狂野欧美白嫩少妇大欣赏| 国产免费视频播放在线视频| 色哟哟·www| 亚洲综合色惰| 国产午夜精品久久久久久一区二区三区| 全区人妻精品视频| 在线精品无人区一区二区三| 免费人成在线观看视频色| 高清毛片免费看| 性色av一级| 国产淫片久久久久久久久| 大片免费播放器 马上看| 美女xxoo啪啪120秒动态图| 女人精品久久久久毛片| 精品久久久久久久久亚洲| freevideosex欧美| 寂寞人妻少妇视频99o| 黄色日韩在线| 亚洲在久久综合| 亚洲人成网站在线观看播放| 又黄又爽又刺激的免费视频.| 亚洲国产成人一精品久久久| 少妇熟女欧美另类| 老司机亚洲免费影院| 国产成人免费无遮挡视频| 成人漫画全彩无遮挡| 精品久久久久久久久av| 少妇人妻久久综合中文| 久久国产乱子免费精品| 22中文网久久字幕| 成年美女黄网站色视频大全免费 | 赤兔流量卡办理| 国产色爽女视频免费观看| 国产精品蜜桃在线观看| 免费在线观看成人毛片| 大香蕉97超碰在线| 亚洲成人手机| 色哟哟·www| 日韩三级伦理在线观看| 国产极品粉嫩免费观看在线 | 赤兔流量卡办理| 亚洲三级黄色毛片| videos熟女内射| av线在线观看网站| 美女视频免费永久观看网站| 韩国av在线不卡| 精品人妻熟女av久视频| 老司机影院成人| 黑丝袜美女国产一区| 不卡视频在线观看欧美| 亚洲欧美成人精品一区二区| 久久这里有精品视频免费| 高清毛片免费看| 纯流量卡能插随身wifi吗| 成人美女网站在线观看视频| 午夜免费观看性视频| 欧美精品人与动牲交sv欧美| 各种免费的搞黄视频| 一级二级三级毛片免费看| 精品一区在线观看国产| 一级毛片电影观看| videossex国产| 久久久久久久大尺度免费视频| 国产精品不卡视频一区二区| 赤兔流量卡办理| 少妇的逼水好多| 18+在线观看网站| 亚洲精华国产精华液的使用体验| 18+在线观看网站| 日日摸夜夜添夜夜添av毛片| 亚洲欧美成人精品一区二区| 亚洲男人天堂网一区| 深夜精品福利| 国产亚洲欧美在线一区二区| 人人妻人人澡人人爽人人夜夜| 亚洲人成77777在线视频| 精品国产一区二区三区久久久樱花| 国产精品av久久久久免费| 黄频高清免费视频| 天天操日日干夜夜撸| 中国美女看黄片| 国产一区二区三区av在线| 波多野结衣一区麻豆| 久久这里只有精品19| 首页视频小说图片口味搜索| 大片电影免费在线观看免费| 午夜成年电影在线免费观看| 欧美大码av| 窝窝影院91人妻| 91九色精品人成在线观看| 热99国产精品久久久久久7| 曰老女人黄片| 午夜两性在线视频| 久久久国产欧美日韩av| 日韩电影二区| 欧美老熟妇乱子伦牲交| 在线天堂中文资源库| 午夜福利在线免费观看网站| 99国产综合亚洲精品| 大型av网站在线播放| 18在线观看网站| 亚洲天堂av无毛| 日韩制服骚丝袜av| 久久ye,这里只有精品| 777米奇影视久久| 国产人伦9x9x在线观看| 亚洲欧美日韩高清在线视频 | 各种免费的搞黄视频| 久久午夜综合久久蜜桃| 丰满迷人的少妇在线观看| 淫妇啪啪啪对白视频 | 亚洲精品在线美女| 欧美日韩亚洲高清精品| 韩国精品一区二区三区| 国产亚洲av片在线观看秒播厂| 啦啦啦免费观看视频1| 久久精品熟女亚洲av麻豆精品| 人妻 亚洲 视频| 另类亚洲欧美激情| 欧美黄色片欧美黄色片| 欧美日韩成人在线一区二区| 成人手机av| 日韩有码中文字幕| 午夜老司机福利片| 各种免费的搞黄视频| 777米奇影视久久| 欧美日韩亚洲国产一区二区在线观看 | 一级,二级,三级黄色视频| 欧美精品人与动牲交sv欧美| 一级,二级,三级黄色视频| 嫁个100分男人电影在线观看| 欧美精品一区二区免费开放| 别揉我奶头~嗯~啊~动态视频 | 国产精品国产三级国产专区5o| 久久久久久免费高清国产稀缺| 三上悠亚av全集在线观看| 国产无遮挡羞羞视频在线观看| 成人av一区二区三区在线看 | 18在线观看网站| 麻豆国产av国片精品| 亚洲av男天堂| 一级毛片电影观看| 久久国产亚洲av麻豆专区| 在线 av 中文字幕| 黑人巨大精品欧美一区二区蜜桃| 日韩熟女老妇一区二区性免费视频| 亚洲天堂av无毛| 午夜免费鲁丝| 精品人妻熟女毛片av久久网站| 满18在线观看网站| 成人国语在线视频| 男男h啪啪无遮挡| 成人av一区二区三区在线看 | 国产亚洲欧美精品永久| 大片免费播放器 马上看| 下体分泌物呈黄色| 极品少妇高潮喷水抽搐| 美女高潮喷水抽搐中文字幕| 亚洲av日韩在线播放| 欧美日韩视频精品一区| 国产精品欧美亚洲77777| 一本综合久久免费| 国产一级毛片在线| 少妇被粗大的猛进出69影院| 99精国产麻豆久久婷婷| 在线观看免费午夜福利视频| 脱女人内裤的视频| 国产日韩欧美亚洲二区| 国产有黄有色有爽视频| 国产麻豆69| 男女高潮啪啪啪动态图| 窝窝影院91人妻| 国产精品久久久久久精品电影小说| 80岁老熟妇乱子伦牲交| 亚洲精华国产精华精| 午夜久久久在线观看| 91精品伊人久久大香线蕉| 97精品久久久久久久久久精品| 午夜免费成人在线视频| 人妻一区二区av| 国产成人精品久久二区二区91| 免费在线观看黄色视频的| 精品一区在线观看国产| av电影中文网址| 亚洲精品国产区一区二| 一区福利在线观看| 国产男女超爽视频在线观看| 大香蕉久久成人网| 国产精品一区二区在线不卡| 18禁观看日本| 这个男人来自地球电影免费观看| 黄色a级毛片大全视频| 新久久久久国产一级毛片| 69av精品久久久久久 | 亚洲avbb在线观看| 一边摸一边抽搐一进一出视频| 亚洲欧美精品综合一区二区三区| 久热这里只有精品99| 久久九九热精品免费| 91九色精品人成在线观看| 精品亚洲成a人片在线观看| 天堂中文最新版在线下载| 爱豆传媒免费全集在线观看| 亚洲精品第二区| 男女无遮挡免费网站观看| 国产在线视频一区二区| 国产精品成人在线| 国产成人欧美在线观看 | 欧美激情极品国产一区二区三区| 中文字幕人妻丝袜一区二区| 国内毛片毛片毛片毛片毛片| 亚洲人成电影观看| 亚洲欧洲精品一区二区精品久久久| 日韩免费高清中文字幕av| 久久久久网色| 免费在线观看视频国产中文字幕亚洲 | 80岁老熟妇乱子伦牲交| 亚洲av美国av| 俄罗斯特黄特色一大片| 久久久精品国产亚洲av高清涩受| 欧美大码av| 一级毛片电影观看| 欧美国产精品va在线观看不卡| 成年动漫av网址| 成年美女黄网站色视频大全免费| 亚洲熟女精品中文字幕| 亚洲欧美一区二区三区久久| 亚洲国产看品久久| 久久香蕉激情| 免费日韩欧美在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国内毛片毛片毛片毛片毛片| 丰满人妻熟妇乱又伦精品不卡| 久热爱精品视频在线9| 狠狠精品人妻久久久久久综合| 免费观看人在逋| 久久影院123| 亚洲精品国产av蜜桃| 欧美老熟妇乱子伦牲交| 视频在线观看一区二区三区| 啦啦啦在线免费观看视频4| 婷婷色av中文字幕| 岛国毛片在线播放| 国产激情久久老熟女| 久久久久久久大尺度免费视频| 欧美激情高清一区二区三区| av片东京热男人的天堂| 美国免费a级毛片| 一本一本久久a久久精品综合妖精| 国产色视频综合| 久久人人爽av亚洲精品天堂| 99国产极品粉嫩在线观看| 亚洲精品中文字幕一二三四区 | 精品卡一卡二卡四卡免费| 天天躁夜夜躁狠狠躁躁| 人人妻人人爽人人添夜夜欢视频| 亚洲激情五月婷婷啪啪| xxxhd国产人妻xxx| 午夜视频精品福利| 欧美午夜高清在线| 丝袜人妻中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 狠狠精品人妻久久久久久综合| 交换朋友夫妻互换小说| 精品少妇内射三级| 在线观看www视频免费| 深夜精品福利| 99热国产这里只有精品6| 侵犯人妻中文字幕一二三四区| 久久av网站| 日韩中文字幕视频在线看片| 国内毛片毛片毛片毛片毛片| 国产区一区二久久| 大片免费播放器 马上看| 亚洲国产成人一精品久久久| 国产在线免费精品| 操美女的视频在线观看| 妹子高潮喷水视频| 狠狠狠狠99中文字幕| 熟女少妇亚洲综合色aaa.| 亚洲欧美一区二区三区久久| 精品一品国产午夜福利视频| 亚洲欧美一区二区三区黑人| 别揉我奶头~嗯~啊~动态视频 | 99re6热这里在线精品视频| 最近最新免费中文字幕在线| 中国国产av一级| 国产淫语在线视频| 菩萨蛮人人尽说江南好唐韦庄| 妹子高潮喷水视频| 精品国产一区二区三区久久久樱花| svipshipincom国产片| 亚洲精品一区蜜桃| 国产成人精品在线电影| 一区二区av电影网| 欧美 亚洲 国产 日韩一| 国产精品久久久人人做人人爽| 99精品久久久久人妻精品| 少妇猛男粗大的猛烈进出视频| 精品国产一区二区三区四区第35| 91av网站免费观看| 久久影院123| 韩国高清视频一区二区三区| h视频一区二区三区| 亚洲av日韩在线播放| 成人手机av| 搡老乐熟女国产| 下体分泌物呈黄色| 久久久国产欧美日韩av| 国产成人精品久久二区二区91| 免费不卡黄色视频| 女人被躁到高潮嗷嗷叫费观| 免费高清在线观看视频在线观看| 国产在视频线精品| 黑丝袜美女国产一区| 丁香六月欧美| 成年人免费黄色播放视频| 欧美久久黑人一区二区| 国产一卡二卡三卡精品| 777米奇影视久久| 久久久久久人人人人人| 69av精品久久久久久 | 婷婷成人精品国产| 老鸭窝网址在线观看| 精品卡一卡二卡四卡免费| 欧美老熟妇乱子伦牲交| 免费av中文字幕在线| 亚洲欧美色中文字幕在线| 午夜老司机福利片| 久久久精品区二区三区| 亚洲国产欧美一区二区综合| 国产老妇伦熟女老妇高清| 日韩欧美免费精品| 午夜91福利影院| 亚洲精品国产精品久久久不卡| 国产亚洲午夜精品一区二区久久| 自线自在国产av| 国产精品成人在线| 亚洲中文字幕日韩| 中亚洲国语对白在线视频| 欧美国产精品va在线观看不卡| 少妇裸体淫交视频免费看高清 | 国产伦理片在线播放av一区| 亚洲伊人久久精品综合| 亚洲专区国产一区二区| 丝袜美腿诱惑在线| 欧美国产精品一级二级三级| 国产精品久久久久成人av| 亚洲成人免费av在线播放| 国产一区二区激情短视频 | 老汉色∧v一级毛片| 亚洲中文av在线| 桃花免费在线播放| 亚洲专区字幕在线| 日韩电影二区| 香蕉国产在线看| 一级毛片精品| 欧美av亚洲av综合av国产av| 亚洲av欧美aⅴ国产| av天堂久久9| 丝袜在线中文字幕| 十八禁网站网址无遮挡| 男人操女人黄网站| 亚洲欧洲精品一区二区精品久久久| 久久精品亚洲熟妇少妇任你| av天堂在线播放| 国产亚洲欧美在线一区二区| 亚洲熟女精品中文字幕| 免费在线观看日本一区| 亚洲第一青青草原| 日韩中文字幕视频在线看片| 国产亚洲av片在线观看秒播厂| 亚洲欧美成人综合另类久久久| 每晚都被弄得嗷嗷叫到高潮| 麻豆国产av国片精品| 麻豆乱淫一区二区| 久久久国产精品麻豆| 法律面前人人平等表现在哪些方面 | 丝袜在线中文字幕| 午夜福利免费观看在线| 精品一区在线观看国产| 啦啦啦啦在线视频资源| 搡老岳熟女国产| 亚洲成人免费电影在线观看| 一区二区日韩欧美中文字幕| 色精品久久人妻99蜜桃| 人人澡人人妻人| 日韩欧美国产一区二区入口| 国产精品自产拍在线观看55亚洲 | 日韩欧美一区视频在线观看| 欧美一级毛片孕妇| 久久久久久免费高清国产稀缺| 国产福利在线免费观看视频| 国产成人精品久久二区二区免费| 亚洲精品美女久久久久99蜜臀| 色婷婷av一区二区三区视频| 丝瓜视频免费看黄片| 69精品国产乱码久久久| 亚洲专区国产一区二区| av电影中文网址| 人妻 亚洲 视频| 极品人妻少妇av视频| 亚洲熟女毛片儿| 亚洲av美国av| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 老鸭窝网址在线观看| 在线观看免费高清a一片| 久久久欧美国产精品| 亚洲国产欧美在线一区| 男女午夜视频在线观看| 操出白浆在线播放| 老熟妇乱子伦视频在线观看 | 亚洲av男天堂| 黄网站色视频无遮挡免费观看| 香蕉丝袜av| av电影中文网址| 国产伦人伦偷精品视频| 免费日韩欧美在线观看| 永久免费av网站大全| 777久久人妻少妇嫩草av网站| 啪啪无遮挡十八禁网站| 九色亚洲精品在线播放| 久久精品成人免费网站| 日韩一区二区三区影片| 久久九九热精品免费| 飞空精品影院首页| 人人妻人人澡人人看| 丝瓜视频免费看黄片| 两人在一起打扑克的视频| 亚洲成人免费av在线播放| 91麻豆av在线| 午夜福利一区二区在线看| 久久影院123| 捣出白浆h1v1| 91精品伊人久久大香线蕉| 19禁男女啪啪无遮挡网站| 亚洲精品久久午夜乱码| 精品乱码久久久久久99久播| 精品一区二区三区av网在线观看 | 一个人免费看片子| 中文字幕色久视频| 欧美在线黄色| 亚洲国产中文字幕在线视频| av视频免费观看在线观看| 国产片内射在线| 欧美激情 高清一区二区三区| 国产精品一区二区在线不卡| 欧美激情 高清一区二区三区| 久久综合国产亚洲精品| 久久久久精品人妻al黑| 国精品久久久久久国模美| 国产精品一区二区在线观看99| 99久久99久久久精品蜜桃| 国产成人精品无人区| 亚洲一区中文字幕在线| 一级片免费观看大全| 秋霞在线观看毛片| 男女高潮啪啪啪动态图| 精品视频人人做人人爽| 青春草亚洲视频在线观看| 美女高潮喷水抽搐中文字幕| 久久99热这里只频精品6学生| 亚洲精品美女久久av网站| 精品国内亚洲2022精品成人 | 久久午夜综合久久蜜桃| 日韩欧美一区二区三区在线观看 | 新久久久久国产一级毛片| 在线av久久热| 日韩视频在线欧美| 欧美中文综合在线视频| 蜜桃在线观看..| 91老司机精品| 久久女婷五月综合色啪小说| 中文字幕精品免费在线观看视频| 后天国语完整版免费观看| 午夜影院在线不卡| 一级毛片精品| 中国美女看黄片| videos熟女内射| 我的亚洲天堂| 啪啪无遮挡十八禁网站| 国产在线一区二区三区精| 久久影院123| 51午夜福利影视在线观看| 婷婷色av中文字幕| 亚洲avbb在线观看| 国产一区二区三区在线臀色熟女 | 国产精品秋霞免费鲁丝片| 国产不卡av网站在线观看| 国产成+人综合+亚洲专区| 国产精品偷伦视频观看了| 搡老岳熟女国产| 丁香六月欧美| 最新的欧美精品一区二区| 国产亚洲精品第一综合不卡| 精品少妇久久久久久888优播| 最近最新中文字幕大全免费视频| 18禁观看日本| 国产精品一区二区在线不卡| 51午夜福利影视在线观看| 999久久久国产精品视频| 久久ye,这里只有精品| 免费观看av网站的网址| 午夜视频精品福利| 亚洲欧洲精品一区二区精品久久久| 少妇裸体淫交视频免费看高清 | 青草久久国产| 日韩电影二区| 午夜91福利影院| 亚洲熟女精品中文字幕| 久久精品aⅴ一区二区三区四区| 欧美日韩国产mv在线观看视频| 国产精品偷伦视频观看了| 久久狼人影院| 99九九在线精品视频| 日本av免费视频播放| avwww免费| 丰满少妇做爰视频| 欧美人与性动交α欧美精品济南到| 国产av一区二区精品久久| 日韩免费高清中文字幕av| 国产黄频视频在线观看| 亚洲伊人色综图| 欧美少妇被猛烈插入视频| 久久中文看片网| 亚洲 欧美一区二区三区| 国产精品av久久久久免费| 手机成人av网站| 亚洲色图 男人天堂 中文字幕|