技术领域
[0001] 本
发明涉及土工测试测量技术领域，具体为一种基于土壤水盐变化确定冻土上限变化规律方法及装置。
背景技术
[0002] 据统计，全球陆地面积70％存在冻土，其中14％为永久冻土，56％为季节冻土；我国是世界上第三冻土大国，其中
多年冻土分布面积为2.068×106km2，占我国国土面积的21.5％，季节性冻土分布面积很广，为5.137×106km2，占我国国土面积的53.5％，两者合计约占全国总面积的3/4左右，冻土路基的上限变化主要是气温变化致使土壤反复冻融而引起的，而土壤的冻融是一个非常复杂的过程，它伴随物理、化学及
力学现象和子过程，最主要包括上限
位置变化、水分迁移和盐分积聚等。
[0003] 因此可以了解到影响多年冻土地区路
稳定性的主要部位是上限附近及其上部的季节融化层。而季节融化层比较活跃，随着
气候和地质条件的改变所发生的融冻变化直接危及路基稳定和道路建筑的安全，而上限的位置及其变化受各种自然条件综合作用，直接反映季节融化层的特性，因此确定多年冻土上限的位置及其变化是冻土地区工程勘察的一个重要内容。
发明内容
[0004] (一)解决的技术问题
[0005] 针对
现有技术的不足，本发明提供了一种基于土壤水盐变化确定冻土上限变化规律方法及装置，解决了确定多年冻土上限的位置及其变化规律是冻土地区工程勘察困难的问题。
[0006] (二)技术方案
[0007] 为实现以上目的，本发明通过以下技术方案予以实现：一种基于土壤水盐变化确定冻土上限变化规律方法，经过冻土制样取样及在试验箱进行冻融试验后，公式①⑤的上边界条件可以分别由季节冻土试样最上表面插入式
传感器固定环中的插入式土壤水分与盐分传感器测量的土壤含水率θ、土壤含盐量w数据求得的函数确定，传感器探针长短不一，数据取平均值后再计算函数，公式①⑤的下边界条件可以分别由多年冻土试样最下表面，即试样槽可移动
滑板与多年冻土试样
接触面处的插入式传感器固定环中的插入式土壤水分与盐分传感器测得的土壤含水率θ、土壤含盐量w数据求得的函数确定，公式①⑤的初始条件由季节冻土和多年冻土试样中由下至上多个插入式传感器固定环中的插入式土壤水分与盐分传感器测得土壤含水率θ、土壤含盐量w数据求得的函数确定；
[0008] 公式①⑤的初始条件及边界条件确定后再根据分离变量法求解解析解得式②⑥，最后结合土壤水分、盐分判断指标由式②⑥计算判断冻土上限变化过程如下：
[0009] (1)由土壤冻融水分传导方程
[0010] 式中：θ为土壤体积含水率(cm3·cm-3)，D为土中水的扩散系数(cm3·s-1)，t为时间(s)，z为土体z轴方向的尺寸分量(cm)，l为季节冻土试样高度(cm)，l1为多年冻土试样高度(cm)；
[0011] 初始条件为t＝0时，θ(z,0)＝θ1(z)；
[0012] 边界条件为z＝0时，θ(0,t)＝θ2(t)；
[0013] z＝l+l1时，θ(l+l1,t)＝θ3(t)；
[0014] 由分离变量法得解析解为：
[0015]
[0016] 其中
[0017]
[0018] t0为测试试验结束时间，n＝1,2,3···；θi为插入式土壤水分与盐分传感器测试数据确定的含水率函数；
[0019] 根据公式②能得到任意深度z任意时间t的含水率θ，以多年冻土与季节冻土的界面上下含水率差值较大为判断依据，由计算机将时间t，t可根据具体试验要求规定调节，代入公式②进行赋值试算，直至某两个相近z处含水率差值Δθ出现突变，Δθ明显激增，即确定为多年冻土上限位置z1，z1为两个相近z值得平均值，若出现多个z1，以最大值为准；
[0020] 不断调节时间t可以在试验时间t1模拟得到现场时间t1'内多年冻土上限位置变化规律；
[0021] (2)由土壤冻融盐分传导方程
[0022] 式中：w为土壤含盐率(g·cm-3)，Dsh为浓度梯度下盐分的水动力弥散系数(cm2·s-1)；
[0023] 初始条件为t＝0时，w(z,0)＝w1(z)；
[0024] 边界条件为z＝0时，w(0,t)＝w2(t)；
[0025] z＝l+l1时，w(l+l1,t)＝w3(t)；
[0026] 由分离变量法得解析解为：
[0027]
[0028] 其中
[0029]
[0030] t0为测试试验结束时间，n＝1,2,3···；wi为插入式土壤水分与盐分传感器测试数据确定的含盐率函数；
[0031] 根据公式⑥能得到任意深度z任意时间t的含盐率w，以多年冻土与季节冻土的界面上下含盐率差值较大为判断依据，由计算机将时间t，代入公式⑥进行赋值试算，直至某两个相近两个z处含盐率差值Δw发生突变，Δw明显激增，即确定为多年冻土上限位置z2，z2为两个相近z值得平均值，若试算出现多个z2，以最大值为准；
[0032] 不断调节时间t可以在试验时间t1模拟得到现场时间t1'内多年冻土上限位置变化规律；
[0033] (3)综上，根据不同公式解析解求得的多年冻土上限位置为
[0034] 根据相似理论，试验所观测到的试样模型冻融变化量与现场
原型冻土上限冻融变化量之比也为1:a，故现场原型冻土上限位置为az，根据此公式可以求得现场冻土上限位置变化规律。
[0035] 一种基于土壤水盐变化确定冻土上限变化规律装置，包括外
箱体、内箱体、
数据采集模
块和系统控制及数据采集计算机，所述内箱体位于外箱体内部，所述外箱体和内箱体的上端分别设置有外箱盖和内箱盖，所述内箱体的左侧外壁和外箱体的内壁之间的位置安装有加热模组和制冷模组，所述加热模组和制冷模组分别与内箱体连通，所述内箱体的右侧外壁和外箱体的内壁之间的位置设置有补水装置，所述内箱体左侧内壁上方安装有日光模拟灯，所述内箱体内壁底部固定连接有金属试样台底座，所述金属试样台底座内部安装有转动
电机，所述转动电机
输出轴安装有金属
转轴，所述金属转轴顶部固定连接有筒状磁
铁旋转台，所述筒状
磁铁旋转台顶部设置有试样槽，所述试样槽的顶端安装有红外测距传感器，所述内箱盖下方安装有
压实装置,所述试样槽设置有若干个，所述试样槽包括试样槽外壁、试样槽内壁和试样槽底座，所述试样槽底座底部安装有升降
气缸，所述升降气缸外侧设置有升降气缸导气管孔和升降气缸
开关，所述升降气缸顶部固定连接有可移动滑板，所述试样槽内壁固定连接有保温层，所述可移动滑板两侧与试样槽的保温层内壁相接触,所述可移动滑板上放置有多年冻土试样，所述多年冻土试样的底部以及四周的位置均放置有
半导体制冷片，所述多年冻土试样顶部放置有季节冻土试样，所述多年冻土试样与季节冻土试样内部设置有插入式传感器固定环，所述插入式传感器固定环设置有多个，且由下至上均匀设置，所述插入式传感器固定环包括固定环内壁和固定环外壁，所述固定环内壁表面设置有传感器固定卡槽，所述传感器固定卡槽均匀设置，所述传感器固定卡槽内部安装有插入式土壤水分与盐分传感器，所述插入式土壤水分与盐分传感器内侧设置有土壤水分和盐分传感器感应探针，所述土壤水分和盐分传感器感应探针位于插入式传感器固定环内侧,且土壤水分和盐分传感器感应探针设置为不同长度。
[0036] 优选的，所述t可根据具体试验要求规定调节，所述t1远远小于t1'。
[0037] 优选的，所述压实装置包括驱动气缸，所述驱动气缸顶部通过加压装置固定
螺栓与内箱盖固定连接，所述驱动气缸的伸缩
杆底部安装有
支撑架，所述支撑架中部通过支撑架连接螺栓与伸缩杆底部连接，所述支撑架底部固定连接有压头，所述压头的数量和位置均与试样槽的数量和位置一一对应。
[0038] 优选的，所述外箱盖内箱盖上分别设有
拉环，所述外箱体、内箱体、外箱盖和内箱盖均使用
隔热金属材料制成，且内壁均涂抹保温材料层,所述日光模拟灯控制端与系统控制及数据采集计算机连接,所述试样槽内壁与保温层之间的位置设置有隔
热层。
[0039] 优选的，所述插入式土壤水分与盐分传感器和红外测距传感器均通过试验数据传输
导线与数据采集模块连接，所述数据采集模块输出端与系统控制及数据采集计算机输入端连接，所述加热模组和制冷模组均通过试验数据传输导线与系统控制及数据采集计算机连接。
[0040] 优选的，所述补水装置包括
马氏瓶、蒸馏水传输
导管、补水装置喷头、蒸馏水传输导
管接头、可拆卸蒸馏水传输导管和筒状旋转台积水回收瓶，所述马氏瓶位于内箱体的右侧外壁和外箱体的内壁之间的位置，所述马氏瓶通过蒸馏水传输导管、蒸馏水传输导管接头和可拆卸蒸馏水传输导管与补水装置喷头连接，所述补水装置喷头位于内箱体内部筒状磁铁旋转台的顶部，所述补水装置喷头的控制端与系统控制及数据采集计算机电性连接。
[0041] (三)有益效果
[0042] 本发明提供了一种基于土壤水盐变化确定冻土上限变化规律方法及装置。具备以下有益效果：
[0043] (1)、本发明在试验模拟时，不仅可以监测季节冻土在冻融过程中的水分变化、盐分迁移变化状况，还可以根据传感器测量确定土壤水分传导、盐分传导微分方程的初始条件及边界条件，并根据分离变量法求解出解析解，依据以土壤含水率、土壤含盐率作为指标的判别标准，计算冻土上限变化规律，试验
进程及试验数据监测由
计算机系统控制，试验结果更加准确。
[0044] (2)、该基于土壤水盐变化确定冻土上限变化规律方法及装置，现场取样不仅适用于原状土，还能为了研究多因素对冻土上限的影响，室内根据拟研究的土质类型自行调配制备试样。
[0045] (3)、该基于土壤水盐变化确定冻土上限变化规律方法及装置，试验箱可以较准确模拟冻土现场实际的
温度、降水、光照变化情况，变化进程由计算机系统控制，可以在最大程度上模拟自然环境变化。筒状磁铁旋转台带动试样槽旋转，使模拟的温度、降水、日照等环境变化更加均匀。
[0046] (4)、该基于土壤水盐变化确定冻土上限变化规律方法及装置，试样槽可以通过调节
底板高度来达到试验要求的高度，既避免了因为试样槽太高，致使土体试样距离试样槽上表面距离太大，而出现土体温度降水光照模拟误差，还可以避免出现因试样尺寸需求而要准备多个不同规格试样槽的情况。
[0047] (5)、该冻土上限变化规律测试方法，利用相似准则将大型冻土试验原型模型化和简易化，模型试验结果可用于反推原型。根据建立的比例关系，调整模型尺寸和试验时间，大大
加速了试验进程，能在短时间内模拟现场长时间的环境变化，大大缩短了监测时间。
附图说明
[0048] 图1为本发明整体的结构示意图；
[0049] 图2为本发明试样槽的结构示意图；
[0050] 图3为本发明插入式传感器固定环示意图；
[0051] 图4为本发明筒状磁铁旋转台及试样槽放置的结构示意图；
[0052] 图5为本发明可拆卸加压装置的结构示意图；
[0053] 图6为本发明试样槽尺寸的结构示意图。
[0054] 图中：1外箱体、2内箱体、3外箱盖、4内箱盖、5金属试样台底座、6金属转轴、7筒状磁铁旋转台、8试样槽外壁、9试样槽内壁、10试样槽底座、11保温层、12隔热层、13升降气缸、14升降气缸导气管孔、15升降气缸开关、16可移动滑板、17多年冻土试样、18季节冻土试样、
19半导体制冷片、20插入式土壤水分与盐分传感器、21插入式传感器固定环、22固定环内壁、23传感器固定卡槽、24土壤水分和盐分传感器感应探针、25红外测距传感器、26制冷模组、27加热模组、28补水装置喷头、29蒸馏水传输导管、30马氏瓶、31日光模拟灯、32试验数据传输导线、33数据采集模块、34系统控制及数据采集计算机、35加压装置固定螺栓、36驱动气缸、37伸缩杆、38支撑架连接螺栓、39支撑架、40压头、41蒸馏水传输导管接头、42可拆卸蒸馏水传输导管、43筒状旋转台积水回收瓶。
具体实施方式
[0055] 下面将结合本发明
实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0056] 请参阅图1-6，本发明提供一种技术方案：一种基于土壤水盐变化确定冻土上限变化规律装置，包括外箱体1、内箱体2、数据采集模块33和系统控制及数据采集计算机34，内箱体2位于外箱体1内部，外箱体1和内箱体2的上端分别设置有外箱盖3和内箱盖4，内箱体2的左侧外壁和外箱体1的内壁之间的位置安装有加热模组27和制冷模组26，加热模组27和制冷模组26分别与内箱体2连通，内箱体2的右侧外壁和外箱体1的内壁之间的位置设置有补水装置，内箱体2左侧内壁上方安装有日光模拟灯31，内箱体2内壁底部固定连接有金属试样台底座5，金属试样台底座5内部安装有转动电机，转动电机输出轴安装有金属转轴6，金属转轴6顶部固定连接有筒状磁铁旋转台7，筒状磁铁旋转台7顶部设置有试样槽，试样槽的顶端安装有红外测距传感器25，内箱盖4下方安装有压实装置,试样槽设置有若干个，试样槽包括试样槽外壁8、试样槽内壁9和试样槽底座10，试样槽底座10底部安装有升降气缸13，升降气缸13外侧设置有升降气缸导气管孔14和升降气缸开关15，升降气缸13顶部固定连接有可移动滑板16，试样槽内壁9固定连接有保温层11，可移动滑板16两侧与试样槽的保温层11内壁相接触,可移动滑板16上放置有多年冻土试样17，多年冻土试样17的底部以及四周的位置均放置有半导体制冷片19，多年冻土试样17顶部放置有季节冻土试样18，多年冻土试样17与季节冻土试样18内部设置有插入式传感器固定环21，插入式传感器固定环21设置有多个，且由下至上均匀设置，插入式传感器固定环21包括固定环内壁22和固定环外壁，固定环内壁22表面设置有传感器固定卡槽23，传感器固定卡槽23均匀设置，传感器固定卡槽23内部安装有插入式土壤水分与盐分传感器20，插入式土壤水分与盐分传感器20内侧设置有土壤水分和盐分传感器感应探针24，土壤水分和盐分传感器感应探针24位于插入式传感器固定环21内侧,且土壤水分和盐分传感器感应探针24设置为不同长度。
[0057] 压实装置包括驱动气缸36，驱动气缸36顶部通过加压装置固定螺栓35与内箱盖4固定连接，驱动气缸36的伸缩杆37底部安装有支撑架39，支撑架39中部通过支撑架连接螺栓38与伸缩杆37底部连接，支撑架39底部固定连接有压头40，压头40的数量和位置均与试样槽的数量和位置一一对应。
[0058] 外箱盖3内箱盖4上分别设有拉环，外箱体1、内箱体2、外箱盖3和内箱盖4均使用隔热金属材料制成，且内壁均涂抹保温材料层,日光模拟灯31控制端与系统控制及数据采集计算机34连接,试样槽内壁9与保温层11之间的位置设置有隔热层12。
[0059] 插入式土壤水分与盐分传感器20和红外测距传感器25均通过试验数据传输导线32与数据采集模块33连接，数据采集模块33输出端与系统控制及数据采集计算机34输入端连接，加热模组27和制冷模组26均通过试验数据传输导线32与系统控制及数据采集计算机
34连接。
[0060] 补水装置包括马氏瓶30、蒸馏水传输导管29、补水装置喷头28、蒸馏水传输导管接头41、可拆卸蒸馏水传输导管42和筒状旋转台积水回收瓶43，马氏瓶30位于内箱体2的右侧外壁和外箱体1的内壁之间的位置，马氏瓶30通过蒸馏水传输导管29、蒸馏水传输导管接头41和可拆卸蒸馏水传输导管42与补水装置喷头28连接，补水装置喷头28位于内箱体2内部筒状磁铁旋转台7的顶部，补水装置喷头28的控制端与系统控制及数据采集计算机34电性连接。
[0061] 使用时，在可移动滑板16上方放置多年冻土试样17，多年冻土试样17上方是季节冻土试样18，多年冻土试样17的底部及内部四周放置半导体制冷片19用于多年冻土制样及保持土样的低温条件，季节冻土试样18与多年冻土试样17中由下至上排列有多个插入式传感器固定环21，每个插入式传感器固定环21有多个传感器固定卡槽23，插入式土壤水分与盐分传感器20固定于传感器固定卡槽23内，土壤水分和盐分传感器感应探针24测得的
信号通过试验数据传输导线32传输至外部的数据采集模块33和系统控制及数据采集计算机34中，制冷模组26和加热模组27加热与降温速率由计算机PID控制
软件控制，PID温控软件安装在系统控制及数据采集计算机34中，可以同步显示当前温度，并自动绘制温度曲线，软件还可以自动设定温度变化，如正弦、余弦、线性等，温控系统适用于-40℃—+80℃，数字显示设置温度和现实温度，主要控制制冷模组26和加热模组27工作；
[0062] 系统控制及数据采集计算机34内置湿度控制系统，补水装置采用蒸馏水喷淋方式，根据试验需要可以安装和拆卸补水装置，由储存蒸馏水的马氏瓶30与补水装置喷头28通过蒸馏水传输导管29、蒸馏水传输导管接头41和可拆卸蒸馏水传输导管42连接，蒸馏水喷淋速率及时间由计算机控制软件控制，补水装置的作用是在试验环境中土样与马氏瓶30中的蒸馏水保持水力联系，以模拟土体在冻融过程中的环境水分补给条件，内箱体2内部安装的日光模拟灯31主要模拟土体在冻融过程中的太阳光照情况，使试验环境更加接近自然环境，光照强度及时间由计算机控制软件控制；
[0063] 插入式土壤水分与盐分传感器20固定在传感器固定卡槽23中形成测量网络埋于土样，插入式传感器固定环21由下至上均匀排布在季节冻土试样18与多年冻土试样17中，插入式土壤水分与盐分传感器20不仅可以测量确定水分传导方程与盐分传导方程的初始条件，季节冻土最上表面的插入式土壤水分与盐分传感器20可以测量确定水分传导方程与盐分传导方程的上边界条件，多年冻土最下表面的插入式土壤水分与盐分传感器20可以测量确定水分传导方程与盐分传导方程的下边界条件；
[0064] 加压装置可以通过加压装置固定螺栓35安装在内箱盖4上，可以随内箱盖4打开安装与拆卸，试样制作完成后拆卸下来；压头40可以竖直压入试样槽中，压头40的压实速率和进程可以由系统控制及数据采集计算机34内置的控制系统进行控制，系统控制及数据采集计算机34可以控制实验进程，系统控制及数据采集计算机34包含试验控制界面、土壤水分数据采集界面、土壤盐分数据采集界面；
[0065] 在试验中需要采用模拟结构及方法需要如下理论进行模拟试验时间加速：
[0066] 由于现场冻土试验耗时较长，本发明为减少时间成本，达到在较短时间内模拟长时间环境变化的效果，基于相似原理，建立冻土模型尺寸与试验时间比例关系达到加速试验时间的目的。本发明依托三维水分传输微分方程、盐分传输微分方程进行冻土模型试验相似准则推导过程如下：
[0067] (1)水分传输微分控制方程：
[0068]
[0069] 0≤x≤d，0≤y≤d，0≤z≤l+l1，t≥0
[0070] 式中：θ为土壤体积含水率(cm3·cm-3)；D为土中水的扩散系数(cm3·s-1)；t为时间(s)；x，y，z分别为土体x轴y轴z轴方向的尺寸分量(cm)，l为季节冻土试样高度(cm)，l1为多年冻土试样高度(cm)，d为试样槽内径(cm)。
[0071] 初始条件为t＝0时，θ(x,y,z,0)＝θ4(x,y,z)；
[0072] 边界条件为z＝0时，θ(x,y,0,t)＝θ5(x,y,t)；z＝l+l1时，θ(x,y,l+l1,t)＝θ6(x,y,t)
[0073] x＝0，θ(0,y,z,t)＝θ7(y,z,t)；x＝d，θ(d,y,z,t)＝θ8(y,z,t)；
[0074] y＝0，θ(x,0,z,t)＝θ9(x,z,t)；y＝d，θ(x,d,z,t)＝θ10(x,z,t)；
[0075] 根据相似理论可以化简得：
[0076] 式中CD、Ct、Cl分别为土体的土中水的扩散系数、时间、几何的相似常数；θi为插入式土壤水分与盐分传感器20测试数据确定的含水率函数。
[0077] (2)盐分传输微分控制方程：
[0078]
[0079] 0≤x≤d，0≤y≤d，0≤z≤l+l1，t≥0
[0080] 式中：w为土壤含盐率(g·cm-3)；Dsh为浓度梯度下盐分的水动力弥散系数(cm2·s-1)。
[0081] 初始条件为t＝0时，w(x,y,z,0)＝w4(x,y,z)；
[0082] 边界条件为z＝0时，w(x,y,0,t)＝w5(x,y,t)，z＝l+l1时，w(x,y,l+l1,t)＝w6(x,y,t)；
[0083] x＝0，w(0,y,z,t)＝w7(y,z,t)；x＝d，w(d,y,z,t)＝w8(y,z,t)；
[0084] y＝0，w(x,0,z,t)＝w9(x,z,t)；y＝d，w(x,d,z,t)＝w10(x,z,t)；
[0085] 根据相似理论可以化简得：
[0086] 式中CDsh、Cθ、Ct、Cl分别为土体的浓度梯度下盐分的水动力弥散系数、体积
含水量、时间、几何的相似常数；wi为插入式土壤水分与盐分传感器20测试数据确定的含盐率函数。
[0087] 公式⑨⑩的初始条件及边界条件通过安置在试样上表面、下表面及整体高度的插入式土壤水分与盐分传感器20测量数据后求解函数得到。
[0088] 综上，当采用原状土为制作模型时，可得：
[0089] CD＝CDsh＝1
[0090] 令：Cθ＝Cw＝1(即试验湿度、盐分浓度均采用与天然条件同样值，相似常数均为1)[0091] 当原状土作为试样模
型材料时，试验水分、盐分与现场原型的相似比例为Cθ＝Cw＝1时，CD＝CDsh＝1，可得时间的比例系数为模型几何尺寸比例系数的平方，故假设试验模型与现场原型几何比例为Cl＝1:a(a为常数)时，则试验时间与现场比例为Ct＝1:a2，且试验所观测到的试样模型冻融变化量与现场原型冻土上限冻融变化量之比也为1:a。
[0092] 一种基于土壤水盐变化确定冻土上限变化规律方法的测试步骤如下：
[0093] (1)试验制样：根据上述相似理论进行试样尺寸计算：由具体试验要求确定试验时间t1及要模拟的现场时间t1'(t1远远小于t1')求得模型时间比例为Ct＝t1:t1'＝1:a2(a为常数)；所以可得模型几何比例为Cl＝1:a，再根据现场要模拟的季节冻土高度L与多年冻土高度L1求得模型的季节冻土试样高度为 多年冻土试样高度为
[0094] 试验制样可分两种情况：第一种原状土制样，将取自模拟现场的季节冻土与多年冻土按照试验计算尺寸进行切割制成试样；多年冻土试样取自原状冻土多年冻土层，取土深度在-5～-10m范围内；季节冻土试样取自原状冻土季节冻融循环层，取土深度在-1～-5m范围内；
[0095] 第二种室内制备土样，为了研究多因素对冻土上限的影响，可室内制样；将取自现场的土样进行烘干洗盐，再根据模拟的土壤含水量、含盐率数据进行室内重新制样。在冻土制备时，首先根据计算多年冻土试样17高度l1和季节冻土试样18高度l选择相应的模具，调节试样槽升降气缸13以改变试样槽高度，试样槽高度调节由红外测距传感器25测量确定。根据计算出的多年冻土与季节冻土两种试样高度l1与l、试样槽内径d及土样干
密度ρ再计
2
算出所需的经多次烘干洗盐后的季节冻土试样18干
质量m1＝ρπ(d/2) l和多年冻土试样17干质量m2＝ρπ(d/2)2l1。再根据已知数据的冻土体积含水率θ、含盐率w准备相应质量的去离子水与NaCl、Na2SO4、NaHCO3，将NaCl、Na2SO4、NaHCO3按一定比例混合均匀后溶解在相应质量的去离子水中，待盐分完全溶解至溶液中后，将盐水溶液加入干土中，边加入溶液边搅拌，使得盐水溶液与土混合均匀，密封至塑料桶中保持12小时，目的是使土体内水分盐分分布均匀。在模型箱
侧壁上涂一层防冻
硅油，以减小侧壁的
摩擦力；先分层装多年冻土，装入时分层埋入插入式传感器固定环21，土样底部以及四周的位置均放置半导体制冷片19，用压实装置对每层填土进行夯实，直至达到规定高度完成，然后将半导体制冷片19调节到最低温度制冷12小时，上层刮平，多年冻土试样17制作完成后再分层填入季节冻土，装入时也分层埋入插入式传感器固定环21，同样用压实装置对每层填土进行夯实直至达到规定高度完成。压实装置使用完毕拆卸下来。
[0096] (2)试验进程控制：本试验的温度、湿度、光照条件应与现场原状土保持一致，取样或制样结束后，将试样槽安放在试验内箱体2的筒状磁铁旋转台7，再启动电机带动金属转轴6和筒状磁铁旋转台7转动，从而带动试样槽在箱体内部均匀转动。
[0097] 根据试验时间t1及要模拟的现场时间t1'(t1远远小于t1')，冻土地区现场t1'时间的温度、降水、光照现场数据，将现场t1'时间的温度变化数据输入系统控制及数据采集计算机34中，计算机系统控制制冷模组26和加热模组27进行温度调控，将这些数据在t1时间内重现于试验箱。同样根据调查的降水数据，输入计算机后用计算机软件控制补水装置喷头28的喷淋速率，使现场t1'时间的降水量在t1时间均匀降于筒状磁铁旋转台7。同理调节日光模拟灯31的照射
频率与时间，使现场t1'时间的光照
辐射量在t1时间均匀照射于试验箱。
[0098] (3)试验结果测量：所述的土壤水分、盐分监测系统，为插入式土壤水分与盐分传感器20，插入式土壤水分与盐分传感器20安置在插入式传感器固定环21中并埋设在土样中，试样槽季节冻土试样18自上表面至多年冻土试样17与可移动滑板16接触面的高度l+l1中均匀分布多个插入式传感器固定环21，传感器测得数据会由数据采集模块采集输入计算机系统保存输出。
[0099] 冻土上限计算过程如下：
[0100] 经过冻土制样取样及在试验箱进行冻融试验后，公式①⑤的上边界条件可以分别由季节冻土试样18最上表面插入式传感器固定环21中的插入式土壤水分与盐分传感器20测量的土壤含水率θ、土壤含盐量w数据求得的函数确定，传感器探针长短不一，数据取平均值后再计算函数，公式①⑤的下边界条件可以分别由多年冻土试样17最下表面，即试样槽可移动滑板16与多年冻土试样17接触面处的插入式传感器固定环21中的插入式土壤水分与盐分传感器20测得的土壤含水率θ、土壤含盐量w数据求得的函数确定，公式①⑤的初始条件由季节冻土试样18和多年冻土试样17中由下至上多个插入式传感器固定环21中的插入式土壤水分与盐分传感器20测得土壤含水率θ、土壤含盐量w数据求得的函数确定；
[0101] 公式①⑤的初始条件及边界条件确定后再根据分离变量法求解解析解得式②⑥，最后结合土壤水分、盐分判断指标由式②⑥计算判断冻土上限变化过程如下：
[0102] (1)由土壤冻融水分传导方程
[0103] 式中：θ为土壤体积含水率(cm3·cm-3)，D为土中水的扩散系数(cm3·s-1)，t为时间(s)，z为土体z轴方向的尺寸分量(cm)，l为季节冻土试样高度(cm)，l1为多年冻土试样高度(cm)；
[0104] 初始条件为t＝0时，θ(z,0)＝θ1(z)；
[0105] 边界条件为z＝0时，θ(0,t)＝θ2(t)；
[0106] z＝l+l1时，θ(l+l1,t)＝θ3(t)；
[0107] 由分离变量法得解析解为：
[0108]
[0109] 其中
[0110]
[0111] t0为测试试验结束时间，n＝1,2,3···；θi为插入式土壤水分与盐分传感器20测试数据确定的含水率函数；
[0112] 根据公式②能得到任意深度z任意时间t的含水率θ，由于多年冻土、季节冻土的含水率具有一定差值，多年冻土的水为冻结状态，季节冻土的液态水含量远远高于多年冻土，故以多年冻土与季节冻土的界面上下含水率差值较大为判断依据，由计算机将时间t，t可根据具体试验要求规定调节，代入公式②进行赋值试算，直至相近z处含水率差值Δθ出现突变，Δθ明显激增，即确定多年冻土上限位置z1，z1为两个相近z值得平均值，若出现多个z1，以最大值为准；
[0113] 不断调节时间t可以在试验时间t1模拟得到现场时间t1'内多年冻土上限位置变化规律；
[0114] (2)由土壤冻融盐分传导方程
[0115] 式中：w为土壤含盐率(g·cm-3)，Dsh为浓度梯度下盐分的水动力弥散系数(cm2·s-1)；
[0116] 初始条件为t＝0时，w(z,0)＝w1(z)；
[0117] 边界条件为z＝0时，w(0,t)＝w2(t)；
[0118] z＝l+l1时，w(l+l1,t)＝w3(t)；
[0119] 由分离变量法得解析解为：
[0120]
[0121] 其中
[0122]
[0123] t0为测试试验结束时间，n＝1,2,3···；wi为插入式土壤水分与盐分传感器20测试数据确定的含盐率函数；
[0124] 根据公式⑥能得到任意深度z任意时间t的含盐率w，由于多年冻土、季节冻土的含盐率具有一定差值，土壤盐分以离子形式存在于液态水中，多年冻土的水为冻结状态，季节冻土的液态水含量远远高于多年冻土，因此季节冻土的含盐率也高于多年冻土，故以多年冻土与季节冻土的界面上下含盐率差值较大为判断依据，由计算机将时间t，代入公式⑥进行赋值试算，直至某相近两个z处含盐率差值Δw发生突变，Δw明显激增，即确定为多年冻土上限位置z2，z2为两个相近z值得平均值，若试算出现多个z2，以最大值为准；
[0125] 不断调节时间t可以在试验时间t1模拟得到现场时间t1'内多年冻土上限位置变化规律；
[0126] (3)综上，根据不同公式解析解求得的多年冻土上限位置为
[0127] 根据相似理论，试验所观测到的试样模型冻融变化量与现场原型冻土上限冻融变化量之比也为1:a，故现场原型冻土上限位置为az，根据此公式可以求得现场冻土上限位置变化规律，不断调节时间t可以在试验时间t1模拟得到现场时间t1'内多年冻土上限位置变化情况。
[0128] t可根据具体试验要求规定调节，t1远远小于t1'。
[0129] 需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下。由语句“包括一个......限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素”。
[0130] 尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、
修改、替换和变型，本发明的范围由所附
权利要求及其等同物限定。
