基于LSTM实时校正的WRF_WRF-Hydro耦合径流预报.pdf

年 月水利学报 第 卷第 期文章编号：（）收稿日期：；网络首发日期：网络首发地址：?基金项目：国家自然科学基金项目（）；长江生态环境保护修复联合研究（第二期）（）作者简介：刘昱辰（），博士生，主要从事数值大气模拟、流域水文预报研究。：通信作者：刘佳（），正高级工程师，主要从事气象水文耦合模拟预报、雷达遥感数据同化研究。：基于 实时校正的?耦合径流预报刘昱辰，刘佳，刘录三，李传哲，王瑜（中国环境科学研究院，北京 ；中国水利水电科学研究院 流域水循环模拟与调控国家重点实验室，北京 ）摘要：为改进?陆气耦合系统的径流预报效果，减小耦合系统在峰现时间、洪峰流量预报上的误差，本文在使用变分数据同化技术充分降低预报降雨误差水平的基础上，采用长短期记忆人工神经网络 对?耦合系统的径流预报过程开展了实时校正研究，并与自回归滑动平均模型 实时校正结果进行对比。研究结果表明，通过数据同化技术可有效提升 模式降雨预报精度，降低 模式的输入误差，但径流预报准确性仍有待提升。对比 和 两种实时校正模型对耦合径流预报结果的实时校正：在前 预见期，两种模型在中国北方半湿润、半干旱地区山区小流域 场典型洪水预报中的表现基本接近，除场次 外，和 两种模型在 预见期的衰减速率分别为 和 ，随着预见期的延长，径流预报精度的衰减速度在整体上慢于 模型，预报效果优于 模型。关键词：；实时校正；?耦合系统；径流预报；数据同化中图分类号：文献标识码：?研究背景在气候变化的大背景下，极端降雨过程会更加频繁和强烈地发生，气象预报与水文预报的结合，有助于探索大气过程对陆面水文过程的影响，提高降雨和径流的模拟、预报能力，是水文预报发展的必然趋势 。相比其他陆气耦合系统，中尺度数值大气模式 （）与其陆面水文模式 的耦合可以实现大气与陆面水文过程的单向与双向交互，模拟水热通量在大气 陆面界面的运动与转化过程，其应用前景十分广泛 。但受到模型结构体系、初始状态和输入信息误差、参数辨识能力等限制，该耦合系统的不确定性较为显著且难以量化 ，模拟结果难以直接应用于实际预报，对洪水预警的指导意义有限。由于降雨和径流形成的物理机制较为复杂，现有各类陆气耦合预报系统在模拟和预报精度方面往往都会存在一定的偏差。研究表明 ，对陆气耦合系统的径流预报结果进行实时校正，能够有效提高耦合系统的降雨 径流预报能力。随着计算机技术的发展，人工智能与大数据技术发展迅速，其中，作为人工智能核心技术之一的人工神经网络，广泛应用于水文预报的实时校正领域 。带有记忆单元的长短时记忆神经网络（，）的提出，使时间序列上的记忆信息变得可控，也使 在模拟和预报具有时间序列属性的过程方面发展的愈发成熟。等 研究表明，将 结合 可以改进径流预报结果，有效提高 在径流预报方面的潜力。等 建立了基于 的网格地表产流模型，结果证明 在我国北方山区小流域地表径流预报中可表现出良好的性能。在 应用于水文预报领域的研究中，大多是采用 直接对径流等特征量进行预测 ，缺乏与具有物理机制的水文模型相结合。在传统误差实时校正模型中，自回归滑动平均模型（，）是一种较为常用的函数模型，相对人工智能误差校正方法，模型参数简单，不易受环境变量的影响。这种方法在水文预报领域已被普遍认可，得到了广泛的应用研究 。为充分提高耦合系统的径流预报能力，本研究在使用 模式数据同化改进预报降雨的同时，针对耦合系统仍然存在的径流预报偏差，采用 和 模型分别对预报径流进行实时校正。首先，通过同化雷达反射率与传统气象监测数据 （），改进 模式的降雨预报精度，在此基础上构建出?耦合系统；其次，配合 模型，实现预报径流的实时校正，提高陆气耦合系统的综合预报能力；最后，与传统实时校正模型 结果进行对比验证。研究成果可为中小尺度流域的洪水风险管理与预报预警提供技术支撑。研究区域概况本研究选取位于大清河流域南支的阜平流域和北支的紫荆关流域作为研究区，流域面积分别为 和 ，研究区地处太行山东麓，全境为山地，地势高低落差明显，陡峭的地形、裸露的地表和复杂的产汇流机制，极易造成严重的洪水灾害。流域面积研究区地理位置、地形地势以及雨量站的分布情况见图 。图 研究区位置图根据阜平流域和紫荆关流域 年以来发生的历史洪水情况，依据同流域相关研究 ，按照不同的降雨类型及径流特征，选取了 场典型降雨过程开展陆气耦合径流预报研究。选取的降雨场次在气象成因、强度、历时、落区上均具有一定的代表性，其分别具有雨强小且持续时间较长、雨强较小但引发了明显的流量涨落、极短时间内雨强大引发洪水陡涨陡落等多种不同特点。为方便分析和比较，所选取 场降雨的降雨时长均为 ，其降雨历时等详细信息见表 。其中，场次 在统计时段降雨发生前，已发生过较多降雨场次，自 年 月 日至 月 日，阜平流域发生的累积面雨量达 ，并且在 月 日至 月 日，出现过一个降雨强度大且持续时间长的洪水过程，其中洪峰流量达到了?。因此，在场次 发生时，流域内存在较大基流，且土壤湿度较为饱和，极易引发洪水，月 日的最高洪峰流量达到?。表 研究流域 场典型降雨场次信息降雨场次所在流域降雨历时 面雨量?最大流量?（?）阜平?阜平?阜平?阜平?紫荆关?紫荆关?降雨场次的降雨强度、落区、历时的不同，会产生不同的径流形成结果。对流域而言，降雨时空分布的均匀程度是影响流域出口断面流量过程的最重要因素。为定量的描述降雨过程的时空分布情况，采用变差系数 ，分别从空间和时间尺度对典型降雨过程的均匀程度进行评价，并将其作为划分降雨类型的指标：珋()槡（）式中：计算空间分布均匀性时，为第 个雨量站的累积降雨量，珋 为 个雨量站的累积降雨量的平均值，为雨量站的数量；计算时间分布均匀性时，为流域第 个小时的面雨量，珋 为 个小时面雨量的平均值，为降雨持续时间总小时数。通过设置 临界值划分降雨类型，降雨空间分布均匀性的临界值 为 ，降雨时间分布均匀性的临界值 为 ，当降雨场次的 值小于 临界值时，判断降雨在时间或空间上是相对均匀的 ，。各降雨场次空间分布和时间分布的 值及降雨类型见表 。表 场降雨的 值及降雨类型降雨场次空间分布 空间分布是否均匀是是否是否否时间分布 时间分布是否均匀是是否否否否 实验方案研究区域的地面观测数据包括来自阜平流域 个、紫荆关流域 个雨量站的实测小时降雨量和阜平水文站、紫荆关水文站的实测小时径流量，空间分布见图 ，均由河北省水文水资源勘测局提供。本研究针对雷达反射率和 数据进行同化。雷达数据采用了位于石家庄市的 波段多普勒雷达系统，该系统的初始数据来自中国气象局。雷达扫描范围见图 ，覆盖半径为 ，空间分辨率为 。是传统的地面与高空气象数据的集合，由世界气象组织（）提供，本研究通过程序仅截取覆盖研究区的 数据作为同化资料，并只针对该区域的数据转化为可被 同化的格式。（）是由国家环境预报中心提供的实时数据，用于驱动 模式的初始场和边界场，空间分辨率为 ，每隔 更新一次。?耦合预报系统构建考虑到两个研究流域位于山区，降雨的不确定性更为突出，模式预报结果仍会受到空间误差的影响 。为了获取可靠的预报降雨，前期利用天气雷达对云雨系统的探测能力和获取瞬时降雨信息的优势，结合传统气象监测数据，基于 三维变分数据同化系统，开展了数据同化频率对 模式降雨预报的影响研究。前期成果 为本文研究提供了可靠的降雨输入。为使设置的嵌套网格与雷达数据的空间分辨率尽可能相近，并提高优化效率，采用三层嵌套网格的方式逐层增加区域的水平分辨率，最外层网格范围的设置利于同化大尺度的 数据，使模式运行更加稳定。最内层网格完全覆盖研究区，相邻嵌套层的嵌套比例设为 ，网格尺寸分别为 、和 。全球预报系统为预报提供了初始和横向边界条件。由于 没有经过同化分析处理，因此在 模式中常作为预报降雨的初始驱动 ，流域的嵌套网格和雷达观测范围见图 。图 阜平流域和紫荆关流域的网格嵌套方案表 模式主要物理参数化方案序号物理参数化方案选择微物理过程 长波辐射 （）短波辐射 陆面方案 边界层 （）积云对流方案 （）图 基于数据同化的?陆气耦合预报系统示意参数化方案的选择借鉴了 等 在相同研究流域的参数化方案敏感性试验结果。对同化试验中使用的降雨产生有重大影响的参数化细节见表 。其中，和 网格的积云对流方案被关闭，模式仅在 最内层嵌套网格（即 网格）上运行，并采用 分辨率的次网格进行水文过程的模拟。陆面模式方案选择 ，驱动数据输入步长和输出数据的时间间隔均为 。研究发现，在两个研究流域，数据同化对 模式预报降雨量及其时空分布的改进明显：数据同化前，模式的预报降雨量远低于实际观测，数据同化后，各小时的预报降雨量均有所增加，同化雷达反射率结合 传统气象监测数据对 预报降雨有明显的改进效果，且 同化时间间隔和 分辨率的降雨输出可以使 模式更好地预报降雨的时间和空间分布情况。因此，在本研究中，的输入由逐小时循环同化后的 预报降雨提供，基于数据同化的陆气耦合流程见图 。径流预报偏差校正方法 实时校正模型 是一种特殊的循环神经网络，其允许先前输入的信息保存在网络的内部，从而影响网络的输出，相比于传统的循环神经网络，在隐藏层中加入了记忆细胞单元，可以较好地记住较长的历史信息。图 模型结构示意图 表 设置参数表参数名称变量设置输出层个数 隐含层神经元数 学习率 学习率衰减因子 批处理量 时间步大小 迭代更新次数 隐含层激活函数 记忆细胞单元可以对输入的数据进行选择性的记忆和遗忘，通过引入遗忘门、输入门、输出门，使神经网络长期传输记忆信息，结构见图 。在本研究中，模型建模在 下运行，研究所使用的 及其变体模块，均采用 （版本 ）框架中的 模块设计。根据文献 ，当流域较小或可获取的数据量不充足时，采用单层隐含层可获得更好的训练效率及运行结果，因而本研究选择了 层作为 的隐含层数。通过试算得出，隐藏神经元个数 时可以使模型均方误差最小。模型训练前，还需要确定模型的学习率、批处理量参数和选代次数等参数。其中，学习率决定了模型权重更新的幅度，但目前还缺少确定最佳学习率的最优方法。在本研究依据经验将学习率设定为 ，同时将学习率衰减因子设定为 ，使学习率随着迭代次数增加缓慢变小，则模型收缩时震荡能得到缓解。批处理量参数即一次输入模型中的样本个数，参照相关研究 ，批处理样本数与隐含层神经元数相同时处理效果较好，故将批处理量定为 。是参数选代更新次数，根据本研究中的数据量大小，将 定为 。主要参数设置见表 。在本研究中，基于 的逐小时径流预报构建了逐小时预报径流误差的历史序列，并将其作为 输入变量，以 为时间步长预报不同预见期下的径流误差，最后将预测的误差序列叠加至径流预报结果中，完成耦合系统水文预报的实时校正。实时校正模型 模型由自回归模型（模型）与移动平均模型（模型）为基础“混合”构成，是可以反映水文变量在时序变化上统计特性的数学模型。本研究针对研究流域内实测径流量和 预报径流量之间的逐小时误差序列 （），通过构建具有时间特征的误差序列，建立基于误差的自回归方程。设 时刻的预报值为 （），实测值为（），则 （）（）（）即为误差序列。假定误差序列具有一定的自相关性，则可以利用过去的误差预报未来的误差。（，）的结构取决于阶数 和 。在本研究流域，已有相关研究表明当 ，时可以取得较好的校正效果 ，因此，本研究选用 （，）模型对耦合系统的预报径流过程进行实时校正。评价指标利用纳什效率系数（）、均方根误差（）、洪峰流量的相对误差（）、洪量的相对误差（）以及峰现时间绝对误差（）对 模拟的径流量进行评估。对于评价指标，用以验证水文模式模拟结果的可信程度。用来衡量预报值同观测值之间的偏差。为洪峰的相对误差，反映模拟洪峰的可信程度，为洪量的相对误差，反映模拟洪量的可信程度。为实际峰现时间和预报峰现时间的时差，用于表示峰现时间预报的误差。（）（珔）（）（）槡（）（）?（）（）?（）（）式中：为洪水过程中第 小时的径流预报值；为洪水过程中第 小时的径流实测值；珔为 的平均值；为时间长度；为场次洪水洪峰流量的预报值；为场次洪水洪峰流量的实测值；为场次洪水洪量的预报值；为场次洪水洪量的实测值；为实际峰现时间；为预报峰现时间。结果分析 校正前的?耦合预报结果同化后预报降雨驱动下的紫荆关流域和阜平流域的径流预报结果如图 所示。为更好的显示降雨引发的流量起涨过程，在 雨洪过程的基础上，对实测径流时间进行了延长。图 同化后预报降雨驱动下的?耦合系统径流预报结果总体上，通过 同化系统，可以使 模式预报降雨量增大，降雨预报整体与实际降雨量更为接近，其预报径流的起涨过程与实际流量过程线也较为吻合。同时，模式对量级小、陡涨陡落、退水较快的洪水过程有更好的模拟效果，但对量级较大且退水较慢的洪水过程流量过程线吻合程度较差。表 给出了 径流预报结果。表 同化前后预报降雨动下的 径流预报结果场次空间 时间 同化前的预报流量同化后的预报流量 洪峰误差?洪量误差?峰现误差 洪峰误差?洪量误差?峰现误差 从评价指标结果可以看出，径流预报精度主要取决于 模式的降雨预报误差，与驱动数据质量和降雨分布类型密切相关。如强对流降雨场次 ，降雨在空间和时间尺度的变差系数 值分别达到了 和 ，短时集中的强降雨过程导致其预报结果存在较大误差，不理想的降雨预报间接影响了洪水的预报精度。虽然数据同化技术为水文过程的模拟与预报提供了更精确的降雨输入，但径流与降雨并非简单的线性关系，模型结构的复杂程度和参数率定过程都会对径流预报产生影响，多因素导致 在径流预报结果上仍存在一定的不确定性。实时校正径流预报结果采用 模型配合?耦合预报系统，对 场径流的预报结果进行实时校正，预见期为 、和 的径流预报结果见图 和表 。图 校正后的?径流预报结果由图 可以看到，模型对?耦合预报系统径流预报的实时校正具有较好的改进效果，尤其是对耦合系统在 预见期的预报结果，改进最为明显。预见期超过 后，模型对耦合系统径流预报的改进程度总体呈下降趋势，预报径流过程的波动增大，但洪峰流量和洪量没有明显的增加或降低。分析其原因，可能是由于随着预见期的延长，预报径流误差明显增加，径流序列数据之间的关联性降低。模型此时已较难准确学习到时间序列的数据特征，因而无法对预报径流序列进行有效的校正，具体体现在径流预报的波动也就越大。表 校正下?耦合系统的径流预报指标场次预见期?场次预见期?场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 将 模型引入?耦合预报系统之后（评价指标见表 ），场次洪水的预报精度有了明显的提高。整体来看，场雨洪过程 预见期的 在 之间，平均 为 。随着预见期的延长，逐渐降低。当预见期小于 时，预报精度下降较为缓慢，但仍有较好的预报校正效果；预见期大于 时，预报精度下降较快，其中表现最好的雨洪过程（场次 ）在 预见期的 仍可达到 。从洪峰相对误差指标来看，预见期的洪峰预报精度最高，径流过程的洪峰相对误差位于区间范围为 ，预见期大多数场次洪峰相对误差小于 。但由于预报流量过程线随着预见期的延长波动较大，洪峰流量和洪量在较长预见期内不呈现明显的规律。实时校正径流预报结果以数据同化后的预报降雨驱动 ，采用 模型对场典型雨洪过程的径流预报结果进行实时校正，预见期为 、和 的径流预报结果见图 和表 。图 校正后的?径流预报结果经 模型校正后，在 预见期，耦合系统能准确预报出所有径流过程的发展趋势，在 之间。预见期的预报径流比较稳定，平均 为 ，但 预见期后的 开始出现较明显的下降。除预见期外，的实时校正效果与校正前的径流预报结果也有一定关系。例如，校正后场次 在 预见期的径流预报 为 ，预见期的径流预报 仍能达到 ，实时校正结果十分理想。对比之下场次 的预报结果和校正效果则相对较差，预报精度衰减速度也相对较快。校正前，场次 的径流预报 为 ，洪峰误差为 ，场次 的 为 ，洪峰误差为 。由于校正前场次 的径流预报精度较低，作为误差序列的实时校正模型，是根据以往的历史径流误差和扰动来推测未来时刻的误差值，所以当误差序列相关性较差时，误差校正效果也不明显。实时校正结果对比将 与 模型校正下的耦合系统径流预报结果进行深度对比，进一步分析 对?陆气耦合系统的径流预报改进效果。图 给出了校正前后 场典型雨洪过程不同预见期的 指标变化情况。由图 可以看出，与校正前的径流预报结果相比，实时校正模型使预报径流的精度明显提高。根据不同预见期与 预见期的 指标变化百分比计算预报径流的 指标衰减速度，可以发现，校正下耦合系统预报径流的衰减速度明显慢于 模型。例如，场次 在 模型校正下，预见期径流预报的衰减速度为 ，预见期的衰减速度为 ，而 模型校正下预见期的衰减速度为 ，预见期的衰减速度为 。除场次 以外，其他场次的径流预报都显示了相似的衰减规律。这与校正前耦合系统的径流预报结果有关，误差序列之间的相关性越强，径流误差的校正效果则更为显著。表 校正下?耦合系统的径流预报指标场次预见期?场次预见期?场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 场次 校正前 预报值 预报值 预报值 预报值 预报值 预报值 讨论与结论 讨论为了获得较高的水文预报精度，通常采用实时校正方法，针对历史预报误差的相关性对预报径流进行实时校正，但该方法受预见期的影响较大。就本研究流域而言，由于流域汇流时间较短，实时校正后的预报径流过程线出现了不同程度的波动（如图 、图 ），且波动随着预见期的延长而增大。可见，与 模型对中小流域在 以上的预见期下预报径流的实时校正能力有限。考虑到大尺度流域汇水时间较长，流量起涨平稳缓慢，预报径流误差即使在 或 预见期以上也具有较高的相关性，与 模型的校正潜力则可以体现在更长预见期的径流预报。结论为进一步提高?耦合系统的径流预报能力，在数据同化预报降雨的同时，针对耦合系统仍然存在的径流预报偏差，研究采用 模型对预报径流进行实时校正，并与 模型的径流校正结果进行了对比，结论如下：（）（，）与 都能进一步改进?耦合系统的径流预报能力，尤其在预见期对预报径流的改进效果最好。校正后，研究区 场雨洪场次的 预见期 区间分别为 （）、（），但随着预见期的延长，两个模型对径流预报的改进程度逐渐降低。（）随着预见期的延长，模型对耦合系统径流预报的整体改进效果优于 模型，耦合系统预报精度的衰减速度慢于 模型。当预见期达到 时，在研究区大部分雨洪场次的 指标在 以上。图 和 对 场典型雨洪过程实时校正后的 指标对比图本研究采用 构建了预报径流的误差校正模型，实现了对?耦合系统预报径流的实时校正。研究还存在一些不足，如 模型在预见期超出 时的预报精度仍有待提高等。在进一步的研究中，将采用 与小波变换或 神经网络等方法结合，利用径流序列的频率域误差特征，结合多种校正方法，开展多模型联合优化方法，以更好地提高陆气耦合预报系统的准确性。参考文献：，?（），（）：，（）：，：，：，（?）：孙明坤，李致家，刘志雨，等 模型与新安江模型在陈河流域的应用对比 湖泊科学，（）：顾天威，陈耀登，高玉芳，等基于 模式的清江流域洪水模拟研究 水文，（）：，郭生练，钟逸轩，吴旭树，等水库洪水概率预报和汛期运行水位动态控制 中国防汛抗旱，（）：雷晓辉，王浩，廖卫红，等变化环境下气象水文预报研究进展 水利学报，（）：熊怡，周建中，孙娜，等基于自适应变分模态分解和长短期记忆网络的月径流预报 水利学报，（）：，（）：，：，（）：姜淞川，陆建忠，陈晓玲，等基于 网络鄱阳湖抚河流域径流模拟研究 华中师范大学学报（自然科学版），（）：，（），（）：崔巍，顾冉浩，陈奔月，等 与 神经网络在福建小流域水文预报中的应用对比 人民珠江，（）：焦伟杰，龙海峰基于自回归模型的分布式水文模型预报校正 水资源与水工程学报，（）：田济扬天气雷达多源数据同化支持下的陆气耦合水文预报 北京：中国水利水电科学研究院，谢平，霍竞群，桑燕芳，等基于 模型的水文序列相依变异分级方法及验证 水利学报，（）：，（）：，（）：王维基于?陆气双向耦合的混合产流区水文模拟研究 北京：中国水利水电科学研究院，：，（）：，（）：，：刘艳基于深度学习的水文要素模拟与预测研究 天津：天津大学，王菲菲基于长短期记忆神经网络的城市降雨径流预测研究 重庆：重庆大学，徐乾顺基于机器学习的中长期径流预报模型研究及系统集成 武汉：华中科技大学，（）：，（）：，（）：?，（，；，）：?，?（），（），：；?；（责任编辑：韩昆）（上接第 页），（，）：，（），（）：（），（），：（），；（），：；（责任编辑：于福亮）