水文地质学中的反演问题.docVIP

地下水反演是一个综合性问题，研究重点是对含水层特征进行刻画。许多建模人员已经对概念模型的不确定性（显著的时间和空间变化）、尺度依赖性、许多类型的未知参数（导水系数、补给和边界条件等）、非线性以及对含水层性质经常表现出低灵敏性的状态变量（特别是水头和浓度）进行了研究。由于这些难题的存在，因此不可能把模型校准与模拟过程分割开来，而是应当将其作为认识含水层特征的重要一步。事实上，真实参数的估计方法尽管在计算细节上会有一些差异，但在本质上没有什么区别。就是否还有充足的空间能够改进地下水反演问题存在争议：开发用户友好的代码、通过地质统计学方法调整变异性、将地质学信息和其它类型的数据（温度、同位素浓度和地下水年龄等）相结合、正确计算不确定性等。尽管如此，即使利用现有的代码进行自动校准，也仍然极大地推动了模拟工作的发展。因此，主张将这一方法的应用进行标准化。 一、概　述 广义上，反演模拟是指利用模拟参数的测量值来收集模型信息的一个过程，其中包括两个相关的概念：模型识别和参数估计，在本文中，后者与校准具有相同的含义。“模型识别”是应用一些方法来寻找模型的特征，如控制方程、边界条件、时间范围或非均质性等。而“参数估计”是对那些刻画特征的性质进行赋值。 地下水模拟中的“反演”与上述定义大致相同。事实上，一些水文地质学反演问题已经建立了标准规程，例如，水文地质学家采用抽水试验来确定含水层的特征。解释这些试验结果的过程包括：首先确定最合适的模型（根据实际的测量值和地质条件），然后根据曲线拟合来估计模型参数，利用自动程序可以很容易完成这项工作。在解析解很少的领域中，如非饱和带，自动化的校准程序已经标准化了（Kool和Parker，1988；Hollenbeck和Jensen，1998）。在通过适当的测试刻画空间变化的领域中，如医学（扫描或X线断层摄影术）（Rudin等，1999；Rhli等，2002）和地球物理学（Sambridge和Mosegaard，2000）等等，自校准已经成为标准。然而，地下水模拟显示出了能够区别于其它模拟的几个特性： ???? 成本 地下水模拟费用相对较高。模型每次运行时，都要求解大型方程组（以获取系统准确而真实的特征）。另外，计算机科学的发展可以大大提高模型的性能； ???? 时间依赖性 状态变量，如水位和浓度具有时间依赖性。许多流动问题本质上都处于稳定流条件下，但即使如此，仍然会包含水位随时间波动的信息，这样就需要进行非稳定流的模拟； ???? 非均质性 水力传导系数K通常是最重要的水力学指标，其值会相差几个数量级；导水系数也是如此，相当于二维模型中的K（在下文中，这两个参数将会交替使用）； ???? 不同类型的参数。当重点考虑导水系数时，其它参数（补给、边界流量等）同样具有不确定性和相关性； ???? 尺度依赖性 野外测量的参数通常只能表示含水层的一小部分，这样，在性质和数量上就与模型所需的参数具有很大的差异； ???? 模型不确定性 含水层的几何特征和非均质性受地质条件的控制，而地质条件很难得到确切的认识。 ???? 灵敏度低 在不同的问题中，状态变量可能会对模拟参数表现出很低的灵敏度（即信息量不足）。特别是水头（最常用的，有时是唯一的测量指标）有时几乎不包含关于水力传导系数的信息。 由于以上特征，含水层模型预测具有极高的不确定性。另外，参数估计不像其它领域那样可以明确地表达出来，在解释抽水试验结果时，就可能遇到这种情况。研究人员会将试验获得的水头值直接输入到程序中，经过适当的定性模型分析后，得到参数值。反演是含水层模拟过程中的一步重要工作。地下水的姐妹学科——地表水文学也具有许多类似的特征，其中的很多问题已经得到解决。 地下水反演的另一个重要特点是其相对独立性。目前，已经开发出了很多独立于其它领域的反演方法。最早的方法是将简单的等效水头（假定已知）代入流动方程中，这样可以得到关于导水系数的一阶偏微分方程（Stallman，1956）。Neuman（1973）将这一方法称为“直接法”，这种方法相对简单、易于理解，在Nelson（1960，1961）之后，得到了广泛应用。事实上，根据测量水头画出流网后，可以得到导水系数（Bennet和Meyer，1952）。然而，这一方法存在几个不足之处：首先，整个研究范围内（时间和空间）的水头值（以及补给、贮水系数和边界条件）必须是已知的，这样才能进行插值工作，但同时也会引入误差，使最终估计的导水系数受到了一些人为影响；其次，这一方法不稳定，水头微小的误差会造成导水系数很大的误差。为了解决第一个问题，最近大部分的反演方法采用了Neuman（1973）所称为的“间接法”，这一方法承认测量存在误差，并寻找可以将这些误差最小化的水力性质。也就是说，通过对目标函数的最小化来寻找参