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

    一種改進(jìn)的遺傳算法在間歇化工過程設(shè)計中的應(yīng)用

    2020-08-22 09:41:36魏良霄黨樂平
    化學(xué)工業(yè)與工程 2020年4期
    關(guān)鍵詞:間歇平行遺傳算法

    魏良霄,張 婷,黨樂平,郝 琳

    (天津大學(xué)化工學(xué)院,天津 300072)

    為了在全球市場上保持競爭力,優(yōu)化生產(chǎn)流程的設(shè)計和管理是重中之重。 間歇化工過程以其技術(shù)密集性、動態(tài)性、多樣性、靈活性和不確定性等特點在化工及相關(guān)行業(yè)中發(fā)揮著重要作用[1]。 雖然在間歇過程的設(shè)計和優(yōu)化方面均有了很大的進(jìn)展,但由于這是一個離散變量和連續(xù)變量共存的混合整數(shù)非線性規(guī)劃MINLP (Mixed-Integer Nonlinear Programming)問題[2],求解過程中目標(biāo)函數(shù)的復(fù)雜度隨著變量的數(shù)量變化呈指數(shù)增長趨勢。 所以對問題進(jìn)行建模和求解仍然相當(dāng)困難[3]。

    差分運算是間歇過程設(shè)計中最早使用的方法[4],但解決方案的復(fù)雜性使得該方法很快被分支定界法[5-6]、數(shù)學(xué)規(guī)劃法和其他啟發(fā)式方法所取代[7-8]。 然而,這些方法即使是在解決簡單問題的過程中就非常麻煩,顯然不適合解決大規(guī)模問題。人工智能算法或通過多種方法相結(jié)合的方式來解決上述問題獲得了成功的應(yīng)用。 Loonkar 和Robinson[9]使用迭代方法來最小化設(shè)備投資。 Yeh等[10]使用啟發(fā)式算法討論單一產(chǎn)品化學(xué)過程的綜合問題。 鄭國文等[11]基于Yeh 的研究通過添加中間存儲單元,進(jìn)一步研究了這一過程。 Kirkpatrick等[12]提出了模擬退火算法(SAA),Yuan 等[13]在間歇過程的設(shè)計和調(diào)度中使用了該算法。 楊志才等[14]提出了一種分解策略,然后引入了基于該策略和隨機(jī)算法的組合算法來解決間歇過程優(yōu)化的MINLP 問題。 正如Moniz 等[15]所指出的,我們?nèi)匀恍枰碌姆椒?它應(yīng)該對啟發(fā)式方法和智能優(yōu)化算法進(jìn)行整合,以解決現(xiàn)實中更加復(fù)雜的問題。

    本研究提出了一種基于遺傳算法(GA)的MINLP 模型,并用回溯法對該模型進(jìn)行了補充。 該算法在保持遺傳算法優(yōu)點的前提下結(jié)合回溯法指導(dǎo)搜索過程,消除劣質(zhì)基因;并且在建模過程中,考慮了設(shè)備利用率最大化問題。 最后通過實例與利用枚舉法和回溯法相結(jié)合的方法建立的模型進(jìn)行了比較。

    1 遺傳算法的基本原理

    遺傳算法是由美國密歇根大學(xué)的Holland[16]提出的。 它是解決復(fù)雜系統(tǒng)問題的框架,具有良好的適用性和魯棒性。 一個典型的遺傳算法的基本要求是:1)解域的遺傳表示。 2)用于評估解域的適應(yīng)度函數(shù)。

    每個候選解決方案的標(biāo)準(zhǔn)表示形式是位數(shù)組。其他類型和結(jié)構(gòu)的數(shù)組基本上可以以相同的方式使用。 如圖1 所示,遺傳算法的過程包括以下步驟[17]:1)創(chuàng)建一組編碼方案。 2)使用適應(yīng)度函數(shù)評估不同的字符串。 3)字符串根據(jù)其評估的適應(yīng)度大小進(jìn)行排序。 4)然后創(chuàng)建2 個子節(jié)點,對2 個字符串進(jìn)行交叉和編譯操作,然后替換上一代的2個字符串。 5) 最后重復(fù)該循環(huán),直到達(dá)到停止標(biāo)準(zhǔn)。

    2 數(shù)學(xué)模型

    2.1 假設(shè)

    圖1 遺傳算法典型流程圖Fig.1 Typical structure of genetic algorithm

    為降低模型建立時的復(fù)雜程度,我們做如下假設(shè):1)生產(chǎn)工藝和配方已知;2)總生產(chǎn)時間和生產(chǎn)需求已經(jīng)明確;3)有充足的原料供應(yīng);4)設(shè)備在整個時間段內(nèi)無故障;5)設(shè)備為非標(biāo)設(shè)備;6)設(shè)備數(shù)量不加以限制;7)處理時間已經(jīng)確定,不會改變;8)處理時間包括進(jìn)料時間、出料時間、反應(yīng)時間和清洗時間。

    2.2 約束方程

    為了提高設(shè)備的利用率,規(guī)定可以將前一間歇級的產(chǎn)品分為xi次生產(chǎn),即設(shè)置中間存儲,如圖2 間歇級2 所示。 很明顯如果間歇級1 中的物料尚未被間歇級2 處理完,間歇級1 是不能接受新的物料的投入的,否則會造成中間存儲中物料的積壓。 于是可以將這個分次生產(chǎn)過程等效為一個處理時間較長的過程,因此有:

    式(1)中:ti為產(chǎn)品在間歇級i 的處理時間,h; t′i為產(chǎn)品在間歇級i 的等效處理時間,h;xi為上一級產(chǎn)品在在下一級的分批生產(chǎn)次數(shù)。

    間歇化工過程通常包含1 個或多個瓶頸階段。通常這些階段涉及到更長的處理時間或者資本密集型設(shè)備[18]。 將與瓶頸階段相對應(yīng)的時間稱為循環(huán)時間:

    式(2)中:ni為間歇級i 的異步平行單元數(shù);T 為限制循環(huán)時間,h。

    圖2 分段生產(chǎn)示意圖Fig.2 Schematic diagram of sectional production

    在1 臺設(shè)備上每批應(yīng)完成的產(chǎn)品量為:

    式(3)中:mi為間歇級i 的同步平行單元數(shù);Bi為間歇級i 每1 批完成的產(chǎn)品量,kg;Qi為間歇級i 的生產(chǎn)總量,kg;H 為年生產(chǎn)時間,h。

    設(shè)備的體積與設(shè)備加工的產(chǎn)品量之間的關(guān)系為:

    式(4)中:Si為間歇級i 的尺寸因子,kg/m3;Vi為設(shè)備體積,m3。

    并且有:

    引入平行單元以降低這些瓶頸階段的主導(dǎo)地位[19]。同步平行單元的最大數(shù)值被指定為在非覆蓋式下操作時需要設(shè)置的同步平行單元的數(shù)。

    式(6)中: Vmax為設(shè)備最大體積,m3。

    反應(yīng)設(shè)備總數(shù)為:

    式(8)中: nmax為間歇級的最大異步平行單元數(shù)。

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

    該模型以設(shè)備投資最低為目標(biāo)函數(shù)[20],在計算過程中設(shè)備成本的估算通常以體積的指數(shù)函數(shù)表示[21]:

    式(9)中:αi為間歇級i 的成本系數(shù);βi為間歇級i的成本指數(shù)。

    故總投資為:

    3 實例計算

    3.1 問題描述

    以實際項目為例,該項目包含7 個間歇級,每個間歇級的出料都是固體,包含過濾和干燥2 種操作。因此,我們可以將每個間歇級的設(shè)備等效為1 個具有過濾和干燥功能的反應(yīng)器。 為了確保轉(zhuǎn)化率不受較大影響,將設(shè)備最大體積設(shè)為5 m3。 年生產(chǎn)時間為330 d,每天24 h 不間斷生產(chǎn),故總生產(chǎn)時間為7 920 h。 其他所有輸入數(shù)據(jù)在表1 中給出。

    表1 產(chǎn)品加工要求清單Table 1 List of product processing requirements

    3.2 算法實現(xiàn)

    該模型的主體為遺傳算法,有3 組變量:同步平行單元數(shù)、異步平行單元數(shù)和連續(xù)生產(chǎn)數(shù)。 其中連續(xù)生產(chǎn)是指在2 個連續(xù)批次的同一間歇級生產(chǎn)時間間隔可能較大的情況下前一間歇級中該步驟的產(chǎn)物分為幾次投入到此間歇級中進(jìn)行生產(chǎn),以減少設(shè)備的體積,提高設(shè)備的利用率。

    由于變量較多,生成的染色體較為復(fù)雜。 我們決定直接以十進(jìn)制數(shù)生成染色體,并且取消編碼和解碼過程。 該模型的每個染色體由14 個基因組成。分為2 部分,前7 個基因代表每個間歇級的同步平行單元數(shù),后7 個基因代表每個間歇級的異步平行單元數(shù)。 在突變過程中,應(yīng)根據(jù)基因所處的位置使用相應(yīng)的函數(shù)進(jìn)行突變,否則極大可能會產(chǎn)生劣質(zhì)基因。

    異步平行單元數(shù)最大值由經(jīng)驗指定,在本模型中指定為4。 然而當(dāng)使用這些基因計算反應(yīng)器的體積時,發(fā)現(xiàn)反應(yīng)器的體積超過了最大體積的極限。于是采用懲罰操作對進(jìn)化過程進(jìn)行優(yōu)化,簡單的流程如圖3 所示。

    圖3 初始模型流程圖Fig.3 Flow chart of initial model

    程序運行后最優(yōu)值為9.146×105,但是在重新運行模型后發(fā)現(xiàn)最優(yōu)值發(fā)生了變化,圖4a 是模型運行100 次后最優(yōu)值的波動情況。 從圖4a 中可以看出,最優(yōu)值一直在變化,其中最低值為8.998×105,最高值為10.580×105,增加了17.6%。

    圖4 最優(yōu)值波動曲線Fig.4 Fluctuation curve of optimal solution after optimization

    經(jīng)過計算,發(fā)現(xiàn)隨機(jī)種群生成后,許多染色體不僅不能滿足限制條件的約束,甚至出現(xiàn)嚴(yán)重錯誤。 比如在這些染色體的指導(dǎo)下會出現(xiàn)上一批物料尚未加工完畢便又投入新的物料的現(xiàn)象,這明顯不符合實際情況。 并且這些染色體在種群中占很高的比例,這使得在種群在進(jìn)化過程中很難獲得更優(yōu)秀的基因,并且覆蓋范圍變得狹窄,這使得GA 的特性無法正常發(fā)揮。 圖5a 顯示了100 代種群進(jìn)化過程中每一代中有缺陷基因的比例。 從圖5a 中我們可以看出,在幾乎每一代中,有缺陷基因的比例都超過80%甚至接近100%。

    圖5 每一代中有缺陷基因的比例Fig.5 The proportion of faulty genes in each generation

    因此,在計算適應(yīng)度前必須對基因進(jìn)行優(yōu)化。于是建立了一種自尋優(yōu)機(jī)制,每次形成新的種群時都對染色體進(jìn)行判斷,并在出現(xiàn)錯誤時對染色體進(jìn)行優(yōu)化。 簡單的流程如圖6 所示。

    再次運行模型100 次后發(fā)現(xiàn)最優(yōu)值幾乎穩(wěn)定在8.998×105,如圖4b 所示。 即使在最優(yōu)值發(fā)生變化的地方,最大值也只有9.072×105,只比8.998×105高0.822 %。 并且大部分錯誤基因出現(xiàn)的比例下降到了50%以下,如圖5b 所示。 當(dāng)然自優(yōu)化機(jī)制的作用不僅僅是這樣,它真正的作用是在錯誤基因出現(xiàn)后對其進(jìn)行優(yōu)化,而錯誤基因出現(xiàn)的比例下降是由于上一代錯誤的基因減少造成的,這是一種額外效果。

    3.3 結(jié)果

    該程序在MATLAB2017b 中運行,運行時間為0.634 s。 經(jīng)過一定代數(shù)的繁衍后收斂到8.998 ×105這一值。 限制循環(huán)時間為120 h,其他數(shù)據(jù)如表2 所示。

    將計算結(jié)果繪制成甘特圖如圖7 所示。 從圖7中可以看出,間歇級1 有7 個同步平行單元,間歇級2 有2 個同步平行單元,整個過程中沒有異步平行單元。 間歇級1 是時間限制級。 中間存儲器設(shè)置在間歇級2 中,并且調(diào)用間歇級2 的2 個批次生產(chǎn)時間間隔從原來的120 h 變?yōu)楫?dāng)前的24 h。 每個階段反應(yīng)器的利用率如圖8 所示。

    圖6 自優(yōu)化機(jī)制流程圖Fig.6 Flow chart of self-optimization mechanism

    表2 遺傳算法計算結(jié)果Table 2 Calculation results of genetic algorithm

    圖7 利用遺傳算法得到的甘特圖Fig.7 Gantt chart obtained by genetic algorithm

    圖8 利用遺傳算法得到的設(shè)備利用率Fig.8 Reactor utilization obtained by genetic algorithm

    4 與一般啟發(fā)式方法進(jìn)行比較

    本模型參考鄭國文等[11]所建立的模型,利用枚舉法和回溯法進(jìn)行建模。 首先計算所有滿足約束方程的x 的值,然后將它們帶入圖9 所示的流程圖中來搜索局部最優(yōu)值。 最后從所有局部最適宜值中選擇最小的結(jié)果作為最優(yōu)結(jié)果。

    圖9 啟發(fā)式算法流程圖Fig.9 Flow chart of heuristic method

    該程序在MATLAB 2017b 中運行并運行14 s。程序運行后,設(shè)備投資局部最優(yōu)值分布折線如圖10所示,其中橫坐標(biāo)代表不同的中間存儲方案,縱坐標(biāo)代表確定存儲方案后模型求得的最優(yōu)解。 從圖10 中可以看出,設(shè)備投資成本的分布范圍較大,最小值出現(xiàn)在3 個位置。 它們分別位于橫坐標(biāo)57、88和96 上,所有與最優(yōu)值相關(guān)的計算結(jié)果列在表3 中。

    圖10 局部最優(yōu)值分布曲線Fig.10 Distribution curve of optimal solution

    從表3 中可以看出,3 個結(jié)果的同步平行單元數(shù)和異步平行單元數(shù)相同,但是連續(xù)生產(chǎn)次數(shù)不同。 顯然,實際生產(chǎn)中在間歇級1 前設(shè)置中間存儲是沒有意義的,因此我們只選取列表中x1=1 的情況,即僅保留表3 中的第一種情況。

    表3 啟發(fā)式算法計算結(jié)果Table 3 Calculation results of heuristic algorithm

    這種生產(chǎn)方法以甘特圖的形式表示,見圖11。從圖11 中可以看出,間歇級1 為過程中的時間限制級,包含4 個同步平行單元4 個異步平行單元。 異步平行單元使限制時間由原來的120 h 減少到了60 h。 在間歇級2 中設(shè)置中間存儲,使得上一間歇級的物料分兩批加入此間歇級,從而減小后續(xù)階段反應(yīng)器的體積,增大反應(yīng)器利用率。 可以得出每個階段的反應(yīng)器利用率,如圖12 所示。

    比較以上2 個結(jié)果,可以看出,遺傳算法獲得的設(shè)備投資結(jié)果比啟發(fā)式算法獲得的結(jié)果低2. 66%。 這是因為同步平行單元數(shù)和異步平行單元數(shù)是隨機(jī)生成的,并且在后續(xù)的交叉變異等操作過程中使得每個階段的同步平行單元數(shù)和異步平行單元數(shù)更加具有多樣性,因此優(yōu)化過程中遺傳算法的覆蓋范圍比一般啟發(fā)式算法的覆蓋范圍更廣。 并且一般啟發(fā)式算法的同步平行單元數(shù)和異步平行單元數(shù)都是按照一定的規(guī)律進(jìn)行變化的,這使得在優(yōu)化過程中易于陷入局部最優(yōu)解。用遺傳算法建立的模型在計算速度上是用一般啟發(fā)式算法建立的模型的22 倍。 但是,當(dāng)中間存儲的容量變小時,2 者之間的時間差距逐漸減小,這表明要解決的問題越復(fù)雜,遺傳算法的優(yōu)勢越明顯。 但是,設(shè)備的利用率沒有變化。 這是因為用一般啟發(fā)式算法建立的模型使用枚舉法來優(yōu)化設(shè)備的利用率,并且所有可能性都包含在指定范圍內(nèi),因此無法進(jìn)一步優(yōu)化。

    圖11 利用一般啟發(fā)式方法得到的甘特圖Fig.11 Gantt chart obtained using general heuristic method

    圖12 利用一般啟發(fā)式方法得到的設(shè)備利用率Fig.12 Reactor utilization obtained using general heuristic method

    5 結(jié)論

    介紹了一種改進(jìn)的遺傳算法在解決單產(chǎn)品間歇化工過程設(shè)計問題中的應(yīng)用。 通過引入自我優(yōu)化機(jī)制和懲罰操作在不影響遺傳算法特性的前提下對基因進(jìn)行人為干預(yù),從而減少種群中劣質(zhì)染色體的出現(xiàn),使種群得以順利穩(wěn)定的繁衍下去。 最后與一般啟發(fā)式算法建立的模型進(jìn)行比較。 結(jié)果表明該模型在計算結(jié)果、收斂速度和計算速度上均優(yōu)于一般啟發(fā)式方法。

    猜你喜歡
    間歇平行遺傳算法
    間歇供暖在散熱器供暖房間的應(yīng)用
    煤氣與熱力(2022年4期)2022-05-23 12:44:46
    向量的平行與垂直
    平行
    逃離平行世界
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    基于遺傳算法和LS-SVM的財務(wù)危機(jī)預(yù)測
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    再頂平行進(jìn)口
    汽車觀察(2016年3期)2016-02-28 13:16:36
    基于改進(jìn)的遺傳算法的模糊聚類算法
    中文字幕av成人在线电影| 联通29元200g的流量卡| 一进一出抽搐gif免费好疼| av专区在线播放| 很黄的视频免费| 亚洲最大成人手机在线| 69人妻影院| 天美传媒精品一区二区| 天堂动漫精品| 欧美潮喷喷水| 欧美xxxx黑人xx丫x性爽| 久久久久久国产a免费观看| 亚洲真实伦在线观看| 欧美日韩国产亚洲二区| 亚洲人与动物交配视频| 亚洲av一区综合| 丰满的人妻完整版| 神马国产精品三级电影在线观看| 亚洲中文日韩欧美视频| 女人被狂操c到高潮| 久久久色成人| 国产亚洲精品久久久com| 午夜老司机福利剧场| 51国产日韩欧美| 欧美高清性xxxxhd video| 最近最新中文字幕大全电影3| 伊人久久精品亚洲午夜| 又粗又爽又猛毛片免费看| 午夜影院日韩av| 亚洲av日韩精品久久久久久密| 亚洲aⅴ乱码一区二区在线播放| 久9热在线精品视频| 久久久久久久亚洲中文字幕| 老熟妇乱子伦视频在线观看| 成人三级黄色视频| 丝袜美腿在线中文| 日本黄色片子视频| 欧美精品啪啪一区二区三区| 村上凉子中文字幕在线| 美女黄网站色视频| 久久午夜亚洲精品久久| 亚洲avbb在线观看| 国产伦一二天堂av在线观看| 最好的美女福利视频网| 床上黄色一级片| 国产精品嫩草影院av在线观看 | 国产久久久一区二区三区| 亚洲最大成人手机在线| 国产又黄又爽又无遮挡在线| 国产毛片a区久久久久| 97人妻精品一区二区三区麻豆| 亚洲av成人av| 99在线视频只有这里精品首页| 99久久成人亚洲精品观看| netflix在线观看网站| 亚洲精品亚洲一区二区| 国产av不卡久久| 亚洲aⅴ乱码一区二区在线播放| 色av中文字幕| 成人性生交大片免费视频hd| 一进一出抽搐动态| 欧美日韩中文字幕国产精品一区二区三区| 久久99热这里只有精品18| 亚洲精品日韩av片在线观看| 中文资源天堂在线| 亚洲人成网站在线播| 午夜激情欧美在线| a在线观看视频网站| 色精品久久人妻99蜜桃| 69av精品久久久久久| 亚洲国产精品成人综合色| 国产精品无大码| 欧美日本亚洲视频在线播放| 久久欧美精品欧美久久欧美| 麻豆久久精品国产亚洲av| 国产伦在线观看视频一区| 看片在线看免费视频| 免费黄网站久久成人精品| 又黄又爽又刺激的免费视频.| 国产免费男女视频| 91av网一区二区| а√天堂www在线а√下载| 国产精品一区www在线观看 | 国产 一区 欧美 日韩| 久久久久久久久久成人| 国产精品乱码一区二三区的特点| 五月伊人婷婷丁香| 欧美高清成人免费视频www| 久久99热6这里只有精品| 午夜精品久久久久久毛片777| 一区二区三区免费毛片| 亚洲国产欧美人成| 国产精品自产拍在线观看55亚洲| 免费观看人在逋| av在线天堂中文字幕| 中文字幕久久专区| 国产单亲对白刺激| 麻豆国产av国片精品| 超碰av人人做人人爽久久| 少妇被粗大猛烈的视频| 亚洲av中文字字幕乱码综合| 精品久久久久久久人妻蜜臀av| 精品国内亚洲2022精品成人| 久久久久久久久久久丰满 | 亚州av有码| 久久精品国产自在天天线| 国产精品综合久久久久久久免费| av天堂在线播放| 成人毛片a级毛片在线播放| 此物有八面人人有两片| 日本色播在线视频| 日本五十路高清| 久久久久国产精品人妻aⅴ院| netflix在线观看网站| 亚洲无线观看免费| netflix在线观看网站| 亚洲天堂国产精品一区在线| 级片在线观看| 亚洲精品乱码久久久v下载方式| 色视频www国产| 干丝袜人妻中文字幕| 色视频www国产| 一级av片app| 精品日产1卡2卡| 国产av一区在线观看免费| a级一级毛片免费在线观看| 麻豆av噜噜一区二区三区| 国产成人影院久久av| 欧美日韩中文字幕国产精品一区二区三区| 国内毛片毛片毛片毛片毛片| 欧美一区二区亚洲| a级毛片免费高清观看在线播放| 精品久久久久久久久av| 欧美一区二区亚洲| 岛国在线免费视频观看| 狠狠狠狠99中文字幕| 黄色视频,在线免费观看| 九九热线精品视视频播放| 亚洲av日韩精品久久久久久密| 观看美女的网站| 亚洲内射少妇av| 婷婷六月久久综合丁香| 色噜噜av男人的天堂激情| 热99在线观看视频| 日本色播在线视频| 国产精品无大码| 一a级毛片在线观看| 国内揄拍国产精品人妻在线| 免费电影在线观看免费观看| 久久久久久伊人网av| 久久久久久久久久成人| 天天一区二区日本电影三级| 国产黄片美女视频| 亚洲电影在线观看av| av中文乱码字幕在线| 又粗又爽又猛毛片免费看| 在线观看免费视频日本深夜| 国产高清不卡午夜福利| 成年版毛片免费区| 免费无遮挡裸体视频| 伦理电影大哥的女人| 精品福利观看| 亚洲精品国产成人久久av| 欧美zozozo另类| 国产精品久久久久久av不卡| 夜夜爽天天搞| 床上黄色一级片| www.色视频.com| 国产精品福利在线免费观看| avwww免费| 男人舔女人下体高潮全视频| 日韩欧美三级三区| 免费不卡的大黄色大毛片视频在线观看 | 国产三级中文精品| 日本五十路高清| 12—13女人毛片做爰片一| 99国产极品粉嫩在线观看| 久久久久久伊人网av| 女生性感内裤真人,穿戴方法视频| 欧美黑人巨大hd| 一本久久中文字幕| 国模一区二区三区四区视频| 免费观看人在逋| 我要看日韩黄色一级片| 欧美成人性av电影在线观看| 一个人观看的视频www高清免费观看| 中亚洲国语对白在线视频| 他把我摸到了高潮在线观看| 88av欧美| 极品教师在线免费播放| 综合色av麻豆| 久久国内精品自在自线图片| 日韩亚洲欧美综合| 热99re8久久精品国产| 老司机福利观看| 亚洲在线自拍视频| 网址你懂的国产日韩在线| 亚洲av.av天堂| 亚洲国产色片| 韩国av一区二区三区四区| 国产av在哪里看| 成人无遮挡网站| 国产在线精品亚洲第一网站| 中文字幕av成人在线电影| 99久久无色码亚洲精品果冻| 啦啦啦观看免费观看视频高清| 一区二区三区激情视频| 夜夜爽天天搞| 人妻丰满熟妇av一区二区三区| 日本-黄色视频高清免费观看| 欧美bdsm另类| 亚洲国产色片| 免费看光身美女| 波多野结衣高清作品| 亚洲精品久久国产高清桃花| 国产精品人妻久久久久久| 日韩一本色道免费dvd| 少妇裸体淫交视频免费看高清| 国产人妻一区二区三区在| 亚洲精品乱码久久久v下载方式| 一区二区三区激情视频| 免费大片18禁| 国产 一区 欧美 日韩| 免费在线观看影片大全网站| 男女之事视频高清在线观看| 成人一区二区视频在线观看| 亚洲成人中文字幕在线播放| 精品免费久久久久久久清纯| 赤兔流量卡办理| 精品久久久久久,| 99riav亚洲国产免费| 中文字幕久久专区| 国产 一区精品| 最近在线观看免费完整版| 久久久久久久久中文| av中文乱码字幕在线| 久久国内精品自在自线图片| 日日摸夜夜添夜夜添小说| 国内精品美女久久久久久| 久久国产乱子免费精品| а√天堂www在线а√下载| 日本在线视频免费播放| 97碰自拍视频| 精品午夜福利在线看| 久久草成人影院| av天堂中文字幕网| 国产欧美日韩精品亚洲av| 99热这里只有精品一区| 亚洲一区二区三区色噜噜| 91久久精品国产一区二区成人| 午夜久久久久精精品| 欧美日韩黄片免| 成人特级av手机在线观看| 99热这里只有是精品在线观看| 国产毛片a区久久久久| 国产精品国产三级国产av玫瑰| 国产成人av教育| 免费看光身美女| 亚洲av.av天堂| 蜜桃亚洲精品一区二区三区| 日本黄色视频三级网站网址| 国产精品自产拍在线观看55亚洲| 精品一区二区三区av网在线观看| 亚洲国产精品sss在线观看| 给我免费播放毛片高清在线观看| av在线亚洲专区| 亚洲美女视频黄频| 极品教师在线视频| 麻豆成人av在线观看| 久久人人精品亚洲av| 日韩亚洲欧美综合| 成人高潮视频无遮挡免费网站| 亚洲国产精品sss在线观看| 日韩高清综合在线| 亚洲精华国产精华液的使用体验 | 非洲黑人性xxxx精品又粗又长| 久久香蕉精品热| 有码 亚洲区| 尤物成人国产欧美一区二区三区| 欧美日本亚洲视频在线播放| 国产亚洲精品久久久com| 网址你懂的国产日韩在线| 国产亚洲91精品色在线| 色av中文字幕| 午夜福利18| 俺也久久电影网| 少妇熟女aⅴ在线视频| 精品乱码久久久久久99久播| 成年女人毛片免费观看观看9| 国产爱豆传媒在线观看| 最近最新免费中文字幕在线| 国产精品不卡视频一区二区| 天堂√8在线中文| 欧美色欧美亚洲另类二区| 中文在线观看免费www的网站| 久久精品人妻少妇| 日日干狠狠操夜夜爽| 欧美日本视频| 久久中文看片网| 精品欧美国产一区二区三| 日本成人三级电影网站| 我要搜黄色片| 国产69精品久久久久777片| 亚洲成人久久性| 日本爱情动作片www.在线观看 | 精品久久久久久久久久免费视频| 日韩欧美精品v在线| 精品欧美国产一区二区三| 国产精品乱码一区二三区的特点| 成人特级av手机在线观看| 久久精品国产亚洲av天美| 男人舔奶头视频| 亚洲精华国产精华精| 久久久久久久久久成人| АⅤ资源中文在线天堂| 午夜精品一区二区三区免费看| 亚洲成人中文字幕在线播放| 久久亚洲真实| 老熟妇乱子伦视频在线观看| 国国产精品蜜臀av免费| 美女被艹到高潮喷水动态| 少妇人妻一区二区三区视频| 亚洲四区av| 免费av不卡在线播放| 亚洲av.av天堂| 欧美日韩亚洲国产一区二区在线观看| 1000部很黄的大片| 春色校园在线视频观看| 成人国产综合亚洲| 久久久久久久精品吃奶| 色5月婷婷丁香| 日韩人妻高清精品专区| 99riav亚洲国产免费| 亚洲无线观看免费| 国产高清有码在线观看视频| 啦啦啦观看免费观看视频高清| 国产精品免费一区二区三区在线| 中出人妻视频一区二区| 搡老岳熟女国产| 国产淫片久久久久久久久| 亚洲中文字幕日韩| 亚洲aⅴ乱码一区二区在线播放| 国产探花极品一区二区| 九九热线精品视视频播放| 99精品在免费线老司机午夜| 精品乱码久久久久久99久播| 久久精品国产99精品国产亚洲性色| 久久欧美精品欧美久久欧美| 亚洲av免费高清在线观看| 成年女人看的毛片在线观看| 欧美激情在线99| 高清在线国产一区| 在线观看一区二区三区| 能在线免费观看的黄片| 91在线观看av| 琪琪午夜伦伦电影理论片6080| 一本精品99久久精品77| 美女免费视频网站| 狂野欧美白嫩少妇大欣赏| 日本成人三级电影网站| 日日干狠狠操夜夜爽| 日韩,欧美,国产一区二区三区 | 精品99又大又爽又粗少妇毛片 | 日本撒尿小便嘘嘘汇集6| 九九在线视频观看精品| 成人无遮挡网站| 不卡视频在线观看欧美| 亚洲18禁久久av| 精品久久久噜噜| 少妇被粗大猛烈的视频| 欧美极品一区二区三区四区| 久久欧美精品欧美久久欧美| 在线看三级毛片| 国产视频一区二区在线看| 亚洲av熟女| 淫妇啪啪啪对白视频| 在线观看一区二区三区| 天堂影院成人在线观看| 搡老岳熟女国产| 女的被弄到高潮叫床怎么办 | 日本一二三区视频观看| 一个人看视频在线观看www免费| 搡老妇女老女人老熟妇| 国产在线男女| 精品国产三级普通话版| 亚洲第一电影网av| 91av网一区二区| 亚洲av电影不卡..在线观看| 亚洲成人精品中文字幕电影| 欧美3d第一页| 人人妻,人人澡人人爽秒播| 无人区码免费观看不卡| 日本免费一区二区三区高清不卡| 69av精品久久久久久| 黄色一级大片看看| 午夜亚洲福利在线播放| 国产精品乱码一区二三区的特点| 黄色视频,在线免费观看| 亚洲,欧美,日韩| 国产av不卡久久| 窝窝影院91人妻| 91麻豆精品激情在线观看国产| 97超级碰碰碰精品色视频在线观看| 亚洲乱码一区二区免费版| 久久午夜亚洲精品久久| 欧美+亚洲+日韩+国产| 高清毛片免费观看视频网站| 三级国产精品欧美在线观看| 日本一二三区视频观看| 一进一出抽搐gif免费好疼| 少妇的逼好多水| 简卡轻食公司| h日本视频在线播放| 女的被弄到高潮叫床怎么办 | 99久久成人亚洲精品观看| 国产高清不卡午夜福利| 久久久精品欧美日韩精品| 联通29元200g的流量卡| 色哟哟哟哟哟哟| 舔av片在线| 日本撒尿小便嘘嘘汇集6| 国语自产精品视频在线第100页| 99热只有精品国产| x7x7x7水蜜桃| 精品不卡国产一区二区三区| 九九热线精品视视频播放| 又爽又黄无遮挡网站| 国产精品98久久久久久宅男小说| 88av欧美| 毛片女人毛片| 成人国产综合亚洲| 在线免费观看的www视频| 日韩 亚洲 欧美在线| 两性午夜刺激爽爽歪歪视频在线观看| 男女边吃奶边做爰视频| 国内精品久久久久久久电影| 国产伦一二天堂av在线观看| 九九久久精品国产亚洲av麻豆| 狂野欧美白嫩少妇大欣赏| 亚洲精华国产精华液的使用体验 | a级毛片免费高清观看在线播放| 日本在线视频免费播放| 12—13女人毛片做爰片一| 久久精品综合一区二区三区| 亚洲人与动物交配视频| 欧美日韩综合久久久久久 | 村上凉子中文字幕在线| 国产高清有码在线观看视频| 男人狂女人下面高潮的视频| 免费一级毛片在线播放高清视频| 神马国产精品三级电影在线观看| 亚洲七黄色美女视频| 国产高清三级在线| 婷婷色综合大香蕉| 俺也久久电影网| 男女那种视频在线观看| 日本-黄色视频高清免费观看| 成人毛片a级毛片在线播放| 久久久成人免费电影| 麻豆精品久久久久久蜜桃| 国产午夜精品论理片| 长腿黑丝高跟| 别揉我奶头~嗯~啊~动态视频| 日日夜夜操网爽| www.色视频.com| 久久久色成人| 在线播放国产精品三级| 永久网站在线| 欧美黑人欧美精品刺激| 观看美女的网站| 日日摸夜夜添夜夜添小说| 男人舔奶头视频| 观看美女的网站| 亚洲在线自拍视频| 国产高清激情床上av| 禁无遮挡网站| 国产成年人精品一区二区| 色吧在线观看| 最新在线观看一区二区三区| 在线观看av片永久免费下载| 日日干狠狠操夜夜爽| 久久久久久久亚洲中文字幕| 免费av不卡在线播放| 国产成年人精品一区二区| 简卡轻食公司| 九九在线视频观看精品| 有码 亚洲区| 免费看日本二区| 全区人妻精品视频| 久久6这里有精品| 欧美一级a爱片免费观看看| 一区二区三区激情视频| 婷婷色综合大香蕉| 俺也久久电影网| 99久国产av精品| 婷婷精品国产亚洲av| 美女免费视频网站| 一边摸一边抽搐一进一小说| 成人美女网站在线观看视频| 2021天堂中文幕一二区在线观| 天美传媒精品一区二区| 人妻丰满熟妇av一区二区三区| 搡女人真爽免费视频火全软件 | 亚洲精品一卡2卡三卡4卡5卡| 一区福利在线观看| 最新在线观看一区二区三区| 99视频精品全部免费 在线| 五月玫瑰六月丁香| 国产成人av教育| 午夜福利在线观看免费完整高清在 | 久久久久精品国产欧美久久久| 久久久久免费精品人妻一区二区| 精品乱码久久久久久99久播| 日韩亚洲欧美综合| 日韩大尺度精品在线看网址| 日本 av在线| 午夜激情福利司机影院| 黄色一级大片看看| 99久久久亚洲精品蜜臀av| 最近中文字幕高清免费大全6 | 日韩欧美国产在线观看| 久久久午夜欧美精品| 国产淫片久久久久久久久| 白带黄色成豆腐渣| АⅤ资源中文在线天堂| 国产在视频线在精品| av中文乱码字幕在线| or卡值多少钱| 免费av观看视频| av天堂在线播放| 国产欧美日韩精品亚洲av| 国产色爽女视频免费观看| 欧美日韩精品成人综合77777| 好男人在线观看高清免费视频| 亚洲经典国产精华液单| 狂野欧美激情性xxxx在线观看| 女人被狂操c到高潮| 九九爱精品视频在线观看| 国产午夜福利久久久久久| 欧美xxxx黑人xx丫x性爽| 亚洲国产日韩欧美精品在线观看| 精品乱码久久久久久99久播| 午夜精品在线福利| 熟女电影av网| 一夜夜www| 窝窝影院91人妻| 亚洲精华国产精华液的使用体验 | 久久香蕉精品热| 亚洲国产精品sss在线观看| 99久久精品国产国产毛片| 一个人看视频在线观看www免费| 免费观看人在逋| 中文字幕高清在线视频| 国产一区二区在线av高清观看| 99热这里只有精品一区| 午夜福利在线在线| 欧美成人一区二区免费高清观看| 国产精品99久久久久久久久| АⅤ资源中文在线天堂| 亚洲avbb在线观看| 一区福利在线观看| 国产中年淑女户外野战色| 全区人妻精品视频| 亚洲最大成人手机在线| 伦理电影大哥的女人| 色综合站精品国产| 免费一级毛片在线播放高清视频| 他把我摸到了高潮在线观看| 欧美日韩亚洲国产一区二区在线观看| 老熟妇仑乱视频hdxx| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品久久久久久毛片| 制服丝袜大香蕉在线| 欧美潮喷喷水| 中文字幕熟女人妻在线| 国产色婷婷99| 三级国产精品欧美在线观看| 亚洲国产色片| 国产熟女欧美一区二区| 国内毛片毛片毛片毛片毛片| 两性午夜刺激爽爽歪歪视频在线观看| 麻豆国产97在线/欧美| av黄色大香蕉| 村上凉子中文字幕在线| 韩国av在线不卡| 我的老师免费观看完整版| 99riav亚洲国产免费| 国产极品精品免费视频能看的| 亚洲五月天丁香| 亚洲无线观看免费| 亚洲av一区综合| 亚洲中文日韩欧美视频| 欧美日韩精品成人综合77777| 亚洲美女搞黄在线观看 | 丝袜美腿在线中文| 看片在线看免费视频| 国产精品精品国产色婷婷| 成人综合一区亚洲| www日本黄色视频网| 国内精品久久久久久久电影| 成人av一区二区三区在线看| 亚洲欧美日韩高清专用| 黄色日韩在线| 两个人视频免费观看高清| 国产精品1区2区在线观看.| 女人十人毛片免费观看3o分钟| 欧美高清性xxxxhd video| 国产精品免费一区二区三区在线| 亚洲最大成人手机在线| 欧美日韩黄片免| 国内精品久久久久久久电影| 少妇丰满av| 亚洲图色成人| 免费高清视频大片| 亚洲无线观看免费|