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

    三價(jià)鐵對(duì)有機(jī)物存在下厭氧氨氧化脫氮的影響

    2021-04-29 14:28:48王海月毛念佳耿金菊任洪強(qiáng)
    中國(guó)環(huán)境科學(xué) 2021年4期
    關(guān)鍵詞:三價(jià)亞硝酸鹽硝化

    王海月,彭 玲,毛念佳,耿金菊,任洪強(qiáng),許 柯

    三價(jià)鐵對(duì)有機(jī)物存在下厭氧氨氧化脫氮的影響

    王海月,彭 玲,毛念佳,耿金菊,任洪強(qiáng),許 柯*

    (南京大學(xué)環(huán)境學(xué)院,污染控制與資源再利用國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210023)

    考察了三價(jià)鐵(2.24~7.84mg/L)存在下厭氧氨氧化系統(tǒng)對(duì)有機(jī)物的耐受性能,并通過16SrRNA高通量測(cè)序技術(shù)和定量PCR探究其機(jī)理.結(jié)果表明,進(jìn)水COD濃度為50和100mg/L時(shí),4個(gè)反應(yīng)器的氨氮和總氮去除率均較高(>90%),三價(jià)鐵的強(qiáng)化作用不明顯;進(jìn)水COD濃度繼續(xù)升高(150和200mg/L),厭氧氨氧化受到抑制,三價(jià)鐵的強(qiáng)化作用逐漸增加;COD濃度為200mg/L時(shí),添加三價(jià)鐵(7.84mg/L)可將氨氮和總氮去除率由61.3%和79.8%(對(duì)照組)提升至71.2%和84.7%.16SrRNA高通量測(cè)序技術(shù)表明,有機(jī)物存在下,污泥微生物群落結(jié)構(gòu)出現(xiàn)變化,主要表現(xiàn)為厭氧氨氧化菌豐度的降低及反硝化菌群的大量增殖,進(jìn)水添加三價(jià)鐵可提高浮霉菌(Planctomycetes)的豐度.定量PCR結(jié)果表明, 三價(jià)鐵能夠提高厭氧氨氧化菌16S rRNA及功能基因的豐度.

    三價(jià)鐵;生物脫氮;厭氧氨氧化;厭氧顆粒污泥膨脹床;COD

    厭氧氨氧化(Anammox)可在厭氧條件下實(shí)現(xiàn)氨氮和亞硝酸鹽的同時(shí)去除,相較于傳統(tǒng)硝化反硝化工藝,具有節(jié)省60%曝氣消耗,無需外加有機(jī)碳源,減少90%污泥產(chǎn)生量等優(yōu)勢(shì),近年來得到廣泛關(guān)注[1-2].

    迄今為止,全世界已建立了200多個(gè)基于Anammox工藝的全規(guī)模污水處理廠[3].Anammox工藝主要適用于處理低碳氮比的高氨氮廢水[4],而實(shí)際廢水中含有濃度,種類不同的有機(jī)物質(zhì)[5],通常認(rèn)為有機(jī)物的存在會(huì)對(duì)厭氧氨氧化菌產(chǎn)生負(fù)面影響[6-8].在Anammox系統(tǒng)中,脫氮性能取決于厭氧氨氧化和反硝化的競(jìng)爭(zhēng)共存,但是,由于反硝化反應(yīng)在熱力學(xué)上比厭氧氨氧化反應(yīng)更有利,以及厭氧氨氧化菌的低生長(zhǎng)速率和低產(chǎn)率,在有足夠化學(xué)需氧量 (COD)的系統(tǒng)中很難維持厭氧氨氧化菌的高豐度,從而影響系統(tǒng)的穩(wěn)定性[9].相關(guān)研究發(fā)現(xiàn)進(jìn)水COD濃度約為280mg/L時(shí)會(huì)完全抑制厭氧氨氧化過程[6].因此,為避免有機(jī)負(fù)荷過大,通常在主流厭氧氨氧化工藝中采取有機(jī)碳分離預(yù)處理措施.例如,以混合厭氧反應(yīng)器(HAR)作為預(yù)處理捕集COD并進(jìn)行資源回收[10],但低溫和低底物濃度會(huì)阻礙厭氧預(yù)處理的應(yīng)用.由此可見,當(dāng)前的研究重點(diǎn)是提高Anammox工藝自身耐受有機(jī)物的性能,開發(fā)可行有效的新技術(shù)至關(guān)重要.

    鐵作為含鐵蛋白的組成部分,在微生物的代謝和生長(zhǎng)中起著關(guān)鍵作用[11].鐵元素可能影響厭氧氨氧化菌的代謝,Van Niftrik等[12]發(fā)現(xiàn)厭氧氨氧化菌體內(nèi)含有鐵儲(chǔ)存顆粒.研究表明,在合理的Fe2+和Fe3+濃度下,Anammox的活性可以通過維持較高的氧化還原電位和聚集?;?高絲氨酸內(nèi)酯(AHLs)來增強(qiáng),而AHLs被認(rèn)為是細(xì)胞間通訊的重要物質(zhì)[13].Wang等[14]發(fā)現(xiàn)添加到聚乙烯醇/殼聚糖/鐵凝膠珠中的鐵離子與帶負(fù)電的細(xì)胞外聚合物之間的相互作用可以提升Anammox生物質(zhì)顆粒的致密性,從而有效地提高反應(yīng)器中的生物質(zhì)保留能力,保護(hù)細(xì)胞抵御外部壓力.Yin等[15]報(bào)道添加Fe3+(5mg/L) 促進(jìn)了海洋厭氧氨氧化菌的生長(zhǎng)速度,從而實(shí)現(xiàn)了更高的總氮去除率.此外,研究發(fā)現(xiàn)Fe3+濃度從2.24mg/L增加至5.6mg/L時(shí),厭氧氨氧化菌生長(zhǎng)速率從0.1343d-1增加到0.1709d-1[16].然而,關(guān)于Fe3+對(duì)有機(jī)物存在下厭氧氨氧化系統(tǒng)的脫氮規(guī)律變化及污泥微生物群落動(dòng)態(tài)變化的研究報(bào)道尚少,且Fe3+對(duì)相關(guān)氮素轉(zhuǎn)化功能基因的長(zhǎng)期影響仍不清楚.

    本文系統(tǒng)考察了三價(jià)鐵存在下厭氧氨氧化系統(tǒng)對(duì)有機(jī)物的耐受性能,闡述了三價(jià)鐵對(duì)厭氧氨氧化工藝運(yùn)行中氮素轉(zhuǎn)化,污泥活性的影響,分析了微生物群落結(jié)構(gòu)的變化及相關(guān)氮素轉(zhuǎn)化功能基因的表達(dá).通過這些研究,為三價(jià)鐵強(qiáng)化厭氧氨氧化工藝提供理論基礎(chǔ).

    1 材料與方法

    1.1 實(shí)驗(yàn)裝置

    實(shí)驗(yàn)采用厭氧膨脹床反應(yīng)器(EGSB),EGSB反應(yīng)器能使微生物有效持留,并使污泥流態(tài)化,促進(jìn)基質(zhì)傳遞及產(chǎn)氣釋放.此外,維持污泥膨脹所進(jìn)行的出水回流,可削弱基質(zhì)自抑制效應(yīng),非常適于厭氧氨氧化工藝的啟動(dòng)及運(yùn)行.實(shí)驗(yàn)采用的ESGB反應(yīng)器由有機(jī)玻璃制成,其高徑比為6,有效容積為0.8L.模擬廢水由蠕動(dòng)泵泵入反應(yīng)器,在升流運(yùn)動(dòng)中完成氮素轉(zhuǎn)化,氣體由頂部的氣室排出,出水由溢流堰排放.反應(yīng)器出水回流比為8,液面上升速度為0.1m/h,為防止光的負(fù)面影響,反應(yīng)器外裹錫紙.實(shí)驗(yàn)裝置如圖1所示.

    實(shí)驗(yàn)采用4個(gè)完全一致的EGSB反應(yīng)器,編號(hào)R1,R2,R3,R4,其中R1為對(duì)照組,R2,R3,R4為實(shí)驗(yàn)組,進(jìn)水添加三價(jià)鐵濃度分別為2.24,4.48,7.84mg/L,以氯化鐵形式提供.模擬廢水pH值調(diào)節(jié)在6.7~6.8,曝高純氮?dú)?>99.5%)至溶解氧(DO)低于0.5mg/L后,進(jìn)入密封進(jìn)水箱.反應(yīng)器連續(xù)進(jìn)水,溫度通過夾套層控制在(35 ± 1)℃.

    圖1 實(shí)驗(yàn)裝置示意

    1.2 接種污泥與實(shí)驗(yàn)進(jìn)水

    試驗(yàn)污泥取自于南京大學(xué)污染控制與資源再利用國(guó)家重點(diǎn)實(shí)驗(yàn)室的具有良好厭氧氨氧化性能的污水處理反應(yīng)器,污泥含量為3.65g VSS/L.穩(wěn)定運(yùn)行后反應(yīng)器引入有機(jī)物(葡萄糖),并加入三價(jià)鐵.實(shí)驗(yàn)采用模擬廢水,氮素以硫酸銨及亞硝酸鈉提供,濃度按需配制.模擬廢水組分如表1所示,為滿足微生物生長(zhǎng)要求,模擬廢水中加入1.25mL/L微量元素濃縮液,組分參考已發(fā)表文獻(xiàn)[17].

    表1 模擬廢水組分

    1.3 分析項(xiàng)目與檢測(cè)方法

    1.3.1 常規(guī)指標(biāo)測(cè)定與分析 采集進(jìn)水和出水水樣,經(jīng)0.45μm玻璃纖維膜過濾后進(jìn)行分析.NH4+-N, NO2--N,NO3--N,COD采用標(biāo)準(zhǔn)方法測(cè)定(APHA 1998).溶解氧(DO) ,pH值采用便攜式溶氧儀(HQ40d, HACH, America and PB-10, Sartorius, Germany)和pH計(jì)(FE20, METTLER TOLEDO Inc., USA)測(cè)定.

    氨氮,亞硝酸鹽及總氮去除率參照式1,2,3計(jì)算.

    式中:(Inf.NH4+),(Eff.NH4+),(Inf.NO2-),(Eff.NO2-),(Inf.TN),(Eff.TN)分別代表進(jìn)出水中氨氮,亞硝酸鹽及總氮濃度,mg/L.

    1.3.2 微生物特性指標(biāo)檢測(cè) (1)比厭氧氨氧化活性(SAA)測(cè)定:從反應(yīng)器中取一定量污泥加入100mL血清瓶中,使瓶?jī)?nèi)VSS為2.0g/L,加入NH4+-N及NO2--N濃度分別為70和92mg/L的模擬廢水80mL,通入高純氮(>99.5%)15min,保證其無氧環(huán)境.血清瓶置于恒溫振蕩器中,溫度設(shè)定為35℃.用注射器間隔取樣,測(cè)定血清瓶中NH4+-N,NO2--N和NO3--N濃度以確定比厭氧氨氧化活性.(2)氮素轉(zhuǎn)化功能基因及功能菌豐度測(cè)定:取定量污泥,使用試劑盒FastDNATMSpin Kit for soil提取DNA.定量PCR分析選取細(xì)菌16S rRNA,厭氧氨氧化菌16S rRNA及相關(guān)氮素轉(zhuǎn)化功能基因.引物選取參考已發(fā)表文獻(xiàn)[17],并由南京金斯瑞生物科技有限公司合成.反應(yīng)體系為25μL,包括:12.5μL的2×EasyTaq?PCR SuperMix,2μL DNA模板(稀釋至20ng/μL),正反向引物各1μL(20μM)及8.5μL的ddH2O,擴(kuò)增程序循環(huán)數(shù)為25.定量PCR擴(kuò)增在實(shí)時(shí)熒光定量PCR儀AB7500(Life Technologies,美國(guó))上進(jìn)行.(3)微生物群落結(jié)構(gòu)分析:取適量污泥樣品,用FastDNA SPIN Kit for soil試劑盒進(jìn)行DNA提取.采用細(xì)菌通用引物對(duì)污泥樣品DNA進(jìn)行16S rRNA擴(kuò)增,引物序列及擴(kuò)增體系參考已發(fā)表文獻(xiàn)[17].混合4份PCR擴(kuò)增產(chǎn)物,用E.Z.N.A.TMCycle-Pure Kit試劑盒進(jìn)行純化,送至江蘇中宜金大分析檢測(cè)有限公司進(jìn)行分析.所有測(cè)序數(shù)據(jù)均已上傳至Figshare平臺(tái),DOI為10.6084/ m9.figshare.12850043,各樣品的序列信息見表2.

    1.3.3 數(shù)據(jù)分析 運(yùn)用SPSS 22.0軟件對(duì)數(shù)據(jù)進(jìn)行單因素方差分析(ANOVA,置信水平95%),<0.05 視為具有顯著性影響.采用Canoco 軟件進(jìn)行冗余分析(RDA),進(jìn)一步分析細(xì)菌群落相對(duì)豐度和功能基因及功能菌豐度與環(huán)境因子間的相關(guān)性.

    表2 樣品序列信息

    2 結(jié)果與討論

    2.1 有機(jī)物存在下的進(jìn)出水水質(zhì)

    反應(yīng)器中污泥含量為3.65g VSS/L.運(yùn)行前20d,各組反應(yīng)器平行運(yùn)行.20d達(dá)到穩(wěn)定,進(jìn)水添加三價(jià)鐵和有機(jī)物,有機(jī)物按50,100,150,200mg/L COD梯度同步提高.0~60d進(jìn)水氨氮及亞硝酸鹽濃度分別穩(wěn)定為152.5,200mg/L.出水氨氮,亞硝酸鹽,硝酸鹽濃度及其去除率如圖2所示.

    由圖2可知,反應(yīng)器R1在20~30d時(shí),平均氨氮,亞硝酸鹽及總氮去除率分別提高為95.5%,98.1%及95.1%,相較于未添加有機(jī)物,總氮去除率明顯提高;30~40d時(shí),COD濃度提高為100mg/L,反應(yīng)器氮素去除率仍較高;進(jìn)水COD提高為150mg/L時(shí)(40~50d),反應(yīng)器中氨氮去除率降低至83.2%,平均亞硝酸鹽及總氮去除率分別為97.5%及89.6%;而當(dāng)進(jìn)水COD進(jìn)一步提升至200mg/L 時(shí)(50~60d),氨氮去除率下降為61.3%,亞硝酸鹽平均去除率為97.3%,總氮去除率下降為79.8%,出水pH值為8.3 ± 0.2,幾乎超過厭氧氨氧化菌最適生長(zhǎng)pH值范圍.

    有機(jī)物存在下,實(shí)驗(yàn)組R2,R3,R4呈現(xiàn)出與R1類似的運(yùn)行特征,但由于進(jìn)水三價(jià)鐵添加量的不同,反應(yīng)器在不同階段運(yùn)行效能存在一定差異.R2,R3, R4中在同一時(shí)段氮素去除率高于R1.當(dāng)COD濃度提升至200mg/L時(shí), R2,R3,R4的氨氮去除率分別為63.0%,68.2%,71.2%;總氮去除率分別為80.5%, 82.9%,84.7%.

    圖2 有機(jī)物存在下出水氮素濃度及其去除率

    (a) R1; (b) R2; (c) R3; (d) R4;上方濃度為COD梯度,下同

    厭氧氨氧化反應(yīng)器為混菌系統(tǒng).胡勇有等[18]發(fā)現(xiàn)葡萄糖存在下厭氧氨氧化菌與反硝化菌協(xié)同脫氮適宜的COD與氨氮的比值在0~1.57;Winkler等[19]認(rèn)為在進(jìn)水的C/N比小于0.5時(shí),厭氧氨氧化菌能勝過異養(yǎng)反硝化菌.在本實(shí)驗(yàn)中,進(jìn)水COD濃度為50,100mg/L時(shí),各反應(yīng)器中氮素去除率提升明顯,此時(shí)厭氧氨氧化作用起主導(dǎo)作用,反硝化菌協(xié)同降解厭氧氨氧化反應(yīng)產(chǎn)生的硝酸鹽,提高了總氮去除率,其他研究也觀察到了類似現(xiàn)象[6].而COD濃度提升為150mg/L時(shí),出水中氨氮濃度明顯提升,這一現(xiàn)象,在進(jìn)水COD濃度升高為200mg/L時(shí)更明顯,此時(shí)COD與氨氮比值為1.32,出水氨氮濃度顯著升高,同時(shí)亞硝酸鹽及硝酸鹽去除效果明顯,主要的脫氮途徑從厭氧氨氧化轉(zhuǎn)變?yōu)榉聪趸?厭氧氨氧化活性降低的一種可能解釋是,在高有機(jī)物含量下,異養(yǎng)反硝化菌將以更快的速率生長(zhǎng),厭氧氨氧化菌在與反硝化菌競(jìng)爭(zhēng)生存空間和亞硝酸鹽(電子受體)時(shí)占劣勢(shì),同時(shí)反硝化過程使pH值升高,可能超出厭氧氨氧化菌生長(zhǎng)范圍[20].朱葛夫等[21]采用 CSTR 反應(yīng)器研究了以異養(yǎng)反硝化污泥啟動(dòng)厭氧氨氧化系統(tǒng)的可行性,并指出系統(tǒng)的最適C/N比為2.14,與本研究結(jié)果有明顯不同,這可能與采用的反應(yīng)器載體,接種污泥的來源有關(guān),推測(cè)是電子受體競(jìng)爭(zhēng)而不是COD濃度本身導(dǎo)致了厭氧氨氧化活性降低.

    在進(jìn)水存在有機(jī)物的條件下,由于實(shí)驗(yàn)組中三價(jià)鐵的加入,反應(yīng)器在整體運(yùn)行中的氮素轉(zhuǎn)化呈現(xiàn)出一定的不同.在同一時(shí)段,實(shí)驗(yàn)組中氨氮,亞硝酸鹽及總氮去除率略高于對(duì)照組R1,這在R3,R4中體現(xiàn)得更為明顯.有研究表明,三價(jià)鐵作為酶的常見活性劑,有可能提高微生物代謝活性,促進(jìn)其生長(zhǎng)[22],實(shí)驗(yàn)組中較高的氮素去除率,可能因?yàn)槿齼r(jià)鐵的添加,提高了反硝化菌及厭氧氨氧化菌的活性.

    2.2 污泥比厭氧氨氧化活性

    為探究有機(jī)物存在下污泥活性的變化,選取各反應(yīng)器20,30,40,50,60d時(shí)的污泥,測(cè)定其比厭氧氨氧化活性,其結(jié)果如圖3所示.

    各反應(yīng)器運(yùn)行至20d時(shí),由于接種污泥相同且同步運(yùn)行, SAA幾乎一致.20d時(shí)進(jìn)水加入葡萄糖,COD濃度為50mg/L.運(yùn)行至30d時(shí),各反應(yīng)器污泥SAA分別為106.5,115.2,118.6,118.2mg N/(g VSS·d).由出水水質(zhì)可得,由于反應(yīng)器中反硝化作用,總氮去除率有一定升高,同時(shí)厭氧氨氧化活性并未受到抑制,氮素去除效果較好. 40d時(shí), R2,R3,R4污泥SAA分別比R1高7.9,11.9,13.8mg N/(g VSS·d),說明添加三價(jià)鐵提高了污泥厭氧氨氧化活性.

    圖3 有機(jī)物存在時(shí)各反應(yīng)器污泥比厭氧氨氧化活性的變化

    當(dāng)進(jìn)水COD濃度為150和200mg/L時(shí),各反應(yīng)器污泥SAA明顯下降,實(shí)驗(yàn)組中污泥SAA高于對(duì)照組,60d時(shí),R1,R2,R3,R4污泥SAA分別為46.8,48.5, 50.3,52.4mg N/(g VSS·d),可推斷出在高有機(jī)物濃度下,三價(jià)鐵的添加可以提高厭氧氨氧化污泥的活性,增強(qiáng)污泥的有機(jī)物耐受能力.

    由圖3可知,無論進(jìn)水COD濃度如何,同一時(shí)間點(diǎn)各反應(yīng)器中污泥SAA均與三價(jià)鐵添加量呈現(xiàn)出一定的正相關(guān), SAA高低順序?yàn)?R4 > R3 > R2 > R1,與各反應(yīng)器氮去除率結(jié)果一致.

    2.3 微生物群落變化

    由于厭氧氨氧化反應(yīng)器為混菌系統(tǒng),研究有機(jī)物存在下污泥微生物群落的變化,可從微觀層面上解釋三價(jià)鐵對(duì)厭氧氨氧化運(yùn)行過程的影響.本研究選取16S rRNA高通量測(cè)序技術(shù)分析微生物群落結(jié)構(gòu),比較有機(jī)物存在時(shí)不同三價(jià)鐵添加量對(duì)微生物群落結(jié)構(gòu)的影響,探討微生物群落結(jié)構(gòu)變化與脫氮效能之間的關(guān)系.

    2.3.1 有機(jī)物存在下細(xì)菌門水平群落結(jié)構(gòu)變化 圖4為有機(jī)物存在下反應(yīng)器中微生物群落結(jié)構(gòu)門水平分布圖.由圖4得,20d時(shí)各反應(yīng)器中污泥微生物群落結(jié)構(gòu)接近一致,變形菌(Proteobacteria)和綠彎菌(Chloroflexi)為優(yōu)勢(shì)菌門,其豐度分別為34.72%~ 35.25%及19.60%~20.38%.據(jù)報(bào)道,主要的反硝化微生物屬于變形菌(Proteobacteria)[23].綠彎菌(Chloroflexi)可支持生物膜結(jié)構(gòu)的完整性,并且可以在缺氧條件下清除源自厭氧氨氧化菌的代謝產(chǎn)物[24].放線菌(Actinobacteria)及擬桿菌(Bacteroidetes)在各反應(yīng)器中也有一定比例.此外,各反應(yīng)器中與厭氧氨氧化菌相關(guān)的浮霉菌(Planctomycetes)豐度接近,為5.52%~5.54%,這與其他厭氧氨氧化反應(yīng)器的報(bào)道相似[4].

    圖4 有機(jī)物存在下污泥中門水平優(yōu)勢(shì)菌分布

    經(jīng)歷有機(jī)負(fù)荷沖擊后,60d時(shí)反應(yīng)器中門水平群落結(jié)構(gòu)變化明顯. R1,R2,R3,R4中變形菌(Proteobacteria)豐度分別提升至40.16%,39.69%, 44.25%,43.67%;綠彎菌(Chloroflexi)則分別下降至13.53%,10.66%,6.21%,6.23%;放線菌(Actinobacteria)及擬桿菌(Bacteroidetes)均有一定增長(zhǎng).門水平上浮霉菌(Planctomycetes)明顯減少,這種變化也類似于處理高濃度COD廢水的其他厭氧氨氧化菌反應(yīng)器[4].但各反應(yīng)器中浮霉菌(Planctomycetes)減少的程度不同,R1,R2,R3,R4中浮霉菌(Planctomycetes)豐度分別降至1.83%,2.16%,2.24%,2.45%.結(jié)果表明,浮霉菌(Planctomycetes)豐度與所添加三價(jià)鐵濃度呈正相關(guān)性(2=0.962,<0.05),而所有已確認(rèn)的厭氧氨氧化菌均屬于浮霉菌(Planctomycetes)門,其豐度可反映厭氧氨氧化菌的生存狀態(tài)[5],可以推斷三價(jià)鐵的添加促進(jìn)了厭氧氨氧化菌的富集生長(zhǎng).

    2.3.2 有機(jī)物存在下細(xì)菌屬水平群落結(jié)構(gòu)變化 圖5為有機(jī)物存在下污泥屬水平下優(yōu)勢(shì)菌群分布變化.是具有厭氧氨氧化能力的優(yōu)勢(shì)屬,但其豐度隨著進(jìn)水COD的升高而降低,這在R1,R2中體現(xiàn)的更為明顯. 60d時(shí),各反應(yīng)器中其豐度分別為0.83%,0.93%,1.03%,1.13%,可見三價(jià)鐵添加的實(shí)驗(yàn)組中菌富集效果更好,這與反應(yīng)器SAA及運(yùn)行效能一致.可以推斷,在一定濃度范圍內(nèi),厭氧氨氧化菌的富集效果與進(jìn)水鐵離子添加量呈正相關(guān)(2=0.994,<0.01).

    圖5 有機(jī)物存在下污泥中屬水平優(yōu)勢(shì)菌分布(> 0.1%)

    圖6 環(huán)境因子與屬水平優(yōu)勢(shì)細(xì)菌種群間RDA分析

    有機(jī)物存在下,反硝化菌屬,豐度明顯增加.在20d時(shí),幾乎難以檢測(cè)到屬,而在60d時(shí),各反應(yīng)器中屬豐度分別提升為2.43%,2.52%,2.89%,3.12%,與對(duì)照組R1相比,添加三價(jià)鐵的反應(yīng)器中屬豐度較高.反硝化菌屬于-變形菌, 20d時(shí),各反應(yīng)器中其豐度為0.23%~0.53%,而運(yùn)行至60d時(shí),其在各反應(yīng)器中豐度分別提升至4.28%, 4.53%,4.59%,5.34%.說明COD的存在使反硝化菌得以生存和占據(jù)優(yōu)勢(shì)地位.

    運(yùn)用Canoco 4.5軟件,將細(xì)菌群落結(jié)構(gòu)的相對(duì)豐度與環(huán)境因子進(jìn)行冗余分析(圖6).

    由圖6可看出,20d時(shí),各反應(yīng)器無明顯差異,而運(yùn)行60d時(shí),表現(xiàn)出明顯差異,表明三價(jià)鐵的添加對(duì)有機(jī)物存在下細(xì)菌群落影響較大.厭氧氨氧化菌和一些反硝化菌屬,與SAA以及氨氮,亞硝酸鹽及總氮去除率呈正相關(guān),與進(jìn)水COD呈負(fù)相關(guān).而其余反硝化菌屬,,,則相反,即與進(jìn)水COD呈正相關(guān),與SAA以及氨氮,亞硝酸鹽及總氮去除率呈負(fù)相關(guān).表明COD抑制了Anammox過程,而三價(jià)鐵可以通過促進(jìn)豐度的增加,進(jìn)而促進(jìn)氨氮,亞硝酸鹽和總氮的去除.

    2.4 功能基因及功能菌豐度變化

    本研究中,有機(jī)物存在下反應(yīng)器污泥細(xì)菌16S rRNA,厭氧氨氧化菌16S rRNA,氮素轉(zhuǎn)化功能基因及豐度如圖7所示.

    圖7 有機(jī)物存在下amoA基因,hzsB基因,細(xì)菌及厭氧氨氧化菌16S rRNA 豐度變化

    (a)細(xì)菌16S rNRA;(b)厭氧氨氧化菌16S rRNA;(c)基因;(d)基因; “*”表示實(shí)驗(yàn)組與對(duì)照組之間差異顯著(<0.05)

    運(yùn)行至60d時(shí),各反應(yīng)器中細(xì)菌16S rRNA含量分別為3.67×107,3.80×107,4.01×107,4.12×107拷貝數(shù)/ng DNA,反應(yīng)器R3,R4中豐度顯著高于對(duì)照組R1(<0.05), R3,R4中較多的微生物量或是因?yàn)槿齼r(jià)鐵的添加.

    20d時(shí),各反應(yīng)器中厭氧氨氧化菌16S rRNA基因豐度幾乎一致;運(yùn)行至40d時(shí),由于低濃度有機(jī)物的促進(jìn)作用,各反應(yīng)器中厭氧氨氧化菌16S rRNA基因豐度有所提升,R3,R4中豐度顯著高于R1(<0.05),且豐度與三價(jià)鐵添加量呈現(xiàn)出一定的正相關(guān).隨著進(jìn)水COD進(jìn)一步提高,由于較高濃度有機(jī)物下反硝化菌的競(jìng)爭(zhēng)作用,厭氧氨氧化菌含量均有不同程度下降,實(shí)驗(yàn)組中豐度均高于對(duì)照組R1(<0.05).可以發(fā)現(xiàn),厭氧氨氧化菌的豐度與厭氧氨氧化菌的活性有很強(qiáng)的相關(guān)性,這與前人的研究一致[25],因此實(shí)現(xiàn)良好的厭氧氨氧化菌保留率可能是在污水處理系統(tǒng)中增強(qiáng)Anammox的關(guān)鍵.

    厭氧氨氧化功能基因豐度變化與厭氧氨氧化菌16S rRNA基因變化趨勢(shì)一致,均先增加而后減少.同一時(shí)間,R2,R3,R4中的豐度顯著高于R1(<0.05).

    由于反應(yīng)器運(yùn)行過程中嚴(yán)格的條件控制,功能基因在運(yùn)行過程中含量幾乎保持穩(wěn)定,反應(yīng)器R1中其豐度為5.00×102~8.60×102拷貝數(shù)/ng DNA, R2中其豐度保持在5.50×102~1.13×103拷貝數(shù)/ng DNA,R3,R4中其豐度則分別為4.50×102~8.00×102拷貝數(shù)/ng DNA,5.00×102~6.50×102拷貝數(shù)/ng DNA.

    如前所述,有機(jī)物添加使得厭氧氨氧化反應(yīng)器中反硝化菌增長(zhǎng),而整體氮素脫除,依賴于厭氧氨氧化菌及反硝化菌的共同作用.為探究此過程中反硝化作用的變化,選取4種反硝化功能基因,,,測(cè)定,其結(jié)果如圖8所示.

    由圖8得,存在有機(jī)物時(shí),反硝化功能基因,,,的豐度隨著COD的增加而增加,這與前人的研究一致[26].

    運(yùn)行至60d時(shí),基因的豐度在各反應(yīng)器中分別提高至2.13×106,2.65×106,2.76×106,3.01×106拷貝數(shù)/ng DNA,實(shí)驗(yàn)組R2,R3,R4中豐度顯著高于R1(< 0.05).及也呈現(xiàn)類似的變化規(guī)律,說明三價(jià)鐵的添加顯著提高了反硝化功能基因的豐度.

    功能基因及功能菌豐度與環(huán)境因子的冗余分析見圖9.

    圖8 有機(jī)物存在下反硝化功能基因豐度變化

    (a)基因;(b)基因;(c)基因;(d)基因; “*”表示實(shí)驗(yàn)組與對(duì)照組之間差異顯著(<0.05)

    圖9 環(huán)境因子與功能基因及功能菌豐度間RDA分析

    由圖9可看出,20d時(shí),R1,R2,R3,R4無明顯差異,而運(yùn)行60d時(shí)差異明顯,表明三價(jià)鐵的添加使各反應(yīng)器之間功能基因及功能菌豐度的差異變大.厭氧氨氧化菌16S rRNA和厭氧氨氧化功能基因的豐度均與SAA以及氨氮,亞硝酸鹽及總氮去除率呈正相關(guān),與進(jìn)水COD呈負(fù)相關(guān),這也進(jìn)一步驗(yàn)證了當(dāng)COD濃度較高時(shí),厭氧氨氧化菌及厭氧氨氧化功能基因受到抑制,氨氮,亞硝酸鹽及總氮的去除率下降.而功能基因的豐度與SAA以及氨氮,亞硝酸鹽及總氮去除率,進(jìn)水COD之間均無明顯相關(guān)性.反硝化功能基因,,,的豐度與進(jìn)水COD呈正相關(guān),與SAA以及氨氮,亞硝酸鹽及總氮去除率呈負(fù)相關(guān),這也進(jìn)一步驗(yàn)證了有機(jī)物的添加可以提高反硝化功能基因的豐度.

    3 結(jié)論

    3.1 在高濃度的有機(jī)物(200mg/L COD)下,與空白組相比,投加7.84mg/L Fe3+可使氨氮和總氮去除率分別提高9.9%和4.9%.

    3.2 有機(jī)物存在下,污泥微生物群落結(jié)構(gòu)出現(xiàn)變化,表現(xiàn)為厭氧氨氧化菌豐度的降低及反硝化菌屬,豐度的提升.同一有機(jī)物濃度下,添加三價(jià)鐵組的污泥微生物的豐度更高.

    3.3 厭氧氨氧化菌16S rRNA及基因豐度在一定范圍內(nèi)與三價(jià)鐵添加量呈正相關(guān).

    [1] Miao L, Yang G, Tao T, et al. Recent advances in nitrogen removal from landfill leachate using biological treatments - A review [J]. Journal of Environmental Management, 2019,235:178-185.

    [2] 呂 愷,王康舟,姚雪薇,等.基于氨氮,硝氮及乙酸條件下Anammox菌的培養(yǎng) [J]. 中國(guó)環(huán)境科學(xué), 2020,40(10):4348-4353. Lü K, Wang K Z, Yao X W, et al. Enrichment of Anammox under conditions of ammonium, nitrate and acetic acid as substrates [J]. China Environmental Science, 2020,40(10):4348-4353.

    [3] Ma W, Li G, Huang B, et al. Advances and challenges of mainstream nitrogen removal from municipal wastewater with anammox-based processes [J]. Water Environment Research, 2020,92(11):1899-1909.

    [4] Pijuan M, Ribera-Guardia A, Balcazar J L, et al. Effect of COD on mainstream anammox: Evaluation of process performance, granule morphology and nitrous oxide production [J]. Science of the Total Environment, 2020:712:136372.

    [5] Li J, Li J, Gao R, et al. A critical review of one-stage anammox processes for treating industrial wastewater: Optimization strategies based on key functional microorganisms [J]. Bioresource Technology, 2018,265:498-505.

    [6] Chen C, Sun F, Zhang H, et al. Evaluation of COD effect on anammox process and microbial communities in the anaerobic baffled reactor (ABR) [J]. Bioresource Technology, 2016,216:571-578.

    [7] Qin Y, Cao Y, Ren J, et al. Effect of glucose on nitrogen removal and microbial community in anammox-denitrification system [J]. Bioresource Technology, 2017,244:33-39.

    [8] 陳重軍,汪瑤琪,姜 瀅,等.有機(jī)物對(duì)連續(xù)流Anammox脫氮及微生物群落影響 [J]. 中國(guó)環(huán)境科學(xué), 2019,39(12):5049-5055. Chen C J, Wang Y Q, Jiang Y, et al. Effects of organic matter on nitrogen removal and microbial community in anammox reactor [J]. China Environmental Science, 2019,39(12):5049-5055.

    [9] Zhang M, Wang S, Ji B, et al. Towards mainstream deammonification of municipal wastewater: Partial nitrification-anammox versus partial denitrification-anammox [J]. Science of the Total Environment, 2019, 692:393-401.

    [10] Li X, Sun S, Yuan H, et al. Mainstream upflow nitritation-anammox system with hybrid anaerobic pretreatment: Long-term performance and microbial community dynamics [J]. Water Research, 2017,125: 298-308.

    [11] Ferousi C, Lindhoud S, Baymann F, et al. Iron assimilation and utilization in anaerobic ammonium oxidizing bacteria [J]. Current Opinion in Chemical Biology, 2017,37:129-136.

    [12] Van Niftrik L, Geerts W J C, Van Donselaar E G, et al. Combined structural and chemical analysis of the anammoxosome: A membrane- bounded intracytoplasmic compartment in anammox bacteria [J]. Journal of Structural Biology, 2008,161:401-410.

    [13] Liu H, Chen Z, Guan Y, et al. Role and application of iron in water treatment for nitrogen removal: A review [J]. Chemosphere, 2018,204: 51-62.

    [14] Wang J, Liang J, Sun L, et al. Granule-based immobilization and activity enhancement of anammox biomass via PVA/CS and PVA/ CS/Fe gel beads [J]. Bioresource Technology, 2020,309:123448- 123448.

    [15] Yin S, Li J, Dong H, et al. Enhanced nitrogen removal through marine anammox bacteria (MAB) treating nitrogen-rich saline wastewater with Fe(III) addition: Nitrogen shock loading and community structure [J]. Bioresource Technology, 2019:287:121405.

    [16] Wang X, Shu D, Yue H. Taxonomical and functional microbial community dynamics in an Anammox-ASBR system under different Fe (III) supplementation [J]. Applied Microbiology and Biotechnology, 2016,100:10147-10163.

    [17] 毛念佳,嵇鴻民,許 柯,等.鐵離子對(duì)厭氧氨氧化污泥富集培養(yǎng)的影響 [J]. 凈水技術(shù), 2018,37(1):22-29. Mao N J, Ji H M, Xu K, et al. Effect of ferric iron on enrichment cultivation process of anammox sludge [J]. Water Purification Technolugy, 2018,37(1):22-29.

    [18] 胡勇有,梁輝強(qiáng),朱靜平,等.有機(jī)碳源環(huán)境下的厭氧氨氧化批式實(shí)驗(yàn) [J]. 華南理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2007,35(6):116-119.Hu Y Y, Liang H Q, Zhu J P, et al. Batch experiment of anammox under organic carbon source environment [J]. Journal of South China University of Technology (Natural Science Edition), 2007,35(6):116-119.

    [19] Winkler M, Yang J, Kleerebezem R, et al. Nitrate reduction by organotrophic Anammox bacteria in a nitritation/anammox granular sludge and a moving bed biofilm reactor [J]. Bioresource Technology, 2012,114:217-223.

    [20] Du R, Cao S, Zhang H, et al. Flexible nitrite supply alternative for mainstream anammox: Advances in enhancing process stability [J]. Environmental Science & Technology, 2020,54:6353-6364.

    [21] 朱葛夫,張凈瑞,劉超翔,等.厭氧氨氧化工藝的啟動(dòng)及有機(jī)物濃度對(duì)其影響研究 [J]. 環(huán)境工程, 2016,34(2):27-32. Zhu G F, Zhang J R, Liu C X, et al. Study on start-up of anaerobic ammonia oxidation and the effect of COD concentration [J]. Environmental Engineering, 2016,34(2):27-32.

    [22] 丁幫璟,徐夢(mèng)珊,李丹丹,等.不同菜地土壤的鐵氨氧化脫氮過程探究 [J]. 中國(guó)環(huán)境科學(xué), 2020,40(8):3506-3511. Ding B J, Xu M S, Li D D, et al. Investigation of nitrogen removal by feammox in soils from different vegetable fields [J]. China Environmental Science, 2020,40(8):3506-3511.

    [23] Zhang X, Damacharla D, Ma D, et al. Quantitative proteomics reveals novel protein interaction partners of PP2A catalytic subunit in pancreatic beta-cells [J]. Molecular and Cellular Endocrinology, 2016,424:1-11.

    [24] Kindaichi T, Yuri S, Ozaki N, et al. Ecophysiological role and function of uncultured Chloroflexi in an anammox reactor [J]. Water Science and Technology, 2012,66:2556-2561.

    [25] Wang S, Peng Y, Ma B, et al. Anaerobic ammonium oxidation in traditional municipal wastewater treatment plants with low-strength ammonium loading: Widespread but overlooked [J]. Water Research, 2015,84:66-75.

    [26] Wang, Yang R, Guo Y, et al. Investigation of COD and COD/N ratio for the dominance of anammox pathway for nitrogen removal via isotope labelling technique and the relevant bacteria [J]. Journal of Hazardous Materials, 2019,366:606-614.

    Effect of Fe3+on nitrogen removal of Anammox in the presence of organic matter.

    WANG Hai-yue, PENG Ling, MAO Nian-jia, GENG Jin-ju, REN Hong-qiang, XU Ke*

    (State Key Laboratory of Pollution Control and Resource Reuse, School of the Environment, Nanjing University, Nanjing 210023, China)., 2021,41(4):1672~1680

    The tolerance of Anammox to organic matter in the presence of Fe3+(2.24~7.84mg/L) was investigated by running four reactors (R1control, R2 2.24mg/L Fe3+, R3 4.48mg/L Fe3+, R4 7.84mg/L Fe3+). The mechanism was clarified by using 16SrRNA high-throughput sequencing technology and qPCR. The results showed that high NH4+-N and TN removal rates (> 90%) of all four reactors could be achieved at the influent COD concentration of 50 and 100mg/L, and the presence of Fe3+did not show the obvious positive effect. As the influent COD concentration increased to 150 and 200mg/L, the performance of Anammox was inhibited and the positive effect of Fe3+increased. When the influent COD concentration was 200mg/L, the addition of Fe3+(7.84mg/L) increased the NH4+-N and TN removals from 61.3% and 79.8% (R1) to 71.2% and 84.7%, respectively. 16SrRNA high-throughput sequencing results indicated the decrease of Anammox bacteria and the proliferation of denitrifying bacteria in the presence of organic matter. The presence of Fe3+could increase the abundance of Planctomycetes. Fe3+had promoting effects on the abundance of Anammox 16S rRNA and functional geneby qPCR data analysis.

    Fe3+;biological nitrogen removal;Anammox;expanded granular sludge bed (EGSB);COD

    X703

    A

    1000-6923(2021)04-1672-09

    王海月(1997-),女,新疆伊犁人,南京大學(xué)碩士研究生,研究方向?yàn)槲鬯锩摰?

    2020-07-27

    水體污染控制與治理科技重大專項(xiàng)(2017ZX07202003)

    * 責(zé)任作者, 副教授, kexu@nju.edu.cn

    猜你喜歡
    三價(jià)亞硝酸鹽硝化
    羊亞硝酸鹽中毒的病因、臨床表現(xiàn)、診斷與防治措施
    高位池亞硝酸鹽防控
    冬棚養(yǎng)殖需警惕亞硝酸鹽超標(biāo)!一文為你講解亞硝酸鹽過高的危害及處理方法
    MBBR中進(jìn)水有機(jī)負(fù)荷對(duì)短程硝化反硝化的影響
    鋅-鎳合金三價(jià)鉻鈍化工藝的研究
    熱鍍鋅板三價(jià)鉻鈍化膜的耐蝕性分析
    家畜硝酸鹽和亞硝酸鹽中毒的診斷、鑒別和防治
    厭氧氨氧化與反硝化耦合脫氮除碳研究Ⅰ:
    采用12種密度泛函理論方法表征三種三價(jià)鈾復(fù)合物
    海水反硝化和厭氧氨氧化速率同步測(cè)定的15N示蹤法及其應(yīng)用
    国产成人精品婷婷| 丝瓜视频免费看黄片| 嫩草影院精品99| 亚洲伊人久久精品综合| 男人舔奶头视频| 精品国产一区二区三区久久久樱花 | 免费大片黄手机在线观看| 女人被狂操c到高潮| 国产精品嫩草影院av在线观看| 天美传媒精品一区二区| 18禁裸乳无遮挡动漫免费视频 | av女优亚洲男人天堂| 性色avwww在线观看| 一区二区三区四区激情视频| 中国三级夫妇交换| 亚洲一区二区三区欧美精品 | 性插视频无遮挡在线免费观看| 精品久久久久久久久av| 亚洲av.av天堂| 国产极品天堂在线| 我的老师免费观看完整版| 国产精品99久久久久久久久| 青春草视频在线免费观看| 中文字幕免费在线视频6| 韩国高清视频一区二区三区| 国产中年淑女户外野战色| 日韩免费高清中文字幕av| 亚洲国产最新在线播放| 久久久久久久精品精品| 黄色怎么调成土黄色| 搡老乐熟女国产| 久久久精品94久久精品| 国产 一区精品| 1000部很黄的大片| 亚洲成人av在线免费| 少妇高潮的动态图| 精品少妇黑人巨大在线播放| 可以在线观看毛片的网站| 在线观看人妻少妇| 建设人人有责人人尽责人人享有的 | 大片免费播放器 马上看| 日本wwww免费看| 2022亚洲国产成人精品| 神马国产精品三级电影在线观看| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久久丰满| 欧美3d第一页| 乱系列少妇在线播放| 午夜爱爱视频在线播放| 国内精品宾馆在线| 一区二区三区乱码不卡18| 啦啦啦中文免费视频观看日本| 精品一区在线观看国产| 日韩三级伦理在线观看| 视频区图区小说| 国产精品成人在线| 亚洲精品亚洲一区二区| 1000部很黄的大片| 亚洲电影在线观看av| h日本视频在线播放| 天天一区二区日本电影三级| 久久女婷五月综合色啪小说 | 国产免费一级a男人的天堂| 精品国产露脸久久av麻豆| 国产男女超爽视频在线观看| 看十八女毛片水多多多| 亚洲av二区三区四区| 国产成人精品福利久久| 国产成人精品福利久久| 国产精品一区二区在线观看99| 欧美日韩综合久久久久久| 少妇高潮的动态图| 看十八女毛片水多多多| 伦理电影大哥的女人| 高清午夜精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 成人欧美大片| 日韩欧美一区视频在线观看 | 噜噜噜噜噜久久久久久91| 一边亲一边摸免费视频| 久久人人爽人人爽人人片va| 夫妻午夜视频| 久久ye,这里只有精品| 国产精品女同一区二区软件| 欧美zozozo另类| 嫩草影院入口| 国语对白做爰xxxⅹ性视频网站| 亚洲综合精品二区| 亚洲精品日本国产第一区| 狂野欧美激情性bbbbbb| 婷婷色av中文字幕| 午夜精品国产一区二区电影 | 亚洲国产日韩一区二区| 欧美最新免费一区二区三区| 亚洲色图av天堂| 免费观看av网站的网址| 我的老师免费观看完整版| 亚洲av欧美aⅴ国产| 午夜日本视频在线| freevideosex欧美| 日本午夜av视频| 日韩中字成人| 久久久久精品久久久久真实原创| 日韩不卡一区二区三区视频在线| 亚洲怡红院男人天堂| 国产伦精品一区二区三区视频9| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品乱久久久久久| 自拍欧美九色日韩亚洲蝌蚪91 | 久久99热6这里只有精品| 一级av片app| 直男gayav资源| videos熟女内射| 涩涩av久久男人的天堂| 在线观看一区二区三区激情| 午夜精品一区二区三区免费看| 欧美变态另类bdsm刘玥| 中国美白少妇内射xxxbb| 久久亚洲国产成人精品v| 久久亚洲国产成人精品v| 亚洲国产精品999| 毛片一级片免费看久久久久| 欧美性感艳星| 黄色怎么调成土黄色| 大片免费播放器 马上看| 亚洲最大成人手机在线| 简卡轻食公司| 女人被狂操c到高潮| 欧美成人精品欧美一级黄| 国产在线一区二区三区精| 观看免费一级毛片| 国产成人精品一,二区| 亚洲丝袜综合中文字幕| 欧美高清成人免费视频www| 欧美+日韩+精品| 国产毛片a区久久久久| 欧美最新免费一区二区三区| 两个人的视频大全免费| 男女边摸边吃奶| 国产 一区精品| 国产亚洲5aaaaa淫片| 国产欧美亚洲国产| 各种免费的搞黄视频| 禁无遮挡网站| 久久精品熟女亚洲av麻豆精品| 在线观看美女被高潮喷水网站| 国产男人的电影天堂91| 亚洲自偷自拍三级| 精品一区在线观看国产| 观看美女的网站| 欧美激情在线99| 男插女下体视频免费在线播放| 日韩 亚洲 欧美在线| 人妻少妇偷人精品九色| 久久精品夜色国产| 一级av片app| 亚洲图色成人| 成年女人在线观看亚洲视频 | 国产成年人精品一区二区| 在线观看国产h片| 欧美日韩视频高清一区二区三区二| 久久综合国产亚洲精品| 久久久久久久精品精品| 欧美成人一区二区免费高清观看| 在现免费观看毛片| 三级国产精品片| 如何舔出高潮| 国产亚洲av片在线观看秒播厂| videos熟女内射| 青春草视频在线免费观看| 插逼视频在线观看| 91午夜精品亚洲一区二区三区| 国产精品偷伦视频观看了| 少妇人妻 视频| 久久人人爽人人爽人人片va| 在线观看美女被高潮喷水网站| 久久国内精品自在自线图片| 欧美丝袜亚洲另类| 中文字幕亚洲精品专区| 老师上课跳d突然被开到最大视频| 小蜜桃在线观看免费完整版高清| 国产一区二区三区综合在线观看 | 青青草视频在线视频观看| 中国国产av一级| 久久国内精品自在自线图片| a级毛色黄片| 成人二区视频| 直男gayav资源| 赤兔流量卡办理| 99精国产麻豆久久婷婷| 久久久久久九九精品二区国产| 日韩欧美 国产精品| 22中文网久久字幕| 婷婷色麻豆天堂久久| 汤姆久久久久久久影院中文字幕| 欧美最新免费一区二区三区| 一个人看视频在线观看www免费| 波多野结衣巨乳人妻| 国产精品99久久久久久久久| 亚洲色图av天堂| 午夜福利网站1000一区二区三区| 国产黄频视频在线观看| 精华霜和精华液先用哪个| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品午夜福利在线看| 久久精品综合一区二区三区| 久久久久久久精品精品| 禁无遮挡网站| 国产69精品久久久久777片| 国产在线一区二区三区精| 欧美高清成人免费视频www| 最近中文字幕2019免费版| 三级国产精品欧美在线观看| 成人毛片a级毛片在线播放| 日韩av在线免费看完整版不卡| 哪个播放器可以免费观看大片| 免费看光身美女| 亚洲国产精品国产精品| av在线app专区| 国产色婷婷99| 特级一级黄色大片| 国产精品久久久久久久电影| 欧美日韩在线观看h| 国产国拍精品亚洲av在线观看| 69av精品久久久久久| 免费黄色在线免费观看| 亚洲不卡免费看| 91狼人影院| 日韩人妻高清精品专区| 久久久久久伊人网av| 免费大片18禁| 欧美精品国产亚洲| 精品人妻视频免费看| 亚洲精品色激情综合| 一二三四中文在线观看免费高清| av女优亚洲男人天堂| 亚洲国产精品国产精品| 亚洲国产日韩一区二区| 成人毛片a级毛片在线播放| 最新中文字幕久久久久| 亚洲精品国产成人久久av| 日韩免费高清中文字幕av| 晚上一个人看的免费电影| 看黄色毛片网站| 大陆偷拍与自拍| 日韩中字成人| 午夜亚洲福利在线播放| 熟女人妻精品中文字幕| 高清在线视频一区二区三区| 在线精品无人区一区二区三 | 色5月婷婷丁香| 纵有疾风起免费观看全集完整版| 少妇人妻久久综合中文| av在线播放精品| 性色av一级| 久热这里只有精品99| 老司机影院毛片| av免费观看日本| 欧美性猛交╳xxx乱大交人| 最近的中文字幕免费完整| 亚洲天堂国产精品一区在线| 九草在线视频观看| 亚洲一区二区三区欧美精品 | tube8黄色片| 九草在线视频观看| 欧美bdsm另类| 毛片一级片免费看久久久久| 美女高潮的动态| 波多野结衣巨乳人妻| 精品人妻偷拍中文字幕| 免费不卡的大黄色大毛片视频在线观看| 亚洲,一卡二卡三卡| 欧美3d第一页| 欧美老熟妇乱子伦牲交| 日韩欧美 国产精品| 日韩三级伦理在线观看| 丰满少妇做爰视频| 亚洲人成网站在线播| www.av在线官网国产| 亚洲激情五月婷婷啪啪| 大码成人一级视频| 毛片女人毛片| 久久久久久久亚洲中文字幕| av天堂中文字幕网| 国产成人91sexporn| 亚洲精品影视一区二区三区av| 特大巨黑吊av在线直播| 国产成人一区二区在线| 日韩av不卡免费在线播放| 亚洲国产欧美人成| 午夜福利视频1000在线观看| 新久久久久国产一级毛片| 欧美日韩综合久久久久久| 亚洲国产精品专区欧美| 毛片一级片免费看久久久久| 精品99又大又爽又粗少妇毛片| 国产精品不卡视频一区二区| 汤姆久久久久久久影院中文字幕| 嫩草影院精品99| 韩国高清视频一区二区三区| 国产成人一区二区在线| 日本与韩国留学比较| 一级二级三级毛片免费看| 韩国av在线不卡| 国产精品不卡视频一区二区| 欧美亚洲 丝袜 人妻 在线| 欧美另类一区| 51国产日韩欧美| 国产午夜精品一二区理论片| 日韩av在线免费看完整版不卡| 国产黄a三级三级三级人| 久久久久精品性色| 免费人成在线观看视频色| 日本免费在线观看一区| 久久久久久久久大av| 性插视频无遮挡在线免费观看| 尤物成人国产欧美一区二区三区| 精品酒店卫生间| av国产免费在线观看| 久久久久久久精品精品| 一级av片app| 亚洲最大成人中文| 亚洲精品日韩av片在线观看| 偷拍熟女少妇极品色| 天堂俺去俺来也www色官网| 久久久欧美国产精品| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 在线免费十八禁| 成人一区二区视频在线观看| 边亲边吃奶的免费视频| 精品久久久久久久人妻蜜臀av| 国产成人精品久久久久久| 日本wwww免费看| 一区二区三区免费毛片| 亚洲图色成人| 最近中文字幕2019免费版| 99热网站在线观看| 好男人视频免费观看在线| 亚洲天堂国产精品一区在线| 亚洲第一区二区三区不卡| 丰满人妻一区二区三区视频av| 亚洲四区av| 国产综合懂色| 全区人妻精品视频| 又粗又硬又长又爽又黄的视频| 国产精品久久久久久精品电影小说 | 2018国产大陆天天弄谢| 精品一区二区三区视频在线| 成年女人在线观看亚洲视频 | 成人无遮挡网站| 乱码一卡2卡4卡精品| 黄色日韩在线| 成人国产av品久久久| 天堂俺去俺来也www色官网| 日韩一区二区视频免费看| 国产午夜精品久久久久久一区二区三区| 国产高清有码在线观看视频| 欧美成人午夜免费资源| 99视频精品全部免费 在线| 777米奇影视久久| 午夜爱爱视频在线播放| 欧美日本视频| 亚洲熟女精品中文字幕| 久久午夜福利片| av专区在线播放| 亚洲天堂av无毛| 18禁裸乳无遮挡免费网站照片| 亚洲天堂国产精品一区在线| 亚洲av日韩在线播放| 国产精品爽爽va在线观看网站| 国产精品久久久久久久久免| 国产黄片美女视频| 亚洲一级一片aⅴ在线观看| 老司机影院成人| 欧美老熟妇乱子伦牲交| 亚洲真实伦在线观看| av国产精品久久久久影院| 夜夜爽夜夜爽视频| 男人舔奶头视频| 午夜福利视频精品| 国产免费又黄又爽又色| 久久99精品国语久久久| 卡戴珊不雅视频在线播放| videos熟女内射| 久久精品人妻少妇| 天美传媒精品一区二区| 伊人久久国产一区二区| 大话2 男鬼变身卡| 我要看日韩黄色一级片| 国产91av在线免费观看| 欧美精品人与动牲交sv欧美| 在线观看免费高清a一片| 晚上一个人看的免费电影| 国产老妇女一区| 一级片'在线观看视频| 国产乱来视频区| 白带黄色成豆腐渣| 亚洲av.av天堂| 国产黄色视频一区二区在线观看| 免费av观看视频| 久久久久精品久久久久真实原创| 男插女下体视频免费在线播放| 成人午夜精彩视频在线观看| 哪个播放器可以免费观看大片| 精品久久久精品久久久| 国产又色又爽无遮挡免| 久久国产乱子免费精品| 插阴视频在线观看视频| 18+在线观看网站| 亚洲美女视频黄频| 精品人妻熟女av久视频| 九草在线视频观看| 亚洲精品色激情综合| 热re99久久精品国产66热6| 亚洲经典国产精华液单| 久久人人爽人人爽人人片va| 国产精品久久久久久久久免| 免费观看av网站的网址| 中国国产av一级| 国产在线一区二区三区精| 中文字幕av成人在线电影| 一区二区三区四区激情视频| 国产 一区 欧美 日韩| 国产黄色免费在线视频| 男男h啪啪无遮挡| 欧美变态另类bdsm刘玥| 久久99蜜桃精品久久| 性色avwww在线观看| 精品久久久久久久人妻蜜臀av| 看非洲黑人一级黄片| 在线 av 中文字幕| 亚洲欧美日韩东京热| 日韩三级伦理在线观看| 热re99久久精品国产66热6| 一级毛片黄色毛片免费观看视频| 国产亚洲最大av| 伊人久久国产一区二区| 日日啪夜夜撸| 春色校园在线视频观看| kizo精华| 亚洲综合色惰| 麻豆乱淫一区二区| 波野结衣二区三区在线| 国产视频首页在线观看| 精品人妻一区二区三区麻豆| 免费高清在线观看视频在线观看| 国产高清三级在线| 国产爽快片一区二区三区| 国产永久视频网站| 日韩欧美精品免费久久| 人人妻人人看人人澡| 免费av毛片视频| 天堂俺去俺来也www色官网| 高清午夜精品一区二区三区| 国产高清国产精品国产三级 | 久热久热在线精品观看| 永久网站在线| 午夜亚洲福利在线播放| 男人添女人高潮全过程视频| 精品国产乱码久久久久久小说| 精品一区二区三区视频在线| 看免费成人av毛片| 欧美国产精品一级二级三级 | 少妇人妻一区二区三区视频| 春色校园在线视频观看| 内射极品少妇av片p| 国产av不卡久久| 久久久久久久国产电影| 亚洲精品乱码久久久久久按摩| 天天躁日日操中文字幕| 久久亚洲国产成人精品v| 高清日韩中文字幕在线| 国产精品蜜桃在线观看| 男的添女的下面高潮视频| 大片免费播放器 马上看| 久久久欧美国产精品| 国产成年人精品一区二区| 晚上一个人看的免费电影| 边亲边吃奶的免费视频| videossex国产| 国产精品不卡视频一区二区| 777米奇影视久久| 免费黄网站久久成人精品| 国产亚洲5aaaaa淫片| 毛片女人毛片| 白带黄色成豆腐渣| 久久影院123| 青春草视频在线免费观看| 国内揄拍国产精品人妻在线| 韩国高清视频一区二区三区| 三级国产精品欧美在线观看| 97热精品久久久久久| 黄色一级大片看看| 18禁在线播放成人免费| 国产欧美日韩精品一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av成人精品一区久久| 欧美日韩精品成人综合77777| 国产乱来视频区| av.在线天堂| 日本av手机在线免费观看| 亚洲图色成人| 日本-黄色视频高清免费观看| 国产精品人妻久久久影院| 黄色视频在线播放观看不卡| 又黄又爽又刺激的免费视频.| 国产黄片美女视频| 菩萨蛮人人尽说江南好唐韦庄| 最近中文字幕高清免费大全6| 久热久热在线精品观看| 国产免费一级a男人的天堂| av在线播放精品| 国产精品人妻久久久久久| 97超视频在线观看视频| 熟妇人妻不卡中文字幕| freevideosex欧美| 日韩免费高清中文字幕av| 国产一区二区三区综合在线观看 | 久久精品国产亚洲av涩爱| 国产精品人妻久久久影院| 九九爱精品视频在线观看| 精品人妻视频免费看| 免费黄色在线免费观看| 交换朋友夫妻互换小说| 欧美日韩国产mv在线观看视频 | 欧美变态另类bdsm刘玥| 亚洲自偷自拍三级| 日日啪夜夜爽| 视频区图区小说| 男人狂女人下面高潮的视频| 亚洲激情五月婷婷啪啪| 老司机影院成人| 国产有黄有色有爽视频| 男男h啪啪无遮挡| 国产探花极品一区二区| 99九九线精品视频在线观看视频| 国产在线男女| 一级a做视频免费观看| 亚洲精品国产av蜜桃| 欧美日韩在线观看h| a级毛色黄片| 国产黄色免费在线视频| 亚洲国产色片| 三级国产精品欧美在线观看| videossex国产| 黄色日韩在线| 性插视频无遮挡在线免费观看| 人妻制服诱惑在线中文字幕| 久久人人爽av亚洲精品天堂 | 成人综合一区亚洲| 亚州av有码| 国产一区有黄有色的免费视频| 日本wwww免费看| 精品99又大又爽又粗少妇毛片| 亚洲欧美日韩无卡精品| 99热这里只有是精品在线观看| 国产精品国产三级专区第一集| 秋霞在线观看毛片| 22中文网久久字幕| 亚洲成人久久爱视频| 在线观看av片永久免费下载| 好男人在线观看高清免费视频| 国产午夜福利久久久久久| 男插女下体视频免费在线播放| 高清午夜精品一区二区三区| 王馨瑶露胸无遮挡在线观看| 黄色日韩在线| av卡一久久| 成人亚洲精品一区在线观看 | 免费播放大片免费观看视频在线观看| 国产精品久久久久久久久免| 欧美bdsm另类| 女的被弄到高潮叫床怎么办| 精品久久久久久久久亚洲| 中文字幕制服av| 亚洲欧美成人综合另类久久久| 久久99热6这里只有精品| 免费看a级黄色片| 建设人人有责人人尽责人人享有的 | 欧美日韩国产mv在线观看视频 | 伊人久久精品亚洲午夜| 我要看日韩黄色一级片| 国产成人一区二区在线| 超碰av人人做人人爽久久| 高清视频免费观看一区二区| 亚洲最大成人av| 高清视频免费观看一区二区| 亚洲精品久久午夜乱码| 欧美激情久久久久久爽电影| 91在线精品国自产拍蜜月| 岛国毛片在线播放| 69av精品久久久久久| 看十八女毛片水多多多| 18禁在线无遮挡免费观看视频| 99视频精品全部免费 在线| 亚洲综合精品二区| 亚洲电影在线观看av| 午夜福利视频精品| 插逼视频在线观看| 亚洲最大成人手机在线| 男的添女的下面高潮视频| 精品国产三级普通话版| 啦啦啦啦在线视频资源| 欧美日韩在线观看h| 欧美 日韩 精品 国产| 精品国产三级普通话版| 极品少妇高潮喷水抽搐| 五月玫瑰六月丁香| 日韩电影二区| 男女国产视频网站| 免费看光身美女| 色播亚洲综合网| 人体艺术视频欧美日本| 亚洲最大成人手机在线|