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

    微塑料對(duì)結(jié)冰過(guò)程中環(huán)境因子遷移的影響

    2022-12-19 12:53:12王志超竇雅嬌康延秋楊文煥敬雙怡李衛(wèi)平
    中國(guó)環(huán)境科學(xué) 2022年11期
    關(guān)鍵詞:冰體結(jié)冰塑料

    王志超,竇雅嬌,康延秋,周 鑫,楊文煥,敬雙怡,李衛(wèi)平

    微塑料對(duì)結(jié)冰過(guò)程中環(huán)境因子遷移的影響

    王志超,竇雅嬌,康延秋,周 鑫,楊文煥,敬雙怡,李衛(wèi)平*

    (內(nèi)蒙古科技大學(xué)能源與環(huán)境學(xué)院,內(nèi)蒙古 包頭 014010)

    為探究結(jié)冰過(guò)程中微塑料賦存對(duì)典型環(huán)境因子(總氮、總磷、鹽度、化學(xué)需氧量、懸浮物)等分布及遷移規(guī)律的影響,采用室內(nèi)模擬的方式研究不同條件下(結(jié)冰比例、結(jié)冰溫度及方式、初始濃度)各典型環(huán)境因子在冰-水相間的分布及遷移規(guī)律,采用物質(zhì)分配系數(shù)()表征環(huán)境因子的遷移能力.結(jié)果表明,結(jié)冰過(guò)程中冰體對(duì)環(huán)境因子具有排斥作用,不同條件下微塑料賦存對(duì)環(huán)境因子的分布規(guī)律均有一定程度的影響,由于微塑料自身特性致使原本應(yīng)遷移至冰下水體中的環(huán)境因子部分滯留于冰體中,導(dǎo)致冰體中環(huán)境因子濃度的提高,為對(duì)照組冰體中環(huán)境因子濃度的1.13~1.49倍;同時(shí)使冰下水體中環(huán)境因子濃度下降,為對(duì)照組冰下水體中環(huán)境因子濃度的0.73~0.93倍;微塑料的賦存同時(shí)導(dǎo)致環(huán)境因子分配系數(shù)值提高0.04~0.18,致使環(huán)境因子向冰下水體遷移的能力下降,但微塑料的賦存并未改變結(jié)冰過(guò)程中環(huán)境因子由冰體向冰下水體遷移的趨勢(shì),微塑料對(duì)于環(huán)境因子的攜帶作用小于冰體對(duì)環(huán)境因子的排斥作用;同時(shí)不同條件下結(jié)冰過(guò)程中微塑料對(duì)環(huán)境因子分布與遷移機(jī)理的影響均可以從結(jié)晶學(xué)角度與共晶理論得到解釋.

    結(jié)冰過(guò)程;微塑料;環(huán)境因子;遷移規(guī)律;作用機(jī)理

    冰封是全球中高緯度地區(qū)地表水存在的一種自然現(xiàn)象及水文特征,全球大約5億個(gè)湖泊會(huì)在冬季時(shí)結(jié)冰[1].而我國(guó)地處北半球,絕大部分地區(qū)的湖泊均有不同程度的冰情出現(xiàn),特別是黃河以北、青藏高原等地區(qū)的湖泊,每年冰封時(shí)間長(zhǎng)達(dá)4~6個(gè)月[2],冰封現(xiàn)象影響范圍大,持續(xù)時(shí)間長(zhǎng),已成為我國(guó)特別是我國(guó)北方地區(qū)水體的重要特征.以位于內(nèi)蒙古高原的典型寒區(qū)湖泊岱海湖為例,冬季嚴(yán)寒漫長(zhǎng),全年平均氣溫僅5℃,每年11月左右進(jìn)入冰封期,一直持續(xù)到次年4月開(kāi)始融化,年冰封期可達(dá)6個(gè)月[3].

    冰蓋為覆蓋某水域全部或部分的冰,一般常見(jiàn)于高原地區(qū).湖泊結(jié)冰初期湖體表面先由冰晶進(jìn)一步形成冰核,冰核在水平及垂直方向逐漸生長(zhǎng),最終宏觀表現(xiàn)為冰蓋的形成[4].因此冰封期內(nèi)冰蓋的存在不僅會(huì)降低光穿透率、減弱外界與湖體間的能量及其物質(zhì)交換等[5],也使得多數(shù)污染物由于冷凍濃縮效應(yīng)等在結(jié)冰時(shí)不斷被排入水體中,進(jìn)而打破冰-水相間濃度平衡,使得污染物均勻的擴(kuò)散到整個(gè)冰下水體,最終導(dǎo)致冰下水體污染物濃度明顯升高,水質(zhì)進(jìn)一步惡化;研究發(fā)現(xiàn),受冰封作用的影響,有機(jī)農(nóng)藥、重金屬、營(yíng)養(yǎng)鹽等常規(guī)污染物在自然水體中的濃度均表現(xiàn)出冰少水多的特性,如烏梁素海冰封時(shí)不同形態(tài)的氮、磷從冰體遷移到冰下水體中,使水體中各項(xiàng)指標(biāo)均大于冰體,以總氮(TN)為例,冰下水體中TN濃度就為冰體的3.16倍[6];但新型污染物微塑料在冰-水相間的遷移規(guī)律卻與之相反,呈現(xiàn)出冰多水少的特點(diǎn),在對(duì)北極海冰的研究中發(fā)現(xiàn)冰樣中微塑料豐度高于水體幾個(gè)數(shù)量級(jí)[7],同時(shí)研究發(fā)現(xiàn)我國(guó)烏梁素海湖冰中微塑料豐度為水體的10~100倍[8].另一方面,微塑料因比表面積大、吸附污染能力強(qiáng)等特性,已成為水體中TN、總磷(TP)等湖泊污染物的重要載體,進(jìn)而與污染物表現(xiàn)出相互影響的關(guān)系,已有研究表明[9]冰封作用下鹽度與微塑料之間存在顯著的正相關(guān)性.

    目前對(duì)于結(jié)冰過(guò)程中環(huán)境因子的遷移影響多數(shù)集中于重金屬、有機(jī)農(nóng)藥等傳統(tǒng)污染物且研究指標(biāo)單一,對(duì)于微塑料大量賦存情況下對(duì)冰-水環(huán)境中典型污染物的遷移效應(yīng)研究較少.針對(duì)此特性,本文前期曾對(duì)我國(guó)北方寒區(qū)湖泊岱海做過(guò)相關(guān)研究[3],發(fā)現(xiàn)冰蓋中微塑料豐度約為同一采樣點(diǎn)水體中微塑料豐度的4~9倍,其中微塑料豐度與岱海典型環(huán)境因子存在一定相關(guān)性,但野外實(shí)驗(yàn)易受水生生物或底泥沉積物物質(zhì)交換等各種因素的影響,難以定性定量得出結(jié)冰過(guò)程中微塑料對(duì)環(huán)境因子分布及遷移規(guī)律的影響,因此較多學(xué)者選用室內(nèi)模擬實(shí)驗(yàn)的方法避免不足.呂宏洲等[10]為研究冰封期烏梁素海中單一因素對(duì)營(yíng)養(yǎng)鹽等遷移的影響進(jìn)行了室內(nèi)結(jié)冰模擬實(shí)驗(yàn),研究表明單一因素不同,營(yíng)養(yǎng)鹽等在冰-水相中的分布、遷移特征不同,并且室內(nèi)模擬實(shí)驗(yàn)的結(jié)果與烏梁素海冰封期實(shí)際情況較為接近.本研究針對(duì)前期野外實(shí)驗(yàn)結(jié)果[3]及岱海冰封時(shí)間長(zhǎng)、微塑料豐度大等特性為背景開(kāi)展室內(nèi)模擬實(shí)驗(yàn),定性并定量評(píng)估結(jié)冰過(guò)程中微塑料賦存對(duì)湖泊冰-水環(huán)境中典型環(huán)境因子的影響,以期為冰封作用下微塑料對(duì)環(huán)境因子的遷移影響提供參考,為微塑料賦存條件下寒區(qū)湖泊季節(jié)性污染研究奠定理論基礎(chǔ).

    1 材料與方法

    1.1 實(shí)驗(yàn)裝置與儀器

    因天然水體的結(jié)冰過(guò)程是自上而下的,為了使室內(nèi)模擬實(shí)驗(yàn)與自然水體的結(jié)冰效果相似,本研究采用一種自制敞口單向?qū)Ю浣Y(jié)冰模擬器(如圖1(a)所示)[11],實(shí)驗(yàn)裝置采用兩組圓柱形不銹鋼儀器組合而成(內(nèi)圓柱20cm、50cm;外圓柱30cm、60cm),在內(nèi)、外圓柱間利用石棉保溫材料包裹及填充,阻斷桶壁和桶底與外界的熱量傳遞,以確保能量自頂而下傳遞,進(jìn)而模擬野外湖泊從表面開(kāi)始的結(jié)冰過(guò)程,其中內(nèi)圓柱可獨(dú)立取出以便結(jié)冰后冰柱的取出及后續(xù)的清洗工作,該套裝置可一定程度上避免野外工作時(shí)由于外源水體污染物等干擾所帶來(lái)的誤差.

    圖1 實(shí)驗(yàn)裝置與儀器

    本實(shí)驗(yàn)采用自制電阻絲測(cè)厚儀(圖2(b))測(cè)量結(jié)冰比例[12],實(shí)驗(yàn)開(kāi)始前將總長(zhǎng)度為的鎳鉻電阻絲提前放入裝置內(nèi),利用12V蓄電池對(duì)電阻絲加熱直至融化其周圍冰體,后上拉電阻絲使下?lián)醢謇帘w下表面,此時(shí)所測(cè)量的冰體上部電阻絲長(zhǎng)度為1,從而可得結(jié)冰厚度2=-1.結(jié)冰比例可通過(guò)冰厚與內(nèi)圓柱高度之比求得.整個(gè)結(jié)冰過(guò)程在設(shè)有溫控的試驗(yàn)箱中完成,低溫試驗(yàn)箱溫度最低可達(dá)-30℃,溫度上下波動(dòng)在±1℃.

    1.2 檢測(cè)指標(biāo)及方法

    本研究選取TN、TP、鹽度、COD、懸浮物(SS)為典型環(huán)境因子.為排除天然湖泊復(fù)雜的水質(zhì)及水動(dòng)力條件對(duì)研究的影響,實(shí)驗(yàn)選用超純水及各物質(zhì)標(biāo)準(zhǔn)藥品進(jìn)行標(biāo)準(zhǔn)溶液的配置,鹽度的配置采用《海洋調(diào)查規(guī)范》[13]中的配方;配制SS所用泥樣采自岱海湖底泥,取回后放置于60℃烘箱中烘干至恒重,冷卻研磨后使用200目不銹鋼篩子過(guò)濾并置于干燥器保存[14];其余環(huán)境因子的配置方法按照《水和廢水監(jiān)測(cè)分析方法(第四版)》[15]配置.樣品中環(huán)境因子TN、TP、鹽度、COD、SS的檢測(cè)方法均按國(guó)家標(biāo)準(zhǔn)方法及相關(guān)文獻(xiàn)測(cè)定[16-19].

    1.3 實(shí)驗(yàn)設(shè)計(jì)

    以2021年結(jié)冰前岱海湖中5種環(huán)境因子及微塑料實(shí)測(cè)的質(zhì)量濃度及豐度為依據(jù),設(shè)定TN、TP、鹽度、COD、SS質(zhì)量濃度和微塑料豐度分別為4.14mg/L、0.15mg/L、14.84′10-12、107mg/L、32.72mg/L、280n/L;其中微塑料類型采用岱海湖常見(jiàn)的高密度塑料顆粒PVC與低密度塑料顆粒PP,微塑料粒徑經(jīng)岱海湖實(shí)測(cè)后選取粒徑占比最大的950mm.

    冰核形成與冰體生長(zhǎng)階段為結(jié)冰的兩個(gè)主要過(guò)程,結(jié)冰比例、初始結(jié)冰溫度、溶液初始濃度都是影響結(jié)冰過(guò)程的重要因素[20],為研究結(jié)冰過(guò)程中微塑料賦存對(duì)岱海典型環(huán)境因子遷移規(guī)律的影響,采用不同結(jié)冰比例下(即結(jié)冰厚度,不同的結(jié)冰厚度決定著冰晶體捕獲雜質(zhì)含量的多少)、不同結(jié)冰溫度及方式下(即初始溶液所處環(huán)境的溫度,其決定著溶液的冷卻速度)、不同初始濃度下(即配置溶液時(shí)所設(shè)置的各環(huán)境因子的質(zhì)量濃度大小,不同的溶液初始濃度影響著冰的生長(zhǎng)速度)實(shí)驗(yàn)組與對(duì)照組的方式研究微塑料對(duì)環(huán)境因子分布及遷移規(guī)律的影響,實(shí)驗(yàn)組中使用超純水配置環(huán)境因子與微塑料混合溶液,對(duì)照組使用超純水僅配置環(huán)境因子溶液,其余變量與實(shí)驗(yàn)組均一致,所有實(shí)驗(yàn)均設(shè)置3組平行,檢測(cè)誤差均控制在5%以內(nèi),具體實(shí)驗(yàn)設(shè)置如下:

    1.3.1 不同結(jié)冰比例下微塑料對(duì)環(huán)境因子遷移規(guī)律的影響 使用超純水配置如上溶液后放入結(jié)冰模擬器中,在-15℃低溫實(shí)驗(yàn)箱中冷凍,經(jīng)本團(tuán)隊(duì)長(zhǎng)期對(duì)冰封期岱海湖檢測(cè)發(fā)現(xiàn)其結(jié)冰比例為8%~20%,故將溶液冷凍至結(jié)冰比例分別為8%、14%、20%時(shí)取出,實(shí)驗(yàn)完成后冰樣采用鋁箔紙密封并避光融化,融化后與冰下水樣共同進(jìn)行環(huán)境因子質(zhì)量濃度的檢測(cè).

    1.3.2 不同結(jié)冰溫度及方式下微塑料對(duì)環(huán)境因子遷移規(guī)律的影響 配置與上述濃度相同的溶液并放入結(jié)冰模擬器中,在-5,-15,-25℃下放入低溫實(shí)驗(yàn)箱冷凍,除此之外,為探究自然環(huán)境中溫度晝夜變化對(duì)結(jié)冰過(guò)程產(chǎn)生的影響,基于冬季岱海實(shí)況進(jìn)行-5℃與-25℃時(shí)間等長(zhǎng)變溫實(shí)驗(yàn),4組實(shí)驗(yàn)結(jié)冰比例為14%時(shí)取出,冰樣采用鋁箔紙密封并避光融化,融化后與冰下水樣共同進(jìn)行環(huán)境因子質(zhì)量濃度的檢測(cè).

    1.3.3 不同初始濃度下微塑料對(duì)環(huán)境因子遷移規(guī)律的影響 配置3種不同濃度分別為上述濃度及其上下變化30%的溶液并放入結(jié)冰模擬器中,濃度由低至高變化,實(shí)驗(yàn)分別命名為C1、C2、C3,在-15℃的低溫實(shí)驗(yàn)箱中冷凍,結(jié)冰比例為14%時(shí)取出,冰樣采用鋁箔紙密封并避光融化,融化后與冰下水樣共同進(jìn)行環(huán)境因子質(zhì)量濃度的檢測(cè).

    系統(tǒng)效應(yīng)不同于個(gè)體效應(yīng)。個(gè)體效應(yīng)是指每一個(gè)具體標(biāo)準(zhǔn)都有其特定的功能,在標(biāo)準(zhǔn)實(shí)施中產(chǎn)生特定的效應(yīng)。系統(tǒng)效應(yīng)則是指由若干個(gè)具有內(nèi)在聯(lián)系的標(biāo)準(zhǔn)個(gè)體組成的標(biāo)準(zhǔn)系統(tǒng)具有的特定功能,在實(shí)施中產(chǎn)生的特定效應(yīng)。但這種效應(yīng)絕不是互不相干的個(gè)體效應(yīng)的簡(jiǎn)單相加,而是從要素量的集合達(dá)到整體質(zhì)的飛躍中產(chǎn)生的,是相關(guān)標(biāo)準(zhǔn)之間相互作用產(chǎn)生的相干效應(yīng),比各個(gè)標(biāo)準(zhǔn)效應(yīng)的簡(jiǎn)單總和大得多。[2]

    1.4 數(shù)據(jù)處理與分析方法

    數(shù)據(jù)處理利用Microsoft Excel2018完成,圖片的繪制利用Origin2019、Soliworks及Microsoft Office PowerPoint2018完成.

    用分配系數(shù)表征結(jié)冰過(guò)程中各環(huán)境因子的排出效應(yīng),即:

    =i/w(1)

    式中:i為冰體中環(huán)境因子濃度,mg/L或10-12;w為對(duì)應(yīng)冰下水體中環(huán)境因子濃度,mg/L或10-12.

    2 結(jié)果與討論

    2.1 不同結(jié)冰比例下微塑料賦存對(duì)環(huán)境因子遷移的影響

    如圖2所示,無(wú)論實(shí)驗(yàn)中是否存在微塑料,不同結(jié)冰比例下實(shí)驗(yàn)組與對(duì)照組中各環(huán)境因子的分布均呈現(xiàn)出原水樣濃度高于冰體中濃度,小于冰下水體中濃度的變化規(guī)律,即冰樣中濃度<原水樣濃度<冰下水樣濃度.實(shí)驗(yàn)組中各環(huán)境因子在冰體中的濃度是對(duì)照組冰體中濃度的1.19~1.39倍,而各環(huán)境因子在冰下水體中的濃度則是對(duì)照組冰下水體濃度的0.77~0.93倍,表明微塑料的賦存致使原本應(yīng)遷移至冰下水體中的環(huán)境因子而滯留于冰體中,其原因可能是微塑料與大多數(shù)污染物遷移規(guī)律不同,微塑料在冰體中的豐度高于相應(yīng)水體中微塑料的豐度,表現(xiàn)出冰多水少的特點(diǎn)[3],而環(huán)境中的微塑料又因其比表面積較小等特點(diǎn)表現(xiàn)出易吸附環(huán)境中其他物質(zhì)的特征,因此實(shí)驗(yàn)組中冰體內(nèi)部由于微塑料的大量賦存導(dǎo)致一定程度上攜帶其他物質(zhì),致使實(shí)驗(yàn)組中冰體內(nèi)部環(huán)境因子濃度的提高,故表現(xiàn)出實(shí)驗(yàn)組與對(duì)照組中各環(huán)境因子分布規(guī)律的不同.

    隨著結(jié)冰比例的增加實(shí)驗(yàn)組與對(duì)照組中環(huán)境因子的分配系數(shù)值均呈現(xiàn)下降趨勢(shì),各物質(zhì)在結(jié)冰比例為20%時(shí)值達(dá)到最小值,此時(shí)各物質(zhì)在冰體中的濃度逐漸減少,冰下水體中的濃度逐漸增大,表明水體中無(wú)論是否有微塑料的存在,隨著結(jié)冰比例的升高各環(huán)境因子向冰下水體的遷移能力提高.但由于微塑料的賦存,在不同的結(jié)冰比例下實(shí)驗(yàn)組中各環(huán)境因子值較對(duì)照組中值提高了0.04~ 0.17,表明微塑料的賦存導(dǎo)致環(huán)境因子向冰下水體遷移的能力下降,其中結(jié)冰比例為8%時(shí)實(shí)驗(yàn)組與對(duì)照組中值差距為0.07~0.17,表明結(jié)冰比例越低各環(huán)境因子在冰-水間的遷移能力受微塑料賦存影響最大,因此這一現(xiàn)象出現(xiàn)的原因可能是結(jié)冰初期冰體所俘獲的物質(zhì)由于冷凍濃縮效應(yīng)大量通過(guò)遷移通道向冰下水體遷移,而微塑料的存在導(dǎo)致一部分的環(huán)境因子與微塑料結(jié)合,減少各環(huán)境因子在冰-水間的遷移量;同時(shí)作為固體污染物的微塑料會(huì)占據(jù)冰體內(nèi)部的晶體位置,致使冰體內(nèi)部的遷移通道減少[20];此外雖然結(jié)冰比例的提高致使冰體中微塑料的豐度相對(duì)增加,但由于微塑料豐度與冰厚呈現(xiàn)出負(fù)相關(guān),因此新生成冰體中微塑料的增加量較之前減少,則對(duì)于環(huán)境因子遷移影響變小,而出現(xiàn)結(jié)冰初期微塑料對(duì)各環(huán)境因子的遷移能力影響最大這一現(xiàn)象.

    除此之外,實(shí)驗(yàn)組中各環(huán)境因子值的提高并不相同,表明微塑料吸附各環(huán)境因子的能力不同,之間的相關(guān)關(guān)系表現(xiàn)出差異,袁海英等[21]在對(duì)滇池中微塑料污染與富營(yíng)養(yǎng)化的相關(guān)性研究中發(fā)現(xiàn),微塑料豐度與TN表現(xiàn)出極顯著正相關(guān)關(guān)系,而與TP等指標(biāo)間并未有顯著關(guān)系出現(xiàn).

    綜上所述,不同結(jié)冰比例下微塑料的賦存僅能一定程度上改變環(huán)境因子在冰-水相的分布規(guī)律,而環(huán)境因子的遷移規(guī)律并未發(fā)生改變,在結(jié)冰過(guò)程中仍然表現(xiàn)出由冰體向冰下水體遷移的規(guī)律,表明微塑料對(duì)于環(huán)境因子的攜帶作用小于冰體對(duì)環(huán)境因子的排斥作用,不足以改變物質(zhì)遷移的總體趨勢(shì).

    2.2 不同結(jié)冰溫度及方式下微塑料賦存對(duì)環(huán)境因子遷移的影響

    在微塑料賦存的情況下冰體中各環(huán)境因子濃度為對(duì)照組的1.13~1.41倍,實(shí)驗(yàn)組冰下水體中環(huán)境因子濃度為對(duì)照組的0.80~0.93倍,在不同結(jié)冰溫度條件下微塑料的賦存仍然對(duì)環(huán)境因子濃度分布產(chǎn)生影響,隨著溫度的降低各環(huán)境因子冰體與冰下水體中物質(zhì)濃度相差逐漸減小;其中變溫結(jié)冰方式同樣會(huì)對(duì)環(huán)境因子的分布產(chǎn)生影響,變溫結(jié)冰方式和恒溫結(jié)冰方式冰體中環(huán)境因子濃度關(guān)系為:-25℃的恒溫結(jié)冰方式>-15℃的恒溫結(jié)冰方式>-5℃與-25℃時(shí)間等長(zhǎng)變溫結(jié)冰方式>-5℃的恒溫結(jié)冰方式,時(shí)間等長(zhǎng)變溫的方式下相較于-15℃、-25℃的恒溫結(jié)冰方式,結(jié)冰速率較低,各物質(zhì)更易遷移到冰下水體中,因此時(shí)間等長(zhǎng)變溫結(jié)冰方式下冰體中物質(zhì)的含量更少.

    不同結(jié)冰溫度及方式下實(shí)驗(yàn)組中環(huán)境因子的值相較于對(duì)照組值升高范圍為0.05~0.18,各環(huán)境因子的遷移規(guī)律受微塑料賦存的影響大小不同,其中TP的遷移規(guī)律受微塑料影響值下降范圍最大為0.12~0.18.這是由于結(jié)冰溫度是決定結(jié)冰速率的因素之一,溫度越低,外界與冰體之間的熱交換量變多,冰體中產(chǎn)生更多的冰體顆粒及晶核數(shù)等,冰體的結(jié)冰速率加快,從而提高了水分子向冰-水界面移動(dòng)的速度,當(dāng)物質(zhì)向冰下水體移動(dòng)的速度小于水分子移動(dòng)的速度時(shí),溶液中的物質(zhì)來(lái)不及“逃逸”,進(jìn)而被俘獲在冰體中,從而導(dǎo)致冰體中物質(zhì)的濃度增加,而微塑料的賦存使一部分的環(huán)境因子聚集在微塑料上,導(dǎo)致其受冰體的排斥作用減弱,向冰下水體的移動(dòng)速度下降,更多的滯留于冰體中[22].

    另一方面結(jié)冰溫度越低,冰體需要更多的面積來(lái)釋放熱量,因此冰晶體的生長(zhǎng)變?yōu)闃?shù)枝狀以增大表面積,故冰體產(chǎn)生樹(shù)枝狀的分支進(jìn)一步捕獲物質(zhì)[23],使冰體的純度降低,環(huán)境因子的濃度增大,則對(duì)物質(zhì)的排斥作用就越低,實(shí)驗(yàn)組中微塑料的加入會(huì)使冰體的純度進(jìn)一步降低進(jìn)而加劇此現(xiàn)象的產(chǎn)生.在不同結(jié)冰溫度及方式條件下,微塑料的賦存對(duì)環(huán)境因子分布及遷移規(guī)律的影響與不同結(jié)冰比例相比并未發(fā)生較大變化,總體上仍然表現(xiàn)出實(shí)驗(yàn)組中冰體中環(huán)境因子的濃度與分配系數(shù)值較高,而冰下水體中環(huán)境因子濃度較低的規(guī)律,表明微塑料賦存對(duì)于物質(zhì)分布及遷移規(guī)律的影響并不取決于結(jié)冰條件的變化.

    2.3 不同結(jié)冰濃度下微塑料賦存對(duì)環(huán)境因子遷移的影響

    如圖4所示,在結(jié)冰比例與結(jié)冰溫度相同時(shí),初始濃度的增加伴隨著冰體中各環(huán)境因子濃度的增加,但分配系數(shù)值的變化趨勢(shì)卻與之相反,冰體中物質(zhì)濃度占總濃度的比例卻在逐漸降低.初始濃度提高后冰體對(duì)環(huán)境因子的排斥作用就越明顯,冰體中的環(huán)境因子濃度的增加量遠(yuǎn)低于物質(zhì)的遷移量,冰體對(duì)環(huán)境因子的俘獲作用低于冰體的排斥作用,所以環(huán)境因子向冰下水體遷移的能力增強(qiáng),遷移到冰下水中的環(huán)境因子變多,導(dǎo)致雖然冰體中環(huán)境因子的濃度在增加而值在降低,其原因一方面是因?yàn)槌跏紳舛鹊奶岣邔?dǎo)致了溶液凝固點(diǎn)的下降,使冰-水界面的穩(wěn)定性下降,出現(xiàn)了較多的冰晶體分支,使其純度降低;另一方面則是因?yàn)槌跏紳舛鹊奶岣呤沟萌芤褐袧撛诰Ш肆孔兇?晶核間相互撞擊的頻率和能量提高,再次成核的幾率再次增加[24],冰晶體生長(zhǎng)速率提高的同時(shí)值下降.

    但值得注意的是相同濃度下微塑料的賦存使實(shí)驗(yàn)組冰體中各環(huán)境因子濃度為對(duì)照組的1.17~ 1.49倍,冰下水體中環(huán)境因子為對(duì)照組的0.73~0.93倍,分配系數(shù)值下降范圍為0.07~0.18,這是由于實(shí)驗(yàn)組相較于對(duì)照組中加入微塑料相當(dāng)于提高了溶液的初始濃度,因此出現(xiàn)實(shí)驗(yàn)組中各環(huán)境因子在冰體中的濃度與分配系數(shù)值等與對(duì)照組不同的趨勢(shì).這與Gao等[25]利用冷凍工藝去除石油精煉廠和紙漿廠污水中有機(jī)污染物的結(jié)果相一致,隨著初始濃度的提高冰體中COD與TOC等的濃度就越高.

    綜上所述,不同結(jié)冰比例、結(jié)冰溫度及方式、初始濃度下微塑料賦存對(duì)于環(huán)境因子的分布及遷移規(guī)律的影響程度不同,微塑料的賦存對(duì)環(huán)境因子遷移的能力造成了影響,但并不能改變環(huán)境因子由冰體向冰下水體遷移的整體趨勢(shì).此外,微塑料作為一種吸附力較強(qiáng)的顆粒態(tài)物質(zhì),與其他懸浮顆粒相比除了具有改變環(huán)境中其他污染物濃度及遷移能力的性質(zhì)外,更值得注意的是,微塑料由于自身特性及毒性等對(duì)于環(huán)境中其他污染物的影響以及生物體造成的毒害作用更大,故在今后的研究中,建議在此基礎(chǔ)上更加注重微塑料后續(xù)所帶來(lái)的實(shí)際影響.

    2.4 結(jié)冰過(guò)程中微塑料對(duì)環(huán)境因子遷移機(jī)理的影響

    結(jié)晶指物質(zhì)從溶液、熔融物或蒸氣中以晶體狀態(tài)析出的過(guò)程,故結(jié)冰過(guò)程中環(huán)境因子在冰-水相中的遷移過(guò)程可以從結(jié)晶學(xué)角度解釋.冰晶體的形成主要包括成核與生長(zhǎng)2個(gè)過(guò)程,在水溫下降至冰點(diǎn)溫度以下溶液系統(tǒng)會(huì)自動(dòng)形成冰晶,當(dāng)水溫繼續(xù)降低時(shí),溶液內(nèi)部污染物的溶解度隨之下降,此時(shí)一部分物質(zhì)會(huì)被冰體析出,導(dǎo)致冰-水相間壓力的不平衡,提高了臨界尺寸冰核的可能性,冰層形成的速度加快.冰晶體的生長(zhǎng)速率與水分子在冰核上聚集的速度及冰-水界面的狀態(tài)有關(guān),冰-水界面的水分子僅通過(guò)界面躍遷就可以在冰核表面附著[26].如圖5所示,在冰-水界面周圍水分子由于氫鍵的作用形成冰核,并由于密度差的原因上浮附著在冰體底部,此時(shí)由于冰-水界面與液相間水分子及物質(zhì)濃度差的作用下,水分在液相內(nèi)向上遷移至冰-水界面處,其余物質(zhì)則向下遷移至液相中.即本文所研究的各環(huán)境因子由冰體遷移至水體中,表明在結(jié)冰過(guò)程中冰體有排斥現(xiàn)象出現(xiàn),冰晶體析出導(dǎo)致多數(shù)溶質(zhì)在冰下水體中濃縮,而微塑料的存在導(dǎo)致原本因密度差作用向液相擴(kuò)散的一部分環(huán)境因子分子聚集在微塑料上,進(jìn)而滯留在冰體內(nèi)部,并未向水體擴(kuò)散,即在結(jié)冰過(guò)程中微塑料的存在對(duì)環(huán)境因子向冰下水體遷移產(chǎn)生了一定阻礙作用,使冰體的凈化程度下降.

    圖5 結(jié)冰過(guò)程中水分子與環(huán)境因子遷移示意

    在不同的領(lǐng)域中共晶理論的含義不同,結(jié)冰過(guò)程中物質(zhì)的遷移同樣可以采取共晶理論進(jìn)行解釋,通常利用共晶理論與溶液的溫度-濃度之間的平衡曲線解釋.如圖6所示,曲線TED、曲線DA與直線TFDB這3條線將坐標(biāo)軸分為H、I、J、K 4個(gè)區(qū)域,H區(qū)間內(nèi)狀態(tài)表現(xiàn)為冰體與溶液共存;I區(qū)間狀態(tài)為冰、溶液與析出溶質(zhì)共存;J區(qū)間狀態(tài)為飽和溶液與析出溶質(zhì)共存;K區(qū)間內(nèi)狀態(tài)為溶質(zhì)及溶劑以溶液的狀態(tài)共存.本研究中各環(huán)境因子溶液未達(dá)到飽和狀態(tài),且所處溫度均低于溶液的冰點(diǎn),所以此時(shí)溶液的狀態(tài)在H區(qū)間,為固相的冰與環(huán)境因子溶液共存的狀態(tài).D點(diǎn)為共晶點(diǎn),此點(diǎn)為各態(tài)變化的臨界點(diǎn),溶劑與溶液之間以溶液狀態(tài)共存[11],若溫度繼續(xù)下降至TF以下,溶液中將有水與環(huán)境因子一同析出.因此環(huán)境因子不斷在新生長(zhǎng)的冰-水界面處濃縮,故一定量的微塑料賦存在溶液中相當(dāng)于提高溶液濃度,會(huì)一定程度上影響遷移能力的大小,卻不能改變冰體中各環(huán)境因子逐漸向冰下水體遷移.

    圖6 結(jié)冰過(guò)程中共晶理論示意

    3 結(jié)論

    3.1 結(jié)冰過(guò)程中冰體對(duì)環(huán)境因子具有排斥作用,無(wú)論是否存在微塑料,各環(huán)境因子均呈現(xiàn)出冰樣中濃度<原水樣濃度<冰下水樣濃度的規(guī)律,且各環(huán)境因子均由冰體遷移至冰下水體.

    3.2 不同結(jié)冰比例、結(jié)冰溫度及方式、初始濃度下微塑料賦存對(duì)環(huán)境因子的分布規(guī)律均有一定程度的影響,微塑料賦存導(dǎo)致冰體中環(huán)境因子濃度與分配系數(shù)K值的提高,致使原本應(yīng)遷移至冰下水體下的環(huán)境因子部分滯留于冰體中,進(jìn)而導(dǎo)致環(huán)境因子向冰下水體遷移的能力下降.但不同結(jié)冰條件下微塑料的賦存并未改變結(jié)冰過(guò)程中環(huán)境因子由冰體向冰下水體遷移的整體趨勢(shì),表明微塑料對(duì)于環(huán)境因子的攜帶作用小于冰體對(duì)環(huán)境因子的排斥作用.

    3.3 不同條件下結(jié)冰過(guò)程中微塑料對(duì)環(huán)境因子分布與遷移機(jī)理的影響均可以從結(jié)晶學(xué)角度與共晶理論得到解釋.

    [1] Verpoorter C, Kutser T, Seekell D A, et al. A global inventory of lakes based on high-resolution satellite imagery [J]. Geophysical Research Letters, 2014,41(18):6396-6402.

    [2] 楊婷婷.烏梁素海冰封期底泥氮形態(tài)變化機(jī)制研究 [D]. 北京:中央民族大學(xué), 2019.

    Yang T T. Mechanisms of nitrogen morphology changes in the bottom sediment of the wuliangsuhai Lake during the ice-out period [D]. Beijing: Minzu University of China, 2019.

    [3] 王志超,竇雅嬌,周 鑫,等.岱海冰封期微塑料與環(huán)境因子的關(guān)系及風(fēng)險(xiǎn)評(píng)價(jià)[J]. 中國(guó)環(huán)境科學(xué), 2022,42(2):889-896.

    Wang Z C, Dou Y J, Zhou X, et al. Relationship between microplastics occurrence and environmental factors and risk assessment during ice- covered period of the daihai Lake [J]. China Environmental Science, 2022,42(2):889-896.

    [4] 趙萬(wàn)里.結(jié)冰和融冰過(guò)程中馬拉硫磷的遷移規(guī)律研究 [D]. 煙臺(tái):煙臺(tái)大學(xué), 2021.

    Zhao W L. Study on the migration law of malathion in the processes of freezing and melting [D]. Yantai: Yantai University, 2021.

    [5] Cooper M G, Smith L C, Rennermalm A K, et al. Spectral attenuation coefficients from measurements of light transmission in bare ice on the greenland ice sheet [J]. The Cryosphere, 2021,15(4):1931-1953.

    [6] 許冬雪,李 興,王 勇,等.冰封期烏梁素海不同形態(tài)氮,磷和葉綠素a的空間分布特征及其響應(yīng)關(guān)系[J]. 生態(tài)環(huán)境學(xué)報(bào), 2021,30(9): 1855-1864.

    Xu D X, Li X, Wang Y, et al. Spatial distribution characteristics and the response of different forms of nitrogen, phosphorus and chlorophyll-a in lake ulansuhai during the frozen period [J]. Ecology and Environment Sciences.

    [7] Kataoka T, Nihei Y, Kudou K, et al. Assessment of the sources and inflow processes of microplastics in the river environments of japan [J]. Environmental pollution, 2019,244:958-965.

    [8] 王志超,楊建林,楊 帆,等.烏梁素海冰蓋中微塑料的分布特征及其與鹽度、葉綠素a的響應(yīng)關(guān)系[J]. 環(huán)境科學(xué), 2021,42(2):673-680.

    Wang Z C, Yang J L, Yang F, et al. Distribution characteristics of microplastics in ice sheets and its response to salinity and chlorophyll a in the lake wuliangsuhai [J]. Environmental Science, 2021,42(2): 673-680.

    [9] Hoffmann L, Eggers S L, Allhusen E, et al. Interactions between the ice algae fragillariopsis cylindrus and microplastics in sea ice [J]. Environment International, 2020,139:105697.

    [10] 呂宏洲,李暢游,史小紅,等.不同條件下烏梁素海污染物在冰-水體系中分布規(guī)律的模擬[J]. 湖泊科學(xué), 2015,27(6):1151-1158.

    Lv H Z, Li C Y, Shi X H, et al. Pollutant distribution under different conditions in lake ulansuhai ice-water system [J]. Journal of Lake Sciences, 2015,27(6):1151-1158.

    [11] 張 巖,任方云,唐元慶,等.融冰過(guò)程中鐵離子和錳離子的遷移規(guī)律[J]. 中國(guó)環(huán)境科學(xué), 2021,41(5):2391-2398.

    Zhang Y, Ren F Y, Tang Y Q, et al. Migration of iron and manganese ions during ice melting [J]. China Environmental Science, 2021,41(5): 2391-2398.

    [12] 唐元慶.水體結(jié)冰和融冰過(guò)程中鐵、錳、鈣、鎂的遷移規(guī)律研究 [D]. 煙臺(tái):煙臺(tái)大學(xué), 2020.

    Tang Y Q. The migration law of iron, manganese, calcium and magnesium in water icing and melting processes [D]. Yan Tai: Yantai University, 2020.

    [13] GB/T 12763.4-2007 海洋調(diào)查規(guī)范第4部分:海水化學(xué)要素調(diào)查[S].

    GB/T 12763.4-2007 Specifications for oceanographic survey part: survey of chemical parameters in sea water [S].

    [14] 喬延龍,殷小亞,肖廣俠,等.懸浮物脅迫中國(guó)對(duì)蝦幼體的急性毒性研究[J]. 漁業(yè)科學(xué)進(jìn)展, 2019,40(3):50-56.

    Qiao Y L, Yin X Y, Xiao G X, et al. Acute toxicity effects of suspended solids stress on fenneropenaeus chinensis larvae [J]. Progress in Fishery Sciences, 2019,40(3):50-56.

    [15] 國(guó)家環(huán)境保護(hù)總局《水和廢水監(jiān)測(cè)分析方法》編委會(huì).水和廢水監(jiān)測(cè)分析方法(第四版) [M]. 北京:中國(guó)環(huán)境科學(xué)出版社, 2002.

    State Environmental Protection Administration. “Analysis methods for water and wastewater” Editorial Board. Analysis Methods for Water and Wastewater(4th edition) [M]. Beijing: China Environmental Science Press, 2002.

    [16] HJ636-2012 水質(zhì)總氮的測(cè)定堿性過(guò)硫酸鉀消解紫外分光光度法[S].

    HJ636-2012 Water quality-determination of total nitrogen-alkaline potassium persulfate digestion uv spectrophotometric method [S].

    [17] GB 11893-89 水質(zhì)總磷的測(cè)定鉬酸銨分光光度法[S].

    GB 11893-89 Water quality-determination of total phosphorus- ammonium molybdate spectrophotometric method [S].

    [18] HJ 828-2017 水質(zhì)化學(xué)需氧量的測(cè)定重鉻酸鹽法[S].

    HJ 828-2017 Water quality-determination of the chemical oxygen demand-dichromate method [S].

    [19] GB 11901-89 水質(zhì)懸浮物的測(cè)定重量法[S].

    GB 11901-89 Water quality-determination of suspended substance- gravimetric method [S].

    [20] 于愛(ài)鑫.模擬結(jié)冰和融冰過(guò)程中阿特拉津的遷移規(guī)律 [D]. 煙臺(tái):煙臺(tái)大學(xué), 2020.

    Yu A X. Migration law of atrazine during simulated freezing and melting processes [D]. Yantai: Yantai University, 2020.

    [21] 袁海英,侯 磊,梁?jiǎn)⒈?等.滇池近岸水體微塑料污染與富營(yíng)養(yǎng)化的相關(guān)性 [J]. 環(huán)境科學(xué), 2021,42(7):3166-3175.

    Yuan H Y, Hou L, Liang Q B, et al. Correlation between microplastics pollution and eutrophication in the near shore waters of dianchi lake [J]. Environmental Science, 2021,42(7):3166-3175.

    [22] Yu T, Ma J Z, Li Q. Factors affecting ice crystal purity during freeze concentration process for urine treatment [J]. Journal of Harbin Institute of Technology, 2007,14(5):593-597.

    [23] Salonen K, Lepp?ranta M, Viljanen M, et al. Perspectives in winter limnology: closing the annual cycle of freezing lakes [J]. Aquatic Ecology, 2009,43(3):609-616.

    [24] Luo C, Chen W, Han W. Experimental study on factors affecting the quality of ice crystal during the freezing concentration for the brackish water [J]. Desalination, 2010,260(1-3):231-238.

    [25] Gao W, Habib M, Smith D W. Removal of organic contaminants and toxiciy from industrial effluents using freezing processes [J]. Desalination, 2009,245(1-3):108-119.

    [26] 王赫偉.冰封期河流污染物變化規(guī)律及機(jī)理研究-以太子河本溪城區(qū)段為例 [D]. 沈陽(yáng):遼寧大學(xué), 2021.

    Wang H W. Study on the variation law and mechanism of pollutants in river during freezing period-taking benxi urban section of taizi river as an example [D]. Shenyang: Liaoning University, 2021.

    Effect of microplastics on the transport of environmental factors during icing and the mechanism of action.

    WANG Zhi-chao, DOU Ya-jiao, KANG Yan-qiu, ZHOU Xin, YANG Wen-huan, JING Shuang-yi, LI Wei-ping*

    (School of Energy and Environment, Inner Mongolia University of Science and Technology, Baotou 014010, China)., 2022,42(11):5369~5377

    In order to explore the influence of microplastics reserve on the distribution and migration of typical environmental factors (total nitrogen, total phosphorus, salinity, chemical oxygen demand and suspended solids) during the icing process, the indoor simulations were adopted to study the distribution and migration of typical environmental factors in the ice-water phase under different conditions (icing ratio, freezing temperature and mode and initial concentration). The migration ability of environmental was characterized by material distribution coefficient(). The results showed that the ice had a repulsive effect on environmental factors during the freezing process, and the occurrence of microplastics under different conditions had a certain degree of influence on the distribution of environmental factors. The environmental factors that should have migrated to the water under the ice were partially retained in the ice due to the characteristics of microplastics, resulting in an increase in the concentration of environmental factors in the ice, which was 1.13~1.49 times that of the control group. At the same time, the concentration of environmental factors in the subglacial water decreased, which was 0.73~0.93 times that of the control group. The presence of microplastics also increased the partition coefficientof environmental factors by 0.04~0.18 and decreased the ability of environmental factors to migrate to subglacial waters. However, the occurrence of microplastics did not change the trend of environmental factors transferring from ice to subglacial water during the freezing process. The carrying effect of microplastics on environmental factors was less than the repulsion effect of ice on environmental factors. At the same time, the influence of microplastics on the distribution and migration mechanism of environmental factors during freezing under different conditions could be explained from the perspective of crystallography and cocrystal theory.

    icing process;microplastics;environmental factors;transport pattern;mechanism of action

    X524

    A

    1000-6923(2022)11-5369-09

    王志超(1988-),男,山東德州人,副教授,博士,從事環(huán)境污染控制與生態(tài)修復(fù).發(fā)表論文30余篇.

    2022-04-25

    國(guó)家自然科學(xué)基金資助項(xiàng)目(42007119);內(nèi)蒙古自治區(qū)高等學(xué)校青年科技英才支持計(jì)劃項(xiàng)目(NJYT22066)

    * 責(zé)任作者, 教授, sjlwp@163.com

    猜你喜歡
    冰體結(jié)冰塑料
    彈體高速侵徹冰體研究
    通體結(jié)冰的球
    高速?gòu)楏w侵徹冰材料過(guò)程數(shù)值模擬研究
    冰體質(zhì)量和撞擊角度對(duì)船首結(jié)構(gòu)碰撞性能的影響
    基于船-水-冰耦合技術(shù)的撞擊參數(shù)對(duì)船冰碰撞性能的影響
    冬天,玻璃窗上為什么會(huì)結(jié)冰花?
    魚缸結(jié)冰
    塑料也高級(jí)
    Coco薇(2016年8期)2016-10-09 16:58:11
    塑料
    塑料和聚合物的超聲焊接
    日韩伦理黄色片| 国内揄拍国产精品人妻在线| 久久久久网色| 18禁裸乳无遮挡动漫免费视频| 爱豆传媒免费全集在线观看| 久久久久久久大尺度免费视频| 成人毛片60女人毛片免费| 亚洲精品456在线播放app| 婷婷色综合www| 久久久久网色| 伦精品一区二区三区| 18禁在线无遮挡免费观看视频| 亚洲欧美精品自产自拍| 午夜老司机福利剧场| 天堂俺去俺来也www色官网| 亚洲高清免费不卡视频| 亚洲欧美日韩无卡精品| 内射极品少妇av片p| 男女国产视频网站| 91久久精品国产一区二区三区| 最黄视频免费看| 久久国产精品大桥未久av | 日韩一本色道免费dvd| 国产伦在线观看视频一区| 国产日韩欧美在线精品| 最近的中文字幕免费完整| 亚洲欧美清纯卡通| 中文欧美无线码| 全区人妻精品视频| kizo精华| 欧美日韩国产mv在线观看视频 | 又爽又黄a免费视频| 91久久精品电影网| 日韩电影二区| 成人影院久久| 久久女婷五月综合色啪小说| 国产 精品1| 国产av国产精品国产| 亚洲精品色激情综合| 精品国产一区二区三区久久久樱花 | 夜夜爽夜夜爽视频| 大片免费播放器 马上看| 日韩制服骚丝袜av| 最黄视频免费看| 极品少妇高潮喷水抽搐| 麻豆乱淫一区二区| 日韩 亚洲 欧美在线| 日韩av免费高清视频| 久久热精品热| 免费人成在线观看视频色| 少妇精品久久久久久久| 亚洲精品成人av观看孕妇| 国产精品熟女久久久久浪| 六月丁香七月| 2021少妇久久久久久久久久久| 亚洲四区av| 韩国av在线不卡| a级毛色黄片| 亚洲精品456在线播放app| 观看美女的网站| 久久精品人妻少妇| 久久精品国产亚洲av涩爱| 国模一区二区三区四区视频| 亚洲精品一区蜜桃| 国产男人的电影天堂91| 99热国产这里只有精品6| av.在线天堂| 大陆偷拍与自拍| 在线免费观看不下载黄p国产| 老女人水多毛片| 国产成人a∨麻豆精品| 丰满乱子伦码专区| 午夜福利高清视频| 精品亚洲成国产av| 国产午夜精品久久久久久一区二区三区| 少妇高潮的动态图| av视频免费观看在线观看| 亚洲四区av| 啦啦啦中文免费视频观看日本| 在线观看免费日韩欧美大片 | 久久久久久伊人网av| 制服丝袜香蕉在线| 国产淫片久久久久久久久| 美女福利国产在线 | 日本av手机在线免费观看| 丰满迷人的少妇在线观看| 国产精品欧美亚洲77777| 日韩精品有码人妻一区| 综合色丁香网| 午夜免费鲁丝| 成人无遮挡网站| 日本黄色片子视频| 成人无遮挡网站| 好男人视频免费观看在线| 免费av中文字幕在线| 精品久久国产蜜桃| 亚洲欧洲国产日韩| 国产精品麻豆人妻色哟哟久久| 只有这里有精品99| 免费大片黄手机在线观看| 久久久久久久亚洲中文字幕| 中文字幕免费在线视频6| 日韩一区二区三区影片| 一级毛片我不卡| 综合色丁香网| www.色视频.com| 免费播放大片免费观看视频在线观看| av播播在线观看一区| 日本vs欧美在线观看视频 | 亚洲精品一二三| h视频一区二区三区| 久久婷婷青草| 国产成人a区在线观看| 亚洲欧美一区二区三区国产| 天美传媒精品一区二区| 亚洲美女视频黄频| 亚洲熟女精品中文字幕| 美女主播在线视频| 国产一区二区三区av在线| 成人国产麻豆网| 亚洲精品,欧美精品| 欧美变态另类bdsm刘玥| 青青草视频在线视频观看| 三级国产精品欧美在线观看| 国产在线视频一区二区| 80岁老熟妇乱子伦牲交| 国产精品蜜桃在线观看| 国产在线一区二区三区精| 99热这里只有精品一区| 联通29元200g的流量卡| 欧美+日韩+精品| 久久久久久久久久久免费av| av在线蜜桃| 男女免费视频国产| 只有这里有精品99| 高清日韩中文字幕在线| 夫妻午夜视频| 精品一区在线观看国产| 亚洲综合精品二区| 亚洲国产色片| 妹子高潮喷水视频| 亚洲精品国产av蜜桃| 免费黄网站久久成人精品| 国产久久久一区二区三区| 免费看av在线观看网站| 日韩不卡一区二区三区视频在线| 久久久久久久久久久免费av| 人妻夜夜爽99麻豆av| 伊人久久国产一区二区| 五月伊人婷婷丁香| 亚洲一区二区三区欧美精品| .国产精品久久| 亚洲精品乱码久久久久久按摩| 九九久久精品国产亚洲av麻豆| 日本wwww免费看| 91aial.com中文字幕在线观看| 一个人免费看片子| 国产精品欧美亚洲77777| 亚洲精品色激情综合| 久久99蜜桃精品久久| 国产女主播在线喷水免费视频网站| 午夜福利网站1000一区二区三区| 国产av码专区亚洲av| 久久毛片免费看一区二区三区| 成人特级av手机在线观看| 久久久久久人妻| 国产视频内射| 日韩 亚洲 欧美在线| 老司机影院成人| 国产成人91sexporn| 免费播放大片免费观看视频在线观看| 成年美女黄网站色视频大全免费 | 七月丁香在线播放| 久久久久久久国产电影| 免费不卡的大黄色大毛片视频在线观看| 天天躁夜夜躁狠狠久久av| 亚洲av福利一区| 18禁在线无遮挡免费观看视频| 在线亚洲精品国产二区图片欧美 | 男女国产视频网站| 欧美3d第一页| 内射极品少妇av片p| 久久精品久久久久久噜噜老黄| 一区二区三区免费毛片| 亚洲国产成人一精品久久久| 一个人看视频在线观看www免费| 男人舔奶头视频| 日本与韩国留学比较| 激情五月婷婷亚洲| videossex国产| 91精品国产九色| .国产精品久久| 亚洲精品乱久久久久久| 最近手机中文字幕大全| 亚洲人成网站高清观看| 永久免费av网站大全| 国产精品三级大全| 国产中年淑女户外野战色| 联通29元200g的流量卡| 97超碰精品成人国产| 欧美一区二区亚洲| 黄色怎么调成土黄色| 日韩 亚洲 欧美在线| av在线老鸭窝| 欧美最新免费一区二区三区| 亚洲真实伦在线观看| 国产精品免费大片| 熟女人妻精品中文字幕| 人人妻人人看人人澡| 免费人妻精品一区二区三区视频| 国产视频首页在线观看| 亚洲精品日本国产第一区| 99久久人妻综合| 在线观看人妻少妇| 春色校园在线视频观看| 一级黄片播放器| 高清日韩中文字幕在线| 国产成人免费无遮挡视频| 国产视频内射| 男人添女人高潮全过程视频| 日韩av在线免费看完整版不卡| 午夜免费鲁丝| 欧美3d第一页| 亚洲国产精品999| 国产精品国产三级专区第一集| 深夜a级毛片| 尾随美女入室| 亚洲中文av在线| 成人黄色视频免费在线看| 欧美日韩视频高清一区二区三区二| 亚洲精品国产av成人精品| 精品99又大又爽又粗少妇毛片| 午夜福利在线观看免费完整高清在| 国产国拍精品亚洲av在线观看| a 毛片基地| 欧美人与善性xxx| 九九久久精品国产亚洲av麻豆| 大香蕉97超碰在线| 成年女人在线观看亚洲视频| 在线看a的网站| 成人18禁高潮啪啪吃奶动态图 | 欧美最新免费一区二区三区| 最黄视频免费看| 亚洲色图综合在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美日韩另类电影网站 | 色5月婷婷丁香| 九九在线视频观看精品| xxx大片免费视频| 亚洲精品国产av成人精品| 伊人久久国产一区二区| 大陆偷拍与自拍| 亚洲精品乱久久久久久| 国产老妇伦熟女老妇高清| 日韩一区二区视频免费看| 日韩av在线免费看完整版不卡| 搡老乐熟女国产| 青春草国产在线视频| 久久鲁丝午夜福利片| 少妇丰满av| 精品久久久精品久久久| 亚洲人成网站高清观看| 亚洲精品视频女| 免费大片18禁| 亚洲图色成人| 99视频精品全部免费 在线| 成人特级av手机在线观看| av黄色大香蕉| 街头女战士在线观看网站| av不卡在线播放| 亚洲色图av天堂| 欧美性感艳星| 亚洲欧美成人精品一区二区| 啦啦啦在线观看免费高清www| 内地一区二区视频在线| 观看av在线不卡| videossex国产| 午夜免费观看性视频| 你懂的网址亚洲精品在线观看| 日韩精品有码人妻一区| 精品久久久久久久久av| 久久久久久久久久久免费av| freevideosex欧美| 欧美少妇被猛烈插入视频| 韩国av在线不卡| 中文资源天堂在线| 日本av免费视频播放| 欧美成人a在线观看| 精品人妻一区二区三区麻豆| 国产成人精品久久久久久| 国产精品国产三级国产专区5o| 80岁老熟妇乱子伦牲交| 91午夜精品亚洲一区二区三区| 国产精品爽爽va在线观看网站| 干丝袜人妻中文字幕| 国产精品一区二区三区四区免费观看| 一区二区av电影网| 亚洲综合精品二区| 国产精品人妻久久久久久| 少妇人妻一区二区三区视频| 18禁在线无遮挡免费观看视频| 美女中出高潮动态图| 这个男人来自地球电影免费观看 | 久久久久精品性色| 久久久久久人妻| 嫩草影院入口| 国产精品国产三级国产专区5o| 亚洲最大成人中文| 国产精品av视频在线免费观看| 97热精品久久久久久| 高清在线视频一区二区三区| freevideosex欧美| 男人和女人高潮做爰伦理| 亚洲欧美日韩东京热| 免费av中文字幕在线| 晚上一个人看的免费电影| 亚洲成人av在线免费| 国产精品一区二区在线不卡| 涩涩av久久男人的天堂| 久久久亚洲精品成人影院| 久久6这里有精品| av在线app专区| 我的老师免费观看完整版| 美女xxoo啪啪120秒动态图| 国产日韩欧美亚洲二区| 亚洲熟女精品中文字幕| 在线免费十八禁| 国产在线视频一区二区| 日韩伦理黄色片| 午夜老司机福利剧场| 中文在线观看免费www的网站| 黄色配什么色好看| 男人舔奶头视频| 国产精品一区二区在线不卡| 美女cb高潮喷水在线观看| 亚洲精品乱码久久久v下载方式| 天堂中文最新版在线下载| 成人特级av手机在线观看| 干丝袜人妻中文字幕| 日本免费在线观看一区| 免费久久久久久久精品成人欧美视频 | av在线观看视频网站免费| 97热精品久久久久久| 国产成人午夜福利电影在线观看| av黄色大香蕉| 色综合色国产| 亚洲精品色激情综合| 日本黄色片子视频| 免费观看av网站的网址| 多毛熟女@视频| 免费黄频网站在线观看国产| 日本av免费视频播放| 最近2019中文字幕mv第一页| 国产熟女欧美一区二区| 久久99热这里只频精品6学生| 黄色欧美视频在线观看| 国产成人91sexporn| 亚洲欧美日韩另类电影网站 | 亚洲精品日韩在线中文字幕| 日本爱情动作片www.在线观看| 成人美女网站在线观看视频| 哪个播放器可以免费观看大片| 日韩 亚洲 欧美在线| 精品国产露脸久久av麻豆| 天堂8中文在线网| 女的被弄到高潮叫床怎么办| 寂寞人妻少妇视频99o| 搡老乐熟女国产| 亚洲国产欧美人成| 国产 精品1| 免费av不卡在线播放| 国产精品三级大全| 成年av动漫网址| 天天躁夜夜躁狠狠久久av| 亚洲婷婷狠狠爱综合网| 午夜视频国产福利| 久久久成人免费电影| 老熟女久久久| 欧美三级亚洲精品| 日日摸夜夜添夜夜爱| 久久鲁丝午夜福利片| 欧美zozozo另类| 女人久久www免费人成看片| av在线播放精品| 亚洲最大成人中文| 亚洲婷婷狠狠爱综合网| 免费观看av网站的网址| a 毛片基地| 久久精品人妻少妇| 啦啦啦啦在线视频资源| 亚洲精品日韩av片在线观看| 啦啦啦啦在线视频资源| 在线看a的网站| 好男人视频免费观看在线| 欧美zozozo另类| videos熟女内射| 久久精品久久久久久噜噜老黄| 成年免费大片在线观看| 日本黄色片子视频| 亚洲av中文av极速乱| 国产 精品1| 成人毛片60女人毛片免费| 国产免费一区二区三区四区乱码| 中国美白少妇内射xxxbb| 啦啦啦在线观看免费高清www| 新久久久久国产一级毛片| 在线精品无人区一区二区三 | 在线亚洲精品国产二区图片欧美 | 欧美高清成人免费视频www| 日本猛色少妇xxxxx猛交久久| 久久久久久久亚洲中文字幕| www.色视频.com| 极品教师在线视频| 免费高清在线观看视频在线观看| 成年人午夜在线观看视频| 性高湖久久久久久久久免费观看| 在线观看国产h片| 欧美97在线视频| 欧美精品国产亚洲| 亚洲av成人精品一区久久| 亚洲怡红院男人天堂| 嘟嘟电影网在线观看| 91aial.com中文字幕在线观看| 国产成人aa在线观看| 两个人的视频大全免费| 春色校园在线视频观看| 精品国产一区二区三区久久久樱花 | 国产伦精品一区二区三区视频9| 我要看黄色一级片免费的| 99国产精品免费福利视频| 一级毛片电影观看| 黄色欧美视频在线观看| 一级毛片 在线播放| 成人二区视频| 亚洲一区二区三区欧美精品| 国产成人aa在线观看| 熟女电影av网| 天堂8中文在线网| 国产av一区二区精品久久 | 男人和女人高潮做爰伦理| 性高湖久久久久久久久免费观看| 一级二级三级毛片免费看| 日本欧美视频一区| 欧美bdsm另类| 啦啦啦啦在线视频资源| 亚洲一区二区三区欧美精品| 在现免费观看毛片| 夫妻性生交免费视频一级片| 欧美另类一区| 观看免费一级毛片| 最黄视频免费看| 免费看av在线观看网站| 男女边摸边吃奶| 日韩亚洲欧美综合| 三级经典国产精品| a级毛片免费高清观看在线播放| 久久久久久久亚洲中文字幕| 亚洲国产av新网站| 国产v大片淫在线免费观看| 日韩欧美 国产精品| 久久午夜福利片| 亚洲人成网站在线观看播放| 久久久久精品久久久久真实原创| 人妻少妇偷人精品九色| 国产精品久久久久久久电影| 国精品久久久久久国模美| 国产伦理片在线播放av一区| 亚洲自偷自拍三级| 成年av动漫网址| 国产极品天堂在线| 下体分泌物呈黄色| 成年美女黄网站色视频大全免费 | 中文乱码字字幕精品一区二区三区| 亚洲中文av在线| 嫩草影院新地址| 国产伦精品一区二区三区四那| 草草在线视频免费看| 国产精品伦人一区二区| 91久久精品电影网| 国产在视频线精品| 伊人久久精品亚洲午夜| 欧美少妇被猛烈插入视频| 中文在线观看免费www的网站| 日本色播在线视频| 亚洲人与动物交配视频| 色网站视频免费| 色视频在线一区二区三区| 欧美变态另类bdsm刘玥| 亚洲精品色激情综合| 精品亚洲成国产av| 久久久久久久久大av| 亚洲av不卡在线观看| 亚洲激情五月婷婷啪啪| 色婷婷久久久亚洲欧美| 九草在线视频观看| 久久ye,这里只有精品| 欧美变态另类bdsm刘玥| 夫妻午夜视频| 国产伦精品一区二区三区四那| 成人综合一区亚洲| www.av在线官网国产| 亚洲色图综合在线观看| 精品酒店卫生间| 久久午夜福利片| 美女脱内裤让男人舔精品视频| 超碰av人人做人人爽久久| 国产片特级美女逼逼视频| 精品熟女少妇av免费看| 色婷婷久久久亚洲欧美| 两个人的视频大全免费| 狂野欧美白嫩少妇大欣赏| 日日啪夜夜撸| 欧美丝袜亚洲另类| 婷婷色麻豆天堂久久| 国产69精品久久久久777片| av免费观看日本| av线在线观看网站| 七月丁香在线播放| 欧美97在线视频| 男女国产视频网站| kizo精华| 91在线精品国自产拍蜜月| 看十八女毛片水多多多| 中文字幕精品免费在线观看视频 | 久久精品人妻少妇| 日韩一区二区三区影片| 男女边摸边吃奶| 国产乱来视频区| 欧美日韩视频精品一区| 欧美成人午夜免费资源| 亚洲一级一片aⅴ在线观看| 免费少妇av软件| 免费av中文字幕在线| 色5月婷婷丁香| 久久6这里有精品| 亚洲国产高清在线一区二区三| 精品人妻熟女av久视频| 汤姆久久久久久久影院中文字幕| 国产片特级美女逼逼视频| 国产有黄有色有爽视频| 老司机影院成人| 交换朋友夫妻互换小说| 久久精品夜色国产| 欧美日韩精品成人综合77777| 亚洲精品一区蜜桃| 亚洲精品色激情综合| 久久久精品免费免费高清| 99九九线精品视频在线观看视频| 看免费成人av毛片| 美女内射精品一级片tv| 国产精品久久久久久久电影| 我要看黄色一级片免费的| 亚洲电影在线观看av| 简卡轻食公司| 亚洲精品,欧美精品| 国产精品久久久久久久久免| 视频区图区小说| 国产深夜福利视频在线观看| 七月丁香在线播放| 在线观看一区二区三区激情| 久久久精品94久久精品| 亚洲国产日韩一区二区| 插阴视频在线观看视频| 免费观看在线日韩| 色婷婷久久久亚洲欧美| 成人毛片a级毛片在线播放| 欧美精品亚洲一区二区| 中文字幕av成人在线电影| 美女高潮的动态| 亚洲欧洲国产日韩| 精品人妻一区二区三区麻豆| 我的女老师完整版在线观看| 一级爰片在线观看| 777米奇影视久久| 国产精品国产av在线观看| 久久人人爽av亚洲精品天堂 | 日本爱情动作片www.在线观看| 亚洲色图av天堂| 高清欧美精品videossex| 中文字幕av成人在线电影| a级毛色黄片| 欧美日韩精品成人综合77777| av在线app专区| 啦啦啦在线观看免费高清www| 精品视频人人做人人爽| 亚洲国产精品一区三区| www.av在线官网国产| 夜夜骑夜夜射夜夜干| 国产永久视频网站| 国产精品国产三级国产专区5o| 婷婷色综合大香蕉| 亚洲精品一区蜜桃| 黑人高潮一二区| 久久久欧美国产精品| 日本av手机在线免费观看| .国产精品久久| 成年人午夜在线观看视频| 免费高清在线观看视频在线观看| 在线亚洲精品国产二区图片欧美 | 国产又色又爽无遮挡免| 最近中文字幕高清免费大全6| 一个人看的www免费观看视频| 五月天丁香电影| 国产极品天堂在线| 欧美成人午夜免费资源| 大片电影免费在线观看免费| 91精品一卡2卡3卡4卡| 日日摸夜夜添夜夜爱| 免费观看在线日韩| 久久久久性生活片| 亚洲欧美日韩无卡精品| 老师上课跳d突然被开到最大视频| 国产国拍精品亚洲av在线观看| 日韩在线高清观看一区二区三区| 男女边吃奶边做爰视频|