地理学报 ›› 1999, Vol. 66 ›› Issue (6): 526-533.doi: 10.11821/xb199906006

• 论文 • 上一篇    下一篇

河流与含水层水力耦合模型及其应用

蒋业放, 张兴有   

  1. 中国科学院地理研究所,北京100101
  • 收稿日期:1999-04-03 修回日期:1999-06-23 出版日期:1999-11-15 发布日期:1999-11-15
  • 基金资助:
    中国科学院A 类重大科研项目(KZ2951-A1-203)

An Integrated ModelForm ulation ofStream - aquifer Interaction System and ItsApplication

JIANGYe-fang, ZHANGXing-you   

  1. Institute of Geography, Chinese Academy of Sciences, Beijing 100101
  • Received:1999-04-03 Revised:1999-06-23 Online:1999-11-15 Published:1999-11-15
  • Supported by:
    the key project of Chinese Academy of Sciences,No.KZ2951-A1-203

摘要: 本文在分析河水与地下水相互作用规律的基础上,提出了河流-含水层相互作用水力耦合模型的建模思路、数值解法及应用实例。河流模型采用忽略河槽调蓄作用的圣维南连续方程,含水层模型为潜水二维渗流方程, 二者通过动态水量交换机制实现耦合; 模型的水量交换机制可以分别模拟含水层顶托排泄、河流压力渗漏和淋滤渗漏等不同方式的水量转换过程。河水模型和地下水模型分别用欧拉法和三角网格差分法来求解。

关键词: 河流-含水层相互作用, 水力耦合模型, 三角网格差分法

Abstract: Based on theanalysisofhydraulicm echanism ofstream - aquiferinteraction sys- tem, an integrated m odelis developed to sim ulate the interaction process. To allow accu- rate sim ulation ofthesystem response, stream flow and ground waterflow aresim ultane- ously m odeled and a water exchange m echanism which can dealwith continuous flow and noncontinuous flow is designed to integrate the stream flow m odeland ground waterflow m odel. Mathem atically, the stream flow issim ulated by using Saint- Venant’s continuouse-quation w ith stream storage neglected and Euler’sApproach isem ployed to generate solu- tions to the equations. The stream system is discretized by im posing a grid ofretangular cells that are assigned num ericalvalues for selected hydrologic conditions and hydraulic properties. Ground w aterflow m odel, w hich isin the form ofphreaticaquiferequation, is solved with finite difference m ethod ofirregular triangulated netw orks. To integrate the m odels, the stream nodes overlap the aquifernodes where stream spass. The volum e of w ater exchange betw een stream s and aquifers is sim ulated by using Darcy’slaw, w hich isbased on thedifference betw een stream stage and watertable in the cell, and theverticalhydraulicconductivity ofthestream bottom, and itreachesthem axi- m um degree when noncontinuous flow occurs beneath the stream bottom . It can be ex- pressed as: (E= m in[γ·(Hr - (H) ·η,Em,qu] (H < (Hr γ·((Hr - (H) (H ≥(Hr (1-3) where E(m /d) is the exchangerate, γ= kzm, k(m /d) and m (m ) are theverticalhydraulic conductivity and thethicknessofthepoorly perm eablestrataon the stream bed respective- ly, (Hr(m ) is the stream stage, (H(m ) is the ground w atertable, Em (m /d) is the m axi- m um perm eability capacity ofthepoorly perm eablestrata and qu= Q/·1/W, Q,Land W are the stream flow, stream coordination and stream w idth respectively, and ηisa dim ension- less factorwhich isinversely proportionalto sand contentofthe stream w ater. Acasestudy showsthatthe integrated m odelcan accurately sim ulatestream flow and ground w ater m ovem ent and the calculated volum e of each stream section flow m atches properly w ith theobserved values. Thefactsprovethatthem odelresultsareaccurateand reliable.

Key words: stream -aquiferinteraction, integrated m odel, finitedifferencem ethod ofirreg- ulartriangulated netw orks

中图分类号: 

  • P343.1