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

    發(fā)動機油中水分的近紅外光譜分析

    2023-02-22 06:25:54商芷萱全宇軒
    光譜學(xué)與光譜分析 2023年2期
    關(guān)鍵詞:油樣波長含水量

    劉 閣, 陳 彬, 商芷萱, 全宇軒

    1. 河北省危險化學(xué)品安全與控制技術(shù)重點實驗室, 華北科技學(xué)院化學(xué)與環(huán)境工程學(xué)院, 河北 燕郊 065201 2. 華北科技學(xué)院機電工程學(xué)院, 河北省礦山設(shè)備安全監(jiān)測重點實驗室, 河北 燕郊 065201

    引 言

    發(fā)動機潤滑油用于減少移動部件的摩擦, 并保持不同部件的清潔, 可以用作清潔劑和分散劑。 發(fā)動機在正常工作中極易在機油中產(chǎn)生污垢、 鹽、 水、 燃燒的不完全產(chǎn)物等污染物, 不但影響發(fā)動機油理化性能, 而且加速了添加劑、 基礎(chǔ)油的化學(xué)反應(yīng), 產(chǎn)生腐蝕性酸和水分等其他有害物質(zhì)[1]。 機油在發(fā)動機中的攪動、 熱和壓力會隨著時間的推移使水污染乳化。 因而發(fā)動機油中的水分主要以溶解或乳化狀態(tài)存在, 會導(dǎo)致腐蝕、 添加劑消耗、 氫致開裂和微點蝕等問題[2], 因此對發(fā)動機油中的水分進行監(jiān)測和量化是至關(guān)重要的。

    國內(nèi)外采用不同的技術(shù)對發(fā)動機油的水污染進行了大量的研究[3]。 開發(fā)了多種測定油中水分含量的方法。 由于添加劑會干擾卡爾-費歇爾法對這類油中水分的測定, 會導(dǎo)致高估水分含量。 核磁共振波譜用于發(fā)動機油的降解, 檢測水中的油污染[4]等各種研究, 對發(fā)動機油中區(qū)分相對較低的水濃度(0%, 1%, 2%V/V)之間有顯著的統(tǒng)計學(xué)差異[5]。

    由于在紅外光的作用下水的偶極矩發(fā)生變化, 氫原子具有最大的振動, 紅外光譜測試油中水分可能是傳統(tǒng)測水的有效替代方法, 得到了極大的關(guān)注。 雖然O—H健的光譜信號可能與其他功能基團的信號重疊, 但使用化學(xué)計量信號處理技術(shù)可以獲得準確的結(jié)果。 特別是自布格-蘭伯特-比爾提出光譜吸收定量評價的數(shù)學(xué)規(guī)律以來, 光譜分析已應(yīng)用于潤滑油性能檢測方面[6-7], 得到了迅速的發(fā)展。 Liu等[8]開發(fā)了一種基于可見-近紅外光譜的小型含水率監(jiān)測系統(tǒng), 采用反射光學(xué)探頭, 具有實現(xiàn)潤滑油狀態(tài)實時監(jiān)測的潛力。 采用偏最小二乘法(partial least squares, PLS)和反向傳播神經(jīng)網(wǎng)絡(luò)(back-propagation neural network, BPNN)算法建立處理模型。 采用馬爾滕斯不確定性檢驗、 區(qū)間偏最小二乘(iPLS)和遺傳算法(genetic algorithm, GA)變量選擇技術(shù)選擇最顯著的預(yù)測變量。 Holland等[9]采用傅里葉變換紅外光譜(Fourier transform infrared spectroscopy, FTIR)研究了SAE 15W-40柴油機潤滑油在不同污染水平的水分含量, 用方差分析建模和檢測極限計來預(yù)測區(qū)分污染水平隨時間變化的能力。 近紅外光譜的信息源是分子內(nèi)部原子間振動的倍頻與合頻, 不同物質(zhì)在近紅外區(qū)域有豐富的吸收光譜, 每種成分都有特定的吸收特征, 采用近紅外光譜對含水潤滑油光譜特征的分析也取得了一些進展[10], Holzki等[11]研發(fā)了一種用于潤滑油含水率原位測定的光纖倏變場吸收傳感器, 測量無包層光纖周圍消光場在約2 950 nm光譜范圍內(nèi)的吸光度。 并采用了介質(zhì)電泳技術(shù)來提高消失場區(qū)域內(nèi)水滴的濃度, 檢測靈敏度平均提高了1.7倍。 陳彬等[12]分析了含水汽輪機油的近紅外光譜, 在1 200, 1 400和1 700 nm波長附近有明顯的吸收峰存在, 其中1 400 nm附近是對應(yīng)O—H鍵的1倍頻; 結(jié)合連續(xù)投影算法(successive projections algorithm, SPA), 建立油中含水量預(yù)測的全波段偏最小二乘法(PLS)模型, 實現(xiàn)了油中含水量的定量分析。 并采用無信息變量消除法提取了光譜的有效波長, 選擇的波長主要集中在O—H鍵的合頻附近1 600~1 700 nm、 1倍頻1 400~1 600 nm和2倍頻900~1 000 nm處, 具有含水油液光譜的物理特征意義。

    因而擬采用近紅外光譜結(jié)合偏最小二乘回歸的方法檢測發(fā)動機油中水分, 在水分分析建模過程中, 基于油液近紅外光譜復(fù)雜、 非線性和難以用明確數(shù)學(xué)模型表達的特點, 采用正交信號校正和幾種傳統(tǒng)的光譜預(yù)處理方法去除光譜的基線偏移、 儀器誤差和測試噪聲。 為了保障所建模型的魯棒性和普適性, 根據(jù)回歸系數(shù)的重要程度, 進行了關(guān)鍵光譜波長的選擇。 將決定系數(shù)R2、 相對預(yù)測性能(relative prediction performance, RPP)、 均方根誤差(root mean squared error, RMSEP)作為模型優(yōu)劣評判的性能指標。 根據(jù)這些指標對所建模型進行比較分析, 獲取發(fā)動機油中水分的近紅外檢測的可靠有效的方法, 為發(fā)動機的安全高效運行提供理論研究和技術(shù)支撐。

    1 實驗部分

    1.1 儀器與油樣

    實驗采用海洋光學(xué)NIRQUEST近紅外光譜儀(型號512-1.9, 波長范圍900~1 700 nm, InGaAs檢測器), 選用SAE10W-30型發(fā)動機油(適用于外部環(huán)境溫度范圍為-25~30 ℃, 殼牌(中國)有限公司生產(chǎn))作為原始油樣, 用移液管將10 mL原始油樣放入70個直徑15 mm、 高度150 mm的玻璃試管中, 配置不同的含水比例(0~2.5%V/V)污染油樣。 配置方法是通過使用移液管用等量的蒸餾水代替適當體積的油, 如除去10 μL的機油, 添加10 μL的蒸餾水, 從而得到0.1%的水污染樣品。 分別用10~250 μL的蒸餾水得到0.1%~2.5%共34個污染級別的油樣, 每個污染級別油樣配置2個樣品和2個0%污染的原始油樣共70個樣品。 每個油樣采用超聲震蕩器進行充分攪拌2 h(50 Hz, 25 ℃), 以保證油水有效的混合。 混合均勻后的油樣采用831型庫侖水分儀(瑞士萬通)進行水分的精確檢測, 作為油樣的最終含水量, 并將它們儲存在一個封閉的試管中, 并直立放置在一個封閉的盒子中, 以防止光降解。

    1.2 光譜測試方法

    近紅外光譜儀是基于比爾-朗伯定律獲取油樣的光譜, 光譜的波長范圍為900~1 700 nm, 光學(xué)分辨率(FWHM)為3.1 nm, 波長重復(fù)性為0.1 nm。 采用光源為50 W石英鹵素?zé)簦?使用的檢測器類型為InGaAs, 色散原件為25 μm光柵, 所有樣品均在室溫下(25 ℃)加入1 cm比色皿中進行測試。 光譜采集軟件為OCEANVIEW, 掃描次數(shù)為20次。 化學(xué)計量學(xué)分析中使用的光譜為每個樣品的三次重復(fù)的平均值。

    然后將光譜儀采集的各個油樣光譜數(shù)據(jù)采用UnscramblerX化學(xué)計量學(xué)軟件進行光譜預(yù)處理和化學(xué)計量學(xué)分析, 采用MATLAB軟件對光譜數(shù)據(jù)根據(jù)波長回歸系數(shù)進行重要性分析和選擇。 光譜的預(yù)處理采用Savitzky-Golay(15點平滑窗, 二階多項式)、 標準正態(tài)變量變換(standard normal variable transformation, SNV)、 多元散射校正(multiplicative scatter correction, MSC)、 正交信號校正(orthogonal signal correction, OSC)以及使用(Savitzky-Golay, SG)算法的一階或二階導(dǎo)數(shù)等方法, 采用非線性迭代偏最小二乘(non-iterative PLS, NIPALS)算法構(gòu)建不同油樣和光譜的非線性關(guān)系, 以消除近紅外光譜中不相關(guān)或正交變量。 再進行主成分分析(principal component analysis, PCA)選擇并消除對回歸模型產(chǎn)生干擾的離群油樣。 并將70個油樣劃分為訓(xùn)練集和驗證集, 訓(xùn)練模型采用完整交叉驗證(留一法)進行構(gòu)建。 將不同含水量的油樣與相應(yīng)的光譜數(shù)據(jù)作為回歸模型的輸入輸出量, 同樣采用交叉驗證的最小均方根(RMSECV)確定回歸模型的最佳因子數(shù)。 采用預(yù)測決定系數(shù)(R2)、 預(yù)測均方根誤差(RMSEP)和RMSECV對所建立的PLS回歸模型進行評價。 采用殘差預(yù)測比(RPD)表示擬合優(yōu)度。 RPD是一種用于確定模型預(yù)測精度的無量綱統(tǒng)計量, 為數(shù)據(jù)集的標準偏差與RMSECV的比率。 一般地, RPD小于1.5表示校正模型不適合預(yù)測, 1.5~2表示有可能區(qū)分高低水平, 2.0~2.5表示有可能進行近似定量預(yù)測, 在2.5~3以及3以上表示預(yù)測效果是好的和較優(yōu)的模型[13]。

    最后為了降低回歸模型的運算量, 以及獲取光譜數(shù)據(jù)中的重要波長變量, 根據(jù)所建模型的回歸系數(shù)選擇不相關(guān)的波長, 降低模型的維數(shù), 重新建立重要波長變量的回歸模型, 并與全波長的回歸模型進行比較。

    2 結(jié)果與討論

    2.1 油樣光譜的特征以及異常值的處理

    配置的油樣經(jīng)過超聲震蕩器充分混合處理后, 油中的水分會發(fā)生微小變化, 因而采用庫侖水分儀進行精確測定含水量, 有助于在獲取分析發(fā)動機油中近紅外光譜數(shù)據(jù)后, 采用化學(xué)計量分析中回歸模型的響應(yīng)具有穩(wěn)健型和可靠性。 配置了經(jīng)過水分儀檢測后的含水量的油樣70個, 所配置的油樣可分為35個級別。 其中34個0.1%~2.5%水污染級別的油樣中, 每個污染級別有2個相同含水量級別的油樣, 共68個油樣。

    2.1.1 油樣光譜的特征

    所有油樣的原始光譜如圖1(a)所示, 在波長931, 1 195~1 212和1 391~1 430 nm處出現(xiàn)了較強的吸收峰, 對應(yīng)于O—H鍵伸縮振動頻率引起的2倍頻、 伸縮振動和扭轉(zhuǎn)振動的合頻、 伸縮振動1倍頻光譜[8]。 在1 391, 1 413和1 430 nm處觀察到三個峰值, 在1 195和1 212 nm處有兩個峰值, 在931 nm處有一個小峰。 這些峰值的強度隨含水量的增加而有明顯的增加。

    圖1 不同含水量油樣的原始光譜和預(yù)處理后光譜(a): 原始光譜; (b): SG一階導(dǎo)數(shù)處理后的光譜Fig.1 Original and pretreated spectra of oil samples with different water contents(a): Original spectrum; (b): Spectra after SG first derivative processing

    采用SG卷積求取光譜的一階導(dǎo)數(shù)如圖1(b)所示, 消除了原始光譜的基線漂移以及背景干擾的影響后, O—H鍵的伸縮振動和扭轉(zhuǎn)振動的合頻、 伸縮振動1倍頻光譜引起峰值的導(dǎo)數(shù)幅值非常明顯, 可見采用近紅外光譜分析發(fā)動機油中水分具有較強的分辨能力。 由圖1的原始光譜和預(yù)處理后光譜可見, 有兩條光譜曲線的變化明顯與其他光譜曲線存在差異, 可能是異常油樣, 在建立回歸模型前需要進行油樣異常值的剔除。

    2.1.2 異常值的剔除

    由于油樣在近紅外測試過程中會存在漏光、 光源強度發(fā)生變化以及異常油樣干擾等因素, 影響光譜數(shù)據(jù)的異常, 將對發(fā)動機油的含水量回歸模型的預(yù)測性能、 魯棒性與適用性產(chǎn)生較大地影響。 因而在建立回歸模型前需要對異常光譜進行檢測和剔除。

    首先對所有油樣進行PCA分析, 得到圖2(a)所示的兩個主成分分數(shù)圖, 可以看到隨著油中水分的增加, 油樣從左到右分布在主成分PC1的兩側(cè), 即含水量相近的油樣對兩個主成分的相關(guān)性相似, 低含水量的油樣(<0.30%)和高含水量(>0.80%)在兩個主成分的相關(guān)性方面具有較大的差異。 表明配置的油樣用于預(yù)測具有有效性。 然后用置信度為97%的Hotelling’s T2統(tǒng)計橢圓檢測油樣, 發(fā)現(xiàn)有1.00%和2.50%兩個污染級別的4個油樣位于橢圓外, Hotelling’s T2統(tǒng)計量是通過測量在主成分方向上的變化到模型中心的距離來反應(yīng)油樣在主成分子空間中偏離模型的程度, 說明這兩個污染級別的油樣是異常值(圖中以圓圈表示)。

    圖2 原始光譜PCA分析的分數(shù)圖(a)和影響圖(b)Fig.2 Fraction diagram (a) and influence Diagram (b) of PCA analysis of original spectrum

    在PCA分析的影響圖[圖2(b)]中, 表示的是油樣的F殘差和Hotelling’s T2統(tǒng)計量的控制限情況, 其中F殘差(置信度為97%)表示油樣在殘差子空間中偏離主成分模型的程度。 可見1.00%和2.50%兩個污染級別的4個油樣位于兩個控制限的右上方, 表示這4個油樣的F殘差較大、 且Hotelling’s T2統(tǒng)計量也較大, 是需要剔除的異常值。 另外在圖2(b)F殘差控制限的左上方出現(xiàn)了0.80%, 0.90%和0.33%三個污染級別的6個油樣, 表明這三個級別的油樣F殘差較大, 在主成分方向上的變化到模型中心的距離較小。 在實驗中基于發(fā)動機油中水分含量一般情況下較小, 因而為了所建模型的魯棒性以及預(yù)測的準確性, 將這三個污染級別的油樣不作為異常值剔除(圖中以圓圈表示)。

    剔除1.00%和2.50%兩個污染等級的油樣異常值后, 剩余的66個油樣中52個油樣用于構(gòu)建回歸模型的訓(xùn)練集, 剩余的14個油樣作為回歸模型的驗證集。

    2.2 構(gòu)建回歸模型

    偏最小二乘回歸(partial least squares regression, PLSR)結(jié)合了多元線性回歸、 典型相關(guān)分析和主成分分析的特點, 是最常用的化學(xué)計量學(xué)方法, 在光譜建模分析中得到廣泛應(yīng)用。 將剔除異常值后不同含水量的油樣近紅外光譜數(shù)據(jù)與含水量構(gòu)建 PLSR進行光譜分析。 PLSR數(shù)據(jù)矩陣X(n,m)由光譜儀器采集的油樣的近紅外光譜組成, 響應(yīng)變量y(n, 1)是水分含量的值, 其中n為訓(xùn)練數(shù)據(jù)集中的樣本數(shù)52個;m為近紅外光譜的波長數(shù)為512個。

    圖3 原始光譜的預(yù)測值和實測值的相關(guān)圖Fig.3 Correlation between predicted and measured values for the original spectrum

    2.2.1 光譜的預(yù)處理

    原始光譜的光散射、 噪聲和基線漂移等影響因素對原始近紅外光譜的疊加或相乘效應(yīng), 致使一些不相關(guān)的信息添加到構(gòu)建的預(yù)測模型中, 減弱了由于樣品的物理特性造成的一些差異, 因而在構(gòu)建回歸模型前需要采用不同的光譜預(yù)處理方法對光譜的數(shù)據(jù)進行處理, 以消除干擾信息, 提高模型的穩(wěn)定性。

    表1 光譜數(shù)據(jù)不同預(yù)處理后的PLS結(jié)果Table 1 PLS modeling results of spectral data after different pretreatments

    圖4 OSC預(yù)處理后全光譜的預(yù)測值和實測值的相關(guān)圖Fig.4 Correlation between predicted and measured values for full spectrum after OSC pretreatment

    2.2.2 光譜特征波長的選擇

    雖然全近紅外光譜經(jīng)過OSC預(yù)處理后的模型預(yù)測效果較好, 但是模型中包含了大量的不相關(guān)信息, 這會增加模型變量的數(shù)量和復(fù)雜性, 同時降低PLS模型的泛化性能。 因此, 很有必要對建模中的近紅外光譜波長變量進行選擇, 提取特征波長, 代表原始光譜信息中有用的部分, 有效地消除冗余變量, 進一步簡化模型, 從而提高PLS模型的穩(wěn)健性。

    根據(jù)PLS模型回歸系數(shù)是回歸模型光譜波長變量的權(quán)重系數(shù), 表征了所建模型的x變量和給定y響應(yīng)之間的關(guān)系。 通過將歸一化的回歸系數(shù)與回歸模型的預(yù)測值進行復(fù)雜的數(shù)學(xué)處理, 所得的系數(shù)矩陣中大于等于1的系數(shù)對應(yīng)的波長為模型影響較大的變量, 以此就可以對OSC預(yù)處理光譜進行特征波長的選擇。

    2.3 油樣的預(yù)測

    根據(jù)所建的OSC預(yù)處理后的全譜PLS模型以及特征波長選擇的PLS模型, 對油樣的預(yù)測集進行驗證, 得到如圖5所示的特征波長選擇前后的預(yù)測結(jié)果比較。 可見基于近紅外全譜的PLS模型以及特征波長選擇后的PLS模型都能較好地對訓(xùn)練集以外的油樣預(yù)測。 其中全譜PLS模型的預(yù)測集中14個油樣的預(yù)測的最大標準差為0.000 7, 主要是0.9%和2%兩個油樣偏離引起的; 通過特征波長選擇后的PLS模型對14個油樣的預(yù)測的最大標準差為0.000 6, 也是0.9%和2%兩個油樣引起; 這主要是由于兩個PLS模型的訓(xùn)練集中的油樣含水量主要分布在較低水平上, 含水量較大的油樣訓(xùn)練樣本較少, 從而造成PLS模型對較高含水量的預(yù)測能力下降。 另外由圖5可見, 特征波長選擇后的PLS模型對預(yù)測集的預(yù)測效果較優(yōu), 每個油樣的預(yù)測值更接近實測值。 說明經(jīng)過特征波長選擇后建立的PLS模型不僅沒有降低模型的精度和預(yù)測能力, 反而由于消除了不相關(guān)變量的信息, 更具有泛化性能。

    圖5 特征波長選擇前后的預(yù)測比較Fig.5 Comparison of predicted values before and after characteristic wavelength selection

    3 結(jié) 論

    猜你喜歡
    油樣波長含水量
    模擬工況條件下改性烴類航空潤滑油綜合性能對比分析*
    潤滑與密封(2023年7期)2023-07-28 03:28:24
    HPLC-PDA雙波長法同時測定四季草片中沒食子酸和槲皮苷的含量
    基礎(chǔ)油對發(fā)動機油TEOST 33C沉積物的影響
    蠟碳數(shù)分布和蠟含量對結(jié)蠟層分層的影響研究
    無硫磷鉬酸酯作為抗氧添加劑與胺類抗氧劑的協(xié)同性能
    結(jié)合Sentinel- 1B和Landsat8數(shù)據(jù)的針葉林葉片含水量反演研究
    森林工程(2018年4期)2018-08-04 03:23:16
    一次冰雹過程成雹機理的數(shù)值模擬
    雙波長激光治療慢性牙周炎的療效觀察
    日本研發(fā)出可完全覆蓋可見光波長的LED光源
    中國照明(2016年4期)2016-05-17 06:16:15
    基于兩相混合流理論PEMFC含水量特性分析
    亚洲一码二码三码区别大吗| 久久热在线av| 中文字幕av电影在线播放| 国产精品一区二区精品视频观看| 久久久国产欧美日韩av| 黄色一级大片看看| 日韩一本色道免费dvd| 日本wwww免费看| 久久国产精品人妻蜜桃| 久久国产精品影院| av欧美777| 亚洲成色77777| 超色免费av| 爱豆传媒免费全集在线观看| 午夜91福利影院| 超碰97精品在线观看| 久久久国产欧美日韩av| 精品国产乱码久久久久久男人| 成人黄色视频免费在线看| 亚洲精品国产一区二区精华液| 中文字幕高清在线视频| 国产精品久久久久久精品古装| 国产淫语在线视频| 欧美日韩综合久久久久久| 麻豆国产av国片精品| 日韩av免费高清视频| 天天躁日日躁夜夜躁夜夜| 久久久久久人人人人人| www.999成人在线观看| 母亲3免费完整高清在线观看| 18禁黄网站禁片午夜丰满| 欧美久久黑人一区二区| 国产精品一区二区在线观看99| 一级,二级,三级黄色视频| 午夜激情久久久久久久| a级毛片黄视频| 无遮挡黄片免费观看| 日韩熟女老妇一区二区性免费视频| 天天躁夜夜躁狠狠久久av| 啦啦啦在线观看免费高清www| 一本久久精品| 免费少妇av软件| 中文字幕另类日韩欧美亚洲嫩草| 777米奇影视久久| 我要看黄色一级片免费的| 国产野战对白在线观看| 大陆偷拍与自拍| 一区二区av电影网| 啦啦啦视频在线资源免费观看| 飞空精品影院首页| 岛国毛片在线播放| 国产一区二区在线观看av| 美女福利国产在线| av在线app专区| 久久国产精品大桥未久av| 精品久久久久久电影网| 五月天丁香电影| 国产在线观看jvid| 日韩av不卡免费在线播放| 国产精品免费大片| 亚洲中文字幕日韩| 国产三级黄色录像| 久久久久视频综合| 亚洲精品成人av观看孕妇| 久热这里只有精品99| 亚洲黑人精品在线| 久久久国产一区二区| 高潮久久久久久久久久久不卡| 99re6热这里在线精品视频| 午夜免费男女啪啪视频观看| 日韩av在线免费看完整版不卡| 2018国产大陆天天弄谢| 国产精品一二三区在线看| 国产伦理片在线播放av一区| av福利片在线| 无限看片的www在线观看| 欧美av亚洲av综合av国产av| 亚洲精品自拍成人| 人人妻人人添人人爽欧美一区卜| 99久久人妻综合| 色网站视频免费| 99热网站在线观看| 最近手机中文字幕大全| 另类亚洲欧美激情| 18禁黄网站禁片午夜丰满| netflix在线观看网站| 国产亚洲精品第一综合不卡| a级毛片在线看网站| 国产精品三级大全| 建设人人有责人人尽责人人享有的| 亚洲视频免费观看视频| 国产熟女午夜一区二区三区| 精品久久久精品久久久| 精品人妻在线不人妻| 1024视频免费在线观看| 高清不卡的av网站| 美女午夜性视频免费| 亚洲国产欧美一区二区综合| 中文字幕最新亚洲高清| 91字幕亚洲| 成人亚洲欧美一区二区av| 操美女的视频在线观看| 日本欧美国产在线视频| 另类亚洲欧美激情| 一区二区av电影网| 色婷婷久久久亚洲欧美| 欧美变态另类bdsm刘玥| 最新的欧美精品一区二区| 日韩av免费高清视频| 中文字幕色久视频| 久久久久网色| 精品少妇黑人巨大在线播放| 桃花免费在线播放| 色婷婷av一区二区三区视频| 亚洲欧美一区二区三区国产| 国产欧美亚洲国产| 欧美人与性动交α欧美精品济南到| 亚洲成人手机| 精品亚洲乱码少妇综合久久| 波野结衣二区三区在线| 久久99精品国语久久久| www.精华液| 久久久精品区二区三区| 国产成人一区二区在线| 国产成人免费观看mmmm| 人妻人人澡人人爽人人| 国产不卡av网站在线观看| 国产成人精品久久久久久| 亚洲国产中文字幕在线视频| 两个人看的免费小视频| 少妇粗大呻吟视频| 日韩制服丝袜自拍偷拍| 大陆偷拍与自拍| 欧美日韩福利视频一区二区| 亚洲精品国产av成人精品| 在线亚洲精品国产二区图片欧美| av天堂在线播放| 久久女婷五月综合色啪小说| 亚洲国产欧美一区二区综合| 多毛熟女@视频| 久久久久网色| 国产高清视频在线播放一区 | 一级片'在线观看视频| 在线观看免费高清a一片| 妹子高潮喷水视频| 亚洲成色77777| 国产国语露脸激情在线看| 色视频在线一区二区三区| 狂野欧美激情性bbbbbb| 国产精品久久久人人做人人爽| 一本一本久久a久久精品综合妖精| 国产色视频综合| 国产精品一区二区在线不卡| 韩国高清视频一区二区三区| 91精品三级在线观看| 国产免费福利视频在线观看| 一级毛片女人18水好多 | 欧美xxⅹ黑人| 成人18禁高潮啪啪吃奶动态图| 中文字幕精品免费在线观看视频| 亚洲男人天堂网一区| 日韩一本色道免费dvd| 脱女人内裤的视频| 一区二区三区激情视频| 亚洲国产看品久久| 在线观看免费午夜福利视频| 亚洲一区中文字幕在线| 国产真人三级小视频在线观看| 免费看十八禁软件| 亚洲欧美日韩高清在线视频 | 日韩大片免费观看网站| 老司机亚洲免费影院| 精品一品国产午夜福利视频| 亚洲av欧美aⅴ国产| 狂野欧美激情性xxxx| 少妇 在线观看| 80岁老熟妇乱子伦牲交| av网站免费在线观看视频| 天堂8中文在线网| 一级黄片播放器| 精品人妻在线不人妻| 高清黄色对白视频在线免费看| 日韩制服骚丝袜av| 热99国产精品久久久久久7| 丁香六月天网| 满18在线观看网站| 久久久精品94久久精品| 精品国产超薄肉色丝袜足j| 宅男免费午夜| 国产精品 国内视频| 国产精品欧美亚洲77777| 国产精品99久久99久久久不卡| 999精品在线视频| 天堂中文最新版在线下载| 国产成人免费观看mmmm| 最近手机中文字幕大全| 亚洲欧美精品综合一区二区三区| 午夜激情av网站| 9191精品国产免费久久| 午夜福利影视在线免费观看| 亚洲精品在线美女| 热99久久久久精品小说推荐| 国产成人精品在线电影| 国产成人精品无人区| 欧美日韩成人在线一区二区| 免费av中文字幕在线| 黑人猛操日本美女一级片| 亚洲欧美中文字幕日韩二区| 悠悠久久av| 日本黄色日本黄色录像| 91麻豆av在线| av天堂在线播放| 人体艺术视频欧美日本| 亚洲av欧美aⅴ国产| 国产高清国产精品国产三级| 国产一区有黄有色的免费视频| 香蕉丝袜av| 亚洲精品国产色婷婷电影| 韩国高清视频一区二区三区| 啦啦啦 在线观看视频| 日韩av免费高清视频| 777久久人妻少妇嫩草av网站| 国产在线视频一区二区| 黄色片一级片一级黄色片| 一二三四在线观看免费中文在| 一区二区日韩欧美中文字幕| 中文欧美无线码| 制服人妻中文乱码| 欧美日韩av久久| 巨乳人妻的诱惑在线观看| 成人免费观看视频高清| 七月丁香在线播放| 中文字幕精品免费在线观看视频| 亚洲av日韩精品久久久久久密 | 99国产精品免费福利视频| 大香蕉久久成人网| 啦啦啦视频在线资源免费观看| 悠悠久久av| 母亲3免费完整高清在线观看| 国产精品欧美亚洲77777| 一级a爱视频在线免费观看| 久久精品aⅴ一区二区三区四区| 亚洲一区二区三区欧美精品| 亚洲欧美清纯卡通| av欧美777| 亚洲激情五月婷婷啪啪| 青春草亚洲视频在线观看| 中文字幕亚洲精品专区| 国产精品av久久久久免费| 七月丁香在线播放| 中文乱码字字幕精品一区二区三区| 日韩电影二区| 麻豆国产av国片精品| 97人妻天天添夜夜摸| 男人添女人高潮全过程视频| 亚洲av成人精品一二三区| 在线观看免费高清a一片| kizo精华| 国产极品粉嫩免费观看在线| 热re99久久精品国产66热6| 久久亚洲国产成人精品v| xxx大片免费视频| 美女视频免费永久观看网站| 久久人妻福利社区极品人妻图片 | 国产又爽黄色视频| 婷婷色综合大香蕉| 人人澡人人妻人| 我要看黄色一级片免费的| 欧美97在线视频| 日本一区二区免费在线视频| 免费日韩欧美在线观看| 成在线人永久免费视频| 男女边吃奶边做爰视频| 久久精品熟女亚洲av麻豆精品| 又大又爽又粗| av电影中文网址| 男人添女人高潮全过程视频| 国产一区二区在线观看av| 欧美精品啪啪一区二区三区 | 2021少妇久久久久久久久久久| 亚洲成av片中文字幕在线观看| 国产熟女午夜一区二区三区| 欧美乱码精品一区二区三区| kizo精华| 久久ye,这里只有精品| 少妇粗大呻吟视频| 国产伦人伦偷精品视频| 欧美大码av| 成年动漫av网址| 麻豆av在线久日| 国产精品一区二区免费欧美 | 高清不卡的av网站| 亚洲第一av免费看| 日本wwww免费看| 日本欧美国产在线视频| a级片在线免费高清观看视频| 国产精品成人在线| 天堂8中文在线网| 国产一级毛片在线| 性色av乱码一区二区三区2| 丁香六月欧美| 精品一区二区三卡| 久久九九热精品免费| 国精品久久久久久国模美| 久久av网站| 亚洲天堂av无毛| 美女福利国产在线| 999久久久国产精品视频| avwww免费| 亚洲天堂av无毛| 一本大道久久a久久精品| 一个人免费看片子| 久久免费观看电影| 9色porny在线观看| www.自偷自拍.com| 大型av网站在线播放| 人人澡人人妻人| 国产黄色视频一区二区在线观看| 一本一本久久a久久精品综合妖精| 精品一区二区三区四区五区乱码 | 亚洲精品中文字幕在线视频| 亚洲国产精品一区二区三区在线| 国产精品久久久久成人av| 18禁国产床啪视频网站| 精品一区二区三区四区五区乱码 | 十八禁高潮呻吟视频| 婷婷色综合大香蕉| 女人爽到高潮嗷嗷叫在线视频| 天堂俺去俺来也www色官网| 18禁观看日本| 国产av一区二区精品久久| 国产福利在线免费观看视频| 七月丁香在线播放| 欧美日韩国产mv在线观看视频| 欧美成人午夜精品| 欧美97在线视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产av影院在线观看| 亚洲,一卡二卡三卡| 亚洲熟女毛片儿| kizo精华| 这个男人来自地球电影免费观看| 欧美精品av麻豆av| 欧美日韩亚洲高清精品| 亚洲欧美一区二区三区黑人| 9色porny在线观看| 欧美精品高潮呻吟av久久| 亚洲免费av在线视频| 九色亚洲精品在线播放| 精品久久久久久久毛片微露脸 | 欧美亚洲 丝袜 人妻 在线| 中文字幕av电影在线播放| 亚洲国产欧美日韩在线播放| 叶爱在线成人免费视频播放| 国产精品 欧美亚洲| 91字幕亚洲| 色精品久久人妻99蜜桃| 国产1区2区3区精品| 性色av一级| 国产精品成人在线| 亚洲午夜精品一区,二区,三区| 亚洲av美国av| 2018国产大陆天天弄谢| 日韩 欧美 亚洲 中文字幕| 最黄视频免费看| 欧美 亚洲 国产 日韩一| 久久久欧美国产精品| 亚洲国产中文字幕在线视频| 日本五十路高清| 天天躁夜夜躁狠狠久久av| 亚洲色图综合在线观看| 岛国毛片在线播放| 一级毛片电影观看| 日本五十路高清| 成年美女黄网站色视频大全免费| 欧美性长视频在线观看| www.999成人在线观看| 午夜激情久久久久久久| 成年av动漫网址| 人人妻人人澡人人看| 国产在视频线精品| 激情视频va一区二区三区| 欧美日韩亚洲综合一区二区三区_| 欧美大码av| 91成人精品电影| 99精国产麻豆久久婷婷| 一本大道久久a久久精品| 亚洲,欧美,日韩| 女人精品久久久久毛片| 久热爱精品视频在线9| 日韩一本色道免费dvd| 欧美日韩国产mv在线观看视频| 久久久精品免费免费高清| 人人妻人人澡人人看| 亚洲伊人色综图| 又大又爽又粗| 一边摸一边抽搐一进一出视频| 伊人亚洲综合成人网| 美女福利国产在线| 欧美精品高潮呻吟av久久| 亚洲综合色网址| 日韩 欧美 亚洲 中文字幕| 国产成人av教育| 日韩制服丝袜自拍偷拍| 国产精品一国产av| 国产一级毛片在线| 国产色视频综合| 国产高清视频在线播放一区 | 中文字幕另类日韩欧美亚洲嫩草| 日韩免费高清中文字幕av| 国产伦理片在线播放av一区| 精品少妇久久久久久888优播| 久久亚洲精品不卡| 亚洲av电影在线观看一区二区三区| 日韩电影二区| 国产高清国产精品国产三级| 午夜福利乱码中文字幕| 在线观看国产h片| 国产野战对白在线观看| 午夜激情久久久久久久| 国产成人系列免费观看| 又粗又硬又长又爽又黄的视频| 久久精品亚洲熟妇少妇任你| 成在线人永久免费视频| 欧美国产精品va在线观看不卡| 国产精品九九99| 黄片播放在线免费| 老汉色∧v一级毛片| 男人操女人黄网站| 女性被躁到高潮视频| 国产男女超爽视频在线观看| 国产精品一区二区在线不卡| 精品一品国产午夜福利视频| av又黄又爽大尺度在线免费看| 天天影视国产精品| 亚洲国产av新网站| 婷婷色av中文字幕| 电影成人av| 欧美精品亚洲一区二区| 最近中文字幕2019免费版| 亚洲欧美日韩高清在线视频 | 老司机在亚洲福利影院| 国产亚洲欧美在线一区二区| 久久热在线av| 久热这里只有精品99| 两个人看的免费小视频| 国产精品免费大片| 亚洲,欧美,日韩| 亚洲欧洲日产国产| 亚洲男人天堂网一区| 2018国产大陆天天弄谢| 十八禁高潮呻吟视频| 午夜激情久久久久久久| 国产在线观看jvid| 精品久久久久久电影网| 国产高清视频在线播放一区 | 日韩伦理黄色片| 97精品久久久久久久久久精品| 在线看a的网站| 麻豆乱淫一区二区| 国产成人精品久久二区二区91| 午夜福利视频精品| 老汉色av国产亚洲站长工具| 咕卡用的链子| 久久久国产精品麻豆| 日本vs欧美在线观看视频| 国产在线观看jvid| av片东京热男人的天堂| 国产爽快片一区二区三区| 亚洲免费av在线视频| 日韩制服丝袜自拍偷拍| 亚洲 国产 在线| a级毛片在线看网站| 色婷婷久久久亚洲欧美| 51午夜福利影视在线观看| 日韩 亚洲 欧美在线| 天堂中文最新版在线下载| 真人做人爱边吃奶动态| 美女高潮到喷水免费观看| 一级毛片我不卡| 五月开心婷婷网| 人人澡人人妻人| 80岁老熟妇乱子伦牲交| 啦啦啦在线免费观看视频4| 黑人巨大精品欧美一区二区蜜桃| 国产欧美亚洲国产| 国产97色在线日韩免费| 一边摸一边抽搐一进一出视频| 亚洲,欧美精品.| 国产国语露脸激情在线看| 欧美国产精品va在线观看不卡| 午夜免费男女啪啪视频观看| 亚洲人成电影免费在线| 手机成人av网站| 天天躁狠狠躁夜夜躁狠狠躁| 黑丝袜美女国产一区| 日本vs欧美在线观看视频| 曰老女人黄片| 久久天堂一区二区三区四区| 免费观看a级毛片全部| 亚洲欧美清纯卡通| 日本黄色日本黄色录像| 精品一区二区三区av网在线观看 | 你懂的网址亚洲精品在线观看| 中文精品一卡2卡3卡4更新| 考比视频在线观看| 激情五月婷婷亚洲| 国产又爽黄色视频| 99香蕉大伊视频| av线在线观看网站| 亚洲伊人久久精品综合| 久久久欧美国产精品| 国产精品三级大全| 看十八女毛片水多多多| 国产黄色视频一区二区在线观看| 国产免费福利视频在线观看| 美女视频免费永久观看网站| xxx大片免费视频| 国产亚洲av片在线观看秒播厂| 国产精品 欧美亚洲| 人妻 亚洲 视频| 啦啦啦中文免费视频观看日本| 亚洲精品美女久久av网站| 国产精品熟女久久久久浪| 一级毛片黄色毛片免费观看视频| 丰满饥渴人妻一区二区三| 婷婷色麻豆天堂久久| av天堂在线播放| 免费不卡黄色视频| 亚洲成人国产一区在线观看 | 建设人人有责人人尽责人人享有的| 交换朋友夫妻互换小说| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 亚洲国产成人一精品久久久| 日本a在线网址| 老司机靠b影院| av福利片在线| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美精品永久| 日本黄色日本黄色录像| 丁香六月欧美| 交换朋友夫妻互换小说| 国产视频一区二区在线看| 久久久精品国产亚洲av高清涩受| av在线app专区| 久久热在线av| 亚洲精品美女久久av网站| 亚洲一区中文字幕在线| 大型av网站在线播放| 三上悠亚av全集在线观看| 国产精品一国产av| 日韩电影二区| 各种免费的搞黄视频| 人人妻人人爽人人添夜夜欢视频| 婷婷丁香在线五月| 日本黄色日本黄色录像| 女性被躁到高潮视频| 国产精品久久久久久人妻精品电影 | 在线观看一区二区三区激情| 一级片'在线观看视频| 国产三级黄色录像| 国精品久久久久久国模美| 老鸭窝网址在线观看| 亚洲国产精品成人久久小说| 亚洲欧洲国产日韩| 欧美精品亚洲一区二区| 99热网站在线观看| 午夜免费成人在线视频| 又粗又硬又长又爽又黄的视频| 伊人久久大香线蕉亚洲五| 高清av免费在线| 中文字幕人妻丝袜一区二区| 亚洲精品久久午夜乱码| 国产成人精品久久二区二区91| 曰老女人黄片| 精品久久久久久电影网| 亚洲欧美一区二区三区国产| 国产成人免费观看mmmm| 男女床上黄色一级片免费看| 久久 成人 亚洲| 国产精品 国内视频| 国产精品熟女久久久久浪| 国产欧美日韩综合在线一区二区| h视频一区二区三区| 男女国产视频网站| 肉色欧美久久久久久久蜜桃| 老司机亚洲免费影院| 男人操女人黄网站| 日韩制服丝袜自拍偷拍| 两个人看的免费小视频| 亚洲,欧美,日韩| 老司机在亚洲福利影院| 操美女的视频在线观看| 男人添女人高潮全过程视频| 色播在线永久视频| 欧美xxⅹ黑人| 嫩草影视91久久| 日本av免费视频播放| 在线观看免费视频网站a站| 咕卡用的链子| 纵有疾风起免费观看全集完整版| 亚洲 国产 在线| 国产成人免费观看mmmm| 精品亚洲乱码少妇综合久久| 日本猛色少妇xxxxx猛交久久| 青春草视频在线免费观看| 国产免费一区二区三区四区乱码| 亚洲七黄色美女视频| 999久久久国产精品视频| 18禁裸乳无遮挡动漫免费视频| 国产又爽黄色视频| 精品卡一卡二卡四卡免费| 黄频高清免费视频| 久久精品久久久久久久性|