修改稿雅鲁藏布江流域径流变化特征分析

雅鲁藏布江流域径流变化特征及趋势分析基金项目：国家自然科学基金（40830639）第一作者简介：洛珠尼玛基金项目：国家自然科学基金（40830639）第一作者简介：洛珠尼玛（1966－），男（藏族），西藏那曲人，工程师，河海大学西藏特培生，主要研究方向为流域水文过程模拟及预报。洛珠尼玛1，王建群2，徐幸仪2（1.西藏自治区水文局，西藏拉萨850000；2.河海大学水文水资源学院，江苏南京210098）摘要：本文采用Mann-Kendall秩次相关检验法等分析方法对雅鲁藏布江流域径流变化特征及年降水量、年蒸发皿观测系列进行了详细研究。雅鲁藏布江上中下游控制站径流上世纪60年代最丰，80年代最枯，年平均流量过程总体上呈略有下降趋势，但不显著o1960〜1981年，中水年和丰水年份较多；1982〜1997年，中水年和枯水年份较多；自1998年以来，丰水年份的机会较多。径流量主要集中在汛期（7〜9月），其中，8月份径流量最大，2月〜3月份径流量最小，径流量月最大值是月最小值的约12倍；径流年内分配不均匀系数和7〜9月占全年比例年际间较稳定。雅鲁藏布江流域降水量和径流变化不显著，气温和蒸发呈明显增加趋势。关键词：气候变化；青藏高原；水循环；径流；降水；蒸发1引言雅鲁藏布江流域界于东经82°00^〜97°07‘和北纬28°00‘〜31°16‘之间，流域面积24.2万km2；东西最大长度约1500km,南北最大宽度约290km,平均宽度约166km；雅鲁藏布江流域发源于喜马拉雅山北麓杰马央宗冰川，源头海拔5590m。青藏高原因其高海拔、地表物理性质复杂等特点，在全球气候变化中具有其自身特殊的演变规律［1，2］，是全球变化研究的热点区域。研究位于青藏高原腹地的雅鲁藏布江流域径流变化特征及趋势，对全球变化研究有着重要的意义。关于西藏高原地区气候变化趋势和雅鲁藏布江流域径流变化特征，有着大量的研究。文献［3-6］的研究一致认为，20世纪60年代以来西藏高原年平均气温全面上升。文献［7］利用西藏1971-2000年月降水量资料分析了近30年高原降水的变化趋势，结果发现西藏大部分地区年降水量变化为正趋势，降水倾向率为1.4-66.6mm/10a,而阿里地区呈较为明显的减少趋势。文献［8］利用西藏高原近40年来的逐月气象数据分析了近40年来西藏高原气候变化的主要特征，认为降水趋势变率空间分布上的基本规律是，其大小由东往西逐渐减小,藏中和藏东为上升趋势，藏西为下降趋势。文献［9,10］分别采用Morlet小波分析了雅鲁藏布江干流径流年际变化和季节性变化的多时间尺度演变规律,揭示了不同尺度下径流变化周期。文献［11］利用雅鲁藏布江流域1956-2000年近45年的实测径流资料对雅鲁藏布江流域径流变化特征进行了研究，认为雅鲁藏布江流域径流的年际变化较稳定,年内分配极不平衡。文献［12］选取雅鲁藏布江流域1956-2008年各干支流10个代表站的年径流实测资料，分析年径流的变化特征及其内部相依性，并进行趋势分析和周期成分分析，认为雅鲁藏布江流域径流的年际变化较小，年内分配极不平衡，径流受季节影响较大。文献［13］运用Mann-Kendall非参数检验方法,分析了雅鲁藏布江流域1956-2000年的年径流演变规律及其驱动因子，认为雅鲁藏布江流域径流总体上呈减少的趋势,演变过程中表现出明显的阶段性和突变性。文献［14］采用Hilbert-Huang变换对雅鲁藏布江干支流上的10个水文站1956-2000年天然径流资料进行分析,探讨了雅鲁藏布江流域年径流变化的近似周期及其演变趋势，发现雅鲁藏布江干支流在20世纪60年代左右径流量较大,但在进入60年代中后期之后,年径流量呈逐渐减小的趋势,这种趋势在80年代附近达到极小值,此后径流量逐渐回升,并在90年代后逐渐进入相对丰水期。关于雅鲁藏布江流域水循环特征和演变趋势尽管有着大量的研究，但是，基于近期最新水文气象观测，研究雅鲁藏布江流域气候变化背景下的径流变化特征和趋势，仍然是一个十分值得关注的研究课题。本文将利用雅鲁藏布江干流上中下游三个水文站的数据，对雅鲁藏布江流域径流变化特征和演变趋势进行分析研究。2径流年际变化特征雅鲁藏布江干流上自上而下分布着奴各沙、羊村、奴下等水文站，其中奴各沙站在年楚河入雅江口的下游，积水面积106060km2；羊村站在拉萨河入雅江口的下游，积水面积151507km2；奴下站在尼洋河入雅江口的下游，积水面积191235km2。根据主要干流控制站奴各沙、羊村、奴下水文站，1960〜2009年年平均流量观测资料统计分析得雅鲁藏布江流域年径流量年际变化统计见表1。表1 雅鲁藏布江流域年径流量年际变化分析Table1AnalysisoftheInterannualvariationofannualrunoff站名1960〜1969平均年径流量（10sm3）1970〜1979平均年径流量（10sm3）1980〜1989平均年径流量（10sm3）1990〜1999平均年径流量（10sm3）2000〜2009平均年径流量（10sm3）1960〜2009平均年径流量（10sm3）奴各沙188.84160.01136.08151.72174.74162.28羊村340.59280.51249.70291.14336.74298.90由表1可以看出，雅江奴下站多年平均径流量597.96X108/3。，雅江流域上世纪60年代最丰、80年代最枯。采用Mann-Kendall秩次相关检验法［⑸、线性回归分析法对雅鲁藏布江主要控制站1960〜2009的年平均流量系列进行统计分析，结果见表2和图1〜图3。在表1中，M-KU表示Mann-Kendall秩次相关检验统计量U,Mann-Kendall倾斜度0，具体定义如下。对时间序列x,x,x(n为样本长度)，定义Mann-Kendall秩次相关检验12n统计量U如下：S-1丸Var(S)1)TOC\o"1-5"\h\zU=\ 01)S+1JVar(S)S=艺工Sgn(x一x) (2)kii=1k=i+1Var(S)=[n(n-1)(2n+5)-£t(t-1)(2t+5)]/18 (3)jjjj=1其中m为“结”的总数，t.,j=1,2,…,m表示第j个“结”的宽度。当n增大时,jU很快收敛于标准正态分布。给定显著性水平«=0.05，其双边正态分位数值U=1.96，当|U|>U时，时间序列变化趋势显著。称统计量U为a/2 a/2Mann-Kendall秩次相关系数。还定义Mann-Kendall倾斜度0如下:(x-x)0=Median一^——j ,Vj<i (4)Ii_j丿其中1<j<i<n。当0>0时，表示上升趋势；当0<0时，表示下降趋势。表2 雅鲁藏布江主要控制站年平均流量系列趋势分析Table2Long-termchangestrendofannualaveragedischarge水文站M-KUM-KB线性趋势((mVs)/a)统计年份奴各沙-0.93-1.67-1.461960-20091-s量流年份1-s量流年份年平均流量—多年平均值一线性（年平均流量）羊村-0.16-0.50-0.211960-2009奴下-0.47-0.67-1.321960-2009图1奴各沙站年平均流量系列变化趋势Fig.1ChangestrendofannualaveragedischargeofNugesha1s量流年份—年平均流量—多年平均值—1s量流年份—年平均流量—多年平均值—线性（年平均流量）图2羊村站年平均流量系列变化趋势Fig.2ChangestrendofannualaveragedischargeofYangcun

年份丨^^年平均流量 多年平均值—线性（年平均流量）图3奴下站年平均流量系列变化趋势Fig.3ChangestrendofannualaveragedischargeofNuxia由表2和图1〜图3可知，雅鲁藏布江主要控制站年平均流量过程总体上略有下降趋势。M-K统计量U=-0.16〜-0.93，其绝对值没有超过显著性水平为a=0.05的双边正态分位数值1.96，因此主要控制站年平均流量过程下降趋势不显著。为进一步探讨年径流丰枯变化的规律，首先作出最下游奴下站年径流量P-III型理论频率曲线，设小于25%为丰水年，大于75%为枯水年，频率在25〜75%为中水年，计算得出丰水年与中水年的分界流量分别为2100m3/s和1610m3/s；年径流系列进行丰、中、枯分类，并绘图见图4。Fig.4TheprocessofwetyearanddryyearinNuxia

由图4可以看出，1960〜1981年，中水年和丰水年份较多；1982〜1997年，中水年和枯水年份较多；自1998年以来，出现丰水年份的机会较多，出现中水年份的机会较少，只有2009年是枯水年。3径流年内变化特征雅鲁藏布江主要控制站1960-2009年多年月平均流量见图5。由图5可以看出，径流量主要集中在汛期（7〜9月），8月份径流量最大，上游代表奴各沙站2月份径流量最小，中游和下游代表站羊村、奴下站3月份径流量最小，三站径流量月最大值与月最小值的倍数分别为12.3、12.6和12.2。1 2 3 4 56 7 8 9 10 11 12月 份务—星流径均平月年多1 2 3 4 56 7 8 9 10 11 12月 份务—星流径均平月年多奴各沙图5雅鲁藏布江主要控制站多年月平均流量Fig.5Theaverageofmanyyearsofmonthlyrunoff为了进一步分析径流年内分配特征变化规律，采用不均匀系数Cu和7〜9月占全年比例Cx来刻画径流年内分配特征，采用Mann-Kendall统计量U和倾斜度B、线性趋势提取发来分析Cu和Cx的变化趋势。其中年径流年内分配不均匀系数Cu定义如下：Cu=-E5）c2=I£(W-E)2

i5）i=1E=丄县W12ii=1其中W.为各月径流量，i=1,2,-12。对下游奴下站径流年内分配不均匀系i数系列和径流7〜9月占全年比例系列进行统计分析，结果见表3、表4。由表

3、表4可以看出，径流年内分配不均匀系数和7~9月占全年比例年际间变化较稳定。表3雅鲁藏布江主要控制站径流年内分配不均匀系数系列趋势分析Table3Long-termchangestrendofCu水文站M-KUM-KB线性趋势统计年份奴各沙-0.6859-0.0013-0.00061960-2009羊村-1.3016-0.0019-0.00171960-2009奴下-0.7499-0.0007-0.00071960-2009表4雅鲁藏布江主要控制站径流7〜9月占全年比例系列趋势分析Table4Long-termchangestrendofCx水文站M-KUM-KB线性趋势统计年份奴各沙-0.1004-0.0001-0.000051960-2009羊村-0.8534-0.0005-0.00051960-2009奴下-0.5603-0.0003-0.00031960-20094年降水和蒸发变化趋势采用Mann-Kendall秩次相关检验法、线性回归分析法对雅鲁藏布江上中下游代表站，1960〜2009年年降水量、20cm蒸发皿蒸发量系列进行统计分析，结果见表5、表6。表5雅鲁藏布江代表站降水量变化趋势分析Table5Long-termchangestrendofannualprecipitation水文站M-KUM-KB线性趋势统计年份奴各沙-1.79-1.94-2.161960-2009羊村-0.38-0.53-0.671960-2009奴下1.070.990.991960-2009表6雅鲁藏布江代表站蒸发量变化趋势分析Table6Long-termchangestrendofannualpotentialevaporation水文站M-KUM-KB线性趋势统计年份奴各沙6.231.5527.891960-2009羊村1.065.075.171960-2009奴下3.514.606.301960-2009由表5、表6可知，上、中游代表站年降水量呈下降趋势，下游代表站呈增加趋势，但均不显著；年蒸发量为显著增加趋势。5结语本文采用Mann-Kendall秩次相关检验法、线性回归分析法、丰中枯分类、不均匀性度量等方法对雅鲁藏布江上中下游控制站的年径流特性和年内分配特征以及年降水量、年蒸发皿观测系列进行了详细研究，结论如下：雅鲁藏布江上中下游控制站径流变化上世纪60年代最丰，80年代最枯，下游控制奴下站多年平均径流量597.96X108/3。上中下游控制站的年平均流量过程总体上略有下降趋势，但不显著。下游控制奴下站为例，1960〜1981年，中水年和丰水年份较多；1982〜1997年，中水年和枯水年份较多；自1998年以来，出现丰水年份的机会较多，出现中水年的机会较少，只有2009年是枯水年。上中下游控制站平均径流量主要集中在汛期(7〜9月)，8月份径流量最大，2月~3月份径流量最小，三站径流量月最大值与月最小值的倍数分别为12.3、12.6和12.2。径流年内分配不均匀系数和7~9月占全年比例年际间较稳定。上、中游代表站年降水量呈下降趋势，下游代表站呈增加趋势，但均不显著；年蒸发量呈显著增加趋势。雅鲁藏布江流域降水量和径流变化不显著，气温和蒸发呈明显增加趋势。参考文献：郑度.青藏高原自然地域系统研究[J].中国科学D辑，1999,26(4)：336-341.(ZhengDu.StudyonRegionalPhysio-geographicalSysteminQinghai-TibetanPlateau[J].ScienceinChina(SeriesD),1999.(inChinese))李潮流，康世昌.青藏高原不同时段气候变化的研究综述[J].地理学报，2006，61(3)：327-335.(LIChao-liu,KANGShi-chang.ReviewofStudiesinClimateChangeovertheTibetanPlateau[J].ActaGeographicaSinica,2006.(inChinese))3周顺武，张勇.近42年雅鲁藏布江中游四季气温变化特征分析[J].成都气象学院学报，2000,15(4)：307-311.(ZH0UShun-wu,ZHANGYong.AnalysisofseasonaltemperaturevariationovermiddlereachesofYalutsangpoRiverduringpast42years[J].JournalofChengduInstituteofMeteorology,2000.(inChinese))杜军.西藏高原近40年的气温变化[J].地理学报，2001，56(6)：682-690.(DUJun.ChangeofTemperatureinTibetanPlateauFrom1961to2000[J].ActaGeographicaSinica,2001.(inChinese))盛文萍，高清竹，李玉娥，等.藏北地区气候变化特征及其影响分析[J].高原气象，2008，27(3)：509-516.(SHENGWen-ping,GAOQing-zhu,LIYu-e,WANYun-fan,YANGKai.CharacteristicofClimateChangeinNorthernTibetandItsImpact[J].PlateauMeteorology,2008.(inChinese))李发鹏，徐综学.西藏高原气象要素长期变化趋势及其持续性特征分析[A].夏军，贾绍凤，刘苏峡.第七届中国水论坛论文集-水系统与水资源可持续管理[C].中国水利水电出版社，2010，4：104-109.(LiFa-peng,XuZong-xue.Long-teamTrendandDurativeCharacteristicsforMeteorologyElementsofQinghai-TibetanPlateau[A].XiaJun,JiaShao-Feng,LiuSu-xia.SustainableManagementofWaterSystemandWaterResource[C].ChinaWaterPowerPress,2010.(inChinese))杜军，马玉才.西藏高原降水变化趋势的气候分析[J].地理学报，2004，59(3)：375-382.(DUJun,MAYu-cai.ClimaticTrendofRainfalloverTibetanPlateaufrom1971to2000[J].ActaGeographicaSinica,2004.(inChinese))徐宗学，巩同梁，赵芳芳.近40年来西藏高原气候变化特征分析[J].亚热带资源与环境学报，2006，1(1)：24-32.(XUZong-xue,GONGTong-liang,ZHAOFang-fang.AnalysisofClimateChangeinTibetanPlateauoverthePast40Years[J].JournalofSubtropicalResourcesandEnvironment,2006.(inChinese))李占玲，徐宗学，巩同梁•雅鲁藏布江流域径流特性变化分析 [J].地理研究,2008,27(2):353-361.(LIZhanTing,XUZong-xue,GONGTong-liang.CharacteristicsofrunoffvariationintheYarlungZangboRiverbasin[J].GeographicalResearch,2008.(inChinese))钱晓燕,袁鹏,邵骏,颜志衡.雅鲁藏布江河川径流变化的季节性规律研究[J].水资源与水工程学报,2010,21(l):29-33.(QIANXiao—yan,YUANPeng,SHAOJun,YANZhi-heng.RuleoftheSeasonalRunoffVariationintheYarlungZangboRiver[J].JournalofWaterResourcesandWaterEngineering,2010.(inChinese))刘剑，姚治君，陈传友•雅鲁藏布江径流变化趋势及原因分析[J].自然资源学报,2007,22(3):471—477.(LIUJian,YAOZhi—jun1,CHENChuan—you.EvolutionTrendandCausationAnalysisoftheRunoffEvolutionintheYarlungZangboRiverBasin[J].JournalofNaturalResources,2007.(inChinese))黄艳，袁鹏，颜志衡.雅鲁藏布江流域天然径流特征及规律分析[J].红水河，2010,29(1):65—69.(HUANGYan,YUANPeng,YANZhi—heng.AnalysisonCharacteristicsandLawsofYarlungZangboRiverBasin'sNaturalRunoff[J].HongshuiRiver,2010.(inChinese))黄俊雄，徐宗学，巩同梁.雅鲁藏布江径流演变规律及其驱动因子分析[J].水文，2007,27(5):31-35. (HUANGJun—xiong,XUZong—xue,GONGTong—liang.CharacteristicsandDrivingFactorsoftheRunoffVariationsintheYarlungZangboRiver[J].JournalofChinaHydrology,2007.(inChinese))邵骏，袁鹏，颜志衡，张文江•基于HHT的雅鲁藏布江径流变化周期及趋势分析.中山大学学报(自然科学版)，2010，49(1)：125—130.(SHAOJun,YUANPeng,YANZhi—heng,ZHANGWen—jiang.AnalysisonPeriodicVariationsandTrendofAnnualRunoffintheYarlungZangboRiverwithHilbert—HuangTransform.ActaScientiarumNaturaliumUniversitatisSunyatseni,2010.(inChinese))魏凤英.现代气候统计诊断与预测技术[M].北京：气象出版社，1999:69—72.(WeiFeng—ying.Statisticaldiagnosisandpredictiontechnologyofmodernclimate[M].Beijing:MeteorologicalPress,1999.(inChinese))ChangeCharacterandTrendofRunoffinYarlungZangboRiverBasinLuoZhuNi-ma1(HydrologyBureauofTibetAutonomousRegion,Lhasa850000,China)WangJian-qun,XuXing-yi(CollegeofHydrologyandWaterResources,HohaiUniversity,Nanjing

210098,China)Abstract:Inthispaper,thechangecharacterofrunoff,annualprecipitationandtheseriesofannualevaporationobservedbyevaporatingdishareanalyzedbyusingMann-KendalltestandothermethodsintheYarlungZangboRiverBasin.Themostabundantannualrunoffhappenedinthe196'0s,the1980'sannualrunoffisthemostinadequate/insufficient.Thereislittledown