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

    牛糞與玉米秸稈配比對混合厭氧發(fā)酵速率的影響

    2016-10-28 09:19:25姜宗姍平佳芃郭延凱陳興華郭建博
    河北科技大學(xué)學(xué)報 2016年4期

    姜宗姍,平佳芃,郭延凱,陳興華,郭建博,廉 靜

    (1.河北科技大學(xué)環(huán)境科學(xué)與工程學(xué)院,河北石家莊 050018;2.河北省污染防治生物技術(shù)實驗室,河北石家莊 050018;3.杭州能源環(huán)境工程有限公司,浙江杭州 310000;4.天津城建大學(xué)環(huán)境與市政工程學(xué)院,天津  300384)

    ?

    牛糞與玉米秸稈配比對混合厭氧發(fā)酵速率的影響

    姜宗姍1,2,3,平佳芃1,2,郭延凱1,2,陳興華3,郭建博1,4,廉靜1,2

    (1.河北科技大學(xué)環(huán)境科學(xué)與工程學(xué)院,河北石家莊050018;2.河北省污染防治生物技術(shù)實驗室,河北石家莊050018;3.杭州能源環(huán)境工程有限公司,浙江杭州310000;4.天津城建大學(xué)環(huán)境與市政工程學(xué)院,天津 300384)

    為了研究牛糞與玉米秸稈的不同配比(以下簡稱糞稈配比,以揮發(fā)性固體(VS)計算)對厭氧發(fā)酵速率的影響,初步確定影響發(fā)酵產(chǎn)氣限速的因素,在中溫(38±1)℃條件下,以厭氧顆粒污泥為接種物,以不同配比的牛糞與玉米秸稈為底物進行混合厭氧發(fā)酵,測定日產(chǎn)氣量、沼氣中CH4與CO2的濃度、沼液中的總碳(TC)、揮發(fā)性脂肪酸(VFAs)濃度和pH值。結(jié)果表明,糞稈配比為2∶1時,水解速率常數(shù)、單位基質(zhì)產(chǎn)甲烷量和生物轉(zhuǎn)化產(chǎn)甲烷效率最大,其值分別為0.043 7d-1,271.93mL/g,71.59%。發(fā)酵初始(第1日)VFAs中乙酸濃度與牛糞比例成正比,發(fā)酵中期(第5日)丙酸積累濃度與秸稈比例成正比。發(fā)酵周期內(nèi)糞稈配比與限速階段的關(guān)系:第1日,秸稈比例越大,產(chǎn)氣限制階段越傾向于水解酸化階段,第2-15日,秸稈比例越大,產(chǎn)氣限速階段越傾向于產(chǎn)氫產(chǎn)乙酸階段;第16-30日,各組發(fā)酵產(chǎn)氣限制階段均為水解酸化階段。該試驗重點對糞稈配比與產(chǎn)氣限速因素的關(guān)系進行研究,為進一步提高混合厭氧發(fā)酵產(chǎn)氣效率提供了理論和試驗依據(jù)。

    生物能;混合厭氧發(fā)酵;水解;揮發(fā)性脂肪酸;甲烷;限速階段

    中國是農(nóng)業(yè)大國,同時也是世界上農(nóng)業(yè)廢棄物產(chǎn)生量較大的國家之一[1],其農(nóng)業(yè)廢棄物的產(chǎn)生量每年約有40多億t,其中農(nóng)作物秸稈7.0億t(折合標(biāo)準(zhǔn)煤量3.53億t),畜禽糞便排放量26.1億t,其余的農(nóng)業(yè)垃圾約7.5億t[2]。但目前中國農(nóng)村大部分的農(nóng)業(yè)廢棄物未經(jīng)合理處置便被直接排入環(huán)境,如牲畜糞便隨意排放、大面積秸稈直接焚燒等,對環(huán)境造成嚴重的污染和破壞[3]。因此,綜合利用農(nóng)業(yè)廢棄物資源對于節(jié)約資源、保護環(huán)境、增加農(nóng)民收入、促進循環(huán)農(nóng)業(yè)的發(fā)展都具有重要的現(xiàn)實意義[4]。

    厭氧發(fā)酵生物技術(shù)可以將農(nóng)作物秸稈、牲畜糞便、市政廢水等有機廢棄物作為原料生產(chǎn)沼氣,這在許多國家都有廣泛的應(yīng)用[5-6]。許多研究發(fā)現(xiàn)混合發(fā)酵比單一發(fā)酵有更高的基質(zhì)利用率,能夠增加80%~400%的產(chǎn)氣量[7-8],且混合發(fā)酵可以避免單一發(fā)酵的底物利用率低、易酸化等問題。目前,對于不同原料為底物的混合厭氧發(fā)酵已有較多的研究,如徐鑫等[9]利用響應(yīng)曲面法研究得到,在39 ℃條件下,原料配比(秸稈∶廚余)為1.29∶1,碳氮比(C/N)為20.42,總固體(TS)質(zhì)量分數(shù)為10.07%時,CH4的總產(chǎn)氣量達到最高值。李偉等[10]的研究結(jié)果表明,在35 ℃條件下,雞糞與玉米秸稈、豬糞與玉米秸稈均以2∶1配比的產(chǎn)氣效果最好。但這些研究大多集中在優(yōu)化配比及發(fā)酵產(chǎn)氣特性的研究上,而對不同原料配比條件下產(chǎn)氣性能差異的原因研究得不夠深入,對發(fā)酵過程限速階段的研究也較少。

    因此,本文通過對厭氧發(fā)酵中水解酸化、產(chǎn)氫產(chǎn)乙酸、產(chǎn)甲烷過程的關(guān)鍵參數(shù)進行監(jiān)測,從底物中微生物的引入、C/N、發(fā)酵基質(zhì)、理論產(chǎn)氣量等角度討論了牛糞與玉米秸稈的不同配比(以下簡稱糞稈配比,以揮發(fā)性固體(VS)計算,下同)對發(fā)酵特性的影響,初步嘗試尋找糞稈配比與發(fā)酵產(chǎn)氣限速因素之間的關(guān)系,以期在獲得最優(yōu)糞稈配比的基礎(chǔ)上,確定產(chǎn)氣限速因素,為進一步解決產(chǎn)氣限速問題、提高混合厭氧發(fā)酵的產(chǎn)氣效率提供支撐。

    1 材料與方法

    1.1試驗材料

    發(fā)酵原料為玉米秸稈、牛糞,分別取自河北科技大學(xué)新校區(qū)附近的農(nóng)田和某養(yǎng)牛場。自然風(fēng)干后的玉米秸稈使用微型植物粉碎機粉碎至30目(約0.613mm)后備用。接種物取自河北省某淀粉廢水處理廠厭氧反應(yīng)器中的顆粒污泥。原料、接種物基本特性見表1。

    1.2試驗裝置

    厭氧發(fā)酵試驗裝置由150mL錐形瓶和250mL量筒組成,分別作為發(fā)酵罐和集氣裝置。試驗裝置如圖1所示。

    表1 原料和接種物的基本特性

    注:表中TS表示總固體;VS表示揮發(fā)性固體。

    圖1 發(fā)酵試驗裝置示意圖Fig.1 Schematic of experimental equipment

    1.3試驗設(shè)計與方法

    將牛糞與粉碎過篩后的玉米秸稈按不同配比組合試驗,共7組,試驗設(shè)計如表2所示。按照表2中的配比將牛糞、玉米秸稈、接種物、水加入150mL錐形瓶后,塞好橡膠塞,再用封口膜密封,置于恒溫水浴鍋中(溫度為(38±1)℃)進行發(fā)酵,量筒倒置于飽和NaCl水槽中。試驗周期為30天。

    試驗每日測量產(chǎn)氣量與沼氣成分。在發(fā)酵的第1,2,5,10,16,30日使用橡膠塞上配備的取樣管取出沼液約3mL,測定揮發(fā)性脂肪酸(VFAs)和pH值。試驗設(shè)置2組平行樣。

    表2 試驗設(shè)計

    注:表中w(牛糞)和w(玉米秸稈)均以揮發(fā)性固體(VS)計算。

    1.4 測試項目與方法

    采用向下排水集氣法測定日沼氣產(chǎn)量;CH4和CO2的濃度采用氣相色譜儀(GC-7890T,上海天美科學(xué)儀器有限公司提供)測定,分析條件:載氣為氬氣,流量為21mL/min,柱溫為40 ℃,進樣口溫度為110 ℃,TCD檢測器溫度為140 ℃,進樣量為0.1mL;采用精密pH值測定儀(DELTA-320)測定pH值;VFAs采用氣相色譜儀(GC-7900,上海天美科學(xué)儀器有限公司提供)測定[11],分析條件:FID檢測器溫度為240 ℃,進樣器溫度為150 ℃,分流比為10∶1,柱流量為1.1mL/min,進樣量為2μL;沼液中TC采用總有機碳分析儀測定;采用凱氏定氮法測量總氮;采用重鉻酸鉀法測量總有機碳(TOC);采用烘干法測量原料中總固體(TS)和揮發(fā)性固體(VS)的含量[12]。

    1.5數(shù)據(jù)處理

    本文對未溶解有機碳的水解過程采用一級動力學(xué)方程進行擬合[13],即

    (1)

    式(1)中:S0為未溶解有機碳的濃度,mol/L;S為t時刻未溶解有機碳的濃度,mol/L;Kd為不同糞稈配比中未溶解有機碳的水解速率常數(shù),d-1。式(1)為對顆粒性有機物水解過程有較好擬合效果的一級動力學(xué)模型,在2.1中對混合發(fā)酵過程中的未溶解性有機碳數(shù)據(jù)進行擬合。

    發(fā)酵階段VFAs中游離態(tài)、分子態(tài)單酸的濃度通過式(2)計算得到[14]:

    (2)

    式(2)中:[A-]表示解離態(tài)單酸的濃度;[HA]表示分子態(tài)單酸的濃度;Ka表示各種單酸的解離常數(shù);H+表示氫離子濃度;本試驗測定的VFAs為[HA]和[A-]的總濃度。式(2)用于2.4中未解離態(tài)丙酸濃度的計算。

    采用Modified Gompertz方程[15],對厭氧發(fā)酵生物化學(xué)甲烷勢進行測定:

    (3)

    式(3)中:M(t)表示時間t時的累積產(chǎn)甲烷量,mL;Pm表示最大累積產(chǎn)甲烷量,mL;Rm表示最大產(chǎn)甲烷速率,mL/d;λ表示滯留時間。將擬合后的Pm和Rm分別除以初始投加的VS質(zhì)量,得到單位VS最大產(chǎn)甲烷量Ps(mL/g)以及最大比產(chǎn)甲烷速率Rs[mL/(g·d)]的數(shù)值。通過分析Ps和Rs等指標(biāo)來比較甲烷生物轉(zhuǎn)化效率。式(3)在2.5中作為產(chǎn)氣性能模型對產(chǎn)甲烷量數(shù)據(jù)進行擬合。

    理論產(chǎn)氣量可以通過化學(xué)計量學(xué)的方法根據(jù)經(jīng)驗的化學(xué)方程式來估算,若有機質(zhì)以分子式CaHbOcNd表示,并假設(shè)它完全降解轉(zhuǎn)換為CO2和CH4,則有機質(zhì)完全降解過程可用式(4)表示[16]:

    (4)

    從而理論甲烷產(chǎn)率Yt為[17]

    (5)

    式(4)、式(5)將用于表5中“理論CH4產(chǎn)量值和“理論CO2產(chǎn)量值”的計算。

    2 結(jié)果與討論

    2.1不同糞稈配比對未溶解有機碳水解動力學(xué)參數(shù)的影響

    厭氧發(fā)酵中未溶解有機碳的水解動力學(xué)模型采用對顆粒性有機物水解過程有較好模擬效果的一級動力學(xué)模型(因80%的產(chǎn)氣量在前15日產(chǎn)生,所以對各組試驗中未溶解有機碳前15日、后15日(即第16—30日)的數(shù)據(jù)分別進行一級動力學(xué)擬合),所得動力學(xué)參數(shù)如表3所示。

    表3 不同糞稈配比混合發(fā)酵未溶解有機碳水解速率常數(shù)

    擬合后所得的線性相關(guān)系數(shù)R2皆大于0.97,表明該線性模型可以很好地反映不同糞稈配比厭氧發(fā)酵時未溶解有機碳的實際水解過程。發(fā)酵過程中,在發(fā)酵前15日內(nèi),隨玉米秸稈比例的增加,未溶解有機碳的水解速率常數(shù)(Kd)從純糞發(fā)酵的0.032 3d-1增加至糞稈配比為2∶1發(fā)酵的0.043 7d-1,而后降低至純秸稈發(fā)酵的0.030 2d-1,水解速率常數(shù)在0.030 2~0.043 7d-1范圍內(nèi)變化,與趙云飛等[18]研究所得水解速率常數(shù)(0.035~0.094d-1)的范圍相近。在發(fā)酵后15日內(nèi),未溶解有機碳的水解速率常數(shù)在0.004 4~0.009 0d-1范圍內(nèi)變化。

    在發(fā)酵前15日內(nèi),隨秸稈比例增加,Kd呈先增加后降低的趨勢,糞稈配比為2∶1時Kd最大為0.043 7d-1(見表3)。推測其原因,一方面可能由于適合的糞稈配比為水解微生物提供了適宜生存的C/N(糞稈配比為2∶1時,原料C/N為20.23,見表2),從而促進了水解過程。另一方面,一定比例牛糞的添加既可以引入一定量瘤胃中厭氧纖維素分解菌促進水解,而牛糞中的碳酸氫根、氨氮又可以增加發(fā)酵體系的堿度,從而緩解水解產(chǎn)酸過程導(dǎo)致的pH值不穩(wěn)定。另外,后15日純秸稈發(fā)酵的Kd為各組最大值(0.009 0 d-1),表明秸稈單獨發(fā)酵的水解過程更慢、周期更長。

    前15日內(nèi)比后15日內(nèi)的未溶解有機碳水解速率常數(shù)大10倍左右(見表3),這與前15日內(nèi),VFAs出現(xiàn)積累(見圖2)、各組累計產(chǎn)氣量能達到總量的80%以上(見表4),而后15日內(nèi)VFAs、日產(chǎn)氣量大大降低的現(xiàn)象相一致。有研究表明[19],固體有機物(特別是如麥稈和玉米秸稈等高纖維含量的有機質(zhì))的水解產(chǎn)酸速率低于產(chǎn)甲烷速率,是固體有機物厭氧發(fā)酵的限速步驟。本研究結(jié)果表明,當(dāng)以牛糞和秸稈作為基質(zhì)時,可能由于秸稈中的部分半纖維素、纖維素和木質(zhì)素在短時間內(nèi)難以被水解菌群代謝成有機酸,從而導(dǎo)致可被產(chǎn)甲烷菌群代謝的底物不足,因而后15日較慢的水解過程導(dǎo)致的產(chǎn)氣限速現(xiàn)象較前15日更為嚴重。

    表4 不同糞稈配比混合厭氧發(fā)酵基于Gompertz模型產(chǎn)氣性能比較

    2.2不同糞稈配比對發(fā)酵液pH值的影響

    圖2 糞稈配比對發(fā)酵液中pH值的影響Fig.2 Effect of various ratios of cow manure and corn straw on pH value in anaerobic fermentation

    不同糞稈配比的發(fā)酵液pH值隨發(fā)酵時間推移,而經(jīng)歷了先上升后降低,最后又緩慢升高的過程(見圖2)。在發(fā)酵液配好的第1日,不同糞稈配比發(fā)酵液的pH值為7.6~7.8,糞稈配比對發(fā)酵起始時的pH值影響不大;但在第2日,糞稈配比中糞的比例越大,其發(fā)酵液的pH值越高;且在隨后的第5,10,16,30日,發(fā)酵液pH值的變化趨勢也遵循著此基本規(guī)律。

    圖3 糞稈配比對發(fā)酵沼液中VFAs及COD質(zhì)量濃度的影響Fig.3 Effect of various ratios of cow manure and corn straw on VFAs in anaerobic fermentation

    隨發(fā)酵時間的推移,發(fā)酵液pH值的變化規(guī)律在一定程度上可以響應(yīng)發(fā)酵液中VFAs濃度的變化趨勢,二者變化趨勢大致相反。發(fā)酵液pH值升高時,VFAs濃度在降低,第5日時,尤為明顯,發(fā)酵液pH值降低時,VFAs濃度在升高。發(fā)酵液pH值隨時間的變化趨勢會受發(fā)酵液中的酸、堿、氨氮和碳酸等物質(zhì)的濃度綜合影響,因此不像VFAs值一樣,會因糞稈配比的不同呈現(xiàn)出明顯的規(guī)律變化性。

    2.3不同糞稈配比對沼液中VFAs濃度的影響

    糞稈配比對發(fā)酵沼液中乙酸、丙酸和丁酸的濃度變化趨勢影響較大(如圖3所示)。當(dāng)w(牛糞)∶w(玉米秸稈)(以下簡化為w(糞)∶w(稈))≥2∶1(即糞稈配比為純糞,6∶1,4∶1,2∶1)時,VFAs濃度(以COD計)自高峰值先迅速下降,而后逐漸升高,隨后隨著發(fā)酵時間的延長逐漸下降。當(dāng)w(糞)∶w(稈)<2∶1(即糞稈配比為1∶2,1∶4,純秸稈)時,VFAs呈現(xiàn)自較低值逐步上升,隨后隨著發(fā)酵時間的延長逐漸下降的趨勢(見圖3d))。w(糞)∶w(稈)≥2∶1的各組VFAs(乙酸為主)峰值出現(xiàn)在第1日,質(zhì)量濃度分別為1 466.72,1 233.44,1 132.93,946.20mg/L(見圖3a))。w(糞)∶w(稈)<2∶1的各組VFAs(丙酸為主)峰值出現(xiàn)在第5日,質(zhì)量濃度分別為993.85,1 260.32,2 262.60mg/L(見圖3c))。

    第1日的VFAs濃度與糞添加量成正比(見圖3d)),是因為牛糞含有大量以乙酸、丙酸為主的VFAs。在整個發(fā)酵周期內(nèi),秸稈比例越大,VFAs濃度越高(第5日達到最大值),丙酸的積累越嚴重,丙酸完全降解所需時間越長(見圖3c))。丙酸積累程度隨秸稈比例的增加而增加,可能有以下幾個原因:1)糞稈配比不同導(dǎo)致C/N不同,從而影響發(fā)酵體系中微生物體內(nèi)的生物氧化過程,并進一步導(dǎo)致NADH/NAD+(NADH為還原態(tài)煙酰胺腺嘌呤二核苷酸;NAD+為氧化態(tài)煙酰胺腺嘌呤二核苷酸)的值、發(fā)酵產(chǎn)物產(chǎn)量、酸化功能菌群的變化,而丙酸型發(fā)酵代謝途徑有利于維持NADH/NAD+的穩(wěn)定,因此C/N在某一范圍內(nèi)變化時,微生物可能更傾向于選擇有利于維持系統(tǒng)穩(wěn)定的丙酸型發(fā)酵[20-22]。2)一定量牛糞的添加可引入代謝丙酸的產(chǎn)氫產(chǎn)乙酸菌,可強化乙酸化過程,從而減輕了丙酸的積累程度。瘤胃微生物菌群可將丙酸降解成乙酸直至分解成為沼氣[23]。3)纖維素含量較高的厭氧發(fā)酵過程較容易出現(xiàn)丙酸型發(fā)酵,因此,推測丙酸積累可能與玉米秸稈的纖維結(jié)構(gòu)有關(guān),本研究與陳廣銀等[24]的研究結(jié)論相一致:稻草與豬糞混合發(fā)酵中稻草比例愈高,發(fā)酵液中丙酸含量愈高。

    圖4 糞稈配比對厭氧發(fā)酵日產(chǎn)氣量的影響Fig.4 Effect of various ratios of cow manure and corn straw on biogas production rate in anaerobic fermentation

    VFAs的濃度變化、積累情況,均可體現(xiàn)水解酸化、產(chǎn)氫產(chǎn)乙酸、產(chǎn)甲烷各階段的底物轉(zhuǎn)化速率情況,再結(jié)合產(chǎn)氣等情況,可以初步判斷整個發(fā)酵過程限制產(chǎn)氣速率的階段。第1日時,隨著糞比例的增加,乙酸濃度、日產(chǎn)氣量和CH4濃度均增加(見圖3a)、圖4和圖5a)),可以推測此時,糞比例越小,發(fā)酵產(chǎn)氣限制階段越傾向于水解酸化階段。第2—5日時,隨秸稈比例的增加,VFAs的消耗量小于產(chǎn)生量,發(fā)生了以丙酸為主的酸積累(見圖3c)、圖3d)),第6—15日時,VFAs的消耗量逐漸大于產(chǎn)量,VFAs濃度緩慢下降(見圖3d)),因此可以推測第2—15日時,秸稈比例越大,發(fā)酵產(chǎn)氣限速階段越傾向于產(chǎn)氫產(chǎn)乙酸階段。但值得注意的是,第2—15日時純糞發(fā)酵過程VFAs的濃度較低,未發(fā)生積累(見圖3d)),可見其限速階段應(yīng)該為水解階段。第16—30日時,VFAs的濃度都保持較低水平(見圖3d)),此時最大日均產(chǎn)氣量僅為48.56mL(見圖4),水解速率常數(shù)僅為第1—15日的1/10左右(見表3),可知這時各組發(fā)酵產(chǎn)氣限制階段均為水解酸化階段。

    2.4不同糞稈配比對厭氧發(fā)酵產(chǎn)氣特性的影響

    糞稈配比對厭氧發(fā)酵日產(chǎn)氣量變化趨勢影響較大(如圖4所示)。w(糞)∶w(稈)≥2∶1時,日產(chǎn)氣量在第1日達到最大峰值,因為糞添加比例越大VFAs(乙酸為主)含量越高,產(chǎn)甲烷菌可直接利用的基質(zhì)則越多,利于產(chǎn)氣,之后日產(chǎn)氣量波動下降至第5日較低水平。隨后,隨著VFAs利用速率大于產(chǎn)生速率,日產(chǎn)氣量逐漸上升至第8日的高峰,之后緩慢降低,并在第12日累計產(chǎn)氣量達到總產(chǎn)氣量的80%。當(dāng)w(糞)∶w(稈)<2∶1時,發(fā)酵開始后,在第3日才達到日產(chǎn)氣最大值,推測是因為秸稈比例越大初始VFAs越低,初始沒有足夠的甲烷菌可利用的基質(zhì),水解酸化階段就更容易成為VFAs產(chǎn)生的限制階段,從而成為產(chǎn)氣的限制階段。而第3日時的VFAs(丙酸為主)有一定的積累,使得產(chǎn)氣突增。日產(chǎn)氣量自第3日開始下降直至第7日,隨后逐漸上升第10日的高峰,之后緩慢降低,在第15日累計產(chǎn)氣量達到總產(chǎn)氣量的80%。

    隨秸稈比例的增加,第5日左右的VFAs(丙酸為主)濃度保持在較高水平(如圖3所示),這可能是引起第7日w(糞)∶w(稈)<2∶1的試驗組出現(xiàn)產(chǎn)氣量極低值的原因。有文獻報道丙酸質(zhì)量濃度大于900 mg/L時,發(fā)酵過程即出現(xiàn)嚴重的酸抑制[25],但發(fā)酵中酸抑制的情況還需要通過對解離態(tài)酸和未解離態(tài)酸的濃度大小來判定。對厭氧微生物菌群,特別是產(chǎn)甲烷菌群活性造成明顯抑制作用是未解離態(tài)脂肪酸的濃度[26]。DUARTER等[27]發(fā)現(xiàn)未解離的揮發(fā)性脂肪酸濃度達到0.17~0.42 mmol/L時,厭氧反應(yīng)過程受到抑制。KROEKER等[28]提出未解離的揮發(fā)脂肪酸濃度達到0.50~1.00 mmol/L時,產(chǎn)甲烷菌活性完全喪失。本試驗中,以丙酸為主的酸積累與產(chǎn)氣降低的相關(guān)性有2種可能:其一,由各組沼液此時的pH值計算得出,各組中未解離態(tài)丙酸最大濃度為0.074 mmol/L(純秸稈發(fā)酵),未達到抑制的最低限值0.17 mmol/L,并且純秸稈發(fā)酵體系第5日丙酸質(zhì)量濃度為3 677.47 mg/L(最大值)時(見圖3),依然出現(xiàn)日產(chǎn)氣高峰279 mL,可說明第7日日產(chǎn)氣量下降的主要原因不是丙酸抑制甲烷菌群作用,而是僅有少量乙酸供甲烷菌群代謝,同時,雖然發(fā)酵液中有大量丙酸,但其降解成乙酸的“產(chǎn)氫產(chǎn)乙酸過程”速率較低所致[29-30]。其二,由于試驗未在第7日取沼液分析VFAs,因此也有可能w(糞)∶w(稈)<2∶1的試驗組在第7日產(chǎn)生了丙酸積累高峰,并超過了酸抑制的最低限值,發(fā)生了丙酸抑制產(chǎn)甲烷菌活性的情況。

    圖5 糞稈配比對發(fā)酵沼氣中CH4,CO2質(zhì)量分數(shù)的影響Fig.5 Effect of various ratios of cow manure and corn straw on CH4 and CO2 in methane

    糞稈配比不同,發(fā)酵過程沼氣中CH4濃度變化趨勢也不同(見圖5a))。前6日各組CH4濃度均呈逐漸上升趨勢,糞添加比例由高到低的發(fā)酵體系的日均CH4質(zhì)量分數(shù)(前6日)分別為56.45%,55.76%,54.53%,53.98%,49.77%,45.64%,38.23%。可以看出隨牛糞比例增大,CH4濃度越大。此后直至第20日,受VFAs濃度、組成和甲烷菌活性等因素影響,各組CH4質(zhì)量分數(shù)在45%~67%之間波動。第20—30日,CH4濃度基本隨著秸稈比例增大而增加。發(fā)酵初期,隨著糞比例增加CH4濃度也較高的原因:一方面,初期VFAs(乙酸為主)濃度隨著糞比例增大而增加,甲烷菌相對充足、可直接利用乙酸,另一方面,是因為牛糞里瘤胃微生物中大量的產(chǎn)甲烷菌可以強化甲烷化過程。值得注意的是發(fā)酵第16—30日時,雖然VFAs濃度較低,日產(chǎn)氣量隨之減少,但對CH4濃度影響不大,說明發(fā)酵過程中甲烷菌活性較穩(wěn)定,可能是因為接種物比例較大(占總體積53%),菌群活性較高且接種物為發(fā)酵前后顆粒形態(tài)未發(fā)生變化的顆粒污泥。

    糞稈配比不同,發(fā)酵過程沼氣中CO2濃度變化趨勢不同(見圖5b))。各組發(fā)酵體系從第1日起CO2濃度逐漸上升至第3日的高峰,糞添加比例由高到低的發(fā)酵體系CO2質(zhì)量分數(shù)峰值分別為33.82%,37.61%,38.33%,42.79%,45.08%,47.10%,55.39%??梢钥闯?,秸稈比例越大CO2濃度越大。第3—20日CO2濃度受水解酸化、甲烷化程度的影響逐步波動下降。第20—30日,糞稈配比由高到低的發(fā)酵體系日均CO2質(zhì)量分數(shù)分別為11.70%,13.18%,13.94%,15.24%,14.56%,16.67%,19.28%?;旧显谡麄€發(fā)酵周期內(nèi),秸稈比例越大CO2濃度越高,可能有2個方面原因:一方面,秸稈比例較大的發(fā)酵體系,其理論產(chǎn)CO2量較大。另一方面,秸稈比例越大,初期未溶解有機質(zhì)濃度較高,經(jīng)歷水解階段亦能貢獻15%左右的CO2。此外,有研究表明,CO2的濃度隨C/N的增加而增加。第20—30日各組發(fā)酵體系的CO2濃度降低到20%以下,可能是因為同型產(chǎn)乙酸菌不斷消耗CO2產(chǎn)出CH4所致。此外,各組的理論CH4與CO2產(chǎn)量相差不大(見表5),但CH4和CO2實際質(zhì)量分數(shù)的變化范圍分別為50%~70%,10%~50%,可見實際中CO2的濃度遠小于CH4,主要是因為CO2的溶解度遠高于CH4(常壓下,40 ℃,CO2和CH4在水中的溶解度分別[31]為0.53mL/mL,0.023 7mL/mL),因此發(fā)酵液中溶解了一定量的CO2,使得實際檢測到的CO2濃度偏低。

    2.5不同糞稈配比對厭氧發(fā)酵單位VS累計產(chǎn)氣特性的影響

    采用CurveExpert軟件中的ModifiedGompertz模型對不同糞稈配比單位VS累計產(chǎn)CH4量的數(shù)據(jù)分別進行擬合,獲得相關(guān)參數(shù)(見表4)。

    所得相關(guān)系數(shù)R2值皆大于0.99,表明該模型可以很好地反映不同糞稈配比發(fā)酵產(chǎn)CH4的過程。最大累計產(chǎn)CH4量(Ps)、最大產(chǎn)CH4速率(Rs)隨秸稈比例的增加,呈現(xiàn)先升高后降低趨勢,當(dāng)w(糞)∶w(稈)=2∶1時,Ps和Rs獲得最大值。停滯時間(λ)隨秸稈比例的增加,基本呈現(xiàn)先降后升的趨勢,w(糞)∶w(稈)=4∶1時,λ獲得最小值0.010天,純糞發(fā)酵時,λ獲得最大值1.71天。當(dāng)w(糞)∶w(稈)≥2∶1時,各組發(fā)酵周期(即85%Ps)相差不大(最大相差0.91天),當(dāng)w(糞)∶w(稈)<2∶1時,發(fā)酵周期隨秸稈比例增大而延長。

    上述結(jié)果表明,糞稈過高、過低都不利于發(fā)酵,當(dāng)w(糞)∶w(稈)=2∶1時,Ps,Rs,λ和發(fā)酵周期分別為271.93 mL/g,21.71 mL/(g·d),0.094天,13.58天。綜合4個參數(shù)考慮,確定該比例為發(fā)酵最佳糞稈配比,與牛糞或玉米秸稈單獨發(fā)酵相比,單位VS累計CH4量分別提高了21.45%和18.21%。

    2.6不同糞稈配比厭氧發(fā)酵產(chǎn)氣性能比較

    厭氧發(fā)酵產(chǎn)氣的性能與發(fā)酵原料組成成分、形態(tài)等因素有關(guān)。根據(jù)化學(xué)計量學(xué)的方法,結(jié)合牛糞、玉米秸稈各元素含量,計算各試驗配比下理論CH4與CO2的產(chǎn)量(見表5)。

    從表5中可以看出,每克VS發(fā)酵原料的CH4和CO2理論產(chǎn)量都隨秸稈比例的增加而增加。每克VS發(fā)酵原料的實際CH4產(chǎn)量、實際CO2產(chǎn)量、生物轉(zhuǎn)化產(chǎn)CH4效率、生物轉(zhuǎn)化產(chǎn)CO2效率,都隨著秸稈比例的增加呈現(xiàn)先增加后降低的趨勢,最大值在w(糞)∶w(稈)=2∶1時出現(xiàn),分別為274.16mL/g,160.04mL/g,71.59%,40.75%,這說明牛糞與玉米秸稈一定比例的混合發(fā)酵比單一發(fā)酵有更高的生物轉(zhuǎn)化效率。

    研究發(fā)現(xiàn),w(糞)∶w(稈)=2∶1發(fā)酵體系的理論產(chǎn)CH4,CO2量都小于秸稈單獨發(fā)酵,但其實際發(fā)酵產(chǎn)CH4,CO2量都大于純秸稈發(fā)酵,可能原因有以下幾點:1)纖維素本身的結(jié)晶結(jié)構(gòu)及木質(zhì)素的物理屏障作用,使得以木質(zhì)纖維素為主要成分的玉米秸稈難以被降解,雖然牛糞主要成分也為木質(zhì)纖維素,但牛糞中含有大量瘤胃微生物,主要包括水解細菌、產(chǎn)氫產(chǎn)乙酸菌群和產(chǎn)甲烷菌群等能夠分泌的各種酶類的菌群[32]。因此牛糞也充當(dāng)了“接種物”的角色,其添加強化了水解、酸化、甲烷化等階段,使得木質(zhì)纖維素降解能夠高效進行。2)根據(jù)本試驗的研究,牛糞的添加可以減輕丙酸積累的情況。3)牛糞的添加調(diào)整了發(fā)酵體系的C/N及堿度,更適于微生物發(fā)酵。此外,牛糞中的成分相對于玉米秸稈更加復(fù)雜,能夠提供發(fā)酵過程中所需的微量元素等物質(zhì)。

    表5 不同糞稈配比厭氧發(fā)酵產(chǎn)氣性能比較

    3 結(jié) 論

    1)隨玉米秸稈比例的增加,未溶解有機碳的水解速率常數(shù)(Kd)呈先增后減的規(guī)律,w(糞)∶w(稈)=2∶1時Kd得到最大值;初始VFAs(乙酸為主)的濃度與牛糞添加比例成正比;丙酸積累(第5日左右)的程度與玉米秸稈添加比例成正比。

    2)結(jié)合水解速率常數(shù)、VFAs、日產(chǎn)氣量變化情況來判斷發(fā)酵周期內(nèi)糞稈配比與限速階段的關(guān)系:第1日時,發(fā)酵產(chǎn)氣的限制階段越傾向于水解酸化階段。第2—15日時,發(fā)酵產(chǎn)氣限速階段越傾向于產(chǎn)氫產(chǎn)乙酸階段。第16—30日時,各組發(fā)酵產(chǎn)氣限制階段均為水解酸化階段。

    3)整個發(fā)酵周期內(nèi),秸稈比例越大CO2濃度越高;當(dāng)w(糞)∶w(稈)=2∶1時,水解速率常數(shù)(前15日)、最大產(chǎn)甲烷量、實際生物轉(zhuǎn)化產(chǎn)CH4效率分別為0.043 7d-1,271.93mL/g,71.59%。

    /< class="emphasis_italic">References

    :

    [1]阮建雯,蔡宗濤,張霞.云南省農(nóng)業(yè)廢棄物資源化利用狀況及對策[J].中國沼氣,2008,26(2):46-51.

    RUANJianwen,CAIZongtao,ZHANGXia.ResourceutilizationofagriculturalresiduesandstrategiesinYunnanProvince[J].ChinaBiogas, 2008, 26(2): 46-51.

    [2]張廣裕.農(nóng)業(yè)廢棄物處理及環(huán)境治理研究[J].甘肅農(nóng)業(yè),2012(23):52-53.

    ZHANGGuangyu.Wastedisposalandenvironmentalmanagementofagricultural[J].JournalofGansuAgriculturalUniversity, 2012(23): 52-53.

    [3]許智,葉小梅,常志洲,等.稻秸、餐廚垃圾及人糞尿混合厭氧發(fā)酵[J].環(huán)境工程學(xué)報,2012,6(7):2447-2453.

    XUZhi,YEXiaomei,CHANGZhizhou,etal.Co-digestionofricestraw,kitchenwasteandtoiletwasteforbiogasproduction[J].ChineseJournalofEnvironmentalEngineering, 2012, 6(7): 2447-2453.

    [4]曹穩(wěn)根,高貴珍,方雪梅,等.我國農(nóng)作物秸稈資源及其利用現(xiàn)狀[J].宿州學(xué)院學(xué)報,2007,22(6):110-113.

    CAOWen’gen,GAOGuizhen,FANGXuemei,etal.Resourcesandutilizationofcropstalksstatus[J].JournalofSuzhouCollege, 2007, 22(6): 110-113.

    [5]孟祥林.循環(huán)農(nóng)業(yè)與“循環(huán)鏈”構(gòu)建:經(jīng)典案例、存在問題與發(fā)展選擇[J].河北科技大學(xué)學(xué)報(社會科學(xué)版),2014,14(2):15-21.

    MENGXianglin.Constructionofcircularagricultureand"circularchain":Classiccase,problemsanddevelopmentchoice[J].JournalofHebeiUniversityofScienceandTechnology(SocialSciences), 2014, 14(2): 15-21.

    [6]HAMEDM,EL-MASHAD,ZHANGRuihong.Biogasproductionfromco-digestionofdairymanureandfoodwaste[J].BioresourceTechnologyBioresourceTechnology, 2010, 101(11): 4021-4028.

    [7]LUOQingming,LIXiujin,ZHUBaoning.AnaerobicbiogasificationofNaOH-treatedcornstalk[J].TransactionsoftheCSAE, 2005, 21(2):111-115.

    [8]WEILANDP.AnaerobicwastedigestioninGermany-statusandrecentdevelopment[J].Biodegradation, 2004, 11(6): 415-421.

    [9]徐鑫,寇巍,張大雷,等.以秸稈和廚余垃圾為原料發(fā)酵產(chǎn)甲烷工藝的優(yōu)化[J].可再生能源,2010,28(6):76-83.

    XUXin,KOUWei,ZHANGDalei,etal.Technologicaloptimizationonstalkandkitchenresiduesfermentationformethaneproduction[J].RenewableEnergyResources, 2010, 28(6): 76-83.

    [10]李偉,李文靜,張彤,等.不同配比雞糞、豬糞與玉米稈混合發(fā)酵的產(chǎn)氣效果[J].西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版),2011,39(5):79-84.

    LIWei,LIWenjing,ZHANGTong,etal.Effectofdifferentmixturesonanaerobicfermentationofdungandcornstraw[J].JournalofNorthwestA&FUniversity(NaturalScienceEdition), 2011, 39(5):79-84.

    [11]羅毅,李曉麗,劉暉,等. 超聲波、堿、超聲波/堿預(yù)處理污泥對污泥水解效果的影響[J].河北工業(yè)科技,2015,32(6):552-556.

    LUOYi,LIXiaoli,LIUHui,etal.Effcetofultrasonic,alkaliandultrasonic/alkalipretreatmentonsludgehydrolysis[J].HebeiJournalofIndustrialScienceandTechnology, 2015, 32(6): 552-556.

    [12]中國科學(xué)院成都生物研究所《沼氣發(fā)酵常規(guī)分析》編寫組.沼氣發(fā)酵常規(guī)分析[M].北京:北京科學(xué)技術(shù)出版社,1984:67-80.

    [13]SOSNOWSKIP,KLEPACZ-SMOLKAA,KACEOREKK,etal.Kineticinvestigationsofmethaneco-fermentationofsewagesludgeandorganicfractionofmunicipalsolidwastes[J].BioresourceTechnology, 2008, 99(13): 5731-5737.

    [14]GINKELSV,LOGANBE.Inhibitionofbiohydrogenproductionbyundissociatedaceticandbutyricacids[J].EnvironmentalScience&Technology, 2005, 39(23): 9351-9356.

    [15]ZWIETERINGMH,JONGENBURGERI,ROMBOUTSFM,etal.Modelingofthebacterialgrowthcurve[J].ApplicationofEnvironmentalMicrobiology,1990, 56(6): 1875-1881.

    [16]HEON,PARKS,KANGH.Effectsofmixtureratioandhydraulicretentiontimeonsingle-stageanaerobicco-digestionoffoodwasteandwasteactivatedsludge[J].JournalofEnvironmentalScienceandHealth, 2004, 39(7): 1739-1756.

    [17]MLLERH,SOMMERS,AHRINGB.Methaneproductivityofmanure,strawandsolidfractionsofmanure[J].BiomassandBioenergy, 2004, 26(5): 485-495.

    [18]趙云飛,劉曉玲,李十中,等.有機成分比例對高固體濃度厭氧發(fā)酵產(chǎn)甲烷的影響[J].中國環(huán)境科學(xué),2012,32(6):1110-1117.

    ZHAOYunfei,LIUXiaoling,LIShizhong,etal.Effectsoforganicsubstancemixingratiosonmethanebioconversionthroughhigh-solidsanaerobicco-fermentation[J].ChinaEnvironmentalScience, 2012, 32(6): 1110-1117.

    [19]陳廣銀,曹杰,葉小梅,等.pH值調(diào)控對秸稈兩階段厭氧發(fā)酵產(chǎn)沼氣的影響[J].生態(tài)環(huán)境學(xué)報,2015,24(2):336-342.

    CHENGuangyin,CAOJie,YEXiaomei,etal.EffectofpHvaluescontrolonbiogasproductionofagriculturalstrawwithtwo-phaseanaerobicfermentation[J].EcologyandEnvironmentalSciences, 2015, 24(2): 336-342.

    [20]劉和,劉曉玲,邱堅,等.C/N對污泥厭氧發(fā)酵產(chǎn)酸類型及代謝途徑的影響[J].環(huán)境科學(xué)學(xué)報,2010,30(2):340-346.

    LIUHe,LIUXiaoling,QIUJian,etal.TheeffectsofC/Nratioontheproductionofvolatilefattyacidsandthemetabolicpathwayofanaerobicfermentationofsewagesludge[J].JournalofEnvironmentalSciences, 2010, 30(2): 340-346.

    [21]任南琪,趙丹,陳曉蕾,等.厭氧生物處理丙酸產(chǎn)生和積累的原因及控制對策[J].中國環(huán)境科學(xué),2001,32(1):83-89.

    RENNanqi,ZHAODan,CHENXiaolei,etal.Anaerobicbiologicaltreatmentofacidproducedandaccumulatedcausesandcontrolmeasures[J].ChinaEnvironmentalScience, 2001, 32(1): 83-89.

    [22]劉曉玲.城市污泥厭氧發(fā)酵產(chǎn)酸條件優(yōu)化及其機理研究[D].無錫:江南大學(xué),2008.

    LIUXiaoling.TheConditionOptimizationofSewageSludgeforProducingVolatileFattyAcidsandtheInvestigationofAcidogenicMechanism[D].Wuxi:JiangnanUniversity, 2008.

    [23]郭嫣秋.瘤胃產(chǎn)甲烷菌定量檢測與微生物菌群調(diào)控研究[D].杭州:浙江大學(xué),2008.

    GUOYanqiu.QuantificationofMethanogensandManipulationofMicrobialCommunityintheRumen[D].Hangzhou:ZhejiangUniversity, 2008.

    [24]陳廣銀,鄭正,鄒星星,等.稻草與豬糞不同比例混合的厭氧消化特性[J].環(huán)境化學(xué),2009,28(1):21-25.

    CHENGuangyin,ZHENGZheng,ZOUXingxing,etal.Characteristicsofanaerobicco-digestionofricestrawandswinefecesinvariousratios[J].EnvironmentalChemistry, 2009, 28(1): 21-25.

    [25]WANGYY,ZHANGYL,WANGJB.Effectsofvolatilefattyacidconcentrationsonmethaneyieldandmethanogenicbacteria[J].BiomassandBioenergy, 2009, 33(5): 848-853.

    [26]BOEK.OnlineMonitoringandControloftheBiogasProcess[D].Copenhagen:TechnicalUniversityofDenmark, 2006.

    [27]DUARTERAC,ANDERSONGK.Inhibitionmodelinginanaerobicdigestion[J].WaterScienceandTechnology, 1982, 14(6/7): 749-763.

    [28]KROEKEREJ,SCHULTEDD,SPARLINGAB,etal.Anaerobictreatmentprocessstability[J].WaterPollutionControlFederation, 1979, 51(4): 718-727.

    [29]李紅麗,曹霏霏,王巖.揮發(fā)性脂肪酸對厭氧干式發(fā)酵產(chǎn)甲烷的影響[J].環(huán)境工程學(xué)報,2014,8(6):2572-2578.

    LIHongli,CAOFeifei,WANGYan.Effectofvolatilefattyacidsonmethaneproductionofdryanaerobicdigestion[J].ChineseJournalofEnvironmentalEngineering, 2014, 8(6): 2572-2578.

    [30]裴占江,劉杰,王粟,等.餐廚垃圾與牛糞聯(lián)合厭氧消化效率研究[J].中國沼氣,2014,32(4):3-7.

    PEIZhanjiang,LIUJie,WANGSu,etal.Anaerobicdigestionefficiencyofmixedfoodwastesandcattlemanure[J].ChinaBiogas, 2014, 32(4): 3-7.

    [31]李夢龍.化學(xué)數(shù)據(jù)速查手冊[M].北京:化學(xué)工業(yè)出版社,2003:169-170.

    [32]陳慶今,劉煥彬,胡勇有.瘤胃微生物對纖維素降解機理[J].微生物學(xué)雜志,2002,22(2):44-46.

    CHENQingjin,LIUHuanbin,HUYongyou.Cellulosedegradationmechanismofrumenmicrobes[J].JournalofMicrobiology, 2002, 22(2): 44-46.

    Effectofdifferentratiosofcowmanureandcornstrawonthemixedanaerobicfermentationrate

    JIANGZongshan1, 2,3,PINGJiapeng1,2,GUOYankai1,2,CHENXinghua3,GUOJianbo1,4,LIANJing1,2

    (1.SchoolofEnvironmentalScienceandEngineering,HebeiUniversityofScienceandTechnology,Shijiazhuang,Hebei050018,China; 2.PollutionPreventionBiotechnologyLaboratoryofHebeiProvince,Shijiazhuang,Hebei050018,China; 3.HangzhouEnergy&EnvironmentalEngineeringCompanyLimited,Hangzhou,Zhejiang310000,China; 4.SchoolofEnvironmentalandMunicipalEngineering,TianjinChengjianUniversity,Tianjin300384,China)

    Inordertostudytheeffectofthedifferentratiosontheanaerobicfermentationrateisinvestigated,andtherate-limitingfactorsarepreliminarilydetermined,atmesophilic(38±1)℃condition,withanaerobicgranularsludgeasinoculums,differentratiosofcowmanureandcornstrawareusedassubstrateformixedanaerobicfermentation.Bymeasuringdailybiogasproduction,theconcentrationsofCH4andCO2inthemarshgas,TC,theconcentrationofVFAsandpHvalue,Theresultsshowthatunderthemixtureratioof2∶1,thehydrolysisrateconstants,cumulativebiogasyieldandbiodegradabilityCH4reachtheirhighlimits,whichare0.043 7d-1, 271.93mL/gand71.59%,respectively.Moreover,itisfoundthattheconcentrationofaceticacidisproportionaltotheamountofcowmanureatthebeginning(thefirstday)ofmixedfermentation,andtheconcentrationofpropionicacidisproportionaltotheamountofcornstrawinmediumfermentationstage(thefifthday).Inaddition,rate-limitingstepofbiogasproductionisrelatedtotheratioofcowmanureandcorninfermentationmaterial.Withtheincreasingofcornstrawproportion,onthe1stday,ittendstohydrolysisacidogenesis;fromthe2thdayto15thday,ittendstohydrogen-productionacetogenisis;andfromthe16thdayto30thday,itishydrolysisacidogenesis.Thepaperfocusesontherelationshipbetweentheratioofcowmanureandcornstrawandtherate-limitingstepforbiogasproduction,whichcouldprovideatheoreticalandexperimentalsupportforimprovingtheefficiencyofbiogasproductioninmixedfermentation.

    biomassenergy;mixedanaerobicfermentation;hydrolyzation;volatilefattyacids;methane;rate-limitingstep

    1008-1542(2016)04-0396-10

    10.7535/hbkd.2016yx04013

    2015-12-15;

    2016-04-06;責(zé)任編輯:王海云

    國家自然科學(xué)基金(51208170);河北省高校百名優(yōu)秀創(chuàng)新人才支持計劃(BR2-211)

    姜宗姍(1985—),女(滿族),河北承德人,碩士,主要從事生物質(zhì)能利用方面的研究。

    廉靜副教授。E-mail:ningjingzhiyuan98@163.com

    X712

    A

    姜宗姍,平佳芃,郭延凱,等.牛糞與玉米秸稈配比對混合厭氧發(fā)酵速率的影響[J].河北科技大學(xué)學(xué)報,2016,37(4):396-405.

    JIANGZongshan,PINGJiapeng,GUOYankai,etal.Effectofdifferentratiosofcowmanureandcornstrawonthemixedanaerobicfermentationrate[J].JournalofHebeiUniversityofScienceandTechnology,2016,37(4):396-405.

    高清黄色对白视频在线免费看| 成人漫画全彩无遮挡| 久久久久久久久久久久大奶| 亚洲在久久综合| 亚洲国产精品国产精品| a级毛片黄视频| 日本av手机在线免费观看| 搡女人真爽免费视频火全软件| 男女免费视频国产| av视频免费观看在线观看| 天天躁夜夜躁狠狠躁躁| 一级毛片我不卡| 午夜老司机福利剧场| 高清不卡的av网站| 亚洲综合色网址| 色网站视频免费| 亚洲三级黄色毛片| 久久精品国产亚洲av天美| 国产片特级美女逼逼视频| 亚洲人与动物交配视频| 99久国产av精品国产电影| 亚洲精品久久午夜乱码| 自拍欧美九色日韩亚洲蝌蚪91| 免费高清在线观看视频在线观看| 国产男女超爽视频在线观看| 捣出白浆h1v1| 高清在线视频一区二区三区| 性色av一级| 国产精品久久久久久av不卡| 日韩av免费高清视频| 亚洲国产av新网站| 婷婷色综合大香蕉| 91午夜精品亚洲一区二区三区| 全区人妻精品视频| 性高湖久久久久久久久免费观看| 亚洲国产日韩一区二区| 成年女人在线观看亚洲视频| 中文字幕亚洲精品专区| 国产片特级美女逼逼视频| 9色porny在线观看| 少妇的丰满在线观看| 亚洲精品乱码久久久久久按摩| 99国产精品免费福利视频| 亚洲欧美色中文字幕在线| 性色av一级| 99久国产av精品国产电影| 久久综合国产亚洲精品| 亚洲伊人久久精品综合| 久久久亚洲精品成人影院| 美女大奶头黄色视频| 18在线观看网站| 久久99一区二区三区| 国产一区二区激情短视频 | a级毛片在线看网站| 欧美激情国产日韩精品一区| 国产xxxxx性猛交| 久久国产精品男人的天堂亚洲 | 最新的欧美精品一区二区| 午夜福利在线观看免费完整高清在| 高清视频免费观看一区二区| 亚洲欧美色中文字幕在线| 十八禁高潮呻吟视频| 免费人成在线观看视频色| 免费大片黄手机在线观看| 亚洲av日韩在线播放| 久久久久精品久久久久真实原创| 精品久久久久久电影网| 欧美日韩视频精品一区| 伊人亚洲综合成人网| 99视频精品全部免费 在线| videosex国产| 国产一区二区激情短视频 | 日韩精品有码人妻一区| 永久网站在线| 看免费成人av毛片| 久久亚洲国产成人精品v| 亚洲精华国产精华液的使用体验| 十八禁高潮呻吟视频| 成人二区视频| 下体分泌物呈黄色| 亚洲五月色婷婷综合| 香蕉精品网在线| 香蕉精品网在线| 亚洲熟女精品中文字幕| 这个男人来自地球电影免费观看 | 久久99热这里只频精品6学生| av国产精品久久久久影院| 看免费av毛片| 亚洲高清免费不卡视频| 性色avwww在线观看| 国产精品熟女久久久久浪| 日韩 亚洲 欧美在线| 天天躁夜夜躁狠狠躁躁| 国产亚洲一区二区精品| 国产成人欧美| 激情视频va一区二区三区| 亚洲欧美一区二区三区国产| 久久国内精品自在自线图片| 亚洲内射少妇av| 丰满迷人的少妇在线观看| av.在线天堂| 亚洲av电影在线进入| 性高湖久久久久久久久免费观看| 看免费成人av毛片| 最近中文字幕2019免费版| 国产男女超爽视频在线观看| 90打野战视频偷拍视频| 欧美日韩综合久久久久久| 女性生殖器流出的白浆| 亚洲精品美女久久久久99蜜臀 | 午夜精品国产一区二区电影| 国产爽快片一区二区三区| 国产69精品久久久久777片| 国产日韩欧美视频二区| a 毛片基地| 水蜜桃什么品种好| 国产老妇伦熟女老妇高清| 好男人视频免费观看在线| 在线观看免费高清a一片| 婷婷成人精品国产| 晚上一个人看的免费电影| 另类精品久久| 久久久精品区二区三区| 亚洲国产成人一精品久久久| 曰老女人黄片| 韩国精品一区二区三区 | 免费观看性生交大片5| 黄色 视频免费看| 捣出白浆h1v1| 少妇被粗大猛烈的视频| 亚洲,一卡二卡三卡| 欧美日韩亚洲高清精品| 少妇高潮的动态图| 搡女人真爽免费视频火全软件| 精品卡一卡二卡四卡免费| 男人舔女人的私密视频| 久久久国产精品麻豆| 中文字幕人妻熟女乱码| 日日撸夜夜添| 日本av手机在线免费观看| 中文字幕精品免费在线观看视频 | 两性夫妻黄色片 | 国产免费一区二区三区四区乱码| 久久青草综合色| 中文字幕精品免费在线观看视频 | 美女脱内裤让男人舔精品视频| 少妇高潮的动态图| 人成视频在线观看免费观看| 亚洲色图 男人天堂 中文字幕 | 国产成人免费观看mmmm| 最新中文字幕久久久久| 99久久中文字幕三级久久日本| av免费观看日本| 国产黄频视频在线观看| 日产精品乱码卡一卡2卡三| 久久久久久久久久久久大奶| 男女下面插进去视频免费观看 | 最近2019中文字幕mv第一页| 男人舔女人的私密视频| 久久毛片免费看一区二区三区| 插逼视频在线观看| 亚洲天堂av无毛| 日韩大片免费观看网站| 免费大片18禁| 精品一区二区三区视频在线| 亚洲精品中文字幕在线视频| 亚洲三级黄色毛片| 你懂的网址亚洲精品在线观看| 这个男人来自地球电影免费观看 | 国产欧美日韩一区二区三区在线| a级片在线免费高清观看视频| 国产免费一级a男人的天堂| 少妇人妻 视频| 高清av免费在线| 久久热在线av| 一区二区三区乱码不卡18| 一边摸一边做爽爽视频免费| 久久久久久久国产电影| 国产色婷婷99| 91精品三级在线观看| 久久人妻熟女aⅴ| xxx大片免费视频| a级毛片黄视频| 久久久精品区二区三区| 国产又爽黄色视频| 国产免费福利视频在线观看| 免费播放大片免费观看视频在线观看| 97在线视频观看| 久久国产亚洲av麻豆专区| 美女国产高潮福利片在线看| 国产精品久久久av美女十八| 人人妻人人澡人人爽人人夜夜| 一区二区三区精品91| 最新中文字幕久久久久| 国产成人91sexporn| 一区二区av电影网| 亚洲av.av天堂| 最近最新中文字幕免费大全7| 中文精品一卡2卡3卡4更新| 十分钟在线观看高清视频www| 日韩中文字幕视频在线看片| 日韩三级伦理在线观看| av天堂久久9| 美女福利国产在线| 美女内射精品一级片tv| 精品久久久精品久久久| 下体分泌物呈黄色| 精品亚洲乱码少妇综合久久| 美女中出高潮动态图| 在线天堂中文资源库| 中文天堂在线官网| 最黄视频免费看| 少妇的逼好多水| 免费观看性生交大片5| 一二三四在线观看免费中文在 | 欧美3d第一页| videos熟女内射| 午夜视频国产福利| 日本黄色日本黄色录像| 国产精品久久久久久久电影| 日韩av不卡免费在线播放| 日本wwww免费看| 男女国产视频网站| 美女内射精品一级片tv| 黄色视频在线播放观看不卡| 内地一区二区视频在线| 久久精品aⅴ一区二区三区四区 | www日本在线高清视频| 9热在线视频观看99| 中文乱码字字幕精品一区二区三区| 欧美人与性动交α欧美软件 | 亚洲精品美女久久av网站| 黄色配什么色好看| 国产有黄有色有爽视频| av网站免费在线观看视频| 亚洲欧美成人综合另类久久久| 男女国产视频网站| 日韩免费高清中文字幕av| 飞空精品影院首页| 另类亚洲欧美激情| 国产av精品麻豆| 香蕉国产在线看| 国国产精品蜜臀av免费| 久久热在线av| 日本欧美国产在线视频| 欧美 日韩 精品 国产| 国产精品久久久av美女十八| 欧美+日韩+精品| 久久99一区二区三区| 亚洲av成人精品一二三区| 亚洲av电影在线观看一区二区三区| av国产久精品久网站免费入址| a级毛色黄片| 亚洲三级黄色毛片| 国产欧美亚洲国产| 亚洲精品日韩在线中文字幕| av又黄又爽大尺度在线免费看| 视频中文字幕在线观看| 久久99热6这里只有精品| 亚洲国产看品久久| 久久久久久久久久久免费av| 边亲边吃奶的免费视频| 18禁裸乳无遮挡动漫免费视频| 亚洲精品久久午夜乱码| 精品卡一卡二卡四卡免费| 在线免费观看不下载黄p国产| 欧美成人精品欧美一级黄| 国产亚洲av片在线观看秒播厂| 亚洲人成网站在线观看播放| 欧美日韩一区二区视频在线观看视频在线| 亚洲熟女精品中文字幕| 日本爱情动作片www.在线观看| 国产成人欧美| 不卡视频在线观看欧美| 国产国语露脸激情在线看| 亚洲色图综合在线观看| 欧美精品亚洲一区二区| 国产日韩欧美在线精品| 亚洲av国产av综合av卡| 久久ye,这里只有精品| 咕卡用的链子| 一级爰片在线观看| 天天躁夜夜躁狠狠久久av| 一级毛片电影观看| 国产精品久久久久久精品古装| 亚洲中文av在线| 自拍欧美九色日韩亚洲蝌蚪91| 热re99久久精品国产66热6| 捣出白浆h1v1| 国产成人91sexporn| 亚洲国产色片| 伦精品一区二区三区| 最新中文字幕久久久久| 9191精品国产免费久久| 黑人巨大精品欧美一区二区蜜桃 | 一边亲一边摸免费视频| 51国产日韩欧美| 人人妻人人爽人人添夜夜欢视频| 精品久久国产蜜桃| 男女无遮挡免费网站观看| www.色视频.com| 中文字幕制服av| 亚洲精品,欧美精品| 亚洲欧美中文字幕日韩二区| 边亲边吃奶的免费视频| 国产精品久久久久久精品电影小说| 大陆偷拍与自拍| 成人漫画全彩无遮挡| 亚洲精品久久成人aⅴ小说| 一本色道久久久久久精品综合| 老熟女久久久| 又大又黄又爽视频免费| 熟妇人妻不卡中文字幕| 少妇人妻 视频| 欧美成人精品欧美一级黄| 91国产中文字幕| 侵犯人妻中文字幕一二三四区| 永久免费av网站大全| 99久久精品国产国产毛片| 久久国产精品大桥未久av| 黄色毛片三级朝国网站| 成人免费观看视频高清| 热re99久久国产66热| 国产欧美亚洲国产| 精品少妇久久久久久888优播| a级片在线免费高清观看视频| 欧美另类一区| 亚洲精华国产精华液的使用体验| 制服人妻中文乱码| kizo精华| 黄片播放在线免费| 日韩成人伦理影院| 观看美女的网站| 久久亚洲国产成人精品v| 亚洲精品国产色婷婷电影| 免费人成在线观看视频色| 大码成人一级视频| 精品一区二区三卡| 有码 亚洲区| 精品一区二区三卡| 男女边吃奶边做爰视频| 午夜视频国产福利| 国产高清国产精品国产三级| 国产一区有黄有色的免费视频| 黄片无遮挡物在线观看| 亚洲欧美日韩卡通动漫| 晚上一个人看的免费电影| 久久久久国产网址| 亚洲国产av新网站| 一二三四在线观看免费中文在 | 亚洲精品久久久久久婷婷小说| 伦理电影免费视频| 欧美丝袜亚洲另类| 不卡视频在线观看欧美| 日产精品乱码卡一卡2卡三| 亚洲av电影在线观看一区二区三区| 久久精品久久久久久久性| 九草在线视频观看| 亚洲欧美成人综合另类久久久| 99九九在线精品视频| 乱人伦中国视频| 在线看a的网站| 亚洲欧洲日产国产| 麻豆精品久久久久久蜜桃| 在线精品无人区一区二区三| 美女主播在线视频| 国产亚洲精品久久久com| 人妻 亚洲 视频| 亚洲av欧美aⅴ国产| 91久久精品国产一区二区三区| 欧美+日韩+精品| 亚洲一级一片aⅴ在线观看| 人妻少妇偷人精品九色| 国产淫语在线视频| 母亲3免费完整高清在线观看 | 99re6热这里在线精品视频| 国产熟女午夜一区二区三区| 黄色毛片三级朝国网站| 国产日韩欧美亚洲二区| 日本-黄色视频高清免费观看| 亚洲欧洲精品一区二区精品久久久 | 大片电影免费在线观看免费| 国产精品久久久久久精品电影小说| 久久精品熟女亚洲av麻豆精品| 日韩av不卡免费在线播放| 2018国产大陆天天弄谢| a 毛片基地| 久久久久久久久久久免费av| 大片电影免费在线观看免费| 老熟女久久久| 亚洲精品aⅴ在线观看| 国国产精品蜜臀av免费| 亚洲五月色婷婷综合| 高清视频免费观看一区二区| 26uuu在线亚洲综合色| 久久久久人妻精品一区果冻| 国产精品久久久久久精品古装| 在线天堂中文资源库| 人人妻人人澡人人看| 啦啦啦视频在线资源免费观看| 久久99热6这里只有精品| 久热这里只有精品99| 最新的欧美精品一区二区| 欧美性感艳星| 丝瓜视频免费看黄片| 亚洲在久久综合| 哪个播放器可以免费观看大片| 国产精品欧美亚洲77777| 久久久精品区二区三区| 免费看光身美女| 日韩制服骚丝袜av| 最近的中文字幕免费完整| 国产成人精品福利久久| 免费av中文字幕在线| 香蕉精品网在线| 人成视频在线观看免费观看| 黄色毛片三级朝国网站| 亚洲精品久久成人aⅴ小说| 国产一区二区在线观看av| 欧美3d第一页| 精品人妻偷拍中文字幕| 成年人午夜在线观看视频| 国产成人精品在线电影| 两个人看的免费小视频| 亚洲精品久久成人aⅴ小说| 最近最新中文字幕免费大全7| 中文字幕亚洲精品专区| 亚洲av中文av极速乱| 99久国产av精品国产电影| 成人亚洲欧美一区二区av| 精品少妇久久久久久888优播| 成年动漫av网址| 一本—道久久a久久精品蜜桃钙片| 男人舔女人的私密视频| 亚洲色图 男人天堂 中文字幕 | 久久久久久久大尺度免费视频| freevideosex欧美| 天堂俺去俺来也www色官网| 五月伊人婷婷丁香| 久久久国产欧美日韩av| 一本色道久久久久久精品综合| 精品视频人人做人人爽| 午夜激情久久久久久久| 午夜福利视频在线观看免费| 免费播放大片免费观看视频在线观看| 下体分泌物呈黄色| 亚洲av电影在线进入| 欧美性感艳星| 精品一区二区三卡| 一本色道久久久久久精品综合| 国产爽快片一区二区三区| 午夜精品国产一区二区电影| 视频在线观看一区二区三区| 五月天丁香电影| 日本免费在线观看一区| 国产成人a∨麻豆精品| 亚洲国产日韩一区二区| 在线观看美女被高潮喷水网站| 国产成人午夜福利电影在线观看| 国产激情久久老熟女| 国产一区有黄有色的免费视频| 亚洲av福利一区| 亚洲欧美成人精品一区二区| 91精品国产国语对白视频| 狠狠婷婷综合久久久久久88av| 黄色视频在线播放观看不卡| 免费观看a级毛片全部| 男女午夜视频在线观看 | 国产精品国产av在线观看| 免费大片18禁| 国产av码专区亚洲av| av女优亚洲男人天堂| av播播在线观看一区| 99久久人妻综合| 一级片'在线观看视频| 久久av网站| 美女内射精品一级片tv| 国产永久视频网站| 中文字幕亚洲精品专区| 亚洲av福利一区| 久久av网站| 国产免费一区二区三区四区乱码| 天天躁夜夜躁狠狠久久av| 久热久热在线精品观看| 国产精品国产三级国产av玫瑰| 国产xxxxx性猛交| 搡老乐熟女国产| 成人亚洲精品一区在线观看| 亚洲欧美中文字幕日韩二区| 观看av在线不卡| 国产在线一区二区三区精| 国产激情久久老熟女| 免费看光身美女| 日韩欧美精品免费久久| freevideosex欧美| 亚洲国产精品国产精品| 精品一区在线观看国产| 亚洲成色77777| 国产 精品1| 99香蕉大伊视频| 天天影视国产精品| 两个人看的免费小视频| 精品国产露脸久久av麻豆| 色婷婷久久久亚洲欧美| 成人手机av| 伊人久久国产一区二区| 国产亚洲最大av| 亚洲成色77777| 国产精品99久久99久久久不卡 | 免费观看av网站的网址| 99热这里只有是精品在线观看| 亚洲欧洲国产日韩| 成人午夜精彩视频在线观看| 少妇猛男粗大的猛烈进出视频| 午夜福利网站1000一区二区三区| 欧美 日韩 精品 国产| 亚洲一级一片aⅴ在线观看| 欧美老熟妇乱子伦牲交| 精品国产国语对白av| 天天躁夜夜躁狠狠久久av| 精品99又大又爽又粗少妇毛片| 国产成人免费观看mmmm| 成年av动漫网址| 成人国产av品久久久| 建设人人有责人人尽责人人享有的| 日本wwww免费看| 久久精品国产亚洲av天美| 日韩av免费高清视频| 黄色配什么色好看| 国产精品一区www在线观看| 精品国产乱码久久久久久小说| 男人添女人高潮全过程视频| 国精品久久久久久国模美| 亚洲欧美精品自产自拍| 香蕉精品网在线| 久久久国产一区二区| 97超碰精品成人国产| 精品少妇久久久久久888优播| videosex国产| 国语对白做爰xxxⅹ性视频网站| 国产熟女欧美一区二区| 美国免费a级毛片| 色94色欧美一区二区| 99久久综合免费| 久久女婷五月综合色啪小说| 国产免费一区二区三区四区乱码| 少妇高潮的动态图| 久久99热这里只频精品6学生| 99re6热这里在线精品视频| 国精品久久久久久国模美| 国产欧美亚洲国产| 99久久综合免费| 一级毛片黄色毛片免费观看视频| 婷婷色综合www| 精品福利永久在线观看| 啦啦啦啦在线视频资源| 99热国产这里只有精品6| 飞空精品影院首页| 日本午夜av视频| 999精品在线视频| 观看av在线不卡| 亚洲精品中文字幕在线视频| 亚洲精品久久久久久婷婷小说| 久久精品夜色国产| 伊人亚洲综合成人网| 色婷婷久久久亚洲欧美| 午夜福利网站1000一区二区三区| 黄色视频在线播放观看不卡| 这个男人来自地球电影免费观看 | 日韩在线高清观看一区二区三区| 欧美少妇被猛烈插入视频| 青春草视频在线免费观看| 啦啦啦在线观看免费高清www| 国产精品麻豆人妻色哟哟久久| 色婷婷av一区二区三区视频| 亚洲美女搞黄在线观看| 国产麻豆69| 乱码一卡2卡4卡精品| 少妇的丰满在线观看| 国产高清三级在线| 18禁裸乳无遮挡动漫免费视频| 一级爰片在线观看| 国产精品人妻久久久影院| 又黄又爽又刺激的免费视频.| av片东京热男人的天堂| 亚洲精品一区蜜桃| 国产成人91sexporn| 夜夜爽夜夜爽视频| 一级a做视频免费观看| 亚洲,欧美,日韩| 精品久久久精品久久久| 精品熟女少妇av免费看| 大香蕉97超碰在线| 男女下面插进去视频免费观看 | videosex国产| 精品卡一卡二卡四卡免费| 久久国产精品大桥未久av| videossex国产| 国产av精品麻豆| 日本欧美国产在线视频| 精品人妻偷拍中文字幕| 天堂俺去俺来也www色官网| 亚洲天堂av无毛| 九九爱精品视频在线观看| av片东京热男人的天堂| 国产日韩欧美在线精品| xxx大片免费视频| 久久这里有精品视频免费| 欧美精品av麻豆av| 久久亚洲国产成人精品v| 亚洲国产精品一区三区| 国产在线视频一区二区| 免费观看无遮挡的男女| 五月开心婷婷网| 成人无遮挡网站|