• <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)展
    国产蜜桃级精品一区二区三区| 亚洲,欧美精品.| 18禁国产床啪视频网站| 1024香蕉在线观看| 丰满的人妻完整版| 国产高清激情床上av| 在线观看一区二区三区| 18禁黄网站禁片免费观看直播| 国产精品香港三级国产av潘金莲| 日本精品一区二区三区蜜桃| 18美女黄网站色大片免费观看| 国产成人av教育| 麻豆成人午夜福利视频| 在线观看免费午夜福利视频| 精品欧美国产一区二区三| 给我免费播放毛片高清在线观看| 久久久久久国产a免费观看| 一进一出好大好爽视频| 99久久精品热视频| 免费av毛片视频| 岛国视频午夜一区免费看| 免费高清视频大片| 成人特级av手机在线观看| 日本免费a在线| 亚洲欧美日韩卡通动漫| 麻豆av在线久日| 可以在线观看的亚洲视频| 欧美丝袜亚洲另类 | 日本黄色视频三级网站网址| 国产91精品成人一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美激情综合另类| 国产亚洲精品一区二区www| 成人av一区二区三区在线看| 国产欧美日韩精品亚洲av| 精品久久久久久久末码| 99久久成人亚洲精品观看| 日韩国内少妇激情av| 2021天堂中文幕一二区在线观| 19禁男女啪啪无遮挡网站| 午夜福利视频1000在线观看| 真实男女啪啪啪动态图| 久久精品91蜜桃| 丰满的人妻完整版| 丝袜人妻中文字幕| 黄色丝袜av网址大全| 少妇熟女aⅴ在线视频| 2021天堂中文幕一二区在线观| 淫秽高清视频在线观看| 亚洲欧美激情综合另类| 在线国产一区二区在线| 在线视频色国产色| 国产成人精品久久二区二区91| 欧美丝袜亚洲另类 | 欧美zozozo另类| 成年免费大片在线观看| 曰老女人黄片| 黄色女人牲交| 日本一二三区视频观看| 一级毛片女人18水好多| 欧美日韩黄片免| 亚洲18禁久久av| 黄片大片在线免费观看| 精品国产乱子伦一区二区三区| 亚洲av五月六月丁香网| 国产精品,欧美在线| 男女视频在线观看网站免费| 精品99又大又爽又粗少妇毛片 | 又大又爽又粗| 一级毛片女人18水好多| 午夜福利在线在线| 免费观看人在逋| 日本黄色片子视频| 国产精品99久久99久久久不卡| 中文字幕人成人乱码亚洲影| 这个男人来自地球电影免费观看| 久久国产精品人妻蜜桃| 女同久久另类99精品国产91| 亚洲激情在线av| 色精品久久人妻99蜜桃| 亚洲人成网站高清观看| 欧美日韩黄片免| 亚洲精品美女久久av网站| 亚洲色图 男人天堂 中文字幕| 国产在线精品亚洲第一网站| 大型黄色视频在线免费观看| 在线观看美女被高潮喷水网站 | 国产成人精品久久二区二区免费| 国产精品野战在线观看| 少妇熟女aⅴ在线视频| 国产一区二区三区视频了| 久久国产乱子伦精品免费另类| 精品国产美女av久久久久小说| 久久久久久大精品| 男女视频在线观看网站免费| 99久久国产精品久久久| 免费看日本二区| 免费在线观看日本一区| 国产成人啪精品午夜网站| 成人亚洲精品av一区二区| 亚洲最大成人中文| 一进一出好大好爽视频| 18禁黄网站禁片午夜丰满| 亚洲 欧美一区二区三区| 一进一出抽搐gif免费好疼| АⅤ资源中文在线天堂| 亚洲国产色片| 无限看片的www在线观看| 国产野战对白在线观看| 日日干狠狠操夜夜爽| 欧美另类亚洲清纯唯美| 日本成人三级电影网站| 视频区欧美日本亚洲| 亚洲欧美激情综合另类| 免费av不卡在线播放| 久久婷婷人人爽人人干人人爱| 最近最新免费中文字幕在线| 可以在线观看毛片的网站| 国产高清视频在线播放一区| a级毛片a级免费在线| 在线看三级毛片| 日本熟妇午夜| 日韩欧美三级三区| 国产亚洲欧美98| 少妇熟女aⅴ在线视频| 国产主播在线观看一区二区| 亚洲精品久久国产高清桃花| 国产成人精品久久二区二区免费| 天天躁狠狠躁夜夜躁狠狠躁| 欧美zozozo另类| 国产av一区在线观看免费| 女同久久另类99精品国产91| a级毛片a级免费在线| 成人鲁丝片一二三区免费| 岛国在线观看网站| 亚洲av成人av| 久久这里只有精品19| 亚洲av片天天在线观看| 热99在线观看视频| 成人永久免费在线观看视频| 午夜福利欧美成人| 好男人在线观看高清免费视频| 丁香六月欧美| 久久精品aⅴ一区二区三区四区| 国模一区二区三区四区视频 | 欧美性猛交黑人性爽| 欧美色欧美亚洲另类二区| 亚洲成av人片在线播放无| 久久国产乱子伦精品免费另类| 啪啪无遮挡十八禁网站| 午夜精品一区二区三区免费看| 久久精品人妻少妇| 夜夜爽天天搞| 啦啦啦韩国在线观看视频| 亚洲真实伦在线观看| 亚洲精品中文字幕一二三四区| 丰满人妻熟妇乱又伦精品不卡| 亚洲一区二区三区色噜噜| 久久香蕉国产精品| 99热这里只有精品一区 | 国产高清有码在线观看视频| 中文字幕久久专区| 亚洲成人中文字幕在线播放| 俺也久久电影网| 看免费av毛片| 中文字幕人成人乱码亚洲影| 在线观看免费午夜福利视频| 欧美日韩精品网址| 亚洲av日韩精品久久久久久密| 亚洲精品456在线播放app | 午夜福利视频1000在线观看| 亚洲专区字幕在线| 黄片小视频在线播放| 女警被强在线播放| 老司机深夜福利视频在线观看| 国产精品av视频在线免费观看| 午夜福利视频1000在线观看| 亚洲欧美日韩东京热| 国产高清激情床上av| 1024手机看黄色片| 757午夜福利合集在线观看| 国产精品,欧美在线| 国产成人av教育| 亚洲欧美日韩东京热| 精品乱码久久久久久99久播| 在线看三级毛片| 黄色日韩在线| 国产精品日韩av在线免费观看| 夜夜夜夜夜久久久久| 成年版毛片免费区| 亚洲精品国产精品久久久不卡| 日本黄色片子视频| 精品不卡国产一区二区三区| 黄色视频,在线免费观看| 亚洲国产精品久久男人天堂| 婷婷精品国产亚洲av在线| 免费搜索国产男女视频| 夜夜躁狠狠躁天天躁| 亚洲美女黄片视频| 欧美成人一区二区免费高清观看 | 免费观看人在逋| 国产人伦9x9x在线观看| 日本成人三级电影网站| 久久精品人妻少妇| 亚洲av中文字字幕乱码综合| 国内精品一区二区在线观看| 12—13女人毛片做爰片一| 国语自产精品视频在线第100页| 欧美成人一区二区免费高清观看 | 99热精品在线国产| netflix在线观看网站| 亚洲 国产 在线| 日韩三级视频一区二区三区| 1024手机看黄色片| 怎么达到女性高潮| 亚洲国产高清在线一区二区三| 午夜福利欧美成人| 国产精品久久久久久人妻精品电影| 日本黄色片子视频| 国产成人精品久久二区二区免费| 18禁观看日本| 国产高清视频在线观看网站| 久久久久久久久久黄片| 国产精品亚洲av一区麻豆| 黄色成人免费大全| 18禁美女被吸乳视频| 人妻丰满熟妇av一区二区三区| 日韩成人在线观看一区二区三区| 久久久国产成人免费| 少妇熟女aⅴ在线视频| 亚洲av成人av| 伦理电影免费视频| 色在线成人网| 别揉我奶头~嗯~啊~动态视频| 欧美精品啪啪一区二区三区| 成年女人永久免费观看视频| 国产综合懂色| 国产一区二区在线观看日韩 | 免费在线观看亚洲国产| 亚洲欧美精品综合一区二区三区| 久久中文字幕人妻熟女| 啦啦啦观看免费观看视频高清| 99re在线观看精品视频| 99精品久久久久人妻精品| 国产毛片a区久久久久| 免费高清视频大片| 国产精品久久电影中文字幕| 少妇人妻一区二区三区视频| 好男人在线观看高清免费视频| 国产亚洲精品一区二区www| 亚洲美女黄片视频| 国产视频一区二区在线看| 亚洲av五月六月丁香网| www日本黄色视频网| 最近最新中文字幕大全免费视频| 国产成人av教育| 黄色丝袜av网址大全| 国产伦一二天堂av在线观看| 国产高清视频在线播放一区| 亚洲av免费在线观看| 亚洲黑人精品在线| 成人18禁在线播放| 精品不卡国产一区二区三区| 亚洲,欧美精品.| 好男人在线观看高清免费视频| 免费观看精品视频网站| 黄色 视频免费看| 久久精品影院6| 欧美日韩一级在线毛片| 成年人黄色毛片网站| 免费在线观看亚洲国产| 国产一区在线观看成人免费| 婷婷精品国产亚洲av| 日韩欧美在线二视频| 在线观看免费午夜福利视频| 十八禁人妻一区二区| 91av网一区二区| 美女黄网站色视频| 99国产精品一区二区蜜桃av| 成年人黄色毛片网站| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久国内视频| 久久香蕉国产精品| 热99在线观看视频| 此物有八面人人有两片| 亚洲欧美日韩高清在线视频| 欧美不卡视频在线免费观看| 色噜噜av男人的天堂激情| 嫩草影院入口| 国产成人av激情在线播放| 日本黄色视频三级网站网址| 黄色 视频免费看| 亚洲专区国产一区二区| 亚洲欧美日韩东京热| 香蕉久久夜色| 亚洲国产日韩欧美精品在线观看 | 色视频www国产| 美女高潮的动态| 黄色丝袜av网址大全| 精品久久久久久久久久久久久| 欧美另类亚洲清纯唯美| 色综合欧美亚洲国产小说| 精品久久蜜臀av无| 国产v大片淫在线免费观看| 首页视频小说图片口味搜索| 久久精品国产99精品国产亚洲性色| 免费看日本二区| 欧洲精品卡2卡3卡4卡5卡区| 国产人伦9x9x在线观看| 看片在线看免费视频| 成人永久免费在线观看视频| 欧美精品啪啪一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 成年版毛片免费区| 熟女电影av网| 999久久久精品免费观看国产| 国产精品综合久久久久久久免费| 在线a可以看的网站| 亚洲人成网站在线播放欧美日韩| 国产精品免费一区二区三区在线| 悠悠久久av| 三级男女做爰猛烈吃奶摸视频| 久久久久免费精品人妻一区二区| 欧美最黄视频在线播放免费| 精品日产1卡2卡| 性欧美人与动物交配| 老司机午夜福利在线观看视频| 国产单亲对白刺激| 老司机福利观看| 亚洲第一欧美日韩一区二区三区| 91麻豆av在线| 欧洲精品卡2卡3卡4卡5卡区| 国产高清videossex| 高清在线国产一区| 看片在线看免费视频| 亚洲色图av天堂| 嫩草影视91久久| 男女下面进入的视频免费午夜| 国产成人福利小说| 麻豆av在线久日| 成年女人毛片免费观看观看9| 黑人欧美特级aaaaaa片| 天天一区二区日本电影三级| 又粗又爽又猛毛片免费看| 欧美午夜高清在线| 岛国在线免费视频观看| 日本撒尿小便嘘嘘汇集6| 欧美国产日韩亚洲一区| 岛国在线观看网站| 91av网一区二区| 极品教师在线免费播放| 狂野欧美白嫩少妇大欣赏| 制服人妻中文乱码| 不卡av一区二区三区| 老熟妇仑乱视频hdxx| 亚洲熟妇熟女久久| 久久午夜亚洲精品久久| 特大巨黑吊av在线直播| 久久久久久久久中文| 亚洲专区国产一区二区| 国产精品影院久久| 欧美性猛交黑人性爽| 美女 人体艺术 gogo| 国产精品 欧美亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 村上凉子中文字幕在线| 偷拍熟女少妇极品色| 成年免费大片在线观看| 中文字幕高清在线视频| 日韩人妻高清精品专区| 国产亚洲欧美在线一区二区| 国产成人av教育| 亚洲在线观看片| 日本成人三级电影网站| 国产一区二区在线观看日韩 | 欧美成人性av电影在线观看| 欧美成狂野欧美在线观看| 免费一级毛片在线播放高清视频| 国产av在哪里看| 欧美在线一区亚洲| 国产日本99.免费观看| 免费观看精品视频网站| 又紧又爽又黄一区二区| 美女大奶头视频| 制服丝袜大香蕉在线| 最新在线观看一区二区三区| 美女被艹到高潮喷水动态| 一进一出抽搐gif免费好疼| 免费电影在线观看免费观看| 日日干狠狠操夜夜爽| 色综合欧美亚洲国产小说| 国产精品一区二区三区四区免费观看 | 老司机深夜福利视频在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲av电影在线进入| 国产精品亚洲av一区麻豆| 日韩国内少妇激情av| 天天躁狠狠躁夜夜躁狠狠躁| 听说在线观看完整版免费高清| 国产亚洲欧美98| 香蕉久久夜色| 国产精品久久久av美女十八| 两个人视频免费观看高清| 色视频www国产| 99国产精品一区二区三区| 日韩欧美精品v在线| 啪啪无遮挡十八禁网站| 级片在线观看| 国产午夜福利久久久久久| 国产成人福利小说| 亚洲欧美激情综合另类| 亚洲自偷自拍图片 自拍| 亚洲自拍偷在线| 欧美成人性av电影在线观看| 精品久久久久久久久久久久久| 夜夜看夜夜爽夜夜摸| 国产成年人精品一区二区| 日韩有码中文字幕| netflix在线观看网站| 97人妻精品一区二区三区麻豆| 99国产精品一区二区蜜桃av| 国产成人啪精品午夜网站| 久久国产精品影院| 国产精品久久久久久精品电影| 欧美xxxx黑人xx丫x性爽| 最近最新中文字幕大全电影3| 一区二区三区激情视频| 亚洲熟妇熟女久久| 美女 人体艺术 gogo| 老司机深夜福利视频在线观看| 日韩欧美精品v在线| 日日干狠狠操夜夜爽| 久久草成人影院| 国产精品野战在线观看| 国产极品精品免费视频能看的| 丰满人妻一区二区三区视频av | 真人做人爱边吃奶动态| 久久久久久人人人人人| 亚洲精品美女久久久久99蜜臀| 给我免费播放毛片高清在线观看| 亚洲精品乱码久久久v下载方式 | 欧美不卡视频在线免费观看| 亚洲 欧美 日韩 在线 免费| 热99在线观看视频| 久久中文字幕一级| 精品一区二区三区四区五区乱码| 欧美最黄视频在线播放免费| 国内精品美女久久久久久| 国产一区二区激情短视频| 精品一区二区三区四区五区乱码| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 一个人免费在线观看电影 | 免费无遮挡裸体视频| 免费av毛片视频| 国产欧美日韩一区二区精品| 精品国内亚洲2022精品成人| 岛国视频午夜一区免费看| 岛国在线免费视频观看| 天堂影院成人在线观看| 亚洲中文字幕一区二区三区有码在线看 | 精品一区二区三区视频在线观看免费| 日韩三级视频一区二区三区| 国产精品亚洲av一区麻豆| 欧美激情久久久久久爽电影| 国产一区在线观看成人免费| 禁无遮挡网站| 国产伦一二天堂av在线观看| 亚洲,欧美精品.| 成人鲁丝片一二三区免费| 麻豆成人av在线观看| 日韩三级视频一区二区三区| 特大巨黑吊av在线直播| 全区人妻精品视频| 露出奶头的视频| 熟女少妇亚洲综合色aaa.| 欧美高清成人免费视频www| 国产精品av视频在线免费观看| 99久久国产精品久久久| 亚洲狠狠婷婷综合久久图片| 无人区码免费观看不卡| 欧美激情久久久久久爽电影| 亚洲精品国产精品久久久不卡| 中出人妻视频一区二区| 欧美不卡视频在线免费观看| 久久久久久久午夜电影| 91老司机精品| 午夜免费观看网址| 欧美成狂野欧美在线观看| 91老司机精品| 欧美日本亚洲视频在线播放| 好男人在线观看高清免费视频| 欧美日韩瑟瑟在线播放| 十八禁网站免费在线| 国产精品 欧美亚洲| 日韩精品青青久久久久久| 99riav亚洲国产免费| 又黄又粗又硬又大视频| 19禁男女啪啪无遮挡网站| 国产成人一区二区三区免费视频网站| 亚洲av片天天在线观看| 在线免费观看的www视频| 精品午夜福利视频在线观看一区| 国产高清视频在线播放一区| 亚洲aⅴ乱码一区二区在线播放| 日日摸夜夜添夜夜添小说| 国产免费av片在线观看野外av| 免费在线观看成人毛片| 小说图片视频综合网站| 欧美日韩乱码在线| 国产av麻豆久久久久久久| 热99在线观看视频| 两人在一起打扑克的视频| 中文亚洲av片在线观看爽| 两个人视频免费观看高清| 真人做人爱边吃奶动态| 亚洲激情在线av| 中出人妻视频一区二区| 又粗又爽又猛毛片免费看| 男人的好看免费观看在线视频| 在线十欧美十亚洲十日本专区| 亚洲国产欧美人成| 亚洲国产中文字幕在线视频| 国产乱人视频| 亚洲成av人片免费观看| 欧美乱码精品一区二区三区| 久久久水蜜桃国产精品网| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人精品中文字幕电影| 国产精品美女特级片免费视频播放器 | 午夜日韩欧美国产| 亚洲激情在线av| 美女扒开内裤让男人捅视频| 久久人人精品亚洲av| 午夜精品久久久久久毛片777| 国产黄片美女视频| 我要搜黄色片| 国产成人精品久久二区二区免费| 日韩欧美三级三区| 亚洲国产中文字幕在线视频| 最好的美女福利视频网| 高清毛片免费观看视频网站| 亚洲欧美日韩东京热| 久久久久久国产a免费观看| 99热这里只有是精品50| 九九久久精品国产亚洲av麻豆 | 色尼玛亚洲综合影院| 午夜两性在线视频| 啦啦啦免费观看视频1| 色播亚洲综合网| 久久精品国产99精品国产亚洲性色| 少妇的逼水好多| 午夜视频精品福利| 久久中文字幕人妻熟女| 精品日产1卡2卡| 男人舔女人的私密视频| 最近在线观看免费完整版| 老熟妇乱子伦视频在线观看| 亚洲精品美女久久av网站| 久久伊人香网站| 97人妻精品一区二区三区麻豆| 老司机午夜十八禁免费视频| 999久久久国产精品视频| 亚洲欧美日韩高清专用| 亚洲中文字幕日韩| 母亲3免费完整高清在线观看| 老司机深夜福利视频在线观看| a级毛片a级免费在线| 精品国内亚洲2022精品成人| 色吧在线观看| 天堂影院成人在线观看| 国产高清激情床上av| 欧美3d第一页| 免费看a级黄色片| 久久久久性生活片| 日韩人妻高清精品专区| 国产欧美日韩一区二区精品| 国产99白浆流出| 禁无遮挡网站| 亚洲 欧美 日韩 在线 免费| 中国美女看黄片| 亚洲精品一卡2卡三卡4卡5卡| 一个人免费在线观看电影 | 国产亚洲av高清不卡| 国产黄色小视频在线观看| 欧美绝顶高潮抽搐喷水| 国产伦人伦偷精品视频| 99精品在免费线老司机午夜| 亚洲自拍偷在线| 禁无遮挡网站| 动漫黄色视频在线观看| 99国产综合亚洲精品| 欧美国产日韩亚洲一区| 国产不卡一卡二| 久久草成人影院| 在线观看一区二区三区| 少妇的丰满在线观看| 欧美av亚洲av综合av国产av| 亚洲第一欧美日韩一区二区三区| 欧美av亚洲av综合av国产av| 国产97色在线日韩免费| 村上凉子中文字幕在线| 欧美色欧美亚洲另类二区| 日韩欧美国产一区二区入口| 午夜激情欧美在线| 男人舔奶头视频| 欧美丝袜亚洲另类 | 此物有八面人人有两片| 久久久久性生活片| 亚洲乱码一区二区免费版| 国产v大片淫在线免费观看| 美女免费视频网站| 法律面前人人平等表现在哪些方面| 男人舔女人下体高潮全视频| 高清在线国产一区|