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

    基于SOM 模型的珠江河口河網(wǎng)鉈生態(tài)風(fēng)險(xiǎn)評(píng)估

    2021-06-03 09:25:24藍(lán)璇利鋒張超董漢英楊清書(shū)余明輝文汝兵楊玉潔
    熱帶海洋學(xué)報(bào) 2021年3期
    關(guān)鍵詞:小潮馬口石龍

    藍(lán)璇, 利鋒, 張超, 董漢英, 楊清書(shū), 余明輝, 文汝兵, 楊玉潔

    1. 華南理工大學(xué)土木與交通學(xué)院, 廣東 廣州 510641;

    2. 中山大學(xué)大氣科學(xué)學(xué)院, 廣東 珠海 519082;

    3. 中山大學(xué)海洋工程與技術(shù)學(xué)院, 廣東 珠海 519082;

    4. 武漢大學(xué)水資源與水電工程國(guó)家重點(diǎn)實(shí)驗(yàn)室, 湖北 武漢 430072

    珠江, 為中國(guó)第二大河流(按流量), 流經(jīng)云南、貴州、廣西、廣東、湖南、江西6 個(gè)省(區(qū)), 流域總?cè)丝跒?980 萬(wàn), 其主要支流為西江、東江、北江。珠江在下游三角洲漫流成網(wǎng)河區(qū), 經(jīng)由八大口門(mén)流入南海。珠江河口河網(wǎng)地區(qū), 是世界上最為復(fù)雜的河口區(qū)之一, 也是粵港澳大灣區(qū)的核心區(qū)域。隨著近20 年來(lái)的經(jīng)濟(jì)發(fā)展, 珠江流域礦產(chǎn)被大量開(kāi)采冶煉, 導(dǎo)致很多毒害重金屬元素排入珠江, 對(duì)流域的水環(huán)境安全造成嚴(yán)重的影響和威脅(劉娟 等,2008)。加強(qiáng)生態(tài)文明建設(shè)和生態(tài)環(huán)境保護(hù), 既是大灣區(qū)建設(shè)世界一流灣區(qū)的重要支撐(萬(wàn)軍 等, 2019),也是重要目標(biāo)要求, 維持良好的水環(huán)境質(zhì)量是建設(shè)綠色生態(tài)灣區(qū)的重要任務(wù)之一(唐天均, 2018)。以綠色發(fā)展推動(dòng)粵港澳大灣區(qū)建設(shè)成為廣東省當(dāng)前面臨的重大課題, 優(yōu)良的生態(tài)環(huán)境是支撐粵港澳大灣區(qū)經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展的先決條件 (孫加龍, 2018)。《粵港澳大灣區(qū)發(fā)展規(guī)劃綱要》中要求強(qiáng)化珠三角水資源安全保障, 因此水安全保障任務(wù)更為迫切(畢建培 等, 2019)。

    鉈(thallium)是最具毒性的稀有分散元素之一,在食物鏈中具有富集性(劉敬勇 等, 2010)。當(dāng)水中鉈的濃度為1~60mg·L-1時(shí), 可致使水體中的魚(yú)死亡。鉈對(duì)成人的最小致死劑量?jī)H為12mg·kg-1體重,對(duì)兒童則為8.8~15mg·kg-1體重(陳學(xué)鵬,2009)。然而目前對(duì)珠江鉈的研究主要集中于西江和北江中下游,針對(duì)珠江河口河網(wǎng)鉈的研究較少。有研究者分析了北江嚴(yán)重鉈污染的產(chǎn)生原因, 發(fā)現(xiàn)由于鉈的化學(xué)活動(dòng)性, 石灰中和過(guò)程的沉降作用甚微, 金屬冶煉廠排放口廢水的鉈含量高達(dá)600~700μg·L-1以上(陳永亨 等 2013)。對(duì)廣東西江流域的某超大型含鉈硫鐵礦和硫酸廠附近的河流沉積物中鉈污染的狀況的研究結(jié)果表明, 研究區(qū)域內(nèi)整條河流都呈現(xiàn)出明顯的鉈污染現(xiàn)象, 含量高達(dá)1.30~17.31mg·kg-1, 即使是進(jìn)入西江入口附近, 在大量西江水體的稀釋作用下,沉積物中鉈的含量仍然高達(dá)我國(guó)淺海沉積物背景值的三倍; 除礦區(qū)附近2~4km 內(nèi)的河流表層水體的鉈濃度有異常值外, 鉈在這些河水中的濃度并不高,為0.01~0.10μg·L-1, 均低于或接近于我國(guó)飲用水的最大允許濃度值(解小凡 等, 2015; 劉娟 等, 2015)。

    目前, 自組織映射神經(jīng)網(wǎng)絡(luò)(self-organizing maps, SOM)方法主要應(yīng)用于數(shù)據(jù)挖掘、統(tǒng)計(jì)學(xué)等方面, 尚無(wú)將此方法應(yīng)用于污染物尤其重金屬分析的公開(kāi)報(bào)道。傳統(tǒng)評(píng)價(jià)方法大多都是運(yùn)用數(shù)理統(tǒng)計(jì)和經(jīng)驗(yàn)等方法將污染物所測(cè)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)整理, 根據(jù)一定的規(guī)律(公式)對(duì)其進(jìn)行分類, 最后得到污染物的分類情況和污染水平。傳統(tǒng)評(píng)價(jià)方法一般受人為因素影響較大, 且容易造成誤差累積。而SOM 方法是自主學(xué)習(xí)方法, 由計(jì)算機(jī)根據(jù)程序自主對(duì)數(shù)據(jù)進(jìn)行分析, 將基因相似的數(shù)據(jù)歸為一類, 從而得到分類結(jié)果, 能夠精細(xì)準(zhǔn)確地對(duì)像鉈這種靈敏度高、閾值低的重金屬進(jìn)行聚類分析, 極大限度地排除了人為經(jīng)驗(yàn)和各種統(tǒng)計(jì)誤差造成的不準(zhǔn)確影響。

    本文在珠江河口與河網(wǎng)處總共布設(shè)11 個(gè)斷面進(jìn)行同步采樣, 研究分析鉈的暴露水平以及鉈含量的時(shí)空分布情況, 并基于潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)法, 運(yùn)用SOM 基因表達(dá)聚類分析模型得到鉈的濃度分布結(jié)果和珠江河口河網(wǎng)潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果, 為珠江河口河網(wǎng)重金屬鉈的污染防治提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 SOM 基因表達(dá)聚類分析

    SOM 是基于神經(jīng)網(wǎng)絡(luò)方式的數(shù)據(jù)矩陣和可視化方式(Kohonen, 2001), 由Kohonen 在1982 年提出,是一種無(wú)監(jiān)督學(xué)習(xí)的神經(jīng)元網(wǎng)絡(luò)模型。其分析過(guò)程是先找到一組中心點(diǎn)(又稱為codebook vector), 然后根據(jù)最相似原則把數(shù)據(jù)集的每個(gè)對(duì)象映射到對(duì)應(yīng)的中心點(diǎn)。在神經(jīng)網(wǎng)絡(luò)術(shù)語(yǔ)中, 每個(gè)神經(jīng)元對(duì)應(yīng)于一個(gè)中心點(diǎn)。

    SOM 數(shù)據(jù)集中的每個(gè)對(duì)象每次處理一個(gè), 判斷最近的中心點(diǎn), 然后更新中心點(diǎn)。中心點(diǎn)之間存在拓?fù)湫螤铐樞? 在更新一個(gè)中心點(diǎn)的同時(shí), 鄰近的中心點(diǎn)也會(huì)隨著更新, 直到達(dá)到設(shè)定的閾值或中心點(diǎn)不再有顯著變化。最終獲得一系列的中心點(diǎn)(codes)隱式地定義多個(gè)簇, 與這個(gè)中心點(diǎn)最近的對(duì)象歸為同一個(gè)簇。SOM 強(qiáng)調(diào)簇中心點(diǎn)之間的鄰近關(guān)系, 相鄰的簇之間相關(guān)性更強(qiáng), 更有利于解釋結(jié)果, 常用于可視化網(wǎng)絡(luò)數(shù)據(jù)或基因表達(dá)數(shù)據(jù)。

    輸出的結(jié)果“SOM neighbours distances”顯示了到所有鄰居的距離, 這種可視化結(jié)果也稱為U 型矩陣圖(U-matrix)。用不同顏色表示神經(jīng)元的高度, 即相鄰的神經(jīng)元權(quán)值向量之間的距離。顏色越深表示距離越遠(yuǎn), 顏色淺則表示兩個(gè)神經(jīng)元之間的距離近(Kaski et al,1992)。

    1.2 潛在生態(tài)評(píng)價(jià)法

    Hakanson(1980)提出了沉積物潛在生態(tài)危害評(píng)價(jià)方法, 計(jì)算公式如下:

    式中Er為潛在生態(tài)風(fēng)險(xiǎn)參數(shù), Cr為鉈的污染參數(shù), Cx為標(biāo)稱沉積物鉈含量的實(shí)測(cè)值, C0為參比值。目前各國(guó)學(xué)者對(duì)參比值的選擇各不相同, 本文以中國(guó)陸沉積物背景值0.55(解小凡 等, 2015)為參比值。Tr為鉈的毒性響應(yīng)系數(shù), 本文選取的Tr值為40(高博 等, 2008)。

    潛在生態(tài)風(fēng)險(xiǎn)程度分為5 個(gè)等級(jí): 低度、中度、較高、高、極高。

    1.3 樣品采集與處理

    本文以珠江河口與河網(wǎng)中各大支流和出海口門(mén)為研究對(duì)象, 共設(shè)立11 個(gè)研究站位, 采樣站位監(jiān)測(cè)斷面示意圖見(jiàn)圖1。其中, 在西江、北江和東江上分別設(shè)立馬口、三水、石龍南和石龍北等4 個(gè)研究站位, 在西江的主要出海口門(mén)設(shè)立磨刀門(mén)站位, 在東江和北江的出??陂T(mén)設(shè)立虎門(mén)站位和5 個(gè)伶仃洋斷面。

    采樣時(shí)間: 枯水期大潮(2016 年11月30 日—12 月1日); 枯水期小潮(2016 年12 月5 日—12 月6 日); 豐水期大潮(2017 年7 月8 日—7 月9 日); 豐水期小潮(2017 年7 月17 日—7 月18 日)。大小潮期間分別采樣4 次。

    根據(jù)GPS 定位儀準(zhǔn)確定位研究站位監(jiān)測(cè)斷面,使用有機(jī)玻璃采水器于每個(gè)采樣點(diǎn)分表層、中層、底層三層采集水樣, 樣品采集后迅速轉(zhuǎn)移至聚乙烯塑料瓶中加相應(yīng)的保存劑進(jìn)行預(yù)處理; 所有樣品暫存于有冰的專用采樣箱內(nèi), 運(yùn)回實(shí)驗(yàn)室4℃保存并進(jìn)行測(cè)試。

    圖1 珠江河口與河網(wǎng)研究區(qū)域及研究站位置示意圖該圖基于國(guó)家測(cè)繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號(hào)為GS(2019)4342 的標(biāo)準(zhǔn)地圖制作, 底圖無(wú)修改Fig. 1 Study area and sections in Pearl River Estuary and network. The map is based on the standard map GS(2019)4342 downloaded from the website of the Standard Map Service of the State Administration of Surveying,Mapping and Geoinformation, and the base map has not been modified

    重金屬元素測(cè)試所采用的儀器為電感耦合等離子體質(zhì)譜聯(lián)用儀(PE 公司NexION 350D 型)。為保證實(shí)驗(yàn)的準(zhǔn)確性, 每批測(cè)試樣品選取3~5 個(gè)樣品做加標(biāo)回收實(shí)驗(yàn), 向樣品中等量加入標(biāo)準(zhǔn)物質(zhì), 按實(shí)驗(yàn)步驟進(jìn)行檢測(cè), 測(cè)定加標(biāo)回收率, 所得回收率均在90%~110%之間。同時(shí)通過(guò)平行樣方法檢測(cè)實(shí)驗(yàn)的分析精度, 平行測(cè)定樣品數(shù)為每批測(cè)試樣品總數(shù)的20%, 每個(gè)樣品做五組平行樣, 所得相對(duì)標(biāo)準(zhǔn)偏差在5%以內(nèi)(RSD<5%)。

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

    通過(guò)運(yùn)用Rstudio 建立SOM 模型, 輸入鉈濃度數(shù)據(jù), 并輸出聚類分析運(yùn)算結(jié)果。

    2 結(jié)果與討論

    2.1 珠江河口與河網(wǎng)地區(qū)鉈的暴露水平

    通過(guò)對(duì)各斷面鉈的含量數(shù)據(jù)進(jìn)行分析, 將11 個(gè)采樣斷面合并為6 個(gè)斷面, 分別代表珠江河口與河網(wǎng)典型區(qū)域, 得到了6 個(gè)斷面的四個(gè)時(shí)段(枯水期大、小潮和豐水期大、小潮)表、中、底層的最大值、最小值及平均值(表1、2)。所有數(shù)據(jù)檢出率均為100%。

    表1 珠江河口與河網(wǎng)典型區(qū)域枯水期各斷面鉈含量Tab. 1 Thallium contents in typical sections of Pearl River Estuary and network during withered water period

    表2 珠江河口與河網(wǎng)典型區(qū)域豐水期各斷面鉈含量Tab. 2 Thallium contents in typical sections of Pearl River Estuary and network during high water period

    2.1.1 河口區(qū)域

    河口區(qū)域共設(shè)有 7 個(gè)采樣斷面, 分別代表伶仃洋、虎門(mén)和磨刀門(mén)三大河口。由表1 和表2 可知, 伶仃洋斷面鉈平均含量在枯水期大潮為0.016μg·L-1, 小潮為 0.113μg·L-1; 豐水期大潮為0.063μg·L-1, 小潮為0.063μg·L-1?;㈤T(mén)斷面鉈平均含 量 在 枯 水 期 大 潮 為 0.017μg·L-1, 小 潮 為0.016μg·L-1; 豐水期大潮為 0.069μg·L-1, 小潮為0.077μg·L-1。磨刀門(mén)斷面鉈平均含量在枯水期大潮為0.126μg·L-1, 小潮為0.019μg·L-1; 豐水期大潮為0.134μg·L-1, 小潮為0.105μg·L-1。整體而言,伶仃洋在枯水期小潮時(shí)鉈的平均含量超過(guò)了國(guó)家飲用水標(biāo)準(zhǔn) 0.1μg·L-1(以下簡(jiǎn)稱國(guó)家標(biāo)準(zhǔn)); 而磨刀門(mén)除了枯水期小潮時(shí)不超標(biāo)之外, 其余三個(gè)時(shí)段鉈含量均超過(guò)了國(guó)家標(biāo)準(zhǔn)。

    進(jìn)一步分析發(fā)現(xiàn), 伶仃洋枯水期小潮的中層的鉈含量超標(biāo), 最大值達(dá)1.270μg·L-1, 導(dǎo)致這個(gè)時(shí)期的鉈含量超過(guò)國(guó)家標(biāo)準(zhǔn)。磨刀門(mén)枯水期大潮情況與之相似, 中層鉈含量最大為1.302μg·L-1; 而豐水期的大小潮都是因?yàn)檎w的鉈含量均處于較高的水平,故超過(guò)了國(guó)家標(biāo)準(zhǔn)。

    2.1.2 河網(wǎng)上游區(qū)域

    河網(wǎng)區(qū)域共設(shè)有4 個(gè)斷面, 分別為西江的馬口斷面、北江的三水?dāng)嗝嬉约皷|江的石龍北、石龍南斷面, 4 個(gè)斷面均位于研究區(qū)域上游。由表1 可知,馬口斷面鉈的平均含量在枯水期大潮為0.228μg·L-1, 小潮為 0.020μg·L-1; 豐水期大潮為0.082μg·L-1, 小潮為0.086μg·L-1。三水?dāng)嗝驺B的平均 含 量 在 枯 水 期 大 潮 為 0.015μg·L-1, 小 潮 為0.124μg·L-1; 豐水期大潮為 0.108μg·L-1, 小潮為0.026μg·L-1。石龍斷面鉈的平均含量在枯水期大潮為0.068μg·L-1, 小潮為0.191μg·L-1; 豐水期大潮為0.116μg·L-1, 小潮為0.080μg·L-1。其中, 馬口枯水期大潮、三水枯水期小潮和豐水期大潮、石龍枯水期大潮和豐水期小潮鉈含量均未超出國(guó)家標(biāo)準(zhǔn)),其余均在國(guó)家標(biāo)準(zhǔn)限值以下。

    進(jìn)一步分析發(fā)現(xiàn), 枯水期大潮期間馬口和三水兩斷面上層的鉈含量最高分別達(dá)1.320、1.330μg·L-1,從而使這兩個(gè)斷面的鉈含量平均值大大升高。三水?dāng)嗝尕S水期大潮鉈平均含量超過(guò)國(guó)家標(biāo)準(zhǔn)是因?yàn)殂B含量整體均處于較高的水平, 且下層的鉈平均含量特別高。石龍斷面枯水期大潮和豐水期小潮的鉈濃度超標(biāo)同樣是因?yàn)樵搮^(qū)域整體的鉈含量均處于較高的水平。

    2.1.3 河口與河網(wǎng)區(qū)域鉈的暴露特征及原因

    從表1 和表2 可知, 除了豐水期小潮時(shí)鉈的含量為河網(wǎng)區(qū)域小于河口區(qū)域, 其余三個(gè)時(shí)段均為河網(wǎng)區(qū)域大于河口區(qū)域。其原因是在枯水期時(shí), 上游來(lái)水量小, 大潮時(shí)候頂托上溯增加, 此時(shí)河口區(qū)域的鉈向河網(wǎng)區(qū)域輸送, 導(dǎo)致了此時(shí)段河網(wǎng)區(qū)域鉈含量大于河口區(qū)域; 小潮時(shí), 雖然頂托上溯減弱,但上溯能力仍然大于上游來(lái)水輸送能力(胡溪 等,2012; 孔蘭 等, 2011), 因此河口區(qū)域的鉈向河網(wǎng)區(qū)域輸送, 導(dǎo)致了此時(shí)段河網(wǎng)區(qū)域鉈含量大于河口區(qū)域。

    當(dāng)豐水期時(shí), 上游來(lái)水量大, 大潮時(shí)頂托上溯增加, 但上溯能力小于上游來(lái)水輸送能力, 與此同時(shí)外海向河口區(qū)域涌入, 稀釋了河口鉈的濃度。因此, 雖然此時(shí)段鉈從河網(wǎng)區(qū)域向河口區(qū)域輸送, 鉈含量仍然是河網(wǎng)區(qū)域>河口區(qū)域; 而小潮時(shí), 頂托上溯減少, 鉈向河口輸送, 導(dǎo)致了此時(shí)段河網(wǎng)區(qū)域鉈含量小于河口區(qū)域。

    2.1.4 研究區(qū)域與其他區(qū)域?qū)Ρ惹闆r

    本文共獲得監(jiān)測(cè)數(shù)據(jù)528 個(gè), 其中48 個(gè)超過(guò)了國(guó)家飲用水標(biāo)準(zhǔn)(0.1μg·L-1), 超標(biāo)個(gè)數(shù)占總數(shù)的9.09%, 但監(jiān)測(cè)數(shù)據(jù)總體均未超過(guò)廣東省工業(yè)排放標(biāo)準(zhǔn)。表3 列舉了部分國(guó)內(nèi)外天然水體中鉈的含量,在天然水體中鉈的含量一般極低, 相比之下, 珠江河網(wǎng)河口區(qū)域鉈的暴露水平處于較高水平。

    2.2 珠江河口河網(wǎng)鉈的空間分布

    為方便分析珠江河口河網(wǎng)區(qū)域鉈的空間分布情況, 消除極值對(duì)平均數(shù)的影響, 采用幾何平均的方法, 分別計(jì)算各采樣時(shí)段各斷面上、中、下三個(gè)層面的幾何平均值。上、中、下層的關(guān)系和位置以圖2 中的石龍北斷面的柱狀圖為例; 將國(guó)家限值0~0.1μg·L-1平均分成三等, 認(rèn)為0~0.033μg·L-1為低水 平; 0.033~0.066μg·L-1為 中 等 水 平; 0.066~0.1μg·L-1為高水平。

    計(jì)算結(jié)果表明, 每個(gè)斷面的鉈濃度幾何平均值均處于國(guó)家限值之內(nèi)(圖2)。珠江河口與河網(wǎng)鉈的濃度基本處于中等偏下的水平, 個(gè)別斷面?zhèn)€別層面含量突出。河網(wǎng)區(qū)域, 馬口斷面在空間上的鉈分布從高到低為上層、下層、中層, 整體含量為中等偏下水平, 每層差別較大; 三水整體含量為低水平, 每層差距不大; 石龍北為中層、下層、上層, 整體含量為中等偏低水平, 每層差別較為明顯; 石龍南為中層、下層、上層, 整體含量為中等偏上水平, 每層差別較大。河口上, 伶仃洋的5 個(gè)斷面除了斷面1 和斷面4(下面簡(jiǎn)稱1—5), 其余鉈含量相差不大, 基本處于低水平; 伶仃洋1 的鉈分布從高到低為中層、上層和下層差別不大, 整體含量為中等偏下水平;伶仃洋4 情況與1 相似, 只是整體含量比1 低; 虎門(mén)整體含量為低水平, 每層差距不大; 磨刀門(mén)為上層、中層和下層差別不大, 整體含量為中等偏上水平。珠江河口河網(wǎng)鉈的濃度在空間上整體表現(xiàn)為河網(wǎng)區(qū)域>河口區(qū)域; 而河網(wǎng)區(qū)域中, 表現(xiàn)為石龍>馬口>三水; 河口區(qū)域中, 則為磨刀門(mén)>伶仃洋>虎門(mén)。

    圖2 珠江河口河網(wǎng)鉈的空間分布該圖基于國(guó)家測(cè)繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號(hào)為GS(2019)4342 的標(biāo)準(zhǔn)地圖制作, 底圖無(wú)修改Fig. 2 Spatial distribution of thallium in the estuary and network of the Pearl River. The map is based on the standard map GS(2019)4342 downloaded from the website of the Standard Map Service of the State Administration of Surveying, Mapping and Geoinformation, and the base map has not been modified

    在采樣周期內(nèi), 馬口斷面的鉈含量偏高, 推測(cè)原因是西江作為珠江水系的最大干流, 上游來(lái)流量大, 馬口斷面上游分布著許多礦廠, 礦廠排放的含鉈廢水隨流而來(lái), 導(dǎo)致了馬口斷面的鉈的含量較高。此外, 馬口斷面上層的鉈含量較高, 說(shuō)明此處的鉈大多呈離子態(tài)溶于水里, 部分鉈為化合物成為沉積物淤積在底部。造成石龍北、石龍南鉈含量高的原因可能是東江下游河網(wǎng)區(qū)的支柱行業(yè)排放導(dǎo)致的,石龍南所受支柱行業(yè)影響較大, 故石龍南的鉈含量比石龍北高。此外, 該兩處斷面的鉈濃度均為中下兩層濃度較高, 推測(cè)原因是此處鉈為化合態(tài)。伶仃洋1 主要受到周邊島上和岸上的鉈排放影響, 而伶仃洋4 處在珠江河口出??谔? 受外海的高鉈含量污染物和沉積物的影響, 故4 比1 的含量高。磨刀門(mén)存在嚴(yán)重的淤積問(wèn)題和咸潮上溯問(wèn)題, 受工業(yè)污染和生活污水影響, 是珠江八大口門(mén)重金屬污染最為嚴(yán)重的區(qū)域之一, 此處的鉈含量最高, 鉈污染較為嚴(yán)重(宋美英, 2014)。

    2.3 珠江河口與河網(wǎng)區(qū)域鉈的時(shí)間分布特征

    為分析珠江河口與河網(wǎng)鉈的時(shí)間分布情況, 采用同空間分布的方法, 分別計(jì)算出采樣斷面上每個(gè)層次在枯水期大小潮、豐水期大小潮四個(gè)時(shí)期的幾何平均值。結(jié)果發(fā)現(xiàn), 雖然個(gè)別層次鉈的濃度處于較高水平, 但是整體上每個(gè)斷面的計(jì)算結(jié)果都在國(guó)家限值內(nèi), 詳見(jiàn)圖3。珠江河口河網(wǎng)鉈的濃度在時(shí)間從高到低分別為豐水期小潮>豐水期大潮>枯水期大潮>枯水期小潮。

    圖3 珠江河口河網(wǎng)鉈的時(shí)間分布圖a、b、c、d 分別為枯水期大潮、枯水期小潮、豐水期大潮、豐水期小潮; 該圖基于國(guó)家測(cè)繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號(hào)為GS(2019)4342 的標(biāo)準(zhǔn)地圖制作, 底圖無(wú)修改Fig. 3 Temporal distribution of thallium concentration in the estuary and network of the Pearl River. a, b, c, and d are spring tide of withered water period, neap tide of withered water period spring tide of high water period and neap tide of high water period, respectively. The map is based on the standard map GS(2019)4342 downloaded from the website of the Standard Map Service of the State Administration of Surveying, Mapping and Geoinformation, and the base map has not been modified

    2.3.1 珠江河口與河網(wǎng)區(qū)域枯水期鉈的分布

    如圖3a 所示, 枯水期大潮時(shí), 馬口斷面鉈含量水平很高; 其次是石龍北的中層, 為中等偏上的水平; 馬口中層、石龍南上層、伶仃洋2 下層、伶仃洋4 中下層和磨刀門(mén)中上層鉈含量較相似, 為中等偏下; 其余均為較低含量水平。

    圖3b 中, 枯水期小潮時(shí), 石龍北的上中下層面鉈的濃度均處于中等水平; 其次是伶仃洋2 的濃度也較高, 為中等偏下水平; 而馬口、石龍南、三水和磨刀門(mén)的鉈濃度為中等到低水平之間; 其余均為較低含量水平。

    綜合來(lái)看, 枯水期河口區(qū)的鉈污染小于河網(wǎng)區(qū)。這兩個(gè)時(shí)期的磨刀門(mén)上層和中層鉈濃度含量較高, 結(jié)合馬口上層鉈含量較高的情況, 可推斷枯水期時(shí)西江流域上的來(lái)沙量不多, 且此時(shí)磨刀門(mén)的泥沙淤積問(wèn)題不嚴(yán)重, 此時(shí)水體中的鉈多為離子態(tài)存在。而石龍北在枯水期的鉈含量高, 可見(jiàn)東江下游河網(wǎng)區(qū)的工業(yè)生產(chǎn)和排放影響較大。而石龍北斷面的鉈在中下層含量多, 此處的沉積物容易堆積, 推測(cè)鉈主要賦存在沉積物中。一般情況下, 大潮時(shí)海水涌入, 水體污染物濃度應(yīng)該降低, 小潮時(shí)陸源污染進(jìn)入, 水體污染物濃度升高。然而, 伶仃洋4 在大潮的鉈濃度卻比小潮時(shí)候高, 推測(cè)是因?yàn)榭菟跁r(shí)候外海的鉈濃度較高, 大潮時(shí)海水進(jìn)入伶仃洋造成此處鉈濃度升高, 而小潮時(shí)上游鉈濃度較少的水體輸入稀釋了此處的鉈。

    結(jié)合水動(dòng)力特征, 枯水期于枯季徑流微弱, 珠江河網(wǎng)區(qū)為潮汐控制, 潮水填充河道, 潮區(qū)界、潮流界可達(dá)三水、馬口以上(朱金格 等, 2010), 故河口的鉈濃度整體比河網(wǎng)低。而西江河網(wǎng)的潮汐傳播速度大于北江河網(wǎng)(陳小齊 等, 2020), 故地理位置相鄰的三水的鉈濃度比馬口小; 而東江河網(wǎng)上的石龍南、北斷面距離河口區(qū)域較近, 鉈隨潮汐傳播上溯所需時(shí)間較馬口和三水短, 故在同一采樣時(shí)間內(nèi),此處濃度較馬口和三水高。

    2.3.2 珠江河口與河網(wǎng)區(qū)域豐水期鉈的分布

    圖3c 中可看出, 豐水期大潮時(shí), 采樣斷面整體鉈的濃度都升高, 為中等偏上水平。此時(shí)伶仃洋口門(mén)區(qū)5 個(gè)斷面和虎門(mén)的鉈含量較其他斷面低; 鉈濃度最高的斷面是磨刀門(mén)。而在豐水期小潮時(shí), 除了三水的鉈濃度較低以外, 其他斷面的鉈含量都偏高,為中等偏上, 濃度最高的為磨刀門(mén)(圖3d)。

    整體上豐水期的鉈濃度都處于較高的水平, 大潮時(shí)河口區(qū)的鉈污染小區(qū)河網(wǎng)區(qū)。磨刀門(mén)的鉈濃度在豐水期大小潮都很高, 且上層和下層均出現(xiàn)最高濃度。豐水期時(shí)候西江不僅上游流量大, 且含沙量也不小, 在磨刀門(mén)的淤積問(wèn)題嚴(yán)重, 推測(cè)此時(shí)磨刀門(mén)斷面上鉈不僅以游離態(tài)存在水中, 還有大量的鉈與沉積物存在。河口處的伶仃洋口門(mén)區(qū)和虎門(mén)在大潮時(shí)由于海水涌入稀釋而使鉈濃度降低, 而小潮時(shí)陸源鉈污染進(jìn)入而導(dǎo)致鉈濃度升高。至于三水和馬口, 雖然它們?cè)诘乩砦恢蒙舷嘟? 但在大潮時(shí)三水的鉈污染較馬口嚴(yán)重, 小潮時(shí)截然相反。馬口大小潮時(shí)鉈濃度差別不大, 大潮略低于小潮, 與其他斷面上的鉈時(shí)間分布規(guī)律大致相同, 而三水可能受到人為干預(yù)的影響較大, 在小潮時(shí)候北江上游的鉈排放減少, 導(dǎo)致此時(shí)三水鉈含量偏低。

    同樣結(jié)合河口河網(wǎng)水動(dòng)力特征進(jìn)行分析, 豐水期三角洲河道整體行洪能力增強(qiáng), 地形下切引起珠江三角洲河網(wǎng)內(nèi)潮動(dòng)力增強(qiáng), 而下游河段徑流受阻(陳小齊 等, 2020), 導(dǎo)致河網(wǎng)內(nèi)潮差增大, 引起河網(wǎng)上游水位雍高, 故大潮時(shí)河網(wǎng)的鉈濃度大于河口;小潮時(shí)潮汐動(dòng)力減弱, 河網(wǎng)內(nèi)潮差較小, 上游流量順利下排, 故此時(shí)河口的鉈濃度大于河口。

    2.4 SOM 模型分析結(jié)果

    文中采用 Rstudio 平臺(tái)進(jìn)行建模及聚類分析,采用R 語(yǔ)言自帶的Kohonen 軟件包對(duì)528 個(gè)鉈濃度進(jìn)行帶有基因表達(dá)的SOM 分析, 所得結(jié)果如圖4a所示。根據(jù)潛在生態(tài)評(píng)價(jià)法, 按照公式和參數(shù)對(duì)鉈濃度數(shù)據(jù)進(jìn)行處理, 運(yùn)用SOM 模型進(jìn)行分析, 最后得到結(jié)果如圖4b 所示。顏色深淺表示神經(jīng)元之間的遠(yuǎn)近距離, 同樣亦表示鉈濃度高低之間的遠(yuǎn)近距離。顏色越紅即數(shù)據(jù)之間的差距越大, 濃度越高; 顏色越藍(lán)表示數(shù)據(jù)之間的差距越小, 含量越小。生態(tài)風(fēng)險(xiǎn)聚類結(jié)果同理。鉈濃度和潛在生態(tài)風(fēng)險(xiǎn)聚類結(jié)果表見(jiàn)表4。

    如圖4a 所示, 磨刀門(mén)和石龍南下層鉈的含量最高, 其次為石龍南下層、馬口上層、磨刀門(mén)中層、馬口中層、磨刀門(mén)上層, 其余各斷面的含量均處于中等偏下至較低水平。SOM 結(jié)果與幾何平均值計(jì)算結(jié)果(圖2)存在些許差別, 但是大體上它們均表現(xiàn)出一致的空間分布特征規(guī)律。兩種結(jié)果產(chǎn)生差別的原因可能由幾何平均值統(tǒng)計(jì)出結(jié)果的受到干擾而SOM 模型結(jié)果是直接根據(jù)528 個(gè)數(shù)據(jù)計(jì)算, 受到干擾較少。

    圖4 SOM 基因表達(dá)聚類結(jié)果Fig. 4 Clustering results of SOM gene expression. Left and right are clustering results of concentration and potential ecological risks of thallium, respectively

    表4 鉈濃度和潛在生態(tài)風(fēng)險(xiǎn)聚類結(jié)果表Tab. 4 Clustering results of concentration and potential ecological risks of thallium

    如圖4b 所示, 珠江河網(wǎng)河口有5 個(gè)站位水層的潛在生態(tài)風(fēng)險(xiǎn)程度處于極高風(fēng)險(xiǎn), 它們分別是磨刀門(mén)上層、馬口上層、磨刀門(mén)中層、馬口中層、石龍南下層; 處于高風(fēng)險(xiǎn)程度的站位水層有4 個(gè), 分別是石龍南上中層、磨刀門(mén)下層、石龍北下層; 處于較高風(fēng)險(xiǎn)程度的有2 個(gè), 為石龍北上、中層; 處于中度風(fēng)險(xiǎn)的有3 個(gè), 為三水下層、伶仃洋5 下層和馬口下層; 其余點(diǎn)位均處于低度生態(tài)風(fēng)險(xiǎn)。

    綜上所述, 磨刀門(mén)、馬口和石龍南、石龍北的鉈含量和潛在生態(tài)風(fēng)險(xiǎn)程度較其他站點(diǎn)高(但都未超過(guò)國(guó)家標(biāo)準(zhǔn))。磨刀門(mén)的鉈含量最高且潛在生態(tài)風(fēng)險(xiǎn)程度也極高, 這主要受到人類活動(dòng)和工業(yè)排放的影響且與磨刀門(mén)的淤積和咸潮上溯問(wèn)題有關(guān)。馬口的鉈含量和潛在生態(tài)風(fēng)險(xiǎn)也很高, 該處的鉈主要是西江上游帶來(lái)的。石龍北的鉈污染情況比石龍南輕,由此推斷東江上游河網(wǎng)排放的鉈含量較高。

    SOM 基因表達(dá)聚類分析是一種無(wú)監(jiān)督學(xué)習(xí)的模式, 不受人為因素干擾, 應(yīng)用于污染物聚類分析是合適的。尤其是像鉈這種毒性大, 靈敏度高, 且自然水體中的鉈含量閾值低的重金屬, 分析時(shí)細(xì)微誤差可能就導(dǎo)致結(jié)果千差萬(wàn)別。SOM 聚類結(jié)果顯示,在只輸入鉈含量數(shù)據(jù)的情況下, 不設(shè)置其他干擾參數(shù), 模型結(jié)果仍然能夠準(zhǔn)確輸出鉈含量和潛在生態(tài)風(fēng)險(xiǎn)程度的高低區(qū)域。雖然SOM 結(jié)果與數(shù)理統(tǒng)計(jì)結(jié)果些許差別, 但是整體上所呈現(xiàn)出的規(guī)律是一致的,證明SOM 是一種可取的污染物聚類分析方法。

    3 結(jié)論

    1) 本研究所有站點(diǎn)的鉈濃度均未超過(guò)廣東省工業(yè)排放標(biāo)準(zhǔn), 表明整體上珠江河口河網(wǎng)的鉈含量處于較低水平, 鉈的污染情況不嚴(yán)重。

    2) 珠江河口河網(wǎng)的鉈含量時(shí)間分布特征為豐水期小潮>豐水期大潮>枯水期大潮>枯水期小潮,且除了豐水期小潮時(shí)段外, 其余水文時(shí)段河網(wǎng)區(qū)域的鉈含量均大于河口區(qū)域; 而珠江河口河網(wǎng)鉈濃度的空間分布特征為磨刀門(mén)>石龍>馬口>伶仃洋>虎門(mén)>三水。

    3) SOM 聚類分析結(jié)果表明, 珠江河口河網(wǎng)上有4 個(gè)斷面具有較高的潛在生態(tài)風(fēng)險(xiǎn), 其程度高低依次為磨刀門(mén)>馬口>石龍南>石龍北, 而其他斷面的潛在生態(tài)風(fēng)險(xiǎn)程度均處于低的水平。

    4) 運(yùn)用SOM 模型對(duì)水體中的鉈進(jìn)行聚類分析,能夠在最大程度上排除人為干預(yù)、統(tǒng)計(jì)分析等各種誤差, 因此是一種合適的污染物聚類分析方法。

    猜你喜歡
    小潮馬口石龍
    馬口魚(yú)
    垂釣(2023年11期)2024-01-21 16:07:04
    非遺視野下湖北馬口窯的保護(hù)、繼承與開(kāi)發(fā)
    流行色(2019年5期)2019-12-13 18:38:50
    鷓鴣天·游石龍水庫(kù)
    力量哪里去了
    十月(2019年4期)2019-07-29 17:33:30
    曲靖市檔案學(xué)會(huì)傾囊相助 力推石龍村脫貧出列
    云南檔案(2019年5期)2019-06-20 09:16:36
    希 望
    希望
    新一季流行色已發(fā)布?快來(lái)入手同色系數(shù)碼3C小潮物!
    來(lái)自馬口窯的對(duì)話
    ——馬口窯文獻(xiàn)與當(dāng)代陶藝創(chuàng)作研究展
    馬口煤礦分層開(kāi)采工作面上覆采空區(qū)自燃治理技術(shù)研究
    亚洲一区二区三区欧美精品| 亚洲第一区二区三区不卡| 99国产精品免费福利视频| 国产色婷婷99| 在线亚洲精品国产二区图片欧美| 色播在线永久视频| 黄片播放在线免费| 欧美精品av麻豆av| 免费不卡黄色视频| 亚洲国产av影院在线观看| 欧美日韩成人在线一区二区| av.在线天堂| 国产午夜精品一二区理论片| 精品亚洲成a人片在线观看| 丰满饥渴人妻一区二区三| 亚洲精品国产av成人精品| 热re99久久国产66热| 我的亚洲天堂| 国产激情久久老熟女| 日本vs欧美在线观看视频| 99精品久久久久人妻精品| 男的添女的下面高潮视频| 久久精品国产亚洲av涩爱| 欧美黑人精品巨大| 少妇精品久久久久久久| 免费黄网站久久成人精品| 99热全是精品| 99久国产av精品国产电影| 国产精品女同一区二区软件| 中文字幕最新亚洲高清| 熟女少妇亚洲综合色aaa.| 色吧在线观看| 精品一区二区免费观看| 亚洲国产av新网站| 精品国产一区二区久久| 一区在线观看完整版| 大香蕉久久成人网| 18在线观看网站| 天天操日日干夜夜撸| 亚洲成人av在线免费| 国产精品一区二区在线观看99| 日韩熟女老妇一区二区性免费视频| 性高湖久久久久久久久免费观看| 美女高潮到喷水免费观看| 亚洲精品久久成人aⅴ小说| 男人添女人高潮全过程视频| 欧美97在线视频| 国产精品麻豆人妻色哟哟久久| 久久久国产欧美日韩av| 三上悠亚av全集在线观看| 一级黄片播放器| 最近2019中文字幕mv第一页| 老司机在亚洲福利影院| 欧美精品亚洲一区二区| 一本久久精品| av电影中文网址| 国产成人a∨麻豆精品| 色播在线永久视频| 亚洲精品国产av成人精品| 热99国产精品久久久久久7| av国产久精品久网站免费入址| 黄色怎么调成土黄色| 国精品久久久久久国模美| 午夜福利影视在线免费观看| 欧美人与性动交α欧美软件| 老司机亚洲免费影院| 日本欧美国产在线视频| 国产黄色视频一区二区在线观看| 国产 精品1| 欧美日韩av久久| 免费黄色在线免费观看| 亚洲精品中文字幕在线视频| 中文字幕亚洲精品专区| 国产亚洲av片在线观看秒播厂| 久久精品熟女亚洲av麻豆精品| 悠悠久久av| 熟妇人妻不卡中文字幕| 熟妇人妻不卡中文字幕| 婷婷成人精品国产| 如何舔出高潮| 亚洲国产成人一精品久久久| tube8黄色片| 久久精品国产a三级三级三级| 女的被弄到高潮叫床怎么办| 国产爽快片一区二区三区| 欧美 日韩 精品 国产| 国产深夜福利视频在线观看| 国产欧美日韩一区二区三区在线| 十八禁高潮呻吟视频| 午夜影院在线不卡| 在线观看人妻少妇| 久久精品aⅴ一区二区三区四区| 亚洲第一青青草原| 欧美黑人欧美精品刺激| 黑丝袜美女国产一区| 国产免费福利视频在线观看| 在线看a的网站| 天天躁夜夜躁狠狠久久av| 色精品久久人妻99蜜桃| 免费人妻精品一区二区三区视频| 国产欧美日韩综合在线一区二区| 亚洲av成人精品一二三区| 午夜老司机福利片| 中文精品一卡2卡3卡4更新| 久久综合国产亚洲精品| 亚洲成人手机| 美女国产高潮福利片在线看| 一级a爱视频在线免费观看| 亚洲精品中文字幕在线视频| 国产在线一区二区三区精| 久久精品aⅴ一区二区三区四区| 久久久欧美国产精品| 看免费成人av毛片| 精品国产露脸久久av麻豆| 满18在线观看网站| av在线观看视频网站免费| 日韩精品有码人妻一区| 成年女人毛片免费观看观看9 | 麻豆精品久久久久久蜜桃| 18禁动态无遮挡网站| 亚洲图色成人| 97在线人人人人妻| 久久精品aⅴ一区二区三区四区| 五月开心婷婷网| 香蕉丝袜av| 天堂中文最新版在线下载| 在线精品无人区一区二区三| 免费在线观看完整版高清| 在线观看免费视频网站a站| 国产伦理片在线播放av一区| 国产欧美亚洲国产| 在线天堂最新版资源| 最新的欧美精品一区二区| 亚洲熟女精品中文字幕| 18禁裸乳无遮挡动漫免费视频| 丰满乱子伦码专区| 久久99热这里只频精品6学生| 嫩草影院入口| 美女国产高潮福利片在线看| 亚洲精品在线美女| 韩国精品一区二区三区| 国产亚洲av片在线观看秒播厂| 亚洲人成网站在线观看播放| 99国产精品免费福利视频| 久久精品人人爽人人爽视色| 男的添女的下面高潮视频| 日本色播在线视频| 色婷婷av一区二区三区视频| 亚洲av福利一区| 制服人妻中文乱码| xxxhd国产人妻xxx| 亚洲精品国产av成人精品| 亚洲欧美激情在线| 午夜福利视频在线观看免费| 国产亚洲av片在线观看秒播厂| 午夜免费观看性视频| 99国产精品免费福利视频| 欧美日韩综合久久久久久| a级毛片在线看网站| 在线观看免费日韩欧美大片| 日本av免费视频播放| 操美女的视频在线观看| 成人漫画全彩无遮挡| 在现免费观看毛片| 秋霞在线观看毛片| 黑人巨大精品欧美一区二区蜜桃| 亚洲第一av免费看| 免费看不卡的av| 国产高清不卡午夜福利| 欧美精品一区二区免费开放| 国产成人免费观看mmmm| 一级片'在线观看视频| 人人妻人人添人人爽欧美一区卜| 国语对白做爰xxxⅹ性视频网站| 日韩一区二区三区影片| 亚洲少妇的诱惑av| 国产成人一区二区在线| 欧美日韩亚洲综合一区二区三区_| 在线天堂中文资源库| av有码第一页| 亚洲精品日韩在线中文字幕| 十分钟在线观看高清视频www| 天美传媒精品一区二区| 亚洲精品日本国产第一区| 欧美成人精品欧美一级黄| 免费少妇av软件| 亚洲成色77777| 欧美xxⅹ黑人| 在线观看免费日韩欧美大片| 亚洲国产欧美在线一区| 免费观看性生交大片5| 国产黄色视频一区二区在线观看| av在线app专区| 黄色视频不卡| 日韩人妻精品一区2区三区| 性少妇av在线| 久久人人爽av亚洲精品天堂| www日本在线高清视频| 国产精品蜜桃在线观看| 国产一卡二卡三卡精品 | 波多野结衣一区麻豆| 久久影院123| 五月天丁香电影| 成年动漫av网址| 成年人免费黄色播放视频| av在线观看视频网站免费| 在线观看三级黄色| 精品国产国语对白av| 欧美精品亚洲一区二区| 亚洲国产精品国产精品| 精品少妇黑人巨大在线播放| 亚洲一码二码三码区别大吗| 日本wwww免费看| 午夜免费鲁丝| 国产精品免费视频内射| 午夜影院在线不卡| 国产精品一区二区在线观看99| 波多野结衣一区麻豆| 久久热在线av| 中文乱码字字幕精品一区二区三区| 久久久久久久国产电影| 久久毛片免费看一区二区三区| 亚洲欧美一区二区三区黑人| 成人国产麻豆网| 亚洲国产精品一区三区| 亚洲av电影在线观看一区二区三区| 一本色道久久久久久精品综合| 国产日韩欧美在线精品| 在线观看免费午夜福利视频| 精品久久久精品久久久| 男女免费视频国产| 一级片免费观看大全| 国产精品无大码| 两个人看的免费小视频| 一边摸一边抽搐一进一出视频| 国产深夜福利视频在线观看| kizo精华| 中文天堂在线官网| 国产男人的电影天堂91| 好男人视频免费观看在线| 免费黄网站久久成人精品| 大香蕉久久网| 在线观看免费日韩欧美大片| 美女扒开内裤让男人捅视频| 亚洲国产精品成人久久小说| 国产成人系列免费观看| 老司机影院毛片| 日韩不卡一区二区三区视频在线| 亚洲精品久久午夜乱码| 午夜激情av网站| 精品亚洲成国产av| 19禁男女啪啪无遮挡网站| 大香蕉久久成人网| 欧美精品av麻豆av| 国产 精品1| 日韩电影二区| 中文字幕人妻丝袜制服| 七月丁香在线播放| 日日爽夜夜爽网站| 亚洲精品国产区一区二| 日本欧美视频一区| 亚洲七黄色美女视频| 日韩不卡一区二区三区视频在线| 热re99久久国产66热| 亚洲美女搞黄在线观看| 久久久久久久久久久免费av| 国产极品粉嫩免费观看在线| 啦啦啦在线免费观看视频4| 18禁动态无遮挡网站| 亚洲精品美女久久久久99蜜臀 | 99国产综合亚洲精品| 99香蕉大伊视频| 青春草亚洲视频在线观看| 国产亚洲av高清不卡| 久久久欧美国产精品| av在线观看视频网站免费| 国产人伦9x9x在线观看| 少妇 在线观看| 亚洲第一av免费看| 欧美日韩亚洲综合一区二区三区_| 中文字幕人妻熟女乱码| 男女边吃奶边做爰视频| 午夜福利影视在线免费观看| 91成人精品电影| 久久人妻熟女aⅴ| av网站在线播放免费| av女优亚洲男人天堂| 精品亚洲成a人片在线观看| 成人免费观看视频高清| 亚洲欧美成人综合另类久久久| 久久 成人 亚洲| 久久精品国产a三级三级三级| 国产成人精品在线电影| 亚洲熟女精品中文字幕| 日韩成人av中文字幕在线观看| 亚洲av综合色区一区| av在线老鸭窝| 国产在线视频一区二区| 日本色播在线视频| 国产av国产精品国产| 亚洲美女搞黄在线观看| 亚洲国产日韩一区二区| 看非洲黑人一级黄片| 高清不卡的av网站| 国产乱人偷精品视频| 欧美日韩一区二区视频在线观看视频在线| 国产极品粉嫩免费观看在线| 最近中文字幕高清免费大全6| 18禁动态无遮挡网站| 五月开心婷婷网| 免费不卡黄色视频| 超碰成人久久| 久久国产精品男人的天堂亚洲| 国产福利在线免费观看视频| 国产男女内射视频| 亚洲精品成人av观看孕妇| 熟女av电影| 日本午夜av视频| 99久久99久久久精品蜜桃| 国产又爽黄色视频| 亚洲欧美一区二区三区黑人| 大片电影免费在线观看免费| 九色亚洲精品在线播放| 久久天堂一区二区三区四区| a级毛片在线看网站| 亚洲av电影在线观看一区二区三区| 建设人人有责人人尽责人人享有的| 国产精品久久久人人做人人爽| 国产精品成人在线| 91老司机精品| 纯流量卡能插随身wifi吗| 看十八女毛片水多多多| 精品久久久久久电影网| 国语对白做爰xxxⅹ性视频网站| 午夜影院在线不卡| av在线观看视频网站免费| 久久99热这里只频精品6学生| 免费黄网站久久成人精品| 精品国产一区二区三区久久久樱花| 桃花免费在线播放| 曰老女人黄片| 女性生殖器流出的白浆| 日本猛色少妇xxxxx猛交久久| 一区二区三区激情视频| 欧美日韩综合久久久久久| 国产成人欧美在线观看 | 国产欧美日韩综合在线一区二区| 大片免费播放器 马上看| 国产日韩欧美在线精品| 亚洲国产欧美一区二区综合| 国产成人啪精品午夜网站| 侵犯人妻中文字幕一二三四区| 九草在线视频观看| 国产97色在线日韩免费| 亚洲欧美激情在线| 亚洲第一区二区三区不卡| 亚洲精品久久午夜乱码| 午夜免费观看性视频| 波多野结衣一区麻豆| 欧美日韩精品网址| 亚洲人成77777在线视频| 丰满迷人的少妇在线观看| 日韩制服丝袜自拍偷拍| 国产免费福利视频在线观看| 国产亚洲午夜精品一区二区久久| 飞空精品影院首页| 国产一区二区三区综合在线观看| 91国产中文字幕| 看免费av毛片| 老熟女久久久| 免费女性裸体啪啪无遮挡网站| 色综合欧美亚洲国产小说| 亚洲四区av| 欧美激情 高清一区二区三区| av在线观看视频网站免费| 国产精品女同一区二区软件| 婷婷成人精品国产| 天天添夜夜摸| 亚洲精品日韩在线中文字幕| 亚洲欧美成人综合另类久久久| 五月天丁香电影| 久久毛片免费看一区二区三区| 中文字幕人妻熟女乱码| 无限看片的www在线观看| 国产一区有黄有色的免费视频| 日韩大码丰满熟妇| 丁香六月欧美| avwww免费| 一区在线观看完整版| xxxhd国产人妻xxx| 女人高潮潮喷娇喘18禁视频| 亚洲一码二码三码区别大吗| 两个人看的免费小视频| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久久久久| 国产男女内射视频| 美女视频免费永久观看网站| 精品国产露脸久久av麻豆| 日本猛色少妇xxxxx猛交久久| 国产亚洲av高清不卡| 国产一区二区三区综合在线观看| 日本wwww免费看| 9热在线视频观看99| 国产在线视频一区二区| 色婷婷av一区二区三区视频| 最近手机中文字幕大全| 一区二区日韩欧美中文字幕| 一级,二级,三级黄色视频| 中国国产av一级| 亚洲专区中文字幕在线 | 精品午夜福利在线看| 国产日韩欧美视频二区| 不卡视频在线观看欧美| 观看av在线不卡| 国产精品久久久久成人av| 精品视频人人做人人爽| 亚洲精品自拍成人| 中文字幕高清在线视频| h视频一区二区三区| 99精品久久久久人妻精品| 中文字幕精品免费在线观看视频| 亚洲精品久久成人aⅴ小说| 欧美日韩成人在线一区二区| 亚洲婷婷狠狠爱综合网| 亚洲精品久久成人aⅴ小说| 欧美日韩视频精品一区| 9色porny在线观看| 丁香六月天网| 少妇人妻精品综合一区二区| 岛国毛片在线播放| 在线观看三级黄色| 极品少妇高潮喷水抽搐| 最新的欧美精品一区二区| 国产精品久久久久成人av| 十八禁网站网址无遮挡| 女人精品久久久久毛片| 一级a爱视频在线免费观看| 国产在线免费精品| 天堂8中文在线网| 高清在线视频一区二区三区| 久久久亚洲精品成人影院| 中文字幕av电影在线播放| 国产成人午夜福利电影在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美日韩另类电影网站| 国产精品.久久久| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成人av在线免费| 哪个播放器可以免费观看大片| 最近手机中文字幕大全| 91精品三级在线观看| 日本猛色少妇xxxxx猛交久久| 999精品在线视频| 亚洲精品aⅴ在线观看| 国产成人免费无遮挡视频| 亚洲,欧美精品.| 大码成人一级视频| 纯流量卡能插随身wifi吗| 天天操日日干夜夜撸| 大片免费播放器 马上看| 午夜福利,免费看| 男女边吃奶边做爰视频| 国产野战对白在线观看| 青草久久国产| 五月天丁香电影| 国产av精品麻豆| 国产日韩一区二区三区精品不卡| 欧美97在线视频| 国产麻豆69| 黄色 视频免费看| 又大又爽又粗| 婷婷色综合www| 男女床上黄色一级片免费看| 日韩一卡2卡3卡4卡2021年| 欧美日韩av久久| 丝袜美腿诱惑在线| 精品第一国产精品| 美女中出高潮动态图| 日本vs欧美在线观看视频| 国产日韩欧美视频二区| 亚洲欧美一区二区三区黑人| 亚洲精品av麻豆狂野| 一区二区日韩欧美中文字幕| 操美女的视频在线观看| 婷婷成人精品国产| 亚洲精品日韩在线中文字幕| 欧美日韩av久久| 精品国产乱码久久久久久男人| 在线观看免费高清a一片| 老汉色av国产亚洲站长工具| 免费av中文字幕在线| 一区福利在线观看| av又黄又爽大尺度在线免费看| 亚洲av在线观看美女高潮| 男女边摸边吃奶| 国语对白做爰xxxⅹ性视频网站| 中文字幕人妻丝袜一区二区 | 国产野战对白在线观看| 日韩精品有码人妻一区| 婷婷色av中文字幕| 日韩人妻精品一区2区三区| 热99国产精品久久久久久7| 色网站视频免费| 国产成人a∨麻豆精品| 亚洲国产日韩一区二区| 丰满乱子伦码专区| 在线观看免费视频网站a站| 中文精品一卡2卡3卡4更新| 亚洲欧美精品综合一区二区三区| 国产成人精品久久二区二区91 | 亚洲国产毛片av蜜桃av| 亚洲欧洲国产日韩| 国产欧美日韩一区二区三区在线| 高清黄色对白视频在线免费看| 欧美乱码精品一区二区三区| 成人亚洲精品一区在线观看| 国产成人免费无遮挡视频| 国产一区二区三区综合在线观看| 天堂8中文在线网| av视频免费观看在线观看| 国产成人啪精品午夜网站| 亚洲国产看品久久| 又粗又硬又长又爽又黄的视频| 波多野结衣av一区二区av| 亚洲av男天堂| 国产日韩欧美亚洲二区| 午夜老司机福利片| 在线亚洲精品国产二区图片欧美| 亚洲色图 男人天堂 中文字幕| 啦啦啦啦在线视频资源| av在线观看视频网站免费| 国产精品国产三级国产专区5o| 亚洲av欧美aⅴ国产| 久久久久精品性色| 制服丝袜香蕉在线| 国产老妇伦熟女老妇高清| 中文字幕人妻丝袜制服| 中文字幕精品免费在线观看视频| 精品国产一区二区三区久久久樱花| 亚洲婷婷狠狠爱综合网| 人人妻人人澡人人爽人人夜夜| 美女福利国产在线| 夫妻性生交免费视频一级片| 中文字幕高清在线视频| www.av在线官网国产| 天美传媒精品一区二区| 亚洲精品美女久久久久99蜜臀 | 黄片播放在线免费| 宅男免费午夜| 国产视频首页在线观看| 两个人免费观看高清视频| 2018国产大陆天天弄谢| 国产精品香港三级国产av潘金莲 | 亚洲欧洲精品一区二区精品久久久 | 丝瓜视频免费看黄片| 在线观看免费午夜福利视频| www.自偷自拍.com| 亚洲av男天堂| 欧美精品高潮呻吟av久久| av网站在线播放免费| 一区福利在线观看| 欧美另类一区| 赤兔流量卡办理| 国产视频首页在线观看| 亚洲,欧美精品.| 校园人妻丝袜中文字幕| 中文字幕人妻丝袜一区二区 | 亚洲人成网站在线观看播放| 男女高潮啪啪啪动态图| 人妻 亚洲 视频| 国产成人午夜福利电影在线观看| 多毛熟女@视频| 天天影视国产精品| 一级毛片 在线播放| www.自偷自拍.com| 最近最新中文字幕免费大全7| av在线老鸭窝| 大片电影免费在线观看免费| 午夜福利乱码中文字幕| 青春草亚洲视频在线观看| 欧美黑人欧美精品刺激| 亚洲 欧美一区二区三区| 国产欧美亚洲国产| 精品国产国语对白av| 免费高清在线观看日韩| 精品卡一卡二卡四卡免费| 18禁国产床啪视频网站| 男女高潮啪啪啪动态图| 久久天堂一区二区三区四区| 亚洲国产精品一区三区| 国产在线视频一区二区| 丝袜美足系列| 51午夜福利影视在线观看| 日韩 亚洲 欧美在线| 亚洲精品久久久久久婷婷小说| 美女视频免费永久观看网站| 久久精品久久久久久噜噜老黄| 精品亚洲成a人片在线观看| 成人三级做爰电影| h视频一区二区三区| 亚洲国产精品一区三区| 伦理电影大哥的女人| 咕卡用的链子| 午夜福利网站1000一区二区三区| 97精品久久久久久久久久精品| 亚洲熟女精品中文字幕| 亚洲国产精品国产精品| 色综合欧美亚洲国产小说| 蜜桃在线观看..| 国产av国产精品国产| 国产av码专区亚洲av| 亚洲国产精品成人久久小说| 亚洲精品,欧美精品| 黑人巨大精品欧美一区二区蜜桃| 日韩中文字幕视频在线看片|