水位恢复法求渗透系数公式
水位恢复法测渗透系数,水位渗透说白了就是恢复往井里灌水,然后看水位怎么降回去。法求这方法比抽水试验省事,系数不用长时间抽水,公式对含水层扰动也小,水位渗透野外应急时特别管用。恢复但很多人只记住了公式,法求没搞懂公式背后的系数物理过程,一换条件就套错。公式
先看最基本的水位渗透公式。假设井是恢复完整井,含水层是法求承压的,稳定注水后突然停注,系数水位恢复过程满足:
K = (Q / 4πT) × ln(t / t′)
这里Q是公式注水流量,T是导水系数,t是从停注开始算的总时间,t′是停注后某一时刻的剩余恢复时间。实际工程里常用简化形式,用直线段斜率求K,即把ln(t/t′)对水位降深S作图,得直线斜率m,则K = 2.303Q / (4πm)。
但别急着套。这公式有个前提——水位恢复必须满足对数线性关系。我见过有人拿非稳定流数据硬拟合,结果K值偏大两倍。原因是含水层边界效应或者井壁阻力影响了早期数据。你画图时,前几分钟的点往往偏离直线,那叫井储效应,得舍弃。
举个例子。去年在河北某粉质粘土场地做注水试验,井径0.1m,注水流量稳定在0.5 L/s,停了之后每30秒记录一次水位。前5个点明显弯曲,我直接删掉,从第6个点开始拟合,R²到了0.998。算出的K是3.2×10⁻⁴ cm/s,跟室内渗透试验结果差不到10%。要是把弯曲段也算进去,K会变成5.8×10⁻⁴,那就离谱了。
还有个坑——潜水含水层和承压含水层公式不一样。潜水井的水位恢复受重力排水影响,恢复曲线不是纯对数型,得用修正公式。有人图省事,拿承压公式硬算潜水井,结果K偏高。实际做法是,当水位恢复高度远小于含水层厚度时(比如小于10%),承压公式还能凑合用;否则得用Boulton延迟排水模型,那计算量就上去了。
公式里的Q怎么取?不是取停注前的瞬时流量,而是取整个注水阶段的平均流量。我见过有人用最后5分钟的流量,如果注水设备压力不稳,那误差就大了。最好在注水过程中每10分钟记一次流量,取算术平均。要是流量波动超过15%,就别用这个法了,换恒定水头试验更靠谱。
最后说个细节。ln(t/t′)里的t′不是从停注开始算的剩余时间,而是从停注时刻到观测时刻的时间差。很多人搞混,把t′当成观测时刻,那对数项就全错了。我习惯在记录表上专门列一列“t′ = 当前时间 - 停注时间”,防止现场手忙脚乱算错。
水位恢复法求K,核心就一句话:画对图,选对段,用对公式。剩下的就是耐心等水位恢复,别急着收仪器。有时候恢复过程要持续几个小时,但数据质量比速度重要。你要是赶时间,至少也得等到恢复量达到总降深的70%以上,不然斜率段太短,拟合出来没意义。









