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

    城市內(nèi)澇風(fēng)險分析模型系統(tǒng)開發(fā)與應(yīng)用

    2018-05-14 11:47:06周玉文婁富豪劉子龍楊偉明劉原吳獻平王中昌
    河北科技大學(xué)學(xué)報 2018年1期
    關(guān)鍵詞:徑流系數(shù)

    周玉文 婁富豪 劉子龍 楊偉明 劉原 吳獻平 王中昌

    摘要:為了全面客觀評估城市排水管網(wǎng)的排水能力、預(yù)防城市內(nèi)澇災(zāi)害的發(fā)生,基于EPA SWMM與ArcEngine技術(shù),結(jié)合等流時線法與水量平衡原理,開發(fā)了城市內(nèi)澇風(fēng)險分析模型系統(tǒng)(USRAMS)。闡述了USRAMS的原理與功能,說明了利用USRAMS進行管網(wǎng)排水能力評估與內(nèi)澇風(fēng)險分析的方法,并驗證了USRAMS的適用性;對滄州市進行了現(xiàn)狀管網(wǎng)排水能力評估和城市內(nèi)澇風(fēng)險分析,以主題圖的形式準(zhǔn)確直觀地表達了瓶頸管段位置以及城市內(nèi)澇風(fēng)險的分布情況。實例證明,USRAMS地表產(chǎn)匯流模塊只需徑流系數(shù)與匯水時間2個參數(shù),較好地適應(yīng)了中國當(dāng)前情況,二維地表淹水模擬采用水量平衡原理,計算穩(wěn)定、快速,其模型結(jié)果可為城市內(nèi)澇防治與雨水管網(wǎng)系統(tǒng)的規(guī)劃改造提供一定參考。

    關(guān)鍵詞:城市給水排水工程; USRAMS;等流時線;徑流系數(shù);排水能力評估;內(nèi)澇風(fēng)險分析

    中圖分類號:TU992文獻標(biāo)志碼:A

    Development and application of urban stormwater

    risk analysis modeling system

    ZHOU Yuwen1,2, LOU Fuhao1,2, LIU Zilong1,3, YANG Weiming1,2, LIU Yuan1,2,

    WU Xianping1,2, WANG Zhongchang4

    (1.College of Architecture and Civil Engineering, Beijing University of Technology, Beijing 100124, China; 2.Key Laboratory of Beijing for Water Quality Science and Water Environment Recovery Engineering, Beijing 100124, China; 3.Beijing Municipal Institute of City Planning & Design, Beijing 100037, China; 4.Cangzhou Urban Planning & Design Institute, Cangzhou, Hebei 061000, China)

    Abstract:In order to evaluate the pipeline drainage capacity and urban flood risk objectively, urban storm water risk analysis modeling system (USRAMS) is developed based on EPA SWMM and ArcEngine combined with the principle of equal flow time-line method and water balance. The principle and function of USRAMS and the method of assessing pipeline drainage capacity and urban flood risk with USRAMS are described, and the suitability of USRAMS is verified. Taking Cangzhou as an example, USRAMS is used to assess the pipeline drainage capacity and urban flood risk, and the position of bottleneck pipe and the waterlogging risk distribution are expressed in thematic map accurately and directly. The facts show that USRAMS only needs two parameters, namely the runoff coefficient and the catchment time, and adapts to the current situation in china better; two-dimensional earth surface waterlogging simulation adopts hydrologic budget theory, so its calculation is stable and fast. The results of USRAMS may provide reference for urban waterlogging prevention and control and the planning and reforming of stormwater pipe network system.

    Keywords:urban water supply and drainage engineering; USRAMS; isochronous stream line; runoff coefficient; drainage capacity assessment; local flooding risk analysis

    中國城市內(nèi)澇災(zāi)害較嚴重[1-3]。為應(yīng)對城市內(nèi)澇災(zāi)害,2013年國務(wù)院辦公廳發(fā)布了《關(guān)于做好城市排水防澇設(shè)施建設(shè)工作的通知》,住建部編制了《城市排水(雨水)防澇綜合規(guī)劃編制大綱》,提出了使用水力模型對城市管網(wǎng)排水能力進行評估的要求[4-5]。2014版《室外排水設(shè)計規(guī)范》大大提高了雨水管渠設(shè)計標(biāo)準(zhǔn),明確提出內(nèi)澇防治標(biāo)準(zhǔn),并要求匯水面積超過2 km2時采用數(shù)學(xué)模型法計算設(shè)計流量[6]。2013年以來各地陸續(xù)開展城市排水防澇規(guī)劃修編工作,而評估雨水排水管網(wǎng)排水能力和城市淹水風(fēng)險是做好規(guī)劃的基本前提。本文采用城市內(nèi)澇風(fēng)險分析模型系統(tǒng)(USRAMS),對滄州市現(xiàn)狀雨水管網(wǎng)排水能力及城市內(nèi)澇風(fēng)險進行評估分析。

    1城市內(nèi)澇風(fēng)險分析模型系統(tǒng)開發(fā)情況

    城市內(nèi)澇風(fēng)險分析模型系統(tǒng)(urban stormwater risk analysis modeling system)(簡稱USRAMS)是北京工業(yè)大學(xué)給排水系統(tǒng)研究室開發(fā)的[7]。USRAMS主要由一維管網(wǎng)模型模擬計算部分與二維地表淹水模擬計算部分組成,基于ArcGIS提供的嵌入式組件庫ArcEngine進行開發(fā)。其地圖數(shù)據(jù)管理平臺為ArcGIS,基礎(chǔ)層數(shù)據(jù)庫為MicrosoftAccess[8]。USRAMS采用目前國內(nèi)常用的SWMM作為一維管網(wǎng)模型水力計算的引擎,添加了等流時線模塊;二維地表漫流采用自主研發(fā)的內(nèi)澇淹水計算引擎,基于路徑追蹤[9]以及水量平衡原理進行二維地表淹水模擬分析。軟件系統(tǒng)基本構(gòu)架如圖1所示。USRAMS借助ArcGIS進行管網(wǎng)模型基礎(chǔ)數(shù)據(jù)的處理。

    2基本功能

    USRAMS的基本功能包括:基礎(chǔ)數(shù)據(jù)編輯與查看(降雨數(shù)據(jù)、管網(wǎng)數(shù)據(jù)、匯水區(qū)數(shù)據(jù)),數(shù)據(jù)檢查,上下游追蹤,管網(wǎng)縱剖面圖查看,主題圖繪制,數(shù)據(jù)導(dǎo)出及圖形的縮放、移動、查詢等。

    USRAMS主要利用ArcGIS與CADTableConvert完成模型數(shù)據(jù)的處理工作,也可對模型數(shù)據(jù)進行編輯與查看。軟件提供了降雨編輯器、入流編輯器、調(diào)蓄池容量編輯器、水泵性能曲線編輯器等;為保證模型數(shù)據(jù)質(zhì)量,USRAMS提供了多種數(shù)據(jù)檢查功能,對管網(wǎng)拓撲結(jié)構(gòu)關(guān)系、模型數(shù)據(jù)完整性、模型參數(shù)合理性進行檢查,包括孤立點檢查、管道逆坡檢查、管道錯位檢查、節(jié)點高程檢查、字段缺失檢查等 。

    上游追蹤功能可快速選擇與指定節(jié)點相連的所有上游管道并以高亮形式顯示。下游追蹤功能可快速選擇指定節(jié)點到最下游節(jié)點之間的管道,并以高亮形式顯示。通過上下游分析可以檢查管網(wǎng)連通性及查看管道路徑;通過縱剖面圖繪制功能可查看管道縱剖面圖,了解管道鋪設(shè)情況。USRAMS可繪制多種主題圖,包括管道充滿度、管道超載能力、管段瓶頸狀況、節(jié)點洪水量、流向、淹水分布圖。軟件還可以根據(jù)用戶需要自定義漸變主題圖。

    3模擬功能

    3.1 設(shè)計降雨生成

    降雨過程線是雨水管網(wǎng)水力模型的輸入條件,它有降雨量、降雨歷時、重現(xiàn)期、降雨峰值以及雨峰位置5個參數(shù)[10]。常見雨型主要有SCS雨型[11]、Keifer & Chu[12]雨型(芝加哥雨型)、均勻雨型、Huff[13]雨型、Yen & Chow[14]雨型(三角雨型)和24 h設(shè)計暴雨過程線[15-16]等,國內(nèi)多采用120 min或180 min的芝加哥雨型。USRAMS提供芝加哥雨型生成器(如圖2所示),還可以輸入實測降雨過程。

    3.2地表產(chǎn)匯流模型

    EPA SWMM采用前損后損法與非線性水庫模型進行地表產(chǎn)流和匯流過程模擬。此外,等流時線法、瞬時單位線法也可進行地表產(chǎn)匯流計算[17]。中國在進行雨水管網(wǎng)設(shè)計時多采用推理公式法,其主要設(shè)計參數(shù)為設(shè)計重現(xiàn)期、降雨歷時、匯水面積和徑流系數(shù),與之相對應(yīng)的計算方法為等流時線法[7]。

    USRAMS在SWMM計算引擎的基礎(chǔ)上添加了等流時線模塊,其降雨-徑流關(guān)系可用卷積分過程表達[18],流量過程線為

    式中:Q(ti)為ti時刻流量;F為匯水區(qū)面積;tc為地面匯水時間,取值范圍為5~15 min;qti-j為ti-j時刻地表凈雨強度,當(dāng)i-j≤0時,取ti-j=0;α為單位換算系數(shù),當(dāng)Q,q,F(xiàn)和t分別以L/s,mm/min,m2和min為單位計算時,α=60。

    3.3一維管網(wǎng)模型

    對于管網(wǎng)匯流,USRAMS采用與SWMM相同的計算方法,包括恒定流、運動波、動力波方法。USRAMS根據(jù)管道充滿度繪制結(jié)果主題圖,并據(jù)此評估管網(wǎng)排水能力,步驟如圖3所示。

    1)收集基礎(chǔ)資料,包括管道、檢查井、水泵、蓄水池、降雨等。

    2)利用ArcGIS構(gòu)建管網(wǎng)拓撲結(jié)構(gòu)與模型數(shù)據(jù)庫。

    3)利用芝加哥雨型生成器構(gòu)建降雨過程線,進行管網(wǎng)排水能力評估時一般選用120 min或180 min的芝加哥雨型。

    4)利用動力波對一維管網(wǎng)模型進行模擬計算。

    5)根據(jù)計算結(jié)果進行管道排水能力評估,主要指標(biāo)有節(jié)點溢流情況、管道最大充滿度、管段超載系數(shù)等。

    6)根據(jù)管道水力狀態(tài)分析產(chǎn)生管道超載(管道運行水面高于管頂)的原因是自身過流能力不足還是下游水位頂托。當(dāng)管段發(fā)生超載時,管段存在2種水力狀態(tài)如圖4所示。由水力學(xué)理論可知:當(dāng)超載管段屬于狀態(tài)a)時,管道實際過流能力小于管道滿流設(shè)計流量,發(fā)生管道超載的原因為下游水位頂托;當(dāng)超載管段屬于狀態(tài)b)時,管道實際過流能力大于管道滿流設(shè)計流量,發(fā)生管道超載的原因為自身過流能力不足。對于受下游水位頂托影響的管段,首先考慮改善下游水力狀態(tài);對于自身過流能力不足的管段,首先考慮增大設(shè)計斷面或坡度。因此,將水力狀態(tài)屬于狀態(tài)b)的超載管段定義為瓶頸管段。

    3.4二維地表淹水模擬計算

    USRAMS利用城市數(shù)字地形模型(DEM),基于路徑追蹤算法劃分排水流域,根據(jù)水量平衡原理進行二維地表淹水模擬計算,得到積水深度、積水量、積水時間等結(jié)果[7],將最大淹水結(jié)果、淹水時間與風(fēng)險等級以.asc格式文件輸出。風(fēng)險等級根據(jù)淹水深度與淹水時間按表1進行劃分。淹水深度與淹水時間的閾值可根據(jù)需要設(shè)定。

    輸出的.asc格式文件可利用ArcGIS進行渲染,制作最大淹水深度、淹水時間與風(fēng)險等級主題圖。城市內(nèi)澇風(fēng)險評估時其主要步驟如圖5所示。

    1)構(gòu)建降雨模型、管網(wǎng)模型、數(shù)字高程模型。暴雨可采用當(dāng)?shù)匕l(fā)生過的較有代表性的大暴雨歷程,或根據(jù)當(dāng)?shù)氐牡湫陀晷屯ㄟ^同頻率放大或同倍比法求得[19]。

    2)搭建USRAMS內(nèi)澇淹水模型。

    3)運行模型得到最大淹水深度、淹水時間與風(fēng)險等級,并輸出為.asc格式文件。

    4)制作淹水深度、淹水時間、風(fēng)險等級圖。

    4實例研究

    4.1研究區(qū)域概況

    研究區(qū)域為河北省滄州市,京杭大運河由南至北穿過市區(qū),將滄州市分為運東和運西兩大排水區(qū)[20]。滄州市現(xiàn)狀排水管網(wǎng)長307.5 km,匯水區(qū)面積為78.6 km2,設(shè)計重現(xiàn)期大部分在0.5 a以下。城區(qū)地勢西高東低,地面平坦。城區(qū)綜合徑流系數(shù)為0.6,地表匯水時間為15 min。滄州市地形分布如圖6所示。

    4.2現(xiàn)狀管網(wǎng)排水能力分析

    1)收集并輸入降雨、管網(wǎng)資料,構(gòu)建一維管網(wǎng)模型,對數(shù)據(jù)進行檢查。

    2)結(jié)合滄州市暴雨強度公式(見式(2))生成芝加哥雨型。

    式(2)中,取雨峰系數(shù)r=0.5,降雨歷時取120 min,分別生成重現(xiàn)期為0.5,1,2,3 a的芝加哥雨型,其中重現(xiàn)期為3 a的芝加哥雨型如圖7所示。

    3)分別進行重現(xiàn)期0.5,1,2,3 a情景下的降雨排水模擬計算,地表產(chǎn)匯流模型選擇等流時線法,管網(wǎng)匯流模型選擇動力波法,模擬時間取2 h,時間步長取1 s。根據(jù)模型結(jié)果得到管網(wǎng)最大充滿度,并繪制最大充滿度結(jié)果主題圖(如圖8所示)。

    由圖8可知,隨著重現(xiàn)期的提高,充滿度為1的管段數(shù)量增加。統(tǒng)計重現(xiàn)期為<0.5,0.5~1,1~2,2~3,≥3 a管道未充滿的長度,見表2。

    由表2和圖9可知,大部分管道的排水標(biāo)準(zhǔn)都在0.5 a以下,且主要分布在老城區(qū);運西新區(qū)的管道標(biāo)準(zhǔn)普遍高于老城區(qū),即滄州市老城區(qū)管道排水重現(xiàn)期多在0.5 a以下。滄州市規(guī)劃管網(wǎng)設(shè)計標(biāo)準(zhǔn)要提高到三年一遇,則現(xiàn)狀管道中有92.8%未達到標(biāo)準(zhǔn)。對于管道產(chǎn)生超載的原因需進一步分析。

    4.3瓶頸管段分析

    正如第3.3小節(jié)中所述,管道發(fā)生超載的原因有2種,分別為管道自身過流能力不足與受下游水位頂托影響。根據(jù)模型結(jié)果,USRAMS將管道分為3種類型,分別為未超載、瓶頸管段、受頂托管段。其中,瓶頸管段定義為管道超載時,管道最大水力坡度大于等于管道坡度;受頂托管段定義為管道超載時,管道最大水力坡度小于管道坡度。以1 a重現(xiàn)期模型情景為例,其瓶頸管段分布情況如圖10所示。

    通過USRAMS明確了管道超載的原因,并通過主題圖展示了瓶頸管段與受頂托管段的分布情況,為進一步的管網(wǎng)改造與規(guī)劃提供參考。

    4.4現(xiàn)狀內(nèi)澇風(fēng)險評估

    USRAMS城市內(nèi)澇模型的DEM的精度為10 m×10 m。選取一場24 h典型實測降雨進行風(fēng)險評估,降雨記錄時間間隔為5 min,總雨量為196.4 mm,降雨過程線如圖11所示。

    產(chǎn)匯流模型選用等流時線法,一維管網(wǎng)匯流模型選用動力波法,時間步長取1 s,地表水量平衡計算步長采用1 min,報告步長取5 min。根據(jù)模型計算結(jié)果,將淹水深度輸出為.asc格式文件,利用內(nèi)澇淹水結(jié)果統(tǒng)計工具統(tǒng)計得到淹水時間與風(fēng)險等級并輸出為.asc格式文件。最后利用ArcGIS分別制作最大淹水深度、最大淹水時間、風(fēng)險等級主題圖,如圖12—圖14所示。

    對于風(fēng)險等級為中和高的區(qū)域,應(yīng)盡快采取措施加以治理,主要措施有管網(wǎng)改造、增建泵站、源頭控制(LID)、增加地表排水通道等。

    5結(jié)語

    USRAMS軟件在SWMM計算引擎的基礎(chǔ)上添加了等流時線模塊,通過二維地表淹水模擬可計算得到淹水深度、淹水時間以及內(nèi)澇風(fēng)險等級的分布情況,在此過程計算中采用水量平衡原理,加快了計算速度,且具有良好的穩(wěn)定性??傊?,USRAMS軟件可為內(nèi)澇治理措施的制定與管網(wǎng)系統(tǒng)的規(guī)劃設(shè)計提供參考。

    目前,國內(nèi)主流的模型軟件主要有MIKE URBAN,Inforworks,PCSWMM等,都集成或開發(fā)了獨立的LID模塊,而USRAMS的針對性較強,主要用于城市雨水管網(wǎng)排水能力評估與城市內(nèi)澇風(fēng)險分析,系統(tǒng)中并未集成LID模塊,但USRAMS參數(shù)需求少,地表產(chǎn)匯流模型只需徑流系數(shù)與匯水時間2個參數(shù),并且二維地表淹水模擬計算速度均優(yōu)于上述軟件。為了更好地適應(yīng)中國海綿城市建設(shè)與內(nèi)澇防治工作,USRAMS將逐步集成LID、大排水系統(tǒng)設(shè)計等模塊。

    參考文獻/References:

    [1]侯玉棟,李樹平,周巍巍,等. 2011年國內(nèi)媒體報道城市暴雨事件分析[J]. 給水排水, 2012,38(sup):44-49.

    HOU Yudong,LI Shuping,ZHOU Weiwei,et al. Analysis of Chinese media reports on urban storm events in 2011[J]. Water & Wastewater Engineering, 2012,38(sup):44-49.

    [2]王軍, 馬洪濤. 城市排澇規(guī)劃有關(guān)問題探討[J]. 給水排水, 2014,40(3):9-12.

    WANG Jun, MA Hongtao. Discussion on some issues about urban logging prevention planning[J]. Water & Wastewater Engineering, 2014,40(3):9-12.

    [3]周玉文. 城市排水(雨水)防澇工程的系統(tǒng)架構(gòu)[J]. 給水排水, 2015,41(12):1-5.

    ZHOU Yuwen.System architecture of waterlogging prevention on water drainage (rainwater)system[J]. Water & Wastewater Engineering, 2015,41(12):1-5.

    [4]謝映霞. 排水防澇 重任在肩——寫給《城市排水(雨水)防澇綜合規(guī)劃編制大綱》頒布一周年[J]. 給水排水, 2014,40(6):1-3.

    XIE Yingxia. Drainage and waterlogging prevention, heavy burden on shoulder——for the anniversary of “City drainage (rainwater)and waterlogging prevention program outline” [J]. Water & Wastewater Engineering, 2014,40(6):1-3.

    [5]車伍,葛裕坤,唐磊,等.我國城市排水(雨水)防澇綜合規(guī)劃剖析[J]. 中國給水排水, 2016,32(10):15-21.

    CHE Wu,GE Yukun,TANG Lei,et al. Analysis of integrated plan for urban drainage ( stormwater) and flood control[J].China Water & Wastewater, 2016,32(10):15-21.

    [6]GB 50014—2006,室外排水設(shè)計規(guī)范[S].

    [7]劉子龍. 基于模型的城市內(nèi)澇災(zāi)害防治理論與技術(shù)研究[D]. 北京: 北京工業(yè)大學(xué), 2015.

    LIU Zilong. The Model-Based Techniques and Theories for Urban Local Flooding Prevention[D].Beijing: Beijing University of Technology,2015.

    [8]姚宇. 基于GeoDatabase的城市排水管網(wǎng)建模的應(yīng)用研究[D]. 上海:同濟大學(xué), 2007.

    YAO Yu. Study on The Application of Urban Drainage Networks Modeling Based on GeoDatabase[D].Shanghai: Tongji University,2007.

    [9]MAKSIMOVIC' , PRODANOVIC' D, BOONYA-AROONNET S, et al. Overland flow and pathway analysis for modelling of urban pluvial flooding[J]. Journal of Hydraulic Research, 2009,47(4):512-523.

    [10]高琳. 鄭州市排水防澇系統(tǒng)設(shè)計降雨研究[D]. 北京: 北京工業(yè)大學(xué), 2016.

    GAO Lin. Research on the Design Rainfall of Urban Drainage and Flooding Prevention System in Zhengzhou[D]. Beijing: Beijing University of Technology,2016.

    [11]REILLY J, PIECHOTA T C. Actual storm events outperform synthetic design storms: A review of SCS curve number reliability[C]// Proceedings of World & Environmental Resources Congress 2005.Washington D C: American Society of Civil Engineers, 2005:456-463.

    [12]KEIFER C J,CHU H H. Synthetic storm pattern for drainage design[J]. Journal of the Hydraulics Division, 1957,83(4):1-25.

    [13]HUFF F A. Time distributions of heavy rainstorms in Illinois[J]. Water Resources Research, 1967,3(4):1007-1019.

    [14]YEN B C, CHOW V T. Design hyetographs for small drainage structures[J]. Journal of the Hydraulics Division, 1980,106(6):1055-1076.

    [15]周玉文,高琳,楊偉明,等. 基于水利雨型峰值的城市長歷時設(shè)計暴雨雨型制作方法[P].中國專利: CN201510781439.7,2016-03-23.

    [16]朱勇年. 設(shè)計暴雨雨型的選用——以杭州市為例[J]. 中國給水排水, 2016,32(1):94-96.

    ZHU Yongnian. Selection of design rainfall pattern: Case of Hangzhou[J].China Water & Wastewater, 2016,32(1):94-96.

    [17]周玉文, 孟昭魯, 宋軍. 城市雨水管網(wǎng)非線性運動波法模擬技術(shù)[J]. 給水排水, 1995,21(4):9-11.

    ZHOU Yuwen, MENG Zhaolu, SONG Jun.Simulation of urban storm sewer system by non-linear kinematic wave method[J]. Water & Wastewater Engineering, 1995,21(4):9-11.

    [18]周玉文, 孟昭魯. 城市雨水口流域等流時線法降雨徑流模擬模型[J]. 沈陽建筑工程學(xué)院學(xué)報, 1994,10(4):339-344.

    ZHOU Yuwen,MENG Zhaolu. Precipitation runoff simulation modeling by isochronal method for urban inlet basin watershed[J].Journal of Shenyang Architectural and Civil Engineering Institute,1994,10(4):339-344.

    [19]GB 51222—2017,城鎮(zhèn)內(nèi)澇防治技術(shù)規(guī)范[S].

    [20]周勁松. 城市總體規(guī)劃中的城市特色塑造——以滄州市城市總體規(guī)劃為例[J]. 中華建設(shè), 2011,20(4):74-75.

    ZHOU Jinsong. Urbanization in urban master planning:A case study of Cangzhou citys master plan[J].China Construction, 2011,20(4):74-75.

    猜你喜歡
    徑流系數(shù)
    降雨特征對半透水道路徑流系數(shù)的影響
    水資源保護(2022年3期)2022-05-31 03:06:36
    天鵝湖子湖治理對地表水環(huán)境影響的研究
    綠色科技(2020年2期)2020-06-24 12:50:52
    小流域截洪溝洪峰流量計算方法研究與比較
    無資料山丘小流域徑流系數(shù)的計算研究
    低影響開發(fā)技術(shù)在市政道路排水中的應(yīng)用分析
    低影響開發(fā)(LID)模式技術(shù)在城市排水中應(yīng)用研究
    基于物理模型的城市綜合徑流系數(shù)特性研究
    大清河流域流域性暴雨洪水演變規(guī)律淺析
    地下水(2016年3期)2016-06-27 06:32:16
    關(guān)于綠色雨水排水系統(tǒng)的研究進展
    坡度與降雨強度對狗牙根草地產(chǎn)流的影響
    好男人视频免费观看在线| 爱豆传媒免费全集在线观看| 秋霞在线观看毛片| 天天躁夜夜躁狠狠久久av| а√天堂www在线а√下载| 毛片女人毛片| 亚洲中文字幕日韩| 乱人视频在线观看| 国产精品无大码| 亚洲在久久综合| 麻豆国产av国片精品| 免费不卡的大黄色大毛片视频在线观看 | h日本视频在线播放| 久久精品国产鲁丝片午夜精品| 六月丁香七月| 亚州av有码| 欧美一区二区精品小视频在线| 国产精品不卡视频一区二区| 特大巨黑吊av在线直播| 国产成人影院久久av| 综合色丁香网| 国产成人精品婷婷| 亚洲欧美日韩无卡精品| 两性午夜刺激爽爽歪歪视频在线观看| 别揉我奶头 嗯啊视频| 国内精品宾馆在线| 听说在线观看完整版免费高清| 国产 一区精品| 日韩精品青青久久久久久| 欧美日本视频| 日日干狠狠操夜夜爽| 村上凉子中文字幕在线| 天堂av国产一区二区熟女人妻| 久久久久九九精品影院| 黄片wwwwww| 插逼视频在线观看| 成人综合一区亚洲| 精品无人区乱码1区二区| 久久久久久伊人网av| 国产高清三级在线| 一边亲一边摸免费视频| 岛国毛片在线播放| 桃色一区二区三区在线观看| 亚洲真实伦在线观看| 九九爱精品视频在线观看| 一边亲一边摸免费视频| 久久久久免费精品人妻一区二区| 久久韩国三级中文字幕| 悠悠久久av| 麻豆精品久久久久久蜜桃| 成人国产麻豆网| 亚洲天堂国产精品一区在线| 午夜a级毛片| 男人狂女人下面高潮的视频| 久久精品夜夜夜夜夜久久蜜豆| 免费无遮挡裸体视频| 亚洲在线观看片| 高清在线视频一区二区三区 | 亚洲av中文av极速乱| 久久精品人妻少妇| 精品人妻偷拍中文字幕| 在线天堂最新版资源| 色综合色国产| 国产伦理片在线播放av一区 | 波野结衣二区三区在线| av福利片在线观看| 白带黄色成豆腐渣| 成人亚洲欧美一区二区av| 熟女电影av网| 熟女电影av网| 全区人妻精品视频| 不卡视频在线观看欧美| 青春草视频在线免费观看| 亚洲欧美日韩高清在线视频| 亚洲不卡免费看| 成年av动漫网址| 日韩 亚洲 欧美在线| 国产美女午夜福利| 久久久国产成人免费| 亚洲欧美中文字幕日韩二区| 亚洲国产精品成人久久小说 | 久久亚洲国产成人精品v| 国内精品久久久久精免费| 欧美成人免费av一区二区三区| 麻豆av噜噜一区二区三区| 精品无人区乱码1区二区| 亚洲精品久久国产高清桃花| 亚洲av一区综合| 狠狠狠狠99中文字幕| 免费看美女性在线毛片视频| 亚洲精品色激情综合| 婷婷六月久久综合丁香| 亚洲成人精品中文字幕电影| 国产精品福利在线免费观看| 国产女主播在线喷水免费视频网站 | 欧美激情国产日韩精品一区| www.av在线官网国产| 搞女人的毛片| 日韩欧美在线乱码| 中文资源天堂在线| 日本熟妇午夜| .国产精品久久| 成人亚洲欧美一区二区av| 大香蕉久久网| 国产午夜精品论理片| 久久亚洲国产成人精品v| 国产精品永久免费网站| 国产中年淑女户外野战色| 精品无人区乱码1区二区| 日韩亚洲欧美综合| 在现免费观看毛片| 男人的好看免费观看在线视频| a级毛片a级免费在线| 亚洲无线在线观看| 欧美精品国产亚洲| 国产三级中文精品| 成年女人看的毛片在线观看| 国产一级毛片七仙女欲春2| 美女脱内裤让男人舔精品视频 | 熟女电影av网| 久久精品国产亚洲av涩爱 | 久久国产乱子免费精品| 日韩欧美一区二区三区在线观看| 中文字幕av成人在线电影| 亚洲激情五月婷婷啪啪| 又粗又爽又猛毛片免费看| 蜜臀久久99精品久久宅男| 此物有八面人人有两片| 九九久久精品国产亚洲av麻豆| 99久久成人亚洲精品观看| 亚洲欧美精品综合久久99| a级毛色黄片| 女人被狂操c到高潮| 国产麻豆成人av免费视频| 亚洲欧美清纯卡通| 日日摸夜夜添夜夜爱| 亚洲精品国产av成人精品| 久久国内精品自在自线图片| 久久久久久九九精品二区国产| 日本一本二区三区精品| 嘟嘟电影网在线观看| 99热这里只有是精品在线观看| 久久人人爽人人片av| 女的被弄到高潮叫床怎么办| 国产91av在线免费观看| 国产精品一二三区在线看| 婷婷色av中文字幕| 亚洲不卡免费看| 国产av麻豆久久久久久久| 久久久久久久久中文| 国产色婷婷99| 日韩在线高清观看一区二区三区| 一个人看视频在线观看www免费| 成人午夜精彩视频在线观看| 国产成人一区二区在线| 99久久九九国产精品国产免费| 可以在线观看的亚洲视频| 亚洲va在线va天堂va国产| 能在线免费看毛片的网站| a级毛片a级免费在线| 少妇被粗大猛烈的视频| 中文亚洲av片在线观看爽| 欧美最黄视频在线播放免费| a级毛色黄片| 欧美zozozo另类| 99久久中文字幕三级久久日本| 国产老妇女一区| 国产伦一二天堂av在线观看| 晚上一个人看的免费电影| 在线观看午夜福利视频| 美女国产视频在线观看| 神马国产精品三级电影在线观看| 天堂av国产一区二区熟女人妻| 免费观看的影片在线观看| 九九热线精品视视频播放| 五月伊人婷婷丁香| 三级经典国产精品| 伊人久久精品亚洲午夜| 日韩国内少妇激情av| 国产欧美日韩精品一区二区| 大香蕉久久网| 一进一出抽搐gif免费好疼| 波多野结衣高清无吗| 一进一出抽搐gif免费好疼| 亚洲激情五月婷婷啪啪| 亚洲欧美成人综合另类久久久 | 国产精品久久久久久久电影| 男的添女的下面高潮视频| 深夜a级毛片| 一边亲一边摸免费视频| 免费搜索国产男女视频| 成人漫画全彩无遮挡| 成人特级av手机在线观看| 国产一区二区三区在线臀色熟女| 国产一区二区三区在线臀色熟女| 日韩av在线大香蕉| 麻豆国产av国片精品| 蜜桃久久精品国产亚洲av| 国产av麻豆久久久久久久| 婷婷亚洲欧美| 永久网站在线| av免费观看日本| 久久久久久大精品| 日韩av在线大香蕉| 午夜免费激情av| 一进一出抽搐动态| 国产亚洲av片在线观看秒播厂 | 久久久欧美国产精品| 欧美日韩国产亚洲二区| 国产亚洲精品久久久久久毛片| 在线a可以看的网站| 日产精品乱码卡一卡2卡三| 国内揄拍国产精品人妻在线| 哪里可以看免费的av片| 国产成人a∨麻豆精品| 少妇人妻精品综合一区二区 | 久久久久久久久久久丰满| 国产激情偷乱视频一区二区| 久久精品国产亚洲av涩爱 | 男人舔女人下体高潮全视频| 波多野结衣高清无吗| 久久99热这里只有精品18| 99九九线精品视频在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 麻豆成人午夜福利视频| 欧洲精品卡2卡3卡4卡5卡区| 午夜久久久久精精品| 在现免费观看毛片| 亚洲一区高清亚洲精品| 99久国产av精品| 晚上一个人看的免费电影| 日韩欧美国产在线观看| 最近最新中文字幕大全电影3| 亚洲美女搞黄在线观看| 特级一级黄色大片| 日韩人妻高清精品专区| 国产不卡一卡二| 五月玫瑰六月丁香| 欧美xxxx黑人xx丫x性爽| 久久久久性生活片| 精品久久久噜噜| av福利片在线观看| 日韩精品有码人妻一区| 欧美性猛交╳xxx乱大交人| 国产精品嫩草影院av在线观看| 美女内射精品一级片tv| 麻豆精品久久久久久蜜桃| 色尼玛亚洲综合影院| 可以在线观看的亚洲视频| 禁无遮挡网站| 麻豆成人午夜福利视频| 99热这里只有精品一区| 好男人在线观看高清免费视频| 一区福利在线观看| 嫩草影院入口| 午夜福利视频1000在线观看| av女优亚洲男人天堂| 美女脱内裤让男人舔精品视频 | 亚洲国产色片| 国产一级毛片在线| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩高清在线视频| 啦啦啦观看免费观看视频高清| 不卡视频在线观看欧美| 亚洲丝袜综合中文字幕| 久久精品国产鲁丝片午夜精品| 51国产日韩欧美| 老司机福利观看| 国产亚洲av片在线观看秒播厂 | ponron亚洲| 亚洲中文字幕一区二区三区有码在线看| 九草在线视频观看| 边亲边吃奶的免费视频| 一区二区三区高清视频在线| 久久欧美精品欧美久久欧美| 自拍偷自拍亚洲精品老妇| 亚洲av成人av| 床上黄色一级片| 能在线免费看毛片的网站| 国产精品蜜桃在线观看 | 午夜老司机福利剧场| 1024手机看黄色片| 亚洲久久久久久中文字幕| 尾随美女入室| 成年女人看的毛片在线观看| 美女大奶头视频| 久久精品国产自在天天线| 亚洲无线在线观看| 久久午夜亚洲精品久久| 99热这里只有精品一区| 国产毛片a区久久久久| 久久人人精品亚洲av| 国产男人的电影天堂91| 日韩亚洲欧美综合| 免费观看精品视频网站| 免费观看在线日韩| 日日撸夜夜添| 97热精品久久久久久| 国产精品伦人一区二区| 国产精品99久久久久久久久| 亚洲成人久久性| 精品一区二区免费观看| 亚洲五月天丁香| 蜜桃亚洲精品一区二区三区| 少妇熟女aⅴ在线视频| 男女啪啪激烈高潮av片| 天堂√8在线中文| 国产黄色小视频在线观看| 美女 人体艺术 gogo| 国产在线精品亚洲第一网站| 成人亚洲欧美一区二区av| 国产精品一及| 老司机影院成人| 看黄色毛片网站| 身体一侧抽搐| 高清午夜精品一区二区三区 | 国产精品久久久久久久久免| 黄色配什么色好看| 菩萨蛮人人尽说江南好唐韦庄 | 男人舔女人下体高潮全视频| 欧美在线一区亚洲| 国产久久久一区二区三区| 亚洲国产精品合色在线| 亚洲乱码一区二区免费版| 99九九线精品视频在线观看视频| 99久久精品国产国产毛片| 免费在线观看成人毛片| 精品人妻一区二区三区麻豆| 欧美激情久久久久久爽电影| 国产日韩欧美在线精品| 高清午夜精品一区二区三区 | 国产亚洲欧美98| 99精品在免费线老司机午夜| 国产伦在线观看视频一区| 小说图片视频综合网站| 此物有八面人人有两片| 国产成年人精品一区二区| 亚洲国产欧美在线一区| 亚洲av一区综合| 欧美成人免费av一区二区三区| 十八禁国产超污无遮挡网站| 在线播放无遮挡| 亚洲成人av在线免费| 亚洲精品色激情综合| 一夜夜www| 精品免费久久久久久久清纯| 久久久精品欧美日韩精品| av在线播放精品| 日韩 亚洲 欧美在线| 国产大屁股一区二区在线视频| 欧美激情国产日韩精品一区| 99久久九九国产精品国产免费| 国产爱豆传媒在线观看| 欧美性感艳星| 久久久色成人| 欧美成人免费av一区二区三区| 精品一区二区三区视频在线| 好男人视频免费观看在线| 亚洲av一区综合| 国产精华一区二区三区| 国产一区二区激情短视频| 在线观看一区二区三区| 成人特级黄色片久久久久久久| 九草在线视频观看| 亚洲成av人片在线播放无| 特级一级黄色大片| 成年女人看的毛片在线观看| av在线蜜桃| 久久欧美精品欧美久久欧美| 99视频精品全部免费 在线| 欧美高清成人免费视频www| 国产成人精品一,二区 | 嫩草影院新地址| 99在线人妻在线中文字幕| 久久精品国产清高在天天线| 国产一区亚洲一区在线观看| 在线观看美女被高潮喷水网站| 亚洲精品国产成人久久av| 久久精品国产亚洲网站| 国产精品电影一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品成人久久久久久| 久久精品影院6| 日日摸夜夜添夜夜爱| 欧美色视频一区免费| 亚洲国产精品国产精品| 成人永久免费在线观看视频| 岛国毛片在线播放| 久久精品国产自在天天线| 国产亚洲av嫩草精品影院| 99久久人妻综合| 久久久久久久久久黄片| 国产av不卡久久| a级毛片a级免费在线| 色视频www国产| 国产精品国产高清国产av| 精品熟女少妇av免费看| 亚洲无线在线观看| 国产单亲对白刺激| 成人欧美大片| 欧美潮喷喷水| 国内少妇人妻偷人精品xxx网站| 日韩欧美一区二区三区在线观看| 色综合站精品国产| 亚洲国产精品成人综合色| 99久久人妻综合| 欧美+日韩+精品| 国产三级在线视频| 一夜夜www| 能在线免费观看的黄片| 久久久久久久久久久丰满| 少妇的逼水好多| 国产高清激情床上av| 偷拍熟女少妇极品色| 男人和女人高潮做爰伦理| 在线播放无遮挡| 久久久午夜欧美精品| 偷拍熟女少妇极品色| 高清在线视频一区二区三区 | 国产精品乱码一区二三区的特点| 亚洲精品乱码久久久久久按摩| 男人舔女人下体高潮全视频| 有码 亚洲区| 日本黄色片子视频| av专区在线播放| 我的老师免费观看完整版| 日日撸夜夜添| 成人漫画全彩无遮挡| 国产精品久久电影中文字幕| 天堂网av新在线| 亚洲性久久影院| 麻豆成人av视频| 黄片wwwwww| 黄色视频,在线免费观看| 草草在线视频免费看| 人妻久久中文字幕网| 午夜精品国产一区二区电影 | 亚洲中文字幕一区二区三区有码在线看| АⅤ资源中文在线天堂| 国产精品嫩草影院av在线观看| 国产又黄又爽又无遮挡在线| 欧美日韩在线观看h| 亚洲最大成人av| 日本一二三区视频观看| 天天一区二区日本电影三级| 12—13女人毛片做爰片一| 美女内射精品一级片tv| 久久久a久久爽久久v久久| 五月伊人婷婷丁香| а√天堂www在线а√下载| 熟女人妻精品中文字幕| 欧美最新免费一区二区三区| 欧美日韩乱码在线| 国产成人freesex在线| 国产精品嫩草影院av在线观看| 国产av麻豆久久久久久久| 我要看日韩黄色一级片| 成人午夜精彩视频在线观看| 一级毛片久久久久久久久女| 亚洲av中文字字幕乱码综合| 精品免费久久久久久久清纯| 亚洲国产精品国产精品| 哪里可以看免费的av片| 中文在线观看免费www的网站| 国产蜜桃级精品一区二区三区| 久久久久网色| 天堂av国产一区二区熟女人妻| 久久精品国产自在天天线| 午夜激情福利司机影院| ponron亚洲| 久久久久性生活片| 此物有八面人人有两片| 午夜免费激情av| 亚洲在久久综合| videossex国产| 一个人看视频在线观看www免费| 亚洲av第一区精品v没综合| 欧美成人精品欧美一级黄| 高清日韩中文字幕在线| 国产69精品久久久久777片| 哪个播放器可以免费观看大片| 禁无遮挡网站| 国内精品美女久久久久久| 亚洲欧美日韩卡通动漫| 国产一级毛片七仙女欲春2| 乱人视频在线观看| 尾随美女入室| 最近视频中文字幕2019在线8| 欧美日韩国产亚洲二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久精品国产鲁丝片午夜精品| 男女下面进入的视频免费午夜| 国产日本99.免费观看| 精品国产三级普通话版| 日本三级黄在线观看| 草草在线视频免费看| 亚洲av电影不卡..在线观看| 国产成人91sexporn| 蜜桃亚洲精品一区二区三区| 久久精品国产清高在天天线| 成人鲁丝片一二三区免费| 热99在线观看视频| 又爽又黄无遮挡网站| 高清日韩中文字幕在线| 日本黄色片子视频| 亚洲精品日韩av片在线观看| 亚洲国产精品久久男人天堂| 国产av在哪里看| 国产黄色视频一区二区在线观看 | .国产精品久久| 国产成人影院久久av| 国语自产精品视频在线第100页| 欧美成人一区二区免费高清观看| 国产伦一二天堂av在线观看| 99久久无色码亚洲精品果冻| h日本视频在线播放| 99久久精品国产国产毛片| 国产亚洲5aaaaa淫片| 久久精品夜色国产| 男女下面进入的视频免费午夜| 精品无人区乱码1区二区| 欧美日本亚洲视频在线播放| 99热6这里只有精品| 色噜噜av男人的天堂激情| 久久99热6这里只有精品| 免费观看a级毛片全部| 精品久久久久久久末码| 变态另类成人亚洲欧美熟女| 精品人妻熟女av久视频| 国产高潮美女av| 国产一级毛片在线| 精品久久久久久久久亚洲| 日韩欧美国产在线观看| 国产真实伦视频高清在线观看| 少妇熟女欧美另类| 全区人妻精品视频| a级毛色黄片| 一区福利在线观看| 高清毛片免费看| 久久久久久国产a免费观看| av在线老鸭窝| 久久久久久久午夜电影| 日韩视频在线欧美| 一夜夜www| 国产 一区 欧美 日韩| 欧美性感艳星| 看黄色毛片网站| 99在线人妻在线中文字幕| 欧美人与善性xxx| 精品免费久久久久久久清纯| 国产老妇女一区| 亚洲精品久久国产高清桃花| 天堂中文最新版在线下载 | 久久99热6这里只有精品| www.色视频.com| 欧美日本视频| 精品午夜福利在线看| 中国美白少妇内射xxxbb| 国模一区二区三区四区视频| 国产亚洲精品久久久久久毛片| 少妇人妻精品综合一区二区 | 日韩强制内射视频| 亚洲av中文字字幕乱码综合| 一区二区三区免费毛片| 看免费成人av毛片| 国产老妇伦熟女老妇高清| 看非洲黑人一级黄片| 亚洲国产精品久久男人天堂| 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区| 国产亚洲91精品色在线| 在线播放无遮挡| 观看美女的网站| 精品熟女少妇av免费看| 亚洲最大成人手机在线| 亚洲国产精品sss在线观看| 在线免费十八禁| 色噜噜av男人的天堂激情| 桃色一区二区三区在线观看| 日日啪夜夜撸| 午夜久久久久精精品| 免费观看a级毛片全部| 美女被艹到高潮喷水动态| 人妻久久中文字幕网| 日韩精品有码人妻一区| 免费av观看视频| 成人特级av手机在线观看| 熟女电影av网| 婷婷色av中文字幕| 日韩中字成人| 99久久精品一区二区三区| 中文精品一卡2卡3卡4更新| 国产精品人妻久久久久久| 99久久精品国产国产毛片| 国产精品久久久久久久电影| 干丝袜人妻中文字幕| 亚洲熟妇中文字幕五十中出| 成人美女网站在线观看视频| а√天堂www在线а√下载| 日韩欧美在线乱码| 毛片女人毛片| 欧美日韩在线观看h| 蜜桃久久精品国产亚洲av| 听说在线观看完整版免费高清| 亚洲av电影不卡..在线观看| 亚洲国产欧美在线一区| 国产又黄又爽又无遮挡在线| а√天堂www在线а√下载| 亚洲欧美精品自产自拍| 狂野欧美白嫩少妇大欣赏| 在线免费观看的www视频| 尾随美女入室| 色播亚洲综合网| 久久热精品热| 午夜爱爱视频在线播放|