CN101843476B - Calibration facility of non-invasive core temperature measuring device and use method thereof - Google Patents
Calibration facility of non-invasive core temperature measuring device and use method thereof Download PDFInfo
- Publication number
- CN101843476B CN101843476B CN2010101228197A CN201010122819A CN101843476B CN 101843476 B CN101843476 B CN 101843476B CN 2010101228197 A CN2010101228197 A CN 2010101228197A CN 201010122819 A CN201010122819 A CN 201010122819A CN 101843476 B CN101843476 B CN 101843476B
- Authority
- CN
- China
- Prior art keywords
- core temperature
- temperature
- value
- container
- variable quantity
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
The invention relates to a calibration facility of a non-invasive core temperature measuring device and a use method thereof. The calibration facility of the invention is provided with a container (9) filled with liquid medium (12), wherein an electric heater (13) is arranged in the container, the bottom of the container is provided with a water pump blender (11), and a core temperature sensor (8) is arranged in the liquid; and the output end of the core temperature sensor extends out of the container and is connected with a data acquisition device. The facility can be used for calibrating the transient term coefficient gamma<calibration> depending on the heat resistance Re and the thermal capacity Cs of the heat insulating material layer of the non-invasive core temperature measuring device, and can also be used for correcting the error compensation coefficient beta<calibration> of the temperature variation value estimated by using the facility, wherein the gamma<calibration> = Re*Cs. By adopting the calibration facility and the use method thereof, the measuring accuracy of human body core temperature can be improved.
Description
This case is that application number is 200710106796.9, and the applying date is: on June 27th, 2007, name is called dividing an application of " method of non-invasive core temperature measuring, device and calibration facility and scaling method thereof "
Technical field
The invention belongs to a kind of caliberating device and scaling method that obtains human body dynamic physiological parameter device, specifically is a kind of caliberating device and method for using thereof of non-invasive core temperature measuring device.。
Background technology
The dynamic variation property of human body core temperature is comprising the information of thermoregulation system very horn of plenty, is important vital sign, and it can be used as a kind of intermediate parameters is that the comprehensive diagnos of the multiple disease of human body provides reference frame.Continuous, the kinetic measurement branch of core temperature have wound and noinvasive two class methods, have the wound method to bring certain damage to human body, for the occasion of the long-term continuous monitoring of needs, have the wound method that very big inconvenience is arranged.Noinvasive measurement of bldy temperature method has heat insulation method, zero heat flow method etc.After heat insulation method hypothesis body surface was done insulation, skin temperature equaled core temperature.The heat insulation method testing time is long, and precision receives very big restriction because of the difference of adiabator performance.The basic ideas of zero heat flow method are to send hot-fluid with the heat radiation of compensation body surface at body surface to human body, when hot-fluid is zero, think that skin temperature equals core temperature.Zero heat flow method needs heating, chiller to come the balance hot-fluid, and device is comparatively complicated.Also have the scholar to carry out such as infrared method microwave method, ultrasonic method, the research of nuclear magnetic resonance, method measuring body nuclear temperature in addition.These methods all have certain applications condition separately.
Summary of the invention
The technical problem that the present invention will solve provides a kind of calibration facility and method for using thereof of non-invasive core temperature measuring device, adopts present device and method can improve the accuracy of measuring the human body core temperature with non-invasive core temperature measuring device.
The method that core temperature of the present invention is measured comprises following content:
1, at body surface insulation material layer is set; The shell temperature pick off is set between insulation material layer inboard and the body surface; Be used for human body shell temperature (Ts); At insulation material layer arranged outside environment temperature sensor, be used for testing environment temperature (Te), and be provided for producing the electroheat member of simulated environment temperature;
2, carry out the analysis and the calculating of diabatic process with analog circuit, comprise following process:
With electroheat member simulated environment intensification (Δ Te); Record skin temperature variable quantity (Δ Ts) this moment; As personalized correction coefficient (A), and the ratio that obtains variation of ambient temperature amount (Δ Te) and skin temperature variable quantity (Δ Ts) is as the ambient temperature correction coefficient (α) of ambient temperature to the skin temperature influence with this ratio that obtains body surface tissue equivalent thermal resistance (Rc) and said insulation material layer equivalent thermal resistance (Re);
Get the core temperature value (Tc in a certain moment of human body with the electric body-temperature instrumentation
0), record shell temperature value (Ts this moment simultaneously
0) and ambient temperature value (Te
0), datum mark as a reference;
Promptly begin to measure estimation after accomplishing said process: measure body surface temperature value (Ts) and ambient temperature (Te); Equivalent thermal resistance (Re), thermal capacity (Cs) according to said insulation material layer; With aforementioned individuation correction coefficient of obtaining (A) and ambient temperature positive coefficient (α); Be available transient simulation circuit analysis method estimation human body core temperature (Tc), in estimation, proofread and correct body surface tissue equivalent's thermal resistance (Rc) with individuation correction coefficient (A); Remove the influence of ambient temperature (Te) with ambient temperature correction coefficient (α) to skin temperature (Ts), in the estimation with the core temperature value (Tc of said reference point
0) calibrate.
The device of non-invasive core temperature measuring of the present invention is: be provided with the binding strip (1) that can be fixed on the human body; Insulation material layer (2) is set on the binding strip; The position that body surface is being pasted in insulation material layer inside is provided with shell temperature pick off (3); The insulation material layer outside is provided with environment temperature sensor (4) and is used to produce the electroheat member (5) of simulated environment temperature, and described shell temperature pick off is connected with the signals collecting display device with the environment temperature sensor outfan.
Do an explanation in the face of the principle of the inventive method down, the inventive method adopts the analysis and the calculating of breadboardin diabatic process: referring to
Fig. 2 the inventive method is simplified thermal resistance analysis figure with the steady heat conduction of breadboardin diabatic process
Simplification breadboardin figure during steady heat conduction that Fig. 3 the inventive method is adopted
Simplification breadboardin figure when the increase transient state transient term heat conduction that Fig. 4 the inventive method is adopted changes.
Can analyze according to the foregoing circuit simulation drawing: skin temperature transducer is measured body surface temperature value Ts and is received four kinds of factor affecting, i.e. core temperature Tc, ambient temperature Te, adiabator physical property (thermal resistance Re, thermal capacity Cs) and the equivalent thermal resistance Rc of body surface organization.The simplification mathematical relationship of Tc and above-mentioned four factors, figure can be derived as (according to Kirchhoff's current law (KCL)) by transient state heat conduction breadboardin:
Tc＝Rc×Cs×dTs/dt+Rc/Re×(Ts-Te)+Ts (1)
Tc and the simplification mathematical relationship of above-mentioned four factors in the steady heat conduction stage are:
Tc＝Rc/Re×(Ts-Te)+Ts (2)
Above-mentioned formula (1), (2) promptly are the basic representations that the inventive method adopts transient simulation circuit analysis method estimation human body core temperature (Tc); Can find out from above-mentioned formula (1): body surface temperature T s and ambient temperature Te can record through pick off; Insulation material layer thermal resistance Re, thermal capacity Cs (depending primarily on density p and specific heat capacity c) can artificially design, and have only the thermal resistance Rc of body surface tissue equivalent to be difficult for obtaining, and Rc also can have more different because of the individual variation of factors such as the body surface fat of different human body, muscle; In addition; The variation of body surface temperature T s is that human body core temperature Tc and ambient temperature Te change coefficient result, if ambient temperature Te is constant, and the then just directly dynamic change of antimer nuclear temperature Tc of the variation of Ts; Therefore, need from the variation of Ts, remove the influence of variation of ambient temperature.
Therefore the inventive method is before measuring estimation, and first simulated environment variations in temperature Δ Te records the skin temperature variation delta Te after stablizing, and can try to achieve according to formula (2)
Rc/Re＝ΔTs/(ΔTe-ΔTs)＝A
With the personalized correction coefficient of A as the embodiment body surface thermal resistance Rc of tissue equivalent;
The Δ Ts/ Δ Te=α that can try to achieve simultaneously is as embodying the ambient temperature correction coefficient of from skin temperature changes, removing ambient temperature effect.
Can realize like this (TC) with increasing transient term Analog Circuits Analysis method estimation human body core temperature.
The calibration facility structure that the present invention measures human body core temperature device is: establish the container (9) that liquid medium (12) are housed; Be provided with electric heater (13) in the container; Container bottom is provided with water pump agitator (11); Be provided with core temperature pick off (8) in the liquid, its outfan is drawn outside the container, is connected with the data acquisition display device.
The scaling method of calibration facility of the present invention is: be used to demarcate the transient term coefficient gamma by apparatus of the present invention non-invasive core temperature measuring device insulation material layer thermal resistance Re and thermal capacity Cs decision
Mark=Re * Cs and the error compensation factor beta that is used to calibrate the core temperature changing value of estimating with the inventive method
Mark, detailed process is:
1, it is temperature-resistant to keep calibration facility container inner core, non-invasive core temperature measuring device of the present invention is strapped on the outer wall of calibration facility container of the present invention, measures simulation shell temperature (Ts) after stable;
2, it is temperature-resistant to keep calibration facility container inner core, and with little electric heating piece simulated environment intensification Δ Teh, detecting simulation body surface variations in temperature is Δ Tsh;
3, little electric heating piece stops heating, after stablizing, certain point as datum mark, is write down this some simulation body surface temperature value, ambient temperature value and container inner core temperature value, is designated as respectively (Tso, Teo, Tco);
4, be starting point with the datum mark, the temperature simulation human body core temperature in the calibration facility container is changed, and continuous record and storage simulation shell temperature (Ts), ambient temperature (Te) and core temperature value (Tc);
5, repeatedly repeat above-mentioned 1-4 step, and one group of measurement data that will at every turn obtain is preserved for use; Comprise following content then:
6, aforementioned every group of measurement data handled respectively, obtained the individuation correction coefficient A=Rc/Re=Δ Tsh/ (Δ Teh-Δ Tsh) and the environment correction factor α=Δ Tsh/ Δ The of every group of data respectively;
7, according to every group of above-mentioned data; With the analogue body nuclear temperature Tc in the inventive method estimation container; The Tc value of estimation is done contrast with the calibration facility container core temperature value of actual measurement; Observe the error between the two, simultaneously artificial Re * Cs value of regulating in the transient term coefficient is got and is made the transient state coefficient γZhi of transient error Re * Cs value hour as this group measurement data;
Ask the meansigma methods γ that respectively organizes transient state coefficient γZhi
Mark, with the calibration value of its transient term coefficients R e * Cs as with apparatus of the present invention and method estimation human body core temperature the time;
8,, obtain each measurement data Tc, Ts, Te variable quantity with respect to its datum mark Tco, Tso, Teo according to every group of above-mentioned data;
The shell temperature variable quantity Tse that the computing environment variations in temperature causes:
Tse＝α×(Te-Teo)
Calculating is changed the shell temperature variable quantity that causes by simple core temperature:
dTsc＝(Ts-Tso)-α×(Te-Teo)；
With the inventive method estimation core temperature variable quantity dTc, the core temperature variable quantity dTc that estimates is compared with the corresponding core temperature variable quantity of actual measurement, obtain error compensation factor beta by shell temperature variable quantity estimation core temperature variable quantity:
β=core temperature changes estimated value/core temperature and changes measured value=dTc/ (Tc
Real-Tco);
The error compensation factor beta of the meansigma methods of asking each grouping error penalty coefficient β as the time with the inventive method estimation human body core temperature
Mark
The invention provides a kind of noinvasive human body core temperature new measurement method and new equipment; Adopt the inventive method and device to carry out long-term, continuous, dynamic, personalized non-invasive core temperature measuring to human body; Easy to use, manufacturing and operating cost are low, and accuracy is high.Adopt calibration facility of the present invention and its scaling method can improve the accuracy that changes with the inventive method and measurement device human body core temperature.
Description of drawings
Fig. 1, heat insulation method measuring body nuclear temperature block diagram
Fig. 2, the inventive method are simplified thermal resistance analysis figure with the steady heat conduction of breadboardin diabatic process
Simplification breadboardin figure when Fig. 3, steady heat conduction that the inventive method adopted
Fig. 4, the simplification breadboardin figure when the increase transient state heat conduction that the inventive method adopted changes.
Fig. 5, apparatus of the present invention frame structure sketch map
Fig. 6, embodiment 1 apparatus of the present invention structural representation
Fig. 7, embodiment 3 calibration facility structural representations of the present invention
Fig. 8, embodiment 5 transient state coefficient gamma
MarkDemarcate sketch map
Fig. 9, embodiment 5 error compensation factor beta are demarcated sketch map
1-binding strip 2-insulation material layer 3-shell temperature pick off 4-environment temperature sensor 5-heating 6-elastic webbing 7-locker 8-core temperature pick off 9-container 10-lid 11-stirs pump 12-medium 13-electrical heating elements 14-measuring disk
Specific embodiments:
Referring to Fig. 6, be an embodiment of core temperature measuring device of the present invention, this device is provided with the measuring disk 14 of band binding strip 1, and binding strip is processed by soft cloth, and flexible material 6 and locker 7 are to guarantee the degree of frapping to human body on it.Measuring disk 14 is provided with insulation material layer 2; Shell temperature pick off 3 is equipped with at the position that human body skin is being pasted in its inboard, is embedded with environment temperature sensor 4 and the heating that is connected with battery 5 near the insulation material layer in the outside, and measuring disk is cuboid; Be of a size of 40*40*15mm3 (the wide * of long * is high); Heavy 10g, heating 5 is the soft thin slice of little electric heating of 35*35*1mm, is embedded in the adiabator.Isolated by insulation material layer 2 between little electric heating piece and the shell temperature pick off 3; The two vertical dimension is 10mm; Shell temperature pick off 3 can pass through good thermal conductor and closely contact with skin surface; Environment temperature sensor 4 contacts with little electric heating piece 5 good heat conduction, and shell temperature pick off 3, environment temperature sensor 4, little electric heating piece heating 5 overlap with the measuring disk center, and shell temperature pick off 3 is connected with the temperature signal collection device respectively with environment temperature sensor 4 outfans.
This example is fixed on embodiment 1 said device on the human body for an embodiment of core temperature measuring method of the present invention, can carry out the core temperature estimation, and detailed process is following:
1. the measuring point is arrived in the measuring device bondage, stablized 30 minutes;
2, with little electric heating piece heating 5 simulated environment intensification Δ Teh, stablize 20 minutes after, the detection of skin variations in temperature is Δ Tsh, then can get:
Individuation correction coefficient: A=Rc/Re=Δ Tsh/ (Δ Teh-Δ Tsh)
Ambient temperature correction coefficient: α=Δ Tsh/ Δ Teh
3. after stopping heating, write down the skin temperature Ts and the ambient temperature Te in a certain moment, use electric body-temperature instrumentation amount core temperature (rectum or axillary fossa) simultaneously and note, be designated as Tso, Teo, Tco respectively, as a reference datum mark;
4, begin to carry out the estimation of core temperature:
Gather different the shell temperature Ts and the ambient temperature Te of human body constantly, adopt formula 3 to carry out the estimation of body nuclear degree:
Tc＝A×Re×Cs×[(dTs-α×dTe)/dt]+
A×[(Tso+dTs-α×dTe)-(Teo+dTe)]+[(Tso+dTs)-α×dTe]
(3)
In the formula: Tc is a human body core temperature estimated value
Cs is the thermal capacity of insulation material layer
The body surface temperature value of Ts for recording immediately
The ambient temperature value of Te for recording immediately
Re is the thermal resistance of insulation material layer
The personalized correction coefficient of A for having obtained through environment temperature rise Δ Teh
A-Rc/Re＝ΔTsh/(ΔTeh-ΔTsh)
Ambient temperature correction coefficient alpha=Δ Tsh/ Δ the Teh of α for obtaining through environment temperature rise Δ Teh
Core temperature value Tco with reference point does single-point calibration to the value of estimating with formula 3 then; Be about to the Tco value that the above-mentioned by formula initial point value of human body core temperature Tc and the datum mark of 3 estimations survey and do the difference comparison, the difference that obtains is as the correcting value of follow-up core temperature estimated value.
In above-mentioned estimation equation; Heat insulator thermal resistance Re, thermal capacity Cs can artificially design (for example through changing the density p and the specific heat capacity c of insulation material layer); A and α also obtain; (dTs-α * dTe) for to have removed the skin temperature variable quantity after the ambient temperature effect, so can realize according to Ts that records and Te estimation core temperature Tc.
This example is another embodiment of core temperature measuring method of the present invention
The difference of this method and embodiment 2 is that the evaluation method of the 4th step core temperature is different; This example is tried to achieve individuation correction coefficient A; Behind the Tso of ambient temperature correction coefficient alpha and reference point, Teo, the Tco, carry out the estimation of core temperature Tc by following formula:
The steady-state quantity estimation:
Tc
Surely=Tco+ (dTs-α * dTe)/(1-α) (4)
The transient volume estimation:
Tc
Wink=A * Re * Cs * (dTs-α * dTe)/dt (5)
Then the core temperature estimated value is: Tc=Tc
Surely+ Tc
Wink
In the formula
Tc is a human body core temperature estimated value
Cs is the thermal capacity of insulation material layer
The body surface temperature value of Ts for recording immediately
The ambient temperature value of Te for recording immediately
Re is the thermal resistance of insulation material layer
The personalized correction coefficient of A for obtaining through environment temperature rise Δ Teh
A＝Rc/Re＝ΔTsh/(ΔTeh-ΔTsh)
The ambient temperature correction coefficient of α for obtaining through environment temperature rise Δ Teh
α＝ΔTsh/ΔTeh
Core temperature value Tco with reference point does single-point calibration to the value of estimating with formula 3 then; Be about to the initial point value of human body core temperature Tc of above-mentioned estimation and Tco value that datum mark is surveyed and do difference relatively, the difference that obtains is as the correcting value of follow-up core temperature estimated value.
The estimation principle of said method is:
As aforementioned: individuation correction coefficient A=Rc/Re=dTsh/ (dTeh-dTsh)
Ambient temperature correction coefficient alpha=dTsh/dTeh.
So: the relation of A and α is: A=α/(1-α)
When ambient temperature is constant, formula 1
Tc＝Rc×Cs×dTs/dt+Rc/Re×(Ts-Te)+Ts
＝A×Re×Cs×[(dTs-α×dTe)/dt]+A×(Ts-Te)+Ts
Being varied to of core temperature is divided in the middle steady-state quantity:
dTc＝A×dTs+dTs＝α/(1-α)×dTs+dTs＝dTs/(1-α)
Remove the influence of ambient temperature, promptly get skin temperature:
dTc＝(dTs-α×dTe)/(1-α)(6)
So can get formula 5 steady-state quantities
Tc
Surely=Tco+ (dTs-α * dTe)/(1-α) (4)
And according to formula (3), the transient volume estimation:
Tc
Wink=A * Re * Cs * [(dTs-α * dTe)/dt] (5)
Wherein Re * Cs can be according to the transient state coefficient of in advance the measuring device insulation material layer being demarcated
γ
Mark=Re * Cs value is calculated
Thereby can obtain the core temperature estimated value is steady-state value and instantaneous value sum:
Tc=Tc
Surely+ Tc
Wink
In order to improve the accuracy of above-mentioned estimated value, further scheme is to use the error compensation factor beta
MarkIt is compensated, and the evaluation method that compensates is:
After aforesaid steady-state quantity being compensated, get through the error compensation factor beta of demarcating
(Tc
Surely)
Mend=Tco+ (dTs-α * dTe)/[β
Mark* (1-α)]
The estimation of transient volume is with aforementioned identical;
Then the core temperature estimated value is:
Tc=((Tc
Surely)
Mend+ Tc
Wink
This example is the embodiment of non-invasive core temperature measuring device calibration facility of the present invention, and this calibration facility is a body temperature simulative generator, and concrete structure is referring to Fig. 7; Container simulation trunk is drum shape; Overall dimension is decided to be the high 250mm of bucket, internal diameter 80mm, and 25mm is thick for the ABS base.Top cover 10 is the thick ABS material of 10mm, leaves cable hole and steam vent on the top cover, and is provided with the seal that supports heating 13.Electrical heating elements is an alternating current 220V, 500W type cylinder steel pipe heating rod, and long 200mm, diameter 10mm is used for the analogue body nuclear heat source.Medium is a water in the bucket, and the immersible pump of 12V4W is fixed on the base, stirs current, makes the inner bucket water temperature even.
The bucket wall of body temperature simulative generator is used for analogue body and examines the thermal resistance of body surface, and its thickness is confirmed as follows.Human body skin heat conductivity=0.47W/ (m.k), fatty heat conductivity=0.21W/ (m.k), normal abdominal part male skin fold (skin+subcutaneous fat) thickness is 5-15mm, the female skin pleat is 12-20mm.This example equal skin depth=10mm that makes even, its mesocuticle is thick=2mm, and fat is thick=8mm, according to the heat conduction steady-state process, entire thermal resistance equals the series resistances sum, that is:
It is λ a=0.236W/ (m.k) that the substitution data computation gets the average heat conductivity of rhicnosis.ABS material thermal conductivity λ abs=0.2512W/ (m.k), the principle that then equates according to above-mentioned thermal resistance equally, the ABS material thickness of 10mm rhicnosis thickness of equal value does
So ABS bucket wall thickness is decided to be 11mm in the model.
Be water temperature in the control body temperature simulative generator, other has designed temperature control panel, and temperature-controlled precision is at ± 0.1 ℃, so that realize the constant of temperature in the body temperature simulative generator and change according to requirement of experiment.
Be provided with analogue body nuclear temperature pick off 8 in the bucket, be used to detect mimic core temperature.
This example is the embodiment of non-invasive core temperature measuring device calibration facility scaling method of the present invention.
This example adopts embodiment 4 described calibration facilities and embodiment 3 described measuring disk, and its temperature sensor outfan is connected with the temperature signal collection panel, is used to gather the temperature signal of calibration facility and measuring disk, and the temperature of control calibration facility container kernel.
For the error that exists in correctly estimating the core temperature transient volume and reduce to estimate; New measuring disk will be demarcated its parameter with calibration facility before use; This method promptly can be used for demarcating non-invasive core temperature measuring device of the present invention, can be used for simultaneously calibrating with non-invasive core temperature measuring device of the present invention and measuring method carrying out the estimated value that the side amount obtains.
The purpose one of this method is the transient state coefficient gamma=Re * Cs that calibrates embodiment 1 said measuring disk, so that make its available embodiment 2 or embodiment 3 said methods carry out the estimation of human body core temperature transient components; The 2nd, obtain with the present invention and carry out the error compensation factor beta that the human body core temperature changes the estimation generation, to compensate the error that estimation exists.
Detailed process is:
1, room temperature is 25 ℃, and simulation body temperature generator core temperature temperature control is at 37 ℃.Measuring disk is strapped on the simulation body temperature generator by shown in Figure 7, stablizes 30 minutes.At this moment, under the little situation of environmental temperature fluctuation (dTe＜1 ℃), shell temperature is near a stationary value;
2, simulation body temperature generator core temperature still temperature control 37 ℃ constant, little electric heating piece simulated environment intensification Δ Teh, the detection of skin variations in temperature is Δ Tsh;
3, little electric heating piece stops heating, stablizes 30 minutes, makes little electric heating piece temperature fall back to the temperature spot before heating up.Then, write down some shell temperature values, ambient temperature value and core temperature value simultaneously, be designated as Tso, Teo, Tco respectively;
4, simulation body temperature generator core temperature with 1 ℃ gradient change, kept 1 hour on each thermograde in 36-41 ℃ of scope.Continuous record body surface, environment and core temperature value are recorded in the computer with Tst, Tet and Tct respectively simultaneously.
5, repeat above-mentioned 1-4, measure 10 times, be kept in the different catalogues the one group of result who at every turn measures for use.
6, according to the above-mentioned data scaling transient state coefficient gamma of respectively organizing:
Coefficient in the middle core temperature of formula (1) among transient term Rc * Cs * dTs/dt should not directly be obtained, but through little electric heating piece simulated environment variations in temperature Teh, can obtain:
A＝Rc/Re＝ΔTsh/(ΔTeh-ΔTsh)
So Rc=A * Re=Δ Tsh/ (Δ Teh-Δ Tsh) * Re
Then transient term is expressed as:
A×Re×Cs×dTs/dt＝ΔTsh/(ΔTeh-ΔTsh)×Re×Cs×dTs/dt
Because A=Δ Tsh/ (Δ Teh-Δ Tsh) part can be tried to achieve by little electric heating piece simulated environment variations in temperature Δ Tsh/ of aforementioned the 2nd step; And Re * Cs partly be one with the relevant amount of the thermophysical property of measuring disk own; In case the physical property and the shape of measuring disk have confirmed that Re * Cs has just confirmed.So Re * Cs also be one quantitatively, as the transient state coefficient gamma, concrete calibration process is with Re * Cs for we:
To 10 groups of experimental datas, do following processing:
According to little electric heating piece simulated environment variations in temperature data, obtain the individuation correction coefficient:
A＝Rc/Re＝ΔTsh/(ΔTeh-ΔTsh)
Obtain environment correction factor
α＝dTsh/dTeh
With formula 3:
Tc＝A×Re×Cs×[(dTs-α×dTe)/dt]+
A×[(Tso+dTs-α×dTe)-(Teo+dTe)]+[(Tso+dTs)-α×dTe]
(3)
Estimation simulation body temperature generator core temperature; And with the actual measurement simulation body temperature generator core temperature do contrast; The artificial simultaneously Re * Cs value of regulating; Observe transient error size between core temperature estimated value and measured value, as shown in Figure 8, get the transient state coefficient γZhi that transient error Re * Cs value hour tests for this group (Tc1 that obtains during Re among the figure * Cs=80 is near measured value).
Such 10 groups of calibration experiment data just obtain 10 transient state coefficient γZhi, and are as shown in table 1.Then the transient state coefficient of this measuring disk is that 10 γ average, γ '=87.1, the calibration coefficient of used formula (3) transient term when estimating when this value promptly can be used as with embodiment 1 measuring disk measurement human body core temperature.
Calibration experiment transient state coefficient table
| |
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | Average |
| γ＝A×Re×Cs | 90 | 85 | 88 | 91 | 86 | 88 | 90 | 84 | 80 | 89 | 87.1 |
7, according to the above-mentioned data scaling error compensation factor beta of respectively organizing
Find in the measurement; In process, because there is certain error in operations such as the choosing of datum mark, ambient temperature effect deduction, so cause having certain error between estimated value and the actual value of core temperature variable quantity dTc by shell temperature variable quantity dTs derivation core temperature variable quantity dTc; For this error is recompensed; This example is to embodiment 1 said measuring disk, and the measuring method of embodiment 2,3 said estimation core temperatures is demarcated an error compensation factor beta, and method is following:
10 groups of experimental datas to aforementioned 1-5 obtains process as follows:
1. according to little electric heating piece simulated environment intensification Δ Teh value and the corresponding shell temperature changing value Δ Tsh of record, can calculate:
The individuation correction coefficient:: A=Rc/Re=Δ Tsh/ (Δ Teh-Δ Tsh)
Ambient temperature correction coefficient: α=Δ Tsh/ Δ Teh
2. treating excess syndrome is tested step 3 data recorded value Tco, and Tso and Teo be respectively as core temperature, the benchmark of shell temperature and variation of ambient temperature (referring to Fig. 9);
3. experiment with computing step 4 experimental data is with respect to the variable quantity (referring to Fig. 9) of benchmark
Core temperature Tc is the variable quantity of its benchmark relatively: Tc-Tco;
Shell temperature Ts is the variable quantity of its benchmark relatively: Ts-Tso;
Ambient temperature Te is the variable quantity of its benchmark relatively: Te-Teo.
4. the shell temperature variable quantity that causes of computing environment variations in temperature: α * (Te-Teo);
5. because the variation of shell temperature is the coefficient result of core temperature and ambient temperature; That is: shell temperature variable quantity=core temperature changes actuating quantity+variation of ambient temperature actuating quantity, changes the shell temperature variable quantity that causes so calculate core temperature here:
ΔTsc＝(Ts-Tso)-α×(Te-Teo)；
6. change the stable state composition that the shell temperature variable quantity dTsc that causes estimates the core temperature variable quantity by simple core temperature: adopt embodiment 3 described estimation formulas 6 to estimate here:
dTc＝(dTs-α×dTe)/(1-α)(6)
The actual measurement stable state core temperature variable quantity that obtains is compared in the core temperature variable quantity dTc that 7. will 6. estimate and aforementioned the 4th step, just obtains the error compensation factor beta by shell temperature variable quantity estimation core temperature variable quantity, promptly
β=core temperature changes estimator/core temperature and changes actual measured amount
=dTc/ (Tc
Real-Tco)
8., 10 groups of calibration experiment data just obtain 10 β values, as shown in table 2,10 β get its meansigma methods and get 0.8835.This is worth as the error compensation factor beta with shell temperature variable quantity estimation core temperature variable quantity
Mark
To use the error compensation factor beta with the core temperature variable quantity that the inventive method estimates
MarkAfter the correction, do contrast with actual measurement core temperature variable quantity again, its error greatly reduces (referring to table 2).
The errors table of table 2 during with error compensation coefficient estimate core temperature variable quantity
| Experiment | β Mark | Use β MarkTc variable quantity after the compensation (℃) | Actual measurement Tc variable quantity (℃) | Use β MarkAverage compensation back Tc variable quantity estimation error (℃) |
| 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 |
| Average | 0.8835 | -0.0047 |
Through above transient state coefficient gamma
MarkWith the error compensation factor beta
MarkDemarcation, calibration process just is through with, and carries out actual measurement estimation and with the error compensation coefficient estimated value is compensated improving the accuracy of estimating the result with calibrated measuring disk.
Claims (2)
1. the calibration facility of non-invasive core temperature measuring device; It is characterized in that: establish the container (9) that liquid medium (12) are housed; Be provided with electric heater (13) in the container, container bottom is provided with water pump agitator (11), is provided with core temperature pick off (8) in the liquid; Its outfan is drawn outside the container, is connected with data acquisition unit.
2. the method for using of the calibration facility of non-invasive core temperature measuring device according to claim 1 is characterized in that: the transient term coefficient gamma that is used to demarcate non-invasive core temperature measuring device insulation material layer thermal resistance Re and thermal capacity Cs decision
Mark=Re * Cs and being used to calibrates the error compensation factor beta with this device estimation core temperature changing value
Mark, detailed process is:
(1), to keep calibration facility container inner core temperature-resistant, non-invasive core temperature measuring device is strapped on the outer wall of calibration facility container, measures simulation body surface temperature T s after stable;
(2), to keep calibration facility container inner core temperature-resistant, with little electric heating piece simulated environment intensification Δ Teh of non-invasive core temperature measuring device, detecting simulation body surface variations in temperature is Δ Tsh;
(3), little electric heating piece stops heating, stable after, certain point as datum mark, is write down this point and simulates body surface temperature value, ambient temperature value and container inner core temperature value, be designated as Tso, Teo, Tco respectively;
(4), be starting point with the datum mark, the temperature simulation human body core temperature in the calibration facility container is changed, and continuous record is simulated body surface temperature T s, ambient temperature Te and core temperature value Tc with storage;
(5), repeatedly repeat above-mentioned (1)-(4) step, and one group of measurement data that will at every turn obtain is preserved for use;
Comprise following content then:
(6), aforementioned every group of measurement data handled respectively, obtain respectively every group of data individuation correction coefficient A=Rc/Re=Δ Tsh/ (Δ Teh-Δ Tsh) and
Environment correction factor α=Δ Tsh/ Δ Teh;
Wherein Rc is a body surface tissue equivalent thermal resistance
(7) according to every group of above-mentioned data; Core temperature value Tc in the estimation container; The Tc value of estimation is done contrast with the calibration facility container core temperature value of actual measurement; Observe the error between the two, simultaneously artificial Re * Cs value of regulating in the transient term coefficient is got and is made the transient term coefficient γZhi of transient error Re * Cs value hour as this group measurement data;
Ask the meansigma methods γ that respectively organizes transient term coefficient γZhi
Mark, with the calibration value of its transient term coefficients R e * Cs as with non-invasive core temperature measuring device estimation human body core temperature the time;
(8), according to every group of above-mentioned data, obtain each measurement data Tc, Ts, Te variable quantity with respect to its datum mark Tco, Tso, Teo;
The shell temperature variable quantity Tse that the computing environment variations in temperature causes:
Tse＝α×(Te-Teo)
Calculating is changed the shell temperature variable quantity that causes by simple core temperature:
dTsc＝(Ts-Tso)-α×(Te-Teo)；
With non-invasive core temperature measuring device estimation core temperature variable quantity dTc; The core temperature variable quantity dTc of estimation is compared with the corresponding core temperature variable quantity Tc reality-Tco of actual measurement, obtains error compensation factor beta by shell temperature variable quantity estimation core temperature changing value:
β=core temperature variable quantity estimated value/core temperature variable quantity measured value=dTc/ (Tc
Real-Tco);
The error compensation factor beta of the meansigma methods of asking each grouping error penalty coefficient β as the time with non-invasive core temperature measuring device estimation human body core temperature
Mark
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010101228197A CN101843476B (en) | 2007-06-27 | 2007-06-27 | Calibration facility of non-invasive core temperature measuring device and use method thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010101228197A CN101843476B (en) | 2007-06-27 | 2007-06-27 | Calibration facility of non-invasive core temperature measuring device and use method thereof |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2007101067969A Division CN101112306B (en) | 2007-06-27 | 2007-06-27 | Method and equipment for non-invasive core temperature measuring and the calibration equipment and the calibration method thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101843476A CN101843476A (en) | 2010-09-29 |
| CN101843476B true CN101843476B (en) | 2012-03-28 |
Family
ID=42768441
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010101228197A Expired - Fee Related CN101843476B (en) | 2007-06-27 | 2007-06-27 | Calibration facility of non-invasive core temperature measuring device and use method thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101843476B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102011114620B4 (en) * | 2011-09-30 | 2014-05-08 | Dräger Medical GmbH | Apparatus and method for determining body core temperature |
| EP2997886B1 (en) * | 2014-09-13 | 2018-08-08 | Arc Devices Limited | Non-touch detection of body core temperature |
| WO2019133469A1 (en) * | 2017-12-28 | 2019-07-04 | Robert Bosch Gmbh | Core temperature detection system and method |
| CN109008989A (en) * | 2018-06-14 | 2018-12-18 | 杭州感到科技有限公司 | The measurement method and equipment of abdomen core temperature |
| CN109932503B (en) * | 2019-03-25 | 2020-02-18 | 江苏精策医疗科技有限公司 | Simulation experiment and noninvasive blood glucose detection method based on body heat production-heat dissipation model |
| CN112386231A (en) * | 2019-08-15 | 2021-02-23 | 安徽华米信息科技有限公司 | In-vivo temperature measuring method, device, equipment and storage medium |
| CN113551808B (en) * | 2021-08-24 | 2024-05-14 | 中国人民解放军总医院 | Core body temperature monitoring experimental facility |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1249691A1 (en) * | 2001-04-11 | 2002-10-16 | Omron Corporation | Electronic clinical thermometer |
| CN2776291Y (en) * | 2005-02-08 | 2006-05-03 | 广东威尔医学科技股份有限公司 | Body internal temperature measuring device |
| CN2879180Y (en) * | 2006-03-16 | 2007-03-14 | 东营研博电子科技有限公司 | Programmable crude oil water measuring apparatus |
-
2007
- 2007-06-27 CN CN2010101228197A patent/CN101843476B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1249691A1 (en) * | 2001-04-11 | 2002-10-16 | Omron Corporation | Electronic clinical thermometer |
| CN2776291Y (en) * | 2005-02-08 | 2006-05-03 | 广东威尔医学科技股份有限公司 | Body internal temperature measuring device |
| CN2879180Y (en) * | 2006-03-16 | 2007-03-14 | 东营研博电子科技有限公司 | Programmable crude oil water measuring apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101843476A (en) | 2010-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101112306B (en) | Method and equipment for non-invasive core temperature measuring and the calibration equipment and the calibration method thereof | |
| CN101843476B (en) | Calibration facility of non-invasive core temperature measuring device and use method thereof | |
| KR102630649B1 (en) | Apparatus, systems and methods for non-invasive thermal irradiation | |
| CN104596667B (en) | The sensitivity method of ultrasonic listening interior of articles transient state non-uniform temperature field | |
| Patel et al. | A self-heated thermistor technique to measure effective thermal properties from the tissue surface | |
| CN101736570B (en) | Device and method for testing fabric contact coldness | |
| Alvarado et al. | A hot-wire method based thermal conductivity measurement apparatus for teaching purposes | |
| CN101788513A (en) | Measurement device of thermal conductivity of materials and method thereof | |
| Kun et al. | Accurate flexible temperature sensor based on laser-induced graphene material | |
| CN105548246A (en) | Heat-conductivity-coefficient measuring experiment system through steady-state method and measuring method | |
| CN108107072A (en) | A kind of soil body Determination of conductive coefficients method and experimental rig | |
| CN109758149B (en) | Electrical impedance imaging method combining time difference imaging and frequency difference imaging | |
| Valvano et al. | Analysis of the Weinbaum-Jiji model of blood flow in the Canine Kidney cortex for self-heated thermistors | |
| Alkhwaji et al. | New mathematical model to estimate tissue blood perfusion, thermal contact resistance and core temperature | |
| Saurabh et al. | Continuous core body temperature estimation via SURFACE temperature measurements using wearable sensors-is it feasible? | |
| CN203502367U (en) | Device for testing heat conductivity coefficient of material by transient plane heat source method | |
| CN114526844B (en) | Thermal parameter self-testing method for thermopile sensor | |
| Fang et al. | Optimization of a wearable device for core body temperature monitoring based on the dual-heat-flux model | |
| Oh et al. | Detection of temperature distribution via recovering electrical conductivity in MREIT | |
| Frim et al. | Heat flux transducer measurement error: a simplified view | |
| CN205426383U (en) | Temperature measurement probe and system | |
| CN1254215C (en) | Method of non destructive measuring biological active tissue heat parameter | |
| CN220932223U (en) | Core body temperature monitoring experimental facility | |
| Chudzik | Measuring system with a dual needle probe for testing the parameters of heat-insulating materials | |
| CN113551808B (en) | Core body temperature monitoring experimental facility |
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: 20120328 Termination date: 20160627 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |