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

    從水質(zhì)變化趨勢看常州市武南區(qū)域重點(diǎn)河流治理成效

    2019-10-23 09:20:48康若熙謝文理趙文星龐曉斯
    中國環(huán)境科學(xué) 2019年10期
    關(guān)鍵詞:內(nèi)梅羅高錳酸鹽永安

    康若熙,謝文理,趙文星,龐曉斯,余 剛*

    從水質(zhì)變化趨勢看常州市武南區(qū)域重點(diǎn)河流治理成效

    康若熙1,謝文理2,趙文星1,龐曉斯1,余 剛1*

    (1.清華大學(xué)環(huán)境學(xué)院,新興有機(jī)污染物控制北京市重點(diǎn)實(shí)驗(yàn)室,北京 100084;2.常州市環(huán)境監(jiān)測中心,江蘇 常州 213001)

    為探究武南區(qū)域重點(diǎn)河流水質(zhì)的變化規(guī)律、驅(qū)動(dòng)因素以及河流治理成效,基于2006~2018年連續(xù)水質(zhì)監(jiān)測數(shù)據(jù),綜合分析了4條重點(diǎn)河流(太滆運(yùn)河、武宜運(yùn)河、武進(jìn)港和永安河)水質(zhì)演變趨勢,并對(duì)污染較重的永安河各項(xiàng)水質(zhì)指標(biāo)進(jìn)行了季節(jié)性分析和相關(guān)性分析.結(jié)果表明:2006~2018年,4條河流水質(zhì)整體呈好轉(zhuǎn)趨勢,修正內(nèi)梅羅指數(shù)分別下降36.2%,31.5%,56.4%,48.7%,受河流清淤工程影響,永安河2017年水質(zhì)有所下降;4條河流氨氮濃度、總氮濃度與高錳酸鹽指數(shù)下降趨勢明顯(<0.05),總磷濃度則存在一定波動(dòng);永安河的總氮、氨氮和高錳酸鹽指數(shù)間或存在同源關(guān)系,氨氮和總氮季節(jié)性變化明顯,雨季濃度低于旱季,總磷和高錳酸鹽指數(shù)沒有明顯季節(jié)性變化趨勢;城鎮(zhèn)化發(fā)展與產(chǎn)業(yè)結(jié)構(gòu)由傳統(tǒng)工業(yè)、農(nóng)業(yè)向第三產(chǎn)業(yè)的轉(zhuǎn)變均對(duì)區(qū)域水環(huán)境改善有積極作用.

    武南區(qū)域;水質(zhì);變化趨勢;驅(qū)動(dòng)因素

    太湖是我國第三大淡水湖泊,地處長江下游的太湖流域,是我國人口最密集、工農(nóng)業(yè)最發(fā)達(dá)的地區(qū)之一[1].然而,隨著經(jīng)濟(jì)的高速發(fā)展,流域內(nèi)水質(zhì)下降、水體富營養(yǎng)化加劇的問題日益突出[2].相關(guān)研究表明,近20年來,太湖流域水體氮、磷營養(yǎng)元素污染嚴(yán)重,湖體夏季水華藍(lán)藻暴發(fā)等情況已嚴(yán)重影響到流域水環(huán)境質(zhì)量及生活用水安全[3-5],因入湖河流污染負(fù)荷重,流速低,自凈能力差所造成的流域生態(tài)風(fēng)險(xiǎn)上升、水質(zhì)性缺水情況增多等問題層出不窮[6-7],流域水環(huán)境污染與水生態(tài)受到了廣泛關(guān)注.

    已有部分學(xué)者分析研究了太湖流域水質(zhì)變化趨勢.于東升[8]等采用Mann-Kendall法對(duì)太湖西岸水質(zhì)變化趨勢進(jìn)行了分析,結(jié)果表明,2009~2015年太湖西岸水質(zhì)呈逐年改善趨勢,影響區(qū)域水質(zhì)的主要驅(qū)動(dòng)因子為氨氮和COD;左一鳴等[9]對(duì)太湖各水質(zhì)指標(biāo)進(jìn)行了聚類分析,表明總氮、氨氮與亞硝酸鹽氮,高錳酸鹽與總磷具有相同的變化趨勢;朱偉等[10]對(duì)太湖梅梁灣和貢湖灣的水質(zhì)變化趨勢分析情況顯示,2010年以來,監(jiān)測區(qū)域總氮濃度呈波動(dòng)下降趨勢,總磷濃度在2015、2016年有所反彈.然而,已有研究多側(cè)重于水質(zhì)指標(biāo)變化情況分析,對(duì)影響水質(zhì)變化的深層次因素探討較少,且研究對(duì)象多為太湖湖體或流域內(nèi)單條河流,尚缺乏針對(duì)太湖流域西北部區(qū)域長時(shí)間跨度水質(zhì)變化趨勢及影響因素的深入分析.

    本研究以太湖流域西北部武南區(qū)域的4條重點(diǎn)河流為對(duì)象,系統(tǒng)分析了河流上、中、下游點(diǎn)位在2006~2018年間的水質(zhì)監(jiān)測數(shù)據(jù),采用修正內(nèi)梅羅指數(shù)法探討了各河流水質(zhì)的總體情況,并對(duì)細(xì)分水質(zhì)指標(biāo)的變化趨勢進(jìn)行了Mann-Kendall檢驗(yàn),同時(shí)通過主成分分析與相關(guān)性分析探討了水質(zhì)變化的內(nèi)在驅(qū)動(dòng)因素.此外,本研究對(duì)污染較重河流永安河的各水質(zhì)指標(biāo)進(jìn)行了季節(jié)性和相關(guān)性分析,研究結(jié)果對(duì)科學(xué)評(píng)價(jià)武南區(qū)域水環(huán)境治理進(jìn)展和完善治理方案具有重要意義.

    1 材料與方法

    1.1 研究區(qū)域概況

    常州市武進(jìn)區(qū)地處江蘇省南部,太湖流域西北部,東鄰無錫市,西毗金壇區(qū)和丹陽市,南接宜興市,北靠常州天寧區(qū)、鐘樓區(qū)和新北區(qū).區(qū)域內(nèi)地勢低平,河網(wǎng)繁復(fù)稠密,屬于典型平原河網(wǎng)結(jié)構(gòu).

    武南區(qū)域是指武進(jìn)區(qū)內(nèi)京杭運(yùn)河以南的部分,面積約800km2,是多條河流的入湖處.武南區(qū)域內(nèi)“兩縱(武進(jìn)港、武宜運(yùn)河)一橫(太滆運(yùn)河)”是骨干河流,接納了區(qū)域內(nèi)主要污染負(fù)荷,其中,武進(jìn)港與太滆運(yùn)河為太湖西北部主要入湖河流,武宜運(yùn)河為太湖湖西區(qū)主要入湖河流[11].此外,永安河是區(qū)域內(nèi)污染較重的一條河流.該區(qū)域由于環(huán)境容量有限、本地和外來污染負(fù)荷高,多個(gè)監(jiān)測斷面水質(zhì)超標(biāo)時(shí)常發(fā)生.

    1.2 點(diǎn)位選擇與數(shù)據(jù)來源

    武南區(qū)域內(nèi)共包含5個(gè)國考斷面和9個(gè)省控點(diǎn)位,覆蓋湖區(qū)及8條入太湖河流.因本研究主要針對(duì)區(qū)域內(nèi)骨干河流(武進(jìn)港、武宜運(yùn)河、太滆運(yùn)河)與污染較重河流(永安河),考慮到各點(diǎn)位已有監(jiān)測數(shù)據(jù)間的充分性與可比性,選取4條河流上、中、下游共12個(gè)國考、省考和重點(diǎn)監(jiān)測點(diǎn)位,點(diǎn)位名稱及布設(shè)位置見圖1.

    本研究數(shù)據(jù)為2006~2018年常州市環(huán)境監(jiān)測站水質(zhì)例行監(jiān)測結(jié)果,各點(diǎn)位每月均有至少1次水質(zhì)監(jiān)測數(shù)據(jù),指標(biāo)包括高錳酸鹽指數(shù)、氨氮、總磷、流速、溫度等.武進(jìn)區(qū)城鎮(zhèn)化與產(chǎn)業(yè)結(jié)構(gòu)數(shù)據(jù)來源為2007~2018年常州統(tǒng)計(jì)年鑒[12].

    圖1 武南區(qū)域重點(diǎn)河流研究點(diǎn)位分布圖

    1.3 評(píng)價(jià)方法

    1.3.1 修正內(nèi)梅羅指數(shù) 常用的水質(zhì)評(píng)價(jià)方法包括單因子指數(shù)評(píng)價(jià)法[13]、綜合污染指數(shù)評(píng)價(jià)法[14]、模糊綜合評(píng)價(jià)法[15]、層次分析法[16]等.其中,綜合污染指數(shù)法具有操作簡便,可直觀反映多個(gè)水質(zhì)參數(shù)與相應(yīng)標(biāo)準(zhǔn)間的綜合對(duì)應(yīng)關(guān)系[17].內(nèi)梅羅指數(shù)是一種兼顧極值或突出最大值的計(jì)權(quán)型多因子環(huán)境質(zhì)量指數(shù),其相較于SWQI法與均值法更加科學(xué)合理[18].修正的內(nèi)梅羅法相較傳統(tǒng)內(nèi)梅羅法更多地考慮了各污染指標(biāo)間的權(quán)重,修正了極端值對(duì)指數(shù)的影響,其主要計(jì)算公式如下[19]:

    本研究采用修正后的內(nèi)梅羅指數(shù)對(duì)研究區(qū)域水質(zhì)進(jìn)行衡量與分析.

    1.3.2 Mann-Kendall統(tǒng)計(jì)檢驗(yàn) Mann-Kendall統(tǒng)計(jì)檢驗(yàn)是環(huán)境數(shù)據(jù)時(shí)間序列趨勢分析的一種方法,也是檢驗(yàn)水文水質(zhì)指標(biāo)時(shí)間序列單調(diào)趨勢的有效工具[20].該方法能有效區(qū)分某一自然過程是處于自然波動(dòng)還是有確定的變化趨勢[21].本研究所涉及的一系列水質(zhì)指標(biāo)為時(shí)間變量,故用Mann-Kendall統(tǒng)計(jì)檢驗(yàn)對(duì)各水質(zhì)指標(biāo)進(jìn)行趨勢分析.

    1.3.3 主成分分析 主成分分析法是一種將多維因子納入同一系統(tǒng)進(jìn)行降維定量化研究的多元統(tǒng)計(jì)分析方法,它可從諸要素中找出影響水質(zhì)狀況的主導(dǎo)因子,逐級(jí)劃分合并以達(dá)到綜合評(píng)價(jià)水環(huán)境的目的[22].本研究采用主成分分析法與相關(guān)性分析法歸納分析各水質(zhì)指標(biāo)的深層次影響因素.

    2 結(jié)果與討論

    2.1 修正內(nèi)梅羅指數(shù)分析結(jié)果

    修正內(nèi)梅羅指數(shù)的評(píng)價(jià)因子選擇主要從水資源利用角度出發(fā),評(píng)價(jià)結(jié)果盡可能客觀反映區(qū)域內(nèi)的水資源用途與污染類型[23].武南片區(qū)污染源主要由工業(yè)源、農(nóng)業(yè)源和生活源三部分構(gòu)成,其中,生活源與畜禽養(yǎng)殖源排放造成的氮、磷污染倍受關(guān)注.綜合考慮研究區(qū)域污染特征,選取高錳酸鹽指數(shù)、氨氮、總氮、總磷4項(xiàng)具有代表性評(píng)價(jià)因子進(jìn)行水質(zhì)評(píng)價(jià).因修正內(nèi)梅羅指數(shù)法需在同一標(biāo)準(zhǔn)下比較,根據(jù)各河流點(diǎn)位已有的水質(zhì)標(biāo)準(zhǔn)值,選取《國家地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)GB3838-2002》[24]中IV類水質(zhì)標(biāo)準(zhǔn)作為基準(zhǔn).以河流上、中、下游水質(zhì)指標(biāo)年均值代表河流整體水質(zhì)狀況,所得修正內(nèi)梅羅指數(shù)計(jì)算結(jié)果見圖2.

    根據(jù)修正內(nèi)梅羅指數(shù)計(jì)算結(jié)果,2006~2018年,4條河流總體水質(zhì)呈好轉(zhuǎn)趨勢.太滆運(yùn)河整體呈上游水質(zhì)較好,中游污染較嚴(yán)重,下游水質(zhì)有所好轉(zhuǎn)的特點(diǎn),中游運(yùn)村點(diǎn)位沿岸分布的化工廠、機(jī)械廠、水泥廠、農(nóng)業(yè)面源以及錫溧漕河來水為河流主要污染來源[25-26].2006~2010年間,得益于“十一五”期間所采取的污染源監(jiān)管、河道污染攔截、上游湖蕩凈化等措施,太滆運(yùn)河水質(zhì)改善明顯,修正內(nèi)梅羅指數(shù)由1.64下降至0.94;2010~2018年,太滆運(yùn)河河口生態(tài)攔截工程對(duì)上游紅湖大橋水質(zhì)改善明顯,中游水污染控制工程持續(xù)推進(jìn),全河修正內(nèi)梅羅指數(shù)維持在0.90左右,水質(zhì)狀況較好.武進(jìn)港是入太湖梅梁灣的主要河流,受京杭運(yùn)河來水及沿岸污染源排放影響[27],上游慈瀆橋點(diǎn)位水質(zhì)污染較中、下游嚴(yán)重.“十一五”期間,武進(jìn)港治理工程對(duì)入太湖污染負(fù)荷量削減達(dá)20%,河流水質(zhì)提升,2006~2010年修正內(nèi)梅羅指數(shù)下降21.1%.2013~2016年,受沿線污染源排放和支流匯水的影響,武進(jìn)港水質(zhì)略有下降,但隨著河流整治工程的竣工,2017~2018年,武進(jìn)港各點(diǎn)位水質(zhì)有所好轉(zhuǎn),修正內(nèi)梅羅指數(shù)達(dá)0.95.武宜運(yùn)河流經(jīng)武南區(qū)域工業(yè)、農(nóng)業(yè)產(chǎn)業(yè)聚集地,承接區(qū)域內(nèi)較多工業(yè)、農(nóng)業(yè)、生活污染負(fù)荷,河流中、下游水質(zhì)狀況較差.受河流綜合治理工程影響,2006~2018年,武宜運(yùn)河水質(zhì)持續(xù)好轉(zhuǎn),修正內(nèi)梅羅指數(shù)下降56.4%.永安河流量僅為1.1m3/s,不及另三條河流的十分之一,河流流經(jīng)前黃鎮(zhèn)內(nèi)農(nóng)業(yè)區(qū)、工業(yè)區(qū)與集鎮(zhèn)生活區(qū),沿線農(nóng)業(yè)面源、工業(yè)源與生活源眾多,小流量與高沿線污染排放共同導(dǎo)致了河流水體污染嚴(yán)重.盡管2006~ 2018年,永安河存在明顯好轉(zhuǎn)趨勢,修正內(nèi)梅羅指數(shù)降幅達(dá)48.7%,但與其他3條河流仍存在較大差距.因當(dāng)?shù)卣掷m(xù)推動(dòng)“三河工程”,永安河拓浚工程于2017年10月底竣工,工程投資16億元,全長約20km,涉及項(xiàng)目包括河道清淤、河道拓寬、建設(shè)生態(tài)保護(hù)岸等.因河道清淤過程對(duì)河流水質(zhì)的不利影響,2017年永安河水質(zhì)有所下降,但清淤工程完成后,2018年河流水質(zhì)狀況明顯好轉(zhuǎn).

    圖2 2006~2018年4條重點(diǎn)河流修正內(nèi)梅羅指數(shù)

    2.2 4項(xiàng)水質(zhì)指標(biāo)時(shí)間跨度分析

    圖3是4條重點(diǎn)河流2006~2018年水質(zhì)指標(biāo)變化曲線,可見,太滆運(yùn)河、武宜運(yùn)河和武進(jìn)港的氨氮濃度、總氮濃度和高錳酸鹽指數(shù)均波動(dòng)下降,得益于“十一五”、“十二五”期間各項(xiàng)整治措施,2006~2010年3項(xiàng)水質(zhì)指標(biāo)濃度下降尤為明顯,2010~2015年各指標(biāo)維持低位波動(dòng).各河流高錳酸鹽指數(shù)基本達(dá)標(biāo),氨氮與總氮變化趨勢相似,武進(jìn)港各項(xiàng)水質(zhì)指標(biāo)相對(duì)較好.目前,區(qū)域內(nèi)3條主干河流氨氮濃度基本達(dá)到IV類水質(zhì)標(biāo)準(zhǔn).永安河氨氮污染仍較重,盡管自2006年至今,河流氨氮濃度得到了一定控制,但河流水質(zhì)仍處于劣V類水平,氨氮為永安河最應(yīng)引起重視的污染物.太滆運(yùn)河、武宜運(yùn)河和武進(jìn)港總磷濃度均符合IV類水質(zhì)標(biāo)準(zhǔn),且變化較平穩(wěn),永安河總磷濃度在2006~2017年均維持在相對(duì)較高的劣V類水平, 為重點(diǎn)關(guān)注的污染指標(biāo).自永安河拓浚工程施工完成后,河流水體磷污染得到了有效控制,2018年全年平均總磷濃度達(dá)到IV類水質(zhì)標(biāo)準(zhǔn).

    表1 Mann-Kendall水質(zhì)變化趨勢檢驗(yàn)

    注:*,**和***分別表示通過90%,95%和99%的顯著性檢驗(yàn).

    表1是Mann-Kendall法水質(zhì)變化趨勢檢驗(yàn)結(jié)果,可見4條河流氨氮濃度、總氮濃度和高錳酸鹽指數(shù)在2006~2018年下降趨勢顯著(<0.05).武進(jìn)港和武宜運(yùn)河總磷濃度下降趨勢也較明顯,而永安河和太滆運(yùn)河總磷濃度增減趨勢不明顯,存在不規(guī)則波動(dòng).

    2.3 水質(zhì)變化驅(qū)動(dòng)因素分析

    “十一五”至“十三五”期間,武進(jìn)區(qū)政府采取了系列措施治理當(dāng)?shù)厮h(huán)境,本研究重點(diǎn)分析主要措施對(duì)武南區(qū)域水質(zhì)變化的影響.在2006~2017年期間,城鎮(zhèn)化和產(chǎn)業(yè)結(jié)構(gòu)調(diào)整是推動(dòng)武進(jìn)區(qū)發(fā)展的兩大主要因素.武進(jìn)區(qū)政府相繼采取了系列產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化與經(jīng)濟(jì)結(jié)構(gòu)調(diào)整措施,如鼓勵(lì)產(chǎn)業(yè)集聚,培育新產(chǎn)業(yè),扶持服務(wù)業(yè)等[28],并通過民營經(jīng)濟(jì)帶動(dòng)、特色園區(qū)建設(shè)、創(chuàng)新創(chuàng)業(yè)發(fā)展、產(chǎn)教融合互促等做法實(shí)現(xiàn)了區(qū)域經(jīng)濟(jì)的高質(zhì)量發(fā)展[29].

    本研究采用SPSS 22統(tǒng)計(jì)軟件,對(duì)城鎮(zhèn)化推進(jìn)和產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化相關(guān)指標(biāo)分別進(jìn)行主成分分析.城鎮(zhèn)化發(fā)展可體現(xiàn)在經(jīng)濟(jì)與收入增長、人口與城市規(guī)模擴(kuò)大等方面[30],基于上述考慮,選取總產(chǎn)值、城鎮(zhèn)人口占比、人均綠地面積、人均道路面積和人均GDP共5項(xiàng)指標(biāo)作為城鎮(zhèn)化推進(jìn)的主要驅(qū)動(dòng)因素;區(qū)域產(chǎn)業(yè)結(jié)構(gòu)可由各產(chǎn)業(yè)產(chǎn)值以及其他特征產(chǎn)業(yè)數(shù)據(jù)體現(xiàn),選取農(nóng)作物播種面積、淡水養(yǎng)殖面積、工業(yè)產(chǎn)值占比、第三產(chǎn)業(yè)產(chǎn)值占比和農(nóng)林牧漁業(yè)產(chǎn)值占比5項(xiàng)指標(biāo)作為產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化的主要驅(qū)動(dòng)因素.主成分分析結(jié)果見表2,KMO值分別為0.719和0.679,Bartlett的球形檢驗(yàn)值均為0.000(<0.05),適合進(jìn)行因子分析.

    城鎮(zhèn)化推進(jìn)因素可提取1個(gè)有效因子,對(duì)方差的貢獻(xiàn)率為92.747%,基本涵蓋了全部信息.其中,總產(chǎn)值、城鎮(zhèn)人口占比、人均綠地面積、人均道路面積和人均GDP均占有較高正向權(quán)重,即地區(qū)經(jīng)濟(jì)發(fā)展越快,城市建設(shè)越好,該因子值越大,將該有效因子定義為城鎮(zhèn)化因子.產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化因素可提取出1個(gè)有效因子,因子對(duì)方差貢獻(xiàn)率為90.643%,其中,農(nóng)作物播種面積、淡水養(yǎng)殖面積、工業(yè)產(chǎn)值占比、和農(nóng)林牧漁業(yè)產(chǎn)值占比的載荷值均較高且為正,第三產(chǎn)業(yè)產(chǎn)值占比比值較高,但載荷值為負(fù),即區(qū)域產(chǎn)業(yè)越由工業(yè)、農(nóng)業(yè)向第三產(chǎn)業(yè)轉(zhuǎn)型,該因子值越小,將此因子命名為產(chǎn)業(yè)結(jié)構(gòu)因子.

    分別計(jì)算了2006~2017年城鎮(zhèn)化因子和產(chǎn)業(yè)結(jié)構(gòu)因子,并將兩因子分別與4條河流監(jiān)測點(diǎn)位年均綜合污染指數(shù)進(jìn)行Pearson相關(guān)性分析,結(jié)果見表3.在0.05(雙尾)水平下,太滆運(yùn)河、武進(jìn)港和武宜運(yùn)河這3條區(qū)內(nèi)主要河流綜合污染指數(shù)均與城鎮(zhèn)化因子顯著負(fù)相關(guān),與產(chǎn)業(yè)結(jié)構(gòu)因子顯著正相關(guān).武南區(qū)域內(nèi)城鎮(zhèn)化發(fā)展的穩(wěn)步推進(jìn)對(duì)重點(diǎn)河流水質(zhì)改善具積極影響,產(chǎn)業(yè)結(jié)構(gòu)由工業(yè)、農(nóng)業(yè)向第三產(chǎn)業(yè)的轉(zhuǎn)變也對(duì)區(qū)內(nèi)河流水質(zhì)有正向作用.政府所推進(jìn)的新型城鎮(zhèn)化建設(shè)與產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化理念對(duì)區(qū)域內(nèi)主干河流水質(zhì)改善成效較顯著.從污染源排放角度探究內(nèi)在原因可知,區(qū)內(nèi)城鎮(zhèn)化的推進(jìn)可大幅減少農(nóng)村生活源的直接排放,并減少部分農(nóng)業(yè)面源向地表徑流的排放,從而降低河流污染負(fù)荷,提升河流水質(zhì);產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化可使區(qū)域內(nèi)工業(yè)企業(yè)數(shù)量減少,產(chǎn)業(yè)聚集程度提高,從而降低工業(yè)源的入河負(fù)荷,促進(jìn)區(qū)域內(nèi)水質(zhì)改善.

    表2 各因子主成分分析結(jié)果

    注:*表示相關(guān)系數(shù)30.75.

    表3 內(nèi)梅羅指數(shù)與因子Pearson相關(guān)性分析結(jié)果

    注:*表示在0.05水平下顯著相關(guān),**表示在0.01水平下顯著相關(guān).

    同時(shí),各因子與永安河年均綜合污染指數(shù)不存在顯著的相關(guān)關(guān)系,這主要是因?yàn)橛腊埠恿髁枯^小,監(jiān)測數(shù)據(jù)波動(dòng)較大,且河流水質(zhì)也受相關(guān)治理工程影響.因此,本研究將單獨(dú)分析永安河水質(zhì).

    2.4 重污染河流(永安河)各指標(biāo)季節(jié)性變化及相關(guān)性分析

    永安河位于武進(jìn)區(qū)內(nèi),北連采菱港,南接太滆運(yùn)河.流經(jīng)湖塘、南夏墅、禮嘉和前黃鎮(zhèn),屬太湖流域武澄錫虞區(qū)骨干調(diào)節(jié)河道.由于河流沿岸化工、電鍍、印染企業(yè)繁多,且排放管控難,永安河及其主要支浜的水質(zhì)黑臭問題一度非常嚴(yán)重,其所排放的污染份額曾占據(jù)太滆運(yùn)河的80%[31].由前述分析可知,永安河各項(xiàng)水質(zhì)指標(biāo)濃度遠(yuǎn)高于其它幾條河流,且水質(zhì)變幅較大.根據(jù)當(dāng)?shù)亟邓?將一年分為雨季(4~9月)和旱季(1~3月,10~12月)兩個(gè)時(shí)段,并分別分析永安河氨氮、總氮、總磷、高錳酸鹽指數(shù)變化趨勢如圖4所示.

    分析結(jié)果顯示,永安河氨氮和總氮季節(jié)性變化較明顯,雨季濃度明顯低于旱季,這主要是由于降水對(duì)各污染指標(biāo)有一定的稀釋作用.就時(shí)間跨度而言,2006~2018年,2項(xiàng)指標(biāo)呈好轉(zhuǎn)趨勢,但在2017年有明顯反彈.河流氨氮和總氮濃度的降低主要是由于當(dāng)?shù)卣畬?duì)永安河采取了治理措施,政府對(duì)河流沿岸高污染企業(yè)進(jìn)行了專項(xiàng)整治,規(guī)范了工業(yè)廢水處理要求,并關(guān)停了不合格企業(yè).而2017年的反彈則是永安河清淤工程所致暫時(shí)現(xiàn)象.永安河總磷濃度和高錳酸鹽指數(shù)的季節(jié)性變化不明顯,雨季和旱季濃度相差不大.就時(shí)間跨度而言,總磷濃度和高錳酸鹽指數(shù)在2006~2018年間變化趨勢不明顯,說明2項(xiàng)污染指標(biāo)與工廠點(diǎn)源排放關(guān)系不大,更多來自河流沿岸面源污染.

    永安河沿岸,工業(yè)點(diǎn)源與農(nóng)業(yè)面源分布錯(cuò)綜復(fù)雜,所流經(jīng)的前黃鎮(zhèn)為武南區(qū)域最大的農(nóng)業(yè)區(qū),經(jīng)估算,前黃鎮(zhèn)所排農(nóng)業(yè)面源占比高達(dá)全區(qū)的44%.降水既能對(duì)河流產(chǎn)生稀釋效應(yīng),也會(huì)因地表徑流將更多污染源帶入河流.因此要使河流水質(zhì)進(jìn)一步好轉(zhuǎn),在進(jìn)行工業(yè)點(diǎn)源整改與河道清淤的同時(shí),也應(yīng)注重河流沿岸農(nóng)業(yè)面源管理.

    為進(jìn)一步探究各水質(zhì)污染指標(biāo)間的內(nèi)部聯(lián)系,對(duì)高錳酸鹽指數(shù)、氨氮、總磷、總氮4項(xiàng)指標(biāo)進(jìn)行Pearson相關(guān)性分析,結(jié)果見表4.高錳酸鹽指數(shù)與總氮、氨氮間存在顯著相關(guān)性,而總磷與其他3項(xiàng)指標(biāo)相關(guān)性均不顯著.該結(jié)果說明河流中高錳酸鹽指數(shù)、總氮和氨氮相互影響,指標(biāo)間或存在同源關(guān)系.對(duì)比區(qū)域內(nèi)污染源估算結(jié)果,永安河流經(jīng)地主要污染源為農(nóng)業(yè)源和生活源,各污染源對(duì)高錳酸鹽指數(shù)、氨氮和總氮濃度均有貢獻(xiàn),為3項(xiàng)指標(biāo)的共同污染源;而水產(chǎn)養(yǎng)殖源、種植源和畜禽養(yǎng)殖源為影響總磷濃度的主要污染源,其他污染源對(duì)總磷排放的貢獻(xiàn)較小.

    表4 各水質(zhì)指標(biāo)Pearson相關(guān)性分析結(jié)果

    注:*表示在0.05水平下顯著相關(guān),**表示在0.01水平下顯著相關(guān).

    3 結(jié)論

    3.1 在2006~2018年間太滆運(yùn)河、武宜運(yùn)河、武進(jìn)港、永安河4條河流水質(zhì)總體好轉(zhuǎn),修正內(nèi)梅羅指數(shù)降幅分別達(dá)36.2%,31.5%,56.4%,48.7%.受河流清淤工程影響,永安河水質(zhì)在2017年有所反彈.

    3.2 2006~2018年間,4條河流氨氮濃度、總氮濃度和高錳酸鹽指數(shù)下降趨勢顯著(<0.05),武進(jìn)港和武宜運(yùn)河總磷濃度下降明顯,太滆運(yùn)河和永安河不同年份總磷濃度存在不規(guī)則波動(dòng).

    3.3 永安河氨氮濃度和總氮濃度季節(jié)性變化明顯,雨季濃度明顯低于旱季,總磷濃度和高錳酸鹽指數(shù)不存在明顯季節(jié)性變化趨勢.氨氮、總氮和高錳酸鹽指數(shù)間相關(guān)性較顯著,各污染物或存在同源關(guān)系,而總磷與其他3項(xiàng)指標(biāo)間不存在顯著相關(guān)關(guān)系.

    3.4 武南區(qū)域內(nèi)城鎮(zhèn)化發(fā)展對(duì)區(qū)域水環(huán)境改善有積極影響,產(chǎn)業(yè)結(jié)構(gòu)由傳統(tǒng)工業(yè)、農(nóng)業(yè)向第三產(chǎn)業(yè)的轉(zhuǎn)變有助于河流水質(zhì)好轉(zhuǎn).

    [1] Su L, Xue Y, Li L, et al. Microplastics in taihu lake, China [J]. Environmental pollution, 2016,216:711-719.

    [2] 蔡履冰.太湖流域水體富營養(yǎng)化成因及防治對(duì)策的初步研究[J]. 中國環(huán)境監(jiān)測, 2003,3(19):52-55. Cai L B. Preliminary study on the cause and the preventing and curing of the water body eutrophication of the Taihu Basin [J]. Environmental Monitoring in China, 2003,3(19):52-55.

    [3] 李 未,秦伯強(qiáng).太湖梅梁灣富營養(yǎng)化主要驅(qū)動(dòng)因子的多時(shí)間尺度分析 [J]. 湖泊科學(xué), 2012,24(6):865-874. Li W, Qin B Q. Multiple temporal scale analysis of main driving factors of eutrophication in Meiliang Bay, Lake Taihu [J]. Journal of Lake Sciences, 2012,24(6):865-874.

    [4] 毛新偉,徐 楓,徐 彬,等.太湖水質(zhì)及富營養(yǎng)化變化趨勢分析 [J]. 水資源保護(hù), 2009,25(1):48-51. Mao X W, Xu F, Xu B, et al. Changes of water quality and eutrophication in Taihu Lake [J]. Water Resources Protection, 2009, 25(1):48-51.

    [5] 朱廣偉.太湖富營養(yǎng)化現(xiàn)狀及原因分析 [J]. 湖泊科學(xué), 2008,(1): 21-26. Zhu G W. Eutrophic status and causing factors for a large, shallow and subtropical Lake Taihu, China [J]. Journal of Lake Sciences, 2008,(1): 21-26.

    [6] 李娟英,曹宏宇,崔 昱,等.太湖流域主要水系水環(huán)境特征分析與富營養(yǎng)化評(píng)價(jià)[J]. 水生態(tài)學(xué)雜志, 2012,33(4):7-13. Li J Y, Cao H Y, Cui Y, et al. Water environmental characteristics analysis and eutrophication assessment on lakes in Taihu Basin [J]. Journal of Hydroecology, 2012,33(4):7-13.

    [7] Xu X, Yang G, Tan Y, et al. Ecological risk assessment of ecosystem services in the Taihu Lake Basin of China from 1985 to 2020 [J]. Science of the Total Environment, 2016,554:7-16.

    [8] 于東升,袁宏林,張 穎,等.太湖西岸水質(zhì)變化趨勢及主要驅(qū)動(dòng)因子 [J]. 環(huán)境污染與防治, 2017,39(10):1063-1066+1070. Yu D S, Yuan H L, Zhang Y, et al. Trend and dominating factors of water quality change in the west bank of Taihu Lake [J]. Environmental Pollution & Control, 2017,39(10):1063-1066+1070.

    [9] 左一鳴,崔廣柏,馮 健.太湖水質(zhì)指標(biāo)相關(guān)關(guān)系的研究 [J]. 人民長江, 2005,(10):29-30. Zuo Y M, Cui G B, Feng J. Study on the correlation of water quality indicators in Taihu Lake [J]. Yangtze River, 2005,(10):29-30.

    [10] 朱 偉,談?dòng)狼?王若辰,等.太湖典型區(qū)2010~2017年間水質(zhì)變化趨勢及異常分析 [J]. 湖泊科學(xué), 2018,30(2):296-305. Zhu W, Tan Y Q, Wang R C, et al. The trend of water quality variation and analysis in typical area of Lake Taihu, 2010~2017 [J]. Journal of Lake Sciences, 2018,30(2):296-305.

    [11] 張 宇,凌 虹,吳俊鋒,等.武宜運(yùn)河小流域平原河網(wǎng)地區(qū)水環(huán)境質(zhì)量現(xiàn)狀分析及改善對(duì)策研究 [J]. 環(huán)境科學(xué)與管理, 2015,40(12): 20-25. Zhang Y, Ling H, Wu J F, et al. Present situation and countermeasures of water environment quality in plain river network region in small watershed of Wuyi River [J]. Environmental Science and Management, 2015,40(12): 20-25.

    [12] 常州市統(tǒng)計(jì)局,國家統(tǒng)計(jì)局常州調(diào)查隊(duì).常州統(tǒng)計(jì)年鑒 [J]. 中國統(tǒng)計(jì)出版社, 2007~2018. Changzhou Statistical Burea, National Bureau of Statistics, Changzhou Division. Changzhou Statistical Yearbook [J]. China Statistics Press, 2007~2018.

    [13] 郭 晶,王丑明,黃代中,等.洞庭湖水污染特征及水質(zhì)評(píng)價(jià)[J]. 環(huán)境化學(xué), 2019,38(1):152-160. Guo J, Wang C M, Huang D Z, et al. Pollution characterization and water quality assessment of Dongting Lake [J]. Environmental Chemistry, 2019,38(1):152-160.

    [14] Lumb A, Halliwell D, Sharma T. Application of CCME water quality index to monitor water quality: a case study of the Mackenzie River Basin, Canada [J]. Environmental Monitoring and Assessment, 2006, 113(1/2/3):411-429.

    [15] 付 寧,任學(xué)昌,朱文萍,等.模糊數(shù)學(xué)法在蘭州市地表水環(huán)境質(zhì)量評(píng)價(jià)中的應(yīng)用分析研究[J]. 環(huán)境科學(xué)與管理, 2016,41(12):168-172+ 186. Fu N, Ren X C, Zhu W P, et al. Application analysis of fuzzy mathematics in water quality evaluation of surface water in Lanzhou [J]. Environmental Science and Management, 2016,41(12):168-172+ 186.

    [16] 李 霄,柴 璐,王曉光,等.基于層次分析法的丹東地區(qū)地下水污染防治區(qū)劃[J]. 地質(zhì)與資源, 2018,27(4):396-405. Li X, Chai L, Wang X G, et al. Regionalization of groundwater pollution prevention in Dandong area based on analytic hierarchy process [J]. Geology and Resources, 2018,27(4):396-405.

    [17] 張亞麗,周 揚(yáng),程 真,等.不同水質(zhì)評(píng)價(jià)方法在丹江口流域水質(zhì)評(píng)價(jià)中應(yīng)用比較[J]. 中國環(huán)境監(jiān)測, 2015,31(3):58-61. Zhang Y L, Zhou Y, Cheng Z, et al. Comparison and analysis of different evaluation methods for water quality in the Danjiangkou Valley [J]. Environmental Monitoring in China, 2015,31(3):58-61.

    [18] 高紅杰,鄭利杰,嵇曉燕,等.典型城市地表水質(zhì)綜合評(píng)價(jià)方法研究[J]. 中國環(huán)境監(jiān)測, 2017,33(2):55-60. Gao H J, Zheng L J, Ji X Y, et al. Study on comprehensive evaluation methods for surface water of typical cities [J]. Environmental Monitoring in China, 2017,33(2):55-60.

    [19] 吳喜軍,董 穎,張亞寧.改進(jìn)的內(nèi)梅羅污染指數(shù)法在黃河干流水質(zhì)評(píng)價(jià)中的應(yīng)用[J]. 節(jié)水灌溉, 2018,(10):51-53+58. Wu X J, Dong Y, Zhang Y N. Application of improved nemerow pollution index method in water quality evaluation of Yellow River [J]. Water Saving Irrigation, 2018,(10):51-53+58.

    [20] 王順久.長江上游川江段氣溫、降水及徑流變化趨勢分析[J]. 資源科學(xué), 2009,31(7):1142-1149. Wang S J. Changing pattern of the temperature, precipitation and runoff in Chuanjiang section of the Yangtze River [J]. Resources Science, 2009,31(7):1142-1149.

    [21] 賈立斌,吳偉宏,袁國華.基于Mann-Kendall的中國近岸海域海洋生態(tài)環(huán)境承載力評(píng)價(jià)與預(yù)警[J]. 生態(tài)經(jīng)濟(jì), 2019,35(2):208-213+224. Jia L B, Wu W H, Yuan G H. Carrying capacity assessment and warning for China nearshore marine ecological environment based on Mann–Kendall [J]. Ecological Economy, 2019,35(2):208-213+224.

    [22] 商 博,于光金,王桂勛,等.基于PCA的區(qū)域環(huán)境質(zhì)量綜合評(píng)價(jià)及應(yīng)用實(shí)例研究[J]. 中國環(huán)境監(jiān)測, 2013,29(5):12-15. Shang B, Yu G J, Wang G X, et al. Research-based on the principal component analysis in regional environmental quality evaluation and application examples [J]. Environmental Monitoring in China, 2013,29(5):12-15.

    [23] 侯佳辛.基于內(nèi)梅羅污染指數(shù)法的區(qū)域水環(huán)境評(píng)價(jià)——以瑪納斯河流域?yàn)槔齕J]. 內(nèi)蒙古水利, 2012,(5):8-9. Hou J X. Regional water environment assessment based on Nemero Pollution Index Method: A case study of Manas River Basin [J]. Inner Mongolia Water Resources, 2012,(5):8-9.

    [24] GB3838-2002 地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)[S]. GB3838-2002 Surface Water Quality Standard [S].

    [25] 徐憲根,曹慧敏,顧禮明,等.太滆運(yùn)河沿程水質(zhì)變化及原因解析[J]. 廣東化工, 2017,44(12):197-198+194. Xu X G, Cao H M, Gu L M, et al. The water quality changes and cause analysis of Taige Canal [J]. Guangdong Chemical Industry, 2017, 44(12):197-198+194.

    [26] 謝文理,田 穎,祁紅娟,等.太滆運(yùn)河氨氮污染時(shí)空特征及來源研究 [J]. 安徽農(nóng)業(yè)科學(xué), 2018,46(22):55-57. Xie W L, Tian Y, Qi H J, et al. Study on spatial-temporal features and source of ammonia-nitrogen pollution of Taige Canal [J]. Journal of Anhui Agricultural Sciences, 2017,44(12):197-198+194.

    [27] 張大偉,李楊帆,孫 翔,等.入太湖河流武進(jìn)港的區(qū)域景觀格局與河流水質(zhì)相關(guān)性分析 [J]. 環(huán)境科學(xué), 2010,31(8):1775-1783. Zhang D W, Li Y F, Sun X, et al. Relationship between landscape pattern and river water quality in Wujingang Region,Taihu Lake Watershed [J]. Environmental Science, 2010,31(8):1775-1783.

    [28] 李菁紅.促進(jìn)武進(jìn)產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化的財(cái)政政策研究[J]. 當(dāng)代經(jīng)濟(jì), 2014,(18):78-79. Li J H. Study on the fiscal policy to promote the industrial structure optimization of Wujin [J]. Contemporary Economics, 2014,(18):78- 79.

    [29] 邱 靈,韓曉旭.創(chuàng)新驅(qū)動(dòng)縣域經(jīng)濟(jì)高質(zhì)量發(fā)展——以江蘇省常州市武進(jìn)區(qū)為例 [J]. 全球化, 2019,(3):73-81+135. Qiu L, Han X X. Innovation drives the high quality development of county economy——Taking Wujin District, Changzhou City, Jiangsu Province as an example [J]. Globalization, 2019,(3):73-81+135.

    [30] 尚 娟,童 瑤.供給側(cè)要素投入對(duì)新型城鎮(zhèn)化發(fā)展影響的測算 [J/OL]. 統(tǒng)計(jì)與決策, 2019,(12):50-53. Shang J, Tong Y. Measurement of the impact of supply-side factor input on the development of new urbanization [J/OL]. Statistics & Decision, 2019,(12):50-53.

    [31] 中國環(huán)境報(bào).江蘇常州治理永安河-控源截污為入湖河流減負(fù) [EB/OL]. http://www.er-china.com/PowerLeader/html/2011/06/20110 623090746.shtml China Environmental News. Governance of Yongan River in Changzhou, Jiangsu Province, Controlling source and pollution for load reduction of river inflowing into the lake [EB/OL]. http:// www.er-china.com/PowerLeader/html/2011/06/20110623090746.shtml

    Spatiotemporal trends in water quality and the effectiveness of water treatment projects for Wunan District, Changzhou.

    KANG Ruo-xi1, XIE Wen-li2, ZHAO Wen-xing1, PANG Xiao-si1, YU Gang1*

    (1.Beijing Key Laboratory for Emerging Organic Contaminants Control, School of Environment, Tsinghua University, Beijing 100084, China;2.Changzhou Environmental Monitoring Center, Changzhou 213001, China)., 2019,39(10):4186~4193

    The aims of this study were to explore spatiotemporal differences in water quality of 4 main rivers (Taige River, Wuyi River, Wujin River and Yongan River) in Wunan District based on 2006~2018 water quality monitoring data, to identify the contributions of different activities on water quality parameters, to evaluate the effects of treatment projects on water quality, and to conduct seasonal and correlation analysis for the most heavily polluted river (i.e., Yongan River). The results showed that water quality of 4 rivers improved from 2006 to 2018, with the modified Nemerow index declining by 36.2%, 31.5%, 56.4% and 48.7%. Water quality of Yongan River deteriorated in 2017 due to dredging activities. The ammonia nitrogen (NH3-N) concentration, total nitrogen (TN) concentration and permanganate index significantly decreased year by year (<0.05), but the total phosphorus (TP) concentrations fluctuated across the same period. The correlation analysis indicated that NH3-N, TN and permanganate index may stem from the same sources. NH3-N and TN exhibited an obvious seasonal variation, but no clear trend was observed for TP and permanganate index. Overall, increased urbanization, decreased traditional industry, and decreased agricultural activity had positive effects on improving water quality in Wunan District.

    Wunan District;water quality;variation tendency;driving force

    X522

    A

    1000-6923(2019)10-4186-08

    康若熙(1993-),女,四川成都人,清華大學(xué)環(huán)境學(xué)院碩士研究生,主要從事水質(zhì)數(shù)據(jù)分析研究.

    2019-03-25

    國家水體污染控制與治理科技重大專項(xiàng)(2017ZX07202006)

    * 責(zé)任作者, 教授, yg-den@tsinghua.edu.cn

    猜你喜歡
    內(nèi)梅羅高錳酸鹽永安
    長江經(jīng)濟(jì)帶水質(zhì)監(jiān)測斷面嘉興段評(píng)價(jià)分析
    綠色科技(2022年10期)2022-06-23 10:14:26
    內(nèi)鏡止血在急性非靜脈曲張性上消化道出血治療中的應(yīng)用
    長安永安
    從多元函數(shù)角度探討高錳酸鹽指數(shù)準(zhǔn)確度的影響因素
    云南化工(2021年6期)2021-12-21 07:31:16
    全自動(dòng)高錳酸鹽指數(shù)分析儀在水環(huán)境監(jiān)測中的應(yīng)用
    中國測試(2021年4期)2021-07-16 07:49:04
    一種新內(nèi)梅羅指數(shù)法的修正方法及應(yīng)用
    永安橋 百年永安
    中國公路(2017年18期)2018-01-23 03:00:49
    內(nèi)梅羅指數(shù)法在環(huán)境質(zhì)量評(píng)價(jià)中的適用性與修正原則
    高錳酸鹽指數(shù)測定中若干問題的探討
    竹鄉(xiāng)永安
    国产成人精品无人区| videossex国产| 一区二区三区免费毛片| 免费看日本二区| 久久精品熟女亚洲av麻豆精品| 黄色日韩在线| 视频区图区小说| 久久久久久久精品精品| 性色av一级| 99久久精品一区二区三区| 久久久国产精品麻豆| 中文欧美无线码| 色哟哟·www| 婷婷色综合www| av线在线观看网站| 人人妻人人澡人人爽人人夜夜| 91aial.com中文字幕在线观看| 各种免费的搞黄视频| 六月丁香七月| 久久久久久久大尺度免费视频| 日日摸夜夜添夜夜添av毛片| 国产极品天堂在线| 亚洲欧洲精品一区二区精品久久久 | 99热网站在线观看| 91精品国产国语对白视频| 国产无遮挡羞羞视频在线观看| 纵有疾风起免费观看全集完整版| 国产免费福利视频在线观看| 少妇人妻一区二区三区视频| 女人精品久久久久毛片| 亚洲精品国产成人久久av| 国产av国产精品国产| 777米奇影视久久| 搡老乐熟女国产| 日日摸夜夜添夜夜爱| 丰满少妇做爰视频| 丰满饥渴人妻一区二区三| 久久精品国产亚洲av天美| 秋霞伦理黄片| 岛国毛片在线播放| 亚洲精品,欧美精品| 精品久久久久久久久亚洲| 精品午夜福利在线看| 18禁在线播放成人免费| 久久青草综合色| 自拍欧美九色日韩亚洲蝌蚪91 | 国产黄频视频在线观看| 免费黄网站久久成人精品| 国产极品粉嫩免费观看在线 | 妹子高潮喷水视频| 国产免费又黄又爽又色| 色94色欧美一区二区| 国产免费一级a男人的天堂| 蜜臀久久99精品久久宅男| 成人国产av品久久久| 极品教师在线视频| 精品久久久噜噜| 欧美日本中文国产一区发布| 伊人亚洲综合成人网| 国产熟女午夜一区二区三区 | 亚洲精品日韩av片在线观看| 日日撸夜夜添| 亚洲欧洲精品一区二区精品久久久 | 欧美 日韩 精品 国产| 精品久久久噜噜| 欧美激情极品国产一区二区三区 | 三级国产精品片| 80岁老熟妇乱子伦牲交| 久久久欧美国产精品| 美女脱内裤让男人舔精品视频| 亚洲精品乱码久久久久久按摩| 视频区图区小说| 欧美另类一区| 日韩欧美精品免费久久| 国产精品蜜桃在线观看| 国产av一区二区精品久久| www.av在线官网国产| 国产免费一级a男人的天堂| 国产 一区精品| 日韩中文字幕视频在线看片| 中文字幕av电影在线播放| 欧美 亚洲 国产 日韩一| 丝袜在线中文字幕| 人妻 亚洲 视频| 国产白丝娇喘喷水9色精品| av网站免费在线观看视频| 日韩中字成人| 国产精品欧美亚洲77777| 高清不卡的av网站| 亚洲人与动物交配视频| 国产一区有黄有色的免费视频| 国产伦在线观看视频一区| 欧美日韩综合久久久久久| 日产精品乱码卡一卡2卡三| 亚洲欧美成人综合另类久久久| 国产精品免费大片| 嫩草影院入口| 嘟嘟电影网在线观看| 成人美女网站在线观看视频| 国产黄片美女视频| 2022亚洲国产成人精品| 亚洲精品色激情综合| 久久ye,这里只有精品| 视频区图区小说| 日日啪夜夜撸| 人人妻人人看人人澡| 色5月婷婷丁香| kizo精华| 91成人精品电影| 欧美3d第一页| 老司机亚洲免费影院| 插阴视频在线观看视频| 99久久综合免费| 热re99久久国产66热| √禁漫天堂资源中文www| 大码成人一级视频| 精品国产一区二区久久| 日韩欧美 国产精品| 国产精品一区二区性色av| 高清午夜精品一区二区三区| 亚洲国产精品一区二区三区在线| 三级经典国产精品| 亚洲无线观看免费| 99久久精品一区二区三区| 午夜视频国产福利| 日本黄色片子视频| 日韩一区二区三区影片| 老司机影院成人| 国语对白做爰xxxⅹ性视频网站| 久久99热这里只频精品6学生| 久久久久视频综合| 美女主播在线视频| 成人午夜精彩视频在线观看| 三级国产精品片| 在线观看国产h片| 久久久久视频综合| 国产一区亚洲一区在线观看| 大片电影免费在线观看免费| 久久精品国产自在天天线| freevideosex欧美| av.在线天堂| 欧美日韩一区二区视频在线观看视频在线| 国产欧美亚洲国产| 亚洲,欧美,日韩| 国产色婷婷99| 丝瓜视频免费看黄片| 晚上一个人看的免费电影| 极品教师在线视频| 国产日韩欧美视频二区| 日韩欧美精品免费久久| 大片免费播放器 马上看| 一二三四中文在线观看免费高清| 我要看日韩黄色一级片| 大片电影免费在线观看免费| 国产精品人妻久久久影院| 国产成人精品无人区| 欧美区成人在线视频| 最近中文字幕2019免费版| 伦精品一区二区三区| 寂寞人妻少妇视频99o| 一级,二级,三级黄色视频| 国产永久视频网站| 纵有疾风起免费观看全集完整版| 亚洲欧美成人精品一区二区| 91精品国产九色| 亚洲精品456在线播放app| 精品一区二区三区视频在线| 色婷婷久久久亚洲欧美| 最近2019中文字幕mv第一页| 久久久久国产精品人妻一区二区| 一级av片app| 免费看av在线观看网站| 亚洲精品中文字幕在线视频 | 精品人妻熟女av久视频| 国产在视频线精品| 少妇 在线观看| 日本黄色片子视频| 亚洲国产精品成人久久小说| 五月天丁香电影| 亚洲欧美一区二区三区国产| 99热这里只有是精品在线观看| 日本av手机在线免费观看| 春色校园在线视频观看| 人妻少妇偷人精品九色| 麻豆乱淫一区二区| 最后的刺客免费高清国语| 成年人午夜在线观看视频| 国产中年淑女户外野战色| 国产69精品久久久久777片| 99热全是精品| 国产中年淑女户外野战色| 亚洲精品成人av观看孕妇| av国产精品久久久久影院| 高清欧美精品videossex| 日本欧美国产在线视频| av国产精品久久久久影院| 国产男人的电影天堂91| 成人国产av品久久久| 18禁动态无遮挡网站| 最后的刺客免费高清国语| 精品少妇内射三级| 黄色毛片三级朝国网站 | 国产免费视频播放在线视频| 久久99热这里只频精品6学生| 国产成人精品婷婷| 夫妻午夜视频| 在线观看av片永久免费下载| 亚洲欧美精品专区久久| 久久99精品国语久久久| 亚洲欧美成人综合另类久久久| 91久久精品国产一区二区三区| 久久人人爽av亚洲精品天堂| 日本av免费视频播放| 国产日韩欧美视频二区| 亚洲精品中文字幕在线视频 | 桃花免费在线播放| 日日撸夜夜添| av免费在线看不卡| 国产精品一区二区在线不卡| 国语对白做爰xxxⅹ性视频网站| 国内精品宾馆在线| a级片在线免费高清观看视频| 国产在线男女| 国产熟女欧美一区二区| 日韩欧美 国产精品| 国产精品一区二区在线观看99| 国产高清不卡午夜福利| 永久网站在线| 九草在线视频观看| 国产成人精品久久久久久| 伊人亚洲综合成人网| 精华霜和精华液先用哪个| 亚洲精品自拍成人| 久久久久精品久久久久真实原创| 2018国产大陆天天弄谢| 成人18禁高潮啪啪吃奶动态图 | 亚洲天堂av无毛| 亚洲av.av天堂| 99热网站在线观看| 精品久久久久久久久av| 久久午夜综合久久蜜桃| 黄色毛片三级朝国网站 | a级毛色黄片| 少妇熟女欧美另类| 高清毛片免费看| 99久久人妻综合| 熟妇人妻不卡中文字幕| 国产精品久久久久久久久免| 91精品国产九色| 日韩欧美一区视频在线观看 | 蜜桃在线观看..| 99re6热这里在线精品视频| 伦理电影大哥的女人| 插逼视频在线观看| 中文天堂在线官网| 最近中文字幕2019免费版| 亚洲一区二区三区欧美精品| 亚洲欧美成人综合另类久久久| 亚洲综合色惰| 国国产精品蜜臀av免费| 97在线视频观看| 久久影院123| 午夜激情久久久久久久| 亚洲自偷自拍三级| 人人妻人人添人人爽欧美一区卜| 成人无遮挡网站| 另类精品久久| 国产淫片久久久久久久久| 国产成人免费无遮挡视频| 中文字幕精品免费在线观看视频 | 在线观看免费高清a一片| 久久97久久精品| 人妻系列 视频| 乱码一卡2卡4卡精品| 人妻系列 视频| 中文字幕久久专区| 一个人看视频在线观看www免费| 国产精品一区二区三区四区免费观看| 日韩一本色道免费dvd| 综合色丁香网| 十八禁高潮呻吟视频 | 美女xxoo啪啪120秒动态图| 亚洲精品第二区| 中国美白少妇内射xxxbb| 亚洲美女黄色视频免费看| 久久6这里有精品| 免费黄频网站在线观看国产| 这个男人来自地球电影免费观看 | 欧美bdsm另类| 最近的中文字幕免费完整| 日韩在线高清观看一区二区三区| 亚洲欧洲国产日韩| 夜夜骑夜夜射夜夜干| 欧美97在线视频| 亚洲一区二区三区欧美精品| 久久精品久久精品一区二区三区| 成年女人在线观看亚洲视频| 国产熟女午夜一区二区三区 | 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 最新中文字幕久久久久| 免费看日本二区| 免费观看av网站的网址| 99九九线精品视频在线观看视频| 免费黄网站久久成人精品| 自拍偷自拍亚洲精品老妇| 免费观看无遮挡的男女| 丁香六月天网| av福利片在线| 美女cb高潮喷水在线观看| 国产淫语在线视频| 欧美日韩精品成人综合77777| 自拍偷自拍亚洲精品老妇| 日产精品乱码卡一卡2卡三| 亚洲高清免费不卡视频| 自线自在国产av| 天堂中文最新版在线下载| 久久影院123| 欧美 日韩 精品 国产| 一级毛片黄色毛片免费观看视频| 天堂8中文在线网| 十八禁高潮呻吟视频 | 久久久久久伊人网av| 成年人午夜在线观看视频| 日韩成人伦理影院| 亚洲伊人久久精品综合| 一级二级三级毛片免费看| 特大巨黑吊av在线直播| 国产女主播在线喷水免费视频网站| 成人漫画全彩无遮挡| 三上悠亚av全集在线观看 | 男女啪啪激烈高潮av片| 久久久久久久精品精品| 国产精品不卡视频一区二区| 99热网站在线观看| 亚洲成人一二三区av| 欧美日韩一区二区视频在线观看视频在线| 国产高清国产精品国产三级| 国产日韩欧美亚洲二区| 人妻一区二区av| 日本wwww免费看| 丝袜脚勾引网站| 精品人妻熟女毛片av久久网站| 七月丁香在线播放| 下体分泌物呈黄色| 岛国毛片在线播放| 伦理电影免费视频| 蜜桃久久精品国产亚洲av| 高清黄色对白视频在线免费看 | 久久精品熟女亚洲av麻豆精品| 少妇人妻一区二区三区视频| 女性被躁到高潮视频| 欧美日本中文国产一区发布| 国产精品免费大片| 亚洲av不卡在线观看| 全区人妻精品视频| 美女xxoo啪啪120秒动态图| 女的被弄到高潮叫床怎么办| 成人亚洲欧美一区二区av| 欧美另类一区| 久久精品国产鲁丝片午夜精品| 日韩人妻高清精品专区| 啦啦啦视频在线资源免费观看| 国产精品麻豆人妻色哟哟久久| 午夜日本视频在线| 最近手机中文字幕大全| 亚洲在久久综合| 日韩三级伦理在线观看| 欧美3d第一页| 日本-黄色视频高清免费观看| 久久久国产一区二区| 另类精品久久| 国内少妇人妻偷人精品xxx网站| 99九九线精品视频在线观看视频| 精品亚洲成a人片在线观看| 国内少妇人妻偷人精品xxx网站| 嫩草影院新地址| 久久久久久久亚洲中文字幕| 如日韩欧美国产精品一区二区三区 | 免费黄网站久久成人精品| 国产女主播在线喷水免费视频网站| 色视频在线一区二区三区| 在线观看人妻少妇| 亚洲精品久久久久久婷婷小说| 夜夜骑夜夜射夜夜干| 80岁老熟妇乱子伦牲交| 美女大奶头黄色视频| av在线播放精品| 又黄又爽又刺激的免费视频.| 亚洲精品亚洲一区二区| 最后的刺客免费高清国语| 中国美白少妇内射xxxbb| 国产永久视频网站| 久久6这里有精品| 国产伦精品一区二区三区四那| 欧美日韩视频精品一区| 国产淫片久久久久久久久| 人妻制服诱惑在线中文字幕| 亚洲av福利一区| 欧美丝袜亚洲另类| 中文天堂在线官网| 我要看黄色一级片免费的| 亚洲久久久国产精品| 精品午夜福利在线看| 欧美日韩亚洲高清精品| 国产黄色免费在线视频| 欧美日韩国产mv在线观看视频| 亚洲精品久久午夜乱码| freevideosex欧美| 丰满少妇做爰视频| 久热这里只有精品99| 国产精品嫩草影院av在线观看| 日韩视频在线欧美| 日本av免费视频播放| 下体分泌物呈黄色| 国产精品福利在线免费观看| 精品99又大又爽又粗少妇毛片| 视频区图区小说| 老熟女久久久| 又粗又硬又长又爽又黄的视频| 两个人免费观看高清视频 | 亚洲国产精品999| 9色porny在线观看| 国产在线免费精品| 最近手机中文字幕大全| 精品人妻偷拍中文字幕| 超碰97精品在线观看| 69精品国产乱码久久久| 久久精品夜色国产| 欧美精品一区二区大全| 人妻人人澡人人爽人人| 国产免费又黄又爽又色| 国产一区二区在线观看av| 国产亚洲一区二区精品| 一级毛片黄色毛片免费观看视频| 黑人猛操日本美女一级片| 久久久久久久大尺度免费视频| 国产一区有黄有色的免费视频| 日韩中字成人| 建设人人有责人人尽责人人享有的| 久久人妻熟女aⅴ| 国产精品国产三级国产av玫瑰| 青春草视频在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 国产69精品久久久久777片| 国产在线免费精品| 久久午夜综合久久蜜桃| 久久亚洲国产成人精品v| 日韩不卡一区二区三区视频在线| 日本午夜av视频| 日韩亚洲欧美综合| a级片在线免费高清观看视频| 少妇的逼好多水| 热99国产精品久久久久久7| 亚洲图色成人| 3wmmmm亚洲av在线观看| 亚洲av不卡在线观看| 日本与韩国留学比较| 国产成人aa在线观看| a级片在线免费高清观看视频| 免费观看性生交大片5| 国产av一区二区精品久久| 视频中文字幕在线观看| 国产成人午夜福利电影在线观看| 丰满人妻一区二区三区视频av| 亚洲三级黄色毛片| 黑人高潮一二区| 亚洲美女搞黄在线观看| 国产精品国产三级国产av玫瑰| 精品99又大又爽又粗少妇毛片| 午夜av观看不卡| 久久精品久久久久久久性| 久久热精品热| 欧美日韩av久久| 国产黄色免费在线视频| 最近的中文字幕免费完整| 国产午夜精品久久久久久一区二区三区| 下体分泌物呈黄色| 中文字幕人妻丝袜制服| a 毛片基地| 亚洲第一区二区三区不卡| 中文字幕av电影在线播放| 99视频精品全部免费 在线| 综合色丁香网| 成人免费观看视频高清| 亚洲婷婷狠狠爱综合网| 亚洲性久久影院| 一区二区三区免费毛片| 免费av不卡在线播放| 美女国产视频在线观看| 丝袜喷水一区| 亚洲综合色惰| 久久国产精品男人的天堂亚洲 | 亚洲性久久影院| 成人18禁高潮啪啪吃奶动态图 | videos熟女内射| 日韩伦理黄色片| 高清在线视频一区二区三区| 色吧在线观看| 午夜免费男女啪啪视频观看| 免费久久久久久久精品成人欧美视频 | 婷婷色av中文字幕| 在线天堂最新版资源| 好男人视频免费观看在线| 免费观看av网站的网址| 免费大片18禁| 亚洲国产日韩一区二区| 插阴视频在线观看视频| 午夜福利在线观看免费完整高清在| 男男h啪啪无遮挡| 在线看a的网站| 日韩亚洲欧美综合| 伊人亚洲综合成人网| 欧美xxxx性猛交bbbb| 亚洲欧美精品自产自拍| 啦啦啦啦在线视频资源| 国产黄频视频在线观看| 午夜老司机福利剧场| 亚洲精品日韩在线中文字幕| 国产日韩欧美视频二区| 午夜福利网站1000一区二区三区| 欧美成人精品欧美一级黄| 精品久久久久久久久亚洲| 婷婷色av中文字幕| 色吧在线观看| 亚洲第一av免费看| 一本—道久久a久久精品蜜桃钙片| 国产午夜精品一二区理论片| av天堂久久9| 波野结衣二区三区在线| 欧美日韩国产mv在线观看视频| 黄色日韩在线| 欧美日韩av久久| 亚洲av国产av综合av卡| 午夜老司机福利剧场| 国产亚洲91精品色在线| 午夜老司机福利剧场| 国产熟女欧美一区二区| 午夜老司机福利剧场| 高清午夜精品一区二区三区| videossex国产| 九九在线视频观看精品| 青春草亚洲视频在线观看| 欧美高清成人免费视频www| 各种免费的搞黄视频| 丁香六月天网| 爱豆传媒免费全集在线观看| 久久久久久久久久成人| 国产亚洲精品久久久com| 亚洲欧美精品专区久久| 国产精品麻豆人妻色哟哟久久| av播播在线观看一区| 亚洲精品第二区| 纵有疾风起免费观看全集完整版| 精品一区二区三区视频在线| 交换朋友夫妻互换小说| 在线 av 中文字幕| 午夜福利在线观看免费完整高清在| 夜夜骑夜夜射夜夜干| 亚洲经典国产精华液单| 18禁在线无遮挡免费观看视频| 成人免费观看视频高清| 美女福利国产在线| 免费av不卡在线播放| 国精品久久久久久国模美| 亚洲经典国产精华液单| 99九九在线精品视频 | 国产男人的电影天堂91| 黄色欧美视频在线观看| 日本猛色少妇xxxxx猛交久久| 国产精品国产av在线观看| 激情五月婷婷亚洲| 熟女av电影| 欧美日韩视频高清一区二区三区二| 人人妻人人爽人人添夜夜欢视频 | 日韩三级伦理在线观看| 偷拍熟女少妇极品色| 你懂的网址亚洲精品在线观看| 国产精品伦人一区二区| 亚洲欧洲国产日韩| 蜜桃在线观看..| 国产精品女同一区二区软件| 欧美三级亚洲精品| 久久精品夜色国产| av在线观看视频网站免费| 丝袜喷水一区| 丰满乱子伦码专区| 麻豆成人av视频| 国产黄片视频在线免费观看| 在线亚洲精品国产二区图片欧美 | 大香蕉久久网| 丝袜脚勾引网站| 亚洲真实伦在线观看| 精品人妻一区二区三区麻豆| 国产视频首页在线观看| 精品人妻熟女毛片av久久网站| 亚洲国产精品成人久久小说| 亚洲精品,欧美精品| 热re99久久精品国产66热6| 五月玫瑰六月丁香| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩卡通动漫| 大陆偷拍与自拍| 妹子高潮喷水视频| www.色视频.com| 久久久国产欧美日韩av| 亚洲,一卡二卡三卡| 波野结衣二区三区在线| 亚洲综合精品二区| 亚洲国产日韩一区二区| 在线免费观看不下载黄p国产| 一区在线观看完整版| 高清毛片免费看| 亚洲国产欧美在线一区| 亚洲三级黄色毛片|