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

    典型PPCPs在河流沉積物中的吸附特性

    2019-04-28 12:52:36余綿梓李適宇胡嘉鏜
    中國環(huán)境科學(xué) 2019年4期
    關(guān)鍵詞:三氯氯霉素卡馬西平

    余綿梓,袁 嘯,李適宇,2*,胡嘉鏜,2**

    ?

    典型PPCPs在河流沉積物中的吸附特性

    余綿梓1,袁 嘯1,李適宇1,2*,胡嘉鏜1,2**

    (1.中山大學(xué)環(huán)境科學(xué)與工程學(xué)院,廣東 廣州 510275;2.廣東省環(huán)境污染控制與修復(fù)技術(shù)重點實驗室,廣東 廣州510275)

    選取了咖啡因、氯霉素、卡馬西平、磺胺甲噁唑和三氯生等5種PPCPs,在實驗室條件下,近似模擬了自然河流的水/沉積物界面,應(yīng)用中心復(fù)合實驗設(shè)計,考察了溫度、pH值、有機(jī)質(zhì)含量和流速對PPCPs吸附比例的影響;利用多元回歸方程擬合實驗數(shù)據(jù)得到PPCPs吸附模型,并利用獨立數(shù)據(jù)對模型進(jìn)行驗證,以建立適用于自然河流的PPCPs吸附模型.實驗結(jié)果表明,咖啡因和卡馬西平的吸附過程是放熱反應(yīng),而磺胺甲噁唑、氯霉素和三氯生的吸附為吸熱反應(yīng);pH值的升高對磺胺甲噁唑和三氯生的吸附能力產(chǎn)生抑制,但能促進(jìn)咖啡因吸附,而對氯霉素和卡馬西平的影響不大;有機(jī)質(zhì)含量、流速和初始濃度對5種PPCPs的吸附比例變化趨勢影響一致,隨著各因素濃度或速度升高,PPCPs的吸附比例均隨之增大,但影響程度有所不同.另外,吸附模型結(jié)果表明,5種PPCPs的吸附比例擬合值與實測值相關(guān)系數(shù)均達(dá)到0.8以上;獨立數(shù)據(jù)驗證結(jié)果表明,擬合值與實測值相關(guān)系數(shù)均達(dá)到0.8以上;因此,多元回歸方程能較好地擬合環(huán)境因素與PPCPs吸附之間的關(guān)系,同時,在所考察的環(huán)境因素濃度范圍內(nèi),多元回歸方程較好地預(yù)測了PPCPs在自然河流沉積物中的吸附行為.

    藥物和個人護(hù)理品(PPCPs);吸附;沉積物;影響因素;多元回歸方程

    藥物和個人護(hù)理品(PPCPs)的概念于1999年提出[1],由于使用量大、使用次數(shù)多,導(dǎo)致PPCPs持續(xù)不斷進(jìn)入河流、土壤和地下水等環(huán)境介質(zhì)中,尤其是在河流中的檢出率最高[2],近年來受到廣泛關(guān)注.相關(guān)研究表明,PPCPs在水環(huán)境中的暴露會對水生生物產(chǎn)生毒害作用[3],并通過食物鏈累積和人類食用行為進(jìn)入人體[4],對人類健康和生態(tài)系統(tǒng)構(gòu)成極大的威脅[5].

    PPCPs進(jìn)入水環(huán)境中會發(fā)生擴(kuò)散、降解和吸附等一系列遷移轉(zhuǎn)化過程,其中PPCPs在底層水和沉積物間的吸附過程是影響PPCPs最終歸趨的關(guān)鍵過程[6],直接影響到PPCPs在水體中的分布及其生物可利用度[7].當(dāng)沉積物受到水流沖刷時,水/沉積物界面的PPCPs則會通過擴(kuò)散、再懸浮等途徑,重新進(jìn)入上覆水體,引起水體二次污染[8].另外,PPCPs吸附到沉積物之后會發(fā)生垂向遷移過程,不僅會影響水流和溶質(zhì)的運移,還對深層沉積物和地下水造成危害[9].因此,了解 PPCPs在沉積物中的吸附行為和吸附特性是一項重要的研究內(nèi)容,對研究水環(huán)境中PPCPs的遷移轉(zhuǎn)化和歸趨有著重要意義.

    研究表明,pH值[10]、有機(jī)質(zhì)含量[11]和水體擾動程度[12]等因素對PPCPs的吸附結(jié)果均會產(chǎn)生影響.然而,目前對環(huán)境因素影響PPCPs吸附的研究多依據(jù)經(jīng)濟(jì)合作與發(fā)展組織化學(xué)品測試準(zhǔn)則OECD guideline 106批實驗方法,且通過單因素實驗進(jìn)行[13],該方法無法模擬真實河流中的水/沉積物界面,而且單因素實驗難以反映自然環(huán)境下多因素構(gòu)成的復(fù)雜體系,從而導(dǎo)致實驗研究與實際測定的結(jié)果往往出現(xiàn)偏差[14].另一方面,關(guān)于PPCPs吸附模型,目前大多采用吸附動力學(xué)和吸附等溫式模型,但上述模型均未考慮環(huán)境因素的作用[15-16],因此在實際應(yīng)用中擬合結(jié)果與實際測定結(jié)果存在較大誤差[17].為了考慮不同環(huán)境因素對PPCPs吸附的影響,多元回歸方程的應(yīng)用越來越得到重視.例如,Franco等[18]采用多元回歸方程較好地擬合了有機(jī)物吸附量與辛醇水分配系數(shù)及電離常數(shù)之間的關(guān)系;Al- Khazrajy等[19]同樣利用多元線性回歸擬合了5種PPCPs在不同土壤中的吸附,均取得較好的擬合效果.然而,上述研究中所采用的擬合數(shù)據(jù)僅限于所研究區(qū)域的沉積物,對于環(huán)境條件不同的其它區(qū)域,模型的應(yīng)用效果仍未可知.

    總體而言,目前對PPCPs在河流沉積物中的吸附行為研究尚有待深入,特別是模擬真實河流環(huán)境中多種因素影響下的吸附過程,另外,適用于描述自然水體的PPCPs吸附過程的數(shù)學(xué)模型亦有待進(jìn)一步建立.因此,本文選取卡馬西平、三氯生、氯霉素、磺胺甲噁唑和咖啡因等5種典型PPCPs(理化性質(zhì)如表1)為研究對象,在實驗室條件下近似模擬自然河流的水/沉積物界面,應(yīng)用中心復(fù)合實驗設(shè)計(CCD),考察pH值、溫度、有機(jī)質(zhì)含量和流速等4種環(huán)境因素對不同性質(zhì)PPCPs吸附的影響;在此基礎(chǔ)上,利用多元回歸方程對實驗數(shù)據(jù)進(jìn)行擬合,并結(jié)合獨立數(shù)據(jù)對擬合模型進(jìn)行預(yù)測和驗證,以建立適用于所考察因素濃度水平內(nèi)的PPCPs在河流沉積物中的吸附模型.

    表1 目標(biāo)物的理化性質(zhì)

    1 材料與方法

    1.1 沉積物樣品與實驗用品

    沉積物樣品:樣品采集自廣州市石井河的表層沉積物(0~5cm),采集的沉積物樣品經(jīng)風(fēng)干后過2mm 篩,在550℃下灼燒1h以去除有機(jī)質(zhì)及咖啡因,保存在4℃下待用.

    實驗藥品:咖啡因(CAF)純度99.5%、氯霉素(CAP)純度99.5%、卡馬西平(CBZ)純度99.5%、磺胺甲噁唑(SMX)純度99.5%和三氯生(TCS)純度99.5%均購自德國Dr. Ehrenstorfe公司.甲醇和乙腈(色譜純)購自德國Merck公司;腐植酸(HA)購自美國Aladdin公司.

    PPCPs儲備液:取一定量PPCPs標(biāo)準(zhǔn)品,用甲醇溶解配成濃度為1g/L的儲備液,于-18℃下保存待用.

    PPCPs標(biāo)準(zhǔn)使用液:取一定體積儲備液于50mL容量瓶中,用0.01mol/L CaCl2超純水定容,加入10mmol/L NaN3以抑制微生物作用.

    實驗儀器:三重四極桿液相色譜-質(zhì)譜聯(lián)用儀(LC-MS/MS,Thermo Fisher Scientific)、SPE固相萃取裝置(Automatic Science, China),pH計(Mettler Toledo, Swiss)、超純水儀(Millipore, Germany)、氮吹儀(Anpel Scientific, China).

    1.2 實驗方法

    吸附實驗在OECD Guideline 106[13]基礎(chǔ)上改進(jìn),以在實驗室條件下近似模擬自然河流的水/沉積物界面,裝置如圖1所示.實驗在一系列500mL燒杯中進(jìn)行,取20g純凈沉積物樣品平鋪于干凈燒杯中,加入20mL超純水在暗處浸潤24h使沉積物飽和[9].將上清液用針筒抽取去除.再用針筒緩慢加入100mL含100ng/L PPCPs溶液.溶液中含0.01mol/L CaCl2及100mg/L疊氮化鈉,CaCl2用以保持PPCPs的電離度,疊氮化鈉用以去除微生物影響;溶液中甲醇控制在0.1%以下以避免發(fā)生共溶現(xiàn)象.

    圖1 吸附實驗裝置

    實驗前,潔凈沉積物樣品在121℃下高壓滅菌15min 以去除微生物[24],并調(diào)節(jié)不同的pH值、溫度、有機(jī)質(zhì)含量、初始濃度和流速等條件,考察各因素對咖啡因在河流沉積物中吸附的影響.各因素的濃度水平根據(jù)響應(yīng)面分析軟件Design Expert 8.0.6中的中心復(fù)合實驗設(shè)計確定,取值范圍和大小如表2所示,得到實驗設(shè)計共86組.其中,pH值通過0.1mol/L HCl或0.1mol/L NaOH調(diào)節(jié).有機(jī)質(zhì)含量通過添加HA進(jìn)行調(diào)節(jié);由于自然河流沉積物中有機(jī)質(zhì)的成分復(fù)雜,在實驗室條件下難以模擬,而前人研究表明,底泥中天然有機(jī)質(zhì)的主要成分為HA[25-26],因此本實驗中有機(jī)質(zhì)含量的模擬通過在沉積物樣品中添加HA進(jìn)行,將一定量HA加入沉積物中,加入20mL超純水?dāng)嚢杈鶆蚝笤诎堤幗?4h,使HA與沉積物充分吸附同時使沉積物飽和[9].流速與溫度通過水浴恒溫振蕩器控制;流速大小的主要影響是改變水體的攪動狀況,由于在實驗室條件下難以模擬真實河流的水況,有學(xué)者將流速轉(zhuǎn)換為振蕩器轉(zhuǎn)速以模擬自然條件下不同的水體攪動狀況[12,27].本文采用王華等[28]通過大規(guī)模同步監(jiān)測及泥沙懸浮實驗建立的振蕩器轉(zhuǎn)速與水體流速之間的等量轉(zhuǎn)化關(guān)系,其公式為

    式中:為流速, m/s,為振蕩器轉(zhuǎn)速, r/min.

    實驗以不含底泥樣品作為對照以考察水解、壁吸附等因素的影響,樣品外裹錫紙以遮光去除光解作用.

    取樣后在5000r/min下離心10min,取上清液經(jīng)0.22μm濾膜過濾,濾液經(jīng)固相萃取和氮吹定容后進(jìn)行LC-MS/MS測定.咖啡因的吸附特性通過吸附比例K表示,可通過式(2)計算:

    式中:K為PPCPs的吸附比例;0為PPCPs的初始濃度, ng/L;C為取樣點小時實驗組上清液PPCPs濃度, ng/L;C為取樣時間點時對照組中咖啡因濃度, ng/L.

    表2 影響因素的大小及范圍

    1.3 PPCPs測定

    采用LC-MS/MS對實驗樣品的PPCPs濃度進(jìn)行測定.流動相為0.1%甲酸水溶液:乙腈=85:15,流速為0.3mL/min,進(jìn)樣量為2μL,咖啡因、卡馬西平、氯霉素、磺胺甲噁唑和三氯生的保留時間分別為1.55, 5.04,3.96,3.37和8.00min.

    2 結(jié)果與討論

    2.1 環(huán)境因素對PPCPs吸附的影響

    2.1.1 溫度對PPCPs吸附的影響 由圖2可知,隨著溫度的升高,咖啡因和卡馬西平的吸附比例出現(xiàn)下降,磺胺甲噁唑、氯霉素和三氯生的吸附比例則上升.在pH值為7.5、初始濃度為500ng/L、有機(jī)質(zhì)含量為2.5%、流速為0.5m/s、吸附時長為36h的條件下,隨著溫度從15℃升高到40℃,咖啡因的吸附比例下降了約10%,卡馬西平的吸附比例則下降了5%;而其他3種物質(zhì)的吸附比例均隨之增加約5%.

    從溫度對PPCPs吸附比例的影響可以看出,咖啡因和卡馬西平的吸附過程是放熱反應(yīng),吸附程度隨著溫度的升高而降低;這是由于溫度升高時,吸附平衡會向逆反應(yīng)方向移動,同時物質(zhì)的溶解度增大,從而抑制吸附反應(yīng)[29];此外,溫度升高還一定程度上促進(jìn)了解吸過程的發(fā)生,因此,溫度的升高使咖啡因和卡馬西平的吸附比例降低[30].在天然河流的溫度范圍內(nèi)(~20℃),溫度變化對咖啡因吸附的影響較小,當(dāng)溫度大于30℃并向上遞增時,咖啡因的吸附比例才出現(xiàn)較明顯的下降趨勢.由于溫度的升高同樣會促進(jìn)分子的擴(kuò)散,使得咖啡因與沉積物之間接觸增多,增加吸附反應(yīng)發(fā)生的機(jī)會,抵消一部分對吸附反應(yīng)的抑制作用,因此溫度較高時,升高溫度對咖啡因吸附的抑制作用相比溫度較低時更明顯.研究表明,大部分PPCPs的吸附反應(yīng)是放熱過程,吸附程度隨著溫度的升高而降低.如Wang等[31]研究表明卡馬西平在石墨上的吸附過程為放熱反應(yīng).本研究中磺胺甲噁唑、氯霉素和三氯生的吸附比例隨著溫度的升高而略有增大,表明其吸附過程為吸熱反應(yīng),也可能是由于分子間無規(guī)則運動即分子擺脫束縛引起的[32].王凱等[33]的研究也表明,磺胺甲噁唑的吸附是吸熱反應(yīng),吸附平衡量隨著溫度升高而升高,但變化幅度不大.

    圖2 溫度對PPCPs吸附的影響

    2.1.2 pH值對PPCPs吸附的影響 如圖3所示.磺胺甲噁唑和三氯生的吸附比例隨著pH值升高而降低,咖啡因的吸附比例在pH值約為10左右達(dá)到最高值,而氯霉素和卡馬西平則基本不受pH值變化的影響.

    pH值對PPCPs吸附的影響主要通過改變PPCPs的分子形態(tài),進(jìn)而影響其與沉積物顆粒的吸附機(jī)理從而影響吸附的程度.當(dāng)pH值低于物質(zhì)的p值時,酸性物質(zhì)大部分保持質(zhì)子態(tài)[34],比電離態(tài)更容易吸附,堿性物質(zhì)則帶有正電而易于吸附;隨著pH值升高,酸性物質(zhì)更加趨向電離帶負(fù)電,吸附能力減弱,而堿性物質(zhì)則水解趨向帶正電,疏水性有所增大,吸附能力增強(qiáng)對于兩性離子,pH值較低時,陽離子基團(tuán)占主導(dǎo),隨著pH值升高,去質(zhì)子化作用增強(qiáng),吸附能力下降.從實驗結(jié)果可以看出,與沉積物之間的靜電反應(yīng)是三氯生、磺胺甲噁唑和咖啡因的吸附機(jī)理之一,而卡馬西平和氯霉素則基本不通過靜電反應(yīng)進(jìn)行吸附.

    三氯生的分子結(jié)構(gòu)上含有一個—OH和多個Cl 基團(tuán),p= 7.8,呈弱酸性.在pH為7左右時,三氯生分子中83%~99%帶正電,因此吸附能力較強(qiáng)[21];當(dāng)pH>7.8時,三氯生主要以陰離子形式存在,因此表面帶有負(fù)電荷的沉積物對其吸附效果較差[33].磺胺甲噁唑同樣是一種酸性物質(zhì),在天然pH值條件下同時存在離子態(tài)和質(zhì)子態(tài)[35],隨著pH值變化,不同形態(tài)的磺胺甲噁唑所占比例也隨之變化.離子態(tài)分子由于帶有負(fù)電,與沉積物之間存在靜電排斥作用而吸附能力較弱,因此磺胺甲噁唑的吸附主要依靠質(zhì)子態(tài)分子的憎水反應(yīng)[36].

    咖啡因是一種堿性PPCPs,p值為10.4,因此在本研究pH值范圍內(nèi)(4~11),當(dāng)pH值小于p值時,咖啡因主要呈現(xiàn)電離態(tài),在水體中帶正電,而沉積物顆粒帶負(fù)電,因此咖啡因與沉積物顆粒之間可以通過靜電反應(yīng)進(jìn)行吸附[37].從結(jié)果可以看出,pH值的升高咖啡因趨向于帶正電,與帶負(fù)電荷的沉積物顆粒之間相互作用力增強(qiáng).但pH值大于pKa時,對咖啡因的水解產(chǎn)生抑制作用,從而降低吸附比例.

    卡馬西平和氯霉素的pKa值均較高(分別為13.9和11.03),為中性物質(zhì),憎水性吸附是主要的吸附機(jī)理,在本研究pH值范圍內(nèi),質(zhì)子態(tài)是卡馬西平和氯霉素的主要形態(tài),pH值對其影響很小[38].

    圖3 pH值對PPCPs吸附的影響

    圖4 有機(jī)質(zhì)含量對PPCPs吸附的影響

    2.1.3 有機(jī)質(zhì)含量對PPCPs吸附的影響 如圖4所示,有機(jī)質(zhì)含量增大,在不同程度上促進(jìn)PPCPs的吸附,其中影響最大的為卡馬西平,在吸附進(jìn)行13h未達(dá)到平衡時,有機(jī)質(zhì)含量從0%升高到5%,卡馬西平的吸附比例提高約40%,但隨著接近吸附平衡,有機(jī)質(zhì)含量變化的影響逐漸減弱.但對于其他4種PPCPs,有機(jī)質(zhì)含量變化的影響程度較卡馬西平小.

    圖5 流速對PPCPs吸附的影響

    圖6 初始濃度對PPCPs吸附的影響

    有機(jī)質(zhì)可以通過配體交換、陽離子橋接與陰離子交換等機(jī)制吸附至沉積物顆粒[39],同時一部分以溶解態(tài)形式溶于沉積物間隙水中.沉積物中的有機(jī)質(zhì)與PPCPs進(jìn)行結(jié)合是PPCPs吸附的途徑之一,吸附態(tài)的有機(jī)質(zhì)可以給PPCPs提供特異性吸附點位;此外,由于多數(shù)PPCPs帶有一定的疏水性,可通過疏水反應(yīng)與有機(jī)質(zhì)進(jìn)行吸附,然后通過有機(jī)質(zhì)與沉積物顆粒的吸附間接吸附在沉積物中[40-41].從有機(jī)質(zhì)含量對卡馬西平吸附的影響可以看出,卡馬西平和沉積物中有機(jī)質(zhì)之間的反應(yīng)是其主要的吸附機(jī)理,隨著有機(jī)質(zhì)含量的增大,能夠提供更多吸附點位與卡馬西平進(jìn)行憎水反應(yīng)[41],極性官能團(tuán)數(shù)量也隨之增多,可以作為π受體與卡馬西平進(jìn)行吸附[42].從有機(jī)質(zhì)含量對其他4種PPCPs吸附的影響可以看出,其他4種物質(zhì)不僅通過有機(jī)質(zhì)進(jìn)行吸附,還可通過范德華力、氫鍵作用等途徑與沉積物顆粒進(jìn)行吸附[10],但有機(jī)質(zhì)含量的增大仍有一定的促進(jìn)作用,表明與有機(jī)質(zhì)進(jìn)行吸附也是這幾種PPCPs進(jìn)行吸附的途徑.

    2.1.4 流速對PPCPs吸附的影響 如圖5所示,隨著流速的增大,5種PPCPs的吸附比例均隨之增大,并且增大幅度隨著反應(yīng)時間的增長而減小;除卡馬西平外,其余4種PPCPs吸附比例在57h后均達(dá)到最大值,而卡馬西平吸附比例在57h后仍隨著流速增大而有所增加.

    流速是自然河流重要的水動力指標(biāo),在不同時間、河流不同區(qū)域的流速均會存在差異,改變沉積物上覆水體的攪動狀況,水體的攪動會使沉積物顆粒再懸浮,增大顆粒與PPCPs之間的接觸面積,更有利于吸附過程的進(jìn)行[43].由于在實驗室條件下難以模擬真實河流的水況,因此關(guān)于流速對吸附的影響還鮮有研究.Hajj-Mohamad等[12]通過設(shè)置靜態(tài)和振蕩條件,近似模擬了水體的攪動和靜置狀態(tài),并對比了靜態(tài)和振蕩條件下4種藥物的吸附差異,發(fā)現(xiàn)振蕩條件下4種藥物的吸附能力均大于靜態(tài)條件;但該研究只是設(shè)置了2種條件,而沒有研究不同的攪動狀態(tài)對藥物吸附的影響.由于流速的模擬在實驗室條件下難以進(jìn)行,本研究將流速大小轉(zhuǎn)化為振蕩器轉(zhuǎn)速,有效地模擬了不同攪動水況對PPCPs吸附的影響.水體的攪動能促進(jìn)PPCPs的吸附,攪動幅度越大,PPCPs的吸附比例越高;但這種促進(jìn)作用隨著反應(yīng)時間的增長而減弱,這是由于隨著反應(yīng)時間增長,吸附反應(yīng)已經(jīng)接近平衡,且攪動情況下沉積物顆粒與水體接觸,也有利于解吸的進(jìn)行,因此可以從結(jié)果看到,達(dá)到吸附平衡時,流速的變化對最終吸附比例沒有影響.

    2.1.5 初始濃度對PPCPs吸附的影響 如圖6所示,在實驗條件范圍內(nèi),5種PPCPs的吸附比例均隨著初始濃度的增加而增大,增大幅度約為5%,但均呈現(xiàn)出隨著吸附時間增大而幅度降低的趨勢.

    前人研究多用吸附等溫式對不同濃度下PPCPs的吸附進(jìn)行模擬,如Wang等[31]研究表明卡馬西平和磺胺甲噁唑的吸附等溫式更符合Langmiur模型,磺胺甲噁唑則更符合Freundlich模型;Xu等[34]通過實驗分別對三氯生和磺胺甲噁唑的吸附等溫式進(jìn)行擬合,發(fā)現(xiàn)兩者的吸附等溫式均符合Freundlich公式,吸附能力隨著初始濃度的增大而增強(qiáng),但公式擬合結(jié)果中值不為1,說明2種物質(zhì)的吸附等溫式并不是線性的,不僅僅受到初始濃度的影響,還與其他環(huán)境因素有關(guān).

    本研究為模擬真實河流中PPCPs的濃度(ng/L級別),初始濃度的取值比一般實驗室研究的濃度(mg/L)較小,因此吸附點位對于PPCPs而言相對充足[22],雖然初始濃度的增大能促進(jìn)PPCPs的吸附,但其影響較小.

    2.2 多元回歸方程

    PPCPs吸附模型的構(gòu)建利用多元回歸方程進(jìn)行.采用多元逐步回歸方法[44],考察PPCPs吸附比例與各因素之間的關(guān)系,并應(yīng)用主成分分析方法篩選出各PPCPs的主要影響因素,以建立較為簡化的PPCPs吸附模型.本研究模型的構(gòu)建在SPSS 22.0中完成,將pH值、溫度、有機(jī)質(zhì)、初始濃度、流速和時間作為輸入量,吸附比例K作為輸出量.實驗數(shù)據(jù)共86組,隨機(jī)抽取其中9組作為獨立驗證數(shù)據(jù),另外77組作為擬合數(shù)據(jù);模型的可靠性通過擬合值與實測值的相關(guān)系數(shù)分析.

    得到簡化的PPCPs吸附模型如表3所示,式中參數(shù)乘積項表示2個因素對吸附比例的交互作用.從結(jié)果可以看出,擬合值與實測值的相關(guān)系數(shù)均大于0.8,相關(guān)性較強(qiáng),表明5種PPCPs吸附模型均可以較好地擬合出環(huán)境因素與吸附比例之間的關(guān)系.從主成分分析結(jié)果可以看出,5種PPCPs吸附模型的主要影響因素分別為:咖啡因(溫度、時間和流速),氯霉素(溫度、初始濃度和時間),卡馬西平(時間、有機(jī)質(zhì)含量和流速),磺胺甲噁唑(時間和流速)及三氯生(時間和流速).可以看出,時間因素對5種PPCPs的吸附均有較大影響;除氯霉素外,其他4種PPCPs吸附的主要影響因素均包含流速;此外,溫度、初始濃度和有機(jī)質(zhì)含量也為主要的因素,而pH值對5種PPCPs的吸附影響程度均較小.

    為檢驗?zāi)P偷膽?yīng)用效果,利用所建立的模型對獨立驗證數(shù)據(jù)進(jìn)行驗證,結(jié)果如圖7所示.從結(jié)果可以看出,5種PPCPs吸附模型對驗證數(shù)據(jù)的擬合值與實測值相關(guān)系數(shù)均達(dá)到0.8以上,相關(guān)性較強(qiáng),表明吸附模型可以很好地應(yīng)用于PPCPs吸附結(jié)果的預(yù)測.

    多元回歸方程近年來開始被應(yīng)用于吸附模型的建立.如Radka等[45]采用多元回歸方程擬合7種有機(jī)物的Freundlich系數(shù),擬合值與實測值的相關(guān)系數(shù)可達(dá)到0.97;He等[46]同樣采用多元回歸方程,建立五氯酚在不同土壤中的吸附模型,擬合值與實測值相關(guān)系數(shù)均達(dá)到0.8以上.但上述研究中,均未利用獨立數(shù)據(jù)對模型進(jìn)行驗證,因而吸附模型的應(yīng)用效果未知.本研究中,多元回歸方程能較好地擬合環(huán)境因素與PPCPs吸附之間的關(guān)系,且利用獨立數(shù)據(jù)進(jìn)行驗證,結(jié)果表明了多元回歸方程應(yīng)用于實際預(yù)測中的可行性.

    注::溫度(℃);:初始濃度(ng/L);:時間(h),:pH值;:有機(jī)質(zhì)含量(%);:流速(m/s).

    3 結(jié)論

    3.1 在實驗室條件下近似模擬了自然河流中的水/沉積物界面,綜合考察了不同環(huán)境因素對咖啡因、氯霉素、卡馬西平、磺胺甲噁唑和三氯生等5種PPCPs在沉積物中吸附的影響,并通過擬合和驗證多元回歸方程,得到適用于自然河流的PPCPs吸附模型.實驗結(jié)果表明,咖啡因和卡馬西平的吸附過程是放熱反應(yīng),而磺胺甲噁唑、氯霉素和三氯生的吸附為吸熱反應(yīng);pH值升高抑制磺胺甲噁唑和三氯生的吸附能力,而促進(jìn)咖啡因的吸附過程,但對氯霉素和卡馬西平的影響不大.隨著有機(jī)質(zhì)含量、流速和初始濃度值的增大,5種PPCPs的吸附比例均隨之增大,但增大程度有所不同.

    3.2 通過對77組實驗數(shù)據(jù)進(jìn)行多元回歸方程擬合,5種PPCPs的吸附比例擬合值與實測值相關(guān)系數(shù)均達(dá)到0.8以上.利用模型對9組獨立驗證數(shù)據(jù)進(jìn)行預(yù)測,擬合值與實測值相關(guān)系數(shù)也均達(dá)到0.8以上;結(jié)果表明,多元回歸方程能較好地擬合環(huán)境因素與PPCPs吸附之間的關(guān)系,驗證結(jié)果表明,在所考察的環(huán)境因素取值范圍內(nèi),將其應(yīng)用于預(yù)測自然河流中PPCPs在沉積物中的吸附是可行的.

    [1] Daughton C G, Ternes T A. Pharmaceuticals and personal care products in the environment: agents of subtle change?[J]. Environmental Health Perspectives, 1999,107(Suppl 6):907-938.

    [2] Dai G, Wang B, Huang J, et al. Occurrence and source apportionment of pharmaceuticals and personal care products in the Beiyun River of Beijing, China [J]. Chemosphere, 2015,119:1033-1039.

    [3] Foran C M, Bennett E R, Benson W H. Developmental evaluation of a potential non-steroidal estrogen: triclosan [J]. Marine Environment Research., 2000,50(1-5):153-156.

    [4] Lopes C, Persat H, Babut M. Transfer of PCBs from bottom sediment to freshwater river fish: a food-web modelling approach in the Rh?ne River (France) in support of sediment management. [J]. Ecotoxicology & Environmental Safety, 2012,81:17-26.

    [5] Kemper N. Veterinary antibiotics in the aquatic and terrestrial environment [J]. Ecological Indicators, 2008,8(1):1-13.

    [6] Zhou J, Broodbank N. Sediment-water interactions of pharmaceutical residues in the river environment [J]. Water Research, 2014,48:61-70.

    [7] Maeng S K, Sharma S K, Abel C D T, et al. Role of biodegradation in the removal of pharmaceutically active compounds with different bulk organic matter characteristics through managed aquifer recharge: Batch and column studies [J]. Water Research, 2011,45(16):4722- 4736.

    [8] Ying G, Zhao J, Zhou L, et al. Fate and occurrence of pharmaceuticals in the aquatic environment (surface water and sediment [J]. Comprehensive Analytical Chemistry, 2013,(62):453-557.

    [9] Fenet H, Mathieu O, Mahjoub O, et al. Carbamazepine, carbamazepine epoxide and dihydroxycarbamazepine sorption to soil and occurrence in a wastewater reuse site in Tunisia [J]. Chemosphere, 2012,88(1):49-54.

    [10] Karnjanapiboonwong A, Morse A N, Maul J D, et al. Sorption of estrogens, triclosan, and caffeine in a sandy loam and a silt loam soil [J]. Journal of Soils and Sediments, 2010,10(7):1300-1307.

    [11] Lin A Y, Lin C, Tung H, et al. Potential for biodegradation and sorption of acetaminophen, caffeine, propranolol and acebutolol in lab-scale aqueous environments [J]. Journal of Hazardous Materials, 2010,183(1-3):242-250.

    [12] Hajj-Mohamad M, Darwano H, Duy S V, et al. The distribution dynamics and desorption behaviour of mobile pharmaceuticals and caffeine to combined sewer sediments [J]. Water Research, 2017,108: 57-67.

    [13] Oecd. Adsorption - Desorption Using a Batch Equilibrium Method [M]. OECD Guideline for the Testing of Chemicals, Paris: Organization for economic Cooperation and Development, 2000:1-44.

    [14] Zhou J, Broodbank N. Sediment-water interactions of pharmaceutical residues in the river environment [J]. Water Res. 2014,48:61-70.

    [15] Ho Y S. Review of Second-Order Models for Adsorption Systems [J]. Journal of Hazardous Materials, 2006,136:681-689.

    [16] Samaraweera M, Jolin W, Vasudevan D, et al. Atomistic prediction of sorption free energies of cationic aromatic amines on montmorillonite: a linear interaction energy method [J]. Rapid Communications in Mass Spectrometry, 1993,7(7):163-166.

    [17] Nguyen T H, Goss K, Ball W P. Polyparameter linear free energy relationships for estimating the equilibrium partition of organic compounds between water and the natural organic matter in soils and sediments [J]. Environmental Science & Technology, 2005,39(4): 913-924.

    [18] Franco A, Trapp S. Estimation of the soil-water partition coefficient normalized to organic carbon for ionizable organic chemicals [J]. Environmental Toxicology & Chemistry, 2008,27(10):1995.

    [19] Al-Khazrajy O S A, Boxall A. Impacts of compound properties and sediment characteristics on the sorption behaviour of pharmaceuticals in aquatic systems [J]. Journal of Hazardous Materials, 2016,317:198- 209.

    [20] Zhang W, Ding Y, Boyd S A, et al. Sorption and desorption of carbamazepine from water by smectite clays [J]. Chemosphere, 2010, 81(7):954-960.

    [21] Revitt D M, Balogh T, Jones H. Sorption behaviours and transport potentials for selected pharmaceuticals and triclosan in two sterilised soils [J]. Journal of Soils and Sediments, 2015,15(3):594-606.

    [22] Martínez-Hernández V, Meffe R, Herrera S, et al. Sorption/desorption of non-hydrophobic and ionisable pharmaceutical and personal care products from reclaimed water onto/from a natural sediment [J]. Science of The Total Environment, 2014,472:273-281.

    [23] Liao P, Zhan Z, Dai J, et al. Adsorption of tetracycline and chloramphenicol in aqueous solutions by bamboo charcoal: a batch and fixed-bed column study [J]. Chemical Engineering Journal, 2013, 228:496-505.

    [24] 余綿梓,袁 嘯,李適宇,等.咖啡因在河流沉積物中吸附的影響因素及模擬研究[J]. 環(huán)境科學(xué)學(xué)報, 2018,38(2):560-569. Yu M Z, Yuan X, Li S Y, et al. Laboratory and simulation study on the adsorption of caffeine onto river sediments and the influencing factors [J]. Acta Scientiae Circumstantiae, 2018,38(2):560-569.

    [25] Nguyen L M. Organic matter composition, microbial biomass and microbial activity in gravel-bed constructed wetlands treating farm dairy wastewaters [J]. Ecological Engineering, 2000,16(2):199-221.

    [26] Stevenson F J. Humus chemistry: genesis, composition, reactions. [J]. Soil Science, 1994,135(2):129-130.

    [27] Gaullier C, Dousset S, Billet D, et al. Is pesticide sorption by constructed wetland sediments governed by water level and water dynamics? [J]. Environmental Science and Pollution Research, 2017, 25(4):1-12.

    [28] 王 華,逄 勇.藻類生長的水動力學(xué)因素影響與數(shù)值仿真[J]. 環(huán)境科學(xué), 2008,29(4):884-889. Wang H, Pang Y. Numerical simulation on hydrodynamic character for algae growth [J]. Environmental Science, 2008,29(4):884-889.

    [29] Zeng G, Zhang C, Huang G, et al. Adsorption behavior of bisphenol A on sediments in Xiangjiang River, Central-south China [J]. Chemosphere, 2006,65(9):1490-1499.

    [30] Cheng D, Liu X, Wang L, et al. Seasonal variation and sediment-water exchange of antibiotics in a shallower large lake in North China [J]. Science of The Total Environment, 2014,476-477:266-275.

    [31] Wang C, Li H, Liao S, et al. Coadsorption, desorption hysteresis and sorption thermodynamics of sulfamethoxazole and carbamazepine on graphene oxide and graphite [J]. Carbon, 2013,65:243-251.

    [32] 茶麗娟.卡馬西平在云南典型土壤上的吸附研究[D]. 昆明:昆明理工大學(xué), 2012. Cha M J. Sorption of carbamazepine in yunnan soils [D]. Kunming: Kunming University of Science and Technology, 2012.

    [33] 王 凱,李侃竹,周亦圓,等.河流沉積物對典型PPCPs的吸附特性及其影響因素[J]. 環(huán)境科學(xué), 2015,36(3):847-854. Wang K, Li K Z, Zhou Y Y, et al. Adsorption characteristics of typical PPCPs onto river dediments and its influencing factors [J]. Environmental Science, 2015,36(3):847-854.

    [34] Xu J, Wu L, Chang A C. Degradation and adsorption of selected pharmaceuticals and personal care products (PPCPs) in agricultural soils [J]. Chemosphere, 2009,77(10):1299-1305.

    [35] Hebig K H, Groza L G, Sabourin M J, et al. Transport behavior of the pharmaceutical compounds carbamazepine, sulfamethoxazole, gemfibrozil, ibuprofen, and naproxen, and the lifestyle drug caffeine, in saturated laboratory columns [J]. Science of The Total Environment, 2017,590-591:708-719.

    [36] Kay P, Blackwell P A, Boxall A B. Column studies to investigate the fate of veterinary antibiotics in clay soils following slurry application to agricultural land [J]. Chemosphere, 2005,60(4):497-507.

    [37] Fairbairn D J, Karpuzcu M E, Arnold W A, et al. Sediment-water distribution of contaminants of emerging concern in a mixed use watershed [J]. Science of the Total Environment, 2015,505:896-904.

    [38] Schaffer M, Boxberger N, B?rnick H, et al. Sorption influenced transport of ionizable pharmaceuticals onto a natural sandy aquifer sediment at different pH [J]. Chemosphere, 2012,87(5):513-520.

    [39] 王 磊,應(yīng)蓉蓉,石佳奇,等.土壤礦物對有機(jī)質(zhì)的吸附與固定機(jī)制研究進(jìn)展[J]. 土壤學(xué)報, 2017,54(4):805-818. Wang L, Ying R R, Shi J Q, et al. Advancement in study on adsorption of organic matter on soil minerals and its mechanism [J]. Acta Pedologica Sinica. 2017,54(4):805-818.

    [40] 鮑艷宇,周啟星,萬 瑩,等.土壤有機(jī)質(zhì)對土霉素在土壤中吸附-解吸的影響[J]. 中國環(huán)境科學(xué), 2009,29(6):651-655. Bao Y Y, Zhou Q X, Wang Y, et al. Effect of soil organic matter on adsorption and desorption of oxytetracycline in soils [J]. China Environmental Science, 2009,29(6):651-655.

    [41] Yamamoto H, Nakamura Y, Moriguchi S, et al. Persistence and partitioning of eight selected pharmaceuticals in the aquatic environment: Laboratory photolysis, biodegradation, and sorption experiments [J]. Water Research, 2009,43(2):351-362.

    [42] Ahmed M B, Zhou J L, Ngo H H, et al. Chloramphenicol interaction with functionalized biochar in water: sorptive mechanism, molecular imprinting effect and repeatable application [J]. Science of The Total Environment, 2017,609:885-895.

    [43] 王 鵬,馮 燕,蔡赟杰.水動力條件對太湖底泥吸附苯胺性能的影響[J]. 環(huán)境科學(xué)與技術(shù), 2013,36(5):117-123. Wang P, Feng Y, Cai Y J. Aniline adsorption properties on taihu lake sediment under different hydrodynamic conditions [J]. Environmental Science & Technology, 2013,36(5):117-123.

    [44] Arulsudar N, Subramanian N, Muthy R S. Comparison of artificial neural network and multiple linear regression in the optimization of formulation parameters of leuprolide acetate loaded liposomes [J]. Journal of Pharmacy Sciences, 2005,8(2):243-258.

    [45] Radke M, Lauwigi C, Heinkele G, et al. Fate of the antibiotic sulfamethoxazole and its two major human metabolites in a water sediment test [J]. Environmental Science & Technology, 2009,43(9): 3135-3141.

    [46] He Y, Xu J, Wang H, et al. Detailed sorption isotherms of pentachlorophenol on soils and its correlation with soil properties [J]. Environmental Research. 2006,101(3):362-372.

    致謝:中山大學(xué)環(huán)境科學(xué)與工程學(xué)院樊青娟老師對本研究中PPCPs的測定提供了幫助,劉廣州和徐闖同學(xué)在現(xiàn)場采樣過程中提供了幫助,在此一并表示感謝!

    Experimental and modeling study of sorption characteristics of selected PPCPs onto river sediments.

    YU Mian-zi1, YUAN Xiao1, LI Shi-yu1,2*, HU Jia-tang1,2**

    (1.School of Environmental Science and Engineering, Sun Yat-sun University, Guangzhou 510275, China;2.Guangdong Province Key Laboratory of Environmental Pollution Control and Remediation Technology, Guangzhou 510275, China)., 2019,39(4):1724~1733

    A water/sediment interface was roughly simulated in laboratory, and a central composite design was applied to investigate the influence of temperature, pH, organic content, and velocity on the sorption proprotion of five selected PPCPs (caffeine, chloramphenicol, carbamazepine, sulfamethoxazole and triclosan). Based on the experimental data, multiple regression equations were fitted and calibrated to establish favorable sorption models for PPCPs in the river systems. The results showed that the sorption process of caffeine and carbamazepine were exothermic reaction, and the sorption of sulfamethoxazole, chloramphenicol and triclosan were endothermic reaction. The increase of pH inhibited the sorption capacity of sulfamethoxazole and triclosan, but it could promote the sorption of caffeine, while it had little effect on the sorption of chloramphenicol and carbamazepine. The organic content, velocity and initial concentration had consistent effect on the sorption ratio of PPCPs. With the increase of the values of the factors, the sorption ratio of the five PPCPs increased with different degrees. The fitting and calibration results showed that the correlation coefficients between fitted and measured values of PPCPs sorption proportions were both above 0.8. Therefore, within the ranges of concentration of the investigated factors, the multiple regression equations were able to reasonably model and predict the sorption of PPCPs onto rivers sediment.

    PPCPs;sorption;sediment;parameters;multiple regression equation

    X522

    A

    1000-6923(2019)04-1724-10

    2018-08-03

    廣東省自然科學(xué)基金資助項目(2018A030313135);中央高?;究蒲袠I(yè)務(wù)費(20133800031650007);國家重點研發(fā)計劃項目(2016YFC0502803)

    *責(zé)任作者, 教授, eeslsy@mail.sysu.edu.cn; 副教授, hujtang@mail.sysu. edu.cn

    余綿梓(1993-),男,廣東潮州人,中山大學(xué)碩士研究生,主要從事PPCPs吸附模擬研究.發(fā)表論文1篇.

    猜你喜歡
    三氯氯霉素卡馬西平
    漲瘋了!碘漲50%,三氯漲超30%,溶劑漲超250%……消毒劑要漲價了
    一種氯霉素高靈敏消線法檢測試紙條的制備
    歐盟食品安全局:三氯蔗糖無致癌風(fēng)險
    三氯生對4種水生生物的急性毒性研究
    嬰幼兒慎用氯霉素眼藥水
    卡馬西平與加巴噴丁對復(fù)發(fā)性三叉神經(jīng)痛的治療價值比較
    HPLC法同時測定氯柳酊中氯霉素和水楊酸的含量
    卡馬西平與丙戊酸鈉緩釋片治療腦炎后癲癇的療效觀察
    氯霉素氧化鋅乳膏的制備及質(zhì)量標(biāo)準(zhǔn)
    添加有機(jī)酸加速2,4,6-三氯酚的生物降解
    肉色欧美久久久久久久蜜桃 | 男女边摸边吃奶| 97超视频在线观看视频| 免费观看无遮挡的男女| 色视频www国产| 国产探花在线观看一区二区| 一级毛片黄色毛片免费观看视频| 亚洲成人精品中文字幕电影| 国产男人的电影天堂91| 亚洲高清免费不卡视频| 18禁动态无遮挡网站| 国产伦理片在线播放av一区| 婷婷色麻豆天堂久久| 97超视频在线观看视频| 一本久久精品| 99热这里只有是精品在线观看| 久久精品熟女亚洲av麻豆精品 | 国产黄片视频在线免费观看| 久久精品久久久久久久性| 成人午夜精彩视频在线观看| 边亲边吃奶的免费视频| 中文精品一卡2卡3卡4更新| 国产在视频线精品| 色尼玛亚洲综合影院| 岛国毛片在线播放| av国产免费在线观看| 美女xxoo啪啪120秒动态图| 国产片特级美女逼逼视频| 男女那种视频在线观看| 高清午夜精品一区二区三区| 日韩制服骚丝袜av| 性插视频无遮挡在线免费观看| 女的被弄到高潮叫床怎么办| 亚洲最大成人中文| 91av网一区二区| 美女xxoo啪啪120秒动态图| 人妻系列 视频| 久久久久久久久中文| 一区二区三区高清视频在线| 国产成人精品婷婷| 国内精品宾馆在线| 日韩视频在线欧美| 三级国产精品欧美在线观看| 亚洲av福利一区| 毛片女人毛片| 99热这里只有是精品50| 夜夜看夜夜爽夜夜摸| 免费电影在线观看免费观看| 日韩欧美精品免费久久| 伊人久久国产一区二区| 三级男女做爰猛烈吃奶摸视频| 国产伦精品一区二区三区视频9| 在线观看免费高清a一片| 亚洲av电影不卡..在线观看| 大香蕉久久网| 热99在线观看视频| 性插视频无遮挡在线免费观看| 国产久久久一区二区三区| 少妇高潮的动态图| 最新中文字幕久久久久| 一个人观看的视频www高清免费观看| 中文字幕亚洲精品专区| 亚洲欧美日韩东京热| 免费观看a级毛片全部| 日本免费a在线| 啦啦啦啦在线视频资源| 搞女人的毛片| 少妇被粗大猛烈的视频| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久久久免费av| 日韩制服骚丝袜av| 春色校园在线视频观看| 亚洲人成网站高清观看| 大话2 男鬼变身卡| 2018国产大陆天天弄谢| 欧美xxxx性猛交bbbb| 中文字幕av在线有码专区| 男女边摸边吃奶| 搡老乐熟女国产| 国产精品一二三区在线看| 国产伦在线观看视频一区| av播播在线观看一区| 在线观看美女被高潮喷水网站| 婷婷色av中文字幕| 少妇高潮的动态图| 免费av毛片视频| 精品久久久精品久久久| 国产成人a区在线观看| 欧美xxⅹ黑人| 九色成人免费人妻av| 久99久视频精品免费| 小蜜桃在线观看免费完整版高清| 人人妻人人澡欧美一区二区| 最近中文字幕高清免费大全6| 如何舔出高潮| 国产高清不卡午夜福利| 亚洲自偷自拍三级| 国产精品美女特级片免费视频播放器| 国产成人午夜福利电影在线观看| 国产真实伦视频高清在线观看| 国产精品精品国产色婷婷| 国产69精品久久久久777片| av在线亚洲专区| 五月伊人婷婷丁香| 六月丁香七月| 观看免费一级毛片| 人妻夜夜爽99麻豆av| 国产亚洲午夜精品一区二区久久 | 国产亚洲精品av在线| 国产午夜精品久久久久久一区二区三区| 日韩成人av中文字幕在线观看| 天堂俺去俺来也www色官网 | 国产高清国产精品国产三级 | 日韩 亚洲 欧美在线| 久久久欧美国产精品| 自拍偷自拍亚洲精品老妇| 日韩成人av中文字幕在线观看| 午夜福利成人在线免费观看| 精品国产一区二区三区久久久樱花 | av黄色大香蕉| 国产免费一级a男人的天堂| 国产一区二区亚洲精品在线观看| 久久精品国产亚洲网站| 亚洲最大成人中文| 美女国产视频在线观看| 五月天丁香电影| 日韩一区二区视频免费看| 免费高清在线观看视频在线观看| 嫩草影院新地址| 亚洲美女搞黄在线观看| 亚洲第一区二区三区不卡| 久久久国产一区二区| 一级黄片播放器| 久久精品国产亚洲av涩爱| 日韩不卡一区二区三区视频在线| 毛片女人毛片| 亚洲精品一二三| 国产 一区精品| 亚洲精品影视一区二区三区av| 国产综合懂色| 久久热精品热| 69人妻影院| 又爽又黄a免费视频| 日本黄大片高清| 淫秽高清视频在线观看| 欧美 日韩 精品 国产| 三级男女做爰猛烈吃奶摸视频| 又粗又硬又长又爽又黄的视频| 亚洲精品成人av观看孕妇| 日韩国内少妇激情av| 国产成人91sexporn| 日韩亚洲欧美综合| 麻豆成人av视频| 热99在线观看视频| 成人av在线播放网站| av在线亚洲专区| 亚洲怡红院男人天堂| 好男人在线观看高清免费视频| 国产在线男女| 国产不卡一卡二| 亚洲四区av| 日韩在线高清观看一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 免费观看的影片在线观看| 可以在线观看毛片的网站| 黄片无遮挡物在线观看| 欧美xxxx黑人xx丫x性爽| 免费少妇av软件| 亚洲精品久久久久久婷婷小说| 两个人视频免费观看高清| 天堂影院成人在线观看| 亚洲精品亚洲一区二区| 国产精品久久久久久av不卡| 联通29元200g的流量卡| 国产 一区精品| 久久久久九九精品影院| 小蜜桃在线观看免费完整版高清| 免费播放大片免费观看视频在线观看| 成年人午夜在线观看视频 | 色哟哟·www| 国产精品伦人一区二区| 全区人妻精品视频| 最近视频中文字幕2019在线8| 熟妇人妻不卡中文字幕| av在线蜜桃| 综合色丁香网| 国产高清三级在线| 亚洲自拍偷在线| 久久久久九九精品影院| 色5月婷婷丁香| 亚洲伊人久久精品综合| 亚洲精品成人av观看孕妇| 国产黄片美女视频| 国产精品三级大全| 99热这里只有精品一区| 中文天堂在线官网| 一二三四中文在线观看免费高清| 精品国内亚洲2022精品成人| 欧美日韩精品成人综合77777| 免费观看的影片在线观看| 国产伦理片在线播放av一区| 免费不卡的大黄色大毛片视频在线观看 | 久久久久久久久久久免费av| 97在线视频观看| 欧美极品一区二区三区四区| 精品一区二区三区人妻视频| 欧美zozozo另类| 青春草视频在线免费观看| 91av网一区二区| 好男人在线观看高清免费视频| 日本熟妇午夜| 亚洲内射少妇av| 特大巨黑吊av在线直播| 亚洲精品aⅴ在线观看| 久久精品久久精品一区二区三区| 欧美一级a爱片免费观看看| 国产精品福利在线免费观看| 一本一本综合久久| 国产精品一区二区三区四区免费观看| 久久久精品94久久精品| 精品亚洲乱码少妇综合久久| 亚洲精品视频女| 大片免费播放器 马上看| 99久久精品一区二区三区| 日韩不卡一区二区三区视频在线| 三级毛片av免费| 国产精品1区2区在线观看.| 色综合站精品国产| 高清在线视频一区二区三区| 亚洲自拍偷在线| 五月天丁香电影| 看免费成人av毛片| 少妇的逼好多水| 99热这里只有精品一区| 婷婷色麻豆天堂久久| 亚洲精品国产av成人精品| 亚洲,欧美,日韩| 国产精品人妻久久久影院| 亚洲精品456在线播放app| 永久网站在线| 免费观看无遮挡的男女| 国产精品人妻久久久久久| 丰满人妻一区二区三区视频av| 熟女人妻精品中文字幕| 国产老妇女一区| 亚洲av一区综合| 国产成人精品一,二区| 国产精品麻豆人妻色哟哟久久 | 精品人妻一区二区三区麻豆| 亚洲国产精品sss在线观看| 久久久色成人| 波野结衣二区三区在线| 一级毛片久久久久久久久女| 精品一区二区三卡| 搡女人真爽免费视频火全软件| 久久久久性生活片| 亚洲第一区二区三区不卡| 日日摸夜夜添夜夜爱| 午夜免费男女啪啪视频观看| 亚洲一级一片aⅴ在线观看| 亚洲一区高清亚洲精品| 亚洲精品aⅴ在线观看| 性插视频无遮挡在线免费观看| 亚洲国产日韩欧美精品在线观看| 国产精品一区二区三区四区久久| 精品亚洲乱码少妇综合久久| 3wmmmm亚洲av在线观看| 麻豆精品久久久久久蜜桃| 国产免费福利视频在线观看| av卡一久久| 亚洲国产日韩欧美精品在线观看| 国产黄a三级三级三级人| 尤物成人国产欧美一区二区三区| 国产 一区精品| 51国产日韩欧美| 国产精品久久久久久久电影| 免费观看无遮挡的男女| 国产伦在线观看视频一区| 日韩av免费高清视频| 精品亚洲乱码少妇综合久久| 日韩伦理黄色片| 国产视频内射| 日韩av不卡免费在线播放| 插阴视频在线观看视频| 国产毛片a区久久久久| 乱人视频在线观看| 秋霞伦理黄片| 免费不卡的大黄色大毛片视频在线观看 | 天天躁夜夜躁狠狠久久av| 国产在视频线精品| 黄色日韩在线| 最近中文字幕高清免费大全6| 久久久久久久久中文| 777米奇影视久久| 精品亚洲乱码少妇综合久久| 男女下面进入的视频免费午夜| 日本欧美国产在线视频| 久久亚洲国产成人精品v| 国产永久视频网站| 亚洲第一区二区三区不卡| 亚洲自偷自拍三级| 亚洲精品日本国产第一区| 在线播放无遮挡| 亚洲欧美成人综合另类久久久| 国产精品熟女久久久久浪| 国产伦在线观看视频一区| 日日撸夜夜添| 乱系列少妇在线播放| 国产黄片美女视频| 久久久久久久午夜电影| 欧美极品一区二区三区四区| 91午夜精品亚洲一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 亚洲人成网站在线观看播放| 一级毛片 在线播放| 亚洲国产精品成人久久小说| 午夜免费男女啪啪视频观看| 国产 亚洲一区二区三区 | av在线播放精品| 老师上课跳d突然被开到最大视频| 日日啪夜夜撸| 好男人视频免费观看在线| 国产成人午夜福利电影在线观看| 亚洲av成人精品一二三区| 非洲黑人性xxxx精品又粗又长| 日韩成人伦理影院| 热99在线观看视频| 国产精品美女特级片免费视频播放器| 人人妻人人澡人人爽人人夜夜 | 亚洲精品第二区| 精品酒店卫生间| 午夜激情久久久久久久| 欧美极品一区二区三区四区| 日韩欧美 国产精品| 免费看美女性在线毛片视频| 亚洲久久久久久中文字幕| 寂寞人妻少妇视频99o| 久久99热这里只频精品6学生| 精品人妻一区二区三区麻豆| 男的添女的下面高潮视频| 最近中文字幕高清免费大全6| 老师上课跳d突然被开到最大视频| 黄片无遮挡物在线观看| 午夜免费激情av| 内射极品少妇av片p| 精品午夜福利在线看| 久热久热在线精品观看| ponron亚洲| 国产成人a∨麻豆精品| 啦啦啦韩国在线观看视频| 亚洲美女视频黄频| 能在线免费观看的黄片| 日韩国内少妇激情av| 国产一区亚洲一区在线观看| 极品教师在线视频| 两个人的视频大全免费| 久久韩国三级中文字幕| 欧美变态另类bdsm刘玥| 精品一区二区免费观看| 99热6这里只有精品| av.在线天堂| 国精品久久久久久国模美| 久久这里有精品视频免费| 丰满人妻一区二区三区视频av| 女的被弄到高潮叫床怎么办| 成人毛片60女人毛片免费| 69av精品久久久久久| 中文字幕av在线有码专区| 又粗又硬又长又爽又黄的视频| 99久久人妻综合| 精品久久久久久久久亚洲| 国产精品久久久久久精品电影小说 | 嫩草影院精品99| 精品一区二区三卡| 日本wwww免费看| 麻豆精品久久久久久蜜桃| 亚洲国产高清在线一区二区三| 亚洲aⅴ乱码一区二区在线播放| 欧美日本视频| 国产 亚洲一区二区三区 | 国产乱来视频区| 欧美三级亚洲精品| 国产免费视频播放在线视频 | 校园人妻丝袜中文字幕| 国产伦精品一区二区三区四那| 免费播放大片免费观看视频在线观看| 日韩亚洲欧美综合| 少妇高潮的动态图| a级毛片免费高清观看在线播放| 乱码一卡2卡4卡精品| 日韩伦理黄色片| 女人被狂操c到高潮| 日韩av在线免费看完整版不卡| 嫩草影院新地址| 中文字幕久久专区| 久久久久久久久大av| 秋霞在线观看毛片| 亚洲欧美成人精品一区二区| 国产一区二区在线观看日韩| 免费在线观看成人毛片| 亚洲第一区二区三区不卡| 久久韩国三级中文字幕| 欧美变态另类bdsm刘玥| 国产欧美另类精品又又久久亚洲欧美| 精品亚洲乱码少妇综合久久| 精华霜和精华液先用哪个| 久久精品国产亚洲网站| 国产成人精品婷婷| 最近最新中文字幕大全电影3| 国产精品av视频在线免费观看| 日韩一区二区三区影片| 在线观看av片永久免费下载| 99久国产av精品国产电影| 久久久久久久久久人人人人人人| 激情 狠狠 欧美| 干丝袜人妻中文字幕| 高清午夜精品一区二区三区| 激情 狠狠 欧美| 内地一区二区视频在线| 亚洲国产高清在线一区二区三| 激情五月婷婷亚洲| 国产v大片淫在线免费观看| 成年版毛片免费区| 亚洲国产欧美在线一区| 狂野欧美白嫩少妇大欣赏| 精品久久久久久久人妻蜜臀av| 久久久午夜欧美精品| 校园人妻丝袜中文字幕| 成年av动漫网址| 成人毛片a级毛片在线播放| av一本久久久久| 成年免费大片在线观看| 亚洲色图av天堂| 亚洲精品乱久久久久久| 久久精品国产亚洲网站| 国产激情偷乱视频一区二区| 91av网一区二区| 春色校园在线视频观看| 久久久久国产网址| 中国国产av一级| 中文字幕久久专区| 97超视频在线观看视频| 亚洲国产精品国产精品| 亚洲精品日韩av片在线观看| 麻豆精品久久久久久蜜桃| 在线观看av片永久免费下载| 搡老乐熟女国产| 成年免费大片在线观看| 26uuu在线亚洲综合色| 国产 一区 欧美 日韩| 精品亚洲乱码少妇综合久久| 亚洲欧美日韩无卡精品| 亚洲四区av| 免费av不卡在线播放| 亚洲国产欧美人成| av在线老鸭窝| 搡女人真爽免费视频火全软件| 一级av片app| 成人毛片a级毛片在线播放| 禁无遮挡网站| 色网站视频免费| 一区二区三区四区激情视频| 97超碰精品成人国产| 国产国拍精品亚洲av在线观看| 天堂影院成人在线观看| 国产不卡一卡二| 男人狂女人下面高潮的视频| 久久久久久九九精品二区国产| 一个人观看的视频www高清免费观看| 亚洲,欧美,日韩| 噜噜噜噜噜久久久久久91| 波多野结衣巨乳人妻| 国产精品三级大全| 久久精品熟女亚洲av麻豆精品 | av免费在线看不卡| 午夜久久久久精精品| 国产黄a三级三级三级人| 麻豆成人av视频| 欧美xxxx黑人xx丫x性爽| 身体一侧抽搐| 成人二区视频| 免费不卡的大黄色大毛片视频在线观看 | 日日啪夜夜爽| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 在线播放无遮挡| 国产免费视频播放在线视频 | 美女高潮的动态| 亚洲av成人精品一二三区| 免费少妇av软件| 欧美一区二区亚洲| 亚洲成人av在线免费| 熟女人妻精品中文字幕| 日韩伦理黄色片| 亚洲在线观看片| 深爱激情五月婷婷| 九九在线视频观看精品| 国产精品国产三级国产专区5o| 久久精品夜夜夜夜夜久久蜜豆| 黄色欧美视频在线观看| 超碰av人人做人人爽久久| 亚洲成人av在线免费| 国产麻豆成人av免费视频| 婷婷色综合大香蕉| 99热这里只有是精品50| 在线观看美女被高潮喷水网站| 三级男女做爰猛烈吃奶摸视频| 久99久视频精品免费| 蜜臀久久99精品久久宅男| 99久久人妻综合| 九色成人免费人妻av| 日韩,欧美,国产一区二区三区| 午夜福利视频1000在线观看| 国产黄频视频在线观看| 午夜久久久久精精品| or卡值多少钱| 最近中文字幕2019免费版| 久久精品夜色国产| 日本色播在线视频| 久久久久久久亚洲中文字幕| av在线天堂中文字幕| 最近中文字幕2019免费版| 国产中年淑女户外野战色| 成人国产麻豆网| 成人欧美大片| 亚洲aⅴ乱码一区二区在线播放| av在线蜜桃| 欧美成人精品欧美一级黄| 777米奇影视久久| 欧美+日韩+精品| 国语对白做爰xxxⅹ性视频网站| 91午夜精品亚洲一区二区三区| 九九爱精品视频在线观看| 日本猛色少妇xxxxx猛交久久| 日本熟妇午夜| 国产精品久久久久久精品电影小说 | 99re6热这里在线精品视频| 99久久精品国产国产毛片| av播播在线观看一区| 国产亚洲av片在线观看秒播厂 | 成人一区二区视频在线观看| 免费观看性生交大片5| av在线蜜桃| 18禁动态无遮挡网站| 日韩av免费高清视频| 国产成人精品福利久久| 国产精品美女特级片免费视频播放器| 亚洲欧洲国产日韩| 国产精品人妻久久久久久| 午夜老司机福利剧场| 婷婷色麻豆天堂久久| 中文字幕人妻熟人妻熟丝袜美| 一级二级三级毛片免费看| 熟女电影av网| 一级av片app| 国产男人的电影天堂91| 大香蕉久久网| 少妇人妻精品综合一区二区| 免费av毛片视频| 91精品一卡2卡3卡4卡| 成人无遮挡网站| 久久99蜜桃精品久久| 久久久久久久亚洲中文字幕| 日日干狠狠操夜夜爽| 久久久久久九九精品二区国产| 丝瓜视频免费看黄片| a级毛片免费高清观看在线播放| 三级男女做爰猛烈吃奶摸视频| 麻豆精品久久久久久蜜桃| 爱豆传媒免费全集在线观看| 特级一级黄色大片| 国产精品久久久久久精品电影小说 | 岛国毛片在线播放| 最近2019中文字幕mv第一页| 久久久久久久久大av| 亚洲,欧美,日韩| 国产淫片久久久久久久久| 精品国产三级普通话版| 亚洲人成网站在线观看播放| 插阴视频在线观看视频| 日本色播在线视频| 久久久精品免费免费高清| 成人无遮挡网站| 日韩欧美精品v在线| 日韩电影二区| 亚洲精品影视一区二区三区av| 亚洲精品,欧美精品| 亚洲在线自拍视频| 国精品久久久久久国模美| 日韩av在线免费看完整版不卡| 国产 一区 欧美 日韩| 免费观看性生交大片5| 国产乱人视频| 国产 一区 欧美 日韩| av在线播放精品| 亚洲av中文av极速乱| 亚洲人成网站高清观看| 亚洲aⅴ乱码一区二区在线播放| 国产黄色免费在线视频| 观看免费一级毛片| 亚洲欧美日韩卡通动漫| 2022亚洲国产成人精品| 亚洲人成网站高清观看| 丝袜美腿在线中文| 汤姆久久久久久久影院中文字幕 | 久久99热这里只有精品18| 日日摸夜夜添夜夜添av毛片| 国产精品嫩草影院av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲无线观看免费| 亚洲国产精品成人久久小说| 美女国产视频在线观看| 少妇熟女aⅴ在线视频| 亚洲精品成人av观看孕妇| 国产精品一区二区三区四区久久| 三级经典国产精品|