• <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
    欧美+亚洲+日韩+国产| 午夜a级毛片| 99久久99久久久精品蜜桃| 超碰成人久久| 亚洲 欧美 日韩 在线 免费| 欧美一区二区精品小视频在线| 亚洲精品在线美女| 97碰自拍视频| 国内揄拍国产精品人妻在线| 99在线人妻在线中文字幕| 国产伦人伦偷精品视频| 一本久久中文字幕| 日本撒尿小便嘘嘘汇集6| 欧美极品一区二区三区四区| 人妻丰满熟妇av一区二区三区| www日本在线高清视频| 黑人巨大精品欧美一区二区mp4| 亚洲av成人不卡在线观看播放网| 亚洲成人国产一区在线观看| 欧美日韩亚洲综合一区二区三区_| 母亲3免费完整高清在线观看| 在线观看免费午夜福利视频| 久久久久性生活片| 91大片在线观看| 国产激情久久老熟女| 男女那种视频在线观看| 国产精品一区二区免费欧美| 成人精品一区二区免费| 国产亚洲精品一区二区www| 99国产精品一区二区三区| 国产午夜福利久久久久久| 亚洲国产高清在线一区二区三| 全区人妻精品视频| 动漫黄色视频在线观看| ponron亚洲| 国产精品久久久人人做人人爽| 90打野战视频偷拍视频| 国产亚洲精品久久久久5区| tocl精华| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲精品久久久久久毛片| 午夜福利在线在线| 国产私拍福利视频在线观看| 精品国产乱码久久久久久男人| 久久香蕉激情| 在线看三级毛片| www.熟女人妻精品国产| 黑人操中国人逼视频| 国产精品98久久久久久宅男小说| 午夜福利18| 床上黄色一级片| 久久亚洲真实| 日韩有码中文字幕| tocl精华| 最新在线观看一区二区三区| 久99久视频精品免费| 欧美日韩一级在线毛片| 18禁黄网站禁片午夜丰满| 国产精品亚洲av一区麻豆| 级片在线观看| 亚洲国产日韩欧美精品在线观看 | 大型黄色视频在线免费观看| 色噜噜av男人的天堂激情| 久久精品影院6| 两个人视频免费观看高清| 后天国语完整版免费观看| 色av中文字幕| 国产成人精品久久二区二区免费| 777久久人妻少妇嫩草av网站| 男人舔奶头视频| 国产精品,欧美在线| 国产高清视频在线观看网站| 这个男人来自地球电影免费观看| 免费观看精品视频网站| av天堂在线播放| 久久中文字幕人妻熟女| 国产精品永久免费网站| 麻豆久久精品国产亚洲av| 日韩三级视频一区二区三区| 久久久精品国产亚洲av高清涩受| 特级一级黄色大片| avwww免费| 国产精品 国内视频| 国产高清videossex| 日本成人三级电影网站| 男女床上黄色一级片免费看| 国产成人av教育| 色综合婷婷激情| 国产三级在线视频| 搡老熟女国产l中国老女人| 一级黄色大片毛片| 男女下面进入的视频免费午夜| 俄罗斯特黄特色一大片| 一夜夜www| 人人妻,人人澡人人爽秒播| 99国产极品粉嫩在线观看| 波多野结衣高清无吗| 国产精品,欧美在线| а√天堂www在线а√下载| 麻豆国产av国片精品| 成人三级做爰电影| 波多野结衣巨乳人妻| 午夜两性在线视频| 99re在线观看精品视频| 亚洲av电影在线进入| 午夜福利在线观看吧| 欧美中文日本在线观看视频| 一区二区三区国产精品乱码| 亚洲成人免费电影在线观看| 日韩欧美精品v在线| 久久中文看片网| 91av网站免费观看| 亚洲av第一区精品v没综合| 桃色一区二区三区在线观看| 久久草成人影院| 国产亚洲精品综合一区在线观看 | 久久香蕉激情| 国产aⅴ精品一区二区三区波| 亚洲午夜理论影院| 精品国产亚洲在线| 国产真人三级小视频在线观看| 丰满的人妻完整版| 国产黄色小视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国内久久婷婷六月综合欲色啪| 国产精品亚洲美女久久久| 亚洲成人久久性| 久久精品夜夜夜夜夜久久蜜豆 | 国产区一区二久久| 中文字幕久久专区| 在线观看免费视频日本深夜| 少妇人妻一区二区三区视频| 老汉色av国产亚洲站长工具| 久久久久九九精品影院| 国产激情久久老熟女| 小说图片视频综合网站| 一个人免费在线观看的高清视频| 国产成人啪精品午夜网站| 一区二区三区激情视频| 国内毛片毛片毛片毛片毛片| 日韩欧美在线二视频| 999久久久精品免费观看国产| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美在线二视频| 精品日产1卡2卡| 99国产综合亚洲精品| 校园春色视频在线观看| 国产精品综合久久久久久久免费| 亚洲av熟女| 亚洲天堂国产精品一区在线| 亚洲专区中文字幕在线| 一级片免费观看大全| 精品一区二区三区av网在线观看| 搡老岳熟女国产| 日韩高清综合在线| 国产伦一二天堂av在线观看| 十八禁网站免费在线| 国内精品久久久久久久电影| 一边摸一边做爽爽视频免费| 俄罗斯特黄特色一大片| 亚洲在线自拍视频| 午夜福利欧美成人| 国产成人一区二区三区免费视频网站| 老司机在亚洲福利影院| 久久久精品欧美日韩精品| 老司机福利观看| 午夜两性在线视频| 一夜夜www| 欧美久久黑人一区二区| 999久久久国产精品视频| 久久人妻福利社区极品人妻图片| 精品久久久久久久久久免费视频| 日韩成人在线观看一区二区三区| 国产熟女午夜一区二区三区| 国产黄片美女视频| 亚洲人成77777在线视频| 成人亚洲精品av一区二区| 动漫黄色视频在线观看| 久久人妻福利社区极品人妻图片| 制服人妻中文乱码| 日韩免费av在线播放| 欧美黑人欧美精品刺激| 高清毛片免费观看视频网站| 熟妇人妻久久中文字幕3abv| 日本五十路高清| 亚洲精品美女久久久久99蜜臀| 亚洲欧美激情综合另类| 黄色成人免费大全| 制服丝袜大香蕉在线| 午夜久久久久精精品| 好男人电影高清在线观看| 极品教师在线免费播放| 人妻丰满熟妇av一区二区三区| 亚洲自拍偷在线| 久久中文字幕人妻熟女| 啪啪无遮挡十八禁网站| 国产男靠女视频免费网站| 国产精品久久久人人做人人爽| 桃红色精品国产亚洲av| 精品久久久久久久人妻蜜臀av| 亚洲真实伦在线观看| 欧美在线一区亚洲| 欧美久久黑人一区二区| 日韩免费av在线播放| а√天堂www在线а√下载| 亚洲精品国产一区二区精华液| 国产99久久九九免费精品| 99国产精品一区二区三区| 黄色 视频免费看| 亚洲18禁久久av| 久久九九热精品免费| 久久午夜综合久久蜜桃| 精品不卡国产一区二区三区| avwww免费| 亚洲国产欧美一区二区综合| 在线播放国产精品三级| 国产视频一区二区在线看| 丰满人妻一区二区三区视频av | 精品国产乱子伦一区二区三区| 久久香蕉精品热| 国产蜜桃级精品一区二区三区| 国产成人影院久久av| 久久久久久久久免费视频了| 精品国产亚洲在线| av欧美777| 757午夜福利合集在线观看| 午夜福利高清视频| 国产蜜桃级精品一区二区三区| 天堂动漫精品| 99久久无色码亚洲精品果冻| 很黄的视频免费| 中文字幕久久专区| 18禁观看日本| 亚洲av日韩精品久久久久久密| 99热6这里只有精品| 亚洲av中文字字幕乱码综合| 国内少妇人妻偷人精品xxx网站 | 在线永久观看黄色视频| av超薄肉色丝袜交足视频| 日本三级黄在线观看| 国产精品久久久久久久电影 | 国产免费男女视频| 国产av在哪里看| 十八禁人妻一区二区| 欧美乱色亚洲激情| 淫秽高清视频在线观看| 怎么达到女性高潮| 免费电影在线观看免费观看| 欧美另类亚洲清纯唯美| 男人舔女人的私密视频| 最近在线观看免费完整版| 在线十欧美十亚洲十日本专区| 国产三级黄色录像| 最近最新中文字幕大全电影3| av天堂在线播放| 久久久久九九精品影院| 亚洲欧美精品综合一区二区三区| 听说在线观看完整版免费高清| 淫妇啪啪啪对白视频| 男人舔女人下体高潮全视频| 女人被狂操c到高潮| 三级国产精品欧美在线观看 | 日韩欧美在线二视频| 国产男靠女视频免费网站| 精品欧美国产一区二区三| 老司机在亚洲福利影院| 欧美zozozo另类| 色综合婷婷激情| 中文字幕人妻丝袜一区二区| 久久午夜综合久久蜜桃| 国产免费av片在线观看野外av| 两人在一起打扑克的视频| 成人高潮视频无遮挡免费网站| 亚洲精品色激情综合| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av| 亚洲精品美女久久av网站| 国产麻豆成人av免费视频| svipshipincom国产片| 亚洲五月天丁香| 又黄又粗又硬又大视频| 国产精品精品国产色婷婷| 免费观看精品视频网站| 免费av毛片视频| 狂野欧美激情性xxxx| 久久精品人妻少妇| 欧美日韩福利视频一区二区| 欧美日韩黄片免| 中国美女看黄片| 日本熟妇午夜| 一个人免费在线观看电影 | 国产精品久久久久久亚洲av鲁大| 日韩有码中文字幕| 老熟妇仑乱视频hdxx| 久久精品成人免费网站| 日韩欧美在线乱码| 亚洲av电影在线进入| 精品国内亚洲2022精品成人| 老汉色av国产亚洲站长工具| 成人18禁高潮啪啪吃奶动态图| 久久精品夜夜夜夜夜久久蜜豆 | 嫩草影视91久久| 午夜精品在线福利| 精品一区二区三区视频在线观看免费| 久久久久九九精品影院| 中文亚洲av片在线观看爽| 欧美在线一区亚洲| 成年版毛片免费区| 老鸭窝网址在线观看| 无限看片的www在线观看| 久久久久久久午夜电影| 一区二区三区国产精品乱码| 久久久国产精品麻豆| 国产99久久九九免费精品| 久久久久久久久久黄片| 久久热在线av| 美女黄网站色视频| av片东京热男人的天堂| av中文乱码字幕在线| www.www免费av| 欧美另类亚洲清纯唯美| 香蕉国产在线看| 丁香六月欧美| 亚洲 国产 在线| 一夜夜www| 男女那种视频在线观看| 亚洲精品一区av在线观看| 精品久久久久久久末码| 真人做人爱边吃奶动态| 五月伊人婷婷丁香| 欧美色视频一区免费| 日本熟妇午夜| 国产精品一区二区三区四区久久| 日本撒尿小便嘘嘘汇集6| 久久国产精品影院| 露出奶头的视频| 久久伊人香网站| 亚洲 国产 在线| 亚洲精华国产精华精| 看黄色毛片网站| 国产精品一及| 国产69精品久久久久777片 | 午夜久久久久精精品| 精品日产1卡2卡| 一二三四在线观看免费中文在| 国产真人三级小视频在线观看| 日韩高清综合在线| 夜夜躁狠狠躁天天躁| 欧美av亚洲av综合av国产av| 俺也久久电影网| 久久久久免费精品人妻一区二区| 国产精品久久久人人做人人爽| 99久久国产精品久久久| 777久久人妻少妇嫩草av网站| 视频区欧美日本亚洲| 国产激情久久老熟女| 伊人久久大香线蕉亚洲五| 日韩av在线大香蕉| 精品国产乱子伦一区二区三区| 97人妻精品一区二区三区麻豆| 欧美色视频一区免费| 国内精品一区二区在线观看| 日韩三级视频一区二区三区| 亚洲avbb在线观看| 亚洲成人久久性| 亚洲成av人片免费观看| 老熟妇仑乱视频hdxx| 亚洲人成电影免费在线| 神马国产精品三级电影在线观看 | 每晚都被弄得嗷嗷叫到高潮| 黄色片一级片一级黄色片| 国产一级毛片七仙女欲春2| 欧美在线一区亚洲| 国产成人av激情在线播放| 午夜a级毛片| 十八禁人妻一区二区| 日韩高清综合在线| 999久久久精品免费观看国产| 制服诱惑二区| 成人av一区二区三区在线看| 久久这里只有精品19| 久久久久亚洲av毛片大全| 精品一区二区三区四区五区乱码| 亚洲成人精品中文字幕电影| 人妻久久中文字幕网| 久久精品夜夜夜夜夜久久蜜豆 | 日韩成人在线观看一区二区三区| 成人亚洲精品av一区二区| 久久精品夜夜夜夜夜久久蜜豆 | 18禁黄网站禁片免费观看直播| 中出人妻视频一区二区| 欧美黑人精品巨大| 最新美女视频免费是黄的| 99精品欧美一区二区三区四区| 亚洲欧美日韩高清在线视频| 久久精品综合一区二区三区| 香蕉久久夜色| 两个人视频免费观看高清| 91在线观看av| 久久久久久久久久黄片| 欧美成人免费av一区二区三区| 欧美精品啪啪一区二区三区| 亚洲片人在线观看| 日本 欧美在线| 91九色精品人成在线观看| 久久99热这里只有精品18| 亚洲男人的天堂狠狠| 巨乳人妻的诱惑在线观看| 久久久久性生活片| 国产一区二区在线观看日韩 | 人成视频在线观看免费观看| 在线观看66精品国产| 午夜福利视频1000在线观看| 999久久久国产精品视频| 老司机午夜十八禁免费视频| 国产成人系列免费观看| 高清在线国产一区| 天堂√8在线中文| 中文字幕高清在线视频| av在线天堂中文字幕| 在线十欧美十亚洲十日本专区| 亚洲黑人精品在线| 亚洲成人免费电影在线观看| 女警被强在线播放| 亚洲国产精品sss在线观看| 美女扒开内裤让男人捅视频| 亚洲av美国av| 成人精品一区二区免费| 狠狠狠狠99中文字幕| 777久久人妻少妇嫩草av网站| a在线观看视频网站| 我的老师免费观看完整版| 无人区码免费观看不卡| 五月伊人婷婷丁香| 久久国产精品人妻蜜桃| 哪里可以看免费的av片| 国产精品亚洲美女久久久| 69av精品久久久久久| 久久香蕉精品热| 别揉我奶头~嗯~啊~动态视频| 久久欧美精品欧美久久欧美| 欧美成人免费av一区二区三区| 国产精品久久久人人做人人爽| 精品国内亚洲2022精品成人| www.999成人在线观看| 黄色毛片三级朝国网站| 国产精品精品国产色婷婷| 亚洲最大成人中文| 欧美成人一区二区免费高清观看 | 在线免费观看的www视频| 特大巨黑吊av在线直播| 国产男靠女视频免费网站| 国产成人av教育| 一进一出抽搐gif免费好疼| 欧美人与性动交α欧美精品济南到| 伊人久久大香线蕉亚洲五| 女人被狂操c到高潮| 两个人免费观看高清视频| 50天的宝宝边吃奶边哭怎么回事| 听说在线观看完整版免费高清| 国内精品一区二区在线观看| 毛片女人毛片| 国产精品久久久久久亚洲av鲁大| 欧美日韩亚洲国产一区二区在线观看| 欧美黑人欧美精品刺激| 久久久久免费精品人妻一区二区| 一级片免费观看大全| 身体一侧抽搐| 国产亚洲欧美98| 男女做爰动态图高潮gif福利片| 久久人妻av系列| 九九热线精品视视频播放| 99国产精品一区二区三区| 又大又爽又粗| 好男人电影高清在线观看| 免费看日本二区| 国产私拍福利视频在线观看| 色综合婷婷激情| cao死你这个sao货| 麻豆久久精品国产亚洲av| 国产亚洲欧美98| 黄色女人牲交| 真人做人爱边吃奶动态| 亚洲精品中文字幕一二三四区| 欧美性长视频在线观看| 日本免费a在线| 久久精品亚洲精品国产色婷小说| 露出奶头的视频| 丁香六月欧美| 床上黄色一级片| 精品国内亚洲2022精品成人| 波多野结衣高清作品| 精品久久久久久久久久免费视频| 黄色女人牲交| svipshipincom国产片| 久久婷婷人人爽人人干人人爱| 午夜成年电影在线免费观看| 色噜噜av男人的天堂激情| 特大巨黑吊av在线直播| 国产一区在线观看成人免费| 婷婷丁香在线五月| 日韩 欧美 亚洲 中文字幕| 精品欧美一区二区三区在线| 国产在线观看jvid| 可以免费在线观看a视频的电影网站| 亚洲美女黄片视频| 免费在线观看亚洲国产| 嫩草影视91久久| 国产麻豆成人av免费视频| 亚洲真实伦在线观看| √禁漫天堂资源中文www| 99在线人妻在线中文字幕| 天天一区二区日本电影三级| 精品欧美国产一区二区三| 嫩草影视91久久| 久久久久性生活片| 日韩成人在线观看一区二区三区| 成人18禁高潮啪啪吃奶动态图| 国产一区二区激情短视频| 又黄又爽又免费观看的视频| 一级作爱视频免费观看| 欧美在线黄色| www.熟女人妻精品国产| 亚洲,欧美精品.| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩乱码在线| 黑人巨大精品欧美一区二区mp4| 久久久久九九精品影院| 麻豆久久精品国产亚洲av| 国产日本99.免费观看| 亚洲真实伦在线观看| 久久精品国产99精品国产亚洲性色| 美女免费视频网站| 校园春色视频在线观看| 一区福利在线观看| √禁漫天堂资源中文www| 2021天堂中文幕一二区在线观| 国产97色在线日韩免费| 国产三级在线视频| 很黄的视频免费| 天天躁狠狠躁夜夜躁狠狠躁| 国产v大片淫在线免费观看| 国产一区二区在线观看日韩 | 我的老师免费观看完整版| 一本一本综合久久| 一级黄色大片毛片| 禁无遮挡网站| 日韩欧美 国产精品| 亚洲人成伊人成综合网2020| 日韩成人在线观看一区二区三区| 精品午夜福利视频在线观看一区| 亚洲片人在线观看| 久久午夜综合久久蜜桃| 一二三四社区在线视频社区8| 亚洲av中文字字幕乱码综合| 国产精品精品国产色婷婷| 欧美 亚洲 国产 日韩一| 日本a在线网址| 婷婷精品国产亚洲av| 黄片大片在线免费观看| 国产高清视频在线播放一区| 美女 人体艺术 gogo| 男人舔奶头视频| www国产在线视频色| 男人舔女人下体高潮全视频| 久久久久精品国产欧美久久久| 99精品欧美一区二区三区四区| 可以在线观看毛片的网站| 露出奶头的视频| 国产精品1区2区在线观看.| 亚洲欧美激情综合另类| 国产精品,欧美在线| 免费在线观看视频国产中文字幕亚洲| 毛片女人毛片| 色噜噜av男人的天堂激情| 精品高清国产在线一区| 俺也久久电影网| 婷婷精品国产亚洲av在线| 美女扒开内裤让男人捅视频| 免费在线观看完整版高清| 亚洲 欧美一区二区三区| 欧美日韩亚洲综合一区二区三区_| 亚洲色图av天堂| 国产av又大| 女人高潮潮喷娇喘18禁视频| svipshipincom国产片| 嫩草影院精品99| 一边摸一边做爽爽视频免费| 欧美黄色淫秽网站| 19禁男女啪啪无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 97超级碰碰碰精品色视频在线观看| 99精品欧美一区二区三区四区| 激情在线观看视频在线高清| 中文字幕精品亚洲无线码一区| 香蕉av资源在线| 欧美高清成人免费视频www| 高清毛片免费观看视频网站| 搡老岳熟女国产| 亚洲中文日韩欧美视频| 成人手机av| 国内揄拍国产精品人妻在线| 亚洲中文日韩欧美视频| 国产精华一区二区三区| 亚洲国产看品久久| 中文字幕av在线有码专区| 久久久久国内视频| 亚洲国产精品sss在线观看| 国产激情偷乱视频一区二区| 免费看十八禁软件| av国产免费在线观看| 99国产精品99久久久久| 国产亚洲av嫩草精品影院| 国产又色又爽无遮挡免费看|