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

    邯鋼燒結(jié)增配YD礦粉的優(yōu)化方法比較

    2016-08-25 05:59:28賈文君
    河南冶金 2016年2期
    關(guān)鍵詞:邯鋼堿度礦粉

    賈文君

    (河北鋼鐵集團(tuán)邯鋼公司)

    ?

    邯鋼燒結(jié)增配YD礦粉的優(yōu)化方法比較

    賈文君

    (河北鋼鐵集團(tuán)邯鋼公司)

    通過(guò)運(yùn)用Excel規(guī)劃求解和MATLAB的線性規(guī)劃法、MTKL法、遺傳算法對(duì)邯鋼燒結(jié)增配YD礦粉的燒結(jié)配料問(wèn)題進(jìn)行了優(yōu)化研究,并從最終結(jié)果、求解效率、易用性等方面比較了4種方法的優(yōu)劣。研究結(jié)果表明: MTKL法和遺傳算法都適用于解決燒結(jié)優(yōu)化配料問(wèn)題,但遺傳算法求解效率相對(duì)要高,更適合推廣。

    燒結(jié)優(yōu)化配料Excel規(guī)劃求解linprog函數(shù)蒙特卡洛法遺傳算法

    0 前言

    目前,鋼鐵工業(yè)在不斷發(fā)展,降本增效已成為鋼鐵企業(yè)發(fā)展的主流。隨著世界鐵礦石資源的不斷變化,相應(yīng)地會(huì)對(duì)鐵礦石的種類、配比不斷做出調(diào)整。當(dāng)燒結(jié)用到多種含鐵原料時(shí),以什么樣的配礦方案進(jìn)行設(shè)計(jì)使燒結(jié)礦性能最好,成本最低,對(duì)于鋼鐵工業(yè)的降本增效具有重要的意義[1]。

    國(guó)內(nèi)過(guò)去多應(yīng)用Excel規(guī)劃求解優(yōu)化燒結(jié)配料,后來(lái)發(fā)展為利用MATLAB的線性規(guī)劃法,優(yōu)化的內(nèi)容也僅限于燒結(jié)礦的化學(xué)成分。隨著研究的不斷深入,燒結(jié)礦的物理性能、冶金性能也被加入到了優(yōu)化模型中來(lái)。隨之而來(lái)的問(wèn)題是變量、約束條件數(shù)目增多,優(yōu)化問(wèn)題規(guī)模變大, MATLAB的MTKL法和遺傳算法也被應(yīng)用于燒結(jié)優(yōu)化配料[2]。采用上述四種方法對(duì)邯鋼燒結(jié)增配YD礦粉進(jìn)行了優(yōu)化研究,得到了不同約束條件下的優(yōu)化燒結(jié)配礦方案,通過(guò)對(duì)各方案的比較分析,認(rèn)為MATLAB的MTKL法和遺傳算法比較適合于現(xiàn)在的燒結(jié)優(yōu)化配料。

    1 邯鋼增配YD礦粉的配料方案

    燒結(jié)礦的配料既要滿足高爐生產(chǎn)的需要(如鐵品位、二元或者三元堿度、硫、磷,還有原料庫(kù)存情況、原料品種),還要使生產(chǎn)成本最小。實(shí)現(xiàn)燒結(jié)配料的優(yōu)化是降低燒結(jié)礦成本的重要措施[3]。

    邯鋼某一時(shí)期增配YD礦粉的燒結(jié)杯方案及含鐵料單價(jià)見(jiàn)表1,各原料的化學(xué)成分見(jiàn)表2。邯鋼燒結(jié)配料一般沿用傳統(tǒng)方式根據(jù)經(jīng)驗(yàn)擬定配比方案,只包含含鐵料的配比,生產(chǎn)單位根據(jù)生產(chǎn)部下達(dá)含鐵料的配比按一定的堿度在進(jìn)行配料。表1中方案1~3為生產(chǎn)部門根據(jù)經(jīng)驗(yàn)擬定的增配YD粉的方案,各方案的混勻料成本也在表1中列出。

    表1 增配YD礦粉的燒結(jié)杯試驗(yàn)方案及混勻料成本比較

    注:表中單價(jià)為邯鋼某一時(shí)期的原料價(jià)格。

    表2 各原料的化學(xué)成分  / %

    2 邯鋼增配YD礦粉的配料優(yōu)化方法比較

    2.1利用線性規(guī)劃優(yōu)化燒結(jié)配料[2-3]

    2.1.1利用Excel規(guī)劃求解優(yōu)化燒結(jié)配料

    Excel“規(guī)劃求解”的方法已經(jīng)成功應(yīng)用于燒結(jié)配料中,可求出燒結(jié)生產(chǎn)過(guò)程成本最小的配料方案。利用Excel“規(guī)劃求解”的方法對(duì)邯鋼增配YD礦粉的含鐵料配比方案進(jìn)行優(yōu)化,從而達(dá)到降成本的目的。

    2.1.1.1變量

    以表1中的11種含鐵料作為決策變量,設(shè)變量Xi代表第i種原料在總配料中的比例,i指本期配料使用的原料種類,i取值從1到11。

    2.1.1.2約束條件

    (1)一般約束

    0≤∑Xi≤1,i=1,2,3,…11,i代表第i種原料。

    (2)資源量約束

    考慮邯鋼燒結(jié)實(shí)際生產(chǎn)工藝、資源量限制和循環(huán)利用,設(shè)計(jì)各原料的資源量約束條件見(jiàn)表3,并且忽略各方案燒損差異。

    表3 優(yōu)化YD礦粉資源量限制   / %

    上述兩項(xiàng)約束條件在規(guī)劃求解參數(shù)窗體中設(shè)定,設(shè)定和修改都比較繁瑣。

    (3)化學(xué)成分約束

    燒結(jié)配料追求鐵品位高,硅含量低,堿度適中,爐渣則希望適度提高氧化鎂含量,硫磷則越低越好。由于Excel“規(guī)劃求解”只能解決目標(biāo)函數(shù)和約束條件都是線性方程的問(wèn)題,故其應(yīng)用受到一定限制,在處理燒結(jié)礦堿度約束時(shí),只有先把堿度的計(jì)算公式處理成線性形式,才能應(yīng)用其求解。本例只考慮含鐵料配比優(yōu)化,對(duì)堿度暫不做約束。設(shè)Fei、Sii、Cai、Mgi、Ali、Tii、Si、Pi分別代表第i種原料的鐵、硅、鈣、鎂、硫、磷含量。則目標(biāo)配料的鐵品位、硅、鈣、鎂、硫、磷含量分別可以用下面公式表示,其中i=1,2,3,…,11(注:用Si 、Ca 、Mg 、Al 、Ti分別代表SiO2、CaO、MgO、Al2O3、TiO2含量。)

    Fe=∑Fei·Xi,本例要求Fe≥60

    Si=∑Sii·Xi,本例要求Si ≤5.5

    Ca=∑Cai·X,本例要求Ca ≥1.6

    Mg=∑Mgi·Xi,本例要求Mg ≥0.5

    Al =∑Ali·Xi,本例要求Al ≤2

    Ti=∑Tii·Xi,本例要求Ti ≤0.1

    S=∑Si·Xi,本例要求S ≤0.3

    P=∑Pi·Xi ,本例要求P ≤0.9

    2.1.1.3目標(biāo)

    使目標(biāo)配料在滿足化學(xué)成分的條件下成本最小。

    (1)成本計(jì)算

    設(shè)Pi代表了第i種原料的價(jià)格,i取值從1到11,i指本期配料使用的原料種類。計(jì)算1t配料的成本:

    第i種原料的成本是:Ci=Pi·Xi·1

    總成本是:

    C=∑Ci=∑Pi·Xi,其中i取值為1,2,…,11。

    (2)目標(biāo)函數(shù)

    MinC=∑Ci=∑Pi·Xi,其中i取值為1,2,…,11。

    在Excel上建立規(guī)劃模型,如圖1所示。在可行域內(nèi)進(jìn)行求解,如果無(wú)解,說(shuō)明約束條件太多,沒(méi)有可行域,需要修改約束條件,重新求解,最終求出最優(yōu)解,找出此時(shí)變量的值,即為最佳方案。

    圖1 電子表格模型

    2.1.1.4結(jié)果

    2.1.2應(yīng)用MATLAB的linprog函數(shù)優(yōu)化燒結(jié)配料

    解決燒結(jié)配料優(yōu)化問(wèn)題,MATLAB的線性規(guī)劃法是最普遍被應(yīng)用的方法,運(yùn)用MATLAB提供的linprog函數(shù)編程可以求得最優(yōu)解。

    利用同Excel規(guī)劃求解一樣的約束條件,運(yùn)用MATLAB的linprog函數(shù)編程,進(jìn)行邯鋼增配YD礦粉方案的優(yōu)化配料,相應(yīng)的程序如圖2所示。

    圖2MATLAB程序

    從圖2可以看出,運(yùn)用linprog函數(shù)編程,所有約束條件設(shè)定比較簡(jiǎn)單直接,修改起來(lái)也比較方便。linprog函數(shù)不適于解決堿度等非線性約束問(wèn)題。

    優(yōu)化出的結(jié)果見(jiàn)表1的方案5,方案5和和使用Excel “規(guī)劃求解”優(yōu)化的方案4 結(jié)果一致。

    由此可見(jiàn),運(yùn)用Excel規(guī)劃求解和MATLAB的linprog函數(shù)優(yōu)化燒結(jié)配料,如果約束條件相同,則優(yōu)化結(jié)果相同;Excel規(guī)劃求解參數(shù)設(shè)定和修改約束條件較為繁瑣,利用MATLAB的linprog函數(shù)較為快捷方便;兩種方法求解效率都高,但不適于解決約束條件是非線性函數(shù)的問(wèn)題。線性規(guī)劃方法得到的配礦方案單純地考慮了成本因素,所需價(jià)格低廉的YD粉數(shù)量較高,方案的質(zhì)量指標(biāo)是否滿足要求,需做進(jìn)一步的驗(yàn)證。

    2.1.3兩種優(yōu)化方案和經(jīng)驗(yàn)擬定方案的混合料成本比較

    利用上述兩種線性規(guī)劃方法優(yōu)化出的含鐵料配比方案,需要按一定的堿度進(jìn)行配料計(jì)算,才能應(yīng)用于生產(chǎn)現(xiàn)場(chǎng),本例中按堿度1.95進(jìn)行燒結(jié)配料,配料計(jì)算每個(gè)方案混合料成本,即含鐵料與燃料、熔劑混合后成本,見(jiàn)表4。

    1.指導(dǎo)自學(xué),注重預(yù)習(xí)。預(yù)習(xí)實(shí)質(zhì)上是學(xué)生自學(xué)的開(kāi)始,小學(xué)階段的學(xué)習(xí)主要依賴?yán)蠋煹陌才牛瑢W(xué)生只要完成老師下達(dá)的簡(jiǎn)單“學(xué)習(xí)指令”就行。但初中則要求學(xué)生自覺(jué)主動(dòng)并且有計(jì)劃地學(xué)習(xí)。一般的學(xué)生只是單純完成學(xué)校和老師交給的作業(yè)就覺(jué)得完成了學(xué)習(xí)任務(wù),而優(yōu)秀的學(xué)生基本上都有預(yù)習(xí)課本的學(xué)習(xí)習(xí)慣。從小學(xué)升入初中要求學(xué)生的學(xué)習(xí)態(tài)度實(shí)現(xiàn)從“要我學(xué)”到“我要學(xué)”的轉(zhuǎn)變。

    表4 增配YD礦粉的燒結(jié)杯試驗(yàn)方案及混合料成本比較

    由表4可知,利用Excel “規(guī)劃求解”和MATLAB的linprog函數(shù)優(yōu)化出增配YD礦粉的方案,每噸混合料成本也和混勻料成本一樣比前三個(gè)方案分別降低了18元、14元、10元。

    2.2應(yīng)用MTKL法和遺傳算法優(yōu)化燒結(jié)配料[4-5]

    2.2.1MTKL法和遺傳算法約束條件

    MTKL法和遺傳算法化學(xué)成分及物理性能的約束條件見(jiàn)表5。

    表5 MTKL法和遺傳算法化學(xué)成分及物理性能的約束條件

    為了探索除了原料化學(xué)成分約束外的物理性能的約束,對(duì)邯鋼2012年度燒結(jié)原料配比與轉(zhuǎn)鼓指數(shù)進(jìn)行統(tǒng)計(jì),對(duì)原料配比與原料粒度組成、平均粒徑進(jìn)行統(tǒng)計(jì)及測(cè)定,并應(yīng)用minitab軟件,采用逐步回歸的方法,對(duì)燒結(jié)轉(zhuǎn)鼓與配比關(guān)系進(jìn)行回歸分析,鐵礦石平均粒徑與鐵礦石配比的回歸分析得出如下回歸分析方程:

    Y(轉(zhuǎn)鼓)= 80.3 + 0.111 BD - 0.0785 BC + 0.0622BT + 0.0245 AF+ 0.0943 SHIH(回歸系數(shù):85)

    Y(平均粒徑) = 35.2 +0.152GZH+0.0733NF-0.553MK - 0.662HX- 0.463YD-0.308BD(回歸系數(shù):90)

    另外由于進(jìn)口礦的比例增大及低價(jià)料的應(yīng)用,一些混勻料粒度組成>10mm含量越限(廠控標(biāo)準(zhǔn)),在優(yōu)化計(jì)算中將>10mm百分含量作為一個(gè)約束條件;由于生產(chǎn)現(xiàn)場(chǎng)采用了厚料層燒結(jié),因此原料的制粒效果對(duì)燒結(jié)影響不容忽視,在試驗(yàn)室對(duì)邯鋼2012年度所用51種礦粉進(jìn)行了濕容量的測(cè)定,以濕容量的平均值作為使用低價(jià)礦的一個(gè)約束條件,應(yīng)用MTKL法和遺傳算法的化學(xué)成分及物理性能的約束條件見(jiàn)表5,資源量的約束和線性規(guī)劃的約束相同。

    2.2.2MTKL法和遺傳算法的優(yōu)化計(jì)算

    在完成上述優(yōu)化準(zhǔn)備工作后,應(yīng)用MTKL法和遺傳算法兩種方法分別編程,對(duì)配加YD礦粉進(jìn)行了最優(yōu)化計(jì)算,可以直接優(yōu)化出堿度為1.95,滿足各項(xiàng)約束條件的配比方案,無(wú)需象線性規(guī)劃那樣再經(jīng)過(guò)配料計(jì)算。

    運(yùn)用MTKL法解決最優(yōu)化問(wèn)題,要設(shè)置兩個(gè)初始參數(shù)L和P。L是初始的目標(biāo)函數(shù),即為燒結(jié)礦的配料成本,一般都把這個(gè)值設(shè)置得較大。P為樣本的空間,即為產(chǎn)生的配料方案數(shù),一般來(lái)說(shuō)P越大,目標(biāo)函數(shù)為最小值的概率越大,計(jì)算所花的時(shí)間也越多。當(dāng)設(shè)置L=70000,P=10000時(shí),最優(yōu)方案和控制參數(shù)分別見(jiàn)表6、表7。計(jì)算結(jié)果表明MTKL法程序簡(jiǎn)單可控制性強(qiáng),目標(biāo)函數(shù)和約束可以根據(jù)具體情況選擇,適用于對(duì)目標(biāo)函數(shù)和約束負(fù)責(zé),不宜用線性規(guī)劃解決的問(wèn)題。

    運(yùn)用英國(guó)謝菲爾德大學(xué)開(kāi)發(fā)的遺傳算法工具箱編程求解最優(yōu)配料方案,需要設(shè)置的參數(shù)更多。初始化種群時(shí),設(shè)置種群數(shù)規(guī)模為40,二進(jìn)制編碼,每個(gè)變量染色體的長(zhǎng)度為10,最大進(jìn)化代數(shù)為1000,選擇方式為輪盤選擇,交叉率0.7,變異率0.005,求解得到的最優(yōu)方案和控制參數(shù)見(jiàn)表6、表7。計(jì)算結(jié)果表明遺傳算法計(jì)算精度高,收斂速度快,約束條件表達(dá)式形式靈活,約束函數(shù)既可以是線性的也可以是非線性的,既可以是顯函數(shù)也可以是隱函數(shù)的形式。由于遺傳算法是模擬自然界“優(yōu)勝劣汰”的規(guī)律,所以收斂速度和求解效率比蒙特卡洛法要高,更適合解決燒結(jié)優(yōu)化配料問(wèn)題。

    表6 利用MTKL和遺傳算法優(yōu)化的配比方案   / %

    表7 利用 MTKL和遺傳算法優(yōu)化的控制參數(shù)

    優(yōu)化結(jié)果分析:

    (1)利用MTKL和遺傳算法優(yōu)化燒結(jié)混合料配比,可以直接計(jì)算出符合一定堿度要求的配比,無(wú)需象線性規(guī)劃那樣再經(jīng)過(guò)配料計(jì)算。

    (2)比較表5、表6和表4,MTKL和遺傳算法不只以化學(xué)成分作為約束條件,另外增加了物理性能的約束條件如轉(zhuǎn)鼓、濕容量、平均粒徑、>10 mm含量, 所以和表3中方案4、方案5比成本增加6~7元/t,但比三個(gè)經(jīng)驗(yàn)擬定方案成本低3~12元。

    (3)在滿足化學(xué)和物理性能約束條件下,比較由遺傳算法及MTKL算法得出的成本,差值僅為1元/t,但MTKL編程相對(duì)于遺傳算法較易掌握。

    (4)通過(guò)蒙特卡洛和遺傳算法得出的配料方案在化學(xué)成分、轉(zhuǎn)鼓指數(shù)方面都滿足了生產(chǎn)要求,降低了成本,并且與公司擬定的配料方案比,各礦種配比變化不大,生產(chǎn)中的可實(shí)施性很強(qiáng)。

    3 結(jié)語(yǔ)

    (1)運(yùn)用Excel規(guī)劃求解、MATLAB的線性規(guī)劃、MTKL和遺傳算法對(duì)邯鋼燒結(jié)增配YD礦粉進(jìn)行了優(yōu)化,在保證燒結(jié)礦化學(xué)成分、物理性能等綜合指標(biāo)符合入爐要求的前提下,每噸混合料的的成本降低了3~18元不等,達(dá)到了降成本的目的。

    (2)運(yùn)用Excel規(guī)劃求解和MATLAB的linprog函數(shù)優(yōu)化燒結(jié)配料,如果約束條件相同,則優(yōu)化結(jié)果相同;Excel規(guī)劃求解參數(shù)設(shè)定和修改約束條件較為繁瑣,利用MATLAB的linprog函數(shù)較為快捷方便;兩種方法求解效率都高,但不適于解決約束條件是非線性函數(shù)的問(wèn)題。

    (3)在解決燒結(jié)配礦最優(yōu)化問(wèn)題時(shí),一旦約束條件數(shù)目增多、形式復(fù)雜, MTKL法和遺傳算法將會(huì)比線性規(guī)劃能更有效地解決問(wèn)題。MTKL相對(duì)遺傳算法編程簡(jiǎn)單,遺傳算法求解效率比蒙特卡洛法要高,更適合解決燒結(jié)優(yōu)化配料問(wèn)題。

    (4)線性規(guī)劃方法得到的配礦方案單純地考慮了成本因素,方案的質(zhì)量指標(biāo)是否滿足要求,需做進(jìn)一步的驗(yàn)證;通過(guò)MTKL法和遺傳算法得出的配料方案在化學(xué)成分、轉(zhuǎn)鼓指數(shù)方面都滿足了生產(chǎn)要求,降低了成本,并且與當(dāng)前的配料方案相比,各礦種的配比變化不大,生產(chǎn)中的可實(shí)施性很強(qiáng)。

    [1]袁曉麗. 燒結(jié)優(yōu)化配礦綜合技術(shù)系統(tǒng)的研究[D/OL].長(zhǎng)沙:中南大學(xué),2007:1-10.

    [2]呂學(xué)偉,白晨光,邱貴寶,等.三種優(yōu)化燒結(jié)配料方法的比較[J].燒結(jié)球團(tuán), 2006,31(2):11-15.

    [3]王道群. 線性規(guī)劃在燒結(jié)礦配料中的應(yīng)用[J].冶金叢刊, 2005,157(3):4-7.

    [4]呂學(xué)偉,白晨光. 基于遺傳算法的燒結(jié)配料綜合優(yōu)化研究[J].鋼鐵,2007,42(4):12-15.

    [5]黃玉明,重鋼燒結(jié)配礦優(yōu)化研究[D/OL].重慶,重慶大學(xué),2005:1-67.

    COMPARISON OF OPTIMIZATION METHODS FOR ADDING WITH YD ON HANSTEEL'S SINTERING

    Jia Wenjun

    (Hansteel Group of HBIS)

    In this paper, researches on sintering optimization with methods of excel solver, Linprog function, MTKL, genetic algorithms are carried out and the four ways are compared through final results , solving efficiency, usability. Results show: MTKL and genetic algorithms are effective for solving the sintering burden optimization, but because of its higher efficiency , the genetic algorithm is suitable for popularization.

    sintering optimizing burdeningexcel solverLinprog functionMTKLgenetic algorithms

    2016—2—19

    聯(lián)系人:賈文君,高級(jí)工程師,河北.邯鄲(056015),河北鋼鐵集團(tuán)公司邯鄲分公司技術(shù)中心鐵前研究室;

    猜你喜歡
    邯鋼堿度礦粉
    礦粉摻量及細(xì)度對(duì)水工混凝土性能的影響研究
    邯鋼高爐INBA水渣系統(tǒng)改造實(shí)踐
    山東冶金(2022年4期)2022-09-14 09:00:14
    基于強(qiáng)度特性的超細(xì)礦粉水泥土配合比設(shè)計(jì)研究
    北方交通(2021年9期)2021-09-16 02:11:58
    河鋼邯鋼5項(xiàng)成果榮獲年度冶金科學(xué)技術(shù)獎(jiǎng)
    四川冶金(2019年5期)2019-12-23 09:04:54
    石粉對(duì)水泥-礦粉混凝土性能的影響
    高堿度銅精礦濾布的再生清洗
    金屬礦山(2014年7期)2014-03-20 14:19:58
    礦粉對(duì)瀝青膠漿的性質(zhì)影響及作用機(jī)理
    2013年邯鋼汽車板產(chǎn)量計(jì)劃逾100萬(wàn)噸
    [知識(shí)小貼士]
    金屬世界(2013年6期)2013-02-18 16:30:47
    有機(jī)肥的堿度及其減緩?fù)寥浪峄臋C(jī)制
    男人的好看免费观看在线视频| 悠悠久久av| 91久久精品国产一区二区成人| 成人精品一区二区免费| 国产私拍福利视频在线观看| 久久久国产成人精品二区| 悠悠久久av| 国产亚洲av嫩草精品影院| 伊人久久精品亚洲午夜| 中国国产av一级| 一级毛片aaaaaa免费看小| 国内久久婷婷六月综合欲色啪| 午夜激情福利司机影院| 免费观看人在逋| av福利片在线观看| 99国产精品一区二区蜜桃av| 搞女人的毛片| 18禁在线播放成人免费| 亚洲丝袜综合中文字幕| 深夜精品福利| 一本精品99久久精品77| avwww免费| 精品日产1卡2卡| 国产三级在线视频| 国产中年淑女户外野战色| 国产精品一区二区性色av| 一边摸一边抽搐一进一小说| 免费看光身美女| 在线观看美女被高潮喷水网站| 看黄色毛片网站| 国产男人的电影天堂91| 国产又黄又爽又无遮挡在线| 国产伦精品一区二区三区视频9| 亚洲精品456在线播放app| 人人妻人人澡欧美一区二区| 亚洲成人久久爱视频| 国产免费一级a男人的天堂| 高清毛片免费观看视频网站| 亚洲精品一卡2卡三卡4卡5卡| 两个人视频免费观看高清| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区三区在线臀色熟女| 啦啦啦韩国在线观看视频| 国产男人的电影天堂91| 成人综合一区亚洲| 国产一区亚洲一区在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产高清激情床上av| 极品教师在线视频| 一级a爱片免费观看的视频| 性欧美人与动物交配| 亚洲成人精品中文字幕电影| 久久久国产成人精品二区| 国产真实伦视频高清在线观看| 欧美3d第一页| 97人妻精品一区二区三区麻豆| 亚洲性夜色夜夜综合| 亚洲丝袜综合中文字幕| 欧美性猛交黑人性爽| 亚洲成av人片在线播放无| 国产av麻豆久久久久久久| 国产成人福利小说| 亚洲av中文字字幕乱码综合| 最近视频中文字幕2019在线8| 波多野结衣巨乳人妻| 激情 狠狠 欧美| 欧美日韩国产亚洲二区| 日本a在线网址| 国产白丝娇喘喷水9色精品| 午夜福利在线观看吧| 全区人妻精品视频| 99在线视频只有这里精品首页| 午夜福利高清视频| 国产中年淑女户外野战色| 国产一区二区三区av在线 | 亚洲真实伦在线观看| 亚洲成a人片在线一区二区| 欧美日韩在线观看h| 一卡2卡三卡四卡精品乱码亚洲| 精品国产三级普通话版| 中文在线观看免费www的网站| 亚洲精品色激情综合| 99久久精品一区二区三区| 淫秽高清视频在线观看| 日韩国内少妇激情av| 性插视频无遮挡在线免费观看| 国产精品亚洲一级av第二区| 欧美色视频一区免费| 亚洲av中文av极速乱| 亚洲成av人片在线播放无| 老熟妇乱子伦视频在线观看| 亚洲丝袜综合中文字幕| 如何舔出高潮| 午夜精品在线福利| 国产黄色小视频在线观看| 欧美一区二区国产精品久久精品| 99久久九九国产精品国产免费| 1024手机看黄色片| 久久国产乱子免费精品| av在线老鸭窝| av在线天堂中文字幕| 在线观看66精品国产| 国产黄色视频一区二区在线观看 | 色哟哟·www| 亚洲真实伦在线观看| 国产精品av视频在线免费观看| 在线观看av片永久免费下载| 婷婷亚洲欧美| 成人精品一区二区免费| 天堂影院成人在线观看| 啦啦啦啦在线视频资源| 一级黄片播放器| 久久中文看片网| 欧美日韩国产亚洲二区| 午夜免费激情av| 日韩av不卡免费在线播放| 亚洲图色成人| 国产一区二区三区在线臀色熟女| 变态另类成人亚洲欧美熟女| 秋霞在线观看毛片| 亚洲综合色惰| 中文资源天堂在线| 插阴视频在线观看视频| 男人舔奶头视频| 国产91av在线免费观看| 色综合站精品国产| 三级男女做爰猛烈吃奶摸视频| 一进一出抽搐动态| 丰满的人妻完整版| 日日摸夜夜添夜夜添av毛片| 国产成人影院久久av| 欧美日韩综合久久久久久| 能在线免费观看的黄片| 久久精品国产亚洲av香蕉五月| 精品久久久久久久人妻蜜臀av| 99久久无色码亚洲精品果冻| 欧美+亚洲+日韩+国产| 级片在线观看| 久久综合国产亚洲精品| 亚洲第一区二区三区不卡| 国产精品伦人一区二区| 国产伦精品一区二区三区视频9| 中文字幕av在线有码专区| 国模一区二区三区四区视频| 午夜福利18| 婷婷六月久久综合丁香| 久久久久久大精品| 国产真实乱freesex| 少妇熟女欧美另类| 在线免费观看的www视频| 在线国产一区二区在线| 午夜福利视频1000在线观看| 亚洲欧美日韩无卡精品| 啦啦啦韩国在线观看视频| 一边摸一边抽搐一进一小说| 成人三级黄色视频| 赤兔流量卡办理| 久久九九热精品免费| 搡女人真爽免费视频火全软件 | 国产成人freesex在线 | 成人高潮视频无遮挡免费网站| 国产麻豆成人av免费视频| 久久久久久久久久黄片| 国产精品一区二区三区四区免费观看 | 久久久久久国产a免费观看| 在线天堂最新版资源| 校园春色视频在线观看| 我要看日韩黄色一级片| 午夜福利高清视频| 亚洲欧美日韩高清在线视频| 一边摸一边抽搐一进一小说| www日本黄色视频网| 久久久久免费精品人妻一区二区| 能在线免费观看的黄片| 五月玫瑰六月丁香| 变态另类成人亚洲欧美熟女| 中文字幕免费在线视频6| 天美传媒精品一区二区| av在线观看视频网站免费| 国产av麻豆久久久久久久| 国产精品嫩草影院av在线观看| 一进一出好大好爽视频| 精品福利观看| 成人亚洲精品av一区二区| 一夜夜www| 久久人人精品亚洲av| 亚洲中文日韩欧美视频| 麻豆成人午夜福利视频| 我的女老师完整版在线观看| 婷婷精品国产亚洲av| 蜜桃久久精品国产亚洲av| а√天堂www在线а√下载| 亚洲精品日韩av片在线观看| 精品国内亚洲2022精品成人| 久久午夜亚洲精品久久| 亚洲av免费在线观看| 国产黄片美女视频| 高清毛片免费看| 六月丁香七月| 最新在线观看一区二区三区| 如何舔出高潮| 欧美成人精品欧美一级黄| 欧美成人精品欧美一级黄| 99久国产av精品| 一卡2卡三卡四卡精品乱码亚洲| 日本一二三区视频观看| 国产乱人视频| 亚洲无线观看免费| 欧美成人精品欧美一级黄| 少妇猛男粗大的猛烈进出视频 | 久久精品国产亚洲av香蕉五月| 小说图片视频综合网站| 亚洲在线观看片| 国产成人aa在线观看| 中文字幕精品亚洲无线码一区| 老司机福利观看| 亚洲欧美精品综合久久99| 国产精品久久久久久久久免| 老司机午夜福利在线观看视频| 蜜桃亚洲精品一区二区三区| 女的被弄到高潮叫床怎么办| av在线老鸭窝| 午夜激情欧美在线| 久久6这里有精品| 日日撸夜夜添| 国产精品综合久久久久久久免费| 午夜日韩欧美国产| 精品午夜福利视频在线观看一区| 国产成人a∨麻豆精品| 97超视频在线观看视频| ponron亚洲| 欧美+日韩+精品| 国内久久婷婷六月综合欲色啪| 日本一二三区视频观看| 国产一区二区亚洲精品在线观看| 亚洲真实伦在线观看| 亚洲av五月六月丁香网| 十八禁国产超污无遮挡网站| 欧美xxxx性猛交bbbb| 露出奶头的视频| 搡老岳熟女国产| 变态另类成人亚洲欧美熟女| 欧美最黄视频在线播放免费| 久久久久性生活片| 亚洲av熟女| 亚洲av美国av| 久久久午夜欧美精品| 欧美激情久久久久久爽电影| 观看免费一级毛片| 国产中年淑女户外野战色| 丝袜美腿在线中文| 日韩三级伦理在线观看| 午夜福利在线在线| 网址你懂的国产日韩在线| 国内精品久久久久精免费| 在线看三级毛片| av卡一久久| 床上黄色一级片| 亚洲成人中文字幕在线播放| 全区人妻精品视频| 国产成人一区二区在线| 国产aⅴ精品一区二区三区波| 国产麻豆成人av免费视频| 丝袜美腿在线中文| 精品久久国产蜜桃| 国产不卡一卡二| 日韩欧美 国产精品| 国产成人freesex在线 | 神马国产精品三级电影在线观看| 少妇高潮的动态图| 日日摸夜夜添夜夜添小说| 精品午夜福利视频在线观看一区| 久久九九热精品免费| 日韩欧美免费精品| 国产精品人妻久久久影院| 少妇人妻一区二区三区视频| 国产午夜福利久久久久久| 国产免费男女视频| 黄色欧美视频在线观看| 精品日产1卡2卡| 国产人妻一区二区三区在| 久久6这里有精品| a级毛片免费高清观看在线播放| 一a级毛片在线观看| 欧美+亚洲+日韩+国产| 三级男女做爰猛烈吃奶摸视频| 国产精品人妻久久久久久| 国产大屁股一区二区在线视频| 夜夜夜夜夜久久久久| 亚洲最大成人手机在线| 91在线精品国自产拍蜜月| 久久人人爽人人爽人人片va| av.在线天堂| 搞女人的毛片| 麻豆一二三区av精品| 级片在线观看| 午夜福利成人在线免费观看| 精品人妻熟女av久视频| 国产精品久久视频播放| 亚洲av一区综合| 色哟哟·www| 国产精品,欧美在线| 97人妻精品一区二区三区麻豆| 国产精品综合久久久久久久免费| 女同久久另类99精品国产91| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品电影一区二区三区| 国产在视频线在精品| 欧美成人精品欧美一级黄| 一级黄片播放器| 简卡轻食公司| 一本久久中文字幕| 欧美zozozo另类| 黄色视频,在线免费观看| 国产亚洲91精品色在线| 免费av毛片视频| 看非洲黑人一级黄片| 日产精品乱码卡一卡2卡三| 12—13女人毛片做爰片一| 亚洲人成网站在线播放欧美日韩| 天天一区二区日本电影三级| 亚洲精华国产精华液的使用体验 | 亚洲av五月六月丁香网| 欧美又色又爽又黄视频| 久久久久久久久久黄片| 狠狠狠狠99中文字幕| 人妻少妇偷人精品九色| 亚洲第一电影网av| 三级经典国产精品| 精品福利观看| 波野结衣二区三区在线| 伦理电影大哥的女人| 国产精品国产高清国产av| 伊人久久精品亚洲午夜| 欧美+日韩+精品| 国产精品99久久久久久久久| 久久亚洲国产成人精品v| 午夜日韩欧美国产| 日韩欧美在线乱码| 少妇的逼水好多| 在线国产一区二区在线| 中文字幕av在线有码专区| 天堂av国产一区二区熟女人妻| 天天躁夜夜躁狠狠久久av| 99在线人妻在线中文字幕| 国产成人影院久久av| 午夜精品国产一区二区电影 | 精品人妻偷拍中文字幕| 男人舔奶头视频| 国产精品无大码| 欧美日韩精品成人综合77777| 在线观看一区二区三区| 草草在线视频免费看| 草草在线视频免费看| 亚洲综合色惰| 国产一区亚洲一区在线观看| 久久人人爽人人爽人人片va| 国产精品三级大全| a级毛片a级免费在线| 成人漫画全彩无遮挡| 欧美日韩乱码在线| 69av精品久久久久久| 成人无遮挡网站| 给我免费播放毛片高清在线观看| 午夜福利视频1000在线观看| 麻豆国产av国片精品| 51国产日韩欧美| 麻豆av噜噜一区二区三区| 国产精品伦人一区二区| 国产91av在线免费观看| 天堂影院成人在线观看| 在线播放国产精品三级| 成人二区视频| 蜜桃久久精品国产亚洲av| 老女人水多毛片| av专区在线播放| 91在线观看av| 色综合色国产| 中文亚洲av片在线观看爽| 国产高潮美女av| 亚洲精品一区av在线观看| 黄片wwwwww| 美女黄网站色视频| 久久久久久久午夜电影| 欧美性猛交╳xxx乱大交人| а√天堂www在线а√下载| 精品久久久久久久末码| 俺也久久电影网| 成人性生交大片免费视频hd| 亚洲成a人片在线一区二区| 亚洲av电影不卡..在线观看| 国产在线男女| 内射极品少妇av片p| 婷婷精品国产亚洲av| 麻豆久久精品国产亚洲av| 亚洲av中文av极速乱| 少妇被粗大猛烈的视频| 黑人高潮一二区| 露出奶头的视频| 亚洲国产精品成人综合色| 给我免费播放毛片高清在线观看| 午夜激情欧美在线| 欧美另类亚洲清纯唯美| 观看免费一级毛片| 久久亚洲国产成人精品v| 一个人看的www免费观看视频| 色哟哟哟哟哟哟| 亚洲美女黄片视频| 99热网站在线观看| 狂野欧美激情性xxxx在线观看| 欧美在线一区亚洲| 老司机影院成人| 极品教师在线视频| 亚洲在线自拍视频| 国产伦一二天堂av在线观看| 国产视频一区二区在线看| 亚洲精品日韩在线中文字幕 | 国产毛片a区久久久久| 寂寞人妻少妇视频99o| 日韩欧美精品v在线| 亚洲av免费在线观看| 日韩欧美免费精品| 午夜免费激情av| 美女高潮的动态| 国产成人a区在线观看| 日韩欧美在线乱码| 国产精品亚洲一级av第二区| 国产精品久久久久久亚洲av鲁大| 久久这里只有精品中国| 欧美性猛交黑人性爽| 一边摸一边抽搐一进一小说| 麻豆乱淫一区二区| 国产色爽女视频免费观看| 国产精品99久久久久久久久| 成人特级黄色片久久久久久久| 2021天堂中文幕一二区在线观| 免费av毛片视频| 我要搜黄色片| 男人和女人高潮做爰伦理| 精品午夜福利视频在线观看一区| 国产探花在线观看一区二区| 色综合色国产| 婷婷精品国产亚洲av在线| 亚州av有码| 欧美激情在线99| 久久人人爽人人片av| 男女那种视频在线观看| 亚洲国产精品合色在线| 91午夜精品亚洲一区二区三区| 九九久久精品国产亚洲av麻豆| 黄色视频,在线免费观看| 欧美xxxx黑人xx丫x性爽| 又爽又黄无遮挡网站| 国产av麻豆久久久久久久| 99视频精品全部免费 在线| 成人三级黄色视频| av黄色大香蕉| 免费av观看视频| 久久精品国产自在天天线| 最近手机中文字幕大全| 日本三级黄在线观看| 中文字幕人妻熟人妻熟丝袜美| 级片在线观看| 男女那种视频在线观看| 国产在线精品亚洲第一网站| 国产色爽女视频免费观看| 久久人妻av系列| 丝袜喷水一区| 成人综合一区亚洲| 日韩一区二区视频免费看| 一级毛片aaaaaa免费看小| 亚洲丝袜综合中文字幕| 看黄色毛片网站| 日韩成人伦理影院| 国产精品亚洲美女久久久| 熟妇人妻久久中文字幕3abv| 国产探花极品一区二区| 全区人妻精品视频| 成人毛片a级毛片在线播放| 亚洲乱码一区二区免费版| 国产亚洲91精品色在线| 午夜福利在线观看吧| 婷婷色综合大香蕉| 国国产精品蜜臀av免费| 最近的中文字幕免费完整| 老熟妇乱子伦视频在线观看| or卡值多少钱| av专区在线播放| 色吧在线观看| 春色校园在线视频观看| 久久久久久久亚洲中文字幕| av.在线天堂| 免费不卡的大黄色大毛片视频在线观看 | 久久久久性生活片| 国产精品人妻久久久久久| 又粗又爽又猛毛片免费看| 国产精品三级大全| 国产成人aa在线观看| a级一级毛片免费在线观看| 天天躁夜夜躁狠狠久久av| 久久久久久国产a免费观看| 亚洲成人久久爱视频| a级一级毛片免费在线观看| 日韩人妻高清精品专区| 久久久久久国产a免费观看| 日韩 亚洲 欧美在线| 亚洲av二区三区四区| 99久久无色码亚洲精品果冻| 国产探花极品一区二区| 久久久精品欧美日韩精品| 亚洲第一电影网av| 成人亚洲欧美一区二区av| 色吧在线观看| 久久午夜亚洲精品久久| 天堂网av新在线| 成人毛片a级毛片在线播放| 精品久久国产蜜桃| 少妇人妻精品综合一区二区 | 亚洲内射少妇av| 三级毛片av免费| 国产精品不卡视频一区二区| 亚洲自拍偷在线| 91久久精品国产一区二区成人| 欧美高清性xxxxhd video| 女同久久另类99精品国产91| 日韩三级伦理在线观看| 精品人妻熟女av久视频| 亚洲国产日韩欧美精品在线观看| 国产真实伦视频高清在线观看| 亚洲aⅴ乱码一区二区在线播放| 别揉我奶头 嗯啊视频| 欧美区成人在线视频| 亚洲熟妇中文字幕五十中出| 别揉我奶头~嗯~啊~动态视频| 天堂av国产一区二区熟女人妻| 国产男靠女视频免费网站| 岛国在线免费视频观看| 精品人妻一区二区三区麻豆 | 嫩草影视91久久| 日本免费a在线| 高清毛片免费看| 寂寞人妻少妇视频99o| or卡值多少钱| 欧美极品一区二区三区四区| 在线观看免费视频日本深夜| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲网站| 大香蕉久久网| 深夜精品福利| 国产视频内射| 亚洲经典国产精华液单| 久久人人爽人人爽人人片va| 麻豆精品久久久久久蜜桃| 成年版毛片免费区| 国产男人的电影天堂91| 99热这里只有是精品在线观看| 亚洲国产精品合色在线| 亚洲一级一片aⅴ在线观看| 三级毛片av免费| 色综合色国产| 九色成人免费人妻av| 男人狂女人下面高潮的视频| 九九在线视频观看精品| 淫秽高清视频在线观看| 久久精品国产亚洲av天美| 91久久精品电影网| 国产黄色小视频在线观看| 国产aⅴ精品一区二区三区波| 97碰自拍视频| 免费搜索国产男女视频| 亚洲丝袜综合中文字幕| eeuss影院久久| 成人av在线播放网站| 国产精品一区二区三区四区免费观看 | 色综合亚洲欧美另类图片| 一进一出抽搐动态| 国产精品一区二区三区四区免费观看 | 久久婷婷人人爽人人干人人爱| 成人午夜高清在线视频| 国产精品伦人一区二区| 亚洲性久久影院| 欧美日本亚洲视频在线播放| 露出奶头的视频| 非洲黑人性xxxx精品又粗又长| 亚洲丝袜综合中文字幕| 两个人的视频大全免费| 天堂√8在线中文| 精品欧美国产一区二区三| 乱码一卡2卡4卡精品| 欧美最新免费一区二区三区| 国产精品av视频在线免费观看| 91久久精品国产一区二区成人| 色av中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 91午夜精品亚洲一区二区三区| 色综合亚洲欧美另类图片| aaaaa片日本免费| 不卡一级毛片| 美女内射精品一级片tv| 成人永久免费在线观看视频| 长腿黑丝高跟| 又爽又黄a免费视频| 亚洲成人中文字幕在线播放| 最近的中文字幕免费完整| 国产精品久久久久久久电影| 91在线精品国自产拍蜜月| 亚洲人成网站在线观看播放| 亚洲av免费在线观看| 精品久久久久久成人av| 最好的美女福利视频网| 在线免费十八禁| 成人特级黄色片久久久久久久| 欧美激情在线99| 欧美+日韩+精品| 成人特级av手机在线观看|