基于有限层法的承压水开采地面沉降简明分析.pdf

书书书第 卷第 期烟台大学学报（自然科学与工程版）年 月 （）文章编号：（）：收稿日期：基金项目：国家自然科学基金资助项目（）；江苏省青蓝工程；南通市科技局计划项目（）。通信作者：徐进（），副教授，博士，主要研究方向为岩土工程数值计算。基于有限层法的承压水开采地面沉降简明分析徐进，张豪，王斯海（烟台大学土木工程学院，山东 烟台 ；江苏工程职业技术学院建筑工程学院，江苏 南通 ）摘要：针对承压含水层地下水开采问题，考虑地下水的非稳定流和不同开采条件的影响，基于地下水流的有限层方法，给出了开采条件下的地下水位降深。结合弹性半空间明德林解和二维卷积公式，建立了承压水开采引起上覆土层地面沉降的计算公式。在此基础上，利用 实现了地面沉降的简明分析，通过与已有解的对比，验证了所建立计算方法的正确性。通过算例分析进一步详细分析了覆盖层厚度、渗透系数、贮水率、弹性模量等因素对地面沉降发展的影响，探讨了本文方法应用于水平井、群井等复杂抽水问题的适用性。关键词：地面沉降；承压水；地下水开采；简明计算；二维卷积中图分类号：文献标志码：对含水层中地下水进行开采会破坏原有的水土平衡状态，在地下水位发生降深的同时，引发土层变形，导致地面沉降。同时，人类工程活动中对潜水和承压含水层的降水减压也会导致工程性沉降，对周围建筑物和市政设施造成不利影响。如何正确评估上述环境地质问题日益得到学术界的重视 。从机理来说，地面沉降是一个复杂的流固耦合过程 。基于水 土耦合理论，和 较早提出了地面沉降有限元计算模型，并运用于威尼斯的地面沉降问题 。等考虑了潜水面变化，研究了抽水过程中潜水面变动对地面沉降的影响 。金小荣等利用二维有限元模型分析了工程降水作用下基坑周围土体沉降性状的影响因素 。骆祖江等系统建立了地下水开采引起地面沉降三维全耦合有限元模型，引入了渗透系数动态模型和邓肯 张非线性模型 。郑刚等利用有限元软件建立三维流固耦合模型研究了承压层长期局部降压的土体分层沉降机理和规律 。等基于半解析数值原理给出了抽水作用下三维流固耦合问题的求解格式 。等建立了基坑降水地面沉降的耦合模型，分析了多层地下水系统基坑降水引起的支护结构和周围土体响应 。但是，基于水 土耦合的地面沉降计算模型，由于计算工作量大、参数难以获取，在实际应用中仍受到一定限制。鉴于此，便于实际应用的简化计算方法得到较多关注。骆冠勇等基于一维太沙基固结方程，推导了承压水降压沉降的固结度计算公式 。龚晓南等结合地下水动力学公式与弹性半空间理论，给出了承压水稳定流情形下的最大沉降公式 。为了进一步反映抽水时间的影响，王春波等将其推广到了非稳定流情形 。地下水流计算是地面沉降分析的基础环节，如何实现三维地下水流的高效计算一直备受关注。有限层法是一种典型的混合算法，兼具解析方法和数值方法的优点，其精确高效的特点在地下水流问题中已得到充分体现 。为此，基于地下三维流高烟台大学学报（自然科学与工程版）第 卷效有限层方法，针对承压水开采引起的地面沉降问题，结合明德林解答和二维卷积公式，得到地面沉降的简明计算公式。借助 实现求解，验证计算公式的正确性，并对承压水降压下地面沉降的影响因素进行分析。承压水三维流有限层法图 所示为典型的下卧承压含水层，在其中利用抽水井开采地下水，在抽水井附近将出现地下水位降深，形成地下水的渗流运动。根据达西定律和质量守恒定律，抽水作用下引起的地下水三维非稳定流控制方程如下：（）（，），（）式中：（，）为降深；为含水层渗透系数；（，）为汇源项；为贮水率。式（）中 （，）表示单位时间从单位体积含水层中抽取的地下水量，表达式取决于抽水井的具体形式，对于图 所示的竖向非完整井，假如井中心位置为（，），则取如下形式：（，）（）（）（），（）式中：为 函数，（），其余，为井滤管长度，为抽水量。?图 承压含水层抽水井开采示意图 有限层法是一种介于解析法和数值法之间的混合算法。区别于计算域的全局离散，有限层法沿 方向将含水层离散成 个层元，水位降深在 、平面采用正交完备的解析函数系 （，）表示，结合 方向的有限元形函数（），得到级数形式的降深试探函数 ：（，）（）（，）（）。（）根据标准的有限元格式，得到关于未知系数的线性方程组：，（）式中：为未知系数向量；渗透性矩阵 、贮水矩阵 和流量矢量 可以由具有解析形式的公式计算得到，具体形式参见文献 。在传统的数值方法中，数值积分是进行单元分析的必要环节，占到相当大的计算工作量。相比来说，有限层法不涉及数值积分，在地下水三维非稳定流分析时，计算工作量更低。同时，由于试探函数中解析函数的正交性，求解格式具有天然的解耦性，通过并行计算可以进一步提高计算效率 。承压水开采引起的地面沉降 水位降深引起的附加应力在抽水作用下，承压含水层中形成的降落漏斗将打破原有的平衡状态，导致上覆土层浮托力的减少，上覆土层底部（承压含水层顶板）形成竖向的附加应力，这种由水位下降导致的竖向附加分布力为（，）（，），（）式中，是水的单位重度。地面沉降计算公式承压含水层顶板为隔水层，不用考虑含水层水位下降引起的潜水层越流和上覆土层的固结，地面沉降主要由含水层顶板处的附加作用力（）引起。根据弹性力学的明德林解答可知，无限半空间体内，单位竖向集中力作用下，地表的竖向位移响应为 （，）（）（）（），（）式中：和 分别为弹性模量和泊松比；是上覆土层的厚度。基于二维卷积公式，结合单位集中力位移响应（）和级数降深函数（），可以得到分布力（）作用下的地表沉降计算公式（，）（）（，）（）（，）。（）利用地下水流的有限层求解格式（），将得到的降深级数计算解代入式（）进行积分，可以分析抽水作用下的地表沉降过程。如果需要对土层内部的位移、应力等情况进行分析，可以用相应物理量的第 期徐进，等：基于有限层法的承压水开采地面沉降简明分析明德林解在式（）中对地表沉降 进行替换。式（）难以进一步展开成解析形式，本文利用 进行求解。对比验证与影响因素分析 对比验证为了验证该方法的正确性，选取典型算例进行了分析，其中下卧含水层中抽水井为完整井，地下水流为稳定流。图 给出了不同弹性模量取值（，）下，地表中心沉降 随覆盖层厚度变化的计算结果，并将其与已有解进行了对比。可以看出，三种情况下本文给出的计算结果都与已有解 吻合良好，验证了本文方法的正确性和合理性。同时，图 计算结果也反映了上覆土层厚度 对地面沉降的影响，在特定弹性模量取值下，随着 的增加，含水层埋深越大，地下水开采的影响越小，所引起的地表沉降越小。当土层厚度 ，泊松比 ，时，在承压水降深 的条件下，地面中心沉降接近 ，这与文献 中江南工作井的计算和实测结果相吻合。?图 地表中心沉降随覆盖层厚度的变化 地下水开采是地面沉降发生的诱因，本文研究了开采量对地面沉降的影响，见图 。可以看出，开采量对地面沉降具有直接影响，随着开采量的增加（从 增加至 ），地面沉降明显加剧，中心沉降从 增大到接近 ，两者呈正相关，保持良好的线性关系。因此，结合当地水文地质条件，特别是开采含水层的补给条件，制定合理地下水开采量是地面沉降灾害预防的关键。文献 的地面沉降计算公式可以分析中心沉降，不能给出地面沉降的空间分布。利用本文方法给出了抽水井附近的地面沉降分布，如图 所示。可以看出，竖向单井抽水作用引起的地面沉降呈漏斗状分布，与水位降深分布保持一致，最大沉降出现在抽水井中心，随距井径向距离 的增加逐渐减小。同时，由于本文方法中考虑了地下水的非稳定流，可以给出地面沉降的发展过程，图 给出了不同时刻（，）的沉降曲线，可以看出，随着抽水时间的延续，井附近的地面沉降不断加剧。?图 抽水量对地面沉降的影响?图 地面沉降随抽水时间变化 影响因素分析地面沉降涉及地下水运动和土体变形的相互作用，影响因素很多。本文计算公式建立过程中，虽然采用了一定简化（忽略上覆土层的固结变形），却考虑了一些重要水文地质参数（如渗透系数 和贮水率）与土力学参数（如弹性模量 和泊松比 ）的影响。图 给出了渗透系数 对沉降的影响，可以看出，数值越大，含水层内水力联系越强，井附近的烟台大学学报（自然科学与工程版）第 卷水位降深越均匀，中心沉降和差异沉降越小。图 表明含水层贮水率 对沉降也有显著影响，随着 的增加，含水层中的地下水储量越大，抽水作用下含水层内的水位降深越小，所引起的沉降值也越小。?图 渗透系数对地面沉降的影响?图 贮水率对地面沉降的影响 图 和图 给出了上覆土层力学参数的影响，可以看出，其他条件保持不变的情况下，覆盖土层刚度越大，土体变形越小，沉降量越小，当弹性模量 从 增加到 时，最大沉降量下降了 ；相对而言，泊松比 对沉降的影响较小，当 在从 增加到 时，最大沉降值减少了 。多种抽水形式下的地面沉降近年来，水平井凭借比常规竖向井更高的抽水效率，在 环 境 和 水 文 方 面 得 到 越 来 越 多 的 应用 。与普通的井网结构不同，水平井的井过滤器水平放置（图 ）。在地下水运动特征上，一般来说，水平井引起的水流模式比传统竖向井更加复杂，具有明显的三维流，并缺乏轴对称性 。?图 弹性模量对地面沉降的影响?图 泊松比对地面沉降的影响?图 含水层中水平抽水井示意图 在有限层求解格式中，水平井过滤器可以完全嵌入剖分节面，从而只需对求解格式（）中 进行适当修改，无需离散就可以实现求解 。基于本文方法，水平井附近的地表沉降曲线如图 ，为了便于对比，图中也给出了相同抽水量下竖向井的地表沉降分布。可以看出，相比于竖向井，水平井开采引起的地表沉降不存在轴对称性，沿过滤器方向的沉降和垂直方向的沉降并不一致，沿井轴方向沉降第 期徐进，等：基于有限层法的承压水开采地面沉降简明分析值更大。此外，总体上，水平井引起的沉降分布比竖向井更为平缓，这与水平井的水力学特征相吻合 ，。?图 水平井和竖向井引起的地面沉降分布曲线对比 为了优化开采方案，在含水层中通常会设置群井系统。在传统有限元数值方法中，为了获得精确的近井水流模拟结果，每一个抽水井附近，都需要更精细的网格划分，这会带来额外的工作量，甚至造成网格不规则和离散误差。有限层法中，由于采用具有解析表达式的流量矢量 代替网格划分，不存在上述问题，因此在处理多井问题时更加方便。为了说明本文方法处理多井问题的能力，设计了由三口竖向井和一口水平井组成的群井系统，具体计算参数见表 。图 给出了群井开采下地面沉降的等值线图，可以看出，在抽水井附近形成了四个局部沉降锥。三口竖向井附近的等高线近似于同心圆，水平井附近的等高线为椭圆，正确反映了竖向井和水平井的水力特性，说明了本文方法应用于群井问题的适用性。表 群井问题计算参数 参数取值含水层厚度 含水层渗透系数 含水层贮水率 上覆土层厚度 弹性模量 泊松比 抽水量 计算区域尺寸?图 多井抽水地面沉降等值线 结论针对开采承压水诱发的地面沉降问题，基于地下水三维非稳定流的有限层方法，利用二维卷积公式和弹性半空间明德林解，提出了一种简明的地面沉降计算方法。通过与已有解的对比，验证了本文方法的正确性。本文简化计算公式可以反映地面沉降的空间分布及其时间演化过程。通过算例分析，研究了多种因素对抽水地面沉降的影响规律，结果表明：（）上覆土层厚度越大，地面沉降量越大；弹性模量对地面沉降影响较大，覆盖土层刚度越大，地面沉降越小。相对而言，泊松比的影响较小；（）含水层渗透系数和贮水率对地面沉降也有一定影响，渗透系数越大，含水层内水力联系越强，差异沉降越小；含水层的储水量越大，地面沉降越小；（）对水平井、群井问题进行了分析，所得结果符合相应的水力学特征，说明该方法在处理复杂抽水问题时也具有一定的适用性。参考文献：龚士良国务院审批同意 全国地面沉降防治规划（年）实施四大工程，有效控制地面沉降恶化趋势 上海国土资源，（）：吴林高，等工程降水设计施工与基坑渗流理论 北京：人民交通出版社，张云，薛禹群抽水地面沉降数学模型的研究现状与展望 中国地质灾害与防治学报，（）：，（）：烟台大学学报（自然科学与工程版）第 卷 ，（）：金小荣，俞建霖，祝哨晨，等基坑降水引起周围土体沉降性状分析 岩土力学，（）：骆祖江，王琰，田小伟，等沧州市地下水开采与地面沉降地裂缝模拟预测 水利学报，（）：郑刚，曾超峰，薛秀丽承压含水层局部降压引起土体沉降机理及参数分析 岩土工程学报，（）：，（）：，：骆冠勇，潘泓，曹洪，等承压水减压引起的沉降分析 岩土力学，（）：龚晓南，张杰承压水降压引起的上覆土层沉降分析 岩土工程学报，（）：王春波，丁文其，刘文军，等非稳定承压水降水引起土层沉降分布规律分析 同济大学学报（自然科学版），（）：，（）：，（）：王少伟，徐进，杨伟涛地下水流并行有限层方法及同伦反演研究 计算力学学报，（）：，（）：，：，：，（）：，（，；，）：，：；（责任编辑李春梅）