首页 / 专利库 / 土壤科学 / 土壤学 / 土壤学 / 土壤特性 / 土壤水分张力 / 一种基于多重分形的土壤水分特征曲线模拟方法

一种基于多重分形的土壤分特征曲线模拟方法

阅读:457发布:2020-05-15

专利汇可以提供一种基于多重分形的土壤分特征曲线模拟方法专利检索,专利查询,专利分析的服务。并且本 发明 属于农业技术领域,具体说,涉及一种基于多重分形的 土壤 水 分特征曲线模拟方法,包括:获取土壤数字图像;对所述土壤数字图像进行处理,获取土壤孔隙分布特征数据,并计算土壤孔隙分布的多重分形参数;根据土壤孔隙直径与土壤毛管吸 力 方程确定土壤水分特征曲线中压力 水头 与土壤孔隙半径之间的关系;建立压力水头与土壤孔隙半径多重分形参数的奇异性指数函数,并对所述函数进行积分变换,获得土壤水分特征曲线模型。本发明基于多重分形方法,充分考虑考虑土壤结构对水分特征曲线模拟的影响,提高了土壤水分特征曲线的模拟 精度 。,下面是一种基于多重分形的土壤分特征曲线模拟方法专利的具体信息内容。

1.一种基于多重分形的土壤分特征曲线模拟方法,其特征在于,包括:
步骤1,获取土壤数字图像;
步骤2,对所述土壤数字图像进行处理,获取土壤孔隙分布特征数据,并计算土壤孔隙分布的多重分形参数;
步骤3,根据土壤孔隙直径与土壤毛管吸Young-Laplace方程确定土壤水分特征曲线中压力水头与土壤孔隙半径之间的关系,具体包括:
假设土壤的孔隙分布分数Pi包含在孔隙当中,并用孔隙半径ri进行估计:
式中c1为经验常数,αi为奇异性指数;
孔隙半径ri、奇异性指数为αi的孔隙数目Ni之间满足下式:
式中c2为经验常数,f(αi)为具有奇异性指数αi的孔隙分形维数,设定多孔介质的所有孔隙中经验常数c1和c2相等,Pi的值只随αi值变化,则包含在所有孔隙当中的部分孔隙可以表示为:
对(3)式积分有:
式中αmax和αmin分别为奇异性指数α的最大值和最小值;
确定α和f(α)之间的关系,用下面的公式来描述他们之间的关系:
式中α0=(αmax+αmin)/2;
根据土壤毛管吸力Young-Laplace方程获取土壤孔隙半径与压力水头之间的关系,所述土壤毛管吸力Young-Laplace方程为:
式中:h为压力水头;r为土壤孔隙半径;σ为表面张力;ω为接触,在25℃时,水的表面张力σ取0.075N/m,接触角ω取0°,则:
式中h是压力水头;
步骤4,建立压力水头与土壤孔隙半径多重分形参数的奇异性指数函数,并对所述函数进行积分变换,获得土壤水分特征曲线模型,具体包括:
将式(5)、(7)代入式(4):得:
通过对式(8)的数值积分得到相应于某压力水头h的土壤含水量θ(h),即:
式中:hmin是相应于最小孔隙的压力水头,a1和a2表示为:
假定土壤孔隙半径r达到最大孔隙半径rmax时的土壤含水量θ(r)等于饱和含水量θs,相应于最大孔隙半径rmax的压力水头hmax为进气口压力,饱和含水量用下式估计:
由式(9)和(11)可以得到由参数θs,αmax,αmin,hmin和hmax表示的土壤水分特征曲线模型:
式(12)即为基于多重分形的土壤水分特征曲线模型。

说明书全文

一种基于多重分形的土壤分特征曲线模拟方法

技术领域

[0001] 本发明属于农业技术领域,具体说,涉及一种基于多重分形的土壤水分特征曲线模拟方法。

背景技术

[0002] 土壤水分特征曲线是描述土壤含水量与吸(基质势)之间的关系曲线。它反映了土壤水能量与土壤水含量的函数关系,是进行土壤水分和溶质运移模拟的重要物理参数,土壤水分特征曲线的模拟精度直接影响了土壤水分和溶质运移模拟预测的精度,因此它是表示土壤基本水力特性最重要的指标之一。
[0003] 土壤水分特征曲线是一个高度非线性的函数,受到多种环境和土壤本身因素的制约,在空间上和时间上都具有强烈的变异性,准确测定非常困难。目前确定土壤水分特征曲线的方法主要有直接测定法和间接模拟法两种,直接方法确定土壤水分特征曲线的过程费时耗力,且测定结果具有较大不确定性;间接方法由于较少的考虑土壤结构对水分特征曲线模拟的影响,导致土壤水分模拟的精度不高。

发明内容

[0004] 本发明的目的在于基于多重分形理论,提供一种间接确定土壤水分特征曲线的方法,以提高土壤水分特征曲线的模拟精度。
[0005] 本发明的实施例提供了一种基于多重分形的土壤水分特征曲线模拟方法,包括:
[0006] 步骤1,获取土壤数字图像;
[0007] 步骤2,对该土壤数字图像进行处理,获取土壤孔隙分布特征数据,并计算土壤孔隙分布的多重分形参数;
[0008] 步骤3,根据土壤孔隙直径与土壤毛管吸力Young-Laplace方程确定土壤水分特征曲线中压力水头与土壤孔隙半径之间的关系;
[0009] 步骤4,建立压力水头与土壤孔隙半径多重分形参数的奇异性指数函数,并对该函数进行积分变换,获得土壤水分特征曲线模型。
[0010] 进一步,步骤3具体包括:
[0011] 假设土壤的孔隙分布分数Pi包含在孔隙当中,并用孔隙半径ri进行估计:
[0012]
[0013] 式中c1为经验常数,αi为奇异性指数;
[0014] 孔隙半径ri、奇异性指数为αi的孔隙数目Ni之间满足下式:
[0015]
[0016] 式中c2为经验常数,f(αi)为具有奇异性指数αi的孔隙分形维数,设定多孔介质的所有孔隙中经验常数c1和c2相等,Pi的值只随αi值变化,则包含在所有孔隙当中的部分孔隙可以表示为:
[0017]
[0018] 对(3)式积分有:
[0019]
[0020] 式中αmax和αmin分别为奇异性指数α的最大值和最小值;
[0021] 确定α和f(α)之间的关系,用下面的公式来描述他们之间的关系:
[0022]
[0023] 式中α0=(αmax+αmin)/2;
[0024] 据土壤毛管吸力Young-Laplace方程获取土壤孔隙半径与压力水头之间的关系,土壤毛管吸力Young-Laplace方程为:
[0025]
[0026] 式中:h为压力水头;r为土壤孔隙半径;σ为表面张力;ω为接触,在25℃时,水的表面张力σ取0.075N/m,接触角ω取0°,则:
[0027]
[0028] 式中h是压力水头。
[0029] 进一步,步骤4具体包括:
[0030] 将式(5)、(7)代入式(4):得:
[0031]
[0032] 通过对式(8)的数值积分得到相应于某压力水头h的土壤含水量θ(h),即:
[0033]
[0034] 式中:hmin是相应于最小孔隙的压力水头,a1和a2表示为:
[0035]
[0036] 假定土壤孔隙半径r达到最大孔隙半径rmax时的土壤含水量θ(r)等于饱和含水量θs,相应于最大孔隙半径rmax的压力水头hmax为进气口压力,饱和含水量用下式估计:
[0037]
[0038] 由式(9)和(11)可以得到由参数θs,αmax,αmin,hmin和hmax表示的土壤水分特征曲线模型:
[0039]
[0040] 式(12)即为基于多重分形的土壤水分特征曲线模型与现有技术相比本发明的有益效果是:基于多重分形方法,充分考虑考虑土壤结构对水分特征曲线模拟的影响,提高了土壤水分特征曲线的模拟精度。附图说明
[0041] 图1是本发明一种基于多重分形的土壤水分特征曲线模拟方法的流程图

具体实施方式

[0042] 下面结合附图所示的各实施方式对本发明进行详细说明,但应当说明的是,这些实施方式并非对本发明的限制,本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代,均属于本发明的保护范围之内。
[0043] 参图1所示,图1是本发明一种基于多重分形的土壤水分特征曲线模拟方法的流程图。
[0044] 本实施例提供了一种基于多重分形的土壤水分特征曲线模拟方法,包括:
[0045] 步骤S1,获取土壤数字图像;
[0046] 步骤S2,对土壤数字图像进行处理,获取土壤孔隙分布特征数据,并计算土壤孔隙分布的多重分形参数;
[0047] 步骤S3,根据土壤孔隙直径与土壤毛管吸力Young-Laplace方程确定土壤水分特征曲线中压力水头与土壤孔隙半径之间的关系;
[0048] 步骤S4,建立压力水头与土壤孔隙半径多重分形参数的奇异性指数函数(将土壤水分特征曲线中的压力水头表示为土壤孔隙多重分形参数的函数),并对该函数进行积分变换,获得(基于土壤孔隙分布多重分形奇异性指数的)土壤水分特征曲线模型。
[0049] 本实施例提供的土壤水分特征曲线模拟方法,基于多重分形方法,充分考虑考虑土壤结构对水分特征曲线模拟的影响,提高了土壤水分特征曲线的模拟精度。
[0050] 下面详述本发明的具体实施过程:
[0051] 1)选定研究区域,田间采集土壤样本,通过CT图像扫描获得土壤数字图像;采用计算机数字图像处理技术和ImageJ图像处理软件获得土壤孔隙分布状况,并计算土壤孔隙分布的多重分形参数。
[0052] 2)假设土壤的孔隙分布分数Pi包含在孔隙当中可以用孔隙半径ri来估计:
[0053]
[0054] 式中c1为经验常数,αi为奇异性指数,孔隙半径ri、奇异性指数为αi的孔隙数目Ni之间满足下式:
[0055]
[0056] 式中c2为经验常数,f(αi)为具有奇异性指数αi的孔隙分形维数,我们认为在所研究的多孔介质的所有孔隙中经验参数c1和c2相等,Pi的值只随αi值变化,因此包含在所有孔隙当中的部分孔隙可以表示为:
[0057]
[0058] 对(3)式积分有:
[0059]
[0060] 式中αmax和αmin分别为奇异性指数α的最大值和最小值。
[0061] 3)确定α和f(α)之间的关系,用下面的公式来描述他们之间的关系:
[0062]
[0063] 式中α0=(αmax+αmin)/2。
[0064] 4)孔隙半径与压力水头之间的关系可以用毛管上升理论方程表示,根据土壤毛管吸力Young-Laplace方程:
[0065]
[0066] 式中:h为土壤负压(基质势);r为孔隙尺寸;σ为表面张力;ω为接触角。在25℃时,水的表面张力σ约为0.075N/m,一般取接触角ω为0°。则:
[0067]
[0068] 式中h是压力水头。
[0069] 将式(5)、(7)代入式(4)可得:
[0070]
[0071] 5)相应于某压力水头h的土壤含水量θ(h)可以通过对式(8)的数值积分得到,即:
[0072]
[0073] 式中:hmin是相应于最小孔隙的压力水头,a1和a2可表示为:
[0074]
[0075] 6)假定土壤孔隙半径r达到最大孔隙半径rmax时的土壤含水量θ(r)等于饱和含水量θs,相应于最大孔隙半径rmax的压力水头hmax为进气口压力,饱和含水量可以用下式估计:
[0076]
[0077] 由式(9)和(11)可以得到由参数θs,αmax,αmin,hmin和hmax表示的土壤水分特征曲线模型:
[0078]
[0079] a1和a2可表示为:
[0080]
[0081]
[0082] 式(12)即为基于多重分形的土壤水分特征曲线模型。
[0083] 本实施例提供的土壤水分特征曲线模拟方法,采用多重分形方法定量描述土壤结构的非均质特征,并应用于土壤水分特征曲线的模拟中,提高了土壤水分特征曲线的模拟精度。
[0084] 上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
[0085] 对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈