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

    拉薩河全氟化合物的時(shí)空分布特征研究

    2018-11-28 03:46:50孫殿超王小萍王傳飛
    中國環(huán)境科學(xué) 2018年11期
    關(guān)鍵詞:氟化合物拉薩河全氟

    孫殿超,龔 平,王小萍,3*,王傳飛

    ?

    拉薩河全氟化合物的時(shí)空分布特征研究

    孫殿超1,3,龔 平1,2,王小萍1,2,3*,王傳飛1,2

    (1.中國科學(xué)院青藏高原研究所,中國科學(xué)院環(huán)境與地表過程實(shí)驗(yàn)室,北京 100101;2.中國科學(xué)院青藏高原地球科學(xué)卓越創(chuàng)新中心,北京 100101;3.中國科學(xué)院大學(xué),北京 100049)

    對(duì)拉薩河27個(gè)采樣點(diǎn)、拉薩污水處理廠分別在豐水期和枯水期進(jìn)行水樣的采集,并用雙三元二維液相色譜-串聯(lián)質(zhì)譜法(DGLC-MS /MS)首次對(duì)水樣中的17種全氟化合物(PFASs)進(jìn)行檢測(cè)分析,結(jié)果表明,PFASs在拉薩普遍存在,拉薩河表層水中平均Σ13PFASs濃度為322pg/L,范圍從60到1724pg/L.拉薩河PFASs的濃度水平低于其他大城市流域,但高于南北極表層水體.拉薩河中的主要化合物為全氟丁酸(PFBA)(25%)、全氟戊酸(PFPeA)(22%)、全氟辛酸(PFOA)(14%)、全氟辛烷磺酸(PFOS)(14%)和全氟丁烷磺酸(PFBS) (13%).空間分布特征表明,拉薩河中全氟化合物的排放還存在一些點(diǎn)排放源,拉薩市產(chǎn)生的污水并未完全匯入到拉薩市污水處理廠中處理.基于PFBA和PFBS的全氟產(chǎn)品使用比較廣泛,而基于PFOA和PFOS的全氟產(chǎn)品只在某些區(qū)域或行業(yè)有相對(duì)多的使用量或排量.PFASs在豐水期和枯水期的濃度水平和組成特征有顯著的差異,拉薩全氟化合物產(chǎn)品的生產(chǎn)和使用模式受季節(jié)影響較大.污水處理廠是PFASs排放到拉薩河中的一個(gè)重要點(diǎn)源.污水處理廠每年排放到拉薩河中的全氟化合物約為37mg,拉薩河排放到雅魯藏布江的物質(zhì)荷載約為292t/a.在未來,所有的廢水都應(yīng)嚴(yán)格地納入污水處理廠經(jīng)達(dá)標(biāo)處理后再排放.盡管PFOA和PFOS在拉薩的生產(chǎn)和使用進(jìn)行了一定的轉(zhuǎn)型,但是基于全氟化合物的產(chǎn)品的限制應(yīng)進(jìn)一步系統(tǒng)地執(zhí)行和強(qiáng)有力地落實(shí).

    全氟烷基酸;持久性有機(jī)污染物;空間分布;季節(jié)變異;物質(zhì)荷載

    全氟化合物(PFASs)由于具有持久性、生物富集傾向和不斷出現(xiàn)的證據(jù)表明的生物有害性,因而日益?zhèn)涫荜P(guān)注.全氟化合物被廣泛應(yīng)用于織物、地毯、紙張、航空液壓油、泡沫滅火劑、不粘鍋、包裝材料等生產(chǎn)制造領(lǐng)域,因而其被大量地排放到環(huán)境中[1].根據(jù)全氟化合物化學(xué)鍵的類型,其分為離子型和分子型,其中,離子型的全氟化合物具有水溶性,故水環(huán)境是全氟化合物存在的一個(gè)重要載體.水又是生命體賴以生存的一種必需物質(zhì),從而全氟化合物能夠隨著水侵入到生物體和人體中,因此,全氟化合物在水環(huán)境中的遷移、轉(zhuǎn)化和歸趨需要引起足夠的重視.

    2009年5月,全氟辛烷磺酸(PFOS)和其鹽類以及全氟磺酰氟(PFOSF)被列入《斯德哥爾摩公約》持久性有機(jī)污染物名錄中.2015年10月,持久性有機(jī)污染物審查委員會(huì)提議將全氟辛酸(PFOA)列入持久性有機(jī)污染物審查名單.至此,全氟化合物的消減拉開了帷幕.全氟化合物的生產(chǎn)和使用也由長鏈(碳原子數(shù)大于8)向短鏈的化合物進(jìn)行轉(zhuǎn)變,如全氟丁酸(PFBA),全氟丁烷磺酸(PFBS).全氟化合物的生產(chǎn)中心也由歐美等發(fā)達(dá)國家向發(fā)展中國家轉(zhuǎn)移,尤其是亞洲國家[2].自此,全氟化合物的生產(chǎn)和排放模式均發(fā)生了變化.

    河流是全氟化合物生物地球化學(xué)循環(huán)的重要途徑.一般而言,全氟化合物通過干濕沉降、地表徑流、冰雪融水、污水處理廠中水的排放等過程進(jìn)入河流.因此,河流中全氟化合物的分布特征和物質(zhì)荷載可以反映全氟化合物的遷移特征以及生產(chǎn)排放信息.然而,目前的研究主要集中在空間分布上,對(duì)季節(jié)變異的研究還相對(duì)缺乏[3-4].如果全氟化合物的季節(jié)動(dòng)態(tài)沒有被考慮到,則它的遷移、傳輸和物質(zhì)荷載評(píng)估的不確定性將會(huì)增加.而關(guān)于流域中全氟化合物季節(jié)變異的規(guī)律還存在一定的爭議,在歐洲易北河的研究中,全氟己酸(PFHxA)和PFOA表現(xiàn)出了顯著的季節(jié)變異特征[5];在日本流域的研究中,只有PFOS和PFOA在部分流域表現(xiàn)出了季節(jié)變異[6].也有研究發(fā)現(xiàn)全氟化合物的濃度和組成特征表現(xiàn)出了不同的季節(jié)特征,然而此結(jié)論還缺乏統(tǒng)計(jì)上的驗(yàn)證[7].此外,在西班牙的研究發(fā)現(xiàn),流域中全氟化合物在整年并沒有出現(xiàn)季節(jié)變異[8].因此,同時(shí)研究河流中全氟化合物的空間分布和季節(jié)變異對(duì)于了解全氟化合物的遷移以及物質(zhì)荷載評(píng)估有重要意義.

    青藏高原被稱為“亞洲水塔”,其為周邊國家和區(qū)域的數(shù)億人口提供了水資源.離子型的全氟化合物又具有水溶性和移動(dòng)性,然而目前關(guān)于青藏高原水體中全氟化合物濃度水平、時(shí)空分布特征和物質(zhì)荷載的研究還是空白.因此,本研究采集了拉薩河和拉薩市污水處理廠的樣品,旨在:調(diào)查拉薩河中全氟化合物的濃度水平;確定拉薩河中全氟化合物的時(shí)空分布特征;估算污水處理廠排放到拉薩河和拉薩河排放到雅魯藏布江的全氟化合物的物質(zhì)荷載.這將為評(píng)估拉薩水環(huán)境質(zhì)量提供數(shù)據(jù)支持,為保障亞洲水塔環(huán)境安全提供決策依據(jù).

    1 材料和方法

    1.1 樣品的采集

    分別在2016年豐水期(6月)和枯水期(12月)采集了拉薩河27個(gè)點(diǎn)的表層水以及拉薩污水處理廠的中水(圖1).拉薩河是世界上海拔最高的河流之一,為3600m~5500m,全長568km,流域面積31760km2,最小流量8.95m3/s,最大流量2830m3/s,年平均流量288m3/s.它發(fā)源于念青唐古拉山脈中段北側(cè),在拉薩市南郊匯入雅魯藏布江.拉薩污水處理廠承接了拉薩市污水處理的功能,一期工程日均處理量5.5t.實(shí)際上,污水處理廠屬于超負(fù)荷運(yùn)載,實(shí)際產(chǎn)生的污水量要大于處理量,污水處理量也沒有季節(jié)上的差異.處理后的中水直接排入了拉薩河.

    圖1 采樣點(diǎn)示意

    W、N、E、S分別代表西、北、東、南支的點(diǎn).數(shù)字從小到大表示從上游到下游

    每個(gè)采樣點(diǎn)用聚丙烯塑料瓶(PP瓶)采集2.5L水,并送往中國科學(xué)院青藏高原研究所避光存儲(chǔ)在4°C的冷藏室中待處理.PP瓶在采樣前均用甲醇、超純水和采樣點(diǎn)處的水預(yù)先清洗3次,然后再進(jìn)行樣品的采集.采樣前3d均沒有發(fā)生降水事件.

    1.2 樣品的前處理

    采集的水樣先用高溫除雜后的玻璃纖維濾膜(沃特曼,英國)進(jìn)行過濾,然后向過濾好的2L水樣中加入回收率物質(zhì)PFAC-MAX(惠靈頓實(shí)驗(yàn)室,加拿大),再用WAX柱(沃特世,6cc,150mg,30μm,美國)以1~2滴/s的速率進(jìn)行固相萃取(色譜科,SPE-12,美國).在上樣前先分別用5mL甲醇和5mL超純水對(duì)WAX柱進(jìn)行活化.待上完樣,WAX柱用4mL醋酸銨緩沖溶液進(jìn)行除雜,然后串聯(lián)一個(gè)WAX柱抽真空干燥30min.此過程與其他研究中的處理過程相似[9].處理好的WAX柱用20mL 0.1%NH4OH的甲醇溶液進(jìn)行洗脫,最后將洗脫液氮吹濃縮到200μL,待上機(jī)分析測(cè)試.

    1.3 儀器分析

    本研究分析測(cè)試的目標(biāo)化合物共有17種,分別為:全氟丁酸(PFBA)、全氟戊酸(PFPeA)、全氟己酸(PFHxA)、全氟庚酸(PFHpA)、全氟辛酸(PFOA)、全氟任酸(PFNA)、全氟癸酸(PFDA)、全氟十一烷酸(PFUnA)、全氟十二烷酸(PFDoA)、全氟十三烷酸(PFTrDA)、全氟十四烷酸(PFTeDA)、全氟十六烷酸(PFUnDA)、全氟十八烷酸(PFODA)、全氟丁烷磺酸(PFBS)、全氟戊烷磺酸(PFHxS)、全氟辛烷磺酸(PFOS)、全氟癸烷磺酸(PFDS).

    本研究所用的儀器為雙三元二維液相色譜-質(zhì)譜聯(lián)用儀(UltiMate 3000DGLC/(?)ESI-MS/MS) (賽默飛,美國).該儀器有較高的靈敏度和較低的化學(xué)噪音,定性定量結(jié)果可靠穩(wěn)定、便捷快速,且檢測(cè)范圍廣.該分析方法可以較好地應(yīng)用到水樣中全氟化合物的檢測(cè)分析.儀器左泵輸送的流動(dòng)相為(A)0.1%的甲酸溶液、(B)乙腈、(C)蒸餾水;右泵輸送的流動(dòng)相為(A)30mmol/L甲酸銨緩沖溶液(pH=4)、(B)乙腈、(C)蒸餾水.流動(dòng)相流速為2mL/min.色譜柱為Acclaim 120C18(賽默飛,4.6×150mm,5.0μm,美國).進(jìn)樣量為10μl,進(jìn)樣速率為0.8mL/min.質(zhì)譜為電噴霧離子源(HESI源)負(fù)離子模式.質(zhì)譜的鞘氣設(shè)置為40units;輔助氣為12units;源電壓2500V;汽化溫度350℃;離子傳輸管溫度400℃;掃描時(shí)間0.01s.

    1.4 質(zhì)量控制

    為了進(jìn)行質(zhì)量控制,萃取樣品前在樣品中均加入了同位素標(biāo)記的回收率物質(zhì).加入樣品中的回收率如下:13C4-PFBA為40%±27%(平均值±標(biāo)準(zhǔn)偏差),13C2-PFHxA為75%±27%,13C4-PFOA為95%±31 %,18O2-PFHxS為88%±28 %,13C5-PFNA為92%±32%,13C2-PFDA為91%±35%,13C4-PFOS為90%±33%,13C2-PFUnA為73%±31%,13C2-PFDoA為70%±32%.為了進(jìn)行全氟化合物分布特征在空間和時(shí)間上的比較,所有的樣品均用回收率進(jìn)行校正.

    此外,用超純水設(shè)置了空白.空白的處理分析流程和樣品處理流程相一致.方法檢出限(MDLs)定義為空白濃度的算數(shù)平均值加3倍的標(biāo)準(zhǔn)偏差.空白樣品中全氟化合的濃度和方法檢出限如表1所示.在進(jìn)樣過程中,每分析10個(gè)樣品進(jìn)1次溶劑空白,以檢測(cè)可能的污染.

    表1 空白樣品中全氟化合的濃度和方法檢出限(pg/L)

    1.5 數(shù)據(jù)分析

    樣品的濃度均為扣除空白后的值.扣除空白后濃度為負(fù)值的點(diǎn)按0計(jì)算,大于0小于檢測(cè)限的值按MDL/2計(jì)算.統(tǒng)計(jì)檢驗(yàn)用SPSS 22.0(IBM公司,美國)進(jìn)行,顯著性水平設(shè)定為=0.05,雙尾檢驗(yàn).圖形制作用Origin 9.0(Origin Lab,美國)和ArcMap 10.4(ESRI,美國).聚類分析用R 3.3.3.

    1.6 物質(zhì)荷載評(píng)估

    物質(zhì)荷載用式(1)來計(jì)算:

    ML=

    C

    ×

    F

    (1)

    式中ML為PFASs的年均荷載,單位t/a;為拉薩河水體中PFASs的平均濃度,單位pg/L;為拉薩河的水流量,單位m3/s.

    2 結(jié)果與討論

    2.1 濃度水平和組成特征

    本研究共收集了27個(gè)采樣點(diǎn)的54個(gè)河流水樣,所有河流樣品中PFTeDA、PFHxDA、PFODA和PFDS均在MDLs以下,PFUdA、PFDoA和PFTrDA分別有61%、61%和98%的樣品在MDLs以下,其他化合物均有較高的檢出頻率.河流樣品中Σ13PFASs的平均值為322pg/L,濃度范圍從60到1724pg/L,中位數(shù)是270pg/L,表示總平均濃度水平呈偏左分布.各化合物的平均水平從大到小依次為(pg/L):PFBA (80)>PFPeA(71)>PFOS>(46)>PFOA(44)>PFBS(33)>PFHxA(16)>PFHxS(12)>PFHpA(10)>PFNA(4)>PFDA (3)>PFUdA(1)>PFDoA(0.4)>PFTrDA(0.04).

    拉薩河流中全氟化合物的濃度水平相比其他城市流域的低,但比南北極表面水體的濃度高.拉薩河平均Σ10PFASs的濃度比美國、歐洲、亞洲等地區(qū)的河流的低1個(gè)數(shù)量級(jí)(表2).全氟化合物的濃度水平的高低與工業(yè)發(fā)展水平和人口密度有直接的關(guān)系[4].拉薩河表層水較低的濃度水平可能歸因于較低的工業(yè)發(fā)展水平和較少的人口密度.根據(jù)國務(wù)院發(fā)展規(guī)劃,到2020年,拉薩市人口不超過500000,這將決定相對(duì)少量全氟產(chǎn)品的需求.

    拉薩河全氟化合物的濃度水平與同樣處于偏遠(yuǎn)地區(qū)的表層水體相比,它略高于南北極等偏遠(yuǎn)地區(qū)表層水體(表2).南北極等偏遠(yuǎn)地區(qū)的全氟化合物主要來源于大氣或者洋流的遠(yuǎn)距離傳輸.拉薩深處內(nèi)陸,遠(yuǎn)離海洋,不存在洋流的遠(yuǎn)距離傳輸.拉薩河中的全氟化合物除了本地排放外,可能還有部分與大氣的遠(yuǎn)距離傳輸有關(guān)系.

    表2 不同流域全氟化合物的濃度水平和組成特征

    注: “–”表示化合物沒有被分析; “n.d.” 表示未檢出; “

    組成特征是評(píng)價(jià)全氟化合物污染狀態(tài)的另一個(gè)指標(biāo).根據(jù)各全氟化合物濃度的平均水平,拉薩河中短鏈的全氟化合物(C£8)占了總濃度的97%,這說明短鏈的全氟產(chǎn)品在拉薩的生產(chǎn)和使用量多于長鏈的.也有可能是由于溶水性的長鏈全氟化合物更傾向于分配到沉積物和懸浮顆粒物中[18],因而溶水性的短鏈化合物在水相占有較大的比例.PFBA、PFPeA、PFOA、PFBS和PFOS是主要化合物.PFOA和PFOS的比例均為14%,而PFOA和PFOS的替代品PFBA和PFBS的比例分別為25%和13%.另一種主要化合物是PFPeA,其占比為22%.其余的全氟化合物占比則低于10%,分別為PFHxA(5%)> PFHxS(4%)>PFHpA(3%)>PFNA(1.3%)>PFDA(1.1%) >PFUdA(0.2%)>PFDoA(0.1%)>PFTrDA(0.01%).說明拉薩PFOA和PFOS正逐漸地被短鏈的化合物取代,C4和C5類的短鏈化合物成為了主流,這可能是由于區(qū)域消減PFOS和PFOA的結(jié)果,也可能是由于西藏全氟產(chǎn)品,例如織物、印染材料等生產(chǎn)使用的轉(zhuǎn)型.全氟化合物的消減已經(jīng)有了實(shí)質(zhì)性的進(jìn)展,這與其他地方最新的研究結(jié)果[3,7]相一致.

    2.2 時(shí)空分布特征

    根據(jù)拉薩河各采樣點(diǎn)不同時(shí)期全氟化合物濃度的平均值,其空間分布特征如圖2所示.拉薩河?xùn)|、北、西、南支的平均Σ13PFASs濃度分別為218、667、323、380pg/L,北、西和南支的濃度均高于東支.圖2-h中,拉薩河北支的濃度最高,東支的濃度最低.南支在污水處理廠的下游,污水廠的中水是一個(gè)點(diǎn)排放源,而北支的濃度要高于南支,說明北支上存在其他的排放源.此外,西支的濃度高于東支,這說明拉薩河西支存在其他的排放源.西支沿線的堆龍德慶區(qū)內(nèi)有廠房、養(yǎng)雞場等,這可能是一個(gè)影響因素.

    如果存在面源的排放,下游的濃度特征應(yīng)和上游基本保持一致,如果存在一個(gè)點(diǎn)源排放,下游的濃度會(huì)隨著稀釋擴(kuò)散作用而降低[2].本研究中,下游點(diǎn)的濃度并非都隨著稀釋作用而降低.如圖2-a、b圖中達(dá)孜區(qū)內(nèi),E4 PFBA和PFBS的濃度比E3、E2、E1的高,E2 PFBA和PFBS的濃度比E1高.圖2-e堆龍德慶區(qū)內(nèi),W2 PFHxA的濃度比W1高,W4比W3高.圖2-c、d、f、g說明存在明顯的點(diǎn)源排放.c中N3 PFOA的濃度明顯高于其他點(diǎn),d中S2 PFOS的濃度明顯高于其他點(diǎn),f中S1 PFHxS有較高的濃度,g中N2 PFPeA有較高的濃度.N3、N2分別位于生活區(qū),S2位于污水廠下游.對(duì)圖2-a、b和c、d進(jìn)行對(duì)比,PFBA和PFBS的空間分布相對(duì)均勻,而PFOA和PFOS在個(gè)別點(diǎn)處的濃度明顯較高,說明基于PFBA和PFBS的全氟產(chǎn)品的使用比較廣泛,而基于PFOA和PFOS的全氟產(chǎn)品只在某些區(qū)域或行業(yè)有相對(duì)多的使用量或排量.在天津大沽河道發(fā)現(xiàn)較高濃度的全氟化合物,這可能來源于污水或是工業(yè)廢水的點(diǎn)源排放,而山東濰坊河流中全氟化合物較高的濃度水平可能是由于受到非點(diǎn)源的降水和地表徑流的影響[19].影響拉薩河全氟化合物空間分布差異的因素并非是單一的,可能還受到多種因素的影響,需要進(jìn)一步分析.

    圖2 拉薩河全氟化合物的空間分布特征

    拉薩河表層水中全氟化合物的濃度水平在豐水期和枯水期存在季節(jié)性的差異.拉薩河表層水中平均Σ13PFASs的濃度在豐水期和枯水期分別為479和164pg/L,豐水期是枯水期的2.9倍.拉薩河?xùn)|支平均Σ13PFASs的濃度在豐水期和枯水期分別為318和119pg/L,豐水期是枯水期的2.7倍.拉薩河北支平均Σ13PFASs的濃度在豐水期和枯水期分別為911和423pg/L,豐水期是枯水期的2.2倍.拉薩河西支平均Σ13PFASs的濃度在豐水期和枯水期分別為517和129pg/L,其豐水期的濃度是枯水期的4倍.同樣的,南支平均Σ13PFASs的濃度在豐水期和枯水期分別為602和166pg/L,豐水期的是枯水期的3.6倍.

    圖3 拉薩河表層水中全氟化合物的季節(jié)變異

    **表示不同季節(jié)間差異極顯著,<0.01;*表示不同季節(jié)間差異顯著,<0.05

    為了進(jìn)一步確定全氟化合物的季節(jié)變異,對(duì)不同化合物在枯水期和豐水期進(jìn)行獨(dú)立樣本T檢驗(yàn).結(jié)果表明PFBA、PFHxA、PFHpA、PFBS、PFHxS和PFOS的濃度水平在豐水期顯著(<0.05)大于枯水期(圖3),這意味著季節(jié)性的人類活動(dòng)對(duì)這幾種全氟化合物的排放有著重要的影響.對(duì)采樣點(diǎn)在不同的季節(jié)進(jìn)行配對(duì)T檢驗(yàn),結(jié)果表明E2、E5、E6、E7、E8、E9、E12、E13、W4、W6、S1、和S5全氟化合物的濃度在豐水期和枯水期差異顯著(<0.05).河流中全氟化合物的季節(jié)動(dòng)態(tài)主要受市場上全氟產(chǎn)品季節(jié)性需求的影響,也受地表徑流和突發(fā)性降水事件的影響[5,20].由于在拉薩河采樣的前幾天均無降水事件發(fā)生,因而拉薩河一些點(diǎn)全氟化合物排放特征在豐水期和枯水期的顯著差異是由其生產(chǎn)和使用的季節(jié)性差異引起的.這也反映出拉薩的生產(chǎn)生活活動(dòng)受季節(jié)的影響較大.歐洲易北河的研究也發(fā)現(xiàn)PFHxA和PFOA的濃度水平在不同季節(jié)有顯著的差異[5],然而西班牙河流中全氟化合物的濃度水平在全年都比較穩(wěn)定[8].因而,全氟化合物的季節(jié)變異還存在區(qū)域上的差異性.全氟化合物區(qū)域存在差異的因素還需要進(jìn)一步分析.

    2.3 影響因素分析

    影響全氟化合物時(shí)空分布差異的因素是多元的,人口則是諸多因素中的一個(gè)因素.拉薩河全氟化合物的分布特征與人口數(shù)量之間的關(guān)系如圖4所示.總濃度和人口數(shù)量之間無規(guī)律.PFOA和PFOS的濃度水平隨著人口數(shù)量的增多而增大,而PFBA的濃度水平隨著人口數(shù)量的增多而降低.PFPeA在人口最多的區(qū)域濃度比較低,而在人口最少的區(qū)域濃度最大.城關(guān)區(qū)主要為人口聚居的區(qū)域,而達(dá)孜區(qū)主要為工業(yè)發(fā)展區(qū),說明在拉薩PFOA和PFOS主要與人們的生活活動(dòng)有關(guān),在生活活動(dòng)中還有較大的使用量,而PFBA和PFPeA主要與工業(yè)生產(chǎn)有關(guān).中國南部沿海地區(qū)的研究發(fā)現(xiàn)人口數(shù)量、人類活動(dòng)和國民生產(chǎn)總值對(duì)該區(qū)全氟化合物的污染有極大的影響[21].而在中國經(jīng)濟(jì)發(fā)達(dá)的地區(qū),人口數(shù)量與全氟羧酸的濃度水平有顯著的正相關(guān)關(guān)系,而與總?cè)衔锏臐舛葲]有顯著的相關(guān)性[22].

    具有相近全氟化合物排放特征或來源的點(diǎn)可以聚為一類[20].對(duì)河流不同采樣點(diǎn)在豐水期和枯水期的濃度分別進(jìn)行聚類分析,如圖5所示.在豐水期,S2和N2各為一類,其他點(diǎn)被聚為一類.說明S2和N2的全氟化合物有各自的來源,S2和N2與其他點(diǎn)全氟化合物的來源也不一致.S2在污水處理廠的下游,污水廠的中水對(duì)其有較大的影響.N2在居民區(qū)內(nèi),該段水域可能存在生活污水直接排入的情況.在枯水期,N3為一類,E4、S2和W2被聚為一類,其他點(diǎn)被聚為一類,而這些點(diǎn)又分布在不同的區(qū)域,說明E4、S2和W2存在點(diǎn)源排放的情況.S2位于污水處理廠的下游,E4位于達(dá)孜工業(yè)園的中游,可能這兩個(gè)點(diǎn)會(huì)受到污水廠和工業(yè)園較大的影響.W2沿線存在建筑施工活動(dòng),可能這是影響W2的一個(gè)主要因素.對(duì)美國羅德島和紐約表層水樣中的全氟化合物進(jìn)行聚類和主成分分析,結(jié)果發(fā)現(xiàn)一些點(diǎn)排放源,如:機(jī)場滅火劑的使用和紡織廠廢水的排放,可以聚為一類,一些采樣點(diǎn)可能受到不同來源的全氟化合物的影響,而各類源的貢獻(xiàn)度不同[23].

    圖4 拉薩河全氟化合物的分布特征與人口之間的關(guān)系

    豐水期和枯水期聚類的結(jié)果不一致,說明拉薩全氟化合物的生產(chǎn)和排放模式在季節(jié)上不固定,也就是說該區(qū)域生產(chǎn)、生活活動(dòng)受到季節(jié)的影響較大.拉薩是一個(gè)旅游城市,拉薩獨(dú)特的自然風(fēng)光和人文景觀吸引了大量的游客前來參觀,在豐水期,旅游活動(dòng)比枯水期興盛,從而豐水期全氟化合物的生產(chǎn)、生活活動(dòng)和排放量也比枯水期多.

    圖5 拉薩河全氟化合物的聚類分析

    A,豐水期; B,枯水期

    拉薩河全氟化合物季節(jié)變異的另一個(gè)重要原因是水量的變異.拉薩河流域?qū)儆诩撅L(fēng)影響的區(qū)域,在豐水期,由于冰雪融化和降雨致使拉薩河有大量的水源補(bǔ)給,從而其流量較大.拉薩河流量最大的時(shí)候達(dá)到2830m3/s.而冰雪和降雨中都有可能攜帶了全氟化合物,這對(duì)全氟化合物的傳輸有重要的影響[24].與日本琵琶湖的研究相比較,除了2個(gè)點(diǎn)的PFOS和PFOA存在季節(jié)變異,其他12個(gè)點(diǎn)均不存在顯著的季節(jié)變異,這可能與流量的變化較小有關(guān)[6].而易北河全氟化合物季節(jié)變異可能是由河流流量的季節(jié)性差異造成的,易北河的流量在2月份達(dá)到最大值,1528m3/s,而PFHxA和PFOA的濃度相應(yīng)在2月份最低[5].

    污水處理廠是全氟化合物一個(gè)重要的點(diǎn)排放源.污水入水中全氟化合物的濃度水平與經(jīng)濟(jì)發(fā)展水平顯著呈正相關(guān),而且入水中全氟化合物的濃度對(duì)中水中的濃度有著絕對(duì)影響[3].2016年,拉薩國民生產(chǎn)總值排名在27個(gè)中國省會(huì)城市中占第27位,且拉薩的產(chǎn)業(yè)以第三產(chǎn)業(yè)為主.較低的經(jīng)濟(jì)發(fā)展水平直接決定了拉薩污水處理廠入水中全氟化合物濃度水平不高.拉薩污水處理廠入水中全氟化合物的濃度為56ng/L,低于上海的134ng/L[22].入水中的主要化合物為PFOA和PFOS,這與化工、電子、造紙、紡織等不同類型的污水處理廠中的研究相一致[25].經(jīng)過污水處理過程,可能一些前體化合物發(fā)生了轉(zhuǎn)變,或者長鏈的化合物降解成了短鏈的化合物,部分化合物被吸附到活性污泥中,部分被分配到中水中,隨后中水從污水廠排放到河流中.拉薩污水處理廠中水的排放對(duì)拉薩河的物質(zhì)荷載有重要地影響.

    2.4 物質(zhì)荷載

    污水處理廠會(huì)排放出PFOA和PFOS以及其他的全氟化合物.拉薩污水處理廠的中水中全氟化合物的濃度在豐水期和枯水期分別為34和2.8ng/L.本研究計(jì)算了污水處理廠排放的物質(zhì)荷載,結(jié)果表明,污水處理廠排放到拉薩河的全氟化合物的年均荷載約為37mg/a.拉薩河排放到雅魯藏布江的全氟化合物的物質(zhì)荷載為9~2869t/a,平均為292t/a.相較而言,通過河流排放到渤海灣的全氟化合物的物質(zhì)荷載約為873t/a[26].從黃河到渤海灣的物質(zhì)荷載約為3.9t/a[18].由于較大的流量,拉薩河排放到雅魯藏布江的全氟化合物的量與其他流域的物質(zhì)荷載相比并不低.實(shí)際上,本研究中的物質(zhì)荷載被低估了,因?yàn)樗袡z測(cè)到的是溶解相的全氟化合物,還有少量粒子相的全氟化合吸附在懸浮顆粒物上[18],本研究只計(jì)算了水體中溶解相的全氟化合物.在冰雪中的研究發(fā)現(xiàn),青藏高原全氟化合物有較長時(shí)間的沉積記錄[24],因此,由于氣候變暖導(dǎo)致的冰雪融化將會(huì)釋放出更多的全氟化合物到拉薩河中.氣候變暖會(huì)影響持久性有機(jī)污染物包括全氟化合物的遠(yuǎn)距離傳輸,青藏高原有可能會(huì)變成全氟化合物的一個(gè)匯區(qū),這個(gè)問題應(yīng)當(dāng)被重視,進(jìn)一步的研究亟需開展.青藏高原被稱為“亞洲水塔”,其為周邊國家和地區(qū)提供了大量的水資源,因此,保護(hù)源頭的水環(huán)境質(zhì)量安全舉足輕重.全氟化合物的全球消減應(yīng)該達(dá)成共識(shí),而且這個(gè)共識(shí)需要被強(qiáng)有力地和系統(tǒng)地執(zhí)行.

    3 結(jié)論

    3.1 拉薩地區(qū)普遍受到全氟化合物的污染,拉薩河污染水平遠(yuǎn)低于其他大城市的流域,略高于南北極地區(qū)表面水體.主要化合物為PFBA、PFPeA、PFOA、PFBS和PFOS.

    3.2 在空間分布上,拉薩河北支的濃度最高,東支的濃度最低.由于擴(kuò)散效應(yīng),下游化合物的濃度理論上應(yīng)當(dāng)降低,而拉薩河下游全氟化合物的濃度并非都低于上游,說明拉薩市內(nèi)的污水并未完全納入污水處理廠處理,拉薩河還受到其他源排放的影響.PFBA和PFBS的空間分布相對(duì)均勻,而PFOA和PFOS在個(gè)別點(diǎn)處濃度明顯較高,說明基于PFBA和PFBS這兩種化合物的全氟產(chǎn)品使用比較廣泛,而基于PFOA和PFOS的全氟產(chǎn)品只在某些區(qū)域或行業(yè)有相對(duì)多的使用量或排放量.

    3.3 在季節(jié)特征上,全氟化合物在豐水期和枯水期的濃度水平和組成特征有顯著的差異,表明拉薩全氟化合物的生產(chǎn)和排放模式在季節(jié)上不固定,即言之,拉薩生產(chǎn)、生活活動(dòng)受到季節(jié)的影響較大.

    3.4 污水處理廠是拉薩河受全氟化合物污染的一個(gè)重要的點(diǎn)排放源,污水處理廠排放到拉薩河中的全氟化合物的年排量約為37mg,拉薩河排放到雅魯藏布江的年排量約為292t.

    [1] Wang Z, Cousins I T, Scheringer M, et al. Global emission inventories for C4-C14perfluoroalkyl carboxylic acid (PFCA) homologues from 1951 to 2030, part II: the remaining pieces of the puzzle [J]. Environ Int, 2014,69:166-76.

    [2] Wang T, Vestergren R, Herzke D, et al. Levels, Isomer Profiles, and Estimated Riverine Mass Discharges of Perfluoroalkyl Acids and Fluorinated Alternatives at the Mouths of Chinese Rivers [J]. Environmental Science & Technology, 2016,50(21):11584-11592.

    [3] 周 珍,胡宇寧,史亞利,等.武漢地區(qū)水環(huán)境中全氟化合物污染水平及其分布特征[J]. 生態(tài)毒理學(xué)報(bào), 2017,(3):425-433.

    [4] 樸海濤,陳 舒,焦杏春,等.大運(yùn)河豐水期水體中全氟化合物的分布[J]. 中國環(huán)境科學(xué), 2016,36(10):3040-3047.

    [5] Zhao Z, Xie Z Y, Tang J H, et al. Seasonal variations and spatial distributions of perfluoroalkyl substances in the rivers Elbe and lower Weser and the North Sea [J]. Chemosphere, 2015,129:118-125.

    [6] Tsuda T, Inoue A, Igawa T, et al. Seasonal Changes of PFOS and PFOA Concentrations in Lake Biwa Water [J]. Bulletin of Environmental Contamination and Toxicology, 2010,85(6):593-597.

    [7] Zhu Z Y, Wang T Y, Meng J, et al. Perfluoroalkyl substances in the Daling River with concentrated fluorine industries in China: seasonal variation, mass flow, and risk assessment [J]. Environmental Science and Pollution Research, 2015,22(13):10009-10018.

    [8] Pignotti E, Casas G, Llorca M, et al. Seasonal variations in the occurrence of perfluoroalkyl substances in water, sediment and fish samples from Ebro Delta (Catalonia, Spain) [J]. Science of the Total Environment, 2017,607:933-943.

    [9] Boiteux V, Bach C, Sagres V, et al. Analysis of 29per- and polyfluorinated compounds in water, sediment, soil and sludge by liquid chromatography-tandem mass spectrometry [J]. International Journal of Environmental Analytical Chemistry, 2016,96(8):705-728.

    [10] Pan C G, Ying G G, Zhao J L, et al. Spatiotemporal distribution and mass loadings of perfluoroalkyl substances in the Yangtze River of China [J]. Science of the Total Environment, 2014,493:580-587.

    [11] Yeung L W Y, Yamashita N, Taniyasu S, et al. A survey of perfluorinated compounds in surface water and biota including dolphins from the Ganges River and in other waterbodies in India [J]. Chemosphere, 2009,76(1):55-62.

    [12] Habibullah-Al-Mamun M, Ahmed M K, Raknuzzaman M, et al. Occurrence and distribution of perfluoroalkyl acids (PFAAs) in surface water and sediment of a tropical coastal area (Bay of Bengal coast, Bangladesh) [J]. Science of the Total Environment, 2016,571:1089- 1104.

    [13] Zushi Y, Ye F, Motegi M, et al. Spatially Detailed Survey on Pollution by Multiple Perfluorinated Compounds in the Tokyo Bay Basin of Japan [J]. Environmental Science & Technology, 2011,45(7):2887- 2893.

    [14] Lam N H, Cho C R, Lee J S, et al. Perfluorinated alkyl substances in water, sediment, plankton and fish from Korean rivers and lakes: A nationwide survey [J]. Science of the Total Environment, 2014,491: 154-162.

    [15] Ahrens L, Felizeter S, Sturm R, et al. Polyfluorinated compounds in waste water treatment plant effluents and surface waters along the River Elbe, Germany [J]. Mar Pollut Bull, 2009,58(9):1326-33.

    [16] Levengood J M, Soucek D J, Sass G G, et al. Interspecific and Spatial Comparisons of Perfluorinated Compounds in Bighead and Silver Carp in the Illinois River, Illinois, USA [J]. Bulletin of Environmental Contamination and Toxicology, 2015,95(5):561-566.

    [17] Zhao Z, Xie Z Y, Moller A, et al. Distribution and long-range transport of polyfluoroalkyl substances in the Arctic, Atlantic Ocean and Antarctic coast [J]. Environmental Pollution, 2012,170:71-77.

    [18] Zhao P J, Xia X H, Dong J W, et al. Short- and long-chain perfluoroalkyl substances in the water, suspended particulate matter, and surface sediment of a turbid river [J]. Science of the Total Environment, 2016,568:57-65.

    [19] Yao Y M, Zhu H K, Li B, et al. Distribution and primary source analysis of per- and poly-fluoroalkyl substances with different chain lengths in surface and groundwater in two cities, North China [J]. Ecotoxicology and Environmental Safety, 2014,108:318-328.

    [20] Zhang Y Z, Wang B, Wang W, et al. Occurrence and source apportionment of Per- and poly-fluorinated compounds (PFCs) in North Canal Basin, Beijing [J]. Sci Rep, 2016,6:36683.

    [21] Kwok K Y, Wang X H, Ya M L, et al. Occurrence and distribution of conventional and new classes of per- and polyfluoroalkyl substances (PFASs) in the South China Sea [J]. Journal of Hazardous Materials, 2015,285:389-397.

    [22] Zhang W, Zhang Y T, Taniyasu S, et al. Distribution and fate of perfluoroalkyl substances in municipal wastewater treatment plants in economically developed areas of China [J]. Environmental Pollution, 2013,176:10-17.

    [23] Zhang X M, Lohmann R, Dassuncao C, et al. Source Attribution of Poly- and Perfluoroalkyl Substances (PFASs) in Surface Waters from Rhode Island and the New York Metropolitan Area [J]. Environmental Science & Technology Letters, 2016,3(9):316-321.

    [24] Wang X P, Halsall C, Codling G, et al. Accumulation of Perfluoroalkyl Compounds in Tibetan Mountain Snow: Temporal Patterns from 1980 to 2010 [J]. Environmental Science & Technology, 2014,48(1):173- 181.

    [25] Kim H Y, Seok H W, Kwon H O, et al. A national discharge load of perfluoroalkyl acids derived from industrial wastewater treatment plants in Korea [J]. Science of the Total Environment, 2016,563:530- 537.

    [26] Chen H, Wang X M, Zhang C, et al. Occurrence and inputs of perfluoroalkyl substances (PFASs) from rivers and drain outlets to the Bohai Sea, China [J]. Environmental Pollution, 2017,221:234-243.

    Special distribution and seasonal variation of perfluoroalkyls substances in Lhasa River Basin, China.

    SUN Dian-chao1,3, GONG Ping1,2, WANG Xiao-ping1,2,3*, WANG Chuan-fei1,2

    (1.Key Laboratory of Tibetan Environment Changes and Land Surface Processes, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing 100101, China;2.Chinese Academy of Sciences Center for Excellence in Tibetan Plateau Earth Sciences, Beijing 100101, China;3.University of Chinese Academy of Sciences, Beijing 100049, China)., 2018,38(11):4298~4306

    The levels and compositions of 17 Perfluoroalkyls Substances (PFASs) were determined by UltiMate 3000DGLC-MS/MS for the first time in river water, wastewater of Lhasa, one of the highest plateau city (average elevation 3700m), during wet and dry season in 2016. Lhasa was commonly contaminated by PFASs. The average Σ13PFASs concentration of surface water of Lhasa River was 322pg/L, ranged from 60 to 1724pg/L. The predominant PFASs were short-chain pefluorobutanoate (PFBA) (25%), perfluoropentanoate (PFPeA) (22%), perfluorooctanoic acid (PFOA) (14%), perfluorooctanesulfonate (PFOS) (14%) and perfluorobutanesulfonate (PFBS) (13%). The spatial distribution of PFASs in Lhasa River revealed that the wastewater not totally assimilated into wastewater treatment plant, and Lhasa River basin was mainly polluted by point sources emission. There was a limited transform of production and use from PFOA and PFOS in Tibet. Seasonal variation of profiles and levels of PFASs in both wet and dry seasons suggested that the use and discharge schema of PFASs in Lhasa were influenced by season. The annual discharge load of PFASs from WWTP to Lhasa River was estimated to be 37mg/a, input amount of PFASs from The Lhasa River basin into YarlungZangbo River were estimated to be 292t/a. In the future, all the wastewater should be strictly incorporated into WWTP and PFASs-based product should be restrained.

    perfluoroalkyl acids;POPs;spatial distribution;seasonal variation;mass loading

    X522

    A

    1000-6923(2018)11-4298-09

    孫殿超(1988-),男,甘肅張掖人,博士研究生,主要研究方向?yàn)榄h(huán)境污染與氣候變化.

    2018-04-18

    國家自然科學(xué)基金資助項(xiàng)目(41671480)

    * 責(zé)任作者, 研究員, wangxp@itpcas.ac.cn

    猜你喜歡
    氟化合物拉薩河全氟
    全氟烷基化合物暴露與成年人抑郁癥間的關(guān)系:基于NHANES 2005~2018
    有關(guān)全氟化合物(PFCs)的來源及影響
    談河流水體污染特征及與浮游細(xì)菌群落多樣性響應(yīng)研究進(jìn)展
    拉薩河流域浮游植物生態(tài)學(xué)研究現(xiàn)狀
    拉薩河,從古遠(yuǎn)流來
    布達(dá)拉(2018年4期)2018-05-14 08:53:30
    拉薩市拉薩河城區(qū)段綜合整治工程4~#閘設(shè)計(jì)
    基于概率統(tǒng)計(jì)和模糊綜合評(píng)價(jià)法的水質(zhì)評(píng)價(jià)模型及其應(yīng)用——以拉薩河水質(zhì)評(píng)價(jià)為例
    1種制備全氟聚醚羧酸的方法
    1種制備全氟烯醚磺酰氟化合物的方法
    典型全氟化合物(PFCs)的降解與控制研究進(jìn)展
    91午夜精品亚洲一区二区三区| av国产免费在线观看| 日韩欧美一区视频在线观看 | 99热网站在线观看| 精品久久久久久成人av| 白带黄色成豆腐渣| 中文乱码字字幕精品一区二区三区 | 一级爰片在线观看| 2021少妇久久久久久久久久久| 亚洲国产欧美在线一区| 两个人的视频大全免费| 97超碰精品成人国产| 韩国高清视频一区二区三区| 免费观看无遮挡的男女| 菩萨蛮人人尽说江南好唐韦庄| 你懂的网址亚洲精品在线观看| 久久这里有精品视频免费| 亚洲精品456在线播放app| 搡老妇女老女人老熟妇| 日韩制服骚丝袜av| 日本黄大片高清| 成人一区二区视频在线观看| 99热6这里只有精品| 一级毛片久久久久久久久女| 只有这里有精品99| 午夜精品一区二区三区免费看| 亚洲人与动物交配视频| 久久99热这里只频精品6学生| 日韩中字成人| 国产女主播在线喷水免费视频网站 | 免费黄频网站在线观看国产| 成人性生交大片免费视频hd| 欧美另类一区| 精品久久久久久久久亚洲| av福利片在线观看| 别揉我奶头 嗯啊视频| 国语对白做爰xxxⅹ性视频网站| 亚洲av成人av| 国产精品无大码| 在线a可以看的网站| 国产淫语在线视频| 国产黄色视频一区二区在线观看| 欧美另类一区| 成人国产麻豆网| 国产精品一区二区性色av| 人人妻人人澡人人爽人人夜夜 | 成人午夜精彩视频在线观看| 免费看日本二区| 亚洲四区av| 国产精品爽爽va在线观看网站| 免费播放大片免费观看视频在线观看| 国产精品久久久久久久电影| 熟妇人妻久久中文字幕3abv| 十八禁国产超污无遮挡网站| 男插女下体视频免费在线播放| 偷拍熟女少妇极品色| 亚洲av男天堂| 久久99精品国语久久久| 国产精品一区二区在线观看99 | 亚洲电影在线观看av| 国产精品人妻久久久久久| 丝瓜视频免费看黄片| 国产精品久久久久久精品电影小说 | 观看美女的网站| 亚洲国产精品专区欧美| 日韩av在线免费看完整版不卡| 三级男女做爰猛烈吃奶摸视频| 97超碰精品成人国产| 亚洲人成网站高清观看| 三级国产精品片| 2021少妇久久久久久久久久久| 一个人免费在线观看电影| 亚洲欧洲日产国产| 精品久久久久久久久亚洲| 日韩人妻高清精品专区| 国产男女超爽视频在线观看| 尾随美女入室| 国产一区有黄有色的免费视频 | 国产激情偷乱视频一区二区| 99九九线精品视频在线观看视频| 青青草视频在线视频观看| 一区二区三区四区激情视频| 日日摸夜夜添夜夜爱| 免费观看无遮挡的男女| 国语对白做爰xxxⅹ性视频网站| 免费黄色在线免费观看| 又爽又黄无遮挡网站| 99久久人妻综合| 联通29元200g的流量卡| 久久久久免费精品人妻一区二区| 久久久久久久久久人人人人人人| 男女啪啪激烈高潮av片| 少妇丰满av| 街头女战士在线观看网站| 亚洲18禁久久av| 国产亚洲av嫩草精品影院| 丝袜喷水一区| 国产黄色免费在线视频| 午夜福利视频精品| 国产色婷婷99| 欧美日韩在线观看h| 日本爱情动作片www.在线观看| 又爽又黄无遮挡网站| 日本熟妇午夜| 国产不卡一卡二| 国产免费视频播放在线视频 | 精品午夜福利在线看| 大香蕉97超碰在线| 看黄色毛片网站| 免费高清在线观看视频在线观看| 日本免费a在线| 久久久久久久国产电影| 国产一级毛片七仙女欲春2| xxx大片免费视频| 亚洲精品日本国产第一区| 午夜免费观看性视频| 日韩av不卡免费在线播放| 噜噜噜噜噜久久久久久91| 久久人人爽人人片av| 国产午夜精品论理片| 亚洲乱码一区二区免费版| 欧美一区二区亚洲| 如何舔出高潮| 成人亚洲欧美一区二区av| 国产精品一及| 国产成人aa在线观看| 日韩成人伦理影院| 精品久久久久久电影网| 水蜜桃什么品种好| 高清视频免费观看一区二区 | 午夜久久久久精精品| 国产精品久久久久久久电影| 国产亚洲午夜精品一区二区久久 | 免费av观看视频| 自拍偷自拍亚洲精品老妇| 日韩成人av中文字幕在线观看| 欧美成人精品欧美一级黄| 国产毛片a区久久久久| 国产精品av视频在线免费观看| 高清欧美精品videossex| 亚洲丝袜综合中文字幕| 亚洲精品视频女| 国产精品女同一区二区软件| 中文字幕人妻熟人妻熟丝袜美| 人妻夜夜爽99麻豆av| 国产一区二区亚洲精品在线观看| 天堂俺去俺来也www色官网 | 狂野欧美激情性xxxx在线观看| 精品久久久精品久久久| 亚洲av免费高清在线观看| 国产三级在线视频| 国产一区二区三区av在线| 免费电影在线观看免费观看| 久久久精品欧美日韩精品| 久久人人爽人人片av| 国产午夜精品一二区理论片| 纵有疾风起免费观看全集完整版 | 97在线视频观看| 一区二区三区免费毛片| 麻豆精品久久久久久蜜桃| 日韩欧美国产在线观看| 久久鲁丝午夜福利片| 国产av在哪里看| 一边亲一边摸免费视频| 国产在视频线精品| 国产精品一二三区在线看| 国产乱人偷精品视频| 国产伦一二天堂av在线观看| 午夜日本视频在线| 精品久久久久久久久av| 日本欧美国产在线视频| 天堂影院成人在线观看| 岛国毛片在线播放| 国产精品福利在线免费观看| 国产亚洲午夜精品一区二区久久 | 成年免费大片在线观看| 久久亚洲国产成人精品v| 男女啪啪激烈高潮av片| 黑人高潮一二区| 乱人视频在线观看| 好男人在线观看高清免费视频| 天堂影院成人在线观看| 亚洲av二区三区四区| 日韩伦理黄色片| 69人妻影院| 一级av片app| 午夜免费男女啪啪视频观看| 亚洲欧美日韩东京热| 欧美精品一区二区大全| 久久这里只有精品中国| av在线播放精品| 国产麻豆成人av免费视频| 成人一区二区视频在线观看| 国产免费一级a男人的天堂| 国产成人91sexporn| 免费观看无遮挡的男女| 午夜激情福利司机影院| 国产淫语在线视频| 精品国产露脸久久av麻豆 | 国产毛片a区久久久久| 老司机影院毛片| 一区二区三区高清视频在线| 26uuu在线亚洲综合色| 乱人视频在线观看| 在线免费观看不下载黄p国产| 精品久久久久久成人av| 麻豆成人av视频| 日韩av在线免费看完整版不卡| 国产成年人精品一区二区| 久久韩国三级中文字幕| 亚洲精品中文字幕在线视频 | 亚洲人成网站在线观看播放| av网站免费在线观看视频 | 亚洲精品一二三| 在线观看一区二区三区| 亚洲精品第二区| 色综合色国产| 激情 狠狠 欧美| 又爽又黄a免费视频| 神马国产精品三级电影在线观看| 春色校园在线视频观看| 久久久久国产网址| 国产精品爽爽va在线观看网站| 直男gayav资源| 永久免费av网站大全| 成人漫画全彩无遮挡| 国产综合懂色| 国产 一区精品| 久久久色成人| 国产精品一区二区在线观看99 | 欧美一级a爱片免费观看看| av播播在线观看一区| 国产亚洲av嫩草精品影院| 搞女人的毛片| 有码 亚洲区| 亚洲av不卡在线观看| av在线蜜桃| 搡老妇女老女人老熟妇| 国产精品人妻久久久久久| 国产成人免费观看mmmm| 精品久久久久久成人av| 日韩欧美精品v在线| 精品欧美国产一区二区三| 久久精品国产亚洲av涩爱| 高清日韩中文字幕在线| 特大巨黑吊av在线直播| 国产黄色免费在线视频| 国产精品久久久久久精品电影小说 | 久久久久性生活片| 1000部很黄的大片| 丰满少妇做爰视频| 亚洲性久久影院| 99久久中文字幕三级久久日本| 国产精品国产三级国产专区5o| 国产有黄有色有爽视频| av卡一久久| 国产亚洲精品av在线| 免费不卡的大黄色大毛片视频在线观看 | 国产有黄有色有爽视频| 毛片一级片免费看久久久久| 日韩一区二区视频免费看| 在线观看一区二区三区| 久久草成人影院| av网站免费在线观看视频 | 国产午夜精品论理片| 亚洲av福利一区| 晚上一个人看的免费电影| 在线观看av片永久免费下载| 男人和女人高潮做爰伦理| 国产精品人妻久久久久久| 99热网站在线观看| a级毛色黄片| 欧美三级亚洲精品| 欧美zozozo另类| 亚洲av中文字字幕乱码综合| 黑人高潮一二区| 精品久久久久久成人av| 亚洲成人中文字幕在线播放| 欧美丝袜亚洲另类| 我的老师免费观看完整版| 两个人视频免费观看高清| 午夜视频国产福利| 2018国产大陆天天弄谢| 国产色婷婷99| 亚洲精品乱久久久久久| 亚洲图色成人| 国产一区二区在线观看日韩| 久久草成人影院| 午夜福利在线观看免费完整高清在| 国产av国产精品国产| 午夜免费观看性视频| av免费在线看不卡| 久久精品国产亚洲av天美| 亚洲精品,欧美精品| av女优亚洲男人天堂| 免费av不卡在线播放| 美女xxoo啪啪120秒动态图| 青春草视频在线免费观看| 亚洲欧洲国产日韩| 如何舔出高潮| 黄片无遮挡物在线观看| 三级国产精品片| 边亲边吃奶的免费视频| 好男人在线观看高清免费视频| 国产男女超爽视频在线观看| 大片免费播放器 马上看| 特级一级黄色大片| 欧美日韩视频高清一区二区三区二| 大又大粗又爽又黄少妇毛片口| 久久这里只有精品中国| 免费黄网站久久成人精品| 深夜a级毛片| 久久精品夜色国产| 成人毛片a级毛片在线播放| 午夜免费激情av| 亚洲色图av天堂| 欧美性感艳星| 午夜免费男女啪啪视频观看| 精品一区在线观看国产| 免费人成在线观看视频色| 亚洲av国产av综合av卡| 丰满少妇做爰视频| 免费av不卡在线播放| 两个人视频免费观看高清| 国产探花在线观看一区二区| 亚洲乱码一区二区免费版| 国产精品久久久久久av不卡| 精品人妻一区二区三区麻豆| 国产精品人妻久久久久久| 国产一区二区在线观看日韩| 日本一二三区视频观看| 欧美高清成人免费视频www| 日本午夜av视频| 亚洲最大成人手机在线| 国产成人午夜福利电影在线观看| 国内揄拍国产精品人妻在线| 在线 av 中文字幕| 人妻制服诱惑在线中文字幕| 五月玫瑰六月丁香| 99热全是精品| 亚洲成色77777| 精品久久国产蜜桃| 久久草成人影院| 日本黄大片高清| 久久久久九九精品影院| av.在线天堂| 国语对白做爰xxxⅹ性视频网站| 人人妻人人澡人人爽人人夜夜 | 赤兔流量卡办理| 国产亚洲午夜精品一区二区久久 | 国产永久视频网站| 乱系列少妇在线播放| 最近中文字幕2019免费版| 欧美激情在线99| 国产一区二区三区综合在线观看 | 成人特级av手机在线观看| 欧美潮喷喷水| 91久久精品国产一区二区三区| 久久久精品免费免费高清| 国产精品人妻久久久影院| av在线老鸭窝| 久久人人爽人人爽人人片va| 日韩欧美一区视频在线观看 | av在线亚洲专区| 激情 狠狠 欧美| 欧美高清成人免费视频www| 国产亚洲最大av| 精品国内亚洲2022精品成人| 久久久久久久午夜电影| 日日摸夜夜添夜夜添av毛片| 18禁在线播放成人免费| 十八禁网站网址无遮挡 | 欧美三级亚洲精品| 久久久久网色| 久久精品人妻少妇| 久久久久精品性色| av网站免费在线观看视频 | 夜夜看夜夜爽夜夜摸| 欧美成人一区二区免费高清观看| 老司机影院毛片| 欧美成人一区二区免费高清观看| 亚洲电影在线观看av| 99久久精品国产国产毛片| 欧美三级亚洲精品| 能在线免费观看的黄片| 亚洲激情五月婷婷啪啪| 久久久久性生活片| 美女大奶头视频| 国产精品久久久久久久久免| 亚洲av成人精品一二三区| 色尼玛亚洲综合影院| www.色视频.com| 中文字幕人妻熟人妻熟丝袜美| 高清日韩中文字幕在线| 亚洲三级黄色毛片| 亚洲精品国产成人久久av| a级一级毛片免费在线观看| 男女国产视频网站| 国产伦理片在线播放av一区| 国产黄色小视频在线观看| 免费大片黄手机在线观看| 91久久精品国产一区二区成人| 国产毛片a区久久久久| 男女视频在线观看网站免费| 人体艺术视频欧美日本| 婷婷色综合www| 身体一侧抽搐| 插阴视频在线观看视频| 亚洲欧美日韩东京热| 免费大片18禁| 九草在线视频观看| 一区二区三区免费毛片| 日本黄大片高清| 国产在视频线精品| 白带黄色成豆腐渣| 成人午夜高清在线视频| 成人亚洲欧美一区二区av| 一二三四中文在线观看免费高清| 国产伦理片在线播放av一区| a级毛片免费高清观看在线播放| 国产午夜精品一二区理论片| 观看美女的网站| 黄色一级大片看看| 好男人视频免费观看在线| 国产熟女欧美一区二区| 国产精品一区www在线观看| 天天躁夜夜躁狠狠久久av| 日韩一区二区三区影片| 欧美bdsm另类| av在线老鸭窝| 久久久久久久午夜电影| 久久精品国产亚洲av涩爱| 亚洲aⅴ乱码一区二区在线播放| 三级男女做爰猛烈吃奶摸视频| 午夜福利视频1000在线观看| 99久久精品一区二区三区| 七月丁香在线播放| 成人一区二区视频在线观看| 偷拍熟女少妇极品色| www.av在线官网国产| 欧美xxxx黑人xx丫x性爽| 精品熟女少妇av免费看| 亚洲成人久久爱视频| 免费看a级黄色片| 特大巨黑吊av在线直播| 日韩大片免费观看网站| 一级毛片aaaaaa免费看小| 97在线视频观看| 色5月婷婷丁香| 69av精品久久久久久| 国产精品1区2区在线观看.| 久久国内精品自在自线图片| 久久热精品热| 亚洲精品456在线播放app| 亚洲国产av新网站| 国产探花在线观看一区二区| 偷拍熟女少妇极品色| 成年版毛片免费区| 精品欧美国产一区二区三| 亚洲av中文av极速乱| 国产一级毛片在线| 我的女老师完整版在线观看| 成人欧美大片| 日本爱情动作片www.在线观看| 欧美最新免费一区二区三区| 日本欧美国产在线视频| 国产av在哪里看| 久久久久久久国产电影| 精品午夜福利在线看| 18+在线观看网站| 久久99热这里只频精品6学生| 老司机影院毛片| 中文乱码字字幕精品一区二区三区 | 国产在视频线精品| 亚洲精品亚洲一区二区| 亚洲av不卡在线观看| 色吧在线观看| 久久久精品免费免费高清| 观看美女的网站| 成人亚洲欧美一区二区av| 简卡轻食公司| 亚洲精品成人久久久久久| 高清在线视频一区二区三区| 97在线视频观看| 国产一区二区三区综合在线观看 | 欧美zozozo另类| 夜夜看夜夜爽夜夜摸| 国产美女午夜福利| 91av网一区二区| 午夜久久久久精精品| 国产一区二区亚洲精品在线观看| 免费大片黄手机在线观看| 毛片女人毛片| 亚洲欧美精品专区久久| av线在线观看网站| 久久久久久国产a免费观看| 亚洲精华国产精华液的使用体验| 色哟哟·www| 欧美潮喷喷水| 日本猛色少妇xxxxx猛交久久| 亚洲美女搞黄在线观看| 看非洲黑人一级黄片| 久久草成人影院| 中文欧美无线码| 成年免费大片在线观看| av在线观看视频网站免费| 国产成人aa在线观看| av网站免费在线观看视频 | 午夜视频国产福利| 亚洲欧洲日产国产| 亚洲av中文字字幕乱码综合| 国产精品日韩av在线免费观看| 成人欧美大片| 国产一区二区亚洲精品在线观看| 国产伦精品一区二区三区四那| 国产成人福利小说| 国产欧美日韩精品一区二区| 国产亚洲91精品色在线| 免费少妇av软件| 一级片'在线观看视频| 欧美97在线视频| 国产男人的电影天堂91| 日本一本二区三区精品| 国产精品人妻久久久影院| 久久久久久久久中文| 毛片一级片免费看久久久久| 在线观看av片永久免费下载| 亚洲人成网站在线播| 2018国产大陆天天弄谢| 国产成人freesex在线| 99久国产av精品国产电影| 国产午夜精品久久久久久一区二区三区| 日本三级黄在线观看| 国产精品日韩av在线免费观看| 欧美一区二区亚洲| 久久久久久久久久人人人人人人| 非洲黑人性xxxx精品又粗又长| 26uuu在线亚洲综合色| 久久精品熟女亚洲av麻豆精品 | 成人综合一区亚洲| 床上黄色一级片| 国产在线男女| 美女xxoo啪啪120秒动态图| 内射极品少妇av片p| 69av精品久久久久久| 国产男女超爽视频在线观看| 国产亚洲精品久久久com| 国产女主播在线喷水免费视频网站 | 亚洲av中文字字幕乱码综合| 我要看日韩黄色一级片| 日本熟妇午夜| 久久午夜福利片| 欧美xxⅹ黑人| 狂野欧美白嫩少妇大欣赏| 免费看a级黄色片| 欧美 日韩 精品 国产| eeuss影院久久| 亚洲精品一区蜜桃| 精品午夜福利在线看| 亚洲av中文字字幕乱码综合| 欧美一区二区亚洲| 少妇高潮的动态图| 丰满少妇做爰视频| 成人毛片60女人毛片免费| 午夜激情欧美在线| 亚洲精品日韩av片在线观看| 国产久久久一区二区三区| 国产精品一及| 亚洲av中文字字幕乱码综合| 一个人免费在线观看电影| 黄色日韩在线| 毛片女人毛片| 汤姆久久久久久久影院中文字幕 | 国产白丝娇喘喷水9色精品| 午夜福利在线观看吧| 久久久久久九九精品二区国产| 99久久人妻综合| 99九九线精品视频在线观看视频| 国产一级毛片七仙女欲春2| 蜜桃亚洲精品一区二区三区| 亚洲欧洲日产国产| 欧美日韩综合久久久久久| a级毛色黄片| 精品午夜福利在线看| 欧美不卡视频在线免费观看| 99热全是精品| 亚洲国产欧美在线一区| 三级经典国产精品| a级毛色黄片| 欧美精品国产亚洲| 国产精品1区2区在线观看.| 色视频www国产| 亚洲自偷自拍三级| 国产白丝娇喘喷水9色精品| 三级国产精品片| 国产男人的电影天堂91| 国产黄片美女视频| 日韩一区二区三区影片| 国产男人的电影天堂91| 亚洲成色77777| 亚洲一区高清亚洲精品| 精品久久国产蜜桃| 国产av码专区亚洲av| 欧美日韩亚洲高清精品| 97在线视频观看| av.在线天堂| 成年av动漫网址| 偷拍熟女少妇极品色| 欧美精品一区二区大全| 午夜久久久久精精品| 日日干狠狠操夜夜爽| 中文字幕av成人在线电影| 最近最新中文字幕大全电影3| 欧美zozozo另类| 只有这里有精品99|