地下水模拟不确定性分析

引言：最近在做一个二维剖面地下水模型，奈何当地地质、水文地质资料少的可怜，在建模型时发现模型概化及后续的校正都存在很大的问题，具有强烈的不确定性，不知如何下手。就想着能不能通过敏感性分析多分析下，通过分析多种可能性，尝试得出结果。然后就查找不确定性的论文，突然就发现了这么一个新大陆，原来不确定性分析还有这么多道道，原谅我之前只知道参数的敏感性分析[笑哭]。
在此将昨天看的几篇文献简单整理下，一是笔记二是供大家参考学习。不当之处还请指正。

什么是地下水模拟不确定性？

地下水数值模拟被广泛用于多个方面，如水资源评价、场地污染、水资源管理和规划等。但是存在多种因素影响地下水模型的准确性，如气候变化和人类活动对地下水系统的影响、自然水文过程的复杂性、理论技术水平以及资料缺乏等，这些因素导致地下水模型和实际的地下水系统之间始终存在差异，造成模拟结果和实际情况的偏差，影响模型的可靠性。这种差异即为地下水模拟的不确定性。

地下水模拟不确定性分类？

按不确定来源可分为：

1. 参数不确定性
2. 模型不确定性
3. 资料不确定性

参数不确定性

参数不确定性是指模型输入参数，如渗透系数、储水系数、弥散度、孔隙度、化学反应速率等水文地质参数和化学参数，具有空间变异性，这使得地下水模拟具有不确定性。如渗透系数随深度逐渐减小，是普遍性的非均质性；再如弥散度具有明显的尺度效应。
对于这种不确定性是目前研究相对较多的。几乎大部分地下水数值模型研究都会有参数敏感性分析的内容。敏感性分析即属不确定分析方法中的一种。

模型不确定性

模型不确定性又可分为概念模型和数值模型不确定性。

概念模型不确定性

概念模型不确定性是指因缺乏水文地质条件或对水文地质条件理解不够透彻，使所建概念模型不能准确反映真实情况，在此类模型基础上所得出的结果就有很大的不确定性。概念模型是地下水数值模拟的基础，因其不确定性所造成结果的不确定性或不可靠性要远大于其他两种不确定性所造成的，且易被忽视，因此应当引起重视。
概念模型不确定性主要表现在：

1. 模型边界位置和类型不准确
2. 地质单元厚度不准确
3. 参数分区不够精细
4. 对含水层类型判断错误（如浅层微承压当做潜水）
5. 对水力联系判断错误等

针对概念模型不确定分析一般采用多模型模拟的方法，即先确定一个理想或假设真实的基本实例模拟，然后改变其中的一个或者几个边界条件（如隔水变为定水头）、或者改变一下模型结构（如忽略弱透水层、承压水当做潜水），设置2-3个对比模型，然后评估不同对比模型对模拟结果的影响。
推荐2篇模型不确定性分析实例：

1.吴吉春,陆乐.地下水模拟不确定性分析[J].南京大学学报(自然科学版),2011,47(03):227-234.

简介：利用贝叶斯方法定量分析了模型边界条件偏差对模拟结果不确定性的影响。

2.曾献奎,王栋,吴吉春.地下水流概念模型的不确定性分析[J].南京大学学报(自然科学版),2012,48(06):746-752.
简介：基于多模型模拟方法利用蒙特卡洛模拟分析了不同模型结构对模型输出的影响，并研究了模型结构导致的不确定被模型参数不确定所弥补或抵消的可能。

数值模型不确定性

数值模型不确定性是在概念上和求解过程中造成的误差。如现有的对流弥散方程式建立在典型单元体（REV）概念基础上的，无法在足够细的尺度上刻画介质的非均质性，因此难以适应非均质介质中的溶质运移问题。此外，大多数溶质运移模型的数值解法都很难完全避免截断误差、解振荡和数值弥散问题。

资料不确定性

最后一种模拟不确定性为观测资料不确定性。
地下水模拟过程中需用到大量的观测资料，如初始水头（溶质）分布、用于模型校正的水头（溶质）观测序列、抽注水井和泉流量、降雨量和蒸发量等。当观测资料时空分布较稀疏时，就需要插值进行推姑；另一方面，观测资料采集整理过程中也可能产生误差，这些观测资料中的误差就造成了地下水模拟的不确定性。

不确定分析的方法

根据不确定性给分析方法求解原理可分为：

1. 蒙特卡洛法
2. 矩方程法
3. 贝叶斯法
4. 条件模拟法
5. 敏感度分析法
6. 一次二阶矩法等

具体的分析方法有大量的资料，在此不做介绍。

参考资料：

1. 吴吉春,陆乐.地下水模拟不确定性分析[J].南京大学学报(自然科学版),2011,47(03):227-234.
2. 曾献奎,王栋,吴吉春.地下水流概念模型的不确定性分析[J].南京大学学报(自然科学版),2012,48(06):746-752.
3. 薛禹群, 谢春红. 地下水数值模拟[M]. 北京: 科学出版社, 2007.