天天操天天玩_19国产精品麻豆免费观看_欧美一级一级做性视频_国产欧美在线播放

歡迎訪問科曼環保網~

免費發布 登錄 注冊 退出

環保行業信息
資訊服務平臺

MBR中膜污染的人工神經網絡預測研究進展

摘要:

膜生物反應器(MBR)作為一種新型廢水處理技術在污水處理方面具有廣闊的應用前景。但是,膜污染是制約MBR進一步發展的瓶頸性問題。近年來,隨著數學算法及計算機技術的發展,將人工神經網絡(ANN)等機器學習算法應用于MBR的膜污染預測成為研究的熱點。總結了膜污染的影響因素,探討了基于經典數學模型膜污染預測的優缺點,綜述了近年來國內外學者運用簡單ANN、優化算法ANN和深度學習ANN對MBR膜污染預測的研究,提出優化算法ANN與深度學習ANN在面對復雜環境下更具優勢。此外,還探討了當前運用ANN機器學習算法進行膜污染預測存在的缺陷,指出ANN模型在中試和工業化規模的MBR膜污染預測中應用較少,并對其未來的發展進行了展望。

膜生物反應器(Membrane bioreactor,MBR)是一種將活性污泥法與膜分離技術相結合的廢水處理技術,具有出水水質好、負荷高和占地面積小等優點〔1-2〕。目前,大型(10 000 m3/d)和超大型(100 000 m3/d)MBR污水處理廠在包括中國、美國和歐洲等世界各地相繼建成并投入使用。據預測,中國的大型MBR污水處理廠將超過300座,總污水處理規模可達1 500萬t/d〔2〕。但是,膜污染會導致跨膜壓差(Transmembrane pressure,TMP)的上升和膜通量的下降,增加MBR的運行成本,成為限制MBR廣泛應用的瓶頸性問題〔3〕。

機器學習(Machine learning,ML)是一種基于以往經驗或數據進行分類與回歸的技術,不僅被運用于農業、氣候、安全、教育、醫學等行業〔4〕,也被廣泛運用于污水處理領域。人工神經網絡(Artificial neural network,ANN)作為一種機器學習模型,不僅具有參數少、預測性能好、泛化能力強等優點〔5〕,而且具有極強的非線性映射能力、學習能力和黑箱建模能力。因此,應用ANN模型進行膜污染預測無需考慮復雜的膜污染中間過程,只需考慮與膜污染相關的指標即可。石寶強等〔6〕早在2006年就綜述了ANN在MBR膜污染預測中的應用研究,但是隨著對MBR膜污染機理的深入研究和ANN模型理論與技術的發展,亟需新的綜述對近年來的研究進行總結。而F. SCHMITT等〔7-8〕雖然綜述了人工智能用于預測膜污染的研究,但是欠缺對傳統數學模型與神經網絡模型差異的深入分析。基于此,筆者簡述了影響膜污染的因素,對比了傳統數學模型與人工神經網絡模型的差異與優劣,總結了近年來ANN模型在膜污染預測中的應用及優缺點,并對其未來的發展進行了展望。

1 膜污染的影響因素

膜材料、反應器操作條件和污泥混合液特性是與膜污染密切相關的主要因素〔9〕,其涉及到的具體參數見圖1。

1.jpg

由圖1可知,膜污染的影響因素眾多且相互影響,因此很難厘清各自的具體作用。例如,進水水質、污泥停留時間(Sludge retention time,SRT)和水力停留時間(Hydraulic retention time,HRT)等操作條件的改變不僅會直接影響膜污染,而且會導致污泥混合液特性發生變化,進而改變膜污染速率〔10〕。

早期研究人員建立了關于操作條件、膜通量、過濾濾阻等參數的數學模型〔11-15〕。此后,有研究建立了表達污泥混合液特性〔例如混合液懸浮固體濃度(Mixed liquid suspend solids,MLSS)、顆粒粒徑(Particle size distribution,PSD)、溶解性微生物產物(Soluble microbial products,SMP)等〕與膜污染關系的數學模型〔16-18〕。但是,這些傳統的數學模型一般為指數式經驗模型,往往需要通過諸多假設來簡化計算,且各參數物理意義不明確,存在預測精度低和通用性差等問題〔6-7〕。因此,開發能夠更加精準預測膜污染的模型對膜污染控制及其機理研究有重要意義。

膜污染模型的構建需要選擇合適的影響因子作為輸入變量,并尋找變量與膜污染之間的關系。因此,厘清各膜污染因素之間的相互作用及其對膜污染的影響是構建膜污染預測模型的前提條件。

1.1 膜材料性質

膜材料的性質包括膜的材質、親疏水性和粗糙程度等〔19〕,這些性質均會對污染物在膜表面的遷移轉化產生顯著影響,進而影響膜污染。例如,Y. JEONG等〔20〕比較了Al2O3陶瓷膜與聚偏氟乙烯(Polyvinylidene fluoride,PVDF)膜污染狀況,結果發現陶瓷膜具有更低的膜污染速率。膜的親疏水性也與膜污染密切相關,大量研究通過表面改性和共混改性法提高了膜表面的親水性能,從而減緩了膜污染〔21〕。此外,膜表面粗糙程度也與膜污染有關,例如有學者發現粗糙的膜表面有利于膠體顆粒在膜上的積聚,引發膜孔堵塞,從而導致嚴重的膜污染〔19〕。

1.2操作參數

不同的MBR系統操作條件也會影響膜污染,主要包括進水水質、運行溫度、SRT、HRT、錯流流速和曝氣強度等。SRT與HRT是影響MBR處理效能、污泥混合液特性和膜污染的重要因素。一般認為SRT的增加會降低污泥活性和絮體大小,此外較長的SRT促進了內源性衰變和細胞裂解,導致SMP的增加〔22〕。例如,在厭氧MBR的研究中,Zhi HUANG等〔23〕比較了不同SRT(30 d、60 d、∞)對膜污染的影響,結果表明在較長的SRT下厭氧MBR中會積累較多的蛋白質與多糖,MLSS增加,從而增加了膜污染速率。對于HRT而言,較短的HRT會導致有機負荷的增加,從而導致微生物分泌更多的胞外聚合物(Extracellular polymeric substances,EPS)和SMP,從而提升膜污染速率。例如,N. FALLAH等〔24〕比較了不同HRT對膜污染的影響,結果表明隨著HRT的降低,EPS和SMP顯著增加,污泥混合液的黏度也顯著增加,導致了較高的膜污染速率。此外,通過在膜表面增加錯流流速和曝氣強度等物理剪切方式,減少顆粒物在膜表面的沉積也是減緩膜污染的常見方法。Yingyu AN等 〔25〕在厭氧MBR中發現,一定范圍內增加含甲烷生物氣曝氣強度能夠減緩膜污染,但是過高的曝氣強度會破壞污泥絮體,導致污泥顆粒平均粒徑降低,從而更易形成致密的膜污染層,同時膜污染速率也會提高。

1.3 污泥混合液特性

對膜污染存在影響的污泥混合液特性主要包括MLSS、EPS和SMP等。盡管MLSS并不是影響膜污染的主要因素,但是過高的MLSS(>15 g/L)會增加混合液的黏度,導致膜污染速率的提高,過低的MLSS(<6 g/L)則會降低SRT并增加食物微生物比率,同樣導致較高的膜污染速率〔26〕。EPS和SMP主要是由微生物生長代謝和衰敗過程中所釋放的蛋白質、多糖、核酸、脂類物質以及腐殖酸等有機物組成,其也被證明與膜污染密切相關,有研究表明EPS中多糖和蛋白質的含量與膜污染速率呈正相關關系〔27〕。此外,SMP不僅可以吸附在膜表面,而且可以堵塞膜孔,從而降低膜面積并提高膜污染速率。

2 膜污染預測模型的相關研究

2.1 經典數學模型

在過去幾十年里,研究人員相繼構建了表達污泥混合液特性和運行條件等與膜污染關系的經驗公式,并建立了經典數學模型,其部分代表性成果見表1。

1.jpg

經典串聯阻力模型是一種使用達西定律來描述不同膜污染層的過濾濾阻模型,數十年來被廣泛用于膜通量的預測〔11〕。基于該模型,B. TANSEL等〔12〕構建了通量模型,在死端過濾與錯流過濾初期成功地預測了通量下降到60%所需的時間。S. CHANG等〔13〕結合達西定律與Hagen-Poisenille公式〔14〕構建了膜纖維半徑與膜通量關系的模型,描述了中空纖維膜中的膜污染情況。J. HERMIA等〔15〕基于膜通量下降的孔隙阻塞機制,根據膜孔堵塞情況選取不同n值(n=0,1,1.5,2),建立了通量預測模型以解釋膜通量與膜污染間的關系。

研究人員還構建了污泥混合液特性與膜污染間的數學模型。例如M. F. R. ZUTHI等〔16〕以濾餅層模型與阻塞模型為基礎,構建了SMP與膜污染間的半經驗數學模型,對TMP的變化進行了精準的預測。T. SATO等〔17〕建立了MLSS、COD與過濾濾阻之間的關系公式,其預測誤差為±24%。此外,Fangang MENG等〔18〕基于MLSS、PSD和EPS,構建了膜污染預測模型,用以闡明污泥混合液特性與膜污染的關系。

盡管經典數學模型在膜污染預測方面取得了一定的成果,但是膜污染是一個受進水特性、微生物活動、過濾模式等多因素影響的復雜動態過程,這些經典數學模型往往為了簡化僅考慮部分影響因素或在建立模型過程中設置諸多假設,例如,基于串聯阻力定律的預測模型雖然有廣泛的應用,但并沒有考慮生物膜的生長和不可逆污染物的形成等作用對膜污染的影響〔28〕,因而在實際應用中不可避免會存在一定的誤差。

2.2 機器學習常見算法

機器學習能從大量歷史數據中挖掘其隱藏的規律,并做出分類或回歸預測。機器學習常見算法及優缺點見表2〔29-33〕。

1.jpg

由表2可知,機器學習能夠彌補傳統的數學模型預測精度低、實際應用效果差等缺點〔29-33〕。目前機器學習已被運用到MBR膜污染的研究中,不少學者使用隨機森林、SVM和多元回歸等來建立與膜污染相關的預測模型。隨機森林是一種以決策樹為基學習器的集成算法,李威威等〔31〕采用隨機森林算法預測膜通量,通過模型計算得到的預測值與真實值的平均誤差僅3.98%。SVM具有的核函數可以解決非線性問題,梁楷〔32〕基于SVM建立的膜污染預測模型有效地預測了膜通量,平均誤差僅為2.63%。多元回歸能選取自變量的最優組合共同預測因變量,Zhan WANG等〔33〕通過多元回歸建立了化學清洗模型,通過模型預測通量恢復率,對減緩膜污染起到了指導性作用。盡管以上機器學習算法能夠在膜污染預測方面獲得較高的預測精度,但是仍然存在一定的局限性,例如,雖然隨機森林對數據集的適應能力很強并具有較快的運算速度,但其往往忽略數據集之間的相關性,導致預測精度下降。

鑒于膜污染過程的復雜性和非線性,簡單的機器學習只能粗略表示變量之間的關系,無法取得較高的預測精度。近年來,由于ANN能夠通過學習非線性的復雜關系來構建數學模型,并且具有較強的泛化能力,對未知數據具有較好的預測效果,因此,部分學者逐漸將其應用于膜污染的預測。對Web of Science上2006年至2020年運用ANN算法研究膜污染的文獻數量進行統計,結果見圖2。

1.jpg

由圖2可知,從2016年開始,關于ANN在膜污染方面研究的文獻數量迅速增加,可見ANN已經成為研究膜污染預測及控制的重要工具。

2.3 人工神經網絡算法

ANN是一種模仿動物神經網絡行為特征進行分布式并行信息處理的算法。其簡單拓撲圖見圖3。

1.jpg

如圖3所示,ANN信息處理功能的實現主要依靠神經節點的輸入輸出、神經元的閾值及連接權值大小的調整。ANN的預測性能可以通過均方根誤差(Root mean square error,RMSE)、決定系數R2、相對誤差Error等參數來評估,其計算分別見公式(1)~(3)。

1.jpg

公式(1)~(3)中,yt(i)、y(i)yt(i)、y(i)分別表示第i個預測和輸出值,y?y?表示所有n個數據所預測的平均值。RMSE被用來衡量模型預測值與真實值的偏離程度,其值越接近于0則表示預測值越接近真實值,決定系數R2越接近于1則表示模型的擬合效果越好。

人工神經網絡在MBR的應用中具有諸多功能,其中一種功能是對MBR中水質進行預測。例如,Yuhang CAI等〔34〕運用小波神經網絡(Wavelet neural networks,WNN)對船舶廢水COD與氨氮去除率進行模擬,其R2分別為0.999與0.997,可為反應器的運行與調控提供良好的參考。此外,A. R. PENDASHTEH等〔35〕成功構建了一種處理含油高鹽廢水的水質預測模型,采用該模型對水質進行預測,結果表明,廢水中COD、總有機碳(Total organic carbon,TOC)和油脂的預測值與真實值較為接近,R2達到0.982。

ANN預測模型的另一個功能是預測膜污染并運用敏感性分析來識別影響膜污染的關鍵因素。例如A. R. ALKMIM等〔36〕運用多層感知器神經網絡(Mutilayer perception neural networks,MLPNN)對膜透水率進行預測,R2可以達到0.964,且計算得出混合液揮發性懸浮固體濃度(Mixed liquid volatile suspended solids,MLVSS)對膜透水率的變化貢獻最大。H. HAZRATI等〔37〕利用反向傳播神經網絡(Back propagation neural networks,BPNN)模型對TMP進行預測,其R2=0.999,且發現MLSS、HRT和運行時間3個變量均對TMP具有較大的影響。

2.4 人工神經網絡算法的分類及在膜污染中的應用

對圖2所統計的運用ANN算法研究膜污染的文獻進行分類并統計,結果見圖4。

1.jpg

圖4(a)根據不同用法將ANN算法劃分為簡單人工神經網絡、結合優化算法的神經網絡和深度學習神經網絡(Deep learning neural networks,DLNN)3類。圖4(b)展示了不同神經網絡文獻占比。由圖4(b)可知,在統計的82篇文獻中,82.93%的研究采用簡單神經網絡構建輸入變量與目標參數的模型,另有13.41%的研究采用了結合優化算法的神經網絡,而基于深度學習神經網絡的研究所占比例較小,不到5%。

筆者匯總了2016—2020年在MBR中運用不同人工神經網絡算法模擬膜污染的部分研究〔36-43〕,結果見表3。

1.jpg

目前在MBR膜污染預測中常用的簡單人工神經網絡主要有MLPNN、BPNN和徑向基函數神經網絡(Radical basis function neural networks,RBFNN)等。這些簡單神經網絡均包括輸入層、隱藏層、輸出層,但是又各具特點。MLPNN的不同層之間采用全連接的方式傳遞結果;BPNN的輸出結構采用前向傳播方式,誤差采用反向傳播方式進行;RBFNN則采用徑向基函數作為激活函數,它的輸出是輸入的徑向基函數和神經元參數的線性組合。大量研究證明基于簡單ANN的預測模型具有較高的預測精度,例如F. SCHMITT等〔39〕在AO-MBR中利用MLPNN預測TMP,并構建進水水質(總氮、總磷、硝酸鹽)與TMP的模型關系,其R2為0.850。盡管簡單人工神經網絡在膜污染預測方面具有較高的精度,但是仍存在擬合時間長、參數眾多且無法解決某些最優化問題等缺點。

在簡單神經網絡基礎上結合優化算法,能夠通過改善訓練方式來最小/最大化損失函數,具有加快神經網絡訓練速度、快速尋找參數最優值以及提高模型預測精度等優點。應用于膜污染預測的常見優化算法包括遺傳算法(Genetic algorithm,GA)和粒子群算法(Particle swarm optimization,PSO)等。其中,GA是基于生物體進化規律模型化問題的求解過程,具有優越的優化性能,被廣泛運用于機器學習和信號處理等領域〔44〕。S. A. MIRBAGHERI等〔45〕采用MLPNN預測TMP及膜透水率,并用GA算法對權值進行優化,結果證實優化后的模型具有更高的預測精度。PSO則是一種基于智能體的仿生優化算法,在求解全局最優解方面具有較大的優勢,可以完成神經網絡連接權值的訓練、結構設計和特征選擇等。劉志峰等〔46〕研究了一種PSO-BPNN,用粒子群算法代替傳統的梯度下降算法,結果表明優化后的模型具有更高的預測準確度,平均誤差從2.35%下降至0.83%。但是,結合優化算法模型的建立與參數優化過程需要耗費大量的時間,限制了其的廣泛應用。

此外,簡單ANN和結合優化算法的ANN的預測性能高度依賴于樣本量,并且在訓練過程中往往會陷入局部最優狀態,從而導致模型的預測穩定性較差。近年來,DLNN由于能夠有效彌補簡單ANN和結合優化算法ANN的不足,并挖掘數據的深度信息而被廣泛研究,其主要包括深度置信網絡(Deep belief network,DBN)、卷積神經網絡(Convolutional neural networks,CNN)、堆疊式自動編碼器(Stacked auto-encoders,SAE)網絡等。通常,DLNN是一種可以有效而靈活地表達高度變化的非線性的算法,具有強大的學習能力和非線性映射能力〔47〕。Honggui HAN等〔43〕開發了一種基于自組織深度置信網絡(Self-organizing deep belief network,SDBN)的膜透水率的預測方法,這種方法可以通過調整模型結構及參數提高模型預測精度,與簡單ANN相比,該模型具有更高的精度(RMSE為0.872)。在一些情況下,深度神經網絡的預測效果可能比結合優化算法的ANN更佳。Shuai SHI等〔48〕基于堆疊降噪自動編碼器(Stacked denoising auto-encoders,SDAE)深度學習神經網絡建立了MBR的性能預測模型,結果發現與GA-BPNN模型相比,其具有更高的R2及更低的均方誤差(Mean square error,MSE)。盡管深度學習神經網絡具有更高的預測精度,但是目前使用深度神經網絡在MBR中進行建模的研究相對較少,其原因可能是一方面深度學習神經網絡的復雜性較高,模型構建的時間遠超簡單ANN,另一方面使用簡單ANN已經可以滿足大多數的研究任務。

綜上所述,目前在MBR膜污染預測模型中,簡單神經網絡因高預測精度與高效性,在短時間內仍然是研究人員運用的主流工具。但是,隨著優化算法的運用以及深度神經網絡的開發,結合優化算法的ANN與深度學習ANN在面對更復雜環境下更具優勢,有望替代簡單神經網絡成為未來研究膜污染預測模型的主要工具。

2.5 ANN在中試MBR膜污染預測中的應用

雖然研究已經證實ANN在MBR膜污染預測上具有較高的性能,但是由于技術設備規模與廢水特性的差異,實驗室小試規模反應器的操作參數往往與實際的大規模反應器有明顯的區別,因此,推動ANN在中試規模MBR膜污染中的預測研究意義重大。Y. J. CHOI等〔49〕收集了中試規模的浸沒式MBR的長期運行數據,并構建了膜化學清洗后的通量、膜過濾時間及進水水質與膜透水率間的ANN模型,該模型R2可達0.950,實現了在中試規模下對MBR膜污染的精準預測。盡管如此,目前相關研究仍大多停留于小試階段,在針對實際廢水的中試MBR以及污水處理廠中應用依然較少。因此亟需根據中試及實際污水處理廠長期運行的數據建立ANN膜污染預測模型,從而更加精準地預測實際污水處理中MBR的膜污染狀況。

3 結論與展望

目前,膜污染仍然是阻礙MBR廣泛應用的瓶頸性問題,近年來,ANN模型因變量簡單、預測精度高等優點被用于MBR膜污染的預測。但是研究表明,ANN模型用于MBR膜污染預測研究仍存在以下缺陷,是未來研究的重要方向:

(1)ANN本質上屬于“黑箱”模型,目前由于神經網絡結構的特殊性難以對中間過程進行解釋。因此,未來的研究可以側重于通過解析神經網絡的復雜結構來厘清影響因子對膜污染的具體作用。

(2)ANN模型存在權值多、建模過程長等缺陷。因此未來膜污染預測模型的發展趨勢是優先使用更多優化的線性和非線性數據分析方法對數據進行降維或特征提取,從而減少數據的復雜程度,進而提高模型預測的精度并減少建模時間。

(3)目前ANN模型應用于中試和工業化規模的MBR膜污染預測較少,亟需基于長期實際運行的數據建立ANN膜污染預測模型,構建完備的預測體系,為調整污水處理廠運行操作條件提供技術支持,進而提高其運行效率、減緩膜污染并降低污水處理成本。




以上是小編收集整理的MBR中膜污染的人工神經網絡預測研究進展部分內容來自網絡,如有侵權請聯系刪除:153045535@qq.com;
本文地址:http://m.mirasorganicspa.com/shuichuli/1073.html



您可能還會對下面的文章感興趣:

天天操天天玩_19国产精品麻豆免费观看_欧美一级一级做性视频_国产欧美在线播放

        亚洲国产婷婷| 亚洲免费一级电影| 欧美在线视频一区| 欧美日韩国产色站一区二区三区| 国产亚洲成av人在线观看导航| 99亚洲视频| 美女视频黄a大片欧美| 国产亚洲精品aa| 亚洲午夜av| 欧美视频第二页| 亚洲精品日韩在线观看| 久久一区二区三区四区五区| 国产日韩欧美不卡| 国产美女搞久久| 9色porny自拍视频一区二区| 欧美不卡福利| 在线看成人片| 美女图片一区二区| 激情综合色综合久久| 久久成人免费网| 国产一区99| 久久er99精品| 国产视频精品xxxx| 久久国产精品久久国产精品| 国产日韩av一区二区| 久久成人资源| 韩日成人av| 老司机成人网| 亚洲福利久久| 欧美日韩不卡在线| 中国成人亚色综合网站| 欧美日韩在线观看一区二区| 一区二区欧美激情| 国产精品久久久久免费a∨| 亚洲在线视频网站| 国产精品日韩专区| 久久www免费人成看片高清| 国产日韩亚洲| 美女久久一区| aa级大片欧美| 国产美女一区二区| 久久久免费观看视频| 亚洲国产专区| 欧美色图一区二区三区| 性久久久久久久久| 在线不卡免费欧美| 欧美久久久久| 欧美伊人久久| 亚洲黄色有码视频| 欧美人与禽猛交乱配| 亚洲欧美一区二区原创| 在线成人小视频| 国产精品高潮视频| 久久综合九色综合欧美就去吻| 亚洲精品一级| 国产模特精品视频久久久久| 欧美1区2区3区| 正在播放欧美一区| 在线观看91精品国产入口| 欧美美女bb生活片| 久久精品国产亚洲高清剧情介绍| 亚洲国产美女精品久久久久∴| 欧美日韩中文在线| 久久精品在线播放| 国产精品99久久99久久久二8 | 欧美一级一区| 91久久综合亚洲鲁鲁五月天| 国产精品无码永久免费888| 狂野欧美一区| 欧美在线视频一区二区三区| 99re在线精品| 亚洲国产成人在线播放| 国产日韩一区二区| 欧美色图麻豆| 欧美理论在线播放| 美日韩精品免费| 欧美影院精品一区| 亚洲午夜激情网站| 亚洲人体偷拍| 一区二区三区在线免费播放| 国产精品一香蕉国产线看观看| 欧美激情精品久久久久久变态| 久久精品盗摄| 久久成人18免费观看| 亚洲在线免费观看| 一区二区三区精品久久久| 91久久精品网| 伊人久久婷婷色综合98网| 国产精品捆绑调教| 欧美成人网在线| 日韩视频在线你懂得| 黄色成人av在线| 国产精品久久二区| 免播放器亚洲一区| 久久精品亚洲一区二区三区浴池| 一区二区精品在线观看| 狠狠色狠色综合曰曰| 国内精品免费在线观看| 国产精品欧美久久| 欧美精品入口| 浪潮色综合久久天堂| 亚洲欧美中文日韩在线| 一本色道**综合亚洲精品蜜桃冫 | 性伦欧美刺激片在线观看| 日韩一级黄色av| 影音欧美亚洲| 亚洲第一精品影视| 一区二区三区无毛| 国产一区二区成人| 国产精自产拍久久久久久| 国产精品vip| 欧美日韩一区国产| 欧美人与禽猛交乱配视频| 麻豆成人在线观看| 免费不卡在线视频| 蜜臀a∨国产成人精品| 久久久国产成人精品| 欧美在线黄色| 午夜日韩激情| 久久免费少妇高潮久久精品99| 欧美综合国产精品久久丁香| 欧美一区国产二区| 久久久久久夜精品精品免费| 欧美综合国产精品久久丁香| 欧美一区免费| 久久激情网站| 久久综合综合久久综合| 免费人成精品欧美精品| 农夫在线精品视频免费观看| 欧美理论大片| 国产精品久久久久久福利一牛影视 | 伊人婷婷久久| 在线观看亚洲精品| 亚洲精品国产无天堂网2021| 9久草视频在线视频精品| 一区二区三区国产| 亚洲免费在线视频| 久久se精品一区精品二区| 另类天堂av| 欧美另类极品videosbest最新版本| 欧美喷潮久久久xxxxx| 欧美午夜精品理论片a级按摩| 国产精品一区二区三区乱码| 黑人巨大精品欧美黑白配亚洲| 激情亚洲网站| 中文在线不卡视频| 久久精品一区二区三区中文字幕| 久久免费黄色| 欧美性jizz18性欧美| 国产日韩欧美二区| 国产一区二区久久久| 一区二区三区四区五区精品| 欧美一区在线看| 欧美国产精品劲爆| 国产麻豆成人精品| 亚洲黄色在线看| 亚洲欧美精品| 欧美揉bbbbb揉bbbbb| 含羞草久久爱69一区| 一本久久综合亚洲鲁鲁| 久久嫩草精品久久久精品| 欧美日韩亚洲91| 在线播放国产一区中文字幕剧情欧美 | 亚洲一二三区在线观看| 久久久99久久精品女同性 | 欧美日韩国产小视频| 激情亚洲网站| 中文在线不卡| 欧美国产欧美亚洲国产日韩mv天天看完整| 欧美无乱码久久久免费午夜一区| 在线成人av| 亚洲欧美日韩国产精品| 欧美精品日韩精品| 在线欧美不卡| 久久福利一区| 国产精品免费久久久久久| 欧美日韩精品一区二区天天拍小说| 亚洲大片精品永久免费| 欧美亚洲综合网| 欧美日韩在线免费| 亚洲精品一区二区三区婷婷月 | 亚洲先锋成人| 欧美激情国产日韩精品一区18| 国产婷婷成人久久av免费高清| 性感少妇一区| 国产精品女同互慰在线看| 日韩视频亚洲视频| 麻豆久久精品| 激情小说另类小说亚洲欧美| 美国三级日本三级久久99| 黑人巨大精品欧美黑白配亚洲| 欧美一区日本一区韩国一区| 欧美三区视频| 日韩天堂在线观看| 欧美日韩综合视频网址| 一区二区三区黄色| 欧美视频一区二区三区| 日韩视频免费在线| 欧美大片一区二区三区| av成人老司机|