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

    1990—2013年湖南四水入洞庭湖汛期徑流量的變化規(guī)律

    2015-06-06 12:40:43劉培亮毛德華鄧曉磊
    水資源保護(hù) 2015年4期
    關(guān)鍵詞:四水桃江湘潭

    劉培亮,毛德華,周 慧,2,鄧曉磊

    (1.湖南師范大學(xué)資源與環(huán)境科學(xué)學(xué)院,湖南 長沙 410081; 2.湖南省水文水資源勘測局,湖南 長沙 410007)

    ?

    1990—2013年湖南四水入洞庭湖汛期徑流量的變化規(guī)律

    劉培亮1,毛德華,周 慧1,2,鄧曉磊1

    (1.湖南師范大學(xué)資源與環(huán)境科學(xué)學(xué)院,湖南 長沙 410081; 2.湖南省水文水資源勘測局,湖南 長沙 410007)

    根據(jù)湖南四水(湘江、資水、沅江、澧水)主要代表水文站1990—2013年汛期月徑流量實(shí)測數(shù)據(jù),應(yīng)用滑動平均法、小波分析法和Mann-Kendall非參數(shù)趨勢檢驗(yàn)法、分配不均勻系數(shù)和變化幅度等方法和指標(biāo),分析四水流域汛期徑流量的年際和年內(nèi)變化規(guī)律,對汛期內(nèi)分配指標(biāo)與湖南省水旱災(zāi)害狀況進(jìn)行了灰色關(guān)聯(lián)度分析。結(jié)果表明:①湖南四水入洞庭湖汛期徑流量年際分配極不均勻,入湖汛期徑流量呈現(xiàn)逐漸減少的趨勢。②湘潭站和桃江站的第一主周期均為6 a,桃源站和石門站的第一主周期均為4 a。③桃江站和桃源站分別于2003年和2004年發(fā)生了突變,湘潭站和石門站均未發(fā)生過突變。④四水入湖徑流量在汛期內(nèi)各月之間分配極不均勻,但隨時間的推移,逐漸趨向均勻;汛期內(nèi)各月徑流量呈現(xiàn)波動遞減的趨勢;4站各項(xiàng)分配指標(biāo)值總體上呈現(xiàn)下降的趨勢。⑤四水入湖徑流量汛期內(nèi)分配特征指標(biāo)與流域水旱災(zāi)害的灰色關(guān)聯(lián)度極高,并在一定程度上反應(yīng)流域汛期水旱災(zāi)害情況。

    汛期徑流量;水旱災(zāi)害;小波分析;Mann-Kendall;灰色關(guān)聯(lián)度分析;湖南四水

    河川徑流量的變化規(guī)律是近幾十年水文水資源學(xué)研究的熱點(diǎn)[1],也是河流水文過程和區(qū)域水資源可持續(xù)利用研究的關(guān)鍵問題[2]。在氣候變化以及人類活動的驅(qū)動下,河川徑流量會發(fā)生年內(nèi)和年際分配的變化,這種變化必然給水資源管理利用、農(nóng)業(yè)以及生態(tài)系統(tǒng)帶來一系列的影響[3-4]。特別是近幾十年來,全球氣溫不斷升高,中國各種極端氣候頻繁出現(xiàn),發(fā)生強(qiáng)降水、洪澇和干旱事件的頻率增加[5],由此帶來的河川徑流量的變化,對經(jīng)濟(jì)和社會的可持續(xù)發(fā)展產(chǎn)生了嚴(yán)重影響。

    湖南四水(湘江、資水、沅江、澧水)作為湖南省主要河流和長江重要支流,其流域覆蓋湖南省大部分地區(qū),徑流量占洞庭湖入湖徑流量的54.6%。近幾十年來,由于自然因素和人類活動的影響,湖南省水旱災(zāi)害發(fā)生的頻率呈加劇趨勢[6],因此研究四水流域徑流量的變化規(guī)律對整個湖南省和洞庭湖流域的區(qū)域經(jīng)濟(jì)和社會的可持續(xù)發(fā)展有著重要意義。國內(nèi)學(xué)者也紛紛對湖南四水流域或洞庭湖出入徑流的規(guī)律做了大量的探討[7-10],但他們大多是對湖南四水流域或洞庭湖出入徑流的年內(nèi)分配或年際變化規(guī)律進(jìn)行分析,專門針對汛期徑流量變化規(guī)律方面的研究很少。筆者認(rèn)為南方河川徑流量汛期變化特征能夠很好地反映流域內(nèi)的水旱災(zāi)害情況。

    湖南省受季風(fēng)環(huán)流和地形因素的影響,年降水集中在4—9月(即汛期),降水量約占全年降水總量的69%,并且以局地性或全省性暴雨形式出現(xiàn),促使山嶺地區(qū)發(fā)生山洪災(zāi)害,并在洞庭湖平原區(qū)匯集成巨大的災(zāi)害性洪水[11]。湖南省農(nóng)業(yè)地域類型屬于水稻種植業(yè)類型,7—9月正值農(nóng)作物需水高峰期,平均氣溫均在35℃以上,而降水量一般只占全年降水量的20%,極易產(chǎn)生夏秋連旱的現(xiàn)象。特別是20世紀(jì)80年代以來,年年水旱災(zāi)害交替演進(jìn)的趨勢日趨明顯[12]。根據(jù)相關(guān)研究成果[13-15],20世紀(jì)90年代,長江流域來水整體偏豐,進(jìn)入21世紀(jì)以來,長江流域來水整體偏枯,而近十多年的來水可能是水文系列“豐枯交替”中的枯水系列,因此,本研究以影響流域水旱災(zāi)害的汛期徑流量為研究對象,采用長江流域河流近期一次“豐枯交替”的時間序列,這對湖南省水資源的開發(fā)和管理、經(jīng)濟(jì)和社會的可持續(xù)發(fā)展以及水旱災(zāi)害的防治有重要的意義。

    利用湖南四水主要水文控制站湘潭站(湘江)、桃江站(資水)、桃源站(沅江)和石門站(澧水)4個水文站1990—2013年24 a汛期(4—9月)月徑流量實(shí)測資料,應(yīng)用5年滑動平均法、小波分析法和Mann-Kendall非參數(shù)趨勢檢驗(yàn)法,分析四水汛期徑流量的年際變化規(guī)律,同時采用分配不均勻系數(shù)和變化幅度等指標(biāo)分析四水汛期徑流量的年內(nèi)變化規(guī)律,并對這些指標(biāo)與湖南省水旱災(zāi)害狀況進(jìn)行了灰色關(guān)聯(lián)度分析,以期認(rèn)識湖南四水入湖汛期徑流量變化規(guī)律,為流域水資源開發(fā)利用和水旱災(zāi)害的防治提供依據(jù)。

    圖1 湖南四水代表性控制站位置

    1 數(shù)據(jù)來源和研究方法

    1.1 數(shù)據(jù)來源

    本研究以湘潭水文站、桃江水文站、桃源水文站、石門水文站1990—2013年的實(shí)測數(shù)據(jù)作為湖南四水入洞庭湖汛期徑流量變化規(guī)律分析的基礎(chǔ)(圖1),其中長時間序列汛期徑流量數(shù)據(jù)和月徑流量數(shù)據(jù)均來源于湖南省水文水資源勘測局,湖南省水旱災(zāi)害情況來源于各年《湖南統(tǒng)計年鑒》。

    1.2 研究方法

    1.2.1 徑流年際變化規(guī)律研究方法

    a.滑動平均法。當(dāng)時間序列的數(shù)值由于受周期變動和隨機(jī)波動的影響,起伏較大,不易顯示出事件的發(fā)展趨勢時,使用滑動平均法可以消除這些因素的影響,顯示出事件的發(fā)展方向與趨勢(即趨勢線)。具體計算方法參照文獻(xiàn)[16]。

    b.小波分析。小波分析是一種時頻窗均可改變、表征信號特征能力的時頻局域化分析方法,借助小波變換的時頻局部化優(yōu)勢可在多時間尺度下尋找周期,具體介紹和方法計算參考文獻(xiàn)[17-18]。

    圖2 四水流域汛期徑流量變化趨勢

    c.Mann-Kendall(M-K)非參數(shù)趨勢檢驗(yàn)法。 M-K非參數(shù)趨勢檢驗(yàn)法很少受異常值干擾,也不受數(shù)據(jù)分布特征的影響,只要滿足時間序列隨機(jī)獨(dú)立即可,因此被廣泛應(yīng)用于水文和氣象時間序列的變化趨勢分析,具體計算參照文獻(xiàn)[19-20]。

    1.2.2 年內(nèi)分配變化規(guī)律研究方法

    a.不均勻性。采用不均勻系數(shù)(Cv)和完全調(diào)節(jié)系數(shù)(Cr)來衡量四水汛期內(nèi)各月徑流量分配變化規(guī)律。具體計算方法參照文獻(xiàn)[21],但文獻(xiàn)[21]中分析的是全年徑流量年內(nèi)變化特征,采用12個月份的資料,而筆者分析汛期徑流量,采用4—9月共計6個月的資料,因此對公式進(jìn)行了修改。

    徑流量汛期內(nèi)分配不均勻系數(shù)(Cv)計算公式為

    (1)

    徑流量汛期內(nèi)分配完全調(diào)節(jié)系數(shù)(Cr)計算公式為

    (2)

    從公式中可以看出,Cv和Cr值越大,汛期內(nèi)各月徑流量相差越懸殊,汛期內(nèi)徑流量分配越不均勻。

    b.變化幅度。采用相對變化幅度和絕對變化幅度兩個指標(biāo)來衡量四水汛期徑流量變化幅度,具體計算方法參照文獻(xiàn)[22]。

    1.2.3 徑流量汛期內(nèi)分配特征與水旱災(zāi)害關(guān)系的分析

    灰色關(guān)聯(lián)度分析法主要是用灰色系統(tǒng)模型對系統(tǒng)發(fā)展態(tài)勢的相互關(guān)系進(jìn)行分析,實(shí)質(zhì)上是將無限空間問題利用有限數(shù)列取代,其基本思路是根據(jù)曲線間的相似程度來判斷其關(guān)聯(lián)程度,具體計算方法參照文獻(xiàn)[23-24]。

    2 結(jié)果分析和討論

    2.1 四水流域汛期徑流量年際變化規(guī)律

    2.1.1 年際趨勢分析

    為分析湖南四水入洞庭湖汛期徑流量年際變化趨勢,分別對流域代表站湘潭站、桃江站、桃源站、石門站的1990—2013年汛期徑流量進(jìn)行了線性趨勢分析和5年的滑動平均處理(圖2)。

    從圖2可以看出:線性趨勢線和5年滑動平均曲線均明顯顯示出4站(湘潭站、桃江站、桃源站、石門站)汛期徑流量呈現(xiàn)下降的趨勢,說明湖南四水入洞庭湖汛期徑流量逐漸減少。這主要受近些年四水流域興建的大量水利工程和流域水土保持工程等對汛期徑流量的調(diào)節(jié)所致[7]。從趨勢線方程看出,桃源站遞減速率最快(7.006 1億m3/a),湘潭站次之(4.473 6億m3/a),桃江站緊隨其后(2.957 4億m3/a),石門站最慢(1.346 0億m3/a)。湘潭站和桃江站汛期徑流量隨時間波動變化的趨勢基本一致,桃源站和石門站變化趨勢一致。前兩站相對后兩站的波動頻率和幅度較小,滑動曲線多成“U”字型,后兩站多成“V”字型。湘潭站和桃源站分別在2010年和2004年前后汛期徑流量發(fā)生了較大幅度的變化;湘潭站在2010年前后幾年汛期徑流量均處于平均值之下,而在2010年迅速超過了平均值;桃源站在2004年之后汛期徑流量基本處于平均值之下。這種變化特征主要是受湖南省地形地勢的影響所致。湖南省東、西、南3面山嶺相連,中部多為丘陵和盆地,起伏和緩;北部為平原,地勢低平。從總體上看,全省恰似一個3面高、中部低、朝北部敞開的“馬蹄形”盆地。湘江和資水多流經(jīng)丘陵平原區(qū),而沅江和澧水多流經(jīng)山區(qū),并且沅江和澧水流域位于夏季風(fēng)的迎風(fēng)坡,夏季降水的幅度變化大,導(dǎo)致了這種變化特征的形成。

    4站最大汛期徑流量分別出現(xiàn)在1994年、1994年、2002年和1998年,最小汛期徑流量前3站均出現(xiàn)在2011年,石門站出現(xiàn)在2006年;前3站最大汛期徑流量和最小汛期徑流量之間均相差3~4倍,石門站甚至相差8倍,說明1990—2013年湖南四水入湖汛期徑流量的年際變化幅度較大,其年際分配極不均勻。4站汛期徑流量變化曲線基本位于多年平均徑流量之上與之下的變化時間分別是2007年、2003年、2005年和2004年,說明在這個時間點(diǎn)之前四水流域處于豐水期,在這時間之后處于枯水期,而目前正處于枯水期。

    2.1.2 小波周期分析

    為了分析湖南四水入洞庭湖汛期徑流量的周期變化,分別對流域代表站湘潭站、桃江站、桃源站、石門站的1990—2013年汛期徑流量進(jìn)行了小波周期分析(圖3,圖4)。

    從圖3和圖4可以看出:4站汛期徑流量在演化過程中存在多時間尺度特征。對湘潭站來說,在徑流演化過程中存在著3 a、5~6 a、10~16 a和22~24 a的4類尺度的周期變化規(guī)律。其中,在5~6 a時間尺度上出現(xiàn)了5次明顯的豐枯交替振蕩,在2010年處有個明顯的豐水周期;在10~16 a時間尺度上出現(xiàn)了2次振蕩;22~24 a時間尺度上出現(xiàn)了1次振蕩。同時可以看出,5~6 a和10~16 a的豐枯交替較為明顯;10~16 a和22~25 a時間尺度的周期變化具有全局性。從圖4(a)看出:湘潭站在周期上表現(xiàn)出4個很明顯的峰值,依次對應(yīng)著3 a、6 a、15 a和24 a的時間尺度。其中,最大峰值對應(yīng)著6 a的時間尺度,說明6 a左右的周期振蕩最強(qiáng),為湘江汛期徑流變化的第一主周期。

    圖3 4站小波系數(shù)等值線圖

    圖4 4站小波方差圖

    對桃江站來說,在徑流演化過程中存在著2~3 a、6~8 a、12~16 a和20~24 a的4類尺度的周期變化規(guī)律。其中,在6~8 a尺度上出現(xiàn)了豐枯交替的6次振蕩,在1999年處出現(xiàn)了明顯的豐水周期,前后也有很明顯的枯水周期;在12~16 a時間尺度上出現(xiàn)了2次振蕩;在20~24 a時間尺度上出現(xiàn)了1次振蕩。6~8 a和22~24 a的豐枯交替較為明顯,12~16 a和20~24 a時間尺度具有全局性。從圖4(b)看出:桃江站在周期上表現(xiàn)出4個很明顯的峰值,依次對應(yīng)著3 a、6 a、14 a和24 a的時間尺度,其中,最大峰值對應(yīng)著6 a的時間尺度,說明6 a左右的周期振蕩最強(qiáng),為資水汛期徑流變化的第一主周期。

    對桃源站來說,在徑流演化過程中存在著4~16 a和22~24 a的2類尺度的周期變化規(guī)律。其中,在4~16 a尺度上出現(xiàn)了豐枯交替的5次振蕩,在2004年處出現(xiàn)了明顯的豐水周期,前后也有很明顯的枯水周期;在22~24 a時間尺度上出現(xiàn)了1次振蕩。4~16 a的豐枯交替較為明顯,2類時間尺度都具有全局性。從圖4(c)看出:桃江站在周期上表現(xiàn)出4個很明顯的峰值,依次對應(yīng)著4 a、8 a、15 a和24 a的時間尺度,其中,最大峰值對應(yīng)著4 a的時間尺度,說明4 a左右的周期振蕩最強(qiáng),為沅江汛期徑流變化的第一主周期。

    對石門站來說,在徑流演化過程中存在著4~7 a、7~12 a和20~24 a的3類尺度的周期變化規(guī)律。其中,在4~7 a尺度上出現(xiàn)了豐枯交替的2次振蕩;在7~12 a時間尺度上出現(xiàn)了4次振蕩,在2003年處出現(xiàn)了明顯的豐水周期,前后也有很明顯的枯水周期;在20~24 a時間尺度上出現(xiàn)了1次振蕩。7~12 a的豐枯交替較為明顯,20~24 a時間尺度具有全局性。從圖4(d)看出:石門站在周期上表現(xiàn)出3個很明顯的峰值,依次對應(yīng)著4 a、9 a和24 a的時間尺度,其中,最大峰值對應(yīng)著4 a的時間尺度,說明4 a左右的周期振蕩最強(qiáng),為澧水汛期徑流變化的第一主周期。

    湖南四水這些變化周期主要是受集水流域降水的影響,出現(xiàn)11 a和22 a左右的周期變化主要是與本文選取的資料恰好是流域一個“豐—枯”周期有關(guān),即與太陽黑子活動的準(zhǔn)11 a和準(zhǔn)22 a周期相吻合;出現(xiàn)6 a左右的周期則與ENSO事件的5~7 a的準(zhǔn)周期相聯(lián)系;出現(xiàn)2~3 a周期振蕩與平流層大氣環(huán)流的準(zhǔn)2 a周期振蕩(QBO)存在一致性關(guān)系[25]。

    2.1.3 突變點(diǎn)分析

    為進(jìn)一步分析湖南四水入洞庭湖汛期徑流量年際變化規(guī)律,分別對流域代表站湘江湘潭站、資水桃江站、沅江桃源站和澧水石門站汛期徑流量進(jìn)行了Mann-Kendall趨勢檢驗(yàn)和突變分析(圖5)。

    圖5 四水流域汛期徑流量M-K檢測統(tǒng)計值

    根據(jù)M-K檢驗(yàn)結(jié)果,湘潭站在2003年之前、桃江站在2000年之前、桃源站在2004年之前和石門站在1995~1999年之間的UF統(tǒng)計量為正值,即這些年份各站汛期徑流量呈現(xiàn)上升趨勢,但均在95%的置信區(qū)間內(nèi),所以增加的趨勢不明顯。4站其他年份UF統(tǒng)計量為負(fù)值,表明這些年份呈現(xiàn)下降的變化趨勢。只有桃江站在2007—2013年和桃源站在2011—2013年突破了95%負(fù)值置信區(qū)間,該時間內(nèi)汛期年徑流量呈現(xiàn)顯著下降趨勢。其他兩站UF統(tǒng)計量均處于95%置信區(qū)間內(nèi)變化,說明該時間內(nèi)汛期年徑流量上升或下降的趨勢不明顯。湘潭站的UF和UB曲線在2003年發(fā)生了相交,但在其后幾年內(nèi)UF曲線均未突破95%置信區(qū)間,說明湘江入湖汛期徑流量未發(fā)生突變。桃江站的UF和UB曲線在2003年發(fā)生了相交,UF曲線隨后在2007年突破了95%置信區(qū)間,說明資水入湖汛期徑流量在2003年發(fā)生了突變。桃源站的UF和UB曲線在2004年發(fā)生了相交,UF曲線隨后在2011年突破了95%置信區(qū)間,說明沅江入湖汛期徑流量在2004年發(fā)生了突變。石門站的UF和UB曲線先后發(fā)生了4次相交,但在其后幾年內(nèi)UF曲線均未突破95%置信區(qū)間,說明,澧水入湖汛期徑流量在1990—2013年間沒有發(fā)生過突變。

    從20世紀(jì)90年代以來,四水流域水利水保工程建設(shè)逐漸增多,人類活動對徑流量的干擾的作用逐漸積累加強(qiáng),加之氣候變化,降水減少這另一主導(dǎo)因素的共同作用,資水在2003年發(fā)了突變,沅江在2004年發(fā)生了突變。湘江由于支流多(2 157條),流域面積廣(94 660 km2),并實(shí)行了48個梯級開發(fā),所以入湖汛期徑流量未發(fā)生突變。澧水由于流經(jīng)山區(qū)和20余座大中型水庫的調(diào)節(jié),也未發(fā)生突變。

    2.2 四水流域汛期徑流量年內(nèi)變化規(guī)律

    2.2.1 汛期徑流量月變化特征

    湖南四水屬于亞熱帶季風(fēng)氣候區(qū)河流,夏季降水多,冬季降水少,河流以降水補(bǔ)給為主。受降水條件、蒸發(fā)、流域自然地理特征和人類活動等的影響,四水流域入湖徑流量在汛期內(nèi)各月之間分配極不均勻[7,11]。從湘潭站、桃江站、桃源站和石門站4站汛期內(nèi)各月徑流量變化曲線圖(圖6)可以看出:

    圖6 四水流域徑流量汛期內(nèi)各月變化趨勢

    總體上,4站各月徑流量呈現(xiàn)波動中遞減的趨勢,其中7月和8月遞減的趨勢比較明顯。各月徑流量的波動變化的幅度不一,其中6月、7月、8月波動幅度較大,特別是7月的最小徑流量和最大徑流量相差最懸殊,年際分配最不均勻;4月徑流量變化曲線與4—9月平均徑流量曲線基本吻合,9月徑流量變化曲線基本位于4—9月平均徑流量曲線之下,可見4月和9月波動幅度較小,年際分配最為均勻,5月次之。24年間,湘潭站月最大徑流量出現(xiàn)在2010年的6月,桃江站最大徑流量出現(xiàn)在1998年6月,桃源站最大徑流量出現(xiàn)在1996年的7月,石門站最大徑流量出現(xiàn)在1991年7月。

    從表1統(tǒng)計情況可以看出,4站4—7月徑流量之和均占汛期徑流量之和的78%以上,說明4—7月是湖南四水入洞庭湖的主汛期。其中湘潭站4—5月所占比重最大,但與6—7月的相差不大,為39.96%,其他3站均為6—7月所占比重最大,分別達(dá)到42.95%、48.47%、53.17%。從4站多年平均徑流量來看,近24年湖南四水入洞庭湖年均汛期徑流量的總量為1 186.88億m3,其中湘江占總量的39.31%,資水為13.46%,沅江為39.77%,澧水為7.46%。

    表1 四水流域多年平均汛期徑流量分配

    2.2.2 徑流量的汛期內(nèi)分配特點(diǎn)

    為進(jìn)一步分析1990—2013年湖南四水入洞庭湖的汛期徑流量年內(nèi)變化特征,對流域代表站湘潭站、桃江站、桃源站和石門站的不均勻系數(shù)、完全調(diào)節(jié)系數(shù)、相對變化幅度和絕對變化幅度等指標(biāo)進(jìn)行統(tǒng)計,見圖7和表2。

    圖7 四水流域流量汛期內(nèi)分配特征指標(biāo)值

    表2 四水流域徑流量汛期內(nèi)分配特征多年平均指標(biāo)值

    水文站不均勻系數(shù)完全調(diào)節(jié)系數(shù)相對變化幅度絕對變化幅度/億m3湘潭0.520.204.96101.45桃江0.500.194.8135.38桃源0.590.226.05122.62石門0.770.2917.7431.19

    從圖7中不均勻系數(shù)和完全調(diào)節(jié)系數(shù)可以看出:4站各自的不均勻系數(shù)和完全調(diào)節(jié)系數(shù)變化趨勢基本相似,湘潭站、桃江站和桃源站之間的不均勻系數(shù)和完全調(diào)節(jié)系數(shù)變化規(guī)律也基本相似,石門站與其他3站相差很大。湘潭站、桃江站和桃源站的不均勻系數(shù)和完全調(diào)節(jié)系數(shù)在1991—1994年、1999—2002年和2006—2009年數(shù)值相對較小,而石門站只是在1994年、2001年和2012年數(shù)值較小,即這些年份徑流量在汛期各月之間的分配較為均勻。從圖7中相對變化幅度和絕對變化幅度可以看出:湘潭站、桃江站和桃源站相對變化幅度值大多在3~10之間變化,波動幅度較為平穩(wěn);石門站有些年份值很小,有些年份相當(dāng)大,甚至超過了50,波動幅度極大。由于四水各自入湖流量的巨大差異,4站絕對變化幅度差別較大,但從波動趨勢來看,湘潭站和桃江站較為相似,沅江和澧水較為相似。

    對比圖7中的各個指標(biāo)可以看出:4站各項(xiàng)指標(biāo)波動變化幅度均較大,變化曲線均由多個“V”字組成,說明四水入洞庭湖徑流量汛期內(nèi)分配有些年份極不均勻;4站各項(xiàng)指標(biāo)總體呈現(xiàn)出下降的趨勢,說明四水入洞庭湖徑流量汛期內(nèi)分配隨時間增加逐漸向均勻化轉(zhuǎn)變,這主要是由于隨著社會經(jīng)濟(jì)的發(fā)展要求,加強(qiáng)了對流域水資源的管理,以及水利工程的建設(shè)等,都使得流域徑流汛期內(nèi)分配趨向于均勻化。

    從表2多年平均值和圖6變化曲線來看,不均勻系數(shù)、完全調(diào)節(jié)系數(shù)和相對變化幅度,石門站數(shù)值最高,桃源站次之,湘潭站隨后,桃江站最低,說明石門站的汛期內(nèi)徑流量分配最不均勻,湘潭站和桃江站的最較為均勻。受澧水地形和地勢條件的影響,經(jīng)常產(chǎn)生氣旋雨和鋒面雨,且降雨強(qiáng)度大,在時空分布上也不均勻,加上澧水河流較短,上中游河道坡度陡,洪水暴漲暴落,極易發(fā)生山洪,所以澧水汛期徑流量的汛期內(nèi)分配極不均勻。其他3條河流由于河流較長,流域范圍較大,干流上多較大型水庫和水電站,所以徑流量汛期內(nèi)分配比石門站較為均勻。

    2.2.3 徑流量汛期內(nèi)分配特征與水旱災(zāi)害的關(guān)系分析

    為分析湖南四水入湖汛期徑流量變化特征與流域的水旱災(zāi)害的關(guān)系,采用湘潭站、桃江站、桃源站和石門站的不均勻系數(shù)、完全調(diào)節(jié)系數(shù)、相對變化幅度和絕對變化幅度等徑流量汛期內(nèi)分配指標(biāo)值與湖南省受水災(zāi)和旱災(zāi)面積進(jìn)行灰色關(guān)聯(lián)度分析,分析結(jié)果見表3。

    表3 灰色關(guān)聯(lián)度統(tǒng)計結(jié)果

    從表3中可以看出流域的水災(zāi)、旱災(zāi)、水旱災(zāi)害之和狀況與4站的各項(xiàng)分配特征指標(biāo)的灰色關(guān)聯(lián)度分別均高于0.78、0.81、0.84,說明湖南四水入湖徑流量汛期內(nèi)分配特征指標(biāo)值與流域水旱災(zāi)害的面積關(guān)聯(lián)度極高,變化曲線相似,其中旱災(zāi)比水災(zāi)的關(guān)聯(lián)度高,水旱災(zāi)害之和又比前兩者高,這與南方地區(qū)經(jīng)常發(fā)生的干旱為伏旱有一定的關(guān)系。4站各項(xiàng)指標(biāo)與水旱災(zāi)害的關(guān)聯(lián)度之間的差別大都在0~1.5,說明4者指標(biāo)值的大小一定程度上均能反映流域的水旱災(zāi)害狀況。

    從圖7中得到1994年和2012年是湖南四水入洞庭湖徑流量汛期內(nèi)分配較為均勻的年份(4站4項(xiàng)分配指標(biāo)值均較小),這兩年湖南省受水旱災(zāi)害的面積較小,1994年受水旱災(zāi)面積153.393萬hm2,2012年受水旱災(zāi)面積106.462萬hm2。四水入洞庭湖徑流量汛期內(nèi)分配極不均勻的年份(4站4項(xiàng)分配指標(biāo)值均較大)為:1990年、1996年、1998年、2003年、2005年、2010年。經(jīng)查各年《湖南統(tǒng)計年鑒》,1996年、1998年和2010年湖南省出現(xiàn)了較大的洪澇災(zāi)害,1996年受水災(zāi)面積123.666萬hm2,1998年受水災(zāi)面積213.41萬hm2,2010年受水災(zāi)面積201.133萬hm2;1990年和2003年湖南省出現(xiàn)了較為嚴(yán)重的干旱災(zāi)害,1990年受旱災(zāi)面積148.533萬hm2,2003年受旱災(zāi)面積189.508萬hm2;2005年《湖南統(tǒng)計年鑒》只有水災(zāi)受災(zāi)面積數(shù)據(jù),為59.413萬hm2,無旱災(zāi)受災(zāi)面積數(shù)據(jù),但是4站7—9月的徑流量均不足正常年份的1/2,說明湖南省發(fā)生了嚴(yán)重的旱情。因此,湖南四水入洞庭湖徑流量汛期內(nèi)分配特征能較好地反映當(dāng)年湖南省水旱災(zāi)害情況,如若反映徑流量汛期內(nèi)分配不均勻性的4項(xiàng)指標(biāo)出現(xiàn)異常偏大值情況,說明當(dāng)年發(fā)生水旱災(zāi)害的可能性大。由于缺少湖南省具體的4—9月水旱災(zāi)害情況資料,此結(jié)論有待進(jìn)一步研究。

    3 結(jié) 論

    a.1990—2013年湖南四水入洞庭湖汛期徑流量年際分配極不均勻,4站汛期徑流量變化曲線呈現(xiàn)下降的趨勢,湖南四水入湖汛期徑流量逐漸減少。湘潭站和桃江站汛期徑流量的隨時間的波動變化趨勢基本一致,而桃源站和石門站的變化趨勢一致。

    b.小波周期分析顯示4站汛期徑流量在演化過程中存在多時間尺度特征。湘潭站的變化周期為3 a、5~6 a、10~16 a和22~24 a,第一主周期為6 a;桃江站的變化周期為2~3 a、6~8 a、12~16 a和20~24 a,第一主周期為6 a;桃源站的變化周期為4~16 a和22~24 a,第一主周期為4 a;石門站的變化周期為4~7 a、7~12 a和20~24 a,第一主周期為4 a。

    c.桃江站在2007—2013年和桃源站在2011—2013年年突破了95%負(fù)值置信區(qū)間,該時間內(nèi)汛期年徑流量呈現(xiàn)顯著下降趨勢。其他兩站UF統(tǒng)計量均處于95%置信區(qū)間內(nèi)變化,該時間內(nèi)汛期年徑流量上升或下降的趨勢不明顯。桃江站和桃源站分別于2003年和2004年發(fā)生了突變。湘潭站和石門站沒有發(fā)生過突變。

    d.1990—2013年湖南四水入湖徑流量汛期內(nèi)分布極不均勻。4—7月是湖南四水入洞庭湖的主汛期;總體上,4站各月徑流量呈現(xiàn)波動中逐年遞減的趨勢。24 a湖南四水入洞庭湖年平均汛期徑流總量為1186.88億m3,其中湘江占總量的39.31%,資水為13.46%,沅江為39.77%,澧水為7.46%。

    e.4站各自的不均勻系數(shù)和完全調(diào)節(jié)系數(shù)變化趨勢基本相似,湘潭站、桃江站和桃源站在1991—1994年、1999—2002年和2006—2009年的徑流量汛期分配較為均勻,而石門站只是在1994年、2001年和2012年的分配較為均勻。湘潭站、桃江站和桃源站相對變化幅度值大都在3~10之間變化,波動幅度較為平穩(wěn);石門站有時超過50,波動幅度極大;石門站汛期內(nèi)分配最不均勻。4站各項(xiàng)指標(biāo)總體呈現(xiàn)下降的趨勢,說明四水入洞庭湖徑流量汛期內(nèi)分配隨時間增加逐漸向均勻化轉(zhuǎn)變。

    f.湖南四水入湖徑流量汛期內(nèi)分配特征指標(biāo)與湖南省水旱災(zāi)害的灰色關(guān)聯(lián)度極高。汛期內(nèi)分配特征指標(biāo)能較好反映當(dāng)年流域汛期水旱災(zāi)害情況,如若四水入湖徑流量汛期內(nèi)分配不均勻性的4項(xiàng)指標(biāo)出現(xiàn)異常偏大,則當(dāng)年發(fā)生水旱災(zāi)害的可能性較大。

    [1] 王浩,嚴(yán)登華,賈仰文,等.現(xiàn)代水文水資源學(xué)科體系及研究前沿和熱點(diǎn)問題[J].水科學(xué)進(jìn)展,2010,21(4): 479-489.(WANG Hao,YAN Denghua,JIA Yangwen,et al.Subject system of modern hydrology and water resources and research frontiers and hot issues[J].Advances in Water Science,2010,21(4): 479-489.(in Chinese))

    [2] 胡四一.現(xiàn)代水文學(xué)及水資源的學(xué)科前沿和研究需求[J].水文,2006,26(3):1-4.(HU Siyi.The academic foreland and research demand of hydrology and water resources[J].Journal of China Hydrology,2006,26(3):1-4.(in Chinese))

    [3] 鄢波,夏自強(qiáng),周艷先,等.黑龍江哈巴羅夫斯克站徑流變化規(guī)律[J].水資源保護(hù),2013,29(3):29-33.((YAN Bo,XIA Ziqiang,ZHOU Yanxian,et al.Variation of runoff at Khabarovsk Station on Heilongjiang River[J].Water Resources Protection,2011,27(4):48-54.(in Chinese))

    [4] 王浩,賈仰文,王建華,等.人類活動影響下的黃河流域水資源演化規(guī)律初探[J].自然資源學(xué)報,2005,20(2):157-162.(WANG Hao,JIA Yangwen,WANG jianhua,et al.Evolutionary laws of Yellow River Basin’s water resources under the impact of human activities[J].Journal of Natural Resources,2005,20(2):157-162.(in Chinese))

    [5] 任國玉,姜彤,李維京,等.氣候變化對中國水資源情勢影響綜合分析[J].水科學(xué)進(jìn)展,2008,19(6):772-779.(REN Guiyu,JIANG Tong,LI Weijing,et al.An integrated assessment of climate change impacts on China’s water resources[J].Advances in Water Science,2008,19(6):772-779.(in Chinese))

    [6] 謝煉,朱毅.湖南水旱災(zāi)害頻繁成因淺析[J].湖南水利水電,2005(5):11-12.(XIE Lian,ZHU Yi.Frequent flood and drought disasters’ cause analyses of Hunan Province[J].Water Resources & Hydropower of Hunan Province,2005(5):11-12.(in Chinese))

    [7] 覃紅燕,謝永宏,鄒冬生.湖南四水入洞庭湖水沙演變及成因分析[J].地理科學(xué),2012,32(5):609-615.(QIN Hongyan,XIE Yonghong,ZOU Dongsheng.Changes of runoff and sediment discharge into Dongting Lake from the four rivers in Hunan Province[J].Scientia Geographica Sinica,2012,32(5):609-615.(in Chinese))

    [8] 李 暉,尹輝,白旸,等.近60年洞庭湖區(qū)水沙演變特征及趨勢預(yù)測[J].水土保持研究,2013,20(3):139-142.(LI Hui,YIN Hui,BAI Yang,et al.Research on evolution characteristics and tendency prediction of runoff and sediment in recent 60 years in Dongting Lake[J].Research of Soil and Water Conservation,2013,20(3):139-142.(in Chinese))

    [9] 毛北平,梅軍亞,張金輝,等.洞庭湖三口洪道水沙輸移變化分析[J].人民長江,2010,41(2):38-46.(MAO Beiping,MEI Junya,ZHANG Jinhui,et al.Analysis of water and sediments transportation of three river-outlets flood channels from Yangtze River to Dongting Lake[J].Yangtze River,2010,41(2): 38-46.(in Chinese))

    [10] 胡光偉,毛德華,張旺,等.60年來湘江干流徑流泥沙過程變化及驅(qū)動力分析[J].湖南師范大學(xué)自然科學(xué)學(xué)報,2014,37(2):12-37.(HU Guangwei,MAO Dehua,ZHANG Wang,et al.Analysis on variation characteristics and driving forces of annual runoff and sediment in the Xiangjiang River over the past 60 years[J].Journal of Natural Science of Hunan Normal University,2014,37(2):12-37.(in Chinese))

    [11] 毛德華,李景保,龔重惠,等.湖南省洪澇災(zāi)害研究[M].長沙:湖南師范大學(xué)出版社,2000:96-133.

    [12] 李景保,王克林,朱寧,等.湖南省水旱災(zāi)害與暴雨徑流資源的調(diào)控[J].自然資源學(xué)報,2004,19(6):716-724.(LI Jingbao,WANG Kelin,ZHU Ning,et al.Rlood & drought disaster and regulation of storm-runoff resources in Hunan Province[J].Journal of Natural Resources,2004,19(6):716-724.(in Chinese))

    [13].戴仕寶,楊世倫.近50年來長江水資源特征變化分析[J].自然資源學(xué)報,2006,21(4):501-506.(DAI Shibao,YANG Shilun.Variations in water resources of the Yangtze River over the last five decades[J].Journal of Natural Resources,2006,21(4):501-506.(in Chinese))

    [14] 張建云,章四龍,王金星,等.近50年來中國六大流域年際徑流變化趨勢研究[J]水科學(xué)進(jìn)展,2007,18 (2): 230-234.(ZHANG Jianyun,ZHANG Silong,WANG Jinxing,et al.Study on runoff trends of the six larger basins in China over the past 50 years[J] Advances in Water Science,2007,18 (2):230-234.(in Chinese))

    [15] 胡光偉,毛德華,李正最,等.三峽工程建設(shè)對洞庭湖的影響研究綜述[J].自然災(zāi)害學(xué)報,2013,22(5): 44-52.(HU Guangwei,MAO Dehua,LI Zhengzui,et al.Research overview on the influence of the Three Gorges Project construction on Dongting Lake[J].Journal of Natural Disasters,2013,22(5): 44-52.(in Chinese))

    [16] 赫振平,趙春.桑干河流域徑流年際變化規(guī)律[J].水資源保護(hù),2011,27(3):54-57.(HE Zhenping,ZHAO Chun.Study on interannual variation of runoff in Sanggan River Basin[J].Water Resources Protection,2011,27(3):54-57.(in Chinese))

    [17] 劉宇峰,孫虎,原志華.基于小波分析的汾河河津站徑流與輸沙的多時間尺度特征[J].地理科學(xué),2012,32(6): 764-770.(LIU Yufeng,SUN Hu,YUAN Zhihua.Multitime scale features of runoff and sediment discharge on wavelet analysis at Hejin Station of Fenhe River Basin[J].Scientia Geographica Sinica,2012,32(6): 764-770.(in Chinese))

    [18] 劉二佳,張曉萍,張建軍,等.1956—2005年窟野河徑流變化及人類活動對徑流的影響分析[J].自然資源學(xué)報,2013,28(7): 1159-1168.(LIU Erjia,ZHANG Xiaoping,ZHANG Jianjun,et al.Variation of annual streamflow and the effect of human activity in the Kuye River during 1956 to 2005[J].Journal of Natural Resources,2013,28(7): 1159-1168.(in Chinese))

    [19] 葉許春,張奇,劉健,等.鄱陽湖流域天然徑流變化特征與水旱災(zāi)害[J].自然災(zāi)害學(xué)報,2012,21(1): 140-147.(YE Xuchun,ZHANG Qi,LIU Jian,et al.Natural runoff change characteristics and flood/drought disasters in Poyang Lake catchment basin[J].Journal of Natural Disasters,2012,21(1): 140-147.(in Chinese))

    [20] 馬聰,魯春霞,李亦秋,等.安康水庫下游徑流演變及水庫運(yùn)行對徑流影響分析[J].長江流域資源與環(huán)境,2013,22(11):1433-1440.(MA Cong,LU Chunxia,LI Yiqiu,et al.Runoff change in the lower reaches of Ankang Reservoir and the influence of Ankang Reservoir on its downstream[J].Resources and Environment in the Yangtze Basin,2013,22(11):1433-1440.(in Chinese))

    [21] 郭巧玲,楊云松,暢祥生,等.1957—2008年黑河流域徑流年內(nèi)分配變化[J].地理科學(xué)科展,2011,30(5): 550-556.(GUO Qiaoling,YANG Yunsong,CHANG Xiangsheng,et al.Annual variation of Heihe River runoff during 1957—2008[J].Progress in Geography,2011,30(5): 550-556.(in Chinese))

    [22] 陸志華,夏自強(qiáng),于嵐嵐,等.松花江干流中游段徑流年內(nèi)分配變化規(guī)律[J].河海大學(xué)學(xué)報:自然科學(xué)版,2012,40 (1):63-69.(LU Zhihua,XIA Ziqiang,YU Lanlan et al.Changes of annual runoff distribution in middle reaches of Songhua River[J].Journal of Hohai University:Natural Sciences,2012,40 (1):63-69.(in Chinese))

    [23] 李雯,趙寶峰,王杰,等.基于灰色關(guān)聯(lián)度法的涇河徑流量變化因素分析[J].水資源與水工程學(xué)報,2009,20(2):92-94.(LI Wen,ZHAO Baofeng,WANG Jie,et al.Analysis on the change factors of runoff in Jinghe River based on grey correlation degree[J].Journal of Water Resources & Water Engineering,2009,20(2):92-94.(in Chinese))

    [24] 朱金峰,梁忠民,湯曉芳.湖南省農(nóng)村水資源保護(hù)現(xiàn)狀評價與趨勢預(yù)測[J].南水北調(diào)與水利科技,2013,11(3): 6-11.(ZHU Jinfeng,LIANG Zhongming,TANG Xiaofang.Evauation of current status and trend prediction of rural water resources protection in Hunan Province[J].South-to-North Water Transfers and Water Science & Technology,2013,11(3): 6-11.(in Chinese))

    [25] 張劍明,章新平,黎祖賢,等.湖南省 46年來降水的氣候特征[J].熱帶氣象學(xué)報,2008,24(5):512-518.(ZHANG Jianming,ZHANG Xinping,LI Zuxian,et al.Climatological features of annual precipitation in Hunan Province during 46 years[J].Journal of Tropical Meteorology,2008,24(5):512-518.(in Chinese))

    ·簡訊·

    第五屆全國水工抗震防災(zāi)學(xué)術(shù)交流會將在湖北宜昌召開

    2015年11月6-8日,第五屆全國水工抗震防災(zāi)學(xué)術(shù)會議將在湖北省宜昌市舉行。會議主辦單位:中國水力發(fā)電工程學(xué)會抗震防災(zāi)專業(yè)委員會;會議承辦單位:三峽大學(xué),水資源高效利用與工程安全國家工程研究中心,湖北省水力發(fā)電工程學(xué)會;會議協(xié)辦單位:河海大學(xué),清華大學(xué),武漢大學(xué)等。

    我國是一個多地震國家,目前全國已建或在建的高于百米的大壩數(shù)量已經(jīng)一百多座,一批在建和擬建的200~300 m級戰(zhàn)略性水電工程,規(guī)模巨大、地震烈度高,抗震防災(zāi)成為工程建設(shè)中的關(guān)鍵??偨Y(jié)水電工程抗震的關(guān)鍵技術(shù)和經(jīng)驗(yàn)教訓(xùn)具有重要意義。本次水工抗震防災(zāi)學(xué)術(shù)會議將就水工建筑物抗震設(shè)計、分析理論,地震震后水工建筑物震損調(diào)查、震害分析、災(zāi)害評估及修復(fù)加固,水工結(jié)構(gòu)現(xiàn)場檢測與監(jiān)測、健康診斷與鑒定,水庫引發(fā)地震、地震預(yù)測和預(yù)警,水電工程抗震措施、抗震安全評價體系,水工結(jié)構(gòu)工程振動、爆炸與沖擊等關(guān)鍵科學(xué)技術(shù)問題,進(jìn)行廣泛的交流探討。

    (本刊編輯部供稿)

    Variation law of runoff in flood seasons into Dongting Lake from Four Rivers in Hunan Province during 1990—2013

    LIU Peiliang1, MAO Dehua1, ZHOU Hui1,2, DENG Xiaolei1

    (1.CollegeofResourcesandEnvironmentalScience,HunanNormalUniversity,Changsha410081,China;2.HydrologyandWaterResourcesSurveyBureauofHunanProvince,Changsha410007,China)

    According to the data of runoff of the Four Rivers(Xiangjiang River, Zishui River, Yuanjiang River, Lishui River) measured by main hydrological station during 1990—2013, we analyze the annual variation law of runoff in flood season of the Four River and conduct a gray correlation degree analysis over the distribution of indicators in the flood season and floods and droughts situation in Hunan Province with the methods of moving average, wavelet analysis, Mann-Kendall mutation test, non-uniformity coefficient and amplitude of variation, etc.The results indicate that: (1) the distribution of inter-annual flood season’s runoff into Dongting Lake from the Four Rivers was quite uneven, while the flood season’s runoff into Dongting Lake presented a trend of decreasing in general.(2) For Xiangtan station and Taojiang station, the first main cycle was 6 years; for Taoyuan station and Shimen station, the first main cycle was 4 years.(3) At Taojiang station and Taoyuan station, mutations occurred in 2003 and 2004, respectively.For Xiangtan station and Shimen station, there had been no mutation.(4) During the flood season, the distribution of runoff into Dongting Lake from the Four Rivers was quite uneven, however, the runoff distribution tended to be even over time.The monthly runoff showed a decreasing trend with fluctuation.As a whole, the indicators from the four stations presented a decreasing trend.(5) The grey correlation degree of the runoff distribution characteristics index value of Four Rivers into Dongting Lake during the flood season and river basin’s flood and drought disasters was extremely high, reflecting the situation of the basin’s flood and drought disasters.

    runoff in flood season; flood and drought disasters; wavelet analysis; Mann-Kendall; grey correlation degree analysis; Four Rivers in Hunan Province

    10.3880/j.issn.1004-6933.2015.04.010

    教育部2012年高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(20124306110009);湖南省重點(diǎn)學(xué)科建設(shè)項(xiàng)目(2011-001);湖南省自然科學(xué)基金重點(diǎn)項(xiàng)目(12JJ2026)

    劉培亮(1989—),男,碩士研究生,研究方向?yàn)樗Y源與水環(huán)境,水旱災(zāi)害。E-mail: lpl198938@163.com

    毛德華,教授。E-mail: mdh408122@sohu.com

    TV121+.4

    A

    1004-6933(2015)04-0052-10

    2014-11-27 編輯:彭桃英)

    猜你喜歡
    四水桃江湘潭
    桃江醫(yī)改值得學(xué)習(xí)的寶貴之處
    蹤影組詩
    滇池(2022年5期)2022-04-30 21:44:36
    桃江鳳凰山:屈原發(fā)出《天問》絕唱
    開封市“四水”存在的問題與解決思路
    四水同治
    湘潭是個好地方
    湘潭紅色文化軟實(shí)力的提升研究
    活力(2019年21期)2019-04-01 12:16:10
    湘潭大學(xué)藝術(shù)學(xué)院作品選
    流行色(2017年12期)2017-10-26 03:08:22
    河?xùn)|·河西
    遼河(2015年2期)2016-04-06 04:40:20
    湘潭高新區(qū)兩大特色產(chǎn)業(yè)園躋身“湖南隊(duì)”
    精品久久久久久久久亚洲| 午夜a级毛片| 精品99又大又爽又粗少妇毛片| а√天堂www在线а√下载| 五月伊人婷婷丁香| 黄色日韩在线| 久久久久久久亚洲中文字幕| 日韩 亚洲 欧美在线| av在线天堂中文字幕| 国产精品福利在线免费观看| 国产精品亚洲美女久久久| 成人亚洲精品av一区二区| 性色avwww在线观看| 久久久久久国产a免费观看| 久久国产乱子免费精品| 精品不卡国产一区二区三区| 国内精品久久久久精免费| 国产中年淑女户外野战色| 超碰av人人做人人爽久久| 成人国产麻豆网| 亚洲第一区二区三区不卡| 亚洲精华国产精华液的使用体验 | 九色成人免费人妻av| 午夜影院日韩av| 国语自产精品视频在线第100页| av在线老鸭窝| 久久久久性生活片| 99国产精品一区二区蜜桃av| 在线免费观看的www视频| 国产精品亚洲美女久久久| 国产欧美日韩一区二区精品| 超碰av人人做人人爽久久| 国产精品国产三级国产av玫瑰| 97碰自拍视频| h日本视频在线播放| 免费观看在线日韩| 婷婷亚洲欧美| 亚洲国产精品久久男人天堂| 1024手机看黄色片| 日本a在线网址| 91麻豆精品激情在线观看国产| 此物有八面人人有两片| 亚洲经典国产精华液单| 人人妻人人看人人澡| 亚洲成av人片在线播放无| 欧美+日韩+精品| 日日摸夜夜添夜夜添小说| 一本精品99久久精品77| 国产麻豆成人av免费视频| 嫩草影院精品99| 黄片wwwwww| 精品熟女少妇av免费看| 成人精品一区二区免费| 晚上一个人看的免费电影| 午夜免费激情av| 欧美性猛交╳xxx乱大交人| 久久亚洲精品不卡| 寂寞人妻少妇视频99o| 亚洲国产日韩欧美精品在线观看| 99热这里只有精品一区| 国产淫片久久久久久久久| 日韩欧美免费精品| 久久精品国产清高在天天线| 欧美日韩精品成人综合77777| 女同久久另类99精品国产91| 久久久久久久久中文| 中文字幕久久专区| 国产综合懂色| 亚洲精品乱码久久久v下载方式| 欧美极品一区二区三区四区| 卡戴珊不雅视频在线播放| 久久久久久大精品| 亚洲精品粉嫩美女一区| 亚洲av成人av| 国产精品一及| 久久久午夜欧美精品| 伦精品一区二区三区| 亚洲中文日韩欧美视频| 国产伦一二天堂av在线观看| 真人做人爱边吃奶动态| 日韩大尺度精品在线看网址| 精品一区二区三区av网在线观看| 蜜桃久久精品国产亚洲av| 别揉我奶头~嗯~啊~动态视频| 亚洲婷婷狠狠爱综合网| 一级毛片久久久久久久久女| 国国产精品蜜臀av免费| 最新中文字幕久久久久| avwww免费| 久久精品国产亚洲av香蕉五月| 精品人妻视频免费看| 亚洲激情五月婷婷啪啪| 三级国产精品欧美在线观看| av天堂中文字幕网| 久久久久国产精品人妻aⅴ院| 久久鲁丝午夜福利片| 国产午夜精品论理片| 亚洲精品乱码久久久v下载方式| 一边摸一边抽搐一进一小说| 一区二区三区四区激情视频 | 成人亚洲欧美一区二区av| 国产 一区 欧美 日韩| 美女高潮的动态| 国产片特级美女逼逼视频| 亚洲aⅴ乱码一区二区在线播放| 久久久精品94久久精品| 日日撸夜夜添| 亚洲欧美日韩无卡精品| 日韩欧美 国产精品| 亚洲人成网站在线播放欧美日韩| 伊人久久精品亚洲午夜| 小蜜桃在线观看免费完整版高清| 真实男女啪啪啪动态图| 国产久久久一区二区三区| 欧美丝袜亚洲另类| 乱系列少妇在线播放| 亚洲精品国产av成人精品 | 赤兔流量卡办理| 亚洲人成网站高清观看| 国产色婷婷99| 色噜噜av男人的天堂激情| 国产三级在线视频| 老熟妇仑乱视频hdxx| 99视频精品全部免费 在线| 欧美激情久久久久久爽电影| 天天躁夜夜躁狠狠久久av| 美女免费视频网站| 日韩一本色道免费dvd| 我的老师免费观看完整版| 一级黄色大片毛片| 日韩精品中文字幕看吧| 国产伦在线观看视频一区| 在线国产一区二区在线| 中文字幕久久专区| 国产高清三级在线| 久久精品国产鲁丝片午夜精品| 老司机福利观看| 成年免费大片在线观看| 老司机影院成人| 一边摸一边抽搐一进一小说| 天天躁日日操中文字幕| 久久6这里有精品| 深爱激情五月婷婷| 国产成人福利小说| 国产精品人妻久久久影院| 国产成人91sexporn| 久久这里只有精品中国| 一区福利在线观看| 亚洲精品成人久久久久久| 97人妻精品一区二区三区麻豆| 蜜臀久久99精品久久宅男| 91午夜精品亚洲一区二区三区| 午夜福利高清视频| 国内少妇人妻偷人精品xxx网站| 1000部很黄的大片| 永久网站在线| 不卡一级毛片| 两个人的视频大全免费| 精品一区二区三区人妻视频| 你懂的网址亚洲精品在线观看 | 91在线精品国自产拍蜜月| 日本免费a在线| 少妇裸体淫交视频免费看高清| 午夜亚洲福利在线播放| 国产黄色小视频在线观看| 亚洲av一区综合| 精品免费久久久久久久清纯| 俄罗斯特黄特色一大片| 亚洲av一区综合| 毛片女人毛片| 中文字幕精品亚洲无线码一区| 成人亚洲精品av一区二区| 老女人水多毛片| 精品国内亚洲2022精品成人| 国产色婷婷99| 精品国内亚洲2022精品成人| 人妻丰满熟妇av一区二区三区| 精品久久久久久久久久免费视频| 中文在线观看免费www的网站| 亚洲精品乱码久久久v下载方式| 老司机福利观看| 亚洲av不卡在线观看| 国产精品久久久久久久久免| 女人被狂操c到高潮| 九九爱精品视频在线观看| 校园春色视频在线观看| 久久人妻av系列| 免费观看精品视频网站| 免费av观看视频| 成人亚洲欧美一区二区av| 日本免费一区二区三区高清不卡| 精华霜和精华液先用哪个| 天堂√8在线中文| 亚洲av美国av| 亚洲精品色激情综合| 中文字幕熟女人妻在线| 男女视频在线观看网站免费| 给我免费播放毛片高清在线观看| 久久午夜亚洲精品久久| 毛片一级片免费看久久久久| 夜夜爽天天搞| 国产成人aa在线观看| 极品教师在线视频| 日本在线视频免费播放| 成人综合一区亚洲| 免费观看精品视频网站| 小蜜桃在线观看免费完整版高清| 熟妇人妻久久中文字幕3abv| 大香蕉久久网| 成人特级黄色片久久久久久久| 日韩欧美 国产精品| 亚洲精品成人久久久久久| 99久久中文字幕三级久久日本| 国产精品三级大全| 国产精品女同一区二区软件| 18禁在线无遮挡免费观看视频 | 国产aⅴ精品一区二区三区波| 日韩大尺度精品在线看网址| 欧美+日韩+精品| 少妇丰满av| 精品国产三级普通话版| 男女下面进入的视频免费午夜| 国产高清激情床上av| 国产私拍福利视频在线观看| 欧美性猛交╳xxx乱大交人| 一本久久中文字幕| 国产69精品久久久久777片| 精品一区二区三区av网在线观看| 欧美最黄视频在线播放免费| 18+在线观看网站| 久久鲁丝午夜福利片| h日本视频在线播放| 亚洲成人av在线免费| av.在线天堂| 日韩av不卡免费在线播放| 18禁黄网站禁片免费观看直播| 99久久中文字幕三级久久日本| 国产男人的电影天堂91| 俄罗斯特黄特色一大片| 嫩草影院入口| 国产精品美女特级片免费视频播放器| 色视频www国产| 国产aⅴ精品一区二区三区波| 中文字幕av在线有码专区| 老熟妇乱子伦视频在线观看| 寂寞人妻少妇视频99o| 国产一区二区三区av在线 | 亚洲激情五月婷婷啪啪| 欧美日韩精品成人综合77777| 国产精品人妻久久久影院| 免费看美女性在线毛片视频| 欧美成人a在线观看| 麻豆成人午夜福利视频| 国产精品久久久久久久久免| 久久久精品大字幕| 国产探花极品一区二区| 寂寞人妻少妇视频99o| 久久国产乱子免费精品| 黑人高潮一二区| 日本免费a在线| 午夜精品在线福利| 成人鲁丝片一二三区免费| 国产又黄又爽又无遮挡在线| 日韩在线高清观看一区二区三区| 欧美最新免费一区二区三区| 男人舔奶头视频| 国产蜜桃级精品一区二区三区| 久久6这里有精品| 欧美丝袜亚洲另类| 欧美高清性xxxxhd video| 国产精品亚洲美女久久久| 国产成人精品久久久久久| 国产午夜精品久久久久久一区二区三区 | 国产综合懂色| 亚洲成av人片在线播放无| 人人妻人人澡人人爽人人夜夜 | 一卡2卡三卡四卡精品乱码亚洲| 国产av麻豆久久久久久久| 国产精品99久久久久久久久| 22中文网久久字幕| 亚州av有码| 日韩av在线大香蕉| 美女高潮的动态| 久久久欧美国产精品| 人妻夜夜爽99麻豆av| 18禁黄网站禁片免费观看直播| 国产欧美日韩精品亚洲av| 一区二区三区免费毛片| 女的被弄到高潮叫床怎么办| 国产激情偷乱视频一区二区| 亚洲av五月六月丁香网| 国产精品福利在线免费观看| 国产成年人精品一区二区| 在线观看美女被高潮喷水网站| 免费观看的影片在线观看| 亚洲国产精品国产精品| 日韩,欧美,国产一区二区三区 | 91久久精品国产一区二区三区| 久久午夜福利片| 国产亚洲欧美98| 亚洲精华国产精华液的使用体验 | 色5月婷婷丁香| 日本免费一区二区三区高清不卡| 日韩欧美三级三区| 夜夜爽天天搞| 欧美成人精品欧美一级黄| 国产欧美日韩精品一区二区| 欧美区成人在线视频| 亚洲欧美日韩无卡精品| 综合色av麻豆| 国产精品美女特级片免费视频播放器| 国产成人a区在线观看| 亚洲成人av在线免费| 黄片wwwwww| 国产一区二区激情短视频| 亚洲aⅴ乱码一区二区在线播放| 三级毛片av免费| 99热精品在线国产| 色哟哟哟哟哟哟| 免费电影在线观看免费观看| 狠狠狠狠99中文字幕| 精品人妻熟女av久视频| 一本精品99久久精品77| 久久久久国产网址| av在线观看视频网站免费| 在线观看66精品国产| 亚洲欧美精品自产自拍| 成人国产麻豆网| av在线播放精品| 美女高潮的动态| 成人精品一区二区免费| 亚洲一区二区三区色噜噜| 国产蜜桃级精品一区二区三区| 人妻少妇偷人精品九色| 两性午夜刺激爽爽歪歪视频在线观看| 美女大奶头视频| 嫩草影院精品99| 最近2019中文字幕mv第一页| 精品久久久久久久久久免费视频| 亚洲经典国产精华液单| 3wmmmm亚洲av在线观看| 天美传媒精品一区二区| 亚洲人与动物交配视频| 亚洲在线自拍视频| 在线天堂最新版资源| 久久这里只有精品中国| 日日干狠狠操夜夜爽| 久久久久性生活片| 亚洲人成网站在线观看播放| 日韩欧美 国产精品| 国产一区亚洲一区在线观看| 亚洲va在线va天堂va国产| 99热这里只有精品一区| 精品少妇黑人巨大在线播放 | 69av精品久久久久久| 中出人妻视频一区二区| 国产伦在线观看视频一区| 久久久午夜欧美精品| 女生性感内裤真人,穿戴方法视频| 久久国内精品自在自线图片| 老师上课跳d突然被开到最大视频| 一个人看的www免费观看视频| 97人妻精品一区二区三区麻豆| 日韩欧美免费精品| 麻豆av噜噜一区二区三区| 麻豆一二三区av精品| 麻豆久久精品国产亚洲av| 国产精品人妻久久久影院| 禁无遮挡网站| 精品久久久久久久久久免费视频| 夜夜看夜夜爽夜夜摸| 99久久中文字幕三级久久日本| 日韩欧美精品免费久久| 亚洲第一电影网av| .国产精品久久| 丰满的人妻完整版| 欧美成人一区二区免费高清观看| 精品99又大又爽又粗少妇毛片| 色综合色国产| 一a级毛片在线观看| 色播亚洲综合网| 亚洲国产精品久久男人天堂| 黄色日韩在线| 色综合亚洲欧美另类图片| 成人无遮挡网站| av黄色大香蕉| 亚洲中文字幕日韩| 黑人高潮一二区| 麻豆国产av国片精品| 99精品在免费线老司机午夜| 老熟妇仑乱视频hdxx| 久久精品综合一区二区三区| 狂野欧美白嫩少妇大欣赏| 免费av观看视频| 亚洲乱码一区二区免费版| 亚洲av一区综合| 特大巨黑吊av在线直播| 亚洲av五月六月丁香网| 国产蜜桃级精品一区二区三区| 在线观看午夜福利视频| 老女人水多毛片| 色哟哟哟哟哟哟| 国国产精品蜜臀av免费| 免费在线观看成人毛片| 欧美激情国产日韩精品一区| 婷婷色综合大香蕉| 18禁在线无遮挡免费观看视频 | eeuss影院久久| 欧美成人a在线观看| 老女人水多毛片| 悠悠久久av| 激情 狠狠 欧美| 成年女人毛片免费观看观看9| 婷婷亚洲欧美| 91麻豆精品激情在线观看国产| 亚洲av免费高清在线观看| 亚洲av中文字字幕乱码综合| 日本与韩国留学比较| 亚洲aⅴ乱码一区二区在线播放| 亚洲图色成人| 国产成人精品久久久久久| 成年女人永久免费观看视频| 成年免费大片在线观看| 伦理电影大哥的女人| 丰满人妻一区二区三区视频av| 乱系列少妇在线播放| 国产精品一区二区三区四区免费观看 | 国产精品精品国产色婷婷| 欧美极品一区二区三区四区| 亚洲一区二区三区色噜噜| 日韩强制内射视频| 亚洲人成网站在线播| 亚洲经典国产精华液单| 亚洲一级一片aⅴ在线观看| 国产精品嫩草影院av在线观看| 日日摸夜夜添夜夜添小说| 亚洲精品一区av在线观看| 日韩欧美免费精品| 成年av动漫网址| 国产亚洲精品久久久久久毛片| 国产欧美日韩精品亚洲av| 中文字幕精品亚洲无线码一区| 熟妇人妻久久中文字幕3abv| 亚洲美女视频黄频| 俺也久久电影网| 国产av麻豆久久久久久久| 欧美区成人在线视频| 亚洲欧美精品自产自拍| 观看美女的网站| 国产伦精品一区二区三区四那| av国产免费在线观看| 国产成人aa在线观看| 日本成人三级电影网站| 九九热线精品视视频播放| 夜夜夜夜夜久久久久| 亚洲精品一卡2卡三卡4卡5卡| 97超视频在线观看视频| 国产视频内射| 欧美区成人在线视频| 嫩草影院精品99| 99热这里只有精品一区| 人妻夜夜爽99麻豆av| 搡老妇女老女人老熟妇| 久久久久久久久久久丰满| 日韩欧美 国产精品| 欧美zozozo另类| 成人一区二区视频在线观看| 人人妻人人澡欧美一区二区| 成人无遮挡网站| 亚洲欧美成人精品一区二区| 国产私拍福利视频在线观看| 亚洲成人久久爱视频| 成人高潮视频无遮挡免费网站| 91麻豆精品激情在线观看国产| 欧美色欧美亚洲另类二区| 国产视频一区二区在线看| 国产精品日韩av在线免费观看| 免费黄网站久久成人精品| 欧美又色又爽又黄视频| 一级毛片我不卡| 日日摸夜夜添夜夜添小说| 亚洲欧美成人综合另类久久久 | 91狼人影院| 免费不卡的大黄色大毛片视频在线观看 | 91av网一区二区| 久久国内精品自在自线图片| 欧美人与善性xxx| 国产高清视频在线观看网站| 成人亚洲精品av一区二区| 亚洲成人av在线免费| 国产一区二区亚洲精品在线观看| 性欧美人与动物交配| 亚洲国产高清在线一区二区三| 超碰av人人做人人爽久久| 亚洲国产精品国产精品| 亚洲真实伦在线观看| 最近手机中文字幕大全| 欧美不卡视频在线免费观看| 有码 亚洲区| 内地一区二区视频在线| 久久精品人妻少妇| 亚洲精品久久国产高清桃花| 亚洲成av人片在线播放无| 99久久中文字幕三级久久日本| 国产精品一二三区在线看| 久久精品91蜜桃| 狠狠狠狠99中文字幕| 国产欧美日韩精品一区二区| 久久精品夜夜夜夜夜久久蜜豆| 51国产日韩欧美| 在线a可以看的网站| 久久综合国产亚洲精品| 国内精品久久久久精免费| 男人舔女人下体高潮全视频| 可以在线观看的亚洲视频| 精品不卡国产一区二区三区| 午夜亚洲福利在线播放| 国产欧美日韩精品一区二区| 久久久久久久久大av| 中出人妻视频一区二区| 大型黄色视频在线免费观看| 一级毛片我不卡| 看免费成人av毛片| 激情 狠狠 欧美| 免费观看精品视频网站| 男女做爰动态图高潮gif福利片| 国产亚洲欧美98| 欧美一区二区精品小视频在线| 天天躁日日操中文字幕| 亚洲,欧美,日韩| 成年女人永久免费观看视频| 观看免费一级毛片| 日韩精品有码人妻一区| 夜夜夜夜夜久久久久| ponron亚洲| 给我免费播放毛片高清在线观看| 亚洲av五月六月丁香网| 综合色丁香网| 国产亚洲精品综合一区在线观看| 欧美精品国产亚洲| 级片在线观看| 一区福利在线观看| 免费看av在线观看网站| 亚洲精品日韩在线中文字幕 | 丝袜美腿在线中文| 搡老熟女国产l中国老女人| 麻豆国产97在线/欧美| 久久久久久久亚洲中文字幕| 草草在线视频免费看| 老司机福利观看| 国产真实伦视频高清在线观看| 非洲黑人性xxxx精品又粗又长| 悠悠久久av| 日本欧美国产在线视频| 99视频精品全部免费 在线| 卡戴珊不雅视频在线播放| 成人永久免费在线观看视频| 淫秽高清视频在线观看| 日韩一区二区视频免费看| 国产高清视频在线观看网站| 国产在视频线在精品| 亚洲美女视频黄频| 女的被弄到高潮叫床怎么办| 亚洲乱码一区二区免费版| 秋霞在线观看毛片| 国产精品综合久久久久久久免费| 给我免费播放毛片高清在线观看| 欧美日韩综合久久久久久| 天天躁日日操中文字幕| 欧美不卡视频在线免费观看| 亚洲av第一区精品v没综合| 国产伦在线观看视频一区| 亚洲性夜色夜夜综合| 国产精品人妻久久久久久| 91久久精品国产一区二区三区| 亚洲丝袜综合中文字幕| 久久人人爽人人爽人人片va| 久久久a久久爽久久v久久| 我的女老师完整版在线观看| 久久久久久久久久久丰满| 99久久精品热视频| 欧美国产日韩亚洲一区| 男女做爰动态图高潮gif福利片| 精品久久久久久成人av| 美女免费视频网站| 老司机福利观看| 色噜噜av男人的天堂激情| 特大巨黑吊av在线直播| 国产一区二区在线观看日韩| 成熟少妇高潮喷水视频| 久久精品国产自在天天线| 特大巨黑吊av在线直播| 3wmmmm亚洲av在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产一区二区在线观看日韩| 午夜福利成人在线免费观看| 最新中文字幕久久久久| 欧美高清性xxxxhd video| 久久这里只有精品中国| 别揉我奶头 嗯啊视频| 久久久成人免费电影| 一本精品99久久精品77| 免费不卡的大黄色大毛片视频在线观看 | 一级毛片电影观看 | 22中文网久久字幕| 欧美xxxx性猛交bbbb| 神马国产精品三级电影在线观看| 啦啦啦啦在线视频资源| 男女下面进入的视频免费午夜| 观看免费一级毛片| 男女边吃奶边做爰视频| 日日摸夜夜添夜夜添小说| 欧美性猛交黑人性爽| 国产亚洲精品av在线| 午夜亚洲福利在线播放|