基于Monte Carlo方法的地表水地下水耦合模拟模型不确定分析

为分析参数的不确定性对地表水地下水耦合模拟模型输出结果的影响,本文运用Monte Carlo方法对地表水地下水耦合模拟模型进行不确定性分析,并根据不确定性分析结果进行生态恶化风险评估。以石头口门水库上游饮马河汇水流域为例,建立地表水地下水耦合模拟模型,并运用HydroGeosphere(HGS)软件求解。利用局部灵敏度分析方法甄选出耦合模拟模型中灵敏度较高的参数,作为随机变量。最后,运用Monte Carlo方法对耦合模拟模型进行不确定性分析。为了减少不确定性分析过程中的计算负荷,应用Kriging方法,建立耦合模拟模型的替代模型。结果表明:地表水地下水耦合模拟模型中灵敏度较高的参数为渗透系数、孔隙度和曼宁粗糙系数,其中渗透系数的变化不仅对耦合模型中地下水水位产生影响,也对地表水流量产生影响;Kriging替代模型可以在保证一定精度的条件下大幅度减少计算负荷;风险评估结果表明,在当前的水资源开发利用条件状态下,研究区地下水生态环境恶化风险的概率为6%,地表水生态环境恶化风险概率为15%。