999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于自適應粒子群算法的燃煤-捕碳機組熱力系統優化設計

2016-10-12 07:37:34付文鋒侯艷峰王藍婧李嘉華楊勇平
動力工程學報 2016年9期
關鍵詞:優化系統設計

付文鋒, 侯艷峰, 王藍婧, 李嘉華, 楊勇平

(1.華北電力大學 電站設備狀態監測與控制教育部重點實驗室,河北保定 071003;2.華北電力大學 控制與計算機工程學院,河北保定 071003)

?

基于自適應粒子群算法的燃煤-捕碳機組熱力系統優化設計

付文鋒1,侯艷峰1,王藍婧2,李嘉華1,楊勇平1

(1.華北電力大學 電站設備狀態監測與控制教育部重點實驗室,河北保定 071003;2.華北電力大學 控制與計算機工程學院,河北保定 071003)

以某1 000 MW超超臨界燃煤機組為例,提出了一種新型燃煤-捕碳機組熱力系統設計方案,建立了該方案下機組的熱經濟性計算框架及回熱系統參數優化模型,并引入自適應粒子群算法進行優化計算.結果表明:新設計燃煤-捕碳機組的熱經濟性顯著改善,循環熱效率比捕碳改造機組相對提高10.7%;自適應粒子群算法收斂快、穩定性好,其優化結果明顯優于其他方法,能夠適用于燃煤-捕碳機組的熱力系統優化設計.

燃煤-捕碳機組; 熱力系統; 自適應粒子群算法; 優化設計

全球變暖是人類共同面對的環境問題.政府間氣候變化專門委員會(IPCC)在《第五次評估報告》中明確指出,化石燃料燃燒和工業生產中產生的CO2是全球溫室氣體增長的主要來源[1].我國從2007年開始成為CO2排放第一大國,面臨著巨大的減排壓力[2].其中,燃煤電廠作為CO2排放大戶,其CO2排放量約占總量的50%,必須承擔起減排和緩解氣候變化的責任.《國家能源科技十二五規劃》中明確提出,到2020年中國的碳排放強度要在2005年的基礎上下降40%~45%[3].因此,開展針對燃煤機組CO2減排技術的研究具有重要意義.

適用于燃煤機組捕碳的方法主要有:燃燒前捕集[4]、富氧燃燒[5]和燃燒后捕集[6].其中,采用前2種技術時,燃燒系統需要進行大幅調整,投資成本較高,在現役機組中應用較少.燃燒后捕集技術適應性強,發展相對成熟,可以廣泛應用于現役機組和新建機組.燃燒后捕集CO2的方法主要有化學吸收法[7]、物理吸附法[8]、低溫蒸餾法[9]和膜分離法等[10].其中,化學吸收法具有捕碳效率高、吸收劑可循環使用以及適合大規模捕碳等優點,吸引眾多學者開展了廣泛研究.在吸收劑方面,Léonard等[8]研究了單乙醇胺(MEA)溶液與CO2的反應機理及其吸收特性;馬雙忱等[9]研究了氨水溶液與MEA溶液捕碳效率的影響因素.在捕碳系統耗能方面,Fu等[10]以仿真方法研究了MEA溶液再生能耗數值;Abu-Zahra等[11-12]研究了捕碳系統操作參數對系統能耗的影響.在系統供能方面,劉芳等[13]和Mokhtar等[14]分別提出了以機組抽汽和太陽能為捕碳系統供能的方案.在機組能耗變化方面,Freguia等[15]和Wagener 等[16]分別研究了吸收劑采用氨水溶液和MEA溶液時,捕碳系統對機組能耗的影響.

以現有的研究報道來看,捕碳系統與熱力系統的集成通常僅是簡單地將捕碳裝置加裝在燃煤機組的某些特定位置,而原熱力系統并未隨捕碳裝置的引入同時進行優化.由于新增的捕碳系統破壞了原熱力系統結構,導致系統內部各參數偏離最優,影響了機組的熱經濟性,因此應當對燃煤-捕碳機組的整體熱力系統重新進行優化設計.筆者基于吸收劑采用MEA溶液、再生能源取自機組抽汽的捕碳方法,提出了一種新型燃煤-捕碳機組耦合熱力系統,并引入自適應粒子群(AWPSO)算法[17]對其熱力系統參數進行優化設計,同時對不同設計方法進行了對比分析.

1 捕碳系統工藝流程

捕碳系統見圖1,主要包括吸收塔、再生塔、貧富液換熱器、再沸器、冷凝器和汽水分離器.由脫硫脫硝系統排出的煙氣經過增壓風機從吸收塔底部自下而上流動,與此同時,用于吸收CO2的MEA貧液從吸收塔頂部自上而下噴淋.煙氣與MEA貧液在吸收塔中逆流接觸并發生化學反應,絕大部分CO2被吸收,貧液變為富液,被處理過的煙氣經煙囪排向大氣.富液經過貧富液換熱器后從再生塔頂部自上而下噴淋入塔,之后進入再沸器進行解析,富液轉變為貧液.解析出的CO2經再生塔進入塔頂冷凝器和汽水分離器進行提純,其中冷凝水回流到再生塔中,高純度的CO2經過壓縮封存以備后期利用.再生得到的貧液由再生塔底流出,經貧富液換熱器和冷凝器進行熱量交換后與補充的MEA溶液混合再次進入吸收塔;MEA溶液在貧、富液泵的推動下,往復循環于吸收塔與再生塔之間,構成了連續吸收和解析CO2的工藝流程.

圖1 煙氣CO2捕集系統Fig.1 Flow chart of the flue gas CO2 capture system

MEA溶液具有較強的堿性,溫度在20~50 ℃時可與煙氣中的CO2迅速反應生成較穩定的氨基甲酸鹽,使煙氣中的CO2得以脫除.氨基甲酸鹽在再沸器中被加熱,溫度高于105 ℃時發生分解,釋放出CO2,溶劑MEA得以再生[18].解析過程需要大量耗能,該部分能量由再沸汽(一般為飽和水蒸氣或微過熱水蒸氣)提供.當再沸器內溫度高于122 ℃時,再生得到的MEA容易發生熱降解,因此需要選擇合適的再沸汽源并嚴格控制再沸汽參數.捕碳系統主要操作參數見表1.

2 新型燃煤-捕碳機組熱力系統結構設計

由于MEA再生能耗巨大,為了滿足這部分熱量要求,需要從機組抽取大量的蒸汽.經過初步估算,當捕碳效率為85%時,捕碳系統耗汽量可達到低壓缸入口蒸汽量的一半左右[19-20].在汽輪機中若抽取如此大量的蒸汽,最適宜的抽汽點應選在中、低壓缸的連接管道上.抽汽在再沸器中放熱后凝結成水,為了回收工質和熱量,此部分水應返回凝結水系統.同時,塔頂冷凝器具有較大的放熱量,可以用來加熱凝結水.經過Aspen Plus軟件模擬,當捕碳效率達到65%以上時,塔頂冷凝器的放熱量足以取代最后一臺低壓加熱器來加熱凝結水.

表1 捕碳系統主要操作參數

基于以上分析,筆者提出了如圖2所示的燃煤-捕碳機組的熱力系統結構,該系統具有如下特點:(1)回熱系統采用“3高3低1除氧”的7級加熱系統.其中,高壓缸設置2段抽汽,中壓缸設置3段抽汽,低壓缸設置2段抽汽;(2)捕碳系統的再沸汽源取自中壓缸排汽,并使中壓缸排汽壓力在額定工況時恰好匹配再沸汽要求;(3)充分利用捕碳系統余熱,其中塔頂冷凝器余熱用于預熱凝結水,再沸器余熱回收于除氧器.

圖2 新型燃煤-捕碳機組熱力系統結構Fig.2 Thermodynamic system of the newly-built coal-fired unit with carbon capture

3 新型燃煤-捕碳機組熱力系統參數設計

3.1自適應粒子群算法

自適應粒子群算法在標準粒子群(PSO)算法的基礎上對慣性權重系數進行了改進,使該算法在收斂速度和穩定性上均得到了提升,其數學描述為

(2)

ω由下式定義:

(3)

式中:ωmax、ωmin分別為權重ω的最大值和最小值;f為粒子當前的適應度值;favg、fmin分別為當前所有粒子的平均適應度值和最小適應度值.

AWPSO算法流程見圖3.

圖3 自適應粒子群算法流程圖Fig.3 Flow chart of the AWPSO algorithm

3.2熱經濟性計算及尋優模型

根據汽輪機汽態膨脹線,加熱器各級抽汽焓hj可以表示為各級抽汽壓力pj的函數:

(4)

式中:j=1,2,…,z,z為加熱器級數.

加熱器的給水焓升τ、蒸汽放熱量q和疏水放熱量γ均可表示為抽汽壓力pj的函數:

(5)

根據各級加熱器熱量平衡和質量平衡建立燃煤-捕碳機組熱力系統汽水分布矩陣方程:

(6)

式中:A為熱力系統結構系數矩陣;α為各級抽汽流量份額矩陣;qf為主系統外的能量流矩陣,包括捕碳系統能耗和給水泵焓升等;τ為各加熱器給水焓升矩陣.

通過求解式(6),各級抽汽流量份額αj也可以表示為抽汽壓力pj的函數:

(7)

機組循環熱效率η可按下式計算:

(8)

式中:w為單位工質的循環做功量;q0為單位工質的循環吸熱量;h0為新蒸汽焓;hc為汽輪機排汽焓;hfw為鍋爐給水焓;qrh為工質再熱吸熱量;n為再熱前的抽汽流數.

結合式(4)和式(7),循環熱效率最終可以表示為抽汽壓力pj的函數:

(9)

以機組循環熱效率為尋優目標并考慮約束條件,最終可以確定燃煤-捕碳機組熱力系統優化數學模型為

(10)

約束條件:pc

式中:p0和pc分別為汽輪機的主蒸汽壓力和排汽壓力.

3.3應用實例

以某1 000 MW超超臨界燃煤機組作為參考機組,機組額定參數和回熱系統參數見表2和表3.

表2 某1 000 MW燃煤機組額定參數

參考上述機組確定新設計燃煤-捕碳機組的邊界參數.其中,主蒸汽壓力、主蒸汽溫度、再熱蒸汽壓力、再熱蒸汽溫度和排汽壓力保持不變.新設計的回熱抽汽個數變為7,由塔頂冷凝器首先加熱凝結水并取消了8號加熱器.由于1號加熱器和2號加熱器的抽汽壓力會影響到鍋爐給水溫度和再熱蒸汽壓力,因此在優化計算時未改變1號和2號加熱器的抽汽壓力.

表3 某1 000 MW燃煤機組回熱系統參數

應用AWPSO算法,構造適應度函數F=-η.參數設置如下:種群規模N=20;學習因子c1=c2=2;慣性權重wmax=0.9,wmin=0.5;迭代步數t=500.經過50次重復優化計算,得到了各級抽汽壓力的設計參數.在相同的熱力系統結構下,分別采用平均分配法、遺傳(GA)算法和PSO算法對熱力系統參數進行設計.不同方法的比較結果見表4和表5.3種智能尋優方法對應的適應度值進化曲線見圖4.

表4不同方法的優化計算結果

Tab.4 Calculation results of different optimization methods MPa

4 對比分析

4.1傳統設計方法與智能尋優方法的對比

由表5可見,在對燃煤-捕碳機組回熱系統進行設計時,智能尋優方法明顯優于以平均分配法為代表的傳統設計方法.平均分配法由于忽略了蒸汽放熱量隨加熱器給水焓升的變化,與機組實際情況差別較大,導致優化效果較差.而在應用GA算法、PSO算法和AWPSO算法進行優化時,所采用的數學模型充分考慮了各種實際因素的影響,沒有對熱力系統進行簡化和假設,所以取得了良好的優化效果.由于燃煤-捕碳機組熱力系統自身的特殊性,目標函數與優化變量之間的耦合性和非線性更強,智能尋優方法更適用于這種復雜系統的優化設計.

表5不同捕碳方案下的給水焓升分配及循環熱效率的比較

Tab.5Comparison of feed-water enthalpy rise distribution and cycle thermal efficiency under different CO2capture schemes

對比參數改造方案1)新設計燃煤-捕碳機組方案平均分配法AWPSO算法PSO算法GA算法τ1/(kJ·kg-1)175.54179.95179.95179.95179.95τ2/(kJ·kg-1)239.81136.30189.62189.71191.30τ3/(kJ·kg-1)120.64136.30222.27217.74210.81τ4/(kJ·kg-1)247.56136.30161.6166.03171.38τ5/(kJ·kg-1)94.01136.3082.6783.9482.19τ6/(kJ·kg-1)80.99136.3080.1580.2680.91τ7/(kJ·kg-1)83.01136.3081.3980.0081.11τ8/(kJ·kg-1)75.72凝結水預熱焓升τ/(kJ·kg-1)119.46119.46119.46119.46循環熱效率η0.36920.40070.40870.40860.4066

注:1)以文獻[19]的方案計算.

圖4 最佳適應度值進化曲線Fig.4 Best fitness trends of various optimization algorithms

4.2智能尋優方法的對比

由圖4可見,GA算法的收斂速度最慢,尋優效果最差,而AWPSO算法收斂速度最快,尋優效果最好,這主要是算法自身尋優機制的不同所致.

GA算法在尋優過程中需要通過“選擇”、“交叉”和“變異”來實現“個體”與“種群”的更新,計算復雜度高,在求解到一定范圍后,進化方向偏離最優方向,做了大量無為而冗余的迭代;尋找最優值的過程中比較均勻地向最優區域移動,并且GA算法過早收斂于局部最優,使得尋優效果較差.AWPSO算法與PSO算法通過隨機速度來改變“個體”,計算復雜度相對較低,“個體”與“種群”在更新時具有高度的隨機性,粒子具有“記憶”特性,同時關注自身最優值與種群最優值,整個尋優過程中始終跟隨當前最優解,從而能夠在較短的時間內找到最優解.AWPSO算法又在PSO算法的基礎上對慣性權重描述進行了改進,增強了算法前期的收斂能力和后期的搜索能力,因而具有更高的求解速度和精度.

4.3設計方案的主要優點

當主蒸汽質量流量相同時,參考機組、捕碳改造機組和新設計燃煤-捕碳機組的發電功率和循環熱效率的對比見圖5.

圖5 機組發電功率及循環熱效率的對比Fig.5 Comparison of unit power output and cycle thermal efficiency

由圖5可以看出,若在參考機組的基礎上直接進行捕碳改造,并未對系統參數進行優化,發電功率和循環熱效率分別下降了21.8%和25%;對燃煤-捕碳熱力系統進行重新優化設計后,發電功率和循環熱效率分別比改造機組提高了7.8%和10.7%.這是由于新設計的燃煤-捕碳熱力系統有以下幾方面優點:(1)再沸汽源參數更接近捕碳系統的要求,減少了參數調節過程中的能量損失;(2)回收了塔頂冷凝器的余熱,代替了最后一級低壓加熱器;(3)對各級加熱器的給水焓升重新進行了最優化分配.以上各因素使得新設計燃煤-捕碳機組的熱經濟性更具優勢.經計算,若機組年利用小時數以5 000 h計,新設計的燃煤-捕碳機組比改造機組年多發電約3.05×108kW·h,年減排CO2約3.265×106t,具有可觀的經濟效益和社會效益.

5 結 論

(1) 重新設計了一種燃煤-捕碳機組熱力系統,可以有效地減少再沸汽源的能量損失并能利用塔頂冷凝器的余熱.

(2) 建立了熱力系統參數優化模型,引入自適應粒子群算法對某1 000 MW超超臨界燃煤-捕碳機組進行尋優計算,結果表明:新系統在優化后,機組的循環熱效率比未優化的捕碳改造機組提高了10.7%.

(3) 自適應粒子群算法在應用于燃煤-捕碳熱力系統優化時,收斂快、精度高,適用于捕碳系統與燃煤機組的集成設計,并可以推廣應用于其他復雜能量系統的參數優化問題.

[1]王偉光,鄭國光.氣候變化綠皮書:應對氣候變化(2014)[M].北京:社會科學文獻出版社,2014:96-124.

[2]徐鋼,田龍虎,劉彤,等.中國電力工業CO2減排戰略分析[J].中國電機工程學報,2011,31(17):1-8.

XU Gang, TIAN Longhu, LIU Tong,etal. Strategic analysis of CO2mitigation in Chinese power industry[J]. Proceedings of the CSEE, 2011, 31(17): 1-8.

[3]楊勇平,楊志平,徐鋼,等.中國火力發電能耗狀況及展望[J].中國電機工程學報,2013,33(23):1-11.

YANG Yongping, YANG Zhiping, XU Gang,etal. Situation and prospect of energy consumption for China's thermal power generation [J].Proceedings of the CSEE, 2013,33(23):1-11.

[4]陳新明,史紹平,閆姝,等.燃燒前CO2捕集技術在IGCC發電中的應用[J].化工學報,2014,65(8):3193-3201.

CHEN Xinming, SHI Shaoping, YAN Shu,etal. Application of CO2capture technology before burning in IGCC power generation system [J].Journal of Chemical Industry and Engineering(China), 2014,65(8):3193-3201.

[5]閻維平,董靜蘭,任海鋒,等.增壓富氧燃燒與捕集CO2電站的經濟性分析[J].動力工程學報,2012,32(9):712-717.

YAN Weiping, DONG Jinglan, REN Haifeng,etal. Economic analysis of the power station with pressurized oxy-fuel combustion and CO2capture [J].Journal of Chinese Society of Power Engineering, 2012,32(9):712-717.

[6]羅聰,鄭瑛,丁寧,等.納米復合鈣基高溫CO2吸收劑的合成與性能[J].中國電機工程學報,2011,31(8):45-50.

LUO Cong, ZHENG Ying, DING Ning,etal.Synthesis and performance of a nano synthetic Ca-based sorbent for high temperature CO2capture[J].Proceedings of the CSEE, 2011,31(8):45-50.

[7]劉練波,郜時旺,許世森,等.燃煤煙氣CO2捕集系統與電廠系統集成分析[J].中國電機工程學報,2014,34(23):3843-3848.

LIU Lianbo,GAO Shiwang,XU Shisen,etal.Analysis on system integration between carbon capture and coal-fired power plant [J].Proceedings of the CSEE,2014,34(23):3843-3848.

[8]LéONARD G, TOYE D, HEYEN G. Relevance of accelerated conditions for the study of monoethanolamine degradation in post-combustion CO2capture[J].Canadian Journal of Chemical Engineering, 2015, 93(2):348-355.

[9]馬雙忱,王夢璇,孫云雪,等.氨水與MEA的CO2捕集對比研究[J].動力工程學報,2012,32(1):52-58.

MA Shuangchen, WANG Mengxuan, SUN Yunxue,etal. Efficiency comparison of CO2removal respectively with ammonia solution and MEA [J].Journal of Chinese Society of Power Engineering, 2012,32(1):52-58.

[10]FU K Y, RONGWONG W, LIANG Z W,etal. Experimental analyses of mass transfer and heat transfer of post-combustion CO2absorption using hybrid solvent MEA-MeOH in an absorber[J]. Chemical Engineering Journal, 2015,260:11-19.

[11]ABU-ZAHRA M R M, SCHNEIDERS L H J, NIEDERER J P M,etal. CO2capture from power plants: Part I. a parametric study of the technical performance based on monoethanolamine[J]. International Journal of Greenhouse Gas Control, 2007, 1(1):37-46.

[12]ABU-ZAHRA M R M, NIEDERER J P M, FERON P H M,etal. CO2capture from power plants: Part II. a parametric study of the economical performance based on mono-ethanolamine[J]. International Journal of Greenhouse Gas Control, 2007, 1(2):135-142.

[13]劉芳,王淑娟,陳昌和,等.電廠煙氣氨法脫碳技術研究進展[J].化工學報,2009,60(2):269-278.

LIU Fang, WANG Shujuan, CHEN Changhe,etal. Research progress of CO2capture by using ammonia from flue gas of power plant [J]. Journal of Chemical Industry and Engineering(China), 2009, 60 (2): 269-278.

[14]MOKHTAR M, ALI M T, KHALILPOUR R,etal. Solar-assisted post-combustion carbon capture feasibility study[J]. Applied Energy, 2012, 92(4):668-676.

[15]FREGUIA S, ROCHELLE G T. Modeling of CO2capture by aqueous monoethanolamine[J]. Aiche Journal, 2003, 49(7):1676-1686.

[16]WAGENER D H V, LIEBENTHAL U, PLAZA J M,etal. Maximizing coal-fired power plant efficiency with integration of amine-based CO2capture in greenfield and retrofit scenarios[J]. Energy, 2014, 72(7):824-831.

[17]付文鋒,王藍婧,李飛,等.基于改進動態自適應粒子群優化算法的汽輪機組給水回熱分配優化方法[J].動力工程學報,2014,34(8):653-659.

FU Wenfeng, WANG Lanjing, LI Fei,etal. Distribution optimization on feedwater enthalpy rise in steam turbine unit based on modified PSO algorithm [J].Journal of Chinese Society of Power Engineering, 2014,34(8):653-659.

[18]韓中合,王繼選.CO2減排機理及與燃煤電廠集成特性的研究[J].動力工程學報,2013,33(10):808-814.

HAN Zhonghe, WANG Jixuan. Study on mechanism of CO2capture process and its integration with power generation systems [J]. Journal of Chinese Society of Power Engineering, 2013, 33(10):808-814.

[19]劉文毅,唐寶強,徐鋼,等.大規模碳捕集電站貧液CO2負載率優化研究[J].動力工程學報,2015,35(1):76-82.

LIU Wenyi, TANG Baoqiang, XU Gang,etal. Optimization on lean solvent CO2loading in coal-fired power plants with large-scale CO2capture [J]. Journal of Chinese Society of Power Engineering, 2015, 35 (1): 76-82.

[20]OEXMANN J, HENSEL C, KATHER A. Post-combustion CO2-capture from coal-fired power plants:preliminary evaluation of an integrated chemical absorption process with piperazine-promoted potassium carbonate[J]. International Journal of Greenhouse Gas Control, 2008, 2(4):539-552.

Optimized Design for Thermodynamic System of Coal-fired Power Plants with CO2Capture Based on AWPSO Algorithm

FUWenfeng1,HOUYanfeng1,WANGLanjing2,LIJiahua1,YANGYongping1

(1. MOE's Key Laboratory of Condition Monitoring and Control for Power Plant Equipment,North China Electric Power University, Baoding 071003, Hebei Province, China;2. School of Control and Computer Engineering, North China Electric Power University, Baoding 071003, Hebei Province, China)

Taking the 1 000 MW ultra-supercritical coal-fired unit as an example, a new design scheme was proposed for the thermal system of power plant with CO2capture, based on which a computing framework for the thermal efficiency and a parameter optimization model for the regenerative system were set up, while the adaptive weighted particle swarm optimization (AWPSO) algorithm was applied for relevant calculations. Results show that the thermal economy of the newly-designed coal-fired power unit with CO2capture has been significantly improved, with its cycle thermal efficiency 10.7% higher than the retrofitted unit before optimization. The AWPSO algorithm is characterized by quick convergence and high stability, and its optimization results are obviously better than other methods, which therefore may be used for design optimization on thermodynamic system of coal-fired carbon-captured power units.

coal-fired power plant with CO2capture; thermodynamic system; AWPSO algorithm; design optimization

2015-10-09

國家自然科學基金煤炭聯合基金資助項目(U1261210);中央高校基本科研業務費專項資金資助項目(2014MS109, 2014MS135)

付文鋒(1982-),男,河北唐山人,講師,博士研究生,主要從事節能原理與技術方面的研究.電話(Tel.):13673222032;

E-mail:fwf_1982@163.com.

1674-7607(2016)09-0746-07

TK284.1

A學科分類號:470.20

猜你喜歡
優化系統設計
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
WJ-700無人機系統
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
瞞天過海——仿生設計萌到家
藝術啟蒙(2018年7期)2018-08-23 09:14:18
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
設計秀
海峽姐妹(2017年7期)2017-07-31 19:08:17
主站蜘蛛池模板: 成人自拍视频在线观看| 精品一区二区三区视频免费观看| 国产精品区网红主播在线观看| 亚洲视频三级| 国产美女无遮挡免费视频网站| 欧美a√在线| 国产凹凸视频在线观看| 狠狠干综合| 97青青青国产在线播放| 国产又粗又爽视频| 国产美女视频黄a视频全免费网站| 色香蕉影院| 日本欧美午夜| 亚洲中文字幕23页在线| 国产99精品久久| 青青草国产免费国产| 亚洲无码在线午夜电影| 国产精品网拍在线| 亚洲综合18p| 久一在线视频| 国产无套粉嫩白浆| 欧洲欧美人成免费全部视频| 久久99国产综合精品1| 最新国产在线| 四虎影视国产精品| 欧美亚洲日韩不卡在线在线观看| 国产高清在线精品一区二区三区| 91久久偷偷做嫩草影院免费看 | 国产精品毛片在线直播完整版| 国产欧美亚洲精品第3页在线| 99伊人精品| 黄色网在线| 在线综合亚洲欧美网站| 久久亚洲日本不卡一区二区| 国内精品手机在线观看视频| 免费看久久精品99| 天天干天天色综合网| 日韩a在线观看免费观看| a国产精品| 欧美日本视频在线观看| 欧美日韩第二页| 性视频久久| 免费一级毛片完整版在线看| 99精品国产电影| 乱人伦99久久| 国产高清国内精品福利| 亚洲热线99精品视频| 亚洲男人天堂久久| 成人在线不卡视频| 波多野结衣一区二区三区四区视频| 色综合五月婷婷| 国产精品偷伦视频免费观看国产| 亚洲国产无码有码| 亚洲色图在线观看| 在线中文字幕日韩| 特级精品毛片免费观看| 亚洲综合18p| 国产精品成人不卡在线观看| 亚洲无码精品在线播放| 丰满人妻久久中文字幕| 国产欧美日韩一区二区视频在线| 国产99视频精品免费视频7| 欧美性色综合网| 五月婷婷综合在线视频| 精品综合久久久久久97超人| 狠狠五月天中文字幕| 久久中文字幕2021精品| AV无码一区二区三区四区| 亚洲第一视频区| 亚洲日韩在线满18点击进入| 国产中文一区a级毛片视频| 免费视频在线2021入口| www.youjizz.com久久| 中文字幕 欧美日韩| 亚洲一区黄色| 被公侵犯人妻少妇一区二区三区| 精品1区2区3区| 亚洲无码精彩视频在线观看| 国产精品极品美女自在线网站| 日韩高清在线观看不卡一区二区| 免费日韩在线视频| 国产精品分类视频分类一区|