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

    基于SWMM模型的城市暴雨洪水模擬與分析

    2017-11-01 01:24:46,趙,
    水利信息化 2017年5期
    關(guān)鍵詞:龍洞恩施市清江

    楊 斌 ,趙 睿 ,

    (1. 湖北工業(yè)大學(xué),湖北 武漢 430068;2. 長(zhǎng)江中游水文水資源勘測(cè)局,湖北 武漢 430010)

    基于SWMM模型的城市暴雨洪水模擬與分析

    楊 斌1,趙 睿1,2

    (1. 湖北工業(yè)大學(xué),湖北 武漢 430068;2. 長(zhǎng)江中游水文水資源勘測(cè)局,湖北 武漢 430010)

    為了進(jìn)行城市暴雨洪水模擬分析,以湖北省恩施市中心城區(qū)為研究區(qū)域,在河流水系分布和洪水來(lái)源分析的基礎(chǔ)上,提出洪水分析計(jì)算方案,構(gòu)建一、二維耦合 SWMM 模型計(jì)算典型頻率(10,20,50,100 年一遇)河道洪水演進(jìn)和漫堤淹沒(méi)過(guò)程,獲得較為準(zhǔn)確的洪水淹沒(méi)范圍、水深、淹沒(méi)歷時(shí)和洪峰到達(dá)時(shí)間等信息,并結(jié)合 GIS 進(jìn)行洪水影響評(píng)價(jià)。結(jié)果表明,清江干流主要河段基本滿足 20 年一遇的防洪標(biāo)準(zhǔn),清江干流和龍洞河部分河段現(xiàn)狀防洪能力較低。SWMM 模型是一種成熟的城市雨洪徑流模型,適宜于城市暴雨洪水模擬與分析應(yīng)用。

    暴雨洪水管理模型;一、二維耦合;暴雨洪水模擬;洪水影響分析

    0 引言

    伴隨著我國(guó)城市化進(jìn)程的不斷加快,城市人口急劇增加導(dǎo)致城市下墊面環(huán)境狀況發(fā)生顯著變化。在全球變化的背景下,城市化進(jìn)程直接影響了流域產(chǎn)匯流過(guò)程,以及區(qū)域水資源供、用、耗、排水條件,使得城市小尺度氣候發(fā)生變化。受到下墊面環(huán)境變化和氣候變化等多因素的共同作用,我國(guó)城市洪澇災(zāi)害發(fā)生次數(shù)不斷增加,危害程度不斷提升[1]。2012 年 7 月 21 日,北京市出現(xiàn)特大暴雨,全市受災(zāi)人口 160.2 萬(wàn)人,緊急轉(zhuǎn)移 9.7 萬(wàn)人,損毀公共設(shè)施448 處,因?yàn)?zāi)直接經(jīng)濟(jì)損失高達(dá) 116.4 億元。因此,在城市重點(diǎn)區(qū)域內(nèi),針對(duì)可能發(fā)生的暴雨洪水淹沒(méi)水深、淹沒(méi)范圍、洪峰到達(dá)時(shí)間等洪水演進(jìn)過(guò)程要素進(jìn)行模擬計(jì)算,并結(jié)合 GIS 進(jìn)行洪水影響分析,以標(biāo)示城市范圍內(nèi)受洪水災(zāi)害的危險(xiǎn)程度,對(duì)于城市防洪除澇減災(zāi)、人口轉(zhuǎn)移疏散、防洪預(yù)案編制等具有重要的意義[2]。

    為了解決城市內(nèi)澇和排水等問(wèn)題,世界各國(guó)研究了城市雨洪徑流模型,其中比較典型的是暴雨洪水管理模型(SWMM)、蓄水、處理和溢流模型(STORM)和沃林福德模型(Wallinford)[3-4]。SWMM 模型是在美國(guó)環(huán)境保護(hù)署資助下,由美國(guó)水資源有限公司、梅特卡夫有限公司和佛羅里達(dá)大學(xué)共同研發(fā)的,提供了地表產(chǎn)流、河道徑流、管網(wǎng)匯流和水污染負(fù)荷的模擬分析,其中地表徑流采用非線性水庫(kù)模型,管網(wǎng)匯流采用圣維南方程,水污染負(fù)荷采用累積—沖刷模型[5]。STORM 模型是由美國(guó)陸軍工程兵團(tuán)開(kāi)發(fā)的,主要用于降雨過(guò)程中水量、水質(zhì)的模擬和城市蓄水、地表漫流等計(jì)算[6]。Wallinford 模型是由英國(guó)環(huán)境部開(kāi)發(fā)的,可進(jìn)行城市水循環(huán)過(guò)程的模擬分析、城市雨水、污水排水系統(tǒng)能力模擬[7]。由于 SWMM 模型是提出最早、應(yīng)用最廣泛的城市雨洪徑流模型,本研究選用了該模型進(jìn)行城市暴雨洪水模擬計(jì)算。

    國(guó)內(nèi)外學(xué)術(shù)界對(duì) SWMM 模型參數(shù)敏感性分析、模型適用范圍、模型擴(kuò)展應(yīng)用等相關(guān)研究較為深入,特別是基于 GIS 的 SWMM 模型應(yīng)用成為一大熱點(diǎn)。Zaghloul 等[8]運(yùn)用神經(jīng)網(wǎng)絡(luò)算法對(duì)徑流參數(shù)、匯水區(qū)域參數(shù)等進(jìn)行分析,初步提出了 SWMM 模型參數(shù)適用范圍。Smith 等[9]將 GIS 和 SWMM 相結(jié)合,研究了土地利用類型對(duì)降雨徑流的影響。董欣等[10]在運(yùn)用不確定分析的 HSY 算法和 Monte Carlo采樣法,在實(shí)測(cè)數(shù)據(jù)的基礎(chǔ)上,對(duì) SWMM 模型參數(shù)進(jìn)行驗(yàn)證,為模型參數(shù)的確定提供依據(jù)。Krebs 等[11]利用基因遺傳算法研究了加入低影響開(kāi)發(fā)(LID)措施的現(xiàn)代城市區(qū)域,提出了加入 LID 措施的模型參數(shù)適用范圍。周力寧[12]提出暴雨內(nèi)澇風(fēng)險(xiǎn)識(shí)別技術(shù),進(jìn)行了不同暴雨強(qiáng)度下內(nèi)澇風(fēng)險(xiǎn)等級(jí)分析,并基于 SWMM 模型計(jì)算內(nèi)澇發(fā)生時(shí)積水區(qū)域的積水深度和內(nèi)澇風(fēng)險(xiǎn)等級(jí)。申邵洪等[13]選擇湖北省??悼h為研究區(qū)域,利用 TUFLOW 軟件構(gòu)建暴雨洪水計(jì)算和淹沒(méi)分析模型,對(duì)城區(qū)受淹地區(qū)、淹沒(méi)水深、淹沒(méi)面積等進(jìn)行了深入分析。

    恩施市位于高山河流沖刷形成的峽谷盆地中,中心城區(qū)地勢(shì)低,境內(nèi)河流眾多,水系形態(tài)復(fù)雜,涉河工程多,下游出口斷面窄,城市排水不暢,暴雨山洪引發(fā)的城市內(nèi)澇積水頻發(fā),在山區(qū)城市暴雨洪水分析方面具有一定的代表性。因此,本文選擇恩施市中心城區(qū)為研究區(qū)域,在河流水系分布和洪水來(lái)源分析的基礎(chǔ)上,提出洪水分析計(jì)算方案,構(gòu)建了一、二維耦合 SWMM 模型,計(jì)算典型頻率(10,20,50,100 年一遇)洪水演進(jìn)過(guò)程中淹沒(méi)范圍、水深、淹沒(méi)歷時(shí)和洪峰到達(dá)時(shí)間等信息,并結(jié)合 GIS 技術(shù)進(jìn)行洪水影響評(píng)價(jià),為城市暴雨洪水內(nèi)澇分析和防洪預(yù)案編制提供范例。鑒于恩施市排水管網(wǎng)資料匱乏,為了簡(jiǎn)化暴雨洪水過(guò)程分析,主要使用 SWMM 模型的地表產(chǎn)流、河道徑流模塊。

    1 研究區(qū)概況

    恩施市是恩施土家族苗族自治州政府所在地,位于湖北省西南部,清江中上游。全市總面積 3 967 km2,轄 3 個(gè)街道、4 個(gè)鎮(zhèn)、9 個(gè)鄉(xiāng),1 個(gè)辦事處,172 個(gè)行政村和 34 個(gè)社區(qū)。據(jù) 2013 年統(tǒng)計(jì),全市戶籍總?cè)丝?81.43 萬(wàn)人,其中城鎮(zhèn)人口(含鄉(xiāng)鎮(zhèn)居委會(huì))36.76 萬(wàn)人,鄉(xiāng)村人口 44.67 萬(wàn)人。

    恩施市城區(qū)主要河流是清江干流,發(fā)源于利川西北部齊躍山東麓的龍洞溝,自北向南從流經(jīng)利川市城區(qū)、騰龍洞、雪照河、恩施市龍鳳鎮(zhèn)、紅廟開(kāi)發(fā)區(qū)、小渡船、舞陽(yáng)壩、六角亭進(jìn)入峽口。清江流域面積 2 928.0 km2,在恩施境內(nèi)流長(zhǎng) 127.0 km,其中城區(qū)段 18.0 km,主河道坡降 7.0‰。清江干流城區(qū)河段中下游左岸有支流龍洞河、帶水河、蔡家河匯入,下游段峽口右岸有支流高橋河、巴公溪匯入。其中,清江干流進(jìn)入恩施市區(qū)時(shí),流經(jīng)了一個(gè)大型水利樞紐工程(大龍?zhí)端畮?kù))。

    龍洞河是清江左岸一級(jí)支流,發(fā)源于恩施市大樹(shù)埡鄉(xiāng)魚(yú)仙洞,流至恩施市舞陽(yáng)壩官坡與清江交匯,流域面積 38.5 km2,河長(zhǎng) 15.3.0 km,主河道坡降 1.8‰。研究區(qū)主要河流水系分布如圖 1 所示。

    圖 1 恩施市城區(qū)主要河流水系分布

    恩施市所在的清江流域?qū)儆趤啛釒Ъ撅L(fēng)性氣候,降雨充沛。多年平均降水量 1 470.1 mm,最大年降水量 1 918.4 mm。降雨年內(nèi)分配主要集中在 4—10 月,占全年水量的 70% 以上,暴雨多集中在 7—9 月,最大 24 h 降水 326.9 mm。

    在恩施市防汛抗旱指揮部辦公室的支持下,系統(tǒng)性搜集了城區(qū)地形、高分辨率遙感影像、水文氣象、社會(huì)經(jīng)濟(jì)等數(shù)據(jù),并進(jìn)行了河道斷面測(cè)量工作。其中,城區(qū)地形和水文數(shù)據(jù)來(lái)源于恩施州水文局,高分辨率遙感影像數(shù)據(jù)來(lái)源于 Google Earth 下載的 Worldview 影像,氣象數(shù)據(jù)來(lái)源于恩施州氣象局,社會(huì)經(jīng)濟(jì)數(shù)據(jù)來(lái)源于恩施市統(tǒng)計(jì)年鑒(2014 版)。為了獲得更為精細(xì)的河道斷面觀測(cè)資料,在湖北省連續(xù)運(yùn)行衛(wèi)星定位服務(wù)綜合系統(tǒng)(CORS)的支持下,使用 GPS RTK 設(shè)備,在研究區(qū)范圍內(nèi)清江干流和龍洞河河段進(jìn)行了控制斷面測(cè)量工作。

    2 洪水分析模型

    根據(jù)恩施市河流水系分析,城區(qū)主要洪水來(lái)源為清江干流。在城區(qū)中東部(舞陽(yáng)壩,恩施州人民政府所在地)附近,有支流龍洞河匯入清江。因此,恩施市主要洪水來(lái)源為:1)清江干流洪水;2)龍洞河洪水。

    2.1 洪水分析計(jì)算方案

    由于恩施市城市防洪標(biāo)準(zhǔn)定為 20 年一遇,因此主要進(jìn)行 10,20,50 ,100 年一遇 4 種洪水量級(jí)的洪水分析計(jì)算,主要考慮清江干流、龍洞河及相應(yīng)洪水組合,共設(shè)計(jì)了 8 種計(jì)算方案,如表 1 所示。

    表 1 恩施市洪水分析計(jì)算方案

    2.2 洪水分析方法

    采用水動(dòng)力學(xué)方法進(jìn)行暴雨洪水分析計(jì)算。根據(jù)恩施城區(qū)地形、水文氣象及河道斷面觀測(cè)資料,采用一維水動(dòng)力學(xué)模型進(jìn)行河道洪水演進(jìn)分析,采用二維水動(dòng)力學(xué)模型進(jìn)行河道外城區(qū)洪水淹沒(méi)分析。其中,一維水動(dòng)力學(xué)方法選用圣維南方程,二維水動(dòng)力學(xué)模型選用深度平均二維淺水波方程。構(gòu)建一、二維耦合模型開(kāi)展暴雨洪水分析能夠計(jì)算出暴雨洪水淹沒(méi)范圍和水深、淹沒(méi)歷時(shí)及洪鋒到達(dá)時(shí)間等信息,滿足洪水影響評(píng)價(jià)的需求。

    1)模型范圍。一、二維耦合模型的范圍如圖 1所示,其中藍(lán)色線條代表河道一維水動(dòng)力學(xué)模型范圍,紅色線條范圍內(nèi)的面狀區(qū)域是二維模型區(qū)域。一維河道在橫向的寬度由沿著左右河岸劃定的線條確定,河岸線內(nèi)是一維河道,河岸線外是二維網(wǎng)格區(qū)域。

    模型邊界:清江干流的上邊界范圍為大龍?zhí)端畮?kù)溢洪道出口,下邊界范圍為恩施連珠大橋下游約1 000 m 處;龍洞河的上邊界范圍為恩施東陽(yáng)酒店上游約 700 m 處,下邊界為龍洞河匯入清江的入口處。

    2)一、二維耦合模型。按照 SWMM 模型的設(shè)置要求,進(jìn)行下墊面的概化。將研究區(qū)域下墊面進(jìn)行二維網(wǎng)格劃分,構(gòu)建一、二維耦合 SWMM 模型。其中,一維河道屬于二維網(wǎng)格的非活動(dòng)區(qū)域,即位于一維河道的網(wǎng)格不參與二維運(yùn)算,河道結(jié)點(diǎn)處水位由一維河道斷面控制;二維區(qū)域通過(guò)網(wǎng)格進(jìn)行運(yùn)算,模型可將地形(.txt 格式)、初始水位、土地利用等數(shù)據(jù)賦在網(wǎng)格上。

    一、二維連接線連接河道結(jié)點(diǎn)和一、二維交界面上,兩控制斷面間河道水位由斷面處水位線性內(nèi)插獲得。當(dāng)河道內(nèi)水位超過(guò)一、二維交界面處所在網(wǎng)格的地形高程,則該處網(wǎng)格水位等于河道水位,并以該水位作為二維運(yùn)算的邊界條件;同時(shí),將二維運(yùn)算得到的流量作為一維河道的邊界條件。

    3)模型邊界條件。對(duì)于清江干流洪水分析,上游邊界條件為大龍?zhí)端畮?kù)各頻率設(shè)計(jì)下泄洪水過(guò)程線,中間加入龍洞河相應(yīng)洪水匯入。在清江最下游斷面給定實(shí)測(cè)的河道坡降,以該斷面處的水位流量關(guān)系曲線作為下邊界條件。

    對(duì)于龍洞河洪水分析,龍洞河上邊界條件為根據(jù)《湖北省暴雨徑流查算圖表》(2012 年)計(jì)算的設(shè)計(jì)洪水過(guò)程線,下邊界條件為清江 20 年一遇洪水水位頂托。

    2.3 模型參數(shù)設(shè)置

    1)模型糙率。通過(guò)實(shí)地調(diào)研過(guò)程中拍攝的照片獲取地面植被情況,參考水力學(xué)計(jì)算手冊(cè)[14]選取洪水分析模型范圍的糙率系數(shù),其參考糙率值如表 2所示。

    2)模型計(jì)算參數(shù)。水力學(xué)模型計(jì)算時(shí),清江干流各方案計(jì)算時(shí)間 72 h,二維模型的計(jì)算時(shí)間步長(zhǎng)為 3.0 s,一維模型的計(jì)算時(shí)間步長(zhǎng)為 1.5 s。龍洞河各方案計(jì)算時(shí)間為 40 h,洪峰值出現(xiàn)在第 10~12 h,二維模型的計(jì)算時(shí)間步長(zhǎng)為 3.0 s,一維模型的計(jì)算時(shí)間步長(zhǎng)為 1.5 s。

    表 2 恩施市水力學(xué)計(jì)算參考糙率

    3 洪水計(jì)算與分析

    3.1 設(shè)計(jì)洪水過(guò)程推求

    1)清江干流設(shè)計(jì)洪水過(guò)程。在一維水力學(xué)模型起點(diǎn)斷面以上的清江干流上,有大龍?zhí)端畮?kù)。因此,進(jìn)入研究區(qū)的清江干流洪水是經(jīng)過(guò)大龍?zhí)端畮?kù)調(diào)蓄后的下泄洪水。根據(jù)《清江大龍?zhí)端麡屑~初步設(shè)計(jì)報(bào)告》,不同設(shè)計(jì)洪水削峰流量如表 3 所示。

    表 3 大龍?zhí)端畮?kù)削峰流量表

    因此,經(jīng)過(guò)調(diào)洪演算之后,清江干流設(shè)計(jì)洪水過(guò)程如圖 2 所示。

    圖 2 經(jīng)過(guò)大龍?zhí)端畮?kù)調(diào)蓄后清江干流設(shè)計(jì)洪水過(guò)程線

    2)龍洞河設(shè)計(jì)暴雨洪水過(guò)程。龍洞河流域內(nèi)無(wú)實(shí)測(cè)洪水資料,且龍洞河流域面積小于 1 000 km2,適合采用《湖北省暴雨徑流查算圖表》(2002 年)推求龍洞河城區(qū)河段控制斷面10,20,50,100 年一遇的設(shè)計(jì)洪水過(guò)程線。設(shè)計(jì)洪水過(guò)程結(jié)果如圖 3所示。

    圖 3 龍洞河設(shè)計(jì)暴雨洪水過(guò)程線

    3.2 洪水計(jì)算結(jié)果

    以 100 年一遇為例。在清江入口處給定 100 年一遇洪水經(jīng)過(guò)大龍?zhí)端畮?kù)調(diào)蓄后的下泄洪水,龍洞河入口處給定 20 年一遇設(shè)計(jì)洪水過(guò)程。下游邊界條件給定實(shí)測(cè)的河道坡降,換算成下游出口的水位流量關(guān)系曲線。

    1)洪水淹沒(méi)水深、歷時(shí)。如圖 4 所示,在 100 年一遇洪水條件下,清江干流洪水淹沒(méi)范圍主要集中在恩施市小渡船、黃泥壩、施州大道、東門(mén)及南門(mén),平均淹沒(méi)水深約為 0.8~1.5 m,平均淹沒(méi)歷時(shí)約為 3~6 h;龍洞河洪水淹沒(méi)范圍主要集中在官坡、棲鳳橋及土橋壩,平均淹沒(méi)水深約為 0.5~0.8 m,平均淹沒(méi)歷時(shí)約為 2~4 h。

    圖 4 100 年一遇洪水淹沒(méi)計(jì)算結(jié)果

    2)洪峰到達(dá)時(shí)間。以洪水通過(guò)恩施城區(qū)的典型控制斷面(恩施市水文站,位于施州大橋上游 100 m處)為基礎(chǔ),計(jì)算洪峰到達(dá)時(shí)間。

    如圖 5 所示,在 100 年一遇洪水條件下,清江干流洪峰到達(dá)時(shí)間約為 6~12 h,龍洞河洪峰到達(dá)時(shí)間約為 3~6 h。

    圖 5 100 年一遇洪峰到達(dá)時(shí)間計(jì)算結(jié)果

    3.3 洪水影響評(píng)價(jià)

    1)洪水淹沒(méi)分析。如圖 6 所示,利用 ArcGIS軟件,結(jié)合高分辨率遙感影像數(shù)據(jù)和社會(huì)經(jīng)濟(jì)數(shù)據(jù),進(jìn)行恩施市洪水淹沒(méi)分析。在 10~100 年一遇洪水條件下,恩施市主要受淹區(qū)域集中在小渡船、黃泥壩、施州大道、東門(mén)、南門(mén)、官坡、棲鳳橋及土橋壩一帶,受淹區(qū)域總面積約為 2.69 km2,詳細(xì)計(jì)算結(jié)果如表 4 所示。

    2)現(xiàn)狀防洪能力評(píng)價(jià)。如圖 7 所示,利用ArcGIS 軟件結(jié)合高分辨率遙感影像和社會(huì)經(jīng)濟(jì)數(shù)據(jù),進(jìn)行恩施市現(xiàn)狀防洪能力評(píng)價(jià)。根據(jù)《山洪災(zāi)害分析評(píng)價(jià)技術(shù)要求(試行)》文件[15],將現(xiàn)狀防洪能力低于 10 a 的河段定義為極高危險(xiǎn)區(qū),現(xiàn)狀防洪能力在 10~20 a 之間的河段為高危險(xiǎn)區(qū),現(xiàn)狀防洪能力大于 20 a 的河段為危險(xiǎn)區(qū)。清江干流紅旗大橋以上基本上是裸露河床,無(wú)河道堤防工程,實(shí)際現(xiàn)狀防洪能力低于 20 a(高危險(xiǎn)區(qū));紅旗大橋至黃石大橋全部河段均修筑漿砌石堤防,堤防高度在 4~6 m 左右,受到龍洞河洪水頂托作用,實(shí)際現(xiàn)狀防洪能力在 20~50 a 之間(危險(xiǎn)區(qū));黃石大橋至連珠大橋部分河段修筑漿砌石堤防,堤防高度在 4~6 m 左右,但受到巴公溪、高橋河等洪水頂托作用,實(shí)際現(xiàn)狀防洪能力低于 10 a(極高危險(xiǎn)區(qū))。龍洞河只有部分下游河段修筑漿砌石堤防,實(shí)際現(xiàn)狀防洪能力在10~20 a 之間(高危險(xiǎn)區(qū))。

    圖 6 恩施市城區(qū)主要洪水淹沒(méi)區(qū)

    圖 7 恩施市城區(qū)現(xiàn)狀防洪能力評(píng)價(jià)

    4 結(jié)語(yǔ)

    1)SWMM 模型是一種成熟的城市雨洪徑流分析模型,能夠較為方便地構(gòu)建一、二維耦合水力學(xué)模型,快速計(jì)算河道洪水過(guò)程和洪水淹沒(méi)分布,適宜于城市暴雨洪水模擬與分析。

    2)根據(jù)恩施市城區(qū)河流水系、洪水來(lái)源、地形地貌等情況,設(shè)計(jì)了一套洪水分析計(jì)算方案,并基于 SWMM 模型進(jìn)行雨洪過(guò)程模擬分析,取得了較為準(zhǔn)確的洪水淹沒(méi)水深、淹沒(méi)歷時(shí)、洪峰到達(dá)時(shí)間等結(jié)果,為城市雨洪淹沒(méi)分析和現(xiàn)狀防洪能力評(píng)價(jià)提供依據(jù)。結(jié)果表明,清江干流主要河段基本滿足了 20 年一遇的防洪標(biāo)準(zhǔn),清江干流和龍洞河部分河段現(xiàn)狀防洪能力較低。這一結(jié)果與恩施市防汛抗旱指揮部辦公室發(fā)布的《恩施市城區(qū)防洪應(yīng)急預(yù)案》[16]基本吻合。

    表 4 恩施市雨洪淹沒(méi)面積統(tǒng)計(jì)表

    3)下一步準(zhǔn)備通過(guò)重點(diǎn)單位走訪和實(shí)地調(diào)研,盡可能地搜集完善恩施市排水管網(wǎng)空間分布、管徑大小等數(shù)據(jù)資料,利用 SWMM 排水管網(wǎng)模塊進(jìn)一步分析城市洪水內(nèi)澇消退過(guò)程。

    [1] 張冬冬,嚴(yán)登華,王義成,等. 城市內(nèi)澇災(zāi)害風(fēng)險(xiǎn)評(píng)估及綜合應(yīng)對(duì)研究進(jìn)展[J]. 災(zāi)害學(xué),2014 (1): 144-149.

    [2] 戴晶晶,劉增賢,陸沈釣. 基于數(shù)值模擬的城市內(nèi)澇風(fēng)險(xiǎn)評(píng)估研究—以蘇州市城市中心區(qū)為例[J]. 中國(guó)水利,2015 (13): 20-23.

    [3] 胡偉賢,何文華,黃國(guó)如,等. 城市雨洪模擬技術(shù)研究進(jìn)展[J]. 水科學(xué)進(jìn)展,2010 (1): 137-144.

    [4] 秦語(yǔ)涵,王紅武,張一龍. 城市雨洪徑流模型研究進(jìn)展[J]. 環(huán)境科學(xué)與技術(shù),2016 (1): 13-19.

    [5] CHIRISTOPHER Z. Review of urban storm watermodels[J].Enviornmental Modeling & Software, 2001 (16):195-231.

    [6] GORDON S I. Stormwater Runoff Models [M]. U. S. :Springer, 1985: 145-162.

    [7] HUANG G, WU S Y. Simulation study on effect of rainwater utilization measures on urban stormwater based on infoworks CS[J]. Water Resources and Power, 2013 (5): 65-71.

    [8] ZAGHLOUL N A, M A Abu Kiefa. Neural network solution of inverse parameters used in the sensitivity-calibration analysis of the SWMM model simulations[J]. Advances in Engineering Software, 2001, 32 (7): 87-595.

    [9] SMITH D, LI J,BANTING D. A PCSWMM/GIS-based water balance model for the reesor Creek watershed[J].Atmospheric Research, 2005, 77 (1): 338-406.

    [10] 董欣,杜鵬飛,李志一,等. SWMM 模型在城市不透水區(qū)地表徑流模擬中的參數(shù)識(shí)別與驗(yàn)證[J]. 環(huán)境科學(xué),2008 (6): 1495-1501.

    [11] KREBS G, KOKKONEN T, VALTANEN M, et al. A high resolution application of a storm water management model(SWMM) using genetic parameter optimization[J]. Urban Water Journal, 2013,10 (6): 394-410.

    [12] 周力寧. 基于 SWMM 的城市內(nèi)澇風(fēng)險(xiǎn)識(shí)別研究[M]. 成都:西南交通大學(xué),2016: 35-44.

    [13] 申邵洪,江炎生,蔡斌,等. 城鎮(zhèn)區(qū)域洪水風(fēng)險(xiǎn)分析研究-以湖北省保康縣為例[J]. 長(zhǎng)江科學(xué)院院報(bào),2016,33 (11): 109-115.

    [14] 李煒. 水力計(jì)算手冊(cè) [M]. 2 版. 北京:中國(guó)水利水電出版社,2006: 53-65.

    [15] 全國(guó)山洪災(zāi)害防治組. 山洪災(zāi)害分析評(píng)價(jià)技術(shù)要求(試行)[M]. 北京:國(guó)家防汛抗旱總指揮部辦公室,2014:79-112.

    [16] 恩施州水文水資源勘測(cè)局. 恩施市城區(qū)防洪應(yīng)急預(yù)案[M].恩施:恩施市防汛抗旱指揮部辦公室,2009: 33-54.

    Urban storm flood simulation and analysis based on storm water management model

    YANG Bin1, ZHAO Rui1,2

    (1. Hubei University of Technology, Wuhan 430068, China;2. Middler Changjiang River Bureau of Hydrological and Water Resources Survey, Wuhan 430010, China)

    In order to conduct urban storm flood simulation and analysis, the paper starts with river system distribution and flood source analysis in case study of Enshi City, and a flood calculation and analysis program are proposed. On the basis of storm water management model (SWMM), an one dimension and two dimension coupled hydraulics model is constructed, which is applied to calculate river flood evolution and diffuse embankment process under the designed flood conditions of 0.01, 0.02, 0.05 and 0.10 (probability). The comprehensive information of flood area, water depth, duration and peak arrival time is also retrieved, which is used to calculate flood effect analysis.Results reveal that the main parts of Qingjiang River basically meet the flood prevention standard of 0.05 (probability)in Enshi City, while a small section of Qingjiang River and Longdong River has lower flood prevention ability than its flood prevention standard. And results also suggest that SWMM is a reliable urban storm water model, which is suitable for urban storm flood simulation and analysis applications.

    storm water management model (SWMM); one dimension and two dimension couple; storm flood monitoring; flood effect analysis

    TV223

    A

    1674-9405(2017)05-0056-07

    10.19364/j.1674-9405.2017.05.011

    2017-06-14

    楊 斌(1969-),男,湖南常德人,教授,主要從事管理科學(xué)與工程、智慧城市等方面研究。

    猜你喜歡
    龍洞恩施市清江
    蘇元老撰并楷書(shū)《龍洞記碑》
    A Review of Studies since the 1980’s on the Tieto-urman Song of the White Wolf
    清江引
    影劇新作(2018年1期)2018-05-26 09:00:52
    魚(yú)躍清江 廣場(chǎng)舞
    文化交流(2017年8期)2017-09-14 22:03:27
    孽龍洞玄想
    星火(2016年5期)2016-02-11 03:31:47
    淺談古村落的保護(hù)與開(kāi)發(fā)——以湖北恩施市盛家壩鄉(xiāng)小溪村為例
    民族大家庭(2015年2期)2015-08-21 07:59:54
    同飲清江水 共護(hù)母親河——首個(gè)“清江保護(hù)日”在長(zhǎng)陽(yáng)舉行
    恩施市基層就業(yè)創(chuàng)業(yè)工作的對(duì)策
    生態(tài)環(huán)境保護(hù)的實(shí)踐與探索——以恩施市屯堡鄉(xiāng)為例
    恩施市中心醫(yī)院:真心服務(wù)體現(xiàn)品牌價(jià)值
    18禁观看日本| 欧美日韩精品网址| 无遮挡黄片免费观看| 国产一区二区激情短视频| 一边摸一边抽搐一进一小说| 国产激情久久老熟女| 国产亚洲精品一区二区www| 妹子高潮喷水视频| 亚洲18禁久久av| 久久亚洲真实| 国产人伦9x9x在线观看| 久久国产精品影院| 国产一级毛片七仙女欲春2| 无限看片的www在线观看| 国产精品 国内视频| 欧美性猛交╳xxx乱大交人| av福利片在线| 国产99久久九九免费精品| 动漫黄色视频在线观看| 国产欧美日韩精品亚洲av| 在线观看免费日韩欧美大片| 亚洲av成人精品一区久久| 久热爱精品视频在线9| 一a级毛片在线观看| 搡老妇女老女人老熟妇| 国产精品久久电影中文字幕| 欧美日韩精品网址| 日韩欧美一区二区三区在线观看| 母亲3免费完整高清在线观看| 久久草成人影院| 亚洲成av人片在线播放无| 亚洲专区中文字幕在线| 舔av片在线| 三级毛片av免费| 一级毛片女人18水好多| 在线观看www视频免费| 亚洲激情在线av| 欧美黑人巨大hd| 婷婷丁香在线五月| 亚洲欧美精品综合一区二区三区| 美女黄网站色视频| 成人av在线播放网站| 中文字幕久久专区| 午夜福利在线在线| 淫秽高清视频在线观看| 首页视频小说图片口味搜索| 欧美激情久久久久久爽电影| 亚洲国产精品sss在线观看| 免费观看人在逋| 成年免费大片在线观看| 欧美精品啪啪一区二区三区| 婷婷精品国产亚洲av在线| 草草在线视频免费看| 婷婷亚洲欧美| 女警被强在线播放| 国产精品久久电影中文字幕| 白带黄色成豆腐渣| 88av欧美| 免费在线观看亚洲国产| 日本熟妇午夜| 免费av毛片视频| 国产精品久久视频播放| 免费在线观看日本一区| 日韩精品青青久久久久久| 久久精品综合一区二区三区| 国产成人欧美在线观看| 精品不卡国产一区二区三区| 亚洲美女视频黄频| 亚洲国产欧美一区二区综合| 最新美女视频免费是黄的| 午夜两性在线视频| 亚洲18禁久久av| 欧美人与性动交α欧美精品济南到| 精品国内亚洲2022精品成人| 亚洲乱码一区二区免费版| 精品国产乱码久久久久久男人| 国产视频内射| 精品福利观看| 亚洲av成人av| e午夜精品久久久久久久| 国产成+人综合+亚洲专区| 亚洲成人精品中文字幕电影| 成人午夜高清在线视频| 最新美女视频免费是黄的| 亚洲精品在线美女| 国产精品自产拍在线观看55亚洲| 一卡2卡三卡四卡精品乱码亚洲| 老司机在亚洲福利影院| 国产欧美日韩精品亚洲av| 久久精品成人免费网站| 免费看日本二区| 九色成人免费人妻av| 久久精品人妻少妇| 一进一出好大好爽视频| 99精品欧美一区二区三区四区| 女同久久另类99精品国产91| 久久久久久人人人人人| 午夜免费成人在线视频| 国内揄拍国产精品人妻在线| 国产精品免费一区二区三区在线| 欧美日韩亚洲综合一区二区三区_| 在线观看午夜福利视频| 最新在线观看一区二区三区| 国产高清videossex| 真人做人爱边吃奶动态| 亚洲自偷自拍图片 自拍| 国语自产精品视频在线第100页| 一区二区三区高清视频在线| 在线观看一区二区三区| 久久久久久国产a免费观看| 一个人观看的视频www高清免费观看 | x7x7x7水蜜桃| 黑人欧美特级aaaaaa片| 特级一级黄色大片| 亚洲熟女毛片儿| av天堂在线播放| 丰满的人妻完整版| a级毛片a级免费在线| 狠狠狠狠99中文字幕| 国产欧美日韩一区二区精品| 欧美乱码精品一区二区三区| 丰满的人妻完整版| 亚洲人成伊人成综合网2020| 亚洲专区国产一区二区| 亚洲中文av在线| 久久久国产成人精品二区| 激情在线观看视频在线高清| 亚洲熟妇中文字幕五十中出| 免费看日本二区| 国产精品免费视频内射| 黄色视频不卡| 特大巨黑吊av在线直播| 日韩欧美 国产精品| 国产精品精品国产色婷婷| 两人在一起打扑克的视频| 高潮久久久久久久久久久不卡| 黄片小视频在线播放| 视频区欧美日本亚洲| 亚洲色图 男人天堂 中文字幕| 看片在线看免费视频| 男女午夜视频在线观看| 99riav亚洲国产免费| 国产激情欧美一区二区| 国产一区二区在线观看日韩 | 精品人妻1区二区| 亚洲精品中文字幕一二三四区| 熟女少妇亚洲综合色aaa.| 国产精品久久久久久久电影 | 亚洲片人在线观看| 精品熟女少妇八av免费久了| 激情在线观看视频在线高清| 国产三级黄色录像| 国产久久久一区二区三区| 999精品在线视频| 舔av片在线| 欧美中文综合在线视频| 成人国产一区最新在线观看| 日本 欧美在线| 中文字幕久久专区| 久久精品夜夜夜夜夜久久蜜豆 | 婷婷精品国产亚洲av| 男女床上黄色一级片免费看| 成年女人毛片免费观看观看9| 一二三四在线观看免费中文在| 亚洲18禁久久av| 夜夜夜夜夜久久久久| 日韩欧美国产一区二区入口| 日本 欧美在线| 日日夜夜操网爽| 女人高潮潮喷娇喘18禁视频| 亚洲七黄色美女视频| 最近最新中文字幕大全电影3| 一区二区三区高清视频在线| 亚洲精华国产精华精| 一级毛片女人18水好多| 成人三级黄色视频| 久久香蕉激情| 国产熟女xx| 又粗又爽又猛毛片免费看| 国产精品一区二区三区四区免费观看 | 精品一区二区三区四区五区乱码| 日本一本二区三区精品| 小说图片视频综合网站| 国产成人系列免费观看| 久久中文字幕人妻熟女| 一进一出抽搐动态| 国产探花在线观看一区二区| 老汉色∧v一级毛片| 中文字幕人妻丝袜一区二区| 18美女黄网站色大片免费观看| 校园春色视频在线观看| 伊人久久大香线蕉亚洲五| 狠狠狠狠99中文字幕| 午夜精品一区二区三区免费看| 亚洲真实伦在线观看| 母亲3免费完整高清在线观看| 窝窝影院91人妻| 国产一区在线观看成人免费| 亚洲成av人片免费观看| 搞女人的毛片| 国产一区二区激情短视频| 校园春色视频在线观看| 国产亚洲精品第一综合不卡| 真人一进一出gif抽搐免费| 国产伦在线观看视频一区| 宅男免费午夜| 91麻豆av在线| 午夜视频精品福利| 这个男人来自地球电影免费观看| 亚洲va日本ⅴa欧美va伊人久久| 波多野结衣高清作品| 中出人妻视频一区二区| 91字幕亚洲| 亚洲国产精品成人综合色| 啦啦啦观看免费观看视频高清| 成人三级黄色视频| 日本一区二区免费在线视频| 国产精品99久久99久久久不卡| 久久久久久久午夜电影| 国产久久久一区二区三区| 日韩欧美国产一区二区入口| 最近在线观看免费完整版| 无限看片的www在线观看| 亚洲18禁久久av| 国产v大片淫在线免费观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av熟女| 人妻丰满熟妇av一区二区三区| www日本黄色视频网| 精品熟女少妇八av免费久了| 久久精品91蜜桃| 91av网站免费观看| 国产1区2区3区精品| 草草在线视频免费看| 免费在线观看成人毛片| 黄色 视频免费看| 最近最新中文字幕大全免费视频| 国产熟女xx| 国产精品一区二区三区四区久久| 亚洲国产精品久久男人天堂| 精品人妻1区二区| 国产亚洲av嫩草精品影院| 亚洲欧美日韩东京热| 久久精品影院6| √禁漫天堂资源中文www| 天天躁夜夜躁狠狠躁躁| 十八禁人妻一区二区| 欧美绝顶高潮抽搐喷水| 国产av不卡久久| 蜜桃久久精品国产亚洲av| 久久久久亚洲av毛片大全| 嫁个100分男人电影在线观看| 欧美 亚洲 国产 日韩一| 十八禁网站免费在线| 欧美中文综合在线视频| 啦啦啦免费观看视频1| 在线观看免费午夜福利视频| 欧美日韩黄片免| 男女视频在线观看网站免费 | 男女视频在线观看网站免费 | 亚洲精品久久国产高清桃花| 无限看片的www在线观看| 亚洲熟妇熟女久久| 国产高清videossex| 亚洲国产精品sss在线观看| 久99久视频精品免费| 小说图片视频综合网站| 一级片免费观看大全| 亚洲av片天天在线观看| 男女床上黄色一级片免费看| 最新美女视频免费是黄的| 五月伊人婷婷丁香| 亚洲成a人片在线一区二区| 老司机靠b影院| 18美女黄网站色大片免费观看| 久久精品亚洲精品国产色婷小说| 最近视频中文字幕2019在线8| 欧美一级a爱片免费观看看 | 欧美中文日本在线观看视频| 首页视频小说图片口味搜索| 丝袜人妻中文字幕| 两人在一起打扑克的视频| 欧美精品啪啪一区二区三区| 亚洲精品久久成人aⅴ小说| 国产精华一区二区三区| 日日夜夜操网爽| 国产三级在线视频| 69av精品久久久久久| 最新美女视频免费是黄的| 日本免费一区二区三区高清不卡| 国产精品av久久久久免费| 亚洲avbb在线观看| 精品一区二区三区视频在线观看免费| 琪琪午夜伦伦电影理论片6080| 国产精品亚洲av一区麻豆| 三级男女做爰猛烈吃奶摸视频| 亚洲精品美女久久av网站| x7x7x7水蜜桃| 亚洲精品在线美女| 国产亚洲av嫩草精品影院| 精品福利观看| 男女视频在线观看网站免费 | 国产精品av久久久久免费| 一级毛片精品| 欧美黑人欧美精品刺激| 久久人妻av系列| 国产精品av视频在线免费观看| 可以免费在线观看a视频的电影网站| 精品电影一区二区在线| 操出白浆在线播放| 欧美三级亚洲精品| 视频区欧美日本亚洲| 2021天堂中文幕一二区在线观| 丰满人妻一区二区三区视频av | 精品久久久久久久久久久久久| 悠悠久久av| 久久中文字幕人妻熟女| 在线播放国产精品三级| 大型黄色视频在线免费观看| 国产伦人伦偷精品视频| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 俄罗斯特黄特色一大片| 夜夜夜夜夜久久久久| 亚洲人成网站在线播放欧美日韩| 91麻豆精品激情在线观看国产| 麻豆成人午夜福利视频| 日韩大尺度精品在线看网址| 午夜福利高清视频| 久久中文字幕人妻熟女| 亚洲第一欧美日韩一区二区三区| 熟女电影av网| 久久久精品欧美日韩精品| 蜜桃久久精品国产亚洲av| 亚洲一区二区三区色噜噜| 欧美性猛交╳xxx乱大交人| 色尼玛亚洲综合影院| 国产三级中文精品| 免费搜索国产男女视频| 国产av麻豆久久久久久久| 国产在线观看jvid| 禁无遮挡网站| 一本大道久久a久久精品| 天堂动漫精品| 中出人妻视频一区二区| 亚洲成人久久性| 亚洲一码二码三码区别大吗| 欧美乱色亚洲激情| 制服丝袜大香蕉在线| 在线国产一区二区在线| 嫩草影院精品99| 成人精品一区二区免费| av有码第一页| 国产亚洲精品av在线| 久久久久久久午夜电影| 91在线观看av| 在线观看美女被高潮喷水网站 | 黄色视频不卡| 国产精品永久免费网站| 欧美日本亚洲视频在线播放| 欧美黑人巨大hd| 日日摸夜夜添夜夜添小说| 午夜福利18| 亚洲黑人精品在线| 国产精品 欧美亚洲| 日本黄大片高清| 国产高清有码在线观看视频 | 男人舔女人的私密视频| 欧美在线一区亚洲| 亚洲精品av麻豆狂野| 国产高清激情床上av| bbb黄色大片| 国产精品99久久99久久久不卡| 一个人免费在线观看的高清视频| 蜜桃久久精品国产亚洲av| 亚洲av成人不卡在线观看播放网| 在线观看午夜福利视频| 欧美乱色亚洲激情| 国产欧美日韩精品亚洲av| 亚洲男人的天堂狠狠| 久久欧美精品欧美久久欧美| a在线观看视频网站| 欧美性长视频在线观看| 天堂√8在线中文| 亚洲国产中文字幕在线视频| 神马国产精品三级电影在线观看 | 午夜两性在线视频| cao死你这个sao货| 国产亚洲精品一区二区www| 久久精品综合一区二区三区| 制服人妻中文乱码| 男女下面进入的视频免费午夜| 99精品久久久久人妻精品| 亚洲中文字幕日韩| av视频在线观看入口| 精品国内亚洲2022精品成人| 日本 av在线| 99久久精品国产亚洲精品| 夜夜看夜夜爽夜夜摸| 欧美成人性av电影在线观看| 日日夜夜操网爽| 亚洲男人天堂网一区| 成人三级做爰电影| 午夜福利在线在线| 国产真实乱freesex| 亚洲中文字幕日韩| 精品免费久久久久久久清纯| 中文字幕高清在线视频| 亚洲,欧美精品.| 美女扒开内裤让男人捅视频| 在线观看免费视频日本深夜| 日韩欧美免费精品| 高潮久久久久久久久久久不卡| 国产蜜桃级精品一区二区三区| 国产日本99.免费观看| 黄色 视频免费看| 麻豆av在线久日| 欧美性长视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲第一电影网av| 亚洲欧美日韩无卡精品| 热99re8久久精品国产| 国产亚洲精品久久久久久毛片| 亚洲七黄色美女视频| 波多野结衣巨乳人妻| 免费在线观看完整版高清| 亚洲av电影不卡..在线观看| 欧美日本亚洲视频在线播放| 美女扒开内裤让男人捅视频| 午夜免费观看网址| 一级毛片精品| 一个人免费在线观看电影 | 亚洲人成电影免费在线| 狂野欧美白嫩少妇大欣赏| 777久久人妻少妇嫩草av网站| 国产亚洲欧美在线一区二区| 亚洲欧洲精品一区二区精品久久久| 国产野战对白在线观看| 国产真实乱freesex| 亚洲中文字幕日韩| 一级片免费观看大全| 特级一级黄色大片| 天天添夜夜摸| 亚洲av中文字字幕乱码综合| 黄片小视频在线播放| 国产久久久一区二区三区| 国产在线观看jvid| 国产亚洲欧美98| 日韩欧美国产一区二区入口| 国产亚洲av嫩草精品影院| 亚洲成av人片在线播放无| 亚洲欧美一区二区三区黑人| 日本成人三级电影网站| 最近在线观看免费完整版| 欧美又色又爽又黄视频| 熟女少妇亚洲综合色aaa.| 国产精品,欧美在线| 怎么达到女性高潮| 久久中文字幕人妻熟女| 亚洲专区国产一区二区| 欧美日韩福利视频一区二区| 亚洲欧美日韩无卡精品| 亚洲天堂国产精品一区在线| 91麻豆精品激情在线观看国产| 日本在线视频免费播放| 国产成人系列免费观看| 五月伊人婷婷丁香| 非洲黑人性xxxx精品又粗又长| 久久精品91蜜桃| www.www免费av| 亚洲av五月六月丁香网| 桃红色精品国产亚洲av| 欧美又色又爽又黄视频| 国产久久久一区二区三区| 十八禁网站免费在线| 韩国av一区二区三区四区| 男女视频在线观看网站免费 | 日韩大码丰满熟妇| 天天躁夜夜躁狠狠躁躁| 色av中文字幕| 国产精品一区二区三区四区久久| 欧美乱码精品一区二区三区| 久久精品91无色码中文字幕| 91字幕亚洲| 欧美极品一区二区三区四区| 十八禁网站免费在线| 国产伦在线观看视频一区| 午夜福利在线观看吧| 国内少妇人妻偷人精品xxx网站 | 欧美日韩精品网址| 别揉我奶头~嗯~啊~动态视频| 人人妻,人人澡人人爽秒播| 美女扒开内裤让男人捅视频| 亚洲国产欧洲综合997久久,| 国产熟女xx| 中国美女看黄片| 亚洲 国产 在线| 亚洲精品国产一区二区精华液| 国产精品一区二区三区四区久久| 国产免费男女视频| 毛片女人毛片| 天天躁狠狠躁夜夜躁狠狠躁| 国产在线精品亚洲第一网站| 久久久精品欧美日韩精品| 国产精品98久久久久久宅男小说| 日本一区二区免费在线视频| 在线观看一区二区三区| 国产激情欧美一区二区| 99国产精品一区二区蜜桃av| 悠悠久久av| 精品午夜福利视频在线观看一区| 免费高清视频大片| 757午夜福利合集在线观看| tocl精华| 五月伊人婷婷丁香| 精品国内亚洲2022精品成人| av视频在线观看入口| 老汉色∧v一级毛片| 精品久久蜜臀av无| 亚洲精品色激情综合| 亚洲专区字幕在线| 高清毛片免费观看视频网站| 午夜视频精品福利| 亚洲男人的天堂狠狠| 日韩三级视频一区二区三区| 国产精品国产高清国产av| 好男人在线观看高清免费视频| 久久久国产欧美日韩av| 淫妇啪啪啪对白视频| 国产精品精品国产色婷婷| 国产视频一区二区在线看| xxxwww97欧美| 午夜免费观看网址| 国产伦在线观看视频一区| 在线观看免费日韩欧美大片| 最好的美女福利视频网| 中文字幕最新亚洲高清| 成在线人永久免费视频| 一边摸一边做爽爽视频免费| 99热6这里只有精品| 丰满的人妻完整版| 中文资源天堂在线| 黄片大片在线免费观看| 久久人人精品亚洲av| 99热这里只有是精品50| 日本免费a在线| 亚洲电影在线观看av| 少妇的丰满在线观看| 一二三四社区在线视频社区8| 成人高潮视频无遮挡免费网站| 亚洲一卡2卡3卡4卡5卡精品中文| 18禁黄网站禁片免费观看直播| 国内精品久久久久精免费| 国产伦人伦偷精品视频| 99久久久亚洲精品蜜臀av| 99热这里只有精品一区 | cao死你这个sao货| 欧美+亚洲+日韩+国产| 国产高清videossex| 叶爱在线成人免费视频播放| 亚洲中文字幕一区二区三区有码在线看 | 中出人妻视频一区二区| 变态另类丝袜制服| 黄色成人免费大全| 免费在线观看成人毛片| 亚洲激情在线av| 欧美乱色亚洲激情| 午夜精品久久久久久毛片777| 日本 av在线| 亚洲av片天天在线观看| 男插女下体视频免费在线播放| 2021天堂中文幕一二区在线观| 亚洲av成人精品一区久久| 两个人视频免费观看高清| 欧美一区二区精品小视频在线| 久久伊人香网站| 俄罗斯特黄特色一大片| 麻豆久久精品国产亚洲av| 亚洲男人的天堂狠狠| 国产一区二区三区在线臀色熟女| 一边摸一边做爽爽视频免费| 久久久国产成人精品二区| 99在线人妻在线中文字幕| 777久久人妻少妇嫩草av网站| 国产日本99.免费观看| 国产精品免费一区二区三区在线| 欧美性长视频在线观看| 一进一出抽搐动态| 国产精品,欧美在线| 亚洲最大成人中文| 99精品久久久久人妻精品| 久久99热这里只有精品18| 香蕉久久夜色| 久久香蕉激情| 他把我摸到了高潮在线观看| 香蕉久久夜色| 一级作爱视频免费观看| 草草在线视频免费看| 亚洲专区字幕在线| 免费在线观看黄色视频的| 国产91精品成人一区二区三区| 亚洲最大成人中文| 非洲黑人性xxxx精品又粗又长| 欧美黑人欧美精品刺激| 国产黄片美女视频| 亚洲一区高清亚洲精品| 女生性感内裤真人,穿戴方法视频| 一二三四在线观看免费中文在| 欧美丝袜亚洲另类 | 日韩有码中文字幕| 香蕉av资源在线| 欧美日韩瑟瑟在线播放| 美女 人体艺术 gogo| 精品乱码久久久久久99久播| 婷婷亚洲欧美|