地下工程渗流场有限元模拟研究

在水电、石油、交通等部门的地下工程中,经常会遇到地下水的渗流问题,如果处理不当,可能诱发应力恶化而影响结构稳定。本文以非饱和介质的渗流模型为基础,实现了地下洞室开挖过程中的非稳定渗流场解析。
按分期施工方案,在保持渗流场变化连续性的前提下,对各施工阶段进行了准稳定渗流分析,并在此基础上建立了地下工程非稳定渗流数值模型。其中,对于洞室系统,基于水量平衡原理,采用了子结构模型;对于防渗帷幕,采用了隔水层模型;对于排水井列,采用了杆单元模型。
根据上述模型,采用有限差分法作为数值模拟方法,编写了地下工程非稳定渗流场分析的有限元程序,实现了复杂渗流区域的渗流场的动态模拟。工程实例计算分析表明,用该模型计算的结果是合理的。      

一种真空预压区地下水位计算方法

根据抽真空时加固区内真空渗流场的形成机制和变化特点,把真空渗流场的形成过程分为非稳定渗流和稳定渗流两个阶段。利用渗流量计算公式,计算加固区内渗流量随时间变化过程曲线。结合两个阶段加固区地下渗流场的变化特性,分别利用Theis和Dupuit计算模型计算不同时刻的地下水位,并用实测进行验证。     

排水砂槽渗流试验与Signorini型变分不等式方法验证

排水是岩土体及工程构筑物渗流控制的主要措施之一,其实质是通过在渗流域内形成潜在溢出边界或低水位边界而实现渗流控制的。在排水渗控条件下,渗流场往往具有强烈的边界非线性特征,Signorini型变分不等式方法从理论上为稳定和非稳定排水渗流问题提供了有效的分析方法,但其实际效果还需要得到试验的验证。通过开展含5个排水廊道的排水砂槽模型试验,研究了复杂排水条件下砂土渗流的基本规律,并通过试验数据与数值计算成果的对比分析,论证了Signorini型变分不等式方法的有效性和正确性。试验结果表明,在稳定渗流条件下,排水砂槽上游侧3个排水廊道对渗流控制起主导作用,而下游侧2个廊道则失去排水功能,数值计算与试验成果吻合较好;在非稳定渗流条件下,受测压管精度、砂样均匀性和毛细效应的影响,数值计算与试验成果存在一定偏差,但也较好地揭示了复杂排水条件下砂槽中的非稳定渗流特征。排水砂槽试验结果验证了Signorini型变分不等式方法的有效性和正确性,为复杂排水条件下岩土体及工程构筑物的渗控结构优化设计提供了有效的分析手段。     

三维稳定渗流的 p 型自适应有限元分析

建立了三维稳定渗流的p型自适应有限元分析的误差估计器。充分利用p型有限元前处理少、计算精度高和收敛速度快的优点，分析了p型有限元法求解渗流问题的技术特点，给出了求解三维稳定渗流的自适应升阶策略。通过具体算例，研究了单元阶次及网格尺度对计算结果的影响，验证了三维p型有限元法在求解渗流问题时的可行性。     

泄洪雾雨作用下边坡饱和非饱和非稳定渗流研究

泄洪雾雨作用下的边坡岩体饱和非饱和渗流问题是目前我国西南地区在建和已建高坝所面临的一个共同问题。从理论上分析了雾雨作用下饱和非饱和介质中水分运动特征;根据饱和渗流和非饱和渗流的数学模型,建立起统一的岩体饱和非饱和非稳定渗流数学模型;编制了有地表入渗作用的饱和非饱和非稳定渗流三维有限元计算程序;对有限元求解中的容水度问题、非饱和水力参数问题以及初始水头场问题作了相应的优化处理;以实际工程为例进行了饱和非饱和渗流分析,计算结果表明,比天然降雨雨强大很多的雾雨入渗会形成对边坡稳定不利的暂态饱和区,并引起地下水位的抬高;暂态饱和区的大小和地下水位的增幅取决于边坡表面的护坡效应和排水、防渗帷幕等措施,因此应该加强边坡表面的保护以及合理的设置排水、帷幕等措施,以减少入渗和地下水位的抬高。     

基于非饱和土固结理论的有限元强度折减法

库水位下降使库岸坡体内产生饱和－非饱和非稳定渗流,非稳定渗流作用是边坡失稳的重要原因之一。为了分析非稳定渗流对边坡稳定性的影响,首先推导了固、液两相孔隙介质的固结方程,考虑了土体渗流与变形的耦合作用。在此基础上结合有限元强度折减法求解边坡稳定安全系数,将渗流、变形及稳定分析采用一套统一的有限元方法。并通过算例分析了库水位骤降情况下,坡体的渗透系数、水位降落比对稳定安全系数的影响,计算结果表明,所提的理论和方法是有效可行的,为饱和-非饱和非稳定渗流作用下边坡稳定问题的分析提供了实用工具。     

台兰河地下水库辐射井抽水过程的非稳定渗流场的有限元分析

针对台兰河地下水库大口径辐射井抽水效果的数值模拟问题,在分析了大口径辐射井渗流行为基础上,提出了应用辐射井子结构法和子结构开关器等辐射井精细模拟技术,并联用改进的截止负压法、迭代增量法和求解大型稀疏矩阵的预处理共轭梯度算法进行有辐射井影响的地下水非稳定渗流场的有限元分析。应用自主研发的三维可视化渗流有限元计算软件GWSS对台兰河地下水库辐射井抽水试验过程中的非稳定渗流场进行了数值模拟研究。研究结果表明：所有观测井计算水位的平均绝对误差为0.22 m,单井水位的平均误差最大值为 0.40 m,最小值为0.02 m,各观测井地下水位的模拟值与实测值的变化趋势吻合较好。经验证算法的合理性和程序的可靠性较好。提出的辐射井子结构法可以精细模拟辐射井的渗流行为和局部精细渗流场,可用于分析辐射井抽水过程中的地下水非稳定渗流变化及取水效率评价。     

祁连山疏勒河上游基流变化及其影响因素分析

基流是河川径流的重要组成部分,也是干旱半干旱地区最重要的枯季水文特征之一. 基流在维持河川径流和维护河流湖泊生态等方面具有重要意义. 基于祁连山疏勒河上游出山口控制水文站昌马堡站1954-2009年的日径流资料,采用递归数字滤波法和平滑最小值法对日径流进行基流分割,对比分析基流和基流指数BFI的变化特征,并从气温和降水两方面探讨基流变化的主要影响因素. 结果表明：两种方法的计算结果有差异,但平滑最小值法更适用于研究区. 疏勒河上游多年平均基流量为6.07×10⁸m³,基流变化总体上呈上升趋势,进入21世纪以后基流量增加最为显著. 研究区多年平均基流指数为0.66,基流指数的年际变化相对比较稳定,但年内波动幅度较大. 在长时间尺度上,降水是基流变化的主要影响因素,降水量的增加在一定程度上影响了基流和基流指数的变化. 在全球变暖的背景下,气温升高导致冰雪融水的增加和冻土的退化,间接的影响了基流的变化.     

岩体三维非稳定渗流模型及数值模拟

地下水是导致地质灾害发生的重要因素之一。岩体本身的复杂性及地下水在岩体内随时间而变化造成地下水在岩体内运动机理的复杂性。文章采用数值模拟方法探索地下水在岩体内的运动规律,了解其致灾机理。以岩体结构力学为基础,运用现代分形理论,建立主干裂隙分形网络;根据裂隙发育规模与工程尺度关系将岩体看作拟连续介质与块裂介质混合介质,并根据两类介质接触处水头相等及节点流量相等建立合理的三维非稳定渗流模型;给出渗流模型的有限元解法,开发出相应的有限元软件;给出算例,计算结果体现了主干裂隙在渗流中的强导水作用及网络状裂隙的贮水功能与渗流滞后效应;强调根据裂隙发育规模与工程尺度关系确定合理渗流模型的重要性。     

非饱和堤岸的渗流特征及其稳定性研究

通过试验确定了非饱和黏土层的土-水特征曲线和强度参数,通过理论分析建立了非饱和堤岸渗流-应力耦合模型。在此基础上,利用耦合计算程序分析了非饱和土堤岸在河水位变动时的非稳定渗流场特征,结合强度折减有限元法分析了河水位反复升降后非饱和堤岸稳定性的变化。结果表明：水位快速上升时,堤脚渗透流速下降后又逐步上升并趋向稳定,河水位快速上升使堤岸边坡的稳定性降低,随着渗流场中孔隙水压力的调整,堤岸边坡的安全系数又有所回升;河水位骤降时,黏性土层饱和度变化相对滞后,水位骤降加大了堤岸的渗透流速,边坡的稳定性迅速降低,水位下降约120 h后堤岸边坡进入较危险时段;非饱和土的基质吸力提高了堤岸边坡整体稳定性,河水位反复升降降低了堤岸边坡整体稳定性,且河水位越低,水位反复升降对堤岸稳定性产生的影响越大     

复杂水力路径下非饱和流动的数值分析

多孔介质中的流动问题,与孔隙介质的特性,含水量状态以及含水量的变化历史密切相关。基于毛细循环滞回理论模型,考虑含水量变化历史对土水特征关系的影响,在开发的U-DYSAC2有限元程序中进行了相应的数值实施。在试验给定的初边值条件下进行了非饱和渗流模拟分析,并将模拟结果与实测数据比较,表明在压力边界条件反复变化下,考虑滞回效应能获得更接近实测的结果,证实该模型在模拟各种循环变化条件下非饱和土渗流初边值问题的适用性与必要性。对入渗重分布反复变化条件下非饱和土柱流动的数值模拟表明,考虑滞回与不考虑滞回条件下,含水量、孔隙水压力和湿峰的迁移的预测在入渗后的重分布过程差异较大。考虑滞回效应时,土柱上部的脱湿速率、下部的吸湿速率比不考虑滞回时要低。从而证实了非饱和多孔介质中的土水状态依赖于含水量变化,而且强烈依赖于土体的水力路径变化。因此,循环边界条件变化下,毛细滞回效应在非饱和渗流模拟中的影响显著,必须加以考虑。     

考虑气相影响的降雨入渗过程分析研究

降雨入渗过程是水在下渗的过程中驱替空气的水-气二相流过程,对这一过程的精确模拟一直是渗流计算的难点,目前的处理方法通常是忽略孔隙气压力变化的影响。根据多相流理论,结合质量守恒定律和达西定律,建立了水-气二相流模型,模型的求解采用积分有限差分法和Newton-Raphson迭代方法,通过变换主要变量来表达相态的变化,实现了水相、气相边界条件及降雨入渗边界的精确模拟。利用上述模型对一土柱试验进行模拟,从而验证了模型的正确性,研究了一均质土层的降雨入渗过程,得到了孔隙水压力、孔隙气压力和毛细压力及含水率的变化过程。根据入渗率与地表孔隙气压力的变化关系,验证了孔隙气压力的增大对入渗水流产生阻滞作用。在求解非稳定渗流问题中,考虑空气压力变化的影响是值得研究的。     

堤防工程三维随机渗流场理论及研究

给出分析各向异性非均质稳定随机渗流场问题的三维有限元模型;结合实际工程问题,统计分析长江荆南干堤士性参数的分布特征,通过Kolomogorov-Smirnov统计检验表明,渗透系数呈高斯分布假设可以接受;通过对长江荆南干堤随机渗流场的三维有限元统计模拟的数值分析,研究长江荆南干堤渗流场的各种随机特性,并进一步对随机模拟结果进行统计检验,验证模拟结果的合理性;在实际的分析研究中把上下游水位的随机波动引入三维有限元的随机分析模型,分析上下游水位的变异性对渗流场矢量的随机干扰和边界条件的随机性对随机渗流场分析结果变异性的影响。在此基础上进一步考虑施加诸如垂直截渗墙、下游导渗沟等抗渗措施后,它们作为复杂边界条件的扰动,在与场内土性参数的变异性共同影响下,对渗流场水头势分布的随机干扰特性,并与相应的确定性稳定渗流场问题的结果对比,证实随机渗流场研究的必要性、可行性及实用性。实现了对长江荆南干堤的三维渗流场的较为全面的随机场模拟及特性分析,分析得到的结论通过统计检验并结合实测工程数据对照证明是可靠的,所研制的程序是适用的。     

孔隙水压力对岩石力学参数的影响

多孔多相的岩石介质中,孔隙水压力的变化会引起岩石的变形,其力学参数随之改变,而岩石的力学参数变化反过来会引起渗流场发生改变。本文利用岩石力学参数在其中的作用,推导出孔隙水压力变化时,岩石变形与孔隙度、压缩系数和体积弹性模量间的关系式。通过理论计算与分析表明,在渗流场和变形场的非线性程度不很高,参数关系不是很复杂的情况下,利用岩石力学参数的更新来进行耦合分析的方法,对于近期工程应用行之有效。     

考虑渗流条件下开挖面失稳离心试验研究

当土压平衡盾构穿越高水位地层（如穿越江河）时,地下水与土舱之间的高水压差会产生过大的渗透力,导致开挖面失稳。为了研究渗流条件下开挖面失稳问题,开发了一套隧道离心模型试验装置,主要包括刚性模型箱、模型盾构、开挖面伺服加载系统、水位控制系统、储水箱。针对饱和砂质粉土地层,开展了一系列不同水位高度的稳态渗流开挖面失稳模型试验。结果显示,开挖面失稳过程中随着开挖面位移的增加,有效支护压力迅速下降;在达到最小值 之后缓慢回升并趋于稳定;极限有效支护压力 与水头压力 呈线性关系。     

Stability analysis of the downstream slope of a levee

Consideration of the average pore water pressure and average unit weight of soil in the slope of a levee, ignoring the position of phreatic line for stability analysis is not correct. In addition, the influence of a rock-filled portion on the upstream slope, ignored in Kozeny analysis, plays a significant role in the stability of the downstream slope. In this study, levees resting on an impervious base, without a filter and with a filter, taking into consideration a realistic pore water pressure and corresponding seepage force, with unsaturated and saturated unit weights of soil above and below the phreatic line, respectively, by taking into account the capillary effect, location of the inverted filter and rock-filled portion of the upstream slope have been analysed. A numerical technique has been applied to compute the elevation of the phreatic line and a computer program was developed. The appropriate unit weight of soil above and below the phreatic line, taking into consideration the capillary effect, and actual pore water pressure ensure greater safety of the downstream slope. The factor of safety has been found to be 16.5% more for a levee with a filter than a levee without a filter for an upstream pondage depth of 18 m, this is increased by 22.0% in cases where there is an upstream rock-filled portion.     

低渗透煤层气-水流固耦合数学模型及数值模拟

基于有效应力原理,结合弹塑性几何方程、本构方程及平衡方程建立了应力场控制方程。依据流体力学中的质量守恒原理建立了煤层气、水和煤岩体固体颗粒的渗流场方程;以上方程再配以辅助方程和定解条件构成了应力作用下煤层 气-水流固耦合的数学模型,并对简化模型进行了单相流流-固耦合数值模拟分析,得到了压力动态分布曲线,分析了耦合和非耦合情况对压力分布的影响,进一步为煤层气流-固耦合分析提供了一定的理论基础。     

A new parabolic variational inequality formulation of Signorini's condition for non‐steady seepage problems with complex seepage control systems

By extending Darcy's law to the dry domain above the free surface and specifying the boundary condition on the potential seepage surfaces as Signorini's type, a partial differential equation (PDE) defined in the entire domain of interest is formulated for non‐steady seepage flow problems with free surfaces. A new parabolic variational inequality (PVI) formulation equivalent to the PDE formulation is then proposed, in which the flux part of the complementary condition of Signorini's type in the PDE formulation is transformed into natural boundary condition. Consequently, the singularity at the seepage points is eliminated and the difficulty in selecting the trial functions is significantly reduced. By introducing an adaptive penalized Heaviside function in the finite element analysis, the numerical stability of the discrete PVI formulation is well guaranteed. The proposed approach is validated by the existing laboratory tests with sudden rise and dropdown of water heads, and then applied to capture the non‐steady seepage flow behaviors in a homogeneous rectangular dam with five drainage tunnels during a linear dropdown of upstream water head. The non‐steady seepage flow in the surrounding rocks of the underground powerhouse in the Shuibuya Hydropower Project is further modeled, in which a complex seepage control system is involved. Comparisons with the in situ monitoring data show that the calculation results well illustrate the non‐steady seepage flow process during impounding and the operation of the reservoir as well as the seepage control effects of the drainage hole arrays and drainage tunnels. Copyright © 2010 John Wiley & Sons, Ltd.     

库水升降条件下考虑饱和渗透系数空间变异性的白水河滑坡渗流变形分析

堆积层滑坡的岩土体渗透系数具有一定的不确定性,且渗透系数是饱和-非饱和渗流分析的重要参数,开展考虑其空间变异性的库岸堆积层滑坡渗流变形分析具有重要意义。以三峡库区中的白水河滑坡为研究对象,基于地面核磁共振技术获取的岩土体渗透系数,分析滑坡体渗透系数的空间变异特征,采用半变异函数方法求得滑坡体渗透系数的竖直波动范围,在此基础上建立渗透系数的非平稳随机场模型。以非侵入式随机有限元的方式开展库水升降两种工况下不确定模型与确定模型的流固耦合模拟,分析两种模型的渗流场、位移变形特征及其差异。结果表明：相比于确定模型,不确定模型孔压改变的滞后性更为明显,且库水下降工况下整体的变形更大,若忽略滑体渗透系数的非平稳空间变异特征将会低估滑坡的实际变形。