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

    南方季節(jié)性河流河道內(nèi)生態(tài)需水研究
    ——以長江荊南三口為例

    2017-07-07 14:57:18呂殿青
    環(huán)境科學(xué)導(dǎo)刊 2017年4期
    關(guān)鍵詞:生態(tài)

    何 蒙,呂殿青,代 穩(wěn),2

    (1.湖南師范大學(xué)資源與環(huán)境科學(xué)學(xué)院,湖南 長沙 410006;2.六盤水師范學(xué)院,貴州 六盤水 553004)

    南方季節(jié)性河流河道內(nèi)生態(tài)需水研究
    ——以長江荊南三口為例

    何 蒙1,呂殿青1,代 穩(wěn)1,2

    (1.湖南師范大學(xué)資源與環(huán)境科學(xué)學(xué)院,湖南 長沙 410006;2.六盤水師范學(xué)院,貴州 六盤水 553004)

    考慮南方季節(jié)性河流年內(nèi)徑流分布嚴(yán)重不均的問題,依據(jù)1951—2015年長江荊南三口5站實(shí)測原型年徑流量序列,采用Mann-Kendall等方法檢測其徑流序列的突變年份,通過GEV概率密度最大流量、汛期最小輸沙量等方法分別計(jì)算了荊南三口河道內(nèi)生態(tài)需水量、輸沙需水量和水質(zhì)凈化需水量。結(jié)果表明:①水文序列的突變年份判別為1970年,由此將水文序列劃分為變異前(1951—1970年)和變異后(1971—2015年)兩段。②水文變異前,河道內(nèi)年生態(tài)需水量、輸沙需水量和水質(zhì)凈化需水量分別為1239.27×108m3、910.01×108m3、425.70×108m3;水文變異后,河道內(nèi)年生態(tài)需水量、輸沙需水量和水質(zhì)凈化需水量分別為563.32×108m3、501.13×108m3、111.54×108m3。③在季節(jié)上,為保障季節(jié)性河流河道內(nèi)全年均滿足生態(tài)流量,1、2、3、4、11、12月份應(yīng)滿足的生態(tài)需水量為1647.28m3/s,5—10月份應(yīng)滿足的生態(tài)需水量分別為873.87m3/s、2499.59m3/s 、5812.76m3/s 、4346.89m3/s 、3901.18m3/s 、1721.70m3/s。④從綜合角度考慮,水文變異下長江荊南三口季節(jié)性河流河道內(nèi)年生態(tài)需水量為752.71×108m3,年輸沙需水量為910.01×108m3,年水質(zhì)凈化需水量為425.70×108m3。

    生態(tài)需水量;輸沙需水量;水質(zhì)凈化需水量;季節(jié)性河流;荊南三口

    隨著人類活動(dòng)的加劇以及其對生態(tài)環(huán)境影響程度的增強(qiáng),河流生態(tài)需水問題成為當(dāng)前研究的熱點(diǎn)。在全球范圍內(nèi),生態(tài)環(huán)境需水的計(jì)算方法有200多種[1],主要包括水文設(shè)定法、水力參數(shù)法和棲息地模擬法等。胡波等人[2]采用生態(tài)需水系數(shù)-水文參數(shù)耦合模型分別計(jì)算了瀾滄江與紅河的河道內(nèi)生態(tài)需水量與生態(tài)需水流量;倪晉仁等[4]結(jié)合觀測資料與河流生態(tài)需水確定原則對黃河下游河流最小生態(tài)環(huán)境需水量進(jìn)行了初步研究;于魯冀等人[5]利用斜率法、曲率法和多目標(biāo)評價(jià)法3種改進(jìn)濕周法計(jì)算分析了賈魯河不同頻率年河道內(nèi)生態(tài)需水量,但在生態(tài)需水量的界定方面至今尚無統(tǒng)一的概念。許炯心將天然河流劃分為常流河流與季節(jié)河流,常流河流是指河道內(nèi)常年維持一定流量的河流,而季節(jié)性河流是指一年中某一季或某一段較長的時(shí)間干涸無水的河流[6]。依據(jù)地理學(xué)的定義,河流是由河道和河道中足夠流動(dòng)的水構(gòu)成,而足夠流動(dòng)的水則是河流各項(xiàng)功能的基礎(chǔ),也是生態(tài)環(huán)境功能的基礎(chǔ)。一般而言,季節(jié)性河流多出現(xiàn)在干旱地區(qū),但在我國南方豐水地區(qū),也存在季節(jié)性缺水河流。我國南方季節(jié)性缺水河流,水資源總量豐富,但時(shí)空分布跨度大,許多河流在枯水時(shí)期出現(xiàn)了諸如斷流、水環(huán)境容量驟降等導(dǎo)致生態(tài)系統(tǒng)功能損壞的現(xiàn)象[7]。20世紀(jì)50年代以來,長江荊南三口5站在枯水期除松滋河西支的新江口有較少的流量外,其余4站每年均發(fā)生不同程度的斷流現(xiàn)象,尤其是特枯年景,如2006年的沙道觀和管家鋪斷流時(shí)間均達(dá)半年以上[8],使生態(tài)需水量極度貧乏,同時(shí)河道水體含沙量大,嚴(yán)重?fù)p害河道水生態(tài)系統(tǒng)。在氣候變化和人類活動(dòng)的共同作用下,地表水文過程發(fā)生了顯著的變化,而生態(tài)需水不僅是水資源配置的一個(gè)基本理論問題,更是水資源管理中一個(gè)實(shí)際應(yīng)用問題,為保障南方季節(jié)性河流生態(tài)系統(tǒng)的良性發(fā)展,在水文變異基礎(chǔ)上對河流系統(tǒng)的河道內(nèi)生態(tài)需水進(jìn)行研究很有必要。本文在運(yùn)用Mann-Kendall、累積距平法等方法判別水文突變年份的基礎(chǔ)上,從水文變異前后尺度上,計(jì)算長江荊南三口年與月的河道內(nèi)生態(tài)需水量、年輸沙需水量和年水質(zhì)凈化需水量,為制定該地區(qū)水資源管理、水資源開發(fā)利用和優(yōu)化水利工程生態(tài)水調(diào)度方案等提供參考。

    1 研究區(qū)概況與基礎(chǔ)數(shù)據(jù)來源

    荊江是指長江干流枝城至城陵磯河段的別稱,全長約360km。荊南三口(荊江南岸的松滋河、虎渡河、藕池河口)是指分泄長江水進(jìn)入洞庭湖北部地區(qū)所組成的水網(wǎng)[9]。其在豐水季節(jié),水流豐沛,甚至洪水奔涌,且含沙量大;而枯水季節(jié),河水?dāng)嗔?,河床裸露:為典型的南方季?jié)性缺水河流。

    本文選取位于荊南三口河系的5個(gè)主要水文控制站點(diǎn)(即松滋口的新江口站和沙道觀、虎渡口的彌陀寺站、藕池河的管家鋪站和康家崗站)的總年徑流量與總輸沙量數(shù)據(jù)作為判別變異點(diǎn)、計(jì)算年生態(tài)需水量、輸沙需水量的原始數(shù)據(jù),采用其月徑流量觀測數(shù)據(jù)作為逐月計(jì)算生態(tài)需水量的原始數(shù)據(jù)。數(shù)據(jù)來自于湖南省水利水電勘測設(shè)計(jì)研究總院和湖南省水文資源勘測網(wǎng)站上的水情日報(bào)表。部分缺失數(shù)據(jù)通過與相鄰的水文站水文序列建立回歸關(guān)系進(jìn)行插補(bǔ)(相關(guān)系數(shù)達(dá)0.8以上)。

    2 研究方法

    在氣候和人類活動(dòng)的雙重影響下,水文序列往往發(fā)生變異,使得前后水文序列的總體分布不一致,不符合水文資料的一致性要求。于是選用長江荊南三口五站1951—2015年總徑流量作為考察水文變異的指標(biāo),通過水文變異點(diǎn)來劃分水文變異前后水文序列,并據(jù)此利用荊南三口總年徑流量、月均流量、月輸沙量計(jì)算水文變異前后長江荊南三口河道內(nèi)生態(tài)需水量、輸沙需水量和水質(zhì)凈化需水量,以便作比較分析。

    2.1 水文突變年份的檢驗(yàn)

    近60多年來長江中上游地區(qū)修建了荊江系統(tǒng)裁彎、葛洲壩水利樞紐、三峽水庫等眾多水利工程,從而改變了荊江河段的水沙輸送量及分配時(shí)間。受其影響,荊南三口一方面水沙量減少,另一方面在枯水季節(jié),除松滋河西支的新江口外,其余4站每年均發(fā)生不同程度的斷流現(xiàn)象,且斷流天數(shù)呈逐期增加趨勢(表1)。因此需要檢驗(yàn)水文序列的異點(diǎn),便于比較水文變異前后荊南三口河道內(nèi)生態(tài)需水量、輸沙需水量和水質(zhì)凈化需水量。

    考慮單一檢測方法有其適用流域的局限性,本研究分別運(yùn)用M-K檢驗(yàn)法[10-11]、累計(jì)距平法[9]及滑動(dòng)t檢驗(yàn)法[12]聯(lián)合判斷長江荊南三口5個(gè)控制站(新江口站、沙道觀、彌陀寺、管家鋪和康家崗)近65年的總年徑流(5站年徑流量之和)突變年份。然后將以上三種方法的檢測結(jié)果進(jìn)行綜合分析,確定荊南三口河系徑流序列的變異點(diǎn)。

    表1 荊南三口各時(shí)期平均斷流天數(shù)[8]

    2.2 河道內(nèi)生態(tài)需水量計(jì)算方法

    對于河流生態(tài)系統(tǒng)而言,流量、水位、流速、水質(zhì)等都是限制生物發(fā)展的重要生態(tài)因子。其中流量不僅與河流生態(tài)系統(tǒng)的限制性因子有關(guān),而且是其系統(tǒng)的重要限制因子,因此可將流量作為計(jì)算河道內(nèi)生態(tài)需水量的原始數(shù)據(jù)。在生物生長與發(fā)展的過程中會(huì)存在其最適生存環(huán)境,而生物在最適環(huán)境中,具有個(gè)數(shù)最多、生長速度最快、生長得最好的特征,隨著生態(tài)環(huán)境偏離最優(yōu)狀態(tài),該生物雖然可以出現(xiàn)生長,但生長速度減慢,且當(dāng)生態(tài)因子繼續(xù)偏離最適值,超過耐受上限,該物種就不再生長,個(gè)數(shù)減少,直到消失。而長期的自然選擇使得生物適應(yīng)出現(xiàn)頻率較高的環(huán)境因子,所以頻率最大的流量最適合河道內(nèi)生態(tài)系統(tǒng)[13]。因此,本文采用廣義極值(GEV)分布的頻率最大處流量作為河道內(nèi)生態(tài)需水量,并運(yùn)用P-Ⅲ型分布90%頻率流量與其結(jié)果進(jìn)行比較。

    2.2.1 GEV概率密度最大流量法

    GEV概率密度分布函數(shù)為:

    (1)

    (2)

    式中:σ是尺度參數(shù),影響分布尺度;μ是位置參數(shù),影響分布水平位置;ξ是形狀參數(shù),影響分布尾部情況,ξ=0、ξ>0和ξ<0分別對應(yīng)Gumbel分布、Weibull分布和Frechet分布。

    運(yùn)用線形矩法對GEV分布進(jìn)行參數(shù)估計(jì)。將觀測到的n個(gè)樣本值的變量X從小到大排列。xj為序列中第j個(gè)值。前三階的L-矩l1、l2、l3及L-偏態(tài)系數(shù)t3為:

    (3)

    (4)

    (5)

    (6)

    由式(3)~(6)可估算出GEV參數(shù):

    (7)

    (8)

    (9)

    式中:Γ(1+ξ)為伽馬函數(shù)。

    對(1)式進(jìn)行最值分析,得頻率最大處流量(即各月河道內(nèi)生態(tài)流量)為:

    (10)

    2.2.2P-Ⅲ型90%保證率法

    P-Ⅲ型曲線與水文現(xiàn)象符合較好,在水文計(jì)算中廣為采用,P-Ⅲ型頻率分配(密度)曲線有很復(fù)雜的數(shù)學(xué)形式,此處不做詳細(xì)贅述,僅將簡易計(jì)算方式列出:將觀測到的n個(gè)樣本值的變量X從大到小排列,xi為序列中第i個(gè)值:

    (11)

    (12)

    經(jīng)過數(shù)學(xué)推導(dǎo),得90%頻率各月河道內(nèi)生態(tài)流量為:

    (13)

    式中:Cv為離差系數(shù),表示系列的相對離散程度;Ki為模比系數(shù);Cs為偏差系數(shù),它反映隨機(jī)變量偏離均值的情況;Ф為離均系數(shù),此處根據(jù)90%頻率和偏差系數(shù)Cs查附表Ф值(福斯特-雷布京)表,再用(13)式計(jì)算x90%值即各月河道內(nèi)生態(tài)流量。

    2.3 河道內(nèi)輸沙需水和水質(zhì)凈化需水量計(jì)算方法

    水流和泥沙既是塑造河流形態(tài)的物質(zhì)基礎(chǔ),又是水生動(dòng)植物生長繁殖的必備水文條件。由于長江荊南三口具有斷流期長、含沙量大、水污染嚴(yán)重等水文特征,因此在研究該地區(qū)河道內(nèi)生態(tài)需水量時(shí),還應(yīng)當(dāng)考慮河流輸沙需水量和水質(zhì)凈化需水量。輸沙需水量的計(jì)算采用汛期最小輸沙量法[14],計(jì)算公式如下:

    (14)

    式中:W為輸沙需水量(108m3),St為多年平均輸沙量(t),Cmax為多年最大月平均含沙量的平均值,Cij為第i年j月的月平均含沙量(t/m3);n為統(tǒng)計(jì)年數(shù)。

    水質(zhì)凈化生態(tài)需水量就是改善河道內(nèi)水質(zhì)所需要的基本水量,我國一般采用近10a最枯月平均流量或90%最枯月平均流量法[15]??紤]出現(xiàn)水文變異情況時(shí),水文變異前后水文序列長度未必相等,且荊南三口枯水期長期斷流,全年徑流基本集中在汛期,因此采用90%保證率汛期最小月平均流量法計(jì)算荊南三口的自凈需水量。

    3 結(jié)果與分析

    3.1 水文序列突變年份

    將1951—2015年長江荊南三口5站總年徑流量的Mann-Kendall檢驗(yàn)統(tǒng)計(jì)量分為兩段,并將結(jié)果與α=0.05對應(yīng)的正態(tài)分布區(qū)間Uα∈(-1.96,1.96)進(jìn)行對比分析。結(jié)果表明,1951—2015年荊南三口徑流變化過程在整體上呈現(xiàn)減少趨勢,可信度達(dá)到95%。UF、UB統(tǒng)計(jì)量如圖2,在1951—2015年荊南三口徑流變化過程中,出現(xiàn)突變年份是1969年和1988年。

    為驗(yàn)證用M-K檢驗(yàn)法檢驗(yàn)的荊南三口徑流量變化過程中的突變年份,本文采用徑流量累計(jì)距平法和滑動(dòng)t檢驗(yàn)法作進(jìn)一步檢驗(yàn)。由累計(jì)距平法所得結(jié)果可知,荊南三口徑流過程在1951—2015年發(fā)生兩次突變(圖3),其突變年份分別是1970年和1980年。由滑動(dòng)t檢驗(yàn)法所檢驗(yàn)的結(jié)果為1970年(圖4)。

    將以上三種方法檢測出的突變年份進(jìn)行比較分析,結(jié)果表明:用M-K檢測的1969年前后為突變點(diǎn)的可信度達(dá)到95%,1980年為突變年份的可信度不顯著,而累計(jì)距平法和滑動(dòng)t檢驗(yàn)法檢測出的突變年份均包含1970年,與M-K法檢測出的突變年份基本無異。由此,可以判定長江荊南三口徑流序列在1970年前后發(fā)生突變。而在此期間下荊江正在實(shí)施系統(tǒng)裁彎取直工程,說明此工程對荊南三口徑流量改變的影響較為顯著。

    3.2 河道內(nèi)生態(tài)需水量

    根據(jù)以上變異年份將所研究水文序列分為變異前(1951—1970年)與變異后(1971—2015年)兩時(shí)段。一般而言,計(jì)算河道內(nèi)生態(tài)需水時(shí),只考慮水文變異點(diǎn)前水文序列,但由于荊南三口的河流為季節(jié)性河流,全年徑流量大部集中在5—10月份,而1、2、3、4、11及12月份,除松滋河西支的新江口外,其余4站每年均發(fā)生不同程度的斷流現(xiàn)象。河流斷流則說明河道內(nèi)沒有流動(dòng)的水,根本無法滿足河道內(nèi)生物生存、繁衍所需要的水量。更為重要的是,以荊江系統(tǒng)裁彎工程、葛洲壩水利樞紐和三峽水庫等水利工程建設(shè)為代表的人類活動(dòng)已對荊南三口水沙環(huán)境造成了不可逆轉(zhuǎn)的影響。從這一角度出發(fā),本研究課題組認(rèn)為,在對比分析水文變異前后生態(tài)需水量變化的基礎(chǔ)上,將水文變異后的計(jì)算結(jié)果作為荊南三口季節(jié)性河流河道內(nèi)生態(tài)需水量更有現(xiàn)實(shí)意義。由于水文現(xiàn)象既具有必然性, 又具有隨機(jī)性, 因此計(jì)算河道內(nèi)生態(tài)流量除運(yùn)用GEV概率密度最大流量法外, 還要與P-Ⅲ型90%保證率法進(jìn)行比較。

    3.2.1 GEV概率密度最大流量法計(jì)算結(jié)果

    運(yùn)用GEV概率密度最大流量法,對荊南三口水文變異前(1951—1970年)1月平均流量序列進(jìn)行矩參數(shù)估計(jì),得到變異前1月份參數(shù):

    ξ= -0.16173,σ= 54.9334,μ=100.7246

    由(10)式計(jì)算變異前1月份頻率最大流量為:

    x=σ/ξ[(ξ+1)-ξ-1]+μ=110.27(m3/s)

    以同樣的方法可分別計(jì)算出變異前后其它月份頻率最大流量(表2)及頻率最大年徑流量。

    水文變異前后頻率最大年徑流量分別為1239.27×108m3和563.32×108m3。如圖5顯示,水文變異前(1951—1970年的20個(gè)年份)年徑流量未滿足河道內(nèi)生態(tài)需水量(1239.27×108m3)的年份有1953、1959、1960、1966、1967、1969與1970等7個(gè)年份;滿足河道內(nèi)生態(tài)需水量的有13個(gè)年份,且這13年的年均徑流量為1514.40×108m3。水文變異后(1971—2015年的45個(gè)年份),年徑流量未滿足河道內(nèi)生態(tài)需水量(563.32×108m3)的年份有1986、1992、1994、1997、2001、2002、2003、2004、2006、2007、2008、2009、2010、2011、2013、2014及2015年等17個(gè)年份;而滿足河道內(nèi)生態(tài)需水量的有28個(gè)年份,其年均徑流量為771.60×108m3,遠(yuǎn)大于其生態(tài)需水量563.32×108m3。水文變異后,生態(tài)徑流減少量較大,但水利工程對荊南三口水環(huán)境的影響不可逆轉(zhuǎn),由此課題組認(rèn)為,在按照年徑流量對荊南三口實(shí)施水利工程調(diào)控時(shí),應(yīng)在滿足河道內(nèi)生態(tài)需水量563.32×108m3的基礎(chǔ)上,再盡可能達(dá)到771.60×108m3的水量, 這將對維持河道水生態(tài)系統(tǒng)健康更為有利。

    表2 長江荊南三口各月河道內(nèi)生態(tài)需水量及頻率滿足表

    相比于北方缺水地區(qū)的季節(jié)性河流而言,南方豐水地區(qū)的季節(jié)性河流總量豐富、需水量要求大,但時(shí)空分布跨度大,許多河流在枯水時(shí)期出現(xiàn)不同程度的斷流現(xiàn)象,季節(jié)性缺水特征明顯。考慮荊南三口全年徑流高度集中在5—10月份,且枯水期長期斷流等因素,在計(jì)算整個(gè)序列的河道內(nèi)生態(tài)需水量時(shí)還應(yīng)當(dāng)逐月進(jìn)行討論。如果在水文序列中某年某月平均流量大于或等于對應(yīng)月的生態(tài)需水流量,則認(rèn)為該年該月滿足了河道內(nèi)生態(tài)需水量,否則該月未滿足河道內(nèi)生態(tài)需水量。于是將各月滿足生態(tài)需水量的年份除以研究的總年數(shù),求出滿足生態(tài)流量的頻率(表2)。如表2所示,水文變異前后生態(tài)需水量的滿足頻率基本保持在50%~70%,且各月份滿足頻率在水文變異前后差別不大,甚至某些月份的滿足頻率在水文變異后還略大于變異前,如6、8、9等月份。這是由于水文變異后各月份月均流量減少,河道內(nèi)生態(tài)需水量的標(biāo)準(zhǔn)不同。如7月份的生態(tài)流量從變異前的10992.67m3/s下降至變異后的5812.76m3/s,變異前后相差5000m3/s以上,可見水文變異后,荊南三口河道內(nèi)的生態(tài)流量大幅度減少,生態(tài)環(huán)境嚴(yán)重惡化。

    從圖6長江荊南三口逐月生態(tài)流量曲線可以看出,水文變異前后1、2、3、4、12月份的生態(tài)流量變化較小,5—11月份生態(tài)流量變化較大,尤其是7—9月份,減少量均在4500m3/s以上,這主要是三峽水庫調(diào)洪調(diào)度所致。1—5月份和11—12月份的生態(tài)流量在水文變異前后都未達(dá)到其平均生態(tài)流量,且1、2、3、4、12月份在水文變異前的生態(tài)流量還未及變異后的平均生態(tài)流量(1647.28m3/s)。可見,1、2、3、4、12月份的流量無論在變異前還是在變異后都無法達(dá)到生態(tài)需水的要求。因此,為保障季節(jié)性河流河道內(nèi)全年均滿足生態(tài)流量,應(yīng)將水文變異后的平均生態(tài)流量作為枯水期1、2、3、4、11、12月份的生態(tài)需水量,水文變異后5—10月份的頻率最大流量作為相對應(yīng)月份的生態(tài)需水量。由以上分析可知,在按照月均流量對荊南三口實(shí)施水利工程調(diào)控時(shí),1、2、3、4、11、12月份應(yīng)滿足的生態(tài)需水量為1647.28m3/s,5—10月份應(yīng)滿足的生態(tài)需水量依次為873.87m3/s、2499.59m3/s、5812.76m3/s、4346.89m3/s、3901.18m3/s、1721.70m3/s。

    依照上述生態(tài)流量數(shù)據(jù),將其轉(zhuǎn)化為生態(tài)徑流量得全年的生態(tài)需水量為752.71×108m3,其結(jié)果比按頻率最大年徑流量所得結(jié)果563.32×108m3更大,且與水文變異后28個(gè)滿足生態(tài)需水量年份的年均徑流量771.60×108m3基本無異。由此可見,對于南方河流徑流分布不均的季節(jié)性河流而言,按月考慮河道內(nèi)生態(tài)需水量更為合理。

    3.2.2 與P-Ⅲ型90%保證率法比較

    運(yùn)用(11)、(12)式得到水文變異前1月份的模比Ki,求得離差系數(shù)Cv的值為0.53,運(yùn)用適線法(圖7)選取合適的Cs值,并查表得出Ф值,運(yùn)用(13)式求出變異前1月份生態(tài)流量x90%=95.68m3/s。以同樣的方式計(jì)算出P-Ⅲ型90%保證率法中變異前后各月份的生態(tài)流量(表3)。通過比較GEV概率密度最大流量法與P-Ⅲ型90%保證率法的計(jì)算結(jié)果可以看出,無論是水文變異前還是變異后,GEV概率密度最大流量法的計(jì)算結(jié)果皆大于90%保證率法的計(jì)算結(jié)果,說明其計(jì)算結(jié)果也適用于90%保證率下的徑流。通過采用Kolmogorov-Smirnov(K-S)法檢驗(yàn)兩種方法計(jì)算結(jié)果的擬合程度,水文變異前后兩者的漸近顯著性皆為0.996,遠(yuǎn)遠(yuǎn)大于0.05,說明兩樣本都服從于同一分布。此外,將P-Ⅲ型90%保證率法、GEV概率密度最大流量法與Tennant法[16]進(jìn)行對比,發(fā)現(xiàn)兩種方式所得河道內(nèi)生態(tài)需水量都在Tennant法“好”等級以上。因此,采用GEV概率密度最大流量法計(jì)算河道內(nèi)生態(tài)需水量是可行的。

    表3 兩種生態(tài)需水量計(jì)算方法對比 (m3/s)

    3.3 輸沙需水量和水質(zhì)凈化需水量

    長江荊南三口河系豐水期平均流量較大,輸沙也主要靠豐水期較大的流量完成。根據(jù)公式(14)計(jì)算得出水文變異前年輸沙需水量為910.01×108m3,小于其生態(tài)需水量(即1239.27×108m3);水文變異后年輸沙需水量為501.13×108m3,也小于其生態(tài)需水量(563.32×108m3)。據(jù)統(tǒng)計(jì)水文變異前僅有1969年的年徑流量未滿足輸沙需水量(910.01×108m3),其它19個(gè)年份均滿足輸沙需水,且滿足年份的年均徑流量為1389.84×108m3;變異后,有1994、1997、2001、2006、2009、2011、2013、2015等8個(gè)年份無法滿足輸沙需水(501.13×108m3),其它37個(gè)年份均滿足輸沙需水量,且這37個(gè)年份的年均徑流量為715.47×108m3。根據(jù)90%汛期最枯月平均流量法可得荊南三口水文變異前后年水質(zhì)凈化需水量分別為425.70×108m3、111.54×108m3。據(jù)統(tǒng)計(jì),變異前后所有年份的徑流量皆能滿足水質(zhì)凈化需水,且變異前后各年份年均徑流量分別為1365.52×108m3、652.74×108m3。

    由表4可以看出,在生態(tài)需水量、輸沙需水量和水質(zhì)凈化需水量這三類需水量中,無論是變異前還是變異后生態(tài)需水量的值都最大。但從輸沙需水量的角度來看,由于泥沙是固體物質(zhì),輸沙對水量的需求也應(yīng)當(dāng)更大;而更重要是,水文變異后荊南三口水量的減少趨勢明顯,但沙量的減少卻并不穩(wěn)定,且變異前后其含沙量都較大,使得水文變異后泥沙淤積更為嚴(yán)重。因此,為保障河道內(nèi)具有良好的輸沙能力,將水文變異前的輸沙需水量(910.01×108m3)作為水文變異下荊南三口季節(jié)性河流河道內(nèi)的輸沙需水量更適宜。從水質(zhì)凈化需水量的角度考慮,隨著荊南三口斷流天數(shù)逐期增加,其水質(zhì)標(biāo)準(zhǔn)也大大降低。改革開放以來,隨著工業(yè)的迅速發(fā)展以及城市化進(jìn)程的不斷推進(jìn),使得水體污染日益嚴(yán)重,所需凈化的水量大,在水文變異后的水資源分布條件下,實(shí)現(xiàn)荊南三口水質(zhì)保護(hù)目標(biāo)的可能性極低。所以,為達(dá)到水質(zhì)保護(hù)目標(biāo),本文將水文變異前的水質(zhì)凈化需水量(425.70×108m3)作為水文變異下荊南三口季節(jié)性河流河道內(nèi)的輸沙需水量。

    表4 長江荊南三口三類需水量比較 (108m3)

    4 結(jié)論與討論

    根據(jù)長江荊南三口季節(jié)性河流河道內(nèi)徑流演變特征以及生態(tài)需水、輸沙需水、水質(zhì)凈化需水的要求,在水文變異條件下長江荊南三口季節(jié)性河流河道內(nèi)年生態(tài)需水量為752.71×108m3,輸沙需水量為910.01×108m3,水質(zhì)凈化需水量為425.70×108m3。

    (1)由于水文現(xiàn)象既具有必然性, 又具有隨機(jī)性, 為保障河道內(nèi)生態(tài)系統(tǒng)的健康發(fā)展,其生態(tài)流量應(yīng)要適應(yīng)不同流量保證率下的環(huán)境,因此在計(jì)算河道內(nèi)生態(tài)流量時(shí),運(yùn)用GEV概率密度最大流量法與P-Ⅲ型90%保證率法進(jìn)行比較,發(fā)現(xiàn)二者的擬合程度較高,且計(jì)算結(jié)果都在Tennant法好等級以上,從而認(rèn)為運(yùn)用GEV概率密度最大流量法計(jì)算河道內(nèi)生態(tài)需水量可行。

    (2)從年際上看,長江荊南三口水文變異前(1951—1970年的20個(gè)年份)河道內(nèi)年生態(tài)需水量為1239.27×108m3,年徑流量未滿足河道內(nèi)生態(tài)需水量的年份有7個(gè);水文變異后(1971—2015年的45個(gè)年份),河道內(nèi)生態(tài)需水量為563.32×108m3,年徑流量未滿足河道內(nèi)生態(tài)需水量的有17個(gè)年份。通過前后對比分析發(fā)現(xiàn),兩者的生態(tài)徑流相差較大,但以水利工程為代表的人類活動(dòng)已對荊南三口河道內(nèi)生態(tài)環(huán)境造成了不可逆轉(zhuǎn)的影響。鑒于此,本課題組認(rèn)為把水文變異后的生態(tài)徑流量作為水文變異下荊南三口季節(jié)性河流的生態(tài)需水量更為合理。從季節(jié)上看, 1—5月份以及11、12月份的生態(tài)流量在水文變異前后都未達(dá)到其平均生態(tài)流量,且變異前1、2、3、4、12月份的生態(tài)流量還未及水文變異后的平均生態(tài)流量(1647.28m3/s)。由此可見,1、2、3、4、12月份的流量無論在變異前還是在變異后都無法達(dá)到生態(tài)需水的要求。為保障季節(jié)性河流河道內(nèi)全年均滿足生態(tài)流量,可將水文變異后的平均生態(tài)流量作為枯水期1、2、3、4、11、12月份的生態(tài)需水量,水文變異后5—10月份的頻率最大流量作為相對應(yīng)月份的生態(tài)需水量。據(jù)此,1、2、3、4、11、12月份應(yīng)滿足的生態(tài)需水量為1647.28 m3/s,5—10月份應(yīng)滿足的生態(tài)需水量分別為873.87m3/s、2499.59m3/s、5812.76m3/s、4346.89m3/s、3901.18m3/s、1721.70m3/s。將月均生態(tài)流量轉(zhuǎn)化為生態(tài)徑流量得全年的生態(tài)需水量為752.71×108m3,其結(jié)果比按頻率最大年徑流量所得結(jié)果563.32×108m3更大,且與28個(gè)滿足生態(tài)需水量年份的年均徑流量771.60×108m3基本無異。因此,通過綜合考慮認(rèn)為,將752.71×108m3作為長江荊南三口季節(jié)性河流河道內(nèi)生態(tài)需水量更具合理性。

    (3)據(jù)計(jì)算,長江荊南三口水文變異前年生態(tài)需水量、輸沙需水量、水質(zhì)凈化需水量分別為1239.27×108m3、910.01×108m3、425.70×108m3;水文變異后三種需水量分別為563.32×108m3、501.13×108m3、111.54×108m3。通過比較荊南三口水文變異前后年生態(tài)需水量、輸沙需水量、水質(zhì)凈化需水量可知,各需水量前后相差都較大,生態(tài)需水量的值最高。從輸沙需水量來看,一方面,水文變異后荊南三口水量的減少趨勢明顯,但沙量的減少卻并不穩(wěn)定,且變異前后其含沙量都較大,使得水文變異后泥沙淤積更為嚴(yán)重;另一方面由于泥沙是固體物質(zhì),其對水量的需求也更大。因此,將水文變異前的輸沙需水量(910.01×108m3)作為水文變異下荊南三口季節(jié)性河流河道內(nèi)的輸沙需水量更為合理。從水質(zhì)凈化需水量的角度分析,荊南三口隨著斷流天數(shù)的增加,水質(zhì)標(biāo)準(zhǔn)也大大降低,且水文變異后的水質(zhì)凈化需水量結(jié)果較小,水污染嚴(yán)重,實(shí)現(xiàn)荊南三口水質(zhì)保護(hù)目標(biāo)的可能性極低。在這種狀況下,本文將水文變異前的水質(zhì)凈化需水量(425.70×108m3)作為水文變異下荊南三口季節(jié)性河流河道內(nèi)的輸沙需水量。

    (4)水流是河道內(nèi)生態(tài)環(huán)境的物質(zhì)基礎(chǔ),又是河道內(nèi)生物繁殖的必備水文條件,而河道內(nèi)生態(tài)需水量是受多種環(huán)境、生物和生態(tài)等因子共同影響的,由于資料獲取困難等原因,本文僅用流量對荊南三口河道內(nèi)生態(tài)需水量、輸沙需水量、水質(zhì)凈化需水量進(jìn)行了計(jì)算。如何有效獲取與量化多種環(huán)境、生物和生態(tài)等因子,并將其融入河道內(nèi)生態(tài)需水量的計(jì)算當(dāng)中,還需要進(jìn)一步研究。

    [1]趙長森,劉昌明,夏軍,等. 閘壩河流河道內(nèi)生態(tài)需水研究--以淮河為例[J]. 自然資源學(xué)報(bào),2008(3):400-411.

    [2]胡波,鄭艷霞,翟紅娟, 等. 生態(tài)需求流量與河道內(nèi)生態(tài)需水量計(jì)算研究——以瀾滄江、紅河為例[J]. 長江科學(xué)院院報(bào),2015(3):99-106.

    [4]倪晉仁,金玲,趙業(yè)安.黃河下游河流最小生態(tài)環(huán)境需水量初步研究[J].水利學(xué)報(bào),2002(10):1-7.

    [5]張強(qiáng),李劍鋒,陳曉宏, 等. 水文變異下的黃河流域生態(tài)流量[J]. 生態(tài)學(xué)報(bào),2011,31(17):4826-4834.

    [6]高志勇,馮佳睿,周川閩,等. 干旱氣候環(huán)境下季節(jié)性河流沉積特征-以庫車河剖面下白堊統(tǒng)為例[J]. 沉積學(xué)報(bào),2014(6):1060-1071.

    [7]陽書敏,邵東國,沈新平.南方季節(jié)性缺水河流生態(tài)環(huán)境需水量計(jì)算方法[J].水利學(xué)報(bào),2005,36(11):1341-1346.

    [8]帥紅,李景保,何霞,劉雯, 等. 環(huán)境變化下長江荊南三口徑流變化特征檢測與歸因分析[J]. 水土保持學(xué)報(bào),2016,30(1):83-88.

    [9]李景保,羅中海,葉亞亞, 等. 三峽水庫運(yùn)行對長江荊南三口水文和生態(tài)的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào),2016,27(4):1285-1293.

    [10]Mann H B.Non-parametric test against trend[J].econometica,1945,13(3):245-259.

    [11]Kendall M G,Gibbons J D.Rank Correlation Methods [M].Londan:Edward Arnold Oxford Universitu Press,1975.

    [12]謝平,陳廣才,雷紅富, 等. 水文變異診斷系統(tǒng)[J]. 水力發(fā)電學(xué)報(bào),2010,29(1):85-91.

    [13]李劍鋒,張強(qiáng),陳曉宏,等. 考慮水文變異的黃河干流河道內(nèi)生態(tài)需水研究[J]. 地理學(xué)報(bào),2011,66(1):99-110.

    [14]李麗娟,鄭紅星. 海灤河流域河流系統(tǒng)生態(tài)環(huán)境需水量計(jì)算[J]. 地理學(xué)報(bào),2000,55(4):495-500.

    [15]海河流域水生態(tài)回復(fù)研究課題組.海河流域生態(tài)環(huán)境需水量計(jì)算方法與標(biāo)準(zhǔn)[R]. 2003.

    [16]Tennant DL.Instream flow regimens for fish,wildlife,recreation,and related environmental resources[J].American Fisheries Society,1976,1(4):6-10.

    [17]彭濤,陳曉宏,王高旭, 等. 基于濕周法的河道最小生態(tài)需水量多目標(biāo)評價(jià)模型[J]. 水利水電科技進(jìn)展,2012,(5):6-10.

    [18]李嘉,王玉蓉,李克鋒, 等. 計(jì)算河段最小生態(tài)需水的生態(tài)水力學(xué)法[J]. 水利學(xué)報(bào),2006,37(10):1169-1174.

    [19]劉劍宇,張強(qiáng),顧西輝. 水文變異條件下鄱陽湖流域的生態(tài)流量[J]. 生態(tài)學(xué)報(bào),2015,35(16):5477-5485.

    [20]馮夏清,章光新,尹雄銳. 基于生態(tài)保護(hù)目標(biāo)的太子河下游河道生態(tài)需水量計(jì)算[J]. 環(huán)境科學(xué)學(xué)報(bào),2010,30(7):1466-1471.

    [21]謝永宏,李峰,陳心勝. 洞庭湖最小生態(tài)需水量研究[J]. 長江流域資源與環(huán)境,2012(1):64-70.

    [22]張強(qiáng),李劍鋒,陳曉宏, 等. 水文變異下的黃河流域生態(tài)流量[J]. 生態(tài)學(xué)報(bào),2011,31(17):4826-4834.

    Research on Ecological Water Demand within the Main Rivers of Seasonal Streams at the ThreeDiversion Entrances in ChangjiangJingjiang River

    HE Meng1,LU Dian-qing1,DAI Wen1'2

    (1.College of Resources and Environment Science, Hunan Normal University, Changsha Hunan 410006, China)

    Considering the yearly uneven runoff distribution in the southern seasonal streams, this article adopted Mann-Kendall method to identify the years when runoff sequence suddenly changes based on the annual runoff time series of prototype tests coming from the three divisions and five stations of ChangjiangJingjiang River from 1951 to 2015. Through GEV′s maximum flow of probability density, the minimum sediment runoff of flood season and other methods, ecological water demand, sediment water demand, and water purification water demand within the rivers of Jingjiang three diversions were separately calculated. The result showed that the sudden change of hydrologic sequence was in 1970. Therefore, hydrologic sequence could be divided into ex ante variation from 1951 to 1970 and ex post variation from 1971 to 2015. Before the hydrologic variation, annual ecological water demand, sediment water demand and water purification water demand were respectively 1 239.27×108 m3, 910.01×108 m3, and 425.70×108 m3. After the hydrologic variation, annual ecological water demand, sediment water demand, and water purification water demand are respectively 563.32×108 m3, 501.13×108 m3, and 111.54×108 m3. In seasons, in order to make ecological water demand within seasonal streams meet the whole year's ecological flow, with the ecological flow of 1 647.28 m3/s in January, February, March, April, November, and December must be ensured. From May to October the ecologic flows were separately 873.87 m3/s, 2 499.59 m3/s, 5812.76 m3/s, 4346.89 m3/s, 3901.18 m3/s, and 1721.70 m3/s. From the general perspectives, ecological water demand within the main rivers of seasonal streams at the three diversion entrances in ChangjiangJingjiang River under hydrologic variation was 752.71×108 m3, the sediment water demand was 910.01×108 m3, and water purification water demand was 425.70×108 m3

    ecological flow requirement; water requirements for transporting sediment; water requirement for self-purification of the stream; seasonal water-deficient river; three diversion entrances in ChangjiangJingjiang River

    2017-01-17

    國家自然科學(xué)基金項(xiàng)目(41571100);湖南省重點(diǎn)學(xué)科建設(shè)項(xiàng)目(2011001);湖南省教育廳重點(diǎn)項(xiàng)目(16A129)。

    何蒙(1992-),女,湖南瀏陽人,碩士研究生,主要從事生態(tài)需水研究。

    呂殿青(1975-),女,博士,教授,博士生導(dǎo)師,主要從事土壤生態(tài)與環(huán)境變化研究。

    X52

    A

    1673-9655(2017)04-0022-09

    猜你喜歡
    生態(tài)
    “生態(tài)養(yǎng)生”娛晚年
    住進(jìn)呆萌生態(tài)房
    生態(tài)之旅
    生態(tài)之旅
    生態(tài)之旅
    大營鎮(zhèn)生態(tài)雞
    貴茶(2019年3期)2019-12-02 01:46:32
    生態(tài)之旅
    潤豐達(dá) 微平衡生態(tài)肥
    茶葉通訊(2017年2期)2017-07-18 11:38:40
    生態(tài)保護(hù) 有你有我
    “知”與“信”:《逃逸行為》的生態(tài)自我
    女人爽到高潮嗷嗷叫在线视频| 精品国产一区二区三区四区第35| 成人免费观看视频高清| 18禁黄网站禁片午夜丰满| 亚洲精品国产色婷婷电影| 女人爽到高潮嗷嗷叫在线视频| 精品电影一区二区在线| 欧美激情极品国产一区二区三区| 十八禁人妻一区二区| 免费久久久久久久精品成人欧美视频| 亚洲av成人一区二区三| 韩国精品一区二区三区| 成年人免费黄色播放视频| 亚洲人成电影观看| 人人妻人人添人人爽欧美一区卜| 老司机靠b影院| 免费在线观看完整版高清| 亚洲五月天丁香| 91麻豆av在线| 色综合婷婷激情| 99久久久亚洲精品蜜臀av| 午夜福利一区二区在线看| 香蕉丝袜av| 日本wwww免费看| 丰满迷人的少妇在线观看| 一二三四在线观看免费中文在| 99re在线观看精品视频| 电影成人av| 精品日产1卡2卡| 一本大道久久a久久精品| 日本五十路高清| 91成年电影在线观看| 欧美日韩一级在线毛片| 91九色精品人成在线观看| 天天躁夜夜躁狠狠躁躁| 在线免费观看的www视频| 18美女黄网站色大片免费观看| 亚洲男人的天堂狠狠| 热99国产精品久久久久久7| 国产三级在线视频| 久久精品国产亚洲av高清一级| 成人国语在线视频| 欧美精品啪啪一区二区三区| 日韩 欧美 亚洲 中文字幕| 香蕉丝袜av| 欧美黄色淫秽网站| 女性被躁到高潮视频| 成人三级黄色视频| 久久中文字幕人妻熟女| 久久久久国产一级毛片高清牌| 亚洲午夜理论影院| 国产精品久久久av美女十八| 欧美在线一区亚洲| 级片在线观看| 身体一侧抽搐| 激情在线观看视频在线高清| 日韩欧美一区视频在线观看| 日韩大码丰满熟妇| 成人18禁高潮啪啪吃奶动态图| 麻豆成人av在线观看| 婷婷六月久久综合丁香| av超薄肉色丝袜交足视频| 一级片'在线观看视频| 老司机亚洲免费影院| 一级a爱片免费观看的视频| 日本免费a在线| 97超级碰碰碰精品色视频在线观看| 在线观看免费高清a一片| 亚洲全国av大片| 露出奶头的视频| 精品久久久久久电影网| 久久精品亚洲熟妇少妇任你| 一a级毛片在线观看| cao死你这个sao货| 国产人伦9x9x在线观看| 成人av一区二区三区在线看| 俄罗斯特黄特色一大片| 国产精品久久电影中文字幕| 变态另类成人亚洲欧美熟女 | 久久99一区二区三区| 色综合欧美亚洲国产小说| 国产精品乱码一区二三区的特点 | 在线看a的网站| 久久草成人影院| www日本在线高清视频| 国产真人三级小视频在线观看| 亚洲国产精品一区二区三区在线| 成人18禁在线播放| 欧美亚洲日本最大视频资源| 黑丝袜美女国产一区| 精品久久久久久电影网| 一边摸一边抽搐一进一出视频| 亚洲片人在线观看| 午夜精品在线福利| 国产熟女xx| 成人三级黄色视频| 久久久久九九精品影院| 黄色毛片三级朝国网站| 国产在线观看jvid| 丝袜人妻中文字幕| 麻豆一二三区av精品| 亚洲精品在线观看二区| 久久亚洲真实| 国产高清国产精品国产三级| 男人舔女人的私密视频| а√天堂www在线а√下载| 欧美另类亚洲清纯唯美| 麻豆久久精品国产亚洲av | 欧美日韩瑟瑟在线播放| 88av欧美| 十分钟在线观看高清视频www| 免费在线观看视频国产中文字幕亚洲| 精品国产美女av久久久久小说| 十八禁人妻一区二区| 怎么达到女性高潮| 最好的美女福利视频网| 一边摸一边抽搐一进一小说| 精品一区二区三区av网在线观看| av中文乱码字幕在线| 亚洲av第一区精品v没综合| 亚洲精品国产一区二区精华液| av免费在线观看网站| 亚洲,欧美精品.| 波多野结衣一区麻豆| 久久国产亚洲av麻豆专区| 91成人精品电影| 日本a在线网址| √禁漫天堂资源中文www| 国产视频一区二区在线看| www.999成人在线观看| 久99久视频精品免费| 欧美日韩国产mv在线观看视频| 免费少妇av软件| 中亚洲国语对白在线视频| 丰满迷人的少妇在线观看| 国产精品一区二区精品视频观看| 在线观看一区二区三区激情| 国产精品九九99| 18禁裸乳无遮挡免费网站照片 | 欧美日韩av久久| 久久人人97超碰香蕉20202| 国产蜜桃级精品一区二区三区| 叶爱在线成人免费视频播放| 亚洲色图 男人天堂 中文字幕| 精品日产1卡2卡| netflix在线观看网站| 欧美亚洲日本最大视频资源| 精品一区二区三区视频在线观看免费 | 午夜91福利影院| 国产成人欧美| 国产91精品成人一区二区三区| 99国产精品一区二区蜜桃av| 国产单亲对白刺激| 动漫黄色视频在线观看| 欧美日韩乱码在线| 亚洲精品久久午夜乱码| 国产一区在线观看成人免费| av福利片在线| 久久国产精品影院| 丝袜美腿诱惑在线| 久久久久久久久免费视频了| 欧美在线黄色| 操出白浆在线播放| 变态另类成人亚洲欧美熟女 | 一进一出抽搐动态| 国产精品98久久久久久宅男小说| 91成年电影在线观看| 成人亚洲精品av一区二区 | 很黄的视频免费| 手机成人av网站| 中文字幕人妻丝袜一区二区| 婷婷丁香在线五月| 999精品在线视频| 丝袜美腿诱惑在线| 脱女人内裤的视频| 午夜福利一区二区在线看| 精品熟女少妇八av免费久了| 久热爱精品视频在线9| 亚洲熟妇熟女久久| 高清在线国产一区| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产成人免费| 一本大道久久a久久精品| 国产91精品成人一区二区三区| 午夜老司机福利片| 国产精品亚洲一级av第二区| 国产免费现黄频在线看| 搡老乐熟女国产| 欧美久久黑人一区二区| 在线观看一区二区三区激情| 日本一区二区免费在线视频| 精品福利观看| 五月开心婷婷网| 久久草成人影院| 不卡一级毛片| 国产精品久久视频播放| 久久精品亚洲精品国产色婷小说| 国产欧美日韩一区二区三| 精品国产一区二区久久| 精品一区二区三区视频在线观看免费 | 69精品国产乱码久久久| 欧美精品亚洲一区二区| 琪琪午夜伦伦电影理论片6080| 婷婷丁香在线五月| ponron亚洲| 久久久久久久精品吃奶| 色播在线永久视频| 国产主播在线观看一区二区| 久久久久久亚洲精品国产蜜桃av| 黄网站色视频无遮挡免费观看| 久久久精品欧美日韩精品| 久久伊人香网站| ponron亚洲| 一个人观看的视频www高清免费观看 | 免费观看精品视频网站| 另类亚洲欧美激情| 极品人妻少妇av视频| 国产精品1区2区在线观看.| 成人免费观看视频高清| 人人澡人人妻人| 无遮挡黄片免费观看| 高清毛片免费观看视频网站 | 国产一区二区三区综合在线观看| 亚洲精品国产一区二区精华液| 999久久久国产精品视频| 男女下面插进去视频免费观看| 亚洲 欧美一区二区三区| 精品一区二区三卡| 国产国语露脸激情在线看| 欧美精品亚洲一区二区| 51午夜福利影视在线观看| 夜夜躁狠狠躁天天躁| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 又黄又爽又免费观看的视频| 俄罗斯特黄特色一大片| 真人做人爱边吃奶动态| 88av欧美| 国产人伦9x9x在线观看| 精品一品国产午夜福利视频| 99久久精品国产亚洲精品| 丰满饥渴人妻一区二区三| 51午夜福利影视在线观看| 国产激情欧美一区二区| 深夜精品福利| 熟女少妇亚洲综合色aaa.| 韩国精品一区二区三区| 69精品国产乱码久久久| 久久精品91蜜桃| 高清黄色对白视频在线免费看| 国产成人免费无遮挡视频| 国产精品久久久人人做人人爽| 极品人妻少妇av视频| 日韩中文字幕欧美一区二区| 99国产精品一区二区蜜桃av| 亚洲专区字幕在线| 老汉色av国产亚洲站长工具| 最近最新中文字幕大全电影3 | 岛国在线观看网站| videosex国产| 可以在线观看毛片的网站| 丰满的人妻完整版| 校园春色视频在线观看| 久久人人精品亚洲av| 如日韩欧美国产精品一区二区三区| www日本在线高清视频| 亚洲av日韩精品久久久久久密| 亚洲欧洲精品一区二区精品久久久| 村上凉子中文字幕在线| 午夜老司机福利片| 妹子高潮喷水视频| 这个男人来自地球电影免费观看| 精品福利观看| 国产成人精品无人区| 看黄色毛片网站| 99国产精品免费福利视频| 精品午夜福利视频在线观看一区| 久久人人97超碰香蕉20202| 国产精品偷伦视频观看了| 欧美激情久久久久久爽电影 | 国产成人欧美在线观看| 国产亚洲精品第一综合不卡| 精品无人区乱码1区二区| 国产99白浆流出| 一级,二级,三级黄色视频| 国产精品久久久久成人av| videosex国产| 美女福利国产在线| 成人免费观看视频高清| 国产成人免费无遮挡视频| 欧美黄色片欧美黄色片| 夫妻午夜视频| 一区福利在线观看| av网站在线播放免费| 波多野结衣高清无吗| 中文欧美无线码| 国产乱人伦免费视频| 久久人妻熟女aⅴ| 看黄色毛片网站| 老汉色av国产亚洲站长工具| 午夜免费鲁丝| 深夜精品福利| 淫妇啪啪啪对白视频| 黄色女人牲交| 夫妻午夜视频| 婷婷六月久久综合丁香| √禁漫天堂资源中文www| 丁香欧美五月| 国产成人欧美在线观看| 亚洲国产精品999在线| 极品人妻少妇av视频| 午夜视频精品福利| 久99久视频精品免费| 国产一区二区三区视频了| 国产黄色免费在线视频| 亚洲欧美精品综合一区二区三区| 一级片免费观看大全| 欧美日韩国产mv在线观看视频| 亚洲精品国产精品久久久不卡| 久久热在线av| 宅男免费午夜| 91大片在线观看| av视频免费观看在线观看| 成人亚洲精品av一区二区 | 精品第一国产精品| 美女高潮喷水抽搐中文字幕| 精品久久久久久久久久免费视频 | 夜夜爽天天搞| 国产亚洲欧美在线一区二区| 日本wwww免费看| 手机成人av网站| 18禁国产床啪视频网站| a在线观看视频网站| 国产激情欧美一区二区| 中出人妻视频一区二区| 在线观看www视频免费| 亚洲五月色婷婷综合| 亚洲成人免费av在线播放| 在线天堂中文资源库| bbb黄色大片| 欧美精品亚洲一区二区| 91精品三级在线观看| 久久久精品国产亚洲av高清涩受| 中文字幕高清在线视频| 日韩 欧美 亚洲 中文字幕| 啪啪无遮挡十八禁网站| 在线观看午夜福利视频| 亚洲精品国产色婷婷电影| www.自偷自拍.com| 如日韩欧美国产精品一区二区三区| 国产亚洲精品一区二区www| av在线播放免费不卡| 精品熟女少妇八av免费久了| 亚洲专区字幕在线| 精品国产一区二区三区四区第35| 嫁个100分男人电影在线观看| 搡老乐熟女国产| av中文乱码字幕在线| 一级片免费观看大全| 久久 成人 亚洲| 午夜福利欧美成人| 在线观看日韩欧美| 久久久久国产精品人妻aⅴ院| 男女床上黄色一级片免费看| 国产成人欧美在线观看| 超碰97精品在线观看| 久久精品国产亚洲av高清一级| xxxhd国产人妻xxx| 一a级毛片在线观看| 人人妻人人澡人人看| 国产精品98久久久久久宅男小说| 久久久久久久午夜电影 | 国产欧美日韩一区二区三区在线| 久久久久久久久久久久大奶| 美女福利国产在线| 国产欧美日韩一区二区三| 男女高潮啪啪啪动态图| 美女高潮喷水抽搐中文字幕| 啦啦啦在线免费观看视频4| 19禁男女啪啪无遮挡网站| 日日爽夜夜爽网站| 亚洲性夜色夜夜综合| 搡老熟女国产l中国老女人| 99国产综合亚洲精品| 美国免费a级毛片| 香蕉丝袜av| 两性午夜刺激爽爽歪歪视频在线观看 | 在线播放国产精品三级| 免费在线观看完整版高清| 91老司机精品| 成人免费观看视频高清| 香蕉久久夜色| 人成视频在线观看免费观看| 亚洲一区中文字幕在线| 一个人观看的视频www高清免费观看 | 老司机深夜福利视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲午夜精品一区,二区,三区| 国产精品久久久久久人妻精品电影| 一边摸一边抽搐一进一出视频| 亚洲国产欧美一区二区综合| 精品国产超薄肉色丝袜足j| 后天国语完整版免费观看| 12—13女人毛片做爰片一| 免费在线观看黄色视频的| 波多野结衣高清无吗| 女人精品久久久久毛片| av在线播放免费不卡| 高清在线国产一区| 天堂√8在线中文| 一级毛片高清免费大全| 久久人妻熟女aⅴ| 激情在线观看视频在线高清| 黄片大片在线免费观看| 日本免费一区二区三区高清不卡 | 美国免费a级毛片| 欧美成人性av电影在线观看| 麻豆成人av在线观看| 中文欧美无线码| 国产精品综合久久久久久久免费 | 成人18禁在线播放| 99国产精品免费福利视频| 国产1区2区3区精品| 久久久精品欧美日韩精品| 一本大道久久a久久精品| 99国产精品一区二区蜜桃av| 午夜两性在线视频| 18禁黄网站禁片午夜丰满| 国产野战对白在线观看| www.精华液| 国产91精品成人一区二区三区| 精品国产乱子伦一区二区三区| 亚洲av日韩精品久久久久久密| 好看av亚洲va欧美ⅴa在| 高清黄色对白视频在线免费看| 少妇 在线观看| 久久九九热精品免费| 亚洲专区中文字幕在线| 老司机午夜福利在线观看视频| 亚洲片人在线观看| 夜夜躁狠狠躁天天躁| 午夜福利免费观看在线| 一区二区日韩欧美中文字幕| 亚洲欧洲精品一区二区精品久久久| 美女高潮喷水抽搐中文字幕| 丝袜美足系列| 亚洲国产精品一区二区三区在线| 亚洲中文字幕日韩| 国产精品秋霞免费鲁丝片| 一级a爱片免费观看的视频| 法律面前人人平等表现在哪些方面| 欧美老熟妇乱子伦牲交| 男男h啪啪无遮挡| 久久人妻福利社区极品人妻图片| 国产精品一区二区精品视频观看| 一区在线观看完整版| 国产亚洲精品第一综合不卡| 国产精品久久视频播放| 国产精品98久久久久久宅男小说| 巨乳人妻的诱惑在线观看| 国产成人精品无人区| 色精品久久人妻99蜜桃| 在线观看66精品国产| 黄片大片在线免费观看| 亚洲精品国产区一区二| 丝袜人妻中文字幕| 一边摸一边抽搐一进一出视频| 中文字幕人妻丝袜制服| 好看av亚洲va欧美ⅴa在| 叶爱在线成人免费视频播放| 国产1区2区3区精品| www.999成人在线观看| 亚洲欧洲精品一区二区精品久久久| 一a级毛片在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产xxxxx性猛交| 欧美成人免费av一区二区三区| 一a级毛片在线观看| 久久久久精品国产欧美久久久| 午夜福利一区二区在线看| 国产aⅴ精品一区二区三区波| 国产野战对白在线观看| 国产成人欧美在线观看| 久久欧美精品欧美久久欧美| 欧美色视频一区免费| 丰满饥渴人妻一区二区三| 老鸭窝网址在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久午夜电影 | 久热爱精品视频在线9| 日本a在线网址| 手机成人av网站| 久久精品成人免费网站| 亚洲欧美日韩另类电影网站| 国产精华一区二区三区| 妹子高潮喷水视频| 国产精品偷伦视频观看了| 亚洲精品国产区一区二| 日韩高清综合在线| 久久久久久免费高清国产稀缺| 午夜精品在线福利| 少妇粗大呻吟视频| 极品人妻少妇av视频| 亚洲自拍偷在线| 精品日产1卡2卡| 婷婷六月久久综合丁香| 国产片内射在线| 国产一卡二卡三卡精品| 免费av中文字幕在线| 一区福利在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲精品第一综合不卡| 日本vs欧美在线观看视频| 新久久久久国产一级毛片| 亚洲欧美精品综合久久99| 亚洲第一av免费看| 久久久精品国产亚洲av高清涩受| 成人av一区二区三区在线看| 国产精品成人在线| 妹子高潮喷水视频| 制服诱惑二区| 少妇粗大呻吟视频| 欧美日韩一级在线毛片| 黄色 视频免费看| 午夜精品久久久久久毛片777| 欧美日韩黄片免| 一区在线观看完整版| ponron亚洲| 热99re8久久精品国产| 午夜免费成人在线视频| 国产97色在线日韩免费| 欧美av亚洲av综合av国产av| www国产在线视频色| 啦啦啦免费观看视频1| 国产精品久久电影中文字幕| avwww免费| 国产高清videossex| 欧美老熟妇乱子伦牲交| 亚洲在线自拍视频| 久久香蕉国产精品| 黑人操中国人逼视频| 99国产综合亚洲精品| 高清欧美精品videossex| 亚洲精品中文字幕一二三四区| 九色亚洲精品在线播放| 精品久久久久久电影网| 欧美老熟妇乱子伦牲交| 免费在线观看视频国产中文字幕亚洲| 后天国语完整版免费观看| 欧美在线黄色| 久久国产亚洲av麻豆专区| 亚洲五月天丁香| www.自偷自拍.com| 搡老乐熟女国产| 日韩欧美免费精品| 丝袜在线中文字幕| 亚洲男人的天堂狠狠| 久久久久久大精品| 亚洲欧美激情在线| 在线av久久热| 99在线人妻在线中文字幕| 女性被躁到高潮视频| 久久香蕉激情| 久久精品国产亚洲av香蕉五月| 精品国产国语对白av| a在线观看视频网站| 天天躁狠狠躁夜夜躁狠狠躁| av国产精品久久久久影院| 交换朋友夫妻互换小说| 国产精品影院久久| 国产在线精品亚洲第一网站| 亚洲欧美一区二区三区黑人| 天天躁夜夜躁狠狠躁躁| 成人影院久久| 18禁美女被吸乳视频| 999久久久国产精品视频| 久久久久精品国产欧美久久久| 国产有黄有色有爽视频| 欧美日韩黄片免| 国产激情久久老熟女| 成人亚洲精品一区在线观看| 午夜两性在线视频| 精品午夜福利视频在线观看一区| 琪琪午夜伦伦电影理论片6080| 精品少妇一区二区三区视频日本电影| 亚洲一区高清亚洲精品| 亚洲熟妇中文字幕五十中出 | 香蕉国产在线看| 精品国产一区二区三区四区第35| 看黄色毛片网站| 91国产中文字幕| av在线播放免费不卡| 亚洲午夜理论影院| 国产av又大| 搡老岳熟女国产| 一区二区三区国产精品乱码| 大香蕉久久成人网| 精品一区二区三卡| 国产亚洲精品综合一区在线观看 | 无遮挡黄片免费观看| 亚洲色图综合在线观看| 久久精品成人免费网站| 亚洲第一av免费看| 亚洲中文日韩欧美视频| 日韩大码丰满熟妇| 国产一区在线观看成人免费| 麻豆一二三区av精品| 91国产中文字幕| 亚洲精品久久成人aⅴ小说| 好看av亚洲va欧美ⅴa在| 别揉我奶头~嗯~啊~动态视频| 他把我摸到了高潮在线观看| 欧美性长视频在线观看| 国产麻豆69| 老汉色av国产亚洲站长工具|