• <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
    .国产精品久久| 免费av观看视频| 97超碰精品成人国产| 亚洲国产高清在线一区二区三| 亚洲精品日韩av片在线观看| 国产黄a三级三级三级人| 午夜老司机福利剧场| 97精品久久久久久久久久精品| 久久韩国三级中文字幕| av专区在线播放| 精品少妇黑人巨大在线播放| 日韩国内少妇激情av| 高清视频免费观看一区二区| 尤物成人国产欧美一区二区三区| 欧美成人精品欧美一级黄| 亚洲国产精品国产精品| 在线观看一区二区三区激情| 97超碰精品成人国产| 夫妻性生交免费视频一级片| 国产男女超爽视频在线观看| 人妻夜夜爽99麻豆av| 国产午夜精品久久久久久一区二区三区| 久久久久久伊人网av| 一级a做视频免费观看| 国产高清有码在线观看视频| 性色avwww在线观看| 日韩 亚洲 欧美在线| 久久久久久久国产电影| 男女那种视频在线观看| 日韩国内少妇激情av| 午夜福利在线观看免费完整高清在| 国产成人91sexporn| 一二三四中文在线观看免费高清| 久久精品国产a三级三级三级| 国产精品.久久久| 国产一区二区在线观看日韩| 能在线免费看毛片的网站| 丝瓜视频免费看黄片| 日韩国内少妇激情av| 亚洲熟女精品中文字幕| 日日摸夜夜添夜夜爱| 91精品国产九色| 国产精品99久久久久久久久| 精华霜和精华液先用哪个| 在线观看美女被高潮喷水网站| 三级国产精品片| 午夜免费男女啪啪视频观看| 日本与韩国留学比较| 18禁在线播放成人免费| 国产欧美日韩一区二区三区在线 | 精品一区在线观看国产| 美女高潮的动态| 欧美日本视频| 少妇人妻一区二区三区视频| 久久人人爽av亚洲精品天堂 | 人妻 亚洲 视频| 永久免费av网站大全| 亚洲国产欧美在线一区| 亚洲人与动物交配视频| 人人妻人人爽人人添夜夜欢视频 | 亚洲性久久影院| 欧美亚洲 丝袜 人妻 在线| 91久久精品国产一区二区三区| 国产成人aa在线观看| av专区在线播放| 男人舔奶头视频| 男女国产视频网站| 免费看日本二区| 欧美激情国产日韩精品一区| 国产伦精品一区二区三区视频9| 欧美zozozo另类| 成人亚洲精品av一区二区| 黄色一级大片看看| av免费在线看不卡| 在线 av 中文字幕| 在线播放无遮挡| 久久久午夜欧美精品| 久久ye,这里只有精品| 久久人人爽人人片av| 日韩在线高清观看一区二区三区| eeuss影院久久| 久久6这里有精品| 国产欧美日韩一区二区三区在线 | 国产成人aa在线观看| 免费看不卡的av| 汤姆久久久久久久影院中文字幕| 卡戴珊不雅视频在线播放| 成年av动漫网址| 成年版毛片免费区| 免费不卡的大黄色大毛片视频在线观看| 男插女下体视频免费在线播放| 老女人水多毛片| av在线蜜桃| 国产男女内射视频| 超碰97精品在线观看| 免费看光身美女| 久久影院123| 狂野欧美白嫩少妇大欣赏| 国产精品国产av在线观看| 午夜日本视频在线| 亚洲国产精品成人综合色| 高清在线视频一区二区三区| 亚洲欧美精品自产自拍| 舔av片在线| 啦啦啦啦在线视频资源| 久久久久久久久久人人人人人人| 国产成人福利小说| 女的被弄到高潮叫床怎么办| 肉色欧美久久久久久久蜜桃 | 久久久久久久久久成人| 国产精品秋霞免费鲁丝片| 午夜福利视频1000在线观看| 少妇的逼水好多| 亚洲真实伦在线观看| 午夜亚洲福利在线播放| 国产日韩欧美亚洲二区| 亚洲精品乱码久久久久久按摩| 中国国产av一级| 超碰97精品在线观看| 国产精品.久久久| 我的女老师完整版在线观看| 国产中年淑女户外野战色| 精品少妇黑人巨大在线播放| 日韩不卡一区二区三区视频在线| 水蜜桃什么品种好| 国产成人精品福利久久| 精品久久久久久久久av| 亚洲精品中文字幕在线视频 | 中文字幕人妻熟人妻熟丝袜美| 日本黄大片高清| 国产精品三级大全| 成人二区视频| 1000部很黄的大片| 日本av手机在线免费观看| 亚洲图色成人| 亚洲精品视频女| 少妇裸体淫交视频免费看高清| 成人亚洲精品av一区二区| 一级毛片黄色毛片免费观看视频| 舔av片在线| 久久久久久久国产电影| 国产永久视频网站| 国产精品久久久久久av不卡| 少妇的逼水好多| av免费观看日本| 国产免费一区二区三区四区乱码| 亚洲欧美成人精品一区二区| 国产精品久久久久久精品电影| 久久久久网色| 久久久久国产网址| 亚洲精品乱久久久久久| 成人二区视频| 亚洲精品aⅴ在线观看| 舔av片在线| 成年av动漫网址| 亚洲欧美精品自产自拍| 日日摸夜夜添夜夜爱| 国产精品麻豆人妻色哟哟久久| 亚洲精品亚洲一区二区| 国产免费一区二区三区四区乱码| 亚洲国产精品999| 日韩精品有码人妻一区| 狂野欧美激情性bbbbbb| 免费av不卡在线播放| 国产男人的电影天堂91| 极品教师在线视频| 日韩欧美 国产精品| 国国产精品蜜臀av免费| 嘟嘟电影网在线观看| 99视频精品全部免费 在线| 国产亚洲午夜精品一区二区久久 | 免费电影在线观看免费观看| 热re99久久精品国产66热6| 99久久九九国产精品国产免费| 日韩大片免费观看网站| 一级av片app| 欧美日韩综合久久久久久| 九色成人免费人妻av| 欧美3d第一页| 国产av码专区亚洲av| 亚洲国产欧美人成| 国产 一区 欧美 日韩| 久久精品人妻少妇| 欧美日韩国产mv在线观看视频 | 亚洲精品自拍成人| 色播亚洲综合网| 男女无遮挡免费网站观看| 国产精品成人在线| 在线 av 中文字幕| 五月开心婷婷网| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品第二区| 国产成人精品福利久久| 搡老乐熟女国产| 亚洲不卡免费看| 国产乱人视频| 亚洲va在线va天堂va国产| 久热这里只有精品99| 精华霜和精华液先用哪个| 五月天丁香电影| 日本熟妇午夜| 欧美另类一区| 99热这里只有是精品50| 在线观看一区二区三区| 午夜激情福利司机影院| 国产精品久久久久久久电影| 久久久久久久午夜电影| 国产日韩欧美在线精品| 黄片无遮挡物在线观看| 麻豆成人午夜福利视频| 老司机影院毛片| 综合色av麻豆| 男女那种视频在线观看| 最新中文字幕久久久久| 日韩一区二区三区影片| 国产高潮美女av| 免费观看性生交大片5| 亚洲在线观看片| 久久国产乱子免费精品| 99热全是精品| 亚洲欧美精品自产自拍| 国模一区二区三区四区视频| 国产欧美亚洲国产| 久久久精品欧美日韩精品| 最近最新中文字幕大全电影3| 国产精品无大码| 亚洲最大成人手机在线| 久久久久久久久久久免费av| 免费大片18禁| 国产又色又爽无遮挡免| 大码成人一级视频| 久久精品国产鲁丝片午夜精品| 亚洲一区二区三区欧美精品 | 午夜福利在线在线| 在线 av 中文字幕| 日本免费在线观看一区| 午夜福利视频精品| 国产精品99久久久久久久久| 国产高清国产精品国产三级 | 欧美日韩一区二区视频在线观看视频在线 | 男人舔奶头视频| 亚洲va在线va天堂va国产| 欧美人与善性xxx| 精品人妻视频免费看| 毛片女人毛片| 看免费成人av毛片| 久久影院123| 国产日韩欧美在线精品| 99久久中文字幕三级久久日本| 黄色视频在线播放观看不卡| 人人妻人人看人人澡| 深夜a级毛片| 亚洲欧美日韩卡通动漫| 中文字幕免费在线视频6| 人体艺术视频欧美日本| 九九久久精品国产亚洲av麻豆| 久久久久久久精品精品| 热re99久久精品国产66热6| 人妻制服诱惑在线中文字幕| 精华霜和精华液先用哪个| 卡戴珊不雅视频在线播放| 久久久久久久久久久丰满| 亚洲电影在线观看av| 大码成人一级视频| 亚洲成人久久爱视频| 亚洲欧美成人精品一区二区| 免费高清在线观看视频在线观看| 亚洲精品国产色婷婷电影| 亚洲欧美日韩另类电影网站 | 好男人在线观看高清免费视频| 欧美+日韩+精品| 亚洲av福利一区| 国产91av在线免费观看| 在线观看av片永久免费下载| 久久久久九九精品影院| 欧美日韩综合久久久久久| 日本爱情动作片www.在线观看| 青春草视频在线免费观看| 26uuu在线亚洲综合色| 麻豆成人午夜福利视频| 日本黄大片高清| 精品久久久久久久末码| 岛国毛片在线播放| 免费观看性生交大片5| 欧美高清成人免费视频www| 国产成人免费观看mmmm| 亚洲精品乱码久久久v下载方式| av在线观看视频网站免费| 春色校园在线视频观看| 成人免费观看视频高清| 啦啦啦中文免费视频观看日本| 波多野结衣巨乳人妻| 国产精品国产三级国产av玫瑰| 欧美 日韩 精品 国产| 国产乱人视频| 51国产日韩欧美| 免费观看av网站的网址| 精品久久久精品久久久| 最近中文字幕2019免费版| 嫩草影院新地址| 国产伦精品一区二区三区四那| 我的女老师完整版在线观看| videos熟女内射| 精品国产乱码久久久久久小说| 韩国高清视频一区二区三区| 亚洲av中文av极速乱| 一本一本综合久久| 亚洲电影在线观看av| 久久午夜福利片| 日韩国内少妇激情av| 听说在线观看完整版免费高清| 久久这里有精品视频免费| 精品午夜福利在线看| 久久精品国产鲁丝片午夜精品| eeuss影院久久| 久久99热6这里只有精品| 久久久亚洲精品成人影院| 大香蕉97超碰在线| 身体一侧抽搐| 亚洲av免费高清在线观看| 最近中文字幕2019免费版| 一个人观看的视频www高清免费观看| 欧美+日韩+精品| 嫩草影院精品99| 亚洲丝袜综合中文字幕| 一边亲一边摸免费视频| 成年女人在线观看亚洲视频 | 日韩在线高清观看一区二区三区| 搡老乐熟女国产| 小蜜桃在线观看免费完整版高清| 日韩强制内射视频| 亚洲精华国产精华液的使用体验| 我要看日韩黄色一级片| 国产成人一区二区在线| 免费大片18禁| 日韩,欧美,国产一区二区三区| 一级毛片久久久久久久久女| 亚洲在久久综合| 偷拍熟女少妇极品色| 在线播放无遮挡| 91久久精品国产一区二区三区| 午夜免费观看性视频| 国产男人的电影天堂91| 欧美潮喷喷水| 亚洲精品国产色婷婷电影| 人妻 亚洲 视频| 精品久久国产蜜桃| 久久久欧美国产精品| 国产黄a三级三级三级人| 99热这里只有是精品50| 久久久国产一区二区| 国产免费又黄又爽又色| 午夜激情久久久久久久| 精品视频人人做人人爽| 精品国产乱码久久久久久小说| 精品人妻熟女av久视频| av在线播放精品| 久久精品夜色国产| 别揉我奶头 嗯啊视频| 人妻夜夜爽99麻豆av| 91久久精品国产一区二区成人| 日韩一区二区视频免费看| 又大又黄又爽视频免费| 毛片一级片免费看久久久久| 国产色婷婷99| 国产一区二区三区av在线| 一级毛片 在线播放| av女优亚洲男人天堂| 色综合色国产| 自拍欧美九色日韩亚洲蝌蚪91 | av在线观看视频网站免费| 国语对白做爰xxxⅹ性视频网站| 国产精品国产av在线观看| 亚洲av男天堂| 天堂俺去俺来也www色官网| 色网站视频免费| 久久久久久久精品精品| 中文字幕久久专区| 日韩不卡一区二区三区视频在线| 亚洲精品第二区| 亚洲第一区二区三区不卡| 天美传媒精品一区二区| 亚洲av成人精品一区久久| 少妇人妻一区二区三区视频| 水蜜桃什么品种好| 久久99精品国语久久久| 又爽又黄无遮挡网站| 97人妻精品一区二区三区麻豆| 国产精品久久久久久久久免| 精品久久久久久久久亚洲| 欧美 日韩 精品 国产| 日韩伦理黄色片| 在线亚洲精品国产二区图片欧美 | 在线观看一区二区三区激情| 一级片'在线观看视频| 夜夜爽夜夜爽视频| 国国产精品蜜臀av免费| 麻豆成人av视频| 午夜亚洲福利在线播放| 亚洲经典国产精华液单| 亚洲最大成人手机在线| 男女啪啪激烈高潮av片| 成年免费大片在线观看| 日韩中字成人| 国产91av在线免费观看| 少妇人妻久久综合中文| 又黄又爽又刺激的免费视频.| 黑人高潮一二区| 成人欧美大片| 免费大片18禁| 一个人看的www免费观看视频| 日韩欧美精品免费久久| 国产人妻一区二区三区在| 午夜福利网站1000一区二区三区| 秋霞在线观看毛片| 久久国产乱子免费精品| 国产乱人偷精品视频| 各种免费的搞黄视频| 日韩视频在线欧美| 久久久精品欧美日韩精品| 久久精品国产亚洲av涩爱| 久久99热这里只频精品6学生| 国产欧美亚洲国产| 日韩免费高清中文字幕av| 国产精品一区二区性色av| 国产美女午夜福利| 久久国产乱子免费精品| 极品少妇高潮喷水抽搐| 又粗又硬又长又爽又黄的视频| 成人欧美大片| 午夜亚洲福利在线播放| 久久久久久久久久人人人人人人| 男人爽女人下面视频在线观看| 午夜日本视频在线| 亚洲国产成人一精品久久久| 国产精品国产三级国产av玫瑰| 联通29元200g的流量卡| 麻豆国产97在线/欧美| 久久久欧美国产精品| 视频中文字幕在线观看| 亚洲成人av在线免费| 免费观看的影片在线观看| 两个人的视频大全免费| 97精品久久久久久久久久精品| 国产一区二区三区av在线| 久久久久性生活片| 男的添女的下面高潮视频| 亚洲国产精品成人综合色| 免费看不卡的av| 成人国产av品久久久| 国产v大片淫在线免费观看| 国产探花在线观看一区二区| 免费av毛片视频| 麻豆精品久久久久久蜜桃| 亚洲av欧美aⅴ国产| videossex国产| 少妇高潮的动态图| 欧美日韩视频高清一区二区三区二| 人人妻人人看人人澡| 色哟哟·www| 在线观看人妻少妇| 久久久久久伊人网av| 九九爱精品视频在线观看| 午夜日本视频在线| 国模一区二区三区四区视频| 欧美高清成人免费视频www| 亚洲精品亚洲一区二区| 国产成人精品婷婷| 在线观看人妻少妇| 一级毛片 在线播放| 国产日韩欧美亚洲二区| 国产精品秋霞免费鲁丝片| 国产爽快片一区二区三区| 亚洲天堂av无毛| 干丝袜人妻中文字幕| 亚洲av中文av极速乱| 国产精品嫩草影院av在线观看| 精品国产一区二区三区久久久樱花 | 久久精品久久久久久噜噜老黄| tube8黄色片| 国产亚洲午夜精品一区二区久久 | 六月丁香七月| 久久久久久九九精品二区国产| 国产一级毛片在线| 日韩三级伦理在线观看| 夜夜爽夜夜爽视频| 搡老乐熟女国产| 国产白丝娇喘喷水9色精品| 禁无遮挡网站| 国产精品偷伦视频观看了| 精品99又大又爽又粗少妇毛片| 亚洲一级一片aⅴ在线观看| 国产午夜精品一二区理论片| 日日摸夜夜添夜夜爱| 亚洲电影在线观看av| 狂野欧美白嫩少妇大欣赏| 少妇被粗大猛烈的视频| 天堂网av新在线| 免费大片18禁| 亚洲国产高清在线一区二区三| 免费在线观看成人毛片| 亚洲欧美日韩无卡精品| 极品教师在线视频| 少妇丰满av| 女的被弄到高潮叫床怎么办| 亚洲真实伦在线观看| 国产人妻一区二区三区在| 寂寞人妻少妇视频99o| 国精品久久久久久国模美| 色综合色国产| 能在线免费看毛片的网站| 国产欧美亚洲国产| 噜噜噜噜噜久久久久久91| 久久久欧美国产精品| 国产男女内射视频| 日韩欧美一区视频在线观看 | 国产精品一区二区三区四区免费观看| 色视频在线一区二区三区| 不卡视频在线观看欧美| 成人漫画全彩无遮挡| 免费看不卡的av| 深爱激情五月婷婷| 人妻一区二区av| 亚洲av电影在线观看一区二区三区 | 亚洲真实伦在线观看| 啦啦啦啦在线视频资源| 免费观看a级毛片全部| 亚洲成人一二三区av| 成人黄色视频免费在线看| 香蕉精品网在线| 亚洲欧美成人综合另类久久久| 久久久久久久久久久丰满| 成人特级av手机在线观看| 一个人看的www免费观看视频| 亚洲欧美一区二区三区国产| 亚洲怡红院男人天堂| 人人妻人人看人人澡| 好男人视频免费观看在线| 国产永久视频网站| 高清毛片免费看| 激情五月婷婷亚洲| 99精国产麻豆久久婷婷| 日本黄色片子视频| 精品午夜福利在线看| 国产探花极品一区二区| 精品一区二区三区视频在线| 国产一区二区在线观看日韩| 99久久精品一区二区三区| 免费观看无遮挡的男女| 日韩av不卡免费在线播放| 99久久人妻综合| 亚洲精品国产色婷婷电影| 欧美日韩综合久久久久久| 天堂俺去俺来也www色官网| 亚洲精品乱久久久久久| 免费看光身美女| 日韩亚洲欧美综合| 国产一区二区在线观看日韩| 人妻 亚洲 视频| 国产免费一区二区三区四区乱码| 亚洲真实伦在线观看| 亚洲人成网站高清观看| 国产高清有码在线观看视频| 国内少妇人妻偷人精品xxx网站| 蜜桃久久精品国产亚洲av| 国产视频首页在线观看| 国产老妇女一区| 国产视频内射| 777米奇影视久久| 少妇高潮的动态图| 99热6这里只有精品| 中国三级夫妇交换| 久久鲁丝午夜福利片| 精品人妻熟女av久视频| 免费av毛片视频| 高清视频免费观看一区二区| 2021少妇久久久久久久久久久| 亚洲精品成人久久久久久| 国产淫语在线视频| 国产成人a∨麻豆精品| 男人舔奶头视频| 国产色爽女视频免费观看| 亚洲av一区综合| 三级国产精品欧美在线观看| 欧美成人午夜免费资源| 亚洲av电影在线观看一区二区三区 | 91午夜精品亚洲一区二区三区| 丝瓜视频免费看黄片| 成人漫画全彩无遮挡| 三级国产精品欧美在线观看| 欧美成人午夜免费资源| 欧美性猛交╳xxx乱大交人| 在线看a的网站| 熟妇人妻不卡中文字幕| 欧美性猛交╳xxx乱大交人| 亚洲av日韩在线播放| 日本三级黄在线观看| 免费大片18禁| 午夜福利高清视频| 亚洲av在线观看美女高潮| 成人漫画全彩无遮挡| 1000部很黄的大片| 国产乱人偷精品视频| 亚洲av电影在线观看一区二区三区 | 亚洲,一卡二卡三卡| 日韩人妻高清精品专区| av在线app专区| 观看美女的网站| 成人毛片a级毛片在线播放| 成年av动漫网址| 免费大片18禁| 在线观看av片永久免费下载| 观看免费一级毛片| 国产免费视频播放在线视频| 亚洲精品久久久久久婷婷小说|