CN105809680A - 一种基于多重分形的土壤水分特征曲线模拟方法 - Google Patents
一种基于多重分形的土壤水分特征曲线模拟方法 Download PDFInfo
- Publication number
- CN105809680A CN105809680A CN201610125382.XA CN201610125382A CN105809680A CN 105809680 A CN105809680 A CN 105809680A CN 201610125382 A CN201610125382 A CN 201610125382A CN 105809680 A CN105809680 A CN 105809680A
- Authority
- CN
- China
- Prior art keywords
- alpha
- soil
- max
- formula
- fractal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Abstract
本发明属于农业技术领域,具体说,涉及一种基于多重分形的土壤水分特征曲线模拟方法,包括:获取土壤数字图像;对所述土壤数字图像进行处理,获取土壤孔隙分布特征数据,并计算土壤孔隙分布的多重分形参数;根据土壤孔隙直径与土壤毛管吸力方程确定土壤水分特征曲线中压力水头与土壤孔隙半径之间的关系;建立压力水头与土壤孔隙半径多重分形参数的奇异性指数函数,并对所述函数进行积分变换,获得土壤水分特征曲线模型。本发明基于多重分形方法,充分考虑考虑土壤结构对水分特征曲线模拟的影响,提高了土壤水分特征曲线的模拟精度。
Description
技术领域
本发明属于农业技术领域,具体说,涉及一种基于多重分形的土壤水分特征曲线模拟方法。
背景技术
土壤水分特征曲线是描述土壤含水量与吸力(基质势)之间的关系曲线。它反映了土壤水能量与土壤水含量的函数关系,是进行土壤水分和溶质运移模拟的重要物理参数,土壤水分特征曲线的模拟精度直接影响了土壤水分和溶质运移模拟预测的精度,因此它是表示土壤基本水力特性最重要的指标之一。
土壤水分特征曲线是一个高度非线性的函数,受到多种环境和土壤本身因素的制约,在空间上和时间上都具有强烈的变异性,准确测定非常困难。目前确定土壤水分特征曲线的方法主要有直接测定法和间接模拟法两种,直接方法确定土壤水分特征曲线的过程费时耗力,且测定结果具有较大不确定性;间接方法由于较少的考虑土壤结构对水分特征曲线模拟的影响,导致土壤水分模拟的精度不高。
发明内容
本发明的目的在于基于多重分形理论,提供一种间接确定土壤水分特征曲线的方法,以提高土壤水分特征曲线的模拟精度。
本发明的实施例提供了一种基于多重分形的土壤水分特征曲线模拟方法,包括:
步骤1,获取土壤数字图像;
步骤2,对该土壤数字图像进行处理,获取土壤孔隙分布特征数据,并计算土壤孔隙分布的多重分形参数;
步骤3,根据土壤孔隙直径与土壤毛管吸力Young-Laplace方程确定土壤水分特征曲线中压力水头与土壤孔隙半径之间的关系;
步骤4,建立压力水头与土壤孔隙半径多重分形参数的奇异性指数函数,并对该函数进行积分变换,获得土壤水分特征曲线模型。
进一步,步骤3具体包括:
假设土壤的孔隙分布分数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)即为基于多重分形的土壤水分特征曲线模型与现有技术相比本发明的有益效果是:基于多重分形方法,充分考虑考虑土壤结构对水分特征曲线模拟的影响,提高了土壤水分特征曲线的模拟精度。
附图说明
图1是本发明一种基于多重分形的土壤水分特征曲线模拟方法的流程图。
具体实施方式
下面结合附图所示的各实施方式对本发明进行详细说明,但应当说明的是,这些实施方式并非对本发明的限制,本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代,均属于本发明的保护范围之内。
参图1所示,图1是本发明一种基于多重分形的土壤水分特征曲线模拟方法的流程图。
本实施例提供了一种基于多重分形的土壤水分特征曲线模拟方法,包括:
步骤S1,获取土壤数字图像;
步骤S2,对土壤数字图像进行处理,获取土壤孔隙分布特征数据,并计算土壤孔隙分布的多重分形参数;
步骤S3,根据土壤孔隙直径与土壤毛管吸力Young-Laplace方程确定土壤水分特征曲线中压力水头与土壤孔隙半径之间的关系;
步骤S4,建立压力水头与土壤孔隙半径多重分形参数的奇异性指数函数(将土壤水分特征曲线中的压力水头表示为土壤孔隙多重分形参数的函数),并对该函数进行积分变换,获得(基于土壤孔隙分布多重分形奇异性指数的)土壤水分特征曲线模型。
本实施例提供的土壤水分特征曲线模拟方法,基于多重分形方法,充分考虑考虑土壤结构对水分特征曲线模拟的影响,提高了土壤水分特征曲线的模拟精度。
下面详述本发明的具体实施过程:
1)选定研究区域,田间采集土壤样本,通过CT图像扫描获得土壤数字图像;采用计算机数字图像处理技术和ImageJ图像处理软件获得土壤孔隙分布状况,并计算土壤孔隙分布的多重分形参数。
2)假设土壤的孔隙分布分数Pi包含在孔隙当中可以用孔隙半径ri来估计:
式中c1为经验常数,αi为奇异性指数,孔隙半径ri、奇异性指数为αi的孔隙数目Ni之间满足下式:
式中c2为经验常数,f(αi)为具有奇异性指数αi的孔隙分形维数,我们认为在所研究的多孔介质的所有孔隙中经验参数c1和c2相等,Pi的值只随αi值变化,因此包含在所有孔隙当中的部分孔隙可以表示为:
对(3)式积分有:
式中αmax和αmin分别为奇异性指数α的最大值和最小值。
3)确定α和f(α)之间的关系,用下面的公式来描述他们之间的关系:
式中α0=(αmax+αmin)/2。
4)孔隙半径与压力水头之间的关系可以用毛管上升理论方程表示,根据土壤毛管吸力Young-Laplace方程:
式中:h为土壤负压(基质势);r为孔隙尺寸;σ为表面张力;ω为接触角。在25℃时,水的表面张力σ约为0.075N/m,一般取接触角ω为0°。则:
式中h是压力水头。
将式(5)、(7)代入式(4)可得:
5)相应于某压力水头h的土壤含水量θ(h)可以通过对式(8)的数值积分得到,即:
式中:hmin是相应于最小孔隙的压力水头,a1和a2可表示为:
6)假定土壤孔隙半径r达到最大孔隙半径rmax时的土壤含水量θ(r)等于饱和含水量θs,相应于最大孔隙半径rmax的压力水头hmax为进气口压力,饱和含水量可以用下式估计:
由式(9)和(11)可以得到由参数θs,αmax,αmin,hmin和hmax表示的土壤水分特征曲线模型:
a1和a2可表示为:
式(12)即为基于多重分形的土壤水分特征曲线模型。
本实施例提供的土壤水分特征曲线模拟方法,采用多重分形方法定量描述土壤结构的非均质特征,并应用于土壤水分特征曲线的模拟中,提高了土壤水分特征曲线的模拟精度。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
Claims (3)
1.一种基于多重分形的土壤水分特征曲线模拟方法,其特征在于,包括:
步骤1,获取土壤数字图像;
步骤2,对所述土壤数字图像进行处理,获取土壤孔隙分布特征数据,并计算土壤孔隙分布的多重分形参数;
步骤3,根据土壤孔隙直径与土壤毛管吸力Young-Laplace方程确定土壤水分特征曲线中压力水头与土壤孔隙半径之间的关系;
步骤4,建立压力水头与土壤孔隙半径多重分形参数的奇异性指数函数,并对所述函数进行积分变换,获得土壤水分特征曲线模型。
2.根据权利要求1所述的一种基于多重分形的土壤水分特征曲线模拟方法,其特征在于,所述步骤3具体包括:
假设土壤的孔隙分布分数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是压力水头。
3.根据权利要求2所述的一种基于多重分形的土壤水分特征曲线模拟方法,其特征在于,所述步骤4具体包括:
将式(5)、(7)代入式(4):得:
通过对式(8)的数值积分得到相应于某压力水头h的土壤含水量θ(h),即:
式中:hmin是相应于最小孔隙的压力水头,a1和a2表示为:
假定土壤孔隙半径r达到最大孔隙半径rmax时的土壤含水量θ(r)等于饱和含水量θs,相应于最大孔隙半径rmax的压力水头hmax为进气口压力,饱和含水量用下式估计:
由式(9)和(11)可以得到由参数θs,αmax,αmin,hmin和hmax表示的土壤水分特征曲线模型:
式(12)即为基于多重分形的土壤水分特征曲线模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610125382.XA CN105809680B (zh) | 2016-03-06 | 2016-03-06 | 一种基于多重分形的土壤水分特征曲线模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610125382.XA CN105809680B (zh) | 2016-03-06 | 2016-03-06 | 一种基于多重分形的土壤水分特征曲线模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105809680A true CN105809680A (zh) | 2016-07-27 |
CN105809680B CN105809680B (zh) | 2019-02-26 |
Family
ID=56466699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610125382.XA Active CN105809680B (zh) | 2016-03-06 | 2016-03-06 | 一种基于多重分形的土壤水分特征曲线模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105809680B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106290796A (zh) * | 2016-09-13 | 2017-01-04 | 中国水利水电科学研究院 | 一种基于多重分形的土壤盐分空间变异性表征方法 |
CN106769629A (zh) * | 2017-01-04 | 2017-05-31 | 塔里木大学 | 一种土壤质量含水量检测方法 |
CN108444895A (zh) * | 2018-06-14 | 2018-08-24 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
CN109211756A (zh) * | 2018-10-31 | 2019-01-15 | 湖北工业大学 | 基于vg模型的变形土进气值预测方法 |
CN109711103A (zh) * | 2019-01-28 | 2019-05-03 | 中国科学院地理科学与资源研究所 | 全范围土壤水分特征曲线拟合方法及拟合系统 |
CN110261583A (zh) * | 2019-07-23 | 2019-09-20 | 石河子大学 | 精确测量低持水能力多孔介质水分特征曲线的实验方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104182647A (zh) * | 2014-09-03 | 2014-12-03 | 塔里木大学 | 获得土壤水分特征曲线Van Genchten模型参数的方法 |
CN204439487U (zh) * | 2015-03-13 | 2015-07-01 | 西北农林科技大学 | 一种利用离心机法测量土壤水分特征曲线的系统 |
JP2015176485A (ja) * | 2014-03-17 | 2015-10-05 | 国立大学法人豊橋技術科学大学 | 三次元モデル特徴抽出方法、及び三次元モデル・アノテーションシステム |
-
2016
- 2016-03-06 CN CN201610125382.XA patent/CN105809680B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015176485A (ja) * | 2014-03-17 | 2015-10-05 | 国立大学法人豊橋技術科学大学 | 三次元モデル特徴抽出方法、及び三次元モデル・アノテーションシステム |
CN104182647A (zh) * | 2014-09-03 | 2014-12-03 | 塔里木大学 | 获得土壤水分特征曲线Van Genchten模型参数的方法 |
CN204439487U (zh) * | 2015-03-13 | 2015-07-01 | 西北农林科技大学 | 一种利用离心机法测量土壤水分特征曲线的系统 |
Non-Patent Citations (2)
Title |
---|
刘建立 等: "估计土壤水分特征曲线的间接方法研究进展", 《水利学报》 * |
黄冠华 等: "土壤水分特性曲线的分形模拟", 《水科学进展》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106290796A (zh) * | 2016-09-13 | 2017-01-04 | 中国水利水电科学研究院 | 一种基于多重分形的土壤盐分空间变异性表征方法 |
CN106769629A (zh) * | 2017-01-04 | 2017-05-31 | 塔里木大学 | 一种土壤质量含水量检测方法 |
CN106769629B (zh) * | 2017-01-04 | 2019-02-19 | 塔里木大学 | 一种土壤质量含水量检测方法 |
CN108444895A (zh) * | 2018-06-14 | 2018-08-24 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
CN108444895B (zh) * | 2018-06-14 | 2020-08-25 | 长安大学 | 一种高效黄土体非饱和渗透参数获取方法 |
CN109211756A (zh) * | 2018-10-31 | 2019-01-15 | 湖北工业大学 | 基于vg模型的变形土进气值预测方法 |
CN109211756B (zh) * | 2018-10-31 | 2020-11-03 | 湖北工业大学 | 基于vg模型的变形土进气值预测方法 |
CN109711103A (zh) * | 2019-01-28 | 2019-05-03 | 中国科学院地理科学与资源研究所 | 全范围土壤水分特征曲线拟合方法及拟合系统 |
CN110261583A (zh) * | 2019-07-23 | 2019-09-20 | 石河子大学 | 精确测量低持水能力多孔介质水分特征曲线的实验方法 |
CN110261583B (zh) * | 2019-07-23 | 2022-03-04 | 石河子大学 | 精确测量低持水能力多孔介质水分特征曲线的实验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105809680B (zh) | 2019-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105809680A (zh) | 一种基于多重分形的土壤水分特征曲线模拟方法 | |
Tiwari et al. | A diffuse interface model with immiscibility preservation | |
CN103198494B (zh) | 一种基于建筑外轮廓线及激光点云的建筑立面提取方法 | |
BRPI0503982A (pt) | método, sistema, e meio fìsico legìvel por computador para criar e desenhar um gráfico dentro de um programa de aplicação | |
CN104954185B (zh) | 一种基于深度置信网的云计算负载预测方法 | |
CN108090262B (zh) | 一种基于sph理论直升机着水载荷计算的粒子对构建方法 | |
JP2012164314A5 (zh) | ||
CN103390274A (zh) | 一种基于区域相关信息熵的图像分割质量评价方法 | |
CN108242962B (zh) | 一种基于测量报告的室内信号传播损耗计算方法及装置 | |
EP2650806A3 (en) | Simulation method for a tire | |
CN106202002A (zh) | 一种用于检测水文系列参数是否变异的方法 | |
CN105045958A (zh) | 基于bp神经网络的gps高程拟合的实现系统及方法 | |
CN104298737A (zh) | 一种基于抽样最大核密度稳健模型的离散点云拟合方法 | |
RU2017140845A (ru) | Сегментация анатомической структуры на основе модели | |
Nasonova et al. | Climate change impact on streamflow in large-scale river basins: projections and their uncertainties sourced from GCMs and RCP scenarios | |
WO2017201813A1 (zh) | 地面灌溉地表水流运动过程模拟方法 | |
SG10201809536XA (en) | Systems and methods for quantification of postural balance of users in an augmented reality environment | |
CN107292100B (zh) | 一种岩石结构面分形维数的计算方法 | |
CN103268423A (zh) | 一种地理现象多点模拟空间尺度选择的方法 | |
CN108280265A (zh) | 基于新型径向基函数的网格变形插值方法、系统及介质 | |
CN106815607A (zh) | 一种基于反距离权重插值反函数的等值线图像数据提取方法 | |
CN103246636B (zh) | 一种汉字字形视觉重心的确定方法 | |
CN107506572B (zh) | 获取目标点的高度的方法和装置 | |
CN103142212B (zh) | 基于一致球的身体对称性检测方法 | |
CN104572575A (zh) | 一种特大变形动态网格生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |