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

    基于近紅外光譜的沼液揮發(fā)性脂肪酸含量快速檢測

    2020-12-03 02:24:46劉金明郭坤林張鴻瓊李文哲許永花
    農(nóng)業(yè)工程學(xué)報 2020年18期
    關(guān)鍵詞:厭氧發(fā)酵總酸丙酸

    劉金明,郭坤林,甄 峰,張鴻瓊,李文哲,許永花

    ·農(nóng)業(yè)生物環(huán)境與能源工程·

    基于近紅外光譜的沼液揮發(fā)性脂肪酸含量快速檢測

    劉金明1,2,3,郭坤林2,甄 峰1,3,張鴻瓊1,4,李文哲1,4,許永花5※

    (1. 東北農(nóng)業(yè)大學(xué)工程學(xué)院,哈爾濱 150030;2. 黑龍江八一農(nóng)墾大學(xué)電氣與信息學(xué)院,大慶 163319;3. 中國科學(xué)院可再生能源重點(diǎn)實驗室,廣州 510640;4. 黑龍江省寒地農(nóng)業(yè)可再生資源利用技術(shù)及裝備重點(diǎn)實驗室,哈爾濱 150030;5. 東北農(nóng)業(yè)大學(xué)電氣與信息學(xué)院,哈爾濱 150030)

    揮發(fā)性脂肪酸(Volatile Fatty Acids,VFA)作為厭氧發(fā)酵過程的重要中間產(chǎn)物,其在厭氧反應(yīng)器中的累積能夠反映出產(chǎn)甲烷菌的不活躍狀態(tài)或厭氧發(fā)酵條件的惡化。為了實現(xiàn)對農(nóng)牧廢棄物厭氧發(fā)酵進(jìn)行過程分析和狀態(tài)監(jiān)控,將近紅外光譜(Near Infrared Spectroscopy,NIRS)與偏最小二乘(Partial Least Squares,PLS)相結(jié)合構(gòu)建玉米秸稈和畜禽糞便厭氧發(fā)酵液乙酸、丙酸和總酸含量快速檢測模型。將競爭自適應(yīng)重加權(quán)采樣法(Competitive Adaptive Reweighted Sampling,CARS)與遺傳模擬退火(Genetic Simulated Annealing,GSA)算法相結(jié)合構(gòu)建CARS-GSA算法對沼液中的乙酸、丙酸和總酸進(jìn)行特征波長優(yōu)選,原始光譜數(shù)據(jù)1 557個波長點(diǎn)經(jīng)預(yù)處理和波長優(yōu)選后,得到乙酸、丙酸和總酸特征波長變量分別為135、101和245個,建立的回歸模型驗證決定系數(shù)分別為0.988、0.923和0.886,預(yù)測均方根誤差(Root Mean Squared Error of Prediction,RMSEP)分別為0.111、0.120和0.727,相對分析誤差分別為9.685、3.685和3.484,與全譜建模相比RMSEP分別減少了17.78%、15.49%和1.22%,能夠滿足農(nóng)牧廢棄物厭氧發(fā)酵過程發(fā)酵液中乙酸和丙酸含量的快速檢測需求,基本滿足總酸的檢測需求。結(jié)果表明,通過構(gòu)建CARS-GSA算法優(yōu)選乙酸、丙酸和總酸的敏感波長變量,參與建模的波長點(diǎn)數(shù)量顯著減少,有效降低了變量維度和模型復(fù)雜度,提升了回歸模型檢測精度和預(yù)測能力,為快速準(zhǔn)確檢測沼液VFA提供了新途徑。

    厭氧發(fā)酵;揮發(fā)性脂肪酸;快速檢測;近紅外光譜;偏最小二乘;遺傳模擬退火算法;競爭自適應(yīng)重加權(quán)采樣

    0 引 言

    揮發(fā)性脂肪酸(Volatile Fatty Acids,VFA)作為厭氧發(fā)酵過程的重要中間產(chǎn)物,為產(chǎn)甲烷階段提供了底物[1]。產(chǎn)甲烷菌主要利用VFA形成甲烷,只有少部分甲烷由二氧化碳和氫氣生成,但二氧化碳和氫氣生成甲烷時也經(jīng)過高分子有機(jī)物形成VFA的中間過程[2]。VFA在厭氧反應(yīng)器中的積累能反映出產(chǎn)甲烷菌的不活躍狀態(tài)或厭氧發(fā)酵條件的惡化,較高的VFA濃度對產(chǎn)甲烷菌有抑制作用,過高的VFA濃度甚至?xí)?dǎo)致厭氧發(fā)酵發(fā)生“酸敗”[3]。在反應(yīng)器運(yùn)行過程中,發(fā)酵液的VFA濃度常用作厭氧發(fā)酵過程的重要監(jiān)控指標(biāo)[4]。通過監(jiān)測發(fā)酵液中VFA的變化情況,可以很好地了解有機(jī)物的降解過程以及產(chǎn)甲烷菌的活性和系統(tǒng)的運(yùn)行情況[5]。為了對厭氧發(fā)酵狀態(tài)進(jìn)行有效監(jiān)控,有必要對發(fā)酵液的VFA含量進(jìn)行快速、準(zhǔn)確測定。

    傳統(tǒng)的VFA檢測方法主要有精餾法、高效液相色譜法、氣相色譜法和各種滴定技術(shù)[5-6]。但傳統(tǒng)檢測方法存在處理時間長、設(shè)備操作復(fù)雜等問題,難以滿足厭氧發(fā)酵過程中通過快速測定VFA實現(xiàn)發(fā)酵過程狀態(tài)監(jiān)測的需求。針對厭氧發(fā)酵過程狀態(tài)監(jiān)測對VFA快速檢測的需求,相關(guān)學(xué)者深入研究了快速滴定法[7]、新型色譜技術(shù)[8]、電化學(xué)傳感器[9]、生物傳感器[10]和光譜分析技術(shù)[11]在VFA快速檢測方面的應(yīng)用。光譜分析技術(shù)因其簡便、快捷、無損、低成本的優(yōu)勢,已在發(fā)酵液VFA檢測方面得到了廣泛應(yīng)用[12-13],其中以近紅外光譜(Near Infrared Spectroscopy,NIRS)定量分析技術(shù)的應(yīng)用最為廣泛[14-16]。在應(yīng)用NIRS對水體中的VFA含量進(jìn)行快速檢測方面,主要以光譜預(yù)處理方法和多元定量校正方法的研究為主,在VFA特征波長優(yōu)選方面尚需進(jìn)一步拓展,以消除不相關(guān)和非線性波長點(diǎn)對模型精度的影響。

    當(dāng)前,NIRS特征波長優(yōu)選方法正朝著多種特征波長優(yōu)選方法相結(jié)合的方向發(fā)展,將區(qū)間偏最小二乘法(Interval Partial Least Squares,iPLS)[17]、協(xié)同區(qū)間偏最小二乘法(Synergy iPLS,SiPLS)[18]、反向區(qū)間偏最小二乘法(Backward iPLS,BiPLS)[19]、連續(xù)投影算法[20]、競爭自適應(yīng)重加權(quán)采樣法(Competitive Adaptive Reweighted Sampling,CARS)[21]等特征波長優(yōu)選算法與遺傳算法(Genetic Algorithm,GA)[22]、模擬退火算法(Simulated Annealing Algorithm,SA)[23]、粒子群優(yōu)化算法[24]等智能優(yōu)化算法相結(jié)合進(jìn)行NIRS特征波長變量優(yōu)選已成為重要研究方向[25-27]。在NIRS定量分析方面,GA因其強(qiáng)大的特征波長優(yōu)選能力已得到廣泛應(yīng)用[28-29],但GA存在早熟問題且進(jìn)化后期搜索效率低。

    遺傳模擬退火(Genetic Simulated Annealing,GSA)算法是對GA的改進(jìn),通過結(jié)合SA的溫度參數(shù)設(shè)計適應(yīng)度函數(shù),引入Metropolis擾動解選擇復(fù)制策略,有效利用GA強(qiáng)大搜索能力的同時解決了GA的兩點(diǎn)不足,在NIRS特征波長優(yōu)選方面取得了較好的應(yīng)用效果[30]。GSA在與iPLS、SiPLS和BiPLS相結(jié)合進(jìn)行特征波長優(yōu)選方面的性能也顯著優(yōu)于GA[30-32],但在使用iPLS、SiPLS和BiPLS進(jìn)行特征譜區(qū)優(yōu)選時,難以避免譜區(qū)內(nèi)部存在冗余波長點(diǎn)。這些不相關(guān)和非線性的冗余波長點(diǎn)導(dǎo)致GSA編碼過長,嚴(yán)重影響了GSA特征波長點(diǎn)二次搜索的性能。

    因此,本文針對以秸稈和糞便為原料的厭氧發(fā)酵過程中,沼液VFA快速檢測的需求,提出基于NIRS建立VFA快速檢測模型,并將GSA與CARS相結(jié)合構(gòu)建CARS-GSA算法進(jìn)行VFA特征波長優(yōu)選,有效解決iPLS、SiPLS和BiPLS敏感波段初步定位過程中存在冗余波長點(diǎn)的問題,進(jìn)而獲取滿足實際檢測需求的有效特征波長變量,以提高VFA快速檢測模型的效率和精度。

    1 材料與方法

    1.1 樣品采集與制備

    試驗用玉米秸稈取自東北農(nóng)業(yè)大學(xué)校內(nèi)試驗田,豬糞取自哈爾濱市三元畜產(chǎn)實業(yè)公司,牛糞取自哈爾濱市宇峰奶牛養(yǎng)殖農(nóng)民專業(yè)合作社,接種物取自黑龍江省寒地農(nóng)業(yè)可再生資源利用技術(shù)及裝備重點(diǎn)實驗室。將采集的玉米秸稈自然風(fēng)干后一部分經(jīng)鍘草機(jī)切成10 mm的秸稈段備用,另一部分經(jīng)錘片式粉碎機(jī)(10 mm篩網(wǎng))粉碎成秸稈粉備用。分別以秸稈段、秸稈粉、牛糞、豬糞、秸稈粉豬糞混合物(按總固體(Total Solid,TS)比1∶1)為厭氧發(fā)酵原料,以實驗室500 L發(fā)酵罐常年馴化正常產(chǎn)氣的牛糞厭氧發(fā)酵液為接種物,進(jìn)行批式厭氧發(fā)酵試驗。秸稈、牛糞、豬糞和接種物的TS濃度分別為86.02%、26.62%、31.22%和4.76%,按TS接種比1∶1,調(diào)整厭氧發(fā)酵原料和接種物添加量,使5種原料對應(yīng)的發(fā)酵系統(tǒng)起始TS濃度分別為7%、6%、8%、7%和7%。在中溫(36±1)℃恒溫水浴槽中,分別采用5和10 L下口瓶作為反應(yīng)器,進(jìn)行2個批次的厭氧發(fā)酵試驗,有效發(fā)酵容積分別為3.5和7 L。試驗過程中每天定時對厭氧發(fā)酵反應(yīng)器進(jìn)行手搖攪拌2次,混勻料液的同時避免浮渣結(jié)殼。為了獲取有代表性的VFA濃度數(shù)據(jù)樣本,采集發(fā)酵液樣品主要在批式厭氧發(fā)酵前半程進(jìn)行。5 L發(fā)酵罐從裝樣后第2天開始,每天8:00采集發(fā)酵液樣品40 mL存放于3個15 mL離心管中,共計采樣16次。為防止料液TS濃度變高,對厭氧發(fā)酵過程產(chǎn)生不良影響,于第8天補(bǔ)水300 mL。10 L發(fā)酵罐從裝樣后第2天開始采樣,共計采樣15次,不需補(bǔ)水;共計采集與制備發(fā)酵液樣品155個,于-20℃冰箱冷凍保存。

    1.2 光譜數(shù)據(jù)采集

    發(fā)酵液冷凍樣品溶解后在冷凍離心機(jī)中以12 000 r/min離心10 min后,取上清液待測。使用Nicolet公司的Antaris II型傅里葉近紅外光譜儀對采集樣品進(jìn)行透射光譜掃描,光譜采集范圍4 000~10 000 cm-1(1 000~2 500 nm),分辨率為8.0 cm-1,樣品掃描32次,數(shù)據(jù)保存格式為lg(1/T),背景每小時掃描一次,裝樣方式為1 mm光程石英比色皿前置通道掃描。在保持室內(nèi)溫濕度基本穩(wěn)定的情況下,每個樣品裝樣3次,取3次掃描平均值作為樣品的原始光譜。原始光譜的波長數(shù)量為1 557個,數(shù)據(jù)點(diǎn)間距為3.86 cm-1,起始波數(shù)為10 001.03 cm-1,結(jié)束波數(shù)為3 999.64 cm-1。

    1.3 VFA測定

    使用安捷倫GC-6890N氣相色譜儀測定厭氧發(fā)酵過程中沼液的VFA濃度。采用外標(biāo)法建立VFA標(biāo)準(zhǔn)曲線,先制備乙酸、丙酸、丁酸、異丁酸和異戊酸的混合標(biāo)準(zhǔn)溶液,再使用去離子水稀釋至6種不同濃度,并測定不同濃度標(biāo)準(zhǔn)溶液各成分對應(yīng)的出峰時間和積分面積。將混合溶液的保留時間與單品的保留時間進(jìn)行比較,根據(jù)已知標(biāo)準(zhǔn)溶液中各物質(zhì)的濃度和積分時間繪制標(biāo)準(zhǔn)曲線。對溶解、離心并采集透射光譜數(shù)據(jù)后的厭氧發(fā)酵液樣品上清液進(jìn)行VFA含量測定。將其與25%偏磷酸溶液按體積比10∶1進(jìn)行混合,然后再以12 000 r/min離心10 min后取上清液,將上清液使用0.45m超濾膜過濾,取濾液進(jìn)行VFA濃度測定。

    1.4 波長優(yōu)選方法

    1.4.1 CARS算法

    CARS算法基于“適者生存”的原則,將蒙特卡洛采樣(Monte-Carlo Sampling,MCS)、指數(shù)衰減函數(shù)和自適應(yīng)加權(quán)采樣(Adaptive Reweighted Sampling,ARS)相結(jié)合獲取波長子集,基于偏最小二乘(Partial Least Squares,PLS)回歸系數(shù)絕對值的大小獲取一系列變量組合,并選擇交叉驗證均方根誤差(Root Mean Squared Error of Cross Validation,RMSECV)值最小的子集作為特征波長。CARS在迭代過程中引入MCS和ARS 2個隨機(jī)因素,難以保證每次優(yōu)選結(jié)果的一致性??梢圆捎枚啻芜\(yùn)行CARS算法,每次都選中的波長點(diǎn)代表著光譜數(shù)據(jù)中與待測目標(biāo)屬性相關(guān)性高的波長點(diǎn),選定這些多次都選中波長作為特征波長,能夠建立高性能的回歸模型。

    1.4.2 CARS-GSA算法

    CARS-GSA算法以CARS優(yōu)選后的特征波長為輸入,采用GSA算法對CARS優(yōu)選結(jié)果進(jìn)行再優(yōu)化,以剔除CARS優(yōu)選結(jié)果中相關(guān)性較差的波長點(diǎn),從而進(jìn)一步提高建模性能。CARS-GSA以CARS優(yōu)選后特征波長點(diǎn)數(shù)為碼長,以PLS回歸模型的折RMSECV為目標(biāo)函數(shù),按初始種群個數(shù)約為碼長的三分之一進(jìn)行二進(jìn)制編碼和種群初始化。“1”和“0”分別表示該波長點(diǎn)對應(yīng)的數(shù)據(jù)“是”、“否”選中參與運(yùn)算。在確定初始溫度、退溫操作,并計算適應(yīng)度函數(shù)值后,執(zhí)行多個輪次的GSA選擇、交叉、變異和Metropolis選擇復(fù)制進(jìn)化操作,完成NIRS特征波長點(diǎn)的優(yōu)選。多次執(zhí)行GSA算法對CARS優(yōu)選結(jié)果進(jìn)行再優(yōu)化,并選擇多次重復(fù)選中的波長點(diǎn)作為特征波長變量建立PLS回歸模型,能夠得到較高的回歸模型性能。

    1.5 回歸模型建立及評價

    本文算法包括光譜預(yù)處理、樣本集劃分、特征波長優(yōu)選及回歸模型構(gòu)建等全部在Matlab R2012b軟件平臺中實現(xiàn)。

    2 結(jié)果與分析

    2.1 采集數(shù)據(jù)分析

    在采用安捷倫GC-6890N氣相色譜儀測定155個發(fā)酵液樣本的VFA濃度時,得到81個乙酸濃度有效數(shù)據(jù)、78個丙酸濃度有效數(shù)據(jù)和87個總酸濃度有效數(shù)據(jù)(總酸濃度為乙酸、丙酸、丁酸、異丁酸和異戊酸質(zhì)量分?jǐn)?shù)之和)。對獲得的VFA樣本有效濃度數(shù)據(jù)進(jìn)行四分位數(shù)分析,并繪制箱線圖如圖1所示。

    圖1 樣本VFA濃度箱線圖

    由圖1可知,乙酸樣本在低濃度區(qū)域占比較大,丙酸樣本略微偏向低濃度區(qū)域,總酸樣本分布比較均勻。乙酸樣本偏離嚴(yán)重的原因在于厭氧發(fā)酵產(chǎn)乙酸、產(chǎn)甲烷平衡期產(chǎn)甲烷菌能夠及時將生成的乙酸轉(zhuǎn)化為甲烷和二氧化碳,進(jìn)而使平衡期階段(在整個發(fā)酵周期中時間占比較大)的乙酸濃度偏低。

    為消除光譜區(qū)域中平頂峰對建模結(jié)果的影響,先剔除原始光譜數(shù)據(jù)中波數(shù)4 933.02~5 295.57 cm-1的95個波長點(diǎn),再用剩余的1 462個有效波長點(diǎn)建立乙酸、丙酸和總酸回歸模型,并對不同光譜預(yù)處理方法下的回歸模型性能進(jìn)行評測。經(jīng)計算比較后確定乙酸濃度回歸模型采用的光譜預(yù)處理方法為MSC+SG,丙酸回歸模型采用的光譜預(yù)處理方法為SG+MSC,總酸回歸模型采用的光譜預(yù)處理方法為FD+SNV+SG。樣品原始光譜及預(yù)處理后的乙酸、丙酸和總酸光譜數(shù)據(jù)的平均光譜如圖2所示。

    圖2 樣品光譜數(shù)據(jù)

    對81個乙酸樣品的原始光譜依次進(jìn)行MSC和SG平滑處理后,使用SPXY法劃分為60個校正集樣本和21個驗證集樣本;對78個丙酸樣品的原始光譜數(shù)據(jù)依次進(jìn)行SG平滑和MSC處理后,使用SPXY法劃分為60個校正集樣本和18個驗證集樣本;對87個總酸樣品的原始光譜數(shù)據(jù)依次進(jìn)行FD、SNV和SG平滑處理后,使用SPXY法劃分為70個校正集樣本和17個驗證集樣本。乙酸、丙酸和總酸濃度值如表1所示。

    表1 樣品VFA濃度

    注:SD是Standard deviation的縮寫,NS是Number of sample的縮寫。

    Note: SD is short for standard deviation, NS is short for number of sample.

    2.2 特征波長優(yōu)選

    2.2.1 CARS特征波長優(yōu)選

    在使用CARS優(yōu)選乙酸回歸模型特征波長時,先執(zhí)行500輪次CARS算法,再按重復(fù)選中次數(shù)遞增的方式選取RMSEP最小時對應(yīng)的特征波長優(yōu)選結(jié)果作為CARS的特征波長(記為CARS500)。執(zhí)行500次CARS算法共得到乙酸特征波長383個以波數(shù)表示,下同,選中次數(shù)最多的特征波長波數(shù)為4 416.19 cm-1,對應(yīng)著乙酸-CH3基團(tuán)的組合頻,選中次數(shù)為457次。選中次數(shù)較多的特征波長點(diǎn)主要分布在4 000~4 600、4 750~4 930、5 300~5 500、5 750~6 050、6 750~7 100和7 500~7 800 cm-1區(qū)域。其中4 000~4 600 cm-1對應(yīng)著乙酸-CH3基團(tuán)的組合頻,4 750~4 930 cm-1對應(yīng)著C=O和-OH基團(tuán)的組合頻,5 300~5 500 cm-1對應(yīng)著-COOH基團(tuán)的一級倍頻,5 750~6 050 cm-1對應(yīng)著-CH3基團(tuán)的一級倍頻,6 750~7 100 cm-1對應(yīng)著C=O和-OH基團(tuán)的二級倍頻,7 500~7 800 cm-1對應(yīng)著-CH3基團(tuán)的二級倍頻。CARS500優(yōu)選特征波長與乙酸平均光譜如圖3所示。

    圖3 CARS500優(yōu)選乙酸特征波長

    為分析不同重復(fù)選中次數(shù)下,CARS500優(yōu)選特征波長的建模性能,建立RMSECV、RMSEP和波長點(diǎn)個數(shù)隨重復(fù)選中次數(shù)的變化關(guān)系,如圖4所示。

    由圖4可知,RMSECV隨著選中波長點(diǎn)個數(shù)的減少整體上呈先迅速減少、再波浪狀向前、最后跳躍式快速上升的形式,其中波長點(diǎn)數(shù)為120時,RMSECV得到最小值0.163,對應(yīng)重復(fù)選中次數(shù)為39次。RMSEP隨選中波長點(diǎn)個數(shù)減少整體呈鋸齒型變化并逐漸增加的形式,其中重復(fù)選中次數(shù)為30、選中波長數(shù)量為142時,所建PLS回歸模型的RMSEP獲得最小值為0.116。采用RMSEP最小時對應(yīng)的142個波長點(diǎn)作為CARS500優(yōu)選的乙酸特征波長。

    圖4 RMSE、波長數(shù)量和選中次數(shù)間的關(guān)系

    2.2.2 CARS-GSA特征波長優(yōu)選

    在使用CARS-GSA優(yōu)選發(fā)酵液乙酸特征波長時,以CARS500優(yōu)選的142個波長點(diǎn)為碼長隨機(jī)生成50個染色體構(gòu)建初始種群,執(zhí)行GSA算法進(jìn)行特征波長點(diǎn)二次優(yōu)選。GSA算法的初溫確定系數(shù)取100,退溫系數(shù)取0.9,進(jìn)化代數(shù)取100,交叉概率取0.7,變異概率取0.01,鄰域解擾動位數(shù)取10。連續(xù)執(zhí)行算法50次,優(yōu)選的乙酸特征波長中選中35次以上的波長共計14個。其中,波數(shù)4 057.49、4 319.77、4 354.48、4 358.33、4 362.19、4 366.05、4 408.48、4 412.33、4 416.19、4 420.05、4 531.90和4 539.61 cm-1對應(yīng)著-CH3基團(tuán)的組合頻,波數(shù)4 925.30 cm-1對應(yīng)著C=O基團(tuán)的組合頻,波數(shù)5 311.00 cm-1對應(yīng)著-COOH基團(tuán)的一級倍頻。CARS-GSA優(yōu)選特征波長與乙酸平均光譜如圖5所示。

    圖5 CARS-GSA優(yōu)選乙酸特征波長

    為分析CARS-GSA優(yōu)選特征波長的建模性能,建立RMSECV、RMSEP與波長點(diǎn)個數(shù)間的對應(yīng)關(guān)系,如圖6所示。由圖6可知,RMSECV和RMSEP隨選中波長點(diǎn)個數(shù)增加整體上呈先迅速減少、再趨于平緩、最后略有上升的趨勢,但RMSECV的最小值要早于RMSEP出現(xiàn)。RMSECV最小值對應(yīng)的波長點(diǎn)數(shù)為54、重復(fù)選中次數(shù)為26,RMSEP最小值對應(yīng)的波長點(diǎn)數(shù)為135、重復(fù)選中次數(shù)為10,說明僅以RMSECV最小確定特征波長的方式容易導(dǎo)致回歸模型產(chǎn)生過擬合的問題。因此,選擇RMSEP最小時對應(yīng)的135個選中波長作為CARS-GSA優(yōu)選的乙酸特征波長。由圖3和圖5中RMSECV和RMSEP最小值的對比可知,CARS-GSA優(yōu)選特征波長的建模性能優(yōu)于CARS-500的建模性能。

    圖6 RMSE與波長數(shù)量間的關(guān)系

    2.2.3 特征波長優(yōu)選結(jié)果

    按上述方法執(zhí)行CARS500和CARS-GSA進(jìn)行丙酸和總酸特征波長優(yōu)選,得到101個丙酸特征波長和245個總酸特征波長。乙酸、丙酸和總酸特征波長分布情況如圖7所示。

    由圖7可知,VFA特征波長全部位于8 000 cm-1以下的中低頻區(qū)域,其中4 000~4 933、5 296~5 600和6 600~7 200 cm-1區(qū)域分布的特征波長點(diǎn)最多,這3部分正好對應(yīng)著光譜數(shù)據(jù)中吸收峰較強(qiáng)、分辨率較好的區(qū)域。丙酸特征波長在4 100~4 500 cm-1區(qū)域有53個,對應(yīng)著-CH2和-CH3基團(tuán)的組合頻;在4 000~4 900 cm-1區(qū)域有7個,對應(yīng)著C=O和-OH基團(tuán)的組合頻;在5 300~5 320 cm-1區(qū)域有2個,對應(yīng)著-COOH基團(tuán)的一級倍頻;在5 670~5 700 cm-1區(qū)域有9個,對應(yīng)著-CH2基團(tuán)的一級倍頻;在6 000~6 070 cm-1區(qū)域有13個,對應(yīng)著-CH3基團(tuán)的一級倍頻;在6 860~7 060 cm-1區(qū)域有17個,對應(yīng)著C=O、-CH2和-OH的二級倍頻??偹崽卣鞑ㄩL在4 000~4 720 cm-1區(qū)域有139個,對應(yīng)著C-C、C=O、-CH、-CH2和-CH3基團(tuán)的組合頻;在4 800~4 930 cm-1區(qū)域有27個,對應(yīng)著C=O和-OH基團(tuán)的組合頻;在5 300~5 380 cm-1區(qū)域有19個,對應(yīng)著-COOH基團(tuán)的一級倍頻;在5 930~6 010 cm-1區(qū)域有11個,對應(yīng)著-CH、-CH2和-CH3基團(tuán)的一級倍頻;在6 590~6 600 cm-1區(qū)域有2個,對應(yīng)著C=O基團(tuán)的二級倍頻;在6 730~7 200 cm-1區(qū)域有47個,對應(yīng)著C=O、-CH、-CH2、-CH3和-OH的二級倍頻。通過分析乙酸、丙酸和總酸特征波長可知,CARS-GSA與CARS500優(yōu)選特征波長結(jié)果具有很好的一致性,CARS-GSA只是剔除掉CARS500優(yōu)選特征波長中選中次數(shù)較少的相關(guān)性較差波長點(diǎn)。

    圖7 VFA特征波長優(yōu)選結(jié)果

    2.3 回歸模型評價與分析

    為評測2種波長優(yōu)選算法的性能,以CARS500和CARS-GSA優(yōu)選后的特征波長變量作為PLS回歸模型的輸入,建立沼液VFA定量回歸模型,并與全譜建模結(jié)果(Full-PLS)、單次CARS(運(yùn)行10次取最佳結(jié)果)優(yōu)選特征波長的建模效果進(jìn)行對比,結(jié)果如表2所示。

    表2 VFA PLS回歸模型評價指標(biāo)

    注:PCs是principal components的縮寫。

    Note: PCs is short for principal components.

    由表2可知,在單次CARS優(yōu)選特征波長建立的VFA回歸模型中,乙酸和丙酸CARS回歸模型的性能優(yōu)于全譜建模,而總酸CARS回歸模型的性能弱于全譜建模。原因在于乙酸和丙酸的結(jié)構(gòu)相對簡單,CARS能夠快速定位到相關(guān)性高的特征波長點(diǎn),而總酸的結(jié)構(gòu)相對復(fù)雜,不同基團(tuán)對應(yīng)的特征波長點(diǎn)數(shù)量較多,當(dāng)使用CARS剔除波長點(diǎn)時可能去掉某些相關(guān)性較高的特征波長點(diǎn),導(dǎo)致建模性能受到影響。多次執(zhí)行CARS算法進(jìn)行特征波長優(yōu)選可以解決單次CARS算法優(yōu)選總酸特征波長建模性能較差的問題。

    采用CARS-GSA作為厭氧發(fā)酵過程中發(fā)酵液乙酸、丙酸和總酸濃度的特征波長優(yōu)選方案,以優(yōu)選后的特征波長分別建立乙酸、丙酸和總酸濃度PLS回歸模型并進(jìn)行性能評測,其結(jié)果如圖8所示。

    圖8 VFA實測值與預(yù)測值分布

    3 結(jié) 論

    1)本研究采用NIRS結(jié)合化學(xué)計量學(xué)方法進(jìn)行沼液VFA的快速檢測,構(gòu)建模型進(jìn)行特征波長優(yōu)選,建立的乙酸、丙酸和總酸PLS回歸模型驗證決定系數(shù)分別為0.988、0.923和0.886,預(yù)測均方根誤差分別為0.111、0.120和0.727,且RPD都大于3,能夠滿足農(nóng)牧廢棄物厭氧發(fā)酵過程中對發(fā)酵液乙酸和丙酸濃度的快速檢測需求,基本滿足總酸濃度的檢測需求。

    2)CARS500在發(fā)酵液VFA濃度NIRS特征波長優(yōu)選方面具有良好的性能,通過多次執(zhí)行CARS算法并選取重復(fù)選中波長點(diǎn)作為特征波長的方式能夠有效提高建模精度和效率,并在一定程度上解決了CARS算法優(yōu)選特征波長結(jié)果的隨機(jī)性問題。

    3)CARS-GSA采用GSA對多次CARS優(yōu)選的特征波長進(jìn)行二次優(yōu)化,能夠有效去除CARS500優(yōu)選波長中相關(guān)性較弱的冗余波長點(diǎn),在提高建模精度和檢測效率的同時,確立了乙酸、丙酸和總酸相關(guān)基團(tuán)與其特征波長的對應(yīng)關(guān)系和特征波長在光譜區(qū)間內(nèi)的分布規(guī)律。

    [1]史風(fēng)梅,徐洪濤,盧玢宇,等. 溫度對養(yǎng)殖廢水厭氧發(fā)酵初期產(chǎn)酸的影響及其原因分析[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(增刊):42-47. Shi Fengmei, Xu Hongtao, Lu Binyu, et al. Effects of temperature on production of volatile fatty acids in mesophilic anaerobic fermentation of swine wastewater and its cause analysis[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(Supp.): 42-47. (in Chinese with English abstract)

    [2]于佳動,趙立欣,馮晶,等. 序批式秸稈牛糞混合厭氧干發(fā)酵過程物料理化及滲濾特性[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(20): 228-234. Yu Jiadong, Zhao Lixin, Feng Jing, et al. Physicochemical and percolating characteristics of sequencing batch dry anaerobic digestion of straw-cow manure mixture[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(20): 228-234. (in Chinese with English abstract)

    [3]宋香育,張克強(qiáng),房芳,等. 工藝措施對豬糞秸稈混合厭氧干發(fā)酵產(chǎn)氣性能的影響[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(11):233-239. Song Xiangyu, Zhang Keqiang, Fang Fang, et al. Influences of different technological strategies on performance of anaerobic co-digestion of pig manure with straw in solid-state[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(11): 233-239. (in Chinese with English abstract)

    [4]Sbarciog M, Giovannini G, Chamy R, et al. Control and estimation of anaerobic digestion processes using hydrogen and volatile fatty acids measurements[J]. Water Science and Technology, 2018, 78(10): 2027-2035.

    [5]Chatterjee B, Radhakrishnan L, Mazumder D. New approach for determination of volatile fatty acid in anaerobic digester sample[J]. Environmental Engineering Science, 2018, 35(4): 333-351.

    [6]Lahav O, Morgan B E. Titration methodologies for monitoring of anaerobic digestion in developing countries: A review[J]. Journal of Chemical Technology and Biotechnology, 2004, 79(12): 1331-1341.

    [7]Hey T, Sandstrom D, Ibrahim V, et al. Evaluating 5 and 8 pH-point titrations for measuring VFA in full-scale primary sludge hydrolysate[J]. Water SA, 2013, 39(1): 17-22.

    [8]Ward A J, Bruni E, Lykkegaard M K, et al. Real time monitoring of a biogas digester with gas chromatography, near-infrared spectroscopy, and membrane-inlet mass spectrometry[J]. Bioresource Technology, 2011, 102(5): 4098-4103.

    [9]Peris M, Escuder-Gilabert L. On-line monitoring of food fermentation processes using electronic noses and electronic tongues: A review[J]. Anal Chim Acta, 2013, 804: 29-36.

    [10]Yu Z, Leng X, Zhao S, et al. A review on the applications of microbial electrolysis cells in anaerobic digestion[J]. Bioresource Technology, 2018, 255: 340-348.

    [11]Nguyen Duc, Gadhamshetty Venkataramana, Nitayavardhana Saoharit, et al. Automatic process control in anaerobic digestion technology: A critical review[J]. Bioresource Technology, 2015, 193: 513-522.

    [12]Bruni E, Ward A J, Kocks M, et al. Comprehensive monitoring of a biogas process during pulse loads with ammonia[J]. Biomass & Bioenergy, 2013, 56: 211-220.

    [13]Falk H M, Reichling P, Andersen C, et al. Online monitoring of concentration and dynamics of volatile fatty acids in anaerobic digestion processes with mid-infrared spectroscopy[J]. Bioprocess and Biosystems Engineering, 2015, 38(2): 237-249.

    [14]Stockl Andrea, Lichti Fabian. Near-infrared spectroscopy (NIRS) for a real time monitoring of the biogas process[J]. Bioresource Technology, 2018, 247: 1249-1252.

    [15]Nespeca Maurílio Gustavo, Rodrigues Caroline Varella, Santana Kamili Oliveira, et al. Determination of alcohols and volatile organic acids in anaerobic bioreactors for H2production by near infrared spectroscopy[J]. International Journal of Hydrogen Energy, 2017, 42(32): 20480-20493.

    [16]Li L, Peng X Y, Wang X M, et al. Anaerobic digestion of food waste: A review focusing on process stability[J]. Bioresource Technology, 2018, 248: 20-28.

    [17]Rato Tiago J, Reis Marco S. Multiresolution interval partial least squares: A framework for waveband selection and resolution optimization[J]. Chemometrics and Intelligent Laboratory Systems, 2019, 186: 41-54.

    [18]鄒小波,張俊俊,黃曉瑋,等. 基于音頻和近紅外光譜融合技術(shù)的西瓜成熟度判別[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(9): 301-307. Zou Xiaobo, Zhang Junjun, Huang Xiaowei, et al. Distinguishing watermelon maturity based on acoustic characteristics and near infrared spectroscopy fusion technology[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(9): 301-307. (in Chinese with English abstract)

    [19]Zareef M, Chen Q S, Ouyang Q, et al. Prediction of amino acids, caffeine, theaflavins and water extract in black tea using FT-NIR spectroscopy coupled chemometrics algorithms[J]. Analytical Methods, 2018, 10(25): 3023-3031.

    [20]張亞坤,羅斌,宋鵬,等. 基于近紅外光譜的大豆葉片可溶性蛋白含量快速檢測[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(18):187-193. Zhang Yakun, Luo Bin, Song Peng, et al. Rapid determination of soluble protein content for soybean leaves based on near infrared spectroscopy[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(18): 187-193. (in Chinese with English abstract)

    [21]王巧華,梅璐,馬美湖,等. 利用機(jī)器視覺與近紅外光譜技術(shù)的皮蛋無損檢測與分級[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(24): 314-321. Wang Qiaohua, Mei Lu, Ma Meihu, et al. Nondestructive testing and grading of preserved duck eggs based on machine vision and near-infrared spectroscopy[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(24): 314-321. (in Chinese with English abstract)

    [22]Jiang H, Xu W, Chen Q. Comparison of algorithms for wavelength variables selection from near-infrared (NIR) spectra for quantitative monitoring of yeast (Saccharomyces cerevisiae) cultivations[J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2019, 214: 366-371.

    [23]Ye Dandan, Sun Laijun, Zou Borui, et al. Non-destructive prediction of protein content in wheat using NIRS[J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2018, 189: 463-472.

    [24]Zhu Y W, Chen X Y, Wang S M, et al. Simultaneous measurement of contents of liquirtin and glycyrrhizic acid in liquorice based on near infrared spectroscopy[J]. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2018, 196: 209-214.

    [25]楊瑋,李民贊,鄭立華,等. 冬棗氮素含量預(yù)測模型中特征波長選擇方法的應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(增刊2):164-168. Yang Wei, Li Minzan, Zheng Lihua, et al. Application of spectral screening method on prediction model of nitrogen content of jujube leaves[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(Supp.2): 164-168. (in Chinese with English abstract)

    [26]朱瑤迪,張佳燁,李苗云,等. 肽聚糖對肉制品中產(chǎn)氣莢膜梭菌芽孢萌發(fā)率影響及預(yù)測[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(4): 287-293. Zhu Yaodi, Zhang Jiaye, Li Miaoyun, et al. Effect of different peptidoglycan on clostridium perfringens spore germination and quantitative prediction[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(4): 287-293. (in Chinese with English abstract)

    [27]Yun Y H, Bin J, Liu D L, et al. A hybrid variable selection strategy based on continuous shrinkage of variable space in multivariate calibration[J]. Analytica Chimica Acta, 2019, 1058: 58-69.

    [28]Yang M, Xu D, Chen S, et al. Evaluation of machine learning approaches to predict soil organic matter and ph using vis-NIR spectra[J]. Sensors, 2019, 19(2): 263.

    [29]Song J, Li G, Yang X. Optimizing GA-PLS model of soluble solids content in Fukumoto navel orange based on Vis-NIR transmittance spectroscopy using discrete wavelet transform[J]. Journal of the Science of Food and Agriculture, 2019, 99(11): 4898-4903.

    [30]劉金明,初曉冬,王智,等. 玉米秸稈纖維素和半纖維素NIRS特征波長優(yōu)選[J]. 光譜學(xué)與光譜分析,2019,39(3):743-750. Liu Jinming, Chu Xiaodong, Wang Zhi, et al. Optimization of characteristic wavelength variables of near infrared spectroscopy for detecting contents of cellulose and hemicellulose in corn stover[J]. Spectroscopy and Spectral Analysis, 2019, 39(3): 743-750. (in Chinese with English abstract)

    [31]劉金明,程秋爽,甄峰,等. 基于GSA的厭氧發(fā)酵原料碳氮比NIRS快速檢測[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2019,50(11): 323-330. Liu Jinming, Cheng Qiushuang, Zhen Feng, et al. Rapid determination of C/N ratio for anaerobic fermentation feedstocks using near infrared spectroscopy based on GSA[J]. Transactions of the Chinese Society for Agricultural Machinery, 2019, 50(11): 323-330. (in Chinese with English abstract)

    [32]Liu Jinming, Li Nan, Zhen Feng, et al. Rapid detection of carbon-nitrogen ratio for anaerobic fermentation feedstocks using near-infrared spectroscopy combined with BiPLS and GSA[J]. Applied Optics, 2019, 58(18): 5090-5097.

    Rapid determination of volatile fatty acids in biogas slurry based on near infrared spectroscopy

    Liu Jinming1,2,3, Guo Kunlin2, Zhen Feng1,3, Zhang Hongqiong1,4, Li Wenzhe1,4, Xu Yonghua5※

    (1.,,150030,; 2.,,163319,; 3.,510640,; 4.,150030,; 5.,,150030,)

    Volatile Fatty Acids (VFA), serving as important intermediate products in Anaerobic Digestion (AD), have been considered as the key variables in most AD monitoring strategies, as they respond to incoming imbalances, indicating the buffer capacity of digesters to process disturbance and imminent digester failure that caused by sudden operational changes. In order to ensure efficient operation of AD while improve the utilization rate of raw materials, it is necessary to accurately monitor and evaluate the operation state of biogas engineering, via detecting the concentrations of VFA in the process of biogas production with corn stover and animal manure as feedstocks. Previously, the rapid detection models of Acetic Acid (AA), Propionic Acid (PA) and Total Acid (TA) in biogas slurry have been constructed, using the Near Infrared Spectroscopy (NIRS) technique combined with the Partial Least Squares (PLS), aiming to overcome the time consuming and high-cost in the traditional chemical analysis method. However, a prediction model can trigger the high complexity and low accuracy, due to the spectroscopic data generally includes quantities of invalid redundant information. In this study, an integrated algorithm was presented, based on the Competitive Adaptive Reweighted Sampling (CARS) and genetic simulated annealing algorithm (GSA), to optimize the characteristic wavelength variables of AA, PA, and TA, and thereby to improve the efficiency and precision of NIRS detection models. An AD experiment was carried out with corn stover, pig manure and cow manure as feedstocks, where 155 samples of biogas slurry were collected. The NIRS data of biogas slurry was acquired in a transmittance mode using the AntarisTMII FT-NIR spectrophotometer equipped with a quartz cuvette. A Gas Chromatography (GC) system was used to measure the VFA of biogas slurry, where 81 valid data of AA, 78 valid data of PA, and 87 valid data of TA were obtained to establish the regression model. One segment of the spectrum with 95 wavelength points was removed from 4 933.02 to 5 295.57 cm-1, and 1462 wavelength variables remained, mainly due to the saturation of spectrum can be caused by the strong combination band of -OH from water. The spectral preprocessing methods were selected, according to the mean relative error of calibration set. Correspondingly, the samples were divided into the calibration set and validation set, using Sample Set Portioning based on Joint X-Y Distances (SPXY) algorithm. The number of characteristic wavelength variables for AA, PA, and TA were 135, 101, and 245, respectively. The PLS regression models were established with the characteristic wavelengths of AA, PA, and TA, where the results were the coefficients of multiple determination for prediction is 0.988, root mean squared error of prediction (RMSEP) of 0.111, and the residual predictive deviation (RPD) of 9.685 for AA, coefficients of multiple determination for prediction is 0.922, RMSEP of 0.120, and RPD of 3.685 for PA, coefficients of multiple determination for prediction is 0.886, RMSEP of 0.727, and RPD of 3.484 for TA. Meanwhile, compared with the whole spectrum model, the RMSEP in the CARS-GSA model decreased by 17.78%, 15.49%, and 1.22%, respectively, showing that the number of wavelengths significantly decreased after the optimization, whereas, the performance of regressive model was obviously higher than that of the whole wavelengths. The results demonstrate that the CARS-GSA model can fulfil the requirement of rapid detection for AA and PA concentrations in biogas slurry during anaerobic fermentation with agricultural waste as feedstocks, while basically meet the detection requirement of TA concentration. The CARS-GSA model also can be used to enhance the forecasting capability of the model, while reduce its complexity. The findings can provide a new way to improve the accuracy and robustness of prediction model, base on optimizing sensitive wavelengths for AA, PA, and TA, further for rapid and accurate measurement of VFA concentrations in biogas slurry.

    anaerobic digestion; volatile fatty acids; rapid determination; near infrared spectroscopy; partial least squares; genetic simulated annealing algorithm; competitive adaptive reweighted sampling

    劉金明,郭坤林,甄峰,等. 基于近紅外光譜的沼液揮發(fā)性脂肪酸含量快速檢測[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(18):188-196.doi:10.11975/j.issn.1002-6819.2020.18.023 http://www.tcsae.org

    Liu Jinming, Guo Kunlin, Zhen Feng, et al. Rapid determination of volatile fatty acids in biogas slurry based on near infrared spectroscopy[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(18): 188-196. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.18.023 http://www.tcsae.org

    2020-05-10

    2020-06-28

    中國科學(xué)院可再生能源重點(diǎn)實驗室(Y907k81001);國家重點(diǎn)研發(fā)計劃(2019YFD1100603);黑龍江省博士后面上資助(LBH-Z19087);黑龍江八一農(nóng)墾大學(xué)三橫三縱支持計劃(ZRCQC202007);黑龍江八一農(nóng)墾大學(xué)學(xué)成人才科研啟動計劃(XDB202006)

    劉金明,博士,副教授,主要從事光譜分析技術(shù)在農(nóng)業(yè)領(lǐng)域的應(yīng)用研究。Email:jinmingliu2008@126.com

    許永花,副教授,主要從事光譜分析技術(shù)方面的研究。Email:xyhsy@126.com

    10.11975/j.issn.1002-6819.2020.18.023

    O657.33

    A

    1002-6819(2020)-18-0188-09

    猜你喜歡
    厭氧發(fā)酵總酸丙酸
    響應(yīng)面法優(yōu)化辣椒紅酸湯的配方及接種發(fā)酵工藝
    葡萄皮渣醋制作工藝優(yōu)化及抗氧化活性研究
    餐廚垃圾厭氧發(fā)酵熱電氣聯(lián)供系統(tǒng)優(yōu)化
    初探影響液態(tài)法米酒總酸總酯的因素
    法國梧桐落葉、香樟青葉與豬糞混合厭氧發(fā)酵特性的探究
    滴定法測定食醋中總酸含量及其不確定度評定
    太陽能-地能熱泵耦合系統(tǒng)在沼氣工程厭氧發(fā)酵增溫中的設(shè)計與應(yīng)用
    食品中丙酸鈉、丙酸鈣測定方法的改進(jìn)
    餐廚垃圾厭氧發(fā)酵產(chǎn)沼氣的初步探究
    K/γ-Al2O3催化丙酸甲酯合成甲基丙烯酸甲酯
    精品久久久久久久久久免费视频| 美女高潮的动态| 夜夜爽天天搞| 乱人视频在线观看| 一夜夜www| 免费大片18禁| 超碰av人人做人人爽久久 | 狠狠狠狠99中文字幕| 亚洲国产精品sss在线观看| 女同久久另类99精品国产91| 波多野结衣巨乳人妻| svipshipincom国产片| 很黄的视频免费| 精品久久久久久久末码| 观看美女的网站| 美女黄网站色视频| 2021天堂中文幕一二区在线观| 三级男女做爰猛烈吃奶摸视频| 丁香六月欧美| 亚洲不卡免费看| 很黄的视频免费| 最后的刺客免费高清国语| 色播亚洲综合网| 在线观看午夜福利视频| 欧美性猛交╳xxx乱大交人| 男女做爰动态图高潮gif福利片| 在线视频色国产色| 波多野结衣高清无吗| 日韩欧美 国产精品| 成年女人看的毛片在线观看| 精品一区二区三区视频在线观看免费| 99视频精品全部免费 在线| 亚洲最大成人中文| 特级一级黄色大片| 中文资源天堂在线| 小蜜桃在线观看免费完整版高清| 亚洲国产中文字幕在线视频| 熟女人妻精品中文字幕| 亚洲成av人片免费观看| 黄片小视频在线播放| 日本精品一区二区三区蜜桃| 久久久久久九九精品二区国产| 中文字幕人成人乱码亚洲影| www.999成人在线观看| 日日干狠狠操夜夜爽| 男女之事视频高清在线观看| 亚洲精品国产精品久久久不卡| 精品国产三级普通话版| 搡女人真爽免费视频火全软件 | 亚洲国产精品久久男人天堂| 级片在线观看| 欧美日韩国产亚洲二区| 在线观看66精品国产| 久久久久久人人人人人| 久久精品人妻少妇| 成人特级黄色片久久久久久久| 国产精品国产高清国产av| 成人永久免费在线观看视频| 日日夜夜操网爽| 在线播放国产精品三级| 亚洲一区高清亚洲精品| 女警被强在线播放| 日本一二三区视频观看| 国产高清有码在线观看视频| 90打野战视频偷拍视频| 久久精品国产亚洲av香蕉五月| 男人舔奶头视频| 亚洲成人精品中文字幕电影| 亚洲激情在线av| 午夜亚洲福利在线播放| 国产一区二区亚洲精品在线观看| 亚洲av美国av| 国产欧美日韩一区二区三| eeuss影院久久| 亚洲成人精品中文字幕电影| 亚洲精品亚洲一区二区| 亚洲中文字幕一区二区三区有码在线看| 三级男女做爰猛烈吃奶摸视频| 日韩国内少妇激情av| 亚洲无线观看免费| 成人精品一区二区免费| 一级毛片女人18水好多| 可以在线观看毛片的网站| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产一区在线观看成人免费| 亚洲 国产 在线| 在线天堂最新版资源| 麻豆成人av在线观看| 国产精品国产高清国产av| 香蕉久久夜色| 在线播放无遮挡| 99热精品在线国产| 级片在线观看| 美女被艹到高潮喷水动态| 人妻久久中文字幕网| 国产中年淑女户外野战色| 毛片女人毛片| 99国产精品一区二区蜜桃av| 国产一区二区三区在线臀色熟女| 免费看美女性在线毛片视频| 丰满的人妻完整版| 美女免费视频网站| avwww免费| 欧美黄色片欧美黄色片| 好男人电影高清在线观看| 久久香蕉精品热| 欧美日韩瑟瑟在线播放| 国产亚洲精品一区二区www| 精品久久久久久久人妻蜜臀av| 韩国av一区二区三区四区| 看片在线看免费视频| 亚洲欧美日韩高清在线视频| 国产精品永久免费网站| 欧美一级毛片孕妇| 欧美黄色淫秽网站| 一进一出抽搐动态| 国产真实乱freesex| 久久久久国内视频| 亚洲国产精品成人综合色| 51午夜福利影视在线观看| 露出奶头的视频| 国产成人福利小说| 国产精品免费一区二区三区在线| 午夜影院日韩av| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av电影不卡..在线观看| 一级作爱视频免费观看| 日韩有码中文字幕| 亚洲成人久久爱视频| 男插女下体视频免费在线播放| 亚洲国产精品成人综合色| 成人三级黄色视频| 午夜免费成人在线视频| 亚洲精品在线观看二区| 亚洲美女黄片视频| 精品久久久久久,| 国产男靠女视频免费网站| 日韩精品青青久久久久久| 老熟妇仑乱视频hdxx| 两个人的视频大全免费| 成人av一区二区三区在线看| 国内精品美女久久久久久| av在线蜜桃| 色哟哟哟哟哟哟| 国产一区二区在线观看日韩 | 男女做爰动态图高潮gif福利片| avwww免费| 久久婷婷人人爽人人干人人爱| 国产伦精品一区二区三区视频9 | 国产精品精品国产色婷婷| 欧美一区二区亚洲| а√天堂www在线а√下载| 最近最新中文字幕大全免费视频| 久久久久久久精品吃奶| 日本 欧美在线| 久久久久九九精品影院| 欧美乱妇无乱码| 午夜精品久久久久久毛片777| bbb黄色大片| 97超级碰碰碰精品色视频在线观看| 成人鲁丝片一二三区免费| 久久草成人影院| 老司机午夜福利在线观看视频| 国产精品亚洲av一区麻豆| 国产 一区 欧美 日韩| 国产av麻豆久久久久久久| 757午夜福利合集在线观看| 亚洲一区二区三区色噜噜| 97人妻精品一区二区三区麻豆| 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产清高在天天线| 亚洲av电影在线进入| 母亲3免费完整高清在线观看| 国产精品嫩草影院av在线观看 | 中文字幕av在线有码专区| 国产探花极品一区二区| 老司机在亚洲福利影院| 不卡一级毛片| 国产不卡一卡二| 亚洲精品色激情综合| 欧美成人a在线观看| 久99久视频精品免费| 中出人妻视频一区二区| 怎么达到女性高潮| 中文字幕高清在线视频| 香蕉久久夜色| 麻豆成人av在线观看| 久久久久久人人人人人| 3wmmmm亚洲av在线观看| www.999成人在线观看| 综合色av麻豆| 一二三四社区在线视频社区8| 日日干狠狠操夜夜爽| 可以在线观看的亚洲视频| 亚洲成人久久爱视频| 91麻豆av在线| 国产激情偷乱视频一区二区| 国产毛片a区久久久久| 人妻夜夜爽99麻豆av| 国产真人三级小视频在线观看| 九色成人免费人妻av| 女人十人毛片免费观看3o分钟| 国产亚洲精品久久久com| 国产精品美女特级片免费视频播放器| 国内精品久久久久久久电影| 国产精品久久电影中文字幕| 亚洲人成网站高清观看| 国产伦一二天堂av在线观看| 欧美日韩精品网址| 午夜福利高清视频| 草草在线视频免费看| 国产一区二区亚洲精品在线观看| 俄罗斯特黄特色一大片| 在线观看午夜福利视频| 一进一出抽搐动态| 国产精品 欧美亚洲| 精品久久久久久久久久免费视频| 韩国av一区二区三区四区| 精品不卡国产一区二区三区| 午夜免费成人在线视频| 国产免费av片在线观看野外av| 丰满乱子伦码专区| 免费观看的影片在线观看| 校园春色视频在线观看| 亚洲av成人精品一区久久| 国产免费男女视频| 国内精品美女久久久久久| 热99在线观看视频| 国产伦人伦偷精品视频| www.www免费av| 国产成人欧美在线观看| 午夜视频国产福利| 日韩国内少妇激情av| 亚洲专区中文字幕在线| 黄片小视频在线播放| 十八禁人妻一区二区| 国内毛片毛片毛片毛片毛片| 精品久久久久久久末码| 欧美黄色淫秽网站| 黑人欧美特级aaaaaa片| 又黄又粗又硬又大视频| 91久久精品电影网| 国产私拍福利视频在线观看| 亚洲av电影在线进入| 国产亚洲精品久久久com| 久久伊人香网站| 搞女人的毛片| av天堂在线播放| 19禁男女啪啪无遮挡网站| 最近最新免费中文字幕在线| 夜夜夜夜夜久久久久| 在线播放无遮挡| 神马国产精品三级电影在线观看| 91麻豆精品激情在线观看国产| 18禁在线播放成人免费| www国产在线视频色| 欧美日韩精品网址| 宅男免费午夜| 国产麻豆成人av免费视频| 国产色婷婷99| 熟女电影av网| 岛国视频午夜一区免费看| 成人无遮挡网站| 精品午夜福利视频在线观看一区| 窝窝影院91人妻| 日韩中文字幕欧美一区二区| 国产精品一区二区三区四区久久| 亚洲五月天丁香| 国产精品野战在线观看| 悠悠久久av| 国产伦人伦偷精品视频| 国产又黄又爽又无遮挡在线| 亚洲成人久久性| 深爱激情五月婷婷| 精品无人区乱码1区二区| 国产毛片a区久久久久| 婷婷六月久久综合丁香| 69av精品久久久久久| e午夜精品久久久久久久| 人妻丰满熟妇av一区二区三区| 国产精品国产高清国产av| 久久国产乱子伦精品免费另类| 欧美一区二区亚洲| 日韩有码中文字幕| 国产 一区 欧美 日韩| 成年版毛片免费区| 少妇人妻一区二区三区视频| 18禁美女被吸乳视频| 午夜激情欧美在线| 欧美一级a爱片免费观看看| 99久久九九国产精品国产免费| 亚洲美女黄片视频| 一个人免费在线观看的高清视频| 成人av在线播放网站| 国产在线精品亚洲第一网站| 网址你懂的国产日韩在线| 亚洲国产精品sss在线观看| 听说在线观看完整版免费高清| 在线看三级毛片| 亚洲成av人片在线播放无| 此物有八面人人有两片| 成年女人永久免费观看视频| 国产真实乱freesex| 精品一区二区三区av网在线观看| 99久久九九国产精品国产免费| 少妇人妻一区二区三区视频| 极品教师在线免费播放| 最新在线观看一区二区三区| 国产精品99久久久久久久久| 亚洲国产精品999在线| xxxwww97欧美| 日本精品一区二区三区蜜桃| 在线观看美女被高潮喷水网站 | 亚洲色图av天堂| 国产精品影院久久| 19禁男女啪啪无遮挡网站| 欧美成狂野欧美在线观看| 蜜桃久久精品国产亚洲av| 老司机深夜福利视频在线观看| 中文字幕av成人在线电影| 99精品欧美一区二区三区四区| 18美女黄网站色大片免费观看| 国产成人欧美在线观看| 日韩欧美在线乱码| 悠悠久久av| 久久中文看片网| 中亚洲国语对白在线视频| 色综合欧美亚洲国产小说| 香蕉久久夜色| 国产三级中文精品| 变态另类丝袜制服| 成年人黄色毛片网站| 国产精品永久免费网站| 久久久久国产精品人妻aⅴ院| 日本与韩国留学比较| 亚洲国产精品999在线| 九色成人免费人妻av| 给我免费播放毛片高清在线观看| 天堂影院成人在线观看| 老熟妇乱子伦视频在线观看| 在线视频色国产色| 天美传媒精品一区二区| 日本撒尿小便嘘嘘汇集6| 国产精品电影一区二区三区| 免费在线观看影片大全网站| 久久久久久久久久黄片| www.www免费av| 尤物成人国产欧美一区二区三区| 亚洲色图av天堂| 少妇丰满av| 啪啪无遮挡十八禁网站| 亚洲男人的天堂狠狠| av天堂中文字幕网| 久久婷婷人人爽人人干人人爱| 一个人观看的视频www高清免费观看| 欧美成人免费av一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美精品v在线| 香蕉丝袜av| 免费在线观看成人毛片| 色综合婷婷激情| 91字幕亚洲| 乱人视频在线观看| xxxwww97欧美| 久久久久免费精品人妻一区二区| 亚洲精品一区av在线观看| 99国产综合亚洲精品| 日本撒尿小便嘘嘘汇集6| 欧美乱码精品一区二区三区| 成人三级黄色视频| 两个人视频免费观看高清| 亚洲av电影不卡..在线观看| 亚洲色图av天堂| 精品乱码久久久久久99久播| 桃红色精品国产亚洲av| 国产在线精品亚洲第一网站| 欧美日韩黄片免| 亚洲熟妇中文字幕五十中出| 精品欧美国产一区二区三| 啦啦啦观看免费观看视频高清| 午夜福利高清视频| 国产一区二区三区在线臀色熟女| 最新在线观看一区二区三区| 欧美+日韩+精品| 国产一区在线观看成人免费| 嫩草影院入口| 免费av观看视频| 在线免费观看的www视频| 欧美中文日本在线观看视频| 九九在线视频观看精品| 在线免费观看的www视频| 国产中年淑女户外野战色| 99国产精品一区二区三区| 色吧在线观看| 亚洲一区高清亚洲精品| 黄色日韩在线| 一二三四社区在线视频社区8| 黄色片一级片一级黄色片| 好看av亚洲va欧美ⅴa在| 热99在线观看视频| 亚洲国产色片| 在线播放无遮挡| 精品国产三级普通话版| 桃红色精品国产亚洲av| 一区二区三区国产精品乱码| 久久亚洲精品不卡| 99久久精品国产亚洲精品| 成人特级黄色片久久久久久久| 天堂√8在线中文| 国内精品久久久久久久电影| 欧美绝顶高潮抽搐喷水| 美女免费视频网站| 啪啪无遮挡十八禁网站| 观看免费一级毛片| 欧美一区二区精品小视频在线| 乱人视频在线观看| 日韩欧美在线二视频| 久久久久久人人人人人| 欧美激情在线99| 国产精品一区二区免费欧美| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲精品一区二区www| 女人高潮潮喷娇喘18禁视频| 午夜福利免费观看在线| 特大巨黑吊av在线直播| 亚洲欧美日韩东京热| 色综合亚洲欧美另类图片| 中文字幕精品亚洲无线码一区| 在线观看免费视频日本深夜| 久久精品综合一区二区三区| 特大巨黑吊av在线直播| 精品一区二区三区人妻视频| 亚洲av一区综合| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区三区视频了| 天天躁日日操中文字幕| 香蕉久久夜色| 男女床上黄色一级片免费看| 很黄的视频免费| 白带黄色成豆腐渣| 2021天堂中文幕一二区在线观| 在线观看美女被高潮喷水网站 | 真人做人爱边吃奶动态| 国产一级毛片七仙女欲春2| 国产精品自产拍在线观看55亚洲| av黄色大香蕉| 午夜福利18| 久久性视频一级片| 欧美绝顶高潮抽搐喷水| a在线观看视频网站| 亚洲国产色片| 国产精品国产高清国产av| 亚洲性夜色夜夜综合| 美女被艹到高潮喷水动态| 91麻豆精品激情在线观看国产| 99久久综合精品五月天人人| 亚洲成人精品中文字幕电影| 国产麻豆成人av免费视频| 麻豆久久精品国产亚洲av| 亚洲av免费高清在线观看| 亚洲电影在线观看av| 午夜福利视频1000在线观看| 欧美区成人在线视频| 十八禁人妻一区二区| 九九在线视频观看精品| 国模一区二区三区四区视频| 小说图片视频综合网站| 亚洲五月天丁香| 99久久九九国产精品国产免费| 国产黄色小视频在线观看| 国产伦在线观看视频一区| 国产精品精品国产色婷婷| 亚洲在线观看片| 深夜精品福利| av黄色大香蕉| 色噜噜av男人的天堂激情| x7x7x7水蜜桃| 久久欧美精品欧美久久欧美| 最近在线观看免费完整版| 少妇高潮的动态图| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕人妻熟人妻熟丝袜美 | 国产伦人伦偷精品视频| 精品久久久久久久久久免费视频| eeuss影院久久| 亚洲一区高清亚洲精品| 国产午夜精品久久久久久一区二区三区 | 特级一级黄色大片| 男女床上黄色一级片免费看| 免费看日本二区| 两个人的视频大全免费| 女人十人毛片免费观看3o分钟| 可以在线观看毛片的网站| 午夜a级毛片| 在线免费观看不下载黄p国产 | av视频在线观看入口| 精品国内亚洲2022精品成人| 国产一区在线观看成人免费| 女警被强在线播放| 在线观看免费视频日本深夜| 18禁在线播放成人免费| 国产探花在线观看一区二区| 男女床上黄色一级片免费看| 国产又黄又爽又无遮挡在线| 国产主播在线观看一区二区| 国产视频内射| 身体一侧抽搐| 中亚洲国语对白在线视频| 9191精品国产免费久久| 国产精品一及| 午夜福利成人在线免费观看| 大型黄色视频在线免费观看| 久久天躁狠狠躁夜夜2o2o| 91在线观看av| 长腿黑丝高跟| www.熟女人妻精品国产| 亚洲真实伦在线观看| 亚洲av免费高清在线观看| 波野结衣二区三区在线 | 精品熟女少妇八av免费久了| 变态另类成人亚洲欧美熟女| 欧美一区二区亚洲| 美女被艹到高潮喷水动态| 欧美午夜高清在线| 99久久综合精品五月天人人| 午夜福利在线在线| 麻豆国产av国片精品| 一级黄色大片毛片| 狠狠狠狠99中文字幕| 亚洲国产欧美网| 波多野结衣高清作品| svipshipincom国产片| 日韩欧美精品v在线| 18禁黄网站禁片午夜丰满| 国产午夜福利久久久久久| 亚洲五月天丁香| 日本熟妇午夜| 听说在线观看完整版免费高清| 免费在线观看日本一区| 日韩欧美一区二区三区在线观看| www国产在线视频色| 男人舔奶头视频| 麻豆成人午夜福利视频| 女同久久另类99精品国产91| 欧美zozozo另类| 国产成+人综合+亚洲专区| 亚洲av日韩精品久久久久久密| 国产色婷婷99| 在线免费观看的www视频| 成人18禁在线播放| 色噜噜av男人的天堂激情| 欧美3d第一页| av中文乱码字幕在线| 国产美女午夜福利| 国产高潮美女av| 一本一本综合久久| 在线天堂最新版资源| 亚洲欧美日韩高清专用| 亚洲专区中文字幕在线| 男人和女人高潮做爰伦理| 亚洲精品在线美女| 中国美女看黄片| 国语自产精品视频在线第100页| 狠狠狠狠99中文字幕| 久久精品亚洲精品国产色婷小说| 亚洲av日韩精品久久久久久密| 国产色婷婷99| 五月玫瑰六月丁香| 国产高清三级在线| 国产精品一及| 久久久久精品国产欧美久久久| 一区福利在线观看| 一卡2卡三卡四卡精品乱码亚洲| 在线视频色国产色| 欧美中文日本在线观看视频| 人人妻人人澡欧美一区二区| 日本撒尿小便嘘嘘汇集6| 国产精品久久久久久亚洲av鲁大| 别揉我奶头~嗯~啊~动态视频| 身体一侧抽搐| 一本一本综合久久| 亚洲精品在线美女| 亚洲成av人片免费观看| 国产99白浆流出| 久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 日韩高清综合在线| 九九在线视频观看精品| 18禁裸乳无遮挡免费网站照片| 小说图片视频综合网站| 免费看十八禁软件| 久久国产精品影院| 日韩欧美精品免费久久 | 黄色成人免费大全| 久久精品91无色码中文字幕| 亚洲专区国产一区二区| 欧美xxxx黑人xx丫x性爽| av福利片在线观看| 色精品久久人妻99蜜桃| 国产亚洲精品久久久久久毛片| 十八禁人妻一区二区| 精品一区二区三区视频在线观看免费| 青草久久国产| 人妻夜夜爽99麻豆av| 欧美日韩福利视频一区二区| 久久欧美精品欧美久久欧美| 亚洲成a人片在线一区二区| 无遮挡黄片免费观看| 国产v大片淫在线免费观看| 色噜噜av男人的天堂激情| 日本在线视频免费播放| 丰满的人妻完整版| 最近最新中文字幕大全电影3| 99精品欧美一区二区三区四区| 啦啦啦韩国在线观看视频|