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

    正滲透過程中汲取質(zhì)反向滲透研究進(jìn)展

    2015-07-25 03:35:32謝朋張忠國(guó)孫濤吳月吳秋燕李繼定李珊
    化工進(jìn)展 2015年10期
    關(guān)鍵詞:錯(cuò)流滲透壓通量

    謝朋,張忠國(guó),孫濤,吳月,吳秋燕,李繼定,李珊

    (1 輕工業(yè)環(huán)境保護(hù)研究所,北京 100089;2 清華大學(xué)化學(xué)工程系化學(xué)工程國(guó)家重點(diǎn)聯(lián)合實(shí)驗(yàn)室,北京 100084)

    正滲透(forward osmosis,簡(jiǎn)稱FO)是近年膜技術(shù)領(lǐng)域的研究熱點(diǎn)之一,我國(guó)目前有很多學(xué)者從事相關(guān)研究工作。由于具有低能耗、低污染和高回收率等優(yōu)勢(shì)[1-3],正滲透技術(shù)在海水淡化[4-8]、廢水處理[9-12]、飲料加工[13-15]、發(fā)電[2,16-18]等領(lǐng)域均具有良好的應(yīng)用前景。

    正滲透是以選擇性半透膜兩側(cè)溶液的滲透壓差為驅(qū)動(dòng)力,溶液中的水從低滲透壓側(cè)向高滲透壓側(cè)滲透的傳遞過程[19]。根據(jù)常用的非對(duì)稱半透膜的取向不同,正滲透技術(shù)可以分為正滲透模式和壓力阻尼滲透(pressure retarded osmosis,簡(jiǎn)稱PRO)模式[19]。膜的活性層(或分離層)朝向原料液,而支撐層朝向汲取液的操作模式稱為FO 模式;膜的活性層朝向汲取液,而支撐層朝向原料液的操作模式稱為PRO 模式。

    理想的半透膜只允許水分子通過,溶質(zhì)和其他雜質(zhì)則被完全截留。但事實(shí)上,由于膜兩側(cè)溶質(zhì)濃度的差異和膜材料截留能力有限,汲取液的溶質(zhì)(或稱汲取質(zhì))會(huì)透過半透膜進(jìn)入原料液,發(fā)生反向滲透(以下簡(jiǎn)稱“反滲”)現(xiàn)象,造成汲取質(zhì)的流失,并產(chǎn)生嚴(yán)重的濃差極化,影響水通量和污染物截留效果,從而制約了正滲透技術(shù)的應(yīng)用[20-22]。另外,在實(shí)際應(yīng)用過程中,汲取質(zhì)反滲還可能會(huì)對(duì)后續(xù)的一些工藝產(chǎn)生影響。因此,研究汲取質(zhì)反滲現(xiàn)象對(duì)于正滲透技術(shù)的發(fā)展和應(yīng)用尤為重要。本文重點(diǎn)介紹了正滲透過程中FO 模式下溶質(zhì)濃差極化和汲取質(zhì)反滲模型,歸納了汲取質(zhì)反滲過程的影響因素及影響規(guī)律,以供有關(guān)科研人員參考。

    1 FO 過程汲取質(zhì)反滲模型研究進(jìn)展

    在正滲透技術(shù)的發(fā)展歷程中,關(guān)于汲取質(zhì)反滲的研究起步比較晚。1997年,Loeb 等[23]首先提出了正滲透中汲取質(zhì)反滲的概念,并建立了初步的相關(guān)模型。隨后,Tang 等[24]、McCutcheon 等[25]在Loeb 模型基礎(chǔ)上建立了比較完善的濃差極化模型。2010年,Phillip 等[20]建立了更加完善的汲取質(zhì)反滲模型,并通過實(shí)驗(yàn)驗(yàn)證了其有效性,但這一模型僅考慮了內(nèi)濃差極化,而忽略了外濃差極化。2012年,Phillip 等[26]改進(jìn)了其提出的汲取質(zhì)反滲模型,考慮了原料液側(cè)的外濃差極化,但依舊沒有考慮汲取液側(cè)的外濃差極化。2013年,Suh 等[27]進(jìn)一步完善了該模型,全面考慮膜兩側(cè)的外濃差極化,但該模型只適用于FO 模式,不適用于PRO 模式。

    1.1 FO 過程的濃差極化現(xiàn)象

    汲取質(zhì)濃差極化現(xiàn)象不僅影響滲透驅(qū)動(dòng)膜過程的水通量,而且影響汲取質(zhì)的反滲過程。濃差極化模型是汲取質(zhì)反滲模型不可或缺的理論基礎(chǔ)之一。

    滲透驅(qū)動(dòng)膜過程存在兩種不同形式的濃差極化現(xiàn)象,即外濃差極化(external concentration polarization,簡(jiǎn)稱ECP)和內(nèi)濃差極化(internal concentration polarization,簡(jiǎn)稱ICP)[22]。對(duì)于FO過程,原料液中溶質(zhì)被截留在活性層表面附近,形成濃縮的外濃差極化(concentrated ECP,簡(jiǎn)稱cECP);在滲透壓驅(qū)動(dòng)下,水透過活性層,支撐層中的汲取液被稀釋,形成稀釋的內(nèi)濃差極化(diluted ICP,簡(jiǎn)稱dICP)。對(duì)于PRO 過程,原料液中溶質(zhì)被截留在支撐層附近,形成濃縮的內(nèi)濃差極化(concentrated ICP,簡(jiǎn)稱cICP);水透過活性層進(jìn)入汲取液側(cè),活性層表面附近的汲取液被稀釋,形成稀釋的外濃差極化(diluted ECP,簡(jiǎn)稱dECP)。

    圖1 中,(a)和(b)分別為FO 模式和PRO 模式下,汲取質(zhì)在非對(duì)稱性膜中的濃度分布情況,其中考慮了支撐層的內(nèi)濃差極化和膜兩側(cè)表面的外濃差極化。

    對(duì)于圖1(a)所示的FO 過程,非對(duì)稱性膜的活性層朝向原料液,原料液側(cè)的水在滲透壓差的作用下透過活性層進(jìn)入支撐層,稀釋了支撐層中的汲取液,形成稀釋的ICP。ICP 的存在降低了有效滲 透壓差,是FO 過程中影響水通量的最主要因 素[22-23,28-33]。

    在壓力驅(qū)動(dòng)過程中, ECP 是造成膜阻力增大的主要因素;而在正滲透過程中,ECP 對(duì)FO 過程的影響相對(duì)于ICP 要小得多[6,31],因此常被忽 略[20,23-25]。但最近研究發(fā)現(xiàn),在錯(cuò)流速度較低或者水通量較高的情況下,ECP 對(duì)FO 過程影響較大,應(yīng)予以考慮[27]。

    1.2 FO 過程的汲取質(zhì)反滲模型

    1997年,Loeb 等[23]首先提出了正滲透中汲取質(zhì)反滲的概念,并在PRO 模式下建立了汲取質(zhì)反滲的數(shù)學(xué)模型,如式(1)所示。由于未考慮支撐層內(nèi)濃差極化和膜兩側(cè)的外濃差極化現(xiàn)象,所以這一模型同樣適用于FO 模式,如式(1)。

    圖1 非對(duì)稱膜正滲透過程中汲取質(zhì)在膜內(nèi)分布情況

    式中,Js為汲取質(zhì)反滲通量,mol/(m2·h);B 為汲取質(zhì)的滲透系數(shù),m/d。

    FO 過程支撐層中溶質(zhì)反向滲透通量由方向相反的兩部分決定:一部分是鹽濃度梯度引起的分子擴(kuò)散,另一部分是水滲透引起的對(duì)流擴(kuò)散[34]。2010年,Phillip 等[20]在Loeb 模型的基礎(chǔ)上對(duì)汲取質(zhì)反滲模型進(jìn)行了完善,即考慮了支撐層內(nèi)的濃差極化,建立了公式,如式(2)。

    式中,Jw為水通量,L/( m2·h);D 為汲取質(zhì)在主體溶液中的擴(kuò)散系數(shù),m2/s;S 為支撐層的結(jié)構(gòu)參數(shù),可由式(3)表示[20],其測(cè)試方法可參考文獻(xiàn)[35]。

    式中,A 為水的滲透系數(shù),m/(d·Pa)。

    2012年,Phillip 等[26]對(duì)其提出的汲取質(zhì)反滲模型進(jìn)行改進(jìn),將原料液側(cè)的外濃差極化考慮進(jìn)來(lái),但仍未考慮汲取液側(cè)的外濃差極化。改進(jìn)后的公式如式(4)。

    式中,kF為原料液側(cè)汲取質(zhì)的傳質(zhì)系數(shù),m/s;kF與舍伍德數(shù)相關(guān)[26,32],可由式(5)計(jì) 算[5,32,36]。

    式中,Sh 為舍伍德數(shù);δF為原料液側(cè)外濃差極化厚度,m;dh為水力直徑,m。

    舍伍德數(shù)可由式(6)、式(7)計(jì)算。

    式中,Re 為雷諾數(shù);Sc 為施密特?cái)?shù); L 為流動(dòng)通道的長(zhǎng)度,m。

    施密特?cái)?shù)Sc 可由式(8)計(jì)算:

    式中,μ 為水的黏度,Pa·s;ρ 為水的密度,kg/m3。

    2013年,Suh 等[27]考慮了支撐層的內(nèi)濃差極化和膜兩側(cè)的外濃差極化,建立了目前最完善的汲取質(zhì)反滲模型,如式(9)所示。

    式中,K 是溶質(zhì)在支撐層內(nèi)擴(kuò)散的阻力系 數(shù)[17,32];Jsw為選擇性反向汲取質(zhì)通量,mol/L;kD為汲取液側(cè)汲取質(zhì)的傳質(zhì)系數(shù),與舍伍德數(shù)相關(guān)。

    式中,Ds為汲取質(zhì)在支撐層內(nèi)的有效擴(kuò)散系數(shù),m2/s;ε 為支撐層的孔隙率;τ 為支撐層的彎度。其中Ds與D 的關(guān)系式如式(11)[37]。

    由于同時(shí)考慮了內(nèi)外濃差極化的影響,Suh 等建立的FO 模式下的汲取質(zhì)反滲模型更加接近真實(shí)情況,但其公式也更加復(fù)雜。汲取質(zhì)反滲現(xiàn)象主要數(shù)學(xué)模型的公式詳見表1。

    表1 汲取質(zhì)反滲現(xiàn)象的主要數(shù)學(xué)模型

    Tang 等[24]在總結(jié)Loeb 等[23]正滲透模型的基礎(chǔ)上將Jw/Js定義為選擇性反向汲取質(zhì)通量Jsw,即水通量與汲取質(zhì)逆向通量的比值,表示每產(chǎn)生定量的水所損失的汲取質(zhì)的量。選擇性反向汲取質(zhì)通量由活性層的透水性和透鹽性決定,與錯(cuò)流速度和汲取質(zhì)濃度無(wú)關(guān)[20,37-38]。Jsw可用式(12)表示[20]。

    式中,n 為由汲取質(zhì)產(chǎn)生的溶解物種類(如NaCl,n=2);R 為理想氣體常數(shù),J/(mol·K);T 為絕對(duì)溫度,K。該式對(duì)于理解正滲透過程中的汲取質(zhì)反向滲透非常重要,它表明膜的支撐層結(jié)構(gòu)參數(shù)S 和汲取液濃度是影響汲取質(zhì)反向滲透的主要 因素。

    2 汲取質(zhì)反滲過程的影響因素

    依據(jù)前面得出的汲取質(zhì)反滲通量公式(9),可知汲取質(zhì)反滲與汲取質(zhì)的滲透系數(shù)B、汲取質(zhì)傳質(zhì)系數(shù)kD(=ShD/dh)、汲取質(zhì)在支撐層內(nèi)擴(kuò)散的阻力系數(shù)K(=S/D)密切相關(guān)。B 由汲取質(zhì)本身特性決定。膜的結(jié)構(gòu)參數(shù)S,因膜的種類不同而不同,同一種膜一般認(rèn)為S 值穩(wěn)定在一定范圍內(nèi)。K 和kD與汲取質(zhì)擴(kuò)散系數(shù)D 和黏度μ 相關(guān)。綜合各相關(guān)條件可知,汲取質(zhì)的反滲通量受汲取質(zhì)本身滲透系數(shù)、操作溫度、膜材料、初始滲透壓、膜結(jié)構(gòu)和所帶電荷等條件影響。

    2.1 滲透壓差或汲取液濃度的影響

    由式(1)可知,F(xiàn)O 過程中汲取質(zhì)的反滲通量與汲取質(zhì)的滲透系數(shù)、膜活性層兩側(cè)的有效濃度差成正比。同一種汲取質(zhì)的滲透系數(shù)維持在穩(wěn)定范圍內(nèi),但其濃度與滲透壓正相關(guān),所以汲取質(zhì)的反向滲透通量應(yīng)隨滲透壓差的增大而增加。如Suh 等[27]以去離子水為原料液,以NaCl 為汲取質(zhì),采用FO模式進(jìn)行實(shí)驗(yàn),發(fā)現(xiàn)當(dāng)汲取質(zhì)濃度由1 mol/L 逐漸增加到4 mol/L 時(shí),汲取質(zhì)的反滲通量也逐漸 增大。

    大量研究同樣發(fā)現(xiàn),在濃度較低時(shí)汲取質(zhì)的反滲通量隨著滲透壓差的升高而近似線性上升,但是在汲取質(zhì)濃度較高時(shí)汲取質(zhì)的反滲通量上升速度減緩[20,39-40]。這是因?yàn)樵跐B透壓差較高時(shí),水通量較大,由于水分子和汲取質(zhì)分子在活性層內(nèi)的流動(dòng)方向是相反的,高水通量一定程度上限制了汲取質(zhì)的反向滲透[6,24]。

    滲透壓差較低時(shí),雖然能夠保證汲取質(zhì)反滲通量保持在很小的范圍內(nèi),但是也限制了水通量;過高的滲透壓差雖然能夠提高水通量,但也會(huì)使汲取質(zhì)反滲通量增加,從而增加運(yùn)行成本。所以實(shí)際應(yīng)用中,應(yīng)該選擇合理的滲透壓差,既能保證水通量,又不會(huì)造成汲取質(zhì)的大量流失。

    2.2 膜表面錯(cuò)流速度的影響

    汲取液和原料液在膜表面的錯(cuò)流速度對(duì)于汲取質(zhì)反滲也很重要,因其影響膜的外濃差極化[37,41-43]。錯(cuò)流速度與汲取質(zhì)反滲通量呈正相關(guān)性,即錯(cuò)流速度較低時(shí)汲取質(zhì)反滲通量較低,錯(cuò)流速度較高時(shí)汲取質(zhì)反滲通量也較高[21]。以FO 操作模式為例,其原因是當(dāng)錯(cuò)流速度較低時(shí),處于原料液側(cè)邊界層內(nèi)的汲取質(zhì)(自汲取液側(cè)擴(kuò)散而來(lái))受到的錯(cuò)流流動(dòng)擾動(dòng)很小,從而不斷聚集,濃度逐漸升高,而汲取液側(cè)汲取液的稀釋效應(yīng)依然存在,導(dǎo)致膜兩側(cè)滲透壓差降低,進(jìn)而降低水通量和鹽通量;相反,膜兩側(cè)錯(cuò)流速度較高時(shí),活性層表面汲取質(zhì)隨水流走,原料液側(cè)的ECP 減弱,膜兩側(cè)維持較高的滲透壓差,水通量和鹽通量均增大。

    此外,Hancock 等[21]還發(fā)現(xiàn),為了減少汲取質(zhì)因反滲而產(chǎn)生的損耗,膜兩側(cè)宜保持較低的流速,但此時(shí)水通量也較低;汲取液流速較高,而原料液流速較低的操作對(duì)于控制汲取質(zhì)的反滲最為不利;為了降低選擇性半透膜反向汲取質(zhì)通量,膜兩側(cè)宜保持相同的錯(cuò)流速度。膜表面錯(cuò)流速度對(duì)水通量、汲取質(zhì)反滲通量的影響情況詳見圖2[21]。由圖2 可以看出,當(dāng)錯(cuò)流速度較低時(shí),汲取質(zhì)反滲通量和水通量也較低;隨著錯(cuò)流速度的增大,汲取質(zhì)反滲通量和水通量增大。

    2.3 膜結(jié)構(gòu)與膜材料的影響

    目前,常用的正滲透膜為非對(duì)稱性的三乙酸纖維素(cellulose triacetate,簡(jiǎn)稱CTA)膜和薄壁復(fù)合(thin-film composite,簡(jiǎn)稱TFC)膜。以美國(guó)HTI公司生產(chǎn)的正滲透膜為例[44],其CTA 膜由支撐層和活性層構(gòu)成,支撐層為聚酯篩網(wǎng)或聚酯無(wú)紡布結(jié)構(gòu),活性層為乙酸纖維素,膜的厚度為 50~80μm;而TFC 膜結(jié)構(gòu)比較復(fù)雜,作為支撐層的聚酯篩網(wǎng)被膜層包在中間,活性層為聚酰胺或聚砜。受膜材料和結(jié)構(gòu)的影響,這兩種膜對(duì)鹽的截留率不同。在相同的操作條件下,TFC 膜的汲取質(zhì)反滲通量Js較CTA 大,但由于TFC 膜具有更高的水通量,其Js/Jw較CTA 膜小得多,如表2 所示。因此,TFC 膜成為目前非對(duì)稱性正滲透膜最主要的發(fā)展方向。不過,為了降低內(nèi)濃差極化,提高水 通量,一些研究人員嘗試開發(fā)對(duì)稱性正滲透 膜[45,48-49]。但由于不同研究工作之間的正滲透操作條件(如采用的汲取質(zhì),詳見表2)不同,目前 尚難以確定對(duì)稱性正滲透膜在選擇性反向汲取質(zhì)通量方面是否具有優(yōu)勢(shì)。

    圖2 錯(cuò)流速度的影響[21]

    2.4 溫度的影響

    溫度是汲取質(zhì)反滲過程中一個(gè)不可忽視的因素。溫度會(huì)影響膜結(jié)構(gòu)、汲取質(zhì)的擴(kuò)散系數(shù)、溶液黏度等。由FO 模型汲取質(zhì)反滲模型可知,汲取質(zhì)反滲通量主要受支撐層內(nèi)汲取質(zhì)擴(kuò)散的阻力系數(shù)K和汲取質(zhì)的傳質(zhì)系數(shù)k 影響。一般溶液的擴(kuò)散系數(shù)D 隨著溫度的升高而增大且呈指數(shù)變化,由式(10)、式(11)可知擴(kuò)散系數(shù)增大會(huì)導(dǎo)致汲取質(zhì)在支撐層內(nèi)的擴(kuò)散阻力系數(shù)K 減小;溶液的黏度隨溫度的升高而減小,且擴(kuò)散系數(shù)的變化大于黏度的變化,即μ/D 不斷減小,由式(5)、式(8)可知汲取質(zhì)傳質(zhì)系數(shù)k 隨溫度增大而減小。

    表2 膜結(jié)構(gòu)和膜材料對(duì)水通量及汲取質(zhì)反滲的影響

    以NaCl、MgCl2、Na2SO4和MgSO4作為汲取質(zhì),在低濃度時(shí),汲取質(zhì)反滲通量隨溫度的升高基本保持恒定,因此可忽略溫度對(duì)汲取質(zhì)反滲的影響;但在高濃度時(shí),隨著溫度的升高,汲取質(zhì)反滲通量先上升再下降[40]。這是因?yàn)榈蜐舛葧r(shí),溫度對(duì)擴(kuò)散系數(shù)和溶液黏度影響較小,因此阻力系數(shù)K和傳質(zhì)系數(shù)k 變化不大,汲取質(zhì)反滲通量基本恒定。當(dāng)濃度較高時(shí),隨著溫度的升高,汲取質(zhì)的擴(kuò)散系數(shù)增大,溶液黏度降低,更利于汲取質(zhì)的反向滲透;當(dāng)溫度增大到一定程度,汲取質(zhì)的傳質(zhì)系數(shù)受溫度影響降低,膜結(jié)構(gòu)發(fā)生變化,傳遞能力受限,繼續(xù)升高溫度時(shí),汲取質(zhì)反滲通量逐漸降低。不過,在高溫下操作時(shí),膜的壽命會(huì)受到影響。因此,應(yīng)根據(jù)具體情況選擇適宜的操作溫度,以便得到較高的水通量,同時(shí)控制較低的汲取質(zhì)反滲 通量。

    2.5 汲取質(zhì)種類的影響

    以電解質(zhì)為汲取質(zhì)時(shí),水通量和汲取質(zhì)反滲通量都隨著汲取質(zhì)濃度的增加而增大。當(dāng)汲取質(zhì)濃度較低時(shí),水通量和汲取質(zhì)反滲通量與汲取質(zhì)濃度呈近似線性關(guān)系;但汲取質(zhì)濃度過高時(shí),水通量和汲取質(zhì)反滲通量與汲取質(zhì)濃度不再呈線性關(guān)系。這是因?yàn)殡娊赓|(zhì)分子量小,很容易進(jìn)入支撐層中,形成內(nèi)濃差極化,影響水通量和汲取質(zhì)反滲通量。

    以中性溶質(zhì)為汲取質(zhì)時(shí),若其分子比較小,如甘油和赤蘚糖醇[50],水通量和汲取質(zhì)反滲通量隨汲取質(zhì)濃度的變化趨勢(shì)與以電解質(zhì)為汲取質(zhì)時(shí)類似;若是葡萄糖、木糖等[44]分子量比較大的汲取質(zhì)時(shí),在PRO 模式中活性層能夠有效地截留大分子物質(zhì),只有極少數(shù)的汲取質(zhì)通過活性層進(jìn)入支撐層,形成的內(nèi)濃差極化極小,而且此時(shí)水通量較低,外濃差極化不明顯[19],反滲通量隨濃度變化呈線性關(guān)系。

    電解質(zhì)類汲取質(zhì)的擴(kuò)散系數(shù)明顯高于中性汲取質(zhì),但其汲取質(zhì)反滲通量卻更低[39]。這是因?yàn)橐噪娊赓|(zhì)溶液為汲取液時(shí),產(chǎn)生的滲透壓高,水通量大,極大阻礙了汲取質(zhì)的反向滲透。

    在實(shí)際應(yīng)用中,汲取液不可能均是單一的電解質(zhì)溶液或中性汲取質(zhì)溶液,大多時(shí)候是多種電解質(zhì)的混合液、電解質(zhì)與中性汲取質(zhì)的混合液或多種中性汲取質(zhì)的混合液。研究發(fā)現(xiàn)[39],由NaCl 與葡萄糖、甘油等組成的汲取液,增加某一汲取質(zhì)的濃度,另一汲取質(zhì)的反向滲透通量也會(huì)隨著增加。這說明汲取質(zhì)反向傳遞存在著耦合效應(yīng)。

    在不同操作條件下,一些汲取質(zhì)在不同濃度或滲透壓下的反滲通量如表3 所示。表中同時(shí)給出了反滲通量與汲取質(zhì)濃度或滲透壓之間的非線性擬合結(jié)果,其中NaCl、KCl、MgCl2的滲透壓是依據(jù)文獻(xiàn)[51],由其濃度計(jì)算得到的??梢钥闯觯橙≠|(zhì)的反滲通量可以通過汲取液的濃度或滲透壓的一元多項(xiàng)式表達(dá),即

    式中,a1、a2、a3、b1、b2、b3為系數(shù),C 為汲取質(zhì)的濃度,Δπ 為汲取液側(cè)的滲透壓。不過,該式的有效性尚需更多的研究進(jìn)行驗(yàn)證。

    表3 不同汲取質(zhì)不同條件下的反滲通量

    2.6 膜取向的影響

    膜的取向也會(huì)影響汲取質(zhì)的反滲通量。當(dāng)以去離子水或低濃度鹽溶液為原料液時(shí),PRO 模式比FO模式的水通量和汲取質(zhì)反滲通量都要高[19,40,52]。以NaCl 汲取液為例[24],當(dāng)以10 mg/L 對(duì)氨基馬尿酸和0.01 mol/L NaCl 的混合溶液為原料液時(shí),在常溫下,PRO 模式的NaCl 汲取質(zhì)反滲通量要高于FO 模式,如圖3 (a)所示。這是因?yàn)閷?duì)于正滲透過程,影響汲取質(zhì)反滲的主要因素是內(nèi)濃差極化,而PRO 模式中低鹽溶液或去離子水形成的內(nèi)濃差極化要小的多,甚至可忽略,所以PRO 過程主要受外濃差極化影響,外濃差極化形成的膜污染相對(duì)于內(nèi)濃差極化要小很多,因此PRO 模式水通量和汲取質(zhì)反滲通量要高于FO 模式。不過,PRO 模式的Js/Jw要小于FO 模式,如圖3 (b)所示。這是由于PRO 模式的水通量較FO 模式高,而水和汲取質(zhì)的滲透方向相反,因此PRO 模式下相對(duì)較高的水通量會(huì)對(duì)汲取質(zhì)滲透過程產(chǎn)生更大的阻礙作用,從而導(dǎo)致其Js/Jw較FO 模式低。

    圖3 膜取向?qū)橙≠|(zhì)(NaCl)反滲的影響[24]

    實(shí)際應(yīng)用一般宜采用FO 模式。這是因?yàn)檎5脑弦阂话悴豢赡苁侨ルx子水或低濃度鹽溶液,大多會(huì)含有無(wú)機(jī)鹽、有機(jī)物、懸浮物、膠體等物質(zhì),這些物質(zhì)會(huì)在PRO 過程中的支撐層內(nèi)累積,很難通 過錯(cuò)流作用使之疏離,從而導(dǎo)致嚴(yán)重的膜污染,影響水通量和汲取質(zhì)反滲通量[53];而FO 模式下,活性層能有效截留這些物質(zhì),只允許水通過,膜污染很輕,通量穩(wěn)定[10-11,54],從而提高膜的使用壽命。

    2.7 離子水力半徑的影響

    同反滲透、納濾等采用半透膜的膜過程一樣,離子的水力半徑也會(huì)影響其在正滲透膜內(nèi)的遷移過程。Hancock 等[21]以NaCl、MgCl2為汲取質(zhì),分別以CaSO4、K2SO4、MgSO4、Ba(NO3)2為原料液,發(fā)現(xiàn)Na+、Cl-、K+、NH4+、HCO3-等水力半徑較小的離子反向滲透通量均比較大,而 Mg2+、Ba2+、Ca2+等水力半徑較大的離子反向滲透通量相對(duì)較低,其中Na+和Cl-的反滲通量如表4 所示。此外,吳晴[40]分別用NaCl 和MgCl2研究離子反滲規(guī)律,發(fā)現(xiàn)在初始階段,MgCl2溶液中的Cl-通量要大于NaCl 中的Cl-通量,但經(jīng)過一段時(shí)間后,MgCl2溶液中的Cl-通量明顯小于NaCl 中的Cl-通量。這說明水力半徑大的離子對(duì)其他離子的反滲具有抑制作用。根據(jù)唐南膜平衡理論[55-56],溶液體系中總體呈電中性,陰、陽(yáng)離子的反滲是相互作用的,如果陰離子反滲通量大,即使陽(yáng)離子的離子半徑很大,反滲通量也會(huì)相應(yīng)變大,反之亦然。

    離子的水力半徑除了影響汲取質(zhì)的反滲通量,還會(huì)破壞汲取質(zhì)陰、陽(yáng)離子的同步逆向擴(kuò)散過程,導(dǎo)致汲取質(zhì)陰、陽(yáng)離子不按其化學(xué)式中的比例進(jìn)行同步逆向擴(kuò)散。由表4 可知,當(dāng)以NaCl 為汲取質(zhì),以 CaSO4、K2SO4、MgSO4、NH4HCO3、SO2、H3BO3為原料液時(shí),陰、陽(yáng)離子反滲通量相差不大;但當(dāng)以Ba(NO3)2為原料液時(shí),Cl-的反滲通量明顯高于Na+。這主要是因?yàn)樵弦褐蠳O3-離子半徑較其他原料液中的陰離子半徑小,較容易透過半透膜進(jìn)入汲取液側(cè),為維持電離平衡,汲取液側(cè)的Cl-透過半透膜進(jìn)入原料液側(cè)。

    總之,水力半徑較小的一價(jià)離子的反向滲透通量比較高,但由于其滲透性強(qiáng),濃差極化效應(yīng)較小,所以能夠產(chǎn)生較高的水通量;而水力半徑較大的二價(jià)離子的反向滲透通量則比較低,但由于其滲透性相對(duì)較差,濃差極化效應(yīng)較大,影響水通量,且Mg2+、Ca2+等離子易于結(jié)垢。因此,應(yīng)結(jié)合實(shí)際需求,在綜合考慮水通量和汲取質(zhì)反向滲透通量等因素的基礎(chǔ)上,選擇適宜的汲取質(zhì)。

    3 結(jié) 語(yǔ)

    汲取質(zhì)反向滲透是以溶解態(tài)物質(zhì)作為汲取質(zhì)的正滲透過程中難以避免的現(xiàn)象,對(duì)于正滲透技術(shù)的推廣應(yīng)用具有重要的影響。汲取質(zhì)反向滲透通量的數(shù)學(xué)模型雖然不斷發(fā)展,但仍然存在很多缺陷,如數(shù)學(xué)公式表達(dá)形式不夠直觀,需迭代求解等。此外,目前尚缺乏較完善的適用于壓力阻尼滲透(PRO)模式的數(shù)學(xué)模型。通過擬合一些文獻(xiàn)的實(shí)驗(yàn)數(shù)據(jù),發(fā)現(xiàn)無(wú)論是正滲透(FO)模式還是PRO模式,汲取質(zhì)的反滲通量均可以通過其濃度或汲取液滲透壓的一元多項(xiàng)式表達(dá),但各項(xiàng)系數(shù)的物理意義仍不明確。這些均需進(jìn)一步地完善,以便更好地預(yù)測(cè)汲取質(zhì)的反滲通量,從而科學(xué)、有效地評(píng)估其影響。

    表4 FO 操作模式下NaCl 汲取質(zhì)中的Na+和Cl-在不同 條件下的反滲通量(18℃)[21]

    汲取質(zhì)反向滲透過程的影響因素很多,如汲取液濃度、溫度、膜結(jié)構(gòu)與膜材料、汲取質(zhì)的性質(zhì)、離子半徑、膜表面流速、膜取向,等等。這些因素的影響規(guī)律及影響機(jī)制的研究對(duì)于控制汲取質(zhì)反向滲透、選用汲取質(zhì)以及開發(fā)新型汲取質(zhì)和膜材料均具有重要的指導(dǎo)作用,因此亟需進(jìn)一步開展有關(guān) 工作。

    總之,汲取質(zhì)反滲過程的深入研究,將會(huì)極大推動(dòng)新型汲取質(zhì)和膜材料的開發(fā),以及正滲透技術(shù)的應(yīng)用,從而利于解決我國(guó)水環(huán)境污染和水資源短缺等重大問題。

    符 號(hào) 說 明

    A——水滲透系數(shù),m/(d·Pa)

    B——汲取質(zhì)滲透系數(shù),m/d

    C——汲取質(zhì)濃度,mol/L

    CDb——汲取液側(cè)主體溶液中汲取質(zhì)濃度,mol/L

    CDm——汲取液側(cè)膜表面汲取質(zhì)濃度,mol/L

    CFb——原料液側(cè)主體溶液中汲取質(zhì)濃度,mol/L

    CFm——原料液側(cè)膜表面汲取質(zhì)濃度,mol/L

    Ci——活性層與支撐層界面處汲取質(zhì)濃度,mol/L

    D——汲取質(zhì)在主體溶液中的擴(kuò)散系數(shù),m2/s

    Ds——汲取質(zhì)在支撐層內(nèi)的擴(kuò)散系數(shù),m2/s

    dh——水力半徑,m

    Js——汲取質(zhì)反滲通量,mol/(m2·h)

    Jsw——汲取質(zhì)反滲通量與水通量的比率,mol/L

    Jw——水通量,L/( m2·h)

    K——支撐層內(nèi)汲取質(zhì)的阻力系數(shù)

    kD——汲取液側(cè)汲取質(zhì)的傳質(zhì)系數(shù),m/s

    kF——原料液側(cè)汲取質(zhì)的傳質(zhì)系數(shù),m/s

    L——流動(dòng)通道的長(zhǎng)度,m

    n——汲取質(zhì)產(chǎn)生的溶解物種類

    R ——理想氣體常數(shù),J/(mol·K)

    Re ——雷諾數(shù)

    S ——支撐層的結(jié)構(gòu)參數(shù)

    Sc ——施密特?cái)?shù)

    Sh ——舍伍德數(shù)

    T ——絕對(duì)溫度,K

    ts——支撐層厚度,m

    δD——汲取液側(cè)外濃差極化的厚度,m

    δF——原料液側(cè)外濃差極化的厚度,m

    ε ——支撐層的孔隙率

    μ ——黏度,Pa·s

    ρ ——水的密度,kg/m3

    τ ——支撐層的彎度

    Δπ ——汲取液側(cè)的滲透壓,Pa

    [1] Bharadwaj R,Singh D,Mahapatra A. Seawater desalination technologies[J]. International Journal of Nuclear Desalination,2008,3(2):151-159.

    [2] Paul D. Microelectronics water treatment system overview[J]. Ultrapure Water,2009,26(5):38-40.

    [3] Melo M,Schluter H,F(xiàn)erreira J,et al. Advanced performance evaluation of a reverse osmosis treatment for oilfield produced water aiming reuse[J]. Desalination,2010,250(3):1016-1018.

    [4] McCutcheon J R,McGinnis R L,Elimelech M. A novel ammonia-carbon dioxide forward (direct) osmosis desalination process[J]. Desalination,2005,174(1):1-11.

    [5] McCutcheon J R,McGinnis R L,Elimelech M. Desalination by ammonia-carbon dioxide forward osmosis:Influence of draw and feed solution concentrations on process performance[J]. Journal of Membrane Science,2006,278(1):114-123.

    [6] Kravath R E,Davis J A. Desalination of sea water by direct osmosis[J]. Desalination,1975,16(2):151-155.

    [7] Tan C H,Ng H Y. A novel hybrid forward osmosis-nanofiltration (FO-NF) process for seawater desalination:Draw solution selection and system configuration[J]. Desalination and Water Treatment,2010,13(1-3):356-361.

    [8] Choi J S,Kim H,Lee S,et al. Theoretical investigation of hybrid desalination system combining reverse osmosis and forward osmosis[J]. Desalination and Water Treatment,2010,15(1-3):114-120.

    [9] Achilli A,Cath T Y,Marchand E A,et al. The forward osmosis membrane bioreactor:A low fouling alternative to MBR processes[J]. Desalination,2009,239(1):10-21.

    [10] Holloway R W,Childress A E,Dennett K E,et al. Forward osmosis for concentration of anaerobic digester centrate[J]. Water Research,2007,41(17):4005-4014.

    [11] Cornelissen E R,Harmsen D,De Korte K F,et al. Membrane fouling and process performance of forward osmosis membranes on activated sludge[J]. Journal of Membrane Science,2008,319(1):158-168.

    [12] Cartinella J L,Cath T Y,F(xiàn)lynn M T,et al. Removal of natural steroid hormones from wastewater using membrane contactor processes[J]. Environmental Science & Technology,2006,40(23):7381-7386.

    [13] Petrotos K B,Quantick P,Petropakis H. A study of the direct osmotic concentration of tomato juice in tubular membrane-module configuration. I. The effect of certain basic process parameters on the process performance[J]. Journal of Membrane Science,1998,150(1):99-110.

    [14] Petrotos K B,Quantick P C,Petropakis H. Direct osmotic concentration of tomato juice in tubular membrane-module configuration. II. The effect of using clarified tomato juice on the process performance[J]. Journal of Membrane Science,1999,160(2):171-177.

    [15] Garcia-Castello E M,McCutcheon J R,Elimelech M. Performance evaluation of sucrose concentration using forward osmosis[J]. Journal of membrane science,2009,338(1):61-66.

    [16] Lee K L,Baker R W,Lonsdale H K. Membranes for power generation by pressure-retarded osmosis[J]. Journal of Membrane Science,1981,8(2):141-171.

    [17] Seppaèlaè A,Lampinen M J. Thermodynamic optimizing of pressure-retarded osmosis power generation systems[J]. Journal of Membrane Science,1999,161(1):115-138.

    [18] Achilli A,Cath T Y,Childress A E. Power generation with pressure retarded osmosis:An experimental and theoretical investigation[J]. Journal of Membrane Science,2009,343(1):42-52.

    [19] Cath T Y,Childress A E,Elimelech M. Forward osmosis:Principles,applications,and recent developments[J]. Journal of Membrane Science,2006,281(1):70-87.

    [20] Phillip W A,Yong J S,Elimelech M. Reverse draw solute permeation in forward osmosis:modeling and experiments[J]. Environmental Science & Technology,2010,44(13):5170-5176.

    [21] Hancock N T,Cath T Y. Solute coupled diffusion in osmotically driven membrane processes[J]. Environmental Science & Technology,2009,43(17):6769-6775.

    [22] Cath T Y,Childress A E,Elimelech M. Forward osmosis:principles,applications,and recent developments[J]. Journal of Membrane Science,2006,281(1):70-87.

    [23] Loeb S,Titelman L,Korngold E,et al. Effect of porous support fabric on osmosis through a Loeb-Sourirajan type asymmetric membrane[J]. Journal of Membrane Science,1997,129(2):243-249.

    [24] Tang C Y,She Q,Lay W C L,et al. Coupled effects of internal concentration polarization and fouling on flux behavior of forward osmosis membranes during humic acid filtration[J]. Journal of Membrane Science,2010,354(1):123-133.

    [25] Gray G T,McCutcheon J R,Elimelech M. Internal concentration polarization in forward osmosis:role of membrane orientation[J]. Desalination,2006,197(1):1-8.

    [26] Yong J S,Phillip W A,Elimelech M. Coupled reverse draw solute permeation and water flux in forward osmosis with neutral draw solutes[J]. Journal of Membrane Science,2012,392:9-17.

    [27] Suh C,Lee S. Modeling reverse draw solute flux in forward osmosis with external concentration polarization in both sides of the draw and feed solution[J]. Journal of Membrane Science,2013,427:365-374.

    [28] Song X,Liu Z,Sun D D. Nano gives the answer:breaking the bottleneck of internal concentration polarization with a nanofiber composite forward osmosis membrane for a high water productionrate[J]. Advanced Materials,2011,23(29):3256-3260.

    [29] Mehta G D,Loeb S. Performance of permasep B-9 and B-10 membranes in various osmotic regions and at high osmotic pressures[J]. Journal of Membrane Science,1979,4:335-349.

    [30] Zou S Z A L. Relating solution physicochemical properties to internal concentration polarization in forward osmosis[J]. Jornal of Mmbrane Science,2011,(1-2):459-467.

    [31] Gray G T,McCutcheon J R,Elimelech M. Internal concentration polarization in forward osmosis:Role of membrane orientation[J]. Desalination,2006,197(1):1-8.

    [32] McCutcheon J R,Elimelech M. Influence of concentrative and dilutive internal concentration polarization on flux behavior in forward osmosis[J]. Journal of Membrane Science,2006,284(1):237-247.

    [33] Ng H Y,Tang W,Wong W S. Performance of forward (direct) osmosis process:Membrane structure and transport phenomenon[J]. Environmental Science & Technology,2006,40(7):2408-2413.

    [34] Chanukya B S,Patil S,Rastogi N K. Influence of concentration polarization on flux behavior in forward osmosis during desalination using ammonium bicarbonate[J]. Desalination,2013,312:39-44.

    [35] Tiraferri A,Yip N Y,Straub A P,et al. A method for the simultaneous determination of transport and structural parameters of forward osmosis membranes[J]. Journal of Membrane Science,2013,444:523-538.

    [36] Gu B,Kim D Y,Kim J H,et al. Mathematical model of flat sheet membrane modules for FO process:Plate-and-frame module and spiral-wound module[J]. Journal of Membrane Science,2011,379(1):403-415.

    [37] Xiao D,Tang C Y,Zhang J,et al. Modeling salt accumulation in osmotic membrane bioreactors:Implications for FO membrane selection and system operation[J]. Journal of Membrane Science,2011,366(1):314-324.

    [38] Tang C Y,She Q,Lay W C L,et al. Modeling double-skinned FO membranes[J]. Desalination,2011,283:178-186.

    [39] 邊麗霞,方彥彥,王曉琳. 正滲透過程中水與溶質(zhì)的傳遞現(xiàn)象[J]. 化工學(xué)報(bào),2014,65(7):2813-2820.

    [40] 吳晴. 正滲透過程中驅(qū)動(dòng)溶質(zhì)反滲規(guī)律研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué),2013

    [41] Gruber M F,Johnson C J,Tang C Y,et al. Computational fluid dynamics simulations of flow and concentration polarization in forward osmosis membrane systems[J]. Journal of Membrane Science,2011,379(1):488-495.

    [42] Sagiv A,Semiat R. Finite element analysis of forward osmosis process using NaCl solutions[J]. Journal of Membrane Science,2011,379(1):86-96.

    [43] Su J,Chung T S. Sublayer structure and reflection coefficient and their effects on concentration polarization and membrane performance in FO processes[J]. Journal of Membrane Science,2011,376(1):214-224.

    [44] Lu X,Boo C,Ma J,et al. Bidirectional diffusion of ammonium and sodium cations in forward osmosis:Role of membrane active layer surface chemistry and charge[J]. Environmental Science & Technology,2014,48(24):14369-14376.

    [45] Boo C,Khalil Y F,Elimelech M. Performance evaluation of trimethylamine-carbon dioxide thermolytic draw solution for engineered osmosis[J]. Journal of Membrane Science,2015,473:302-309.

    [46] 宗同強(qiáng),趙玉宏,李娜. 正滲透膜成分與結(jié)構(gòu)分析[J]. 廣州化學(xué),2013,38(2):37-41.

    [47] 李海軍,張捍民,楊鳳林. 對(duì)稱正滲透膜制備及膜性能表征[J]. 環(huán)境科學(xué)與技術(shù),2013,36(12):10-14.

    [48] Wang K Y,Ong R C,Chung T S. Double-skinned forward osmosis membranes for reducing internal concentration polarization within the porous sublayer[J]. Industrial & Engineering Chemistry Research,2010,49(10):4824-4831.

    [49] Yang Q,Wang K Y,Chung T S. Dual-layer hollow fibers with enhanced flux as novel forward osmosis membranes for water production[J]. Environ. Sci. Technol.,2009,43:2800-2805.

    [50] Fang Y,Bian L,Wang X. Understanding membrane parameters of a forward osmosis membrane based on nonequilibrium thermodynamics[J]. Journal of Membrane Science,2013,437:72-81.

    [51] 吳秋燕,張忠國(guó),楊柳,等.正滲透汲取液的研究進(jìn)展[J]. 環(huán)境科學(xué)與技術(shù),2015,38(6):139-145.

    [52] Tan C H,Ng H Y. Modified models to predict flux behavior in forward osmosis in consideration of external and internal concentration polarizations[J]. Journal of Membrane Science,2008,324:209-219.

    [53] Mi B,Elimelech M. Chemical and physical aspects of organic fouling of forward osmosis membranes[J]. Journal of Membrane Science,2008,320:292-302.

    [54] Qin J J,Kekre K A,Oo M H,et al. Preliminary study of osmotic membrane bioreactor:Effects of draw solution on water flux and air scouring on fouling[J]. Water Science & Technology,2010,62(6):1353-1360.

    [55] Tang C Y,She Q,Lay W C L,et al. Coupled effects of internal concentration polarization and fouling on flux behavior of forward osmosis membranes during humic acid filtration[J]. Journal of Membrane Science,2010,354(1):123-133.

    [56] 張朝章. 關(guān)于唐南(Donnan)平衡的討論[J]. 安徽大學(xué)學(xué)報(bào):自然科學(xué)版,1986,(1):35-44.

    猜你喜歡
    錯(cuò)流滲透壓通量
    冬小麥田N2O通量研究
    高考生物問答復(fù)習(xí)之滲透壓
    錯(cuò)流過濾技術(shù)在中小工廠的應(yīng)用
    內(nèi)壓式超濾工藝處理山區(qū)水庫(kù)水的試驗(yàn)
    化基本概念為源頭活水
    ——2017年滲透壓相關(guān)高考真題賞析
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    錯(cuò)流旋轉(zhuǎn)填料床脫除細(xì)顆粒物研究
    心搏驟停后綜合征患者血漿滲透壓測(cè)定的臨床意義
    蜂蜜中耐高滲透壓酵母菌的分離與鑒定
    春、夏季長(zhǎng)江口及鄰近海域溶解甲烷的分布與釋放通量
    午夜精品在线福利| 亚洲va日本ⅴa欧美va伊人久久| 十八禁国产超污无遮挡网站| 久久午夜亚洲精品久久| 一区二区三区免费毛片| 亚洲成人久久性| 国产免费一级a男人的天堂| 日韩精品青青久久久久久| 男女啪啪激烈高潮av片| 欧美bdsm另类| 88av欧美| 亚洲精品国产成人久久av| 性插视频无遮挡在线免费观看| 性色avwww在线观看| 99精品久久久久人妻精品| 桃色一区二区三区在线观看| 国产久久久一区二区三区| 人妻丰满熟妇av一区二区三区| 窝窝影院91人妻| 午夜日韩欧美国产| 国产探花极品一区二区| h日本视频在线播放| 九九在线视频观看精品| av在线观看视频网站免费| 99久久精品热视频| 天堂√8在线中文| 精品免费久久久久久久清纯| 午夜福利视频1000在线观看| 午夜福利视频1000在线观看| 丝袜美腿在线中文| 最近中文字幕高清免费大全6 | 十八禁网站免费在线| 美女黄网站色视频| 高清日韩中文字幕在线| 欧美bdsm另类| 成人国产一区最新在线观看| 91午夜精品亚洲一区二区三区 | 国产精品爽爽va在线观看网站| 制服丝袜大香蕉在线| 91在线精品国自产拍蜜月| 给我免费播放毛片高清在线观看| 亚洲四区av| 欧美丝袜亚洲另类 | 精品人妻偷拍中文字幕| 亚洲欧美日韩东京热| 国产午夜精品论理片| 日本熟妇午夜| 网址你懂的国产日韩在线| 国产精品98久久久久久宅男小说| 夜夜看夜夜爽夜夜摸| 两个人的视频大全免费| 国内少妇人妻偷人精品xxx网站| 欧美又色又爽又黄视频| 国产三级在线视频| 婷婷六月久久综合丁香| 97超视频在线观看视频| 观看美女的网站| 别揉我奶头~嗯~啊~动态视频| 欧美另类亚洲清纯唯美| 午夜日韩欧美国产| 日本色播在线视频| 天天躁日日操中文字幕| 亚洲图色成人| 欧美区成人在线视频| 九九久久精品国产亚洲av麻豆| 日日撸夜夜添| 一本久久中文字幕| 久久精品国产亚洲av天美| 日韩一区二区视频免费看| 两个人的视频大全免费| 色综合站精品国产| 中出人妻视频一区二区| 亚洲经典国产精华液单| 国产男人的电影天堂91| 波野结衣二区三区在线| 一级a爱片免费观看的视频| 麻豆国产97在线/欧美| 俄罗斯特黄特色一大片| 88av欧美| 欧美激情国产日韩精品一区| 91在线观看av| 大又大粗又爽又黄少妇毛片口| 啦啦啦韩国在线观看视频| 人妻丰满熟妇av一区二区三区| 日韩,欧美,国产一区二区三区 | 老熟妇仑乱视频hdxx| 欧美日韩瑟瑟在线播放| 伦理电影大哥的女人| 女同久久另类99精品国产91| 欧美一区二区精品小视频在线| 亚洲精品一区av在线观看| 日韩中字成人| 他把我摸到了高潮在线观看| 免费电影在线观看免费观看| 美女高潮喷水抽搐中文字幕| 久久久久久久午夜电影| 69av精品久久久久久| 在线免费十八禁| 美女大奶头视频| bbb黄色大片| 国产真实伦视频高清在线观看 | 少妇人妻一区二区三区视频| 久久久国产成人精品二区| 久久午夜亚洲精品久久| 禁无遮挡网站| 午夜精品在线福利| 少妇高潮的动态图| 国产精品乱码一区二三区的特点| 小蜜桃在线观看免费完整版高清| 国语自产精品视频在线第100页| 日本五十路高清| 久久99热这里只有精品18| 国产高清有码在线观看视频| 免费看av在线观看网站| 欧美+日韩+精品| 国产精品久久久久久久电影| 日本色播在线视频| 91在线精品国自产拍蜜月| 在线免费观看不下载黄p国产 | 国产乱人伦免费视频| 成人国产麻豆网| 91久久精品国产一区二区成人| 最新在线观看一区二区三区| 色5月婷婷丁香| 91午夜精品亚洲一区二区三区 | 日韩一区二区视频免费看| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品色激情综合| 两个人视频免费观看高清| 非洲黑人性xxxx精品又粗又长| 久久人妻av系列| 成年版毛片免费区| 男人舔女人下体高潮全视频| av在线天堂中文字幕| 夜夜爽天天搞| 一区二区三区四区激情视频 | 哪里可以看免费的av片| 老熟妇仑乱视频hdxx| 久久精品夜夜夜夜夜久久蜜豆| 日本精品一区二区三区蜜桃| 男女之事视频高清在线观看| 好男人在线观看高清免费视频| 国产精品久久久久久久久免| 精品人妻一区二区三区麻豆 | 在线看三级毛片| 校园春色视频在线观看| 亚洲av中文字字幕乱码综合| 国产黄片美女视频| 国产三级中文精品| 香蕉av资源在线| 亚洲国产精品成人综合色| 无人区码免费观看不卡| 特大巨黑吊av在线直播| 欧美国产日韩亚洲一区| 少妇裸体淫交视频免费看高清| 国国产精品蜜臀av免费| 性色avwww在线观看| 久久99热6这里只有精品| 成年免费大片在线观看| 日本在线视频免费播放| 日韩一区二区视频免费看| 一级黄色大片毛片| 午夜免费成人在线视频| 欧美黑人欧美精品刺激| 九色成人免费人妻av| 在线观看美女被高潮喷水网站| 国产人妻一区二区三区在| 久久久久精品国产欧美久久久| 搡老熟女国产l中国老女人| 动漫黄色视频在线观看| 一级a爱片免费观看的视频| 露出奶头的视频| 国产亚洲精品久久久com| 午夜精品一区二区三区免费看| 日韩人妻高清精品专区| 嫩草影院入口| 亚洲精华国产精华液的使用体验 | 国产免费男女视频| 久久婷婷人人爽人人干人人爱| 亚洲av免费高清在线观看| 在线观看舔阴道视频| 久久精品夜夜夜夜夜久久蜜豆| 精品欧美国产一区二区三| 午夜免费男女啪啪视频观看 | 性插视频无遮挡在线免费观看| 啦啦啦观看免费观看视频高清| 亚洲真实伦在线观看| 免费搜索国产男女视频| 成人综合一区亚洲| 男女下面进入的视频免费午夜| 真人一进一出gif抽搐免费| 国产精品不卡视频一区二区| 高清日韩中文字幕在线| 很黄的视频免费| 国产黄片美女视频| 免费看美女性在线毛片视频| 国产精品一区二区免费欧美| 91av网一区二区| a在线观看视频网站| 午夜免费男女啪啪视频观看 | 嫁个100分男人电影在线观看| 国产伦精品一区二区三区四那| 亚洲成人中文字幕在线播放| 久99久视频精品免费| 老司机深夜福利视频在线观看| 国产精品国产三级国产av玫瑰| aaaaa片日本免费| 国产高清不卡午夜福利| 免费大片18禁| 国产欧美日韩精品一区二区| 国产真实伦视频高清在线观看 | 久久人人精品亚洲av| 一本久久中文字幕| 欧美日韩国产亚洲二区| 免费电影在线观看免费观看| 男人狂女人下面高潮的视频| 欧美+亚洲+日韩+国产| 亚洲精华国产精华精| 日韩在线高清观看一区二区三区 | 如何舔出高潮| 欧美又色又爽又黄视频| 最近在线观看免费完整版| 天天一区二区日本电影三级| 有码 亚洲区| 两人在一起打扑克的视频| 国产精品亚洲美女久久久| 日韩欧美三级三区| 日韩欧美在线二视频| 大型黄色视频在线免费观看| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| 欧美日韩乱码在线| 午夜视频国产福利| 久久99热这里只有精品18| 精品无人区乱码1区二区| 女同久久另类99精品国产91| 欧美在线一区亚洲| 中文资源天堂在线| 有码 亚洲区| 免费看日本二区| 日韩欧美 国产精品| 日韩av在线大香蕉| 熟妇人妻久久中文字幕3abv| 有码 亚洲区| 国内久久婷婷六月综合欲色啪| 嫩草影视91久久| 久久亚洲真实| 亚洲,欧美,日韩| АⅤ资源中文在线天堂| 国产69精品久久久久777片| 午夜福利在线观看免费完整高清在 | 国产 一区 欧美 日韩| 免费看光身美女| 欧美极品一区二区三区四区| 国内精品久久久久精免费| 一级黄色大片毛片| 最好的美女福利视频网| 观看美女的网站| 色精品久久人妻99蜜桃| 日本一二三区视频观看| 精品久久久久久久久久免费视频| 亚洲电影在线观看av| 亚洲国产欧美人成| 最近在线观看免费完整版| 成熟少妇高潮喷水视频| 一级黄色大片毛片| 久久人人精品亚洲av| 俄罗斯特黄特色一大片| 中文字幕免费在线视频6| 国产精品98久久久久久宅男小说| 亚洲av美国av| 亚洲人与动物交配视频| 欧美+日韩+精品| 自拍偷自拍亚洲精品老妇| 制服丝袜大香蕉在线| 免费看日本二区| 伦理电影大哥的女人| 草草在线视频免费看| 国产精品久久视频播放| 国产大屁股一区二区在线视频| 又紧又爽又黄一区二区| 国产精品,欧美在线| 日本欧美国产在线视频| 中文字幕久久专区| 日本黄色视频三级网站网址| 床上黄色一级片| 麻豆成人午夜福利视频| av专区在线播放| 91久久精品电影网| 久久人人爽人人爽人人片va| 午夜激情欧美在线| 亚洲无线观看免费| 美女大奶头视频| 天美传媒精品一区二区| 日本撒尿小便嘘嘘汇集6| 日本熟妇午夜| 亚洲一级一片aⅴ在线观看| 男人和女人高潮做爰伦理| 91狼人影院| 久久婷婷人人爽人人干人人爱| 亚洲欧美日韩高清专用| 91在线观看av| 久久这里只有精品中国| 午夜激情欧美在线| 免费人成在线观看视频色| 国产精品人妻久久久久久| 日本成人三级电影网站| 99久久久亚洲精品蜜臀av| 男女做爰动态图高潮gif福利片| 88av欧美| 我要看日韩黄色一级片| 国产日本99.免费观看| 国产v大片淫在线免费观看| 老司机深夜福利视频在线观看| 美女免费视频网站| 性插视频无遮挡在线免费观看| 又黄又爽又刺激的免费视频.| 人妻夜夜爽99麻豆av| 国产一区二区亚洲精品在线观看| 热99re8久久精品国产| 中国美白少妇内射xxxbb| 黄片wwwwww| 午夜老司机福利剧场| 一本精品99久久精品77| 成人av一区二区三区在线看| 淫秽高清视频在线观看| 欧美成人免费av一区二区三区| 变态另类成人亚洲欧美熟女| 亚洲av美国av| 午夜激情福利司机影院| 白带黄色成豆腐渣| 99久久中文字幕三级久久日本| 欧美3d第一页| 国产精品98久久久久久宅男小说| 欧美成人免费av一区二区三区| 天堂av国产一区二区熟女人妻| 国产 一区精品| 日本欧美国产在线视频| 欧美精品啪啪一区二区三区| 99热网站在线观看| 男人和女人高潮做爰伦理| 国产伦精品一区二区三区四那| 夜夜看夜夜爽夜夜摸| 九色成人免费人妻av| 日本黄色片子视频| 黄色一级大片看看| 久久久久久久久大av| 一个人观看的视频www高清免费观看| 久久国产精品人妻蜜桃| 一进一出抽搐gif免费好疼| 噜噜噜噜噜久久久久久91| 欧美极品一区二区三区四区| 在线观看免费视频日本深夜| 国产一区二区激情短视频| 99热这里只有是精品在线观看| 麻豆久久精品国产亚洲av| а√天堂www在线а√下载| 精品久久久久久久久久免费视频| 日日撸夜夜添| 日本一本二区三区精品| 国产精品人妻久久久影院| 午夜久久久久精精品| 99久久精品热视频| 国产综合懂色| 国产精品一区二区三区四区久久| 久久久久久久亚洲中文字幕| 婷婷精品国产亚洲av在线| 成年人黄色毛片网站| 一区二区三区激情视频| 国产欧美日韩一区二区精品| 久久热精品热| aaaaa片日本免费| 午夜福利在线在线| 很黄的视频免费| 日日啪夜夜撸| 欧美色欧美亚洲另类二区| 99国产极品粉嫩在线观看| 一级黄片播放器| 日韩精品青青久久久久久| 无人区码免费观看不卡| 亚洲精品在线观看二区| 亚洲一区高清亚洲精品| 五月伊人婷婷丁香| 久久草成人影院| 五月伊人婷婷丁香| 国产久久久一区二区三区| 草草在线视频免费看| 国内久久婷婷六月综合欲色啪| www.色视频.com| 久久精品国产鲁丝片午夜精品 | 亚洲欧美激情综合另类| 午夜影院日韩av| 久久久国产成人精品二区| 深夜a级毛片| 免费看av在线观看网站| 一a级毛片在线观看| 校园春色视频在线观看| 久久久久久伊人网av| 久久久久久久久久久丰满 | 一级av片app| 亚洲国产精品sss在线观看| 91精品国产九色| 少妇丰满av| 国产久久久一区二区三区| 亚洲综合色惰| 欧美一级a爱片免费观看看| 国产av不卡久久| 欧美一级a爱片免费观看看| 亚洲熟妇熟女久久| 日韩在线高清观看一区二区三区 | 少妇的逼水好多| or卡值多少钱| 国产探花极品一区二区| 国产精品久久视频播放| 色哟哟哟哟哟哟| 亚洲人与动物交配视频| 成人三级黄色视频| 国产av一区在线观看免费| 亚洲av免费在线观看| 大又大粗又爽又黄少妇毛片口| 国产精品日韩av在线免费观看| 欧美日本视频| 欧美zozozo另类| 欧美日韩综合久久久久久 | 少妇熟女aⅴ在线视频| 天天躁日日操中文字幕| 老熟妇乱子伦视频在线观看| 日本爱情动作片www.在线观看 | 国产av在哪里看| 最近在线观看免费完整版| 亚洲电影在线观看av| 能在线免费观看的黄片| 亚洲欧美日韩无卡精品| 成熟少妇高潮喷水视频| 97超视频在线观看视频| 人人妻人人看人人澡| 校园人妻丝袜中文字幕| av福利片在线观看| 国产真实乱freesex| 亚洲真实伦在线观看| 老司机福利观看| 久久精品国产清高在天天线| 国产激情偷乱视频一区二区| 亚洲av第一区精品v没综合| h日本视频在线播放| 免费观看的影片在线观看| 两个人的视频大全免费| 欧美3d第一页| 久久久久久伊人网av| 99热精品在线国产| 日韩欧美免费精品| 精品不卡国产一区二区三区| 麻豆成人午夜福利视频| 精品人妻视频免费看| 欧美国产日韩亚洲一区| 99在线人妻在线中文字幕| 男女视频在线观看网站免费| 国产精品98久久久久久宅男小说| 久9热在线精品视频| 亚洲真实伦在线观看| 校园春色视频在线观看| 国产男人的电影天堂91| 欧美日韩国产亚洲二区| 成人无遮挡网站| 国产女主播在线喷水免费视频网站 | 成年女人毛片免费观看观看9| 成人国产麻豆网| 午夜精品一区二区三区免费看| a级毛片免费高清观看在线播放| 欧美激情国产日韩精品一区| 99riav亚洲国产免费| av在线蜜桃| bbb黄色大片| 欧美+亚洲+日韩+国产| 91精品国产九色| 精品国产三级普通话版| 国产色爽女视频免费观看| 热99在线观看视频| 国产久久久一区二区三区| 国产欧美日韩精品亚洲av| 欧美高清成人免费视频www| 日韩强制内射视频| 欧美日韩瑟瑟在线播放| 男人和女人高潮做爰伦理| 国产一级毛片七仙女欲春2| 亚洲人成伊人成综合网2020| 日韩精品有码人妻一区| 九九爱精品视频在线观看| 天堂√8在线中文| 免费观看精品视频网站| 国产黄色小视频在线观看| 小说图片视频综合网站| 高清毛片免费观看视频网站| 午夜免费成人在线视频| 成人鲁丝片一二三区免费| www日本黄色视频网| 国产色爽女视频免费观看| 国产精品人妻久久久影院| 尾随美女入室| 97热精品久久久久久| 日本成人三级电影网站| 精品无人区乱码1区二区| 日本 av在线| 高清在线国产一区| 综合色av麻豆| 波野结衣二区三区在线| 18禁黄网站禁片午夜丰满| 日本五十路高清| 久久久久久久久久久丰满 | x7x7x7水蜜桃| 草草在线视频免费看| 天堂av国产一区二区熟女人妻| 男女边吃奶边做爰视频| 成年版毛片免费区| 国产人妻一区二区三区在| 婷婷精品国产亚洲av在线| 在线播放国产精品三级| 国产精品人妻久久久影院| 久久久久久久久久成人| 日韩高清综合在线| 日本-黄色视频高清免费观看| 成年免费大片在线观看| av国产免费在线观看| 精品一区二区三区av网在线观看| 欧美日韩乱码在线| 精品一区二区三区视频在线观看免费| 国产人妻一区二区三区在| 久久精品国产自在天天线| 精品国内亚洲2022精品成人| 久久久久久大精品| 欧美最新免费一区二区三区| 国产麻豆成人av免费视频| 午夜影院日韩av| 真实男女啪啪啪动态图| 18禁在线播放成人免费| 99热这里只有是精品50| 欧美黑人欧美精品刺激| 精品久久久噜噜| 男女做爰动态图高潮gif福利片| 老司机午夜福利在线观看视频| 国产爱豆传媒在线观看| 日韩亚洲欧美综合| 亚洲av二区三区四区| 国产精品久久视频播放| 欧美绝顶高潮抽搐喷水| 美女xxoo啪啪120秒动态图| 国产午夜精品论理片| 欧美性猛交黑人性爽| 免费av观看视频| 久久人人爽人人爽人人片va| 日日摸夜夜添夜夜添小说| 在线免费观看的www视频| 国内精品久久久久精免费| 国产精品一及| 免费观看精品视频网站| 久久久久久伊人网av| 天堂网av新在线| 日本a在线网址| 国产精品不卡视频一区二区| 国产一区二区三区在线臀色熟女| 欧美人与善性xxx| 欧美日韩乱码在线| 久久精品国产亚洲av天美| 欧美日韩精品成人综合77777| 内地一区二区视频在线| 欧美一区二区精品小视频在线| 久久久国产成人免费| 嫩草影视91久久| 日韩 亚洲 欧美在线| 国产免费av片在线观看野外av| 精品人妻一区二区三区麻豆 | 久久久久久久久大av| 国产午夜精品论理片| 亚洲av熟女| 午夜福利高清视频| 黄色丝袜av网址大全| 亚洲精品粉嫩美女一区| .国产精品久久| 99热这里只有是精品在线观看| 久久精品久久久久久噜噜老黄 | 99在线人妻在线中文字幕| 少妇高潮的动态图| 成人无遮挡网站| 国产真实伦视频高清在线观看 | 精品午夜福利在线看| 亚洲,欧美,日韩| av女优亚洲男人天堂| 久久这里只有精品中国| 欧美最黄视频在线播放免费| 性插视频无遮挡在线免费观看| 亚洲精品在线观看二区| 亚洲美女搞黄在线观看 | 一本一本综合久久| 国产不卡一卡二| 亚洲欧美激情综合另类| 精品午夜福利在线看| 精品福利观看| xxxwww97欧美| 一区福利在线观看| 久久久久免费精品人妻一区二区| 1024手机看黄色片| 亚洲真实伦在线观看| 免费观看精品视频网站| 一区二区三区免费毛片| 亚洲成人精品中文字幕电影| 乱码一卡2卡4卡精品| 国产伦人伦偷精品视频| 国内少妇人妻偷人精品xxx网站| 99久久成人亚洲精品观看| 男女边吃奶边做爰视频| 在线天堂最新版资源| 人妻久久中文字幕网| 国产一区二区亚洲精品在线观看| 久久久久国产精品人妻aⅴ院| 免费在线观看成人毛片|