• <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
    饑餓的磊編
    国产白丝娇喘喷水9色精品| 女的被弄到高潮叫床怎么办| 韩国精品一区二区三区| 中文字幕人妻丝袜一区二区 | 日韩一本色道免费dvd| 日韩一卡2卡3卡4卡2021年| 校园人妻丝袜中文字幕| 飞空精品影院首页| 亚洲精品美女久久久久99蜜臀 | 在线天堂最新版资源| 免费播放大片免费观看视频在线观看| 久久久久人妻精品一区果冻| 国产亚洲一区二区精品| 日韩视频在线欧美| 成人毛片60女人毛片免费| 色吧在线观看| 91精品国产国语对白视频| 国产精品久久久久久精品古装| 精品一品国产午夜福利视频| 亚洲精品视频女| 国产成人精品婷婷| 国产精品免费视频内射| 午夜福利一区二区在线看| 人体艺术视频欧美日本| 成人国语在线视频| 五月天丁香电影| 国产一区有黄有色的免费视频| 久久精品人人爽人人爽视色| 99热网站在线观看| 99久国产av精品国产电影| av在线app专区| 国产精品.久久久| 婷婷色综合大香蕉| 高清不卡的av网站| videos熟女内射| 久久女婷五月综合色啪小说| 中国三级夫妇交换| a级毛片黄视频| 免费少妇av软件| 免费观看av网站的网址| 精品一区二区三区四区五区乱码 | 亚洲精品美女久久久久99蜜臀 | 精品国产露脸久久av麻豆| 久久毛片免费看一区二区三区| 国产成人午夜福利电影在线观看| 人人妻人人爽人人添夜夜欢视频| 国产精品蜜桃在线观看| av在线老鸭窝| 精品亚洲成a人片在线观看| 成人午夜精彩视频在线观看| 丝袜美腿诱惑在线| 亚洲色图 男人天堂 中文字幕| 日韩三级伦理在线观看| 卡戴珊不雅视频在线播放| 美女国产高潮福利片在线看| 国产一级毛片在线| 国产精品蜜桃在线观看| 亚洲,欧美精品.| 午夜福利网站1000一区二区三区| 久久久国产一区二区| 免费观看性生交大片5| 校园人妻丝袜中文字幕| 国产精品99久久99久久久不卡 | 欧美 亚洲 国产 日韩一| 自线自在国产av| 99国产精品免费福利视频| 男女边摸边吃奶| 成人国语在线视频| 毛片一级片免费看久久久久| 亚洲精品在线美女| 日韩视频在线欧美| 黄色一级大片看看| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美在线一区| av国产精品久久久久影院| 亚洲欧洲精品一区二区精品久久久 | 美女高潮到喷水免费观看| 精品视频人人做人人爽| 精品视频人人做人人爽| 亚洲成色77777| 亚洲av国产av综合av卡| 国产熟女午夜一区二区三区| 国产一区亚洲一区在线观看| 中国国产av一级| 男男h啪啪无遮挡| 最近2019中文字幕mv第一页| 老女人水多毛片| 国产又色又爽无遮挡免| 久久精品人人爽人人爽视色| www.av在线官网国产| 欧美精品一区二区大全| 亚洲欧洲国产日韩| 亚洲综合色惰| 最近中文字幕高清免费大全6| 美女国产高潮福利片在线看| 老女人水多毛片| 成年人免费黄色播放视频| 超色免费av| av网站在线播放免费| 亚洲情色 制服丝袜| 日韩免费高清中文字幕av| 免费黄色在线免费观看| 精品视频人人做人人爽| 99国产综合亚洲精品| 只有这里有精品99| 亚洲美女视频黄频| 只有这里有精品99| 黄频高清免费视频| 亚洲国产精品999| 精品亚洲成国产av| 性色av一级| 亚洲精品国产av成人精品| 亚洲情色 制服丝袜| 国产精品人妻久久久影院| 日本av免费视频播放| 亚洲欧美一区二区三区黑人 | 26uuu在线亚洲综合色| 天天操日日干夜夜撸| 90打野战视频偷拍视频| 国产成人午夜福利电影在线观看| 乱人伦中国视频| 国产日韩欧美视频二区| 超碰成人久久| 国精品久久久久久国模美| 97在线视频观看| 桃花免费在线播放| 日本av免费视频播放| 男的添女的下面高潮视频| 亚洲成色77777| 欧美激情 高清一区二区三区| 欧美av亚洲av综合av国产av | 少妇人妻精品综合一区二区| 亚洲欧美一区二区三区久久| 国产精品一区二区在线不卡| 可以免费在线观看a视频的电影网站 | 国产成人精品久久久久久| 成人黄色视频免费在线看| 搡老乐熟女国产| 午夜激情av网站| 天天躁狠狠躁夜夜躁狠狠躁| 赤兔流量卡办理| 国产乱人偷精品视频| 精品国产超薄肉色丝袜足j| 99久国产av精品国产电影| 久久 成人 亚洲| 久久热在线av| 有码 亚洲区| 七月丁香在线播放| 欧美日韩精品成人综合77777| 日韩av在线免费看完整版不卡| 国产有黄有色有爽视频| 久久久久国产网址| 日韩在线高清观看一区二区三区| 国产精品麻豆人妻色哟哟久久| 美女xxoo啪啪120秒动态图| 午夜av观看不卡| 美女大奶头黄色视频| 国产有黄有色有爽视频| 黄色一级大片看看| 人妻一区二区av| 一区二区日韩欧美中文字幕| 少妇精品久久久久久久| 亚洲av在线观看美女高潮| 婷婷成人精品国产| 久久人人爽人人片av| a级毛片黄视频| 久久国内精品自在自线图片| 毛片一级片免费看久久久久| 日韩成人av中文字幕在线观看| 男女边吃奶边做爰视频| 亚洲一级一片aⅴ在线观看| 亚洲av成人精品一二三区| 日本vs欧美在线观看视频| 国产免费视频播放在线视频| 国产一区二区 视频在线| 日本猛色少妇xxxxx猛交久久| 日韩不卡一区二区三区视频在线| 国产精品无大码| 国产精品国产三级专区第一集| 中文天堂在线官网| 亚洲国产精品一区二区三区在线| 免费高清在线观看日韩| 哪个播放器可以免费观看大片| 2018国产大陆天天弄谢| 色视频在线一区二区三区| 久久精品久久久久久久性| 久久毛片免费看一区二区三区| 国产av精品麻豆| 久久久久久久久久久免费av| 亚洲国产日韩一区二区| 亚洲成av片中文字幕在线观看 | 丝袜美腿诱惑在线| 午夜福利在线免费观看网站| 久久精品国产亚洲av涩爱| 亚洲第一青青草原| 亚洲,欧美精品.| 熟女av电影| 熟女av电影| 尾随美女入室| 免费在线观看完整版高清| 亚洲精品一二三| 午夜免费男女啪啪视频观看| av视频免费观看在线观看| 午夜福利乱码中文字幕| 在线观看免费视频网站a站| 男女无遮挡免费网站观看| 日韩精品免费视频一区二区三区| 亚洲三级黄色毛片| 日韩精品免费视频一区二区三区| 美女午夜性视频免费| 亚洲三级黄色毛片| 一区福利在线观看| 考比视频在线观看| 蜜桃在线观看..| 亚洲第一青青草原| 欧美+日韩+精品| 日本黄色日本黄色录像| 黄色一级大片看看| 色94色欧美一区二区| 免费高清在线观看日韩| 精品国产一区二区三区久久久樱花| 欧美日韩综合久久久久久| 99国产综合亚洲精品| 日本-黄色视频高清免费观看| 制服丝袜香蕉在线| 亚洲精品aⅴ在线观看| 久久久久久久久久久久大奶| 男女啪啪激烈高潮av片| 这个男人来自地球电影免费观看 | 久久久久久久久久久久大奶| 麻豆av在线久日| 一区二区三区精品91| 国产精品一国产av| 视频区图区小说| 国产97色在线日韩免费| 国产成人精品婷婷| 国产亚洲最大av| 人人澡人人妻人| 黄片播放在线免费| 在线观看免费高清a一片| 亚洲国产最新在线播放| 高清欧美精品videossex| 欧美日韩av久久| 青春草国产在线视频| 久热久热在线精品观看| 狠狠精品人妻久久久久久综合| 亚洲精品一区蜜桃| 国产成人精品一,二区| 午夜福利一区二区在线看| 亚洲,欧美,日韩| 国产免费现黄频在线看| 免费观看av网站的网址| 亚洲激情五月婷婷啪啪| 两性夫妻黄色片| 国产人伦9x9x在线观看 | 久久精品aⅴ一区二区三区四区 | 久久久久久久久久久久大奶| 在线观看免费日韩欧美大片| 人人妻人人澡人人看| 日韩中文字幕视频在线看片| 国产亚洲最大av| 国产成人免费无遮挡视频| 高清欧美精品videossex| 最近2019中文字幕mv第一页| 超碰97精品在线观看| 日日啪夜夜爽| 精品人妻在线不人妻| 久久精品国产亚洲av涩爱| 色网站视频免费| 电影成人av| av福利片在线| 爱豆传媒免费全集在线观看| 国产亚洲最大av| 黄色毛片三级朝国网站| 在现免费观看毛片| 精品国产超薄肉色丝袜足j| 亚洲少妇的诱惑av| av视频免费观看在线观看| 欧美激情极品国产一区二区三区| 一级毛片黄色毛片免费观看视频| 1024视频免费在线观看| 少妇人妻 视频| av不卡在线播放| 男女无遮挡免费网站观看| 母亲3免费完整高清在线观看 | 国产精品久久久久成人av| 亚洲精品久久成人aⅴ小说| 男女免费视频国产| 国产麻豆69| 五月天丁香电影| 日韩在线高清观看一区二区三区| 2022亚洲国产成人精品| 国产成人免费观看mmmm| 色哟哟·www| 一区二区三区四区激情视频| 热99国产精品久久久久久7| 另类精品久久| 欧美精品高潮呻吟av久久| 亚洲男人天堂网一区| 久久人妻熟女aⅴ| 亚洲色图 男人天堂 中文字幕| 少妇的丰满在线观看| 观看美女的网站| 日本欧美国产在线视频| 美女视频免费永久观看网站| 久久国产精品男人的天堂亚洲| 亚洲精品美女久久久久99蜜臀 | 国产免费视频播放在线视频| 亚洲欧美一区二区三区国产| 老熟女久久久| 亚洲第一青青草原| 亚洲第一区二区三区不卡| av免费观看日本| 亚洲精品国产一区二区精华液| 亚洲精品av麻豆狂野| 国产精品免费视频内射| 十八禁高潮呻吟视频| 国产伦理片在线播放av一区| 日本-黄色视频高清免费观看| 香蕉国产在线看| 日本猛色少妇xxxxx猛交久久| 久久韩国三级中文字幕| 精品少妇一区二区三区视频日本电影 | 免费高清在线观看视频在线观看| 成人国产av品久久久| 久久久久久久国产电影| 考比视频在线观看| 精品国产露脸久久av麻豆| 日本色播在线视频| 色哟哟·www| 美国免费a级毛片| 晚上一个人看的免费电影| 秋霞伦理黄片| 国产成人精品福利久久| 2022亚洲国产成人精品| 精品久久蜜臀av无| 国产精品久久久久久av不卡| 激情视频va一区二区三区| 亚洲国产看品久久| 国产高清不卡午夜福利| 成人黄色视频免费在线看| av一本久久久久| 在线 av 中文字幕| 久久精品国产亚洲av天美| 人人妻人人爽人人添夜夜欢视频| 天天操日日干夜夜撸| 国产成人免费观看mmmm| 亚洲精品av麻豆狂野| 精品少妇一区二区三区视频日本电影 | 国产在线免费精品| 汤姆久久久久久久影院中文字幕| 亚洲精品自拍成人| 一级毛片 在线播放| 久久韩国三级中文字幕| 国产精品久久久久久精品古装| 午夜福利一区二区在线看| 大话2 男鬼变身卡| 久久久久精品久久久久真实原创| 国产精品女同一区二区软件| 亚洲欧美精品综合一区二区三区 | 黄色视频在线播放观看不卡| 免费黄网站久久成人精品| 99国产精品免费福利视频| 在线观看免费日韩欧美大片| √禁漫天堂资源中文www| 国产老妇伦熟女老妇高清| 丝袜喷水一区| 99精国产麻豆久久婷婷| 国产欧美日韩综合在线一区二区| xxxhd国产人妻xxx| 国产亚洲一区二区精品| 国产不卡av网站在线观看| 中文天堂在线官网| 中文精品一卡2卡3卡4更新| 性色avwww在线观看| av线在线观看网站| 久久久久久久久久久久大奶| av卡一久久| kizo精华| 亚洲av男天堂| 啦啦啦啦在线视频资源| 亚洲成色77777| 国产欧美日韩综合在线一区二区| 国产成人av激情在线播放| 日韩一区二区三区影片| 少妇被粗大的猛进出69影院| 精品福利永久在线观看| 大片电影免费在线观看免费| 2021少妇久久久久久久久久久| 国产日韩欧美亚洲二区| 免费高清在线观看日韩| 精品一区二区三区四区五区乱码 | 久久久精品区二区三区| 久久av网站| 中文欧美无线码| 国产成人av激情在线播放| 国产又爽黄色视频| av福利片在线| 久久精品夜色国产| 色94色欧美一区二区| 亚洲国产精品一区三区| 午夜影院在线不卡| 国产精品无大码| 午夜老司机福利剧场| 男女午夜视频在线观看| 欧美激情高清一区二区三区 | 久久精品国产自在天天线| 国产av国产精品国产| 在线亚洲精品国产二区图片欧美| 久久久久国产网址| 亚洲欧洲国产日韩| av国产久精品久网站免费入址| 女人高潮潮喷娇喘18禁视频| 国产有黄有色有爽视频| 亚洲国产欧美在线一区| 欧美日韩av久久| 亚洲人成网站在线观看播放| 午夜日本视频在线| 久久人人爽av亚洲精品天堂| 欧美亚洲日本最大视频资源| 国产欧美亚洲国产| 一区二区av电影网| 性色avwww在线观看| 免费少妇av软件| 老鸭窝网址在线观看| 中文字幕最新亚洲高清| 成人二区视频| 国产一区二区激情短视频 | 最新中文字幕久久久久| 色播在线永久视频| 亚洲国产毛片av蜜桃av| 欧美精品一区二区大全| 久久午夜综合久久蜜桃| 丝袜在线中文字幕| 夫妻性生交免费视频一级片| 亚洲精品乱久久久久久| 天天躁夜夜躁狠狠躁躁| 亚洲人成电影观看| 欧美精品亚洲一区二区| 看十八女毛片水多多多| 一区二区三区四区激情视频| 深夜精品福利| 国产视频首页在线观看| 一级爰片在线观看| 欧美日韩视频精品一区| a 毛片基地| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 欧美 日韩 精品 国产| 成人毛片60女人毛片免费| 夫妻午夜视频| 免费黄频网站在线观看国产| 视频在线观看一区二区三区| 久久精品国产综合久久久| 久久这里有精品视频免费| 1024香蕉在线观看| 丝袜美腿诱惑在线| 夫妻午夜视频| 侵犯人妻中文字幕一二三四区| 久久精品国产亚洲av涩爱| 亚洲综合色惰| 欧美精品亚洲一区二区| 亚洲在久久综合| 女人精品久久久久毛片| 777久久人妻少妇嫩草av网站| 国产精品二区激情视频| 精品一区二区免费观看| 国产 精品1| 亚洲国产日韩一区二区| 少妇的丰满在线观看| 国产日韩欧美视频二区| 成人亚洲精品一区在线观看| 啦啦啦在线观看免费高清www| 国产亚洲一区二区精品| 亚洲人成77777在线视频| 色哟哟·www| 久久久久久人人人人人| 天堂中文最新版在线下载| 久久影院123| 极品少妇高潮喷水抽搐| 亚洲一区二区三区欧美精品| 国产精品久久久久久av不卡| 国产精品久久久久久精品电影小说| 熟女电影av网| 9191精品国产免费久久| 久久女婷五月综合色啪小说| xxxhd国产人妻xxx| 天天影视国产精品| 波多野结衣一区麻豆| 2021少妇久久久久久久久久久| 国产毛片在线视频| 成人国语在线视频| 啦啦啦啦在线视频资源| 你懂的网址亚洲精品在线观看| 成人18禁高潮啪啪吃奶动态图| 久久精品国产鲁丝片午夜精品| 天天躁夜夜躁狠狠久久av| 中文字幕av电影在线播放| 久久精品熟女亚洲av麻豆精品| 一区在线观看完整版| 久久久久网色| 欧美激情极品国产一区二区三区| 中文精品一卡2卡3卡4更新| 飞空精品影院首页| 成人毛片60女人毛片免费| 黄色毛片三级朝国网站| 久久精品aⅴ一区二区三区四区 | 日韩伦理黄色片| 777久久人妻少妇嫩草av网站| 性高湖久久久久久久久免费观看| √禁漫天堂资源中文www| 99久久人妻综合| 日产精品乱码卡一卡2卡三| 性色avwww在线观看| 欧美+日韩+精品| 熟女少妇亚洲综合色aaa.| 熟妇人妻不卡中文字幕| 亚洲av中文av极速乱| 婷婷色综合www| 亚洲国产日韩一区二区| 色哟哟·www| videos熟女内射| 69精品国产乱码久久久| 日韩 亚洲 欧美在线| 丝袜美腿诱惑在线| 久久99热这里只频精品6学生| 在线观看免费视频网站a站| 色婷婷av一区二区三区视频| 免费黄网站久久成人精品| 国产麻豆69| 涩涩av久久男人的天堂| 十分钟在线观看高清视频www| 韩国精品一区二区三区| 精品国产一区二区久久| 亚洲国产精品国产精品| 午夜av观看不卡| 999精品在线视频| 制服诱惑二区| 少妇人妻精品综合一区二区| 在线天堂最新版资源| 久久99热这里只频精品6学生| 赤兔流量卡办理| 爱豆传媒免费全集在线观看| 天天躁夜夜躁狠狠躁躁| 久久精品亚洲av国产电影网| www日本在线高清视频| 老司机影院毛片| 啦啦啦中文免费视频观看日本| 少妇的逼水好多| 欧美精品国产亚洲| 新久久久久国产一级毛片| 观看av在线不卡| 99热全是精品| 一级片免费观看大全| 欧美人与性动交α欧美软件| 国产免费现黄频在线看| 人人妻人人添人人爽欧美一区卜| 少妇被粗大的猛进出69影院| 爱豆传媒免费全集在线观看| 亚洲av在线观看美女高潮| av福利片在线| 超碰97精品在线观看| 2022亚洲国产成人精品| 黄色视频在线播放观看不卡| 成年美女黄网站色视频大全免费| av卡一久久| 精品酒店卫生间| 青青草视频在线视频观看| 国产成人午夜福利电影在线观看| 国产精品免费大片| 国产日韩欧美在线精品| 国产97色在线日韩免费| 久久精品国产亚洲av天美| 久久精品国产自在天天线| 国产一级毛片在线| 亚洲av中文av极速乱| 欧美日韩av久久| 男女下面插进去视频免费观看| 免费在线观看视频国产中文字幕亚洲 | 亚洲欧美一区二区三区国产| 纵有疾风起免费观看全集完整版| 99久国产av精品国产电影| 中文字幕另类日韩欧美亚洲嫩草| 久久午夜福利片| 国产免费现黄频在线看| 少妇猛男粗大的猛烈进出视频| 日本黄色日本黄色录像| 精品少妇久久久久久888优播| 人人妻人人爽人人添夜夜欢视频| 青春草亚洲视频在线观看| 1024香蕉在线观看| 制服人妻中文乱码| a级片在线免费高清观看视频| 一级毛片电影观看| 国产精品国产av在线观看| 男女高潮啪啪啪动态图| 午夜免费男女啪啪视频观看| 精品一区在线观看国产| 亚洲人成网站在线观看播放| 一级毛片电影观看| 国产成人午夜福利电影在线观看| 精品一区在线观看国产| 18+在线观看网站| 国产视频首页在线观看| 1024香蕉在线观看| 9色porny在线观看| 中文字幕人妻丝袜一区二区 | 精品一区在线观看国产| kizo精华| 青青草视频在线视频观看| 91精品三级在线观看| 亚洲 欧美一区二区三区| 黄色配什么色好看| 热re99久久国产66热| 国产精品av久久久久免费|