CN106198341A - 降水入渗土柱模拟***及非饱和渗透系数测定方法 - Google Patents
降水入渗土柱模拟***及非饱和渗透系数测定方法 Download PDFInfo
- Publication number
- CN106198341A CN106198341A CN201610480183.0A CN201610480183A CN106198341A CN 106198341 A CN106198341 A CN 106198341A CN 201610480183 A CN201610480183 A CN 201610480183A CN 106198341 A CN106198341 A CN 106198341A
- Authority
- CN
- China
- Prior art keywords
- water
- soil sample
- rainfall
- undisturbed soil
- computer
- 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
- 238000001556 precipitation Methods 0.000 title claims abstract description 108
- 238000001764 infiltration Methods 0.000 title claims abstract description 101
- 230000008595 infiltration Effects 0.000 title claims abstract description 101
- 238000004088 simulation Methods 0.000 title claims abstract description 83
- 238000003556 assay Methods 0.000 title claims abstract description 11
- 239000002689 soil Substances 0.000 claims abstract description 389
- 238000002386 leaching Methods 0.000 claims abstract description 85
- 239000005413 snowmelt Substances 0.000 claims abstract description 47
- 230000000694 effects Effects 0.000 claims abstract description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 358
- 239000000523 sample Substances 0.000 claims description 286
- 238000012544 monitoring process Methods 0.000 claims description 71
- 238000012360 testing method Methods 0.000 claims description 70
- 238000000034 method Methods 0.000 claims description 45
- 239000004033 plastic Substances 0.000 claims description 44
- 229920003023 plastic Polymers 0.000 claims description 44
- HPNSNYBUADCFDR-UHFFFAOYSA-N chromafenozide Chemical compound CC1=CC(C)=CC(C(=O)N(NC(=O)C=2C(=C3CCCOC3=CC=2)C)C(C)(C)C)=C1 HPNSNYBUADCFDR-UHFFFAOYSA-N 0.000 claims description 43
- 230000008569 process Effects 0.000 claims description 41
- 238000002474 experimental method Methods 0.000 claims description 36
- 239000002985 plastic film Substances 0.000 claims description 36
- 230000035939 shock Effects 0.000 claims description 35
- 238000012384 transportation and delivery Methods 0.000 claims description 32
- 229920006395 saturated elastomer Polymers 0.000 claims description 31
- 230000005540 biological transmission Effects 0.000 claims description 26
- 230000008859 change Effects 0.000 claims description 23
- 239000004927 clay Substances 0.000 claims description 17
- 238000004891 communication Methods 0.000 claims description 15
- 229920001971 elastomer Polymers 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000003628 erosive effect Effects 0.000 claims description 9
- 230000003993 interaction Effects 0.000 claims description 9
- 238000000626 liquid-phase infiltration Methods 0.000 claims description 9
- 238000005303 weighing Methods 0.000 claims description 9
- 238000009434 installation Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 8
- 229910000831 Steel Inorganic materials 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000007423 decrease Effects 0.000 claims description 6
- 229910052751 metal Inorganic materials 0.000 claims description 6
- 230000003020 moisturizing effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 239000010959 steel Substances 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims description 6
- 238000004078 waterproofing Methods 0.000 claims description 6
- 230000003247 decreasing effect Effects 0.000 claims description 5
- 230000000149 penetrating effect Effects 0.000 claims description 5
- 230000006735 deficit Effects 0.000 claims description 4
- 239000002196 Pyroceram Substances 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 210000002615 epidermis Anatomy 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 239000010813 municipal solid waste Substances 0.000 claims description 3
- 230000003204 osmotic effect Effects 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 238000005057 refrigeration Methods 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- 238000005527 soil sampling Methods 0.000 claims description 3
- 239000007921 spray Substances 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 238000003780 insertion Methods 0.000 claims 2
- 230000037431 insertion Effects 0.000 claims 2
- 238000004080 punching Methods 0.000 claims 1
- 238000009533 lab test Methods 0.000 abstract description 6
- 238000000205 computational method Methods 0.000 abstract description 3
- 238000012956 testing procedure Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 15
- 230000007547 defect Effects 0.000 description 5
- 239000002184 metal Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004140 cleaning Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 239000008400 supply water Substances 0.000 description 2
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000005194 fractionation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 239000002910 solid waste Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
- G01N5/04—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder
- G01N5/045—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder for determining moisture content
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
降水入渗土柱模拟***及非饱和渗透系数测定方法,模拟***包括底部构件、连接在底部构件上的一个或多个串联的土柱实验标准构件以及土柱实验标准构件顶部的多种形式降水模拟装置;模拟***装土构件为多个标准构件且由计算机自动化控制,并且基于该***实现降水入渗系数和非饱和渗透系数的测定,并给出了测试步骤及计算方法,测定时采用大直径的原状土柱进行室内实验,运用多种形式的降水模拟装置模拟四种水力边界条件,即低强度降雨全入渗边界条件、中强度降雨有径流入渗边界条件、高强度降雨稳定入渗边界条件和融雪入渗边界条件;本发明具有实用性强,使用效果好,便于推广使用的特点。
Description
技术领域
本发明属于非饱和土技术领域,特别涉及一种降水入渗土柱模拟***及非饱和渗透系数测定方法。
背景技术
降水主要是指降雨和降雪,水分以各种形式从大气到达地面,其它形式的降水还包括露、霜、雹等。降水是水文循环的重要环节,也是人类用水的基本来源。降水资料是分析合理洪枯水情、流域旱情的基础,也是水资源的开发利用如防洪、发电、灌溉等的规划设计与管理运用的基础。降水入渗补给地下水的过程是大气水到土壤水到地下水“三水”相互转换关系中最基本的环节之一,降水入渗对地下水的补给量即为降水补给量,它是地下水的主要补给方式,同时,也是区域水均衡计算中的一个重要均衡要素。
土柱实验被广泛应用于农业、林业、地质、土木和环境等研究领域。应用土柱实验可以在实验室内模拟土壤水分和污染物迁移规律。土柱通常分为原状土柱和扰动土柱两种。原状土柱能够用来测试土体本身的结构及其物理性质;当前土柱实验一般采用扰动土柱,扰动土柱是经过筛分形成,或者按照一定比例混合填装所形成的,其不能用来测试土体本身的结构特性。不管采用哪种土柱进行试验,都存在以下问题:①土柱实验装土构件单节长度大,便捷性差,以至于土柱安装困难,对仪器的清洗造成不便;②监测设备安装繁琐,需要对土柱进行钻孔,破坏其结构性;③监测设备一般为人工观测,使得实验精度低、人为性强。
非饱和渗透系数与含水量或基质势的关系是描述非饱和土壤中水分运移和溶质输送的重要函数关系之一,是分析降水条件下土坡稳定性、固体废物填埋场、地下污水的迁移和填土工程等问题的重要参数。由于在非饱和土中有基质吸力的存在,不能用常规的饱和渗透实验方法确定其渗透系数,使得非饱和土渗透系数的确定具有较大的难度。非饱和渗透系数的测定既可在实验室,也可在现场进行。室内实验的优点是可在设定的水力边界条件下测试原状试样和重塑试样的非饱和渗透性能,同时与原位实测相比,室内实验周期短、经济且更成熟。但通常由于非饱和土的成分、结构以及赋存环境的差异,导致了室内实验所用小试件难以代表真正土体的实际情况。
发明内容
为了克服上述现有技术的缺陷,本发明的目的在于提供一种降水入渗土柱模拟***及非饱和渗透系数测定方法,基于大直径的原状土柱进行室内实验,采用多种形式的降水模拟装置,模拟***装土构件为多个标准构件且由计算机自动化控制,并且基于该***实现非饱和渗透系数的测定,并给出了测试步骤及计算方法,具有实用性强,使用效果好,便于推广使用的特点。
为了达到上述目的,本发明的技术方案为:
降水入渗土柱模拟***,包括底部构件(1)、连接在底部构件(1)上的一个或多个串联的土柱实验标准构件(4)以及土柱实验标准构件(4)顶部的多种形式降水模拟装置(C);
所述的底部构件(1)包括位于最下方的底座(1-1),底座(1-1)上的集水点通过塑料软管(2)接入出渗量量杯(3),塑料软管(2)上设置有第三流量传感器(3-1),第三流量传感器(3-1)接入计算机(7);底座(1-1)的上方设置有承力柱(1-2),承力柱(1-2)的上部设置有高进气值陶土板(1-3),高进气值陶土板(1-3)的四周边沿均与底部构件(1)管件(1-7)的内壁水平紧贴,高进气值陶土板(1-3)的顶部设置有滤纸(1-4),滤纸(1-4)的上表面与原状土样(12)接触;管件(1-7)的顶端设置有外螺纹连接段(1-6),外螺纹连接段(1-6)通过法兰(6)与土柱实验标准构件(4)连接;
所述土柱实验标准构件(4)由两个相同的半圆柱体经卡箍(4-30)通过土柱实验标准构件(4)管壁(4-1)的卡箍凹槽(4-3)处连接成一个圆柱体,土柱实验标准构件(4)的管壁(4-1)上设置有圆形小孔(4-4),圆形小孔(4-4)与橡胶塞(5-7)配合使用;多个土柱实验标准构件(4)通过法兰(6)将上下端的螺纹连接段(4-2)进行连接;***件(5)通过橡胶塞(5-7)***原状土样(12)内,***件(5)内传感器所采集的数据都实时传输给计算机(7),土壤热传导吸力探头(4-7)经圆形小孔(4-4)***原状土样(12)内;土壤热传导吸力探头(4-7)内传感器所采集的数据都实时传输给计算机(7);土柱实验标准构件(4)上固定有多个测压管(4-9),多个测压管(4-9)的每个入水口经圆形小孔(4-4)***原状土样(12)内;所述的***件(5)在土柱上按照同一列布置,土壤热传导吸力探头(4-7)在土柱上按照同一列布置,测压管(4-9)在土柱上按照同一列布置;
所述的承力柱(1-2)包括承力柱支座(1-22)以及固定在其上的承力柱主体(1-21),所述承力柱支座(1-22)与底座(1-1)为一体成型,承力柱支座(1-22)在土柱竖向投影按照“一个圆心+以底座(1-1)半径1/2为半径的圆周向五等份”方式布置;承力柱主体(1-21)长度不同使得底座(1-1)呈现坡度;
所述的土柱底部构件(1)、土柱实验标准构件(4)及圆柱构件(10)均由耐高温玻璃钢制成;
所述的***件(5)的最前端设置有***针头(5-6),***件(5)内部管道的转角处设置有橡胶垫片(5-1),内部管道中设置有弹出件(5-2),弹出件(5-2)包括温湿度传感器探头(4-5)或冷热一体金属元件(4-6),弹出件(5-2)后端与导线(9)进行连接,弹出件(5-2)的尾部套设有轻质弹簧(5-4),轻质弹簧(5-4)的末端设置有探针控制器(5-5),***件(5)的外侧中段设置有挡板(5-3);
所述圆形小孔(4-4)形状大小与橡胶塞(5-7)相一致且结合紧密,排布方式为:纵向上相隔排列且遵循上密下疏原则,横向上绕土柱实验标准构件(4)外圆周长8等分排布;
所述的卡箍(4-30)由两半圆环钢圈(4-34)组成并通过一端的铆钉(4-35)铆接,卡箍(4-30)的另一端接头(4-33)通过螺丝杆(4-31)和螺母(4-32)来调节卡箍(4-30)的松紧程度,使标准圆柱构件紧密结合;
所述的法兰(6)内侧设有法兰螺纹(6-1),法兰(6)的两端设置有转动把手(6-2);
所述的计算机(7)的信号端与微控制器(7-0)的信号端相连,微控制器(7-0)设置有温度探测输入端(7-1)、温湿度探测输出端(7-2)、基质吸力探测输出端(7-3);温度探测输入端(7-1)经导线(9)连接冷热一体金属元件(4-6),温湿度探测输出端(7-2)经导线(9)连接温湿度传感器探头(4-5),基质吸力探测输出端(7-3)经导线(9)连接土壤热传导吸力探头(4-7);
所述的多种形式降水模拟装置(C)包括低强度降雨全入渗模拟装置(16)、中强度降雨有径流入渗模拟装置(17)、高强度降雨稳定入渗模拟装置(14)和融雪入渗模拟装置(8);
所述的低强度降雨全入渗模拟装置(16)包括外侧刻有刻度的输水水桶(16-1),输水水桶(16-1)的下方通过输水管(16-12)与降雨喷头(16-9)连通,降雨喷头(16-9)设置于土柱实验标准构件(4)上部的圆柱构件(10)顶部,圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;U型水头控制管(10-2)内的水面上设置有轻质塑料片(16-8),所述轻质塑料片(16-8)形式为薄圆片,所述轻质塑料片(16-8)上方的圆形凹槽内设置有永久磁铁(16-7),永久磁铁(16-7)的正上方设置有拉线(16-3)吊挂的电线圈(16-6),电线圈(16-6)外接有导线(9),拉线(16-3)的上端缠于转轮(16-2)上,转轮(16-2)上设置有把手(16-4),转轮(16-2)上端用拉线(16-3)通过着力构件(16-11)的小孔与止水阀(16-5)下部相连接,止水阀(16-5)下端设置有轻质弹簧(16-10),所述圆柱构件(10)与U型水头控制管(10-2)相连接,轻质塑料片(16-8)、永久磁铁(16-7)和电线圈(16-6)都在U型水头控制管(10-2)的滑槽(16-13)内运动;所述滑槽(16-13)嵌于U型水头控制管(10-2)内壁,所述电线圈(16-6)按照边沿处三等分设置有滑轨(16-14),所述滑轨(16-14)可在滑槽(16-13)内自由上下运动;所述的低强度降雨全入渗模拟装置(16)控制原状土样(12)表面水位高度低于1mm;
所述的中强度降雨有径流入渗模拟装置(17)包括设置在圆柱构件(10)内的降雨容器(17-9)和设置在降雨容器(17-9)底部的圆形降雨孔(17-10),圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;所述降雨容器(17-9)的顶部设置有降雨容器内压力控制管(17-8)和与外部水源连接的进水管(17-7);所述进水管(17-7)上设置有进水电磁阀(17-2)、进水水泵(17-1)和用于对降雨量进行实时监测的第一流量传感器(17-3);所述降雨容器内压力控制管(17-8)上设置有压力控制电磁阀(17-5)和压力传感器(17-6),所述降雨容器内压力控制管(17-8)的尾部连接有空气压缩机(17-4);所述降雨容器(17-9)的顶部内壁上设置有用于对降雨容器(17-9)的水位进行实时监测的水位传感器(17-11);所述圆柱构件(10)内的原状土样(12)上表面外侧开有出水口(10-1),出水口(10-1)通过塑料软管(2)接入径流量量杯(13),所述塑料软管(2)上设置有第二流量传感器(13-1),所述第二流量传感器(13-1)通过导线(9)接入计算机(7);
所述的高强度降雨稳定入渗模拟装置(14)包括外侧刻有刻度的输水水桶(14-1),输水水桶(14-1)的下方通过输水管(14-12)与降雨喷头(14-9)连通,降雨喷头(14-9)设置于土柱实验标准构件(4)上部的圆柱构件(10)顶部,圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;U型水头控制管(10-2)内的水面上设置有轻质塑料片(14-8),所述轻质塑料片(14-8)形式为薄圆片,所述轻质塑料片(14-8)上方的圆形凹槽内设置有永久磁铁(14-7),永久磁铁(14-7)的正上方设置有拉线(14-3)吊挂的电线圈(14-6),电线圈(14-6)外接有导线(9),拉线(14-3)的上端缠于转轮(14-2)上,转轮(14-2)上设置有把手(14-4),转轮(14-2)上端用拉线(14-3)通过着力构件(14-11)的小孔与止水阀(14-5)下部相连接,止水阀(14-5)下端设置有轻质弹簧(14-10),所述圆柱构件(10)与U型水头控制管(10-2)相连接,轻质塑料片(14-8)、永久磁铁(14-7)和电线圈(14-6)都在U型水头控制管(10-2)的滑槽(14-13)内运动;所述滑槽(14-13)嵌于U型水头控制管(10-2)内壁,所述电线圈(14-6)按照边沿处三等分设置有滑轨(14-14),所述滑轨(14-14)可在滑槽(14-13)内自由上下运动;所述的高强度降雨稳定入渗模拟装置(14)控制原状土样(12)表面水位高度高于1cm以上;
所述的融雪入渗模拟装置(8)包括设置在土柱外侧的温度控制器(8-1)以及通过导线(9)与其所连接的调温元件(8-3),所述调温元件(8-3)位于顶盖(8-2)下方,顶盖(8-2)的顶部设置有超声波测距传感器(8-4),顶盖(8-2)位于土柱实验标准构件(4)上部的圆柱构件(10)的上方且紧密接触,圆柱构件(10)内的原状土样(12)上表面外侧开有出水口(10-1),所述出水口(10-1)通过塑料软管(2)接入径流量量杯(13),塑料软管(2)上设置有第二流量传感器(13-1),所述第二流量传感器(13-1)通过导线(9)接入计算机(7),原状土样(12)上设置有碎冰(11)。
基于降水入渗土柱模拟***的非饱和渗透系数测定方法,包括以下步骤:
步骤一、组装土柱实验标准构件
分别对土柱实验标准构件(4)的两块半圆柱体管壁(4-1)进行拼接,对土柱实验标准构件(4)的接缝进行密封及防水处理,然后把卡箍(4-30)套在卡箍凹槽(4-3)上,并通过扳手上紧套在螺丝杆件(4-31)上的螺母(4-32),使卡箍(4-30)牢固地套在卡箍凹槽(4-3)上,然后将多个土柱实验标准构件(4)通过法兰(6)串联实现纵向拼接;
步骤二、安装原状土样
选取预先准备好的原状土样(12),将原状土样(12)竖立在地面上,将已经连接好的多个土柱实验标准构件(4)套住原状土样(12),将原状土样(12)与土柱实验标准构件(4)之间的缝隙采用密封及防水处理保证后续实验时水不从缝隙直接流下;
步骤三、组装土柱实验仪器设备
先将底座(1-1)放置于水平地面上,然后将承力柱主体(1-21)套放在对应的承力柱支座(1-22)上,将高进气值陶土板(1-3)水平放置于承力柱主体(1-21)上方,所述高进气值陶土板(1-3)上表面铺设有滤纸(1-4),所述承力柱主体(1-21)、高进气值陶土板(1-3)、滤纸(1-4)均位于管件(1-7)内部;在集水处最低位置设置一个出水孔外接塑料软管(2),将所述塑料软管(2)的另一端接入出渗量量杯(3),其中所述塑料软管(2)上安装第三流量传感器(3-1),所述第三流量传感器(3-1)通过导线(9)接入计算机(7);将法兰(6)的法兰螺纹(6-1)对准外螺纹连接段(1-6),通过转动把手(6-2)将法兰(6)紧密地安装在底部构件(1)上方,然后将拼接组装而成的土柱实验标准构件(4)连同其套住的原状土样(12)通过土柱实验标准构件(4)的底部螺纹连接段(4-2)与所述底部构件(1)上端的法兰(6)的法兰螺纹(6-1)进行组装,使拼接组装而成的土柱实验标准构件(4)连同其套住的原状土样(12)位于底部构件(1)滤纸(1-4)的正上方,最后通过法兰(6)将圆柱构件(10)进行连接;
步骤四、测定原状土样的初始状态
原状土样(12)的初始含水率:多个温湿度传感器探头(4-5)分别对原状土样(12)的湿度进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的湿度信号,并将各个测试点处原状土样(12)的湿度信号记录为各个测试点处原状土样(12)的初始含水率θc;
原状土样(12)的初始温度:多个温湿度传感器探头(4-5)分别对原状土样(12)的温度进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的温度信号,并将各个测试点处原状土样(12)的温度信号记录为各个测试点处原状土样(12)的初始温度Tc;
原状土样(12)的基质吸力:多个土壤热传导吸力探头(4-7)分别对原状土样(12)的基质吸力进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的基质吸力信号,并将各个测试点处原状土样(12)的基质吸力信号记录为各个测试点处原状土样(12)的初始基质吸力Fac;
原状土样(12)的水头高度:多个测压管(4-9)分别对原状土样(12)的水头高度进行监测得到各个测试点初始阶段对应的水头高度hc;
原状土样(12)的饱和含水率:将原状土样(12)的取样地点带回的其它土样进行饱和含水率测定,作为原状土样(12)的饱和含水率;取土样放入称量盒内,为其注水直至水面浸没土样,浸没10分钟之后将多余的水清除,称质量为m,之后将土样和称量盒放入烘箱内,进行烘干,之后将其置于天平上进行称重得质量为ms,之后利用公式计算求得原状土样(12)的饱和含水率θsat;
步骤五、模拟四种水力边界条件
根据测试要求,多种形式降水模拟装置(C)配合土柱实验标准构件(4)工作来模拟四种水力边界条件,即低强度降雨全入渗边界条件、中强度降雨有径流入渗边界条件、高强度降雨稳定入渗边界条件和融雪入渗边界条件;
步骤六、模拟过程中的各参数监测
多个温湿度传感器探头(4-5)分别对模拟过程中的原状土样(12)的温湿度进行监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的温度、湿度信号,并将各个测试点处原状土样(12)的温度、湿度信号记录为各个测试点记录时刻对应的温度Ti、含水率θi;多个土壤热传导吸力探头(4-7)分别对模拟过程中的原状土样(12)的基质吸力进行监测,将监测到的基质吸力信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的基质吸力信号,并将各个测试点处原状土样(12)的基质吸力信号记录为各个测试点记录时刻对应的基质吸力Fa;多个测压管(4-9)分别对原状土样(12)的水头高度进行监测得到各个测试点记录时刻对应的水头高度hi;
上述所有的监测,其监测时间频率设置如下:降水入渗5分钟内,记录时间间隔为5秒,降水入渗5-15分钟内,记录时间间隔为10秒,降水入渗15-30分钟内,记录时间间隔为15秒,降水入渗30-60分钟内,记录时间间隔为20秒,降水入渗60分钟以后记录时间间隔为60秒,直到实验达到稳定后2-4小时以上;
步骤七、模拟过程中监测结果分析
对所记录的某一时刻的温度Ti、基质吸力Fa监测数据做插值处理,得到某一时刻原状土样(12)所对应的温度云图和基质吸力云图;同时对所记录的某一时刻原状土样(12)所对应的含水率θi监测数据做插值处理,得到某一时刻原状土样(12)对应的含水率云图;
依据含水率云图的变化规律,找出每个时刻所对应的湿润前锋的位置,所述湿润前锋的位置指的是湿润带的边缘,与下部未湿润带之间含水率存在明显突变的部分,各位置连线形成湿润前锋线,从而观察湿润前锋位置随时间t的变化规律;根据含水率云图,依据含水率的大小,找出每个时刻所对应的饱和含水率θsat的等值线,从而确定完全饱和带,所述完全饱和带定义是土柱上表面以下一定深度内出现水分完全饱和的部分;当湿润前锋与完全饱和带在同一监测时刻出现时,所述湿润前锋线与饱和含水率θsat等值线之间的区域被定义为降水入渗过渡带;
绘制基质吸力Fa和含水率θi的关系图,从而分别得到各土层的土-水特征曲线;
步骤八、降水入渗系数计算
根据公式Qr=Qz-Qj,计算得到降水入渗量Qr,单位为cm3;其中Qz为总降水量,单位为cm3;Qj为降水径流量,单位为cm3;
根据公式△S=Qr-Qc,计算得到降水入渗的水分亏损量△S,单位为cm3;其中Qr为降水入渗量,单位为cm3;Qc为降水出渗量,单位为cm3;
根据公式Vr=Qr/t计算得到降水入渗率Vr,单位为cm3/s;其中Qr为降水入渗量,单位为cm3;t为实验测试时间,单位为s;
根据公式Vc=Qc/t计算得到降水出渗率Vc,单位为cm3/s;其中Qc为降水出渗量,单位为cm3;t为实验测试时间,单位为s;
根据公式α=Qc/Qz计算得到降水入渗系数α,单位为无量纲;其中Qc为降水出渗量,单位为cm3;Qz为总降水量,单位为cm3;
步骤九、非饱和渗透系数计算
基于土柱实验的降水入渗,可以概化为一维垂向入渗,其数学模型如下
将Darcy定律代入以上方程(1)得
对于非饱和土,渗透系数k与含水率存在函数关系,所以方程(2)可化为下式
根据Fredlund&Morgenstern所提出的理论,试样所受的法向应力(σ-ua)和基质吸力值的变化将会引起体积含水率θw的变化,即
式中:σ——总应力;
m1 w——与法向应力(σ-ua)变化有关的水的体积变化系数;
——与基质吸力(ua-uw)变化有关的水的体积变化系数;
将方程(4)对时间微分,同时,在非稳定渗流过程中土体单元上并没有外荷载作用,假定在非饱和区气相连续不变,得
其中:即土-水特征曲线的斜率;
由(3)、(5)式可得
其中:k为非饱和渗透系数,单位为无量纲;h为总水头高度,单位为cm;y为土柱测点高度,单位为cm;γw为水的重度,单位为N/m3;为土-水特征曲线的斜率,单位为无量纲;t为时间,单位为s;
依据方程(6)求得任意时刻土柱任意高度处的非饱和渗透系数k,利用插值法绘制其分布图。
步骤五所述的四种水力边界条件,其低强度降雨全入渗的具体模拟步骤为:
(1)、向输水水桶(16-1)内加水(15),此时使止水阀(16-5)处于关闭状态,然后打开止水阀(16-5),水(15)经输水管(16-12)从降雨喷头(16-9)流至圆柱构件(10)内的抗水压冲击板(A)上,再到达原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用;位于U型水头控制管(10-2)内的轻质塑料片(16-8)与其上方的永久磁铁(16-7)在水(15)的浮力作用下上升,达到与圆柱构件(10)内水位高度相同,当上部的圆柱构件(10)内水位上升到高出原状土样(12)上表面1mm时,通过导线(9)给电线圈(16-6)供电形成电磁铁,同时通过把手(16-4)带动转轮(16-2)转动,所述转轮(16-2)上缠有拉线(16-3),所述电线圈(16-6)随着拉线(16-3)在滑槽(16-13)内垂直运动,转动把手(16-4)调整拉线(16-3)长度,使止水阀(16-5)在轻质弹簧(16-10)的作用下接近于关闭状态来进行粗调,然后改变电线圈(16-6)的外部导线(9)所接电源电压的大小来进行精调,使此时止水阀(16-5)刚好处于关闭状态;当U型水头控制管(10-2)内的水位随着上部的圆柱构件(10)内的水位下降时,所述轻质塑料片(16-8)和永久磁铁(16-7)也随之下降,永久磁铁(16-7)与电线圈(16-6)之间的磁力减小,所述轻质弹簧(16-10)受到的荷载增加产生压缩形变,位于轻质弹簧(16-10)上的止水阀(16-5)打开,水(15)继续经输水管(16-12)从降雨喷头(16-9)流入圆柱构件(10)内进行补给;当补给的水位高度再次大于1mm时,止水阀(16-5)就会处于关闭状态,当补给的水位高度再次小于1mm时,止水阀(16-5)打开继续补给,从而模拟出了低强度降雨全入渗边界条件;
(2)、在模拟低强度降雨全入渗过程中,渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在低强度降雨全入渗过程中,观测到输水水桶(16-1)的补水间隔为2小时以上时,说明达到入渗稳定,关闭供水***,停止实验;此时,输水水桶(16-1)内水的损失量与圆柱构件(10)内未渗入的水量之差记录为总降雨量Qz,在低强度降雨全入渗条件下,总降雨量Qz数值上等于入渗量Qr;此时,查看出渗量测量量杯(3)内渗出原状土样(12)水的量,并将该读数记录为低强度降雨全入渗条件下的出渗量Qc。
步骤五所述的四种水力边界条件,其中强度降雨有径流入渗的具体模拟步骤为:
(1)、计算机(7)发出信号接通空气压缩机(17-4)的供电回路,然后驱动压力控制电磁阀(17-5)打开,空气压缩机(17-4)启动,产生的压缩空气经过压力控制电磁阀(17-5)调节压力后通过降雨容器内降雨压力控制管(17-8)进入降雨容器(17-9)内,将降雨容器(17-9)内的水从圆形降雨孔(17-10)中喷出,之后到达圆柱构件(10)内的抗水压冲击板(A)上,再到原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用,模拟出了中强度降雨有径流入渗边界条件;中强度降雨模拟过程中,压力传感器(17-6)对经过降雨容器(17-9)内降雨压力控制管(17-8)的压缩空气压力进行实时监测并将监测到的信号输出给计算机(7),计算机(7)将其接收到的压力监测值与实验设定的压力阈值进行比对,当压力监测值小于压力阈值时,计算机(7)驱动压力控制电磁阀(17-5)打开的开度增大,当压力监测值大于压力阈值时,计算机(7)驱动压力控制电磁阀(17-5)打开的开度减小,从而使经过降雨容器(17-9)内降雨压力控制管(17-8)的压缩空气压力维持在压力阈值范围内;同时,所述水位传感器(17-11)对降雨容器(17-9)内的水位进行实时监测并将监测到的信号输出给计算机(7),计算机(7)将其接收到的水位监测值与实验设定的水位设定值进行比对,当水位监测值小于水位设定值时,计算机(7)控制接通进水水泵(17-1)的供电回路,并驱动进水电磁阀(17-2)打开,外部水源流出的水通过进水管(17-7)流入降雨容器(17-9)内,当水位监测值大于水位设定值时,计算机(7)控制断开进水水泵(17-1)的供电回路,并驱动进水电磁阀(17-2)关闭,停止往降雨容器(17-9)内注水,从而使降雨容器(17-9)内水位维持稳定,保证了降雨过程的持续进行;
(2)、在模拟中强度降雨有径流入渗过程中,未渗入原状土样(12)内的水从出水口(10-1)内流出并经过塑料软管(2)流入径流量量杯(13),渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);模拟中强度降雨过程中,径流量量杯(13)对未渗入原状土样(12)内的水的量进行监测并记录,当相邻的两个采样时刻的未渗入原状土样(12)内的水的量的差值小于等于1cm3时,说明达到了入渗稳定,关闭供水***,停止实验;此时,第一流量传感器(17-3)对降雨量进行实时监测并将监测到的信号实时输出给计算机(7)并记录为总降雨量Qz;查看径流量测量量杯(13)内未渗入原状土样(12)内的水的量,并将该读数记录为降雨径流量Qj;查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为降雨出渗量Qc。
步骤五所述的四种水力边界条件,其高强度降雨稳定入渗的具体模拟步骤为:
(1)、向输水水桶(14-1)内加水(15),此时使止水阀(14-5)处于关闭状态,然后打开止水阀(14-5),所述水(15)经输水管(14-12)从降雨喷头(14-9)流至圆柱构件(10)内的抗水压冲击板(A)上,再到达原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用;位于U型水头控制管(10-2)内的轻质塑料片(14-8)与其上方的永久磁铁(14-7)在水(15)的浮力作用下上升,达到与圆柱构件(10)内水位高度相同,当上部的圆柱构件(10)内水位上升到预设实验高度时,通过导线(9)给电线圈(14-6)供电形成电磁铁,同时通过把手(14-4)带动转轮(14-2)转动,所述转轮(14-2)上缠有拉线(14-3),所述电线圈(14-6)随着拉线(14-3)在滑槽(14-13)内垂直运动,转动把手(14-4)调整拉线(14-3)长度使止水阀(14-5)在轻质弹簧(14-10)的作用下接近于关闭状态来进行粗调,然后改变电线圈(14-6)的外部导线(9)所接电源电压的大小来进行精调,使此时止水阀(14-5)刚好处于关闭状态;当U型水头控制管(10-2)内的水位随着上部的圆柱构件(10)内的水位下降时,所述轻质塑料片(14-8)和永久磁铁(14-7)也随之下降,永久磁铁(14-7)与电线圈(14-6)之间的磁力减小,所述轻质弹簧(14-10)受到的荷载增加产生压缩形变,位于轻质弹簧(14-10)上的止水阀(14-5)打开,水(15)继续经输水管(14-12)从降雨喷头(14-9)流入圆柱构件(10)内进行补给;当补给的水位高度再次大于预设实验高度时,止水阀(14-5)就会处于关闭状态,当补给的水位高度再次小于预设实验高度时,止水阀(14-5)打开继续补给,模拟出了高强度降雨稳定入渗边界条件;
(2)、在模拟高强度降雨稳定入渗过程中,渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在高强度降雨稳定入渗过程中,观测到输水水桶(14-1)的补水间隔为2小时以上时,说明达到入渗稳定;此时,输水水桶(14-1)内水(15)的损失量记录为总降雨量Qz,同时,将输水水桶(14-1)内水(15)的损失量与圆柱构件(10)内未渗入的水的量之差记录为降雨入渗量Qr,查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为降雨出渗量Qc。
步骤五所述的四种水力边界条件,其融雪入渗的具体模拟步骤为:
(1)、计算机(7)发出信号接通温度控制器(8-1)的供电回路,从而控制调温元件(8-3)加热到预设实验温度T1,温度控制器(8-1)对调温元件(8-3)进行实时监测并将监测到的信号实时传输给计算机(7),计算机(7)将其接收到的监测值与调温元件(8-3)预设实验温度T1相比对,当监测值达到调温元件(8-3)预设实验温度T1时,计算机(7)向温度控制器(8-1)发出信号断开调温元件(8-3)的供电回路,调温元件(8-3)停止加热,当监测值低于调温元件(8-3)预设实验温度T1时,计算机(7)经微控制器(7-0)向温度控制器(8-1)发出信号接通调温元件(8-3)的供电回路,调温元件(8-3)开始加热,从而使调温元件(8-3)保持为预设实验温度T1;预设实验温度T1的范围为0℃~80℃;
(2)、计算机(7)发出信号经微控制器(7-0)从温度探测输入端(7-1)经导线(9)到达冷热一体金属元件(4-6),控制冷热一体金属元件(4-6)开始工作,使原状土样(12)达到预设实验温度T2;融雪入渗模拟过程中,埋置于原状土样(12)的温湿度传感器探头(4-5)对附近土壤的温度进行实时监测并将监测到的信号经导线(9)通过微控制器(7-0)的温湿度输出端(7-2)实时传输给计算机(7),计算机(7)将其接收到的温度监测值与冷热一体金属元件(4-6)预设实验温度T2相比对,当温度监测值高于冷热一体金属元件(4-6)预设实验温度T2时,计算机(7)经微控制器(7-0)向温度探测输入端(7-1)发出信号控制冷热一体金属元件(4-6)开始制冷,当温度监测值低于冷热一体金属元件(4-6)实验温度T2时,计算机(7)经微控制器(7-0)向温度探测输入端(7-1)发出信号控制冷热一体金属元件(4-6)开始加热,从而使冷热一体金属元件(4-6)的温度保持为预设实验温度T2,模拟出了融雪低温入渗条件;预设实验温度T2的范围为-20℃~20℃;
(3)、在模拟融雪入渗的过程中,未渗入原状土样(12)内的水从所述出水口(10-1)内流出并经过塑料软管(2)流入径流量量杯(13);渗出原状土样(12)内的水从所述通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在融雪入渗过程中,当相邻的两个采样时刻的未渗入原状土样(12)内的水的量的差值小于等于1cm3时,说明达到了融雪入渗稳定,停止实验;此时,查看径流量测量量杯(13)内未渗入原状土样(12)内的水的量,并将该读数记录为融雪径流量Qj;查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为融雪出渗量Qc;根据原状土样(12)入渗稳定时上表层处的含水率θ,结合原状土样(12)上表层处的初始含水率θc,从而换算出融雪入渗量Qr;融雪入渗量Qr与融雪径流量Qj的和即为总融雪量Qz。
本发明的优点:
1.本发明的土柱实验标准构件(4)由两个相同的半圆柱体连接成一个圆柱体,可以方便地进行原状样和重塑样的实验,同时方便实验仪器的清洗工作,克服目前土柱实验难以采用原状样进行实验以及后期清洗工作复杂的缺陷。
2.橡胶塞(5-7)与***件(5)的配合使用,克服现有土柱实验测试元件***时溢水及水分亏损的缺陷,具有测试结果精确、可靠性高、可重复使用的优点。
3.土柱实验标准构件(4)可根据实验要求多个串联,克服了现有技术中土柱实验中装土部件单节长度大、操作不便的缺陷;且整个***的各部件采用标准构件,使得实验标准化,实验效率大大提高。
4.将模拟***与计算机连接,实现自动化控制,具有精确性、标准性高的特点,克服了现阶段大部分的仪器设备人工操作容易出现误差的缺陷。
5.本发明提出了多种形式降水模拟装置(C),涉及多种强度降雨入渗模拟装置和融雪入渗模拟装置;多种强度降雨入渗模拟装置包括低强度降雨全入渗模拟装置(16)、中强度降雨有径流入渗模拟装置(17)和高强度降雨稳定入渗模拟装置(14);低强度降雨全入渗模拟装置(16)和高强度降雨稳定入渗模拟装置(14)利用永久磁铁和电线圈磁性相斥的原理,从而控制降雨的速度;中强度降雨有径流入渗模拟装置(17)通过空压机实现了降雨强度的控制;融雪入渗模拟装置(8)通过调温元件(8-3)和冷热一体金属元件(4-6)实现融雪低温入渗模拟过程。
6、基于多种形式下的降水入渗土柱模拟***实现非饱和渗透系数的测定,并给出了测试步骤及计算方法,具有实用性强,使用效果好,便于推广使用的特点。
7、本发明基于大直径(不小于50cm)的原状土柱进行室内实验,与传统的室内实验小试件相比,能够更好地代表真正土体的实际情况,与原位实测相比,本发明的大直径原状土柱实验也具备周期短、花费少、技术手段成熟的特点。
8、本发明首次提出多种形式下的降水入渗土柱模拟***及非饱和渗透系数的测定,本发明结构新颖合理,能够有效模拟降水入渗的多种工况,具有准确度高、可靠性强。
附图说明
图1A是基于低强度降雨全入渗模拟装置的***总图。
图1B是低强度降雨全入渗模拟装置(16)的示意图。
图1C是图1B中B处放大图。
图2A是基于中强度降雨有径流模拟装置的***总图。
图2B是中强度降雨有径流入渗模拟装置(17)的示意图。
图3A是基于高强度降雨稳定入渗装置的***总图。
图3B是高强度降雨稳定入渗模拟装置(14)的示意图。
图3C是图3B中B处放大图。
图4A是基于融雪入渗模拟装置的***总图。
图4B是融雪入渗模拟装置(8)的示意图。
图5是土柱实验标准构件(4)与底部构件(1)的装配图。
图6是土柱实验标准构件(4)中管壁(4-1)示意图,其中图6A是主视图,图6B是后视图。
图7是图5中A处的放大图。
图8是底座(1-1)的俯视图。
图9是***件(5)的示意图。
图10是***件(5)与橡胶塞(5-7)的配合示意图。
图11是本发明的卡箍(4-30)的拆分示意图。
图12是法兰(6)的示意图。
图13是本发明降水入渗土柱模拟***总示意图。
具体实施方式
下面参照附图对本发明做详细叙述。
参照图13,降水入渗土柱模拟***,包括底部构件(1)、连接在底部构件(1)上的一个或多个串联的土柱实验标准构件(4)以及土柱实验标准构件(4)顶部的多种形式降水模拟装置(C);所述的多种形式降水模拟装置(C)包括低强度降雨全入渗模拟装置(16)、中强度降雨有径流入渗模拟装置(17)、高强度降雨稳定入渗模拟装置(14)和融雪入渗模拟装置(8);
参照图5,所述的底部构件(1)包括位于最下方的底座(1-1),底座(1-1)上的集水点通过塑料软管(2)接入出渗量量杯(3),塑料软管(2)上设置有第三流量传感器(3-1),第三流量传感器(3-1)接入计算机(7);底座(1-1)的上方设置有承力柱(1-2),承力柱(1-2)的上部设置有高进气值陶土板(1-3),高进气值陶土板(1-3)的四周边沿均与底部构件(1)管件(1-7)的内壁水平紧贴,高进气值陶土板(1-3)的顶部设置有滤纸(1-4),滤纸(1-4)的上表面与原状土样(12)接触;管件(1-7)的顶端设置有外螺纹连接段(1-6),外螺纹连接段(1-6)通过法兰(6)与土柱实验标准构件(4)连接;
参照图6,所述土柱实验标准构件(4)由两个相同的半圆柱体经卡箍(4-30)通过土柱实验标准构件(4)管壁(4-1)的卡箍凹槽(4-3)处连接成一个圆柱体,土柱实验标准构件(4)的管壁(4-1)上设置有圆形小孔(4-4),圆形小孔(4-4)与橡胶塞(5-7)配合使用;多个土柱实验标准构件(4)通过法兰(6)将上下端的螺纹连接段(4-2)进行连接;***件(5)通过橡胶塞(5-7)***原状土样(12)内,***件(5)内传感器所采集的数据都实时传输给计算机(7),土壤热传导吸力探头(4-7)经圆形小孔(4-4)***原状土样(12)内;土壤热传导吸力探头(4-7)内传感器所采集的数据都实时传输给计算机(7);土柱实验标准构件(4)上固定有多个测压管(4-9),多个测压管(4-9)的每个入水口经圆形小孔(4-4)***原状土样(12)内;所述的***件(5)在土柱上按照同一列布置,土壤热传导吸力探头(4-7)在土柱上按照同一列布置,测压管(4-9)在土柱上按照同一列布置;
参照图7、图8,所述的承力柱(1-2)包括承力柱支座(1-22)以及固定在其上的承力柱主体(1-21),所述承力柱支座(1-22)与底座(1-1)为一体成型,承力柱支座(1-22)在土柱竖向投影按照“一个圆心+以底座(1-1)半径1/2为半径的圆周向五等份”方式布置;承力柱主体(1-21)长度不同使得底座(1-1)呈现坡度;
所述的土柱底部构件(1)、土柱实验标准构件(4)及圆柱构件(10)均由耐高温玻璃钢制成;
参照图9、图10,所述的***件(5)的最前端设置有***针头(5-6),***件(5)内部管道的转角处设置有橡胶垫片(5-1),内部管道中设置有弹出件(5-2),弹出件(5-2)包括温湿度传感器探头(4-5)或冷热一体金属元件(4-6),弹出件(5-2)后端与导线(9)进行连接,弹出件(5-2)的尾部套设有轻质弹簧(5-4),轻质弹簧(5-4)的末端设置有探针控制器(5-5),***件(5)的外侧中段设置有挡板(5-3);
所述圆形小孔(4-4)形状大小与橡胶塞(5-7)相一致且结合紧密,排布方式为:纵向上相隔排列且遵循上密下疏原则,横向上绕土柱实验标准构件(4)外圆周长8等分排布;
参照图11,所述的卡箍(4-30)由两半圆环钢圈(4-34)组成并通过一端的铆钉(4-35)铆接,卡箍(4-30)的另一端接头(4-33)通过螺丝杆(4-31)和螺母(4-32)来调节卡箍(4-30)的松紧程度,使标准圆柱构件紧密结合;
参照图12,所述的法兰(6)内侧设有法兰螺纹(6-1),法兰(6)的两端设置有转动把手(6-2);
所述的计算机(7)的信号端与微控制器(7-0)的信号端相连,微控制器(7-0)设置有温度探测输入端(7-1)、温湿度探测输出端(7-2)、基质吸力探测输出端(7-3);温度探测输入端(7-1)经导线(9)连接冷热一体金属元件(4-6),温湿度探测输出端(7-2)经导线(9)连接温湿度传感器探头(4-5),基质吸力探测输出端(7-3)经导线(9)连接土壤热传导吸力探头(4-7)。
实施例一
本实施例中与本发明配套使用的多种强度降雨模拟装置是低强度降雨全入渗模拟装置(16);
参照图1,所述的低强度降雨全入渗模拟装置(16)包括外侧刻有刻度的输水水桶(16-1),输水水桶(16-1)的下方通过输水管(16-12)与降雨喷头(16-9)连通,降雨喷头(16-9)设置于土柱实验标准构件(4)上部的圆柱构件(10)顶部,圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;U型水头控制管(10-2)内的水面上设置有轻质塑料片(16-8),所述轻质塑料片(16-8)形式为薄圆片,所述轻质塑料片(16-8)上方的圆形凹槽内设置有永久磁铁(16-7),永久磁铁(16-7)的正上方设置有拉线(16-3)吊挂的电线圈(16-6),电线圈(16-6)外接有导线(9),拉线(16-3)的上端缠于转轮(16-2)上,转轮(16-2)上设置有把手(16-4),转轮(16-2)上端用拉线(16-3)通过着力构件(16-11)的小孔与止水阀(16-5)下部相连接,止水阀(16-5)下端设置有轻质弹簧(16-10),所述圆柱构件(10)与U型水头控制管(10-2)相连接,轻质塑料片(16-8)、永久磁铁(16-7)和电线圈(16-6)都在U型水头控制管(10-2)的滑槽(16-13)内运动;所述滑槽(16-13)嵌于U型水头控制管(10-2)内壁,所述电线圈(16-6)按照边沿处三等分设置有滑轨(16-14),所述滑轨(16-14)可在滑槽(16-13)内自由上下运动;所述的低强度降雨全入渗模拟装置(16)控制原状土样(12)表面水位高度低于1mm。
本实施例是低强度降雨全入渗模拟***的非饱和渗透系数测定方法包括以下步骤:
步骤一、组装土柱实验标准构件
分别对土柱实验标准构件(4)的两块半圆柱体管壁(4-1)进行拼接,对土柱实验标准构件(4)的接缝进行密封及防水处理,然后把卡箍(4-30)套在卡箍凹槽(4-3)上,并通过扳手上紧套在螺丝杆件(4-31)上的螺母(4-32),使卡箍(4-30)牢固地套在卡箍凹槽(4-3)上,然后将多个土柱实验标准构件(4)通过法兰(6)串联实现纵向拼接;
步骤二、安装原状土样
选取预先准备好的原状土样(12),将原状土样(12)竖立在地面上,将已经连接好的多个土柱实验标准构件(4)套住原状土样(12),将原状土样(12)与土柱实验标准构件(4)之间的缝隙采用密封及防水处理保证后续实验时水不从缝隙直接流下;
步骤三、组装土柱实验仪器设备
先将底座(1-1)放置于水平地面上,然后将承力柱主体(1-21)套放在对应的承力柱支座(1-22)上,将高进气值陶土板(1-3)水平放置于承力柱主体(1-21)上方,所述高进气值陶土板(1-3)上表面铺设有滤纸(1-4),所述承力柱主体(1-21)、高进气值陶土板(1-3)、滤纸(1-4)均位于管件(1-7)内部;在集水处最低位置设置一个出水孔外接塑料软管(2),将所述塑料软管(2)的另一端接入出渗量量杯(3),其中所述塑料软管(2)上安装第三流量传感器(3-1),所述第三流量传感器(3-1)通过导线(9)接入计算机(7);将法兰(6)的法兰螺纹(6-1)对准外螺纹连接段(1-6),通过转动把手(6-2)将法兰(6)紧密地安装在底部构件(1)上方,然后将拼接组装而成的土柱实验标准构件(4)连同其套住的原状土样(12)通过土柱实验标准构件(4)的底部螺纹连接段(4-2)与所述底部构件(1)上端的法兰(6)的法兰螺纹(6-1)进行组装,使拼接组装而成的土柱实验标准构件(4)连同其套住的原状土样(12)位于底部构件(1)滤纸(1-4)的正上方,最后通过法兰(6)将圆柱构件(10)进行连接;
步骤四、测定原状土样的初始状态
原状土样(12)的初始含水率:多个温湿度传感器探头(4-5)分别对原状土样(12)的湿度进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的湿度信号,并将各个测试点处原状土样(12)的湿度信号记录为各个测试点处原状土样(12)的初始含水率θc;
原状土样(12)的初始温度:多个温湿度传感器探头(4-5)分别对原状土样(12)的温度进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的温度信号,并将各个测试点处原状土样(12)的温度信号记录为各个测试点处原状土样(12)的初始温度Tc;
原状土样(12)的基质吸力:多个土壤热传导吸力探头(4-7)分别对原状土样(12)的基质吸力进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的基质吸力信号,并将各个测试点处原状土样(12)的基质吸力信号记录为各个测试点处原状土样(12)的初始基质吸力Fac;
原状土样(12)的水头高度:多个测压管(4-9)分别对原状土样(12)的水头高度进行监测得到各个测试点初始阶段对应的水头高度hc;
原状土样(12)的饱和含水率:将原状土样(12)的取样地点带回的其它土样进行饱和含水率测定,作为原状土样(12)的饱和含水率;取土样放入称量盒内,为其注水直至水面浸没土样,浸没10分钟之后将多余的水清除,称质量为m,之后将土样和称量盒放入烘箱内,进行烘干,之后将其置于天平上进行称重得质量为ms,之后利用公式计算求得原状土样(12)的饱和含水率θsat;
步骤五、模拟低强度降雨全入渗
(1)、向输水水桶(16-1)内加水(15),此时使止水阀(16-5)处于关闭状态,然后打开止水阀(16-5),水(15)经输水管(16-12)从降雨喷头(16-9)流至圆柱构件(10)内的抗水压冲击板(A)上,再到达原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用;位于U型水头控制管(10-2)内的轻质塑料片(16-8)与其上方的永久磁铁(16-7)在水(15)的浮力作用下上升,达到与圆柱构件(10)内水位高度相同,当上部的圆柱构件(10)内水位上升到高出原状土样(12)上表面1mm时,通过导线(9)给电线圈(16-6)供电形成电磁铁,同时通过把手(16-4)带动转轮(16-2)转动,所述转轮(16-2)上缠有拉线(16-3),所述电线圈(16-6)随着拉线(16-3)在滑槽(16-13)内垂直运动,转动把手(16-4)调整拉线(16-3)长度使止水阀(16-5)在轻质弹簧(16-10)的作用下接近于关闭状态来进行粗调,然后改变电线圈(16-6)的外部导线(9)所接电源电压的大小来进行精调,使此时止水阀(16-5)刚好处于关闭状态;当U型水头控制管(10-2)内的水位随着上部的圆柱构件(10)内的水位下降时,所述轻质塑料片(16-8)和永久磁铁(16-7)也随之下降,永久磁铁(16-7)与电线圈(16-6)之间的磁力减小,所述轻质弹簧(16-10)受到的荷载增加产生压缩形变,位于轻质弹簧(16-10)上的止水阀(16-5)打开,水(15)继续经输水管(16-12)从降雨喷头(16-9)流入圆柱构件(10)内进行补给;当补给的水位高度再次大于1mm时,止水阀(16-5)就会处于关闭状态,当补给的水位高度再次小于1mm时,止水阀(16-5)打开继续补给,从而模拟出了低强度降雨全入渗边界条件;
(2)、在模拟低强度降雨全入渗过程中,渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在低强度降雨全入渗过程中,观测到输水水桶(16-1)的补水间隔为2小时以上时,说明达到入渗稳定,关闭供水***,停止实验;此时,输水水桶(16-1)内水的损失量与圆柱构件(10)内未渗入的水量之差记录为总降雨量Qz,在低强度降雨全入渗条件下,总降雨量Qz数值上等于入渗量θi;此时,查看出渗量测量量杯(3)内渗出原状土样(12)水的量,并将该读数记录为低强度降雨全入渗条件下的出渗量θi;
步骤六、模拟过程中的各参数监测
多个温湿度传感器探头(4-5)分别对模拟过程中的原状土样(12)的温湿度进行监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的温度、湿度信号,并将各个测试点处原状土样(12)的温度、湿度信号记录为各个测试点记录时刻对应的温度Ti、含水率θi;多个土壤热传导吸力探头(4-7)分别对模拟过程中的原状土样(12)的基质吸力进行监测,将监测到的基质吸力信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的基质吸力信号,并将各个测试点处原状土样(12)的基质吸力信号记录为各个测试点记录时刻对应的基质吸力Fa;多个测压管(4-9)分别对原状土样(12)的水头高度进行监测得到各个测试点记录时刻对应的水头高度hi;
上述所有的监测,其监测时间频率设置如下:降水入渗5分钟内,记录时间间隔为5秒,降水入渗5-15分钟内,记录时间间隔为10秒,降水入渗15-30分钟内,记录时间间隔为15秒,降水入渗30-60分钟内,记录时间间隔为20秒,降水入渗60分钟以后记录时间间隔为60秒,直到实验达到稳定后2-4小时以上;
步骤七、模拟过程中监测结果分析
对所记录的某一时刻的温度Ti、基质吸力Fa监测数据做插值处理,得到某一时刻原状土样(12)所对应的温度云图和基质吸力云图;同时对所记录的某一时刻原状土样(12)所对应的含水率θi监测数据做插值处理,得到某一时刻原状土样(12)对应的含水率云图;
依据含水率云图的变化规律,找出每个时刻所对应的湿润前锋的位置,所述湿润前锋的位置指的是湿润带的边缘,与下部未湿润带之间含水率存在明显突变的部分,各位置连线形成湿润前锋线,从而观察湿润前锋位置随时间t的变化规律;根据含水率云图,依据含水率的大小,找出每个时刻所对应的饱和含水率θsat的等值线,从而确定完全饱和带,所述完全饱和带定义是土柱上表面以下一定深度内出现水分完全饱和的部分;当湿润前锋与完全饱和带在同一监测时刻出现时,所述湿润前锋线与饱和含水率θsat等值线之间的区域被定义为降水入渗过渡带;
绘制基质吸力Fa和含水率θi的关系图,从而分别得到各土层的土-水特征曲线;
步骤八、降水入渗系数计算
根据公式Qr=Qz-Qj,计算得到降水入渗量Qr,单位为cm3;其中Qz为总降水量,单位为cm3;Qj为降水径流量,单位为cm3;
根据公式△S=Qr-Qc,计算得到降水入渗的水分亏损量△S,单位为cm3;其中Qr为降水入渗量,单位为cm3;Qc为降水出渗量,单位为cm3;
根据公式Vr=Qr/t计算得到降水入渗率Vr,单位为cm3/s;其中Qr为降水入渗量,单位为cm3;t为实验测试时间,单位为s;
根据公式Vc=Qc/t计算得到降水出渗率Vc,单位为cm3/s;其中Qc为降水出渗量,单位为cm3;t为实验测试时间,单位为s;
根据公式α=Qc/Qz计算得到降水入渗系数α,单位为无量纲;其中Qc为降水出渗量,单位为cm3;Qz为总降水量,单位为cm3;
步骤九、非饱和渗透系数计算
基于土柱实验的降水入渗,可以概化为一维垂向入渗,其数学模型如下
将Darcy定律代入以上方程(1)得
对于非饱和土,渗透系数k与含水率存在函数关系,所以方程(2)可化为下式
根据Fredlund&Morgenstern所提出的理论,试样所受的法向应力(σ-ua)和基质吸力值的变化将会引起体积含水率θw的变化,即
式中:σ——总应力;
m1 w——与法向应力(σ-ua)变化有关的水的体积变化系数;
——与基质吸力(ua-uw)变化有关的水的体积变化系数;
将方程(4)对时间微分,同时,在非稳定渗流过程中土体单元上并没有外荷载作用,假定在非饱和区气相连续不变,得
其中:即土-水特征曲线的斜率;
由(3)、(5)式可得
其中:k为非饱和渗透系数,单位为无量纲;h为总水头高度,单位为cm;y为土柱测点高度,单位为cm;γw为水的重度,单位为N/m3;为土-水特征曲线的斜率,单位为无量纲;t为时间,单位为s;
依据方程(6)求得任意时刻土柱任意高度处的非饱和渗透系数k,利用插值法绘制其分布图。
实施例二
本实施例中与本发明配套使用的多种强度降雨模拟装置是中强度降雨有径流入渗模拟装置(17);
参照图2,所述的中强度降雨有径流入渗模拟装置(17)包括设置在圆柱构件(10)内的降雨容器(17-9)和设置在降雨容器(17-9)底部的圆形降雨孔(17-10),圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;所述降雨容器(17-9)的顶部设置有降雨容器内压力控制管(17-8)和与外部水源连接的进水管(17-7);所述进水管(17-7)上设置有进水电磁阀(17-2)、进水水泵(17-1)和用于对降雨量进行实时监测的第一流量传感器(17-3);所述降雨容器内压力控制管(17-8)上设置有压力控制电磁阀(17-5)和压力传感器(17-6),所述降雨容器内压力控制管(17-8)的尾部连接有空气压缩机(17-4);所述降雨容器(17-9)的顶部内壁上设置有用于对降雨容器(17-9)的水位进行实时监测的水位传感器(17-11);所述圆柱构件(10)内的原状土样(12)上表面外侧开有出水口(10-1),出水口(10-1)通过塑料软管(2)接入径流量量杯(13),所述塑料软管(2)上设置有第二流量传感器(13-1),所述第二流量传感器(13-1)通过导线(9)接入计算机(7)。
本实施例是基于中强度降雨有径流模拟***的非饱和渗透系数测定方法包括以下步骤:
步骤一、组装土柱实验标准构件
具体步骤同实施例一;
步骤二、安装原状土样
具体步骤同实施例一;
步骤三、组装土柱实验仪器设备
具体步骤同实施例一;
步骤四、测定原状土样的初始状态
具体步骤同实施例一;
步骤五、模拟中强度降雨有径流入渗
(1)、计算机(7)发出信号接通空气压缩机(17-4)的供电回路,然后驱动压力控制电磁阀(17-5)打开,空气压缩机(17-4)启动,产生的压缩空气经过压力控制电磁阀(17-5)调节压力后通过降雨容器内降雨压力控制管(17-8)进入降雨容器(17-9)内,将降雨容器(17-9)内的水从圆形降雨孔(17-10)中喷出,之后到达圆柱构件(10)内的抗水压冲击板(A)上,再到原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用,模拟出了中强度降雨有径流入渗边界条件;中强度降雨模拟过程中,压力传感器(17-6)对经过降雨容器(17-9)内降雨压力控制管(17-8)的压缩空气压力进行实时监测并将监测到的信号输出给计算机(7),计算机(7)将其接收到的压力监测值与实验设定的压力阈值进行比对,当压力监测值小于压力阈值时,计算机(7)驱动压力控制电磁阀(17-5)打开的开度增大,当压力监测值大于压力阈值时,计算机(7)驱动压力控制电磁阀(17-5)打开的开度减小,从而使经过降雨容器(17-9)内降雨压力控制管(17-8)的压缩空气压力维持在压力阈值范围内;同时,所述水位传感器(17-11)对降雨容器(17-9)内的水位进行实时监测并将监测到的信号输出给计算机(7),计算机(7)将其接收到的水位监测值与实验设定的水位设定值进行比对,当水位监测值小于水位设定值时,计算机(7)控制接通进水水泵(17-1)的供电回路,并驱动进水电磁阀(17-2)打开,外部水源流出的水通过进水管(17-7)流入降雨容器(17-9)内,当水位监测值大于水位设定值时,计算机(7)控制断开进水水泵(17-1)的供电回路,并驱动进水电磁阀(17-2)关闭,停止往降雨容器(17-9)内注水,从而使降雨容器(17-9)内水位维持稳定,保证了降雨过程的持续进行;
(2)、在模拟中强度降雨有径流入渗过程中,未渗入原状土样(12)内的水从出水口(10-1)内流出并经过塑料软管(2)流入径流量量杯(13),渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);模拟中强度降雨过程中,径流量量杯(13)对未渗入原状土样(12)内的水的量进行监测并记录,当相邻的两个采样时刻的未渗入原状土样(12)内的水的量的差值小于等于1cm3时,说明达到了入渗稳定,关闭供水***,停止实验;此时,第一流量传感器(17-3)对降雨量进行实时监测并将监测到的信号实时输出给计算机(7)并记录为总降雨量Qz;查看径流量测量量杯(13)内未渗入原状土样(12)内的水的量,并将该读数记录为降雨径流量Qj;查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为降雨出渗量Qc;
步骤六、模拟过程中的各参数监测
具体步骤同实施例一;
步骤七、模拟过程中监测结果分析
具体步骤同实施例一;
步骤八、降水入渗系数计算
具体步骤同实施例一;
步骤九、非饱和渗透系数计算
具体步骤同实施例一。
实施例三
本实施例中与本发明配套使用的多种强度降雨模拟装置是高强度降雨稳定入渗模拟装置(14);
参照图3,所述的高强度降雨稳定入渗模拟装置(14)包括外侧刻有刻度的输水水桶(14-1),输水水桶(14-1)的下方通过输水管(14-12)与降雨喷头(14-9)连通,降雨喷头(14-9)设置于土柱实验标准构件(4)上部的圆柱构件(10)顶部,圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;U型水头控制管(10-2)内的水面上设置有轻质塑料片(14-8),所述轻质塑料片(14-8)形式为薄圆片,所述轻质塑料片(14-8)上方的圆形凹槽内设置有永久磁铁(14-7),永久磁铁(14-7)的正上方设置有拉线(14-3)吊挂的电线圈(14-6),电线圈(14-6)外接有导线(9),拉线(14-3)的上端缠于转轮(14-2)上,转轮(14-2)上设置有把手(14-4),转轮(14-2)上端用拉线(14-3)通过着力构件(14-11)的小孔与止水阀(14-5)下部相连接,止水阀(14-5)下端设置有轻质弹簧(14-10),所述圆柱构件(10)与U型水头控制管(10-2)相连接,轻质塑料片(14-8)、永久磁铁(14-7)和电线圈(14-6)都在U型水头控制管(10-2)的滑槽(14-13)内运动;所述滑槽(14-13)嵌于U型水头控制管(10-2)内壁,所述电线圈(14-6)按照边沿处三等分设置有滑轨(14-14),所述滑轨(14-14)可在滑槽(14-13)内自由上下运动;所述的高强度降雨稳定入渗模拟装置(14)控制原状土样(12)表面水位高度高于1cm以上。
本实施例是基于高强度降雨稳定模拟***的非饱和渗透系数测定方法包括以下步骤:
步骤一、组装土柱实验标准构件
具体步骤同实施例一;
步骤二、安装原状土样
具体步骤同实施例一;
步骤三、组装土柱实验仪器设备
具体步骤同实施例一;
步骤四、测定原状土样的初始状态
具体步骤同实施例一;
步骤五、模拟高强度降雨稳定入渗
(1)、向输水水桶(14-1)内加水(15),此时使止水阀(14-5)处于关闭状态,然后打开止水阀(14-5),所述水(15)经输水管(14-12)从降雨喷头(14-9)流至圆柱构件(10)内的抗水压冲击板(A)上,再到达原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用;位于U型水头控制管(10-2)内的轻质塑料片(14-8)与其上方的永久磁铁(14-7)在水(15)的浮力作用下上升,达到与圆柱构件(10)内水位高度相同,当上部的圆柱构件(10)内水位上升到预设实验高度时,通过导线(9)给电线圈(14-6)供电形成电磁铁,同时通过把手(14-4)带动转轮(14-2)转动,所述转轮(14-2)上缠有拉线(14-3),所述电线圈(14-6)随着拉线(14-3)在滑槽(14-13)内垂直运动,转动把手(14-4)调整拉线(14-3)长度使止水阀(14-5)在轻质弹簧(14-10)的作用下接近于关闭状态来进行粗调,然后改变电线圈(14-6)的外部导线(9)所接电源电压的大小来进行精调,使此时止水阀(14-5)刚好处于关闭状态;当U型水头控制管(10-2)内的水位随着上部的圆柱构件(10)内的水位下降时,所述轻质塑料片(14-8)和永久磁铁(14-7)也随之下降,永久磁铁(14-7)与电线圈(14-6)之间的磁力减小,所述轻质弹簧(14-10)受到的荷载增加产生压缩形变,位于轻质弹簧(14-10)上的止水阀(14-5)打开,水(15)继续经输水管(14-12)从降雨喷头(14-9)流入圆柱构件(10)内进行补给;当补给的水位高度再次大于预设实验高度时,止水阀(14-5)就会处于关闭状态,当补给的水位高度再次小于预设实验高度时,止水阀(14-5)打开继续补给,模拟出了高强度降雨稳定入渗边界条件;
(2)、在模拟高强度降雨稳定入渗过程中,渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在高强度降雨稳定入渗过程中,观测到输水水桶(14-1)的补水间隔为2小时以上时,说明达到入渗稳定;此时,输水水桶(14-1)内水(15)的损失量记录为总降雨量Qz,同时,将输水水桶(14-1)内水(15)的损失量与圆柱构件(10)内未渗入的水的量之差记录为降雨入渗量Qr,查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为降雨出渗量Qc;
步骤六、模拟过程中的各参数监测
具体步骤同实施例一;
步骤七、模拟过程中监测结果分析
具体步骤同实施例一;
步骤八、降水入渗系数计算
具体步骤同实施例一;
步骤九、非饱和渗透系数计算
具体步骤同实施例一。
实施例四
本实施例中与本发明配套使用的多种强度降雨模拟装置是融雪入渗模拟装置(8);
参照图4,所述的融雪入渗模拟装置(8)包括设置在土柱外侧的温度控制器(8-1)以及通过导线(9)与其所连接的调温元件(8-3),所述调温元件(8-3)位于顶盖(8-2)下方,顶盖(8-2)的顶部设置有超声波测距传感器(8-4),顶盖(8-2)位于土柱实验标准构件(4)上部的圆柱构件(10)的上方且紧密接触,圆柱构件(10)内的原状土样(12)上表面外侧开有出水口(10-1),所述出水口(10-1)通过塑料软管(2)接入径流量量杯(13),塑料软管(2)上设置有第二流量传感器(13-1),所述第二流量传感器(13-1)通过导线(9)接入计算机(7),原状土样(12)上设置有碎冰(11)。
本实施例是基于融雪入渗模拟***的非饱和渗透系数测定方法包括以下步骤:
步骤一、组装土柱实验标准构件
具体步骤同实施例一;
步骤二、安装原状土样
具体步骤同实施例一;
步骤三、组装土柱实验仪器设备
具体步骤同实施例一;
步骤四、测定原状土样的初始状态
具体步骤同实施例一;
步骤五、模拟模拟融雪入渗
(1)、计算机(7)发出信号接通温度控制器(8-1)的供电回路,从而控制调温元件(8-3)加热到预设实验温度T1,温度控制器(8-1)对调温元件(8-3)进行实时监测并将监测到的信号实时传输给计算机(7),计算机(7)将其接收到的监测值与调温元件(8-3)预设实验温度T1相比对,当监测值达到调温元件(8-3)预设实验温度T1时,计算机(7)向温度控制器(8-1)发出信号断开调温元件(8-3)的供电回路,调温元件(8-3)停止加热,当监测值低于调温元件(8-3)预设实验温度T1时,计算机(7)经微控制器(7-0)向温度控制器(8-1)发出信号接通调温元件(8-3)的供电回路,调温元件(8-3)开始加热,从而使调温元件(8-3)保持为预设实验温度T1;预设实验温度T1的范围为0℃~80℃;
(2)、计算机(7)发出信号经微控制器(7-0)从温度探测输入端(7-1)经导线(9)到达冷热一体金属元件(4-6),控制冷热一体金属元件(4-6)开始工作,使原状土样(12)达到预设实验温度T2;融雪入渗模拟过程中,埋置于原状土样(12)的温湿度传感器探头(4-5)对附近土壤的温度进行实时监测并将监测到的信号经导线(9)通过微控制器(7-0)的温湿度输出端(7-2)实时传输给计算机(7),计算机(7)将其接收到的温度监测值与冷热一体金属元件(4-6)预设实验温度T2相比对,当温度监测值高于冷热一体金属元件(4-6)预设实验温度T2时,计算机(7)经微控制器(7-0)向温度探测输入端(7-1)发出信号控制冷热一体金属元件(4-6)开始制冷,当温度监测值低于冷热一体金属元件(4-6)实验温度T2时,计算机(7)经微控制器(7-0)向温度探测输入端(7-1)发出信号控制冷热一体金属元件(4-6)开始加热,从而使冷热一体金属元件(4-6)的温度保持为预设实验温度T2,模拟出了融雪低温入渗条件;预设实验温度T2的范围为-20℃~20℃;
(3)、在模拟融雪入渗的过程中,未渗入原状土样(12)内的水从所述出水口(10-1)内流出并经过塑料软管(2)流入径流量量杯(13);渗出原状土样(12)内的水从所述通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在融雪入渗过程中,当相邻的两个采样时刻的未渗入原状土样(12)内的水的量的差值小于等于1cm3时,说明达到了融雪入渗稳定,停止实验;此时,查看径流量测量量杯(13)内未渗入原状土样(12)内的水的量,并将该读数记录为融雪径流量Qj;查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为融雪出渗量Qc;根据原状土样(12)入渗稳定时上表层处的含水率θ,结合原状土样(12)上表层处的初始含水率θc,从而换算出融雪入渗量Qr;融雪入渗量Qr与融雪径流量Qj的和即为总融雪量Qz;
步骤六、模拟过程中的各参数监测
具体步骤同实施例一;
步骤七、模拟过程中监测结果分析
具体步骤同实施例一;
步骤八、降水入渗系数计算
具体步骤同实施例一;
步骤九、非饱和渗透系数计算
具体步骤同实施例一。
Claims (10)
1.降水入渗土柱模拟***,其特征在于,包括底部构件(1)、连接在底部构件(1)上的一个或多个串联的土柱实验标准构件(4)以及土柱实验标准构件(4)顶部的多种形式降水模拟装置(C);
所述的底部构件(1)包括位于最下方的底座(1-1),底座(1-1)上的集水点通过塑料软管(2)接入出渗量量杯(3),塑料软管(2)上设置有第三流量传感器(3-1),第三流量传感器(3-1)接入计算机(7);底座(1-1)的上方设置有承力柱(1-2),承力柱(1-2)的上部设置有高进气值陶土板(1-3),高进气值陶土板(1-3)的四周边沿均与底部构件(1)管件(1-7)的内壁水平紧贴,高进气值陶土板(1-3)的顶部设置有滤纸(1-4),滤纸(1-4)的上表面与原状土样(12)接触;管件(1-7)的顶端设置有外螺纹连接段(1-6),外螺纹连接段(1-6)通过法兰(6)与土柱实验标准构件(4)连接;
所述土柱实验标准构件(4)由两个相同的半圆柱体经卡箍(4-30)通过土柱实验标准构件(4)管壁(4-1)的卡箍凹槽(4-3)处连接成一个圆柱体,土柱实验标准构件(4)的管壁(4-1)上设置有圆形小孔(4-4),圆形小孔(4-4)与橡胶塞(5-7)配合使用;多个土柱实验标准构件(4)通过法兰(6)将上下端的螺纹连接段(4-2)进行连接;***件(5)通过橡胶塞(5-7)***原状土样(12)内,***件(5)内传感器所采集的数据都实时传输给计算机(7),土壤热传导吸力探头(4-7)经圆形小孔(4-4)***原状土样(12)内;土壤热传导吸力探头(4-7)内传感器所采集的数据都实时传输给计算机(7);土柱实验标准构件(4)上固定有多个测压管(4-9),多个测压管(4-9)的每个入水口经圆形小孔(4-4)***原状土样(12)内;所述的***件(5)在土柱上按照同一列布置,土壤热传导吸力探头(4-7)在土柱上按照同一列布置,测压管(4-9)在土柱上按照同一列布置。
2.根据权利要求1所述的降水入渗土柱模拟***,其特征在于,所述的承力柱(1-2)包括承力柱支座(1-22)以及固定在其上的承力柱主体(1-21),所述承力柱支座(1-22)与底座(1-1)为一体成型,承力柱支座(1-22)在土柱竖向投影按照“一个圆心+以底座(1-1)半径1/2为半径的圆周向五等份”方式布置;承力柱主体(1-21)长度不同使得底座(1-1)呈现坡度;
所述的土柱底部构件(1)、土柱实验标准构件(4)及圆柱构件(10)均由耐高温玻璃钢制成。
3.根据权利要求1所述的降水入渗土柱模拟***,其特征在于,所述的***件(5)的最前端设置有***针头(5-6),***件(5)内部管道的转角处设置有橡胶垫片(5-1),内部管道中设置有弹出件(5-2),弹出件(5-2)包括温湿度传感器探头(4-5)或冷热一体金属元件(4-6),弹出件(5-2)后端与导线(9)进行连接,弹出件(5-2)的尾部套设有轻质弹簧(5-4),轻质弹簧(5-4)的末端设置有探针控制器(5-5),***件(5)的外侧中段设置有挡板(5-3)。
4.根据权利要求1所述的降水入渗土柱模拟***,其特征在于,所述圆形小孔(4-4)形状大小与橡胶塞(5-7)相一致且结合紧密,排布方式为:纵向上相隔排列且遵循上密下疏原则,横向上绕土柱实验标准构件(4)外圆周长8等分排布。
5.根据权利要求1所述的降水入渗土柱模拟***,其特征在于,所述的卡箍(4-30)由两半圆环钢圈(4-34)组成并通过一端的铆钉(4-35)铆接,卡箍(4-30)的另一端接头(4-33)通过螺丝杆(4-31)和螺母(4-32)来调节卡箍(4-30)的松紧程度,使标准圆柱构件紧密结合;
所述的法兰(6)内侧设有法兰螺纹(6-1),法兰(6)的两端设置有转动把手(6-2)。
6.根据权利要求1所述的降水入渗土柱模拟***,其特征在于,所述的计算机(7)的信号端与微控制器(7-0)的信号端相连,微控制器(7-0)设置有温度探测输入端(7-1)、温湿度探测输出端(7-2)、基质吸力探测输出端(7-3);温度探测输入端(7-1)经导线(9)连接冷热一体金属元件(4-6),温湿度探测输出端(7-2)经导线(9)连接温湿度传感器探头(4-5),基质吸力探测输出端(7-3)经导线(9)连接土壤热传导吸力探头(4-7)。
7.根据权利要求1所述的降水入渗土柱模拟***,其特征在于,所述的多种形式降水模拟装置(C)包括低强度降雨全入渗模拟装置(16)、中强度降雨有径流入渗模拟装置(17)、高强度降雨稳定入渗模拟装置(14)和融雪入渗模拟装置(8);
所述的低强度降雨全入渗模拟装置(16)包括外侧刻有刻度的输水水桶(16-1),输水水桶(16-1)的下方通过输水管(16-12)与降雨喷头(16-9)连通,降雨喷头(16-9)设置于土柱实验标准构件(4)上部的圆柱构件(10)顶部,圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;U型水头控制管(10-2)内的水面上设置有轻质塑料片(16-8),所述轻质塑料片(16-8)形式为薄圆片,所述轻质塑料片(16-8)上方的圆形凹槽内设置有永久磁铁(16-7),永久磁铁(16-7)的正上方设置有拉线(16-3)吊挂的电线圈(16-6),电线圈(16-6)外接有导线(9),拉线(16-3)的上端缠于转轮(16-2)上,转轮(16-2)上设置有把手(16-4),转轮(16-2)上端用拉线(16-3)通过着力构件(16-11)的小孔与止水阀(16-5)下部相连接,止水阀(16-5)下端设置有轻质弹簧(16-10),所述圆柱构件(10)与U型水头控制管(10-2)相连接,轻质塑料片(16-8)、永久磁铁(16-7)和电线圈(16-6)都在U型水头控制管(10-2)的滑槽(16-13)内运动;所述滑槽(16-13)嵌于U型水头控制管(10-2)内壁,所述电线圈(16-6)按照边沿处三等分设置有滑轨(16-14),所述滑轨(16-14)可在滑槽(16-13)内自由上下运动;所述的低强度降雨全入渗模拟装置(16)控制原状土样(12)表面水位高度低于1mm;
所述的中强度降雨有径流入渗模拟装置(17)包括设置在圆柱构件(10)内的降雨容器(17-9)和设置在降雨容器(17-9)底部的圆形降雨孔(17-10),圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;所述降雨容器(17-9)的顶部设置有降雨容器内压力控制管(17-8)和与外部水源连接的进水管(17-7);所述进水管(17-7)上设置有进水电磁阀(17-2)、进水水泵(17-1)和用于对降雨量进行实时监测的第一流量传感器(17-3);所述降雨容器内压力控制管(17-8)上设置有压力控制电磁阀(17-5)和压力传感器(17-6),所述降雨容器内压力控制管(17-8)的尾部连接有空气压缩机(17-4);所述降雨容器(17-9)的顶部内壁上设置有用于对降雨容器(17-9)的水位进行实时监测的水位传感器(17-11);所述圆柱构件(10)内的原状土样(12)上表面外侧开有出水口(10-1),出水口(10-1)通过塑料软管(2)接入径流量量杯(13),所述塑料软管(2)上设置有第二流量传感器(13-1),所述第二流量传感器(13-1)通过导线(9)接入计算机(7);
所述的高强度降雨稳定入渗模拟装置(14)包括外侧刻有刻度的输水水桶(14-1),输水水桶(14-1)的下方通过输水管(14-12)与降雨喷头(14-9)连通,降雨喷头(14-9)设置于土柱实验标准构件(4)上部的圆柱构件(10)顶部,圆柱构件(10)内高出原状土样(12)上表面2-5厘米处设置有抗水压冲击板(A),抗水压冲击板(A)上设置有筛孔;U型水头控制管(10-2)内的水面上设置有轻质塑料片(14-8),所述轻质塑料片(14-8)形式为薄圆片,所述轻质塑料片(14-8)上方的圆形凹槽内设置有永久磁铁(14-7),永久磁铁(14-7)的正上方设置有拉线(14-3)吊挂的电线圈(14-6),电线圈(14-6)外接有导线(9),拉线(14-3)的上端缠于转轮(14-2)上,转轮(14-2)上设置有把手(14-4),转轮(14-2)上端用拉线(14-3)通过着力构件(14-11)的小孔与止水阀(14-5)下部相连接,止水阀(14-5)下端设置有轻质弹簧(14-10),所述圆柱构件(10)与U型水头控制管(10-2)相连接,轻质塑料片(14-8)、永久磁铁(14-7)和电线圈(14-6)都在U型水头控制管(10-2)的滑槽(14-13)内运动;所述滑槽(14-13)嵌于U型水头控制管(10-2)内壁,所述电线圈(14-6)按照边沿处三等分设置有滑轨(14-14),所述滑轨(14-14)可在滑槽(14-13)内自由上下运动;所述的高强度降雨稳定入渗模拟装置(14)控制原状土样(12)表面水位高度高于1cm以上。
8.根据权利要求7所述的降水入渗土柱模拟***,其特征在于,所述的融雪入渗模拟装置(8)包括设置在土柱外侧的温度控制器(8-1)以及通过导线(9)与其所连接的调温元件(8-3),所述调温元件(8-3)位于顶盖(8-2)下方,顶盖(8-2)的顶部设置有超声波测距传感器(8-4),顶盖(8-2)位于土柱实验标准构件(4)上部的圆柱构件(10)的上方且紧密接触,圆柱构件(10)内的原状土样(12)上表面外侧开有出水口(10-1),所述出水口(10-1)通过塑料软管(2)接入径流量量杯(13),塑料软管(2)上设置有第二流量传感器(13-1),所述第二流量传感器(13-1)通过导线(9)接入计算机(7),原状土样(12)上设置有碎冰(11)。
9.基于上述任一权利要求所述的降水入渗土柱模拟***的非饱和渗透系数测定方法,其特征在于,包括以下步骤:
步骤一、组装土柱实验标准构件
分别对土柱实验标准构件(4)的两块半圆柱体管壁(4-1)进行拼接,对土柱实验标准构件(4)的接缝进行密封及防水处理,然后把卡箍(4-30)套在卡箍凹槽(4-3)上,并通过扳手上紧套在螺丝杆件(4-31)上的螺母(4-32),使卡箍(4-30)牢固地套在卡箍凹槽(4-3)上,然后将多个土柱实验标准构件(4)通过法兰(6)串联实现纵向拼接;
步骤二、安装原状土样
选取预先准备好的原状土样(12),将原状土样(12)竖立在地面上,将已经连接好的多个土柱实验标准构件(4)套住原状土样(12),将原状土样(12)与土柱实验标准构件(4)之间的缝隙采用密封及防水处理保证后续实验时水不从缝隙直接流下;
步骤三、组装土柱实验仪器设备
先将底座(1-1)放置于水平地面上,然后将承力柱主体(1-21)套放在对应的承力柱支座(1-22)上,将高进气值陶土板(1-3)水平放置于承力柱主体(1-21)上方,所述高进气值陶土板(1-3)上表面铺设有滤纸(1-4),所述承力柱主体(1-21)、高进气值陶土板(1-3)、滤纸(1-4)均位于管件(1-7)内部;在集水处最低位置设置一个出水孔外接塑料软管(2),将所述塑料软管(2)的另一端接入出渗量量杯(3),其中所述塑料软管(2)上安装第三流量传感器(3-1),所述第三流量传感器(3-1)通过导线(9)接入计算机(7);将法兰(6)的法兰螺纹(6-1)对准外螺纹连接段(1-6),通过转动把手(6-2)将法兰(6)紧密地安装在底部构件(1)上方,然后将拼接组装而成的土柱实验标准构件(4)连同其套住的原状土样(12)通过土柱实验标准构件(4)的底部螺纹连接段(4-2)与所述底部构件(1)上端的法兰(6)的法兰螺纹(6-1)进行组装,使拼接组装而成的土柱实验标准构件(4)连同其套住的原状土样(12)位于底部构件(1)滤纸(1-4)的正上方,最后通过法兰(6)将圆柱构件(10)进行连接;
步骤四、测定原状土样的初始状态
原状土样(12)的初始含水率:多个温湿度传感器探头(4-5)分别对原状土样(12)的湿度进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的湿度信号,并将各个测试点处原状土样(12)的湿度信号记录为各个测试点处原状土样(12)的初始含水率θc;
原状土样(12)的初始温度:多个温湿度传感器探头(4-5)分别对原状土样(12)的温度进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的温度信号,并将各个测试点处原状土样(12)的温度信号记录为各个测试点处原状土样(12)的初始温度Tc;
原状土样(12)的基质吸力:多个土壤热传导吸力探头(4-7)分别对原状土样(12)的基质吸力进行一次监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的基质吸力信号,并将各个测试点处原状土样(12)的基质吸力信号记录为各个测试点处原状土样(12)的初始基质吸力Fac;
原状土样(12)的水头高度:多个测压管(4-9)分别对原状土样(12)的水头高度进行监测得到各个测试点初始阶段对应的水头高度hc;
原状土样(12)的饱和含水率:将原状土样(12)的取样地点带回的其它土样进行饱和含水率测定,作为原状土样(12)的饱和含水率;取土样放入称量盒内,为其注水直至水面浸没土样,浸没10分钟之后将多余的水清除,称质量为m,之后将土样和称量盒放入烘箱内,进行烘干,之后将其置于天平上进行称重得质量为ms,之后利用公式计算求得原状土样(12)的饱和含水率θsat;
步骤五、模拟四种水力边界条件
根据测试要求,多种形式降水模拟装置(C)配合土柱实验标准构件(4)工作来模拟四种水力边界条件,即低强度降雨全入渗边界条件、中强度降雨有径流入渗边界条件、高强度降雨稳定入渗边界条件和融雪入渗边界条件;
步骤六、模拟过程中的各参数监测
多个温湿度传感器探头(4-5)分别对模拟过程中的原状土样(12)的温湿度进行监测,并将监测到的信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的温度、湿度信号,并将各个测试点处原状土样(12)的温度、湿度信号记录为各个测试点记录时刻对应的温度Ti、含水率θi;多个土壤热传导吸力探头(4-7)分别对模拟过程中的原状土样(12)的基质吸力进行监测,将监测到的基质吸力信号传给微控制器(7-0),微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7),计算机(7)接收并记录多个测试点处原状土样(12)的基质吸力信号,并将各个测试点处原状土样(12)的基质吸力信号记录为各个测试点记录时刻对应的基质吸力Fa;多个测压管(4-9)分别对原状土样(12)的水头高度进行监测得到各个测试点记录时刻对应的水头高度hi;
上述所有的监测,其监测时间频率设置如下:降水入渗5分钟内,记录时间间隔为5秒,降水入渗5-15分钟内,记录时间间隔为10秒,降水入渗15-30分钟内,记录时间间隔为15秒,降水入渗30-60分钟内,记录时间间隔为20秒,降水入渗60分钟以后记录时间间隔为60秒,直到实验达到稳定后2-4小时以上;
步骤七、模拟过程中监测结果分析
对所记录的某一时刻的温度Ti、基质吸力Fa监测数据做插值处理,得到某一时刻原状土样(12)所对应的温度云图和基质吸力云图;同时对所记录的某一时刻原状土样(12)所对应的含水率θi监测数据做插值处理,得到某一时刻原状土样(12)对应的含水率云图;
依据含水率云图的变化规律,找出每个时刻所对应的湿润前锋的位置,所述湿润前锋的位置指的是湿润带的边缘,与下部未湿润带之间含水率存在明显突变的部分,各位置连线形成湿润前锋线,从而观察湿润前锋位置随时间t的变化规律;根据含水率云图,依据含水率的大小,找出每个时刻所对应的饱和含水率θsat的等值线,从而确定完全饱和带,所述完全饱和带定义是土柱上表面以下一定深度内出现水分完全饱和的部分;当湿润前锋与完全饱和带在同一监测时刻出现时,所述湿润前锋线与饱和含水率θsat等值线之间的区域被定义为降水入渗过渡带;
绘制基质吸力Fa和含水率θi的关系图,从而分别得到各土层的土-水特征曲线;
步骤八、降水入渗系数计算
根据公式Qr=Qz-Qj,计算得到降水入渗量Qr,单位为cm3;其中Qz为总降水量,单位为cm3;Qj为降水径流量,单位为cm3;
根据公式△S=Qr-Qc,计算得到降水入渗的水分亏损量△S,单位为cm3;其中Qr为降水入渗量,单位为cm3;Qc为降水出渗量,单位为cm3;
根据公式Vr=Qr/t计算得到降水入渗率Vr,单位为cm3/s;其中Qr为降水入渗量,单位为cm3;t为实验测试时间,单位为s;
根据公式Vc=Qc/t计算得到降水出渗率Vc,单位为cm3/s;其中Qc为降水出渗量,单位为cm3;t为实验测试时间,单位为s;
根据公式α=Qc/Qz计算得到降水入渗系数α,单位为无量纲;其中Qc为降水出渗量,单位为cm3;Qz为总降水量,单位为cm3;
步骤九、非饱和渗透系数计算
基于土柱实验的降水入渗,可以概化为一维垂向入渗,其数学模型如下
将Darcy定律代入以上方程(1)得
对于非饱和土,渗透系数k与含水率存在函数关系,所以方程(2)可化为下式
根据Fredlund&Morgenstern所提出的理论,试样所受的法向应力(σ-ua)和基质吸力值的变化将会引起体积含水率θw的变化,即
式中:σ——总应力;
m1 w——与法向应力(σ-ua)变化有关的水的体积变化系数;
——与基质吸力(ua-uw)变化有关的水的体积变化系数;
将方程(4)对时间微分,同时,在非稳定渗流过程中土体单元上并没有外荷载作用,假定在非饱和区气相连续不变,得
其中:即土-水特征曲线的斜率;
由(3)、(5)式可得
其中:k为非饱和渗透系数,单位为无量纲;h为总水头高度,单位为cm;y为土柱测点高度,单位为cm;γw为水的重度,单位为N/m3;为土-水特征曲线的斜率,单位为无量纲;t为时间,单位为s;
依据方程(6)求得任意时刻土柱任意高度处的非饱和渗透系数k,利用插值法绘制其分布图。
10.根据权利要求9所述的降水入渗土柱模拟***的非饱和渗透系数测定方法,其特征在于:
步骤五所述的四种水力边界条件,其低强度降雨全入渗的具体模拟步骤为:
(1)、向输水水桶(16-1)内加水(15),此时使止水阀(16-5)处于关闭状态,然后打开止水阀(16-5),水(15)经输水管(16-12)从降雨喷头(16-9)流至圆柱构件(10)内的抗水压冲击板(A)上,再到达原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用;位于U型水头控制管(10-2)内的轻质塑料片(16-8)与其上方的永久磁铁(16-7)在水(15)的浮力作用下上升,达到与圆柱构件(10)内水位高度相同,当上部的圆柱构件(10)内水位上升到高出原状土样(12)上表面1mm时,通过导线(9)给电线圈(16-6)供电形成电磁铁,同时通过把手(16-4)带动转轮(16-2)转动,所述转轮(16-2)上缠有拉线(16-3),所述电线圈(16-6)随着拉线(16-3)在滑槽(16-13)内垂直运动,转动把手(16-4)调整拉线(16-3)长度使止水阀(16-5)在轻质弹簧(16-10)的作用下接近于关闭状态来进行粗调,然后改变电线圈(16-6)的外部导线(9)所接电源电压的大小来进行精调,使此时止水阀(16-5)刚好处于关闭状态;当U型水头控制管(10-2)内的水位随着上部的圆柱构件(10)内的水位下降时,所述轻质塑料片(16-8)和永久磁铁(16-7)也随之下降,永久磁铁(16-7)与电线圈(16-6)之间的磁力减小,所述轻质弹簧(16-10)受到的荷载增加产生压缩形变,位于轻质弹簧(16-10)上的止水阀(16-5)打开,水(15)继续经输水管(16-12)从降雨喷头(16-9)流入圆柱构件(10)内进行补给;当补给的水位高度再次大于1mm时,止水阀(16-5)就会处于关闭状态,当补给的水位高度再次小于1mm时,止水阀(16-5)打开继续补给,从而模拟出了低强度降雨全入渗边界条件;
(2)、在模拟低强度降雨全入渗过程中,渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在低强度降雨全入渗过程中,观测到输水水桶(16-1)的补水间隔为2小时以上时,说明达到入渗稳定,关闭供水***,停止实验;此时,输水水桶(16-1)内水的损失量与圆柱构件(10)内未渗入的水量之差记录为总降雨量Qz,在低强度降雨全入渗条件下,总降雨量Qz数值上等于入渗量Qr;此时,查看出渗量测量量杯(3)内渗出原状土样(12)水的量,并将该读数记录为低强度降雨全入渗条件下的出渗量Qc;
步骤五所述的四种水力边界条件,其中强度降雨有径流入渗的具体模拟步骤为:
(1)、计算机(7)发出信号接通空气压缩机(17-4)的供电回路,然后驱动压力控制电磁阀(17-5)打开,空气压缩机(17-4)启动,产生的压缩空气经过压力控制电磁阀(17-5)调节压力后通过降雨容器内降雨压力控制管(17-8)进入降雨容器(17-9)内,将降雨容器(17-9)内的水从圆形降雨孔(17-10)中喷出,之后到达圆柱构件(10)内的抗水压冲击板(A)上,再到原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用,模拟出了中强度降雨有径流入渗边界条件;中强度降雨模拟过程中,压力传感器(17-6)对经过降雨容器(17-9)内降雨压力控制管(17-8)的压缩空气压力进行实时监测并将监测到的信号输出给计算机(7),计算机(7)将其接收到的压力监测值与实验设定的压力阈值进行比对,当压力监测值小于压力阈值时,计算机(7)驱动压力控制电磁阀(17-5)打开的开度增大,当压力监测值大于压力阈值时,计算机(7)驱动压力控制电磁阀(17-5)打开的开度减小,从而使经过降雨容器(17-9)内降雨压力控制管(17-8)的压缩空气压力维持在压力阈值范围内;同时,所述水位传感器(17-11)对降雨容器(17-9)内的水位进行实时监测并将监测到的信号输出给计算机(7),计算机(7)将其接收到的水位监测值与实验设定的水位设定值进行比对,当水位监测值小于水位设定值时,计算机(7)控制接通进水水泵(17-1)的供电回路,并驱动进水电磁阀(17-2)打开,外部水源流出的水通过进水管(17-7)流入降雨容器(17-9)内,当水位监测值大于水位设定值时,计算机(7)控制断开进水水泵(17-1)的供电回路,并驱动进水电磁阀(17-2)关闭,停止往降雨容器(17-9)内注水,从而使降雨容器(17-9)内水位维持稳定,保证了降雨过程的持续进行;
(2)、在模拟中强度降雨有径流入渗过程中,未渗入原状土样(12)内的水从出水口(10-1)内流出并经过塑料软管(2)流入径流量量杯(13),渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);模拟中强度降雨过程中,径流量量杯(13)对未渗入原状土样(12)内的水的量进行监测并记录,当相邻的两个采样时刻的未渗入原状土样(12)内的水的量的差值小于等于1cm3时,说明达到了入渗稳定,关闭供水***,停止实验;此时,第一流量传感器(17-3)对降雨量进行实时监测并将监测到的信号实时输出给计算机(7)并记录为总降雨量Qz;查看径流量测量量杯(13)内未渗入原状土样(12)内的水的量,并将该读数记录为降雨径流量Qj;查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为降雨出渗量Qc;
步骤五所述的四种水力边界条件,其高强度降雨稳定入渗的具体模拟步骤为:
(1)、向输水水桶(14-1)内加水(15),此时使止水阀(14-5)处于关闭状态,然后打开止水阀(14-5),所述水(15)经输水管(14-12)从降雨喷头(14-9)流至圆柱构件(10)内的抗水压冲击板(A)上,再到达原状土样(12)的上表面,所述抗水压冲击板(A)起到保护原状土样(12)不受降雨侵蚀的作用;位于U型水头控制管(10-2)内的轻质塑料片(14-8)与其上方的永久磁铁(14-7)在水(15)的浮力作用下上升,达到与圆柱构件(10)内水位高度相同,当上部的圆柱构件(10)内水位上升到预设实验高度时,通过导线(9)给电线圈(14-6)供电形成电磁铁,同时通过把手(14-4)带动转轮(14-2)转动,所述转轮(14-2)上缠有拉线(14-3),所述电线圈(14-6)随着拉线(14-3)在滑槽(14-13)内垂直运动,转动把手(14-4)调整拉线(14-3)长度使止水阀(14-5)在轻质弹簧(14-10)的作用下接近于关闭状态来进行粗调,然后改变电线圈(14-6)的外部导线(9)所接电源电压的大小来进行精调,使此时止水阀(14-5)刚好处于关闭状态;当U型水头控制管(10-2)内的水位随着上部的圆柱构件(10)内的水位下降时,所述轻质塑料片(14-8)和永久磁铁(14-7)也随之下降,永久磁铁(14-7)与电线圈(14-6)之间的磁力减小,所述轻质弹簧(14-10)受到的荷载增加产生压缩形变,位于轻质弹簧(14-10)上的止水阀(14-5)打开,水(15)继续经输水管(14-12)从降雨喷头(14-9)流入圆柱构件(10)内进行补给;当补给的水位高度再次大于预设实验高度时,止水阀(14-5)就会处于关闭状态,当补给的水位高度再次小于预设实验高度时,止水阀(14-5)打开继续补给,模拟出了高强度降雨稳定入渗边界条件;
(2)、在模拟高强度降雨稳定入渗过程中,渗出原状土样(12)内的水从通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在高强度降雨稳定入渗过程中,观测到输水水桶(14-1)的补水间隔为2小时以上时,说明达到入渗稳定,关闭供水***,停止实验;此时,输水水桶(14-1)内水(15)的损失量记录为总降雨量Qz,同时,将输水水桶(14-1)内水(15)的损失量与圆柱构件(10)内未渗入的水的量之差记录为降雨入渗量Qr,查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为降雨出渗量Qc;
步骤五所述的四种水力边界条件,其融雪入渗的具体模拟步骤为
(1)、计算机(7)发出信号接通温度控制器(8-1)的供电回路,从而控制调温元件(8-3)加热到预设实验温度T1,温度控制器(8-1)对调温元件(8-3)进行实时监测并将监测到的信号实时传输给计算机(7),计算机(7)将其接收到的监测值与调温元件(8-3)预设实验温度T1相比对,当监测值达到调温元件(8-3)预设实验温度T1时,计算机(7)向温度控制器(8-1)发出信号断开调温元件(8-3)的供电回路,调温元件(8-3)停止加热,当监测值低于调温元件(8-3)预设实验温度T1时,计算机(7)经微控制器(7-0)向温度控制器(8-1)发出信号接通调温元件(8-3)的供电回路,调温元件(8-3)开始加热,从而使调温元件(8-3)保持为预设实验温度T1;预设实验温度T1的范围为0℃~80℃;
(2)、计算机(7)发出信号经微控制器(7-0)从温度探测输入端(7-1)经导线(9)到达冷热一体金属元件(4-6),控制冷热一体金属元件(4-6)开始工作,使原状土样(12)达到预设实验温度T2;融雪入渗模拟过程中,埋置于原状土样(12)的温湿度传感器探头(4-5)对附近土壤的温度进行实时监测并将监测到的信号经导线(9)通过微控制器(7-0)的温湿度输出端(7-2)实时传输给计算机(7),计算机(7)将其接收到的温度监测值与冷热一体金属元件(4-6)预设实验温度T2相比对,当温度监测值高于冷热一体金属元件(4-6)预设实验温度T2时,计算机(7)经微控制器(7-0)向温度探测输入端(7-1)发出信号控制冷热一体金属元件(4-6)开始制冷,当温度监测值低于冷热一体金属元件(4-6)实验温度T2时,计算机(7)经微控制器(7-0)向温度探测输入端(7-1)发出信号控制冷热一体金属元件(4-6)开始加热,从而使冷热一体金属元件(4-6)的温度保持为预设实验温度T2,模拟出了融雪低温入渗条件;预设实验温度T2的范围为-20℃~20℃;
(3)、在模拟融雪入渗的过程中,未渗入原状土样(12)内的水从所述出水口(10-1)内流出并经过塑料软管(2)流入径流量量杯(13);渗出原状土样(12)内的水从所述通水孔(1-5)内流出并经过塑料软管(2)流入出渗量测量量杯(3);在融雪入渗过程中,当相邻的两个采样时刻的未渗入原状土样(12)内的水的量的差值小于等于1cm3时,说明达到了融雪入渗稳定,停止实验;此时,查看径流量测量量杯(13)内未渗入原状土样(12)内的水的量,并将该读数记录为融雪径流量Qj;查看出渗量测量量杯(3)内渗出原状土样(12)的水的量,并将该读数记录为融雪出渗量Qc;根据原状土样(12)入渗稳定时上表层处的含水率θ,结合原状土样(12)上表层处的初始含水率θc,从而换算出融雪入渗量Qr;融雪入渗量Qr与融雪径流量Qj的和即为总融雪量Qz。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610480183.0A CN106198341A (zh) | 2016-06-25 | 2016-06-25 | 降水入渗土柱模拟***及非饱和渗透系数测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610480183.0A CN106198341A (zh) | 2016-06-25 | 2016-06-25 | 降水入渗土柱模拟***及非饱和渗透系数测定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106198341A true CN106198341A (zh) | 2016-12-07 |
Family
ID=57461159
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610480183.0A Pending CN106198341A (zh) | 2016-06-25 | 2016-06-25 | 降水入渗土柱模拟***及非饱和渗透系数测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106198341A (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106908368A (zh) * | 2017-03-14 | 2017-06-30 | 上海理工大学 | 模拟降雨土柱入渗试验*** |
CN107421847A (zh) * | 2017-05-17 | 2017-12-01 | 中南大学 | 一种非饱和土层总吸力实时测量装置及其测量方法 |
CN107449706A (zh) * | 2017-06-06 | 2017-12-08 | 湖北工业大学 | 基于分形理论的变形土体饱和、非饱和渗透系数预测方法 |
CN107976393A (zh) * | 2017-11-24 | 2018-05-01 | 金陵科技学院 | 透水混凝土透水系数常水头和变水头测试装置及测试方法 |
CN108760564A (zh) * | 2018-06-29 | 2018-11-06 | 河海大学 | 用于测试非饱和砂土中气体持久性的实验装置 |
CN108760599A (zh) * | 2018-05-16 | 2018-11-06 | 河海大学 | 一种入渗模拟测量装置 |
CN108827495A (zh) * | 2018-08-31 | 2018-11-16 | 瑞安市益华汽车配件有限公司 | 一种汽车废气的温度传感器 |
CN109115972A (zh) * | 2018-08-01 | 2019-01-01 | 陈亚洲 | 一种获取地下水参数一维岩柱模拟水压自记数据实验装置 |
CN109187285A (zh) * | 2018-10-19 | 2019-01-11 | 中国科学院寒区旱区环境与工程研究所 | 土体渗流试验装置 |
CN109325253A (zh) * | 2018-08-01 | 2019-02-12 | 苏州智道势能信息科技有限公司 | 一种密封件密封性模拟测试方法 |
CN109975522A (zh) * | 2019-04-28 | 2019-07-05 | 长沙理工大学 | 一种计算后峰型雨水入渗作用下均质地基基质吸力的方法 |
CN111721688A (zh) * | 2020-06-24 | 2020-09-29 | 湘潭大学 | 一种边界有压渗流试验装置及滑坡模型试验*** |
CN112504337A (zh) * | 2020-11-30 | 2021-03-16 | 同济大学 | 一种基于光栅光纤监测的软土降水与回灌模型试验装置 |
CN113029891A (zh) * | 2019-12-24 | 2021-06-25 | 山西大地华基建材科技有限公司 | 一种反应透水砖透水性能的演示装置 |
CN113435140A (zh) * | 2021-07-22 | 2021-09-24 | 河海大学 | 一种适用于复杂降雨时间序列的降雨入渗边界处理方法 |
CN113435139A (zh) * | 2021-07-22 | 2021-09-24 | 河海大学 | 一种双重渗透介质降雨入渗边界的处理方法 |
CN113466100A (zh) * | 2021-06-21 | 2021-10-01 | 苏州巨力复合材料有限公司 | 一种模拟降雨环境的透水测试方法及装置 |
CN114428033A (zh) * | 2022-03-31 | 2022-05-03 | 深圳市勘察研究院有限公司 | 一种强降雨地表含水量预警监测装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707035A (zh) * | 2012-06-28 | 2012-10-03 | 长安大学 | 一种全自动土柱淋滤实验装置 |
CN102890047A (zh) * | 2012-09-27 | 2013-01-23 | 清华大学 | 一种室内柱状土层降雨入渗自动测定*** |
CN103018421A (zh) * | 2012-11-22 | 2013-04-03 | 中国地质大学(武汉) | 用于模拟石油烃类污染包气带土壤的多节式土柱 |
CN203275257U (zh) * | 2013-05-30 | 2013-11-06 | 长安大学 | 模块化土体渗透性试验仪 |
CN105181533A (zh) * | 2015-08-29 | 2015-12-23 | 西安科技大学 | 基于室内模拟试验的降雪入渗参数测定***及方法 |
CN105181531A (zh) * | 2015-08-29 | 2015-12-23 | 西安科技大学 | 黄土水分迁移规律室内模拟***及特征参数测定方法 |
CN204945001U (zh) * | 2015-09-21 | 2016-01-06 | 河海大学 | 一种室内模拟降雨入渗的试验装置 |
-
2016
- 2016-06-25 CN CN201610480183.0A patent/CN106198341A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707035A (zh) * | 2012-06-28 | 2012-10-03 | 长安大学 | 一种全自动土柱淋滤实验装置 |
CN102890047A (zh) * | 2012-09-27 | 2013-01-23 | 清华大学 | 一种室内柱状土层降雨入渗自动测定*** |
CN103018421A (zh) * | 2012-11-22 | 2013-04-03 | 中国地质大学(武汉) | 用于模拟石油烃类污染包气带土壤的多节式土柱 |
CN203275257U (zh) * | 2013-05-30 | 2013-11-06 | 长安大学 | 模块化土体渗透性试验仪 |
CN105181533A (zh) * | 2015-08-29 | 2015-12-23 | 西安科技大学 | 基于室内模拟试验的降雪入渗参数测定***及方法 |
CN105181531A (zh) * | 2015-08-29 | 2015-12-23 | 西安科技大学 | 黄土水分迁移规律室内模拟***及特征参数测定方法 |
CN204945001U (zh) * | 2015-09-21 | 2016-01-06 | 河海大学 | 一种室内模拟降雨入渗的试验装置 |
Non-Patent Citations (1)
Title |
---|
王红: "非饱和黄土渗透性函数的试验研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106908368B (zh) * | 2017-03-14 | 2020-03-17 | 上海理工大学 | 模拟降雨土柱入渗试验*** |
CN106908368A (zh) * | 2017-03-14 | 2017-06-30 | 上海理工大学 | 模拟降雨土柱入渗试验*** |
CN107421847A (zh) * | 2017-05-17 | 2017-12-01 | 中南大学 | 一种非饱和土层总吸力实时测量装置及其测量方法 |
CN107449706A (zh) * | 2017-06-06 | 2017-12-08 | 湖北工业大学 | 基于分形理论的变形土体饱和、非饱和渗透系数预测方法 |
CN107449706B (zh) * | 2017-06-06 | 2019-11-08 | 湖北工业大学 | 基于分形理论的变形土体饱和、非饱和渗透系数预测方法 |
CN107976393A (zh) * | 2017-11-24 | 2018-05-01 | 金陵科技学院 | 透水混凝土透水系数常水头和变水头测试装置及测试方法 |
CN108760599A (zh) * | 2018-05-16 | 2018-11-06 | 河海大学 | 一种入渗模拟测量装置 |
CN108760564A (zh) * | 2018-06-29 | 2018-11-06 | 河海大学 | 用于测试非饱和砂土中气体持久性的实验装置 |
CN108760564B (zh) * | 2018-06-29 | 2024-04-19 | 河海大学 | 用于测试非饱和砂土中气体持久性的实验装置 |
CN109115972A (zh) * | 2018-08-01 | 2019-01-01 | 陈亚洲 | 一种获取地下水参数一维岩柱模拟水压自记数据实验装置 |
CN109325253A (zh) * | 2018-08-01 | 2019-02-12 | 苏州智道势能信息科技有限公司 | 一种密封件密封性模拟测试方法 |
CN109325253B (zh) * | 2018-08-01 | 2023-05-26 | 苏州智道势能信息科技有限公司 | 一种密封件密封性模拟测试方法 |
CN108827495A (zh) * | 2018-08-31 | 2018-11-16 | 瑞安市益华汽车配件有限公司 | 一种汽车废气的温度传感器 |
CN109187285A (zh) * | 2018-10-19 | 2019-01-11 | 中国科学院寒区旱区环境与工程研究所 | 土体渗流试验装置 |
CN109975522A (zh) * | 2019-04-28 | 2019-07-05 | 长沙理工大学 | 一种计算后峰型雨水入渗作用下均质地基基质吸力的方法 |
CN113029891A (zh) * | 2019-12-24 | 2021-06-25 | 山西大地华基建材科技有限公司 | 一种反应透水砖透水性能的演示装置 |
CN111721688A (zh) * | 2020-06-24 | 2020-09-29 | 湘潭大学 | 一种边界有压渗流试验装置及滑坡模型试验*** |
CN112504337A (zh) * | 2020-11-30 | 2021-03-16 | 同济大学 | 一种基于光栅光纤监测的软土降水与回灌模型试验装置 |
CN112504337B (zh) * | 2020-11-30 | 2022-05-13 | 同济大学 | 一种基于光栅光纤监测的软土降水与回灌模型试验装置 |
CN113466100A (zh) * | 2021-06-21 | 2021-10-01 | 苏州巨力复合材料有限公司 | 一种模拟降雨环境的透水测试方法及装置 |
CN113435140A (zh) * | 2021-07-22 | 2021-09-24 | 河海大学 | 一种适用于复杂降雨时间序列的降雨入渗边界处理方法 |
CN113435139A (zh) * | 2021-07-22 | 2021-09-24 | 河海大学 | 一种双重渗透介质降雨入渗边界的处理方法 |
CN113435139B (zh) * | 2021-07-22 | 2022-03-18 | 河海大学 | 一种双重渗透介质降雨入渗边界的处理方法 |
CN113435140B (zh) * | 2021-07-22 | 2022-03-29 | 河海大学 | 一种适用于复杂降雨时间序列的降雨入渗边界处理方法 |
CN114428033A (zh) * | 2022-03-31 | 2022-05-03 | 深圳市勘察研究院有限公司 | 一种强降雨地表含水量预警监测装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106198341A (zh) | 降水入渗土柱模拟***及非饱和渗透系数测定方法 | |
CN205898792U (zh) | 多状态原状土柱降雨入渗模块化模拟装置 | |
CN105911257A (zh) | 多状态原状土柱入渗模拟***及非饱和渗透系数测定方法 | |
CN106018747A (zh) | 煤矸石淋滤液入渗土柱模拟***及特征参数测定方法 | |
CN105929138A (zh) | 煤矸石淋滤液消融入渗土柱模拟***及特征参数测定方法 | |
CN106093347A (zh) | 多种强度煤矸石淋滤液入渗模拟***及特征参数测定方法 | |
CN105910975A (zh) | 融雪入渗模拟***及非饱和渗透系数测定方法 | |
CN205826656U (zh) | 一种冻土冻融室内试验装置 | |
CN205384202U (zh) | 一种基于水气两相流的降雨入渗规律测试*** | |
CN207717762U (zh) | 一种人工模拟室内降雨蒸发的试验装置 | |
CN110426504A (zh) | 基于示踪水力层析反演的岩溶含水层砂箱试验***及方法 | |
CN109187285A (zh) | 土体渗流试验装置 | |
CN107144523B (zh) | 一种模拟降雨时空变化的边坡试验装置 | |
CN110178595A (zh) | 模拟气候环境变化和降雨时空变化的植被覆盖层模型箱 | |
CN107642352A (zh) | 一种三维模拟石油开发实验装置 | |
CN109064864A (zh) | 一种模拟地热尾水回灌路径的装置及其使用方法 | |
CN106644892A (zh) | 一种可加载的渗透试验装置 | |
CN104459052A (zh) | 一种适用于复杂地表的植被蒸散监测装置及方法 | |
CN102607892A (zh) | 稻田大型原状土柱渗漏水收集测量方法及装置 | |
CN205898790U (zh) | 模块化多种强度煤矸石淋滤液入渗土柱实验装置 | |
CN205982257U (zh) | 基于室内土柱实验的降水入渗模拟装置 | |
CN205982258U (zh) | 基于室内土柱实验的煤矸石淋滤液入渗模拟装置 | |
CN108152476A (zh) | 一种固沟保塬工程措施模拟装置及其模拟方法与应用 | |
CN205898791U (zh) | 煤矸石淋滤液消融入渗室内土柱模拟装置 | |
CN206725404U (zh) | 一种用于连续监测的多功能土柱模拟集成装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161207 |
|
RJ01 | Rejection of invention patent application after publication |