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

    沸石載體恢復(fù)受饑餓影響厭氧氨氧化菌的性能研究

    2021-07-16 07:05:48余道道孫敬起霍唐燃趙云鵬劉思彤
    關(guān)鍵詞:沸石饑餓反應(yīng)器

    余道道 孫敬起 霍唐燃 趙云鵬 劉思彤,,?

    沸石載體恢復(fù)受饑餓影響厭氧氨氧化菌的性能研究

    余道道1孫敬起2霍唐燃2趙云鵬2劉思彤1,2,?

    1.北京大學(xué)深圳研究生院環(huán)境與能源學(xué)院, 深圳 518055; 2.北京大學(xué)環(huán)境科學(xué)與工程學(xué)院, 北京 100871; ?通信作者, E-mail: liusitong@pku.edu.cn

    利用沸石良好的吸氨能力和微生物載體功能, 將受饑餓影響的厭氧氨氧化菌在含沸石的反應(yīng)器內(nèi)進(jìn)行活性恢復(fù), 研究沸石對(duì)長(zhǎng)期饑餓影響下厭氧氨氧化菌恢復(fù)過(guò)程中菌群活性及結(jié)構(gòu)的影響。結(jié)果表明, 經(jīng)過(guò) 82 天的運(yùn)行, 沸石反應(yīng)器和對(duì)照反應(yīng)器(無(wú)沸石添加)內(nèi)厭氧氨氧化菌均獲得恢復(fù)。沸石反應(yīng)器的脫氮效率從第 59 天開始與對(duì)照反應(yīng)器拉開差距, 運(yùn)行第 82 天時(shí), 沸石反應(yīng)器中氮去除負(fù)荷(NRR)為 177.8 mgN/(L·d), 顯著大于對(duì)照反應(yīng)器 (154.3mgN/(L·d))。利用 16SrRNA 高通量測(cè)序方法, 對(duì)反應(yīng)器內(nèi)菌群的微生物組成與功能進(jìn)行分析, 發(fā)現(xiàn)投加沸石載體使得恢復(fù)后的厭氧氨氧化菌的豐度更高, 可以更好地促進(jìn)受饑餓抑制的厭氧氨氧化菌的恢復(fù)。

    厭氧氨氧化菌; 沸石; 細(xì)菌恢復(fù); 群落結(jié)構(gòu)

    厭氧氨氧化是在厭氧條件下, 以氨為電子供體, 亞硝酸鹽為電子受體, 將氨氧化成氮?dú)獾纳镞^(guò)程。與傳統(tǒng)的生物脫氮方法相比, 厭氧氨氧化具有氮去除負(fù)荷高、污泥產(chǎn)量低以及無(wú)需外加碳源等優(yōu)點(diǎn)[1], 受到學(xué)者們的普遍關(guān)注。然而, 由于厭氧氨氧化菌生長(zhǎng)速度極慢, 倍增時(shí)間長(zhǎng), 因此其培養(yǎng)過(guò)程極其緩慢[2–3]。更重要的是, 在反應(yīng)器長(zhǎng)期運(yùn)行之后, 如果種泥保存不當(dāng), 使污泥處于饑餓狀態(tài), 活性會(huì)受到更大的影響。因此, 研究長(zhǎng)期受饑餓影響厭氧氨氧化菌的有效恢復(fù)策略具有重要的意義。

    目前, 恢復(fù)受饑餓影響厭氧氨氧化菌活性的措施主要有改變進(jìn)水負(fù)荷、投加中間產(chǎn)物 N2H4以及添加陶?;蝈缧菭罾w維材料等載體。候曉幫等[4]采用 4 種方式(正常負(fù)荷、降低負(fù)荷、投加 N2H4和投加 NH2OH)復(fù)活厭氧氨氧化菌, 發(fā)現(xiàn)降低負(fù)荷或投加 N2H4是快速恢復(fù)厭氧氨氧化菌活性的有效方法。楊開亮等[5]采用以彗星狀纖維材料為載體(特點(diǎn)是利用中間的小珠將纖維固定起來(lái), 兩端呈放射狀)的序批式反應(yīng)器(sequencing batch reactor, SBR), 對(duì)停止運(yùn)行 100 天的反應(yīng)器中厭氧氨氧化菌進(jìn)行恢復(fù)培養(yǎng)。結(jié)果表明, 生物膜的形成可增強(qiáng)厭氧氨氧化菌的活性, 縮短恢復(fù)周期, 恢復(fù)后厭氧氨氧化菌的脫氮性能良好。但是, 這些方法都存在不足之處: 外源添加藥品增加成本, 且 N2H4有劇毒、極不穩(wěn)定; 采用陶?;蝈缧菭罾w維材料等載體存在成本高、恢復(fù)過(guò)程漫長(zhǎng)等問(wèn)題。因此, 本研究擬用一種高效吸附 NH4+-N、經(jīng)濟(jì)成本低且環(huán)境友好的材料——沸石, 通過(guò)底物吸附耦合生物降解的模式, 對(duì)受饑餓影響的厭氧氨氧化菌進(jìn)行恢復(fù)。

    沸石是分布在自然界中的硅酸鹽礦物, 內(nèi)部具有豐富的多孔結(jié)構(gòu), 孔體積可達(dá)自身體積的一半, 比表面積高達(dá) 400~1000m2/g, 能提供多個(gè)吸附活性位點(diǎn)[6]。沸石具有良好的離子吸附性及交換性, 尤其對(duì) NH4+-N 的吸附性能良好[7]。Miladinovic 等[8]發(fā)現(xiàn)表面有生物膜的沸石比普通沸石脫氮效果好, 沸石生物膜上的硝化作用可以協(xié)同沸石的離子交換作用共同去除 NH4+-N。Qiu 等[9]分別以陶瓷、沸石和碳酸鹽為填料, 研究 3 種曝氣生物濾池(biological aerated filter, BAF)的性能, 發(fā)現(xiàn)沸石 BAF 和碳酸鹽 BAF 的 NH4+-N 去除效率比陶瓷顆粒 BAF 高, 且抗沖擊能力更好。Yapsakli 等[10]利用沸石吸附和生物降解相結(jié)合的方法去除廢水中 NH4+-N, 原理是將微生物的代謝作用與沸石的 NH4+-N 吸附作用耦合, 通過(guò)沸石吸附與解吸附 NH4+-N 離子的功能, 為厭氧氨氧化反應(yīng)提供 NO2?-N 與 NH4+-N 之比相對(duì)穩(wěn)定的進(jìn)水, 促進(jìn)廢水中 NH4+-N 的去除?,F(xiàn)有研究多以沸石為載體來(lái)促進(jìn)廢水中 NH4+-N 的去除, 鮮有研究利用沸石促進(jìn)受饑餓影響厭氧氨氧化菌的恢復(fù)。我們通過(guò)比較在反應(yīng)器內(nèi)加入和不加沸石兩種情況下長(zhǎng)期受饑餓影響厭氧氨氧化菌的活性恢復(fù)情況, 利用 16S rRNA 高通量測(cè)序, 探究恢復(fù)過(guò)程中菌群變化規(guī)律, 為沸石用于恢復(fù)受饑餓影響厭氧氨氧化菌的研究和應(yīng)用奠定基礎(chǔ)。

    1 材料和方法

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

    實(shí)驗(yàn)裝置由進(jìn)水系統(tǒng)、曝氣系統(tǒng)、反應(yīng)器主體、出水系統(tǒng)、水浴循環(huán)系統(tǒng)以及自控系統(tǒng)組成, 反應(yīng)器有效容積為 5L。通過(guò)時(shí)間控制器控制反應(yīng)器的進(jìn)水、攪拌、沉降與出水階段的時(shí)間。兩個(gè)反應(yīng)器頂部均有硅膠墊圈, 以此隔絕氧氣。在反應(yīng)器底部放置轉(zhuǎn)子, 并通過(guò)控制磁力攪拌器進(jìn)行攪拌。為了保證厭氧氨氧化菌所需的恒溫條件, 反應(yīng)器外層設(shè)有溫度控制在 37±0.5oC 的水浴夾層。其中一個(gè)反應(yīng)器為不加沸石的對(duì)照, 另一個(gè)反應(yīng)器底部放入一個(gè)自制三層沸石架(圖 1)。沸石架由 3 只定制不銹鋼架焊接而成, 單只架外圈直徑為 11cm, 高 5cm。將每只不銹鋼架網(wǎng)面的中間掏空, 從而使攪拌過(guò)程中厭氧氨氧化菌與沸石充分接觸。網(wǎng)面被掏空的部分呈圓形, 內(nèi)圈直徑為 4cm, 沸石架高約 18cm, 沸石反應(yīng)器內(nèi)的沸石添加量約為 100g/L。整個(gè)沸石反應(yīng)器裝置如圖 2 所示。

    1.2 實(shí)驗(yàn)廢水和接種污泥

    反應(yīng)器進(jìn)水采用人工模擬廢水, 其物質(zhì)組成如表 1 所示。進(jìn)水前, 用 N2/CO2(95%/5%)對(duì)人工模擬廢水進(jìn)行曝氣, 除去水中的溶解氧, 并將進(jìn)水 pH 控制在 6.8~7.2。

    反應(yīng)器接種的是放置在室溫中試反應(yīng)器內(nèi)、在NH4+-N 和 NO2?-N 濃度接近 0mg/L 的污水中饑餓長(zhǎng)達(dá) 88 天的厭氧氨氧化污泥。接種后, 各反應(yīng)器內(nèi)混合液揮發(fā)性懸浮固體(mixed liquor volatile suspended solid, MLVSS)的濃度均為4500 mg/L。

    1.3 反應(yīng)器運(yùn)行和水質(zhì)檢測(cè)方法

    對(duì)照反應(yīng)器和沸石反應(yīng)器共運(yùn)行 82 天, 第 1~10天為啟動(dòng)階段, 第 11~30 天為負(fù)荷提升階段, 第 31~82 天為穩(wěn)定運(yùn)行階段。反應(yīng)器啟動(dòng)初期, 兩個(gè)反應(yīng)器進(jìn)水的 NH4+-N 和 NO2?-N 濃度均保持在 50±5 mg/L。在負(fù)荷提升階段結(jié)束以前, 當(dāng)連續(xù) 3 個(gè)水力停留時(shí)間(hydraulic retention time, HRT)反應(yīng)器出水中 NO2?-N 濃度都穩(wěn)定在 5mg/L 以下時(shí), 通過(guò)提高進(jìn)水 NH4+-N 和 NO2?-N 濃度的方式來(lái)提升反應(yīng)器進(jìn)水氮負(fù)荷。當(dāng)進(jìn)水的 NH4+-N 和 NO2?-N 濃度提高至 200mg/L 左右時(shí), 即負(fù)荷提升階段結(jié)束后, 不再繼續(xù)提高進(jìn)水中 NH4+-N 和 NO2?-N 的濃度, 改為通過(guò)縮短 HRT 的方式來(lái)加大反應(yīng)器進(jìn)水的氮負(fù)荷。在穩(wěn)定運(yùn)行階段末期, 兩個(gè)反應(yīng)器的 HRT 均由最初的 24 小時(shí) 縮短至 8 小時(shí)。

    圖1 裝滿沸石的沸石架

    圖2 沸石反應(yīng)器裝置示意圖

    表1 合成廢水的組成

    每天取兩個(gè)反應(yīng)器的進(jìn)出水進(jìn)行水質(zhì)檢測(cè)分析。在反應(yīng)器運(yùn)行的 3 個(gè)階段, 分別取兩個(gè)反應(yīng)器內(nèi)懸浮態(tài)污泥樣品以及沸石反應(yīng)器內(nèi)沸石表面附著態(tài)污泥樣品, 用于分析胞外聚合物(extracellular po-lymeric substance, EPS)和微生物的組成結(jié)構(gòu), 樣品采集時(shí)間和編號(hào)如表 2 所示。

    使用梅特勒 FE20 型 pH 計(jì)測(cè)定 pH 值, 使用雷磁 JPB-607A 型便攜式溶解氧儀測(cè)定溶解氧。氨氮(NH4+-N)、亞硝酸鹽氮(NO2?-N)、硝酸鹽氮(NO3?- N)、總氮(total nitrogen, TN)和 MLVSS 等指標(biāo)均按照國(guó)家標(biāo)準(zhǔn)方法[11]測(cè)定。

    表2 對(duì)照反應(yīng)器及沸石反應(yīng)器污泥樣品編號(hào)

    說(shuō)明: “?”表示反應(yīng)器運(yùn)行的第1天, 無(wú)附著態(tài)污泥。

    1.4 EPS 提取及測(cè)定

    EPS 是特定環(huán)境條件下細(xì)菌代謝分泌的、包裹在細(xì)胞壁外的高分子聚合物, 主要由蛋白質(zhì)(protein, PN)和多糖(polysaccharide, PS)構(gòu)成[12]。EPS 是顆粒污泥快速形成與維持穩(wěn)定的關(guān)鍵因素之一[13], 厭氧氨氧化污泥中的 EPS 含量越大, 污泥顆粒穩(wěn)定性越強(qiáng)。本研究采用熱提取法[14]提取 EPS。取 0.5mL 厭氧氨氧化污泥, 靜置 20 分鐘后去除上清液。用0.05%NaCl 溶液將樣品定容至 10mL, 將混合液放入 60oC 水浴中加熱 30 分鐘。在 10000rpm 的轉(zhuǎn)速下將混合液離心 15 分鐘, 吸取上清液, 并用 0.45μm濾膜過(guò)濾, 即可得到 EPS 提取液。采用快速 Lowry法蛋白質(zhì)試劑盒測(cè)定 EPS 中 PN 含量, 采用蒽酮硫酸比色法測(cè)定 EPS 中 PS 含量[15]。采用標(biāo)準(zhǔn)重量法測(cè)定可揮發(fā)性懸浮固體濃度(volatile suspended solids, VSS)[11], 利用公式: (測(cè)得 PN/PS 濃度×提取液體積)/VSS, 將 PN 與 PS 含量相加, 即得最終 EPS濃度, EPS 中 PN 和 PS 以每克 VSS 中的蛋白和多糖含量計(jì)算。

    在整個(gè) 82 天的運(yùn)行周期中, 采集并測(cè)定第 1天、第 15 天、第 30 天、第 45 和第 82 天共 5 個(gè)時(shí)間點(diǎn)的污泥樣本。由于測(cè)定 EPS 需至少 3 個(gè)平行樣, 每個(gè)平行樣取 0.5mL 菌樣, 而起始階段沸石反應(yīng)器內(nèi)采集到的附著態(tài)污泥過(guò)少(低于 1.5mL), 故只測(cè)定對(duì)照反應(yīng)器內(nèi)懸浮樣及沸石反應(yīng)器內(nèi)懸浮污泥的EPS含量。

    1.5 DNA 提取、聚合酶鏈?zhǔn)綌U(kuò)增反應(yīng)(PCR)及Illumina測(cè)序

    采用試劑盒法提取污泥中的 DNA, 所用的試劑盒為 FastDNA? SPIN Kit for Soil 土壤 DNA 快速提取試劑盒, 采用瓊脂糖凝膠電泳對(duì) DNA 提取情況進(jìn)行檢驗(yàn)。提取出的 DNA 樣品置于?20oC 冰箱內(nèi)冷凍保存。利用細(xì)菌通用引物 338F (ACTCCTACGGGA GGCAGCAG)和 806R (GGACTACHVGGGTWTCTA AT)[16]對(duì)合格的 DNA 樣品進(jìn)行 16S rRNA 基因在V3~V4 可變區(qū)的高通量測(cè)序(Illumina Miseq 2000, 上海美吉生物醫(yī)藥科技有限公司)。測(cè)序深度指數(shù)為 Coverage, 每個(gè)樣品測(cè) 40000~70000 條序列, 樣品覆蓋度均在 99.8%以上。通過(guò) Miseq 測(cè)序得到的雙端測(cè)序 reads, 按照重疊關(guān)系進(jìn)行拼接, 將測(cè)序獲得的原始序列數(shù)據(jù)進(jìn)行質(zhì)控過(guò)濾, 得到高質(zhì)量序列, 然后進(jìn)行操作分類單元(operational taxonomic units, OTU)的聚類分析, 并通過(guò) RDP classifier 貝葉斯算法, 在 97%的相似水平上對(duì) OTU 進(jìn)行物種分類學(xué)注釋(80%置信度)?;?OTU 數(shù)據(jù), 進(jìn)行物種分類學(xué)分析, 分析樣品的多樣性和菌群結(jié)構(gòu)。通過(guò)每個(gè)OTU 對(duì)應(yīng)的 greengene id, 獲得 OTU 對(duì)應(yīng)的 clusters of orthologous groups of proteins (COG)家族信息和KEGG orthology(KO)信息, 并計(jì)算各 COG 的豐度和 KO 豐度。根據(jù) COG 數(shù)據(jù)庫(kù)的信息, 可以從 Egg NOG 數(shù)據(jù)庫(kù)中解析得到各個(gè) COG 的描述信息和功能信息, 從而得到 COG 功能豐度譜。

    2 結(jié)果與討論

    2.1 反應(yīng)器的脫氮性能

    兩個(gè)反應(yīng)器進(jìn)出水的 NH4+-N, NO2?-N和NO3?-N濃度變化見圖 3(a)和(b), TN去除率及氮去除負(fù)荷(nitrogen removal rate, NRR)變化見圖 3(c)和(d)。在反應(yīng)器的啟動(dòng)階段, 兩個(gè)反應(yīng)器進(jìn)水 NH4+-N 和NO2?-N 濃度均為 50±5mg/L, HRT 設(shè)置為 24 小時(shí)。兩個(gè)反應(yīng)器出水中 NH4+-N 和 NO2?-N濃度均呈下降趨勢(shì)。啟動(dòng)階段末期(第 9 天), 對(duì)照反應(yīng)器的出水 NH4+-N 降至 11.9mg/L, 沸石反應(yīng)器的出水 NH4+-N 降至 10.2mg/L; 兩反應(yīng)器的出水NO2?-N 均降至 0mg/L。對(duì)照反應(yīng)器和沸石反應(yīng)器的 NRR 分別為18.6mgN/(L·d)和 18.12mgN/(L·d), 無(wú)明顯差別。在負(fù)荷提升階段(第 11~30 天), 厭氧氨氧化菌的脫氮活性迅速提升, TN 去除率呈直線增長(zhǎng), NRR的增長(zhǎng)速率也比啟動(dòng)階段有一定程度的提高。由圖 3(c)可以看出, 在負(fù)荷提升階段, 兩個(gè)反應(yīng)器的 TN 去除率曲線基本上重合。第 11~15 天, TN 去除率由 75%上升至 84%, 第 16~30 天, TN 去除率穩(wěn)定在 80%± 3%。兩個(gè)反應(yīng)器的 NRR 從第 11 天的 26.5mgN/ (L·d)增加到第 30 天的 73 mgN/(L·d)。

    在穩(wěn)定運(yùn)行階段(第 31~82 天), 厭氧氨氧化菌的脫氮活性開始穩(wěn)定, 此時(shí)沸石反應(yīng)器的穩(wěn)定性及脫氮性能開始高于對(duì)照反應(yīng)器。特別是當(dāng) HRT 由12 小時(shí) 縮短為 8 小時(shí)(第 59 天), 兩反應(yīng)器的 TN 去除率均呈斷崖式下降, 對(duì)照反應(yīng)器下降 42.3%, 沸石反應(yīng)器下降 39.1%。兩反應(yīng)器的 NRR 下降幅度不大, 對(duì)照反應(yīng)器下降 25.64mgN/(L·d), 沸石反應(yīng)器僅下降 17.08mgN/(L·d)。由此看出, 在反應(yīng)器內(nèi)部受到較大沖擊時(shí)(突然縮減 HRT 等), 沸石反應(yīng)器具有更好的穩(wěn)定性。TN 去除率的顯著下降可能是由于過(guò)早地縮短 HRT, 使進(jìn)水中 TN 超出反應(yīng)器的TN 去除負(fù)荷, 導(dǎo)致在出水時(shí), 反應(yīng)器內(nèi)還有大量NH4+-N 和 NO2?-N 未去除, 厭氧氨氧化菌被過(guò)高的NO2?-N 濃度抑制, 使得 TN 去除效率大幅度下滑。

    由圖 3(c)和(d)可以看出, 在穩(wěn)定運(yùn)行階段的后期(第 59~82 天), 沸石反應(yīng)器的 NRR 和 TN 去除率開始高于對(duì)照反應(yīng)器, 并逐漸與對(duì)照反應(yīng)器拉開差距。沸石反應(yīng)器的 TN 去除率由第 58 天的 47.6%增加到第 82 天的 64.1%, 增幅為 16.5%; 對(duì)照反應(yīng)器的 TN 去除率在這 25 天里由 43.9%增至 55.2%, 增幅 11.4%; 兩反應(yīng)器 TN 去除率增幅相差5.1%。沸石反應(yīng)器的 NRR 在這 25 天里由 139.3mgN/(L·d)增至 177.8mgN/(L·d), 增幅為 38.5mgN/(L·d); 對(duì)照反應(yīng)器的 NRR 在這 25 天里由 129.1mgN/(L·d)增至154.3mgN/(L·d), 增幅為 25.2mgN/(L·d)。沸石反應(yīng)器的 NRR 增幅比對(duì)照反應(yīng)器提高 52.78%, 說(shuō)明在反應(yīng)器受到較大沖擊后, 與對(duì)照反應(yīng)器相比, 沸石反應(yīng)器可以更高效地恢復(fù)反應(yīng)器的整體脫氮效率。在整個(gè) 82 天的運(yùn)行周期中, 沸石反應(yīng)器與對(duì)照反應(yīng)器的 NRR 經(jīng) t 檢驗(yàn)得到=0.6353(>0.05), 說(shuō)明兩個(gè)反應(yīng)器的整體脫氮效率無(wú)明顯差異。但是, 在反應(yīng)器運(yùn)行的最后 24 天, 即當(dāng) HRT 縮短為 8 小時(shí)的階段, 兩個(gè)反應(yīng)器 NRR 值經(jīng) t 檢驗(yàn)得出= 0.0007(<0.01), 即具有顯著的差異性。這說(shuō)明在穩(wěn)定運(yùn)行期后期, 沸石反應(yīng)器的總體脫氮效率顯著高于對(duì)照反應(yīng)器。綜上所述, 添加沸石可以提高反應(yīng)器的整體脫氮效率, 并增強(qiáng)反應(yīng)器內(nèi)部系統(tǒng)的穩(wěn)定性。

    圖3 對(duì)照反應(yīng)器(a)和沸石反應(yīng)器(b)進(jìn)出水中NH4+-N, NO2?-N和NO3?-N濃度變化、反應(yīng)器總氮(TN)去除率變化(c)以及氮去除負(fù)荷(NRR)變化(d)

    2.2 反應(yīng)器運(yùn)行過(guò)程中 EPS 多糖和蛋白濃度變化分析

    由圖 4 可知, 兩反應(yīng)器內(nèi)懸浮污泥 EPS 含量總體上隨時(shí)間增加而升高。前 3 次采樣 (第 15 天、第30 天和第 45 天)的兩反應(yīng)器內(nèi)懸浮態(tài)污泥 EPS 含量經(jīng) t 檢驗(yàn)得=0.5152(>0.05), 無(wú)顯著性差異。并且, 前 3 次采樣的沸石反應(yīng)器內(nèi)懸浮態(tài)污泥 EPS 含量均小于對(duì)照反應(yīng)器, 可能是由于前期過(guò)多的菌泥在沸石表面形成過(guò)厚的生物膜, 影響內(nèi)層厭氧氨氧化菌的恢復(fù)[17], 同時(shí)沸石會(huì)截留部分細(xì)菌, 使其無(wú)法與反應(yīng)器內(nèi)的基質(zhì)充分反應(yīng)。但是, 到穩(wěn)定運(yùn)行后期(第 82 天), 沸石反應(yīng)器內(nèi)懸浮菌樣 EPS 含量(103.2 mg/gVSS)首次超越對(duì)照反應(yīng)器(97.1mg/gVSS)。陳方敏等[18]發(fā)現(xiàn), 高活性的厭氧氨氧化細(xì)菌 EPS 濃度范圍為 100.1~105.6mg/gVSS, 本研究與之相近。綜上所述, 雖然在第 15 天、第 30 天和第 45 天, 對(duì)照組反應(yīng)器內(nèi)懸浮態(tài)污泥 EPS 含量高于沸石反應(yīng)器, 但兩者差異并不顯著。經(jīng)過(guò)長(zhǎng)時(shí)間(82 天)的運(yùn)行, 沸石反應(yīng)器內(nèi)懸浮態(tài)污泥 EPS 濃度高于對(duì)照反應(yīng)器。結(jié)合 TN 去除率來(lái)看, 在反應(yīng)器運(yùn)行的前 58天內(nèi), 兩反應(yīng)器的 TN 去除曲線幾乎重合, TN 總體去除率并無(wú)明顯的差異(=0.8201>0.05)。沸石反應(yīng)器在懸浮污泥 EPS 含量低于對(duì)照反應(yīng)器的情況下, 還保持與對(duì)照反應(yīng)器幾乎一樣的脫氮效率, 這也說(shuō)明添加沸石載體后, 沸石吸附 NH4+-N, 使得反應(yīng)器的脫氮更加高效和穩(wěn)定。

    圖4 反應(yīng)器內(nèi)懸浮污泥的蛋白濃度(a)和多糖濃度(b)

    2.3 反應(yīng)器內(nèi)微生物群落結(jié)構(gòu)分析

    為了進(jìn)一步探究沸石吸氨對(duì)反應(yīng)器內(nèi)菌群組成的影響, 將對(duì)照反應(yīng)器和沸石反應(yīng)器內(nèi)懸浮污泥及沸石反應(yīng)器內(nèi)沸石表面附著態(tài)污泥進(jìn)行 16S rRNA基因高通量測(cè)序。首先對(duì)采集樣本的細(xì)菌測(cè)序結(jié)果進(jìn)行主成分分析(principal component analysis, PCA), 探究不同時(shí)期兩個(gè)反應(yīng)器內(nèi)微生物群落結(jié)構(gòu)的相似性與差異性。PCA 分析可以對(duì)數(shù)據(jù)進(jìn)行簡(jiǎn)化, 去除干擾元素, 降維后在二維坐標(biāo)圖中表征出樣本間的相似性。樣本組成越相似, 則在 PCA 圖中的距離越近, 出現(xiàn)樣品數(shù)據(jù)點(diǎn)聚集的現(xiàn)象。

    如圖 5 所示, 選取的樣品共聚為 4 類。反應(yīng)器運(yùn)行啟動(dòng)階段, 即第一次采樣的采樣點(diǎn) B_1, B_2 和B_3 聚為一類, 說(shuō)明在該階段, 沸石反應(yīng)器內(nèi)污泥的微生物群落結(jié)構(gòu)沒(méi)有發(fā)生大的改變。隨著反應(yīng)器的運(yùn)行, 從負(fù)荷提升期開始, 沸石反應(yīng)器內(nèi)污泥的群落結(jié)構(gòu)發(fā)生顯著的改變——沸石反應(yīng)器懸浮樣C_2 和沸石反應(yīng)器附著樣 C_3 開始與對(duì)照反應(yīng)器懸浮樣 C_1 拉大距離。但是, C_2 與 C_3 的距離仍然接近, 說(shuō)明此時(shí)沸石反應(yīng)器內(nèi)懸浮態(tài)污泥與附著態(tài)污泥的群落結(jié)構(gòu)相似。到穩(wěn)定運(yùn)行期, 沸石反應(yīng)器內(nèi)懸浮樣(D_2 和 E_2)、沸石反應(yīng)器內(nèi)附著樣(D_3和 E_3)以及對(duì)照反應(yīng)器懸浮樣(D_1 和 E_1)各自形成 3 個(gè)距離較遠(yuǎn)的團(tuán)簇, 說(shuō)明在穩(wěn)定運(yùn)行期, 對(duì)照反應(yīng)器懸浮樣與沸石反應(yīng)器懸浮樣的群落結(jié)構(gòu)差異繼續(xù)變大, 而沸石反應(yīng)器內(nèi)部懸浮態(tài)污泥與附著態(tài)污泥的群落結(jié)構(gòu)也開始出現(xiàn)差異。上述結(jié)果表明, 沸石的加入對(duì)反應(yīng)器內(nèi)污泥微生物群落結(jié)構(gòu)的演替有顯著的影響。

    2.4 菌群分類及豐度變化情況

    由圖 6 和 7 可知, 對(duì)照反應(yīng)器內(nèi)懸浮污泥中厭氧氨氧化菌屬主要是其次是。隨著反應(yīng)器的運(yùn)行,在對(duì)照反應(yīng)器內(nèi)懸浮污泥中的豐度占比為B_1 6.9%, C_1 29%, D_1 30%和 E_1 34%, 呈增加的趨勢(shì), 且在負(fù)荷提升期(第 15~45 天)的幅度較大, 為 22.1%。在穩(wěn)定運(yùn)行期(第 45~82 天), 增幅較小, 僅為 5%。沸石反應(yīng)器內(nèi)懸浮菌樣中厭氧氨氧化菌屬主要是, 其初始占比為0.29%。隨著反應(yīng)器的運(yùn)行,在沸石反應(yīng)器內(nèi)懸浮菌樣中的豐度占比為 B_2 4.9%, C_2 20%, D_2 33%和 E_2 41%, 呈增加的趨勢(shì)。與對(duì)照反應(yīng)器不同的是, 沸石反應(yīng)器內(nèi)懸浮態(tài)菌在穩(wěn)定運(yùn)行時(shí)期,菌屬豐度占比的增幅(21%)較大。其中, 在第 60 天及第 90 天取得的污泥樣品中, 沸石反應(yīng)器內(nèi)懸浮菌樣中菌屬的相對(duì)豐度均超過(guò)對(duì)照反應(yīng)器, 兩反應(yīng)器內(nèi)屬豐度占比差距由 3%(第 60 天)拉大到 7% (第 90 天)。上述結(jié)果說(shuō)明沸石的加入可以提升反應(yīng)器內(nèi)懸浮污泥中厭氧氨氧化菌的相對(duì)豐度和增長(zhǎng)速率。

    屬在沸石反應(yīng)器內(nèi)附著態(tài)污泥中的豐度占比為 B_2 12%, C_2 36%, D_2 34%和 E_2 35%。附著態(tài)污泥內(nèi)菌屬含量在負(fù)荷提升期(第 15 天)及穩(wěn)定運(yùn)行期前期(第 45 天)均比同時(shí)期對(duì)照反應(yīng)器內(nèi)懸浮態(tài)污泥高。但是, 在穩(wěn)定運(yùn)行期的中、后期(第 45, 60 和 90 天), 附著態(tài)污泥內(nèi)屬的豐度占比變化不大, 維持在 35%±1%??赡苁且?yàn)榉惺砻娓街勰鄡?nèi)厭氧氨氧化菌的相對(duì)豐度在這段時(shí)間一直處于較高水平, 到達(dá)“瓶頸”。在反應(yīng)器運(yùn)行前期(第1~45 天), 附著態(tài)污泥內(nèi)屬的豐度占比增長(zhǎng)比懸浮態(tài)污泥快的原因可能是沸石良好的吸附 NH4+-N 能力以及沸石對(duì)厭氧氨氧化菌提供的載體功能。在反應(yīng)器運(yùn)行期間, 沸石反應(yīng)器內(nèi)的沸石表面形成一層厭氧氨氧化生物膜, 通過(guò)沸石自身吸附及解吸 NH4+-N 的能力, 在生物膜微環(huán)境內(nèi)進(jìn)行穩(wěn)定的厭氧氨氧化反應(yīng), 從而能夠更好地促進(jìn)受饑餓抑制厭氧氨氧化菌的恢復(fù)。

    圖5 主成分分析

    屬是一種兼性厭氧反硝化菌屬, 可在微氧條件下利用系統(tǒng)內(nèi) NO3?-N 或 NO2?-N 進(jìn)行反硝化反應(yīng)[19]。從圖 6 和 7 可以看出, 在厭氧氨氧化菌恢復(fù)過(guò)程中,屬的豐度呈逐步上升的趨勢(shì)。以對(duì)照反應(yīng)器為例,屬在懸浮態(tài)污泥中的豐度占比為 B_1 5.8%, C_1 33%, D_1 34%和 E_1 38%, 豐度占比從第 2 次采樣(第15 天)的 5.8%升至第 5 次采樣(第 45 天)的 38%, 增幅達(dá) 27.2%, 與屬在對(duì)照反應(yīng)器懸浮態(tài)污泥中豐度占比的變化情況(第 45 天比第15 天提高 22.1%)類似, 兩者均呈逐步遞增趨勢(shì)。這說(shuō)明在受饑餓影響厭氧氨氧化菌的恢復(fù)過(guò)程中,與屬于協(xié)同共生關(guān)系。在沸石反應(yīng)器內(nèi),屬在懸浮態(tài)污泥中的豐度占比為 B_2 6.3%, C_2 31%, D_2 37%和E_2 36%, 在附著態(tài)污泥中豐度占比為 B_2 6.9%, C_2 32%, D_2 39%和 E_2 36%。沸石反應(yīng)器內(nèi)懸浮態(tài)與附著態(tài)污泥中屬的豐度占比幾乎相同, 且沸石反應(yīng)器內(nèi)污泥中屬的豐度占比與對(duì)照反應(yīng)器內(nèi)污泥差別不大, 說(shuō)明沸石對(duì)反應(yīng)器內(nèi)屬的豐度影響甚微。

    圖6 對(duì)照反應(yīng)器內(nèi)懸浮態(tài)污泥樣品在屬水平上的物種分類及豐度占比

    圖7 沸石反應(yīng)器內(nèi)懸浮態(tài)污泥樣品(a)和附著態(tài)污泥樣品(b)在屬水平上的物種分類及豐度占比

    圖8 微生物群落功能的COG功能豐度變化

    2.5 功能預(yù)測(cè)分析

    16s rDNA 高通量測(cè)序結(jié)果可以用來(lái)預(yù)測(cè)微生物群落的功能。采用 PICRUSt1 功能預(yù)測(cè)法, 對(duì) 4個(gè)采樣時(shí)間點(diǎn)(第 15 天、第 45 天、第 60 天和第 82天), 兩個(gè)反應(yīng)器內(nèi)不同形態(tài)污泥中的微生物群落功能進(jìn)行 COG 功能預(yù)測(cè), 得到 COG 功能豐度箱(圖8)??梢钥闯? 3 種污泥中豐度較高的 COG 功能基因分別為細(xì)胞壁和細(xì)胞膜生物合成功能、能量產(chǎn)生和轉(zhuǎn)化功能、氨基酸轉(zhuǎn)運(yùn)和代謝功能以及信號(hào)轉(zhuǎn)換功能。其中, 豐度排第一位和第二位的功能分別是氨基酸轉(zhuǎn)運(yùn)和代謝以及能量的產(chǎn)生和轉(zhuǎn)化, 說(shuō)明反應(yīng)器中微生物的活性較高, 生長(zhǎng)速率較快, 微生物增加量大。豐度排第三位的是細(xì)胞壁和細(xì)胞膜生物合成功能, 該功能有利于微生物的附著態(tài)生長(zhǎng), 說(shuō)明厭氧氨氧化菌容易附著在沸石載體表面, 通過(guò)沸石的吸氨能力及載體功能, 促進(jìn)厭氧氨氧化反應(yīng)高效、穩(wěn)定地運(yùn)行。結(jié)合污泥中物種分類及豐度占比可以發(fā)現(xiàn), 沸石表面附著污泥內(nèi)厭氧氨氧化菌的豐度比其他兩種污泥更快地達(dá)到峰值(35%)。

    進(jìn)一步地, 通過(guò)對(duì)兩個(gè)反應(yīng)器內(nèi) 3 種不同污泥狀態(tài)的 4 個(gè)采樣時(shí)間點(diǎn)(第 15 天、第 45 天、第 60 天和第 82 天)的 COG 功能基因豐度之間進(jìn)行 t 檢驗(yàn), 尋找<0.05 的 COG 功能分類, 結(jié)果如圖 9 所示。兩種對(duì)照反應(yīng)器懸浮樣以及沸石反應(yīng)器懸浮樣和附著樣均存在顯著性差異的 COG 功能分類: 信號(hào)傳導(dǎo)通路、細(xì)胞防御機(jī)理。對(duì)于沸石表面附著的污泥, 向反應(yīng)器內(nèi)加入沸石可能會(huì)增加污泥內(nèi)細(xì)胞的防御代謝, 增強(qiáng)細(xì)菌自身的排異及適應(yīng)功能。對(duì)于懸浮態(tài)細(xì)菌, 向反應(yīng)器內(nèi)加入沸石可能會(huì)加強(qiáng)細(xì)菌之間的通信。沸石反應(yīng)器內(nèi)懸浮污泥屬的豐度比對(duì)照反應(yīng)器高 7%, 說(shuō)明沸石的加入可使細(xì)菌之間信號(hào)傳導(dǎo)更加密切, 從而促進(jìn)細(xì)菌功能的恢復(fù)。

    圖9 兩反應(yīng)器內(nèi)污泥的兩種不同COG功能豐度及顯著性差異

    3 結(jié)論

    本文通過(guò)研究投加沸石載體對(duì)受饑餓影響厭氧氨氧化菌恢復(fù)過(guò)程的影響, 得到如下主要結(jié)論。

    1)投加沸石載體可使反應(yīng)器內(nèi)厭氧氨氧化反應(yīng)更加穩(wěn)定和高效, NH4+-N 吸附和微生物降解同時(shí)發(fā)生。2)投加沸石載體可使恢復(fù)后的厭氧氨氧化菌()的豐度更高, 更好地促進(jìn)受饑餓抑制厭氧氨氧化菌的長(zhǎng)期恢復(fù)。3)沸石為載體, 運(yùn)行生物膜–活性污泥復(fù)合系統(tǒng)反應(yīng)器, 通過(guò)底物吸附+生物降解耦合手段, 在沸石表面的微環(huán)境中達(dá)到長(zhǎng)期高效穩(wěn)定的 NH4+-N 底物供給, 可以更好地促進(jìn)厭氧氨氧化反應(yīng)進(jìn)行。

    [1] Jetten M S, Wagner M, Fuerst J, et al. Microbiology and application of the anaerobic ammonium oxidation “anammox” process. CurrOpin Biotechnol, 2001, 12 (3): 283–288

    [2] Zhang L, Narita Y, Gao L, et al. Maximum specific growth rate of anammox bacteria revisited. Water Research, 2017, 11(6): 296–303

    [3] Tang C J, Zheng P, Chai L Y, et al. Characterization and quantification of anammox start-up in UASB reac-tors seeded with conventional activated sludge. Inter-national Biodeterioration & Biodegradation, 2013, 82 (8): 141–148

    [4] 侯曉幫, 操家順, 周可為. 厭氧氨氧化抑制試驗(yàn)及復(fù)活策略研究. 中國(guó)給水排水, 2017, 33(23): 55–60

    [5] 楊開亮, 占志強(qiáng), 王瑩, 等. 厭氧氨氧化菌活性的快速恢復(fù)及運(yùn)行. 上海海事大學(xué)學(xué)報(bào), 2017, 38(4): 84–87

    [6] Ban Z, Dave G. Laboratory Studies on Recovery of N and P from Human Urine through Struvite Crystalli-sation and Zeolite Adsorption. Environmental Tech-nology, 2004, 25(1): 111–121

    [7] Inglezakis V J. The concept of “capacity” in zeolite ion-exchange systems. Journal of Colloid&Interface Science, 2005, 281(1): 68–79

    [8] Miladinovic N, Weatherley L R. Intensification of ammonia removal in a combined ion-exchange and nitrification column. Chemical Engineering Journal, 2008, 135(1): 15–24

    [9] Qiu L P, Zhang S B, Wang G W, et al. Performances and nitrification properties of biological aerated filters with zeolite, ceramic particle and carbonate media. Bioresour Technol, 2010, 101(19): 7245–7251

    [10] Yapsakli K, Aktan C K, Mertoglu B. Anammox-zeolite system acting as buffer to achieve stable effluent nitrogen values. Biodegradation, 2017, 28(1): 69–79

    [11] 國(guó)家環(huán)境保護(hù)總局. 水和廢水監(jiān)測(cè)分析方法. 北京: 中國(guó)環(huán)境科學(xué)出版社, 2002

    [12] Sheng G P, Yu H Q, Li X Y. Extracellular polymeric substances (EPS) of microbial aggregates in biolo-gical wastewater treatment systems: a review. Bio-technology Advances, 2010, 28(6): 882–894

    [13] Guo J, Wang S, Lian J, et al. Rapid start-up of the anammox process: effects of five different sludge extracellular polymeric substances on the activity of anammox bacteria. Bioresource Technology, 2016, 220: 641–646

    [14] Abzac P, Bordas F, Joussein E, et al. Characterization of themineral fraction associated to extracellular polymeric substances (EPS) in anaerobic granular slu-dges, Environmental Science & Technology, 2010, 44(1): 412–418

    [15] Shu D, He Y, Yue H, et al. Metagenomic and quantitative insights into microbial communities and functional genes of nitrogen and iron cycling in twelve wastewater treatment systems. Chemical Engineering Journal, 2016, 290: 21–30

    [16] Gaudy A F. Colorimetric determination of protein and carbohydrate. Water Wastes, 1962, 7: 17–22

    [17] Ali M, Oshiki M, Rathnayake L, et al. Rapid and successful start-up of anammox process by immobili-zing the minimal quantity of biomass in PVA-SA gel beads. Water Research, 2015, 79: 147–157

    [18] 陳方敏, 顧澄偉, 胡羽婷, 等. 厭氧氨氧化污泥恢復(fù)過(guò)程中的顆粒特性. 環(huán)境科學(xué), 2018, 39(12): 5605–5611

    [19] Fahrbach M, Kuever J, Meinke R, et al.gen. nov. , sp. nov. , a 17beta-oe-stradiol-degrading, denitrifying betaproteobacte-rium. International Journal of Systematic & Evolutio-nary Microbiology, 2006, 56(7): 1547–1552

    Effect of Zeolite on Recovery of Anammox Bacteria Affected by Starvation

    YU Daodao1, SUN Jingqi2, HUO Tangran2, ZHAO Yunpeng2, LIU Sitong1,2,?

    1. Department of Environment and Energy, Peking University Shenzhen Graduate School, Shenzhen 518055; 2. Department of Environmental Science & Engineering, Peking University, Beijing 100871; ? Corresponding author, E-mail: liusitong@pku.edu.cn

    Zeolite has good ammonium absorption capacity and microbial carrier function. Anammox bacteria affected by starvation were activated in a reactor containing zeolite to study the effect of zeolite on the activity and structure of anammox bacteria after long-term starvation. The results showed that anammox bacteria recovered in the zeolite reactor and the control reactor (without zeolite addition) after 82 days of operation. The nitrogen removal efficiency of zeolite reactor began to lag behind that of the control reactor on the 59th day. On the 82nd day of operation, the nitrogen removal load (NRR) of zeolite reactor was 177.8 mgN/(L·d), significantly higher than that of the control reactor (154.3 mgN/(L·d)). At the same time, 16S rRNA high-throughput sequencing method was used to analyze the microbial composition and function of the microbial community in the reactor. It was found that adding zeolite carrier could increase the abundance of anammox bacteriaafter recovery, and could better promote the recovery of anammox bacteria inhibited by hunger.

    anammox bacteria; zeolite; bacterial recovery; community structure

    10.13209/j.0479-8023.2021.011

    2020–03–07;

    2020–03–24

    猜你喜歡
    沸石饑餓反應(yīng)器
    聚丙烯環(huán)管反應(yīng)器升密操作與控制
    云南化工(2021年11期)2022-01-12 06:06:36
    沸石分子篩發(fā)展簡(jiǎn)述
    云南化工(2021年10期)2021-12-21 07:33:24
    5種沸石分子篩的吸附脫碳對(duì)比實(shí)驗(yàn)
    煤氣與熱力(2021年9期)2021-11-06 05:22:56
    EGSB反應(yīng)器處理阿維菌素廢水
    上旋流厭氧反應(yīng)器在造紙廢水處理中的應(yīng)用
    回憶饑餓
    一張饑餓年代的教師“特供證”
    文史春秋(2016年8期)2016-02-28 17:41:34
    費(fèi)托合成微反應(yīng)器研究進(jìn)展
    沸石再生
    石油化工(2015年9期)2015-08-15 00:43:05
    饑餓的磊編
    国产精品日韩av在线免费观看| 淫秽高清视频在线观看| 大话2 男鬼变身卡| 最近中文字幕2019免费版| 午夜免费观看性视频| 精品久久久噜噜| 97精品久久久久久久久久精品| 久久精品夜色国产| 美女高潮的动态| 国产v大片淫在线免费观看| 18禁动态无遮挡网站| 小蜜桃在线观看免费完整版高清| 午夜福利网站1000一区二区三区| 有码 亚洲区| 亚洲四区av| 乱码一卡2卡4卡精品| 三级经典国产精品| 菩萨蛮人人尽说江南好唐韦庄| av在线播放精品| 韩国高清视频一区二区三区| 欧美变态另类bdsm刘玥| 亚洲真实伦在线观看| 床上黄色一级片| 久久久久九九精品影院| 午夜精品一区二区三区免费看| 国产黄片美女视频| 国产精品一及| 老司机影院成人| 最近最新中文字幕大全电影3| 人人妻人人澡人人爽人人夜夜 | 乱人视频在线观看| 麻豆乱淫一区二区| 伦精品一区二区三区| 国产精品国产三级专区第一集| 成人特级av手机在线观看| 国产精品一区www在线观看| 国产成人精品一,二区| 全区人妻精品视频| 又粗又硬又长又爽又黄的视频| 免费观看的影片在线观看| 高清在线视频一区二区三区| 国产午夜福利久久久久久| 日日撸夜夜添| 成年人午夜在线观看视频 | 又大又黄又爽视频免费| 能在线免费看毛片的网站| 丝瓜视频免费看黄片| 看免费成人av毛片| 日韩制服骚丝袜av| 亚洲欧美清纯卡通| 日韩一区二区视频免费看| 纵有疾风起免费观看全集完整版 | 3wmmmm亚洲av在线观看| 在线观看美女被高潮喷水网站| 麻豆精品久久久久久蜜桃| 小蜜桃在线观看免费完整版高清| 国产真实伦视频高清在线观看| 18禁在线播放成人免费| 午夜福利成人在线免费观看| 国产高清不卡午夜福利| 99久久中文字幕三级久久日本| 成人无遮挡网站| 亚洲精品一区蜜桃| 大陆偷拍与自拍| 久久久成人免费电影| 亚洲精品影视一区二区三区av| 成年av动漫网址| 精品酒店卫生间| 久久久久精品久久久久真实原创| 亚洲欧美清纯卡通| 80岁老熟妇乱子伦牲交| 亚洲精品成人久久久久久| 国产高清有码在线观看视频| 精品久久久久久成人av| 寂寞人妻少妇视频99o| 亚洲一级一片aⅴ在线观看| 婷婷色综合www| 日本一本二区三区精品| 在线观看av片永久免费下载| 美女xxoo啪啪120秒动态图| 成年免费大片在线观看| 精品少妇黑人巨大在线播放| 国内揄拍国产精品人妻在线| 国产精品久久久久久av不卡| av在线老鸭窝| 亚洲电影在线观看av| xxx大片免费视频| 国产精品嫩草影院av在线观看| 亚洲av成人精品一区久久| 寂寞人妻少妇视频99o| 九色成人免费人妻av| 国产精品一区www在线观看| 中文字幕av在线有码专区| 国产黄色免费在线视频| 国产免费视频播放在线视频 | 午夜激情久久久久久久| 韩国av在线不卡| 国产成人一区二区在线| 亚洲自拍偷在线| 国精品久久久久久国模美| 高清毛片免费看| 欧美区成人在线视频| 爱豆传媒免费全集在线观看| 一本一本综合久久| kizo精华| 婷婷六月久久综合丁香| 男人和女人高潮做爰伦理| 免费人成在线观看视频色| 亚洲精品视频女| 美女xxoo啪啪120秒动态图| 91精品国产九色| 少妇高潮的动态图| 精品午夜福利在线看| 久久久久精品性色| 日韩国内少妇激情av| 少妇丰满av| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产精品国产精品| av在线蜜桃| 又粗又硬又长又爽又黄的视频| 日本免费a在线| 国产v大片淫在线免费观看| 亚洲丝袜综合中文字幕| 在线观看av片永久免费下载| 亚洲aⅴ乱码一区二区在线播放| 精品人妻一区二区三区麻豆| 亚洲精品乱码久久久久久按摩| 亚洲av电影在线观看一区二区三区 | 国产精品一二三区在线看| 看免费成人av毛片| 又粗又硬又长又爽又黄的视频| 狂野欧美白嫩少妇大欣赏| 91狼人影院| 男人舔女人下体高潮全视频| 男人爽女人下面视频在线观看| 午夜免费男女啪啪视频观看| 久久精品夜夜夜夜夜久久蜜豆| 久久久久九九精品影院| 亚洲国产日韩欧美精品在线观看| 久久久成人免费电影| 男人和女人高潮做爰伦理| 国产又色又爽无遮挡免| 国产精品av视频在线免费观看| 大片免费播放器 马上看| 亚洲精品456在线播放app| 亚洲av中文字字幕乱码综合| 国产av码专区亚洲av| 日韩av在线免费看完整版不卡| 最近最新中文字幕免费大全7| 亚洲精品久久午夜乱码| av国产久精品久网站免费入址| 如何舔出高潮| 亚洲性久久影院| 麻豆久久精品国产亚洲av| 久久久久久久久久黄片| 亚洲va在线va天堂va国产| 国产精品伦人一区二区| 欧美潮喷喷水| 国产黄片美女视频| 亚洲精品自拍成人| 国产亚洲91精品色在线| 成人无遮挡网站| 又爽又黄无遮挡网站| 国产 一区 欧美 日韩| 国产成人freesex在线| 亚洲精品aⅴ在线观看| 搡女人真爽免费视频火全软件| 国产视频首页在线观看| 老司机影院成人| 黄色欧美视频在线观看| 日韩欧美 国产精品| 国产男人的电影天堂91| 高清欧美精品videossex| 精品久久久久久久久久久久久| 国产毛片a区久久久久| 国产午夜精品久久久久久一区二区三区| 亚洲熟妇中文字幕五十中出| 免费观看的影片在线观看| 久久久久久久大尺度免费视频| kizo精华| 日韩欧美 国产精品| 日韩不卡一区二区三区视频在线| 日韩成人av中文字幕在线观看| 男女那种视频在线观看| av在线老鸭窝| 日本-黄色视频高清免费观看| 又大又黄又爽视频免费| 亚洲国产精品成人综合色| 久久精品国产亚洲网站| 欧美潮喷喷水| 午夜福利高清视频| 国产亚洲91精品色在线| 国产精品久久视频播放| 成年免费大片在线观看| 国产一区有黄有色的免费视频 | 国产一区亚洲一区在线观看| 亚洲精品国产av蜜桃| 欧美三级亚洲精品| 中国国产av一级| 国产在视频线在精品| 日韩精品青青久久久久久| 天天躁日日操中文字幕| 久久久久久久久久成人| 又粗又硬又长又爽又黄的视频| 亚洲成人一二三区av| 你懂的网址亚洲精品在线观看| 综合色av麻豆| 高清在线视频一区二区三区| 午夜福利在线在线| 久久久色成人| 亚洲成色77777| 激情五月婷婷亚洲| 成人美女网站在线观看视频| 亚洲美女视频黄频| 少妇人妻精品综合一区二区| 国产免费福利视频在线观看| 精品久久久久久久末码| 久久精品国产亚洲av涩爱| 国产男人的电影天堂91| 精品亚洲乱码少妇综合久久| 亚洲欧美精品专区久久| 国产精品熟女久久久久浪| 日本熟妇午夜| 狂野欧美激情性xxxx在线观看| 亚洲av在线观看美女高潮| 亚洲电影在线观看av| av播播在线观看一区| 国产大屁股一区二区在线视频| 亚洲性久久影院| 免费高清在线观看视频在线观看| 天堂√8在线中文| 在线 av 中文字幕| 成人一区二区视频在线观看| 边亲边吃奶的免费视频| 99热网站在线观看| 3wmmmm亚洲av在线观看| 色播亚洲综合网| 精品人妻一区二区三区麻豆| 欧美xxⅹ黑人| 久久精品综合一区二区三区| 好男人在线观看高清免费视频| 中文字幕久久专区| 国产永久视频网站| 国产男人的电影天堂91| 午夜爱爱视频在线播放| 免费人成在线观看视频色| 精品一区二区三区人妻视频| 久久人人爽人人爽人人片va| 大香蕉97超碰在线| 午夜免费观看性视频| 国产三级在线视频| 国产精品久久久久久精品电影小说 | 午夜福利视频精品| 免费看光身美女| 亚洲国产欧美人成| 一级爰片在线观看| 日韩制服骚丝袜av| 99热这里只有精品一区| 久久99精品国语久久久| 熟女人妻精品中文字幕| 最近中文字幕2019免费版| 一个人免费在线观看电影| 国产欧美日韩精品一区二区| 97超视频在线观看视频| 高清毛片免费看| 亚洲av中文av极速乱| 欧美日韩视频高清一区二区三区二| 国产精品蜜桃在线观看| 成人亚洲欧美一区二区av| 国产成人免费观看mmmm| 最后的刺客免费高清国语| 最近最新中文字幕大全电影3| av天堂中文字幕网| 国产成人91sexporn| 婷婷色综合www| av在线亚洲专区| 久久人人爽人人爽人人片va| av在线老鸭窝| 我的女老师完整版在线观看| 亚洲av在线观看美女高潮| 中文欧美无线码| 久久久久精品性色| 精品国产一区二区三区久久久樱花 | 国产精品99久久久久久久久| 国国产精品蜜臀av免费| 国产淫语在线视频| 日韩伦理黄色片| 2021少妇久久久久久久久久久| 天堂av国产一区二区熟女人妻| 男人舔奶头视频| 综合色av麻豆| 欧美3d第一页| 久久鲁丝午夜福利片| 精华霜和精华液先用哪个| 午夜老司机福利剧场| 国产一区有黄有色的免费视频 | 草草在线视频免费看| 亚洲av免费高清在线观看| 80岁老熟妇乱子伦牲交| 只有这里有精品99| 国产黄色小视频在线观看| 日本三级黄在线观看| 欧美激情国产日韩精品一区| 日韩一区二区三区影片| 看十八女毛片水多多多| av在线观看视频网站免费| 国模一区二区三区四区视频| 日韩电影二区| 2022亚洲国产成人精品| 国产高清有码在线观看视频| 国产精品蜜桃在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产精品美女特级片免费视频播放器| 天堂av国产一区二区熟女人妻| 中文欧美无线码| 成人午夜高清在线视频| 老师上课跳d突然被开到最大视频| 精华霜和精华液先用哪个| 建设人人有责人人尽责人人享有的 | 国产麻豆成人av免费视频| 久久久国产一区二区| 久久精品国产亚洲av天美| 免费观看a级毛片全部| 亚洲精品成人av观看孕妇| 国产免费一级a男人的天堂| 丰满少妇做爰视频| 久久久久久久久久人人人人人人| 国产成人a∨麻豆精品| 日韩,欧美,国产一区二区三区| av一本久久久久| 日韩视频在线欧美| 国产欧美另类精品又又久久亚洲欧美| 狠狠精品人妻久久久久久综合| 91aial.com中文字幕在线观看| 少妇熟女欧美另类| 国产极品天堂在线| 久久精品久久久久久久性| 日韩人妻高清精品专区| 能在线免费观看的黄片| 午夜福利高清视频| 亚洲欧美一区二区三区国产| 亚洲精品亚洲一区二区| 秋霞在线观看毛片| 午夜福利在线在线| 免费在线观看成人毛片| 日日啪夜夜撸| 久久99热6这里只有精品| 亚洲精品自拍成人| 国产精品国产三级国产专区5o| 日韩人妻高清精品专区| 简卡轻食公司| 亚洲精品乱码久久久久久按摩| 简卡轻食公司| 男人舔奶头视频| 亚洲欧美一区二区三区国产| 午夜福利在线在线| 国产亚洲一区二区精品| 亚洲av中文av极速乱| 丰满少妇做爰视频| 成人美女网站在线观看视频| 久久久久久久大尺度免费视频| 人妻一区二区av| 美女大奶头视频| 欧美xxxx黑人xx丫x性爽| 欧美极品一区二区三区四区| 欧美区成人在线视频| av又黄又爽大尺度在线免费看| 麻豆乱淫一区二区| 精品久久久久久久久av| 九色成人免费人妻av| 美女内射精品一级片tv| 三级毛片av免费| 久久鲁丝午夜福利片| 搡女人真爽免费视频火全软件| 久久综合国产亚洲精品| 久热久热在线精品观看| 小蜜桃在线观看免费完整版高清| 1000部很黄的大片| 亚洲18禁久久av| 亚洲最大成人手机在线| 成年av动漫网址| 国产探花极品一区二区| 深夜a级毛片| 久久韩国三级中文字幕| 中文字幕人妻熟人妻熟丝袜美| 看十八女毛片水多多多| 99久国产av精品国产电影| 男人爽女人下面视频在线观看| 免费人成在线观看视频色| 免费av不卡在线播放| 精品久久久久久久久亚洲| 亚洲av电影在线观看一区二区三区 | www.av在线官网国产| 国产成人精品久久久久久| av又黄又爽大尺度在线免费看| 三级毛片av免费| 午夜免费观看性视频| 国产高潮美女av| 亚州av有码| 能在线免费看毛片的网站| 97精品久久久久久久久久精品| 国产成人午夜福利电影在线观看| 国产av在哪里看| 日韩成人av中文字幕在线观看| 成年免费大片在线观看| 久久久久久久久久成人| 亚洲欧美成人综合另类久久久| 一级毛片我不卡| 成年女人在线观看亚洲视频 | 午夜福利网站1000一区二区三区| 久久久精品94久久精品| 亚洲成人一二三区av| 久久久久国产网址| 又粗又硬又长又爽又黄的视频| 午夜精品一区二区三区免费看| 一边亲一边摸免费视频| 18+在线观看网站| av一本久久久久| 91精品一卡2卡3卡4卡| 超碰av人人做人人爽久久| 欧美人与善性xxx| 成人午夜高清在线视频| 国产高清不卡午夜福利| 亚洲最大成人手机在线| 日韩av在线免费看完整版不卡| 一级黄片播放器| 自拍偷自拍亚洲精品老妇| 免费看av在线观看网站| 男女下面进入的视频免费午夜| 国产在线一区二区三区精| 色吧在线观看| 99久久精品一区二区三区| 搡老乐熟女国产| 免费少妇av软件| 男女边摸边吃奶| 久久精品夜夜夜夜夜久久蜜豆| 22中文网久久字幕| 搡老乐熟女国产| 一级a做视频免费观看| 午夜福利在线观看吧| 国产色爽女视频免费观看| 精品国产露脸久久av麻豆 | 七月丁香在线播放| 亚洲成人久久爱视频| 亚洲无线观看免费| 黄片wwwwww| 一本久久精品| 尾随美女入室| 亚洲国产精品专区欧美| 午夜爱爱视频在线播放| 欧美激情久久久久久爽电影| 777米奇影视久久| 久久精品久久久久久久性| 欧美日韩综合久久久久久| 久久久色成人| 亚洲欧美精品自产自拍| 三级毛片av免费| 亚洲成人精品中文字幕电影| 成年av动漫网址| 亚洲一级一片aⅴ在线观看| 青春草视频在线免费观看| av一本久久久久| 亚洲在久久综合| 欧美成人a在线观看| 免费大片黄手机在线观看| av国产久精品久网站免费入址| 国产午夜福利久久久久久| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区av在线| 亚洲精品乱码久久久v下载方式| 日韩大片免费观看网站| 人人妻人人看人人澡| 亚洲av一区综合| 少妇高潮的动态图| 国产精品嫩草影院av在线观看| 亚洲精品乱码久久久久久按摩| 国产精品伦人一区二区| 国产精品一及| 看黄色毛片网站| 亚洲精品久久午夜乱码| 熟女电影av网| 成人性生交大片免费视频hd| 国产午夜精品久久久久久一区二区三区| 欧美高清性xxxxhd video| 日韩三级伦理在线观看| 中国美白少妇内射xxxbb| 日韩av在线大香蕉| 寂寞人妻少妇视频99o| 国产乱来视频区| 又大又黄又爽视频免费| 欧美97在线视频| 乱人视频在线观看| 免费av不卡在线播放| 久久久久网色| 国产真实伦视频高清在线观看| 久久久久久久久久久丰满| 中文精品一卡2卡3卡4更新| 亚洲精品一二三| 女人十人毛片免费观看3o分钟| 成年女人看的毛片在线观看| 国产精品国产三级国产专区5o| 亚洲欧美成人综合另类久久久| 亚洲在线观看片| 深爱激情五月婷婷| 欧美日韩综合久久久久久| 国产美女午夜福利| av免费观看日本| 日本一本二区三区精品| 欧美性猛交╳xxx乱大交人| h日本视频在线播放| 成年人午夜在线观看视频 | 国产精品熟女久久久久浪| 蜜臀久久99精品久久宅男| 欧美高清性xxxxhd video| 一边亲一边摸免费视频| 亚洲在久久综合| 高清欧美精品videossex| 汤姆久久久久久久影院中文字幕 | 亚洲精品,欧美精品| 久久精品国产亚洲av天美| 有码 亚洲区| 久久鲁丝午夜福利片| 一本一本综合久久| 欧美日韩在线观看h| 观看免费一级毛片| 狂野欧美激情性xxxx在线观看| 久久久久久久久久久免费av| 国产精品一区www在线观看| 日本午夜av视频| 天堂网av新在线| 欧美日韩视频高清一区二区三区二| 99视频精品全部免费 在线| 免费黄色在线免费观看| 亚洲欧洲国产日韩| 成人欧美大片| 国产免费视频播放在线视频 | 成人鲁丝片一二三区免费| 高清日韩中文字幕在线| 97超碰精品成人国产| 97热精品久久久久久| 91午夜精品亚洲一区二区三区| 内地一区二区视频在线| 国产精品日韩av在线免费观看| 大香蕉久久网| 熟女人妻精品中文字幕| 国产精品av视频在线免费观看| 99久久精品一区二区三区| 一个人观看的视频www高清免费观看| 国产免费一级a男人的天堂| 特大巨黑吊av在线直播| 少妇的逼好多水| 亚洲国产日韩欧美精品在线观看| av卡一久久| 久久热精品热| av线在线观看网站| 久久久精品欧美日韩精品| 大陆偷拍与自拍| 韩国高清视频一区二区三区| 在线免费观看的www视频| 国产伦在线观看视频一区| 春色校园在线视频观看| 亚洲熟女精品中文字幕| 午夜精品一区二区三区免费看| 免费观看性生交大片5| 国产女主播在线喷水免费视频网站 | .国产精品久久| 非洲黑人性xxxx精品又粗又长| www.色视频.com| 人人妻人人看人人澡| 欧美日韩国产mv在线观看视频 | 纵有疾风起免费观看全集完整版 | 国产精品福利在线免费观看| 男女那种视频在线观看| 日日干狠狠操夜夜爽| 午夜福利网站1000一区二区三区| 精品久久久久久成人av| 久久久精品94久久精品| 美女脱内裤让男人舔精品视频| 久久99热6这里只有精品| 黑人高潮一二区| 少妇熟女aⅴ在线视频| 国产成人aa在线观看| 日韩制服骚丝袜av| 尾随美女入室| 在线 av 中文字幕| 日韩精品有码人妻一区| 人人妻人人看人人澡| 亚洲va在线va天堂va国产| 亚洲国产精品sss在线观看| 国产成人a区在线观看| 久久精品久久精品一区二区三区| 91久久精品电影网| 久久精品国产鲁丝片午夜精品| 国产精品日韩av在线免费观看| 精品国产一区二区三区久久久樱花 | 日韩精品有码人妻一区| 久久久久久久久中文| 高清欧美精品videossex| 97人妻精品一区二区三区麻豆| 午夜免费观看性视频| 丰满少妇做爰视频| 亚洲精品日韩av片在线观看| 午夜免费观看性视频| 欧美日韩国产mv在线观看视频 | 日韩精品青青久久久久久| 天堂影院成人在线观看| 亚洲av一区综合| 午夜福利高清视频| 极品少妇高潮喷水抽搐| 久久99热这里只频精品6学生| 中文天堂在线官网| 午夜久久久久精精品| 两个人视频免费观看高清| 亚洲精品乱码久久久久久按摩| 国产高清三级在线|