奶牛泌乳曲线的应用

泌乳牛在整个泌乳期的变化规律2012-01-16 13:51:38 来源：牛业科技之光浏览：14次泌乳牛在整个泌乳期内，其产奶量、体重、干物质进食量的变化规律如下：（一）产奶量的变化母牛产犊后，在正常情况下产奶量逐日增加，到产后6～10周出现最高日产奶量，其产奶量比初乳期可提高30％。

若按泌乳月计算，则第二个泌乳月为产奶量最高峰。

此后3～5个月为平稳期，高产奶牛其曲线基本不变，有些母牛每月可下降3％～5％，第六至十个月为泌乳后期，高产奶牛下降得很少或不下降，有些牛每月平均下降3％～5％，另一些牛则下降7％～8％，直到停奶后进入干奶期。

（二）体重的变化据测定，饲养得好的、体重较稳定的母牛，产后0～70 天体重共减少35千克，每天平均约减500克；其减重时间分布，前30 天平均每日每头减重1.8～2.2千克，后40天仅为30～100克。

产后71～150天阶段，一般母牛和高产奶牛可保持体重不增不减，但中低产母牛此时日增重可达150克。

产后151～305 天阶段，一般奶牛均能明显恢复原来体重，包括胎儿的生长发育，综合起来日增重可达400～500克。

以后60天的干乳期，其日增重为350～600克。

(三)干物质进食量的变化母牛产犊后,随着产奶量的逐日增加，日粮干物质进食量也在增加,但在产后0～70天或0～90天阶段,母牛进食的养分不能满足泌乳养分的需要,故发生能量和蛋白质的负平衡,母牛减重用于产奶,一般于产后90 天左右母牛干物质进食量达最高峰,以后基本保持平衡,养分能满足维持和泌乳的需要,且有部分养分用于修补体组织和胎儿的需要.一般在产后7～8 个月母牛对于干物质的进食量开始下降,但高产奶牛一直保持到第九至十个泌乳月也不会下降. 干奶期母牛对于干物质的进食量下降十分明显,高产奶牛可下降40%～50% ，到临产前5～7 天食欲最差，干物质进食量仅为体重的1%。

奶牛生产性能测定（DHI）（一）体细胞计数的应用牛奶体细胞通常由巨噬细胞、淋巴细胞和多形核嗜中性白细胞（PMN）等组成。

正常情况下，牛奶中的体细胞数一般在20 万个/mL－30 万个/mL。

当乳房受到外伤或者发生疾病（如乳房炎等）时体细胞数就会迅速增加。

如果体细胞数超过50 万个/mL，就导致产奶量下降。

测量牛奶体细胞数的变化有助于及早发现乳房损伤或感染、预防治疗乳腺炎，同时还可降低治疗费用，减少牛只的淘汰，增加产奶能力。

因此，体细胞数反映了牛奶产量、质量以及牛只的健康状况，也是奶牛乳房健康水平的重要标志。

1、奶牛理想的体细胞数：第1 胎≤15 万/ml 第2 胎≤25 万/ml 第3 胎≤30 万/ml 影响体细胞数变化的主要因素有：病原微生物对乳腺组织感染、应激、环境、气候、泌乳天数、遗传、胎次等，其中致病菌影响最大。

临床乳房炎的发生，将会损失20％～70％的奶量，个别牛只甚至会无乳汁分泌。

例如：一个牛场有泌乳牛300 头，体细胞数平均40 万/ml，一年仅奶产量损失的费用就可达4.5 万元（假设头胎牛占25%，奶价2.2 元/千克)，这其中还不包括因乳腺炎造成的其它损失，如乳房永久性破坏、牛只间相互传染、头胎牛过早干奶与淘汰、兽药费、抗生素残留奶、原料奶质量下降等，约占总费用的36%。

2、体细胞数对奶牛乳房健康及牛奶品质的影响：测定牛奶体细胞，是判断乳房炎轻重的有力手段，特别是能预示隐性乳房炎。

奶牛一旦患有乳房炎，产奶量、奶的质量都会有相应的变化。

患乳房炎的奶牛其乳腺组织的泌乳能力下降，达不到遗传潜力的产奶峰值，并对干奶牛的治疗花费较大。

如果能有效的避免乳房炎，就可达到高的产奶峰值，获得巨大的经济回报。

患乳房炎的奶牛所分泌的牛奶与正常牛奶的主要区别，是干物质含量减少及各种乳成分的含量比例发生变化。

如乳房炎达到很重的程度，牛奶将接近血液成分。

所以，牛奶体细胞数与产奶量是成反比关系，高体细胞数牛奶中脂肪、蛋白、乳糖等成分都将发生变化。

实验四乳牛泌乳曲线的绘制与分析一、实验目的和要求通过累计各泌乳月份的产乳量及绘制泌乳曲线图，学习乳牛生产力水平的分析方法。

二、实验材料两头泌乳母牛各一个泌乳期的产乳原始记录。

（1）2011号奶牛第二个泌乳期的产乳量统计表（表1）（2）2015号奶牛第五个泌乳期的产乳量统计表（表2）（3）计算器三、实验方法及步骤1．根据表1、表2日产乳量统计表计算出全期实际产乳天数、实际产乳量、全期平均产乳量，并查出全期最高日产乳量。

2.累计各泌乳月的产乳量（自产犊日起，每30天为一个泌乳月份）。

3.计算出各泌乳月份的日平均产乳量（最后一个泌乳月不足30天时，按实际天数计算其平均值）。

4.用曲线绘制成图。

5.填写乳牛泌乳性能分析表。

6.对比两头牛的泌乳曲线和有关数据，分析其不同特点及具体的优缺点，以确定其生产力水平的高低。

表1 、2011号牛（中国荷斯坦牛）第二个泌乳期日产乳量统计表 kg41828211410741 5202823211410741 62024211491 7192491 8212824951 9222825951 10222725853 1122853 122323262013853 132******** 14242219138653 15252226138653 162514261913863 17251913865 182514251913865 192512251912620251225191262126121812622131862312241862414241811252714241242表2 、2015号牛（中国荷斯坦牛）第二个泌乳期日产乳量统计表 kg表3乳牛泌乳性能分析表表4各泌乳月产乳量累计（kg ）55平均日产乳量504540353252151（kg）5123456789101112泌乳月份泌乳曲线图四、作业根据表1、表2数据绘制泌乳曲线及分析泌乳性能报告。

实验五牛日粮配方的设计一、目的要求熟悉牛精料补充料配制原则与方法，初步掌握不同泌乳阶段奶牛TMR(日粮)和肉牛日粮配制技术。

荷斯坦奶牛泌乳曲线与生产周期阶段划分52AnimalScience&V eterinaryMedicineV o1.18No.2(TotalNo.74)荷斯坦奶牛泌乳曲线与生产周期阶段划分苏建华(扬州大学实验农牧场,江苏扬州225009)摘要以353头泌乳期产奶超过270d以上的荷斯坦奶牛产奶量记录为依据,按5产奶量百分比间隔,由微机分别逐头计算45～100产奶量所达到的泌乳天数,泌乳持续时间,求出平均值,用推导方式描述奶牛泌乳曲线图.根据泌乳曲线,建议将奶牛生产周期各阶段划分为:围产期即产前,产后各3周,泌乳盛期从泌乳第4周开始至17周;泌乳中期由I8周至3O周;泌乳后期由31周至于奶;干奶期至少需5周.关t词:荷斯坦奶牛;泌乳曲线;奶牛生产周期中圈分类号;$823.91文献标识码:B文章编号:1008—9381(2001)O2—0052—02荷斯坦奶牛从产犊泌乳开始到泌乳高峰,再逐步用中文VISUALF0XPR06.0自行编程,由计降低直至干奶,这个过程是奶牛生产的一个周期.了算机按所给程序自动对353头奶牛每日的产奶量分解奶牛泌乳过程的产奶规律,合理划分生产周期各个析,求出每头奶牛的最高日产奶量.从最高日产奶量阶段,将有利于针对奶牛的泌乳和生产性能,有目的的45比例开始,间隔5产奶量百分比,分别逐头地加强各阶段的饲养管理,采取有效的饲喂,防病等计算45～i00产奶量所达到的泌乳天数,泌乳持技术措施,充分发挥奶牛的泌乳潜力,促进整个泌乳续时间.剔除因病异常牛只的产奶量记录和两头极端期产奶量的提高.为此,笔者对奶牛泌乳过程中不同数据,求出平均值,用生物统计"从单组数据的样本估产奶量百分比情况下的产奶天数,持续时间等关系进计总体"的方法分析泌乳曲线各个统计点的数据,以行了初步探讨,现将有关情况简要介绍如下:推导方式描述泌乳曲线全部数据点.数据源奶牛的均1数据来源与处理方法值情况:头数,353头;平均胎次,2.85;平均最高日1.1数据来源产,36.65kg;平均达最高日产天数48.95d;平均数据资料来源于扬州大学实验农牧场牧业公司305d产量,7492.3kg.和扬州市扬农乳业有限公司等奶牛场,1998年1月2泌乳曲线数据与泌乳曲线圈至2000年1月泌乳期产奶超过270d以上的全部荷2.1泌乳曲线数据斯坦奶牛产奶量记录.按1.3方法收集的各种产奶量百分比例情况下1.2产奶量记录方式的产奶天数,持续时间的泌乳曲线统计点数据,见表采用每日记录法.日三次挤奶,并用电子称准确1.由表1推导出的泌乳曲线数据见表2.称重记载,累加三次奶量为每头奶牛的日产奶量,每2.2泌乳曲线图.日将各头牛的日产奶数据输入电子计算机的奶牛产利用产奶量百分比与产奶天数的关系数据,在计奶量数据表中.算机上应用MicrosoftGraph软件绘制出以下奶牛1.3产奶量记录的处理泌乳曲线图(图1).E要墅.注:X横座标表示产奶天数(d);Y纵座标表示产奶量百分比(oA)圈1奶牛泌乳曲线圈一'¨¨I-～牧产～,●,●～畜生～-t..^.....-',∞∞∞∞矗;柏∞幻∞o动物科学与动物医学2001年3月第18卷第2期(总第74期)53 表1泌乳曲线统计点数据产奶量百分比()到达该比例天数(d)95%置信区间本比例持续时间(d)95置信区间产奶量百分比(%)到达该比例天数(d)95置信区间本比例持续时间(d)95置信区间囊2泌乳曲线推导数据产奶量百分比(%)到达该比例天数产奶量百分比(%)刭达该比例天数产奶量百分比(%)刭达该比例天数3.744.535.396.578.8611.6115.O32O.14Z5.119095lOO9590858O757031.8140.7248.9559.0968.5089.07116.43139.93166.106560555O45189.49215.27231.77244.88258.443涝乳曲线分析从上图我们可以看出,荷斯坦奶牛产犊后1d的产奶量一般低于其最高日产量的45,具体比例与各奶牛场第一次是否挤净初奶的要求有关.产犊后第3～4d,产奶量为45左右;第20~-21d产奶量可迅速达到8OoA;第48-~50d产奶量达到最大值;维持2 ～4d后开始较大幅度降低,到116d后,产奶量回落到8O,以后平稳下降至190d时,产奶量已只有65,到260d时产奶量平均降低至45.这时,奶牛的产奶量变化较大,有的迅速降低或停奶,有的仍平稳产奶,也有的干奶期已到,人工进行干奶.所以,最后阶段的产奶量与奶牛的生产性能,怀孕,人工干奶等影响有关,数据的95置性区间相对误差较大.4讨论4.1在奶牛泌乳曲线研究中,自Gains(1927)提出产奶量数学模型以来,主要有wood模型,王一兰一一丁模型,改良分室模型等.这些研究均用数学模型以函数关系来预测不同胎次,不同季节产犊的奶牛个体产奶量.笔者试图用推导方式描述泌乳曲线图型, 并引入产奶量百分比的概念来揭示荷斯坦奶牛群体产奶量与时间关系的一般规律.4.2从产犊开始到再次产犊是一个生产周期.根据中华人民共和国专业标准一《高产奶牛饲养管理规范》(ZBB43002—1985)规定:将一个生产周期划分为:围产期,产奶盛期,产奶中期,产奶后期,干奶期等五个阶段.一般认为,围产期是指母牛分娩前后各15 d以内的时间;泌乳盛期是指母牛分娩16d到100d 以内l泌乳中期是指产后101d~-200dl泌乳后期是指201d至于奶以前的阶段l干奶期是指停止挤奶到分娩前15d的一段时向.由于荷斯坦奶牛产前产后生理状况和代谢的特殊性,需要精心的护理,笔者完全赞同不少学者提出将产前三周与产后三周鲻分为围产前期和围产后期的观点.从泌乳曲线图上我们可以看出,产后21d～119d时,产奶量在最高日产奶量的8O%以上,一般又可以把这一时期称为泌乳高峰期,我们将这14周划分为泌乳盛期是名列其实的. 产后120d～210d内产奶量由8O逐步降低到6O,将这13周定义为泌乳中期l把这以后到干奶前的时间定义为泌乳后期.另外,将干奶期控{II在5周时间是必须的,这有利于母牛的体况恢复,营养蓄积, 胎儿生长,特别是对乳腺细胞的恢复生长起到了较大的作用,为下一个生产周期打下了物质基础.因此,根据泌乳曲线图以产奶量百分比划分生产周期的各个阶段是完全合理的.4.3就奶牛各个个体的泌乳曲线来说,并非完全相同,这与产犊胎次,产犊季节,产奶性能等都有着密切的联系,本次所绘奶牛泌乳曲线图只是为了揭示奶牛整体泌乳情况的普遍规律,用以合理划分奶牛生产周期的各个阶段.虽然数据的运算量较大,但数据的样本数还受到一定的限制,有待今后进一步论证.(收稿日期;2001—01—22)"曲o盯扎掘圳挑∞韶H一一一一一～一H鳃∞一～一∞一一一一盯玲蛇一一～踮～一一一0凹鹃朋moH=:三8撒捣毗蛇∞扎;2扎m。

实验四乳牛泌乳曲线的绘制与分析一、实验目的和要求通过累计各泌乳月份的产乳量及绘制泌乳曲线图，学习乳牛生产力水平的分析方法。

二、实验材料两头泌乳母牛各一个泌乳期的产乳原始记录。

（1）2011号奶牛第二个泌乳期的产乳量统计表（表1）（2）2015号奶牛第五个泌乳期的产乳量统计表（表2）（3）计算器三、实验方法及步骤1．根据表1、表2日产乳量统计表计算出全期实际产乳天数、实际产乳量、全期平均产乳量，并查出全期最高日产乳量。

2.累计各泌乳月的产乳量（自产犊日起，每30天为一个泌乳月份）。

3.计算出各泌乳月份的日平均产乳量（最后一个泌乳月不足30天时，按实际天数计算其平均值）。

4.用曲线绘制成图。

5.填写乳牛泌乳性能分析表。

6.对比两头牛的泌乳曲线和有关数据，分析其不同特点及具体的优缺点，以确定其生产力水平的高低。

1表3乳牛泌乳性能分析表5平均日产乳量544332211（kg）泌乳月份泌乳曲线图四、作业根据表1、表2数据绘制泌乳曲线及分析泌乳性能报告。

实验五牛日粮配方的设计一、目的要求熟悉牛精料补充料配制原则与方法，初步掌握不同泌乳阶段奶牛TMR(日粮)和肉牛日粮配制技术。

二、实习材料计算器、电脑、奶牛和肉牛饲养标准、常规饲料营养成分表。

三、实习方法与步骤(一)试差法设计平衡奶牛日粮试差法配合日粮的步骤如下：（1）首先根据不同品种、年龄、类型、生产水平等，参照适当的饲料标准，确定所需各种营养物质的数量或比率。

（2）选择当地的常用饲料并确定其数量。

（3）试配日粮按饲料成分计算出配合日粮的营养物质含量，并与饲养标准相比较，一般来说，首先考虑能量和蛋白质两项，然后再考虑其他。

（4）修正日粮，如果所配日粮与确定的需要量不符合时，则应渐增某种饲料的用量，以求最后与确定的需要量符合。

为体重为600 kg、日产奶量为30 kg、乳脂率为3.5％的奶牛配制日粮。

可用饲料为全株玉米青贮、羊草、玉米、麸皮、豆粕、棉子饼、磷酸氢钙、石粉、食盐等。

如何评价奶牛的产奶性能？
奶牛产奶性能的评价指标主要有产奶量、乳脂率、乳蛋白率、泌乳速度、泌乳持续性（泌乳曲线）、体细胞数等。

通过DHI测定，可以实现对奶牛生产性能比较完整而系统的评价。

DHI每月采奶样1次，通过DHI分析仪的测定和软件的统计分析，形成DHI报告，其内容主要包括生产性能测定月报表（个体牛成绩）、体细胞分布一览表、体细胞走势表、各泌乳阶段生产性能汇总表、泌乳曲线、牛群分布统计一览表等。

这些信息，可以用以指导牛群选育、防治乳房炎，改进饲养管理、提高生产性能等工作。

因此，作为奶牛遗传育种T作的基础，DHI测定是奶牛生产中需要重视的一项工作，推广应用DHI测定体系是大势所趋。

问题是如何克服DHI测定仪器昂贵且依赖进口、饲养场对DHI测定意义的认识不足等困难。

1。

奶牛的泌乳规律
一、母牛不同泌乳月的泌乳规律母牛在同一泌乳期内不同泌乳月的产奶量，常随其品种、个体、健康和年龄状况，以及饲养管理不同而异。

通常是奶牛分娩后，产奶量逐步上升，约在第2个或第3个泌乳月时达到高峰阶段，然后逐步下降，一般高产牛上升幅度大，曲线在高峰期平稳，下降缓慢，所以牛产后泌乳量能很快上升并在高峰期维持较长时间，总产奶量较高。

而那些即使具有**产奶量，但下降速度快、幅度大的牛总产奶量低。

二、不同年龄和胎次的奶牛泌乳规律母牛的产奶量随年龄和胎次的变化而不同。

一般初产母牛，由于身体尚未成熟，泌乳能力较差，产奶量低于高峰胎次产奶量的80%～85%，随着年龄、胎次的增加及身体的成熟，第三或第四胎便可达到**产奶胎次，有些牛在饲养管理较好、饲料全价的情况下，产奶高峰胎次可持续至8胎，但多数母牛在5～6胎时产奶量已下降。

三、不同挤奶时间和挤奶间隔时间的泌乳规律挤奶时间的不同也会影响母牛的产奶量。

母牛乳腺在夜间活动频繁，乳汁生成和分泌较多，所以早晨产奶量高于其他时间，而且早餐前挤奶（早上5~6时）较好。

挤奶间隔时间对产奶量也有一定影响。

实验证明，挤奶间隔时间相等，产奶**多，如每隔8小时挤奶一次，或每隔12小时挤奶一次。

2012-11-29 05:27 荷斯坦奶农俱乐部网访问量：427 [ 字号：大中小 ]酮含量下降，刺激泌乳，泌乳的维持通过挤奶刺激和激素调节两个方面来完成的。

挤奶可以促进催乳、ACTH、肾上腺糖皮质类固醇和催产素的分泌，从而维持泌乳过程。

一般情况下，由于产奶量随着时间的推移而呈现出规律性的变化，表现为在泌乳初期产奶量迅速增加，经过一段时间后达到泌乳高峰期，以后便逐渐下降，直至干奶，这种规律性即构成泌乳曲线。

如果以时间为横坐标，以产奶量为纵坐标，所得到的泌乳期产奶量随时间变化的曲线就叫做泌乳曲线，这是产奶量受体内神经激素调节的结果，也是反应奶牛泌乳情况既直观又方便的形式。

在研究泌乳曲线时主要应考虑准确的取样方法和简单易行、准确度高又有高度产奶遗传代表性的数学模型，并结合实践制定出有实际生产意义的选择指数，以提高选种的准确性。

本文以新疆呼图壁种牛场荷斯坦牛不同胎次各泌乳月平均产奶量为基础，利用Wilmink 模型、Nelder逆多项式模型、Wood不完全伽玛函数模型等三种泌乳曲线模型对各泌乳月平均产奶量进行了曲线拟合，以及为新疆呼图壁种牛场产奶量的预测、品种选育和遗传改良、提高饲养管理水平等方面提供参考。

一、材料与方法1.数据来源取自新疆呼图壁种牛场2004-2008年共1345头奶牛的泌乳记录。

为了使制定的校正系数尽可能反映客观实际，在Excel 表里对数据资料进行了如下处理：(1)删除泌乳月少于4个月的记录;(2)删除月产奶量小于100千克的记录;(3)删除母牛在泌乳期患病的记录。

将数据整理后最后得到13336个泌乳期记录，其中头胎记录数为整体的，30.4%，第二胎占29.4%，第三胎为15.8%，第四胎为11.1%，第四胎以上的为13.3%，见表1。

表1 泌乳期产奶量统计(单位：千克)2.分析方法将整理好的数据导入SAS(8.1)软件，利用PROC- NLIN分析过程中的DUD法来进行泌乳曲线拟合，估计3个模型中的参数，根据拟合度(R2)越大越好、残差平方和(RSS，Residual Sum of Squares)和估计误差标准误(SEE，Standardized Estimated Error)越小越好的原则确定最佳泌乳曲线模型。

实验四乳牛泌乳曲线的绘制与分析一、实验目的和要求通过累计各泌乳月份的产乳量及绘制泌乳曲线图，学习乳牛生产力水平的分析方法。

二、实验材料两头泌乳母牛各一个泌乳期的产乳原始记录。

（1）2011号奶牛第二个泌乳期的产乳量统计表（表1）（2）2015号奶牛第五个泌乳期的产乳量统计表（表2）（3）计算器三、实验方法及步骤1．根据表1、表2日产乳量统计表计算出全期实际产乳天数、实际产乳量、全期平均产乳量，并查出全期最高日产乳量。

2.累计各泌乳月的产乳量（自产犊日起，每30天为一个泌乳月份）。

3.计算出各泌乳月份的日平均产乳量（最后一个泌乳月不足30天时，按实际天数计算其平均值）。

4.用曲线绘制成图。

5.填写乳牛泌乳性能分析表。

6.对比两头牛的泌乳曲线和有关数据，分析其不同特点及具体的优缺点，以确定其生产力水平的高低。

表1 、2011号牛（中国荷斯坦牛）第二个泌乳期日产乳量统计表 kg112226.525.520.513.58 6.553 1223232620138 6.553 1323.52226.520138 6.553 14242227.519138653 1525222619.5138653 16251426191386 4.53 172514.525.51913865 2.5 182514251913865 2.5 1925122519127.56 4.5 2.5 2025122519127.56 4.5 2.5 21261224.518127.56 4.5 2.5 2225.51324.51810.57.56 4.5 2.5表2 、2015号牛（中国荷斯坦牛）第二个泌乳期日产乳量统计表 kg表3乳牛泌乳性能分析表表4各泌乳月产乳量累计（kg）5 5平均日产乳量54 5 4 0 3 5 3 0 2 5 2 0 15（kg）5123456789101112泌乳月份泌乳曲线图四、作业根据表1、表2数据绘制泌乳曲线及分析泌乳性能报告。

第29卷第1期黄牛杂志Vo1.29No.1 2003年1月Journal of Yellow Cattle Science Jan.2003文章编号:1001-9111(2003)01-0019-03高产荷斯坦奶牛泌乳曲线的特征及相关分析程郁昕1~陈宏权2~王治华1~许春香3(1.安徽技术师范学院动物科学系~安徽凤阳233100;2.安徽农业大学畜牧水产学院~安徽合肥230036;3.安徽淮南乳品公司奶牛三场~安徽淮南232000)摘要:选取285胎次第三~四胎305d产奶量超过8000kg的安徽荷斯坦牛的泌乳资料~利用Wood模型(Y t=at b-ct e)进行高产荷斯坦奶牛泌乳曲线特征及相关的分析O结果表明:按照达到泌乳高峰周时间的不同对选取奶牛进行分组~11组Wood模型拟合泌乳曲线的平均拟合精度(R2)为0.939(幅度为0.899-0.962)O随着达到泌乳高峰周时间的推迟~规模因子a具有逐渐降低的趋势;而产奶量上升率b~产奶量下降因子c~泌乳持久力表现与a相反的变化趋势O在相关分析中~MY(实际产奶量)与EMY(估计产奶量)~PY(高峰月产奶量)为强正相关~r(相关系数)分别为0.999和0.776~与持久力为中等正相关(r=0.379)~与a~b~c的相关程度较弱;持久力与MY~EMY均为中等正相关~r分别为0.379及0.382~与b~a~c的r分别是0.905~-0.875~0.731~与PY的相关程度较弱(r=-0.140)O关键词:高产荷斯坦奶牛;Wood模型;泌乳曲线;特征分析;相关分析中图分类号:S823文献标识码:A自Wood首次提出用不完全伽玛函数作为描述奶牛泌乳曲线的数学模型以来[1]~世界各国的畜牧科技工作者针对这个课题进行了大量的研究~相继提出了多种用于描述泌乳曲线的数学模型[2]O我国于1980年以来也陆续开展了有关此领域的探讨~其主要研究方向在于:针对我国奶牛业的发展状况探索新的泌乳曲线模型[3-5]~同一模型下不同取样方法的比较[6~7]~以及曲线参数遗传特性的分析[8~9]等O 而无论进行哪一种研究~其最根本的宗旨就是要找出泌乳曲线波型与泌乳量高低之间的内在联系~从而为奶牛的育种工作提供理论上的依据O本次研究的目的旨在利用安徽荷斯坦奶牛的泌乳资料~通过Wood模型对高产奶牛的泌乳曲线进行研究~探索高产荷斯坦奶牛泌乳曲线的特征O1材料和方法1.1样本来源及分组方法资料取自安徽省淮南乳品公司1998~2001年间285胎次第三~四胎305d产奶量超过8000kg 的荷斯坦牛的泌乳资料~每个胎次的产奶记录按周统计(kg/周)O根据达到泌乳高峰周的时间对所取产奶胎次进行分组~即4周前~4周~5周~ ~12周~12周之后共11组O1.z数学模型本次研究采用已被国内外学者公认为估计效果较好的Wood模型进行分析[1]~模型结构为:Y t=at b e-ct(1)式中Yt为第t泌乳周的平均日产奶量(kg)~a 为规模因子~b为产奶量上升率参数~c为产奶量下降因子~t是时间变量(t=1~2~ ~43.5周)~e是自然对数的底O上式所描述的曲线波型由a~b~c三个参数共同决定~改变参数值可使曲线波型符合个体不同波型的泌乳曲线O当然~事实上这三个参数可以由不同奶牛群体或奶牛个体在一定时期内的实际泌乳资料估计得到;同时~通过这三个参数可以进一步估计出一些次级参数O泌乳持久力:P=c-(b+1)高峰泌乳周:WPY=b/c高峰泌乳量:PY=a(b/c)b e-b305d估计泌乳量:EMY=Zt1Y t收稿日期:2002-07-03基金项目:安徽省教育厅自然科学基金项目(2000JL057)资助.作者简介:程郁昕(1969-)~女~安徽技术师范学院动物科学系讲师~研究方向为数量遗传和动物育种.1.3模型参数的估计及拟合效果的分析估计模型参数时可先将其转变为线性形式,即:Ln(Y t)=Ln(a)+bLn(t)-ct(2)令Ln(a)=A,Ln(t)=Bi,Ln(Y t)=Y i为了使(2)式的估计误差最小,根据最小二乘原理,可使误差最小的a~b~c值作为估计参数,即使下式中的O值为最小:Z[Y t-(A+bB i-ct)]2=O也就是说,要使O对a~b~c的偏导数为零:8O 8A =O,8O8b=O,8O8c=O用矩阵方法稍加整理得:F L1JAbcn Z B i-Z t iZ B i Z B2i Z t i B iZ t i Z t i B i-Z tFL1J2i-1Z Y tZ Y i B iZ Y i tFL1J i求出A~b~c后,进一步对A值作指数转换求出a值G利用EXCEL程序进一步分析Wood模型对高产荷斯坦奶牛泌乳曲线的拟合效果,得出MY与EMY的相关指数(R2)G2结果与分析2.1拟合效果按不同泌乳高峰周分组估计的Wood模型参数及拟合效果见表1G就拟合效果而言,11组的平均R2值为O.939(幅度为O.899~O.962)G王治华等[1O]研究报道,当按照公牛家系分组利用Wood模型估计北京荷斯坦奶牛泌乳曲线参数时,22个公牛家系的平均R2值为O.943(幅度为O.854~O.983)G郭智群等[4]按间隔3O d取样,Wood模型拟合北京荷斯坦牛泌乳曲线参数的R2值为O.942G徐慧如等[6]曾经按照不同胎次分组并利用Wood模型进行参数估计,结果表明R2均值的变异幅度为O.93~O.99之间G上述结果与本文对安徽荷斯坦牛的参数估计结果十分相近G2.2高产奶牛泌乳曲线的特征由表1可见,所有高产奶牛几乎均在产犊12周之内达到了泌乳高峰,就MY和EMY而言,似乎从表观上看不出高峰周与泌乳量高低之间的内在联系,但可发现约2/3(65%)高产牛的高峰泌乳周集中在6~1O周G进一步对曲线参数分析不难看出,随着到达泌乳高峰周时间的推迟,规模因子a具有逐渐降低的趋势;而参数b则随着高峰周的推迟而逐渐增加,意味着奶牛到达泌乳高峰越早,产奶量上升的速度越慢,反之亦然;参数c同样随着高峰周的推迟而逐渐上升,但由于它在模型中以负指数的形式出现,所以不难推知,其实际涵义是奶牛到达泌乳高峰的时间越短,产奶量下降的速率越慢G泌乳持久力具有与参数b相似的变化趋势G由上述分析可见,对于高产奶牛而言,若到达泌乳高峰的时间较早,则往往具有较高的规模因子a,但同时伴随着产奶量下降较快~泌乳持久力较低的不足;若达到泌乳高峰的时间较晚,则往往具有较高的泌乳持久力,但同时伴随着规模因子a较低的缺陷G此结果预示着这些高产奶牛的泌乳量有着进一步提高的遗传潜力G表1高产荷斯坦牛泌乳曲线的参数估计高峰周n a b c MY(kg)EMY(kg)R2PY(kg)持久力4周前2643.6O.O19O.O189591.239597.2O O.92342.8259.96 4周154O.6O.112O.O289424.5O9428.23O.95143.4453.3O 5周1636.8O.145O.O289324.2O19329.52O.9434O.4159.98 6周4334.6O.175O.O299318.429322.49O.96239.7864.O7 7周2231.6O.267O.O389149.379152.12O.9574O.7263.O1 8周4328.7O.295O.O379132.289136.55O.96O39.4271.48 9周3624.3O.425O.O479O44.119O49.3O O.9424O.5O78.O3 1O周4121.6O.429O.O438788.2O8796.O7O.88937.7389.7O 11周2318.9O.488O.O458689.248697.44O.9O137.131OO.93 12周1616.7O.589O.O499374.129379.O8O.9574O.O912O.58 12周后415.3O.63O O.O461O258.211O265.86O.94442.38151.25 X28.5O.325O.O379281.289286.71O.9394O.4O82.94Shanks等和王治华等分别研究报道,荷斯坦奶牛的泌乳曲线波型具有遗传性,其中高峰泌乳量遗传力(h2)的变化幅度为(O.23O.O2)~(O.43 O.O1);高峰时期的h2为(O.18O.O1)~(O.27O.O2);参数a~b~c的h2分别为(O.O9O.O1)~ (O.27O.O1)~(O.O6O.O1)~(O.17O.O2)~ (O.17 O.O1)~(O.23 O.O2)[8~11]G可以设想,鉴于曲线波型具有遗传性,考虑选取高峰期较早的高O2黄牛杂志第29卷产家系与高峰期较迟的高产家系进行交配 然后通过选择 可望培育出产奶量更高~泌乳曲线波型更为理想的奶牛品系2.3相关分析高产荷斯坦奶牛泌乳曲线参数的相关系数见表2表2高产荷斯坦牛泌乳曲线参数的相关系数abcMY EMY PY 持久力a 1-0.989-0.9460.0670.0640.516-0.875b 10.9510.0290.031-0.3960.905c 1-0.192-0.191-0.4830.731MY 10.9990.7760.379EMY 10.7740.382PY 1-0.140持久力1由表2可见 MY 与EMY ~PY 为强正相关 r 分别为0.999和0.776 与持久力为中等正相关(r =0.379) 与参数a ~b ~c 表现出较弱的相关关系 持久力与MY 呈中等正相关(r =0.379) 与b ~a ~c 的相关程度较强 r 分别达到0.905~-0.875~0.731 与PY 表现为弱负相关(r =-0.140) 这些结果与陈幼春[3]~王治华[8]~Shanks [11]等的研究结果有一定的差异 可能是受到样本来源等许多因素的限制 有关泌乳曲线相关性的研究尚需进一步的深入 参考文献:[1]Wood P D P .Algebraic model of lactation curve in cattle [J ].Nature 1967 216(10):164-165.[2]王治华.泌乳曲线数学模型的研究进展[J ].中国奶牛 1993 (4):44-46.[3]陈幼春 等.利用伽玛曲线估测奶牛产奶量的研究[J ].畜牧兽医学报 1988 19(2):80-90.[4]郭智群 等.奶牛泌乳曲线数学模型拟合和早期预报[J ].畜牧兽医学报 1990 21(2):115-120.[5]储明星 等.Wood 模型和多元线性回归对奶牛个体泌乳曲线的拟合~预测及分析[J ].中国奶牛 1991 (2):29-31.[6]徐慧如 等.北京黑白花奶牛泌乳曲线的分析及应用[J ].中国畜牧杂志 1988 24(4):33-35.[7]王金良 等.用不完全伽玛函数泌乳曲线估计奶牛产奶量[J ].中国奶牛 1990 (5):31-34.[8]王治华 等.黑白花奶牛泌乳曲线参数的方差组分和遗传特性分析[J ].畜牧兽医学报 1994 25(5):417-422.[9]郭智群 等.奶牛泌乳曲线数学模型的遗传分析[J ].畜牧兽医学报 1993 24(2):136-142.[10]王治华 等.公牛家系泌乳曲线的研究[J ].黄牛杂志 1994 20(2):8-13.[11]Shanks R et al .Genetic aspects of lactation curves [J ].J .Dairy Sci . 1981 64:1852-1860.The Character and Correlation analysis on lactation CurVe ofHigh -Yield Holstein CowsG~ENG Yu -xin 1 C~EN ~ong -guan 2 WANG Zhi -hua 1 XU Chun -xiang3(1.An z Tec nz cCl N07?Cl C0llege Fe n gyC n g An z 233100 C zn C ;2.An z A g7z c lt 7Cl U nz ze7s z ty Hefe z An z 230036 C zn C ;3.An z H C zn C n DC z 7y C0?PC n y H C zn C n An z 232000 C zn C )a b stract :T he lactation period records of 285parities of Anhui ~olstein co W s W hose 305-day milk yield (the third parity and the fourth parity )W ere beyond 8000kg W ere selected to characteristically and correlatively analysis lactation curve of high -producing ~olsteins based on Wood model (Yt =atbe -ct ).T he results sho W ed that the average fitting accuracy (R 2)of lactation curve W as 0.939(the range W as 0.889~0.962)based on Wood model .With the postponement of attaining peak W eek of lactation parameter a W ent do W n gradually ;parameters b and c and lactation persistency W ent up little by little .T he correlation analysis sho W ed that milk yield (MY )had high positive correlation W ith estimated milk yield (EMY )and peak yield (PY ) correlation coefficient (r )W as 0.999and 0.776respectively ;had medium positive correlation W ith persistency (r =0.379);and had lo W correlation W ith parameters a b c respectively .Persistency had medium positive correlation W ith MY (r =0.379)and EMY (r =0.382);had high correlation W ith parameters b (r =0.905) a (r =-0.875)and c (r =0.731)respectively ;and had lo W negative correlation W ith PY (r =-0.140).K ey words :~igh -producing ~olstein co W s ;Wood model ;L actation curve ;Character analysis ;Correlation analysis12第1期程郁昕等:高产荷斯坦奶牛泌乳曲线的特征及相关分析。