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

    基于EFDC模型的河道型水庫(kù)藻類(lèi)生長(zhǎng)對(duì)流域污染負(fù)荷削減的響應(yīng)*
    ——以廣東長(zhǎng)潭水庫(kù)為例

    2015-05-10 01:34:54李一平王靜雨
    湖泊科學(xué) 2015年5期
    關(guān)鍵詞:點(diǎn)位藻類(lèi)葉綠素

    李一平,王靜雨,滑 磊

    (1:河海大學(xué)環(huán)境學(xué)院,南京 210098)

    (2:浙江省水利河口研究院,杭州 310020)

    基于EFDC模型的河道型水庫(kù)藻類(lèi)生長(zhǎng)對(duì)流域污染負(fù)荷削減的響應(yīng)*
    ——以廣東長(zhǎng)潭水庫(kù)為例

    李一平1,王靜雨1**,滑 磊2

    (1:河海大學(xué)環(huán)境學(xué)院,南京 210098)

    (2:浙江省水利河口研究院,杭州 310020)

    將環(huán)境(EFDC)模型應(yīng)用于廣東梅州的河道型水庫(kù)——長(zhǎng)潭水庫(kù),經(jīng)過(guò)模型的驗(yàn)證,建立長(zhǎng)潭水庫(kù)水生態(tài)數(shù)學(xué)模型,模擬出水庫(kù)在現(xiàn)狀污染負(fù)荷下的藻類(lèi)生長(zhǎng)情況,以葉綠素a濃度為指標(biāo),研究水庫(kù)在豐水年與枯水年富營(yíng)養(yǎng)化改善程度對(duì)不同流域污染削減方案的響應(yīng)關(guān)系.結(jié)果表明:現(xiàn)狀污染負(fù)荷削減10%,在豐水年庫(kù)尾、庫(kù)中和庫(kù)首點(diǎn)位的葉綠素a平均濃度分別降低13.99%、12.00%和10.35%;枯水年3個(gè)代表點(diǎn)位葉綠素a平均濃度分別下降8.42%、5.63%和2.10%.污染削減20%的情況下,豐水年3個(gè)代表點(diǎn)位的葉綠素a濃度分別降低26.78%、19.25%和17.04%,枯水年對(duì)應(yīng)的降低幅度分別為11.72%、7.97%和5.12%;庫(kù)尾地區(qū)表現(xiàn)出河道的特征,在污染負(fù)荷削減的情況下藻類(lèi)生長(zhǎng)能夠受到有效遏制,庫(kù)首地區(qū)則表現(xiàn)出湖泊的特性,枯水年這一特征表現(xiàn)更加突出,污染削減20%的情況下,葉綠素a的平均濃度僅降低5.12%.

    河道型水庫(kù);長(zhǎng)潭水庫(kù);EFDC模型;污染削減;藻類(lèi)生長(zhǎng)響應(yīng)

    隨著我國(guó)三峽水庫(kù)、小浪底水庫(kù)等水利工程的修建完成,對(duì)于水庫(kù)的保護(hù)與管理工作也日益迫切.河道型水庫(kù)形成后,其水動(dòng)力學(xué)特征介于湖泊和河道之間[1],水體滯留時(shí)間變長(zhǎng),水流變緩慢,比河流更容易發(fā)生富營(yíng)養(yǎng)化[1],加之水庫(kù)與河流的不可分割性,水庫(kù)水質(zhì)變差必然影響到下游河道的水質(zhì).因此需要研究河道型水庫(kù)的富營(yíng)養(yǎng)化特征與控制對(duì)策.

    近年來(lái),EFDC模型被廣泛地應(yīng)用于湖泊、水庫(kù)、河流、海灣等地表水的模擬與分析,可以定量模擬環(huán)境特征、污染負(fù)荷與湖泊水質(zhì)間的動(dòng)態(tài)響應(yīng)關(guān)系,亦可用作湖泊、河道等水體富營(yíng)養(yǎng)化控制的決策工具.李一平等利用EFDC模型建立太湖的水環(huán)境動(dòng)力學(xué)模型,研究“引江濟(jì)太”對(duì)于太湖水動(dòng)力的調(diào)控效果[2-3].衛(wèi)志宏等利用EFDC模型建立洱海湖泊及湖灣的三維水動(dòng)力水質(zhì)模型,研究洱海主要污染物允許排放量的控制分配[4].李杰等利用EFDC模型,研究美國(guó)北卡羅來(lái)納州Falls Lake水庫(kù)內(nèi)物質(zhì)輸運(yùn)過(guò)程,模擬結(jié)果顯示水庫(kù)結(jié)構(gòu)和徑流分布是影響污染物滯留時(shí)間的重要因素[5].陳異暉利用EFDC模型建立了滇池水質(zhì)模型,研究了人工調(diào)水工程對(duì)于滇池的水質(zhì)改善效果,但由于基礎(chǔ)資料不足,建立的模型對(duì)于滇池水質(zhì)模擬并不夠精確[6].李一平等利用EFDC模型研究了河道型水庫(kù)水動(dòng)力特征與氣候條件的響應(yīng)關(guān)系[7],相關(guān)文獻(xiàn)報(bào)道[2-4]顯示,EFDC模型具有通用性好、數(shù)值計(jì)算能力強(qiáng)、數(shù)據(jù)輸出應(yīng)用范圍廣等特點(diǎn),尤其對(duì)于水體水動(dòng)力的模擬精度已達(dá)到相當(dāng)高的水平.但目前國(guó)內(nèi)還沒(méi)有關(guān)于河道型水庫(kù)藻類(lèi)生長(zhǎng)與污染負(fù)荷削減響應(yīng)關(guān)系的研究.針對(duì)上述問(wèn)題,本文選用環(huán)境流體動(dòng)力學(xué)模型,以葉綠素a為指標(biāo),來(lái)模擬不同污染削減方案下長(zhǎng)潭水庫(kù)藻類(lèi)生長(zhǎng)情況,利用EFDC健全的模擬結(jié)果可視化功能及后處理功能,研究分析水庫(kù)不同位置藻類(lèi)生長(zhǎng)的特征,探討河道型水庫(kù)控制藻類(lèi)生長(zhǎng)的對(duì)策.

    1 長(zhǎng)潭水庫(kù)水環(huán)境數(shù)學(xué)模型構(gòu)建

    1.1 EFDC模型簡(jiǎn)介

    長(zhǎng)潭水庫(kù)水動(dòng)力學(xué)模型采用EFDC模型構(gòu)建,EFDC是在美國(guó)國(guó)家環(huán)保署資助下由維吉尼亞海洋研究所(VIMS, Virginia Institute of Marine Science at the College of William and Mary)的John Hamrick等根據(jù)多個(gè)數(shù)學(xué)模型集成開(kāi)發(fā)研制的綜合模型,被用于模擬水系統(tǒng)一維、二維和三維流場(chǎng),物質(zhì)輸運(yùn)(包括溫度、鹽度和泥沙的輸運(yùn)),生態(tài)過(guò)程以及淡水入流等.與POM、ECOM模型相比,EFDC模型的優(yōu)勢(shì)為:靈活的邊界處理技術(shù),通用的文件輸入格式,能快速耦合水動(dòng)力、泥沙和水質(zhì)模塊,省略了不同模型接口程序的研發(fā)過(guò)程.同時(shí)EFDC開(kāi)發(fā)有完整的前、后處理軟件EFDC-Explorer,采用可視化的界面操作,能快速生成網(wǎng)格數(shù)據(jù)和處理圖像文件.EFDC自1992年發(fā)布以來(lái),不斷更新完善,隨著EFDC-Explorer等前、后處理工具的出現(xiàn),EFDC得到了更加廣泛的應(yīng)用.

    EFDC模型垂向上采用σ坐標(biāo)變換,能較好地?cái)M和近岸復(fù)雜的岸線和地形;采用Gelperin等修正的Mellor-Yamada 2.5階湍封閉模式較客觀地提供垂向混合系數(shù),避免其人為選取造成的誤差.動(dòng)量方程、連續(xù)方程及狀態(tài)方程為:

    ?t(mxmyHu)+?x(myHuu)+?y(mxHvu)+?z(mxmywu)-femxmyHv=

    (1)

    ?t(mxmyHv)+?x(myHuv)+?y(mxHvv)+?z(mxmywv)+femxmyHu=

    (2)

    mxmyfe=mxmyf-u?ymx+v?xmy

    (3)

    水質(zhì)變量的質(zhì)量守恒控制方程為:

    (4)

    式中,u、v、w分別為邊界擬合正交曲線坐標(biāo)x、y、z方向上的速度分量;mx和my分別為度量張量對(duì)角元素的平方根,m=mxmy為度量張量行列式的平方根;Av為垂向紊動(dòng)黏滯系數(shù);Qu、Qv分別為動(dòng)量源、匯項(xiàng);C為水質(zhì)變量濃度;Kx、Ky和Kz分別為x、y、z方向的擴(kuò)散系數(shù);SC為單位體積源匯項(xiàng).

    水質(zhì)變量的質(zhì)量守恒方程包括物理輸運(yùn)、平流擴(kuò)散以及生態(tài)動(dòng)力學(xué)過(guò)程,式(4)等號(hào)左邊后3項(xiàng)為平流輸運(yùn)項(xiàng),等號(hào)右邊前3項(xiàng)為擴(kuò)散輸運(yùn)項(xiàng).這6項(xiàng)和物理輸運(yùn)項(xiàng)類(lèi)似,數(shù)值解法和水動(dòng)力模型中鹽度質(zhì)量守恒方程相似.等號(hào)右邊最后1項(xiàng)表示每個(gè)狀態(tài)變量的動(dòng)力學(xué)過(guò)程和外部負(fù)荷.

    圖1 長(zhǎng)潭水庫(kù)水動(dòng)力模型(底部高程)以及代表點(diǎn)位置Fig.1 The bottom elevation and representative points of Changtan Reservoir water dynamic model

    圖2 長(zhǎng)潭水庫(kù)水動(dòng)力模型(網(wǎng)格局部放大圖)Fig.2 Enlarged view of local grids of Changtan Reservoir

    藻類(lèi)是模型模擬的一個(gè)重要變量,模型考慮藍(lán)藻、硅藻和綠藻3種藻類(lèi).下標(biāo)x表示藻的種類(lèi).模型考慮的源匯項(xiàng)包括:生長(zhǎng)、基礎(chǔ)代謝、捕食、沉降以及外部負(fù)荷.描述這3種藻的動(dòng)力學(xué)方程基本相同,只是方程中的參數(shù)取值不同.動(dòng)力學(xué)方程為:

    (5)

    式中,Bx為x種類(lèi)藻的生物量(g/m3),Px為x種類(lèi)藻的生產(chǎn)率(d-1),BMx為x種類(lèi)藻的基礎(chǔ)代謝率(d-1),PRx為x種類(lèi)藻的被捕食率(d-1),WSx為x種類(lèi)藻的沉降速率(m/d),WBx為x種類(lèi)藻的外部負(fù)荷(g/d),V為模擬單元.

    1.2 研究區(qū)域概況

    長(zhǎng)潭水庫(kù)(24°42′2″~24°50′15″N,116°4′8″~116°8′1″E)位于廣東省梅州蕉嶺縣蕉城鎮(zhèn)西北約6km的石窟河長(zhǎng)潭峽谷中,是以發(fā)電、供水、防洪為主要功能的大(二)型水庫(kù).長(zhǎng)潭水庫(kù)長(zhǎng)興電站壩址至長(zhǎng)潭水庫(kù)壩址之間全長(zhǎng)約22km,寬100~350m,平均水深30m,庫(kù)容1.4×108m3,水面面積4.3km2,屬于典型的河道型水庫(kù).

    長(zhǎng)潭水庫(kù)入庫(kù)河流有石窟河干流(中山河)、中赤河、查干河.石窟河是韓江一級(jí)支流,又名石窟溪、蕉嶺河,源于福建武平縣東留,流經(jīng)平遠(yuǎn)、蕉嶺、梅縣等縣市,于梅縣丙村鎮(zhèn)東州壩注入梅江.石窟河干流長(zhǎng)178km,流域面積3777km2.

    長(zhǎng)潭水庫(kù)于1978年3月開(kāi)始建設(shè),1987年開(kāi)始蓄水發(fā)電,1991年3月竣工.水庫(kù)設(shè)計(jì)洪水位為151.50m,正常庫(kù)容為1.145×108m3,死庫(kù)容為0.5985×108m3,設(shè)計(jì)正常高水位為148.00m,汛期防洪限制水位為144.00m,死水位為136.50m,屬季節(jié)調(diào)節(jié)水庫(kù).長(zhǎng)潭水庫(kù)多年平均流量為55.59m3/s,多年平均徑流量為17.85×108m3,多年平均徑流系數(shù)為0.52.

    目前,廣東省梅州市城區(qū)的供水水源來(lái)自清涼山水庫(kù),根據(jù)梅州城區(qū)需水預(yù)測(cè)結(jié)果,到2020年,梅州城區(qū)需水總量將達(dá)9.305×107m3,大大超過(guò)清涼山水庫(kù)供水工程現(xiàn)有的年供水能力,顯然不能滿足梅州城區(qū)的需水要求.因此,長(zhǎng)潭水庫(kù)作為梅州城區(qū)的后備水源地成為必然的趨勢(shì).

    然而,據(jù)當(dāng)?shù)赜嘘P(guān)部門(mén)反映,長(zhǎng)潭水庫(kù)一年中發(fā)生數(shù)次類(lèi)似藍(lán)藻水華的事件.根據(jù)水環(huán)境調(diào)查情況[8],長(zhǎng)潭水庫(kù)上游及閩粵省界水體富營(yíng)養(yǎng)化污染較重,浮游植物以藍(lán)、綠藻為優(yōu)勢(shì)種,長(zhǎng)潭水庫(kù)總氮、總磷超標(biāo)嚴(yán)重,目前水庫(kù)水質(zhì)為Ⅳ~Ⅴ類(lèi)水,不滿足集中式生活飲用水地表水源地水質(zhì)所需要滿足的Ⅲ類(lèi)水(GB 3838-2002標(biāo)準(zhǔn)規(guī)定)的要求.水庫(kù)主要是受到上游福建省境內(nèi)大量養(yǎng)殖廢水和農(nóng)業(yè)、生活面源污染,為了保證后備水源地的水質(zhì)安全,治理污染勢(shì)在必行,因此有必要開(kāi)展長(zhǎng)潭水庫(kù)的水體富營(yíng)養(yǎng)化與污染負(fù)荷削減之間的響應(yīng)關(guān)系研究,為制定科學(xué)合理的治污措施提供參考.

    1.3 模型的構(gòu)建

    1.3.1 網(wǎng)格生成 EFDC是基于有限差分求解水動(dòng)力方程的數(shù)值模擬系統(tǒng),因此本研究需要進(jìn)行網(wǎng)格概化.模型上邊界始于長(zhǎng)興電站,下邊界為長(zhǎng)潭水庫(kù)大壩壩址,共22km,模型范圍內(nèi)在平面上共劃分為11214個(gè)正方形網(wǎng)格,網(wǎng)格距為20m,模型在垂向上采用σ坐標(biāo),平均分為10層,用庫(kù)底和表層水體厚度來(lái)定義垂向網(wǎng)格的高度,垂向上分為10層,各層所占的水深比例均為0.1.長(zhǎng)潭水庫(kù)水動(dòng)力模型底部高程及局部放大圖見(jiàn)圖1、圖2.

    1.3.2 初始條件 本次模擬的初始條件設(shè)置為模擬時(shí)段第1d長(zhǎng)興電站壩下水位的觀測(cè)值,水質(zhì)的初始參數(shù)設(shè)置為2011年5月長(zhǎng)潭水庫(kù)水質(zhì)監(jiān)測(cè)的實(shí)測(cè)值,水溫設(shè)置為20℃,總磷、總氮、銨態(tài)氮和硝態(tài)氮的濃度分別為0.06、2.17、0.10和1.58mg/L.

    1.3.3 邊界條件 模型邊界條件包括水動(dòng)力邊界和氣象邊界條件,水動(dòng)力邊界條件包括1個(gè)主河道和7條主要入庫(kù)支流,長(zhǎng)興電站的調(diào)度資料構(gòu)成水動(dòng)力上邊界,長(zhǎng)潭電站的水位資料構(gòu)成水動(dòng)力下邊界.氣象邊界條件包括氣壓、氣溫、相對(duì)濕度、降雨量、風(fēng)向、風(fēng)速、云量等逐日數(shù)據(jù),來(lái)源于梅州氣象站.

    1.4 參數(shù)選取與模型驗(yàn)證

    表1為長(zhǎng)潭水庫(kù)水生態(tài)動(dòng)力學(xué)模型中主要參數(shù)的含義與取值,EFDC對(duì)水動(dòng)力的模擬已經(jīng)十分成熟,水動(dòng)力的模擬中,大部分物理參數(shù)都未作改變[6].以長(zhǎng)潭水庫(kù)水位、葉綠素a、總氮、總磷濃度實(shí)測(cè)資料為對(duì)照,利用2011年的水質(zhì)監(jiān)測(cè)數(shù)據(jù)作為該模型的驗(yàn)證資料,比較模型的計(jì)算值與實(shí)測(cè)值之間的誤差.長(zhǎng)潭水庫(kù)水生態(tài)動(dòng)力模型水質(zhì)參數(shù)的驗(yàn)證結(jié)果見(jiàn)圖3.可以看出,EFDC模型模擬出的總氮濃度總體偏高,但與實(shí)測(cè)數(shù)據(jù)的趨勢(shì)一致,模擬的總磷濃度在水庫(kù)中上游地區(qū)偏低,下游擬合程度較好.模擬的葉綠素a濃度誤差在20%以內(nèi),且與實(shí)測(cè)數(shù)據(jù)變化趨勢(shì)一致.整體而言,建立的模型能夠反映出水質(zhì)變化的趨勢(shì),可以用來(lái)模擬預(yù)測(cè)長(zhǎng)潭水庫(kù)藻類(lèi)生長(zhǎng)情況.

    表1 長(zhǎng)潭水庫(kù)水生態(tài)動(dòng)力模型參數(shù)取值

    圖3 長(zhǎng)潭水庫(kù)水位、總氮、總磷和葉綠素a濃度驗(yàn)證圖Fig.3 Verification of water level, total nitrogen, total phosphorus and chlorophyll-a concentrations of Changtan Reservoir

    2 長(zhǎng)潭水庫(kù)藻類(lèi)生長(zhǎng)對(duì)污染負(fù)荷削減的響應(yīng)關(guān)系

    2.1 模擬方案

    根據(jù)長(zhǎng)潭水庫(kù)污染源調(diào)查結(jié)果[8],核算出現(xiàn)狀污染負(fù)荷,長(zhǎng)潭水庫(kù)的氮、磷總負(fù)荷分別為4324.18和65.53t/a,其中長(zhǎng)興電站攜帶的上游流域氮、磷負(fù)荷分別為2563.2和40.05t/a,分別占水庫(kù)污染總負(fù)荷的59.3%和61.1%.依據(jù)《廣東省“十二五”主要污染物總量控制規(guī)劃》制定出兩種污染削減方案,即現(xiàn)狀流域污染負(fù)荷分別削減10%和20%,研究在不同營(yíng)養(yǎng)鹽輸入的邊界條件下水庫(kù)藻類(lèi)生長(zhǎng)的變化情況.

    國(guó)內(nèi)有研究表明,流量增加引起的流速加快對(duì)于水體中藻類(lèi)生長(zhǎng)作用顯著[9-11],因此有必要研究在不同入庫(kù)流量下水動(dòng)力條件對(duì)污染削減方案的配合作用.長(zhǎng)潭水庫(kù)的來(lái)水主要由上游長(zhǎng)興電站的下泄流量和長(zhǎng)潭水庫(kù)流域降雨匯流組成,其中長(zhǎng)興電站上游來(lái)水占水庫(kù)來(lái)水總量的85%以上,根據(jù)長(zhǎng)潭水庫(kù)上游長(zhǎng)興電站的調(diào)度資料逐日變化圖,選取1992和2004年分別作為豐水年和枯水年.相應(yīng)的逐月平均流量見(jiàn)表2.

    表2 長(zhǎng)潭水庫(kù)典型年份逐月入庫(kù)平均流量(m3/s)

    2.2 不同污染削減情景下長(zhǎng)潭水庫(kù)葉綠素a濃度模擬

    根據(jù)長(zhǎng)潭水庫(kù)的調(diào)度資料,查得長(zhǎng)潭水庫(kù)典型年分別為:1992年(豐水年)和2004年(枯水年).由此制定相應(yīng)的2種計(jì)算方案.以2011年污染源現(xiàn)狀調(diào)查結(jié)果作為水質(zhì)模型的邊界條件,以現(xiàn)狀污染負(fù)荷削減10%和20%為原則分別制定出相應(yīng)計(jì)算方案,模擬出長(zhǎng)潭水庫(kù)在現(xiàn)狀污染負(fù)荷的不同削減方案下不同水文年水質(zhì)的改善情況,考慮到藍(lán)藻暴發(fā)多發(fā)生在夏季,制定的方案計(jì)算時(shí)間均為每年的6-8月,共92d.長(zhǎng)潭水庫(kù)兼具河道和湖泊的特性,因此,分別選取庫(kù)尾、庫(kù)中和庫(kù)首3個(gè)代表點(diǎn)位(圖1),定性、定量地研究其藻類(lèi)生長(zhǎng)情況.

    2.3 模擬結(jié)果分析

    藻類(lèi)生長(zhǎng)受氣溫影響,普遍認(rèn)為氣溫在26℃時(shí)藻類(lèi)生物量增長(zhǎng)最快,在水庫(kù)中最易形成水華[12],本文模擬的是夏季長(zhǎng)潭水庫(kù)藻類(lèi)生長(zhǎng)情況,當(dāng)?shù)貧鉁厥诌m宜藻類(lèi)生長(zhǎng),加之進(jìn)入水庫(kù)的氮、磷等污染源超標(biāo),在污染負(fù)荷不削減的情況下,水體內(nèi)藻類(lèi)迅速生長(zhǎng),豐水年長(zhǎng)潭水庫(kù)3個(gè)代表點(diǎn)位的葉綠素a濃度分別在模擬的第18、25和33d達(dá)到峰值,枯水年3個(gè)代表點(diǎn)位的葉綠素a濃度分別在第22、29和30d達(dá)到峰值.

    在污染負(fù)荷削減之后,長(zhǎng)潭水庫(kù)3個(gè)代表點(diǎn)位處總氮、總磷濃度均呈現(xiàn)下降趨勢(shì),庫(kù)尾、庫(kù)中點(diǎn)位水流速度較快,污染物質(zhì)不易停留,總氮、總磷濃度降低的幅度明顯大于庫(kù)首點(diǎn)位.在豐水年,污染負(fù)荷削減10%之后,長(zhǎng)潭水庫(kù)3個(gè)代表點(diǎn)位的總氮濃度降低15%左右,總磷濃度下降幅度更大,流速最快的庫(kù)尾點(diǎn)位達(dá)到了25.24%.在枯水年,污染負(fù)荷削減20%后,庫(kù)首點(diǎn)位總氮、總磷濃度分別為1.50和0.100mg/L(表3、4),仍然面臨富營(yíng)養(yǎng)化的風(fēng)險(xiǎn)[13].

    表3 不同污染負(fù)荷削減方案下長(zhǎng)潭水庫(kù)3個(gè)代表點(diǎn)位總氮濃度平均值對(duì)比

    表4 不同污染負(fù)荷削減方案下長(zhǎng)潭水庫(kù)3個(gè)代表點(diǎn)位總磷濃度平均值對(duì)比

    在現(xiàn)狀污染負(fù)荷不削減的情況下,長(zhǎng)潭水庫(kù)庫(kù)尾、庫(kù)中和庫(kù)首點(diǎn)位兩種水文典型年葉綠素a的峰值濃度均超過(guò)50μg/L,說(shuō)明目前長(zhǎng)潭水庫(kù)接納的污染物已經(jīng)超過(guò)環(huán)境容量,即使是在水動(dòng)力條件較好的豐水年,仍有藻類(lèi)暴發(fā)的可能;污染負(fù)荷削減10%,枯水年水庫(kù)葉綠素a的平均濃度僅有小幅度的降低,說(shuō)明目前長(zhǎng)潭水庫(kù)接納的污染物超標(biāo)較多,加之在枯水年,夏季的月平均入庫(kù)流量?jī)H為豐水年的35%,不良的水動(dòng)力條件下,污染削減10%對(duì)于遏制藻類(lèi)生長(zhǎng)幾乎沒(méi)有作用.污染負(fù)荷削減20%對(duì)于水質(zhì)的改善作用仍然有限,庫(kù)首點(diǎn)位葉綠素a的平均濃度僅降低5.12%.相應(yīng)地,在豐水年污染負(fù)荷削減10%后,水庫(kù)3個(gè)代表點(diǎn)位的葉綠素a濃度均可以降低10%以上,污染削減力度增大至20%后,葉綠素a平均濃度普遍可以下降20%左右(表5).

    表5 不同污染負(fù)荷削減方案下長(zhǎng)潭水庫(kù)3個(gè)代表點(diǎn)位葉綠素a濃度平均值對(duì)比

    長(zhǎng)潭水庫(kù)是典型的河道型水庫(kù),庫(kù)尾區(qū)域水深較淺,水流較快,表現(xiàn)出河道的特性,較大的流量產(chǎn)生的沖刷和稀釋作用使藻類(lèi)不能集中暴發(fā)[14],因此污染負(fù)荷削減后,庫(kù)尾區(qū)域的水質(zhì)改善程度較好.但在豐水年污染負(fù)荷不削減的情況下,庫(kù)尾區(qū)域葉綠素a的平均濃度最高,其原因是豐水年降雨量增多,庫(kù)尾地區(qū)上游的養(yǎng)殖業(yè)廢水隨著地表徑流直接進(jìn)入水庫(kù),為藻類(lèi)生長(zhǎng)提供大量營(yíng)養(yǎng)鹽,也說(shuō)明在營(yíng)養(yǎng)鹽過(guò)量的情況下水動(dòng)力條件對(duì)藻類(lèi)生長(zhǎng)的作用不明顯.相應(yīng)地,庫(kù)首區(qū)域水深變深,流速變慢,表現(xiàn)出湖泊的特性,污染負(fù)荷削減20%后,在枯水年葉綠素a平均濃度僅降低5.12%,說(shuō)明緩慢的流態(tài)是藻類(lèi)大量生長(zhǎng)的主要環(huán)境條件之一.

    污染負(fù)荷削減20%后,在豐水年的水動(dòng)力條件下,水庫(kù)3個(gè)代表點(diǎn)位的葉綠素a濃度峰值表現(xiàn)出明顯下降,但是在流量大幅度降低的枯水年,除庫(kù)尾點(diǎn)位,庫(kù)中和庫(kù)首點(diǎn)位葉綠素a峰值濃度沒(méi)有明顯變化.其原因是水庫(kù)庫(kù)尾區(qū)域表現(xiàn)出河道的特性,流速相對(duì)較快,污染物質(zhì)不易停留,大多積累在庫(kù)首區(qū)域,造成了庫(kù)首點(diǎn)位的葉綠素a峰值濃度相對(duì)庫(kù)尾和庫(kù)中點(diǎn)位始終偏高(表6).

    表6 不同污染負(fù)荷削減方案下長(zhǎng)潭水庫(kù)3個(gè)代表點(diǎn)位葉綠素a濃度峰值對(duì)比

    3 結(jié)論

    本文利用EFDC模型建立了河道型水庫(kù)——長(zhǎng)潭水庫(kù)的水動(dòng)力水生態(tài)數(shù)學(xué)模型,利用水位、水質(zhì)實(shí)測(cè)資料對(duì)模型進(jìn)行校核,實(shí)現(xiàn)了模型的本地化,模擬出了長(zhǎng)潭水庫(kù)藻類(lèi)生長(zhǎng)對(duì)于流域污染削減方案的響應(yīng),對(duì)于長(zhǎng)潭水庫(kù)防治富營(yíng)養(yǎng)化提供了指導(dǎo)性的對(duì)策,對(duì)于河道型水庫(kù)的管理與水質(zhì)改善也具有參考意義.

    1) 長(zhǎng)潭水庫(kù)是典型的河道型水庫(kù),庫(kù)尾區(qū)域具有河道的特性,水體交換程度好,污染物不易停留.庫(kù)首區(qū)域表現(xiàn)出湖泊的特性,水體流動(dòng)緩慢,污染物質(zhì)容易滯留,在適宜的環(huán)境條件下容易暴發(fā)藍(lán)藻水華.

    2) 現(xiàn)狀污染負(fù)荷削減10%可以對(duì)藻類(lèi)生長(zhǎng)產(chǎn)生一定程度的抑制作用,但在水動(dòng)力條件較差的枯水年,庫(kù)首點(diǎn)位葉綠素a的峰值濃度僅降低1.06%;增大污染負(fù)荷削減的力度并且通過(guò)水庫(kù)調(diào)度增加入庫(kù)流量,可以明顯改善庫(kù)尾區(qū)域水質(zhì)情況,葉綠素a的平均濃度可以降低26.78%,峰值濃度可以降低16.17%.

    3) 長(zhǎng)潭水庫(kù)上游養(yǎng)殖業(yè)產(chǎn)生的廢水對(duì)污染負(fù)荷貢獻(xiàn)最大,營(yíng)養(yǎng)鹽超標(biāo)的情況下即使增加長(zhǎng)興電站的下泄流量也難以改善庫(kù)首地區(qū)的富營(yíng)養(yǎng)化狀況,因此大力治理養(yǎng)殖業(yè)排污才是保障長(zhǎng)潭水庫(kù)水質(zhì)安全的長(zhǎng)遠(yuǎn)措施.

    4) EFDC模型對(duì)水質(zhì)基礎(chǔ)資料的完整性要求較高,由于沒(méi)有環(huán)保監(jiān)測(cè)斷面,長(zhǎng)潭水庫(kù)缺少長(zhǎng)系列的水質(zhì)實(shí)測(cè)資料,若要更加精確地預(yù)測(cè)水庫(kù)富營(yíng)養(yǎng)化,還需開(kāi)展更深入的資料收集工作.

    [1] 張 遠(yuǎn),鄭丙輝,富 國(guó)等.河道型水庫(kù)基于敏感性分區(qū)的營(yíng)養(yǎng)狀態(tài)標(biāo)準(zhǔn)與評(píng)價(jià)方法研究.環(huán)境科學(xué)學(xué)報(bào),2006,26(6):1016-1021.

    [2] 郝文彬,唐春燕,滑 磊等.引江濟(jì)太調(diào)水工程對(duì)太湖水動(dòng)力的調(diào)控效果.河海大學(xué)學(xué)報(bào):自然科學(xué)版,2012,40(2):129-133.

    [3] Li YP, Achary K, Yu ZB. Modeling impacts of Yangtze River water transfer on water ages in Lake Taihu, China.EcologicalEngineering, 2011, 37: 325-334.

    [4] 衛(wèi)志宏,唐雄飛,楊振祥等.洱海主要污染物允許排放量的控制分配.湖泊科學(xué),2013,25(5):665-673.DOI 10.18307/2013.0507.

    [5] 李 杰,林 晶,吳增茂.FallsLake水庫(kù)內(nèi)溶解物輸運(yùn)模擬.水科學(xué)進(jìn)展,2011,22(3):413-420.

    [6] 陳異暉.基于EFDC模型的滇池水質(zhì)模擬.云南環(huán)境科學(xué),2005,24(4):28-30.

    [7] 李一平,滑 磊,王沛芳等.河道型水庫(kù)水動(dòng)力特征與氣候條件的響應(yīng)關(guān)系.湖泊科學(xué),2013,25(3):317-323.DOI 10.18307/2013.0301.

    [8] 王 超,高超越,王沛芳等.廣東長(zhǎng)潭水庫(kù)富營(yíng)養(yǎng)化與浮游植物分布特征.湖泊科學(xué),2013,25(5):749-755.DOI 10.18307/2013.0518.

    [9] 李錦秀,杜 斌,孫以三.水動(dòng)力條件對(duì)富營(yíng)養(yǎng)化影響規(guī)律探討.水利水電技術(shù),2005,36(5):15-18.

    [10] Lei Z, Yu ZL, Rui Z. A Three-dimensional water quality modeling approach for exploring the eutrophication responses to load reduction scenarios in Lake Yilong(China).EnvironmentalPollution, 2013, 177: 13-21.

    [11] Guo ZW, Zong XX. Prediction of algal blooming using EFDC model: case study in the Daoxiang Lake.EcologicalModelling, 2011, 222:1245-1252.

    [12] 華錦彪.洋河水庫(kù)“水華”發(fā)生的實(shí)驗(yàn)研究[學(xué)位論文].北京:北京大學(xué),1991.

    [13] 王煥校.污染生態(tài)學(xué).北京:高等教育出版社,2000:147.

    [14] 梁培瑜,王 烜,馬芳冰等.水動(dòng)力條件對(duì)水體富營(yíng)養(yǎng)化的影響.湖泊科學(xué),2013,25(4):455-462.DOI 10.18307/2013.0401.

    Response of algae growth to pollution reduction of drainage basin based on EFDC model for channel reservoirs: A case of Changtan Reservoir, Guangdong Province

    LI Yiping1, WANG Jingyu1& HUA Lei2

    (1:CollegeofEnvironment,HohaiUniversity,Nanjing210098,P.R.China)

    (2:ZhejiangInstituteofHydraulics&Estuary,Hangzhou310020,P.R.China)

    The EFDC model has been used to study response of water quality to varying degrees of pollution reduction in Changtan Reservoir. After verification of the model, a water quality model of Changtan Reservoir was built. Algae growth was simulated under the present situation of pollution load. Pollution reduction scenarios based on two hydrologic years were simulated using the water quality model. The results indicated that overmuch nutrient could lead to eutrophication even in wet years. Under 10% reduction of pollution load, average chlorophyll-a concentrations in three areas (reservoir tail, middle and head areas) decreased by 13.99%, 12.00% and 10.35% in wet year and by 8.42%, 5.63% and 2.10% in dry year, respectively. Under the condition of pollution reduced by 20%, chlorophyll-a concentrations in three areas decreased by 26.78%, 19.25% and 17.04% in wet year and by 11.72%, 7.97% and 5.12% in dry year, respectively. The water in the reservoir tail area behaved flowing-water similar to rivers, while the water in the reservoir head area also behaved still-water similar to lakes, which has been more dramatic in dry year and the average chlorophyll-a concentration decreased by 5.12%.

    Channel reservoir; Changtan Reservoir; EFDC model; pollution load reduction; response of algae growth

    *水利部公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目(201101020)資助.2014-07-28收稿;2014-12-29收修改稿.李一平(1978~),男,博士,副教授;E-mail:liyiping@hhu.edu.cn.

    J.LakeSci.(湖泊科學(xué)), 2015, 27(5): 811-818

    DOI 10.18307/2015.0507

    ?2015 byJournalofLakeSciences

    **通信作者;E-mail:wjyhhu@126.com.

    猜你喜歡
    點(diǎn)位藻類(lèi)葉綠素
    1年期國(guó)債收益率跌至1%
    銀行家(2025年1期)2025-02-08 00:00:00
    藻類(lèi)水華控制技術(shù)及應(yīng)用
    機(jī)器人快速示教方法及示教點(diǎn)位姿變換的研究
    細(xì)菌和藻類(lèi)先移民火星
    軍事文摘(2020年20期)2020-11-16 00:31:40
    提取葉綠素
    桃樹(shù)葉綠素含量與SPAD值呈極顯著正相關(guān)
    吃蔬菜有個(gè)“321模式” 三兩葉菜類(lèi),二兩其他類(lèi),一兩菌藻類(lèi)
    葉綠素家族概述
    浮游藻類(lèi)在水體PAHs富集與降解中的研究
    垂直面內(nèi)建立基線的特殊點(diǎn)位高程測(cè)量法
    亚洲成人久久爱视频| 熟女人妻精品中文字幕| 免费观看的影片在线观看| 国产精华一区二区三区| 人妻少妇偷人精品九色| 在线看三级毛片| 蜜桃亚洲精品一区二区三区| 精品久久久久久成人av| 日本三级黄在线观看| 免费av观看视频| 我的老师免费观看完整版| 免费av毛片视频| 日韩欧美免费精品| 干丝袜人妻中文字幕| 久久久精品大字幕| 一区二区三区激情视频| 亚洲18禁久久av| 18禁裸乳无遮挡免费网站照片| 九九爱精品视频在线观看| 国产精品99久久久久久久久| 18禁黄网站禁片免费观看直播| 国产精品自产拍在线观看55亚洲| 在现免费观看毛片| 99精品在免费线老司机午夜| 日韩欧美 国产精品| 99久久精品热视频| 一区二区三区四区激情视频 | 亚洲精品久久国产高清桃花| 老女人水多毛片| 婷婷精品国产亚洲av| 日韩中字成人| 国产精品久久久久久av不卡| 精品久久久久久久久久久久久| 伦理电影大哥的女人| 搡老岳熟女国产| 亚洲性夜色夜夜综合| 午夜爱爱视频在线播放| 成人一区二区视频在线观看| 久久久久久九九精品二区国产| 国产精品一区二区三区四区久久| 无人区码免费观看不卡| 国产成人一区二区在线| 久久精品久久久久久噜噜老黄 | 亚洲精品一卡2卡三卡4卡5卡| 97人妻精品一区二区三区麻豆| 俄罗斯特黄特色一大片| 91麻豆av在线| 九九在线视频观看精品| 国产男人的电影天堂91| 欧美xxxx性猛交bbbb| 99精品在免费线老司机午夜| 国产三级在线视频| 嫩草影院精品99| 麻豆精品久久久久久蜜桃| 欧美最黄视频在线播放免费| 内射极品少妇av片p| 欧美不卡视频在线免费观看| 亚洲成人久久性| 久久欧美精品欧美久久欧美| 久久久久久九九精品二区国产| 色综合亚洲欧美另类图片| 狂野欧美激情性xxxx在线观看| 91麻豆精品激情在线观看国产| 亚洲精品成人久久久久久| 国产精品久久久久久精品电影| 哪里可以看免费的av片| 校园人妻丝袜中文字幕| 久久久色成人| 九九在线视频观看精品| 日韩欧美 国产精品| eeuss影院久久| 亚洲18禁久久av| 精品免费久久久久久久清纯| 真人一进一出gif抽搐免费| 97超级碰碰碰精品色视频在线观看| 久久久国产成人精品二区| 久久国产精品人妻蜜桃| 免费在线观看影片大全网站| 嫩草影院精品99| 桃色一区二区三区在线观看| 黄色一级大片看看| 一a级毛片在线观看| 欧美在线一区亚洲| 久久久精品大字幕| 午夜影院日韩av| 天美传媒精品一区二区| 亚洲精品456在线播放app | 嫩草影院入口| 人妻丰满熟妇av一区二区三区| 色哟哟哟哟哟哟| 女生性感内裤真人,穿戴方法视频| 国产人妻一区二区三区在| 国产av麻豆久久久久久久| 欧美+亚洲+日韩+国产| a级毛片a级免费在线| 一级黄片播放器| netflix在线观看网站| 最好的美女福利视频网| 嫩草影院入口| 天堂网av新在线| 免费在线观看影片大全网站| 天美传媒精品一区二区| 日本一本二区三区精品| 欧美最新免费一区二区三区| 欧美人与善性xxx| 欧美性猛交╳xxx乱大交人| 亚洲无线在线观看| 亚洲欧美日韩卡通动漫| 国产精品久久久久久av不卡| 1024手机看黄色片| 国产精品久久久久久久久免| 久久精品国产亚洲av天美| 波多野结衣巨乳人妻| 国产亚洲精品av在线| 国产综合懂色| 欧美国产日韩亚洲一区| 看黄色毛片网站| 全区人妻精品视频| 精品久久久久久久久亚洲 | 老司机午夜福利在线观看视频| 成人综合一区亚洲| 国产乱人伦免费视频| 日韩精品青青久久久久久| 久久精品国产清高在天天线| 日韩欧美免费精品| 国产精品一区www在线观看 | 色综合亚洲欧美另类图片| 日本爱情动作片www.在线观看 | 99热这里只有是精品50| 国产三级在线视频| 男插女下体视频免费在线播放| av专区在线播放| 午夜影院日韩av| 一个人观看的视频www高清免费观看| 老女人水多毛片| ponron亚洲| 国产精品精品国产色婷婷| 蜜桃亚洲精品一区二区三区| 亚洲中文字幕日韩| 午夜亚洲福利在线播放| 中文字幕人妻熟人妻熟丝袜美| 亚洲中文字幕一区二区三区有码在线看| 床上黄色一级片| 美女高潮的动态| 99久久九九国产精品国产免费| 香蕉av资源在线| 国产欧美日韩一区二区精品| 伊人久久精品亚洲午夜| 国产高清视频在线观看网站| 91精品国产九色| 亚洲专区国产一区二区| 久久精品影院6| 日本黄色片子视频| 亚洲电影在线观看av| 国产精品三级大全| 九九爱精品视频在线观看| 亚洲中文字幕一区二区三区有码在线看| 综合色av麻豆| 日韩欧美国产在线观看| 精品人妻视频免费看| 亚洲精华国产精华精| 精品国内亚洲2022精品成人| 色噜噜av男人的天堂激情| 亚洲va在线va天堂va国产| 免费一级毛片在线播放高清视频| 男人舔女人下体高潮全视频| 一本一本综合久久| 亚洲精华国产精华液的使用体验 | 免费看a级黄色片| 丝袜美腿在线中文| 日韩在线高清观看一区二区三区 | 国产视频内射| 白带黄色成豆腐渣| 乱人视频在线观看| 九色国产91popny在线| 精品人妻偷拍中文字幕| 在线观看午夜福利视频| 小蜜桃在线观看免费完整版高清| 黄色女人牲交| 亚洲成a人片在线一区二区| 国内精品一区二区在线观看| 国产精品国产三级国产av玫瑰| 一级黄片播放器| 黄色日韩在线| 国产精品精品国产色婷婷| 久久久久久久久久久丰满 | 美女大奶头视频| 欧美高清成人免费视频www| 日韩中文字幕欧美一区二区| 在线观看一区二区三区| 最近视频中文字幕2019在线8| 成年人黄色毛片网站| 日本一二三区视频观看| 亚洲avbb在线观看| 欧美三级亚洲精品| 国产真实伦视频高清在线观看 | 中出人妻视频一区二区| netflix在线观看网站| 91午夜精品亚洲一区二区三区 | 亚洲欧美日韩高清在线视频| 亚洲人成网站在线播放欧美日韩| 午夜福利欧美成人| 99国产精品一区二区蜜桃av| 老司机午夜福利在线观看视频| 欧美人与善性xxx| 午夜福利在线观看免费完整高清在 | 中文字幕av在线有码专区| 少妇被粗大猛烈的视频| 亚洲va在线va天堂va国产| 亚洲中文日韩欧美视频| 国内揄拍国产精品人妻在线| 最新在线观看一区二区三区| 国产免费一级a男人的天堂| 欧美日韩乱码在线| 免费一级毛片在线播放高清视频| 国产主播在线观看一区二区| xxxwww97欧美| 日本-黄色视频高清免费观看| 极品教师在线免费播放| 国产精品野战在线观看| 国产精品国产三级国产av玫瑰| 色视频www国产| 免费观看在线日韩| 久久天躁狠狠躁夜夜2o2o| 亚洲午夜理论影院| 亚洲一级一片aⅴ在线观看| 啦啦啦啦在线视频资源| 国产高清视频在线播放一区| 午夜福利在线观看吧| 久久久久免费精品人妻一区二区| 亚洲av免费在线观看| 久久久国产成人精品二区| 国产熟女欧美一区二区| 制服丝袜大香蕉在线| 亚洲av成人精品一区久久| 国产精品福利在线免费观看| 亚洲人与动物交配视频| 日本黄大片高清| 亚洲av第一区精品v没综合| 国产aⅴ精品一区二区三区波| 国产精品1区2区在线观看.| 国产不卡一卡二| videossex国产| 亚洲av成人精品一区久久| 婷婷色综合大香蕉| 日本黄色片子视频| 国产三级中文精品| av女优亚洲男人天堂| 一个人看的www免费观看视频| 亚洲无线在线观看| 窝窝影院91人妻| 国产主播在线观看一区二区| 亚洲精品456在线播放app | 国产麻豆成人av免费视频| 网址你懂的国产日韩在线| 九九热线精品视视频播放| 国产大屁股一区二区在线视频| 国产高潮美女av| 国产精品女同一区二区软件 | a级毛片a级免费在线| 乱码一卡2卡4卡精品| 精品久久久久久久久亚洲 | eeuss影院久久| 非洲黑人性xxxx精品又粗又长| 欧美日韩综合久久久久久 | 精品国内亚洲2022精品成人| 久久久久久久久久黄片| 内地一区二区视频在线| 久久国产乱子免费精品| 免费在线观看成人毛片| 老熟妇仑乱视频hdxx| 久久久久九九精品影院| .国产精品久久| 日本免费a在线| 美女被艹到高潮喷水动态| 久久婷婷人人爽人人干人人爱| 在线看三级毛片| 成人av一区二区三区在线看| 久久久色成人| 女同久久另类99精品国产91| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av涩爱 | 99热网站在线观看| 直男gayav资源| 国产真实伦视频高清在线观看 | 男女下面进入的视频免费午夜| 国产日本99.免费观看| 亚洲中文字幕一区二区三区有码在线看| 免费人成视频x8x8入口观看| 免费看光身美女| 淫秽高清视频在线观看| 免费人成在线观看视频色| 久久中文看片网| 久久精品国产99精品国产亚洲性色| 亚洲专区中文字幕在线| 嫁个100分男人电影在线观看| 亚洲最大成人中文| 91在线精品国自产拍蜜月| 亚洲午夜理论影院| 国产一区二区三区在线臀色熟女| 在线国产一区二区在线| 露出奶头的视频| 淫妇啪啪啪对白视频| 搡女人真爽免费视频火全软件 | 日本一本二区三区精品| 97超视频在线观看视频| 级片在线观看| 国产91精品成人一区二区三区| 欧美极品一区二区三区四区| 日本色播在线视频| 男人狂女人下面高潮的视频| 国产精品无大码| 久久精品国产亚洲网站| 此物有八面人人有两片| 成人av在线播放网站| 亚洲电影在线观看av| 97热精品久久久久久| 老熟妇乱子伦视频在线观看| 日日干狠狠操夜夜爽| av在线亚洲专区| 淫秽高清视频在线观看| netflix在线观看网站| 亚洲内射少妇av| 男女那种视频在线观看| 丰满的人妻完整版| 黄色丝袜av网址大全| 无遮挡黄片免费观看| 在线看三级毛片| 亚洲精品亚洲一区二区| 国内精品久久久久久久电影| 婷婷精品国产亚洲av在线| 99精品久久久久人妻精品| 熟女电影av网| 欧美成人一区二区免费高清观看| av天堂中文字幕网| 日本 欧美在线| 黄色配什么色好看| 99国产精品一区二区蜜桃av| 免费无遮挡裸体视频| 天美传媒精品一区二区| 国产免费男女视频| 天美传媒精品一区二区| 18+在线观看网站| 精品国产三级普通话版| 桃色一区二区三区在线观看| 国产av在哪里看| 97热精品久久久久久| 欧美日本视频| 五月玫瑰六月丁香| 免费看日本二区| 国语自产精品视频在线第100页| 久久中文看片网| 久久99热6这里只有精品| 伦精品一区二区三区| 波野结衣二区三区在线| 69人妻影院| 麻豆一二三区av精品| 久久久久久久久久黄片| 九九爱精品视频在线观看| 久久人人精品亚洲av| 日本欧美国产在线视频| 日韩,欧美,国产一区二区三区 | 久久精品夜夜夜夜夜久久蜜豆| 最近最新免费中文字幕在线| 啦啦啦韩国在线观看视频| 国产真实伦视频高清在线观看 | 欧美日韩综合久久久久久 | 久久久国产成人精品二区| 国产精品久久久久久久电影| 18禁在线播放成人免费| 老师上课跳d突然被开到最大视频| 久久久午夜欧美精品| 天天躁日日操中文字幕| 国产欧美日韩精品一区二区| 亚洲精品成人久久久久久| 久久午夜亚洲精品久久| 伊人久久精品亚洲午夜| 男女边吃奶边做爰视频| 国产av一区在线观看免费| 一进一出抽搐gif免费好疼| 久久久久国内视频| 日本黄色视频三级网站网址| 亚洲国产色片| 婷婷精品国产亚洲av在线| 97超级碰碰碰精品色视频在线观看| 最好的美女福利视频网| 特大巨黑吊av在线直播| 99久久精品热视频| 十八禁国产超污无遮挡网站| 国产91精品成人一区二区三区| 亚洲天堂国产精品一区在线| 亚洲自偷自拍三级| 麻豆一二三区av精品| 国产一区二区三区在线臀色熟女| 欧美人与善性xxx| 日日干狠狠操夜夜爽| 99久久九九国产精品国产免费| 欧美zozozo另类| 看十八女毛片水多多多| 男插女下体视频免费在线播放| 国产极品精品免费视频能看的| 久久久久国内视频| 国产女主播在线喷水免费视频网站 | 在线免费观看不下载黄p国产 | 亚洲不卡免费看| 男女做爰动态图高潮gif福利片| 欧洲精品卡2卡3卡4卡5卡区| 99久久中文字幕三级久久日本| 精品久久久久久久久久免费视频| 国产精品三级大全| 国产精品不卡视频一区二区| 我要搜黄色片| 一个人看的www免费观看视频| 日韩欧美精品免费久久| 99久久久亚洲精品蜜臀av| 无人区码免费观看不卡| 久久国产乱子免费精品| 国产av一区在线观看免费| 丰满人妻一区二区三区视频av| 成人美女网站在线观看视频| 久久久久久大精品| 少妇熟女aⅴ在线视频| 日韩av在线大香蕉| 在线观看免费视频日本深夜| 好男人在线观看高清免费视频| 大型黄色视频在线免费观看| 日本黄色片子视频| 国产高清激情床上av| 老师上课跳d突然被开到最大视频| 波多野结衣高清作品| 赤兔流量卡办理| 99九九线精品视频在线观看视频| 少妇熟女aⅴ在线视频| 毛片一级片免费看久久久久 | 我要搜黄色片| 免费av不卡在线播放| 美女 人体艺术 gogo| 男女啪啪激烈高潮av片| 欧美丝袜亚洲另类 | 中文在线观看免费www的网站| 国产成年人精品一区二区| 日韩大尺度精品在线看网址| av在线天堂中文字幕| 国产精品国产高清国产av| 99riav亚洲国产免费| 永久网站在线| 亚洲成a人片在线一区二区| 婷婷色综合大香蕉| 精品人妻一区二区三区麻豆 | 国产精品女同一区二区软件 | 99精品久久久久人妻精品| 一进一出抽搐gif免费好疼| xxxwww97欧美| 婷婷六月久久综合丁香| 欧美一级a爱片免费观看看| 狂野欧美激情性xxxx在线观看| 日本熟妇午夜| 国产精华一区二区三区| 午夜福利在线在线| 亚洲美女搞黄在线观看 | 国产真实乱freesex| av中文乱码字幕在线| 欧美另类亚洲清纯唯美| 精品乱码久久久久久99久播| 亚洲av中文字字幕乱码综合| 少妇裸体淫交视频免费看高清| 免费av不卡在线播放| 最近在线观看免费完整版| 此物有八面人人有两片| 国产精品久久电影中文字幕| 长腿黑丝高跟| 老司机深夜福利视频在线观看| 精品久久久久久久末码| 国内少妇人妻偷人精品xxx网站| 又爽又黄无遮挡网站| 国产精品久久久久久精品电影| 成人一区二区视频在线观看| 精品免费久久久久久久清纯| 一区二区三区激情视频| 国产日本99.免费观看| 午夜a级毛片| 国产色婷婷99| 校园人妻丝袜中文字幕| 亚洲电影在线观看av| x7x7x7水蜜桃| 男女做爰动态图高潮gif福利片| 啦啦啦观看免费观看视频高清| 欧美性感艳星| 亚洲人成网站在线播放欧美日韩| 非洲黑人性xxxx精品又粗又长| 欧美bdsm另类| 淫妇啪啪啪对白视频| 听说在线观看完整版免费高清| 国产精品国产高清国产av| 一边摸一边抽搐一进一小说| 免费无遮挡裸体视频| 欧美xxxx黑人xx丫x性爽| 成年女人永久免费观看视频| 欧美精品啪啪一区二区三区| 我要看日韩黄色一级片| 真人做人爱边吃奶动态| 国产乱人视频| 88av欧美| 欧美极品一区二区三区四区| 精品久久久久久,| 精品人妻熟女av久视频| 在线观看免费视频日本深夜| 亚洲va在线va天堂va国产| 久久久久国产精品人妻aⅴ院| 有码 亚洲区| av在线老鸭窝| 国产私拍福利视频在线观看| 国产综合懂色| 国产亚洲91精品色在线| 天堂av国产一区二区熟女人妻| 黄色视频,在线免费观看| 色5月婷婷丁香| 嫩草影院精品99| 免费av观看视频| 九九热线精品视视频播放| 赤兔流量卡办理| 观看美女的网站| 黄色配什么色好看| 搡女人真爽免费视频火全软件 | 国产亚洲精品久久久久久毛片| 国产精品嫩草影院av在线观看 | 国产成年人精品一区二区| 亚洲三级黄色毛片| 九九热线精品视视频播放| 免费一级毛片在线播放高清视频| 中文字幕熟女人妻在线| 欧美日韩国产亚洲二区| 亚洲av一区综合| 色综合色国产| 久久午夜福利片| 欧美三级亚洲精品| 国产麻豆成人av免费视频| 人人妻,人人澡人人爽秒播| 无人区码免费观看不卡| 99国产极品粉嫩在线观看| 久久99热6这里只有精品| 国产一区二区在线观看日韩| 乱码一卡2卡4卡精品| 琪琪午夜伦伦电影理论片6080| 日韩在线高清观看一区二区三区 | 欧美日韩亚洲国产一区二区在线观看| 亚洲美女黄片视频| 国产精品久久电影中文字幕| 日本精品一区二区三区蜜桃| 日韩在线高清观看一区二区三区 | 午夜福利成人在线免费观看| 最新中文字幕久久久久| 尤物成人国产欧美一区二区三区| 88av欧美| 免费一级毛片在线播放高清视频| 嫩草影院精品99| 国产精品乱码一区二三区的特点| 一个人观看的视频www高清免费观看| 小说图片视频综合网站| 久久午夜亚洲精品久久| 欧美又色又爽又黄视频| 精品午夜福利视频在线观看一区| 男女做爰动态图高潮gif福利片| 老女人水多毛片| 有码 亚洲区| 在线观看舔阴道视频| 男女之事视频高清在线观看| 少妇熟女aⅴ在线视频| 夜夜爽天天搞| av.在线天堂| 三级毛片av免费| 色噜噜av男人的天堂激情| 国产在视频线在精品| 97超视频在线观看视频| 成人一区二区视频在线观看| 一个人免费在线观看电影| 一级黄色大片毛片| 可以在线观看的亚洲视频| www日本黄色视频网| 色尼玛亚洲综合影院| 亚洲精品色激情综合| 国产伦精品一区二区三区四那| 天堂动漫精品| 尾随美女入室| 午夜福利在线观看吧| 国产精品亚洲一级av第二区| 亚洲av免费高清在线观看| 亚洲专区中文字幕在线| 亚洲欧美精品综合久久99| 深夜精品福利| 少妇的逼水好多| 国内毛片毛片毛片毛片毛片| 国产精品一区二区三区四区免费观看 | 成年女人永久免费观看视频| 久久久久免费精品人妻一区二区| 91久久精品国产一区二区成人| 国产免费一级a男人的天堂| 干丝袜人妻中文字幕| 97超视频在线观看视频| 韩国av一区二区三区四区| 日本 欧美在线| av.在线天堂| 午夜激情福利司机影院| 亚洲国产精品sss在线观看| 亚洲男人的天堂狠狠| 亚洲18禁久久av| 国产精品不卡视频一区二区| 日韩中字成人| 狠狠狠狠99中文字幕| 久久久久国产精品人妻aⅴ院| 精品久久久久久久末码| 1000部很黄的大片| 很黄的视频免费| 精品国产三级普通话版| 国产色爽女视频免费观看|