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

    冰封對烏梁素海水質(zhì)污染特征的影響

    2023-04-19 02:52:32趙同國李衛(wèi)平田業(yè)蓬原凌飛
    灌溉排水學報 2023年3期
    關(guān)鍵詞:烏梁素海冰封類水

    張 巖,王 楠,趙同國,李衛(wèi)平,田業(yè)蓬,原凌飛

    冰封對烏梁素海水質(zhì)污染特征的影響

    張 巖1,王 楠1,趙同國1,李衛(wèi)平2,田業(yè)蓬3,原凌飛4

    (1.煙臺大學 土木工程學院,山東 煙臺 264005;2.內(nèi)蒙古科技大學 能源與環(huán)境學院,內(nèi)蒙古 包頭 014010;3.煙臺大學 特里爾可持續(xù)發(fā)展技術(shù)學院,山東 煙臺 264005;4.煙臺市套子灣污水處理有限公司,山東 煙臺 264013)

    【目的】明確烏梁素海冰封期的水質(zhì)污染特征,揭示冰封對湖泊水質(zhì)等級的影響。【方法】于2021年1月在烏梁素海布設(shè)10個采樣點,采集冰層及冰下水樣品。參照《地表水環(huán)境質(zhì)量標準》(GB 3838—2002)和烏梁素海的污染特征,選取CODMn、NH3-N、TP、TN、氟化物和銅6項水質(zhì)指標作為目標污染物,采用單因子水質(zhì)標識指數(shù)法、綜合水質(zhì)標識指數(shù)法和主成分分析法對烏梁素海冰封期與非冰封期的水質(zhì)進行綜合評價?!窘Y(jié)果】各項目標污染物的單因子水質(zhì)標識指數(shù)為:冰封期冰下水≥非冰封期≥冰封期冰融水;6項目標污染物的污染程度依次為:TN>TP>CODMn>NH3-N>銅>氟化物,其中CODMn、NH3-N、TN質(zhì)量濃度都發(fā)生了水質(zhì)等級的變化;各采樣點在冰封期的綜合水質(zhì)標識指數(shù)結(jié)果全部劣于非冰封期,其中采樣點I12由Ⅱ類惡化為Ⅲ類,其余采樣點水質(zhì)等級雖未發(fā)生變化,但其質(zhì)量濃度距離相應的水質(zhì)等級下限值更加接近;烏梁素海的整體水質(zhì)受NH3-N、TN和TP水質(zhì)因子的影響較大,較高的氮、磷值將為烏梁素海水體富營養(yǎng)化帶來潛在風險;各采樣點的水質(zhì)優(yōu)劣順序在空間上呈出水區(qū)>湖中區(qū)>進水區(qū)的變化?!窘Y(jié)論】冬季結(jié)冰過程中,冰對目標污染物具有排斥效應,污染物由冰體向冰下水遷移致使冰下水中污染物質(zhì)量濃度增大,加劇了烏梁素海的水質(zhì)惡化,這為冰封期湖泊水環(huán)境管理帶來嚴峻考驗。

    烏梁素海;結(jié)冰效應;冰封期;主成分分析法;水質(zhì)標識指數(shù)

    0 引 言

    【研究意義】湖泊是地球上淡水資源的重要組成部分,是維護水域生態(tài)系統(tǒng)健康發(fā)展的重要保障。在高緯度、高海拔地區(qū),有5 000萬個以上湖泊存在冰封期[1],形成的冰層阻擋了冰下水體與外界的物質(zhì)和能量交換[2],降低了各種生化反應速率[3],使冰下水體自凈能力大大減弱[4-5]。由于結(jié)冰過程中冰層的排斥作用,污染物向冰下水遷移,在冰下水中甚至是沉積物中聚集,使冰下水環(huán)境受到嚴重破壞,湖泊冰封期水污染嚴重[6-7]。在春季冰層融化過程中,冰體內(nèi)的少量污染物會在早期集中釋放到冰下水中,也會對冰下水環(huán)境再次產(chǎn)生影響[8]。

    【研究進展】近年來,針對冰封期湖泊國內(nèi)外學者開展了越來越多的研究,這些研究大多集中在污染物、營養(yǎng)鹽和浮游動植物的分布及污染特征以及結(jié)冰對湖泊冰下水污染影響。孫悅等[9]研究了白洋淀冬季冰封期的水污染特征,結(jié)果表明TN、TP是白洋淀冬季冰封期冰下水水質(zhì)的限制因子;王司陽等[10]以烏梁素海、奈倫湖為研究對象,對黃河流域上游干旱地區(qū)冰封期湖泊營養(yǎng)狀態(tài)特征進行了分析,發(fā)現(xiàn)冰蓋的形成使湖泊富營養(yǎng)化程度加劇;劉曉旭等[11]以內(nèi)蒙古達里諾爾湖為研究對象,對不同相態(tài)下湖水中的營養(yǎng)鹽濃度分布特征和定量關(guān)系進行了研究,結(jié)果表明冬季結(jié)冰過程中達里湖水體中的氮和磷存在濃縮效應,導致冰封期湖內(nèi)營養(yǎng)鹽濃度超過國家Ⅴ類水質(zhì)標準;李佳等[12]分析了烏梁素海冰封期浮游藻類群落結(jié)構(gòu)以及與營養(yǎng)物質(zhì)之間的關(guān)系,揭示了營養(yǎng)物質(zhì)在冰-水介質(zhì)的分布規(guī)律,并對冰封期的水質(zhì)進行了評價;盧興順等[13]分析了烏梁素海冰封期污染物的分布遷移,揭示了烏梁素海冰封期營養(yǎng)物質(zhì)分布規(guī)律,為湖泊富營養(yǎng)化治理等研究提供理論依據(jù)。孫馳[14]基于第一性原理分析了冰的物理特點對營養(yǎng)元素分布的影響,揭示了重金屬在冰-水中的分布特征?!厩腥朦c】目前對非冰封和冰封條件下的湖泊水質(zhì)污染特征對比以及結(jié)冰效應對水質(zhì)等級變化的影響卻鮮有報道?!緮M解決的關(guān)鍵問題】本研究將烏梁素海進行網(wǎng)格化處理,設(shè)置10個代表性采樣點,參照《地表水環(huán)境質(zhì)量標準》(GB 3838—2002),結(jié)合了烏梁素海的環(huán)境條件,選取較為典型的CODMn、NH3-N、TP、TN、氟化物和銅6項水質(zhì)指標,運用單因子水質(zhì)標識指數(shù)、綜合水質(zhì)標識指數(shù)以及主成分分析對該水體非冰封期和冰封期的水質(zhì)進行分析,旨在從不同角度明確冰封對烏梁素海水質(zhì)等級和水污染特征的影響。

    1 材料與方法

    1.1 研究區(qū)概況

    烏梁素海(40°36'—41°03'N,108°43'—108°57'E)位于內(nèi)蒙古自治區(qū)巴彥淖爾市烏拉特前旗境內(nèi),河套灌區(qū)最東部,是黃河流域內(nèi)典型的寒旱區(qū)湖泊濕地,水域面積為293 km2,湖區(qū)水深在0.5~1.5 m之間,最大水深4 m,流域基本概況如圖1所示[15]。烏梁素海是河套灌區(qū)工農(nóng)業(yè)退水的主要承泄渠道[16],隨著工農(nóng)業(yè)迅速發(fā)展,大量污染物如氮、磷經(jīng)由各排干進入烏梁素海,對烏梁素海水質(zhì)的影響嚴重,某些污染物的質(zhì)量濃度遠超出了《地表水環(huán)境質(zhì)量標準》中的限定值,甚至惡化到Ⅴ類水標準,對湖泊的生態(tài)環(huán)境造成嚴重破壞[17]。湖泊于每年11月中下旬開始結(jié)冰,冰封期長達5個月,冰層厚度為0.3~0.6 m,約占水體的1/3[18]。冰蓋形成后整個冰下水生態(tài)系統(tǒng)處于“靜止”狀態(tài)[2],大大削弱了湖泊水體的自凈能力[4]。

    圖1 烏梁素海流域概況和采樣點布設(shè)[15,18]

    1.2 樣品采集和數(shù)據(jù)處理

    根據(jù)河套灌區(qū)各排干與烏梁素海的相對位置,依據(jù) 2 km×2 km 的梅花型方式,將研究區(qū)域網(wǎng)格化處理,布設(shè)了10個采樣點(圖1),包括進水區(qū):J11、I12、西大灘L11、北部明水區(qū)域L15,湖中區(qū):大卜洞O10、蘆葦區(qū)域N13、旅游區(qū)Q10、Q8,出水區(qū):海壕V3、二點T5,于2021年1月進行冰層和冰下水的采集:用冰鉆采集冰芯樣品,將其放置在切割板上,用冰鋸切分為3段冰柱,并將其分別放在塑料瓶中,在室溫下自然融化;用“注射器”式采水器通過冰孔抽取各采樣點的冰下水樣。參照《水和廢水監(jiān)測分析方法》[20],對采集的冰樣以及冰下水樣進行檢測,用以表征烏梁素海的水質(zhì)特征。

    為避免研究過程中外源污染物以及降水等原因?qū)Ω黜椝|(zhì)指標的影響,根據(jù)物質(zhì)守恒原理,非冰封期各指標質(zhì)量濃度1計算式為(不考慮冰和水的密度差異):

    式中:ice為冰層中各項指標的質(zhì)量濃度(mg/L);ice為冰層的體積(m3);water為冰下水體中各項指標的質(zhì)量濃度(mg/L);water為冰下水的體積(m3);1為計算所得的非冰封期各項水質(zhì)指標的質(zhì)量濃度(mg/L)。

    分配系數(shù)用以表征污染物在結(jié)冰過程中由冰層向冰下水層遷移的能力,計算式為:

    1.3 水質(zhì)評價方法

    選用單因子水質(zhì)標識指數(shù)、綜合水質(zhì)標識指數(shù)和主成分分析對烏梁素海冰封期和非冰封期水質(zhì)進行評價。通過單因子水質(zhì)標識指數(shù)判斷烏梁素海各采樣點的水質(zhì)指標是否達到《地表水環(huán)境質(zhì)量標準》(GB 3838—2002)中的水域功能分類,將檢測數(shù)據(jù)在同一類水質(zhì)類別里進行定量比較,得出縱向空間水質(zhì)變化趨勢;運用主成分分析能夠?qū)⒏呔S度數(shù)據(jù)化為低維度數(shù)據(jù),以減少冗余信息所造成的誤差,同時篩選出能囊括絕大部分信息的水質(zhì)指標,對烏梁素海的綜合水質(zhì)進行分析,并對各個采樣點的水質(zhì)進行排序,定性、定量的描述烏梁素海不同時期的綜合水質(zhì);通過綜合水質(zhì)標識指數(shù)對單因子水質(zhì)標識指數(shù)的結(jié)果進行進一步分析,結(jié)合《地表水環(huán)境質(zhì)量標準》(GB 3838—2002)評價湖泊的綜合水質(zhì)類別,將冰封期與非冰封期的水質(zhì)等級進行對比,分析冰封對湖泊水質(zhì)等級的影響。

    1.3.1 單因子水質(zhì)標識指數(shù)

    單因子水質(zhì)指數(shù)P[21]由1位整數(shù)、小數(shù)點后2位或3位有效數(shù)字組成,計算式為:

    式中:1代表水質(zhì)指標的水質(zhì)類別;2代表水質(zhì)數(shù)據(jù)在1類水質(zhì)變化區(qū)間中所處的位置,根據(jù)公式按四舍五入的原則計算確定;3代表水質(zhì)類別與功能區(qū)劃設(shè)定類別的比較結(jié)果,視評價指標的污染程度,3為1位或2位有效數(shù)字。

    1.3.2 主成分分析法

    通過主成分分析法將烏梁素海復雜的試驗數(shù)據(jù)進行降維處理,分析各項污染物所攜帶的信息,篩選出比原始污染物數(shù)目少,但能影響總體水質(zhì)的代表性污染物[9],再進行成分得分系數(shù)矩陣等相關(guān)方程和公式的計算,得到10個采樣點的主成分分析的綜合得分,定量描述各采樣點的污染程度,并對各個采樣點的水質(zhì)優(yōu)劣進行排序,從空間上確定湖泊不同區(qū)域的污染程度。

    1.3.3 綜合水質(zhì)標識指數(shù)

    綜合水質(zhì)標識指數(shù)wq[22]由整數(shù)位和3位或4位小數(shù)組成,其計算式為:

    式中:1、2由計算獲得,3和4根據(jù)比較結(jié)果得到。其中,1為綜合水質(zhì)類別;2為綜合水質(zhì)在1類水質(zhì)變化區(qū)間內(nèi)所處位置;3為參與綜合水質(zhì)評價的水質(zhì)指標中,劣于水環(huán)境功能區(qū)目標的單項指標個數(shù);4為綜合水質(zhì)類別與水體功能區(qū)類別的比較結(jié)果,視綜合水質(zhì)的污染程度,4為1位或2位有效數(shù)字。

    2 結(jié)果與分析

    2.1 烏梁素海冰封期湖泊水深與冰厚分布

    為了更好地表征各污染物向冰下水的遷移程度,對10個采樣點的水深以及冰厚進行了現(xiàn)場測量,結(jié)果見圖2。烏梁素海全湖的平均水深為2.19 m,西北和南部的水深較深,位于南部的采樣點T5最深,為3.13 m;冰厚分布較均勻,在0.35~0.6 m之間,平均值為0.49 m,湖區(qū)中部地區(qū)的冰厚值多大于平均值。

    圖2 烏梁素海各采樣點冰厚和水深變化

    2.2 烏梁素海冰封期各項水質(zhì)因子特征

    由圖3可知,6項目標污染物在非冰封期湖水、冰融水和冰下水的含量呈“凹”字形分布,具體表現(xiàn)為:冰下水>非冰封期湖水>冰融水,其中CODMn和TP在各采樣點質(zhì)量濃度分布較均勻,但遷移系數(shù)在采樣點L11較大,分別為0.52和1.04,表明該采樣點的CODMn和TP向冰下水遷移能力較弱,更容易被冰晶捕獲在冰體中,導致冰融水中CODMn、TP質(zhì)量濃度較大。J11采樣點的NH3-N和TN質(zhì)量濃度遠大于其他采樣點,遷移系數(shù)卻遠小于其他采樣點,分別為0.03和0.02,表明NH3-N和TN的遷移能力較強,加重了冰下水的氮源污染,春季冰融期易引起水體富營養(yǎng)化。采樣點L11的TN遷移系數(shù)較高,為0.23,污染物向冰下水遷移的量較少,導致冰融水中的TN質(zhì)量濃度明顯大于其他采樣點。氟化物的質(zhì)量濃度在整個湖區(qū)分布較均勻,但其遷移系數(shù)在J11采樣點明顯大于其他采樣點,為0.30,遷移能力較弱,所以J11采樣點的氟化物冰中質(zhì)量濃度稍高。采樣點I12的銅質(zhì)量濃度明顯低于其他點,但遷移系數(shù)卻是最大,為0.06,其他均在0.02以下,說明I12采樣點低質(zhì)量濃度的銅遷移能力較弱,可能是由于采樣點I12的銅質(zhì)量濃度較低,進入冰封期后向冰下水的遷移量相對較少,導致遷移系數(shù)較大。

    2.3 烏梁素海冰封期單因子水質(zhì)標識指數(shù)評價

    依據(jù)《地表水質(zhì)量標準》(GB 3838—2002)中地表水水域的水環(huán)境功能與保護目標的功能區(qū)對象劃分,烏梁素海的水域功能分類為Ⅲ類,本研究以地表Ⅲ類水為限值,對烏梁素海非冰封期和冰封期冰下水和冰中的6項指標進行單因子水質(zhì)標識指數(shù)計算,結(jié)果如圖4所示。

    圖3 烏梁素海各采樣點目標污染物質(zhì)量濃度

    圖4 烏梁素海各采樣點目標污染物的單因子水質(zhì)標識指數(shù)

    6項水質(zhì)指標的單因子水質(zhì)標識指數(shù)結(jié)果表現(xiàn)為:冰下水≥非冰封期湖水≥冰融水。各采樣點冰融水中污染物的質(zhì)量濃度未超過Ⅲ類水標準,但冰封期冰下水和非冰封期結(jié)果明顯大于冰融水,其中TN和TP污染最為嚴重,所有采樣點均超出了地表Ⅲ類水標準,其次為CODMn和NH3-N,污染較輕的是銅和氟化物。湖泊從非冰封期進入冰封期后,有6個采樣點的CODMn下降了1個水質(zhì)等級,其中采樣點J11、L11、Q8、Q10和V3由Ⅱ類水降低到Ⅲ類水,I12采樣點由Ⅲ類水下降至Ⅳ類水;采樣點J11和采樣點Q8的NH3-N下降了1個水質(zhì)等級,其中采樣點J11從Ⅲ類水惡化到Ⅳ類水,其余采樣點NH3-N質(zhì)量濃度變化未達到等級變化程度;采樣點I12、L15和Q8的TN質(zhì)量濃度惡化了1個水質(zhì)等級,其中采樣點I12惡化為Ⅴ類水,其余2個采樣點為Ⅳ類水,其余采樣點雖未發(fā)生水質(zhì)等級的變化,但TN質(zhì)量濃度距離相應水質(zhì)等級區(qū)間的下限值接近了21%~70%,甚至采樣點J11和采樣點L11的TN單因子水質(zhì)標識指數(shù)為7.54和6.43,已經(jīng)為劣Ⅴ類水。所有采樣點的TP雖未發(fā)生水質(zhì)等級的下降,但都超出了Ⅲ類水標準。氟化物非冰封期和冰封期冰下水的單因子水質(zhì)標識指數(shù)結(jié)果在1.20~1.50之間,整體變化幅度不大。銅的非冰封期和冰封期冰下水的結(jié)果完全重合,這是因為烏梁素海水體中的銅質(zhì)量濃度較低,雖然冰封期會受到冰蓋的排斥作用而向冰下水遷移,但其質(zhì)量濃度遠遠達不到水質(zhì)等級單因子水質(zhì)標識指數(shù)的變化幅度。

    2.4 烏梁素海冰封期綜合水質(zhì)標識指數(shù)評價

    為對冰封期湖泊的橫向空間的污染特征進行比較,采用綜合水質(zhì)標識指數(shù)能夠進一步明確冰封對各采樣點水質(zhì)變化的影響。由圖5可知,湖泊冰封期綜合水質(zhì)標識指數(shù)明顯大于非冰封期的。采樣點I12綜合水質(zhì)下降了1個等級,由Ⅱ類水下降為Ⅲ類水。其余采樣點水質(zhì)等級雖未發(fā)生變化,但這些采樣點各項指標的質(zhì)量濃度距離相應的水質(zhì)等級下限值接近10%~31%:差值最大的是采樣點L15、Q8,表明進入冰封期后,污染最為嚴重;差值最小的是采樣點L11、T5,表明污染較輕。冰封期各個采樣點參與評價的水質(zhì)指標未達到水環(huán)境功能區(qū)目標的數(shù)量相較于非冰封期增加了5個。

    圖5 烏梁素海各采樣點綜合水質(zhì)標識指數(shù)

    2.5 烏梁素海冰封期污染主成分分析

    為了進一步確定冰封期湖泊水質(zhì)的重要影響因子,采用主成分分析對冰封期湖泊空間上水質(zhì)優(yōu)劣進行排序。由表1的因子荷載矩陣可知,3個主成分的累計貢獻率達到86.703%,表明這3個主成分可反映數(shù)據(jù)的大部分信息。貢獻率最大的是第1個主成分為45.296%,最小的是第3個主成分為13.931%。與第1主成分關(guān)系密切的相關(guān)指標有NH3-N和TN;與第2主成分關(guān)系密切的指標有TP。由此可知,影響烏梁素海冰封期綜合水質(zhì)的主要因子為NH3-N、TN和TP。由表2主成分分析的因子荷載矩陣可知,3個主成分的累計貢獻率達到84.376%,表明這3個主成分可反映數(shù)據(jù)的大部分信息。3個主成分中貢獻率最大為44.244%,最小是24.028%。與第1主成分關(guān)系密切的相關(guān)指標有NH3-N和TN,因子負載為0.921和0.917;與第2主成分關(guān)系密切的指標有氟化物;與第3主成分關(guān)系密切的指標有CODMn。因此,影響烏梁素海非冰封期綜合水質(zhì)的主要因子為NH3-N、TN、氟化物。

    表1 冰封期主成分分析因子荷載矩陣

    表2 非冰封期主成分分析因子荷載矩陣

    通過因子荷載矩陣的分析可以反映出烏梁素海水質(zhì)的主要影響因子是NH3-N、TN和TP。根據(jù)成分得分系數(shù)矩陣等相關(guān)方程和計算式得出10個采樣點的主成分分析的綜合得分(圖6),綜合得分越高,水質(zhì)越差。各采樣點水質(zhì)優(yōu)劣順序為:J11

    圖6 烏梁素海主成分分析綜合得分

    3 討 論

    3.1 結(jié)冰過程中污染物的遷移機理

    由烏梁素海冰封期和非冰封期的單因子水質(zhì)標識指數(shù)、綜合水質(zhì)標識指數(shù)和主成分分析評價的評價結(jié)果可知,結(jié)冰過程中CODMn、NH3-N、TP、TN、氟化物和銅6項水質(zhì)指標均呈現(xiàn)出由冰體向冰下水體遷移的規(guī)律,這可從湖冰的形成過程進行解釋:當冬季湖泊表面溫度下降到0 ℃時,表層水達到過冷狀態(tài)后水分子開始有冰核產(chǎn)生并形成細小、平整且不含雜質(zhì)的冰晶,這些冰晶相互連接便可形成柔性冰,柔性冰經(jīng)凍結(jié)變厚,形成堅硬冰蓋,成核形成冰晶,并逐漸形成冰水界層[23],此時空氣中的冷量能只能通過湖泊表面向下傳遞,由于氮、磷等污染物的凝固點較低,冷能僅能夠使在冰-水界面處的水分子達到凝固點而平穩(wěn)析出,烏梁素海冰封期6種污染物由冰體向冰下水遷移,水分子形成冰晶的過程中會將液態(tài)的污染物質(zhì)排斥到冰水界層附近;同時湖泊結(jié)冰過程中也會產(chǎn)生熱通量,進而產(chǎn)生溶質(zhì)通量,溶質(zhì)通量會比冰更快的脫離結(jié)冰區(qū)域,即污染物質(zhì)會比冰晶更快的脫離結(jié)冰區(qū),進而進入到冰水界層。溶質(zhì)通量的驅(qū)動因素為溫度梯度和質(zhì)量濃度梯度:溫度梯度會驅(qū)動物質(zhì)從低溫區(qū)域遷移到高溫區(qū)域(索雷效應)[24],即隨著結(jié)冰過程的進行,各污染物會由冰體遷移到冰下水體,只有少部分污染物以“污染物胞”的形式存在于冰體中(圖7);質(zhì)量濃度梯度則會驅(qū)動物質(zhì)從高質(zhì)量濃度遷移至低質(zhì)量濃度區(qū)域(菲克擴散定律)[24],即污染物質(zhì)會從質(zhì)量濃度較高的冰水界層遷移至冰下水中,導致原本質(zhì)量濃度較為均勻的水體發(fā)生變化,冰下水中污染物質(zhì)量濃度增加。

    圖7 冰結(jié)構(gòu)剖面示意圖

    3.2 湖泊冰封的環(huán)境效應和啟示

    結(jié)冰過程中污染物的遷移效應使得烏梁素海冰封期冰下水體中污染物的質(zhì)量濃度增加,冰下水環(huán)境質(zhì)量的惡化會隨著結(jié)冰厚度的逐漸增加而進一步加劇,特別是TN、NH3-N、TP的單因子水質(zhì)標識指數(shù)較高,表明冰封對烏梁素海富營養(yǎng)元素的環(huán)境效應更為明顯,這與楊文煥等[25]的研究結(jié)果一致。而冰下水中較高的氮、磷量將為翌年春季水體的富營養(yǎng)化帶來潛在風險,特別是在進水區(qū)域,隨著河套灌區(qū)春灌排水的流入,富營養(yǎng)化的風險系數(shù)還將提高,這需要引起湖泊水管理部門的重視。同時,由于烏梁素海水深較淺,冰體所占比例較高,冰下水中污染物的濃縮系數(shù)較大,因此可從淺水湖泊冰封期污染特征的特殊性入手,探索淺水湖泊冰封期水環(huán)境治理的措施。

    冰封對湖泊水質(zhì)污染特征的影響也提示我們需要尤為關(guān)注水體冰封期冰下水環(huán)境質(zhì)量的變化,盡管冰體的覆蓋使得冰下水環(huán)境的變化更加“隱蔽”,但冰下水環(huán)境的變化絕非“靜止”。特別是對于飲用水水源地而言,隨著結(jié)冰過程的進行,冰下水質(zhì)會逐漸變差,加之低溫條件下水處理廠的化學處理工藝的效率會降低,因此如何能在水體冰封期確保出廠水達到《生活飲用水衛(wèi)生標準》(GB 5749—2022)的要求帶來新的課題和挑戰(zhàn)。

    4 結(jié) 論

    1)烏梁素海冰封期的綜合水質(zhì)標識指數(shù)結(jié)果全部劣于非冰封期,甚至出現(xiàn)了水質(zhì)等級的下降,表明結(jié)冰過程對污染物具有濃縮效應,致使污染物被排斥到冰下水體。

    2)烏梁素海各區(qū)域的水質(zhì)優(yōu)劣順序為:進水區(qū)<湖中區(qū)<出水區(qū),表明湖區(qū)進水區(qū)的水質(zhì)污染較嚴重,應加強對匯入湖區(qū)的農(nóng)業(yè)退水和農(nóng)業(yè)廢水等外源污染物排放的治理。

    3)烏梁素海的水質(zhì)受CODMn、NH3-N、TN、TP水質(zhì)因子的影響較大,較高的營養(yǎng)元素將為湖泊春季融冰期的富營養(yǎng)化帶來潛在風險。

    4)湖泊冰封期冰下水環(huán)境質(zhì)量的惡化,為冰封期湖泊的管理帶來嚴峻的考驗,也可為冰封期湖泊的治理提供新的思路,但需要進一步研究湖泊冰封期水環(huán)境治理的具體措施。

    [1] VERPOORTER C, KUTSER T, SEEKELL D A, et al. A global inventory of lakes based on high-resolution satellite imagery[J]. Geophysical Research Letters, 2014, 41(18): 6 396-6 402.

    [2] FANG X, STEFAN H G. Simulations of climate effects on water temperature, dissolved oxygen, and ice and snow covers in lakes of the contiguous US under past and future climate scenarios[J]. Limnology and Oceanography, 2009, 54(6): 2 359-2 370.

    [3] BELZILE C, GIBSON J A E, VINCENT W F. Colored dissolved organic matter and dissolved organic carbon exclusion from lake ice: Implications for irradiance transmission and carbon cycling[J]. Limnology and Oceanography, 2002, 47(5): 1 283-1 293.

    [4] 黃文峰, 韓紅衛(wèi), 牛富俊, 等. 季節(jié)性冰封熱融淺湖水溫原位觀測及其分層特征[J]. 水科學進展, 2016, 27(2): 280-289.

    HUANG Wenfeng, HAN Hongwei, NIU Fujun, et al. Field observations on water temperature and stratification in a seasonally ice-covered shallow thermokarst lake[J]. Advances in Water Science, 2016, 27(2): 280-289.

    [5] 李衛(wèi)平, 滕飛, 楊文煥, 等. 烏梁素海冰封期冰-水中污染物空間分布特征及污染評價[J]. 灌溉排水學報, 2020, 39(2): 122-128, 144.

    LI Weiping, TENG Fei, YANG Wenhuan, et al. Spatial distribution characteristics and pollution assessment of ice-water pollutants during the ice-sealing period in wuliangsuhai[J]. Journal of Irrigation and Drainage, 2020, 39(2): 122-128, 144.

    [6] 姜忠峰. 烏梁素海綜合需水分析及生態(tài)系統(tǒng)健康評價[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學, 2011.

    JIANG Zhongfeng. Analysis of ecological water requirement and ecosystem health assessment for lake wuliangsuhai[D]. Hohhot: Inner Mongolia Agricultural University, 2011.

    [7] 周石磊, 孫悅, 岳哿丞, 等. 雄安新區(qū)-白洋淀冬季冰封期水體好氧反硝化菌群落空間分布特征及驅(qū)動因素[J]. 環(huán)境科學, 2020, 41(5): 2 177-2 187.

    ZHOU Shilei, SUN Yue, YUE Gecheng, et al. Spatial distribution characteristics and driving factors of aerobic denitrification bacterial community structure from Baiyangdian Lake in Xiongan new area during the winter freezing period[J]. Environmental Science, 2020, 41(5): 2 177-2 187.

    [8] XUE S, CHEN J, TIE M. Release of dissolved organic matter from melting ice[J]. Environmental Progress & Sustainable Energy, 2016, 35(5): 1 458-1 467.

    [9] 孫悅, 李再興, 張藝冉, 等. 雄安新區(qū)-白洋淀冰封期水體污染特征及水質(zhì)評價[J]. 湖泊科學, 2020, 32(4): 952-963.

    SUN Yue, LI Zaixing, ZHANG Yiran, et al. Water pollution characteristics and water quality evaluation during the freezing period in Lake Baiyangdian of Xiongan New Area[J]. Journal of Lake Sciences, 2020, 32(4): 952-963.

    [10] 王司陽, 張笑欣, 田世民, 等. 黃河流域干旱區(qū)湖泊冰封期浮游植物群落結(jié)構(gòu)特征及影響因子研究[J]. 水利學報, 2020, 51(9): 1 070-1 079.

    WANG Siyang, ZHANG Xiaoxin, TIAN Shimin, et al. Study on phytoplankton community structure characteristics and its influencing factors of lakes in arid regions of the Yellow River Basin during ice-sealing period[J]. Journal of Hydraulic Engineering, 2020, 51(9): 1 070-1 079.

    [11] 劉曉旭, 李暢游, 李文寶, 等. 冰封期達里諾爾湖同位素與營養(yǎng)鹽分布特征及關(guān)系的定量分析[J]. 湖泊科學, 2015, 27(6): 1 159-1 167.

    LIU Xiaoxu, LI Changyou, LI Wenbao, et al. The distribution and relationship of isotope and nutrient during freeze-up period in the Lake Dalinuoer: A quantitative approach[J]. Journal of Lake Sciences, 2015, 27(6): 1 159-1 167.

    [12] 李佳, 侯俊青, 趙子聞, 等. 烏梁素海冰封期浮游藻類分布特征研究及水質(zhì)評價[J]. 環(huán)境科學與技術(shù), 2019, 42(9): 61-67.

    LI Jia, HOU Junqing, ZHAO Ziwen, et al. Distribution characteristics of phytoplankton and water quality evaluation in the ice-sealing period of lake ulansuhai[J]. Environmental Science & Technology, 2019, 42(9): 61-67.

    [13] 盧興順, 丁曉宇, 趙子聞, 等. 烏梁素海冰封期污染物分布及遷移特征研究[J]. 水生態(tài)學雜志, 2021, 42(6): 41-48.

    LU Xingshun, DING Xiaoyu, ZHAO Ziwen, et al. Distribution and exchange of pollutants during the icebound season of ulansuhai lake[J]. Journal of Hydroecology, 2021, 42(6): 41-48.

    [14] 孫馳. 基于第一性原理的烏梁素海冰、水介質(zhì)中重金屬遷移特征研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學, 2019.

    SUN Chi. Transport characteristics of heavy metals based on first principle in ice and water medium in the lake ulansuhai[D]. Hohhot: Inner Mongolia Agricultural University, 2019.

    [15] 翟佳倫, 史小紅, 劉禹, 等. 烏梁素海冰封期水溫與溶解氧濃度變化研究[J]. 干旱區(qū)研究, 2021, 38(3): 629-639.

    ZHAI Jialun, SHI Xiaohong, LIU Yu, et al. Change law of water temperature and dissolved oxygen concentration of Wuliangsu Sea in icebound period[J]. Arid Zone Research, 2021, 38(3): 629-639.

    [16] 盧興順. 植物碳源強化人工濕地處理烏梁素海農(nóng)田退水效能研究[D]. 呼和浩特: 內(nèi)蒙古大學, 2022.

    LU Xingshun. Study on the efficiency of plant carbon sourceenhanced constructed wetlands for treating agricultural drainage water in the ulansuhai lake[D]. Hohhot: Inner Mongolia University, 2022.

    [17] 朱鵬航. 基于PCLake模型的烏梁素海水質(zhì)突變及閾值識別研究[D]. 呼和浩特: 內(nèi)蒙古大學, 2022.

    ZHU Penghang. A study on water quality mutation and threshold identification in Wuliangsu lake based on PCLake model[D]. Hohhot: Inner Mongolia University, 2022.

    [18] 全棟, 史小紅, 趙勝男, 等. 2006—2017年烏梁素海夏季水體營養(yǎng)狀態(tài)及影響因子[J]. 湖泊科學, 2019, 31(5): 1 259-1 267.

    QUAN Dong, SHI Xiaohong, ZHAO Shengnan, et al. Eutrophication of Lake Ulansuhai in 2006—2017 and its main impact factors[J]. Journal of Lake Sciences, 2019, 31(5): 1 259-1 267.

    [19] 史小紅. 烏梁素海營養(yǎng)元素及其存在形態(tài)的數(shù)值模擬分析[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學, 2007.

    SHI Xiaohong. Modelling analysis of eutrophic elements and their speciation for Wuliangsuhai lake in Inner Mongolia[D]. Hohhot: Inner Mongolia Agricultural University, 2007.

    [20] 國家環(huán)境保護總局. 水和廢水監(jiān)測分析方法[M]. 北京: 中國環(huán)境科學出版社, 2002.

    [21] 徐祖信. 我國河流單因子水質(zhì)標識指數(shù)評價方法研究[J]. 同濟大學學報(自然科學版), 2005, 33(3): 321-325.

    XU Zuxin. Single factor water quality identification index for environmental quality assessment of surface water[J]. Journal of Tongji University, 2005, 33(3): 321-325.

    [22] 徐祖信. 我國河流綜合水質(zhì)標識指數(shù)評價方法研究[J]. 同濟大學學報(自然科學版), 2005, 33(4): 482-488.

    XU Zuxin. Comprehensive water quality identification index for Environmental Quality assessment of surface water[J]. Journal of Tongji University (Natural Science Edition), 2005, 33(4): 482-488.

    [23] YU Tao, MA Jun, ZHANG Liqiu. Factors affecting ice crystal purity during freeze concentration process for urine treatment[J]. Journal of Harbin Institute of Technology, 2007, 14(5): 593-597.

    [24] ABDELMOATY H M, MAHGOUB A U, ABDELDAYEM A W. Performance analysis of salt reduction levels in indirect freeze desalination system with and without magnetic field exposure[J]. Desalination, 2021, 508: 115 021.

    [25] 楊文煥, 杜璇, 姚植, 等. 基于VAR模型的烏梁素海富營養(yǎng)化影響因子評價[J]. 環(huán)境科學與技術(shù), 2021, 44(2): 25-32.

    YANG Wenhuan, DU Xuan, YAO Zhi, et al. Evaluation of the influencing factors of wuliangsuhai eutrophication based on VAR model[J]. Environmental Science & Technology, 2021, 44(2): 25-32.

    Change in Water Quality of Ulansuhai Lake after Its Surface Freezing

    ZHANG Yan1,WANG Nan1, ZHAO Tongguo1, LI Weiping2, TIAN Yepeng3, YUAN Lingfei4

    (1. Civil Engineering College, Yantai University, Yantai 264005, China; 2. School of Energy and Environment, Inner Mongolia University of Science and Technology, Baotou 014010, China; 3. Trier Institute of Technology for Sustainable Development,Yantai University, Yantai 264005, China; 4. Yantai Taoziwan Wastewater Treatment Co, Yantai 264013, China)

    【Objective】Chemical quality of surface water bodies is affected by many factors but most of them are driven by oxygen. Bioavailable oxygen in water decreases after the water surface is frozen, and this paper investigates how this impacts water quality in Ulansuhai lake in Inner Mongolia.【Method】Ice and water samples were taken in January 2021 from 10 sites in the lake. Contents of CODMn, nitrate, total phosphorus (TP), total nitrogen (TN), fluoride and copper in each sample were analyzed. They were compared to those measured before the lake was frozen. Change in these chemicals due to freezing was assessed using the single-factor and the comprehensive water quality identification index, as well as the principal component analysis.【Result】①Single factor water quality identification index analysis indicated that the content of all chemicals in fluid water was ranked in the order of water underneath ice≥water before freezing≥melted ice. The pollution degree of the six chemicals was ranked in the order of TN>TP>CODMn>NH3-N>copper>fluoride. The chemicals that resulted in water quality classification change were CODMn, NH3-N and TN. ② Comprehensive water quality identification index analysis showed that freezing worsened the water quality in all 10 sampling sites, though the deterioration varied with sampling sites. Nitrate, TN and TP had the greatest influence on water quality, as increasing nitrogen and phosphorus risked eutrophication. Spatially, water quality in the lake area was ranked in the order of outlet area>middle lake area>inlet area.【Conclusion】Surface freezing in Ulansuhai Lake increased chemical concentrations in the water beneath the ice, leading to water quality deterioration. To ameliorate water pollution, discharge of pollutants to the lake, especially for TN and TP in autumn and winter, should be reduced.

    Ulansuhai Lake; icing effect; water quality evaluation; principal component analysis method; water quality identification index

    張巖, 王楠, 趙同國, 等. 冰封對烏梁素海水質(zhì)污染特征的影響[J]. 灌溉排水學報, 2023, 42(3): 120-127.

    ZHANG Yan, WANG Nan, ZHAO Tongguo, et al. Change in Water Quality of Ulansuhai Lake after Its Surface Freezing[J]. Journal of Irrigation and Drainage, 2023, 42(3): 120-127.

    1672 - 3317(2023)03 - 0120 - 08

    X524

    A

    10.13522/j.cnki.ggps.2022384

    2022-07-13

    山東省優(yōu)秀中青年科學家科研獎勵基金項目(BS2014HZ021);國家自然科學基金項目(51609207);內(nèi)蒙古自治區(qū)自然科學基金項目(2019MS05006)

    張巖(1984-),男。教授,主要研究方向為水資源利用與保護。E-mail: zhangyan-992@163.com

    責任編輯:白芳芳

    猜你喜歡
    烏梁素海冰封類水
    神秘的冰封之灣
    感受極地冰封之美
    冰封奇觀
    奇聞怪事(2021年2期)2021-04-20 12:12:01
    冰封在“空中”的湖魚
    浙江麗水:打好劣V類水剿滅戰(zhàn)
    十字路口的烏梁素海
    中國三峽(2017年4期)2017-06-06 10:44:22
    熒光類水滑石的細菌檢測及殺菌性能
    烏梁素海地區(qū)近50年氣候變化特征分析
    十四烷酸插層稀土類水滑石的合成及其對PVC的熱穩(wěn)定作用
    中國塑料(2015年6期)2015-11-13 03:03:11
    PA6/類水滑石納米復合材料的制備與性能研究
    中國塑料(2015年5期)2015-10-14 00:59:49
    色94色欧美一区二区| 国产成人91sexporn| 欧美国产精品一级二级三级| 亚洲欧美清纯卡通| 国产深夜福利视频在线观看| 亚洲精品美女久久av网站| 精品人妻一区二区三区麻豆| 蜜桃在线观看..| 成人综合一区亚洲| videosex国产| 久久免费观看电影| 久久人妻熟女aⅴ| 永久免费av网站大全| 伊人亚洲综合成人网| 中文字幕另类日韩欧美亚洲嫩草| 国产精品 国内视频| 男女免费视频国产| 免费在线观看黄色视频的| 亚洲一级一片aⅴ在线观看| 男女边吃奶边做爰视频| 国产免费视频播放在线视频| 亚洲欧美中文字幕日韩二区| 久久精品国产鲁丝片午夜精品| 国产日韩欧美在线精品| 国产一区二区三区综合在线观看 | 国产亚洲最大av| 午夜久久久在线观看| 午夜免费鲁丝| 欧美日韩av久久| 久热久热在线精品观看| 国产又色又爽无遮挡免| 狠狠婷婷综合久久久久久88av| 母亲3免费完整高清在线观看 | 考比视频在线观看| 看十八女毛片水多多多| 久久久久国产网址| 午夜日本视频在线| 亚洲伊人久久精品综合| 国产日韩一区二区三区精品不卡| 亚洲欧美清纯卡通| 热re99久久精品国产66热6| 极品少妇高潮喷水抽搐| 日日啪夜夜爽| 熟女av电影| 老司机影院毛片| 久久久亚洲精品成人影院| 看非洲黑人一级黄片| 日韩人妻精品一区2区三区| 69精品国产乱码久久久| 男女高潮啪啪啪动态图| 亚洲精品国产av成人精品| 欧美精品亚洲一区二区| 99久久人妻综合| 18禁国产床啪视频网站| 成人毛片60女人毛片免费| 成人黄色视频免费在线看| 成人国语在线视频| 久久午夜综合久久蜜桃| 观看美女的网站| 一二三四中文在线观看免费高清| 高清不卡的av网站| 侵犯人妻中文字幕一二三四区| 女的被弄到高潮叫床怎么办| 国产又色又爽无遮挡免| 边亲边吃奶的免费视频| 少妇的丰满在线观看| 亚洲情色 制服丝袜| 亚洲欧美清纯卡通| 三上悠亚av全集在线观看| 美女中出高潮动态图| 在线天堂中文资源库| 亚洲精品aⅴ在线观看| av免费观看日本| 又粗又硬又长又爽又黄的视频| 日本爱情动作片www.在线观看| 婷婷色综合www| 超碰97精品在线观看| 在线观看人妻少妇| 国产一区二区激情短视频 | 免费黄网站久久成人精品| 国产又爽黄色视频| 日本免费在线观看一区| 国产麻豆69| 免费av中文字幕在线| 久久97久久精品| 欧美日韩av久久| 熟妇人妻不卡中文字幕| 久久这里有精品视频免费| 在线精品无人区一区二区三| 国产福利在线免费观看视频| 日本wwww免费看| 五月开心婷婷网| 亚洲精品乱码久久久久久按摩| 日韩不卡一区二区三区视频在线| 国国产精品蜜臀av免费| 国产精品嫩草影院av在线观看| 晚上一个人看的免费电影| 80岁老熟妇乱子伦牲交| 亚洲精品美女久久久久99蜜臀 | 男的添女的下面高潮视频| 香蕉国产在线看| 又粗又硬又长又爽又黄的视频| 国产精品久久久久久av不卡| 久久 成人 亚洲| av国产精品久久久久影院| 久久久精品区二区三区| 国产成人精品在线电影| 大片电影免费在线观看免费| 麻豆精品久久久久久蜜桃| 高清视频免费观看一区二区| 国产麻豆69| 亚洲欧美清纯卡通| 中国美白少妇内射xxxbb| 街头女战士在线观看网站| 免费观看无遮挡的男女| 精品亚洲成国产av| 午夜福利网站1000一区二区三区| 免费av中文字幕在线| 激情视频va一区二区三区| 国产探花极品一区二区| 色吧在线观看| 精品一区二区三区视频在线| 亚洲,欧美,日韩| 汤姆久久久久久久影院中文字幕| 一二三四中文在线观看免费高清| 亚洲国产看品久久| 久久99蜜桃精品久久| 啦啦啦中文免费视频观看日本| 亚洲精品国产av蜜桃| 大陆偷拍与自拍| 国产精品人妻久久久影院| 9热在线视频观看99| 中国美白少妇内射xxxbb| 美女主播在线视频| 日本欧美国产在线视频| 午夜福利,免费看| 中文字幕人妻熟女乱码| 老司机亚洲免费影院| 成人午夜精彩视频在线观看| 看非洲黑人一级黄片| 一本久久精品| 国产不卡av网站在线观看| 国产精品一国产av| 我要看黄色一级片免费的| 成年女人在线观看亚洲视频| 69精品国产乱码久久久| 曰老女人黄片| 欧美+日韩+精品| 亚洲国产av影院在线观看| 欧美人与性动交α欧美软件 | 大片电影免费在线观看免费| 午夜精品国产一区二区电影| 天天操日日干夜夜撸| 制服人妻中文乱码| 在线亚洲精品国产二区图片欧美| 插逼视频在线观看| 成人漫画全彩无遮挡| 纯流量卡能插随身wifi吗| 免费女性裸体啪啪无遮挡网站| 黑人猛操日本美女一级片| 飞空精品影院首页| 22中文网久久字幕| 精品国产国语对白av| 男女无遮挡免费网站观看| 国产极品天堂在线| 国语对白做爰xxxⅹ性视频网站| 精品一品国产午夜福利视频| 日本av免费视频播放| 国产成人a∨麻豆精品| 又黄又爽又刺激的免费视频.| 日韩av不卡免费在线播放| 在线观看美女被高潮喷水网站| 伊人久久国产一区二区| 久久久久久久久久久久大奶| 午夜福利视频精品| 国产精品国产av在线观看| 一区在线观看完整版| 国产免费一区二区三区四区乱码| 国产精品国产av在线观看| 最新中文字幕久久久久| 欧美xxxx性猛交bbbb| 日本黄色日本黄色录像| 亚洲内射少妇av| 国产精品久久久久久精品古装| 黄色一级大片看看| 国产欧美另类精品又又久久亚洲欧美| 婷婷成人精品国产| 亚洲精品aⅴ在线观看| 亚洲精品久久午夜乱码| 免费人成在线观看视频色| 国产精品成人在线| 午夜福利影视在线免费观看| 国产午夜精品一二区理论片| 亚洲综合色惰| 亚洲国产欧美日韩在线播放| 插逼视频在线观看| www.熟女人妻精品国产 | 成人国产av品久久久| xxxhd国产人妻xxx| 亚洲丝袜综合中文字幕| 国产麻豆69| 少妇人妻精品综合一区二区| 欧美日韩综合久久久久久| 多毛熟女@视频| 欧美日韩亚洲高清精品| a级毛片在线看网站| 一区二区三区四区激情视频| 久久精品aⅴ一区二区三区四区 | 日韩av免费高清视频| 18禁动态无遮挡网站| 插逼视频在线观看| 亚洲精品国产色婷婷电影| 久久久久久人妻| 永久免费av网站大全| 高清不卡的av网站| 国产毛片在线视频| 2018国产大陆天天弄谢| 日日啪夜夜爽| 国产一区二区在线观看日韩| 97精品久久久久久久久久精品| 国产一区二区激情短视频 | 欧美精品高潮呻吟av久久| 成人毛片a级毛片在线播放| 在线观看美女被高潮喷水网站| 久久久久久久亚洲中文字幕| 搡女人真爽免费视频火全软件| 在现免费观看毛片| 捣出白浆h1v1| 色哟哟·www| 18禁观看日本| 22中文网久久字幕| 黄色毛片三级朝国网站| 日韩成人伦理影院| 一本—道久久a久久精品蜜桃钙片| 麻豆精品久久久久久蜜桃| 一级毛片黄色毛片免费观看视频| 久久青草综合色| av在线app专区| 亚洲精品久久久久久婷婷小说| 国产在视频线精品| 国产一区二区激情短视频 | 亚洲人与动物交配视频| 中文字幕另类日韩欧美亚洲嫩草| 少妇猛男粗大的猛烈进出视频| 少妇被粗大的猛进出69影院 | 亚洲精品中文字幕在线视频| 天天操日日干夜夜撸| 男男h啪啪无遮挡| 97超碰精品成人国产| 色婷婷av一区二区三区视频| 中文字幕亚洲精品专区| 欧美日韩成人在线一区二区| 国产老妇伦熟女老妇高清| 国产在线视频一区二区| 两个人看的免费小视频| 午夜免费观看性视频| 亚洲美女视频黄频| 久久久久久人人人人人| 如日韩欧美国产精品一区二区三区| 久久国内精品自在自线图片| 国产免费现黄频在线看| 久久97久久精品| 久久毛片免费看一区二区三区| 日本vs欧美在线观看视频| 妹子高潮喷水视频| 午夜日本视频在线| 九色亚洲精品在线播放| 国产免费福利视频在线观看| 丰满少妇做爰视频| 亚洲av电影在线进入| 国产日韩一区二区三区精品不卡| 赤兔流量卡办理| 亚洲精品色激情综合| 男人舔女人的私密视频| 日韩中字成人| 国产无遮挡羞羞视频在线观看| 交换朋友夫妻互换小说| 在线 av 中文字幕| 日本免费在线观看一区| 精品久久久精品久久久| 国产免费现黄频在线看| 国产1区2区3区精品| 国产成人午夜福利电影在线观看| 大片电影免费在线观看免费| 国产日韩欧美在线精品| 日韩三级伦理在线观看| 久久精品国产鲁丝片午夜精品| 如日韩欧美国产精品一区二区三区| 九色亚洲精品在线播放| 美女大奶头黄色视频| 中文字幕av电影在线播放| 亚洲人成网站在线观看播放| 多毛熟女@视频| 欧美精品高潮呻吟av久久| 少妇人妻 视频| 建设人人有责人人尽责人人享有的| 天天操日日干夜夜撸| av电影中文网址| 国产欧美亚洲国产| 国产av国产精品国产| 天天影视国产精品| 亚洲一级一片aⅴ在线观看| 熟妇人妻不卡中文字幕| 日韩精品有码人妻一区| 飞空精品影院首页| 日本午夜av视频| 侵犯人妻中文字幕一二三四区| 满18在线观看网站| 中文乱码字字幕精品一区二区三区| 亚洲经典国产精华液单| 一二三四中文在线观看免费高清| 国产一区二区激情短视频 | 蜜桃国产av成人99| 欧美日韩视频精品一区| 热re99久久精品国产66热6| 国产成人精品久久久久久| 咕卡用的链子| 精品国产露脸久久av麻豆| 亚洲av福利一区| 国产亚洲av片在线观看秒播厂| 国产一区二区在线观看日韩| 国产欧美日韩一区二区三区在线| 日本-黄色视频高清免费观看| 热re99久久国产66热| 久久精品国产亚洲av涩爱| av一本久久久久| h视频一区二区三区| 亚洲一级一片aⅴ在线观看| 亚洲图色成人| 日韩不卡一区二区三区视频在线| 黑人巨大精品欧美一区二区蜜桃 | 精品亚洲成国产av| 亚洲精品乱久久久久久| 日本91视频免费播放| 亚洲欧洲国产日韩| 精品99又大又爽又粗少妇毛片| 久久99蜜桃精品久久| 99热全是精品| 2022亚洲国产成人精品| 99热全是精品| av播播在线观看一区| 三级国产精品片| 亚洲av国产av综合av卡| 亚洲精品第二区| 中文天堂在线官网| 下体分泌物呈黄色| 一级,二级,三级黄色视频| 五月天丁香电影| 日本黄大片高清| 热99国产精品久久久久久7| 中文字幕另类日韩欧美亚洲嫩草| 亚洲四区av| av.在线天堂| 亚洲欧美清纯卡通| 亚洲国产精品一区二区三区在线| av女优亚洲男人天堂| 亚洲第一区二区三区不卡| 亚洲国产最新在线播放| 你懂的网址亚洲精品在线观看| 18禁在线无遮挡免费观看视频| 大香蕉久久网| 久久精品国产综合久久久 | 人妻少妇偷人精品九色| 亚洲一区二区三区欧美精品| 99热6这里只有精品| 2022亚洲国产成人精品| 婷婷色av中文字幕| 日本午夜av视频| 欧美成人精品欧美一级黄| 日本色播在线视频| 18禁动态无遮挡网站| 少妇被粗大猛烈的视频| 久久久久久人妻| 少妇人妻 视频| 老熟女久久久| 少妇的丰满在线观看| 一级爰片在线观看| 久久人人爽人人爽人人片va| 久久婷婷青草| 制服人妻中文乱码| 捣出白浆h1v1| 亚洲婷婷狠狠爱综合网| 中文字幕另类日韩欧美亚洲嫩草| 国产片特级美女逼逼视频| 成人毛片a级毛片在线播放| 欧美精品亚洲一区二区| 一区二区三区精品91| 国产日韩欧美视频二区| 91国产中文字幕| 国产成人欧美| 最新中文字幕久久久久| 丁香六月天网| 日韩av免费高清视频| 成人午夜精彩视频在线观看| 欧美性感艳星| 国产在线免费精品| 久久精品国产综合久久久 | 一二三四在线观看免费中文在 | 日本欧美国产在线视频| a级毛片在线看网站| 秋霞在线观看毛片| 欧美人与性动交α欧美精品济南到 | 夜夜骑夜夜射夜夜干| 国产有黄有色有爽视频| 老司机亚洲免费影院| 在线看a的网站| www日本在线高清视频| 国产福利在线免费观看视频| 一级黄片播放器| 日产精品乱码卡一卡2卡三| 日韩一区二区三区影片| 日本色播在线视频| 伊人久久国产一区二区| 久久人人爽人人爽人人片va| 亚洲欧美成人精品一区二区| 人人妻人人添人人爽欧美一区卜| 欧美激情极品国产一区二区三区 | 一级片免费观看大全| 免费日韩欧美在线观看| 日韩人妻精品一区2区三区| 日韩一区二区三区影片| 欧美bdsm另类| 极品少妇高潮喷水抽搐| 日本猛色少妇xxxxx猛交久久| 国产精品国产av在线观看| 国产精品女同一区二区软件| 黄色一级大片看看| 婷婷色综合www| 人成视频在线观看免费观看| 亚洲精品成人av观看孕妇| 欧美日韩成人在线一区二区| 狠狠精品人妻久久久久久综合| 国产成人91sexporn| 巨乳人妻的诱惑在线观看| 欧美另类一区| 欧美+日韩+精品| 国产免费视频播放在线视频| 欧美xxⅹ黑人| 麻豆乱淫一区二区| 久久久久精品久久久久真实原创| 欧美精品高潮呻吟av久久| 色哟哟·www| 日韩成人伦理影院| 精品少妇内射三级| 成人毛片60女人毛片免费| 99热全是精品| 亚洲欧洲国产日韩| 精品国产一区二区三区四区第35| 久久热在线av| 性高湖久久久久久久久免费观看| 最近手机中文字幕大全| 国产熟女欧美一区二区| 久久韩国三级中文字幕| 久久久久精品人妻al黑| 亚洲欧美一区二区三区黑人 | 黄色一级大片看看| 精品一区二区三区视频在线| 人妻少妇偷人精品九色| 亚洲美女搞黄在线观看| 26uuu在线亚洲综合色| 在线亚洲精品国产二区图片欧美| 69精品国产乱码久久久| 国产精品久久久久久精品古装| 一区二区三区精品91| 亚洲精品成人av观看孕妇| av在线观看视频网站免费| 国产成人aa在线观看| 亚洲国产日韩一区二区| 久久 成人 亚洲| 国产免费现黄频在线看| 99九九在线精品视频| 成人亚洲欧美一区二区av| 日本爱情动作片www.在线观看| 精品人妻熟女毛片av久久网站| 亚洲欧洲精品一区二区精品久久久 | 熟女av电影| av又黄又爽大尺度在线免费看| 免费少妇av软件| 亚洲综合色网址| 日韩电影二区| 日韩av免费高清视频| 日韩电影二区| 男女边吃奶边做爰视频| 日韩精品免费视频一区二区三区 | 美女视频免费永久观看网站| 两性夫妻黄色片 | 欧美日韩一区二区视频在线观看视频在线| 夜夜骑夜夜射夜夜干| 亚洲第一区二区三区不卡| 日韩欧美精品免费久久| 黄色配什么色好看| 天美传媒精品一区二区| 少妇高潮的动态图| 久久97久久精品| 成人免费观看视频高清| 91精品国产国语对白视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av在线观看视频网站免费| 久久久精品区二区三区| 大香蕉久久成人网| 国产精品99久久99久久久不卡 | 在线观看一区二区三区激情| 精品一区二区三区四区五区乱码 | 久久精品久久久久久久性| 亚洲精品,欧美精品| 国产日韩欧美视频二区| 蜜臀久久99精品久久宅男| 亚洲精品第二区| 欧美日韩成人在线一区二区| 激情视频va一区二区三区| 黄片无遮挡物在线观看| 国产乱来视频区| 久久99蜜桃精品久久| 亚洲av综合色区一区| av免费观看日本| 18禁动态无遮挡网站| 久久精品久久久久久噜噜老黄| 久久人人97超碰香蕉20202| 亚洲国产av新网站| videos熟女内射| 国产麻豆69| 老女人水多毛片| 亚洲综合色惰| 国产成人精品一,二区| 国产亚洲欧美精品永久| 人妻 亚洲 视频| 国产亚洲最大av| 高清视频免费观看一区二区| 久久99热这里只频精品6学生| 亚洲国产日韩一区二区| 91精品国产国语对白视频| 久久久国产精品麻豆| 久久人人爽av亚洲精品天堂| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 不卡视频在线观看欧美| 国产av国产精品国产| 日韩熟女老妇一区二区性免费视频| 成人影院久久| 日韩成人伦理影院| 亚洲av电影在线观看一区二区三区| 国产毛片在线视频| 亚洲,一卡二卡三卡| 如日韩欧美国产精品一区二区三区| 美女视频免费永久观看网站| 久久久久精品久久久久真实原创| 黑人巨大精品欧美一区二区蜜桃 | 新久久久久国产一级毛片| 大码成人一级视频| 青春草视频在线免费观看| 午夜免费鲁丝| 99re6热这里在线精品视频| 黄色毛片三级朝国网站| 色网站视频免费| 婷婷色综合www| 成人毛片a级毛片在线播放| 99热全是精品| av免费观看日本| 亚洲少妇的诱惑av| 久久久国产一区二区| 久久精品aⅴ一区二区三区四区 | 高清av免费在线| 久久亚洲国产成人精品v| 亚洲欧美日韩卡通动漫| 老司机影院毛片| 大香蕉97超碰在线| 亚洲一区二区三区欧美精品| 999精品在线视频| 日日撸夜夜添| 亚洲丝袜综合中文字幕| 久久人人爽av亚洲精品天堂| 久久久久久久久久久免费av| av天堂久久9| 国产在线免费精品| 9191精品国产免费久久| 日本免费在线观看一区| 亚洲人与动物交配视频| 26uuu在线亚洲综合色| 亚洲国产欧美日韩在线播放| 七月丁香在线播放| 亚洲国产精品一区三区| av有码第一页| 春色校园在线视频观看| 91午夜精品亚洲一区二区三区| 欧美日韩国产mv在线观看视频| 制服丝袜香蕉在线| 男女边摸边吃奶| 97人妻天天添夜夜摸| 99视频精品全部免费 在线| 18禁在线无遮挡免费观看视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 看免费成人av毛片| 制服人妻中文乱码| 欧美激情极品国产一区二区三区 | 最近中文字幕2019免费版| 日本欧美国产在线视频| 亚洲第一区二区三区不卡| 欧美xxxx性猛交bbbb| 欧美日韩精品成人综合77777| 中文字幕亚洲精品专区| 国产精品国产三级国产av玫瑰| 精品亚洲成国产av| 精品久久国产蜜桃| 久久久久精品人妻al黑| 一本久久精品| 亚洲av中文av极速乱| 青春草国产在线视频| 一区二区av电影网| 欧美97在线视频| 在线观看美女被高潮喷水网站| 欧美最新免费一区二区三区| 美女国产高潮福利片在线看| 亚洲美女黄色视频免费看| 涩涩av久久男人的天堂| 国产国拍精品亚洲av在线观看| 午夜福利,免费看|