径向基函数配点法在分析承压和非承压地下水混合流动计算中的应用.pdfVIP

　 2015年3 月 地下水 Mar. ,2015 第37卷　 第2期　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Groundwater　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Vol．37　 NO．2 径向基函数配点法在分析承压和非 承压地下水混合流动计算中的应用 李永莲,周德亮,牛忠星 (辽宁师范大学数学学院,辽宁 大连 116029) [摘　 要]　 采用径向基函数配点法建立承压和非承压地下水混合流动问题的数值模拟模型。 模拟过程中针对 非承压含水层出现的非线性的配点方程组,利用迭代法将其转化为线性的配点方程组,然后再计算求解。 结果表 明,迭代法对于非线性的配点方程组也是有效的;同时也说明了径向基函数配点法求解该问题不仅实施简单,计算 效率较高,而且计算结果令人满意。 [关键词]　 径向基函数;配点法;承压;非承压;地下水;非线性 [中图分类号]　 P641.2　 　 [文献标识码]　 B　 　 [文章编号]　 1004 -1184(2015)02 -0047 -02 　 　 目前,水文地质领域内使用的数值模拟方法主要有有限 沉积含水层补给;河谷地区上游地段,为一接受灌溉、降雨补 [1] 单元法、有限差分法、边界元法等 。 这些方法对网格具有 给率为ε 的地区;中、下游地段的冲积、沉积含水层顶部,覆 很强的依赖性,且网格剖分容易造成局部区域不守恒,计算 盖有厚8 cm左右的不透水覆盖层,因存在不透水覆盖层,灌 [2] 量大,成本高等缺陷。 而无网格方法 (无单元方法)则不需 溉、降雨的余水,表现为地表径流不参加含水层内部的地下 要预先定义网格或单元,避免了网格生成的复杂过程,具有 水流动。 图3所示的剖面A-A,假定平行于西断层方向,剖 节点配置灵活,计算量小,计算程序易实现,计算效率较高等 面形状不变,以此来反映全区域冲积、沉积含水层的变化形 [3] 优点 ,是求解微分方程定解问题的一类新的有效的数值计 态。 [4] 算方法。 径向基函数(Radial Basic Function,RBF) 配点法 作为无网格方法的典型代表,具有实施简单、不需要积分的 [5] [6] 特点,在一维土壤水分运动 、多层介质热传导反问题 、突 [7] 发性水污染事故的反演模拟 等问题的微分方程求解中已 获得了成功的应用。 本文采用径向基函数配点法解决平面二维承压和非承压 地下水混合流动问题,建立具体的数值模拟模型,用 MAT- LAB程序实现该算法,并对结果进行实际意义的分析。 模拟 过程中基于承压和非承压地下水方程的差异,对每个节点进 行承压和非承压的判断。 而对于非承压含水层出现的非线 性的配点方程组,为简化计算,采用迭代法将其转化为线性 的配点方程组来求解。 1　 径向基函数配点法求解承压和非承 图1　 河谷地区平面图 压地下水混合流动问题 1.1　 研究区概况及概念模型 [8] 考虑承压和非承压地下水混合流动问题 :如图 1 所 示,两河相会,两断层相交,其间围成一地堑形河谷地区。 该 地区存在两条断层,南断层为一下降断层,倾向北东,其南侧 为大体积花岗岩墙,视为不透水边界;西断层也为一下降断 层,倾向南东,其西侧又为一大体积花岗岩墙。 该地区东侧 和北