CN110286144A - 一种粮堆结露预测方法 - Google Patents
一种粮堆结露预测方法 Download PDFInfo
- Publication number
- CN110286144A CN110286144A CN201910626913.7A CN201910626913A CN110286144A CN 110286144 A CN110286144 A CN 110286144A CN 201910626913 A CN201910626913 A CN 201910626913A CN 110286144 A CN110286144 A CN 110286144A
- Authority
- CN
- China
- Prior art keywords
- temperature
- moisture condensation
- grain
- point
- line
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/56—Investigating or analyzing materials by the use of thermal means by investigating moisture content
- G01N25/66—Investigating or analyzing materials by the use of thermal means by investigating moisture content by investigating dew-point
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种粮堆结露预测方法,包括:步骤一、选取粮堆内竖直方向的测温平面,并且获取测温平面的粮温数组;步骤二、计算露点温度,将粮温数组中每个点位的温度分别与露点温度进行对比,并记录结露点的坐标;步骤三、制作测温平面的温度场云图,并将结露点的坐标在云图中依次连接形成结露线;步骤四、计算结露线的长度L及结露线的移动速度V,并根据结露线的长度及结露线的移动速度判断粮堆是否有结露风险。本发明提供的粮堆结露预测方法,仅需要粮堆内的粮温数据以及粮食的水分数据即可预测出粮堆是否有结露风险,需要的粮情数据种类较少,便于实施,结露判断结果生成快速。
Description
技术领域
本发明属于粮情监控技术领域,特别涉及一种粮堆结露预测方法。
背景技术
粮堆结露是指粮堆中的水汽达到饱和状态,开始在粮粒、器材、仓墙表面出现凝结水的现象。粮堆结露是威胁粮食安全的重要隐患之一。当粮堆内温度骤降到一定程度,籽粒间孔隙中所含水汽达到饱和状态时,粮食籽粒表面开始结露。
粮堆结露发生是一个逐步累积的过程,与粮堆温湿度长之间的耦合密切相关。中国专利,申请号“201410222195.4”根据粮食平衡水分原理和WU模型,通过对粮堆各个测量平面虚拟测量点温湿度的重建,根据温湿度梯度、水汽分压压力梯度实现对粮堆结露位置与实践的预测,避免发生结露和霉变。但该方法判断结露过程中,判断步骤较为复杂,需计算粮堆内的温度场、湿度场以及水蒸气分压等因素,综合判断。据了解,现今仍有大部分粮情监控系统仍没有装备湿度检测传感器,粮情数据中缺失粮堆的湿度数据。因此目前上述结露判断预警方法难以实施。
而据经验知,正常储粮过程中,粮堆内粮食水分一般变化幅度不大,可以用入仓时的水分替代粮堆内粮食籽粒的水分。由此可见,可以获得粮堆内较为准确的温度与水分数据。因此,需要一种能够依据粮堆温度和水分等变量对粮堆是否结露进行准确、快速的判断的新方法,该新方法能够依据粮堆温度、水分等变量实现对结露位置的快速判断,避免因粮情数据中缺失湿度数据而无法判断的问题,确保储粮安全。
发明内容
本发明的目的是提供一种粮堆结露预测方法,以历史储藏过程中或实时的粮情数据为依据,根据粮堆温度与水分,实现对粮堆结露时间与结露位置的快速判断,简化了粮堆结露判断过程和方法。
本发明提供的技术方案为:
一种粮堆结露预测方法,包括:
步骤一、选取粮堆内竖直方向的测温平面,并且获取所述测温平面的粮温数组后制作温度场云图;
步骤二、计算露点温度,将所述粮温数组中每个点位的温度分别与所述露点温度进行对比,并记录结露点的坐标;
其中,所述结露点为温度低于所述露点温度的点位;以及
所述露点温度为:td=ERHr×(A+B×t)+C×t-19.2;
式中,t为干球温度;ERHr为粮食平衡相对湿度,A、B、C分别为公式参数;
步骤三、将所述结露点的坐标在所述温度场云图中依次连接形成结露线;
步骤四、计算所述结露线的长度及结露线的移动速度,
当L≥L0,并且V≥V0时,判断所述测温平面有结露风险;
其中,L为所述结露线的长度,L0为结露线临界长度,V为所述结露线的移动速度,V0为结露线临界移动速度。
优选的是,在所述步骤二中,粮食平衡相对温度为:
式中,ERHr为粮食平衡相对湿度,M为粮食含水率,T为粮食温度,A1、A2、B1、B2和D分别为公式参数。
优选的是,在所述步骤四中,所述结露线的长度为:
l=C×R;
式中,R为温度场云图中结露线所占的像素点个数,C为每个像素点对应的实际长度。
优选的是,所述每个像素点对应的实际长度为:
式中,N为所述测温平面实际长度,n为所述测温平面的温度云场图的长度所占的像素点个数。
优选的是,所述结露线的移动速度为:
式中,Δd为采样周期,相邻两个周期的所有结露点在温度云场的高度方向上的平均移动距离所占的像素点个数,G为每个像素点对应的实际高度。
优选的是,所述相邻两个采样周期的所有结露点在温度云场的高度方向上的平均移动距离所占的像素点个数为:
式中,p为结露点的个数,Si为相邻采样周期中的同一个结露点移动的距离。
本发明的有益效果是:
本发明提供的粮堆结露预测方法,仅需要粮堆内的粮温数据以及粮食的水分数据即可预测出粮堆是否有结露风险,需要的粮情数据种类较少,便于实施,结露判断结果生成快速;同时本方法能够实现结露位置的图形化表达,从而更加直观的体现出结露情况。
附图说明
图1为本发明所述的粮堆结露预测方法流程图。
图2为本发明实施例中的2018年1月11日某个测温平面的温度场云图。
图3为本发明实施例中的2018年1月15日同一个测温平面的温度场云图。
具体实施方式
下面结合附图对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
如图1所示,本发明提供了一种粮堆结露预测方法,包括如下步骤:
S1:选择任意竖直测温平面,通过插值的方式补全缺失的粮温数据,形成粮温数组。
S2:计算粮食平衡相对湿度;根据粮食平衡相对湿度以及温度,利用露点温度计算模型,计算空气的结露温度。
其中,所述粮食平衡相对湿度ERHr利用CAE模型进行计算;
所述CAE模型为:
式中,ERHr为粮食平衡相对湿度,M为粮食含水率(%),T为粮食温度(℃),A1、A2、B1、B2及D为CAE方程的5个参数,根据粮食种类和吸着类型确定。其中,当粮食籽粒间空气温湿度高于粮食籽粒的温湿度时,粮食处于解析状态;当粮食籽粒间空气温湿度高于粮食籽粒的温湿度时,粮食籽粒处于吸附状态。当表层冷温粮食温度低于计算得出的高温空气露点时,高温空气中的多余的水气将会冷凝在粮食表面,形成结露现象,因此结露区域粮食一般处于吸附状态;CAE方程参数取值如表1所示。
表1 CAE方程参数表
其中,所述露点计算模型为:
td=ERHr×(A+B×t)+C×t-19.2;
式中,A=0.198,B=0.0017,C=0.8400;td为露点温度(℃);t为干球温度(℃);ERHr为粮食平衡相对湿度(%)。
S3:从粮温数组(对应的坐标)的最下面一行依次向上扫描,将每个点位的温度与露点温度进行对比,若某一点位的温度小于露点温度,则预测该点存在结露风险,记录该点的坐标位置。
S4:利用云图生成工具,结合插值完成的该测温平面内所有的粮温数据数组,生成该平面的温度场云图。
S5:将记录的结露坐标位置在云图中连接,实现云图中结露位置预测的图形化表达。
S6:根据当前采样周期结露线的长度,以及相邻两个采样周期的移动速度,判断测温平面是是否有结露风险,将所有测温平面的数据汇总,即可判断出哪些测温平面有结露风险,同时可以直接确定结露位置。
所述结露线长度计算方法为:
统计云图中结露线所占的像素点个数为R,云图的像素尺寸为n×m,其中,n表示云图长度的像素尺寸(云图长度对应的像素点个数),m表示云图高度的像素尺寸(云图长度对应的像素点个数),云图所代表的测温平面的实际尺寸为N×M(长度×高度),则每个像素点的实际尺寸长度C和高度G分别为:
那么云图中结露线长度l为:
l=C×R。
所述结露线移动速度计算方法为:
某一张云图的像素尺寸为n×m,云图中的结露点位置坐标为(xi,yj1),其中(i=0,……,n;j=0,……,m),前一时刻云图中同一列结露点位置坐标为(xi,yj2),则该点的移动距离为:
Si=yj1-yj2;
两张云图所有结露点的平均移动距离Sp为:
式中,p为结露点的个数。
将温度场云图中的结露线移动距离转换到实际空间的移动距离,公式为:
Ss=G×Sp;
则结露线的移动速率为V:
式中,Δd为相邻两张云图间隔的时间。
其中,若某一列不存在结露点,则该列结露点坐标为(xi,0),如果某一列中存在多个结露点,则以位于最上方的结露点坐标为准。
其中,当L≥L0,并且V≥V0时,判断粮堆有结露风险;
L0表示结露线临界长度,V0表示结露线临界移动速度。
实施例
以某储备库的粮仓(平房仓库)粮情数据为例,该粮仓尺寸为46×26m,粮堆高度6m,储藏作物为小麦,入仓时粮食水分为12%。仓内布置66条测温电缆,呈11×6长方形排布,每根电缆上有4个测温点。
结露预测与图形化表达的实施步骤为:
S1:选择测温平面中测温点6×4的平面,通过插值的方式补全缺失的粮温数据,形成粮温数组。
S2:选择CAE模型参数,根据粮情数据和CAE模型参数,计算平衡相对湿度;根据平衡相对湿度以及温度,利用露点温度计算模型,计算空气的结露温度。
其中,M=12%,T=12.5℃,A1、A2、B1、B2及D分别为4.393、4.845、7.843、3.858、203.892,带入下式计算平衡相对湿度,
之后,计算露点温度:
td=ERHr×(A+B×t)+C×t-19.2=10.3;
式中,A=0.198,B=0.0017,C=0.8400;t=18℃;ERHr=63%。
S3:从粮温数组的最下面一行每一列依次向上扫描,判断每个点位的温度与露点温度,若小于露点温度,则表示该点存在结露风险,记录该点的坐标位置。其中,若某一列不存在结露点,则该列结露点坐标为(xi,0)。如果某一列中存在多个结露点,则以位于最上方的结露点坐标为准。
S4:利用云图生成工具,结合插值完成的该平面内所有的粮温数据数组,生成该平面的温度场云图。
S5:将记录的坐标在云图中连接,从而实现云图中结露位置的图形化表达。如图2所示,为该粮仓2018年1月11日某个测温平面的温度场云图,图2中的线L1为云图中的结露线。如图3所示,为该粮仓2018年1月15日同一个测温平面的温度场云图,图3中的线L2为云图中的结露线。
S6:计算结露线长度,分析其移动方向。经计算,图3中的结露线长度为16.5424米,大于设定的结露线临界长度5m,对比2018年1月11日,结露线沿高度方向向下移动速度为0.001848m/s,大于设定的临界移动速度0.001m/s。因此判断该测温平面有结露风险。按照上述方法,可以逐个判断粮堆中哪些测温平面具有结露风险,并且锁定结露点的位置。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
Claims (6)
1.一种粮堆结露预测方法,其特征在于,包括如下步骤:
步骤一、选取粮堆内竖直方向的测温平面,并且获取所述测温平面的粮温数组后制作温度场云图;
步骤二、计算露点温度,将所述粮温数组中每个点位的温度分别与所述露点温度进行对比,并记录结露点的坐标;
其中,所述结露点为温度低于所述露点温度的点位;以及
所述露点温度为:td=ERHr×(A+B×t)+C×t-19.2;
式中,t为干球温度;ERHr为粮食平衡相对湿度,A、B、C分别为公式参数;
步骤三、将所述结露点的坐标在所述温度场云图中依次连接形成结露线;
步骤四、计算所述结露线的长度及结露线的移动速度,
当L≥L0,并且V≥V0时,判断所述测温平面有结露风险;
其中,L为所述结露线的长度,L0为结露线临界长度,V为所述结露线的移动速度,V0为结露线临界移动速度。
2.根据权利要求1所述的粮堆结露预测方法,其特征在于,在所述步骤二中,粮食平衡相对温度为:
式中,ERHr为粮食平衡相对湿度,M为粮食含水率,T为粮食温度,A1、A2、B1、B2和D分别为公式参数。
3.根据权利要求2所述的粮堆结露预测方法,其特征在于,在所述步骤四中,所述结露线的长度为:
l=C×R;
式中,R为温度场云图中结露线所占的像素点个数,C为每个像素点对应的实际长度。
4.根据权利要求3所述的粮堆结露预测方法,其特征在于,所述每个像素点对应的实际长度为:
式中,N为所述测温平面实际长度,n为所述测温平面的温度云场图的长度所占的像素点个数。
5.根据权利要求1或4所述的粮堆结露预测方法,其特征在于,所述结露线的移动速度为:
式中,Δd为温度场云图的采样周期,相邻两个采样周期的所有结露点在温度云场的高度方向上的平均移动距离所占的像素点个数,G为每个像素点对应的实际高度。
6.根据权利要求5所述的粮堆结露预测方法,其特征在于,所述相邻两个采样周期的所有结露点在温度云场的高度方向上的平均移动距离所占的像素点个数为:
式中,p为结露点的个数,Si为相邻采样周期中的同一个结露点移动的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910626913.7A CN110286144B (zh) | 2019-07-12 | 2019-07-12 | 一种粮堆结露预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910626913.7A CN110286144B (zh) | 2019-07-12 | 2019-07-12 | 一种粮堆结露预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110286144A true CN110286144A (zh) | 2019-09-27 |
CN110286144B CN110286144B (zh) | 2020-07-21 |
Family
ID=68021384
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910626913.7A Active CN110286144B (zh) | 2019-07-12 | 2019-07-12 | 一种粮堆结露预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110286144B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111149523A (zh) * | 2019-12-30 | 2020-05-15 | 黑龙江八一农垦大学 | 集装箱存放玉米的保存方法 |
CN111307866A (zh) * | 2020-03-20 | 2020-06-19 | 河南工业大学 | 一种粮食结露临界参数判断方法 |
CN112763399A (zh) * | 2020-12-28 | 2021-05-07 | 西安特种设备检验检测院 | 一种锅炉低温受热面烟气腐蚀风险区域的检测方法 |
CN113237918A (zh) * | 2021-05-20 | 2021-08-10 | 北京盈丰利泰科贸有限公司 | 一种烟叶密封堆垛结露提前预警方法、装置及电子设备 |
CN114113216A (zh) * | 2020-08-26 | 2022-03-01 | 宝山钢铁股份有限公司 | 带钢板卷库结露预警方法 |
CN116818835A (zh) * | 2023-08-30 | 2023-09-29 | 中储粮成都储藏研究院有限公司 | 一种基于粮仓温度的结露判定方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03195883A (ja) * | 1989-12-26 | 1991-08-27 | Sanyo Electric Co Ltd | 高湿度低温貯蔵庫の露受け装置 |
CN102156007A (zh) * | 2011-02-14 | 2011-08-17 | 国家粮食局科学研究院 | 粮情预测系统、粮情控制系统及粮堆安全状态判别方法 |
CN104007776A (zh) * | 2014-05-23 | 2014-08-27 | 国家粮食局科学研究院 | 一种基于温度湿度场耦合的粮堆结露预警防控系统及方法 |
CN106980332A (zh) * | 2017-05-23 | 2017-07-25 | 安徽博微长安电子有限公司 | 一种储粮状态的检测方法 |
CN108333216A (zh) * | 2018-04-27 | 2018-07-27 | 中储粮成都储藏研究院有限公司 | 粮堆结露预警检测仪 |
CN109060873A (zh) * | 2018-03-16 | 2018-12-21 | 河南工业大学 | 一种粮食质热传递和结露参数检测装置及方法 |
-
2019
- 2019-07-12 CN CN201910626913.7A patent/CN110286144B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03195883A (ja) * | 1989-12-26 | 1991-08-27 | Sanyo Electric Co Ltd | 高湿度低温貯蔵庫の露受け装置 |
CN102156007A (zh) * | 2011-02-14 | 2011-08-17 | 国家粮食局科学研究院 | 粮情预测系统、粮情控制系统及粮堆安全状态判别方法 |
CN104007776A (zh) * | 2014-05-23 | 2014-08-27 | 国家粮食局科学研究院 | 一种基于温度湿度场耦合的粮堆结露预警防控系统及方法 |
CN106980332A (zh) * | 2017-05-23 | 2017-07-25 | 安徽博微长安电子有限公司 | 一种储粮状态的检测方法 |
CN109060873A (zh) * | 2018-03-16 | 2018-12-21 | 河南工业大学 | 一种粮食质热传递和结露参数检测装置及方法 |
CN108333216A (zh) * | 2018-04-27 | 2018-07-27 | 中储粮成都储藏研究院有限公司 | 粮堆结露预警检测仪 |
Non-Patent Citations (2)
Title |
---|
王东华: "粮堆的结露和预防处理", 《仓储科技》 * |
章铖: "温湿度对稻谷粮堆结露的影响及实仓结露预警", 《中国粮油学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111149523A (zh) * | 2019-12-30 | 2020-05-15 | 黑龙江八一农垦大学 | 集装箱存放玉米的保存方法 |
CN111307866A (zh) * | 2020-03-20 | 2020-06-19 | 河南工业大学 | 一种粮食结露临界参数判断方法 |
CN111307866B (zh) * | 2020-03-20 | 2022-07-12 | 河南工业大学 | 一种粮食结露临界参数判断方法 |
CN114113216A (zh) * | 2020-08-26 | 2022-03-01 | 宝山钢铁股份有限公司 | 带钢板卷库结露预警方法 |
CN114113216B (zh) * | 2020-08-26 | 2023-10-13 | 宝山钢铁股份有限公司 | 带钢板卷库结露预警方法 |
CN112763399A (zh) * | 2020-12-28 | 2021-05-07 | 西安特种设备检验检测院 | 一种锅炉低温受热面烟气腐蚀风险区域的检测方法 |
CN112763399B (zh) * | 2020-12-28 | 2021-09-21 | 西安特种设备检验检测院 | 一种锅炉低温受热面烟气腐蚀风险区域的检测方法 |
CN113237918A (zh) * | 2021-05-20 | 2021-08-10 | 北京盈丰利泰科贸有限公司 | 一种烟叶密封堆垛结露提前预警方法、装置及电子设备 |
CN116818835A (zh) * | 2023-08-30 | 2023-09-29 | 中储粮成都储藏研究院有限公司 | 一种基于粮仓温度的结露判定方法 |
CN116818835B (zh) * | 2023-08-30 | 2023-11-21 | 中储粮成都储藏研究院有限公司 | 一种基于粮仓温度的结露判定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110286144B (zh) | 2020-07-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110286144A (zh) | 一种粮堆结露预测方法 | |
CN104165978B (zh) | 提高水质监测数据可靠性的方法 | |
An et al. | Quantifying time-series of leaf morphology using 2D and 3D photogrammetry methods for high-throughput plant phenotyping | |
CN104040311B (zh) | 用于房间内温度监视的系统及方法 | |
CN104007776B (zh) | 一种基于温度湿度场耦合的粮堆结露预警防控系统及方法 | |
US11624825B2 (en) | Object tracking | |
CN103344742B (zh) | 一种粮仓储粮虫霉活动危害度预测的方法 | |
CN106885846B (zh) | 树木缺陷检测装置及检测方法 | |
de Moraes Frasson et al. | Rainfall interception by maize canopy: Development and application of a process-based model | |
CN105224714B (zh) | 气象数据的处理方法及装置 | |
CN106802164A (zh) | 一种适用于粮情监测的绝对水势计算方法及云图生成方法 | |
CN106017546B (zh) | 一种原位测量粮仓单元装置 | |
CN102870137A (zh) | 用于基于2d运动估计相对3d速度和加速度投影函数的数据挖掘方法和系统 | |
CN104729097B (zh) | 一种热水系统及其显示控制方法、装置 | |
CN106247935A (zh) | 一种用于体积测量系统的探测物体积计算方法 | |
Gleicher et al. | Interpreting three-dimensional spore concentration measurements and escape fraction in a crop canopy using a coupled Eulerian–Lagrangian stochastic model | |
CN111310317B (zh) | 基于大数据与插值预测的粮仓空时温度场预测方法及装置 | |
CN113052399A (zh) | 温度监控系统、方法、装置、存储介质及电子设备 | |
CN109579827A (zh) | 一种基于弧形阵列的磁性目标探测和定位方法 | |
CN108764749A (zh) | 一种适于储粮宜存评价的3i策略方法 | |
CN105137041B (zh) | 土壤参数空间分布的监测方法和系统 | |
CN110810027A (zh) | 一种基于粮温的粮仓内环流控温方法 | |
CN110348520A (zh) | 一种基于气象数据预测仓储粮堆表层温度的方法和系统 | |
Liang et al. | Evaluation of ECOSTRESS evapotranspiration estimates over heterogeneous landscapes in the continental US | |
CN108268646A (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 |