西农萨能奶山羊泌乳规律研究*

黄耀辉，王永军，胡秀芝

(西北农林科技大学 动物科技学院，陕西 杨凌 712100)



科学研究

西农萨能奶山羊泌乳规律研究*

黄耀辉，王永军*，胡秀芝

(西北农林科技大学 动物科技学院，陕西 杨凌 712100)

研究以50只同期产羔且具有完整的240 d泌乳记录的泌乳母羊为研究对象，采集其每天泌乳量，绘制泌乳曲线，利用wood模型对泌乳曲线进行拟合，对其泌乳规律进行了研究。结果表明：西农萨能奶山羊240 d平均泌乳量为(544.39±111.56)kg；泌乳高峰出现在泌乳第39天，其泌乳高峰值平均为2.76 kg/d；西农萨能奶山羊240 d泌乳曲线符合Wood模型，其泌乳曲线可用方程Y=1.6997x0.1885e-0.0048x(R2=0.9297)；胎次对参数c值、泌乳峰值、泌乳持续力和总泌乳量有显著影响(P<0.05)，其中第4胎次总泌乳量最大，为627.2 kg，而第2胎次最小，为464.1 kg，第4胎次母羊的泌乳峰值平均为3.25 kg/d，依次为5胎次及以上的2.86 kg/d，3胎次的2.75 kg/d及2胎次的2.32 kg/d；以第3胎次持续力最好为6.65，且泌乳持续力在第3胎以后有下降的趋势；在机械挤奶条件下，产羔数对泌乳曲线不存在显著影响(P>0.05)。综合分析，wood模型拟合的泌乳曲线能客观反映西农萨能奶山羊的泌乳规律。

西农萨能奶山羊；泌乳量；Wood模型；泌乳曲线

西农萨能奶山羊作为一个优质奶山羊品种，具有产奶量高、遗传性能稳定等特点，但是有关泌乳规律的详细研究报道很少。泌乳曲线能够客观、详尽的反映泌乳规律，已被大多数学者认可，而且能够反映泌乳曲线的数学模型不仅在奶牛上得到了广泛运用，在绵羊、山羊应用效果也较好。现有研究发现，依据家畜的泌乳规律对其实施科学化管理是充分发挥其泌乳性能的必然选择。目前关于泌乳规律的描述主要采用不同的泌乳曲线方程进行描述，其中Wood泌乳曲线方程是目前应用最为广泛的泌乳曲线方程之一。Wood模型估计的三个参数a、b、c与泌乳曲线特征有很大关联性，因此Wood模型成为了描述泌乳曲线运用最广泛的模型[1-3]。

因此本研究以西农萨能奶山羊泌乳母羊为研究主体，利用机器挤奶方式测定羊群泌乳期的泌乳量，并运用wood模型进行泌乳数据拟合，得到反映泌乳规律的泌乳曲线方程。通过wood方程推导的二级公式，与实际各值进行比较分析，揭示西农萨能奶山羊产后240 d的泌乳规律。这对于充分发挥母羊泌乳潜力、合理日常营养、指导羊群遗传选育等均具有重要意义。

1 材料与方法

1.1 试验动物饲养管理

泌乳母羊按照西北农林科技大学教学试验基地种羊场母羊饲养管理办法进行饲养与管理。羔羊生后第 2 日即与母羊隔离，之后母羊每天进行机械挤奶，每天早晚各挤奶一次。母羊日粮包括精料、青贮玉米和苜蓿草粉，自由饮水，并注意保持补饲料及饮水清洁。

1.2 数据采集与整理

本研究以240 d泌乳量记录完整的50只西农莎能奶山羊母羊作为研究对象，母羊养殖于西北农林科技大学教学试验基地种羊场，按照人为延长黄体期技术进行同期发情处理并配种。以2009年12月底顺产的母羊作为备选试验用羊，以每只母羊为一个重复进行每天泌乳量测定，连续测定240 d的泌乳量。其中二胎母羊10只，三胎母羊11只，4胎母羊9只，5胎以上母羊20只；产单羔母羊15只，双羔母羊23只，三羔以上母羊12只。

表1 试验用羊精料补充料配方及营养成分

注：①预混料为每千克饲粮提供：维生素A 25 000 IU，维生素D312 500 IU，维生素E 250 mg，钴3 mg，铜112.5 mg，碘6.5 mg，铁87.5 mg，锰132.5 mg，硒5.5 mg，锌170 mg。②营养水平为计算值。

Notes:①Premix provides in per kg of diet:vitamin A 25 000 IU，vitamin D312 500 IU，vitamin E 250 mg，Co 3 mg，Cu 112.5 mg，I 6.5 mg，Fe 87.5 mg，Mn 132.5 mg，Se 5.5 mg，Zn 170 mg.②Nutrient levels are calculated values.

1.3 分析统计指标

1.3.1 泌乳量 按下式进行计算

(1)

式中：Yi表示第i(i:1… ,50)母羊204 d泌乳量；t 是泌乳天数(t:1…,240)),Yit是第i母羊第t天的泌乳量。

1.3.2 泌乳曲线拟合 为了准确描述西农萨能奶山羊240 d 泌乳曲线，本研究选用Wood[4]泌乳曲线方程对母羊每天泌乳量进行了拟合。Wood于 1967 年提出动物的泌乳曲线方程的一般模式为：

Y=atbe(-ct)

(2)

式中:Y表示母羊产后第t天泌乳量，t 为母羊产后泌乳天数，a、b、c为描述泌乳曲线形状的3个特征参数：其中a值非常接近母羊产后第一天的泌乳量，因此也常被认为代表母羊产后第一天的泌乳量，b描述了泌乳高峰出现前泌乳曲线上升的速度，c描述了泌乳高峰出现后泌乳曲线下降的速度[5]。

通过对 Wood方程进行数学推导，在运用非线性回归分析以及最小二乘法，估计a，b，c三个系数的值。可获得预测峰值泌乳量Ymax(peak milk yield)和泌乳高峰出现日期Dmax(lactation peaked date，产后第几泌乳日)以及泌乳持续力P(persistence)的估计公式分别为

Ymax=a(b/c)be-b

(3)

Dmax=b/c

(4)

P=- (b+1) ln c

(5)

1.4 数据处理

试验数据用Excel 2007、SPSS 17.0等软件进行处理，用Duncan法进行多重比较，结果用“平均值±标准差”表示，P<0.05表示差异显著。

2 结果与分析

2.1 西农萨能奶山羊泌乳量测定结果

50只西农萨能奶山羊240 d泌乳量测定结果见表2。由表2可知，50只西农萨能奶山羊个体间泌乳量存在较大差异，其240 d泌乳量为247.1～786.6 kg，平均泌乳量为(544.39±111.56)kg。

2.2 西农萨能奶山羊泌乳曲线

根据50只奶山羊240 d泌乳期内每天泌乳量测定结果，绘制西农萨能奶山羊240 d 的泌乳曲线，如图1所示。由图1可知，西农萨能奶山羊泌乳曲线总体呈现泌乳前期泌乳量快速上升期、泌乳高峰期、到达泌乳高峰后泌乳量缓慢下降的变化趋势。

对图1所示泌乳曲线与Wood 泌乳曲线方程进行计算机拟合后所获得的描述泌乳曲线形状的3个特征参数a，b，c的值分别为1.6997±0.3098，0.1885±0.0457和0.0048±0.0011。据此，西农萨能奶山羊泌乳曲线可用如下方程描述：

Y=1.6997x0.1885e-0.0048x(R2=0.9297)

(6)

利用泌乳曲线方程可对母羊泌乳期间的每日泌乳量进行估计，将每日泌乳量的估计值与实测值进行配对t检验，结果差异不显著(P>0.05)。不仅如此，利用泌乳曲线方程还可进一步预测泌乳高峰出现日期、峰值泌乳量和总泌乳量，把a，b，c的值分别代入公式3和公式4，得出泌乳峰值及峰值出现的时间分别为2.81 kg和39.34 d，分别相当于实际观测值的101.81%和100.88%；240天母羊的总泌乳量预测值为544.02 kg，而实际值为544.43 kg，预测值略低于实际值，相当于实际观测值的99.93%。上述结果表明，西农萨能奶山羊的泌乳曲线完全符合Wood模型。这与Torres-Hemandez等[6]，Gipson等[7]和Williams等[8]的研究结果一致。

表2 西农萨能奶山羊240 d泌乳量测定结果

图1 西农萨能奶山羊产后240 d泌乳曲线

2.3 胎次对泌乳量的影响

50只母羊的不同胎次泌乳曲线参数表明：胎次对泌乳曲线有显著影响(P<0.05)(见表3和图2)。

表3 泌乳高峰出现日期、峰值泌乳量和总泌乳量预测结果与实际观测值的比较

表4 不同胎次泌乳参数

注:同行不同上标字母表示差异显著(P<0.05)。

Notes:varions superscript in the same line show significant difference(P<0.05).

图2 不同胎次羊泌乳曲线

和其他胎次相比，第4胎次的母羊参数a最高(2.0268)，这与赵志华等[9]的研究结果相一致；而第二胎次母羊的a值最低，为1.5341。另外，不同胎次的参数c值、泌乳峰值、持续力和总泌乳量差异显著(P<0.05)。

从图2可以看出，第4胎次的母羊泌乳量总是高于其他胎次的泌乳量，通过显著性检验我们还得出，第4胎次母羊的泌乳峰值明显高于其他胎次，但是其他胎次间的泌乳峰值差异不显著。第4胎次母羊的最大泌乳量平均为为3.25 kg/d，其他胎次泌乳峰值依次为5胎次及以上的2.86 kg/d，3胎次的2.75 kg/d及2胎次的2.32 kg/d。和泌乳峰值的结果相对应，240 d各胎次的泌乳总量有相似的趋势：第4胎次泌乳总量最大，为627.22 kg，而第2胎次最小，为464.06 kg。以上结果显示了第4胎次母羊泌乳性能优于其他胎次，而第2胎次相比其他胎次泌乳性能较低。Carnicella等[10]研究得出，山羊第3、第4胎次母羊产奶量明显高于(9%～14%)第1、第2胎次母羊。因此可推测第2胎次的母羊还不够成熟，因而泌乳性能还未达到最佳。

试验结果表明，第3胎次持续力最好，为6.65，且持续力在第3胎以后有下降的趋势(第4胎次6.31，5胎及以上6.27)。Portolano等[11]、Susmel[12]和Giaccone[13]研究表明，持续力和胎次呈负相关关系，即胎次越高，持续力越低，这与本研究结论相一致。

2.4 产羔数对泌乳量的影响

50只母羊的不同产羔数泌乳曲线参数表明：产羔数对泌乳曲线没有显著影响(P>0.05)(见表5和图3)。由表5可知：不同产羔数的参数a、b、c值、峰值时间、泌乳峰值、持续力和总泌乳量差异均不显著(P>0.05)。

产双羔或多羔的母羊与单羔母羊相比产奶量提高，这是很多学者均得出的结论[14-15]，多产的原因主要在于多羔母羊的乳腺排空频率和排空速度提高[16]。本研究中，产羔数对泌乳参数无明显影响(P>0.05，表5)，因此对泌乳曲线的形状也基本无影响(图3)。这与Carnicella等[17]的研究结果不一致，因为两个研究中母羊的挤奶方式不同，Carnicella等采用的是母羊与羔羊不隔离，自由哺乳，而本研究采用的是机械挤奶，母羊与羔羊隔离，因此母羊每日泌乳量与羔羊数无关。

表5 不同产羔数泌乳参数

图3 不同产羔数母羊的泌乳曲线

3 讨 论

通过西农萨能奶山羊与陕北白绒山羊、国外杂种羊、努比羊泌乳曲线方程的比较分析[18-20]，见表6，西农萨能奶山羊产后泌乳规律如下：(1)与陕北白绒山羊、国外杂种羊、努比羊比较而言，西农萨能奶山羊具有很好的泌乳性能，产后第一天泌乳量大(a值大)；(2)西农萨能奶山羊泌乳初期每日泌乳量上升的速度最慢(b 值最小)，但泌乳峰值最大，而且泌乳峰值出现的时间晚，在本试验研究中发现：西农萨能奶山羊泌乳峰值在产后39 天左右，而陕北白绒山羊泌乳峰值出现在母羊产后第21天左右，国外杂种羊、努比羊的泌乳峰值分别在产后第57天左右和第18～19天；(3)泌乳持续力强，泌乳后期日泌乳量下降最慢(c 值最小)。

西农萨能奶山羊与其他品种山羊泌乳曲线比较见图4，由图4可以知道，西农萨能奶山羊的泌乳曲线处在最上面，表明泌乳性能很好；泌乳曲线前期上升速度较快、泌乳高峰值最大且出现较迟、泌乳高峰期后曲线下降速度较慢，表明泌乳持续力强。出现这个结果的原因与本品种的长期选育有关系，众所周知，西农萨能奶山羊作为特有的奶山羊品种，其主要目的就是产奶，这就决定了在选育过程中以泌乳性能作为第一选育因素。一般高产母羊，都具有a值大，b值和c值都相对较小，即泌乳曲线起点高，上升和下降速度都较小，泌乳曲线高峰期维持时间长，西农萨能奶山羊泌乳方程也符合此特点。

表6 西农萨能奶山羊、陕北白绒山羊、国外杂种羊和努比羊泌乳曲线方程的比较

图4 西农萨能奶山羊和其它品种羊的泌乳曲线

4 结 论

(1)西农萨能奶山羊240 d平均泌乳量为(547.15±102.85)kg；泌乳高峰出现在泌乳第39天，其泌乳峰值平均为2.76 kg/d。

(2)西农萨能奶山羊240 d泌乳曲线符合Wood模型，其泌乳曲线可用方程Y=1.7095x0.1854e-0.0048x(R2=0.9297)。

(3)胎次对参数c值、泌乳峰值、持续力和总泌乳量有显著影响(P<0.05)，其中第四胎次总泌乳量最大，为623.5 kg，而第二胎次最小，为461.3 kg；第四胎次母羊的泌乳峰值平均为3.25 kg/d，依次为五胎次及以上的2.86 kg/d，三胎次的2.75 kg/d及二胎次的2.32 kg/d；以第三胎次持续力最好，为6.65，且持续力在第三胎以后有下降的趋势。

(4)在机械挤奶条件下，产羔数对泌乳曲线不存在显著影响(P>0.05)。

[1] Cappio-Borlino A,Pulina G,Rossi G.A non-linear modification of Wood's equation fitted to lactation curves of Sardinian dairy ewes[J]. Small Ruminant Research,1995,18(1):75-79.

[2] 张 荣,王小玲,王爱华.奶山羊泌乳性能的研究[J].畜禽业,2004(9):8-9.

[3] 高腾云,安立龙.两种杂种牛的产奶性能及泌乳曲线[J].畜牧兽医学报,2000,31(5):411-415.

[4] Wood P D P, Algebraic model of the lactation curve in cattle[J]. Nature,1967,216(10): 164-165.

[5] 周卫东.中国荷斯坦奶牛泌乳曲线的数学模型及其遗传特性研究[D].浙江 杭州:浙江大学,2001.

[6] Tones-Hemandez G，Hohenboken W D.Biometric properties of lactations in ewes raising single or twin lambs[J].Anim Prod,1980，30(3): 431-436.

[7] Gipson T A,Grossman M.Diphasic analysis of lactation curves in dairy goats[J]. J Dairy Sci,1989,72(4):1 035-1 044.

[8] Williams J C.An empirical model for the lactation curve of white British dairy goats[J]. Anim Prod，1993,57(1): 91-97.

[9] 赵智华,彭顺清,李凤翔,等.用Wood模型拟合四川地区黑白花奶牛泌乳曲线的研究[J].四川畜牧兽医,2000,27(113)：19-20.

[10] Carnicella D, Dario M, Ayres M C C.et al. The effect of diet, parity, year and number of kids on milk yield and milk composition in Maltese goat[J]. Small Rumin Res,2008，77(1)：71-74.

[11] Portolano B,Spatafora F,Bono G, et al.Application of the Wood model to lactation curves of Comisana sheep[J]. Small Ruminant Research,1997,24(1):7- 13.

[12] Susmel P, Spanghero M, Mills C R,et al. The application of the wood model to the lactation curves of Italian Simmental cows[J]. ASPA National Congress,1993,10: 161-168.

[13] Giaccone P,Portolano B,Bomumo A, et al. Milk production and quality in the Derivata di Siria goat[J]. Zoot Nutr Anim,1995,21: 97-109.

[14] Alexander G，Davies H L.Relationship of milk production to number of lambs born or suckled[J]. Australian J Agr Res,1959,10(5):720-724.

[15] Loerch S C,McClure K E,Parker C F. Effects of number of lambs suckled and supplemental protein source on lactating ewe performance[J].J Anim Sci,1985,60(1):6-13.

[16] Gardner R W,Hogue D E.Effects of energy intake and number of lambs suckled on milk yield, milk composition and energetic efficiency of lactating ewes[J]. J Anim Sci,1964,23(4):935-974.

[17] Carnicella D, Dario M, Ayres M C C ,et al. The effect of diet, parity, year and number of kids on milk yield and milk composition in Maltese goat[J]. Small Rumin Res,2008,77(1):71-74.

[18] 梁铁刚.陕北白绒山羊泌乳规律及哺乳期羔羊体重生长规律研究[D].陕西 杨凌:西北农林科技大学,2012.

[19] 朱仁俊, 楮 柯, 冷 静, 等.努比羊和云岭黑山羊初产泌乳性能的研究. 安徽农业科学,2010,38(6):2 956-2 958.

[20] Montaldo H,Almanza A,Juárez A.Genetic group, age and season effects on lactation curve shape in goats[J]. Small Ruminant Research,1997,24(3):195 - 202.

Study on Lactation Rules of Xinong Saanen Dairy Goats

HUANG Yao-hui1, WANG Yong-jun1*,HU Xiu-zhi1

(CollegeofAnimalScienceandTechnology,NorthwestA&FUniversity,Yangling,Shaanxi,712100)

To study the lactation rules of Xinong Saanen dairy goat,a total of 50 Xinong Saanen dairy goat ewes which were lambed at the same period and had a full lactation record of 240 days were chosen to collecte the data of daily milk yield of every ewe, draw the lactation curve, us the wood model to fit the lactation curve so as to research the rules of ewe lactation . The results showed that:the average milk yield of Xinong Saanen dairy goat in 240 days was (544.39±111.56)kg； the daily milk yield saw a peak on 39th day and the average lactation peak was 2.76 kg/d；Wood model for the lactation curve of dairy goat showed a reasonable fit to the milk production data of Saanen ewes, and the lactation curve of Xinong Saanen dairy goat was described by the equation below: Y=1.6997x0.1885e-0.0048x(R2=0.9297); There were a strong influence of parity on the parameter c, peak production, persistency and total milk yield (P<0.05) and the forth ewes had the highest milk yield 627.2 kg,while the second ewes had the least yield of 464.1 kg;the peak production of the forth lambing ewes on average was 3.25 kg/d,and that of 5th2.86 kg/d,3rdlawbing 2.75 kg/d lambing or more and ,2ndlawbing 2.32 kg/d ,respectively; the third lambing ewes had the best persistency,6.65, and persistency values tended to decrease after the third lambing;(4)in conditions of mechanical milking, the type of lambing did not significantly influence the shape of the lactation curve(P>0.05). In conclusion, the lactation curve fitting wood model can objectively reflect the lactation rules of Xinong Saanen dairy goats.

Xinong Saanen dairy goats; milk yield; wood model; lactation curve

2015-04-27

2015-05-13

公益性行业科研专项(K312020806)

黄耀辉(1990-)，男，陕西安康人，在读硕士，研究方向：动物营养与饲料科学。E-mail:huangyaohui0624@126.com

*[通讯作者] 王永军(1964-)，男，陕西礼泉人，副教授，硕士生导师，研究方向：动物营养与饲料科学。E-mail:dkxywyj@163.com

S811.6

A

1005-5228(2015)08-0011-06