第二章 土壤水分运动基本方程2
第二章 土壤水分运动基本方程2
第二章 土壤水分运动基本方程 如前所述，达西定律是由达西（Darcy ，Henry 1856）通过饱和砂柱渗透试验得出，后由Richards （1931）将其扩伸至非饱和水流中，并规定导水率为土壤负压h 的函数，即 ()H h k q ?= （2-2-1） 式中：H ?——为水势梯度； k （h ）——为导水率，是土壤负压h 的函数； q ——为水流通量或流速。 Richards 方程垂向一维方程为 ) 1)(( ) (±??-=??-=z h k z H k q z θθ 注意：H=h ±z ，垂直坐标向上为“+”；向下时为“–”。 由于k （h ）受滞后影响较大，上式仅适用于单纯的吸湿或脱湿过程。若将导水率作为容积含水率函数，即以k （θ）代替人k （h ），则可避免滞后作用的影响。 一般说来达西定律对饱和与非饱和水流均可适用，即水流通量与势能梯度成正比。但在饱和土壤中，压力为正值，其总水头包括了由该点在地下水面以下深度来确定的静水压力（正值）和相对于基准面高度来确定的位置水头，总水头为压力水头和位置水头之和，水由总水头高处向低处流动。在非饱和土壤中，基质势为负值，土水势在不考虑溶质势、温度势及气压势时，只包括重力势和基质势。因此，总水头常以负压水头和位置水头之和来表示。 一维Richards 方程的几种形式： 根据() ()θθ θD h k =??（K=C ×D ）得： x h k q x ??-=)(θ x D q x ??-=θ θ)( y h k q y ??-=) (θ y D q y ??-=θθ)( )1)( (±??-=z h k q z θ )]()([θθθk z D q z ±??-=
土壤入渗速度测定实验
实验一 土壤入渗速度的测定实验 一、实验目的 1.测定特土壤的垂直入渗特性曲线。 2.掌握测定土壤吸渗与入渗速度的操作方法。 二、实验原理 考斯加可夫公式:i t ＝i 1t -a ---------------------------- (1) i t ——入渗开始后时间t 的入渗速度; i 1——在第一个单位时间土壤的渗透系数,相当于t ＝l 时的土壤下渗速度; a —指数。 对公式(1)取对数得 lgi t ＝lgi 1-a·lgt ----------------------- (2) 实测的lgi t ,lgt 点应成直线关系,取t=1时的i 值,极为i 1,该直线的斜率为a 值。 计算时t a ,t b 时刻对应i a ,i b ,代入下式得 b a b a t t i i a lg lg lg lg --= ----------------------- (3) 若已知i 1,a 值也可以按下述方法推求,有式(1)积分得 a t a t t a i dt t i idt I ---= ==??110 10 1 ----------------------- (4) I 为时间t 内总入渗量(累积入渗量),由实测数据得出,由于i 1已知,故a 可以求出。该法的缺点时很难测定第一个单位时间的入渗强度。 三、实验设备 1.土壤入渗仪:一套; 2.秒表:一只 3.量筒、滤纸、烧杯 4.排水管 5.接渗瓶 四、实验步骤 1.装土:将玻璃管从入渗仪上取下,底部放入一片滤纸,然后装土,在装土期间,
用木棒稍捣,要求土样均匀,装土至玻璃管即可,再在土样上部放入一张滤纸,把玻璃管与入渗仪连接好。 2.加水:关闭水阀,打开排气阀,用烧杯向加水槽加水,使量桶里的水位到达到一定刻度处,然后关闭排气阀。 3.建立水头开始实验:用烧杯迅速向玻璃管加水至玻璃管上标线,水头建立后,立即打开供水阀,同时打开秒表计时,三者要求同时进行,动作要迅速、准确、细心。 4.记数:实验开始后秒表不能中断,要求每隔1分钟1次,共读10次,再每隔2分钟读1次,共读10次,再每隔3分钟读1次,共读5次,以后每隔5分钟读1次,直到两相邻时段内,读数差值相等,说明土壤入渗已经达到稳定,即停止实验,记录项目为记录表中的第l项与第2项。 土壤非饱与垂直入渗率测定表 日期: 土质: 垂直入渗仪横断面面积(mm2): 马氏瓶横断面面积(mm2): 五、实验资料整理 1.根据实验数据,将记录的马氏瓶读数算为毫升,再计算为水层深度。 2.计算时段平均入渗速度。
土壤的入渗特性及渗吸速度测定_灌排工程学
第一部分 课程实验及指导 实验一：土壤的入渗特性及渗吸速度测定 一、实验目的 土壤渗吸速度是反映土壤透水性能的重要指标，它是农田水量平衡计算的重要依据。旱田在进行地面灌溉时，灌溉水在重力作用下自地表逐渐向下湿润。为保证最有效地利用灌溉水，既要使计划湿润层得到均匀的灌溉 水，又不产生多余的水量向深层渗漏，必须了解水向土中入渗的规律。 二、实验设备 渗吸速度测试仪、量杯、秒表等。 三、实验过程 1．取自然风干土碾碎过筛，要求碎块不大于2毫米，测筒底铺滤纸，装土至给定深度，适当沉实，再盖滤纸。 2．在量杯内灌水，并关闭放水管和通气管（如图所示），放在支架上。 3．实验开始时同时完成：掀动计时秒表，迅速使测试仪中土样上建立水层2厘 米。 图1-1-1土壤入渗特性实验装置 4．实验开始后，定时记载量杯中水量读数，时间间隔初期较短，以后逐渐加大。并填写表1-1-1： 表1-1-1 土壤入渗特性测定记录表 四、实验原理 在地面形成一定水层的入渗称为有压入渗，对于均质土的入渗强度，已有若干计算公 式，菲利普根据严格的数学推导，求的解析解为： f i t s i += -2/12 （1-1-1）
i —t 时刻的入渗强度； s —与土壤初始含水率有关的特性常数，称为吸水率； i f —稳定入渗率，即饱和土壤渗透系数。 考斯加可夫根据野外实测资料分析，发现入渗强度（渗吸速度）与时间之间呈指数关系，其形式为： α-=t i i 1 （1-1-2） 式中 i 1—第一个单位时间的入渗强度； α—反映土壤性质与入渗初始时土壤含水率的经验常数。 饱和与非饱和土壤水分运动均服从达西定律，所不同者，在饱和情况下，认为渗透系数是常数；而在非饱和情况下，渗透系数是变量，其值随土壤含水率而异，含水率越低，渗透系数越大。 五、实验要求 1．根据水室断面和测筒断面，求出△t 时间内测筒下渗的水量。 2．求出各时段平均入渗速度v 。 3．用坐标纸点绘渗吸速度随时间变化过程线。 4．分析确定供水开始时土壤渗吸速度i f 、渗吸系数及透水指数α值。 5．填写实验报告。 六、思考题 利用菲利普公式和考斯加可夫公式求s 或i 1时，讲选取第一个单位时刻的i 值，如何理解这第一个单位时刻的意思？它是根据i 的取值单位还是绘图时的取值单位？
土壤入渗理论与方法(改)
西南林业大学 硕士研究生文献综述 论文题目：土壤入渗理论与方法 学院：环境科学与工程学院 年级： 2014级 成员：冯晓月阮书鹏曹向文 指导教师：宋维峰 2015年4 月25 日
摘要 入渗是水文学中重要的基本概念，定量确定土壤入渗性能对认识水循环及水利用具有重要的理论意义和实践价值。当然，从不同角度出发去探讨土壤入渗也有不同的科研意义。本文试图通过对目前国内外对土壤入渗的研究做一个系统性归纳与对比，从而为下一步的学术论文打下基础。 关键词 土壤入渗；方法；模型；影响因素
Abstract Infiltration is a vital basic concepts in hydrology, qualitatively analysis soil infiltration capability has important theoretical significance and practical value in the water cycle and use. Of course, using different angle of view to discuss soil infiltration also have different research significance. This article attempts to do a systematic induction and comparison through the study of soil water infiltration at home and abroad, which lays the foundation for the next academic paper. Key words soil infiltration;methods; the influence of factors
数学物理方程第三版第一章答案(全)
数学物理方程第三版答案 第一章． 波动方程 §1 方程的导出。定解条件 1．细杆（或弹簧）受某种外界原因而产生纵向振动，以u(x,t)表示静止时在x 点处的点在时刻t 离开原来位置的偏移，假设振动过程发生的张力服从虎克定律，试证明),(t x u 满足方程 ()?? ? ??????=??? ??????x u E x t u x t ρ 其中ρ为杆的密度，E 为杨氏模量。 证：在杆上任取一段，其中两端于静止时的坐标分别为 x 与+x x ?。现在计算这段杆在时刻t 的相对伸长。在时刻t 这段杆两端的坐标分别为： ),();,(t x x u x x t x u x ?++?++ 其相对伸长等于 ),()],([)],([t x x u x x t x u x t x x u x x x ?+=??-+-?++?+θ 令 0→?x ，取极限得在点x 的相对伸长为x u ),(t x 。由虎克定律，张力),(t x T 等于 ),()(),(t x u x E t x T x = 其中)(x E 是在点x 的杨氏模量。 设杆的横截面面积为),(x S 则作用在杆段),(x x x ?+两端的力分别为 x u x S x E )()(x u x x S x x E t x )()();,(?+?+).,(t x x ?+ 于是得运动方程 tt u x x s x ???)()(ρx ESu t x =),(x x x x x ESu x x |)(|)(-?+?+ 利用微分中值定理，消去x ?，再令0→?x 得 tt u x s x )()(ρx ?? = x ESu () 若=)(x s 常量，则得 22)(t u x ??ρ=))((x u x E x ????
数学物理方程谷超豪版第二章课后答案
第 二 章 热 传 导 方 程 §1 热传导方程及其定解问题的提 1. 一均匀细杆直径为l ，假设它在同一截面上的温度是相同的，杆的表面和周围介质发生热交换，服从于规律 dsdt u u k dQ )(11-= 又假设杆的密度为ρ，比热为c ，热传导系数为k ，试导出此时温度u 满足的方程。 解：引坐标系：以杆的对称轴为x 轴，此时杆为温度),(t x u u =。记杆的截面面积4 2 l π为S 。 由假设，在任意时刻t 到t t ?+内流入截面坐标为x 到x x ?+一小段细杆的热量为 t x s x u k t s x u k t s x u k dQ x x x x ????=???-???=?+221 杆表面和周围介质发生热交换，可看作一个“被动”的热源。由假设，在时刻t 到t t ?+在截面为 x 到x x ?+一小段中产生的热量为 ()()t x s u u l k t x l u u k dQ ??-- =??--=11 1124π 又在时刻t 到t t ?+在截面为x 到x x ?+这一小段内由于温度变化所需的热量为 ()()[]t x s t u c x s t x u t t x u c dQ t ????=?-?+=ρρ,,3 由热量守恒原理得： ()t x s u u l k t x s x u k t x s t u c x t ??-- ????=????11 2 24ρ 消去t x s ??，再令0→?x ，0→?t 得精确的关系： ()11 224u u l k x u k t u c -- ??=??ρ 或 ()()11 22 2112244u u l c k x u a u u l c k x u c k t u --??=--??=??ρρρ 其中 ρ c k a =2 2. 试直接推导扩散过程所满足的微分方程。 解：在扩散介质中任取一闭曲面s ，其包围的区域 为Ω，则从时刻1t 到2t 流入此闭曲面的溶质，由dsdt n u D dM ??-=，其中D 为扩散系数，得 ?????= 2 1 t t s dsdt n u D M 浓度由u 变到2u 所需之溶质为 ()()[]???????????ΩΩΩ ??=??=-=2 12 1121,,,,,,t t t t dvdt t u C dtdv t u C dxdydz t z y x u t z y x u C M 两者应该相等，由奥、高公式得： ????????Ω Ω??==????????? ??????+???? ??????+??? ??????=2 12 11t t t t dvdt t u C M dvdt z u D z y u D y x u D x M 其中C 叫做孔积系数=孔隙体积。一般情形1=C 。由于21,,t t Ω的任意性即得方程： ?? ? ??????+???? ??????+??? ??????=??z u D z y u D y x u D x t u C 3. 砼(混凝土)内部储藏着热量，称为水化热，在它浇筑后逐渐放出，放热速度和它所储藏的 水化热成正比。以()t Q 表示它在单位体积中所储的热量，0Q 为初始时刻所储的热量，则 Q dt dQ β-=，其中β为常数。又假设砼的比热为c ，密度为ρ，热传导系数为k ，求它在浇后温度u 满足的方程。 解： 可将水化热视为一热源。由Q dt dQ β-=及00Q Q t ==得()t e Q t Q β-=0。由假设，放 热速度为 t e Q ββ-0 它就是单位时间所产生的热量，因此，由原书71页，(1.7)式得 ??? ? ??-=+??? ? ????+??+??=??-ρρββc k a e c Q z u y u x u a t u t 20222222 2 4. 设一均匀的导线处在周围为常数温度0u 的介质中，试证:在常电流作用下导线的温度满足微分方程 ()2201224.0ρω ρωρc r i u u c P k x u c k t u +--??=?? 其中i 及r 分别表示导体的电流强度及电阻系数，表示横截面的周长，ω表示横截面面积，而k 表示导线对于介质的热交换系数。 解：问题可视为有热源的杆的热传导问题。因此由原71页(1.7)及(1.8)式知方程取形式为
土壤入渗实验指导书
《水文学原理》实验指导书 天津农学院水利工程系 2006.9
实验一土壤渗透系数的测定 [实验目的]： １．掌握土壤下渗的物理过程及下渗机理； ２．测量土壤渗透系数Ｋ； ３．学习正确使用渗透筒。 [实验原理]： 下渗过程一般划分为三个阶段。第一阶段为渗润阶段，这阶段，土壤含水量较小，分子力和毛管力均很大，再加上重力的作用，所以此时土壤吸收水分的能力特别大，以致初始下渗容量很大，而且由于分子力和毛管力随土壤含水量增加快速减小，使得下渗容量迅速递减。第二阶段为渗漏阶段，土壤颗粒表面已形成水膜，因此分子力几乎趋于零，这时水主要在毛管力和重力作用下向土壤入渗，下渗容量比渗润阶段明显减小，而且由于毛管力随土壤含水量增加趋于减小阶段，所以这阶段下渗容量的递减速度趋缓。第三阶段为渗透阶段，在这一阶段，土壤含水量已达到田间持水量以上，这时不仅分子力早已不起作用，毛管力也不再起作用了。控制这一阶段下渗的作用力仅为重力。与分子力和毛管力相比，重力只是一个小而稳定的作用力，所以在渗透阶段，下渗容量必达到一个稳定的极小值，称为稳定下渗率。 [实验仪器]： 1.渗透筒（渗透环）一套——渗透筒是用金属做的一套无底同心圆柱筒，筒底 具刀口，同心环内管的横截面积为1000cm2,内径35.8cm，高30-50cm，外筒内径60cm（亦可用土埂围堰代替外筒）； 2.量筒500ml和1000ml各一个； 3.水桶2个；温度计1支（刻度0-50℃）；秒表（普通钟表）1块；量水测针或 木制厘米尺一个；席片或塑料薄膜（灌水时防止冲刷用）。 [实验步骤]： 1.选取具有代表性的地块，把渗透筒的内筒插入土中，深度10cm左右，同时插 好外筒。如无外筒，可筑埂围堰，高度和内筒高相平，埂顶宽20cm，并捣实之。 2.同内外插入量水测针或木制厘米尺各一支，筒内水层厚度一般保持5cm。 3.把席子或塑料薄膜放入筒底，同时把温度计插入筒内。在开始灌水时，土壤 吸水速度较快，为使筒内达到一定水层，第一次灌水要快，同时视水层下降
推求Gardner-Russo持水曲线模型参数的简单入渗法
基金项目 作者简介 湖北监利人博士生 主要研究方向为地面水资源与地下水资源及环境 推求持水曲线模型参数的简单入渗法 薛绪掌张仁铎 武汉大学水资源与水电工程科学国家重点实验室湖北武汉 国家农业信息化工程技术研究中心 北京 中山大学环境科学与工程学院广东 广州 摘要本文基于水平一维非饱和土壤水分运动规律 推求了 用模拟的结果进行拟合其决定系数 为 利用数值模拟数据和实验数据检验该方法将用此方法 结果表明本研究所求得的参数有较高的精度关键词土壤水分渗流运动 参数数值模拟 等和 直接测量土壤水分特征曲线和非饱和土壤导水率的方法 土壤水力特性土壤质地资料被成功地用来预测非饱和土壤水力特性 和 描述非饱和土壤水力特性模型中的参数该方法是在假 等模型中的参 法来推求更多描述土壤水分运动模型中的参数基本理论 水平一维非饱和土壤水分运动建立在 其表达式为
其表达式如下 是土壤饱和体积 方程描述如下 式中 其初始和边界条件其中为土壤初始体积含水率 式中为任意位置 则湿润峰位置的土壤水基质势 很低有 式中为湿润峰 由于 可得 变化的函数表达式 该式右边须乘以一个 参数 式中 当时入渗通量 和湿润峰厚度
此公式相似于表征的水分入渗模型 当土壤湿润峰为时其相对应的土壤水累积入渗量 其中 其中 其中 土壤饱和水 力传导度和土壤饱和体积含水率取 风干土含水量和 根据入渗率和湿润峰之间的关系 ?利用迭代法求得可得参数 和
本研究应用程序 模拟中用到列出了土柱长度为 用到 分别为和 体积含水率和土壤饱和导水率和后借助求解 和为了验证所推导的计算 算得到的参数值和输入的参数值进行比较和参数敏感性分析并将参数估计值代入模型中得到的 表土壤类型水力特性参数 土壤??? 实验方法年 验室温度控制在土壤为风干散装土系采自北京昌平小汤山国家精准农业基地个土壤剖面层 次和土壤样品自然风干且过表 表供试土壤的基本性质 土壤剖面层次深度有机质团粒结构状况 将各种供试土壤按照装土容重分成 在实验前取自然分干土样利用烘干法测定供试土壤的重 和初始体积含水率 在实验室内进行了传统的的有机 试验土柱是界面直径为将供试土壤按设 计容重分层均匀装入圆筒在实验过程中 结果和讨论 数值分析湿润峰 为了验证其结果将模拟结果点绘在二维坐标中和图 图和图分别描述了 其拟合结果见表
数学物理方程第二版答案解析(平时课后知识题作业任务)
数学物理方程第二版答案 第一章． 波动方程 §1 方程的导出。定解条件 4. 绝对柔软逐条而均匀的弦线有一端固定，在它本身重力作用下，此线处于铅垂平衡位置，试导出此线的微小横振动方程。 解：如图2，设弦长为l ，弦的线密度为ρ，则x 点处的张力)(x T 为 )()(x l g x T -=ρ 且)(x T 的方向总是沿着弦在x 点处的切线方向。仍以),(t x u 表示弦上各点在时刻t 沿垂直于x 轴方向的位移，取弦段),,(x x x ?+则弦段两端张力在u 轴方向的投影分别为 )(sin ))(();(sin )(x x x x l g x x l g ?+?+--θρθρ 其中)(x θ表示)(x T 方向与x 轴的夹角 又 . sin x u tg ??=≈θθ 于是得运动方程 x u x x l t u x ???+-=???)]([22ρ∣x u x l g x x ??--?+][ρ∣g x ρ 利用微分中值定理，消去x ?，再令0→?x 得 ])[(2 2x u x l x g t u ??-??=??。 5. 验证 2 221),,(y x t t y x u --= 在锥2 22y x t -->0中都满足波动方程 222222y u x u t u ??+??=??证：函数2221),,(y x t t y x u --=在锥2 22y x t -->0内对变量t y x ,,有
二阶连续偏导数。且 t y x t t u ?---=??-2 3 222)( 22 52222 32222 2) (3) (t y x t y x t t u ?--+---=??- - )2()(2 2223 222y x t y x t ++?--=- x y x t x u ?--=??- 23 222)( ()() 225222232222 23x y x t y x t x u - ---+--=?? ( )()222 252222y x t y x t -+- -=- 同理 ()()222 25 2222 22y x t y x t y u +---=??- 所以 ()() .22 22 2225222222 2t u y x t y x t y u x u ??=++--=??+ ??- 即得所证。 §2 达朗贝尔公式、 波的传抪 3.利用传播波法，求解波动方程的特征问题（又称古尔沙问题） ??? ? ???==??=??=+=-).()(0022222x u x u x u a t u at x at x ψ? ())0()0(ψ?= 解：u(x,t)=F(x-at)+G(x+at) 令 x-at=0 得 )(x ?=F （0）+G （2x ） 令 x+at=0 得 )(x ψ=F （2x ）+G(0)
土壤入渗实验报告
一、实验目的 1．加深对土壤渗吸速度变化的一般规律的了解。 2．了解土壤质地对土壤渗吸速度的影响。 3．掌握土壤渗吸速度的常规测定方法及装置原理。 二、实验设备 水在土壤中入渗分为有压入渗和无压入渗。如漫灌、畦灌和沟灌都属于有压入渗。喷灌、滴灌属于无压入渗。本试验是模拟有压入渗条件下，土壤渗吸速度的测定。 本试验为室内试验，试验装置如图4-1-1。试验仪器大体分为由两部分，即试样渗吸桶和供水马氏瓶。双环入渗试验的外环外径为15cm，内径14cm；内环的外径直径10cm，内径直径9cm，高15cm。安装后要求内环环顶端与渗吸筒齐平，下端插入土内10cm。试验桶正上方为自动供 水箱（即为马氏瓶），使内环保持稳定的水层深度。供水马氏瓶外径6cm，内 径5cm。此外再配备秒表、水桶、水勺和刮土板等试验用具。 三、实验方法及步骤 1．实验准备工作 a．人员分工 每组实验人员3～5人，其中一人计时兼指挥，一人读取供水水位数 值，一人加水，其余人员做记录和观察渗吸规律。 b．准备工作 和内环一并称重， （1）测量试样桶容积V，按欲模拟土壤干容重 干 M。 计算出干土重' （2）将筛网贴紧桶底铺好，然后开始填装。土样一般分5～6次填装， 均匀夯实，层间要“打毛”。土样全部装好后用刮板刮平表面，最后将马 氏瓶安装好待用。 (3) 关闭供水箱（马氏瓶）的出水口，向水箱内注水，然后用胶塞密 封注水进水口。图4-1-1 试验装置示意图 (4) 在试样图环内表层铺塑料薄膜，向环内注入约5cm深的水层，打 开供水箱开关，用注射器抽水，直至马氏瓶能正常供水（目的是调节马氏瓶）。 (5) 检查秒表是否正常及回零位。 (6) 记录供水箱原始水位读数。 2. 实验方法及步骤 试验人员必须精力集中，认真负责，在统一指挥下，分工协作，作好记录。 a．迅速抽取塑料薄膜，并开始记时水位数值。 b．读取第一分钟末供水箱的水位，按试验要求读取水位数值。 c．实验至渗吸速度稳定后(即每两次水位读数差相同)，实验结束。 3. 注意事项 a．供水箱出水口必须淹没在内环水面以下0.5～1.0cm。 b．水位读数要读取每分钟末的数值，该数是计算渗吸规律重要的参数之一。 c．试验开始时迅速向外环加水至0.5～1.0cm时，使内外环水位大致保持相同水深，但外环加水不计入总量。 d．内环的供水量，由水箱上的标尺读数换算获取。 四、试验原理及资料分析整理
降雨和灌水入渗条件下土壤水分运动2.docx
第五章降雨和灌水入渗条件下土壤水分运 动 第一节水向土中入渗过程 一、概述 降雨和灌水入渗是田间水循环的重要环节，与潜水蒸发一样，是水资源评价和农田水分 状况调控的重要依据。 水渗入土壤的强度主要取决于降雨或灌水的方式和强度以及土壤渗水性能。如果土壤渗水性能较强，大于外界供水强度，则入渗强度主要决定于外界供水强度，在入渗过程中土壤表面含水率随入渗而逐渐提高，直至达到某一稳定值。如果降雨或灌水强度较大，超过了土壤渗水能力，入渗强度就决定于土壤的入渗性能，这样就会形成径流或地表积水。这两种情况可能发生在入渗过程的不同阶段，如在稳定灌溉强度（例如喷灌）下，开始时灌溉强度小于土壤入渗能力，入渗率等于灌溉强度；但经过一定时间后，土壤入渗能力减少，灌水强度大于土壤入渗能力，于是产生余水，如图2－5－ 1所示的降雨或灌水条件下的入渗过程。开始时入渗速率较高，以后逐渐减小。土壤的入渗能力随时间而变化，与土壤原始湿度和土壤 水的吸力有关，同时也与土壤剖面上土质条件、结构等因素有关。一般来说，开始入渗阶段，土壤入渗能力较高，尤其是在入渗初期，土壤 比较干燥的情况，然后随土壤水的入渗速率逐 渐减小，最后接近于一常量，而达到稳定入渗 阶段。 在较干旱的条件下，土壤表层的水势梯度 较陡。所以，入渗速率较大，但随着入渗水渗 入土中，土壤中基模吸力下降。湿润层的下移 使基模吸力梯度减小。在垂直入渗情况下，如 供水强度较大，使土壤剖面上达到饱和，当入 渗强度等于土壤饱和水力传导度时，将达到稳 定入渗阶段。如供水强度较小，小于饱和土壤 水力传导度时，达到稳定入渗阶段的入渗强度将等于该湿度条件下的非饱和土壤水力传导 度。 入渗过程中，土壤剖面上水分分布与土表入渗条件有关。根据 Coleman和Bodman 的研 究， 当均质土壤地表有积水入渗时，典型含水率分布剖面可分为四个区，即表层有一薄层为饱和带，以下是含水率变化较大的过渡带，其下是含水率分布较均匀的传导层，以下是湿润程度随深度减小的湿润层，该层湿度梯度越向下越陡，直到湿润锋。随着入渗时间延续，传导层 会不断向深层发展，湿润层和湿润锋也会下移，含水率分布曲线逐渐变平缓。
二维吸渗与入渗条件下土壤水力特性参数反演方法研究
二维吸渗与入渗条件下土壤水力特性参数反演方法研究 土壤水力特性参数取值是影响非饱和土壤水运动数值计算精度的关键。采用数值模拟、理论分析和室内试验对比相结合的技术路线,综合运用土壤水动力学、数值模拟与数值反演、多目标优化、代理模型和多种计算机语言综合集成技术,开展土壤二维负压吸渗、积水入渗水分运动参数的反演方法研究,取得以下主要结果:（1）提出了一种新的土壤水力特性参数反演方法,即“两步法”。第一步,以吸渗/入渗结束时刻的土壤含水率(θfinal),即ψ(θ final)最小作为目标函数,采用遗传算法反演饱和含水率;第二步, 以累积吸渗/入渗量ψ（Q）和吸渗/入渗速率ψ（v）最小作为目标函数,采用由多向量遗传算法和粒子群算法所构建的混合算法反演水力特性参数α、n和 Ks;与传统的加权和多目标反演方法相比,所提方法能够有效解决不同目标函数权重系数难以确定的问题,且具有高的求解效率和强的稳健性。（2）以所提“两步法”为基础,分别对二维吸渗和积水入渗条件下多种典型土壤、不同初始含水量条件下的van Genuchten–Mualem模型中水力特性参数进行了反演。 结果表明所得土壤水力特性参数反演值与典型土壤参考值（以RETC软件给出的典型值为比较时的参考值）具有好的一致性,说明所提反演方法具有高的可靠性;采用反演所得土壤水力特性参数分别绘制土壤水分特征曲线和导水率曲线,并与参考值绘制的曲线进行比较,结果表明两者具有高的一致性,说明反演所得 参数可较为精确的估算土壤水分特征曲线和导水率曲线;量化比较了考虑土壤含水率和累积入渗量存在测量误差条件下反演所得水力特性参数估算土壤水分特 征曲线和导水率曲线和参考值曲线,结果表明两者间具有小的差异和满意的估算精度,说明了所提反演方法具有强的稳健性;对积水入渗土壤垂直剖面含水率非 均一分布条件下水力特性参数进行了反演,结果表明典型土壤不同含水率分布模式下所得水力特性参数估算值与参考值差异较小,且采用反演结果绘制的土壤水分特征曲线和导水率曲线与参考值绘制的土壤水分特征曲线和导水率曲线基本 一致,说明所提反演范围具有较为广泛的使用范围,可用于生产实践。（3）建立了基于Kriging代理模型的土壤水力特性参数反演模型。根据土壤积水入渗的累积入渗量和最终含水率对土壤水力特性参数进行了反演估算,结果表明反演结果与典型土壤参考值具有高的一致性;量化比较了考虑土壤含水率和累积入渗量存在
第五章 降雨和灌水入渗条件下土壤水分运动2
第五章降雨和灌水入渗条件下土壤水分运动 第一节水向土中入渗过程 一、概述 降雨和灌水入渗是田间水循环的重要环节，与潜水蒸发一样，是水资源评价和农田水分状况调控的重要依据。 水渗入土壤的强度主要取决于降雨或灌水的方式和强度以及土壤渗水性能。如果土壤渗水性能较强，大于外界供水强度，则入渗强度主要决定于外界供水强度，在入渗过程中土壤表面含水率随入渗而逐渐提高，直至达到某一稳定值。如果降雨或灌水强度较大，超过了土壤渗水能力，入渗强度就决定于土壤的入渗性能，这样就会形成径流或地表积水。这两种情况可能发生在入渗过程的不同阶段，如在稳定灌溉强度（例如喷灌）下，开始时灌溉强度小于土壤入渗能力，入渗率等于灌溉强度；但经过一定时间后，土壤入渗能力减少，灌水强度大于土壤入渗能力，于是产生余水，如图2－5－1所示的降雨或灌水条件下的入渗过程。开始时入渗速率较高，以后逐渐减小。土壤的入渗能力随时间而变化，与土壤原始湿度和土壤水的吸力有关，同时也与土壤剖面上土质条件、结构等因素有关。一般来说，开始入渗阶段，土壤入渗能力较高，尤其是在入渗初期，土壤 比较干燥的情况，然后随土壤水的入渗速率逐 渐减小，最后接近于一常量，而达到稳定入渗 阶段。 在较干旱的条件下，土壤表层的水势梯度 较陡。所以，入渗速率较大，但随着入渗水渗 入土中，土壤中基模吸力下降。湿润层的下移 使基模吸力梯度减小。在垂直入渗情况下，如 供水强度较大，使土壤剖面上达到饱和，当入 渗强度等于土壤饱和水力传导度时，将达到稳 定入渗阶段。如供水强度较小，小于饱和土壤 水力传导度时，达到稳定入渗阶段的入渗强度将等于该湿度条件下的非饱和土壤水力传导度。 入渗过程中，土壤剖面上水分分布与土表入渗条件有关。根据Coleman和Bodman的研究，当均质土壤地表有积水入渗时，典型含水率分布剖面可分为四个区，即表层有一薄层为饱和带，以下是含水率变化较大的过渡带，其下是含水率分布较均匀的传导层，以下是湿润程度随深度减小的湿润层，该层湿度梯度越向下越陡，直到湿润锋。随着入渗时间延续，传导层会不断向深层发展，湿润层和湿润锋也会下移，含水率分布曲线逐渐变平缓。
土壤入渗实验报告
土壤入渗实验报告 一、实验目的： 进行土壤入渗试验，对土壤入渗规律有大致了解，并且利用测的数据绘出土壤累积入渗量和时间的关系曲线，利用该曲线求出入渗强度和时间的关系。 二、实验仪器： 直径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和时间 的关系。 四、实验数据： 五、数据分析：
01 2 3 4 5 6 5 10 15 20 25 30 35 4045入渗时间t/min 土壤累积入渗量I 和时间t 的关系曲线图 根据土壤累积入渗量I 和时间t 的拟合函数求导得到土壤入渗强度i 和时间t 的关系 00.05 0.10.150.20.250.30.350.40.45 10 20 30 40 50 入渗时间t/min 土壤入渗强度 i/cm·min 土壤入渗强度i 和时间t 的关系曲线图 六、 体会： 1. 土壤渗流试验看似简单，不过做成功很难。实验室中大多数仪器损坏。实验条件艰 难，但是还是尽量利用实验室中的其他条件勉强把实验做完了。
降雨和灌水入渗条件下土壤水分运动2
第五章降雨和灌水入渗条件下土壤水分运 动 第一节水向土中入渗过程 、概述 降雨和灌水入渗是田间水循环的重要环节， 与潜水蒸发一样，是水资源评价和农田水分 状况调控的重要依据。 水渗入土壤的强度主要取决于 降雨或灌水的方式 和强度以及土壤渗水性能。如果土壤渗 水性能较强，大于外界供水强度，则入渗强度主要决定于外界供水强度， 在入渗过程中土壤 表面含水率随入渗而逐渐提高， 直至达到某一稳定值。如果降雨或灌水强度较大，超过了土 壤渗水能力，入渗强度就决定于土壤的入渗性能， 这样就会形成径流或地表积水。 这两种情 况可能发生在入渗过程的不同阶段， 如在稳定灌溉强度（例如喷灌）下， 开始时灌溉强度小 于土壤入渗能力，入渗率等于灌溉强度； 但经过一定时间后，土壤入渗能力减少，灌水强度 大于土壤入渗能力，于是产生余水，如图 2-5- 1所示的降雨或灌水条件下的入渗过程。 开 始时入渗速率较高， 以后逐渐减小。土壤的入渗能力随时间而变化， 与土壤原始湿度和土壤 水的吸力有关，同时也与土壤剖面上土质条件、 土壤入渗能力较高，尤其是在入渗初期，土壤 比较干燥的情况，然后随土壤水的入渗速率逐 渐减小，最后接近于一常量，而达到稳定入渗 阶段。 在较干旱的条件下，土壤表层的水势梯度 较陡。所以，入渗速率较大，但随着入渗水渗 入土中，土壤中基模吸力下降。湿润层的下移 使基模吸力梯度减小。在垂直入渗情况下，如 供水强度较大，使土壤剖面上达到饱和，当入 渗强度等于土壤饱和水力传导度时，将达到稳 定入渗阶段。如供水强度较小，小于饱和土壤 水力传导度时，达到稳定入渗阶段的入渗强度将等于该湿度条件下的非饱和土壤水力传导 度。 入渗过程中，土壤剖面上水分分布与土表入渗条件有关。 根据Coleman 和Bodman 的研究, 当均质土壤地表有积水入渗时， 典型含水率分布剖面可分为四个区， 即表层有一薄层为饱和 带，以下是含水率变化较大的过渡带， 其下是含水率分布较均匀的传导层， 以下是湿润程度 随深度减小的湿润层，该层湿度梯度越向下越陡， 直到湿润锋。随着入渗时间延续，传导层 会不断向深层发展，湿润层和湿润锋也会下移，含水率分布曲线逐渐变平缓。 结构等因素有关。一般来说，开始入渗阶段, 图2-5-1入滲率对时间的茏祭曲线
第三章 蒸发条件下土壤水分运动2
第三章 蒸发条件下土壤水分运动 土壤水分蒸发可以发生在土壤表面和植物体上。植物体的蒸发一般称为蒸腾，土壤表面蒸发称为土面蒸发。本章主要讨论土面蒸发条件下土壤水分运动。 土面蒸发所消耗的水分来自两部分，一部分是指直接消耗地下水面以上土层中水分，一部分消耗地下水，消耗地下水部分称为潜水蒸发。土壤水分蒸发有稳定蒸发和不稳定蒸发两种状态，当土壤水分的蒸发量与地下水补给量相平衡时为土壤水分的稳定蒸发状态，一般在连续干旱期，且地下水有侧向补给时，会出现这种情况；当土壤水蒸发量不等于地下水补给量时，土壤水分为不稳定蒸发状态，在降雨或灌水后的蒸发初期或地下水无侧向补给时，常处于这种状态。 土壤水蒸发一方面决定于外界（大气）蒸发能力（常以水面蒸发表示），另一方面决定于土层从地下水面向地表输水的能力，其输水能力大小一方面取决于土质条件，同时也决定于表土含水率。由于土壤水蒸发的水分是从土表散失，因此为了研究土壤水蒸发问题首先必须了解表土蒸发规律。 第一节 表土蒸发 一、形成干土层前的表土蒸发过程 由于表土蒸发的主要影响因素不同，其蒸发过程可以分为以下两个阶段。 1．表土蒸发保持稳定阶段 表层土壤水分的蒸发主要是由于土壤水汽压力与地表大气中水汽压力有一定差值，在压力梯度作用下，土壤中水汽向大气中扩散而产生的。压力差越大，土壤中水汽扩散的水量越大。在这一阶段，表土含水率越高（在某一定值以上），土壤水汽压力基本不随含水率的变化而改变，其数值趋近于饱和水汽压力。在这种情况下，土壤水分蒸发主要取决于外界条件（温度、湿度、风速等）。在外界条件不变的情况下，土壤水分的蒸发将不随含水率降低而变化，因此这一阶段称为稳定蒸发阶段，蒸发强度可用下式表示： )(0101P P -=βε （2-3-1） 式中：1ε––––稳定蒸发阶段土壤水分蒸发强度（m 3/d ·m 2 或 m/d ）； 0β––––质量交换系数，与外界条件有关； 1P ––––土壤表层的水汽压力（Pa ）； 0P ––––大气中的水汽压力（Pa ）。 在这一阶段内土壤水蒸发强度接近水面蒸发强度，而与土壤含水率无关。此阶段土壤含 水率的下限（临界含水率）即是蒸发强度与土壤毛管输水能力保持平衡之点；临界含水率即
土壤入渗实验报告
河海大学 土壤入渗实验报告班级：农水2班
土壤入渗实验报告 一、实验目的： 进行土壤入渗试验，对土壤入渗规律有大致了解，并且利用测的数据绘出土壤累积入渗量和时间的关系曲线，利用该曲线求出入渗强度和时间的关系。 二、实验仪器： 直径4.5cm的土柱圆筒、宽5cm长4cm的马氏瓶、托盘、烧杯、天平、直尺、滤纸、秒表、疏松土壤、自来水。 三、实验步骤： 1)把准备好的土装入圆筒中，每装5cm就夯实一次，直至土的的顶端位于圆筒 进水孔的下缘。 2)检测马氏瓶是否漏气。如果漏气就换试验设备或者用凡士林涂抹。 3)将水灌入马氏瓶，把马氏瓶的出水孔和圆筒的进水口用橡胶管连接好，调节 圆 4)筒和马氏瓶的相对高度使得马氏瓶的出水孔刚好出水。 5)读取马氏瓶中水的高度。 6)实验开始，量取土柱量筒中土壤稳定下渗时土柱的淹水的深度，分别读取实 验 7)开始后第1、3、5、7、10、15、20、25、30、35、40min钟时马氏瓶中水的 高度。 8)根据马氏瓶的面积和圆筒的面积求出圆筒中从开始到不同时刻的累积入渗 量I，画出土壤累积入渗量I和时间的关系曲线，并利用该曲线求出入渗强度i和时间的关系。 四、实验数据： 见后附页。 五、数据分析： I f=0.20 cm/min ； S=2.44 ； I1=1.22 cm/min ； ?=0.51 。
入渗强度-时间曲线图 入渗强度-时间双对数曲线图 湿润锋深度与时间关系图时间t/min 湿 润 锋 深 度 /mm F=38.5×t0.563
随着入渗时间的延长，土壤累积的入渗深度和入渗水量不断增大，初始时刻增加较快，随着时间的不断延长，土壤累积的入渗深度和入渗水量增长速度减缓。 土壤入渗强度随时间的延长而迅速减小，在初始时刻减小很快，随着时间的推移最后趋近于某一稳定值，形成这种现象的主要因素是入渗路径的不断加长，从水柱面到入渗锋面的水势梯度逐渐减小所以入渗强度也在不断减小，最后接近于该种土壤的渗透系数。 六、注意事项： 1.土壤入渗要现实土壤足够湿润，然后才能开始记录试验数据，否则会出错。 2.土壤入渗强度开始较强，后来变小，逐渐稳定。 3.土壤渗流试验看似简单，不过做成功很难，有许多需要值得注意的地方。 4.读数时秒表不能停，从打开止水阀门时按下秒表，直到入渗完全。 5.土壤入渗要现实土壤足够湿润，然后才能开始记录试验数据，否则会出错。 6.装土时不能把土压的太紧也不能不压，必须保持每次加入土的量和夯实的次 数一致，尽量保证土壤入渗的均匀度。 7.土壤入渗强度开始较强，后来变小，逐渐稳定，需要大致了解这个规律 8.往马氏瓶注水时需要打开阀门，注完后需要关闭才能开始实验。 9.
数学物理方程课程
《数学物理方程》课程 教学大纲 课程代码：B0110040 课程名称：数学物理方程／equation of mathematic physics 课程类型：学科基础课 学时学分：64学时/4学分 适用专业：地球物理学 开课部门：基础课教学部 一、课程的地位、目的和任务 课程的地位：数学物理方程是地球物理学专业的一门重要的专业（或技术）基础课。数学物理方程是反应自然中物理现象的基本模型，也是一种基本的数学工具，与数学其他学科和其他科学技术领域诸如数值分析、优化理论、系统工程、物理、化学、生物等学科都有广泛联系。对于将来从事工程地震技术工作及自然科学研究的学生来说是必不可少的。期望学生通过该门课程的学习，能深刻地理解数学物理方程的不同定解问题所反应的物理背景。 课程的目的与任务：使学生了解数学物理方程建立的依据和过程，认识这门学科与物理学、力学、化学、生物学等自然科学和社会科学以及工程技术的极密切的广泛的联系。掌握经典数学物理方程基本定解问题的提法和相关的基本概念和原理，重点掌握求解基本线性偏微分方程定解问题的方法和技巧。使学生掌握与本课程相关的重要理论的同时，注意启发和训练学生联系自己的专业，应用所学知识来处理和解决实际问题的能力。 二、课程与相关课程的联系与分工 学生在进入本课程学习之前，应修课程包括：大学物理、高等数学、线性代数、复变函数、场论与向量代数。这些课程的学习，为本课程奠定了良好的数学基础。本课程学习结束后，可进入下列课程的学习：四大力学、电磁场与微波技术、近代物理实验等。且为进一步选修偏微分方程理论、数值计算、控制理论与几何分析等课程打下基础。
三、教学内容与基本要求 第一章绪论 1.教学内容 第一节偏微分方程的基本概念 第二节弦振动方程及定解条件 第三节热传导方程及定解条件 第四节拉普拉斯方程及定解条件 第五节二阶线性偏微分方程的分类 第六节线性算子 2.重点难点 重点：物理规律“翻译”成数学物理方程的思路和步骤，实际问题近似于抽象为理想问题 难点：数学物理方程的数学模型建立及数学物理方程的解空间是无限维的函数空间 3.基本要求 （1）了解数学物理方程研究的基本内容，偏微分方程的解、阶、维数、线性与非线性、齐次与非齐次的概念；了解算子的定义。了解三类典型方程的建立及其定解问题（初值问题、边值问题和混合问题）的提法，定解条件的物理意义。 （2）掌握微分算子的运算规律，理解线性问题的叠加原理 （3）了解二阶线性方程的特征理论 （4）掌握两个变量二阶线性偏微分方程分类方法及化简方法 （5）掌握三类方程的标准形式及其化简过程，会三类方程的比较，并能通过标准形式求得某些方程的通解。 第二章分离变量法 1.教学内容 第一节有界弦的自由振动。 第二节有界长杆的热传导问题。 第三节二维拉普拉斯方程的边值问题。 第四节非齐次方程得求解问题。