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

    基于SWMM模型的LID設施的雨洪控制效果分析

    2019-06-25 08:57:20李思祎張建豐計寶鑫
    中國農(nóng)村水利水電 2019年6期
    關鍵詞:排放口徑流系數(shù)過程線

    李思祎,張建豐,李 濤,計寶鑫

    (西安理工大學水利水電學院,西安 710048)

    隨著高速的城市化,城市的不透水路面日益增加。近幾十年來,西安市城區(qū)內(nèi)澇災害頻繁,城市內(nèi)澇問題得到廣泛關注。如今,城市居民生產(chǎn)和生活中受內(nèi)澇災害影響嚴重[1,2]。研究表明,暴雨洪水管理模型SWMM(Storm Water Management Model)是目前研究較為深入,在國內(nèi)外應用最為廣泛的降雨-徑流模擬模型[3]。SWMM模型已經(jīng)更新至SWMM5.0版本,較之前增加了對低影響開發(fā)措施的支持,包括生物滯留池、連續(xù)多孔透水鋪裝路面、雨水桶、入滲池、植草溝5種措施。目前SWMM模型成為國內(nèi)外最主要的低影響開發(fā)措施的計算模型之一[4]。

    國內(nèi)學者運用SWMM模型在城市防洪排澇方面做了大量深入的研究,如熊麗君等選取上海市中心城區(qū),應用SWMM模型分析研究了該研究區(qū)域的降雨和產(chǎn)流特征,得出低影響開發(fā)措施對削減小雨量或低強度降雨下的徑流效果明顯,同時建立該地區(qū)徑流量與降雨因子的回歸方程,為該地區(qū)低影響開發(fā)方案提供參考[5]。柳浩林等運用SWMM模型分析兩種不同LID措施透水磚和下沉式綠地對降雨徑流過程的影響,得出LID措施對雨洪控制效果顯著[6]。目前,國內(nèi)在應用SWMM模型研究模擬中還存在若干問題,模型參數(shù)選取、加設LID措施時照搬國外案例經(jīng)驗不符合國內(nèi)實情等[7]。

    本文將準確的排水管網(wǎng)和下墊面微地形信息結(jié)合,劃分若干個降雨徑流匯水區(qū)分區(qū)?;赟WMM模型探求每一個獨立排水管網(wǎng)系統(tǒng)所在下墊面的降雨徑流特征規(guī)律和加設4種不同LID組合措施對雨洪的控制效果和規(guī)律。尤其針對城市化程度高,不透水率大的城區(qū)的防洪管理,提供改進建議。同時對西安市以及國內(nèi)其他城市的雨水管理以及內(nèi)澇防治工作提供一些理論參考與依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域概況

    以西安市南二環(huán)與環(huán)城南路之間、雁塔北路以西、長安北路以東為研究對象。該研究區(qū)域多年平均氣溫為13.3 ℃,在1月份氣溫最低,平均在-0.5~1.3 ℃,8月氣溫最高,平均氣溫為26.4~26.9 ℃。多年平均降雨量629.9 mm,7-10月的降水占全年降水量的55.9%。該區(qū)域降雨時空分布不均,夏季容易出現(xiàn)暴雨天氣,因此加重了城市內(nèi)澇發(fā)生的頻率。

    該區(qū)域建筑基本建于20世紀八九十年代,城市化水平高,建筑屋面面積占該區(qū)域總面積41.8%,不透水路面占該區(qū)域面積44.8%,不透水率高達86.6%,排水系統(tǒng)相對較老,極易造成城市內(nèi)澇災害,是西安市內(nèi)澇頻發(fā)的區(qū)域之一。研究區(qū)域面積211.25 hm2,區(qū)域大體地勢是西北高,東南低,該雨水管網(wǎng)系統(tǒng)沿文藝路截流建設路、友誼路、建西街等管道系統(tǒng)的雨水,最終排入護城河。

    1.2 研究區(qū)域概化

    根據(jù)西安市最新的排水管網(wǎng)分布圖,對子匯水區(qū)進行劃分,首先根據(jù)研究區(qū)域雨水管網(wǎng)圖以及雨水井的位置,對研究區(qū)域的主干道進行劃分。其次根據(jù)該區(qū)域的土地利用類型圖、數(shù)字高程圖,對已經(jīng)粗略劃分好的子匯水區(qū)進行詳細的劃分。這樣可以盡量減小模型模擬的誤差,使得劃分好的每一個子匯水區(qū)域真實的降雨徑流過程與模型模擬的降雨徑流過程保持相似。如圖1所示,劃分結(jié)果為:共22個子匯水區(qū),28個匯水節(jié)點,27段排水管道和1個排水總出口,圖中每個黑點代表一個子匯水區(qū)中心。

    圖1 子匯水區(qū)劃分與排水管網(wǎng)圖Fig.1 Sub-catchment area division and drainage pipe network diagram

    1.3 模型參數(shù)確定

    1.3.1 子匯水區(qū)特征參數(shù)的選取

    通過西安市地形圖,提取需要輸入模型中各子匯水區(qū)的具體參數(shù),如表1所示。

    1.3.2 水力參數(shù)選取

    根據(jù)模型的用戶手冊和相關文獻[8,9]中,提供的模型參數(shù)的選取辦法,如表2所示,對水力參數(shù)進行選取。

    1.3.3 模型驗證

    對上述選取的水力參數(shù),需要進行校驗,但在本研究中,研究區(qū)域?qū)儆诜比A的城市街道,并且排放口的位置處于地下數(shù)十米,現(xiàn)狀排放口處的實際流量過程線不易測得。所以,采用一種無實測資料的模型驗證方法,具體方法是參照實際的下墊面的情況,依據(jù)《室外排水設計規(guī)范》對不同下墊面的徑流參數(shù)進行取值。再計算不同下墊面的面積占研究區(qū)域總面積的面積權(quán)重,計算得到該區(qū)域的實際綜合徑流系數(shù)。最后通過模型模擬確定的徑流系數(shù)與實際綜合徑流系數(shù)進行比較,檢驗模型的誤差。表3是不同下墊面徑流系數(shù)取值范圍表。表4是城市地區(qū)不同區(qū)域的綜合徑流系數(shù)取值范圍表。

    表1 研究區(qū)域基本特征參數(shù)表Tab.1 Table of characteristic parameters of sub-catchment area

    表2 SWMM模型中主要水力參數(shù)取值范圍表Tab.2 The value range of main hydraulic parameters in SWMM model is shown in table

    表3 不同下墊面徑流系數(shù)取值范圍表Tab.3 A variety of underlying surface runoff coefficient value table

    表4 不同城市區(qū)域綜合徑流系數(shù)取值范圍表Tab.4 Urban area comprehensive runoff coefficient values range table

    因為研究區(qū)域?qū)儆谖靼彩卸h(huán)以內(nèi)的老城區(qū),該區(qū)域中底層建筑較多,城市化程度高,屬于城市建筑密集區(qū)。所以依據(jù)表4所示,綜合徑流系數(shù)在0.6~0.85之間,結(jié)合不同的下墊面類型,不透水路面徑流系數(shù)、綠地徑流系數(shù)、徑流系數(shù)分別選取0.9、0.1、0.9。 表5所示是實際綜合徑流系數(shù)加權(quán)計算結(jié)果。

    表5 實際綜合徑流系數(shù)計算表Tab.5 Comprehensive runoff coefficient calculation table

    SWMM模型模擬采用兩個不同降雨重現(xiàn)期1年一遇和2年一遇,運用SWMM模型進行模擬確定徑流系數(shù),計算結(jié)果如表6所示。

    表6 模擬徑流系數(shù)與計算綜合徑流系數(shù)對比表Tab.6 Simulate runoff coefficient and runoff coefficient calculation comprehensive contrast table

    如表6所示,在1年一遇(P=1)條件下模擬徑流系數(shù)為0.791與實際綜合徑流系數(shù)0.793差異不大。為了判別在2年一遇(P=2)條件下的差異性,引入模擬綜合徑流系數(shù)與實際綜合徑流系數(shù)的變異系數(shù)CV(ψ)。如果變異系數(shù)小于5%,則證明變異性小。

    (1)

    經(jīng)公式(1)計算,變異系數(shù)為3.3%,由于徑流系數(shù)是一個變量,會隨著降雨重現(xiàn)期的增大而增大。又因為此次模擬的降雨徑流系數(shù)的變異系數(shù)小于5%,基本判斷出在2年一遇的強度降雨下模擬得到的降雨徑流系數(shù)變異性小。所以SWMM模型可以準確的模擬研究區(qū)域的實際情況,可以用于下述的情景模擬。

    1.4 模擬情景設計

    1.4.1 設計雨型和雨量

    城市暴雨雨型設計的方法主要有芝加哥法、Pilgrim&Cordery 法、Huff法[10]等。國內(nèi)的水文手冊和室外排水設計規(guī)范的設計暴雨公式一般都采用芝加哥法[11]。根據(jù)《西安市水文手冊》和西安市多年的暴雨資料確定設計暴雨公式如下:

    <20年)

    (2)

    (3)

    式中:i為暴雨強度,mm/min;N為設計重現(xiàn)期,a;t為設計暴雨歷時,min。

    不同重現(xiàn)期設計降雨過程線如圖2所示,取降雨歷時取120 min,雨峰系數(shù)r=0.4,設計4種不同降雨重現(xiàn)期為1年一遇、2年一遇、5年一遇、10年一遇。

    圖2 不同降雨強度下降雨過程線Fig.2 Design different frequency rainfall process line

    1.4.2 設計不同雨峰系數(shù)下降雨過程線

    通過公式(4)、(5)可知,雨峰系數(shù)r(峰值時間與降雨歷時的比值)是芝加哥雨型的唯一參數(shù)。所以通過公式可得到隨著降雨歷時變化的設計降雨過程線。

    (4)

    (5)

    式中:i前為峰前時間的雨強,mm/min;i后為峰后時間的雨強,mm/min;tb為峰前時間,min;ta為峰后時間,min;r為峰值系數(shù);a,b,c為地方參數(shù)。

    如圖3所示是設計重現(xiàn)期P=1年時,雨峰系數(shù)r=0.3、0.4、0.5、0.6、0.7的設計降雨過程線。

    圖3 P=1設計不同雨峰系數(shù)下降雨過程線Fig.3 Design frequency P=1 different rainfall peak rainfall process line

    1.4.3 不同LID措施設置

    情景設計選取在海綿城市建設中應用較為廣泛的LID設施:透水鋪裝路面和綠色屋頂[12,13]。由于考慮到研究區(qū)域為不透水面積較大的老城區(qū),透水路面改造率有限并且有一定的改造難度,所以結(jié)合實際的下墊面情況,在設計不同LID措施情景中,保持透水路面面積不變,設計為15%透水鋪裝路面,15%是指透水路面面積占總道路面積為15%。建筑屋頂由于易于改造,根據(jù)不同屋頂?shù)母脑炻试O計以下4種情景。其中25%、50%、75%綠色屋頂是指綠色屋頂占總建筑屋頂面積為25%、50%、75%。加裝不同類型LID措施的情景設計如表7所示。

    表7 不同情景設計表Tab.7 Design tables for different scenarios

    2 結(jié)果與分析

    2.1 現(xiàn)狀條件下降雨徑流過程模擬

    2.1.1 不同設計雨峰系數(shù)條件下降雨徑流過程

    圖4和表8分別是不同雨峰系數(shù)下排放口處流量過程線和不同雨峰系數(shù)下水文特征值表。

    圖4 不同雨峰系數(shù)條件下流量過程線Fig.4 Flow process line under different rain-peak coefficients

    雨峰系數(shù)洪水總量/萬m3洪峰值/(m3·s-1)峰值時間/min徑流系數(shù)r=0.33.592 36.429510.792r=0.43.593 26.951610.791r=0.53.598 37.148710.790r=0.63.589 97.316820.791

    從圖4和表8可以看出:隨著雨峰系數(shù)的增大,對于排放口產(chǎn)生的洪量和下墊面的徑流系數(shù)基本沒有發(fā)生變化。是因為在相同的下墊面類型和固定的設計降雨量,洪量和徑流系數(shù)是不發(fā)生變化的。但當雨峰系數(shù)每增加0.1,洪峰出現(xiàn)的時間推遲10 min左右。并且洪峰流量也隨之增大,當r=0.3時,Q=6.429 m3/s;r=0.6時,Q=7.316 m3/s,洪峰流量增大13.8%。根據(jù)排水管網(wǎng)的實際情況和模型運行的過程分析可知,在降雨初期時排水管道存在大量空間,能有效排走雨峰值出現(xiàn)越早的雨量,所以此時的洪峰流量值越小。但當雨峰值出現(xiàn)較晚,雨水管網(wǎng)的空間較小,導致洪峰流量增大,使得管道負荷增大。所以也可以看出在雨水管網(wǎng)設計時,城市的雨峰系數(shù)選取是不可忽視的。

    2.1.2 不同設計降雨強度下降雨徑流過程

    圖5以雨峰系數(shù)r=0.4設計P=1、P=2、P=5、P=10 四個不同重現(xiàn)期排放口流量過程線。表9是不同設計降雨強度下降雨徑流過程中水文特征值表。

    圖5 不同降雨強度流量過程線Fig.5 Different rainfall intensity flow process line

    重現(xiàn)期/a洪水總量/萬m3峰值時間/min徑流系數(shù)P=13.654 8600.774P=25.878 0540.803P=57.706 1500.834P=109.867 5480.861

    從圖5和表9可以看出:在降雨重現(xiàn)期P=2年時的洪峰出現(xiàn)時間比P=1年時提前6 min;P=5年時的洪峰出現(xiàn)時間比P=2年時提前4 min;P=10年時的洪峰出現(xiàn)時間比P=5年時提前2 min。得出隨著降雨強度的增大,洪峰出現(xiàn)時間會提前,但洪峰出現(xiàn)時間的速度會變緩。根據(jù)對市政管網(wǎng)建設分析可知,雨水管網(wǎng)的設計規(guī)范采用的是一年一遇的設計標準,當降雨強度增大時會導致一定程度的管網(wǎng)滿流溢流。當降雨重現(xiàn)期P=1、2、5、10年時流域出口產(chǎn)生的洪水總量分別為3.654 8、5.878 0、7.706 1、9.867 5 萬m3,徑流系數(shù)在分別為0.774、0.803、0.834、0.861。是因為降雨強度的增大,排放口管道滿流的時間會提前,持續(xù)時間變長,不能持續(xù)排走雨水,導致洪量值增大。管道滿流時間增長,城市區(qū)域內(nèi)澇積水時間越長,模擬的規(guī)律與城市區(qū)域?qū)嶋H內(nèi)澇中特點基本一致。

    2.2 不同情景下降雨徑流過程模擬

    2.2.1 不同情景下降雨徑流過程

    根據(jù)上述模擬情景設計,運用SWMM模型模擬現(xiàn)狀條件下和不同情境下研究區(qū)域的降雨徑流過程。圖6為設計降雨重現(xiàn)期P=1的情況下研究區(qū)域排放口的流量過程線。表10是不同情景下排放口所對應的水文特征值。

    從圖6和表10可以看出:加裝LID措施,明顯有減少洪水流量,減小洪峰值,推遲洪峰出現(xiàn)時間的效果。加設不同LID措施較現(xiàn)狀條件下洪量分別削減5.4%、17.7%、26.8%、36.6%,洪峰分別削減2.0%、12.4%、17.7%、28.5%。表明隨著LID布設面積比例的增加,排放口洪峰、洪量大幅下降,洪水歷時縮短,洪峰向后延遲。從相關文獻中,模擬的結(jié)果可以得到驗證,王書敏等[14]對綠色屋頂對城市區(qū)域降雨徑流量模擬,得出相對于未加裝綠色屋頂,綠色屋頂?shù)木G化率達到100%,50%時研究區(qū)域峰值徑流量分別降低了5.3%和1.8%。表明綠色屋頂可以儲蓄洪水、削減洪峰,有利于緩解城市內(nèi)澇。

    圖6 不同情景下P=1的排放口流量過程線Fig.6 The discharge outlet flow process line in different scenarios with the recurrence period P=1

    情景設計洪量/萬m3洪量削減率/%洪峰/(m3·s-1)洪峰削減率/%洪峰時間/min洪峰推遲時間/min現(xiàn)狀3.654 807.010580情景13.458 15.46.872.0602情景23.006 917.76.1412.4624情景32.675 426.85.7717.7646情景42.317 636.65.0128.5668

    2.2.2 加設LID措施在不同設計降雨重現(xiàn)期下降雨徑流過程

    為了研究在不同設計降雨重現(xiàn)期,現(xiàn)狀條件和加裝LID措施的降雨徑流規(guī)律,選取LID布設情景4在P=1、P=2、P=5、P=10四個降雨重現(xiàn)期下的流量過程線如圖7所示,不同降雨強度下洪量削減對比表如表11所示。

    圖7 不同設計降雨強度下排放口處流量過程線Fig.7 Discharge process line at discharge outlet under different design rainfall intensity

    從圖7和表11可以看出:隨著設計降雨強度的增大,洪水流量和洪峰值逐漸增大,洪峰出現(xiàn)的時間提前,徑流過程歷時增加。洪量值與現(xiàn)狀條件相比,在重現(xiàn)期 1、2、5、10年時洪量削減量分別為1.337 2、1.862 9、2.352 4、2.211 1 萬m3,洪量衰減率分別為 36.6%,33.7%,29.2%,22.7%,可以看出洪水總量的削減量逐漸增加,但是s洪量的削減幅度減小。表明隨著降雨強度的增加,雨洪削減的效果減小。從相關文獻中,上述的研究規(guī)律可以從印證:翟丹丹等[15]北方典型城市的設計降雨對綠色屋頂裝置進行模擬降雨試驗,發(fā)現(xiàn)在設計降雨重現(xiàn)期在3年一遇時,加裝綠色屋頂措施對降雨徑流的削減效果稍高于重現(xiàn)期為5年的。兩者均在29.0%~33.1% 之間。說明綠色屋頂蓄水容量是有限,隨著雨量的增大,對洪峰洪量削減效果越不明顯。

    表11 不同降雨強度下洪量削減對比表Tab.11 Comparison table of flood reduction under different rainfall intensity

    同時,從圖7排放口流量過程線可以看出,加裝LID措施導致洪峰出現(xiàn)的時間有提前的趨勢。所以進行進一步的探究,以1 min為模擬時間步長,模擬不同重現(xiàn)期下現(xiàn)狀條件下和情景4的降雨徑流過程。分析在不同降雨強度的不同情景下,洪峰出現(xiàn)時間和推遲時間,如表12所示。

    表12 不同降雨強度下洪峰推遲時間對比表Tab.12 Comparison table of flood peak delay time under different rainfall intensities

    從表12可以得到以下結(jié)論:在相同降雨強度下,加裝LID措施明顯能推遲洪峰出現(xiàn)時間,有滯洪錯峰的效果。在重現(xiàn)期 1、2、5、10 時洪峰推遲時間分別為7、6、4、3 min,隨著降雨強度增大,滯洪錯峰的效果明顯減弱。

    3 結(jié) 語

    近年來,隨著城市化的發(fā)展,城市不透水面積激增,內(nèi)澇頻繁發(fā)生。為了解決面臨的難題,提出了海綿城市理念。針對改造困難的老城區(qū),運用低影響開發(fā)理念,改善老城區(qū)排水系統(tǒng)現(xiàn)狀,從而有效控制降雨徑流過程。本文基于城市雨洪模擬技術(shù)方法,運用SWMM模型對獨立排水區(qū)域的現(xiàn)狀條件下和加裝LID措施的雨洪過程進行模擬和分析,得到以下結(jié)論:

    (1)雨峰位置的不同,排水管網(wǎng)的洪峰值變化較大。所以在城市管道設計中,準確的得到雨峰系數(shù),對雨水管道設計標準的選取有較大的參考價值。

    (2)加裝LID措施,能有效改善內(nèi)澇現(xiàn)象。隨著布設LID措施的面積比例增加,排放口處洪量、洪峰值大幅度減小,洪峰出現(xiàn)的時間延后,洪水歷時縮短。

    (3)加裝LID措施,滯洪錯峰效果明顯,但是隨著降雨強度的增大,LID措施對削減雨洪效果明顯減弱。

    猜你喜歡
    排放口徑流系數(shù)過程線
    水下航行器熱射流多孔排放口結(jié)構(gòu)仿真分析
    遮流罩排放口對水下航行器熱射流溫度場影響的數(shù)值仿真
    降雨特征對半透水道路徑流系數(shù)的影響
    水資源保護(2022年3期)2022-05-31 03:06:36
    煤制油項目有組織排放口VOCs排放源強分析
    化工管理(2021年19期)2021-08-02 01:25:38
    基于Excel繪制改正系數(shù)過程線淺析
    智能城市(2018年7期)2018-07-10 08:30:00
    基于青山水庫洪水調(diào)節(jié)論述給排水系統(tǒng)設計
    科學與財富(2018年7期)2018-05-21 08:46:30
    無資料山丘小流域徑流系數(shù)的計算研究
    基于青山水庫論述調(diào)洪編程計算過程
    第一類污染物廢水總排放口“超標”應如何處理
    Estimating runoff coefficient for quantity assessment of roof rainwater harvesting system
    国产乱来视频区| 秋霞伦理黄片| 成年人午夜在线观看视频| av国产久精品久网站免费入址| 我要看日韩黄色一级片| 久久人人爽人人爽人人片va| av国产免费在线观看| 深夜a级毛片| 国产精品国产三级专区第一集| 色网站视频免费| 看免费成人av毛片| 91狼人影院| 午夜免费观看性视频| 亚洲人成网站高清观看| 精品亚洲成a人片在线观看 | 精品视频人人做人人爽| 丝瓜视频免费看黄片| 嘟嘟电影网在线观看| 少妇人妻精品综合一区二区| 免费看日本二区| 国产v大片淫在线免费观看| 干丝袜人妻中文字幕| 下体分泌物呈黄色| 久久午夜福利片| 在线观看免费高清a一片| 天堂8中文在线网| 欧美日韩视频高清一区二区三区二| 多毛熟女@视频| 大片电影免费在线观看免费| 亚洲伊人久久精品综合| 色婷婷av一区二区三区视频| 一区二区av电影网| 久久99热这里只有精品18| 3wmmmm亚洲av在线观看| 一区在线观看完整版| av国产免费在线观看| 亚洲天堂av无毛| 国产精品成人在线| 黑人高潮一二区| 少妇人妻久久综合中文| 国产精品一及| 免费不卡的大黄色大毛片视频在线观看| 国产伦精品一区二区三区视频9| 精华霜和精华液先用哪个| 亚洲熟女精品中文字幕| 欧美日韩视频精品一区| 亚洲欧美日韩无卡精品| 亚洲天堂av无毛| 精品久久久久久电影网| 国产伦精品一区二区三区视频9| 亚洲欧美成人综合另类久久久| av国产免费在线观看| 日韩不卡一区二区三区视频在线| 一区二区av电影网| 日韩制服骚丝袜av| 日韩伦理黄色片| 日韩中字成人| 亚洲欧美成人综合另类久久久| 久久久成人免费电影| 爱豆传媒免费全集在线观看| av免费在线看不卡| 亚洲怡红院男人天堂| 久久久久久伊人网av| 一本—道久久a久久精品蜜桃钙片| 日韩三级伦理在线观看| 在线天堂最新版资源| 亚洲欧美中文字幕日韩二区| 99精国产麻豆久久婷婷| a级一级毛片免费在线观看| 插阴视频在线观看视频| 国产伦精品一区二区三区视频9| 欧美精品一区二区免费开放| 熟女av电影| 人妻 亚洲 视频| 亚洲无线观看免费| 国产欧美亚洲国产| .国产精品久久| 成年美女黄网站色视频大全免费 | 97热精品久久久久久| 国模一区二区三区四区视频| 能在线免费看毛片的网站| 能在线免费看毛片的网站| 99久国产av精品国产电影| 国产美女午夜福利| 美女高潮的动态| 特大巨黑吊av在线直播| 伊人久久国产一区二区| 少妇精品久久久久久久| 日韩电影二区| 伊人久久精品亚洲午夜| 伦精品一区二区三区| 成人高潮视频无遮挡免费网站| 亚洲成色77777| 免费观看在线日韩| 国产爱豆传媒在线观看| 国产爱豆传媒在线观看| 日本av免费视频播放| 又大又黄又爽视频免费| 日韩制服骚丝袜av| 大片电影免费在线观看免费| 中国国产av一级| 六月丁香七月| av国产久精品久网站免费入址| 国产淫片久久久久久久久| 成年免费大片在线观看| 老熟女久久久| 中文字幕久久专区| 在线精品无人区一区二区三 | 免费大片18禁| 男男h啪啪无遮挡| 中文乱码字字幕精品一区二区三区| 熟女人妻精品中文字幕| 免费黄色在线免费观看| 香蕉精品网在线| 亚洲人成网站高清观看| av视频免费观看在线观看| 国产精品不卡视频一区二区| 人妻制服诱惑在线中文字幕| 国产av码专区亚洲av| 国产永久视频网站| 日产精品乱码卡一卡2卡三| 视频区图区小说| 自拍欧美九色日韩亚洲蝌蚪91 | 伦理电影免费视频| 久久97久久精品| av不卡在线播放| 免费看av在线观看网站| 国产在线一区二区三区精| 91狼人影院| 色网站视频免费| 久久久午夜欧美精品| 久久国产乱子免费精品| 99热网站在线观看| 久久影院123| 国产精品秋霞免费鲁丝片| 国产欧美另类精品又又久久亚洲欧美| 久久ye,这里只有精品| 久久国产亚洲av麻豆专区| 2021少妇久久久久久久久久久| 国产 精品1| 18禁动态无遮挡网站| 国产日韩欧美在线精品| 国产又色又爽无遮挡免| 欧美激情极品国产一区二区三区 | 伦理电影大哥的女人| 大陆偷拍与自拍| 欧美精品一区二区大全| 亚洲精品乱码久久久久久按摩| 又大又黄又爽视频免费| 国产成人精品婷婷| 色视频在线一区二区三区| 插阴视频在线观看视频| 午夜激情福利司机影院| 中文字幕免费在线视频6| 精品亚洲成a人片在线观看 | 中文精品一卡2卡3卡4更新| 久久久久久久精品精品| 一本一本综合久久| 香蕉精品网在线| 精品视频人人做人人爽| 看十八女毛片水多多多| 午夜福利网站1000一区二区三区| 欧美亚洲 丝袜 人妻 在线| 亚洲怡红院男人天堂| 在线观看三级黄色| 少妇裸体淫交视频免费看高清| 丰满少妇做爰视频| 在线免费十八禁| 国产精品伦人一区二区| 一级爰片在线观看| 免费看光身美女| 制服丝袜香蕉在线| 街头女战士在线观看网站| 国产精品av视频在线免费观看| 一级a做视频免费观看| 日日摸夜夜添夜夜爱| 亚洲av免费高清在线观看| 香蕉精品网在线| av专区在线播放| 亚洲av二区三区四区| 成人毛片a级毛片在线播放| 哪个播放器可以免费观看大片| 色哟哟·www| 超碰av人人做人人爽久久| 中文乱码字字幕精品一区二区三区| 一级av片app| 一级毛片aaaaaa免费看小| 人妻制服诱惑在线中文字幕| 中文字幕av成人在线电影| 久久久成人免费电影| 激情 狠狠 欧美| 在现免费观看毛片| 欧美成人精品欧美一级黄| 国产成人a区在线观看| 欧美激情极品国产一区二区三区 | 成人午夜精彩视频在线观看| 国产伦在线观看视频一区| 亚洲成人中文字幕在线播放| 午夜免费鲁丝| 中文字幕制服av| 大香蕉久久网| 菩萨蛮人人尽说江南好唐韦庄| 黄色怎么调成土黄色| 国产69精品久久久久777片| 丝袜喷水一区| 一个人看的www免费观看视频| 国产精品成人在线| 久久人人爽人人爽人人片va| av在线播放精品| 少妇人妻精品综合一区二区| 亚洲国产精品国产精品| 人人妻人人爽人人添夜夜欢视频 | 成人一区二区视频在线观看| 欧美最新免费一区二区三区| 只有这里有精品99| 舔av片在线| 色婷婷久久久亚洲欧美| 久久久久久久久久久丰满| 成年女人在线观看亚洲视频| .国产精品久久| 免费观看a级毛片全部| 久久影院123| 五月天丁香电影| 伊人久久国产一区二区| 久久鲁丝午夜福利片| 国产精品国产三级专区第一集| 国产亚洲91精品色在线| 成人无遮挡网站| 麻豆乱淫一区二区| 国产色爽女视频免费观看| 国产白丝娇喘喷水9色精品| 亚洲精品久久久久久婷婷小说| www.色视频.com| 好男人视频免费观看在线| 国产在线一区二区三区精| 日日摸夜夜添夜夜添av毛片| 精品熟女少妇av免费看| 日韩精品有码人妻一区| 18禁在线播放成人免费| 成年人午夜在线观看视频| 精品久久久噜噜| 国产精品久久久久久久电影| 99热国产这里只有精品6| 色视频在线一区二区三区| 在线看a的网站| 在线观看免费日韩欧美大片 | 青春草视频在线免费观看| 亚洲国产日韩一区二区| 国产亚洲5aaaaa淫片| 国产国拍精品亚洲av在线观看| 黄色怎么调成土黄色| 三级国产精品片| 久久这里有精品视频免费| av在线蜜桃| 亚洲国产精品专区欧美| 欧美日韩视频高清一区二区三区二| 啦啦啦啦在线视频资源| 一级毛片aaaaaa免费看小| 午夜福利在线观看免费完整高清在| 爱豆传媒免费全集在线观看| 一级毛片我不卡| 中文字幕人妻熟人妻熟丝袜美| 在线精品无人区一区二区三 | 成人无遮挡网站| 国产91av在线免费观看| 五月开心婷婷网| 亚洲欧美日韩无卡精品| 伦理电影免费视频| 最近2019中文字幕mv第一页| 男女无遮挡免费网站观看| 中文天堂在线官网| 日韩av免费高清视频| 亚洲国产精品成人久久小说| 99久久人妻综合| 老司机影院毛片| 日韩伦理黄色片| 超碰av人人做人人爽久久| 韩国av在线不卡| 亚洲三级黄色毛片| av.在线天堂| 高清日韩中文字幕在线| 91精品国产国语对白视频| 亚洲av不卡在线观看| 亚洲在久久综合| 久久国产乱子免费精品| 天堂8中文在线网| 亚洲av综合色区一区| 欧美人与善性xxx| 亚洲美女黄色视频免费看| 人人妻人人添人人爽欧美一区卜 | 成人午夜精彩视频在线观看| 久久久午夜欧美精品| 免费av中文字幕在线| a 毛片基地| av视频免费观看在线观看| 2022亚洲国产成人精品| 精品一区在线观看国产| 国产精品女同一区二区软件| 久久久久精品性色| 成人午夜精彩视频在线观看| 亚洲精品视频女| 国产真实伦视频高清在线观看| 亚洲欧洲日产国产| 国产精品久久久久久av不卡| 久久久久精品性色| 日韩在线高清观看一区二区三区| 蜜臀久久99精品久久宅男| 成人影院久久| 欧美丝袜亚洲另类| 国产在线男女| 男人舔奶头视频| 我要看日韩黄色一级片| 久久精品国产亚洲网站| 岛国毛片在线播放| 色吧在线观看| 女人十人毛片免费观看3o分钟| 成年av动漫网址| 伦理电影大哥的女人| 少妇人妻 视频| av国产久精品久网站免费入址| 国产深夜福利视频在线观看| 亚洲成人av在线免费| 99热网站在线观看| 亚洲av免费高清在线观看| 久久精品熟女亚洲av麻豆精品| 交换朋友夫妻互换小说| 九色成人免费人妻av| 国产精品99久久99久久久不卡 | 国产成人freesex在线| 街头女战士在线观看网站| 国产综合精华液| 国产精品久久久久久久久免| 九草在线视频观看| 亚洲国产精品专区欧美| 婷婷色综合www| 青春草视频在线免费观看| 精品人妻偷拍中文字幕| 成人特级av手机在线观看| 亚洲精品成人av观看孕妇| 色5月婷婷丁香| 黄片wwwwww| videossex国产| 丝袜喷水一区| 在线观看av片永久免费下载| 我要看日韩黄色一级片| 婷婷色av中文字幕| 亚洲最大成人中文| 精品久久久久久电影网| 色婷婷久久久亚洲欧美| 国产精品人妻久久久久久| 一区在线观看完整版| 嫩草影院入口| 夜夜爽夜夜爽视频| 亚洲激情五月婷婷啪啪| 在线观看美女被高潮喷水网站| 中文字幕制服av| 天堂中文最新版在线下载| 男女国产视频网站| 免费看av在线观看网站| 五月开心婷婷网| 国产精品人妻久久久影院| 97精品久久久久久久久久精品| 99国产精品免费福利视频| 国产精品av视频在线免费观看| 成人无遮挡网站| 亚洲av成人精品一区久久| 男女边摸边吃奶| 久久人人爽av亚洲精品天堂 | 欧美日韩一区二区视频在线观看视频在线| 国产精品免费大片| 中文乱码字字幕精品一区二区三区| 日韩在线高清观看一区二区三区| 少妇丰满av| 一本久久精品| 免费看不卡的av| 国产亚洲5aaaaa淫片| 99九九线精品视频在线观看视频| 中文字幕制服av| 1000部很黄的大片| 超碰av人人做人人爽久久| av线在线观看网站| 亚洲精品视频女| 啦啦啦中文免费视频观看日本| 国产日韩欧美亚洲二区| 午夜福利在线观看免费完整高清在| 久久ye,这里只有精品| 久久久久性生活片| 国产探花极品一区二区| 亚洲伊人久久精品综合| 亚洲欧美一区二区三区国产| 国产在线免费精品| 大片电影免费在线观看免费| 超碰97精品在线观看| 99视频精品全部免费 在线| 久久ye,这里只有精品| 18禁动态无遮挡网站| 欧美xxⅹ黑人| 人人妻人人添人人爽欧美一区卜 | 黄色配什么色好看| 精品国产乱码久久久久久小说| 插阴视频在线观看视频| 免费看不卡的av| 欧美bdsm另类| 久久久久精品性色| 人妻制服诱惑在线中文字幕| 黑人高潮一二区| 国产v大片淫在线免费观看| 日本-黄色视频高清免费观看| 久久人人爽av亚洲精品天堂 | 男女下面进入的视频免费午夜| 六月丁香七月| 久久精品熟女亚洲av麻豆精品| 中文欧美无线码| av不卡在线播放| 国产成人一区二区在线| 色婷婷av一区二区三区视频| 亚洲欧美一区二区三区国产| 性高湖久久久久久久久免费观看| 中文精品一卡2卡3卡4更新| 免费看av在线观看网站| 你懂的网址亚洲精品在线观看| 欧美一区二区亚洲| 极品少妇高潮喷水抽搐| 永久免费av网站大全| 亚洲国产精品专区欧美| 男人和女人高潮做爰伦理| 日韩一区二区视频免费看| 精品国产一区二区三区久久久樱花 | 久久婷婷青草| 精品一区二区免费观看| 乱系列少妇在线播放| 伦理电影大哥的女人| 久久久国产一区二区| 欧美日韩亚洲高清精品| 在线观看av片永久免费下载| 一二三四中文在线观看免费高清| 精品亚洲成a人片在线观看 | 毛片一级片免费看久久久久| 日韩一区二区视频免费看| 免费播放大片免费观看视频在线观看| 老司机影院成人| 九草在线视频观看| 免费久久久久久久精品成人欧美视频 | 中文天堂在线官网| 亚洲熟女精品中文字幕| 免费黄色在线免费观看| 永久免费av网站大全| 春色校园在线视频观看| 亚洲国产高清在线一区二区三| 好男人视频免费观看在线| 成人高潮视频无遮挡免费网站| 激情 狠狠 欧美| 久久久精品免费免费高清| 久久人妻熟女aⅴ| 高清午夜精品一区二区三区| 日本vs欧美在线观看视频 | 免费观看无遮挡的男女| 80岁老熟妇乱子伦牲交| 成人国产麻豆网| 色视频www国产| 日本免费在线观看一区| 亚洲欧美日韩卡通动漫| 日韩欧美一区视频在线观看 | 久久精品国产自在天天线| 精品一区二区三区视频在线| 日韩大片免费观看网站| 91aial.com中文字幕在线观看| 日韩亚洲欧美综合| 网址你懂的国产日韩在线| 欧美性感艳星| 观看免费一级毛片| 国产一区二区在线观看日韩| 亚洲欧美一区二区三区国产| av福利片在线观看| 成人18禁高潮啪啪吃奶动态图 | 日本色播在线视频| 少妇人妻精品综合一区二区| 99视频精品全部免费 在线| 黄色怎么调成土黄色| 精品久久久噜噜| 日本wwww免费看| 性色av一级| 男女啪啪激烈高潮av片| 美女cb高潮喷水在线观看| 亚洲经典国产精华液单| 一级毛片久久久久久久久女| h视频一区二区三区| 一本色道久久久久久精品综合| 热re99久久精品国产66热6| 边亲边吃奶的免费视频| 一区二区三区免费毛片| 国产高清不卡午夜福利| 成人影院久久| 99热这里只有是精品50| 亚洲国产av新网站| 热re99久久精品国产66热6| 亚洲av日韩在线播放| 午夜老司机福利剧场| 免费大片黄手机在线观看| 我要看日韩黄色一级片| 免费人成在线观看视频色| 少妇人妻久久综合中文| 久久久久久久国产电影| 国产午夜精品久久久久久一区二区三区| 欧美xxxx性猛交bbbb| 国产免费一级a男人的天堂| 一区二区三区乱码不卡18| 女性生殖器流出的白浆| 少妇裸体淫交视频免费看高清| 人人妻人人看人人澡| 亚洲精品色激情综合| 丝瓜视频免费看黄片| 国产乱人视频| av女优亚洲男人天堂| 亚洲欧美一区二区三区黑人 | 亚洲激情五月婷婷啪啪| 国产精品熟女久久久久浪| 日韩,欧美,国产一区二区三区| 亚洲国产精品一区三区| h视频一区二区三区| 日韩av免费高清视频| 成人综合一区亚洲| 丰满人妻一区二区三区视频av| 91精品国产九色| 在线观看免费日韩欧美大片 | 久久精品国产自在天天线| 日本黄色片子视频| 亚洲精品,欧美精品| 韩国av在线不卡| 久久久午夜欧美精品| 亚洲精品第二区| 久久久久久九九精品二区国产| 女性被躁到高潮视频| 亚洲国产高清在线一区二区三| 男男h啪啪无遮挡| 成人免费观看视频高清| 国产黄片视频在线免费观看| 亚洲在久久综合| 校园人妻丝袜中文字幕| 内射极品少妇av片p| 国产精品女同一区二区软件| kizo精华| 国产亚洲午夜精品一区二区久久| 国产色婷婷99| 最近最新中文字幕免费大全7| 国产乱人偷精品视频| 亚洲欧洲国产日韩| 欧美极品一区二区三区四区| 大话2 男鬼变身卡| 国产精品精品国产色婷婷| 欧美日韩亚洲高清精品| 久久鲁丝午夜福利片| 国产永久视频网站| 午夜激情久久久久久久| 下体分泌物呈黄色| 十八禁网站网址无遮挡 | 高清毛片免费看| 精品国产一区二区三区久久久樱花 | 七月丁香在线播放| 亚洲第一av免费看| 国产国拍精品亚洲av在线观看| av专区在线播放| 97在线视频观看| 国产精品人妻久久久久久| 国产精品精品国产色婷婷| 中文资源天堂在线| 亚洲国产精品一区三区| 99国产精品免费福利视频| 超碰av人人做人人爽久久| 久久精品国产亚洲av涩爱| 高清不卡的av网站| 国产精品一区二区在线观看99| 大片免费播放器 马上看| 精品国产一区二区三区久久久樱花 | 国产精品久久久久久av不卡| 亚洲经典国产精华液单| 插阴视频在线观看视频| 免费观看在线日韩| 亚洲四区av| 人妻夜夜爽99麻豆av| 最近的中文字幕免费完整| 中文欧美无线码| 亚洲精品456在线播放app| 国产精品福利在线免费观看| 免费在线观看成人毛片| 亚洲精品,欧美精品| 亚洲精品456在线播放app| 日韩伦理黄色片| 极品少妇高潮喷水抽搐| 国产亚洲91精品色在线| 日韩欧美一区视频在线观看 | 青青草视频在线视频观看| www.av在线官网国产| 国精品久久久久久国模美| 亚洲av二区三区四区| 能在线免费看毛片的网站| 男女免费视频国产| 黄色视频在线播放观看不卡| 亚洲激情五月婷婷啪啪| 亚洲人与动物交配视频| 免费高清在线观看视频在线观看| 亚洲国产成人一精品久久久| 又大又黄又爽视频免费| 又黄又爽又刺激的免费视频.| 九色成人免费人妻av| 国产精品不卡视频一区二区| 久久97久久精品| 国产有黄有色有爽视频| 国产人妻一区二区三区在| 国产黄频视频在线观看| 亚洲图色成人| 午夜激情久久久久久久| 少妇的逼水好多| 久久久久久久精品精品| 亚洲av二区三区四区|