复杂函数边界控制下的潜水非稳定流模型及解的应用-水利学报.pdf

水 利 学 报 2018年6月 SHUILI XUEBAO 第49卷 第6期 文章编号：0559-9350（2018）06-0725-07 复杂函数边界控制下的潜水非稳定流模型及解的应用 吴 丹，陶月赞，林 飞 （合肥工业大学 土木与水利工程学院，安徽 合肥 230009） 摘要：为解决复杂河渠水位边界影响下的潜水非稳定流模型难以求解的问题，建立不依赖边界函数的变换过程的 Fourier变换方法，利用卷积定义和卷积的微分性质，给出模型的理论解；对实际河渠水位过程采用Lagrange线性 插值，将插值函数代入理论解，可简便地获得问题的实际解。研究表明：（1）该方法求解过程比较简明，且解是 由形式较简单的常用函数组成；（2）依据潜水位变动速度的时间过程计算模型参数的配线法，方法简便；（3）边界 水位变化过程，对河渠与潜水之间的水量交换作用，有2倍于边界水位变幅的累积效应。 关键词：潜水非稳定流；Lagrange插值；Fourier变换；卷积微分性质；累积效应 中图分类号：P641.132 文献标识码：A doi ：10.13243/ki.slxb 1 研究背景 将完全切割潜水含水层的河渠水位作为第一类边界，边界一侧潜水非稳定渗流问题，是地下水 ［1-3］ 渗流力学中的最经典问题之一 。该问题也属带自由面渗流问题，所以一直是地下水渗流力学领域 ［4］ ［1］ 中的研究难点与热点之一 。上述问题的解，可为灌溉渠系设计提供重要的计算依据 ：通过调控渠 ［5-6］ 系水位将影响区地下水水位控制在合理区间，是现代许多灌区的基本要求 。 就边界条件而言，河渠水位变化过程不仅受河渠自身水力特征和区域水文气象条件控制，而且 ［7-8］ 受到明显的人工调控影响，具强烈的随机性 ；对这种复杂多变的时变函数所构成的边界条件，随 着对其概化与处理的方法不同，模型的求解方法以及所获得的解也有所差别。 这类问题的经典模型——J.G.Ferris模型（1950年），假定河渠水位变化瞬时完成之后长期 保持稳定［1，3 ，9-10］，即下文模型（Ⅰ）中的边界条件H （t ）为常数H0；对河渠水位复杂多变条件下的模 型，利用常用积分变换方法求解时，按对边界处理方式与求解方法不同，现行方法主要可分为三大 类： 方法一。H （t ）是如时间线性函数等相对简单的函数，可通过常用的积分变换方法对模型直接求 解［1，3，9-10］；该方法对一些特定情形（如经典模型所假设的情形）有较好的适用性，但不具普遍应用意 义。 方法二。H （t ）采用水位或水位变速“分段等值”方法离散化处理，依据经典模型的解，利用叠加 ［1］ 原理构建问题的解 。这类方法就求解过程而言，方法比较简便；但当H （t ）比较复杂时，划分的时 段数较多，计算偏繁琐。 方法三。H （t ）形式复杂，包括一些初始条件较复杂问题，在求解过程中或解的表达式中，需借 助不常用的特殊函数或是构建专用的特殊函数［4-6，11-14］，这为应用带来不便。综合上述方法优点，针对 收稿日期：2017-10-19；网络出版日期：2018-06-08 网络出版地址：http：///kcms/detail/11.1882.TV1510.001.html 基金项目：国家自然科学基金项目 作者简介：吴丹（1986-），女，安徽池州人，博士生，主要从事水资源研究。E-mail：wudanyt@163.com 通讯作者：陶月赞（1964-），男，博士，教授，博士生导师，主要从事水资源与水环境研究。E-mail：taoyuezan@163.com