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

    陽澄淀泖區(qū)動態(tài)洪水風險圖的編制及其管理系統(tǒng)的開發(fā)

    2017-11-10 07:55:16鐘桂輝劉曙光張梟鳴胡子琛俞燕彥
    水利水電科技進展 2017年6期
    關(guān)鍵詞:圩區(qū)泵站洪水

    鐘桂輝,劉曙光,張梟鳴,胡子琛,俞燕彥

    (1.同濟大學土木工程學院,上海 200092; 2.北京艾森思科技有限公司,北京 100085)

    陽澄淀泖區(qū)動態(tài)洪水風險圖的編制及其管理系統(tǒng)的開發(fā)

    鐘桂輝1,劉曙光1,張梟鳴1,胡子琛1,俞燕彥2

    (1.同濟大學土木工程學院,上海 200092; 2.北京艾森思科技有限公司,北京 100085)

    基于大量的實測和調(diào)研資料,建立了能反映陽澄淀泖區(qū)現(xiàn)狀河網(wǎng)、圩區(qū)分布及水工建筑物調(diào)度的水動力模型,經(jīng)率定和驗證后能夠用于該區(qū)域的洪水模擬?;赪eb技術(shù)、GIS技術(shù)和數(shù)據(jù)庫技術(shù),研發(fā)了嵌套水動力模型的動態(tài)洪水風險圖管理系統(tǒng)。針對任意降雨、邊界條件、水工建筑物參數(shù)及調(diào)度參數(shù),系統(tǒng)能快速調(diào)用模型進行模擬計算并生成動態(tài)洪水風險圖。與在線預(yù)報數(shù)據(jù)結(jié)合,可實現(xiàn)洪水預(yù)報預(yù)警,為洪水風險圖在防洪減災(zāi)、洪水管理、洪水預(yù)報等方面提供技術(shù)支撐和應(yīng)用借鑒。

    洪水風險圖; 洪水管理系統(tǒng); 水動力模型; 陽澄淀泖區(qū)

    1968年美國會通過《全國洪水保險法》推動洪水風險圖誕生以來,洪水風險圖在世界各國得到了廣泛的應(yīng)用[1-6]。我國自2004年開始洪水風險圖編制試點工作,至今已經(jīng)在全國很多地區(qū)編制了包括防洪保護區(qū)、中小流域、蓄滯洪區(qū)及城區(qū)等不同類型的洪水風險圖,為我國洪水風險圖的繪制和洪水風險管理探索出了經(jīng)驗、理論和方法[7-13]。洪水風險圖作為防洪減災(zāi)的重要技術(shù)支撐,是制定流域防洪規(guī)劃、規(guī)范和約束經(jīng)濟社會的發(fā)展行為、部署防洪工程及非工程措施、開展洪水保險和防汛搶險救災(zāi)等工作的重要依據(jù)。

    隨著全球氣候變暖、海平面上升及城鎮(zhèn)化的快速發(fā)展,洪水風險已經(jīng)并將繼續(xù)發(fā)生顯著變化,洪水風險圖的靈活運用及動態(tài)更新已成為流域及各級洪水管理部門重點關(guān)注的問題。特別是太湖流域,人口眾多、經(jīng)濟發(fā)達,特殊的地理環(huán)境和氣候特點,使得該地區(qū)洪澇災(zāi)害頻繁,損失嚴重。全球氣候變暖使得太湖流域臺風和暴雨趨勢增強,海平面的上升又降低了流域洪水外排的能力。人口、資產(chǎn)密度的急速增長,土地利用方式的急劇改變,流域內(nèi)水面率與河網(wǎng)通達率的下降,洪水風險的不確定性及危害性在不斷增加[14-16]。另外,太湖流域水系復(fù)雜,河湖密布,近50%的面積位于圩區(qū)或城市大包圍之中,圩區(qū)及城市大包圍是平原低洼地區(qū)防洪除澇的有效措施,但大量的澇水外排,導致圩外河網(wǎng)水位迅速大幅上漲,如2015年6月,太湖流域遭受3場強降雨侵襲,圩區(qū)及城市大包圍的澇水外排使得大運河水位迅速上升,大運河全線水位超警戒,無錫段水位達到5.18 m(1991年歷史高水位4.88 m),區(qū)域及流域防洪壓力大大增加。防洪減災(zāi)及洪水風險管理需要洪水風險圖的支撐,也需要探索動態(tài)和實時洪水風險圖,預(yù)見未來洪水風險,推廣洪水風險圖在其他領(lǐng)域的應(yīng)用,并從而提高洪水風險的管理能力[17]。

    計算機技術(shù)及信息化平臺的快速發(fā)展,使集成水文、氣象、模型計算、數(shù)據(jù)同化、GIS展示等為一體的動態(tài)洪水圖編制成為可能。本文以太湖流域陽澄淀泖區(qū)為研究對象,基于大量的基礎(chǔ)實測和調(diào)研資料,建立了高精度的一、二維耦合的水動力模型,結(jié)合計算機技術(shù),將模型嵌套于洪水風險圖管理系統(tǒng)中,實現(xiàn)了動態(tài)洪水風險圖的編制。洪水風險的增長是氣候變化和社會經(jīng)濟易損性增長共同作用的結(jié)果[16]。因此,開展針對氣候變化、人類活動和防洪工程體系綜合影響下的洪水模擬,動態(tài)生成洪水風險圖進行洪水風險分析和預(yù)測,能有效增強區(qū)域的洪水綜合管理與應(yīng)急響應(yīng)的能力。

    1 水動力模型

    1.1 區(qū)域概況

    陽澄淀泖區(qū)位于太湖流域東北部,屬太湖下游的一個水利分區(qū),西至望虞河,北倚長江,南至太浦河,面積4 314 km2。區(qū)內(nèi)河網(wǎng)縱橫交錯,湖泊眾多,有大小河道2萬余條,湖蕩146個。

    區(qū)域外圍有長江控制線、環(huán)太湖控制線、望虞河東岸控制線、太浦河北岸控制線、運河澹臺湖段東岸控制線、淀山湖昆山堤段控制線六大防洪工程。這些控制線沿線設(shè)有各種閘門及泵站,并按照聯(lián)合調(diào)度的控制規(guī)則來控制和調(diào)節(jié)外圍洪水。區(qū)域內(nèi)部除了有通江骨干河道、大小湖蕩進行洪水調(diào)節(jié)外,還建有大小聯(lián)圩366個,面積占總研究區(qū)域面積的50%。每個圩區(qū)都有多個泵站和節(jié)制閘進行調(diào)節(jié),實現(xiàn)各自圩區(qū)的防洪排澇。

    陽澄淀泖區(qū)地勢低平,河道水面比降小,水流流速緩慢,上游受太湖洪水下泄影響,下游又受長江潮位頂托,排水難度大,洪澇滯蓄時間長。洪澇特點表現(xiàn)為:①洪澇災(zāi)害頻繁。根據(jù)歷史資料統(tǒng)計,平均每4~5 a發(fā)生一次洪澇災(zāi)害。新中國成立以來,相繼發(fā)生了1954、1991、1999年流域性大洪水。②平原河網(wǎng)地區(qū)洪澇不分。洪澇災(zāi)害主要以平原河網(wǎng)地區(qū)受災(zāi)為主,每遇大范圍持久降雨或局部大暴雨時,難以區(qū)分洪災(zāi)和澇災(zāi)。③洪澇災(zāi)害經(jīng)濟損失巨大。④城市化發(fā)展導致區(qū)域防洪形勢更趨嚴峻。城市基礎(chǔ)設(shè)施的大規(guī)模建設(shè)使不透水地面面積擴大、水面率減少,對產(chǎn)匯流特性產(chǎn)生了較大影響;城市防洪除澇能力的增強和區(qū)域大規(guī)模的圩區(qū)建設(shè),又造成地區(qū)排澇量進一步加大,致使圩外河道水位上漲加快,高水位持續(xù)時間延長,加重了流域區(qū)域防洪除澇壓力[17]。

    1.2 水動力模型的建立與驗證

    選擇丹麥水動力研究所開發(fā)的MIKE FLOOD軟件建立一、二維耦合模型。

    一維水動力模型的基本方程是圣維南方程組,包括水流連續(xù)方程和動量方程:

    (1)

    (2)

    式中:x為距離;t為時間;A為過水斷面面積;Q為流量;h為水位;q為旁側(cè)入流量;C為謝才系數(shù);R為水力半徑;g為重力加速度。

    二維水動力模型基于三維不可壓縮的N-S(Navier-Stockes)方程,服從Bousinesq渦黏假定和靜水壓力假定,采用非結(jié)構(gòu)網(wǎng)格與蛙跳格式(leap frog)、通過控制體積法求解連續(xù)方程和動量方程。

    耦合模型是將一維河網(wǎng)兩岸的堤防與二維地面進行側(cè)向連接,實現(xiàn)水量、動量的實時交換,保證數(shù)學模型在空間、時間上的計算精度[18-19]。

    1.2.1 水動力模型的建立

    a. 數(shù)據(jù)來源及處理。數(shù)據(jù)主要包括基礎(chǔ)地理數(shù)據(jù)、水文數(shù)據(jù) 、河道斷面數(shù)據(jù)。河流水系數(shù)據(jù)、DEM高程數(shù)據(jù)、道路交通及土地利用等基礎(chǔ)地理數(shù)據(jù)來自江蘇省測繪局提供的2013年全要素矢量地圖,比尺為1∶10 000。水文數(shù)據(jù)為歷史觀測數(shù)據(jù),包括1951—2015年23個水文站點的降水資料,典型年主要測站的水位資料,水文數(shù)據(jù)進行了數(shù)據(jù)的檢驗和校正。骨干河道的斷面為實測數(shù)據(jù),測量間距為2 km,斷面變化較大及轉(zhuǎn)折點處進行了加密。其他河道斷面為全國水利普查數(shù)據(jù)。模型中的斷面進行了插值處理,一般每500 m一個計算斷面。

    b. 河網(wǎng)概化。區(qū)域內(nèi)河網(wǎng)復(fù)雜,湖蕩連通,且50%的區(qū)域面積為圩區(qū)控制。為了保證模型的精度,建模時盡量真實再現(xiàn)原有河網(wǎng)狀況。將河道分為圩內(nèi)、圩外兩大類。圩外河道盡量按實際情況全面放入模型中,半高地和高地處的河道(主要集中在常熟、昆山和太倉地區(qū))按調(diào)蓄水面積進行概化。圩內(nèi)河道盡量保留與圩外河網(wǎng)直接由閘泵相連的骨干河道,考慮洪水期水閘基本關(guān)閉,模型只考慮泵站的布置和抽排能力。調(diào)查每個圩區(qū)的水面率,根據(jù)實際水面率對圩區(qū)增加了適量的調(diào)蓄水面。湖泊處理成虛擬河道加上調(diào)蓄水面,湖泊實際水面面積與虛擬河道水面面積的差作為附加調(diào)蓄面積添加在調(diào)蓄節(jié)點處。模型共概化河道949條(其中圩外河道483條,圩內(nèi)河道466條),湖泊50個,計算水位點6 988個、計算流量點4 699個,模型河網(wǎng)如圖1所示。

    圖1 陽澄淀泖區(qū)水系概化

    c. 水工建筑物處理。區(qū)域內(nèi)水工建筑物主要包括水閘和泵站兩種。望虞河、太湖、長江、太浦河沿線閘門按照實際閘門的位置、規(guī)模及調(diào)度情況加入模型,邊界河道和澹臺湖以上(非城區(qū)包圍)沿運河東岸上的城區(qū)包圍的水工建筑物也按實際閘門的位置、規(guī)模及調(diào)度情況模擬。因圩區(qū)內(nèi)有大量的水閘和泵站,在保證圩區(qū)排澇總量不變的基礎(chǔ)上,根據(jù)圩區(qū)規(guī)模、圩區(qū)排澇能力和圩區(qū)實際泵站分布進行概化處理。洪水期間圩區(qū)水閘關(guān)閉。經(jīng)概化處理,模型包含泵站480個、水閘86個。

    d. 城市大包圍的處理。城市大包圍相當于一個大圩區(qū),研究區(qū)域內(nèi)主要有蘇州城區(qū)、昆山城區(qū)、吳江城區(qū)、常熟城區(qū)和太倉城區(qū)。城市大包圩考慮了所有的水閘和泵站,水閘和泵站的位置、規(guī)模及調(diào)度按實際情況處理。

    e. 網(wǎng)格剖分及耦合連接。為保證模型的精度,將每個圩區(qū)單獨網(wǎng)格剖分,沿江的半高地與高地區(qū)域統(tǒng)一進行網(wǎng)格剖分,局部地區(qū)進行加密,共生成網(wǎng)格82 904個,平均網(wǎng)格面積0.04 km2。一、二維耦合采用側(cè)向連接,將一維河道的堤防與二維陸域進行連接。

    f. 邊界條件、初始條件及參數(shù)選取。陽澄淀泖區(qū)是太湖流域一個獨立的水利分區(qū),南、北和西邊分別以長江、望虞河、太浦河和太湖口門為開邊界。東部邊界主要包括青松大包圍與嘉寶大包圍,根據(jù)調(diào)度實際,陽澄淀泖區(qū)難以向東泄洪,從區(qū)域防洪不利情況考慮,不考慮通過青松、嘉寶大包圍向下游排水,僅考慮淀山湖通過攔路港向下游泄洪以及吳淞江與鹽鐵塘與下游的水量交換,因此東部邊界基本為閉邊界,僅泖甸站與黃渡站為水位控制的開邊界。京杭運河參照2015年實測最大值,給定流量188 m3/s。初始流量設(shè)為零,初始水位設(shè)為河道多年平均高水位。河道糙率的初始值選擇文獻[20]中各河道的給定值,陸域的糙率按《洪水風險圖編制技術(shù)細則(試行)》賦予初值,率定時局部進行調(diào)整。

    g. 設(shè)計暴雨。本文涉及的降雨有20年一遇、50年一遇、100年一遇最大24 h設(shè)計面暴雨;50年一遇、100年一遇最大15 d設(shè)計面暴雨;“91北”、“99南”最大30 d設(shè)計面暴雨。最大24 h設(shè)計暴雨選取了楓橋、平望、直塘、花橋30年以上長系列資料適線法計算,并進行點面折算,利用《江蘇省暴雨洪水圖集》[21]最大24 h設(shè)計暴雨概化雨型進行時段雨量分配,并同時控制1 h、6 h降雨量。蒸發(fā)和滲透的影響采用初損后損法扣損,不考慮初損,后損為1 mm/h。長歷時設(shè)計暴雨中,蒸發(fā)和滲透的影響較大,不能單純用扣損法。選擇區(qū)域內(nèi)23個雨量站1951—2012年長系列資料進行了長歷時設(shè)計暴雨的計算,并利用泰森多邊形進行點面關(guān)系折算,設(shè)計典型年為1999年型,最大15 d同頻控制。最后再利用文獻[20]的模型,按照土地類型進行蒸發(fā)和入滲計算,得到分片后的設(shè)計面凈雨。

    1.2.2 模型率定和驗證

    選取2009年7月21日至8月10日洪水進行率定,率定后各河道的曼寧系數(shù)介于0.0225~0.035之間。選取2012年8月6—15日、2013年10月6—13日兩場典型洪水進行驗證。湘城站、陳墓站和金家壩站是陽澄淀泖區(qū)最具代表性的水位站點,兩段驗證期內(nèi)各站點計算值與實測值相差較小(圖2),絕大部分水位誤差控制在5 cm以內(nèi), NASH確定性系數(shù)介于0.83~0.98之間(表1),表明模型精度良好。

    圖2 部分水位站實測水位與計算水位對比

    水位站名年份最高水位/m出現(xiàn)日期實測計算實測計算水位最大誤差/m確定性系數(shù)湘城陳墓巴城金家壩20123.673.6008-1008-100.070.8320133.843.7710-0910-090.070.9220123.753.6608-1008-100.090.9320133.823.8010-0910-100.020.9620123.593.5308-1008-100.060.8420133.843.7410-0910-090.100.9120123.763.6908-1008-100.070.9120134.023.9810-0910-100.040.98

    2 動態(tài)洪水風險圖

    所謂動態(tài)洪水風險圖,是指能任意設(shè)置不同頻率的降雨,更改模型邊界條件、任意設(shè)置水工建筑物參數(shù)和調(diào)度規(guī)則、預(yù)設(shè)不同潰口位置和潰口發(fā)展過程,實現(xiàn)快速調(diào)用模型計算并生成洪水風險圖。同時,洪水風險的過程能動態(tài)展示;不同的方案可以調(diào)用對比;重點地區(qū)的水位和流量可以進行關(guān)注、預(yù)測并報警。

    為完成上述目的,首先,基于Web技術(shù)、GIS技術(shù)和數(shù)據(jù)庫技術(shù)研發(fā)了B/S (Browser/Server)結(jié)構(gòu)動態(tài)洪水風險圖管理系統(tǒng)。考慮到該系統(tǒng)主要用于防洪規(guī)劃、風險分析、防洪決策,系統(tǒng)采用面向服務(wù)(SOA)的體系架構(gòu),將應(yīng)用程序的不同功能單元(數(shù)據(jù)庫、模型等)通過定義好的接口和契約進行鏈接,實現(xiàn)應(yīng)用程序之間的數(shù)據(jù)交換和功能實現(xiàn)。然后,開發(fā)了MIKEF LOOD模型相關(guān)模塊數(shù)據(jù)交互組件的接口,包括降雨更改接口、在線數(shù)據(jù)調(diào)用接口、水工建筑物更改及調(diào)用接口,預(yù)設(shè)的潰口接口等。系統(tǒng)可以靈活調(diào)用模型模塊,更改相關(guān)參數(shù),從而實現(xiàn)動態(tài)洪水風險圖在系統(tǒng)中的編制。

    2.1 動態(tài)洪水風險圖管理系統(tǒng)

    系統(tǒng)的開發(fā)基于JAVAEE技術(shù)[22],將數(shù)據(jù)庫訪問、企業(yè)級Java組件、命名和目錄服務(wù)、動態(tài)頁面生成、XML、JSON、事務(wù)管理等有機集成在一起,并且提供集群等高級特性,使之特別適合構(gòu)建復(fù)雜的大型應(yīng)用系統(tǒng),并且保證系統(tǒng)具有很好的可擴展性。采用輕量級的構(gòu)建方法,在POJO的基礎(chǔ)上,封裝了Spring、Spring MVC和Hibernate,實現(xiàn)了簡單的構(gòu)架、強大的功能、良好的擴展性以及在不同的應(yīng)用服務(wù)器之間的高可移植性。系統(tǒng)建立的主要目標是實現(xiàn)在線數(shù)據(jù)的及數(shù)據(jù)庫數(shù)據(jù)的存儲和傳遞、水動力模型的封裝和調(diào)用、洪水風險圖結(jié)果的動態(tài)展示、動態(tài)洪水風險圖的編制及部分關(guān)注點水位的預(yù)報預(yù)警。系統(tǒng)采用分層的設(shè)計架構(gòu),整個架構(gòu)分為4層:應(yīng)用層、認證層、組件層和數(shù)據(jù)層,各層都有一套定義好的接口和構(gòu)件,具體見圖3。

    系統(tǒng)的主界面如圖4所示,主要由菜單欄,圖層區(qū)和顯示區(qū)。菜單欄顯示了系統(tǒng)的主要功能,包括區(qū)域概況、靜態(tài)風險圖、動態(tài)計算、實時預(yù)報及結(jié)果查看5個部分。

    區(qū)域概況是區(qū)域的基本信息的概述,包括地理信息、社會經(jīng)濟、行政區(qū)劃及歷史洪水等。

    圖3 系統(tǒng)架構(gòu)

    圖4 系統(tǒng)主界面

    靜態(tài)洪水風險圖是對洪水風險圖成果進行分類展示,任意計算方案的洪水風險圖都可以保存入庫并顯示。點擊“靜態(tài)風險圖”菜單,系統(tǒng)會以分類表格的形式羅列出該區(qū)域的風險圖研究成果,包括名稱、類型、相關(guān)的說明以及風險圖圖件(圖4)。已完成的靜態(tài)洪水風險圖有3類:①歷史洪水風險圖,主要是2009年、2012年、2013年典型年洪水風險圖。②不同頻率設(shè)計暴雨下內(nèi)澇洪水風險圖,包括20年、50年、100年一遇最大24 h設(shè)計暴雨;50年、100年一遇最大15 d設(shè)計暴雨;“91北”、“99南”最大30 d設(shè)計暴雨下的洪水風險圖。③潰堤洪水風險圖,經(jīng)調(diào)研和論證后假定部分堤防發(fā)生潰堤,從而編制的潰堤洪水風險圖。

    動態(tài)計算即為改變降雨及邊界等參數(shù)的動態(tài)洪水風險方案計算,計算結(jié)果保存在“結(jié)果查看”中,并能生成動態(tài)洪水風險圖。

    實時預(yù)報是動態(tài)計算的延伸,將降雨和邊界條件與實時水雨情數(shù)據(jù)庫連接,系統(tǒng)一方面能對實時的降雨、水文站點水位、河道流量等數(shù)據(jù)進行查詢分析;另外一方面,可以用實時水雨情為變動參數(shù)調(diào)用模型滾動計算,展現(xiàn)當前水雨情下水動力模型模擬預(yù)測的洪水風險信息。系統(tǒng)還設(shè)置了研究區(qū)域重要關(guān)注點(斷面、河道水位、流量),并給關(guān)注點設(shè)置預(yù)警閥值,當模型滾動計算中關(guān)注點出現(xiàn)超出預(yù)警閥值時,能進行洪水預(yù)報預(yù)警。

    結(jié)果查看是以模型計算的結(jié)果為依據(jù),動態(tài)展示一維水位流量變化過程、二維洪水演進過程,自動生成洪水淹沒風險圖、洪水淹沒歷時風險圖及洪水到達時間洪水風險圖。

    2.2 動態(tài)洪水風險圖編制

    動態(tài)洪水風險圖通過調(diào)用系統(tǒng)中的水動力模型,進行降雨條件的任意更改,如改變成其他設(shè)計頻率的降雨,或者依據(jù)天氣預(yù)報情況更改區(qū)域內(nèi)水文站點的降雨,同時還可以改變其他邊界條件,如改變邊界的水位過程、潮位過程、水閘泵站的參數(shù)、潰口的位置等,實現(xiàn)對任意洪水方案的快速模擬,圖5為動態(tài)洪水風險圖更改參數(shù)設(shè)置的界面。

    動態(tài)風險圖參數(shù)的更改包括5個部分:方案信息、降雨條件、邊界條件、泵站、水閘及潰口。“方案信息”是對新計算方案基本信息的設(shè)置和保存,方案信息包括:方案名稱、計算開始時間、計算結(jié)束時間、輸出步長和方案描述信息?!敖涤陾l件”提供降雨時序數(shù)據(jù)的查看及批量修改,降雨時序數(shù)據(jù)柱狀圖展示,雨量站的地圖定位?!斑吔鐥l件”提供水位時序數(shù)據(jù)的查看和批量修改,水位過程線展示,邊界點的地圖定位?!氨谜尽碧峁┍谜緟?shù)的更改和設(shè)置,包括:泵站的起排水位、終排水位和排水能力,可分別按行政區(qū)劃分片批量設(shè)置,也可以單獨設(shè)置每個泵站的參數(shù)?!八l”提供水閘參數(shù)的更改和設(shè)置,與泵站類似,可以批量修改也可以單獨設(shè)置?!皾⒖凇笨梢约僭O(shè)潰口,進行潰口、漫堤的動態(tài)洪水分析計算。

    參數(shù)更改設(shè)置完成后,通過模型模擬計算接口驅(qū)動MIKE11水動力模塊、MIKE21水動力模塊和MIKE FLOOD耦合模型,迅速進行洪水風險分析計算,計算結(jié)果導出到數(shù)據(jù)庫接口,并按照一定的表結(jié)構(gòu)導入到數(shù)據(jù)庫中,形成一個新方案洪水風險圖數(shù)據(jù),本文模型較為精細,82 904個計算網(wǎng)格,模擬24 h暴雨內(nèi)澇的時間約為20 min。以新方案計算數(shù)據(jù)為基礎(chǔ),通過風險圖配圖、添加輔助信息、說明信息和洪水風險統(tǒng)計信息,在結(jié)果查看模塊中生成新的洪水風險圖,從而完成一個動態(tài)洪水風險圖的編制。洪水風險的時空模擬結(jié)果可以在“結(jié)果查看”中進行回放、暫停及動畫展示。動態(tài)洪水風險圖可以用JPG等不同格式的不同圖幅、比例尺和分辨率查看和輸出。

    3 結(jié) 語

    本文基于大量實測數(shù)據(jù),采用一、二維耦合技術(shù)建立了陽澄淀泖區(qū)平原河網(wǎng)水動力模型,該模型包括河道949條,泵站480個,水閘86個,計算網(wǎng)格82 904個。模型的建立充分考慮了平原復(fù)雜河網(wǎng)間水量的交換、水工建筑物的影響及調(diào)度、圩區(qū)及城市大包圍的影響等。采用歷史洪澇數(shù)據(jù)對模型進行了率定和驗證,絕大部分水位誤差小于5 cm,精度較高,能夠用于該地區(qū)其他洪水方案的洪水風險分析和洪水風險圖的編制。

    通過多學科交叉、多模型耦合及先進的計算機軟件,研發(fā)了基于Web技術(shù)、GIS技術(shù)和數(shù)據(jù)庫技術(shù)的動態(tài)洪水風險圖編制系統(tǒng)。該系統(tǒng)將水動力模型通過接口嵌入其中,實現(xiàn)已完成方案的洪水風險圖查詢及動態(tài)展示;對任意組合洪水方案進行快速模擬計算,生成動態(tài)洪澇風險圖;與實時水雨工情數(shù)據(jù)庫連接,滾動計算并預(yù)報預(yù)警。系統(tǒng)能為防洪減災(zāi)、洪水管理、城市規(guī)劃避洪搶險提供決策依據(jù)。當然,該系統(tǒng)還需要進一步完善,如需要進一步提高模型的運行速度,今后考慮將其他水文模型及水動力模型入庫,實現(xiàn)任意模型的選擇、模擬及洪水預(yù)報預(yù)警。

    [1] VOUSDOUKAS M,VOUKOUVALAS E,MENTASCHI L,et al.Developments in large-scale coastal flood hazard mapping[J].Natural Hazard and Earth System,2016,16(8): 1841-1853.

    [2] GUSYEV M,GADEKE A,CULLMANN J,et al.Connecting global and local-scale flood risk assessment: a case study of the Rhine River basin flood hazard[J].Flood Risk Management,2014,9(4):343-354.

    [3] MENTZAFOU A,MARKGIANNI V,DIMITRIOU E.The use of geospatial technologies in flood hazard mapping and assessment: case study from River Evros[J].Pure and Applied Geophysics,2017,174(2):679-700.

    [4] BUCHELE B,KREIBICH H,KRON A,et al.Flood-risk mapping: contributions towards an enhanced assessment of extreme events and associated risks[J].Natural Hazards & Earth System Sciences,2006,6(4):485-503.

    [5] MOTEVALLI A,VAFAKHAH M.Flood hazard mapping using synthesis hydraulic and geomorphic properties at watershed scale[J].Stochastic Environmental Research and Risk Assessment,2016,30(7): 1889-1900.

    [6] BALDASSARRE G D,CASTELLARIN A,MONTANARI A,et al.Probability-weighted hazard maps for comparing different flood risk management strategies: a case study[J].Natural Hazards,2009,50(3): 479-496.

    [7] ZHONG Guihui,LIU Shuguang,HAN Chao,et al.Urban flood mapping for Jiaxing City based on hydrodynamic modeling and GIS analysis[J].Coastal Research,2014,68: 168-175.

    [8] 李娜,程曉陶,邱紹偉.上海市城區(qū)洪水風險圖制作[J].中國防汛抗旱,2009,19(6):38-41.(LI Na,CHENG Xiaotao,QIU Shaowei.Flood risk map making in Shanghai City[J].China Flood & Drought Managemen,2009,19(6): 38-41.(in Chinese))

    [9] 王煒.二維洪水數(shù)值模擬在太浦河洪水風險圖編制中的應(yīng)用[D].上海:同濟大學,2008.

    [10] 李昌志,黃金池,何曉燕,等.丹江口水庫洪水風險圖編制[J].中國防汛抗旱,2008,18(4): 63-66.(LI Changzhi,HUANG Jinchi,HE Xiaoyan,et al.Flood risk map making in Danjiangkou Reservoir[J].China Flood & Drought Management,2008,18(4): 63-66.(in Chinese))

    [11] 李航,于德萬,劉心玲.松花江洪水風險圖編制及管理應(yīng)用[J].東北水利水電,2016(3): 43-44.(LI Hang,YU Wande,LIU Xinling.Flood risk mapping in Songhua River and its management application[J].Water Resources & Hydropower of Northeast China,2016(3): 43-44.(in Chinese))

    [12] 王靜,喻朝慶,程曉陶,等.太湖流域大尺度洪水分析中圩區(qū)影響洪澇分布的模型[J].水利水電技術(shù),2010,41(9): 91-96.(WANG Jing,YU Chaoqing,CHENG Xiaotao,et al.Simulating the impacts of polder areas on distribution of flooding in Taihu Basin for broad scale flood analysis[J].Water Resources and Hydropower Engineering,2010,41(9): 91-96.(in Chinese))

    [13] 葉愛民,劉曙光,韓超,等.MIKEFLOOD耦合模型在杭嘉湖流域嘉興地區(qū)洪水風險圖編制工作中的應(yīng)用[J].中國防汛抗旱,2016,26(2): 56-60.(YE Aimin,LIU Shuguang,HAN Chao,et al.Application of MIKE FLOOD coupling model in flood risk mapping at Jiaxing area of Zhejiang Province[J].China Flood & Drought Management,2016,26(2): 56-60.(in Chinese))

    [14] 徐宗學,劉瀏.太湖流域氣候變化檢測與未來氣候變化情景預(yù)估[J].水利水電科技進展,2012,32(1): 1-7.(XU Zongxue,LIU Liu.Detection of climate change and projection of future climate change scenarios in Taihu Lake Basin[J].Advances in Science and Technology of Water Resources,2012,32(1): 1-7.(in Chinese))

    [15] 關(guān)穎慧.長江流域極端氣候變化及其未來趨勢預(yù)測[D].楊凌: 西北農(nóng)林科技大學,2015.

    [16] HALL J W,EWANS E P,PENNING-ROWSELL E C,et al.Quantified scenarios analysis of drivers and impacts of changing flood risk in England and Wales: 2030—2100[J].Environmental Hazards,2003,5(3/4):51-65.

    [17] 葉劍春,章杭慧.太湖流域洪水風險管理實踐與思考[J].水利水電科技進展,2015,35(5):136-141.(YE Jianchun,ZHANG Hanghui.Practices and thinking of flood risk management in Taihu Lake Basin[J].Advances in Science and Technology of Water Resources,2015,35(5):136-141.(in Chinese))

    [18] Danish Hydraulic Institute.A modelling system for rivers and channels,MIKE 11 reference manual[R] Danmark: Danish Hydraulic Institute,2014.

    [19] Danish Hydraulic Institute.MIKEFLOOD 1D-2D modelling user manual[R] Danmark: Danish Hydraulic Institute,2014.

    [20] 程文輝,王船海,朱琰.太湖流域模型[M].南京: 河海大學出版社,2006.

    [21] 江蘇省水文總站.江蘇省暴雨洪水圖集[G].南京:江蘇省水文總站,1984.

    [22] 曹鳴鵬,趙偉,許林英.J2EE技術(shù)及其實現(xiàn)[J].計算機應(yīng)用,2001,10(21): 20-23.(CAO Mingpeng,ZHAO Wei,XU Linying.The J2EE technology and its implementation[J].Computer Applications,2001,10(21): 20-23.(in Chinese))

    DynamicfloodriskmappinganddevelopmentofamanagementsystemforYangcheng&DianmaoDistrict//

    ZHONG Guihui1, LIU Shuguang1, ZHANG Xiaoming1,HU Zicheng1,YU Yanyan2

    (1.CollegeofCivilEngineering,TongjiUniversity,Shanghai200092,China; 2.BeijingAisensiScienceandTechnologyCo.,Ltd.,Beijing100085,China)

    A hydrodynamic model integrated with river networks, ponler areas and the controlling of hydraulic structures was established. After validation and calibration with sufficient measured data and current situation investigation, this numerical model can successfully simulate the flood. Based on Web, GIS and database technology, the management system of dynamic flood risk maps coupledwith the hydrodynamic model was developed. Dynamic flood risk mapping can be produced after the employment of the hydrodynamic model,which can deal with multiple customized data, including precipitation, boundary conditions, parameters of hydraulic structuresand scheduling disciplines. Combined with online forecasting data, the management system can provide technical support and references for the application of flood risk maps on flood control and disaster relief, as well as flood management and forecasting in China.

    dynamic flood risk mapping; flood management system; hydrodynamic model; Yangcheng & Dianmao District

    “十二五”國家科技支撐計劃 (2014BAL05B02)

    鐘桂輝(1971—),女,副教授,博士研究生,主要從事防洪減災(zāi)研究。E-mail:guihui_zhong@#edu.cn

    劉曙光(1962—),男,教授,博士,主要從事防洪減災(zāi)研究。E-mail:liusgliu@#edu.cn

    10.3880/j.issn.1006-7647.2017.06.011

    TV877

    A

    1006-7647(2017)06-0062-07

    2017-01-09 編輯:鄭孝宇)

    猜你喜歡
    圩區(qū)泵站洪水
    基于GIS 空間分析的上海市青浦圩區(qū)現(xiàn)狀研究
    張家邊涌泵站建設(shè)難點及技術(shù)創(chuàng)新實踐
    洪水時遇到電線低垂或折斷該怎么辦
    又見洪水(外二首)
    天津詩人(2017年2期)2017-11-29 01:24:25
    2016年河南省己建成泵站數(shù)量
    圩區(qū)治理對區(qū)域防洪排澇影響之我見
    全省已建成泵站數(shù)量
    洪水來了
    幼兒畫刊(2016年8期)2016-02-28 21:00:52
    上海奉賢莊行地區(qū)圩區(qū)除澇能力研究
    河南省2014年已建成泵站數(shù)量
    免费黄网站久久成人精品| 亚洲国产欧美日韩在线播放| 人妻 亚洲 视频| 久久久久久久久久人人人人人人| 精品国产乱码久久久久久小说| 亚洲精品国产色婷婷电影| 一区二区三区乱码不卡18| 久久久久久人人人人人| 色5月婷婷丁香| 国产一区二区在线观看av| 亚洲精品国产色婷婷电影| 亚洲欧美中文字幕日韩二区| 巨乳人妻的诱惑在线观看| 夫妻性生交免费视频一级片| 七月丁香在线播放| 高清欧美精品videossex| 亚洲精品一区蜜桃| 最近中文字幕2019免费版| 精品国产一区二区久久| 午夜老司机福利剧场| 成年av动漫网址| 色5月婷婷丁香| 汤姆久久久久久久影院中文字幕| 国产精品女同一区二区软件| 99热全是精品| 少妇被粗大的猛进出69影院 | 国产熟女午夜一区二区三区| 中文字幕最新亚洲高清| 国产69精品久久久久777片| 免费看不卡的av| 22中文网久久字幕| 成人影院久久| 亚洲五月色婷婷综合| 亚洲丝袜综合中文字幕| 这个男人来自地球电影免费观看 | av免费在线看不卡| 亚洲av免费高清在线观看| 午夜福利视频在线观看免费| 永久免费av网站大全| 精品久久蜜臀av无| 国产片特级美女逼逼视频| 午夜影院在线不卡| a级毛片在线看网站| 一边亲一边摸免费视频| 中文乱码字字幕精品一区二区三区| 久久久国产欧美日韩av| 人成视频在线观看免费观看| 亚洲欧洲精品一区二区精品久久久 | 水蜜桃什么品种好| 色5月婷婷丁香| 日本与韩国留学比较| 久久影院123| 少妇熟女欧美另类| 亚洲性久久影院| 在线观看免费日韩欧美大片| 校园人妻丝袜中文字幕| 婷婷成人精品国产| 波野结衣二区三区在线| 啦啦啦在线观看免费高清www| 精品久久久精品久久久| 免费少妇av软件| 国产成人精品久久久久久| av又黄又爽大尺度在线免费看| 国内精品宾馆在线| 久久久欧美国产精品| 亚洲一区二区三区欧美精品| av国产精品久久久久影院| 欧美成人精品欧美一级黄| 在线天堂中文资源库| av福利片在线| 一本久久精品| 欧美变态另类bdsm刘玥| av播播在线观看一区| 青青草视频在线视频观看| 婷婷色av中文字幕| 夫妻性生交免费视频一级片| 99久久中文字幕三级久久日本| 亚洲欧美清纯卡通| 亚洲人成网站在线观看播放| 国产熟女午夜一区二区三区| 久久午夜福利片| 国产在线一区二区三区精| 99精国产麻豆久久婷婷| 国产精品人妻久久久影院| a级毛片黄视频| 日韩精品免费视频一区二区三区 | 最近手机中文字幕大全| 夜夜爽夜夜爽视频| 免费人妻精品一区二区三区视频| 美女内射精品一级片tv| 亚洲精品久久成人aⅴ小说| 日韩av不卡免费在线播放| 男女高潮啪啪啪动态图| 亚洲成人一二三区av| 黄色配什么色好看| 亚洲综合色惰| 少妇 在线观看| 免费在线观看完整版高清| 香蕉国产在线看| 亚洲 欧美一区二区三区| 精品99又大又爽又粗少妇毛片| 免费高清在线观看日韩| 国产精品欧美亚洲77777| 99热6这里只有精品| 宅男免费午夜| 国产一区二区在线观看av| 99久久精品国产国产毛片| 精品一区二区三卡| 国产黄频视频在线观看| 精品久久久精品久久久| 精品人妻在线不人妻| 丝袜人妻中文字幕| 美女内射精品一级片tv| 成人午夜精彩视频在线观看| 丝袜人妻中文字幕| 免费观看无遮挡的男女| 精品国产一区二区三区四区第35| 天天躁夜夜躁狠狠久久av| 亚洲av电影在线观看一区二区三区| 国产精品久久久久久久久免| 亚洲欧美一区二区三区黑人 | 秋霞在线观看毛片| 岛国毛片在线播放| 亚洲国产欧美在线一区| 天堂俺去俺来也www色官网| 精品人妻熟女毛片av久久网站| 免费播放大片免费观看视频在线观看| 亚洲欧美中文字幕日韩二区| 深夜精品福利| 一本色道久久久久久精品综合| 国产精品久久久久久精品电影小说| 宅男免费午夜| 色视频在线一区二区三区| 麻豆乱淫一区二区| 日韩熟女老妇一区二区性免费视频| 美女国产视频在线观看| 日韩一区二区三区影片| 免费大片黄手机在线观看| 色婷婷av一区二区三区视频| 一边摸一边做爽爽视频免费| 三级国产精品片| 熟女av电影| 日韩中字成人| 国语对白做爰xxxⅹ性视频网站| 男人添女人高潮全过程视频| av福利片在线| 曰老女人黄片| 国产av一区二区精品久久| 久久99一区二区三区| 国产成人91sexporn| 老司机影院成人| 另类亚洲欧美激情| 国产成人91sexporn| 欧美日韩av久久| 亚洲欧美成人精品一区二区| 日韩不卡一区二区三区视频在线| 2021少妇久久久久久久久久久| 另类亚洲欧美激情| 亚洲av男天堂| 99热这里只有是精品在线观看| 国产一区二区激情短视频 | av不卡在线播放| 午夜视频国产福利| www.色视频.com| 亚洲国产色片| 久久精品人人爽人人爽视色| 狂野欧美激情性xxxx在线观看| 狠狠精品人妻久久久久久综合| 熟女电影av网| 91在线精品国自产拍蜜月| 亚洲av在线观看美女高潮| 精品福利永久在线观看| 亚洲成人手机| 日韩成人伦理影院| 少妇的丰满在线观看| 一边摸一边做爽爽视频免费| 欧美xxⅹ黑人| 亚洲国产精品一区三区| 国产色婷婷99| 一区二区av电影网| 日韩一区二区三区影片| 99精国产麻豆久久婷婷| 日本爱情动作片www.在线观看| 日日啪夜夜爽| 亚洲精品自拍成人| 午夜91福利影院| 国产av一区二区精品久久| 26uuu在线亚洲综合色| 欧美国产精品va在线观看不卡| 18禁国产床啪视频网站| xxx大片免费视频| 欧美3d第一页| 香蕉丝袜av| 久久久久人妻精品一区果冻| 亚洲国产精品国产精品| 成人黄色视频免费在线看| 夜夜骑夜夜射夜夜干| 亚洲欧美成人精品一区二区| 国产成人精品婷婷| 精品亚洲乱码少妇综合久久| 日本av免费视频播放| 午夜福利影视在线免费观看| 高清视频免费观看一区二区| 国产精品免费大片| 国产亚洲av片在线观看秒播厂| 亚洲av中文av极速乱| 大片免费播放器 马上看| 欧美变态另类bdsm刘玥| 日韩中文字幕视频在线看片| 熟妇人妻不卡中文字幕| av在线观看视频网站免费| 一本久久精品| av有码第一页| 黑人欧美特级aaaaaa片| 精品视频人人做人人爽| 成人国产麻豆网| 中文乱码字字幕精品一区二区三区| 欧美日本中文国产一区发布| 丰满乱子伦码专区| 天天影视国产精品| 久久这里只有精品19| 侵犯人妻中文字幕一二三四区| 久久久国产精品麻豆| 免费av中文字幕在线| 久久人人97超碰香蕉20202| 人成视频在线观看免费观看| 婷婷色综合www| 亚洲精品美女久久久久99蜜臀 | 国产成人午夜福利电影在线观看| 日韩三级伦理在线观看| 老熟女久久久| 日本色播在线视频| 99热全是精品| 一级片免费观看大全| 精品卡一卡二卡四卡免费| 国产午夜精品一二区理论片| 久久久国产一区二区| 午夜精品国产一区二区电影| 免费人妻精品一区二区三区视频| 中文字幕精品免费在线观看视频 | 又黄又粗又硬又大视频| 色94色欧美一区二区| 少妇人妻 视频| 国产精品人妻久久久影院| 国产精品久久久久久精品古装| 男女免费视频国产| 久久国产亚洲av麻豆专区| 精品人妻在线不人妻| 色婷婷av一区二区三区视频| 国产欧美日韩一区二区三区在线| 久久99蜜桃精品久久| 99久国产av精品国产电影| www.av在线官网国产| 国产亚洲午夜精品一区二区久久| 一级爰片在线观看| 桃花免费在线播放| 国产爽快片一区二区三区| 午夜日本视频在线| 水蜜桃什么品种好| 午夜福利网站1000一区二区三区| 美国免费a级毛片| 91精品三级在线观看| 亚洲欧洲日产国产| 巨乳人妻的诱惑在线观看| 美女国产高潮福利片在线看| 亚洲精品久久成人aⅴ小说| av在线观看视频网站免费| 欧美 亚洲 国产 日韩一| 亚洲国产日韩一区二区| 国产综合精华液| 午夜免费观看性视频| 日日撸夜夜添| 男人舔女人的私密视频| 日本午夜av视频| 国产成人aa在线观看| 人体艺术视频欧美日本| 热re99久久国产66热| 欧美人与善性xxx| 多毛熟女@视频| 91成人精品电影| 国产一区二区三区av在线| 人妻少妇偷人精品九色| 国产成人精品无人区| 王馨瑶露胸无遮挡在线观看| 欧美激情极品国产一区二区三区 | 热99久久久久精品小说推荐| 精品久久国产蜜桃| 国产xxxxx性猛交| 精品一区二区三区四区五区乱码 | 男女啪啪激烈高潮av片| a 毛片基地| 中文字幕另类日韩欧美亚洲嫩草| 国产免费一区二区三区四区乱码| 欧美激情 高清一区二区三区| 亚洲精品乱久久久久久| 内地一区二区视频在线| 亚洲激情五月婷婷啪啪| 性色avwww在线观看| 国产乱人偷精品视频| 日韩av在线免费看完整版不卡| 99久久综合免费| 国产又爽黄色视频| 日韩不卡一区二区三区视频在线| 天天躁夜夜躁狠狠躁躁| 精品国产乱码久久久久久小说| 人妻系列 视频| 国产女主播在线喷水免费视频网站| 一本久久精品| 全区人妻精品视频| 国产在视频线精品| 亚洲成色77777| 亚洲精品色激情综合| av卡一久久| 成人国产av品久久久| 国产综合精华液| 蜜桃国产av成人99| 亚洲成色77777| 亚洲精品色激情综合| 青春草国产在线视频| 免费在线观看完整版高清| 日本av免费视频播放| 国产成人免费无遮挡视频| 少妇熟女欧美另类| 久久久久国产网址| 日韩伦理黄色片| 最近手机中文字幕大全| 国产色爽女视频免费观看| 一级片免费观看大全| 久久影院123| 国产深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 五月伊人婷婷丁香| 丰满少妇做爰视频| 国产av码专区亚洲av| 91精品伊人久久大香线蕉| 看免费成人av毛片| 久久国产亚洲av麻豆专区| 免费看光身美女| 两个人看的免费小视频| 色5月婷婷丁香| 亚洲av电影在线观看一区二区三区| 另类亚洲欧美激情| 在线观看国产h片| 久久久欧美国产精品| 欧美少妇被猛烈插入视频| 精品国产一区二区三区四区第35| 国产激情久久老熟女| 国产亚洲精品第一综合不卡 | 久久狼人影院| 国产精品欧美亚洲77777| 夜夜爽夜夜爽视频| 亚洲欧洲日产国产| 热99久久久久精品小说推荐| 精品一区在线观看国产| 国产毛片在线视频| 人人妻人人添人人爽欧美一区卜| 九草在线视频观看| 男男h啪啪无遮挡| 黄色配什么色好看| 亚洲婷婷狠狠爱综合网| 全区人妻精品视频| 国产老妇伦熟女老妇高清| 高清黄色对白视频在线免费看| 一级毛片我不卡| 亚洲av电影在线进入| 大片免费播放器 马上看| 大香蕉久久网| 好男人视频免费观看在线| 久久久国产精品麻豆| 精品少妇内射三级| 99精国产麻豆久久婷婷| 热99国产精品久久久久久7| 人体艺术视频欧美日本| 久久精品国产综合久久久 | 97人妻天天添夜夜摸| 母亲3免费完整高清在线观看 | 飞空精品影院首页| 亚洲欧美成人综合另类久久久| 成年av动漫网址| 国产精品秋霞免费鲁丝片| 国产在线一区二区三区精| 美女国产视频在线观看| 国产xxxxx性猛交| 一边亲一边摸免费视频| 国产不卡av网站在线观看| 青春草亚洲视频在线观看| 日韩在线高清观看一区二区三区| 久久影院123| 欧美+日韩+精品| 啦啦啦中文免费视频观看日本| 午夜福利,免费看| 国产精品免费大片| 亚洲国产毛片av蜜桃av| av片东京热男人的天堂| 夜夜爽夜夜爽视频| 亚洲欧洲精品一区二区精品久久久 | 日韩av免费高清视频| 精品人妻一区二区三区麻豆| 中文精品一卡2卡3卡4更新| 国产色婷婷99| 成人影院久久| 日韩,欧美,国产一区二区三区| 久久97久久精品| 久久久久国产网址| 性色avwww在线观看| 亚洲第一区二区三区不卡| 在现免费观看毛片| 交换朋友夫妻互换小说| 久久99热6这里只有精品| 街头女战士在线观看网站| 黄色 视频免费看| 26uuu在线亚洲综合色| 久久精品国产a三级三级三级| av免费在线看不卡| 日韩一区二区视频免费看| 国产精品欧美亚洲77777| 少妇被粗大的猛进出69影院 | 精品一品国产午夜福利视频| 中国国产av一级| 欧美+日韩+精品| 观看美女的网站| 99国产精品免费福利视频| 啦啦啦中文免费视频观看日本| 国产成人精品福利久久| 咕卡用的链子| 男女下面插进去视频免费观看 | 五月伊人婷婷丁香| 日日啪夜夜爽| 国产亚洲欧美精品永久| 久久精品aⅴ一区二区三区四区 | 国产成人av激情在线播放| 岛国毛片在线播放| 欧美人与性动交α欧美软件 | 国产成人欧美| 99热网站在线观看| 国产亚洲av片在线观看秒播厂| 午夜激情久久久久久久| 亚洲国产精品999| 国精品久久久久久国模美| 99精国产麻豆久久婷婷| 国产片特级美女逼逼视频| 男女边摸边吃奶| 国产精品久久久久久精品古装| 看非洲黑人一级黄片| 99久久人妻综合| 国产精品国产三级国产av玫瑰| 天天影视国产精品| www.熟女人妻精品国产 | 综合色丁香网| 大香蕉久久网| 夜夜骑夜夜射夜夜干| 22中文网久久字幕| 久热久热在线精品观看| 十分钟在线观看高清视频www| 日韩精品免费视频一区二区三区 | 丁香六月天网| 一边亲一边摸免费视频| 久久久欧美国产精品| 午夜免费鲁丝| 最近中文字幕2019免费版| 美女国产高潮福利片在线看| 如日韩欧美国产精品一区二区三区| 亚洲精品国产av蜜桃| 各种免费的搞黄视频| 一区二区日韩欧美中文字幕 | 伊人久久国产一区二区| 国产亚洲欧美精品永久| 热99久久久久精品小说推荐| 久久 成人 亚洲| 黑丝袜美女国产一区| 亚洲国产精品成人久久小说| 黄片无遮挡物在线观看| 国产亚洲一区二区精品| 国产精品国产三级专区第一集| 美女主播在线视频| 一本久久精品| 99视频精品全部免费 在线| av国产久精品久网站免费入址| 精品酒店卫生间| 看免费av毛片| 久久久久久久久久成人| 成人漫画全彩无遮挡| 亚洲三级黄色毛片| 女人精品久久久久毛片| 午夜福利视频在线观看免费| 十分钟在线观看高清视频www| 狠狠精品人妻久久久久久综合| 亚洲精品视频女| 2021少妇久久久久久久久久久| 国产白丝娇喘喷水9色精品| 午夜免费鲁丝| 深夜精品福利| 久久综合国产亚洲精品| 亚洲精华国产精华液的使用体验| 国产亚洲欧美精品永久| 国产成人精品一,二区| 爱豆传媒免费全集在线观看| 欧美精品高潮呻吟av久久| 美女大奶头黄色视频| 国产亚洲一区二区精品| 亚洲色图 男人天堂 中文字幕 | 久热这里只有精品99| 成人18禁高潮啪啪吃奶动态图| 高清视频免费观看一区二区| 精品一区在线观看国产| 亚洲,欧美,日韩| 午夜免费鲁丝| videosex国产| 高清欧美精品videossex| 男女下面插进去视频免费观看 | 中文欧美无线码| 日韩一区二区三区影片| 亚洲精品日韩在线中文字幕| 亚洲欧美一区二区三区国产| 欧美精品一区二区大全| 下体分泌物呈黄色| 亚洲精品视频女| 亚洲精品456在线播放app| 免费播放大片免费观看视频在线观看| 亚洲五月色婷婷综合| 两性夫妻黄色片 | 香蕉丝袜av| 亚洲精品日本国产第一区| 2022亚洲国产成人精品| 欧美日韩亚洲高清精品| 日韩制服骚丝袜av| 天天影视国产精品| 国产片内射在线| 久久久久国产网址| 午夜91福利影院| 中文字幕另类日韩欧美亚洲嫩草| 亚洲综合精品二区| 日本免费在线观看一区| 亚洲精品日本国产第一区| 国产精品国产三级专区第一集| 香蕉精品网在线| 亚洲伊人久久精品综合| 涩涩av久久男人的天堂| 啦啦啦在线观看免费高清www| 国产xxxxx性猛交| 91久久精品国产一区二区三区| 日韩,欧美,国产一区二区三区| 一级毛片我不卡| 国产永久视频网站| 国产精品一二三区在线看| 日韩不卡一区二区三区视频在线| 新久久久久国产一级毛片| 国产男人的电影天堂91| 午夜免费男女啪啪视频观看| 亚洲色图综合在线观看| 亚洲成色77777| 亚洲精品日本国产第一区| 亚洲精品,欧美精品| 亚洲精品国产av蜜桃| av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 欧美另类一区| 午夜福利,免费看| 18禁在线无遮挡免费观看视频| 亚洲精品aⅴ在线观看| 99热这里只有是精品在线观看| 搡女人真爽免费视频火全软件| 日韩视频在线欧美| 精品第一国产精品| av有码第一页| 亚洲国产最新在线播放| 人成视频在线观看免费观看| 婷婷成人精品国产| 精品国产一区二区久久| 久久狼人影院| 亚洲精品久久午夜乱码| 欧美日韩一区二区视频在线观看视频在线| 丝袜脚勾引网站| 国产熟女午夜一区二区三区| 国产在线一区二区三区精| 一边亲一边摸免费视频| 欧美老熟妇乱子伦牲交| 啦啦啦视频在线资源免费观看| 寂寞人妻少妇视频99o| 欧美3d第一页| 国产一区二区在线观看av| 黄色 视频免费看| 国产精品国产三级国产av玫瑰| 18+在线观看网站| 久久青草综合色| 寂寞人妻少妇视频99o| 侵犯人妻中文字幕一二三四区| 成人黄色视频免费在线看| 亚洲图色成人| 国国产精品蜜臀av免费| 免费大片黄手机在线观看| 一本大道久久a久久精品| 中文天堂在线官网| 国产精品久久久久久精品古装| 色吧在线观看| 亚洲精华国产精华液的使用体验| 亚洲 欧美一区二区三区| 美女脱内裤让男人舔精品视频| 国产片特级美女逼逼视频| 一区在线观看完整版| av在线老鸭窝| 综合色丁香网| 一级毛片黄色毛片免费观看视频| 麻豆精品久久久久久蜜桃| 美女国产视频在线观看| 18禁国产床啪视频网站| 热re99久久国产66热| 国产免费福利视频在线观看| 在线免费观看不下载黄p国产| 精品国产乱码久久久久久小说| 天天影视国产精品| 校园人妻丝袜中文字幕| 插逼视频在线观看| 欧美日韩av久久| 日本av免费视频播放|