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

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

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

    高志強 丁偉 唐榕 李昱 張弛

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

    關(guān)鍵詞:豐滿水庫;水文情勢;魚類;生態(tài)完整度;生態(tài)調(diào)度

    中圖分類號:TV697.1文獻標(biāo)志碼:A文章編號:

    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研究背景

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

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

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

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

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

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

    2.2漂浮性洄游魚類調(diào)查

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

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

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

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

    首先根據(jù)出流過程對水文情勢改變程度進行評價,采用Richer[7]所提出的水文變化指標(biāo)體系(Indicators of Hydrologic Alteration,IHA)來衡量水文情勢的變化。

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

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

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

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

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

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

    (2)魚類需求指標(biāo)確定。

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

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

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

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

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

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

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

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

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

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

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

    (3)模型求解。

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

    3結(jié)果與討論

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

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

    進一步分析各生態(tài)指標(biāo)多年平均改善程度,選生態(tài)完整度最大的調(diào)度方案與現(xiàn)有調(diào)度方案相比較,考慮生態(tài)完整度調(diào)度方案中的落水率、漲水率、高流量脈沖持續(xù)時間、高流量脈沖數(shù)、最大流量出現(xiàn)時間、7月平均流量指標(biāo)隸屬度分別增加0083、0112、0397、0261、0273、0153,見圖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)度方案落水過快所導(dǎo)致的魚卵擱淺問題;(2)漲水率指標(biāo)隸屬度由0364增加到0476,天然狀態(tài)下1173 m3/s,由現(xiàn)有調(diào)度方案的558 m3/s增加到704 m3/s,從而增強了對魚類產(chǎn)卵刺激,加大了魚類的產(chǎn)卵規(guī)模;(3)高脈沖持續(xù)時間指標(biāo)隸屬度由0486增加到0883,天然狀態(tài)下為91 d,由現(xiàn)有調(diào)度方案68 d增加到87 d,從而有利于加大魚類的產(chǎn)卵規(guī)模;(4)高流量脈沖數(shù)指標(biāo)隸屬度由0498增加到0759,天然狀態(tài)下發(fā)生247次,由現(xiàn)有調(diào)度方案119次增加到185次,增加魚類產(chǎn)卵的可能性;(5)最大流量出現(xiàn)時間指標(biāo)隸屬度由0548增加到0821,天然狀態(tài)下出現(xiàn)時間為第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,為魚類的產(chǎn)卵以及魚類的漂浮提供了良好的流量環(huán)境。

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

    4結(jié)論

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

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

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

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

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

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

    參考文獻(References):

    [1]李翀,彭靜,廖文根.長江中游四大家魚發(fā)江生態(tài)水文因子分析及生態(tài)水文目標(biāo)確定[J].中國水利水電科學(xué)研究院學(xué)報,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]康玲,黃云燕,楊正祥,等.水庫生態(tài)調(diào)度模型及其應(yīng)用[J].水利學(xué)報,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]李清清,覃暉,陳廣才,等.長江中游水文情勢變化及對魚類的影響分析[J].人民長江,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].地理科學(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é)報,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].中國水利水電科學(xué)研究院學(xué)報,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(#]徐天寶,彭靜,李翀.葛洲壩水利工程對長江中游生態(tài)水文特征的影響[J].長江流域資源與環(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]張洪波,黃強,彭少明,等.黃河生態(tài)水文評估指標(biāo)體系構(gòu)建及案例研究[J].水利學(xué)報,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]楊揚.考慮生態(tài)需水分析的水庫調(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)的水庫引水與供水聯(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]王海霞,張弛,周惠成.引水條件下水庫生態(tài)調(diào)度方法研究[J].水資源與水工程學(xué)報,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]楊娜,梅亞東,于樂江.考慮天然水流模式的多目標(biāo)水庫優(yōu)化調(diào)度模型及應(yīng)用[J].河海大學(xué)學(xué)報(自然科學(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]左其亭,梁士奎.基于水文情勢分析的閘控河流生態(tài)需水調(diào)控模型研究[J].水力發(fā)電學(xué)報,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]潘曉健.考慮魚類多樣性的豐滿水庫生態(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]楊娜,梅亞東,尹志偉.建壩對下游河道水文情勢影響RVA評價方法的改進[J].長江流域資源與環(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]張陵蕾,吳宇雷,張志廣,等.基于魚類棲息地生態(tài)水文特征的生態(tài)流量過程研究[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]駱輝煌,楊青瑞,李倩,等.長江上游珍稀特有魚類保護區(qū)魚類生境特征初步研究[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| 久久天躁狠狠躁夜夜2o2o| 一级av片app| 可以在线观看的亚洲视频| 久久精品国产99精品国产亚洲性色| 麻豆一二三区av精品| 一级黄色大片毛片| 欧美日韩国产亚洲二区| 人人妻,人人澡人人爽秒播| 97碰自拍视频| 国产一区二区三区在线臀色熟女| 午夜福利视频1000在线观看| 天堂影院成人在线观看| 亚洲成av人片在线播放无| 1000部很黄的大片| АⅤ资源中文在线天堂| 日韩 亚洲 欧美在线| 乱码一卡2卡4卡精品| 免费无遮挡裸体视频| 不卡一级毛片| 少妇的逼好多水| 午夜亚洲福利在线播放| 午夜精品国产一区二区电影 | 亚洲最大成人中文| 精品久久久久久久久亚洲| 黄色视频,在线免费观看| 国产精品一区www在线观看| 精品久久久久久久久久免费视频| 国产免费一级a男人的天堂| 国产私拍福利视频在线观看| 午夜福利在线观看吧| 最新在线观看一区二区三区| 91在线观看av| 国产一区二区三区av在线 | 最新在线观看一区二区三区| 日韩,欧美,国产一区二区三区 | 两个人的视频大全免费| 精品乱码久久久久久99久播| 好男人在线观看高清免费视频| 观看免费一级毛片| 免费在线观看影片大全网站| 九九热线精品视视频播放| 久久国产乱子免费精品| 男女视频在线观看网站免费| 国产av一区在线观看免费| 亚洲经典国产精华液单| 丰满人妻一区二区三区视频av| 国产精品久久视频播放| 久久人妻av系列| 在线观看66精品国产| 国产精品无大码| a级毛片免费高清观看在线播放| 麻豆国产av国片精品| 久久久久久久亚洲中文字幕| 久久热精品热| 免费看美女性在线毛片视频| 精品99又大又爽又粗少妇毛片| 内射极品少妇av片p| 亚洲欧美日韩无卡精品| 丝袜美腿在线中文| 午夜久久久久精精品| 三级男女做爰猛烈吃奶摸视频| 国产不卡一卡二| 99国产精品一区二区蜜桃av| 中文在线观看免费www的网站| 欧美区成人在线视频| 一进一出好大好爽视频| 国模一区二区三区四区视频| 国产色爽女视频免费观看| 国产熟女欧美一区二区| 变态另类丝袜制服| 欧美色欧美亚洲另类二区| 久久亚洲国产成人精品v| 2021天堂中文幕一二区在线观| 欧美性猛交╳xxx乱大交人| a级毛色黄片| 国产精品电影一区二区三区| 深夜a级毛片| 日韩成人av中文字幕在线观看 | 国产精品,欧美在线| 中国美女看黄片| 亚洲天堂国产精品一区在线| 一区二区三区高清视频在线| 麻豆国产97在线/欧美| 能在线免费观看的黄片| 国产免费男女视频| 黄片wwwwww| 夜夜爽天天搞| 人妻丰满熟妇av一区二区三区| 亚洲成人久久性| 国产一区二区激情短视频| 51国产日韩欧美| 久久久久久久久中文| 国产三级在线视频| 欧美成人一区二区免费高清观看| av专区在线播放| 亚洲精品在线观看二区| 亚洲高清免费不卡视频| 亚洲丝袜综合中文字幕| av天堂在线播放| 国产精品女同一区二区软件| av黄色大香蕉| 12—13女人毛片做爰片一| 色综合站精品国产| 午夜福利成人在线免费观看| 日本-黄色视频高清免费观看| 哪里可以看免费的av片| 最新在线观看一区二区三区| 欧美日韩综合久久久久久| 免费电影在线观看免费观看| 午夜激情欧美在线| 国语自产精品视频在线第100页| 久久精品国产清高在天天线| 51国产日韩欧美| 欧美+日韩+精品| 少妇丰满av| 久久午夜福利片| 国产又黄又爽又无遮挡在线| 国产三级中文精品| 亚洲精品影视一区二区三区av| 国产又黄又爽又无遮挡在线| av在线观看视频网站免费| 狠狠狠狠99中文字幕| 国产av一区在线观看免费| 亚洲中文字幕一区二区三区有码在线看| 少妇的逼好多水| 麻豆乱淫一区二区| 亚洲欧美日韩东京热| www.色视频.com| 亚洲四区av| 成人精品一区二区免费| 亚洲美女视频黄频| 午夜激情欧美在线| 欧美激情久久久久久爽电影| 国产一区二区在线av高清观看| 秋霞在线观看毛片| 亚洲专区国产一区二区| 亚洲最大成人av| 欧美成人精品欧美一级黄| 日本成人三级电影网站| a级毛片a级免费在线| 久久草成人影院| 国产高清不卡午夜福利| 舔av片在线| 人妻久久中文字幕网| 国产黄片美女视频| 色综合站精品国产| 精品日产1卡2卡| 亚洲av五月六月丁香网| a级毛片免费高清观看在线播放| 成年女人毛片免费观看观看9| 波野结衣二区三区在线| 露出奶头的视频| 亚洲av电影不卡..在线观看| 午夜福利在线观看免费完整高清在 | 中文字幕av成人在线电影| 在线观看美女被高潮喷水网站| 内地一区二区视频在线| 婷婷精品国产亚洲av在线| 女人十人毛片免费观看3o分钟| 精品人妻一区二区三区麻豆 | 在线观看午夜福利视频| av在线亚洲专区| 国产精品av视频在线免费观看| 国产精华一区二区三区| 免费看光身美女| 久久精品国产清高在天天线| 成人二区视频| 精品不卡国产一区二区三区| 黄色视频,在线免费观看| 亚洲图色成人| 人人妻人人澡欧美一区二区| 久久九九热精品免费| 日韩欧美在线乱码| 菩萨蛮人人尽说江南好唐韦庄 | 久久精品综合一区二区三区| 国产乱人偷精品视频| 99久久久亚洲精品蜜臀av| 搡老妇女老女人老熟妇| 99精品在免费线老司机午夜| 久久这里只有精品中国| 在线看三级毛片| 亚洲专区国产一区二区| 日本免费a在线| 久久久精品94久久精品| 精品人妻一区二区三区麻豆 | 给我免费播放毛片高清在线观看| 日韩国内少妇激情av| 国产精品,欧美在线| aaaaa片日本免费| 亚洲av二区三区四区| 九色成人免费人妻av| 极品教师在线视频| av天堂中文字幕网| 一个人免费在线观看电影| 中文字幕人妻熟人妻熟丝袜美| 久久热精品热| 午夜福利视频1000在线观看| 成人高潮视频无遮挡免费网站| 波多野结衣高清作品| 久久久国产成人精品二区| 中文亚洲av片在线观看爽| 黄色视频,在线免费观看| 级片在线观看| 18禁黄网站禁片免费观看直播| 美女cb高潮喷水在线观看| 久久久久免费精品人妻一区二区| 男女做爰动态图高潮gif福利片| 久99久视频精品免费| 国产毛片a区久久久久| 欧美精品国产亚洲| 久久精品国产亚洲av天美| 大型黄色视频在线免费观看| 国产亚洲精品久久久com| 亚洲成av人片在线播放无| 香蕉av资源在线| 一级黄片播放器| 国产视频一区二区在线看| 中国国产av一级| 又爽又黄a免费视频| 欧美又色又爽又黄视频| 99久国产av精品国产电影| 天天一区二区日本电影三级| 舔av片在线| 国产大屁股一区二区在线视频| 日韩精品青青久久久久久| 日韩欧美精品免费久久| 少妇熟女aⅴ在线视频| 国产高清视频在线观看网站| 高清毛片免费看| 久久久午夜欧美精品| 久久天躁狠狠躁夜夜2o2o| 91在线观看av| 国产一区二区激情短视频| 蜜桃久久精品国产亚洲av| 日本免费一区二区三区高清不卡| 久久久久国内视频| 日韩精品中文字幕看吧| av在线天堂中文字幕| 十八禁国产超污无遮挡网站| 久久亚洲精品不卡| 成人鲁丝片一二三区免费| 免费在线观看影片大全网站| 春色校园在线视频观看| 三级经典国产精品| 一区二区三区四区激情视频 | 日本免费一区二区三区高清不卡| 久久国内精品自在自线图片| 亚洲经典国产精华液单| 岛国在线免费视频观看| 国产午夜精品论理片| 天堂动漫精品| 乱码一卡2卡4卡精品| 亚洲国产精品sss在线观看| 在线免费十八禁| 亚洲精品一卡2卡三卡4卡5卡| 乱人视频在线观看| 久久久久久久久中文| 精品乱码久久久久久99久播| 99久久成人亚洲精品观看| 午夜福利18| 国产人妻一区二区三区在| 欧美最新免费一区二区三区| 国产黄色小视频在线观看| 啦啦啦韩国在线观看视频| 久久精品国产亚洲av香蕉五月| 国产精品一二三区在线看| 精品熟女少妇av免费看| 成人无遮挡网站| 国产精品乱码一区二三区的特点| 亚洲性夜色夜夜综合| 成人美女网站在线观看视频| 丰满人妻一区二区三区视频av| 老司机福利观看| 在线观看av片永久免费下载| 别揉我奶头~嗯~啊~动态视频| .国产精品久久| 欧美日韩乱码在线| 国产精品不卡视频一区二区| 成人特级av手机在线观看| 又黄又爽又刺激的免费视频.| 亚洲精品成人久久久久久| 国产麻豆成人av免费视频| 国产精品久久久久久久电影| 嫩草影院入口| 色综合色国产| 日韩av在线大香蕉| 日本熟妇午夜| 乱人视频在线观看| 97超碰精品成人国产| 亚洲自偷自拍三级| av在线老鸭窝| 在线免费观看的www视频| 亚洲欧美精品综合久久99| 神马国产精品三级电影在线观看| 人妻少妇偷人精品九色| 免费在线观看影片大全网站| 好男人在线观看高清免费视频| 亚洲成人av在线免费| 女人十人毛片免费观看3o分钟| 午夜免费男女啪啪视频观看 | 成人永久免费在线观看视频| 舔av片在线| 欧美国产日韩亚洲一区| 嫩草影视91久久| 日韩欧美免费精品| 久久久久精品国产欧美久久久| 国产一级毛片七仙女欲春2| 男女边吃奶边做爰视频| 校园春色视频在线观看| 欧美一区二区精品小视频在线| 欧美性感艳星| 又黄又爽又刺激的免费视频.| 国产乱人视频| 一区二区三区免费毛片| 成人精品一区二区免费| 99久久精品热视频| 非洲黑人性xxxx精品又粗又长| 天天一区二区日本电影三级| 久久人人爽人人爽人人片va| 极品教师在线视频| 日韩欧美 国产精品| a级毛片a级免费在线| 晚上一个人看的免费电影| 九九在线视频观看精品| 国产一区二区三区在线臀色熟女| 99久久成人亚洲精品观看| 免费在线观看影片大全网站| 一夜夜www| 国产乱人偷精品视频| 亚洲成人久久性| 日韩欧美免费精品| 亚洲在线自拍视频| 91av网一区二区| 高清毛片免费观看视频网站| 欧美日本视频| 亚洲av美国av| 91精品国产九色| 国产熟女欧美一区二区| 少妇裸体淫交视频免费看高清| 成人美女网站在线观看视频| 国产男人的电影天堂91| 精品一区二区三区av网在线观看| 免费电影在线观看免费观看| 日本一本二区三区精品| 亚洲av美国av| 日韩大尺度精品在线看网址| 变态另类丝袜制服| 看片在线看免费视频| 亚洲国产欧洲综合997久久,| 99热精品在线国产| 国产一区二区激情短视频| 成人欧美大片| 哪里可以看免费的av片| 国产激情偷乱视频一区二区| 淫秽高清视频在线观看| 亚洲欧美成人精品一区二区| 国产精品一区二区三区四区久久| 免费人成在线观看视频色| 亚洲自偷自拍三级| 99热精品在线国产| 丰满乱子伦码专区| 日韩成人伦理影院| 国产精品女同一区二区软件| 国产探花极品一区二区| 日韩av在线大香蕉| 一级黄色大片毛片| 日本一二三区视频观看| 在线天堂最新版资源| 国模一区二区三区四区视频| 欧美三级亚洲精品| 好男人在线观看高清免费视频| 寂寞人妻少妇视频99o| 久久国内精品自在自线图片| 99久久久亚洲精品蜜臀av| 国产淫片久久久久久久久| 麻豆久久精品国产亚洲av| 1024手机看黄色片| 国内少妇人妻偷人精品xxx网站| 成人性生交大片免费视频hd| 亚洲av中文av极速乱| 久久人人爽人人片av| 人人妻,人人澡人人爽秒播| 国产中年淑女户外野战色| 在线观看av片永久免费下载| 免费在线观看成人毛片| 亚洲乱码一区二区免费版| 亚洲成a人片在线一区二区| 一区二区三区高清视频在线| 精品无人区乱码1区二区| 99国产极品粉嫩在线观看| 97超级碰碰碰精品色视频在线观看| 最好的美女福利视频网| 精品午夜福利视频在线观看一区| 婷婷六月久久综合丁香| 十八禁国产超污无遮挡网站| 亚洲内射少妇av| 可以在线观看毛片的网站| 精品久久国产蜜桃| 中文字幕av在线有码专区| 九九热线精品视视频播放| 18+在线观看网站| 深爱激情五月婷婷| 成人高潮视频无遮挡免费网站| 日日摸夜夜添夜夜爱| 国产精品1区2区在线观看.| 亚洲精品久久国产高清桃花| 国内精品美女久久久久久| 麻豆成人午夜福利视频| 蜜臀久久99精品久久宅男| 夜夜夜夜夜久久久久| 国产精品福利在线免费观看| 老司机福利观看| 久久亚洲精品不卡| 国产 一区 欧美 日韩| 97超视频在线观看视频| 可以在线观看毛片的网站| 欧美性猛交黑人性爽| 国产精品不卡视频一区二区| 一区二区三区四区激情视频 | 99久久精品热视频| 成人美女网站在线观看视频| 亚洲精品国产成人久久av| 日日摸夜夜添夜夜添av毛片| 一区二区三区四区激情视频 | 国产蜜桃级精品一区二区三区| 老师上课跳d突然被开到最大视频| 精品乱码久久久久久99久播| 国产激情偷乱视频一区二区| 少妇裸体淫交视频免费看高清| 久久久久国产网址| 性色avwww在线观看| 亚洲一区二区三区色噜噜| 99在线视频只有这里精品首页| 嫩草影视91久久| 欧美3d第一页| 久久99热这里只有精品18| 久久人妻av系列| 国产国拍精品亚洲av在线观看| 色视频www国产| 麻豆久久精品国产亚洲av| 国产麻豆成人av免费视频| 国产男靠女视频免费网站| 日本黄色片子视频| 久久精品国产亚洲av香蕉五月| 真实男女啪啪啪动态图| 欧美一区二区国产精品久久精品| 久久精品国产亚洲av涩爱 | 99热这里只有精品一区| 亚洲av免费在线观看| 成人欧美大片| 日本免费a在线| 亚洲av二区三区四区| 国产亚洲精品久久久com| 久久久久久九九精品二区国产| 亚洲成a人片在线一区二区| 看非洲黑人一级黄片| av黄色大香蕉| 97在线视频观看| 搞女人的毛片| 亚洲精品久久国产高清桃花| ponron亚洲| 久久精品夜色国产| 久久精品国产自在天天线| 中文在线观看免费www的网站| 亚洲熟妇熟女久久| 欧美人与善性xxx| 天天躁日日操中文字幕| 亚洲成人久久性| 97超级碰碰碰精品色视频在线观看| 看非洲黑人一级黄片| 国产 一区 欧美 日韩| 国产一区二区亚洲精品在线观看| 精品少妇黑人巨大在线播放 | 欧美性感艳星| av天堂中文字幕网| 精品人妻熟女av久视频| 婷婷色综合大香蕉| 精品一区二区三区av网在线观看| 亚洲性夜色夜夜综合| 看十八女毛片水多多多| av女优亚洲男人天堂| 蜜臀久久99精品久久宅男| 精品久久久久久成人av| 久久人人爽人人爽人人片va| 熟妇人妻久久中文字幕3abv| 高清毛片免费观看视频网站| 久久欧美精品欧美久久欧美| 99热这里只有精品一区| 国产av不卡久久| 亚洲,欧美,日韩| 亚洲精品久久国产高清桃花| 欧美+亚洲+日韩+国产| 欧美成人免费av一区二区三区| 亚洲一级一片aⅴ在线观看| 麻豆一二三区av精品| 一夜夜www| 中文亚洲av片在线观看爽| 热99re8久久精品国产| 国产不卡一卡二| 国内精品久久久久精免费| 偷拍熟女少妇极品色| 国产高清不卡午夜福利| 久久精品夜色国产| 久久这里只有精品中国| 国产伦精品一区二区三区四那| 亚洲七黄色美女视频| 欧美三级亚洲精品| 成人鲁丝片一二三区免费| 国产私拍福利视频在线观看| 亚洲成av人片在线播放无| 嫩草影院入口| 又粗又爽又猛毛片免费看| 国产v大片淫在线免费观看| 黑人高潮一二区| 尾随美女入室| 久久久久九九精品影院| av天堂在线播放| 在线播放国产精品三级| 淫秽高清视频在线观看| 欧美最新免费一区二区三区| 久久精品夜色国产| 99久久精品热视频| 午夜免费男女啪啪视频观看 | 一进一出抽搐gif免费好疼| 成人鲁丝片一二三区免费| 草草在线视频免费看| 看片在线看免费视频| 亚洲国产精品成人综合色| 99视频精品全部免费 在线| 秋霞在线观看毛片| 一进一出好大好爽视频| 欧美日本视频| 亚洲精品国产av成人精品 | 亚洲国产精品国产精品| 综合色av麻豆| av天堂中文字幕网| 久久久久性生活片| 99久久中文字幕三级久久日本| 中文资源天堂在线| 男插女下体视频免费在线播放| 2021天堂中文幕一二区在线观| 亚洲天堂国产精品一区在线| 欧美性猛交╳xxx乱大交人| 女同久久另类99精品国产91| 97超级碰碰碰精品色视频在线观看| 亚洲av.av天堂| 欧美高清成人免费视频www| 欧美成人a在线观看| 精品人妻偷拍中文字幕| 99热这里只有精品一区| 成人特级av手机在线观看| av福利片在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品,欧美在线| 国产伦一二天堂av在线观看| 免费观看的影片在线观看| 国产 一区精品| 亚洲在线自拍视频| 岛国在线免费视频观看| 22中文网久久字幕| 内地一区二区视频在线| 精品久久国产蜜桃| 老熟妇乱子伦视频在线观看| 亚洲成人久久爱视频| 搡女人真爽免费视频火全软件 | 精品午夜福利在线看| 一个人看视频在线观看www免费| 老司机影院成人| 久久精品国产自在天天线| 麻豆一二三区av精品| 在线观看午夜福利视频| 亚洲美女黄片视频| 啦啦啦观看免费观看视频高清| 日韩av在线大香蕉| 久久欧美精品欧美久久欧美| 国产又黄又爽又无遮挡在线| 搡老妇女老女人老熟妇| 国产精品综合久久久久久久免费| 99久国产av精品| 亚洲国产欧美人成| 亚洲av五月六月丁香网| 简卡轻食公司| 日韩欧美一区二区三区在线观看| 色5月婷婷丁香| 国产欧美日韩一区二区精品| 最新中文字幕久久久久| 久久久久久久久中文| 春色校园在线视频观看| 黄片wwwwww| 18+在线观看网站| 国产精华一区二区三区| 黄片wwwwww| 久久久欧美国产精品| 99久久精品国产国产毛片| 精品久久久久久久久亚洲| 波多野结衣巨乳人妻| 男人的好看免费观看在线视频| 18禁黄网站禁片免费观看直播| 亚洲av中文字字幕乱码综合| 欧美成人一区二区免费高清观看| 日韩欧美精品v在线| 免费人成在线观看视频色| 欧美xxxx黑人xx丫x性爽| 久久久午夜欧美精品| 欧美高清性xxxxhd video| АⅤ资源中文在线天堂| 欧美人与善性xxx| 性色avwww在线观看|