• <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)永安
    欧美不卡视频在线免费观看| 亚洲av成人av| 夜夜爽天天搞| 老司机福利观看| 亚洲国产日韩欧美精品在线观看 | 欧美又色又爽又黄视频| 波多野结衣高清作品| 最近最新中文字幕大全免费视频| 亚洲欧美激情综合另类| 国产伦精品一区二区三区四那| 日本成人三级电影网站| 婷婷六月久久综合丁香| 亚洲精品中文字幕一二三四区| 亚洲人成电影免费在线| 国产欧美日韩精品一区二区| 亚洲一区二区三区色噜噜| 久久精品国产亚洲av香蕉五月| 叶爱在线成人免费视频播放| 亚洲人成电影免费在线| 天堂av国产一区二区熟女人妻| 久久精品91蜜桃| 久久亚洲精品不卡| 欧美zozozo另类| 国产单亲对白刺激| 手机成人av网站| 美女免费视频网站| 久久久久久久精品吃奶| 美女扒开内裤让男人捅视频| 白带黄色成豆腐渣| 男女做爰动态图高潮gif福利片| 日日干狠狠操夜夜爽| 亚洲一区高清亚洲精品| 韩国av一区二区三区四区| 男人舔女人的私密视频| 欧美一区二区国产精品久久精品| 日韩有码中文字幕| 天堂av国产一区二区熟女人妻| 亚洲中文字幕一区二区三区有码在线看 | 国产精品女同一区二区软件 | 伊人久久大香线蕉亚洲五| 美女高潮喷水抽搐中文字幕| 桃红色精品国产亚洲av| 夜夜躁狠狠躁天天躁| 免费av不卡在线播放| av国产免费在线观看| 国内毛片毛片毛片毛片毛片| 热99re8久久精品国产| 国内精品久久久久精免费| 熟女电影av网| 老汉色av国产亚洲站长工具| 免费观看的影片在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩东京热| 欧美日本亚洲视频在线播放| 男女之事视频高清在线观看| 亚洲av电影不卡..在线观看| 麻豆成人av在线观看| 国产又色又爽无遮挡免费看| 国产精品久久视频播放| 中文字幕最新亚洲高清| 欧美三级亚洲精品| 床上黄色一级片| 精品一区二区三区视频在线 | av在线蜜桃| 亚洲,欧美精品.| 一本久久中文字幕| 91老司机精品| 人妻夜夜爽99麻豆av| 亚洲国产欧洲综合997久久,| 久久久国产成人免费| 精品一区二区三区av网在线观看| 99久久久亚洲精品蜜臀av| 美女高潮喷水抽搐中文字幕| 久久久国产成人免费| 精品国产乱子伦一区二区三区| 成年女人看的毛片在线观看| 亚洲精品色激情综合| 俺也久久电影网| 亚洲欧洲精品一区二区精品久久久| 国产精品1区2区在线观看.| avwww免费| 欧美激情在线99| 99久久综合精品五月天人人| 国产精品乱码一区二三区的特点| 一级a爱片免费观看的视频| 噜噜噜噜噜久久久久久91| 18禁裸乳无遮挡免费网站照片| 在线观看日韩欧美| 免费在线观看成人毛片| 欧美黄色淫秽网站| 级片在线观看| 岛国在线观看网站| 亚洲av成人不卡在线观看播放网| 亚洲国产精品sss在线观看| 国产私拍福利视频在线观看| 欧美日韩乱码在线| 看免费av毛片| 在线观看舔阴道视频| www.自偷自拍.com| 午夜福利高清视频| 国产高清激情床上av| 久久久久久大精品| 欧美zozozo另类| 国产亚洲精品久久久com| 日韩有码中文字幕| 国产欧美日韩精品一区二区| 国产乱人伦免费视频| 悠悠久久av| 国产一区在线观看成人免费| 色综合亚洲欧美另类图片| 中文字幕熟女人妻在线| 久久中文字幕人妻熟女| 亚洲av第一区精品v没综合| 在线看三级毛片| 国产 一区 欧美 日韩| 久久久久久久久中文| 91九色精品人成在线观看| 午夜福利高清视频| 国产精品一及| 一个人免费在线观看的高清视频| 男女视频在线观看网站免费| 婷婷丁香在线五月| 成人精品一区二区免费| 婷婷精品国产亚洲av| 999精品在线视频| 国产一区二区在线观看日韩 | 日本成人三级电影网站| 热99在线观看视频| 九九热线精品视视频播放| 嫁个100分男人电影在线观看| x7x7x7水蜜桃| 一级毛片精品| 99视频精品全部免费 在线 | 日本黄色片子视频| 精品久久久久久久久久免费视频| 欧美日韩国产亚洲二区| 国产精品国产高清国产av| 美女免费视频网站| 亚洲av免费在线观看| 国产黄色小视频在线观看| 男人和女人高潮做爰伦理| 久久久久精品国产欧美久久久| 成人永久免费在线观看视频| 老汉色∧v一级毛片| 此物有八面人人有两片| 在线永久观看黄色视频| АⅤ资源中文在线天堂| 色视频www国产| 99久久精品一区二区三区| 亚洲av成人不卡在线观看播放网| 免费搜索国产男女视频| 成人一区二区视频在线观看| 日韩欧美精品v在线| 欧美在线一区亚洲| 男女那种视频在线观看| av天堂在线播放| 91在线精品国自产拍蜜月 | 搞女人的毛片| 久久久久久久久久黄片| 亚洲欧美激情综合另类| 一个人看的www免费观看视频| 丰满人妻一区二区三区视频av | 成人鲁丝片一二三区免费| 国产精品国产高清国产av| 国产精品久久久久久精品电影| 男人的好看免费观看在线视频| 中文字幕久久专区| 国产一区二区在线观看日韩 | 欧美日韩福利视频一区二区| 国产亚洲精品综合一区在线观看| 亚洲美女视频黄频| 国产精品美女特级片免费视频播放器 | 国产伦精品一区二区三区视频9 | 国产成人一区二区三区免费视频网站| 日韩精品中文字幕看吧| 国产一区二区在线av高清观看| av片东京热男人的天堂| 在线十欧美十亚洲十日本专区| 久久这里只有精品中国| 女警被强在线播放| 日韩国内少妇激情av| 国产不卡一卡二| 18禁裸乳无遮挡免费网站照片| 人妻丰满熟妇av一区二区三区| 精品人妻1区二区| 欧美zozozo另类| 国产 一区 欧美 日韩| 亚洲专区中文字幕在线| 中亚洲国语对白在线视频| 亚洲人成网站高清观看| 成人国产一区最新在线观看| 国产精品女同一区二区软件 | 国产欧美日韩精品亚洲av| 国产成人av教育| 亚洲最大成人中文| 夜夜看夜夜爽夜夜摸| 国产一区二区在线av高清观看| 欧美中文日本在线观看视频| 日本黄大片高清| av在线天堂中文字幕| 法律面前人人平等表现在哪些方面| 亚洲午夜精品一区,二区,三区| 国产极品精品免费视频能看的| 久久久久久久久免费视频了| 亚洲自偷自拍图片 自拍| 婷婷精品国产亚洲av在线| 色综合站精品国产| 免费在线观看成人毛片| 免费在线观看视频国产中文字幕亚洲| 亚洲成人久久性| 1024香蕉在线观看| 国产亚洲av嫩草精品影院| svipshipincom国产片| 亚洲国产精品合色在线| 久久国产乱子伦精品免费另类| 狠狠狠狠99中文字幕| 久久精品国产99精品国产亚洲性色| 在线十欧美十亚洲十日本专区| 亚洲片人在线观看| 欧美激情在线99| 成熟少妇高潮喷水视频| 免费无遮挡裸体视频| 亚洲片人在线观看| 最近最新免费中文字幕在线| 午夜福利成人在线免费观看| 成人欧美大片| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 在线观看日韩欧美| 黄色日韩在线| 国内毛片毛片毛片毛片毛片| 亚洲 欧美一区二区三区| av国产免费在线观看| 香蕉国产在线看| 黑人操中国人逼视频| 久久婷婷人人爽人人干人人爱| 操出白浆在线播放| 黑人巨大精品欧美一区二区mp4| 床上黄色一级片| 成人亚洲精品av一区二区| 欧美日韩黄片免| 欧美日韩乱码在线| 波多野结衣巨乳人妻| 色噜噜av男人的天堂激情| 久久久色成人| 超碰成人久久| 露出奶头的视频| 成人18禁在线播放| 女人高潮潮喷娇喘18禁视频| 禁无遮挡网站| 国产精品98久久久久久宅男小说| 18禁美女被吸乳视频| 怎么达到女性高潮| 久久久成人免费电影| aaaaa片日本免费| 亚洲中文字幕日韩| av视频在线观看入口| 窝窝影院91人妻| 男人舔女人下体高潮全视频| xxxwww97欧美| 老汉色∧v一级毛片| 欧美激情在线99| 免费一级毛片在线播放高清视频| 91麻豆精品激情在线观看国产| 日本黄色片子视频| 18禁观看日本| 黄片大片在线免费观看| h日本视频在线播放| 欧美日韩中文字幕国产精品一区二区三区| 高清在线国产一区| 国产一区在线观看成人免费| 99久久精品国产亚洲精品| 亚洲专区国产一区二区| 亚洲av第一区精品v没综合| 99热精品在线国产| 黄片小视频在线播放| 国产69精品久久久久777片 | 欧美在线一区亚洲| 一本综合久久免费| 90打野战视频偷拍视频| 怎么达到女性高潮| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| av中文乱码字幕在线| 1024手机看黄色片| 国语自产精品视频在线第100页| 人人妻人人澡欧美一区二区| 亚洲五月婷婷丁香| 黄色日韩在线| 69av精品久久久久久| 麻豆一二三区av精品| 精品国产美女av久久久久小说| 99久久精品国产亚洲精品| 不卡一级毛片| 99riav亚洲国产免费| 亚洲自偷自拍图片 自拍| 国产成人福利小说| www.熟女人妻精品国产| 中文字幕人成人乱码亚洲影| 夜夜夜夜夜久久久久| 成人av一区二区三区在线看| 又大又爽又粗| 成年版毛片免费区| 亚洲七黄色美女视频| 看片在线看免费视频| 丰满人妻一区二区三区视频av | 亚洲国产精品久久男人天堂| 97人妻精品一区二区三区麻豆| 欧美色视频一区免费| 又紧又爽又黄一区二区| 老鸭窝网址在线观看| 在线观看舔阴道视频| 久久精品影院6| 国产高清videossex| 亚洲欧美日韩无卡精品| 免费在线观看亚洲国产| 国产爱豆传媒在线观看| 亚洲欧美日韩高清专用| 我的老师免费观看完整版| 久久天躁狠狠躁夜夜2o2o| 久久精品亚洲精品国产色婷小说| 最近视频中文字幕2019在线8| 亚洲在线自拍视频| 岛国视频午夜一区免费看| 国产精品久久电影中文字幕| 巨乳人妻的诱惑在线观看| 99热6这里只有精品| 精品一区二区三区视频在线观看免费| 日韩高清综合在线| 久久国产精品影院| 午夜福利在线观看吧| 日本熟妇午夜| 香蕉丝袜av| 国产激情偷乱视频一区二区| 免费高清视频大片| 精品久久久久久久久久免费视频| 亚洲,欧美精品.| 蜜桃久久精品国产亚洲av| 村上凉子中文字幕在线| 老司机午夜十八禁免费视频| 91字幕亚洲| 两性夫妻黄色片| 国产精品1区2区在线观看.| 亚洲七黄色美女视频| 国产伦一二天堂av在线观看| 亚洲av第一区精品v没综合| 亚洲专区字幕在线| 搡老妇女老女人老熟妇| 国产精品国产高清国产av| 国产高清有码在线观看视频| 老熟妇乱子伦视频在线观看| 日本 欧美在线| 国产午夜精品论理片| 国产精品精品国产色婷婷| 中文字幕av在线有码专区| 琪琪午夜伦伦电影理论片6080| 日日夜夜操网爽| 十八禁网站免费在线| 国产高清三级在线| 国产高清视频在线观看网站| 啦啦啦韩国在线观看视频| 国产av不卡久久| 亚洲 欧美一区二区三区| 天堂√8在线中文| 看免费av毛片| 综合色av麻豆| 搡老熟女国产l中国老女人| 嫩草影院精品99| 欧美日本视频| 午夜亚洲福利在线播放| 亚洲va日本ⅴa欧美va伊人久久| 制服丝袜大香蕉在线| 一本综合久久免费| 成人午夜高清在线视频| 一个人看的www免费观看视频| 男女视频在线观看网站免费| 国产视频一区二区在线看| 俺也久久电影网| 日韩中文字幕欧美一区二区| 精品一区二区三区视频在线 | 国产三级黄色录像| 亚洲国产欧美网| 亚洲狠狠婷婷综合久久图片| 国产精品99久久久久久久久| 亚洲欧美日韩卡通动漫| 三级男女做爰猛烈吃奶摸视频| 亚洲精品久久国产高清桃花| 精品一区二区三区视频在线观看免费| 亚洲人成电影免费在线| 婷婷精品国产亚洲av在线| av片东京热男人的天堂| 国产精品久久久av美女十八| 午夜福利成人在线免费观看| 五月伊人婷婷丁香| 欧美国产日韩亚洲一区| 无遮挡黄片免费观看| 成人av在线播放网站| 熟妇人妻久久中文字幕3abv| 99久久无色码亚洲精品果冻| 亚洲欧美精品综合一区二区三区| 两个人视频免费观看高清| 免费无遮挡裸体视频| 少妇裸体淫交视频免费看高清| 一级毛片女人18水好多| 视频区欧美日本亚洲| 观看免费一级毛片| 女生性感内裤真人,穿戴方法视频| 国产精品永久免费网站| 一区福利在线观看| av黄色大香蕉| 国产精品久久久久久精品电影| 九色国产91popny在线| 国内揄拍国产精品人妻在线| 国产高清三级在线| 欧美午夜高清在线| 精品99又大又爽又粗少妇毛片 | 精品久久久久久久末码| 久久久久国产精品人妻aⅴ院| 天天躁日日操中文字幕| 老司机深夜福利视频在线观看| 欧美成人性av电影在线观看| 老司机深夜福利视频在线观看| 久久久久久久午夜电影| 亚洲自偷自拍图片 自拍| 久久久精品欧美日韩精品| 波多野结衣巨乳人妻| 97超视频在线观看视频| 真实男女啪啪啪动态图| 国产又色又爽无遮挡免费看| 母亲3免费完整高清在线观看| 无人区码免费观看不卡| 首页视频小说图片口味搜索| 伦理电影免费视频| 男人舔奶头视频| 成人av在线播放网站| 岛国在线观看网站| 啦啦啦免费观看视频1| 亚洲精品久久国产高清桃花| 又黄又粗又硬又大视频| 久久久久久久午夜电影| 熟女少妇亚洲综合色aaa.| 19禁男女啪啪无遮挡网站| 国产一区二区三区在线臀色熟女| 亚洲avbb在线观看| 超碰成人久久| 午夜福利欧美成人| 国产亚洲av嫩草精品影院| 欧美乱色亚洲激情| 男人舔女人下体高潮全视频| 日本一二三区视频观看| 麻豆成人av在线观看| 青草久久国产| 国产亚洲精品久久久久久毛片| 亚洲国产色片| 老熟妇仑乱视频hdxx| 男女视频在线观看网站免费| 熟女电影av网| 久久性视频一级片| 亚洲精品在线观看二区| 三级毛片av免费| av在线蜜桃| 亚洲欧美激情综合另类| 在线视频色国产色| 午夜久久久久精精品| 国产欧美日韩一区二区三| 99热精品在线国产| 99在线人妻在线中文字幕| 久久午夜综合久久蜜桃| 1024香蕉在线观看| 日本撒尿小便嘘嘘汇集6| 国产黄片美女视频| 亚洲熟妇熟女久久| 国产69精品久久久久777片 | 亚洲欧美日韩卡通动漫| 窝窝影院91人妻| 久久精品影院6| 成年免费大片在线观看| 亚洲中文字幕一区二区三区有码在线看 | 国产午夜精品久久久久久| 亚洲精品456在线播放app | 久久久久久久久免费视频了| 欧美黄色片欧美黄色片| 精品午夜福利视频在线观看一区| 日韩中文字幕欧美一区二区| 日韩av在线大香蕉| 国产熟女xx| 欧美黑人欧美精品刺激| 97超级碰碰碰精品色视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 男女做爰动态图高潮gif福利片| 嫁个100分男人电影在线观看| 成人一区二区视频在线观看| 伊人久久大香线蕉亚洲五| 制服丝袜大香蕉在线| 国产欧美日韩精品一区二区| 亚洲乱码一区二区免费版| 91在线精品国自产拍蜜月 | 国产三级黄色录像| 久久伊人香网站| www.自偷自拍.com| 69av精品久久久久久| 搡老熟女国产l中国老女人| 男女做爰动态图高潮gif福利片| 俺也久久电影网| 欧美绝顶高潮抽搐喷水| 中文字幕人妻丝袜一区二区| 亚洲欧美激情综合另类| 久久国产精品影院| 欧美在线一区亚洲| 精品久久久久久久人妻蜜臀av| 午夜福利欧美成人| 亚洲人成电影免费在线| 熟女电影av网| 观看美女的网站| 91老司机精品| 精华霜和精华液先用哪个| 午夜激情欧美在线| 亚洲成人免费电影在线观看| 嫩草影院精品99| 天堂网av新在线| 亚洲av美国av| 99久国产av精品| 悠悠久久av| 免费看十八禁软件| 成年版毛片免费区| 精品免费久久久久久久清纯| 啦啦啦韩国在线观看视频| 国产精品av视频在线免费观看| x7x7x7水蜜桃| av天堂中文字幕网| 日韩成人在线观看一区二区三区| 桃红色精品国产亚洲av| 国产午夜精品论理片| 日韩欧美三级三区| 亚洲欧洲精品一区二区精品久久久| 人人妻人人澡欧美一区二区| 18禁裸乳无遮挡免费网站照片| 一个人看视频在线观看www免费 | 亚洲 国产 在线| 青草久久国产| 亚洲国产高清在线一区二区三| 免费大片18禁| 亚洲av成人精品一区久久| www.999成人在线观看| 无限看片的www在线观看| 两个人视频免费观看高清| 两个人的视频大全免费| 午夜激情福利司机影院| 精品国产美女av久久久久小说| 日本一本二区三区精品| 制服丝袜大香蕉在线| 一个人看的www免费观看视频| 无人区码免费观看不卡| 久久国产乱子伦精品免费另类| 淫妇啪啪啪对白视频| 国产三级黄色录像| 国产精品1区2区在线观看.| 久久性视频一级片| 人妻夜夜爽99麻豆av| 久久久国产成人精品二区| 美女免费视频网站| 啦啦啦韩国在线观看视频| 90打野战视频偷拍视频| 亚洲人成网站在线播放欧美日韩| 久9热在线精品视频| 午夜福利欧美成人| 久久国产精品人妻蜜桃| 午夜亚洲福利在线播放| 久久国产精品影院| 欧美黄色淫秽网站| 国产精品一及| 亚洲欧美日韩高清在线视频| 亚洲aⅴ乱码一区二区在线播放| 一区二区三区激情视频| 99久国产av精品| 国产精品久久久久久精品电影| 亚洲无线观看免费| 国产高清有码在线观看视频| 日本黄大片高清| 国产成人aa在线观看| 最好的美女福利视频网| 亚洲av成人精品一区久久| 久久久国产精品麻豆| 熟女电影av网| 中文字幕精品亚洲无线码一区| 丰满人妻一区二区三区视频av | 最新美女视频免费是黄的| 日本精品一区二区三区蜜桃| 久久久久国产精品人妻aⅴ院| 丝袜人妻中文字幕| 一级a爱片免费观看的视频| 国产aⅴ精品一区二区三区波| 一夜夜www| 精品电影一区二区在线| 欧美大码av| 欧美成人性av电影在线观看| 欧美+亚洲+日韩+国产| 真实男女啪啪啪动态图| 99国产综合亚洲精品| 免费av不卡在线播放| 女人被狂操c到高潮| 搡老岳熟女国产| 亚洲av中文字字幕乱码综合| 亚洲黑人精品在线| 波多野结衣高清无吗| 真实男女啪啪啪动态图| cao死你这个sao货| 性欧美人与动物交配| 久久香蕉国产精品| 国产亚洲精品久久久久久毛片| 国产成人一区二区三区免费视频网站| 99热只有精品国产| 日韩欧美在线乱码| 真人一进一出gif抽搐免费| av欧美777|