基于改进蒙特卡洛算法的河流水质概率区间时空变化预测研究.pdfVIP

2017年第5期 水利技术监督 理论研究 DOI：10．3969／j．issn．1008—1305．2017．05．024 基于改进蒙特卡洛算法的河流水质 概率区间时空变化预测研究 曹一梅，胡嘉佳，李显鸿 (云南省水文水资源局昆明分局，云南昆明650051) 摘要：文章引入交叉熵抽样方法对蒙特卡洛随机抽样进行优化和改进，并将改进算法结合水动力一水质模型对云 南中部某河流水质概率区间时空变化进行预测。结果表明：改进算法在河流主要水质指标趋势的预测精度上好于 累积概率最高达O．646；改进算法可估算不同概率下各水质指标随时闽变化的波动曲线，可以河流水质预警提供 关键依据；研究河流中上游出现高浓度氨氮和低浓度溶解氧的概率高于下游。研究成果对于河流水质预警系统构 建提供重要依据。 关键词：交叉熵抽样方法；改进蒙特卡洛算法；水质概率区间时空预测；云南中部河流 中图分类号：TV213．4 文献标识码：A 文章编号：1008—1305(2017)05—0072一03 水质概率的预测可以为区域水质预警提供重 计算状态均通过方程(1)进行确定。 要的依据，水质预警可对区域水质恶化进行有效 在确定不同状态变量的失效概率后，需要各变 控制，近些年来，有许多学者展开区域河流水质 量参数风险值进行计算，计算方程为： Ⅳ 1 概率的研究，取得一定研究成果¨巧1，在这些成 (2) 果中蒙特卡洛算法可对水质数据进行随机抽样计 式中，y(，)一参数风险值；Js一计算仿真状态变量 算，在河流水质概率应用较为广泛。但是传统蒙 参数；x一随机状态变量；E(F)一均方差值。 特卡洛算法存在仿真计算和抽样时段过长的局限， 在参数风险值确定后，需要对算法的收敛精度 使得随机抽样计算效率较低。当前，有学者引入 进行判定分析，判定方程为： 交叉熵抽样方法对对传统算法进行过优化和改进， 并在水资源领域得到相关应用n引，但是在河流 ㈩ p=掣=瑶铲 水质概率预测中应用较少，特别是对河流沿程水 式中，卢一模型计算收敛均差系数；¨E(F)]一期 质概率预测还未得到具体应用，为此本文引入改 望最小均差值，经变化得到判定方程为： 进的蒙特卡洛算法，以云南红河上游流域某河流 Ⅳ=_』j圭鲁 ㈩、(4)’ 为实例，对该河流水质概率区间进行时空变化 肚嵩备[口E(F)]2 式中，在收敛精度确定的条件下，可以减少随机变 预测。 量的抽样均方差来优化抽样次数，从而提高模型收 1 改进的蒙特卡洛算法原理 敛计算精度。 本文引入交叉熵抽样方法对模型随机抽