中图分类号：S822 文献标志码：A 文章编号：1001-4330（2025）02-0484-10
0 引言
【研究意义】驴奶具有独有的营养成分和药用价值[1]，驴乳具有低致敏性[2]。驴乳汁较其他家畜乳更清、更稀，虽然产量不及奶牛，但与其他奶源相比，驴乳含有更多的溶菌酶和免疫球蛋白[3]，其酪蛋白与乳清的比例较低似乎是快速消化和低过敏特征的主要原因[4]。驴乳是一种免疫增强剂[5]驴乳可以作为一种有益的饮品成分来研究[6-7] 。而在驴乳中分化出的粘膜菌株SRV5和SRV10已证明了2种菌株的降胆固醇能力[8]。驴奶产量有限，但有独特的品质和药用价值[9]。驴奶富含蛋白质、钙、铁以及多种维生素。此外，驴奶还被用于制作高端的化妆品和保健品[10]。因此新疆驴泌乳方面的选种选育尤为必要。【前人研究进展】史鹏等[]通过多种模型对南疆驴的泌乳数据进行拟合发现，六次多项式模型的拟合效果最佳，且新疆的乳高峰期为第16周左右。李金霞等[的研究结果表明，高产驴与低产驴的泌乳曲线走势大致相同。【本研究切入点】新疆驴泌乳曲线的相关研究基本局限于同一地区，但对东疆地区新疆驴和南疆地区新疆驴不同饲养管理方式差异方面鲜有研究。【拟解决的关键问题】研究新疆东疆地区新疆驴和南疆地区新疆驴全泌乳期日产奶量，测定其相关体尺性状指标，分析东疆地区新疆驴和南疆地区新疆驴产奶量及泌乳规律差异，以及东疆地区半舍饲半放牧新疆驴产奶量与体尺指标间的关系，为生产中不同生产方向泌乳驴的选育提供参考。
材料与方法
1.1 材料
驴选自新疆哈密市伊吾县山南开发区塞外骞驴养殖合作社65头，南疆岳普湖驴场的新疆驴70头，选取产驹1个月之后正常泌乳健康新疆驴作为研究对象，收集其自产驹1月之后的150d早中晚3次的产奶量，每3d测定1次得到全期实际产奶量。精料组成成分为玉米秸秆、豆粕、麸皮和玉米等。每日07：00、13：00、19：00（依据季节调整挤奶时间前后不超 
）由挤奶工使用机器依次挤奶，并对称重记录。挤奶结束将母驴和驴驹混圈2h，刺激母驴产奶，然后将驴驹分圈饲养。所有驴进行统一的饲养管理，自由饮水，
1.2 方法
在挤奶前待测泌乳母驴保定于围栏内，保持水平站立，使用不同工具分别测量体高、体长、管围、胸围、胸宽、前段、尻宽度、尻角度、腰强度、乳房深度、乳房附着高度、乳头长度和乳头直径13个体尺数据。并收集其自产驹1月之后的 1 5 0 d 早中晚3次的产奶量，每3d测定1次得到全期实际产奶量。
1.3 数据处理
数据采用Excel2010和SPSS21.0进行处理，利用GraphpadPrism8.4软件中非线性回归模块套入六次多项式式对不同地区及不同胎次的泌乳驴进行产奶量的拟合，并比较其差异。
使用ANOVA方差分析进行方差齐性检验，Pearson相关系数判断13个变量之间的相关性。
通过对体尺数据与产奶量之间双变量相关分析，判断不同泌乳水平奶驴体尺指标之间的差异；通过对体尺数据进行回归分析，建立对应泌乳阶段的最优预测模型。
2 结果与分析
2.1 泌乳曲线拟合效果
研究表明，模型I为逆多项式模型，模型公式为 
;模型 I 为立方模型公式为 
;模型Ⅲ为六次多项式模型其模型公式为 
。其中 Y 为第 X d 的产奶量， X 为产奶天数， a， b， c， d ， e ， f ， g 为模型参数系数，其中拟合效果最佳的为六次多项式模型。第3胎次泌乳曲线拟合优度 
为0.8882；第4胎次泌乳曲线拟合优度 
；第5胎次泌乳曲线 
0.946 7。
3个模型对南疆地区新疆驴与东疆地区新疆驴拟合效果为模型 Igt; 模型 Igt; 模型I，故六次多项式模型对生产数据的拟合效果最佳，在后续不同胎次泌乳曲线拟合中选取六次多项式模型作为最佳拟合模型。表 1 ～ 3
表1新疆驴逆多项式模型拟合注： Y 是泌乳期第 X 天的测定日产奶量， X 为时间变量为泌乳日； 
为拟合度， S S E 为剩余平方和， R M S E 为剩余标准差，下同 Notes： Y is the measured dailymilk production on the X thdayof lactation，the timevariable X is the lactation day， 
is thedegree of fit，SSEis heresidual sumof squares，RMSE is the remaining standard deviation，the same asbelow
表2 新疆驴立方模型拟合Tab.2 Fitting of cubemodelsof Xinjiang Donkeys2.2 新疆驴最优拟合泌乳量变化
研究表明，东疆地区新疆驴的拟合优度为0.973，南疆地区新疆驴泌乳曲线的拟合优度为0.942，并绘制曲线。表4
六次多项式模型拟合南疆地区新疆驴与东疆地区新疆驴的泌乳曲线，从泌乳第30d开始记录，东疆地区新疆驴泌乳前期泌乳量上升速率相似，但东疆地区新疆驴上升起点更为提前一些。
且泌乳高峰的持续时间，东疆地新疆区驴相比南疆地区新疆驴更佳，在到达泌乳高峰之后，两者的
下降速率相似。图1
表3 东疆地区和南疆地区新疆驴六次多项式模型拟合
表4东疆地区新疆驴和南疆地区新疆驴泌乳日均量统计2.4 东疆地区新疆驴3、4和5胎次泌乳曲线差异
研究表明，在胎次为3的新疆驴，东疆地区新疆驴到达泌乳高峰的时间与南疆地区新疆驴相比更短，在到达泌乳高峰后，两者下降速率相似。在高峰期东疆地区新疆驴曲线亦宽于南疆地区新疆驴，东疆地区新疆驴产奶高峰的维持时间比南疆地区新疆驴更久。表6，图3
东疆地区新疆驴日均产奶量高于南疆地区新疆驴。图2
2.3 南疆地区新疆驴3、4和5胎次泌乳曲线差异
研究表明，南疆地区新疆驴3、4和5胎次的拟合优度分别为0.917、0.895和0.925；东疆地区
泌乳前期，南疆地区和东疆地区新疆驴第4胎次的新疆驴前期泌乳水平的差距并不明显，在泌乳量增加速率方面东疆地区新疆驴的4胎次高于南疆地区。泌乳高峰期的奶产量东疆地区新疆驴的4胎次新疆驴高于南疆地区新疆驴，再到达泌乳高峰期后两者下降速率相等。图4
南疆地区新疆驴和东疆地区新疆驴5胎次泌乳曲线总体均呈先缓慢上升后迅速下降的趋势，且二者均在100d左右到达泌乳高峰期，高峰期产奶量东疆地区新疆驴高于南疆地区新疆驴。到达泌乳高峰期后两者均缓慢下降，且下降速率相似。图5
新疆驴的泌乳规律相似，从产驹后1月开始记录，虽然地区不同，但泌乳曲线呈现相同的特点，在整个哺乳期间，呈缓慢上升后迅速下降趋势，均在产驹后第100d左右达到泌乳高峰，且3个胎次中泌乳水平由高到低依次是第3、第4、第5胎次。虽然地区不同，但第3胎次的新疆驴泌乳水平较高。但是也存在一定的差异，东疆地区（半舍饲）新疆驴泌乳初期泌乳水平低于南疆地区（全舍饲）新疆驴，但在泌乳中期东疆半舍饲半放牧的新疆驴泌乳水平高于南疆地区新疆驴。在达到泌乳高峰后南疆地区新疆驴和东疆地区新疆驴泌乳曲线均呈迅速下降趋势。表7
表7 南疆地区新疆驴泌乳量统计东疆地区新疆驴第3胎次的泌乳量差异显著高于南疆地区新疆驴，但东疆地区和南疆地区4、5胎次的新疆驴泌乳量差异不显著。半舍饲半放牧的饲养方式对产奶量的影响并不大，反而在改善动物饲养成本降低有显著作用。表8，图6
2.5 体尺数据与产奶量
2.5.1 体尺数据与矫正产奶量的描述性统计
研究表明，利用Excel软件对新疆驴的产奶量数据与泌乳量数据数据进行整理，计算平均值和标准差以及最大值最小值。表9
表8 东疆地区新疆驴泌乳量统计
表9基础参数2.5.2 泌乳量与体尺数据相关性
研究表明，泌乳量与胸围、胸宽、尻宽度、腰强度、乳房深度、乳头长度之间呈极显著正相关（ Plt; 0 . 0 1 ），与体高、体长、管围、前段、尻角度、乳头 直径相关性不显著。表10
表10 参数估计值2.5.3 体尺数和泌乳量回归
研究表明，其模型方程为 
一
。
建立的多元线性回归模型不显著，对模型进行优化，依据参数估计值将对泌乳量影响小的相关体尺指标进行剔除，选择胸围、胸宽、尻宽度、腰强度、乳房深度、乳房附着高度、乳头长度6个差异显著的体尺指标与月均产奶量进行线性回归，最优线性回归方程为Y=-9.83+0.42X，-
。表 1 1 ～ 1 3
最优线性回归方程预测泌乳量的拟合度测试，其 
，体尺指标估测产奶量的准确度可达0.939。表14
3讨论
3.1 不同地区驴泌乳曲线差异
Muhatai[13]研究表明，冬秋季产驴驹的日常产奶量高于夏春季产驴驹，且驴的泌乳曲线的一般模型与大部分反刍动物的标准模型大致相同。 
Alessandro等[14]对意大利驴的研究表明，Wood模型能够从汇总数据中识别出密度曲线的形状（ 
= 0 . 1 1 ， P lt; 0 . 0 1 ⟩ ，泌乳曲线高峰出现在 4 8 d 。栗敏杰等[15]在对荷斯坦奶牛的研究中发现，AS模型为1、3胎最佳拟合模型，wood模型为2胎牛最佳拟合模型，头胎牛和泌乳牛的泌乳曲线差异较大。贾先波等在对中国荷斯坦牛泌乳曲线数学模型拟合研究中，分别使用wood模型、Ali-Schaeffer模型、六次多项式模型等分别进行拟合后发现六次多项式模型拟合效果最佳，试验结果与其结果一致。
史鹏等[17]在对新疆驴泌乳曲线的研究表明，新疆驴在第5＼～15周产奶量随时间变化呈缓慢上升趋势，在15周到第16周左右达到泌乳高峰，在第16周以后呈迅速下降趋势，且六次多项式模型为最佳拟合模型，对2、3、4胎次拟合优度分别为0 . 9 1 9 ， 0 . 8 5 8 ， 0 . 9 2 0 ，且刘玲玲等[18]对新疆哈萨克马泌乳曲线拟合的研究结果也表明，六次多项式模型拟合效果最为理想。试验研究与其结果一致，六次多项式模型为拟合3、4和5胎次的泌乳驴最优数学模型，其拟合优度分别为0.9467、0.9761和0.8884。泌乳曲线呈先缓慢上升后迅速下降的趋势，在其开始泌乳100d左右达到泌乳高峰，且泌乳量水平越高其后期下降趋势越迅速。与史鹏等的研究结果相似但稍有不同，可能由于品种相同但饲养管理方式和环境略有差异导致。
rorresernerareen（ee erreerrenrta
表12 多元回归方程显著性检验 Tab.12 Significancetestofmultiple regressionequation
表13 多元回归方程显著性检验
表14 最优回归方程拟合优度李金霞等[12]的研究也表明，疆岳驴的泌乳高峰期为第2～4月，在达到第5个泌乳月后大幅下降，疆岳驴产乳量个体差异大。研究中泌乳高峰期为3月左右与其结果相似，东疆地区新疆驴在泌乳期长度以及泌乳高峰期峰值产奶量方面优于南疆地区新疆驴。
3.2 体尺与乳产量相关性
陈斌等[19的研究表明，南疆地区新疆驴的胸围、前段、尻宽度、尻角度、乳房深度、乳头直径与泌乳量呈显著正相关（ 
），泌乳量与腰强度、乳房附着高度、乳头长度呈极显著正相关，与其他性状相关性不显著。研究中泌乳量与胸围、胸宽、尻宽度、腰强度、乳房深度、乳房附着高度之间呈极显著正相关（ 
），与其他性状相关性不显著，研究结果与其存在一定差异，可能是饲养管理方式存在一定差异导致。
刘玲玲[20在对新疆伊犁马的研究中发现，其胸围、尻角度、乳房围与产奶量之间呈极显著正相关，与研究结果相似，可以说明马属动物在泌乳特征方面存在相似性。赵开森等[21在对奶牛的研究中发现，胸围对产奶量有显著的相关关系，与试验研究结果一致，可能是物种间存在差异导致，但也受季节，饲养条件等因素影响。Lin，CY等[22]发现，在荷斯坦奶牛中更长更宽的奶牛其产奶量更多，且产奶量与乳头长度乳头宽度存在极显著正相关。
4结论
六次多项式模型对新疆驴泌乳曲线拟合效果最优，东疆地区新疆驴泌乳高峰更长，不仅表现出更长的泌乳高峰期，而且其泌乳峰值产奶量也明显高于南疆地区新疆驴。此外，在泌乳量的前期，东疆地区新疆驴表现出更快的泌乳量增速。东疆地区新疆驴半舍饲半放牧的新疆驴的泌乳表现优于南疆地区。东疆地区新疆驴体尺与产奶量的方程可作为放牧条件下驴泌乳水平的潜在模型之一。
参考文献（References）
[1]陈英，杨红．新疆驴产业发展现状与对策研究[J]．山西农 经，2022，（2）：154-156. CHENYing，YANGHong.Researchondevelopment statusand countermeasuresofdonkey industryin Xinjiang[J].Shanxi Agricultural Economy，2022，（2） ：154-156.
[2]范小雪，韦海涛，高兴明，等．特色乳的功能活性研究进展 [J].乳业科学与技术，2021，44（3）：31-36. FAN Xiaoue，WEI Haiao，GAO Xinging，et al.Research progresson functional activity ofcharacteristicmilk[J].Dairy Science and Technology，2019，44（3） ：31-36.
[3]Aspri M，Economou N，PapademasP.Donkeymilk：an overview on functionality，technology，and futureprospects[J]. FoodReviewsInternational，2017，33（3）：316-333.
「4]Chiofalo B.DugoP.Bonaccorsi IL.etal.Comparison of maior
lipid components in human and donkey milk：new perspectives for ahypoallergenic diet in humans[J]. Immunopharmacology and Immunotoxicology，2011， 3 3 （ 4 ） ： 6 3 3 - 6 4 4 ：
[5]Jirillo F，Jirillo E，Magrone T. Donkeys and goats milk consumption and benefits to human health with special reference to theinflammatory status[J].Current Pharmaceutical Design， 2010，16（7） ： 859-863.
[6]LiY，Ma QS，Liu GQ，et al.Effects of donkey milk on oxidative stress and inflammatory response[J].Journal ofFood Biochemistry，2022，46（4）：e13935.
[7]LiY，FanYM，Shaikh A S，etal.Dezhou donkey（Equus asinus）milk a potential treatment strategy for type 2 diabetes[J]. Journal ofEthnopharmacology，2020，246：112221.
[8]Rastogi S，Mittal V，Singh A. In vitro evaluation of probiotic potential and safety assessment of Lactobacillus mucosae strains isolated from donkeyslactation[J].ProbioticsandAntimicrobial Proteins，2020，12（3）：1045-1056.
[9]翟宝菊，肖辉，李莉．鲜驴奶体内抑瘤作用的实验研究［J]. 新疆医科大学学报，2011，34（3）：274-278. ZHAI Baoju，XIAO Hui，LILi.Experimental study on tumor inhibitionof fresh donkey milk in vivo[J].Journal of Xinjiang Medical University，201，34（3）：274278.
[10］杨浩峰．新疆驴奶开发利用的初步研究[D]．乌鲁木齐：新 疆医科大学，2006. YANG Haofeng. Preliminary study on development and utilization of donkey milk in Xinjiang[D]. Urumqi： Xinjiang Medical University，2006.
[11］史鹏，陈斌，何海迎，等．新疆驴不同胎次泌乳曲线的拟合 研究[J]．新疆农业大学学报，2022，45（6）：447-450. SHIPeng，CHEN Bin，HE Haiying，et al．Fitting analysis of lactationcurveof differentparity of Xinjiangdonkey[J]. Journal of Xinjiang Agricultural University，2022，45（6）： 447-450.
[12]李金霞，陆东林，吾布力·沙地克，等．疆岳驴泌乳曲线分 析[J]．新疆农业科学，2014，51（10）：1923-1927. LI Jinxia，LU Donglin，Wubuli Shadike，et al．Analysis of Xinjiang Jiangyue donkeys location curve[J]. Xinjiang Agricultural Sciences，2014，51（10）：1923-1927.
[13]Muhatai G，Cheng L，Rugoho I，et al.Effect of parity，milking timeandstage of lactation onmilkyieldof Jiangyue donkey （Equusasinus）inNorth West China[J].Journal ofDairy Research，2017， 8 4 （ 1 ） ： 2 3 - 2 6
[14]DAlessandro A G，Martemucci G.Lactation curve and ffects of milking regimen on milk yield and quality，and udder health in Martina Franca jennies （Equus asinus）[J]. Journal of Animal Science，2012，90（2）：669-681.
[15]栗敏杰，任小丽，闫磊，等．荷斯坦奶牛泌乳曲线拟合分析 [J]．家畜生态学报，2020，41（9）：34-38. LI Minjie，REN Xiaoli，YAN Lei，etal. Lactation curve fitting analysisofHolsteindairycows[J].Journal ofDomesticAnimal Ecology，2019，41（9） ：34 -38.
[16］贾先波，陈仕毅，王杰，等．中国荷斯坦牛泌乳曲线数学模 型拟合及应用[J]．南京农业大学学报，2016，39（1）：133－ 138. JIA Xianbo，CHENShiyi，WANG Jie，et al.Application of different modelsto the lactation curvesin Chinese Holstein[J]. Journal of Nanjing Agricultural University，2016，39（1）：13- 138.
[17］史鹏．新疆驴泌乳曲线的拟合研究[D]．乌鲁木齐：新疆农 业大学，2022. SHI Peng. Study on Fitting lactation curve of donkeys in Xinjiang [D].Urumqi：XinjiangAgricultural University，2022.
[18]刘玲玲，吕士鹏，方超，等．哈萨克马泌乳曲线的拟合研究 [J].中国畜牧杂志，2020，56（8）：91-94. LIU Linging，LYU Shieng，FANG Chao，et al. Study on fitting lactationcurveof Hasakma[J].Chinese Journal of Animal Science，2019，56（8） ：91-94.
[19]陈斌，何海迎，史鹏，等．新疆驴体型结构与泌乳性状的相 关性分析[J]．新疆农业大学学报，2021，44（6）：412-416. CHEN Bin，HE Haiying，SHI Peng，et al. Correlation analysis of body structure and lactation traits of donkeys in Xinjiang[J]. Journal of Xinjiang Agricultural University，2021，44（6） ：412- 416.
[20］刘玲玲．新疆伊犁马早期泌乳性状及体型外貌与泌乳性能 关系的研究[D]．乌鲁木齐：新疆农业大学，2014. LIULingling.Studyon the relationship between early lactation traits，body shapeand lactation performance of Yili horses in Xinjiang[D]. Urumqi： Xinjiang Agricultural University，2014.
[21］赵开森，王杰，孙文强，等．西门塔尔牛体尺指标与产奶量 的相关性分析[J]．四川农业大学学报，2023，41（5）：901- 905. ZHAO Kaisen，WANG Jie， SUN Wenqiang， et al. Correlation analysis betweenbodysize indexesand milk yield of Simmental cattle[J]．Journal of Sichuan Agricultural University，2023，41 （5）：901-905.
[22]Lin CY，LeeAJ，McAllisterAJ，et al.Intercorrelations among milk production traits and body and udder measurements in Holstein heifers1[J]. Journal of Dairy Science，1987，70（11）： 2385 -2393.
Comparative analysis of lactation patterns of Xinjiang Donkey in Eastern Xinjiang Region and Xinjiang Donkey in Southern Xinjiang Region
ZHANG Guowei，CHEN Bin，LIU Wujun （College of Animal Science， Xinjiang Agricultural University，Urumqi 83Oo52，China）
Abstract：【Objective】To study the difference of milk yield and lactation regularity of Comparative analysis of lactation patterns of Xinjiang Donkey in Eastern Xinjiang Region and Xinjiang Donkey in Southern Xinjiang Region.【Methods】 In this study， 65 healthy Xinjiang Donkey in Eastern Xinjiang Region were selected from half house feeding and half grazing.7O farmed Xinjiang Donkey in Southorn Xinjiang whose lactation curve was fited by inverse polynomial model，cubic model and sixth -degree polynomial model to explore the diffrence of lactation regularity of donkeys in diferent areas of Xinjang.The optimal regression model was established basedon body size traitsand average monthly milk yield，which was usedas the prediction model of milk level under semi- grazing and semi-house feeding conditions.【Results】Theresults showed that the sixth-degree polynomial model had the best fiting effect on the lactation curve of Xinjiang donkeys.The lactation peak was longer in eastern Xinjiang，and the average daily milk yield and peak milk yield of Xinjiang Donkeyin Eastern Xinjiang Region were betterthan thatof Xinjiang Donkeyin Southorn Xinjiang.The regression model of body size and milk yield in eastern Xinjiang could be used as a prediction model of lactation level under semi - grazing and semi- house feeding condition【Conclusion】The lactation performance of semi - herding and semi-herding Xinjiang Donkey in Eastern Xinjiang Region is beter than thatof the donkeys in southern Xinjiang.
Key Words：Xinjiang Donkey in Eastern Xinjiang Region; Xinjiang Donkey in Southorn Xinjiang； lactationregularity；feeding and management；body measure
