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

    鍋爐爐膛熱負(fù)荷分布擬合算法研究

    2014-06-07 10:02:40劉福國楊晨旭
    計(jì)量學(xué)報(bào) 2014年5期
    關(guān)鍵詞:輻射熱水冷壁爐膛

    劉福國, 楊晨旭

    (1.國網(wǎng)山東省電力公司電力科學(xué)研究院,山東濟(jì)南 250002;2.華能沾化熱電有限公司,山東濱州 256800)

    鍋爐爐膛熱負(fù)荷分布擬合算法研究

    劉福國1, 楊晨旭2

    (1.國網(wǎng)山東省電力公司電力科學(xué)研究院,山東濟(jì)南 250002;2.華能沾化熱電有限公司,山東濱州 256800)

    鍋爐爐膛輻射熱負(fù)荷測(cè)量及分布規(guī)律的研究可為蒸發(fā)受熱面及節(jié)流孔設(shè)計(jì)提供必需的數(shù)據(jù)。提出了爐膛熱負(fù)荷測(cè)量數(shù)據(jù)的處理方法,將測(cè)得的熱負(fù)荷數(shù)據(jù)轉(zhuǎn)換成相對(duì)于蒸發(fā)受熱面平均熱負(fù)荷的分配不均系數(shù),該系數(shù)是爐膛相對(duì)高度的函數(shù),它的多項(xiàng)式回歸模型不僅要使誤差平方和達(dá)到最小,而且根據(jù)平均熱負(fù)荷的定義,多項(xiàng)式系數(shù)還要滿足給定的等式約束。對(duì)于采用普通燃燒技術(shù)和空氣分級(jí)燃燒技術(shù)的兩類燃煤鍋爐,利用該回歸模型擬合爐膛輻射熱負(fù)荷數(shù)據(jù),得到不同的分布規(guī)律,并與某300 MW鍋爐爐膛熱負(fù)荷測(cè)量數(shù)據(jù)進(jìn)行了對(duì)比;這種以最小二乘法為基礎(chǔ)的改進(jìn)回歸模型既滿足了回歸變量的內(nèi)在聯(lián)系,又對(duì)測(cè)量數(shù)據(jù)進(jìn)行了最佳逼近。

    計(jì)量學(xué);鍋爐;爐膛輻射熱流量;數(shù)據(jù)處理;改進(jìn)的數(shù)據(jù)回歸模型;空氣分級(jí)燃燒

    1 前 言

    爐膛熱負(fù)荷是影響燃煤鍋爐安全和經(jīng)濟(jì)運(yùn)行的關(guān)鍵參數(shù)之一,該參數(shù)通常根據(jù)設(shè)計(jì)經(jīng)驗(yàn)選取,用于確定受熱面大小或進(jìn)行受熱管的管壁溫度計(jì)算[1~3]。在實(shí)際運(yùn)行中,爐膛輻射熱負(fù)荷分布極不均勻,某些局部區(qū)域的熱負(fù)荷有時(shí)會(huì)達(dá)到平均熱負(fù)荷的數(shù)倍以上,高熱負(fù)荷區(qū)是安全運(yùn)行的薄弱環(huán)節(jié)。受熱面輻射熱負(fù)荷測(cè)量及分布規(guī)律研究是改進(jìn)爐膛設(shè)計(jì)、提高運(yùn)行可靠性的基礎(chǔ)工作,爐膛熱負(fù)荷測(cè)量研究已得到廣泛開展[4~7],一些研究成果被提煉成經(jīng)驗(yàn)設(shè)計(jì)曲線,在前蘇聯(lián)頒布的鍋爐熱力計(jì)算標(biāo)準(zhǔn)中以熱負(fù)荷分配不均系數(shù)的形式給出[8]。近年來,鍋爐燃用煤質(zhì)惡化[9],爐膛受熱面結(jié)焦腐蝕和超溫爆管事故增加[10];同時(shí),低NOx燃燒技術(shù)和超臨界壓力直流鍋爐大量投用[11,12],導(dǎo)致爐膛輻射熱負(fù)荷分布發(fā)生變化,因此,采用新燃燒方式的超高參數(shù)直流鍋爐,爐膛熱負(fù)荷分布研究需要進(jìn)一步開展[13~15]。

    在鍋爐熱力計(jì)算或水動(dòng)力計(jì)算中,需要確定爐膛不同位置的熱負(fù)荷,因此,熱負(fù)荷分配不均系數(shù)常擬合成爐膛高度的函數(shù),目前,在進(jìn)行回歸分析時(shí)沒有考慮熱負(fù)荷不均系數(shù)的定義對(duì)模型的要求,這會(huì)造成計(jì)算結(jié)果不合理。本文給出一種新的爐膛熱負(fù)荷測(cè)量數(shù)據(jù)的處理方法,對(duì)于采用普通燃燒技術(shù)和空氣分級(jí)燃燒技術(shù)的兩類鍋爐,利用該方法處理了爐膛輻射熱負(fù)荷數(shù)據(jù),得到不同的分布規(guī)律;本文給出的回歸模型既滿足誤差平方和達(dá)到最小,又兼顧了熱負(fù)荷不均系數(shù)的物理意義,使熱負(fù)荷計(jì)算更為合理。

    2 爐膛熱負(fù)荷測(cè)量數(shù)據(jù)預(yù)處理

    2.1 爐膛輻射熱負(fù)荷及分配不均系數(shù)

    在鍋爐設(shè)計(jì)中,爐膛高度上的輻射熱負(fù)荷q計(jì)算如下[1]

    式中,q為爐內(nèi)某一高度上的輻射熱負(fù)荷,kW/m2;η為爐膛高度上輻射熱負(fù)荷分配不均系數(shù);φ為鍋爐的保熱系數(shù);Bj為實(shí)際燃燒產(chǎn)生煙氣的燃料量,kg/s;Qa為每千克計(jì)算燃料的燃燒產(chǎn)物所擁有的總熱量,kJ/kg;I″l為每千克燃料在爐膛出口生成的煙氣的焓,kJ/kg;∑H為爐膛內(nèi)輻射受熱面總面積,m2。

    式(1)右邊除η之外的項(xiàng)可看作爐膛輻射受熱面的平均熱負(fù)荷q0,即

    大容量電廠鍋爐爐膛內(nèi)輻射受熱面主要是水冷壁蒸發(fā)受熱面,有時(shí)還布置少量的輻射式再熱器或過熱器,式(2)中的∑H為所有輻射受熱面面積之和。

    由式(1)和式(2)可知,η為某一高度上的熱負(fù)荷q與爐膛輻射受熱面平均熱負(fù)荷q0的比值。在爐膛輻射熱負(fù)荷測(cè)量研究時(shí),測(cè)點(diǎn)常布置在水冷壁蒸發(fā)受熱面,計(jì)算η值時(shí),采用水冷壁蒸發(fā)受熱面平均熱負(fù)荷q′0代替爐膛輻射受熱面平均熱負(fù)荷q0,即

    當(dāng)輻射受熱面全部為水冷壁蒸發(fā)受熱面時(shí),爐膛輻射受熱面平均熱負(fù)荷q0即為水冷壁的平均熱負(fù)荷q′0,即q0=q′0。

    式(4)是通過熱負(fù)荷測(cè)量數(shù)據(jù)的回歸分析得到。當(dāng)式(4)已知時(shí),水冷壁蒸發(fā)受熱面平均熱負(fù)荷q′0采用下式計(jì)算

    式中,L為爐膛橫截面的周長,對(duì)于實(shí)際爐膛,L為高度的線段函數(shù);X為爐膛高度。

    由式(3)得q=ηq′0,代入式(5),并結(jié)合式(4)得到

    式(6)積分整理后變?yōu)?/p>

    因此,將水冷壁蒸發(fā)受熱面的熱負(fù)荷數(shù)據(jù)整理成式(4)所示的多項(xiàng)式時(shí),系數(shù)ci不僅要滿足最小二乘原理,還要滿足式(7)的要求。

    2.2 水冷壁蒸發(fā)受熱面平均熱負(fù)荷q′0的計(jì)算

    將水冷壁蒸發(fā)受熱面輻射熱負(fù)荷q數(shù)據(jù)整理成相對(duì)高度x=X/H的多項(xiàng)式

    q′0即為根據(jù)測(cè)量數(shù)據(jù)求得的水冷壁蒸發(fā)受熱面平均輻射熱負(fù)荷。

    2.3 熱負(fù)荷測(cè)量數(shù)據(jù)預(yù)處理

    綜上,輻射熱負(fù)荷測(cè)量數(shù)據(jù)預(yù)處理方法是先將輻射熱負(fù)荷q整理成式(8)所示的多項(xiàng)式函數(shù),然后,根據(jù)式(9)計(jì)算受熱面平均熱負(fù)荷q′0,q′0確定后,利用式(3)計(jì)算熱負(fù)荷分配不均系數(shù)η,進(jìn)而得到熱負(fù)荷分配不均系數(shù)η沿相對(duì)高度x的分布數(shù)據(jù)。

    式(8)代入式(5),積分得到

    3 一種新的熱負(fù)荷分配不均系數(shù)擬合方法

    3.1有約束條件的最小二乘回歸模型

    得到熱負(fù)荷分配不均系數(shù)η數(shù)據(jù)后,將η回歸成式(4)所示的多項(xiàng)式函數(shù),多項(xiàng)式的系數(shù)還要滿足式(7)的要求。

    若有m組分配不均系數(shù)η的數(shù)據(jù)(xi,ηi),i=1,2…,m,根據(jù)最小二乘原理,應(yīng)使誤差的平方和達(dá)到最小,上述回歸模型就是在滿足方程(7)的約束下,求誤差函數(shù)達(dá)到極小值時(shí)的系數(shù)ci的值,即

    3.2 熱負(fù)荷分配不均系數(shù)的實(shí)際數(shù)據(jù)

    表1給出了3組輻射熱負(fù)荷分配不均系數(shù)數(shù)據(jù),其中第一組適用于采用普通燃燒技術(shù)的鍋爐,該組數(shù)據(jù)來自文獻(xiàn)[8]提供的η分布曲線;第2組為針對(duì)采用空氣分級(jí)燃燒技術(shù)的爐膛,該組數(shù)據(jù)是利用設(shè)計(jì)廠家提供的輻射熱負(fù)荷分布曲線轉(zhuǎn)化而來;第3組數(shù)據(jù)是利用一臺(tái)300 MW鍋爐蒸發(fā)受熱面的熱負(fù)荷測(cè)量數(shù)據(jù),采用前文所述的預(yù)處理方法得到。

    表1 熱負(fù)荷分配不均系數(shù)

    3.3 回歸分析方法及結(jié)果分析

    利用前面提出的有約束的最小二乘法回歸模型,分別對(duì)上述3組數(shù)據(jù)進(jìn)行擬合分析,可以得到式(4)中的系數(shù)ci。

    求解是在Matlab 2008a平臺(tái)上進(jìn)行,需要編寫2個(gè)代碼文件,一是誤差函數(shù)文件“myfun.m”,另一個(gè)是約束等式文件“myequ.m”;選擇6次多項(xiàng)式進(jìn)行回歸分析,即式(4)中的n=7;對(duì)于表1中的3組數(shù)據(jù),式(10)中的m分別為m=12,40,11;式(10)表示的有約束條件的函數(shù)最小值問題可能存在多個(gè)解,因此,確定的搜索策略是,先求解在沒有任何約束的情況下,誤差函數(shù)達(dá)到最小時(shí)的系數(shù)ci值,這是純粹的最小二乘法回歸問題,采用函數(shù)“fminsearch”進(jìn)行搜索,語句為

    x′為無約束搜索的初值。這次搜索得到系數(shù)ci的解為[c01,c02,…,c07],該解作為上述有約束搜索的初值,即x0=[c01,c02,…,c07],進(jìn)而用函數(shù)“fmincon”求解有約束的函數(shù)極值問題,語句為

    [x,fval,exitflag,output,lambda]=fmincon(@m(xù)yfun,x0,[],[],[],[],lb,ub,@m(xù)yequ)

    這條語句中的lb、ub為搜索的下限和上限,規(guī)定此次搜索在x0附近進(jìn)行,即

    式中,ε為搜索寬度,可取為ε=0.1~0.5。

    采用有約束的最小二乘回歸模型對(duì)上述3組數(shù)據(jù)進(jìn)行擬合分析,得到式(4)的多項(xiàng)式系數(shù)見表2,表中最后一項(xiàng)為這些系數(shù)代入式(7)左邊的計(jì)算結(jié)果,可見,多項(xiàng)式系數(shù)基本滿足約束等式(7)的要求。

    表2 不同條件下多項(xiàng)式(7)的擬合系數(shù)

    圖1、2為擬合多項(xiàng)式和擬合采用的數(shù)據(jù)的對(duì)比,擬合多項(xiàng)式的系數(shù)滿足式(7)所示的條件約束,且使誤差函數(shù)達(dá)到最小,這是一種改進(jìn)型的最小二乘法。

    圖1 普通燃燒技術(shù)和空氣分級(jí)燃燒技術(shù)的爐膛熱負(fù)荷分布

    由圖1可以看出,與文獻(xiàn)[8]給出的典型的普通燃燒技術(shù)相比,空氣分級(jí)配風(fēng)燃燒技術(shù)爐膛熱負(fù)荷的最高點(diǎn)的位置明顯上移,上移的距離接近1/3的相對(duì)爐膛高度,這是因?yàn)榉旨?jí)配風(fēng)改變了燃燒的燃燒分布。

    圖2給出了某鍋爐在200 MW負(fù)荷下爐膛熱負(fù)荷分布的測(cè)量結(jié)果,以及與典型的普通燃燒技術(shù)爐膛熱負(fù)荷分布的對(duì)比,由于該鍋爐未采用分級(jí)燃燒,爐膛熱負(fù)荷分布較接近文獻(xiàn)[8]給出的數(shù)據(jù)。

    圖2 實(shí)測(cè)鍋爐爐膛熱負(fù)荷分布與文獻(xiàn)[8]的數(shù)據(jù)對(duì)比

    4 結(jié) 語

    (1)當(dāng)燃煤鍋爐爐膛內(nèi)未布置再熱器或過熱器等輻射受熱面時(shí),水冷壁蒸發(fā)受熱面平均熱負(fù)荷等于爐膛輻射受熱面平均熱負(fù)荷。

    (2)將爐膛水冷壁蒸發(fā)受熱面的輻射熱負(fù)荷整理成分配不均系數(shù),然后進(jìn)行分配不均系數(shù)的多項(xiàng)式回歸分析時(shí),多項(xiàng)式系數(shù)不僅要滿足最小二乘原理,還要滿足等式約束;給出一個(gè)有約束條件的最小二乘法回歸模型,該模型也可應(yīng)用于其它類似的工程問題中。

    (3)有約束的最小二乘法回歸方法就是求誤差函數(shù)在約束條件下的極小值問題,給出了在matlab平臺(tái)上的搜索方法。

    (4)由于空氣分級(jí)配風(fēng)改變了燃料的燃燒分布,采用該燃燒技術(shù)的鍋爐爐膛火焰中心明顯上移,與普通燃燒技術(shù)相比,火焰中心上移的距離約為1/3的相對(duì)爐膛高度。

    [1] 馮俊凱,沈幼庭.鍋爐原理及計(jì)算[M].北京:科學(xué)出版社,2003.

    [2] JB/Z 201-83電站鍋爐水動(dòng)力計(jì)算方法[S].

    [3] Mohammad A,Shamshirgaran S R.A case study:the effects of the design factors on the thermal profile of shahid rajaiee boiler[J].AppliedThermalEngineering,2008,28(8):955-961.

    [4] 曹漢鼎,朱智士,丁士發(fā),等.1025t/h亞臨界鍋爐爐膛流動(dòng)、燃燒、傳熱數(shù)學(xué)模型耦合計(jì)算及爐膛熱負(fù)荷分布[J].動(dòng)力工程,1992,12(2):1-8.

    [5] Xu M,Azevedo J L T,Carvalho M G.Modelling of the combustion process and NOxemission in a utility boiler[J].Fuel,2000,70(13):1611-1619.

    [6] 郭宏生,徐通模,惠世恩,等.四角布置切向燃燒鍋爐爐膛內(nèi)的熱負(fù)荷分布[J].動(dòng)力工程,1993,13(4):1-6.

    [7] 劉福國,宗世田.煤焦顆粒燃燒模擬在W型火焰鍋爐爐膛設(shè)計(jì)熱力計(jì)算中的應(yīng)用[J].燃燒科學(xué)與技術(shù),2005,11(5):400-406.

    [8] 北京鍋爐廠.鍋爐機(jī)組熱力計(jì)算標(biāo)準(zhǔn)方法[M].北京:機(jī)械工業(yè)出版社,1976.

    [9] Kouprianov V I.Modeling of thermal characteristics for a furnace of a500 MW boiler fired with high-ash coal[J].Energy,2001,26(9):839-853.

    [10] 蔣翀,張明杰,許家偉,等.600 MW超超臨界機(jī)組鍋爐下水冷壁出口工質(zhì)溫度的分布特性[J].動(dòng)力工程,2009,29(4):335-340.

    [11] Taniguchi M,Yamamoto K,Kobayashi H,etal.A reduced NOxreaction model for pulverized coal combustion under fuel-rich conditions[J].Fuel,2002,81(3):363-371.

    [12] Pan J,Yang D,Yu H,etal.Mathematical modeling and thermal-hydraulic analysis of vertical water wall in an ultra supercritical boiler[J].AppliedThermal Engineering,2009,29(11):2500-2507.

    [13] 張永福.鍋爐參數(shù)對(duì)爐內(nèi)傳熱的影響[J].熱力發(fā)電,2000,29(1):15-16.

    [14] Zhang M,Bie R S,Yu Z Z,etal.Heat flux profile of the furnace wall of a 300 MWe CFB boiler[J].Powder Technology,2010,203(3):548-554.

    [15] Li Z Q,Jing JP,Liu G K,etal.Measurement of gas species,temperatures,char burnout,and wall heat fluxes in a200MWe lignite-fired boiler at different loads[J].AppliedEnergy,2010,87(4):1217-1230.

    Study on Radiant Heat Fluxes Fitting Algorithm in Furnace of Utility Boiler

    LIU Fu-guo1, YANG Chen-xu2
    (1.State Grid Shandong Electric Power Research Institute,Jinan,Shandong 250002,China;
    2.Huaneng Zhanhua Cogeneration Co Ltd,Binzhou,Shandong 256800,China)

    Design of evaporative heating surface and restriction orifices benefit from measurements and investigations of furnace radiantheat fluxes in coal-fired boilers.A newmethod of heat fluxes data processing is proposed,theheat fluxesdata are converted into fluxes factors basing on the average heat flux of the furnace surface,the heat fluxes factor is a function of relative height of furnace,and in it's polynomial regressionmodel,ensuring the sum of squared errors to theminimum is as important asmeeting the equality constraintof polynomial coefficients according to the definition of the average heat flux.An improved version of polynomial regression model based on the least squares method is presented accordingly,engineers constantly have similar needs in other problems.The improved regressionmodelwas used to analyze furnace heat fluxes data for two types of coal-fired boiler with ordinary combustion technology and air-staged combustion technology.This new regressionmodel based on the least squaresmethod fit the data optimally tomeet the intrinsic link between the regression variables.

    Metrology;Boiler;Furnace radiant heat flux;Data processing;Improved data regression model;Air staged combustion

    TB941

    A

    1000-1158(2014)05-0440-05

    10.3969/j.issn.1000-1158.2014.05.07

    2012-08-24;

    2014-06-21

    劉福國(1969-),男,江蘇邳州人,山東電力科學(xué)研究院高級(jí)工程師,主要從事鍋爐性能優(yōu)化研究。lephico@163.com

    猜你喜歡
    輻射熱水冷壁爐膛
    小小觀察家
    超超臨界循環(huán)流化床鍋爐膜式水冷壁管 溫度與應(yīng)力分析
    1000MW超超臨界鍋爐水冷壁開裂原因分析
    二維爐膛氣液兩相對(duì)沖流動(dòng)數(shù)值模擬
    層燃型垃圾焚燒鍋爐的爐膛與爐膛溫度簡(jiǎn)析
    水冷壁高溫腐蝕研究進(jìn)展
    船用燃油輔鍋爐爐膛爆燃分析
    水上消防(2020年2期)2020-07-24 09:27:06
    某廠鍋爐水冷壁管兩次典型爆管失效案例分析及預(yù)防
    肯尼亞項(xiàng)目大型池火輻射熱模型選擇及安全間距確定
    基于CMOS工藝鎢微測(cè)輻射熱計(jì)陣列集成芯片的設(shè)計(jì)與制作*
    另类精品久久| 十八禁网站网址无遮挡| 夜夜爽天天搞| 色综合欧美亚洲国产小说| netflix在线观看网站| 在线观看www视频免费| 两个人免费观看高清视频| 国产伦人伦偷精品视频| 国产区一区二久久| 99九九在线精品视频| 纯流量卡能插随身wifi吗| 欧美黄色片欧美黄色片| 麻豆av在线久日| 可以免费在线观看a视频的电影网站| 国产99久久九九免费精品| 欧美成人午夜精品| 久久av网站| 中亚洲国语对白在线视频| 桃花免费在线播放| 成人影院久久| 黄色毛片三级朝国网站| 天堂8中文在线网| 亚洲情色 制服丝袜| 亚洲欧洲日产国产| 亚洲avbb在线观看| 色综合欧美亚洲国产小说| 在线播放国产精品三级| av天堂久久9| 两个人看的免费小视频| 黄色怎么调成土黄色| 热re99久久国产66热| 午夜福利影视在线免费观看| 麻豆乱淫一区二区| 一二三四社区在线视频社区8| a在线观看视频网站| 免费av中文字幕在线| 欧美激情久久久久久爽电影 | 午夜两性在线视频| 亚洲色图av天堂| 少妇猛男粗大的猛烈进出视频| 青青草视频在线视频观看| 50天的宝宝边吃奶边哭怎么回事| 人人妻人人添人人爽欧美一区卜| 国产aⅴ精品一区二区三区波| 91老司机精品| 啦啦啦在线免费观看视频4| 欧美乱码精品一区二区三区| 在线观看一区二区三区激情| 夜夜爽天天搞| 人妻久久中文字幕网| 丁香六月欧美| 欧美人与性动交α欧美软件| 国产1区2区3区精品| 精品午夜福利视频在线观看一区 | 操美女的视频在线观看| 久久久精品国产亚洲av高清涩受| 99久久国产精品久久久| 99热网站在线观看| 成人国语在线视频| 黄色a级毛片大全视频| 大片免费播放器 马上看| 在线观看免费高清a一片| 老司机福利观看| 亚洲成a人片在线一区二区| 少妇的丰满在线观看| 亚洲成国产人片在线观看| 69精品国产乱码久久久| aaaaa片日本免费| 亚洲欧洲精品一区二区精品久久久| 日日摸夜夜添夜夜添小说| 天天躁日日躁夜夜躁夜夜| 免费高清在线观看日韩| 成人av一区二区三区在线看| 夜夜骑夜夜射夜夜干| a级毛片黄视频| 国产一区二区三区视频了| 午夜精品国产一区二区电影| 丁香六月天网| 国产91精品成人一区二区三区 | 啦啦啦在线免费观看视频4| 激情视频va一区二区三区| 美女高潮喷水抽搐中文字幕| 无人区码免费观看不卡 | cao死你这个sao货| 国产精品久久电影中文字幕 | 亚洲午夜精品一区,二区,三区| 三上悠亚av全集在线观看| 妹子高潮喷水视频| 日本五十路高清| 搡老岳熟女国产| 人妻一区二区av| 黑丝袜美女国产一区| tocl精华| 亚洲精品久久午夜乱码| 水蜜桃什么品种好| 国产欧美日韩一区二区三区在线| 久久久久久久大尺度免费视频| 一区二区三区国产精品乱码| 两个人看的免费小视频| 国产免费视频播放在线视频| 大片免费播放器 马上看| 黄网站色视频无遮挡免费观看| 91成人精品电影| 宅男免费午夜| 亚洲精品粉嫩美女一区| 热99久久久久精品小说推荐| 在线观看一区二区三区激情| 热99久久久久精品小说推荐| 国产91精品成人一区二区三区 | 亚洲精品成人av观看孕妇| 少妇精品久久久久久久| 成年人免费黄色播放视频| 亚洲中文av在线| 男女无遮挡免费网站观看| 国产真人三级小视频在线观看| 午夜福利,免费看| 精品国产一区二区三区久久久樱花| 视频区图区小说| 黑人操中国人逼视频| 国产精品久久久av美女十八| 国产成人影院久久av| 久久热在线av| 水蜜桃什么品种好| 精品亚洲成a人片在线观看| 久久人人97超碰香蕉20202| 国产精品一区二区精品视频观看| 久久久久网色| 手机成人av网站| av免费在线观看网站| 色在线成人网| 岛国毛片在线播放| 亚洲国产精品一区二区三区在线| 欧美日韩亚洲综合一区二区三区_| 90打野战视频偷拍视频| 久久久久久久精品吃奶| 女警被强在线播放| 人妻久久中文字幕网| 国产老妇伦熟女老妇高清| 亚洲精品一二三| 在线观看免费午夜福利视频| www日本在线高清视频| 波多野结衣av一区二区av| 久久中文字幕人妻熟女| 亚洲av国产av综合av卡| 午夜福利在线观看吧| 男女边摸边吃奶| 麻豆国产av国片精品| 亚洲欧洲日产国产| 欧美精品亚洲一区二区| 伊人久久大香线蕉亚洲五| 波多野结衣av一区二区av| 亚洲av电影在线进入| 18禁裸乳无遮挡动漫免费视频| 一边摸一边抽搐一进一出视频| 1024香蕉在线观看| 国产成人一区二区三区免费视频网站| 窝窝影院91人妻| 一边摸一边做爽爽视频免费| 国产精品久久久久久精品古装| 在线观看舔阴道视频| 亚洲午夜理论影院| 国产亚洲av高清不卡| 久久久国产欧美日韩av| 精品久久久久久久毛片微露脸| 日本一区二区免费在线视频| 久久青草综合色| 咕卡用的链子| 一级毛片女人18水好多| 黑人欧美特级aaaaaa片| 99久久99久久久精品蜜桃| 亚洲欧美精品综合一区二区三区| 久久久久久久国产电影| 亚洲人成77777在线视频| 久久久久国产一级毛片高清牌| 久热爱精品视频在线9| 欧美激情极品国产一区二区三区| 免费av中文字幕在线| 亚洲国产成人一精品久久久| 1024视频免费在线观看| 亚洲av成人一区二区三| 国产一卡二卡三卡精品| 最新的欧美精品一区二区| 亚洲avbb在线观看| 天天躁日日躁夜夜躁夜夜| 在线 av 中文字幕| 午夜日韩欧美国产| 美女福利国产在线| 在线观看66精品国产| 国产欧美日韩一区二区三区在线| 国内毛片毛片毛片毛片毛片| 青青草视频在线视频观看| 精品亚洲成国产av| 汤姆久久久久久久影院中文字幕| svipshipincom国产片| 女性被躁到高潮视频| 午夜两性在线视频| 女同久久另类99精品国产91| 悠悠久久av| h视频一区二区三区| 老司机影院毛片| 一区二区av电影网| 大片免费播放器 马上看| 国产精品久久久av美女十八| 性色av乱码一区二区三区2| 午夜精品久久久久久毛片777| 美女扒开内裤让男人捅视频| 国产老妇伦熟女老妇高清| 嫩草影视91久久| 在线观看人妻少妇| 在线观看免费视频网站a站| 国产精品九九99| 超色免费av| 法律面前人人平等表现在哪些方面| 色婷婷久久久亚洲欧美| 国产在线视频一区二区| 国产高清国产精品国产三级| 97在线人人人人妻| 91国产中文字幕| 国产精品久久久久成人av| 久久天堂一区二区三区四区| 香蕉国产在线看| 一级毛片女人18水好多| 亚洲中文字幕日韩| 嫩草影视91久久| 考比视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 一区二区av电影网| 热re99久久国产66热| 丰满饥渴人妻一区二区三| 国产一卡二卡三卡精品| 一区二区av电影网| 老司机靠b影院| 亚洲国产精品一区二区三区在线| 精品国内亚洲2022精品成人 | 亚洲国产看品久久| 久久久精品国产亚洲av高清涩受| 欧美日韩黄片免| 国产福利在线免费观看视频| 午夜福利视频精品| 天天躁狠狠躁夜夜躁狠狠躁| 69av精品久久久久久 | 国产1区2区3区精品| 亚洲国产欧美网| 在线永久观看黄色视频| 狠狠婷婷综合久久久久久88av| 久久午夜综合久久蜜桃| 亚洲精品久久成人aⅴ小说| 女警被强在线播放| 婷婷成人精品国产| 高清黄色对白视频在线免费看| 久久 成人 亚洲| 他把我摸到了高潮在线观看 | 在线观看一区二区三区激情| 日韩一卡2卡3卡4卡2021年| 激情视频va一区二区三区| 满18在线观看网站| 美国免费a级毛片| 国产精品二区激情视频| 亚洲国产中文字幕在线视频| 国产又爽黄色视频| 最新在线观看一区二区三区| 国产精品影院久久| 国产国语露脸激情在线看| 国产日韩欧美视频二区| 男女边摸边吃奶| 国产精品久久久久成人av| 国产在线免费精品| 我的亚洲天堂| 美女高潮喷水抽搐中文字幕| 国产精品一区二区在线不卡| 精品一品国产午夜福利视频| 成人三级做爰电影| 亚洲国产毛片av蜜桃av| 99热网站在线观看| 人人妻,人人澡人人爽秒播| 国产av又大| 国产精品免费视频内射| 交换朋友夫妻互换小说| 国产精品熟女久久久久浪| 色94色欧美一区二区| 最近最新中文字幕大全电影3 | 午夜福利,免费看| 美女扒开内裤让男人捅视频| 日韩视频在线欧美| 久久久久视频综合| 少妇精品久久久久久久| 在线av久久热| 日本五十路高清| 99精品久久久久人妻精品| 黑人巨大精品欧美一区二区mp4| 久久人妻熟女aⅴ| 国产成人精品久久二区二区免费| 亚洲中文av在线| 国产在线精品亚洲第一网站| 免费在线观看完整版高清| 亚洲欧洲日产国产| 国产主播在线观看一区二区| www.精华液| 亚洲欧美色中文字幕在线| 久久久精品国产亚洲av高清涩受| 水蜜桃什么品种好| 国产精品九九99| 狂野欧美激情性xxxx| 亚洲欧美日韩另类电影网站| 久久狼人影院| 欧美激情极品国产一区二区三区| 日本vs欧美在线观看视频| 一本大道久久a久久精品| 国产不卡一卡二| 18在线观看网站| 看免费av毛片| 十八禁高潮呻吟视频| 男女下面插进去视频免费观看| 国产成人啪精品午夜网站| 蜜桃国产av成人99| 欧美精品一区二区大全| 欧美国产精品va在线观看不卡| 亚洲av日韩精品久久久久久密| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看免费日韩欧美大片| 国产精品 欧美亚洲| 黄色 视频免费看| 黑丝袜美女国产一区| 18禁裸乳无遮挡动漫免费视频| 99精品在免费线老司机午夜| 亚洲精品在线美女| 国产精品久久久人人做人人爽| 免费观看a级毛片全部| 亚洲国产欧美网| 精品一区二区三卡| 男女免费视频国产| 国产片内射在线| 悠悠久久av| 18禁美女被吸乳视频| 精品卡一卡二卡四卡免费| 日本黄色日本黄色录像| 日日摸夜夜添夜夜添小说| 捣出白浆h1v1| 国精品久久久久久国模美| 亚洲国产看品久久| 麻豆国产av国片精品| 国产1区2区3区精品| 午夜精品国产一区二区电影| 免费高清在线观看日韩| 丰满人妻熟妇乱又伦精品不卡| 12—13女人毛片做爰片一| 亚洲一码二码三码区别大吗| 免费一级毛片在线播放高清视频 | 热99国产精品久久久久久7| 涩涩av久久男人的天堂| 欧美日韩亚洲综合一区二区三区_| 熟女少妇亚洲综合色aaa.| 999久久久精品免费观看国产| 国产一区二区三区在线臀色熟女 | 少妇被粗大的猛进出69影院| 欧美日韩亚洲国产一区二区在线观看 | svipshipincom国产片| 精品视频人人做人人爽| 国产精品 欧美亚洲| 日韩精品免费视频一区二区三区| 成人18禁在线播放| 欧美午夜高清在线| 国产人伦9x9x在线观看| 深夜精品福利| 久久99热这里只频精品6学生| 亚洲精品美女久久久久99蜜臀| 淫妇啪啪啪对白视频| 热99久久久久精品小说推荐| 成人国产av品久久久| 亚洲成av片中文字幕在线观看| 亚洲五月婷婷丁香| av天堂久久9| 狠狠精品人妻久久久久久综合| 亚洲av日韩精品久久久久久密| 亚洲第一欧美日韩一区二区三区 | 国产又爽黄色视频| 久9热在线精品视频| 精品国产超薄肉色丝袜足j| 交换朋友夫妻互换小说| netflix在线观看网站| 蜜桃在线观看..| 黄色片一级片一级黄色片| 99久久99久久久精品蜜桃| 亚洲伊人色综图| 啪啪无遮挡十八禁网站| 丰满人妻熟妇乱又伦精品不卡| 欧美亚洲日本最大视频资源| 国内毛片毛片毛片毛片毛片| 亚洲国产中文字幕在线视频| 老司机靠b影院| cao死你这个sao货| 黄色视频在线播放观看不卡| 亚洲一区二区三区欧美精品| 超碰97精品在线观看| 欧美成狂野欧美在线观看| 不卡一级毛片| 女警被强在线播放| 午夜福利视频在线观看免费| 熟女少妇亚洲综合色aaa.| 午夜福利乱码中文字幕| 777米奇影视久久| 丝瓜视频免费看黄片| 亚洲av国产av综合av卡| 精品久久久精品久久久| 免费黄频网站在线观看国产| 香蕉丝袜av| 欧美日韩亚洲国产一区二区在线观看 | 90打野战视频偷拍视频| av国产精品久久久久影院| 美女高潮到喷水免费观看| 精品少妇内射三级| 亚洲欧美一区二区三区黑人| 侵犯人妻中文字幕一二三四区| 日韩人妻精品一区2区三区| 高清欧美精品videossex| av有码第一页| 精品第一国产精品| 涩涩av久久男人的天堂| 午夜福利在线观看吧| 99精品在免费线老司机午夜| 丝瓜视频免费看黄片| 亚洲国产欧美一区二区综合| 少妇猛男粗大的猛烈进出视频| 久久亚洲精品不卡| 欧美黄色淫秽网站| 乱人伦中国视频| 国产精品亚洲av一区麻豆| 国产不卡一卡二| 午夜成年电影在线免费观看| 色在线成人网| 精品久久久久久电影网| 久久婷婷成人综合色麻豆| 亚洲av日韩在线播放| 久久ye,这里只有精品| 国产欧美日韩一区二区三区在线| 满18在线观看网站| 精品少妇久久久久久888优播| 丝袜喷水一区| 波多野结衣一区麻豆| 久热爱精品视频在线9| 精品福利永久在线观看| 亚洲精品在线美女| 精品久久久久久电影网| 久久精品熟女亚洲av麻豆精品| 国产精品一区二区在线不卡| 亚洲中文av在线| 国产成人精品久久二区二区免费| 在线天堂中文资源库| 国产aⅴ精品一区二区三区波| 午夜福利欧美成人| 国产成人一区二区三区免费视频网站| 777久久人妻少妇嫩草av网站| 亚洲精品在线美女| 黄色视频,在线免费观看| √禁漫天堂资源中文www| 极品少妇高潮喷水抽搐| 91国产中文字幕| 黄色视频在线播放观看不卡| 国产精品一区二区在线观看99| 成人18禁高潮啪啪吃奶动态图| 国产真人三级小视频在线观看| 欧美精品啪啪一区二区三区| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美精品综合一区二区三区| 日本av免费视频播放| 国产在线一区二区三区精| 自拍欧美九色日韩亚洲蝌蚪91| 久久天堂一区二区三区四区| 久久国产精品人妻蜜桃| 少妇被粗大的猛进出69影院| 亚洲综合色网址| 人人妻,人人澡人人爽秒播| 成人黄色视频免费在线看| 国产精品99久久99久久久不卡| 大香蕉久久网| 国产精品久久电影中文字幕 | 50天的宝宝边吃奶边哭怎么回事| 久久国产精品男人的天堂亚洲| 国产av又大| 变态另类成人亚洲欧美熟女 | 欧美 日韩 精品 国产| 国产成人免费观看mmmm| 两个人免费观看高清视频| 欧美另类亚洲清纯唯美| tube8黄色片| 亚洲国产av影院在线观看| 老汉色∧v一级毛片| 一夜夜www| 日韩免费高清中文字幕av| 精品少妇黑人巨大在线播放| 欧美精品av麻豆av| 国产精品香港三级国产av潘金莲| 国产极品粉嫩免费观看在线| 又紧又爽又黄一区二区| 久久中文字幕人妻熟女| 免费一级毛片在线播放高清视频 | 精品一品国产午夜福利视频| 99热网站在线观看| 91字幕亚洲| 老司机午夜福利在线观看视频 | 亚洲精品国产一区二区精华液| 深夜精品福利| 欧美人与性动交α欧美软件| 久久性视频一级片| 乱人伦中国视频| 亚洲成a人片在线一区二区| 精品亚洲乱码少妇综合久久| 国产极品粉嫩免费观看在线| 国产精品久久久久久精品古装| 国产精品久久久久久人妻精品电影 | 无限看片的www在线观看| 亚洲男人天堂网一区| 黑人猛操日本美女一级片| 亚洲国产毛片av蜜桃av| 亚洲国产av影院在线观看| av线在线观看网站| 免费久久久久久久精品成人欧美视频| 黄色视频不卡| 欧美日韩av久久| 在线观看免费午夜福利视频| 成年动漫av网址| 国产精品一区二区免费欧美| 国产欧美日韩一区二区三| 国产一区二区 视频在线| 亚洲成人国产一区在线观看| 精品久久蜜臀av无| 精品亚洲成a人片在线观看| 久9热在线精品视频| 亚洲欧美激情在线| 一级片免费观看大全| 免费在线观看视频国产中文字幕亚洲| 一个人免费在线观看的高清视频| 99精品在免费线老司机午夜| 欧美日本中文国产一区发布| 99精品在免费线老司机午夜| 久久久久久久精品吃奶| 纯流量卡能插随身wifi吗| 美女扒开内裤让男人捅视频| 性高湖久久久久久久久免费观看| 老熟妇乱子伦视频在线观看| 99久久国产精品久久久| 国产av一区二区精品久久| 大片免费播放器 马上看| 91av网站免费观看| 欧美激情久久久久久爽电影 | 91成年电影在线观看| 99精品在免费线老司机午夜| 99re6热这里在线精品视频| 精品一区二区三区av网在线观看 | a级毛片在线看网站| 黄色视频在线播放观看不卡| tocl精华| 国产精品香港三级国产av潘金莲| 国产免费视频播放在线视频| 1024视频免费在线观看| 男女午夜视频在线观看| 黄色a级毛片大全视频| 午夜福利在线免费观看网站| 成人永久免费在线观看视频 | www日本在线高清视频| 国产伦人伦偷精品视频| 五月天丁香电影| 免费女性裸体啪啪无遮挡网站| 黑人操中国人逼视频| 亚洲熟妇熟女久久| 日日爽夜夜爽网站| 国产在线免费精品| 日日摸夜夜添夜夜添小说| 欧美黄色片欧美黄色片| 久久国产精品人妻蜜桃| 国产一区二区三区综合在线观看| 手机成人av网站| 国产麻豆69| 欧美激情 高清一区二区三区| 午夜福利影视在线免费观看| 久久精品亚洲av国产电影网| 在线亚洲精品国产二区图片欧美| 老司机亚洲免费影院| 久久中文字幕一级| 国产成人av激情在线播放| 国产免费av片在线观看野外av| 丰满人妻熟妇乱又伦精品不卡| 国产不卡av网站在线观看| 淫妇啪啪啪对白视频| 亚洲av电影在线进入| 久久精品亚洲精品国产色婷小说| 黄频高清免费视频| 欧美日韩精品网址| 国产黄频视频在线观看| 久久人人97超碰香蕉20202| 正在播放国产对白刺激| 免费女性裸体啪啪无遮挡网站| 午夜视频精品福利| 国产高清videossex| 性色av乱码一区二区三区2| 午夜精品国产一区二区电影| 中文字幕最新亚洲高清| 日本a在线网址| 伊人久久大香线蕉亚洲五| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人国产一区在线观看| 在线观看免费高清a一片| 精品免费久久久久久久清纯 | 9191精品国产免费久久| 国产男女超爽视频在线观看| 美国免费a级毛片| 国产福利在线免费观看视频| 国产精品免费视频内射| 制服诱惑二区| 91成人精品电影| 男女高潮啪啪啪动态图| 两个人免费观看高清视频| 久久国产精品人妻蜜桃| 人人妻人人澡人人爽人人夜夜| 五月天丁香电影|