一种适于储粮宜存评价的3I策略方法与流程

本发明涉及粮仓粮情监测领域，具体涉及一种适于储粮宜存评价的3i策略方法。

背景技术

利用粮仓储备粮食，对保障国家粮食的供给平衡以及保持粮食品质至关重要。储粮监管是利用人工以及信息化技术手段，监察储粮是否按计划进行仓储、调用、轮换的管理过程。目前，对储粮的监管方法主要形式是以人工为主的清仓查库，一些有实力的粮食仓储企业开始通过在粮仓内安装视频监控设备辅助监管，行业内尚缺乏对粮仓内部质量实时宜存评价和质量追溯的方法。

技术实现要素：

本发明设计开发了一种适于储粮宜存评价的3i策略方法，本发明的发明目的是利用粮情监测系统采集的粮仓温度信息进行不同周期的等效积温计算，根据实时滚动积温直接对粮仓微生物危害、昆虫危害和品质劣变的发生进行宜存评价和质量追溯。

本发明提供的技术方案为：

一种适于储粮宜存评价的3i策略方法，包括如下步骤：

步骤一、确定所述储粮仓储的采样范围、测温点和积温周期长度；其中，所述积温周期长度包括短周期、中周期和长周期；

步骤二、在所述采样范围和所述采样周期内进行实时温度计算；

步骤三、通过如下公式进行等效积温计算：

式中，为短周期等效积温，t0s为短周期等效积温的温度起点，ts,d为粮仓位置的实际温度，t为温度ts,d持续时间长度；

式中，为中周期等效积温，t0m为中周期等效积温的温度起点，ts,d为粮仓位置的实际温度，t表示温度ts,d持续时间长度；

式中：ct^pl为长周期等效积温，t0l为长周期等效总积温的温度起点，ts,d为粮仓位置的实际温度，tθ表示温度ts,d持续时间长度，δ表示时间tθ的个数；

步骤四、根据所述短周期等效积温进行微生物发生判断与预测；

根据所述中周期等效积温进行昆虫发生判断与预测；以及

根据所述长周期等效积温进行粮食质量宜存评价。

优选的是，在所述步骤四中，根据所述短周期等效积温进行微生物发生判断，包括：

如果且则判断所述短周期内测试点产生微生物；

如果或者则判断所述短周期内测试点未产生微生物；

其中，ctw为所述短周期微生物发生的积温阈值，为粮仓所述短周期内平均湿度，为粮仓所述短周期内微生物产生的湿度阈值；

根据所述短周期等效积温进行微生物发生预测，包括：

如果ts,d≥tcw,d，则预测所述短周期内测试点可能产生大量微生物；

如果ts,d＜tcw,d，则预测所述短周期内测试点不可能产生大量微生物；

其中，ts,d为测试点的粮食温度，tcw,d为微生物发生温度。

优选的是，在所述步骤四中，根据所述中周期等效积温进行昆虫发生判断，包括：

如果且则判断所述中周期内测试点可能产生昆虫；

如果或者则判断所述中周期内测试点不可能产生昆虫；

其中，ctq为所述中周期昆虫发生的积温阈值，为粮仓所述中周期内平均湿度，为粮仓所述中周期内昆虫产生的湿度阈值；

根据所述中周期等效积温进行昆虫发生预测，包括：

如果tempd≥ts,d，则预测所述中周期内测试点可能产生昆虫；

如果tempd＜ts,d，则预测所述中周期内测试点不可能产生昆虫；

其中，ts,d为测试点的粮食温度，tempd为昆虫发生温度。

优选的是，在所述步骤四中，根据所述长周期等效积温进行粮食质量宜存评价，包括：

如果则判断测试点及周边粮食处于优质态；

如果则判断测试点及周边粮食处于宜存态；

如果则判断测试点及周边粮食处于轻度陈化态；

如果则判断测试点及周边粮食处于陈化态；

其中，cts1＝(t1-t0l)×α，cts2＝(t2-t0l)×α，cts3＝(t3-t0l)×α，cts1，cts2，cts3分别宜存、轻度不宜存和重度不宜存等效积温，t0l为长期等效积温的温度起点，t1为谷物长期储藏宜存温度，t2为谷物长期储藏轻度不宜存温度，t3为谷物长期储藏重度不宜存温度，α为储藏时间；

当p4>0时，预测全仓粮食变陈化的时间t＝τ×(1-p4)；

当p4＝0时，预测全仓粮食变陈化的时间t＞τ×(1-p4)；

其中，τ为从储藏开始到粮食出现陈化的时间，p4为陈化率，v4为处于陈化态的粮食体积，v为全仓内粮食总体积。

优选的是，所述短周期微生物发生的积温阈值ctw计算过程包括：

采集粮食温度ts,d、粮食湿度rhs,d，通过emc模型计算得到粮食水分ws,d后，再通过微生物模型计算得到此时的温度tcw,d，对所述短周期的时间内的温度tcw,d进行积温计算得到所述积温阈值ctw；

所述短周期微生物发生的湿度阈值计算过程包括：

采集粮食温度ts,d，通过微生物模型计算得到此时的水分ws,d后，再通过emc模型计算得到此时的湿度rhcw,d，进而得到所述短周期内的时间内的所述湿度阈值

其中，所述emc模型为emc＝[ln(1-erh)-lnerh-a-ct]/b；式中，erh为谷物的相对湿度；emc为谷物湿基含水率，即谷物水分；t为环境温度，a，b，c为拟合系数；

所述微生物模型为式中，a1、b1、c1、d1、e1为常数，dd表示间隔天数。

优选的是，所述中周期昆虫发生的积温阈值ctq计算过程包括：

采集粮仓湿度rhs,d，选择昆虫模型，确定影响粮食品质的昆虫个数sv，通过昆虫模型，计算此时粮食昆虫发生温度为tempd，对所述中周期的时间内的温度tempd进行积温计算得到所述积温阈值ctq；

采集粮仓温度ts,d，选择昆虫模型，确定影响粮食品质的昆虫个数sv，通过昆虫模型，计算此时粮食昆虫发生的湿度为rhcq,d，进而得到所述中周期内的时间内的所述湿度阈值

优选的是，所述昆虫模型中的玉米象鼻虫害发生模型为

式中，sv为玉米象从虫卵到成虫的存活数量，temp为粮食温度，rh为粮食湿度，a、b、c、d、e、f、g为常数。

优选的是，在所述步骤二中，在所述采样范围内进行实时温度计算过程包括如下步骤：

确定采样范围内的第一方向、第二方向和第三方向，由所述第一方向和所述第二方向组成第一平面，由所述第一方向和所述第二方向组成第一平面，所述第一方向和所述第三方向组成第二平面，由所述第二方向和所述第三方向组成第三平面；其中，所述第一方向、所述第二方向和所述第三方向分别相互垂直，所述第一平面、所述第二平面和所述第三平面分别相互垂直；

计算过程如下：

所述第一方向上平均温度为：

0＜i≤n或0＜j≤m；

所述第二方向上平均温度为：

0＜i≤n或0＜k≤l；

所述第三方向上平均温度为：

0＜j≤m或0＜k≤l；

所述第一平面上平均温度为：

所述第二平面上平均温度为：

所述第三平面上平均温度为：

所述粮仓整体平均温度为：

优选的是，所述短周期为3～14天，所述中周期为21～42天，所述长周期为以粮食储藏开始为起点持续到评测分析结束。

优选的是，在所述步骤二中，在所述采样范围和所述采样周期内进行实时温度计算时，所述短周期和所述中周期的测试采用周期滚动测试的方式。

本发明与现有技术相比较所具有的有益效果：利用短周期(滚动)等效积温实现了储粮霉变的微生物爆发情况的预测和判断，利用中周期(滚动)等效积温实现了储粮昆虫发生情况的预测和判断，以长周期等效总积温实现储粮宜存评价与质量追溯，解决了仓储谷物实时快捷评价方法缺失的问题，降低了粮情品质检测的操作难度、减轻了劳动力。

附图说明

图1为本发明所述的流程图。

图2为本发明所述的实施例中粮仓东西向横剖面与测温点布置图。

图3为本发明所述的实施例中粮仓南北向横剖面与测温点布置图。

图4为本发明所述的实施例中粮仓坐标系示意图。

图5为本发明所述的实施例中各平面等效积温结果示意图。

图6为本发明所述的实施例中微生物发生的储粮温度与水分关系示意图。

图7为本发明所述的实施例中xoy的四个平面的日均等效积温线示意图。

具体实施方式

下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。

如图1所示，3i(积分，integration)是以三个不同周期对储粮温度进行积分，基于上述理论，本发明提供了一种适于储粮宜存评价的3i策略方法，包括如下过程：

步骤一、根据粮仓粮情监测系统，建立测温点三维矩阵坐标系。平房仓以粮仓的一个底角为坐标系原点，构成底角三边分别为x、y、z轴正方向，x、y、z轴正方向测温点个数分别为n、m和l；圆筒仓以粮仓外切正方体的底边某一顶点为坐标系原点，构成底角三边分别为x、y、z轴正方向，x、y、z轴正方向测温点个数分别为n、m和l；

步骤二、采集粮情数据；

步骤三、去除干扰数据；

步骤四、选择微生物与昆虫模型；

步骤五、确定积温周期长度，短周期为pes天、中周期为pem天、长周期为pel天；在本实施例中，作为一种优选，3<pes<14，21<pem<42，长周期以储藏开始为起点持续到现在；

步骤六、实时温度计算；根据所述步骤一中建立的测温点坐标系，仓内各测温点的坐标表示为(x,y,z)＝(i,j,k)，其中，0<i≤n，0<j≤m，0<k≤l，n、m、l分别为粮仓x、y、z正方向的测温点数，将粮仓分为三种平面，分别为平行于xoy平面，点集表示为(1～n，1～m，k)，平行于xoz平面，点集表示为(1～n，j，1～l)，平行于yoz平面，点集表示为(i，1～m，1～l)；将粮仓第d天测温点(i,j,k)处的温度表示为t(xi,yj,zk,τd)，d表示时间，取值1～τ，τ为从储藏开始到现在的储藏天数；

第d天粮仓内各条线(测温点组成的各条线)平均温度计算公式为：

0＜i≤n或0＜j≤m；

0＜i≤n或0＜k≤l；

0＜j≤m或0＜k≤l；

第d天粮仓内各平面(测温点组成的平面)平均温度计算公式为：

步骤七、第d天粮仓整体均温计算公式为：

其中，上述点、线、面、粮仓整体的实际温度统称为温度ts,d；

步骤七、短期(滚动)、中期(滚动)等效积温计算；设置短、中、长周期温度起点分别为t0s，t0m，t0l；从储藏开始，以每一天为起点，向后分别截取短周期内温度、中周期内温度，与时间相乘后累加作为该位置这个周期内(滚动)等效积温；

各位置短周期(滚动)等效积温计算公式如：

式中，为短周期(滚动)等效积温，t0s为短周期(滚动)等效积温的温度起点；t为温度ts,d持续时间长度，单位为天；d取值为1～pes(从开始到储藏的日期pes)；

各位置中周期(滚动)等效积温计算公式如：

式中，为中周期(滚动)等效积温，t0m为中周期(滚动)等效积温的温度起点；t表示温度ts,d持续时间；d取值为1～pem(从开始到储藏的日期pem)；

步骤八、对数据进行评测分析；

(1)、微生物发生判断与预测：

微生物发生判断方法为：若且则判断短周期pes天内该位置产生了微生物，若或则判断短周期pes天内该位置没有产生微生物；其中，ctw为短周期微生物发生的积温阈值，为粮仓pes天内平均湿度，为粮仓短周期pes天内微生物产生的湿度阈值，ctw、由微生物发生模型与emc模型求得；仓内点、线、平面或粮仓整体的判断方法相同；记录微生物发生日期；

微生物发生预测方法为：根据粮食湿度rhs,d与微生物模型，求得微生物发生温度tcw,d；若ts,d≥tcw,d，则可预测短周期pes天内该位置可能产生大量微生物；若ts,d＜tcw,d，则可预测短周期pes天内该位置不可能产生微生物；仓内点、线、平面或粮仓整体的预测方法相同；

其中，湿度阈值计算方法为：根据当前粮食温度ts,d与微生物发生公式求得微生物发生的水分ws,d，将温度与水分带入emc模型中，求得每天的该位置粮食湿度rhcw,d，依次求得pes天内的该位置粮食湿度均值为作为粮仓短周期pes天内微生物产生的湿度阈值；

积温阈值ctw计算方式为；根据粮食的emc模型与粮食温度ts,d、粮食湿度rhs,d，计算得到此时的粮食水分ws,d，将粮食水分ws,d带入微生物模型中求得此时的该位置微生物发生温度tcw,d，对短周期pes天内的该位置微生物发生温度tcw,d积温得到ctw；

谷物三参数emc模型为：

emc＝[ln(1-erh)-lnerh-a-cte]/b；

式中：erh为谷物相对湿度，单位为％；emc为谷物湿基含水率，也即谷物水分，单位为％；te环境温度，单位为℃，a，b，c是不同粮种解吸和吸附不同状态下的拟合系数；

粮食的微生物发生温度、水分关系模型为：

式中，a1、b1、c1、d1、e1为常数，dd表示间隔天数，时间与短周期天数相同；

(2)、昆虫发生判断与预测：

昆虫发生判断方法为：进行中周期(滚动)等效积温计算，判断若且则该位置附近中周期pem天内可能产生昆虫；若或则该位置附近中周期pem天内没有可能产生昆虫，ctq为昆虫发生的积温阈值，为粮仓中周期pem天内平均湿度，为粮仓pem天内昆虫产生的平均湿度，也即湿度阈值；ctq、由昆虫发生模型求得；记录昆虫产生的日期；

昆虫产生预测方法为：选择昆虫产生模型与昆虫个数sv只/kg(sv<4)，计算当前湿度下的昆虫产生的温度为tempd，若tempd≥ts,d，则粮食在该位置附近容易产生昆虫，若tempd＜ts,d，则粮食在该位置附近不容易产生昆虫；

其中，积温阈值ctq计算方法为：选择昆虫模型，确定影响粮食品质的昆虫个数sv只/kg(sv<4)，然后根据昆虫模型与粮仓湿度rhs,d，计算此时粮食发生昆虫的温度为tempd(tempd可代表点、线、面或粮仓的温度)，对温度tempd进行中周期pem天积温得到积温阈值ctq；

湿度阈值计算方法为：根据选择的昆虫模型和昆虫个数sv，利用粮仓温度ts,d求得昆虫发生的湿度rhcq,d，中周期pem天内的湿度求均值得到

(3)、粮情宜存评价与质量追溯的方法为：

选择粮仓某一类位置(点、线、面或整个粮仓)进行长周期等效总积温计算，即从储藏开始至当前时间的总积温ct^pl，计算公式为：

式中：t0l为长周期等效总积温的温度起点；ts,d为粮仓位置的实际温度，tθ为该位置粮食温度ts,d持续时间，δ表示时间tθ的个数，ct^pl为长周期等效总积温；

设定粮食质量评价积温阈值；根据粮食(玉米、稻谷、小麦)储存品质判定规则现行国家标准(gb/t20570-202015、gb/t20569-2006、gt/b20571-2006)中设定的粮食宜存、轻度不宜存、重度不宜存的品质指标，结合参考文献(华北地区实仓储藏稻谷的品质变化及控制研究)中储藏时间、温度与谷物品质的关系，设定谷物长期储藏宜存温度为t1，轻度不宜存温度为t2，重度不宜存温度为t3，储藏时间α天(以t1、t2、t3温度储藏α天后粮食品质达到国家标准宜存、轻度不宜存和重度不宜存线)，则α天的宜存、轻度不宜存、重度不宜存积温线分别为：

cts1＝(t1-t0l)×α；

cts2＝(t2-t0l)×α；

cts3＝(t3-t0l)×α；

式中：t0l为长期等效积温的温度起点，cts1，cts2，cts3分别宜存、轻度不宜存和重度不宜存等效积温；

根据设定积温阈值线划分粮食评价时区：

若为优质时区，该位置及周边粮食处于优质态；

若为宜存时区，该位置及周边粮食处于宜存态；

若为轻度陈化时区，该位置及周边粮食处于轻度陈化态；

若为陈化时区，该位置及周边粮食处于陈化态；

点、线、面以及粮仓的粮食宜存评价时区同上；

粮仓粮食质量宜存评价与质量追溯；根据上述时区评价各位置粮食所处状态，根据粮仓尺寸结构及各位置的分布，计算各位置(点、线、面或粮仓中的某一类)代表粮食的体积，统计各位置粮食所处时区，各时区粮食体积分别为v1，v2，v3，v4，粮堆总体积为v，则优质率p1、宜存率p2、轻度陈化率p3、陈化率p4计算公式如下所示；若p4>0，利用(1-p4)与已储藏时间τd的乘积(τd为从储藏开始到粮食出现陈化的时间)，预测全仓粮食变陈化的时间为τ×(1-p4)，若p4＝0，则预测陈化时间大于τ；

步骤九、生成粮情宜存评价与质量追溯报告；以文档形式生成报告，该报告记录操作人员进行粮食质量宜存评价与追溯的具体日期，以及所在粮仓的基本信息(粮仓尺寸、所在地、储藏时间等信息)，同时记录微生物、昆虫发生的次数以及具体日期，粮食质量变化的日期(包括粮情质量变化日期，以及预测粮情质量发生变化的时间)。

在另一种实施例中，中周期中的昆虫模型选取玉米象鼻虫害发生模型为

式中，sv为玉米象从虫卵到成虫的存活数量，依据参考文献(玉米象取食对新稻谷储藏品质的影响)，sv取值范围为(sv<4只/kg)时玉米象对谷物品质影响较小，a＝-0.124731，b＝1.14602，c＝25.04，d＝6.58875，e＝2.78128，f＝68.3030，g＝0.314531，temp为粮食温度，rh为粮食湿度，粮食湿度以粮仓湿度代替。

实施例

本发明提供了一种基于温度积分的仓储粮食质量评测分析方法，包括如下过程：

步骤一、如图2～4所示，根据现有粮情系统建立粮仓测温点坐标系；粮仓东西方向长47.5m，南北方向宽26m，仓房高8m，粮面高6m；内置测温电缆78根，东西方向13根，间距3.75m，两侧离墙面1.25m；南北方向6根，间距5m，两侧离墙面0.5m；第一层测温点距底面0.5m，第四层测温点距粮面高度0.3m，每根测温电缆上点均匀分布，间隔约1.7m；以粮仓西北方向底顶点作为坐标原点，竖直向上为z轴正向，向东为x轴正向，向南为y轴正向。根据上述坐标系，建立传感器测温点xyz数据三维，矩阵x方向n个测温点，y方向m个测温点，z方向l个测温点；

步骤二、采集各测温点温度数据，建立温度数据三维矩阵；采集时间为2014年5月6日-2015年5月6日；

步骤三、去除干扰及非正常数据；

步骤四、选择微生物模型与昆虫模型，确定粮食质量判断阈值；

步骤五、确定积温周期的长度，短周期为pes天、中周期为pem天、长周期为pel天，其中：3<pes<14，21<pem<42，长周期积温以2014年5月6日为积温第一天至查询日期截止。

步骤六、实时温度计算；根据s1建立的传感器测温点坐标系，仓内各测温点的坐标表示为(x,y,z)＝(i,j,k),其中0<i≤13，0<j≤6，0<k≤4，13、6、4分别为粮仓x、y、z方向的测温点数。将粮仓分为三类平面，分别为平行于xoy平面，点集表示为(1～13，1～6，k)，平行于xoz平面，点集表示为(1～13，j，1～4)，平行于yoz平面，点集表示为(i，1～6，1～4)；分别求各平面所有测温点采集温度的均值作为该平面的温度，又因为粮仓每天均采集温度，因此当天温度表示为t(xi,yj,zk,τd)，d表示时间，取值1～τ，τ为从储藏开始到现在的储藏天数(一天采集一次数据)，该粮库储粮365天，则τ取值为365；选择短周期为7天，中周期为21天；特别的：某一天各平面均温计算公式为：

以上各平面均分统称为ts,d；

步骤七、短(7天)、中(21天)周期(滚动)等效积温计算；设置短、中、长周期积温起点分别为t0s，t0m，t0l；从储藏开始，各平面以每一天为起点，向后分别截取7、21天温度进行温度累加，与时间相乘后作为该平面这个周期内的7天(滚动)等效积温和21天(滚动)等效积温；各平面7天(滚动)等效积温计算公式如：

式中，t0s为长周期等效总积温的温度起点，t取值为1天，d表示以储藏开始之后的每天，取值为1～358；

选择五个平面进行7天(滚动)等效积温，结果如图5所示；

各平面21天(滚动)等效积温计算公式如：

式中，t0m为长周期等效总积温的温度起点，t取值为1天，d取值为1～344；

步骤八、对数据进行评测分析；

(1)、微生物发生判断与预测：

进行短周期(滚动)等效积温计算，若且则判断7天内该xoy平面产生了微生物，若或则判断7天内该xoy平面没有产生微生物，ctw为短周期微生物发生的等效积温，为粮仓7天内平均湿度，为粮仓7天内微生物产生的平均湿度，ctw、由微生物发生模型与emc模型求得；仓内点、线、平面或粮仓整体的判断方法相同；记录微生物产生的日期；

所述微生物发生预测方法为：

根据粮食湿度rhs,d与微生物模型，求得微生物发生温度tcw,d。若ts,d≥tcw,d，则可预测7天内该位置可能产生大量微生物；若ts,d＜tcw,d，则可预测7天内该位置不可能产生微生物；仓内点、线、平面或粮仓整体的预测方法相同；

湿度阈值计算方法为：根据当前粮食温度txoy,k,d与微生物发生公式求得微生物发生的水分ws，将粮食温度txoy,k,d与水分ws带入emc方程中，求得每天的该位置粮食湿度rhcw,d，依次求得7天内的湿度rhcw,d的均值为

微生物发生的积温阈值ctw计算方式为，根据粮食的emc模型与粮食温度txoy,k,d、湿度rhs,d，计算得到此时的粮食水分ws，将粮食水分ws带入微生物模型中求得此时的粮食微生物发生的温度tcxoy,k,d，对7天内的tcxoy,k,d积温得到ctw；

谷物三参数emc模型为：emc＝[ln(1-erh)-lnerh-a-ct]/b；式中：erh为谷物相对湿度，单位为％；emc为谷物湿基含水率，即谷物水分，单位为％；t为环境温度，单位为℃；a，b，c是不同粮种解吸和吸附不同状态下的拟合系数，如表2所示；

表2拟合系数

所述粮食的微生物发生温湿度、水分关系模型为：

式中，a1、b1、c1、d1、e1为常数，dd表示间隔天数，同pes时间；

特别的，粮食为玉米时，a1＝-19.644，b1＝2.514，c1＝1.13×10⁹，d1＝5.398，e1＝126.609，dd表示间隔天数，这里取值为7；

选取dd为7天计算微生物发生的储粮温度与水分关系如图6所示。

昆虫发生判断方法为：进行中周期(滚动)等效积温计算，判断若且则该位置附近这21天内可能产生昆虫；若且则该位置附近这21天内没有可能产生昆虫；ctq为昆虫发生的积温阈值，为粮仓pes天内平均湿度，为粮仓pes天内昆虫产生的平均湿度，也即湿度阈值。ctq、由昆虫发生模型求得。记录昆虫产生的日期；

昆虫产生预测方法为：选择昆虫产生模型与昆虫个数sv只/kg(sv<4)，计算当前湿度下的昆虫产生的温度为tempxoy,k,d，比较若tempxoy,d≥txoy,k,d，则粮食在该位置附近容易产生昆虫，若tempxoy,k,d＜txoy,k,d，则粮食在该位置附近不容易产生昆虫；

积温阈值ctq计算方法为：选择昆虫模型，确定影响粮食品质的昆虫个数sv只/kg(sv<4)，然后根据昆虫模型与粮仓湿度rhs,d，计算此时粮食发生昆虫的温度为tempxoy,k,d(tempxoy,k,d可代表点、线、面或粮仓的温度),对温度tempxoy,k,d进行21积温得到积温阈值ctq；

湿度阈值计算方法为：根据选择的昆虫模型，以及昆虫个数sv＝3，利用粮仓温度ts,d求得昆虫发生的湿度rhcq,d，21天内的湿度求均值得到xoz平面、yoz平面预测方法相同；

玉米象鼻虫害发生公式：

式中，sv为玉米象虫卵到成虫的存活数量，依据参考文献(玉米象取食对新稻谷储藏品质的影响)，sv取值范围为(sv<4只/kg)时玉米象虫对谷物品质影响较小；在本实施例中，a＝-0.124731，b＝1.14602，c＝25.04，d＝6.58875，e＝2.78128，f＝68.3030，g＝0.314531；temp为粮食温度，rh为粮食湿度，粮食湿度以粮仓湿度代替；

(3)、所述粮食宜存评价和质量追溯的方法为：

选择平面进行周期等效总积温，即从储藏开始至当前时间的总积温，各平面的长周期等效总积温计算公式为：

式中：tθ选取为1，δ取值为τ。

设定粮食质量宜存评价的积温阈值；根据粮食(玉米、稻谷、小麦)储存品质判定规则现行国家标准中设定的粮食宜存、轻度不宜存、重度不宜存的品质指标，结合参考文献(华北地区实仓储藏稻谷的品质变化及控制研究)中储藏时间、温度与谷物品质的关系，设定谷物长期储藏宜存温度为t1，轻度不宜存温度为t2，重度不宜存温度为t3，储藏时间1095天(以t1、t2、t3温度储藏1095天后粮食品质达到国家标准宜存、轻度不宜存和重度不宜存线)，则1095天宜存、轻度不宜存、重度不宜存等效积温线分别为：

cts1＝(t1-t0l)×1095；

cts2＝(t2-t0l)×1095；

cts3＝(t3-t0l)×1095；

式中：t0l为长周期等效积温的温度起点，cts1，cts2，cts3分别宜存、轻度不宜存和重度不宜存等效积温线；

日均等效积温计算；将宜存、轻度不宜存、重度不宜存等效积温线换算成粮食质量评价日均等效积温线，ctr1＝cts1/β，ctr2＝cts2/β，ctr3＝cts3/β，β取值为1～1095之间的整数。同时将平面长期等效总积温换算到日均等效总积温，即

根据上述粮食质量评价日均等效积温线划分粮食质量宜存评价时区。若ctrxoy,k,d＜ctr1，为优质时区，该层及周边粮食处于优质态，若tr1＜trxoy,k,d＜tr2，为宜存时区，该层及周边粮食处于宜存态，若ctr2＜ctrxoy,k,d＜ctr3，为轻度陈化时区，该层及周边粮食处于轻度陈化态，若ctr3＜ctrxoy,k,d，为陈化时区，该层及周边粮食处于陈化态；利用xoz或yoz平面进行质量评价时，方法同上，且质量评价只选择三类平面中的一类；如图7所示，小于ctr1为优质时区，位于ctr1与ctr2之间为宜存时区，位于ctr2与ctr3之间为轻度陈化时区，大于ctr3为陈化时区；图中显示x0y的四个平面的日均等效积温线，分别为ctrxoy,1,d、ctrxoy,2,d、ctrxoy,3,d、ctrxoy,4,d；

根据上述设定的评价时区，统计各平面粮食日均等效总积温的时区分布，计算各平面所代表的粮食体积，如：选择为xoy平面，根据测温点排布以及粮仓尺寸结构计算xoy平面各层代表的体积，分别为v1＝1668.25m³，v2＝3356.25m³，v3＝3356.25m³，v4＝1420.25m³粮仓总储粮体积为9880m³，储藏时间到达a日期时，第四层由优质态变为宜存态，则粮食宜存率优质率为85.62％。储藏时间到b日期时第四层由宜存时区变为宜存时区，则粮食轻度陈化率为14.38％，优质率为85.62％；若p4>0,利用(1-p4)与已储藏时间τd的乘积(τd为从储藏开始到粮食出现陈化的时间长度)，预测全仓粮食变陈化的时间为τd×(1-p4)，若p4＝0，则预测陈化时间大于τd。由图可知τd约为330天，陈化率为14.38％，则预测整仓粮食陈化的时间为282天；

步骤九、生成粮情宜存评价与质量追溯报告；以文档形式生成报告，该报告记录操作人员进行粮食质量宜存评价与追溯的具体日期，以及所在粮仓的基本信息(粮仓尺寸、所在地、储藏时间等信息)，同时记录微生物、昆虫发生的次数以及具体日期，粮食质量变化的日期(包括粮情质量变化日期，以及预测粮情质量发生变化的时间)。

尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。