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

    耦合水文情勢(shì)及魚(yú)類(lèi)需求的生態(tài)調(diào)度研究

    2018-07-05 08:54:00高志強(qiáng)丁偉唐榕李昱張弛
    南水北調(diào)與水利科技 2018年2期

    高志強(qiáng) 丁偉 唐榕 李昱 張弛

    摘要:針對(duì)豐滿(mǎn)水庫(kù)下游魚(yú)類(lèi)多樣性減少的問(wèn)題,在滿(mǎn)足最小生態(tài)流量需求的基礎(chǔ)上,開(kāi)展了營(yíng)造適宜魚(yú)類(lèi)的水文情勢(shì)研究。通過(guò)篩選改變度大且適宜魚(yú)類(lèi)的水文情勢(shì)指標(biāo)構(gòu)建生態(tài)完整度函數(shù),將此作為生態(tài)目標(biāo)建立水庫(kù)多目標(biāo)優(yōu)化調(diào)度模型。調(diào)度結(jié)果表明,所建立的生態(tài)調(diào)度模型在滿(mǎn)足發(fā)電、供水和最小生態(tài)需水量的基礎(chǔ)上,可將生態(tài)完整度至少提高48%,其中,漲水率、高流量脈沖數(shù)指標(biāo)隸屬度分別增加了0.112、0.261,且漲水率指標(biāo)在枯水年改善程度好,高流量脈沖數(shù)指標(biāo)在枯水年改善程度大、豐水年改善程度小。本文所建立的模型顯著提高了適宜魚(yú)類(lèi)水文情勢(shì)的生態(tài)完整度,可為有魚(yú)類(lèi)修復(fù)任務(wù)的流域提供參考。

    關(guān)鍵詞:豐滿(mǎn)水庫(kù);水文情勢(shì);魚(yú)類(lèi);生態(tài)完整度;生態(tài)調(diào)度

    中圖分類(lèi)號(hào):TV697.1文獻(xiàn)標(biāo)志碼:A文章編號(hào):

    16721683(2018)02001407

    Abstract:

    To solve the problem of the dwindling fish resources in the lower reaches of the Fengman reservoir,we conducted the study on hydrological regime suitable for fishes on the basis of satisfying the minimum ecological water demand.We selected hydrological regime indicators that changed severely and were suitable for fishes,and used them to establish an ecological integrity function.Using the proposed function as the ecological objective,we established a multiobjective optimal scheduling model.Results showed that this model could improve ecological integrity by at least 48%.The memberships of flow rising rate and high flow pulse number increased by 0112 and 0261 respectively.The flow rising rate was well improved in dry years.The high flow pulse number was improved greatly in dry years and slightly in wet years.This model can significantly improve the ecological integrity,and can provide a reference for the watershed with fish restoration tasks.

    Key words:

    Fengman reservoir;hydrological regime;fish;ecological integrity;reservoir scheduling considering ecology

    1研究背景

    自然條件下,河流水文情勢(shì)的豐枯變化維持著水生態(tài)系統(tǒng)的穩(wěn)定[1]。然而水庫(kù)的“蓄豐補(bǔ)枯”調(diào)節(jié)坦化了天然徑流,極大地改變了天然水文情勢(shì)的變化,產(chǎn)生了河流連通能力下降、下游濕地萎縮、生物棲息地喪失等一系列的生態(tài)問(wèn)題[23],因此在水庫(kù)調(diào)度中必須要考慮生態(tài)問(wèn)題。

    目前水庫(kù)生態(tài)調(diào)度主要是保證泄流滿(mǎn)足最小生態(tài)需水量[4],然而這種調(diào)度方式?jīng)]有反映河流天然[HJ1.98mm]狀態(tài)下的峰現(xiàn)時(shí)間、洪水脈沖、漲水次數(shù)等水文情勢(shì)信息,對(duì)生物的洄游、產(chǎn)卵、生長(zhǎng)等過(guò)程帶來(lái)不利影響[5],因此生態(tài)調(diào)度除了滿(mǎn)足水量的需求外,還需營(yíng)造適宜的水文情勢(shì)[6]。水文情勢(shì)可通過(guò)生態(tài)指標(biāo)量化,在調(diào)度時(shí)進(jìn)行調(diào)控滿(mǎn)足,由于指標(biāo)眾多[711],如何合理地篩選關(guān)鍵指標(biāo)進(jìn)行控制十分重要。當(dāng)生態(tài)監(jiān)測(cè)資料比較充分時(shí),可以采用數(shù)據(jù)挖掘法篩選指標(biāo),如YiChen E.Yang[12]利用GP挖掘法篩選了IHA指標(biāo),并建立了魚(yú)類(lèi)豐度和IHA指標(biāo)之間的函數(shù)關(guān)系;當(dāng)生態(tài)監(jiān)測(cè)資料缺乏時(shí),一般按個(gè)人經(jīng)驗(yàn)[1316]和生態(tài)指標(biāo)改變度大小的順序[1718]進(jìn)行篩選,如Suen[13]根據(jù)個(gè)人經(jīng)驗(yàn)挑選代表石門(mén)水庫(kù)下游生態(tài)系統(tǒng)的指標(biāo),研究了生態(tài)指標(biāo)改善對(duì)興利目標(biāo)的影響;楊娜[17]通過(guò)篩選改變度最大的前6個(gè)水文情勢(shì)指標(biāo),研究了丹江口水庫(kù)在水量分配和天然模式兩種調(diào)度情境下各指標(biāo)的改善程度。但這兩種方法都沒(méi)有結(jié)合流域的生態(tài)保護(hù)對(duì)象篩選指標(biāo),可能導(dǎo)致所修復(fù)水文情勢(shì)指標(biāo)對(duì)保護(hù)對(duì)象的作用并不明顯,達(dá)不到修復(fù)目的。

    本文針對(duì)第二松花江流域豐滿(mǎn)水庫(kù)下游河段產(chǎn)漂浮性卵洄游魚(yú)類(lèi)資源嚴(yán)重衰退問(wèn)題,篩選水文改變度大且適合魚(yú)類(lèi)需求的生態(tài)指標(biāo)構(gòu)建生態(tài)完整度函數(shù),在保證發(fā)電和供水任務(wù)的基礎(chǔ)上,建立考慮生態(tài)完整度的水庫(kù)多目標(biāo)優(yōu)化調(diào)度模型,探討調(diào)度結(jié)果對(duì)發(fā)電和供水的影響,分析各生態(tài)指標(biāo)多年平均隸屬度的變化情況以及年際變化情況。

    2.1豐滿(mǎn)水庫(kù)特征及流量數(shù)據(jù)

    豐滿(mǎn)水庫(kù)位于我國(guó)第二松花江流域,壩址在吉林省吉林市城區(qū)上游24 km處,是一座以發(fā)電為主,兼有防洪、供水、養(yǎng)殖的綜合水利樞紐工程,其電站裝機(jī)為1 002.5 MW,保證出力為166 MW;水庫(kù)下游的主要供水任務(wù)包括:長(zhǎng)春市引水11 m3/s,下游生態(tài)下泄150 m3/s,吉林市區(qū)以下灌溉期(5月-6月)取水350 m3/s,且生態(tài)下泄水量可兼作為下游的灌溉用水。

    豐滿(mǎn)水庫(kù)始建于1937年,沒(méi)有建庫(kù)前河道的天然水文數(shù)據(jù)序列,只在1946年以后有連續(xù)的出庫(kù)流量資料,為此本文在判別水文指標(biāo)改變度時(shí),以1946-1990年(1991年豐滿(mǎn)水庫(kù)上游修建了白山水庫(kù))反推的入庫(kù)流量數(shù)據(jù)序列作為下游河道天然水文數(shù)據(jù)序列;在水庫(kù)調(diào)節(jié)計(jì)算時(shí),采用1946-2009年豐滿(mǎn)水庫(kù)的入庫(kù)流量數(shù)據(jù)。

    2.2漂浮性洄游魚(yú)類(lèi)調(diào)查

    豐滿(mǎn)水庫(kù)壩址至飲馬河口區(qū)段是產(chǎn)漂浮性卵的[HJ]洄游魚(yú)類(lèi)如“四大家魚(yú)”的主要產(chǎn)卵區(qū),分布有鰲龍河口、榆樹(shù)、龍王廟、飲馬河口等產(chǎn)卵場(chǎng),見(jiàn)圖1。根據(jù)文獻(xiàn)[19],豐滿(mǎn)水庫(kù)下游魚(yú)類(lèi)多樣性不斷降低,主要是由于水庫(kù)修建及調(diào)度管理、環(huán)境污染、私捕濫撈、河道采砂等造成的,其中水庫(kù)的調(diào)度管理影響較大,對(duì)魚(yú)類(lèi)的洄游產(chǎn)卵、覓食生長(zhǎng)過(guò)程帶來(lái)不利影響,為此,豐滿(mǎn)水庫(kù)在滿(mǎn)足發(fā)電和供水任務(wù)的基礎(chǔ)上,還急需開(kāi)展?fàn)I造適宜魚(yú)類(lèi)水文情勢(shì)的研究。

    2.3生態(tài)目標(biāo)函數(shù)確定方法

    本文所提出的生態(tài)目標(biāo)確定方法,首先計(jì)算生態(tài)指標(biāo)改變度,然后選擇改變度大且影響?hù)~(yú)類(lèi)的指標(biāo)作為待修復(fù)對(duì)象,最后利用二元比較法確定生態(tài)目標(biāo)函數(shù)。

    (1)水文情勢(shì)指標(biāo)改變度計(jì)算。

    首先根據(jù)出流過(guò)程對(duì)水文情勢(shì)改變程度進(jìn)行評(píng)價(jià),采用Richer[7]所提出的水文變化指標(biāo)體系(Indicators of Hydrologic Alteration,IHA)來(lái)衡量水文情勢(shì)的變化。

    為了便于水庫(kù)實(shí)際調(diào)度的控制管理,這里定義漲水(落水)率指標(biāo)的單位為流量,即首先求每次漲水(落水)前后的流量差值,然后將年內(nèi)各次漲水(落水)前后流量差值的均值作為該年漲水率(落水)率指標(biāo)。

    在判斷水文情勢(shì)變化程度時(shí),主要評(píng)價(jià)調(diào)度后各水文情勢(shì)指標(biāo)相對(duì)于天然狀態(tài)下的改變程度,本文采用高斯隸屬度函數(shù)(式(1))反映水文情勢(shì)指標(biāo)相對(duì)于天然狀態(tài)的隸屬程度[15],采用改變度函數(shù)評(píng)價(jià)各指標(biāo)的改變程度(式(2))。鑒于豐滿(mǎn)水庫(kù)沒(méi)有建庫(kù)前的水文資料,本文以1946-1990年反推的入庫(kù)流量序列作為下游河道天然流量序列。

    μ(xi)=e-[SX(](xi-mi)2[]2σ2i[SX)][JY](1)

    v(xi)=1-μ(xi)[JY](2)

    式中:μ(xi)、v(xi)分別為第i個(gè)水文指標(biāo)的隸屬度函數(shù)和改變度函數(shù);xi為調(diào)度后第i個(gè)水文指標(biāo)年值;mi、σi分別為天然狀態(tài)下第i個(gè)水文指標(biāo)的均值和標(biāo)準(zhǔn)差。

    在天然水文情勢(shì)變化評(píng)估時(shí),IHA指標(biāo)改變度大于067的,認(rèn)為是高度改變;在033~067之間,認(rèn)為是中度改變;小于033的認(rèn)為是低度改變;其中高度和中度改變的指標(biāo)需要修復(fù)[20]。根據(jù)豐滿(mǎn)水庫(kù)各指標(biāo)的改變程度,如表1所示,需要修復(fù)的指標(biāo)有26個(gè),說(shuō)明豐滿(mǎn)水庫(kù)的調(diào)度極大地改變了下游的天然水文情勢(shì)。

    (2)魚(yú)類(lèi)需求指標(biāo)確定。

    由于豐滿(mǎn)水庫(kù)壩址至飲馬河口區(qū)段產(chǎn)漂浮性卵的洄游魚(yú)類(lèi)處于嚴(yán)重衰退狀態(tài),洄游產(chǎn)卵規(guī)模減少,本文以產(chǎn)漂浮性卵的洄游魚(yú)類(lèi)為保護(hù)對(duì)象,選擇水文情勢(shì)改變度大且對(duì)洄游產(chǎn)卵及魚(yú)卵發(fā)育有重要作用的指標(biāo)進(jìn)行修復(fù)。根據(jù)洄游產(chǎn)卵以及魚(yú)卵發(fā)育時(shí)所需要的洪水脈沖、漲落水率以及漲水持續(xù)時(shí)間等水文情勢(shì)刺激[2122],確定魚(yú)類(lèi)需求指標(biāo)如下:a.高流量脈沖發(fā)生次數(shù)以及高流量脈沖持續(xù)時(shí)間指標(biāo),許多魚(yú)類(lèi)已經(jīng)適應(yīng)了在高流量脈沖時(shí)節(jié)產(chǎn)卵,高流量脈沖作為魚(yú)類(lèi)洄游產(chǎn)卵的信號(hào),決定魚(yú)類(lèi)產(chǎn)卵的次數(shù)與規(guī)模;b.最大1日、3日平均流量指標(biāo),一般完成一次產(chǎn)卵過(guò)程需要2~3日,這就要求有足夠長(zhǎng)時(shí)間的高流量閾值來(lái)保證產(chǎn)卵的進(jìn)行;c.最大流量出現(xiàn)時(shí)間指標(biāo),該指標(biāo)表征了魚(yú)類(lèi)產(chǎn)卵進(jìn)入高峰期;d.漲水率和落水率指標(biāo),魚(yú)類(lèi)產(chǎn)卵的完成還需要一定的漲、落水率,漲水時(shí)流速加快刺激魚(yú)類(lèi)排卵,而排卵后需要一定的落水過(guò)程才能使魚(yú)卵順利漂向下游;e.6月-8月份平均流量指標(biāo),第二松花江流域魚(yú)類(lèi)的洄游產(chǎn)卵集中在6月-8月份,產(chǎn)卵后魚(yú)卵的漂浮以及發(fā)育都需要一定的流量環(huán)境。

    據(jù)此需求,可從改變度比較大的IHA指標(biāo)中選擇跟魚(yú)類(lèi)相關(guān)的指標(biāo):其中IHA高流量脈沖次數(shù)指標(biāo)改變度為0502,高流量脈沖持續(xù)時(shí)間指標(biāo)改變度為0514,均為中度改變,需要修復(fù);平均最大1日、3日流量指標(biāo)改變度為0292、0302,屬于低度改變,不用考慮修復(fù);最大流量出現(xiàn)時(shí)間指標(biāo)改變度為0452,屬中度改變,需要修復(fù);漲、落水率指標(biāo)的改變度分別為0636、0633,屬于中度改變,需要修復(fù);6月-8月份中只有7月平均流量指標(biāo)改變度為0397,屬于中度改變,需要修復(fù)。經(jīng)以上分析,按改變度大小排序,需要修復(fù)的指標(biāo)有落水率、漲水率、高脈沖流量持續(xù)時(shí)間、高流量脈沖數(shù)、最大流量出現(xiàn)時(shí)間及7月平均流量。

    (3)生態(tài)完整度函數(shù)的構(gòu)建。

    由于所要修復(fù)的指標(biāo)比較多,都作為調(diào)度目標(biāo)難以實(shí)現(xiàn),為此需要由所篩選的指標(biāo)構(gòu)建生態(tài)完整度函數(shù)(式(3)),作為水庫(kù)生態(tài)調(diào)度的目標(biāo)。

    Pe=∑[DD(]6[]i=1[DD)]ωi×μi[JY](3)

    式中:Pe為年生態(tài)完整度;ωi為第i個(gè)生態(tài)指標(biāo)權(quán)重;μi為第i個(gè)生態(tài)指標(biāo)隸屬度值。

    在確定指標(biāo)權(quán)重時(shí),目前的研究均按等權(quán)重處理[510],但各生態(tài)指標(biāo)值的改變程度不同,其修復(fù)時(shí)的重要程度就不同,在確定生態(tài)目標(biāo)時(shí)各指標(biāo)權(quán)重也應(yīng)不同。為此,本文根據(jù)各指標(biāo)值多年平均改變度大小的不同,利用二元比較法確定落水率、漲水率、高流量脈沖持續(xù)時(shí)間、高流量脈沖數(shù)、最大流量出現(xiàn)時(shí)間及7月平均流量指標(biāo)的權(quán)重ωi分別為0244、0244、0163、0163、0105及0081。

    2.4水庫(kù)調(diào)度模型構(gòu)建

    (1)目標(biāo)函數(shù)。

    本文將生態(tài)改變度函數(shù)作為生態(tài)目標(biāo)加入調(diào)度中,由于豐滿(mǎn)水庫(kù)以發(fā)電為主,模型考慮以生態(tài)完整度多年平均值、多年平均發(fā)電量和發(fā)電保證率作為調(diào)度目標(biāo),各目標(biāo)越大越好:

    式中:Vt+1、Vt分別為t+1、t時(shí)段末庫(kù)容;It、Dt、SUt分別為t時(shí)段入庫(kù)水量、泄流量、蒸發(fā)滲漏損失量;Zt為t時(shí)段的水位;Zd為死水位;Zl為汛限水位;Zn為正常高水位;m、n、l為長(zhǎng)春市供水、下游生態(tài)流量、灌溉不發(fā)生破壞的時(shí)段數(shù);Ps、Pt、Pa分別為長(zhǎng)春市供水、下游生態(tài)流量、灌溉的保證率。

    (3)模型求解。

    以豐滿(mǎn)水庫(kù)實(shí)際調(diào)度圖作為初始解,調(diào)度圖中防棄水線、加大出力線和保證出力線上的水位作為模型優(yōu)化變量,采用NSGAⅡ算法對(duì)上述多目標(biāo)優(yōu)化調(diào)度模型進(jìn)行編程求解。

    3結(jié)果與討論

    以豐滿(mǎn)水庫(kù)1946-2009年入庫(kù)流量作為輸入資料,通過(guò)對(duì)模型求解得到108個(gè)帕累托(Pareto)解集,即108個(gè)生態(tài)調(diào)度方案,為了便于分析現(xiàn)有調(diào)度(實(shí)際出流)方案和考慮生態(tài)完整度調(diào)度方案的結(jié)果變化,選出使各目標(biāo)及供水保證率取得最大值和最小值作為代表性方案,對(duì)應(yīng)結(jié)果見(jiàn)表2。

    與現(xiàn)有調(diào)度方案所得結(jié)果相比,本文所建立的生態(tài)調(diào)度模型可將生態(tài)完整度由0437增大到0650以上,至少提高48%,生態(tài)效益改善顯著。為使水庫(kù)泄流的水文情勢(shì)更貼近天然狀態(tài),在非汛期(9月-次年5月),為貼近入流小的情勢(shì),與現(xiàn)有調(diào)度圖中各調(diào)度線相比,考慮生態(tài)完整度所得的調(diào)度圖中各條調(diào)度線在非汛期比較高,使得泄流減少;在汛期(6月-8月),為貼近汛期入流大的情勢(shì),各條調(diào)度線比較低,使得泄流增大。調(diào)度圖變化使發(fā)電量增加07%~26%,這是由于考慮生態(tài)完整度的調(diào)度線在非汛期3月-5月份比較高,即在較高水位時(shí),水庫(kù)仍進(jìn)行限制供水,從而使得水庫(kù)水位抬高,增加發(fā)電量;發(fā)電保證率降低24%~51%,但滿(mǎn)足75%的約束,這主要是3月-5月份泄流受到限制,出力減少造成發(fā)電保證率降低;供水保證率都滿(mǎn)足約束,但均有所降低:長(zhǎng)春市用水保證率降低05%~08%,最小生態(tài)保證率降低32%~48%,灌溉保證率降低18%~80%,同樣由于調(diào)度圖在3月-5月份較高,泄流受到限制所致??傮w來(lái)看,本文所建立的生態(tài)調(diào)度模型,在保證發(fā)電和供水任務(wù)的基礎(chǔ)上,可顯著增大生態(tài)完整度,具有可行性。

    進(jìn)一步分析各生態(tài)指標(biāo)多年平均改善程度,選生態(tài)完整度最大的調(diào)度方案與現(xiàn)有調(diào)度方案相比較,考慮生態(tài)完整度調(diào)度方案中的落水率、漲水率、高流量脈沖持續(xù)時(shí)間、高流量脈沖數(shù)、最大流量出現(xiàn)時(shí)間、7月平均流量指標(biāo)隸屬度分別增加0083、0112、0397、0261、0273、0153,見(jiàn)圖2。各指標(biāo)變化及其生態(tài)意義如下:(1)落水率指標(biāo)隸屬度由0367增加到0450,天然狀態(tài)下落水率為1098 m3/s,由現(xiàn)有調(diào)度方案的175 m3/s降低到1241 m3/s,減緩了現(xiàn)有調(diào)度方案落水過(guò)快所導(dǎo)致的魚(yú)卵擱淺問(wèn)題;(2)漲水率指標(biāo)隸屬度由0364增加到0476,天然狀態(tài)下1173 m3/s,由現(xiàn)有調(diào)度方案的558 m3/s增加到704 m3/s,從而增強(qiáng)了對(duì)魚(yú)類(lèi)產(chǎn)卵刺激,加大了魚(yú)類(lèi)的產(chǎn)卵規(guī)模;(3)高脈沖持續(xù)時(shí)間指標(biāo)隸屬度由0486增加到0883,天然狀態(tài)下為91 d,由現(xiàn)有調(diào)度方案68 d增加到87 d,從而有利于加大魚(yú)類(lèi)的產(chǎn)卵規(guī)模;(4)高流量脈沖數(shù)指標(biāo)隸屬度由0498增加到0759,天然狀態(tài)下發(fā)生247次,由現(xiàn)有調(diào)度方案119次增加到185次,增加魚(yú)類(lèi)產(chǎn)卵的可能性;(5)最大流量出現(xiàn)時(shí)間指標(biāo)隸屬度由0548增加到0821,天然狀態(tài)下出現(xiàn)時(shí)間為第202 d(7月下旬),由現(xiàn)有調(diào)度方案第175 d(6月下旬)推遲到第195 d(7月中旬);(6)7月平均流量指標(biāo)隸屬度由0603增加到0756,天然狀態(tài)下平均流量為9142m3/s,由現(xiàn)有調(diào)度方案下的5397m3/s增加到7783 m3/s,為魚(yú)類(lèi)的產(chǎn)卵以及魚(yú)類(lèi)的漂浮提供了良好的流量環(huán)境。

    兩方案下生態(tài)完整度年際變化過(guò)程見(jiàn)圖3,考慮生態(tài)完整度調(diào)度方案大部分年份的生態(tài)完整度比現(xiàn)有調(diào)度方案大,尤其在枯水年及連續(xù)枯水年更為顯著,生態(tài)完整度增加到04左右。以下簡(jiǎn)要分析對(duì)魚(yú)類(lèi)產(chǎn)卵影響最大的漲水率和高流量脈沖數(shù)的年際變化過(guò)程:(1)如圖4所示,總體上漲水率在豐水年和平水年的隸屬度比枯水年大,需要對(duì)枯水年漲水率指標(biāo)進(jìn)行改善,該模型使?jié)q水率在枯水年由現(xiàn)有調(diào)度方案下的1024 m3/s 增加到3112 m3/s,使魚(yú)類(lèi)在枯水年產(chǎn)卵刺激增強(qiáng),這是由于考慮生態(tài)完整度的調(diào)度線在汛期相對(duì)較低,使得泄流更加容易;(2)由于在汛期調(diào)度線比較低,使得水庫(kù)棄水量及棄水次數(shù)增加,高流量脈沖次數(shù)也隨之增加;但對(duì)于汛期來(lái)水比較大的豐水年,見(jiàn)圖5,在1954年、1956年等特豐水年,汛期泄流次數(shù)多,使得高流量脈沖次數(shù)比天然狀態(tài)下的次數(shù)多,反而造成其相對(duì)于天然狀態(tài)下的指標(biāo)隸屬度降低。

    4結(jié)論

    本文主要針對(duì)豐滿(mǎn)水庫(kù)下游天然水文情勢(shì)改變,產(chǎn)漂流性卵洄游魚(yú)類(lèi)資源衰退問(wèn)題,在滿(mǎn)足發(fā)電和供水的基礎(chǔ)上,進(jìn)一步營(yíng)造了適宜魚(yú)類(lèi)的水文情勢(shì)。通過(guò)篩選水文情勢(shì)改變度大且和魚(yú)類(lèi)洄游產(chǎn)卵相關(guān)的指標(biāo)構(gòu)建生態(tài)完整度函數(shù),建立了考慮生態(tài)完整度的水庫(kù)多目標(biāo)優(yōu)化調(diào)度模型,探討了調(diào)度結(jié)

    果對(duì)發(fā)電和供水的影響,分析了各生態(tài)指標(biāo)多年平均隸屬度的變化情況以及下游生態(tài)流量年際變化情況,主要有以下結(jié)論。

    (1)豐滿(mǎn)水庫(kù)壩址至飲馬河口區(qū)段營(yíng)造適宜水文情勢(shì)的研究表明,該調(diào)度模型在滿(mǎn)足發(fā)電和用水的基礎(chǔ)上,可將下游水文情勢(shì)的生態(tài)完整度提高48%,但會(huì)降低發(fā)電量以及長(zhǎng)春市、下游生態(tài)下泄和灌溉供水保證率。

    (2)各指標(biāo)多年平均隸屬度的分析表明,生態(tài)完整度中落水率、漲水率、高流量脈沖持續(xù)時(shí)間、高流量脈沖數(shù)、最大流量出現(xiàn)時(shí)間、7月平均流量指標(biāo)隸屬度分別增加0083、0112、0397、0261、0273、0153,為魚(yú)類(lèi)的洄游產(chǎn)卵以及魚(yú)卵發(fā)育創(chuàng)造了更有利的條件。

    (3)生態(tài)完整度和指標(biāo)隸屬度年際變化分析表明,生態(tài)完整度在枯水年及連續(xù)枯水年改善程度更佳,生態(tài)完整度達(dá)到04左右;漲水率指標(biāo)在枯水年改善程度比較好,高流量脈沖數(shù)指標(biāo)在枯水年改善程度大、豐水年改善程度小。

    由于第二松花江流域的生態(tài)資料所限,本文僅從流量角度為下游河段營(yíng)造適宜的水流下泄模式,后期[JP3]研究還應(yīng)進(jìn)一步考慮水溫、營(yíng)養(yǎng)物質(zhì)分布等生態(tài)要素。[HJ1.55mm]

    參考文獻(xiàn)(References):

    [1]李翀,彭靜,廖文根.長(zhǎng)江中游四大家魚(yú)發(fā)江生態(tài)水文因子分析及生態(tài)水文目標(biāo)確定[J].中國(guó)水利水電科學(xué)研究院學(xué)報(bào),2006,4(3):170176.(LI C,PENG J,LIAO W G.Study on the ecohydrological factors and flow regime requirement on spawning of four major Chinese carps in the middle reaches of Yangtze River[J].Journal of China Institute of Water Resources and Hydropower Research,2006,4(3):170176.(in Chinese))DOI:10.13244/j.cnki.jiwhr.2006.03.002.

    [2]康玲,黃云燕,楊正祥,等.水庫(kù)生態(tài)調(diào)度模型及其應(yīng)用[J].水利學(xué)報(bào),2010,41(2):134141.(KANG L,HUANG Y Y,YANG Z X.Reservoir ecological operation model and its application[J].Journal of Hydraulic Engineering,2010,41(2):134141.(in Chinese))DOI:10.13243/j.cnki.slxb.2010.02.008.

    [3]李清清,覃暉,陳廣才,等.長(zhǎng)江中游水文情勢(shì)變化及對(duì)魚(yú)類(lèi)的影響分析[J].人民長(zhǎng)江,2012(11):8689.(LI Q Q,QIN H,CHEN G C.Analysis on hydrologic alteration of middle Yangtze River and its impact on fish[J].Yangtze River,2012,43(11):8689.(in Chinese))DOI:10.16232/j.cnki.10014179.2012.11.003.

    [4]湯潔,佘孝云,林年豐,等.生態(tài)需水的理論和方法研究進(jìn)展[J].地理科學(xué),2005,25(3):33673373.(TANG J,SHE X Y,LIN N F.Advances in researches on the theories and methods of ecoenvironmental water demand[J].Scientia Geograph Ica Sinica,2005,25(3):33673373.(in Chinese))

    [5]董哲仁,張晶.洪水脈沖的生態(tài)效應(yīng)[J].水利學(xué)報(bào),2009,40(3):281288.(DONG Z R,ZHANG J.Ecological effect of flood pulses[J].Journal of Hydraulic Engineering,2009,40(3):281288.(in Chinese))

    [6]李翀,廖文根.河流生態(tài)水文學(xué)研究現(xiàn)狀[J].中國(guó)水利水電科學(xué)研究院學(xué)報(bào),2009,7(2):301306.(LI C,LIAO W G.Research status of riverine ecohydrology[J].Journal of China Institute of Water Resources and Hydropower Research,2009,7(2):301306.(in Chinese))DOI:10.13244/j.cnki.jiwhr.2009.02.001.

    [7]RICHTER B D,BAUMGARTNER J V,POWELL J,et al.A method for assessing hydrologic alteration within ecosystems[J].Conservation Biology,1996,10(4):11631174.DOI:10.1046/j.15231739.1996.10041163.x.

    [8]GROWNS J,MARSH N.Characterisation of flow in regulated and unregulated streams in eastern Australia[M].Cooperative Research Centre for Freshwater Ecology,2000.

    [9]JIANPING S,EDWIN H E,WAYLAND E J.Ecohydrologic indicators for rivers of Northern Taiwan[C].ASCE/EWRI World Water and Environmental Resources Congress,Salt Lake City,UT,USA:ASCE,2004:143151.[ZK)]

    [10][ZK(#]徐天寶,彭靜,李翀.葛洲壩水利工程對(duì)長(zhǎng)江中游生態(tài)水文特征的影響[J].長(zhǎng)江流域資源與環(huán)境,2007,16(1):7275.(XU T B,PEN J,LI C.Influence of Gezhouba dam on the ecohydrological characteristics in the middle reach of the Yangtze River[J].Resources and Environment in the Yangtze Basin,2007,16(1):7275.(in Chinese))

    [11]張洪波,黃強(qiáng),彭少明,等.黃河生態(tài)水文評(píng)估指標(biāo)體系構(gòu)建及案例研究[J].水利學(xué)報(bào),2012,43(6):675683.(ZHANG H B,HUANG Q,PENG S M.An ecohydrologic index system and its application on the Yellow River[J].Journal of Hydraulic Engineering,2012,43(6):675683.(in Chinese))DOI:10.13243/j.cnki.slxb.2012.06.005.

    [12] YICHEN E.YANG,XIMING C,EDWIN E H.Identification of hydrologic indicators related to fish diversity and abundance:A data mining approach for fish community analysis[J].Water Resources Research,2008,44(4).DOI:10.1029/2006WR005764.

    [13]SUEN J P,EHEART J W.Reservoir management to balance ecosystem and human needs:Incorporating the paradigm of the ecological flow regime[J].Water Resources Research,2006,42(W034173).DOI:10.1029/2005WR004314,2006.

    [14]楊揚(yáng).考慮生態(tài)需水分析的水庫(kù)調(diào)度研究[D].大連:大連理工大學(xué),2012.(YANG Y.Research on reservoir operation considering ecological water demand analysis[D].Dalian:Dalian University of Technology,2012.(in Chinese))

    [15]王海霞.考慮生態(tài)目標(biāo)的水庫(kù)引水與供水聯(lián)合調(diào)度研究[D].大連:大連理工大學(xué),2015.(WANG H X.Research on the joint operation of water diversion and water supply considering ecological objective[D].Dalian:Dalian University of Technology,2015.(in Chinese))

    [16]王海霞,張弛,周惠成.引水條件下水庫(kù)生態(tài)調(diào)度方法研究[J].水資源與水工程學(xué)報(bào),2016,27(1):16.(WANG H X,ZHANG C,ZHOU H C.Research on ecological operation method under water diversion condition[J].Journal of Water Resources & Water Engineering,2016,27(1):16.(in Chinese))DOI:10.11705 /j.issn.1672643X.2016.01.01.

    [17]楊娜,梅亞?wèn)|,于樂(lè)江.考慮天然水流模式的多目標(biāo)水庫(kù)優(yōu)化調(diào)度模型及應(yīng)用[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,41(1):8589.(YANG N,MEI Y D,YU L J.Multiobjective reservoir optimal operation model considering natural flow regime and its application[J].Journal of Hohai University(Natural Sciences),2013,41(1):8589.(in Chinese))DOI:10.3876 /j.issn.10001980.2013.01.016.

    [18]左其亭,梁士奎.基于水文情勢(shì)分析的閘控河流生態(tài)需水調(diào)控模型研究[J].水力發(fā)電學(xué)報(bào),2016(12):7076.(ZUO Q T,LIANG S K.Regulation model of ecological water demands by sluicecontrolled rivers based on hydrological regime analysis[J].Journal of Hydroelectric Engineering,2016,35(12):7076.(in Chinese))

    [19]潘曉健.考慮魚(yú)類(lèi)多樣性的豐滿(mǎn)水庫(kù)生態(tài)調(diào)度研究[D].大連:大連理工大學(xué),2016.(PAN X J.Study on fengman reservoir ecological operation based on fish diversity[D].Dalian :Dalian University of Technology,2016.(in Chinese))

    [20]楊娜,梅亞?wèn)|,尹志偉.建壩對(duì)下游河道水文情勢(shì)影響RVA評(píng)價(jià)方法的改進(jìn)[J].長(zhǎng)江流域資源與環(huán)境,2010,19(5):560565.(YANG N,MEI Y D,YIN Z W.Impact assessment of dams om the flow regime of lower river by improved RVA[J].Resources and Environment in the Yangtze Basin,2010,19(5):560565.(in Chinese))

    [21]張陵蕾,吳宇雷,張志廣,等.基于魚(yú)類(lèi)棲息地生態(tài)水文特征的生態(tài)流量過(guò)程研究[J].水電能源科學(xué),2015,33(3):1013.(ZHANG L L,WU Y L,ZHANG Z G,et al.Determination of ecological flow regime based on ecohydrological characteristics of fish habitat[J].Water Resources and Power,2015,33(3):1013.(in Chinese))

    [22]駱輝煌,楊青瑞,李倩,等.長(zhǎng)江上游珍稀特有魚(yú)類(lèi)保護(hù)區(qū)魚(yú)類(lèi)生境特征初步研究[J].淡水漁業(yè),2014,44(6):4448.(LUO H H,YANG Q R,LI Q.A preliminary study on fish habitat characteristics in the national nature reserve for the rare and endemic fishes in the upper reaches of the Yangtze River[J].Freshwater Fisheries,2014,44(6):4448.(in Chinese))

    中文字幕av在线有码专区| 色哟哟哟哟哟哟| 亚洲精品在线观看二区| 精品久久久久久久久av| 亚洲av免费高清在线观看| 精品日产1卡2卡| 99精品在免费线老司机午夜| 欧美一区二区精品小视频在线| 在线免费观看的www视频| 卡戴珊不雅视频在线播放| 欧美3d第一页| 国产在视频线在精品| 免费看日本二区| 一区福利在线观看| www.色视频.com| 男女做爰动态图高潮gif福利片| 欧美精品国产亚洲| 一级毛片久久久久久久久女| 欧美xxxx黑人xx丫x性爽| 亚洲精品久久国产高清桃花| 欧美成人一区二区免费高清观看| 好男人在线观看高清免费视频| 久久人人爽人人片av| 精品日产1卡2卡| 国产 一区 欧美 日韩| 午夜视频国产福利| 少妇的逼水好多| 搡老妇女老女人老熟妇| 国国产精品蜜臀av免费| 日本免费a在线| 九九热线精品视视频播放| 露出奶头的视频| 综合色av麻豆| 麻豆精品久久久久久蜜桃| 午夜福利在线观看吧| 禁无遮挡网站| 国产中年淑女户外野战色| 亚洲欧美精品综合久久99| 欧美日韩国产亚洲二区| 国产不卡一卡二| 在线观看美女被高潮喷水网站| 韩国av在线不卡| 真人做人爱边吃奶动态| 亚洲欧美日韩无卡精品| 夜夜爽天天搞| 国产又黄又爽又无遮挡在线| 九九在线视频观看精品| 亚洲熟妇熟女久久| 97超级碰碰碰精品色视频在线观看| 成人精品一区二区免费| 日本与韩国留学比较| 乱人视频在线观看| 97超级碰碰碰精品色视频在线观看| 97超碰精品成人国产| 欧美+亚洲+日韩+国产| a级毛片免费高清观看在线播放| 美女cb高潮喷水在线观看| 国产成人影院久久av| 日本黄色视频三级网站网址| 国产欧美日韩精品一区二区| 久久久久免费精品人妻一区二区| 欧美不卡视频在线免费观看| 大又大粗又爽又黄少妇毛片口| 青春草视频在线免费观看| 亚洲无线观看免费| 久久精品91蜜桃| 欧美日韩精品成人综合77777| 精品国产三级普通话版| 十八禁网站免费在线| 日韩欧美免费精品| 男人和女人高潮做爰伦理| 别揉我奶头 嗯啊视频| 在线国产一区二区在线| 亚洲国产高清在线一区二区三| 国产av在哪里看| avwww免费| 久久久久久久久久黄片| 久久人人爽人人爽人人片va| 国产一区二区亚洲精品在线观看| 亚洲七黄色美女视频| 欧美色视频一区免费| 一级av片app| 天堂av国产一区二区熟女人妻| 亚洲av五月六月丁香网| 亚州av有码| 久久午夜亚洲精品久久| 日韩欧美在线乱码| 精华霜和精华液先用哪个| 最近中文字幕高清免费大全6| 女的被弄到高潮叫床怎么办| 国产精品亚洲一级av第二区| 有码 亚洲区| 精品人妻熟女av久视频| 搞女人的毛片| 日韩欧美一区二区三区在线观看| 久久久色成人| 国产午夜精品论理片| 精品人妻视频免费看| 一区二区三区四区激情视频 | 亚洲人成网站在线播| 精品福利观看| 一进一出好大好爽视频| 简卡轻食公司| 搞女人的毛片| 色综合亚洲欧美另类图片| 禁无遮挡网站| 久久久久免费精品人妻一区二区| 国产一区二区在线观看日韩| 午夜福利视频1000在线观看| 亚洲成av人片在线播放无| 国产精品无大码| 亚洲av第一区精品v没综合| 精品午夜福利视频在线观看一区| 亚洲成人av在线免费| 黑人高潮一二区| 午夜爱爱视频在线播放| 国产成人freesex在线 | 欧美潮喷喷水| 国产 一区精品| 91久久精品国产一区二区成人| 日本在线视频免费播放| 国产极品精品免费视频能看的| 长腿黑丝高跟| 欧美中文日本在线观看视频| 高清午夜精品一区二区三区 | 国产成人aa在线观看| 18+在线观看网站| 欧美日韩在线观看h| 99久久无色码亚洲精品果冻| 99久久精品一区二区三区| 最好的美女福利视频网| 97热精品久久久久久| 天堂网av新在线| 国产精品久久久久久亚洲av鲁大| 日日摸夜夜添夜夜爱| 精品久久久久久久久久久久久| 国产午夜福利久久久久久| 国产一区二区在线av高清观看| 最近2019中文字幕mv第一页| 婷婷亚洲欧美| 国产精品三级大全| 国产v大片淫在线免费观看| 在线a可以看的网站| 国产一区亚洲一区在线观看| 非洲黑人性xxxx精品又粗又长| 日本熟妇午夜| 久久精品影院6| 国产精品一二三区在线看| 亚洲最大成人中文| 有码 亚洲区| 最新在线观看一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 99久久精品热视频| 51国产日韩欧美| 熟妇人妻久久中文字幕3abv| 国产精品人妻久久久影院| 在现免费观看毛片| 国内精品久久久久精免费| 日日摸夜夜添夜夜添av毛片| 午夜亚洲福利在线播放| 久久久久久久久久黄片| 嫩草影院入口| 最近的中文字幕免费完整| 99久久成人亚洲精品观看| 两性午夜刺激爽爽歪歪视频在线观看| 97碰自拍视频| 99在线人妻在线中文字幕| 国内精品宾馆在线| 干丝袜人妻中文字幕| АⅤ资源中文在线天堂| 淫秽高清视频在线观看| 能在线免费观看的黄片| 日韩一区二区视频免费看| 国产一级毛片七仙女欲春2| avwww免费| 最近中文字幕高清免费大全6| 国产单亲对白刺激| 欧美+日韩+精品| 亚洲精品一区av在线观看| 淫秽高清视频在线观看| 三级毛片av免费| 国产精品人妻久久久久久| 寂寞人妻少妇视频99o| 五月玫瑰六月丁香| 夜夜爽天天搞| 国产单亲对白刺激| 看黄色毛片网站| 99久久成人亚洲精品观看| 亚洲国产日韩欧美精品在线观看| 亚洲精品亚洲一区二区| 欧美一区二区国产精品久久精品| 在现免费观看毛片| 亚洲精品乱码久久久v下载方式| 国产精品不卡视频一区二区| 国产爱豆传媒在线观看| 国产精品一区二区三区四区久久| 久久久久九九精品影院| 国产综合懂色| 亚洲成人av在线免费| 国产精品日韩av在线免费观看| 中出人妻视频一区二区| 在线天堂最新版资源| 久久久久久大精品| 欧美一级a爱片免费观看看| 色在线成人网| 最近2019中文字幕mv第一页| 成年女人看的毛片在线观看| 色播亚洲综合网| 日日摸夜夜添夜夜添小说| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美日韩高清专用| 国产精品99久久久久久久久| 97在线视频观看| 国产久久久一区二区三区| 精品人妻偷拍中文字幕| 春色校园在线视频观看| 亚洲色图av天堂| 久久久久久久午夜电影| 亚洲久久久久久中文字幕| 大型黄色视频在线免费观看| 国内精品宾馆在线| 又黄又爽又免费观看的视频| 有码 亚洲区| 国产69精品久久久久777片| 国内精品久久久久精免费| 在现免费观看毛片| 简卡轻食公司| 乱码一卡2卡4卡精品| 国产探花极品一区二区| 少妇人妻一区二区三区视频| 亚洲国产欧美人成| 国产男靠女视频免费网站| 床上黄色一级片| 久久九九热精品免费| 青春草视频在线免费观看| 成人亚洲欧美一区二区av| 少妇被粗大猛烈的视频| 69av精品久久久久久| 中文字幕免费在线视频6| 网址你懂的国产日韩在线| 久久婷婷人人爽人人干人人爱| 日日摸夜夜添夜夜爱| 国产精品野战在线观看| 亚洲av免费在线观看| 男人狂女人下面高潮的视频| 人人妻,人人澡人人爽秒播| 欧美日韩乱码在线| 最近的中文字幕免费完整| 久久久精品94久久精品| 免费看美女性在线毛片视频| 国产黄色视频一区二区在线观看 | 亚洲一区二区三区色噜噜| 国产不卡一卡二| 91在线精品国自产拍蜜月| 热99re8久久精品国产| 看黄色毛片网站| 婷婷精品国产亚洲av在线| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 欧美日韩国产亚洲二区| 2021天堂中文幕一二区在线观| 禁无遮挡网站| 国产视频一区二区在线看| 欧美性感艳星| 久久久久久久亚洲中文字幕| 久久久成人免费电影| 免费看av在线观看网站| 国产男靠女视频免费网站| 亚洲欧美日韩高清专用| 亚洲国产高清在线一区二区三| 成人午夜高清在线视频| 成人永久免费在线观看视频| 91久久精品国产一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美清纯卡通| 国产白丝娇喘喷水9色精品| 露出奶头的视频| 国产 一区 欧美 日韩| 一a级毛片在线观看| 99热6这里只有精品| 日本黄大片高清| 亚洲婷婷狠狠爱综合网| 久久久久久久久久成人| 听说在线观看完整版免费高清| 国产探花在线观看一区二区| 久久久久久久久大av| 免费看日本二区| 亚洲av中文av极速乱| 给我免费播放毛片高清在线观看| 国产精品久久久久久久电影| 一区二区三区免费毛片| 国产白丝娇喘喷水9色精品| 免费观看人在逋| 欧美潮喷喷水| 精品久久久久久久久av| 啦啦啦韩国在线观看视频| 色综合亚洲欧美另类图片| 欧美成人免费av一区二区三区| 日产精品乱码卡一卡2卡三| 简卡轻食公司| 婷婷亚洲欧美| 少妇被粗大猛烈的视频| 乱码一卡2卡4卡精品| 国产色婷婷99| 看黄色毛片网站| 日本五十路高清| 亚洲五月天丁香| 日韩欧美国产在线观看| 在线a可以看的网站| 亚洲内射少妇av| 男人狂女人下面高潮的视频| 国产精品亚洲一级av第二区| 丰满乱子伦码专区| 老司机午夜福利在线观看视频| 国产高清不卡午夜福利| 亚洲av成人av| 亚洲精品456在线播放app| 日韩一区二区视频免费看| 国产又黄又爽又无遮挡在线| 成人毛片a级毛片在线播放| 欧美精品国产亚洲| 中文字幕精品亚洲无线码一区| 国内精品一区二区在线观看| 国产成年人精品一区二区| 热99re8久久精品国产| 伦理电影大哥的女人| 亚洲av美国av| 成年女人看的毛片在线观看| 免费观看的影片在线观看| 欧美日韩综合久久久久久| 夜夜爽天天搞| 亚洲精品一卡2卡三卡4卡5卡| 欧美激情在线99| 亚洲人成网站在线观看播放| 综合色av麻豆| 亚洲国产日韩欧美精品在线观看| 深爱激情五月婷婷| 不卡视频在线观看欧美| 精品免费久久久久久久清纯| 在线国产一区二区在线| 听说在线观看完整版免费高清| 小蜜桃在线观看免费完整版高清| 淫秽高清视频在线观看| 中国美白少妇内射xxxbb| 成人特级av手机在线观看| 国产探花在线观看一区二区| 男女边吃奶边做爰视频| 日韩欧美在线乱码| 我要搜黄色片| 又黄又爽又免费观看的视频| 亚洲精品影视一区二区三区av| 白带黄色成豆腐渣| 欧美成人a在线观看| 老熟妇乱子伦视频在线观看| 激情 狠狠 欧美| 十八禁国产超污无遮挡网站| 亚洲精品一卡2卡三卡4卡5卡| 国产毛片a区久久久久| 少妇被粗大猛烈的视频| 日日干狠狠操夜夜爽| 亚洲图色成人| 亚洲国产高清在线一区二区三| a级毛片免费高清观看在线播放| 蜜桃久久精品国产亚洲av| 1000部很黄的大片| 亚洲七黄色美女视频| 国产精品电影一区二区三区| 亚洲欧美成人综合另类久久久 | 欧美一区二区亚洲| 亚洲最大成人中文| 亚洲高清免费不卡视频| 全区人妻精品视频| 久久九九热精品免费| 寂寞人妻少妇视频99o| 日韩亚洲欧美综合| 日韩一本色道免费dvd| 少妇人妻精品综合一区二区 | 日本免费a在线| 国内精品宾馆在线| 国产色婷婷99| avwww免费| 欧美丝袜亚洲另类| 中文亚洲av片在线观看爽| 成人特级av手机在线观看| 成人鲁丝片一二三区免费| 亚洲专区国产一区二区| 成人特级av手机在线观看| 亚洲第一电影网av| 亚洲精品日韩av片在线观看| 在线免费十八禁| 三级男女做爰猛烈吃奶摸视频| 成人漫画全彩无遮挡| 99热这里只有是精品50| 亚洲精品成人久久久久久| 亚洲四区av| 国产伦在线观看视频一区| av免费在线看不卡| 成年免费大片在线观看| 成年女人看的毛片在线观看| 日本爱情动作片www.在线观看 | 免费av毛片视频| 麻豆精品久久久久久蜜桃| 永久网站在线| 亚洲最大成人av| 一a级毛片在线观看| 最新在线观看一区二区三区| av在线蜜桃| 久久6这里有精品| 人人妻人人看人人澡| 国产精品人妻久久久影院| 一本精品99久久精品77| 嫩草影视91久久| 99riav亚洲国产免费| 伦精品一区二区三区| 久久精品综合一区二区三区| 亚洲国产精品国产精品| 2021天堂中文幕一二区在线观| 乱系列少妇在线播放| 三级男女做爰猛烈吃奶摸视频| 超碰av人人做人人爽久久| 国产综合懂色| 一区福利在线观看| 男人的好看免费观看在线视频| 欧美区成人在线视频| 18禁在线无遮挡免费观看视频 | 国产国拍精品亚洲av在线观看| 亚洲不卡免费看| 久久欧美精品欧美久久欧美| 一本一本综合久久| 亚洲av五月六月丁香网| 久久久久久国产a免费观看| 男女那种视频在线观看| 午夜a级毛片| 狂野欧美白嫩少妇大欣赏| 国产熟女欧美一区二区| 91久久精品国产一区二区成人| 日日摸夜夜添夜夜添av毛片| 综合色丁香网| 国产白丝娇喘喷水9色精品| 最近手机中文字幕大全| 嫩草影视91久久| 成人二区视频| 高清日韩中文字幕在线| 99精品在免费线老司机午夜| 日韩欧美免费精品| 欧美日韩综合久久久久久| 99在线视频只有这里精品首页| 99久久久亚洲精品蜜臀av| 99热全是精品| 亚洲欧美成人综合另类久久久 | 欧洲精品卡2卡3卡4卡5卡区| 免费高清视频大片| 国产成人一区二区在线| 少妇被粗大猛烈的视频| 亚洲第一电影网av| 欧美精品国产亚洲| 色5月婷婷丁香| 欧美激情在线99| 日本黄色视频三级网站网址| 少妇熟女欧美另类| 亚洲美女视频黄频| 联通29元200g的流量卡| 偷拍熟女少妇极品色| 国产白丝娇喘喷水9色精品| 给我免费播放毛片高清在线观看| 日韩av不卡免费在线播放| 人妻少妇偷人精品九色| 久久精品91蜜桃| 国产三级在线视频| 免费人成视频x8x8入口观看| 欧美日韩国产亚洲二区| 久久久精品94久久精品| 国产色婷婷99| 91麻豆精品激情在线观看国产| 精品日产1卡2卡| 国产午夜精品久久久久久一区二区三区 | 国产成人福利小说| 久久久精品大字幕| 99久久中文字幕三级久久日本| 亚洲精品乱码久久久v下载方式| 亚洲av免费在线观看| 免费大片18禁| 亚洲精品一卡2卡三卡4卡5卡| 禁无遮挡网站| 国产三级在线视频| 亚洲天堂国产精品一区在线| 亚洲精品国产成人久久av| 欧美高清成人免费视频www| 91久久精品国产一区二区成人| 午夜福利在线观看免费完整高清在 | 老司机影院成人| 亚洲精品在线观看二区| 国产久久久一区二区三区| 天堂√8在线中文| 国产精品99久久久久久久久| 一本久久中文字幕| 久久欧美精品欧美久久欧美| 色尼玛亚洲综合影院| 天堂√8在线中文| 内射极品少妇av片p| 国产成人影院久久av| 波多野结衣高清无吗| 成熟少妇高潮喷水视频| 九色成人免费人妻av| 国产精品99久久久久久久久| 毛片女人毛片| 精品人妻视频免费看| 欧美另类亚洲清纯唯美| 日韩精品有码人妻一区| 日韩高清综合在线| 性欧美人与动物交配| 91久久精品电影网| 国产成人aa在线观看| 国产免费一级a男人的天堂| 1024手机看黄色片| 午夜福利在线观看免费完整高清在 | 国产亚洲精品av在线| 舔av片在线| 91在线观看av| 嫩草影院精品99| 国产成人一区二区在线| 欧美一区二区亚洲| 天天一区二区日本电影三级| 亚洲av熟女| 亚洲色图av天堂| 能在线免费观看的黄片| 国产高清激情床上av| 真人做人爱边吃奶动态| 亚洲精品影视一区二区三区av| 久久热精品热| 国内久久婷婷六月综合欲色啪| 久久久久久大精品| 精品不卡国产一区二区三区| 国产中年淑女户外野战色| 菩萨蛮人人尽说江南好唐韦庄 | 桃色一区二区三区在线观看| av免费在线看不卡| 国产伦一二天堂av在线观看| 日本一二三区视频观看| 亚洲欧美日韩无卡精品| 国产色婷婷99| 日韩av在线大香蕉| 午夜福利高清视频| 欧洲精品卡2卡3卡4卡5卡区| 综合色av麻豆| 全区人妻精品视频| 一级av片app| 久久精品影院6| 99久久九九国产精品国产免费| 熟女人妻精品中文字幕| 日韩制服骚丝袜av| 中国美白少妇内射xxxbb| 春色校园在线视频观看| 美女大奶头视频| 免费观看精品视频网站| 一区二区三区四区激情视频 | 欧美性猛交黑人性爽| 网址你懂的国产日韩在线| 非洲黑人性xxxx精品又粗又长| 久久久久久伊人网av| 成人永久免费在线观看视频| 1024手机看黄色片| 精品国产三级普通话版| 国产亚洲精品综合一区在线观看| 97热精品久久久久久| 亚洲乱码一区二区免费版| 最近最新中文字幕大全电影3| 少妇裸体淫交视频免费看高清| 久久精品91蜜桃| 亚洲欧美精品综合久久99| 国产成人福利小说| av福利片在线观看| 欧美3d第一页| 变态另类丝袜制服| 人人妻人人澡人人爽人人夜夜 | 日韩av不卡免费在线播放| 亚洲第一电影网av| 亚洲av一区综合| 寂寞人妻少妇视频99o| 中文字幕免费在线视频6| 亚洲无线在线观看| 你懂的网址亚洲精品在线观看 | 麻豆乱淫一区二区| 深夜精品福利| 白带黄色成豆腐渣| 99热这里只有是精品在线观看| 日本精品一区二区三区蜜桃| 亚洲欧美日韩东京热| 国产麻豆成人av免费视频| 国产精品久久久久久精品电影| 久久久久久久久久成人| 国产人妻一区二区三区在| 日韩欧美一区二区三区在线观看| 伦精品一区二区三区| 国产人妻一区二区三区在| 99热网站在线观看| 国产精品乱码一区二三区的特点| 午夜爱爱视频在线播放| 91麻豆精品激情在线观看国产| 国产成人freesex在线 | 久久精品夜夜夜夜夜久久蜜豆| 亚洲成av人片在线播放无| 精品欧美国产一区二区三| 精品国产三级普通话版| 久久久久免费精品人妻一区二区| 99在线视频只有这里精品首页| 99久久精品热视频| 精品少妇黑人巨大在线播放 | 日韩亚洲欧美综合| 九九久久精品国产亚洲av麻豆| 搡老熟女国产l中国老女人| 国产真实乱freesex| 高清日韩中文字幕在线| 又黄又爽又免费观看的视频|