CN1055829C - Brain willis circulation dynamics analysing method and instrument - Google Patents
Brain willis circulation dynamics analysing method and instrument Download PDFInfo
- Publication number
- CN1055829C CN1055829C CN 96116339 CN96116339A CN1055829C CN 1055829 C CN1055829 C CN 1055829C CN 96116339 CN96116339 CN 96116339 CN 96116339 A CN96116339 A CN 96116339A CN 1055829 C CN1055829 C CN 1055829C
- Authority
- CN
- China
- Prior art keywords
- centerdot
- artery
- right sides
- tau
- integral
- 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 Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
The present invention belongs to a cerebral circulation blood dynamics analyzing method and an instrument. A method for rapidly analyzing and calculating the dynamical parameter of cerebral circulation is established according to the anatomical model of a cerebral circulating arterial vascular bed, an equivalent network model thereof and a corresponding control equation, by using a Doppler to detect the blood flow speed of main branch blood vessels in the skull and by using an imaging method to obtain the geometric parameter of cerebral vessels. According to the analytical method, the analytical instrument is designed and is composed of a detection system, a collection analysis and calculation system and a storage and output system. The present invention has the advantages of great enhancement of a detection analysis speed and has important significance for the basic research and the clinical application of cerebral circulation.
Description
The invention belongs to field of medical technology, is a kind of cerebro-blood circulation dynamics analytical method and instrument.
Cerebrovascular disease is serious threat human life's commonly encountered diseases and a frequently-occurring disease.Cerebrovascular disease at first causes disturbance of cerebral circulation, makes some regional area cerebral tissue ischemia, and then causes brain cell death.
The cerebral circulation Arterial system, mainly by carotid artery, vertebra-basilar artery and basis cranii WILLIS ring etc. constitute blood supply net in the brain.Blood flows into this network from left and right sides internal carotid artery and left and right sides vertebral artery respectively.Analyze the interior blood motion rule of such multi-source network and the mechanical characteristic of each bifurcated artery blood vessel, crucial clinical and theory value is arranged preventing, diagnosing and treat cerebrovascular disease.
Because brain blood bed is enclosed in the skull deeply, therefore accurate the measurement and the dynamics of analyzing cerebral circulation, still technical theoretically all very difficult.
The eighties is ground into the TCD of merit mid-term abroad, and the using ultrasound Doppler technology can detect some endarterial blood flow rate of intracranial, but still powerless with analysis to the analysis, the particularly detection of cerebrovascular dynamic characteristic of whole cerebral circulation.
External late nineteen seventies, domestic late nineteen eighties has been invented method and the instrument of measuring the carotid artery system vessel properties from common carotid artery in succession.But this quasi-instrument can only be done quantitative description and analysis to the carotid artery system blood vessel, but can't to whole cerebral circulation particularly vertebra one basilar artery system and WILLIS ring etc. detect and analyze.
Therefore, development is used to detect and analyze the method and the instrument of whole cerebral circulation characteristic, becomes present clinical urgency problem to be solved.
The objective of the invention is to design a kind of method and instrument that can carry out the quick nondestructive check and analysis to whole cerebral circulation characteristic.
Because no matter cerebral circulation vascular bed is geometrical property or dynamics, all circulation has essential difference with body.Therefore, the kinetic model of extensive use once, analytical method etc. all can not be indiscriminately imitated in the research of cerebral circulation simply in the body circulating research.The present invention considers the particularity of brain Willis ring structure, in the kinetic parameter analytical method of having set up the suitable clinical practice of a cover on the theoretical model basis of experimental verification.
According to anatomical structure and data, the letter of cerebral circulation arteries can be become as Fig. 1.Shown anatomical model (Hillen, B et al, Analysis of flow andvascular resistance in a model of the circle of Willis, J.Biomech., 21,807-814,1988), wherein, each artery segment of cerebral circulation dissects title and the symbol contrast is as follows:
R
A, R
M, R
PBefore representing brain respectively, in, the vascular bed Peripheral resistance in artery downstream, 1 and 2 are left and right sides internal carotid artery inlet, 3 and 4 are left and right sides vertebral artery inlet.
Tremulous pulse title pipeline section symbol | Tremulous pulse title pipeline section symbol |
Internal carotid artery c 1,c 2Basilar artery b vertebral artery v 1,v 2Posterior communicating artery pc 1,pc 2Posterior cerebral artery I p 11,p 22 | Anterior cerebral artery I a 11,a 21Anterior communicating artery ac middle cerebral artery m 1,m 2Posterior cerebral artery II p 12,p 22Anterior cerebral artery II a 12,a 22 |
According to above-mentioned cerebrovascular anatomical model, the present invention has set up the cerebro-blood circulation dynamics lumped parameter model, and its equivalent circuit network model as shown in Figure 2.Wherein, we adopt one seven model of element simulation carotid artery subsystem, and in this subsystem, we have not only considered vascular resistance and compliance but also have considered the induction reactance (L of blood vessel
A1And L
A2).With a flow resistance R
AcDescribe the anterior communicating artery, the left and right sides carotid artery system is joined together.To vertebra one basilar artery system, respectively with the simulation of R-C two model of element, basilar artery are with a flow resistance R with every vertebral artery
bSimulation, its elasticity and induction reactance equivalence are in vertebral artery and posterior cerebral artery; The both sides posterior cerebral artery is adopted one 4 model of element simulation respectively; Between carotid artery system and vertebral-basilar artery system, adopt a R-L two model of elements simulations arteriae communicans posterior cerebri of connecting.Adopt above-mentioned aggregative model can simulate the cerebral circulation physiological conditions more truly, and aspect Mathematical treatment, be unlikely to undue complexity.
To cerebral circulation model shown in Figure 1, in the middle of concrete the application, also can suitably simplify according to actual needs, be constant even be zero as some parameters (R, C, L etc.) that can get in the network model.The implication of each the parameter representative in the network model specifically is listed as follows: parameter meaning P
1, P
2, P
3, P
4, four inputs (carotid artery and vertebral artery) inlet pressure C
C1, C
C2Left and right sides common carotid artery compliance R
M1, R
M2Left and right sides middle cerebral artery bed resistance R
A12, R
A22Left and right sides anterior cerebral artery bed resistance C
A1, C
A2Left and right sides anterior cerebral artery compliance R
A11, R
A21Left and right sides anterior cerebral artery section resistance L
A1, L
A2Left and right sides anterior cerebral artery induction reactance R
AcArteriae communicans anterior cerebri resistance R
C1, R
C2Left and right sides carotid artery resistance R
V1, R
V2Left and right sides vertebral artery resistance R
bBasilar artery resistance C
V1, C
V2Left and right sides vertebral artery compliance R
PC1, R
PC2Left and right sides arteriae communicans posterior cerebri resistance L
PC1, L
PC2Left and right sides arteriae communicans posterior cerebri induction reactance C
P1, C
P2Left and right sides posterior cerebral artery compliance L
P1, L
P2Left and right sides posterior cerebral artery induction reactance R
P12, R
P22Left and right sides posterior cerebral artery bed resistance R
P11, R
P21Left and right sides posterior cerebral artery resistance Q
C1, Q
C2Flow Qv1 in the carotid artery of the left and right sides, Qv2 left and right sides vertebra power arteries and veins flow Qm1, Qm2 left and right sides middle cerebral artery flow Q
A1, Q
A2Flow Q in the anterior cerebral artery
PC1, Q
PC2Flow Q in the posterior cerebral artery
Lp1, Q
Lp2Flow is concluded Pi representative pressure as can be known, Q to above-mentioned parameter in the arteriae communicans posterior cerebri
iRepresent flow, R
iRepresent resistance of blood flow (also claiming flow resistance), C
iExpression vascular compliance (also claiming fluid capacitance), L represents blood vessel induction reactance (also claiming influenza), subscript is represented the symbol of each corresponding vessel segment.
In the R-C-L analogue model that every blood vessel adopted
L
i=1.341
i/D
i 2 (1)
R
i=1.631
i/ D
i 4(2) wherein 1
iAnd D
iBe the length and the diameter of corresponding vessel segment, unit is cm.
According to this equivalent network model, the present invention has set up the governing equation (also claiming state equation) of following brain circulation system:
Here D=[d
Ij], E=[e
Ij], be 13 rank matrixes, X=[x
i], b=[b
i], be 13 rank vectors,
Expression is to the derivative of time, wherein X=(P
C1, P
C2, P
A1, P
A2, P
v, P
P1, P
P2, Q
A1, Q
A2, Q
Pc1, Q
Pc2, Q
Lp1, Q
Lp2)
T(4) b=(b
1, b
2, 0,0, b
5, 0,0,0,0,0,0,0,0)
T The matrix D element is d
11=C
C1, d
22=C
C2, d
33=C
A1, d
44=C
A1, d
55=C
v, d
66=C
P1, d
77=C
P2, d
88=L
A1, d
99=L
A1, d
1010=L
PC1, d
1111=L
PC2, d
1212=L
P1, d
1313=L
P2, d
Ij=0 (i ≠ j) (6) matrix E element is
, e
18=-1, e
110=-1,
, e
29=-1, e
211=-1,
, e
38=1,
, e
49=1,
, e
512=-1, e
513=-1,
, e
612=1,
, e
713=1, e
88=-R
A11, e
81=1, e
83=-1, e
99=-R
A21, e
92=1, e
94=-1,
, e
101=1,
, e
112=1,
e
1212=-R
b, e
125=1, e
126=-1, e
1213=-R
b, e
1213=-R
b, e
135=1, e
137=-1, e
1312=-R
b, all the other elements are 0.(7)
Use the result of equation (3), can calculate blood flow in other blood vessel by following formula:
Q
B=Q
Lp1+ Q
Lp2,
Q
P11=Q
P12+Q
PC1,Q
P21=Q
P22+Q
PC2, (8)
Governing equation (3) has provided the physical characteristic of network (as flow resistance R
i, fluid capacitance C
i, induction reactance L
i) wait with blood flowing characteristic (as pressure P
i, flow Q
i) between quantitative relationship.Can solve 13 unknown quantitys by equation (3).If known whole R
i, L
i, C
i, and know the pressure P of four arrival ends
1-P
4, can accurately solve every endovascular flow Q so
iOtherwise, if known four end points pressure P
1-P
4And flow Q in some blood vessel that can measure
i(as carotid artery, vertebral artery, middle cerebral artery, the anterior communicating artery, or the like), just can use equation (3) and adopt approximating method to solve the physical characteristic of corresponding pipeline section (as resistance R
i, compliance C
iDeng), be called indirect problem.Concrete grammar is seen another patent application (No95111513.8) of inventor.But this method is in clinical practice, and a patient's of analysis situation takes longer, generally wants more than 2 hours.The present invention utilizes the ultransonic state-of-the-art technology of current intracranial on above-mentioned model based, a kind of new analysis calculation method is provided, and the development corresponding instrument.The present invention adopts transcranial doppler to detect blood flow rate in the main branch vessel of intracranial, adopt the iconography method to obtain the cerebrovascular geometric parameter, and, reflect that brain Willis encircles the kinetic parameter of blood circulation state thereby can calculate rapidly and accurately to the simplification that equation is made a match and managed.
According to above-mentioned cerebrovascular dynamic model and state equation, can obtain the cerebral hemodynamic parameter value quickly and accurately particularly by following step.
1. use physiological pressure transducer or ultrasonic blood vessel diameter detecting instrument, detect patient's carotid artery and vertebrarterial pressure pulsation, and use brachial artery pressure numerical value, its waveform is demarcated, thereby obtain P1 (t), P2 (t), P3 (t), waveform and the numerical value of P4 (t) P1 (t).
2. use B ultrasound probe and doppler ultrasound velocity of blood flow probe, directly detect the interior blood-flow waveform of carotid artery and vertebral artery and numerical value: Qc1 (t), Qc2 (t), Qv1 (t), Qv2 (t).
3. use TCD,transcranial Doppler doppler transducer (TCD) and detect the intracranial middle cerebral artery, flow velocity waveforms and numerical value in anterior cerebral artery and the posterior cerebral artery.By patient DSA or/and the MRA blood vessel image is learned data, obtain the proportionate relationship of outer tremulous pulse (as common carotid artery, the vertebral artery etc.) diameter of entocranial artery diameter and cranium, the numerical value of the outer artery diameter of the cranium that application Type B ultrasonic probe or blood vessels caliber ultrasonic probe record, can obtain the caliber of intracranial vessel, further calculate endarterial blood-flow waveform and numerical value, i.e. Qa1 (t), Qa2 (t), Qm1 (t), Qm2 (t), Qlp1 (t) and Qlp2 (t) before, during and after the brain with formula group (8).
4. use DSA and MRA data, obtain the length and the diameter of relevant vessel segment, using formula (1), (2) calculate the resistance and the induction reactance of some single vessel of intracranial.Be the numerical value of Ra11, Ra12, Rpc1, Rpc2, Rac, Rb, La1, La2, Lpc1, Lpc2, Lp1, Lp2.
5. calculate carotid artery and vertebral artery resistance Rc1, Rc2, Rv1, Rv2.
(9)
Wherein, i=1,2, k desirable 10,11 or 12, and Zc1, Zc2 and Zv1, Zv2 are respectively carotid artery and vertebrarterial input impedance mould.They can be analyzed by the Fouriev of carotid artery and vertebral artery pressure and flow waves and obtain, promptly
i=1,2...,4 (10)
I=1,2 (11) p here
Ij, φ
Pij, Q
Cij, φ
QcijAnd Q
Vij, φ
QvijBe respectively mould and argument that pressure and flow Fourier change
6. adopt permanent theoretical method, can obtain the resistance of each blood vessel terminal of intracranial by following formula:
Here P and Q represent P (t) and Q (t) meansigma methods at a cardiac cycle, and P
i=Q
CiR
Ci+ (Q
Ai2+ (1)
I+1Qac) R
Ai1+ Q
Ai2R
Ai2(18) Q
Ac=(Q
A12R
A12-Q
A22R
A22)/R
Ac(19) Q
Pi2=Q
Lpi+ Q
Ci-Q
Mi-Q
Ai2+ (1)
I+1Q
Ac(20) f
i=P
i-Q
CiR
Ci-Q
Ai2R
Ai1(21) h
Ii=(1+R
Ai1/ R
Ac) Q
Ai2(22) h
12=-(R
A11/ R
Ac) Q
A22, h
21=-(R
A21/ R
Ac) Q
A12(23) i=1,2
7. use equation group (8) and can further obtain the interior pressure P of each pipeline section
C1=P
1-Q
C1R
C1., P
C2=P
2-Q
C2R
C2(24)
P
a1=Q
a12·R
a12, P
a2=Q
a22·R
a22 (26)
8. flow Q in other blood vessels of intracranial
A1, Q
A2, Q
Pc1And Q
Pc2For:
Wherein
i=1,2
Wherein
i=1,2
9. the intracranial compliance is:
Here, To and Ts are that heart contraction begins and the time that finishes.
According to above-mentioned analytical method, the present invention has also developed corresponding instrument, i.e. brain Willis ring circulation dynamics analytical tool.This instrument is divided into 3 parts: detection system, acquisition analysis system, storage output system.Wherein detection system has two parts, first is intracranial vessel blood flow rate test section, and it is a ultrasonic Doppler probe, and supersonic frequency is 2-3MHZ, adopt pulse to wear the cranium Doppler technology, can detect the blood flow rate in middle cerebral artery, prerolandic artery Rolando and artery and the basilar artery.This part also can adopt the TCD technology of present comparative maturity.Second portion is the outer vascular pressure flow detection part of cranium, comprises continuous ultrasound ripple doppler transducer, physiological pressure transducer or ultrasonic blood vessel diameter probe that a frequency is 5-8MHZ.Continuous ultrasound ripple probe is used to detect blood flow rate and the blood flow in cranium arteria carotis externa or the vertebral artery.Pressure transducer or ultraphonic pipe diameter probe are used to detect the blood vessel wall waveform of beating, and use brachial arterial pressure numerical value again it is demarcated, and then obtain pressure waveform numerical value.The collection computing system is reached by a master computer and A/D converter, communication microcomputer to be formed corresponding to the analytical calculation software of said method.By A/D converter aforementioned detected analog signal conversion is become digital signal, and adopt in the computer initial data as analytical calculation.If above-mentioned signal has converted digital signal (as data among the TCD) in collection, then can adopt the microcomputer communication mode, import data into master computer.Analytical calculation software is worked out according to the inventive method, and software block diagram is seen Fig. 5.It is according to above-mentioned blood flow and pressure information, and the cerebrovascular geometric parameter that obtains from imaging examination, and analysis meter is calculated the every hematodinamics index of brain Willis ring circulation.Data are stored with output and are made up of hard disk, floppy disk, tape, printer etc., the normal parameter value contrast of storing in patient's brain recirculation dynamic index that it can obtain the analytical calculation system and the main frame is judged, output analyzing and testing result, and prompting patient's hematodinamics state, for the clinical analysis cerebrovascular function provides quantitative target and reference proposition.
The apparatus structure block diagram as shown in Figure 6 and Figure 7.The present invention can provide two kinds of forms.A kind of distributing, promptly the test section utilizes existing equipment, and as TCD, CVA etc., promptly detecting the intracranial vessel blood flow rate is a unit, and detecting outer vascular flow speed of cranium and pressure is a unit, adopts communication modes to connect with main frame.Another kind is an integrated form, wherein all probes all is installed on the main frame, forms an integral body.
Clinical practice is for example: certain male patient, the applying pressure pick off records the left carotid artery pressure waveform, as shown in Figure 8, with ultrasonic technique detect flow waveform, as shown in Figure 9.The vertebral artery flow waveform, as shown in phantom in Figure 10.Use the inventive method, on Compaq486 DX-100 computer, moved about 30 seconds, promptly calculate cerebral hemodynamic parameter, list in the following table: the cerebral circulation power numerical value resistance (dyns/cm that Theoretical Calculation obtains
5) induction reactance (dyns
2/ cm
5) compliance (cm
5/ dyn) Rm
12.08 * 10
4La1 1.04 * 10
2Cc
17.53 * 10
-5Ra
124.11 * 10
4Lpc1 3.16 * 10
2Ca
18.26 * 10
-6Rp
122.22 * 10
4Lp1 2.47 * 10
2Cp
13.89 * 10
-6Rc
11.84 * 10
3Cv1 5.23 * 10
-5Rv
12.37 * 10
3These value inverse iterations are gone into the aforementioned calculation formula, calculate theoretical carotid artery and vertebral artery flow waveform and see solid line among Fig. 9 and Figure 10, contrast finds that Theoretical Calculation result is very approaching with the actual measurement value.This proves absolutely that the method that the present invention proposes is successful.It is feasible that clinical practice is used, and analytical calculation speed improves greatly.Fig. 1 is a cerebral circulation vascular bed anatomical model sketch map; Fig. 2 is the equivalent circuit network model corresponding to Fig. 1 anatomical model; Fig. 3 is a lumped parameter model of element sketch map; Fig. 4 is data flow block diagram among the present invention; Fig. 5 is analytical calculation software block diagram among the present invention; Fig. 6 is an instrument distributing block diagram of the present invention; Fig. 7 is the centralized block diagram of instrument of the present invention; Fig. 8 is detected carotid artery pressure waveform; Fig. 9 is carotid artery stream waveform, and wherein dotted line is a detection waveform, and solid line is the Theoretical Calculation waveform; Figure 10 is the vertebral artery flow waveform, and wherein dotted line is a detection waveform, and solid line is the Theoretical Calculation waveform;
Claims (2)
1 one kinds of brain willis circulation dynamics analysing methods according to the anatomical model of cerebrovascular bed, are set up equivalent circuit network model and control corresponding equation, it is characterized in that calculating the acquisition cerebral hemodynamic parameter by following step:
(1) uses physiological pressure transducer or ultrasonic blood vessel diameter detecting instrument, detect patient's carotid artery and vertebrarterial pressure pulsation, and use brachial artery pressure numerical value, its waveform is demarcated, thereby obtain P1 (t), P2 (t), P3 (t), waveform and the numerical value of P4 (t);
(2) use super ripple probe of Type B and doppler ultrasound velocity of blood flow probe, directly detect the interior blood-flow waveform of carotid artery and vertebral artery and numerical value: Qc1 (t), Qc2 (t), Qv1 (t), Qv2 (t);
(3) use TCD,transcranial Doppler doppler transducer (TCD) and detect the intracranial middle cerebral artery, flow velocity waveforms in anterior cerebral artery and the posterior cerebral artery and numerical value pass through patient DSA or/and the MRA blood vessel image is learned data, obtain the outer tremulous pulse of entocranial artery diameter and cranium (as common carotid artery, vertebral artery etc.) proportionate relationship of diameter, the numerical value of the outer artery diameter of the cranium that application Type B ultrasonic probe or blood vessels caliber ultrasonic probe record, can obtain the caliber of intracranial vessel, and before further calculating brain, in, blood-flow waveform and numerical value in the artery, i.e. Qa1 (t), Qa2 (t), Qm1 (t), Qm2 (t), Qlp1 (t) and Qlp2 (t);
(4) use DSA and MRA data, obtain the length and the diameter of relevant vessel segment, using formula
L
i=1.341
i/D
i 2 (1)
R
i=1.631
i/ D
i 4(2) resistance and the induction reactance that calculates some single vessel of intracranial is the numerical value of Ral1, Ral2, Rapc1, Rpc2, Rac, Rb, La1, La2, Lpc1, Lpc2, Lp1, Lp2;
(5) calculate carotid artery and vertebral artery resistance Rc1, Rc2, Rv1, Rv2:
(9)
Wherein, i=1,2, k desirable 10,11 or 12, and Zc1, Zc2 and Zv1, Zv2 are respectively carotid artery and vertebrarterial input impedance mould, and they can be analyzed by the Fouriev of carotid artery and vertebral artery pressure and flow waves and obtain, promptly
i=1,2...,4 (10)
I=1,2 (11) p here
Ij, φ
Pij, Q
Cij, φ
QcijAnd Q
Vij, φ
QvijBe respectively the mould argument that pressure and flow Fourier change
(6) adopt permanent theoretical method, can obtain the resistance of each blood vessel terminal of intracranial by following formula:
Here P and Q represent P (t) and Q (t) meansigma methods at a cardiac cycle, and P
i=Q
CiR
Ci+ (Q
Ai2+ (1)
I+1Q
Ac) R
Ai1+ Q
Ai2R
Ai2(18) Q
Ac=(Q
A12R
A12-Q
A22R
A22)/R
Ac(19) Q
Pi2=Q
Lpi+ Q
Ci-Q
Mi-Q
Ai2+ (1)
I+1Q
Ac(20) f
i=P
i-Q
CiR
Ci-Q
Ai2R
Ai1(21) h
Ii=(1+R
Ai1/ R
Ac) Q
Ai2(22) h
12=-(R
A11/ R
Ac) Q
A22, h
21=-(R
A21/ R
Ac) Q
A12(23) i=1,2 (7) use equation (8) can further obtain the interior pressure of pipeline section: P fully
C1=P
1-Q
C1R
C1, P
C2=P
2-Q
C2R
C2(24)
P
a1=Q
a12·R
a12, P
a2=Q
a22·R
a22 (26)
(8) flow Q in other blood vessels of intracranial
A1, Q
A2, Q
Pc1And Q
Pc2For:
i=1,2
I=1,2 (9) intracranial compliances are:
Here, To and Ts are that heart contraction begins and the time that finishes, and have the implication of related parameter to be listed as follows: parameter meaning P
1, P
2, P
3, P
4Four inputs (carotid artery and vertebral artery) inlet pressure C
C1, C
C2Left and right sides common carotid artery compliance R
M1, R
M2Left and right sides middle cerebral artery bed resistance R
A12, R
A22Left and right sides anterior cerebral artery bed resistance C
A1, C
A2Left and right sides anterior cerebral artery compliance R
A11, R
A21Left and right sides anterior cerebral artery section resistance L
A1, L
A2Left and right sides anterior cerebral artery induction reactance R
AcArteriae communicans anterior cerebri resistance R
C1, R
C2Left and right sides carotid artery resistance R
V1, R
V2Left and right sides vertebral artery resistance R
bBasilar artery resistance C
V1, C
V2Left and right sides vertebral artery compliance R
PC1, R
PC2Left and right sides arteriae communicans posterior cerebri resistance L
PC1, L
PC2Left and right sides arteriae communicans posterior cerebri induction reactance C
P1, C
P2Left and right sides posterior cerebral artery compliance L
P1, L
P2Left and right sides posterior cerebral artery induction reactance R
P12, R
P22Left and right sides posterior cerebral artery bed resistance R
P11, R
P21Left and right sides posterior cerebral artery resistance Q
C1, Q
C2Flow Qv1 in the carotid artery of the left and right sides, Qv2 left and right sides vertebra power arteries and veins flow Qm1, Qm2 left and right sides middle cerebral artery flow Q
A1, Q
A2Flow Q in the anterior cerebral artery
PC1, Q
PC2Flow Q in the posterior cerebral artery
Lp1, Q
Lp2Flow in the arteriae communicans posterior cerebri
2. a brain Willis encircles the circulation power analytical tool, by detection system, acquisition analysis system and storage output system are formed, wherein detection system comprises that a frequency that is used to detect intracranial vessel blood speed is the impulse ultrasound doppler transducer of 2-3MHZ, a frequency that is used to detect outer vascular flow speed of cranium and blood flow is the continuous ultrasound doppler transducer of 5-8MHZ, a physiological pressure transducer or a ultraphonic pipe diameter probe that is used to detect the blood vessel wall waveform, the storage output system is by being used to store the brain recirculation dynamic achievement data that cerebral circulation normal parameter numerical value and analytical calculation system obtain, and the hard disk of output comparing result, floppy disk, printer is formed, the collection analysis computing system comprises master computer, the analog signal conversion that the test section is obtained becomes the A/D converter of digital signal, it is characterized in that the analytical calculation system also comprise the digital signal that will measure import into the communication microcomputer of master computer and method according to claim 1 establishment be used for analytical calculation relatively the software of cerebral hemodynamic parameter form.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 96116339 CN1055829C (en) | 1996-04-25 | 1996-04-25 | Brain willis circulation dynamics analysing method and instrument |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 96116339 CN1055829C (en) | 1996-04-25 | 1996-04-25 | Brain willis circulation dynamics analysing method and instrument |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1146319A CN1146319A (en) | 1997-04-02 |
CN1055829C true CN1055829C (en) | 2000-08-30 |
Family
ID=5123452
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 96116339 Expired - Fee Related CN1055829C (en) | 1996-04-25 | 1996-04-25 | Brain willis circulation dynamics analysing method and instrument |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1055829C (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7650181B2 (en) | 2005-09-14 | 2010-01-19 | Zoll Medical Corporation | Synchronization of repetitive therapeutic interventions |
CN101991436A (en) * | 2010-10-30 | 2011-03-30 | 徐州市宝兴医疗设备有限公司 | Central arterial pressure and AI (Artificial Intelligence) measuring device |
US20140323928A1 (en) | 2013-04-30 | 2014-10-30 | Zoll Medical Corporation | Compression Depth Monitor with Variable Release Velocity Feedback |
CN112641467B (en) * | 2020-12-21 | 2022-03-01 | 西安交通大学 | Non-invasive transcranial ultrasonic blood pressure measuring device and method based on SonoVue |
-
1996
- 1996-04-25 CN CN 96116339 patent/CN1055829C/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN1146319A (en) | 1997-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Caballero et al. | A review on computational fluid dynamics modelling in human thoracic aorta | |
Yang et al. | Computation of flow pressure fields from magnetic resonance velocity mapping | |
Willemet et al. | Arterial pressure and flow wave analysis using time-domain 1-D hemodynamics | |
Cullen et al. | The use of system identification techniques in the analysis of oculomotor burst neuron spike train dynamics | |
CN101172042B (en) | Brain blood vessel circulation kinetic analysis method and apparatus | |
AU2010210731A1 (en) | Calculating cardiovascular parameters | |
Urchuk et al. | MR measurement of pulsatile pressure gradients | |
CN1055829C (en) | Brain willis circulation dynamics analysing method and instrument | |
AU2010210791A1 (en) | Detection of vascular conditions using arterial pressure waveform data | |
Bertaglia et al. | Computational hemodynamics in arteries with the one-dimensional augmented fluid-structure interaction system: viscoelastic parameters estimation and comparison with in-vivo data | |
CN1055828C (en) | Analysis method and instrument for dynamics of cerebral blood circulation | |
Liu et al. | Analysis of dynamic cerebral autoregulation using an ARX model based on arterial blood pressure and middle cerebral artery velocity simulation | |
Creswell et al. | An experimental method for evaluating constitutive models of myocardium in in vivo hearts | |
Keunen et al. | Nonlinear transcranial Doppler analysis demonstrates age-related changes of cerebral hemodynamics | |
Vakorin et al. | Inferring neural activity from BOLD signals through nonlinear optimization | |
Barschdorff et al. | Neural network based multi sensor heart sound analysis | |
Özbay et al. | Effects of window types on classification of carotid artery Doppler signals in the early phase of atherosclerosis using complex-valued artificial neural network | |
Swiercz et al. | Neural network technique for detecting emergency states in neurosurgical patients | |
Avanzolini et al. | Electrical analogs for monitoring vascular properties in artificial heart studies | |
Klebæk et al. | Neural network for sonogram gap filling | |
Zhang et al. | An in vitro circulatory device for studying blood flow electrical impedance in human common carotid arteries | |
Siewnicka et al. | Hybrid cardiovascular simulator as a tool for physical reproduction of the conditions prevailing in the apex of the heart | |
CN112541303B (en) | Method, device and storage medium for measuring residual blood percentage | |
JP2005502147A (en) | Device having an artificial neuron for describing the transmission characteristics of excitatory neurons, and a method for obtaining the transmission characteristics of excitatory neurons using artificial neurons | |
CN203290891U (en) | Noninvasive detecting system for intracranial blood circulation obstacle based on network topology analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1044688 Country of ref document: HK |