CN108981960A - A kind of AD strategy process suitable for grain storage supervision - Google Patents

A kind of AD strategy process suitable for grain storage supervision Download PDF

Info

Publication number
CN108981960A
CN108981960A CN201810566299.5A CN201810566299A CN108981960A CN 108981960 A CN108981960 A CN 108981960A CN 201810566299 A CN201810566299 A CN 201810566299A CN 108981960 A CN108981960 A CN 108981960A
Authority
CN
China
Prior art keywords
grain
temperature
silo
days
layer
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.)
Pending
Application number
CN201810566299.5A
Other languages
Chinese (zh)
Inventor
吴子丹
吴文福
韩峰
崔宏伟
朱浩天
陈中旭
徐岩
刘哲
张亚秋
陈龙
秦骁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JIDA SCIENCE APPARATUS CO Ltd CHANGCHUN
Jilin University
Original Assignee
JIDA SCIENCE APPARATUS CO Ltd CHANGCHUN
Jilin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by JIDA SCIENCE APPARATUS CO Ltd CHANGCHUN, Jilin University filed Critical JIDA SCIENCE APPARATUS CO Ltd CHANGCHUN
Priority to CN201810566299.5A priority Critical patent/CN108981960A/en
Publication of CN108981960A publication Critical patent/CN108981960A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K13/00Thermometers specially adapted for specific purposes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology

Landscapes

  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Storage Of Harvested Produce (AREA)

Abstract

The invention discloses a kind of AD strategy process suitable for grain storage supervision, temperature average value and grain depot location universe temperature threshold value in grain depot location one day including step 1: are determined according to silo geography information;Step 2: determining each layer grain temperature threshold value of silo;Step 3: according to the grain temperature in sampling period acquisition silo, constituting three-dimensional grain temperature matrix;Step 4: when grain temperature is not in every layer of grain temperature threshold range or when in the universe temperature threshold range of grain depot location in silo, grain temperature abnormality in silo.AD strategy process of the present invention suitable for grain storage supervision, can judge whether grain temperature is abnormal in silo according to grain grain temperature variation amplitude inside annual temperature period and silo, can also judge whether the temperature rise of grain or temperature drop are abnormal in silo according to the differential variation rate of grain grain temperature inside silo, moreover it is possible to judge whether Grain Quantity is abnormal in silo according to adjacent two days grain temperature change rate.

Description

A kind of AD strategy process suitable for grain storage supervision
Technical field
The present invention relates to grain reserves technical fields, and more particularly, the present invention relates to a kind of AD suitable for grain storage supervision Strategy process.
Background technique
In China, stored grain system carries China's grain strategic materials reserve, and plays the important task for adjusting market, and Ensure the important foundation of national stability, and the supervision and verification to grain material object are an important rings during grain reserves Section.As China's Agricultural Level continues to develop, the market supply and increase in demand, the scale of grain storage constantly expand.But with The increase of grain depot quantity, grain depot management is not good at, fraud is increasing.Grain depot mainly pass through disengaging storehouse charge book and Artificial investigation is to judge silo inventory information, and this mode takes time and effort, while the carelessness in management work will cause to inventory The false judgment of information is primarily present following defect:
(1) it needs to increase new hardware in silo, at high cost, installation is inconvenient
(2) grain feelings are caused once occurring will cause the loss of key message the problems such as power failure, plant maintenance in use process The interruption of monitoring
(3) there is artificial pause, destroy supervisory systems, allow invalidation of government's supervising phenomenon.
Therefore, grain temperature changing rule is studied, improving grain storage regulation technique has Ensuring Food Safety positive meaning.
Summary of the invention
It, can be according to annual temperature the purpose of the present invention is having designed and developed a kind of AD strategy process suitable for grain storage supervision Grain grain temperature changes amplitude to judge whether grain temperature is abnormal in silo inside period and silo.
The present invention can also judge the temperature rise or temperature of grain in silo according to the differential variation rate of grain grain temperature inside silo Whether drop is abnormal.
The present invention can also judge whether Grain Quantity is abnormal in silo according to adjacent two days grain temperature change rate.
Technical solution provided by the invention are as follows:
A kind of AD strategy process suitable for grain storage supervision, includes the following steps:
Step 1: temperature average value and grain depot location universe in grain depot location one day are determined according to silo geography information Temperature threshold value:
AT=C1+C2×WD+C3×JD+C4×HD+(C5+C6×WD+C7×JD+C8×HD)
×sin(2×3.14/12×MOTH+C9+C10×WD+C11×JD+C12×HD)
ATmax=C1+C5+ (C2+C6) × WD+ (C3+C7) × JD+ (C4+C8) × HD
ATmin=C1+C5+ (C2-C6) × WD+ (C3-C7) × JD+ (C4-C8) × HD
Wherein, AT is temperature average value in grain depot location one day, ATmaxFor grain depot location universe temperature maximum value, ATminFor grain depot location universe temperature minimum value, WD indicates that grain depot location latitude, JD indicate grain depot location longitude, HD table Show that grain depot location height above sea level, MOTH indicate month, C1~C12 is constant;
Step 2: determine each layer grain temperature threshold value of silo:
Wherein, AkminFor the grain temperature minimum value of silo kth layer, AkmaxFor the grain temperature maximum value of silo kth layer, hkFor silo Kth layer and grain face distance, H are grain face height in silo;
Step 3: according to the grain temperature in sampling period acquisition silo, constituting three-dimensional grain temperature matrixWherein,It is temperature measuring point in silo (i, j, k) in the grain temperature of the d days kth layers, i is the position of X-direction in XOY plane, j XOY The position of Y-direction in plane, XOY plane are parallel to horizontal plane, and X-direction is vertical with Y-direction, and k is grain where temperature measuring point (i, j, k) Clad number, d are the time;
Step 4: when grain temperature meets in siloWith One or more of when, grain temperature abnormality in silo.
Preferably, further includes:
The grain temperature in silo is acquired, and determines practical grain temperature change rate:
Wherein, Di,j,k,dFor temperature measuring point in silo (i, j, k) the d days kth layers practical grain temperature change rate,It is temperature measuring point in silo (i, j, k) in the grain temperature of d days kth layers of d- Δ, Δ d is time interval;
When | Di,j,k,d| > | Dai,j,k,d| when, then the temperature rise of grain or temperature drop are abnormal in silo, wherein Dai,j,k,dFor Prediction grain temperature change rate of the temperature measuring point (i, j, k) in the d days kth layers in silo.
Preferably, the prediction grain temperature change rate of every layer of each temperature measuring point is identical in the silo.
Preferably, further includes:
According to the grain temperature change rate of the temperature ecology curve of silo and each layer of silo and the range prediction difference grain layer of grain face:
Wherein, Dak,dFor the prediction grain temperature change rate of the d days kth layers, G is proportionality coefficient,For d+ θkIt Variation of ambient temperature rate, θkNumber of days is lagged for kth layer grain temperature.
Preferably, the variation of ambient temperature rate are as follows:
Wherein, DHadFor the d days variation of ambient temperature rates, ATdFor the d days environment temperatures, ATd-δFor d- δ days environment Temperature, δ are the number of days of temperature interval.
Preferably, the kth layer grain temperature lags number of days are as follows:
Wherein, Q is each layer grain temperature retardation coefficient.
Preferably,
When | Di,j,k,d| > β × | Di,j,k,d-1| when, wherein Di,j,k,d-1It is temperature measuring point in silo (i, j, k) at the d-1 days The practical grain temperature change rate of kth layer, β be adjacent two days grain temperature change rate coefficient, then in silo the d days kth layers Grain Quantity Abnormal or generation ventilation, stifling operation.
Preferably, adjacent two days grain temperature change rate factor beta >=1.5.
Preferably, in the step 3, improper in data or interference data need to be removed by acquiring the grain temperature in silo.
It is of the present invention the utility model has the advantages that
(1) the AD strategy process of the present invention suitable for grain storage supervision, can be according in annual temperature period and silo Portion's grain grain temperature changes amplitude to judge whether grain temperature is abnormal in silo;It can also be according to the differential variation of grain grain temperature inside silo Rate judges whether the temperature rise of grain or temperature drop are abnormal in silo;Grain in silo can also be judged according to adjacent two days grain temperature change rate Whether eclipse number amount is abnormal or occurs abnormal to go out to put in storage with ventilating etc..
(2) hold disconnected property, since grain grain temperature keeps the continuity of trend, the acquisition of grain feelings data interrupt front and back grain feelings data according to So remain high correlation and continuity.Since hardware fault or human factor cause data discontinuous, grain feelings cannot be destroyed The continuity of data supervision is since abnormal operation causes in case of big jump.
(3) Rong Quexing, grain feelings information cause to lack due to hardware fault, can carry out abnormal determination, and pass through numerical value Method completion.
(4) fault-tolerance, grain feelings information lead to mistake due to hardware or communication failure, can carry out abnormal determination, and lead to Numerical method is crossed to be modified.
Detailed description of the invention
Fig. 1 is the flow chart of the AD strategy process of the present invention suitable for grain storage supervision.
Fig. 2 is barn temperature measuring point arrangement schematic diagram described in the embodiment of the present invention 1.
Fig. 3 is the annual ecological temperature profile in silo location described in the embodiment of the present invention 1.
Fig. 4 is the grain temperature change curve of second layer temperature measuring point (3,3,2) described in the embodiment of the present invention 1.
Fig. 5 is barn temperature measuring point arrangement schematic diagram described in the embodiment of the present invention 2.
Fig. 6 is the grain temperature change curve of second layer temperature measuring point (7,3,2) described in the embodiment of the present invention 2.
Specific embodiment
Present invention will be described in further detail below with reference to the accompanying drawings, to enable those skilled in the art referring to specification text Word can be implemented accordingly.
(for AD strategy, intend from grain temperature as shown in Figure 1, the present invention provides a kind of AD strategy process suitable for grain storage supervision Amplitude (grain temperature amplitude of variation, A), the grain temperature difference divide change rate (D) to be analyzed), include the following steps:
Step 1: temperature average value and grain depot location universe in grain depot location one day are determined according to silo geography information Temperature threshold value:
AT=C1+C2×WD+C3×JD+C4×HD+(C5+C6×WD+C7×JD+C8×HD)
×sin(2×3.14/12×MOTH+C9+C10×WD+C11×JD+C12×HD)
ATmax=C1+C5+ (C2+C6) × WD+ (C3+C7) × JD+ (C4+C8) × HD
ATmin=C1+C5+ (C2-C6) × WD+ (C3-C7) × JD+ (C4-C8) × HD
Wherein, AT is temperature average value in grain depot location one day, ATmaxFor grain depot location universe temperature maximum value, ATminFor grain depot location universe temperature minimum value, WD indicates that grain depot location latitude, JD indicate grain depot location longitude, HD table Show that grain depot location height above sea level, MOTH indicate month, C1~C12 is constant;
Step 2: determine each layer grain temperature threshold value of silo:
Wherein, AkminFor the grain temperature minimum value of silo kth layer, AkmaxFor the grain temperature maximum value of silo kth layer, hkFor silo Kth layer and grain face distance, H are grain face height in silo;
Step 3: according to the grain temperature in sampling period acquisition silo, and removing improper in data or interference data, constitute Three-dimensional grain temperature matrixWherein,Grain temperature for temperature measuring point in silo (i, j, k) in the d days kth layers, i For the position of X-direction in XOY plane, j is the position of Y-direction in XOY plane, and XOY plane is parallel to horizontal plane, X-direction and the side Y To vertical, k is the silo number of plies where temperature measuring point (i, j, k), and d is the time;
In the present embodiment, three-dimensional system of coordinate is established, horizontal warehouse is constituted using a base angle of silo as coordinate system original temperature measuring point Three side of base angle is respectively X, Y, Z axis positive direction, and X, Y, Z axis positive direction temperature measuring point number is respectively n, m and l.Tower silo is with silo Interior outermost cable for measuring temperature forms cuboid, using a certain base angle of cuboid as coordinate system original temperature measuring point, constitutes three side of base angle difference For X, Y, Z axis positive direction, X, Y, Z axis positive direction temperature measuring point number is respectively n, m and l, i.e., had both been suitble to horizontal warehouse or had been suitble to cylinder Storehouse.
Step 4: when grain temperature meets in siloWithOne or more of when, grain temperature abnormality in silo.
I.e. in annual periodicity variation, a sinusoidal variation rule is presented in ambient temperature curve.During grain storage, silo The interior mild grain depot location atmospheric temperature of grain grain is similar in annual cycles inner curve, and a similar sinusoidal week is also presented Phase property variation tendency, i.e., amplitude (grain temperature value) variation of grain temperature temperature measuring point should be in annual ambient temperature curve maximum value and minimum in silo Between value.
Step 5: the grain temperature in acquisition silo, and determine practical grain temperature change rate:
Wherein, Di,j,k,dFor temperature measuring point in silo (i, j, k) the d days kth layers practical grain temperature change rate,It is temperature measuring point in silo (i, j, k) in the grain temperature of d days kth layers of d- Δ, Δ d is time interval;
The prediction grain temperature change rate of every layer of each temperature measuring point is identical in the silo, i.e. Dai,j,k,d=Dak,d, wherein Dak,d For the prediction grain temperature change rate of the d days kth layers.
According to the grain temperature change rate of the temperature ecology curve of silo and each layer of silo and the range prediction difference grain layer of grain face:
Wherein, G is proportionality coefficient,For d+ θkIt variation of ambient temperature rate, θkIt is stagnant for kth layer grain temperature Number of days afterwards.
The variation of ambient temperature rate are as follows:
Wherein, DHadFor the d days variation of ambient temperature rates, ATdFor the d days environment temperatures, ATd-δFor d- δ days environment Temperature, δ are the number of days of temperature interval.
The kth layer grain temperature lags number of days are as follows:
Wherein, Q is each layer grain temperature retardation coefficient, is empirical value.
When | Di,j,k,d| > | Dai,j,k,d| when, wherein Dai,j,k,dIt is temperature measuring point in silo (i, j, k) in the d days kth layers Prediction grain temperature change rate, then the temperature rise of grain or temperature drop are abnormal in silo.
Step 6: when | Di,j,k,d| > β × | Di,j,k,d-1| when, wherein Di,j,k,d-1Exist for temperature measuring point in silo (i, j, k) The practical grain temperature change rate of the d-1 days kth layers, β be adjacent two days grain temperature change rate coefficient, then in silo the d days kth layers grain Eclipse number amount is abnormal or ventilation, stifling operation occurs.In the present embodiment, adjacent two days grain temperature change rate factor beta >=1.5.
Design parameter determines as shown in Table 1.
One design parameter decision table of table
Step 7: comprehensive judgement being carried out to analyzed grain temperature temperature measuring point, determines whether the exception detected belongs to improper behaviour Grain temperature abnormality caused by making, if the exception of abnormal operation, then recording exceptional type and date.The ruler of the composite record silo Silos essential information and the grain storages such as very little structure, location, longitude and latitude, height above sea level, stored up crop, time of putting in storage, export time Number, the proportion, position of journey China Oil and Food Import and Export Corporation temperature abnormality type, and its specific abnormal date occurs, eventually form anomaly analysis As a result report file.
Embodiment 1 (by taking horizontal warehouse as an example)
S1. the Minitype granary established using laboratory carries out storage test in Changchun City Chaoyang District high hill village.Grain Storehouse is having a size of 3.8m × 3.8m × 6m, grain face height 4.5m in storehouse, and grain reservation is corn, and kind is first jade 33, water when putting in storage It is divided into 15%.Arrange 25 cable for measuring temperature in storehouse, it is horizontal, indulge each 5,0.8 meter of cable for measuring temperature spacing, away from bulkhead 0.3m, Mei Ge electricity 4 temperature measuring points are arranged on cable, first layer temperature measuring point is away from orlop 0.5m (close to orlop), and the 4th layer away from grain face 0.4m, two layers of centre Temperature measuring point is evenly arranged, as shown in Figure 2.Silo was aerated on June 30th, 2017, was locally generated heat in grain heap after ventilation Area, hot zone reduce after being gradually increased, and disappear in August 4th or so.
S2. grain feelings data are called, the acquisition time of grain feelings data is on August 15th, 11 days 1 May in 2017.
S3. grain depot location is Changchun City, and climate model is as follows:
AT=C1+C2×WD+C3×JD+C4×HD+(C5+C6×WD+C7×JD+C8×HD)
×sin(2×3.14/12×MOTH+C9+C10×WD+C11×JD+C12×HD)
The maximum value AT of universe temperaturemax, minimum value ATminCalculation formula are as follows:
ATmax=C1+C5+ (C2+C6) × WD+ (C3+C7) × JD+ (C4+C8) × HD
ATmin=C1+C5+ (C2-C6) × WD+ (C3-C7) × JD+ (C4-C8) × HD
AT indicates that temperature on average in one day, WD indicate that latitude, JD indicate that longitude, HD indicate that height above sea level, MOTH indicate in model Month, C1~C12 are constant.C1=444.2059, C2=-6.916006, C3=-0.581316, C4=-0.002506, C5= 114.0762, C6=-6.199186, C7=-0.333406, C8=0.000456, C9=1.240849, C10=0.003992, C11 =-0.002969, C12=-1.14 × 10-6, grain depot location latitude is 43.55 °, and longitude is 125.29 °, and height above sea level is 209 meters, The annual ecological temperature curve of generation on annual January 1 to December 31 is as shown in Figure 3.
Determine that grain temperature AD strategy supervision threshold value table two is as follows using universe climate model:
Two universe climate model of table determines that grain temperature AD strategy supervises threshold value table
S4. removal interference and improper data.
S5. data are recombinated according to silo inner sensor distribution situation, constitutes three-dimensional grain temperature matrix or array.
Using silo direction northwest bottom top temperature measuring point as coordinate original temperature measuring point, it is straight up Z axis forward direction, is eastwards X-axis Forward direction, it is positive for Y-axis southwards.According to above-mentioned coordinate system, sensor temperature measuring point XYZ data three-dimensional, matrix X-direction 5 surveys are established Warm spot, 5 temperature measuring points of Y-direction, 4 temperature measuring points of Z-direction, each layer temperature measuring point of silo can indicate respectively0 < i≤5,0 < j≤5, d indicates date, i.e., in the May, 2017 since grain storage On August 15th, 11 days 1.
S6. the judgement of A strategy is carried out to every layer of Z axis positive direction each temperature measuring point grain temperature data respectively and D strategy determines.
S7. comprehensive judgement is carried out to analyzed grain temperature temperature measuring point, whether the exception for determining that S6 is detected belongs to improper behaviour Grain temperature abnormality caused by making, if the exception of abnormal operation, then recording exceptional type and date.Respectively to every layer of silo each grain Warm temperature measuring point is analyzed in the grain temperature variation of time series.Grain temperature change curve is generated such as to second layer temperature measuring point (3,3,2) Shown in Fig. 4.
The region F grain temperature abnormality is judged by A strategy, method particularly includes: obtaining grain temperature maximum in the region F is 34 DEG C, this When be July, grain temperature > ATmax=25 DEG C, the temperature measuring point grain temperature abnormality known to inquiry threshold value table.In G temperature measuring point decline change rate drop It is low.
The grain temperature abnormality temperature rise of the region E is judged by D strategy, method particularly includes:
Known by step 1 barn structure, the second layer is away from grain face h2For 2.8m, grain face height is 4.5m, and the Q time delays experience system Number is selected as 25 days, and the 2nd layer of grain temperature lag number of days θ is calculated2It is 15 days.
In formula: Q is each layer grain temperature retardation coefficient, this silo empirical value is 25 days.
According to the temperature ecology curve of the silo and each layer of silo away from grain relation of plane, the grain temperature prediction of different grain layers is determined Change rate Da2,d, calculation formula is
In formula, h2Away from grain face height be 2.8m for this 2 layers, H is grain face height 4.5m, G be proportionality coefficient (with where grain depot Ground, storehouse type etc. are information-related, and 0.95) empirical value is taken as, Da2,dFor the d days 2 layers of grain temperature change rates,For d- θ2It Variation of ambient temperature rate, θ2Number of days is lagged for the 2nd layer of grain temperature, is 15 days.A-quadrant is on June 30 or so, then seeks gas on June 15 The change rate DHa of temperatured-15=0.93, then Da2,d=0.57.
The practical grain temperature change rate D3,3,2,dAre as follows: according to the current 2 layers of grain temperature data of reading, the region E grain temperature change rate meter Calculate formula are as follows:
D3,3,2,d=(28-18)/15
The practical grain temperature change rate D in the region E is calculated according to formula3,3,2,d=0.67, compare D3,3,2,d> DHa2,d, inquire threshold value The region E can determine that the temperature measuring point overheating by D strategy known to table two.
The region G is judged according to AD strategy, although the judgement of D strategy is normal, A strategy judges that the temperature measuring point is still abnormal.
Thus temperature measuring point grain temperature variation judgement is completed, and successively judges that barn temperature measuring orders grain using plane as class using AD strategy Whether grain temperature is abnormal.The date and position that recording exceptional occurs, while calculating the shared survey of temperature measuring point of abnormal heating, overheating The percentage of warm spot sum.
S8. it is analyzed according to step S7, dimensional structure, location, longitude and latitude, height above sea level, the locating work of the composite record silo The number of the silos essential informations such as object, time of putting in storage, export time and grain storage process China Oil and Food Import and Export Corporation temperature abnormality type, proportion, Position, and the date being specifically abnormal, and form anomaly analysis report file.
Embodiment 2 (horizontal warehouse)
No. 13 silos of Heilongjiang Province's grain depot of grain feelings data from third Grain Reserve Ecology area of China.The storehouse is horizontal warehouse, Stored up crop is corn, and corn moisture is about 13.8% when putting in storage, and puts in storage and finishes on May 6th, 2014,2 sunrise April in 2016 Storehouse.The silo east-west direction long 47.5m, the wide 26m of North and South direction, the high 8m in warehouse, the high 6m of grain face.Built-in cable for measuring temperature 78, east West is to 13, and spacing 3.75m, two sides are from metope 1.25m;North and South direction 6, spacing 5m, two sides are from metope 0.5m.4th layer Temperature measuring point is away from grain face height 0.4m, and first layer temperature measuring point is away from bottom surface 0.5m, and temperature measuring point is uniformly distributed on every cable for measuring temperature, is spaced About 1.7m.Barn temperature measuring point arrangement schematic diagram is as shown in Figure 5.
Second layer temperature measuring point (7,3,2) are selected to generate grain temperature change curve as shown in Figure 6.
The region U grain temperature abnormality is judged by A strategy, method particularly includes: obtaining grain temperature maximum value in the region F is about -15 DEG C, Minimum value is about -17 DEG C, regional environment temperature ATmaxIt is -5 DEG C, ATminIt is -13 DEG C, is at this time -2 months December, grain temperature is maximum ValueKnow that the temperature measuring point grain temperature is too low.It is inquired grain depot historical operation record and learns the grain on December 10 Library is aerated, and cold wind is passed through grain depot, leads to grain temperature abnormality.
AD strategy process of the present invention suitable for grain storage supervision, can be according to grain inside annual temperature period and silo Grain temperature changes amplitude to judge whether grain temperature is abnormal in silo;It can also be sentenced according to the differential variation rate of grain grain temperature inside silo Whether the run out of grain temperature rise of grain or temperature drop in storehouse be abnormal;Grain number in silo can also be judged according to adjacent two days grain temperature change rate Whether amount is abnormal or occurs abnormal to go out to put in storage with ventilating etc..The judgment method is simple to operation, passes through grain grain temperature Variation characteristic and foundation silo where region air temperature model, the normal grain temperature threshold value and grain temperature for different geographical that you can get it become Rate threshold value, to carry out abnormal judgement and abnormity early warning to grain temperature situation of change.Compared to adding hardware, this method in silo With efficient, at low cost, data acquisition is simple, is not easy to crack, have it is fault-tolerant, hold and lack, hold the special temperature measuring point such as disconnected.
Although the embodiments of the present invention have been disclosed as above, but its is not only in the description and the implementation listed With it can be fully applied to various fields suitable for the present invention, for those skilled in the art, can be easily Realize other modification, therefore without departing from the general concept defined in the claims and the equivalent scope, the present invention is simultaneously unlimited In specific details and legend shown and described herein.

Claims (9)

1. a kind of AD strategy process suitable for grain storage supervision, which comprises the steps of:
Step 1: temperature average value and grain depot location universe temperature in grain depot location one day are determined according to silo geography information Threshold value:
AT=C1+C2×WD+C3×JD+C4×HD+(C5+C6×WD+C7×JD+C8×HD)×sin(2×3.14/12×MOTH+ C9+C10×WD+C11×JD+C12×HD)
ATmax=C1+C5+ (C2+C6) × WD+ (C3+C7) × JD+ (C4+C8) × HD
ATmin=C1+C5+ (C2-C6) × WD+ (C3-C7) × JD+ (C4-C8) × HD
Wherein, AT is temperature average value in grain depot location one day, ATmaxFor grain depot location universe temperature maximum value, ATminFor Grain depot location universe temperature minimum value, WD indicate that grain depot location latitude, JD indicate that grain depot location longitude, HD indicate grain depot Location height above sea level, MOTH indicate month, and C1~C12 is constant;
Step 2: determine each layer grain temperature threshold value of silo:
Wherein, AkminFor the grain temperature minimum value of silo kth layer, AkmaxFor the grain temperature maximum value of silo kth layer, hkFor silo kth layer With grain face distance, H is grain face height in silo;
Step 3: according to the grain temperature in sampling period acquisition silo, constituting three-dimensional grain temperature matrixWherein,It is temperature measuring point in silo (i, j, k) in the grain temperature of the d days kth layers, i is the position of X-direction in XOY plane, j XOY The position of Y-direction in plane, XOY plane are parallel to horizontal plane, and X-direction is vertical with Y-direction, and k is grain where temperature measuring point (i, j, k) Clad number, d are the time;
Step 4: when grain temperature meets in silo With One or more of when, grain temperature abnormality in silo.
2. the AD strategy process suitable for grain storage supervision as described in claim 1, which is characterized in that further include:
The grain temperature in silo is acquired, and determines practical grain temperature change rate:
Wherein, Di,j,k,dFor temperature measuring point in silo (i, j, k) the d days kth layers practical grain temperature change rate,For grain For temperature measuring point (i, j, k) in the grain temperature of d days kth layers of d- Δ, Δ d is time interval in storehouse;
When | Di,j,k,d| > | Dai,j,k,d| when, then the temperature rise of grain or temperature drop are abnormal in silo, wherein Dai,j,k,dFor silo Prediction grain temperature change rate of the interior temperature measuring point (i, j, k) in the d days kth layers.
3. the AD strategy process suitable for grain storage supervision as claimed in claim 2, which is characterized in that every layer of each survey in the silo The prediction grain temperature change rate of warm spot is identical.
4. the AD strategy process suitable for grain storage supervision as claimed in claim 3, which is characterized in that further include:
According to the grain temperature change rate of the temperature ecology curve of silo and each layer of silo and the range prediction difference grain layer of grain face:
Wherein, Dak,dFor the prediction grain temperature change rate of the d days kth layers, G is proportionality coefficient,For d+ θkIt ring Border rate of temperature change, θkNumber of days is lagged for kth layer grain temperature.
5. the AD strategy process suitable for grain storage supervision as claimed in claim 4, which is characterized in that the variation of ambient temperature rate Are as follows:
Wherein, DHadFor the d days variation of ambient temperature rates, ATdFor the d days environment temperatures, ATd-δFor d- δ days environment temperature, δ is the number of days of temperature interval.
6. the AD strategy process suitable for grain storage supervision as claimed in claim 5, which is characterized in that the kth layer grain temperature lag Number of days are as follows:
Wherein, Q is each layer grain temperature retardation coefficient.
7. the AD strategy process for being suitable for grain storage supervision as described in any one of claim 2-6, which is characterized in that
When | Di,j,k,d| > β × | Di,j,k,d-1| when, wherein Di,j,k,d-1It is temperature measuring point in silo (i, j, k) in the d-1 days kth layers Practical grain temperature change rate, β is that adjacent two days grain temperature change rate coefficient, then in silo the Grain Quantity of the d days kth layers extremely or Ventilation, stifling operation occur for person.
8. the AD strategy process suitable for grain storage supervision as claimed in claim 7, which is characterized in that adjacent two days grain temperature become Rate factor beta >=1.5.
9. the AD strategy process for being suitable for grain storage supervision as described in claim 1-6 or 8, which is characterized in that in the step 3, Grain temperature in acquisition silo need to remove improper in data or interference data.
CN201810566299.5A 2018-06-05 2018-06-05 A kind of AD strategy process suitable for grain storage supervision Pending CN108981960A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810566299.5A CN108981960A (en) 2018-06-05 2018-06-05 A kind of AD strategy process suitable for grain storage supervision

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810566299.5A CN108981960A (en) 2018-06-05 2018-06-05 A kind of AD strategy process suitable for grain storage supervision

Publications (1)

Publication Number Publication Date
CN108981960A true CN108981960A (en) 2018-12-11

Family

ID=64540040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810566299.5A Pending CN108981960A (en) 2018-06-05 2018-06-05 A kind of AD strategy process suitable for grain storage supervision

Country Status (1)

Country Link
CN (1) CN108981960A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110608816A (en) * 2019-09-04 2019-12-24 安徽省粮油信息中心(安徽粮食批发交易市场管委会) Method for abnormal investigation and result analysis and identification during grain temperature detection
CN110631721A (en) * 2019-09-20 2019-12-31 辽宁省粮食科学研究所 Granary heat insulation judgment method based on grain condition big data
CN110686371A (en) * 2019-10-24 2020-01-14 辽宁省粮食科学研究所 Granary air conditioner automatic temperature control method based on temperature field cloud picture
CN110704512A (en) * 2019-10-22 2020-01-17 吉林大学 Granary ventilation time interval judgment method based on historical grain situation data
CN117666447A (en) * 2024-01-31 2024-03-08 泰州市衡顺电控科技有限公司 Digital intelligent granary data supervision system and method of multi-control scheme

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110608816A (en) * 2019-09-04 2019-12-24 安徽省粮油信息中心(安徽粮食批发交易市场管委会) Method for abnormal investigation and result analysis and identification during grain temperature detection
CN110631721A (en) * 2019-09-20 2019-12-31 辽宁省粮食科学研究所 Granary heat insulation judgment method based on grain condition big data
CN110704512A (en) * 2019-10-22 2020-01-17 吉林大学 Granary ventilation time interval judgment method based on historical grain situation data
CN110704512B (en) * 2019-10-22 2022-05-24 吉林大学 Granary ventilation time interval judgment method based on historical grain situation data
CN110686371A (en) * 2019-10-24 2020-01-14 辽宁省粮食科学研究所 Granary air conditioner automatic temperature control method based on temperature field cloud picture
CN110686371B (en) * 2019-10-24 2021-03-26 辽宁省粮食科学研究所 Granary air conditioner automatic temperature control method based on temperature field cloud picture
CN117666447A (en) * 2024-01-31 2024-03-08 泰州市衡顺电控科技有限公司 Digital intelligent granary data supervision system and method of multi-control scheme
CN117666447B (en) * 2024-01-31 2024-04-12 泰州市衡顺电控科技有限公司 Digital intelligent granary data supervision system and method of multi-control scheme

Similar Documents

Publication Publication Date Title
CN108981960A (en) A kind of AD strategy process suitable for grain storage supervision
Luo et al. Development of an IoT-based big data platform for day-ahead prediction of building heating and cooling demands
Buffa et al. Advanced control and fault detection strategies for district heating and cooling systems—A review
Gunay et al. Inverse blackbox modeling of the heating and cooling load in office buildings
CN105958940B (en) The performance monitor system and its monitoring method of solar cell module
Shi et al. Prediction of indoor temperature and relative humidity based on cloud database by using an improved BP neural network in Chongqing
US11048247B2 (en) Building management system to detect anomalousness with temporal profile
Crevoisier et al. Regional US carbon sinks from three-dimensional atmospheric CO2 sampling
Esrafilian-Najafabadi et al. Occupancy-based HVAC control using deep learning algorithms for estimating online preconditioning time in residential buildings
Li et al. Response of energy consumption for building heating to climatic change and variability in Tianjin City, China
CN108764802A (en) A kind of ABC strategy process suitable for the supervision of grain storage number
Chung Estimating solar insolation and power generation of photovoltaic systems using previous day weather data
Kalogirou et al. Development of a neural network-based fault diagnostic system for solar thermal applications
Wang et al. Responses of extreme high temperatures to urbanization in the Beijing–Tianjin–Hebei urban agglomeration in the context of a changing climate
Bersani et al. Model predictive control versus traditional relay control in a high energy efficiency greenhouse
CN105373849A (en) Stepped temperature prediction method for photovoltaic battery assembly
Zsiborács et al. Assessing shading losses of photovoltaic power plants based on string data
CN108764749A (en) A kind of 3I strategy process for preferably depositing evaluation suitable for grain storage
Succetti et al. Challenges and perspectives of smart grid systems in Islands: a real case study
Miller et al. A seed dataset for a public, temporal data repository for energy informatics research on commercial building performance
JP5658597B2 (en) Energy-saving behavior evaluation system
Xie et al. Simulation of water temperature in paddy fields by a heat balance model using plant growth status parameter with interpolated weather data from weather stations
Gao et al. Evolutionary characteristics of the atmospheric circulations for frequent and infrequent dust storm springs in northern China and the detection of potential future seasonal forecast signals
Muhammad Lawan et al. Wind energy assessment and mapping using terrain nonlinear autoregressive neural network (TNARX) and wind station data
Quinting et al. The intensity and motion of hybrid cyclones in the Australian region in a composite potential vorticity framework

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20181211

RJ01 Rejection of invention patent application after publication