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

    微藻培養(yǎng)光生物反應(yīng)器內(nèi)傳遞現(xiàn)象的研究進(jìn)展

    2012-10-19 03:35:48陳智杰姜澤毅張欣欣張欣茹
    化工進(jìn)展 2012年7期
    關(guān)鍵詞:生物

    陳智杰,姜澤毅,2,張欣欣,3,張欣茹

    (1北京科技大學(xué)機(jī)械工程學(xué)院,北京 100083;2北京科技大學(xué),北京市高校節(jié)能與環(huán)保工程研究中心,北京100083; 3北京科技大學(xué),冶金工業(yè)節(jié)能減排北京市重點(diǎn)實(shí)驗(yàn)室,北京 100083)

    進(jìn)展與述評(píng)

    微藻培養(yǎng)光生物反應(yīng)器內(nèi)傳遞現(xiàn)象的研究進(jìn)展

    陳智杰1,姜澤毅1,2,張欣欣1,3,張欣茹1

    (1北京科技大學(xué)機(jī)械工程學(xué)院,北京 100083;2北京科技大學(xué),北京市高校節(jié)能與環(huán)保工程研究中心,北京100083;3北京科技大學(xué),冶金工業(yè)節(jié)能減排北京市重點(diǎn)實(shí)驗(yàn)室,北京 100083)

    微藻規(guī)?;囵B(yǎng)過程中光生物反應(yīng)器內(nèi)傳遞現(xiàn)象是影響微藻的生長(zhǎng)及產(chǎn)量的重要因素。本文重點(diǎn)綜述了光生物反應(yīng)器內(nèi)傳遞現(xiàn)象(光傳遞、傳熱、傳質(zhì)和傳動(dòng)量傳遞)及其數(shù)學(xué)模型研究進(jìn)展,分析了光生物反應(yīng)器結(jié)構(gòu)和尺寸對(duì)光傳遞和傳質(zhì)的影響,總結(jié)影響各傳遞現(xiàn)象的重要參數(shù),如光吸收系數(shù)、體積傳質(zhì)系數(shù)等,為高效光生物反應(yīng)器的設(shè)計(jì)、優(yōu)化及放大提供了參考依據(jù)。

    微藻;光生物反應(yīng)器;傳熱;傳質(zhì);光衰減

    微藻(microalgae)是一種能有效利用光能、二氧化碳和水生長(zhǎng)的低等植物,其種類繁多、分布廣泛;具有生物量大、生長(zhǎng)周期短、易培養(yǎng)及脂質(zhì)含量高等特點(diǎn),被國際上認(rèn)為是一種最有潛力替代石油的生物資源[1],其規(guī)?;囵B(yǎng)是解決當(dāng)前能源短缺和環(huán)境污染的有效手段。

    微藻規(guī)?;庾责B(yǎng)培養(yǎng)主要在光生物反應(yīng)器內(nèi)進(jìn)行,以提供微藻生長(zhǎng)必要的營養(yǎng)物、光照、CO2、溫度、鹽度和通氣條件等。微藻光生物反應(yīng)器包括開放式和封閉式[2-4],高效光生物反應(yīng)器的設(shè)計(jì)、優(yōu)化和放大是能源微藻規(guī)?;钠款i問題之一[1]。高效光生物反應(yīng)器研究主要涉及4個(gè)方面[5]:一是增大反應(yīng)器比表面積;二是增強(qiáng)氣液傳質(zhì)效率;三是提供高效光源;四是提高光傳遞效率。因此,深入研究光生物反應(yīng)器內(nèi)流動(dòng)、傳熱、傳質(zhì)和光傳遞現(xiàn)象對(duì)于實(shí)現(xiàn)能源微藻規(guī)?;哂兄匾F(xiàn)實(shí)意義。有鑒于此,本文綜述了光生物反應(yīng)器內(nèi)的各傳遞現(xiàn)象及其數(shù)理模型,闡述了光生物反應(yīng)器結(jié)構(gòu)和尺寸變化對(duì)各傳遞現(xiàn)象的影響,并分析了主要參數(shù)對(duì)各相關(guān)傳遞現(xiàn)象的影響。

    1 光生物反應(yīng)器內(nèi)光傳遞現(xiàn)象

    1.1 現(xiàn)象描述

    植物生長(zhǎng)離不開光合作用,光的可獲得性和光照強(qiáng)度是控制光合作用影響微藻細(xì)胞生長(zhǎng)的主要因素[6-7]。然而,光生物反應(yīng)器結(jié)構(gòu)多樣,且微藻細(xì)胞自身存在相互遮蔽和光吸收效應(yīng),使入射光穿過藻液時(shí)不斷衰減。因此,光在培養(yǎng)體系內(nèi)的分布存在時(shí)空的非線性變化規(guī)律,研究發(fā)現(xiàn)隨著培養(yǎng)細(xì)胞密度和透光距離的增加光衰減現(xiàn)象呈上升趨勢(shì)[8-10]。

    對(duì)于微藻戶外培養(yǎng),不僅存在著光衰減現(xiàn)象,而且光照強(qiáng)度和光照時(shí)間均隨時(shí)間呈現(xiàn)周期性變化。在光生物反應(yīng)器研究中,有學(xué)者將培養(yǎng)液的流動(dòng)和光衰減結(jié)合,采用光暗循環(huán)周期[7,11]來體現(xiàn)光生物反應(yīng)器內(nèi)光分布,這一參數(shù)直接影響微藻的生產(chǎn)率和CO2固定率。此外,過量的光照強(qiáng)度還會(huì)引起光抑制作用[12-13]。

    因此,建立光生物反應(yīng)器光傳遞數(shù)學(xué)模型有效分析光分布是設(shè)計(jì)高效微藻光生物反應(yīng)器的關(guān)鍵之一。

    1.2 數(shù)學(xué)模型

    目前,學(xué)者們通常假設(shè)[8,10]光在培養(yǎng)液中傳播方向不變,忽略固體顆粒的光散射,輻射為單色光,采用經(jīng)典光學(xué)理論中Lambert-Beer定律[14]來描述光生物反應(yīng)器內(nèi)的細(xì)胞濃度和光徑對(duì)光衰減的一般規(guī)律,如式(1)。

    式中,I為入射光被藻類吸收后剩余的光強(qiáng)度,W/m2;I0為入射光強(qiáng)度,W/m2;α為光吸收系數(shù),m-1;L為光徑,m;ka為比消光系數(shù),m2/g;X為藻細(xì)胞濃度,g/L。

    近年來,國內(nèi)外有不少學(xué)者利用Lambert-Beer定律的指數(shù)衰減關(guān)系對(duì)多種藻細(xì)胞培養(yǎng)液中的光衰減進(jìn)行研究,得到光衰減經(jīng)驗(yàn)公式,如式(2)~式(5)。

    一般情況下,采用Lambert-Beer定律指數(shù)關(guān)系已足夠準(zhǔn)確描述培養(yǎng)液中光分布。然而,由于光散射和可選擇性吸收,隨著藻細(xì)胞濃度的增加,藻細(xì)胞間距逐漸減小,某些情況下細(xì)胞濃度與光衰減的變化關(guān)系會(huì)偏離Lambert-Beer定律[19,20],呈非線性相關(guān)。

    式中,I0為入射光光照強(qiáng)度,W/m2;I為光程L位置的光照強(qiáng)度,W/m2;OD為光學(xué)密度;L為光程,m。

    Fernandez等[8]在研究三角褐指藻時(shí)發(fā)現(xiàn),當(dāng)藻細(xì)胞濃度超過1.3 g/L時(shí),藻液的光衰減特性不符合Lambert-Beer定律,提出描述光衰減現(xiàn)象較為簡(jiǎn)單的雙曲模型,見式(8)。

    式中,At為光衰減系數(shù);Atmax為最大光衰減度;Kat為光衰減常數(shù),g/m;K'a為雙曲吸收系數(shù),m-1;Laverage為平均光徑,m;c為細(xì)胞濃度,g/L。

    華東理工大學(xué)朱篤等[21]利用雙曲模型研究聚球藻7942培養(yǎng)液中的光衰減,得到較Lambert-Beer定律更好的光分布。暨南大學(xué)徐明芳等[22]在研究螺旋藻時(shí)發(fā)現(xiàn),采用雙曲模型能在較大藻生物量濃度范圍內(nèi)很好的描述光衰減現(xiàn)象。

    此外,江西師范大學(xué)張志斌[23]基于藻細(xì)胞中葉綠素對(duì)光的吸收存在的“包裹效應(yīng)”理論,假設(shè)比消光系數(shù)隨藻細(xì)胞濃度增大趨于一個(gè)定值,從而得到光衰減指數(shù)模型,見式(9)。實(shí)驗(yàn)證明該模型在較寬藻細(xì)胞濃度下對(duì)不同藻培養(yǎng)液光衰減都有較好的描述。

    式中,A(X)為光衰減系數(shù),m-1;為細(xì)胞濃度X趨于無窮大時(shí)的比消光系數(shù),m2/g;X為藻細(xì)胞濃度,g/L;C為常數(shù),C=;ε為經(jīng)驗(yàn)常數(shù)。

    1.3 光吸收系數(shù)

    光吸收系數(shù)是 Lambert-Beer定律中的一個(gè)參數(shù),它與光強(qiáng)度無關(guān),主要用于表征光衰減程度。華東理工大學(xué)劉晶璘等[24]分析螺旋藻時(shí)發(fā)現(xiàn),培養(yǎng)基的光吸收系數(shù)與藻液自身的光吸收系數(shù)之比小于1%,藻液的光吸收系數(shù)主要由藻體決定。不同的藻種由于藻體細(xì)胞的結(jié)構(gòu)、化學(xué)組成及光合色素的種類與含量不同,光譜吸收特性不同。

    在低培養(yǎng)密度時(shí),目前普遍采用以特征波長(zhǎng)表示的藻細(xì)胞濃度擬合得到光吸收系數(shù),見式(10)。

    式中,α為總光吸收系數(shù),m-1;ODi為以消光度表示的藻細(xì)胞密度;A,B為常數(shù),是光衰減系數(shù);Ex為比消光系數(shù),m2/g;cx為藻細(xì)胞濃度,g/L。

    高細(xì)胞濃度時(shí),消光系數(shù)可采用雙曲模型來定義[21],見式(11)。

    式中,A(X)為光衰減系數(shù),m―1;Amax=670 m-1,b=2.25 g/L;X為藻細(xì)胞濃度,g/L。

    清華大學(xué)吳良柏[25]在Bricaud的基礎(chǔ)上,考慮不同波長(zhǎng)光波的能量對(duì)平均吸收系數(shù)的影響,將太陽光照近似為溫度5800 K的黑體輻射,得到小球藻液接受太陽光照的平均光吸收系數(shù),見式(12)。

    式中,α(λ)為隨波長(zhǎng)變化的光吸收系數(shù),m-1;λ為波長(zhǎng),m;c1為第一輻射常量,其值為3.7419×10-16,W/m2;c2為第二輻射常量,其值為 1.4388×10-2,W/m2;Ts為 5000K;σ為玻爾茲曼常數(shù),5.67× 10-8W/(m2/K4)。

    華東理工大學(xué)劉晶璘等[24]研究了光吸收系數(shù)的波長(zhǎng)分布特征,得到隨著波長(zhǎng)的增加光吸收系數(shù)基本呈現(xiàn)平穩(wěn)下降趨勢(shì)。同時(shí),基于藻體大小對(duì)光吸收系數(shù)的影響規(guī)律,提出了藻體比消光率的概念,得到藻體顆粒越大的藻種其藻液的比消光系數(shù)越小,在相同幾何尺寸的光生物反應(yīng)器中藻液的平均光強(qiáng)就越高。Cornet[24]提出,當(dāng)藻細(xì)胞濃度較低時(shí),比消光系數(shù)(ka或k)是一個(gè)常數(shù),隨著藻液濃度的增加,比消光系數(shù)不再是一個(gè)常數(shù),會(huì)隨著藻液濃度的增加發(fā)生變化。張志斌[23]實(shí)驗(yàn)研究得到比消光系數(shù)隨著光程的增大而下降,隨細(xì)胞濃度的增大而減少。

    2 光生物反應(yīng)器內(nèi)傳熱現(xiàn)象

    2.1 現(xiàn)象描述

    培養(yǎng)液溫度是微藻細(xì)胞生長(zhǎng)的主要限制因子,同時(shí)影響培養(yǎng)液中的組分?jǐn)U散。光生物反應(yīng)器內(nèi)的溫度分布[26-27]涉及氣液流動(dòng)、輻射換熱、蒸發(fā)散熱、生物反應(yīng)熱等復(fù)雜的物理化學(xué)過程。一般來說,每個(gè)藻種都有其最佳的生長(zhǎng)溫度范圍。Richmod[28]研究發(fā)現(xiàn),當(dāng)溫度低于最佳生長(zhǎng)溫度15 ℃時(shí)微藻仍能生存,但當(dāng)溫度超過其最佳生長(zhǎng)溫度2~4 ℃時(shí)就會(huì)造成藻細(xì)胞的大量死亡。

    對(duì)于封閉式光生物反應(yīng)器來說,培養(yǎng)液溫度很容易達(dá)到 55 ℃。當(dāng)反應(yīng)器尺寸增大反應(yīng)器體積增大,單位體積的表面積迅速減小,而冷卻效率卻正比于反應(yīng)器表面積,因此傳熱問題成為光生物反應(yīng)器規(guī)?;南拗埔蛩豙3]。

    2.2 數(shù)學(xué)模型

    Gutierrez 等[29]忽略太陽與培養(yǎng)液、培養(yǎng)液與外界環(huán)境間的輻射換熱,以及培養(yǎng)液內(nèi)上下液面間的對(duì)流換熱,考慮培養(yǎng)液與壁面的傳熱,建立數(shù)學(xué)模型,如式(13)、式(14)。

    式中,Qs、Qst、Qc、Qe、Qr、Qk分別指培養(yǎng)液獲得的太陽熱、壁面獲得的太陽熱、對(duì)流熱損失量、蒸發(fā)熱損失量、培養(yǎng)液的輻射熱流量、導(dǎo)熱損失熱流量,W;Mw、Mt分別指培養(yǎng)液和壁面的質(zhì)量,kg;Cpw、Cpt分別指培養(yǎng)液和壁面的比熱容,kJ/(kg?K);Tw、Tt分別指培養(yǎng)液和壁面的溫度,K;Uta、Uwt分別指培養(yǎng)液與環(huán)境和壁面之間的對(duì)流換熱系數(shù),

    W /(m2?K)。

    Quentin等[27]綜合考慮了太陽輻射、大氣輻射、培養(yǎng)液蒸發(fā)等因素,建立數(shù)學(xué)模型如式(15),模擬結(jié)果與實(shí)驗(yàn)結(jié)果都表明,對(duì)于開放池式光生物反應(yīng)器,反應(yīng)器內(nèi)培養(yǎng)液的溫度與環(huán)境溫度有密切關(guān)系,培養(yǎng)液溫度峰值在每天的13時(shí)到16時(shí)之間。

    式中,Qra,r、Qra,d、Qra,D、Qre,s、Qra,a、Qre,a、

    Qra,g、Qc、Qev、Qb、Qcond分別指反應(yīng)器的輻射熱、太陽擴(kuò)散輻射熱、直接太陽輻射熱、地面的太陽反輻熱、空氣輻射熱、地面輻射熱、對(duì)流換熱、蒸發(fā)熱、氣泡的換熱量和導(dǎo)熱,W;ρw為培養(yǎng)液的密度,kg/m;Vr為反應(yīng)器體積,m3;Cpw為培養(yǎng)液比熱容,kJ/(kg?℃);Tr為培養(yǎng)液溫度,℃。

    目前,關(guān)于光生物反應(yīng)器內(nèi)溫度分布研究主要針對(duì)開放式且采用零維模型[27,29],仍無法準(zhǔn)確描述微藻培養(yǎng)過程中培養(yǎng)液上下溫度差異。與此同時(shí),為了簡(jiǎn)化傳熱問題 Gotez[30]在無藻情況下分析藻液的熱平衡,建立了封閉式光生物反應(yīng)器的零維數(shù)學(xué)模型。由于封閉式光生物反應(yīng)器內(nèi)溫度分布非均勻性[6],應(yīng)綜合考慮培養(yǎng)液光學(xué)特性、流動(dòng)狀態(tài)及藻細(xì)胞生物反應(yīng)熱建立沿培養(yǎng)液高度方向的溫度模型。

    2.3 相關(guān)參數(shù)

    2.3.1 熱導(dǎo)率

    熱導(dǎo)率是藻液流體的一個(gè)重要熱物性參數(shù),是建立光生物反應(yīng)器內(nèi)傳熱數(shù)學(xué)模型的必要參數(shù),目前罕有對(duì)藻液熱導(dǎo)率測(cè)量報(bào)道。清華大學(xué)吳良柏[25]采用 3ω法在 20~38 ℃之間測(cè)量小球藻濃度為C(C=2.7×107個(gè)/m L)和C/2流體的熱導(dǎo)率,得到熱導(dǎo)率在0.58 W/(m?K)到0.62 W/(m?K)之間。

    2.3.2 比熱容

    比熱容是材料的物性參數(shù)之一。上海交通大學(xué)王爽等[31]利用NETZSCH DSC404型差示量熱掃描儀測(cè)量了3種海藻(紅藻、綠藻和褐藻)粉末在40~550 ℃溫度范圍內(nèi)的比熱容。目前,仍末見藻液流體比熱容測(cè)量報(bào)道,實(shí)際計(jì)算中較多采用近似于水的比熱容。

    3 光生物反應(yīng)器內(nèi)氣液傳質(zhì)現(xiàn)象

    3.1 現(xiàn)象描述

    微藻培養(yǎng)系統(tǒng)是典型的多相反應(yīng)體系,存在明顯的氣相、液相和固相。微藻培養(yǎng)光生物反應(yīng)器內(nèi)空氣或CO2以氣泡形式存在于液相中,以提供光合作用所需的無機(jī)碳源。光合作用的產(chǎn)物O2必須及時(shí)排除,以避免由于溶氧濃度[3]過高而導(dǎo)致藻細(xì)胞死亡。由于光生物反應(yīng)器內(nèi)多相流動(dòng)與反應(yīng)器結(jié)構(gòu)和尺寸有密切聯(lián)系,隨著光生物反應(yīng)器結(jié)構(gòu)和尺寸發(fā)生變化,反應(yīng)器內(nèi)無機(jī)碳和溶氧濃度分布將發(fā)生明顯的改變。因此,深入分析供氣條件、光生物反應(yīng)器結(jié)構(gòu)對(duì)CO2和O2體積傳質(zhì)速率的影響規(guī)律,建立無機(jī)碳和溶氧濃度數(shù)學(xué)模型對(duì)優(yōu)化光生物反應(yīng)器結(jié)構(gòu)和提高氣液傳質(zhì)效率具有重要的實(shí)際意義。

    光生物反應(yīng)器內(nèi)氣液傳質(zhì)包括氣體主流中CO2擴(kuò)散到藻細(xì)胞光合反應(yīng)部位,以及光合作用產(chǎn)物O2從液相主體中排除。二氧化碳從氣泡向藻細(xì)胞內(nèi)光合作用部位的主要擴(kuò)散過程如圖1所示[32]。

    圖1 微藻培養(yǎng)過程的二氧化碳傳遞過程示意圖

    3.2 數(shù)理模型

    Rubio等[33]對(duì)氣升式管狀生物反應(yīng)器進(jìn)行研究,認(rèn)為微藻培養(yǎng)過程是一個(gè)偽均勻的氣液兩相系統(tǒng),并假設(shè):①在極小的時(shí)間間隔內(nèi),光合作用率(溶氧濃度和無機(jī)碳濃度)只隨時(shí)間變化;②根據(jù)流體靜力學(xué),假設(shè)液相的體積流量為常數(shù);③認(rèn)為與培養(yǎng)液同溫度下,CO2和O2呈飽和狀態(tài);④認(rèn)為大部分的氣升區(qū)是柱塞流,除了脫氣裝置。根據(jù)組分的質(zhì)量守恒得到無機(jī)碳和溶氧的控制方程,如式(16)。

    該模型能夠以小時(shí)為單位,同時(shí)預(yù)測(cè)上升管和下降管內(nèi)的溶氧水平、CO2的供給量以及pH值變化情況。

    浙江大學(xué)程桂林等[34]對(duì)10 L氣升式柱型光生物反應(yīng)器進(jìn)行研究,假設(shè):①液相連續(xù)流動(dòng),并忽略氣液間相互作用;②氣泡在液相中分布均勻,且軸向?qū)ΨQ;③實(shí)驗(yàn)過程認(rèn)為溫度、蒸汽壓恒定。根據(jù)質(zhì)量守恒定律運(yùn)用集中參數(shù)法建立了無機(jī)碳和溶氧的控制方程,如式(17)、式(18)。

    式中,[O2]、[CT]為溶氧、溶解的總碳濃度,mol/m3;[O2]*、[CO2]*為液相中氧氣和二氧化碳的飽和濃度,mol/m3;FO2、FCO2為氣相中氧氣和二氧化碳的摩爾流率,mol/s;QL為液體的體積流率,m3/s;KLaLO2、KLaLCO2分別為氧氣和二氧化碳的體積傳質(zhì)系數(shù),s-1;RO2、RCO2分別為氧氣的生成率和二氧化碳的消耗率,mol?m3/s;DO2、DCO2分別為氧氣和二氧化碳的擴(kuò)散系數(shù),m2?s;s為截面面積,m2;UL為液體的流速,m/s;ε為氣含率。

    研究結(jié)果表明上述模型能針對(duì)各自光生物反應(yīng)器結(jié)構(gòu)較好的模擬溶氧濃度、pH值隨時(shí)間的變化情況,不僅可用于微藻光生物反應(yīng)器培養(yǎng)去除密閉空氣中微量CO2的模擬和預(yù)測(cè),且對(duì)光生物反應(yīng)器的優(yōu)化設(shè)計(jì)和微藻的高密度培養(yǎng)有指導(dǎo)意義。但是,上述模型均為零維模型,尚末能精確的描述無機(jī)碳和溶氧濃度的分布,應(yīng)結(jié)合CFD技術(shù)采用多相流動(dòng)數(shù)學(xué)模型,綜合分析不同光生物反應(yīng)器結(jié)構(gòu)對(duì)無機(jī)碳和溶氧濃度分布的影響。

    3.3 體積傳質(zhì)系數(shù)

    化學(xué)工程中,常采用體積傳質(zhì)系數(shù)來代表反應(yīng)器內(nèi)的總括傳質(zhì)系數(shù)。體積傳質(zhì)系數(shù)是用來評(píng)價(jià)光生物反應(yīng)器性能的常用參數(shù)[2],它與許多因素相關(guān),如表觀氣速[35]、黏度[35]、曝氣量及曝氣方式[36]、擴(kuò)散系數(shù)[37]、氣含率[33-34]等。

    體積傳質(zhì)系數(shù)(KLa)是兩個(gè)參數(shù)的乘積,實(shí)際應(yīng)用中由于測(cè)量傳質(zhì)系數(shù)(KL)和有效接觸面積(a)比較困難,因此常將它們并為一項(xiàng)稱為體積傳質(zhì)系數(shù),單位為 h-1。根據(jù)氣液傳質(zhì)理論,描述體積傳質(zhì)系數(shù)的模型有4種[38]:① 結(jié)合尺寸分析和實(shí)驗(yàn)數(shù)據(jù)的與現(xiàn)象相關(guān)的模擬;② 采用雙膜理論的空間模型;③ 采用Higbie滲透理論的時(shí)間模型;④ 結(jié)合雙膜理論和滲透理論的模型。關(guān)于體積傳質(zhì)系數(shù)的測(cè)量方法[32]有3種:亞硫酸鹽法、動(dòng)態(tài)法和定態(tài)法,且以測(cè)量O2的體積傳質(zhì)系數(shù)為主。對(duì)于光生物反應(yīng)器內(nèi)CO2體積傳質(zhì)系數(shù)的確定常采用經(jīng)驗(yàn)公式法[33-34]。

    4 光生物反應(yīng)器內(nèi)多相流動(dòng)

    4.1 現(xiàn)象描述

    微藻細(xì)胞光自養(yǎng)培養(yǎng)時(shí)需保證藻細(xì)胞在培養(yǎng)液中均勻懸浮,通常采用氣升、鼓泡或轉(zhuǎn)槳攪拌的方式實(shí)現(xiàn)。培養(yǎng)液的流動(dòng)狀態(tài)會(huì)影響培養(yǎng)液混合時(shí)間、剪切(混合)速率和微藻細(xì)胞的光暗循環(huán)周期,較佳的混合時(shí)間可以保證高的細(xì)胞濃度、使藻細(xì)胞懸浮、消除熱分層、幫助營養(yǎng)物均勻分布、提高氣液傳質(zhì)和降低相互遮蔽程度;合理的培養(yǎng)液剪切(混合)速率可以提高氣液傳質(zhì)效率促進(jìn)細(xì)胞生長(zhǎng),但過高的剪切(混合)速率會(huì)損傷藻細(xì)胞從而導(dǎo)致藻細(xì)胞死亡;光暗循環(huán)周期的長(zhǎng)短是影響微藻生產(chǎn)量的重要因素。

    此外,根據(jù)生物反應(yīng)器的放大研究可知[39-40],隨著應(yīng)器結(jié)構(gòu)和尺寸的變化,反應(yīng)器內(nèi)部分的流動(dòng)與傳質(zhì)過程會(huì)發(fā)生明顯的變化。因此,深入分析光生物反應(yīng)器結(jié)構(gòu)和尺寸對(duì)多相流動(dòng)的影響,不僅有助于提高微藻生長(zhǎng)速度[3],而且對(duì)于光生物反應(yīng)器的優(yōu)化和放大具有十分重要的意義。

    4.2 CFD模型

    關(guān)于生物反應(yīng)器內(nèi)流動(dòng)的研究,前人多采用理論結(jié)合實(shí)驗(yàn)研究的方法。隨著計(jì)算機(jī)技術(shù)的發(fā)展,計(jì)算流體力學(xué)(CFD)軟件廣泛應(yīng)用于流場(chǎng)模擬。近年來,有學(xué)者將其用于模擬不同結(jié)構(gòu)光生物反應(yīng)器內(nèi)的流動(dòng)過程,表1列舉了近年CFD技術(shù)在光生物反應(yīng)器中的應(yīng)用。

    表1 CFD技術(shù)在光生物反應(yīng)器內(nèi)流動(dòng)的應(yīng)用

    4.3 相關(guān)參數(shù)

    4.3.1 表觀氣/液速

    對(duì)于氣動(dòng)力推動(dòng)的氣升式光生物反應(yīng)器來說,液體的循環(huán)時(shí)間強(qiáng)烈依賴于表觀氣速,隨著表觀氣速的提高指數(shù)下降[47]。提高表觀氣速不僅有利于氣液的混合、提高液體循環(huán)速度和增加氣液傳質(zhì),而且還會(huì)阻止反應(yīng)器溶氧的積聚。表2給出了不同光生物反應(yīng)器內(nèi)表觀氣速和體積傳質(zhì)系數(shù)的比較。從表中可以看出,在相同表觀氣速下,板式光生物反應(yīng)器具有較高的體積傳質(zhì)系數(shù);采用垂直布置比水平布置更有利于提高體積傳質(zhì)系數(shù)。一般情況下,氣液兩相的表觀速度可表示為式(20)。

    式中,Ugs、Ugs分別為氣相和液相的表觀速度,m/s;Vg、V1分別為氣相和液相速度,m/s;A為截面面積,m2。

    表2 不同光生物反應(yīng)器內(nèi)表觀氣速和體積傳質(zhì)系數(shù)的比較

    4.3.2 氣含率

    氣含率是指反應(yīng)器內(nèi)氣相所占的體積分?jǐn)?shù),是表征氣液兩相體系最重要的參數(shù)。影響氣含率的因素有氣液界面張力、氣泡尺寸、表觀氣速、黏度等,其中表觀氣速對(duì)氣含率的影響最為顯著[53]。測(cè)量氣含率的方法有光纖法、電導(dǎo)法、壓差法和體積膨脹法。一般情況下,多采用壓差法和體積膨脹法。

    壓差法根據(jù)不同氣含率流體產(chǎn)生的靜壓不同,測(cè)量局部區(qū)域的平均氣含率,見式(21)。

    式中,Δh為壓差計(jì)的讀數(shù);ΔH為相鄰測(cè)點(diǎn)的間距,m。

    4.3.3 黏度

    黏度是影響藻液流動(dòng)狀態(tài)的重要參數(shù),通常認(rèn)為光自養(yǎng)情況下黏度與水接近[54]。清華大學(xué)吳良柏等[25]分析小球藻藻液濃度為C(C=2.7×107個(gè)/m L)和C/2的黏度,得到黏度隨濃度的增加而增加,都比純水大,但差異較小,認(rèn)為小球藻流體為牛頓流體。Petkov等[55]分析尖柵列藻屬、綠球藻和porphyridium sordidum 在不同濃度和溫度下的黏度,得到的黏度范圍在0.82~2.62 mPa·s。華東理工大學(xué)劉晶璘等[56]分析鈍頂螺旋藻的流變特性,發(fā)現(xiàn)其藻液屬冪律流體,隨著藻濃度的增加,藻液的稠度系數(shù)增大,藻液變黏稠。

    5 結(jié) 語

    微藻培養(yǎng)光生物反應(yīng)器內(nèi)傳遞現(xiàn)象涉及光學(xué)、生物學(xué)、化學(xué)和熱物理等多個(gè)學(xué)科,問題尤為復(fù)雜。作者認(rèn)為今后對(duì)于規(guī)模化微藻培養(yǎng)光生物反應(yīng)器的研究應(yīng)從以下幾個(gè)方面進(jìn)行開展。

    (1)結(jié)合最佳微藻生長(zhǎng)影響因子,分析培養(yǎng)液流動(dòng)、傳質(zhì)及光生物反應(yīng)器結(jié)構(gòu)對(duì)各傳遞現(xiàn)象參數(shù)的影響,提煉培養(yǎng)液氣液傳質(zhì)及流動(dòng)特性對(duì)微藻生長(zhǎng)的影響機(jī)制。

    (2)基于現(xiàn)代實(shí)驗(yàn)測(cè)試技術(shù),系統(tǒng)認(rèn)識(shí)不同細(xì)胞濃度下藻液流體的相關(guān)物性參數(shù)。

    (3)探索各傳遞現(xiàn)象在規(guī)模放大中變化規(guī)律,設(shè)計(jì)適用于微藻規(guī)模培養(yǎng)高效光生物反應(yīng)器。

    [1] 李元廣,譚天偉,黃英明.微藻生物柴油產(chǎn)業(yè)化技術(shù)中的若干科學(xué)問題及其分析[J].中國基礎(chǔ)科學(xué),2009,11(5):64-70.

    [2] Ugwu C U,Aoyagi H,Uchiyama H. Phtobioreactors for mass cultivation of algae[J].Bioresource Technology,2008,99(10):4021-4028.

    [3] Xu L,Weathers P J,Xiong X R,el al. M icroalgal bioreactors:Challenges and opportunities[J].Engineering in Life Sciences,2009,9(3):178-189.

    [4] 劉晶璘,張嗣良.封閉式光生物反應(yīng)器研究進(jìn)展[J].生物工程學(xué)報(bào),2000,16(2):119-123.

    [5] 程麗華,張林,陳歡林,等.微藻固定CO2研究進(jìn)展[J].生物工程學(xué)報(bào),2005,21(2):177-181.

    [6] Grima E M,Acien Fernandez F G,Camacho F G,et al. Photobioreactor:Light regime,mass transfer and scale up[J].Journal of Biotechnology,1999,70(1-3):231-247.

    [7] 劉春朝,劉瑞,王鋒.微藻培養(yǎng)過程的光特性研究進(jìn)展[J]. 生物加工過程,2011,9(6):69-76,

    [8] Acien Fernandez F G,Camacho F G,Sanchez Perez J A,et al. A model for light distribution and average solar irradiance inside outdoor tubular photobioreactors for the m icroalgal mss culture[J].Biotechnology and Bioengineering,1997,55(5):701-714.

    [9] Cornet J F,Dussap C G,Dubertret G. A structure model for simulation of cultures of the cyanobacterium spirulina platereis in photobioreactorsⅠ:Coupling between light transfer and grow th kinetics[J].Biotechnology and Bioengineering,1992,40(7):817-825.

    [10] Hsieh C H,Wu W T. A novel photobioreactor w ith transparent rectangular chambers for cultivation of m icroalgae[J].Biochemical Engineering Journal,2009,46(3):300-305..

    [11] Eduardo J L,Carlos Henrique G S,Lucy F L,et al. Effect of light cycles(night/day) on CO2fixation and biomass production by m icroalgae in photobioreactor[J].Chemical Engineering and Processing:Process Intensification,2009,48(1):306-310.

    [12] Meireles L A. On-line control of light intensity in a m icroalgal bioreactor using a novel automatic system[J].Enzyme and Microbial Technology,2008,7(5):554-559.

    [13] Richmode A. Principles for attaining maximal m icroalgal productivity in photobioreactors:An overview[J].Hydrobiologia,2004,512(1-3):33-37.

    [14] Lee C G. Calculation of light penetration depth in photobioreactor[J].Biotechnology and Bioprocess Engineering,1999,4(1) :78-81.

    [15] 康瑞娟,蔡昭鈴,施定基.光強(qiáng)在魚腥藻7120培養(yǎng)液中的衰減及其對(duì)藻細(xì)胞生長(zhǎng)的影響[J].化工冶金,2000,21(4):384-388.

    [16] 康瑞娟,周文齊,蔡昭鈴,等.聚球藻7002在光生物反應(yīng)器中的光自養(yǎng)培養(yǎng)[J].生物工程學(xué)報(bào),2000,16(5):618-622.

    [17] 陳燁,駱其君,馬斌,等.角毛藻光衰減特性及其在光生物反應(yīng)器培養(yǎng)中的應(yīng)用[J].寧波大學(xué)學(xué)報(bào):理工版,2004,17(2):143-146.

    [18] 孫穎穎,王懷忠,閻斌倫.海洋微藻培養(yǎng)中的光衰減研究[J].淮海工學(xué)院學(xué)報(bào):自然科學(xué)版,2008,17(2):62-65.

    [19] 溫少紅.螺旋藻培養(yǎng)液光衰減規(guī)律的研究[J].海洋通報(bào),2001,20(4):92-96.

    [20] 溫少紅,王長(zhǎng)海,鞠寶,等.紫球藻培養(yǎng)液光衰減規(guī)律的研究[J].海洋通報(bào),2001,20(2):35-39.

    [21] 朱篤,李元廣,葉勤,等.轉(zhuǎn)基因聚球藻7942培養(yǎng)中光傳遞及對(duì)細(xì)胞生長(zhǎng)的影響[J].過程工程學(xué)報(bào),2003,3(3):231-237.

    [22] 徐明芳,周遠(yuǎn)志,區(qū)德洪.光生物反應(yīng)器中光衰減特征與螺旋藻生長(zhǎng)動(dòng)力學(xué)研究[J].海洋學(xué)報(bào),2001,25(11):32-37.

    [23] 張志斌.微藻光自養(yǎng)培養(yǎng)中光傳遞及其生長(zhǎng)動(dòng)力學(xué)研究[D].南昌:江西師范大學(xué),2008,15-24.

    [24] 劉晶璘,張嗣良.光生物反應(yīng)器微藻生長(zhǎng)與比消光系數(shù)的關(guān)系[J].化學(xué)反應(yīng)工程與工藝,1999,15(4):407-414.

    [25] 吳良柏.小球藻流體熱物性及其培養(yǎng)條件控制的研究[D].北京:清華大學(xué),2010,15-20.

    [26] Shang H,Scott J A,Shepherd S H,Ross G M. A dynamic thermal model for heating microalgae incubator ponds using off-gas[J].Chemical Engineering Science.,2010,65(16):4591-4597.

    [27] Quentin B,Andy S,F(xiàn)ringer O B,et al. Mechanistic modeling of broth temperature in outdoor photobioreactors[J].Environ.Sci.Technol.,2010,44(6):2197-2203.

    [28] Richmod A. Principles for attaining maximal m icroalgal productivity in photobioreactors:An overview[J].Hydrobiologia,2004,512(1-3):33-37.

    [29] Gutierrez J,Porta-Gandara M A,et al. Passive temperature solar control of an outdoor photobioreactor[J].Renewable Energy,2008,33(8):1892-1903.

    [30] Goetz V,F(xiàn)Beorqne F L,Pruvost J,et al. A generic temperature model for solar photobioreactors[J].Chemical Engineering Journal,2011, 175(15):443-449.

    [31] 王爽,姜秀民.大型海藻生物質(zhì)的比熱容實(shí)驗(yàn)研究[J].太陽能學(xué)報(bào),2009,30(10):1356-1360.

    [32] Chisti M Y. Airlift Bioreactors[M]. Elsevier Appiled Science,1989:13-15.

    [33] Rubio F C,Acien Fernandz F G,et al. Prediction of dissolved oxygen and carbon dioxide concentration profiles in tubular photobioreactors for m icroalgal culture[J].Biotechnology and Bioengineering,1999,62(1):71-86.

    [34] 程桂林,程麗華,周成旭,等.光生物反應(yīng)器脫除空氣中CO2的模型研究[J].生物工程學(xué)報(bào),2006,22(5):821-828.

    [35] 趙東勝,劉桂敏,趙艷麗,等. 氣升式反應(yīng)器研究進(jìn)展[J].化工進(jìn)展,2007,26(6):810-813.

    [36] 張峰,景文珩,邢衛(wèi)紅,等.外循環(huán)氣升式膜反應(yīng)器中的體積傳質(zhì)系數(shù)[J].化學(xué)工程,2008,36(11):5-16.

    [37] Doucha J,Livansky K. Productivity CO2/O2exchange and hydraulics in outdoor open high density m icroalgal (Chlorella. Sp) photobioreactors operated in a M iddle and South European climate[J].Journal of Appllied Phycology,2006,18(6):811-826.

    [38] Huang Q S,Chao Y,Mao G Z,et al. CFD simulation of hydrodynam ics and mass trnasfer in an internal airlift loop reactor using a steady two-fluid model[J].Chemical Engineering Science,2010,65(20):5527-5536.

    [39] 張嗣良.工業(yè)生物過程優(yōu)化與放大研究中的科學(xué)問題[J].中國基礎(chǔ)科學(xué),2009,11(5):27-30.

    [40] 銀建中,程紹杰,賈凌云,等.生物反應(yīng)器放大因素與方法研究[J].化工裝備技術(shù),2009,30(1):22-27.

    [41] Mudde R F,Van Den Akker H E A. 2D and 3D simulation of an internal airlift loop reactor on the basis of a two-fluid model[J].Chemical Engineering Science,2001,56(21-22):6351-6358.

    [42] Kerdouss F,Bnnari A,Proulx P. CFD modeling of gas dispersion and bubble size in a double turbine stirred tank[J].Chemical Engineering Science,2006,61(10):3313-3322.

    [43] Pruvost J,Pottier L,Legrand J. Numerical investigation of hydrodynam ic and m ixing conditions in a torus photobioreactor[J].Chemical Engineering Science,2006,61(14):4476-4489.

    [44] Iris P N,Clemens P. Simulations of light intensity variation in photobioreactors[J].Journal of Biotechnology,2007,131(3):276-285.

    [45] 林晨,李元廣,等.一種新型多節(jié)隔板-平板式光生物反應(yīng)器的數(shù)值和實(shí)驗(yàn)研究[J].高校化工程學(xué)報(bào),2009,23(2):263-269.

    [46] Sato T,Yamada D,Hirabayashi S. Development of virtual photobioreactor for m icroalgae culture considering turbulent flow and flashing light effect[J].Energy Conversion and Management,2010,51(6):1196-1201.

    [47] 蔡昭鈴,片岡廣.氣升式反應(yīng)器內(nèi)高粘性介質(zhì)中氣液傳質(zhì)的研究[J].化學(xué)反應(yīng)工程與工藝,1991,7(2):136-141.

    [48] Velarde R R,Urbina E C,et al. Hydrodynam ic and mass transfer characterization of a flat-panel airlift photobioreactor w ith high light path[J].Chemical Engineering and Processing,2010,49(1):97-103.

    [49] M iron A S,Camacho F G,et al. Bubble-column and airlift photobioreactors for algal culture[J].AIChE Journal,2000,46(9):1872-1887.

    [50] Sierra E,Acien F G,et al. Characterization of a flat plate photobioreactor for the production of m icroalgae[J].Chem.Eng.Journal,2008,138(1-3):136-147.

    [51] Contreras A,Garcia F,et al. Interaction between CO2-mass transfer,light availability,and bhdrodynam ic stress in the grow th of Phaeodactylum tricornum tum in a concentric tube airlift photobioreactor[J].Biotechnol.Bioeng.,1998,60(3):317-325.

    [52] M iyamoto K,Wable O,Benenmann J R. Vertical tubular reactor for microalgae cultivation[J].Biotechnology Letters,1999,10(10):703-708.

    [53] 何廣湘,楊索和,靳海波.鼓泡床反應(yīng)器內(nèi)流動(dòng)與傳質(zhì)行為的研究進(jìn)展[J].化學(xué)工業(yè)與工程,2007,24(1):75-80.

    [54] 孔維利,陳劍佩,李元廣.能源微藻敞開式光生物反應(yīng)器增設(shè)內(nèi)構(gòu)件CFD研究[J].化工進(jìn)展,2010,29(s1):107-112.

    [55] Petkov G D,Bratkova S G. Viscosity of algal cultures and estimation of turbulency in devices for the mass culture of m icroalgae[J].Algological Studies,1996,114(81):99-104.

    [56] 劉晶璘,李元廣,張嗣良.用管式光生物反應(yīng)器培養(yǎng)螺旋藻的研究[J].生物工程學(xué)報(bào),1999,15(4):525-528.

    Research progress of transport phenomena within photobioreactor for m icroalgae culturing

    CHEN Zhijie1,JIANG Zeyi1,2,ZHANG Xinxin1,3,ZHANG Xinru1
    (1School of Mechanical Engineering,2Beijing Engineering Research Center for Energy Saving and Environmental Protection,3Beijing Key Laboratory for Energy Saving and Em ission Reduction of Metallurgical Industry,University of Science and Technology Beijing,Beijing 100083,China)

    Transport phenomena w ithin photobioreactor has a great impact on the grow th rate and biomass production for mass cultivation of m icroalgae. This paper focuses on the transport phenomena w ithin photobioreactor,including light distribution,heat transfer,mass and momentum transfer,and the research progress of mathematical models of each phenomenon. In this review,the effect of the type and scale of photobioreactor on light distribution,mass transfer and other phenomena is discussed,and the important parameters of each phenomenon,such as light absorption coefficient and volumetric mass transfer coefficient are summarized. The review would provide reference for the design,optim ization and scale-up of high efficiency photobioreactor suitable for mass cultivation of microalgae.

    m icroalgae; photobioreactor; heat transfer; mass transfer; light attenuation

    Q 819;TK 6

    A

    1000–6613(2012)–07–1407–08

    2012-02-06;修改稿日期:2012-02-20。

    中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(FRF-AS-10-005B)。

    陳智杰(1980—),男,博士研究生。聯(lián)系人:姜澤毅,副教授。E-mail zyjiang@ustb.edu.cn。

    猜你喜歡
    生物
    生物多樣性
    生物多樣性
    上上生物
    發(fā)現(xiàn)不明生物
    史上“最黑暗”的生物
    軍事文摘(2020年20期)2020-11-28 11:42:50
    第12話 完美生物
    航空世界(2020年10期)2020-01-19 14:36:20
    最初的生物
    自然生物被直銷
    清晨生物初歷直銷
    生物的多樣性
    曰老女人黄片| 婷婷色综合大香蕉| 制服人妻中文乱码| av卡一久久| 欧美 亚洲 国产 日韩一| 亚洲美女搞黄在线观看| 国产xxxxx性猛交| 成人免费观看视频高清| 国产一区二区激情短视频 | 又黄又粗又硬又大视频| 青春草国产在线视频| 热re99久久精品国产66热6| 午夜免费观看性视频| 中文字幕人妻丝袜制服| 久久久久视频综合| 啦啦啦视频在线资源免费观看| 日韩一区二区三区影片| 欧美人与善性xxx| 色94色欧美一区二区| 波多野结衣一区麻豆| 久久影院123| 水蜜桃什么品种好| 亚洲精品一二三| 精品人妻偷拍中文字幕| 国产亚洲最大av| 午夜激情av网站| 韩国精品一区二区三区| 欧美激情高清一区二区三区 | 免费观看无遮挡的男女| 九九爱精品视频在线观看| 国产在线视频一区二区| 亚洲精品日韩在线中文字幕| 亚洲av电影在线观看一区二区三区| 国产精品99久久99久久久不卡 | 美女大奶头黄色视频| 亚洲国产欧美在线一区| 丝袜脚勾引网站| 高清在线视频一区二区三区| 亚洲人成77777在线视频| 26uuu在线亚洲综合色| 久久青草综合色| 国产无遮挡羞羞视频在线观看| 欧美老熟妇乱子伦牲交| 亚洲美女黄色视频免费看| 国产熟女午夜一区二区三区| 1024视频免费在线观看| 中国国产av一级| 国产视频首页在线观看| 免费观看无遮挡的男女| 免费观看无遮挡的男女| 丰满饥渴人妻一区二区三| 欧美97在线视频| 黑丝袜美女国产一区| 亚洲人成电影观看| 色哟哟·www| 国产片特级美女逼逼视频| 在线观看人妻少妇| 在线天堂中文资源库| 久久免费观看电影| 大香蕉久久成人网| 成人国语在线视频| 18在线观看网站| 午夜福利在线免费观看网站| 80岁老熟妇乱子伦牲交| 在线观看一区二区三区激情| 性高湖久久久久久久久免费观看| 亚洲一区二区三区欧美精品| 日韩大片免费观看网站| 香蕉国产在线看| 欧美精品一区二区免费开放| 美女国产视频在线观看| 黄色怎么调成土黄色| 午夜老司机福利剧场| 美女国产高潮福利片在线看| 欧美亚洲日本最大视频资源| 水蜜桃什么品种好| 国产精品蜜桃在线观看| 久久精品国产亚洲av天美| 久久婷婷青草| 老汉色av国产亚洲站长工具| 久久精品人人爽人人爽视色| 成人亚洲精品一区在线观看| 成人毛片a级毛片在线播放| 香蕉国产在线看| 国产一区亚洲一区在线观看| 男女国产视频网站| 久久久久人妻精品一区果冻| 亚洲国产最新在线播放| 国产深夜福利视频在线观看| 成人毛片a级毛片在线播放| 国产精品 欧美亚洲| 亚洲精品第二区| 国产精品不卡视频一区二区| 亚洲美女视频黄频| 五月天丁香电影| 一区二区日韩欧美中文字幕| 欧美bdsm另类| 日韩av免费高清视频| 黄片播放在线免费| 久久久久久久大尺度免费视频| 91精品三级在线观看| 国产熟女午夜一区二区三区| 国产精品久久久久成人av| 国产又爽黄色视频| 久久久亚洲精品成人影院| 国产人伦9x9x在线观看 | 寂寞人妻少妇视频99o| av网站在线播放免费| a 毛片基地| 久久久欧美国产精品| 91国产中文字幕| av在线播放精品| 乱人伦中国视频| 精品国产乱码久久久久久男人| 在线观看三级黄色| 老汉色∧v一级毛片| 久久久国产精品麻豆| 国产一区二区在线观看av| 高清欧美精品videossex| 成年女人在线观看亚洲视频| 五月伊人婷婷丁香| 色婷婷久久久亚洲欧美| a级片在线免费高清观看视频| 波多野结衣一区麻豆| 丝袜美足系列| 80岁老熟妇乱子伦牲交| 一区二区av电影网| 777久久人妻少妇嫩草av网站| 国产伦理片在线播放av一区| 精品国产一区二区久久| 水蜜桃什么品种好| 国产一区二区激情短视频 | 免费看不卡的av| 日本av免费视频播放| 91精品伊人久久大香线蕉| 午夜福利视频在线观看免费| 一区二区三区精品91| 国精品久久久久久国模美| 高清在线视频一区二区三区| 超色免费av| 一二三四中文在线观看免费高清| 精品一区二区三卡| 哪个播放器可以免费观看大片| 在线观看免费高清a一片| 久久久a久久爽久久v久久| 男女啪啪激烈高潮av片| 国产精品久久久久久av不卡| 考比视频在线观看| 国产极品粉嫩免费观看在线| 最近2019中文字幕mv第一页| 性高湖久久久久久久久免费观看| 激情五月婷婷亚洲| 精品人妻一区二区三区麻豆| 黄片无遮挡物在线观看| 久久国内精品自在自线图片| 99re6热这里在线精品视频| 国产一区二区在线观看av| 男人爽女人下面视频在线观看| 18+在线观看网站| 欧美人与性动交α欧美精品济南到 | 婷婷色综合www| 精品人妻在线不人妻| 国产精品一国产av| 五月天丁香电影| 午夜福利,免费看| 一区二区av电影网| 青春草国产在线视频| 久久久久久久久久久免费av| 自拍欧美九色日韩亚洲蝌蚪91| 一边摸一边做爽爽视频免费| 欧美激情高清一区二区三区 | 黄网站色视频无遮挡免费观看| 日韩熟女老妇一区二区性免费视频| 中文精品一卡2卡3卡4更新| 麻豆乱淫一区二区| 在线天堂中文资源库| 中文字幕色久视频| 晚上一个人看的免费电影| 丝袜美腿诱惑在线| 欧美少妇被猛烈插入视频| 纯流量卡能插随身wifi吗| 免费观看在线日韩| 日日摸夜夜添夜夜爱| 欧美精品高潮呻吟av久久| 日产精品乱码卡一卡2卡三| av卡一久久| 午夜免费观看性视频| 精品国产乱码久久久久久小说| av不卡在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 最近的中文字幕免费完整| 亚洲精品美女久久久久99蜜臀 | 国产极品粉嫩免费观看在线| 久久久久精品性色| 免费高清在线观看日韩| 青青草视频在线视频观看| 精品人妻一区二区三区麻豆| 五月伊人婷婷丁香| 国产日韩欧美视频二区| 久久女婷五月综合色啪小说| 久久综合国产亚洲精品| 欧美精品高潮呻吟av久久| 国产精品无大码| 卡戴珊不雅视频在线播放| 国产精品蜜桃在线观看| 最新的欧美精品一区二区| 色播在线永久视频| 高清视频免费观看一区二区| 丰满乱子伦码专区| 妹子高潮喷水视频| av线在线观看网站| 黄频高清免费视频| 制服诱惑二区| 好男人视频免费观看在线| 久久人人爽人人片av| 尾随美女入室| 精品久久蜜臀av无| 韩国高清视频一区二区三区| 免费高清在线观看日韩| 多毛熟女@视频| 欧美 日韩 精品 国产| 美女脱内裤让男人舔精品视频| 亚洲精品国产一区二区精华液| 亚洲精品成人av观看孕妇| 国产一区二区激情短视频 | 女性生殖器流出的白浆| 亚洲内射少妇av| 成人二区视频| 男的添女的下面高潮视频| 久久久久久久精品精品| 国产精品香港三级国产av潘金莲 | 国产免费福利视频在线观看| 啦啦啦在线免费观看视频4| 亚洲欧美色中文字幕在线| 欧美亚洲 丝袜 人妻 在线| 大片免费播放器 马上看| h视频一区二区三区| 欧美人与性动交α欧美精品济南到 | 国产亚洲最大av| 天天躁夜夜躁狠狠躁躁| av在线观看视频网站免费| 9191精品国产免费久久| 18禁国产床啪视频网站| 精品午夜福利在线看| 丰满饥渴人妻一区二区三| 免费高清在线观看日韩| 9色porny在线观看| 美女午夜性视频免费| 只有这里有精品99| 中文乱码字字幕精品一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 97精品久久久久久久久久精品| 久久久久国产网址| 亚洲精品国产色婷婷电影| 韩国精品一区二区三区| 青草久久国产| 日韩免费高清中文字幕av| 亚洲欧美成人精品一区二区| 在线观看三级黄色| 亚洲美女搞黄在线观看| 大香蕉久久成人网| 亚洲精品日本国产第一区| 捣出白浆h1v1| 欧美精品av麻豆av| 国产免费又黄又爽又色| 久久久精品国产亚洲av高清涩受| 亚洲精品第二区| 亚洲人成77777在线视频| 黑人猛操日本美女一级片| 免费黄网站久久成人精品| 母亲3免费完整高清在线观看 | a级片在线免费高清观看视频| 亚洲av中文av极速乱| 天天躁夜夜躁狠狠久久av| 叶爱在线成人免费视频播放| 狠狠精品人妻久久久久久综合| 精品一品国产午夜福利视频| 亚洲成色77777| 超色免费av| 亚洲国产欧美在线一区| 亚洲色图综合在线观看| 美国免费a级毛片| 最近2019中文字幕mv第一页| 国产精品蜜桃在线观看| 在线观看人妻少妇| 九草在线视频观看| 丝袜美腿诱惑在线| 五月伊人婷婷丁香| 伦理电影免费视频| 不卡av一区二区三区| 亚洲国产成人一精品久久久| 国产淫语在线视频| 亚洲欧美成人精品一区二区| 丝袜美腿诱惑在线| 成人黄色视频免费在线看| 日韩精品有码人妻一区| 哪个播放器可以免费观看大片| 男女免费视频国产| 欧美精品一区二区免费开放| 少妇人妻久久综合中文| 91久久精品国产一区二区三区| 我的亚洲天堂| 丝袜脚勾引网站| 在线亚洲精品国产二区图片欧美| 天天影视国产精品| 欧美+日韩+精品| 99国产精品免费福利视频| 亚洲精品美女久久久久99蜜臀 | 日本av免费视频播放| 亚洲美女黄色视频免费看| 色网站视频免费| 母亲3免费完整高清在线观看 | 国产精品女同一区二区软件| 啦啦啦啦在线视频资源| 人人妻人人爽人人添夜夜欢视频| 欧美精品一区二区免费开放| 亚洲 欧美一区二区三区| 久久精品久久久久久久性| a级片在线免费高清观看视频| 免费人妻精品一区二区三区视频| 激情五月婷婷亚洲| 一级a爱视频在线免费观看| 女的被弄到高潮叫床怎么办| 国产又色又爽无遮挡免| 欧美精品一区二区大全| 五月伊人婷婷丁香| 久久久精品国产亚洲av高清涩受| 美女大奶头黄色视频| 久久久亚洲精品成人影院| videosex国产| 国产黄频视频在线观看| xxx大片免费视频| 中文字幕最新亚洲高清| 久久国产亚洲av麻豆专区| 夜夜骑夜夜射夜夜干| 韩国高清视频一区二区三区| 老女人水多毛片| 亚洲人成网站在线观看播放| 美女高潮到喷水免费观看| 亚洲成人av在线免费| 精品久久蜜臀av无| 在线观看一区二区三区激情| 亚洲av.av天堂| 中文字幕制服av| 91国产中文字幕| 亚洲欧美一区二区三区国产| a级毛片在线看网站| 在现免费观看毛片| 免费观看在线日韩| 热99国产精品久久久久久7| 国产高清不卡午夜福利| 丝袜美足系列| 波野结衣二区三区在线| 王馨瑶露胸无遮挡在线观看| 成人亚洲精品一区在线观看| 精品少妇内射三级| 伦理电影大哥的女人| 午夜免费鲁丝| 久久久久精品人妻al黑| 久久97久久精品| 国产在线视频一区二区| av在线播放精品| 在线精品无人区一区二区三| 91在线精品国自产拍蜜月| 午夜福利,免费看| 中国国产av一级| 久久久欧美国产精品| 国产精品免费大片| 午夜福利一区二区在线看| 高清av免费在线| 久久久久久久亚洲中文字幕| 精品久久蜜臀av无| 成年女人在线观看亚洲视频| 有码 亚洲区| 精品一品国产午夜福利视频| 丰满乱子伦码专区| 欧美激情 高清一区二区三区| 国产精品欧美亚洲77777| 精品国产一区二区三区四区第35| 午夜激情久久久久久久| 美女中出高潮动态图| av又黄又爽大尺度在线免费看| 精品亚洲成国产av| 一本色道久久久久久精品综合| 日日啪夜夜爽| 国产精品亚洲av一区麻豆 | 你懂的网址亚洲精品在线观看| 丰满少妇做爰视频| 叶爱在线成人免费视频播放| 春色校园在线视频观看| 高清黄色对白视频在线免费看| av不卡在线播放| 叶爱在线成人免费视频播放| 亚洲av免费高清在线观看| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久成人av| 精品国产露脸久久av麻豆| av电影中文网址| 叶爱在线成人免费视频播放| 免费日韩欧美在线观看| 天天躁夜夜躁狠狠久久av| 亚洲欧美精品综合一区二区三区 | 超碰成人久久| 国产成人精品久久二区二区91 | 热re99久久精品国产66热6| 国产精品国产三级专区第一集| 侵犯人妻中文字幕一二三四区| 精品人妻偷拍中文字幕| 亚洲伊人久久精品综合| 爱豆传媒免费全集在线观看| 国产福利在线免费观看视频| 国产极品天堂在线| 爱豆传媒免费全集在线观看| 国产成人a∨麻豆精品| 亚洲国产精品999| 男女高潮啪啪啪动态图| 国产野战对白在线观看| 久久久久久伊人网av| 高清在线视频一区二区三区| 亚洲美女视频黄频| 欧美日韩亚洲高清精品| 日本午夜av视频| 三级国产精品片| 男女下面插进去视频免费观看| 精品卡一卡二卡四卡免费| 一区福利在线观看| 免费播放大片免费观看视频在线观看| 香蕉精品网在线| 亚洲成色77777| 亚洲成人手机| 国产在线视频一区二区| 色视频在线一区二区三区| 国产又色又爽无遮挡免| 国产精品免费大片| 精品人妻偷拍中文字幕| 精品少妇内射三级| 国产白丝娇喘喷水9色精品| 国产综合精华液| 亚洲欧美成人精品一区二区| 欧美中文综合在线视频| 韩国高清视频一区二区三区| 成人亚洲欧美一区二区av| 黄色怎么调成土黄色| 日韩视频在线欧美| 亚洲国产av影院在线观看| 另类精品久久| 天天操日日干夜夜撸| 欧美日韩亚洲高清精品| 日韩熟女老妇一区二区性免费视频| 欧美国产精品一级二级三级| 国产野战对白在线观看| 又大又黄又爽视频免费| 超碰97精品在线观看| 久久免费观看电影| 春色校园在线视频观看| 性高湖久久久久久久久免费观看| 天堂8中文在线网| 蜜桃在线观看..| 国产精品免费大片| 免费少妇av软件| 国产成人免费观看mmmm| 国产亚洲av片在线观看秒播厂| 国产男女内射视频| 亚洲,欧美,日韩| 精品少妇内射三级| 人成视频在线观看免费观看| 黄网站色视频无遮挡免费观看| 久久精品久久久久久噜噜老黄| 一个人免费看片子| 国产亚洲最大av| 国产一区二区激情短视频 | 久久久久久久久久人人人人人人| www日本在线高清视频| 日韩av在线免费看完整版不卡| 国产无遮挡羞羞视频在线观看| 久久久精品国产亚洲av高清涩受| 一区福利在线观看| 免费观看性生交大片5| 国产熟女欧美一区二区| 少妇熟女欧美另类| 成人国产麻豆网| 国产色婷婷99| 在线观看美女被高潮喷水网站| 老司机影院毛片| 亚洲精品日本国产第一区| 91国产中文字幕| 国产乱人偷精品视频| 日韩av免费高清视频| 亚洲国产精品999| 精品卡一卡二卡四卡免费| 夫妻性生交免费视频一级片| 久久韩国三级中文字幕| 日韩,欧美,国产一区二区三区| 亚洲精品国产一区二区精华液| 在线亚洲精品国产二区图片欧美| 哪个播放器可以免费观看大片| 国产精品蜜桃在线观看| 五月伊人婷婷丁香| 飞空精品影院首页| 亚洲欧美精品自产自拍| 精品少妇内射三级| 国产探花极品一区二区| 777米奇影视久久| 午夜影院在线不卡| 最近最新中文字幕免费大全7| 天天躁夜夜躁狠狠躁躁| 哪个播放器可以免费观看大片| 老汉色∧v一级毛片| 国产一区二区三区av在线| 少妇熟女欧美另类| 大香蕉久久成人网| 亚洲精品在线美女| 在线观看免费日韩欧美大片| 人妻系列 视频| 久久精品亚洲av国产电影网| 麻豆精品久久久久久蜜桃| 亚洲精品久久久久久婷婷小说| 久久人人爽人人片av| 欧美亚洲 丝袜 人妻 在线| 麻豆精品久久久久久蜜桃| 免费人妻精品一区二区三区视频| 建设人人有责人人尽责人人享有的| av网站免费在线观看视频| 最近2019中文字幕mv第一页| 日韩av在线免费看完整版不卡| 麻豆乱淫一区二区| 亚洲av国产av综合av卡| 高清欧美精品videossex| 精品国产一区二区三区久久久樱花| 国产1区2区3区精品| 80岁老熟妇乱子伦牲交| 少妇人妻 视频| 精品国产超薄肉色丝袜足j| 久久久国产精品麻豆| 免费av中文字幕在线| 国产亚洲av片在线观看秒播厂| 性少妇av在线| 久久精品国产a三级三级三级| 久久精品亚洲av国产电影网| 水蜜桃什么品种好| 人成视频在线观看免费观看| a级毛片黄视频| 纵有疾风起免费观看全集完整版| 亚洲精品日本国产第一区| 精品国产乱码久久久久久男人| 国产成人免费无遮挡视频| 精品人妻熟女毛片av久久网站| 男女边吃奶边做爰视频| 天天影视国产精品| 七月丁香在线播放| 伊人久久国产一区二区| 又大又黄又爽视频免费| 妹子高潮喷水视频| 如日韩欧美国产精品一区二区三区| 777久久人妻少妇嫩草av网站| 免费高清在线观看日韩| 美国免费a级毛片| 成人亚洲精品一区在线观看| 久久国产亚洲av麻豆专区| 91精品三级在线观看| 亚洲国产精品一区三区| 国产不卡av网站在线观看| 国产视频首页在线观看| 美女中出高潮动态图| 一本—道久久a久久精品蜜桃钙片| 18禁观看日本| 街头女战士在线观看网站| 久久久久久免费高清国产稀缺| 街头女战士在线观看网站| 乱人伦中国视频| 免费不卡的大黄色大毛片视频在线观看| 国产精品不卡视频一区二区| 国产免费一区二区三区四区乱码| 国产高清不卡午夜福利| 亚洲婷婷狠狠爱综合网| 99国产综合亚洲精品| 免费黄色在线免费观看| 久久精品亚洲av国产电影网| 国产无遮挡羞羞视频在线观看| av免费观看日本| 一级爰片在线观看| 丰满迷人的少妇在线观看| 国产成人精品福利久久| 欧美人与性动交α欧美软件| 男女边摸边吃奶| 十八禁高潮呻吟视频| 97在线视频观看| 国产白丝娇喘喷水9色精品| 人妻少妇偷人精品九色| 久久精品国产自在天天线| 巨乳人妻的诱惑在线观看| 精品第一国产精品| 大陆偷拍与自拍| 热99国产精品久久久久久7| 亚洲欧美一区二区三区黑人 | 久久av网站| 亚洲情色 制服丝袜| 欧美+日韩+精品| 日日摸夜夜添夜夜爱| videosex国产| 人妻人人澡人人爽人人| 你懂的网址亚洲精品在线观看| 中文乱码字字幕精品一区二区三区| 90打野战视频偷拍视频| 欧美日韩国产mv在线观看视频| 老司机影院毛片| 最近手机中文字幕大全| 你懂的网址亚洲精品在线观看| 国产精品久久久久成人av| 最近最新中文字幕大全免费视频 | 精品亚洲成a人片在线观看| 又大又黄又爽视频免费| 久久久久久人妻| 久久久精品国产亚洲av高清涩受| 少妇人妻久久综合中文|