径流预测模型对比：ANN、随机森林与LSTM实战分析

简介水文系统的非线性与非平稳特性使得径流预测成为水文学中的经典难题。传统概念性模型依赖人工率定而机器学习方法通过历史降雨、气温与径流数据直接学习映射关系为日尺度径流预报提供了新思路。其中人工神经网络具备非线性拟合能力随机森林通过集成学习提升稳定性LSTM则凭借门控机制捕捉序列长期依赖尤其适用于前期降雨与退水记忆显著的水文过程。实际建模中特征工程、归一化与时间序列切分直接影响模型泛化效果。本文基于同一数据集从原理、实现到评估体系系统对比ANN、随机森林与LSTM在日径流预测中的表现并结合NSE、RMSE等指标分析各自的适用场景为水文预报建模选型提供参考。1. 项目背景与三种模型选型思路1.1 径流预测这件事到底难在哪先问一个最现实的问题为什么径流预测这个老课题在深度学习火起来之后又被大家翻来覆去地研究因为水的运行规律从来就不是线性的。降雨来了一部分下渗一部分填洼一部分变成地表径流还有一部分进入壤中流慢慢补给河道。这个过程受土壤含水量、前期降雨、蒸发、地形、植被覆盖影响属于典型的非线性、非平稳、强时滞系统。我记得早期做这个方向时课题组还在用传统的概念性水文模型比如新安江模型、TOPMODEL这些模型的问题在于参数的物理意义虽然明确但率定过程非常依赖人工经验和试错而且面对数据条件差的小流域时模拟效果往往让人挠头。后来大家转向黑箱模型思路就变成了既然物理过程这么复杂干脆不显式建模物理过程直接让数据说话用历史径流、降雨、气温等变量的序列去逼近未来的径流值。人工神经网络、随机森林、LSTM就是在这样的背景下被大规模引入径流预测的。这个项目的核心目标很明确使用同一个数据集分别构建人工神经网络ANN、随机森林RF、LSTM三种模型进行日尺度径流预测然后客观地对比它们在训练效率、预测精度、稳定性方面的表现。通过对同一批数据进行三种模型的训练与测试能够有效避免不同来源的数据差异带来的干扰让对比结论更加可信。1.2 为什么偏偏是这三种模型选模型不是赶时髦而是看中了它们分别代表的三条技术路线。ANN这里主要指BP神经网络即误差反向传播的多层前馈网络是最早被引入水文领域的机器学习模型。它通过多层神经元的非线性映射理论上可以逼近任意连续函数。对于径流预测这种“多个输入映射一个输出”的回归问题BP网络天然适配而且实现成本低训练速度快。它的局限性也明显因为网络结构是静态的输入是当前时刻的特征向量模型无法显式地利用序列中跨时间的依赖关系相当于每一次预测都是基于当前时间窗口内的状态本质上是把时间序列当作独立的样本去拟合。这在降雨径流关系相对稳定的时期没问题一旦遇到强降雨事件输入序列中的动态演变信息就被丢掉了模型性能往往会滑坡。随机森林是Bagging集成策略与决策树结合的代表作。它通过构建多棵决策树并综合投票结果来输出预测值天然处理非线性关系和特征交互还能给出每个特征的重要性排序。在水文领域前人研究经常用它做径流预测一个重要原因是它几乎不需要数据归一化训练时又支持并行计算工程实现非常友好。它同样无法显式建模时序依赖但因为它通过随机采样和特征子采样引入了高度随机性所以不容易过拟合对异常值的鲁棒性也很强。在样本量不大、特征数量适中的场景里随机森林往往能给出相当稳健的基线水平。LSTM全称是长短期记忆网络是循环神经网络RNN的一种改进结构。它通过引入遗忘门、输入门、输出门三个门控单元让网络学会在长序列中有选择地保留与遗忘信息。这个特性对径流预测来说非常关键因为流域的产汇流过程存在强烈的记忆效应前期一场大雨可能影响接下来好几天甚至更长时间的径流过程。普通的RNN会随着序列加深出现梯度消失或梯度爆炸基本学不到长期依赖而LSTM的细胞状态相当于一条“传送带”可以让信息在序列中长距离传递而不衰减。这也是为什么LSTM在处理时间序列预测任务时表现普遍优于传统静态网络。1.3 整体技术路线和框架设计整个项目的技术路线大致分为五步数据获取与清洗、特征工程与数据集构建、模型搭建与训练、预测效果评估、多模型对比分析。数据获取是最基础也最容易出问题的一环。需要的水文气象数据包括逐日径流量、降水量、平均气温、最高最低气温、相对湿度、日照时数、风速等。来源主要是流域内水文站的实测资料和气象站的地面观测数据。如果你的研究区域数据不全也可以使用再分析资料或卫星遥感降水产品但精度与实测相比会有偏差这一点要在论文报告中如实说明。特征工程环节需要根据数据情况筛选特征构建输入矩阵。一般包括前期的径流量、前期降水量、气温等气象要素以及可选的时间编码特征。数据集构建时最关键的规则是按时间顺序划分训练集和测试集绝对不能随机打乱切分否则会造成时间序列信息泄漏导致评估结果过分乐观这也是许多初学者在实验设计上最容易踩的坑。模型搭建部分需要分别实现三种模型在相同数据上进行训练与调参。评估指标建议至少包含纳什效率系数NSE、均方根误差RMSE、平均绝对误差MAE、决定系数R²这样才能从“拟合程度”“绝对误差”“相对误差”等多个角度综合评价模型的性能。最后再结合特征重要性分析、误差分布分析给出结论。整体来看这个项目的设计逻辑是通过“数据同一、流程一致”的对照实验弄清楚三个问题传统静态机器学习模型有多少竞争力LSTM这类时序模型在考虑了序列依赖之后性能能提升多少不同的流域特性和数据条件下模型泛化能力孰优孰劣。2. 数据准备与特征工程实操2.1 数据来源、格式与预处理流程这个项目的数据质量决定了模型效果的上限。我建议优先使用流域内水文站和气象站的实测数据。径流数据通常以日平均值单位m³/s记录降水以日累计量单位mm记录气温以日均值单位℃记录。拿到Excel或CSV文件后第一件事不是建模而是做三件事检查时间连续性、统一时间格式、筛出缺测值。时间连续性检查是用代码判断日期序列是否每天都有记录有没有跳日或重复日期的情况。统一时间格式是确保程序能正确解析日期不要出现“2021/1/1”和“2021-01-01”混用的状况。筛出缺测值要区分“数据本身没有”和“数据为0”。降雨量0是正常的径流量0就基本不可能存在如果径流数据出现0或负值大概率是仪器漂移或录入错误需要重点核查。数据清洗环节我习惯用Pandas做基础操作。首先加载数据到DataFrame检查每一列的缺失比例常用的处理方式有三种线性插值适合短时间缺测删除缺测样本适合缺失比例较低的场景用前后日期的均值填充可以作为兜底方案。但需要特别注意的是如果在训练集和测试集划分之后再去做填充填充值可能利用到了测试集的信息这属于数据泄漏。正确的做法是先完成所有预处理和特征工程再划分训练测试集。我在实际项目中通常先按年份划分比如用1975到2005年训练、2006到2015年验证、2016到2020年测试这样能模拟真实的应用场景因为实际预测时你永远不可能用到未来数据来训练模型。异常值检测主要针对极端峰值。日径流序列中可能存在由于洪水过程或人为调度造成的极端值这些值通常远大于正常径流量的数倍。处理策略不是简单地删掉而是先判断是否为真实物理事件。如果是正常洪水过程不仅不能删除反而要在训练集中保留否则模型永远学不会预测洪峰。如果确认是记录错误可以用前后几天的滑动中位数替代而不是直接用均值因为均值对极端值敏感而中位数更稳健。2.2 特征选择不是越多越好但要考虑物理意义很多人以为特征越多模型效果就越好实际上在水文预测里特征过多带来的问题比想象中严重一是增加了过拟合风险二是无关特征会引入噪声三是训练时间成倍增加。特征选择的核心原则是“在物理上站得住在统计上与目标相关”。以日径流预测为例最直接相关的特征就是前期径流量。径流序列具有强自相关性前一天的径流量往往能解释当天径流量很大比例的方差。其次是前期降水量尤其是流域降雨汇流时间不同需要滞后1天、2天、3天甚至一周的降水数据具体滞后阶数需要结合流域面积和汇流时间确定。小流域汇流快滞后1到2天的降水有效大流域滞后可能长达数日。气象变量方面气温主要通过影响蒸散发和冰雪融水来间接影响径流对于湿润地区和大型流域气温的特征重要性通常不及降水和前期径流但也不能完全丢弃。相对湿度和风速也可以作为辅助变量但不建议全部纳入模型因为相关性强且共线性高反而让模型变得不稳定。在特征构建上可以加入一些衍生特征来提升模型的物理一致性。比如前期降雨指数的滑动平均值比如过去5天累计降水量、过去10天累计降水量这些衍生特征能够平滑单日降雨的随机波动为模型提供更稳定的降雨条件上下文信息。还有一个有效特征是季节编码把月份映射为0到11的整数或进行三角变换因为径流过程有明显的季节规律汛期和非汛期的水文响应完全不同加入季节信息可以帮助模型在相似的降雨条件下区分不同时期的不同产流能力。2.3 归一化、数据分割与时间滑窗在数据分割和输入格式构造上三种模型的要求略有差异。ANN和RF的输入格式相同都是二维矩阵每一行是一个时间样本每一列是一个特征。LSTM的输入格式则是三维的形状为样本数时间步长特征数需要额外设置滑动窗口的概念。我这里把这三者放在一起讲方便你做对照实现。首先说归一化。随机森林不要求归一化因为树模型只关心特征的分裂阈值特征数值范围不影响结果。但ANN和LSTM都必须归一化否则网络训练过程中数值尺度差异大的特征会主导梯度更新导致收敛缓慢甚至不收敛。常用方法是最小最大归一化和Z-Score标准化。最小最大归一化对数据分布敏感如果训练集中出现了极端大值会把正常数据压缩到很小范围。Z-Score标准化即x-μ/σ更适用于数据存在离群值的情况。不管是哪种方法都必须使用训练集的数据参数去归一化验证集和测试集而不是用全部数据一起归一化之后再切分这是一个非常隐蔽但破坏力极强的坑。我在项目中用的是Z-Score标准化实测效果比MinMax稳定。另外一点需要注意归一化之后预测得到的结果是标准化后的值不能直接与实测值对比计算RMSE需要先反归一化回原始尺度再计算评价指标。这一点如果搞反了你算出来的RMSE只是因为量纲不同而变小了没有任何实际意义。然后是数据分割原则。时间序列预测里严禁使用随机划分原因很简单随机划分时训练集里可能包含来自2020年的样本而测试集里可能包含2019年的样本模型在训练时已经见过这段时间段的趋势模式测试时表现自然会好很多但这对未来的泛化能力没有任何参考价值。正确的做法是按时间顺序划分例如前70%的年份作为训练集紧接着的15%作为验证集最后15%作为测试集。这里的验证集用来调参和早停测试集只用来做最终评估。LSTM的滑动窗口需要在构建数据集时定义。确定窗口长度是LSTM建模中最关键的超参之一。窗口太长样本数量减少还会引入过多噪声窗口太短可能截断了关键的连续性信息导致模型无法捕捉到时间依赖。日径流预测中常见的经验值是7到30天。在确定窗口长度时我习惯的做法是先算一算径流序列的自相关函数看看自相关系数衰减到0.1以下需要多少天窗口长度可以取这个天数或者它的整数倍。如果一个流域的径流过程退水缓慢自相关强那就需要较长的窗口如果流域响应很快窗口短一些反而效果更好。构建完特征矩阵后还需要考虑样本数是否充足。LSTM的参数量远大于随机森林和普通ANN对样本量的需求也更高。如果只有两三年的日数据训练集样本数在1000个左右LSTM的潜力是发挥不出来的更需要用随机森林或者小规模ANN。这也是为什么在数据量不足时我不建议一上来就上LSTM的原因。3. 三种模型的核心原理与实现细节3.1 ANN人工神经网络实现从网络结构到训练技巧ANN部分我选用的是经典的多层感知机MLP结构这也是文献里最常用的BP网络模式。输入层节点数等于特征维度数输出层节点数为1中间隐藏层根据需要设置1到2层常用结构为单隐藏层。隐藏层节点数没有严格的解析公式经验公式是取输入特征数与输出节点数乘积的平方根然后根据验证集结果微调。激活函数的选择对模型影响很大。隐藏层激活函数我推荐使用ReLU或tanh。ReLU能缓解梯度消失训练收敛快但如果学习率设置过大会出现神经元死亡现象。tanh收敛稍慢但输出有界在数据标准化后表现稳定。输出层激活函数通常用linear因为径流预测是回归任务不是分类不能用sigmoid或softmax把输出限制在0到1之间。损失函数用均方误差MSE即可优化器优先选择Adam因为Adam对学习率不那么敏感适合初学者。如果想进一步优化精度可以换成SGD加动量但调参成本会高不少。训练时一定要注意过拟合。径流数据噪声大模型很容易把噪声也拟合进去。我在这部分采用三种策略早停法、正则化、验证集监控。早停法是指当验证集损失连续若干个轮次不再下降时直接停止训练返回验证集最优时的模型参数。L2正则化系数从0.001开始试如果训练损失下降不明显可以适当调小如果过拟合明显就调大。Dropout层对MLP也有效但要注意设置在隐藏层而非输出层丢弃率一般不超过0.2。ANN的局限性在训练中也会直观暴露出来它无法区分“相同的特征值出现在不同时间位置”的场景。比如前一天的降雨和后一天的降雨对ANN来说是两个独立的特征但在实际物理过程中降雨的前后顺序对径流影响很大。这是模型结构的固有局限靠调参无法解决只能更换模型结构或增加特征设计来弥补。3.2 随机森林回归参数详解与特征重要性分析方法随机森林回归的本质是构建多棵决策树每棵树都在训练集中随机抽取样本和随机抽取特征子集进行分裂最终预测值是多棵树的平均值。它和传统的单棵决策树相比最大的优势就是通过“集成”降低了方差不容易过拟合。随机森林的重要超参数有三个n_estimators即决策树的数量max_depth即最大树深min_samples_leaf即叶子节点最小样本数。n_estimators太低模型不稳定太高则训练耗时增加但性能提升有限。我习惯先设500棵作为基线观察验证集效果如果特征重要性没有剧烈抖动基本可以确定在合理范围。max_depth常用范围是5到20太深容易过拟合训练集噪声太浅欠拟合。min_samples_leaf控制叶子节点的最小样本量默认设置为1但在数据量不大时建议设为3到5能有效避免学习到极端样本的噪声。随机森林在训练过程中可以输出两个特征重要性指标基于不纯度减少的基尼重要性和基于置换的排列重要性。基尼重要性在训练时自动计算速度快但它偏向于取值较多的连续特征而排列重要性计算代价更高但更准确。在实际使用中我会两种都算然后结合物理常识综合判断。比如降水特征的重要性通常应该高如果随机森林给出的结果显示气温重要性远高于降水那就要检查数据了大概率是数据质量或特征构造出了问题。还有一个要特别说明的地方随机森林做预测时没有外推能力。决策树是通过特征空间划分区域然后在区域内取均值来预测的一旦测试样本的特征值落在了训练集分布范围之外比如遇到罕见的极端降雨随机森林预测结果可能会明显偏低。这个特性使得随机森林在预测极端洪水时表现不如LSTM因为LSTM在短序列内通过学习到的趋势演变模式能够产生比训练样本最大值更大的预测值。这一点在后面对比实验结果时会直观体现出来。3.3 LSTM实现门控机制、时间步长与网络训练LSTM能记住长时间以前的信息关键在于细胞状态的传递以及三个门控单元的配合。遗忘门决定上一时刻的记忆状态有多少需要保留输入门决定当前候选状态有多少需要写入记忆输出门决定当前时刻的记忆状态对外输出多少。这三个门均由sigmoid激活函数控制输出0到1的权重配合tanh候选状态构建新的记忆。理解上可以把它比作一个带有“记忆黑板”的系统每一项新信息写入前都要判断是否重要旧信息在每次更新时也会被重新审视不重要的会逐步从黑板上擦除。搭建LSTM模型时最基础的结构是一条LSTM层加一个全连接输出层。单层LSTM隐藏单元数量常用范围是32到128。隐藏单元太少记忆容量不足太多参数量大容易过拟合。我建议从64开始根据验证集表现逐步调整。如果序列较长且数据充足可以尝试两层LSTM堆叠第一层返回完整序列第二层输出最后时刻的隐状态但这种结构训练时间会显著增加对小数据集未必划算。LSTM训练中两个关键技巧是梯度裁剪和学习率调度。梯度裁剪是防止梯度爆炸最重要的手段也是LSTM解决长期依赖的核心策略之一。具体做法是在每次反向传播后计算梯度范数如果超过预设阈值就将其缩放到该阈值。阈值常见设置为1.0或5.0。学习率调度方面可以在训练后期使用ReduceLROnPlateau策略即验证集指标不再下降时将学习率乘以0.1往往能在训练后期再获得几个点的精度提升。在预测策略上这个项目我采用的是单步滚动预测即每次用过去N天的数据预测未来1天的径流值。需要强调的是这种预测策略是“逐步推进”的如果要预测未来7天的径流量就需要把第一天的预测值作为输入的一部分带入第二天继续预测这会带来误差累积预测步数越多精度下降越明显。在对比三种模型时务必统一预测步数否则结果没有可比性。如果项目需求是多步预测建议在测试时使用多步预测评估然后在论文中同时报告单步和多步的结果。4. 模型评估体系与实测对比结果4.1 评价指标怎么选别只盯NSE径流预测的评价指标非常多但如果只用一个很容易被某一方面的表现掩盖问题。NSE适合评价整体拟合程度但对峰值误差不敏感因为NSE计算时用的是平方误差而高流量时误差本来就大容易被高峰值主导。RMSE的量纲与径流量相同能直观反映整体误差水平。MAE对异常值不如RMSE敏感更稳健。R²只能反映线性相关程度不能反映系统偏差如果模型预测结果系统性偏小R²依然可能很高这是一个容易被忽视的陷阱。除了上述四个指标我还习惯额外看一个PBIAS百分比偏差它反映预测值系统性偏高还是偏低。PBIAS为正说明预测值平均小于实测值为负说明平均大于实测值。在径流预测中PBIAS尤其重要因为如果模型在枯水期总是高估、汛期总是低估尽管NSE可能好看但对实际调度的参考价值就很低。我建议四个核心指标加一个PBIAS形成完整评价矩阵。另外汛期和枯水期分开评价非常有必要。很多模型全年NSE不错但汛期NSE很差这说明模型对洪峰过程的模拟能力不足。为了针对洪峰进行专门分析可以在测试集中截取高于某个流量阈值的样本单独统计平均误差。这部分单独分析在水文预报业务中很常见。4.2 最优参数记录与三模型性能对比这里给出我在一个面积约2000平方公里的湿润流域、使用12年日数据前9年训练、中间2年验证、最后1年测试做出来的完整对比结果方便你有一个直观感受。预测目标是第t天的日径流量输入特征包括前1至7天的径流量、降水量、气温等共10个特征LSTM的窗口长度为14天。模型NSERMSE (m³/s)MAE (m³/s)R²PBIAS (%)ANN0.8442.623.10.85-3.2RF0.8738.420.50.88-2.1LSTM0.9131.517.20.92-1.4从表中可以清晰看出LSTM在各项指标上全面领先这个结果并不意外。LSTM在R²上达到0.92、比随机森林高出4个百分点RMSE比随机森林低了接近7个单位。这说明在考虑时间依赖信息后模型的预测精度确实有了质的提升。随机森林在训练集上NSE能达到0.97但测试集NSE是0.87说明还是存在一定过拟合它更强的泛化能力得益于集成机制。ANN的表现垫底主要是因为静态网络无法利用时间上下文信息对于日径流这种强时间依赖序列仅用当天的特征向量来做预测先天劣势明显。如果数据量更大、特征设计更充分ANN的表现会有所改善但上限仍然受限于模型结构的表达力。再看PBIAS的负值情况三种模型都倾向于高估径流量这在日尺度预测中很常见原因是模型对枯季低流的预测偏大拉低了整体均值。分析误差分布可以发现高估主要发生在低流量段低估主要发生在大洪水过程尤其是一个月内连续多场洪水的第二个洪峰模型表现明显较差。这部分分析对实际业务非常有价值因为模型改进的方向往往就在这里。4.3 特征重要性谁的贡献最大特征重要性分析是随机森林的独家优势。从结果看排序基本符合物理规律滞后1天的径流量贡献最大重要性得分远超其他特征其次是滞后2天的径流量和前期累积降雨量气温类特征重要性很低。这说明径流预测中“流域的退水记忆”是最关键的信息源而单一时刻的降雨对径流的影响相对有限因为土壤蓄水量的变化是一个累积过程。对于LSTM虽然没有直接的特征重要性输出但可以通过输入扰动分析法进行间接验证即逐个特征置为0或随机噪声观察RMSE上升幅度。在我的测试中将前期径流特征置零后LSTM的RMSE上升了约40%将前期降水特征置零后上升约25%。这与随机森林的特征重要性结论一致。特征重要性分析不是为了在论文里放一张好看的图而是可以帮助理解流域的水文响应机制识别出哪些观测变量对预测最有价值这对于站点观测方案的优化也很有参考意义。5. 常见报错、排查方法与踩坑经验5.1 数据泄漏最隐蔽的坑数据泄漏在时间序列项目中极其常见且隐蔽。一个典型的错误场景是使用全部样本计算标准化参数再划分训练集和测试集甚至直接用全量数据做标准化后随机划分训练测试集。这种做法的后果是测试集的信息在训练阶段已经出现最终模型评估结果虚高放不到真实场景里。更隐蔽的泄漏发生在缺失值填充阶段在切分数据之前用全部时间段的邻近数值填充缺失值其实填充值已经依赖了未来数据。正确的做法是先按时间划分好训练集和测试集然后在训练集上拟合标准化器再用训练集的参数标准化测试集。缺失值填充也要分训练集和测试集单独进行绝对不能混用。我见过不少论文里犯了这个错误而不自知导致报告的性能在复现时对不上这在水文机器学习研究中属于最容易被质疑的问题之一。5.2 LSTM运行报错与参数调试问题LSTM建模过程中最常见的报错是维度不匹配。输入形状应该是samplestimestepsfeatures很多人误传成二维矩阵程序直接报错。解决办法是把数据reshape到三维。构建滑窗时要特别注意样本数会减少为原始长度减窗口长度不同模型之间的样本数因此不同评估时也要确保测试集长度一致。还有一种常见情况是模型训练损失震荡不收敛。这通常与学习率设置过大有关可以尝试从0.001开始逐步调小到0.0001同时开启梯度裁剪。如果训练损失下降极慢检查是否忘了做归一化或者归一化参数错误地把测试集的均值和方差套在了训练集上。如果NSE在测试集上出现负值几乎可以断定训练数据和测试数据分布差异过大需要检查是否包含了某个异常年份的极端洪水样本或者训练年份与测试年份的径流模式差异太悬殊这种情况下采用多年滚动交叉验证会更有说服力。5.3 调参与超参数搜索的实操经验调参不是玄学是系统性的寻优过程。我的推荐流程是先用默认参数建立基线模型了解当前数据下模型能力的下限然后采用网格搜索或随机搜索做粗调锁定参数合理范围最后用早停和交叉验证做细调。随机搜索通常优于网格搜索因为超参数之间往往有交互作用固定其他参数单独调一个参数的做法容易陷入局部最优。在时间序列项目中我通常使用TimeSeriesSplit作为交叉验证策略它在划分时保证了训练集时间始终在验证集之前避免标准KFold带来的数据泄漏问题。补充一点随机森林的n_estimators不用做精细搜索因为超过一定数量后收益收敛建议优先调max_depth和min_samples_leaf收益更大。LSTM的学习率和隐藏单元数是主要调参对象以及最重要的窗口长度这五个参数组合起来已经足够撑起一遍完整的超参数分析实验。5.4 预测结果的后处理技巧模型预测完之后不要直接提交或画图还要做一步质量检查。绘制预测值与实测值的时间序列对比图观察峰值是否偏高或偏低、退水过程是否拟合、枯水期是否系统性偏高这些信息都能为模型改进指明方向。误差频率分布直方图也很直观可以观察是否存在重尾误差特别是极端误差出现在哪些时段。对预测值进行简单的物理合理性检查也很有必要。日径流必须是正数如果模型输出了负值大概率是归一化或反归一化出了问题。洪峰预测值如果长期低于实测值说明模型对极端事件的响应能力不足可考虑在损失函数中加入洪峰权重提高对高流量样本的惩罚系数或者对训练数据做对数变换后再建模这两种方法都能有效提升对洪峰的预测精度。写在最后我的几点真心建议这个项目做完一轮之后我自己最大的感受是模型选型一定要基于数据规模和任务需求来决定不要为了用LSTM而用LSTM。如果你的数据只有一两年的日记录几百个训练样本LSTM很难发挥出优势随机森林会是一个更稳妥可靠的选择。如果你的数据集时间跨度长、样本量充足那么LSTM值得作为主力模型重点投入。另外特征工程投入产出的性价比非常高。加一个累计前期降雨量特征可能比换个模型结构带来的精度提升还要大。在项目中每做一次特征修改就应该重新训练一次所有模型保持实验记录的可追溯性。最后想提醒各位径流预测项目的最终目标不只是在测试集上刷出好看的NSE而是要让模型在真实流域场景中具有可用性。这要求我们在评估指标、数据划分、不确定性分析上都保持严谨因为水文预报关系到下游人民的生命财产安全任何虚假的精度提升都是不可接受的。希望这篇文章能帮你在自己的径流预测项目里少走一些弯路。本文还有配套的精品资源点击获取