弱透水层瞬态水力参数的地质统计反演识别研究

弱透水层释水变形伴随着水力参数(渗透系数K和贮水率Ss)的瞬态变化，传统的“移动窗口”式反演方法假定每个反演窗口的初始状态是静态的，这与弱透水层释水变形存在明显滞后效应的实际情况不符。本文提出了在高维参数反演优化情境下，时间域上模型参数的动态变化可与空间域上模型参数的一维非均质性等价，因而可用地质统计反演方法可用于识别弱透水层瞬态K和Ss的新思路，实现了该方法在水文地质领域的创新性应用。以上海越流含水层系统为例，运用似线性地质统计反演算法，高精度地拟合了F16分层标处的长观变形数据，获取了原位透镜体状弱透水层K和Ss的瞬态特征。反演结果显示，观测时段内K的变幅为1. 02×10-5~5. 89×10-4 m/a，平均值为1. 38×10-4 m/a，Ss的变幅为4. 34×10-4~6. 02×10-3 m-1，平均值为2. 65×10-3 m-1。基于反演所得参数瞬态特征，分析原位弱透水层在1994年前后分别以弹性和塑性变形为主，弱透水层K呈现与含水层水位相反的季节性变化特征，该现象可能与渗流力作用下细粒物质在含水层与弱透水层透镜体交界面附近的对孔隙的阻塞和疏通作用有关。该方法可获取完整观测时段内水文地质参数的动态变化轨迹，在计算评价地表水 地下水交互领域的河床渗透性、原位修复领域的渗透反应墙渗透性等的动态演化特征方面具有潜在的应用价值，并可作为构建水文地质领域、岩土工程等领域有关本构模型的有效辅助手段。

Characterization of transient aquitard hydraulic parameters using the geostatistical inversion approach