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

    基于核偏最小二乘法的濕法冶金萃取過程建模

    2021-06-02 13:09:06
    濕法冶金 2021年3期
    關(guān)鍵詞:傳質(zhì)水相數(shù)據(jù)模型

    陳 溥

    (柳州鐵道職業(yè)技術(shù)學(xué)院,廣西 柳州 545616)

    目前,國內(nèi)外有關(guān)銅溶劑萃取過程中的建模研究較少,建立使用模型對于萃取過程中的自動控制有重要意義[1-5]。建立關(guān)鍵生產(chǎn)指標(biāo)預(yù)測模型,對萃取過程的控制和優(yōu)化至關(guān)重要,但影響過程操作的重要動力學(xué)傳質(zhì)速率無法測量。研究提出了一種適用于銅萃取過程中的結(jié)合機理模型和數(shù)據(jù)模型的混合模型。該模型由參數(shù)未知(傳質(zhì)速率)的機理模型和參數(shù)未知的數(shù)據(jù)預(yù)測模型組合而成。利用核偏最小二乘法(kernel partial least square,KPLS)建立傳質(zhì)速率預(yù)測模型,通過仿真與BP(back propagation)神經(jīng)網(wǎng)絡(luò)的混合模型進行對比分析。

    1 銅溶劑萃取過程概述

    銅萃取過程可分解為萃取、洗滌、反萃取3個階段[6-7]。

    萃取階段:含有銅離子的水相與有機相混合,銅離子被萃入到有機相[8]。洗滌階段:負(fù)載有機相用洗滌劑洗滌,將夾帶的金屬離子洗出有機相。反萃取階段:用反萃取劑對洗滌后有機相進行反萃取,使銅離子返回水相。

    影響銅溶劑萃取過程的因素有相比、溫度、水相pH、萃取劑濃度、被萃取金屬離子濃度等[9-10]。

    2 混合建模

    銅萃取過程機理模型的未知參數(shù)實際上無法測量,所以,提出了一種混合模型。該模型將機理模型和數(shù)據(jù)模型相結(jié)合,由參數(shù)(傳質(zhì)速率)未知的機理模型和參數(shù)未知的數(shù)據(jù)預(yù)測模型組成[11]。用KPLS法分別建立數(shù)據(jù)預(yù)測模型。

    2.1 動態(tài)機理模型

    動態(tài)機理模型描述萃取過程中組分濃度隨時間的變化,是溶劑萃取過程控制系統(tǒng)的依據(jù),有助于實現(xiàn)組分濃度優(yōu)化控制。

    在建立動態(tài)模型之前,對于多級逆流萃取過程做出以下假設(shè):

    1)兩相完全混合,并且僅在混合器中進行傳質(zhì);

    2)兩相不混溶,澄清器沒有返混和逆流現(xiàn)象;

    3)傳質(zhì)系數(shù)、萃取平衡等溫線參數(shù)和效率參數(shù)不恒定,根據(jù)離線測量估算[12];

    4)萃取過程中兩相體積不變;

    5)萃取過程中各組分濃度隨時間變化,是時間的函數(shù)。

    根據(jù)物料平衡關(guān)系,建立第i階段萃取過程的動力學(xué)模型,見式(1)~(3)[13]。

    (1)

    (2)

    ρoi(t)=f(ρa(i+1)(t-ti),ρa(i+1)(t),qa(t),qo(t),A,B)。

    (3)

    式中:Vi—i級混合澄清器中混合相體積,L;qa、qo—水相、有機相流量,m3/h;ρa(i+1)為i+1級澄清器出口水相中Cu2+質(zhì)量濃度,g/L;ρo(i-1)為i-1級澄清器出口的有機相Cu2+質(zhì)量濃度,g/L;t、ti—萃取、滯后時間,h;Ki為第i級傳質(zhì)系數(shù),L/s;ρoi(t)-ρo*i(t)為傳質(zhì)的驅(qū)動力,kN/m3;ρo*i為第i級達到平衡時有機相中Cu2+質(zhì)量濃度(理想狀態(tài)),g/L,與上層Cu2+進入水相的濃度、下一階段進入有機相的濃度、水相和有機相的流量、效率參數(shù)(α)和平衡萃取線常數(shù)A、B之間存在未知的函數(shù)關(guān)系,當(dāng)i=N時,Ki[ρoi(t)-ρo*i(t)]為傳質(zhì)速率,g/(L·h)。

    用McCabe-Thiele圖計算平衡關(guān)系確定萃取級數(shù),萃取部分理論平衡值根據(jù)水相和有機相中Cu2+及水相和有機相流量比確定。萃取平衡值x*計算公式[14]為

    (4)

    式中:a=-qa/qo;b=ρo(i-1)-aρa(i+1),g/L。

    α受有機相中試劑濃度和水相料液酸度影響。根據(jù)理論濃度平衡值x*(水相)和y*(有機相)、輸入的水相濃度和有機相中Cu2+質(zhì)量濃度ρin及實際Cu2+質(zhì)量濃度平衡值ρout估算每個萃取單元的萃取效率。見式(5),該值接近于1[15]。

    (5)

    當(dāng)萃取達到平衡時,萃取過程處于穩(wěn)定狀態(tài),見式(6)[16]:

    q(tend)[ρin(tend)-ρout(tend)]-
    KiV[ρout(tend)-ρ*out(tend)]=0,

    (6)

    一般來說,ρout(tend)-ρ*out(tend)非常小。用ε表示,ρout(t)-ρ*out(t),導(dǎo)出傳質(zhì)系數(shù)K,見式(7)[17]:

    (7)

    傳質(zhì)系數(shù)受多種因素影響,如操作條件、反應(yīng)器尺寸等。在交叉參考微反應(yīng)器中,銅的傳質(zhì)系數(shù)為31.12~214.34 L/s[18];混合澄清器的K=50~400 L/s。

    2.2 BP神經(jīng)網(wǎng)絡(luò)模型

    BP神經(jīng)網(wǎng)絡(luò)有3層,輸入層、輸出層和隱藏層,每層有多個節(jié)點[19]。節(jié)點是神經(jīng)元,神經(jīng)元之間的連接線為權(quán)重?;綛P網(wǎng)絡(luò)如圖1所示。

    圖1 基本BP神經(jīng)網(wǎng)絡(luò)模型

    假設(shè)輸入層、輸出層、隱含層節(jié)點數(shù)分別為r、m、n;xi為輸入層i神經(jīng)元的輸入值;zk、sk分別為隱含層第k個節(jié)點的輸出值和凈輸入值;yj、sj分別為輸出層j節(jié)點的實際輸出和凈輸入值;vik、wkj分別為輸入層i節(jié)點與隱含層k節(jié)點和隱含層k節(jié)點與輸出層j節(jié)點之間的連接權(quán)值。

    正向傳輸過程中,輸出層和隱含層的輸出計算公式見(8)~(11)[20]。

    zk=f(sk),k=1,2,……,n;

    (8)

    (9)

    yj=f(sj),j=1,2,……,m;

    (10)

    (11)

    在式(8)~(11)中,選擇激勵函數(shù)f(x)=1/(1+e-x),其導(dǎo)數(shù)為f’=f(1-f)。通常,誤差函數(shù)(也稱性能函數(shù))E可以評估網(wǎng)絡(luò)性能,E越小,網(wǎng)絡(luò)性能越好。見式(12),性能函數(shù)取均方誤差函數(shù)。

    (12)

    誤差反向傳播將誤差定義擴展到隱含層至輸入層,見式(13)。

    (13)

    標(biāo)準(zhǔn)BP算法用梯度下降法來調(diào)整層間權(quán)重,見式(14)、(15)。

    (14)

    (15)

    如果η較大,則修正速度快且易振蕩;η較小時,則會導(dǎo)致網(wǎng)絡(luò)訓(xùn)練時間長且收斂速度慢。η取值范圍為0.01~0.8。BP算法流程如圖2所示。

    圖2 BP算法流程

    2.3 核偏最小二乘法

    偏最小二乘法(partial least square,PLS)是一種多元統(tǒng)計數(shù)據(jù)分析方法,可以解決普通多元回歸方法無法解決的許多問題[21]。偏最小二乘方法假定感興趣區(qū)域由一些潛在向量決定,這些向量是觀測值的線性組合。假設(shè)有向量個因變量{y1,y2,…,yq}和p個自變量{x1,x2,…,xp},觀測n個樣本點。X={x1,x2,…,xp}為輸入變量的n×p矩陣;Y={y1,y2,…,yq}為相應(yīng)輸出變量的(n×p)矩陣。

    首先,從X和Y中提取得分向量t1和u1(t1為x1,x2,…,xp的線性組合,u1為y1,y2,…,yq的線性組合),要求t1和u1盡可能代表矩陣X和Y。

    t1=Xw1;

    (16)

    u1=Yc1。

    (17)

    然后,建立X對t1、Y對u1的回歸方程,見式(18)(19),且滿足式(20)。

    X=t1p1T+E1;

    (18)

    Y=u1q1T+F1;

    (19)

    u1=f1(t1)+r1,

    (20)

    式中,p1、q1分別為得分向量對應(yīng)的載荷向量。

    殘差計算公式為(21)(22)。

    E1=X-t1p1T;

    (21)

    F1=Y-u1q1T。

    (22)

    提取E1和F1得分向量,重復(fù)上述算法,直到得到第a個得分向量。殘差Ea和Fa幾乎沒有任何關(guān)鍵信息。

    由于銅萃取過程的動力學(xué)反應(yīng)為非線性,PLS算法無法解決該問題,因此引進核偏最小二乘法。該方法是在核希爾伯特空間(Reproducing Kemel Hilbert Space,RKHS)中建立的機器學(xué)習(xí)算法之一。用該方法建立模型的目的是最小化結(jié)構(gòu)風(fēng)險,即在(RKHS)中尋找最優(yōu)模型。KPLS遵循PLS基本思想,它們的區(qū)別在于PLS建立從屬變量(y)和自變量(x)之間的相關(guān)性,而KPLS僅將自變量(x)投影到RKHS上,并且在RKHS中找到Ф(x)和y的相關(guān)性。KPLS的本質(zhì)是在特征空間中構(gòu)建PLS模型,有效地獲得非線性回歸模型。

    同樣,觀測n個樣本,所謂“核技術(shù)”,實際上是Ф(xi)TФ(xj)=K(xi,xj)??梢钥吹?,ФФT為所有映射的輸入數(shù)據(jù)點{Ф(xi)}n(i=1)之間交叉內(nèi)積運算的核Gram矩陣K(n×n),因此,可以用核函數(shù)來代替非線性映射。通常,當(dāng)提取t分量后,矩陣見式(23)。

    K←(I-ttT)K(I-ttT)=
    K-ttTK-KttT+ttTKttT,

    (23)

    式中,I為n維單位矩陣。對PLS算法的非線性校正也可以得到類似的KPLS算法。修正的結(jié)果是從KYYT和YYT矩陣中提取得分向量t、u,步驟如下:

    步驟1,隨機設(shè)定初始值u(可設(shè)定為u等于Y列中的任意一列);

    步驟2,根據(jù)t=ФФTu=Ku,計算得分向量;

    步驟3,根據(jù)c=YTt計算系數(shù)向量;

    步驟4,計算得分向量;

    步驟5,重復(fù)2~5步,直到收斂;

    步驟6,矩陣K和Y進行縮并,K←(I-ttT)K(I-ttT),Y←Y-ttTY;

    步驟7,返回步驟1,計算下一個t、u。

    回歸系數(shù)矩陣B通過公式(24)計算得到,

    B=ФTU(TTKU)-1TTY。

    (24)

    預(yù)測的訓(xùn)練數(shù)據(jù)見式(25),

    Y*=ФB=TTTY。

    (25)

    KPLS算法流程如圖3所示。

    圖3 KPLS算法流程

    3 仿真結(jié)果與分析

    3.1 仿真參數(shù)

    混合模型由參數(shù)未知(傳質(zhì)速率)的機理模型和參數(shù)未知的數(shù)據(jù)預(yù)測模型組成。通過BP神經(jīng)網(wǎng)絡(luò)和KPLS建立傳質(zhì)速率預(yù)測模型,并與機理模型形成混合模型。數(shù)據(jù)模型是黑盒模型。在未知傳質(zhì)速率模型結(jié)構(gòu)和相關(guān)參數(shù)情況下,可以獲得傳質(zhì)速率背后的過程變量之間的復(fù)雜關(guān)系,從而降低建模難度。在考慮過程的物理特性的同時,機理模型和數(shù)據(jù)模型組合有效地利用了相關(guān)的數(shù)據(jù)信息。這種混合模式具有很高的精度,適合于實際工業(yè)生產(chǎn)過程,為銅溶劑萃取工藝先進控制系統(tǒng)的開發(fā)和在線優(yōu)化奠定了基礎(chǔ)。銅萃取混合模型如圖4所示。

    圖4 銅萃取過程混合模型

    計算機為CPU i5-2450m,內(nèi)存為8 Gb。仿真軟件為matlab2018。實驗室實際采集200組萃取過程穩(wěn)態(tài)數(shù)據(jù),100組用于訓(xùn)練,100組用于測試。為了驗證和比較上述預(yù)測模型性能,利用以下3個性能指標(biāo)對各算法性能進行評價。

    1)均方根誤差(RMSE)

    (26)

    2)平均絕對誤差(MAE)

    (27)

    3)平均絕對百分比誤差(MAPE)

    (28)

    式中:xi和x*i分別為銅傳質(zhì)速率實際值和預(yù)測值,M為預(yù)測樣品總數(shù)。

    3.2 仿真分析

    建立傳質(zhì)速率黑箱模型,輸出傳質(zhì)速率預(yù)測值。黑箱模型的輸入對應(yīng)于傳輸速率估計值的影響因素。相同條件下,利用BP神經(jīng)網(wǎng)絡(luò)和KPLS對傳輸速率進行建模。圖5為BP神經(jīng)網(wǎng)絡(luò)及KPLS黑箱模型預(yù)測效果??梢钥闯觯築P神經(jīng)網(wǎng)絡(luò)的訓(xùn)練效果良好,表明BP神經(jīng)網(wǎng)絡(luò)黑箱模型具有準(zhǔn)確的預(yù)測能力,非線性映射能力較好;但訓(xùn)練中收斂速度較慢,并且在訓(xùn)練集比較小時,容易發(fā)生過擬合;KPLS將初始輸入映射到高維特征空間,有效解決了非線性回歸問題。

    a—BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型;b—KPLS預(yù)測模型。圖5 黑箱預(yù)測結(jié)果

    2個數(shù)據(jù)模型性能指標(biāo)見表1。

    表1 BP神經(jīng)網(wǎng)絡(luò)和KPLS數(shù)據(jù)模型預(yù)測性能指標(biāo)

    數(shù)據(jù)模型的正確性在混合模型構(gòu)建中起重要作用。以上2種算法能夠準(zhǔn)確預(yù)測傳質(zhì)速率。由表1性能指標(biāo)對比結(jié)果看出,KPLS方法效果較好。數(shù)據(jù)模型的準(zhǔn)確性解決了傳質(zhì)速率無法測量的問題,這為建立更接近實際過程的混合模型奠定了基礎(chǔ)。

    對比單一模型和混合模型預(yù)測結(jié)果可以說明上述銅萃取過程混合模型的有效性,對比結(jié)果如圖6所示,各模型的性能指標(biāo)見表2。

    a—BP神經(jīng)網(wǎng)絡(luò);b—KPLS。

    表2 單一模型和混合模型預(yù)測性能指標(biāo)

    由圖6看出:機理模型能夠很好預(yù)測過程趨勢;然而,由于萃取過程較為復(fù)雜,預(yù)測結(jié)果與實際值誤差較大;單一數(shù)據(jù)模型相對于機理模型需要考慮的因素更多,預(yù)測效果優(yōu)于機理模型。由于數(shù)據(jù)模型不使用已知的經(jīng)驗知識,僅依賴數(shù)據(jù)獲得輸入輸出之間的關(guān)系,純數(shù)據(jù)模型的泛化能力差,所以預(yù)測值和實際值誤差較大?;旌夏P蛯?shù)據(jù)模型和機理模型優(yōu)點相結(jié)合,精度較高,泛化能力強,有良好的預(yù)測效果,對銅溶液萃取工藝先進控制系統(tǒng)的開發(fā)和操作優(yōu)化有一定參考性?;贙PLS的混合模型性能較優(yōu),RMSE為0.066 52,MAE為0.043 12,MAPE為0.019 82;對銅溶液萃取工藝先進控制系統(tǒng)的開發(fā)和操作優(yōu)化有一定參考性。

    4 結(jié)論

    提出了一種混合模型,該模型由存在未知參數(shù)(傳質(zhì)速率)的機理模型和未知參數(shù)的數(shù)據(jù)預(yù)測模型組成,通過數(shù)據(jù)模型預(yù)測機理模型的未知參數(shù)(傳質(zhì)速率度)?;旌夏P托阅艿姆治鼋Y(jié)果表明,基于KPLS的混合模型性能較優(yōu),均方根誤差、平均絕對誤差、平均絕對百分比誤差都較低,可作為銅萃取過程中自動控制系統(tǒng)開發(fā)的依據(jù)??紤]到實驗室設(shè)備和數(shù)據(jù)規(guī)模,建模方法還處于初級階段,還有待進一步完善。

    猜你喜歡
    傳質(zhì)水相數(shù)據(jù)模型
    海上中高滲透率砂巖油藏油水相滲曲線合理性綜合分析技術(shù)
    更 正
    面板數(shù)據(jù)模型截面相關(guān)檢驗方法綜述
    加熱爐爐內(nèi)跟蹤數(shù)據(jù)模型優(yōu)化
    電子測試(2017年12期)2017-12-18 06:35:36
    地下水流速與介質(zhì)非均質(zhì)性對于重非水相流體運移的影響
    氨基酸鹽吸收二氧化碳過程的傳質(zhì)特性
    PTFE膜吸收CO2的工藝參數(shù)對傳質(zhì)性能的影響
    用三辛胺和磷酸三丁酯萃取、銨溶液反萃取鉬的研究
    濕法冶金(2014年3期)2014-04-08 01:04:51
    清潔轉(zhuǎn)向酸H+表面?zhèn)髻|(zhì)行為實驗研究
    煤顆粒熱解的傳熱傳質(zhì)分析
    計算物理(2014年1期)2014-03-11 17:00:36
    国精品久久久久久国模美| 色哟哟·www| www.色视频.com| 国产精品一区二区在线不卡| 日日撸夜夜添| 国产一区亚洲一区在线观看| 激情五月婷婷亚洲| 最近中文字幕高清免费大全6| 黄色怎么调成土黄色| 亚洲综合色网址| 爱豆传媒免费全集在线观看| 国产极品粉嫩免费观看在线| 久久久久久伊人网av| 99热网站在线观看| 丰满迷人的少妇在线观看| 国产综合精华液| 黄片播放在线免费| 免费观看无遮挡的男女| 亚洲精品av麻豆狂野| 人人妻人人添人人爽欧美一区卜| 日韩av免费高清视频| 国产成人免费观看mmmm| 亚洲av男天堂| 色婷婷av一区二区三区视频| 国产成人免费观看mmmm| 一边亲一边摸免费视频| 亚洲欧洲日产国产| 男女高潮啪啪啪动态图| 亚洲欧洲精品一区二区精品久久久 | 日韩精品免费视频一区二区三区 | 亚洲国产欧美在线一区| 热99国产精品久久久久久7| 天天躁夜夜躁狠狠躁躁| 色网站视频免费| 少妇人妻 视频| 亚洲,一卡二卡三卡| 丝袜在线中文字幕| 成人免费观看视频高清| 香蕉国产在线看| 黑人欧美特级aaaaaa片| 欧美亚洲日本最大视频资源| 亚洲中文av在线| 观看av在线不卡| 青青草视频在线视频观看| 亚洲五月色婷婷综合| 亚洲三级黄色毛片| 交换朋友夫妻互换小说| 国产国拍精品亚洲av在线观看| 国产在视频线精品| 国产精品久久久av美女十八| 国产福利在线免费观看视频| 国产视频首页在线观看| 蜜桃国产av成人99| 国产免费现黄频在线看| 18禁裸乳无遮挡动漫免费视频| 人妻人人澡人人爽人人| 久久久久久久亚洲中文字幕| 99国产综合亚洲精品| 欧美日韩精品成人综合77777| 国产精品女同一区二区软件| 女的被弄到高潮叫床怎么办| 亚洲中文av在线| 亚洲美女视频黄频| 亚洲高清免费不卡视频| 亚洲国产精品国产精品| 亚洲av电影在线进入| 欧美人与性动交α欧美精品济南到 | av电影中文网址| 一本久久精品| av国产精品久久久久影院| av在线播放精品| 久久热在线av| 男女国产视频网站| 精品少妇久久久久久888优播| 亚洲情色 制服丝袜| 亚洲,一卡二卡三卡| 中文字幕精品免费在线观看视频 | 全区人妻精品视频| 母亲3免费完整高清在线观看 | 国产在线免费精品| 免费观看av网站的网址| 五月伊人婷婷丁香| 青春草亚洲视频在线观看| 亚洲欧美清纯卡通| 秋霞伦理黄片| 天天操日日干夜夜撸| 在线观看三级黄色| 青春草国产在线视频| 亚洲国产精品999| 街头女战士在线观看网站| 丰满迷人的少妇在线观看| 成人二区视频| 亚洲经典国产精华液单| 精品久久蜜臀av无| 国产精品国产三级国产专区5o| 久久毛片免费看一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲欧美成人综合另类久久久| 咕卡用的链子| 久热这里只有精品99| 一级,二级,三级黄色视频| 国产片内射在线| 啦啦啦啦在线视频资源| 成年人免费黄色播放视频| 日本午夜av视频| 国产成人免费无遮挡视频| 男女午夜视频在线观看 | 久久人人爽人人片av| 久久久久久久久久成人| 七月丁香在线播放| 97精品久久久久久久久久精品| 18禁国产床啪视频网站| 菩萨蛮人人尽说江南好唐韦庄| 丰满少妇做爰视频| 精品亚洲乱码少妇综合久久| 久久久久久久久久久久大奶| 国产精品欧美亚洲77777| a级毛色黄片| 黄色一级大片看看| 有码 亚洲区| 99热6这里只有精品| 亚洲性久久影院| 国产国语露脸激情在线看| 午夜91福利影院| 国产成人精品一,二区| 亚洲av欧美aⅴ国产| 亚洲欧洲精品一区二区精品久久久 | 汤姆久久久久久久影院中文字幕| 99久久综合免费| 高清毛片免费看| 成年美女黄网站色视频大全免费| 精品第一国产精品| 天天躁夜夜躁狠狠躁躁| 男女边摸边吃奶| 日本爱情动作片www.在线观看| 亚洲精品一区蜜桃| 在线看a的网站| 男女国产视频网站| 极品少妇高潮喷水抽搐| 亚洲一区二区三区欧美精品| 久久97久久精品| 亚洲图色成人| 日韩av在线免费看完整版不卡| 久久人人爽av亚洲精品天堂| 亚洲国产精品国产精品| 大香蕉久久成人网| 亚洲欧美清纯卡通| 91精品三级在线观看| 国产国语露脸激情在线看| 我要看黄色一级片免费的| 精品一品国产午夜福利视频| 久久毛片免费看一区二区三区| 中文欧美无线码| 久久 成人 亚洲| 亚洲第一av免费看| 国产精品成人在线| 日韩免费高清中文字幕av| 男人添女人高潮全过程视频| 蜜臀久久99精品久久宅男| 国产成人免费无遮挡视频| 国产成人精品一,二区| 亚洲精品久久久久久婷婷小说| xxxhd国产人妻xxx| 久久久a久久爽久久v久久| 成人毛片a级毛片在线播放| 多毛熟女@视频| 中文乱码字字幕精品一区二区三区| 久久人人97超碰香蕉20202| 精品人妻一区二区三区麻豆| 涩涩av久久男人的天堂| 午夜福利,免费看| 久久免费观看电影| av又黄又爽大尺度在线免费看| 日日摸夜夜添夜夜爱| a 毛片基地| 亚洲国产精品国产精品| 国产极品天堂在线| 2021少妇久久久久久久久久久| 国产午夜精品一二区理论片| 亚洲综合色惰| 国产亚洲一区二区精品| 欧美日韩一区二区视频在线观看视频在线| 9色porny在线观看| 天天影视国产精品| 国产黄频视频在线观看| 在现免费观看毛片| 亚洲中文av在线| 少妇的逼好多水| 超碰97精品在线观看| 一级毛片黄色毛片免费观看视频| av国产精品久久久久影院| 免费黄频网站在线观看国产| 午夜av观看不卡| 99久久精品国产国产毛片| 久久av网站| 亚洲av.av天堂| 韩国av在线不卡| 色94色欧美一区二区| 2018国产大陆天天弄谢| 黑人巨大精品欧美一区二区蜜桃 | 激情视频va一区二区三区| 免费高清在线观看日韩| xxxhd国产人妻xxx| 午夜免费男女啪啪视频观看| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美一区二区三区国产| 精品国产一区二区三区久久久樱花| 国产亚洲最大av| 亚洲精品456在线播放app| 国产免费现黄频在线看| 国产成人av激情在线播放| 一级爰片在线观看| 精品一区二区三区四区五区乱码 | a级毛色黄片| 日韩欧美一区视频在线观看| 国产av码专区亚洲av| 又大又黄又爽视频免费| 青青草视频在线视频观看| 在线免费观看不下载黄p国产| 中文字幕免费在线视频6| 日韩在线高清观看一区二区三区| 国产免费一级a男人的天堂| 亚洲成国产人片在线观看| 最近中文字幕高清免费大全6| 99热这里只有是精品在线观看| 少妇被粗大的猛进出69影院 | 欧美另类一区| 黄色视频在线播放观看不卡| 乱人伦中国视频| 免费大片黄手机在线观看| 国产黄频视频在线观看| 天美传媒精品一区二区| 亚洲成人一二三区av| 宅男免费午夜| 一区在线观看完整版| 女人精品久久久久毛片| 中文乱码字字幕精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 亚洲一码二码三码区别大吗| 热re99久久国产66热| 国产精品不卡视频一区二区| 我的女老师完整版在线观看| 成人二区视频| 2022亚洲国产成人精品| 久久97久久精品| 另类亚洲欧美激情| 久久人人爽人人爽人人片va| 男人舔女人的私密视频| 1024视频免费在线观看| a 毛片基地| 国产精品一区二区在线观看99| 精品第一国产精品| 午夜福利乱码中文字幕| 18禁裸乳无遮挡动漫免费视频| 欧美3d第一页| 纵有疾风起免费观看全集完整版| 啦啦啦中文免费视频观看日本| 男女无遮挡免费网站观看| 观看av在线不卡| 人妻少妇偷人精品九色| 你懂的网址亚洲精品在线观看| 午夜免费鲁丝| 中国美白少妇内射xxxbb| 少妇的丰满在线观看| 亚洲精品国产av蜜桃| 人成视频在线观看免费观看| 亚洲四区av| 捣出白浆h1v1| 久久ye,这里只有精品| 高清在线视频一区二区三区| 成人手机av| 亚洲四区av| av有码第一页| 欧美性感艳星| 另类亚洲欧美激情| 国产激情久久老熟女| 亚洲欧美日韩另类电影网站| 欧美激情国产日韩精品一区| 黄色配什么色好看| 亚洲美女视频黄频| 丰满少妇做爰视频| 免费少妇av软件| 亚洲精品乱码久久久久久按摩| 美女xxoo啪啪120秒动态图| 激情视频va一区二区三区| 亚洲,欧美,日韩| 亚洲美女黄色视频免费看| 人人澡人人妻人| 久热久热在线精品观看| 日韩,欧美,国产一区二区三区| 国产一区二区在线观看日韩| 一二三四中文在线观看免费高清| 999精品在线视频| 捣出白浆h1v1| 极品人妻少妇av视频| 王馨瑶露胸无遮挡在线观看| av在线app专区| 亚洲国产精品一区二区三区在线| 欧美另类一区| 啦啦啦啦在线视频资源| 亚洲经典国产精华液单| 亚洲欧洲日产国产| 夜夜爽夜夜爽视频| 中国三级夫妇交换| 欧美 日韩 精品 国产| 美女大奶头黄色视频| 国语对白做爰xxxⅹ性视频网站| 亚洲国产日韩一区二区| 校园人妻丝袜中文字幕| 久久av网站| 精品一区在线观看国产| 一级爰片在线观看| 国产精品国产三级国产av玫瑰| 亚洲国产欧美日韩在线播放| 制服诱惑二区| 99热这里只有是精品在线观看| 欧美精品一区二区大全| 午夜福利影视在线免费观看| 精品人妻在线不人妻| 香蕉国产在线看| 大片免费播放器 马上看| 91成人精品电影| 在线看a的网站| xxx大片免费视频| 国内精品宾馆在线| 国产亚洲av片在线观看秒播厂| 99视频精品全部免费 在线| 亚洲av福利一区| 久久99热6这里只有精品| 成人亚洲精品一区在线观看| 久久久久久伊人网av| 97精品久久久久久久久久精品| 国产在线视频一区二区| 日韩不卡一区二区三区视频在线| 日本欧美国产在线视频| 97超碰精品成人国产| 久久精品久久久久久噜噜老黄| 欧美人与性动交α欧美软件 | 美女国产视频在线观看| 国产精品久久久久久精品电影小说| 国产伦理片在线播放av一区| 国产麻豆69| 狂野欧美激情性bbbbbb| 国产女主播在线喷水免费视频网站| 99久国产av精品国产电影| 精品国产一区二区三区久久久樱花| 高清不卡的av网站| 不卡视频在线观看欧美| 国产成人av激情在线播放| 国产精品国产av在线观看| 黑人高潮一二区| 免费黄网站久久成人精品| 大香蕉久久成人网| 黄色 视频免费看| 日韩av不卡免费在线播放| 女性被躁到高潮视频| 国产高清不卡午夜福利| 又黄又粗又硬又大视频| 色网站视频免费| 美女福利国产在线| 国产探花极品一区二区| 啦啦啦中文免费视频观看日本| 国产探花极品一区二区| 最近中文字幕高清免费大全6| 尾随美女入室| 青春草亚洲视频在线观看| 亚洲欧美成人精品一区二区| 国产亚洲最大av| 人妻 亚洲 视频| av片东京热男人的天堂| 免费少妇av软件| 黑人高潮一二区| 国产亚洲精品第一综合不卡 | 欧美变态另类bdsm刘玥| 日本猛色少妇xxxxx猛交久久| 丝袜喷水一区| 国产男女内射视频| 精品一区二区三卡| √禁漫天堂资源中文www| 国产淫语在线视频| 午夜久久久在线观看| 成人午夜精彩视频在线观看| www.熟女人妻精品国产 | 久久97久久精品| 免费看光身美女| 亚洲,一卡二卡三卡| 在线精品无人区一区二区三| 国产成人精品婷婷| 亚洲精品国产av蜜桃| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久成人av| 欧美老熟妇乱子伦牲交| 在线观看三级黄色| 国产高清三级在线| 国产欧美日韩综合在线一区二区| 亚洲五月色婷婷综合| 插逼视频在线观看| 午夜福利视频精品| 性高湖久久久久久久久免费观看| 丝袜美足系列| 大陆偷拍与自拍| 一本大道久久a久久精品| 精品熟女少妇av免费看| a级毛片黄视频| 91午夜精品亚洲一区二区三区| 亚洲人与动物交配视频| 在线观看免费日韩欧美大片| 狠狠精品人妻久久久久久综合| 最新中文字幕久久久久| av在线app专区| 亚洲精品色激情综合| 巨乳人妻的诱惑在线观看| 午夜激情久久久久久久| av视频免费观看在线观看| 午夜久久久在线观看| 亚洲人成网站在线观看播放| 午夜免费鲁丝| videosex国产| 亚洲av免费高清在线观看| 满18在线观看网站| 国产日韩欧美视频二区| 日韩伦理黄色片| 欧美成人午夜免费资源| 国产一区亚洲一区在线观看| 午夜免费男女啪啪视频观看| 啦啦啦中文免费视频观看日本| 久久精品人人爽人人爽视色| 99久久精品国产国产毛片| 中文天堂在线官网| 嫩草影院入口| 久久久久久久久久成人| 国产色爽女视频免费观看| av播播在线观看一区| 男人爽女人下面视频在线观看| 人人妻人人爽人人添夜夜欢视频| 男女高潮啪啪啪动态图| 男人添女人高潮全过程视频| 国产成人免费观看mmmm| 亚洲精品国产av蜜桃| 制服诱惑二区| 久久精品久久久久久久性| 2022亚洲国产成人精品| 激情视频va一区二区三区| 欧美日韩成人在线一区二区| 99热全是精品| 男女啪啪激烈高潮av片| 亚洲精品日本国产第一区| 美女脱内裤让男人舔精品视频| 亚洲国产最新在线播放| 国产欧美亚洲国产| 美女xxoo啪啪120秒动态图| 另类亚洲欧美激情| 亚洲成国产人片在线观看| 黄色 视频免费看| 最新的欧美精品一区二区| 一级毛片我不卡| 岛国毛片在线播放| 成人漫画全彩无遮挡| 欧美精品高潮呻吟av久久| 亚洲av成人精品一二三区| 国产精品一二三区在线看| 天天操日日干夜夜撸| 久久国内精品自在自线图片| 亚洲av电影在线进入| 久久99热这里只频精品6学生| a级片在线免费高清观看视频| 国产男人的电影天堂91| 91成人精品电影| 美女国产视频在线观看| 纯流量卡能插随身wifi吗| freevideosex欧美| 亚洲欧洲精品一区二区精品久久久 | 欧美精品一区二区免费开放| 国产精品99久久99久久久不卡 | 国产在视频线精品| 国产在线免费精品| 丰满迷人的少妇在线观看| 秋霞伦理黄片| 亚洲精品色激情综合| 黑人巨大精品欧美一区二区蜜桃 | 女性被躁到高潮视频| www.av在线官网国产| 全区人妻精品视频| 久久99热这里只频精品6学生| 欧美变态另类bdsm刘玥| 大片电影免费在线观看免费| videos熟女内射| freevideosex欧美| 久久韩国三级中文字幕| 久久久a久久爽久久v久久| 视频区图区小说| 欧美3d第一页| 亚洲精品自拍成人| 最近手机中文字幕大全| 深夜精品福利| 欧美少妇被猛烈插入视频| 一级爰片在线观看| 我的女老师完整版在线观看| 在线亚洲精品国产二区图片欧美| 视频在线观看一区二区三区| 在线免费观看不下载黄p国产| 欧美精品人与动牲交sv欧美| 老司机影院成人| 欧美日韩一区二区视频在线观看视频在线| 久久国产亚洲av麻豆专区| 亚洲国产日韩一区二区| 久久热在线av| 18+在线观看网站| 国产视频首页在线观看| 国产成人午夜福利电影在线观看| 国产亚洲欧美精品永久| 又黄又爽又刺激的免费视频.| 免费播放大片免费观看视频在线观看| 亚洲五月色婷婷综合| 超色免费av| 国产精品久久久久久久久免| 国产精品成人在线| 超碰97精品在线观看| 亚洲色图 男人天堂 中文字幕 | 天堂中文最新版在线下载| 韩国av在线不卡| 一区二区av电影网| 久久毛片免费看一区二区三区| 精品少妇久久久久久888优播| 天堂中文最新版在线下载| a级毛色黄片| 精品一区二区免费观看| 少妇被粗大的猛进出69影院 | 亚洲国产精品999| av国产久精品久网站免费入址| 男女边摸边吃奶| 女的被弄到高潮叫床怎么办| 精品午夜福利在线看| 国产在线视频一区二区| 免费日韩欧美在线观看| 日韩中字成人| 天天影视国产精品| 久久97久久精品| 大话2 男鬼变身卡| 中文字幕亚洲精品专区| 最近最新中文字幕免费大全7| 国产毛片在线视频| 春色校园在线视频观看| 在线观看免费日韩欧美大片| 在线观看人妻少妇| 黄色怎么调成土黄色| 欧美xxⅹ黑人| 中文字幕制服av| 成人国产av品久久久| 亚洲欧美中文字幕日韩二区| 一区二区三区精品91| videos熟女内射| 免费大片18禁| 久久久久久久久久久免费av| 日韩 亚洲 欧美在线| 美女主播在线视频| 两个人免费观看高清视频| 美女福利国产在线| 亚洲精品久久久久久婷婷小说| 亚洲精品一区蜜桃| 亚洲情色 制服丝袜| 人人妻人人澡人人爽人人夜夜| 亚洲人与动物交配视频| 七月丁香在线播放| 国产不卡av网站在线观看| 精品国产国语对白av| 亚洲av中文av极速乱| 男女午夜视频在线观看 | 国产av一区二区精品久久| 精品少妇内射三级| 中文天堂在线官网| 精品亚洲乱码少妇综合久久| 久久综合国产亚洲精品| 精品国产露脸久久av麻豆| 国产69精品久久久久777片| 啦啦啦在线观看免费高清www| 五月开心婷婷网| 久久久久网色| 国产一区二区在线观看av| 国产毛片在线视频| 免费日韩欧美在线观看| 亚洲欧洲国产日韩| 亚洲三级黄色毛片| 18禁在线无遮挡免费观看视频| 2022亚洲国产成人精品| 日本欧美视频一区| 亚洲精品av麻豆狂野| freevideosex欧美| 久久久a久久爽久久v久久| 亚洲色图综合在线观看| 精品人妻偷拍中文字幕| 日日啪夜夜爽| xxxhd国产人妻xxx| 一级片免费观看大全| 在线观看人妻少妇| 校园人妻丝袜中文字幕| 久久久亚洲精品成人影院| 国产爽快片一区二区三区| 久久影院123| 制服丝袜香蕉在线| 亚洲图色成人| 欧美性感艳星| 成人毛片a级毛片在线播放| av线在线观看网站| 大片电影免费在线观看免费| 女性生殖器流出的白浆| 国产精品熟女久久久久浪| 国产欧美日韩综合在线一区二区| 久久精品aⅴ一区二区三区四区 | 美女福利国产在线| 国产精品国产三级专区第一集| 欧美精品高潮呻吟av久久| 毛片一级片免费看久久久久| 欧美亚洲 丝袜 人妻 在线| 日韩电影二区|