DREAM算法分析地下水数值模拟不确定性的影响因素.PDFVIP

第６２卷　 第２期　　 地　质　论　评 　　　　 ＧＥＯＬＯＧＩＣＡＬＲＥＶＩＥＷ　　 Ｖｏｌ６２　Ｎｏ２ ２０１６年 ３月 Ｍａｒ．　２０１６ ＤＲＥＡＭ算法分析地下水数值模拟 不确定性的影响因素 １，２） １） ３） １） 杨运 ，吴吉春 ，骆乾坤 ，吴剑锋 １）南京大学地球科学与工程学院水科学系，南京，２１００９３； ２）水利部淮河水利委员会，安徽蚌埠，２３３００１；　３）合肥工业大学资源与环境工程学院，合肥，２３０００９ 内容提要：地下水数值模拟常受到模型不确定性、观测资料不确定性等多种不确定性因素的影响，对这些影响 因素进行定量分析十分必要。将差分进化自适应梅特罗波利斯（ＤｉｆｆｅｒｅｎｔｉａｌＥｖｏｌｕｔｉｏｎＡｄａｐｔｉｖｅＭｅｔｒｏｐｏｌｉｓ，ＤＲＥＡＭ） 算法与ＭＯＤＦＬＯＷ模型结合应用于地下水数值模拟不确定性的定量分析。以模型结构概化、水位观测资料误差这 两种重要不确定性来源为例，通过一个理想地下水流模型，系统分析两者对模型参数及模型输出结果不确定性的影 响。研究结果表明：模型结构概化及水位观测资料误差共同引起了地下水数值模拟的不确定性，但模型结构概化起 到了主控作用。模型结构概化合理时，模型参数及模型输出结果的不确定性较小，并随观测资料误差不确定性的增 大而增大；模型结构概化不合理时，模型参数及模型输出结果主要受控于模型结构概化带来的影响，且不确定性显 著增大；观测资料误差相同情况下，模型结构概化越接近于真实的水文地质条件，模型参数及模型输出结果的不确 定性越小。 关键词：ＤＲＥＡＭ；模型结构概化；观测资料误差；数值模拟；不确定性；地下水 　　地下水数值模拟已成为刻画地下水流和溶质运 即对复杂的含水层参数分布进行分区简化，与真实 移行为的重要手段。然而地下水系统复杂的特性加 的水文地质条件就可能产生差异。而另一方面，在 之人类认知的局限性，使得地下水模型和实际的地 搜集、测量、记录水位观测数据过程中，由于人为或 下水系统之间往往存在着差异，这种差异一般归结 系统原因，不可避免又会带来水位观测资料的误差。 为地下水模拟的不确定性（姚足金等，２００９；吴吉春 模型结构概化与水位观测资料误差对地下水模拟不 等， 确定性的影响不容忽视。 ２０１１）。一般而言，地下水模拟不确定性按照来 源可以分为参数不确定性、模型不确定性和观测资 差分进化 自适应梅特罗波利斯（Ｄｉｆｆｅｒｅｎｔｉａｌ 料不确定性（Ｈａｓｓａｎｅｔａｌ．，２００８）。这些不确定性因 ＥｖｏｌｕｔｉｏｎＡｄａｐｔｉｖｅＭｅｔｒｏｐｏｌｉｓ，ＤＲＥＡＭ）算法（Ｖｒｕｇｔ 素的存在严重影响了地下水模拟结果的可靠性，因 ｅｔａｌ．，２００９）是近年来提出的具有较高求解效率的 此，对其进行定量分析十分必要。 马尔科夫链蒙特卡洛采样方法之一（Ｙｕｓｔｒｅｓｅｔａｌ．， 近十多年来，关于地下水数值模拟不确定性的 ２０１２），已被成功用于多种不确定性分析问题（曹飞 研究多集中于参数不确定性的影响研究（Ｂｅｖｅｎｅｔ 凤，２０１０；ＬａｌｏｙａｎｄＶｒｕｇｔ，２０１２；ＳｈｉＸｉａｏｑｉｎｇｅｔａｌ．， ａｌ．，２００１；Ｈａｓｓａｎｅｔａｌ．，２００９；ＬｕＤａｎｅｔａｌ．，２０１４； ２０１４）。本文将其与ＭＯＤＦＬＯＷ模型结合应用于定 Ｌｉａｎｇｅｔａｌ．，２０１５），如渗透系数、孔隙度等参数的空 量分析地下水数值模拟的不确定性。以模型结构概 间变异性对地下水模拟不确定性带来了不同程度的 化、水位观测资料误差这两种重要不确定性来源为 影响（陈彦等，２００５；施小清等，２０１２；骆乾坤等， 例，通过一个理想的地下水流模型，系统探讨两者对