不同铺设形式土工膜对土石坝防渗效果
在土石坝运行过程中，上下游水位差形成的渗透动力会驱动水流在坝体及地基岩土体内部发生渗流。水流在渗流过程中会对岩土体颗粒产生渗透力作用，这种力的作用会引起岩土体发生变形，进而改变土体的孔隙率、渗透系数等核心水力参数，同时也会导致渗流边界形状发生改变[6]。而岩土体孔隙结构及边界条件的变化又会反作用于渗流场，使得渗流路径、渗透速度等渗流特性参数发生相应调整，形成'渗流作用-应力变形-渗流特性改变'的动态耦合过程。因此，若忽略渗流场与应力场之间的耦合作用，仅采用单一渗流场分析方法计算得到的结果会与工程实际存在较大偏差，无法真实反映土石坝的渗流特性。基于此，本文在研究土石坝渗流场时，必须考虑应力场与渗流场的双向耦合作用，采用渗流-应力耦合计算方法开展数值模拟分析。
为精准构建渗流-应力耦合计算模型，需基于连续介质力学理论及达西定律，建立耦合控制方程。其中，应力场控制方程采用弹塑性本构关系描述岩土体变形特性，考虑土体的非线性应力-应变关系；渗流场控制方程基于质量守恒定律，结合达西定律描述水流运动规律，并引入孔隙水压力与有效应力的耦合关系[7]。模型参数选取方面，坝体填土及地基岩土体的物理力学参数（如弹性模量、泊松比、黏聚力、内摩擦角）及水力参数（如初始孔隙率、初始渗透系数）均通过室内土工试验及现场原位测试获取，确保参数的真实性与代表性。
通过耦合计算模型分析发现，土工膜的铺设形式对堤坝和地基内的压力水头分布具有显著调控作用：仅在地基铺设土工膜时，由于土工膜阻断了地基内部的竖向渗流通道，地基内的压力水头呈现显著降低趋势，降幅可达30%~50%，但坝体内部渗流路径未发生明显改变，因此堤坝内部的压力水头分布几乎与未铺设土工膜时保持一致；仅在上游坡面铺设土工膜时，土工膜直接阻断了上游坡面的入渗通道，使得坝体内部渗流动力减弱，坝体各深度处的压力水头均有所降低，平均降幅约20%~35%，但由于地基渗流通道未被阻断，地基内的压力水头分布几乎不受影响。
基于上述计算结果，结合不同工程工况需求，可提出土工膜铺设形式的优化设计思路：当工程核心需求为降低堤坝内部压力水头、避免坝体发生渗透变形时，可采用仅在上游面铺设土工膜的方案，该方案施工难度较低且经济性较好；当堤坝坝体填土渗透系数较低（小于1×10⁻⁶ cm/s），而地基岩土体渗透系数较大时，地基渗流是导致工程渗漏的主要原因，此时单独在地基处铺设土工膜是最优选择，可实现以最低成本达到最佳防渗效果的目标；当堤坝坝体填土与地基岩土体的渗透系数均较大（大于1×10⁻⁵ cm/s），即坝体与地基均存在严重渗漏风险时，需采用上游坡面与地基复合铺设土工膜的方案，通过构建'坡面阻渗-地基截渗'的双重防渗体系，实现对总渗流量的有效控制。
