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

    淠河水質(zhì)時(shí)空變化特征及驅(qū)動(dòng)因素分析

    2024-07-04 10:49:58朱沉靜譚茜郭志陽汪可欣吳師
    人民長江 2024年6期
    關(guān)鍵詞:淠河新安渡口

    朱沉靜 譚茜 郭志陽 汪可欣 吳師

    收稿日期:2023-11-23;接受日期:2024-04-03

    基金項(xiàng)目:引江濟(jì)淮工程(安徽段)輸水沿線水源保護(hù)區(qū)劃分、重大突發(fā)性水污染事件預(yù)警系統(tǒng)和水源地保護(hù)規(guī)劃等課題研究項(xiàng)目(YJJH-ZT-ZX-20230706545);國家水體污染控制與治理科技重大專項(xiàng)(2017X07603-005)

    作者簡介:朱沉靜,女,助理工程師,碩士,研究方向?yàn)樗Y源開發(fā)利用與保護(hù)。E-mail:chenj_zhu@cug.edu.cn

    Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

    文章編號(hào):1001-4179(2024) 06-0030-07

    引用本文:朱沉靜,譚茜,郭志陽,等.

    淠河水質(zhì)時(shí)空變化特征及驅(qū)動(dòng)因素分析

    [J].人民長江,2024,55(6):30-36.

    摘要:為全面了解近年來淠河流域水質(zhì)污染狀況及時(shí)空變化規(guī)律,基于2018~2022年淠河干流水質(zhì)考核斷面pH、溶解氧(DO)、化學(xué)需氧量(COD)、五日生化需氧量(BOD5)、氨氮(NH3-N)、總氮(TN)及總磷(TP)監(jiān)測數(shù)據(jù),利用水質(zhì)指數(shù)(WQI)法對(duì)河流水質(zhì)進(jìn)行綜合評(píng)價(jià),并結(jié)合相關(guān)性分析識(shí)別影響水質(zhì)變化的關(guān)鍵水質(zhì)指標(biāo)。結(jié)果表明:① 淠河干流水質(zhì)總體較好,WQI均值為(75.94±7.47),處于“良好”水平,但非汛期氨氮濃度為0.020~1.480 mg/L,負(fù)荷較高。② 空間分布上水質(zhì)優(yōu)劣排序?yàn)闄M排頭>大店崗>窯崗嘴>新安渡口。③ 在時(shí)間上干流汛期(6~9月)水質(zhì)明顯優(yōu)于非汛期(1~5月和10~12月)。④ 結(jié)合相關(guān)性分析得出,研究區(qū)內(nèi)主要污染指標(biāo)為NH3-N、TN和COD;橫排頭—窯崗嘴段及新安渡口—大店崗段流域土地利用類型耕地占比高,水體受化肥施用和畜禽養(yǎng)殖業(yè)產(chǎn)生的COD、TN和NH3-N污染;窯崗嘴—新安渡口段水質(zhì)受城鎮(zhèn)生活污水直排和污水廠尾水排放影響,NH3-N污染嚴(yán)重;汛期污染程度伴隨集中降雨強(qiáng)度逐步減弱。為保護(hù)淠河水生態(tài)環(huán)境,建議從增強(qiáng)流域自凈能力、加強(qiáng)農(nóng)業(yè)面源污染防控、完善污水截污管網(wǎng)系統(tǒng)等方面加強(qiáng)措施。

    關(guān)? 鍵? 詞:水質(zhì); WQI法; 變化趨勢; 影響因素; 淠河

    中圖法分類號(hào): X824

    文獻(xiàn)標(biāo)志碼: A

    DOI:10.16232/j.cnki.1001-4179.2024.06.005

    0? 引 言

    河流是孕育城市文明的搖籃,其水質(zhì)安全與否直接關(guān)系到居民的生存發(fā)展。隨著城市擴(kuò)張和人口增長,各類水體污染源帶來的威脅日益嚴(yán)峻[1-3]。針對(duì)河流開展水質(zhì)評(píng)估有助于了解水質(zhì)污染狀況以及時(shí)空變化規(guī)律,對(duì)河流的保護(hù)及管理工作具有重要意義。利用水質(zhì)指標(biāo)定量評(píng)價(jià)水質(zhì)可以采用多種方法,包括主成分分析法、人工神經(jīng)網(wǎng)絡(luò)評(píng)價(jià)法、單因子評(píng)價(jià)法以及水質(zhì)指數(shù)(Water Quality Index,WQI)法等。上述方法中,水質(zhì)指數(shù)法不僅能夠?qū)⒋罅繌?fù)雜的水質(zhì)指標(biāo)轉(zhuǎn)變?yōu)橐粋€(gè)反映水質(zhì)整體狀況的WQI評(píng)分,還能客觀顯示污染因子對(duì)水體污染的貢獻(xiàn)程度。相較其他方法,WQI法具有計(jì)算過程簡單、評(píng)價(jià)結(jié)果直觀等優(yōu)點(diǎn),因而被廣泛應(yīng)用于評(píng)估地表水水質(zhì),尤其是在確定水質(zhì)時(shí)空變化趨勢方面[4-6]。

    淠河是淮河中游南岸的一條較大支流,承擔(dān)著重要的區(qū)域供水及灌溉功能。近年來受強(qiáng)烈人類活動(dòng)影響,僅2019年淠河下游新安渡口斷面水質(zhì)就有8次超標(biāo)[7],這給當(dāng)?shù)爻青l(xiāng)生活和工農(nóng)業(yè)生產(chǎn)帶來極大不便。由于淠河水質(zhì)受到污染物入河量、降水頻率、降水強(qiáng)度、水流交換速度及水體周邊環(huán)境等諸多因素的影

    響,其在空間和時(shí)間尺度上均呈現(xiàn)較強(qiáng)的特異性,難以制定針對(duì)性的水質(zhì)保障措施。目前,針對(duì)淠河水質(zhì)水期變化趨勢[8-9]和淠河城區(qū)段污染物入河結(jié)構(gòu)[7]已有一定的研究,但這些研究中使用的評(píng)價(jià)方法多為GB 3838-2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中的單因子評(píng)價(jià)方法,因其“單項(xiàng)否決”的特點(diǎn),難以識(shí)別不同污染因子對(duì)水質(zhì)的影響。此外,針對(duì)多因素驅(qū)動(dòng)下的淠河水質(zhì)演變機(jī)制的研究也較少。

    因此,本研究在探明淠河重要水質(zhì)指標(biāo)時(shí)空分布特征的基礎(chǔ)上,利用WQI法對(duì)干流水質(zhì)整體狀況進(jìn)行評(píng)價(jià),并結(jié)合相關(guān)性分析識(shí)別不同時(shí)空條件下影響WQI指數(shù)的主要參數(shù),以期豐富對(duì)淠河水質(zhì)演變規(guī)律及主要驅(qū)動(dòng)因素的認(rèn)識(shí),為當(dāng)?shù)厮|(zhì)保障和水污染治理提供科學(xué)依據(jù)。

    1? 研究區(qū)概況

    淠河發(fā)源于岳西和金寨縣境內(nèi)的大別山北麓,由南向北流經(jīng)岳西、霍山、金寨、六安、霍邱、壽縣等縣(市)后,在正陽關(guān)匯入淮河,全長253 km,流域總面積6 000 km2。本研究淠河河段范圍為干流橫排頭樞紐到大店崗國考斷面,長度約為118 km。其中淠河橫排頭—窯崗嘴斷面為上游農(nóng)村段,河道長度約29 km;窯崗嘴—新安渡口斷面為城區(qū)段,河道長度約19 km;新安渡口—大店崗斷面為下游農(nóng)村段,河道長度約70 km(圖1)。

    淠河流域?qū)賮啛釒П辈窟吘壍臇|亞季風(fēng)氣候,年平均氣溫14.2~15.4 ℃,夏季高溫多雨,冬季低溫干燥。多年平均降雨量1 334 mm,且具有明顯的季節(jié)性,降雨主要集中分布在6~9月,約占全年的67%[10],其中6~7月多為渦切變型暴雨(俗稱梅雨)形成主汛期,8~9月為臺(tái)風(fēng)型暴雨形成伏汛。流域暴雨中心多發(fā)生在佛子嶺、響洪甸水庫上游。

    2? 數(shù)據(jù)來源和研究方法

    2.1? 數(shù)據(jù)來源

    以6~9月作為汛期,10月至次年5月作為非汛期。淠河水質(zhì)數(shù)據(jù)來自六安市生態(tài)環(huán)境局,包含橫排頭、窯崗嘴2個(gè)省級(jí)考核斷面及新安渡口、大店崗2個(gè)國家級(jí)考核斷面(斷面位置如圖1所示)。2018年1月至2022年12月逐月實(shí)測酸堿度(pH)、溶解氧(DO)、化學(xué)需氧量(COD)、五日生化需氧量(BOD5)、氨氮(NH3-N)、總磷(TP)、總氮(TN),共7項(xiàng)指標(biāo),監(jiān)測頻次為每月1次。研究中涉及的降雨量數(shù)據(jù)來源于2019~2023年《安徽統(tǒng)計(jì)年鑒》。

    2.2? 研究方法

    2.2.1? 水質(zhì)評(píng)價(jià)方法

    采用Pesce等提出的水質(zhì)指數(shù)(WQI)計(jì)算方法[11],如公式(1)所示:WQI=ni=1CiPini=1Pi

    (1)

    式中:n為水質(zhì)評(píng)價(jià)過程中所有指標(biāo)的數(shù)量;Ci為指標(biāo)i實(shí)測濃度的標(biāo)準(zhǔn)化賦分?jǐn)?shù)值,參考GB 3838-2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》,從0~100劃分為11個(gè)級(jí)別,根據(jù)水質(zhì)指標(biāo)的濃度進(jìn)行賦值;Pi為指標(biāo)i的權(quán)重,參考水質(zhì)指標(biāo)對(duì)水生生物的重要性進(jìn)行劃分,其中最大權(quán)重為4,最小值為1[12-15],如表1所列。

    WQI的取值范圍為0~100,值越高代表水質(zhì)越好。根據(jù)WQI得分,將水質(zhì)分為5個(gè)等級(jí):優(yōu)秀(90~100)、良好(70~90)、中等(50~70)、一般(25~ 50)和差(0~25)[16]。

    2.2.2? 相關(guān)性分析

    相關(guān)性分析是一種統(tǒng)計(jì)分析方法,用于衡量兩個(gè)或多個(gè)變量之間的關(guān)系強(qiáng)度和方向。這種關(guān)系可以用相關(guān)系數(shù)來量化,最常見的是Pearson相關(guān)系數(shù)(線性相關(guān)),其值介于-1與1之間,負(fù)值表示負(fù)相關(guān),正值表示正相關(guān)。在顯著性的前提下,絕對(duì)值越大,相關(guān)性越強(qiáng)。本研究利用SPSS 26.0軟件計(jì)算淠河考核斷面水質(zhì)指標(biāo)和WQI指數(shù)之間的Pearson相關(guān)系數(shù)。

    3? 結(jié)果與討論

    3.1? 水質(zhì)指標(biāo)時(shí)空變化特征

    3.1.1? 水質(zhì)指標(biāo)空間差異特征

    如表2所列,2018~2022年,各監(jiān)測斷面COD年均濃度值從大到小依次為窯崗嘴>新安渡口>大店崗>橫排頭。各斷面BOD5、TN及NH3-N年均濃度值從大到小則均表現(xiàn)為新安渡口>大店崗>窯崗嘴>橫排頭。已有研究表明,河流水質(zhì)與流域土地利用方式密切相關(guān)[17]。根據(jù)《六安市十四五水土保持規(guī)劃》,2020年六安市縣(區(qū))耕地占土地總面積比例分別為霍山縣(8.16%)、金安區(qū)(49.97%)、裕安區(qū)(46.27%)、霍邱縣(67.39%)。橫排頭—窯崗嘴屬金安、裕安區(qū)境內(nèi),流域范圍內(nèi)分布有大量耕地。土壤中未被吸收的殘留無機(jī)物和農(nóng)藥等,易隨降雨形成的徑流進(jìn)入附近水體;此外,耕地灌溉也增加了土壤水分,易于土壤氮及化肥進(jìn)入地表水體,最終導(dǎo)致窯崗嘴斷面COD濃度較高。窯崗嘴—新安渡口段兩岸為六安主城區(qū),該段沿線布有多個(gè)入河排污口,水體承接了大量的生活生產(chǎn)廢水及城鎮(zhèn)地表徑流,故新安渡口斷面TN及NH3-N濃度偏高。TP濃度平均值各斷面則無明顯變化。已有研究結(jié)果表明,磷和其他磷離子均與泥沙具有較強(qiáng)的親和力[18]。橫排頭以下河道進(jìn)入丘陵和平原區(qū),河面變寬,比降變緩,水流速度較小,利于懸移質(zhì)泥沙沉淀,進(jìn)而提升了河道吸附和固定磷的能力,因此整體上淠河干流總磷濃度波動(dòng)幅度小。

    3.1.2? 水質(zhì)指標(biāo)時(shí)間變化特征

    由于淠河流域降水年內(nèi)分布不均勻,汛期與非汛期水質(zhì)指標(biāo)存在明顯差異,根據(jù)統(tǒng)計(jì),六安市2020年汛期降雨量占年降雨總量的70.46%,因此,需對(duì)水質(zhì)變化情況以年內(nèi)時(shí)間尺度加以分析。

    從圖2(a)可以看出,汛期(6~9月)DO相較非汛期(10月至次年5月)呈下降趨勢,推測可能是由于汛期水溫偏高,導(dǎo)致溶解氧飽和值下降;同時(shí),汛期較大的徑流量和降雨量易于有機(jī)物等還原性物質(zhì)進(jìn)入水體,隨著溫度的升高,微生物降解有機(jī)質(zhì)消耗溶解氧的活動(dòng)更加劇烈。前人研究結(jié)果指出大氣入滲和水生植物光合作用釋放是水體DO的主要來源[19]。汛期COD濃度普遍高于非汛期,呈現(xiàn)出與DO相反的變化趨勢(圖2(c)),進(jìn)一步佐證:相較非汛期,汛期水體受有機(jī)物污染和還原性無機(jī)物污染程度更為嚴(yán)重。

    水體NH3-N、TN、TP及BOD5濃度則普遍表現(xiàn)為非汛期高于汛期,且非汛期氨氮指標(biāo)負(fù)荷較高(圖2(b)~(f)),表明非汛期硝化細(xì)菌作用下的水中含氮類有機(jī)物的分解活躍。這可能是由于橫排頭樞紐等攔河構(gòu)筑物的建成攔截了大部分上游來水,使得水體流動(dòng)性變小,不利于營養(yǎng)鹽類的輸出。此外,汛期集中降水產(chǎn)生的大量地表徑流注入水體,在稀釋污染物的同時(shí),也在一定程度上改善了河道的水動(dòng)力條件,加快水體更新,抑制了富營養(yǎng)化。

    3.2? WQI水質(zhì)評(píng)價(jià)

    研究區(qū)2018~2022年各監(jiān)測斷面WQI水質(zhì)指數(shù)評(píng)價(jià)結(jié)果如圖3所示。

    總體而言,淠河WQI值處于“良好”(70~90)及“中等”(50~70)兩個(gè)等級(jí)區(qū)間,多年WQI均值為(75.94±7.47)。由圖3(a)可知,2019年1~12月各考核斷面WQI“良好” 占比相較其余年份最低。近兩年來,淠河水質(zhì)持續(xù)提升,2022年WQI“良好”占比相較前4 a最高。

    從時(shí)間層面上看(圖3(b)),汛期各斷面WQI等級(jí)“良好”占比,相較非汛期有所上升。由于部分非汛期時(shí)段存在WQI等級(jí)為“中等”的情況,因此若以WQI等級(jí)“良好”作為判斷該監(jiān)測斷面水質(zhì)是否達(dá)標(biāo)的依據(jù),則汛期的水質(zhì)達(dá)標(biāo)率要明顯高于非汛期。這可能是由于汛期伴隨的大量雨水的稀釋作用,使得水體中污染因子濃度降低。從空間層面上看(圖3(b)),各考核斷面WQI值由高到低依次為橫排頭>大店崗>窯崗嘴>新安渡口。橫排頭斷面WQI

    值相較其他斷面最高。這是由于橫排頭斷面以上為東西淠河交匯口,承接上游佛子嶺水庫及響洪甸水庫來水,水庫水質(zhì)常年保持Ⅱ類標(biāo)準(zhǔn),因此該斷面水質(zhì)較好。而新安渡口斷面WQI值為4個(gè)斷面中最低,推測是由于城鎮(zhèn)污水處理廠尾水排放及大雁河、鳳凰河、蘇大堰河等5條水質(zhì)較差的支流匯入,使大量污染物進(jìn)入窯崗嘴—新安渡口段;加之區(qū)段內(nèi)新安橡膠壩和城

    北橡膠壩的設(shè)立,使得壩體之間水流速變小,對(duì)污染物稀釋作用較弱。上述多種因素共同作用最終導(dǎo)致新安渡口斷面水質(zhì)較差。

    3.3? 水質(zhì)時(shí)空變化驅(qū)動(dòng)因素分析

    由前文分析可知,淠河干流橫排頭、窯崗嘴、新安渡口及大店崗斷面的水質(zhì)指標(biāo)及WQI指數(shù)具有較大的差異性。因此,將2018年1月至2022年12月上述各考核斷面WQI值與7項(xiàng)水質(zhì)指標(biāo)進(jìn)行Pearson相關(guān)性分析很有必要。通過比較其相關(guān)系數(shù)值大小,以查明影響水質(zhì)的主要污染因子。

    結(jié)果如圖4所示:總體來看,氮磷營養(yǎng)鹽指標(biāo)(TN、NH3-N、TP)及有機(jī)化合物指標(biāo)(COD、BOD5)與WQI之間呈現(xiàn)顯著性負(fù)相關(guān)(P<0.05);溶解氧指標(biāo)與WQI呈顯著性正相關(guān)(P<0.05)。故影響淠河水質(zhì)的污染物類型可大致分為2大類:營養(yǎng)鹽及有機(jī)物質(zhì)。

    淠河橫排頭—窯崗嘴段內(nèi)建有安徽省六安市淠河國家濕地公園,濕地系統(tǒng)中的水生植物及微生物對(duì)氮磷等營養(yǎng)物質(zhì)具有良好的去除作用[20-21],因此水質(zhì)與其他段相比整體偏好。如圖4(a)~(b)所示,NH3-N和COD與WQI指數(shù)之間呈現(xiàn)顯著負(fù)相關(guān)性,說明該段影響水質(zhì)的主要污染指標(biāo)為NH3-N及COD。前人研究結(jié)果證實(shí),農(nóng)村化肥施用和畜禽糞便,是淮河流域氮污染的主要來源[22]。根據(jù)《安徽統(tǒng)計(jì)年鑒》,2020年六安市化肥施用共計(jì)為162 547 t,其中氮肥46 248 t,磷肥14 962 t。橫排頭—窯崗嘴段耕地居多,且沿河分布有少量的排污口(圖1),由此推測除內(nèi)源輸入外,COD、NH3-N主要來源為隨地表徑流進(jìn)入水體的化肥、畜禽糞便及農(nóng)田中殘留有機(jī)肥升溫發(fā)酵形成有機(jī)質(zhì),其次為管網(wǎng)排水。

    相關(guān)研究表明,面源污染物濃度受降水量和降水強(qiáng)度影響,在相同的污染物累積量條件下,表現(xiàn)出隨降雨量及降雨強(qiáng)度的增強(qiáng),徑流中污染物濃度逐漸降低的趨勢[23]。從降雨量與WQI的關(guān)系圖(圖5)可以發(fā)現(xiàn),汛期水質(zhì)考核斷面WQI指數(shù)值較高的點(diǎn)通常伴隨較大的降雨量,而非汛期存在部分時(shí)段WQI值較低且降雨量少,說明汛期降雨在一定程度上改善了淠河水質(zhì)。

    淠河窯崗嘴—新安渡口段內(nèi)共設(shè)置有5座污水處理廠及多處生活污水入河排污口(圖1),水質(zhì)較差。2020年六安市淠河水環(huán)境綜合治理項(xiàng)目在該河段范圍內(nèi)啟動(dòng),包含灘岸整治、濱岸帶植被修復(fù)及支流水環(huán)境治理等工程措施,水質(zhì)整體呈現(xiàn)回升向好趨勢。如圖4(c)所示,新安渡口斷面NH3-N和TN與WQI指數(shù)之間呈現(xiàn)顯著負(fù)相關(guān)性,說明該段影響水質(zhì)的主要污染指標(biāo)為NH3-N及TN。已有研究結(jié)果表明,窯崗嘴—新安渡口段水體氮的主要來源為污水廠尾水及城鎮(zhèn)污水,其直接入河的NH3-N含量約為275.06 t/a,約占77.10%[7]。因此,城區(qū)污水入河量的控制是改善該段水環(huán)境質(zhì)量的重要任務(wù)。

    如圖4(d)所示,大店崗斷面NH3-N、TN和COD與WQI指數(shù)之間呈現(xiàn)顯著負(fù)相關(guān)性,說明NH3-N、TN和COD為影響新安渡口—大店崗段的主要水質(zhì)指標(biāo)。淠河新安渡口—大店崗段屬霍邱及壽縣境內(nèi),河勢平坦,其中霍邱縣2020年土地利用類型耕地占比約為67.39%;淮南市2019年耕地33.37萬hm2(500.56萬畝)(第三次全國國土調(diào)查),約88.34%的耕地分布于壽縣、鳳臺(tái)縣、潘集區(qū),與上游橫排頭—窯崗嘴段相似流域內(nèi)耕地占比高。該段沿線也布置有入河排污口(圖1),因此農(nóng)業(yè)面源污染和城鎮(zhèn)管網(wǎng)點(diǎn)源污染是該段水體污染物的主要外源輸入。

    4? 結(jié)論與建議

    本文重點(diǎn)分析了2018~2022年淠河干流各考核斷面水質(zhì)指標(biāo)和WQI值的時(shí)空變化特征,以及它們之間的相關(guān)性。研究表明:① 近年淠河干流WQI均值為(75.94±7.47),水質(zhì)處于“良好”水平,各水質(zhì)指標(biāo)總體較好。② 研究區(qū)WQI指數(shù)時(shí)空變異性大,時(shí)間上汛期(6~9月)水質(zhì)明顯優(yōu)于非汛期(1~5月和10月至次年12月),空間上水質(zhì)優(yōu)劣排序?yàn)闄M排頭>大店崗>窯崗嘴>新安渡口。③ 結(jié)合相關(guān)性分析,各考核斷面NH3-N、TN和COD與WQI指數(shù)之間呈現(xiàn)顯著負(fù)相關(guān)性,表明研究區(qū)內(nèi)主要污染指標(biāo)為NH3-N、TN和COD。④ 橫排頭—窯崗嘴段及新安渡口—大店崗段流域土地利用類型耕地占比高,水體主要受化肥施用和畜禽養(yǎng)殖業(yè)產(chǎn)生的NH3-N、TN和COD污染;窯崗嘴—新安渡口段水質(zhì)受城鎮(zhèn)生活污水直排和污水廠尾水排放影響,NH3-N污染嚴(yán)重;汛期污染程度伴隨集中降雨強(qiáng)度逐步減弱。

    因此,為保障淠河流域水生態(tài)環(huán)境,建議從以下幾個(gè)方面加強(qiáng)措施:① 加強(qiáng)濕地保護(hù)修復(fù),增強(qiáng)流域自凈能力。針對(duì)窯崗嘴—新安渡口段氨氮高負(fù)荷問題,建議因地制宜,在充分依托淠河現(xiàn)有濕地資源的同時(shí),對(duì)受河道采砂等活動(dòng)破壞的濕地植被區(qū)進(jìn)行恢復(fù);在支流入干流處等關(guān)鍵區(qū)域建設(shè)生態(tài)濕地。利用水生植物修復(fù)污染水體。② 加強(qiáng)農(nóng)業(yè)面源污染防控。首先,應(yīng)科學(xué)確定高效施肥方式,推進(jìn)有機(jī)肥替代化肥,指導(dǎo)農(nóng)民科學(xué)合理、科學(xué)使用農(nóng)藥;其次,要連通農(nóng)田內(nèi)部溝、塘,優(yōu)化布設(shè)生態(tài)溝、生態(tài)塘和植被過濾帶,提高對(duì)徑流污染的攔截能力。③ 完善污水截污管網(wǎng)系統(tǒng)。加強(qiáng)對(duì)沿河排口、暗涵內(nèi)排口、沿河截留干管的監(jiān)管,嚴(yán)厲打擊偷排漏排;在重點(diǎn)入河排污口實(shí)現(xiàn)數(shù)字化自動(dòng)監(jiān)控,推進(jìn)雨污分流改造。

    參考文獻(xiàn):[1]? CARPENTER S R,CARACO N F,CORRELL D L,et al.Nonpoint pollution of surface waters with phosphorus and nitrogen[J].Ecological Applications,1998,8(3):559-568.

    [2]? WANG J Y,DA L J,SONG K,et al.Temporal variations of surface water quality in urban,suburban and rural areas during rapid urbanization in Shanghai,China[J].Environmental Pollution,2008,152(2):387-393.

    [3]? 王平,趙琦琳,張妍,等.撫仙湖水體磷的時(shí)空分布特征分析[J].人民長江,2023,54(增2):126-130.

    [4]? MOHEBBI M R,SAEEDI R,MONTAZERI A,et al.Assessment of water quality in groundwater resources of Iran using a modified drinking water quality index (DWQI)[J].Ecological Indicators,2013,30:28-34.

    [5]? WU Z S,LAI X J,LI K Y.Water quality assessment of rivers in Lake Chaohu Basin (China) using water quality index[J].Ecological Indicators,2021,121∶107021.

    [6]? 齊家蕙,謝崇寶,楊麗原.沂河水質(zhì)評(píng)價(jià)模型研究及其應(yīng)用[J].中國農(nóng)村水利水電,2023(8):103-110.

    [7]? 黃志,吳蕾,海子彬,等.干旱年份淠河水污染綜合治理研究[J].安徽農(nóng)學(xué)通報(bào),2021,27(21):135-137,168.

    [8]? 汪萬芬,譚綠貴,劉曉升.基于模糊綜合評(píng)價(jià)法的淠河水環(huán)境質(zhì)量評(píng)價(jià)[J].資源開發(fā)與市場,2008(5):410-412,416.

    [9]? 汪萬芬,戈海玉,劉曉升.六安市地表水環(huán)境質(zhì)量評(píng)價(jià)及變化趨勢分析[J].安徽建筑工業(yè)學(xué)院學(xué)報(bào)(自然科學(xué)版),2009,17:59-62.

    [10]李遠(yuǎn)平,楊太保,馬建國.淠河流域汛期降水集中度和集中期的變化特征[J].資源科學(xué),2012,34(3):418-423.

    [11]PESCE S F,WUNDERLIN D A.Use of water quality indices to verify the impact of Córdoba City (Argentina) on Suquía River[J].Water Research,2000,34(11):2915-2926.

    [12]KOER M A T,SEVGILI H.Parameters selection for water quality index in the assessment of the environmental impacts of land-based trout farms[J].Ecological Indicators,2014,36:672-681.

    [13]劉琦,田雨露,劉洋,等.兩類水質(zhì)綜合評(píng)價(jià)方法的特點(diǎn)及其在河流水環(huán)境管理中的作用[J].生態(tài)學(xué)報(bào),2019,39(20):7538-7546.

    [14]吳利,張雁,陳亞軍,等.巢湖流域河湖系統(tǒng)水環(huán)境因子分布特征及營養(yǎng)狀態(tài)評(píng)價(jià)[J].水生態(tài)學(xué)雜志,2023,44(1):65-72.

    [15]范澤宇,袁林,周亦,等.基于WQI法的滴水湖及入湖河道水質(zhì)研究[J].上海海洋大學(xué)學(xué)報(bào),2021,30(2):276-283.

    [16]JONNALAGADDA S B,MHERE G.Water quality of the odzi river in the eastern highlands of zimbabwe[J].Water Research,2001,35(10):2371-2376.

    [17]代孟均,張兵,杜倩倩,等.不同緩沖區(qū)的土地利用方式對(duì)地表水水質(zhì)的影響:以海河流域天津段為例[J].環(huán)境科學(xué),2024,45(3):1512-1524.

    [18]MLLER B,F(xiàn)INGER D,STURM M,et al.Present and past bio-available phosphorus budget in the ultra-oligotrophic Lake Brienz[J].Aquatic Sciences,2007,69:227-239.

    [19]張航,陳亞.2016~2020年長江下游干流水質(zhì)變化特征分析[J].水利水電快報(bào),2023,44(4):123-127,135.

    [20]吳一帆,唐洋博,李瑋,等.生態(tài)系統(tǒng)優(yōu)先保護(hù)區(qū)識(shí)別及生態(tài)治理成效分析:以六安主城區(qū)為例[J].人民長江,2023,54(5):112-120.

    [21]周欣.環(huán)巢湖濕地生態(tài)修復(fù)研究:以蘆溪濕地為例[J].治淮,2024(1):40-42.

    [22]宋大平,左強(qiáng),劉本生,等.農(nóng)業(yè)面源污染中氮排放時(shí)空變化及其健康風(fēng)險(xiǎn)評(píng)價(jià)研究:以淮河流域?yàn)槔跩].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2018,37(6):1219-1231.

    [23]馬騰飛,石林,吳超群,等.華南沿海飲用水水源水庫水質(zhì)變化及影響因素[J].人民長江,2021,52(10):46-54.

    (編輯:劉 媛)

    Spatial and temporal variation characteristics of water quality and its influencing factors in Pihe River

    ZHU Chenjing,TAN Qian,GUO Zhiyang,WANG Kexin,WU Shi

    (Anhui Survey & Design Institute of Water Resources & Hydropower Co.,Ltd.,Hefei 230088,China)

    Abstract:

    To comprehensively reveal water quality pollution and the spatiotemporal variation characteristics of the Pihe River Basin in recent years,the water quality index (WQI) method was adopted to evaluate the river water quality based on monitoring data of pH,dissolved oxygen (DO),chemical oxygen demand (COD),5-day biochemical oxygen demand (BOD5),ammonia nitrogen (NH3-N),total nitrogen (TN),and total phosphorus (TP) from the mainstream of the Pihe River from 2018 to 2022.Additionally,correlation analysis was considered to identify key water quality indicators influencing water quality.The results were as follows:① The mainstream of the Pihe River has generally favorable water quality,with WQI averaging at (75.94±7.47),ranking in the “good” range.However,during non-flood seasons,the concentration of NH3-N ranged from 0.020~1.480 mg/L,signifying a relatively high load.②The spatial distribution of water quality ranked as follows:Hengpaitou>Dadiangang>Yaogangzui>Xin′an Ferry.③ At the time level,the water quality of mainstream during the flood season (June to September) was markedly superior to that during the non-flood season (January to May and October to December).④ The correlation analysis showed that the water was mainly polluted by NH3-N,TN,and COD in the study area.In the Hengpaitou-Yaogangzui section and Xin′an Ferry-Dadiangang section with arable land as a major land use type,the water was mainly affected by COD,TN,and NH3-N generated from fertilizer application and livestock farming.In the Yaogangzui-Xin′an Ferry section,the water quality was severely impacted by the direct discharge of urban domestic sewage and effluent from sewage treatment plants,causing serious NH3-N pollution.During the flood season,water pollution gradually alleviates as concentrated rainfalls happen.From the perspective of protecting the ecological environment of the Pihe River,this study proposed to boost the self-purification capacity of the watershed,strengthen the prevention of agricultural non-point source pollution,and improve the intercepting and treating sewage systems.

    Key words:

    water quality; WQI; variation trend; influencing factors; Pihe River

    猜你喜歡
    淠河新安渡口
    為君守衛(wèi)
    界首渡口緬懷紅軍
    西江月(2021年3期)2021-12-21 06:34:18
    新安灘
    漢江渡口——兩代艄公的“價(jià)值觀”
    淠河水域淡水魚類資源調(diào)查與分析
    每一個(gè)彷徨的人生渡口,唯有自渡
    文苑(2020年6期)2020-06-22 08:42:04
    張新安教授簡介
    六安市淠河河道治理的城市生態(tài)效應(yīng)研究
    新安灘
    程?hào)|斌的詩
    詩歌月刊(2018年4期)2018-09-18 02:11:56
    亚洲欧洲日产国产| 天天一区二区日本电影三级| 女的被弄到高潮叫床怎么办| 99热这里只有是精品在线观看| 国产乱来视频区| 少妇人妻精品综合一区二区| 日本一二三区视频观看| 一本一本综合久久| 桃色一区二区三区在线观看| 美女高潮的动态| 亚洲人与动物交配视频| 国产视频首页在线观看| 综合色av麻豆| 亚洲自拍偷在线| 2022亚洲国产成人精品| 免费无遮挡裸体视频| 99久久精品一区二区三区| 3wmmmm亚洲av在线观看| 国产高清三级在线| 国产精品永久免费网站| 国产美女午夜福利| 免费看日本二区| 晚上一个人看的免费电影| 日韩亚洲欧美综合| videossex国产| 久久综合国产亚洲精品| 国产乱人视频| 99久久九九国产精品国产免费| 欧美激情久久久久久爽电影| 午夜视频国产福利| 中文字幕av在线有码专区| 国产久久久一区二区三区| 最近2019中文字幕mv第一页| 久久久久久国产a免费观看| 亚洲国产精品sss在线观看| 看黄色毛片网站| 亚洲美女搞黄在线观看| 久久人人爽人人片av| 成人性生交大片免费视频hd| 午夜福利视频1000在线观看| 一夜夜www| 哪个播放器可以免费观看大片| 欧美日韩国产亚洲二区| 午夜精品国产一区二区电影 | 美女xxoo啪啪120秒动态图| 偷拍熟女少妇极品色| 国产精品女同一区二区软件| 69人妻影院| 简卡轻食公司| 99热这里只有是精品在线观看| 欧美潮喷喷水| 国产精品综合久久久久久久免费| 亚州av有码| 免费观看的影片在线观看| 亚洲精品乱码久久久v下载方式| 中文资源天堂在线| 只有这里有精品99| 99九九线精品视频在线观看视频| 久久99热6这里只有精品| 国产欧美日韩精品一区二区| 日韩人妻高清精品专区| 狂野欧美激情性xxxx在线观看| 亚洲国产精品成人久久小说| 国产精品电影一区二区三区| 午夜福利网站1000一区二区三区| 国产一区二区亚洲精品在线观看| 少妇熟女欧美另类| 赤兔流量卡办理| 日韩一本色道免费dvd| 国产精品一区二区三区四区久久| 在线免费十八禁| 天堂网av新在线| 人妻夜夜爽99麻豆av| 日韩欧美三级三区| 最近2019中文字幕mv第一页| 久久久久免费精品人妻一区二区| 91av网一区二区| 亚洲av一区综合| 久久草成人影院| 中文在线观看免费www的网站| 91狼人影院| 夜夜看夜夜爽夜夜摸| 亚洲真实伦在线观看| 日本猛色少妇xxxxx猛交久久| 99久国产av精品国产电影| 国产在视频线在精品| 亚洲av二区三区四区| 3wmmmm亚洲av在线观看| 欧美xxxx黑人xx丫x性爽| 国产91av在线免费观看| 精品国产一区二区三区久久久樱花 | 国产精华一区二区三区| 校园人妻丝袜中文字幕| 久久韩国三级中文字幕| 毛片女人毛片| 亚洲av电影不卡..在线观看| 久久久a久久爽久久v久久| 中文字幕av在线有码专区| 久久久久久九九精品二区国产| 亚洲18禁久久av| 在线免费十八禁| kizo精华| 日韩欧美精品v在线| 一级二级三级毛片免费看| 一级黄片播放器| 久久久久免费精品人妻一区二区| 免费播放大片免费观看视频在线观看 | 人妻系列 视频| 亚洲在久久综合| 五月伊人婷婷丁香| 丰满少妇做爰视频| 国产成人a区在线观看| 国产一级毛片七仙女欲春2| 爱豆传媒免费全集在线观看| 观看美女的网站| 亚洲五月天丁香| 久久久久久久久久黄片| 亚洲av成人av| 最后的刺客免费高清国语| 麻豆av噜噜一区二区三区| 国产一级毛片在线| av专区在线播放| 精品国产露脸久久av麻豆 | 一级二级三级毛片免费看| 久久精品国产亚洲网站| 少妇熟女aⅴ在线视频| 精华霜和精华液先用哪个| 亚洲欧美日韩卡通动漫| 精品久久久久久成人av| 精品国产露脸久久av麻豆 | 欧美性感艳星| 老司机福利观看| 夫妻性生交免费视频一级片| 一边亲一边摸免费视频| 久久99热6这里只有精品| 国产精品99久久久久久久久| 色综合亚洲欧美另类图片| 国产在视频线精品| 国产激情偷乱视频一区二区| 国产 一区精品| 国产成人a∨麻豆精品| 天堂影院成人在线观看| 又粗又爽又猛毛片免费看| a级毛片免费高清观看在线播放| 51国产日韩欧美| 欧美+日韩+精品| 久久这里只有精品中国| 午夜福利在线观看吧| 级片在线观看| 在线观看av片永久免费下载| 亚洲av熟女| 国产精品久久久久久久久免| 晚上一个人看的免费电影| 免费电影在线观看免费观看| 精品久久久久久电影网 | 99热6这里只有精品| 精品人妻一区二区三区麻豆| 国产伦在线观看视频一区| 亚洲第一区二区三区不卡| 97人妻精品一区二区三区麻豆| 国产精品三级大全| 免费黄网站久久成人精品| 在线天堂最新版资源| 欧美精品国产亚洲| 精品99又大又爽又粗少妇毛片| 波野结衣二区三区在线| 有码 亚洲区| 蜜桃亚洲精品一区二区三区| av.在线天堂| 国产亚洲最大av| 精品久久久久久成人av| 欧美一区二区精品小视频在线| 亚洲最大成人av| 亚洲久久久久久中文字幕| 国产av码专区亚洲av| 99在线人妻在线中文字幕| 欧美xxxx性猛交bbbb| 国产成人一区二区在线| 久久久久久久久大av| 国产老妇伦熟女老妇高清| 午夜免费男女啪啪视频观看| 亚洲美女视频黄频| 中文字幕精品亚洲无线码一区| 欧美人与善性xxx| 高清日韩中文字幕在线| 成人av在线播放网站| 国产精品不卡视频一区二区| 3wmmmm亚洲av在线观看| 中文乱码字字幕精品一区二区三区 | 亚洲欧美成人精品一区二区| 3wmmmm亚洲av在线观看| 成人性生交大片免费视频hd| 黄色配什么色好看| a级毛片免费高清观看在线播放| 嘟嘟电影网在线观看| 免费av不卡在线播放| 国产男人的电影天堂91| 男人狂女人下面高潮的视频| 麻豆久久精品国产亚洲av| 中文字幕熟女人妻在线| 国产亚洲av片在线观看秒播厂 | 亚洲精品日韩av片在线观看| 午夜福利在线观看吧| 久久精品熟女亚洲av麻豆精品 | www.色视频.com| 久久久色成人| 久久鲁丝午夜福利片| 丰满人妻一区二区三区视频av| 老司机影院成人| 精品少妇黑人巨大在线播放 | 99久久精品一区二区三区| 国产午夜精品久久久久久一区二区三区| 亚洲五月天丁香| 99热精品在线国产| 日韩av不卡免费在线播放| 亚洲电影在线观看av| 亚洲欧美一区二区三区国产| 麻豆av噜噜一区二区三区| 国产日韩欧美在线精品| 网址你懂的国产日韩在线| 国产精品福利在线免费观看| 长腿黑丝高跟| 午夜福利成人在线免费观看| av在线蜜桃| 欧美zozozo另类| 国产欧美日韩精品一区二区| 99久久精品热视频| 极品教师在线视频| 国产精品国产三级专区第一集| 观看美女的网站| 啦啦啦观看免费观看视频高清| 免费看a级黄色片| 久久久久九九精品影院| 特级一级黄色大片| 日韩强制内射视频| av专区在线播放| www.色视频.com| 国产亚洲一区二区精品| 国内精品美女久久久久久| 青青草视频在线视频观看| 在线播放国产精品三级| 2021天堂中文幕一二区在线观| 好男人视频免费观看在线| 成人欧美大片| 久久久成人免费电影| 午夜爱爱视频在线播放| 在线观看66精品国产| 亚洲国产欧美人成| 大又大粗又爽又黄少妇毛片口| 亚洲欧美日韩无卡精品| 久99久视频精品免费| 国产毛片a区久久久久| 亚洲国产欧洲综合997久久,| 国产精品久久电影中文字幕| 国产三级在线视频| 午夜激情福利司机影院| 你懂的网址亚洲精品在线观看 | 中文字幕久久专区| 日韩一区二区视频免费看| 3wmmmm亚洲av在线观看| 少妇熟女欧美另类| 中文字幕av成人在线电影| 毛片一级片免费看久久久久| 在线播放无遮挡| 国产精品国产三级国产av玫瑰| 三级国产精品片| 亚洲第一区二区三区不卡| 天美传媒精品一区二区| 一级毛片aaaaaa免费看小| 日韩精品有码人妻一区| 国产乱来视频区| 人妻夜夜爽99麻豆av| 国产高清有码在线观看视频| 亚洲人成网站在线观看播放| 亚洲一级一片aⅴ在线观看| 久久久久性生活片| 午夜福利高清视频| 1024手机看黄色片| 日本-黄色视频高清免费观看| 亚洲av一区综合| 99久久精品一区二区三区| 亚洲怡红院男人天堂| 一级黄片播放器| 国产免费福利视频在线观看| 69人妻影院| 97超碰精品成人国产| 最近最新中文字幕免费大全7| 亚洲一区高清亚洲精品| 在线观看av片永久免费下载| 免费观看的影片在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 一级毛片久久久久久久久女| 国产大屁股一区二区在线视频| 寂寞人妻少妇视频99o| 亚洲精品乱久久久久久| 九九在线视频观看精品| 免费观看a级毛片全部| 日韩av在线大香蕉| 国产大屁股一区二区在线视频| 久久精品人妻少妇| 天天躁日日操中文字幕| 在线播放国产精品三级| 亚洲中文字幕日韩| 亚洲av免费在线观看| 国产麻豆成人av免费视频| 亚洲成av人片在线播放无| 亚洲精品日韩在线中文字幕| 国产精品爽爽va在线观看网站| 99热6这里只有精品| 亚洲,欧美,日韩| 菩萨蛮人人尽说江南好唐韦庄 | 欧美另类亚洲清纯唯美| 日本黄色视频三级网站网址| 国产精品人妻久久久影院| 如何舔出高潮| av播播在线观看一区| 在线免费十八禁| 国产午夜福利久久久久久| 色播亚洲综合网| 伦精品一区二区三区| АⅤ资源中文在线天堂| 久久久久久久亚洲中文字幕| 精品久久久久久电影网 | 国产成年人精品一区二区| 午夜福利成人在线免费观看| 2022亚洲国产成人精品| 欧美xxxx性猛交bbbb| 国产午夜福利久久久久久| 免费搜索国产男女视频| 搞女人的毛片| 一级毛片久久久久久久久女| 日韩在线高清观看一区二区三区| 精品久久国产蜜桃| 天堂av国产一区二区熟女人妻| 国产成人免费观看mmmm| 精品久久久久久成人av| 精华霜和精华液先用哪个| 久久人人爽人人片av| 国产精品一区二区性色av| 欧美成人一区二区免费高清观看| 看免费成人av毛片| 一二三四中文在线观看免费高清| 色5月婷婷丁香| 中文字幕熟女人妻在线| 午夜免费男女啪啪视频观看| 有码 亚洲区| 国产淫语在线视频| 精品久久国产蜜桃| 日韩欧美在线乱码| 乱人视频在线观看| 内射极品少妇av片p| 欧美性猛交黑人性爽| 亚洲av成人精品一区久久| 久久这里只有精品中国| 亚洲久久久久久中文字幕| 久久99蜜桃精品久久| 日韩 亚洲 欧美在线| 亚洲电影在线观看av| 嫩草影院入口| 在线观看66精品国产| av视频在线观看入口| 国产成人一区二区在线| 一级二级三级毛片免费看| 天堂av国产一区二区熟女人妻| 国产欧美另类精品又又久久亚洲欧美| 免费观看精品视频网站| 国产男人的电影天堂91| 久久久精品大字幕| 亚洲内射少妇av| 可以在线观看毛片的网站| 亚洲自偷自拍三级| 国产精品野战在线观看| 有码 亚洲区| 国产激情偷乱视频一区二区| 亚洲国产欧美人成| a级一级毛片免费在线观看| 爱豆传媒免费全集在线观看| 国产一区有黄有色的免费视频 | 成年女人永久免费观看视频| 我要搜黄色片| 麻豆av噜噜一区二区三区| 成年版毛片免费区| 国产一区二区亚洲精品在线观看| 久久草成人影院| 纵有疾风起免费观看全集完整版 | 久久久久精品久久久久真实原创| 日本一二三区视频观看| 色5月婷婷丁香| АⅤ资源中文在线天堂| 亚洲欧美精品综合久久99| 午夜福利在线观看免费完整高清在| 久久久精品欧美日韩精品| 99久久成人亚洲精品观看| 亚洲精品亚洲一区二区| 一夜夜www| 国产激情偷乱视频一区二区| 看非洲黑人一级黄片| 国产老妇伦熟女老妇高清| 国产精品乱码一区二三区的特点| 亚洲国产欧美人成| 日日摸夜夜添夜夜添av毛片| 成人漫画全彩无遮挡| 免费观看精品视频网站| 麻豆国产97在线/欧美| 青春草视频在线免费观看| 国产精品人妻久久久影院| 男女那种视频在线观看| 亚洲精品色激情综合| 又黄又爽又刺激的免费视频.| 高清日韩中文字幕在线| 又黄又爽又刺激的免费视频.| 免费电影在线观看免费观看| 男女那种视频在线观看| 又粗又硬又长又爽又黄的视频| 观看免费一级毛片| 国产精品三级大全| 特级一级黄色大片| 一区二区三区免费毛片| 国产精品久久久久久久电影| 日本欧美国产在线视频| av在线蜜桃| 亚洲人与动物交配视频| 国产av在哪里看| 如何舔出高潮| 美女国产视频在线观看| 91精品国产九色| 岛国在线免费视频观看| 亚洲va在线va天堂va国产| 午夜视频国产福利| 国产精品一区二区三区四区久久| 特大巨黑吊av在线直播| 欧美一区二区精品小视频在线| 最后的刺客免费高清国语| 网址你懂的国产日韩在线| 最近中文字幕高清免费大全6| 国产亚洲精品av在线| 国产一区二区亚洲精品在线观看| 亚洲av成人精品一区久久| 亚洲av成人精品一二三区| 桃色一区二区三区在线观看| 全区人妻精品视频| 视频中文字幕在线观看| 亚洲av成人精品一区久久| 少妇人妻一区二区三区视频| 国产精品嫩草影院av在线观看| 成人高潮视频无遮挡免费网站| 有码 亚洲区| 国产成人精品一,二区| 精品免费久久久久久久清纯| 精品一区二区三区视频在线| 赤兔流量卡办理| 欧美+日韩+精品| 日本午夜av视频| 啦啦啦啦在线视频资源| av专区在线播放| 国产免费又黄又爽又色| 国产不卡一卡二| 中文字幕亚洲精品专区| 91aial.com中文字幕在线观看| 久久久国产成人免费| 亚洲成色77777| 亚洲av熟女| 国产免费男女视频| 激情 狠狠 欧美| 国产精品不卡视频一区二区| 国产男人的电影天堂91| 欧美成人午夜免费资源| 三级国产精品欧美在线观看| 深爱激情五月婷婷| 久久久久久久久中文| 亚洲国产精品专区欧美| 欧美xxxx性猛交bbbb| 国产亚洲精品av在线| 久久婷婷人人爽人人干人人爱| 久久久国产成人免费| 免费人成在线观看视频色| 久久精品影院6| 日本一本二区三区精品| 亚洲av男天堂| 一级爰片在线观看| 亚洲av成人av| 日韩一区二区视频免费看| 国产国拍精品亚洲av在线观看| 别揉我奶头 嗯啊视频| 国产极品天堂在线| 一级黄片播放器| 中文字幕av成人在线电影| 精品人妻熟女av久视频| 永久免费av网站大全| 国内精品宾馆在线| 麻豆久久精品国产亚洲av| 1024手机看黄色片| 两个人的视频大全免费| 少妇猛男粗大的猛烈进出视频 | 国产成人a∨麻豆精品| 亚洲成av人片在线播放无| 中文精品一卡2卡3卡4更新| 国产综合懂色| 国产私拍福利视频在线观看| 国产91av在线免费观看| 自拍偷自拍亚洲精品老妇| 久久精品国产亚洲网站| 久久精品人妻少妇| 少妇熟女欧美另类| 午夜福利在线在线| 麻豆精品久久久久久蜜桃| 日韩一本色道免费dvd| 一个人观看的视频www高清免费观看| 岛国毛片在线播放| 亚洲国产精品成人综合色| 久久久久国产网址| 亚洲婷婷狠狠爱综合网| 久久鲁丝午夜福利片| 一级毛片久久久久久久久女| 国产伦理片在线播放av一区| 我要看日韩黄色一级片| 中文资源天堂在线| 91在线精品国自产拍蜜月| 国产精品国产高清国产av| 99久久无色码亚洲精品果冻| 欧美激情久久久久久爽电影| 欧美一区二区亚洲| 日日摸夜夜添夜夜添av毛片| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 国产一区二区在线观看日韩| 国产精品蜜桃在线观看| 少妇人妻一区二区三区视频| 久久人人爽人人片av| 精品人妻偷拍中文字幕| 国产老妇伦熟女老妇高清| 中文天堂在线官网| 久久精品国产亚洲网站| 午夜日本视频在线| 久久人人爽人人片av| 久久99精品国语久久久| 麻豆av噜噜一区二区三区| 村上凉子中文字幕在线| 韩国av在线不卡| 国产精品国产高清国产av| 国产欧美另类精品又又久久亚洲欧美| 乱码一卡2卡4卡精品| 午夜a级毛片| 免费黄网站久久成人精品| 久久99精品国语久久久| 简卡轻食公司| 日韩av在线免费看完整版不卡| 国产麻豆成人av免费视频| 一区二区三区高清视频在线| 麻豆成人午夜福利视频| 身体一侧抽搐| 一卡2卡三卡四卡精品乱码亚洲| av免费在线看不卡| 中文资源天堂在线| 亚洲av免费高清在线观看| 日本五十路高清| 天堂影院成人在线观看| 狂野欧美激情性xxxx在线观看| 成年版毛片免费区| 免费黄色在线免费观看| 亚洲欧美成人综合另类久久久 | 亚洲av电影在线观看一区二区三区 | 国产乱来视频区| 日韩大片免费观看网站 | 久久精品综合一区二区三区| 欧美激情久久久久久爽电影| 国产精品嫩草影院av在线观看| 亚洲欧美精品专区久久| 国产av一区在线观看免费| 一级毛片电影观看 | 免费观看性生交大片5| 免费观看在线日韩| 亚洲国产成人一精品久久久| 久久人人爽人人爽人人片va| 日本熟妇午夜| 亚洲成人久久爱视频| 最新中文字幕久久久久| 啦啦啦韩国在线观看视频| 久久精品久久精品一区二区三区| 亚洲欧美精品自产自拍| 韩国高清视频一区二区三区| 高清日韩中文字幕在线| 日韩成人伦理影院| 欧美色视频一区免费| av国产免费在线观看| 中文字幕久久专区| 欧美成人a在线观看| 久久精品人妻少妇| 国产综合懂色| 两个人视频免费观看高清| 插逼视频在线观看| 欧美成人精品欧美一级黄| 长腿黑丝高跟| 免费av不卡在线播放| 久久人人爽人人片av| 亚洲av男天堂| 亚洲最大成人中文| 水蜜桃什么品种好| 七月丁香在线播放| 久久鲁丝午夜福利片| 亚洲国产色片| 国产精品av视频在线免费观看| 午夜视频国产福利| 免费观看的影片在线观看| 少妇人妻一区二区三区视频| 久久久精品大字幕| 六月丁香七月| 欧美日韩在线观看h| 国产色婷婷99| 一个人观看的视频www高清免费观看| 天天躁夜夜躁狠狠久久av| 日本三级黄在线观看| 国产高清国产精品国产三级 | 99热6这里只有精品| 亚洲av男天堂|