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

    環(huán)太湖各水資源分區(qū)入出湖總氮負荷與濃度變化

    2023-06-09 04:16:10李順章李思昊陳彥彤余海翔王昱璋
    水利水電快報 2023年4期
    關(guān)鍵詞:湖西入湖藍藻

    李順章 李思昊 陳彥彤 余海翔 王昱璋

    摘要:

    為明確環(huán)太湖各水資源分區(qū)入出湖總氮負荷與濃度的時空變化規(guī)律,基于2012~2018年環(huán)太湖各水資源分區(qū)水質(zhì)水量監(jiān)測資料,計算分析了各水資源分區(qū)年、季(藍藻暴發(fā)期)尺度入出湖總氮負荷與濃度的年際變化及其負荷貢獻率時空變化規(guī)律。結(jié)果表明:環(huán)太湖入湖總氮多年平均濃度明顯高于出湖,各水資源分區(qū)年、季尺度入出湖總氮濃度的年際變化特征相似,總體而言,湖西區(qū)、浙西區(qū)、武澄錫虞區(qū)年、季尺度入湖總氮濃度明顯高于陽澄淀柳區(qū)和杭嘉湖區(qū),西太湖出湖高于東太湖出湖,入出湖水質(zhì)均未出現(xiàn)顯著改善。環(huán)太湖入湖總氮多年平均負荷遠大于出湖,各水資源分區(qū)年、季尺度入出湖總氮負荷年際變化規(guī)律相似,湖西區(qū)、浙西區(qū)在環(huán)太湖入湖負荷中的貢獻排名前兩位,在出湖負荷中太浦河的貢獻最大,藍藻暴發(fā)期武澄錫虞區(qū)入湖總氮負荷對全年的貢獻最大,藍藻暴發(fā)期太浦河和望虞河出湖總氮負荷對全年的貢獻受豐枯年份變化影響較大。

    關(guān)鍵詞:

    水資源分區(qū); 總氮; 水量加權(quán)平均法; 藍藻; 太湖

    中圖法分類號:X524

    文獻標志碼:A

    DOI:10.15974/j.cnki.slsdkb.2023.04.019

    文章編號:1006-0081(2023)04-0115-08

    0 引 言

    太湖位于中國長三角地區(qū),也是上海、蘇州、無錫等大中城市的重要水源地。太湖時常暴發(fā)不同強度的藍藻水華,嚴重威脅著供水安全,受到環(huán)太湖地區(qū)各級政府的高度重視和社會各界的廣泛關(guān)注[1]。現(xiàn)階段太湖營養(yǎng)過剩的狀況還未得到根本扭轉(zhuǎn),較高的氮磷營養(yǎng)鹽濃度導(dǎo)致藍藻水華頻發(fā)[2-3]。氮作為湖泊富營養(yǎng)化的重要控制因子之一,控制湖體總氮濃度一直都是治理太湖富營養(yǎng)化的重點[4-5]。程聲通等[6]研究表明,太湖總氮污染物主要來自入湖河道輸入,入湖河道總氮污染物輸入量超過太湖總?cè)牒摵傻?0%,因此,嚴格控制入湖氮負荷與濃度是解決太湖水體富營養(yǎng)化問題的關(guān)鍵[7],為此,亟需從不同時空尺度探明入出太湖總氮負荷與濃度的變化規(guī)律。已有學(xué)者開展了相關(guān)研究:呂文等[8]研究了太湖入湖總氮負荷的空間變化特征,認為西部、西北部湖區(qū)對應(yīng)的入湖河道污染物濃度較大,湖心、東部、東南部湖區(qū)濃度較低。馬倩等[9]探明了太湖入湖河道總氮濃度的年際變化特征,發(fā)現(xiàn)湖西區(qū)入湖河道總氮濃度與望虞河的總氮濃度之比為1.49~3.46,而望虞河入湖總氮濃度明顯小于湖西區(qū)水體總氮濃度,表明“引江濟太”工程調(diào)水引流的水質(zhì)明顯優(yōu)于太湖湖西區(qū)入湖河道水質(zhì)。呂文[8]、徐彬[10]等揭示了太湖入湖總氮負荷變化的影響因素,前者發(fā)現(xiàn)總氮入湖通量受水量影響,其低值和高值分別發(fā)生在2013年特枯年和2016年特豐年;后者的研究表明日均入湖水量與太湖氨氮自凈能力呈現(xiàn)出一定的正相關(guān)關(guān)系,長期引調(diào)水可以提高太湖的氨氮自凈能力。然而,值得一提的是,已有的相關(guān)研究較少以水資源分區(qū)為單位研究入出湖總氮負荷與濃度的變化,而且在時間尺度上也鮮少關(guān)注季尺度。因此,本文基于2012~2018年環(huán)太湖入出湖水量、水質(zhì)資料,從年、季尺度深入解析了環(huán)太湖各水資源分區(qū)入出湖河流總氮負荷與濃度變化。

    1 數(shù)據(jù)來源與研究方法

    1.1 數(shù)據(jù)來源

    本文中所采用的2012~2018年入出太湖水量資料依據(jù)江蘇省水文水資源勘測局提供的環(huán)太湖水文巡測數(shù)據(jù)。入出湖河流水量的計算參考李瓊芳等[11]的計算方法,規(guī)定入湖為正、出湖為負。本文中所采用的2012~2018年入出太湖主要河流的總氮資料分別由江蘇省水文水資源勘測局無錫分局、常州分局、蘇州分局提供,水質(zhì)監(jiān)測斷面布設(shè)見圖1。巡測段、站歸類見表1。

    1.2 研究方法

    參考文獻[11],將環(huán)太湖區(qū)域劃分浙西區(qū)、湖西區(qū)、武澄錫虞區(qū)、陽澄淀泖區(qū)、杭嘉湖區(qū)等5個水資源分區(qū)以及望虞河、太浦河兩條主要河流(圖1),并采用文獻中污染物負荷和濃度計算方法,計算不同年份環(huán)太湖和不同水資源分區(qū)的入出湖總氮負荷與濃度,分析其年際變化規(guī)律;計算不同年份各水資源分區(qū)入出湖負荷在總?cè)氤龊恼急?,解析不同水資源分區(qū)總氮貢獻的年際變化;計算不同年份各水資源分區(qū)藍藻暴發(fā)期(5~10月)入出湖總氮負荷與濃度以及其入出湖負荷在全年入出湖的占比,探明藍藻暴發(fā)期各水資源分區(qū)入出湖總氮負荷與濃度以及對全年貢獻的年際變化。

    2 結(jié)果分析

    2.1 年尺度

    2.1.1 環(huán)太湖與不同水資源分區(qū)入出湖總氮負荷變化特征

    環(huán)太湖和不同水資源分區(qū)入湖總氮負荷變化過程見圖2(a)。從圖2(a)可以看出,環(huán)太湖和不同分區(qū)入湖總氮年負荷雖均隨年份變化而變化,但其年際變化特征不盡相同。環(huán)太湖、湖西區(qū)、浙西區(qū)總氮負荷年際變化較其他水資源分區(qū)的

    更加顯著,而且湖西區(qū)和浙西區(qū)的各年入湖總氮負荷要遠高于其他分區(qū),尤其是湖西區(qū),更為突出。環(huán)太湖的入湖總氮負荷呈現(xiàn)出的年際變化規(guī)律與湖西區(qū)、浙西區(qū)相似,其年負荷第一大來源是湖西區(qū),其次是浙西區(qū)。2012~2018年期間,環(huán)太湖、湖西區(qū)和浙西區(qū)的最大入湖年負荷值均出現(xiàn)在2016年,分別為53 278.10,34 983.30 t和14 729.20 t,2016年流域發(fā)生特大洪水,湖西區(qū)和浙西區(qū)的入湖水量顯著增加[12],導(dǎo)致這兩區(qū)進入太湖的總氮負荷也顯著增加;環(huán)太湖和湖西區(qū)的最小入湖年負荷值均出現(xiàn)在2013年,分別為33 600.16 t和22 752.05 t,而其他水資源分區(qū)卻并非如此,這進一步表明環(huán)太湖負荷年際變化規(guī)律在很大程度上取決于湖西區(qū)的變化。對于具有引江濟太功能的望虞河,通過其流入太湖的總氮年負荷在各水資源分區(qū)中是明顯偏小的,表明各年由望虞河帶入的總氮負荷非常有限。

    環(huán)太湖和不同水資源分區(qū)出湖總氮負荷變化過程見圖2(b)。從圖2(b)可以看出,環(huán)太湖和不同分區(qū)出湖總氮負荷呈現(xiàn)出不一致的變化規(guī)律。環(huán)太湖、太浦河流出太湖的總氮負荷年際變化明顯比其他水資源分區(qū)顯著,而且兩者的波動規(guī)律具有高度一致性。從太浦河流出太湖的年總氮負荷除了在2013年低于浙西區(qū)、陽澄淀泖區(qū)、杭嘉湖區(qū),在其他年份均高于從其他分區(qū)流出太湖的年總氮負荷。2012~2018年期間,環(huán)太湖與太浦河的最大和最小出湖年負荷值(28 960.99 ,12 701.36 t和15 375.16 ,1 463.63 t)也分別出現(xiàn)在2016年和2013年,表明環(huán)太湖負荷年際變化規(guī)律在很大程度上取決于從太浦河流出的總氮變化。對于同時具有引排功能的望虞河,其年際波動特征與太浦河相似,其最大值與最小值出現(xiàn)的年份與太浦河完全一致。在流域降水量較小的年份,其排出太湖的總氮負荷也比較小,但在大水年份,其排出太湖的總氮負荷明顯增加,如2016年,其出湖負荷僅低于太浦河。表明望虞河在大水年份通過排水有利于降低湖區(qū)總氮負荷。

    環(huán)太湖、太浦河最大和最小出湖年負荷值出現(xiàn)的時間與環(huán)太湖、湖西區(qū)最大和最小入湖年負荷值出現(xiàn)的時間一致。根據(jù)林鵬等[12]研究,2012~2018年期間入出湖年水量最大和最小值分別發(fā)生在2016年和2013年,表明入出湖水量的大小對入出湖總氮負荷具有決定性的作用:對于入湖負荷,湖西區(qū)入湖水量具有決定性作用;對于出湖負荷,太浦河出湖水量具有決定性作用。對比環(huán)太湖入出湖總氮負荷,發(fā)現(xiàn)在2012~2018年間每年由入湖水流攜帶進入太湖的總氮負荷始終高于出湖水流帶出的總氮負荷。

    2.1.2 不同水資源分區(qū)入出湖總氮負荷對環(huán)太湖的貢獻變化

    不同水資源分區(qū)入湖總氮負荷對環(huán)太湖的貢獻變化過程見圖3(a)。從圖3(a)可以看出,不同分區(qū)入湖總氮負荷對環(huán)太湖的貢獻率雖均隨年份變化而變化,但總體而言變化幅度不大。在所有年份中湖西區(qū)的貢獻最大,平均而言,為70%~80%;浙西區(qū)的貢獻位居第二,貢獻率的平均值為20%~30%;其他水資源分區(qū)各年的貢獻率明顯小于湖西區(qū)和浙西區(qū),均低于10%。2012~2018年期間,雖然湖西區(qū)2016年的入湖總氮負荷最大,但該年湖西區(qū)對環(huán)太湖的貢獻率卻低于其他年份,表明2016年其他水資源分區(qū)的貢獻有所增加;浙西區(qū)最大年總氮負荷和貢獻率出現(xiàn)的時間一致,均為2016年。湖西區(qū)和浙西區(qū)兩區(qū)一共超過90%的貢獻率表明,削減這兩區(qū)的入湖總氮負荷是降低湖體氮負荷的關(guān)鍵。對于望虞河,通過其流入太湖的總氮年負荷對環(huán)太湖的貢獻明顯偏小,特別是2016年大水年,由望虞河帶入的總氮負荷的貢獻更是偏低。

    不同水資源分區(qū)出湖總氮負荷對環(huán)太湖的貢獻變化過程見圖3(b)。從圖3(b)可以看出,與不同分區(qū)入湖總氮負荷對環(huán)太湖的貢獻率年際變化相比,不同分區(qū)出湖總氮負荷對環(huán)太湖的貢獻率年際變化明顯加劇,其中變化最劇烈的是太浦河,最大值(43.86%)出現(xiàn)在2016年,最小值(9.52%)在2013年;其次是浙西區(qū),貢獻率在8%~33%之間變化;再者是杭嘉湖區(qū),其貢獻率的變化范圍為10%~30%;排名第4的是望虞河,2016年貢獻率最大(19.01%),2013年最?。?.35%);剩余的其他兩個水資源分區(qū)貢獻率變化相對平穩(wěn)得多,其中武澄錫虞區(qū)各年的貢獻率總體明顯偏小。2012~2018年期間,太浦河和望虞河的貢獻率最大值和最小值出現(xiàn)的年份一致,但浙西區(qū)和杭嘉湖區(qū)的貢獻率最大值和最小值出現(xiàn)的年份正好與太浦河和望虞河相反。表明太浦河和望虞河作為重要排水通道,在大水年2016年將水排出太湖的同時也增加了排出總氮負荷的貢獻率。在小水年份2013年,浙西區(qū)和杭嘉湖區(qū)的出湖水量增加,也導(dǎo)致兩區(qū)隨水流排出的總氮負荷貢獻率最大。由此可見,不同水資源分區(qū)年出湖總氮負荷貢獻率的變化與來水情況緊密相關(guān)。

    2.1.3 環(huán)太湖和不同水資源分區(qū)入出湖總氮濃度變化特征

    不同水資源分區(qū)河流入出湖總氮平均濃度變化如圖4所示。不同水資源分區(qū)河流入湖總氮平均濃度變化如圖4(a)所示。從圖4(a)可以看出,所有水資源分區(qū)各年入湖水質(zhì)均較差,屬于Ⅴ類或劣Ⅴ類,其中武澄錫虞區(qū)的入湖水質(zhì)總體而言最差,杭嘉湖區(qū)的最好。從年際變化來看,除了杭嘉湖區(qū)的入湖水質(zhì)未呈現(xiàn)明顯改善外,其他水資源分區(qū)和環(huán)太湖均呈現(xiàn)出朝水質(zhì)改善的方面發(fā)展,但改善效果仍有待大力提升。特別是入湖水量最大和次大的湖西區(qū)和浙西區(qū),其總氮濃度一直維持在劣Ⅴ類,水質(zhì)改善工作任重道遠。湖西區(qū)和浙西區(qū)入湖水量大,再疊加其水質(zhì)差,導(dǎo)致環(huán)太湖的平均總氮濃度也一致保持在劣Ⅴ類。從變化幅度來說,武澄錫虞區(qū)和陽澄淀泖區(qū)明顯比環(huán)太湖和其他水資源分區(qū)顯著。這兩區(qū)在2014年出現(xiàn)特高值,這是因為:2014年武澄錫虞區(qū)7~8月入湖水量大,而其他月份入湖水量為零,再加上主要入湖河道武進港7月總氮濃度出現(xiàn)特高值(6.66 mg/L),導(dǎo)致水量加權(quán)法計算得到的總氮平均濃度大;陽澄淀泖區(qū)由于其主要入湖河流滸光運河4月總氮濃度出現(xiàn)特高值,同時該月降水較常年同期多45%,入湖水量驟增,最終導(dǎo)致水量加權(quán)平均的年入湖總氮濃度較高。在空間分布上,不同分區(qū)2012~2018年入湖總氮多年平均濃度由高至低排序依次為湖西區(qū)(3.79 mg/L)>浙西區(qū)(3.62 mg/L)>武澄錫虞區(qū)(3.59 mg/L)>陽澄淀泖區(qū)(2.57 mg/L)>望虞河(2.41 mg/L)>杭嘉湖區(qū)(2.30 mg/L)。從入湖河道總氮濃度的空間分布差異來看,湖西區(qū)、浙西區(qū)以及武澄錫虞區(qū)對應(yīng)的太湖西部的總氮濃度高于太湖東部水資源分區(qū),與文獻[8]的分析結(jié)果一致。

    不同水資源分區(qū)河流出湖總氮平均濃度變化如圖4(b)所示。從圖4(b)可以看出,與入湖水質(zhì)相比,總體而言出湖水質(zhì)明顯偏好。除了浙西區(qū)各年出湖水質(zhì)屬于劣Ⅴ類,其他水資源分區(qū)各年出湖水質(zhì)在Ⅲ~Ⅴ類之間變化。其中武澄錫虞區(qū)的出湖水質(zhì)最好,陽澄淀泖區(qū)僅次于武澄錫虞區(qū)。從年際變化來看,除了陽澄淀泖區(qū)的出湖水質(zhì)有所變差外,其他水資源分區(qū)和環(huán)太湖的出湖濃度變化都比較穩(wěn)定。從變化幅度來說,武澄錫虞區(qū)和陽澄淀泖區(qū)明顯比環(huán)太湖和其他水資源分區(qū)的顯著。在空間分布上,不同分區(qū)2012~2018年出湖總氮多年平均濃度由高至低排序依次為浙西區(qū)(3.21 mg/L)>杭嘉湖區(qū)(2.00 mg/L)>太浦河(1.92 mg/L)>望虞河(1.68 mg/L)>陽澄淀泖區(qū)(1.23 mg/L)>武澄錫虞區(qū)(0.92 mg/L)。

    2.2 季尺度

    2.2.1 不同水資源分區(qū)藍藻暴發(fā)期入出湖總氮負荷變化特征

    不同水資源分區(qū)藍藻暴發(fā)期(5~10月)入出湖總氮負荷變化過程見圖5。不同水資源分區(qū)藍藻暴發(fā)期(5~10月)入湖總氮負荷變化過程見圖5(a)。從圖5(a)可以看出,湖西區(qū)和浙西區(qū)藍藻暴發(fā)期總氮負荷年際變化呈現(xiàn)出了與其年尺度的相似特征:年際變化遠較其他分區(qū)顯著,且其值遠高于其他分區(qū),湖西區(qū)尤其突出;2012~2018年期間最大值和最小值出現(xiàn)的時間也與年尺度的一致。表明這兩個區(qū)的年尺度總氮負荷在很大程度上取決于其藍藻暴發(fā)期的總氮負荷。對于望虞河,除了2013年小水年份在藍藻暴發(fā)期通過加大從長江向太湖引水導(dǎo)致其總氮負荷有所增加外,其他年份的總氮負荷都很小。

    不同水資源分區(qū)藍藻暴發(fā)期(5~10月)出湖總氮負荷變化過程見圖5(b)。從圖5(b)可以看出,藍藻暴發(fā)期內(nèi)不同水資源分區(qū)出湖總氮負荷同樣呈現(xiàn)出不一致的變化規(guī)律。太浦河的年際變化最為顯著,其次是望虞河,再者是陽澄淀泖區(qū),而浙西區(qū)、杭嘉湖區(qū)和武澄錫虞區(qū)的年際變化相對穩(wěn)定。太浦河和望虞河藍藻暴發(fā)期總氮負荷呈現(xiàn)出與年尺度相似的年際變化規(guī)律:太浦河和望虞河藍藻暴發(fā)期總氮負荷的年際波動一致。2012~2014年期間太浦河藍藻暴發(fā)期出湖的總氮負荷低于浙西區(qū)和杭嘉湖區(qū),但在其他年份均高于其余分區(qū)出湖。太浦河和望虞河作為主要排水通道,2013年小水年排出太湖的總氮負荷最?。ǚ謩e為1 126.68 t和332.10 t),但在2016年大水年,兩者排出太湖的總氮負荷最大(分別為7 248.17 t和3 917.38 t),且該年望虞河僅次于太浦河。太浦河和望虞河的年際變化特征表明,其出湖總氮負荷的大小與不同來水量緊密相關(guān)。此外,望虞河藍藻暴發(fā)期出湖的總氮負荷在一般年份通常比浙西區(qū)、杭嘉湖區(qū)和陽澄淀泖區(qū)小。對比藍藻暴發(fā)期入出湖總氮負荷可以看出,在2012~2018年間每年由入湖水流攜帶進入太湖的總氮負荷始終高于出湖水流帶出的總氮負荷。

    2.2.2 不同水資源分區(qū)藍藻暴發(fā)期入出湖總氮負荷對年負荷的貢獻變化

    不同水資源分區(qū)藍藻暴發(fā)期入出湖總氮負荷對年負荷的貢獻變化過程見圖6。不同水資源分區(qū)藍藻暴發(fā)期入湖總氮負荷對年負荷的貢獻變化過程見圖6(a)。從圖6(a)可以看出,不同分區(qū)藍藻暴發(fā)期入湖總氮負荷對年負荷貢獻率的年際變化存在明顯差異:望虞河最顯著,浙西區(qū)、杭嘉湖區(qū)和陽澄淀泖區(qū)的年際變化顯著程度相似,武澄錫虞區(qū)和湖西區(qū)相對穩(wěn)定,武澄錫虞區(qū)各年的貢獻率均高于其他分區(qū)。藍藻暴發(fā)期入湖總氮負荷最大的湖西區(qū)對全年負荷的貢獻率維持在50%左右,表明非藍藻暴發(fā)期入湖總氮負荷對全年負荷的貢獻與藍藻暴發(fā)期的同等重要。杭嘉湖區(qū)、陽澄淀泖區(qū)、浙西區(qū)和望虞河藍藻暴發(fā)期入湖總氮負荷對年負荷貢獻率最大值均發(fā)生在2013年枯水年,說明2013年這些分區(qū)其他月份入湖負荷量有所減少。常規(guī)年份武澄錫虞區(qū)在非藍藻暴發(fā)期閘門通常是關(guān)閉的,幾乎沒有總氮負荷流入太湖,僅在2016年大水年少量總氮負荷流入太湖。在2016年大水年,杭嘉湖區(qū)、陽澄淀泖區(qū)、浙西區(qū)、湖西區(qū)的貢獻率為42%~73%,但望虞河小于30%,表明大水年藍藻暴發(fā)期望虞河的貢獻率明顯低于枯水年份。2014年和2015年望虞河藍藻暴發(fā)期從長江向太湖引水非常有限,其貢獻率基本為零。

    不同水資源分區(qū)藍藻暴發(fā)期(5~10月)出湖總氮負荷對年負荷貢獻的變化過程見圖6(b)。從圖6(b)可以看出,兩個主要排水通道太浦河和望虞河的貢獻率年際變化最顯著,其次是武澄錫虞區(qū),剩余3個水資源分區(qū)的年際變化比較平穩(wěn)。在常規(guī)年份,太浦河藍藻暴發(fā)期出湖總氮負荷對全年負荷的貢獻率通常都小于50%,表明非藍藻暴發(fā)期的出湖負荷貢獻占主導(dǎo)作用;在2013年枯水年,其在藍藻暴發(fā)期的貢獻率達到最大(76.98%),表明非藍藻暴發(fā)期的出湖負荷貢獻隨出湖水量的減少明顯降低;在2016年大水年,其在藍藻暴發(fā)期的貢獻率接近60%,表明藍藻暴發(fā)期的出湖負荷貢獻隨出湖水量的增加有所增加。與太浦河相反,望虞河藍藻暴發(fā)期出湖總氮負荷對全年負荷的貢獻率除了2013年在50%左右外,其他年份都明顯高于50%,表明在非小水年份望虞河在藍藻暴發(fā)期具有重要的排水功能,其在藍藻暴發(fā)期排水量的增加導(dǎo)致其出湖總氮負荷對年出湖負荷的貢獻率也提高。武澄錫虞區(qū)在藍藻暴發(fā)期的貢獻率在2016年最高,表明其在2016年大水年通過增加藍藻暴發(fā)期排水量從而提高了貢獻率。剩余3個水資源分區(qū)各年的貢獻率在50%左右變化,表明不同來水年對其貢獻率影響不大,藍藻暴發(fā)期和非藍藻暴發(fā)期對全年出湖負荷的影響比較均衡。

    2.2.3 不同水資源分區(qū)藍藻暴發(fā)期入出湖總氮濃度的變化

    不同水資源分區(qū)藍藻暴發(fā)期入出湖總氮平均濃度變化如圖7所示。不同水資源分區(qū)藍藻暴發(fā)期入湖總氮平均濃度變化如圖7(a)所示。從圖7(a)可以看出,所有水資源分區(qū)各年藍藻暴發(fā)期入湖水質(zhì)均較差,屬于Ⅳ類至劣Ⅴ類,其中武澄錫虞區(qū)的水質(zhì)總體而言最差,望虞河最好。從年際變化來看,武澄錫虞區(qū)和陽澄淀泖區(qū)的波動最大,其他水資源分區(qū)的變化相對穩(wěn)定。值得一提的是,總氮負荷最大的湖西區(qū)和浙西區(qū)的濃度呈現(xiàn)穩(wěn)中有降的趨勢,但仍在劣Ⅴ類范圍,水質(zhì)提升仍然面臨嚴峻挑戰(zhàn)。自2015年,陽澄淀泖區(qū)的水質(zhì)有所改善,但由于其對全湖的負荷貢獻偏小,對降低入湖總氮負荷貢獻有限。與其他水資源分區(qū)的水質(zhì)相比,望虞河和杭嘉湖區(qū)的水質(zhì)雖然偏優(yōu),但仍有提升空間。由于望虞河2015年藍藻暴發(fā)期內(nèi)沒有入湖水量,圖7(a)中望虞河的總氮濃度實際上是實測的總氮濃度的算術(shù)平均值(1.75 mg/L)。就不同水資源分區(qū)2012~2018年入湖總氮多年平均濃度而言,由高至低依次排序為武澄錫虞區(qū)(3.62 mg/L)>浙西區(qū)(3.37 mg/L)>湖西區(qū)(3.18 mg/L)>陽澄淀泖區(qū)(2.26 mg/L)>杭嘉湖區(qū)(2.17 mg/L)>望虞河(2.06 mg/L)。

    不同水資源分區(qū)藍藻暴發(fā)期出湖總氮平均濃度變化如圖7(b)所示。從圖7(b)可以看出,與入湖水質(zhì)相比,總體而言藍藻暴發(fā)期出湖水質(zhì)明顯偏好。除了浙西區(qū)和杭嘉湖區(qū)各年藍藻暴發(fā)期出湖水質(zhì)屬于劣Ⅴ類,其他水資源分區(qū)的出湖水質(zhì)在Ⅲ類至Ⅴ類之間變化,其中武澄錫虞區(qū)最好。在年際變化方面,除了浙西區(qū)的波動稍大一些,其他各區(qū)的變化相對穩(wěn)定,但都沒有改善的趨勢。望虞河藍藻暴發(fā)期出湖總氮濃度基本維持在Ⅴ類水質(zhì)。太浦河的出湖總氮濃度總體上要低于望虞河但高于陽澄淀泖區(qū)。藍藻暴發(fā)期不同分區(qū)的2012~2018年出湖總氮多年平均濃度由高至低依次排序為浙西區(qū)(2.96 mg/L)>杭嘉湖區(qū)(1.99 mg/L)>望虞河(1.56 mg/L)>太浦河(1.49 mg/L)>陽澄淀泖區(qū)(1.27 mg/L)>武澄錫虞區(qū)(0.82 mg/L)。

    3 結(jié) 論

    本文研究了2012~2018年環(huán)太湖各水資源分區(qū)入出湖總氮負荷與濃度變化,得到主要結(jié)論如下:

    (1) 2012~2018年間環(huán)太湖入湖總氮多年平均濃度(3.62 mg/L)明顯高于出湖的總氮多年平均濃度(1.78 mg/L);湖西區(qū)、浙西區(qū)以及武澄錫虞區(qū)入湖總氮濃度明顯高于太湖東南部的陽澄淀泖區(qū)和杭嘉湖區(qū),西太湖出湖總氮濃度高于東太湖;與環(huán)太湖入湖總氮多年平均濃度相比,望虞河入湖總氮多年平均濃度(2.41 mg/L)明顯偏低。

    (2) 入湖總氮多年平均負荷(42 021.66 t)遠大于出湖總氮多年平均負荷值(19 488.19 t),且各年環(huán)太湖入湖總氮負荷始終大于出湖;各年湖西區(qū)入湖總氮負荷對環(huán)太湖總?cè)牒摵傻呢暙I最大,占比高達70%以上,其次是浙西區(qū),占比在13%~28%之間變化;在出湖總氮負荷中,平均而言太浦河的貢獻最大,占比達到32.10%;望虞河出入湖總氮負荷在豐枯年份差異較大,在2013年枯水年,其入湖的總氮負荷比其他年份的略高,而在2016年豐水年,其出湖的總氮負荷則明顯高于其他年份。

    (3) 藍藻暴發(fā)期入湖總氮負荷所呈現(xiàn)出的特征與年尺度相似,其中湖西區(qū)與浙西區(qū)排名前兩位,但藍藻暴發(fā)期湖西區(qū)入湖負荷對全年的貢獻并不突出,而武澄錫虞區(qū)對全年的貢獻率最大。藍藻暴發(fā)期太浦河和望虞河出湖總氮負荷對全年的貢獻率年際變化較大,主要受年際水量變化的影響。

    (4) 藍藻暴發(fā)期入出湖總氮平均濃度變化特征與年尺度相似,其中湖西區(qū)和浙西區(qū)入湖水質(zhì)有所改善,但依然較差,其余水資源分區(qū)入湖水質(zhì)改善不明顯;所有出湖水質(zhì)雖隨年份波動,但沒有明顯改善勢頭。

    參考文獻:

    [1] 房振南,金科,王雪姣,等.長三角一體化生態(tài)綠色發(fā)展示范區(qū)主要河湖水質(zhì)變化趨勢分析[J].水利水電快報,2021,42(4):68-74.

    [2] 朱偉,談永琴,王若辰,等.太湖典型區(qū)2010-2017年間水質(zhì)變化趨勢及異常分析[J].湖泊科學(xué),2018,30(2):296-305.

    [3] 蔡梅,王元元,龔李莉,等.新一輪太湖流域水環(huán)境綜合治理的思考和建議[J].水利規(guī)劃與設(shè)計,2021(2):4-6,45.

    [4] 代倩子,張坤,徐彬.2014-2018年太湖梅梁湖和東太湖水質(zhì)變化特征分析[J].中國農(nóng)村水利水電,2020(7):82-84.

    [5] 秦伯強.淺水湖泊湖沼學(xué)與太湖富營養(yǎng)化控制研究[J].湖泊科學(xué),2020,32(5):1229-1243.

    [6] 程聲通,錢益春,張紅舉.太湖總磷、總氮宏觀水環(huán)境容量的估算與應(yīng)用[J].環(huán)境科學(xué)學(xué)報,2013,33(10):2848-2855.

    [7] 秦伯強,楊柳燕,陳非洲,等.湖泊富營養(yǎng)化發(fā)生機制與控制技術(shù)及其應(yīng)用[J].科學(xué)通報,2006(16):1857-1866.

    [8] 呂文,楊惠,楊金艷,等.環(huán)太湖江蘇段入湖河道污染物通量與湖區(qū)水質(zhì)的響應(yīng)關(guān)系[J].湖泊科學(xué),2020,32(5):1454-1462.

    [9] 馬倩,田威,吳朝明.望虞河引長江水入太湖水體的總磷、總氮分析[J].湖泊科學(xué),2014,26(2):207-212.

    [10] 徐彬,顧蘇莉,毛新偉.引江濟太調(diào)水與太湖氨氮自凈能力的相關(guān)性分析[C]∥中國水利學(xué)會.科技創(chuàng)新與水利改革——中國水利學(xué)會2014學(xué)術(shù)年會論文集(上冊).北京:中國水利學(xué)會,2014.

    [11] 李瓊芳,許樹洪,陳啟慧,等.環(huán)太湖各水資源分區(qū)入出湖河流總磷濃度與負荷變化分析[J].湖泊科學(xué),2022,34(1):74-89.

    [12] 林鵬,陳啟慧,李瓊芳,等.環(huán)太湖各水資源分區(qū)入出湖水量及貢獻分析[J].水資源保護,2021,37(3):66-73.

    (編輯:江 文)

    Variation of inflow and outflow TN loads and concentrations of different water resources zones around Taihu Lake

    LI Shunzhang,LI Sihao,CHEN Yantong,YU Haixiang ,WANG Yuzhang

    (Dayu College,Hohai University,Nanjing 210098,China)Abstract:

    To clarify the spacial-temporal variations of the total nitrogen(TN) loads and concentrations of Taihu Lake,based on the observed water quality and quantity data from 2012 to 2018,the TN loads and flow-weighted mean concentrations into and out of Taihu Lake were calculated,the annual and seasonal (cyanobacteria outbreak period) scale TN loads and concentrations of inflows and outflows and their contribution rates were calculated by water resources zones,then the paper analyzed the spacial-temporal variations of the TN loads and concentrations.The results showed that the average annual TN concentrations of the inflow rivers were significantly higher than those of the outflow rivers,and the inter-annual variation of TN concentrations of the inflow and outflow rivers were similar in all water resource zones.The inter-annual average TN load of the inflow rivers was much larger than those out of the lake,and the interannual variation patterns of annual and seasonal loads of TN into and out of the lake were similar in all water resource zones,the annual and seasonal TN concentrations of inflow rivers from Huxi,Zhexi and Wuchengxiyu zone were significantly higher than those from Yangchengdianmao zone and Hangjiahu zone,and the TN concentrations of flows out of west lake was higher than those from east lake,and the water quality of river around the lake was not significantly improved;The annual average loads of TN flowing into the lake was much greater than that flowing out of the lake,and the annual and seasonal variation patterns of TN loads of inflow and outflow rivers were similar in each water resources zone.The proportion of TN loads of inflow rivers from Huxi and Zhexi zone was first two among all zones,and the proportion of TN load from Taipu River was the largest.During the cyanobacteria outbreak period,The contribution of TN load of inflows during the cyanobacteria outbreak period to the whole annual TN load from Wuchengxiyu zone was the largest.The contribution of TN load from Taipu and Wangyu river during the cyanobacteria outbreak period to the whole annual TN load were greatly influenced by the variation of abundant and dry years.

    Key words:

    water resources zone; total nitrogen; weighted average method of water quantity; cyanophyta; Taihu Lake

    收稿日期:

    2022-05-01

    基金項目:

    國家自然科學(xué)基金面上項目(51879069);大學(xué)生創(chuàng)新訓(xùn)練基金項目(202110294045);江蘇省水利科技項目(2018001)

    作者簡介:

    李順章,男,研究方向為水文與水資源工程。E-mail:2824119941@qq.com

    猜你喜歡
    湖西入湖藍藻
    降雨事件下洱海北部主要入湖河流污染物特征分析*
    新時代湖西紅色文化傳承與發(fā)展創(chuàng)新研究*
    環(huán)洱海主要入湖河流水質(zhì)特征及入湖污染負荷估算
    人民長江(2022年1期)2022-04-01 02:44:35
    龍湖·首開湖西星辰美學(xué)館
    南美白對蝦養(yǎng)殖池塘藍藻水華處理舉措
    針對八月高溫藍藻爆發(fā)的有效處理方案
    入湖河口濕地恢復(fù)與重建規(guī)劃設(shè)計初探——以資興市興寧河入湖河口濕地為例
    滇池入湖河流磷負荷時空變化及形態(tài)組成貢獻*
    可怕的藍藻
    巢湖滅藻
    民生周刊(2012年25期)2012-04-29 00:44:03
    色婷婷久久久亚洲欧美| 变态另类成人亚洲欧美熟女 | 变态另类成人亚洲欧美熟女 | 免费搜索国产男女视频| av在线天堂中文字幕| 大型av网站在线播放| 69精品国产乱码久久久| 欧美国产日韩亚洲一区| 欧美性长视频在线观看| 少妇 在线观看| 悠悠久久av| 国产成人系列免费观看| 国产av一区二区精品久久| 国产一卡二卡三卡精品| 日本免费a在线| 黄色片一级片一级黄色片| 久久人人爽av亚洲精品天堂| 亚洲中文字幕一区二区三区有码在线看 | 国产成人系列免费观看| 免费女性裸体啪啪无遮挡网站| 日韩有码中文字幕| 亚洲天堂国产精品一区在线| 搞女人的毛片| 亚洲午夜理论影院| 亚洲成av片中文字幕在线观看| 欧美黑人精品巨大| 日本欧美视频一区| 午夜a级毛片| 嫩草影视91久久| 中文字幕色久视频| 久久人妻福利社区极品人妻图片| 国产精品爽爽va在线观看网站 | 欧美+亚洲+日韩+国产| 美女午夜性视频免费| 国产精品二区激情视频| 欧美黑人欧美精品刺激| 高清在线国产一区| 精品无人区乱码1区二区| 麻豆国产av国片精品| 中国美女看黄片| 性色av乱码一区二区三区2| 久久久久久久午夜电影| 国内毛片毛片毛片毛片毛片| АⅤ资源中文在线天堂| 欧美日本中文国产一区发布| 一进一出好大好爽视频| 黄色女人牲交| 19禁男女啪啪无遮挡网站| 国产亚洲精品久久久久5区| 最近最新免费中文字幕在线| 很黄的视频免费| 真人做人爱边吃奶动态| 精品国产亚洲在线| 黄色视频不卡| 操美女的视频在线观看| 欧美乱色亚洲激情| 一区福利在线观看| a在线观看视频网站| 成人18禁在线播放| 欧美成人午夜精品| 99在线人妻在线中文字幕| 后天国语完整版免费观看| 两个人看的免费小视频| 午夜精品国产一区二区电影| 成人亚洲精品一区在线观看| 桃色一区二区三区在线观看| 亚洲最大成人中文| 亚洲视频免费观看视频| 99在线视频只有这里精品首页| 中文字幕最新亚洲高清| 老司机靠b影院| 欧美黑人欧美精品刺激| 免费在线观看亚洲国产| 国产熟女xx| 老熟妇仑乱视频hdxx| 久久中文看片网| 18禁裸乳无遮挡免费网站照片 | 香蕉国产在线看| 久久人人爽av亚洲精品天堂| 久久精品91蜜桃| 久久中文看片网| 久久久久九九精品影院| 日日爽夜夜爽网站| 91字幕亚洲| АⅤ资源中文在线天堂| 最新在线观看一区二区三区| 女性被躁到高潮视频| 亚洲第一av免费看| 午夜亚洲福利在线播放| 欧美日韩亚洲国产一区二区在线观看| 日本黄色视频三级网站网址| 女性生殖器流出的白浆| 国产成人欧美在线观看| 一二三四在线观看免费中文在| 国产精品久久久人人做人人爽| 一级,二级,三级黄色视频| svipshipincom国产片| 国产乱人伦免费视频| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 中文字幕色久视频| 亚洲熟女毛片儿| 99热只有精品国产| 三级毛片av免费| 国产精品野战在线观看| 精品不卡国产一区二区三区| 精品无人区乱码1区二区| 亚洲五月婷婷丁香| 黄片大片在线免费观看| 人人妻人人澡人人看| 极品人妻少妇av视频| 久久草成人影院| 国产高清videossex| 嫩草影视91久久| 国产一卡二卡三卡精品| 久久精品国产亚洲av高清一级| 一级片免费观看大全| 久久久国产成人精品二区| 99精品在免费线老司机午夜| 欧美国产精品va在线观看不卡| 亚洲成a人片在线一区二区| 成在线人永久免费视频| 又大又爽又粗| bbb黄色大片| 免费少妇av软件| 好男人电影高清在线观看| 午夜福利18| 日本 av在线| 久久天堂一区二区三区四区| 国产精品日韩av在线免费观看 | 丰满人妻熟妇乱又伦精品不卡| 热re99久久国产66热| 少妇 在线观看| 国产又色又爽无遮挡免费看| 亚洲熟妇熟女久久| or卡值多少钱| 亚洲一区中文字幕在线| 女人爽到高潮嗷嗷叫在线视频| av免费在线观看网站| 亚洲 欧美 日韩 在线 免费| 久久草成人影院| 亚洲国产中文字幕在线视频| 久久热在线av| 夜夜躁狠狠躁天天躁| 一边摸一边抽搐一进一出视频| 精品电影一区二区在线| 国产精品免费一区二区三区在线| 免费观看人在逋| 亚洲一卡2卡3卡4卡5卡精品中文| 伦理电影免费视频| 性色av乱码一区二区三区2| 天天躁狠狠躁夜夜躁狠狠躁| 不卡av一区二区三区| 欧美成人一区二区免费高清观看 | av网站免费在线观看视频| 久久天躁狠狠躁夜夜2o2o| 亚洲精华国产精华精| 成人三级黄色视频| 嫩草影视91久久| 男女下面进入的视频免费午夜 | www.熟女人妻精品国产| 香蕉久久夜色| 色哟哟哟哟哟哟| 国产蜜桃级精品一区二区三区| 日本a在线网址| 嫁个100分男人电影在线观看| 91大片在线观看| 欧美日韩黄片免| 午夜福利免费观看在线| 伦理电影免费视频| 一进一出抽搐gif免费好疼| 一个人免费在线观看的高清视频| 99久久综合精品五月天人人| 精品一区二区三区视频在线观看免费| 香蕉丝袜av| 日本三级黄在线观看| 精品一区二区三区视频在线观看免费| 国产精品av久久久久免费| 亚洲国产精品久久男人天堂| av中文乱码字幕在线| 亚洲黑人精品在线| 国内精品久久久久精免费| 涩涩av久久男人的天堂| 黄片大片在线免费观看| 香蕉久久夜色| 亚洲精品美女久久久久99蜜臀| x7x7x7水蜜桃| 久久香蕉精品热| 亚洲中文av在线| 桃色一区二区三区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 岛国视频午夜一区免费看| 精品无人区乱码1区二区| 久久久久国产精品人妻aⅴ院| www.自偷自拍.com| 18美女黄网站色大片免费观看| 国产在线精品亚洲第一网站| 欧美在线黄色| 制服诱惑二区| 日韩精品中文字幕看吧| www.999成人在线观看| 精品欧美国产一区二区三| 50天的宝宝边吃奶边哭怎么回事| 亚洲五月色婷婷综合| 国产主播在线观看一区二区| 满18在线观看网站| 精品午夜福利视频在线观看一区| 丁香六月欧美| 97碰自拍视频| 嫩草影院精品99| 日本vs欧美在线观看视频| 不卡av一区二区三区| 99国产极品粉嫩在线观看| 天堂动漫精品| 亚洲国产欧美一区二区综合| 成人亚洲精品av一区二区| 精品卡一卡二卡四卡免费| 亚洲精品国产区一区二| 国产高清激情床上av| 亚洲成人国产一区在线观看| 69av精品久久久久久| 久久久国产成人免费| www国产在线视频色| 欧美日本视频| 黑人巨大精品欧美一区二区蜜桃| АⅤ资源中文在线天堂| 两性夫妻黄色片| 成人18禁高潮啪啪吃奶动态图| 麻豆久久精品国产亚洲av| 久久人妻熟女aⅴ| 50天的宝宝边吃奶边哭怎么回事| 国产主播在线观看一区二区| 十分钟在线观看高清视频www| 麻豆成人av在线观看| 国产一区二区三区在线臀色熟女| 国产在线精品亚洲第一网站| 在线观看一区二区三区| 免费不卡黄色视频| 悠悠久久av| 亚洲国产精品成人综合色| 久久性视频一级片| 欧美另类亚洲清纯唯美| 一级毛片女人18水好多| av电影中文网址| 99久久久亚洲精品蜜臀av| 亚洲国产精品sss在线观看| videosex国产| 国产麻豆成人av免费视频| 日日摸夜夜添夜夜添小说| 日本黄色视频三级网站网址| 国产真人三级小视频在线观看| 男女之事视频高清在线观看| 首页视频小说图片口味搜索| 精品久久蜜臀av无| 三级毛片av免费| 啦啦啦观看免费观看视频高清 | 久热这里只有精品99| 91av网站免费观看| 亚洲人成伊人成综合网2020| 91九色精品人成在线观看| 欧美日本亚洲视频在线播放| 极品教师在线免费播放| 久久久久久亚洲精品国产蜜桃av| 757午夜福利合集在线观看| 日本 av在线| 日韩一卡2卡3卡4卡2021年| 一区二区三区激情视频| 一级片免费观看大全| 91在线观看av| 亚洲av日韩精品久久久久久密| 久久伊人香网站| 婷婷精品国产亚洲av在线| 亚洲五月婷婷丁香| av免费在线观看网站| 久久中文字幕人妻熟女| 欧美最黄视频在线播放免费| 宅男免费午夜| 亚洲第一电影网av| 美女午夜性视频免费| 非洲黑人性xxxx精品又粗又长| 亚洲五月色婷婷综合| 亚洲狠狠婷婷综合久久图片| 欧美成人免费av一区二区三区| 最近最新中文字幕大全电影3 | 国产三级在线视频| 中文字幕人妻丝袜一区二区| 久9热在线精品视频| 日日摸夜夜添夜夜添小说| 亚洲专区字幕在线| 久久影院123| 无限看片的www在线观看| 在线视频色国产色| 国产成人精品无人区| 性少妇av在线| 欧美国产精品va在线观看不卡| 久久婷婷人人爽人人干人人爱 | 午夜日韩欧美国产| 国产又爽黄色视频| 国产精品久久视频播放| 熟女少妇亚洲综合色aaa.| 一本综合久久免费| 久久国产精品男人的天堂亚洲| 18美女黄网站色大片免费观看| 国内精品久久久久精免费| 成人国产一区最新在线观看| 丝袜美腿诱惑在线| 欧洲精品卡2卡3卡4卡5卡区| 国产免费男女视频| 自拍欧美九色日韩亚洲蝌蚪91| 免费观看精品视频网站| 欧美av亚洲av综合av国产av| 日日干狠狠操夜夜爽| 欧美国产日韩亚洲一区| 欧美日韩黄片免| 国产亚洲精品久久久久5区| 久久影院123| 欧美在线黄色| 国产麻豆69| 变态另类成人亚洲欧美熟女 | 亚洲成国产人片在线观看| 韩国av一区二区三区四区| 国产亚洲精品第一综合不卡| 久久天躁狠狠躁夜夜2o2o| 国产精品野战在线观看| 欧美av亚洲av综合av国产av| 色精品久久人妻99蜜桃| 亚洲avbb在线观看| 久久九九热精品免费| 精品久久蜜臀av无| 国产精品国产高清国产av| 婷婷精品国产亚洲av在线| 久久久国产精品麻豆| 成人国语在线视频| 国产视频一区二区在线看| 欧美日本亚洲视频在线播放| 制服丝袜大香蕉在线| 校园春色视频在线观看| 日本黄色视频三级网站网址| 露出奶头的视频| www.精华液| 不卡一级毛片| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区在线不卡| 久久天堂一区二区三区四区| 制服人妻中文乱码| 在线免费观看的www视频| √禁漫天堂资源中文www| 亚洲国产毛片av蜜桃av| 国产亚洲欧美98| 成人国产一区最新在线观看| 非洲黑人性xxxx精品又粗又长| 纯流量卡能插随身wifi吗| 后天国语完整版免费观看| 在线观看www视频免费| 日韩大码丰满熟妇| 精品国产亚洲在线| 欧美色视频一区免费| 一区二区三区激情视频| 免费看十八禁软件| 亚洲五月色婷婷综合| 一级,二级,三级黄色视频| 亚洲色图av天堂| 狠狠狠狠99中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 91九色精品人成在线观看| 变态另类丝袜制服| 国产成人啪精品午夜网站| 午夜福利高清视频| 亚洲av日韩精品久久久久久密| 日本vs欧美在线观看视频| www日本在线高清视频| 亚洲视频免费观看视频| 国产熟女午夜一区二区三区| 国产成+人综合+亚洲专区| 亚洲av日韩精品久久久久久密| 免费一级毛片在线播放高清视频 | 精品熟女少妇八av免费久了| 国产精品一区二区精品视频观看| 久久精品91无色码中文字幕| 亚洲成人免费电影在线观看| 真人一进一出gif抽搐免费| 侵犯人妻中文字幕一二三四区| 色播亚洲综合网| 亚洲色图av天堂| 国产精品 国内视频| 一区二区三区国产精品乱码| 在线永久观看黄色视频| 制服人妻中文乱码| 美女高潮喷水抽搐中文字幕| 淫妇啪啪啪对白视频| 亚洲成人国产一区在线观看| 国产精品永久免费网站| 日本五十路高清| 国产色视频综合| 亚洲 欧美 日韩 在线 免费| 制服丝袜大香蕉在线| www.熟女人妻精品国产| 黄片小视频在线播放| 两个人看的免费小视频| 变态另类丝袜制服| 精品卡一卡二卡四卡免费| 免费高清在线观看日韩| 搞女人的毛片| 久久亚洲真实| 在线观看66精品国产| 久久精品成人免费网站| 亚洲精品国产区一区二| 国产熟女午夜一区二区三区| 纯流量卡能插随身wifi吗| 黄色片一级片一级黄色片| www.自偷自拍.com| 久久精品国产99精品国产亚洲性色 | 国产午夜精品久久久久久| 中文字幕久久专区| or卡值多少钱| 欧美 亚洲 国产 日韩一| 两个人视频免费观看高清| 亚洲午夜精品一区,二区,三区| 母亲3免费完整高清在线观看| 国产精品亚洲美女久久久| 午夜福利欧美成人| 人人妻人人澡欧美一区二区 | 亚洲第一电影网av| 国产av在哪里看| 制服诱惑二区| 又黄又爽又免费观看的视频| 日韩三级视频一区二区三区| 国内久久婷婷六月综合欲色啪| 日本五十路高清| 久久香蕉激情| 变态另类成人亚洲欧美熟女 | 91老司机精品| 国产在线观看jvid| 久久午夜亚洲精品久久| 色婷婷久久久亚洲欧美| 成人亚洲精品一区在线观看| xxx96com| 他把我摸到了高潮在线观看| 精品午夜福利视频在线观看一区| 精品一品国产午夜福利视频| 久久人人爽av亚洲精品天堂| 成人亚洲精品一区在线观看| 不卡av一区二区三区| 免费搜索国产男女视频| 中文亚洲av片在线观看爽| 男人操女人黄网站| 神马国产精品三级电影在线观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 搡老妇女老女人老熟妇| 美女国产高潮福利片在线看| 嫁个100分男人电影在线观看| 免费在线观看影片大全网站| av片东京热男人的天堂| 精品一品国产午夜福利视频| 久久久久久久久免费视频了| 淫秽高清视频在线观看| 无限看片的www在线观看| 啪啪无遮挡十八禁网站| 国产午夜福利久久久久久| 午夜精品久久久久久毛片777| 大型黄色视频在线免费观看| 可以免费在线观看a视频的电影网站| 欧美黑人精品巨大| 国内精品久久久久精免费| 久久久久久大精品| 岛国在线观看网站| 91国产中文字幕| 午夜福利高清视频| 成人国语在线视频| 色综合亚洲欧美另类图片| 亚洲欧美日韩高清在线视频| 国产成人精品在线电影| 国产精品国产高清国产av| 久久久久九九精品影院| 韩国精品一区二区三区| 校园春色视频在线观看| 亚洲成av片中文字幕在线观看| 久久久久久久久中文| 日日干狠狠操夜夜爽| 国产精品亚洲一级av第二区| 人人妻人人澡人人看| 丁香欧美五月| 欧美久久黑人一区二区| 久久久精品欧美日韩精品| 色老头精品视频在线观看| 嫁个100分男人电影在线观看| 999久久久国产精品视频| 国产成人欧美在线观看| 18禁国产床啪视频网站| 黄色毛片三级朝国网站| 最近最新中文字幕大全免费视频| 亚洲成人免费电影在线观看| 久久这里只有精品19| 一a级毛片在线观看| 91大片在线观看| 亚洲精品粉嫩美女一区| 少妇熟女aⅴ在线视频| 黄片大片在线免费观看| 人人妻,人人澡人人爽秒播| 美女扒开内裤让男人捅视频| 亚洲熟妇中文字幕五十中出| 韩国av一区二区三区四区| 一级片免费观看大全| 亚洲成人久久性| 国产单亲对白刺激| 免费在线观看完整版高清| 在线国产一区二区在线| 国产精品久久久久久人妻精品电影| 国产亚洲精品久久久久久毛片| 国产国语露脸激情在线看| 日韩 欧美 亚洲 中文字幕| 精品久久久久久久久久免费视频| 国产精品av久久久久免费| 一二三四在线观看免费中文在| 在线观看日韩欧美| 美女大奶头视频| 久久热在线av| 欧美日韩亚洲综合一区二区三区_| 久久国产精品影院| 操出白浆在线播放| 午夜a级毛片| 亚洲 欧美一区二区三区| 午夜亚洲福利在线播放| 动漫黄色视频在线观看| 99久久精品国产亚洲精品| 婷婷六月久久综合丁香| 精品久久久久久久人妻蜜臀av | 亚洲av片天天在线观看| a级毛片在线看网站| 一进一出抽搐gif免费好疼| 亚洲欧美精品综合久久99| 亚洲男人天堂网一区| 好看av亚洲va欧美ⅴa在| www.www免费av| 国产欧美日韩一区二区三区在线| 欧美另类亚洲清纯唯美| 国产午夜精品久久久久久| 色婷婷久久久亚洲欧美| 久久久国产成人免费| 美女午夜性视频免费| 免费观看精品视频网站| 最近最新中文字幕大全免费视频| 老汉色∧v一级毛片| 亚洲第一欧美日韩一区二区三区| 欧美黑人欧美精品刺激| 99国产精品一区二区三区| 久久九九热精品免费| 天天躁狠狠躁夜夜躁狠狠躁| 淫妇啪啪啪对白视频| 搡老妇女老女人老熟妇| 黄色 视频免费看| 大码成人一级视频| 亚洲av成人av| 一边摸一边抽搐一进一出视频| 午夜精品久久久久久毛片777| 亚洲美女黄片视频| 少妇的丰满在线观看| 99re在线观看精品视频| svipshipincom国产片| 亚洲欧美精品综合久久99| 丁香欧美五月| 久久伊人香网站| 变态另类丝袜制服| 色播在线永久视频| 亚洲av电影在线进入| 色播在线永久视频| 久久中文字幕人妻熟女| 欧美性长视频在线观看| 热re99久久国产66热| 国产精品爽爽va在线观看网站 | 国产乱人伦免费视频| 精品国产乱码久久久久久男人| 中文字幕高清在线视频| 淫秽高清视频在线观看| 99精品在免费线老司机午夜| 亚洲精品美女久久av网站| 真人做人爱边吃奶动态| 中文字幕人妻丝袜一区二区| 可以在线观看毛片的网站| 亚洲人成77777在线视频| 91精品国产国语对白视频| 亚洲av日韩精品久久久久久密| 久久 成人 亚洲| 91成年电影在线观看| 国产亚洲精品av在线| 国产亚洲精品久久久久5区| 无限看片的www在线观看| 国产av在哪里看| 男女之事视频高清在线观看| 99国产精品一区二区三区| 一夜夜www| av电影中文网址| 大型av网站在线播放| 一二三四在线观看免费中文在| 亚洲成av人片免费观看| ponron亚洲| 97超级碰碰碰精品色视频在线观看| 97人妻精品一区二区三区麻豆 | 亚洲av第一区精品v没综合| 久久香蕉精品热| 日本一区二区免费在线视频| 无遮挡黄片免费观看| 黄色a级毛片大全视频| 一个人观看的视频www高清免费观看 | 亚洲成人免费电影在线观看| 18禁黄网站禁片午夜丰满| 国产成人免费无遮挡视频| 亚洲av美国av| 国产97色在线日韩免费| 国产精品久久久久久亚洲av鲁大| 国产精品久久电影中文字幕| 91字幕亚洲| 亚洲国产中文字幕在线视频| 日本精品一区二区三区蜜桃| 丝袜美足系列|