[0001] 本申请要求2015年7月2日提交的中国专利申请201510383182.X的优先权，将其通过引用并入本文。

[0002] 本发明涉及畜牧养殖领域，具体的涉及一种牛体温的测量方法。

[0003] 体温是机体重要的生理指标，正常情况下，牛体温会在相对恒定的范围内，随生理活动变化而成一定规律性变化(Cooper-Prado et al.2011)，比如发情、排卵、妊娠及分娩等。然而，在牛体温检测的实践中，大部分牛场目前依然采用水银体温计进行直肠测温，该方法测温效率低下，易引起奶牛疾病交叉传播，且不能准确实时获取每头牛的体温数据，需要的人力也较多，耗时费力，远远不能满足现代规模化牧场的管理要求(朱翔宇，2013)。随着生物传感技术的发展，部分发达国家已将无线传感技术应用于奶牛体表测温，并取得了一定的研究成果。Bligh等1974年首次提出将无线遥测技术运用于动物体温测定的构想(Bligh et al.1974)；Kyle等与Suthar等(Kyle et al.1998；Suthar et al.2011)在奶牛阴道植入含温度采集、无线电接收设备的棒状无线遥测系统，对奶牛阴道温度进行实时检测；Raul等与Miranda等(Raul et al.2006；Miranda et al.2009)将与会阴肌肉组织无排斥反应的无线温度传感器，植入奶牛会阴部肌肉，检测会阴部温度，进行发情鉴定；AlZahal等与Edouard等(AlZahal et al.2009,2010；Edouard et al.2011)分别通过食道将无线电瘤胃丸投入瘤胃，检测奶牛瘤胃温度变化，结果显示，无线电瘤胃丸能够有效检测奶牛体内温度变化。

[0004] 上述方法虽不同程度地达到了测温目的，但无论瘤胃丸还是皮下埋置都会对奶牛造成一定伤害，均不符合动物福利要求。近年来，红外测温不断被用于猪体表温度测量，然而效果都不很理想，无法指导生产应用。现代医学的发展，推动了对人体体表温度变化规律的深入研究，并取得了一定阶段性成果，然而，牛体表被毛较多，无毛裸露的体表部位较少，并且检测很不方便，体表温度研究报道较少，仅见到关于热应激对奶牛体表温度的影响(温雅俐，2011)和乳区表面温度变化与乳房炎关系(Metzner et al.2014)等少量研究，这些结果只揭示了部分环境因素及体表个别部位温度与部分病理现象的关系，没有深入揭示体表温度与体温(肛温)间的关系及其变化规律，对牛体表温度变化规律缺乏足够深入系统的研究。对牛体表温度变化规律的研究需要对牛体温进行连续不间断测量，但是奶牛体表被厚厚毛发覆盖，无毛部位较少，且体表接触式测温受外界环境温度影响较大难以反映体温的真实状况，尤其是，检测装置难以固定的问题成为牛体温准确连续测量难以突破的瓶颈。因此，目前还没有找到一种准确的牛体温连续测量方法。

[0005] 本发明的目的在于克服现有技术中存在的缺陷，提供一种能够准确、连续的测量牛体表温度的方法。

[0006] 选择合适的牛体温检测位置对尽可能减少外界环境因素影响，高保真采集体表温度至关重要，而已有研究均没有找到特别合适的温度检测位置。检测位置不同主要关系到两个方面，一方面是关系到被检测部位的体温是否能在一定程度上反映被测动物的真实体温状态；另一方面是关系到该部位是否能够适合进行连续自动的检测，是否能够避免被测动物的活动所引起的体温检测的误差。

[0007] 本发明的发明人在对牛体温状况的长期深入研究过程中发现，牛体表温度与牛的实际体温(直肠温度)之间存在一定程度的差异，而且不同体表部位与真实温度之间的变化趋势也不一致。这为通过测量牛体表温度表征体温带来了难度。发明人经过大量摸索实验，提供了一种牛体温测量方法，该方法包括：步骤一：首先检测牛腿跖部温度和直肠温度，再将检测时间和季节影响作为固定效应，个体影响作为随机效应对牛腿跖部温度和直肠温度数据进行混合线性模型分析；

[0008] 步骤二：在需要测量牛体温时将测量获得的牛腿跖部温度带入所述数学模型计算获得牛直肠温度。

[0009] 可选的，所述牛腿跖部为牛后腿跖部。

[0010] 可选的，所述牛腿跖部为牛腿跖骨下部近籽骨上方。

[0011] 可选的，所述方法包括将体温检测装置固定于牛腿跖部，实时检测牛只温度数据。

[0012] 可选的，述方法包括利用体温检测装置每小时检测2-10次。

[0013] 可选的，所述数学模型是将检测时间和季节影响作为固定效应，个体影响作为随机效应对牛腿跖部温度和直肠温度数据进行混合线性模型分析。

[0014] 本发明还提供了所述牛体温的测量方法在牛机体生理情况检测中的应用。

[0015] 可选的，所述应用包括检测牛腿跖部温度通过数学模型计算得到牛直肠温度，利用所述牛直肠温度作为牛机体生理情况的指标。

[0016] 可选的，牛机体生理情况包括发情、排卵、妊娠、分娩。

[0017] 本发明的优点在于，本发明的测量部位选择为牛腿跖部，有利于温度检测设备的固定，而且可以降低因固定检测设备给动物造成的不适感；另一方面，本发明所提供的方法同以往的直肠手动测温法相比，不但减少了测定的人力物力投入，而且可实现24小时全天候检测，避免了因为人为作息时间制度下的温度检测死角(如夜间、非关注期)，为高产牛生理状态的适时预测和疾病的早期预警提供了技术基础，不但可望实现人力的解放，也为养殖业实现精细化管理提供了基础数据和硬件设施。本发明为深入研究牛发情、妊娠、分娩等生理状态和发热、中毒等伴随温度变化疾病的病理状态体温变化规律提供了技术支撑。本研究自主探索了尾根内侧、颈部、后腿蹄腕部等多个温度探头固定位置，最终通过数据分析，发现后腿跖部采集的温度数据较稳定，同时该位置固定牢固，避免了温度探头脱落现象。附图说明

[0018] 图1为体表测温装置安装示意图；

[0019] 图2为体表测温设备与水银温度计测温一致性分析的Bland—Altman图；

[0020] 图3为不同季节牛一天内体表温度和直肠温度变化曲线图。

[0021] 以下将对本发明的具体实施方式进行详细的说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。

[0022] 本发明提供了一种牛体温的测量方法，该方法包括：步骤一：首先检测牛腿跖部温度和直肠温度，再将检测时间和季节影响作为固定效应，个体影响作为随机效应对牛腿跖部温度和直肠温度数据进行混合线性模型分析，建立牛腿跖部温度和直肠温度关系的数学模型；步骤二：在需要测量牛体温时将测量获得的牛腿跖部温度带入所述数学模型计算获得牛直肠温度。在本发明所提供的方法中，所述牛腿跖部为牛后腿跖部，这是由于发明人在研究过程中发现，牛后腿跖部的温度的变化趋势相比于前腿更接近直肠温度的变化趋势，而且后腿固定设备可尽可能减少牛用嘴巴撕扯导致的设备损害，同时也避免牛只躺卧时采温设备因佩戴在前腿而使牛前腿受伤。

[0023] 在本发明的一种具体实施方式中，所述牛腿跖部为跖骨下部近籽骨上方，该部位由于跖骨下面相连的近籽骨凸起，可使测温设备得以固定。

[0024] 在本发明所提供的方法中，所述方法包括将体温检测装置固定于牛腿跖部，实时检测牛只温度数据，例如可以每小时自动检测2-10次，检测的时间间隔优选为相同的时间。

[0025] 所述数学模型是将检测时间和季节影响作为固定效应，个体影响作为随机效应对牛腿跖部温度和直肠温度数据进行混合线性模型分析。

[0026] 在本发明的一种实施方式中，所述混合线性模型如模型1所示：

[0027] 模型1：Y＝a+X0b0+Xsbs+Xtbt+Xibi+e；

[0028] 模型1中，Y为直肠温度，X0为牛腿跖部温度温度，Xs为季节设计矩阵，bs为Xs的非随机参数向量；Xt为时间设计矩阵，b0为X0的回归系数，bt为Xt的非随机参数向量；Xi为个体设计矩阵，bi为Xi的随机参数向量；a代表截距，e为随机向量误差。本发明所提供的方法对于测量体温所使用的工具没有特别的限制，只要能够准确测量跖部温度即可，例如可以为接触式测温或非接触式测温。优选的情况下，测量体温所使用的工具为能够固定于待测部位的，可以进行实时、连续、自动体温检测的工具。

[0029] 本发明还提供了本发明所述的方法在牛机体生理情况检测中的应用。

[0030] 可选的，所述应用包括检测牛腿跖部温度通过数学模型计算得到牛直肠温度，利用所述牛直肠温度作为牛机体生理情况的指标。

[0031] 可选的，牛机体生理情况包括发情、排卵、妊娠、分娩。

[0032] 具体的，可以将利用本发明所提供的方法计算得到的牛直肠温度应用于了解被测动物的生理情况，从而为指导畜牧生产提供依据。值得注意的是，其直接目的不是为了获得动物疾病的诊断结果或健康状况。

[0033] 下面将通过具体实施方式对本发明进行详细说明。需要理解的是以下实施例的给出仅是为了起到说明的目的，并不是用于对本发明的范围进行限制。本领域的技术人员在不背离本发明的宗旨和精神的情况下，可以对本发明进行各种修改和替换。

[0034] 对比例1

[0035] 选择19头健康牛只，使用手提式红外测温仪(DT8806H，测量精度为±0.1℃)分别测量眼窝、鼻镜、耳蜗、腹部无毛区、乳头、牛角基部和尾根内侧共7个部位的体表温度，同时用水银温度计测量肛门温度，将体表温度与肛门温度进行比较，结果用平均值±标准差表示(如表1所示)：

[0036] 表1

[0037]部位 测量方法 温度
眼窝 红外测温 36.27±0.54
鼻镜 红外测温 35.32±0.49
耳蜗 红外测温 37.78±0.45
腹部无毛 红外测温 35.37±0.49
乳头 红外测温 35.87±0.58
牛角基部 红外测温 35.40±0.47
尾根内侧 红外测温 37.52±0.46
肛门 水银温度计 38.50±0.28

[0038] 结果显示，不同部位温度与肛门温度均存在一定差距，其中只有耳蜗和尾根内侧温度与直肠温度最相近，且在7个测量部位中标准差最小，但耳蜗神经末梢较为丰富，固定温度探头牛只反应比较剧烈。同时考虑到温度探头的固定方法和可能性，分别研究了尾根内侧固定、颈部固定的测温效果。

[0039] 尾根内侧固定试验结果显示，虽然尾根内侧可固定接触式测温探头，但该位点距肛门较近，探头固定不牢，排粪时肛门抽动和牛尾摆动都不利于测温探头固定，且易受粪便影响，与其相连的电源线也不易固定，易导致温度探头丢失，无法正常采集数据，经过1个多月的尝试后不得不改换其他安装位置。

[0040] 颈部固定结果表明：颈部皮肤毛发较密且长，较厚实的毛发有较强的隔热效果，同时，颈部表皮直接暴露于体表，较易受外界环境影响，也不易固定。经试验，颈部固定法可正常采集体表温度数据，但是采集数据与体温偏离较大。表明颈部体表温度测定很难避免环境影响，温度测定结果对真实体温的指示程度较差，颈部测定途径不能达到理想效果。

[0041] 实施例1

[0042] 1材料方法

[0043] 1.1试验时间与地点

[0044] 分别于2014年冬(2月2日至10日；winter，Win.)、夏(6月2日至10日；summer，Sum.)、秋(10月2日至10日；autumn，Aut.)三个季节，在河北石家庄某牛场进行体表温度与直肠温度测定，因为春季与秋季气候条件相似，所以，仅选择秋季进行温度测定分析，以代表春秋季测定效应。

[0045] 1.2试验动物及其饲养管理

[0046] 每个季节分别随机选取7头健康、产后60至90天二胎正常可繁西门塔尔牛作为实验牛；试验牛采用散栏式饲养、全自动TMR日粮饲喂，每天分别在2:00、6:00、12:00、18:00、22:00投放饲料、保证自由采食和清洁饮水。

[0047] 1.3测温设备及其测试与校正

[0048] 市售水银温度计1支，兽用水银体温计2支，体表温度自动测定系统，包括体表温度检测装置7套(该装置采用接触式测温)、数据采集装置1套(该装置用于采集检测装置中存储的体温数据)。如图1所示，该系统采用接触式测温方式，选择牛体表被毛稀疏的后肢蹄腕部为体表测温部位，将蹄腕部肌肉、血管较丰厚部位(即牛腿跖部，牛腿跖骨下部近籽骨上方)作为温度传感器固定位置，并根据该部位形状设计了具有一定弧度的固定外壳，通过外壳与腿部紧密贴合使传感器能够与肌体紧密接触。同时，又在外壳中与肌体贴合一侧的中心设计了一个凹巢，温度传感器固定于凹巢内，进一步降低了温度、风力等外界环境因素对测温效果的影响，检测获得的温度。

[0049] 为使实验结果准确科学，每次试验前将7个体表温度检测装置置于恒温箱中，进行测试实验。首先设定恒温箱温度从30℃升至40℃(每0.5h升温1℃)，测温设备每0.5h自动测温一次。试验期间，使用胶带将水银温度计固定在恒温箱玻璃门内侧，在自动测温时间点，同时读取并记录水银温度计显示的恒温箱内温度。该测试实验重复三次，证明体表温度自动测定系统准确可靠后，进行后续体表温度自动测定实验。

[0050] 1.4体表温度和体内温度的采集

[0051] 预试验：正式试验开始前两天，将体表温度(surface temperature，ST)检测装置紧密固定于实验牛后腿蹄腕部(即牛后腿跖部)，并将测温间隔设定为半小时，每4个小时用数据采集装置收集一次温度数据；根据温度采集结果进行设备调试，同时进行试验牛适应性训练。待设备测温正常稳定、试验牛适应后，进行正式试验。

[0052] 正式试验：自动测温装置继续沿用预试期的半小时测温设置，连续测温3天；同时，采用水银体温计，每2小时测定一次试验牛直肠温度(rectal temperature，RT)。测定直肠温度时，先将试验牛保定，待其休息10分钟后将温度计插入直肠，5分钟后取出温度计，读取数据，每头牛每天共采集12个时间点的直肠温度数据，也连续测温3天。

[0053] 1.5数据处理与分析

[0054] 设备实验室测试数据，使用SAS9.2进行差异显著性分析与线性回归分析，并使用Bland—Altman法进行一致性评价。

[0055] 每头牛每天有48条体表温度数据与12条直肠温度数据。因直肠温度两个小时内变化较小，故取相邻两个时间点的直肠温度平均值作为中间时间点的直肠温度，并舍弃所有非整点体表温度数据，使每头牛每天各有23个相同时间点的体表温度与直肠温度数据，三个季节一共获得1449条相同时间点的体表温度与直肠温度数据。采用SAS 9.2进行体温数据差异显著性分析与协方差分析，同时采用Tablecurve2D/3D进行线性回归分析与非线性回归分析，并进行体表温度校正分析。

[0056] 2结果与分析

[0057] 2.1体表温度检测装置的测试分析

[0058] 实验室测试结果表明，恒温箱温度在30～40℃范围内变化时，7个检测装置的测定2
结果发生了同步变化，且差异不显著，自动测温结果与恒温箱内实际温度呈强相关(R ：
0.9987)。测温结果一致性金标准分析得到的Bland—Altman图显示，7个测温设备间的差异仅为0±0.07℃(95％置信区间为-0.22～0.22℃，图2)，可以忽略不计。说明该装置的准确性和稳定性较好，可用于后续体表温度研究。

[0059] 2.2体表温度与直肠温度的季节及测定日差异分析

[0060] 跖部体表温度与直肠温度统计分析结果(表1、2)显示，直肠温度变化幅度较小，各季节间、各测定日间均无显著差异(p＞0.05)，说明生理条件下，直肠温度基本恒定。各测定日间体表温度(表2)差异也不显著(p＞0.05)，说明自动测温装置性能较为稳定。然而，不同于直肠温度，跖部体表温度测定结果季节差异显著(p＜0.05)，其中，夏季温度最高，冬季温度最低，各季节的测定结果均显著低于直肠温度测定结果(p＜0.05)。

[0061] 表1体表温度与直肠温度统计分析结果

[0062]                                               单位：℃[0063]

[0064] 注：A、B、C表示同行数据的差异显著性，字母相同者表示差异不显著(p＞0.05)，字母不同者表示差异显著(p＜0.05)；

[0065] 表2相同季节三天体表温度与直肠温度情况

[0066]

[0067] 注：A、B、C表示同列数据的差异显著性，字母相同者表示差异不显著(p＞0.05)，字母不同者表示差异显著(p＜0.05)；

[0068] 2.3跖部体表温度和直肠温度的昼夜差异分析

[0069] 表3体表温度与直肠温度昼夜统计分析

[0070]

[0071] 注：A、B表示同列数据的差异显著性，字母相同者表示差异不显著(p＞0.05)，字母不同者表示差异显著(p＜0.05)；

[0072] 实验牛不同季节各测定时间点温度曲线(图3)与统计分析结果(表3)表明，体表温度受时间因素影响显著并呈现明显昼夜变化规律，5:00至7:00为一天中的最低值，随后逐渐升高，到13:00到15:00左右达到一天的最高值，随后逐渐降低，直至最低温度。直肠温度变化幅度较小，但变化规律类似。

[0073] 2.4跖部体表温度与直肠温度回归分析

[0074] 假设体表温度与直肠温度间存在线性关系，将体表温度对直肠温度进行回归分析。

[0075] 2.4.1仅考虑季节为固定效应

[0076] 以季节因素为固定效应，个体为随机效应，对体表温度和直肠温度数据进行混合线性模型分析。建立的混合线性模型A如下：

[0077] Y＝a+X0b0+Xsbs+Xibi+e其中，Y为直肠温度，X0为牛腿跖部温度，Xs为季节设计矩阵，bs为Xs的非随机参数向量，即固定效应；Xi为个体设计矩阵，b0为X0的回归系数，bi为Xi的随机参数向量，即随机效应；a代表截距，e为随机向量误差。

[0078] 采用模型A对体表温度和直肠温度进行混合线性分析，得到模型的参数估计见表4，模型解释的变异为0.399。

[0079] 表4仅考虑季节效应因素为固定效应的混合线性模型结果

[0080]

[0081] 2.4.2同时考虑季节和测量时间为固定效应

[0082] 不考虑时间因素影响得到的决定系数仅为0.399，效果不够理想。考虑到单因素方差分析提示不同测量时间对体表测温结果有显著影响。因此，在前面模型的基础上，加入测量时间作为固定效应，建立的混合线性模型B如下：

[0083] Y＝a+X0b0+Xsbs+Xtbt+Xibi+e其中，在模型A的基础上，Xt为时间设计矩阵，bt为Xt的非随机参数向量。模型B中各参数的估计值见表5，模型解释的变异为0.562。因此，考虑时间因素的回归模型优于不考虑时间因素的回归模型。

[0084] 表5同时考虑季节因素测定时间因素为固定效应的混合线性模型结果[0085]

[0086]

[0087] 2.5直肠温度与校正后体表温度对比

[0088] 运用混合线性模型B得出的回归方程对体表测温数据进行矫正，各校正温度与直肠温度差值绝对数的平均数为0.104℃，90％数据差值的绝对值不超过0.2℃，说明该拟合方程与本测温系统结合，可以取得较理想的测温结果。

[0089] 2.6交叉验证

[0090] 为验证上述线性模型的预测效果，将全部数据随机平均分为10等份，随机选取9份建立混合线性模型，剩下1份使用建立的模型进行预测，依次重复10次，得到每个体表温度对应的预测直肠温度。最后对预测的直肠温度与真实直肠测温数据做相关分析以评价模型的预测效能。考虑时间点因素模型的相关系数为0.644，明显高于不考虑时间点因素模型的相关系数0.569。

[0091] 通过SAS9.3差异显著性分析后发现，校正后的体表温度与真实直肠温度间差异不显著(p＞0.05)，说明该方法校正后得到了理想的效果，同时也说明牛体表温度与直肠温度之间存在良好的曲线关系。可以通过将测量得到的牛腿跖部温度带入混合线性模型B，计算获得牛直肠温度，完成牛体温的测量。即可以通过本发明所提供的方法，利用体表温度直接获得直肠温度数据。

[0092] 虽然，上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。