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

    山地城市主要河流水質(zhì)評價及預(yù)測研究
    ——以四川省綿陽市為例

    2023-11-04 14:23:48王源哲華春林趙麗樊敏梁曉盈周樂樂蔡璨姚婧
    生態(tài)環(huán)境學(xué)報 2023年8期
    關(guān)鍵詞:斷面水質(zhì)污染

    王源哲,華春林,趙麗,樊敏,梁曉盈,周樂樂,蔡璨,姚婧

    1. 西南科技大學(xué)環(huán)境與資源學(xué)院,四川 綿陽 621010;2. 西南科技大學(xué)經(jīng)濟(jì)管理學(xué)院,四川 綿陽 621010;3. 西南科技大學(xué)成都創(chuàng)新研究院,四川 成都 621010

    水是生命之源,是人類生活與社會生產(chǎn)的重要物質(zhì)基礎(chǔ),然而,隨著社會經(jīng)濟(jì)的迅猛發(fā)展,水體水質(zhì)呈現(xiàn)出污染程度各異、污染機(jī)理復(fù)雜的趨勢,水質(zhì)污染已成為制約經(jīng)濟(jì)社會可持續(xù)發(fā)展的主要因素(馬樂寬等,2013)。在水質(zhì)遭受不同程度污染的情景下,對水質(zhì)進(jìn)行綜合評價,確定主要污染物,把握水體水質(zhì)健康及污染程度,進(jìn)而因地制宜地采取精準(zhǔn)治污措施,此外,科學(xué)預(yù)測水質(zhì)變化趨勢,有利于促進(jìn)有效長期的區(qū)域水環(huán)境保護(hù)規(guī)劃的制定(陳潤羊等,2008;高紅杰等,2017)。黨的二十大報告明確提出,我們要統(tǒng)籌水資源、水環(huán)境、水生態(tài)治理,推動重要江河湖庫生態(tài)保護(hù)治理。因此,為了對水資源進(jìn)行合理的開發(fā)和科學(xué)的管理,有必要對水環(huán)境進(jìn)行評價和預(yù)測。

    目前,國內(nèi)外常見的水質(zhì)評價方法包括單因子污染指數(shù)評價法、綜合污染指數(shù)評價法(Ji et al.,2016;Ouyang et al.,2022)、灰色系統(tǒng)評價法(張冉等,2013)、主成分分析法(朱琳等,2018)以及模糊綜合評價法(杜軍凱等,2015;楊浩等,2016)。其中,主成分分析法可明確水質(zhì)的主要影響因素與水質(zhì)變化規(guī)律及差異(蔡廣強(qiáng)等,2018),模糊綜合評價法能夠定量化地處理水環(huán)境功能區(qū)劃及水質(zhì)分級界限具有模糊性的問題(靳會姣等,2014),改進(jìn)的主成分-模糊綜合評價耦合模型充分考慮了水質(zhì)和水量的影響,不僅可以選擇關(guān)鍵指標(biāo),還可以指定水質(zhì)的空間變化和等級(Zhao et al.,2020)。然而,現(xiàn)階段的水質(zhì)評價方法主要側(cè)重于處理評價指標(biāo)以及解決水環(huán)境的不確定性問題,關(guān)于水質(zhì)綜合評價方法的選取的研究尚且不足,并且評價結(jié)果是否科學(xué)合理尚無確切的定義。同時,水質(zhì)預(yù)測作為掌握水環(huán)境的變化趨勢及預(yù)防水質(zhì)惡化的前提與基礎(chǔ),對于及時了解水質(zhì)狀況、預(yù)防水污染、保護(hù)水環(huán)境具有重要的現(xiàn)實意義。常用的水質(zhì)預(yù)測方法有支持向量機(jī)預(yù)測(Vapnik et al.,1995)、水質(zhì)模擬模型(Rani et al.,2010)、混沌理論(Qiao et al.,2016)和神經(jīng)網(wǎng)絡(luò)模型(Zheng et al.,2022)等?,F(xiàn)階段,國內(nèi)外研究人員已將支持向量機(jī)運用到了多種水質(zhì)指標(biāo)的預(yù)測中(Noori et al,2015),支持向量機(jī)預(yù)測法雖然可以解決很多傳統(tǒng)方法解決不了的數(shù)學(xué)問題,但在水質(zhì)預(yù)測中存在支持向量數(shù)多,預(yù)測時間長,無概率性輸出等問題,僅適用于小樣本問題的建模。水質(zhì)模擬模型雖然應(yīng)用范圍廣、適合中長期預(yù)測,但側(cè)重于水環(huán)境系統(tǒng)(劉曉等,2018),需要考慮的參數(shù)較多,模擬的結(jié)果缺乏交互性,F(xiàn)an et al.(2012)通過結(jié)合水文工程中心河流分析系統(tǒng)(HEC-RAS)和S-P(Streeter-Phelps)模型模擬潮汐河水質(zhì),然而受模型參數(shù)的影響,對生物化學(xué)需氧量和溶解氧的模擬結(jié)果預(yù)估過高。混沌理論為河流水質(zhì)變化提供了新的預(yù)測方法,但基于混沌理論建立的預(yù)測模式較為復(fù)雜,需要提供大量時序數(shù)據(jù)資料,適應(yīng)于短期的水質(zhì)預(yù)測(Zaldivar et al.,2000;杜建等,2012)。神經(jīng)網(wǎng)絡(luò)預(yù)測模型具備良好的記憶功能以及精度高、適應(yīng)性強(qiáng)、運算速度快等優(yōu)點(唐亦舜等,2022),在具有不確定性因素和非線性特征的各類水質(zhì)指標(biāo)的水質(zhì)預(yù)測方面(Maier et al.,2010)有著廣闊的應(yīng)用前景。鑒于傳統(tǒng)水質(zhì)預(yù)測方法對水質(zhì)監(jiān)測數(shù)據(jù)的依賴性較強(qiáng),預(yù)測范圍受限以及缺乏考慮參數(shù)的影響,預(yù)測精度較低等問題,神經(jīng)網(wǎng)絡(luò)模型為水質(zhì)預(yù)測研究提供了科學(xué)有效的途徑。

    目前,對綿陽市主要河流的研究主要集中在水環(huán)境容量計算、水污染總量控制及單一的水質(zhì)評價,缺乏水質(zhì)綜合評價及預(yù)測相關(guān)的研究(沈淞濤,2005;戴本林等,2008)。因此,本研究基于綿陽市水質(zhì)監(jiān)測斷面的污染物指標(biāo)監(jiān)測數(shù)據(jù),結(jié)合模糊綜合評價法與主成分分析法進(jìn)行水質(zhì)綜合評價,以獲取水質(zhì)級別并識別水質(zhì)較差的年份和監(jiān)測斷面,并基于水質(zhì)評價結(jié)果構(gòu)建NARX 神經(jīng)網(wǎng)絡(luò)模型對未來區(qū)域水環(huán)境的水質(zhì)狀況進(jìn)行合理的預(yù)判,實現(xiàn)科學(xué)預(yù)防水質(zhì)污染,以期為區(qū)域水環(huán)境的評價與管理、預(yù)測與規(guī)劃提供參考。

    1 研究區(qū)概況

    綿陽市,被稱為中國科技城,是四川第二大經(jīng)濟(jì)體和成渝城市群區(qū)域中心城市,位于四川盆地西北部,介于30°42′-33°03′N、103°45′-105°43′E 之間,總面積約為20 248.4 km2,地形地貌類型復(fù)雜多樣,其中山地區(qū)占61.0%,丘陵區(qū)占20.4%,平原區(qū)占18.6%,地勢起伏大,呈西北部高,東南部低,屬多山多丘陵的山地城市(圖1)。全市受地貌影響,降水豐沛,徑流量大,江河縱橫,水系發(fā)達(dá),境內(nèi)有大小河流及溪溝3 000 余條,均分別注入嘉陵江支流涪江、白龍江與西河,全屬嘉陵江水系。根據(jù)2022年綿陽市環(huán)境質(zhì)量狀況年報,全市境內(nèi)地表水水質(zhì)總體較好,均能達(dá)到或優(yōu)于劃定的水環(huán)境功能類別要求,但芙蓉溪及魯班水庫水質(zhì)為良,呈中營養(yǎng)狀態(tài),通口河水質(zhì)有所下降。綿陽市境內(nèi)的主要河流,不僅直接關(guān)系到全市整體水環(huán)境安全與質(zhì)量,還會影響綿陽市社會經(jīng)濟(jì)可持續(xù)發(fā)展戰(zhàn)略的實施,因此,對綿陽市主要河流的評價及預(yù)測是保護(hù)及管理河流水環(huán)境、推動經(jīng)濟(jì)高質(zhì)量綠色發(fā)展的關(guān)鍵。

    圖1 水質(zhì)監(jiān)測站點分布圖Figure 1 Distribution of water quality monitoring stations

    2 研究方法

    2.1 數(shù)據(jù)來源與處理

    水質(zhì)數(shù)據(jù)來源于綿陽市生態(tài)環(huán)境監(jiān)測站、綿陽市水環(huán)境形式分析報告以及長江生態(tài)環(huán)境保護(hù)修復(fù)城市駐點組長期監(jiān)測。根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)(國家環(huán)境保護(hù)總局,2002),結(jié)合綿陽市水質(zhì)實際污染特征、來源以及歷史監(jiān)測結(jié)果,本研究按照生態(tài)環(huán)境行業(yè)現(xiàn)行標(biāo)準(zhǔn)及技術(shù)規(guī)范對水質(zhì)樣品進(jìn)行采集并分析,篩選出了溶解氧(DO)、高錳酸鹽指數(shù)(CODMn)、五日生化需氧量(BOD5)和氨氮(NH3-N)這4 類對水質(zhì)影響較大并能反映水質(zhì)狀況的指標(biāo)作為主要評價因子,對綿陽市12 個地表水環(huán)境質(zhì)量國控、省控以及市控監(jiān)測斷面2014-2022 年枯水期(1-4 月和12 月)、平水期(5 月和11 月)、豐水期(6-10 月)、年均值水質(zhì)監(jiān)測結(jié)果進(jìn)行評價分析(表1)。

    表1 綿陽市主要監(jiān)測斷面及其水環(huán)境功能標(biāo)準(zhǔn)Table 1 Main monitoring sections and water environment function standards in Mianyang City

    2.2 模糊綜合評價法

    本研究首先根據(jù)篩選出的4 種主要水質(zhì)評價因子,建立評價因子集={DO,CODMn,BOD5,NH3-N},確定評價標(biāo)準(zhǔn)集為={Ⅰ,Ⅱ,Ⅲ,Ⅳ,Ⅴ}。由于水質(zhì)污染程度和水質(zhì)分級標(biāo)準(zhǔn)具有模糊性,故用隸屬度來刻畫分級界限隸屬度(李婧等,2020),通過隸屬函數(shù),按照水質(zhì)分級標(biāo)準(zhǔn),計算出各項因子分別對各級水質(zhì)指標(biāo)的隸屬度,從而確定模糊關(guān)系矩陣。然后,采用超標(biāo)倍數(shù)法對各水質(zhì)指標(biāo)賦予權(quán)重(張婉萍等,2022)。最后,根據(jù)模糊關(guān)系矩陣和評價因子權(quán)重集得到各監(jiān)測斷面水質(zhì)等級(王守坤等,2013)。

    2.3 主成分分析法

    主成分分析法是一種利用降維的思想對多個變量以線性變換形式篩選出重要變量的多元統(tǒng)計分析方法(王麗婧等,2013)。本研究以年份為樣本對選取的監(jiān)測斷面各水期水質(zhì)進(jìn)行主成分分析,分析判斷各污染指標(biāo)之間的相關(guān)性,并根據(jù)各監(jiān)測斷面在2014-2022 年及不同水期的污染綜合得分進(jìn)行排序,以評價水質(zhì)優(yōu)劣程度,即以主成分綜合得分大小反映水質(zhì)的污染程度,即得分越高水質(zhì)越差。

    2.4 NARX 神經(jīng)網(wǎng)絡(luò)模型

    NARX(nonlinear auto regressive model with exogenous inputs)神經(jīng)網(wǎng)絡(luò)模型是通過前期長時間序列的隨機(jī)變量的線性組合來描述未來某一時刻隨機(jī)變量的非線性自回歸模型(Fabio et al.,2022)。NARX 神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)由輸入層、隱含層、輸出層以及輸出到輸入的延時構(gòu)成(圖2)。本研究在2014-2022 年綿陽市水質(zhì)監(jiān)測數(shù)據(jù)及水質(zhì)評價結(jié)果的基礎(chǔ)上,采用MATLAB 2020b 工具箱的NARX神經(jīng)網(wǎng)絡(luò)模型預(yù)測2023-2030 年各監(jiān)測斷面的水質(zhì)指標(biāo)濃度,從而客觀反映未來水質(zhì)變化趨勢。

    圖2 NARX 神經(jīng)網(wǎng)絡(luò)預(yù)測模型結(jié)構(gòu)Figure 2 Structure of NARX neural network prediction model

    3 結(jié)果與討論

    3.1 水質(zhì)綜合評價結(jié)果及分析

    表2 為綿陽市境內(nèi)主要河流2014-2022 年水質(zhì)監(jiān)測斷面各水期4 項水質(zhì)指標(biāo)平均監(jiān)測數(shù)據(jù)的描述性統(tǒng)計結(jié)果。由表可知,從監(jiān)測斷面的水質(zhì)指標(biāo)平均濃度來看,除DO 外,其余3 項水質(zhì)指標(biāo)平均濃度值最高處均為芙蓉溪仙魚橋斷面,說明該斷面水質(zhì)污染較為嚴(yán)重,而平武水文站總體上水質(zhì)則較好。從不同水期各水質(zhì)指標(biāo)的平均濃度可知,BOD5和NH3-N 的平均濃度值在枯水期較高,這與入河污染物受徑流量影響顯著有關(guān)(何子建等,2022),在枯水期內(nèi)BOD5和NH3-N 濃度隨徑流量的減少而增大。而CODMn的平均濃度值在豐水期較高,這與豐水期雨季或氣溫持續(xù)升高,河中水量豐富密切相關(guān),并且徑流量增大引起的CODMn污染物入河的增加量大于徑流量對河流CODMn污染物的稀釋量。從各水期水質(zhì)監(jiān)測斷面的變異系數(shù)來看,NH3-N的變異系數(shù)在各水期之間均較高,說明其濃度值受水期影響波動幅度也較大,而DO 的變異系數(shù)在水期之間均較低,說明其波動變化范圍較為穩(wěn)定。

    表2 不同水期監(jiān)測斷面平均水質(zhì)指標(biāo)濃度的描述性統(tǒng)計Table 2 Descriptive statistics of average water quality index concentration of monitoring sections in different water periods

    為進(jìn)一步明確監(jiān)測斷面各水期的主要污染物,本研究根據(jù)2014-2022 年所選監(jiān)測斷面水質(zhì)指標(biāo)在不同水期內(nèi)的平均值,計算污染因子權(quán)重并確定各水期內(nèi)的主要污染物(表3)。由表可知,研究范圍內(nèi)枯水期的主要污染物為NH3-N,BOD5。同理能夠判斷出平水期和豐水期的主要污染物為BOD5和CODMn。NH3-N 在枯水期的權(quán)重高于平水期和豐水期,而NH3-N 作為城鎮(zhèn)生活污水的主要來源,說明枯水期水質(zhì)受生活污水的影響較大;BOD5在豐水期的權(quán)重高于枯水期與平水期,這也與不同水期下河流斷面徑流和含沙量的差異有關(guān)(伏介雄等,2008),豐水期降雨量大,河流水量增大且含沙量高,從而導(dǎo)致河流有機(jī)污染物濃度的升高。

    表3 2014-2022 年監(jiān)測斷面各水期污染指標(biāo)權(quán)重歸一化結(jié)果Table 3 Pollution index weight normalization results of water periods in the monitoring section from 2014 to 2022

    根據(jù)各斷面水質(zhì)指標(biāo)的描述性統(tǒng)計結(jié)果,選取水質(zhì)較好的平武水文站斷面和水質(zhì)較差的芙蓉溪仙魚橋斷面進(jìn)行分析(圖3)。平武水文站為涪江上游基本控制站,而涪江上游流域地廣人稀、工礦企業(yè)較少,斷面內(nèi)的污染主要來自農(nóng)田、城市地表徑流以及沿岸的畜禽養(yǎng)殖。芙蓉溪仙魚橋斷面位于綿陽市經(jīng)濟(jì)較為發(fā)達(dá)、人口較為密集的游仙區(qū),區(qū)內(nèi)工廠企業(yè)數(shù)量僅次于涪城區(qū),工業(yè)廢水排放總量較大并且該斷面附近支流多,水系密集,是導(dǎo)致水質(zhì)相對較差的重要原因,且NH3-N 為主要超標(biāo)項目,主要來源于以城鎮(zhèn)生活污水及工業(yè)廢水為主的點源污染(胡開明等,2015)。對此綿陽市加強(qiáng)了沿岸主要污染源排查及治理,對沿河的大污染工廠企業(yè)實行關(guān)、停、并、轉(zhuǎn),并采取了嚴(yán)禁規(guī)模養(yǎng)殖戶將廢污水直接排入河道等措施。

    圖3 2014-2022 年各監(jiān)測斷面水質(zhì)指標(biāo)變化對比圖Figure 3 Comparison of water quality indexes of different monitoring sections from 2014 to 2022

    在水質(zhì)評價方面,采用模糊綜合評價法對2014-2022 年12 個監(jiān)測斷面進(jìn)行水質(zhì)評價及分析,以此反映綜合水質(zhì)污染狀況。對于水質(zhì)指標(biāo)的年平均值而言(表4),梓潼垢家渡斷面2016 年的水質(zhì)級別為Ⅲ級;芙蓉溪仙魚橋斷面在2015 年、2016 年水質(zhì)級別分別為Ⅲ級、Ⅱ級,這一水質(zhì)級別評價結(jié)果也與2016 年綿陽市環(huán)境質(zhì)量狀況年報中芙蓉溪仙魚橋斷面水質(zhì)有所好轉(zhuǎn)相符;總體上2014-2022年水質(zhì)級別達(dá)標(biāo)斷面占比為100%,梓潼垢家渡和芙蓉溪仙魚橋斷面年均水質(zhì)水平總體上與其他斷面相比較低,但兩個斷面均達(dá)到了水環(huán)境功能類別劃分(Ⅲ級)的要求,水質(zhì)級別較差的斷面主要出現(xiàn)在2014-2018 年之間,而2019-2022 年水質(zhì)較好,整體水質(zhì)呈現(xiàn)向好趨勢,這與2018 年7 月起施行的《綿陽市水污染防治條例》密切相關(guān)。在枯水期,2015 年的芙蓉溪仙魚橋斷面與2016 年的梓潼垢家渡斷面的水質(zhì)級別均為Ⅳ級,超過了其水環(huán)境功能劃定類別Ⅲ級,凱江老南橋在2016 年水質(zhì)級別波動較大,從Ⅰ級降至Ⅲ級,水質(zhì)為輕度污染,其余斷面水質(zhì)情況均較好,總體水質(zhì)級別達(dá)標(biāo)斷面占比92%(表5)。在平水期,僅芙蓉溪仙魚橋斷面在2017 年水質(zhì)級別為Ⅲ級,水質(zhì)較差,其余斷面水質(zhì)情況較好,水質(zhì)級別達(dá)標(biāo)斷面占比為100%(表6)。在豐水期,水質(zhì)級別較低的斷面出現(xiàn)在2014-2016 年之間,分別為芙蓉溪仙魚橋、梓潼垢家渡、凱江老南橋與豐谷,其余斷面的水質(zhì)級別在2014-2022 年均為Ⅰ級,總體水質(zhì)級別達(dá)標(biāo)率為100%(表7)。

    表4 2014-2022 年各監(jiān)測斷面水質(zhì)級別Table 4 Water quality levels of each monitored section from 2014 to 2022

    表5 2014-2022 年枯水期各監(jiān)測斷面水質(zhì)級別Table 5 Water quality levels of each monitored section during the low water period from 2014 to 2022

    表6 2014-2022 年平水期各監(jiān)測斷面水質(zhì)級別Table 6 Water quality levels of each monitored section during the normal water period from 2014 to 2022

    表7 2014-2022 年豐水期各監(jiān)測斷面水質(zhì)級別Table 7 Water quality levels of each monitored section during the high water period from 2014 to 2022

    3.2 水質(zhì)污染程度變化結(jié)果及分析

    為探究水質(zhì)指標(biāo)之間的相關(guān)性對于水質(zhì)污染的影響,本研究首先采用主成分分析法對研究范圍內(nèi)水質(zhì)指標(biāo)的年平均值進(jìn)行Pearson 相關(guān)性分析,Pearson 相關(guān)性是分析不同指標(biāo)之間同源性的常用方法,相關(guān)性高的指標(biāo)之間具有相似的污染源或遷移特征(王磊等,2020)。由各水質(zhì)指標(biāo)的相關(guān)性系數(shù)可知(表8),CODMn與BOD5在3 個水期內(nèi)均呈極顯著正相關(guān)關(guān)系,表明兩個水質(zhì)指標(biāo)存在同源性,推測水體中的CODMn與BOD5主要來源于工業(yè)及生活污染源,并且在理論上BOD5還與水體中的COD 具備一定的比例關(guān)系,可通過判斷其比例關(guān)系來識別水體污染類型,若COD 與BOD5的比值越大,則還原性無機(jī)物為水體主要污染物,若比值越小,則水體遭受到的有機(jī)物污染較為嚴(yán)重,同時CODMn與BOD5測定值的范圍還可由COD 進(jìn)行推測(張野等,2012),利用水質(zhì)指標(biāo)之間的關(guān)聯(lián)性去預(yù)測或估算其他指標(biāo),從而實現(xiàn)對監(jiān)測數(shù)據(jù)的預(yù)判和數(shù)據(jù)合理性的審核并確保水質(zhì)數(shù)據(jù)監(jiān)測的準(zhǔn)確有效;DO 與NH3-N 在枯水期、平水期為顯著負(fù)相關(guān),與CODMn在平水期、豐水期分別為顯著與極顯著負(fù)相關(guān),DO 與NH3-N、CODMn呈負(fù)相關(guān)這一結(jié)果與宋曉霄等(2012)對水質(zhì)指標(biāo)進(jìn)行數(shù)據(jù)統(tǒng)計和相關(guān)性分析后的結(jié)果相符,CODMn與NH3-N同屬耗氧有機(jī)物,當(dāng)水質(zhì)DO 濃度降低時,有機(jī)物降解速率減緩,CODMn與NH3-N 濃度反而呈升高趨勢,通過分析DO 與其相關(guān)性,可合理判別好氧有機(jī)物的變化趨勢并識別水質(zhì)主要污染物類型。

    表8 2014-2022 年不同水期各水質(zhì)指標(biāo)的相關(guān)性系數(shù)Table 8 Correlation coefficients of water quality indexes in different water periods from 2014 to 2022

    然后,根據(jù)所選監(jiān)測斷面除DO 以外的水質(zhì)數(shù)據(jù)分水期計算綜合主成分值,從而評價各監(jiān)測斷面的水質(zhì)污染程度(圖4),并對水質(zhì)較好與較差的斷面在2014-2022 年之間的主成分得分進(jìn)行綜合評價比較(表9)。結(jié)果表明,在3 個水期之間,芙蓉溪仙魚橋斷面水質(zhì)污染均較差,而平武水文站斷面水質(zhì)均較好,凱江老南橋、梓潼垢家渡及天仙鎮(zhèn)大佛寺斷面受水期影響較大且水質(zhì)得分排序情況在斷面之間不穩(wěn)定,其余斷面在各水期間排序結(jié)果無變化,受水期影響較小。與其他斷面相比,芙蓉溪仙魚橋、平武水文站斷面分別在2015、2016 年水質(zhì)污染較重,但2020-2022 年兩個斷面的整體水質(zhì)狀況有所改善,而這也得益于2020 年四川省印發(fā)了《2020 年全省河湖管理保護(hù)工作要點》以及各地方加強(qiáng)了河道采砂綜合整治以及河湖管理基礎(chǔ)工作的力度等一系列舉措。

    圖4 2014-2022 年不同水期各監(jiān)測斷面水質(zhì)平均主成分得分Figure 4 Average principal component scores of water quality in different sections during different water periods from 2014 to 2022

    表9 2014-2022 年監(jiān)測斷面水質(zhì)主成分得分及排序結(jié)果Table 9 Principal component scores and ranking results of water quality in monitored sections from 2014 to 2022

    3.3 水質(zhì)變化趨勢預(yù)測及分析

    在水質(zhì)預(yù)測方面,本研究選用NARX 神經(jīng)網(wǎng)絡(luò)模型并對2014-2022 年12 個斷面水質(zhì)數(shù)據(jù)分季度進(jìn)行歸一化與反歸一化,為了驗證模型的性能,根據(jù)算法準(zhǔn)則(宋新山,2008;Ayodele et al.,2021;包志家等,2021),將數(shù)據(jù)集隨機(jī)劃分為36 個目標(biāo)時間步,用于訓(xùn)練、驗證和測試。指定隱藏層中有10 個神經(jīng)元,時延為2、4,采用Levenberg-Marquardt算法(Salami et al.,2015)對網(wǎng)絡(luò)進(jìn)行訓(xùn)練,選用均方誤差(MSE)和相關(guān)系數(shù)(r)值指標(biāo)對訓(xùn)練后的模型進(jìn)行定量評價(表10)。模型預(yù)測結(jié)果表明,4 種水質(zhì)指標(biāo)的r均大于0.7,MSE 均小于0.05,說明模型對水質(zhì)的預(yù)測精度良好,具有較強(qiáng)的可靠性和預(yù)測性能。為了使NARX 模型具有較高實用性,預(yù)測對象選取水質(zhì)評價結(jié)果中具代表性的芙蓉溪仙魚橋斷面為例進(jìn)行預(yù)測結(jié)果的進(jìn)一步分析。DO、CODMn、BOD5以及NH3-N 的平均相對誤差分別為5.914%、-0.190%、21.491%以及25.366%(圖5),總體上的平均相對誤差均在可接受范圍內(nèi)。由模型對芙蓉溪仙魚橋斷面2014-2022 年的NH3-N水質(zhì)指標(biāo)的預(yù)測值結(jié)果與實測值進(jìn)行對比可知(圖6),芙蓉溪仙魚橋斷面NH3-N 的相關(guān)系數(shù)在訓(xùn)練、驗證、測試以及整體上分別為0.963、0.415、0.561以及0.868,模型預(yù)測結(jié)果與實測水質(zhì)數(shù)據(jù)的相關(guān)系數(shù)整體上達(dá)到了0.8 以上。通過NARX 神經(jīng)網(wǎng)絡(luò)模型模擬結(jié)果的誤差及相關(guān)系數(shù)分析,模型的預(yù)測性能良好,再次證實了該模型可用于進(jìn)行水質(zhì)預(yù)測,故基于此預(yù)測了芙蓉溪仙魚橋斷面2023-2030年水質(zhì)指標(biāo)的變化情況(圖7)。結(jié)果表明,2023-2030 年DO 與CODMn平均濃度均高于2014-2022 年對應(yīng)指標(biāo)的平均濃度值,而BOD5與NH3-N平均濃度均低于2014-2022 年對應(yīng)指標(biāo)的平均濃度值,說明水體總體改善,但受到有機(jī)物污染的程度有所提高。預(yù)測的水質(zhì)指標(biāo)中,DO 濃度值年際變化不大,較為穩(wěn)定,CODMn濃度值在2024 年及2027 年存在小幅度的波谷,整體上呈波浪形變化,BOD5濃度值呈先上升后下降再回升的變化趨勢,而NH3-N 濃度值整體上呈現(xiàn)逐年下降趨勢,表明NH3-N 對未來水質(zhì)污染的貢獻(xiàn)度有所降低。

    表10 NARX 神經(jīng)網(wǎng)絡(luò)模型預(yù)測的相關(guān)系數(shù)與均方誤差Table 10 Correlation coefficient and mean square error predicted by NARX neural network model

    圖5 2014-2022 年芙蓉溪仙魚橋斷面水質(zhì)指標(biāo)實測值與模擬值結(jié)果對比圖Figure 5 Comparison of measured and simulated water quality indexes in Xianyuqiao section of Furong stream from 2014 to 2022

    圖6 芙蓉溪仙魚橋斷面NH3-N 水質(zhì)指標(biāo)的相關(guān)性回歸圖Figure 6 Correlation regression diagram of NH3-N water quality index in Xianyuqiao section of Furong stream

    圖7 2023-2030 年芙蓉溪仙魚橋斷面水質(zhì)指標(biāo)變化情況Figure 7 Changes of water quality indexes in Xianyuqiao section of Furong stream from 2023 to 2030

    為驗證模型預(yù)測的準(zhǔn)確性,基于水質(zhì)預(yù)測結(jié)果對所選斷面的水質(zhì)再次進(jìn)行水質(zhì)評價。首先,對于2025 及2030 年各監(jiān)測斷面的水質(zhì)指標(biāo)變化情況而言(圖8),豐谷及魯班島斷面2030 年的4 項水質(zhì)指標(biāo)的年均濃度值均高于2025 年,百頃及北川通口斷面2030 年除NH3-N 年均濃度值低于2025 年外,其余3 項水質(zhì)指標(biāo)年均濃度值高于2025 年;梓潼垢家渡及飲馬橋斷面2030 年的DO 及NH3-N年均濃度值高于2025 年,CODMn及BOD5年均濃度值低于2025 年。然后,由所選斷面在2025 及2030年的水質(zhì)指標(biāo)預(yù)測值的主成分得分及排序結(jié)果可知(表11),芙蓉溪仙魚橋斷面的水質(zhì)在所選斷面中仍較差,與其他斷面相比,平武水文站、福田壩斷面的水質(zhì)分別在2025 年、2030 年較好,預(yù)測結(jié)果與歷史時期斷面總體水質(zhì)狀況相符。最后,對所選的12 個斷面的2023-2030 年的水質(zhì)指標(biāo)預(yù)測值進(jìn)行模糊綜合評價,得出各斷面在2023-2030 年的水質(zhì)級別(表12)。結(jié)果表明,天仙鎮(zhèn)大佛寺斷面在2028-2030 年的水質(zhì)級別由2023-2027 年的Ⅰ級降為Ⅱ級,水質(zhì)有所下降;芙蓉溪仙魚橋斷面在2025及2030 年的水質(zhì)級別分別為Ⅱ級與Ⅲ級,其余斷面在未來八年總體水質(zhì)級別均較好,為Ⅰ級,水質(zhì)預(yù)測結(jié)果客觀反映了綿陽市未來的水質(zhì)變化趨勢,有利于識別及掌控區(qū)域水質(zhì)污染程度,可為河流污染預(yù)警及水環(huán)境管理規(guī)劃工作提供參考價值。

    表11 2025 及2030 年各監(jiān)測斷面水質(zhì)平均主成分得分及排序結(jié)果Table 11 Average principal component scores and ranking results of water quality in each monitored section in 2025 and 2030

    表12 2023-2030 年監(jiān)測斷面水質(zhì)評價結(jié)果Table 12 Water quality evaluation results of monitored sections during 2023-2030

    3.4 水環(huán)境管理及規(guī)劃啟示

    分析水體是否受到污染、污染程度、污染原因以及發(fā)展趨勢,對于水污染的控制及預(yù)防、合理規(guī)劃水資源的開發(fā)利用至關(guān)重要。本研究通過對綿陽市2014-2022 年主要水質(zhì)監(jiān)測斷面的水質(zhì)及其指標(biāo)分水期進(jìn)行綜合評價與分析,從水質(zhì)指標(biāo)的權(quán)重分析了各水期的主要污染物,從水質(zhì)指標(biāo)的相關(guān)性分析了各指標(biāo)對水質(zhì)狀況的影響,計算了各斷面在各水期內(nèi)的水質(zhì)級別及得分,從而較為全面客觀地掌握歷史水質(zhì)狀況。由評價結(jié)果可知,從時間上來看,2014-2018 年水質(zhì)相對較差,2019-2022 年水質(zhì)有所好轉(zhuǎn),枯水期的水質(zhì)斷面達(dá)標(biāo)率低于平水期和豐水期,應(yīng)加強(qiáng)此階段未達(dá)標(biāo)水質(zhì)監(jiān)測斷面周圍水域的監(jiān)管及治理;從空間上來看,平武水文站及北川通口斷面水質(zhì)整體較好,均位于北川縣,污染較嚴(yán)重的區(qū)域集中在綿陽市中部水系交匯處(游仙區(qū)的芙蓉溪仙魚橋斷面)及南部(三臺縣的凱江老南橋斷面、梓潼縣的垢家渡斷面以及鹽亭縣的天仙鎮(zhèn)大佛寺斷面)。因此,芙蓉溪、凱江以及梓江內(nèi)這些區(qū)域可以作為未來重點治理對象,并對其附近工廠及沿岸企業(yè)污水排放情況實行嚴(yán)格監(jiān)查。

    在此基礎(chǔ)上,本研究基于歷史實測水質(zhì)數(shù)據(jù)采用了NARX 神經(jīng)網(wǎng)絡(luò)模型對河流水質(zhì)監(jiān)測斷面各項水質(zhì)指標(biāo)進(jìn)行2023-2030 年的水質(zhì)預(yù)測,在滿足水質(zhì)預(yù)測精度要求的前提下,同樣對水質(zhì)監(jiān)測斷面各水質(zhì)指標(biāo)未來的水質(zhì)狀況變化、主成分得分以及水質(zhì)級別進(jìn)行了計算分析,以了解未來水質(zhì)狀況的發(fā)展趨勢及規(guī)律。由水質(zhì)預(yù)測結(jié)果可知,鹽亭縣的天仙鎮(zhèn)大佛寺斷面在2028 年、游仙區(qū)的芙蓉溪仙魚橋斷面在2025 及2030 年均有水質(zhì)下降趨勢。總之,本研究提出的基于模糊綜合評價法-主成分分析法-NARX 神經(jīng)網(wǎng)絡(luò)模型研究框架,能夠在一定程度上客觀真實地反映區(qū)域水環(huán)境的水質(zhì)歷史狀況及未來趨勢。

    然而,鑒于本研究受一些客觀因素的限制,如水質(zhì)數(shù)據(jù)、水文水利參數(shù)、社會經(jīng)濟(jì)統(tǒng)計數(shù)據(jù)等的可獲取性,在水質(zhì)評價方面,目前僅對該區(qū)域主要水質(zhì)監(jiān)測斷面,分別采用模糊綜合評價法和主成分分析法進(jìn)行水質(zhì)等級評定和水質(zhì)優(yōu)劣排序,在未來的研究中,應(yīng)該對水質(zhì)波動較大和水質(zhì)較差的區(qū)域的水質(zhì)監(jiān)測站點進(jìn)行加密,實行更加精細(xì)密集地布控,識別出水質(zhì)污染重點區(qū)域及擴(kuò)散區(qū)域,進(jìn)行空間上由點到面的研究和具有針對性的管控,從而得到區(qū)域整體水環(huán)境更加全面準(zhǔn)確的水質(zhì)評價結(jié)果。在水質(zhì)預(yù)測方面,NARX 神經(jīng)網(wǎng)絡(luò)模型使用單一的水質(zhì)指標(biāo)點數(shù)據(jù)作為輸入預(yù)測了水質(zhì)整個條件分布,但目前尚未清楚這些預(yù)測如何及時向前傳播,且存在單步和單因素預(yù)測的局限性,以及對網(wǎng)絡(luò)時延、隱藏神經(jīng)元數(shù)量和激活函數(shù)的考慮仍有優(yōu)化的余地。因此,后續(xù)研究可將水質(zhì)及其影響水質(zhì)變化的驅(qū)動因素作為神經(jīng)網(wǎng)絡(luò)模型的多輸入指標(biāo),以及考慮將NARX 人工神經(jīng)網(wǎng)絡(luò)與混合動態(tài)神經(jīng)網(wǎng)絡(luò)、遺傳算法和其他機(jī)器學(xué)習(xí)模型進(jìn)行聯(lián)合調(diào)用,進(jìn)一步構(gòu)建更加完整科學(xué)的水質(zhì)預(yù)測模型,提高水質(zhì)指標(biāo)模擬精度,賦予水環(huán)境指標(biāo)物理涵義。

    綜上所述,本研究結(jié)果對于區(qū)域水環(huán)境管理及規(guī)劃工作的啟示如下:一是針對易受水期影響,枯水期水質(zhì)不達(dá)標(biāo)的區(qū)域,應(yīng)通過優(yōu)化超標(biāo)區(qū)域內(nèi)外的水資源配置,提升河流連通性和水環(huán)境容量,從而改善枯水期水質(zhì)狀況;二是針對易受工業(yè)、農(nóng)業(yè)污染等影響的區(qū)域,應(yīng)調(diào)整產(chǎn)業(yè)結(jié)構(gòu),加快產(chǎn)業(yè)轉(zhuǎn)型升級,加大對污染源的控制,從源頭上減少污染物入河量;三是針對水質(zhì)變化復(fù)雜的區(qū)域,應(yīng)建立全過程的水環(huán)境監(jiān)控預(yù)警系統(tǒng)及長效管理機(jī)制,提高污染防治水平。

    4 結(jié)論

    通過對綿陽市主要河流及湖庫水質(zhì)監(jiān)測斷面的水質(zhì)監(jiān)測數(shù)據(jù)進(jìn)行整理分析,利用模糊綜合評價方法及主成分分析法進(jìn)行水質(zhì)評價,并根據(jù)評價結(jié)果,建立NARX 神經(jīng)網(wǎng)絡(luò)模型對未來水質(zhì)變化狀況進(jìn)行預(yù)測,結(jié)論如下:

    1)通過對監(jiān)測斷面各水期污染指標(biāo)平均值進(jìn)行權(quán)重歸一化及模糊綜合評價可知,NH3-N 和BOD5為枯水期主要污染物,平水期和豐水期的主要污染物為BOD5和CODMn;除枯水期外,水質(zhì)監(jiān)測斷面在平水期、豐水期以及年平均狀態(tài)下的達(dá)標(biāo)率均為100%,說明枯水期河流水質(zhì)更易受到污染并超標(biāo)。

    2)根據(jù)主成分分析可知,CODMn與BOD5在3個水期內(nèi)均呈極顯著正相關(guān)關(guān)系,DO 與NH3-N 之間在枯水期、平水期為顯著負(fù)相關(guān),在豐水期則無顯著關(guān)系,DO 與CODMn在平水期、豐水期分別為顯著與極顯著負(fù)相關(guān),應(yīng)基于不同水期各污染指標(biāo)的相關(guān)性對河流水質(zhì)污染進(jìn)行科學(xué)高效地管控。

    3)由NARX 神經(jīng)網(wǎng)絡(luò)模型預(yù)測結(jié)果可知,DO、CODMn、BOD5以及NH3-N 濃度的均方誤差均小于0.05,相關(guān)系數(shù)均大于0.7,構(gòu)建的預(yù)測模型能夠客觀反映未來水質(zhì)變化趨勢。

    因此,本研究提出的基于模糊綜合評價法-主成分分析法-NARX 神經(jīng)網(wǎng)絡(luò)模型研究框架,水質(zhì)綜合評價結(jié)果客觀準(zhǔn)確,可為水環(huán)境功能區(qū)劃及水質(zhì)目標(biāo)的確定提供依據(jù),根據(jù)水質(zhì)評價結(jié)果并結(jié)合歷史水質(zhì)數(shù)據(jù)采用神經(jīng)網(wǎng)絡(luò)模型進(jìn)行水質(zhì)預(yù)測,預(yù)測結(jié)果科學(xué)有效,有利于類似河流的水污染控制規(guī)劃及水環(huán)境綜合管理工作的開展。

    猜你喜歡
    斷面水質(zhì)污染
    水質(zhì)抽檢豈容造假
    環(huán)境(2023年5期)2023-06-30 01:20:01
    ??? ??? ?? ????? ? ??(斷面)
    一月冬棚養(yǎng)蝦常見水質(zhì)渾濁,要如何解決?這9大原因及處理方法你要知曉
    堅決打好污染防治攻堅戰(zhàn)
    堅決打好污染防治攻堅戰(zhàn)
    基于Dijkstra算法的最優(yōu)解列斷面快速搜索方法
    對抗塵污染,遠(yuǎn)離“霾”伏
    都市麗人(2015年5期)2015-03-20 13:33:49
    國內(nèi)首臺全斷面煤巷高效掘進(jìn)機(jī)誕生并交付
    天鋼Ф210mm斷面圓坯的開發(fā)與生產(chǎn)實踐
    天津冶金(2014年4期)2014-02-28 16:52:32
    水質(zhì)總磷測定存在的問題初探
    河南科技(2014年23期)2014-02-27 14:19:07
    大陆偷拍与自拍| 国产成人精品久久二区二区免费| 母亲3免费完整高清在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲成人免费av在线播放| 亚洲少妇的诱惑av| av线在线观看网站| 曰老女人黄片| 在线 av 中文字幕| 亚洲欧美清纯卡通| av电影中文网址| 夜夜骑夜夜射夜夜干| 精品国产国语对白av| 日日摸夜夜添夜夜爱| www.精华液| 多毛熟女@视频| av又黄又爽大尺度在线免费看| 久久这里只有精品19| 国产一区二区在线观看av| 日本午夜av视频| 电影成人av| 大话2 男鬼变身卡| 中文字幕人妻丝袜制服| 免费女性裸体啪啪无遮挡网站| 欧美人与善性xxx| 成年人免费黄色播放视频| 国产成人一区二区在线| 男女床上黄色一级片免费看| 91麻豆精品激情在线观看国产 | 老熟女久久久| 秋霞在线观看毛片| 亚洲欧美清纯卡通| 老司机亚洲免费影院| 久久精品亚洲熟妇少妇任你| 欧美人与性动交α欧美软件| 一本大道久久a久久精品| 2021少妇久久久久久久久久久| netflix在线观看网站| 最近手机中文字幕大全| 欧美人与性动交α欧美软件| 黄色视频在线播放观看不卡| 在线观看www视频免费| 欧美日韩亚洲高清精品| 999精品在线视频| 欧美在线黄色| 青春草视频在线免费观看| 波多野结衣av一区二区av| 国产主播在线观看一区二区 | 亚洲伊人久久精品综合| 精品一区二区三区四区五区乱码 | 91国产中文字幕| 欧美精品一区二区免费开放| 欧美另类一区| 欧美日韩成人在线一区二区| 在线观看免费午夜福利视频| 国产麻豆69| 爱豆传媒免费全集在线观看| 一级片'在线观看视频| 免费av中文字幕在线| 青青草视频在线视频观看| 欧美日韩视频高清一区二区三区二| 人人澡人人妻人| 美女福利国产在线| 高清欧美精品videossex| 亚洲国产精品成人久久小说| 成年人午夜在线观看视频| 秋霞在线观看毛片| 国产无遮挡羞羞视频在线观看| 悠悠久久av| 精品一区二区三卡| 国产主播在线观看一区二区 | 色综合欧美亚洲国产小说| 国产视频首页在线观看| 久久ye,这里只有精品| 黑人欧美特级aaaaaa片| 99久久综合免费| 亚洲午夜精品一区,二区,三区| 日韩 亚洲 欧美在线| 亚洲精品中文字幕在线视频| 国产亚洲av片在线观看秒播厂| 欧美97在线视频| 18禁国产床啪视频网站| 亚洲精品美女久久av网站| 亚洲中文av在线| 久久人妻熟女aⅴ| 汤姆久久久久久久影院中文字幕| 18在线观看网站| 国产国语露脸激情在线看| 熟女少妇亚洲综合色aaa.| 欧美日韩视频高清一区二区三区二| 久久精品久久精品一区二区三区| 午夜91福利影院| av不卡在线播放| 热99国产精品久久久久久7| 精品少妇一区二区三区视频日本电影| 欧美日韩黄片免| 女人高潮潮喷娇喘18禁视频| 日本欧美视频一区| 成人午夜精彩视频在线观看| 午夜福利,免费看| 欧美激情极品国产一区二区三区| 国产在线免费精品| 99国产精品一区二区蜜桃av | 亚洲五月婷婷丁香| 黄色片一级片一级黄色片| 操美女的视频在线观看| 午夜久久久在线观看| 免费一级毛片在线播放高清视频 | 久久国产精品影院| 欧美成人午夜精品| 悠悠久久av| 超碰97精品在线观看| 少妇粗大呻吟视频| 国产又爽黄色视频| 黄片小视频在线播放| 国产av精品麻豆| 亚洲激情五月婷婷啪啪| 欧美人与性动交α欧美精品济南到| 视频区图区小说| 精品欧美一区二区三区在线| 我的亚洲天堂| 人妻人人澡人人爽人人| 少妇 在线观看| 国产在线一区二区三区精| 天天躁狠狠躁夜夜躁狠狠躁| 热99久久久久精品小说推荐| 国产一区二区激情短视频 | 午夜免费男女啪啪视频观看| 19禁男女啪啪无遮挡网站| 国产精品久久久久久人妻精品电影 | 在线观看免费日韩欧美大片| 久久这里只有精品19| 亚洲欧洲日产国产| 亚洲精品久久午夜乱码| 一级毛片我不卡| 日本欧美国产在线视频| 亚洲国产欧美日韩在线播放| 国产在线观看jvid| 精品少妇黑人巨大在线播放| 成人亚洲欧美一区二区av| 亚洲一码二码三码区别大吗| 国产熟女午夜一区二区三区| 日韩中文字幕视频在线看片| 久久人人爽av亚洲精品天堂| 天天躁日日躁夜夜躁夜夜| 日韩中文字幕欧美一区二区 | 久久99精品国语久久久| 久9热在线精品视频| 男女国产视频网站| 国精品久久久久久国模美| 国产精品熟女久久久久浪| 一区二区三区精品91| 亚洲天堂av无毛| 一级片'在线观看视频| 精品少妇久久久久久888优播| 亚洲精品av麻豆狂野| 曰老女人黄片| 亚洲 国产 在线| 国产精品久久久人人做人人爽| 大香蕉久久网| 精品免费久久久久久久清纯 | 超色免费av| 日韩大片免费观看网站| 亚洲欧洲日产国产| 高清不卡的av网站| 欧美人与性动交α欧美软件| 老汉色∧v一级毛片| 欧美在线黄色| 在线精品无人区一区二区三| 老鸭窝网址在线观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产色婷婷电影| 一本大道久久a久久精品| 国产欧美日韩一区二区三区在线| 国产xxxxx性猛交| 国产精品一区二区在线观看99| 国产伦人伦偷精品视频| 老司机午夜十八禁免费视频| 高清不卡的av网站| 一区二区三区激情视频| 在线观看人妻少妇| 国产又色又爽无遮挡免| 99国产综合亚洲精品| 亚洲天堂av无毛| 成人影院久久| 夜夜骑夜夜射夜夜干| 久久久久久人人人人人| 天天影视国产精品| 人人妻人人澡人人爽人人夜夜| 男女之事视频高清在线观看 | 又粗又硬又长又爽又黄的视频| 精品一区二区三区av网在线观看 | av网站免费在线观看视频| 亚洲人成电影免费在线| 日本猛色少妇xxxxx猛交久久| 午夜日韩欧美国产| 18禁观看日本| 午夜免费男女啪啪视频观看| 国产精品一区二区在线观看99| 1024视频免费在线观看| 狂野欧美激情性xxxx| 黄色毛片三级朝国网站| 亚洲欧美日韩另类电影网站| 国产xxxxx性猛交| 国产精品av久久久久免费| 日本av免费视频播放| 免费av中文字幕在线| 2021少妇久久久久久久久久久| 丝袜在线中文字幕| 91精品三级在线观看| 欧美亚洲 丝袜 人妻 在线| 一级片'在线观看视频| 9191精品国产免费久久| 国产精品久久久久久精品古装| 考比视频在线观看| 国产男女超爽视频在线观看| 久久免费观看电影| 亚洲综合色网址| netflix在线观看网站| 日韩制服丝袜自拍偷拍| 欧美人与善性xxx| 1024视频免费在线观看| av一本久久久久| 在线 av 中文字幕| 乱人伦中国视频| 久久久久精品人妻al黑| 色精品久久人妻99蜜桃| 国产精品亚洲av一区麻豆| 中文字幕高清在线视频| 国产一级毛片在线| 午夜免费男女啪啪视频观看| 免费观看人在逋| 亚洲av男天堂| 91老司机精品| 国产视频一区二区在线看| 美女午夜性视频免费| 一级毛片黄色毛片免费观看视频| 一区在线观看完整版| 亚洲国产中文字幕在线视频| 王馨瑶露胸无遮挡在线观看| 欧美日韩黄片免| 中文字幕人妻丝袜一区二区| 久久精品成人免费网站| 国产激情久久老熟女| a级毛片黄视频| 精品国产一区二区三区四区第35| 亚洲精品美女久久久久99蜜臀 | av电影中文网址| av天堂在线播放| 色网站视频免费| 午夜老司机福利片| 日日夜夜操网爽| 国产高清视频在线播放一区 | 波多野结衣av一区二区av| 老司机影院毛片| 国产熟女午夜一区二区三区| 青春草亚洲视频在线观看| 在现免费观看毛片| 乱人伦中国视频| 视频区欧美日本亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲欧美精品永久| 我要看黄色一级片免费的| 久久ye,这里只有精品| 久久精品久久精品一区二区三区| 一级a爱视频在线免费观看| 午夜福利视频在线观看免费| 国产男女超爽视频在线观看| 成人手机av| 国产成人av激情在线播放| 一级黄片播放器| 另类亚洲欧美激情| 久久久久久久久久久久大奶| 国产男女内射视频| 老司机深夜福利视频在线观看 | 成年av动漫网址| 国产欧美亚洲国产| 日本一区二区免费在线视频| 欧美日韩国产mv在线观看视频| av在线app专区| 欧美激情高清一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品影院| 国产在线视频一区二区| 亚洲国产av影院在线观看| 一区二区三区乱码不卡18| 午夜视频精品福利| 免费观看av网站的网址| 亚洲精品乱久久久久久| 黄色毛片三级朝国网站| 乱人伦中国视频| 纵有疾风起免费观看全集完整版| 欧美中文综合在线视频| 午夜激情av网站| 亚洲人成网站在线观看播放| 亚洲av国产av综合av卡| 国产视频首页在线观看| 最近手机中文字幕大全| 手机成人av网站| 久久毛片免费看一区二区三区| 在线精品无人区一区二区三| 两性夫妻黄色片| 精品高清国产在线一区| 狠狠精品人妻久久久久久综合| av有码第一页| 啦啦啦在线观看免费高清www| 欧美日韩精品网址| 亚洲国产精品999| 国产欧美日韩综合在线一区二区| 久久99热这里只频精品6学生| 婷婷成人精品国产| 黄色视频不卡| 中文字幕色久视频| 亚洲熟女毛片儿| 九草在线视频观看| 国产欧美日韩综合在线一区二区| 黄色毛片三级朝国网站| 新久久久久国产一级毛片| 亚洲黑人精品在线| 黄片播放在线免费| av天堂在线播放| 1024视频免费在线观看| 亚洲色图综合在线观看| 国产野战对白在线观看| 亚洲av欧美aⅴ国产| 亚洲国产精品一区三区| 嫩草影视91久久| 国产女主播在线喷水免费视频网站| √禁漫天堂资源中文www| 妹子高潮喷水视频| 色94色欧美一区二区| av网站在线播放免费| 亚洲欧美激情在线| 啦啦啦 在线观看视频| 狂野欧美激情性xxxx| 亚洲欧美激情在线| 国产男女超爽视频在线观看| 国产成人免费观看mmmm| 超碰成人久久| 久久久欧美国产精品| 色精品久久人妻99蜜桃| 看免费av毛片| tube8黄色片| 色精品久久人妻99蜜桃| 丁香六月欧美| 国产主播在线观看一区二区 | 亚洲欧洲日产国产| 亚洲精品在线美女| 久久人妻福利社区极品人妻图片 | 久久狼人影院| 如日韩欧美国产精品一区二区三区| 中文字幕av电影在线播放| 无遮挡黄片免费观看| 十八禁人妻一区二区| 丰满少妇做爰视频| 精品一品国产午夜福利视频| 午夜老司机福利片| 亚洲国产精品成人久久小说| 亚洲视频免费观看视频| 国产一区二区在线观看av| 一边摸一边做爽爽视频免费| 大香蕉久久网| 国产精品人妻久久久影院| 丰满饥渴人妻一区二区三| 国产在线免费精品| 又粗又硬又长又爽又黄的视频| 亚洲三区欧美一区| 天天添夜夜摸| 欧美少妇被猛烈插入视频| 久久人人97超碰香蕉20202| 老鸭窝网址在线观看| 国产1区2区3区精品| 巨乳人妻的诱惑在线观看| 麻豆国产av国片精品| 777久久人妻少妇嫩草av网站| 在线精品无人区一区二区三| xxx大片免费视频| 秋霞在线观看毛片| 国产一区二区三区综合在线观看| 久久久久视频综合| 晚上一个人看的免费电影| 久久精品国产亚洲av高清一级| 欧美精品av麻豆av| 日本91视频免费播放| 亚洲五月婷婷丁香| 老司机亚洲免费影院| 久久鲁丝午夜福利片| 精品欧美一区二区三区在线| 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| 在线观看国产h片| 女人被躁到高潮嗷嗷叫费观| 国产成人精品在线电影| 久久午夜综合久久蜜桃| av在线app专区| 免费观看a级毛片全部| 中文字幕高清在线视频| 巨乳人妻的诱惑在线观看| 国产熟女午夜一区二区三区| 丝袜脚勾引网站| 国产男女超爽视频在线观看| 中国国产av一级| 午夜福利视频在线观看免费| 麻豆国产av国片精品| 9热在线视频观看99| 亚洲人成77777在线视频| 久久人人爽av亚洲精品天堂| 精品一区二区三卡| 天堂俺去俺来也www色官网| 久久女婷五月综合色啪小说| 久久热在线av| avwww免费| 七月丁香在线播放| 精品国产国语对白av| 精品少妇黑人巨大在线播放| 亚洲成国产人片在线观看| 啦啦啦啦在线视频资源| 天天躁夜夜躁狠狠久久av| 丝袜美足系列| 亚洲av美国av| 亚洲国产欧美在线一区| 另类精品久久| 国产男女超爽视频在线观看| 99热国产这里只有精品6| 青春草亚洲视频在线观看| 国精品久久久久久国模美| 美女福利国产在线| 人妻人人澡人人爽人人| 丰满饥渴人妻一区二区三| av天堂在线播放| 欧美黑人欧美精品刺激| 国产色视频综合| 欧美乱码精品一区二区三区| 日韩,欧美,国产一区二区三区| 婷婷色综合大香蕉| 侵犯人妻中文字幕一二三四区| 看免费成人av毛片| 久久久精品94久久精品| 乱人伦中国视频| 制服诱惑二区| 久久国产精品影院| 亚洲欧美一区二区三区黑人| 91精品国产国语对白视频| 国产精品久久久人人做人人爽| 日韩av不卡免费在线播放| 亚洲国产欧美网| 久久鲁丝午夜福利片| 精品一区在线观看国产| 精品亚洲乱码少妇综合久久| 嫁个100分男人电影在线观看 | 国产91精品成人一区二区三区 | 国产成人系列免费观看| 国产熟女午夜一区二区三区| 777久久人妻少妇嫩草av网站| 黄色一级大片看看| 欧美大码av| 国产精品香港三级国产av潘金莲 | 久久久精品区二区三区| 交换朋友夫妻互换小说| netflix在线观看网站| 日日爽夜夜爽网站| 一二三四在线观看免费中文在| 国产成人一区二区在线| av网站在线播放免费| 日本五十路高清| 一区在线观看完整版| 欧美成人午夜精品| 亚洲一区二区三区欧美精品| www日本在线高清视频| 国产三级黄色录像| 精品一区在线观看国产| 激情视频va一区二区三区| 国产精品久久久久久人妻精品电影 | 亚洲少妇的诱惑av| 99精品久久久久人妻精品| 午夜两性在线视频| 大型av网站在线播放| av不卡在线播放| 夫妻性生交免费视频一级片| 免费不卡黄色视频| 999精品在线视频| 两性夫妻黄色片| 又紧又爽又黄一区二区| 国产成人精品久久二区二区91| 久久青草综合色| 另类亚洲欧美激情| 婷婷丁香在线五月| 97在线人人人人妻| 免费在线观看黄色视频的| 久久久精品区二区三区| 一级片免费观看大全| 老司机在亚洲福利影院| 精品免费久久久久久久清纯 | 精品一区二区三区av网在线观看 | 十八禁网站网址无遮挡| 男女国产视频网站| 涩涩av久久男人的天堂| 人人妻人人澡人人看| 妹子高潮喷水视频| 国产精品偷伦视频观看了| 又紧又爽又黄一区二区| 伊人亚洲综合成人网| 美女国产高潮福利片在线看| 国产免费福利视频在线观看| 欧美日韩视频高清一区二区三区二| 国产日韩欧美在线精品| 黄色片一级片一级黄色片| 久久久久久久大尺度免费视频| 日韩一本色道免费dvd| 一级a爱视频在线免费观看| 国产精品 国内视频| 亚洲精品日本国产第一区| 波野结衣二区三区在线| 高清不卡的av网站| 黄片小视频在线播放| 国产精品欧美亚洲77777| 天天影视国产精品| 少妇的丰满在线观看| 男女午夜视频在线观看| 亚洲国产av影院在线观看| 嫩草影视91久久| 日韩一卡2卡3卡4卡2021年| 老熟女久久久| 日本vs欧美在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 美女大奶头黄色视频| 欧美国产精品va在线观看不卡| 亚洲av电影在线进入| 男男h啪啪无遮挡| 亚洲人成77777在线视频| 99久久99久久久精品蜜桃| 国产av精品麻豆| 中文字幕精品免费在线观看视频| 波野结衣二区三区在线| 精品国产国语对白av| 国产免费现黄频在线看| 99久久99久久久精品蜜桃| 大陆偷拍与自拍| 国产免费一区二区三区四区乱码| 国产亚洲欧美在线一区二区| 亚洲成人手机| 亚洲人成77777在线视频| 亚洲中文日韩欧美视频| 亚洲av电影在线进入| 亚洲国产最新在线播放| 亚洲成色77777| 日日摸夜夜添夜夜爱| 国产高清视频在线播放一区 | 一级片免费观看大全| 老司机在亚洲福利影院| 十八禁人妻一区二区| 黄色 视频免费看| 一区二区av电影网| 天天影视国产精品| svipshipincom国产片| 亚洲成av片中文字幕在线观看| 国产99久久九九免费精品| 国产男女超爽视频在线观看| 美国免费a级毛片| cao死你这个sao货| 建设人人有责人人尽责人人享有的| 免费人妻精品一区二区三区视频| 侵犯人妻中文字幕一二三四区| 久久久久久久久久久久大奶| 国产黄色视频一区二区在线观看| 国产成人精品久久二区二区91| av不卡在线播放| 色精品久久人妻99蜜桃| 日韩一区二区三区影片| 色精品久久人妻99蜜桃| 亚洲五月色婷婷综合| 免费人妻精品一区二区三区视频| 国产视频一区二区在线看| 欧美日韩视频精品一区| 亚洲精品国产av蜜桃| 国产极品粉嫩免费观看在线| 亚洲精品一区蜜桃| 菩萨蛮人人尽说江南好唐韦庄| 亚洲伊人色综图| 欧美黄色淫秽网站| 精品一区二区三卡| 色综合欧美亚洲国产小说| 美女中出高潮动态图| 精品一区在线观看国产| 亚洲成人国产一区在线观看 | 男女无遮挡免费网站观看| 欧美国产精品一级二级三级| 18禁国产床啪视频网站| 高清不卡的av网站| 色网站视频免费| av电影中文网址| 亚洲国产欧美一区二区综合| 国产视频一区二区在线看| 一区在线观看完整版| 人人妻,人人澡人人爽秒播 | av国产精品久久久久影院| 波野结衣二区三区在线| 亚洲自偷自拍图片 自拍| 国产女主播在线喷水免费视频网站| 成年人午夜在线观看视频| 一级,二级,三级黄色视频| 国产无遮挡羞羞视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩亚洲高清精品| 黄色片一级片一级黄色片| 国产亚洲精品第一综合不卡| 午夜激情久久久久久久| 精品国产一区二区三区四区第35| www.精华液| 久久av网站| 国产免费一区二区三区四区乱码| 国产有黄有色有爽视频| www日本在线高清视频| 欧美日韩成人在线一区二区| 成年av动漫网址|