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

    水動力作用對污染物在河流水沙兩相中分配的影響研究進展

    2015-06-14 11:13:28洋成浩科唐洪武李志偉
    關(guān)鍵詞:水沙泥沙流速

    肖 洋成浩科唐洪武李志偉

    (1.河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室,江蘇南京 210098;2.河海大學(xué)水利水電學(xué)院,江蘇南京 210098)

    水動力作用對污染物在河流水沙兩相中分配的影響研究進展

    肖 洋1,2,成浩科1,2,唐洪武1,2,李志偉1,2

    (1.河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室,江蘇南京 210098;2.河海大學(xué)水利水電學(xué)院,江蘇南京 210098)

    通過總結(jié)近年來有關(guān)水動力作用對河流中泥沙吸附/釋放重金屬、磷等污染物影響研究成果,分析了水流紊動強度和流速對污染物在河流水沙兩相中分配的影響規(guī)律。當(dāng)水流紊動強度較低時,泥沙主要以床沙形式存在,泥沙對污染物的吸附/釋放速率和強度較低;隨著紊動強度增大,水沙界面切應(yīng)力隨之增大并促使泥沙懸浮,污染物與懸浮泥沙顆粒接觸面積增大,使得泥沙對污染物的吸附/釋放速率和強度顯著增加。流速對污染物在水沙兩相間分配的影響較為復(fù)雜。泥沙靜止時,流速的增大會減小邊界層厚度,增大水體溶解氧含量、氧化還原電位等參數(shù),并使污染物在水體里的擴散由分子擴散轉(zhuǎn)變?yōu)橐晕蓜訑U散為主,增強污染物在水-沙界面的交換通量,對污染物在水沙兩相間的分配產(chǎn)生影響;隨著流速進一步增大,床沙逐漸起動、懸浮,此時除上述因素外,泥沙運動狀態(tài)、懸浮顆粒間碰撞強度、懸浮物的絮凝等均會對泥沙吸附/釋放污染物產(chǎn)生影響,由于影響因素較多,作用機制復(fù)雜,目前關(guān)于泥沙運動對污染物在水沙兩相間分配的影響所得結(jié)論仍存在分歧。相應(yīng)的機理研究及數(shù)學(xué)模型的建立均做了較大簡化。綜觀當(dāng)前的研究成果,水沙運動及其化學(xué)生物過程等對污染物在河流水沙兩相中的分配作用機理和耗氧有機物、有毒有機物等污染物與泥沙的作用機理是這一領(lǐng)域的未來研究重點。

    水動力作用;水流紊動強度;流速;泥沙運動;污染物遷移轉(zhuǎn)化

    長期以來,很多學(xué)者[1-5]通過室內(nèi)靜態(tài)試驗研究了泥沙的粒徑、理化特性、濃度以及水體污染物濃度、溫度、pH、氧化還原電位等因素對泥沙吸附/釋放污染物的影響,建立了上述各因素與泥沙吸附/釋放污染物特性之間的關(guān)系。然而,室內(nèi)靜態(tài)試驗與天然河流特征有較大差別:室內(nèi)靜態(tài)試驗一般是在較小的反應(yīng)器中,通過恒溫震蕩的方法研究泥沙對污染物的吸附/釋放作用,震蕩過程中,泥沙可與水體中污染物充分接觸,對污染物的吸附/釋放速率、總量遠大于天然河流[6];天然河流中,泥沙與污染物之間不是充分接觸,在不同的水動力作用下,泥沙運動形式的改變,會影響污染物在水沙兩相間的分配特征[5,7]。在不同水動力條件下,泥沙的運動形式包括靜止、起動、懸浮、沉降等,這個過程中,泥沙對污染物在河流水沙兩相中分配的影響主要包括3個部分(圖1):(a)污染物在床沙-上覆水界面的擴散、床沙對污染物的吸附和釋放過程;(b)懸沙沉降、床沙起動過程中對污染物吸附/釋放過程和床沙中孔隙水向上覆水的擴散過程;(c)懸沙對污染物的吸附/釋放過程。

    圖1 河流污染物在水沙兩相遷移轉(zhuǎn)化過程示意圖[8]Fig.1 Transport and transformation processes of pollutants in water and sediment[8]

    本文在總結(jié)水流紊動、流速等因素對污染物在水沙兩相間分配的影響試驗研究進展基礎(chǔ)上,討論了相應(yīng)機理和數(shù)學(xué)模型,并對將來的研究方向進行展望,以期揭示水動力對污染物在河流中遷移轉(zhuǎn)化規(guī)律的影響。目前研究中的污染物多選用營養(yǎng)鹽(磷)和重金屬,而對一般有機污染物和有毒有機污染物等方面的研究較少,故本文分析中側(cè)重于磷及重金屬污染物。

    1 水動力作用對污染物在水沙兩相間分配的影響試驗研究

    天然河流中,水流的流動狀態(tài)一般為紊流,水體的紊動強度和流速影響著泥沙的運動形式及其對污染物的吸附/釋放規(guī)律。

    1.1 紊動強度對污染物在水沙兩相間分配的影響

    研究者多采用諧振式紊動模擬裝置(圖2)形成均勻紊流,研究不同紊動強度下泥沙吸附/釋放污染物規(guī)律。如黃歲梁等[6,9]利用諧振式紊動模擬裝置研究了不同紊動強度下泥沙對重金屬的吸附特性,發(fā)現(xiàn)在低紊動強度下,泥沙以底泥形式存在,對污染物吸附速率很?。桓呶蓜訌姸认?,水沙界面切應(yīng)力增大,污染物與泥沙接觸面積、幾率增大,吸附速率、吸附量以及達到吸附平衡的時間均顯著增大。倪晉仁等[10]通過相同的裝置得到了相似的結(jié)論,認為紊動會破壞高含沙懸浮液的絮凝結(jié)構(gòu),使賓漢剪切應(yīng)力減少,促進粗顆粒沉降和底部泥沙濃度增加,進而促進泥沙的吸附作用。王鵬等[11]發(fā)現(xiàn)紊動強度對吸附的影響僅在快吸附階段較為明顯,在快吸附階段,紊動強度的增大會增加苯胺與沉積物的接觸面積,使沉積物對苯胺的吸附量增大;在慢吸附階段,泥沙對苯胺的吸附主要為化學(xué)吸附,水動力條件的影響并不明顯。在諧振式紊動模擬裝置內(nèi)紊動強度較大,水相污染物的垂直濃度分布較為均勻[12],泥沙可與污染物進行充分接觸,與真實河流環(huán)境存在較大差異。

    圖2 諧振式紊動模擬裝置[9](單位:cm)Fig.2 Resonance-type turbulence simulation equipment(units:cm)

    1.2 流速對污染物在水沙兩相間分配的影響

    流速對污染物在水沙兩相間分配的影響比較復(fù)雜。當(dāng)床沙未起動時,流速的增加會增強水沙界面的擾動,促進上覆水與孔隙水的交換,同時會改變水體的物理化學(xué)環(huán)境,如:水體溶解氧濃度、氧化還原電位、邊界層厚度等,對污染物在水沙兩相間的分配產(chǎn)生影響。隨著流速的增大,泥沙受到的水流剪切力作用增大,發(fā)生起動、懸浮、沉降,使污染物在水體中的遷移轉(zhuǎn)化過程更加復(fù)雜。

    1.2.1 床沙靜止時流速對污染物在水沙兩相間分配的影響

    圖3 循環(huán)直水槽示意圖[8](單位:m)Fig.3 Sketch of experimental straight flume(units:m)

    當(dāng)床沙靜止時,污染物主要通過擴散、吸附、解吸作用在河流水沙中遷移轉(zhuǎn)化。李斌等[13]研究底泥污染物釋放特性發(fā)現(xiàn),在床沙靜止不動時,床沙中污染物的釋放以界面釋放為主,釋放量較小,釋放時間長。黃廷林等[14]利用直水槽(圖3)研究了不同流速下鎘的釋放規(guī)律,在靜水中,泥沙中的重金屬釋放過程由Fick擴散控制;在動水條件下,隨著水流速度的提高(泥沙沒有運動),上覆水中的紊動擴散起主要作用,擴散強度大幅增加,使表層泥沙中重金屬釋放,也增大了表層下泥沙孔隙水中重金屬的擴散動力——孔隙水和上覆水之間重金屬的濃度梯度,表層下部泥沙中重金屬的釋放強度也相應(yīng)提高。House等[15-16]對不同流速下床沙對磷的吸附過程進行了研究,認為水槽中流速的增大會使水沙界面邊界層變薄,增大床沙對磷的吸附量。Barlow等[17]對黏性床沙吸附磷的研究中發(fā)現(xiàn)流速及水深與磷的吸附量之間并無顯著關(guān)系,認為在實驗室水槽規(guī)模下水深及水流流速的改變都不足以引起泥沙對磷吸附量的改變,但文中采用連續(xù)試驗的方法,在后續(xù)試驗中,隨著泥沙表面吸附磷的增多,泥沙顆粒的吸附能力降低,試驗結(jié)果同時受到流速和泥沙表面含磷飽和度的影響,不能簡單歸結(jié)為流速對泥沙吸附磷無顯著影響。與室內(nèi)的水槽試驗相比,實際河流并不是一個循環(huán)系統(tǒng),流速對泥沙吸附/釋放污染物的影響還表現(xiàn)在影響了上覆水體與床沙的接觸時間,當(dāng)水流流速較大時會加速污染物在水體中的傳輸,減小上覆水體與床沙的接觸時間,造成床沙對磷吸附量的減?。粶p小水流流速則有利于增加水體和床沙的接觸范圍和時間,增加固著生物的生長速率,進而增加污染物的吸附及顆粒態(tài)污染物的沉降[5]。

    1.2.2 泥沙運動對污染物在水沙兩相間分配的影響

    隨著流速的增大,泥沙的運動形式發(fā)生變化,床沙再懸浮過程及懸沙的沉降過程使污染物在水體中的遷移轉(zhuǎn)化過程更加復(fù)雜。Kalnejais等[18]認為床沙再懸浮過程中污染物的釋放包含床沙中液相污染物的釋放(即孔隙水的釋放)、床沙中固相污染物的釋放和再懸浮泥沙的釋放3個過程;通過螺旋槳的旋轉(zhuǎn)使床沙懸浮,對床沙懸浮過程中污染物的釋放進行研究,發(fā)現(xiàn)在臨界切應(yīng)力以下,床沙未懸浮,重金屬的釋放程度較小,隨著切應(yīng)力的增大,大于臨界切應(yīng)力時,泥沙懸浮,上覆水體中重金屬的濃度隨之增大。Cantwell等[19-20]則認為相對于切應(yīng)力,泥沙的物理、化學(xué)特性在重金屬釋放過程中的作用更加重要。

    王沛芳等[21]利用循環(huán)直水槽研究了不同流速下重金屬在沉積物-孔隙水-上覆水之間的遷移規(guī)律,發(fā)現(xiàn)流速的增加能促進重金屬的釋放,主要原因是孔隙水中的重金屬濃度遠大于上覆水中的濃度,沉積物的再懸浮使上覆水中重金屬濃度增加。黃廷林等[14]發(fā)現(xiàn)隨著流速進一步增大,重金屬的釋放濃度大幅度提高,這部分釋放不僅包括懸移質(zhì)、推移質(zhì)泥沙中重金屬的釋放,還包括靜止狀態(tài)的底部表層及底部以下泥沙中重金屬的釋放。張坤等[22]對床沙中CODCr釋放率隨水流流速的變化規(guī)律進行研究,也得到了CODCr的釋放量隨著流速增大而增大的結(jié)論。

    黃建枝等[23]利用環(huán)形水槽(圖4)研究了床沙再懸浮過程中上覆水中溶解態(tài)的Ni、Cu、Pb的濃度,發(fā)現(xiàn)隨著流速的增加,上覆水中重金屬濃度上升,認為這是由于:(a)床沙再懸浮過程中重金屬由穩(wěn)定態(tài)向可交換態(tài)轉(zhuǎn)換;(b)隨著流速的增加,顆粒間碰撞的幾率增大,使得松散結(jié)合的重金屬更容易從懸浮顆粒上解吸下來;(c)上覆水中膠體顆粒組分增大,而這部分膠體顆粒所攜帶的重金屬含量仍被視為溶解性部分,導(dǎo)致計算的重金屬釋放量增大。

    圖4 環(huán)形水槽示意圖[24]Fig.4 Sketch of annular flume[24]

    Zhang等[24]、朱紅偉等[8,25]對床沙再懸浮過程中孔隙水的作用進行分析發(fā)現(xiàn),對于非吸附性污染物,孔隙水的釋放占主要作用,對于吸附性的磷,釋放分為2個部分:再懸浮初始階段,孔隙水中液相磷的釋放起主導(dǎo)作用,持續(xù)懸浮階段,釋放主要來自于懸浮顆粒上固相磷的釋放。朱紅偉等[25]2013年的研究中出現(xiàn)了不同的現(xiàn)象,在持續(xù)懸浮階段,總磷釋放通量出現(xiàn)波動,認為這是由再懸浮顆粒的吸附及懸浮物的絮凝沉降造成的。

    朱廣偉等[26]也發(fā)現(xiàn)在由波浪引起的泥沙懸浮過程中,水體中總磷濃度有較明顯的增高,而溶解性活性磷(SRP)的濃度在開始階段顯著增高,在后期受以下2個因素的影響:(a)波浪的掀沙作用增加了水體懸浮泥沙的數(shù)量,加劇了水體中懸浮顆粒的運動和顆粒與顆粒、顆粒與離子間的碰撞,增大了水體無機和有機膠體顆粒對溶解性鹽類的吸附作用;(b)水體氧含量隨著波浪強度的增大而增大,使得沉積物再懸浮帶入水體的還原性鐵、錳等離子轉(zhuǎn)化成具有較強離子吸附能力的氧化鐵、氧化錳,水體中的可溶性活性磷又下降到掀沙前的水平。畢春娟等[27]通過泥沙再懸浮過程中上覆水中Cd、Pb、Cr的分配系數(shù)研究,認為再懸浮過程中,Cd、Pb、Cr的遷移轉(zhuǎn)化過程包括:生物膜上吸附重金屬的解吸釋放、孔隙水中溶解態(tài)重金屬的直接再懸浮和擴散釋放、沉積物中硫化物結(jié)合態(tài)重金屬的氧化釋放、上覆水中懸浮顆粒對重金屬的吸附以及沉淀等過程。彭進平等[28]測量了不同流速下水體中總磷的變化情況,發(fā)現(xiàn)隨水體流速的增大,水體中總磷的質(zhì)量濃度變化分為下降期、上升期和突增期3個階段,出現(xiàn)這種情況的原因是,低流速下懸浮顆粒的絮凝沉降造成水體中總磷減小,隨著流速的增大,沉積物顆粒少量懸浮,總磷濃度隨之升高,隨著流速繼續(xù)增大,沉積物普遍運動,總磷濃度出現(xiàn)突增,在此過程中,可溶性活性磷的濃度隨著流速的增大持續(xù)增加,說明沉積物中的磷持續(xù)向上覆水體釋放。Wan等[29]通過環(huán)形水槽研究了在水流突然擴張情況下磷的釋放行為,流速對床沙中磷向上覆水體釋放的影響過程中存在一個7~15cm/s的臨界速度,在臨界速度以下,隨著流速的增大,上覆水中溶解態(tài)總磷濃度增大,而在臨界速度以上,由于擴張段的存在,流態(tài)更為復(fù)雜,流速不是磷釋放的決定性因素,不規(guī)則的水流條件會減弱甚至阻礙磷的釋放過程。

    現(xiàn)場試驗是獲得污染物分布最真實、最有效的途徑,同時現(xiàn)場試驗中影響因素更多,泥沙、水體、生物等因素都會對試驗結(jié)果產(chǎn)生影響。Chen等[30]等對臺灣東北部Daiyujay河洪水過程中河流水體中總磷濃度進行取樣測量,發(fā)現(xiàn)當(dāng)洪峰來臨時,水體中總磷濃度突增,洪峰過后,總磷濃度先是快速降低,然后緩慢降至背景濃度。Zhang等[31]將采集泥沙放到聚氯乙烯托盤內(nèi)置于香港Crooked Island的海岸,以模擬海洋環(huán)境下,沉積物中非離子性疏水有機化合物的釋放過程,結(jié)果發(fā)現(xiàn),污染物在現(xiàn)場試驗中的釋放速率是室內(nèi)模擬試驗結(jié)果的10~100倍,作者認為這是由于室內(nèi)試驗中水的交換和沖洗效果對污染物的釋放影響低于現(xiàn)場試驗。Berg等[32]和Groote等[33]研究了清淤過程中泥沙的再懸浮情況對重金屬在水沙兩相間分配的影響,發(fā)現(xiàn)在泥沙再懸浮過程中,懸浮顆粒中重金屬的含量顯著增加,而溶解態(tài)的重金屬濃度并無顯著增加,認為在泥沙再懸浮過程中,水體氧化還原條件發(fā)生改變,泥沙顆粒表面的金屬硫化物氧化形成鐵錳氧化物/氫氧化物,造成重金屬污染物重新被吸附在泥沙顆粒上,同時藻類等對污染物的吸附也使得污染物向液相轉(zhuǎn)化的增量受限。

    2 污染物在河流水沙兩相中分配模型研究

    動水條件下,污染物與泥沙的相互作用除受泥沙顆粒特性、水環(huán)境物理化學(xué)條件等因素影響之外,還受水動力條件、泥沙運動的影響。由于影響因素眾多,作用復(fù)雜,同時由于觀測手段的限制,定量獲得各種因素的貢獻量的難度較大,因此對其過程的數(shù)學(xué)描述做了較大的簡化。

    2.1 分配系數(shù)法

    分配系數(shù)法是水質(zhì)模型中常用的描述污染物在水沙兩相中分配的方法,反映了達到吸附平衡時,污染物在固液兩相中的分配比。水質(zhì)模型一般分為整體模型和分相模型。整體模型假定推移質(zhì)和底泥對重金屬的吸附/解吸符合平衡狀態(tài),如美國環(huán)境保護署開發(fā)的水質(zhì)模型WASP[34]中,假定重金屬污染物的分配系數(shù)隨時間和空間的變化而變化,但具體如何確定分配系數(shù)的大小,WASP模型沒有很好的建議。分相模型將污染物的遷移轉(zhuǎn)化分為水相、懸移相和底泥相,如高興齋等[35]、禹雪中等[36]、余常昭[37]在研究河流中重金屬遷移轉(zhuǎn)化規(guī)律的過程中,引用平衡關(guān)系Cw=KbSc表達污染物在懸浮相和水相之間的分配關(guān)系,其中,Cw為液相污染物濃度,Sc為固相污染物濃度,Kb為分配系數(shù),在底泥傳質(zhì)中也采用了相同的方法。然而,在不同的水流條件下,泥沙對污染物的吸附/釋放速率存在差異。在泥沙的起動、懸浮、沉降、再懸浮等過程中,會加劇床沙中污染物與水體的交換強度,懸浮作用(懸浮+擴散)造成的上覆水氨氮濃度增加可以達到單純由擴散產(chǎn)生的濃度的數(shù)十倍[38],同時,在底泥的再懸浮過程中,底泥孔隙水中污染物的釋放對水體污染的貢獻也不可忽視。這些因素在上述模型中并沒有考慮,故分配系數(shù)法不能完整地描述泥沙在污染物遷移轉(zhuǎn)化過程中的作用。

    2.2 吸附/釋放動力學(xué)模型

    吸附/釋放動力學(xué)模型最初用于描述靜態(tài)試驗中泥沙對污染物的吸附/釋放過程。黃歲梁等[39]、何用等[40]、褚君達等[41]認為河流中泥沙對污染物的吸附/釋放行為與室內(nèi)靜態(tài)試驗并無顯著區(qū)別,在他們改進的水質(zhì)模型中使用了Langmuir動力學(xué)方程、Freundlich動力學(xué)方程參數(shù)對水質(zhì)模型進行優(yōu)化。

    Elovich方程是一個半經(jīng)驗公式,最初應(yīng)用于土壤及礦物吸附反應(yīng)的動力學(xué)研究,它認為吸附速率隨著吸附量的增加呈指數(shù)減小,在靜態(tài)試驗中獲得了較好的擬合效果。House[15-16]、Barlow等[17]利用Elovich方程對磷在水-床沙中的分配進行分析發(fā)現(xiàn),Elovich方程對試驗結(jié)果的擬合效果比邊界層方程和擴散方程好。然而這些動力學(xué)方程的參數(shù)通常是通過對實測數(shù)據(jù)進行擬合獲得,水動力參數(shù)對吸附/釋放的影響無法獲得。

    也有學(xué)者研究了紊動特性對泥沙吸附/釋放污染物的影響,如周孝德等[42]認為水流紊動強度對重金屬釋放動力學(xué)過程的影響是通過對某一固定時刻重金屬釋放強度的影響來體現(xiàn)的,而水流紊動強度并不隨時間而變化,僅取決于水力條件;他通過室內(nèi)水槽試驗,利用量綱分析法建立了重金屬釋放強度Ct與水流紊動強度之間的相關(guān)關(guān)系Ct=KRe1.3,其中Ct為t時刻的釋放濃度,K為釋放動力學(xué)參數(shù),由靜態(tài)的解吸動力學(xué)試驗得到,Re為雷諾數(shù)。

    2.3 邊界層模型

    停滯邊界層理論最初由Berner[43]提出,他認為在水沙界面存在一個停滯邊界層。在停滯邊界層內(nèi),磷通過擴散作用在水-沙之間遷移,邊界層的厚度與水流流速呈線性負相關(guān)關(guān)系,水流流速越大,邊界層的厚度越小,對磷擴散的阻礙越小,床沙對磷的吸附量越大。該模型[16]的形式如下:

    式中:C(t)——t時刻上覆水中污染物濃度;C0——污染物初始濃度;α——污染物平衡濃度;D0——擴散常數(shù);d——平均水深;z——停滯邊界層厚度,可通過水槽試驗來確定。

    在泥沙未起動情況下,House、Barlow利用式(1)得到了較好的擬合效果,但在有懸沙存在時,懸沙的沉降和再懸浮過程將對試驗結(jié)果產(chǎn)生影響。

    2.4 拋物線擴散模型及擴散方程

    拋物線擴散模型及擴散方程都可用來描述污染物在水沙界面的交換通量,這2個理論均認為造成污染物在上覆水-床沙間的遷移過程的驅(qū)動力為上覆水-床沙中污染物的濃度差,遵循擴散理論,其形式為

    拋物線擴散方程由Cooke[44]提出,形式如下:

    式中:β、R、A、k1——經(jīng)驗系數(shù)。

    House等[16]通過水槽試驗發(fā)現(xiàn),水流的流速越大,拋物線速率常數(shù)越大,床沙對磷的吸附量越大。張坤等[22]利用擴散方程得到了動態(tài)水流條件下底泥污染物中CODcr的釋放強度與時間的擬合關(guān)系:

    拋物線擴散模型及擴散方程主要通過上覆水體中污染物濃度變化,描述污染物在水沙界面的通量變化過程,懸沙對污染物的吸附/釋放作用,表層床沙的吸附/釋放作用對污染物在上覆水-床沙中分配的影響在這2個方程中均通過經(jīng)驗系數(shù)表示。

    2.5 三重區(qū)模型

    含氧量的大小對泥沙吸附/解吸污染物的行為有非常重要的影響。House等[45]于2002年提出了三重區(qū)模型用來描述污染物在水-沙界面及在床沙中的遷移規(guī)律,即污染物在進入床沙時,先進入邊界層,此過程可用邊界層模型進行描述,之后為床沙層,根據(jù)含氧量的區(qū)別,分為有氧層和無氧層兩層,兩層中泥沙對污染物的吸附/釋放參數(shù)存在區(qū)別。有氧層中,使用分配系數(shù)Kb進行表征,無氧層中污染物在顆粒態(tài)和溶解態(tài)的分配使用修正的分配系數(shù)Kb/f(f為大于1的參數(shù),表明泥沙在厭氧條件下的吸附量相對于有氧條件下的衰減程度)進行描述。該模型提供了一種可預(yù)測床沙中污染物垂向分配的方法,動力學(xué)因素的影響體現(xiàn)在邊界層厚度的改變。2012年提出的磷質(zhì)量守恒模型[46]中也將床沙層分為有氧層和無氧層,污染物在有氧層和無氧層中的分配系數(shù)分別使用不同的參數(shù)進行描述,而泥沙運動、水流條件、風(fēng)速等因素對于分配系數(shù)的影響在該模型中并未考慮。

    3 結(jié) 語

    水動力作用是影響污染物在河流水沙兩相中分配的重要因素,水流紊動強度、流速等水動力因素會通過改變水體物理化學(xué)環(huán)境和泥沙運動狀態(tài)對污染物在水體中的遷移轉(zhuǎn)化過程產(chǎn)生影響。由于影響泥沙對污染物的吸附/釋放作用的因素較多,目前關(guān)于不同流速對污染物在水沙兩相間的分配影響以釋放研究為主,吸附研究相對較少,且結(jié)論仍存在分歧。同時對水動力作用影響污染物在水沙兩相間分配的機理研究不夠深入,理論體系不夠完善,常用的理論及數(shù)學(xué)模型均做了較大簡化。因此,筆者認為今后的研究可從以下幾點著手:

    a.影響污染物在河流水沙中分配的因素較多,已有的研究多局限于試驗現(xiàn)象的描述,沒有很好地建立與水動力因素、泥沙運動指標的定量關(guān)系,有待建立污染物在水沙中分配規(guī)律與水流流速、水流紊動強度、懸沙濃度、推移質(zhì)輸沙率、懸移質(zhì)輸沙率、床面形態(tài)等之間的關(guān)系,為水質(zhì)模型提供合理的假設(shè)。

    b.污染物在水沙兩相間的分配包含擴散、沉降、吸附、解吸、化學(xué)反應(yīng)、降解作用、生物活動等多種作用的復(fù)雜過程,確定各種作用在污染物遷移轉(zhuǎn)化過程中的相對貢獻,對于揭示污染物在水體中的遷移轉(zhuǎn)化過程具有十分重要的意義。

    c.目前常用的簡單的、單一的數(shù)學(xué)模型對于描述現(xiàn)象具有一定的價值,但難以真實描述污染物在水沙兩相間分配的內(nèi)在機理,數(shù)學(xué)模型還需在對污染物在水沙兩相分配現(xiàn)象和機理的認識發(fā)展的基礎(chǔ)上,著重加以改進。

    d.河流中的污染物種類眾多,目前的研究多集中在重金屬及富營養(yǎng)化元素(氮磷)等與泥沙有較強親和力的元素上,對耗氧有機物、有毒有機物等污染物的研究較少。

    [1]XIAO Yang,ZHU Xiaolian,CHENG Haoke,et al.Characteristics of phosphorus adsorption by sediment mineral matrices with different particle sizes[J].Water Science and Engineering,2013,6(3):262-271.

    [2]ZHOU Aimin,TANG Hongxiao,WANG Dongsheng.Phosphorus adsorption on natural sediments:modeling and effects of pH and sediment composition[J].Water Research,2005,39(7):1245-1254.

    [3]郭長城,王國祥,喻國華.天然泥沙對富營養(yǎng)化水體中磷的吸附特性研究[J].中國給水排水,2006,22(9):10-13.(GUO Changcheng,WANG Guoxiang,YU Guohua.Study on adsorption of phosphorus in eutrophied water body by natural sediment[J].China Water and Waste Water,2006,22(9):10-13.(in Chinese))

    [4]馮海艷,李文霞,楊忠芳,等.上覆水溶解氧水平對蘇州城市河道底泥吸附/釋放磷影響的研究[J].地學(xué)前緣,2008,15(5):227-234.(FENG Haiyan,LI Wenxia,YANG Zhongfang,et al.Influence of dissolved oxygen in overlying water on releasing/absorption of phosphorus from sediments from Suzhou river,east China[J].Earth Science Frontiers,2008,15(5):227-234.(in Chinese))

    [5]WITHERS P,JARVIE H P.Delivery and cycling of phosphorus in rivers:a review[J].Science of the Total Environment,2008,400(1):379-395.

    [6]HUANG Suiliang,NG Chiuon,GUO Qizhong.Experimental investigation of the effect of flow turbulence and sediment transport patterns on the adsorption of cadmium ions onto sediment particles[J].Journal of Environmental Sciences,2007,19(6):696-703.

    [7]HOUSE W A,DENISON F H.Phosphorus dynamics in a lowland river[J].Water Research,1998,32(6):1819-1830.

    [8]朱紅偉,張坤,鐘寶昌,等.泥沙顆粒和孔隙水在底泥再懸浮污染物釋放中的作用[J].水動力學(xué)研究與進展:A輯,2011,26(5):631-641.(ZHU Hongwei,ZHANG Kun,ZHONG Baochang,et al.Effects of particles and pore water in release of pollutants due to sediment resuspension[J].Chinese Journal of Hydrodynamics,2011,26(5):631-641.(in Chinese))

    [9]黃歲梁.諧振式紊動模擬裝置解吸重金屬試驗及數(shù)學(xué)模擬[J].環(huán)境科學(xué)研究,1995(3):7-12.(HUANG Suiliang.Desorption of sediment and its mathematical model in harmonic vibration device for simulating turbulence[J].Research of Environmental Sciences,1995(3):7-12.(in Chinese))

    [10]倪晉仁,楊軍,鄧寶山.高含沙紊動系統(tǒng)中影響泥沙和銅垂向分布的關(guān)鍵因素分析[J].環(huán)境科學(xué)學(xué)報,2001(1):1-6.(NI Jinren,YANG Jun,DENG Baoshan.Controlling factors on sediment-copper vertical distribution in hyper-concentrated turbulent system[J].Acta Scientiae Circumstantiae,2001(1):1-6.(in Chinese))

    [11]王鵬,馮燕,蔡赟杰.水動力條件對太湖底泥吸附苯胺性能的影響[J].環(huán)境科學(xué)與技術(shù),2013,36(5):117-123.(WANG Peng,F(xiàn)ENG Yan,CAI Yunjie.Aniline adsorption properties on Taihu Lake sediment under different hydrodynamic conditions[J].Environmental Science and Technology,2013,36(5):117-123.(in Chinese))

    [12]楊軍,倪晉仁,鄧寶山.高含沙紊動系統(tǒng)中泥沙銅的相互作用及其垂向濃度分布[J].環(huán)境科學(xué),2001,22(1):90-93.(YANG Jun,NI Jinren,DENG Baoshan.Sediment-copper interaction and vertical distribution in hyper-concentrated turbulent system[J].Environmental Science,2001,22(1):90-93.(in Chinese))

    [13]李彬,張坤,鐘寶昌,等.底泥污染物釋放水動力特性實驗研究[J].水動力學(xué)研究與進展:A輯,2008,23(2):126-133.(LI Bin,ZHANG Kun,ZHONG Baochang,et al.An experimental study on release of pollutants from sediment under hydrodynamic conditions[J].Chinese Journal of Hydrodynamics,2008,23(2):126-133.(in Chinese))

    [14]黃廷林,周孝德.渭河沉積物中重金屬釋放的動態(tài)實驗研究[J].水利學(xué)報,1994,25(11):52-58.(HUANG Tinglin,ZHOU Xiaode.Experimental research on heavy metal release from sediment in Weihe River[J].Journal of Hydraulic Engineering,1994,25(11):52-58.(in Chinese))

    [15]HOUSE W A,DENISON F H,ARMITAGE P D.Comparison of the uptake of inorganic phosphorus to a suspended and stream bed-sediment[J].Water Research,1995,29(3):767-779.

    [16]HOUSE W A,DENISON F H,SMITH J T,et al.An investigation of the effects of water velocity on inorganic phosphorus influxto a sediment[J].Environmental Pollution,1995,89(3):263-271.

    [17]BARLOW K,NASH D,GRAYSON R.Investigating phosphorus interactions with bed sediments in a fluvial environment using a recirculating flume and intact soil cores[J].Water Research,2004,38(14):3420-3430.

    [18]KALNEJAIS L H,MARTIN W R,BOTHNER M H.The release of dissolved nutrients and metals from coastal sediments due to resuspension[J].Marine Chemistry,2010,121(1):224-235.

    [19]CANTWELL M G,BURGESS R M,KING J W.Resuspension of contaminated field and formulated reference sediments part I:Evaluation of metal release under controlled laboratory conditions[J].Chemosphere,2008,73(11):1824-1831.

    [20]CANTWELL M G,BURGESS R M,KESTER D R.Release and phase partitioning of metals from anoxic estuarine sediments during periods of simulated resuspension[J].Environmental Science&Technology,2002,36(24):5328-5334.

    [21]王沛芳,胡燕,王超,等.動水條件下重金屬在沉積物水之間的遷移規(guī)律[J].土木建筑與環(huán)境工程,2012(3):151-158.(WANG Peifang,HU Yan,WANG Chao,et al.Analysis on mobility of heavy metals between sediment-water under different hydrodynamic conditions[J].Journal of Civil,Architectural and Environmental Engineering,2012(3):151-158.(in Chinese))

    [22]張坤,李彬,王道增.動態(tài)水流條件下河流底泥污染物(CODCr)釋放研究[J].環(huán)境科學(xué)學(xué)報,2010(5):985-989.(ZHANG Kun,Li Bin,WANG Daozeng.Contaminant(CODCr)release from river bottom sediments uder flow conditions[J].Acta Scientiae Circumstantiae,2010(5):985-989.(in Chinese))

    [23]黃建枝,葛小鵬,楊曉芳,等.利用環(huán)流槽研究再懸浮條件下涼水河底泥沉積物中重金屬的遷移與分布[J].科學(xué)通報.2012,57(21):2015-2021.(HUANG Jianzhi,GE Xiaopeng,YANG Xiaofang,et al.Remobilization of heavy metals during the resuspension of Liangshui River sediments using an annular flume[J].Chinese Science Bulletin,2012,57(21):2015-2021.(in Chinese))

    [24]ZHANG Kun,CHENG Pengda,ZHONG Baochang,et al.Total phosphorus release from bottom sediments in flowing water[J].Journal of Hydrodynamics:Ser B,2012,24(4):589-594.

    [25]朱紅偉,蔣基安,程鵬達,等.泥沙污染物起動再懸浮釋放機理[J].水科學(xué)進展,2013,24(4):537-542.(ZHU Hongwei,JIANG Ji'an,CHENG Pengda,et al.Machanism of pollutant release due to sediment re-suspension[J].Advances in Water Science,2013,24(4):537-542.(in Chinese))

    [26]朱廣偉,秦伯強,張路,等.太湖底泥懸浮中營養(yǎng)鹽釋放的波浪水槽試驗[J].湖泊科學(xué),2005,17(1):61-68.(ZHU Guangwei,QIN Boqiang,ZHANG Lu,et al.Wave effects on nutrient release of sediments from lake Taihu by flume experiments[J].Journal of Lake Sciences,2005,17(1):61-68.(in Chinese))

    [27]畢春娟,陳振樓,李猛,等.再懸浮作用下長江口近岸沉積物中Cd、Pb和Cr的遷移與釋放[J].環(huán)境科學(xué),2011(9):2512-2521.(BI Chunjuan,CHEN Zhenlou,LI Meng,et al.Transport and release of Cd,Pb and Cr from the Yangtze estuarine sediments during sediment resuspension event[J].Environmental Science,2011(9):2512-2521.(in Chinese))

    [28]彭進平,逄勇,李一平,等.水動力條件對湖泊水體磷素質(zhì)量濃度的影響[J].生態(tài)環(huán)境,2003(4):388-392.(PENG Jinping,PENG Yong,LI Yiping,et al.Effects of hydrodynamic conditions on mass concentration of phosphorus[J].Ecology and Environment,2003(4):388-392.(in Chinese))

    [29]WAN Jun,WANG Ze,LI Zhijie,et al.Critical velocity in phosphorus exchange processes across the sediment-water interface[J].Journal of Environmental Sciences,2013,25(10):1966-1971.

    [30]CHEN Y,LIU J,KUO J,et al.Estimation of phosphorus flux in rivers during flooding[J].Environmental Monitoring and Assessment,2013,185(7):5653-5672.

    [31]ZHANG Y,WU R S,HONG H S,et al.Field study on desorption rates of polynuclear aromatic hydrocarbons from contaminated marine sediment[J].Environmental Toxicology and Chemistry,2000,19(10):2431-2435.

    [32]van den BERG G A,MEIJERS G G,van der Heijdt L M,et al.Dredging-related mobilisation of trace metals:a case study in the Netherlands[J].Water Research,2001,35(8):1979-1986.

    [33]de Groote J,DUMON G,VANGHELUWE M,et al.Environmental monitoring of dredging operations in the Belgian nearshore zone[J].Terra et Aqua,1998(70):19-25.

    [34]AMBROSE R B,WOOL T A,CONNOLLY J P,et al.WASP4,a hydrodynamic and water-quality model-model theory,user's manual,and programmer's guide[R].Athens:United States Environmental Protection Agency,1988.

    [35]高興齋,周佩清,黃雙雙,等.河流重金屬遷移模型的沉降與再懸浮系數(shù)研究[J].環(huán)境科學(xué)研究,1986(1):33-41.(GAO Xingzhai,ZHOU Peiqing,HUANG Shuangshuang,et al.Coefficients of settlement and resuspension in heavy metals migration model[J].Research of Environmental Sciences,1986(1):33-41.(in Chinese))

    [36]禹雪中,楊志峰,鐘德鈺,等.河流泥沙與污染物相互作用數(shù)學(xué)模型[J].水利學(xué)報,2006,37(1):10-15.(YU Xuezhong,YANG Zhifeng,ZHONG Deyu,et al.Numerical model for interaction bet ween sediment and pollutant in river[J].Journal of Hydraulic Engineering,2006,37(1):10-15.(in Chinese))

    [37]余常昭.水環(huán)境中污染物擴散輸移原理與水質(zhì)模型[M].北京:中國環(huán)境科學(xué)出版社,1989.

    [38]REDDY K R,F(xiàn)ISHER M M,IVANOFF D.Resuspension and diffusive flux of nitrogen and phosphorus in a hypereutrophic lake[J].Journal of Environmental Quality,1996,25(2):363-371.

    [39]黃歲梁,萬兆惠,張朝陽.沖積河流重金屬污染物遷移轉(zhuǎn)化數(shù)學(xué)模型研究[J].水利學(xué)報,1995,26(1):47-56.(HUANG Suiliang,WAN Zhaohui,ZHANG Chaoyang.Study on the mathematical model of heavy metal pollutants transport-transformation in fluvial rivers[J].Journal of Hydraulic Engineering,1995,26(1):47-56.(in Chinese))

    [40]何用,李義天.重金屬遷移轉(zhuǎn)化模型研究[J].水科學(xué)進展,2004,15(5):576-583.(HE Yong,LI Yitian.Study of the model of heavy metal pollutants transport[J].Advances in Water Science,2004,15(5):576-583.(in Chinese))

    [41]褚君達,徐惠慈.河流底泥沖刷沉降對水質(zhì)影響的研究[J].水利學(xué)報,1994,25(11):42-47.(CHU Junda,XU Huici.Research on the effects of river sediment scouring sedimentation of water quality[J].Journal of Hydraulic Engineering,1994,25(11):42-47.(in Chinese))

    [42]周孝德,黃廷林.河流底流中重金屬釋放的水流紊動效應(yīng)[J].水利學(xué)報,1994,25(11):22-25.(ZHOU Xiaode,HUANG Tinglin.Effect of water turbulence on the release of heavy metals in bottom stream[J].Journal of Hydraulic Engineering,1994(11):22-25.(in Chinese)))

    [43]BERNER R A.Early diagenesis:a theoretical approach[M].Princeton:Princeton University Press,1980.

    [44]COOKE I J.A kinetic approach to the description of soil phosphate status[J].Journal of Soil Science,1966,17(1):56-64.

    [45]HOUSE W A,DENISON F H.Exchange of inorganic phosphate between river waters and bed-sediments[J].Environmental Science&Technology,2002,36(20):4295-4301.

    [46]CHOMAT C J,WESTPHAL K S.Enhanced understanding of sediment phosphorus dynamics in river systems with a simple supplemental mass balance tool[J].Journal of Environmental Engineering,2012,139(1):34-43.

    Review of influence of hydrodynamic action on distribution of pollutants in water and sediment in river

    XIAO Yang1,2,CHENG Haoke1,2,TANG Hongwu1,2,LI Zhiwei1,2
    (1.State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Hohai University,Nanjing 210098,China;2.College of Water Conservancy and Hydropower Engineering,Hohai University,Nanjing 210098,China)

    Through a review of recent investigations into the effects of hydrodynamic action on the adsorption/release of pollutants(heavy metals,phosphorus,etc.)by sediments,the influence of water turbulence intensity and velocity on the distribution of pollutants in water and sediment in rivers was analyzed.The results show that,when the turbulence intensity is low,sediments occur in a form close to bed sediment,and the rate and extent of adsorption/release of pollutants by sediments are low.As the turbulence intensity increases,the shear stress at the water-sediment interface increases,sediment suspension is promoted,and the contact area between pollutants and suspended sediment particles increases,leading to a significant increase in the rate and extent of adsorption/release of pollutants by sediments.The effect of the water velocity on the distribution of pollutants in water and sediment iscomplex.When the sediment is stationary,with the increase of the velocity,the boundary layer thickness decreases,the dissolved oxygen content in water and redox potential increase,the diffusion of pollutants in the water body changes from molecular diffusion to turbulent diffusion,and the pollutant exchange flux at the watersediment interface increases.With further increase of the velocity,the bed sediment gradually moves and is suspended.At this moment,the state of sediment movement,intensity of collision between suspended matters,and flocculation of sediments also have effects on the adsorption/release of pollutants by sediments.Because of the extensive factors and the complexity of the functional mechanism,there is a disagreement between results regarding the influence of sediment movement on the adsorption/release of pollutants by sediments reported by different researchers,and simplification is carried out through mechanism studies and establishment of mathematical models.It is pointed out that future studies should focus on the mechanism of the influence of the water and sediment movement as well as chemical and biological processes on the distribution of pollutants in water and sediment in rivers and the interaction of oxygen consumption organics and nocuous organics with sediments.

    hydrodynamic action;water turbulent intensity;velocity;sediment movement;transport and transformation of pollutant

    X522;G353.11

    A

    1000-1980(2015)05-0480-09

    10.3876/j.issn.1000-1980.2015.05.013

    20150327

    國家自然科學(xué)基金(51179055,51239003);江蘇高校優(yōu)勢學(xué)科建設(shè)工程項目(SYS1401);中央高校基本科研業(yè)務(wù)費專項(2013/B14020362)

    肖洋(1974—),男,貴州銅仁人,教授,博士,主要從事水力學(xué)及河流動力學(xué)研究。E-mail:sediment_lab@hhu.edu.cn

    猜你喜歡
    水沙泥沙流速
    “流體壓強與流速的關(guān)系”知識鞏固
    泥沙做的父親
    『流體壓強與流速的關(guān)系』知識鞏固
    山雨欲來風(fēng)滿樓之流體壓強與流速
    新疆多泥沙河流水庫泥沙處理措施
    大型水利樞紐下游水沙變異特征
    土壤團聚體對泥沙沉降速度的影響
    愛虛張聲勢的水
    山區(qū)河流上下雙丁壩回流區(qū)水沙特性淺探
    江西建材(2018年1期)2018-04-04 05:26:28
    泥沙滅火
    兒童繪本(2015年2期)2015-05-25 18:10:15
    丰满乱子伦码专区| 午夜日本视频在线| 91精品国产九色| 国产精品久久久久久av不卡| 少妇丰满av| 国产一区二区在线观看日韩| 全区人妻精品视频| 人人妻人人爽人人添夜夜欢视频 | 国国产精品蜜臀av免费| 色吧在线观看| 精品国产乱码久久久久久小说| 国产女主播在线喷水免费视频网站| 免费大片黄手机在线观看| 在线观看免费视频网站a站| 韩国高清视频一区二区三区| 女人十人毛片免费观看3o分钟| 免费不卡的大黄色大毛片视频在线观看| 汤姆久久久久久久影院中文字幕| 国产熟女欧美一区二区| 99国产精品免费福利视频| 久久精品熟女亚洲av麻豆精品| 亚洲国产成人一精品久久久| 大香蕉97超碰在线| 成人影院久久| 一级毛片 在线播放| 中文欧美无线码| 国内少妇人妻偷人精品xxx网站| .国产精品久久| 人人妻人人看人人澡| 老女人水多毛片| 欧美极品一区二区三区四区| 狂野欧美激情性bbbbbb| 日本猛色少妇xxxxx猛交久久| 一区二区三区四区激情视频| 一本一本综合久久| 日本黄色片子视频| 国产爽快片一区二区三区| 又大又黄又爽视频免费| freevideosex欧美| 久久国产乱子免费精品| 久久ye,这里只有精品| 免费观看的影片在线观看| 亚洲av不卡在线观看| 亚洲av中文字字幕乱码综合| 久久久久久久久久久免费av| 久久精品人妻少妇| 97超视频在线观看视频| 各种免费的搞黄视频| 日韩av免费高清视频| 亚州av有码| 国产成人91sexporn| 亚洲精品色激情综合| 91午夜精品亚洲一区二区三区| 简卡轻食公司| 精品久久久久久久末码| 久久99热6这里只有精品| 熟妇人妻不卡中文字幕| 少妇精品久久久久久久| 尾随美女入室| 亚洲色图av天堂| 久久99精品国语久久久| 最新中文字幕久久久久| av国产免费在线观看| 久久99蜜桃精品久久| 3wmmmm亚洲av在线观看| 婷婷色综合大香蕉| 少妇人妻久久综合中文| 国产精品秋霞免费鲁丝片| 国产精品久久久久成人av| 一本久久精品| 日本欧美国产在线视频| 国产深夜福利视频在线观看| 国产毛片在线视频| 九九在线视频观看精品| 身体一侧抽搐| 亚洲欧美一区二区三区黑人 | 少妇猛男粗大的猛烈进出视频| 欧美性感艳星| 日本欧美国产在线视频| 简卡轻食公司| 国产乱来视频区| 中文字幕人妻熟人妻熟丝袜美| 一级爰片在线观看| 久久鲁丝午夜福利片| 两个人的视频大全免费| 国产一区二区三区av在线| 国产男人的电影天堂91| 舔av片在线| videossex国产| 亚洲精品中文字幕在线视频 | 久久国内精品自在自线图片| 妹子高潮喷水视频| 国产欧美日韩一区二区三区在线 | 欧美性感艳星| 亚洲国产色片| 亚洲国产毛片av蜜桃av| 老司机影院毛片| 免费观看a级毛片全部| 男女国产视频网站| 亚洲国产精品成人久久小说| 亚洲国产毛片av蜜桃av| 欧美最新免费一区二区三区| 搡女人真爽免费视频火全软件| 久久97久久精品| 免费黄色在线免费观看| 国产在线男女| 欧美精品人与动牲交sv欧美| 18禁裸乳无遮挡动漫免费视频| 男男h啪啪无遮挡| 国产在视频线精品| av在线app专区| 一个人看的www免费观看视频| 26uuu在线亚洲综合色| 午夜免费鲁丝| 久久99热6这里只有精品| 午夜视频国产福利| 欧美激情极品国产一区二区三区 | 黄色视频在线播放观看不卡| 欧美老熟妇乱子伦牲交| 99精国产麻豆久久婷婷| 伦理电影大哥的女人| 日韩欧美 国产精品| 亚洲伊人久久精品综合| 人人妻人人添人人爽欧美一区卜 | 日韩中字成人| 91久久精品国产一区二区三区| 最黄视频免费看| 欧美一级a爱片免费观看看| 美女脱内裤让男人舔精品视频| 精品国产三级普通话版| 观看美女的网站| 99热全是精品| 丝袜喷水一区| 国产一区二区三区综合在线观看 | h日本视频在线播放| 3wmmmm亚洲av在线观看| 成人漫画全彩无遮挡| 欧美xxxx性猛交bbbb| 国产欧美日韩精品一区二区| av网站免费在线观看视频| 欧美精品国产亚洲| 欧美精品亚洲一区二区| av视频免费观看在线观看| 国产探花极品一区二区| 国产男女超爽视频在线观看| 日韩不卡一区二区三区视频在线| 卡戴珊不雅视频在线播放| 亚洲不卡免费看| 久久鲁丝午夜福利片| 欧美bdsm另类| xxx大片免费视频| 精品一区在线观看国产| 国产美女午夜福利| 99久久综合免费| 亚洲高清免费不卡视频| 亚洲av综合色区一区| 免费在线观看成人毛片| 男人和女人高潮做爰伦理| 国产高清三级在线| 18+在线观看网站| 精品熟女少妇av免费看| 国产成人精品久久久久久| 自拍偷自拍亚洲精品老妇| 久久久久久久久久成人| 免费在线观看成人毛片| 一个人看视频在线观看www免费| 国内少妇人妻偷人精品xxx网站| 午夜激情久久久久久久| 亚洲欧美清纯卡通| 蜜桃亚洲精品一区二区三区| 国产精品一二三区在线看| 成人二区视频| 精品视频人人做人人爽| 欧美高清成人免费视频www| 日韩电影二区| 身体一侧抽搐| 黄片wwwwww| 舔av片在线| 国产视频内射| 久久精品国产亚洲av涩爱| 日日啪夜夜撸| 久久久精品94久久精品| 久久精品夜色国产| 精品一区二区三卡| 夜夜看夜夜爽夜夜摸| 3wmmmm亚洲av在线观看| 26uuu在线亚洲综合色| 国精品久久久久久国模美| 欧美 日韩 精品 国产| 欧美高清成人免费视频www| 日韩av在线免费看完整版不卡| 高清av免费在线| 国产精品久久久久成人av| 在线观看人妻少妇| av.在线天堂| 少妇人妻久久综合中文| 欧美成人a在线观看| 极品教师在线视频| 亚洲,欧美,日韩| 欧美精品人与动牲交sv欧美| 有码 亚洲区| 啦啦啦视频在线资源免费观看| 少妇熟女欧美另类| 久久久精品免费免费高清| 黄色日韩在线| 免费观看a级毛片全部| 国产精品女同一区二区软件| freevideosex欧美| 国产精品av视频在线免费观看| 国产精品欧美亚洲77777| 亚洲va在线va天堂va国产| 高清毛片免费看| 高清视频免费观看一区二区| 日韩人妻高清精品专区| 黑人高潮一二区| 简卡轻食公司| 小蜜桃在线观看免费完整版高清| 国产av精品麻豆| 五月玫瑰六月丁香| 久久久久国产网址| 91精品国产九色| 精品人妻一区二区三区麻豆| 久久久久久久久久成人| 日本-黄色视频高清免费观看| 最黄视频免费看| 毛片一级片免费看久久久久| 国产乱来视频区| 在现免费观看毛片| 久久久国产一区二区| 91精品伊人久久大香线蕉| 久久精品国产亚洲网站| 精品少妇黑人巨大在线播放| 精品一区二区免费观看| 国产淫片久久久久久久久| 午夜视频国产福利| 欧美精品一区二区免费开放| 午夜免费鲁丝| 欧美极品一区二区三区四区| 国产有黄有色有爽视频| 久久精品久久久久久噜噜老黄| 少妇猛男粗大的猛烈进出视频| 国产深夜福利视频在线观看| 久久99热这里只频精品6学生| 人妻少妇偷人精品九色| 秋霞伦理黄片| 亚洲av中文av极速乱| 观看美女的网站| 国产精品久久久久久精品古装| 午夜免费男女啪啪视频观看| 男女边摸边吃奶| 精品一区在线观看国产| 亚洲精品456在线播放app| 99久久精品一区二区三区| 男女边摸边吃奶| 久久久久人妻精品一区果冻| 日韩一区二区三区影片| 小蜜桃在线观看免费完整版高清| 日本色播在线视频| 久久97久久精品| 新久久久久国产一级毛片| 久久热精品热| 欧美日韩视频高清一区二区三区二| 国产欧美日韩一区二区三区在线 | 91aial.com中文字幕在线观看| av播播在线观看一区| 国产欧美亚洲国产| 国产女主播在线喷水免费视频网站| 成人黄色视频免费在线看| 啦啦啦视频在线资源免费观看| 纯流量卡能插随身wifi吗| 国产成人精品一,二区| 国产精品.久久久| 亚洲精品久久久久久婷婷小说| 丰满人妻一区二区三区视频av| 一区二区三区四区激情视频| 国产精品一二三区在线看| 在线观看免费视频网站a站| 99视频精品全部免费 在线| 我的女老师完整版在线观看| 最后的刺客免费高清国语| 在线亚洲精品国产二区图片欧美 | 国产精品国产av在线观看| 精品熟女少妇av免费看| 亚洲精品色激情综合| 婷婷色综合www| 老熟女久久久| 免费观看av网站的网址| 2018国产大陆天天弄谢| 卡戴珊不雅视频在线播放| 午夜精品国产一区二区电影| 欧美激情国产日韩精品一区| 日韩av在线免费看完整版不卡| 免费久久久久久久精品成人欧美视频 | 久久精品国产亚洲av天美| 国产精品99久久久久久久久| 亚洲精品亚洲一区二区| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久久久免| 91精品一卡2卡3卡4卡| 欧美日韩在线观看h| 一二三四中文在线观看免费高清| 天堂俺去俺来也www色官网| 一个人看视频在线观看www免费| 午夜福利在线在线| 亚洲欧洲日产国产| 干丝袜人妻中文字幕| 色5月婷婷丁香| 亚洲真实伦在线观看| 久久综合国产亚洲精品| 日韩一区二区三区影片| 午夜视频国产福利| 26uuu在线亚洲综合色| 国产成人a∨麻豆精品| 国产一级毛片在线| www.av在线官网国产| 26uuu在线亚洲综合色| 欧美成人一区二区免费高清观看| av黄色大香蕉| 网址你懂的国产日韩在线| 日韩不卡一区二区三区视频在线| 久久97久久精品| 亚洲美女搞黄在线观看| 亚洲经典国产精华液单| 欧美bdsm另类| 日韩成人av中文字幕在线观看| 丝瓜视频免费看黄片| 欧美日韩一区二区视频在线观看视频在线| 日韩伦理黄色片| 女人十人毛片免费观看3o分钟| 国产精品久久久久久久久免| 国产在线视频一区二区| 亚洲精品日韩av片在线观看| 校园人妻丝袜中文字幕| 如何舔出高潮| 国产有黄有色有爽视频| 国产精品福利在线免费观看| 国产国拍精品亚洲av在线观看| 日韩人妻高清精品专区| 97在线视频观看| 欧美丝袜亚洲另类| 中国美白少妇内射xxxbb| 91久久精品电影网| 超碰av人人做人人爽久久| 一级爰片在线观看| 日韩人妻高清精品专区| 日韩av免费高清视频| 妹子高潮喷水视频| 国产在线视频一区二区| 免费少妇av软件| 中文字幕制服av| av福利片在线观看| 日本猛色少妇xxxxx猛交久久| 国产v大片淫在线免费观看| 亚洲经典国产精华液单| 麻豆国产97在线/欧美| 国产乱人偷精品视频| av在线观看视频网站免费| 女人十人毛片免费观看3o分钟| 成人美女网站在线观看视频| 婷婷色麻豆天堂久久| 亚洲内射少妇av| 中文欧美无线码| 91在线精品国自产拍蜜月| 欧美极品一区二区三区四区| 天天躁夜夜躁狠狠久久av| 国产黄片视频在线免费观看| 亚洲人成网站高清观看| 久久精品人妻少妇| 日本爱情动作片www.在线观看| 亚洲精品456在线播放app| 黄色日韩在线| 大又大粗又爽又黄少妇毛片口| 日韩成人伦理影院| 亚洲熟女精品中文字幕| 久久国产亚洲av麻豆专区| 尤物成人国产欧美一区二区三区| 日本免费在线观看一区| 高清av免费在线| 国产乱来视频区| 日本黄色日本黄色录像| 菩萨蛮人人尽说江南好唐韦庄| 一级毛片我不卡| 在线观看一区二区三区激情| 日韩av不卡免费在线播放| 偷拍熟女少妇极品色| 国产精品蜜桃在线观看| 国产av精品麻豆| 亚洲内射少妇av| 婷婷色av中文字幕| 久久久久久久久久人人人人人人| 欧美精品国产亚洲| 国产乱来视频区| 三级经典国产精品| 日韩中文字幕视频在线看片 | 日韩成人伦理影院| 3wmmmm亚洲av在线观看| 美女cb高潮喷水在线观看| 美女中出高潮动态图| 久久久午夜欧美精品| 嫩草影院新地址| 精品亚洲成a人片在线观看 | 日韩免费高清中文字幕av| 99久久精品一区二区三区| 观看av在线不卡| 最近手机中文字幕大全| 亚洲精华国产精华液的使用体验| 哪个播放器可以免费观看大片| 国产精品久久久久久久电影| 舔av片在线| 蜜桃久久精品国产亚洲av| 在线观看免费日韩欧美大片 | 亚洲第一区二区三区不卡| 精品一区在线观看国产| 熟妇人妻不卡中文字幕| 日韩欧美 国产精品| 天堂8中文在线网| 国产日韩欧美亚洲二区| 日韩一本色道免费dvd| 国产av国产精品国产| 又黄又爽又刺激的免费视频.| 久久久亚洲精品成人影院| 亚洲色图综合在线观看| 精品人妻一区二区三区麻豆| 亚洲精品日韩在线中文字幕| 国产免费视频播放在线视频| 男女无遮挡免费网站观看| 亚洲欧美日韩无卡精品| 日本一二三区视频观看| kizo精华| 亚洲综合色惰| 亚洲欧美日韩另类电影网站 | 一区在线观看完整版| 日韩免费高清中文字幕av| 日日摸夜夜添夜夜添av毛片| 日韩不卡一区二区三区视频在线| 男人和女人高潮做爰伦理| 亚洲,欧美,日韩| 亚洲精品乱久久久久久| 久久精品久久久久久久性| 国产有黄有色有爽视频| 激情 狠狠 欧美| 国产久久久一区二区三区| 成人毛片a级毛片在线播放| 国产黄片视频在线免费观看| 欧美一级a爱片免费观看看| 亚洲美女黄色视频免费看| 午夜老司机福利剧场| 18+在线观看网站| 久久久成人免费电影| 国产中年淑女户外野战色| 黄色日韩在线| 日韩强制内射视频| 高清在线视频一区二区三区| freevideosex欧美| 久久午夜福利片| 蜜桃在线观看..| 中国三级夫妇交换| 一级av片app| 男女免费视频国产| 一区二区三区四区激情视频| 日韩视频在线欧美| 黄片wwwwww| 欧美精品亚洲一区二区| av福利片在线观看| 国产高清国产精品国产三级 | 最新中文字幕久久久久| 亚洲欧美一区二区三区国产| 亚洲av二区三区四区| 亚洲精品国产成人久久av| 尤物成人国产欧美一区二区三区| 只有这里有精品99| 国产精品久久久久成人av| 啦啦啦中文免费视频观看日本| 亚洲av综合色区一区| 国产亚洲5aaaaa淫片| 男人爽女人下面视频在线观看| 女性生殖器流出的白浆| 深爱激情五月婷婷| 亚洲精华国产精华液的使用体验| 在线观看国产h片| 亚洲,欧美,日韩| 黄色怎么调成土黄色| 不卡视频在线观看欧美| 免费av中文字幕在线| 一级黄片播放器| 亚洲欧美精品自产自拍| 嫩草影院入口| 成人免费观看视频高清| 高清视频免费观看一区二区| 插阴视频在线观看视频| 国产精品欧美亚洲77777| 晚上一个人看的免费电影| 最新中文字幕久久久久| 青春草亚洲视频在线观看| 国产午夜精品一二区理论片| 国产亚洲午夜精品一区二区久久| 一级毛片aaaaaa免费看小| 成人特级av手机在线观看| 91精品一卡2卡3卡4卡| 高清不卡的av网站| 人妻少妇偷人精品九色| 国产爱豆传媒在线观看| 在线 av 中文字幕| 99视频精品全部免费 在线| 国产亚洲最大av| 久久这里有精品视频免费| 亚洲精品国产成人久久av| 不卡视频在线观看欧美| 久久午夜福利片| 人体艺术视频欧美日本| 能在线免费看毛片的网站| 大陆偷拍与自拍| 成人影院久久| 岛国毛片在线播放| 国产在视频线精品| 久久精品人妻少妇| 97精品久久久久久久久久精品| 国产精品偷伦视频观看了| 在线亚洲精品国产二区图片欧美 | 亚洲精品自拍成人| 麻豆成人午夜福利视频| 亚洲熟女精品中文字幕| 一个人免费看片子| 中文精品一卡2卡3卡4更新| 国产精品国产av在线观看| 免费观看性生交大片5| 人体艺术视频欧美日本| 国产精品国产三级国产专区5o| 日韩不卡一区二区三区视频在线| 又黄又爽又刺激的免费视频.| 少妇的逼好多水| 中文资源天堂在线| av国产久精品久网站免费入址| 交换朋友夫妻互换小说| 欧美 日韩 精品 国产| 久久久久久久久大av| 亚洲人成网站在线观看播放| 嫩草影院新地址| 亚洲第一区二区三区不卡| 99热这里只有是精品在线观看| 一区二区av电影网| 成人午夜精彩视频在线观看| 卡戴珊不雅视频在线播放| 国产精品久久久久久久电影| 国产视频内射| 亚洲精品日韩在线中文字幕| 国产美女午夜福利| 最近中文字幕2019免费版| 久久久国产一区二区| 三级国产精品片| 性高湖久久久久久久久免费观看| www.av在线官网国产| 国产成人一区二区在线| 久久久久久伊人网av| 国产深夜福利视频在线观看| 3wmmmm亚洲av在线观看| 麻豆乱淫一区二区| 成人免费观看视频高清| 免费观看性生交大片5| 国产成人一区二区在线| 舔av片在线| 午夜免费观看性视频| 久久久久久伊人网av| 人妻系列 视频| 成人毛片60女人毛片免费| 天堂8中文在线网| 久久久欧美国产精品| 亚洲国产欧美在线一区| 久久久久久人妻| 只有这里有精品99| 下体分泌物呈黄色| 在线看a的网站| 亚洲av免费高清在线观看| 人妻制服诱惑在线中文字幕| 国产综合精华液| 亚洲性久久影院| 午夜激情久久久久久久| 成年av动漫网址| 国产精品福利在线免费观看| 我要看黄色一级片免费的| 亚洲欧美成人精品一区二区| 男女边摸边吃奶| 99九九线精品视频在线观看视频| 国产大屁股一区二区在线视频| 国产爱豆传媒在线观看| 国产免费一区二区三区四区乱码| 少妇的逼水好多| 精品午夜福利在线看| 国产乱来视频区| 欧美另类一区| tube8黄色片| av免费在线看不卡| 欧美高清性xxxxhd video| 亚洲精品日韩在线中文字幕| 97精品久久久久久久久久精品| 日韩人妻高清精品专区| 国产老妇伦熟女老妇高清| 亚洲欧美一区二区三区黑人 | 一级毛片我不卡| 亚洲真实伦在线观看| 国产色爽女视频免费观看| 少妇裸体淫交视频免费看高清| 免费观看av网站的网址| 婷婷色av中文字幕| 啦啦啦在线观看免费高清www| 国产在线免费精品| 亚洲成人手机| 街头女战士在线观看网站| 黄色欧美视频在线观看| 少妇猛男粗大的猛烈进出视频| 日韩人妻高清精品专区| 在线观看人妻少妇| 亚洲欧美日韩卡通动漫| 内射极品少妇av片p| 精品一区二区三区视频在线|