• <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ì)”
    亚洲国产毛片av蜜桃av| 国产激情欧美一区二区| 黄色a级毛片大全视频| 99国产精品99久久久久| ponron亚洲| avwww免费| 天堂动漫精品| 亚洲一码二码三码区别大吗| av网站在线播放免费| 亚洲色图综合在线观看| 女同久久另类99精品国产91| 一进一出好大好爽视频| 男女之事视频高清在线观看| 69精品国产乱码久久久| av天堂久久9| 精品久久久精品久久久| 久久婷婷成人综合色麻豆| 国产高清视频在线播放一区| 国产精品二区激情视频| 国产精华一区二区三区| 91麻豆av在线| 黄色怎么调成土黄色| 国产精品久久久久成人av| 欧美性长视频在线观看| 一a级毛片在线观看| 欧美日韩乱码在线| 久久精品国产99精品国产亚洲性色 | 久久久久久久精品吃奶| 99久久精品国产亚洲精品| 亚洲精品av麻豆狂野| 午夜激情av网站| 国产又爽黄色视频| 久久久久久久久久久久大奶| 亚洲第一欧美日韩一区二区三区| 中文字幕精品免费在线观看视频| 纯流量卡能插随身wifi吗| 亚洲熟妇中文字幕五十中出 | 大片电影免费在线观看免费| 岛国在线观看网站| 国产淫语在线视频| 69av精品久久久久久| e午夜精品久久久久久久| 日日摸夜夜添夜夜添小说| 亚洲五月色婷婷综合| 免费黄频网站在线观看国产| 丝袜美腿诱惑在线| 女人被狂操c到高潮| 咕卡用的链子| 午夜视频精品福利| 亚洲精品美女久久av网站| 久久精品人人爽人人爽视色| aaaaa片日本免费| 黄片播放在线免费| 九色亚洲精品在线播放| 亚洲视频免费观看视频| 制服诱惑二区| 脱女人内裤的视频| av中文乱码字幕在线| 日韩熟女老妇一区二区性免费视频| 亚洲五月婷婷丁香| 黄片播放在线免费| www.自偷自拍.com| 国产亚洲欧美在线一区二区| 免费观看a级毛片全部| 久久久久久久久免费视频了| 少妇粗大呻吟视频| av有码第一页| 午夜精品久久久久久毛片777| 国产成人av激情在线播放| 久久性视频一级片| 亚洲欧美日韩高清在线视频| 亚洲精品国产区一区二| 露出奶头的视频| 欧美成狂野欧美在线观看| 久久久国产一区二区| 制服人妻中文乱码| 亚洲国产欧美一区二区综合| 亚洲黑人精品在线| 欧美日韩成人在线一区二区| 欧美成狂野欧美在线观看| 免费在线观看亚洲国产| 一级作爱视频免费观看| 国产区一区二久久| 欧美成人免费av一区二区三区 | 99香蕉大伊视频| 精品国产一区二区三区四区第35| 久久国产精品男人的天堂亚洲| 一区二区日韩欧美中文字幕| 欧美黄色片欧美黄色片| 亚洲精品美女久久av网站| 99精国产麻豆久久婷婷| 如日韩欧美国产精品一区二区三区| 久久精品国产99精品国产亚洲性色 | 成年动漫av网址| 满18在线观看网站| 搡老乐熟女国产| 国产亚洲精品久久久久久毛片 | 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 99国产精品一区二区蜜桃av | 国产欧美日韩精品亚洲av| 国产欧美日韩一区二区精品| 精品久久久精品久久久| 久久精品国产亚洲av香蕉五月 | 午夜精品国产一区二区电影| 久久中文字幕一级| 国产精品一区二区在线不卡| 一进一出抽搐gif免费好疼 | 国产成人av教育| 男男h啪啪无遮挡| 十八禁高潮呻吟视频| 露出奶头的视频| 国产区一区二久久| 精品少妇一区二区三区视频日本电影| 99精品欧美一区二区三区四区| 欧美成狂野欧美在线观看| 亚洲精品av麻豆狂野| 亚洲精品在线观看二区| 国产伦人伦偷精品视频| 欧美另类亚洲清纯唯美| 丰满饥渴人妻一区二区三| 欧美激情 高清一区二区三区| 成人黄色视频免费在线看| 在线观看免费高清a一片| 人人妻,人人澡人人爽秒播| 久久中文字幕人妻熟女| 免费av中文字幕在线| 亚洲精品成人av观看孕妇| 亚洲成av片中文字幕在线观看| 麻豆国产av国片精品| 日日夜夜操网爽| 午夜免费成人在线视频| 俄罗斯特黄特色一大片| 中文字幕最新亚洲高清| 自线自在国产av| 国产免费现黄频在线看| 久久久国产一区二区| 亚洲专区国产一区二区| 90打野战视频偷拍视频| 亚洲精品国产一区二区精华液| 久久久久国产精品人妻aⅴ院 | 黄色毛片三级朝国网站| 成人精品一区二区免费| 国产又色又爽无遮挡免费看| 欧美成狂野欧美在线观看| 久久国产亚洲av麻豆专区| 日本五十路高清| 丝袜美腿诱惑在线| 久久热在线av| 搡老岳熟女国产| 最近最新免费中文字幕在线| 黄色毛片三级朝国网站| 久热这里只有精品99| 国产精华一区二区三区| 亚洲成人手机| 国产又色又爽无遮挡免费看| 欧美人与性动交α欧美精品济南到| netflix在线观看网站| 波多野结衣av一区二区av| av超薄肉色丝袜交足视频| 成人18禁高潮啪啪吃奶动态图| 女人被躁到高潮嗷嗷叫费观| 69av精品久久久久久| 黄片播放在线免费| 巨乳人妻的诱惑在线观看| 午夜福利在线观看吧| 91老司机精品| 村上凉子中文字幕在线| 在线看a的网站| 欧美精品av麻豆av| 欧美亚洲日本最大视频资源| 亚洲七黄色美女视频| 制服人妻中文乱码| 婷婷精品国产亚洲av在线 | 亚洲中文字幕日韩| 色94色欧美一区二区| 狂野欧美激情性xxxx| 国产精品免费视频内射| 日韩视频一区二区在线观看| 黑人欧美特级aaaaaa片| 亚洲av成人不卡在线观看播放网| 一级片免费观看大全| 夜夜躁狠狠躁天天躁| 亚洲伊人色综图| cao死你这个sao货| xxx96com| 搡老岳熟女国产| 国产色视频综合| 咕卡用的链子| 一进一出好大好爽视频| 日本vs欧美在线观看视频| 午夜激情av网站| 国产午夜精品久久久久久| 久久久久国产一级毛片高清牌| 日韩中文字幕欧美一区二区| 欧美日韩亚洲综合一区二区三区_| 欧美中文综合在线视频| 自线自在国产av| 男女之事视频高清在线观看| 免费一级毛片在线播放高清视频 | 成人精品一区二区免费| 国产欧美日韩一区二区三| 99热国产这里只有精品6| 国产精品乱码一区二三区的特点 | 亚洲片人在线观看| 亚洲免费av在线视频| av免费在线观看网站| 日韩熟女老妇一区二区性免费视频| 成熟少妇高潮喷水视频| 12—13女人毛片做爰片一| 在线看a的网站| 黄色怎么调成土黄色| 丁香欧美五月| 欧美色视频一区免费| 欧美性长视频在线观看| 久久久国产精品麻豆| 欧洲精品卡2卡3卡4卡5卡区| 人人妻,人人澡人人爽秒播| 亚洲一区二区三区欧美精品| 99国产精品一区二区三区| 国产欧美日韩精品亚洲av| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜在线中文字幕| 天天操日日干夜夜撸| 免费一级毛片在线播放高清视频 | 在线视频色国产色| 老鸭窝网址在线观看| 黄色怎么调成土黄色| 久久国产精品男人的天堂亚洲| 一区在线观看完整版| 欧美激情久久久久久爽电影 | 国产不卡一卡二| 欧美日韩视频精品一区| 欧美日韩精品网址| 精品国产乱码久久久久久男人| 我的亚洲天堂| 欧美精品高潮呻吟av久久| 久久精品91无色码中文字幕| 老司机亚洲免费影院| 一二三四社区在线视频社区8| 久久久久视频综合| 好男人电影高清在线观看| 久久久久久久久久久久大奶| 一区二区三区国产精品乱码| 欧美黑人欧美精品刺激| 国产男女内射视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产极品粉嫩免费观看在线| 建设人人有责人人尽责人人享有的| 国产一区二区三区在线臀色熟女 | 亚洲五月色婷婷综合| 国产有黄有色有爽视频| 亚洲中文日韩欧美视频| 欧美激情 高清一区二区三区| 国产单亲对白刺激| 成人国产一区最新在线观看| 午夜精品久久久久久毛片777| 午夜激情av网站| 国产av精品麻豆| 我的亚洲天堂| a级片在线免费高清观看视频| 亚洲免费av在线视频| 老汉色∧v一级毛片| 人成视频在线观看免费观看| 99国产精品一区二区三区| 欧美激情极品国产一区二区三区| 最近最新中文字幕大全电影3 | 欧美日韩福利视频一区二区| 成年版毛片免费区| 国内毛片毛片毛片毛片毛片| 王馨瑶露胸无遮挡在线观看| 久久香蕉国产精品| 看免费av毛片| 亚洲精品美女久久久久99蜜臀| 午夜成年电影在线免费观看| 在线观看一区二区三区激情| 女人精品久久久久毛片| 亚洲av成人不卡在线观看播放网| 久久久久国内视频| 黄片播放在线免费| 村上凉子中文字幕在线| 色婷婷av一区二区三区视频| 激情视频va一区二区三区| 亚洲精品中文字幕一二三四区| 校园春色视频在线观看| 999久久久精品免费观看国产| cao死你这个sao货| 久久青草综合色| 精品人妻1区二区| 国产精品成人在线| 一边摸一边做爽爽视频免费| 午夜福利在线观看吧| 欧美 亚洲 国产 日韩一| 中文亚洲av片在线观看爽 | 亚洲精品成人av观看孕妇| 性色av乱码一区二区三区2| 高潮久久久久久久久久久不卡| 久久国产精品男人的天堂亚洲| 久久这里只有精品19| 不卡一级毛片| 国产精品 国内视频| 老司机深夜福利视频在线观看| 视频区欧美日本亚洲| 久久久久久久久免费视频了| 一进一出好大好爽视频| 久久久久精品国产欧美久久久| 天天影视国产精品| 女人精品久久久久毛片| 自线自在国产av| tocl精华| 黄色成人免费大全| 久久久久久久精品吃奶| 国产精品久久久久成人av| 成人亚洲精品一区在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲久久久国产精品| 法律面前人人平等表现在哪些方面| 精品人妻1区二区| 日韩有码中文字幕| 午夜视频精品福利| 99香蕉大伊视频| 热re99久久精品国产66热6| 久久天躁狠狠躁夜夜2o2o| 亚洲自偷自拍图片 自拍| 视频区欧美日本亚洲| 一级毛片女人18水好多| 欧美黑人欧美精品刺激| 夜夜躁狠狠躁天天躁| 国产91精品成人一区二区三区| 免费久久久久久久精品成人欧美视频| 美女午夜性视频免费| 久久久久久亚洲精品国产蜜桃av| 国产在线一区二区三区精| 久久香蕉激情| 国产欧美日韩一区二区三区在线| 日韩有码中文字幕| 女人被躁到高潮嗷嗷叫费观| 精品熟女少妇八av免费久了| 日韩免费av在线播放| 久久香蕉激情| 不卡一级毛片| 男人操女人黄网站| 999久久久精品免费观看国产| 1024视频免费在线观看| 水蜜桃什么品种好| 午夜福利在线观看吧| 水蜜桃什么品种好| 中国美女看黄片| 美女国产高潮福利片在线看| 91成年电影在线观看| 国产片内射在线| 人人妻,人人澡人人爽秒播| 大型av网站在线播放| 50天的宝宝边吃奶边哭怎么回事| 亚洲第一av免费看| 久久精品亚洲av国产电影网| 中文字幕高清在线视频| 亚洲欧美色中文字幕在线| 十八禁人妻一区二区| 亚洲av成人一区二区三| avwww免费| 欧美色视频一区免费| 欧美性长视频在线观看| 精品乱码久久久久久99久播| 精品国产一区二区久久| 下体分泌物呈黄色| 99在线人妻在线中文字幕 | 极品少妇高潮喷水抽搐| 欧美精品人与动牲交sv欧美| 久久青草综合色| 国产一区有黄有色的免费视频| 国产av又大| 热99国产精品久久久久久7| 人人妻人人爽人人添夜夜欢视频| 欧美久久黑人一区二区| 中文字幕人妻丝袜一区二区| 久久青草综合色| 天堂√8在线中文| 亚洲欧美日韩高清在线视频| 久久香蕉国产精品| 精品国产国语对白av| 日韩视频一区二区在线观看| 中文字幕精品免费在线观看视频| 亚洲一码二码三码区别大吗| 黄色成人免费大全| 18在线观看网站| 亚洲成人免费电影在线观看| 欧美成狂野欧美在线观看| 久久久久久免费高清国产稀缺| 成人手机av| 999久久久精品免费观看国产| 亚洲熟女精品中文字幕| 一进一出抽搐gif免费好疼 | 精品卡一卡二卡四卡免费| 久久久久久久午夜电影 | av线在线观看网站| 免费在线观看亚洲国产| 精品久久久久久,| 成年女人毛片免费观看观看9 | 人人妻人人澡人人看| 日韩免费av在线播放| 精品国产一区二区久久| 精品久久久久久,| 亚洲第一欧美日韩一区二区三区| 久久亚洲精品不卡| 久久国产精品人妻蜜桃| 大香蕉久久网| 香蕉久久夜色| 欧美在线一区亚洲| 久久香蕉精品热| 如日韩欧美国产精品一区二区三区| 美女视频免费永久观看网站| 国产精品秋霞免费鲁丝片| 韩国精品一区二区三区| av片东京热男人的天堂| 精品视频人人做人人爽| 久久影院123| 国产欧美日韩精品亚洲av| 日日爽夜夜爽网站| 国产一区二区三区在线臀色熟女 | 亚洲成人免费av在线播放| 两人在一起打扑克的视频| 国产一区有黄有色的免费视频| 成熟少妇高潮喷水视频| 宅男免费午夜| 怎么达到女性高潮| 久久久国产精品麻豆| 精品久久久精品久久久| 免费人成视频x8x8入口观看| 一级黄色大片毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 十八禁高潮呻吟视频| 可以免费在线观看a视频的电影网站| 欧美在线黄色| 国产aⅴ精品一区二区三区波| 亚洲精品美女久久久久99蜜臀| 一级毛片精品| 丝袜美腿诱惑在线| 777久久人妻少妇嫩草av网站| 满18在线观看网站| 黄色成人免费大全| 欧美日韩视频精品一区| 十八禁高潮呻吟视频| 老司机午夜福利在线观看视频| 无遮挡黄片免费观看| 欧美日韩亚洲国产一区二区在线观看 | 国产一区二区三区在线臀色熟女 | 在线观看一区二区三区激情| 国产精品永久免费网站| 日韩欧美三级三区| 日日摸夜夜添夜夜添小说| 极品教师在线免费播放| 亚洲色图av天堂| 日韩大码丰满熟妇| 少妇粗大呻吟视频| 欧美丝袜亚洲另类 | 久久精品亚洲av国产电影网| 91字幕亚洲| 夜夜夜夜夜久久久久| 国产不卡一卡二| 制服人妻中文乱码| 亚洲av日韩精品久久久久久密| 午夜福利影视在线免费观看| 亚洲五月天丁香| 欧美亚洲日本最大视频资源| 极品人妻少妇av视频| 欧美在线黄色| svipshipincom国产片| 久久中文看片网| 黄色怎么调成土黄色| 日韩熟女老妇一区二区性免费视频| 一本大道久久a久久精品| 国产精品 国内视频| 不卡一级毛片| 色婷婷av一区二区三区视频| 国产一区在线观看成人免费| 亚洲情色 制服丝袜| 90打野战视频偷拍视频| 下体分泌物呈黄色| 成人黄色视频免费在线看| 精品人妻1区二区| 熟女少妇亚洲综合色aaa.| 黄片小视频在线播放| 中文字幕av电影在线播放| 久久人妻福利社区极品人妻图片| 激情在线观看视频在线高清 | av天堂在线播放| 午夜久久久在线观看| 亚洲成av片中文字幕在线观看| 王馨瑶露胸无遮挡在线观看| 国产精品 国内视频| 黄色 视频免费看| 首页视频小说图片口味搜索| 在线永久观看黄色视频| 淫妇啪啪啪对白视频| 成在线人永久免费视频| 国产精品免费视频内射| 黄色片一级片一级黄色片| 国产又爽黄色视频| 亚洲情色 制服丝袜| 久久国产精品大桥未久av| 手机成人av网站| 欧美日韩视频精品一区| 丁香欧美五月| 国产在线精品亚洲第一网站| 丝袜美腿诱惑在线| 欧美成人午夜精品| 国产熟女午夜一区二区三区| 搡老岳熟女国产| 9色porny在线观看| 亚洲av电影在线进入| 亚洲色图av天堂| 久久香蕉激情| 美女午夜性视频免费| 丰满迷人的少妇在线观看| 黑丝袜美女国产一区| 黄色成人免费大全| 99久久国产精品久久久| 久久精品aⅴ一区二区三区四区| 怎么达到女性高潮| 日本a在线网址| 日本黄色视频三级网站网址 | 啦啦啦免费观看视频1| 国产精品一区二区免费欧美| 国产有黄有色有爽视频| 中文亚洲av片在线观看爽 | 91精品三级在线观看| 香蕉久久夜色| 亚洲人成电影观看| 乱人伦中国视频| 黑人猛操日本美女一级片| 成年人免费黄色播放视频| 女性被躁到高潮视频| 欧美色视频一区免费| 免费黄频网站在线观看国产| 制服人妻中文乱码| 叶爱在线成人免费视频播放| 天堂俺去俺来也www色官网| 久99久视频精品免费| 国产又爽黄色视频| 人人澡人人妻人| 十分钟在线观看高清视频www| 91精品国产国语对白视频| av一本久久久久| 乱人伦中国视频| 99re6热这里在线精品视频| 国产精品秋霞免费鲁丝片| 性色av乱码一区二区三区2| 国产精华一区二区三区| 亚洲综合色网址| 妹子高潮喷水视频| 精品乱码久久久久久99久播| 亚洲av欧美aⅴ国产| 成人18禁高潮啪啪吃奶动态图| 777久久人妻少妇嫩草av网站| 欧美黑人欧美精品刺激| 亚洲精品中文字幕一二三四区| 亚洲综合色网址| 日韩免费av在线播放| 搡老岳熟女国产| 亚洲第一青青草原| 日韩一卡2卡3卡4卡2021年| 老司机亚洲免费影院| 99热网站在线观看| xxx96com| 久久久精品区二区三区| 国产亚洲欧美在线一区二区| 18禁裸乳无遮挡动漫免费视频| 下体分泌物呈黄色| 精品乱码久久久久久99久播| 99热国产这里只有精品6| 久久 成人 亚洲| 国产精品香港三级国产av潘金莲| 日韩欧美一区视频在线观看| 欧美日韩福利视频一区二区| 亚洲国产欧美网| 曰老女人黄片| 日韩三级视频一区二区三区| 亚洲免费av在线视频| 12—13女人毛片做爰片一| 一区福利在线观看| 国产成人av教育| 久久香蕉精品热| 欧美av亚洲av综合av国产av| 亚洲人成电影免费在线| 久热爱精品视频在线9| 亚洲美女黄片视频| 99香蕉大伊视频| 不卡av一区二区三区| 啦啦啦免费观看视频1| 99热国产这里只有精品6| 精品一区二区三区视频在线观看免费 | 国产区一区二久久| 欧美+亚洲+日韩+国产| 夜夜爽天天搞| 一本综合久久免费| 午夜福利一区二区在线看| 国产成人免费观看mmmm| 黑人巨大精品欧美一区二区mp4| 亚洲精品美女久久av网站| 人妻一区二区av| av网站在线播放免费| 波多野结衣av一区二区av| 成年人免费黄色播放视频| 成人三级做爰电影| 超碰97精品在线观看| 69av精品久久久久久| 99国产精品一区二区蜜桃av | 成人国语在线视频| 又黄又爽又免费观看的视频| 日韩欧美免费精品| 正在播放国产对白刺激| 99国产精品一区二区三区| 免费高清在线观看日韩| 黄色a级毛片大全视频|