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

    基于高通量測(cè)序的微生物強(qiáng)化污泥減量工藝中微生物群落解析

    2016-10-13 21:37:26宋云龍張金松邵明非鄧仁健
    中國(guó)環(huán)境科學(xué) 2016年7期
    關(guān)鍵詞:氧區(qū)桿菌屬活性污泥

    宋云龍,張金松*,朱 佳,邵明非,鄧仁健

    ?

    基于高通量測(cè)序的微生物強(qiáng)化污泥減量工藝中微生物群落解析

    宋云龍1,張金松1*,朱 佳2,邵明非1,鄧仁健1

    (1.哈爾濱工業(yè)大學(xué)深圳研究生院,土木與環(huán)境工程學(xué)院,廣東 深圳 518055;2.深圳職業(yè)技術(shù)學(xué)院建筑與環(huán)境工程學(xué)院,廣東 深圳 518055)

    利用高通量測(cè)序技術(shù)對(duì)微生物強(qiáng)化原位污泥減量工藝中微生物群落進(jìn)行解析.結(jié)果表明,外源菌劑的投加改變了活性污泥中微生物群落結(jié)構(gòu),菌劑的有效成分乳桿菌屬和醋酸桿屬在強(qiáng)化組中含量顯著增加.科水平上的群落多樣性分析顯示:強(qiáng)化組厭氧區(qū)Siompson指數(shù)、Shannon指數(shù)和Pielou指數(shù)均不同程度升高,群落多樣性增加.聚類分析顯示,微生物群落按時(shí)間順序明顯聚為3簇,強(qiáng)化組厭氧區(qū)微生物群落隨時(shí)間發(fā)生較大演變.主坐標(biāo)分析顯示,強(qiáng)化組和對(duì)照組微生物群落明顯聚為2類,其厭氧區(qū)群落分別按照不同的方向演變,微生物強(qiáng)化和DO降低的協(xié)同作用是強(qiáng)化組群落演變驅(qū)動(dòng)力,DO降低是對(duì)照組群落演變驅(qū)動(dòng)力.典范對(duì)應(yīng)分析顯示顯著影響微生物群落結(jié)構(gòu)的環(huán)境因子依次為pH值、水溫和DO.

    高通量測(cè)序;微生物強(qiáng)化;水解酸化;群落解析

    微生物強(qiáng)化技術(shù)在污水處理系統(tǒng)中主要用于去除難降解有機(jī)物[1-4],近年來(lái)由于我國(guó)污泥問(wèn)題日益嚴(yán)峻[5-7],研究者嘗試將微生物強(qiáng)化技術(shù)用于污泥的處理處置[8-9].微生物強(qiáng)化原位污泥減量作為一種新興、高效的污泥減量技術(shù)受到廣泛關(guān)注,并取得了顯著的污泥減量效果[10-12],但是由于其減量機(jī)理尚不明確,因而尚未得到大范圍內(nèi)推廣應(yīng)用.現(xiàn)有的微生物強(qiáng)化污泥減量機(jī)理分析多停留在理化性質(zhì)層面,而微生物強(qiáng)化對(duì)微生物群落結(jié)構(gòu)的影響方面的研究較少.現(xiàn)有微生物群落分析多是依靠PCR-DGGE技術(shù),僅通過(guò)電泳圖譜中的幾十個(gè)條帶信息對(duì)微生物群落多樣性進(jìn)行定性或半定量的分析, 而DGGE技術(shù)對(duì)于群落中相對(duì)豐度低于1%的非優(yōu)勢(shì)物種檢測(cè)效果甚微[13],難以從多個(gè)分類層對(duì)物種豐度和群落特征進(jìn)行更全面深入的信息挖掘[14-16].以Illumina測(cè)序?yàn)榇淼牡?代高通量測(cè)序具有測(cè)序成本低、信息采集量大、適用性廣等優(yōu)點(diǎn),成為近年來(lái)應(yīng)用最廣泛的測(cè)序技術(shù)[17-20],為深度解析活性污泥微生物群落提供了有力的支持.

    課題組前期通過(guò)投加菌劑的方式對(duì)A/O工藝活性污泥進(jìn)行微生物強(qiáng)化作用,取得了30%左右的污泥減量效果[21].在此基礎(chǔ)上,進(jìn)一步依托高通量測(cè)序平臺(tái)Illumina-HiSeq2000,研究了低污泥濃度工況中外源菌劑的投加對(duì)活性污泥微生物群落結(jié)構(gòu)的影響,從群落演變的層面探討了A/O工藝中微生物強(qiáng)化原位污泥減量的機(jī)理,為微生物強(qiáng)化原位污泥減量技術(shù)的規(guī)?;茝V應(yīng)用提供理論依據(jù).

    1 材料與方法

    1.1 A/O中試試驗(yàn)

    1.1.1 中試裝置 裝置位于深圳市羅芳污水處理廠,整個(gè)系統(tǒng)為兩套完全平行且獨(dú)立運(yùn)行的A/O工藝中試裝置,模擬羅芳污水處理廠一期處理工藝,單臺(tái)設(shè)備處理規(guī)模240m3/d,構(gòu)造如圖1所示,其厭氧1區(qū)、厭氧2區(qū)、好氧1區(qū)、好氧2區(qū)和沉淀池有效容積分別為12.5,25,22,22,30m3.

    1.1.2 試驗(yàn)污水污泥 中試進(jìn)水取自深圳羅芳污水處理廠一期工程的中間沉淀區(qū),具體水質(zhì)如下:pH值為6.25~7.68,均值7.03;COD為102~ 334mg/L,均值155mg/L;BOD5為47~135mg/L,均值67mg/L;總氮為21~52mg/L,均值32mg/L;氨氮為13~42mg/L,均值21mg/L;總磷為1.59~ 6.97mg/L,均值3.25mg/L.試驗(yàn)污泥取自深圳羅芳廠一期AA/O工藝好氧區(qū).

    1.1.3 菌劑投加及運(yùn)行方式 菌劑投加量為日處理水量(240m3/d)的萬(wàn)分之四(96kg),投加點(diǎn)為厭氧2區(qū),一次性投加.投加微生物菌劑的稱強(qiáng)化組,未投加微生物菌劑稱對(duì)照組.中試運(yùn)行條件模擬羅芳污水處理廠的生產(chǎn)工況:反應(yīng)體系MLSS維持在3500~4000mg/L,好氧區(qū)DO維持在3~5mg/L,厭氧區(qū)DO維持在0.05~0.15mg/L,污泥回流比為1~1.2,強(qiáng)化組和對(duì)照組污泥齡分別為28.04和17.99d,兩組HRT均為8.15h,在保證出水達(dá)到GB18918-2002[22]規(guī)定的一級(jí)A標(biāo)準(zhǔn)的前提下盡量減少排泥.

    1.2 試驗(yàn)方法

    1.2.1 取樣 取強(qiáng)化組和對(duì)照組厭氧2區(qū)的第3,10,21d的活性污泥,相應(yīng)標(biāo)記為d3-Y1、d3-Y2、d10-Y1、d10-Y2、d21-Y1、d21-Y2,取強(qiáng)化組和對(duì)照組好氧2區(qū)的第5、21d的活性污泥,相應(yīng)標(biāo)記為d5-H1、d5-H2、d21-H1、d21-H2.活性污泥樣品經(jīng)0.2mm玻璃纖維濾膜過(guò)濾后-80℃保存.

    1.2.2 檢測(cè)方法 用FastDNA?土壤試劑盒(MP Biomedicals,Illkirch,France)進(jìn)行細(xì)菌組DNA提取.16S rRNA基因V6區(qū)序列采用引物對(duì)967F (CAACGCGAAGAACCTTACC)與1046R (CGACAGCCATGCANCACCT)進(jìn)行擴(kuò)增[23]. PCR條件:98℃變性2min,然后28個(gè)循環(huán)(98℃20s,55℃20s, 68℃1min),最后在68℃延伸5min.為減少PCR擴(kuò)增前期的潛在錯(cuò)誤,每DNA樣品同步進(jìn)行2個(gè)平行擴(kuò)增,最后再進(jìn)行混合與后續(xù)分析.擴(kuò)增后采用2.0%的瓊脂糖凝膠電泳對(duì)擴(kuò)增產(chǎn)物進(jìn)行驗(yàn)證.最后,PCR擴(kuò)增產(chǎn)物采用純化試劑盒(天根,北京)進(jìn)行純化,并采用光度法(NanoDrop-1000)測(cè)定產(chǎn)物濃度后混合樣品,保證各樣品核酸濃度基本相同.最后采用Illumina- HiSeq2000(深圳華大基因)高通量測(cè)序平臺(tái)測(cè)序后進(jìn)行群落解析.

    1.2.3 多樣性指數(shù)的計(jì)算 辛普森多樣性指數(shù)(Simpson's diversity index,)[24]:

    式中:p表示種的個(gè)體在群落中的比例;p2表示隨機(jī)取2個(gè)個(gè)體為同種的概率.

    香濃多樣性指數(shù)(Shannon-Wiener,¢)[25]:

    式中:表示總的物種數(shù);p表示第個(gè)種占總數(shù)的比例.

    Pielou均勻度指數(shù)(Pielou's evenness index,)[26]:

    式中:¢為香農(nóng)指數(shù);是¢的最大值.

    2 結(jié)果與討論

    2.1 微生物強(qiáng)化污泥減量效果及污水處理效能

    表1 微生物強(qiáng)化污水處理效果 Table 1 Pollutants removal efficiency by bioaugmentation

    中試裝置有效運(yùn)行了30d,取得了良好的污泥減量效果.強(qiáng)化組和對(duì)照組的MLSS產(chǎn)率分別為0.36,0.49kgMLSS/kgCOD,MLVSS產(chǎn)率分別為0.19,0.26kgMLSS/kgCOD,MLSS和MLVSS減量率分別為26.24%和28.09%.微生物強(qiáng)化污水處理效果見(jiàn)表1,可以看出,投加菌劑后,強(qiáng)化組SS、COD、BOD5、TN、氨氮的去除效果得到不同程度的加強(qiáng),但TP的去除能力略微下降.強(qiáng)化組和對(duì)照組均出水達(dá)到GB18918-2002一級(jí)A標(biāo)準(zhǔn).

    2.2 高通量測(cè)序?qū)ξ⑸锶郝涞蔫b定效率

    本次測(cè)序各樣品原始序列數(shù)分布比較均勻,除D21-H2較高之外,基本在10000條左右(表2).對(duì)原始DNA序列去除了古菌序列,并刪除了錯(cuò)誤序列,最終將各樣品的DNA序列數(shù)統(tǒng)一為7341條.之后采用Uclust方法以及0.97的相似度對(duì)序列進(jìn)行聚類,針對(duì)聚類獲得的OTU采用First方法抽取代表性序列,將代表序列與Greengene數(shù)據(jù)庫(kù)進(jìn)行比對(duì),完成對(duì)各OTU的注釋.從而將各樣本中每條序列鑒定到門、綱、目、科、屬、種等不同的分類水平.

    表2 各樣品測(cè)序深度 Table 2 Sequence depth of each sample

    注:可鑒定序列是指本次測(cè)序測(cè)得序列中與Greengene數(shù)據(jù)庫(kù)已知物種序列相似度不低于97%的序列,低含量物種是指比例低于1%的可鑒定序列所對(duì)應(yīng)的物種.

    本次測(cè)序在門、綱、目、科、屬、種水平上的鑒定率分別是98.76%、94.58%、87.56%、82.65%、51.14%和0.69%,在門、綱、目、科水平鑒定結(jié)果良好,對(duì)屬水平也有一定鑒定效率,對(duì)種水平鑒定能力較差(表3).

    2.3 微生物強(qiáng)化對(duì)群落組成的影響

    表4列舉了A/O中試裝置活性污泥中的優(yōu)勢(shì)微生物,按照相對(duì)豐度降序排列.厭氧區(qū)和好氧區(qū)在門水平上的優(yōu)勢(shì)微生物組成相種類相同,主要有變形菌門(Betaproteobacteria, Alphaproteobacteria, Gammaproteobacteria, Deltaproteobacteria)、放線菌門(Actinobacteria)、擬桿菌門(Bacteroidetes)、浮霉菌門(Planctomycetes)、酸桿菌門(Acidobacteria)和綠彎菌門(Chloroflexi),且按照豐度排序基本一致;在綱、目水平上,厭氧區(qū)和好氧區(qū)優(yōu)勢(shì)微生物組成相同,但按照豐度排序差異較大;在科、屬水平上厭氧區(qū)與好氧區(qū)微生物組成種類基本相似,腐螺旋菌科(Saprospiraceae)和亞硝化單胞菌科(Nitrosomonadaceae)分別是厭氧區(qū)和好氧區(qū)的優(yōu)勢(shì)科,索氏均屬(Thauera)和苯基桿菌屬(Phenylobacterium)分別是厭氧區(qū)和好氧區(qū)的優(yōu)勢(shì)屬;厭氧區(qū)和好氧區(qū)的各科屬微生物豐度差異較大.

    注:( )內(nèi)為相對(duì)豐度.

    微生物菌劑中最多的是厚壁門(Firmicutes),占細(xì)菌總數(shù)的53.48%,其次是變形門(Proteobacteria),占細(xì)菌總數(shù)的46.38%,其中對(duì)污泥水解酸化具有強(qiáng)化作用的主要是乳桿菌屬()和醋酸桿菌屬().乳桿菌屬屬于厚壁門(Firmicutes)、桿菌綱(Bacilli)、乳桿菌目(Lactobacillales)、乳桿菌科(Lactobacillaceae),醋酸桿菌屬屬于厚壁門(Firmicutes)、梭菌綱(Clostridia)、梭菌目(Clostridiales)、真桿菌科(Eubacteriaceae).乳桿菌屬是異養(yǎng)菌,耐酸性,能發(fā)酵分解糖代謝,終產(chǎn)物中50%以上是乳酸[27],在污水污泥處理中有廣泛的應(yīng)用.Vellingiri等[28]利用馴化培養(yǎng)的嗜酸乳桿菌降解生活污水中總懸浮固體,同時(shí)對(duì)COD和BOD5都有很好的降解效果.Oh等[29]通過(guò)向厭氧發(fā)酵罐接種乳酸桿屬?gòu)?qiáng)化初沉污泥水解酸化,SCOD增加59%以上,同時(shí)產(chǎn)生了醋酸、乳酸、丙酸、丁酸和戊酸等多種有機(jī)酸.醋酸桿菌屬也是一類化能異養(yǎng)菌,耐酸性,能將乙醇氧化成醋酸,并可將醋酸和乳酸氧化成CO2、水和H2等[30].兩種菌屬在污水處理過(guò)程中的有機(jī)物水解酸化過(guò)程中都扮演著重要角色,常用作復(fù)合微生物菌劑的主要成分[31-32].對(duì)中試裝置活性污泥高通量分析結(jié)果顯示,第3,10,21d的強(qiáng)化組(對(duì)照組)厭氧區(qū)變形門分別為67.95%(64.71%)、66.96%(77.44%)、64.56%(74.66%),厚壁門分別為0.91%(1.37%)、2.39%(1.42%)、3.04%(2.01%).微生物強(qiáng)化使強(qiáng)化組的厚壁門含量增加了2倍多,變形門的增長(zhǎng)受到抑制.屬水平上第3,10,21d的強(qiáng)化組(對(duì)照組)厭氧區(qū)醋酸桿屬相對(duì)豐度分別為3%(2%)、7%(2%)、8%(3%),乳桿菌屬相對(duì)豐度分別為5%(3%)、10%(4%)、9%(3%).可見(jiàn)菌劑中的功能組分在強(qiáng)化組中都顯著增長(zhǎng),并對(duì)活性污泥水解酸化起到了強(qiáng)化作用.

    群落解析結(jié)果越接近真實(shí),綜合考慮分類層次和鑒定率,對(duì)科水平的數(shù)據(jù)進(jìn)行深入解析.圖2反映了A/O中試裝置厭氧區(qū)和好氧區(qū)不同運(yùn)行階段的科水平的微生物群落組成,活性污泥中科水平微生物群落結(jié)構(gòu)變化明顯.強(qiáng)化組和對(duì)照組的厭氧區(qū)中高含量物種的總比例變化差明顯.所占比例分別是71.80%、72.10%和69.99%,呈現(xiàn)出穩(wěn)中微降的變化趨勢(shì);對(duì)照組高含量物種第3,10,21d所占比例分別是73.75%、75.03%和75.30%,呈現(xiàn)逐步升高的變化趨勢(shì).這說(shuō)明外源菌劑的投加改變了強(qiáng)化組厭氧區(qū)微生物組成和分布.具體的,Comamonadaceae是A/O裝置活性污泥中含量最多的微生物,強(qiáng)化組中的Comamonadaceae由第3d的15.49%降到第21d的12.23%,而對(duì)照組中Comamonadaceae基本不變;強(qiáng)化組中Sinobacteraceae由5.03%降低至3.01%,對(duì)照組Sinobacteraceae基本恒定;同樣還有Alcaligenaceae、Haliangiaceae、Polyangiaceae和Saprospiraceae等在強(qiáng)化組中都明顯減少,在對(duì)照組中基本不變.其中Alcaligenaceae是一類產(chǎn)堿的好氧菌,它在強(qiáng)化組中含量的減少,是由于外加菌劑中的乳酸桿菌和醋酸桿菌等酸性菌在強(qiáng)化污泥水解酸化時(shí)降低了體系的pH值,通過(guò)種間競(jìng)爭(zhēng)的方式抑制了Alcaligenaceae的增殖.Nitrosomonadaceae是亞硝化單胞菌科,在有機(jī)氮轉(zhuǎn)變成硝酸鹽的過(guò)程中起著重要作用.Nitrosomonadaceae變化規(guī)律說(shuō)明強(qiáng)化組對(duì)蛋白質(zhì)等有機(jī)氮分解能力明顯要強(qiáng)于對(duì)照組,這也是強(qiáng)化組的脫氮效果由于對(duì)照組的原因.

    綜上所述,菌劑中的微生物投加到活性污泥后通過(guò)兩方面作用強(qiáng)化污泥減量,一是菌劑中乳桿菌屬和醋酸菌屬等功能微生物自身直接與污泥發(fā)生作用,強(qiáng)化了污泥水解酸化,減少污泥產(chǎn)量,二是通過(guò)分泌胞外物質(zhì)(如抗生素)直接作用或通過(guò)改變生存環(huán)境(如pH值)等方式間接作用,促進(jìn)或抑制反應(yīng)體系中其它微生物的生長(zhǎng),改變?cè)钚晕勰嘀械奈⑸锶郝浣Y(jié)構(gòu),影響了污泥產(chǎn)量.

    2.4 微生物強(qiáng)化對(duì)群落多樣性的影響

    多樣性指數(shù)和均勻度指數(shù)是反映微生物群落變化的綜合量度,Siompson指數(shù)反映了群落內(nèi)微生物多樣性的概率,Shannon指數(shù)反映了群落多樣性的高低,Pielou指數(shù)反映了群落的均勻程度,三者均是指數(shù)越大,微生物分布越均勻,多樣性越高[24-26].圖3反映了科水平上的微生物群落多樣性變化,強(qiáng)化組厭氧區(qū)的Siompson指數(shù)隨著運(yùn)行天數(shù)增加而升高,由第3d的13.07升高到第21d的14.42,增加了10.33%,而對(duì)照組厭氧區(qū)的Siompson指數(shù)則是先降低后回升,到第21d變?yōu)?1.91,各個(gè)時(shí)期Siompson指數(shù)均是強(qiáng)化組高于對(duì)照組;好氧區(qū)Siompson指數(shù)則是強(qiáng)化組降低,對(duì)照組升高,強(qiáng)化組Siompson指數(shù)均低于對(duì)照組.強(qiáng)化組厭氧區(qū)Siompson指數(shù)升高是由于外加菌劑中多種外源微生物的加入增加了其微生物的種類,并改變了各微生物相對(duì)含量;而對(duì)照組Siompson指數(shù)降低則是因?yàn)轶w系中活性污泥取自羅芳污水廠好氧池,進(jìn)入中試裝置厭氧區(qū)后由于pH值、DO和營(yíng)養(yǎng)物質(zhì)等生活條件的差異導(dǎo)致微生物群落發(fā)生小幅度調(diào)整,多樣性降低.好氧區(qū)Siompson指數(shù)變化規(guī)律相反,這是因?yàn)榫鷦┲械墓δ芪⑸锒嗍菂捬趸蚣嫘詤捬跷⑸?在好氧區(qū)難以大量增殖,對(duì)其多樣性影響甚微. Shannon指數(shù)變化規(guī)律與Siompson指數(shù)相似.厭氧區(qū)的Pielou指數(shù)變化規(guī)律與Siompson指數(shù)相似,而兩組好氧區(qū)Pielou指數(shù)波動(dòng)較小,無(wú)明顯差異.本試驗(yàn)微生物菌劑主要投加到厭氧區(qū),其有效微生物經(jīng)過(guò)適當(dāng)時(shí)間的增殖并與土著微生物充分作用,抑制了部分土著微生物的增殖,使得強(qiáng)化組厭氧區(qū)多樣性升高,由此改變微生物群落結(jié)構(gòu),增加了水解酸化菌的含量,并強(qiáng)化反應(yīng)器中活性污泥的水解酸化,從而達(dá)到污泥減量的目的.

    2.5 微生物強(qiáng)化對(duì)群落演變的影響

    通過(guò)聚類分析研究了厭氧區(qū)和好氧區(qū)不同時(shí)期的微生物群落相似性(圖4),發(fā)現(xiàn)活性污泥微生物群落在時(shí)間尺度上發(fā)生明顯的演變,第3d和第5d的樣品聚為一類,第10d和第21d的樣品聚為一類,特別的,第 21d強(qiáng)化組厭氧區(qū)獨(dú)成一類.空間尺度上強(qiáng)化組和對(duì)照組對(duì)比來(lái)看,其厭氧區(qū)變化差異明顯:第3d的強(qiáng)化組和對(duì)照組距離最近,相似性最高;對(duì)照組第10d和第21d位于中間一支,與對(duì)照組第3d相比形成一定分化,d3-Y2與d21-Y2的差異為0.035;強(qiáng)化組厭氧區(qū)第3,10, 21d分別位于3個(gè)不同支,距離最遠(yuǎn),差異最大, d3-Y1與d21-Y1的差異為0.13;第21d時(shí)強(qiáng)化組和對(duì)照組群落差異達(dá)到0.09.這說(shuō)明相對(duì)于對(duì)照組厭氧區(qū),強(qiáng)化組厭氧區(qū)落隨時(shí)間的推移發(fā)生了更為明顯的演變,而外源微生物強(qiáng)化正是這種群落變化的驅(qū)動(dòng)力.

    對(duì)厭氧區(qū)和好氧區(qū)不同時(shí)期的微生物群落進(jìn)行了主坐標(biāo)分析,由圖5可知,強(qiáng)化組和對(duì)照組厭氧區(qū)的群落明顯分為2個(gè)區(qū)域,對(duì)照組群落沿著坐標(biāo)1方向演變,而強(qiáng)化組群落是沿著坐標(biāo)1和坐標(biāo)2對(duì)角線方向演變,并且強(qiáng)化組的演變梯度是對(duì)照組的近2倍.強(qiáng)化組和對(duì)照組污泥和進(jìn)水來(lái)源相同,中試運(yùn)行工況如DO、pH值、MLSS等均保持一致,有無(wú)菌劑投加是二者唯一的差異,因此可知微生物強(qiáng)化是強(qiáng)化組微生物群落沿著不同于對(duì)照組方向進(jìn)行演變的主要因素.對(duì)照組活性污泥由羅芳廠一期AA/O工藝好氧區(qū)引入到中試裝置后,該好氧污泥在對(duì)照組厭氧區(qū)由于DO等工藝參數(shù)的改變,其中微生物群落隨之發(fā)生演變,即DO的降低(由3~5mg/L降低為不足0.5mg/L)是微生物群落沿坐標(biāo)1方向的演變驅(qū)動(dòng)力.而微生物強(qiáng)化是微生物群落沿坐標(biāo)2方向的演變驅(qū)動(dòng)力,強(qiáng)化組微生物群落在2個(gè)驅(qū)動(dòng)力綜合作用沿坐標(biāo)1和坐標(biāo)2對(duì)角線方向演變,并且箭頭更靠近坐標(biāo)2,說(shuō)明微生物強(qiáng)化的驅(qū)動(dòng)作用更強(qiáng).兩組好氧區(qū)的微生物群落在坐標(biāo)1方向上演變較大,對(duì)照組在坐標(biāo)2方向上基本無(wú)變化,強(qiáng)化組在坐標(biāo)2方向上有小幅演變,這是因?yàn)槲⑸锞鷦┑挠行ЫM分是厭氧和兼性厭氧微生物,在好氧區(qū)無(wú)法大量繁殖,難以正常發(fā)揮作用.

    選取水溫、pH值、DO和SCOD、TP、TN、氨氮、硝酸鹽氮8個(gè)指標(biāo),通過(guò)典范對(duì)應(yīng)分析(Canonical Correspondence Analysis, CCA)研究科水平上微生物群落演變的關(guān)鍵影響因子.F1和F2的特征值分別為0.056和0.033,其解釋率分別為48.12%和28.43%,總解釋率為76.55%.pH值是微生物群落演變的最關(guān)鍵影響因子,其次是水溫和DO,而SCOD、TP、TN、氨氮、硝酸鹽氮的影響相對(duì)微弱.pH值的影響最大,說(shuō)明外加菌劑的微生物強(qiáng)化通過(guò)影響強(qiáng)化組活性污泥pH值這一生存環(huán)境來(lái)改變期微生物群落演替方向,使得更多的水解酸化細(xì)菌在較低pH值環(huán)境中充分繁殖,強(qiáng)化了對(duì)活性污泥的水解酸化作用,達(dá)到污泥減量的目的.由于試驗(yàn)是在12月末到1月期間進(jìn)行,水溫波動(dòng)幅度較大(8.9~17.9℃),對(duì)微生物群落演變產(chǎn)生較大影響.厭氧區(qū)DO為0.05~ 0.13mg/L,好氧區(qū)DO為3.09~3.75mg/L,DO的影響主要體現(xiàn)在厭氧區(qū)和好氧區(qū)微生物群落差異上.

    3 結(jié)論

    3.1 微生物強(qiáng)化使A/O中試裝置的MLSS和MLVSS分別減少了26.26%和28.08%,有效減少了剩余污泥產(chǎn)量.強(qiáng)化組的SS、COD、BOD5、TN和氨氮等污染物的去除率增加,但總磷的去除率略微降低,出水水質(zhì)可穩(wěn)定達(dá)到一級(jí)A標(biāo)準(zhǔn).

    3.2 高通量測(cè)序?qū)钚晕勰辔⑸锶郝湓陂T、綱、目、科水平鑒定結(jié)果良好,屬水平上也有一定鑒定效率,種水平的鑒定效果較差.微生物菌劑主要成分是厚壁門的乳桿菌屬和醋酸桿菌屬,微生物強(qiáng)化使強(qiáng)化組活性污泥中厚壁門增加了2倍,乳桿菌屬和醋酸菌屬在活性污泥中屬于低含量菌屬,二者在對(duì)照組中含量基本不變,在強(qiáng)化組含量明顯上升.

    3.3 微生物強(qiáng)化使強(qiáng)化組厭氧區(qū)Siompson指數(shù)、Shannon指數(shù)和Pielou指數(shù)均不同程度升高,群落多樣性增加.聚類分析顯示微生物群落按時(shí)間順序明顯聚為3簇,強(qiáng)化組厭氧區(qū)微生物群落隨時(shí)間發(fā)生較大演變.主坐標(biāo)分析顯示強(qiáng)化組和對(duì)照組群落明顯聚為2類,其厭氧區(qū)群落分別按照不同的方向演變,微生物強(qiáng)化和DO降低的協(xié)同作用是其演變的推動(dòng)力.典范對(duì)應(yīng)分析顯示顯著影響微生物群落結(jié)構(gòu)的環(huán)境因子依次為pH值、水溫和DO.

    3.4 微生物強(qiáng)化使得活性污泥微生物群落結(jié)構(gòu)發(fā)生明顯變化,菌劑通過(guò)直接與污泥作用強(qiáng)化污泥水解酸化,或通過(guò)改變生存環(huán)境(如pH值)等方式間接作用,促進(jìn)或抑制反應(yīng)體系中其它微生物的生長(zhǎng),改變?cè)钚晕勰嘀械奈⑸锶郝浣Y(jié)構(gòu),減少污泥產(chǎn)量.

    -chloroaniline slows bacterial species turnover rate in bioreactors [J]. Water Research, 2013,12,(47):7109-9119.

    陳 雄,韋祎旸,李輝信,等.附著劍菌對(duì)三氯聯(lián)苯的吸附及降解機(jī)制研究 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(5):1423-1428.

    Chang L K, Ibrahim D, Omar I C. A laboratory scale bioremediation of Tapis crude oil contaminated soil by bioaugmentation of Acinetobacter baumannii T30C [J]. African Journal Of Microbiology Research, 2011,9(5):2609-2615.

    侯 穎,王 飛,董維亮,等. Rhodococcus sp. T3-1菌株降解乙草胺的特性 [J]. 中國(guó)環(huán)境科學(xué), 2013,33(10):1785-1790.

    Zhu L L, Qiu L P, Zhang L X. Challenge and Development of Sludge Treatment and Disposal in China [J]. Applied Mechanics and Materials, 2014,8,(675-677):665-668.

    Wang S B. Sludge Management and Its Development in China [J]. Applied Mechanics and Materials, 2013,11,(665-668):1299- 1302.

    Vega E M, Hector G O. Optimizing chemical conditioning for odour removal of undigested sewage sludge in drying processes [J]. Journal Of Environmental Management, 2015,3(150):111- 119.

    Murugesan K, Selvam A. Flocculation and dewaterability of chemically enhanced primary treatment sludge by bioaugmentation with filamentous fungi [A]. International Conference on Solid Waste-Innovation in Technology and Management[C]. 2014,9(168):198-203.

    Wang Y, Qu P C H, Peng D C. Hydrolysis of excess sludge by bioaugmentation and nutrient removal from supernatant with low cost [C]//Proceedings of the 2011 International Conference on Remote Sensing, Environment and Transportation Engineering, 2011,6:4096-9.

    李 俊.微生物菌劑對(duì)污水處理廠污泥減量的影響研究 [D]. 重慶:西南大學(xué), 2008.

    張海珍.MCMP-II多功能復(fù)合微生物制劑用于城市污水處理廠污泥減量的生產(chǎn)性試驗(yàn)研究 [D]. 重慶:重慶大學(xué), 2008.

    黃黎明.微生物制劑IBC強(qiáng)化SBR處理系統(tǒng)的研究 [D]. 廣州:華南理工大學(xué), 2012.

    Muyzer G, Dewaal EC, Uitterlinden AG. Profiling of complex microbial-populations by denaturing gradient gel-electrophoresis analysis of polymerase chain reaction-amplified genes-coding for 16S ribosomal-rna [J]. Applied And Environmental Microbiology, 1993,59(3):695-700.

    趙 鑫,劉一威,李曉東,等.微生物復(fù)合菌劑法污泥原位減量群落動(dòng)態(tài)研究 [J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2015,36(3):839-845.

    薛圓圓,孫寶盛,杜 江,等.貧營(yíng)養(yǎng)條件下IAMBR污泥微生物群落結(jié)構(gòu)的演變 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(3):839-845.

    陳 燕,劉國(guó)華,范 強(qiáng),等.不同溶解氧條件下A/O系統(tǒng)的除碳脫氮效果和細(xì)菌群落結(jié)構(gòu)變化環(huán)境科學(xué), 2015,36(7): 2610-2616.

    Reuter J A, Spacek D V, Snyder M P. High-Throughput Sequencing Technologies [J]. Molecular Cell, 2015,3(58):586- 597.

    Picard C, Fischer A l. Contribution of high-throughput DNA sequencing to the study of primary immunodeficiencies [J]. European Journal of Immunology, 2014,10(44):2854-2861.

    Georgiou G, Ippolito G C. Beausang J. The promise and challenge of high-throughput sequencing of the antibody repertoire [J]. Nature Biotechnology, 2014,2,(32):158-168.

    McAdam P R, Richardson E J, Fitzgerald J. High-throughput sequencing for the study of bacterial pathogen biology [J]. Current Opinion in Microbiology, 2014,6,(19):106-113.

    佳,宋云龍,張金松,等.基于微生物強(qiáng)化的原位污泥減量中試研究 [J]. 給水排水, 2014,40(8):37-43.

    GB18918-2002 城鎮(zhèn)污水處理廠污染物排放標(biāo)準(zhǔn) [S].

    Turnbaugh P J, Hamady M, Yatsunenko T, et al. A core gut microbiome in obese and lean twins [J]. Nature, 2009,457(7228): 480–487.

    Simpson E H. Measurement of diversity [J]. Nature, 1949, 163:688.

    Hill M O. Diversity and evenness: a unifying notation and its consequences [J]. Ecology, 1973,54:427-431.

    Pielou E C. The measurement of diversity in different types of biological collections [J]. Journal of Theoretical Biology, 1966, 13:131-144.

    Zhang, Y X, Vadlani P V. Lactic acid production from biomass-derived sugars via co-fermentation of Lactobacillus brevis and Lactobacillus plantarum [J]. Journal of Bioscience and Bioengineering, 2015,6(119):694-699.

    Vellingiri K, Ramachandran T. Isolation and Purification of Lactobacillus Acidophilus and Analyzing its Influence on Effluent Treatment [J]. International Journal of Engineering and Technology Innovation, 2015,1(5):66-74.

    Oh M H, Lee, S M, Hong S H. Monitoring of Lactobacillus sp inoculated in the reactor to evaluate the solubilization efficiency of primary sludge [J]. International Biodeterioration & Biodegradation, 2013,11(85):603-607.

    Demler M, Weuster-Botz D. Reaction Engineering Analysis of Hydrogenotrophic Production of Acetic Acid by Acetobacterium woodii [J]. Biotechnology and Bioengineering, 2011,2(108):470- 474.

    YIN S, LI S. Ethanol-tolerant acetogenic microorganism composite bacteria used for anaerobic fermentation of medicine residue alcohol comprises Acetobacterium carbinolicum, Clostridium irregularic, Clostridium cellobioparum, and Clostridium glycolicum: China, 103555638-A [P].

    Fan H, Qu H, Zhang L. Complex microbial agent useful for treating wastewater, contains Rhodopseudomonas palustris, Acetobacterium, Bacillus subtilis, Metallosphaera sedula, Metallosphaera prunae, Nitrosomonas, Lactococcus, and Flavobacterium: China, 104195090-A [P].

    * 責(zé)任作者, 教授, zhangjinsong@waterchina.com

    Analysis of microbial community in in-situ sludge reduction process by bioaugmentation using high-throughput sequencing technology

    SONG Yun-long1, ZHANG Jin-song1*, ZHU Jia2, SHAO Ming-fei1, DENG Ren-jian1

    (1.School of Civil and Environment Engineering, Harbin Institute of Technology Shenzhen Graduate School, Shenzhen 518055, China;2.School of Civil and Environmental Engineering, Shenzhen Polytechnic, Shenzhen 518055, China). China Environmental Science, 2016,36(7):2099~2107

    High-throughput sequencing was used to parse the microbial community of activated sludge in an in-situ sludge reduction process. Microbial community structure shifted dramatically as the result of bioaugmentation. The abundance of and increased significantly in bioaugmentation group. Community diversity indices including Simpson index, Shannon index and Pielou index rised at different degrees in anaerobic zone of bioaugmentation group. Cluster analysis indicated that microbial communities were obviously grouped into three clusters in chronological order. In anaerobic zone of bioaugmentation group, the community structure shifted more obviously over time. Through principal coordinate analysis, communities of bioaugmentation group and control group were grouped into two distinct areas. Synergistic efforts of bioaugmentation and DO reduction were the driving forces of community shift in bioaugmentation group, while DO reduction was the driving force in control group. Canonical correspondence analysis showed that pH, water temperature and DO were the dominant environmental factors, which affected the microbial community structure significantly.

    high-throughput sequencing;bioaugmentation;hydrolytic acidification;microbial community analysis

    X172

    A

    1000-6923(2016)07-2099-09

    宋云龍(1983-),男,山東威海人,哈爾濱工業(yè)大學(xué)博士研究生,主要從事污泥減量及其資源化利用的研究.

    2015-11-06

    國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目(2009AA064704)

    猜你喜歡
    氧區(qū)桿菌屬活性污泥
    自動(dòng)清洗高效生物膜反應(yīng)器設(shè)計(jì)與性能
    化工管理(2024年1期)2024-01-17 08:49:20
    潰瘍性結(jié)腸炎患者腸道菌群分布特征分析
    一種污水處理系統(tǒng)及處理方法與應(yīng)用
    能源化工(2021年6期)2021-12-30 15:41:26
    養(yǎng)豬微生物發(fā)酵床芽胞桿菌空間生態(tài)位特性
    活性污泥系統(tǒng)ASM2d模型的自適應(yīng)模糊PID控制
    活性污泥系統(tǒng)ASM2d模型的自適應(yīng)模糊PID控制
    污水活性污泥處理過(guò)程的溶解氧增益調(diào)度控制
    多段多級(jí)AO除磷脫氮工藝
    活性污泥對(duì)管式多孔α-Al2O3陶瓷膜支撐體性能影響的初步研究
    類芽孢桿菌屬β-葡萄糖苷酶在大腸桿菌中可溶性重組表達(dá)的優(yōu)化
    日韩欧美 国产精品| 最近在线观看免费完整版| 亚洲国产精品999在线| 国产色爽女视频免费观看| 最新在线观看一区二区三区| 观看美女的网站| 香蕉丝袜av| 91在线精品国自产拍蜜月 | 亚洲人成电影免费在线| 国产亚洲精品久久久com| 成人三级黄色视频| 日韩成人在线观看一区二区三区| 久久婷婷人人爽人人干人人爱| 深爱激情五月婷婷| www.www免费av| 非洲黑人性xxxx精品又粗又长| 国产v大片淫在线免费观看| 女人被狂操c到高潮| 亚洲精品乱码久久久v下载方式 | 欧美av亚洲av综合av国产av| 久久婷婷人人爽人人干人人爱| 亚洲国产欧美人成| 亚洲五月天丁香| 国产欧美日韩一区二区精品| 无人区码免费观看不卡| 国产久久久一区二区三区| 成熟少妇高潮喷水视频| 尤物成人国产欧美一区二区三区| 成熟少妇高潮喷水视频| 国产综合懂色| 中文字幕高清在线视频| h日本视频在线播放| 身体一侧抽搐| 999久久久精品免费观看国产| svipshipincom国产片| 免费在线观看亚洲国产| 桃色一区二区三区在线观看| 熟女人妻精品中文字幕| 狂野欧美白嫩少妇大欣赏| 啪啪无遮挡十八禁网站| 此物有八面人人有两片| 老熟妇仑乱视频hdxx| 波多野结衣高清无吗| 12—13女人毛片做爰片一| 亚洲国产精品sss在线观看| 男女下面进入的视频免费午夜| 搡女人真爽免费视频火全软件 | 99国产精品一区二区三区| 熟女电影av网| 精品电影一区二区在线| 嫩草影视91久久| 欧美性猛交╳xxx乱大交人| 欧美日韩乱码在线| 国产亚洲精品久久久com| 亚洲一区高清亚洲精品| 国产探花在线观看一区二区| 国产又黄又爽又无遮挡在线| e午夜精品久久久久久久| 丰满人妻熟妇乱又伦精品不卡| 好看av亚洲va欧美ⅴa在| 老司机在亚洲福利影院| 不卡一级毛片| 嫩草影视91久久| 久久精品国产自在天天线| 无遮挡黄片免费观看| 啦啦啦免费观看视频1| 成人欧美大片| 国产麻豆成人av免费视频| 精品熟女少妇八av免费久了| 欧美3d第一页| 国产色婷婷99| 又黄又爽又免费观看的视频| 日韩欧美 国产精品| 国产av在哪里看| 男插女下体视频免费在线播放| 久久精品亚洲精品国产色婷小说| 国内久久婷婷六月综合欲色啪| 欧美成人免费av一区二区三区| 法律面前人人平等表现在哪些方面| 亚洲天堂国产精品一区在线| 91九色精品人成在线观看| 99视频精品全部免费 在线| 久久欧美精品欧美久久欧美| 日韩欧美精品v在线| 婷婷精品国产亚洲av在线| 亚洲男人的天堂狠狠| 十八禁人妻一区二区| 午夜福利视频1000在线观看| 国产单亲对白刺激| 又爽又黄无遮挡网站| 亚洲男人的天堂狠狠| 国产aⅴ精品一区二区三区波| 成人精品一区二区免费| 亚洲中文日韩欧美视频| 一级毛片高清免费大全| 国产视频一区二区在线看| 国产乱人视频| av欧美777| 18禁在线播放成人免费| 国产美女午夜福利| www.熟女人妻精品国产| 日韩国内少妇激情av| 又粗又爽又猛毛片免费看| 久久亚洲精品不卡| 色吧在线观看| 九九久久精品国产亚洲av麻豆| 欧洲精品卡2卡3卡4卡5卡区| xxxwww97欧美| 欧美成狂野欧美在线观看| 男女下面进入的视频免费午夜| 亚洲片人在线观看| 亚洲专区国产一区二区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产高清在线一区二区三| 成人国产一区最新在线观看| 久久这里只有精品中国| 国产精品爽爽va在线观看网站| a在线观看视频网站| 久久久国产成人精品二区| 免费人成在线观看视频色| 亚洲精品一卡2卡三卡4卡5卡| 日本与韩国留学比较| 久久精品影院6| 久久精品91无色码中文字幕| 女警被强在线播放| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 成年免费大片在线观看| 一级毛片女人18水好多| 国产老妇女一区| 亚洲成人中文字幕在线播放| 国产一区二区在线av高清观看| 老鸭窝网址在线观看| 18禁黄网站禁片免费观看直播| 国产亚洲精品一区二区www| 内射极品少妇av片p| 白带黄色成豆腐渣| 1000部很黄的大片| 国产一区在线观看成人免费| 国内揄拍国产精品人妻在线| 十八禁网站免费在线| 狠狠狠狠99中文字幕| 女警被强在线播放| 亚洲av成人精品一区久久| 亚洲精品一卡2卡三卡4卡5卡| 国产野战对白在线观看| 哪里可以看免费的av片| 99riav亚洲国产免费| 蜜桃久久精品国产亚洲av| 熟妇人妻久久中文字幕3abv| 日本免费一区二区三区高清不卡| 中文字幕人妻丝袜一区二区| 国产精华一区二区三区| 久久国产精品人妻蜜桃| 黄色成人免费大全| 男女做爰动态图高潮gif福利片| 麻豆国产97在线/欧美| 老汉色∧v一级毛片| 精品久久久久久久人妻蜜臀av| 精品一区二区三区人妻视频| 一级a爱片免费观看的视频| 国产精品精品国产色婷婷| 狂野欧美白嫩少妇大欣赏| 国产av麻豆久久久久久久| 在线天堂最新版资源| 草草在线视频免费看| 女生性感内裤真人,穿戴方法视频| 天天躁日日操中文字幕| 欧美乱妇无乱码| 黄色丝袜av网址大全| www日本在线高清视频| 免费看a级黄色片| 少妇的逼好多水| 男人和女人高潮做爰伦理| 日韩欧美在线二视频| 又黄又粗又硬又大视频| 亚洲真实伦在线观看| 日本a在线网址| 又粗又爽又猛毛片免费看| 床上黄色一级片| 亚洲无线在线观看| 国产97色在线日韩免费| 手机成人av网站| 国产男靠女视频免费网站| 一级黄色大片毛片| 精品午夜福利视频在线观看一区| 国产单亲对白刺激| 脱女人内裤的视频| 又紧又爽又黄一区二区| 久久精品国产自在天天线| 搡老妇女老女人老熟妇| 日本撒尿小便嘘嘘汇集6| av专区在线播放| 91久久精品电影网| 色吧在线观看| 欧美丝袜亚洲另类 | 高清在线国产一区| 日本与韩国留学比较| 99riav亚洲国产免费| netflix在线观看网站| 国内揄拍国产精品人妻在线| 日本撒尿小便嘘嘘汇集6| 一级黄色大片毛片| а√天堂www在线а√下载| 国产色婷婷99| 国产欧美日韩一区二区三| 黄色片一级片一级黄色片| 欧美不卡视频在线免费观看| 精品久久久久久久毛片微露脸| 国产精品 欧美亚洲| 国产高清视频在线播放一区| 中国美女看黄片| 一区二区三区国产精品乱码| 亚洲av免费在线观看| 亚洲一区高清亚洲精品| 亚洲 国产 在线| 夜夜爽天天搞| 好看av亚洲va欧美ⅴa在| 国产精品,欧美在线| 老司机午夜十八禁免费视频| 久久香蕉精品热| 在线观看美女被高潮喷水网站 | 精品久久久久久久久久久久久| 变态另类成人亚洲欧美熟女| 精品福利观看| 麻豆一二三区av精品| 亚洲最大成人中文| 午夜福利18| 日本免费a在线| 久久亚洲真实| 国产伦精品一区二区三区视频9 | 变态另类丝袜制服| 最后的刺客免费高清国语| 51国产日韩欧美| 最近视频中文字幕2019在线8| 亚洲精品影视一区二区三区av| 在线观看日韩欧美| 亚洲在线观看片| 网址你懂的国产日韩在线| 3wmmmm亚洲av在线观看| 叶爱在线成人免费视频播放| 热99re8久久精品国产| 日日摸夜夜添夜夜添小说| 97碰自拍视频| 久久久精品欧美日韩精品| 国产97色在线日韩免费| 免费看光身美女| 18禁黄网站禁片免费观看直播| 久久久久久久精品吃奶| 精品人妻一区二区三区麻豆 | netflix在线观看网站| 美女大奶头视频| 日本 欧美在线| 久久人妻av系列| 每晚都被弄得嗷嗷叫到高潮| 俄罗斯特黄特色一大片| 欧美中文综合在线视频| 中文亚洲av片在线观看爽| 高清日韩中文字幕在线| 老司机福利观看| 国产爱豆传媒在线观看| 真人一进一出gif抽搐免费| 色老头精品视频在线观看| 亚洲无线在线观看| 国产精品野战在线观看| 国产成年人精品一区二区| e午夜精品久久久久久久| 又黄又粗又硬又大视频| 国产精品永久免费网站| 久久久久免费精品人妻一区二区| 国产精品影院久久| 国产一区二区亚洲精品在线观看| 黄片大片在线免费观看| 脱女人内裤的视频| 嫩草影视91久久| 18禁美女被吸乳视频| 国产精品永久免费网站| 麻豆久久精品国产亚洲av| 床上黄色一级片| 99久久综合精品五月天人人| 午夜福利在线观看免费完整高清在 | av视频在线观看入口| 特大巨黑吊av在线直播| 色老头精品视频在线观看| 国产精品一及| 一进一出抽搐动态| 国产伦一二天堂av在线观看| 亚洲美女视频黄频| 亚洲欧美日韩高清专用| 久久精品国产综合久久久| 搡老岳熟女国产| 国产成人影院久久av| 国产激情欧美一区二区| 亚洲国产精品999在线| www.www免费av| 看免费av毛片| 久久精品国产清高在天天线| 给我免费播放毛片高清在线观看| 亚洲av中文字字幕乱码综合| 搡老岳熟女国产| 亚洲最大成人手机在线| 窝窝影院91人妻| 免费av毛片视频| 校园春色视频在线观看| 亚洲精品成人久久久久久| 三级毛片av免费| 亚洲专区国产一区二区| 又爽又黄无遮挡网站| 狂野欧美白嫩少妇大欣赏| 男人的好看免费观看在线视频| 最新美女视频免费是黄的| 又紧又爽又黄一区二区| 五月伊人婷婷丁香| 久久精品91无色码中文字幕| 天堂影院成人在线观看| 好男人在线观看高清免费视频| 亚洲人与动物交配视频| a在线观看视频网站| 变态另类丝袜制服| 亚洲av日韩精品久久久久久密| 国产淫片久久久久久久久 | 国产精品久久久久久人妻精品电影| 十八禁人妻一区二区| 美女黄网站色视频| 搡老妇女老女人老熟妇| 亚洲国产精品999在线| 亚洲av免费在线观看| 国产亚洲av嫩草精品影院| 亚洲avbb在线观看| 一区二区三区高清视频在线| 综合色av麻豆| 一级黄片播放器| 欧美日韩瑟瑟在线播放| 久久欧美精品欧美久久欧美| 俄罗斯特黄特色一大片| 精品熟女少妇八av免费久了| 日日夜夜操网爽| 又紧又爽又黄一区二区| 欧美黄色片欧美黄色片| 精品久久久久久久末码| 国产精品久久电影中文字幕| 欧美绝顶高潮抽搐喷水| 99精品欧美一区二区三区四区| 禁无遮挡网站| 国产毛片a区久久久久| www.999成人在线观看| 99久久综合精品五月天人人| 国内精品美女久久久久久| 亚洲无线在线观看| 尤物成人国产欧美一区二区三区| 日韩欧美在线乱码| 亚洲精品成人久久久久久| 操出白浆在线播放| 亚洲国产日韩欧美精品在线观看 | 俺也久久电影网| e午夜精品久久久久久久| 日韩中文字幕欧美一区二区| 99精品欧美一区二区三区四区| 久久久久性生活片| 小说图片视频综合网站| 他把我摸到了高潮在线观看| 深爱激情五月婷婷| 蜜桃亚洲精品一区二区三区| 白带黄色成豆腐渣| 午夜福利免费观看在线| 成熟少妇高潮喷水视频| 国产精品综合久久久久久久免费| 天天添夜夜摸| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 久久99热这里只有精品18| 在线观看日韩欧美| 熟妇人妻久久中文字幕3abv| 18禁黄网站禁片免费观看直播| 久久久久国内视频| 国产午夜精品论理片| 久久久久久九九精品二区国产| 我要搜黄色片| 久久久久久久亚洲中文字幕 | 日韩成人在线观看一区二区三区| 99riav亚洲国产免费| 一个人看视频在线观看www免费 | 免费观看精品视频网站| 国产极品精品免费视频能看的| 久久久久久大精品| 最近最新中文字幕大全电影3| 亚洲精品在线美女| 一二三四社区在线视频社区8| 午夜老司机福利剧场| 内射极品少妇av片p| 一a级毛片在线观看| 亚洲无线在线观看| 国产日本99.免费观看| 久久精品国产亚洲av涩爱 | 男女做爰动态图高潮gif福利片| 成人无遮挡网站| 欧美乱码精品一区二区三区| 国产精品久久久人人做人人爽| 美女 人体艺术 gogo| 黄色成人免费大全| 波多野结衣高清作品| 无人区码免费观看不卡| 欧美日韩精品网址| 亚洲熟妇熟女久久| 91麻豆av在线| 日韩精品青青久久久久久| 18美女黄网站色大片免费观看| 欧美日韩国产亚洲二区| 亚洲激情在线av| 成熟少妇高潮喷水视频| 欧美区成人在线视频| 女人十人毛片免费观看3o分钟| 夜夜躁狠狠躁天天躁| 欧美av亚洲av综合av国产av| 国产亚洲精品久久久com| 国产精品99久久久久久久久| 国产不卡一卡二| 欧美乱妇无乱码| 国产精品乱码一区二三区的特点| 精品国产三级普通话版| 日韩有码中文字幕| 国产毛片a区久久久久| eeuss影院久久| 免费在线观看亚洲国产| 在线a可以看的网站| 中文字幕人妻丝袜一区二区| 在线免费观看不下载黄p国产 | 国产精品 欧美亚洲| 亚洲精品日韩av片在线观看 | 99久久精品一区二区三区| 色在线成人网| 久久天躁狠狠躁夜夜2o2o| 成人一区二区视频在线观看| 91av网一区二区| 在线观看午夜福利视频| 久久久久久久亚洲中文字幕 | 欧美三级亚洲精品| 免费av毛片视频| 亚洲av免费高清在线观看| a级一级毛片免费在线观看| 免费av毛片视频| 亚洲第一电影网av| 女人高潮潮喷娇喘18禁视频| 很黄的视频免费| 日韩人妻高清精品专区| 成人无遮挡网站| 国产精品亚洲美女久久久| 一a级毛片在线观看| 国产成人aa在线观看| 国产午夜精品论理片| 一级毛片高清免费大全| 神马国产精品三级电影在线观看| 欧美日韩中文字幕国产精品一区二区三区| 免费无遮挡裸体视频| 欧美zozozo另类| 国内精品一区二区在线观看| 搡老妇女老女人老熟妇| 欧美日韩综合久久久久久 | 午夜福利18| 亚洲久久久久久中文字幕| 国产伦精品一区二区三区视频9 | 男女做爰动态图高潮gif福利片| 中文字幕人妻熟人妻熟丝袜美 | 高清在线国产一区| 美女大奶头视频| 欧美色欧美亚洲另类二区| 悠悠久久av| 最近视频中文字幕2019在线8| 美女被艹到高潮喷水动态| 色噜噜av男人的天堂激情| 久久久国产成人精品二区| 午夜免费成人在线视频| 一级作爱视频免费观看| 99视频精品全部免费 在线| 嫩草影院入口| 少妇熟女aⅴ在线视频| 日本a在线网址| 成人永久免费在线观看视频| 国产成人福利小说| 真实男女啪啪啪动态图| 精华霜和精华液先用哪个| 两人在一起打扑克的视频| 亚洲 国产 在线| 欧美日韩福利视频一区二区| 99久久久亚洲精品蜜臀av| 亚洲成av人片免费观看| 在线观看免费视频日本深夜| 欧美zozozo另类| 日本熟妇午夜| 深爱激情五月婷婷| 黑人欧美特级aaaaaa片| 久久久成人免费电影| 人妻久久中文字幕网| 97人妻精品一区二区三区麻豆| 97超级碰碰碰精品色视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 可以在线观看的亚洲视频| 18禁黄网站禁片免费观看直播| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 十八禁网站免费在线| 久久久久免费精品人妻一区二区| 免费看十八禁软件| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美精品综合久久99| 97超视频在线观看视频| 免费av观看视频| 精品久久久久久久毛片微露脸| 波多野结衣巨乳人妻| 欧美性猛交黑人性爽| 成人一区二区视频在线观看| 男女床上黄色一级片免费看| 国产又黄又爽又无遮挡在线| 深夜精品福利| 精品国内亚洲2022精品成人| 国产av不卡久久| 国产欧美日韩精品一区二区| 噜噜噜噜噜久久久久久91| 欧美+日韩+精品| 日日干狠狠操夜夜爽| 亚洲av二区三区四区| 99国产精品一区二区蜜桃av| 少妇高潮的动态图| 久久久久久久午夜电影| 亚洲国产中文字幕在线视频| 亚洲欧美日韩无卡精品| xxxwww97欧美| 亚洲乱码一区二区免费版| 免费观看的影片在线观看| www国产在线视频色| 久久精品人妻少妇| 精品熟女少妇八av免费久了| 成人特级黄色片久久久久久久| 国产亚洲精品久久久com| 精品一区二区三区人妻视频| 午夜精品在线福利| 不卡一级毛片| 1024手机看黄色片| 精品电影一区二区在线| 非洲黑人性xxxx精品又粗又长| 日本五十路高清| 日本黄色片子视频| 精品久久久久久,| 色视频www国产| 午夜影院日韩av| 国产精品久久久久久久久免 | 9191精品国产免费久久| 综合色av麻豆| 18+在线观看网站| 国产黄色小视频在线观看| 免费看光身美女| 男女视频在线观看网站免费| 亚洲专区国产一区二区| 久久久久九九精品影院| 亚洲精品在线美女| 99在线视频只有这里精品首页| 不卡一级毛片| 亚洲人成电影免费在线| 男人舔奶头视频| xxxwww97欧美| 精品久久久久久,| 亚洲男人的天堂狠狠| 色综合亚洲欧美另类图片| 精品久久久久久久久久免费视频| АⅤ资源中文在线天堂| 亚洲中文日韩欧美视频| 少妇裸体淫交视频免费看高清| 国产精品乱码一区二三区的特点| 成年人黄色毛片网站| 日韩 欧美 亚洲 中文字幕| 国产欧美日韩精品一区二区| 嫩草影院入口| 2021天堂中文幕一二区在线观| 亚洲人成网站高清观看| 天美传媒精品一区二区| 国产老妇女一区| 18禁黄网站禁片午夜丰满| 色综合亚洲欧美另类图片| 99久久99久久久精品蜜桃| 禁无遮挡网站| 午夜福利欧美成人| 人人妻人人看人人澡| 嫁个100分男人电影在线观看| 成人三级黄色视频| 免费观看精品视频网站| 亚洲精品影视一区二区三区av| 男人舔奶头视频| 日韩欧美国产在线观看| 19禁男女啪啪无遮挡网站| 在线播放无遮挡| 国产一区在线观看成人免费| 欧美在线一区亚洲| 国产精品乱码一区二三区的特点| 日韩欧美免费精品| 国产探花极品一区二区| 欧美zozozo另类| 亚洲中文字幕日韩| 一级作爱视频免费观看| 国产野战对白在线观看| 欧美中文综合在线视频| 男女床上黄色一级片免费看| 色综合婷婷激情| 91久久精品国产一区二区成人 | 午夜影院日韩av| 精品欧美国产一区二区三| 日日摸夜夜添夜夜添小说| 亚洲狠狠婷婷综合久久图片| 午夜日韩欧美国产| 一个人免费在线观看电影| 少妇高潮的动态图| 桃红色精品国产亚洲av| 91久久精品电影网| 日本撒尿小便嘘嘘汇集6| 中文字幕人成人乱码亚洲影| 天堂影院成人在线观看| 亚洲国产欧美人成| www日本在线高清视频|