CN114441738A - 确定非线性固结弱透水层水文地质参数的方法 - Google Patents
确定非线性固结弱透水层水文地质参数的方法 Download PDFInfo
- Publication number
- CN114441738A CN114441738A CN202210086911.5A CN202210086911A CN114441738A CN 114441738 A CN114441738 A CN 114441738A CN 202210086911 A CN202210086911 A CN 202210086911A CN 114441738 A CN114441738 A CN 114441738A
- Authority
- CN
- China
- Prior art keywords
- water
- water outlet
- outlet pipe
- experimental
- consolidated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 95
- 239000002689 soil Substances 0.000 claims abstract description 40
- 238000007596 consolidation process Methods 0.000 claims abstract description 24
- 238000002474 experimental method Methods 0.000 claims abstract description 16
- 230000035699 permeability Effects 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 6
- 239000011800 void material Substances 0.000 claims abstract description 5
- 238000009738 saturating Methods 0.000 claims abstract description 3
- 238000012360 testing method Methods 0.000 claims description 8
- 239000011521 glass Substances 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 3
- 229910001220 stainless steel Inorganic materials 0.000 claims description 3
- 239000010935 stainless steel Substances 0.000 claims description 3
- 239000008400 supply water Substances 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 2
- 230000008859 change Effects 0.000 description 6
- 238000007906 compression Methods 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 230000006835 compression Effects 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000003673 groundwater Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000004579 marble Substances 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035755 proliferation Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 229940099259 vaseline Drugs 0.000 description 1
- 238000003911 water pollution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Dispersion Chemistry (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明涉及水文地质参数测定技术领域,尤其是非线性固结弱透水层水文地质参数;本发明所要解决的技术问题是提供一种确定非线性固结弱透水层水文地质参数的方法。确定非线性固结弱透水层水文地质参数的方法,包括以下步骤:A.建立实验模型;B.饱和实验土样;C.开展实验;D.实验数据记录;E.终止实验;F.参数计算:将记录的变形量S在与标准曲线相同模数的双对数坐标系中作S~t实测曲线,与标准曲线配线,任选一匹配点,记下对应的坐标值S、t和计算固结系数Cw、压缩指数Cc和初始渗透系数K(e0); l为弱透水层厚度,e0弱透水层的初始孔隙比;σ′0是弱透水层的初始有效应力;γw是水的重度;f(e0)=1+e0;f(ef)=1+ef,ef是弱透水层的最终孔隙比。
Description
技术领域
本发明涉及水文地质参数测定技术领域,尤其是非线性固结弱透水层水文地质参数。
背景技术
随着社会经济发展和人口数量的激增,水资源短缺和环境恶化成为21世纪人类面临的两个重大问题。地下水资源过度开采引起的地面沉降问题在我国华北平原和长江三角洲地区已普遍存在,其中长三角地区由地下水超采引起的地面沉降大于200mm的面积已占区域面积的十分之一,其中最大累积沉降量达2.80m,造成的经济损失达500亿美元;同时,我国地下水污染呈现出由点到面、由浅到深、由城市到农村的扩展趋势,污染程度日益严重。地下水过度开采和污染相互影响,形成恶性循环。弱透水层的水文地质参数对预测、评价和控制地面沉降至关重要,对地下水资源的评价与管理,以及含水层系统中污染物和热量的迁移转化也非常重要。
弱透水层作为多层含水层系统的重要组成部分,在冲积平原和沉积盆地中广泛分布。由于弱透水层的渗透系数一般小于10–8m/s,比含水层小几个数量级,弱透水层储存的地下水资源和排水引起的变形量常常被忽视。含水层水文地质参数确定主要通过现场抽水(或注水)试验的方法实现,而弱透水层水文地质参数确定方法的研究还比较少,现存的方法也未考虑弱透水层水文地质参数在固结过程中的变化。目前为止,缺少一种确定非线性固结弱透水层水文地质参数计算的有效原理与方法。
发明内容
本发明所要解决的技术问题是提供一种确定非线性固结弱透水层水文地质参数的方法。
本发明解决其技术问题所采用的技术方案是:确定非线性固结弱透水层水文地质参数的方法,包括以下步骤:
A.建立实验模型:包括筒体,在筒体内从下至上依次设有下反滤层、实验土样和上反滤层,马氏瓶通过进水管与筒体顶部连通,量筒通过出水管与下反滤层连通,在筒体顶部固定有百分表,在进水管和出水管上设有阀门;
B.饱和实验土样;
C.开展实验,保持实验土样上部水位不变;
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分钟观测一次;
E.持续观测,直至每隔30分钟百分表的读数不变,终止实验;
l为弱透水层厚度,e0弱透水层的初始孔隙比;σ′0是弱透水层的初始有效应力;γw是水的重度;f(e0)=1+e0;f(ef)=1+ef,ef是弱透水层的最终孔隙比。
进一步的是,对于A步骤,所述土柱实验模型从下到上依次包括下反滤层、试样段、上反滤层;所述下滤层上端的弱透水层柱体底端侧壁上设有出水管,出水管向上延伸至试样段以上,出水管的出水口对应一接水容器,且出水口处设有阀门;所述上反滤层下底面安装有透水有机玻璃板,且通过不锈钢杆连接百分表;所述土柱试验模型上部设有供水的常水头马氏瓶,且供水水头高度大于出水管出水口的水头高度;所述土柱试验模型最上端设有进水口;所述进水口处设有阀门。
进一步的是,对于B步骤,实验土样填充结束后,通过下部出水管缓缓向固结容器加水,至上部接入马氏瓶的出水管溢流。此后静置24小时以上,直至试样厚度稳定。试样内如有气泡,可用真空泵抽出容器内空气。
进一步的是,对于C步骤,打开底部阀门,使水自出水管口溢出,流入量筒;上部进水管利用马氏瓶供水,保持试样上部水位不变。
本发明的有益效果是:本发明所提供的确定非线性固结弱透水层水文地质参数的方法具有以下优点:
①推导的解析解,理论严密,因此方法具有严格的理论依据;
②实验装置及实验过程简单、易操作,也可直接采取现场原装土样进行实验;
③采用配线法确定参数,方法简单易操作;
④一次实验可以同时求得非线性固结弱透水层的固结系数、压缩指数和初始渗透系数,获取的参数多;
⑤由于变形量测量容易实现且误差小,因此该方法求得的参数精度高,有很好的推广应用价值。
附图说明
图1是本发明实验模型的结构示意图;
图2是本发明的方法的流程图;
图3是弱透水层的无量纲变形量变化图;
图5是对应观测时间的弱透水层变形量与实测值对比图;
图中零部件、部位及编号:马氏瓶1、进水管2、百分表3、筒体4、上反滤层5、实验土样6、下反滤层7、出水管8、阀门9、量筒10。
具体实施方式
下面结合附图对本发明作进一步说明。
确定非线性固结弱透水层水文地质参数的方法,首先基于土柱试验模型,推导在弱透水层主体一侧水头降低某一常量条件下弱透水层变形量公式解析解;然后,给出基于变形量随时间变化的实验资料,采用配线法确定弱透水层固结系数、压缩指数和初始渗透系数的方法。
土柱实验模型从下到上依次包括下反滤层7、实验土样6、上反滤层5;所述下反滤层7上端的弱透水层柱体底端侧壁上设有出水管8,出水管8向上延伸至实验土样6以上,出水管8的出水口对应一量筒10,且出水口处设有阀门9;所述上反滤层5下底面安装有透水有机玻璃板,且通过不锈钢杆连接百分表3;所述土柱试验模型上部设有供水的常水头马氏瓶1,且供水水头高度大于出水管8出水口的水头高度。所述土柱试验模型最上端设有进水口;所述进水口处设有阀门9。
非线性固结弱透水层参数确定的原理:
(1)非线性固结弱透水层水流模型的解析解
假设弱透水层柱体饱和承压且测压水头处处相等,弱透水层下侧水头降低某一常量而弱透水层上侧水头不变,测定弱透水层变形量随时间的变化。这时弱透水层为垂向一维流,取如图1坐标系,这样非线性固结弱透水层的水流模型为
s(a,0)=0 0<a<l
式中,s(a,t)为弱透水层a点t时刻的水头变化;为柱体下侧水头变化值;l为弱透水层厚度;为弱透水层的固结系数;K(e0)为弱透水层的初始渗透系数;Cc是弱透水层的压缩指数;e0弱透水层的初始孔隙比;σ′0是弱透水层的初始有效应力;γw是水的重度;f(e0)=1+e0;f(ef)=1+ef;ef是弱透水层的最终孔隙比。
对于非线性固结弱透水层水流模型,经变量变换和分离变量法,得解为:
弱透水层变形量的解析解
t时刻弱透水层的变形量
其中,
为无量纲变形量,而
为无量纲时间。
通过式(3)得弱透水层的无量纲变形量随无量纲时间变化如图3.
弱透水层的参数确定方法:
对式(3)和式(4)两边同时取对数,有
(5)、(6)二式等号右边的第二项都是常数,因此在双对数坐标系内,实验获得的S~t曲线和标准曲线(如图3)在形状上是相同的,只是纵横坐标平移了和cw/l2。采用配线法,将两条曲线重叠,任选一匹配点,记下对应的坐标值S、t和代入(2)、(4)式,得:
固结系数:
初始渗透系数:
压缩指数:
如图1所示,实验步骤:
(1)调试实验装置,保证实验装置密闭、管路畅通。
(2)填充实验土样:把装填弱透水层土样的固结容器擦拭干净,在容器内壁均匀涂抹一层凡士林。测定反滤层(砂土)和现场采取的原状土样(或填充实验土样)的土颗粒密度。填充下反滤层(砂土样),铺设滤网,在固结容器内放置现场采取的原状土样(或填充实验土样),然后铺设滤网和多孔有机玻璃板,填充上反滤层(砂土样),将固结容器用盖板密闭。
(3)饱和实验土样:实验土样填充结束后,通过下部出水管缓缓向固结容器加水,至上部接入马氏瓶的出水管溢流。此后静置24小时以上,直至试样厚度稳定。试样内如有气泡,可用真空泵抽出容器内空气。
(4)开展实验:打开底部阀门,使水自出水管口溢出,流入量筒。上部进水管利用马氏瓶供水,保持试样上部水位不变。
(5)实验数据记录:通过百分表观测弱透水层变形量的时间,宜在打开出水管阀门后第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分钟各观测一次并记录,以后每隔30分钟观测一次。
(6)持续观测,直至每隔30分钟百分表的读数不变,可终止实验;
(7)参数计算:将记录的变形量S在与标准曲线相同模数的双对数坐标系中作S~t实测曲线,与标准曲线配线,使两条曲线重叠最好,任选一匹配点,记下对应的坐标值S、t和代入公式(7)、(8)、(9)计算固结系数Cw、压缩指数Cc和初始渗透系数K(e0)。
本实施例中一种确定非线性固结弱透水层水文地质参数的方法的具体步骤如下:
(1)实验模型
实验模型主体为有机玻璃圆筒,高100cm,外径40cm,内径38.4cm。圆筒内填充物分为三段,下部填充下反滤层,厚度为22cm,中部填充实验用的土样(粉质黏性土)厚度为24cm,上部填充上反滤层,厚度为16.5cm。上、下反滤层的砂颗粒密度是2.65×103kg/m3,上部反滤层的砂土总量是30kg;实验土样的质量是38.76kg,土颗粒的干密度(ρs)是2.58×103kg/m3,干密度(ρd)是1.38×103kg/m3,实验土样的初始孔隙比(e0)是0.869。开始实验前弱透水层和上、下反滤层的水头相等。开始实验时,将下部反滤层中的水头突然降低100cm,并保持不变。观测弱透水层的变形量,观测资料如表1,持续观测了480分钟。
表1变形量观测资料
t | S | t | S | t | S | t | S |
0.00 | 0 | 1.5 | 0.57 | 10 | 1.66 | 80 | 4.65 |
0.07 | 0.12 | 2 | 0.66 | 12 | 1.85 | 100 | 5.01 |
0.10 | 0.14 | 2.5 | 0.76 | 14 | 2.01 | 120 | 5.01 |
0.13 | 0.16 | 3 | 0.85 | 16 | 2.17 | 150 | 5.24 |
0.17 | 0.17 | 3.5 | 0.92 | 18 | 2.33 | 180 | 5.56 |
0.25 | 0.22 | 4 | 0.99 | 20 | 2.47 | 210 | 5.82 |
0.33 | 0.25 | 4.5 | 1.06 | 25 | 2.78 | 240 | 5.94 |
0.42 | 0.28 | 5 | 1.13 | 30 | 3.06 | 270 | 5.99 |
0.50 | 0.31 | 6 | 1.25 | 35 | 3.33 | 300 | 6.08 |
0.67 | 0.36 | 7 | 1.36 | 40 | 3.56 | 360 | 6.28 |
0.83 | 0.40 | 8 | 1.48 | 50 | 3.98 | 420 | 6.40 |
1.00 | 0.45 | 9 | 1.57 | 60 | 4.32 | 480 | 6.48 |
(2)参数计算
将表中变形量S在双对数坐标系中作S~t实测曲线(图4中圆点),与标准曲线(图4中实线)配线,使两条曲线重叠(如图4),任选一匹配点,记下对应的坐标值S=4.651×10-3m,t=4.8×103s和代入式(7)、(8)、(9)得:cw=1.512×10-6m2/s,K(e0)=1.906×10-7m/s和Cc=0.1546。
(3)验证
将配线法求得的参数带入计算公式(2),计算得到对应观测时间的弱透水层变形量与实测值对比如图5,所有实测值与对应时刻的计算值十分吻合,表明计算得到的弱透水层水文地质参数较好的反映了弱透水层的固结性质。
Claims (4)
1.确定非线性固结弱透水层水文地质参数的方法,其特征在于,包括以下步骤:
A.建立实验模型:包括筒体(4),在筒体(4)内从下至上依次设有下反滤层(7)、实验土样(6)和上反滤层(5),马氏瓶(1)通过进水管(2)与筒体(4)顶部连通,量筒(10)通过出水管(8)与下反滤层(7)连通,在筒体(4)顶部固定有百分表(3),在进水管(2)和出水管(8)上设有阀门(9);
B.饱和实验土样;
C.开展实验,保持实验土样上部水位不变;
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分钟观测一次;
E.持续观测,直至每隔30分钟百分表的读数不变,终止实验;
l为弱透水层厚度,e0弱透水层的初始孔隙比;σ′0是弱透水层的初始有效应力;γw是水的重度;f(e0)=1+e0;f(ef)=1+ef,ef是弱透水层的最终孔隙比。
2.如权利要求1所述的确定非线性固结弱透水层水文地质参数的方法,其特征在于:对于A步骤,所述土柱实验模型从下到上依次包括下反滤层、试样段、上反滤层;所述下滤层上端的弱透水层柱体底端侧壁上设有出水管,出水管向上延伸至试样段以上,出水管的出水口对应一接水容器,且出水口处设有阀门;所述上反滤层下底面安装有透水有机玻璃板,且通过不锈钢杆连接百分表;所述土柱试验模型上部设有供水的常水头马氏瓶,且供水水头高度大于出水管出水口的水头高度;所述土柱试验模型最上端设有进水口;所述进水口处设有阀门。
3.如权利要求1所述的确定非线性固结弱透水层水文地质参数的方法,其特征在于:对于B步骤,实验土样填充结束后,通过下部出水管缓缓向固结容器加水,至上部接入马氏瓶的出水管溢流。此后静置24小时以上,直至试样厚度稳定。试样内如有气泡,可用真空泵抽出容器内空气。
4.如权利要求1所述的确定非线性固结弱透水层水文地质参数的方法,其特征在于:对于C步骤,打开底部阀门,使水自出水管口溢出,流入量筒;上部进水管利用马氏瓶供水,保持试样上部水位不变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210086911.5A CN114441738A (zh) | 2022-01-25 | 2022-01-25 | 确定非线性固结弱透水层水文地质参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210086911.5A CN114441738A (zh) | 2022-01-25 | 2022-01-25 | 确定非线性固结弱透水层水文地质参数的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114441738A true CN114441738A (zh) | 2022-05-06 |
Family
ID=81369540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210086911.5A Pending CN114441738A (zh) | 2022-01-25 | 2022-01-25 | 确定非线性固结弱透水层水文地质参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114441738A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116559047A (zh) * | 2023-05-06 | 2023-08-08 | 中国地质大学(武汉) | 渗透实验装置、方法及渗透系数和流态的评估方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102809642A (zh) * | 2012-08-03 | 2012-12-05 | 河海大学 | 一种确定弱透水层水文地质参数的方法 |
-
2022
- 2022-01-25 CN CN202210086911.5A patent/CN114441738A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102809642A (zh) * | 2012-08-03 | 2012-12-05 | 河海大学 | 一种确定弱透水层水文地质参数的方法 |
Non-Patent Citations (1)
Title |
---|
HUAN LUO ETAL: "An Analytical Method to Calculate Groundwater Released From an Aquitard Undergoing Nonlinear Consolidation", 《WATER RESOURCES RESEARCH》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116559047A (zh) * | 2023-05-06 | 2023-08-08 | 中国地质大学(武汉) | 渗透实验装置、方法及渗透系数和流态的评估方法 |
CN116559047B (zh) * | 2023-05-06 | 2024-01-30 | 中国地质大学(武汉) | 渗透实验装置、方法及渗透系数和流态的评估方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107764713B (zh) | 量程可调式双环入渗装置及土壤渗透系数的原位测试方法 | |
CN102809642B (zh) | 一种确定弱透水层水文地质参数的方法 | |
CN111208042B (zh) | 一种反演非饱和垃圾土水力参数的装置及方法 | |
CN110927038A (zh) | 土体土水特征曲线及渗透系数测试的试验装置及试验方法 | |
CN102141503B (zh) | 雨水花园填料渗透系数的测试装置及测试方法 | |
CN110672497A (zh) | 一种多功能渗透管涌测试仪 | |
CN109342295B (zh) | 一种基于流固耦合分析的透水沥青路面渗透性能评价方法 | |
CN108101341A (zh) | 泥浆快速脱水装置 | |
CN111982699A (zh) | 非饱和土压缩特性及渗透特性试验装置 | |
Boersma | Field measurement of hydraulic conductivity above a water table | |
CN114441738A (zh) | 确定非线性固结弱透水层水文地质参数的方法 | |
CN112949118A (zh) | 一种带帽桩复合地基的沉降计算方法 | |
CN113552037B (zh) | 一种测试垃圾双孔隙度渗流参数的装置及方法 | |
CN111175477A (zh) | 饱和粉细砂层诱导注浆实验模型及实验方法 | |
Bouwer | Estimating the ability of the vadose zone to transmit liquids | |
Zhou et al. | Conceptual data model and method of settlement calculation for deformation and water release from saturated soft soil | |
CN110685659A (zh) | 低渗透介质水力压裂改性的三维模拟系统及使用方法 | |
CN115598040A (zh) | 一种孔隙介质两向渗透系数测定装置及方法 | |
CN112782054B (zh) | 原位土体稳定水文特征参数测定实验装置及其实验方法 | |
CN115032135A (zh) | 用于测定超软土固结参数的水力固结试验装置及试验方法 | |
CN213398523U (zh) | 低渗透介质水力压裂实验装置 | |
CN211206497U (zh) | 加卸荷扰动下非饱和土水气二相流运移实验仪 | |
CN114117836A (zh) | 一种同时确定弱透水层水文地质参数和非达西渗流控制参数的方法 | |
CN111141655A (zh) | 一种基于簇井的土壤透气率真空抽提分层测试方法 | |
CN110579586B (zh) | 加卸荷扰动下非饱和土水气二相流运移实验仪及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220506 |
|
RJ01 | Rejection of invention patent application after publication |