降雨入渗补给量的数值模拟研究.pdfVIP

第十章敷学模拟技术 降雨入渗补给量的数值模拟研究+ 张志才陈喜 (河海大学水文水资源与水利工程科学国家重点实验室南京210098) 摘要降雨入渗补给是地下水评价的重要组成部分，该问题的研究对指导地下水告理开发、刺用具有重要的理 论和实际意卫。丰文主要探讨了术通量法在降雨入渗补培计算中的应用，并结合实例舟析研究了利用水通量法与 教值模型耦告椎求降雨^渗朴雏量曲方法。 关键词 降雨^渗补给水通量非饱和渗透系教数值模型 1引言 地下水资源的主要特点是具有可恢复性，它在蒸发、开采过程中不断消耗，同时也不断接受外界补 给?。地下水埋藏较浅的平原地区，降阿人渗补给为地下水(潜水)主要补给来源”’。降雨入渗补给量的 计算，是地下水资源评价的重要组成部分。准确计算入渗补给量对指导人们更合理地开发利用地下水资源． 避免囡过量开采地下水而发生的灾害．具有重要的理论和实际意义。 降雨人渗补给量估算主要有蒸渗仪观测法、水通量法等、示踪方法以及数值模拟方法。其中水通量法与 数值模拟方法研究与应用较为广泛．如Richards(1956)”o，De Boodt(1967)，Wellings(1984)”1利用水通 量法计算土壤水通最以度非饱和渗透系数，Cooper(1990)”1应用水通量法计算土壤下渗量。荆恩春等 (1989)[61利崩大型物理模拟实验和石家庄野外试验场试验，采用定位通量法、纠偏通量法计算降雨人渗补 给量。随着观测技术和计算方法的发展，以土壤水分运动数值模型为基础计算入渗补给量的方法得以迅速发 展。如Ragab(1997)等建立的简单土壤水模型(Sirr-ple Soil Model)以及Finch(1998)对此模型的改进， Gogolev(2002)提出的双层非饱和土壤水运动模型，使得利用数值法计算降雨人渗补给量更可靠。陈喜 (2002)提出的基于土壤水运动数值模拟和参数化方法的“四水”转化模型，成功地应用于淮北平原区降雨 人渗补给量的估算。 本文根据水文实验站实测资料，研究了利用水通量法与数值模型耦合汁算降雨人渗补给量的方法。 2水通量法 2．1 ZFP法 无水量通过的断面称为零通量面(ZFP)，该断面通常近似作为一个平面处理。若记ZIP为Z’，则存在 稳定ZIP的情况下，通过任一断面Z的水量可由公式(I)计算： q(z)=L口(z，l。)dx—L，p(；，b)d； (1) 式中：Q(z)为t。至t：时段内通过断面z的水量；o(z，“)，0(Z，b)分别为计算时段初l。和计算时段末屯断 面Z与Z’之问的平均土壤含水量。 乞2定位通量法 定位通量法的理论基础是达西定律和质量守恒原理。由于靠近潜水面土水势梯度变化剧烈，不宜利用达 西公式，因此可利用定位通量法计算通过潜水面附近任一断面z的水量，如式(2)所示： Q(z)=Q(Z’)+I，．日(z，I。)dz—L日(：，f2)d2 (2) 式中：∥表示某一给定断面；Q(Z)和Q(z’)分别为在t。～t：时段内断面Z和Z’处通过的水量。 · 基金项目：教育部新世纪优秀人才支持计划资助项目(NCET一04—0492)。 第一作者简介：张志才(1980一)，男．河北邯郸人．博士研究生．水文学及木资源专业。E—mail：zhangzhicai一0 @163 com 第=篇技术方法 若能求得通过断面z 7处的水量Q(z’)。由式(2)即可求得断面z处通过的水量Q(z)，断面∥处的水 量Q(Zr)可由达西公式计算。 2．3 ZFP法与定位通量法混合计算 ZFP法与定位通量法均可推求降阿人渗补给量。但由于实际中存在ZFP法失效期(ZFP消失)，因此采 用定位通量法计算连续降雨入渗补给量较好”’。如式(2)所示，应用定位通量法关键是确定p(z’)。因此， 确定非饱和渗透系数K(0)对计算降雨人渗补给量至关重要。根据z胛法，由式(1)、式(3)推求K(p)。 ， K(口)：一q(：)／譬 (3) ， 删 式中：q(：)=口(：)／(f：一f：)；Q(z)由式(1)计算。 3一维土壤水运移数值模型 水通量法方法简单，计算结果准确，但需要较多的实测数据，特别是土壤含水量分布资料，因此一般其 在实验室内进行模拟实验应用。通过模拟实验建立的数值模型通过参数率定和模型验证后可推广到野外较大 区域。因此利用水通量法与数值模型耦合推求降雨人渗补给景具有重要的意义。本文眦非饱和带分层土壤水 计算模型为例”1，探讨水通量法与数值模型耦台推求降雨人渗补给量的方法。 3．1土壤水分运动模拟 根据带有源汇项的一维Richard方程计算非饱和土壤水运动： 等=去(D罢)+iOK+，(枷) (4) 式中：K，D分别为非饱和渗透系数和非饱和扩散率，二者均为土壤音水量的函数。 设D(0)=一K(们(o啦／oo)，廿为土水势”‘，F(f，占