CN101112306A - 无创体核温度测量的方法、装置和标定设备及其标定方法 - Google Patents
无创体核温度测量的方法、装置和标定设备及其标定方法 Download PDFInfo
- Publication number
- CN101112306A CN101112306A CNA2007101067969A CN200710106796A CN101112306A CN 101112306 A CN101112306 A CN 101112306A CN A2007101067969 A CNA2007101067969 A CN A2007101067969A CN 200710106796 A CN200710106796 A CN 200710106796A CN 101112306 A CN101112306 A CN 101112306A
- Authority
- CN
- China
- Prior art keywords
- temperature
- value
- core
- body surface
- core temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
无创体核温度测量的方法、装置和标定设备及其标定方法,本发明在人体体表设置绝热材料层,绝热材料层内侧与人体体表之间设置体表温度传感器,用于检测人体体表温度(Ts),在绝热材料层外侧设置环境温度传感器,用于检测环境温度(Te),并设置用于产生模拟环境温度的电热部件;通过测量得到的体表和环境温度数据,用模拟电路进行传热过程的分析,用增加瞬态项模拟电路分析法估算人体体核温度(Tc)。采用本发明方法和装置可对人体进行长期、连续、动态、个性化的无创体核温度测量,使用方便,制造和运行成本低,准确性高。采用本发明标定设备和其使用方法可提高用本发明方法和装置测量人体体核温度的准确性。
Description
技术领域
本发明属一种获取人体动态生理参数的方法和装置,具体是一种无创体核温度的测量方法、装置和该装置的标定设备及其标定方法。
背景技术
人体体核温度的动态变异性包含着体温调节***极为丰富的信息,是重要的生命体征,它可以作为一种中间参数为人体多种疾病的综合诊断提供参考依据。体核温度的连续、动态测量分有创和无创两类方法,有创法给人体带来一定的损伤,对于需要长期连续监测的场合,有创法有很大的不便。无创体温测量方法有热绝缘法、零热流法等。热绝缘法假设体表做绝热处理后,皮肤温度等于体核温度。热绝缘法测试时间长,精度因绝热材料性能的不同而受到很大的限制。零热流法的基本思路是在体表给人体发送热流以补偿体表散热,当热流为零时认为皮肤温度等于体核温度。零热流法需要加热、冷却装置来平衡热流,装置较为复杂。另外还有学者进行了诸如红外法,微波法,超声法,NMR核磁共振法测量体核温度的研究。这些方法都有各自特定的应用条件。
发明内容
本发明要解决的技术问题是提供一种无创体核温度测量的方法、装置和该装置的标定设备及其标定方法,用本发明方法和装置可实现对人体进行长期、连续、动态、个性化的无创体核温度测量,操作简便易行、精确度高,采用本发明标定设备及其标定方法可提高用本发明方法和装置测量人体体核温度的准确性。
本发明体核温度测量的方法包含下述内容:
1、在人体体表设置绝热材料层,绝热材料层内侧与人体体表之间设置体表温度传感器,用于检测人体体表温度(Ts),在绝热材料层外侧设置环境温度传感器,用于检测环境温度(Te),并设置用于产生模拟环境温度的电热部件;
2、用模拟电路进行传热过程的分析和计算,包含下述过程:
用电热部件模拟环境升温(ΔTe),测得此时皮肤温度变化量(ΔTs),以此获得人体体表组织等效热阻(Rc)与所述绝热材料层等效热阻(Re)的比值作为个性化校正系数(A),并获得环境温度变化量(ΔTe)与皮肤温度变化量(ΔTs)的比值作为环境温度对皮肤温度影响的环境温度校正系数(α);
用电子体温计测得人体某一时刻的体核温度值(Tc0),同时测得此刻的体表温度值(Ts0)和环境温度值(Te0),作为参考基准点;
完成上述过程后即开始测量估算:测量人体体表温度值(Ts)和环境温度(Te),根据所述绝热材料层的等效热阻(Re)、热容量(Cs),和前述获取的个体化校正系数(A)和环境温度正系数(α),即可用瞬态模拟电路分析法估算人体体核温度(Tc),在估算中,用个体化校正系数(A)校正人体体表组织等效热阻(Rc),用环境温度校正系数(α)去除环境温度(Te)对皮肤温度(Ts)的影响,估算中以所述参考基准点的体核温度值(Tc0)作校准。
本发明无创体核温度测量的装置为:设有可固定在人体上的绑缚带(1),绑缚带上设置绝热材料层(2),绝热材料层内部贴着人体体表的部位设有体表温度传感器(3),绝热材料层外部设有环境温度传感器(4)和用于产生模拟环境温度的电热部件(5),所述的体表温度传感器和环境温度传感器输出端与信号采集显示装置连接。
下面对本发明方法的原理作一说明,本发明方法采用电路模拟传热过程的分析和计算:参见
图2本发明方法用电路模拟传热过程的稳态导热简化热阻分析图
图3本发明方法所采用的稳态导热时的简化电路模拟图
图4本发明方法所采用的增加瞬态瞬态项导热变化时的简化电路模拟图。
根据上述电路模拟图可分析出:皮肤温度传感器测量人体体表温度值Ts受四种因素影响,即体核温度Tc,环境温度Te,绝热材料物理性质(热阻Re,热容量Cs)和体表组织等效热阻Rc。Tc与上述四因素的简化数学关系,由瞬态导热电路模拟图可推导为(根据基尔霍夫电流定律):
Tc=Rc×Cs×dTs/dt+Rc/Re×(Ts-Te)+Ts (1)
Tc与上述四因素在稳态导热阶段的简化数学关系为:
Tc=Rc/Re×(Ts-Te)+Ts (2)
上述公式(1)、(2)即是本发明方法采用瞬态模拟电路分析法估算人体体核温度(Tc)的基本表达式,从上述式(1)可看出:人体体表温度Ts和环境温度Te可通过传感器测得,绝热材料层热阻Re、热容量Cs(主要取决于密度ρ和比热容c)可以人为设计,只有人体体表组织等效热阻Rc不易获取,而且Rc还会因不同人体的体表脂肪、肌肉等因素的个体差异而有较大不同,另外,人体体表温度Ts的变化是人体体核温度Tc和环境温度Te变化共同作用的结果,如果环境温度Te恒定,则Ts的变化就能直接反映体核温度Tc的动态变化,因此,需要从Ts的变化中去除环境温度变化的影响。
因此本发明方法在测量估算前,先模拟环境温度变化ΔTe,测得稳定后的皮肤温度变化量ΔTe,根据式(2)即可求得
Rc/Re=ΔTs/(ΔTe-ΔTs)=A
将A作为体现人体体表组织等效热阻Rc的个性化校正系数;
同时将可求得的ΔTs/ΔTe=α作为体现从皮肤温度变化中去除环境温度影响的环境温度校正系数。
这样即可实现用增加瞬态项模拟电路分析法估算人体体核温度(Tc)了。
本发明测量人体体核温度装置的标定设备结构为:设装有液体介质(12)的容器(9),容器内设有电加热器(13),容器底部设有水泵搅拌器(11),液体内设有体核温度传感器(8),其输出端引出容器外,与数据采集显示装置连接。
本发明标定设备的标定方法是:用于标定由本发明装置无创体核温度测量装置绝热材料层热阻Re和热容量Cs决定的瞬态项系数γ标=Re×Cs,和用于校准用本发明方法估算的体核温度变化值的误差补偿系数β标,具体过程为:
1、维持标定设备容器内核心温度不变,将本发明无创体核温度测量装置绑缚在本发明标定设备容器的外壁上,稳定后测出模拟体表温度(Ts);
2、维持标定设备容器内核心温度不变,用微电热片模拟环境升温ΔTeh,检测模拟体表温度变化为ΔTsh;
3、微电热片停止加热,稳定后,将某点作为基准点,记录该点模拟体表温度值、环境温度值和容器内核心温度值,分别记为(Tso、Teo、Tco);
4、以基准点为起始点,使标定设备容器内的温度模拟人体体核温度变化,并连续记录和存储模拟体表温度(Ts)、环境温度(Te)和核心温度值(Tc);
5、多次重复上述1-4步骤,并将每次得到的一组测量数据保存待用;
然后包含下述内容:
6、对前述每组测量数据分别进行处理,分别取得每组数据的
个体化校正系数A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)和
环境校正系数α=ΔTsh/ΔThe ;
7、根据上述的每组数据,用本发明方法估算容器内的模拟体核温度Tc,将估算的Tc值与实测的标定设备容器核心温度值做对比,观察二者之间的误差,同时人为调节瞬态项系数中的Re×Cs值,取使瞬态误差最小时的Re×Cs值作为该组测量数据的瞬态系数γ值;
求各组瞬态系数γ值的平均值γ标,将其作为用本发明装置和方法估算人体体核温度时的瞬态项系数Re×Cs的标定值;
8、根据上述的每组数据,求出各测量数据Tc、Ts、Te相对于其基准点Tco、Tso、Teo的变化量;
计算环境温度变化造成的体表温度变化量Tse:
Tse=α×(Te-Teo)
计算由单纯核心温度变化造成的体表温度变化量:
dTsc=(Ts-Tso)-α×(Te-Teo);
用本发明方法估算核心温度变化量dTc,将估算的核心温度变化量dTc与实测的对应核心温度变化量相比,得到由体表温度变化量估算体核温度变化量的误差补偿系数β:
β=体核温度变化估算值/体核温度变化实测值=dTc/(Tc实-Tco);
求各组误差补偿系数β的平均值作为用本发明方法估算人体体核温度时的误差补偿系数β标。
本发明提供了一种无创人体体核温度测量的新方法和新装置,采用本发明方法和装置可对人体进行长期、连续、动态、个性化的无创体核温度测量,使用方便,制造和运行成本低,准确性高。采用本发明标定设备和其标定方法可提高用本发明方法和装置测量人体体核温度变化的准确性。
附图说明
图1、热绝缘法测量体核温度方框示意图
图2、本发明方法用电路模拟传热过程的稳态导热简化热阻分析图
图3、本发明方法所采用的稳态导热时的简化电路模拟图
图4、本发明方法所采用的增加瞬态导热变化时的简化电路模拟图。
图5、本发明装置方框结构示意图
图6、实施例1本发明装置结构示意图
图7、实施例3本发明标定设备结构示意图
图8、实施例5瞬态系数γ标标定示意图
图9、实施例5误差补偿系数β标定示意图
1-绑缚带 2-绝热材料层 3-体表温度传感器 4-环境温度传感器5-电热元件 6-弹性带 7-锁扣装置 8-体核温度传感器 9-容器10-盖 11-搅拌泵 12-介质 13-电加热元件 14-测量盘
具体实施方案:
实施例1
参见图6,为本发明体核温度测量装置的一个实施例,该装置设有带绑缚带1的测量盘14,绑缚带由软布制成,其上有弹性材料6和锁扣装置7以保证对人体的缚紧度。测量盘14上设有绝热材料层2,其内侧贴着人体皮肤的部位装有体表温度传感器3,靠近外侧的绝热材料层内嵌有环境温度传感器4和与电池连接的电热元件5,测量盘呈长方体,尺寸为40*40*15mm3(长*宽*高),重10g,电热元件5为35*35*1mm的微电热软薄片,镶嵌于绝热材料中。微电热片与体表温度传感器3之间被绝热材料层2隔离,二者垂直距离为10mm,体表温度传感器3可经过良导热体与皮肤表面紧密接触,环境温度传感器4与微电热片5良导热接触,体表温度传感器3、环境温度传感器4、微电热片电热元件5与测量盘中心重合,体表温度传感器3和环境温度传感器4输出端分别与温度信号采集装置连接。
实施例2
本例为本发明体核温度测量方法的一个实施例,将实施例1所述装置固定在人体上,即可进行体核温度估算,具体过程如下:
1.将测量装置绑缚到测量部位,稳定30分钟;
2、用微电热片电热元件5模拟环境升温ΔTeh,稳定20分钟后,检测皮肤温度变化为ΔTsh,则可得:
个体化校正系数:A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)
环境温度校正系数:α=ΔTsh/ΔTeh
3.停止加热后,记录某一时刻的皮肤温度Ts和环境温度Te,同时用电子体温计测量体核温度(直肠或腋窝)并记录下来,分别记为Tso、Teo、Tco,作为参考基准点;
4、开始进行体核温度的估算:
采集不同时刻人体的体表温度Ts和环境温度Te,采用公式3进行体核度估算:
Tc=A×Re×Cs×[(dTs-α×dTe)/dt]+
A×[(Tso+dTs-α×dTe)-(Teo+dTe)]+[(Tso+dTs)-α×dTe]
(3)
式中:Tc为人体体核温度估算值
Cs为绝热材料层的热容量
Ts为即时测得的人体体表温度值
Te为即时测得的环境温度值
Re为绝热材料层的热阻
A为已通过环境温升ΔTeh获取的个性化校正系数
A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)
α为通过环境温升ΔTeh获取的环境温度校正系数α=ΔTsh/ΔTeh
然后用参考基准点的体核温度值Tco对用公式3估算的值做单点校准,即将上述按公式3估算的人体体核温度Tc初始点值和基准点所测的Tco值做差值比较,得到的差值作为后续体核温度估算值的校正量。
在上述估算公式中,热绝缘材料热阻Re、热容量Cs可人为设计(例如通过改变绝热材料层的密度ρ和比热容c),A和α也已求出,(dTs-α×dTe)为去除了环境温度影响后的皮肤温度变化量,所以可以实现根据测得的Ts和Te估算体核温度Tc。
实施例3
本例为本发明体核温度测量方法的又一个实施例
本方法与实施例2的不同点在于第4步体核温度的估算方法有所不同,本例在求得个体化校正系数A,环境温度校正系数α,和参考基准点的Tso、Teo、Tco后,按下面的公式进行体核温度Tc的估算:
稳态量估算:
Tc稳=Tco+(dTs-α×dTe)/(1-α) (4)
瞬态量估算:
Tc瞬=A×Re×Cs×(dTs-α×dTe)/dt (5)
则体核温度估算值为:Tc=Tc稳+Tc瞬
式中
Tc为人体体核温度估算值
Cs为绝热材料层的热容量
Ts为即时测得的人体体表温度值
Te为即时测得的环境温度值
Re为绝热材料层的热阻
A为通过环境温升ΔTeh获取的个性化校正系数
A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)
α为通过环境温升ΔTeh获取的环境温度校正系数
α=ΔTsh/ΔTeh
然后用参考基准点的体核温度值Tco对用公式3估算的值做单点校准,即将上述估算的人体体核温度Tc初始点值和基准点所测的Tco值做差值比较,得到的差值作为后续体核温度估算值的校正量。
上述方法的估算原理是:
如前述:个体化校正系数A=Rc/Re=dTsh/(dTeh-dTsh)
环境温度校正系数α=dTsh/dTeh。
所以:A与α的关系为:A=α/(1-α)
环境温度不变时,公式1
Tc=Rc×Cs×dTs/dt+Rc/Re×(Ts-Te)+Ts
=A×Re×Cs×[(dTs-α×dTe)/dt]+A×(Ts-Te)+Ts
中稳态量中体核温度的变化成分为:
dTc=A×dTs+dTs=α/(1-α)×dTs+dTs=dTs/(1-α)
去除环境温度对皮肤温度的影响,即得:
dTc=(dTs-α×dTe)/(1-α) (6)
故可得公式5稳态量
Tc稳=Tco+(dTs-α×dTe)/(1-α) (4)
而根据公式(3),瞬态量估算:
Tc瞬=A×Re×Cs×[(dTs-α×dTe)/dt] (5)
其中Re×Cs可根据事先对测量装置绝热材料层标定的瞬态系数
γ标=Re×Cs值进行计算
从而可以得到体核温度估算值为稳态值与瞬态值之和:
Tc=Tc稳+Tc瞬
为了提高上述估算值的准确性,进一步方案是用误差补偿系数β标对其进行补偿,进行补偿的估算方法是:
用经过标定的误差补偿系数β对前述的稳态量进行补偿后得
(Tc稳)补=Tco+(dTs-α×dTe)/[β标×(1-α)]
瞬态量的估算与前述相同;
则体核温度估算值为:
Tc=((Tc稳)补+Tc瞬。
实施例4
本例是本发明无创体核温度测量装置标定设备的实施例,本标定设备是一个体温模拟发生器,具体结构参见图7,容器模拟人体躯干呈圆桶形,总体尺寸定为桶高250mm,内径80mm,ABS底座25mm厚。顶盖10为10mm厚ABS材料,顶盖上留有走线孔和排气孔,并设有支撑电热元件13的绝缘装置。电加热元件为交流220V,500W型圆柱钢管加热棒,长200mm,直径10mm,用来模拟体核热源。桶内介质为水,12V4W的潜水泵固定在底座上,搅动水流,使桶内水温均匀。
体温模拟发生器的桶壁用来模拟体核到体表的热阻,其厚度确定如下。人体皮肤导热系数=0.47W/(m.k),脂肪导热系数=0.21W/(m.k),正常腹部男性皮肤皱襞(皮肤+皮下脂肪)厚度为5-15mm,女性皮肤皱襞为12-20mm。本例取平均皮厚=10mm,其中表皮厚=2mm,脂肪厚=8mm,根据导热稳态阶段,总热阻等于串联热阻之和,即:
代入数据计算得皮皱平均导热系数为λa=0.236W/(m.k)。ABS材料导热系数λabs=0.2512W/(m.k),则同样根据上述热阻相等的原理,等价10mm皮皱厚度的ABS材料厚度为
所以模型中ABS桶壁厚度定为11mm。
为控制体温模拟发生器内水温,另设计了温度控制板,控温精度在±0.1℃,以便根据实验要求来实现体温模拟发生器内温度的恒定和变化。
桶内设有模拟体核温度传感器8,用于检测模拟的体核温度。
实施例5
本例为本发明无创体核温度测量装置标定设备标定方法的实施例。
本例采用实施例4所述的标定设备和实施例3所述的测量盘,其温度传感器输出端与温度信号采集控制板连接,用于采集标定设备和测量盘的温度信号,并控制标定设备容器内核的温度。
为了正确估算体核温度瞬态量和减小估算中存在的误差,新测量盘在使用前要用标定设备对其参数进行标定,本方法即可用于标定本发明无创体核温度测量装置,同时可用于校准用本发明无创体核温度测量装置和测量方法进行侧量得到的估算值。
本方法的目的一是标定出实施例1所述测量盘的瞬态系数γ=Re×Cs,以便使其可用实施例2或实施例3所述方法进行人体体核温度瞬态成分的估算;二是求出用本发明进行人体体核温度变化估算产生的误差补偿系数β,以补偿估算存在的误差。
具体过程是:
1、室温25℃,模拟体温发生器核心温度控温在37℃。测量盘按图7所示绑缚在模拟体温发生器上,稳定30分钟。此时,在环境温度波动不大的情况下(dTe<1℃),体表温度接近一个稳定值;
2、模拟体温发生器核心温度仍控温在37℃不变,微电热片模拟环境升温ΔTeh,检测皮肤温度变化为ΔTsh;
3、微电热片停止加热,稳定30分钟,使微电热片温度回落到升温前的温度点。然后,同时记录一点体表温度值、环境温度值和核心温度值,分别记为Tso、Teo、Tco;
4、模拟体温发生器核心温度在36-41℃范围内以1℃的梯度变动,每个温度梯度上保持1小时。同时连续记录体表、环境和核心温度值,分别以Tst、Tet和Tct记录在计算机中。
5、重复上述1-4,测量10次,将每次测量得到的一组结果保存在不同的目录中待用。
6、根据上述各组数据标定瞬态系数γ:
公式(1)中体核温度中瞬态项Rc×Cs×dTs/dt中的系数不宜直接取
得,但通过微电热片模拟环境温度变化Teh,可得到:
A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)
所以Rc=A×Re=ΔTsh/(ΔTeh-ΔTsh)×Re
则瞬态项表示为:
A×Re×Cs×dTs/dt=ΔTsh/(ΔTeh-ΔTsh)×Re×Cs×dTs/dt
因为A=ΔTsh/(ΔTeh-ΔTsh)部分可由前述第2步微电热片模拟环境温度变化ΔTsh/求得,而Re×Cs部分是一个与测量盘本身热物理性质有关的量,一旦测量盘的物理性质和形状确定了,Re×Cs就确定了。所以Re×Cs也是一个定量,我们将Re×Cs作为瞬态系数γ,具体标定过程是:
对10组实验数据,做如下的处理:
根据微电热片模拟环境温度变化数据,取得个体化校正系数:
A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)
取得环境校正系数
α=dTsh/dTeh
用公式3:
Tc=A×Re×Cs×[(dTs-α×dTe)/dt]+
A×[(Tso+dTs-α×dTe)-(Teo+dTe)]+[(Tso+dTs)-α×dTe]
(3)
估算模拟体温发生器核心温度,并和实测的模拟体温发生器核心温度做对比,同时人为调节Re×Cs值,观察核心温度估算值和实测值间瞬态误差大小,如图8所示,取瞬态误差最小时的Re×Cs值为该组实验的瞬态系数γ值(图中Re×Cs=80时得到的Tc1最接近实测值)。
这样10组标定实验数据就得到10个瞬态系数γ值,如表1所示。则该测量盘的瞬态系数为10次γ取平均值,γ’=87.1,该值即可作为用实施例1测量盘测量人体体核温度时所用公式(3)进行估算时瞬态项的标定系数.
标定实验瞬态系数表
实验 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 均值 |
γ=A×Re×Cs | 90 | 85 | 88 | 91 | 86 | 88 | 90 | 84 | 80 | 89 | 87.1 |
7、根据上述各组数据标定误差补偿系数β
测量中发现,在由体表温度变化量dTs推导体核温度变化量dTc的过程中,由于基准点的选取、环境温度影响扣除等操作存在一定的误差,所以导致核心温度变化量dTc的估算值和实际值间存在一定的误差,为对这种误差给以补偿,本例对用实施例1所述测量盘,实施例2、3所述估算体核温度的测量方法标定一个误差补偿系数β,方法如下:
对前述1-5得到的10组实验数据,按如下步骤做处理:
①根据记录的微电热片模拟环境升温ΔTeh值和相应体表温度变化值ΔTsh,可以算得:
个体化校正系数::A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)
环境温度校正系数:α=ΔTsh/ΔTeh
②取实验步骤3记录的数据值Tco,Tso和Teo分别作为核心温度,体表温度和环境温度变化的基准(参见图9);
③计算实验步骤4实验数据相对于基准的变化量(参见图9)
核心温度Tc相对其基准的变化量:Tc-Tco;
体表温度Ts相对其基准的变化量:Ts-Tso;
环境温度Te相对其基准的变化量:Te-Teo。
④计算环境温度变化造成的体表温度变化量:α×(Te-Teo);
⑤由于体表温度的变化是核心温度和环境温度共同作用的结果,即:体表温度变化量=核心温度变化作用量+环境温度变化作用量,所以这里计算出核心温度变化造成的体表温度变化量:
ΔTsc=(Ts-Tso)-α×(Te-Teo);
⑥由单纯核心温度变化造成的体表温度变化量dTsc来估算核心温度变化量的稳态成分:这里采用实施例3所述的估算式6进行估算:
dTc=(dTs-α×dTe)/(1-α) (6)
⑦将⑥估算的核心温度变化量dTc与前述第四步中得到的实测稳态核心温度变化量相比,就得到由体表温度变化量估算体核温度变化量的误差补偿系数β,即
β=体核温度变化估算量/体核温度变化实测量
=dTc/(Tc实-Tco)
⑧、10组标定实验数据就得到10个β值,如表2所示,10个β取其平均值得0.8835。将该值作为用体表温度变化量估算体核温度变化量的误差补偿系数β标。
将用本发明方法估算出的体核温度变化量用误差补偿系数β标修正后,再与实测核心温度变化量做对比,其误差大为降低(参见表2)。
表2用误差补偿系数估算核心温度变化量时的误差表
实验 | β标 | 用β标补偿后的Tc变化量(℃) | 实测Tc变化量(℃) | 用β标均值补偿后Tc变化量估算误差(℃) |
1 | 0.9412 | 2.762 | 2.837 | -0.075 |
2 | 0.8451 | 1.036 | 0.988 | 0.048 |
3 | 0.7852 | 1.326 | 1.213 | 0.113 |
4 | 0.8813 | 1.259 | 1.253 | 0.006 |
5 | 0.938 | 2.758 | 2.827 | -0.069 |
6 | 0.895 | 1.12 | 1.074 | 0.046 |
7 | 0.894 | 0.823 | 0.849 | -0.026 |
8 | 0.923 | 2.175 | 2.235 | -0.06 |
9 | 0.943 | 1.985 | 2.151 | -0.166 |
10 | 0.789 | 1.352 | 1.216 | 0.136 |
均值 | 0.8835 | -0.0047 |
通过以上瞬态系数γ标和误差补偿系数β标的标定,标定过程就结束了,用标定后的测量盘进行实际的测量估算并用误差补偿系数对估算值进行补偿可提高估算结果的准确性。
Claims (7)
1.无创体核温度测量的方法,其特征在于,包含下述内容:
(1)、在人体体表设置绝热材料层,绝热材料层内侧与人体体表之间设置体表温度传感器,用于检测人体体表温度(Ts),在绝热材料层外侧设置环境温度传感器,用于检测环境温度(Te),并设置用于产生模拟环境温度的电热部件;
(2)、用模拟电路进行传热过程的分析和计算,包含下述过程:
用电热部件模拟环境升温(ΔTeh),测得此时皮肤温度变化量(ΔTeh),以此获得人体体表组织等效热阻(Rc)与所述绝热材料层等效热阻(Re)的比值作为个性化校正系数(A),并获得环境温度变化量(ΔTeh)与皮肤温度变化量(ΔTsh)的比值作为环境温度对皮肤温度影响的环境温度校正系数(α);
用电子体温计测得人体某一时刻的体核温度值(TC0),同时测得此刻的体表温度值(TS0)和环境温度值(Te0),作为参考基准点;
完成上述过程后即开始进行测量估算:测量人体体表温度值(Ts)和环境温度(Te),根据所述绝热材料层的等效热阻(Re)、热容量(Cs),和前述获取的个体化校正系数(A)和环境温度校正系数(α),即可用增加瞬态项模拟电路分析法估算人体体核温度(Tc),在估算中,用个体化校正系数(A)校正人体体表组织等效热阻(Rc),用环境温度校正系数(α)去除环境温度(Te)对皮肤温度(Ts)的影响,估算中以所述参考基准点的体核温度值(TC0)作校准。
2.根据权利要求1所述的无创体核温度测量的方法,其特征在于:所述用增加瞬态项模拟电路分析法估算人体体核温度(Tc)的方法是:用下述估算式(3)进行体核温度估算:
Tc=A×Re×Cs×[(dTs-α×dTe)/dt]+
A×[(Tso+dTs-α×dTe)-(Teo+dTe)]+[(Tso+dTs)-α×dTe]
(3)
式中:Tc为人体体核温度估算值
Cs为绝热材料层的热容量
Ts为即时测得的人体体表温度值
Te为即时测得的环境温度值
Re为绝热材料层的热阻
A为通过环境温升ΔTeh获取的个性化校正系数
A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)
Rc为人体体表组织等效热阻
α为通过环境温升ΔTeh获取的环境温度校正系数
α=ΔTsh /ΔTeh
Tso为参考基准点体表温度值
Teo为参考基准点环境温度值
然后将按上述公式估算的人体体核温度(Tc)初始点值和基准点(Tco)值做差值比较,得到的差值作为后续体核温度估算值的校正量。
3.根据权利要求1所述无创体核温度测量的方法,其特征在于:所述用增加瞬态项模拟电路分析法估算人体体核温度(Tc)的方法是:
按下述估算式进行体核温度Tc的估算:
稳态量估算:
Tc稳=Tco+(dTs-α×dTe)/(1-α) (4)
瞬态量估算:
Tc瞬=A×Re×Cs×(dTs-α×dTe)/dt (5)
则体核温度估算值为:Tc=Tc稳+Tc瞬
式中
Tc为人体体核温度估算值
Cs为绝热材料层的热容量
Ts为即时测得的人体体表温度值
Te为即时测得的环境温度值
Re为绝热材料层的热阻
A为通过环境温升ΔTeh获取的个性化校正系数
A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)
α为通过环境温升ΔTeh获取的环境温度校正系数
α=ΔTsh/ΔTeh
然后将按上述公式估算的人体体核温度(Tc)初始点值和基准点(Tco)值做差值比较,得到的差值作为后续体核温度估算值的校正量。
4.根据权利要求3所述的无创体核温度测量的方法,其特征在于:用误差补偿系数β标对估算的体核温度估算值的稳态量进行补偿,用下式进行稳态量估算:
Tc稳=Tco+(dTs-α×dTe)/β标(1-α)所述的误差补偿系数β标为根据标定测量得到的测量估算值与实际值进行比较得出的误差所标定的误差补偿系数,所述估算式中瞬态项系数(Re×Cs)采用标定测量中的标定值(Re×Cs)=γ标。
5.无创体核温度测量的装置,其特征在于,设有可固定在人体上的绑缚带(1),绑缚带上设置绝热材料层(2),绝热材料层内部贴着人体体表的部位设有体表温度传感器(3),绝热材料层外部设有环境温度传感器(4)和用于产生模拟环境温度的电热部件(5),所述的体表温度传感器和环境温度传感器输出端与信号采集装置连接。
6.无创体核温度测量装置的标定设备,其特征在于:设装有液体介质(12)的容器(9),容器内设有电加热器(13),容器底部设有水泵搅拌器(11),液体内设有体核温度传感器(8),其输出端引出容器外,与数据采集装置连接。
7.无创体核温度测量装置标定设备的标定方法,其特征在于:用于标定由本发明装置无创体核温度测量装置绝热材料层热阻Re和热容量Cs决定的瞬态项系数γ标=Re×Cs,和用于校准以本发明方法所述无创体核温度测量方法估算的体核温度变化值的误差补偿系数β标,具体过程为:
(1)、维持标定设备容器内核心温度不变,将本发明无创体核温度测量装置绑缚在本发明标定设备容器的外壁上,稳定后测出模拟体表温度(Ts):
(2)、维持标定设备容器内核心温度不变,用微电热片模拟环境升温ΔTeh,检测模拟体表温度变化为(ΔTsh);
(3)、微电热片停止加热,稳定后,将某点作为基准点,记录该点模拟体表温度值、环境温度值和容器内核心温度值,分别记为Tso、Teo、Tco;
(4)、以基准点为起始点,使标定设备容器内的温度模拟人体体核温度变化,并连续记录和存储模拟体表温度(Ts)、环境温度(Te)和核心温度值(Tc);
(5)、多次重复上述(1)-(4)步骤,并将每次得到的一组测量数据保存待用;
然后包含下述内容:
(6)、对前述每组测量数据分别进行处理,分别取得每组数据的个体化校正系数A=Rc/Re=ΔTsh/(ΔTeh-ΔTsh)和环境校正系数α=ΔTsh/ΔThe;
(7)根据上述的每组数据,用本发明估算方法估算容器内的模拟体核温度Tc,将估算的Tc值与实测的标定设备容器核心温度值做对比,观察二者之间的误差,同时人为调节瞬态项系数中的Re×Cs值,取使瞬态误差最小时的Re×Cs值作为该组测量数据的瞬态系数γ值;
求各组瞬态系数γ值的平均值γ标,将其作为用本发明装置和方法估算人体体核温度时的瞬态项系数Re×Cs的标定值;
(8)、根据上述的每组数据,求出各测量数据Tc、Ts、Te相对于其基准点Tco、Tso、Teo的变化量;
计算环境温度变化造成的体表温度变化量Tse:
Tse=α×(Te-Teo)
计算由单纯核心温度变化造成的体表温度变化量:
dTsc=(Ts-Tso)-α×(Te-Teo);
用本发明方法估算核心温度变化量dTc,将估算的核心温度变化量dTc与实测的对应核心温度变化量相比,得到由体表温度变化量估算体核温度变化量的误差补偿系数β:
β=体核温度变化估算值/体核温度变化实测值=dTc/(Tc实-Tco);
求各组误差补偿系数β的平均值作为用本发明方法估算人体体核温度时的误差补偿系数β标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101067969A CN101112306B (zh) | 2007-06-27 | 2007-06-27 | 无创体核温度测量的方法、装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2007101067969A CN101112306B (zh) | 2007-06-27 | 2007-06-27 | 无创体核温度测量的方法、装置 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101228197A Division CN101843476B (zh) | 2007-06-27 | 2007-06-27 | 无创体核温度测量装置的标定设备及其使用方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101112306A true CN101112306A (zh) | 2008-01-30 |
CN101112306B CN101112306B (zh) | 2010-09-01 |
Family
ID=39020944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007101067969A Expired - Fee Related CN101112306B (zh) | 2007-06-27 | 2007-06-27 | 无创体核温度测量的方法、装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101112306B (zh) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104434048A (zh) * | 2014-12-01 | 2015-03-25 | 电子科技大学 | 一种人体深层组织温度测量仪及测量方法 |
CN105387952A (zh) * | 2015-10-20 | 2016-03-09 | 北京航空航天大学 | 一种间接测量石英振梁加速度计内部石英谐振梁的温度的方法 |
CN106037683A (zh) * | 2015-07-02 | 2016-10-26 | 中国农业科学院北京畜牧兽医研究所 | 一种通过牛体表温度计算牛体温的方法 |
CN106706165A (zh) * | 2015-11-16 | 2017-05-24 | ***通信集团公司 | 一种温度测量的方法及装置 |
CN106798545A (zh) * | 2017-03-03 | 2017-06-06 | 董云鹏 | 体温监测*** |
CN107049253A (zh) * | 2017-04-10 | 2017-08-18 | 深圳市共进电子股份有限公司 | 一种基于人工智能的红外热成像体温检测方法和装置 |
CN108871597A (zh) * | 2018-06-13 | 2018-11-23 | 苏州飞康生物医药有限公司 | 一种智能温度计、智能测温*** |
CN109008989A (zh) * | 2018-06-14 | 2018-12-18 | 杭州感到科技有限公司 | 腹部核温的测量方法和设备 |
CN109186805A (zh) * | 2018-09-12 | 2019-01-11 | 刘凯 | 一种体温计 |
CN109752117A (zh) * | 2017-11-02 | 2019-05-14 | 上海智杏投资管理合伙企业(有限合伙) | 连续测温的体温计 |
CN110381816A (zh) * | 2017-04-04 | 2019-10-25 | 奥尼欧有限公司 | 用于连续且无线监测和分析生物体温度的传感器***和方法 |
CN111060218A (zh) * | 2019-12-20 | 2020-04-24 | 浙江智柔科技有限公司 | 体温测量装置及测量方法 |
CN111141420A (zh) * | 2020-02-04 | 2020-05-12 | 上海申矽凌微电子科技有限公司 | 基于热流法的物体深部温度测量方法及装置 |
CN111693177A (zh) * | 2020-06-23 | 2020-09-22 | 广东小天才科技有限公司 | 体核温度测量方法、装置、介质和智能可穿戴设备 |
CN111741710A (zh) * | 2017-12-28 | 2020-10-02 | 罗伯特·博世有限公司 | 核心温度检测***和方法 |
CN111829665A (zh) * | 2020-07-16 | 2020-10-27 | 深兰自动驾驶研究院(山东)有限公司 | 一种基于用人体作为虚拟黑体的体温测量的方法、装置及存储介质 |
CN112169098A (zh) * | 2020-09-27 | 2021-01-05 | 湖南明康中锦医疗科技发展有限公司 | 一种呼吸支持设备传感器自校准方法和*** |
CN112487692A (zh) * | 2020-12-16 | 2021-03-12 | 中国科学院重庆绿色智能技术研究院 | 一种从额头温度估计体核温度的方法及其应用 |
CN112539839A (zh) * | 2020-10-31 | 2021-03-23 | 中国农业科学院农业信息研究所 | 一种非接触式动物体温监测装置 |
CN113551808A (zh) * | 2021-08-24 | 2021-10-26 | 中国人民解放军总医院 | 核心体温监测实验设备 |
CN113820042A (zh) * | 2021-09-09 | 2021-12-21 | 超越科技股份有限公司 | 一种基于国产化单片机的疫情防控实时测温*** |
CN114206208A (zh) * | 2019-06-11 | 2022-03-18 | F2D医疗公司 | 用于确定人的深层内部温度的温度测量设备和*** |
CN116026493A (zh) * | 2022-06-30 | 2023-04-28 | 荣耀终端有限公司 | 核心体温检测方法及电子设备 |
US11680852B2 (en) | 2020-05-22 | 2023-06-20 | Eaton Intelligent Power Limited | Temperature measurement system |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1249691A1 (en) * | 2001-04-11 | 2002-10-16 | Omron Corporation | Electronic clinical thermometer |
CN2776291Y (zh) * | 2005-02-08 | 2006-05-03 | 广东威尔医学科技股份有限公司 | 体内测温装置 |
-
2007
- 2007-06-27 CN CN2007101067969A patent/CN101112306B/zh not_active Expired - Fee Related
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104434048B (zh) * | 2014-12-01 | 2016-08-24 | 电子科技大学 | 一种人体深层组织温度测量仪及测量方法 |
CN104434048A (zh) * | 2014-12-01 | 2015-03-25 | 电子科技大学 | 一种人体深层组织温度测量仪及测量方法 |
CN106037683A (zh) * | 2015-07-02 | 2016-10-26 | 中国农业科学院北京畜牧兽医研究所 | 一种通过牛体表温度计算牛体温的方法 |
CN105387952A (zh) * | 2015-10-20 | 2016-03-09 | 北京航空航天大学 | 一种间接测量石英振梁加速度计内部石英谐振梁的温度的方法 |
CN105387952B (zh) * | 2015-10-20 | 2018-06-12 | 北京航空航天大学 | 一种间接测量石英振梁加速度计内部石英谐振梁的温度的方法 |
CN106706165B (zh) * | 2015-11-16 | 2019-06-25 | ***通信集团公司 | 一种温度测量的方法及装置 |
CN106706165A (zh) * | 2015-11-16 | 2017-05-24 | ***通信集团公司 | 一种温度测量的方法及装置 |
CN106798545A (zh) * | 2017-03-03 | 2017-06-06 | 董云鹏 | 体温监测*** |
CN106798545B (zh) * | 2017-03-03 | 2019-10-01 | 董云鹏 | 体温监测*** |
CN110381816A (zh) * | 2017-04-04 | 2019-10-25 | 奥尼欧有限公司 | 用于连续且无线监测和分析生物体温度的传感器***和方法 |
WO2018188147A1 (zh) * | 2017-04-10 | 2018-10-18 | 深圳市共进电子股份有限公司 | 基于人工智能的红外热成像体温检测方法、装置和终端设备 |
CN107049253A (zh) * | 2017-04-10 | 2017-08-18 | 深圳市共进电子股份有限公司 | 一种基于人工智能的红外热成像体温检测方法和装置 |
CN109752117A (zh) * | 2017-11-02 | 2019-05-14 | 上海智杏投资管理合伙企业(有限合伙) | 连续测温的体温计 |
CN111741710B (zh) * | 2017-12-28 | 2023-09-15 | 罗伯特·博世有限公司 | 核心温度检测***和方法 |
CN111741710A (zh) * | 2017-12-28 | 2020-10-02 | 罗伯特·博世有限公司 | 核心温度检测***和方法 |
CN108871597A (zh) * | 2018-06-13 | 2018-11-23 | 苏州飞康生物医药有限公司 | 一种智能温度计、智能测温*** |
CN109008989A (zh) * | 2018-06-14 | 2018-12-18 | 杭州感到科技有限公司 | 腹部核温的测量方法和设备 |
CN109008989B (zh) * | 2018-06-14 | 2024-06-11 | 杭州感到科技有限公司 | 腹部核温的测量方法和设备 |
CN109186805A (zh) * | 2018-09-12 | 2019-01-11 | 刘凯 | 一种体温计 |
CN114206208A (zh) * | 2019-06-11 | 2022-03-18 | F2D医疗公司 | 用于确定人的深层内部温度的温度测量设备和*** |
CN111060218A (zh) * | 2019-12-20 | 2020-04-24 | 浙江智柔科技有限公司 | 体温测量装置及测量方法 |
CN111060218B (zh) * | 2019-12-20 | 2021-05-11 | 浙江智柔科技有限公司 | 体温测量装置及测量方法 |
CN111141420A (zh) * | 2020-02-04 | 2020-05-12 | 上海申矽凌微电子科技有限公司 | 基于热流法的物体深部温度测量方法及装置 |
US11680852B2 (en) | 2020-05-22 | 2023-06-20 | Eaton Intelligent Power Limited | Temperature measurement system |
CN111693177A (zh) * | 2020-06-23 | 2020-09-22 | 广东小天才科技有限公司 | 体核温度测量方法、装置、介质和智能可穿戴设备 |
CN111829665A (zh) * | 2020-07-16 | 2020-10-27 | 深兰自动驾驶研究院(山东)有限公司 | 一种基于用人体作为虚拟黑体的体温测量的方法、装置及存储介质 |
CN112169098A (zh) * | 2020-09-27 | 2021-01-05 | 湖南明康中锦医疗科技发展有限公司 | 一种呼吸支持设备传感器自校准方法和*** |
CN112539839A (zh) * | 2020-10-31 | 2021-03-23 | 中国农业科学院农业信息研究所 | 一种非接触式动物体温监测装置 |
CN112487692B (zh) * | 2020-12-16 | 2023-08-22 | 中国科学院重庆绿色智能技术研究院 | 一种从额头温度估计体核温度的方法及其应用 |
CN112487692A (zh) * | 2020-12-16 | 2021-03-12 | 中国科学院重庆绿色智能技术研究院 | 一种从额头温度估计体核温度的方法及其应用 |
CN113551808A (zh) * | 2021-08-24 | 2021-10-26 | 中国人民解放军总医院 | 核心体温监测实验设备 |
CN113551808B (zh) * | 2021-08-24 | 2024-05-14 | 中国人民解放军总医院 | 核心体温监测实验设备 |
CN113820042A (zh) * | 2021-09-09 | 2021-12-21 | 超越科技股份有限公司 | 一种基于国产化单片机的疫情防控实时测温*** |
CN116026493A (zh) * | 2022-06-30 | 2023-04-28 | 荣耀终端有限公司 | 核心体温检测方法及电子设备 |
CN116026493B (zh) * | 2022-06-30 | 2024-04-05 | 荣耀终端有限公司 | 核心体温检测方法及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN101112306B (zh) | 2010-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101112306A (zh) | 无创体核温度测量的方法、装置和标定设备及其标定方法 | |
AU2020300996B2 (en) | Apparatus, systems, and methods for non-invasive thermal interrogation | |
US4386604A (en) | Determination of the basal metabolic rate of humans with a whole-body calorimeter | |
Valvano et al. | Thermal conductivity and diffusivity of biomaterials measured with self-heated thermistors | |
CN101843476B (zh) | 无创体核温度测量装置的标定设备及其使用方法 | |
Patel et al. | A self-heated thermistor technique to measure effective thermal properties from the tissue surface | |
EP3769104B1 (en) | A phantom for multi-parametric calibration in magnetic resonance imaging | |
CN106706165B (zh) | 一种温度测量的方法及装置 | |
GB2474773A (en) | Determination of thermal properties of solid bodies using a flexible membrane serving as both heater and temperature sensor | |
CN110375883B (zh) | 基于主动热流控制的体温计及其测温方法 | |
EP3666179A1 (en) | Core body temperature sensor system based on flux measurement | |
CN109875753A (zh) | 一种缓解儿童腹部不适的中药热敷装置 | |
CN103675018B (zh) | 一种织物热性能测试仪 | |
CN104792439B (zh) | 温度测量方法、装置、探头和*** | |
Saurabh et al. | Continuous core body temperature estimation via SURFACE temperature measurements using wearable sensors-is it feasible? | |
CN111141420A (zh) | 基于热流法的物体深部温度测量方法及装置 | |
Fang et al. | Optimization of a wearable device for core body temperature monitoring based on the dual-heat-flux model | |
CN207832721U (zh) | 一种纤维集合体蓄热性能的测试装置 | |
CN108088871A (zh) | 一种纤维集合体蓄热性能的测试装置及其测试方法 | |
CN212904622U (zh) | 一种测量有限厚度多孔材料中异质含量的装置 | |
CN220932223U (zh) | 核心体温监测实验设备 | |
CN113551808B (zh) | 核心体温监测实验设备 | |
CN113532678A (zh) | 核心体温测量装置 | |
CN102458083B (zh) | 嵌入式保温箱及保温控制板 | |
Deng et al. | Effect of fixing material on skin-contact temperature measurement by wearable sensor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100901 Termination date: 20160627 |
|
CF01 | Termination of patent right due to non-payment of annual fee |