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

    基于PCSWMM模型的城市內(nèi)澇風險評估

    2016-10-19 10:09:42吳海春黃國如
    水資源保護 2016年5期

    吳海春,黃國如, 2

    (1.華南理工大學土木與交通學院,廣東 廣州 510640;2.華南理工大學亞熱帶建筑科學國家重點實驗室,廣東 廣州 510640)

    ?

    基于PCSWMM模型的城市內(nèi)澇風險評估

    吳海春1,黃國如1, 2

    (1.華南理工大學土木與交通學院,廣東 廣州510640;2.華南理工大學亞熱帶建筑科學國家重點實驗室,廣東 廣州510640)

    構(gòu)建了基于PCSWMM的海甸島城市雨洪模型,采用3場次實測暴雨(20081013,20101005和20111005)內(nèi)澇淹沒資料對模型進行驗證,結(jié)果表明所構(gòu)建模型具有良好的精度和可靠性。以1 a、2 a、5 a、10 a和20 a 5種設(shè)計重現(xiàn)期降雨組合同頻率設(shè)計潮位進行模擬計算,對海甸島現(xiàn)狀排水能力進行評估,結(jié)果表明,海甸島84.8%管網(wǎng)排水能力低于1年一遇。以5 a、10 a、20 a、50 a降雨組合對應潮位對海甸島進行內(nèi)澇風險評估,得到了海甸島4種不同重現(xiàn)期暴雨的內(nèi)澇風險等級劃分狀況。關(guān)鍵詞:PCSWMM模型;管網(wǎng)排水能力;內(nèi)澇風險評估;海甸島

    近年來我國城市化快速發(fā)展,加之受全球氣候變化影響,暴雨等極端天氣事件頻發(fā),內(nèi)澇已成為阻礙城市可持續(xù)發(fā)展的重要因素。為了保障人民群眾的生命財產(chǎn)安全,提高城市的防災減災能力和安全保障水平,亟需開展城市內(nèi)澇風險評估工作?!冻鞘信潘?雨水)防澇綜合規(guī)劃編制大綱》(建城[2013]98號)和GB50014—2006《室外排水設(shè)計規(guī)范》(2014版)均推薦使用水力模型進行規(guī)劃模擬,以獲得雨水徑流的流態(tài)、水位變化、淹沒范圍、淹沒深度和淹沒時間等信息,然后采用單一指標或多個指標疊加,綜合評估城市內(nèi)澇災害的危險性[1]。

    在眾多的城市排水模型中,以美國環(huán)保局于20世紀70年代開發(fā)并持續(xù)改進的SWMM模型應用最為廣泛,但SWMM為一維水文水動力模型,無法計算溢出水量在地面的淹沒范圍和淹沒水深[2]。PCSWMM為加拿大水力計算研究所(Computational Hydraulics International,CHI)以SWMM為核心開發(fā)的水文水力學模型,已廣泛應用于一維管道與二維洪泛區(qū)耦合模擬、排水管網(wǎng)設(shè)計和評估、滯洪蓄水設(shè)計與評估、洪水風險分析等領(lǐng)域[3-6]。本文以??谑泻5閸u為研究對象,構(gòu)建PCSWMM城市雨洪模型,并對海甸島進行內(nèi)澇風險評估。

    1 PCSWMM模型基本原理

    1.1PCSWMM模型結(jié)構(gòu)

    PCSWMM是以SWMM為核心開發(fā)的城市雨洪模型,其主要計算原理、計算方法與SWMM基本一致,但相對SWMM而言,PCSWMM增強了前后處理能力和可視化等內(nèi)容[3-6]。

    PCSWMM模型概化思路為大小雙層排水系統(tǒng),即地表二維排水系統(tǒng)(或稱大排水系統(tǒng))和地下一維排水系統(tǒng)(或稱小排水系統(tǒng))。傳統(tǒng)城市雨洪模型只有一套排水系統(tǒng),該系統(tǒng)可以看成是由輸水管道、具有蓄水和銜接功能的節(jié)點以及出口組成,其主要核心是將城市地面概化成一個個“水庫”,“水庫”與“水庫”僅能通過排水管網(wǎng)進行水量交換[7]。后來,Djordjevic等[8]認為道路也具有行洪作用,于是在傳統(tǒng)城市雨洪模型中加入道路排水體系,稱為“只考慮了道路”的雙層排水系統(tǒng)。但“只考慮了道路”的雙層排水系統(tǒng)忽略了雨水在道路之外的其他地表二維的淹沒過程,精度達不到理想的地表二維模擬狀況。大小雙層排水系統(tǒng)將地面概化為“大排水系統(tǒng)”,除了在道路上方建立排水系統(tǒng),在非道路部分也建立一套排水系統(tǒng),另外考慮到建筑物等對水流的阻擋作用,在建筑物上不建立排水系統(tǒng)[3-6]。大小雙層排水系統(tǒng)示意圖見圖1。

    圖1 PCSWMM大小雙層排水系統(tǒng)示意圖

    1.2PCSWMM建模步驟

    a. 繪制阻礙層。創(chuàng)建二維網(wǎng)格前需繪制阻礙層,以考慮建筑物等對水流的阻擋作用,即在遙感影像圖上將這些區(qū)域圈出,不在此區(qū)域內(nèi)創(chuàng)建網(wǎng)格。

    b. 繪制邊界層。在城市區(qū)域,地表類型主要包括道路、草地和河道等。由于道路和草地的糙率不同,河道、道路等對水流的引導作用不同,因此,需繪制不同類型地表邊界,構(gòu)成3種地表類型邊界層。

    c. 基本設(shè)置。①網(wǎng)格類型:PCSWMM網(wǎng)格類型主要有六角形、定向、矩形和自適應等4種類型,其中定向網(wǎng)格主要用于窄長類型邊界,如河流和道路等,而且對定向網(wǎng)格還需繪制該邊界中心線,用于確定導流方向;自適應類型網(wǎng)格主要用于自行添加二維節(jié)點,或直接將CAD提取出的高程點作為二維節(jié)點類型;其他邊界可根據(jù)使用情況選擇六角形或矩形。②采樣因子m:采樣因子m用于生成標高點,每個網(wǎng)格內(nèi)生成m個標高點,取m個標高點的高程平均值作為網(wǎng)格和二維檢查井的井底標高。③網(wǎng)格分辨率、糙率等其他參數(shù)。

    d. 創(chuàng)建二維檢查井。根據(jù)DEM生成二維節(jié)點,在每個二維節(jié)點位置生成一個二維檢查井。為保證二維檢查井不發(fā)生溢流,將二維檢查井深度設(shè)為30 m,取二維檢查井水深作為該網(wǎng)格內(nèi)地表淹沒水深。PCSWMM計算方法與SWMM基本一致,當某檢查井水深大于井深時,多余水量就從檢查井溢出損失,或作為積水儲存于檢查井上方,但此時多余水量只是儲存于該檢查井上方,并未造成該檢查井的水深升高,會使計算的地表淹沒水深偏低。將二維檢查井深度設(shè)定為30 m,僅是為了保證即使發(fā)生內(nèi)澇,淹沒深度也不至于達到30 m,水流也不會從二維檢查井溢出,當然也可以將其設(shè)置為40 m、50 m等足夠大的數(shù)值。

    e. 創(chuàng)建網(wǎng)格及地表二維明渠。相鄰二維檢查井之間以30 m深的明渠管道連接,作為地表二維管道,利用泰森多邊形法繪制明渠管道的垂直平分線,每個檢查井周圍明渠管道的垂直平分線相交即組成了該網(wǎng)格的邊界,并將網(wǎng)格邊界長度作為被該邊界垂直平分的明渠寬度,網(wǎng)格內(nèi)所有二維明渠面積總和作為該網(wǎng)格面積。

    f. 確定網(wǎng)格和二維檢查井底標高。以二維節(jié)點或標高點的地表高程作為該網(wǎng)格和二維檢查井的井底高程。如果選擇以標高點的地表高程作為該網(wǎng)格和二維檢查井的井底高程,還需根據(jù)DEM生成標高點。取該網(wǎng)格內(nèi)若干標高點的標高平均值作為該網(wǎng)格和二維檢查井的井底高程[3-6]。

    g. 一維二維耦合連接。PCSWMM中一維二維耦合有兩種連接方式,分別為“使用底部孔口”和“直接連接到一維檢查井”,其中前者在網(wǎng)格內(nèi)唯一的二維檢查井和網(wǎng)格內(nèi)需要連接的一維檢查井之間建立一個以一維檢查井為起點、以二維檢查井為終點的底部孔口,并以二維檢查井與一維檢查井的下沿標高的高程差作為孔口的入口偏移量,這樣孔口剛好能連接至地面標高,適用于非河網(wǎng)地區(qū);后者將二維檢查井直接移至需要耦合的一維檢查井,并將原一維檢查井的上沿標高增加30 m,下沿標高不變,適用于河流或湖泊的一維二維耦合[3-6]。

    h. 設(shè)置一維檢查井超載深度。使用底部孔口連接一維檢查井和二維檢查井時,還需給一維檢查井設(shè)置不會發(fā)生溢流的超載深度,如20 m等。如前所述,當一維檢查井水深大于井深時,即水剛好到達孔口的入口處時,多余水量從該檢查井溢出損失,或作為積水儲存于該檢查井上方,使該一維檢查井水頭不再升高,水也就不能從一維檢查井通過孔口進入二維檢查井。而超載深度是指水從檢查井溢出前,檢查井內(nèi)水深可以達到的超過檢查井深度的最大深度。設(shè)置超載深度后,一維檢查井便不會發(fā)生溢流,也就能從一維檢查井通過孔口進入二維檢查井[3-6]。

    經(jīng)過以上步驟便可創(chuàng)建二維網(wǎng)格,并將一維、二維耦合連接,具體流程見圖2。網(wǎng)格的計算結(jié)果數(shù)據(jù)如網(wǎng)格最大淹沒深度等,根據(jù)網(wǎng)格內(nèi)二維檢查井及二維明渠等相關(guān)屬性計算得出[3-6]。

    圖2 創(chuàng)建二維網(wǎng)格及一維二維耦合流程

    2 PCSWMM模型構(gòu)建

    2.1研究區(qū)域概況

    海口市地處熱帶濱海地區(qū),海甸島位于??谑斜辈俊⒛隙山龊:涌诙?環(huán)抱于南渡江入??诜帚獾臋M溝河和海甸溪之中,面積13.8 km2。海甸島片區(qū)有海甸五西路明渠、鴨尾溪和白沙河等河道,海甸五西路明渠下游和鴨尾溪下游入海口處分別有一座閘門,對河道防洪排澇起重要作用,水閘凈寬分別為18.5 m和10.0 m,最大開度均為3.0 m。??谑卸嗄昶骄邓? 827 mm,其中5—10月為雨季,降水量占全年降水量的78.1%。海甸島四周臨海,多處路段地勢低洼,島上雨水通過地下排水管道排入河道,最終排入大海。大部分排水管網(wǎng)設(shè)計標準為一年一遇或以下,加之海潮頂托嚴重,每當遇到臺風暴雨或天文風暴潮時,許多路段嚴重積水[9]。

    2004—2012年,造成??谑袨暮Φ臒釒庑?4個,其中2008年10月13日、2010年10月5日、2011年10月5日等強降雨對海口市主城區(qū)造成了嚴重影響,尤其是2011年10月5日8時至10月6日8時,??谑斜辈垦睾5貐^(qū)為暴雨中心區(qū),最大降雨量88.1 mm/h,最大24 h降雨量441.0 mm,強降雨形成內(nèi)澇,海甸島內(nèi)的海甸五西路、人民大道、和平大道、海達路、海甸二東路、海甸三西路等路段嚴重積水,平均淹沒水深0.5 m,造成部分道路交通中斷,經(jīng)濟損失較為慘重[9]。

    2.2模型構(gòu)建

    2.2.1一維模型構(gòu)建

    在構(gòu)建PCSWMM模型前,需對研究區(qū)域排水系統(tǒng)進行概化,主要包括排水管網(wǎng)概化、河道概化。SWMM模型的水力要素包括節(jié)點和管段2種,節(jié)點包括窨井、雨水篦子、探測點、轉(zhuǎn)折點、出水口和蓄水設(shè)施;管段包括排水管道和溝渠。概化較大區(qū)域的排水管網(wǎng)時可把雨水篦子及其相連管道刪除。地表徑流最終通過排水管道匯集到河道中,河道蓄量、水位高低與管道排水能力關(guān)系密切。本模型將河道作為排水系統(tǒng)的一部分,將其分段概化成參數(shù)(河寬、河深、河底高程)各異的明渠。海甸島主要河道沿岸均有排水管道出水口。經(jīng)概化處理后,海甸島概化為3 510根管線、2 675個檢查井、60個出水口和8個蓄水設(shè)施。

    根據(jù)海甸島DEM、街區(qū)和道路分布圖,在ArcGIS中對海甸島進行劃分、調(diào)整和合并,再根據(jù)道路、管線及建筑物分布等對其進行子匯水區(qū)劃分。由于海甸島地勢整體較為平坦,在劃分過程中更加重視街道和社區(qū)單元的分布情況,采用泰森多邊形法將流域劃分為2 925個子匯水區(qū),并為各子匯水區(qū)指定流域出口,該出口可為排水管網(wǎng)的檢查井,也可設(shè)為下游的子流域子匯水區(qū)。

    坡度值為子匯水區(qū)的敏感參數(shù)之一。通過ArcGIS柵格表面坡度計算工具,對海甸島DEM進行坡度計算,獲取各子匯水區(qū)的平均坡度信息。地表漫流寬度直接對子匯水區(qū)匯流時間產(chǎn)生影響。依據(jù)SWMM模型手冊推薦公式,利用子匯水區(qū)面積與地表路徑長度的比值來計算地表漫流寬度。利用地形圖結(jié)合遙感影像圖提取各種用地類型(道路及廣場用地、屋面、城市綠化帶及公園)信息,統(tǒng)計各子匯水區(qū)的不透水率。

    SWMM模型中的水文水動力參數(shù)有些可根據(jù)研究區(qū)域的實際情況予以事先確定,如流域面積、坡度、匯水區(qū)不透水率等,其他一些參數(shù)如匯水區(qū)漫流寬度、不透水區(qū)曼寧系數(shù)、霍頓產(chǎn)流參數(shù)等,則根據(jù)相關(guān)文獻資料選定。本研究參考模型用戶手冊和相關(guān)研究成果取初值,輸入實測降雨和潮位數(shù)據(jù)進行計算,然后根據(jù)計算結(jié)果與實地調(diào)查情況進行對比,再調(diào)整參數(shù)重新計算。調(diào)整參數(shù)至與實地調(diào)查情況接近至誤差允許范圍內(nèi),得到的最終參數(shù)見表1,參數(shù)率定結(jié)果見2.3節(jié)。

    表1 子匯水區(qū)水動力參數(shù)取值

    2.2.2二維模型構(gòu)建

    按前述分析計算步驟構(gòu)建一維二維耦合模型,繪制研究區(qū)邊界、道路邊界、道路中心線圖層、阻礙圖層。其中河道部分在一維模型構(gòu)建時已概化為明渠,故在河道部分不再創(chuàng)建網(wǎng)格,將河道圖層也繪制在阻礙圖層。道路內(nèi)網(wǎng)格類型采用定向網(wǎng)格,分辨率取8 m,采樣因子m取3,曼寧系數(shù)取0.015。道路之外的其他研究區(qū)域網(wǎng)格類型采用六角形網(wǎng)格,分辨率取30 m,采樣因子m取3,曼寧系數(shù)取0.013。最終得到21 732個二維網(wǎng)格,面積從25 m2到1 000 m2不等,平均面積425 m2。鑒于本文所創(chuàng)建網(wǎng)格均處于非河網(wǎng)地區(qū),故選取“使用底部孔口”連接方式進行一維二維耦合。

    2.3模型驗證

    選取2011年10月5—6日暴雨進行參數(shù)率定,結(jié)果見表2。根據(jù)內(nèi)澇實際調(diào)研情況,該場次暴雨造成積水最為嚴重,提取模型計算淹沒點結(jié)果,將其與該場次暴雨實地調(diào)查澇點分布比較,模擬澇點結(jié)果與實際內(nèi)澇發(fā)生位置較為吻合,表明該模型能較為準確地模擬海甸島片區(qū)排水狀況。由表2可知,實測最大水深和模擬最大水深總體誤差較小,模擬得到的積水結(jié)果較為準確地反映了實際內(nèi)澇狀況。分別采用2008年10月13—14日和2010年10月5—6日場次降雨對所構(gòu)建的雨洪模型進行驗證,計算結(jié)果見表2,但由于2008年和2010年暴雨沒有實測淹沒資料,故表2中這兩場次暴雨只有模擬結(jié)果。由表2可看出,3場實測暴雨條件下,易澇點遭受不同程度淹沒,不同場次暴雨造成積水路段有所不同,其中積水最為嚴重的地方均發(fā)生在海甸五西路,與實際情況相同。

    表2 實測暴雨下各澇點最大淹沒深度 m

    3 排水系統(tǒng)能力評估

    3.1評估方法及邊界條件

    利用所構(gòu)建的PCSWMM模型評估海甸島現(xiàn)狀排水管網(wǎng)的排水能力,統(tǒng)計各種降雨條件下一維檢查井和一維管道超載情況。

    3.1.1設(shè)計暴雨

    海南省氣候中心與??谑袣庀缶钟?013年12月修改了??谑斜┯陱姸裙?該公式基于??谑袊覛庀笳?961—2012 年共計52年的自記雨量資料,依據(jù)皮爾遜Ⅲ型分布的最小二乘法修訂而成。設(shè)計暴雨強度計算公式為

    (1)

    式中:q為設(shè)計暴雨強度,L/(s·hm2);t為降雨歷時,min;P為設(shè)計重現(xiàn)期,a。

    采用國際上應用廣泛的芝加哥雨型進行降雨過程分配,降雨歷時為2 h,時間間隔為5 min,雨峰位置r=0.415,得到設(shè)計重現(xiàn)期為1 a、2 a、5 a、10 a和20 a的降雨過程。

    3.1.2設(shè)計潮位分析

    ??谑袑贋I海城市,排水易受潮位頂托影響,故需考慮潮位對排水的作用。近些年有些學者采用Copula函數(shù)建立區(qū)間暴雨和外江洪(潮)水位的聯(lián)合分布,利用聯(lián)合概率密度分析兩者的遭遇幾率,得出“對于不超過10年一遇的暴雨,遭遇同頻率的外江水位的幾率最大;但對超過10年一遇的暴雨,遭遇比其重現(xiàn)期小些的外江水位的幾率大些”[10-11],故對于重現(xiàn)期不超過10年一遇的暴雨采用同頻率設(shè)計潮位進行計算,但對于重現(xiàn)期超過10年一遇的暴雨,考慮其不利情況,本文亦采用同頻率設(shè)計潮位進行計算。

    表3 實測及設(shè)計暴雨條件下一維檢查井溢流和一維管道超載情況統(tǒng)計結(jié)果

    注:一維檢查井共2 675個,一維管道長共87 km。

    3.2排水系統(tǒng)能力評估結(jié)果

    利用PCSWMM模型得到實測及設(shè)計暴雨條件下一維檢查井溢流和一維管道超載情況統(tǒng)計結(jié)果(表3)。由表3可知,針對設(shè)計暴雨,隨著設(shè)計重現(xiàn)期(P)及組合潮位增大,溢流檢查井數(shù)量逐漸增多,超載管道長度逐漸增加,且可得出海甸島84.8%的排水管網(wǎng)排水能力低于1年一遇的結(jié)論。對于實測降雨而言,20101005和20111005場次降雨的降雨歷時和潮位均比較接近,20081013場次降雨的降雨歷時、降雨量及潮位均比20101005場次降雨高,但低于20111005場次的降雨量及雨強。由表3可知,20081013場次降雨的溢流檢查井個數(shù)和管道超載長度比20101005場次大,比20111005場次降雨小??傮w來說,降雨強度越大、潮位越高、歷時越長,溢流檢查井個數(shù)及超載管道長度越大。

    4 內(nèi)澇風險評估

    以??谟炅空?970—2012年43年實測暴雨序列為基礎(chǔ)資料,用皮爾遜Ⅲ型曲線法求得海口24 h不同重現(xiàn)期(5 a、10 a、20 a、50 a)的設(shè)計暴雨量值。選取1996年9月19日17:00至9月20日17:00的24 h暴雨過程為24 h典型暴雨過程(雨型),該場暴雨雨量大(398.7 mm),峰值也大(77.3 mm/h),且雨峰出現(xiàn)在暴雨過程后期,其形成的洪水主峰較大且靠后,對排水防澇安全較為不利。依據(jù)典型暴雨過程,以設(shè)計暴雨量值為控制,采取同頻率放大法推求24 h雨量時程分配。分別以5年一遇、10年一遇、20年一遇、50年一遇24 h設(shè)計降雨組合同頻率潮位進行動態(tài)模擬,對海甸島進行內(nèi)澇風險評估。

    評估標準采用國外常用的洪水風險評價標準,主要考慮洪水淹沒深度和流速來進行風險程度評估,具體計算公式如下:

    (2)

    式中:R為風險等級;d為積水深度,m;v為流速,m/s;n為常數(shù),取0.5;f為危害系數(shù)。

    當積水深度為0~0.25 m時,f取0;當積水深度為0.25~0.75 m時,草地/耕地、森林、城市的f值分別取0、0.5、1;當d>0.75 m或v>2 m/s時,草地/耕地、森林、城市的f值分別取0.5、1、1。根據(jù)計算出的內(nèi)澇風險指數(shù),風險等級劃分如下:R<0.75,風險極低;R值在0.75~1.25之間,風險低;R值在1.25~1.5之間,風險中等;R>1.5,風險高[12]。根據(jù)評估結(jié)果,分別列出不同重現(xiàn)期降雨情景下不同風險等級的面積(表4)。

    由表4可知,針對設(shè)計暴雨,隨著降雨強度增加,風險極低的面積幾乎不變,風險低的面積逐漸減小,風險中等和風險高的面積逐漸增大。針對實測暴雨,20081013場次降雨歷時較長,潮位較高,風險低、風險中等及風險高的面積均比其余兩場大,20111005場次降雨強度遠大于20101005場次降雨,風險低、中等及高風險面積均大于20101005場次降雨。由此可以得出結(jié)論,針對濱海城市,降雨歷時越長,潮位越高;降雨強度越大,內(nèi)澇越嚴重。

    表4 內(nèi)澇風險評估結(jié)果

    5 結(jié) 論

    a. 闡述了PCSWMM模型的基本原理及建模步驟,構(gòu)建了海甸島片區(qū)PCSWMM模型,利用20081013、20101005、20111005這3場實測暴雨及潮位過程計算分析,發(fā)現(xiàn)所構(gòu)建模型具有良好的精度和可靠性。

    b. 以1 a、2 a、5 a、10 a和20 a 5種不同設(shè)計重現(xiàn)期降雨及3場實測降雨組合對應潮位對海甸島進行排水現(xiàn)狀能力評估,結(jié)果表明84.8%管網(wǎng)排水能力僅處于1年一遇以下。針對濱海城市,降雨強度、潮位、降雨歷時越大,檢查井溢流和一維管網(wǎng)超載情況越嚴重。

    c. 以5年一遇、10年一遇、20年一遇和50年一遇24h 4種不同設(shè)計重現(xiàn)期降雨及3場實測降雨組合對應潮位對海甸島進行內(nèi)澇風險評估,結(jié)果表明在50年一遇降雨情景下8%的面積處于風險低、11.1%的面積處于風險中等、6.2%的面積處于風險高。針對濱海城市而言,降雨歷時越長,潮位越高,降雨強度越大,處于內(nèi)澇風險高的面積越大,內(nèi)澇越嚴重。

    [ 1 ] 林琳.合肥市典型沖溝雨水系統(tǒng)排水防澇能力評估研究[J].中國給水排水,2014,30(15):150-154.(LIN Lin.Assessment of flood drainage capacity of typical stormwater gully system in Hefei City [J].China Water & Wastewater,2014,30 (15): 150-154.(in Chinese))

    [ 2 ] 陳小龍,趙冬泉,盛政,等.Digital Water在城市排水防澇規(guī)劃中的應用[J].中國給水排水,2015,31(21):105-108.(CHEN Xiaolong,ZHAO Dongquan,SHENG Zheng,et al.Application of digital water simulation system in urban flood control planning [J].China Water & Wastewater,2015,31 (21): 105-108.(in Chinese))

    [ 3 ] JAMES W,ROSSMAN L A,JAMES W R C.User’s guide to SWMM 5:based on original USEPA SWMM documentation [M].Guelph,Ontario,Canada: Computational Hydraulics International,2010.

    [ 4 ] AHIABLAME L,SHAKYA R.Modeling flood reduction effects of low impact development at a watershed scale[J].Journal of Environmental Management,2016,171:81-91.

    [ 5 ] BECK J.Comparison of three methodologies for quasi-2d river flood modeling with SWMM5[J/OL].Journal of Water Management Modeling,2016,C402.https://www.chijournal.org/Content/Files/C402.pdf

    [ 6 ] KLENZENDORF B,PORESKY A,KELLY M,et al.1D/2D modeling of decentralized stormwater control measures for flood mitigation in Austin,Texas[C]//Proceedings of World Environmental and Water Resources Congress 2015.Alexander,American: American Society of Civil Engineers,2015:204-213.

    [ 7 ] 黃國如,馮杰,劉寧寧,等.城市雨洪模型及應用[M].北京:中國水利水電出版社,2013.

    [ 8 ] DJORDJEVIC S,PRODANOVIC D,MAKSIMOVIC.An approach to simulation of dual drainage[J].Water Science and Technology,1999,39(9): 95-103.

    [ 9 ] 黃國如,黃維,張靈敏,等.基于GIS和SWMM模型的城市暴雨積水模擬[J].水資源與水工程學報,2015,26(4):1-6.(HUANG Guoru,HUANG Wei,ZHANG Lingmin,et al.Simulation of rainstorm waterlogging in urban areas based on GIS and SWMM model [J].Journal of Water Resources & Water Engineering,2015,26 (4): 1-6.(in Chinese))

    [10] 武傳號,黃國如,吳思遠.基于Copula函數(shù)的廣州市短歷時暴雨與潮位組合風險分析[J].水力發(fā)電學報,2014,33(2):33-40.(WU Chuanhao,HUANG Guoru,WU Siyuan.Risk analysis of combinations of short duration rainstorm and tidal level in Guangzhou based on Copula function [J].Journal of Hydroelectric Engineering,2014,33 (2): 33-40.(in Chinese))

    [11] 劉曾美,陳子燊.區(qū)間暴雨和外江洪水位遭遇組合的風險[J].水科學進展,2009,20(5):619-625.(LIU Zengmei,CHEN Zisheng.Risk study of the bivariate encounter of interzone rainstorm and flood level of the outer river [J].Advance in Water Science,2009,20 (5): 619-625.(in Chinese))

    [12] WADE S,RAMSBOTTOM D,FLOYD P,et al.Risks to people: developing new approaches for flood hazard and vulnerability mapping[C]//Proceedings of Defra Flood & Coastal Management Conference.York,UK: HR Wallingford, 2005.

    Risk assessment of urban waterlogging based on PCSWMM model

    WU Haichun1, HUANG Guoru1, 2

    (1.SchoolofCivilEngineeringandTransportation,SouthChinaUniversityofTechnology,Guangzhou510640,China;2.StateKeyLaboratoryofSubtropicalBuildingScience,SouthChinaUniversityofTechnology,Guangzhou510640,China)

    An urban rainstorm model based on PCSWMM for Haidian Island in Haikou City was developed. The model was verified by three measured rainfalls events and the corresponding tidal graph records in 20081013, 20101005, and 20111005. The simulated waterlogging submerged positions and depths are basically consistent with the actual investigation situation, which shows that the model has high precision and reliability. It respectively evaluated the present situation of drainage networks and the waterlogging risks in Haidian Island in Haikou City by using design rainstorms of return periods of 1, 2, 5, 10, 20, and 50 years with the corresponding tidal levels. The result shows that the present situation of 84.8% of the drainage networks in Haidian Island in Haikou City is less than 1-year design return period. Meanwhile, the waterlogging risk assessment of design rainstorms with 5, 10, 20, 50-year return periods and three measured rainstorms with the corresponding tidal levels are analyzed.

    PCSWMM model; water discharge capacity of drainage networks; waterlogging risk assessment; Haidian Island

    10.3880/j.issn.1004-6933.2016.05.003

    水利部公益性行業(yè)科研專項(201301093, 201401048);廣東省科技計劃項目(2016A020223003);華南理工大學亞熱帶建筑科學國家重點實驗室自主研究課題項目(2014ZC09);廣州市水務(wù)科技項目(GZSW-201401)

    吳海春(1988—),男,碩士研究生,研究方向為水文水資源。E-mail: whcscut@163.com

    黃國如,教授,博士生導師。E-mail: huanggr@scut.edu.cn

    TU992

    A

    1004-6933(2016)05-0011-06

    2016-06-22編輯:彭桃英)

    亚洲五月婷婷丁香| 高清视频免费观看一区二区| 日韩欧美免费精品| 日韩制服骚丝袜av| 人人妻人人添人人爽欧美一区卜| 久热爱精品视频在线9| 午夜精品国产一区二区电影| a级片在线免费高清观看视频| 久久人妻熟女aⅴ| 亚洲人成电影免费在线| www.精华液| 亚洲中文日韩欧美视频| 999久久久国产精品视频| 亚洲成人国产一区在线观看| 久久精品成人免费网站| svipshipincom国产片| 免费一级毛片在线播放高清视频 | 精品人妻一区二区三区麻豆| 日韩视频一区二区在线观看| 精品久久久久久电影网| 亚洲av成人一区二区三| 99国产极品粉嫩在线观看| 精品少妇黑人巨大在线播放| 亚洲五月婷婷丁香| 久久热在线av| 婷婷成人精品国产| 热99久久久久精品小说推荐| 久久久久国内视频| av天堂久久9| 丝袜脚勾引网站| 午夜福利视频在线观看免费| 色94色欧美一区二区| 成年人黄色毛片网站| 欧美少妇被猛烈插入视频| 五月天丁香电影| 国产在线观看jvid| 交换朋友夫妻互换小说| 中文精品一卡2卡3卡4更新| 咕卡用的链子| 在线观看一区二区三区激情| av网站免费在线观看视频| 亚洲国产欧美在线一区| 国产精品一区二区精品视频观看| 蜜桃在线观看..| 中国美女看黄片| 美女扒开内裤让男人捅视频| 97精品久久久久久久久久精品| 日韩电影二区| 久久久久国内视频| 我的亚洲天堂| 嫩草影视91久久| 女性生殖器流出的白浆| 日韩一卡2卡3卡4卡2021年| 久久 成人 亚洲| 成年av动漫网址| av有码第一页| 日本撒尿小便嘘嘘汇集6| 国产成人精品久久二区二区免费| 国产日韩欧美亚洲二区| 免费女性裸体啪啪无遮挡网站| e午夜精品久久久久久久| 精品国产国语对白av| 男人舔女人的私密视频| 丝袜美腿诱惑在线| 国产97色在线日韩免费| 黄色视频在线播放观看不卡| 亚洲av美国av| 国产极品粉嫩免费观看在线| 无限看片的www在线观看| 丝袜喷水一区| 电影成人av| 国产高清videossex| 秋霞在线观看毛片| 亚洲中文av在线| 国产极品粉嫩免费观看在线| 一本大道久久a久久精品| 波多野结衣av一区二区av| 肉色欧美久久久久久久蜜桃| 夫妻午夜视频| 天堂俺去俺来也www色官网| 亚洲男人天堂网一区| 久久九九热精品免费| 亚洲情色 制服丝袜| 久久久精品区二区三区| 久久久水蜜桃国产精品网| 国产熟女午夜一区二区三区| 搡老熟女国产l中国老女人| 国产日韩欧美亚洲二区| 久久久久久久久久久久大奶| 高清在线国产一区| 中文字幕制服av| 女性生殖器流出的白浆| 精品久久蜜臀av无| 亚洲熟女毛片儿| 深夜精品福利| 视频区图区小说| 亚洲中文日韩欧美视频| 日本欧美视频一区| 亚洲五月婷婷丁香| 99九九在线精品视频| 涩涩av久久男人的天堂| 亚洲一区二区三区欧美精品| 亚洲 欧美一区二区三区| 男人舔女人的私密视频| 久久精品国产亚洲av高清一级| 亚洲精品国产一区二区精华液| 一级黄色大片毛片| 亚洲国产日韩一区二区| 首页视频小说图片口味搜索| 国产精品久久久久久精品电影小说| 久久精品aⅴ一区二区三区四区| 亚洲精品日韩在线中文字幕| 看免费av毛片| 无遮挡黄片免费观看| 日本a在线网址| 国产成人啪精品午夜网站| 99热全是精品| 日韩视频一区二区在线观看| 男女国产视频网站| 丝瓜视频免费看黄片| 免费观看av网站的网址| 9色porny在线观看| 狠狠精品人妻久久久久久综合| 日韩,欧美,国产一区二区三区| 黄色片一级片一级黄色片| 精品久久久精品久久久| 国产不卡av网站在线观看| 亚洲熟女精品中文字幕| 我要看黄色一级片免费的| 国产淫语在线视频| 国产成人影院久久av| 久久99热这里只频精品6学生| 欧美+亚洲+日韩+国产| 国产福利在线免费观看视频| 日韩 亚洲 欧美在线| 亚洲一区二区三区欧美精品| 成人国产av品久久久| 熟女少妇亚洲综合色aaa.| 搡老熟女国产l中国老女人| 成人国产一区最新在线观看| 99re6热这里在线精品视频| 欧美精品一区二区免费开放| 可以免费在线观看a视频的电影网站| 啦啦啦视频在线资源免费观看| 一边摸一边做爽爽视频免费| 精品一区二区三卡| 天天添夜夜摸| 亚洲人成电影免费在线| 黄色a级毛片大全视频| 亚洲欧美精品综合一区二区三区| 搡老岳熟女国产| 成在线人永久免费视频| av福利片在线| 精品少妇久久久久久888优播| 国产精品一区二区精品视频观看| 国产在线观看jvid| 亚洲九九香蕉| 久久这里只有精品19| 91老司机精品| 午夜两性在线视频| 宅男免费午夜| 最新的欧美精品一区二区| 亚洲精品乱久久久久久| 国产成人av激情在线播放| 国产成+人综合+亚洲专区| 18禁观看日本| 国产精品香港三级国产av潘金莲| 亚洲九九香蕉| 亚洲av男天堂| 王馨瑶露胸无遮挡在线观看| 亚洲 国产 在线| 欧美激情极品国产一区二区三区| 十八禁人妻一区二区| √禁漫天堂资源中文www| 老司机亚洲免费影院| av天堂在线播放| videosex国产| 视频区图区小说| 啦啦啦在线免费观看视频4| 亚洲三区欧美一区| 亚洲视频免费观看视频| 最近最新免费中文字幕在线| 国产精品免费大片| 色老头精品视频在线观看| 美女午夜性视频免费| 51午夜福利影视在线观看| 18禁裸乳无遮挡动漫免费视频| 咕卡用的链子| 精品亚洲成国产av| 大陆偷拍与自拍| 亚洲国产av新网站| 亚洲成人手机| 成人国产一区最新在线观看| videosex国产| 久久这里只有精品19| 91精品伊人久久大香线蕉| 久久久久久久国产电影| 亚洲精品粉嫩美女一区| av福利片在线| 中文字幕高清在线视频| 自线自在国产av| 99精国产麻豆久久婷婷| 国产精品亚洲av一区麻豆| 丝袜美足系列| 亚洲伊人色综图| 丝袜美腿诱惑在线| 国产成人欧美在线观看 | 成年人黄色毛片网站| 国产成人av激情在线播放| 亚洲精品国产区一区二| 一二三四社区在线视频社区8| 成人手机av| 久久青草综合色| 成人av一区二区三区在线看 | 国产成人a∨麻豆精品| 91字幕亚洲| 久久精品人人爽人人爽视色| 国产一区二区在线观看av| 久久精品国产综合久久久| 国产精品影院久久| 国内毛片毛片毛片毛片毛片| 在线av久久热| 老司机影院成人| 免费少妇av软件| 69av精品久久久久久 | 热re99久久国产66热| 自线自在国产av| 精品国产一区二区三区久久久樱花| 国产成人免费观看mmmm| av不卡在线播放| 国产在线一区二区三区精| 啦啦啦 在线观看视频| 免费人妻精品一区二区三区视频| 亚洲精品国产色婷婷电影| 色老头精品视频在线观看| 成年动漫av网址| 中国美女看黄片| 两个人看的免费小视频| 秋霞在线观看毛片| 亚洲男人天堂网一区| tocl精华| 淫妇啪啪啪对白视频 | 欧美黄色片欧美黄色片| 如日韩欧美国产精品一区二区三区| 看免费av毛片| 在线永久观看黄色视频| 欧美日韩亚洲国产一区二区在线观看 | 美女脱内裤让男人舔精品视频| 欧美大码av| 国产精品麻豆人妻色哟哟久久| 欧美日韩视频精品一区| 亚洲欧美日韩另类电影网站| 俄罗斯特黄特色一大片| 脱女人内裤的视频| 免费在线观看完整版高清| 蜜桃在线观看..| 国产欧美日韩一区二区精品| 亚洲熟女毛片儿| 最近最新中文字幕大全免费视频| 咕卡用的链子| 亚洲人成77777在线视频| 黑人欧美特级aaaaaa片| kizo精华| 国产片内射在线| 日本91视频免费播放| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产超薄肉色丝袜足j| 国产免费福利视频在线观看| 亚洲专区字幕在线| 成在线人永久免费视频| 亚洲成人手机| 久久久国产精品麻豆| 精品第一国产精品| 他把我摸到了高潮在线观看 | 99国产精品一区二区三区| 亚洲中文av在线| 两个人看的免费小视频| 欧美黑人精品巨大| 成人三级做爰电影| 午夜精品国产一区二区电影| 免费高清在线观看日韩| 国产成人一区二区三区免费视频网站| 成人18禁高潮啪啪吃奶动态图| 亚洲国产欧美网| 久久热在线av| 欧美 亚洲 国产 日韩一| 少妇精品久久久久久久| 精品国产乱码久久久久久男人| 国产成人欧美在线观看 | 汤姆久久久久久久影院中文字幕| 啦啦啦中文免费视频观看日本| 国产在线一区二区三区精| 国产成人精品久久二区二区91| 欧美精品av麻豆av| 欧美亚洲日本最大视频资源| 国产成人精品在线电影| 大码成人一级视频| 侵犯人妻中文字幕一二三四区| 国产在线视频一区二区| 桃花免费在线播放| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲国产av新网站| 一级a爱视频在线免费观看| 国内毛片毛片毛片毛片毛片| 一区二区三区乱码不卡18| 19禁男女啪啪无遮挡网站| 久久影院123| 久久狼人影院| 亚洲精品乱久久久久久| 丰满饥渴人妻一区二区三| 十八禁网站免费在线| 亚洲五月婷婷丁香| 精品少妇内射三级| 汤姆久久久久久久影院中文字幕| 亚洲国产av影院在线观看| 99国产精品一区二区三区| 久久久欧美国产精品| 久久性视频一级片| 青青草视频在线视频观看| 男人舔女人的私密视频| 成年美女黄网站色视频大全免费| 亚洲欧洲精品一区二区精品久久久| 亚洲 国产 在线| 纵有疾风起免费观看全集完整版| 青青草视频在线视频观看| 日本五十路高清| 最近最新免费中文字幕在线| 久久 成人 亚洲| 男人添女人高潮全过程视频| 国产精品成人在线| 99久久99久久久精品蜜桃| 午夜免费观看性视频| 欧美日韩一级在线毛片| 欧美性长视频在线观看| 日本精品一区二区三区蜜桃| 首页视频小说图片口味搜索| 精品国产超薄肉色丝袜足j| 性色av乱码一区二区三区2| 国产精品 欧美亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成国产人片在线观看| 精品人妻1区二区| 久久久精品国产亚洲av高清涩受| 日韩制服丝袜自拍偷拍| bbb黄色大片| 999久久久国产精品视频| 99香蕉大伊视频| 欧美日韩黄片免| 两性夫妻黄色片| 妹子高潮喷水视频| 美女主播在线视频| 中文字幕制服av| 99精品久久久久人妻精品| 美女视频免费永久观看网站| 中文字幕人妻熟女乱码| 成人国产一区最新在线观看| 99热网站在线观看| 丝袜喷水一区| 亚洲avbb在线观看| 日韩欧美一区视频在线观看| 欧美少妇被猛烈插入视频| 国产av一区二区精品久久| 国产激情久久老熟女| avwww免费| 日本猛色少妇xxxxx猛交久久| avwww免费| 亚洲国产中文字幕在线视频| 少妇被粗大的猛进出69影院| 黑人欧美特级aaaaaa片| 在线天堂中文资源库| 久久青草综合色| 成年人免费黄色播放视频| 亚洲色图综合在线观看| 久久久国产精品麻豆| 丝袜美足系列| 肉色欧美久久久久久久蜜桃| 手机成人av网站| 电影成人av| 97精品久久久久久久久久精品| 久久精品成人免费网站| 欧美日韩黄片免| 麻豆国产av国片精品| 老司机靠b影院| 别揉我奶头~嗯~啊~动态视频 | 在线观看免费日韩欧美大片| 黄色 视频免费看| 国产欧美日韩一区二区精品| 嫩草影视91久久| 亚洲伊人久久精品综合| 免费少妇av软件| 高清av免费在线| 一区二区日韩欧美中文字幕| 精品第一国产精品| 亚洲精品国产精品久久久不卡| 真人做人爱边吃奶动态| 久久国产精品大桥未久av| 精品久久久精品久久久| 69av精品久久久久久 | 他把我摸到了高潮在线观看 | 欧美亚洲日本最大视频资源| 热99re8久久精品国产| 精品人妻1区二区| 亚洲精品乱久久久久久| 一区二区日韩欧美中文字幕| 亚洲情色 制服丝袜| 久久久久精品人妻al黑| 9191精品国产免费久久| 国产成人精品久久二区二区91| 97人妻天天添夜夜摸| 又大又爽又粗| 飞空精品影院首页| 亚洲伊人久久精品综合| 精品一区二区三卡| 91大片在线观看| 国产在线视频一区二区| 丰满少妇做爰视频| 欧美日韩亚洲综合一区二区三区_| 两性午夜刺激爽爽歪歪视频在线观看 | 91老司机精品| 麻豆乱淫一区二区| 亚洲一区二区三区欧美精品| 黄色 视频免费看| 国产成人一区二区三区免费视频网站| 国产一区有黄有色的免费视频| 爱豆传媒免费全集在线观看| 一区二区日韩欧美中文字幕| 淫妇啪啪啪对白视频 | 9色porny在线观看| 亚洲av日韩精品久久久久久密| 91九色精品人成在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲国产精品成人久久小说| 久久性视频一级片| 日韩欧美国产一区二区入口| 亚洲国产欧美在线一区| 亚洲精品一二三| 色视频在线一区二区三区| 亚洲精品粉嫩美女一区| 高清欧美精品videossex| 国产福利在线免费观看视频| 国产精品一二三区在线看| 久9热在线精品视频| 乱人伦中国视频| 亚洲国产欧美在线一区| 亚洲第一青青草原| 热re99久久国产66热| 男女之事视频高清在线观看| 亚洲精品成人av观看孕妇| 麻豆国产av国片精品| 12—13女人毛片做爰片一| 少妇人妻久久综合中文| 美女中出高潮动态图| 国产精品一区二区在线不卡| 性高湖久久久久久久久免费观看| 国产精品久久久久成人av| 亚洲av欧美aⅴ国产| 天堂俺去俺来也www色官网| 国产欧美日韩综合在线一区二区| 另类亚洲欧美激情| 日本av手机在线免费观看| 亚洲av日韩精品久久久久久密| 中文字幕精品免费在线观看视频| 1024视频免费在线观看| 国产精品久久久久久精品电影小说| 肉色欧美久久久久久久蜜桃| 天天躁日日躁夜夜躁夜夜| 宅男免费午夜| 国产成人精品在线电影| 国产日韩一区二区三区精品不卡| 欧美精品高潮呻吟av久久| 精品国产乱码久久久久久小说| 黄色 视频免费看| 51午夜福利影视在线观看| 精品高清国产在线一区| 在线永久观看黄色视频| 亚洲avbb在线观看| 亚洲精品国产av蜜桃| 日韩视频一区二区在线观看| 91大片在线观看| 国产亚洲精品第一综合不卡| 亚洲国产精品一区二区三区在线| 久久99热这里只频精品6学生| 人人妻人人爽人人添夜夜欢视频| 成年女人毛片免费观看观看9 | 亚洲欧美一区二区三区黑人| 天天躁日日躁夜夜躁夜夜| 国产亚洲欧美在线一区二区| 亚洲 欧美一区二区三区| 男女床上黄色一级片免费看| 丰满迷人的少妇在线观看| 国产欧美日韩一区二区三 | 少妇 在线观看| 成年人黄色毛片网站| 日韩欧美国产一区二区入口| 夫妻午夜视频| 成人免费观看视频高清| 9热在线视频观看99| 欧美日韩视频精品一区| 国产一区有黄有色的免费视频| av天堂久久9| 久久女婷五月综合色啪小说| 水蜜桃什么品种好| 亚洲人成电影免费在线| 日韩电影二区| 日日爽夜夜爽网站| 亚洲一码二码三码区别大吗| 飞空精品影院首页| 亚洲综合色网址| 久久av网站| 欧美乱码精品一区二区三区| 免费看十八禁软件| 午夜免费成人在线视频| 中文字幕最新亚洲高清| 国产av精品麻豆| 久久人人爽人人片av| 一区二区av电影网| 国产精品久久久人人做人人爽| 国精品久久久久久国模美| 满18在线观看网站| 悠悠久久av| 老司机福利观看| 国产精品二区激情视频| 欧美人与性动交α欧美软件| 青草久久国产| 午夜影院在线不卡| 高潮久久久久久久久久久不卡| 99国产综合亚洲精品| 一边摸一边做爽爽视频免费| 桃花免费在线播放| 亚洲男人天堂网一区| 国产成+人综合+亚洲专区| 在线观看免费日韩欧美大片| av在线老鸭窝| 我的亚洲天堂| 高清黄色对白视频在线免费看| 12—13女人毛片做爰片一| 成人三级做爰电影| 美女高潮到喷水免费观看| 欧美日韩精品网址| 国产麻豆69| 午夜福利一区二区在线看| 色播在线永久视频| 最新的欧美精品一区二区| 精品人妻一区二区三区麻豆| 欧美久久黑人一区二区| 香蕉丝袜av| 啦啦啦啦在线视频资源| 老司机午夜福利在线观看视频 | 天堂中文最新版在线下载| 亚洲精品国产区一区二| 成人三级做爰电影| 国产亚洲欧美在线一区二区| 久久久久国产精品人妻一区二区| 亚洲伊人久久精品综合| 午夜福利,免费看| 国产精品久久久久成人av| 欧美日韩视频精品一区| 日韩 欧美 亚洲 中文字幕| 国产av国产精品国产| 亚洲精品中文字幕一二三四区 | 久久久国产精品麻豆| 少妇 在线观看| 操美女的视频在线观看| 亚洲成av片中文字幕在线观看| 老司机影院成人| 国产极品粉嫩免费观看在线| 热99久久久久精品小说推荐| 亚洲专区字幕在线| 国产黄色免费在线视频| 老汉色∧v一级毛片| 亚洲一区二区三区欧美精品| 色婷婷av一区二区三区视频| 亚洲专区中文字幕在线| av片东京热男人的天堂| 在线观看人妻少妇| 久久久久精品国产欧美久久久 | 麻豆av在线久日| 少妇粗大呻吟视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美97在线视频| 色视频在线一区二区三区| 最近中文字幕2019免费版| 日韩 亚洲 欧美在线| 精品熟女少妇八av免费久了| 免费少妇av软件| 精品国产乱子伦一区二区三区 | 精品久久久精品久久久| a级毛片在线看网站| 久久久久久免费高清国产稀缺| 99久久精品国产亚洲精品| 国产在线观看jvid| 欧美大码av| 最近中文字幕2019免费版| 久久国产精品影院| 人人妻人人澡人人看| 嫩草影视91久久| 别揉我奶头~嗯~啊~动态视频 | 免费不卡黄色视频| h视频一区二区三区| av国产精品久久久久影院| 亚洲九九香蕉| 亚洲欧美日韩另类电影网站| 成人手机av| 亚洲av成人不卡在线观看播放网 | 人妻 亚洲 视频| 久久99热这里只频精品6学生| 狠狠婷婷综合久久久久久88av| 又黄又粗又硬又大视频| 少妇裸体淫交视频免费看高清 | 99热网站在线观看| 视频区欧美日本亚洲| 中文字幕精品免费在线观看视频| 亚洲精品自拍成人| 男人舔女人的私密视频| 777米奇影视久久|