• <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)流的影響
    国产午夜精品久久久久久| 久久精品亚洲精品国产色婷小说| 日韩欧美国产在线观看| 国产人伦9x9x在线观看| 日本在线视频免费播放| 自拍欧美九色日韩亚洲蝌蚪91| 少妇的丰满在线观看| 女人被狂操c到高潮| 国产日韩一区二区三区精品不卡| 美女国产高潮福利片在线看| 国产精品98久久久久久宅男小说| 欧美日韩瑟瑟在线播放| 精品高清国产在线一区| 757午夜福利合集在线观看| 日韩国内少妇激情av| 欧美国产精品va在线观看不卡| 亚洲最大成人中文| 免费无遮挡裸体视频| 欧美激情久久久久久爽电影 | 精品久久久久久,| 成人三级黄色视频| 久久国产精品男人的天堂亚洲| 一进一出好大好爽视频| 亚洲欧美一区二区三区黑人| 首页视频小说图片口味搜索| 国产日韩一区二区三区精品不卡| 国产97色在线日韩免费| 两性夫妻黄色片| 色av中文字幕| 人妻丰满熟妇av一区二区三区| 久久国产亚洲av麻豆专区| cao死你这个sao货| 可以在线观看的亚洲视频| 人成视频在线观看免费观看| 国产极品粉嫩免费观看在线| 成人三级做爰电影| 成人亚洲精品av一区二区| 日韩高清综合在线| 欧美老熟妇乱子伦牲交| 啦啦啦免费观看视频1| 精品国产亚洲在线| av免费在线观看网站| 国产1区2区3区精品| 日本撒尿小便嘘嘘汇集6| 国产在线观看jvid| 99久久久亚洲精品蜜臀av| 国产精品一区二区免费欧美| 性少妇av在线| 高清毛片免费观看视频网站| 久久国产精品男人的天堂亚洲| 欧美精品啪啪一区二区三区| 真人一进一出gif抽搐免费| 美女午夜性视频免费| 日本免费一区二区三区高清不卡 | 欧美黑人精品巨大| 成人永久免费在线观看视频| 久久国产精品人妻蜜桃| 日本在线视频免费播放| 国产色视频综合| 日日干狠狠操夜夜爽| av中文乱码字幕在线| 黄色成人免费大全| 午夜日韩欧美国产| 我的亚洲天堂| 免费久久久久久久精品成人欧美视频| 操出白浆在线播放| 可以在线观看的亚洲视频| 亚洲一卡2卡3卡4卡5卡精品中文| 18禁美女被吸乳视频| 国产亚洲欧美在线一区二区| 久久人人97超碰香蕉20202| 国产精品一区二区精品视频观看| av天堂在线播放| 美女免费视频网站| 在线十欧美十亚洲十日本专区| 亚洲欧洲精品一区二区精品久久久| 丝袜美足系列| 女人被躁到高潮嗷嗷叫费观| 亚洲第一青青草原| videosex国产| 日日干狠狠操夜夜爽| 久9热在线精品视频| 精品欧美一区二区三区在线| 国产成人免费无遮挡视频| 欧美日韩精品网址| 精品久久久久久,| 91精品三级在线观看| 一夜夜www| 午夜免费激情av| 成人国语在线视频| 夜夜躁狠狠躁天天躁| 日韩国内少妇激情av| av在线播放免费不卡| www.999成人在线观看| tocl精华| 最新美女视频免费是黄的| 亚洲国产高清在线一区二区三 | 多毛熟女@视频| 欧美色视频一区免费| 亚洲五月婷婷丁香| 变态另类丝袜制服| 日韩 欧美 亚洲 中文字幕| 嫩草影视91久久| 久久香蕉精品热| 亚洲自偷自拍图片 自拍| 99久久99久久久精品蜜桃| 日韩大尺度精品在线看网址 | 非洲黑人性xxxx精品又粗又长| 日日干狠狠操夜夜爽| 久久人妻熟女aⅴ| 亚洲久久久国产精品| 在线观看66精品国产| 看片在线看免费视频| 淫妇啪啪啪对白视频| 精品一区二区三区av网在线观看| 麻豆一二三区av精品| 中文字幕高清在线视频| 国产成人av激情在线播放| 别揉我奶头~嗯~啊~动态视频| 中国美女看黄片| 极品教师在线免费播放| 长腿黑丝高跟| 丁香欧美五月| 国产在线观看jvid| avwww免费| 国产av在哪里看| 夜夜看夜夜爽夜夜摸| 欧美日本视频| 日日夜夜操网爽| 久久久久久久久中文| 9色porny在线观看| 91麻豆av在线| 操出白浆在线播放| 日韩一卡2卡3卡4卡2021年| 在线av久久热| 精品午夜福利视频在线观看一区| 操美女的视频在线观看| 免费在线观看亚洲国产| 三级毛片av免费| 亚洲色图 男人天堂 中文字幕| 亚洲人成电影观看| 两性夫妻黄色片| 在线播放国产精品三级| 国产片内射在线| 日日摸夜夜添夜夜添小说| www.精华液| 国产精品久久电影中文字幕| 国产片内射在线| 色尼玛亚洲综合影院| 久久性视频一级片| 久久久精品国产亚洲av高清涩受| 国产精品亚洲av一区麻豆| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲欧美98| 狂野欧美激情性xxxx| 精品欧美国产一区二区三| x7x7x7水蜜桃| 人人澡人人妻人| 中文字幕av电影在线播放| 国内精品久久久久精免费| 日韩欧美国产在线观看| 婷婷丁香在线五月| 国产精品久久电影中文字幕| 正在播放国产对白刺激| 91精品国产国语对白视频| netflix在线观看网站| 一区二区三区国产精品乱码| 日日爽夜夜爽网站| 久久热在线av| 免费久久久久久久精品成人欧美视频| 可以免费在线观看a视频的电影网站| 好男人在线观看高清免费视频 | 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩亚洲国产一区二区在线观看| 色尼玛亚洲综合影院| 精品欧美国产一区二区三| 婷婷丁香在线五月| 久久久久精品国产欧美久久久| 亚洲va日本ⅴa欧美va伊人久久| 日韩精品青青久久久久久| 人人妻人人爽人人添夜夜欢视频| 国产色视频综合| 久久久国产成人免费| 成人三级做爰电影| 9热在线视频观看99| 丁香六月欧美| 美女免费视频网站| 免费在线观看黄色视频的| 久久性视频一级片| 99在线视频只有这里精品首页| 午夜免费激情av| 黄色片一级片一级黄色片| 女人精品久久久久毛片| 精品欧美国产一区二区三| 夜夜爽天天搞| 51午夜福利影视在线观看| 久久精品国产亚洲av香蕉五月| 女生性感内裤真人,穿戴方法视频| 成人亚洲精品一区在线观看| 亚洲aⅴ乱码一区二区在线播放 | 国产欧美日韩一区二区三| 亚洲av成人av| 天堂动漫精品| 久久草成人影院| 成人免费观看视频高清| 国产av精品麻豆| 极品教师在线免费播放| 满18在线观看网站| 在线观看午夜福利视频| av超薄肉色丝袜交足视频| 亚洲aⅴ乱码一区二区在线播放 | 黑人巨大精品欧美一区二区蜜桃| 亚洲国产精品999在线| 精品国产一区二区三区四区第35| 欧美在线黄色| 亚洲成人国产一区在线观看| 丰满的人妻完整版| АⅤ资源中文在线天堂| 日韩欧美免费精品| 99国产精品一区二区蜜桃av| 中亚洲国语对白在线视频| bbb黄色大片| 亚洲av电影不卡..在线观看| 亚洲欧美精品综合久久99| 老司机午夜福利在线观看视频| 女人被狂操c到高潮| 大码成人一级视频| 熟女少妇亚洲综合色aaa.| 正在播放国产对白刺激| 成人永久免费在线观看视频| 国产精品乱码一区二三区的特点 | 午夜免费鲁丝| 国产精品一区二区免费欧美| 99国产精品99久久久久| 自拍欧美九色日韩亚洲蝌蚪91| 免费搜索国产男女视频| 国产精品影院久久| 天天躁狠狠躁夜夜躁狠狠躁| 一a级毛片在线观看| 午夜福利,免费看| 国产精品自产拍在线观看55亚洲| 亚洲av熟女| 91国产中文字幕| 成人手机av| 国产高清有码在线观看视频 | 国产一区二区在线av高清观看| 午夜福利免费观看在线| 宅男免费午夜| 一区二区日韩欧美中文字幕| 长腿黑丝高跟| 免费人成视频x8x8入口观看| 国产精品美女特级片免费视频播放器 | 啦啦啦免费观看视频1| 欧美黄色片欧美黄色片| 在线免费观看的www视频| 亚洲精品美女久久av网站| 日本a在线网址| 久久久久国内视频| 99久久精品国产亚洲精品| 久久国产精品人妻蜜桃| 亚洲最大成人中文| 亚洲成人免费电影在线观看| 久久久国产欧美日韩av| 无遮挡黄片免费观看| 国产免费av片在线观看野外av| 亚洲精品美女久久av网站| 欧美日韩亚洲综合一区二区三区_| 免费人成视频x8x8入口观看| 后天国语完整版免费观看| 国产亚洲欧美98| 大码成人一级视频| 操美女的视频在线观看| 欧美精品亚洲一区二区| 一本综合久久免费| 天天躁夜夜躁狠狠躁躁| 国产91精品成人一区二区三区| 91av网站免费观看| 最近最新免费中文字幕在线| 精品不卡国产一区二区三区| 久久精品人人爽人人爽视色| 久久久精品欧美日韩精品| 久久久久国内视频| 男女床上黄色一级片免费看| 免费看美女性在线毛片视频| 欧美乱妇无乱码| 久久草成人影院| 精品久久久久久久久久免费视频| 国内精品久久久久精免费| 法律面前人人平等表现在哪些方面| 亚洲 欧美一区二区三区| 亚洲色图av天堂| 亚洲欧美激情在线| 精品人妻在线不人妻| 国产一区二区三区在线臀色熟女| 午夜亚洲福利在线播放| 精品熟女少妇八av免费久了| 丁香六月欧美| 亚洲最大成人中文| 欧美精品亚洲一区二区| 亚洲av片天天在线观看| 精品久久久精品久久久| 黄色丝袜av网址大全| 日韩国内少妇激情av| 成人国产综合亚洲| 国产一卡二卡三卡精品| 99国产精品一区二区蜜桃av| 欧洲精品卡2卡3卡4卡5卡区| 麻豆一二三区av精品| 国产亚洲精品一区二区www| 亚洲国产欧美日韩在线播放| 国产av一区二区精品久久| 日本免费a在线| 亚洲人成电影免费在线| 免费在线观看黄色视频的| 人妻久久中文字幕网| 无限看片的www在线观看| av欧美777| 极品人妻少妇av视频| 99国产精品一区二区三区| 人人妻人人澡人人看| 欧美日韩亚洲国产一区二区在线观看| 色哟哟哟哟哟哟| 久久久精品欧美日韩精品| 色播亚洲综合网| 久久精品亚洲精品国产色婷小说| 久热这里只有精品99| 亚洲av熟女| 久久香蕉激情| 一夜夜www| 69av精品久久久久久| 波多野结衣一区麻豆| 免费在线观看亚洲国产| 我的亚洲天堂| 在线观看免费视频日本深夜| 久久香蕉精品热| 97碰自拍视频| 国产不卡一卡二| 一级作爱视频免费观看| 日韩欧美一区视频在线观看| 国产av精品麻豆| 久久伊人香网站| 久久久久久免费高清国产稀缺| 国产一区二区三区综合在线观看| 狂野欧美激情性xxxx| 精品国产国语对白av| 激情在线观看视频在线高清| 久久久久国内视频| 国产av精品麻豆| 性少妇av在线| 久久人人爽av亚洲精品天堂| 不卡av一区二区三区| 久久婷婷成人综合色麻豆| 99re在线观看精品视频| 久久香蕉精品热| 91精品国产国语对白视频| 老司机深夜福利视频在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲美女黄片视频| 国产成人影院久久av| 50天的宝宝边吃奶边哭怎么回事| 不卡av一区二区三区| 精品久久久久久久久久免费视频| 一级毛片高清免费大全| 少妇 在线观看| 国产精品久久久av美女十八| 两个人视频免费观看高清| 亚洲精品国产色婷婷电影| 丝袜美足系列| 欧美+亚洲+日韩+国产| 国产成年人精品一区二区| 最近最新中文字幕大全免费视频| 久久久国产成人免费| 精品久久久久久久久久免费视频| 男男h啪啪无遮挡| 亚洲一区二区三区色噜噜| 亚洲av美国av| 成人手机av| 久久久久久久久免费视频了| 亚洲一区二区三区色噜噜| 国产av精品麻豆| 中文字幕最新亚洲高清| 国产精品一区二区精品视频观看| 怎么达到女性高潮| 亚洲成国产人片在线观看| 亚洲av成人一区二区三| 91大片在线观看| 欧美黄色片欧美黄色片| 老熟妇仑乱视频hdxx| 免费看十八禁软件| netflix在线观看网站| 人人妻,人人澡人人爽秒播| 久久亚洲真实| 成年人黄色毛片网站| 757午夜福利合集在线观看| 国产精品,欧美在线| 99久久综合精品五月天人人| 91九色精品人成在线观看| 久久国产精品影院| 亚洲aⅴ乱码一区二区在线播放 | 黄频高清免费视频| 很黄的视频免费| 麻豆久久精品国产亚洲av| 大陆偷拍与自拍| 欧美一级毛片孕妇| 波多野结衣高清无吗| 一本久久中文字幕| 亚洲一区二区三区色噜噜| 亚洲性夜色夜夜综合| 十分钟在线观看高清视频www| 国产成人免费无遮挡视频| 午夜视频精品福利| 久久香蕉国产精品| 欧美老熟妇乱子伦牲交| 国产精品久久视频播放| 午夜精品久久久久久毛片777| 国产成人av激情在线播放| 999久久久国产精品视频| 久久久久久久午夜电影| 亚洲精华国产精华精| 国产精品 欧美亚洲| 欧美乱妇无乱码| 国产精品一区二区免费欧美| 亚洲一区二区三区色噜噜| 亚洲中文av在线| 在线观看日韩欧美| 国产99久久九九免费精品| 亚洲 国产 在线| 在线天堂中文资源库| 精品国产一区二区久久| 精品国产国语对白av| 99久久99久久久精品蜜桃| 性欧美人与动物交配| 久久天躁狠狠躁夜夜2o2o| 黑人欧美特级aaaaaa片| 激情在线观看视频在线高清| 搡老岳熟女国产| 母亲3免费完整高清在线观看| 亚洲人成伊人成综合网2020| 欧美日本视频| 一级毛片女人18水好多| 啦啦啦韩国在线观看视频| 一区二区三区激情视频| 中文字幕人妻丝袜一区二区| 国产精品秋霞免费鲁丝片| 亚洲免费av在线视频| 首页视频小说图片口味搜索| 大香蕉久久成人网| 99香蕉大伊视频| 国产99久久九九免费精品| 午夜精品国产一区二区电影| 久久久水蜜桃国产精品网| 一边摸一边抽搐一进一小说| 熟女少妇亚洲综合色aaa.| 黄网站色视频无遮挡免费观看| 天堂√8在线中文| av中文乱码字幕在线| 91麻豆精品激情在线观看国产| 久久婷婷人人爽人人干人人爱 | 亚洲欧美日韩另类电影网站| 国产高清有码在线观看视频 | 怎么达到女性高潮| 99久久久亚洲精品蜜臀av| 女人精品久久久久毛片| 成年女人毛片免费观看观看9| 午夜福利高清视频| 两性夫妻黄色片| 变态另类成人亚洲欧美熟女 | 天天躁狠狠躁夜夜躁狠狠躁| 女人精品久久久久毛片| 国产免费av片在线观看野外av| 国产aⅴ精品一区二区三区波| 在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 久久中文看片网| 法律面前人人平等表现在哪些方面| 久久狼人影院| 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区视频在线观看免费| 最近最新中文字幕大全免费视频| 国产午夜福利久久久久久| 国产精品久久视频播放| 一进一出抽搐gif免费好疼| 青草久久国产| 国产野战对白在线观看| 国产麻豆成人av免费视频| 国产一级毛片七仙女欲春2 | 满18在线观看网站| 欧美成人午夜精品| 色精品久久人妻99蜜桃| 亚洲三区欧美一区| 视频在线观看一区二区三区| 免费在线观看影片大全网站| 波多野结衣一区麻豆| 亚洲自拍偷在线| 国产亚洲欧美98| 午夜久久久在线观看| 久久欧美精品欧美久久欧美| 久热这里只有精品99| 窝窝影院91人妻| 国产不卡一卡二| 精品午夜福利视频在线观看一区| 日本五十路高清| 久久天躁狠狠躁夜夜2o2o| 啦啦啦观看免费观看视频高清 | 啦啦啦免费观看视频1| 亚洲欧美一区二区三区黑人| 免费看十八禁软件| 亚洲黑人精品在线| 久久人妻av系列| 午夜影院日韩av| 大香蕉久久成人网| 国产精品亚洲av一区麻豆| 国产亚洲精品综合一区在线观看 | 丝袜美腿诱惑在线| 亚洲天堂国产精品一区在线| 国产日韩一区二区三区精品不卡| 午夜影院日韩av| 一级片免费观看大全| 婷婷丁香在线五月| 国产成人av激情在线播放| 亚洲欧美日韩另类电影网站| 丁香欧美五月| 亚洲熟妇熟女久久| 久久久久国内视频| 欧美绝顶高潮抽搐喷水| 久热这里只有精品99| 欧美人与性动交α欧美精品济南到| 搞女人的毛片| 中文字幕av电影在线播放| 亚洲国产毛片av蜜桃av| 高清毛片免费观看视频网站| 亚洲午夜理论影院| 亚洲av成人不卡在线观看播放网| 亚洲中文字幕一区二区三区有码在线看 | 手机成人av网站| 午夜影院日韩av| 国产精品久久视频播放| 中文字幕av电影在线播放| 国产精品99久久99久久久不卡| 亚洲熟妇中文字幕五十中出| 亚洲狠狠婷婷综合久久图片| 欧美中文日本在线观看视频| 欧美乱码精品一区二区三区| 精品一品国产午夜福利视频| 两个人视频免费观看高清| 午夜久久久久精精品| 女人被躁到高潮嗷嗷叫费观| 国产91精品成人一区二区三区| 很黄的视频免费| 国产精品98久久久久久宅男小说| 亚洲熟女毛片儿| 国产精品久久久人人做人人爽| 在线视频色国产色| 成人亚洲精品一区在线观看| 男女下面插进去视频免费观看| 国产成人av激情在线播放| 国产精品久久久av美女十八| 女性生殖器流出的白浆| 久久久久久亚洲精品国产蜜桃av| 女人爽到高潮嗷嗷叫在线视频| 久久国产精品男人的天堂亚洲| 国产av一区二区精品久久| 久久国产精品人妻蜜桃| 国产野战对白在线观看| 97人妻精品一区二区三区麻豆 | 亚洲一区高清亚洲精品| 亚洲专区字幕在线| 伊人久久大香线蕉亚洲五| 亚洲中文av在线| 在线观看午夜福利视频| 精品久久久精品久久久| 天堂√8在线中文| 夜夜夜夜夜久久久久| 亚洲成人免费电影在线观看| 免费观看人在逋| 可以在线观看毛片的网站| 精品人妻在线不人妻| 两人在一起打扑克的视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美大码av| 精品久久久久久,| 一级,二级,三级黄色视频| 高清在线国产一区| 手机成人av网站| 免费高清视频大片| 日本三级黄在线观看| 韩国精品一区二区三区| 一个人观看的视频www高清免费观看 | 两个人视频免费观看高清| 窝窝影院91人妻| 日韩欧美国产一区二区入口| 免费看十八禁软件| 在线观看免费日韩欧美大片| 亚洲色图综合在线观看| 搞女人的毛片| 午夜福利,免费看| 国产精品一区二区在线不卡| 18禁观看日本| 在线观看午夜福利视频| 在线十欧美十亚洲十日本专区| 999久久久精品免费观看国产| 亚洲精品美女久久av网站| 精品国产乱码久久久久久男人| 88av欧美| 精品欧美国产一区二区三| 曰老女人黄片| 亚洲av熟女| 国产亚洲av高清不卡| 午夜亚洲福利在线播放| 一区在线观看完整版| 变态另类成人亚洲欧美熟女 | 亚洲va日本ⅴa欧美va伊人久久| 十八禁网站免费在线| or卡值多少钱| 国产日韩一区二区三区精品不卡|