地下水动态随机模型的研究进展.docxVIP

地下水动态随机模型的研究进展

1随机理论方法

地下水动态值的模拟是对地下水资源和地下水环境变化的定量分析和处理方法。以前的数值模拟一般采用确定性方法,根据质量守恒定律,建立连续性偏微分方程定量描述。但是多孔介质空隙的几何条件复杂多变,各种水文地质参数都呈空间随机分布,致使地下水在多孔介质中的运动具有很大的不确定性。地下水流的确定性模拟方法由于未考虑地下水流系统的随机性问题,用它来描述地下水水流运动时,其模拟值存在着不准确性。为解决上述问题,人们提出了随机理论方法,它将含水层的水力参数视为随机场,可根据有限的数据模拟地下水水流的空间变异规律。2000年,美国国家研究委员会召开会议,对地下水研究建议了四个方向,其中第二个就是要基于野外实验数据,开发随机模型的研究方法学,来研究具有尺度效应和不确定边界的各向异性含水层。本文重点介绍了随机理论方法在地下水动态模拟研究中的原理、各种随机模型的应用领域及局限性,及近几年来发展的各种随机估参方法。用随机理论方法对地下水动态进行模拟,可以在给出参数估计置信区间和估计精度的基础上,寻求各种随机因素的空间变异性对水流运动的影响,对地下水基础理论研究及地下水资源管理有一定的意义。

2随机动态水模型的基本思想和应用条件

2.1与期望值相关度的把握

在地下水动态模拟预测中,人们真正关心的是区域范围内的地下水动态(考虑到随机因素的影响,应指地下水动态的变化趋势或期望值),以及由随机因素而产生的波动幅度(即与期望值的偏离程度)。随机动态模拟实质上是把确定性模型中的参数看成具有空间分布的随机场,把随机场离散后带入到控制方程中,通过研究随机参数的统计特征(一阶矩和二阶矩)和变异性描述地下水水流变量的宏观规律性,可以给出随机参数及随机变量的估值精度和估值区间,以及由于随机参数和随机变量的波动对水位动态变化的影响程度。

2.2模型假设简化

随机方法在研究地下水水流运动时,必须满足以下条件:渗透系数各向异性、孔隙度呈空间分布;渗透系数的对数在其平均值处随机波动,但整个范围内波动相对较小,无周期性或某种趋势性;水头在其均值附近浮动,水流运动不在某一确定方向上;为便于分析,模型必须通过合理的假设进行简化,假设扰动较小,并进行线性近似;与渗透速率波动的大小相比,当地(实验室)弥散度波动较小。

3随机数值模拟

地下水动态预测最早采用的是比较简单的水均衡法及水文地质比拟法。1905年,E.梅勒第一次用解析法论证了泉水流量的预测方法,20世纪50年代后期,Г.Н.卡门斯基在考虑降水入渗的条件下,通过有限差分法预测了地下水动态变化规律。上个世纪80年代后期,随机方法在国外迅速发展。国内20世纪70年代中后期才开始采用确定性方法预测地下水动态,随机模拟近十几年刚刚发展起来,研究很有限,孙讷正(1989)与陈葆仁(1988)等人分别在其专著里对随机数值模拟进行过简单的介绍,杨金忠(2000)对多孔介质中水分及溶质运移的随机理论进行了总结性介绍。

随机模型按随机因素的性质可分为两类:1)仅含随机边界条件或随机初始条件的模型,即模型中的控制方程为确定性方程,但定解条件是随机的;2)带有附加随机函数和随机系数的模型,即控制方程中的某些参数或函数是随机的,这种分法较易理解。按求解方法可分为:①Monte-Carlo随机模型;②Taylor展开随机模拟;③摄动展开随机模拟;待定系数摄动展开随机模拟;④待定系数摄动展开随机有限元和待定系数摄动展开有限体积法;还有其它一些随机模型,如泊松模型、球状模型、随机分区模型等。

3.1数值计算模型

1965年,Warren和Price最早把Monte-Carlo模型应用于地下水研究,此模型假定随机变量的概率分布函数和协方差函数已知,用伪随机数的生成技术产出多组随机变量,然后把变量带入模型求解水流的统计值。此模拟技术,不受随机变量波动范围的限制,可计算不同节点同一时刻的浓度(水头)协方差,同一节点不同时刻的浓度(水头)协方差等。Freeze(1975)、Sagar和A.Runchal(1982)、El-Kadi和Brutaert(1985)发展了此模型。Yoneda等(1986)根据Monte-Carlo随机模拟方法把任意分布的随机变量转化为正态随机分布,计算了地下水水位。崔秀等(1997)采用积分插值法与Monte-Carlo随机法对二维承压稳定地下水流进行计算和理论研究,其中随机分布数采用混合同余法产生,最后对这两种方法的计算格式及边值条件处理进行了推导、分析。姚磊华(1997)讨论了求解地下水模型的诺伊曼展开、蒙特卡罗随机有限元法,对二维承压稳定流中随机变量的计算时段进行了重点分析。紧接着(1998)又建立二维潜水水流和水质数学模型,根据渭河水位、水质、降水量等随机变量的特征