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

    微生物堵塞過程中生物膜生長(zhǎng)特征對(duì)多孔介質(zhì)滲流特征影響

    2022-06-29 08:57:14王茜丹楊悅鎖武宇輝裴子丞
    中國(guó)環(huán)境科學(xué) 2022年6期
    關(guān)鍵詞:胞外生物膜菌體

    王茜丹,路 瑩,楊悅鎖,武宇輝,徐 喆,裴子丞

    微生物堵塞過程中生物膜生長(zhǎng)特征對(duì)多孔介質(zhì)滲流特征影響

    王茜丹,路 瑩*,楊悅鎖,武宇輝,徐 喆,裴子丞

    (吉林大學(xué)地下水資源與環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,吉林 長(zhǎng)春 130021)

    為研究回灌水刺激下生物膜在砂柱中的生長(zhǎng)規(guī)律,并分析生物膜生長(zhǎng)特征對(duì)堵塞介質(zhì)典型滲流特征的影響,采用室內(nèi)土柱實(shí)驗(yàn)方法,供給營(yíng)養(yǎng)液刺激砂柱內(nèi)生物膜生長(zhǎng),模擬入滲介質(zhì)生物堵塞的過程,監(jiān)測(cè)介質(zhì)滲透系數(shù)的變化,并開展不同回灌時(shí)長(zhǎng)下介質(zhì)水分吸持實(shí)驗(yàn)和彌散試驗(yàn),并對(duì)其內(nèi)生物膜形態(tài)進(jìn)行表征.結(jié)果顯示,在0~5h時(shí),滲透系數(shù)呈先下降后回升趨勢(shì);在回灌至18h前,滲透系數(shù)急劇衰減,水力彌散系數(shù)明顯增大,水分吸持能力的變化并不明顯;18h后滲透系數(shù)降低速率減緩,水分吸持能力明顯升高,溶質(zhì)運(yùn)移過渡到以彌散作用為主導(dǎo)機(jī)制.綜上,可將生物膜在介質(zhì)內(nèi)的生長(zhǎng)分為三個(gè)階段:菌適應(yīng)期,接種菌體并未生長(zhǎng)繁殖,單純由于菌體的流入與流出引起滲透系數(shù)的變化;菌體大量生長(zhǎng)繁殖期,大量生長(zhǎng)繁殖的菌體占據(jù)了部分介質(zhì)孔隙體積;胞外聚合物大量分泌期,逐漸形成具有透水性的生物膜,生物膜的生長(zhǎng)與回灌水提供的養(yǎng)分形成動(dòng)態(tài)平衡.

    地下水人工回灌;生物堵塞;生物膜;滲透系數(shù);彌散試驗(yàn);水分特征曲線

    為了應(yīng)對(duì)地下水的過度開采,利用人工手段將選用的入滲水引滲或注入目標(biāo)含水層中[1],以豐富地下水資源.通過這種人工手段,地下水資源得到了補(bǔ)給,由此引發(fā)的系列環(huán)境問題得到了有效的定期緩解,但是實(shí)踐工程中在回灌水滲入回灌介質(zhì)的滲流路徑上有微生物群落存在[2],這些微生物在水體環(huán)境中能依靠自身特性與介質(zhì)顆粒牢固粘附[3]并在適當(dāng)?shù)臈l件下迅速繁殖,微生物的生物體及其代謝物將粘附或積累在滲透介質(zhì)上,降低了介質(zhì)的滲透能力[5].人工回灌對(duì)地下水資源的管理、儲(chǔ)存和恢復(fù)有重大影響,但這種堵塞問題嚴(yán)重影響著工程的回灌效率、維護(hù)成本以及使用壽命[7],在一定程度上限制了地下水的有效補(bǔ)給和工程的高效實(shí)施[7-8].

    Baveye等[9]、姜桂華等[10]系統(tǒng)歸結(jié)了微生物作用堵塞孔隙的物質(zhì),認(rèn)為主要是因?yàn)槲⑸锏纳L(zhǎng)、衰亡、吸附作用引起細(xì)菌及其新陳代謝物質(zhì)在固體顆粒表面聚集,占據(jù)了介質(zhì)內(nèi)部孔隙通道導(dǎo)致堵塞.Bishop[11]、Flemming等[12]將這種吸附在固體顆粒表面的微生物菌體及其胞外聚合物(EPS)統(tǒng)稱為生物膜,并認(rèn)為大量的呈網(wǎng)狀結(jié)構(gòu)的EPS構(gòu)成其骨架,微生物填充其中.對(duì)于生物膜的兩種組分,Thullner等[13]觀察到微生物發(fā)展成為生物堵塞是以分散的菌落形式存在,而Vandevivere等[14]實(shí)驗(yàn)發(fā)現(xiàn)胞外聚合物也可以使介質(zhì)滲透性能顯著下降.這說明生物膜的兩種組分都會(huì)對(duì)介質(zhì)的滲透性能產(chǎn)生影響,然而EPS是附著于細(xì)胞壁上的大分子有機(jī)物質(zhì)[15],來源于細(xì)胞分泌、細(xì)胞溶解和吸附等[16],與菌體位置和性質(zhì)存在顯著差異.而生物膜在多孔介質(zhì)中的行為模式會(huì)受其自身性質(zhì)影響[17],因此當(dāng)生物膜在介質(zhì)孔隙內(nèi)積聚后,其在不同階段兩種組分的生長(zhǎng)變化會(huì)對(duì)介質(zhì)的透水能力產(chǎn)生不同影響.

    因此,本文采用滲流試驗(yàn)測(cè)定微生物生長(zhǎng)階段含水介質(zhì)的滲透系數(shù)變化,并通過離心機(jī)法和彌散試驗(yàn)對(duì)不同堵塞程度下的水分吸持能力和彌散能力進(jìn)行測(cè)定.通過比對(duì)這些滲流特征變化相對(duì)應(yīng)的生物膜組分變化,進(jìn)一步分析生物膜堵塞過程中生物膜生長(zhǎng)特征引起的介質(zhì)滲透能力的變化規(guī)律.在以往的研究中,多是關(guān)注滲流過程中各項(xiàng)滲透性能指標(biāo)整體的變化,而本研究著重于將各個(gè)滲透性能指標(biāo)的變化與生物膜的生長(zhǎng)水平相對(duì)應(yīng),得到各項(xiàng)滲透性能在不同的生物膜生長(zhǎng)水平影響下發(fā)生的階段性變化特征,為實(shí)際生物堵塞的預(yù)測(cè)與控制提供科學(xué)理論支持.

    1 材料與方法

    1.1 裝置與材料

    本試驗(yàn)采用一維柱裝置,包括層析柱、蠕動(dòng)泵、測(cè)壓管、供水裝置和出水裝置五個(gè)部分(如圖1所示).層析柱由厚度2mm的玻璃制成,內(nèi)徑1cm,高度7cm,實(shí)際裝柱高度5cm.注水口和流出口均設(shè)有400目的篩網(wǎng),防止砂樣流出.在進(jìn)出水口處分別設(shè)置測(cè)壓管監(jiān)測(cè)滲透系數(shù).

    本試驗(yàn)選用粒徑在0.125~0.25mm的細(xì)砂作為供試含水介質(zhì).為了避免物理和化學(xué)因素的影響,將細(xì)砂過篩后洗凈烘干,每次試驗(yàn)前用濃度為37%的濃鹽酸浸泡24h,并清洗至超純水在清洗前后電導(dǎo)率接近一致,滅菌烘干后,冷卻備用[18].選用地下水中優(yōu)勢(shì)菌屬[19]銅綠假單胞菌為接種液,磷酸鹽緩沖液稀釋至目標(biāo)濃度.回灌水的配比為葡萄糖1000mg/L、硝酸鉀289.2mg/L[20].

    圖1 試驗(yàn)裝置示意

    1.2 實(shí)驗(yàn)過程

    將預(yù)處理后的砂樣分層裝入層析柱中,側(cè)面敲擊令其自然沉淀,保證過程中水面始終沒過介質(zhì)一定高度,重復(fù)此過程裝至砂柱頂部,確保砂柱完全飽和[21].用蠕動(dòng)泵以0.25mL/min的恒定流速注入4mL菌液后以相同流速供入營(yíng)養(yǎng)液.每隔2h監(jiān)測(cè)滲流柱的滲透系數(shù),并設(shè)置回灌時(shí)長(zhǎng)分別為5h,18h,28h, 35.5h,44h,50.5h的對(duì)照組,利用掃描電子顯微鏡對(duì)不同回灌時(shí)長(zhǎng)下的生物膜生長(zhǎng)水平進(jìn)行表征,并測(cè)定不同生物膜生長(zhǎng)水平下的水分吸持能力和溶質(zhì)彌散能力.

    在試驗(yàn)運(yùn)行期間,對(duì)不同堵塞程度下的各組進(jìn)行示蹤劑彌散試驗(yàn),示蹤劑選用NaCl溶液,濃度為100mg/L.以0.25mL/min的速度將示蹤液泵入,收集層析柱的流出液,每個(gè)樣品200μL.用水質(zhì)分析儀檢測(cè)流出液中Cl-的濃度.將砂柱放入聚丙烯離心管中,用臺(tái)式高速冷凍離心機(jī)依次在500,1000,1500, 2000r/min的轉(zhuǎn)速下進(jìn)行離心,每種轉(zhuǎn)速下,分別離心5,10,15,25min并取樣.待離心后取出砂柱,用高精度電子天平稱量土樣質(zhì)量.根據(jù)每次離心前后質(zhì)量之差計(jì)算出土壤體積含水率并計(jì)算不同含水率下所對(duì)應(yīng)的土壤水吸力,即可繪制水分特征曲線[22-23].

    待每組實(shí)驗(yàn)運(yùn)行結(jié)束,取出砂柱進(jìn)水口和出水口處的樣品并進(jìn)行預(yù)處理,預(yù)處理方法為:使用濃度為2.5%的戊二醛溶液在4℃下制片固定.磷酸鹽緩沖溶液沖洗后,用乙醇溶液進(jìn)行梯度脫水處理[20].將樣品在真空狀態(tài)下冷凍干燥,標(biāo)號(hào)留用.然后使用掃描電子顯微鏡觀察附著在砂柱上的生物膜形態(tài).

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

    1.3.1 滲透系數(shù) 利用達(dá)西公式計(jì)算初始滲透系數(shù).

    式中:為出水口流量,m3/d;D為出入水口壓強(qiáng)差, m;D為柱高,m;為入滲砂柱橫截面面積,m2.

    采用相對(duì)滲透系數(shù)¢(實(shí)測(cè)滲透系數(shù)與初始滲透系數(shù)0之比)隨時(shí)間的變化規(guī)律反映砂柱滲透能力的變化.

    1.3.2 土壤水吸力 土壤水吸力和水分含量的變化曲線稱為土壤水分特征曲線,是土壤水力參數(shù)中的一個(gè)重要參數(shù).美國(guó)學(xué)者在1980年提出Van Genuchten方程(簡(jiǎn)稱VG方程)[24],因其線型與實(shí)測(cè)數(shù)據(jù)曲線相似且參數(shù)含義明確而在眾多描述方程中最為常用,其表達(dá)式為:

    式中:為土壤含水率,cm3/cm3;θ為土壤殘余含水率, cm3/cm3;θ為土壤飽和含水率,cm3/cm3;為土壤水吸力,cm;為土壤水分特征曲線形狀因子,共有θ、θ、4個(gè)參數(shù)[25].因此,本實(shí)驗(yàn)選用VG模型與實(shí)測(cè)數(shù)據(jù)進(jìn)行擬合,進(jìn)一步求出相關(guān)土壤水力參數(shù)進(jìn)行分析.

    1.3.3 彌散系數(shù) 通過進(jìn)行室內(nèi)彌散試驗(yàn),可以獲得示蹤劑溶液的濃度分布函數(shù)的觀測(cè)數(shù)據(jù),以砂柱孔隙體積PV為橫坐標(biāo),以出水NaCl溶液濃度為縱坐標(biāo),繪制示蹤液在多孔介質(zhì)中的穿透曲線[26].

    求解對(duì)應(yīng)生物膜生長(zhǎng)水平下的縱向彌散系數(shù)L和縱向彌散度L,可用Laplace變換得到該實(shí)驗(yàn)數(shù)學(xué)模型的解析解:

    式中:為示蹤劑溶液的濃度值,mg/L;0為流出液中示蹤劑的濃度值,mg/L;為試驗(yàn)柱內(nèi)液體實(shí)際流速,m/d;exp()為以e為底的指數(shù)函數(shù),其中e=2.718; erfc()為余誤差函數(shù).

    在關(guān)系曲線圖中找到相對(duì)濃度/0=0.159和/0=0.841兩個(gè)點(diǎn),得到與其對(duì)應(yīng)的橫坐標(biāo)PV0.159和PV0.841,其中,PV0.159和PV0.841分別指示蹤劑溶液相對(duì)濃度達(dá)到0.159和0.841時(shí)所對(duì)應(yīng)流經(jīng)的孔隙體積;將PV0.159和PV0.841分別代入式(4)和式(5)中,且令上式相減并同開平方變換,可以得到多孔介質(zhì)縱向彌散系數(shù)的計(jì)算公式:

    同上,多孔介質(zhì)縱向彌散度的計(jì)算公式如下:

    據(jù)此,可以求出縱向彌散系數(shù)L和縱向彌散度L.

    2 結(jié)果與討論

    2.1 生物膜形態(tài)變化特征

    生物膜的生命活動(dòng)在掃描電子顯微鏡下的微觀觀察中變得更為直觀.空白介質(zhì)表面沒有發(fā)現(xiàn)微生物存在的跡象(圖2).接種回灌的砂柱樣品中能明顯觀察到有異于未接種砂柱中桿狀的銅綠假單胞菌菌體.

    進(jìn)出水端的介質(zhì)顆粒表面裹覆的生物膜密度隨回灌時(shí)間的延長(zhǎng)而變高,微生物的生長(zhǎng)繁殖跡象愈加明顯.

    在掃描電鏡裝置4300倍數(shù)下觀察生物膜整體的生長(zhǎng)趨勢(shì),并在更高的倍數(shù)下(×10000)觀察到胞外聚合物的分泌情況(圖3).

    圖3 生物膜形態(tài)掃描電鏡

    根據(jù)圖3可將生物膜的發(fā)展分為三個(gè)時(shí)期,適應(yīng)期:在5h時(shí),菌體處于未進(jìn)行大量生長(zhǎng)繁殖的適應(yīng)期,代謝活性微弱;菌體大量繁殖期:在18h前,大量菌體生長(zhǎng)繁殖,沒有觀察到胞外聚合物跡象;胞外聚合物分泌期:在回灌到18h前開始分泌胞外聚合物,在28h狀態(tài)下,胞外聚合物能夠以絲狀形態(tài)被觀察到,亦稱之為纏枝狀.當(dāng)回灌時(shí)間發(fā)展到35.5h時(shí),胞外聚合物拉絲形態(tài)開始弱化,彼此之間形成帶有孔隙的絲網(wǎng)狀,在此期間菌體生長(zhǎng)繁殖速度減緩.菌體和胞外聚合物以膜狀形態(tài)厚厚裹覆在截止顆粒上.胞外聚合物隨著滲流柱堵塞的進(jìn)一步發(fā)展更加厚密.

    2.2 多孔介質(zhì)相對(duì)滲透系數(shù)變化特征

    回灌模擬試驗(yàn)中,滲流柱的相對(duì)滲透系數(shù)發(fā)生如圖4所示變化.在各組試驗(yàn)過程中,滲透系數(shù)大體上均呈衰減趨勢(shì),并可以根據(jù)相對(duì)滲透系數(shù)的衰減變化分為三個(gè)階段.

    在滲流試驗(yàn)初始階段(0~10h),相對(duì)滲透系數(shù)會(huì)呈現(xiàn)下降至近60%后再回升至接近初始滲透系數(shù)的趨勢(shì).對(duì)應(yīng)前面生物膜形態(tài)特征,此階段,菌體剛剛脫離適應(yīng)期,沒有明顯的生長(zhǎng)繁殖代謝的跡象,因此單純由菌體自身占據(jù)孔隙體積的變化影響多孔介質(zhì)的滲透系數(shù).接種菌體占據(jù)砂柱內(nèi)部一定的有效孔隙空間,滲透介質(zhì)的有效孔隙度隨之降低,導(dǎo)致滲透系數(shù)降低.隨著帶有營(yíng)養(yǎng)物質(zhì)的回灌水代替菌液以恒定流速泵入,未附著的菌體流出含水介質(zhì),被占據(jù)的那部分有效孔隙通道重新恢復(fù),滲透系數(shù)隨著有效孔隙率增加而回升.

    在中期階段(10~26h)入滲介質(zhì)的相對(duì)滲透系數(shù)均顯著下降,降低到20%以下.比對(duì)相應(yīng)階段的生物膜生長(zhǎng)變化發(fā)現(xiàn),這一階段中菌體與介質(zhì)表面接觸面積大而充分且營(yíng)養(yǎng)物質(zhì)充足,迅速繁殖,因此發(fā)生生物堵塞引起砂柱滲透性能急劇衰減,砂柱相對(duì)滲透系數(shù)快速下降.

    在后期階段(26~75h),相對(duì)滲透系數(shù)變化緩慢.此階段菌體開始大量分泌胞外聚合物,菌體的生長(zhǎng)繁殖速率下降,而胞外聚合物促進(jìn)了菌體粘連聚集形成具有透水性的生物膜,但幾乎沒有更進(jìn)一步的占據(jù)孔隙體積.因此這一階段相對(duì)滲透系數(shù)變化緩慢,堵塞程度在衰減至3%以后趨于穩(wěn)定.

    Engesgaard等[27]認(rèn)為生物堵塞后期介質(zhì)滲透性趨于穩(wěn)定是由于生物膜的生長(zhǎng)與回灌水提供的養(yǎng)分形成動(dòng)態(tài)平衡.本實(shí)驗(yàn)在介質(zhì)滲透能力降低了97%以后滲透系數(shù)開始穩(wěn)定,整個(gè)滲流過程中滲透性能衰減了3個(gè)數(shù)量級(jí).Cunningham等[28]通過定水頭滲流模擬試驗(yàn)發(fā)現(xiàn),多孔介質(zhì)由于生物堵塞,引起滲透性能衰減了92%~98%,與本實(shí)驗(yàn)結(jié)果相近.Taylor等[29]曾在室內(nèi)進(jìn)行定流量的生物堵塞滲流柱試驗(yàn),發(fā)現(xiàn)不同的滲流速度和滲流條件下,堵塞發(fā)展的空間尺度有所不同.堵塞可以造成多孔介質(zhì)水力傳導(dǎo)系數(shù)下降3個(gè)數(shù)量級(jí).Brough等[30]實(shí)驗(yàn)中發(fā)現(xiàn),35個(gè)系列試驗(yàn)中生物堵塞引起的的介質(zhì)滲透能力降低了28%~79%.大量的生物堵塞實(shí)驗(yàn)研究結(jié)果都符合Thullner[31]總結(jié)對(duì)比的生物堵塞對(duì)介質(zhì)的滲透性變化影響一般在降低兩到三個(gè)數(shù)量級(jí)的結(jié)論.這與本實(shí)驗(yàn)的實(shí)驗(yàn)結(jié)果相符,說明單純由于生物堵塞引起的滲透性能衰減具有一定的極限性.

    圖4 相對(duì)滲透系數(shù)變化情況

    2.3 介質(zhì)水分吸持能力變化特征

    根據(jù)離心實(shí)驗(yàn)實(shí)測(cè)所得數(shù)據(jù)繪制不同回灌時(shí)長(zhǎng)下的水分特征曲線,見圖5.在賦予砂柱57.647cmH2O的吸力時(shí),由于吸力較低,介質(zhì)內(nèi)的小孔隙排水為主,處在不同堵塞程度下的砂柱離心出來的水量差異并不明顯.在吸力提高到230.588cm H2O時(shí),砂樣含水量開始快速下降,不同砂柱離心所得水量出現(xiàn)差異,體積含水率在堵塞較嚴(yán)重的砂柱中明顯大于堵塞程度輕的砂柱.在更高轉(zhuǎn)速下離心時(shí)砂柱受到的吸力較高,主要是毛管孔隙中水分的變化,介質(zhì)對(duì)水分的吸持能力作用明顯,各砂柱的水分特征曲線差異趨于明顯.

    對(duì)比未注菌的空白砂柱會(huì)發(fā)現(xiàn),在相同的吸力下,回灌時(shí)間越長(zhǎng),砂柱的體積含水率越高,即生物膜在介質(zhì)中的生長(zhǎng)會(huì)提高介質(zhì)的吸持能力.生物膜的生長(zhǎng)特征對(duì)介質(zhì)水分吸持能力的影響可用土水特征參數(shù)刻畫,利用Origin將實(shí)測(cè)數(shù)據(jù)與VG方程進(jìn)行擬合求得參數(shù),如表1.

    圖5 不同回灌時(shí)間節(jié)點(diǎn)的水分特征曲線

    表1 VG方程參數(shù)

    分析表1中求得的土水特征參數(shù)發(fā)現(xiàn),土壤飽和含水率和形狀因子的變化特征并不顯著,而土壤殘余含水率隨著回灌時(shí)間的推移和生物膜生長(zhǎng)繁殖不斷升高,以殘余含水率的變化表征水分吸持能力的變化.在5h、18h、28h、44h、50.5h時(shí)殘余含水率分別為0.00229、0.01477、0.03522、0.04752、0.07633.與空白未注菌的砂柱殘余含水率0.00102相比,在回灌初期殘余含水率的升高幅度小,在18h以后,殘余含水率提高了10倍以上.

    殘余含水率反映了介質(zhì)中以薄膜形式或獨(dú)立水環(huán)的形式存在的水量,與介質(zhì)的結(jié)構(gòu)和成分有關(guān).水分吸持能力越強(qiáng),緊緊吸附在介質(zhì)顆粒表面的水量越多,殘余含水率越高.而殘余含水率與回灌時(shí)間存在明顯相關(guān)關(guān)系,繪制殘余含水率與回灌時(shí)間的關(guān)系如圖6,并比對(duì)生物膜生長(zhǎng)特征階段,將水分吸持能力的變化分區(qū).介質(zhì)的水分吸持能力變化可以以胞外聚合物大量分泌為時(shí)間節(jié)點(diǎn)分為兩個(gè)階段.

    在滲流試驗(yàn)的前18h,出現(xiàn)大量胞外聚合物節(jié)點(diǎn)前的階段,菌體的附著和生長(zhǎng)繁殖雖然引起滲透系數(shù)的急劇下降,但是殘余含水率升高趨勢(shì)并不明顯.在28h的節(jié)點(diǎn)已經(jīng)分泌大量胞外聚合物,相對(duì)于18h時(shí)殘余含水率明顯提高.胞外聚合物的大量分泌對(duì)于滲透系數(shù)衰減的影響微弱,但是對(duì)于介質(zhì)持水能力的升高有明顯的影響.

    圖6 殘余含水率變化

    介質(zhì)的持水能力隨附著菌體及其EPS產(chǎn)量的增多而愈強(qiáng).堵塞越嚴(yán)重,介質(zhì)所受生物膜生長(zhǎng)的影響越顯著,介質(zhì)的水分吸持能力越強(qiáng).而其中胞外聚合物的影響相較于菌體更加顯著.

    2.4 多孔介質(zhì)溶質(zhì)彌散能力變化特征

    經(jīng)過滲流試驗(yàn)初期滲透系數(shù)波動(dòng)階段后,分別在回灌時(shí)間為5h、18h、15h、24.5h的節(jié)點(diǎn)進(jìn)行彌散試驗(yàn).對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行擬合,如圖7所示.在得到各個(gè)節(jié)點(diǎn)下彌散試驗(yàn)的實(shí)測(cè)數(shù)據(jù)后,用Origin軟件進(jìn)行擬合,進(jìn)而分別求出各節(jié)點(diǎn)的縱向彌散系數(shù)及縱向彌散度.

    初始時(shí)介質(zhì)內(nèi)溶質(zhì)運(yùn)移以對(duì)流作用為主導(dǎo)機(jī)制,對(duì)流作用對(duì)溶質(zhì)運(yùn)移的影響大于彌散作用.在回灌時(shí)間0~18h的情況下,菌體的附著與繁殖占據(jù)了介質(zhì)內(nèi)的孔隙通道,引起滲透系數(shù)的衰減和水力彌散系數(shù)的升高.對(duì)流作用開始減弱.Engesgaard等[27]發(fā)現(xiàn)含水介質(zhì)滲透性將隨著介質(zhì)內(nèi)堵塞的發(fā)展呈現(xiàn)明顯的非均質(zhì)性,水動(dòng)力彌散等性能參數(shù)受到堵塞很大的影響,多孔介質(zhì)內(nèi)從均勻流為主導(dǎo)轉(zhuǎn)變?yōu)橄虻牧餍桶l(fā)生從近均勻流模式向非均勻流模式轉(zhuǎn)變.Seifert等[32]也發(fā)現(xiàn)由于進(jìn)水端優(yōu)先接受帶有營(yíng)養(yǎng)物質(zhì)的回灌水補(bǔ)給,所以進(jìn)水端附近附著于介質(zhì)顆粒的生物膜更加有利于形成與發(fā)展,此時(shí)微生物菌落占據(jù)孔隙喉道,導(dǎo)致滲透系數(shù)顯著下降,同時(shí)彌散度呈線性趨勢(shì)增大.

    在回灌至18h以后,胞外聚合物出現(xiàn)大量分泌趨勢(shì),主導(dǎo)機(jī)制轉(zhuǎn)變?yōu)閺浬⒆饔?溶質(zhì)彌散系數(shù)的變化開始穩(wěn)定.

    圖7 不同回灌時(shí)間節(jié)點(diǎn)的穿透曲線

    3 結(jié)論

    3.1 生物膜在介質(zhì)內(nèi)的生長(zhǎng)可以分為三個(gè)階段.在0~5h時(shí),介質(zhì)內(nèi)的接種菌處于接種后的適應(yīng)期,沒有生長(zhǎng)繁殖和分泌胞外聚合物的跡象.在回灌至18h前,介質(zhì)內(nèi)生物膜的生長(zhǎng)以菌體的生長(zhǎng)繁殖為主.在18~28h間分泌纏枝狀的胞外聚合物,逐漸與菌體共同發(fā)展成致密膜狀厚厚裹覆在介質(zhì)顆粒上.

    3.2 介質(zhì)的水分吸持能力在菌體大量生長(zhǎng)繁殖階段的前18h內(nèi)改變并不明顯;18h以后,菌體分泌大量胞外聚合物,介質(zhì)的水分吸持能力顯著增強(qiáng),殘余含水率提高10倍以上.

    3.3 介質(zhì)的彌散能力在菌體大量生長(zhǎng)繁殖階段的前18h內(nèi),由于菌體占據(jù)孔隙空間生長(zhǎng)繁殖,水力彌散系數(shù)明顯增大;18h后分泌大量胞外聚合物,溶質(zhì)運(yùn)移以彌散作用為主導(dǎo)機(jī)制,彌散系數(shù)逐漸穩(wěn)定.

    [1] Todd D K. Groundwater hydrology [M]. New York:Wiley, 1980.

    [2] 北京市地質(zhì)局水文地質(zhì)工程地質(zhì)大隊(duì),等.地下水人工補(bǔ)給[M].北京:地質(zhì)出版社, 1982:170-171.

    Beijing Bureau of Geology hydrogeology engineering Geology Brigade, et al.Artificial recharge of groundwater. Bei Jing: Geological Publishing House, 1982:170-171.

    [3] 宋 偉,劉文虎.透析液生物污染的預(yù)防與控制[J]. 中國(guó)血液化, 2008,(2):113-115.

    Song W, Liu W H. Prevention and control of biological contamination of dialysate [J].Chinese Journal of Blood Purification, 2008,(2):113- 115.

    [4] Schwager A, Boller M. Transport phenomena in intermittent filters [J]. Water Science and Technology, 1997,35(6):13-20.

    [5] Barrett M E, Taylor S. Retrofit of storm water treatment controls in a highway environment [C]//The Fifth International Conference of Sustainable Techniques and Strategies in Urban Water Management, 2004:243-250.

    [6] Jeong HY, Jun S, Cheon J, Park M. A review on clogging mechanisms and managements in aquifer storage and recovery (ASR) applications [J]. Geosciences Journal, 2018,22(4):667-679.

    [7] Barrett M E, Taylor S. Retrofit of storm water treatment controls in a highway environment [C]. The Fifth International Conference of Sustainable Techniques and Strategies in Urban Water Management, 2004:243-250.

    [8] 王世林,龔延風(fēng),陳麗萍,等.地下水人工回灌過程中生物堵塞研究進(jìn)展[J]. 科學(xué)技術(shù)與工程, 2019,19(34):25-31.

    Wang S L, Gong Y F, Chen L P, et al. Advances in bioclogging during Artificial groundwater recharge [J]. Science Technology and Engineering, 2019,19(34):25-31.

    [9] Baveye P, Vandevivere P, Hoyle B L, et al. Environmental impact and mechanisms of the biological clogging of saturated soils and aquifer materials [J]. Critical Reviews in Environmental Science and Technology, 1998,28(2):123.

    [10] 姜桂華,廖資生,徐凌云,等.人工微生物脫氮過程中含水層堵塞問題的實(shí)驗(yàn)研究[J]. 長(zhǎng)春科技大學(xué)學(xué)報(bào), 1998,(2):3-5.

    Jiang G H, Liao Z S, Xv L Y, Sun T. The experimental study on clogging of aquifer during the treatment of nitric pollution by synthetic microorganism [J]. Journal of Jilin University, 1998,(2):3-5.

    [11] Bishop P. Biofilm structure and kinetics [J]. Wat. Sci. Tech., 1997, 36(1):287-294.

    [12] Flemming H, Wingender J. The biofilm matrix [J]. Nature Reviews Microbiology, 2010,8(9): 623-633..

    [13] Thullner M, Zeyer J, Kinzelbach W.Influence of microbial growth on hydraulic properties of pore networks [J]. Transport in Porous Media, 2002,49(1):99-122.

    [14] Vandevivere P, Baveye P. Effect of bacterial extracellular polymers on the saturated hydraulic conductivity of sand columns [J]. Applied and environmental microbiology, 1992,58(5):1690-1698.

    [15] 賀張偉.預(yù)處理方法對(duì)污泥厭氧耦合微生物電解及厭氧消化產(chǎn)能的影響[D]. 哈爾濱工程大學(xué), 2014.

    He Z W. Effects of pretreatment methods of excess sludge on energy yields by anaerobic coupled microbial electrolysis and anaerobic digestion [D] Harbin Engineering University, 2014.

    [16] 陳園園,彭黨聰.生物膜中EPS與微生物的分布及測(cè)定[J]. 中國(guó)給水排水, 2017,33(15):25-30,35.

    Chen Y Y, Peng D C. Distribution and determination of EPS and microorganisms in biofilm [J]. China Water & Wastewater, 2017,33 (15):25-30,35.

    [17] 張文靜,秦運(yùn)琦,劉 丹,等.微生物在多孔介質(zhì)中環(huán)境行為研究 [J]. 中國(guó)環(huán)境科學(xué), 2018,38(10):3975-3984.

    Zhang W J, Qin Y Q, Liu D, et al. Transport behavior of microorganism in the porous media [J]. China Environmental Science, 2018,38(10):3975-3984.

    [18] 徐國(guó)丹.再生水回灌下含水介質(zhì)滲透性變化規(guī)律試驗(yàn)研究[D]. 西安:長(zhǎng)安大學(xué), 2018.

    Xv G D. Experimental reaserch on the variation of permeability of aqueous media in the process of reclaimed water recharge [D].Xi’an: Chang’an University, 2018.

    [19] 李 軍,張翠云,藍(lán)芙寧,等.區(qū)域地下水不同深度微生物群落結(jié)構(gòu)特征[J]. 中國(guó)環(huán)境科學(xué), 2019,39(6):2614-2623.

    Li J, Zhang C Y, Lan F N, et al. Structure characteristics of microbial community at different depths of groundwater [J]. China Environmental Science, 2019,39(6):2614-2623.

    [20] 沈 健,馮 晨,林允照,等.銅綠假單胞菌生物膜形成機(jī)制及耐藥相關(guān)性[J]. 浙江預(yù)防醫(yī)學(xué), 2015,27(1):40-43.

    Shen J, Feng C, Lin Y Z, et al. A study on the mechanism of biofilm formation of pseudomonas aeruginosa and its [J]. Preventive Medicine, 2015,27(1):40-43.

    [21] 趙婧彤,冶雪艷,杜新強(qiáng),等.不同濃度懸浮物顆粒在多孔介質(zhì)中遷移特性研究[J]. 水利水電技術(shù), 2019,50(10):25-31.

    Zhao J T, Ye X Y, Du X Q, et al. Study on the migration characteristics of suspended particles inporous media during different concentration of suspension [J]. Water Resources and Hydropower Engineering, 2019,50(10):25-31.

    [22] 李玉山.測(cè)定土壤水勢(shì)的離心機(jī)法[J]. 土壤, 1981,(4):143-146.

    Li Y S. Centrifuge method for determination of soil water potential [J]. Soils, 1981,(4):143-146.

    [23] 邵明安.不同方法測(cè)定土壤基質(zhì)勢(shì)的差別及準(zhǔn)確性的初步研究[J]. 土壤通報(bào), 1985,(5):222-225.

    Shao M A. A preliminary study on the difference and accuracy of different methods for measuring soil matrix potential [J]. Chinese Journal of Soil Science, 1985,(5):222-225.

    [24] Van Genuchten M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils [J]. Soil Sci Soc Am J, 1980,44(5):892-898.

    [25] 郭向紅,孫西歡,馬娟娟.基于混合遺傳算法估計(jì)van Genuchten方程參數(shù)[J]. 水科學(xué)進(jìn)展, 2009,20(5):677-682.

    Guo X H, Sun X H, Ma J J. Parametric estimation of the van Genuchten's equation based on hybrid genetic [J]. Advances in Water Science, 2009,20(5):677-682.

    [26] 張亞嬌,沈 超,李 韜,等.上海淺層粉性土彌散實(shí)驗(yàn)的室內(nèi)研究[C]// 中國(guó)建筑學(xué)會(huì)工程勘察分會(huì).第十五屆全國(guó)工程物探與巖土工程測(cè)試學(xué)術(shù)大會(huì)論文集. 中國(guó)建筑學(xué)會(huì)工程勘察分會(huì):中國(guó)建筑學(xué)會(huì)工程勘察分會(huì), 2017:5.

    Zhang Y J, Shen C, Li T, et al. Laboratory study on dispersion experiment of shallow silty soil in Shanghai [C]// Engineering Investigation Branch of Architectural Society of China. Proceedings of the 15th National Conference on Engineering Geophysical and Geotechnical Engineering Testing . Engineering Investigation Branch of Architectural Society of China: Engineering Investigation Branch of Architectural Society of China, 2017:5.

    [27] Engesgaard P, Seifert D, Herrera P. Bioclogging in porous media: Tracer Studies [M]. Riverbank Filtration Hydrology. Netherlands: Springer, 2006.

    [28] Cunningham A B, Characklis W G, Abedeen F, et al. Influence of biofilm accumulation on porous media hydrodynamics [J]. Environmental Science and Technology, 1991,25(7):1305–1311.

    [29] Taylor S W, Jaffé P R. Jaffé. Biofilm growth and the related changes in the physical properties of a porous medium: 1. Experimental investigation [J]. Water Resources Research , 1990,26(9):2153-2159.

    [30] Brough M J, Al-Tabbaa A, Martin R J.Active biofilm barriers for waste containment and bioremediation: laboratory assessment [J]. In Situ and On-site Bioremediation, 1997,4(4):233-238.

    [31] Thullner M. Comparison of bioclogging effects in saturated porous media within one- and two-dimensional flow systems [J]. Ecological Engineering, 2008,36(2):176-196.

    [32] Seifert D, Engesgaard P. Use of tracer tests to investigate changes in flow and transport properties due to bioclogging of porous media [J]. Journal of Contaminant Hydrology, 2007,93(1):58-71.

    A study of microbial clogging on the variation of seepage characteristics with biofilm growth in porous medium.

    WANG Qian-dan, LU Ying*, YANG Yue-suo, WU Yu-Hui, XU Zhe, PEI Zi-cheng

    (Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China)., 2022,42(6):2771~2778

    The seepage characteristics of porous media varies with the inner biofilm growth irritated by recharge water during groundwater recharge. Indoor soil column experiments were conducted to stimulate the bio clogging process of porous media by continuous injection with nutrient solutions. Growth morphology biofilm and the corresponding seepage parameters of media (including permeability coefficient, water retention and dispersion coefficient) were measured over a series of time periods during experiments. The results showed that the permeability coefficient slightly decreased and rose in the first 5hours. In 5~18 hours, the permeability coefficient decreased sharply, and the hydraulic dispersion coefficient, and greatly increased, but the change of water holding capacity was not obvious. After 18hours, the decrease rate of the permeability coefficient slowed down, whist the water holding capacity began to increase rapidly, and the solute transport in the porous media transited be recharged by dominant mechanism of dispersion. To sum up, the clogging process by the growth of biofilm in the medium can be divided into three stages: (1) First stage is the adaptation period of bacteria, in which the inoculated bacteria did not dramatically grow and reproduce, and the permeability coefficient changed simply because of the retention and migration of bacteria cells; (2) Second stage is the period of mass growth and reproduction of bacteria, in which the increasing bacteria cells occupy more and more pore volume of medium, resulting the decrease in permeability and the increase in dispersion of the porous media; (3) Last stage is the period of extracellular polymer secretion, in which a permeable biofilm is gradually formed, and the growth and shedding of biofilm forms a dynamic balance with the nutrients provided by recharge water.

    artificial groundwater recharge;bioclogging;biofilm;permeability coefficient;dispersion test;water characteristic curve

    X172

    A

    1000-6923(2022)06-2771-08

    王茜丹(1999-),女,吉林公主嶺人,吉林大學(xué)碩士研究生,主要從事地下水回灌堵塞模擬研究

    2021-11-18

    吉林省自然科學(xué)基金項(xiàng)目(20210101092JC)

    * 責(zé)任作者, 副教授, luying819@jlu.edu.cn

    猜你喜歡
    胞外生物膜菌體
    硝化顆粒污泥胞外聚合物及信號(hào)分子功能分析
    菌體蛋白精養(yǎng)花鰱高產(chǎn)技術(shù)探析
    東北酸菜發(fā)酵過程中菌體的分離與鑒定
    幽門螺桿菌生物膜的研究進(jìn)展
    生物膜胞外聚合物研究進(jìn)展
    水華期間藻類分層胞外聚合物與重金屬的相互作用機(jī)制研究
    菌體蛋白水解液應(yīng)用于谷氨酸發(fā)酵的研究
    黃芩苷對(duì)一株產(chǎn)NDM-1大腸埃希菌體內(nèi)外抗菌作用的研究
    光動(dòng)力對(duì)細(xì)菌生物膜的作用研究進(jìn)展
    NY3菌固定化及生物膜處理含油廢水的研究
    小蜜桃在线观看免费完整版高清| 免费人成视频x8x8入口观看| 韩国av在线不卡| 国产色爽女视频免费观看| 又黄又爽又免费观看的视频| 五月伊人婷婷丁香| 国产精品自产拍在线观看55亚洲| aaaaa片日本免费| av专区在线播放| 欧美不卡视频在线免费观看| 亚洲精品成人久久久久久| 舔av片在线| 在现免费观看毛片| 大又大粗又爽又黄少妇毛片口| 国产毛片a区久久久久| 人人妻人人澡欧美一区二区| 我要看日韩黄色一级片| 精品无人区乱码1区二区| 又黄又爽又刺激的免费视频.| 国语自产精品视频在线第100页| 国产三级中文精品| 久久久久久国产a免费观看| 嫁个100分男人电影在线观看| 一级av片app| 中文资源天堂在线| 亚洲精品在线观看二区| 免费观看人在逋| 桃色一区二区三区在线观看| 日韩国内少妇激情av| 精品人妻1区二区| 成人特级av手机在线观看| 久久久国产成人免费| 男女那种视频在线观看| ponron亚洲| 变态另类丝袜制服| 热99re8久久精品国产| 啦啦啦韩国在线观看视频| 亚洲中文字幕一区二区三区有码在线看| 国产蜜桃级精品一区二区三区| 国产精品久久久久久av不卡| 99热这里只有精品一区| 老熟妇乱子伦视频在线观看| 男人舔女人下体高潮全视频| 女人十人毛片免费观看3o分钟| 女人被狂操c到高潮| 国产精品亚洲美女久久久| 欧美又色又爽又黄视频| 在线观看美女被高潮喷水网站| xxxwww97欧美| 中国美白少妇内射xxxbb| 亚洲av中文av极速乱 | av天堂中文字幕网| 亚洲综合色惰| 最近在线观看免费完整版| 成人综合一区亚洲| 日韩一本色道免费dvd| 日韩欧美三级三区| 国产精华一区二区三区| 日日夜夜操网爽| 老熟妇仑乱视频hdxx| 九色成人免费人妻av| 别揉我奶头~嗯~啊~动态视频| 日日夜夜操网爽| 久久久久久九九精品二区国产| 床上黄色一级片| 简卡轻食公司| 国产午夜精品论理片| 国产精品福利在线免费观看| 亚洲国产精品成人综合色| 国产69精品久久久久777片| 亚洲午夜理论影院| 日韩中字成人| 久久久久国内视频| 99久久精品一区二区三区| ponron亚洲| a级毛片免费高清观看在线播放| 可以在线观看的亚洲视频| 波多野结衣巨乳人妻| 日韩精品青青久久久久久| 成人特级av手机在线观看| 国内精品久久久久久久电影| 国产伦在线观看视频一区| 亚洲欧美日韩无卡精品| 久久天躁狠狠躁夜夜2o2o| 级片在线观看| 麻豆久久精品国产亚洲av| 哪里可以看免费的av片| 两人在一起打扑克的视频| 亚洲性久久影院| av国产免费在线观看| av视频在线观看入口| 校园春色视频在线观看| 九九在线视频观看精品| 久久人人爽人人爽人人片va| 麻豆成人午夜福利视频| 亚洲精华国产精华液的使用体验 | 国产午夜福利久久久久久| 欧美黑人欧美精品刺激| 免费在线观看成人毛片| 91狼人影院| 日本五十路高清| 色在线成人网| 亚洲av免费在线观看| 淫秽高清视频在线观看| 久久亚洲真实| 很黄的视频免费| 国产大屁股一区二区在线视频| 韩国av一区二区三区四区| 欧美日韩黄片免| 午夜福利欧美成人| 欧美一区二区国产精品久久精品| 午夜老司机福利剧场| 久久精品国产99精品国产亚洲性色| 免费在线观看日本一区| 看十八女毛片水多多多| av天堂在线播放| 国产一区二区在线av高清观看| 久久久久国内视频| 亚洲精品乱码久久久v下载方式| 亚洲性久久影院| 大型黄色视频在线免费观看| 亚洲av熟女| 午夜久久久久精精品| 成年女人看的毛片在线观看| 精品久久国产蜜桃| 精品久久久久久成人av| 久久精品91蜜桃| 天堂√8在线中文| 在线免费观看不下载黄p国产 | 91久久精品国产一区二区成人| 亚洲不卡免费看| 免费观看在线日韩| 18禁在线播放成人免费| 精品人妻一区二区三区麻豆 | av天堂在线播放| 岛国在线免费视频观看| 成人一区二区视频在线观看| 婷婷色综合大香蕉| 日韩亚洲欧美综合| 日韩在线高清观看一区二区三区 | 亚洲人成伊人成综合网2020| 黄色日韩在线| 久久久久久大精品| 黄色丝袜av网址大全| 我的老师免费观看完整版| 亚洲精品国产成人久久av| 国产v大片淫在线免费观看| 黄色日韩在线| 亚洲内射少妇av| 国产日本99.免费观看| 伊人久久精品亚洲午夜| 长腿黑丝高跟| 女的被弄到高潮叫床怎么办 | 国产成人av教育| 亚洲,欧美,日韩| 热99在线观看视频| 国产成人av教育| 亚洲欧美日韩东京热| 九九久久精品国产亚洲av麻豆| 午夜精品在线福利| 成人国产综合亚洲| 亚洲av日韩精品久久久久久密| 又爽又黄无遮挡网站| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久| 免费看日本二区| 女同久久另类99精品国产91| av在线天堂中文字幕| 内地一区二区视频在线| 99国产极品粉嫩在线观看| 丝袜美腿在线中文| 美女被艹到高潮喷水动态| 色在线成人网| 成人性生交大片免费视频hd| 成人国产综合亚洲| 淫秽高清视频在线观看| 九色国产91popny在线| 亚洲精品在线观看二区| 成年人黄色毛片网站| 国产亚洲精品综合一区在线观看| 久久中文看片网| 欧美xxxx性猛交bbbb| 国产精品国产高清国产av| 欧美精品国产亚洲| 欧美成人免费av一区二区三区| 熟女电影av网| 国产伦人伦偷精品视频| 窝窝影院91人妻| eeuss影院久久| 国产欧美日韩一区二区精品| 91久久精品国产一区二区三区| 如何舔出高潮| 国内精品美女久久久久久| 国产伦人伦偷精品视频| 亚洲av成人精品一区久久| 婷婷精品国产亚洲av在线| 99国产精品一区二区蜜桃av| 少妇的逼水好多| 在线观看一区二区三区| xxxwww97欧美| 99久久久亚洲精品蜜臀av| 搡老妇女老女人老熟妇| 免费在线观看日本一区| 一进一出好大好爽视频| 91精品国产九色| 动漫黄色视频在线观看| 99国产精品一区二区蜜桃av| 99热这里只有是精品50| 少妇被粗大猛烈的视频| 国产日本99.免费观看| 九九在线视频观看精品| 一个人免费在线观看电影| 国产又黄又爽又无遮挡在线| 97人妻精品一区二区三区麻豆| 桃红色精品国产亚洲av| 日韩国内少妇激情av| 无遮挡黄片免费观看| 嫁个100分男人电影在线观看| 亚洲无线观看免费| 亚洲av成人av| 久久久久久久久久久丰满 | 高清毛片免费观看视频网站| 国产免费男女视频| 久久久久久久午夜电影| 最近最新中文字幕大全电影3| 精品人妻1区二区| 18禁裸乳无遮挡免费网站照片| 在线国产一区二区在线| 久久久久久伊人网av| 桃色一区二区三区在线观看| 国内少妇人妻偷人精品xxx网站| 女人十人毛片免费观看3o分钟| 免费在线观看日本一区| 男人舔女人下体高潮全视频| 麻豆精品久久久久久蜜桃| 亚洲国产高清在线一区二区三| 国产精品女同一区二区软件 | 有码 亚洲区| 九九爱精品视频在线观看| 欧美一区二区精品小视频在线| 免费搜索国产男女视频| 91在线观看av| 国产一区二区在线av高清观看| 久久精品国产亚洲av天美| 日日啪夜夜撸| 美女cb高潮喷水在线观看| 老司机深夜福利视频在线观看| 亚洲18禁久久av| 欧美xxxx性猛交bbbb| 黄色丝袜av网址大全| 亚洲最大成人中文| 国产精品无大码| 性插视频无遮挡在线免费观看| 欧美成人免费av一区二区三区| 91狼人影院| 在线观看av片永久免费下载| 少妇人妻一区二区三区视频| 18禁在线播放成人免费| a在线观看视频网站| 黄色丝袜av网址大全| 一区二区三区激情视频| 国产亚洲91精品色在线| 熟女人妻精品中文字幕| netflix在线观看网站| 欧美日韩精品成人综合77777| 夜夜爽天天搞| 国产成人a区在线观看| 在线观看av片永久免费下载| 亚洲最大成人av| 婷婷丁香在线五月| 小说图片视频综合网站| 少妇裸体淫交视频免费看高清| 亚洲三级黄色毛片| 欧美黑人欧美精品刺激| 国产综合懂色| 久久婷婷人人爽人人干人人爱| x7x7x7水蜜桃| 亚洲专区中文字幕在线| 性色avwww在线观看| 伊人久久精品亚洲午夜| 综合色av麻豆| 国产欧美日韩精品一区二区| 久久久成人免费电影| 又爽又黄无遮挡网站| 国产乱人视频| 麻豆成人午夜福利视频| av在线蜜桃| 成人无遮挡网站| 看免费成人av毛片| 日本爱情动作片www.在线观看 | 动漫黄色视频在线观看| 高清毛片免费观看视频网站| 欧美激情久久久久久爽电影| 少妇裸体淫交视频免费看高清| 成人欧美大片| 久久精品影院6| xxxwww97欧美| 亚洲中文日韩欧美视频| 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 久久精品国产亚洲av香蕉五月| 真人一进一出gif抽搐免费| 亚洲欧美日韩东京热| 国产精品日韩av在线免费观看| 人妻久久中文字幕网| 亚洲成人免费电影在线观看| or卡值多少钱| 亚洲国产精品久久男人天堂| 日韩精品青青久久久久久| 一个人观看的视频www高清免费观看| 天堂av国产一区二区熟女人妻| 美女免费视频网站| 国内精品宾馆在线| 内射极品少妇av片p| 国产精品人妻久久久影院| 乱系列少妇在线播放| 真人做人爱边吃奶动态| 亚洲精品在线观看二区| 日韩av在线大香蕉| 午夜久久久久精精品| 亚洲精品一区av在线观看| 午夜福利在线观看吧| 成人永久免费在线观看视频| 在线观看一区二区三区| 人妻少妇偷人精品九色| 黄色配什么色好看| 88av欧美| 国内精品久久久久久久电影| 少妇高潮的动态图| 久久这里只有精品中国| 日本一二三区视频观看| 亚洲成av人片在线播放无| 国产高清激情床上av| 精品免费久久久久久久清纯| 国产单亲对白刺激| ponron亚洲| or卡值多少钱| 成人鲁丝片一二三区免费| 亚洲色图av天堂| 亚洲av美国av| xxxwww97欧美| 色综合色国产| 麻豆国产97在线/欧美| 亚洲欧美激情综合另类| 久久久久国产精品人妻aⅴ院| 国产一区二区激情短视频| 啪啪无遮挡十八禁网站| 天堂av国产一区二区熟女人妻| 午夜福利欧美成人| 麻豆成人午夜福利视频| 少妇高潮的动态图| 久久精品国产亚洲av香蕉五月| 人妻少妇偷人精品九色| 看免费成人av毛片| 99热网站在线观看| 久久99热6这里只有精品| 免费无遮挡裸体视频| 国产精品国产三级国产av玫瑰| 禁无遮挡网站| 久久久久国产精品人妻aⅴ院| 热99在线观看视频| 午夜爱爱视频在线播放| 色5月婷婷丁香| 18禁黄网站禁片免费观看直播| 国产精华一区二区三区| 国产精品人妻久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 国产成人aa在线观看| 国产成人一区二区在线| 一a级毛片在线观看| 国产欧美日韩一区二区精品| 欧美绝顶高潮抽搐喷水| 亚洲自偷自拍三级| 国产av一区在线观看免费| 久久久久久九九精品二区国产| 搡老妇女老女人老熟妇| 国产单亲对白刺激| 成人国产一区最新在线观看| 欧美黑人巨大hd| 人人妻人人看人人澡| 国产伦精品一区二区三区视频9| 精品人妻熟女av久视频| 伦精品一区二区三区| 我的女老师完整版在线观看| 亚洲欧美激情综合另类| 黄色日韩在线| 我要看日韩黄色一级片| 亚洲av电影不卡..在线观看| 亚洲 国产 在线| 成年人黄色毛片网站| 三级男女做爰猛烈吃奶摸视频| 久久久久久大精品| 免费黄网站久久成人精品| 国产精品福利在线免费观看| 在现免费观看毛片| 欧美日韩国产亚洲二区| 国产国拍精品亚洲av在线观看| 久久草成人影院| 久久久久国产精品人妻aⅴ院| 国产不卡一卡二| 成人av在线播放网站| 亚洲经典国产精华液单| 成人特级av手机在线观看| 日本在线视频免费播放| av在线观看视频网站免费| 少妇裸体淫交视频免费看高清| 又黄又爽又刺激的免费视频.| 亚洲精品亚洲一区二区| 国产精品久久视频播放| 亚洲欧美日韩东京热| 97超视频在线观看视频| 亚洲久久久久久中文字幕| 婷婷精品国产亚洲av在线| 麻豆成人午夜福利视频| h日本视频在线播放| 国产精品精品国产色婷婷| 欧洲精品卡2卡3卡4卡5卡区| 一区二区三区激情视频| 中文字幕熟女人妻在线| 国产在线精品亚洲第一网站| 国产一区二区在线观看日韩| 欧美色欧美亚洲另类二区| 色哟哟哟哟哟哟| 亚洲aⅴ乱码一区二区在线播放| netflix在线观看网站| 国产男人的电影天堂91| 尤物成人国产欧美一区二区三区| 成人三级黄色视频| 在线看三级毛片| 人妻夜夜爽99麻豆av| 国产成人影院久久av| h日本视频在线播放| 69人妻影院| 亚洲无线观看免费| 久久久久久国产a免费观看| 观看免费一级毛片| 国产伦精品一区二区三区视频9| 婷婷丁香在线五月| 色在线成人网| 国产精品一区www在线观看 | 十八禁网站免费在线| 内地一区二区视频在线| av在线观看视频网站免费| 亚洲av第一区精品v没综合| 久久人妻av系列| 亚洲精品日韩av片在线观看| 18禁黄网站禁片午夜丰满| 一区二区三区激情视频| 69人妻影院| 久久久久久久久久成人| 综合色av麻豆| 国产欧美日韩精品一区二区| 国产精品国产三级国产av玫瑰| 热99re8久久精品国产| 亚洲成人免费电影在线观看| 99久久无色码亚洲精品果冻| 伊人久久精品亚洲午夜| 乱人视频在线观看| 成人精品一区二区免费| 精品午夜福利在线看| 国产综合懂色| 在线国产一区二区在线| 在线观看一区二区三区| 婷婷丁香在线五月| 久久国产乱子免费精品| 在线免费观看不下载黄p国产 | 中国美女看黄片| 联通29元200g的流量卡| 丰满的人妻完整版| 欧美高清性xxxxhd video| 精品久久久久久,| 国内精品久久久久久久电影| 我的老师免费观看完整版| 女人十人毛片免费观看3o分钟| 在线国产一区二区在线| 在线免费观看的www视频| 欧美另类亚洲清纯唯美| 亚洲av中文av极速乱 | 国产在线男女| 波多野结衣高清作品| 国产精品美女特级片免费视频播放器| 国内精品美女久久久久久| 久久九九热精品免费| 亚洲av二区三区四区| 亚洲av第一区精品v没综合| 色综合婷婷激情| 国产av不卡久久| 成人美女网站在线观看视频| 免费电影在线观看免费观看| 精品日产1卡2卡| 午夜免费男女啪啪视频观看 | 十八禁网站免费在线| 日韩欧美在线二视频| 精品久久久久久久久亚洲 | 91午夜精品亚洲一区二区三区 | 在线观看免费视频日本深夜| 欧美zozozo另类| 亚洲国产精品合色在线| 俄罗斯特黄特色一大片| 国产在线男女| 狂野欧美激情性xxxx在线观看| 成人欧美大片| 亚洲精品色激情综合| 精品一区二区三区av网在线观看| 色吧在线观看| 久久久久精品国产欧美久久久| 3wmmmm亚洲av在线观看| 麻豆精品久久久久久蜜桃| 一区二区三区激情视频| 日韩高清综合在线| 91狼人影院| 亚洲性夜色夜夜综合| 国产色婷婷99| 亚洲精品色激情综合| 精品人妻一区二区三区麻豆 | 日本与韩国留学比较| 中文字幕免费在线视频6| 久久国产乱子免费精品| 桃色一区二区三区在线观看| 美女大奶头视频| 色尼玛亚洲综合影院| 久久亚洲精品不卡| 亚洲美女搞黄在线观看 | 亚洲av美国av| 午夜久久久久精精品| 精品一区二区免费观看| 国产aⅴ精品一区二区三区波| 国产免费av片在线观看野外av| 国内久久婷婷六月综合欲色啪| 欧美bdsm另类| 久久精品国产鲁丝片午夜精品 | 国产精品永久免费网站| 精品久久久久久久久亚洲 | 久久热精品热| 日韩亚洲欧美综合| 色哟哟哟哟哟哟| 99热6这里只有精品| 少妇人妻一区二区三区视频| 别揉我奶头 嗯啊视频| 十八禁网站免费在线| 看黄色毛片网站| 亚洲av熟女| 日韩精品中文字幕看吧| 天天躁日日操中文字幕| 欧美日韩综合久久久久久 | 国产高清三级在线| 亚洲av成人av| 大又大粗又爽又黄少妇毛片口| 国产欧美日韩精品亚洲av| netflix在线观看网站| 两个人视频免费观看高清| 一个人看的www免费观看视频| 欧美黑人巨大hd| 国产高清不卡午夜福利| 久久久久九九精品影院| 免费观看在线日韩| 别揉我奶头~嗯~啊~动态视频| 亚洲真实伦在线观看| 国产精品久久电影中文字幕| 国产一区二区在线观看日韩| 亚洲男人的天堂狠狠| 我的女老师完整版在线观看| 国产av麻豆久久久久久久| 亚洲精品国产成人久久av| 日日夜夜操网爽| 国产精品精品国产色婷婷| 国产精品乱码一区二三区的特点| 欧美zozozo另类| 国产乱人伦免费视频| 久久久色成人| 大型黄色视频在线免费观看| 国产成人av教育| 无遮挡黄片免费观看| 欧美成人免费av一区二区三区| www日本黄色视频网| 亚洲av电影不卡..在线观看| 国产精品久久久久久精品电影| 又黄又爽又刺激的免费视频.| 亚洲av.av天堂| 成人美女网站在线观看视频| 欧美一区二区国产精品久久精品| 中文资源天堂在线| 日日干狠狠操夜夜爽| 一级黄色大片毛片| 免费观看的影片在线观看| 全区人妻精品视频| 午夜免费激情av| 人妻夜夜爽99麻豆av| 亚洲一区高清亚洲精品| 久久亚洲精品不卡| 真实男女啪啪啪动态图| 国产精品三级大全| 看免费成人av毛片| 精品欧美国产一区二区三| 婷婷精品国产亚洲av在线| 精品久久久久久成人av| 国产精品,欧美在线| 中文资源天堂在线| 欧美性猛交╳xxx乱大交人| 欧美不卡视频在线免费观看| 春色校园在线视频观看| 日韩人妻高清精品专区| 如何舔出高潮| 国产高清激情床上av| 中出人妻视频一区二区| 久久精品国产清高在天天线| 国产精品电影一区二区三区| 最近中文字幕高清免费大全6 | 成人三级黄色视频| 高清在线国产一区| 亚洲av日韩精品久久久久久密| 国产精品亚洲美女久久久| 久久亚洲真实| 国产精品综合久久久久久久免费| 如何舔出高潮|