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

    基于時間序列分析的火電機組運行優(yōu)化研究

    2017-04-12 05:39:38朱龍霞肖明美許云峰
    河北工業(yè)科技 2017年2期
    關(guān)鍵詞:火電廠燃煤預(yù)測值

    朱龍霞,肖明美,栗 濤,許云峰,馬 瑞,劉 濱

    (1.河北科技大學(xué)信息科學(xué)與工程學(xué)院,河北石家莊 050018;2.中國核電工程有限公司河北分公司,河北石家莊 050000;3.河北科技大學(xué)經(jīng)濟管理學(xué)院,河北石家莊 050018)

    基于時間序列分析的火電機組運行優(yōu)化研究

    朱龍霞1,肖明美1,栗 濤2,許云峰1,馬 瑞1,劉 濱3

    (1.河北科技大學(xué)信息科學(xué)與工程學(xué)院,河北石家莊 050018;2.中國核電工程有限公司河北分公司,河北石家莊 050000;3.河北科技大學(xué)經(jīng)濟管理學(xué)院,河北石家莊 050018)

    為了預(yù)測火電機組在變工況狀態(tài)下的運行最優(yōu)值,提高燃煤電廠的燃煤效率,依據(jù)燃煤火電廠3個月的各項數(shù)據(jù),通過分析并挖掘燃煤火電廠與供電煤率密切相關(guān)的可調(diào)控參數(shù)(主蒸汽溫度、主蒸汽壓力)隨時間的變化關(guān)系,利用Holt-Winters時間序列法進行預(yù)測。經(jīng)過計算發(fā)現(xiàn),通過Holt-Winters方法得到的主蒸汽溫度和主蒸汽壓力的相對誤差值分別為0.28%和2.13%。實驗結(jié)果表明:Holt-Winters時間序列法對預(yù)測燃煤火電廠火電機組在變工況下的運行數(shù)據(jù)是合理的,且其誤差小,對燃煤火電廠提高火電機組的燃煤效率具有一定的指導(dǎo)意義。

    煤炭能;Holt-Winters;燃煤火電廠;數(shù)據(jù)挖掘;數(shù)據(jù)分析;預(yù)測

    火力發(fā)電廠在中國電力工業(yè)中占據(jù)十分重要的地位,發(fā)電量可以達到總發(fā)電量的80%,而在重視節(jié)能減排的當今時代,火電廠的煤炭消費仍然達到了中國煤炭總數(shù)的50%以上。燃煤火電廠是依靠汽水系統(tǒng)、燃燒系統(tǒng)、電氣系統(tǒng)和控制系統(tǒng)的共同作用而使燃料中的化學(xué)能量以電能的方式傳輸出去的。按照能量守恒原理,各個系統(tǒng)的能量變化都會波及到下一系統(tǒng),并由于負荷、供電質(zhì)量、經(jīng)濟性等要求,控制系統(tǒng)又會使這些波動反過來影響能量源——燃燒系統(tǒng)[1]。目前,中國通過改善燃煤方式[2-3]和鍋爐改造等方式提高發(fā)電效率,但這主要集中于性能優(yōu)化,而對供電電壓和電網(wǎng)電壓不平衡所造成的燃煤虛耗過多、電力浪費等問題沒有過多考慮。因此,本文從提高火電機組運行效率的角度,尋找當前機組在運行工況條件下所能達到的最佳參數(shù),從而為運行人員提供機組在特定負荷及外部條件下的最優(yōu)運行方式和參數(shù)控制。

    目前確定機組運行優(yōu)化目標值的研究方法可分為3種:1)制造廠所提供的標準值;2)以火電機組經(jīng)濟性為目標函數(shù),采用機組變工況理論優(yōu)化計算得到;3)機組變工況下,以實時運行過程中的動態(tài)優(yōu)化確定[4-6]。國內(nèi)大多采用前兩種方式,然而隨著機組運行時間的延長和機組工作狀態(tài)的動態(tài)變化,并不能反映機組實際運行狀態(tài)下的最優(yōu)化。

    許多學(xué)者從事機組運行優(yōu)化目標值的研究,根據(jù)現(xiàn)代優(yōu)化計算方法,如關(guān)聯(lián)規(guī)則算法[7]、神經(jīng)網(wǎng)絡(luò)算法[8-10]、遺傳算法[11]、時間序列法[12-13]、目標函數(shù)計算法[4-5]等建立相關(guān)模型,從而獲得該機組運行優(yōu)化的最優(yōu)值。燃煤火電廠的任務(wù)之一是為電網(wǎng)用戶提供足夠的滿足質(zhì)量要求的熱負荷,由于各種原因,電網(wǎng)用戶所需的熱負荷是隨時變化的,當熱電廠提供的熱負荷與電網(wǎng)用戶的需求量相平衡時最經(jīng)濟,否則,會造成浪費,或使電網(wǎng)用戶生產(chǎn)難以正常運行。本文利用燃煤火電廠數(shù)據(jù)的時間序列性,對數(shù)以千計的燃煤火電廠某機組的實際運行數(shù)據(jù)進行分析,為運行機組人員提供在各種擾動下的最優(yōu)運行參數(shù)控制,從而提高燃煤電廠的發(fā)電效率。

    1 關(guān)鍵技術(shù)描述

    火電廠運行設(shè)備之間緊密聯(lián)系,各設(shè)備運行狀態(tài)的波動互相關(guān)聯(lián),這種動態(tài)波動的關(guān)聯(lián)情況通過與設(shè)備狀態(tài)相關(guān)的運行數(shù)據(jù)明顯地顯現(xiàn)出來[14]。電廠數(shù)據(jù)具有種類多、數(shù)量大、高維和強耦合的特點,且其在多個子系統(tǒng)的工作過程中,是隨時間的變化而變化的。本文采用socket通信方式并以2 h為一個時間間隔,連續(xù)從分散控制系統(tǒng)(DCS)下載了3個月實時數(shù)據(jù)并進行分析。通過對電廠數(shù)據(jù)進行合法性檢查,剔除冗余數(shù)據(jù)、空白數(shù)據(jù)和錯誤數(shù)據(jù)后,將其根據(jù)控制系統(tǒng)來源進行分組并在組中添加數(shù)據(jù),最后將這些數(shù)據(jù)保存到XML(eXtensible markup language,可擴展標記語言)中。由于XML中的數(shù)據(jù)是從電廠獲得的所有類型數(shù)據(jù),故利用dom4j解析XML,將與該機組供電煤率密切相關(guān)的可調(diào)控參數(shù),包括主蒸汽壓力、主蒸汽溫度、減溫水量、給水溫度、過量空氣系數(shù)、再熱減溫水量等數(shù)據(jù)提取出來并構(gòu)建數(shù)據(jù)倉庫[15]。

    時間序列,也叫時間數(shù)列、歷史復(fù)數(shù)或動態(tài)數(shù)列,是將某種統(tǒng)計指標的數(shù)值按時間先后順序排列所形成的數(shù)列[12,16]。對于平穩(wěn)時間序列分析,指數(shù)平滑法中的Holt雙參數(shù)線性指數(shù)平滑法適用于線性趨勢、無季節(jié)變化的時間序列預(yù)測;Winters線性和季節(jié)性指數(shù)平滑法適用于對含有季節(jié)性因素的時間序列的預(yù)測[13,17]。Holt-Winters模型是通過對歷史數(shù)據(jù)的加權(quán)平均值來跟蹤數(shù)據(jù)的變化,其可分為季節(jié)加性和季節(jié)乘性。

    2 時間序列模型的確立與分析

    任何一個時間序列模型的建立首先都需要對時間序列數(shù)據(jù)作走勢圖,在對其走勢進行分析的基礎(chǔ)上,結(jié)合數(shù)據(jù)信息的先驗知識來進行時間序列模型的確立。以主蒸汽溫度、主蒸汽壓力為例,利用Matlab對其歷史數(shù)據(jù)進行了走勢圖的繪制,得到其時間序列的數(shù)據(jù)概貌,如圖1和圖2所示。

    圖1 主蒸汽溫度隨時間變化關(guān)系圖Fig.1 Relation diagram of main stream temperature with datetime

    圖2 主蒸汽壓力隨時間變化關(guān)系圖Fig.2 Relation diagram of vapor pressure with datetime

    從圖中可以看出,雖然各項數(shù)據(jù)隨時間的變化而變化,且變化程度有強有弱,但總的來說都有趨勢和“季節(jié)”變化特征,故可采用具有趨勢和季節(jié)的時間序列模型對其進行預(yù)測。

    由于無法從圖中看出明顯的遞增或遞減趨勢,本文用季節(jié)性變動分析方法中的移動平均數(shù)比率法對上述幾項時間序列樣本進行計算

    通過計算可以發(fā)現(xiàn),該值在允許誤差值范圍內(nèi)保持恒定,且兩者之間呈比例性而非數(shù)值增加或減少的相互關(guān)系,從而推斷出上述數(shù)據(jù)時間序列中的“季節(jié)效應(yīng)”具有成倍增長或下降的表現(xiàn),故采用Holt-Winters乘法預(yù)測模型對其監(jiān)測指標進行預(yù)測。

    基于以上研究,對該機組下的主蒸汽溫度和主蒸汽壓力構(gòu)建Holt-Winters乘法模型進行預(yù)測,預(yù)測間隔h個時間段的主蒸汽溫度和主蒸汽壓力數(shù)據(jù)的預(yù)測方程為

    F=St-2+h(mt+hrt) ,

    式中:F表示預(yù)測值;mt表示在時間間隔t的估計平均值;St表示適合于該時間段的“季節(jié)”因子;rt表示在時間間隔t的估計趨勢項,即各數(shù)值對現(xiàn)數(shù)據(jù)平均值的期望增加或減少量。這3個系數(shù)確立所依據(jù)的方程組為

    mt=αxt/St-2+(1-a)(mt-1+rt-1) ,

    St=βxt/mt+(1-β)St-2,

    rt=γ(mt-mt-1)+(1-γ)rt-1,

    式中:α,β,γ為Holt-Winters時間序列分析模型的平滑常數(shù),其值在0~1之間變化;“季節(jié)”變量因子為2。該方程通過不斷修正平滑值St和趨勢值rt,從而使預(yù)測值更加接近真實值。

    3 Holt-Winters預(yù)測值計算與回歸分析模型預(yù)測精度對比分析

    EViews(econometrics views)通常稱為計量經(jīng)濟學(xué)軟件包,是對社會經(jīng)濟關(guān)系與經(jīng)濟活動的數(shù)量規(guī)律,采用計量經(jīng)濟學(xué)方法與技術(shù)進行“觀察”。其在科學(xué)數(shù)據(jù)分析與評價、金融分析、經(jīng)濟預(yù)測、銷售預(yù)測和成本分析等領(lǐng)域應(yīng)用非常廣泛[18]。通過采用Holt-Winters乘法預(yù)測模型與回歸分析模型進行預(yù)測,經(jīng)過多次重復(fù)實驗,其中Holt-Winters乘法模型中,主蒸汽溫度權(quán)重系數(shù)α,β,γ分別為0.2,0.5,0.8,主蒸汽壓力權(quán)重系數(shù)α,β,γ分別為0.9,0.1,0.05時效果達到最好;通過對3個月電廠數(shù)據(jù)的整理與組合,主蒸汽溫度的回歸模型為Y= 536.544 666 667 + 0.061 030 303 030 3×X,主蒸汽壓力的回歸模型為Y=16.060 666 666 7+0.025 757 575 757 6×X,其中Y為預(yù)測值,X為次數(shù)。其實際值與分別采用兩種模型進行預(yù)測的預(yù)測值如表1所示。

    表1 X日監(jiān)測燃煤火電廠運行數(shù)據(jù)(實際值/預(yù)測值/使用回歸分析模型預(yù)測值)

    本文根據(jù)表1分別做出主蒸汽溫度的實際值、預(yù)測值與回歸分析模型的預(yù)測值曲線圖,以及主蒸汽壓力的實際值、預(yù)測值與回歸分析模型的預(yù)測值曲線圖,如圖3和圖4所示。

    圖3 主蒸汽溫度的預(yù)測值、實際值與回歸分析預(yù)測值折線圖Fig.3 Line chart of main stream temperature’s predictions actual values, and predictions by regression analysis method

    圖4 主蒸汽壓力的預(yù)測值、實際值和回歸分析預(yù)測值折線圖Fig.4 Line chart of vapor pressure’s predictions, actual values and predictions by regression analysis method

    式中:vi表示實際值;pi表示預(yù)測值;n表示需要預(yù)測項目的數(shù)量;v′表示實際值的平均值;σ即為相對誤差值。

    根據(jù)上述公式,通過Holt-Winters方法得到的主蒸汽溫度和主蒸汽壓力的相對誤差值分別為0.28%,2.13%;通過回歸分析模型得到的主蒸汽溫度和主蒸汽壓力的相對誤差值分別為0.43%,6.36%。通過選取回歸分析模型與本文預(yù)測方法進行對比,以平均相對誤差作為模型預(yù)測效果評價指標,結(jié)果發(fā)現(xiàn)本文預(yù)測方法在預(yù)測精度上明顯高于回歸分析模型預(yù)測方法,更適用于中短期燃煤火電廠機組數(shù)據(jù)預(yù)測。

    4 結(jié) 論

    1)根據(jù)一天之內(nèi)某燃煤火電機組在變工況下的實際運行數(shù)據(jù)的線性關(guān)系可以看出,主蒸汽溫度與主蒸汽壓力在實際運行過程中有增長趨勢,且Holt-Winters算法對其進行預(yù)測是合理的。

    2)根據(jù)文中的誤差分析可以得知,Holt-Winters算法相對于回歸線性模型對預(yù)測燃煤電廠火電機組的運行數(shù)據(jù)具有更高的準確性。

    3)燃煤電廠數(shù)據(jù)在3個月內(nèi)呈現(xiàn)增長趨勢,但在現(xiàn)實生活中燃煤電廠數(shù)據(jù)處于波動狀態(tài),所以線性回歸模型只能呈現(xiàn)燃煤電廠數(shù)據(jù)的總體發(fā)展趨勢,但對預(yù)測具體數(shù)據(jù)將存在很大誤差,而Holt-Winters算法可通過調(diào)整參數(shù)減少預(yù)測數(shù)據(jù)誤差。

    4)Holt-Winters時間序列分析法對進一步提高燃煤電廠火電機組運行的經(jīng)濟性具有一定的指導(dǎo)意義。

    /References:

    [1] 裘國相, 汪思源, 王文標. 數(shù)據(jù)挖掘在火電廠運行優(yōu)化目標值確定中的應(yīng)用[J]. 自動化技術(shù)與應(yīng)用, 2006, 25(3): 6-9. QIU Guoxiang, WANG Siyuan,WANG Wenbiao. Data mining in optimization of the targeted value for thermal power plant[J]. Techniques of Automation and Applications, 2006,25(3): 6-9.

    [2] 姬中國, 郭玉泉. 在建大型火力發(fā)電廠節(jié)能分析與對策[J]. 電力建設(shè), 2007, 28(5): 44-46. JI Zhongguo, GUO Yuquan. Analysis and countermeasures of under-constructed large thermal power plant energy saving[J]. Electric Power Construction, 2007, 28(5): 44-46.

    [3] 蔣敏華, 黃斌. 燃煤發(fā)電技術(shù)發(fā)展展望[J]. 中國電機工程學(xué)報, 2012,32(29): 1-8. JIANG Minhua, HUANG Bin. Prospects on coal-fired power generation technology development[J]. Proceedings of the CSEE, 2012,32(29): 1-8.

    [4] 胡洪華, 黃廷輝, 艾衛(wèi)國, 等. 大型火電機組運行優(yōu)化目標值的研究和確定[J]. 中國電力, 2004,37(9): 22-25. HU Honghua,HUANG Tinghui,AI Weiguo, et al. Study and determination of operation optimization target values of large size fossil-fired power unit[J]. Electric Power, 2004,37(9): 22-25.

    [5] 任浩仁, 李蔚, 盛德仁, 等. 火電機組變工況下運行指標應(yīng)達值的分析[J]. 中國電機工程學(xué)報, 1999,19(9): 50-52. REN Haoren, LI Wei, SHENG Deren, et al.The analyze of operation index for the power unit under different loads[J]. Proceedings of the CSEE, 1999,19(9): 50-52.

    [6] 李建強, 劉吉臻, 張欒英, 等. 基于數(shù)據(jù)挖掘的電站運行優(yōu)化應(yīng)用研究[J]. 中國電機工程學(xué)報, 2006,26(20): 118-123. LI Jianqiang, LIU Jizhen, ZHANG Luanying, et al. The research and application of fuzzy association rule mining in power plant operation optimization[J]. Proceedings of the CSEE, 2006,26(20): 118-123.

    [7] 鄭西西. 基于關(guān)聯(lián)規(guī)則的火電廠優(yōu)化目標值確定的研究[D]. 保定:華北電力大學(xué), 2011. ZHENG Xixi. Reaserch on the Thermal Power Plant Operation Optimization Value Determining Based on Association Rule[D]. Baoding: North China Electric Power University,2011.

    [8] 程偉良, 夏國棟, 徐壽臣, 等. 基于人工神經(jīng)網(wǎng)絡(luò)和模擬進化算法的分級燃燒優(yōu)化[J]. 清華大學(xué)學(xué)報(自然科學(xué)版), 2005, 45(5): 693-696. CHENG Weiliang, XIA Guodong, XU Shouchen, et al. Staged combustion oprimization based on artificial neural network and simulation evolvement algorithm[J]. Journal of Tsinghua University(Science and Technology), 2005, 45(5):693-696.

    [9] 周昊, 朱洪波, 岑可法. 基于人工神經(jīng)網(wǎng)絡(luò)和遺傳算法的火電廠鍋爐實時燃燒優(yōu)化系統(tǒng)[J]. 動力工程學(xué)報, 2003,23(5): 2665-2669. ZHOU Hao, ZHU Hongbo, CEN Kefa. An on-line boiler operating optimization system based on the neural network and the genetic algorithms[J]. Power Engineering,2003,23(5): 2665-2669.

    [10]王子杰, 李健, 孫萬云. 基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的鍋爐燃燒優(yōu)化方法[J]. 華北電力大學(xué)學(xué)報(自然科學(xué)版), 2008, 35(1): 14-17. WANG Zijie, LI Jian, SUN Wanyun. Boiler combusion optimization based on neural network and genetic algorithm[J]. Journal of North China Electric Power University(Natural Science Edition),2008,35(1):14-17.

    [11]蔡杰進, 馬曉茜. 基于遺傳算法的火電機組運行優(yōu)化[J]. 華南理工大學(xué)學(xué)報(自然科學(xué)版), 2006, 34(4): 101-105. CAI Jiejin, MA Xiaoqian. Operation oprimization of steam power unit based on genetic algorithm[J]. Joumal of South China University of Technology(Natural Science Edition), 2006,34(4): 101-105.

    [12]張美英, 何杰. 時間序列預(yù)測模型研究簡介[J]. 江西科學(xué), 2009,27(5): 697-701. ZHANG Meiying, HE Jie. Overview of time series forecasting model[J]. Jiangxi Science, 2009,27(5): 697-701.

    [13]吳濤, 顏輝武, 唐桂剛. 三峽庫區(qū)水質(zhì)數(shù)據(jù)時間序列分析預(yù)測研究[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2006,31(6): 500-502. WU Tao, YAN Huiwu, TANG Guigang. Prediction on time series analysis of water quality in Yangtze Gorger reservoir area[J]. Geomatics and Information Science of Wuhan University,2006,31(6): 500-502.

    [14]鮑文, 于達仁, 王偉, 等. 基于關(guān)聯(lián)規(guī)則的火電廠傳感器故障檢測[J]. 中國電機工程學(xué)報, 2003,23(12): 170-174. BAO Wen, YU Daren, WANG Wei, et al. Sensor fault detection in thermal power plants based on association rule[J]. Proceeding of the CSEE, 2003,23(12): 170-174.

    [15]劉濱, 吳燕, 劉振宇. 數(shù)據(jù)倉庫技術(shù)[J]. 河北工業(yè)科技, 2005,22(5): 305-307. LIU Bin, WU Yan, LIU Zhenyu. Introduction of data warehouse[J]. Hebei Journal of Industrial Science and Techno-logy,2005,22(5): 305-307.

    [16]KALEKAR P S. Time series forecasting using Holt-Winters exponential smoothing[J]. Kanwal Rekhi School of Information Technology, 2004:1-13.

    [17]劉羅曼. 時間序列分析中指數(shù)平滑法的應(yīng)用[J]. 沈陽師范大學(xué)學(xué)報(自然科學(xué)版), 2009,27(4): 416-418. LIU Luoman. Application of exponential smoothing in time series analysis[J]. Journal of Shenyang Normal University(Natural Science Edition), 2009,27(4): 416-418.

    [18]孫紅果, 余星, 王竟竟. EViews在概率統(tǒng)計教學(xué)中應(yīng)用[J]. 邵陽學(xué)院學(xué)報(自然科學(xué)版), 2013,10(3): 24-28. SUN Hongguo, YU Xing, WANG Jingjing. The application of EViews software in probability and statistics[J].Journal of Shaoyang University(Natural Science Edition), 2013,10(3): 24-28.

    Research of operation optimization of thermal power unit based on time series analysis

    ZHU Longxia1, XIAO Mingmei1, LI Tao2, XU Yunfeng1, MA Rui1, LIU Bin3

    (1.School of Information Science and Engineering,Hebei University of Science and Technology, Shijiazhuang, Hebei 050018, China;2.Hebei Branch, China Nuclear Power Engineering Company Limited, Shijiazhuang, Hebei 050000, China;3. School of Economics and Management, Hebei University of Science and Technology, Shijiazhuang, Hebei 050018, China)

    In order to predict the thermal power generating units in operation status under the optimal value, improve the efficiency of coal-fired power plant, according to the data of three months, through the analysis and mining of coal-fired power plants and coal supply rate closely related to the control parameters (this paper mainly introduces the main steam temperature, steam pressure) with time change, it's forecasted by using Holt-Winters time series. It is found that the relative error of the main steam temperature and the main steam pressure obtained by the Holt-Winters method is 0.28% and 2.13%. The experimental results show that the Holt-Winters time series method to predict the operation data of thermal power units in coal-fired power plant under changing condition is reasonable, and the error is small, providing a certain significance to improve the efficiency of coal-fired thermal power units for coal-fired power plants.

    coal energy; Holt-Winters; coal-fired power plant; data mining; data analysis; prediction

    1008-1534(2017)02-0125-05

    2016-11-08;

    2017-03-09;責任編輯:陳書欣

    國家自然科學(xué)基金(71271076)

    朱龍霞(1992—),女,河北石家莊人,碩士研究生,主要從事數(shù)據(jù)分析方面的研究。

    許云峰副教授。E-mail: hbkd_xyf@hebust.edu.cn

    TP391

    A

    10.7535/hbgykj.2017yx02009

    朱龍霞,肖明美,栗 濤,等.基于時間序列分析的火電機組運行優(yōu)化研究[J].河北工業(yè)科技,2017,34(2):125-129. ZHU Longxia, XIAO Mingmei, LI Tao,et al.Research of operation optimization of thermal power unit based on time series analysis[J].Hebei Journal of Industrial Science and Technology,2017,34(2):125-129.

    猜你喜歡
    火電廠燃煤預(yù)測值
    IMF上調(diào)今年全球經(jīng)濟增長預(yù)期
    企業(yè)界(2024年8期)2024-07-05 10:59:04
    加拿大農(nóng)業(yè)部下調(diào)2021/22年度油菜籽和小麥產(chǎn)量預(yù)測值
    ±800kV直流輸電工程合成電場夏季實測值與預(yù)測值比對分析
    法電再次修訂2020年核發(fā)電量預(yù)測值
    國外核新聞(2020年8期)2020-03-14 02:09:19
    全球最大火電廠的清潔改造
    能源(2017年10期)2017-12-20 05:54:15
    歐盟新規(guī)或?qū)?dǎo)致三分之一的燃煤電廠關(guān)閉
    火電廠循環(huán)水處理及系統(tǒng)優(yōu)化探討
    燃煤機組“超低排放”改造中CEMS的選型與應(yīng)用
    “超低”排放技術(shù)在我國燃煤電廠的應(yīng)用
    燃煤電廠節(jié)能管理
    成人亚洲精品av一区二区| 亚洲国产中文字幕在线视频| 国产伦精品一区二区三区四那| 久久国产精品影院| 亚洲片人在线观看| 久久国产乱子伦精品免费另类| 欧美又色又爽又黄视频| 又黄又爽又免费观看的视频| 精品久久久久久久久久免费视频| 最新在线观看一区二区三区| 日本 欧美在线| 性色avwww在线观看| 免费看十八禁软件| 日本a在线网址| 亚洲av电影在线进入| 老鸭窝网址在线观看| aaaaa片日本免费| 成人欧美大片| 午夜免费成人在线视频| 亚洲av成人av| 久久九九热精品免费| 天堂影院成人在线观看| 国产极品精品免费视频能看的| 国产精品久久久久久精品电影| 亚洲欧美日韩卡通动漫| 国内久久婷婷六月综合欲色啪| 岛国在线观看网站| 欧美av亚洲av综合av国产av| 国产成人精品无人区| 亚洲人成伊人成综合网2020| 国产1区2区3区精品| 色哟哟哟哟哟哟| 97超视频在线观看视频| tocl精华| 日韩欧美在线二视频| 99国产精品一区二区三区| 国产爱豆传媒在线观看| 国产一区在线观看成人免费| 国产精品精品国产色婷婷| aaaaa片日本免费| 少妇的逼水好多| 国产精品综合久久久久久久免费| 淫秽高清视频在线观看| e午夜精品久久久久久久| 亚洲成av人片在线播放无| 国产乱人视频| av在线天堂中文字幕| 国内精品久久久久久久电影| 国产免费男女视频| 久久亚洲精品不卡| 国产精华一区二区三区| 天堂av国产一区二区熟女人妻| 身体一侧抽搐| 久久婷婷人人爽人人干人人爱| 色综合亚洲欧美另类图片| 亚洲av成人av| 在线免费观看的www视频| 国产亚洲欧美在线一区二区| 免费搜索国产男女视频| a级毛片在线看网站| tocl精华| 久久精品国产99精品国产亚洲性色| av黄色大香蕉| 免费观看精品视频网站| 国产aⅴ精品一区二区三区波| 亚洲第一电影网av| 久久国产精品人妻蜜桃| 精品久久久久久久久久免费视频| 亚洲av成人精品一区久久| 成人一区二区视频在线观看| 可以在线观看毛片的网站| 91九色精品人成在线观看| 日韩国内少妇激情av| 身体一侧抽搐| 国产一区二区三区在线臀色熟女| 婷婷丁香在线五月| 一个人看的www免费观看视频| 一区福利在线观看| 国产精品久久视频播放| 怎么达到女性高潮| 男女午夜视频在线观看| 老鸭窝网址在线观看| 久久久久精品国产欧美久久久| 国产亚洲精品久久久久久毛片| 色吧在线观看| 熟女人妻精品中文字幕| 日日摸夜夜添夜夜添小说| 麻豆久久精品国产亚洲av| 国产av不卡久久| 国产成人aa在线观看| 黑人操中国人逼视频| av天堂在线播放| 我要搜黄色片| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一卡2卡三卡4卡5卡| 99久久精品热视频| 成人性生交大片免费视频hd| 99国产精品一区二区三区| 少妇丰满av| 在线观看午夜福利视频| 天堂av国产一区二区熟女人妻| 国产精品久久久久久亚洲av鲁大| 国产成人精品久久二区二区91| 蜜桃久久精品国产亚洲av| 三级男女做爰猛烈吃奶摸视频| 色吧在线观看| 欧美色欧美亚洲另类二区| 久久精品亚洲精品国产色婷小说| 免费电影在线观看免费观看| 亚洲色图 男人天堂 中文字幕| 国产激情欧美一区二区| 国产精品影院久久| 岛国在线免费视频观看| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久久电影 | 亚洲国产欧美网| 久久精品人妻少妇| 制服丝袜大香蕉在线| 人妻久久中文字幕网| 三级男女做爰猛烈吃奶摸视频| 88av欧美| 首页视频小说图片口味搜索| 日韩中文字幕欧美一区二区| 国产一区二区三区在线臀色熟女| 色综合欧美亚洲国产小说| 国产精品一区二区精品视频观看| 国产精品日韩av在线免费观看| xxxwww97欧美| 在线a可以看的网站| 黄片小视频在线播放| 熟女人妻精品中文字幕| 精品福利观看| 亚洲国产欧美人成| 欧美成人免费av一区二区三区| h日本视频在线播放| 香蕉av资源在线| 麻豆国产av国片精品| 黑人操中国人逼视频| 一区福利在线观看| 成年人黄色毛片网站| 亚洲中文日韩欧美视频| 久久精品国产亚洲av香蕉五月| 日本一二三区视频观看| 真人做人爱边吃奶动态| 国产蜜桃级精品一区二区三区| 真人一进一出gif抽搐免费| 波多野结衣巨乳人妻| 少妇人妻一区二区三区视频| 国产淫片久久久久久久久 | 三级毛片av免费| 日韩国内少妇激情av| 他把我摸到了高潮在线观看| 国产真实乱freesex| 在线观看舔阴道视频| 男人舔女人下体高潮全视频| 女警被强在线播放| 黄色成人免费大全| 亚洲av日韩精品久久久久久密| 久久久久国产一级毛片高清牌| 九九在线视频观看精品| 国产精品精品国产色婷婷| 欧美日韩亚洲国产一区二区在线观看| www日本在线高清视频| 一进一出好大好爽视频| 国产69精品久久久久777片 | 黄片大片在线免费观看| 中文资源天堂在线| 岛国视频午夜一区免费看| 三级毛片av免费| 香蕉久久夜色| 在线永久观看黄色视频| 亚洲欧美激情综合另类| 国产精品1区2区在线观看.| 亚洲精品美女久久久久99蜜臀| 免费在线观看成人毛片| 精品不卡国产一区二区三区| 琪琪午夜伦伦电影理论片6080| 狂野欧美激情性xxxx| 中文亚洲av片在线观看爽| 琪琪午夜伦伦电影理论片6080| 91麻豆av在线| 欧美+亚洲+日韩+国产| 亚洲av熟女| 欧美一级a爱片免费观看看| 久久欧美精品欧美久久欧美| 午夜免费激情av| 99热精品在线国产| 亚洲精品一卡2卡三卡4卡5卡| 不卡一级毛片| 久久婷婷人人爽人人干人人爱| 我要搜黄色片| 国产单亲对白刺激| 熟妇人妻久久中文字幕3abv| 又粗又爽又猛毛片免费看| 国产欧美日韩一区二区三| 国产成人精品无人区| tocl精华| 欧美黄色淫秽网站| 男女那种视频在线观看| 国产精品99久久久久久久久| 国产极品精品免费视频能看的| 午夜福利18| 亚洲精品美女久久av网站| 欧美日韩乱码在线| 最近最新免费中文字幕在线| 97超级碰碰碰精品色视频在线观看| 精品国内亚洲2022精品成人| 久久香蕉精品热| 国产不卡一卡二| 精品一区二区三区视频在线 | 亚洲精品在线观看二区| 美女cb高潮喷水在线观看 | 免费搜索国产男女视频| 99久久精品一区二区三区| 少妇丰满av| av女优亚洲男人天堂 | 国模一区二区三区四区视频 | 欧美成人性av电影在线观看| 久久久久久久久中文| 国产伦精品一区二区三区视频9 | 蜜桃久久精品国产亚洲av| 久久久久久九九精品二区国产| 欧美成人性av电影在线观看| 亚洲国产日韩欧美精品在线观看| 最近最新中文字幕免费大全7| 美女高潮的动态| 免费观看在线日韩| 亚洲精品影视一区二区三区av| 欧美变态另类bdsm刘玥| 七月丁香在线播放| 国产91av在线免费观看| 欧美一区二区精品小视频在线| 亚洲国产精品成人综合色| 国产激情偷乱视频一区二区| 最新中文字幕久久久久| 国产中年淑女户外野战色| 精品国产一区二区三区久久久樱花 | 又粗又硬又长又爽又黄的视频| 亚洲欧美精品综合久久99| 午夜福利在线观看免费完整高清在| 一个人免费在线观看电影| 真实男女啪啪啪动态图| 成人国产麻豆网| 久久久久性生活片| 国产爱豆传媒在线观看| 欧美成人一区二区免费高清观看| 精品少妇黑人巨大在线播放 | 国产午夜精品论理片| 免费不卡的大黄色大毛片视频在线观看 | 少妇被粗大猛烈的视频| 国产精品国产三级国产专区5o | 99热6这里只有精品| 欧美最新免费一区二区三区| 床上黄色一级片| 亚洲精品一区蜜桃| 色噜噜av男人的天堂激情| 亚洲欧美一区二区三区国产| 老司机影院毛片| 美女xxoo啪啪120秒动态图| 国产成人a区在线观看| 亚洲精品日韩av片在线观看| 国产真实乱freesex| 黄色日韩在线| 岛国在线免费视频观看| 内射极品少妇av片p| 九色成人免费人妻av| 内地一区二区视频在线| 日产精品乱码卡一卡2卡三| 日本欧美国产在线视频| 小说图片视频综合网站| 国产成人a∨麻豆精品| 亚洲无线观看免费| 91久久精品国产一区二区三区| 丝袜喷水一区| 亚洲av男天堂| 99久久人妻综合| 深夜a级毛片| 国产成人aa在线观看| 精品国产露脸久久av麻豆 | videossex国产| 国产乱人偷精品视频| 久久国产乱子免费精品| 国产人妻一区二区三区在| 在线a可以看的网站| 久久99热这里只频精品6学生 | 国产私拍福利视频在线观看| 亚洲人成网站在线播| 亚洲精品久久久久久婷婷小说 | 国产精品国产高清国产av| 国产精品人妻久久久影院| 亚洲欧美一区二区三区国产| 国产午夜精品论理片| 日本色播在线视频| 亚洲欧美日韩卡通动漫| 直男gayav资源| 欧美最新免费一区二区三区| 我的女老师完整版在线观看| 97热精品久久久久久| 天堂影院成人在线观看| 人人妻人人澡人人爽人人夜夜 | 人妻制服诱惑在线中文字幕| 26uuu在线亚洲综合色| 最近中文字幕高清免费大全6| 欧美成人精品欧美一级黄| 中国美白少妇内射xxxbb| 国产三级中文精品| 亚洲av电影在线观看一区二区三区 | 国产高清不卡午夜福利| 亚州av有码| 在线免费十八禁| 精品久久久久久久人妻蜜臀av| 欧美三级亚洲精品| 久久99热这里只有精品18| 性色avwww在线观看| 国产在线一区二区三区精 | 色吧在线观看| 国产在视频线精品| 成人综合一区亚洲| 在线观看av片永久免费下载| 激情 狠狠 欧美| 久久99热6这里只有精品| 九草在线视频观看| 在线免费观看的www视频| av线在线观看网站| 国产伦精品一区二区三区视频9| 极品教师在线视频| 欧美日韩精品成人综合77777| 青春草国产在线视频| 久久久亚洲精品成人影院| 中文乱码字字幕精品一区二区三区 | 插逼视频在线观看| 国模一区二区三区四区视频| 精品不卡国产一区二区三区| 久久午夜福利片| 免费无遮挡裸体视频| 嫩草影院新地址| av在线亚洲专区| 久久久精品欧美日韩精品| 秋霞伦理黄片| 国产av不卡久久| 亚洲av成人精品一区久久| www.色视频.com| 高清毛片免费看| 男女下面进入的视频免费午夜| 免费看日本二区| 日本与韩国留学比较| 国产成人aa在线观看| 一区二区三区高清视频在线| 国产精品久久久久久精品电影小说 | 免费观看性生交大片5| 91狼人影院| 噜噜噜噜噜久久久久久91| 国产精品野战在线观看| 九九热线精品视视频播放| 啦啦啦韩国在线观看视频| 亚洲国产欧洲综合997久久,| 成人综合一区亚洲| 国产一区二区亚洲精品在线观看| 波多野结衣巨乳人妻| 欧美一区二区国产精品久久精品| 中文字幕av在线有码专区| 国产亚洲午夜精品一区二区久久 | 天堂影院成人在线观看| 久久久久久久久大av| 九九爱精品视频在线观看| 国产不卡一卡二| 欧美成人精品欧美一级黄| 五月伊人婷婷丁香| 免费黄网站久久成人精品| 久久精品久久久久久噜噜老黄 | av又黄又爽大尺度在线免费看 | 国产精品av视频在线免费观看| 欧美极品一区二区三区四区| 久久久成人免费电影| 久久精品国产鲁丝片午夜精品| 亚洲av成人精品一二三区| 国产久久久一区二区三区| 婷婷六月久久综合丁香| 高清毛片免费看| 色综合亚洲欧美另类图片| 99久久成人亚洲精品观看| 99久久精品热视频| av在线老鸭窝| 有码 亚洲区| 麻豆一二三区av精品| 丰满乱子伦码专区| 欧美日韩国产亚洲二区| 嫩草影院新地址| 国内精品一区二区在线观看| 国产精品久久久久久av不卡| 国产av码专区亚洲av| 春色校园在线视频观看| 纵有疾风起免费观看全集完整版 | 免费在线观看成人毛片| 91av网一区二区| 18禁动态无遮挡网站| www日本黄色视频网| 欧美极品一区二区三区四区| 男人舔女人下体高潮全视频| 级片在线观看| 亚洲国产欧洲综合997久久,| 国产伦精品一区二区三区四那| 成人特级av手机在线观看| 亚洲精华国产精华液的使用体验| 在线观看美女被高潮喷水网站| a级毛色黄片| 亚洲乱码一区二区免费版| 韩国高清视频一区二区三区| 少妇丰满av| 99久久人妻综合| 99热网站在线观看| 欧美97在线视频| 亚洲18禁久久av| 成人一区二区视频在线观看| 淫秽高清视频在线观看| 搞女人的毛片| 嫩草影院入口| 最近中文字幕2019免费版| 国产精华一区二区三区| 日本免费一区二区三区高清不卡| 久久精品久久久久久久性| 亚洲精品成人久久久久久| 中文字幕亚洲精品专区| 美女国产视频在线观看| 精品久久久噜噜| АⅤ资源中文在线天堂| 看黄色毛片网站| 偷拍熟女少妇极品色| 亚洲欧美日韩东京热| 我要搜黄色片| h日本视频在线播放| 亚洲va在线va天堂va国产| 寂寞人妻少妇视频99o| 最近2019中文字幕mv第一页| 国产真实伦视频高清在线观看| 欧美3d第一页| 久久精品久久久久久噜噜老黄 | 高清av免费在线| 18+在线观看网站| 天美传媒精品一区二区| 老司机影院毛片| 国产亚洲av嫩草精品影院| 成人毛片60女人毛片免费| 免费黄色在线免费观看| 成人三级黄色视频| 午夜久久久久精精品| 天堂√8在线中文| 人人妻人人看人人澡| 深爱激情五月婷婷| 网址你懂的国产日韩在线| 成人无遮挡网站| 免费观看在线日韩| 久久人人爽人人片av| 欧美日韩精品成人综合77777| 免费看av在线观看网站| 久久久久久九九精品二区国产| 精品少妇黑人巨大在线播放 | 中文在线观看免费www的网站| 久久国产乱子免费精品| 国产中年淑女户外野战色| 国产免费福利视频在线观看| 久久午夜福利片| 欧美激情久久久久久爽电影| 深爱激情五月婷婷| 少妇的逼水好多| 一边亲一边摸免费视频| 国产在视频线在精品| 黄片wwwwww| 国产精品伦人一区二区| 欧美又色又爽又黄视频| 麻豆成人午夜福利视频| 国产免费福利视频在线观看| 国产淫片久久久久久久久| av在线蜜桃| 精品一区二区三区视频在线| 亚洲自拍偷在线| 别揉我奶头 嗯啊视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精华一区二区三区| 中文字幕av成人在线电影| 91精品国产九色| 我的女老师完整版在线观看| 亚洲人成网站在线观看播放| 一二三四中文在线观看免费高清| 国产成人精品久久久久久| 少妇人妻精品综合一区二区| 少妇被粗大猛烈的视频| 精品久久久噜噜| 欧美一区二区亚洲| 精品国产一区二区三区久久久樱花 | 丝袜喷水一区| 男女国产视频网站| 国产精品.久久久| 亚洲最大成人手机在线| 亚洲乱码一区二区免费版| 高清在线视频一区二区三区 | 亚洲国产精品国产精品| 国产黄a三级三级三级人| 熟妇人妻久久中文字幕3abv| 午夜福利视频1000在线观看| 日韩av在线免费看完整版不卡| 欧美成人精品欧美一级黄| 国产精品国产三级国产av玫瑰| 波多野结衣巨乳人妻| 亚洲aⅴ乱码一区二区在线播放| 国产午夜福利久久久久久| 麻豆乱淫一区二区| 婷婷色综合大香蕉| www.av在线官网国产| 青春草亚洲视频在线观看| 波多野结衣巨乳人妻| 久久久久久国产a免费观看| 久久久精品大字幕| 亚州av有码| 国产日韩欧美在线精品| 国产一区亚洲一区在线观看| 国产成人精品一,二区| 国产私拍福利视频在线观看| 国产av一区在线观看免费| 成人国产麻豆网| 内地一区二区视频在线| 欧美bdsm另类| 欧美成人免费av一区二区三区| av免费观看日本| 午夜福利成人在线免费观看| 久久精品国产自在天天线| 99热这里只有精品一区| 欧美日韩精品成人综合77777| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 亚洲精品456在线播放app| 美女大奶头视频| 欧美成人精品欧美一级黄| www日本黄色视频网| 一卡2卡三卡四卡精品乱码亚洲| 26uuu在线亚洲综合色| 国产91av在线免费观看| 男女下面进入的视频免费午夜| 精品久久久久久成人av| 丰满人妻一区二区三区视频av| 97超碰精品成人国产| 成人无遮挡网站| 汤姆久久久久久久影院中文字幕 | 国产又黄又爽又无遮挡在线| 久久综合国产亚洲精品| 免费观看在线日韩| 久久草成人影院| 久久精品影院6| www.av在线官网国产| 亚洲欧美中文字幕日韩二区| 久久精品久久久久久噜噜老黄 | 久久精品久久精品一区二区三区| 国产一级毛片七仙女欲春2| 男人狂女人下面高潮的视频| 国产一区亚洲一区在线观看| 国产探花极品一区二区| 国产老妇伦熟女老妇高清| 日韩中字成人| 成人欧美大片| 六月丁香七月| 蜜臀久久99精品久久宅男| 亚洲精品一区蜜桃| 成人漫画全彩无遮挡| 日本黄色视频三级网站网址| 一级av片app| 九色成人免费人妻av| 女人久久www免费人成看片 | 色噜噜av男人的天堂激情| 久久这里有精品视频免费| 亚洲五月天丁香| 午夜日本视频在线| 丝袜喷水一区| 日本免费a在线| 亚洲av日韩在线播放| 亚洲国产成人一精品久久久| 久久精品综合一区二区三区| 成人亚洲欧美一区二区av| 青春草国产在线视频| 国产精品国产三级国产av玫瑰| 久久精品国产亚洲av涩爱| av女优亚洲男人天堂| 国产伦精品一区二区三区视频9| 亚洲成av人片在线播放无| 亚洲最大成人手机在线| 看非洲黑人一级黄片| 小蜜桃在线观看免费完整版高清| 一本久久精品| 日本午夜av视频| 亚洲av中文字字幕乱码综合| 丰满少妇做爰视频| 国产一区二区三区av在线| 成人一区二区视频在线观看| 六月丁香七月| 中文欧美无线码| 精品国内亚洲2022精品成人| 好男人视频免费观看在线| 亚洲欧美精品专区久久| 联通29元200g的流量卡| 97热精品久久久久久| 国产真实伦视频高清在线观看| 亚洲丝袜综合中文字幕| 日韩欧美精品免费久久| 亚洲国产精品sss在线观看| 国产激情偷乱视频一区二区| 一级毛片我不卡| 黄片wwwwww| 美女大奶头视频| 中文字幕av在线有码专区| 国产亚洲av嫩草精品影院| 最新中文字幕久久久久| 中国美白少妇内射xxxbb| 国产老妇伦熟女老妇高清| 99久久精品热视频| 亚洲在久久综合| 51国产日韩欧美|