本研究旨在基于机器人挤奶记录拟合荷斯坦牛各乳区的泌乳曲线,并分析泌乳曲线的群体规律。收集了来自北京地区某规模化牧场2022—2024年1 391头荷斯坦牛771 186条利拉伐挤奶机器人(VMS 2017)挤奶记录,首先根据每日多次挤奶记录整理获得各乳区的日产奶量,使用WOOD不完全伽马模型对奶牛各乳区泌乳曲线进行拟合,获得模型参数a、b、c以及计算得到泌乳持续力、高峰奶量、高峰日3个二级参数和305 d产奶量,并对各乳区泌乳曲线进行分析。采用SAS软件广义混合线性模型过程分析胎次和产犊季节对各乳区泌乳曲线的影响。结果表明:荷斯坦牛左前、右前、左后、右后乳区的日平均产奶量分别为8.53、8.73、10.82、10.65 kg,高峰奶量分别为10.41、10.56、12.86、12.51 kg,高峰日分别为65.45、67.98、75.39、76.57 d,305 d产奶量分别为2 578.85、2 626.33、3 225.13、3 133.45 kg,荷斯坦牛后乳区产奶量高于前乳区,前后乳区指数占比分别为45%与55%。胎次和产犊季节效应对泌乳曲线相关参数a、b、c以及泌乳持续力、高峰奶量、高峰日、305 d产奶量均有显著影响,头胎奶牛4个乳区高峰奶量与305 d产奶量均显著低于经产牛,且高峰日最晚到达,但泌乳持续力显著高于经产牛。左前与右前乳区高峰奶量,305 d产奶量分别在秋季和冬季达到最高,除左前乳区为春季产犊高峰日最晚,其余乳区均在秋季产犊时最晚达到产奶高峰。基于荷斯坦牛的机器人挤奶记录,本研究揭示了北京地区高产奶牛各乳区的泌乳曲线特征,为优化奶牛场饲养管理和牛群遗传改良提供了有用信息。