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

    湖泊水力調(diào)控對(duì)河湖連通伴生洪水風(fēng)險(xiǎn)的緩釋影響

    2022-06-07 05:26:00鄭云鶴劉殷竹于若蘭高學(xué)平
    水利學(xué)報(bào) 2022年3期
    關(guān)鍵詞:沭河駱馬湖河湖

    張 晨,鄭云鶴,劉殷竹,于若蘭,高學(xué)平

    (天津大學(xué) 水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300350)

    1 研究背景

    近幾十年來(lái),我國(guó)水資源與水環(huán)境問(wèn)題日益突出,水資源分布時(shí)空不均與新時(shí)代新發(fā)展的理念不相匹配,嚴(yán)重制約了經(jīng)濟(jì)社會(huì)的發(fā)展。河湖水系連通工程作為國(guó)家治水、用水的重要戰(zhàn)略手段,對(duì)區(qū)域水資源優(yōu)化統(tǒng)籌配置、河湖生態(tài)環(huán)境改善以及抵御水旱災(zāi)害具有重要作用[1]。關(guān)于河湖連通的定義,前人從連通功能[2]、連通特征[3]、連通的時(shí)空性科學(xué)范疇[4]等不同角度詮釋了河湖水系連通的概念與內(nèi)涵,如李原園等[5]闡述了河湖水系連通的定義,即以江河湖庫(kù)等為基礎(chǔ),通過(guò)人為的疏導(dǎo)、調(diào)度,建立江河湖庫(kù)水體之間的水力聯(lián)系。

    隨著河湖水系連通工程的大量興建,對(duì)不同類型和尺度的河湖水系連通進(jìn)行風(fēng)險(xiǎn)識(shí)別、預(yù)測(cè)與評(píng)估,建立應(yīng)對(duì)極端水文事件的河湖水系連通伴生風(fēng)險(xiǎn)管理機(jī)制,制定變化環(huán)境下河湖水系連通伴生風(fēng)險(xiǎn)的適應(yīng)性對(duì)策,成為河湖水系連通伴生風(fēng)險(xiǎn)管控方面的重要研究?jī)?nèi)容。以往研究主要以兩種方法進(jìn)行風(fēng)險(xiǎn)分析:(1)單純的物理過(guò)程模擬。其優(yōu)點(diǎn)是考慮了水流運(yùn)動(dòng)規(guī)律;缺點(diǎn)是風(fēng)險(xiǎn)情景假設(shè)致使工況方案受限,如從引水流量[6]、引調(diào)水時(shí)間[7]、水量分配[8]等單方面研究河湖水系連通伴生風(fēng)險(xiǎn)分析和管控措施,或分別利用HEC-RAS、MIKE、EFDC 等數(shù)值模擬洪水演進(jìn)過(guò)程,進(jìn)而分析洪水風(fēng)險(xiǎn)[9-10]。(2)單獨(dú)的模糊隨機(jī)分析。該方法雖然可以形成更普遍的隨機(jī)方案群,但缺乏物理機(jī)制過(guò)程,且會(huì)受到現(xiàn)有數(shù)據(jù)量和空間尺度的限制,如劉艷麗等[11]基于模糊環(huán)境構(gòu)建的水庫(kù)防洪風(fēng)險(xiǎn)分析模型。該模型結(jié)合了拉丁超立方抽樣與蒙特卡羅算法的優(yōu)點(diǎn),適用于單因素風(fēng)險(xiǎn)與多因素組合風(fēng)險(xiǎn)的分析。

    本文作者等[12]利用Vine Copula函數(shù)與水動(dòng)力模型模擬技術(shù),創(chuàng)建隨機(jī)水情條件下河湖水系連通伴生風(fēng)險(xiǎn)分析模型,探討隨機(jī)水情條件下河湖水系連通伴生風(fēng)險(xiǎn),實(shí)現(xiàn)風(fēng)險(xiǎn)隨機(jī)過(guò)程與物理過(guò)程的結(jié)合;并以沂沭泗水系下游河湖連通工程為例,進(jìn)一步探討湖泊水力調(diào)控對(duì)河湖連通伴生洪水風(fēng)險(xiǎn)的定量影響。

    2 研究區(qū)域與方法

    2.1 研究區(qū)域概況 沂沭泗流域位于淮河流域東北部,流域面積7.96 萬(wàn)km2,占淮河流域面積的29%。沂沭泗水系下游由韓莊運(yùn)河、邳蒼分洪道、中運(yùn)河、沂河下游、老沭河、新沂河以及駱馬湖組成,如圖1 所示。該河湖連通工程中,邳蒼分洪道于1958年開(kāi)挖,連接沂河與中運(yùn)河,防治了沂河上游分洪后的漫流危害;新沂河于1949年開(kāi)挖,承接駱馬湖及老沭河來(lái)水后匯入黃海,在新中國(guó)成立后的沂沭泗流域削減水患中發(fā)揮了巨大作用。沂沭河上中游均位于山區(qū),源短流急,且南四湖流域洪水起落較快,洪水匯集至沂沭河水系下游后,河道比降變緩,行洪不暢,存在發(fā)生洪水災(zāi)害的風(fēng)險(xiǎn)。

    圖1 沂沭泗水系連通狀況及水文站分布

    駱馬湖南北長(zhǎng)20 km,東西寬16 km。湖底高程約19.0~20.0 m,正常蓄水位23.0 m 時(shí),湖面面積375 km2,容積9.0 億m3;設(shè)計(jì)洪水位25.0 m 時(shí),容積15億m3;校核洪水位26.0 m 時(shí),總?cè)莘e19.0億m3。駱馬湖承接沂河下游、中運(yùn)河以及邳蒼分洪道來(lái)水,是南水北調(diào)東線工程第二座調(diào)蓄湖泊,保護(hù)下游1000 萬(wàn)畝耕地和600 萬(wàn)人生命財(cái)產(chǎn)安全[13],在沂沭泗水系下游的防洪體系中具有舉足輕重的作用。

    2.2 風(fēng)險(xiǎn)分析模型控制方程 對(duì)于一維、二維水動(dòng)力數(shù)學(xué)模型,其控制方程為連續(xù)性方程與運(yùn)動(dòng)方程,表述了水體運(yùn)動(dòng)時(shí),重力、壓力、黏滯力和慣性力之間的平衡關(guān)系[14-15]。由于風(fēng)險(xiǎn)過(guò)程發(fā)生的隨機(jī)性,采用改進(jìn)后的連續(xù)性方程與運(yùn)動(dòng)方程:

    式中:u、v 為沿x、y 方向的流速,m/s;g 為重力加速度,m/s2;h 為水深,m;D 為無(wú)量綱經(jīng)驗(yàn)系數(shù),與河床、邊坡形狀相關(guān);R 為水力半徑,m;H 為位置水頭,m;C 為謝才系數(shù),m1/2/s;ω為隨機(jī)源匯項(xiàng),s-1,ω=f (x,y,t) 。等式(2)左側(cè)為時(shí)變加速度項(xiàng)與對(duì)流項(xiàng),等式右側(cè)從左至右依次表示重力產(chǎn)生的壓力梯度、黏滯力項(xiàng)、河床底部摩擦及隨機(jī)源匯項(xiàng)。

    將C-Vine Copula n 維聯(lián)合密度函數(shù)兩邊積分得n 維累積分布函數(shù),求解n 維累計(jì)分布函數(shù)的逆函數(shù),可推導(dǎo)出以流量表述的隨機(jī)源匯項(xiàng)ω具體表達(dá)式為:

    式中:Q 為流量,m3/s;x1,x2,…,xn為關(guān)鍵風(fēng)險(xiǎn)因子;累積分布函數(shù)的逆函數(shù)F-1(x1,x2,…,xn)代表隨機(jī)生成的關(guān)鍵風(fēng)險(xiǎn)因子x1,x2,…,xn序列。具體計(jì)算方法見(jiàn)文獻(xiàn)[12]。

    2.3 風(fēng)險(xiǎn)分析模型建立及驗(yàn)證 建立沂沭泗水系下游河道一維水動(dòng)力模型,即韓莊運(yùn)河、中運(yùn)河、新沂河每隔1 km 劃分一個(gè)斷面,河段長(zhǎng)度總計(jì)251 km,斷面總數(shù)共253 個(gè)。駱馬湖二維水動(dòng)力模型的網(wǎng)格尺寸為100 m×100 m,共劃分46 213 個(gè)網(wǎng)格。

    一維水動(dòng)力模型選用2008年運(yùn)河站洪水過(guò)程率定模型糙率,再以2008年嶂山站洪水過(guò)程校核模型;二維水動(dòng)力模型選用水位庫(kù)容曲線驗(yàn)證模型。其中,運(yùn)河站位于駱馬湖上游的中運(yùn)河段,嶂山站位于駱馬湖下游的新沂河段。經(jīng)模型率定,韓莊運(yùn)河與中運(yùn)河糙率為0.06,新沂河糙率為0.035。模型誤差分析如圖2 所示,結(jié)果表明,平均絕對(duì)誤差(MAE)和均方根誤差(RMSE)均低于0.45 m、納什系數(shù)(NSE)大于0.85。

    圖2 風(fēng)險(xiǎn)分析模型水位模擬值與實(shí)測(cè)值比較

    當(dāng)沂沭泗水系遭遇洪水災(zāi)害時(shí),洪水通過(guò)韓莊運(yùn)河韓莊閘、沂河下游劉道口閘、老沭河勝利堰閘進(jìn)入沂沭泗水系下游及駱馬湖,故選取韓莊運(yùn)河、沂河下游、老沭河邊界斷面洪峰流量作為隨機(jī)變量構(gòu)建Copula 函數(shù)生成邊界條件。

    首先,選用4 種常見(jiàn)水文變量分布(威布爾分布、logistic 分布、正態(tài)分布及對(duì)數(shù)正態(tài)分布)分別對(duì)各變量觀測(cè)序列進(jìn)行邊緣分布擬合,通過(guò)AIC(赤池信息準(zhǔn)則)擬合優(yōu)度檢驗(yàn)比選出最適用的邊緣分布。各分布AIC 值列于表1,其K-S 檢驗(yàn)值均大于0.05,通過(guò)檢驗(yàn)。

    表1 各隨機(jī)變量不同分布類型對(duì)應(yīng)AIC 值

    由表1 結(jié)果可知,韓莊運(yùn)河洪峰流量與沂河下游洪峰流量最優(yōu)擬合分布均為威布爾分布,老沭河洪峰流量最優(yōu)擬合分布為對(duì)數(shù)正態(tài)分布。將韓莊運(yùn)河、沂河下游及老沭河洪峰流量序列用其對(duì)應(yīng)擬合邊緣分布函數(shù)進(jìn)行概率積分轉(zhuǎn)換,使新序列在[0,1]區(qū)間上均勻分布。根據(jù)AIC 準(zhǔn)則,第一層級(jí)從 十種常見(jiàn)的Copula 函 數(shù)(Gaussian Copula、Student-t Copula、Clayton Copula、Gumbel-Hougaard Copula、Frank Copula、Joe Copula、BB 族Copula)選取兩變量間最優(yōu)的Pair Copula 類型,結(jié)果為韓莊運(yùn)河與沂河下游洪峰選用Frank Copula,韓莊運(yùn)河與老沭河洪峰選用Clayton Copula(270 degrees);第二層級(jí)樹(shù)選用Normal Copula。最后通過(guò)貝葉斯估計(jì)方法估參,并利用極大似然估計(jì)對(duì)采用的Copula函數(shù)參數(shù)估計(jì)進(jìn)行修正,得Vine Copula 的p 值為0.55(大于0.05),通過(guò)擬合優(yōu)度檢驗(yàn)。

    2.4 洪水風(fēng)險(xiǎn)源與隨機(jī)邊界條件 結(jié)合沂沭泗流域面積廣、汛期洪水峰高量大、源短流急、來(lái)勢(shì)兇猛等特點(diǎn)[16],空間上劃分3 種風(fēng)險(xiǎn)源:風(fēng)險(xiǎn)源①,南四湖流域未遭遇洪水,沂沭河上游遭遇洪水,南四湖控制下泄,邳蒼分洪道來(lái)水進(jìn)入中運(yùn)河,繼而進(jìn)入駱馬湖,沂河來(lái)水進(jìn)入駱馬湖;風(fēng)險(xiǎn)源②,沂沭河上游未遭遇洪水,南四湖流域遭遇洪水,沂沭河上游來(lái)水東調(diào)進(jìn)入新沂河,南四湖敞泄,南四湖來(lái)水通過(guò)中運(yùn)河進(jìn)入駱馬湖;風(fēng)險(xiǎn)源③,沂沭泗流域全境遭遇洪水,南四湖洪水敞泄至韓莊運(yùn)河,邳蒼分洪道來(lái)水進(jìn)入中運(yùn)河,沂河來(lái)水進(jìn)入駱馬湖。

    根據(jù)歷史水文資料,分別考慮3 種風(fēng)險(xiǎn)源,在韓莊運(yùn)河、沂河下游、老沭河洪峰流量區(qū)間內(nèi)隨機(jī)組合出水情條件,作為風(fēng)險(xiǎn)分析模型的邊界條件,并形成模擬計(jì)算方案,列于表2。

    表2 模擬方案

    2.5 風(fēng)險(xiǎn)等級(jí)計(jì)算方法 本文選取流量、水位和流速作為定量評(píng)估洪水風(fēng)險(xiǎn)的指標(biāo),分別以校核流量、保證水位、不沖流速作為判別發(fā)生風(fēng)險(xiǎn)的限值。根據(jù)風(fēng)險(xiǎn)分析模型模擬結(jié)果,計(jì)算風(fēng)險(xiǎn)概率和風(fēng)險(xiǎn)強(qiáng)度,公式如下:

    式中:P 為風(fēng)險(xiǎn)概率;I 為洪水風(fēng)險(xiǎn)評(píng)估指標(biāo)的模擬值,IL為相應(yīng)指標(biāo)的風(fēng)險(xiǎn)限值,GI>IL為風(fēng)險(xiǎn)評(píng)估指標(biāo)模擬值大于風(fēng)險(xiǎn)限值的模擬組數(shù);GT為總模擬組數(shù);R 為風(fēng)險(xiǎn)強(qiáng)度;I′為大于風(fēng)險(xiǎn)限值的風(fēng)險(xiǎn)評(píng)估指標(biāo)模擬值。借鑒趙勇等[17]依據(jù)風(fēng)險(xiǎn)概率和風(fēng)險(xiǎn)影響劃分風(fēng)險(xiǎn)等級(jí)的方法,本文將式(4)得到的風(fēng)險(xiǎn)概率和風(fēng)險(xiǎn)強(qiáng)度形成風(fēng)險(xiǎn)等級(jí)矩陣,列于表3,進(jìn)而確定風(fēng)險(xiǎn)等級(jí)。

    表3 風(fēng)險(xiǎn)等級(jí)矩陣

    3 沂沭泗水系河湖連通伴生洪水風(fēng)險(xiǎn)分析

    當(dāng)沂沭河上游遭遇洪水,南四湖流域未遭遇洪水時(shí),沂沭泗水系下游連通伴生洪水風(fēng)險(xiǎn)如圖3(a)所示。沂沭河上游遭遇重現(xiàn)期為50年一遇洪水時(shí),南四湖控制下泄,沂河下游流量超過(guò)8000 m3/s時(shí),邳蒼分洪道進(jìn)行分洪,但分洪量較小,故韓莊運(yùn)河及中運(yùn)河段伴生洪水風(fēng)險(xiǎn)較低(P=0,R=0)。新沂河承接駱馬湖嶂山閘下泄洪水及老沭河來(lái)水,新沂河全段流量風(fēng)險(xiǎn)為中風(fēng)險(xiǎn)(P=0.55,R=0.3);大部分河段水位風(fēng)險(xiǎn)處于中風(fēng)險(xiǎn)(P=0.42,R=0.12),??诙尾糠?jǐn)嗝嫠伙L(fēng)險(xiǎn)達(dá)到高風(fēng)險(xiǎn)(P=1,R=0.26),原因是河口淤積造成河床升高;流速風(fēng)險(xiǎn)大部分處于中風(fēng)險(xiǎn)(P=0.54,R=0.12),其中口頭段為新沂河與老沭河匯流處,河床沖蝕較深,流速降低,因此流速低風(fēng)險(xiǎn)(P=0,R=0)。

    圖3 沂沭泗水系下游各風(fēng)險(xiǎn)源下風(fēng)險(xiǎn)等級(jí)分布

    沂沭河上游未遭遇洪水,南四湖流域遭遇重現(xiàn)期為50年一遇洪水時(shí),南四湖敞泄,最大下泄流量不超過(guò)5600 m3/s,沂河下游流量未超過(guò)8000 m3/s,邳蒼分洪道不分洪。由于南四湖下泄流量較小,故韓莊運(yùn)河以及中運(yùn)河段風(fēng)險(xiǎn)較低,駱馬湖總?cè)牒髁枯^小,新沂河段風(fēng)險(xiǎn)亦較低,如圖3(b)所示,沂沭泗水系下游整體處于低風(fēng)險(xiǎn)。

    當(dāng)沂沭泗流域全境遭遇重現(xiàn)期為50年一遇洪水時(shí),沂沭泗水系下游連通伴生風(fēng)險(xiǎn)如圖3(c)所示。韓莊運(yùn)河段處于低風(fēng)險(xiǎn),中運(yùn)河段由于疊加邳蒼分洪道來(lái)水,流量(P=0.5,R=0.25)、水位(P=0.41,R=0.05)、流速(P=0.5,R=0.11)均升至中風(fēng)險(xiǎn);新沂河由于承接駱馬湖下泄洪水及老沭河來(lái)水,全段流量風(fēng)險(xiǎn)升至高風(fēng)險(xiǎn)(P=0.8,R=0.59),水位風(fēng)險(xiǎn)升至高風(fēng)險(xiǎn)(P=0.73,R=0.16),流速風(fēng)險(xiǎn)大部分為高風(fēng)險(xiǎn)(P=0.75,R=0.24);海口段附近為中風(fēng)險(xiǎn)(P=0.65,R=0.27),口頭段低風(fēng)險(xiǎn)(P=0,R=0)。

    當(dāng)洪水重現(xiàn)期為百年一遇時(shí),各風(fēng)險(xiǎn)源洪水風(fēng)險(xiǎn)均有較大提升。對(duì)于風(fēng)險(xiǎn)源①,邳蒼分洪道在大部分時(shí)間行洪,但由于南四湖控制下泄,韓莊運(yùn)河及中運(yùn)河段洪水風(fēng)險(xiǎn)仍處于較低水平;駱馬湖入湖流量和嶂山閘下泄流量顯著增加,新沂河洪水風(fēng)險(xiǎn)顯著提升,由中風(fēng)險(xiǎn)提升至高風(fēng)險(xiǎn)(P=0.77,R=0.43)。對(duì)于風(fēng)險(xiǎn)源②,韓莊運(yùn)河由低風(fēng)險(xiǎn)升至高風(fēng)險(xiǎn)(P=0.56,R=0.41),中運(yùn)河段由低風(fēng)險(xiǎn)升至中風(fēng)險(xiǎn)(P=0.36,R=0.27),新沂河由無(wú)風(fēng)險(xiǎn)升至低風(fēng)險(xiǎn)(P=0.2,R=0.11)。對(duì)于風(fēng)險(xiǎn)源③,當(dāng)沂沭泗流域全境遭遇百年一遇洪水時(shí),沂沭泗水系下游洪水風(fēng)險(xiǎn)均有較大提升,中運(yùn)河由于同時(shí)承接韓莊運(yùn)河和邳蒼分洪道來(lái)水,洪水風(fēng)險(xiǎn)由中風(fēng)險(xiǎn)升至高風(fēng)險(xiǎn);駱馬湖入湖流量巨大,新沂河全段為高風(fēng)險(xiǎn)。

    4 沂沭泗水系下游河湖連通伴生洪水風(fēng)險(xiǎn)管控措施及效果

    由上節(jié)可知,風(fēng)險(xiǎn)源②條件下沂沭泗水系下游整體處于低風(fēng)險(xiǎn),故僅對(duì)風(fēng)險(xiǎn)源①及風(fēng)險(xiǎn)源③進(jìn)行水力調(diào)控,以降低河湖連通伴生洪水風(fēng)險(xiǎn)。圖4 為新沂河調(diào)控前后風(fēng)險(xiǎn)變化對(duì)比圖,因口頭段為新沂河與老沭河匯流點(diǎn),海口段為新沂河入海處,流速較低,調(diào)控前后流速風(fēng)險(xiǎn)概率和強(qiáng)度均接近0,故圖中未表示。新沂河沭陽(yáng)站匯集了駱馬湖下泄洪水及老沭河來(lái)水,洪水風(fēng)險(xiǎn)等級(jí)需重點(diǎn)關(guān)注,該控制站點(diǎn)不同風(fēng)險(xiǎn)源條件下調(diào)控方式及效果列于表4。當(dāng)僅沂沭河上游遭遇重現(xiàn)期為50年一遇洪水時(shí),駱馬湖起調(diào)水位由23.5 降至22.5 m,同時(shí)保證庫(kù)區(qū)汛末蓄水位保持在23 m。經(jīng)駱馬湖調(diào)蓄后,嶂山閘下泄流量減小,因僅降低起調(diào)水位,水庫(kù)削峰效果不明顯;新沂河流量風(fēng)險(xiǎn)概率降低10%左右,風(fēng)險(xiǎn)強(qiáng)度降低16%左右,全段仍處于中風(fēng)險(xiǎn)(P=0.5,R=0.25),如圖4(a)所示。

    表4 不同風(fēng)險(xiǎn)源下駱馬湖調(diào)控方式及效果

    當(dāng)沂沭泗流域全境遭遇重現(xiàn)期為50年一遇洪水時(shí),駱馬湖起調(diào)水位由23.5 降至22.5 m 并進(jìn)行預(yù)泄,盡量使湖區(qū)騰出空間來(lái)承接沂河下游及中運(yùn)河來(lái)水,同時(shí)湖區(qū)汛末蓄水位保持在23 m。調(diào)蓄后,水庫(kù)削峰效果明顯,嶂山閘至口頭流量風(fēng)險(xiǎn)概率由0.63 降至0.5 左右,口頭至??诙瘟髁匡L(fēng)險(xiǎn)概率從0.8 降至0.5,全段風(fēng)險(xiǎn)強(qiáng)度下降約80%,風(fēng)險(xiǎn)等級(jí)由高風(fēng)險(xiǎn)降為中風(fēng)險(xiǎn),如圖4(b)所示。

    當(dāng)風(fēng)險(xiǎn)源①及風(fēng)險(xiǎn)源③洪水重現(xiàn)期提升至百年一遇時(shí),駱馬湖起調(diào)水位由23.5 降至22.5 m 并進(jìn)行預(yù)泄,同時(shí)啟用黃墩滯洪區(qū)進(jìn)行滯洪。對(duì)于風(fēng)險(xiǎn)源①,經(jīng)調(diào)控后,新沂河沭陽(yáng)站洪水風(fēng)險(xiǎn)概率由0.8 降至0.6,風(fēng)險(xiǎn)強(qiáng)度由0.42 降至0.26,由高風(fēng)險(xiǎn)降至中風(fēng)險(xiǎn),如圖4(c)所示;而對(duì)于風(fēng)險(xiǎn)源③,由于駱馬湖入湖流量巨大,遠(yuǎn)超駱馬湖調(diào)蓄能力,調(diào)控幾乎不起作用。

    圖4 新沂河調(diào)控前后風(fēng)險(xiǎn)對(duì)比

    5 水力調(diào)控對(duì)新沂河洪水風(fēng)險(xiǎn)影響規(guī)律

    選取新沂河同一斷面不同水情條件下的水深(本節(jié)中將上文水位變換為水深以消除不同斷面河床高程差異)與流速模擬值,計(jì)算其中位數(shù)作為該斷面水深與流速的代表值,繪制圖5 中散點(diǎn)。文中定義,水深與流速風(fēng)險(xiǎn)概率P 均小于0.4 且風(fēng)險(xiǎn)強(qiáng)度R 均小于0.2 的斷面為安全斷面,水深與流速兩者任一P 大于0.4 或R 大于0.2 的斷面為風(fēng)險(xiǎn)斷面,經(jīng)調(diào)控后由風(fēng)險(xiǎn)斷面變?yōu)榘踩珨嗝娴姆Q之為風(fēng)險(xiǎn)緩釋斷面;安全區(qū)為水深與流速均處于低風(fēng)險(xiǎn)時(shí)(P≤0.4,R≤0.2)對(duì)應(yīng)的區(qū)域;水深與流速的安全閾值分別為低風(fēng)險(xiǎn)上界P=0.4 或R=0.2 對(duì)應(yīng)的數(shù)值,即調(diào)控前水深安全閾值為9.7 m、流速為0.8 m/s,調(diào)控后分別為10.3 m 和0.93 m/s。研究發(fā)現(xiàn),通過(guò)駱馬湖水力調(diào)控,下游河道部分?jǐn)嗝婷撾x風(fēng)險(xiǎn)區(qū),進(jìn)入風(fēng)險(xiǎn)區(qū)與調(diào)控前安全區(qū)之間,文中將其稱為緩沖區(qū)。50年一遇重現(xiàn)期洪水時(shí),綜合風(fēng)險(xiǎn)源①和③,若駱馬湖未進(jìn)行調(diào)控,當(dāng)入湖流量達(dá)到11 046 m3/s 時(shí),下游新沂河河道超過(guò)62%的斷面出現(xiàn)風(fēng)險(xiǎn),為風(fēng)險(xiǎn)斷面;調(diào)控后當(dāng)入湖流量達(dá)到12 672.6 m3/s 時(shí),下游河道仍有約80%的斷面洪水風(fēng)險(xiǎn)為低風(fēng)險(xiǎn),新沂河近42%的斷面進(jìn)入緩沖區(qū),如圖5(a)所示。當(dāng)遭遇百年一遇重現(xiàn)期洪水時(shí),風(fēng)險(xiǎn)源③調(diào)控?zé)o效果,僅考慮風(fēng)險(xiǎn)源①,如圖5(b)所示。若駱馬湖未進(jìn)行調(diào)控,新沂河72%的斷面為風(fēng)險(xiǎn)斷面,調(diào)控后當(dāng)入湖流量達(dá)到12 923 m3/s 時(shí),下游河道斷面風(fēng)險(xiǎn)概率將達(dá)到0.4(風(fēng)險(xiǎn)強(qiáng)度仍小于0.2),新沂河28%的斷面進(jìn)入緩沖區(qū),處于可容忍風(fēng)險(xiǎn)區(qū)域。

    圖5 駱馬湖水力調(diào)控對(duì)新沂河斷面風(fēng)險(xiǎn)影響規(guī)律

    6 湖泊水力調(diào)控的風(fēng)險(xiǎn)緩釋作用

    目前,多數(shù)水系連通風(fēng)險(xiǎn)評(píng)估研究只針對(duì)河河連通,或在河湖連通研究中未考慮湖泊的主動(dòng)調(diào)控作用,相關(guān)研究成果列于表5。如利用圖論相關(guān)方法研究不同連通方案下水系連通程度以及定量評(píng)估水系連通洪水風(fēng)險(xiǎn),選擇區(qū)域河網(wǎng)最優(yōu)連通方案[18];在沂沭河上游的河河連通分析中,以水深、流速風(fēng)險(xiǎn)限值定量計(jì)算水系連通伴生洪水風(fēng)險(xiǎn)概率[12]。在太湖流域的洪水風(fēng)險(xiǎn)研究中,通過(guò)建立大尺度水動(dòng)力模型,對(duì)洪水災(zāi)害進(jìn)行評(píng)估,該研究中加入了湖泊單元,但未研究湖泊在水系連通工程中的主動(dòng)調(diào)控作用[19];采用一二維耦合水動(dòng)力模型的方法分析洞庭湖流域的洪水風(fēng)險(xiǎn),考慮了湖泊的調(diào)控能力,但并未對(duì)其作出定量分析[20]。

    表5 水系連通風(fēng)險(xiǎn)評(píng)估不同研究結(jié)果對(duì)比

    本文通過(guò)湖泊水力調(diào)控,分析發(fā)現(xiàn)河湖連通所伴生的洪水風(fēng)險(xiǎn)存在緩沖區(qū)域。該區(qū)域通過(guò)一定的湖泊調(diào)控方式,提高了下游河道水深、流速安全閾值,突破了原有河河連通研究中單一風(fēng)險(xiǎn)安全閾值的限制,使得風(fēng)險(xiǎn)管控更加靈活。例如,本文通過(guò)湖泊調(diào)控,將下游河道水深安全閾值由9.7 提升至10.3 m,提高6%,流速安全閾值由0.8 提升至0.93 m/s,提高16%。因此,湖泊作為河湖連通工程中的大型儲(chǔ)蓄水單元,對(duì)水系連通工程伴生風(fēng)險(xiǎn)的調(diào)節(jié)、緩釋作用不可忽略。

    7 結(jié)論及建議

    沂沭泗水系下游遭遇50年一遇洪水時(shí),部分河段面臨較高風(fēng)險(xiǎn),利用駱馬湖進(jìn)行水力調(diào)控時(shí),新沂河42%的風(fēng)險(xiǎn)斷面進(jìn)入緩沖區(qū),河湖連通伴生洪水風(fēng)險(xiǎn)降至低風(fēng)險(xiǎn)(P=0.4,R=0.2)。但可預(yù)見(jiàn)的是,駱馬湖調(diào)控能力有限,在遭遇更大重現(xiàn)期洪水時(shí),沂沭泗流域仍面臨著較高風(fēng)險(xiǎn),危及兩岸人民的生命及財(cái)產(chǎn)安全。為此,結(jié)合本文研究成果對(duì)沂沭泗水系下游防洪調(diào)度提出以下建議:(1)駱馬湖起調(diào)水位由23.5 降至22.5 m,并控制駱馬湖入湖流量小于12 923 m3/s,使沂沭泗水系下游洪水風(fēng)險(xiǎn)處于風(fēng)險(xiǎn)緩沖區(qū)域內(nèi)。(2)擴(kuò)大新沂河入海規(guī)模,沂沭河上游來(lái)水優(yōu)先東調(diào),減小沂河下游進(jìn)入駱馬湖流量,同時(shí)盡量避免使用邳蒼分洪道分洪,減輕駱馬湖及中運(yùn)河的防洪壓力。(3)提升中運(yùn)河校核流量至10 000 m3/s,確保南四湖流域遭遇50年一遇洪水時(shí),南四湖超額洪水可以敞泄至韓莊運(yùn)河,保障南四湖防洪安全,同時(shí)在邳蒼分洪道進(jìn)行分洪時(shí),中運(yùn)河發(fā)生洪水的風(fēng)險(xiǎn)概率能夠降至低風(fēng)險(xiǎn)。(4)當(dāng)駱馬湖承接上游超50年一遇洪水時(shí),提升嶂山閘泄洪能力,盡量避免使用黃墩滯洪區(qū)進(jìn)行滯洪,以免造成巨大的經(jīng)濟(jì)損失。同時(shí)提升新沂河入海規(guī)模,減小新沂河口頭至??诙蔚暮樗L(fēng)險(xiǎn)。

    猜你喜歡
    沭河駱馬湖河湖
    美麗的沭河公園
    全面推行河湖長(zhǎng)制打造“三晉”幸福河湖
    山西水利(2022年5期)2022-09-21 02:38:28
    淮河流域省級(jí)河湖長(zhǎng)第一次聯(lián)席會(huì)議召開(kāi)
    治淮(2022年8期)2022-09-03 03:42:26
    我家門(mén)前沭水流
    加強(qiáng)河湖保護(hù)治理 改善河湖生態(tài)環(huán)境
    蘆 葦
    黑龍江:河湖治理保護(hù)成效明顯
    近500年來(lái)駱馬湖演變的驅(qū)動(dòng)力探究
    美麗的駱馬湖
    問(wèn)詢駱馬湖(外一章)
    久久久国产欧美日韩av| 国产精品一区二区在线观看99| 一区二区三区激情视频| 激情在线观看视频在线高清 | 丰满饥渴人妻一区二区三| 国产国语露脸激情在线看| 老司机影院毛片| 丝瓜视频免费看黄片| 免费日韩欧美在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美一区二区三区黑人| 激情在线观看视频在线高清 | 99久久国产精品久久久| 香蕉丝袜av| 久久久久久久大尺度免费视频| 欧美乱码精品一区二区三区| 亚洲伊人久久精品综合| 成在线人永久免费视频| 日本欧美视频一区| 国产日韩欧美视频二区| 日本av手机在线免费观看| 国产国语露脸激情在线看| 国内毛片毛片毛片毛片毛片| 1024香蕉在线观看| 亚洲成人免费电影在线观看| 国产aⅴ精品一区二区三区波| 99精品久久久久人妻精品| 1024视频免费在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲av成人一区二区三| 国产精品亚洲av一区麻豆| 精品福利永久在线观看| 日韩欧美三级三区| 久久久久久人人人人人| 色播在线永久视频| 亚洲综合色网址| 18在线观看网站| 亚洲男人天堂网一区| 精品福利观看| 啦啦啦视频在线资源免费观看| 无人区码免费观看不卡 | 多毛熟女@视频| 亚洲熟女精品中文字幕| 一进一出好大好爽视频| 免费在线观看影片大全网站| 亚洲午夜理论影院| 亚洲精品国产色婷婷电影| 国产精品偷伦视频观看了| 亚洲精品美女久久av网站| 午夜免费成人在线视频| 丰满少妇做爰视频| 十分钟在线观看高清视频www| 亚洲精品中文字幕一二三四区 | 桃红色精品国产亚洲av| 亚洲av电影在线进入| 精品一品国产午夜福利视频| 久久久久久久大尺度免费视频| 亚洲va日本ⅴa欧美va伊人久久| 久久影院123| 亚洲伊人久久精品综合| 免费人妻精品一区二区三区视频| 久久国产精品大桥未久av| 亚洲成人国产一区在线观看| 搡老岳熟女国产| 欧美+亚洲+日韩+国产| 一级黄色大片毛片| 亚洲av片天天在线观看| 久久久精品94久久精品| 夜夜夜夜夜久久久久| 久久久久久久精品吃奶| 岛国毛片在线播放| 久热这里只有精品99| 中文字幕制服av| 国产精品99久久99久久久不卡| 人人妻人人添人人爽欧美一区卜| 最黄视频免费看| 高清毛片免费观看视频网站 | 最近最新中文字幕大全免费视频| 欧美老熟妇乱子伦牲交| 三上悠亚av全集在线观看| 成年人午夜在线观看视频| 窝窝影院91人妻| 亚洲国产中文字幕在线视频| 欧美人与性动交α欧美精品济南到| 天堂中文最新版在线下载| 欧美黄色片欧美黄色片| 国产精品久久久av美女十八| 黄色 视频免费看| 一边摸一边抽搐一进一出视频| 99国产精品99久久久久| svipshipincom国产片| 99热网站在线观看| 午夜精品久久久久久毛片777| 电影成人av| 成人免费观看视频高清| 午夜福利乱码中文字幕| 男人舔女人的私密视频| 精品少妇一区二区三区视频日本电影| 757午夜福利合集在线观看| 窝窝影院91人妻| 国产亚洲欧美精品永久| 99热网站在线观看| 香蕉国产在线看| 久久天躁狠狠躁夜夜2o2o| 国产伦人伦偷精品视频| 国内毛片毛片毛片毛片毛片| 亚洲精品美女久久av网站| 一区二区日韩欧美中文字幕| 午夜免费成人在线视频| 十八禁高潮呻吟视频| 超碰97精品在线观看| 国产视频一区二区在线看| 水蜜桃什么品种好| 天天躁日日躁夜夜躁夜夜| 精品国内亚洲2022精品成人 | 99国产精品免费福利视频| 午夜福利乱码中文字幕| 亚洲伊人色综图| 亚洲国产成人一精品久久久| 美女国产高潮福利片在线看| 国产成人精品久久二区二区免费| 亚洲熟妇熟女久久| 亚洲欧美激情在线| 亚洲va日本ⅴa欧美va伊人久久| 黄色a级毛片大全视频| 国产激情久久老熟女| 国产97色在线日韩免费| 久久久久国产一级毛片高清牌| 别揉我奶头~嗯~啊~动态视频| 欧美激情高清一区二区三区| 国产精品二区激情视频| avwww免费| 国产一区有黄有色的免费视频| www.自偷自拍.com| 成人av一区二区三区在线看| 侵犯人妻中文字幕一二三四区| 热99国产精品久久久久久7| 久久人妻av系列| 欧美大码av| 久久久久精品国产欧美久久久| 久久久久国产一级毛片高清牌| 成人国产av品久久久| 久久人妻av系列| 一级毛片电影观看| 欧美成人免费av一区二区三区 | 18禁国产床啪视频网站| 成人国产一区最新在线观看| 国产无遮挡羞羞视频在线观看| 黄色丝袜av网址大全| 熟女少妇亚洲综合色aaa.| 三级毛片av免费| 欧美人与性动交α欧美精品济南到| 免费人妻精品一区二区三区视频| 久久精品成人免费网站| 久久精品国产亚洲av高清一级| 人人妻人人澡人人爽人人夜夜| 啦啦啦 在线观看视频| 亚洲av美国av| 久久精品国产亚洲av香蕉五月 | 12—13女人毛片做爰片一| 悠悠久久av| 12—13女人毛片做爰片一| 久久久精品区二区三区| 中文亚洲av片在线观看爽 | 一本大道久久a久久精品| 咕卡用的链子| 精品国产国语对白av| 亚洲欧美一区二区三区久久| 午夜福利影视在线免费观看| 欧美日韩国产mv在线观看视频| 天天躁日日躁夜夜躁夜夜| 一本久久精品| 欧美国产精品va在线观看不卡| 免费日韩欧美在线观看| 在线十欧美十亚洲十日本专区| 亚洲va日本ⅴa欧美va伊人久久| 99国产精品一区二区三区| 午夜福利影视在线免费观看| 日韩免费高清中文字幕av| 国产精品久久久久久精品电影小说| 久久久久网色| 欧美国产精品va在线观看不卡| 欧美国产精品va在线观看不卡| 久久久久网色| 亚洲色图av天堂| 黄色成人免费大全| 国产成人欧美| 亚洲av成人不卡在线观看播放网| 亚洲国产成人一精品久久久| 伊人久久大香线蕉亚洲五| 精品亚洲成国产av| 在线观看一区二区三区激情| 欧美黑人精品巨大| 国产不卡av网站在线观看| 亚洲成av片中文字幕在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 最新在线观看一区二区三区| 少妇的丰满在线观看| 中文字幕av电影在线播放| 亚洲精品国产色婷婷电影| 亚洲成a人片在线一区二区| 久久精品91无色码中文字幕| 国产片内射在线| 不卡av一区二区三区| 99久久人妻综合| 久久国产精品大桥未久av| 亚洲伊人久久精品综合| 少妇猛男粗大的猛烈进出视频| 日韩免费av在线播放| 黑丝袜美女国产一区| 日韩大片免费观看网站| 国产精品熟女久久久久浪| 日韩 欧美 亚洲 中文字幕| 国产成人精品久久二区二区91| 欧美日韩亚洲国产一区二区在线观看 | 啦啦啦中文免费视频观看日本| 建设人人有责人人尽责人人享有的| 一区二区日韩欧美中文字幕| 高清黄色对白视频在线免费看| 在线av久久热| 丝袜美足系列| 亚洲三区欧美一区| 久热这里只有精品99| 亚洲七黄色美女视频| 色婷婷av一区二区三区视频| 中文字幕高清在线视频| 精品人妻1区二区| 91av网站免费观看| 夜夜夜夜夜久久久久| 美女福利国产在线| 亚洲专区中文字幕在线| 99国产精品免费福利视频| cao死你这个sao货| 国产成人欧美在线观看 | 一本—道久久a久久精品蜜桃钙片| www.熟女人妻精品国产| 在线观看免费高清a一片| 曰老女人黄片| 少妇粗大呻吟视频| 啦啦啦免费观看视频1| 51午夜福利影视在线观看| 最新在线观看一区二区三区| av国产精品久久久久影院| 麻豆乱淫一区二区| 久久亚洲精品不卡| 亚洲精品美女久久av网站| 国产精品九九99| 国产精品偷伦视频观看了| 91大片在线观看| 青青草视频在线视频观看| 久久热在线av| 正在播放国产对白刺激| 成人影院久久| 欧美日韩av久久| 飞空精品影院首页| 精品一品国产午夜福利视频| 亚洲精品美女久久av网站| 男女之事视频高清在线观看| 久久午夜亚洲精品久久| 亚洲精品乱久久久久久| 国产成人av激情在线播放| 欧美乱码精品一区二区三区| 亚洲精品成人av观看孕妇| 新久久久久国产一级毛片| 亚洲av美国av| 亚洲欧美一区二区三区黑人| 男人操女人黄网站| 亚洲国产欧美在线一区| 国精品久久久久久国模美| 国产亚洲一区二区精品| 欧美黑人欧美精品刺激| 黄片大片在线免费观看| 亚洲国产中文字幕在线视频| 欧美 亚洲 国产 日韩一| 午夜免费成人在线视频| 91av网站免费观看| 成人国产av品久久久| 国产视频一区二区在线看| 欧美日韩亚洲高清精品| 国产在线视频一区二区| av网站免费在线观看视频| 一本色道久久久久久精品综合| 亚洲精品av麻豆狂野| 国产无遮挡羞羞视频在线观看| 91麻豆精品激情在线观看国产 | 国产精品亚洲一级av第二区| 蜜桃在线观看..| 桃红色精品国产亚洲av| 国产日韩一区二区三区精品不卡| 一级a爱视频在线免费观看| 久热这里只有精品99| 丝袜人妻中文字幕| 国产人伦9x9x在线观看| 97人妻天天添夜夜摸| av欧美777| 欧美激情高清一区二区三区| 一本色道久久久久久精品综合| 97人妻天天添夜夜摸| 国产精品免费一区二区三区在线 | 男人舔女人的私密视频| 美女高潮到喷水免费观看| 99香蕉大伊视频| 亚洲色图综合在线观看| 男人舔女人的私密视频| 91精品三级在线观看| 中亚洲国语对白在线视频| 婷婷丁香在线五月| 淫妇啪啪啪对白视频| 婷婷成人精品国产| 国产不卡一卡二| 99在线人妻在线中文字幕 | 少妇裸体淫交视频免费看高清 | 极品少妇高潮喷水抽搐| avwww免费| 99精品欧美一区二区三区四区| 国产福利在线免费观看视频| e午夜精品久久久久久久| 丁香六月欧美| 成年人午夜在线观看视频| 亚洲国产成人一精品久久久| 国产精品.久久久| 夜夜夜夜夜久久久久| 色婷婷久久久亚洲欧美| 91麻豆精品激情在线观看国产 | 午夜激情av网站| 国产亚洲精品第一综合不卡| 肉色欧美久久久久久久蜜桃| 欧美成人免费av一区二区三区 | 狠狠狠狠99中文字幕| 欧美成人免费av一区二区三区 | 中文字幕色久视频| 国产精品熟女久久久久浪| 亚洲国产毛片av蜜桃av| 成人亚洲精品一区在线观看| 菩萨蛮人人尽说江南好唐韦庄| 久久国产亚洲av麻豆专区| 丝瓜视频免费看黄片| 母亲3免费完整高清在线观看| 人妻久久中文字幕网| 97人妻天天添夜夜摸| 欧美黑人精品巨大| 久久精品亚洲av国产电影网| 久久影院123| 中文字幕制服av| 国产三级黄色录像| 日韩有码中文字幕| xxxhd国产人妻xxx| 老司机靠b影院| 国产精品九九99| 一本—道久久a久久精品蜜桃钙片| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品在线观看二区| 性色av乱码一区二区三区2| av一本久久久久| 亚洲色图av天堂| 18禁国产床啪视频网站| 国产一区二区激情短视频| 亚洲人成电影免费在线| 大片免费播放器 马上看| 成人18禁在线播放| 丰满少妇做爰视频| 久久精品国产a三级三级三级| 国产不卡一卡二| 久久人妻福利社区极品人妻图片| 免费不卡黄色视频| netflix在线观看网站| 国产精品久久久久久精品古装| 亚洲第一青青草原| 麻豆国产av国片精品| 亚洲av成人一区二区三| 欧美在线黄色| 欧美亚洲 丝袜 人妻 在线| 黄色毛片三级朝国网站| 自线自在国产av| 一边摸一边抽搐一进一小说 | 99热网站在线观看| 免费在线观看视频国产中文字幕亚洲| 午夜福利一区二区在线看| 一区在线观看完整版| 国产老妇伦熟女老妇高清| 亚洲久久久国产精品| 法律面前人人平等表现在哪些方面| 女人被躁到高潮嗷嗷叫费观| 少妇裸体淫交视频免费看高清 | 国产主播在线观看一区二区| 交换朋友夫妻互换小说| 国产精品免费大片| 女性被躁到高潮视频| 美女午夜性视频免费| 国产免费福利视频在线观看| 亚洲第一欧美日韩一区二区三区 | 蜜桃国产av成人99| 日韩欧美三级三区| 桃红色精品国产亚洲av| 91老司机精品| 每晚都被弄得嗷嗷叫到高潮| 国产三级黄色录像| 无人区码免费观看不卡 | 国产免费av片在线观看野外av| 又大又爽又粗| 黑人欧美特级aaaaaa片| 美国免费a级毛片| 精品亚洲成国产av| 亚洲国产精品一区二区三区在线| 精品欧美一区二区三区在线| 女性被躁到高潮视频| 久久亚洲精品不卡| 午夜成年电影在线免费观看| 国产亚洲精品第一综合不卡| 99精国产麻豆久久婷婷| 精品国产一区二区久久| 咕卡用的链子| 丁香欧美五月| 交换朋友夫妻互换小说| 丰满迷人的少妇在线观看| 亚洲精品国产区一区二| 国产欧美亚洲国产| 国产麻豆69| 免费日韩欧美在线观看| 国产淫语在线视频| 在线观看免费日韩欧美大片| 亚洲va日本ⅴa欧美va伊人久久| 天堂动漫精品| 制服诱惑二区| 91精品国产国语对白视频| 老熟女久久久| 精品一区二区三卡| 日韩人妻精品一区2区三区| 久久人妻av系列| 妹子高潮喷水视频| 日韩三级视频一区二区三区| 亚洲人成伊人成综合网2020| 精品少妇内射三级| 日韩 欧美 亚洲 中文字幕| 久久免费观看电影| 久久青草综合色| 男女午夜视频在线观看| 久久久久久久久久久久大奶| 精品一区二区三区av网在线观看 | 国产亚洲精品第一综合不卡| 天堂中文最新版在线下载| 国产精品久久电影中文字幕 | 18禁裸乳无遮挡动漫免费视频| 热99国产精品久久久久久7| av欧美777| 两性夫妻黄色片| 男男h啪啪无遮挡| 亚洲精品在线观看二区| 亚洲av日韩在线播放| 正在播放国产对白刺激| 12—13女人毛片做爰片一| 99精国产麻豆久久婷婷| 一本色道久久久久久精品综合| 黄色视频在线播放观看不卡| 一本一本久久a久久精品综合妖精| 日本撒尿小便嘘嘘汇集6| 在线十欧美十亚洲十日本专区| 他把我摸到了高潮在线观看 | 久久精品亚洲熟妇少妇任你| 国产主播在线观看一区二区| 青草久久国产| 丰满少妇做爰视频| 亚洲第一av免费看| 汤姆久久久久久久影院中文字幕| 亚洲精品久久午夜乱码| 91成人精品电影| 中国美女看黄片| 美女福利国产在线| 99re在线观看精品视频| 大香蕉久久成人网| 色综合婷婷激情| h视频一区二区三区| 制服诱惑二区| 亚洲自偷自拍图片 自拍| 欧美午夜高清在线| 三级毛片av免费| 黄色成人免费大全| 亚洲中文av在线| 欧美一级毛片孕妇| 热re99久久精品国产66热6| 国产男女超爽视频在线观看| 亚洲欧洲日产国产| 不卡一级毛片| 性色av乱码一区二区三区2| svipshipincom国产片| 狠狠精品人妻久久久久久综合| 亚洲国产精品一区二区三区在线| 亚洲欧美日韩高清在线视频 | 18禁黄网站禁片午夜丰满| 下体分泌物呈黄色| 亚洲精品粉嫩美女一区| 大码成人一级视频| 精品国产超薄肉色丝袜足j| 少妇裸体淫交视频免费看高清 | 性色av乱码一区二区三区2| 国产精品一区二区免费欧美| av天堂在线播放| a在线观看视频网站| 欧美日韩国产mv在线观看视频| 国产精品国产av在线观看| 亚洲国产欧美一区二区综合| 亚洲精品av麻豆狂野| 国产主播在线观看一区二区| 女警被强在线播放| 国产精品一区二区精品视频观看| 成人亚洲精品一区在线观看| 免费久久久久久久精品成人欧美视频| 三上悠亚av全集在线观看| 69av精品久久久久久 | 他把我摸到了高潮在线观看 | 美女视频免费永久观看网站| 日本a在线网址| 午夜免费鲁丝| 久久亚洲精品不卡| 亚洲人成伊人成综合网2020| 国产日韩欧美视频二区| 一夜夜www| 91精品国产国语对白视频| 久久亚洲精品不卡| 两人在一起打扑克的视频| 日本五十路高清| 久久这里只有精品19| 成人特级黄色片久久久久久久 | 男女边摸边吃奶| 黑人操中国人逼视频| 十八禁网站网址无遮挡| 国产精品久久电影中文字幕 | 最近最新中文字幕大全电影3 | 十八禁网站免费在线| 熟女少妇亚洲综合色aaa.| 国产一区二区三区视频了| 欧美激情高清一区二区三区| 成人精品一区二区免费| 欧美老熟妇乱子伦牲交| 久久久久国产一级毛片高清牌| 999久久久国产精品视频| 一级,二级,三级黄色视频| 搡老熟女国产l中国老女人| 午夜福利,免费看| 人妻久久中文字幕网| 97在线人人人人妻| 亚洲国产av影院在线观看| 精品国产亚洲在线| 亚洲av欧美aⅴ国产| 亚洲成人手机| 99在线人妻在线中文字幕 | 少妇的丰满在线观看| 亚洲伊人色综图| av福利片在线| 蜜桃国产av成人99| 一级片免费观看大全| 亚洲专区字幕在线| av有码第一页| xxxhd国产人妻xxx| 日韩三级视频一区二区三区| 中文字幕最新亚洲高清| 成年动漫av网址| 久久精品国产99精品国产亚洲性色 | 男女床上黄色一级片免费看| 久久久久精品国产欧美久久久| 国产不卡一卡二| av又黄又爽大尺度在线免费看| 日韩制服丝袜自拍偷拍| 免费观看a级毛片全部| 国产亚洲一区二区精品| 一级,二级,三级黄色视频| 热99国产精品久久久久久7| 叶爱在线成人免费视频播放| 窝窝影院91人妻| 高清视频免费观看一区二区| 高清在线国产一区| 男女无遮挡免费网站观看| av又黄又爽大尺度在线免费看| 亚洲国产毛片av蜜桃av| 欧美黄色淫秽网站| tube8黄色片| 久久精品国产99精品国产亚洲性色 | 日韩视频一区二区在线观看| 国精品久久久久久国模美| 十八禁网站免费在线| 亚洲精品久久午夜乱码| 亚洲av美国av| 成人18禁高潮啪啪吃奶动态图| 久久av网站| 黄片播放在线免费| 99re6热这里在线精品视频| 欧美久久黑人一区二区| 91麻豆精品激情在线观看国产 | 老熟女久久久| 日本vs欧美在线观看视频| 999久久久精品免费观看国产| 欧美亚洲 丝袜 人妻 在线| 国产精品一区二区精品视频观看| 亚洲黑人精品在线| 老熟妇乱子伦视频在线观看| 亚洲黑人精品在线| 在线十欧美十亚洲十日本专区| 精品免费久久久久久久清纯 | 久久99热这里只频精品6学生| 亚洲成av片中文字幕在线观看| 黄片播放在线免费| 中文字幕人妻丝袜一区二区| 精品一区二区三区av网在线观看 | 满18在线观看网站| 五月天丁香电影| 男女之事视频高清在线观看| 免费在线观看视频国产中文字幕亚洲| 中文字幕另类日韩欧美亚洲嫩草| 成人永久免费在线观看视频 | 精品国内亚洲2022精品成人 | 亚洲精品成人av观看孕妇| 中国美女看黄片| 一边摸一边做爽爽视频免费| 亚洲av国产av综合av卡| 丰满饥渴人妻一区二区三|