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

    磺胺嘧啶在覆鐵石英砂中遷移特征及數(shù)值模擬

    2019-11-28 10:50:18徐紹輝
    中國環(huán)境科學(xué) 2019年11期
    關(guān)鍵詞:流液石英砂離子

    張 惠,林 青,徐紹輝

    磺胺嘧啶在覆鐵石英砂中遷移特征及數(shù)值模擬

    張 惠,林 青,徐紹輝*

    (青島大學(xué)環(huán)境科學(xué)與工程學(xué)院,山東 青島 266071)

    為了探討磺胺類抗生素在土壤中的遷移特征,以磺胺嘧啶(SDZ)為研究對象,以石英砂(QS)為介質(zhì),采用高溫加熱法制備覆鐵石英砂(IOCS),通過土柱出流實驗,研究不同鐵含量和離子強度對SDZ在IOCS中遷移的影響,并利用Hydrus-1D軟件對結(jié)果進行模擬.土柱出流實驗結(jié)果表明,隨著IOCS表面鐵含量增加,其表面粗糙度越大,且表面攜帶的正電荷越多,SDZ越容易吸附在IOCS表面,其相對濃度峰值越低,不利于SDZ的遷移;離子強度的增加,促進了SDZ在QS和IOCS中的遷移,這主要是受靜電作用、表面絡(luò)合作用和競爭吸附作用的控制.模擬結(jié)果顯示,非平衡單點吸附模型(OSM)可以較好地模擬不同鐵含量和離子強度下SDZ的遷移過程,擬合曲線與實測曲線幾乎重合,2>0.97,RMSE<0.032,擬合效果較好;反演所得吸附分配系數(shù)F隨著鐵含量的增加而增加,隨離子強度的增加而減小;阻滯因子d與其呈現(xiàn)相同的變化趨勢;吸附特征參數(shù)>1,而一階動力學(xué)速率常數(shù)均較小.

    磺胺嘧啶;覆鐵石英砂;遷移;Hydrus-1D;鐵含量;離子強度

    抗生素作為一種新型環(huán)境有機污染物,其在環(huán)境中的吸附遷移規(guī)律備受關(guān)注[1-2].進入環(huán)境中的抗生素通過多種途徑匯集于土壤,對土壤中微生物活性、動植物生長產(chǎn)生較大的危害,一些吸附性弱、遷移性強的抗生素則通過地表徑流、淋濾作用進入水體環(huán)境,從而污染地表水與地下水[3-6].磺胺嘧啶(SDZ)是磺胺類抗生素的一種,其在酸性和堿性條件下均易溶解,更易遷移,對環(huán)境的危害更大[7].

    土壤氧化鐵礦物比表面積大,傳輸電子能力強,并且具有較高的氧化還原活性,一些離子型有機污染物更趨于吸附在其表面,且受大量物理、化學(xué)及生物等因素的影響.當(dāng)前已有不少關(guān)于抗生素類化合物在氧化鐵礦物上吸附特征的研究,吸附機制也不盡相同[8-11].為分析單一因素下SDZ在含氧化鐵土壤中的遷移特征,本研究以低反應(yīng)性石英砂為介質(zhì),采用高溫加熱法制備覆鐵石英砂(Iron Oxide Coated Sand,IOCS).據(jù)報道,鐵在IOCS表面以氧化鐵形式存在,比表面積提高5~10倍,且表面主要帶正電荷,吸附容量大大增加[12].目前已有大量關(guān)于納米材料[13-14]、膠體[15-16]、重金屬[17-18]和有機物[19-20]等在IOCS中遷移規(guī)律的研究,而關(guān)于磺胺類抗生素在IOCS中遷移特征的研究鮮有報道.

    Hydrus-1D是用于模擬一維飽和與非飽和多孔介質(zhì)中水流、溶質(zhì)運移以及熱量傳輸?shù)能浖?該軟件用Richards方程表征水流運動過程,用對流-彌散方程描述溶質(zhì)運移過程;采用Galerkin線性有限元法對方程進行求解,并可反演水流、溶質(zhì)在多孔介質(zhì)中的水力傳導(dǎo)率、吸附、解吸、降解和彌散等相關(guān)參數(shù),對水流和溶質(zhì)遷移特征進行量化再現(xiàn),具有靈活的輸入輸出功能,在土壤水分運動和溶質(zhì)運移方面具有廣泛的適用性[21].

    本文以SDZ為研究對象,采用室內(nèi)土柱出流實驗,研究不同鐵含量和離子強度下SDZ在石英砂中的遷移特征,并采用Hydrus-1D對結(jié)果進行模擬,以探討其在多孔介質(zhì)中的遷移轉(zhuǎn)化機理,進而為SDZ的土壤及地下水污染修復(fù)提供理論依據(jù).

    1 材料與方法

    1.1 供試材料

    SDZ,純度99%,購自北京百靈威科技有限公司,石英砂(QS)、HCl、NaOH、FeCl3·6H2O、KBr、CaCl2均為分析純.

    選取20~40目石英砂,首先用0.1mol/L的HCl溶液浸泡24h,用去離子水反復(fù)沖洗,然后用0.1mol/L的NaOH溶液浸泡24h,用去離子水反復(fù)沖洗至pH值為中性,105℃烘干備用.

    IOCS的制備[22].取3份清洗過后的石英砂100g于直徑10cm的培養(yǎng)皿中,分別加入50mL 0.05,0.1, 0.5mol/L的FeCl3·6H2O溶液,邊攪拌邊加入3mL 5mol/L的NaOH溶液,調(diào)pH值為堿性,攪拌均勻后放入烘箱,于105℃高溫加熱,前6h每隔一段時間攪拌一次,防止石英砂顆?;ハ嗾尺B,直至溶液完全蒸干,24h后停止加熱并冷卻至室溫,用去離子水反復(fù)沖洗至顏色不變后于105℃烘干備用,制備的3種IOCS分別為IOCS1、IOCS2、IOCS3,其表面鐵含量分別為1.3,2.1,3.5mg/g.

    1.2 實驗方法

    1.2.1 土柱填裝 實驗土柱為高20cm,內(nèi)徑5cm的有機玻璃柱,分別取QS和不同鐵含量的IOCS約560g分5次進行填裝,每填完一層用玻璃棒輕微攪動并用壓實器壓實,以保證土柱填裝均勻,在土柱兩端接口處分別放置一層脫脂棉,然后放一層內(nèi)徑5cm的濾紙,使溶液均勻流入、流出,并防止砂粒流出堵塞進出水口.

    1.2.2 示蹤實驗 以KBr為示蹤劑,通過Br-穿透曲線來判斷土柱填裝的均勻性,并為SDZ運移模擬提供彌散系數(shù)等水力學(xué)參數(shù).首先用蠕動泵將0.01mol/L的CaCl2溶液以較低流速自下往上輸入分別裝有QS和3種IOCS的砂柱中,使砂柱飽和并排出空氣,然后轉(zhuǎn)變?nèi)肓鞣较?待流場穩(wěn)定后從上往下脈沖輸入1孔隙體積(PV)0.05mol/L的KBr溶液,最后再輸入0.01mol/L的CaCl2溶液,用自動部分收集器收集出流液,直至檢測不到Br-為止,Br-濃度用Br-選擇電極進行測定,實驗完成后繼續(xù)利用該組砂柱進行不同鐵含量下SDZ遷移實驗,每個實驗重復(fù)2次.其中孔隙體積

    PV=(1)

    式中:PV為無量綱時間;為孔隙水流速,cm/min;為時間,min;為土柱長度,cm.

    1.2.3 SDZ遷移實驗 在不同鐵含量實驗中,設(shè)置石英砂表面鐵含量分別為0,1.3,2.1,3.5mg/g,然后將pH=6.14,濃度為0.01mol/L的CaCl2溶液通過蠕動泵自下往上輸入分別裝有這四種砂的砂柱中進行飽和,約40h后轉(zhuǎn)變?nèi)肓鞣较?待流場穩(wěn)定后自上往下分別脈沖輸入30mg/L的SDZ溶液約600mL,然后再用0.01mol/L的CaCl2溶液分別進行洗脫.

    在不同離子強度實驗中,設(shè)置CaCl2濃度分別為0.001,0.005,0.01,0.05mol/L,用IOCS3填裝砂柱, QS砂柱作為對照,然后分別用pH=6.14上述濃度的CaCl2溶液飽和砂柱,SDZ輸入以及洗脫過程同上.

    用自動部分收集器收集出流液.用紫外分光光度計(UV9100C PC,北京萊博泰科儀器有限公司)測定出流液中SDZ的濃度[23],pH計(PHS-3C,上海精科儀器有限公司)測定pH值,ORP復(fù)合電極(雷磁 501,上海精密科技有限公司)測定氧化還原電位(h),便攜式電導(dǎo)率儀(雷磁 DDB-303A,上海精密科技有限公司)測定電導(dǎo)率(c),本實驗重復(fù)2次.

    1.3 運移模型

    局部平衡對流-彌散方程(LEA)描述穩(wěn)定流條件下,只考慮吸附作用時,溶質(zhì)在一維均勻介質(zhì)中的運移過程,且吸附過程瞬時可逆[24].方程如下:

    式中:為溶質(zhì)濃度mg/L;為時間,min;為水動力彌散系數(shù),cm2/min,;為彌散度,cm;為平均孔隙水流速,cm/min,s;為水流通量,cm/min;s為飽和含水量,cm3/cm3;為遷移距離,cm;為容重,g/cm3;為溶質(zhì)在單位質(zhì)量土壤上的吸附量, mg/g;為吸附量隨時間的變化,對于非反應(yīng)性Br-,為0,所以可以利用確定性平衡方程來反演出參數(shù)、s.

    非平衡單點吸附模型(OSM)假定吸附過程為動力學(xué)吸附,不存在瞬時吸附,Freundlich等溫吸附模型可以很好地描述SDZ的吸附過程[25].方程如下:

    式中:為溶質(zhì)濃度,mg/cm3;e為吸附平衡時單位質(zhì)量土壤吸附的溶質(zhì)濃度,mg/g;為動力學(xué)吸附的溶質(zhì)濃度,mg/cm3;為描述動力學(xué)吸附過程的一階速率常數(shù),min-1;F為吸附分配系數(shù),cm3/g.

    利用Freundlich等溫吸附模型線性化方法求線性分配系數(shù),并求得阻滯因子d[26]:

    式中:為線性吸附分配系數(shù),無單位;為吸附特征參數(shù);為容重,g/cm3;0為SDZ初始濃度,mg/L;為體積含水量,cm3/cm3.

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

    采用Microsoft Excel 2007進行數(shù)據(jù)處理,采用Origin 9.0、CorelDRAW X7軟件作圖,采用Hydrus-1D軟件對實驗結(jié)果進行模擬,并反演得到相關(guān)參數(shù).

    2 結(jié)果與討論

    2.1 鐵含量對SDZ在石英砂中遷移的影響

    圖1 不同鐵含量下SDZ的穿透曲線

    圖2 SDZ吸附作用示意

    圖1為不同鐵含量下SDZ的穿透曲線,圖2為SDZ吸附作用示意圖.可以看出,IOCS表面鐵含量越多,SDZ出流時間越晚,穿透曲線峰值越低,越不利于SDZ的遷移.一方面,隨著鐵含量的增加,IOCS表面粗糙度增大,從而比表面積也隨之增大[27],更有利于SDZ的吸附;另一方面,由于表面羥基化作用,氧化鐵在水溶液中可以發(fā)生表面離子化或電離化,從而使IOCS表面出現(xiàn)可變電荷[28].在一定pH值下,介質(zhì)表面正電荷與負電荷數(shù)量相等,表面凈電荷為零,此時的pH值稱為可變電荷的電荷零點(PZC).當(dāng)溶液pH pHPZC時,介質(zhì)表面可變電荷為負電荷[29].一般氧化鐵礦物的電荷零點在8.5左右[30],所以在pH=4~6.5范圍內(nèi),IOCS表面主要帶正電荷,表面鐵主要以水合陽離子形式存在,SDZ主要以中性分子兼少量陰離子形式存在,分子形態(tài)的SDZ具有更強的疏水性,且IOCS表面提供的疏水性吸附點位能更好的與其結(jié)合;由于非專性吸附作用,陰離子形態(tài)的SDZ在靜電引力作用下會與帶正電荷的IOCS產(chǎn)生靜電吸附,且IOCS表面鐵含量越多,所帶的正電荷越多,靜電引力越大,越有利于SDZ的吸附.

    由圖3可以看出,QS和IOCS出流液的pH值不同,且QS>IOCS1>IOCS2>IOCS3.原因可能是IOCS表面的水合氧化鐵易發(fā)生電離生成H+,隨著IOCS表面鐵含量的增加,生成的H+越多,從而溶液的pH值越低.

    出流液的氧化還原電位處于200~360mV之間,屬于氧化條件,整體呈現(xiàn)與pH值相反的規(guī)律,4種情況下氧化還原電位大小順序為:IOCS3>IOCS2> IOCS1>QS.

    出流液電導(dǎo)率整體呈下降的趨勢,電導(dǎo)率變化主要與溶液中Ca2+變化有關(guān),可能是因為遷移過程中,溶液中大量Ca2+吸附在IOCS表面,導(dǎo)致出流液中Ca2+離子濃度有所減少,從而引起電導(dǎo)率下降.

    2.2 離子強度對SDZ在石英砂中遷移的影響

    由圖4可知,同一離子強度下,SDZ在IOCS中遷移的穿透曲線峰值均低于QS,說明氧化鐵的存在抑制了SDZ的遷移.不同離子強度下,隨著CaCl2濃度的增加,SDZ在QS和IOCS中遷移的穿透曲線峰值均增大,說明離子強度的增加促進了SDZ的遷移.QS表面為負電荷,在出流液pH=6.2~7.4范圍內(nèi),SDZ主要以陰離子和中性分子形態(tài)存在,帶負電荷的石英砂與SDZ陰離子產(chǎn)生靜電斥力,不利于SDZ陰離子吸附,中性SDZ分子可能更易吸附在QS表面,但隨著離子強度的增加,越來越多Ca2+與QS結(jié)合,與SDZ中性分子競爭吸附點位,從而促進SDZ遷移;對于IOCS,在出流液pH=3.5~6.5范圍內(nèi),其表面帶正電荷,此時SDZ主要以中性兼少量的陰離子存在,IOCS表面水合氧化鐵的羥基在酸性條件下會產(chǎn)生去質(zhì)子化位點,與SDZ中性分子結(jié)合形成外層表面絡(luò)合物,外層表面絡(luò)合物并不穩(wěn)定,隨著離子強度的增加,大量Ca2+會與水合氧化鐵的羥基結(jié)合吸附在IOCS表面,與SDZ形成競爭吸附,從而促進SDZ遷移;少量SDZ陰離子會與帶正電荷的IOCS產(chǎn)生靜電引力,隨著離子強度的增加,Ca2+與SDZ陰離子作用力加強,從而減少靜電引力作用,促進SDZ遷移.這與不同CaCl2濃度下諾氟沙星在土壤中吸附規(guī)律相似[31].

    圖5為不同離子強度下出流液中pH值、h、c變化特征.SDZ在QS中遷移時,出流液pH值先增大,在1PV左右開始趨于穩(wěn)定,在4PV左右開始降低,與SDZ出流規(guī)律一致,原因可能是在遷移過程中,SDZ陰離子與H+結(jié)合隨出流液流出,從而使出流液中H+數(shù)量減少,pH值先上升并趨于穩(wěn)定,當(dāng)SDZ濃度開始減小時,與H+結(jié)合的SDZ陰離子數(shù)量減少,從而出流液中H+數(shù)量開始增加,pH值隨之降低;SDZ在IOCS中遷移時,出流液pH值呈降低趨勢,且最后均穩(wěn)定在4左右,原因可能是IOCS表面氧化鐵發(fā)生電離,產(chǎn)生大量H+,從而使溶液中H+含量增加,其數(shù)量遠大于與SDZ陰離子結(jié)合的H+,所以出流液pH值持續(xù)降低并趨于穩(wěn)定.

    由氧化還原電位變化可以看出,所有實驗均在氧化性條件下進行,且氧化還原電位整體變化趨勢與pH值變化趨勢相反,SDZ在QS中遷移時,出流液氧化還原電位先降低然后再升高,在IOCS中遷移時,氧化還原電位則呈上升趨勢.

    由電導(dǎo)率曲線可以看出,隨著離子強度的增強,電導(dǎo)率增大,這是因為Ca2+對電導(dǎo)率變化起主導(dǎo)作用,離子強度越高,出流液中Ca2+濃度越大,從而電導(dǎo)率越大.

    2.3 Br-穿透曲線及參數(shù)確定

    圖6 Br-在QS和IOCS中遷移的穿透曲線

    表1 Br-穿透曲線擬合相關(guān)參數(shù)

    圖6為Br-在QS和IOCS中遷移的穿透曲線及擬合曲線,表1為利用LEA擬合Br-穿透曲線得到的相關(guān)參數(shù).Br-穿透曲線對稱性良好,不存在拖尾現(xiàn)象,/0最大值都能達到1,說明土柱填裝比較均勻,不存在大孔隙及優(yōu)先流等現(xiàn)象.利用LEA對Br-運移進行模擬,并反演參數(shù)s、,由于僅有2個未知參數(shù),因此模型解的不唯一性較低,s取實驗中的飽和水流通量.由結(jié)果可知,擬合曲線與實測曲線基本重合,2均在0.98以上,RMSE均小于0.027,擬合效果較好.飽和含水量s與實測值0.46cm3/cm3相差不大,彌散度均在0.13cm左右.

    2.4 不同鐵含量時SDZ擬合曲線及參數(shù)確定

    采用Br-示蹤實驗擬合得到的s、作為固定參數(shù),利用OSM模型對不同鐵含量時SDZ的運移進行模擬,由于模型中的F、、解不是唯一確定的,所以取不同的參數(shù)初始值進行多次模擬計算,直到模擬值與實測值的誤差很小,此時的參數(shù)組合即是最佳參數(shù)組合;阻滯因子d根據(jù)公式(6)求得.由圖7、表2可知,不同鐵含量時SDZ運移的擬合值與實測值接近,2均大于0.97,RMSE均小于0.032,擬合效果較好.OSM模型假定吸附過程為動力學(xué)吸附,不存在瞬時吸附,所以瞬時吸附點位所占分數(shù)=0.吸附分配系數(shù)F隨著IOCS表面鐵含量的增加而增大,說明鐵含量越大,促進了SDZ的吸附.阻滯因子d與F呈相同的變化規(guī)律,說明隨著IOCS表面鐵含量的增大,SDZ遷移受到的阻力作用越大.吸附特征參數(shù)>1,表明不同鐵含量的IOCS表面具有較強的均一性[32].一階動力學(xué)速率常數(shù)數(shù)值較小,說明SDZ在介質(zhì)表面吸附性較弱.

    表2 不同鐵含量時SDZ穿透曲線擬合相關(guān)參數(shù)

    2.5 不同離子強度時SDZ擬合曲線及參數(shù)確定

    不同離子強度下SDZ擬合過程及參數(shù)確定依據(jù)與不同鐵含量情況下類似.由圖8、表3可知,不同離子強度下SDZ穿透曲線與擬合曲線基本吻合,且2均在0.97以上,RMSE均小于0.027,擬合效果較好.在QS和IOCS中,F均隨著離子強度的增大而減小,說明高離子強度促進了SDZ遷移,且同一離子強度下,相對于QS,IOCS中F值均較大.>1,說明QS和IOCS表面均一性較強.數(shù)值較小,說明此條件下吸附作用比較弱.阻滯因子d與F變化趨勢相同,說明離子強度越大,SDZ遷移過程中遇到的阻力作用越小,越容易遷移,這與實驗結(jié)果相符合.

    圖8 不同離子強度下SDZ的擬合曲線

    表3 不同離子強度下SDZ穿透曲線擬合相關(guān)參數(shù)

    續(xù)表3

    3 結(jié)論

    3.1 氧化鐵的存在會抑制SDZ遷移,且鐵含量越多,抑制作用越大.當(dāng)IOCS表面鐵含量為3.5mg/g時,SDZ出流的相對濃度峰值/0最小,約為0.71.

    3.2 離子強度的增加會促進SDZ在QS和IOCS中的遷移,當(dāng)CaCl2濃度為0.05mol/L時,SDZ出流的相對濃度峰值/0最大,分別在0.90和0.76左右.這可能是高離子強度下,Ca2+與SDZ競爭吸附作用較大造成的.

    3.3 運用Hydrus-1D中非平衡單點吸附模型對SDZ遷移過程進行模擬,擬合所得2均大于0.97,RMSE均小于0.032,擬合效果較好;反演得到的吸附分配系數(shù)F、阻滯因子d均隨著鐵含量的增加而增大,隨離子強度的增加而減小,>1,數(shù)值均較小.

    [1] Chen H, Ma L Q, Gao B, et al. Influence of Cu and Ca cations on ciprofloxacin transport in saturated porous media [J]. Journal of hazardous materials, 2013,262:805-811.

    [2] 葛林科,任紅蕾,魯建江,等.我國環(huán)境中新興污染物抗生素及其抗性基因的分布特征[J]. 環(huán)境化學(xué), 2015,34(5):875-883. Ge L K, Ren H L, Lu J J, et al. Occurrence of antibiotics and corresponding resistance genes in the environment of China [J]. Environmental Chemistry, 2015,34(5):875-883.

    [3] Sarmah A K, Meyer M T, Boxall A. A global perspective on the use, sales, exposure pathways, occurrence, fate and effects of veterinary antibiotics (VAs) in the environment [J]. Chemosphere, 2006,65(5): 725-759.

    [4] Diaz-Cruz M S, Maria J, Damià B. Environmental behavior and analysis of veterinary and human drugs in soils, sediments and sludge [J]. Trends in Analytical Chemistry, 2003,22(6):340-351.

    [5] 成玉婷,吳小蓮,向 壘,等.廣州市典型有機蔬菜基地土壤中磺胺類抗生素污染特征及風(fēng)險評價 [J]. 中國環(huán)境科學(xué), 2017,37(3):1154- 1161. Cheng Y T, Wu X L, Xiang L, et al. Distribution and risk assessment of sulfonamide antibiotics in soil from organic vegetable farms in Guangzhou [J]. China Environmental Science, 2017,37(3):1154-1161.

    [6] 張敬沙,龐炳坤,吳 杰,等.獸用抗生素磺胺二甲嘧啶對麥田NH3揮發(fā)的影響[J]. 中國環(huán)境科學(xué), 2018,38(6):2039-2046. Zhang J S, Pang B K, Wu J, et al. Effect of veterinary antibiotics on NH3volatilization from winter-wheat field [J]. China Environmental Science, 2019,(7):1-10.

    [7] 趙方凱,楊 磊,喬 敏,等.土壤中抗生素的環(huán)境行為及分布特征研究進展[J]. 土壤, 2017,49(3):428-436. Zhao F K, Yang L, Qiao M, et al. Occurrence of antibiotics and corresponding resistance genes in the environment of China [J]. Soils, 2017,49(3):428-436.

    [8] 郭學(xué)濤,苗 育,葛建華,等.pH值和離子強度對針鐵礦吸附磺胺二甲基嘧啶的影響機制[J]. 環(huán)境科學(xué)學(xué)報, 2016,36(7):2476-2482. Guo X T, Miao Y, Ge J H, et al. The influence of pH and ionic strength on the sorption of sulfamethazine by goethite [J]. Acta Scientiae Circumstantiae, 2016,36(7):2476-2482.

    [9] Martin S, Shchukarev A, Hanna K, et al. Kinetics and mechanisms of ciprofloxacin oxidation on hematite surfaces [J]. Environmental Science & Technology, 2015,49(20):12197-12205.

    [10] 顧 維,趙 玲,董元華,等.針鐵礦吸附諾氟沙星特征的研究[J]. 中國環(huán)境科學(xué), 2011,31(8):1314-1320. Gu W, Zhao L, Dong Y H, et al. Characteristics of norfloxacin sorption on goethite [J]. China Environmental Science, 2011,31(8):1314-1320.

    [11] Raquel A, Figueroa, Allison A, et al. Sorption of oxytetracycline to iron oxides and iron oxide·rich soils [J]. Environmental Science & Technology, 2005,39:6664-6671.

    [12] 丁春生,盧敬科,張德華,等.鐵鹽改性砂制備及其吸附Zn2+的性能研究[J]. 環(huán)境污染與防治, 2009,31(11):44-52. Ding C S, Lu J K, Zhang D H, et al. Preparztion of ferric salt modified sands for Zn2+adsorption [J]. Environmental Pollution Control, 2009, 31(11):44-52.

    [13] Deng J W, Chong Y S, Yan J, et al. Role of solution chemistry in the retention and release of graphene oxide nanomaterials in uncoated and iron oxide-coated sand [J]. Science of the Total Environment, 2017, 579:776-785.

    [14] Qi Z C, Du T T, Ma P K, et al. Transport of graphene oxide in saturated quartz sand containing iron oxides [J]. Science of the Total Environment, 2019,657:1450-1459.

    [15] Phili, P R, Johnson, Sun N, et al. Colloid transport in geochemically heterogeneous porous media: modeling and measurements [J]. Environmental Science & Technology, 1996,30:3284-3293.

    [16] Bunn R A, Magelky R D, Ryan J N, et al. Mobilization of natural colloids from an iron oxide coated sand aquifer: effect of pH and ionic strength [J]. Environmental Science & Technology, 2002,36:314-322.

    [17] Pierce M L, Moore C B. Adsorption of arsenite and arsenate on amorphous iron hydroxide [J]. Water Research, 1982,16(7):1247- 1253.

    [18] Hsu J C, Lin C J, Liao C H, et al. Removal of As(V) and As(III) by reclaimed iron-oxide coated sands [J]. Journal of Hazardous Materials, 2008,153(1/2):817-826.

    [19] 鄧慧萍,徐迪民,易小萍,等.幾種改性濾料去除水中有機物的性能比較[J]. 同濟大學(xué)學(xué)報, 2001,29(4):444-447. Deng H P, Xu D M, Yi X P, et al. Comparison of performance of several modified filter materials for removing organic matter from water [J]. Journal of Tongji University, 2001,29(4):444-447.

    [20] 趙玉華,賈 瑩,張 旭.不同改性石英砂去除有機物的效能比較[J]. 中國給水排水, 2007,23(19):15-18. Zhao Y H, Jia Y, Zhang X. Comparison of removal efficiency of organic matter by different modified quartz sands [J]. China Water & Wastewater, 2007,23(19):15-18.

    [21] 李 瑋,何江濤,劉麗雅,等.Hydrus-1D軟件在地下水污染風(fēng)險評價中的應(yīng)用[J]. 中國環(huán)境科學(xué), 2013,33(4):639-647.Li W, He J T, Liu L Y, et al. Application of Hydrus-1D software in groundwater contamination risk assessment [J]. China Environmental Science, 2013,33(4):639-647.

    [22] Brigante M, Zanini G, Avena M. Effect of humic acids on the adsorption of paraquat by goethite [J]. Journal of Hazardous Materials, 2010,184(1-3):241-247.

    [23] 尚爾寧,沈惠明,樊 祥.分光光度法測定磺咳合劑中磺胺嘧啶的含量[J]. 抗感染藥學(xué), 2004,1(2):84-85. Shang E N, Shen H M, Fan X. Determination of sulfadiazine content in sulfadiazine anticough mixture by spectrophotometry [J]. Anti-Infection Pharmacy, 2004,1(2):84-85.

    [24] Toride N F, Leij F J, Genuchten M T V. The CXTFIT code for estimating transport parameters from laboratory or field tracer experiments [J]. Research, 1995,137.

    [25] 張步迪,林 青,徐紹輝.磺胺嘧啶在原狀土柱中的運移特征及模擬[J]. 土壤學(xué)報, 2018,55(4):879-888. Zhang B D, Lin Q, Xu S H, et al. Characteristics and simulation of transport of sulfadiazine in undisturbed soil columns [J]. Acta Pedologica Sinica, 2018,55(4):879-888.

    [26] 徐秀姣,林 青,徐紹輝.環(huán)丙沙星在石英砂中的吸附遷移特征及參數(shù)分析[J]. 環(huán)境科學(xué)學(xué)報, 2016,36(6):2085-2094.Xu X J, Lin Q, Xu S H. Adsorption transport and parameter analysis of ciprofloxacin in quartz sand [J]. Acta Scientiae Circumstantiae, 2016, 36(6):2085-2094.

    [27] Qi Z C, Du T T, Ma P K, et al. Transport of graphene oxide in saturated quartz sand containing iron oxides [J]. Science of the Total Environment, 2019,657:1450-1459.

    [28] Se?bastien M, Shchukarev A, Hanna K, et al. Kineticsand mechanisms of ciprofloxacin oxidationon hematite surfaces [J]. Environmental Science & Technology, 2015,49(20):12197-12205.

    [29] 李學(xué)垣.土壤化學(xué)及實驗指導(dǎo)[M]. 北京:農(nóng)業(yè)出版社, 1997:80-81. Li X H. Soil chemistry and experimental guidance [M]. Beijing: Agriculture Press, 1997.80-81.

    [30] Murray J W. The surface chemistry of hydrous manganese dioxide [J]. Journal of Colloid and Interface Science, 1974,46(3):357-371.

    [31] 張勁強,董元華.陽離子強度和陽離子類型對諾氟沙星土壤吸附的影響 [J]. 環(huán)境科學(xué), 2007,28(10):2383-2388. Zhang J Q, Dong Y H. Influence of strength and species of cation on adsorption of norfloxacin in typical soils of China [J]. Environmental Science, 2007,28(10):2383-2388.

    [32] Selim H M. Modeling the transport and retention of inorganics in soils [J]. Advances in Agronomy, 1992,47:331-384.

    Characteristics of sulfadiazine migration in iron oxide coated sand and a comparison with simulation results.

    ZHANG Hui, LIN Qing, XU Shao-hui*

    (College of Environmental Science and Engineering, Qingdao University, Qingdao 266071, China)., 2019,39(11):4712~4721

    In order to investigate the migration characteristics of sulfonamide antibiotics in soil, sulfadiazine (SDZ) was used as a model component, and quartz sand (QS) was used as a medium to obtain iron oxide coated sand (IOCS). The effects of iron content and ionic strength on the migration of SDZ in IOCS were studied through column experiments, and the obtained results were simulated by Hydrus-1D software. The results of column experiments showed that the migration of SDZ was inhabited with an increase in the iron content on the surface of IOCS. The more the iron content, the higher the roughness, and the more the positive charge is carried through, which is beneficial to the adsorption of SDZ. An increase in the ionic strength promoted the migration of SDZ in QS and IOCS, which is mainly controlled by electrostatic interaction, surface complexation, and competitive adsorption. The migration process of SDZ under different iron contents and ionic strength can be better simulated by non-equilibrium one-site sorption model (OSM) in Hydrus-1D. The fitted curve is in agreement with the measured curve, where2>0.97 and RMSE<0.032. The adsorption distribution coefficient,Fincreased with an increase in iron content and decreased with an increase in ionic strength. The retardation factor,dshowed the same trend asF. The adsorption characteristic parameter>1, and the first-order kinetic rate constantwere observed to be small.

    sulfadiazine;iron oxide coated sand;migration;Hydrus-1D;iron content;ionic strength

    X131.1

    A

    1000-6923(2019)11-4712-10

    張 惠(1995-),女,山東聊城人,青島大學(xué)碩士研究生,主要從事地下環(huán)境中水流和溶質(zhì)運移及數(shù)值模擬.

    2019-04-24

    國家自然科學(xué)基金資助項目(41571214,41807010)

    * 責(zé)任作者, 教授, shhxu@qdu.edu.cn

    猜你喜歡
    流液石英砂離子
    高溫高壓下石英砂粒徑對油井水泥石性能的影響
    光學(xué)玻璃電熔窯流液洞的數(shù)值模擬分析
    硅酸鹽通報(2021年1期)2021-02-23 01:49:54
    鳴律
    炎黃地理(2021年12期)2021-01-04 18:52:33
    母豬子宮內(nèi)膜炎和產(chǎn)道惡露綜合征診治
    羅布麻莖傷流液成分研究
    在細節(jié)處生出智慧之花
    小議離子的檢驗與共存
    石英砂巖礦浮選降鐵試驗研究
    鋼渣對亞鐵離子和硫離子的吸附-解吸特性
    鋁離子電池未來展望
    videos熟女内射| 欧美日韩成人在线一区二区| 男女下面插进去视频免费观看 | 97在线人人人人妻| 免费少妇av软件| 国产精品麻豆人妻色哟哟久久| 性色av一级| 天堂俺去俺来也www色官网| 日韩电影二区| 亚洲成av片中文字幕在线观看 | 26uuu在线亚洲综合色| 成人18禁高潮啪啪吃奶动态图| 美女xxoo啪啪120秒动态图| 午夜激情av网站| 日本欧美视频一区| 精品福利永久在线观看| 久久人人爽人人爽人人片va| 亚洲色图综合在线观看| 韩国精品一区二区三区 | 中文字幕精品免费在线观看视频 | 精品久久国产蜜桃| 97精品久久久久久久久久精品| 日韩视频在线欧美| 精品亚洲乱码少妇综合久久| 99视频精品全部免费 在线| 国产熟女午夜一区二区三区| av在线播放精品| 免费看av在线观看网站| 国产毛片在线视频| 国产精品一国产av| 综合色丁香网| 亚洲精品国产av蜜桃| 亚洲久久久国产精品| 日产精品乱码卡一卡2卡三| 欧美少妇被猛烈插入视频| 中文乱码字字幕精品一区二区三区| 精品亚洲成国产av| 久久99一区二区三区| av又黄又爽大尺度在线免费看| 午夜91福利影院| 久久久久久人妻| h视频一区二区三区| 国产免费一级a男人的天堂| 黑人高潮一二区| 最近手机中文字幕大全| 久久狼人影院| 视频在线观看一区二区三区| 国产综合精华液| 国产精品一二三区在线看| 久久毛片免费看一区二区三区| av网站免费在线观看视频| 2018国产大陆天天弄谢| 99热全是精品| 精品少妇黑人巨大在线播放| 国产极品粉嫩免费观看在线| 午夜福利影视在线免费观看| 精品视频人人做人人爽| 日韩一本色道免费dvd| 亚洲av.av天堂| 免费播放大片免费观看视频在线观看| 久久久国产欧美日韩av| av在线老鸭窝| 成人综合一区亚洲| 免费在线观看完整版高清| 伦理电影大哥的女人| 国产av精品麻豆| 亚洲一区二区三区欧美精品| 国产一区二区三区av在线| 一级毛片 在线播放| 午夜久久久在线观看| 最近中文字幕高清免费大全6| 一本大道久久a久久精品| 人妻系列 视频| 老熟女久久久| 国产淫语在线视频| 99国产精品免费福利视频| 色网站视频免费| 只有这里有精品99| 中文字幕av电影在线播放| 三上悠亚av全集在线观看| 黄色一级大片看看| xxx大片免费视频| 在线 av 中文字幕| 午夜福利,免费看| 国产精品无大码| 亚洲精品久久午夜乱码| 新久久久久国产一级毛片| 精品国产露脸久久av麻豆| 色吧在线观看| 久久精品夜色国产| 成年人免费黄色播放视频| 最黄视频免费看| 在线观看免费高清a一片| 亚洲图色成人| 满18在线观看网站| 综合色丁香网| 丝瓜视频免费看黄片| 国产av码专区亚洲av| 亚洲成国产人片在线观看| 国内精品宾馆在线| 亚洲性久久影院| 成人黄色视频免费在线看| av国产精品久久久久影院| 久久久久久人人人人人| 99国产综合亚洲精品| 免费不卡的大黄色大毛片视频在线观看| 又黄又粗又硬又大视频| 日韩在线高清观看一区二区三区| 久久精品久久久久久噜噜老黄| 一二三四中文在线观看免费高清| 老女人水多毛片| 满18在线观看网站| 免费观看性生交大片5| 丝瓜视频免费看黄片| 香蕉丝袜av| 中文精品一卡2卡3卡4更新| 丝袜脚勾引网站| 18+在线观看网站| 如日韩欧美国产精品一区二区三区| 免费观看性生交大片5| 99久久人妻综合| 插逼视频在线观看| av一本久久久久| 丰满饥渴人妻一区二区三| 午夜91福利影院| 国产精品熟女久久久久浪| 亚洲图色成人| 精品福利永久在线观看| 国产成人av激情在线播放| 丝袜在线中文字幕| 成年美女黄网站色视频大全免费| 亚洲国产精品一区三区| av线在线观看网站| 日韩av不卡免费在线播放| 国产成人av激情在线播放| 高清毛片免费看| 大话2 男鬼变身卡| 久久97久久精品| 男女午夜视频在线观看 | 国产成人av激情在线播放| 天天躁夜夜躁狠狠躁躁| 22中文网久久字幕| 中文字幕av电影在线播放| 国产亚洲精品第一综合不卡 | 啦啦啦啦在线视频资源| 香蕉丝袜av| 夫妻午夜视频| 日韩制服丝袜自拍偷拍| 91午夜精品亚洲一区二区三区| 久久青草综合色| 国产乱人偷精品视频| 天天操日日干夜夜撸| 久久午夜综合久久蜜桃| 久久久久久久国产电影| 国产免费一级a男人的天堂| 制服人妻中文乱码| 高清毛片免费看| 久久99一区二区三区| 久久午夜综合久久蜜桃| 丝袜喷水一区| 午夜av观看不卡| 好男人视频免费观看在线| 美女中出高潮动态图| 亚洲激情五月婷婷啪啪| 精品少妇内射三级| 97超碰精品成人国产| 五月伊人婷婷丁香| 男人爽女人下面视频在线观看| 宅男免费午夜| 中国三级夫妇交换| 黄网站色视频无遮挡免费观看| 女性生殖器流出的白浆| 欧美另类一区| 亚洲av在线观看美女高潮| 国产女主播在线喷水免费视频网站| 2018国产大陆天天弄谢| 18禁动态无遮挡网站| 精品熟女少妇av免费看| 看免费av毛片| 久久精品熟女亚洲av麻豆精品| 如日韩欧美国产精品一区二区三区| 成人无遮挡网站| 久久精品熟女亚洲av麻豆精品| 国产精品一国产av| 9热在线视频观看99| 欧美日本中文国产一区发布| 各种免费的搞黄视频| 韩国av在线不卡| 久久精品人人爽人人爽视色| 2018国产大陆天天弄谢| 乱人伦中国视频| 国产女主播在线喷水免费视频网站| 色视频在线一区二区三区| av在线app专区| 男女边摸边吃奶| 中文字幕av电影在线播放| 国产xxxxx性猛交| av播播在线观看一区| 一级毛片我不卡| 另类亚洲欧美激情| 人人妻人人澡人人看| 婷婷色综合www| 欧美老熟妇乱子伦牲交| 日韩中文字幕视频在线看片| 欧美 日韩 精品 国产| 免费观看在线日韩| 日日爽夜夜爽网站| 日韩精品免费视频一区二区三区 | 国产精品国产三级国产专区5o| 国产免费现黄频在线看| 亚洲精品国产色婷婷电影| 乱人伦中国视频| 午夜日本视频在线| 久久国内精品自在自线图片| 女性生殖器流出的白浆| 人人澡人人妻人| 国产亚洲av片在线观看秒播厂| 国产黄色免费在线视频| 亚洲,欧美,日韩| 我的女老师完整版在线观看| 久久久久久人人人人人| 一级片免费观看大全| 欧美日韩视频高清一区二区三区二| videos熟女内射| 国产精品偷伦视频观看了| 又黄又粗又硬又大视频| 亚洲经典国产精华液单| 搡老乐熟女国产| 丁香六月天网| 国产一区二区在线观看日韩| 九色亚洲精品在线播放| 中文字幕人妻丝袜制服| 日本爱情动作片www.在线观看| 国产精品麻豆人妻色哟哟久久| 男人操女人黄网站| 下体分泌物呈黄色| videosex国产| 亚洲欧洲精品一区二区精品久久久 | 日韩人妻精品一区2区三区| 91精品国产国语对白视频| 久久这里有精品视频免费| 久久久久精品久久久久真实原创| 国产精品一二三区在线看| 日本-黄色视频高清免费观看| 天堂8中文在线网| 日本欧美视频一区| 国产国语露脸激情在线看| 国产精品99久久99久久久不卡 | 亚洲国产日韩一区二区| 成人影院久久| 日韩电影二区| 五月玫瑰六月丁香| 伊人久久国产一区二区| 久久精品国产自在天天线| 国产精品一区www在线观看| 夫妻性生交免费视频一级片| 一二三四中文在线观看免费高清| 宅男免费午夜| 欧美另类一区| 亚洲国产毛片av蜜桃av| 亚洲国产av影院在线观看| 久久99热6这里只有精品| 在线看a的网站| 国产1区2区3区精品| 80岁老熟妇乱子伦牲交| xxx大片免费视频| 自线自在国产av| 只有这里有精品99| 中文乱码字字幕精品一区二区三区| 久久精品人人爽人人爽视色| 亚洲精品av麻豆狂野| 亚洲精华国产精华液的使用体验| 亚洲精品色激情综合| 18禁观看日本| 18在线观看网站| 春色校园在线视频观看| 日韩,欧美,国产一区二区三区| 激情五月婷婷亚洲| 日韩电影二区| 又粗又硬又长又爽又黄的视频| 久久精品国产a三级三级三级| 香蕉国产在线看| 黄色视频在线播放观看不卡| 国产女主播在线喷水免费视频网站| 王馨瑶露胸无遮挡在线观看| 国产成人一区二区在线| 中文字幕最新亚洲高清| 午夜免费男女啪啪视频观看| 国产亚洲欧美精品永久| 国产亚洲av片在线观看秒播厂| av福利片在线| 免费看av在线观看网站| 国产精品成人在线| 韩国av在线不卡| 天天操日日干夜夜撸| 人妻人人澡人人爽人人| 大香蕉久久网| 亚洲国产精品一区三区| 日韩成人伦理影院| 亚洲精品国产色婷婷电影| 精品一区二区三卡| 我的女老师完整版在线观看| 久久久久久人妻| 男人操女人黄网站| 国产一级毛片在线| 波多野结衣一区麻豆| 麻豆精品久久久久久蜜桃| 中文天堂在线官网| 五月天丁香电影| 九九在线视频观看精品| 99九九在线精品视频| 大香蕉久久网| 美女福利国产在线| 日韩熟女老妇一区二区性免费视频| 日韩成人av中文字幕在线观看| 国产又色又爽无遮挡免| 激情五月婷婷亚洲| 国产不卡av网站在线观看| 又粗又硬又长又爽又黄的视频| 丁香六月天网| 少妇的丰满在线观看| 99久久人妻综合| 日韩伦理黄色片| 丝袜在线中文字幕| 香蕉国产在线看| 99热这里只有是精品在线观看| 一级片免费观看大全| 亚洲精品日本国产第一区| 插逼视频在线观看| 99久久精品国产国产毛片| 成人亚洲精品一区在线观看| 99久久综合免费| 精品少妇黑人巨大在线播放| 最近中文字幕高清免费大全6| 亚洲性久久影院| 国产高清国产精品国产三级| 中国美白少妇内射xxxbb| 国产亚洲午夜精品一区二区久久| 考比视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 制服人妻中文乱码| 国产伦理片在线播放av一区| 国产日韩一区二区三区精品不卡| 精品人妻在线不人妻| 一边亲一边摸免费视频| 久久精品熟女亚洲av麻豆精品| 欧美精品国产亚洲| 青春草视频在线免费观看| 国产男人的电影天堂91| 亚洲精品久久成人aⅴ小说| 久久久国产一区二区| 国产成人精品一,二区| 亚洲成国产人片在线观看| 在线观看人妻少妇| 亚洲精品乱久久久久久| 一级爰片在线观看| 久久久国产一区二区| 亚洲三级黄色毛片| 免费久久久久久久精品成人欧美视频 | 香蕉丝袜av| 一区二区三区乱码不卡18| 国产成人精品在线电影| 丰满少妇做爰视频| 成年美女黄网站色视频大全免费| 曰老女人黄片| 天天操日日干夜夜撸| 亚洲欧美中文字幕日韩二区| 视频在线观看一区二区三区| 久久精品久久久久久久性| 中文欧美无线码| 天堂俺去俺来也www色官网| 亚洲精品,欧美精品| 激情视频va一区二区三区| 蜜臀久久99精品久久宅男| 少妇 在线观看| 国产精品 国内视频| 久久久久国产精品人妻一区二区| 18在线观看网站| 人妻系列 视频| 久久99蜜桃精品久久| 超碰97精品在线观看| 亚洲熟女精品中文字幕| 狠狠精品人妻久久久久久综合| 日韩精品有码人妻一区| 97在线人人人人妻| 亚洲av福利一区| 精品国产一区二区三区四区第35| 国产精品免费大片| 热99久久久久精品小说推荐| 黄片播放在线免费| 性高湖久久久久久久久免费观看| 人妻人人澡人人爽人人| 国产女主播在线喷水免费视频网站| 亚洲欧美中文字幕日韩二区| 国产极品天堂在线| 如何舔出高潮| 18禁动态无遮挡网站| 精品人妻偷拍中文字幕| 在线免费观看不下载黄p国产| 蜜臀久久99精品久久宅男| 啦啦啦啦在线视频资源| 又大又黄又爽视频免费| 91成人精品电影| 国产精品人妻久久久影院| 久久精品熟女亚洲av麻豆精品| 97在线视频观看| kizo精华| 高清不卡的av网站| 老司机影院毛片| 我要看黄色一级片免费的| 免费女性裸体啪啪无遮挡网站| 亚洲精品日本国产第一区| 国产精品 国内视频| 欧美成人午夜精品| 久久精品国产自在天天线| 亚洲精品一区蜜桃| 国产精品国产三级国产av玫瑰| 精品国产一区二区三区四区第35| 成人影院久久| www.av在线官网国产| 男男h啪啪无遮挡| 国产有黄有色有爽视频| 9热在线视频观看99| 亚洲精品日韩在线中文字幕| 久久精品久久久久久久性| 美女主播在线视频| 丁香六月天网| 欧美少妇被猛烈插入视频| 亚洲国产精品999| 校园人妻丝袜中文字幕| 亚洲成色77777| 国产精品久久久久成人av| 国产成人精品无人区| 亚洲精品成人av观看孕妇| 两性夫妻黄色片 | 最新的欧美精品一区二区| 人妻系列 视频| 99香蕉大伊视频| 国产深夜福利视频在线观看| 国产成人精品婷婷| 尾随美女入室| 99久久综合免费| 午夜福利在线观看免费完整高清在| 午夜激情av网站| 国产精品成人在线| 国产在线一区二区三区精| 五月玫瑰六月丁香| 久久亚洲国产成人精品v| 黄色怎么调成土黄色| 九草在线视频观看| 国产在视频线精品| 国产精品 国内视频| 热re99久久精品国产66热6| 最新中文字幕久久久久| 国产日韩欧美在线精品| 成人漫画全彩无遮挡| 黄色一级大片看看| av不卡在线播放| 午夜日本视频在线| 成年av动漫网址| 国产欧美另类精品又又久久亚洲欧美| 免费看不卡的av| 亚洲精品久久午夜乱码| 久久国产精品大桥未久av| 男人爽女人下面视频在线观看| 久久99精品国语久久久| 久久精品国产亚洲av天美| 中文精品一卡2卡3卡4更新| 欧美激情 高清一区二区三区| 哪个播放器可以免费观看大片| 伊人久久国产一区二区| 97在线视频观看| 欧美日本中文国产一区发布| 我要看黄色一级片免费的| 亚洲精品自拍成人| 久久久久久久大尺度免费视频| 一级毛片 在线播放| 两个人免费观看高清视频| 国产综合精华液| 久久久久久久久久人人人人人人| 美女国产高潮福利片在线看| 中文欧美无线码| 伊人久久国产一区二区| 久久久精品94久久精品| 女性生殖器流出的白浆| 曰老女人黄片| 欧美日韩精品成人综合77777| 十八禁高潮呻吟视频| 日本猛色少妇xxxxx猛交久久| 精品少妇黑人巨大在线播放| 国产精品久久久久成人av| 汤姆久久久久久久影院中文字幕| 老熟女久久久| 男女高潮啪啪啪动态图| videosex国产| 亚洲欧美一区二区三区国产| 又大又黄又爽视频免费| 成年女人在线观看亚洲视频| 一区二区三区四区激情视频| 国产成人精品久久久久久| 少妇熟女欧美另类| 国产精品秋霞免费鲁丝片| 国产免费又黄又爽又色| 青青草视频在线视频观看| 人体艺术视频欧美日本| 熟女人妻精品中文字幕| 国产国语露脸激情在线看| 美女中出高潮动态图| 日本av手机在线免费观看| 中文字幕人妻丝袜制服| 在线观看www视频免费| 七月丁香在线播放| 热99久久久久精品小说推荐| 99久久综合免费| 久久免费观看电影| 韩国精品一区二区三区 | 看免费成人av毛片| 国产片特级美女逼逼视频| 一区在线观看完整版| 久久久国产欧美日韩av| 亚洲情色 制服丝袜| 香蕉丝袜av| 天堂俺去俺来也www色官网| 七月丁香在线播放| 中国美白少妇内射xxxbb| 欧美日韩国产mv在线观看视频| 久久久亚洲精品成人影院| 久久久久久久国产电影| 有码 亚洲区| 日韩,欧美,国产一区二区三区| 亚洲av免费高清在线观看| 日本av免费视频播放| 国产高清三级在线| 亚洲情色 制服丝袜| 中文字幕人妻丝袜制服| 国产精品久久久久久精品古装| 亚洲精品色激情综合| 99国产精品免费福利视频| 99久国产av精品国产电影| 国产精品无大码| 国产精品久久久av美女十八| 最近中文字幕高清免费大全6| 亚洲av日韩在线播放| 欧美日韩亚洲高清精品| 久久久国产一区二区| 久久久久国产网址| 国产片特级美女逼逼视频| 欧美bdsm另类| 亚洲精品久久成人aⅴ小说| 午夜福利视频精品| 亚洲精品视频女| 国产深夜福利视频在线观看| 欧美精品av麻豆av| 黑丝袜美女国产一区| 99热全是精品| 一本—道久久a久久精品蜜桃钙片| 欧美+日韩+精品| 女性被躁到高潮视频| 日韩伦理黄色片| 午夜视频国产福利| av国产久精品久网站免费入址| 超色免费av| 亚洲国产毛片av蜜桃av| 91在线精品国自产拍蜜月| 成人毛片a级毛片在线播放| 久久人人97超碰香蕉20202| 夫妻午夜视频| 99国产综合亚洲精品| 超碰97精品在线观看| 久久久国产一区二区| 久久久久久久久久久免费av| 久久精品夜色国产| 最近最新中文字幕免费大全7| 99国产精品免费福利视频| 亚洲国产色片| 亚洲国产av新网站| 国产伦理片在线播放av一区| 热99国产精品久久久久久7| 欧美日韩精品成人综合77777| 精品人妻在线不人妻| 丰满少妇做爰视频| 亚洲美女黄色视频免费看| 国产xxxxx性猛交| 国产成人a∨麻豆精品| 国产极品天堂在线| 在线观看一区二区三区激情| 国产成人免费无遮挡视频| 国产黄色免费在线视频| 欧美精品高潮呻吟av久久| 一区二区三区四区激情视频| 亚洲精品乱码久久久久久按摩| 亚洲国产精品国产精品| 色视频在线一区二区三区| 男的添女的下面高潮视频| 久久精品国产鲁丝片午夜精品| 青春草亚洲视频在线观看| 亚洲av日韩在线播放| 丝瓜视频免费看黄片| 国产精品女同一区二区软件| 午夜福利网站1000一区二区三区| 精品99又大又爽又粗少妇毛片| 久久ye,这里只有精品| 成年美女黄网站色视频大全免费| 18禁观看日本| 精品第一国产精品| 午夜福利网站1000一区二区三区| 亚洲伊人色综图| 免费大片18禁| 国产永久视频网站| 日本91视频免费播放| 最近中文字幕高清免费大全6| 一级,二级,三级黄色视频| 国产一区有黄有色的免费视频| 亚洲精品av麻豆狂野| av在线观看视频网站免费|