下载文档原格式
土壤入渗实验报告一、实验目的：进行土壤入渗试验，对土壤入渗规律有大致了解，并且利用测的数据绘出土壤累积入渗量和时间的关系曲线，利用该曲线求出入渗强度和时间的关系。
二、实验仪器：直径4.5cm的土柱圆筒、宽5cm长4cm的马氏瓶三、实验步骤：（1）把准备好的土装入圆筒中，每装5cm就夯实一次，直至土的的顶端位于圆筒进水孔的下缘。
（2）检测马氏瓶是否漏气。
如果漏气就换试验设备或者用凡士林涂抹。
（3）将水灌入马氏瓶，把马氏瓶的出水孔和圆筒的进水口用橡胶管连接好，调节圆筒和马氏瓶的相对高度使得马氏瓶的出水孔刚好出水。
（4）读取马氏瓶中水的高度。
（5）实验开始，量取土柱量筒中土壤稳定下渗时土柱的淹水的深度，分别读取实验开始后第1、3、5、7、10、15、20、25、30、35、40min钟时马氏瓶中水的高度。
（6）根据马氏瓶的面积和圆筒的面积求出圆筒中从开始到不同时刻的累积入渗量I，画出土壤累积入渗量I和时间的关系曲线，并利用该曲线求出入渗强度i和时间的关系。
四、实验数据：土壤稳定下渗时土柱的淹水深度为1cm五、数据分析：i/012345651015202530354045入渗时间t/min土壤累积入渗量I 和时间t 的关系曲线图根据土壤累积入渗量I 和时间t 的拟合函数求导得到土壤入渗强度i 和时间t 的关系00.050.10.150.20.250.30.350.40.451020304050入渗时间t/min土壤入渗强度i/cm·min土壤入渗强度i 和时间t 的关系曲线图六、 体会：1. 土壤渗流试验看似简单，不过做成功很难。
实验室中大多数仪器损坏。
实验条件艰难，但是还是尽量利用实验室中的其他条件勉强把实验做完了。
2.土壤入渗要现实土壤足够湿润，然后才能开始记录试验数据，否则会出错。
3.土壤入渗强度开始较强，后来变小，逐渐稳定。
土木工程中的渗透系数测定方法土木工程中的渗透系数测定方法是评估土壤和岩石渗透性的重要手段。
它对于土木工程项目的设计和施工具有重要意义。
准确测定渗透系数可以帮助工程师评估土壤或岩石的渗透性能，以选择合适的建筑材料和采取适当的工程措施，确保工程的稳定和安全性。
一、定负压渗透试验法定负压渗透试验法是一种常用的渗透系数测定方法。
它通过建立负压梯度，测量土壤或岩石中的水渗透速率。
在实验中，首先需要将试验材料制成柱状或土壤墙样品，并保持边界的密封。
接着，在样品上施加一定的负压，测量水渗透过程中的流量并记录时间。
通过对流量和时间的关系进行分析，可以计算出渗透系数。
定负压渗透试验法的优点是操作简单、可控性强，适用于不同类型的土壤和岩石。
二、土柱法土柱法是另一种常用的渗透系数测定方法。
它通过制备土柱样品，浸水并测量土壤中的水渗透速率。
在实验过程中，首先需要选择一定的土壤样品，并保持其湿润，以避免水分蒸发。
接着，在柱状样品的底部设置渗流出口，并将其放置在一个恒定水位的水槽中。
然后观察进水与出水之间的水位差，测量出水的流量和时间。
通过对流量、时间和水位差的关系进行分析，可以计算出土壤的渗透系数。
土柱法适用于渗透系数较小的细颗粒土壤。
三、试井法试井法是一种间接测定渗透系数的方法。
它通过在井孔中进行渗流试验，测量水位变化来推断土壤或岩石的渗透系数。
在试井法中，首先需要选择一个合适的井孔，并在井孔内设置流量计和水位计等仪器。
接着，注入一定数量的水，并观察水位的变化。
通过对水位和时间的变化进行分析，可以间接计算出渗透系数。
试井法适用于渗透系数较大的岩石和土壤。
四、强夯法强夯法是一种特殊的渗透系数测定方法。
它通过制造冲击波，使水分浸透到土壤中，从而估算出土壤的渗透系数。
在强夯法中，首先需要使用特殊的设备进行夯击，使土壤中的颗粒重新排列。
接着，注入一定数量的水分，并观察水分在土壤中的渗透情况。
通过对夯击和渗透过程的观察和分析，可以估算出土壤的渗透系数。
摘要本文通过选取同粒径粘土和砂土装填为不同柱长土柱展开室内低渗透介质中一维土柱溶质运移实验。
土柱按照构成依次编号为：1号.10cm砂；2号.5cm砂+3cm粘土+5cm砂；3号.5cm砂+5cm粘土+5cm砂；4号.5cm砂+8cm粘土+5cm砂。
在土柱顶端持续给予浓度为0.0856mol/L的NaCl溶液，通过稳定每个土柱上下断面的水头，在土柱底端测量土柱渗出液，当渗出液电导率达到峰值时，移除土柱顶端的NaCl溶液，通过自来水淋滤，测量不同时间、不同土柱长度下土柱渗出液的电导率，绘出不同柱长渗出液的电导率和时间关系曲线。
根据观测数据采用图解法（三点公式）利用三个特定浓度点的时间t0.5、t0.84、t0.16，平均孔隙流速u求得弥散系数依次为D L1=2954.34cm2/h、D L2=0.384cm2/h、D L3=0.432cm2/h、D L4=0.577cm2/h；弥散度αL依次为2.653cm、1.825cm、3.470 cm、6.635cm。
然后使用CXTFIT2.1软件拟合观测数据得到弥散系数依次为403cm2/h、0.615m2/h、0.181cm2/h、0.199cm2/h，弥散度依次为0.920cm、2.577cm、2.864cm、4.678 cm。
两种计算方法均显示，随着粘土介质长度的增大，弥散度也随之增大。
然后利用所求参数反演实验数据，在同一坐标系中绘制出C/C0-t曲线。
通过比较三点公式和CXTFIT2.1软件拟合求得的弥散系数，发现后者相对于前者实验精度更高，因为前一种方法由孔隙度或含水量除以渗透流速得到的平均孔隙流速与实际平均孔隙流速存在较大误差，特别对细粒介质的误差很大，不能准确地预测溶质迁移的结果，而后者能充分利用实验观测数据，得到的结果真实可靠。
关键词：溶质运移弥散系数弥散度CXTFIT2.1 穿透曲线ABSTRACTIn this thesis, we have did the solute transport experiments and the leaching experiments in four soil columns with the same particle size and different column lengths. According to the composition of soil columns , soil samples are numbered as: The 1st. 10cm sand; the 2nd. 5cm sand+3cm soil+5cm sand; the 3rd . 5cm sand +5cm soil +5cm sand; the 4th.5cm sand +8cm soil+5cm sand. By stabilizing the soil columns’ hydraulic head , we use the NaCl solution t o simulate the soluble contaminants . Then we get a series of data by measuring the electrical conductivity in the soil columns with the changeable time .By the data we obtained, we draw the BTC.We take two methods to get the dispersion coefficient :The first way is to make the use of three specific concentration point time t0.5，t0.84，t0.16 and the average pore velocity to calculate it. Then we get the dispersion coefficient are D L1=2954.34cm2/h, D L2=0.384cm2/h, D L3=0.432cm2/h, D L4=0.577cm2/h, the dispe rsivity αL are 2.653cm, 1.825cm, 3.470 cm, 6.635 cm. The other way is to enter experimental data into the CXTFIT2.1 programme to obtain dispersion coefficient , and using the inverse technique to get another experimental data. Then we get the diffusion coefficient are 403cm2/h, 0.615m2/h, 0.181cm2/h, 0.199cm2/h, the dispersivity αL are 0.920cm, 2.577cm, 2.864cm, 4.678cm.As shown from the two calculation methods, the dispersivity increases as the soil column length increase. Then analyze the simulation results, the calculated parameters inversion values and the real values. Then we draw the C/C0—t curve in the same coordinate system .For the reason of making full use of the experimental observations, the CXTFIT 2.1 program got a higher accuracy when it is compared with the former way. And we also found that it has a great difference between the D L obtained by the average pore velocity we used and the actual average velocity . The results obtained by porosity and average pore velocity we used were inexact. If we take them as the same one ,we will make a great mistake, especially for fine medium ,so it can’t reflect the results of solute transport accurately.Keywords: s olute transport dispersion coefficient dispersivity CXTFIT2.1 program the break through curve(BTC)目录第一章绪论............................................................................................................................ - 4 - 第二章国内外研究进展........................................................................................................ - 6 - 第一节国内外研究进展.................................................................................................. - 6 - 第二节研究中的问题...................................................................................................... - 7 - 第三节本文的研究思路和内容...................................................................................... - 7 - 第三章实验方案.................................................................................................................... - 8 - 第一节实验仪器.............................................................................................................. - 8 - 第二节实验步骤.............................................................................................................. - 9 - 第三节实验物理参数的测定........................................................................................ - 13 - 第四节实验成果分析方法............................................................................................ - 15 - 第五节实验成果分析.................................................................................................... - 21 - 第四章结论及建议................................................................................................................ - 29 - 第一节本实验得出的结论............................................................................................ - 29 - 第二节实验误差分析及不足........................................................................................ - 29 - 第三节对后续实验的建议............................................................................................ - 30 - 致谢.......................................................................................................................................... - 30 - 参考文献.................................................................................................................................. - 32 -第一章绪论我国是一个水资源相对短缺的国家，人均水资源占有量仅为世界人均占有量的四分之一，水资源己成为我国经济、社会和环境可持续发展的关键制约因素。
弗罗里硅土柱层析分离多氯联苯和有机氯农药的研究1谢飞1 2，周传光2，王艳洁2，姚子伟2，林忠胜2，徐恒振2，吕景才1，赵元凤11大连水产学院，辽宁大连（116023）2国家海洋环境监测中心，辽宁大连（116023）E-mail：xiefei0702@摘要：本文建立了弗罗里硅土(Florisil )柱层析分离多氯联苯（PCB S）与有机氯农药（OCP S）的方法。
比较了不同的Florisil活化时间、活化温度、粒径、用量及层析柱的柱长/柱径对PCB S 与OCP S分离效率的影响。
结果显示，当弗罗里硅土活化时间为30min，活化温度为200℃、粒径为60目～100目、用量为2.0g及其柱长/柱径为230mm/8mm时,分离效果较好，其PCB S 回收率在76.0%～101.1%。
仅有七氯、艾氏剂、DDE、狄氏剂、甲氧-DDT没有分开，但对PCB S干扰较小。
关键词：Florisil层析柱，PCBs和OCPs，GC-ECD，分类分离中图分类号：X551 引言多氯联苯（PCBs)和有机氯农药(OCPs)均属于持久性有机污染物(POPs) ,与常规污染物不同,它们在自然环境中极难降解,并能在全球范围内长距离迁移,被生物体摄入后不易分解,并沿着食物链浓缩放大,对人体危害巨大[1]。
目前环境样品中多氯联苯（PCBs)和有机氯农药(OCPs)的痕量分析，一般采用带电子捕获检测器的毛细管气相色谱法（GC-ECD）。
但是在分析时，这两类污染物之间互相干扰，影响各自的准确测定。
因此，有关PCBs和OCPs的分离，一直是环境分析工作者关注的热点[2-9]。
采用传统的柱层析分离法，操作简单、成本低、重复性好，分离效率还可以。
徐恒振、周传光等[10-12]曾采用Florisil柱层析较好地分离了海洋环境样品中的PCBs和OCPs。
本文在此基础上，对 Florisil柱层析分离PCBs和OCPs的影响因素做进一步的探讨，旨在通过不同的Florisil活化时间、活化温度、粒径、用量及其柱长/柱径等，对PCB S与OCP S分离效率的影响进行研究，以期得到柱层析的最佳分离条件，为海洋环境样品中PCBs和OCPs的检测提供参考资料。
以土的常三轴实验学习Duncan-Chang本构关系模型一、实验过程1、试样制备试验土样取自于南水北调焦作段一处工程，取回后，人工制成含水量15%的土体。
在实验制样过程中，由于含水量较高，所以在通过制样器后，土柱未能成型，于是在原来土样基础上，添加了较干的土，再在制样器侧壁涂抹凡士林。
最后制成高度7厘米，直径3.5厘米的土柱实验样品2、不固结不排水（UU）剪切试验试验是在土木工程学院深部矿井重点实验室进行的，试验装置如图1所示。
图1 常三轴实验仪主要试验步骤为(1)记录体变管的初始读数；(2)对试样加周围压力，并在周围压力下固结。
当孔隙水压力的读数接近零时，说明固结完成，记下排水管的读数；(3)开动马达，合上离合器，按0.0065%/min的剪切应变速率对试样加载。
按百分表读数为0，30，60，90，120，150，180，210，240，300，360，420，480，540，600，660，⋯，的间隙记读排水管读数和量力环量表读数，直到试样破坏为止。
二、邓肯张双曲线模型到目前为止，国内外学者提出的土体本构模型不计其数，但是真正广泛用于工程实际的模型却为数不多，邓肯-张模型为其中之一。
该模型是一种建立在增量广义虎克定律基础上的非线性弹性模型，可经反映应力～应变关系的非线性，模型参数只有８个，且物理意义明确，易于掌握，并可通过静三轴试验全部确定，便于在数值计算中运用，因而，得到了广泛地应用。
1、邓肯-张双曲线模型的本质邓肯-张双曲线模型的本质在于假定土的应力应变之间的关系具有双曲线性质，见图2（a)。
图2(a ） 12()~a σσε- 双曲线图2（b) 1131/()~εσσε-关系图2 三轴试验的应力应变典型关系理论图1963年，康纳（Kondner ）根据大量土的三轴试验的应力应变关系曲线，提出可以用双曲线拟合出一般土的三轴试验13()~a σσε-曲线，即aab a εεσσ+=-31 （1）其中，,a b 为试验常数。
渗透系数测定方法渗透系数是指土壤中水分向地下渗透的能力，也可以理解为土壤对水分传输的难易程度。
渗透系数的测定方法有很多种，下面我将介绍几种常用的测定方法。
1. 饱和渗透法饱和渗透法是一种简单常用的渗透系数测定方法。
首先在一个圆柱形的实验室土柱上布置一个封闭水槽，确定土柱的渗透面积。
然后，将所需测定的土样填充到土柱内，使其饱和。
在土样上方设置一个水箱，控制水头高度，然后记录柱水头变化的时间和水头下降量。
根据达西-利胡特方程，可以通过计算水流速度和水头梯度来计算渗透系数。
2. 不饱和渗透法不饱和渗透法适用于土壤不饱和状态下的渗透系数测定。
实验中需要使用一个均匀细腻的土壤样品，将其压实到一个封闭的渗透仪器中，确保土壤样品不受外界水分输入。
然后，通过应用一系列确定水力头梯度和水流速度的方法，测定不同土壤水分含量下的水流速度，进而计算出渗透系数。
此方法可以模拟不同土层中的不同水分含量下的渗透过程。
3. 稳态渗透法稳态渗透法是一种适用于渗透系数测定的常用方法。
实验中使用一个圆筒形土壤样品，通过施加稳定的水头梯度，在不同的水头条件下测量渗透流量。
根据达西-利胡特方程，通过测量的水流量和水头梯度，可以计算出渗透系数。
此方法适用于可测量渗透流量且渗透系数变化较小的土壤。
4. 双线性渗透法双线性渗透法是一种适用于非饱和土壤的渗透系数测定方法。
该方法要求样品含水量在饱和状态下与在不饱和状态下之间存在明显的差异。
实验中使用两个不同水头梯度的渗透仪器，同时对同一土样进行测试。
通过测量两个试验的渗透流量和水头梯度，可以计算渗透系数。
此方法适用于非饱和土壤和含水量变化较大的土层。
总结起来，渗透系数测定方法包括饱和渗透法、不饱和渗透法、稳态渗透法和双线性渗透法等。
选择合适的方法需要考虑土壤的性质、实验设备的条件以及测定的目的和要求。
通过测定渗透系数可以帮助我们了解土壤的水分运移特性，对于农田水分管理和水文模型的建立都有重要的应用价值。
《土柱实验》教学指导书
一、实验目的
1、加强对黄土特性的学习，掌握土壤弥散系数的概念和方法；
2、观察污染物在土壤中的迁移情况；
3、深入了解土壤对污染物的吸附和解析模型；
4、熟悉模拟实验数据的记录和处理，是将来进行污染物运移数学模拟的基础。
二、实验内容
1、测定土壤的弥散系数
2、观察污染物在土壤中的迁移
3、建立土壤对污染物的等温吸附模型和吸附动力学模型；
4、建立土壤对污染物的等温解析模型和解析动力学模型；
三、仪器相片
活性碳夹层柱身出水口出水或入水孔
四、模拟实验土柱和装置
1、土柱的概述
实验土柱通常分为两种：非扰动的原状土柱和填装的土柱。
原状土柱较好的代表了研究对象的原土壤结构及理化性质。
填装土柱可以是经过筛
分的一种土壤,也可以是按一定比例混合的几种土壤。
填装土柱不能保持原土壤结构,但是,适用于一定目的的专门研究。
例如,通过改变各种参数(土壤粒径、比表面和矿物组分,土壤水的pH值,或土壤的不同混合比例等)，来研究土壤对特定污染物的吸附性能。
2、土柱的尺寸
土柱的尺寸没有硬性的规定，一般认为圆形土柱的长度大于其直径的2.5倍即可。
不过对于要求精确的模拟土柱，其尺寸还要考虑以下几个方面
(1) 污染物通常含有多种物质，在设计实验的过程中应当以迁移能力最强的污染物为标准来设计土柱，一般来说迁移能力上，有机污染物<阳离子污染物<阴离子污染物。
(2) 由于在土柱内建立稳定的非饱和流场的方法不同,当土柱下段含
水率偏高，不能满足实验要求时，必须要增长土柱的长度。
(3) 对于不同的实验目的，必须根据测量方法和测量数据来确定土柱长度与直径的比值。
3、入水方式
入水方式通常有从上往下渗水和从下往上饱水两种。
通常前者用于，地表污染物往土壤中的下渗，可以比较正确的反应该条件下污染物在土壤中的迁移情况，但是无法用于研究土壤对其的吸附作用。
而从下往上饱水，可以让土壤与污染物充分接触，是研究土壤对污染物吸附机理和模型的最佳方式。
一般情况下，土柱模型都采用定流量水，在要求并不是很高的实验中，通常使用双马氏瓶来固定入水水头，但是在精度很高或者对要求很大，很小或很稳定的水流量时，必须使用蠕动泵入水。
4、土柱的注意事项
(1) 砾石层不管入水方式如何，必须在入水口的土壤前铺设1-2cm厚的砾石层，防治水流对土壤层的冲刷，从上往下渗水时尤其必须注意。
而在出水口土壤后也必须铺设1-2cm的砾石层，防治土壤堵塞。
(2) 尼龙丝网任何不同材料的层位之间，都应该放置尼龙丝网，且必须根据流向和材料来选择尼龙丝网的目数，当顺流向且材料由细变粗时，尼龙网需要目数足够大以阻挡细材料往粗材料中的渗透。
(3) 反应渗透墙根据实验的目的和所研究的污染物种类，可以适当的在土柱中添加反应渗透墙。
反应渗透墙的材料组成厚度都必须经过详细的计算以选取合适的数值。
(4) 柱身出水处首先为了防止土壤对柱身处水口的堵塞，必须在出水口前布置合适的尼龙丝网；其次在不同的实验中，此出水口有不同的功能，例如水头的测量和水质电导率的测量。
(5) 土柱内侧必须粗糙，以防止水沿内测管壁的优先流。
5、附属设备
马氏瓶、蠕动泵、橡皮管、示踪剂、电导率仪及探头和收集样品用的容器。
五、实验步骤
1、设计实验
确定所用污染物溶液、土柱尺寸、填充物类型和反应渗透墙情况。
2、装填土柱
装填前，需先测定所填土壤的部分系数，如粒径大小、级配组成、矿物组成、孔隙度、干容重γ和含水量S。
装填时采用干堆法，且要利用土壤的干容重γ和含水量S进行计算。
根据不同的土柱尺寸，每次装填的高度在H须在5-10cm之间，而每次装填的重量为：
W=V·γ(1 +S)
式中：
V——每次装入土体体积(cm3);
γ——天然土体干容重(g/cm3);
S——室内土壤含水量(%);
在每次装入W(g)土壤后,利用压实器进行土壤压实,使其达到规定的土壤土柱高度H。
如此试验土柱的干容重与天然情况下土壤干容重将相等或者接近,由
此来控制土柱筒内土壤的孔隙率与天然情况下的土壤孔隙率基本相符。
3、准备外设
根据设计的入水方式，准备马氏瓶或蠕动泵且将其连通；准备容器收集出水口的样品；对于柱身的出水处，根据计划准备电导率或水压表或直接采样。
4、记录数据
根据不同的实验目的记录数据，具体内容见下章。
5、结果分析
对记录数据进行处理，分析实验结果；然后需要对数据和结果进行验证，最后得到结论及建议。
具体内容见下章。
六、实验结果处理
（一）测定试验土柱的弥散系数
1、示踪剂的选择
一般条件下，试验中采用氯化钠作为示踪物质,以氯离子作为测定土的示踪剂。
虽然氯离子对于土壤胶体来讲属于惰性离子,但由于阳离子的其仍具有一定吸附作用，为了合理地确定出土柱体的动力弥散系数，须首先进行氯离子的静态平衡吸附试验，以确定土壤对它的最大吸附量及平衡液浓度。
2、土壤对氯离子的吸附试验
称量20g风干的土壤盛入三角瓶内,再加入200ml氯离子浓度分别为200ppm，400ppm，600ppm，800ppm，1000ppm，1200mp的溶液(水土比为10：1)，在水浴振荡器上振荡5个小时，静置20小时，用离心机离心澄清，取上清液分析其平衡液浓度，用差减法求得土壤胶体对氯离子(CI-)的吸附量。
绘制土壤对氯离子吸附平衡曲线。
该实验可以得到，土壤对氯离子的最大吸附浓度X。
在下面进行的土柱实验中，先用X浓度的氯离子溶液使土柱中突然达到饱和，然后用2/3X浓度的氯离子溶液进行弥散系数的测定。
3、流体动力弥散系数的测定
在土柱入口连续恒定地注入氯离子浓度为1500mg/l，出口取样测定氯离子的浓度。
以出口处氯离子浓度C为Y轴，相应时间为X轴绘制穿透曲线。
根据穿透曲线查找以下数值：即当C/C0=0.16，0.5，0.84时的时刻：t0.16；t0.5；t0.84，将时间t0.5的数据代入下列公式
v=L/t0.5
V——流速；
L——土柱长度；
得到速度V，然后带入
DL=v3(t0.84-t0.16)2/8L
DL——弥散系数
这样就能得到土壤的弥散系数。
（二）建立土壤对污染物的等温吸附(解析)模型和吸附(解析)动力学模型
1、实验过程
吸附过程土柱先用蒸馏水饱和后，开始加入初始含污染物溶液，调节该溶液流速为Y进行淋滤试验，按一定间隔时间采取淋出液若干毫升进行污染物含量测定.并记录淋入、淋出液的体积。
当淋出液污染物浓度接近或达到初始溶液污染物浓度时，可以近似认为污染物的吸附过程达到饱和状态，停止淋滤。
解吸过程将上述污染物吸附达到饱和状态的试验土柱搁置36 h之后，加入与初始溶液pH值相当的蒸馏水进行污染物的淋溶(流速与上面实验需要一致)解吸试验过程。
同样间隔时间取淋出液样分析污染物含量，并记录当时淋入、淋出液体积。
当淋入液的体积与上述吸附过程停止淋滤时所加初始溶液的体积相同时，即可停止污染物的解吸淋溶试验。
2、建立等温吸附(解析)模型
数据处理的基本思路在于，分析数据之间的对应关系，然后应用适当的方程进行拟合。
对于等温吸附(解析)模型，现在有三种理论可以用来解释，分别是：
Henry曲线：G=K d C
式中：G —— 吸附量；
C —— 平衡时液相离子浓度（mg/L ）；
K d ——吸附平衡常数（L/g ）。
Frendlich 曲线 ：G=KC n
式中：G —— 吸附量；
K —— 常数；
n —— 表示该等温吸附线线性度的常数，介于0与1之间；当液相中被吸附组分浓度很低，或在砂土（CEC 值小）中产生吸附时，n →1；
C —— 平衡时液相离子浓度（mg/L ）；
Langmiur 曲线： KC
1KC S S m ++= 式中：S m —— 某组分的最大吸附浓度（mg/kg ）；
K —— 与键能有关的常数；
实验完成后，将淋出液中溶质的平均浓度C 为X 轴，溶质在土壤中的单位吸附量(单位残留量)G 为Y 轴，绘制等温吸附(解析)曲线，对该曲线分别进行上述三种模型的模拟，最后根据模拟的结果进行分析。
3、建立吸附(解析)动力学模型
吸附(解析)动力学模型，通常有下面四种，分别是：
Elovich 方程：c=a+b ·lnt
双常数速率方程：lnc=a+b ·lnt
抛物线扩散方程：c=a+bt 1/2
一级反应动力学方程：lnc=a+bt
其中：
c 为溶质在单位土体上的吸附量(mg/kg)；
t 为反应时间(h)；
a 为与初始浓度有关的试验常数；
b 为与吸附活化能有关的吸附速率常数。
实验完成后，将淋出液时间t 为X 轴，溶质在土壤中的单位吸附量(单位残留量)c 为Y 轴，绘制吸附(解析)动力学曲线，对该曲线分别进行上述四种模型的模拟，最后根据模拟的结果进行分析。
