CN103070668A - Heart age detector and detection method thereof - Google Patents

Abstract

The invention discloses a heart age detector and a detection method thereof, belonging to the technical field of biomedical engineering. The heart age detector is characterized by comprising a finger artery pulse wave sensor and a sphygmomanometer, wherein the finger artery pulse wave sensor and the sphygmomanometer are connected with a signal acquiring and analyzing meter through connecting cables. According to the heart age detector, by using a non-invasive measurement method, upper arm blood pressure and finger artery pulse wave signals of a human body in the rest and after movement are respectively collected, and heart function parameters are compared and analyzed on the basis of a nonlinear pulse wave communication theory, and thus the purpose of evaluating the heart age is achieved. Through primarily evaluating the heart age and the heart function condition, the heart age detector can provide important reference basis for health, living habits, sports health and the like for people.

Description

Heart age detector and detection method thereof

Technical field:

The invention belongs to the biomedical engineering technology field, relate to a kind of noninvasive detection device for appraiser's systemic heart age.

Background technology:

There is an advertising words to be: " 60 years old heart of 30 years old people, 30 years old heart of 60 years old people ".Its implication behind affords food for thought, and has also proposed a very rigorous problem in science: human heart age and physiological age actually difference have muchly, how to evaluate and test heart age?

Along with the raising of people's living standard, Aging Problem becomes the major issue of global concern.Because increase and other many reasons of life stress, cardiovascular disease become the major disease that threatens human health.The increasing a series of life-threatening diseases such as hypertension, atherosclerosis, heart failure, myocardial infarction that finally cause of the burden of heart.If timely human body heart age so just can be the daily health supervision of people, reference frame for preventing that providing of cardiovascular disease is important.

Summary of the invention:

The present invention combines the noinvasive detection information of pulse wave and blood pressure, utilizes non-linear pulse wave theory to calculate the index of cardiovascular performance parameter, thereby calculates the human heart age, for physical condition provides important references.

A kind of human heart age detector is characterized in that: comprise the upper arm sphygomanometer, refer to affectionately fight wave sensor and cardiac function analysis of information collection instrument, by wire link together (Fig. 1).Can gather simultaneously blood pressure and pulse ripple signal, utilize non-linear pulse wave theory to calculate cardiac functional parameter, measuring and calculating heart age (Fig. 2).

This device respectively human body according to general normal person's physiological parameter scope, is estimated tested heart age under quiescent condition and the blood pressure and pulse ripple signal behind the quantitative exercise.

This device utilizes the noinvasive detection method to carry out heart age evaluation and test, can recuperate etc. for community medicine, rehabilitation, sport and body-building, leisure good reference frame is provided.Instrumentation is simple, and adaptability is large, and broad masses of the people's medical treatment ﹠ health is had good reference significance.

The invention is characterized in, in detector, realize according to the following steps successively:

Step (1), detector initializes:

Input: the radial artery pressure pulse wave that cuff sphygmomanometer detects and measured's pulse pressure value; The pulse waveform that the wave sensor that refers to affectionately to fight gathers;

Measured's personal information comprises sex, age, height, body weight and shoulder breadth at least;

Step (2), detector calculate the average pulse pressure P of measured in the cardiac cycle automatically _m, i.e. the arithmetic mean of instantaneous value of one group of detected value;

Step (3) calculates pressure pulse velocity of wave propagation c (t):

$C\left(t\right)=\sqrt{\frac{P\left(t\right)(1+b*\ln \frac{P\left(t\right)}{\mathrm{Pm}})}{\mathrm{\ρb}}}$ , wherein:

P (t) is pulse pressure value corresponding to each sampled point in the cardiac cycle,

ρ is density of blood, ρ=1.05*10 ^-3Kg/cm ³,

B is the pulse waveform parameter, is calculated as follows:

$b=\frac{a{\left(\frac{0.81}{k+0.48}\right)}^2-1}{\mathrm{\π}}$ ，

$k=\frac{P_m-P_d}{P_s-P_d}$ , wherein:

K is shape factor,

P _d, be diastolic pressure, P _sBe systolic pressure, be measured value,

A is the parameter relevant with velocity profile, a=0.57;

Step (4) calculates barometric gradient F (t):

$F\left(t\right)=\frac{1}{\mathrm{\ρc}\left(t\right)}\frac{\mathrm{dP}\left(t\right)}{\mathrm{dt}}$ , wherein:

, be the derivative of the interior pulse pressure value of a cardiac cycle to the time;

Step (5) is calculated as follows the blood flow ordinary differential equation, obtains blood flow Q (t):

$\frac{\mathrm{dQ}\left(t\right)}{\mathrm{dt}}+\mathrm{\λ}\left(t\right)Q\left(t\right)+\mathrm{\ϵ}\left(t\right)Q^2\left(t\right)=A\left(t\right)F\left(t\right)$ , wherein:

$\mathrm{\λ}\left(t\right)=\frac{8\mathrm{\α\γ}}{R^2\left(t\right)}-\frac{\mathrm{\β}[{\mathrm{\β}}_1^2{\left({\mathrm{\α}}_2{\mathrm{\α}}_{2m}\right)}^4-1]}{P_m-[{\mathrm{\β}}_1^2{\left({\mathrm{\α}}_2{\mathrm{\β}}_{2m}\right)}^4-1](P\left(t\right)-P_m)}\frac{\mathrm{dP}\left(t\right)}{\mathrm{dt}}$ ，

, wherein:

λ (t) is the time dependent first power coefficient of blood flow,

, be the semi-cone angle under the naturalness in the blood vessel,

,

β ₁, be the ratio of the length when body length and the blood vessel stress balance of blood vessel under the normal physiological conditions, β ₁=0.57,

β _2m, vessel radius vessel radius ratio when not being out of shape during for the blood vessel stress balance, β _2m=1.105,

α ₂, be human body correction factor under the physiological condition,

,

γ is blood motion viscosity, is the ratio of hemodynamics viscosity and density of blood, and when first calculating, dynamic viscosity is progressive viscosity, mu ₀, under physiological condition, the desired viscosity when shear rate is enough large, μ ₀=3.6522081, be setting value,

α, β, β ₀, be human body correction factor under the physiological condition, wherein:

α, the non-linear pulse wave propagate coefficient when relating to the interior blood flow temporal evolution of a cardiac cycle during for calculating λ (t),

β, β ₀Non-linear pulse wave propagate coefficient when relating to the interior vessel radius temporal evolution of a cardiac cycle during for calculating ε (t),

α=0.57, β=1.173, β ₀=1.5 β are relevant with the vessel radius elastic dilatation, β ₀It is relevant with the vessel radius elastic shrinkage,

ε (t) is the time dependent secondary power coefficient of blood flow,

A (t) is the time dependent function of vessel radius in the cardiac cycle, describes the dynamic process of blood vessel elasticity shrinkage expansion in the cardiac cycle, represents with following formula:

A (t)=π R ²(t), wherein:

R (t) reacts vessel radius over time,

$R\left(t\right)=R_m*\sqrt{1+b*\ln \frac{P\left(t\right)}{P_m}}$ , wherein:

R _m, the vessel radius value does not represent with following formula when being out of shape,

$R_m=0.042+0.000625(1+0.36\sqrt{G})\sqrt{\mathrm{hl}}$ , wherein:

G is testee's age, age round number in this method,

H is testee's height, and unit is cm in this method,

L is testee's shoulder breadth, and unit is cm in this method;

Step (6) calculates the mean blood flow Q in the cardiac cycle _m, waveform parameter ξ, shear rate D (t) and flow factor S:

$Q_m=\frac{1}{T}{\∫}_0^{t}Q\left(t\right)\mathrm{dt}=\frac{1}{N}{\mathrm{\Σ}}_{n=1}^{n}Q\left(n\right)$ , wherein:

Q (t), the blood flow that measures during for moment t, unit are ml/sec,

Step (9) is calculated as follows heart stroke volume SV and exports value of calculation:

V=Q _m* T, wherein:

T is the time of a cardiac cycle;

SV represents the blood flow that heart is whenever beaten and once exported, and unit is ml/beat;

Step (10) is calculated as follows minute output of heart CO and exports value of calculation:

CO=Q _m*60/1000

CO represents the blood flow of heartbeat output in one minute, and unit is L/min;

Step (11) is calculated as follows power and exports value of calculation:

W _p=P _m*Q _m

W _p, being ventricle potential mean power at any one time, unit is W, is definite value in a cardiac cycle;

For cardiac function, the heart stroke volume is about 60-90ml/beat, and minute output is about 4-8L/min, and mean power is about 0.9-1.8W.According to stroke volume, minute output and mean power, obtain 3 groups of heart age values respectively, 3 values that then will obtain are carried out arithmetic average, obtain final heart age result.

Stroke volume SV〉90 ml/beat, judge that heart age is 15-20 year; 85＜SV＜90 ml/beat judges that heart age is 20-30 year; 80＜SV＜85 ml/beat judges that heart age is 30-40 year; 75＜SV＜80 ml/beat judges that heart age is 40-50 year; 70＜SV＜75 ml/beat judges that heart age is 50-60 year; 65＜SV＜70 ml/beat judges that heart age is 60-70 year; 60＜SV＜65 ml/beat judges that heart age is 70-80 year; SV＜60 ml/beat judges that heart age is more than 90 years old.

Minute output CO〉8L/min, judge that heart age is 15-20 year; 7.33＜CO＜8L/min judges that heart age is 20-30 year; 6.66＜CO＜7.33L/min judges that heart age is 30-40 year; 5.99＜CO＜6.66L/min judges that heart age is 40-50 year; 5.32＜CO＜5.99L/min judges that heart age is 50-60 year; 4.67＜CO＜5.32L/min judges that heart age is 60-70 year; 4＜CO＜4.67L/min judges that heart age is 70-80 year; CO＜4L/min judges that heart age is more than 90 years old.

Mean power W _P1.8W, judge that heart age is 15-20 year; 1.65＜W _P＜1.8W judges that heart age is 20-30 year; 1. 5＜W _P＜1.65W judges that heart age is 30-40 year; 1.35＜W _P＜1.5W judges that heart age is 40-50 year; 1.2＜W _P＜1.05W judges that heart age is 50-60 year; 1.05＜W _P＜1.2W judges that heart age is 60-70 year; 0.9＜W _P＜1.05W judges that heart age is 70-80 year; W _P＜0.9W judges that heart age is more than 90 years old.

Show according to general clinical value, more than three parameters have certain dependency, namely under the generic condition, 3 values that record can be in same heart age scope.Do not get rid of the not situation in same heart age scope, but the difference of 3 values can be very not large, still can obtain the heart age value by arithmetic average.

Description of drawings:

Fig. 1 is the population structure sketch map of heart age detector; Arrow represents to connect cable among the figure.

Fig. 2 utilizes non-linear pulse wave theory to calculate the flow chart of cardiac function and heart age.

The specific embodiment:

The present invention smoothly is wrapped in the upper arm place with upper arm sphygmomanometer and fixes, and the wave sensor that will refer to affectionately fight is fixed on homonymy finger middle finger place.Tested employing sitting posture gets final product, and upper arm and finger keep in the same horizontal line.Under quiescent condition, gather blood pressure and pulse ripple signal.Then, the tested test instrunment that removes is done 1 minute squat action, and number of times is no less than 50.Follow according to the method described above fixing test instrument, again acquisition and recording blood pressure and pulse ripple signal.After finishing, can operate the heart age signal gathering analysis meter, input the information such as tested age, height, body weight, shoulder breadth, instrument automatic analysis cardiac function, reckoning heart age provide result's (can show output, also can connect printer output).

Important technological parameters:

The used heating of this device unit be connected cable and can adopt the heater such as now widely used microwave heat therapeutic, ultrasound thermal therapy, radio-frequency (RF) thermotherapeutic, laserthermia and connect accordingly cable; This installs used conduit can adopt the present clinical Wicresoft's interventional therapy microtubular that has been widely used in; This installs used heating antenna and can improve a little on the basis of existing antenna for thermotherapy and get final product.Therefore, this device has good technological feasibility and clinical practice basis technically.Important technological parameters is as follows:

Supply voltage: 220/380 supply frequency: 50Hz

Conductor length:＞100mm printer interface: USB2.0.

Claims (2)

1. heart age detector, its characteristic is: comprise finger affectionately fight wave sensor, sphygomanometer, the wave sensor that refers to affectionately fight connects by being connected cable with signal gathering analysis meter with sphygomanometer.

2. application rights requires the method for 1 described heart age detector, it is characterized in that step following steps (1), and detector initializes:

Input: the radial artery pressure pulse wave that cuff sphygmomanometer detects and measured's pulse pressure value; The pulse waveform that the wave sensor that refers to affectionately to fight gathers;

Measured's personal information comprises sex, age, height, body weight and shoulder breadth at least;

Step (2), detector calculate the average pulse pressure P of measured in the cardiac cycle automatically _m, i.e. the arithmetic mean of instantaneous value of one group of detected value;

Step (3) calculates pressure pulse velocity of wave propagation c (t):

$C\left(t\right)=\sqrt{\frac{P\left(t\right)(1+b*\ln \frac{P\left(t\right)}{\mathrm{Pm}})}{\mathrm{\ρb}}}$ , wherein:

P (t) is pulse pressure value corresponding to each sampled point in the cardiac cycle,

ρ is density of blood, ρ=1.05*10 ^-3Kg/cm ³,

B is the pulse waveform parameter, is calculated as follows:

$b=\frac{a{\left(\frac{0.81}{k+0.48}\right)}^2-1}{\mathrm{\π}}$ ，

$k=\frac{P_m-P_d}{P_s-P_d}$ , wherein:

K is shape factor,

P _d, be diastolic pressure, P _sBe systolic pressure, be measured value,

A is the parameter relevant with velocity profile, a=0.57;

Step (4) calculates barometric gradient F (t):

$F\left(t\right)=\frac{1}{\mathrm{\ρc}\left(t\right)}\frac{\mathrm{dP}\left(t\right)}{\mathrm{dt}}$ , wherein:

, be the derivative of the interior pulse pressure value of a cardiac cycle to the time;

Step (5) is calculated as follows the blood flow ordinary differential equation, obtains blood flow Q (t):

$\frac{\mathrm{dQ}\left(t\right)}{\mathrm{dt}}+\mathrm{\λ}\left(t\right)Q\left(t\right)+\mathrm{\ϵ}\left(t\right)Q^2\left(t\right)=A\left(t\right)F\left(t\right)$ , wherein:

$\mathrm{\λ}\left(t\right)=\frac{8\mathrm{\α\γ}}{R^2\left(t\right)}-\frac{\mathrm{\β}[{\mathrm{\β}}_1^2{\left({\mathrm{\α}}_2{\mathrm{\α}}_{2m}\right)}^4-1]}{P_m-[{\mathrm{\β}}_1^2{\left({\mathrm{\α}}_2{\mathrm{\β}}_{2m}\right)}^4-1](P\left(t\right)-P_m)}\frac{\mathrm{dP}\left(t\right)}{\mathrm{dt}}$ ，

, wherein:

λ (t) is the time dependent first power coefficient of blood flow,

, be the semi-cone angle under the naturalness in the blood vessel,

,

β ₁, be the ratio of the length when body length and the blood vessel stress balance of blood vessel under the normal physiological conditions, β ₁=0.57,

β _2m, vessel radius vessel radius ratio when not being out of shape during for the blood vessel stress balance, β _2m=1.105,

α ₂, be human body correction factor under the physiological condition,

,

γ is blood motion viscosity, is the ratio of hemodynamics viscosity and density of blood, and when first calculating, dynamic viscosity is progressive viscosity, mu ₀, under physiological condition, the desired viscosity when shear rate is enough large, μ ₀=3.6522081, be setting value,

α, β, β ₀, be human body correction factor under the physiological condition, wherein:

α, the non-linear pulse wave propagate coefficient when relating to the interior blood flow temporal evolution of a cardiac cycle during for calculating λ (t),

β, β ₀Non-linear pulse wave propagate coefficient when relating to the interior vessel radius temporal evolution of a cardiac cycle during for calculating ε (t),

α=0.57, β=1.173, β ₀=-1.5, β is relevant with the vessel radius elastic dilatation, β ₀It is relevant with the vessel radius elastic shrinkage,

ε (t) is the time dependent secondary power coefficient of blood flow,

A (t) is the time dependent function of vessel radius in the cardiac cycle, describes the dynamic process of blood vessel elasticity shrinkage expansion in the cardiac cycle, represents with following formula:

A (t)=π R ²(t), wherein:

R (t) reacts vessel radius over time,

$R\left(t\right)=R_m*\sqrt{1+b*\ln \frac{P\left(t\right)}{P_m}}$ , wherein:

R _m, the vessel radius value does not represent with following formula when being out of shape,

$R_m=0.042+0.000625(1+0.36\sqrt{G})\sqrt{\mathrm{hl}}$ , wherein:

G is testee's age, age round number in this method,

H is testee's height, and unit is cm in this method,

L is testee's shoulder breadth, and unit is cm in this method;

Step (6) calculates the mean blood flow Q in the cardiac cycle _m, waveform parameter ξ, shear rate D (t) and flow factor S:

$Q_m=\frac{1}{T}{\∫}_0^{t}Q\left(t\right)\mathrm{dt}=\frac{1}{N}{\mathrm{\Σ}}_{n=1}^{n}Q\left(n\right)$ , wherein:

Q (t), the blood flow that measures during for moment t, unit are ml/sec,

Step (9) is calculated as follows heart stroke volume SV and exports value of calculation:

SV=Q _m* T, wherein:

T is the time of a cardiac cycle;

SV represents the blood flow that heart is whenever beaten and once exported, and unit is ml/beat;

Step (10) is calculated as follows minute output of heart CO and exports value of calculation:

CO=Q _m*60/1000

CO represents the blood flow of heartbeat output in one minute, and unit is L/min;

Step (11) is calculated as follows power and exports value of calculation:

W _p=P _m*Q _m

W _p, being ventricle potential mean power at any one time, unit is W, is definite value in a cardiac cycle;

For cardiac function, the heart stroke volume is about 60-90ml/beat, and minute output is about 4-8L/min, and mean power is about 0.9-1.8W; According to stroke volume, minute output and mean power, obtain 3 groups of heart age values respectively, 3 values that then will obtain are carried out arithmetic average, obtain final heart age result;

Stroke volume SV〉90 ml/beat, judge that heart age is 15-20 year; 85＜SV＜90 ml/beat judges that heart age is 20-30 year; 80＜SV＜85 ml/beat judges that heart age is 30-40 year; 75＜SV＜80 ml/beat judges that heart age is 40-50 year; 70＜SV＜75 ml/beat judges that heart age is 50-60 year; 65＜SV＜70 ml/beat judges that heart age is 60-70 year; 60＜SV＜65 ml/beat judges that heart age is 70-80 year; SV＜60 ml/beat judges that heart age is more than 90 years old;

Minute output CO〉8L/min, judge that heart age is 15-20 year; 7.33＜CO＜8L/min judges that heart age is 20-30 year; 6.66＜CO＜7.33L/min judges that heart age is 30-40 year; 5.99＜CO＜6.66L/min judges that heart age is 40-50 year; 5.32＜CO＜5.99L/min judges that heart age is 50-60 year; 4.67＜CO＜5.32L/min judges that heart age is 60-70 year; 4＜CO＜4.67L/min judges that heart age is 70-80 year; CO＜4L/min judges that heart age is more than 90 years old;

Mean power W _P1.8W, judge that heart age is 15-20 year; 1.65＜W _P＜1.8W judges that heart age is 20-30 year; 1. 5＜W _P＜1.65W judges that heart age is 30-40 year; 1.35＜W _P＜1.5W judges that heart age is 40-50 year; 1.2＜W _P＜1.05W judges that heart age is 50-60 year; 1.05＜W _P＜1.2W judges that heart age is 60-70 year; 0.9＜W _P＜1.05W judges that heart age is 70-80 year; W _P＜0.9W judges that heart age is more than 90 years old.

Images