Page 107 - 2023年第54卷第10期

图 ４ 渗透系数场 Ｋ ３

图 ５ 基于 ＥＳ － ＭＤＡ方法两次同化反演参数的后验分布
图 ６分别为（ａ）参数 Ｋ 标准场、（ｂ）反演前随机生成的初始参数场、（ｃ）基于 ＥＳ － ＭＤＡ第一次、
３
－ ２
（ ｄ）第二次同化反演的参数场。两次同化反演参数场与参数标准场之间的 ＲＭＳＥ分别为 １．８ × １０ 和 ９．１ ×
－ ３
１０ ，ＳＳＩＭ分别为 ０．９３和 ０．９６。两者都较为准确，且第二次同化反演参数场精度更高。说明增加数据
同化次数可以增加 ＥＳ － ＭＤＡ方法的参数反演精度，第二次同化反演的后验分布更为精确。
３．４ 正演模型验证 在非均质参数场和开采量的校验后，还要进行正演模拟，用来验证预报水头场的
能力。通过将基于两阶段的 ＭＣＭＣ和 ＥＳ － ＭＤＡ所反演出的非均质参数场和开采量输入到地下水数值
模型中，与监测井真实水头值和真实水头场进行对比。如图 ７所示，其中红色实线为真实水头，绿色
虚线和蓝色点划线分别为基于两阶段的 ＭＣＭＣ和基于 ＥＳ － ＭＤＡ参数反演后正演模拟得到的监测井
水头。
根据图 ７可以看出，基于两种方法反演的非均质参数场和开采量进行正演模拟的监测井水头过程
线与真实水头过程线十分接近，两种方法都充分捕捉到了监测井位置的水头变化情况，说明监测井水
头模拟效果较好。

２
— １ ４ ３ —