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

    萬安水庫后汛期實(shí)時(shí)防洪預(yù)報(bào)優(yōu)化調(diào)度策略研究

    2022-11-28 06:38:22程穎新何中政陳濟(jì)天劉孟樺熊芳金
    中國(guó)農(nóng)村水利水電 2022年11期
    關(guān)鍵詞:萬安吉安防洪

    程穎新,何中政,2,陳濟(jì)天,劉孟樺,熊芳金,,尹 恒

    (1.南昌大學(xué)工程建設(shè)學(xué)院,江西南昌 330031;2.南昌大學(xué)鄱陽湖環(huán)境與資源利用教育部重點(diǎn)實(shí)驗(yàn)室,江西南昌 330031;3.江西省水文監(jiān)測(cè)中心,江西南昌 330002;4.江西水利職業(yè)學(xué)院,江西南昌 330013)

    0 引言

    水庫(群)防洪調(diào)度是流域防洪重要的非工程措施[1],也是目前防洪非工程措施研究的熱點(diǎn),國(guó)內(nèi)外諸多學(xué)者已做了大量研究[2]。針對(duì)確定性來水條件下水庫(群)防洪優(yōu)化調(diào)度問題的研究主要集中在優(yōu)化方法,如:線性規(guī)劃、動(dòng)態(tài)規(guī)劃類優(yōu)化方法[3]和智能優(yōu)化算法[4]等。優(yōu)化方法研究對(duì)象從單個(gè)水庫、梯級(jí)水庫群、混聯(lián)水庫群,研究目標(biāo)也從單個(gè)防洪目標(biāo)拓展到多重防洪優(yōu)化目標(biāo)或汛期洪水資源綜合利用目標(biāo)[5]。

    除確定性優(yōu)化方法研究外,來水不確定性條件下水庫(群)防洪調(diào)度隨機(jī)優(yōu)化方法的研究也是熱點(diǎn)之一。而隨著地理信息系統(tǒng)技術(shù)、測(cè)雨雷達(dá)和衛(wèi)星云圖技術(shù)應(yīng)用,各時(shí)效預(yù)報(bào)“定點(diǎn)、定時(shí)、大量”的準(zhǔn)確率都有很大提高[6]。鐘華昱等[7]耦合集合預(yù)報(bào)信息開展了“參數(shù)-模擬-優(yōu)化”框架的水庫多目標(biāo)優(yōu)化調(diào)度規(guī)則研究;宋君輝等[8]根據(jù)凈雨過程與入庫流量過程擬定水庫預(yù)報(bào)調(diào)度規(guī)則,但上述研究主要聚焦于結(jié)合徑流特性和預(yù)報(bào)信息來制定調(diào)度規(guī)則。而在發(fā)電調(diào)度研究領(lǐng)域,長(zhǎng)中短實(shí)時(shí)嵌套預(yù)報(bào)優(yōu)化調(diào)度[9]和短期預(yù)報(bào)優(yōu)化調(diào)度研究[10]已取得了一定成效,而在防洪調(diào)度研究領(lǐng)域涉及實(shí)時(shí)防洪預(yù)報(bào)優(yōu)化調(diào)度相關(guān)內(nèi)容較少。因此,以水庫針對(duì)下游控制點(diǎn)防洪優(yōu)化調(diào)度問題為例,提出了一種水庫實(shí)時(shí)防洪預(yù)報(bào)優(yōu)化調(diào)度策略,結(jié)合預(yù)報(bào)信息開展優(yōu)化調(diào)度攔蓄洪水,并以江西省贛江中游大型控制性水庫萬安對(duì)下游吉安站后汛期防洪調(diào)度開展了相關(guān)實(shí)例研究。

    1 研究方法

    1.1 水庫防洪優(yōu)化調(diào)度模型

    針對(duì)下游控制點(diǎn)的水庫防洪調(diào)度問題中,大壩自身的防洪控制安全目標(biāo)可以轉(zhuǎn)化為運(yùn)行水位、水位變幅、流量和流量變幅等約束進(jìn)行表示??紤]最大削峰準(zhǔn)則的水庫防洪優(yōu)化調(diào)度模型可以表示為:

    式中:T為水庫防洪調(diào)度時(shí)段數(shù);分別表示下游防洪控制點(diǎn)t時(shí)段初、末時(shí)刻流量,m3∕s;此外,還需要考慮一些安全調(diào)度約束:

    (1)水庫水量平衡。

    (2)庫容約束。

    (3)水位變幅約束。

    (4)流量約束。

    (5)流量變幅約束。

    (6)下游防洪控制點(diǎn)安全流量約束。

    式中:It表示在t時(shí)段的入庫流量;Qt-1和Qt分別表示t時(shí)段初、末時(shí)刻的出庫流量,m3∕s;ΔT表示調(diào)度時(shí)間步長(zhǎng),s;Vt、和分別表示在t時(shí)刻的水庫庫容及其上下限,m3;Zt表示在t時(shí)刻運(yùn)行水位,而ΔZt表示在t時(shí)段水位變幅限制,m;分別表示在t時(shí)刻下泄流量上下限,m3∕s;ΔQt表示在t時(shí)段的下泄流量變幅限制,m3∕s;表示水庫下泄流量Qt演進(jìn)到下游防洪控制點(diǎn)流量,m3∕s,研究選用馬斯京根法來進(jìn)行河道流量演算;qt表示水庫到下游防洪控制點(diǎn)的區(qū)間來水,m3∕s;Qfc為下游防洪控制點(diǎn)安全泄量,m3∕s。

    1.2 水庫防洪實(shí)時(shí)預(yù)報(bào)優(yōu)化調(diào)度策略

    不同于洪水過程已知的確定性防洪優(yōu)化調(diào)度,水庫實(shí)時(shí)防洪預(yù)報(bào)優(yōu)化調(diào)度往往只能獲取到存在預(yù)報(bào)誤差的有限預(yù)見期內(nèi)的預(yù)報(bào)洪水過程。實(shí)時(shí)防洪調(diào)度除依據(jù)調(diào)度規(guī)則開展外,還可以通過滾動(dòng)地實(shí)施“預(yù)報(bào)-優(yōu)化-調(diào)度”過程,輔助水庫防洪調(diào)度實(shí)時(shí)運(yùn)行。以洪水預(yù)報(bào)預(yù)見期為τ個(gè)時(shí)段為例,防洪實(shí)時(shí)預(yù)報(bào)調(diào)度的優(yōu)化目標(biāo)需要從式(1)的全時(shí)段優(yōu)化,修改為從當(dāng)前時(shí)段t到時(shí)段t+τ進(jìn)行優(yōu)化,模型輸入也僅采用預(yù)見期時(shí)段t到時(shí)段t+τ內(nèi)的預(yù)報(bào)洪水過程。

    由于水庫針對(duì)下游控制點(diǎn)防洪調(diào)度需要考慮河道洪水演進(jìn),使得上述防洪優(yōu)化調(diào)度問題存在時(shí)滯特性。而POA 主要用于解決多階段的動(dòng)態(tài)決策問題,針對(duì)非線性、時(shí)滯問題仍具有全局最優(yōu)收斂的能力[11],研究采用POA求解上述問題。

    在預(yù)見期內(nèi),首先獲得預(yù)報(bào)信息,然后開展水庫防洪優(yōu)化調(diào)度計(jì)算,依據(jù)優(yōu)化調(diào)度方案指導(dǎo)水庫實(shí)際調(diào)度。綜上,研究提出的水庫防洪實(shí)時(shí)預(yù)報(bào)優(yōu)化調(diào)度策略流程如圖1 所示,具體實(shí)施步驟為:①標(biāo)記當(dāng)前面臨時(shí)段為t,獲得到時(shí)段t+τ-1 預(yù)見期內(nèi)的預(yù)報(bào)洪水過程;②考慮約束條件,根據(jù)式(9)建立預(yù)見期內(nèi)的預(yù)報(bào)優(yōu)化調(diào)度模型,采用POA 求解該模型,獲得最優(yōu)決策序列;③在面臨的時(shí)段t,水庫按照優(yōu)化結(jié)果進(jìn)行下泄;④到下一時(shí)段,令t=t+1,重復(fù)步驟(1)~(3)。

    圖1 實(shí)時(shí)防洪調(diào)度預(yù)報(bào)優(yōu)化調(diào)度策略流程圖Fig.1 Flow chart of o forecast optimization operation strategy for real-time flood control operation

    1.3 洪水預(yù)報(bào)誤差隨機(jī)模擬

    洪水預(yù)報(bào)誤差受不同預(yù)報(bào)模型影響,其分布和特征各不相同。部分學(xué)者將其簡(jiǎn)化視為正態(tài)分布,從預(yù)報(bào)精度評(píng)定等級(jí)出發(fā)[12],考慮服從正態(tài)分布的預(yù)報(bào)誤差,入庫洪水過程及其各時(shí)段標(biāo)準(zhǔn)差可表示為:

    2 研究對(duì)象

    2.1 研究區(qū)域概況

    萬安水庫壩址坐落在贛江中游上段,位于萬安縣城以上約2 km,壩址多年平均徑流量299 億m3,約占贛江總流域面積的44.2%,占贛江入湖水量的43.5%。樞紐主要以發(fā)電為主,兼有防洪、航運(yùn)、灌溉、養(yǎng)殖等綜合效益,大壩在設(shè)計(jì)時(shí)安全系數(shù)較高,壩頂高程105.0 m,汛限水位(死水位)為90.0 m,正常蓄水位為100.0 m。但因移民問題及投資等因素的制約,檢驗(yàn)水庫蓄水后對(duì)贛州市淹沒和泥沙淤積的影響,將初期運(yùn)行時(shí)的汛后蓄水位定為96.0 m,汛限水位和死水位為85.0 m,防洪高水位93.6 m,總庫容11.17 億m3(96.0 m),防洪庫容5.7 億m3,調(diào)洪庫容7.98 億m3,水庫調(diào)節(jié)系數(shù)0.026 7,屬于不完全年調(diào)節(jié)水庫,是贛江流域防洪的主要控制性工程之一,萬安水庫區(qū)域概況如圖2所示。

    圖2 研究對(duì)象萬安水庫區(qū)域概況圖Fig.2 Spatial profile of Wan'an Reservoir

    2.2 萬安水庫現(xiàn)行后汛期洪水調(diào)度方式

    根據(jù)2017年江西省防辦批復(fù)的《萬安水庫汛期調(diào)度運(yùn)行方案》,萬安水庫的防洪任務(wù)是在確保樞紐自身防洪安全條件下,不加重上游贛州市的防洪影響,兼顧下游的防洪安全,其中下游萬安縣城安全泄量為8 800 m3∕s。萬安水庫后汛期為6 月21日至9 月30 日,運(yùn)行水位按93.5~96.0 m 動(dòng)態(tài)控制,在后汛期洪水來臨前,將庫水位預(yù)泄至93.5 m 以下,其洪水調(diào)度方式如表1所示。當(dāng)流域發(fā)生大雨,或預(yù)報(bào)24 h 內(nèi)有暴雨及以上天氣時(shí),其后汛期預(yù)泄調(diào)度方式如表2所示。

    表1 萬安水庫后汛期防洪調(diào)度規(guī)則Tab.1 Flood control regulation operation of Wan′an reservoir in late flood period

    表2 萬安水庫后汛期預(yù)泄調(diào)度規(guī)則Tab.2 Pre-discharge regulation of Wan′an reservoir in late flood season

    3 實(shí)例研究

    建立萬安水庫防洪實(shí)時(shí)預(yù)報(bào)優(yōu)化調(diào)度模型,時(shí)段步長(zhǎng)設(shè)定為6 h,水庫下泄演進(jìn)到下游吉安站的時(shí)滯約18 h。實(shí)例研究設(shè)定了6 種洪水預(yù)見期:24、30、36、48、60、72 h,即τ=4、5、6、8、10、12。分析萬安水庫后汛期預(yù)泄調(diào)度規(guī)則模擬結(jié)果,設(shè)定流量變幅約束為6 000 m3∕(s·6 h),水位變幅約束為1.5 m∕(6 h)。

    3.1 不考慮預(yù)報(bào)誤差的實(shí)例分析

    在本次實(shí)例實(shí)驗(yàn)中不考慮預(yù)報(bào)誤差,預(yù)見期內(nèi)的預(yù)報(bào)信息采用實(shí)際來水進(jìn)行模擬。對(duì)比分析中加入了現(xiàn)行調(diào)度規(guī)則模擬和POA 確定性優(yōu)化結(jié)果。輸入洪水選用了1992 年7 月和2019 年7 月兩場(chǎng)50 年一遇典型設(shè)計(jì)洪水。POA 優(yōu)化結(jié)果使用了全過程的洪水信息,此處僅作為對(duì)比分析中削峰率的理論上限值。表3給了不同方案得到的吉安站洪峰流量和削峰率統(tǒng)計(jì)結(jié)果,實(shí)驗(yàn)結(jié)果表明:依據(jù)現(xiàn)行調(diào)度規(guī)則,萬安水庫對(duì)吉安的攔洪削峰作用較弱,削峰率分別為0.15%、0.64%;以預(yù)見期τ=8為例,實(shí)時(shí)預(yù)報(bào)優(yōu)化調(diào)度策略兩場(chǎng)洪水的削峰率分別為8.43%和13.84%,與POA 優(yōu)化結(jié)果存在差距,但明顯優(yōu)于現(xiàn)行調(diào)度規(guī)則。

    表3 設(shè)計(jì)洪水條件下不同方法的調(diào)度結(jié)果統(tǒng)計(jì)Tab.3 Statistical results of flood control operation by different methods under design flood conditions

    此外,當(dāng)τ≥5 時(shí)預(yù)報(bào)實(shí)時(shí)優(yōu)化調(diào)度策略才會(huì)有削峰效果。這是未來4 個(gè)時(shí)段內(nèi)吉安站流量過程只有第4 個(gè)時(shí)段受萬安水庫第1 個(gè)時(shí)段的下泄影響(萬安水庫到下游吉安站時(shí)滯為3 個(gè)時(shí)段),當(dāng)τ<5 且預(yù)報(bào)未來有洪水時(shí),水庫執(zhí)行預(yù)泄只會(huì)增加未來4 個(gè)時(shí)段內(nèi)下游吉安站的流量過程,使得在預(yù)報(bào)實(shí)時(shí)優(yōu)化調(diào)度策略中水庫無法開展預(yù)泄。這表明水庫對(duì)下游控制點(diǎn)進(jìn)行實(shí)時(shí)防洪預(yù)報(bào)優(yōu)化調(diào)度時(shí),洪水預(yù)見期至少要比洪水演進(jìn)時(shí)滯多2個(gè)時(shí)段。

    圖3 和圖4 給出了1992 年7 月型洪水條件下不同調(diào)度方式的調(diào)洪水位過程和下游吉安站流量過程。已知全過程洪水信息的確定性優(yōu)化過程和考慮預(yù)見期內(nèi)洪水信息的實(shí)時(shí)預(yù)報(bào)優(yōu)化調(diào)度的水位過程趨勢(shì)基本一致,僅在預(yù)泄和部分?jǐn)r蓄洪水時(shí)機(jī)選擇上不同。結(jié)合圖4 吉安站流量過程可知,現(xiàn)行調(diào)度規(guī)則僅針對(duì)面臨時(shí)段洪水進(jìn)行調(diào)度,對(duì)吉安站洪峰的攔蓄缺乏針對(duì)性,使得削峰效果較差。

    圖3 萬安水庫水位過程Fig.3 Water level process of Wan'an Reservoir

    圖4 吉安站流量過程Fig.4 Flow process of Ji'an hydrological station

    3.2 考慮預(yù)報(bào)誤差的實(shí)例分析

    進(jìn)一步考慮預(yù)報(bào)誤差,分析預(yù)報(bào)實(shí)時(shí)優(yōu)化調(diào)度策略的應(yīng)用效果。根據(jù)1.3 小節(jié)中描述的預(yù)報(bào)誤差隨機(jī)模擬方法,運(yùn)用蒙特卡洛方法隨機(jī)模擬10 000 次,針對(duì)2019 年7 月型洪水甲、乙、丙3種不同預(yù)報(bào)精度下的實(shí)時(shí)預(yù)報(bào)優(yōu)化調(diào)度結(jié)果如表4所示。

    表4 不同預(yù)報(bào)精度下預(yù)報(bào)優(yōu)化調(diào)度策略調(diào)度結(jié)果Tab.4 Application effect of forecast optimal dispatching under different forecast accuracy

    當(dāng)τ=5 時(shí)不同預(yù)報(bào)精度下削峰率約為6.60%,當(dāng)τ≥6 時(shí)削峰率維持在12%~13%。而在相同洪水預(yù)見期下,隨著預(yù)報(bào)精度的提高,削峰率呈較小幅度的增大趨勢(shì),吉安站洪峰流量標(biāo)準(zhǔn)差越大。綜上所述,預(yù)報(bào)誤差降低了實(shí)時(shí)預(yù)報(bào)優(yōu)化調(diào)度的削峰效果,并增加了吉安站洪峰流量不確定性,但相比現(xiàn)行調(diào)度規(guī)則模擬結(jié)果仍具有較大優(yōu)勢(shì),研究成果具有實(shí)際意義。

    4 結(jié)論

    (1)萬安水庫現(xiàn)有的后汛期防洪調(diào)度規(guī)則針對(duì)1992 年7 月和2019年7月兩場(chǎng)50年一遇典型洪水過程攔洪削峰效果較弱,僅實(shí)現(xiàn)了不加重下游防洪壓力的效果。

    (2)若要結(jié)合預(yù)報(bào)信息開展防洪預(yù)報(bào)實(shí)時(shí)優(yōu)化調(diào)度,洪水預(yù)報(bào)有效信息至少要比洪水傳播時(shí)滯多2 個(gè)時(shí)段,即萬安水庫針對(duì)下游吉安站開展防洪預(yù)報(bào)優(yōu)化調(diào)度至少需要30 h 有效洪水預(yù)報(bào)信息。

    (3)結(jié)合預(yù)報(bào)信息,研究提出的預(yù)報(bào)優(yōu)化調(diào)度策略指導(dǎo)萬安水庫開展實(shí)時(shí)防洪調(diào)度,與確定性優(yōu)化結(jié)果存在一定差距,但相比現(xiàn)行調(diào)度規(guī)則提高削峰率約10%。

    猜你喜歡
    萬安吉安防洪
    快速組裝防洪擋水墻裝置
    夏季防洪防汛
    金吉安監(jiān)理公司市場(chǎng)開拓取得重大突破
    公益宣傳(防洪進(jìn)行時(shí))
    2019年河南省防洪除澇保護(hù)統(tǒng)計(jì)匯總表(本年達(dá)到)
    浙江萬安科技股份有限公司
    專用汽車(2020年2期)2020-04-08 10:57:50
    蘇萬安
    寶藏(2019年3期)2019-03-28 05:24:22
    蘇萬安 藏石欣賞
    寶藏(2018年12期)2019-01-29 01:51:16
    圣誕節(jié)的那雙手套
    跳高比賽中的意外
    18禁国产床啪视频网站| 最近最新中文字幕大全免费视频 | 国产xxxxx性猛交| tube8黄色片| 国产高清视频在线播放一区 | 在线观看一区二区三区激情| 精品视频人人做人人爽| 亚洲成人国产一区在线观看 | 91精品国产国语对白视频| 美国免费a级毛片| 王馨瑶露胸无遮挡在线观看| 免费在线观看视频国产中文字幕亚洲 | e午夜精品久久久久久久| 久9热在线精品视频| 又大又黄又爽视频免费| 十分钟在线观看高清视频www| 亚洲国产中文字幕在线视频| 免费在线观看黄色视频的| 亚洲国产av影院在线观看| 精品国产乱码久久久久久小说| 国产极品粉嫩免费观看在线| 中文字幕人妻熟女乱码| 久久这里只有精品19| 久久久久久亚洲精品国产蜜桃av| 99热网站在线观看| e午夜精品久久久久久久| 一区在线观看完整版| 最黄视频免费看| 免费在线观看黄色视频的| 大码成人一级视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品一区二区大全| 国产精品av久久久久免费| 欧美激情极品国产一区二区三区| 另类亚洲欧美激情| 免费在线观看影片大全网站 | 美女福利国产在线| 久久ye,这里只有精品| 精品久久久久久电影网| 叶爱在线成人免费视频播放| av不卡在线播放| 热99国产精品久久久久久7| 最近手机中文字幕大全| 国产极品粉嫩免费观看在线| 国产精品三级大全| 我的亚洲天堂| 日本av手机在线免费观看| 国产99久久九九免费精品| 色视频在线一区二区三区| 日本一区二区免费在线视频| 欧美激情高清一区二区三区| 国产精品二区激情视频| 精品人妻熟女毛片av久久网站| 人人妻人人澡人人爽人人夜夜| 国产av精品麻豆| 伊人亚洲综合成人网| 男人操女人黄网站| 日韩中文字幕视频在线看片| 中文字幕人妻丝袜一区二区| 高清视频免费观看一区二区| 男人爽女人下面视频在线观看| 亚洲成人免费电影在线观看 | 国产免费现黄频在线看| 巨乳人妻的诱惑在线观看| tube8黄色片| 久久精品熟女亚洲av麻豆精品| 日韩欧美一区视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久人人人人人| 亚洲 国产 在线| 黄色 视频免费看| 亚洲图色成人| 国产在线一区二区三区精| 乱人伦中国视频| 国产真人三级小视频在线观看| 性高湖久久久久久久久免费观看| 人人妻人人澡人人爽人人夜夜| 永久免费av网站大全| 久久人妻熟女aⅴ| 色94色欧美一区二区| 91国产中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 男女高潮啪啪啪动态图| 精品一区二区三区av网在线观看 | 91九色精品人成在线观看| 99久久人妻综合| 晚上一个人看的免费电影| 波多野结衣av一区二区av| 国产免费又黄又爽又色| 亚洲av美国av| 美女脱内裤让男人舔精品视频| 精品少妇黑人巨大在线播放| 国产高清视频在线播放一区 | 国产免费视频播放在线视频| 成人影院久久| 制服人妻中文乱码| 日韩熟女老妇一区二区性免费视频| 女性被躁到高潮视频| 美女主播在线视频| 性色av乱码一区二区三区2| 国产精品麻豆人妻色哟哟久久| 国产老妇伦熟女老妇高清| 日日摸夜夜添夜夜爱| 人人澡人人妻人| 国产av国产精品国产| 婷婷丁香在线五月| 国产成人免费观看mmmm| 操出白浆在线播放| 国产成人av教育| 一级,二级,三级黄色视频| 国产免费福利视频在线观看| 日韩欧美一区视频在线观看| 在线观看www视频免费| 美女视频免费永久观看网站| 少妇猛男粗大的猛烈进出视频| xxx大片免费视频| 精品国产乱码久久久久久小说| 后天国语完整版免费观看| 亚洲精品日韩在线中文字幕| 亚洲七黄色美女视频| 国产有黄有色有爽视频| 国产有黄有色有爽视频| 亚洲av电影在线进入| 久久中文字幕一级| 99国产综合亚洲精品| 亚洲国产精品国产精品| 久久av网站| 久久狼人影院| 亚洲欧美色中文字幕在线| 美女午夜性视频免费| 欧美 日韩 精品 国产| 成人国产av品久久久| 国产成人一区二区在线| 一二三四社区在线视频社区8| 午夜免费观看性视频| 精品亚洲乱码少妇综合久久| 在线观看www视频免费| 一本大道久久a久久精品| 婷婷色麻豆天堂久久| 国产成人系列免费观看| 黄色 视频免费看| a级毛片黄视频| 欧美 日韩 精品 国产| 久久久久久人人人人人| 午夜两性在线视频| 少妇的丰满在线观看| 可以免费在线观看a视频的电影网站| 久久久久久亚洲精品国产蜜桃av| 国产成人欧美在线观看 | 十八禁高潮呻吟视频| 久久人人97超碰香蕉20202| 亚洲精品久久成人aⅴ小说| 欧美成人精品欧美一级黄| xxxhd国产人妻xxx| 首页视频小说图片口味搜索 | 啦啦啦啦在线视频资源| 天天操日日干夜夜撸| 国产成人系列免费观看| 一区福利在线观看| 国产无遮挡羞羞视频在线观看| 国产成人免费无遮挡视频| 真人做人爱边吃奶动态| 精品国产超薄肉色丝袜足j| 最新在线观看一区二区三区 | 一级毛片 在线播放| 又粗又硬又长又爽又黄的视频| 久久久久久久久免费视频了| www日本在线高清视频| 伊人久久大香线蕉亚洲五| 国产在线观看jvid| 麻豆乱淫一区二区| 久热爱精品视频在线9| 中文字幕亚洲精品专区| 日本av手机在线免费观看| 国产97色在线日韩免费| 国产一区二区三区av在线| 国产在线视频一区二区| 国产精品久久久久久精品古装| 一区二区三区四区激情视频| 高清不卡的av网站| 在线观看免费高清a一片| 中文字幕高清在线视频| 精品少妇内射三级| 欧美精品av麻豆av| 精品国产一区二区久久| 国产精品一区二区在线观看99| 久久国产精品男人的天堂亚洲| 欧美 日韩 精品 国产| 黄色视频在线播放观看不卡| xxxhd国产人妻xxx| 18禁裸乳无遮挡动漫免费视频| 成人三级做爰电影| 亚洲精品乱久久久久久| 日韩伦理黄色片| 狂野欧美激情性xxxx| 男女高潮啪啪啪动态图| 亚洲精品久久午夜乱码| 97人妻天天添夜夜摸| 天天躁夜夜躁狠狠躁躁| 热99久久久久精品小说推荐| 成在线人永久免费视频| 两性夫妻黄色片| 欧美日韩福利视频一区二区| 日韩一卡2卡3卡4卡2021年| 一级,二级,三级黄色视频| 黄色片一级片一级黄色片| 久久久久精品国产欧美久久久 | 亚洲精品一卡2卡三卡4卡5卡 | 久久久精品94久久精品| 天天躁日日躁夜夜躁夜夜| 欧美日韩黄片免| 国产精品久久久人人做人人爽| 久久免费观看电影| 亚洲天堂av无毛| 欧美日韩福利视频一区二区| 亚洲一区二区三区欧美精品| 国产精品成人在线| 午夜av观看不卡| 最近中文字幕2019免费版| 男人舔女人的私密视频| 成人国产av品久久久| 亚洲欧美精品自产自拍| 国产精品二区激情视频| av天堂在线播放| www.自偷自拍.com| 丝瓜视频免费看黄片| 人妻 亚洲 视频| 成人午夜精彩视频在线观看| av在线老鸭窝| 自线自在国产av| 久久毛片免费看一区二区三区| 欧美亚洲日本最大视频资源| 涩涩av久久男人的天堂| 飞空精品影院首页| 美女视频免费永久观看网站| 极品少妇高潮喷水抽搐| 国产精品香港三级国产av潘金莲 | 亚洲国产精品成人久久小说| 国产成人免费无遮挡视频| 高清av免费在线| 91麻豆av在线| 脱女人内裤的视频| 91九色精品人成在线观看| a 毛片基地| 色播在线永久视频| 色视频在线一区二区三区| 国语对白做爰xxxⅹ性视频网站| 丝袜美足系列| 久久精品亚洲熟妇少妇任你| 国产爽快片一区二区三区| 国产亚洲av片在线观看秒播厂| 国产成人一区二区在线| 永久免费av网站大全| 免费观看a级毛片全部| 欧美中文综合在线视频| 久久人人爽av亚洲精品天堂| 黄色毛片三级朝国网站| 91九色精品人成在线观看| 晚上一个人看的免费电影| 无遮挡黄片免费观看| 看免费成人av毛片| 天天躁夜夜躁狠狠久久av| 亚洲欧美成人综合另类久久久| 国产精品一区二区在线不卡| 免费av中文字幕在线| 久久久久久亚洲精品国产蜜桃av| 久久热在线av| 丁香六月欧美| a 毛片基地| 建设人人有责人人尽责人人享有的| a级片在线免费高清观看视频| 久久亚洲国产成人精品v| 男女床上黄色一级片免费看| 99精国产麻豆久久婷婷| 国产伦人伦偷精品视频| 亚洲天堂av无毛| 亚洲色图 男人天堂 中文字幕| 亚洲激情五月婷婷啪啪| 亚洲天堂av无毛| 多毛熟女@视频| 色视频在线一区二区三区| 不卡av一区二区三区| 精品亚洲成a人片在线观看| 老汉色av国产亚洲站长工具| 欧美精品人与动牲交sv欧美| 欧美少妇被猛烈插入视频| 久久国产精品人妻蜜桃| 啦啦啦啦在线视频资源| av视频免费观看在线观看| 日本a在线网址| 国产一区二区三区av在线| 亚洲国产精品一区三区| 亚洲欧美一区二区三区国产| 国产日韩欧美在线精品| 欧美激情高清一区二区三区| 国产男人的电影天堂91| 啦啦啦在线免费观看视频4| 亚洲国产精品一区三区| 50天的宝宝边吃奶边哭怎么回事| 午夜91福利影院| 一区二区三区激情视频| 黄色毛片三级朝国网站| a 毛片基地| 十八禁高潮呻吟视频| 涩涩av久久男人的天堂| 在线天堂中文资源库| 视频区图区小说| 一区二区三区精品91| 亚洲色图 男人天堂 中文字幕| 大型av网站在线播放| 18禁黄网站禁片午夜丰满| 一级毛片 在线播放| 男女边吃奶边做爰视频| 国产高清国产精品国产三级| 精品高清国产在线一区| 晚上一个人看的免费电影| 久久久精品免费免费高清| 大香蕉久久网| 香蕉国产在线看| 女警被强在线播放| 免费久久久久久久精品成人欧美视频| 9191精品国产免费久久| 丁香六月天网| 精品国产一区二区三区久久久樱花| 在线av久久热| 9191精品国产免费久久| 亚洲第一av免费看| 国产伦人伦偷精品视频| 精品一区在线观看国产| 丝袜脚勾引网站| 欧美97在线视频| 欧美xxⅹ黑人| 黑人猛操日本美女一级片| 亚洲欧美中文字幕日韩二区| 欧美97在线视频| 成年美女黄网站色视频大全免费| 国产麻豆69| 日日夜夜操网爽| 老司机亚洲免费影院| 久久久久久久大尺度免费视频| 亚洲国产日韩一区二区| 国产有黄有色有爽视频| bbb黄色大片| 一级毛片女人18水好多 | 久久亚洲精品不卡| 电影成人av| 久久国产精品人妻蜜桃| 嫩草影视91久久| 国产成人精品在线电影| 精品久久蜜臀av无| 黄色视频在线播放观看不卡| 自线自在国产av| 人妻 亚洲 视频| 永久免费av网站大全| 色婷婷av一区二区三区视频| 我要看黄色一级片免费的| 欧美 日韩 精品 国产| 欧美av亚洲av综合av国产av| 免费观看av网站的网址| 日本wwww免费看| 久久青草综合色| 国产精品国产三级国产专区5o| av在线老鸭窝| 亚洲欧美精品自产自拍| 欧美另类一区| 午夜福利影视在线免费观看| 九色亚洲精品在线播放| 色综合欧美亚洲国产小说| 午夜老司机福利片| 丝袜美足系列| 欧美亚洲日本最大视频资源| 男女高潮啪啪啪动态图| 日韩 亚洲 欧美在线| 精品人妻熟女毛片av久久网站| 首页视频小说图片口味搜索 | 1024视频免费在线观看| 18禁裸乳无遮挡动漫免费视频| 校园人妻丝袜中文字幕| √禁漫天堂资源中文www| 男人操女人黄网站| 亚洲欧美日韩高清在线视频 | 日本欧美视频一区| 90打野战视频偷拍视频| 波多野结衣一区麻豆| 日日夜夜操网爽| 久久久久久久国产电影| 国产精品av久久久久免费| 日本午夜av视频| 人人妻人人添人人爽欧美一区卜| 欧美国产精品一级二级三级| 男女边吃奶边做爰视频| 免费日韩欧美在线观看| 男女无遮挡免费网站观看| 大型av网站在线播放| 精品久久蜜臀av无| 曰老女人黄片| 久久毛片免费看一区二区三区| 亚洲一区二区三区欧美精品| 久久这里只有精品19| netflix在线观看网站| 天堂俺去俺来也www色官网| 女性生殖器流出的白浆| 最新在线观看一区二区三区 | 日韩av在线免费看完整版不卡| 99热国产这里只有精品6| 青春草亚洲视频在线观看| 国产老妇伦熟女老妇高清| www.自偷自拍.com| 成在线人永久免费视频| 午夜久久久在线观看| 777久久人妻少妇嫩草av网站| 国产一卡二卡三卡精品| 久久性视频一级片| 日日夜夜操网爽| 黑人巨大精品欧美一区二区蜜桃| 18禁国产床啪视频网站| 亚洲av国产av综合av卡| 免费一级毛片在线播放高清视频 | 国产99久久九九免费精品| 亚洲成人免费电影在线观看 | 一区在线观看完整版| 啦啦啦在线观看免费高清www| 九色亚洲精品在线播放| 80岁老熟妇乱子伦牲交| 成在线人永久免费视频| 各种免费的搞黄视频| 一级a爱视频在线免费观看| 人人妻人人澡人人爽人人夜夜| 超色免费av| av在线播放精品| 超碰成人久久| 777米奇影视久久| 午夜免费鲁丝| 精品卡一卡二卡四卡免费| 久久精品久久久久久久性| 看免费成人av毛片| 香蕉国产在线看| 国产精品 欧美亚洲| 韩国高清视频一区二区三区| www.999成人在线观看| 亚洲色图 男人天堂 中文字幕| 丁香六月欧美| 人人妻,人人澡人人爽秒播 | 女性被躁到高潮视频| 午夜免费成人在线视频| 欧美激情极品国产一区二区三区| 色精品久久人妻99蜜桃| 51午夜福利影视在线观看| 大码成人一级视频| 免费在线观看视频国产中文字幕亚洲 | 五月天丁香电影| 少妇精品久久久久久久| 国产精品二区激情视频| 国产成人av教育| 老司机深夜福利视频在线观看 | 国产精品国产三级国产专区5o| 国产免费现黄频在线看| 国产精品免费大片| 飞空精品影院首页| 亚洲av欧美aⅴ国产| 一本久久精品| 精品国产一区二区三区久久久樱花| 高清不卡的av网站| 男女床上黄色一级片免费看| 欧美精品av麻豆av| 日韩欧美一区视频在线观看| 狂野欧美激情性xxxx| 十八禁高潮呻吟视频| 黄色片一级片一级黄色片| 国产一级毛片在线| 国产精品偷伦视频观看了| 久久国产精品影院| av网站免费在线观看视频| 欧美 亚洲 国产 日韩一| 国产av精品麻豆| 精品少妇内射三级| 欧美激情极品国产一区二区三区| 欧美国产精品一级二级三级| 久久精品aⅴ一区二区三区四区| 伊人亚洲综合成人网| av在线播放精品| 中文字幕亚洲精品专区| 人成视频在线观看免费观看| 久久国产精品人妻蜜桃| 欧美日韩亚洲高清精品| 精品人妻一区二区三区麻豆| 国产免费福利视频在线观看| 国产日韩一区二区三区精品不卡| 国产真人三级小视频在线观看| 国产色视频综合| 咕卡用的链子| 国产深夜福利视频在线观看| 天天躁夜夜躁狠狠久久av| 午夜福利视频精品| 黑人巨大精品欧美一区二区蜜桃| 99精品久久久久人妻精品| 亚洲精品一卡2卡三卡4卡5卡 | 日本色播在线视频| 好男人视频免费观看在线| 国产一区有黄有色的免费视频| 久久性视频一级片| 91麻豆精品激情在线观看国产 | 亚洲国产欧美网| 一区二区三区精品91| 在线观看免费日韩欧美大片| 成人国产一区最新在线观看 | 国产精品亚洲av一区麻豆| 一本久久精品| 国产三级黄色录像| 国产成人欧美在线观看 | 女性生殖器流出的白浆| 桃花免费在线播放| 一本久久精品| 亚洲视频免费观看视频| 欧美日韩亚洲综合一区二区三区_| 制服诱惑二区| 久久精品久久精品一区二区三区| 亚洲 国产 在线| 亚洲视频免费观看视频| 天天躁夜夜躁狠狠久久av| 精品高清国产在线一区| 一区福利在线观看| 一级毛片黄色毛片免费观看视频| 国产成人影院久久av| 日日夜夜操网爽| 久久毛片免费看一区二区三区| 亚洲成国产人片在线观看| 两个人看的免费小视频| 91九色精品人成在线观看| 久热这里只有精品99| 免费高清在线观看视频在线观看| 亚洲av国产av综合av卡| 我的亚洲天堂| 啦啦啦在线免费观看视频4| 久久这里只有精品19| 另类精品久久| 性色av乱码一区二区三区2| 欧美日韩一级在线毛片| 成年人午夜在线观看视频| 少妇被粗大的猛进出69影院| 后天国语完整版免费观看| 欧美人与性动交α欧美软件| 国产精品久久久人人做人人爽| 免费在线观看黄色视频的| 亚洲精品成人av观看孕妇| 人体艺术视频欧美日本| 黄片播放在线免费| 国产真人三级小视频在线观看| 只有这里有精品99| 99久久99久久久精品蜜桃| 天天添夜夜摸| 亚洲欧美日韩高清在线视频 | 久久精品aⅴ一区二区三区四区| 久久精品久久久久久噜噜老黄| 日韩 亚洲 欧美在线| 肉色欧美久久久久久久蜜桃| 我要看黄色一级片免费的| 亚洲少妇的诱惑av| 香蕉国产在线看| 老司机在亚洲福利影院| 久久国产精品大桥未久av| 天堂8中文在线网| 日韩欧美一区视频在线观看| 一区在线观看完整版| 无限看片的www在线观看| 十八禁网站网址无遮挡| av一本久久久久| 日本欧美视频一区| 国产成人免费观看mmmm| 免费少妇av软件| 黑人猛操日本美女一级片| 国产精品国产三级国产专区5o| 久久精品亚洲av国产电影网| 国产野战对白在线观看| 国产精品亚洲av一区麻豆| 国产在线观看jvid| 国产男人的电影天堂91| 国产精品二区激情视频| 亚洲中文日韩欧美视频| 国产高清不卡午夜福利| 丰满饥渴人妻一区二区三| 成年人黄色毛片网站| 亚洲五月色婷婷综合| 国产熟女午夜一区二区三区| 狠狠精品人妻久久久久久综合| 老司机午夜十八禁免费视频| 亚洲欧美日韩另类电影网站| 欧美黄色淫秽网站| av欧美777| 久久热在线av| 大陆偷拍与自拍| 50天的宝宝边吃奶边哭怎么回事| 最近中文字幕2019免费版| a级毛片黄视频| 亚洲成色77777| 丰满饥渴人妻一区二区三| 久久99一区二区三区| 久久久久久久大尺度免费视频| 久久九九热精品免费| 久热爱精品视频在线9| 午夜av观看不卡| 欧美+亚洲+日韩+国产| 精品欧美一区二区三区在线| 老汉色av国产亚洲站长工具| 国产成人精品久久二区二区免费| 男女之事视频高清在线观看 | 久久天堂一区二区三区四区| 啦啦啦中文免费视频观看日本| 国产在线免费精品| 久久人人爽人人片av| 观看av在线不卡| 999精品在线视频| av片东京热男人的天堂| 亚洲国产看品久久| 9191精品国产免费久久|