一种奶牛存栏结构预测方法与流程

本发明涉及畜牧养殖领域，特别涉及一种奶牛存栏结构预测方法。

背景技术：

在畜牧领域中，传统的奶牛存栏结构预测，目前还没有一种明确的算法，仅仅是根据年受胎率等指标进行粗略的简单人为估算，没有一个完整的把影响牛群结构变化的可计算因素都进行参数化并参与相关计算得到指定时间段内牛群结构计算方法及公式，因此偏差比较大。

技术实现要素：

为了有效估算出奶牛群存栏结构，以便牧场做下一年的生产计划，并能跟踪计划与实际对比，及时发现生产异常情况；本发明提供一种奶牛存栏结构预测方法。

其技术方案为，以现有的牛群结构及生产数据为基础，根据奶牛的繁殖逻辑及牛群状态的变化规律，总结出一套系统性的计算公式，并根据积累的牛群档案数据及繁殖登记数据计算得到各群的预测头数；

一种奶牛存栏结构预测方法：

S1建立初始数据库，包括：

S101、根据牧场所有牛群的孕育状态及天数将牛群进行分类；

S102、对牧场牛群进行至少2个月时间的状态统计；包括牛只档案及繁殖登记数据；并指定数据开始日期，根据指定日期统计出相关牛群结构，作为初始数据，形成初始数据库；初始数据包括牛群的期初值、淘汰数及新增数；

S2、调用S1的初始数据库，获取牛群参数：并将获取的参数分配在S1分类的各牛群上，然后再针对各个牛群，将参数分配至1到12月；

S3、采用计算函数针对第一个月进行计算，通过计算，获取未来12个月的数据；计算函数为期初值减去指定月份预计减少数加上新增数；

S4、对S3得到的结果，通过实际数据校验进行纠正；即针对使用的各个参数，结合实际情况用同期最近实际数据对比判断结果是否与实际存在较大偏差，如果偏差大，则以参数实际计算值直接替代；这样预测的结果会更加贴近实际结果。

优选为，所述S1中将牛群进行分类具体为，将牛群分成：未孕泌乳牛、已孕泌乳牛、已孕干奶牛、未孕干奶牛、已孕育成牛、未孕育成牛、各个月龄育成牛、哺乳犊牛、断奶犊牛，再将牛群按泌乳天数、怀孕天数、月龄细分；未孕泌乳牛按泌乳天数细分，分为0~30、30~60、60以上；妊娠泌乳牛及干奶牛及后备牛按妊娠天数细分，按妊娠天数45~75、75~105、105~135、135~165、165~195、195~225、225~255、255以上；未孕干奶牛不再细分；未孕后备牛按出生天数细分，每30天为一个阶段，直到450以及大于450天。

优选为，所述S2中获取的牛群参数包括：泌乳牛年度化死淘率、干奶牛年度化死淘率、青年牛年度化死淘率、犊牛年度化死淘率、泌乳牛怀孕0-120 天流产率、泌乳牛怀孕120-210 天流产率、干奶牛怀孕0-120 天流产率、干奶牛怀孕120-210 天流产率、育成牛怀孕0-120 天流产率、育成牛怀孕120-210 天流产率、成乳牛怀孕210-265天早产率、青年牛怀孕210-265天早产率、出生成活率、母犊率、未孕泌乳牛参配率、未孕干奶牛参配率、青年牛参配率、未孕泌乳牛受胎率、未孕干奶牛受胎率、青年牛受胎率、计划单产。

优选为，所述S3中计算步骤包括：

S301，对初始数据库进行调取，获取各项所述牛群参数，以及预先计算出需要扣减的头数；比如流产率，可以根据该比率可以预算出妊娠牛只不能正常产犊的头数，用期初数值乘以这个比率就可以计算出需要扣减的妊娠牛头数从而得到预测的妊娠牛头数，其他如母犊率、淘汰率等也是这样的计算方法；

S302，针对指定计算时间里的第一个月，以上月期末值为基础，结合系统提供的参数值进行计算，下个月将以上个月计算得到的头数作为期初数结合相关公式及参数进行计算，以此类推，一直到第12个月；

S303，通过积累的实际数据动态更新各个参数，将偏差缩小；各个参数根据实际数据根据前面叙述的计算方法进行定期重新计算，用实际计算得到的值更新参数表，然后用这些参数再次参与计算得到预测的牛只头数。同时本系统较好的解决了类似这种偏差，并根据丰顿奶牛云平台积累的数据进行大数据分析总结出不同区域的相关参数正常值，为预测的准确度打下良好基础。

优选为，所述S302的计算过程为：

预测指定牛群本期死淘头数=∑指定牛群期初头数*年度化死淘率；

预测指定牛群本期转入=∑指定牛群期初头数*指定的年度化指标；

预测指定牛群本期转出=∑指定牛群期初头数*指定的年度化指标；

预测指定牛群本期头数=指定牛群期初头数+预测指定牛群本期转入-预测指定牛群本期转出-预测指定牛群本期死淘头数；

所述年度化死淘率：根据已有数据提取实际淘汰头数除以指定时间范围内每天饲养头数的平均值；包括泌乳牛年度化死淘率、干奶牛年度化死淘率、青年牛年度化死淘率、犊牛年度化死淘率；

所述指定的年度化指标：包括泌乳牛怀孕0-120 天流产率、泌乳牛怀孕120-210 天流产率、干奶牛怀孕0-120 天流产率、干奶牛怀孕120-210 天流产率、育成牛怀孕0-120 天流产率、育成牛怀孕120-210 天流产率、成乳牛怀孕210-265天早产率、青年牛怀孕210-265天早产率、出生成活率、母犊率、未孕泌乳牛参配率、未孕干奶牛参配率、青年牛参配率、未孕泌乳牛受胎率、未孕干奶牛受胎率、青年牛受胎率。

流产率的计算就是指定时间范围内指定妊娠天数的牛只流产发生数除以指定妊娠天数阶段牛只头数；

早产率就是指定时间范围内指定怀孕阶段的牛只流产头数除以指定怀孕阶段的牛只总头数；

出生成活率就是指定时间段内出生牛只成活头数除以出生牛只总数。母犊率就是指定时间范围内出生牛只母犊数除以出生牛只总数；

参配率就是用指定时间段内实际配种头数除以理论发情牛数；

受胎率就是指定时间段内初检怀孕的牛只头数除以配种头数。

根据积累数据，分别计算出牛只分群得到的各个牛群的期初头数。

所有参数需要计算出在各个月份的具体数值。

在期初数据基础上计算出各群牛只期初头数，并设置好个月各群系列参数，先计算指定日期后的第一个月，再以第一个月的得数为初始值计算下个月以此类推。

计算过程可能存在累计偏差，通过调整参数来消除偏差。

积累基础数据，分牛群计算期初牛只头数，再将涉及到的相关计算参数按月份计算或直接设置，再进行计算指定日期后的第一个月的各群牛头数，然后以第一个月的计算结果做第二个月的期初牛头数，再以相同方式进行计算，以此类推。

本发明实施例提供的技术方案带来的有益效果是：该计算方法可以有效预算出未来一年的各个阶段牛只头数，这样为牛场本身可以提供相关决策依据，避免浪费；对乳企可以提前预算处来年各个月可以从各个牧场采购到的牛奶数量，安排好相关生产及销售计划，做到以销定产，为企业避免相关损失，节约成本。

附图说明

图1为本发明实施例的原理框图；

图2为本发明实施例的泌乳牛数量推测线图；

图3为本发明实施例产奶量推测线图；

图4为本发明实施例的奶牛数据展示图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。当然，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

实施例1

参见图1至图4，本发明提供一种奶牛存栏结构预测方法：

S1建立初始数据库，包括：

S101、根据牧场所有牛群的孕育状态及天数将牛群进行分类；

S102、对牧场牛群进行至少2个月时间的状态统计；包括牛只档案及繁殖登记数据；并指定数据开始日期，根据指定日期统计出相关牛群结构，作为初始数据，形成初始数据库；初始数据包括牛群的期初值、淘汰数及新增数；

S2、调用S1的初始数据库，获取牛群参数：并将获取的参数分配在S1分类的各牛群上，然后再针对各个牛群，将参数分配至1到12月；

S3、采用计算函数针对第一个月进行计算，通过计算，获取未来12个月的数据；计算函数为期初值减去指定月份预计减少数加上新增数；

S4、对S3得到的结果，通过实际数据校验进行纠正；即针对使用的各个参数，结合实际情况用同期最近实际数据对比判断结果是否与实际存在较大偏差，如果偏差大，则以参数实际计算值直接替代；这样预测的结果会更加贴近实际结果。

优选为，所述S1中将牛群进行分类具体为，将牛群分成：未孕泌乳牛、已孕泌乳牛、已孕干奶牛、未孕干奶牛、已孕育成牛、未孕育成牛、各个月龄育成牛、哺乳犊牛、断奶犊牛，再将牛群按泌乳天数、怀孕天数、月龄细分；未孕泌乳牛按泌乳天数细分，分为0~30、30~60、60以上；妊娠泌乳牛及干奶牛及后备牛按妊娠天数细分，按妊娠天数45~75、75~105、105~135、135~165、165~195、195~225、225~255、255以上；未孕干奶牛不再细分；未孕后备牛按出生天数细分，每30天为一个阶段，直到450以及大于450天。

优选为，所述S2中获取的牛群参数包括：泌乳牛年度化死淘率、干奶牛年度化死淘率、青年牛年度化死淘率、犊牛年度化死淘率、泌乳牛怀孕0-120 天流产率、泌乳牛怀孕120-210 天流产率、干奶牛怀孕0-120 天流产率、干奶牛怀孕120-210 天流产率、育成牛怀孕0-120 天流产率、育成牛怀孕120-210 天流产率、成乳牛怀孕210-265天早产率、青年牛怀孕210-265天早产率、出生成活率、母犊率、未孕泌乳牛参配率、未孕干奶牛参配率、青年牛参配率、未孕泌乳牛受胎率、未孕干奶牛受胎率、青年牛受胎率、计划单产。

优选为，所述S3中计算步骤包括：

S301，对初始数据库进行调取，获取各项所述牛群参数，以及预先计算出需要扣减的头数；比如流产率，可以根据该比率可以预算出妊娠牛只不能正常产犊的头数，用期初数值乘以这个比率就可以计算出需要扣减的妊娠牛头数从而得到预测的妊娠牛头数，其他如母犊率、淘汰率等也是这样的计算方法；

S302，针对指定计算时间里的第一个月，以上月期末值为基础，结合系统提供的参数值进行计算，下个月将以上个月计算得到的头数作为期初数结合相关公式及参数进行计算，以此类推，一直到第12个月；

S303，通过积累的实际数据动态更新各个参数，将偏差缩小；各个参数根据实际数据根据前面叙述的计算方法进行定期重新计算，用实际计算得到的值更新参数表，然后用这些参数再次参与计算得到预测的牛只头数。同时本系统较好的解决了类似这种偏差，并根据丰顿奶牛云平台积累的数据进行大数据分析总结出不同区域的相关参数正常值，为预测的准确度打下良好基础。

优选为，所述S302的计算过程为：

预测指定牛群本期死淘头数=∑指定牛群期初头数*年度化死淘率；

预测指定牛群本期转入=∑指定牛群期初头数*指定的年度化指标；

预测指定牛群本期转出=∑指定牛群期初头数*指定的年度化指标；

预测指定牛群本期头数=指定牛群期初头数+预测指定牛群本期转入-预测指定牛群本期转出-预测指定牛群本期死淘头数；

所述年度化死淘率：根据已有数据提取实际淘汰头数除以指定时间范围内每天饲养头数的平均值；包括泌乳牛年度化死淘率、干奶牛年度化死淘率、青年牛年度化死淘率、犊牛年度化死淘率；

所述指定的年度化指标：包括泌乳牛怀孕0-120 天流产率、泌乳牛怀孕120-210 天流产率、干奶牛怀孕0-120 天流产率、干奶牛怀孕120-210 天流产率、育成牛怀孕0-120 天流产率、育成牛怀孕120-210 天流产率、成乳牛怀孕210-265天早产率、青年牛怀孕210-265天早产率、出生成活率、母犊率、未孕泌乳牛参配率、未孕干奶牛参配率、青年牛参配率、未孕泌乳牛受胎率、未孕干奶牛受胎率、青年牛受胎率。

流产率的计算就是指定时间范围内指定妊娠天数的牛只流产发生数除以指定妊娠天数阶段牛只头数；

早产率就是指定时间范围内指定怀孕阶段的牛只流产头数除以指定怀孕阶段的牛只总头数；

出生成活率就是指定时间段内出生牛只成活头数除以出生牛只总数。母犊率就是指定时间范围内出生牛只母犊数除以出生牛只总数；

参配率就是用指定时间段内实际配种头数除以理论发情牛数；

受胎率就是指定时间段内初检怀孕的牛只头数除以配种头数。

根据积累数据，分别计算出牛只分群得到的各个牛群的期初头数。

所有参数需要计算出在各个月份的具体数值。

在期初数据基础上计算出各群牛只期初头数，并设置好个月各群系列参数，先计算指定日期后的第一个月，再以第一个月的得数为初始值计算下个月以此类推。

计算过程可能存在累计偏差，通过调整参数来消除偏差。

积累基础数据，分牛群计算期初牛只头数，再将涉及到的相关计算参数按月份计算或直接设置，再进行计算指定日期后的第一个月的各群牛头数，然后以第一个月的计算结果做第二个月的期初牛头数，再以相同方式进行计算，以此类推。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。