成人福利视频在线观看_国产精品日韩久久久久_欧美全黄视频_欧美网色网址

首頁> 資源> 論文>正文

活性污泥數(shù)學(xué)模型的研究應(yīng)用進(jìn)展與問題討論

論文類型 技術(shù)與工程 發(fā)表日期 2002-08-01
來源 《工業(yè)用水與廢水》2002年第4期
作者 李茹瑩,季民,任智勇,胡振苓,馬文杰
關(guān)鍵詞 活性污泥數(shù)學(xué)模型 模擬 參數(shù)校正
摘要 針對國際水質(zhì)協(xié)會(IAWQ)提出的活性污泥1號、2號、3號模型(ASM NO.1-3)的主要特點進(jìn)行了分析比較,指出了模型本身存在的若干問題。根據(jù)對活性污泥數(shù)學(xué)模型應(yīng)用研究的體會,認(rèn)為活性污泥數(shù)學(xué)模型應(yīng)用的難點是水質(zhì)的分析測定、模型簡化和參數(shù)校正。

李茹瑩1,季民1,任智勇1,胡振苓2,馬文杰2
(1.天津大學(xué) 環(huán)境科學(xué)與工程學(xué)院,天津 300072;2.天津市污水處理研究所,天津 300381)

  摘要:針對國際水質(zhì)協(xié)會(IAWQ)提出的活性污泥1號、2號、3號模型(ASM NO.1-3)的主要特點進(jìn)行了分析比較,指出了模型本身存在的若干問題。根據(jù)對活性污泥數(shù)學(xué)模型應(yīng)用研究的體會,認(rèn)為活性污泥數(shù)學(xué)模型應(yīng)用的難點是水質(zhì)的分析測定、模型簡化和參數(shù)校正。
  關(guān)鍵詞:活性污泥數(shù)學(xué)模型;模擬;參數(shù)校正
  中圖分類號:X703
  文獻(xiàn)標(biāo)識碼:A
  文章編號:1009-2455(2002)04-0004-03

Approaches to the Development and Problems of Research
and Application of Activated Sludge Model

LI Ru-ying1,JI Min1,REN Zhi-yong1,HU Zhen-ling2,MA Wen-jie2
(1.School of Environmental Science & Engineering,Tianjin University,Tianjin 300072,China;
2.Tianjin Wastewater Treatment Research Institute,Tianjin 30038l,China)

  Abstract: Analysis and comparison were made to the activated sludge model No.1,2 and 3 (ASM No.1,2
and 3) proposed by International Association on Water Quality (IAWQ).Some problems of the models were
found.According to the practice of applications and researches on activated sludge model, the key issues on modeling are water quality analysis,model simplification and parameter calibration.
  Key words:activated sludge model(ASM);modeling;parameter calibration

  近年來,國外有很多關(guān)于數(shù)學(xué)模型在飲用水處理、污泥處置、各種污水處理等工藝中的應(yīng)用報道。在活性污泥工藝眾多的數(shù)學(xué)模型中,由原國際水質(zhì)協(xié)會(IAWQ)(現(xiàn)已改稱為國際水協(xié)會IWA)推出的活性污泥數(shù)學(xué)模型(ASM)發(fā)展最為成熟,應(yīng)用最為廣泛。但我國在該領(lǐng)域的研究起步較晚,在實際應(yīng)用方面還相當(dāng)滯后。目前國內(nèi)僅有相當(dāng)有限的關(guān)于模型應(yīng)用的報道。本文重點討論IAWQ活性污泥數(shù)學(xué)模型的特點及目前存在的問題,并結(jié)合筆者在研究中的體會,分析其在應(yīng)用中的關(guān)鍵問題及其發(fā)展應(yīng)用前景。

1 IAWQ活性污泥數(shù)學(xué)模型

1.1 IAWQ活性污泥數(shù)學(xué)模型簡介
  為了鼓勵環(huán)境科學(xué)家和工程師更廣泛地把數(shù)學(xué)模型應(yīng)用到廢水生物處理系統(tǒng)的分析設(shè)計和運行管理中去,1983年,原國際水污染控制協(xié)會(IAWPRC)(國際水質(zhì)協(xié)會IAWQ的前身)組織了南非、丹麥、美國等五國專家組成活性污泥工藝模型課題組來完成活性污泥處理系統(tǒng)數(shù)學(xué)模型的研究。ASM課題組于1987年正式發(fā)表了技術(shù)報告,闡明了活性污泥1號模型(ASM No.1)的主要特性。它以矩陣的形式描述了污水中好氧、缺氧條件下所發(fā)生的水解、有機(jī)物降解、微生物生長、衰減等8種反應(yīng),模型中包含13種組分、5個化學(xué)計量系數(shù)和14個動力學(xué)參數(shù)[1]。
  ASM NO.1自推出以來得到廣泛應(yīng)用,但它的缺陷是未包含磷的去除。針對此問題,IAWQ專家組于1995年又推出活性污泥2號模型(ASM No.2),它包含了磷的吸收和釋放,增加了厭氧水解、酵解及與聚磷菌有關(guān)的4個反應(yīng)過程。因為生物除磷機(jī)理很復(fù)雜,所以ASM No.2非常龐大,它包含19種物質(zhì)、19種反應(yīng)、22個化學(xué)計量系數(shù)以及42個動力學(xué)參數(shù)[2]。該模型提出了包含化學(xué)需氧量(COD)、氮和磷去除過程在內(nèi)的綜合性生物處理工藝過程動態(tài)模擬理論,它不是生物除磷模型的最終方案,而是一種折中方案。
  1999年IAWQ專家組經(jīng)過對1號模型應(yīng)用中問題的修正,推出活性污泥3號模型(ASM No.3)。ASM No.3不以水解作用為重點、引入有機(jī)物在微生物體內(nèi)的貯藏及內(nèi)源呼吸,以強(qiáng)調(diào)細(xì)胞內(nèi)部的活動過程。ASM No.3與ASM No.1中主要現(xiàn)象是相關(guān)的,如以城市污水為主的活性污泥系統(tǒng)中的氧消耗、污泥產(chǎn)量、硝化和反硝化,但ASM No.2中的生物除磷在此不予考慮。ASM No.3包括13個組分,12個生化過程。ASM No.3僅考慮微生物轉(zhuǎn)化過程,不包括化學(xué)沉淀,但基于ASM No.2(包含磷的吸收和釋放過程)所提供的信息,很容易加入該過程[3]。
1.2 活性污泥1、2、3號模型(ASM No.1~3)的比較
  ASM諸多版本的共同特點是它們均以Monod方程為基礎(chǔ),都是多維的并包含大量的動力學(xué)參數(shù)和化學(xué)計量系數(shù),均以矩陣的形式描述生物反應(yīng)過程。
  ASM No.1和ASM No.2在應(yīng)用過程中都有一些限制,如他們都要求pH值接近中性并保持恒定;ASM No.1要求系統(tǒng)在恒定溫度下運行,ASM No.2的實用溫度要限制在中等范圍,大概為10~25℃,因為高溫和低溫狀態(tài)下聚磷菌PAOS的特性尚不完全清楚,模型不一定能給出合理的預(yù)測,尤其是對于磷的去除。從ASM No.1到ASM No.2的最突出變化是生物量按反應(yīng)過程進(jìn)行了更為詳細(xì)的劃分,使其濃度不能簡單地用分布參數(shù)XB,M。描述。除了生物除磷過程外,ASM No.2還包含了兩個“化學(xué)過程”,可以用于模擬磷的化學(xué)沉淀。
  ASM No.3與ASM No.1大致相同,只是在廢水特征化這一重要方面作了改動,將重點從水解過程轉(zhuǎn)移到有機(jī)物質(zhì)的儲存過程。在ASM No.1中快速可生物降解基質(zhì)COD必須從呼吸試驗中估算,而對這一試驗的解釋又依賴于異養(yǎng)菌產(chǎn)率YH的值。在ASM No.3中溶解性COD僅由快速可生物降解基質(zhì)SS和惰性可溶性有機(jī)質(zhì)SI組成。從兩個模型所采用的典型廢水組成可以看出:ASM No.3中SS占總COD的40%,而不象ASM No.1中僅為10%。
  與ASM No.1和ASM No.2相比,ASM No.3的一個重要區(qū)別是,通過0.45μm濾膜的過濾作用可以將溶解性成分與顆粒性組分更好地區(qū)分,而前者廢水游液中仍然含有相當(dāng)比例的慢速可生物降解基質(zhì)XS

2 IAWQ活性污泥數(shù)學(xué)模型目前存在的問題

  綜合國外近年來對模型的應(yīng)用研究,已證實了ASM No.1對于描述碳氧化、硝化和反硝化過程具有較強(qiáng)的模擬預(yù)測能力。但是經(jīng)過多年的實際應(yīng)用,也發(fā)現(xiàn)ASM No.1中存在著一些明顯缺點[3]
 ?、貯SM No.1中不包含氮和堿度對異養(yǎng)菌的限制因素,這樣在某些情況下會出現(xiàn)負(fù)濃度(如銨鹽等)。
  ②ASM No.1中包括了可生物降解的溶解性有機(jī)氮(SNS)和顆粒性有機(jī)氮(XNS)。在實踐中,二者都不易測定,應(yīng)將其從模型中刪除。
  ③ASM No.1中氨化反應(yīng)的動力學(xué)難于定量化。在應(yīng)用ASM No.1時,假定所有有機(jī)組分均為恒定組成,即恒定的N與COD質(zhì)量比。
  ④ASM No.1對顆粒性惰性有機(jī)物質(zhì)按其來源將其區(qū)分為兩類,即進(jìn)水所含的XI和生物質(zhì)衰減產(chǎn)生的XD,實際中,難于如此清晰地區(qū)分這兩種組分。因此ASM No.3中將XD和原ASM No.1中的XI統(tǒng)一考慮為XI處理。
  ⑤在模型中細(xì)菌的死亡分解與水解和增長過程分為有機(jī)物儲存,死亡,捕食和分解等階段,從而使得相應(yīng)動力學(xué)參數(shù)的確定變得非常困難。
 ?、奕绻鬯幚韽S中易降解基質(zhì)濃度較高時,在好氧和缺氧條件下可以觀測到聚羥基烷酸PHAS(poly-hydroxy-alkanoates)(有時為糖原質(zhì)-glyco-gen)的儲存現(xiàn)象,但ASM No.1中不包括這一過程。
 ?、逜SM No.1認(rèn)為好氧與缺氧條件下硝化菌的衰減速率是相同的,當(dāng)固體停留時間(SRT)較長或缺氧池體積比例較高時,對最大硝化反應(yīng)速率的預(yù)測就會出現(xiàn)問題。
  除磷過程的加入給ASM No.2模型的使用也帶來了諸多限制。例如模型沒有考慮鉀和鎂對生物除磷的限制作用。但是眾所周知,鉀和鎂是PAOS中構(gòu)成聚磷酸鹽的兩種重要陽離子。這些陽離子的短缺會導(dǎo)致污泥中聚磷菌聚積作用的惡化,從而導(dǎo)致磷去除的明顯降低。據(jù)報道,亞硝酸鹽和一氧化氮對生物除磷有抑制作用,但在模型中沒有考慮這種影響。由于ASM No.2的復(fù)雜性和除磷機(jī)理的不確定性,使得ASM No.2相較于ASM No.1的應(yīng)用限制較多,經(jīng)驗也不如ASM No.1成熟。但它是目前唯一包含磷的去除的較為成功的活性污泥數(shù)學(xué)模型。
  由于ASM No.3提出較晚,目前還沒有經(jīng)過大量的運行數(shù)據(jù)的驗證,所以仍然需要在實踐中對模型進(jìn)行不斷地修正和改進(jìn),特別是對儲存現(xiàn)象的描述。

3 IAWQ活性污泥數(shù)學(xué)模型應(yīng)用過程中幾個問題的討論

  活性污泥數(shù)學(xué)模型的研究和應(yīng)用始終是國際上廢水生物處理領(lǐng)域研究的熱點之一,目前已有許多大型污水處理廠在設(shè)計、改建和運行過程中應(yīng)用了數(shù)學(xué)模型,積累了許多寶貴的經(jīng)驗。從國內(nèi)外大量的應(yīng)用實例中不難看出,IAWQ模型應(yīng)用中問題的焦點在于模型的簡化、水質(zhì)的分析和測定以及參數(shù)的校正。下面筆者結(jié)合國內(nèi)外應(yīng)用報道及自身研究過程中的切身體會對這幾方面問題分別進(jìn)行討論。
3.1 水質(zhì)的分析和測定
  由于模型中涉及的組分較多,而且不能由常規(guī)水質(zhì)分析指標(biāo)直觀表述,需要設(shè)定常規(guī)檢測數(shù)據(jù)與模型參數(shù)的轉(zhuǎn)換方案。如考慮污水除磷脫氮則要測定進(jìn)水構(gòu)成中COD、氮、磷的溶解性和非溶解性、惰性和非惰性物質(zhì)各自的含量,測定的內(nèi)容比目前常規(guī)的水質(zhì)化驗多出很多,而且具體的測定方法還未規(guī)范化,有些還有待于進(jìn)一步研究解決。IAWQ的專家組目前正致力于對進(jìn)水水質(zhì)測定方法簡單化、規(guī)范化方面的研究[4-5]。
3.2 模型的簡化
  由于IAWQ數(shù)學(xué)模型的復(fù)雜性,將其完整的模型直接應(yīng)用于污水廠的設(shè)計、模擬、運行控制是不現(xiàn)實的,這就需要根據(jù)具體條件將其適當(dāng)簡化。但值得環(huán)境工作者注意的是,模型的簡化雖然是必要的,但也是相當(dāng)危險的。如果簡化不當(dāng),將導(dǎo)致應(yīng)用的徹底失敗。因此,在模型應(yīng)用于某一特定污水廠之前,必須充分了解污水水質(zhì)特性、工藝的水力學(xué)特性以及生物反應(yīng)過程和機(jī)理。在此基礎(chǔ)上,分析各變量和參數(shù)對系統(tǒng)模擬的敏感性,忽略不敏感因素,僅考慮對 系統(tǒng)影響較大的因素,從而達(dá)到簡化的目的。
  模型簡化的結(jié)果是矩陣中變量和生物反應(yīng)過程及其相關(guān)參數(shù)的減少,或是非線性模型的線性化。簡化后的模型必須經(jīng)過校正方可使用,而校正過程的難點便集中于參數(shù)的調(diào)整。
3.3 參數(shù)的校正
  簡化的模型必須進(jìn)行參數(shù)校正,對于活性污泥系統(tǒng)這樣的復(fù)雜非線性系統(tǒng),參數(shù)調(diào)整有兩種方法,一種是手動調(diào)整,一種是自動調(diào)整。手動調(diào)整耗時長且結(jié)果重現(xiàn)性差,故僅在無法自動調(diào)整時采用。但是自動調(diào)整也不是可以直接使用,它要求先進(jìn)的在線監(jiān)測儀器和功能很強(qiáng)的數(shù)學(xué)工具。國外曾討論使用界面響應(yīng)法(RSM)進(jìn)行參數(shù)調(diào)整[6]。該方法是用于建立、改進(jìn)、優(yōu)化過程的統(tǒng)計學(xué)和數(shù)學(xué)技術(shù)的結(jié)合。在RSM中,響應(yīng)值(輸出值)與輸入變量的濃度有關(guān)。在活性污泥系統(tǒng)數(shù)學(xué)模型的校正中,RSM用作局部敏感分析工具,其中輸入變量是模型參數(shù),輸出是余差平方和。在荷蘭的Rotterdam污水處理廠的卡羅塞型氧化溝脫氮模型(ASM No.1)中,采用該方法進(jìn)行參數(shù)調(diào)整,獲得了精確可靠的結(jié)果??梢灶A(yù)見,隨著計算機(jī)技術(shù)的飛速發(fā)展,更加快速方便的校正方法必將應(yīng)運而生,從而推動活性污泥數(shù)學(xué)模型更加廣泛的應(yīng)用。

4 結(jié)語

  IAWQ活性污泥數(shù)學(xué)模型在國外的大量應(yīng)用中已取得了寶貴的經(jīng)驗,污水處理工藝的不斷完善和生物處理機(jī)理的不斷明確,以及計算機(jī)技術(shù)的飛速發(fā)展必將使數(shù)學(xué)模型得到更為廣泛的應(yīng)用。
  我國已建及在建的廢水處理廠中大約有80%以上采用活性污泥法,應(yīng)用數(shù)學(xué)模型設(shè)計和管理廢水處理廠是提高我國廢水處理廠設(shè)計、運行和管理水平的必然需要。但數(shù)學(xué)模型的應(yīng)用需要較高的硬件水平,國外污水處理廠大都采用在線監(jiān)測實時控制,而我國污水處理廠水質(zhì)分析大都采用現(xiàn)場取樣、化驗室分析的方法,使得測定結(jié)果引入了大量人為誤差,而且測定值僅是每天某個時刻的瞬時值,不能全面反映水質(zhì)變化情況,因而模擬值與測量值之間誤差較大。相信隨著我國水處理工藝不斷改進(jìn)、硬件水平不斷改善,數(shù)學(xué)模型在我國污水處理中的應(yīng)用必將得到推廣。

參考文獻(xiàn)
[1]M Henze,C P L Grady,W Gujer,G V R Marais,T Matsuo.IAWPRC scientific and technical reports No 1:Activated sludge moddel No1[R].IAWPRC:IAWPRC Task Group on Mathematical Modeling for design and operation of biological wastewater treatment,1986.
[2]W Gujer,M Henze,Takahashi Mino,T Matsuo,et al.Activated sludge model No 2:biological phosphorus removal[J].Wat Sci Tech,1995,3(2):1~11.
[3]Willi Gujer,Mogens Henze,Takahashi Mino,et al.Activated sludeg model No 3[J].Wat Sci Tech,1999,39(1):183~193.
[4]Nowak O,F(xiàn)ranz A,Svardal,Muller V,Kuhn V.Parameter estimation for activated sludge models with the help of mass balances[J].Wat Sci Tech,1999,39(4):113~12O.
[5]王維斌.基于神經(jīng)網(wǎng)絡(luò)的活性污泥建模及其應(yīng)用研究[D].天津;天津大學(xué),2001.
[6]A Abusam,K J Keesman,G van Straten.Parameter estimation procedure for complex non-linear systems:calibration of ASM No.1 for N-removal in a full-scale oxidation ditch[J].Wat Sci Tech,2001,43(7):357~365.


作者簡介
李茹瑩(1977~),女,遼寧鞍山人,天津大學(xué)碩士研究生,主要從事水處理研究,電話(022)27401644,ramyli@eyou.com。

論文搜索

發(fā)表時間

月熱點論文

論文投稿

很多時候您的文章總是無緣變成鉛字。研究做到關(guān)鍵時,試驗有了起色時,是不是想和同行探討一下,工作中有了心得,您是不是很想與人分享,那么不要只是默默工作了,寫下來吧!投稿時,請以附件形式發(fā)至 paper@h2o-china.com ,請注明論文投稿。一旦采用,我們會為您增加100枚金幣。

成人福利视频在线观看_国产精品日韩久久久久_欧美全黄视频_欧美网色网址
美国精品在线观看| 日韩福利电影在线| 精品一区二区三区视频在线观看| 色妞www精品视频| 日韩三级伦理片妻子的秘密按摩| 亚洲欧美另类综合偷拍| 国产成人精品免费在线| 日韩欧美在线观看一区二区三区| 美女mm1313爽爽久久久蜜臀| 日韩美女在线视频| 日韩精品国产欧美| 欧美三电影在线| 一区二区三区产品免费精品久久75| 成人精品国产一区二区4080| 精品成人一区二区三区| 日本色综合中文字幕| 日韩欧美不卡一区| 国产成人亚洲精品青草天美| 国产精品嫩草99a| 成人一区二区三区中文字幕| 久久久久国产精品麻豆| 精品综合久久久久久8888| 久久久精品国产免大香伊| 99久久99精品久久久久久| 国产精品丝袜久久久久久app| 99久久婷婷国产综合精品电影| 亚洲视频在线一区二区| 99精品国产视频| 亚洲成人精品一区二区| 欧美日韩精品欧美日韩精品| 亚洲国产另类精品专区| 在线观看日韩精品| 一区二区三区不卡视频在线观看| 欧美日韩mp4| 国产乱子伦一区二区三区国色天香 | 国产**成人网毛片九色| 久久午夜电影网| 精品一区二区三区免费观看| 欧美国产成人在线| 成人av手机在线观看| 一区二区欧美视频| 欧美mv日韩mv国产网站app| 极品销魂美女一区二区三区| 中文字幕一区二区三区在线观看| 99在线精品观看| 亚洲青青青在线视频| 91精品国产乱码| 另类欧美日韩国产在线| 久久综合久久综合亚洲| 91丨九色丨蝌蚪丨老版| 一区二区激情视频| 精品国产第一区二区三区观看体验 | 在线综合亚洲欧美在线视频| 六月丁香综合在线视频| 国产精品的网站| 色婷婷精品久久二区二区蜜臂av| 奇米影视一区二区三区| 久久综合资源网| 在线观看中文字幕不卡| 国产乱一区二区| 亚洲一区二区三区四区在线观看| 91精品国产一区二区| 成人h精品动漫一区二区三区| 日韩黄色在线观看| 自拍偷在线精品自拍偷无码专区| 欧美不卡一区二区三区四区| 91污在线观看| 经典三级一区二区| 午夜久久久影院| 精品国产乱子伦一区| 色综合夜色一区| 日本中文在线一区| 日韩美女久久久| 久久久午夜精品| 欧美日韩日日摸| 91在线码无精品| 国产一区不卡精品| 日韩成人一区二区三区在线观看| 综合久久给合久久狠狠狠97色| 精品第一国产综合精品aⅴ| 欧美日韩精品是欧美日韩精品| 99re热视频这里只精品| 国产一区二区h| 美女脱光内衣内裤视频久久影院| 亚洲最大的成人av| 中文无字幕一区二区三区| 精品美女一区二区| 在线成人小视频| eeuss国产一区二区三区| 精品在线你懂的| 日韩激情一区二区| 亚洲地区一二三色| 亚洲精品免费在线观看| 精品国产乱码久久久久久夜甘婷婷 | 国产精品国模大尺度视频| xnxx国产精品| 日韩精品影音先锋| 欧美久久久一区| 欧美性生活久久| 国产91综合一区在线观看| 久久99精品国产麻豆婷婷| 日韩成人精品在线| 午夜精品视频一区| 亚洲福利视频导航| 一区二区三区中文字幕| 综合自拍亚洲综合图不卡区| 中文字幕在线观看不卡| 欧美国产视频在线| 国产亚洲一区字幕| 日韩一级片网址| 91 com成人网| 欧美美女一区二区三区| 欧美无砖专区一中文字| 在线视频你懂得一区二区三区| 99精品在线免费| 99国产欧美另类久久久精品| 99免费精品视频| 99国内精品久久| 一本到一区二区三区| 91视频.com| 色婷婷综合中文久久一本| 一本到三区不卡视频| 欧美主播一区二区三区美女| 色综合天天天天做夜夜夜夜做| 91热门视频在线观看| 色综合天天在线| 欧美视频自拍偷拍| 欧美精品日韩综合在线| 欧美一区二区三区在线观看视频 | 国产精品污网站| 国产精品乱码人人做人人爱| 一区在线观看免费| 亚洲乱码中文字幕综合| 亚洲资源在线观看| 天天操天天干天天综合网| 日本成人超碰在线观看| 久久99精品国产麻豆婷婷 | 欧美日韩日日摸| 日韩欧美在线观看一区二区三区| 精品国产污网站| 2021国产精品久久精品| 欧美高清一级片在线观看| 亚洲国产精品黑人久久久 | 亚洲日本丝袜连裤袜办公室| 亚洲与欧洲av电影| 日韩av在线播放中文字幕| 免费观看久久久4p| 国产在线观看一区二区| 国产成人夜色高潮福利影视| 91视频www| 欧美丰满少妇xxxbbb| www欧美成人18+| 国产精品久久久一区麻豆最新章节| 亚洲欧美日韩在线| 日日夜夜免费精品| 国产伦理精品不卡| 91在线视频免费观看| 欧美二区在线观看| 国产视频911| 亚洲男人的天堂在线aⅴ视频| 亚洲成精国产精品女| 精品在线观看免费| 99这里只有精品| 欧美高清视频www夜色资源网| 精品国产一区二区精华| 中文字幕亚洲区| 午夜精品免费在线| 国产成人精品影视| 在线观看国产一区二区| 日韩欧美国产午夜精品| 国产精品久久久一区麻豆最新章节| 午夜电影久久久| 国产成人一区在线| 欧美日韩国产小视频| 国产亚洲综合色| 亚洲午夜久久久久久久久电影网 | 美女任你摸久久| 成人黄色大片在线观看| 欧美二区在线观看| 国产精品福利av| 日本不卡的三区四区五区| 成人午夜视频在线观看| 欧美高清视频一二三区| 国产精品国产精品国产专区不片 | 视频精品一区二区| 粉嫩嫩av羞羞动漫久久久 | 美女视频第一区二区三区免费观看网站| 成av人片一区二区| 欧美一区二区三区精品| 国产精品第13页| 九色综合狠狠综合久久| 91国产视频在线观看| 久久你懂得1024| 亚洲大片在线观看| 99麻豆久久久国产精品免费 | 精品免费国产一区二区三区四区| 亚洲狠狠丁香婷婷综合久久久| 国产乱子伦一区二区三区国色天香| 欧美三区免费完整视频在线观看| 国产日韩欧美一区二区三区综合|