CN105911257A - 多状态原状土柱入渗模拟系统及非饱和渗透系数测定方法 - Google Patents
多状态原状土柱入渗模拟系统及非饱和渗透系数测定方法 Download PDFInfo
- Publication number
- CN105911257A CN105911257A CN201610479162.7A CN201610479162A CN105911257A CN 105911257 A CN105911257 A CN 105911257A CN 201610479162 A CN201610479162 A CN 201610479162A CN 105911257 A CN105911257 A CN 105911257A
- Authority
- CN
- China
- Prior art keywords
- water
- rainfall
- soil sample
- undisturbed soil
- pipe
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 378
- 230000008595 infiltration Effects 0.000 title claims abstract description 76
- 238000001764 infiltration Methods 0.000 title claims abstract description 76
- 238000004088 simulation Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000003204 osmotic effect Effects 0.000 title claims abstract description 6
- 238000012360 testing method Methods 0.000 claims abstract description 65
- 238000002474 experimental method Methods 0.000 claims abstract description 44
- 230000008569 process Effects 0.000 claims abstract description 36
- 230000000694 effects Effects 0.000 claims abstract description 13
- 238000011049 filling Methods 0.000 claims abstract description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 354
- 239000000523 sample Substances 0.000 claims description 253
- 238000002386 leaching Methods 0.000 claims description 73
- 238000012544 monitoring process Methods 0.000 claims description 55
- 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 50
- 239000002985 plastic film Substances 0.000 claims description 36
- 239000004033 plastic Substances 0.000 claims description 33
- 229920003023 plastic Polymers 0.000 claims description 33
- 229920006395 saturated elastomer Polymers 0.000 claims description 32
- 238000012384 transportation and delivery Methods 0.000 claims description 32
- 230000035939 shock Effects 0.000 claims description 30
- 230000008859 change Effects 0.000 claims description 24
- 230000005540 biological transmission Effects 0.000 claims description 18
- 239000004927 clay Substances 0.000 claims description 18
- 238000004891 communication Methods 0.000 claims description 15
- 238000001556 precipitation Methods 0.000 claims description 15
- 229920001971 elastomer Polymers 0.000 claims description 14
- 239000005060 rubber Substances 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 12
- 230000003628 erosive effect Effects 0.000 claims description 9
- 230000033001 locomotion Effects 0.000 claims description 9
- 238000007789 sealing Methods 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000003247 decreasing effect Effects 0.000 claims description 6
- 238000009434 installation Methods 0.000 claims description 6
- 230000003020 moisturizing effect Effects 0.000 claims description 6
- 230000000149 penetrating effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004078 waterproofing Methods 0.000 claims description 6
- 238000005303 weighing Methods 0.000 claims description 6
- 230000001276 controlling effect Effects 0.000 claims description 4
- 230000006735 deficit Effects 0.000 claims description 4
- 229920002430 Fibre-reinforced plastic Polymers 0.000 claims description 3
- 240000002853 Nelumbo nucifera Species 0.000 claims description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 3
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 3
- 229910000831 Steel Inorganic materials 0.000 claims description 3
- 230000033228 biological regulation Effects 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 239000003365 glass fiber Substances 0.000 claims description 3
- 239000004519 grease Substances 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 claims description 3
- 238000005527 soil sampling Methods 0.000 claims description 3
- 239000010959 steel Substances 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 13
- 238000003556 assay Methods 0.000 description 7
- 230000007547 defect Effects 0.000 description 5
- 238000009533 lab test Methods 0.000 description 5
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000004140 cleaning Methods 0.000 description 3
- 210000001215 vagina Anatomy 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 238000004080 punching Methods 0.000 description 2
- 238000012956 testing procedure Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003245 coal Substances 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
- 238000005194 fractionation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 239000011504 laterite Substances 0.000 description 1
- 229910001710 laterite Inorganic materials 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000000203 mixture Substances 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
- 150000003839 salts Chemical class 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
- -1 swelled ground Substances 0.000 description 1
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
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- 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
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/246—Earth materials for water content
-
- 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
- G01N2015/0813—Measuring intrusion, e.g. of mercury
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Environmental & Geological Engineering (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)，弹出件(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-2)、基质吸力探测输出端(7-3)；温湿度探测输出端(7-2)经导线(9)连接温湿度传感器探头(4-5)，基质吸力探测输出端(7-3)经导线(9)连接土壤热传导吸力探头(4-7)；
所述的多种强度降雨入渗模拟装置(C)包括低强度降雨全入渗模拟装置(16)、中强度降雨有径流入渗模拟装置(17)和高强度降雨稳定入渗模拟装置(14)；
所述的低强度降雨全入渗模拟装置(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以上。
基于多状态原状土柱入渗模拟系统的非饱和渗透系数测定方法，包括以下步骤：
步骤一、组装土柱实验标准构件
分别对土柱实验标准构件(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)的湿度信号记录为各个测试点记录时刻对应的含水率θ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小时以上；
步骤七、模拟过程中监测结果分析
对所记录的某一时刻的基质吸力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.本发明的土柱实验标准构件(4)由两个相同的半圆柱体连接成一个圆柱体，可以方便地进行原状样和重塑样的实验，同时方便实验仪器的清洗工作，克服目前土柱实验难以采用原状样进行实验以及后期清洗工作复杂的缺陷。
2.橡胶塞(5-7)与插入件(5)的配合使用，克服现有土柱实验测试元件插入时溢水及水分亏损的缺陷，具有测试结果精确、可靠性高、可重复使用的优点。
3.土柱实验标准构件(4)可根据实验要求多个串联，克服了现有技术中土柱实验中装土部件单节长度大、操作不便的缺陷；且整个系统的各部件采用标准构件，使得实验标准化，实验效率大大提高。
4.将模拟系统与计算机连接，实现自动化控制，具有精确性、标准性高的特点，克服了现阶段大部分的仪器设备人工操作容易出现误差的缺陷。
5.本发明提出了多种强度降雨入渗模拟装置(C)，多种强度降雨入渗模拟装置(C)包括低强度降雨全入渗模拟装置(16)、中强度降雨有径流入渗模拟装置(17)和高强度降雨稳定入渗模拟装置(14)；低强度降雨全入渗模拟装置(16)和高强度降雨稳定入渗模拟装置(14)利用永久磁铁和电线圈磁性相斥的原理，从而控制降雨速度，实现模拟入渗；中强度降雨有径流入渗模拟装置(17)通过空压机实现了降雨强度的控制。
6、基于多状态原状土柱入渗模拟系统实现非饱和渗透系数的测定，并给出了测试步骤及计算方法，具有实用性强，使用效果好，便于推广使用的特点。
7、本发明基于大直径(不小于50cm)的原状土柱进行室内实验，与传统的室内实验小试件相比，能够更好地代表真正土体的实际情况；与原位实测相比，本发明的大直径原状土柱实验也具备周期短、花费少、技术手段成熟的特点。
8、本发明可针对于矿山工程、公路工程、铁路工程和水利水电工程等人类工程活动，尤其是可对比分析人类工程活动扰动作用发生前后不同时期的各种不同类型原状土样(12)的非饱和渗透特性，譬如采煤塌陷扰动作用发生前后不同时期的各种不同类型原状土样(12)的非饱和渗透特性，所述的原状土样(12)包括砾类土、砂类土、细粒土以及诸如淤泥类土、膨胀土、红土、黄土类土、盐渍土、人工填土等特殊土。
9、本发明首次提出多状态原状土柱入渗模拟系统及非饱和渗透系数测定方法，本发明结构新颖合理，能够有效模拟降雨入渗的多种工况，具有准确度高、可靠性强。
附图说明
图1A是基于低强度降雨全入渗模拟装置的系统总图。
图1B是低强度降雨全入渗模拟装置(16)的示意图。
图1C是图1B中B处放大图。
图2A是基于中强度降雨有径流模拟装置的系统总图。
图2B是中强度降雨有径流入渗模拟装置(17)的示意图。
图3A是基于高强度降雨稳定入渗装置的系统总图。
图3B是高强度降雨稳定入渗模拟装置(14)的示意图。
图3C是图3B中B处放大图。
图4是法兰6的示意图。
图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是本发明多状态原状土柱入渗模拟系统的总示意图。
具体实施方式
下面参照附图对本发明做详细叙述。
参照图12，多状态原状土柱入渗模拟系统，包括底部构件(1)、连接在底部构件(1)上的一个或多个串联的土柱实验标准构件(4)以及土柱实验标准构件(4)顶部的多种强度降雨入渗模拟装置(C)；所述的多种强度降雨入渗模拟装置(C)包括低强度降雨全入渗模拟装置(16)、中强度降雨有径流入渗模拟装置(17)和高强度降雨稳定入渗模拟装置(14)；
参照图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)，弹出件(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)的松紧程度，使标准圆柱构件紧密结合；
参照图4，所述的法兰(6)内侧设有法兰螺纹(6-1)，法兰(6)的两端设置有转动把手(6-2)；
所述的计算机(7)的信号端与微控制器(7-0)的信号端相连，微控制器(7-0)设置有温湿度探测输出端(7-2)、基质吸力探测输出端(7-3)；温湿度探测输出端(7-2)经导线(9)连接温湿度传感器探头(4-5)，基质吸力探测输出端(7-3)经导线(9)连接土壤热传导吸力探头(4-7)。
实施例一
本实施例中与本发明配套使用的多种强度降雨入渗模拟装置(C)是低强度降雨全入渗模拟装置(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数值上等于入渗量Qr；此时，查看出渗量测量量杯(3)内渗出原状土样(12)水的量，并将该读数记录为低强度降雨全入渗条件下的出渗量Qc；
步骤六、模拟过程中的各参数监测
多个温湿度传感器探头(4-5)分别对模拟过程中的原状土样(12)的湿度进行监测，并将监测到的信号传给微控制器(7-0)，微控制器(7-0)将监测信号通过串口通信电路实时传输给计算机(7)，计算机(7)接收并记录多个测试点处原状土样(12)的湿度信号，并将各个测试点处原状土样(12)的湿度信号记录为各个测试点记录时刻对应的含水率θ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小时以上；
步骤七、模拟过程中监测结果分析
对所记录的某一时刻的基质吸力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；
步骤六、模拟过程中的各参数监测
具体步骤同实施例一；
步骤七、模拟过程中监测结果分析
具体步骤同实施例一；
步骤八、降雨入渗系数计算
具体步骤同实施例一；
步骤九、非饱和渗透系数计算
具体步骤同实施例一。
实施例三
本实施例中与本发明配套使用的多种强度降雨入渗模拟装置(C)是高强度降雨稳定入渗模拟装置(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以上。
本实施例是基于高强度降雨稳定模拟系统的非饱和渗透系数测定方法包括以下步骤：
步骤一、组装土柱实验标准构件
具体步骤同实施例一；
步骤二、安装原状土样
具体步骤同实施例一；
步骤三、组装土柱实验仪器设备
具体步骤同实施例一；
步骤四、测定原状土样(12)的初始状态
具体步骤同实施例一；
步骤五、模拟高强度降雨稳定入渗边界条件
(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；
步骤六、模拟过程中的各参数监测
具体步骤同实施例一；
步骤七、模拟过程中监测结果分析
具体步骤同实施例一；
步骤八、降雨入渗系数计算
具体步骤同实施例一；
步骤九、非饱和渗透系数计算
具体步骤同实施例一。
Claims (9)
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)在土柱上按照同一列布置；
所述的多种强度降雨入渗模拟装置(C)包括低强度降雨全入渗模拟装置(16)、中强度降雨有径流入渗模拟装置(17)和高强度降雨稳定入渗模拟装置(14)。
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)，弹出件(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-2)、基质吸力探测输出端(7-3)；温湿度探测输出端(7-2)经导线(9)连接温湿度传感器探头(4-5)，基质吸力探测输出端(7-3)经导线(9)连接土壤热传导吸力探头(4-7)。
7.根据权利要求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；
所述的中强度降雨有径流入渗模拟装置(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.基于上述任一权利要求所述的多状态原状土柱入渗模拟系统的非饱和渗透系数测定方法，其特征在于，包括以下步骤：
步骤一、组装土柱实验标准构件
分别对土柱实验标准构件(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)的湿度信号记录为各个测试点记录时刻对应的含水率θ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小时以上；
步骤七、模拟过程中监测结果分析
对所记录的某一时刻的基质吸力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，利用插值法绘制其分布图。
9.根据权利要求8所述的多状态原状土柱入渗模拟系统的非饱和渗透系数测定方法，其特征在于：
步骤五所述的模拟多种强度降雨入渗边界条件，其低强度降雨全入渗的具体模拟步骤为：
(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。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610479162.7A CN105911257A (zh) | 2016-06-25 | 2016-06-25 | 多状态原状土柱入渗模拟系统及非饱和渗透系数测定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610479162.7A CN105911257A (zh) | 2016-06-25 | 2016-06-25 | 多状态原状土柱入渗模拟系统及非饱和渗透系数测定方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105911257A true CN105911257A (zh) | 2016-08-31 |
Family
ID=56759702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610479162.7A Pending CN105911257A (zh) | 2016-06-25 | 2016-06-25 | 多状态原状土柱入渗模拟系统及非饱和渗透系数测定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105911257A (zh) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106324045A (zh) * | 2016-11-01 | 2017-01-11 | 西安理工大学 | 一种黄土的电导率性能及土水特征测试装置 |
| CN106483269A (zh) * | 2016-09-21 | 2017-03-08 | 陕西科技大学 | 一种研究铅/镉在黄土中入渗、迁移及归趋的装置及方法 |
| CN106501158A (zh) * | 2016-12-15 | 2017-03-15 | 中国农业科学院农业资源与农业区划研究所 | 一种室内模拟试验用水分渗漏收集装置 |
| CN107014739A (zh) * | 2017-06-16 | 2017-08-04 | 成都理工大学 | 原状土渗透试验装置及测量渗透系数的方法 |
| CN107255614A (zh) * | 2017-07-27 | 2017-10-17 | 洛阳理工学院 | 重塑非饱和粉质粘土渗透系数的测量及验证方法 |
| CN107957383A (zh) * | 2018-01-17 | 2018-04-24 | 长安大学 | 模拟非饱和带瞬态流-稳态流转化机理的试验装置及方法 |
| CN108444895A (zh) * | 2018-06-14 | 2018-08-24 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
| CN108760599A (zh) * | 2018-05-16 | 2018-11-06 | 河海大学 | 一种入渗模拟测量装置 |
| CN108931623A (zh) * | 2018-04-19 | 2018-12-04 | 中国地质大学（武汉） | 研究降雨入渗下非饱和土的土水特征曲线的模型实验装置 |
| CN109001098A (zh) * | 2018-08-07 | 2018-12-14 | 西南大学 | 一种土壤瞬时渗透特性测试仪及土壤渗透特性计算方法 |
| CN109211756A (zh) * | 2018-10-31 | 2019-01-15 | 湖北工业大学 | 基于vg模型的变形土进气值预测方法 |
| CN109459368A (zh) * | 2018-11-23 | 2019-03-12 | 成都理工大学 | 一种实现多场耦合及原位干湿循环的渗透仪 |
| CN110008302A (zh) * | 2019-04-12 | 2019-07-12 | 长江岩土工程总公司（武汉） | 多影响因子降雨入渗系数定量评价方法 |
| CN110057408A (zh) * | 2019-05-27 | 2019-07-26 | 北京帕克国际工程咨询股份有限公司 | 建筑楼宇主体质量检测装置 |
| CN110296924A (zh) * | 2019-06-28 | 2019-10-01 | 浙江大学 | 一种光热发电储热系统泄漏熔盐在地基中渗流情况的测试系统及方法 |
| CN107607456B (zh) * | 2017-09-11 | 2020-03-24 | 武汉大学 | 基于非线性通量过程控制的非饱和土壤水力传导度的测定方法 |
| CN111238691A (zh) * | 2020-03-09 | 2020-06-05 | 胡增永 | 一种高精度力传感器 |
| CN112432885A (zh) * | 2020-10-26 | 2021-03-02 | 西安理工大学 | 一种测量饱和土体渗流速率的装置和模型选定方法 |
| WO2021128858A1 (zh) * | 2019-12-24 | 2021-07-01 | 中国环境科学研究院 | 浸矿剂迁移测验系统及浸矿剂迁移转化试验系统 |
| CN113435139A (zh) * | 2021-07-22 | 2021-09-24 | 河海大学 | 一种双重渗透介质降雨入渗边界的处理方法 |
| CN116202575A (zh) * | 2023-05-04 | 2023-06-02 | 山东汇杰地理信息科技有限公司 | 一种水土保持土壤流失率监测系统及方法 |
| CN117169089A (zh) * | 2023-11-03 | 2023-12-05 | 东北农业大学 | 一种土壤渗透率检测的张力入渗仪 |
| CN110008302B (zh) * | 2019-04-12 | 2024-06-04 | 长江岩土工程总公司（武汉） | 多影响因子降雨入渗系数定量评价方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080060429A1 (en) * | 2006-08-13 | 2008-03-13 | Noe Stephen A | Runoff Rain Gauge Elevated In Situ Design |
| CN102359084A (zh) * | 2011-07-29 | 2012-02-22 | 中国农业大学 | 河湖包气带渗滤性能的模拟调控系统及其方法 |
| CN102707035A (zh) * | 2012-06-28 | 2012-10-03 | 长安大学 | 一种全自动土柱淋滤实验装置 |
| CN102890047A (zh) * | 2012-09-27 | 2013-01-23 | 清华大学 | 一种室内柱状土层降雨入渗自动测定系统 |
| CN103235018A (zh) * | 2013-04-18 | 2013-08-07 | 中国地质大学（武汉） | 雨水水质监测及入渗模拟综合实验仪 |
| CN203275257U (zh) * | 2013-05-30 | 2013-11-06 | 长安大学 | 模块化土体渗透性试验仪 |
| CN204594833U (zh) * | 2015-05-13 | 2015-08-26 | 深圳市港嘉工程检测有限公司 | 土密度试验灌水法的灌水控制装置 |
| CN204945001U (zh) * | 2015-09-21 | 2016-01-06 | 河海大学 | 一种室内模拟降雨入渗的试验装置 |
-
2016
- 2016-06-25 CN CN201610479162.7A patent/CN105911257A/zh active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080060429A1 (en) * | 2006-08-13 | 2008-03-13 | Noe Stephen A | Runoff Rain Gauge Elevated In Situ Design |
| CN102359084A (zh) * | 2011-07-29 | 2012-02-22 | 中国农业大学 | 河湖包气带渗滤性能的模拟调控系统及其方法 |
| CN102707035A (zh) * | 2012-06-28 | 2012-10-03 | 长安大学 | 一种全自动土柱淋滤实验装置 |
| CN102890047A (zh) * | 2012-09-27 | 2013-01-23 | 清华大学 | 一种室内柱状土层降雨入渗自动测定系统 |
| CN103235018A (zh) * | 2013-04-18 | 2013-08-07 | 中国地质大学（武汉） | 雨水水质监测及入渗模拟综合实验仪 |
| CN203275257U (zh) * | 2013-05-30 | 2013-11-06 | 长安大学 | 模块化土体渗透性试验仪 |
| CN204594833U (zh) * | 2015-05-13 | 2015-08-26 | 深圳市港嘉工程检测有限公司 | 土密度试验灌水法的灌水控制装置 |
| CN204945001U (zh) * | 2015-09-21 | 2016-01-06 | 河海大学 | 一种室内模拟降雨入渗的试验装置 |
Non-Patent Citations (2)
| Title |
|---|
| 李萍: "《黄土中水分迁移规律研究》", 《中国优秀硕士论文全文数据库》 * |
| 王红: "《非饱和黄土渗透性函数的试验研究》", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (32)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106483269A (zh) * | 2016-09-21 | 2017-03-08 | 陕西科技大学 | 一种研究铅/镉在黄土中入渗、迁移及归趋的装置及方法 |
| CN106324045A (zh) * | 2016-11-01 | 2017-01-11 | 西安理工大学 | 一种黄土的电导率性能及土水特征测试装置 |
| CN106501158A (zh) * | 2016-12-15 | 2017-03-15 | 中国农业科学院农业资源与农业区划研究所 | 一种室内模拟试验用水分渗漏收集装置 |
| CN107014739A (zh) * | 2017-06-16 | 2017-08-04 | 成都理工大学 | 原状土渗透试验装置及测量渗透系数的方法 |
| CN107014739B (zh) * | 2017-06-16 | 2023-09-08 | 成都理工大学 | 原状土渗透试验装置及测量渗透系数的方法 |
| CN107255614A (zh) * | 2017-07-27 | 2017-10-17 | 洛阳理工学院 | 重塑非饱和粉质粘土渗透系数的测量及验证方法 |
| CN107607456B (zh) * | 2017-09-11 | 2020-03-24 | 武汉大学 | 基于非线性通量过程控制的非饱和土壤水力传导度的测定方法 |
| CN107957383A (zh) * | 2018-01-17 | 2018-04-24 | 长安大学 | 模拟非饱和带瞬态流-稳态流转化机理的试验装置及方法 |
| CN108931623A (zh) * | 2018-04-19 | 2018-12-04 | 中国地质大学（武汉） | 研究降雨入渗下非饱和土的土水特征曲线的模型实验装置 |
| CN108760599A (zh) * | 2018-05-16 | 2018-11-06 | 河海大学 | 一种入渗模拟测量装置 |
| CN108444895A (zh) * | 2018-06-14 | 2018-08-24 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
| CN108444895B (zh) * | 2018-06-14 | 2020-08-25 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
| CN109001098B (zh) * | 2018-08-07 | 2020-02-21 | 西南大学 | 一种土壤瞬时渗透特性测试仪及土壤渗透特性计算方法 |
| CN109001098A (zh) * | 2018-08-07 | 2018-12-14 | 西南大学 | 一种土壤瞬时渗透特性测试仪及土壤渗透特性计算方法 |
| CN109211756B (zh) * | 2018-10-31 | 2020-11-03 | 湖北工业大学 | 基于vg模型的变形土进气值预测方法 |
| CN109211756A (zh) * | 2018-10-31 | 2019-01-15 | 湖北工业大学 | 基于vg模型的变形土进气值预测方法 |
| CN109459368A (zh) * | 2018-11-23 | 2019-03-12 | 成都理工大学 | 一种实现多场耦合及原位干湿循环的渗透仪 |
| CN109459368B (zh) * | 2018-11-23 | 2023-12-29 | 成都理工大学 | 一种实现多场耦合及原位干湿循环的渗透仪 |
| CN110008302A (zh) * | 2019-04-12 | 2019-07-12 | 长江岩土工程总公司（武汉） | 多影响因子降雨入渗系数定量评价方法 |
| CN110008302B (zh) * | 2019-04-12 | 2024-06-04 | 长江岩土工程总公司（武汉） | 多影响因子降雨入渗系数定量评价方法 |
| CN110057408A (zh) * | 2019-05-27 | 2019-07-26 | 北京帕克国际工程咨询股份有限公司 | 建筑楼宇主体质量检测装置 |
| CN110057408B (zh) * | 2019-05-27 | 2021-10-26 | 北京帕克国际工程咨询股份有限公司 | 建筑楼宇主体质量检测装置 |
| CN110296924A (zh) * | 2019-06-28 | 2019-10-01 | 浙江大学 | 一种光热发电储热系统泄漏熔盐在地基中渗流情况的测试系统及方法 |
| WO2021128858A1 (zh) * | 2019-12-24 | 2021-07-01 | 中国环境科学研究院 | 浸矿剂迁移测验系统及浸矿剂迁移转化试验系统 |
| CN111238691A (zh) * | 2020-03-09 | 2020-06-05 | 胡增永 | 一种高精度力传感器 |
| CN112432885A (zh) * | 2020-10-26 | 2021-03-02 | 西安理工大学 | 一种测量饱和土体渗流速率的装置和模型选定方法 |
| CN113435139A (zh) * | 2021-07-22 | 2021-09-24 | 河海大学 | 一种双重渗透介质降雨入渗边界的处理方法 |
| CN113435139B (zh) * | 2021-07-22 | 2022-03-18 | 河海大学 | 一种双重渗透介质降雨入渗边界的处理方法 |
| CN116202575A (zh) * | 2023-05-04 | 2023-06-02 | 山东汇杰地理信息科技有限公司 | 一种水土保持土壤流失率监测系统及方法 |
| CN116202575B (zh) * | 2023-05-04 | 2023-07-28 | 山东汇杰地理信息科技有限公司 | 一种水土保持土壤流失率监测系统及方法 |
| CN117169089A (zh) * | 2023-11-03 | 2023-12-05 | 东北农业大学 | 一种土壤渗透率检测的张力入渗仪 |
| CN117169089B (zh) * | 2023-11-03 | 2024-01-26 | 东北农业大学 | 一种土壤渗透率检测的张力入渗仪 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105911257A (zh) | 多状态原状土柱入渗模拟系统及非饱和渗透系数测定方法 | |
| CN205898792U (zh) | 多状态原状土柱降雨入渗模块化模拟装置 | |
| CN106198341A (zh) | 降水入渗土柱模拟系统及非饱和渗透系数测定方法 | |
| CN106018747A (zh) | 煤矸石淋滤液入渗土柱模拟系统及特征参数测定方法 | |
| CN106093347A (zh) | 多种强度煤矸石淋滤液入渗模拟系统及特征参数测定方法 | |
| Dong et al. | Aeolian sand transport: a wind tunnel model | |
| Walter et al. | Shear-stress partitioning in live plant canopies and modifications to Raupach’s model | |
| CN105929138A (zh) | 煤矸石淋滤液消融入渗土柱模拟系统及特征参数测定方法 | |
| CN105910975A (zh) | 融雪入渗模拟系统及非饱和渗透系数测定方法 | |
| CN109752303A (zh) | 一种非连续非饱和土壤降雨入渗物理模拟系统及测定方法 | |
| CN107665270B (zh) | 基于Lu模型的大坝下游河岸带热流耦合模拟构建方法 | |
| CN107655653A (zh) | 一种研究积雪漂移的风洞 | |
| CN109064864A (zh) | 一种模拟地热尾水回灌路径的装置及其使用方法 | |
| Western et al. | Soil moisture and runoff processes at Tarrawarra | |
| CN109324170A (zh) | 多场景土体渗流和水盐迁移规律的检测方法 | |
| CN109001098A (zh) | 一种土壤瞬时渗透特性测试仪及土壤渗透特性计算方法 | |
| CN103543095A (zh) | 一种采矿废石堆体降雨入渗实验装置及方法 | |
| CN106353359A (zh) | 一种测试土体试样热物性参数的装置 | |
| CN205898790U (zh) | 模块化多种强度煤矸石淋滤液入渗土柱实验装置 | |
| Zhang et al. | Impact of preferential and lateral flows of water on single‐ring measured infiltration process and its analysis | |
| CN108152476A (zh) | 一种固沟保塬工程措施模拟装置及其模拟方法与应用 | |
| CN205982258U (zh) | 基于室内土柱实验的煤矸石淋滤液入渗模拟装置 | |
| CN205982257U (zh) | 基于室内土柱实验的降水入渗模拟装置 | |
| Wahba et al. | Evaluation of DRAINMOD‐S for simulating water table management under semi‐arid conditions | |
| CN106442603A (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: 20160831 |
|
| RJ01 | Rejection of invention patent application after publication |