基于可变水文地质参数的地下水系统数学模型.pdfVIP

第20卷第3期 水科学进展 Vd．20．No．3 2009年5月 ADVANCESⅡ、『WrATERSCIENCE May，2009 基于可变水文地质参数的地下水系统数学模型 刘玉珍1一，王本德1，姜英震2 (1．大连理工大学土木水利学院，辽宁大连116024；2．辽宁省水利水电科学研究院．辽宁沈阳110003) 摘要：以辽宁省平原区为研究背景，建立了典型区地下水系统数学模型。提出了可变的开采含水层等效渗透系数和 ．平均给水度的水文地质参数模拟方法，对常规分区法进行了改进，并应用于建立的数学模型中．提高了模型精度。 用修正有限元法对数学模型进行数值离散，利用VisllalBasic软件编制计算程序，对地下水系统进行模拟，取得了 较好的模拟结果。建立的模型为进一步研究辽宁省平原区地下水状况和可持续利用提供了条件，研究方法对同类 地区的模型研究具有重要参考价值。 关键词：地下水；数值模拟；水文地质参数；有限元 中图分类号：P641 文献标识码：A 文章编号：1001．6791(2009)03．-0398-05 地下水是一种宝贵的自然资源，同时又是环境的基本要素…。实际的地下水系统具有地形起伏多变、含水 层厚度不均、结构复杂、水文地质参数的时空变异性明显等特征。数学模型对这些特征有相对合理的表征。是 对真实地下水系统的仿真和模拟…2。随着计算机技术的快速发展，数值计算方法在地下水资源分析评价中得到 逐步推广[3-8]。该领域学者对地下水系统数值模拟的工作程序、步骤方面达成一致，强调对水文地质条件合理 概化的重要性。对某系统的地下水进行资源评价，必须在对其水文地质条件进行必要概化的基础上，建立能够 准确描述该系统水流运动特征的数学模型…1。 近年来，虽然有很多数值法用于地下水资源评价的实例，但地下水系统是一个复杂的系统，在应用中。因 评价区域水文地质条件的差异，使得数值模拟方法在边界条件处理、水文地质参数分区，以及计算时段划分等 方面均有不同。水文地质参数的模拟是地下水系统数值模拟中面临的最大问题，在模型实际应用中，参数识别 可概括为分区法和插值法【9—2J，分区法即根据水文地质条件进行分区，每个区域内的水文地质参数被概化为 一个定值。该方法较为简单，在中国广为应用。当研究区面积不大或参数分区较多时，这种方法是有一定精度 的[13J。但也存在一定的局限性，第一，当实际研究区面积很大时，即使参数分区已经很多，精度也会大大降 低；第二，由于含水介质是在特定的历史时期特定的地质条件下沉积形成的，具有一定的成层性[13]，在地下 水开发利用过程中，随着地下水埋深的变化，对应的开采层在发生变化，随之发生变化的含水层水文地质参数 没有被反映出来。基于这种情况，本文结合辽宁省平原区典型实例，建立了地下水系统数学模型，在对其水文 地质参数模拟时，以分区法为前提，充分考虑以上两方面问题，提出了一种新的模拟方法，用动态可变的含水 层等效渗透系数和平均给水度来模拟水文地质参数，试图减少常规的水文地质参数分区法给地下水数值模拟带 来的模拟误差。 辽宁省是中国北方地区水资源严重短缺的省份，中部平原区人口密集、工农业发达，水资源供需矛盾突 出，选取的典型研究区域为秀水河河谷平原区，位于辽宁省中部，属于傍辽河干流发育的冲洪积河谷平原，总 面积995km2。该区域为河谷平原区第四系松散岩类孔隙水，基本以冲积、冲洪积物为主，沉积连续，厚度较 大，呈现出山前冲洪积扇、冲洪积平原构成的冲洪积物地下水水文地质特征，在平原区具有很强的代表性。 收稿日期：20嘀删 作者简介：刘玉珍(1966一)，女。辽宁北票人，教授级高级工程师。博士，主要从事水资源及农田水利方面的研究。 E-mail：hBky驷@163．corn 第3期． 刘玉珍，等：基于可变水文地质参数的地下水系统数学模型 399 1地下水系统数学模型的建立 依据二维地下水运动偏微分方程，在地下水资源系统概化的基础上，根据质量守恒定律及能量守恒定律， 建立系统的数学模型： (善，Y)∈n t0 (1)