基于van Genuchten-Mualem模型的饱和-非饱和介质流动随机数值分析
- 格式:pdf
- 大小:916.22 KB
- 文档页数:7
不同区域土壤水下渗在四水转化中作用研究[摘要] 根据水循环机理及水平衡机理,现有水存在形式主要为大气水、地表水、土壤水、地下水,四种水存在形式的转化具有复杂的联系。
其中土壤水在转化过程中占据中心位置。
土壤水是指地面以下至潜水面以上土壤层中的水分,是联系地表水、地下水的纽带,同时也是物质传输和运移的载体。
本实验基于灌溉下渗实验模型,主要利用时域反射仪(TDR)连续测量土壤剖面水分含量的方法,分别观测北京通州区土壤和昆明土壤在降水入渗与再分配过程。
根据采集的数据经过图像处理与HYDRUS-1D软件模拟图像进行对比得出不同种土壤水入渗过程中不同深度初始含水率以及下渗率、剖面含水率等对实验的影响,据此得出相关结论。
[关键字]土壤下渗灌溉下渗模型时域反射仪初始含水率剖面含水率1 概述“四水”转化是一个庞大的系统工程,对于其中的每一方面而言,都是一个相对独立的子系统,每个子系统各自独立、自成体系,但它们又相互关联、相互依存、相互转化和相互制约,形成了一个独特的循环转化系统。
2 试验条件与方法研究区北京市通州区通州区位于北京市东南部,通州区地处永定河、潮白河冲积洪积平原,地势平坦,自西北向东南倾斜。
其土质多为潮黄土、两合土、沙壤土,土壤肥沃,质地适中。
云南省昆明市位于北纬亚热带,素以”春城”之称而享誉中外。
四是冬无严寒,日照充足,天睛少雨。
两地从地形、气候、降水、土壤等多方面具有不同特点,差异较大,以两地为研究对象进行对比研究实验结果将相对明显。
本实验的土柱灌溉下渗实验模型是由北京土和昆明土的两个土柱以及时域反射仪组成。
从试验地选取优质土进行分装,实验装置如图1所示。
土柱装填在有机玻璃柱中,有机玻璃柱高度65.0cm,内径28.1cm,在不同的深度剖面上水平埋设时域反射仪TDR)探针,每10分钟自动测定一次土壤水分含水率。
时域反射仪的基本原理是,高频电磁脉冲沿传输线在土壤中传播的速度依赖于土壤的介电特性,在一定的电磁波频率范围内(50M~10GHz),矿物质、空气和水的介电特性为常数,因此土体的介电常数主要依赖于土壤容积含水量,这样可以建立土壤容积含水量与土壤介电常数的经验方程,TDR通过测量高频电磁脉冲在土壤中的传播速度求得土壤的介电常数从而计算出土壤的含水量[2]。
中国农业大学科技成果——饱和-非饱和介质中水分和溶质不规则迁移机理与模拟成果简介饱和-非饱和介质中水分与溶质迁移转化理论,是农业水资源可持续利用与水土环境保护的重要理论基础。
本项目主要以认识含水介质中水分与溶质不规则迁移规律,发展相应的模拟理论与方法为目标,以室内外试验、理论分析和数值模拟为手段,系统开展了非饱和土壤水分动态随机与分形模拟理论、含水介质中溶质非费克迁移理论及模拟模型、含水层中抽水井附近非达西流理论与模型、区域地下水动态模拟的理论与方法等方面的研究。
技术内容(1)提出了土壤水力特性参数求解的新方法,建立了非饱和土壤水分动态随机模拟的新理论;(2)探讨了介质非均质性对溶质非费克迁移的影响机理,建立了基于分数微分理论的吸附性溶质迁移的模拟模型,发展了模型求解的新方法;(3)研究了非达西流动的机理,建立了抽水条件下非达西流动的模拟模型,提出了抽水井水力学计算的新方法;(4)结合我国北方地区的实际情况,建立了地下水动态模拟的人工神经网络(ANN)与FEFLOW耦合模型。
该项目解决了含水介质中水分与溶质不规则迁移的一些关键性理论问题,提高了土壤与含水层中水分与溶质迁移模拟与预测的精度和可靠性,揭示了饱和与非饱和介质中水分与溶质不规则迁移的内在机理,深化了水分与溶质迁移不规则现象的认识,丰富了农业资源与环境、农业水利学科的内容,对相关学科的发展有重要促进和推动作用。
该项目在非饱和水分随机运动、非达西流动、非费克迁移的微观机理和宏观规律的研究方面具有明显的创新性,国际著名水文学家Fred Molz教授等认为我们的有关成果达到了国际先进水平。
研究成果在农业水资源高效利用和水土环境保护等领域中具有潜在的应用价值。
投资预算及效益分析研究成果在农业水资源高效利用和水土环境保护等领域中具有潜在的应用价值。
*国家重点基础研究发展规划项目(G1999011803)和ACIAR 项目(LWR1/96/164)资助收稿日期:2000-12-28;收到修改稿日期:2002-03-24封丘地区土壤传递函数的研究*朱安宁 张佳宝(中国科学院南京土壤研究所,南京 210008) 陈效民 陈德立(南京农业大学资源与环境学院,南京 210095) (澳大利亚墨尔本大学)摘 要 以黄淮海平原封丘地区的潮土和风沙土为研究对象,试图寻求解决土壤水动力学特性参数问题的更加实用可行的方法。
根据大量的土壤基本物理性质和土壤持水数据,利用多元逐步回归分析方法,分别建立了van Genuchten 模型的参数H r 、H s 、a 、n 的土壤传递函数模型,并通过统计分析和数学模拟进行模型验证。
研究结果表明,van Genuchten 模型参数与土壤的质地、有机质含量及容重等基本物理性质有一定的线性相关性。
van Genuchten 模型的参数估计模型对粘性土壤的拟合效果较好,而拟合砂性土壤时误差较大。
关键词 土壤物理性质,水分特征曲线,van Genuchten 模型,土壤传递函数,数学模拟中图分类号 S152172当前,区域土壤水分和溶质运移的研究已成为土壤物理学研究的一大热点,然而,用计算机数值模拟土壤水和溶质运移所遇到的主要问题之一是难于获取土壤水分运动参数,如土壤水分特征曲线h (H )、饱和导水率K s 、非饱和导水率K (H )等。
Childs 早在1940年就注意到土壤物理性质影响土壤水分特性,特别是70年代末以来,许多学者在研究土壤水分运动参数与土壤基本理化性质关系方面做了大量的工作[1,2],从而建立了研究土壤水分运动参数的一种新方法)))PTFs 法(Pedo -Transfer Functions,土壤传递函数)。
PTFs 法是用统计模型由土壤基本性质预测土壤水分运动参数的方法[3,4],即由已知的土壤水分运动参数和与之相对应的土壤基本理化性质建立多元回归方程,用此回归方程结合土壤基本理化性质来预测土壤的水分运动参数。
基于HYDRU S模型的盐碱地土壤水盐运移模拟潘延鑫;罗纨;贾忠华;井思媛;李山;武迪【摘要】To determine the movement of salt and water in the saline-alkali flatland of Lubotan,Shaanxi province, based on the saturated-unsaturated soil water and solute transport theory,field monitoring data of local water and salt for many years was applied to simulate the rules of local soil water and salt movement,using the HYDRUS 1D numerical model .The soil water and salt changes were analyzed,and the reasonable field irrigation quota was proposed .The results showed that during the whole reproductive period,soil water content had a similar variation trend under different irrigation quotas .Considering water saving and salt control,farmland irrigation quota of 500 m3·hm-2 is reasonable to control salt accumulation in soil .The simulated results of soil water and salt migration using HYDRUS model are basically consistent with the observed values in field experiment,and the results can be referred for farmland management of water and salt in semblable salinity regions .%为了解陕西卤泊滩盐碱地的水盐运移情况,基于当地2009—2013年田间水盐监测资料,应用饱和-非饱和土壤水分及溶质运移理论,利用HYDRUS-1 D数值模型对当地土壤水分、盐分运移规律进行数值模拟,分析了盐碱地的水盐变化状况,确定合理的田间灌水定额。
非饱和土壤导水率试验计算与模拟分析胡钜鑫;虎胆·吐马尔白;穆丽德尔·托伙加;杨未静【摘要】以非饱和土壤导水率作为研究对象,用瞬时剖面法计算两种土壤非饱和土壤导水率,并与RETC中不同模型的模拟结果进行对比,研究瞬时剖面法计算结果的可靠性.结果表明:两种土壤的K-h与lgK-h模拟曲线和实测值均吻合较好,实测值和不同模型的模拟值均属于高度性相关,且K-θ实测曲线与各模型的模拟曲线变化规律相似,处于各模拟曲线之间.综上所述,瞬时剖面法计算结果与模拟结果相似,具有一定的准确性,可以直接使用在实际生产运用过程中.【期刊名称】《石河子大学学报(自然科学版)》【年(卷),期】2019(037)001【总页数】7页(P105-111)【关键词】非饱和土壤导水率;瞬时剖面法;van Genuchten模型;Mualem模型【作者】胡钜鑫;虎胆·吐马尔白;穆丽德尔·托伙加;杨未静【作者单位】新疆农业大学水利与土木工程学院,新疆乌鲁木齐市,830052;水文水资源与水利工程科学国家重点实验室,江苏南京,210098;新疆农业大学水利与土木工程学院,新疆乌鲁木齐市,830052;新疆农业大学水利与土木工程学院,新疆乌鲁木齐市,830052;新疆农业大学水利与土木工程学院,新疆乌鲁木齐市,830052【正文语种】中文【中图分类】S152.7非饱和土壤土导水率K 是土壤水分参数中的重要参数之一,它反⒊了土壤中的水分在非饱和状态下的运动规律。
非饱和土壤导水率的测定方法包括直接法和间接法,直接法又分为田间测定和室内测定。
田间测定方法包括结壳法[1]、圆盘入渗法[2-4]、双环法[5]等,室内测定方法包括瞬时剖面法、垂直下渗通量法、零通量法[6]等。
其中直接测量法通常耗时耗力,不易测量,因此大部分学者常选⒚间接方法求取非饱和导水率,包括土壤水分再分布法[7-8],或者通过水分特征曲线C 和水平扩散度D 公式推求非饱和土壤导水率K[9],另外通过模拟软件[10],例如Hydrus 和RETC 通过土壤质地资料推求非饱和导水率[11-13]。
自由渗流面及渗流量推求的数值计算方法探讨秦甜甜;丁国辉;程勤波【摘要】自由渗流面具有复杂的非线性,较难确定.本文采用开源地下水数值模拟程序MODFLOW及SUTRA,分别运用MODFLOW模型中干湿单元转化技术、SUTRA模型中单元渗透矩阵调整法以及本文建立的缓变渗透系数矩阵法推求自由渗流面.对比其求解结果表明,采用MODFLOW运用干湿转化技术求解自由渗透面的方法稳定性最好、精度最高,而采用缓变渗透系数矩阵法的SUTRA程序,改善了传统单元渗透矩阵调整法的不稳定性,提高了数值计算精度,避免了MODFLOW必须矩形网格的局限性,是一种实用的计算自由渗流面及估算地下水与河流水量交换量的方法.【期刊名称】《地下水》【年(卷),期】2018(040)004【总页数】3页(P12-14)【关键词】MODFLOW;SUTRA;单元渗透矩阵调整法;缓变渗透系数矩阵法;自由渗流面【作者】秦甜甜;丁国辉;程勤波【作者单位】江苏省地质调查研究院,江苏南京 210018;南京市测绘勘察研究院股份有限公司,江苏南京 210005;河海大学水文水资源学院/河海大学水文水资源与水利工程科学国家重点实验室,江苏南京 210098【正文语种】中文【中图分类】P641.2自由渗流面与地下水潜水水面特征相似,具有复杂的非线性,是河流等地表水入渗补给地下水、大坝渗流等定量化计算的难点,目前,尚未发现确定自由渗流面的解析法,在实际工程应用中常用数值方法求解。
求解该类问题的数值方法可分为两类:变网格法和固定网格法。
由于变网格法操作复杂,容易使计算网格畸变,难以处理存在水平介质分层及各种复杂夹层的情况等缺点,因而研究较少。
为此,国内外学者提出了许多与固定网格法相关的处理技术,主要有剩余流量法、初流量法、节点虚流量法、单元渗透矩阵调整法、复合单元法及截止负压法等[1-6]。
但这些处理技术大多方法复杂,操作困难,应用不便。
本文利用开源地下水数值模拟软件MODFLOW[7]及SUTRA[8],提出了采用干湿单元转化、缓变渗透系数矩阵法推求自由渗流面的方法,成功推求了自由渗流面,有效地提高了目前河流与地下水交互作用计算的求解精度。
不同出流方式下垂直流人工湿地数值模拟范思雨;龙天渝【摘要】建立了二级垂直流人工湿地的二维非恒定微生物动力学模型,针对淹没出流方式、不同的等间隔淹没出流方式和定水头出流方式,采用数值模拟的方法研究二级垂直流人工湿地处理重庆市径流雨水的效果;结果表明,在相同工况下,间隔较大的淹没出流方式能达到最大的有效体积比和水力效率,水流在三层填料中与各种类微生物接触更加充分、持久,在此出流方式下促进了溶解氧的传质,复氧后溶解氧浓度回升更快、浓度峰值更大、高浓度维持时间更长,因此促进了由微生物主导的生物化学降解反应,二级垂直流湿地对COD、氮和磷的去除效果也优于其他出流方式.【期刊名称】《重庆工商大学学报(自然科学版)》【年(卷),期】2017(034)006【总页数】9页(P110-118)【关键词】人工湿地;出流方式;CW2D;数值模拟【作者】范思雨;龙天渝【作者单位】重庆大学三峡库区生态环境教育部重点实验室,重庆,400045;重庆大学三峡库区生态环境教育部重点实验室,重庆,400045;重庆大学低碳绿色建筑国际联合研究中心,重庆,400045【正文语种】中文【中图分类】X703人工湿地污水处理技术因其显著的低能耗、低维护成本、处理范围广和良好的景观效应而被广泛应用[1],又因其较强的抗冲击负荷能力而在蓄纳、处理城市径流雨水上有很好的应用前景[2-6]。
目前大多数人工湿地的设计与运行依然建立在经验的基础上,对其全运行过程中复杂的物理-化学-生物协同作用认知的“黑箱”现象依然存在,而这一作用与人工湿地内的水流规律密切相关[5]。
在模拟研究中将湿地的出流方式抽象成边界条件,边界条件的设定在流体力学的计算中极大地影响了水流规律,因此对人工湿地出流方式的探讨和研究对提高其处理效率、指导工程设计具有重要意义。
目前对于人工湿地进/出水方式的研究主要是两个方面:一是讨论穿孔管在进水处和出水处的配/集水的均匀性;二是讨论进/出水口在垂直方向上的位置(如上、中、下)及其组合搭配对配水均匀性和处理效果的影响。
求van Genuchten模型参数的AM-MCMC方法石晓蕾;徐绍辉;廖凯华【期刊名称】《土壤》【年(卷),期】2012(44)2【摘要】采用基于自适应采样算法的马尔科夫链蒙特卡罗方法(简称AM-MCMC)来估算描述土壤水分特征曲线的van Genuchten模型的参数,并推求出模型参数的后验分布,从而为模型参数的不确定性分析提供依据.结果表明,对于van Genuchten 模型而言,采用AM-MCMC算法能得到模型参数后验均值和方差的分布,并且能推求出模型参数的置信区间,所以用这种算法来求解van Genuchten方程的参数是行之有效的,为求解van Genuchten模型参数提供了一种新的思路.%This paper adopted the adaptive Metropolis algorithm of Markov chain Monte Carlo approach (AM-MCMC) to estimate the parameters of van Genuchten model, which describes the soil water characteristic curve. The posterior statistics of the van Genuchten model parameters were obtained to estimated the model uncertainty. It was found that the distribution of the posterior means and variances of the model parameters could be obtained. Confidence intervals of the model parameters could be ascertained. It was proved effective and novel in solving van Genuchten equation with this algorithm.【总页数】6页(P345-350)【作者】石晓蕾;徐绍辉;廖凯华【作者单位】青岛大学环境科学系,青岛266071;青岛大学环境科学系,青岛266071;南京大学水科学系,南京210093【正文语种】中文【中图分类】S152【相关文献】1.van Genuchten模型参数的物理意义 [J], 陈卫金;程东会;陶伟2.土壤持水曲线van Genuchten模型求参的一种新方法 [J], 刘贤赵;李嘉竹;张振华3.土壤持水曲线van Genuchten模型求参的Matlab实现 [J], 魏义长;刘作新;康玲玲;王云璋;时明立4.利用优化方法求算Van Genuchten方程参数 [J], 李春友;任理;李保国5.土壤水分特征曲线Van-Genuchten模型参数的土壤传输函数比选 [J], 李彬楠;樊贵盛因版权原因,仅展示原文概要,查看原文内容请购买。