摘要：文中用河道来水距平百分率旱涝定义方法对区域内旱涝特征进行了分析。得出该区域丰枯转换周期大约为7～8年。并根据区域内58年径流实测数据建立了中长期自回归预报模型。预报方案精度达72.2%，精度检验达80.0%。从预测值可以确定安宁河下段2012年～2018年为丰水期，2019年～2026年为枯水期。安宁河流域2018年发生大洪水的可能性较大。2019年、2025年安宁河流域发生大旱或重旱的可能性较大。
关键词：等级；周期；径流深；枯水期；丰水期；预测
1 自然概况
安宁河是雅砻江下段左岸最大的一级支流，发源于四川省凉山彝族自治州冕宁县东小相岭，分东西两源：东源柯别河源于冕宁县拖乌乡拖乌伐木场东；西源中江河源于冕宁县大桥镇万年雪山。两源在大桥镇汇合后称安宁河。流域位于东经101°51′～102°48′、北纬26°38′～28°53′之间。两岸多山，支流多短小，比降较大，流域狭长。由北向南流经凉山彝族自治州的冕宁县、西昌市、德昌县及攀枝花市的米易县、盐边县，于盐边县安宁乡大坪地村湾滩处汇入雅砻江。安宁河全长337 km，干流长303 km，落差936m，平均比降3.1‰，流域面积11150km2。
安宁河流域地貌以山地为主，中段沿河两岸多河谷平原和山间盆地。流域地势北高南低，海拔3000m以上山地占流域面积36.5％，海拔1500~3000m山地占流域面积47.4％，海拔1500m以下的河谷平坝、山间盆地和二半山阶地仅占流域面积的16.1％。
安宁河流域面积大于500 km2的支流有4条，即孙水河、海河、茨达河、锦川河。小于500 km2大于100 km2的河流有14条。安宁河冕宁县安宁桥以上为上段，安宁桥至德昌县黄水塘为中段，黄水塘至河口为下段。
2 水文气象
安宁河流域属亚热带季风气候区，是西南季风暖湿气流北上的必经之路。流域内日照时间长，霜冻期短，夏秋多雨，冬春干旱，四季不明显。流域内气候温和，寒暑区别不大。除高寒山区外，多年平均气温均在14℃以上。
流域内降水量的年际变化不大，降水量的年内分配极不均匀，多年平均降水量在1000mm以上。5～10月受来自印度洋的西南季风和来自西太平洋的东南季风交替影响，降水量占全年的90%以上；11～4月降水量大多不足全年降水量的10%。流域北部河源及孙水河流域一带为安宁河流域第一暴雨中心；而德昌至米易小街一带为安宁河流域第二暴雨中心。在垂直方向基本规律是河谷平坝区少雨，山区多雨。流域径流主要由降水补给，时空变化规律同降水基本一致。6～10月径流占年总量的80%左右；11月至翌年5月径流量占年总量的20%左右。
安宁河流域内洪水主要由暴雨形成，大洪水主要出现在6～9月。安宁河流域洪水多发生于局部，总的趋势是冕宁县泸沽镇以上和孙水河流域洪水影响泸沽镇至西昌、德昌一带；而德昌至河口段的洪涝灾害主要受西昌及安宁河下段的茨达河，锦川河影响。
3 旱涝分析
3.1 用河道来水（径流深）距平百分率划分旱涝等级
河道来水距平百分率指河道来水距平值与多年同期平均值的百分比值。计算公式为：
3.2 安宁河流域旱涝分析
（1）从表1、表2可以看出，安宁河中上段和下段旱涝等级并不完全一致。差别最大为1981年，安宁河中上段为枯水年份，而下段为正常年份。58年间完全一致的有29年。
（2）安宁河中上段和下段发生洪水等级也不完全一致。最典型的是1960年、1981年和2012年。1960年孙水河流域发生流量1400m3/s（60年一遇）的洪水，而安宁河下段仅发生最大流量为1480m3/s（2年一遇）的洪水。1981年孙水河流域最大洪水流量仅为270m3/s，而安宁河下段最大洪水流量为3410m3/s（30年一遇）。2012年孙水河流域发生流量2510m3/s（200年一遇）的特大洪水，而安宁河下段仅发生最大流量为2330m3/s（7年一遇）的洪水。
（3）安宁河中上段1954年发生重涝后到1968年再出现重涝14年，1968年到1974年6年，1974年到1998年24年。1975年发生重旱后到1981年再出现重旱6年，1981年到2006年25年，2006年到2011年5年。而安宁河下段1965年发生重涝后到1968年再出现重涝3年，1968年到1974年6年，1974年到1998年24年。1957年至2014年58年间安宁河下段2011年只出现一次重旱。
（4）安宁河下段丰枯变化具有明显的阶段性特征，年丰枯序列大致可划分为3个以丰为主的阶段和4个以枯为主的阶段。1957年～1964年共8年为枯水期，1965年～1974年共10年为丰水期，1975年～1982年共8年为枯水期，1983年～1991年共9年为丰水期，1992年～1997年共6年为枯水期，1998年～2004年共7年为丰水期，2005年～2011年共7年为枯水期。2012年进入丰水期。安宁河下段枯水期持续6～8年，丰水期持续7～10年，丰枯大致以8年为转换周期。
（5）安宁河下段大多枯水期中有偏丰年，丰水期中有偏枯年。1965年～1974年10年中有3年丰水年，1年偏丰年，3年正常年，3年偏枯年；1998年～2004年7年中有3年偏丰年，1年丰水年，3年正常年，明显偏涝。而1992年～1997年6年中有4年为偏枯年，2年正常年，明显偏旱。58年间有几次旱涝急转或交替转换。1967年、1968年、1969年偏枯年急转丰水年再急转偏枯年；1974年、1975年丰水年急转偏枯年；1997年、1998年偏枯年急转丰水年；2011年、2012年枯水年急转偏丰年。
4、径流深预测
现引用安宁河出口控制站湾滩水文站1957年至2014年58年年径流深实测数据对安宁河下段旱涝周期进行预测研究。
先把年径流深实测数据按梯度排列，然后进行回归计算，从而建立自回归预报模型。经过率定，用历年年径流深进行12位梯度排列后进行回归计算得到的预报方程，预测下一年年径流深精度较高。预报方案精度达72.2%，精度检验达80.0%。此方案中许可误差按《GBT 22482-2008 水文情报预报规范》７.５.４.a）规定执行。即：许可误差为多年同期实测变幅的20%：（976-418）×20%=112mm。此预报方案编制率定期为1969年～2004年共36年，检验期为2005年～2014年共10年。具体排列方法见表3。由于篇幅有限表中只列出了头和尾及不合格点数据。根据径流序列 变量前后的相依性，建立自回归模型，表达式为：
5、结语
（1）通过对安宁河旱涝特征分析，得出该区域丰枯转换周期大约为7～8年，14～16年一个轮回。大多枯水期中有偏丰年，丰水期中有偏枯年。
（2）上述方法建立的预报模型完全采用实测数据，没有任何假定，而且只用了径流深一个因子。预报模型对安宁河下段未来径流的估算，虽然有一定误差，但是从一定程度上可反映未来气候变化对安宁河下段径流丰枯变化的影响。
（3）预报方法重点在于数据排列方式。此方法旱涝急转年份时误差较大。年初做当年丰枯预测时请综合分析。笔者用此方法对安宁河湾滩水文站及雅砻江泸宁水文站年最大洪水进行回归计算（许可误差以实测值的20%为准），预报方案精度达70%以上，精度检验也达70%以上。希望从事相关专业的人士深入研究。
（4）从预测值可以确定安宁河下段2012年～2018年为丰水期，2019年～2026年为枯水期。与旱涝特征分析结论一致。
（5）综合分析得出如下结论。安宁河全流域2012年～2018年为丰水期，2019年～2026年为枯水期。安宁河全流域2018年发生大洪水的可能性较大。而安宁河流域2019年、2025年发生大旱或重旱的可能性较大。
参考文献：
[1]钟佩辰.安宁河流域水文特征初步分析[J].城市建设理论研究，2014，（24）.
[2]陈涛，王文彬，陈壮志.四川省孙水河流域旱涝周期分析[J].四川水利，2014，（4）.
[3]王文彬.四川省孙水河流域枯水期丰水期预测[J].地球，2014，（3）：123-124.
[4]王殿武，迟道才，梁凤国，江行久.区域旱涝特征分析及灾害预测技术研究[M].北京：中国水利水电出版社，2011.5：42-50，60-61，79-87.
[5]水利部水文局，长江水利委员会水文局.水文情报预报技术手册[M].北京：中国水利水电出版社，2010.10.
[6]GBT22482-2008，水文情报预报规范[S].
作者简介：
赵洪 男，1965年生，汉族，工程师，四川省西昌水文水资源勘测局。从事水文情报及水资源研究工作。
论文作者:赵洪,吴永键,王文彬,钟佩辰
论文发表刊物:《电力设备》2017年第8期
论文发表时间:2017/7/19
标签：宁河论文; 流域论文; 径流论文; 旱涝论文; 枯水期论文; 偏枯论文; 洪水论文; 《电力设备》2017年第8期论文;