基于贝叶斯推理的环境水力学反问题研究

来源 :浙江大学 | 被引量 : 0次 | 上传用户:dalu008
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着环境水力学正问题模型(预测模型)求解的日益成熟,其有关的反问题的研究越来越受到了的重视。环境水力学反问题解的主要困难在于其不适定性,而这种不适定性(尤其是不唯一性)主要来源于水环境系统的不确定性。水环境系统中含有许多不确定性因素,如测量数据的不确定性、数学模型的不确定性和模型参数的不确定性等,如何正视并解决这些不确定性是环境水力学反问题研究亟需解决的问题。本文从概率论的角度出发采用贝叶斯统计学的方法-贝叶斯推理建立了环境水力学反问题的求解模型,其中模型参数、测量数据、先验信息和最终反问题解的形式等全部采用概率语言描述,这样就较好地解决了环境水力学反问题求解中的不确定性问题。理论研究方面,本文对贝叶斯推理中标准的抽样技术-马尔科夫链蒙特卡罗抽样方法(Markov chain Monte Carlo,MCMC)进行了研究。为了提高马尔科夫链对后验参数空间(解空间)的搜索能力,本文在Metropolis-Hastings算法的核心上提出了一种动态多链的搜索方法(Dynamic Multi-chains Metropolis-Hastings,DMMH)。该算法由于采用多条链的搜索,因而增加了对后验空间的搜索能力,克服了单链搜索能力不足的问题;同时DMMH中马尔科夫链可以根据实时的搜索状态进行调整(减少)链的个数,因而比常规的多链搜索节省了计算时间。环境水力学正问题求解的基础数学模型为对流扩散方程,本文主要采用有限体积法和有限元法作为求解方法。在有限体积法方面,特别研究了基于通量限制器(flux limiter)的有限体积法,该方法具有精度高、稳定性好的优点,比较适用于需要调用大量正演计算的环境水力学反问题的求解。本文提出了一种改进的通量限制器格式M7,数值计算表明M7通量限制器与传统的一些限制器相比,提高了计算精度和显式差分计算的稳定性。应用研究方面,本文采用贝叶斯推理主要求解了两类反问题,即环境水力学的参数估计反问题和污染物源项识别反问题。环境水力学数学模型丰富,为了验证基于贝叶斯推理的反演模型的合理性和可靠性,本文根据模型分类的不同,研究了一维模型、二维模型、稳态模型、非稳态模型、常系数模型、变系数模型、含源模型、非含源模型等各种类型的环境水力学数学模型,大量算例计算表明在正问题求解精度高(数值误差低)的情况下,采用贝叶斯推理和MCMC抽样方法获得的反问题的解具有信息量大(能给出环境水力学参数的后验分布)且估计精度高的优点。在反演结果的可靠性和估计的精度方面,本文重点探讨了测点位置、测量异常值和测点数目对反演结果的影响,研究了反问题求解的非唯一性问题。源项识别反问题主要研究了单个污染源和多个污染源下的位置和强度的识别问题,采用MCMC的贝叶斯推理均能获得较高的识别精度。最优化方法目前是环境水力学反问题求解的主流方法,它和贝叶斯方法两者之间既有区别也有联系。本文也介绍了基于优化方法(遗传算法)的环境水力学反问题求解方法,并提出了基于混合遗传算法结合有限体积法(FVM-HGA)的参数识别方法。该方法与传统遗传算法相比具有较强的全局寻优和局部搜索能力,且反演的精度高、收敛速度快;但其缺点是不能像贝叶斯方法那样获得参数的后验分布。由于MCMC本身也可作为参数识别的工具,本文也对基于MCMC的环境水力学参数识别问题进行了研究,提出了基于贝叶斯理论的以似然函数为收敛准则的参数识别新方法(L-MCMC)。算例表明,L-MCMC在减少反演计算量上效果显著。总而言之,对于诸如环境水力学反问题等一些环境系统反问题,基于MCMC的贝叶斯方法是一种强大的求解工具,既能给出模型参数的分布规律,又能作为优化算法给出模型参数的最优估计,具有重要的研究价值和应用前景。
其他文献
【正】一、看质量证明标志农机产品大多获得国家有关技术鉴定部门颁发的合格证书和证章。如“农业机械推广许可证”、“生产许可证”和“3C(强制性产品认证)标志”等,这样的
长江源区位于青藏高原腹地,具备“高、寒、旱”的特点,其生态系统脆弱。近年来,随着全球气候变化,源区出现了河道断流、冻土退化、植被破坏等一系列的生态环境问题,引起社会各界的
当温家宝在人大闭幕后对英国金融时报记者谈到“今年恐怕是中国经济最困难的一年”时,相信很多人对这样沉重的定调都没有心理准备。那么,作为中国经济大背景下的建材行业,是不是
--、西方诗歌的译介催促了中国新诗的诞生“五四”新文学的诞生,正如鲁迅所说:“一方面是由于社会的要求的,一方面则是受了西方文学的影响。”(鲁迅:《<草鞋脚>小引》,《鲁迅
回 回 产卜爹仇贱回——回 日E回。”。回祖 一回“。回干 肉果幻中 N_。NH lP7-ewwe--一”$ MN。W;- __._——————》 砧叫]们羽 制作:陈恬’#陈川个美食 Back to yield
【正】 把打碎的旧电池掺入积肥坑中(1平方米1节1号电池即可),土粪中含有的有机物在微生物的作用下生成腐殖酸,腐殖酸与旧电池中残锌起化学反应后生成腐殖酸盐。另外,旧电池
目的 观察醒智散对血管性痴呆大鼠学习记忆能力的改善作用,并研究其可能的机制.方法 通过反复夹闭双侧颈总动脉同时腹腔注射硝普钠的方法制备血管性痴呆大鼠模型.然后将大鼠
由冲积河流主槽和河漫滩组成的复式河槽在调节洪水、削减洪峰、贮存泥沙等方面不同于单一河槽。复式河槽随着流量增大、水位抬高,当水位超过滩地高程时洪水漫滩,产生一系列水
<正>距离动用631.7亿元的巨额资本买下万达集团13个文旅项目和76个酒店项目只有月余,"地产狂人"孙宏斌挂帅的融创中国日前又大口吃进了广东广弘房地产有限公司50%的股权,同时
初春是呼吸道感染病多发季节,受气候影响容易引发某些呼吸道传染病的局部流行。2003年春季的“非典”流行,使人们对春季的流行病提心吊胆。鉴于此,本刊特推出一组“初春季节