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

    基于SWMM耦合HEC-RAS 的人工建設(shè)水網(wǎng)防洪排澇模擬研究

    2024-01-01 00:00:00夏銘輝尹世洋
    人民珠江 2024年6期
    關(guān)鍵詞:水網(wǎng)閘門(mén)

    摘要:對(duì)新建城區(qū)的水文、水動(dòng)力過(guò)程進(jìn)行準(zhǔn)確模擬分析是人工水網(wǎng)規(guī)劃建設(shè)及新建城區(qū)內(nèi)澇防治的關(guān)鍵。應(yīng)用 SWMM 分布式水文模型耦合 HEC-RAS 二維水動(dòng)力 SWE-ELW模型,對(duì)安徽省阜陽(yáng)市城南水網(wǎng)進(jìn)行模型模擬及防洪排澇分析。結(jié)果表明,SWMM模型中計(jì)算參數(shù)更多,排澇模數(shù)隨匯水分區(qū)面積增大而減小,且與規(guī)劃泵站排澇模數(shù)接近,計(jì)算結(jié)果更為合理。城南水網(wǎng)排澇標(biāo)準(zhǔn)30年一遇工況下,城南水網(wǎng)水系汛期最高水位為27.89 m,滿足排澇方案要求;各閘門(mén)過(guò)流流量與其河道底寬成正比,各閘門(mén)各過(guò)程線波動(dòng)程度與其河道底寬成反比;水網(wǎng)先蓄后排削峰57.21 m3/s,錯(cuò)峰2.83 h,調(diào)蓄洪量169萬(wàn) m3,削峰率42.45%,防洪排澇效益顯著。該方法具備一定的模擬效果及精度,可為相關(guān)人工建設(shè)水網(wǎng)的防洪排澇研究提供借鑒和指導(dǎo)。

    關(guān)鍵詞:水網(wǎng);防洪排澇;閘門(mén);SWMM;HEC-RAS;二維模型

    中圖分類(lèi)號(hào):TV213.9;TU992 文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):1001-9235(2024)06-0030-09

    Research on Flood Control and Drainage Simulation of Artificial Water Networks Based on SWMM Coupled with HEC-RAS

    XIA Minghui1, YIN Shiyang2

    (1. Beijing Youxinhui Engineering Consulting Corporation, Beijing 100026, China;2. College of Water Resources and HydropowerEngineering, North China Electric Power University, Beijing 102206, China)

    Abstract: Accurate simulation and analysis of hydrology and hydraulic processes in newly constructed urban areas are crucial for the planning and construction of artificial water networks and the prevention of flooding in newly constructed urban areas. A distributed hydrological model of SWMM was coupled with a two-dimensional hydrodynamic model (SWE-ELW) of HEC-RAS to simulate and analyze flood control and drainage in the southern water network of Fuyang City in Anhui Province. The results show that the number of calculation parameters in the SWMM model is more, and the drainage modulus decreases with increasing catchment area, which is close to the drainage modulus in planned pump stations, making the calculation results more reasonable. Under the condition of a 30- year return period of the southern water network drainage standard, the highest water level in the southern water network reaches 27.89 m during flood season, meeting the requirements of the drainage scheme. The flow through each gate is proportional to the river channel′s bottom width, and the fluctuation degree of each process line of each gate is inversely proportional to the river channel′s bottom width. The network first stores water before discharging it, reducing peak flow by 57.21 m3/s, delaying peak flow by 2.83 hours, storing flood by 1.69 million m3, and reaching a peak flow reduction rate of 42.45%, resulting in significant flood control and drainage benefits. This method has certain simulation effects and accuracy, which can provide reference and guidance for relatedstudies on flood control and drainage in artificial water network construction.

    Keywords: water network; flood control and drainage; gate; SWMM; HEC-RAS; two-dimensional model

    近期,中共中央、國(guó)務(wù)院提出加快構(gòu)建國(guó)家水網(wǎng)、建設(shè)現(xiàn)代化高質(zhì)量水利基礎(chǔ)設(shè)施網(wǎng)絡(luò),印發(fā)《國(guó)家水網(wǎng)建設(shè)規(guī)劃綱要》,對(duì)水系連通、綜合治理及引調(diào)排水等工程建設(shè)提出了新的戰(zhàn)略部署和要求。防洪排澇是水網(wǎng)承擔(dān)的重要功能,如何量化評(píng)估水網(wǎng)防洪排澇效益,對(duì)水網(wǎng)設(shè)計(jì)和建設(shè)至關(guān)重要。數(shù)值模型能模擬反映水網(wǎng)水系復(fù)雜的結(jié)構(gòu)及水力特性,為研究防洪排澇問(wèn)題的主要技術(shù)方法。

    部分學(xué)者應(yīng)用數(shù)值模型于天然水網(wǎng)的防洪排澇研究,取得了較好的研究成果。徐馳等[1-2]應(yīng)用MIKE FLOOD模型分析南京市高淳區(qū)河湖水網(wǎng),利用 ArcGIS、MIKE、MATLAB 等工具建立一維二維耦合水網(wǎng)模型,實(shí)現(xiàn)防洪、供水、水環(huán)境治理的統(tǒng)一模擬,構(gòu)建生態(tài)水網(wǎng)評(píng)價(jià)指標(biāo)體系對(duì)南京市高淳區(qū)生態(tài)水網(wǎng)進(jìn)行定量與定性分析。范立柱等[3]應(yīng)用MIKE11一維模型對(duì)九江鎮(zhèn)水系優(yōu)化排澇方案進(jìn)行可行性驗(yàn)算。閆偉偉等[4]構(gòu)建河網(wǎng)規(guī)劃概化圖,建立湖泊調(diào)蓄演算模型來(lái)擬定武漢市新洲區(qū)漲渡湖流域防洪排澇規(guī)劃方案。袁玉等[5]應(yīng)用 HEC-HMS 水文模型對(duì)秦淮河流域降雨徑流模擬,定量計(jì)算圩垸式防洪削減單次洪水洪峰值及洪量。蔣建靈[6]建立杭州市運(yùn)西片河網(wǎng)一維非恒定流水動(dòng)力學(xué)模型,選用2013年10月的“菲特”臺(tái)風(fēng)數(shù)據(jù)對(duì)深隧排洪方案及排澇效果進(jìn)行模擬驗(yàn)證。

    一些城市在城市建設(shè)過(guò)程中,通過(guò)水系連通及綜合治理將原本獨(dú)立的河流連通形成水網(wǎng),同時(shí)以水定城,依水建城。新建城區(qū)的防洪排澇不僅需要考慮開(kāi)發(fā)建設(shè)導(dǎo)致的下墊面變化,而且也要兼顧人工河網(wǎng)、用地變化及市政雨水管網(wǎng)的復(fù)雜水力關(guān)系。

    對(duì)新建城區(qū)的水文、水動(dòng)力過(guò)程進(jìn)行準(zhǔn)確模擬分析是人工水網(wǎng)規(guī)劃建設(shè)及新建城區(qū)內(nèi)澇防治的關(guān)鍵。目前國(guó)內(nèi)外的相關(guān)研究主要集中于 SWMM[7-8]、MIKE 和 Info Works ICM 等模型。符銳[9]通過(guò) SWMM 模型研究石南省昭通市現(xiàn)狀水系及規(guī)劃河道排水防澇規(guī)劃方案。吳沛霖等[10]構(gòu)建張家港城市排水防澇 SWMM模型,運(yùn)用指標(biāo)體系法定義排水防澇風(fēng)險(xiǎn)級(jí)別,繪制不同重現(xiàn)期內(nèi)澇風(fēng)險(xiǎn)圖。劉瑤[11]應(yīng)用 MIKE FLOOD 模型評(píng)估深圳市福田河流域排水管網(wǎng)系排水能力及城市內(nèi)澇風(fēng)險(xiǎn)。李尤等[12]利用 Info Works ICM 構(gòu)建北京市清河內(nèi)典型流域老龍口溝精細(xì)化綜合洪澇耦合模型,分析顯示該流域內(nèi)澇積水問(wèn)題由不透水下墊面比例高、管網(wǎng)排水能力弱、局部地形低洼、局地強(qiáng)降雨易發(fā)等多因素綜合造成。許婷婷等[13]采用 Info Works ICM模型軟件構(gòu)建雄安新區(qū)起步區(qū)洪澇模型進(jìn)行分析,結(jié)果顯示對(duì)河道水位充分預(yù)降可最大化利用規(guī)劃泵站的排水能力。葉陳雷等[14]基于 Info Works ICM 構(gòu)建福州主城區(qū)晉安河片區(qū)水文水動(dòng)力耦合模型,對(duì)水工程聯(lián)合調(diào)度以及低影響開(kāi)發(fā)作用下的洪澇過(guò)程進(jìn)行了量化分析,結(jié)果顯示通過(guò)聯(lián)合調(diào)度對(duì)河道水位調(diào)節(jié)可以避免河道漫溢。

    目前人工建設(shè)河網(wǎng)及新建城區(qū)內(nèi)澇防治研究主要集中于 SWMM、MIKE 和 Info Works ICM 等模型,HEC-RAS 模型相關(guān)應(yīng)用及研究甚少。 HEC- RAS模型為美國(guó)工程兵團(tuán)開(kāi)發(fā)的水動(dòng)力模型,廣泛應(yīng)用于洪水模擬相關(guān)研究[15-18]。該模型在一維二維水動(dòng)力計(jì)算及水工建筑物調(diào)度模擬上性能較為優(yōu)越。

    本文應(yīng)用 SWMM 分布式水文模型耦合 HEC- RAS 二維水動(dòng)力 SWE-ELW 模型,對(duì)安徽省阜陽(yáng)市城南水網(wǎng)進(jìn)行模型模擬及防洪排澇分析,可為相關(guān)人工建設(shè)水網(wǎng)的防洪排澇研究提供借鑒和指導(dǎo)。

    1方法原理

    1.1排澇模數(shù)法

    本文研究區(qū)域各河流無(wú)實(shí)測(cè)流量資料,采用《安徽省淮北地區(qū)除澇水文計(jì)算辦法》排澇模數(shù)法[19]計(jì)算各區(qū)域流量。城南新區(qū)各流域分區(qū)面積都在50 km2以下,采用流域面積在50 km2以下的排澇模數(shù)公式,24 h凈雨24 h排出計(jì)算,見(jiàn)式(1):

    M =(1)

    考慮周邊地區(qū)按規(guī)劃建設(shè)后,用地類(lèi)型改變對(duì)地表徑流的影響,結(jié)合城市防洪排澇特性,城南水網(wǎng)按24 h凈雨10 h排出計(jì)算,見(jiàn)式(2):

    M =(2)

    式中:α為權(quán)重值,其值為0.84;R 為降雨量,mm;M為所求排澇模數(shù),m3/(s·km2)。

    排澇模數(shù)乘以區(qū)域面積,得到區(qū)域排澇流量:

    Q = FM(3)

    式中:F 為區(qū)域面積,km2;M 為排澇模數(shù),m3/(s· km2);Q 為所求排澇流量,m3/s。

    1.2 SWMM

    本文采用 SWMM模型[20]計(jì)算排澇流量,該模型是一個(gè)動(dòng)態(tài)降水徑流模擬模型,選用 GREEN-AMPT入滲模型計(jì)算降雨入滲,見(jiàn)式(4):

    fp = Ks [ ](4)

    式中:Ks 為飽和導(dǎo)水率,mm/h;ψs 為潤(rùn)濕前沿處的毛細(xì)作用力,mm;d 為表面沉積水深,mm;Ls 為飽和層深度,mm。

    其非線性水庫(kù)求解方程見(jiàn)式(5)、(6):

    dt = S dt = Si - Q(5)

    Q = L (h - hp )5/3S0(1)/2(6)

    式中:V=Sh為地表集水量,m3;h 為水深,m;t 為時(shí)間,s;S 為地表面積,m2;i為凈雨強(qiáng)度,mm/s;Q 為出流量,m3/s;L 為子流域?qū)挾?,m;n 為曼寧糙率系數(shù);hp 為地面蓄水深,mm;S0為子流域坡度。

    1.3 HEC-RAS二維模型

    采用 HEC-RAS 二維水動(dòng)力 SWE-ELW 模型[21],其質(zhì)量守恒及動(dòng)量守恒方程見(jiàn)式(7)—(9):

    式中:x、y 分別為笛卡爾直角坐標(biāo)系距離坐標(biāo),m;t為時(shí)間坐標(biāo),s;h 為水深,m;g 為重力加速度,9.8 m/ s2;zs 為水面高程,m;R 為水力半徑,m;fc 為科里奧利參數(shù),1/s;pa 為大氣壓強(qiáng),Pa;q 分別為源匯向量,m3/ s;u、v 分別為 x、y 方向的流速,m/s;vt,xx、vt,yy分別為 x、 y 方向的水平渦黏性系數(shù),m2/s;τb,x、τb,y 分別為 x、y 方向的底部剪應(yīng)力,N/m2;τs,x、τs,y 分別為 x、y 方向的表面風(fēng)應(yīng)力,N/m2。

    SWE-ELW 模型數(shù)值解算采用歐拉-拉格朗日法,見(jiàn)式(10)—(13):

    式中:Ak 為網(wǎng)格垂直面積,m2;AL 為左網(wǎng)格水平面積, m2;AR 為右網(wǎng)格水平面積;Vj為網(wǎng)格平均水流速度矢量,m/s;vt,k 為面渦黏性系數(shù),m2/s;nk為面法向單位向量;nij為網(wǎng)格質(zhì)心i和j 與相鄰網(wǎng)格 k 方向的單元向量,m;?xi,j 為網(wǎng)格質(zhì)心i和j 與相鄰網(wǎng)格 k 方向的距離,m。

    2城南水網(wǎng)

    2.1水系分布

    安徽省阜陽(yáng)市屬淮河流域,地處中國(guó)南方雨量豐沛和北方干旱少雨的過(guò)渡地帶,屬暖溫帶半濕潤(rùn)季風(fēng)氣候區(qū)。春末及夏季,西太平洋副熱帶高壓移近本流域,且比較穩(wěn)定,來(lái)自西南印度洋孟加拉灣和西太平洋的水汽,隨西南和東南季風(fēng)輸入本區(qū)域。阜陽(yáng)市河流均屬淮河水系,水系有外河和內(nèi)河之分,外河主要為潁河、泉河、茨河等6條大河,內(nèi)河為中清河、東清河、華橋溝等64條主干河道。

    城南新區(qū)位于阜陽(yáng)市潁西區(qū),為阜陽(yáng)的金融商務(wù)中心、行政辦公中心、旅游會(huì)展中心、文化體育中心、研發(fā)創(chuàng)意中心和商業(yè)居住中心。城南新區(qū)地面高程為28.5~31.5 m,地形高差較小,整體較平坦。區(qū)域總體地勢(shì)自西北向東南傾斜,北高南低、西高東低,地面平均坡降1/5000~1/10000,為典型平原地貌。根據(jù)《阜陽(yáng)市城市排水(雨水)防澇綜合規(guī)劃(2013—2030)》(阜陽(yáng)市城鄉(xiāng)規(guī)劃設(shè)計(jì)研究院,2014年4月),阜陽(yáng)市規(guī)劃區(qū)被劃分為19個(gè)雨水分區(qū)。

    其中,城南新區(qū)大部分屬于華橋溝分區(qū),極少部分屬于七漁河分區(qū),區(qū)域流域面積18.53 km2。城南新區(qū)澇水經(jīng)城南水網(wǎng)調(diào)蓄后,通過(guò)中清河、東清河外排澇水至下游華僑溝,由華橋溝泵站抽排入潁河,見(jiàn)圖1。

    城南水網(wǎng)為阜陽(yáng)市內(nèi)河水系,位于阜陽(yáng)市潁西區(qū)城南新區(qū)內(nèi),水系面積2.32 km2,占區(qū)域流域面積的12.67%。城南水網(wǎng)是城南新區(qū)重要基礎(chǔ)設(shè)施,負(fù)責(zé)保障區(qū)域排水、水生態(tài)和水環(huán)境。城南水網(wǎng)總布局為8河1灣,東西走向?yàn)榱锖?、五道河和保豐河3條河道,南北走向?yàn)槔衔髑搴?、南北河、西清河、中清河及東清河5條河道,水系中心為雙清灣湖。其中8條河道總長(zhǎng)度25.27 km,雙清灣面積26.10 hm2,見(jiàn)圖2。

    城南水網(wǎng)主要承接城南新區(qū)市政雨水管網(wǎng)排水,周邊用地規(guī)劃見(jiàn)圖3,市政雨水管網(wǎng)分布見(jiàn)圖4。城南水網(wǎng)各河流端口處均設(shè)置有節(jié)制閘,發(fā)生強(qiáng)降雨時(shí),按照調(diào)度安排關(guān)閉水系端口處全部節(jié)制閘,水網(wǎng)形成獨(dú)立調(diào)蓄水庫(kù),先蓄后排。水網(wǎng)排澇時(shí),老西清河為南向北流;南北河、西清河及中清河為北向南流;五道河及保豐河為西向東流;六里河在中清河以西為西向東流,在涵洞前為東向西流,在涵洞后為西向東流;東清河在五道河以北為南向北流,在五道河以南為北向南流。城南水網(wǎng)河流流向整體趨勢(shì)為北向南流,西向東流。

    2.2城南水網(wǎng)排澇方案

    根據(jù)歷史洪水資料,城南新區(qū)區(qū)域澇水與大河洪水幾乎同頻發(fā)生,故當(dāng)發(fā)生排澇標(biāo)準(zhǔn)澇水時(shí),外接大河水位均較高,無(wú)法自排入河,需要通過(guò)泵站進(jìn)行強(qiáng)排。

    根據(jù)《阜陽(yáng)市城市水系規(guī)劃(2013—2030)》(阜陽(yáng)市城鄉(xiāng)規(guī)劃設(shè)計(jì)研究院,2013年10月),城南水網(wǎng)排澇標(biāo)準(zhǔn)為30年一遇,排澇方案采用先蓄后排;汛前水位為26.50 m,發(fā)生強(qiáng)降雨時(shí)關(guān)閉水系內(nèi)所有水工建筑物,充分發(fā)揮河道、湖泊調(diào)蓄功能,將前期澇水蓄存在水系內(nèi);當(dāng)水位到達(dá)27.50 m(華橋溝泵站啟動(dòng)水位)時(shí),打開(kāi)中清河南閘、東清河南閘及東清河北閘,外排澇水至下游華僑溝,由華橋溝泵站抽排入潁河;水系汛期最高水位不高于28.30 m。

    3排澇流量

    阜陽(yáng)市年最大24 h暴雨均值為101.4 mm,CV 為0.472,CS/CV 為2,計(jì)算得到阜陽(yáng)市30年一遇24 h設(shè)計(jì)暴雨雨量為205 mm,按“安徽省84洪水計(jì)算辦法”中24 h 降雨時(shí)程分配成果進(jìn)行分配得到設(shè)計(jì)暴雨過(guò)程線,見(jiàn)圖5。

    根據(jù)圖4排水分區(qū)與管網(wǎng)圖劃分 SWMM匯水分區(qū)并量取面積,匯水分區(qū)面積為0.07~0.75 km2,見(jiàn)圖6;通過(guò)圖3量算加權(quán)得到各匯水分區(qū)綜合不透水面積比;SWMM模型選用 GREEN-AMPT入滲模型,項(xiàng)目區(qū)土壤為黏質(zhì)壤土,模型設(shè)置吸入水頭為210 mm,飽和導(dǎo)水率為1.16 mm/h,初始虧損百分比為0.15;硬化地表糙率采用0.015;非硬化地表采用0.15。

    SWMM 輸入暴雨過(guò)程與模型參數(shù)計(jì)算得到各匯入點(diǎn)排澇流量過(guò)程,見(jiàn)圖7。排澇模數(shù)法公式代入各匯水分區(qū)面積及30年一遇24 h設(shè)計(jì)暴雨雨量205 mm,計(jì)算得到各匯入點(diǎn)排澇流量。

    圖8a所示,各匯入點(diǎn) SWMM模型計(jì)算結(jié)果大于排澇模數(shù)法結(jié)果,流量差值為0.15~1.15 m3/s。

    根據(jù)《阜陽(yáng)市城市排水(雨水)防澇綜合規(guī)劃(2013—2030)》,與項(xiàng)目區(qū)匯水分區(qū)面積接近的排水泵站有儲(chǔ)臺(tái)泵站和河濱路泵站,儲(chǔ)臺(tái)泵站分區(qū)面積為0.43 km2,20年一遇規(guī)劃流量為2.00 m3/s,河濱路泵站分區(qū)面積為0.74 km2,20年一遇規(guī)劃流量為9.00 m3/s。將 SWMM模型計(jì)算結(jié)果、排澇模數(shù)法計(jì)算結(jié)果、儲(chǔ)臺(tái)泵站和河濱路泵站規(guī)劃成果換算為排澇模數(shù)進(jìn)行比較,見(jiàn)圖8b。排澇模數(shù)法排澇模數(shù)均為5.30 m3/(s·km2),接近儲(chǔ)臺(tái)泵站排澇模數(shù),與河濱路泵站排澇模數(shù)數(shù)值相差較大。SWMM 模型計(jì)算結(jié)果排澇模數(shù)介于儲(chǔ)臺(tái)泵站與河濱路泵站排澇模數(shù)之間,且排澇模數(shù)隨分區(qū)面積增大而減小。

    排澇模數(shù)法一般用于農(nóng)田等排澇計(jì)算,公式僅涉及到面積與排水時(shí)間;SWMM模型通過(guò)暴雨過(guò)程線,不透水面積比、糙率等參數(shù)體現(xiàn)開(kāi)發(fā)建設(shè)程度、表面糙率、透水面積等條件的影響,排澇模數(shù)與規(guī)劃泵站排澇模數(shù)接近,計(jì)算結(jié)果更為合理。本次各匯入點(diǎn)30年一遇設(shè)計(jì)排澇流量采用 SWMM計(jì)算成果。

    4二維模型

    4.1二維模型設(shè)置

    本文模型算例模擬范圍為城南水網(wǎng)內(nèi)六里河、五道河、保豐河、老西清河、南北河、西清河、中清河、東清河及雙清灣湖??紤]到水網(wǎng)水系縱橫交錯(cuò),采用 HEC-RAS SWE-ELW二維模型模擬。

    采用2016年《阜陽(yáng)市城南新區(qū)水系綜合治理》設(shè)計(jì)方案用于模型建模,生成1 m×1 m 的模型 DEM, DEM 完全體現(xiàn)設(shè)計(jì)圖紙,實(shí)現(xiàn)河流景觀斷面、河流交匯處平滑順接及雙清灣湖泊島嶼,實(shí)現(xiàn)全部地形細(xì)節(jié)保證 DEM 的高精度,見(jiàn)圖9。

    城南水網(wǎng)內(nèi)計(jì)算河段長(zhǎng)25.27 km,計(jì)算面積為11.44 km2。河流間地塊區(qū)域網(wǎng)格設(shè)置為200 m×200 m,河流及湖泊區(qū)域網(wǎng)格設(shè)置為10 m×10 m,共計(jì)35504個(gè)計(jì)算網(wǎng)格,見(jiàn)圖10a。

    模型上邊界為各匯水分區(qū)匯入點(diǎn),輸入圖7各匯入點(diǎn)排澇流量過(guò)程;根據(jù)排澇方案,排澇過(guò)程中,僅中清河南閘、東清河南閘及東清河北閘會(huì)開(kāi)啟排水,模型內(nèi)邊界只設(shè)置中清河南閘、東清河南閘及東清河北閘。其中,中清河南閘為2扇4 m×5 m 閘門(mén),閘底高程22.85 m,東清河南閘為2扇4 m×4 m 閘門(mén),閘底高程24.20 m,東清河北閘1扇4 m×4 m 閘門(mén),閘底高程24.20 m;模型下邊界為中清河南邊界、東清河南邊界及東清河北邊界,考慮最不利工況,均設(shè)置水位高程27.50 m,見(jiàn)圖10b。

    HEC-RAS 二維模型的糙率采用分區(qū)方式對(duì)不同區(qū)域賦值,河岸為0.08,河底為0.03,湖泊為0.04,綠島為0.10,綠地為0.04,建筑區(qū)為0.10,廣場(chǎng)1區(qū)為0.03,廣場(chǎng)2區(qū)為0.03,涵洞為0.015,河坡為0.05,景觀島嶼為0.10,見(jiàn)圖11。

    模型內(nèi)初始水位為26.50 m;模型內(nèi)邊界各閘門(mén)調(diào)度方案為閘前水位高于27.52 m 開(kāi)啟閘門(mén),閘前水位低于27.51 m關(guān)閉閘門(mén),閘門(mén)初始開(kāi)度為0,閘門(mén)啟閉速度為0.5 m/min,;計(jì)算時(shí)間為48 h,計(jì)算步長(zhǎng)為5 s,存儲(chǔ)步長(zhǎng)為5 min。

    4.2結(jié)果與分析

    圖12所示,城南水系最高水位為27.50~27.89 m;東清河北閘上游水位最高,水位高程為27.89 m;中清河南閘下游、東清河南閘下游及東清河北閘下游水位最低,水位高程為27.50 m;城南水網(wǎng)最高水位分布整體趨勢(shì)為北高南低,西高東低;離閘門(mén)越近水位越低。

    圖13所示,城南水系最大流速為0~1.96 m/s;東清河北閘流速最大,流速為1.96 m/s;老西清河南端、五道河西端、保豐河西端、南北河北端、南北河南端、西清河北端、西清河南端及中清河北端游流速最小,流速為0 m/s;城南水網(wǎng)最大流速分布整體趨勢(shì)為東大西小,;離閘門(mén)越近流速越大。

    圖14所示,城南水網(wǎng)排澇時(shí),老西清河為南向北流;南北河、西清河及中清河為北向南流;五道河及保豐河為西向東流;六里河在中清河以西為西向東流,在涵洞前為東向西流,在涵洞后為西向東流;東清河在五道河以北為南向北流,在五道河以南為北向南流。城南水網(wǎng)河流流向整體趨勢(shì)為北向南流,西向東流。

    圖15所示,城南水網(wǎng)排澇時(shí),各閘門(mén)閘前水位,閘后水位、啟動(dòng)時(shí)間及閘門(mén)開(kāi)度基本一致;中清河南閘最大過(guò)流流量最大,最大流量為35.82 m3/s,東清河北閘最大過(guò)流流量最下,最大流量為25.99 m3/ s;各閘門(mén)閘前水位、閘后水位、閘門(mén)開(kāi)度及閘門(mén)過(guò)流過(guò)程線存在一定的波動(dòng),東清河南閘波動(dòng)最大,中清河南閘波動(dòng)最小。

    對(duì)比圖12、13及圖15可知,各閘門(mén)過(guò)流流量與其河道底寬成正比,各閘門(mén)各過(guò)程線波動(dòng)程度與其河道底寬成反比。

    圖7各匯入點(diǎn)流量過(guò)程線相加得到城南水網(wǎng)入流曲線,圖15d各閘門(mén)過(guò)流過(guò)程線相加得到出流曲線。圖16所示,對(duì)比入流曲線與出流曲線的洪峰、洪量及峰值時(shí)間,城南水網(wǎng)削峰57.21 m3/s,錯(cuò)峰2.83 h,調(diào)蓄洪量169萬(wàn) m3,削峰率42.45%。

    5結(jié)論

    a)城南新區(qū)為無(wú)實(shí)測(cè)水文資料的新建城區(qū),本文采用排澇模數(shù)法及 SWMM 模型計(jì)算排澇流量。結(jié)果顯示,SWMM 模型中計(jì)算參數(shù)更多,排澇模數(shù)隨匯水分區(qū)面積增大而減小,且與規(guī)劃泵站排澇模數(shù)接近,計(jì)算結(jié)果更為合理。

    b)本文耦合 SWMM 分布式水文模型與 HEC- RAS 二維水動(dòng)力模型,同時(shí)發(fā)揮了水文、水動(dòng)力模型各自的優(yōu)勢(shì),能有效模擬閘門(mén)控制下城南水網(wǎng)排澇情景,可為相關(guān)城區(qū)規(guī)劃建設(shè)及水系閘門(mén)調(diào)度提供借鑒和參考。

    c)模擬結(jié)果顯示,城南水網(wǎng)排澇標(biāo)準(zhǔn)30年一遇工況下,城南水網(wǎng)水系汛期最高水位為27.89 m,小于28.30 m,滿足排澇方案要求;各閘門(mén)過(guò)流流量與其河道底寬成正比,各閘門(mén)各過(guò)程線波動(dòng)程度與其河道底寬成反比;城南水網(wǎng)排澇方案下,水系先蓄后排削峰57.21 m3/s,錯(cuò)峰2.83 h,調(diào)蓄洪量169萬(wàn) m3,削峰率42.45%,防洪排澇效益顯著。

    d)本文采用模擬方法具備一定的模擬效果及精度,可為相關(guān)人工建設(shè)水網(wǎng)的防洪排澇研究提供借鑒和指導(dǎo)。

    參考文獻(xiàn):

    [1]徐馳,董慶華,張宏雅,等.南京市高淳區(qū)城市水網(wǎng)防洪及雨洪資源利用模擬研究[J].中國(guó)農(nóng)村水利水電,2020(5):53-57,62.

    [2]徐馳,彭振陽(yáng),黃金鳳,等.多目標(biāo)生態(tài)水網(wǎng)構(gòu)建與評(píng)價(jià)方法研究[J].人民長(zhǎng)江,2022,53(3):79-86.

    [3]范立柱,樊博,李游坤.城鎮(zhèn)河網(wǎng)水系排澇優(yōu)化方案研究[J].水利發(fā)展研究,2023,23(2):65-69.

    [4]閆偉偉,何德舜,廖文俊.武漢市新洲區(qū)漲渡湖流域防洪排澇規(guī)劃方案擬定[J].中國(guó)農(nóng)村水利水電,2019(8):88-92.

    [5]袁玉,高玉琴,吳錫.基于 HEC_HMS 水文模型的秦淮河流域圩垸式防洪模式洪水模擬[J].三峽大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,37(5):34-39.

    [6]蔣建靈.杭州市區(qū)城西洪澇災(zāi)害防治方案研究[J].中國(guó)水利水電科學(xué)研究院學(xué)報(bào),2021,19(4):434-440.

    [7] SMITH M B. Comment on ′Analysis and modeling of flooding inurban drainage systems′[J]. Journal of Hydrology,2006,317(3/4):355-363.

    [8]XU D,LIU Y J. Research on the effect of rainfall flood regulation and control of wetland park based on SWMM model—a case study of wetland park in Yuanjia village , Qishan county, Shaanxi province[J]. Earth and Environmental Science,2018,121(5). DOI:10.1088/1755-1315/121/5/052014.

    [9]符銳.基于 SWMM模型的山前平原區(qū)城市排水防澇計(jì)算方法探究[J].人民珠江,2019,40(2):123-127.

    [10]吳沛霖,俞芳琴,王婷,等.基于 SWMM 的張家港市排水防澇風(fēng)險(xiǎn)評(píng)估[J].水文,2020,40(2):31-37.

    [11]劉瑤.基于模型分析構(gòu)建城市排水防澇方案:以深圳市福田河流域?yàn)槔跠].西安:西安建筑科技大學(xué),2015.

    [12]李尤,邸蘇闖,朱永華,等.基于多模型耦合的城市內(nèi)澇多視角分析[J].中國(guó)農(nóng)村水利水電,2019(5):82-90.

    [13]許婷婷,段夢(mèng),孫建升,等.雄安新區(qū)起步區(qū)排水防澇規(guī)律研究[J].給水排水,2021,47(12):54-59.

    [14]葉陳雷,徐宗學(xué).水工程調(diào)度與低影響開(kāi)發(fā)協(xié)同作用下典型城市片區(qū)洪澇過(guò)程模擬:以福州市為例[J].水利學(xué)報(bào),2022,53(7):833-844.

    [15]徐偉,劉茂,楊杰,等.基于 HEC-RAS 的淮河淮南段洪水漫頂風(fēng)險(xiǎn)分析[J].長(zhǎng)江科學(xué)院院報(bào),2011,28(7):13-18.

    [16]田景環(huán),張科磊,陳猛,等. HEC-RAS模型在洪水風(fēng)險(xiǎn)分析評(píng)估中的應(yīng)用研究[J].水電能源科學(xué),2012,30(4):23-25.

    [17]SURIYA S,MUDGAL B. Impact of urbanization on flooding:The Thirusoolamsub watershed: Acasestudy[J]. Journalof Hydrology,2012,412/413:210-219.

    [18]石振明,熊永峰,彭銘,等.堰塞湖潰壩快速定量風(fēng)險(xiǎn)評(píng)估方法:以2014年魯?shù)榈卣鹦纬傻募t石巖堰塞湖為例[J].水利學(xué)報(bào),2016,47(6):742-751.

    [19]欒巍,劉微.安徽省淮北地區(qū)除澇水文計(jì)算辦法初探[J].水利科技與經(jīng)濟(jì),2014,20(7):63-65.

    [20]ROSSMAN L A. Storm water management model user′s manual, version 5.0[R]. Cincinnati:National Risk Management Research Laboratory,Office of Research and Development,US Environmen? tal Protection Agency,2010.

    [21]US Army Corps of Engineers Hydrologic Engineering Center. HEC-RAS,River Analysis System Hydraulic Reference Manual [M]. Davis CA: US Army Corps of Engineers Hydrologic Engineering Center,2023.

    (責(zé)任編輯:程茜)

    猜你喜歡
    水網(wǎng)閘門(mén)
    水網(wǎng)結(jié)構(gòu)對(duì)多災(zāi)害情景的韌性響應(yīng)研究
    我終于會(huì)過(guò)地鐵閘門(mén)了
    水利部:到2035 年基本形成國(guó)家水網(wǎng)總體格局
    水網(wǎng)藻生長(zhǎng)及鈾對(duì)其生長(zhǎng)影響的模型構(gòu)建
    ●山西加快推進(jìn)縣域小水網(wǎng)配套建設(shè)
    弧形閘門(mén)斜支臂“扭轉(zhuǎn)角”研究
    山西省大水網(wǎng)工程建設(shè)如火如荼
    山西水利(2017年1期)2017-03-07 08:50:31
    把住醫(yī)?;鹚亻l門(mén)
    找準(zhǔn)入口,打開(kāi)思路的閘門(mén)
    山西召開(kāi)大水網(wǎng)決戰(zhàn)暨小水網(wǎng)啟動(dòng)工作推進(jìn)會(huì)
    山西水利(2016年9期)2016-04-08 09:42:52
    大香蕉久久网| 久久久色成人| 亚洲伊人久久精品综合| 永久免费av网站大全| 国产成人精品福利久久| 赤兔流量卡办理| 欧美极品一区二区三区四区| 久久久久性生活片| 亚洲美女搞黄在线观看| 视频中文字幕在线观看| 欧美少妇被猛烈插入视频| 亚洲天堂av无毛| 午夜福利网站1000一区二区三区| 国产成人精品福利久久| 尾随美女入室| 国产91av在线免费观看| 免费看av在线观看网站| 51国产日韩欧美| 热99国产精品久久久久久7| 综合色丁香网| 自拍欧美九色日韩亚洲蝌蚪91 | 国模一区二区三区四区视频| 狂野欧美激情性xxxx在线观看| 国产精品久久久久久精品电影| 3wmmmm亚洲av在线观看| 国产黄频视频在线观看| 91午夜精品亚洲一区二区三区| 一个人看的www免费观看视频| 日韩一区二区视频免费看| 国产精品久久久久久久久免| 国产老妇伦熟女老妇高清| 国产在视频线精品| 噜噜噜噜噜久久久久久91| 欧美精品国产亚洲| 久久久久精品久久久久真实原创| 国内少妇人妻偷人精品xxx网站| 成年版毛片免费区| 亚洲国产精品999| 国产色爽女视频免费观看| 韩国高清视频一区二区三区| 97精品久久久久久久久久精品| 麻豆成人午夜福利视频| 插阴视频在线观看视频| 狂野欧美激情性xxxx在线观看| 久久久久久国产a免费观看| 亚洲国产欧美在线一区| 国产精品一及| 少妇人妻久久综合中文| 精品国产三级普通话版| 美女主播在线视频| 国产一区二区三区av在线| 免费观看在线日韩| 黄色日韩在线| 人体艺术视频欧美日本| 在线a可以看的网站| 三级经典国产精品| 十八禁网站网址无遮挡 | 一级片'在线观看视频| 日韩大片免费观看网站| 久久精品久久久久久噜噜老黄| 2021天堂中文幕一二区在线观| 中文天堂在线官网| 国产免费福利视频在线观看| 日韩av不卡免费在线播放| 边亲边吃奶的免费视频| 日韩一本色道免费dvd| 亚洲精品日韩在线中文字幕| 国产黄色免费在线视频| 神马国产精品三级电影在线观看| 男人舔奶头视频| 91aial.com中文字幕在线观看| 日本与韩国留学比较| 亚洲真实伦在线观看| 欧美日韩视频精品一区| 亚洲精品日韩av片在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲精品视频女| 国内少妇人妻偷人精品xxx网站| 男人舔奶头视频| 日韩av不卡免费在线播放| 欧美xxⅹ黑人| 亚洲一区二区三区欧美精品 | 日本欧美国产在线视频| 人妻少妇偷人精品九色| 国产黄色视频一区二区在线观看| 最近的中文字幕免费完整| 免费观看的影片在线观看| 精品人妻熟女av久视频| 黄色视频在线播放观看不卡| 亚洲欧美中文字幕日韩二区| 亚洲国产精品成人久久小说| 国产人妻一区二区三区在| 性色av一级| 久久久久久久国产电影| 99热这里只有精品一区| 丝瓜视频免费看黄片| 欧美潮喷喷水| 夜夜爽夜夜爽视频| 赤兔流量卡办理| 免费看光身美女| 欧美日韩综合久久久久久| 精品国产露脸久久av麻豆| 国产在线男女| 亚洲丝袜综合中文字幕| 日韩电影二区| 99久久精品一区二区三区| 国产精品国产av在线观看| 国产成人freesex在线| 成人亚洲欧美一区二区av| 日韩av不卡免费在线播放| 欧美日韩在线观看h| 免费播放大片免费观看视频在线观看| www.色视频.com| 中文精品一卡2卡3卡4更新| 久久人人爽人人爽人人片va| 91久久精品国产一区二区三区| 亚洲av一区综合| 卡戴珊不雅视频在线播放| 亚洲va在线va天堂va国产| 熟女人妻精品中文字幕| 久热久热在线精品观看| 成年av动漫网址| av黄色大香蕉| av专区在线播放| 精品少妇黑人巨大在线播放| 国产成人精品福利久久| 最近中文字幕2019免费版| 麻豆久久精品国产亚洲av| 亚洲真实伦在线观看| 日韩精品有码人妻一区| 国产亚洲91精品色在线| 久久久精品欧美日韩精品| 好男人视频免费观看在线| 国产黄色视频一区二区在线观看| 亚洲精品456在线播放app| 亚洲不卡免费看| 国产淫语在线视频| 人妻一区二区av| 久久女婷五月综合色啪小说 | 久久99蜜桃精品久久| 视频区图区小说| 精品久久国产蜜桃| 激情五月婷婷亚洲| 亚洲,一卡二卡三卡| 99久久精品国产国产毛片| 建设人人有责人人尽责人人享有的 | 小蜜桃在线观看免费完整版高清| 日日摸夜夜添夜夜添av毛片| 成人亚洲精品一区在线观看 | 精品久久久久久久人妻蜜臀av| 在线播放无遮挡| 免费观看性生交大片5| videos熟女内射| 国产在线男女| 久久亚洲国产成人精品v| 欧美成人一区二区免费高清观看| 能在线免费看毛片的网站| 中文乱码字字幕精品一区二区三区| 亚洲人与动物交配视频| 久久久久久久大尺度免费视频| 日韩一区二区视频免费看| 精品酒店卫生间| 亚洲欧美精品自产自拍| 亚洲精品久久久久久婷婷小说| 精品久久久久久电影网| a级毛片免费高清观看在线播放| av卡一久久| 少妇熟女欧美另类| 在线观看人妻少妇| 国产探花极品一区二区| 亚洲电影在线观看av| 亚洲av日韩在线播放| 美女主播在线视频| 在线观看三级黄色| 久久99热这里只频精品6学生| 精品人妻熟女av久视频| 亚洲va在线va天堂va国产| 三级国产精品片| 国产综合懂色| 伦精品一区二区三区| 成人毛片60女人毛片免费| 男女下面进入的视频免费午夜| 精品久久国产蜜桃| 一二三四中文在线观看免费高清| 丝袜脚勾引网站| 欧美极品一区二区三区四区| 成人美女网站在线观看视频| 大香蕉久久网| 亚洲精品视频女| 免费不卡的大黄色大毛片视频在线观看| 欧美zozozo另类| 99热全是精品| 日韩视频在线欧美| 三级经典国产精品| 欧美精品人与动牲交sv欧美| 欧美日韩在线观看h| 久久久国产一区二区| av国产精品久久久久影院| 亚洲精品第二区| 日韩视频在线欧美| 免费大片黄手机在线观看| 午夜福利在线在线| 国产亚洲91精品色在线| 六月丁香七月| 日韩av免费高清视频| 亚洲高清免费不卡视频| 韩国av在线不卡| 中文在线观看免费www的网站| 国产精品99久久久久久久久| 中文字幕av成人在线电影| 色吧在线观看| 久久97久久精品| 久久精品久久久久久久性| 成人一区二区视频在线观看| 网址你懂的国产日韩在线| 国产极品天堂在线| 日韩av免费高清视频| 国产亚洲av嫩草精品影院| 九九爱精品视频在线观看| 免费av观看视频| av在线蜜桃| 中国三级夫妇交换| 在线观看免费高清a一片| 亚洲人与动物交配视频| 极品少妇高潮喷水抽搐| 伦理电影大哥的女人| 2022亚洲国产成人精品| 777米奇影视久久| 最后的刺客免费高清国语| 精品99又大又爽又粗少妇毛片| 日韩在线高清观看一区二区三区| 1000部很黄的大片| 精品国产三级普通话版| 亚洲av一区综合| a级毛色黄片| 赤兔流量卡办理| 永久免费av网站大全| 亚洲国产av新网站| 亚洲aⅴ乱码一区二区在线播放| 性色av一级| kizo精华| av网站免费在线观看视频| 18禁裸乳无遮挡动漫免费视频 | 人妻 亚洲 视频| 国产午夜福利久久久久久| 国产精品秋霞免费鲁丝片| 亚洲色图综合在线观看| 日本黄色片子视频| 成人午夜精彩视频在线观看| 免费av不卡在线播放| 亚洲精品日韩av片在线观看| 亚洲av电影在线观看一区二区三区 | 久久99热6这里只有精品| 男人添女人高潮全过程视频| 国产白丝娇喘喷水9色精品| 2021天堂中文幕一二区在线观| 国产亚洲午夜精品一区二区久久 | 亚洲经典国产精华液单| 亚洲精品日韩av片在线观看| 我要看日韩黄色一级片| 黄色视频在线播放观看不卡| 日韩视频在线欧美| 免费av不卡在线播放| 亚洲伊人久久精品综合| 一二三四中文在线观看免费高清| 国产69精品久久久久777片| 日韩在线高清观看一区二区三区| 精品人妻偷拍中文字幕| 有码 亚洲区| 最近2019中文字幕mv第一页| 亚洲丝袜综合中文字幕| 我的老师免费观看完整版| 欧美另类一区| 久久精品夜色国产| 亚洲人成网站在线观看播放| 欧美精品一区二区大全| 国产精品三级大全| 精品人妻熟女av久视频| 亚洲丝袜综合中文字幕| av在线蜜桃| 色哟哟·www| 午夜激情福利司机影院| 麻豆成人av视频| 哪个播放器可以免费观看大片| 观看美女的网站| 九九爱精品视频在线观看| 欧美三级亚洲精品| 又粗又硬又长又爽又黄的视频| 一级毛片我不卡| 涩涩av久久男人的天堂| 三级男女做爰猛烈吃奶摸视频| 久久久久精品久久久久真实原创| videossex国产| 久久久a久久爽久久v久久| 亚洲欧美一区二区三区黑人 | 免费播放大片免费观看视频在线观看| 国产男女内射视频| 2021少妇久久久久久久久久久| 日本一本二区三区精品| 亚洲欧美精品自产自拍| 午夜激情久久久久久久| 97热精品久久久久久| 97在线人人人人妻| 国产老妇女一区| 免费av毛片视频| 性色avwww在线观看| 国产亚洲av片在线观看秒播厂| av福利片在线观看| 国产精品一及| 国产精品久久久久久精品电影小说 | 欧美激情在线99| 国产精品蜜桃在线观看| 成人高潮视频无遮挡免费网站| 26uuu在线亚洲综合色| 欧美日韩视频高清一区二区三区二| 成人二区视频| 岛国毛片在线播放| 草草在线视频免费看| 日韩免费高清中文字幕av| 在线观看人妻少妇| 身体一侧抽搐| 有码 亚洲区| 一级黄片播放器| 麻豆成人av视频| 97人妻精品一区二区三区麻豆| 成人二区视频| 黄色一级大片看看| 亚洲国产精品专区欧美| 性色avwww在线观看| 成人综合一区亚洲| 91久久精品电影网| 成人特级av手机在线观看| 亚洲欧美日韩另类电影网站 | 久久精品久久久久久久性| 国产成人freesex在线| 久久久久精品久久久久真实原创| videossex国产| av在线老鸭窝| 91久久精品电影网| 国产成人福利小说| a级毛色黄片| 大香蕉久久网| 高清日韩中文字幕在线| 精品久久久精品久久久| 一级爰片在线观看| 亚洲精品aⅴ在线观看| 噜噜噜噜噜久久久久久91| 超碰97精品在线观看| 欧美日韩综合久久久久久| 一边亲一边摸免费视频| 禁无遮挡网站| 日韩在线高清观看一区二区三区| 日本与韩国留学比较| 色播亚洲综合网| 国产色爽女视频免费观看| 人妻夜夜爽99麻豆av| 我要看日韩黄色一级片| 极品教师在线视频| 国产精品.久久久| 国内揄拍国产精品人妻在线| 婷婷色综合www| 联通29元200g的流量卡| 黄片无遮挡物在线观看| 毛片一级片免费看久久久久| 国产成人a区在线观看| 91午夜精品亚洲一区二区三区| 国精品久久久久久国模美| 国产乱人偷精品视频| 交换朋友夫妻互换小说| 激情五月婷婷亚洲| 成人鲁丝片一二三区免费| 白带黄色成豆腐渣| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 男人狂女人下面高潮的视频| 成年版毛片免费区| 国产成人91sexporn| 最近2019中文字幕mv第一页| 日韩国内少妇激情av| 免费播放大片免费观看视频在线观看| 一区二区三区免费毛片| 国产日韩欧美在线精品| 2022亚洲国产成人精品| 欧美zozozo另类| 男人添女人高潮全过程视频| 日日啪夜夜撸| 美女内射精品一级片tv| 国产精品麻豆人妻色哟哟久久| freevideosex欧美| 久久ye,这里只有精品| av线在线观看网站| 欧美精品一区二区大全| 日本av手机在线免费观看| 97在线人人人人妻| 亚洲av国产av综合av卡| 一区二区三区四区激情视频| 日本一本二区三区精品| 日韩一本色道免费dvd| 欧美zozozo另类| 免费黄频网站在线观看国产| 国产黄色免费在线视频| 国产成人免费无遮挡视频| 真实男女啪啪啪动态图| 久久99精品国语久久久| 天天躁日日操中文字幕| 可以在线观看毛片的网站| 免费高清在线观看视频在线观看| 在线观看免费高清a一片| 黄色欧美视频在线观看| 久久99蜜桃精品久久| 亚洲欧美清纯卡通| 91狼人影院| 中国三级夫妇交换| 亚洲精品第二区| 久久女婷五月综合色啪小说 | 亚洲四区av| 精品少妇黑人巨大在线播放| 丰满人妻一区二区三区视频av| 美女高潮的动态| 日本免费在线观看一区| 听说在线观看完整版免费高清| 老师上课跳d突然被开到最大视频| 欧美亚洲 丝袜 人妻 在线| 亚洲综合精品二区| 另类亚洲欧美激情| 又大又黄又爽视频免费| 简卡轻食公司| 亚洲国产日韩一区二区| 色哟哟·www| 伊人久久国产一区二区| 汤姆久久久久久久影院中文字幕| 欧美日韩国产mv在线观看视频 | 在线观看一区二区三区激情| 久久热精品热| av在线天堂中文字幕| 亚洲av成人精品一区久久| 日本色播在线视频| 亚洲国产色片| 精品少妇黑人巨大在线播放| 久久人人爽人人爽人人片va| 女的被弄到高潮叫床怎么办| 国产成人精品久久久久久| 亚洲aⅴ乱码一区二区在线播放| 夜夜爽夜夜爽视频| 国产真实伦视频高清在线观看| 热99国产精品久久久久久7| 亚洲最大成人中文| 国产男女内射视频| 国产91av在线免费观看| 午夜福利高清视频| 亚洲国产最新在线播放| 人人妻人人爽人人添夜夜欢视频 | 成年人午夜在线观看视频| 五月开心婷婷网| 亚洲最大成人av| 国产精品麻豆人妻色哟哟久久| av在线观看视频网站免费| 男插女下体视频免费在线播放| 天堂网av新在线| 国产老妇女一区| 久久99精品国语久久久| 国产在视频线精品| 天天一区二区日本电影三级| 久久精品综合一区二区三区| av在线app专区| 日本与韩国留学比较| 国产精品一区www在线观看| 九草在线视频观看| 男女国产视频网站| 日日啪夜夜爽| 只有这里有精品99| 国产在视频线精品| 欧美人与善性xxx| 久久久久久久午夜电影| 亚洲天堂国产精品一区在线| 最新中文字幕久久久久| a级毛色黄片| 午夜福利高清视频| 亚洲精品国产av成人精品| 国产成人精品一,二区| 国产永久视频网站| 男插女下体视频免费在线播放| 搞女人的毛片| 男女国产视频网站| 欧美 日韩 精品 国产| 国产精品一区二区在线观看99| 在线观看美女被高潮喷水网站| 亚洲不卡免费看| 插逼视频在线观看| 免费观看性生交大片5| 联通29元200g的流量卡| 草草在线视频免费看| 欧美丝袜亚洲另类| 大陆偷拍与自拍| 丝袜美腿在线中文| 91狼人影院| av一本久久久久| 成人一区二区视频在线观看| 最近最新中文字幕免费大全7| 亚洲美女搞黄在线观看| 国产极品天堂在线| 国产成人aa在线观看| 两个人的视频大全免费| 一本久久精品| 国产白丝娇喘喷水9色精品| 亚洲av一区综合| 亚洲精品一二三| 有码 亚洲区| 联通29元200g的流量卡| 亚洲av免费在线观看| 新久久久久国产一级毛片| 视频区图区小说| 91精品国产九色| 少妇人妻一区二区三区视频| 麻豆精品久久久久久蜜桃| 亚洲精品乱码久久久久久按摩| 午夜福利视频1000在线观看| 欧美日韩综合久久久久久| 国产精品久久久久久av不卡| 精品酒店卫生间| 亚洲,一卡二卡三卡| 高清视频免费观看一区二区| 在线观看国产h片| 亚洲人成网站在线观看播放| 日韩av免费高清视频| 免费人成在线观看视频色| 各种免费的搞黄视频| 中文字幕人妻熟人妻熟丝袜美| 国产真实伦视频高清在线观看| 建设人人有责人人尽责人人享有的 | 国产精品熟女久久久久浪| 最近中文字幕2019免费版| 国产免费一区二区三区四区乱码| 成人高潮视频无遮挡免费网站| 免费大片18禁| 男女无遮挡免费网站观看| 777米奇影视久久| 黄色欧美视频在线观看| 国产成年人精品一区二区| 国产大屁股一区二区在线视频| 少妇人妻久久综合中文| 亚洲精品一区蜜桃| 少妇人妻久久综合中文| 亚洲av在线观看美女高潮| 草草在线视频免费看| 97在线人人人人妻| 亚洲伊人久久精品综合| 亚洲国产精品国产精品| 精品久久久久久电影网| 免费黄网站久久成人精品| 日韩伦理黄色片| 嘟嘟电影网在线观看| 少妇人妻一区二区三区视频| 亚洲精品456在线播放app| 免费观看a级毛片全部| 亚洲av不卡在线观看| 亚洲天堂av无毛| 晚上一个人看的免费电影| 国产一区二区在线观看日韩| 久久久久久久大尺度免费视频| 欧美bdsm另类| 中文字幕制服av| 最新中文字幕久久久久| 丝瓜视频免费看黄片| 美女国产视频在线观看| 伊人久久国产一区二区| 免费电影在线观看免费观看| 免费少妇av软件| 国产成人福利小说| 精品久久久久久久人妻蜜臀av| 国产高清国产精品国产三级 | 色5月婷婷丁香| 午夜免费鲁丝| 女人久久www免费人成看片| 国产成人a∨麻豆精品| 亚洲成人av在线免费| 日本一二三区视频观看| 在线看a的网站| 亚洲人与动物交配视频| av播播在线观看一区| 极品教师在线视频| 亚洲美女视频黄频| 老师上课跳d突然被开到最大视频| 各种免费的搞黄视频| 一级毛片aaaaaa免费看小| 欧美性猛交╳xxx乱大交人| 黄色配什么色好看| 亚洲怡红院男人天堂| 丰满少妇做爰视频| 亚洲国产高清在线一区二区三| 精品久久久精品久久久| 亚洲欧美一区二区三区黑人 | 亚洲人成网站在线播| 视频区图区小说| 老女人水多毛片| 九色成人免费人妻av| 国产午夜精品一二区理论片| av一本久久久久| 在线精品无人区一区二区三 | 欧美性感艳星| 久久久久久久久久成人| 国产av国产精品国产| 国产老妇女一区| 在线 av 中文字幕| 22中文网久久字幕| 欧美3d第一页| 精品久久久久久久久亚洲| 99久久人妻综合| 18禁裸乳无遮挡动漫免费视频 | 亚洲精品色激情综合| 久久久精品免费免费高清| 狂野欧美激情性xxxx在线观看| 一个人看视频在线观看www免费| 嫩草影院入口| 午夜福利在线观看免费完整高清在| 国产黄色视频一区二区在线观看|