确定非线性固结弱透水层水文地质参数的方法

1.本发明涉及水文地质参数测定技术领域，尤其是非线性固结弱透水层水文地质参数。


背景技术：

2.随着社会经济发展和人口数量的激增，水资源短缺和环境恶化成为21世纪人类面临的两个重大问题。地下水资源过度开采引起的地面沉降问题在我国华北平原和长江三角洲地区已普遍存在，其中长三角地区由地下水超采引起的地面沉降大于200mm的面积已占区域面积的十分之一，其中最大累积沉降量达2.80m，造成的经济损失达500亿美元；同时，我国地下水污染呈现出由点到面、由浅到深、由城市到农村的扩展趋势，污染程度日益严重。地下水过度开采和污染相互影响，形成恶性循环。弱透水层的水文地质参数对预测、评价和控制地面沉降至关重要，对地下水资源的评价与管理，以及含水层系统中污染物和热量的迁移转化也非常重要。
3.弱透水层作为多层含水层系统的重要组成部分，在冲积平原和沉积盆地中广泛分布。由于弱透水层的渗透系数一般小于10
–
8m/s，比含水层小几个数量级，弱透水层储存的地下水资源和排水引起的变形量常常被忽视。含水层水文地质参数确定主要通过现场抽水(或注水)试验的方法实现，而弱透水层水文地质参数确定方法的研究还比较少，现存的方法也未考虑弱透水层水文地质参数在固结过程中的变化。目前为止，缺少一种确定非线性固结弱透水层水文地质参数计算的有效原理与方法。


技术实现要素：

4.本发明所要解决的技术问题是提供一种确定非线性固结弱透水层水文地质参数的方法。
5.本发明解决其技术问题所采用的技术方案是：确定非线性固结弱透水层水文地质参数的方法，包括以下步骤：
6.a.建立实验模型：包括筒体，在筒体内从下至上依次设有下反滤层、实验土样和上反滤层，马氏瓶通过进水管与筒体顶部连通，量筒通过出水管与下反滤层连通，在筒体顶部固定有百分表，在进水管和出水管上设有阀门；
7.b.饱和实验土样；
8.c.开展实验，保持实验土样上部水位不变；
9.d.实验数据记录：通过百分表观测弱透水层变形量的时间，在打开出水管阀门后第0.25、0.5、1、1.5、2、3、4、5、6、7、8、9、10、12、14、16、18、20、25、30、35、40、50、60、80、100、120分钟各观测一次并记录变形量s，以后每隔30分钟观测一次；
10.e.持续观测，直至每隔30分钟百分表的读数不变，终止实验；
11.f.参数计算：将记录的变形量s在与标准曲线相同模数的双对数坐标系中作s～t实测曲线，与标准曲线配线，任选一匹配点，记下对应的坐标值s、t和计算固结
系数cw、压缩指数cc和初始渗透系数k(e0)；
[0012][0013]
l为弱透水层厚度，e0弱透水层的初始孔隙比；σ
′0是弱透水层的初始有效应力；γw是水的重度；f(e0)＝1+e0；f(ef)＝1+ef，ef是弱透水层的最终孔隙比。
[0014]
进一步的是，对于a步骤，所述土柱实验模型从下到上依次包括下反滤层、试样段、上反滤层；所述下滤层上端的弱透水层柱体底端侧壁上设有出水管，出水管向上延伸至试样段以上，出水管的出水口对应一接水容器，且出水口处设有阀门；所述上反滤层下底面安装有透水有机玻璃板，且通过不锈钢杆连接百分表；所述土柱试验模型上部设有供水的常水头马氏瓶，且供水水头高度大于出水管出水口的水头高度；所述土柱试验模型最上端设有进水口；所述进水口处设有阀门。
[0015]
进一步的是，对于b步骤，实验土样填充结束后，通过下部出水管缓缓向固结容器加水，至上部接入马氏瓶的出水管溢流。此后静置24小时以上，直至试样厚度稳定。试样内如有气泡，可用真空泵抽出容器内空气。
[0016]
进一步的是，对于c步骤，打开底部阀门，使水自出水管口溢出，流入量筒；上部进水管利用马氏瓶供水，保持试样上部水位不变。
[0017]
本发明的有益效果是：本发明所提供的确定非线性固结弱透水层水文地质参数的方法具有以下优点：
[0018]
①
推导的解析解，理论严密，因此方法具有严格的理论依据；
[0019]
②
实验装置及实验过程简单、易操作，也可直接采取现场原装土样进行实验；
[0020]
③
采用配线法确定参数，方法简单易操作；
[0021]
④
一次实验可以同时求得非线性固结弱透水层的固结系数、压缩指数和初始渗透系数，获取的参数多；
[0022]
⑤
由于变形量测量容易实现且误差小，因此该方法求得的参数精度高，有很好的推广应用价值。
附图说明
[0023]
图1是本发明实验模型的结构示意图；
[0024]
图2是本发明的方法的流程图；
[0025]
图3是弱透水层的无量纲变形量变化图；
[0026]
图4是s～t实测曲线与标准曲线配线图；
[0027]
图5是对应观测时间的弱透水层变形量与实测值对比图；
[0028]
图中零部件、部位及编号：马氏瓶1、进水管2、百分表3、筒体4、上反滤层5、实验土样6、下反滤层7、出水管8、阀门9、量筒10。
具体实施方式
[0029]
下面结合附图对本发明作进一步说明。
[0030]
确定非线性固结弱透水层水文地质参数的方法，首先基于土柱试验模型，推导在
弱透水层主体一侧水头降低某一常量条件下弱透水层变形量公式解析解；然后，给出基于变形量随时间变化的实验资料，采用配线法确定弱透水层固结系数、压缩指数和初始渗透系数的方法。
[0031]
土柱实验模型从下到上依次包括下反滤层7、实验土样6、上反滤层5；所述下反滤层7上端的弱透水层柱体底端侧壁上设有出水管8，出水管8向上延伸至实验土样6以上，出水管8的出水口对应一量筒10，且出水口处设有阀门9；所述上反滤层5下底面安装有透水有机玻璃板，且通过不锈钢杆连接百分表3；所述土柱试验模型上部设有供水的常水头马氏瓶1，且供水水头高度大于出水管8出水口的水头高度。所述土柱试验模型最上端设有进水口；所述进水口处设有阀门9。
[0032]
非线性固结弱透水层参数确定的原理：
[0033]
(1)非线性固结弱透水层水流模型的解析解
[0034]
假设弱透水层柱体饱和承压且测压水头处处相等，弱透水层下侧水头降低某一常量而弱透水层上侧水头不变，测定弱透水层变形量随时间的变化。这时弱透水层为垂向一维流，取如图1坐标系，这样非线性固结弱透水层的水流模型为
[0035][0036]
s(a，0)＝0
ꢀꢀꢀ
0＜a＜l
[0037][0038][0039]
式中，s(a,t)为弱透水层a点t时刻的水头变化；为柱体下侧水头变化值；l为弱透水层厚度；为弱透水层的固结系数；k(e0)为弱透水层的初始渗透系数；cc是弱透水层的压缩指数；e0弱透水层的初始孔隙比；σ
′0是弱透水层的初始有效应力；γw是水的重度；f(e0)＝1+e0；f(ef)＝1+ef；ef是弱透水层的最终孔隙比。
[0040]
对于非线性固结弱透水层水流模型，经变量变换和分离变量法，得解为：
[0041][0042]
式中，是弱透水层底面的最终有效应力。
[0043]
弱透水层变形量的解析解
[0044]
t时刻弱透水层的变形量
[0045][0046]
其中，
[0047][0048]
为无量纲变形量，而
[0049][0050]
为无量纲时间。
[0051]
通过式(3)得弱透水层的无量纲变形量随无量纲时间变化如图3.
[0052]
从弱透水层的无量纲变形量变化图3看出，起初时刻变形量迅速增大，并趋于稳定，当后，变形量趋于定值
[0053]
弱透水层的参数确定方法：
[0054]
对式(3)和式(4)两边同时取对数，有
[0055][0056][0057]
(5)、(6)二式等号右边的第二项都是常数，因此在双对数坐标系内，实验获得的s～t曲线和标准曲线(如图3)在形状上是相同的，只是纵横坐标平移了和cw/l2。采用配线法，将两条曲线重叠，任选一匹配点，记下对应的坐标值s、t和代入(2)、(4)式，得：
[0058]
固结系数：
[0059][0060]
初始渗透系数：
[0061][0062]
压缩指数：
[0063][0064]
如图1所示，实验步骤：
[0065]
(1)调试实验装置，保证实验装置密闭、管路畅通。
[0066]
(2)填充实验土样：把装填弱透水层土样的固结容器擦拭干净，在容器内壁均匀涂抹一层凡士林。测定反滤层(砂土)和现场采取的原状土样(或填充实验土样)的土颗粒密度。填充下反滤层(砂土样)，铺设滤网，在固结容器内放置现场采取的原状土样(或填充实验土样)，然后铺设滤网和多孔有机玻璃板，填充上反滤层(砂土样)，将固结容器用盖板密闭。
[0067]
(3)饱和实验土样：实验土样填充结束后，通过下部出水管缓缓向固结容器加水，至上部接入马氏瓶的出水管溢流。此后静置24小时以上，直至试样厚度稳定。试样内如有气
6m2/s,k(e0)＝1.906
×
10-7m/s和cc＝0.1546。
[0079]
(3)验证
[0080]
将配线法求得的参数带入计算公式(2)，计算得到对应观测时间的弱透水层变形量与实测值对比如图5，所有实测值与对应时刻的计算值十分吻合，表明计算得到的弱透水层水文地质参数较好的反映了弱透水层的固结性质。