Page 125 - 2023年第54卷第12期

此外，图 ４（ｄ）—（ｆ）三幅图分别为采用 ＣＥＥＭＤＡＮ算法基于华县站 ５月份径流序列计算的时不变阈
值序列、移植阈值序列和时变阈值序列及其干旱识别结果。以时变阈值序列识别结果为基准（图 ４（ｆ）），
３
６０年中有 ２４年处于干旱状态（４０％干旱发生率），累积径流亏缺量为 ３９３．９７ｍ ?ｓ。而基于时不变阈值
序列（图 ４（ｄ））识 别 出的 干旱 频次相 比时 变阈值 序列偏 少约 ２５．００％，累 积 径 流 亏 缺 量 存 在 低 估 约
２９．４２％，如 ２００３—２０２０年均未出现干旱状态，与实际旱涝情况不符，表明该阈值选取不够合理。基
于移植阈值序列（图 ４（ｅ））识别出的干旱频次相比时变阈值序列偏多约 ３３．３３％，累积径流亏缺量存在
高估约 ８７．９１％，如 １９８６—２００２年一直处于干旱状态，也与实际旱涝情况不符，该阈值选取也不够合
理。而时变阈值序列（图 ４（ｆ））识别的干旱状态频次较为适中，累积径流亏缺量也相对正常，能够合理
地刻画出非平稳条件下基准期和变化期的实际旱涝情况。
图 ５分别是标准化径流指数法（图 ５（ａ）—（ｃ））和阈值法（图 ５（ｄ）—（ｆ））对应 ３个方案对华县站旱
涝情况的表征。其中，标准化径流指数序列服从均值为 ０，方差为 １的标准化正态分布，其累积时间
序列（图 ５（ａ）—（ｃ）中右侧纵坐标轴所示）最终应返回 ０值，且上升表示变涝，下降表示变旱。比较结
果表明：时不变标准化径流指数（图 ５（ａ））在基准期数值全部大于 － １．５，对该时期的水文干旱事件烈
度存在低估，其累积 ＳＳＩ － １２ 曲线在基准期呈现整体上升趋势（即持续变涝），趋势斜率为 ０．３４，在
ｉｎｖａｒ
变化期呈现整体下降趋势（即持续变旱），趋势斜率为 － ０．５５，说明此方案计算的标准化指数难以合理
地刻画基准期径流过程的旱涝转化。时不变阈值法（图 ５（ｄ））也同样存在对基准期水文干旱事件径流
亏缺量低估的情况。移植参数标准化径流指数（图 ５（ｂ））对基准期的水文干旱事件表征较为合理，但
变化期出现大量数值小于－ ２．０ ，对该时期的水文干旱事件烈度存在高估，其累积 ＳＳＩ － １２ 曲线在变化
ｔｒａｎ
期持续下降（即持续变旱），趋势斜率为 － １．１４ ，且最终远远偏离 ０值，说明此方案难以合理地刻画变
化期径流过程的旱涝转化过程，移植阈值法（图 ５（ｅ））也存在相同问题。时变标准化径流指数（图 ５（ｃ））
则可以合理地识别基准期和变化期的水文干旱事件，两个时期的 ＳＳＩ － １２数值均有小于 － ２．０的时刻
ｖａｒ

图 ５ 华县站三种设计方案下 ＳＳＩ － １２指数和 ＴＬ － １２阈值方法对旱涝情况的表征结果
５
— １ １ ３ —