CN111898083A - 一种描述干团聚体粒径分布的模型 - Google Patents
一种描述干团聚体粒径分布的模型 Download PDFInfo
- Publication number
- CN111898083A CN111898083A CN202010790137.7A CN202010790137A CN111898083A CN 111898083 A CN111898083 A CN 111898083A CN 202010790137 A CN202010790137 A CN 202010790137A CN 111898083 A CN111898083 A CN 111898083A
- Authority
- CN
- China
- Prior art keywords
- model
- dry
- particle size
- size distribution
- aggregates
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 60
- 239000002245 particle Substances 0.000 title claims abstract description 54
- 239000002689 soil Substances 0.000 claims abstract description 14
- 238000009825 accumulation Methods 0.000 claims abstract description 10
- 238000012216 screening Methods 0.000 claims abstract description 7
- 230000003628 erosive effect Effects 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract description 3
- 230000001186 cumulative effect Effects 0.000 description 5
- 238000000034 method Methods 0.000 description 5
- 230000015556 catabolic process Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000007873 sieving Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
技术领域
本发明属于土壤流失防治技术领域,具体涉及一种描述干团聚体粒径分布的模型。
背景技术
风蚀是导致干旱和半干旱地区肥沃表层土壤颗粒物质流失的主要自然驱动力。土壤风蚀会降低土地生产力,导致土壤退化。土壤表层粗造影响因子和干团聚体粒径分布是影响风蚀的关键因素。干团聚体粒径是影响风蚀模型的重要输入参数,为了提高风蚀模型计算的精确度,需要定量描述土壤干团聚体粒径的分布。平筛(或套筛)是获取干团聚体粒径分布的传统方法,为了减小手动筛分造成的误差,多使用振筛机和一套平筛获取干团聚体粒径分布。为了减小振筛机对干团聚体的破坏,Chepil设计了旋筛来获取干团聚体粒径的分布。干团聚体粒径分布的旋筛方法被广泛应用于WEQ、RWEQ、WEPS等风蚀模型中。
几十年来,科学家一直在探索干团聚体粒径分布参数化的方法。使用数学函数拟合来确定干团聚体粒径分布被认为是一种有效的方法。Marshall和Quirk建立了WD模型(Weibull distribution)来确定干团聚体粒径分布,并通过实验验证了WD模型描述的干团聚体破碎的过程,所述方程为:
其中d为筛孔,P(d)为干团聚体小于筛孔d的累积频率,α为尺度因子,β为形状因子。
Perfect等人介绍了Mandelbrot提出的FD模型(Fractal distribution)来表达干团聚体粒径的分布。实验结果表明,干团聚体的质量大小和数量大小分布均可用FD模型表示,所述方程为:
其中dmax为最大筛孔尺寸,D为分形维数。
Wagner and Ding(1993)进一步改进Lognormal distribution(LD)模型,改进后的模型MLN描述干团聚体粒径分布方程为:
其中d0拟合参数,μ是几何平均粒径(geometric median diameter,GMD),δ是干团聚体粒径分布的标准差(log of the geometric standard deviation,logGSD)。
Perfect等人1993年在不同实验地收集了干团聚体粒径分布数据,在比较FD模型(Fractal distribution)和WD模型(Weibull distribution)后发现,FD模型(Fractaldistribution)表现出了更高的精度。Zobeck等人2003年利用在美国六个州的24个地点采集的5400多个表层土壤样品验证了FD模型和WD模型,结果表明WD模型表现更好。同时发现不同尺寸的筛孔对模型的拟合结果有很大影响,在实践中,不同的实验目的、使用不同的操作程序对干团聚体进行筛分,干团聚体的粒径分布会呈现不同的曲线模式。此外,当前对数正态模型(Lognormal distribution,LD)、分形模型(Fractal distribution,FD)和威布尔模型(Weibull distribution,WD)存在不能很好地表达右侧单峰值和多峰值的干团聚体粒径分布情况的问题,针对此情况,本发明提出了一种结合幂函数和指数函数的新模型来表达干团聚体粒径分布,并且利用几个国家不同出版物的数据进行了验证。
发明内容
本发明的目的就是提供一种描述干团聚体粒径分布的模型,以解决现有干团聚体粒径分布模型不能很好地模拟不同形态的干团聚体粒径分布的问题。
本发明的目的是通过以下技术方案实现的:一种描述干团聚体粒径分布的模型,所述模型结合幂函数和指数函数来表达干团聚体的粒径分布,模型的方程式为:
其中,d为筛孔粒径,P(d)为干团聚体累积物小于筛孔粒径的可能性,λ、γ为拟合参数,dr为干团聚体累积频率等于最大干团聚体累积频率的0.63时的直径,其中,取1-e-1≈0.63。
所述模型可描述不同土壤质地、土地利用方式和筛分方式的干团聚体粒径分布。
较大尺寸的干团聚体粒径分布可由幂函数表达,较小尺寸的干团聚体粒径分布可由幂函数和指数函数联合表达。
本发明针对当前对数正态模型(Lognormal distribution,LD)、分形模型(Fractal distribution,FD)和威布尔模型(Weibull distribution,WD)不能很好的表达右侧单峰值和多峰值的干团聚体粒径分布的情况,通过联合幂函数和指数函数得到新模型后发现,较大尺寸的干团聚体粒径分布可由幂函数表达,较小尺寸的干团聚体粒径分布可由幂函数和指数函数联合表达。
总体上,利用不同土壤质地、不同土地利用方式和不同筛分方式的土壤干团聚体粒径分布数据对新模型和已有模型进行对比分析,新模型模拟精度最高,为提高风蚀模型计算的精确度提供了指导。在使用五个国家的253份数据与其他模型比较后发现,新模型的性能更好,且较大尺寸的干团聚体粒径分布可由幂函数表达,较小尺寸的干团聚体粒径分布可由幂函数和指数函数联合表达。筛分过程中筛孔的选择会影响干团聚体粒径分布的形态,并进一步影响新模型的精度。
附图说明
图1是幂函数和指数函数对新模型影响示意图。其中,d为筛孔,P(d)为干团聚体小于筛孔d的累积频率。
图2是多模态数据和单峰向右偏的干团聚体粒径分布数据图。其中,图(2)和(4)分别以对数形式与图(1)和(3)相对应,d为筛孔,dmax为每个土样的最大筛孔粒径,P(d)为干团聚体小于筛孔d的累积频率。
具体实施方式
实施例1
在对比其他模型后,结合幂函数和指数函数得出新模型,所述方程如下:
其中,d为筛孔粒径,P(d)为干团聚体累积物小于筛孔粒径的可能性,λ,γ为拟合参数,dr(mm)为干团聚体累积频率等于最大干团聚体累积频率的0.63(1-e-1≈0.63)时的直径。
在应用时,通过对比不同类型的数据发现,较大尺寸的干团聚体粒径分布呈现幂函数,较小尺寸的干团聚体粒径分布呈现为幂函数和指数函数联合表达的形式。
实施例2
使用不同国家的253份数据(如表1所示)对FD模型、WD模型、MLN模型和新模型进行比较,用调整的相关系数aR2和均方根误差RMSE来验证三个模型的性能,结果见表2所示,其中aR2的计算公式为:
其中,R2为相关系数,n为干团聚体粒径分布的筛孔数,k为回归模型相关系数的数量。
表1:不同国家干团聚体收集数据
表2:不同干团聚体粒径分布模型的统计特征值和回归参数
表2中Max为最大值,Min为最小值,Ave为平均值,STD为标准差。
从表2中可以看出,FD、WD、MLN和新模型aR2的平均值分别为0.7680、0.9192、0.9082和0.9495,RMSE的平均值分别为0.1067、0.0532、0.0562和0.0350。对于FD、WD和MLN来说,aR2值小于0.9的有145个(72.50%)、73个(31.20%)和105个(41.50%),而对于新模型,aR2值小于0.9的有27个(10.71%)。对于FD、WD和MLN模型,RMSE的值大于0.1的有106个(53.00%)、15个(6.41%)和26个(10.28%),而对于新模型,RMSE的值大于0.1的只有7个(2.79%)。新模型的aR2和RMSE的标准差也很小。这些结果表明,新模型性能比其他模型更好。
实施例3
筛孔数的多少会影响干团聚体粒径分布和模型的精度。对于同一土壤样品,16个筛孔(0.01、0.03、0.045、0.063、0.09、0.125、0.18、0.25、0.355、0.5、0.71、1、1.4、2、4和12.5毫米)的干团聚体粒径分布是多种形态的(图1a),16个筛孔中的7个(0.125、0.25、0.5、1、2、4和12.5毫米)干团聚体粒径分布是左侧单峰的(图1b)。我们进一步总结了FD、WD、MLD和新模型的统计特性,如表3所示,FD、WD、MLD的模型精度显著提高,各模型的不同参数也发生了变化。用16个和7个筛孔的观测数据与新模型的模拟结果比较相似,这表明新模型能够用不同的形态描述干团聚体的粒径分布。且在图2和表4中可以看出FD模型和WD模型在描述多模态数据和单峰向右偏数据时性能比新模型要差。
表3:不同干团聚体粒径分布模型的统计特征值和回归参数
表3中Max为最大值,Min为最小值,Ave为平均值,STD为标准差。
表4:不同干团聚体粒径分布模型的统计特征值和回归参数
表4中Max为最大值,Min为最小值,Ave为平均值,STD为标准差。
Claims (3)
2.根据权利要求1所述的描述干团聚体粒径分布的模型,其特征在于,所述模型可描述不同土壤质地、土地利用方式和筛分方式的干团聚体粒径分布。
3.根据权利要求1所述的描述干团聚体粒径分布的模型,其特征在于,较大尺寸的干团聚体粒径分布可由幂函数表达,较小尺寸的干团聚体粒径分布可由幂函数和指数函数联合表达。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010790137.7A CN111898083B (zh) | 2020-08-07 | 2020-08-07 | 一种描述干团聚体粒径分布的模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010790137.7A CN111898083B (zh) | 2020-08-07 | 2020-08-07 | 一种描述干团聚体粒径分布的模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111898083A true CN111898083A (zh) | 2020-11-06 |
CN111898083B CN111898083B (zh) | 2024-03-22 |
Family
ID=73246711
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010790137.7A Active CN111898083B (zh) | 2020-08-07 | 2020-08-07 | 一种描述干团聚体粒径分布的模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111898083B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0990888A2 (en) * | 1998-09-29 | 2000-04-05 | Horiba, Ltd. | Apparatus and method for measuring a particle size distribution |
US20080208511A1 (en) * | 2005-03-07 | 2008-08-28 | Michael Trainer | Methods and apparatus for determining characteristics of particles |
JP2009036533A (ja) * | 2007-07-31 | 2009-02-19 | Kajima Corp | 地盤材料の粒度計測システム及びプログラム |
JP2012242099A (ja) * | 2011-05-16 | 2012-12-10 | Kajima Corp | 粒状材料の分割式粒度計測方法及びシステム |
JP2013257188A (ja) * | 2012-06-12 | 2013-12-26 | Kajima Corp | 粒状材料の粒度分布計測方法及びシステム |
-
2020
- 2020-08-07 CN CN202010790137.7A patent/CN111898083B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0990888A2 (en) * | 1998-09-29 | 2000-04-05 | Horiba, Ltd. | Apparatus and method for measuring a particle size distribution |
US20080208511A1 (en) * | 2005-03-07 | 2008-08-28 | Michael Trainer | Methods and apparatus for determining characteristics of particles |
US20080204716A1 (en) * | 2005-03-07 | 2008-08-28 | Michael Trainer | Methods and apparatus for determining characteristics of particles |
JP2009036533A (ja) * | 2007-07-31 | 2009-02-19 | Kajima Corp | 地盤材料の粒度計測システム及びプログラム |
JP2012242099A (ja) * | 2011-05-16 | 2012-12-10 | Kajima Corp | 粒状材料の分割式粒度計測方法及びシステム |
JP2013257188A (ja) * | 2012-06-12 | 2013-12-26 | Kajima Corp | 粒状材料の粒度分布計測方法及びシステム |
Non-Patent Citations (2)
Title |
---|
ZHONGLING GUO ET AL.: "Logistic growth models for describing the fetch effect of aeolian sand transport", 《SOIL & TILLAGE RESEARCH 194 (2019) 104306》, pages 1 - 9 * |
黄亚鹏: "不同时间尺度农田风沙流模拟", 《中国沙漠》, pages 1 - 8 * |
Also Published As
Publication number | Publication date |
---|---|
CN111898083B (zh) | 2024-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109815223A (zh) | 一种针对工业监测数据缺失的补全方法及补全装置 | |
CN105445160B (zh) | 一种沥青混合料的空隙特征及其提取方法 | |
CN1132265C (zh) | 用于自动调节介质滤波器的特性的方法和设备 | |
CN110349637B (zh) | 骨料周围界面过渡区体积分数预测方法、装置及终端设备 | |
JP2008523470A5 (zh) | ||
CN111898083A (zh) | 一种描述干团聚体粒径分布的模型 | |
CN109596440A (zh) | 一种沥青混合料温度疲劳性能的评价方法 | |
CN107220946B (zh) | 一种岩石运输带上不良块度图像的实时剔除方法 | |
CN108111852A (zh) | 面向量化分块压缩感知编码的双测量参数率失真控制方法 | |
CN109191571B (zh) | 一种应用3d打印技术制备力学试验标准集料的方法 | |
CN109507971A (zh) | 一种机制砂质量智能监控系统及监控方法 | |
CN116069672B (zh) | 操作系统内核定向模糊测试的种子变异方法及测试方法 | |
CN112557244A (zh) | 用于评价沥青混合料中玄武岩纤维分散均匀性的方法 | |
Ahiakwo et al. | Effects of sieve aperture modification on dewatered cassava mash sieving process | |
CN110362857B (zh) | 岩石结构面各向异性分布特征的理论模型解析方法 | |
CN108996932A (zh) | 一种沥青路面回收料的厂拌热再生级配设计方法 | |
CN108288266A (zh) | 一种混凝土再生骨料边角损失率和砂浆剥落面积的量化评价方法 | |
CN114611298A (zh) | 一种通过喷墨3d打印砂型铺砂过程建模与仿真优化工艺参数的方法和设备 | |
CN108195644A (zh) | 利用石墨烯提高土的无侧限抗压强度的方法 | |
CN1641637A (zh) | 联机分析处理系统及方法 | |
CN103761525A (zh) | 基于显微ct数字图像处理铝电解用阴极炭块组分检测方法 | |
CN114092400B (zh) | 一种基于体视学理论利用二维图像提取沥青混合料细观特征的方法 | |
CN111797730A (zh) | 一种水泥熟料岩相自动分析方法 | |
CN205701369U (zh) | 一种筛分筒 | |
CN111488649A (zh) | 一种联合收割机零部件载荷分布的非参数估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |