CN110443423A - A kind of regional rotation grazing dynamic optimization method based on Task Assignment Model - Google Patents
A kind of regional rotation grazing dynamic optimization method based on Task Assignment Model Download PDFInfo
- Publication number
- CN110443423A CN110443423A CN201910722731.XA CN201910722731A CN110443423A CN 110443423 A CN110443423 A CN 110443423A CN 201910722731 A CN201910722731 A CN 201910722731A CN 110443423 A CN110443423 A CN 110443423A
- Authority
- CN
- China
- Prior art keywords
- cell
- rotation grazing
- livestock
- grazing
- seed stock
- 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
- 238000009304 pastoral farming Methods 0.000 title claims abstract description 426
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000005457 optimization Methods 0.000 title claims abstract description 21
- 244000144980 herd Species 0.000 claims abstract description 196
- 244000144972 livestock Species 0.000 claims abstract description 162
- 241001465754 Metazoa Species 0.000 claims abstract description 45
- 244000144977 poultry Species 0.000 claims description 25
- 241001494479 Pecora Species 0.000 claims description 24
- 239000004459 forage Substances 0.000 claims description 12
- 230000006651 lactation Effects 0.000 claims description 8
- 235000021050 feed intake Nutrition 0.000 claims description 5
- 230000007850 degeneration Effects 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 4
- 238000005192 partition Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000012423 maintenance Methods 0.000 claims description 3
- 244000025254 Cannabis sativa Species 0.000 abstract description 7
- 238000009826 distribution Methods 0.000 abstract description 5
- 230000007613 environmental effect Effects 0.000 abstract description 3
- 238000011161 development Methods 0.000 abstract description 2
- 230000037213 diet Effects 0.000 abstract 1
- 235000005911 diet Nutrition 0.000 abstract 1
- 238000012913 prioritisation Methods 0.000 abstract 1
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 241000287196 Asthenes Species 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 241000283707 Capra Species 0.000 description 1
- 208000034657 Convalescence Diseases 0.000 description 1
- 101100001672 Emericella variicolor andG gene Proteins 0.000 description 1
- 238000009395 breeding Methods 0.000 description 1
- 230000001488 breeding effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 238000010189 synthetic method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Animal Husbandry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Educational Administration (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Housing For Livestock And Birds (AREA)
Abstract
The invention discloses a kind of regional rotation grazing dynamic optimization method based on Task Assignment Model, reasonable regional rotation grazing can guarantee to avoid overgrazing under the premise of livestock normally searches for food.There are Various Complex and dynamic element in regional rotation grazing, the present invention herds number of days by the Types of Grassland in pasture, rotation grazing cycle, cell, seed stock herd carries out the division of rotation grazing cell.Type, the livestock daily diet, rotation grazing cell day grass yield, animal number for considering livestock, establish rotation grazing district dynamic Optimized model based on Task Assignment Model.The matching that different seed stock herds Yu rotation grazing cell are carried out by normalized differential vegetation index and animal number, realizes the distribution of each appointment stage rotation grazing cell, to obtain the regional rotation grazing prioritization scheme within a rotation grazing cycle.It is that regional rotation grazing is planned stage by stage according to progress with dynamic animal number and seed stock herd, under conditions of can satisfy ecologically environmental restriction, grazing time using and space configuration on scientific carry out grazing, improve grassland resources utilization efficiency, promote animal husbandry development.
Description
Technical field
The present invention relates to technical field of animal husbandry, more particularly say, are related to a kind of stroke for considering biology, ecology and economy
Area's rotation grazing optimization method.
Background technique
Regional rotation grazing (Rotational Grazing) is also referred to as rotation grazing on an area basis, can improve grass poultry imbalance, ecological environment
The situation of deterioration is the Grazing system for efficiently using grassland resources.Regional rotation grazing be under the premise of Forage-Livestock Balance, it is planned
Seasonal grazing grassland is divided into several rotation grazing cells (Rotational Grazing Plot), herdsman (Herdsman) is allowed to raise
(feed intake is with sheep unit, 1 weight 45Kg, day to the field-to-field feeding feeding in a certain order of feeding livestock (livestock)
Consume the adult sheep of 1.8Kg standard hay, or with this other comparable domestic animal), the grazing system that utilizes of circulation, it is ensured that
There is certain interval recovery time in the pasture of each rotation grazing cell, avoids and overgrazes and alleviate livestock to grassland
It tramples on and destroys, while can also guarantee that livestock obtains enough herbage within the grazing season, be a kind of putting for adjusting Forage-Livestock Balance
Herd management method.
Assortment is carried out to the livestock (Livestock) of herdsman (Herdsman) raising, with reference to " the People's Republic of China (PRC)
Agricultural industry criteria " NY/T 635-2015, " poultry kind " classification, Ke Yifen in the table 4 of " calculating of the reasonable animal number of natural meadow "
For sheep, goat, ox, buffalo etc..
(Rotational Grazing planning) is planned for unreasonable regional rotation grazing, to livestock
(Livestock) for, it can not meet the daily feeding demand of livestock；For pasture (Pasture), it be easy to cause grassland raw
State system degradation, the vicious circle for causing reasonable animal number (Carrying Capacity) to further decrease.Rationally animal number is
Refer to certain grassland area, a certain using in the period, utilized in conservative grazing (or mowing) and maintain the sustainable production in meadow
Under the premise of, meet the needs of domestic animal normal growth, breeding, production, the most animal numbers that can be carried.Reasonable animal number is again
Claim Carrying.
For in a traditional way, the region that herdsman herds is to determine exist to a certain extent by herdsman's subjective experience
Grass raises unbalanced problem, destroys the ecology in pasture.
Summary of the invention
The invention discloses a kind of regional rotation grazing dynamic optimization method based on Task Assignment Model, for two class seed stock herds one
The rotation grazing cell divided in a rotation grazing cycle is assigned in order.Firstly, on the basis of considering Forage-Livestock Balance, according to grassland
Type, rotation grazing cycle, cell herd number of days, and the two class seed stock herds divided by livestock to herbage demand, take turns to pasture
Herd the division of cell.Secondly, determining the appointment rank in rotation grazing cycle according to two class seed stock herds of division and the number of rotation grazing cell
Section.Then, the feed intake of unified conversion herdsman's variety classes livestock.Finally, establishing rotation grazing district dynamic based on Task Assignment Model
Optimized model is suitable for the rotation grazing cell herded by normalized differential vegetation index (NDVI) screening, according to animal number height, is needed to herbage
The high seed stock herd of the amount of asking preferentially is assigned and each rotation grazing cell is only assigned primary principle within a rotation grazing cycle, every
A appointment stage carries out the matching of seed stock herd and rotation grazing cell, to obtain the regional rotation grazing optimization side within a rotation grazing cycle
Case.It is that can satisfy ecology according to the planning for successively carrying out regional rotation grazing by the appointment stage with dynamic animal number and seed stock herd
Under conditions of environmental constraints, grazing time using and space configuration on scientific carry out grazing, improve grassland resources
Utilization efficiency promotes animal husbandry development.
A kind of regional rotation grazing dynamic optimization method based on Task Assignment Model of the invention, specifically there is the following steps:
Step 1: configuration rotation grazing cell；
Step 11, it according to number of days is herded rotation grazing cycle and partition to any one pasture progress small Division of rotation grazing, obtains
With the M block rotation grazing cell for dividing region；Execute step 12；
RE={ E is expressed as using aggregate form to the rotation grazing cell that division obtains1,E2,…,Ev,…EM}；
The sum M of rotation grazing cell are as follows:
Step 12, preferential treatment seed stock herd RA and conventional seed stock herd RB is respectively divided in the livestock that herdsman raises, executes step
13；
Preferential treatment seed stock herd RA refers to the livestock in growth stage or lactation period；
Conventional seed stock herd RB refers to the livestock in maintenance phase or non-lactation period；
Step 13, the appointment stage g in rotation grazing cycle is determined by rotation grazing cell total number M；Execute step 31；
The preferential treatment seed stock herd RA and routine seed stock herd RB divided according to step 12 is assigned in M rotation grazing cell respectively,
Therefore, the appointment stage sum of rotation grazing cell included in rotation grazing cycle is
Step 2: initialization livestock number；
Step 21, after step 11 processing, fence design is carried out to pasture, the phase of herding is herded beginning, and executes step 22；
Step 22, the livestock number of standard sheep unit is obtained；Execute step 31；
Any one kind livestock of preferential treatment seed stock herd RA is denoted asIt corresponds to describedSheep unit conversion factor
It is denoted asIt corresponds to describedFeeding quantity be denoted as
Any one kind livestock of conventional seed stock herd RB is denoted asIt corresponds to describedSheep unit conversion factor
It is denoted asIt corresponds to describedFeeding quantity be denoted as
For the livestock head number of unified different seed stock herds, rolled over using the quantity of different cultivars multiplied by the sheep unit of different cultivars
Coefficient is calculated, the number of the standard sheep unit of each seed stock herd is obtained；
For the livestock number of the standard sheep unit of preferential treatment seed stock herd RA, it is denoted as RAT, and
For the livestock number of the standard sheep unit of conventional seed stock herd RB, it is denoted as RBT, and
Step 3: the dynamic regional rotation grazing planning based on Task Assignment Model；
Step 31, the small Division of rotation grazing is carried out using normalized differential vegetation index, executes step 32；
From rotation grazing set of cells RE={ E1,E2,…,Ev,…EMIn obtain normalized differential vegetation index be higher than 0.4 rotation grazing
Cell is denoted as qualified rotation grazing set of cells
From rotation grazing set of cells RE={ E1,E2,…,Ev,…EMIn obtain normalized differential vegetation index less than or equal to 0.4
Rotation grazing cell is denoted as degeneration rotation grazing set of cells
In the present invention, the RENDVI_1 ∈ RE, the RENDVI_2 ∈ RE；
Step 32, appointment divided stages are carried out in a rotation grazing cycle, execute step 33；
Step 33, whether seed stock herd is assigned to rotation grazing cell；
Judge assigning stage g, whether seed stock herd is assigned to rotation grazing cell Ev, it is according to rotation grazing cell EvAnimal number
Whether it is more than or equal to the quantity of seed stock herd, indicates are as follows:
IfStage g then is being assigned, giving preferential treatment to seed stock herd RA can be higher than in normalized differential vegetation index
0.4 rotation grazing cellFeeding executes step 34；
If allStage g then is being assigned, preferential treatment seed stock herd RA cannot be in normalized differential vegetation index
Rotation grazing cell higher than 0.4Feeding executes step 35；
IfIt can then be higher than in normalized differential vegetation index in appointment stage g, conventional seed stock herd RB
0.4 rotation grazing cellFeeding executes step 34；
If allStage g then is being assigned, conventional seed stock herd RB cannot be in normalized differential vegetation index
Rotation grazing cell higher than 0.4Feeding executes step 35；
Step 34, result is assigned to obtain；
Any one assigns the appointment result in stage to have: animal number highest rotation grazing cellIt is high with animal number time
Rotation grazing cell
Step 35, judge whether rotation grazing cell is searched for food by livestock；
Assigning the stage to assign after g=1 terminates, and whether record rotation grazing cell is searched for food by livestock；To the wheel that do not searched for food by livestock
Cell is herded, step 32 is repeated to step 34, carries out the appointment that g=2 assigns the stage；Assigning the stage to assign in g=2 terminates,
Whether record rotation grazing cell is searched for food by livestock, to the rotation grazing cell that do not searched for food by livestock, repeats step 32 to step 34, into
The appointment in row g=3 appointment stage；Until theA appointment stage, which assigns, completes, M rotation grazing cell whole seed stock herds
It is searched for food, dynamic regional rotation grazing planning is completed, and step terminates；
If the animal number of all rotation grazing cells is all unsatisfactory for seed stock herd animal number in the stage assigned, merge buffer area
Supplementary feeding is herded or carried out to L and the highest rotation grazing cell of animal number；
It is denoted as by the rotation grazing cell that livestock searches for food:Table
Show that first rotation grazing cell is searched for food by livestock；Indicate that second rotation grazing cell is searched for food by livestock；Indicate v-th of rotation grazing
Cell is searched for food by livestock；Indicate that j-th of rotation grazing cell is searched for food by livestock；Indicate that m-th rotation grazing cell is adopted by livestock
Food；
It is not denoted as by the rotation grazing cell that livestock searches for food:
Indicate that first rotation grazing cell is not searched for food by livestock；Indicate that second rotation grazing cell is not searched for food by livestock；It indicates v-th
Rotation grazing cell is not searched for food by livestock；Indicate that j-th of rotation grazing cell is not searched for food by livestock；Indicate m-th rotation grazing cell not
It is searched for food by livestock.
Regional rotation grazing dynamic optimization method based on Task Assignment Model of the invention, the dynamic optimization of Task Assignment Model are as follows:
It is g-th of appointment stage that v-th of rotation grazing is small under conditions of expression meets normalized differential vegetation index
Area is assigned to preferential treatment seed stock herd RA；
Under conditions of expression meets normalized differential vegetation index, g-th of appointment stage, preferential treatment seed stock herd RA is being taken turns
Herd cell EvIn for forage total amount；
It is g-th of appointment stage that j-th of rotation grazing is small under conditions of expression meets normalized differential vegetation index
Area is assigned to conventional seed stock herd RB；
Under conditions of expression meets normalized differential vegetation index, g-th of appointment stage, conventional seed stock herd RB is being taken turns
Herd cell EjIn for forage total amount；
DgFor the total amount g appointment stage for livestock forage；
The forward direction of livestock productivity is influenced considering livestock feeding, therefore DgIt is maximized are as follows:
Under conditions of expression meets normalized differential vegetation index, g-th of stage, v-th of rotation grazing is small
Area can only be herded for preferential treatment seed stock herd RA；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, seed stock herd is given preferential treatment to
RA can only enter a rotation grazing cell and herd；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, j-th of rotation grazing is small
Area can only be herded for conventional seed stock herd RB；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, conventional seed stock herd
RB can only enter a rotation grazing cell and herd；
IfThen giving preferential treatment to seed stock herd RA can be in rotation grazing cell of the normalized differential vegetation index higher than 0.4Feeding；
IfThen giving preferential treatment to seed stock herd RA cannot be small higher than 0.4 rotation grazing in normalized differential vegetation index
AreaFeeding；
IfThen conventional seed stock herd RB can be in rotation grazing cell of the normalized differential vegetation index higher than 0.4Feeding；
IfThen conventional seed stock herd RB cannot be small higher than 0.4 rotation grazing in normalized differential vegetation index
AreaFeeding；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, gives preferential treatment to seed stock herd RA and arrive
V rotation grazing cell is herded；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, seed stock herd RA is given preferential treatment to
To failing in exam, v rotation grazing cell is herded；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, conventional seed stock herd RB is arrived
Jth rotation grazing cell is herded；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, conventional seed stock herd RB is arrived
Jth rotation grazing cell is not herded.
The present invention is based on the advantages of dynamic regional rotation grazing optimization method of Task Assignment Model to be:
1. the method for the present invention considers the rotation grazing cycle of different Types of Grassland, rotation grazing cell herds number of days, to a pasture
Rotation grazing number of cells and the stage of rotation grazing cycle carry out reasonable disposition.
2. the method for the present invention variety classes livestock is pressed seed stock herd division is carried out to the desirability of herbage on the basis of,
The appointment for dividing seed stock herd to carry out rotation grazing cell, can satisfy demand of the different livestocks to herbage.
3. the method for the present invention is based on Task Assignment Model, it is contemplated that the animal number in rotation grazing cycle in each appointment stage, in animal number
Constraint under rotation grazing cell is assigned, each rotation grazing cell is merely able to be assigned primary within entire rotation grazing cycle, can
Overgrazed is avoided, and each rotation grazing cell includes convalescence, improves grassland resource utilization rate.
Detailed description of the invention
Fig. 1 is the flow chart of the regional rotation grazing dynamic optimization method the present invention is based on Task Assignment Model.
Fig. 2 is the rotation grazing cell schematic diagram in any one pasture.
Fig. 3 is optimization method of the present invention and the comparison diagram that is randomly assigned.
Specific embodiment
Below in conjunction with drawings and examples, the present invention is described in further detail.
Normalized differential vegetation index (NDVI, Normalized Difference Vegetation Index) can be accurately anti-
Reflect vegetation cover situation.It is the annual Chinese vegetation index provided by Resources and environmental sciences data center, the Chinese Academy of Sciences
(NDVI) spatial distribution data collection.Annual Chinese vegetation index (NDVI) spatial distribution data collection is based on continuous time series
SPOT/VEGETATION NDVI satellite remote sensing date, using maximum value synthetic method generate 1998 since annual vegetation refer to
Number data set.The data set effectively reflects vegetative coverage distribution and variation shape of the area in all parts of the country on room and time scale
Condition, to coupling relationship condition monitoring, vegetation resources rationally using and other ecological environment related fieldss research have it is particularly significant
Reference significance.In the present invention, different rotation grazing cell herbage situations can be reacted using normalized differential vegetation index (NDVI).
In the present invention, the rotation grazing cycle, cell herd number of days with reference to People's Republic of China's agricultural industry criteria
" grassland regional rotation grazing technical regulation ", NY/T 1343-2007, the content of the Table A 1 of publication on 04 17th, 2007, in this hair
Table 1 is listed as in bright.
1 different type grassland regional rotation grazing of table designs important technological parameters
| Types of Grassland | Grazing frequency (secondary) | Cell herds number of days (d) | Rotation grazing cycle (d) |
| Warm nature meadow steppe class | 3~4 | 3~5 | 30~40 |
| Warm nature Steppe class | 2~3 | 5~8 | 50~75 |
| Warm nature desert steppe class | 2 | 6~12 | 75~50 |
| High-cold steppe class | 2 | 6~10 | 75~90 |
| Temperate steppe Desert | 2 | 6~10 | 75~90 |
| Lowland meadow class | 3~4 | 3~5 | 30~40 |
| Mountain Meadow class | 3~4 | 3~5 | 30~40 |
| Warm nature thick grass Inertial coupling | 12~9 | 4~6 | 30~40 |
| Alpine meadow class | 2 | 6~12 | 75~80 |
| Hot thick grass, Inertial coupling | 9~12 | 2~4 | 30~40 |
| Warm property thick grass, Inertial coupling | 8~11 | 2~4 | 35~45 |
In the present invention, the sheep unit conversion factor is with reference to People's Republic of China's agricultural industry criteria " natural grass
The calculating of the reasonable animal number in ground ", NY/T 635-2015, the content of the table 4 of publication on May 21st, 2015 arranges in the present invention
It is shown as table 2.
The various adult domestic animals of table 2 are converted into the conversion factor of standard domestic animal unit
In the present invention, any one pasture (Pasture) is divided using network, and it is small to obtain the rotation grazing with region division
Area (Rotational Grazing Plot), as shown in Figure 2.In figure, a pasture is divided into raising area and buffer area, the feeding
There are multiple rotation grazing cells in feeding area.The total block data mesh of the rotation grazing cell is denoted as M.The multiple rotation grazings divided to a pasture are small
Area is expressed as RE={ E using aggregate form1,E2,…,Ev,…,Ej,…EM}；
E1Indicate first rotation grazing cell for belonging to pasture.
E2Indicate second rotation grazing cell for belonging to pasture.
EvIndicate v-th of rotation grazing cell for belonging to pasture.Subscript v indicates the identification number of rotation grazing cell, says for convenience
It is bright, the EvAlso referred to as first any rotation grazing cell.
EMIndicate the last one the rotation grazing cell for belonging to pasture.The sum of subscript M expression rotation grazing cell.
In the present invention, livestock is divided into two class seed stock herds by the demand according to livestock to herbage, and one kind is the domestic animal of lactation period
Poultry or the livestock of growth stage are denoted as preferential treatment seed stock herd RA；Another kind of is the livestock of non-lactation period or the livestock of non-growth stage note
For conventional seed stock herd RB.
In the present invention, the livestock of herdsman (Herdsman) raising generally has multiple kinds (i.e. poultry kind), herdsman's raising
Poultry kind sum be denoted as q, any one kind livestock is denoted as ap, the mark of the kind of subscript p expression herdsman raise livestock.Phase
Corresponding to apConversion factor be denoted as bp, correspond to apFeeding quantity be denoted as cp.The disclosed poultry kind of reference table 2,
The feed intake of each poultry kind is calculated with sheep unit conversion factor.
Each pasture has existing animal number (standing carrying capacity), and existing animal number refers to centainly
The meadow of area utilizes in the period in certain, actually holds feeding standard livestock head number.In the present invention, rotation grazing cell
The appointment stage, which refers to, to be obtained in a rotation grazing cycle according to the total quantity of rotation grazing cell and seed stock herd quantity, referred to as appointment rank
Section；The appointment stage is indicated with letter g.Any one rotation grazing cell EvExisting animal number when assigning stage g, is denoted as
Shown in Figure 1, the present invention is based on the specific steps of the dynamic regional rotation grazing optimization method of Task Assignment Model are as follows:
Step 1: configuration rotation grazing cell；
Step 11, foundation rotation grazing cycle and partition herd number of days and carry out rotation grazing cell to any one pasture (Pasture)
(Rotational Grazing Plot) is divided, and is obtained with the M block rotation grazing cell for dividing region；Execute step 12；
RE={ E is expressed as using aggregate form to the rotation grazing cell that division obtains1,E2,…,Ev,…EM}；
E1Indicate first rotation grazing cell for belonging to pasture；
E2Indicate second rotation grazing cell for belonging to pasture；
EvIndicate any one the rotation grazing cell for belonging to pasture, subscript v indicates the identification number of rotation grazing cell；
EMIndicate the last one the rotation grazing cell for belonging to pasture, subscript M indicates the sum of rotation grazing cell.
In the present invention, rotation grazing cycle, partition herd number of days with reference to People's Republic of China's agricultural industry criteria " grassland
Regional rotation grazing technical regulation ", NY/T 1343-2007 is issued on 04 17th, 2007.That is table 1.
The present invention is allocated rotation grazing cell for two seed stock herds, therefore the sum M of rotation grazing cell are as follows:
X indicates rotation grazing cycle, and Y indicates that cell herds number of days, and L indicates buffer area, is alternative rotation grazing cell.
Such as, the pasture of an existing Alpine meadow class, in corresponding table 1, rotation grazing cycle X=75~80；Cell herds number of days Y
=6~12.It is 1 that buffer area number, which is arranged, and rotation grazing cell herds number of days Y=10, then suitable rotation grazing number of cells M=17.Tool
Body: rotation grazing set of cells isE1Indicate that first rotation grazing for belonging to pasture is small
Area, E2Indicate second rotation grazing cell for belonging to pasture, E3Indicate the third rotation grazing cell for belonging to pasture, E4Expression, which belongs to, herds
4th rotation grazing cell of field, E5Indicate the 5th rotation grazing cell for belonging to pasture, E6Indicate that the 6th rotation grazing for belonging to pasture is small
Area, E7Indicate the 7th rotation grazing cell for belonging to pasture, E8Indicate the 8th rotation grazing cell for belonging to pasture, E9Expression, which belongs to, herds
9th rotation grazing cell of field, E10Indicate the tenth rotation grazing cell for belonging to pasture, E11Indicate the 11st wheel for belonging to pasture
Herd cell, E12Indicate the 12nd rotation grazing cell for belonging to pasture, E13Indicate the 13rd rotation grazing cell for belonging to pasture, E14
Indicate the 14th rotation grazing cell for belonging to pasture, E15Indicate the 15th rotation grazing cell for belonging to pasture, E16Expression, which belongs to, herds
16th rotation grazing cell of field.
In the present invention, due to needing to meet the feeding demand of livestock, when dividing rotation grazing cell, rotation grazing cell number
Amount M cannot be using rounding up, and need to round up calculating.
Step 12, preferential treatment seed stock herd RA and conventional seed stock herd is respectively divided in herdsman (Herdsman) livestock raised
RB, for convenience of the appointment of rotation grazing cell is carried out, preferential treatment seed stock herd RA and the feed intake difference of conventional seed stock herd RB are no more than ± 5%；
Execute step 13；
Preferential treatment seed stock herd RA refers to the livestock in growth stage or lactation period.
Conventional seed stock herd RB refers to the livestock in maintenance phase or non-lactation period.
Step 13, the appointment stage g in rotation grazing cycle is determined by rotation grazing cell total number M；Execute step 31；
The preferential treatment seed stock herd RA and routine seed stock herd RB divided according to step 12 is assigned in M rotation grazing cell respectively,
Therefore, the appointment stage sum of rotation grazing cell included in rotation grazing cycle is
In order to clearly illustrate the appointment stage of rotation grazing cell included in rotation grazing cycle, indicated with g, andG represents the mark that different rotation grazing cells assign the stage.
Such as, Alpine meadow class pasture determines that the appointment stage in rotation grazing cycle is g=by rotation grazing cell total number M=17
8.Specifically: the appointment stage is g=1,2,3,4,5,6,7,8.
Step 2: initialization livestock number；
Step 21, after step 11 processing, fence design is carried out to pasture (Pasture), the phase of herding is herded beginning, and is held
Row step 22；
Step 22, the livestock number of standard sheep unit is obtained；Execute step 31；
Any one kind livestock of preferential treatment seed stock herd RA is denoted asSubscript p indicates the kind of the livestock of herdsman's raising
Mark.It corresponds to describedSheep unit conversion factor be denoted asIt corresponds to describedFeeding quantity be denoted as
Such as, the livestock in the preferential treatment seed stock herd RA of herdsman's raising includes:
Medium-sized sheep 28, is denoted as
Large-scale horse is 12, is denoted asThen kind q=2.
According to table 2, then
Any one kind livestock of conventional seed stock herd RB is denoted asSubscript p indicates the mark of the kind for the livestock that herdsman supports
Know.It corresponds to describedSheep unit conversion factor be denoted asIt corresponds to describedFeeding quantity be denoted as
Such as, the livestock in the conventional seed stock herd RB of herdsman's raising includes:
Medium-sized sheep 20, is denoted as
Large-scale horse is 6, is denoted as
Large-scale ox 5, is denoted asThen kind q=3.
According to table 2, then
In the present invention, for the livestock head number of unified different seed stock herds, using the quantity of different cultivars multiplied by different product
The sheep unit conversion factor of kind, obtains the number of the standard sheep unit of each seed stock herd；
For the livestock number of the standard sheep unit of preferential treatment seed stock herd RA, it is denoted as RAT, and
For the livestock number of the standard sheep unit of conventional seed stock herd RB, it is denoted as RBT, and
Such as, the livestock conversion in the preferential treatment seed stock herd RA of above-mentioned herdsman's raising is the livestock number of standard sheep unit are as follows:
Such as, the livestock conversion in the conventional seed stock herd RB of above-mentioned herdsman's raising is the livestock number of standard sheep unit are as follows:
P indicates the mark of the kind of the livestock of herdsman's raising.
Q indicates the kind sum of the livestock of herdsman's raising.
Indicate the feeding quantity for belonging to the livestock kind of preferential treatment seed stock herd RA.
Indicate the conversion factor for belonging to the livestock kind of preferential treatment seed stock herd RA.
Indicate the feeding quantity for belonging to the livestock kind of conventional seed stock herd RB.
Indicate the conversion factor for belonging to the livestock kind of conventional seed stock herd RB.
Step 3: the dynamic regional rotation grazing planning based on Task Assignment Model；
Step 31, the small Division of rotation grazing is carried out using normalized differential vegetation index, executes step 32；
In the present invention, normalized differential vegetation index selects 0.4 to be used as reference value, is greater than 0.4 with NDVI index or is less than
Equal to 0.4 pair rotation grazing cell is divided into two set.
From rotation grazing set of cells RE={ E1,E2,…,Ev,…EMIn obtain normalized differential vegetation index be higher than 0.4 rotation grazing
Cell is denoted as qualified rotation grazing set of cells
Indicate the rotation grazing cell E that normalized differential vegetation index (NDVI) is higher than 0.41。
Indicate the rotation grazing cell E that normalized differential vegetation index (NDVI) is higher than 0.42。
Indicate the rotation grazing cell E that normalized differential vegetation index (NDVI) is higher than 0.4v。
Indicate the rotation grazing cell E that normalized differential vegetation index (NDVI) is higher than 0.4M。
From rotation grazing set of cells RE={ E1,E2,…,Ev,…EMIn obtain normalized differential vegetation index less than or equal to 0.4
Rotation grazing cell is denoted as degeneration rotation grazing set of cells
Indicate the rotation grazing cell E that normalized differential vegetation index (NDVI) is less than or equal to 0.41。
Indicate the rotation grazing cell E that normalized differential vegetation index (NDVI) is less than or equal to 0.42。
Indicate the rotation grazing cell E that normalized differential vegetation index (NDVI) is less than or equal to 0.4v。
Indicate the rotation grazing cell E that normalized differential vegetation index (NDVI) is less than or equal to 0.4M。
In the present invention, the RENDVI_1 ∈ RE, the RENDVI_2 ∈ RE.
For example, for rotation grazing set of cells RE={ E1,E2,…,Ev,…EMIn any one rotation grazing cell returned
One, which changes vegetation index, is higher than 0.4 screening, and obtained qualified rotation grazing cell grassland has
Then degeneration rotation grazing cell grassland of the normalized differential vegetation index less than or equal to 0.4 has
For example, divided according to Alpine meadow class pasture, including buffer area L=1, total M=17 rotation grazing cell RE=
{E1,E2,E3,E4,E5,E6,E7,E8,E9,E10,E11,E12,E13,E14,E15,E16,ELIn, normalized differential vegetation index is higher than 0.4
Rotation grazing set of cells are as follows:
Then, if rotation grazing cell normalized differential vegetation index is less than or equal to 0.4 rotation grazing cell are as follows:
Step 32, appointment divided stages are carried out in a rotation grazing cycle, execute step 33；
In the present invention, distribution of the stage to rotation grazing cell is assigned by each, to obtain in a rotation grazing cycle
Interior regional rotation grazing optimization processing.
Such as, Alpine meadow class pasture, in g=1 stage, i.e. the 1st~10 of rotation grazing cycle day.
Such as, Alpine meadow class pasture, in g=2 stage, i.e. the 11st~20 of rotation grazing cycle day.
Such as, Alpine meadow class pasture, in g=3 stage, i.e. the 21st~30 of rotation grazing cycle day.
Such as, Alpine meadow class pasture, in g=4 stage, i.e. the 31st~40 of rotation grazing cycle day.
Such as, Alpine meadow class pasture, in g=5 stage, i.e. the 41st~50 of rotation grazing cycle day.
Such as, Alpine meadow class pasture, in g=6 stage, i.e. the 51st~60 of rotation grazing cycle day.
Such as, Alpine meadow class pasture, in g=7 stage, i.e. the 61st~70 of rotation grazing cycle day.
Such as, Alpine meadow class pasture, in g=8 stage, i.e. the 71st~80 of rotation grazing cycle day.
Step 33, whether seed stock herd is assigned to rotation grazing cell；
In the present invention, judge assigning stage g, whether seed stock herd is assigned to rotation grazing cell Ev, it is small according to rotation grazing
Area EvAnimal number whether be more than or equal to the quantity of seed stock herd, indicate are as follows:
It indicates assigning stage g, normalized differential vegetation index (NDVI) is higher than under conditions of 0.4, v-th of rotation grazing
Cell EvAnimal number.
It indicates assigning stage g, normalized differential vegetation index (NDVI) is higher than under conditions of 0.4, j-th of rotation grazing
Cell EjAnimal number.V-th of rotation grazing cell EvWith j-th of rotation grazing cell EjIt is not same rotation grazing cell.
IfStage g then is being assigned, preferential treatment seed stock herd RA can be at normalized differential vegetation index (NDVI)
Rotation grazing cell higher than 0.4Feeding executes step 34.
If allStage g then is being assigned, preferential treatment seed stock herd RA cannot be in normalized differential vegetation index
(NDVI) it is higher than 0.4 rotation grazing cellFeeding executes step 35.
IfStage g then is being assigned, conventional seed stock herd RB can be at normalized differential vegetation index (NDVI)
Rotation grazing cell higher than 0.4Feeding executes step 34.
If allStage g then is being assigned, conventional seed stock herd RB cannot be in normalized differential vegetation index
(NDVI) it is higher than 0.4 rotation grazing cellFeeding executes step 35.
In the present invention, judge to assign the stage in g=1, whether seed stock herd is assigned to rotation grazing cell Ev, it is according to rotation grazing
Cell EvAnimal number whether be more than or equal to the quantity of seed stock herd, indicate are as follows:
It indicates under conditions of normalized differential vegetation index (NDVI) is higher than 0.4, v-th of rotation grazing cell EvLoad
Poultry amount, subscript 1 indicate that first (g=1) assigns the stage.
Indicate that normalized differential vegetation index (NDVI) is higher than under conditions of 0.4, j-th of rotation grazing cell EjLoad poultry
Amount, subscript 1 indicate the appointment stage.
IfThen giving preferential treatment to seed stock herd RA can be in the wheel of normalized differential vegetation index (NDVI) higher than 0.4
Herd cellFeeding executes step 34.
If allThen giving preferential treatment to seed stock herd RA cannot be higher than at normalized differential vegetation index (NDVI)
0.4 rotation grazing cellFeeding executes step 35.
IfThen conventional seed stock herd RB can be in the wheel of normalized differential vegetation index (NDVI) higher than 0.4
Herd cellFeeding executes step 34.
If allThen routine seed stock herd RB cannot be higher than at normalized differential vegetation index (NDVI)
0.4 rotation grazing cellFeeding executes step 35.
Such as, when g=1, i.e. the 1st~10 day of rotation grazing cycle；
ThenMeet RAT and RBT.
Such as, when g=2, the 11st~20 day of rotation grazing cycle:
ThenMeet RAT and RBT.
Such as, when g=3, i.e. the 21st~30 day of rotation grazing cycle:
ThenMeet RAT and RBT.
Such as, when g=4, i.e. the 31st~40 day of rotation grazing cycle:
ThenMeet RAT and RBT.
Such as, when g=5, i.e. the 41st~50 day of rotation grazing cycle:
ThenMeet RAT and RBT.
Such as, when g=6, i.e. the 51st~60 day of rotation grazing cycle:
ThenMeet RAT and RBT.
Such as, when g=7, i.e. the 61st~70 day of rotation grazing cycle:
ThenMeet RAT and RBT.
Such as, when g=8, i.e. the 71st~80 day of rotation grazing cycle:
Since assignable rotation grazing cell is inadequate, then buffer area is added and is herded jointly.
ThenMeet RAT and RBT.
Step 34, result is assigned to obtain；
In the present invention, belonging to g=1 assigns the appointment result in stage to have: animal number highest rotation grazing cellWith
The high rotation grazing cell of animal number time
Under conditions of expression meets normalized differential vegetation index (NDVI), the stage is assigned to take turns v-th in g=1
It herds cell and is assigned to preferential treatment seed stock herd RA.
Under conditions of expression meets normalized differential vegetation index (NDVI), the stage is assigned in g=1, gives preferential treatment to seed stock herd
RA is in rotation grazing cell EvIn for forage total amount.
Under conditions of expression meets normalized differential vegetation index (NDVI), the stage is assigned to take turns j-th in g=1
It herds cell and is assigned to conventional seed stock herd RB.
Under conditions of expression meets normalized differential vegetation index (NDVI), stage, conventional seed stock herd are assigned in g=1
RB is in rotation grazing cell EjIn for forage total amount.
D1To assign the stage for the total amount of livestock forage in g=1.
By animal number highest rotation grazing cellIt distributes to preferential treatment seed stock herd RA to herd, the high rotation grazing of animal number time
CellConventional seed stock herd RB is distributed to herd.Execute step 35.
It indicates first appointment stage, under conditions of meeting normalized differential vegetation index (NDVI), by rotation grazing
Cell EvIt is assigned to preferential treatment seed stock herd RA.
It indicates first appointment stage, under conditions of meeting normalized differential vegetation index (NDVI), by rotation grazing
Cell EjIt is assigned to conventional seed stock herd RB.
Such as, appointment result when g=1 are as follows:It is assigned to preferential treatment seed stock herd RA,It is assigned to conventional poultry kind
Group RB.
Such as, appointment result when g=2 are as follows:It is assigned to preferential treatment seed stock herd RA,It is assigned to conventional poultry kind
Group RB.
Such as, appointment result when g=3 are as follows:It is assigned to preferential treatment seed stock herd RA,It is assigned to conventional poultry kind
Group RB.
Such as, appointment result when g=4 are as follows:It is assigned to preferential treatment seed stock herd RA,It is assigned to conventional poultry kind
Group RB.
Such as, appointment result when g=5 are as follows:It is assigned to preferential treatment seed stock herd RA,It is assigned to conventional poultry kind
Group RB.
Such as, appointment result when g=6 are as follows:It is assigned to preferential treatment seed stock herd RA,It is assigned to conventional poultry kind
Group RB.
Such as, appointment result when g=7 are as follows:It is assigned to preferential treatment seed stock herd RA,It is assigned to conventional poultry kind
Group RB.
Such as, appointment result when g=8 are as follows:It is assigned to preferential treatment seed stock herd RA,It is assigned to conventional poultry kind
Group RB.
Step 35, judge whether rotation grazing cell is searched for food by livestock；
In the present invention, assigning the stage to assign after g=1 terminates, and whether record rotation grazing cell is searched for food by livestock；To not by domestic animal
The rotation grazing cell for raiseeing feeding repeats step 32 to step 34, carries out the appointment that g=2 assigns the stage.The stage is assigned in g=2
Appointment terminates, record rotation grazing cell whether searched for food by livestock, to not by livestock search for food rotation grazing cell, repeat step 32 to
Step 34, the appointment that g=3 assigns the stage is carried out.Until theA appointment stage, which assigns, completes, and M rotation grazing cell is complete
Portion's seed stock herd is searched for food, and dynamic regional rotation grazing planning is completed, and step terminates.
In the present invention, if the animal number of all rotation grazing cells is all unsatisfactory for seed stock herd animal number in the stage assigned,
Then merge buffer area L and supplementary feeding is herded or carried out to the highest rotation grazing cell of animal number.
It is denoted as by the rotation grazing cell that livestock searches for food:Table
Show that first rotation grazing cell is searched for food by livestock；Indicate that second rotation grazing cell is searched for food by livestock；Indicate that v-th of rotation grazing is small
It is searched for food by livestock in area；Indicate that j-th of rotation grazing cell is searched for food by livestock；Indicate that m-th rotation grazing cell is searched for food by livestock.
It is not denoted as by the rotation grazing cell that livestock searches for food:
Indicate that first rotation grazing cell is not searched for food by livestock；Indicate that second rotation grazing cell is not searched for food by livestock；It indicates v-th
Rotation grazing cell is not searched for food by livestock；Indicate that j-th of rotation grazing cell is not searched for food by livestock；Indicate m-th rotation grazing cell not
It is searched for food by livestock.In the present invention, the rotation grazing cell that do not searched for food by livestock is exactly unassigned rotation grazing cell.
If assigning the stage in g=1, the rotation grazing cell searched for food by livestock are as follows:Then not by domestic animal
Raise the rotation grazing cell of feeding are as follows:
The stage such as is assigned in g=1,
The rotation grazing cell searched for food by livestock are as follows:Indicate the 7th rotation grazing cell by livestock
Feeding；Indicate that the 12nd rotation grazing cell is searched for food by livestock.
The rotation grazing cell that do not searched for food by livestock are as follows:ELIt indicates not
The buffer area searched for food by livestock.
The stage is assigned in g=2,
The rotation grazing cell searched for food by livestock are as follows:Indicate that the 4th rotation grazing is small
It is searched for food by livestock in area；Indicate that the 7th rotation grazing cell is searched for food by livestock；Indicate that the tenth rotation grazing cell is adopted by livestock
Food；Indicate that the 12nd rotation grazing cell is searched for food by livestock.
The rotation grazing cell that do not searched for food by livestock are as follows:
The stage is assigned in g=3,
The rotation grazing cell searched for food by livestock are as follows:
The rotation grazing cell that do not searched for food by livestock are as follows:
The stage is assigned in g=4,
The rotation grazing cell searched for food by livestock are as follows:
The rotation grazing cell that do not searched for food by livestock are as follows: REO={ E3,E6,E8,E13,E14,E15,E16,EL}。
The stage is assigned in g=5,
The rotation grazing cell searched for food by livestock are as follows:
The rotation grazing cell that do not searched for food by livestock are as follows: REO={ E6,E8,E13,E14,E15,EL}。
The stage is assigned in g=6,
The rotation grazing cell searched for food by livestock are as follows:
The rotation grazing cell that do not searched for food by livestock are as follows: REO={ E8,E13,E15,EL}。
The stage is assigned in g=7,
The rotation grazing cell searched for food by livestock are as follows:
The rotation grazing cell that do not searched for food by livestock are as follows: REO={ E15,EL}。
The stage is assigned in g=8,
The rotation grazing cell searched for food by livestock are as follows:
Indicate that buffer area is searched for food by livestock.
The rotation grazing cell that do not searched for food by livestock are as follows:Then dynamic divides the small end of extent of rotation grazing.Referring to Fig. 3 institute
Show, with the method for the present invention and be randomly assigned method, the livestock number comparison in the rotation grazing cell in different appointment stages can obtain this
Inventive method can make the livestock number of the small area carrier of rotation grazing more reasonable, can play the effect of protection rotation grazing cell ecology.
Dynamic optimization based on Task Assignment Model are as follows:
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of appointment stage by v-th
Rotation grazing cell is assigned to preferential treatment seed stock herd RA.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of appointment stage gives preferential treatment to seed stock herd
RA is in rotation grazing cell EvIn for forage total amount.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of appointment stage by j-th
Rotation grazing cell is assigned to conventional seed stock herd RB.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of appointment stage, conventional seed stock herd
RB is in rotation grazing cell EjIn for forage total amount.
DgFor the total amount g appointment stage for livestock forage.
Considering livestock feeding in the present invention influences the forward direction of livestock productivity, therefore DgIt is maximized such as formula (3):
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of stage, v-th
Rotation grazing cell can only be herded for preferential treatment seed stock herd RA.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of stage, preferential treatment
Seed stock herd RA can only enter a rotation grazing cell and herd.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of stage, j-th
Rotation grazing cell can only be herded for conventional seed stock herd RB.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of stage, routine
Seed stock herd RB can only enter a rotation grazing cell and herd.
IfThen giving preferential treatment to seed stock herd RA can be in the wheel of normalized differential vegetation index (NDVI) higher than 0.4
Herd cellFeeding.
If0.4 cannot be higher than at normalized differential vegetation index (NDVI) by then giving preferential treatment to seed stock herd RA
Rotation grazing cellFeeding.
IfThen conventional seed stock herd RB can be in the wheel of normalized differential vegetation index (NDVI) higher than 0.4
Herd cellFeeding.
IfThen conventional seed stock herd RB cannot be higher than 0.4 at normalized differential vegetation index (NDVI)
Rotation grazing cellFeeding.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of stage, preferential treatment poultry kind
Group RA is herded to v rotation grazing cell.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of stage, preferential treatment poultry
Population RA is herded to v rotation grazing cell of failing in exam.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of stage, conventional poultry kind
Group RB is herded to jth rotation grazing cell.
Under conditions of expression meets normalized differential vegetation index (NDVI), g-th of stage, conventional poultry kind
Group RB is herded to not jth rotation grazing cell.
Claims (5)
1. a kind of regional rotation grazing dynamic optimization method based on Task Assignment Model, it is characterised in that specifically there is the following steps:
Step 1: configuration rotation grazing cell；
Step 11, it according to number of days is herded rotation grazing cycle and partition to any one pasture progress small Division of rotation grazing, is had
Divide the M block rotation grazing cell in region；Execute step 12；
RE={ E is expressed as using aggregate form to the rotation grazing cell that division obtains1,E2,…,Ev,…EM}；
E1Indicate first rotation grazing cell for belonging to pasture；
E2Indicate second rotation grazing cell for belonging to pasture；
EvIndicate any one the rotation grazing cell for belonging to pasture, subscript v indicates the identification number of rotation grazing cell；
EMIndicate the last one the rotation grazing cell for belonging to pasture, subscript M indicates the sum of rotation grazing cell；
The sum M of rotation grazing cell are as follows:
X indicates rotation grazing cycle, and Y indicates that cell herds number of days, and L indicates buffer area, is alternative rotation grazing cell；
Step 12, preferential treatment seed stock herd RA and conventional seed stock herd RB is respectively divided in the livestock that herdsman raises, executes step 13；
Preferential treatment seed stock herd RA refers to the livestock in growth stage or lactation period；
Conventional seed stock herd RB refers to the livestock in maintenance phase or non-lactation period；
Step 13, the appointment stage g in rotation grazing cycle is determined by rotation grazing cell total number M；Execute step 31；
The preferential treatment seed stock herd RA and routine seed stock herd RB divided according to step 12 is assigned in M rotation grazing cell respectively, therefore,
The appointment stage sum of rotation grazing cell included in rotation grazing cycle is
Step 2: initialization livestock number；
Step 21, after step 11 processing, fence design is carried out to pasture, the phase of herding is herded beginning, and executes step 22；
Step 22, the livestock number of standard sheep unit is obtained；Execute step 31；
Any one kind livestock of preferential treatment seed stock herd RA is denoted asIt corresponds to describedSheep unit conversion factor be denoted asIt corresponds to describedFeeding quantity be denoted as
Any one kind livestock of conventional seed stock herd RB is denoted asIt corresponds to describedSheep unit conversion factor be denoted asIt corresponds to describedFeeding quantity be denoted as
For the livestock head number of unified different seed stock herds, system is converted multiplied by the sheep unit of different cultivars using the quantity of different cultivars
Number, obtains the number of the standard sheep unit of each seed stock herd；
For the livestock number of the standard sheep unit of preferential treatment seed stock herd RA, it is denoted as RAT, and
For the livestock number of the standard sheep unit of conventional seed stock herd RB, it is denoted as RBT, and
P indicates the mark of the kind of the livestock of herdsman's raising；
Q indicates the kind sum of the livestock of herdsman's raising；
Indicate the feeding quantity for belonging to the livestock kind of preferential treatment seed stock herd RA；
Indicate the conversion factor for belonging to the livestock kind of preferential treatment seed stock herd RA；
Indicate the feeding quantity for belonging to the livestock kind of conventional seed stock herd RB；
Indicate the conversion factor for belonging to the livestock kind of conventional seed stock herd RB；
Step 3: the dynamic regional rotation grazing planning based on Task Assignment Model；
Step 31, the small Division of rotation grazing is carried out using normalized differential vegetation index, executes step 32；
From rotation grazing set of cells RE={ E1,E2,…,Ev,…EMIn obtain normalized differential vegetation index be higher than 0.4 rotation grazing cell,
It is denoted as qualified rotation grazing set of cells
Indicate the rotation grazing cell E that normalized differential vegetation index is higher than 0.41；
Indicate the rotation grazing cell E that normalized differential vegetation index is higher than 0.42；
Indicate the rotation grazing cell E that normalized differential vegetation index is higher than 0.4v；
Indicate the rotation grazing cell E that normalized differential vegetation index is higher than 0.4M；
From rotation grazing set of cells RE={ E1,E2,…,Ev,…EMIn obtain normalized differential vegetation index less than or equal to 0.4 rotation grazing it is small
Area is denoted as degeneration rotation grazing set of cells
Indicate the rotation grazing cell E that normalized differential vegetation index is less than or equal to 0.41；
Indicate the rotation grazing cell E that normalized differential vegetation index is less than or equal to 0.42；
Indicate the rotation grazing cell E that normalized differential vegetation index is less than or equal to 0.4v；
Indicate the rotation grazing cell E that normalized differential vegetation index is less than or equal to 0.4M；
In the present invention, the RENDVI_1 ∈ RE, the RENDVI_2 ∈ RE；
Step 32, appointment divided stages are carried out in a rotation grazing cycle, execute step 33；
Step 33, whether seed stock herd is assigned to rotation grazing cell；
Judge assigning stage g, whether seed stock herd is assigned to rotation grazing cell Ev, it is according to rotation grazing cell EvAnimal number whether
More than or equal to the quantity of seed stock herd, indicate are as follows:
It indicates assigning stage g, normalized differential vegetation index is higher than under conditions of 0.4, v-th of rotation grazing cell EvLoad poultry
Amount；
It indicates assigning stage g, normalized differential vegetation index is higher than under conditions of 0.4, j-th of rotation grazing cell EjLoad poultry
Amount；V-th of rotation grazing cell EvWith j-th of rotation grazing cell EjIt is not same rotation grazing cell；
IfStage g then is being assigned, preferential treatment seed stock herd RA can be in wheel of the normalized differential vegetation index higher than 0.4
Herd cellFeeding executes step 34；
If allStage g then is being assigned, giving preferential treatment to seed stock herd RA cannot be higher than in normalized differential vegetation index
0.4 rotation grazing cellFeeding executes step 35；
IfThen it can be higher than 0.4 in normalized differential vegetation index in appointment stage g, conventional seed stock herd RB
Rotation grazing cellFeeding executes step 34；
If allIt cannot then be higher than in normalized differential vegetation index in appointment stage g, conventional seed stock herd RB
0.4 rotation grazing cellFeeding executes step 35；
Step 34, result is assigned to obtain；
Any one assigns the appointment result in stage to have: animal number highest rotation grazing cellThe high rotation grazing with animal number time
Cell
Step 35, judge whether rotation grazing cell is searched for food by livestock；
Assigning the stage to assign after g=1 terminates, and whether record rotation grazing cell is searched for food by livestock；It is small to the rotation grazing that do not searched for food by livestock
Area repeats step 32 to step 34, carries out the appointment that g=2 assigns the stage；Assigning the stage to assign in g=2 terminates, records
Whether rotation grazing cell is searched for food by livestock, to the rotation grazing cell that do not searched for food by livestock, repeats step 32 to step 34, carries out g
=3 assign the appointment in stage；Until theA appointment stage, which assigns, completes, and M rotation grazing cell adopted by whole seed stock herds
Food, dynamic regional rotation grazing planning are completed, and step terminates；
If the animal number of all rotation grazing cells is all unsatisfactory for seed stock herd animal number in the stage assigned, merge buffer area L with
Supplementary feeding is herded or carried out to the highest rotation grazing cell of animal number；
It is denoted as by the rotation grazing cell that livestock searches for food: Indicate first
A rotation grazing cell is searched for food by livestock；Indicate that second rotation grazing cell is searched for food by livestock；Indicate v-th of rotation grazing cell by domestic animal
Poultry feeding；Indicate that j-th of rotation grazing cell is searched for food by livestock；Indicate that m-th rotation grazing cell is searched for food by livestock；
It is not denoted as by the rotation grazing cell that livestock searches for food: It indicates
First rotation grazing cell is not searched for food by livestock；Indicate that second rotation grazing cell is not searched for food by livestock；Indicate v-th of rotation grazing
Cell is not searched for food by livestock；Indicate that j-th of rotation grazing cell is not searched for food by livestock；Indicate m-th rotation grazing cell not by domestic animal
Poultry feeding.
2. the regional rotation grazing dynamic optimization method according to claim 1 based on Task Assignment Model, it is characterised in that based on finger
Send the dynamic optimization of model are as follows:
Under conditions of expression meets normalized differential vegetation index, v-th of rotation grazing cell is referred to g-th of appointment stage
Task preferential treatment seed stock herd RA；
Under conditions of expression meets normalized differential vegetation index, g-th of appointment stage gives preferential treatment to seed stock herd RA in rotation grazing cell
EvIn for forage total amount；
Under conditions of expression meets normalized differential vegetation index, j-th of rotation grazing cell is referred to g-th of appointment stage
Task conventional seed stock herd RB；
Under conditions of expression meets normalized differential vegetation index, in g-th of appointment stage, conventional seed stock herd RB is in rotation grazing cell
EjIn for forage total amount；
DgFor the total amount g appointment stage for livestock forage；
The forward direction of livestock productivity is influenced considering livestock feeding, therefore DgIt is maximized are as follows:
Under conditions of expression meets normalized differential vegetation index, g-th of stage, v-th of rotation grazing cell is only
It can be herded for preferential treatment seed stock herd RA；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, seed stock herd RA is given preferential treatment to only
It can enter a rotation grazing cell to herd；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, j-th of rotation grazing cell is only
It can be herded for conventional seed stock herd RB；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, conventional seed stock herd RB is only
It can enter a rotation grazing cell to herd；
IfThen giving preferential treatment to seed stock herd RA can be in rotation grazing cell of the normalized differential vegetation index higher than 0.4Feeding；
IfThen giving preferential treatment to seed stock herd RA cannot be in rotation grazing cell of the normalized differential vegetation index higher than 0.4Feeding；
IfThen conventional seed stock herd RB can be in rotation grazing cell of the normalized differential vegetation index higher than 0.4Feeding；
IfThen conventional seed stock herd RB cannot be in rotation grazing cell of the normalized differential vegetation index higher than 0.4Feeding；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, seed stock herd RA to v is given preferential treatment to
Rotation grazing cell is herded；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, seed stock herd RA is given preferential treatment to failing in exam
V rotation grazing cell is herded；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, conventional seed stock herd RB to jth wheel
Cell is herded to herd；
Under conditions of expression meets normalized differential vegetation index, g-th of stage, conventional seed stock herd RB is to failing in exam
J rotation grazing cell is herded.
3. the regional rotation grazing dynamic optimization method according to claim 1 or 2 based on Task Assignment Model, it is characterised in that: return
One, which changes vegetation index, selects 0.4 to be used as reference value, is greater than 0.4 with NDVI index or is less than or equal to 0.4 pair of small Division of rotation grazing
Gather for two.
4. the regional rotation grazing dynamic optimization method according to claim 1 or 2 based on Task Assignment Model, it is characterised in that: wheel
Number of cells M is herded to round up.
5. the regional rotation grazing dynamic optimization method according to claim 1 or 2 based on Task Assignment Model, it is characterised in that: excellent
Feed intake difference to seed stock herd RA and conventional seed stock herd RB is no more than ± 5%.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910722731.XA CN110443423B (en) | 2019-08-06 | 2019-08-06 | Dynamic optimization method for grazing in turn in zoning based on assignment model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910722731.XA CN110443423B (en) | 2019-08-06 | 2019-08-06 | Dynamic optimization method for grazing in turn in zoning based on assignment model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110443423A true CN110443423A (en) | 2019-11-12 |
| CN110443423B CN110443423B (en) | 2023-01-17 |
Family
ID=68433699
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910722731.XA Active CN110443423B (en) | 2019-08-06 | 2019-08-06 | Dynamic optimization method for grazing in turn in zoning based on assignment model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110443423B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023036057A1 (en) * | 2021-09-08 | 2023-03-16 | International Business Machines Corporation | Disaster prediction and proactive mitigation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1650699A (en) * | 2004-02-06 | 2005-08-10 | 福建省农业科学院红萍研究中心 | Method for rotational grazing of southern mountain artificial grassland goats in different regions |
| US20080109196A1 (en) * | 2006-11-06 | 2008-05-08 | The Remme Corporation | Tools and Methods for Range Management |
| CN103444625A (en) * | 2013-08-28 | 2013-12-18 | 百绿国际草业(北京)有限公司 | Rotational grazing method |
| CN104765961A (en) * | 2015-04-08 | 2015-07-08 | 兰州大学 | Rotational grazing method for Tan sheep in grassland |
-
2019
- 2019-08-06 CN CN201910722731.XA patent/CN110443423B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1650699A (en) * | 2004-02-06 | 2005-08-10 | 福建省农业科学院红萍研究中心 | Method for rotational grazing of southern mountain artificial grassland goats in different regions |
| US20080109196A1 (en) * | 2006-11-06 | 2008-05-08 | The Remme Corporation | Tools and Methods for Range Management |
| CN103444625A (en) * | 2013-08-28 | 2013-12-18 | 百绿国际草业(北京)有限公司 | Rotational grazing method |
| CN104765961A (en) * | 2015-04-08 | 2015-07-08 | 兰州大学 | Rotational grazing method for Tan sheep in grassland |
Non-Patent Citations (2)
| Title |
|---|
| M.J. CROS等: "Simulating management strategies: the rotational grazing example", 《AGRICULTURAL SYSTEMS》 * |
| 覃琳等: "基于划区轮牧理论的祁连山北麓牧区动态草畜平衡研究", 《农业工程学报》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023036057A1 (en) * | 2021-09-08 | 2023-03-16 | International Business Machines Corporation | Disaster prediction and proactive mitigation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110443423B (en) | 2023-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ryschawy et al. | Paths to last in mixed crop–livestock farming: lessons from an assessment of farm trajectories of change | |
| Putfarken et al. | Site use of grazing cattle and sheep in a large-scale pasture landscape: A GPS/GIS assessment | |
| Mbuza et al. | Characterization of broiler poultry production system in Rwanda | |
| Gelasakis et al. | Description and typology of intensive Chios dairy sheep farms in Greece | |
| Steinfeld et al. | A classification of livestock production systems | |
| Cossins et al. | Options for improvement of the Borana pastoral system | |
| Rada et al. | Productivity growth and the revival of Russian agriculture | |
| Johnson | Agricultural adjustment in China: Problems and prospects | |
| Chen et al. | Chinese integrated fish farming: a comparative bioeconomic analysis | |
| CN110443423A (en) | A kind of regional rotation grazing dynamic optimization method based on Task Assignment Model | |
| Omasaki et al. | Fish production practices of smallholder farmers in western Kenya | |
| CN109329201A (en) | Family domestic animal-meadow configuration technology is herded by a kind of family | |
| Kahi et al. | Influence of production circumstances and economic evaluation criteria on economic comparison of breeds and breed crosses | |
| Ørskov et al. | The role of animals in spreading farmers’ risks: a new paradigm for animal science | |
| Kaur et al. | Cost of milk production in Punjab: A pre-requisite for pricing policy | |
| Cosgrove et al. | High production dairy-beef cattle grazing systems: a review of research in the Manawatu | |
| Tegegne | Transhumance cattle production system in North Gondar, Amhara region, Ethiopia: is it sustainable? | |
| Stanbury | Women and water: Effects of irrigation development in a north Indian village | |
| CN110458368A (en) | A kind of nomadic paths planning method of Forage-Livestock Balance based on network flow | |
| Mbuza et al. | Characterization of layer poultry production in Rwanda | |
| Herrera et al. | Prospective analysis of the feeding basis in a tropical dairy farm with linear programming | |
| Soliman | Economic feed utilization for dairy buffalo under intensive agricultural system | |
| Niu et al. | Towards an Agricultural Economy for China in a New Age: Progress, Problems, Response, and Prospects | |
| Lestari et al. | " Gaduh Culture" as a Model of Capital Access for the Poor Through Cows Cattle in Trenggalek | |
| Haque et al. | An analysis of Profitability and Resource Use Efficiency of Poultry Feed Mills in Bangladesh |
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 |