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

    極端天氣下軌道交通水害致災因子時空分布特征

    2024-12-31 00:00:00江雨潮高燕史天根孫可天施旭
    人民長江 2024年7期
    關鍵詞:風暴潮水害極值

    摘要:

    近年來,在與氣候變化相關的多種極端天氣因素的相互作用下,洪澇事件對城市及軌道交通系統(tǒng)造成了嚴重影響。以粵港澳大灣區(qū)典型沿海城市洪澇事件為例,研究了軌道交通水害的主要致災因子(極端降雨和臺風風暴潮)及相關特征參數(shù),揭示了其時空分布特征。研究表明:① 近年來,廣州市降雨量、降雨強度均明顯增加,但降雨歷時減少,降雨整體呈現(xiàn)增多且趨于集中的極端化趨勢。廣州市降雨時間分布在年內(nèi)存在兩個極值區(qū)間,第一極值區(qū)間與廣州市“龍舟水”相關,第二極值區(qū)間主要受臺風影響。廣州市降雨空間分布不均,主要受季風及地理位置影響。② 近10 a影響廣東省的29場臺風風暴潮災害集中發(fā)生在7~10月,該時間與降雨第二極值區(qū)間重合。在對近年廣東省典型軌道交通水害事件進行總結后發(fā)現(xiàn),水害高發(fā)期與兩個極值區(qū)間重合。雖然廣東省臺風強度存在緩慢降低趨勢,但由于廣東省沿海城市遭受臺風風暴潮災害的頻次不斷上升,且粵中地區(qū)密集分布有廣東省全部的地下軌道交通系統(tǒng),未來廣東省軌道交通系統(tǒng)水害風險將隨著極端天氣時空分布特征的變化進一步加大。研究成果可為未來城市軌道交通水害致災特征研究提供依據(jù),為城市防災減災提供參考。

    關" 鍵" 詞:

    軌道交通; 洪澇災害; 致災因子; 極端降雨; 臺風風暴潮; 特征參數(shù)模型; 廣州市

    中圖法分類號: TV122

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2024.07.002

    收稿日期:

    2023-11-28;接受日期:

    2024-03-12

    基金項目:

    國家重點研發(fā)計劃項目(2022YFC3005203)

    作者簡介:

    江雨潮,男,碩士研究生,研究方向為軌道交通水害及其次生地質災害。E-mail:jiangych37@mail2.sysu.edu.cn

    通信作者:

    高" 燕,女,教授,博士生導師,主要從事地質災害、島礁鈣質砂工程性質等研究。E-mail:gaoyan25@mail.sysu.edu.cn

    Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

    文章編號:1001-4179(2024) 07-0010-10

    引用本文:

    江雨潮,高燕,史天根,等.

    極端天氣下軌道交通水害致災因子時空分布特征

    [J].人民長江,2024,55(7):10-19.

    0" 引 言

    20世紀以來,氣候變化加上高強度的人類活動,改變了區(qū)域或局域尺度的水文過程和氣候過程,致使全球極端強降雨等極端氣候條件突變式發(fā)生[1]。在洪水災害中,一部分水體侵入地下空間,如軌道交通、購物中心等,威脅民眾生命安全[2]。因此,城市軌道交通在給民眾帶來了極大便利的同時,其安全問題也越來越緊迫。近年來,城市內(nèi)強度大、歷時短的極端天氣事件頻發(fā)[3]。根據(jù)IPCC的預測,東亞地區(qū)遭遇強烈熱帶氣旋的頻率和其登陸帶來的降水強度都可能增加[4]。2021年7月河南省遭遇歷史罕見特大暴雨,鄭州市因災死亡失蹤380人,多處地鐵軌道設施遭受水淹無法運行;廣州市2020年5月22日暴雨導致城市嚴重內(nèi)澇,造成4人死亡[5];2018年臺風“山竹”、2017年臺風“天鴿”先后影響珠江口,均造成了破紀錄的高潮位[6]。

    在這種極端復合天氣事件導致復合洪水風險越發(fā)嚴重的背景下,國內(nèi)外學者開始對極端天氣事件因子的時空分布特征展開分析,Konapala等[7]研究表明全球9個區(qū)域中有4個區(qū)域的降雨量呈現(xiàn)變化增加趨勢??梢钥闯觯谌虺叨壬?,極端天氣事件顯著增加的區(qū)域多于明顯減少的區(qū)域;在大洲尺度上,北美洲、歐洲、大洋洲和亞洲東部等地區(qū)甚至在總降水量減少的條件下,極端降雨頻率和強度均有所增加,且增加的趨勢將進一步加?。?-11]。針對國內(nèi)極端天氣的時空分布特征,目前研究多基于地方雨量站[12-15]、水文站[15-19]及驗潮站[20-23]歷年數(shù)據(jù)的統(tǒng)計分析,結果表明全國降雨普遍存在頻次增加、強度增大的趨勢[24-26],空間分布上表現(xiàn)出城市多、城郊少、沿海多、內(nèi)陸少、南部多、北部少的特點[27-30]。中國沿海地區(qū)的風暴潮災害則表現(xiàn)出強度波動幅度較小,但頻次不斷增加,且隨著城市發(fā)展及海平面的不斷升高,其危害有不斷加大的趨勢[31-34]。目前,針對影響軌道交通系統(tǒng)水淹事件的研究相對較少[35-38],且缺乏極端天氣事件對軌道交通系統(tǒng)造成危害的分析[39-42],關于軌道交通水害的研究多為一般暴雨條件下的水害分析[43-45],而對于不同極端天氣條件下,軌道交通水害致災因子的時空分布特征尚不清晰。

    研究軌道交通水害的致災因子,并了解各水害因子的特征參數(shù),對其發(fā)展趨勢及時空分布特征進行數(shù)字化表征是對國家“十四五”計劃的積極響應,也是符合《國務院辦公廳關于加強城市內(nèi)澇治理的實施意見》要求的有效措施[45-46]。本文將分析影響城市、軌道交通洪澇災害的致災因子組成及基本特征,重點從降雨、臺風風暴潮等城市軌道交通水害影響因素進行統(tǒng)計分析,以期為未來城市軌道交通水害致災水體組成、水淹影響過程及致災特征研究提供參考。

    1" 數(shù)據(jù)來源及方法

    本文以廣州市及廣東省沿海地區(qū)為研究區(qū)域,通過數(shù)據(jù)統(tǒng)計和趨勢擬合構成參數(shù)模型,分析廣東省極端天氣條件下軌道交通水害致災因子的發(fā)展趨勢。其中,降雨數(shù)據(jù)來源于“廣州天氣”天氣過程實況數(shù)據(jù)(http:∥www.tqyb.com.cn/gz/weatherLive/tqgc/)、《2007~2022年廣州市氣候公報》以及廣州市各區(qū)縣氣象站2007~2022年逐時降雨數(shù)據(jù);臺風風暴潮和天文潮數(shù)據(jù)來源于廣東省自然資源廳《2013~2015年廣東海洋災害公報》及廣東省沿海18個潮位站(南渡站、廣州站、汕尾站等)的逐時潮位數(shù)據(jù),各潮位站的具體區(qū)位如表1所列。

    2" 軌道交通水害致災因子組成

    本研究為實現(xiàn)當水害形成時快速確定不同致災因子對軌道交通的影響,采用的研究方法包括:① 構建致災因子組成基本單元集;② 確定各致災因子對水害、水體組成的貢獻方式及特征參數(shù);③ 確定各致災因子的數(shù)字化參數(shù)特征[47-48]。根據(jù)GB 50014-2021《室外排水設計標準》中城市雨水設計流量計算規(guī)范,城市軌道交通水害過程中,致災水體主要來源于降雨。同時,結合沿海地區(qū)極端臺風風暴潮天氣對城市洪澇過程的影響,主要選取極端降雨與臺風風暴潮作為軌道交通水害的主要致災因子,并對其特征參數(shù)及時空分布特征展開研究。

    通過實際案例分析及現(xiàn)場、文獻調研表明,降雨是引發(fā)城市軌道交通水淹災害最直接的因素[49],確定其影響水害形成的特征參數(shù)有利于了解極端降雨的發(fā)展趨勢及其對軌道交通水淹災害的影響過程,是實現(xiàn)軌道交通水害預測預報的關鍵環(huán)節(jié)。采用平均降雨強度作為極端降雨影響軌道交通水害的主要特征參量,同時對其他降雨相關影響因素的變化特征進行分析。臺風風暴潮作為引發(fā)沿海城市軌道交通水淹事件的主要驅動因素,具有突發(fā)性強、影響范圍廣以及危害程度大的特點,往往數(shù)小時內(nèi)就能釀成巨大災害,嚴重影響人民群眾生命財產(chǎn)安全。采用風暴增水及最大風暴增水值作為臺風風暴潮的特征參數(shù)[27-31]。主要研究討論的軌道交通致災因子特征參數(shù)及相關影響因素的具體表征如表2所列。

    3" 極端降雨時空分布特征

    廣州市為中國軌道交通系統(tǒng)發(fā)展速度最快、覆蓋范圍最廣的城市之一,以廣州市為例,對其極端降雨時空分布特征進行統(tǒng)計分析。

    3.1" 廣州市降雨時間分布特征

    3.1.1" 降雨相關參數(shù)變化趨勢

    近年來,廣州市每年都會發(fā)生較為嚴重的水淹災害。經(jīng)統(tǒng)計,2007~2022年間廣州市范圍內(nèi)遭遇暴雨及暴雨以上的極端降雨事件共計212次,其中大部分的極端降雨事件都引發(fā)了一定程度的城市內(nèi)澇。本文對廣州市近16 a內(nèi)212次降雨事件的降雨特征參數(shù)進行統(tǒng)計分析,如圖1所示。結果表明:近16 a,廣州市降雨事件的過程降雨量、最大降雨強度以及平均降雨強度都呈增加趨勢,降雨歷時呈現(xiàn)減少趨勢。降雨總體上表現(xiàn)為單次降雨事件的強度增大、歷時減短,降雨次數(shù)減少且時間上趨于集中的態(tài)勢。

    3.1.2" 平均降雨強度分布規(guī)律及年內(nèi)變化趨勢

    為進一步明確降雨特征參數(shù)的數(shù)字化特征,以2007~2022年內(nèi)廣州市極端降雨事件的平均降雨強度為代表,對每年極端降雨事件的平均降雨強度分布規(guī)律進行分析。通過研究過去16 a 廣州市極端降雨事件平均降雨強度的年內(nèi)分布規(guī)律可以看出(圖2),大部分年份的平均降雨強度在時間分布上存在兩個極值,且極值兩側的數(shù)值分布規(guī)律表現(xiàn)出一定的正態(tài)分布規(guī)律,則可應用高斯函數(shù)表述各極值的正態(tài)分布特征?;诟咚购瘮?shù)建立的參數(shù)模型如公式(1)、(2)所示。

    P=P0+(P0-Pc)×e-2×x-xcw2

    (1)

    w=Pc-P02ln4

    (2)

    式中:P表示廣州市發(fā)生極端降雨事件時的平均降雨強度,mm/h;

    x表示時間,d;

    P0表示廣州市極端降雨事件的基礎降雨強度,mm/h;

    Pc表示平均降雨強度極值區(qū)間內(nèi)的極值,mm/h;

    xc表示平均降雨強度極值出現(xiàn)的時間,d;

    w表示鐘狀高斯函數(shù)的寬度。

    對兩個極值出現(xiàn)的時間進行分析后發(fā)現(xiàn),廣州市平均降雨強度的第一個極值出現(xiàn)在4~5月,該時間與廣州市前汛期“龍舟水”的時間相重合[5]。第二個極值出現(xiàn)在7~10月,該時間段內(nèi),廣東省頻繁遭遇臺風災害,廣州市受臺風帶來的風雨影響易形成極端降雨,因此可以推測該時間段內(nèi)的極值與臺風的頻繁登陸有關。由于廣州市極端降雨事件的平均降雨強度在年內(nèi)存在兩個極值區(qū)間,因此,分別將兩個極值Pc1與Pc2代入?yún)?shù)模型中可以得到不同極值區(qū)間內(nèi)平均降雨強度年內(nèi)變化趨勢表達式。

    綜上,為實現(xiàn)對平均降雨強度P時間分布模型的建立,還需要確定的模型參數(shù)包括:基礎降雨強度P0、平均降雨強度極值Pci及各極值出現(xiàn)的時間xci。下文對各模型參數(shù)情況作進一步討論。

    3.1.3" 平均降雨強度極值年際變化趨勢

    為構建廣州市未來遭遇極端降雨事件時的降雨特征參數(shù)模型,分別對兩個極值區(qū)間內(nèi)的模型參數(shù)進行統(tǒng)計分析,即對兩個極值區(qū)段內(nèi)的平均降雨強度極值(Pc1,Pc2)及極值出現(xiàn)的時間(xc1,xc2)進行統(tǒng)計分析。如圖3所示,分析平均降雨強度兩個極值的年際變化趨勢可以看出,各極值區(qū)間內(nèi)極值的年際變化呈現(xiàn)增長趨勢,符合冪增長。且在對第二極值區(qū)間內(nèi)的極值依據(jù)“該年份粵中地區(qū)是否直接遭受臺風風暴潮登陸影響”進行劃分后,發(fā)現(xiàn)第二極值區(qū)間內(nèi)的極值變化趨勢按上限(直接影響)、下限(非直接影響)分別擬合得到的R2值明顯優(yōu)于未劃分時的擬合結果。因此,可以得出結論,第二極值區(qū)間內(nèi)參數(shù)在分布上明顯存在上限與下限,上限的出現(xiàn)受該年廣東省臺風災害是否直接登陸粵中地區(qū)影響。

    根據(jù)擬合結果建立廣州市平均降雨強度極值冪增長模型如公式(3)~(5)所示,對于第一極值區(qū)間內(nèi)的平均降雨強度極值直接使用公式(3)進行模擬。對于第二極值區(qū)間內(nèi)的極值,應分別討論該年7~10月臺風風暴潮是否直接對廣州市產(chǎn)生影響,若臺風直接登陸粵中地區(qū),則選取上限對應的公式(4)進行模擬,若未直接對廣州市產(chǎn)生影響,則應用公式(5)進行模擬。

    Pc1=0.026X2-102.66X+101832.8(3)

    Pc2=0.4X2-1600.21X+1609902.8(4)

    Pc2=0.02X2-78.82X+78526.3" (5)

    式中:Pc1,Pc2分別為廣州市極端降雨事件平均降雨強度第一、第二極值區(qū)間的極值,mm/h;X為時間,a。

    對剩余的模型參數(shù)(xc1,xc2,P0)進行趨勢分析,結果如圖4所示。可見,xc與P0均存在一定的增長趨勢但并不明顯,未來在進行模型建立時可以根據(jù)實際情況取值。

    本節(jié)通過對2007~2022年廣州市極端降雨事件的平均降雨強度模型相關參數(shù)的變化趨勢分析,根據(jù)高斯函數(shù)建立了廣州市極端降雨時間分布特征模型,實現(xiàn)了對廣州市未來極端降雨事件特征參數(shù)發(fā)展態(tài)勢的數(shù)字化模擬,可以為城市軌道交通水淹災害預防、預報提供支撐,對于城市發(fā)展及工程安全建設具有重要意義。

    3.1.4" 平均降雨強度與其他降雨影響因素相關性分析

    為進一步確定未來廣州市極端降雨事件的整體發(fā)展態(tài)勢,對2007~2022年廣州市遭遇的共計212場極端降雨事件的降雨歷時及最大降雨強度與平均降雨強度的相關性進行分析,如圖5所示。由圖5(a)可知,廣州市過去遭遇的極端降雨事件的降雨歷時普遍小于24 h,且隨著平均降雨強度增大,其所對應的降雨歷時減小,減小趨勢符合指數(shù)降低趨勢。依據(jù)指數(shù)函數(shù)構建平均降雨強度-降雨歷時模型,結果如公式(6)所示。形成該趨勢的主要原因是近些年降雨存在極端化的發(fā)展趨勢,部分極端降雨事件表現(xiàn)出降雨強度大、降雨歷時短的特點。即在極短的時間內(nèi)產(chǎn)生大量的降雨量,從而導致降雨歷時短的同時平均降雨強度遠大于正常值的情況出現(xiàn)。

    P=4.59+59.6×e-T4.37

    (6)

    式中:P表示平均降雨強度,mm/h;

    T表示降雨歷時,h。

    平均降雨強度與最大降雨強度的相關性分析如圖5(b)所示,當最大降雨強度小于50 mm/h時,平均降雨強度與最大降雨強度兩者基本滿足1∶1的關系,兩者差異較小,表明該情況下的降雨過程對應的降雨歷時短,較大的降雨量在1 h左右全部降落到城市內(nèi),形成極端降雨天氣;當最大降雨強度增大到50~200 mm/h時,此時的降雨強度已經(jīng)超過了大部分城市的最大防洪能力。而在該范圍內(nèi),平均降雨強度普遍分布在0.09倍的最大降雨強度附近,這表示即使在較長降雨歷時的降雨過程內(nèi),也有可能出現(xiàn)降雨量集中在極短時間內(nèi)降落的極端降雨情況,并且這種情況一旦出現(xiàn),所產(chǎn)生的降雨強度將更大,同時由于降雨過程前期已經(jīng)累計了一定的降雨量,因此對城市產(chǎn)生的排洪壓力也將增大,對城市安全造成嚴重威脅。

    3.2" 廣州市降雨空間分布特征

    為研究廣州市降雨空間變化規(guī)律及其對軌道交通造成的風險,通過2007~2022年廣州市遭遇的212次極端降雨事件的雨量站記錄,對每次降雨事件所影響的主要區(qū)縣進行統(tǒng)計,并應用ArcGIS將廣州市各區(qū)縣遭遇極端降雨事件的總頻次按照顏色深淺以及標記大小的方式標注在廣州市地圖上,如圖6所示。由圖6可知,近16 a廣州市降雨空間分布不均,廣州市降雨主要分布在受氣候和季風控制的市區(qū)北部,如從化區(qū)、增城區(qū)等,以及受臺風風暴潮影響的市區(qū)南部,如南沙區(qū)等。將降雨空間分布圖與軌道交通分布圖對比后可以發(fā)現(xiàn),廣州市極端降雨高頻地區(qū)的軌道交通分布并不密集,廣州市軌道交通水害受極端降雨空間分布影響小。

    因此,廣州市降雨空間分布受地理位置的影響更大,與中心城區(qū)城市化發(fā)展及地下軌道交通建設的相關性低。各主要降雨地區(qū)遭遇極端降雨的頻次隨著大氣候背景下降雨及臺風風暴潮頻次的增加而增加。

    4" 風暴潮時空分布特征

    4.1" 廣東省風暴潮時間分布特征

    4.1.1" 風暴潮年內(nèi)時間分布特征

    經(jīng)統(tǒng)計,2013~2022年明顯影響廣東省的臺風共29場。對臺風登陸的時間進行統(tǒng)計,如圖7所示,結果表明7~10月為廣東省臺風風暴潮高發(fā)期,這與上述的降雨特征參數(shù)第二極值區(qū)間重合。臺風風暴潮的突發(fā)可能帶來極端降雨、超警戒潮位、河流洪水等次生災害,給廣東省帶來洪澇風險。由圖7可以看出,廣東省臺風風暴潮年內(nèi)登陸時間呈現(xiàn)出以8~9月為中心的正態(tài)分布趨勢,為更好地模擬未來廣東省臺風風暴潮年內(nèi)登陸時間及頻次,應用高斯函數(shù)構建廣東省臺風風暴潮時間分布模型,如公式(7)所示。

    N=-0.04+10.53π/2 e(x′-8.52)22.81

    (7)

    式中:N表示廣東省臺風風暴潮登陸頻次;x′表示時間,月。

    據(jù)統(tǒng)計,過去10 a影響廣東省的29場臺風中,恰逢天文大潮的臺風有8場,占總數(shù)的27.6%,可見,無論臺風登陸時是否為天文大潮期,風暴潮都可能產(chǎn)生較高的風暴潮增水從而對城市造成洪澇風險。

    近些年,廣東省在極端天氣的影響下,軌道交通水害事件頻發(fā)。2016年5月,廣州地鐵6號線長湴站發(fā)生積水倒灌;2017年6月,受臺風“苗柏”影響,深圳地鐵車公廟站積水關閉;2020年5月,廣州地鐵13號線官湖站被淹;2021年7月廣州地鐵21號線神舟路站由于出入口擋水墻倒塌導致地鐵站進水;2023年9月,臺風“??钡顷懘鬄硡^(qū),導致廣州地鐵市橋站、深圳地鐵黃閣坑站、榆園站及香港地鐵黃大仙站發(fā)生嚴重積水入侵。通過前文的分析可知,廣東省極端天氣在時間分布上存在兩個極值區(qū)間,分別為每年4~5月的第一極值區(qū)間及7~10月的第二極值區(qū)間??梢钥闯?,以上軌道交通水害事件均發(fā)生在兩個極端天氣極值區(qū)間內(nèi),且受到臺風、降雨的嚴重影響。

    4.1.2" 廣東省風暴潮年際時間分布特征

    為明確廣東省臺風風暴潮未來發(fā)展趨勢,對2013~2022年影響廣東省的29場臺風風暴潮的風暴潮增水及臺風的中心氣壓、中心最大風速進行統(tǒng)計,如圖8所示??梢钥闯?,廣東省在過去10 a中遭受的臺風中心氣壓存在一定的上升趨勢,臺風的中心最大風速存在一定的下降趨勢,臺風風暴潮災害所產(chǎn)生的風暴潮增水存在一定的下降趨勢。以上趨勢均表明廣東省近10 a內(nèi)遭遇的臺風風暴潮災害強度正在緩慢減弱。

    通過分析臺風風暴潮相關參數(shù)的年際變化趨勢,建立臺風風暴潮時間分布模型可以實現(xiàn)對未來廣東省臺風風暴潮災害影響因素、過程的模擬,從而更好地幫助廣東省沿海城市預測、預報、預防由臺風風暴潮誘發(fā)的城市洪澇及軌道交通水淹災害。

    4.2" 廣東省風暴潮空間分布特征

    依據(jù)廣東省沿海各區(qū)市潮位站的潮位變化數(shù)據(jù),對2013~2022年廣東省沿海受臺風風暴潮影響進而產(chǎn)生風暴增水的頻次進行統(tǒng)計,各驗潮站對應的頻次如圖9所示。結果表明:過去10 a對廣東省產(chǎn)生影響的29場臺風,對粵西各潮位站產(chǎn)生影響55次,對粵中產(chǎn)生影響72次,對粵東產(chǎn)生影響25次?;浳魇茱L暴潮影響頻次呈先減小后增大趨勢,粵中以及粵東各區(qū)域受風暴潮影響頻次均呈持續(xù)增加趨勢??梢酝茰y,未來廣東省沿海各區(qū)市受臺風風暴潮威脅的次數(shù)將不

    斷增加。然而,粵中等沿海地區(qū)作為廣東省經(jīng)濟發(fā)展

    最快的地區(qū),其軌道交通網(wǎng)絡的分布密度也遠遠大于

    其他內(nèi)陸城市。若未來出現(xiàn)多場臺風連續(xù)登陸的極端

    天氣情況,將會對廣東省軌道交通系統(tǒng)造成十分嚴重的危害。因此,風暴潮影響頻次的增加將會給廣東省各沿海城市洪澇帶來更大的風險,廣東省軌道交通系統(tǒng)受臺風風暴潮空間分布特征影響將十分顯著。

    5" 極端降雨與風暴潮的耦合關系

    前文分別討論了不同軌道交通水害致災因子各自

    的時空分布特征,并通過對各特征參數(shù)進行分析,構建了不同的特征參數(shù)模型。然而,由于氣候變化的影響,相關因素(如極端降雨、嚴重風暴潮等)的相互作用將

    進一步加強。極端天氣事件期間,風暴潮和強降雨同時發(fā)生,將增加粵港澳大灣區(qū)沿海城市的洪水風險。例如,臺風“天鴿”于2017年8月23日在珠三角登陸,引發(fā)嚴重風暴潮和強降雨。香港、澳門和珠海市受到高海平面、強降雨和地表洪水的嚴重影響。同時,在前文的研究中也發(fā)現(xiàn),廣州市極端降雨第二極值區(qū)間在時間上與風暴潮集中影響廣東省的時間重合,第二極值區(qū)間內(nèi)的上限值直接受臺風影響。因此,應考慮極端天氣條件下多種致災因子之間的相關性來評估其對軌道交通水害的影響程度。

    通過對2013~2022年期間的29場臺風及其對應第二極值區(qū)間內(nèi)的極端降雨事件進行統(tǒng)計,分析臺風中心氣壓及最大中心風速與極端降雨平均降雨強度及風暴潮最大風暴潮增水之間的相關性。應用二元二次方程對各參數(shù)的分布趨勢進行擬合(圖10),構建風暴潮-降雨耦合模型,如公式(8)、(9)所示。從結果中可以看出,平均降雨強度及最大風暴潮增水隨相關因素的變化趨勢相同,與臺風中心氣壓呈負相關,與臺風中心最大風速呈正相關。未來通過耦合模型可以實現(xiàn)在一定臺風強度條件下對臺風造成的極端天氣影響進行模擬。

    P=0.0062p2-12.26p-0.012v2max+0.74vmax+6045.69(8)

    emax=0.025p2-53.28p-0.18v2max+12.47vmax+27764.96(9)

    式中:P為平均降雨強度,mm/h;

    emax為最大風暴潮增水,cm;

    p為臺風中心氣壓,hPa;

    vmax為臺風最大中心風速,m/s。

    6" 結 論

    由于全球范圍內(nèi)氣候變化及城市化快速發(fā)展,極端天氣事件愈加頻繁,極端降雨強度加大以及臺風風暴潮頻發(fā)使軌道交通水淹防范情勢更加嚴峻。本文以廣州市及周邊地區(qū)為代表分析了軌道交通水害致災因子時空分布特征,并建立了特征參數(shù)數(shù)字化模型。主要結論如下:

    (1) 極端降雨和臺風風暴潮災害是引發(fā)城市軌道交通水淹的主要原因。極端降雨和風暴潮的主要致災特征參數(shù)為平均降雨強度及最大風暴潮增水等,其影響因素包括降雨量、降雨歷時、臺風中心氣壓及臺風中心最大風速等。

    (2) 近16 a廣州市212場極端降雨事件的降雨量、平均降雨強度及最大降雨強度都呈增加趨勢,降雨歷時緩慢減少。其中平均降雨強度在年內(nèi)存在兩個極值分布區(qū)間,時間上分別與廣州市“龍舟水”及臺風集中影響時間相符。各極值區(qū)間內(nèi)極值均表現(xiàn)出一定的增長趨勢,且第二極值區(qū)間內(nèi)受臺風影響存在上下限。廣州市極端降雨事件呈現(xiàn)增多且趨于集中的態(tài)勢。對廣東省典型的軌道交通水害事件進行總結后發(fā)現(xiàn),水害事件集中發(fā)生在兩個極值區(qū)間內(nèi)。廣州市軌道交通系統(tǒng)水害受極端降雨時間分布影響顯著,而受極端降雨空間分布影響小,降雨頻次在空間上分布不均,主要受季風和地理位置影響。

    (3) 近29場臺風風暴潮的數(shù)據(jù)表明,廣東省風暴潮增水值與臺風中心氣壓及最大中心風速變化密切相關。臺風風暴潮災害影響廣東省的頻次正在不斷上升,粵中區(qū)域分布有廣東省已建成的全部軌道交通系統(tǒng),軌道交通受風暴潮空間分布的影響顯著。風暴潮引起的平均降雨強度及風暴潮增水隨著臺風中心氣壓的增大而減小,隨著臺風中心最大風速的增大而增大,各特征參數(shù)與影響因素間存在明確的數(shù)字化特征。

    參考文獻:

    [1]" HALLEGATTE S,GREEN C,NICHOLLS R J,et al.Future flood losses in major coastal cities[J].Nature Climate Change,2013,3(9):802-806.

    [2]" OZAKI T,ASANO N,ISHIGAKI T,et al.Characteristics of underground mall inundation by flush flood in short time heavy rainfall[J].Journal of Japan Society of Civil Engineers,2014,70:I_1417-I_1422.

    [3]" QIANG Y,ZHANG L,HE J,et al.Urban flood analysis for Pearl River Delta cities using an equivalent drainage method upon combined rainfall-high tide-storm surge events[J].Journal of Hydrology,2021,597:126293.

    [4]" BERNSTEIN L,BOSCH P,CANZIANI O,et al.IPCC,2007:climate change 2007[M].Geneva:IPCC,2008.

    [5]" 張之琳,邱靜,程濤,等.粵港澳大灣區(qū)城市洪澇問題及其分析[J].水利學報,2022,53(7):823-832.

    [6]" 賈寧,劉強,石先武,等.基于現(xiàn)場調查的臺風“天鴿”(1713)和臺風“山竹”(1822)風暴潮災害影響和致災對比分析[J].海洋預報,2022,39(5):94-99.

    [7]" KONAPALA G,MISHRA A K,WADA Y,et al.Climate change will affect global water availability through compounding changes in seasonal precipitation and evaporation[J].Nature Communications,2020,11(1):3044.

    [8]" SCHUSTER Z T,POTTER K W,LIEBL D S.Assessing the effects of climate change on precipitation and flood damage in Wisconsin[J].Journal of Hydrologic Engineering,2012,17(8):888-894.

    [9]" SHAO Z,GAO F,YANG S L,et al.A new semiparametric and EEMD based framework for mid-term electricity demand forecasting in China:hidden characteristic extraction and probability density prediction[J].Renewable and Sustainable Energy Reviews,2015,52:876-889.

    [10]孔鋒,史培軍,方建,等.全球變化背景下極端降水時空格局變化及其影響因素研究進展和展望[J].災害學,2017,32(2):165-174.

    [11]MIN S K,ZHANG X,ZWIERS F W,et al.Human contribution to more-intense precipitation extremes[J].Nature,2011,470(7334):378-381.

    [12]ZHANG Y,XU Y,DONG W,et al.A future climate scenario of regional changes in extreme climate events over China using the PRECIS climate model[J].Geophysical Research Letters,2006,33(24):L24702.

    [13]顧西輝,張強,孔冬冬.中國極端降水事件時空特征及其對夏季溫度響應[J].地理學報,2016,71(5):718-730.

    [14]殷水清,王楊,謝云,等.中國降雨過程時程分型特征[J].水科學進展,2014,25(5):617-624.

    [15]杜經(jīng)緯,王昊,王浩,等.短歷時雨型對城市排水管網(wǎng)運行狀態(tài)影響研究[J].人民長江,2023,54(4):108-113,162.

    [16]成丹,陳正洪.湖北宜昌市區(qū)暴雨雨型的演變特征[J].干旱氣象,2017,35(2):225-231.

    [17]何玉秀,許有鵬,李子貽,等.城鎮(zhèn)化對極端降水的影響及其貢獻率研究:以太湖平原地區(qū)為例[J].湖泊科學,2022,34(1):262-271.

    [18]崔婷婷,王銀堂,胡慶芳.城鎮(zhèn)化進程中太湖流域降水演變規(guī)律研究[J].江蘇水利,2019(8):1-6.

    [19]于志慧,劉鵬飛,何玉秀,等.太湖平原高度城鎮(zhèn)化下降雨特征時空變化研究[J].湖泊科學,2022,34(6):2134-2143.

    [20]孫周亮,鮑振鑫,舒章康,等.南京市近25年短歷時暴雨雨型特性[J].水文,2019,39(5):78-83.

    [21]李潔,徐強強,王正華.考慮降雨空間異質性的洪水預報模型研究[J].人民長江,2023,54(增1):11-15.

    [22]宋曉猛,張建云,孔凡哲,等.北京地區(qū)降水極值時空演變特征[J].水科學進展,2017,28(2):161-173.

    [23]黃國如,陳易偲,姚芝軍.高度城鎮(zhèn)化背景下珠三角地區(qū)極端降雨時空演變特征[J].水科學進展,2021,32(2):161-170.

    [24]譚暢,孔鋒,郭君,等.1961~2014年中國不同城市化地區(qū)暴雨時空格局變化:以京津冀、長三角和珠三角地區(qū)為例[J].災害學,2018,33(3):132-140.

    [25]于鑫,鄭騰飛,劉顯通,等.廣東省極端降水的時空分布特征[J].廣東氣象,2022,44(2):6-9.

    [26]程遙,紀彭.深圳市“4·11”暴雨洪水分析 [J].水利水電快報,2022,43(5):8-11,19.

    [27]葉榮輝,錢燕,孔俊,等.珠江三角洲大系統(tǒng)風暴潮數(shù)學模型建立與驗證[J].人民長江,2013,44(21):76-80.

    [28]龐古乾,黃玉明,何健,等.珠海市臺風及其次生風暴潮特征分析與評估[J].熱帶氣象學報,2020,36(1):42-50.

    [29]韓晶.臺風山竹和天鴿對珠海沿海風暴潮增水影響[J].吉林水利,2019(447):47-49,53.

    [30]史鍵輝,王永信.珠江口縱深的風暴潮和增水特征[J].海洋預報,2000(4):47-51.

    [31]張敏,羅軍,胡金磊,等.雷州市沿海風暴潮淹沒危險性評估[J].熱帶海洋學報,2019,38(2):1-12.

    [32]李穎,方偉華,林偉,等.可能最大風暴潮風險評估中各等級熱帶氣旋設定方法[J].海洋科學,2014,38(4):71-80.

    [33]周才揚,殷成團,章衛(wèi)勝,等.長江口臺風影響能力研究[J].人民長江,2021,52(1):27-30.

    [34]劉秋興,傅賜福,李明杰,等.“天鴿”臺風風暴潮預報及數(shù)值研究[J].海洋預報,2018,35(1):29-36.

    [35]ISHIGAKI T,TODA K,BABA Y,et al.Experimental study on evacuation from underground space by using real size models[J].Proceedings of Hydraulic Engineering,2006,50:583-588.

    [36]ISHIGAKI T,KAWANAKA R,ONISHI Y,et al.Assessment of safety on evacuating route during underground flooding[C]∥Proceedings of 16th IAHR-APD Congress and 3rd Symposium of IAHR-ISHS,2009:141-146.

    [37]TODA K,KAWAIKE K,YONEYAMA N,et al.Underground inundation analysis by integrated urban flood model[C]∥Proceedings of 16th IAHR-APD Congress and 3rd Symposium of IAHR-ISHS,2009:166-171.

    [38]YONEYAMA N,TODA K,AIHATA S,et al.Numerical analysis for evacuation possibility from small underground space in urban flood[C]∥Proceedings of 16th IAHR-APD Congress and 3rd Symposium of IAHR-ISHS,2009:107-112.

    [39]QUAN R,ZHANG L,LIU M,et al.Risk assessment of rainstorm waterlogging on subway in central urban area of Shanghai,China based on scenario simulation[C]∥2011 19th International Conference on Geoinformatics.Shanghai:IEEE,2011:1-6.

    [40]INOUE T,KAWANAKA R,ISHIGAKI T,et al.Safety evacuation from underground mall during urban flooding[J].Journal of Japan Society of Civil Engineers,2011,67(4):I_973-I_978.

    [41]HAMAGUCHI S,ISHIGAKI T,OZAKI T,et al.Study on vulnerability of mega-underground space to extreme flood[J].Journal of Japan Society of Civil Engineers,2016,72(4):I_1363-I_1368.

    [42]LYU H M,WANG G F,SHEN J,et al.Analysis and GIS mapping of flooding hazards on 10 May 2016,Guangzhou,China[J].Water,2016,8(10):447.

    [43]LYU H M,SUN W J,SHEN S L,et al.Flood risk assessment in metro systems of mega-cities using a GIS-based modeling approach[J].Science of the Total Environment,2018,626:1012-1025.

    [44]LI M,KWAN M P,YIN J,et al.The potential effect of a 100-year pluvial flood event on metro accessibility and ridership:a case study of central Shanghai,China[J].Applied Geography,2018,100:21-29.

    [45]中華人民共和國國民經(jīng)濟和社會發(fā)展第十四個五年規(guī)劃和2035年遠景目標綱要[EB/OL].(2021-03-13)[2023-11-28].https:∥www.gov.cn/xinwen/2021-03/13/content_5592681.html.

    [46]國務院辦公廳.國務院辦公廳關于加強城市內(nèi)澇治理的實施意見:國辦發(fā)〔2021〕11號[Z].北京:國務院辦公廳,2021.

    [47]FORERO-ORTIZ E,MARTNEZ-GOMARIZ E,CAAS PORCUNA M.A review of flood impact assessment approaches for underground infrastructures in urban areas:a focus on transport systems[J].Hydrological Sciences Journal,2020,65(11):1943-1955.

    [48]莫偉麗.地鐵車站水侵過程數(shù)值模擬及避災對策研究[D].杭州:浙江大學,2010.

    [49]MAIDMENT D.水文學手冊[M].張建云,李紀生,等,譯.北京:科學出版社,2002.

    (編輯:謝玲嫻)

    Temporal and spatial distribution characteristics of disaster-causing factors for rail transit flooding under extreme conditions

    JIANG Yuchao1,GAO Yan1,2,SHI Tiangen1,SUN Ketian1,SHI Xu1

    (1.School of Earth Sciences and Engineering,Sun Yat-sen University,Zhuhai 519000,China;" 2.Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai),Zhuhai 519000,China)

    Abstract:

    In recent years,under the interaction of multiple extreme weather factors related to climate change,flooding events have caused serious impacts on cities and rail transit systems.Taking a typical coastal city Guangzhou City in the Guangdong-Hong Kong-Macao Greater Bay Area as a case study,we examined the main contributing factors (extreme rainfall and typhoon storm surge) and related parameters of rail transit flooding,revealing their temporal and spatial distribution characteristics.The study results showed that: ① The amount and intensity of rainfall in Guangzhou City have increased significantly,while the duration of rainfall has decreased,leading to an overall increase in rainfall and a trend toward more extreme precipitation events.There are two extreme intervals of rainfall distribution in Guangzhou City during the year.The first interval is associated with the \"Dragon-Boat-Precipitation,\" while the second interval is mainly influenced by typhoons.The spatial distribution of rainfall in Guangzhou City is uneven,primarily due to the monsoon and geographic location.② The 29 typhoons and storm surges that affected Guangdong Province in the past 10 years occurred between July and October,coinciding with the second extreme rainfall interval.Analyzing typical rail transit water damage events in Guangdong Province in recent years revealed that the high incidence of water damage aligns with these two extreme intervals.Although the intensity of typhoons in Guangdong Province is slowly decreasing,the frequency of typhoons and storm surges in its coastal cities is increasing.Given that all rail transit systems are densely distributed in the central region of Guangdong Province,the risk of water damage to these systems is expected to increase in the future due to changes in the temporal and spatial distribution characteristics of extreme weather.The results of the study can serve as a foundation for future research on the characteristics of urban rail transit water damage and as a reference for urban disaster prevention and mitigation.

    Key words:

    rail transit; flood disasters; disaster-causing factors; extreme rainfall; typhoon storm surge; characteristic parametric model; Guangzhou City

    猜你喜歡
    風暴潮水害極值
    極值點帶你去“漂移”
    2012年“蘇拉”和“達維”雙臺風影響的近海風暴潮過程
    海洋通報(2021年2期)2021-07-22 07:55:24
    極值點偏移攔路,三法可取
    防范未來風暴潮災害的綠色海堤藍圖
    科學(2020年4期)2020-11-26 08:27:00
    基于多變量LSTM神經(jīng)網(wǎng)絡模型的風暴潮臨近預報
    海洋通報(2020年6期)2020-03-19 02:10:18
    一類“極值點偏移”問題的解法與反思
    匹配數(shù)為1的極值2-均衡4-部4-圖的結構
    礦井水害立體防治技術體系
    科技傳播(2014年21期)2014-11-17 17:08:30
    廣東省風暴潮時空分布特征及重點城市風暴潮風險研究
    海洋學報(2014年3期)2014-06-05 14:35:01
    淺談礦井常見水害事故及防治對策
    河南科技(2014年16期)2014-02-27 14:13:33
    日韩欧美一区二区三区在线观看 | 视频在线观看一区二区三区| 久久久久久久大尺度免费视频| 在线十欧美十亚洲十日本专区| 美女国产高潮福利片在线看| 777久久人妻少妇嫩草av网站| 午夜激情av网站| 九色亚洲精品在线播放| 欧美激情久久久久久爽电影 | 黄色a级毛片大全视频| 久久 成人 亚洲| 丝袜美足系列| 亚洲国产精品成人久久小说| 欧美日韩精品网址| 久热爱精品视频在线9| 欧美日韩av久久| 亚洲欧美精品综合一区二区三区| 精品久久久精品久久久| 亚洲一区中文字幕在线| 青青草视频在线视频观看| 国产亚洲欧美精品永久| av又黄又爽大尺度在线免费看| 欧美一级毛片孕妇| 久久久精品国产亚洲av高清涩受| 日韩有码中文字幕| av片东京热男人的天堂| 日韩制服骚丝袜av| 午夜成年电影在线免费观看| 成人18禁高潮啪啪吃奶动态图| 一级毛片精品| 99国产精品一区二区三区| 曰老女人黄片| 老熟妇仑乱视频hdxx| 天天躁夜夜躁狠狠躁躁| 国产亚洲精品第一综合不卡| 精品亚洲乱码少妇综合久久| 手机成人av网站| 99久久国产精品久久久| 日韩有码中文字幕| 久久九九热精品免费| 叶爱在线成人免费视频播放| 久久精品人人爽人人爽视色| 亚洲欧美精品自产自拍| 十八禁人妻一区二区| 国产色视频综合| 亚洲av男天堂| 免费av中文字幕在线| 国精品久久久久久国模美| 又黄又粗又硬又大视频| 国产亚洲一区二区精品| 亚洲伊人色综图| 亚洲精品久久成人aⅴ小说| 99国产精品免费福利视频| 国产一区二区三区在线臀色熟女 | 在线 av 中文字幕| 91大片在线观看| 日本欧美视频一区| 老熟妇乱子伦视频在线观看| 夜夜看夜夜爽夜夜摸| 久久香蕉精品热| 午夜福利18| av福利片在线观看| 国产精品亚洲美女久久久| 精品福利观看| 亚洲美女黄片视频| 亚洲片人在线观看| 老汉色av国产亚洲站长工具| 国产精品精品国产色婷婷| a在线观看视频网站| 久久九九热精品免费| 人妻久久中文字幕网| 国产激情久久老熟女| 最新在线观看一区二区三区| 丁香六月欧美| 午夜激情av网站| 亚洲专区中文字幕在线| 亚洲美女黄片视频| 国产精品久久久久久精品电影| 久久精品影院6| 香蕉丝袜av| a在线观看视频网站| 少妇的丰满在线观看| 天堂影院成人在线观看| 亚洲国产欧美网| 成年版毛片免费区| 亚洲精品在线美女| 老司机靠b影院| 久久人人精品亚洲av| 久久精品夜夜夜夜夜久久蜜豆 | 一个人免费在线观看电影 | 毛片女人毛片| 欧美成人免费av一区二区三区| 国产探花在线观看一区二区| 成人高潮视频无遮挡免费网站| 亚洲中文字幕一区二区三区有码在线看 | svipshipincom国产片| 黄色成人免费大全| 欧美黄色片欧美黄色片| 精品无人区乱码1区二区| 久久久久免费精品人妻一区二区| 色噜噜av男人的天堂激情| www.999成人在线观看| 国产精华一区二区三区| 国内久久婷婷六月综合欲色啪| 国产熟女xx| 亚洲成人中文字幕在线播放| 男女视频在线观看网站免费 | 久久久久久国产a免费观看| 精品一区二区三区四区五区乱码| 熟妇人妻久久中文字幕3abv| 香蕉av资源在线| 真人一进一出gif抽搐免费| 波多野结衣巨乳人妻| 午夜免费成人在线视频| 久久久久久久精品吃奶| 国产精品 欧美亚洲| 熟女电影av网| 亚洲人成77777在线视频| 亚洲免费av在线视频| 又黄又粗又硬又大视频| 男女床上黄色一级片免费看| 欧美色欧美亚洲另类二区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产成人免费| 麻豆一二三区av精品| 午夜影院日韩av| 成人午夜高清在线视频| 中文字幕久久专区| 一级a爱片免费观看的视频| 男女视频在线观看网站免费 | 叶爱在线成人免费视频播放| 欧美大码av| 亚洲色图av天堂| 日韩中文字幕欧美一区二区| 校园春色视频在线观看| 9191精品国产免费久久| 老司机靠b影院| 久久人人精品亚洲av| 亚洲精品久久成人aⅴ小说| 亚洲成人中文字幕在线播放| 十八禁网站免费在线| 一个人免费在线观看的高清视频| 国产成人精品无人区| 午夜精品在线福利| 黄色丝袜av网址大全| 两性午夜刺激爽爽歪歪视频在线观看 | 两个人的视频大全免费| avwww免费| videosex国产| 夜夜爽天天搞| 操出白浆在线播放| 九色国产91popny在线| 亚洲av五月六月丁香网| 精品人妻1区二区| 国产三级黄色录像| 亚洲五月婷婷丁香| 国内久久婷婷六月综合欲色啪| 国语自产精品视频在线第100页| 天堂动漫精品| 成人亚洲精品av一区二区| 久久久精品大字幕| 午夜激情av网站| 成人永久免费在线观看视频| 99热只有精品国产| 婷婷丁香在线五月| 少妇粗大呻吟视频| 1024视频免费在线观看| 欧美乱色亚洲激情| 亚洲精品色激情综合| 免费在线观看成人毛片| 婷婷亚洲欧美| 舔av片在线| 欧美中文综合在线视频| 波多野结衣高清无吗| 亚洲九九香蕉| 久久久久久久久免费视频了| 日韩 欧美 亚洲 中文字幕| 黄片大片在线免费观看| 日本熟妇午夜| 欧美黑人欧美精品刺激| 午夜福利成人在线免费观看| 91字幕亚洲| 两个人视频免费观看高清| 精品久久久久久久久久免费视频| 美女免费视频网站| 最近在线观看免费完整版| 亚洲中文日韩欧美视频| 精品国产超薄肉色丝袜足j| 黑人巨大精品欧美一区二区mp4| 亚洲在线自拍视频| 国产亚洲精品一区二区www| 在线观看www视频免费| 丁香欧美五月| 舔av片在线| 国产成人aa在线观看| 啦啦啦观看免费观看视频高清| 午夜两性在线视频| 国产一区二区在线观看日韩 | 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品一区二区www| 99国产极品粉嫩在线观看| 色综合站精品国产| 国产成人aa在线观看| 看片在线看免费视频| 婷婷精品国产亚洲av在线| 人成视频在线观看免费观看| 午夜日韩欧美国产| 国产精品一区二区三区四区久久| 色av中文字幕| 久久欧美精品欧美久久欧美| 一本一本综合久久| 欧美不卡视频在线免费观看 | 色老头精品视频在线观看| 亚洲熟妇中文字幕五十中出| 成人国语在线视频| 18禁观看日本| a级毛片在线看网站| 男男h啪啪无遮挡| 日韩欧美一区二区三区在线观看| 精品欧美一区二区三区在线| 99国产极品粉嫩在线观看| 俺也久久电影网| 国产成+人综合+亚洲专区| 欧美黄色片欧美黄色片| 日韩欧美在线二视频| 久久久久亚洲av毛片大全| 1024视频免费在线观看| 中文资源天堂在线| 亚洲精品粉嫩美女一区| 久久久久精品国产欧美久久久| 人人妻,人人澡人人爽秒播| 老鸭窝网址在线观看| 精品久久久久久久人妻蜜臀av| 国产高清激情床上av| 无人区码免费观看不卡| 免费在线观看影片大全网站| 一本久久中文字幕| 欧美精品亚洲一区二区| 日日干狠狠操夜夜爽| 免费一级毛片在线播放高清视频| 亚洲男人天堂网一区| www国产在线视频色| 老司机靠b影院| 欧美在线一区亚洲| 国产成+人综合+亚洲专区| 国产欧美日韩精品亚洲av| 美女扒开内裤让男人捅视频| 美女大奶头视频| 成人av在线播放网站| 国产精品98久久久久久宅男小说| 亚洲精品中文字幕一二三四区| 免费一级毛片在线播放高清视频| 丰满人妻一区二区三区视频av | 免费在线观看日本一区| 国产探花在线观看一区二区| 一边摸一边做爽爽视频免费| 国产高清视频在线观看网站| 很黄的视频免费| 白带黄色成豆腐渣| 亚洲,欧美精品.| 亚洲欧美日韩高清专用| 欧美乱色亚洲激情| 2021天堂中文幕一二区在线观| АⅤ资源中文在线天堂| 无人区码免费观看不卡| 国产精品1区2区在线观看.| 欧美一区二区国产精品久久精品 | 欧美 亚洲 国产 日韩一| 他把我摸到了高潮在线观看| 久久久久性生活片| 久久久久国产精品人妻aⅴ院| 天堂√8在线中文| 国产成人av教育| 国产视频一区二区在线看| 久久久久久大精品| 色av中文字幕| 黄色 视频免费看| 久久中文字幕一级| 亚洲男人的天堂狠狠| 每晚都被弄得嗷嗷叫到高潮| 亚洲七黄色美女视频| 两个人看的免费小视频| 老司机靠b影院| 两个人视频免费观看高清| 久久热在线av| 精品久久久久久,| 这个男人来自地球电影免费观看| 高潮久久久久久久久久久不卡| 12—13女人毛片做爰片一| 久久精品国产亚洲av高清一级| 天堂影院成人在线观看| 午夜影院日韩av| 毛片女人毛片| 在线十欧美十亚洲十日本专区| 国产av不卡久久| 美女扒开内裤让男人捅视频| 久久精品人妻少妇| 免费人成视频x8x8入口观看| 精品久久久久久久久久免费视频| 欧美+亚洲+日韩+国产| 欧美日韩亚洲综合一区二区三区_| 久久伊人香网站| 久久精品91蜜桃| 欧美日韩亚洲综合一区二区三区_| 日韩有码中文字幕| 无限看片的www在线观看| 亚洲av成人av| 法律面前人人平等表现在哪些方面| 亚洲熟女毛片儿| 一区福利在线观看| 18美女黄网站色大片免费观看| 国产探花在线观看一区二区| 日韩精品免费视频一区二区三区| 精品熟女少妇八av免费久了| 搞女人的毛片| 久久精品影院6| 日本三级黄在线观看| 91字幕亚洲| 18美女黄网站色大片免费观看| 亚洲av成人一区二区三| 黄色女人牲交| 久久久精品国产亚洲av高清涩受| 欧美 亚洲 国产 日韩一| 久久久水蜜桃国产精品网| 国产精品亚洲av一区麻豆| 999精品在线视频| 色播亚洲综合网| 亚洲av片天天在线观看| 丰满的人妻完整版| 又爽又黄无遮挡网站| 每晚都被弄得嗷嗷叫到高潮| 精品免费久久久久久久清纯| 精品久久久久久久久久免费视频| 五月伊人婷婷丁香| 国产亚洲av嫩草精品影院| 成人一区二区视频在线观看| 18禁黄网站禁片免费观看直播| 精品熟女少妇八av免费久了| 亚洲色图av天堂| 90打野战视频偷拍视频| 久久草成人影院| 国产午夜福利久久久久久| 天天添夜夜摸| 国产精品 欧美亚洲| 天天躁夜夜躁狠狠躁躁| av超薄肉色丝袜交足视频| 成人精品一区二区免费| 99国产极品粉嫩在线观看| 国模一区二区三区四区视频 | 亚洲av片天天在线观看| 两个人免费观看高清视频| 99精品久久久久人妻精品| 国产亚洲av高清不卡| 黄频高清免费视频| 久久人妻av系列| 亚洲国产精品sss在线观看| 国产黄色小视频在线观看| 色综合欧美亚洲国产小说| 曰老女人黄片| 一级a爱片免费观看的视频| 两个人视频免费观看高清| 欧洲精品卡2卡3卡4卡5卡区| 欧美三级亚洲精品| 国产成年人精品一区二区| 99久久国产精品久久久| 在线观看66精品国产| 欧美日韩亚洲综合一区二区三区_| 日本一本二区三区精品| 色在线成人网| 日本一区二区免费在线视频| 亚洲熟女毛片儿| 亚洲精品国产精品久久久不卡| 成人国语在线视频| 黄色成人免费大全| 亚洲av电影在线进入| www.熟女人妻精品国产| 亚洲电影在线观看av| 免费观看人在逋| 高潮久久久久久久久久久不卡| 国产成人系列免费观看| 无人区码免费观看不卡| 全区人妻精品视频| 两个人视频免费观看高清| 窝窝影院91人妻| 在线观看一区二区三区| 久久久精品国产亚洲av高清涩受| 欧美精品亚洲一区二区| 精品久久久久久久毛片微露脸| 狂野欧美激情性xxxx| 不卡一级毛片| а√天堂www在线а√下载| 国产99久久九九免费精品| 国产精品香港三级国产av潘金莲| 久久性视频一级片| 一级毛片高清免费大全| 亚洲成人精品中文字幕电影| av福利片在线| av福利片在线观看| 最近视频中文字幕2019在线8| 高潮久久久久久久久久久不卡| 无遮挡黄片免费观看| 狠狠狠狠99中文字幕| 免费在线观看日本一区| 中文字幕人成人乱码亚洲影| 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| 午夜a级毛片| 欧美色视频一区免费| 老熟妇乱子伦视频在线观看| 床上黄色一级片| 国产人伦9x9x在线观看| 亚洲,欧美精品.| 日韩欧美在线二视频| 99久久精品热视频| 91老司机精品| 亚洲无线在线观看| 观看免费一级毛片| 日本 av在线| 夜夜夜夜夜久久久久| 国产亚洲精品久久久久久毛片| 日韩av在线大香蕉| 欧美日韩国产亚洲二区| 久久久水蜜桃国产精品网| 丰满人妻熟妇乱又伦精品不卡| 精品久久久久久,| 黄色片一级片一级黄色片| 日韩中文字幕欧美一区二区| 俄罗斯特黄特色一大片| 成人欧美大片| 精品乱码久久久久久99久播| 精品日产1卡2卡| 国产野战对白在线观看| bbb黄色大片| 国产精品久久久久久久电影 | 欧美精品啪啪一区二区三区| 亚洲性夜色夜夜综合| 精品欧美一区二区三区在线| 午夜福利在线在线| 日韩欧美精品v在线| 国产精品久久视频播放| 精品不卡国产一区二区三区| 99精品久久久久人妻精品| 久久久久亚洲av毛片大全| 可以免费在线观看a视频的电影网站| 国产不卡一卡二| 成人18禁在线播放| 变态另类丝袜制服| 国产免费男女视频| 1024视频免费在线观看| 久久久久久九九精品二区国产 | 少妇人妻一区二区三区视频| 丰满人妻熟妇乱又伦精品不卡| 中文字幕人妻丝袜一区二区| 久久精品国产综合久久久| 岛国在线免费视频观看| 一级黄色大片毛片| 久久久国产成人免费| 国产av在哪里看| 成熟少妇高潮喷水视频| 日本精品一区二区三区蜜桃| 韩国av一区二区三区四区| 成人午夜高清在线视频| 国产一级毛片七仙女欲春2| 国产精品美女特级片免费视频播放器 | 后天国语完整版免费观看| 亚洲性夜色夜夜综合| 久久国产精品人妻蜜桃| 亚洲片人在线观看| 女人被狂操c到高潮| 亚洲国产精品999在线| 91麻豆av在线| 久久这里只有精品中国| 蜜桃久久精品国产亚洲av| 国产精品久久久久久精品电影| 制服人妻中文乱码| 精品乱码久久久久久99久播| 欧美久久黑人一区二区| 午夜免费观看网址| 日本在线视频免费播放| 此物有八面人人有两片| 国内精品久久久久久久电影| 黄色 视频免费看| 亚洲性夜色夜夜综合| 韩国av一区二区三区四区| 女生性感内裤真人,穿戴方法视频| 久久午夜综合久久蜜桃| 男人舔女人下体高潮全视频| 欧洲精品卡2卡3卡4卡5卡区| 97碰自拍视频| 免费在线观看亚洲国产| 男人舔女人的私密视频| 国产熟女午夜一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩福利视频一区二区| 日本免费一区二区三区高清不卡| av福利片在线| 久久久久久九九精品二区国产 | 亚洲av五月六月丁香网| 夜夜爽天天搞| 国产黄片美女视频| 中国美女看黄片| 国产精品98久久久久久宅男小说| 国产精品久久久久久人妻精品电影| 国产免费男女视频| 天天添夜夜摸| 国产伦在线观看视频一区| 国产一区二区三区在线臀色熟女| 99国产精品99久久久久| 国产一区二区在线av高清观看| 久久午夜综合久久蜜桃| 国产精品一及| 最近最新中文字幕大全电影3| 国产蜜桃级精品一区二区三区| 色尼玛亚洲综合影院| 亚洲片人在线观看| 在线播放国产精品三级| 国产精品,欧美在线| 级片在线观看| 宅男免费午夜| 欧美一级a爱片免费观看看 | 成人18禁高潮啪啪吃奶动态图| 亚洲,欧美精品.| 黄色丝袜av网址大全| 国产主播在线观看一区二区| 欧美国产日韩亚洲一区| avwww免费| 狂野欧美激情性xxxx| 五月伊人婷婷丁香| 国产成人啪精品午夜网站| 麻豆国产97在线/欧美 | 国产午夜精品久久久久久| www日本在线高清视频| 久9热在线精品视频| 亚洲国产精品999在线| 夜夜躁狠狠躁天天躁| 少妇裸体淫交视频免费看高清 | 亚洲一码二码三码区别大吗| 一区二区三区高清视频在线| 精品日产1卡2卡| 国产黄a三级三级三级人| 亚洲国产精品合色在线| 欧美日韩福利视频一区二区| 久久这里只有精品中国| 精品国产乱子伦一区二区三区| 19禁男女啪啪无遮挡网站| 欧美精品啪啪一区二区三区| 别揉我奶头~嗯~啊~动态视频| 亚洲乱码一区二区免费版| 日本一本二区三区精品| 日本a在线网址| 日韩大尺度精品在线看网址| 黄频高清免费视频| 12—13女人毛片做爰片一| 美女扒开内裤让男人捅视频| 久久久久久大精品| 午夜两性在线视频| www国产在线视频色| 亚洲熟妇熟女久久| 婷婷亚洲欧美| 精品无人区乱码1区二区| 免费电影在线观看免费观看| 国产精品亚洲美女久久久| 午夜精品在线福利| av在线播放免费不卡| 一级毛片高清免费大全| 日韩高清综合在线| 亚洲国产精品999在线| 国产探花在线观看一区二区| 欧美日韩乱码在线| 国产精品98久久久久久宅男小说| 桃红色精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 18禁黄网站禁片午夜丰满| 亚洲国产精品999在线| 亚洲国产中文字幕在线视频| 国产精品一及| 校园春色视频在线观看| av天堂在线播放| e午夜精品久久久久久久| 国产男靠女视频免费网站| 免费在线观看黄色视频的| 色哟哟哟哟哟哟| 精品一区二区三区四区五区乱码| 成人三级黄色视频| 黄色丝袜av网址大全| 五月玫瑰六月丁香| 波多野结衣巨乳人妻| 黄频高清免费视频| 精品不卡国产一区二区三区| 99精品久久久久人妻精品| 久久久国产欧美日韩av| 精品久久久久久成人av| 亚洲男人天堂网一区| 亚洲国产欧美网| cao死你这个sao货| 欧美日韩黄片免| 亚洲成a人片在线一区二区| av欧美777| 精品一区二区三区av网在线观看| 美女高潮喷水抽搐中文字幕| 在线国产一区二区在线| 久久精品国产亚洲av高清一级| 熟女电影av网| 免费在线观看亚洲国产| 日本精品一区二区三区蜜桃| 国产乱人伦免费视频| 首页视频小说图片口味搜索| 九色成人免费人妻av| 国产精品亚洲一级av第二区| 亚洲欧美精品综合久久99| a级毛片a级免费在线| 国产99白浆流出| 日本三级黄在线观看| 制服丝袜大香蕉在线| 久久久久久九九精品二区国产 |