基于MCMC和ES-MDA方法的地下水数值模型非均 质参数场及开采量的反演研究

马尔科夫链蒙特卡罗方法(MCMC)和多重数据同化集合平滑器方法(ES-MDA)近年来在地下水参数反演得到广泛应用，但对三维多层非均质含水层参数反演精度和计算效率还缺乏对比分析。本文构建了含有基于Karhunen-Loève展开的非均质参数场的潜水和多层承压水含水层案例，并建立了地下水数值模型和基于Kriging方法的替代模型，模拟含水层分层水头变化，探讨了基于替代模型的MCMC、替代模型和数值模型相结合的两阶段MCMC以及ES-MDA方法反演的含水层渗透系数以及开采量。结果表明，针对本文算例，在非均质参数和开采量的反演中，相比而言，两阶段MCMC反演参数精度更高，ES-MDA方法计算效率更高。本研究为地下水数值模型参数反演方法选择提供参考依据。