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

    Na+和K+共存對(duì)A2/O工藝脫氮除磷效果及污泥性質(zhì)的影響

    2019-07-23 02:18:14張?zhí)m河陳子成郭靜波賈艷萍
    關(guān)鍵詞:氧區(qū)磷率鹽度

    張?zhí)m河,鄭 晶,田 蕊,陳子成,郭靜波,賈艷萍,李 正,張 健

    ?

    Na+和K+共存對(duì)A2/O工藝脫氮除磷效果及污泥性質(zhì)的影響

    張?zhí)m河1,2,鄭 晶1,田 蕊1,陳子成1※,郭靜波3,賈艷萍1,李 正1,張 健1

    (1. 東北電力大學(xué)化學(xué)工程學(xué)院,吉林 132012; 2.吉林建筑大學(xué)松遼流域水環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,長春 130118;3. 東北電力大學(xué)建筑工程學(xué)院,吉林 132012)

    為了揭示多種金屬離子共存的含鹽廢水生物處理系統(tǒng)污染物的去除機(jī)制和污泥特性,考察Na+、K+共存對(duì)A2/O工藝污染物去除率、污泥性質(zhì)和微生物群落的影響,采用高通量測(cè)序技術(shù)分析了厭氧區(qū)、缺氧區(qū)和好氧區(qū)的微生物群落結(jié)構(gòu),結(jié)合脫氮除磷效果和污泥性質(zhì)的變化,探討不同Na+/K+摩爾比下A2/O工藝優(yōu)勢(shì)種群的演替規(guī)律,以期從微生物角度明確Na+、K+共存對(duì)含鹽廢水污染物去除率的影響。結(jié)果表明:當(dāng)進(jìn)水Na+/K+摩爾比分別為2、1和0.5時(shí),A2/O工藝的COD去除率分別為80%、84%和86%,TN去除率分別為73%、77%和80%,K+濃度的提高緩解了Na+對(duì)COD和TN去除率的抑制作用;厭氧區(qū)釋磷率分別為70%、73%和74%,缺氧區(qū)吸磷率分別為53%、55%和58%,好氧區(qū)吸磷率分別為70%、72%和75%。隨著進(jìn)水Na+/K+摩爾比的降低,厭氧區(qū)、缺氧區(qū)和好氧區(qū)微生物群落的豐富度和多樣性降低,微生物群落差異顯著,變形菌門的相對(duì)豐度均升高約30%,擬桿菌門和綠彎菌門相對(duì)豐度逐漸降低。陶氏菌屬和固氮弧菌屬作為優(yōu)勢(shì)菌屬,其相對(duì)豐度逐漸增大,有利于氮磷污染物的去除。通過增加K+的濃度有利于提高氮、磷去除率,增強(qiáng)污泥的生物絮凝性和反硝化聚磷菌的活性。

    金屬離子;廢水;污泥;生物反應(yīng)器;微生物群落結(jié)構(gòu);鹽度

    0 引 言

    食品加工、制革、紡織、水產(chǎn)養(yǎng)殖和石油生產(chǎn)等過程中均會(huì)產(chǎn)生大量的含鹽廢水[1]。含鹽廢水進(jìn)入污水處理廠與活性污泥接觸后,影響活性污泥的沉降性、微生物活性和群落結(jié)構(gòu),導(dǎo)致污染物去除率發(fā)生變化。Na+是污水中最常見的低價(jià)鹽類離子,高Na+鹽度能夠大大降低微生物的豐富度和多樣性,導(dǎo)致酸類、酯類等有機(jī)污染物去除率和活性污泥的沉降性能下降,低Na+鹽度下的微生物群落與高Na+鹽度下的微生物群落差異性顯著[2-4]。低Na+鹽度對(duì)A2/O工藝缺氧區(qū)污泥的絮凝性產(chǎn)生促進(jìn)作用,亞硝酸氧化菌(nitrite-oxidizing bacteria,NOB)對(duì)Na+鹽度的耐受能力低于氨氧化菌(ammonia-oxidizing bacteria,AOB);無論是Na+鹽度升高或降低時(shí),亞硝酸鹽積累率均升高[5-6],Na+鹽度影響自養(yǎng)菌的硝化作用[2]。采用膜生物反應(yīng)器處理水產(chǎn)養(yǎng)殖廢水時(shí)發(fā)現(xiàn),Na+鹽度逐漸增大會(huì)使膜生物反應(yīng)器的處理能力變差,導(dǎo)致短程硝化反硝化現(xiàn)象的發(fā)生,但是有利于提高耐鹽微生物的豐度[7]。

    一些學(xué)者也研究了K+對(duì)活性污泥的影響。在廢水中溫厭氧消化時(shí),K+能減緩Na+鹽度對(duì)厭氧菌的抑制作用,有效提高厭氧消化的甲烷產(chǎn)量和化學(xué)需氧量(chemical oxygen demand,COD)去除率,保持厭氧微生物的形態(tài),提高厭氧微生物的活性[8]。當(dāng)Na+和K+鹽度均為30 g/L時(shí),升流式厭氧污泥床反應(yīng)器中COD降解受到抑制[9-10],K+和Na+降低了顆粒污泥和絮狀污泥的產(chǎn)甲烷活性,Na+的鹽抑制常數(shù)高于K+。低濃度的K+不會(huì)抑制酶活性,有利于污染物去除率的提高[11]。在SBR反應(yīng)器中,當(dāng)K+添加量為40 mg/L時(shí),可以促進(jìn)廢水中COD的去除;當(dāng)K+添加量為100 mg/L時(shí),可以促進(jìn)廢水中氨氮的去除[12]。這些研究主要分析了單一運(yùn)行條件下Na+、K+對(duì)好氧微生物或厭氧微生物活性和污染物去除率的影響,但是關(guān)于不同Na+/K+摩爾比、不同運(yùn)行條件下(厭氧、缺氧和好氧)的微生物群落如何變化,以及Na+、K+共存條件下影響有機(jī)物去除和脫氮作用的主要菌群及其分布規(guī)律尚不明確。

    A2/O工藝由厭氧區(qū)、缺氧區(qū)和好氧區(qū)組成,通過溶解氧的不斷變化可實(shí)現(xiàn)同步脫氮除磷[13]。本研究分別提取不同Na+/K+摩爾比條件下A2/O工藝厭氧、缺氧和好氧區(qū)的活性污泥總DNA,并對(duì)其16S rDNA V4區(qū)進(jìn)行了PCR擴(kuò)增和擴(kuò)增產(chǎn)物的高通量測(cè)序分析,研究不同Na+/K+摩爾比條件下A2/O工藝厭氧區(qū)、缺氧區(qū)和好氧區(qū)微生物群落的變化,分析鹽類金屬離子對(duì)COD去除和脫氮除磷效率的影響,探索含鹽廢水生物脫氮除磷機(jī)理。

    1 材料與方法

    1.1 試驗(yàn)裝置

    活性污泥取自吉林市污水處理廠二沉池。A2/O反應(yīng)器(圖1)采用有機(jī)玻璃制成,有效容積為48 L,其中厭氧區(qū)和缺氧區(qū)有效容積均為12 L,好氧區(qū)有效容積為24 L。厭氧區(qū)和缺氧區(qū)配有攪拌器,好氧區(qū)底部安裝曝氣頭,采用折流方式運(yùn)行。裝置內(nèi)部由隔板分為6部分,其中1個(gè)厭氧區(qū),2個(gè)缺氧區(qū),4個(gè)好氧區(qū),好氧區(qū)的部分混合液回流至缺氧區(qū)。采用豎流式二沉池,容積為5 L。A2/O反應(yīng)器連續(xù)運(yùn)行,污泥回流比為50%,硝化液回流比為100%;混合液懸浮固體(mixed liquor suspended solids,MLSS)為3 000~4 000 mg/L,水力停留時(shí)間(hydraulic retention time,HRT)為8 h,溫度為17~25 ℃。污泥馴化過程中按C:N:P=100:5:1的比例加入營養(yǎng)物質(zhì),以滿足微生物的生長需求。當(dāng)污染物達(dá)到穩(wěn)定去除效果后,采用逐步增加鹽度負(fù)荷的方法進(jìn)行耐鹽污泥的培養(yǎng),以降低鹽度沖擊對(duì)微生物造成的毒害。

    圖1 A2/O工藝流程示意圖

    試驗(yàn)用水采用模擬城市生活污水,實(shí)際城市生活污水COD、TN和TP一般為250~800、20~70和3~12 mg/L,故配水的主要成分為:無水乙酸鈉(1.2 g/L)、氯化銨(0.25 g/L)、磷酸二氫鉀(0.05 g/L)、硫酸鎂(0.05 g/L)和氯化鈣(0.01 g/L)。培養(yǎng)微生物所需的微量元素,主要成分為:氯化鈷(0.42 mg/L)、鉬酸鈉(0.15 mg/L)、硫酸錳(0.13 mg/L)、氯化鐵(0.37 mg/L)和硫酸銅(0.1 mg/L)。本研究進(jìn)水NaCl濃度保持穩(wěn)定,Na+濃度為0.2 mol/L;K+濃度采用KCl調(diào)控,K+濃度分別為0.1、0.2和0.4 mol/L,不同K+濃度下的進(jìn)水pH值均為7.5~8.0。在Na+/K+摩爾比分別為2、1和0.5條件下,考察Na+/K+對(duì)A2/O工藝反硝化除磷效果和生物絮凝性的影響。

    1.2 微生物群落結(jié)構(gòu)測(cè)定與分析

    1.2.1 樣品采集

    反應(yīng)器達(dá)到穩(wěn)定運(yùn)行后,對(duì)污泥進(jìn)行取樣分析。當(dāng)進(jìn)水Na+/K+為2、1和0.5時(shí),YM2、YM1和YM0.5分別代表厭氧區(qū)的污泥樣品,QM2、QM1和QM0.5分別代表缺氧區(qū)的污泥樣品,HM2、HM1和HM0.5分別代表好氧區(qū)的污泥樣品。

    1.2.2 樣品測(cè)定

    采用CTAB法提取活性污泥樣品中的DNA,利用瓊脂糖凝膠電泳檢測(cè)DNA純度和濃度。以稀釋后的樣品基因組DNA作為模板,采用515F-806R對(duì)16S rDNA的V4區(qū)進(jìn)行PCR擴(kuò)增,擴(kuò)增產(chǎn)物采用Illumina Miscq平臺(tái)進(jìn)行測(cè)序。所用酶和緩沖液均采用New England Biolabs公司生產(chǎn)的Phusion? High-Fidelity PCR Master Mix with GC Buffer;反應(yīng)條件:98 ℃預(yù)變性1 min,然后30個(gè)循環(huán)(98 ℃變性10 s,50 ℃退火30 s,72 ℃延伸30 s),最后72 ℃延伸5 min。利用New England Biolabs公司生產(chǎn)的NEB Next? UltraTMDNA Library Prep Kit for Illumina試劑盒進(jìn)行文庫構(gòu)建,再經(jīng)過Qubit定量和文庫檢測(cè)合格后,利用HiSeq進(jìn)行上機(jī)測(cè)序,分析污泥樣品的微生物群落結(jié)構(gòu)。

    1.2.3 數(shù)據(jù)分析

    利用Uparse (V7.0.1001)軟件進(jìn)行聚類分析;利用Qiime(V1.9.1)軟件計(jì)算Chao1、ACE、Shannon、Simpson等指數(shù)[14]。

    Chao1指數(shù)采用式(1)計(jì)算。

    式中1為Chao1指數(shù);obs為測(cè)序分析得到的物種數(shù);1為樣本中數(shù)量為1的數(shù)目;2為樣本中數(shù)量為2的數(shù)目。

    ACE指數(shù)采用式(2)、(3)計(jì)算。

    式中ACE為ACE指數(shù);rare為含有“abund”條序列或少于“abund”條序列的OTU(optical transform unit)數(shù)目;n為含有條序列的OTU數(shù)目;abund為多于“abund”條序列的OTU數(shù)目;abund默認(rèn)為10。

    Shannon指數(shù)采用式(4)計(jì)算。

    式中為Shannon指數(shù);P為樣品中屬于第種個(gè)體的比例。

    Simpson指數(shù)采用式(5)計(jì)算。

    式中為Simposon指數(shù);為物種數(shù)目。

    利用R(V2.15.3)軟件進(jìn)行Beta多樣性指數(shù)組間差異分析。

    1.3 水質(zhì)分析方法

    采用多參數(shù)水質(zhì)測(cè)定儀(北京連華永興科技發(fā)展公司,LH-3BA)分析COD、TN和TP濃度;采用原子吸收光譜儀(日本島津公司,AA-7000)分析污水中Na+、K+濃度;Zeta電位采用微電泳儀(上海中晨數(shù)字技術(shù)設(shè)備有限公司,JS94H2)檢測(cè);脫氫酶活性采用TTC比色法測(cè)定[15];在反應(yīng)器的每個(gè)運(yùn)行周期結(jié)束,利用缺氧區(qū)的活性污泥分析不同條件下反硝化聚磷菌與聚磷菌(DPAOs/PAOs)的比例[16-17],首先取1 L缺氧區(qū)污泥置于三角瓶中,沉降10 min后棄去上清液,污泥再經(jīng)過2.5 h厭氧攪拌后,平均分成2份:①投加KNO3,缺氧條件下運(yùn)行3 h;②好氧條件下運(yùn)行3 h。試驗(yàn)過程中,每隔30 min取樣分析PO43--P和NO3--N的變化,最大缺氧吸磷速率與最大好氧吸磷速率的比值(PURanmax/PURomax),即為DPAOs/PAOs比例。

    厭氧區(qū)釋磷率采用式(6)計(jì)算。

    式中S為厭氧區(qū)釋磷率,%;TPY為厭氧區(qū)總磷濃度,mg/L;TPJ為進(jìn)水總磷濃度,mg/L。

    缺氧區(qū)吸磷率采用式(7)計(jì)算。

    式中X為缺氧區(qū)吸磷率,%;TPQ為缺氧區(qū)總磷濃度,mg/L。

    好氧區(qū)吸磷率采用式(8)計(jì)算。

    式中X為好氧區(qū)吸磷率,%;TPH為好氧區(qū)總磷濃度,mg/L。

    厭氧區(qū)COD去除率計(jì)算方法

    缺氧區(qū)COD去除率計(jì)算方法

    好氧區(qū)COD去除率計(jì)算方法

    2 結(jié)果與分析

    2.1 Na+、K+共存對(duì)污染物去除率的影響

    Na+、K+共存對(duì)A2/O工藝污染物去除率的影響,如圖2所示。隨著進(jìn)水K+濃度的提高,COD去除率增大。當(dāng)進(jìn)水Na+/K+摩爾比為2時(shí),總COD去除率為80%,厭氧區(qū)、缺氧區(qū)和好氧區(qū)COD去除率分別為43%、35%和47%,Na+對(duì)COD去除率產(chǎn)生明顯的抑制作用。Na+在維持細(xì)胞滲透壓和膜電位平衡等方面起著重要作用,過高的鹽濃度會(huì)導(dǎo)致細(xì)胞膜電位的紊亂[18]。外界環(huán)境Na+濃度過高,微生物胞內(nèi)滲透壓低于胞外,細(xì)胞失水,從而影響微生物胞內(nèi)酶的活性,細(xì)胞失活。當(dāng)進(jìn)水Na+/K+摩爾比為1時(shí),總COD去除率為84%,厭氧區(qū)、缺氧區(qū)和好氧區(qū)COD去除率分別為47%、36%和51%,隨著K+濃度的升高,嗜鹽微生物通過吸收K+來平衡外界的Na+鹽度環(huán)境;當(dāng)進(jìn)水Na+/K+摩爾比為0.5時(shí),總COD去除率為86%,厭氧區(qū)、缺氧區(qū)和好氧區(qū)COD去除率繼續(xù)升高,分別為50%、40%和53%。在維持微生物細(xì)胞的滲透壓方面,K+對(duì)細(xì)胞膜通透性的影響大于Na+,一些耐鹽微生物的蛋白質(zhì)和酶需要較高的K+濃度才能維持正常的功能和結(jié)構(gòu)[19]。微生物主要通過細(xì)胞內(nèi)積累大量的K+排出Na+[20],保持細(xì)胞內(nèi)外的滲透壓平衡和離子平衡,減緩高濃度Na+對(duì)細(xì)胞結(jié)構(gòu)和功能的破壞[21],進(jìn)而實(shí)現(xiàn)微生物的耐鹽,K+的增加有利于COD的去除。

    圖2 Na+、K+共存對(duì)A2/O工藝污染物去除率的影響

    當(dāng)進(jìn)水Na+/K+摩爾比為2時(shí),TN去除率為73%,厭氧區(qū)、缺氧區(qū)和好氧區(qū)TN去除率分別為47%、36%和30%;當(dāng)進(jìn)水Na+/K+摩爾比為1時(shí),TN去除率升高至77%,厭氧區(qū)、缺氧區(qū)和好氧區(qū)分別為50%、38%和32%,較高的K+濃度(>50 mg/L)可以減輕Na+對(duì)脫氮微生物的抑制作用[22],當(dāng)進(jìn)水Na+/K+摩爾比為0.5時(shí),TN去除率達(dá)到80%,厭氧區(qū)、缺氧區(qū)和好氧區(qū)的TN去除率升高至53%、39%和35%。本課題組的前期研究表明,Na+濃度的增大(0、10、40 g/L)對(duì)硝化和反硝化過程均產(chǎn)生抑制作用[2]。本研究中,隨著進(jìn)水K+濃度的升高,TN去除率呈增大趨勢(shì),適量的K+可起到促進(jìn)硝化菌和反硝化菌生長的作用,有效提升微生物的脫氮能力,K+的增加有利于TN的去除。

    隨著進(jìn)水Na+/K+摩爾比的降低,A2/O工藝TP去除率逐漸增大。當(dāng)Na+/K+摩爾比為2時(shí),TP去除率為60%,厭氧區(qū)釋磷率為70%,缺氧吸磷率和好氧吸磷率分別為53%和70%;當(dāng)Na+/K+摩爾比為1時(shí),TP去除率為62%,厭氧區(qū)釋磷率為73%,缺氧吸磷率和好氧吸磷率分別為55%和72%;當(dāng)Na+/K+摩爾比為0.5時(shí),TP去除率為65%,厭氧區(qū)釋磷率為74%,缺氧吸磷率和好氧吸磷率分別為58%和75%。有研究表明,隨著鹽度的增加,聚磷菌受到的抑制作用增大[23]。本研究發(fā)現(xiàn),隨著Na+/K+摩爾比的降低,反硝化聚磷菌所占比例增加,反硝化除磷效率提高。K+濃度的提高有利于反硝化聚磷菌數(shù)量的增加,缺氧區(qū)反硝化除磷污泥的活性逐漸增強(qiáng),Na+/K+為2、1、0.5時(shí),DPAOs/PAO分別為43%、45%、和46%。

    2.2 Na+、K+共存對(duì)污泥性質(zhì)的影響

    在不同Na+/K+條件下,厭氧區(qū)、缺氧區(qū)和好氧區(qū)污泥脫氫酶活性的變化,如圖3所示??梢钥闯觯瑫r(shí)投加Na+、K+抑制脫氫酶活性,但隨著進(jìn)水Na+/K+摩爾比的降低,脫氫酶活性呈上升趨勢(shì)。當(dāng)進(jìn)水Na+/K+摩爾比為2時(shí),厭氧區(qū)、缺氧區(qū)和好氧區(qū)脫氫酶活性分別為58.3、58.7和59.4 mg/(L·h),仍低于無Na+、K+運(yùn)行條件下的厭氧區(qū)、缺氧區(qū)和好氧區(qū)脫氫酶活性(分別為62.5、62.3和63.6 mg/(L·h))[24]。當(dāng)進(jìn)水Na+/K+摩爾比為1時(shí),脫氫酶活性升高,厭氧區(qū)、缺氧區(qū)和好氧區(qū)脫氫酶活性分別為56.3、59.1和60.3 mg/(L·h);當(dāng)進(jìn)水Na+/K+摩爾比為0.5時(shí),脫氫酶活性分別為57.1、59.9和60.7 mg/(L·h)。K+濃度的增加有利于缺氧區(qū)和好氧區(qū)脫氫酶活性的提高,污染物去除率上升。

    圖3 不同Na+/K+摩爾下A2/O厭氧區(qū)、缺氧區(qū)和好氧區(qū)污泥的脫氫酶活性

    Na+/K+摩爾比對(duì)Zeta電位的影響,如表2所示。可以看出,隨著進(jìn)水Na+/K+摩爾比的下降,厭氧區(qū)、缺氧區(qū)和好氧區(qū)污泥的Zeta電位均呈上升趨勢(shì)。有研究表明,K+比Na+更容易被吸附,隨著進(jìn)水K+的逐漸增多,污泥絮體所吸附的正電荷增多,引起Zeta電位的升高[25]。圖4表明,隨著進(jìn)水K+濃度的增大,粒徑與重絮凝(flocculation ability,F(xiàn)A)呈增長趨勢(shì),污泥生物絮凝性逐漸增強(qiáng),Na+對(duì)污泥的解絮能力受到抑制[26-27]。

    表2 不同Na+/K+摩爾比條件下Zeta電位的變化

    圖4 Na+/K+摩爾比對(duì)A2/O反應(yīng)器不同區(qū)域污泥粒徑和重絮凝的影響

    2.3 Na+、K+共存對(duì)微生物群落結(jié)構(gòu)的影響

    2.3.1 Na+、K+共存對(duì)微生物多樣性的影響

    本課題組的前期研究表明,Na+質(zhì)量濃度的增大(0、10、40 g/L)導(dǎo)致微生物多樣性逐漸降低,微生物群落生長受到抑制[2]。不同Na+/K+摩爾比下A2/O工藝厭氧區(qū)、缺氧區(qū)和好氧區(qū)微生物多樣性的變化,如表3所示。Chao指數(shù)和ACE指數(shù)可以表示物種豐富度的高低,數(shù)值越高,豐富度越高;Shannon指數(shù)和Simpson指數(shù)可以表示物種多樣性的變化,數(shù)值升高,多樣性升高。隨著進(jìn)水Na+/K+摩爾比的下降,厭氧區(qū)、缺氧區(qū)和好氧區(qū)中Chao指數(shù)、ACE指數(shù)、Shannon指數(shù)和Simpson指數(shù)均降低,群落豐富度和物種多樣性降低,各區(qū)常見物種均減少。這可能是因?yàn)檫M(jìn)水K+濃度升高導(dǎo)致不能適應(yīng)K+鹽度變化的微生物被淘汰,微生物多樣性降低;但是,污染物去除率并未降低(圖2所示)。由此可見,微生物多樣性并不一定是含鹽污水污染物去除的主要制約因素。有研究表明,K+鹽度增加,微生物群落生長受到抑制,微生物多樣性逐漸降低[28-30]。

    表3 微生物多樣性的變化

    注:YM2、YM1、YM0.5、QM2、QM1、QM0.5、HM2、HM1和HM0.5分別代表Na+/K+為2、1、0.5時(shí)采集的厭氧區(qū)、缺氧區(qū)、好氧區(qū)的污泥樣品,下同。

    Note: YM2, YM1, YM0.5, QM2, QM1, QM0.5, HM2, HM1 and HM0.5 represent sludge samples collected in anaerobic, anoxic and oxic zones when the Na+/K+molar ratios are 2, 1 and 0.5, respectively, the same below.

    2.3.2 Na+、K+共存對(duì)A2/O工藝微生物差異性的影響

    當(dāng)進(jìn)水Na+/K+摩爾比分別為2、1和0.5時(shí),厭氧區(qū)、缺氧區(qū)和好氧區(qū)污泥樣品間的微生物群落差異,如圖5所示。通過主坐標(biāo)分析(principal coordinate analysis,PCo A)和非度量多維尺度分析(non-metric multidimensional scaling,NMDS)比較不同生態(tài)系統(tǒng)的Beta多樣性變化,反映單個(gè)樣品間是否存在群落差異[31-32]。每一個(gè)點(diǎn)代表一個(gè)樣本,一個(gè)區(qū)域表示一個(gè)分組,兩點(diǎn)之間距離越近表明兩者的群落結(jié)構(gòu)差異越小。根據(jù)圖5a、b,隨著進(jìn)水Na+/K+摩爾比的降低,樣品間距離加大,樣品分布在不同空間位置上,表明污泥樣品微生物群落間具有顯著差異。

    根據(jù)OTU聚類和微生物豐度對(duì)樣品進(jìn)行聚類分析,結(jié)果如圖5c所示。隨著進(jìn)水Na+/K+摩爾比的降低,微生物群落出現(xiàn)分支,K+的增加導(dǎo)致微生物群落之間的差異性增大,這與Beta多樣性分析結(jié)果相一致;結(jié)合A2/O工藝對(duì)污染物去除率的變化(圖2)可知,進(jìn)水Na+/K+摩爾比對(duì)污染物去除率和微生物群落具有顯著影響。

    2.3.3 Na+、K+共存對(duì)A2/O工藝微生物群落優(yōu)勢(shì)菌的影響

    選取門、屬2個(gè)級(jí)別對(duì)微生物群落結(jié)構(gòu)進(jìn)行分析。A2/O工藝不同區(qū)域活性污泥在門水平上的微生物群落相對(duì)豐度,如圖6a所示。在所有污泥樣品中,變形菌門()均為優(yōu)勢(shì)菌門,相對(duì)豐度最高。當(dāng)進(jìn)水Na+/K+摩爾比為2、1和0.5時(shí),厭氧區(qū)、缺氧區(qū)和好氧區(qū)中變形菌門的相對(duì)豐度升高(升高約30%);擬桿菌門()相對(duì)豐度逐漸降低,厭氧區(qū)擬桿菌門相對(duì)豐度由20.6%降低至9.6%,缺氧區(qū)擬桿菌門相對(duì)豐度由18.2%降低至8.5%,好氧區(qū)擬桿菌門相對(duì)豐度由20.0%降低至7.1%,擬桿菌門不僅在水環(huán)境中廣泛存在[33],還是厭氧消化的主要菌門之一。綠彎菌門()作為具有較強(qiáng)污染物降解能力的主要菌門,其相對(duì)豐度降低,厭氧區(qū)綠彎菌門相對(duì)豐度由9.5%下降至2.0%,缺氧區(qū)綠彎菌門相對(duì)豐度由4.3%下降至3.6%,好氧區(qū)綠彎菌門相對(duì)豐度由7.6%下降至1.9%。為優(yōu)勢(shì)菌門之一,厭氧區(qū)相對(duì)豐度由2.8%下降至1.0%,缺氧區(qū)相對(duì)豐度由1.9%下降至0.9%,好氧區(qū)相對(duì)豐度由2.7%下降至0.7%。

    圖5 微生物群落結(jié)構(gòu)的差異性

    圖6 微生物物種相對(duì)豐度

    圖6b為優(yōu)勢(shì)菌屬水平上的微生物群落分析。隨著進(jìn)水Na+/K+的變化,厭氧區(qū)、缺氧區(qū)和好氧區(qū)陶氏菌屬()的相對(duì)豐度逐漸增大,始終為各區(qū)域的優(yōu)勢(shì)菌屬。當(dāng)Na+/K+摩爾比分別為2、1和0.5時(shí),厭氧區(qū)陶氏菌屬的相對(duì)豐度分別為17.4%、29.3%和43.2%,缺氧區(qū)陶氏菌屬的相對(duì)豐度分別為22.6%、41.3%和42.9%,好氧區(qū)陶氏菌屬的相對(duì)豐度分別為18.8%、37.6%和51.1%。陶氏菌屬作為具有廣泛污染物降解能力的菌屬,其相對(duì)豐度的增大有利于污染物去除率的提高,這與污染物去除率變化相一致。陶氏菌屬和固氮弧菌屬()均為反硝化過程中的重要菌屬[34-35]。隨著進(jìn)水Na+/K+摩爾比的降低,固氮弧菌屬相對(duì)豐度逐漸增大,厭氧區(qū)、缺氧區(qū)和好氧區(qū)固氮弧菌屬的相對(duì)豐度分別由2.1%、2.6%和2.2%增加至4.8%、4.2%和6.0%,從而促進(jìn)反硝化的進(jìn)行。同時(shí),隨著進(jìn)水Na+/K+摩爾比的降低,參與反硝化的陶氏菌屬和固氮弧菌屬增多,未知菌屬所占比例減少,有利于污染物的去除。

    3 結(jié) 論

    1)隨著進(jìn)水Na+/K+摩爾比的下降,化學(xué)需氧量、總氮和總磷去除率升高,A2/O工藝反硝化除磷效率逐漸增強(qiáng)。雖然同時(shí)投加Na+和K+對(duì)A2/O工藝污染物去除產(chǎn)生抑制作用,但進(jìn)水K+濃度的增加緩解了Na+對(duì)微生物的抑制能力。

    2)隨著進(jìn)水Na+/K+摩爾比的下降,厭氧區(qū)、缺氧區(qū)和好氧區(qū)的Zeta電位、污泥粒徑和重絮凝性呈上升趨勢(shì),污泥絮凝性升高,高濃度的K+濃度緩解了Na+的解絮行為。

    3)隨著進(jìn)水Na+/K+摩爾比的下降,厭氧區(qū)、缺氧區(qū)和好氧區(qū)的微生物群落豐富度和多樣性降低,微生物的群落結(jié)構(gòu)差異顯著。優(yōu)勢(shì)菌門始終為變形菌門、擬桿菌門、綠彎菌門、厚壁菌門和綠菌門,陶氏菌屬相對(duì)豐度的增大有利于污染物去除率的提高,固氮弧菌屬相對(duì)豐度的提高有利于促進(jìn)反硝化的進(jìn)行。

    [1] Church J, Hwang J H, Kim K T, et al. Effect of salt type and concentration on the growth and lipid content of Chlorella vulgaris in synthetic saline wastewater for biofuel production[J]. Bioresource Technology, 2017, 243: 147-153.

    [2] 張?zhí)m河,田蕊,陳子成,等. NaCl鹽度對(duì)A2/O工藝去除廢水污染物和系統(tǒng)微生物的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(10):231-237.

    Zhang Lanhe, Tian Rui, Chen Zichen, et al. Effects of NaCl salinity on wastewater pollutants removal and microorganism in A2/O technology process[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2018, 34(10): 231-237. (in Chinese with English abstract)

    [3] 張崇淼,馬宇超,牛全睿,等. 漸增NaCl對(duì)印染廢水處理系統(tǒng)活性污泥微生物的影響[J]. 工業(yè)水處理,2017,37(2):33-37.

    Zhang Chongmiao, Ma Yuchao, Niu Quanrui, et al. Influences of gradual-increase NaCl on the activated sludge microbes in dyeing wastewater treatment systems[J]. Industrial Water Treatment, 2017, 37(2): 33-37. (in Chinese with English abstract)

    [4] Wu G, Guan Y, Zhan X. Effect of salinity on the activity, settling and microbial community of activated sludge in sequencing batch reactors treating synthetic saline wastewater[J].Water Science & Technology, 2008, 58(2): 351-358.

    [5] 趙凱峰,王淑瑩,葉柳,等. NaCl鹽度對(duì)耐鹽活性污泥沉降性能及脫氮的影響[J]. 環(huán)境工程學(xué)報(bào),2010,4(3):570-574.

    Zhao Kaifeng, Wang Shuying, Ye Liu, et al. Effects of salinity on salt-tolerant activated sludge settling performance and nitrogen removal[J]. Chinese Journal of Environmental Engineering, 2010, 4(3): 570-574. (in Chinese with English abstract)

    [6] Mannina G, Capodici M, Cosenza A, et al. Sequential batch membrane bio-reactor for wastewater treatment: The effect of increased salinity[J]. Bioresource Technology, 2016, 209: 205-212.

    [7] Hong J, Li W, Lin B, et al. Deciphering the effect of salinity on the performance of submerged membrane bioreactor for aquaculture of bacterial community[J]. Desalination, 2013, 316(5): 23-30.

    [8] Li J, Shi W, Jiang C, et al. Evaluation of potassium as promoter on anaerobic digestion of saline organic wastewater[J]. Bioresource Technology, 2018, 266: 68-74.

    [9] 鄒小玲,許柯,夏興華,等. NaCl和KCl鹽度對(duì)厭氧污泥的馴化及對(duì)比產(chǎn)甲烷活性的影響[J]. 中國沼氣,2009,27(3):23-25.

    Zou Xiaoling, Xu Ke, Xia Xinghua, et al. NaCl and KCl salinity acclimations of anaerobic sludge and their effects on specific methanogenic activity[J]. China Biogas, 2009, 27(3): 23-25. (in Chinese with English abstract)

    [10] 鄒小玲,許柯,丁麗麗,等. NaCl和KCl對(duì)厭氧污泥抑制的動(dòng)力學(xué)研究[J]. 化工環(huán)保,2009,29(5):394-397.

    Zou Xiaoling, Xu Ke, Ding Lili, et al. Kinetics study of inhibiting action of NaCl and KCl on anaerobic sludge[J]. Environmental Protection of Chemical Industry, 2009, 29(5): 394-397. (in Chinese with English abstract)

    [11] LiuL, Yan H, Tan W, et al. Influence of Na+, K+, Mg2+, Ca2+, and Fe3+on filterability and settleability of drilling sludge[J]. Chinese Journal of Chemical Engineering, 2016, 25(5):658-664.

    [12] 楊紅薇,陳佼,張建強(qiáng). K+、Ca2+、Mg2+對(duì)高鹽肝素廢水處理的影響[J]. 環(huán)境工程學(xué)報(bào),2014,8(10):4267-4272.

    Yang Hongwei, Chen Jiao, Zhang Jianqiang. Effects of K+,Ca2+,Mg2+on high salt heparin wastewater treatment[J]. Chinese Journal of Environmental Engineering, 2014, 8(10):4267-4272. (in Chinese with English abstract)

    [13] Hwang K L, Bang C H, Zoh K D. Characteristics of methane and nitrous oxide emissions from the wastewater treatment plant[J]. Bioresource Technology, 2016, 214: 881-884.

    [14] 徐華勤,肖潤林,鄒冬生,等. 長期施肥對(duì)茶園土壤微生物群落功能多樣性的影響[J].生態(tài)學(xué)報(bào),2007,27(8):3355-3361.

    Xu Huaqin, Xiao Runlin, Zou Dongsheng, et al. Effects of long-term fertilization on functional diversity of soil microbial community of the tea plantation[J]. Acta Ecologica Sinica, 2007, 27(8): 3355-3361. (in Chinese with English abstract)

    [15] 魏民,鄭國臣,李建政,等. 表征水體中生物活性及脫氫酶檢測(cè)方法研究[J]. 東北水利水電,2012,30(8):43-46.

    Wei Min, Zheng Guochen, Li Jianzheng, et al. Study on detection method of biological activity and dehydrogenase in water[J]. Water Resources & Hydropower of Northeast China, 2012, 30(8): 43-46. (in Chinese with English abstract)

    [16] Guo H, Zhou J, Zhang S, et al. Characteristics of nitrogen and phosphorus removal in a sequencing batch reactor[J]. Journal of Environmental Sciences, 2011, 23, Suppl(11): S110-S113.

    [17] Wu D, Ekama G A, Wang H G, et al. Simultaneous nitrogen and phosphorus removal in the sulfur cycle-associated enhanced biological phosphorus removal (EBPR) process[J]. Water Research, 2014, 49(2): 251-264.

    [18] 林兵. 鹽脅迫下好氧顆粒污泥微生物種群變化及生態(tài)響應(yīng)機(jī)制[D]. 鎮(zhèn)江:江蘇科技大學(xué),2014.

    Lin Bing. Change of Microbial Community Structure and Ecology Reponse Mechanism of Aerobic Granular Sludge under Salt Stress[D]. Zhenjiang: Jiangsu University of Science and Technology, 2014. (in Chinese with English abstract)

    [19] 李志華,王曉昌,王耀東. 含鹽量對(duì)好氧顆粒污泥形成過程的影響[J]. 環(huán)境工程學(xué)報(bào),2008,2(9):1228-1230.

    Li Zhihua, Wang Xiaochang, Wang Yaodong. Influence of sodium chloride on formation of aerobic granules[J]. Chinese Journal of Environmental Engineering, 2008, 2(9): 1228-1230. (in Chinese with English abstract)

    [20] 陶衛(wèi)平. 嗜鹽菌的嗜鹽機(jī)制[J]. 生物學(xué)通報(bào),1996,31(1):23-24.

    [21] Shilo, Moshe. Strategies of Microbial Life in Extreme Environments [M]. Verlag Chemie, 1979.

    [22] 何健. 高鹽難降解工業(yè)廢水微生物處理的污泥馴化研究與應(yīng)用[D]. 南京:南京農(nóng)業(yè)大學(xué),2000.

    He Jian. Sludge Acclimation of Microbiological Treatment of Hypersaline Refractory Waste Water[D]. Nanjing: Agricultural University of Nanjing, 2000. (in Chinese with English abstract)

    [23] 李玲玲,周鵬. 活性污泥中功能性菌群抗鹽度沖擊性能研究[J]. 環(huán)境工程學(xué)報(bào),2010,4(1):105-109.

    Li Lingling, Zhou Peng. Effect of salinity shock on specific microorganisms in conventional activated sludge[J]. Chinese Journal of Environmental Engineering, 2010, 4(1): 105-109. (in Chinese with English abstract)

    [24] 張?zhí)m河,田蕊,郭靜波,等. NaCl鹽度對(duì)A2/O工藝缺氧區(qū)胞外聚合物及生物絮凝性的影響[J].環(huán)境科學(xué),2018,39(9):4281-4288.

    Zhang Lanhe, Tian Rui, Guo Jingbo, et al. Effect of NaCl salinity on extracellular polymeric substances and bioflocculation of anoxic sludge in A2/O process[J]. Environmental Sciences, 2018, 39(9): 4281-4288. (in Chinese with English abstract)

    [25] Zhang W, Zheng J, Zheng P, et al. Sludge-Derived biochar for arsenic(iii) immobilization: effects of solution chemistry on sorption behavior[J]. Journal of Environmental Quality, 2015, 44(4): 1119-1126.

    [26] Nouha K, Kumar R S, Balasubramanian S, et al. Critical review of EPS production, synthesis and composition for sludge flocculation[J]. Journal of Environmental Sciences, 2017, 66(4): 225-245.

    [27] Amatore C, Jutand A, Le?Duc G. Mechanistic origin of antagonist effects of usual anionic bases (OH?, CO32?) as modulated by their Countercations (Na+, Cs+, K+) in palladium-catalyzed suzuki-miyaura reactions[J]. Chemistry, 2012, 18(21): 6616-6625.

    [28] Guo X, Miao Y, Wu B, et al. Correlation between microbial community structure and biofouling as determined by analysis of microbial community dynamics[J]. Bioresource Technology, 2015, 197: 99-105.

    [29] 周貴忠,許碩,姚倩,等. 不同鹽度下活性污泥中微生物群落變化規(guī)律及其處理模擬染料廢水[J]. 環(huán)境科學(xué),2017,38(7):2972-2977.

    Zhou Guizhong, Xu Shuo, Yao Qian, et al. Influence of salinity on microbial community in activated sludge and its application in simulated dye wastewater treatment[J]. Environmental Science, 2017, 38(7): 2972-2977. (in Chinese with English abstract)

    [30] Wu G, Guan Y, Zhan X. Effect of salinity on the activity, settling and microbial community of activated sludge in sequencing batch reactors treating synthetic saline wastewater[J]. Water Sci Technol, 2008, 58: 351-358.

    [31] 李墨青. 納米銀對(duì)SBR系統(tǒng)水處理效能及微生物菌群的影響研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué),2014.

    Li Moqing. Impacts of Silver Nanoparticles on Water Treatment Efficiencies of SBR system and Microbial Communities[D]. Harbin: Harbin Institute of Technology, 2014. (in Chinese with English abstract)

    [32] Hallenbeck P C, Liu Y. Recent advances in hydrogen production by photosynthetic bacteria[J]. International Journal of Hydrogen Energy, 2016, 41(7): 4446-4454.

    [33] Lin P Y, Chen H L, Huang C T, et al. Biofilm production, use of intravascular indwelling catheters and inappropriate antimicrobial therapy as predictors of fatality in Chryseobacterium meningosepticum bacteraemia[J]. International Journal of Antimicrobial Agents, 2010, 36(5): 436-440.

    [34] Shen L D, Zheng P H, Ma S J. Nitrogen loss through anaerobic ammonium oxidation in agricultural drainage ditches[J]. Biology & Fertility of Soils, 2016, 52(2): 127-136.

    [35] Tian M, Zhao F, Shen X, et al. The first metagenome of activated sludge from full-scale anaerobic/anoxic/oxic (A2O) nitrogen and phosphorus removal reactor using Illumina sequencing[J]. Acta Scientiae Circumstantiae, 2015, 35(9): 181-190.

    Effect of coexistence of Na+and K+on sludge properties and microbial community structure in A2/O process

    Zhang Lanhe1,2, Zheng Jing1, Tian Rui1, Chen Zicheng1※, Guo Jingbo3, Jia Yanping1, Li Zheng1, Zhang Jian1

    (1.132012,; 2.,130118; 3.,132012)

    A large amount of salty wastewater is produced from industrial processes such as food processing, tanning, textiles, aquaculture and oil production. After the salty wastewater entered the activated sludge system and contacted with the activated sludge, the sedimentation, microbial activity and community structure of activated sludge were influenced, which could result in the change of pollutant removal efficiency. A2/O process could achieve simultaneously nitrogen and phosphorus removal, but it was unclear how Na+and K+influenced the microbial communities and the removal efficiency of pollutants in the anaerobic, anoxic and oxic zones. In order to reveal the removal mechanism of the pollutants and the characteristics of activated sludge in the biological treatment system under the coexisting multiple metal ions, the effect of the coexisting Na+and K+on the removal efficiency of pollutants was investigated and the properties of sludge and microbial community in the anaerobic, anoxic and oxic zones were analyzed by high-throughput sequencing. Combined with the changes of removal efficiency of nitrogen and phosphorus and sludge properties, the succession regulation of dominant populations were explored based on the analysis of microbial community structure under different Na+/K+molar ratios to distinguish the effect of the coexistence of Na+and K+on the removal efficiency of the pollutants in the salty wastewater from the point of microorganisms. The results showed that when the influent Na+/K+molar ratio was 2, 1 and 0.5, the removal efficiencies of COD were 80%, 84% and 86%, respectively. The removal efficiencies of TN were 73%, 77% and 80%, respectively. The increase of K+concentration alleviated the inhibition of Na+on the removal efficiency of COD and TN. The release rates of TP in anaerobic areas were 70%, 73% and 74%, respectively. The phosphorus uptake rates in the anoxic zone were 53%, 55% and 58%, respectively. The phosphorus uptake rates in the oxic zone were 70%, 72% and 75%, respectively. The Zeta potential in anaerobic zone is -25.4, -23.2 and -14.7 mV, respectively. The Zeta potential in the anoxic zone is -33.0, -26.6 and -13.7 mV and the Zeta potential in oxic zone is -30.4, -18.6 and -11.0 mV, respectively. The positive charge adsorbed by sludge floc increased. Moreover, the biological flocculation of sludge increased gradually and the deflocculation ability of Na+was inhibited. The richness and diversity of microbial communities in anaerobic, anoxic and oxic zones decreased and the differences of microbial communities were significant. The relative abundance ofincreased by about 30% and the relative abundance ofanddecreased gradually. As the dominant genera, the relative abundance of nitrogen-fixingandgradually increased and the proportion of unknown bacteria decreased, which were beneficial to the removal of pollutants. Therefore, the increase of K+concentration contributed to increase the removal efficiency of nitrogen and phosphorus under the coexistence of Na+and K+. It was also beneficial to the improvement of dehydrogenase activity of sludge in anoxic zone and oxic zone and the enhancement of particle size and flocculation ability (FA) of sludge.

    metal ions; wastewater; sludge; bioreactor; microbial community structure; salinity

    2018-11-07

    2019-04-07

    國家自然科學(xué)基金(51678119,51508073);吉林省科技發(fā)展計(jì)劃項(xiàng)目(20180201016SF,20180101309JC,20180101079JC);吉林省教育廳科學(xué)技術(shù)研究項(xiàng)目(JJKH20180454KJ,JJKH20180453KJ)

    張?zhí)m河,博士,教授,主要研究方向?yàn)槲鬯锩摰准夹g(shù)。Email:zhanglanhe@163.com

    陳子成,博士,副教授,主要研究方向?yàn)樗幚砑夹g(shù)與理論。Email:chenzicheng@126.com

    10.11975/j.issn.1002-6819.2019.11.024

    X703

    A

    1002-6819(2019)-11-0206-08

    張?zhí)m河,鄭 晶,田 蕊,陳子成,郭靜波,賈艷萍,李 正,張 健. Na+和K+共存對(duì)A2/O工藝脫氮除磷效果及污泥性質(zhì)的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(11):206-213. doi:10.11975/j.issn.1002-6819.2019.11.024 http://www.tcsae.org

    Zhang Lanhe, Zheng Jing, Tian Rui, Chen Zicheng, Guo Jingbo, Jia Yanping, Li Zheng, Zhang Jian. Effect of coexistence of Na+and K+on sludge properties and microbial community structure in A2/O process[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(11): 206-213. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.11.024 http://www.tcsae.org

    猜你喜歡
    氧區(qū)磷率鹽度
    自動(dòng)清洗高效生物膜反應(yīng)器設(shè)計(jì)與性能
    化工管理(2024年1期)2024-01-17 08:49:20
    一種污水處理系統(tǒng)及處理方法與應(yīng)用
    能源化工(2021年6期)2021-12-30 15:41:26
    轉(zhuǎn)爐脫磷工藝中單渣法的研究進(jìn)展
    組合改性石灰石對(duì)農(nóng)村分散性生活污水除磷性能研究
    70 t轉(zhuǎn)爐“留渣+雙渣”脫磷工藝研究與實(shí)踐
    鹽度和pH對(duì)細(xì)角螺耗氧率和排氨率的影響
    鹽度脅迫對(duì)入侵生物福壽螺的急性毒性效應(yīng)
    適用于高鹽度和致密巖層驅(qū)油的表面活性劑
    多段多級(jí)AO除磷脫氮工藝
    馬鋼轉(zhuǎn)爐雙渣法脫磷工藝生產(chǎn)實(shí)踐
    欧美日韩精品成人综合77777| 亚洲综合色惰| 国产日本99.免费观看| 日日摸夜夜添夜夜爱| 色吧在线观看| 久久这里只有精品中国| 看非洲黑人一级黄片| 最好的美女福利视频网| 此物有八面人人有两片| 欧美日本亚洲视频在线播放| 97人妻精品一区二区三区麻豆| 亚洲自拍偷在线| 国产av在哪里看| 亚洲av中文av极速乱| 99国产极品粉嫩在线观看| 亚洲成人精品中文字幕电影| 中文字幕熟女人妻在线| 熟女人妻精品中文字幕| 免费观看精品视频网站| 日本黄大片高清| 最新中文字幕久久久久| 免费观看人在逋| 国产精品久久久久久精品电影小说 | 超碰av人人做人人爽久久| 一进一出抽搐动态| 99热精品在线国产| 夫妻性生交免费视频一级片| 色综合亚洲欧美另类图片| 岛国在线免费视频观看| 亚洲精品456在线播放app| 日韩成人伦理影院| 亚洲精品日韩在线中文字幕 | kizo精华| 国产精品三级大全| 亚洲四区av| 级片在线观看| 波野结衣二区三区在线| 国产成人91sexporn| 欧美在线一区亚洲| 欧洲精品卡2卡3卡4卡5卡区| 岛国毛片在线播放| 午夜a级毛片| 日日干狠狠操夜夜爽| 色视频www国产| 日韩欧美在线乱码| 日本与韩国留学比较| 国产成年人精品一区二区| 亚洲精品粉嫩美女一区| 欧美一区二区国产精品久久精品| АⅤ资源中文在线天堂| 色综合亚洲欧美另类图片| 国产一区亚洲一区在线观看| 波多野结衣高清作品| 变态另类丝袜制服| 久久久久久大精品| 国产黄片视频在线免费观看| 亚洲国产精品成人综合色| 一级毛片我不卡| 精品欧美国产一区二区三| 少妇人妻精品综合一区二区 | 日韩国内少妇激情av| 在线观看免费视频日本深夜| 草草在线视频免费看| 丝袜美腿在线中文| 乱码一卡2卡4卡精品| 波野结衣二区三区在线| 五月玫瑰六月丁香| 麻豆国产av国片精品| 日本一本二区三区精品| 久久久久久九九精品二区国产| 五月玫瑰六月丁香| 最近的中文字幕免费完整| 成人美女网站在线观看视频| 悠悠久久av| 在线免费观看的www视频| av在线播放精品| 一级黄片播放器| 免费人成在线观看视频色| 高清午夜精品一区二区三区 | 欧美潮喷喷水| 天天躁夜夜躁狠狠久久av| 国产亚洲91精品色在线| 久99久视频精品免费| 亚洲欧美日韩高清在线视频| 日本免费一区二区三区高清不卡| АⅤ资源中文在线天堂| 亚洲国产精品成人综合色| 久久精品夜色国产| 又黄又爽又刺激的免费视频.| 久久久久免费精品人妻一区二区| 91aial.com中文字幕在线观看| 国产精品爽爽va在线观看网站| 国产精品嫩草影院av在线观看| 国产高清三级在线| 在线观看av片永久免费下载| 中文字幕熟女人妻在线| 婷婷色av中文字幕| 欧美激情国产日韩精品一区| 欧美日本亚洲视频在线播放| 久久韩国三级中文字幕| 99久久人妻综合| 国语自产精品视频在线第100页| 午夜a级毛片| 夜夜爽天天搞| 国产精品麻豆人妻色哟哟久久 | 国产精品.久久久| 国产不卡一卡二| 欧美在线一区亚洲| 老女人水多毛片| 寂寞人妻少妇视频99o| 亚洲熟妇中文字幕五十中出| 国产麻豆成人av免费视频| 亚洲人成网站在线观看播放| 国产精品福利在线免费观看| 一区二区三区免费毛片| 成人美女网站在线观看视频| av在线观看视频网站免费| 一级二级三级毛片免费看| 干丝袜人妻中文字幕| 中国美白少妇内射xxxbb| 中文欧美无线码| 欧美激情在线99| 性欧美人与动物交配| 日本熟妇午夜| 久久久久久久久久黄片| 国内精品宾馆在线| 伦精品一区二区三区| www日本黄色视频网| 一区二区三区四区激情视频 | 国产男人的电影天堂91| 日本黄色片子视频| 欧美激情久久久久久爽电影| 久久久久久伊人网av| 禁无遮挡网站| 中文资源天堂在线| 国产伦精品一区二区三区视频9| 日韩欧美精品免费久久| 黄色欧美视频在线观看| 神马国产精品三级电影在线观看| 一夜夜www| 国产亚洲欧美98| 中国美白少妇内射xxxbb| 国产 一区精品| 日日摸夜夜添夜夜爱| 国内久久婷婷六月综合欲色啪| 亚洲人成网站在线播放欧美日韩| av又黄又爽大尺度在线免费看 | 国产美女午夜福利| 久久这里有精品视频免费| 国产精品久久久久久亚洲av鲁大| 久久6这里有精品| 精品日产1卡2卡| 婷婷精品国产亚洲av| 亚洲18禁久久av| 久久久久网色| 亚洲成av人片在线播放无| 亚洲成人中文字幕在线播放| 国产精品久久视频播放| 成人美女网站在线观看视频| 久久99精品国语久久久| 一个人观看的视频www高清免费观看| 一级二级三级毛片免费看| 久久久久九九精品影院| 免费av观看视频| 国产午夜精品久久久久久一区二区三区| 国产高潮美女av| 色哟哟·www| 51国产日韩欧美| 久久久久免费精品人妻一区二区| 3wmmmm亚洲av在线观看| 日本一二三区视频观看| 国产亚洲精品av在线| 亚洲成人中文字幕在线播放| 99久国产av精品国产电影| 久久久久久久亚洲中文字幕| 国产亚洲5aaaaa淫片| 一级毛片久久久久久久久女| 国产av一区在线观看免费| 国产一级毛片七仙女欲春2| a级毛片a级免费在线| 久久久久性生活片| 特级一级黄色大片| 亚洲国产日韩欧美精品在线观看| 51国产日韩欧美| 九九爱精品视频在线观看| 特级一级黄色大片| 边亲边吃奶的免费视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 12—13女人毛片做爰片一| www.色视频.com| 久久韩国三级中文字幕| 婷婷色av中文字幕| 亚洲成av人片在线播放无| 大型黄色视频在线免费观看| 麻豆国产97在线/欧美| 99久久精品国产国产毛片| 亚洲欧美成人精品一区二区| 色吧在线观看| 午夜精品国产一区二区电影 | 久久久精品欧美日韩精品| 中文字幕免费在线视频6| 免费观看在线日韩| 少妇猛男粗大的猛烈进出视频 | 日韩欧美三级三区| 国产精品蜜桃在线观看 | 欧美激情在线99| 国内精品久久久久精免费| 国产三级在线视频| 精品欧美国产一区二区三| 两性午夜刺激爽爽歪歪视频在线观看| 成年女人永久免费观看视频| 乱码一卡2卡4卡精品| 亚洲真实伦在线观看| 国产精品久久视频播放| 91精品一卡2卡3卡4卡| 丰满的人妻完整版| 99久国产av精品| 午夜福利在线观看免费完整高清在 | 乱码一卡2卡4卡精品| av在线观看视频网站免费| 日韩高清综合在线| 日韩国内少妇激情av| 又粗又硬又长又爽又黄的视频 | 永久网站在线| 99久久无色码亚洲精品果冻| 欧美变态另类bdsm刘玥| www.色视频.com| 免费在线观看成人毛片| 免费观看在线日韩| 男女视频在线观看网站免费| 日日摸夜夜添夜夜添av毛片| 亚洲熟妇中文字幕五十中出| 麻豆精品久久久久久蜜桃| 国产高清不卡午夜福利| 久久6这里有精品| 又黄又爽又刺激的免费视频.| 免费看光身美女| 小说图片视频综合网站| 99九九线精品视频在线观看视频| 欧美激情在线99| 亚洲av中文字字幕乱码综合| 99热网站在线观看| 男人和女人高潮做爰伦理| 伊人久久精品亚洲午夜| 国产精品人妻久久久影院| 欧美成人精品欧美一级黄| 一个人看的www免费观看视频| 国产亚洲91精品色在线| 国产成人aa在线观看| 欧美不卡视频在线免费观看| 欧美精品国产亚洲| 午夜精品国产一区二区电影 | 国产熟女欧美一区二区| 综合色av麻豆| 欧美一级a爱片免费观看看| 国产午夜福利久久久久久| 搡老妇女老女人老熟妇| 免费黄网站久久成人精品| 亚洲激情五月婷婷啪啪| 校园人妻丝袜中文字幕| 亚洲精华国产精华液的使用体验 | 国产v大片淫在线免费观看| 免费av不卡在线播放| 噜噜噜噜噜久久久久久91| 蜜桃亚洲精品一区二区三区| 天堂中文最新版在线下载 | 中国国产av一级| 能在线免费观看的黄片| 精品久久久噜噜| 男女啪啪激烈高潮av片| 国产成人影院久久av| 99热精品在线国产| 国产精品,欧美在线| 热99re8久久精品国产| 精品不卡国产一区二区三区| 国产视频首页在线观看| 亚洲av成人精品一区久久| 亚洲精品乱码久久久久久按摩| 午夜老司机福利剧场| 免费观看精品视频网站| 久久精品国产清高在天天线| 人妻系列 视频| av在线蜜桃| 嫩草影院精品99| 欧美不卡视频在线免费观看| 九九久久精品国产亚洲av麻豆| 一级毛片我不卡| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲自偷自拍三级| 国产综合懂色| 2021天堂中文幕一二区在线观| 久久国内精品自在自线图片| 国产成人午夜福利电影在线观看| 亚洲国产欧美人成| 久久人妻av系列| 91午夜精品亚洲一区二区三区| 日本一二三区视频观看| 免费看光身美女| 中文字幕熟女人妻在线| 久久人妻av系列| 久久亚洲精品不卡| 国产精品日韩av在线免费观看| 成年免费大片在线观看| 国产激情偷乱视频一区二区| 亚洲人成网站在线观看播放| 日日啪夜夜撸| 麻豆国产97在线/欧美| 可以在线观看的亚洲视频| 天堂影院成人在线观看| 日本熟妇午夜| 国产极品天堂在线| or卡值多少钱| 中文字幕精品亚洲无线码一区| 国产大屁股一区二区在线视频| 日本黄色视频三级网站网址| 午夜a级毛片| 一个人观看的视频www高清免费观看| 51国产日韩欧美| 精品久久久久久久人妻蜜臀av| 国产精品无大码| 能在线免费观看的黄片| av天堂中文字幕网| 欧美性猛交黑人性爽| 国产精品蜜桃在线观看 | 久久久午夜欧美精品| 麻豆成人午夜福利视频| 黄色一级大片看看| 一级毛片电影观看 | 国产淫片久久久久久久久| 黄色日韩在线| 久久久久久久午夜电影| 久久精品国产亚洲av天美| av在线观看视频网站免费| 亚洲欧美清纯卡通| 成人永久免费在线观看视频| 国产精品国产三级国产av玫瑰| 国产精品久久电影中文字幕| 天天躁夜夜躁狠狠久久av| 69av精品久久久久久| 欧美日韩国产亚洲二区| eeuss影院久久| 亚洲最大成人手机在线| 九草在线视频观看| 一个人观看的视频www高清免费观看| 国产精品一二三区在线看| 精品少妇黑人巨大在线播放 | 日本色播在线视频| 久久久色成人| 亚洲av二区三区四区| 亚洲欧美日韩卡通动漫| 国产激情偷乱视频一区二区| 在线免费观看的www视频| 青春草国产在线视频 | www.色视频.com| 最近2019中文字幕mv第一页| 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 国内精品久久久久精免费| 免费观看a级毛片全部| 99热这里只有是精品在线观看| 国产熟女欧美一区二区| 嫩草影院新地址| 国产成人一区二区在线| 美女黄网站色视频| 99久久人妻综合| 欧美又色又爽又黄视频| 日韩高清综合在线| 特级一级黄色大片| 精品久久久久久久久av| 国产精品永久免费网站| 欧美日韩在线观看h| 夜夜夜夜夜久久久久| 成年女人看的毛片在线观看| 欧美zozozo另类| 久久久精品大字幕| 亚洲av免费在线观看| 国产精品蜜桃在线观看 | 不卡视频在线观看欧美| 嫩草影院精品99| 级片在线观看| 免费看a级黄色片| 成年免费大片在线观看| 91午夜精品亚洲一区二区三区| 免费人成视频x8x8入口观看| 少妇高潮的动态图| 激情 狠狠 欧美| 久久久a久久爽久久v久久| 亚洲一级一片aⅴ在线观看| 国产精品伦人一区二区| 国产成人精品久久久久久| 欧美日韩乱码在线| 看免费成人av毛片| 国产精品爽爽va在线观看网站| 97超碰精品成人国产| 日韩三级伦理在线观看| 国内精品久久久久精免费| 久久久成人免费电影| av黄色大香蕉| 亚洲欧美成人综合另类久久久 | 国产v大片淫在线免费观看| 日韩人妻高清精品专区| 久久99热6这里只有精品| 亚洲中文字幕一区二区三区有码在线看| 亚洲熟妇中文字幕五十中出| 国产av一区在线观看免费| 99国产极品粉嫩在线观看| 国产亚洲欧美98| av在线蜜桃| 国产真实乱freesex| 亚洲中文字幕日韩| 亚洲国产精品成人久久小说 | 男女视频在线观看网站免费| 亚洲在久久综合| 69av精品久久久久久| 国产伦一二天堂av在线观看| 狠狠狠狠99中文字幕| 看片在线看免费视频| 淫秽高清视频在线观看| 国产欧美日韩精品一区二区| 一级毛片电影观看 | 老师上课跳d突然被开到最大视频| av免费在线看不卡| 色尼玛亚洲综合影院| 美女黄网站色视频| 别揉我奶头 嗯啊视频| 99热网站在线观看| 直男gayav资源| 高清日韩中文字幕在线| 精品久久久久久久久久久久久| 日韩三级伦理在线观看| 国产乱人偷精品视频| av在线观看视频网站免费| 乱系列少妇在线播放| 亚洲中文字幕日韩| 国产视频内射| 国产一区二区三区av在线 | 桃色一区二区三区在线观看| 2021天堂中文幕一二区在线观| 成人性生交大片免费视频hd| 永久网站在线| 麻豆乱淫一区二区| 国产精品久久久久久亚洲av鲁大| 又爽又黄无遮挡网站| 在线观看av片永久免费下载| 国产淫片久久久久久久久| 亚洲精品日韩在线中文字幕 | 久久久午夜欧美精品| 麻豆成人午夜福利视频| 麻豆久久精品国产亚洲av| 人人妻人人澡人人爽人人夜夜 | 狂野欧美激情性xxxx在线观看| 日韩成人伦理影院| 毛片一级片免费看久久久久| 国产精品一区二区三区四区免费观看| 春色校园在线视频观看| 日韩一本色道免费dvd| 国产精品嫩草影院av在线观看| 亚洲精品亚洲一区二区| 午夜激情福利司机影院| 国产一区二区三区av在线 | 亚洲精品成人久久久久久| 日韩一区二区三区影片| 久久精品国产自在天天线| 夜夜夜夜夜久久久久| 午夜a级毛片| 成人亚洲精品av一区二区| 精品一区二区免费观看| 九九爱精品视频在线观看| 黑人高潮一二区| 国产又黄又爽又无遮挡在线| 国产精品三级大全| 欧美最黄视频在线播放免费| 亚洲内射少妇av| 亚洲三级黄色毛片| 国产精品精品国产色婷婷| 看片在线看免费视频| 亚洲精品粉嫩美女一区| a级一级毛片免费在线观看| 久久久久久久久大av| 99热精品在线国产| 精品久久国产蜜桃| 日韩,欧美,国产一区二区三区 | 一本久久精品| av.在线天堂| 一进一出抽搐动态| 亚洲av中文av极速乱| 欧美极品一区二区三区四区| 26uuu在线亚洲综合色| 国产在视频线在精品| 婷婷六月久久综合丁香| 久久这里有精品视频免费| 亚洲欧洲日产国产| 国产毛片a区久久久久| 亚洲七黄色美女视频| 亚洲电影在线观看av| 一级毛片我不卡| 91午夜精品亚洲一区二区三区| 如何舔出高潮| 国产一区二区在线观看日韩| 精品久久久噜噜| 亚洲av免费在线观看| 在线播放无遮挡| 中国国产av一级| 婷婷精品国产亚洲av| 夜夜爽天天搞| 国产乱人视频| 久99久视频精品免费| 精品人妻视频免费看| 性欧美人与动物交配| 色吧在线观看| 精品久久久噜噜| 我要搜黄色片| 亚洲最大成人中文| 国产一区二区激情短视频| 99久久精品一区二区三区| 九色成人免费人妻av| 精品少妇黑人巨大在线播放 | 免费看美女性在线毛片视频| 精品久久久久久久人妻蜜臀av| 99国产精品一区二区蜜桃av| 人妻夜夜爽99麻豆av| 春色校园在线视频观看| 成人毛片a级毛片在线播放| 亚洲一区高清亚洲精品| 少妇丰满av| 欧美潮喷喷水| 我要看日韩黄色一级片| 亚洲人成网站在线播放欧美日韩| 久久午夜亚洲精品久久| 亚洲欧美精品专区久久| 国产一级毛片在线| 国产精品国产三级国产av玫瑰| 欧美日本视频| 日本免费a在线| 亚洲电影在线观看av| 久久国内精品自在自线图片| 亚洲成人精品中文字幕电影| 99久国产av精品| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 一个人观看的视频www高清免费观看| 白带黄色成豆腐渣| 成人av在线播放网站| 亚洲精品日韩在线中文字幕 | eeuss影院久久| 日韩大尺度精品在线看网址| 久久精品影院6| 国产精品国产高清国产av| 美女内射精品一级片tv| 国产精品综合久久久久久久免费| 午夜福利成人在线免费观看| 国产精品,欧美在线| 麻豆成人av视频| 69av精品久久久久久| 在线观看av片永久免费下载| 少妇熟女欧美另类| 精品一区二区免费观看| 亚洲人成网站高清观看| 狠狠狠狠99中文字幕| 内射极品少妇av片p| 69av精品久久久久久| 两个人的视频大全免费| 热99re8久久精品国产| 免费av毛片视频| 久久久精品大字幕| 欧美zozozo另类| 国产精品乱码一区二三区的特点| 国产男人的电影天堂91| 国产精品久久视频播放| 麻豆av噜噜一区二区三区| 日本黄大片高清| 超碰av人人做人人爽久久| 啦啦啦韩国在线观看视频| 最好的美女福利视频网| 久久99热6这里只有精品| 一进一出抽搐gif免费好疼| 一个人观看的视频www高清免费观看| 免费av不卡在线播放| av卡一久久| 久久精品91蜜桃| 夜夜爽天天搞| 免费看光身美女| 欧美成人a在线观看| 可以在线观看毛片的网站| 男女边吃奶边做爰视频| 最近最新中文字幕大全电影3| 蜜臀久久99精品久久宅男| 久久久国产成人精品二区| 欧美潮喷喷水| 在线免费观看不下载黄p国产| 12—13女人毛片做爰片一| 国产日本99.免费观看| 欧美精品一区二区大全| 最近的中文字幕免费完整| 国产精品久久久久久久电影| 久久久久网色| 麻豆久久精品国产亚洲av| 久久99热6这里只有精品| 国产极品天堂在线| 婷婷色av中文字幕| 波多野结衣高清无吗| 午夜老司机福利剧场| 免费观看在线日韩| 亚洲精品久久国产高清桃花| 天堂中文最新版在线下载 | 91狼人影院| a级毛片a级免费在线| 亚洲欧美日韩东京热| 最近中文字幕高清免费大全6| av天堂中文字幕网| 91久久精品国产一区二区成人| 欧美日本视频| 岛国在线免费视频观看| 一个人看的www免费观看视频| 桃色一区二区三区在线观看|