去年冬天在豫西某水库做现场踏勘，均质计坝顶裂缝从桩号0+320一直延伸到0+480，坝渗最宽处能塞进一根手指。透分这座均质坝是析毕七十年代用附近粉质黏土填筑的，运行四十多年，业设浸润线明显抬高了。均质计回来做渗流计算时发现，坝渗下游坡脚出逸点比设计值高出1.7米——这个数字直接解释了裂缝的透分成因。

均质坝的析毕渗透分析，核心就两件事：确定浸润线位置，业设算出渗流量。均质计听着简单，坝渗真动手做毕业设计才发现坑不少。透分达西定律是析毕基础，但实际填筑土料的业设渗透系数往往比实验室测的大一个数量级，因为现场碾压不可能完全均匀，存在层间薄弱面。我用的GeoStudio做有限元分析，模型建了三种工况：正常蓄水位、设计洪水位、骤降工况。正常蓄水位下计算出的单宽渗流量是0.38 m³/d·m，跟现场量水堰观测的0.41 m³/d·m差得不算离谱，但骤降工况就麻烦了——上游水位从正常蓄水位降到死水位，坝体内部孔隙水压力来不及消散，上游坡面附近产生指向坡外的动水压力，安全系数从1.42掉到1.18。

这里得说个容易被忽略的点。很多同学做毕业设计直接用饱和渗透系数，但均质坝上游面在蓄水后其实存在一个非饱和区。我试过用Fredlund-Xing模型拟合土水特征曲线，把非饱和参数也输进去，浸润线位置比纯饱和计算低了约0.3米。这个差异对下游坡面出逸点位置影响不大，但对坝体上部含水量的判断差很多。指导老师看了之后说，实际工程中这个非饱和区对坝顶裂缝的影响比想象中更明显——干湿循环导致土体反复胀缩，裂缝就这么来的。

网格划分也得留个心眼。一开始我用自动划分，单元尺寸0.5米，结果下游坡脚处渗流速度场出现明显锯齿。后来把坡脚附近局部加密到0.1米，流量结果才稳定下来。有限元这东西，网格太粗会抹平关键的水力梯度，网格太细又容易在迭代时震荡。我的经验是：坝体内部0.3米左右，坡面附近0.1米，边界处再适当放疏。时间步长取1天，骤降工况再缩短到0.5天。

渗透系数各向异性不能忽略。碾压填筑的土料，水平方向渗透系数通常是垂直方向的2到5倍。我按3倍取值算了一组，跟各向同性对比，浸润线在下游坡面的出逸点高度差了0.4米。别小看这0.4米，反滤层设计厚度就得跟着调。有个师兄做类似课题时没考虑各向异性，答辩时被评委追问“现场碾压层面怎么处理的”，直接卡壳。

算完浸润线只是第一步，反滤层设计才是见功夫的地方。按《碾压式土石坝设计规范》，反滤料D15要大于被保护土d85的4到5倍。我取的砂砾料D15=1.2mm，被保护土d85=0.18mm，比值6.7，满足要求。但现场取的原状土做颗分试验，d85其实只有0.12mm，比值变成10——太粗了，细颗粒会流失。最后把反滤料D15调到0.6mm才通过。这种细节，不亲手做一遍试验根本体会不到。

渗流计算报告写完那天，我把结果跟现场观测孔水位对了对。坝轴线处计算水位比实测低0.2米，下游坡脚处反而高0.15米。偏差不大，但趋势反了。后来想明白了：计算时假设坝基不透水，实际坝基有一层2米厚的砂卵石，存在深层渗漏。补做坝基渗流分析后，下游水位才跟实测吻合。均质坝的渗透问题，从来不只是坝体本身的事。