CN107822594A - 一种脑脊液分流流量检测系统及方法 - Google Patents
一种脑脊液分流流量检测系统及方法 Download PDFInfo
- Publication number
- CN107822594A CN107822594A CN201710896666.3A CN201710896666A CN107822594A CN 107822594 A CN107822594 A CN 107822594A CN 201710896666 A CN201710896666 A CN 201710896666A CN 107822594 A CN107822594 A CN 107822594A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- mfrac
- cerebrospinal fluid
- heat
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 210000001175 cerebrospinal fluid Anatomy 0.000 title claims abstract description 134
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000005284 excitation Effects 0.000 claims abstract description 78
- 238000001514 detection method Methods 0.000 claims abstract description 36
- 238000012545 processing Methods 0.000 claims abstract description 24
- 238000012546 transfer Methods 0.000 claims abstract description 16
- 230000008859 change Effects 0.000 claims abstract description 15
- 238000006243 chemical reaction Methods 0.000 claims abstract description 8
- 238000005057 refrigeration Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims description 42
- 239000000203 mixture Substances 0.000 claims description 18
- 238000013178 mathematical model Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 8
- 238000013517 stratification Methods 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 4
- 230000035939 shock Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 7
- 208000003906 hydrocephalus Diseases 0.000 description 5
- 210000004556 brain Anatomy 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000001816 cooling Methods 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 229920001296 polysiloxane Polymers 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000007920 subcutaneous administration Methods 0.000 description 3
- 229910001069 Ti alloy Inorganic materials 0.000 description 2
- 210000000683 abdominal cavity Anatomy 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000002500 effect on skin Effects 0.000 description 2
- 210000005081 epithelial layer Anatomy 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 210000001767 medulla oblongata Anatomy 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- RBTBFTRPCNLSDE-UHFFFAOYSA-N 3,7-bis(dimethylamino)phenothiazin-5-ium Chemical compound C1=CC(N(C)C)=CC2=[S+]C3=CC(N(C)C)=CC=C3N=C21 RBTBFTRPCNLSDE-UHFFFAOYSA-N 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000000576 arachnoid Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000004925 denaturation Methods 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 210000004379 membrane Anatomy 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 229960000907 methylthioninium chloride Drugs 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010412 perfusion Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000741 silica gel Substances 0.000 description 1
- 229910002027 silica gel Inorganic materials 0.000 description 1
- 210000000278 spinal cord Anatomy 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 210000002330 subarachnoid space Anatomy 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- External Artificial Organs (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明涉及一种脑脊液分流流量检测系统及方法。包括MCU模块及与该MCU模块连接的计算机、信号处理模块、热激发驱动模块,所述热激发驱动模块还经制冷模块连接有导热装置,所述信号处理模块还连接有信号采集模块;MCU模块通过热激发驱动模块控制制冷模块工作,通过导热装置传导热激发信号,从而改变热激发点脑脊液温度;而后,信号采集模块采集采样点脑脊液温度信号,信号处理模块进行信号预处理和A/D转换,传输给MCU模块,计算脑脊液分流流量,并上传计算数据至计算机。本发明首创应用热传导技术于脑脊液分流流量检测,可以模拟检测患者在任何条件下的脑脊液分流流量,比现有X‑CT和MRI设备效率高,省钱,无辐射,另外本申请仅需采用一个采样点即可实现脑脊液分流流量的检测。
Description
技术领域
本发明属于热力学和信息科学技术在生物医学上的应用,具体涉及一种脑脊液分流流量检测系统及方法。
背景技术
脑积水,即循环于蛛网膜(granulations)下腔包围着脑部与脊髓的脑脊液异常增加或累积,导致脑脊液压迫脑部,使穿越枕骨大孔的延髓受损的一种疾病。
脑脊液分流术是治疗脑积水效果明显的一种脑神经外科手术,如图1中植入人体内的管线部分所示。在脑室插入硅胶导管,用置入的一组单向阀把流通梗阻于蛛网膜的多余脑脊液排入腹腔,被再吸收。
手术能迅速改善脑积水引起的神经功能损害,但要发挥远期效果及作用,就要随时准确控制及调整脑脊液分流量。调整正确和及时,病情好转;调整错误或延误,则脑延髓受损严重直至危及生命。
控制和调整脑脊液分流量的前提是分流量检测。目前的X-CT检测过程是,先注入美兰等示踪剂到脑脊液循环系统,再核素脑池造影和CT术断层扫描,患者既痛苦又伤身。另一种检测方法是增强型核磁共振脑部检查,每次费用高达千元以上且检测结果不准确。如果不检测、少检测、凭经验,就无法把控治疗进展。据统计,脑脊液分流术目前的平均有效率约50%。分流管阻塞、引流不畅、分流过度,分流失败、分流位置不佳等,其发生率为40%左右,是手术失败的主要原因。
如表1所示,为现成的可检测脑脊液流量的几种技术设备,表中的技术及设备是应用层析扫描成像技术实现的,对使用环境和操作水平均有较高要求,且非常昂贵。
表1可检测脑脊液流速的技术及设备
现有相关技术如下:
1、现有的美国专利的装置,其有一个热激发和2个采样点,是一种用于确定植入式CSF分流器中的连续CSF流速的装置和方法。该系统/方法利用形成在柔性垫上的珀尔帖传感器,放置在患者皮肤上。珀尔帖传感器包括耦合到热电阻器的珀尔帖装置,其在CSF分流位置处与患者的皮肤接触。佩尔贴装置连续操作,由珀尔帖温度传感器控制到低于患者核心温度的预定温度,以形成温度差,使由皮肤温度传感器检测皮肤/CSF流产生的任何热量,并且珀尔帖温度传感器上游和下游温度传感器以及控制温度传感器被用于形成用于校准对应于实时CSF流速的珀尔帖信号的零流量基线。传感器处理装置处理用于生成零流量基线和珀尔帖信号的所有传感器数据[1]。
2、另一美国专利是一种用于定量测量植入皮下的分流管中CSF流量的系统。该系统包括聚集在三个部分中的热敏元件阵列,放置在皮肤表面上的冷却装置和相关联的数据采集和分析装置。将两个传感器部分放置在皮肤上的分流器上,并测量与CSF运动相关的实时温度响应。靠近冷却装置放置的一个阵列收集关于皮肤热性质的数据,包括皮肤热变性,比热,扩散性,灌注和热惯性。该方法包括评估皮肤的热性质并测量分流管中的CSF流量。该方法可用于分流术评估,脑脊液阀调节程序和与流量相关的CSF流量测量。当测量相对CSF流量时,仅需要使用一个传感器部分,而不需要检测皮肤热性质和连续地施加冷却装置[2]。
3、来自美国的基于热传导原理的脑脊液分流量定量检测方法是应用热飞行时间的计算方法,且需要2个采样点[3]。
4、美国相关知识产权装置的测评报告[4]。
5、国内和其他国家没见报道。
参考文献:
[1]huntcheck,Inc..CONTINUOUS REAL-TIME CSF FLOW MONITOR AND METHOD,US2015201882[P].2015-07-23
[2]Shuntcheck,Inc..Real time CSF flow measurement system and method,US8894584[P].2014-11-25
[3]S Rajasekaran,Hongwei Qu,K Zakalik(2015IEEE).Thermal measurementof cerebrospinal fluid flow rate in Hydrocephalus shunt.SENSORS.2015-ieeexplore.ieee.org
[4]U.S.Department of Neurosurgery,Children’s Hospital Boston,HarvardMedical School,Boston,Massachusetts;Department of Neurosurgery,NationwideChildren’s Hospital,Columbus,Ohio;Division of Neurosurgery,Children’sNational Medical Center,Washington,District of Columbia.Evaluation of theShuntCheck Noninvasive Thermal Technique for Shunt Flow Detection inHydrocephalic Patients.[R]NEUROSURGERY 68:198–205,2011。
发明内容
本发明的目的在于提供一种脑脊液分流流量检测系统及方法,采用热传导流量数学模型,可以模拟检测患者在任何条件下的脑脊液分流流量,比现有X-CT和MRI设备效率高,省钱,无辐射,另外本申请仅需采用一个采样点即可实现脑脊液分流流量的检测。
为实现上述技术方案,本发明提供了一种脑脊液分流流量检测系统,包括MCU模块及与该MCU模块连接的计算机、信号处理模块、热激发驱动模块,所述热激发驱动模块还经制冷模块连接有导热装置,所述信号处理模块还连接有信号采集模块;MCU模块通过热激发驱动模块控制制冷模块工作,通过导热装置传导热激发信号,从而改变热激发点脑脊液温度;而后,信号采集模块采集采样点脑脊液温度信号,信号处理模块进行信号预处理和A/D转换,传输给MCU模块,并经由MCU模块传输给计算机,MCU模块和计算机控制检测过程,计算脑脊液分流流量,并上传计算数据至计算机;所述MCU与计算机构成分布式系统。
在本发明一实施例中,所述热激发点、采样点各为一个。
在本发明一实施例中,所述导热装置产生的热激发信号是通过作用于埋设于人体皮肤下的分流管的脑脊液,具体为分流管的热激发点处的脑脊液,从而改变热激发点处的脑脊液温度;该热激发过程可以是直接的,也可以是间接的。
在本发明一实施例中,所述采用热传导流量数学模型计算脑脊液分流流量具体方式如下,
为叙述简便直接,可设热激发点、采样点皮肤层均匀且各项同性,忽略毛细血管和人体热源影响,建立热传导数学模型为齐次偏微分方程:
和
其中,a是皮肤比热ca、密度ka、热传导系数ρa组成的常数;b是脑脊液比热cb、密度kb、热传导系数pb组成的常数;d是导管壁比热cd、密度kd、热传导系数pd组成的常数;u(x,y,z,t)为皮肤解剖结构分层温度场中皮肤任一点G(x,y,z),在时刻t时的温度值;
忽略热能侧面流动部位影响,则方程(1)、(2)、(3)可简化为垂直于皮肤方向和沿分流管线方向的一维热传导方程:
和
通过方程(4)、(5)、(6),即可得皮肤热传导时间τa,分流管热传导时间τd,脑脊液温度变化热传导时间τb,脑脊液流动时间τ0,和热激发开始至采样结束的时间τc之间,有
τc=2τa+τb+2τd+τ0 (7)
已知分流管半径r,热激发与采样点间距离l,即可得出脑脊液流速v和流量m:
v=l/(τc-2τa-2τd-τb) (8)
m=πr2s (9)。
在本发明一实施例中,若考虑毛细血管和人体内部热源,及其热能侧面流动部位影响:
其中,F(x,y,z,t)为单位时间内人体内部单位体积中所产生的热量和外部热交换量之和;
则相应的热传导流量数学模型(1)-(6)应为非齐次热传导方程:
及
及
若考虑导热分层的影响,则以上热传导方程就是一个已知物体在边界上的温度状况或热交换情况,以及已知物体在初始时刻的温度,来求解物体在以后时刻的温度的问题。
本发明还提供了一种脑脊液分流流量检测方法,设置热激发点、采样点,在热激发点产生热激发信号,从而改变热激发点脑脊液温度,而后在采样点采样脑脊液温度信号,最后根据采样的脑脊液温度信号,通过采用热传导流量计算方法计算脑脊液分流流量。
在本发明一实施例中,所述热激发点、采样点各为一个。
在本发明一实施例中,所述热激发信号是通过作用于埋设于人体皮肤下的分流管的脑脊液,具体为分流管的热激发点处的脑脊液,从而改变热激发点处的脑脊液温度;该热激发过程可以是直接的,也可以是间接的。
在本发明一实施例中,所述采用热传导流量计算方法计算脑脊液分流流量的具体方式如下,
为叙述简便直接,可设热激发点、采样点皮肤层均匀且各项同性,忽略毛细血管和人体热源影响,建立热传导数学模型为齐次偏微分方程:
和
其中,a是皮肤比热ca、密度ka、热传导系数ρa组成的常数;b是脑脊液比热cb、密度kb、热传导系数pb组成的常数;d是导管壁比热cd、密度kd、热传导系数pd组成的常数;u(x,y,z,t)为皮肤解剖结构分层温度场中皮肤任一点G(x,y,z),在时刻t时的温度值;
忽略热能侧面流动部位影响,则方程(1)、(2)、(3)可简化为垂直于皮肤方向和沿分流管线方向的一维热传导方程:
和
通过方程(4)、(5)、(6),即可得皮肤热传导时间τa,分流管热传导时间τd,脑脊液温度变化热传导时间τb,脑脊液流动时间τ0,和热激发开始至采样结束的时间τc之间,有
τc=2τa+τb+2τd+τ0 (7)
已知分流管半径r,热激发与采样点间距离l,即可得出脑脊液流速v和流量m:
v=l/(τc-2τa-2τd-τb) (8)
m=πr2s (9)。
在本发明一实施例中,若考虑毛细血管和人体内部热源,及其热能侧面流动部位影响:
其中,F(x,y,z,t)为单位时间内人体内部单位体积中所产生的热量和外部热交换量之和;
则相应的热传导流量数学模型(1)-(6)应为非齐次热传导方程:
及
及
若考虑导热分层的影响,则以上热传导方程就是一个已知物体在边界上的温度状况或热交换情况,以及已知物体在初始时刻的温度,来求解物体在以后时刻的温度的问题。
相较于现有技术,本发明具有以下有益效果:本发明首创应用热传导技术于脑脊液分流流量检测,可以模拟检测患者在任何条件下的脑脊液分流流量,比现有X-CT和MRI设备效率高,省钱,无辐射,另外本申请仅需采用一个采样点即可实现脑脊液分流流量的检测。
附图说明
图1为脑积水分流术及分流量检测示意图。
图2为沿分流管方向的检测原理剖面图。
图3为本发明的脑脊液分流流量检测系统图。
具体实施方式
下面结合附图,对本发明的技术方案进行具体说明。
本发明的一种脑脊液分流流量检测系统,一种脑脊液分流流量检测系统,包括MCU模块及与该MCU模块连接的计算机、信号处理模块、热激发驱动模块,所述热激发驱动模块还经制冷模块连接有导热装置,所述信号处理模块还连接有信号采集模块;MCU模块通过热激发驱动模块控制制冷模块工作,通过导热装置传导热激发信号,从而改变热激发点脑脊液温度;而后,信号采集模块采集采样点脑脊液温度信号,信号处理模块进行信号预处理和A/D转换,传输给MCU模块,并经由MCU模块传输给计算机,MCU模块和计算机控制检测过程,计算脑脊液分流流量,并上传计算数据至计算机;所述MCU与计算机构成分布式系统。
本发明还提供了一种脑脊液分流流量检测方法,通过设置热激发点、采样点,在热激发点产生热激发信号,从而改变热激发点脑脊液温度,而后在采样点采样脑脊液温度信号,最后根据采样的脑脊液温度信号,通过采用热传导流量计算方法计算脑脊液分流流量。
在本发明中,所述热激发点、采样点各为一个。
在本发明中,所述热激发信号是通过作用于埋设于人体皮肤下的分流管(人体皮肤从上到下依次分为上皮层8、真皮层9、皮下层10,分流管1的导管壁11即位于皮下层10下)的脑脊液,具体为分流管的热激发点处的脑脊液,从而改变热激发点处的脑脊液温度;该热激发过程可以是直接的,也可以是间接的。
在本发明中,所述采用热传导流量计算方法计算脑脊液分流流量的具体方式如下,
为叙述简便直接,可设热激发点、采样点皮肤层均匀且各项同性,忽略毛细血管和人体热源影响,建立热传导数学模型为齐次偏微分方程:
和
其中,a是皮肤比热ca、密度ka、热传导系数ρa组成的常数;b是脑脊液比热cb、密度kb、热传导系数pb组成的常数;d是导管壁比热cd、密度kd、热传导系数pd组成的常数;u(x,y,z,t)为皮肤解剖结构分层温度场中皮肤任一点G(x,y,z),在时刻t时的温度值;
忽略热能侧面流动部位影响,则方程(1)、(2)、(3)可简化为垂直于皮肤方向和沿分流管线方向的一维热传导方程:
和
通过方程(4)、(5)、(6),即可得皮肤热传导时间τa,分流管热传导时间τd,脑脊液温度变化热传导时间τb,脑脊液流动时间τ0,和热激发开始至采样结束的时间τc之间,有
τc=2τa+τb+2τd+τ0 (7)
已知分流管半径r,热激发与采样点间距离l,即可得出脑脊液流速v和流量m:
v=l/(τc-2τa-2τd-τb) (8)
m=πr2s (9)
进一步的,若考虑毛细血管和人体内部热源,及其热能侧面流动部位影响:
其中,F(x,y,z,t)为单位时间内人体内部单位体积中所产生的热量和外部热交换量之和;
则相应的热传导流量数学模型(1)-(6)应为非齐次热传导方程:
及
及
若考虑导热分层的影响,则以上热传导方程就是一个已知物体在边界上的温度状况或热交换情况,以及已知物体在初始时刻的温度,来求解物体在以后时刻的温度的问题。这里,本发明上述所使用的人体的各种参数及其函数,系发明人数以万次实验数据取得的,为本发明主要保护点。
以下为本发明的具体实现过程。
图1为脑脊液分流流量检测实例图,图2为沿分流管方向的检测原理剖面图,其中,1为分流管,2为采样点,3为分流管腹腔端,4为热激励点,5为脑室,6为脑脊液(CSF)流向,7为腹腔,8为上皮层、9为真皮层、10为皮下层,11为导管壁,12为脑脊液。
本发明的一种脑脊液分流流量检测系统及方法,采用热传导流量数学模型,可以模拟检测患者在任何条件下的脑脊液分流流量,比现有X-CT和MRI设备效率高,省钱,无辐射,另外本发明仅需采用一个采样点即可实现脑脊液分流流量的检测。
本发明涉及一种脑脊液分流流量检测系统,包括MCU模块及与该MCU模块连接的计算机、信号处理模块、热激发驱动模块,所述热激发驱动模块还经制冷模块连接有导热装置,所述信号处理模块还连接有信号采集模块;MCU模块通过热激发驱动模块控制制冷模块工作,通过导热装置传导热激发信号,从而改变热激发点脑脊液温度,而后信号采集模块采集采样点脑脊液温度信号,并通信号处理模块进行A/D转换后传输给MCU模块,经由MCU模块采用热传导流量计算方法计算脑脊液分流流量,并上传至计算机进行存储处理。
所述技术及设备属于一种基于热传导理论的电子流量传感器,是热力学、流体力学、计算机科学在生物医学上的应用,是一种非侵入式的医疗器械。其用于脑脊液流速检测的传感器在图1、图2中的2、4两个部分。2、4可以是2段代替硅树脂分流管的钛合金管,也可以是依旧不变的硅胶管。使用钛合金管的优点是热阻率远低于硅胶管,易于加温和采样,传感器精度高。
传感器流量检测原理是,在皮肤上施加热激发信号如图2的4处。热激发信号将改变埋在皮肤下的脑脊液温度,分流管里的脑脊液流动会引起沿线皮肤温度变化,随后连续采集2处的温度变化。
利用2、4两处在热传导和脑脊液流作用下的相关时间,再依据流体力学、热力学的热传导偏微分方程(扩散方程),在给定初始条件和设定热能流经8、9、10、11等处边界条件后算出静态热扩散时间;这样,已知分流管内径,采样点距离,各物理点的热扩散时间和流体流动时间差,能算出脑脊液分流流速流量了。因此,脑脊液流量检测系统功能如图3所示。
图3中,闭环控制制冷模块产生热激发信号。信号采集模块将多路模拟信号送入有调理放大、AD转换功能的前端信号处理模块,MCU模块实现现场控制、通信、数据处理功能。
通信是指MCU与AD转换模块、MCU与计算机之间的通信。
MCU数据处理内容有,诊断设备是否正常,逻辑判断所采集的数据是否在某段正常取值范围,查表取值内存单元里经复杂计算后预存的温度补偿、线性插值、初始设定等参数值,转换不同类型数值等简单计算。
计算机从MCU单元获取数据序列,利用多种数学模型等计算脑脊液分流流量并生成数据库文件,计算机也能作为上位机,监控和设定MCU所需的参数值。
本发明具体采用了基于热传导的流量监测原理,具体实现如下:
对具体的热激发点和采样点(见图1,图2),先获取其皮肤层、毛细血管和人体热源的各项参数,作为已知的热传导方程的边界条件。为简便起见,这里设皮肤层均匀且各向同性,忽略毛细血管和人体热源影响,传感器关于热传导的数学模型为齐次偏微分方程:
和
其中,a是皮肤比热ca、密度ka、热传导系数ρa组成的常数;b是脑脊液比热cb、密度kb、热传导系数pb组成的常数;d是导管壁比热cd、密度kd、热传导系数pd组成的常数;u(x,y,z,t)为皮肤解剖结构分层温度场中皮肤任一点G(x,y,z),在时刻t时的温度值;
忽略热能侧面流动部位影响,则方程(1)、(2)、(3)可简化为垂直于皮肤方向和沿分流管线方向的一维热传导方程:
和
通过方程(4)、(5)、(6),即可得皮肤热传导时间τa,分流管热传导时间τd,脑脊液温度变化热传导时间τb,脑脊液流动时间τ0,和热激发开始至采样结束的时间τc之间,有
τc=2τa+τb+2τd+τ0 (7)
已知分流管半径r,热激发与采样点间距离l,即可得出脑脊液流速v和流量m:
v=l/(τc-2τa-2τd-τb) (8)
m=πr2s (9)
进一步的,若考虑毛细血管和人体内部热源,及其热能侧面流动部位影响:
其中,F(x,y,z,t)为单位时间内人体内部单位体积中所产生的热量和外部热交换量之和;
则相应的热传导流量数学模型(1)-(6)应为非齐次热传导方程:
及
及
若考虑导热分层的影响,即图1中所示的8,9,10,11,12的影响,则以上热传导方程就是一个已知物体在边界上的温度状况(或热交换情况)和物体在初始时刻的温度,来求解物体在以后时刻的温度的问题。因此这里的热传导方程最自然的一个定解问题就是在已给的初始条件和边界条件下求问题的解,是公知的。这里,本发明上述所使用的人体的各种参数及其函数,系发明人数以万次实验数据取得的,为本发明主要保护点。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其
他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (10)
1.一种脑脊液分流流量检测系统,其特征在于:包括MCU模块及与该MCU模块连接的计算机、信号处理模块、热激发驱动模块,所述热激发驱动模块还经制冷模块连接有导热装置,所述信号处理模块还连接有信号采集模块;MCU模块通过热激发驱动模块控制制冷模块工作,通过导热装置传导热激发信号,从而改变热激发点脑脊液温度;而后,信号采集模块采集采样点脑脊液温度信号,信号处理模块进行信号预处理和A/D转换,传输给MCU模块;所述MCU与计算机构成分布式系统,系统中MCU是下位机,计算机是上位机。
2.根据权利要求1所述的一种脑脊液分流流量检测系统,其特征在于:所述热激发点、采样点各为一个。
3.根据权利要求1所述的一种脑脊液分流流量检测系统,其特征在于:所述导热装置产生的热激发信号是通过作用于埋设于人体皮肤下的分流管的脑脊液,具体为分流管的热激发点处的脑脊液,从而改变热激发点处的脑脊液温度;该热激发过程可以是直接的,也可以是间接的。
4.根据权利要求1所述的一种脑脊液分流流量检测系统,其特征在于:所述采用热传导流量数学模型计算脑脊液分流流量具体方式如下:
为叙述简便直接,可设热激发点、采样点皮肤层均匀且各项同性,忽略毛细血管和人体热源影响,建立热传导数学模型为齐次偏微分方程:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
和
其中,a是皮肤比热ca、密度ka、热传导系数ρa组成的常数;b是脑脊液比热cb、密度kb、热传导系数pb组成的常数;d是导管壁比热cd、密度kd、热传导系数pd组成的常数;u(x,y,z,t)为皮肤解剖结构分层温度场中皮肤任一点G(x,y,z),在时刻t时的温度值;
忽略热能侧面流动部位影响,则方程(1)、(2)、(3)可简化为垂直于皮肤方向和沿分流管线方向的一维热传导方程:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
和
通过方程(4)、(5)、(6),即可得皮肤热传导时间τa,分流管热传导时间τd,脑脊液温度变化热传导时间τb,脑脊液流动时间τ0,和热激发开始至采样结束的时间τc之间,有
τc=2τa+τb+2τd+τ0 (7)
已知分流管半径r,热激发与采样点间距离l,即可得出脑脊液流速v和流量m:
v=l/(τc-2τa-2τd-τb) (8)
m=πr2s (9)。
5.根据权利要求4所述的一种脑脊液分流流量检测系统,其特征在于:若考虑毛细血管和人体内部热源,及其热能侧面流动部位影响:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>c</mi>
<mi>&rho;</mi>
</mrow>
</mfrac>
</mrow>
其中,F(x,y,z,t)为单位时间内人体内部单位体积中所产生的热量和外部热交换量之和;
则相应的热传导流量数学模型(1)-(6)应为非齐次热传导方程:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
及
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
及
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
若考虑导热分层的影响,则以上热传导方程就是一个已知物体在边界上的温度状况或热交换情况,以及已知物体在初始时刻的温度,来求解物体在以后时刻的温度的问题。
6.一种脑脊液分流流量检测方法,其特征在于:设置热激发点、采样点,在热激发点产生热激发信号,从而改变热激发点脑脊液温度,而后在采样点采样脑脊液温度信号,最后根据采样的脑脊液温度信号,通过采用热传导流量计算方法计算脑脊液分流流量。
7.根据权利要求6所述的一种脑脊液分流流量检测方法,其特征在于:所述热激发点、采样点各为一个。
8.根据权利要求6所述的一种脑脊液分流流量检测方法,其特征在于:所述热激发信号是通过作用于埋设于人体皮肤下的分流管的脑脊液,具体为分流管的热激发点处的脑脊液,从而改变热激发点处的脑脊液温度;该热激发过程可以是直接的,也可以是间接的。
9.根据权利要求6所述的一种脑脊液分流流量检测方法,其特征在于:所述采用热传导流量计算方法计算脑脊液分流流量的具体方式如下,
为叙述简便直接,可设热激发点、采样点皮肤层均匀且各项同性,忽略毛细血管和人体热源影响,建立热传导数学模型为齐次偏微分方程:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
和
其中,a是皮肤比热ca、密度ka、热传导系数ρa组成的常数;b是脑脊液比热cb、密度kb、热传导系数pb组成的常数;d是导管壁比热cd、密度kd、热传导系数pd组成的常数;u(x,y,z,t)为皮肤解剖结构分层温度场中皮肤任一点G(x,y,z),在时刻t时的温度值;
忽略热能侧面流动部位影响,则方程(1)、(2)、(3)可简化为垂直于皮肤方向和沿分流管线方向的一维热传导方程:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
和
通过方程(4)、(5)、(6),即可得皮肤热传导时间τa,分流管热传导时间τd,脑脊液温度变化热传导时间τb,脑脊液流动时间τ0,和热激发开始至采样结束的时间τc之间,有
τc=2τa+τb+2τd+τ0 (7)
已知分流管半径r,热激发与采样点间距离l,即可得出脑脊液流速v和流量m:
v=l/(τc-2τa-2τd-τb) (8)
m=πr2s (9)。
10.根据权利要求9所述的一种脑脊液分流流量检测系统,其特征在于:若考虑毛细血管和人体内部热源,及其热能侧面流动部位影响:
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>c</mi>
<mi>&rho;</mi>
</mrow>
</mfrac>
</mrow>
其中,F(x,y,z,t)为单位时间内人体内部单位体积中所产生的热量和外部热交换量之和;
则相应的热传导流量数学模型(1)-(6)应为非齐次热传导方程:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
及
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
及
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
若考虑导热分层的影响,则以上热传导方程就是一个已知物体在边界上的温度状况或热交换情况,以及已知物体在初始时刻的温度,来求解物体在以后时刻的温度的问题。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710896666.3A CN107822594B (zh) | 2017-09-28 | 2017-09-28 | 一种脑脊液分流流量检测系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710896666.3A CN107822594B (zh) | 2017-09-28 | 2017-09-28 | 一种脑脊液分流流量检测系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107822594A true CN107822594A (zh) | 2018-03-23 |
CN107822594B CN107822594B (zh) | 2020-10-13 |
Family
ID=61644050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710896666.3A Expired - Fee Related CN107822594B (zh) | 2017-09-28 | 2017-09-28 | 一种脑脊液分流流量检测系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107822594B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109979295A (zh) * | 2019-05-07 | 2019-07-05 | 莆田学院 | 一种脊髓型颈椎病的病变演示装置及演示方法 |
CN114366859A (zh) * | 2021-10-29 | 2022-04-19 | 首都医科大学宣武医院 | 一种脑脊液引流及脑温自动监测装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4548516A (en) * | 1984-12-24 | 1985-10-22 | Helenowski Tomasz K | Sequential thermography fluid flow indicator |
US20050204811A1 (en) * | 2004-02-03 | 2005-09-22 | Neff Samuel R | System and method for measuring flow in implanted cerebrospinal fluid shunts |
US20070167867A1 (en) * | 2005-02-24 | 2007-07-19 | Erich Wolf | System for transcutaneous monitoring of intracranial pressure |
US20080125690A1 (en) * | 2006-09-24 | 2008-05-29 | Delaporte Stephen E | Electroactive polymer actuated cerebrospinal fluid shunt |
US20130109998A1 (en) * | 2010-05-28 | 2013-05-02 | Shuntcheck, Inc. | Real time csf flow measurement system & method |
US20150297093A1 (en) * | 2014-04-18 | 2015-10-22 | Vivonics, Inc. | Flow rate sensor system and method for non-invasively measuring the flow rate of a bodily fluid |
CN106595904A (zh) * | 2016-12-16 | 2017-04-26 | 莆田学院 | 已知特征曲线的温度传感器在同性介质中的质量检测方法 |
-
2017
- 2017-09-28 CN CN201710896666.3A patent/CN107822594B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4548516A (en) * | 1984-12-24 | 1985-10-22 | Helenowski Tomasz K | Sequential thermography fluid flow indicator |
US20050204811A1 (en) * | 2004-02-03 | 2005-09-22 | Neff Samuel R | System and method for measuring flow in implanted cerebrospinal fluid shunts |
US20070167867A1 (en) * | 2005-02-24 | 2007-07-19 | Erich Wolf | System for transcutaneous monitoring of intracranial pressure |
US20080125690A1 (en) * | 2006-09-24 | 2008-05-29 | Delaporte Stephen E | Electroactive polymer actuated cerebrospinal fluid shunt |
US20130109998A1 (en) * | 2010-05-28 | 2013-05-02 | Shuntcheck, Inc. | Real time csf flow measurement system & method |
US20150297093A1 (en) * | 2014-04-18 | 2015-10-22 | Vivonics, Inc. | Flow rate sensor system and method for non-invasively measuring the flow rate of a bodily fluid |
CN106595904A (zh) * | 2016-12-16 | 2017-04-26 | 莆田学院 | 已知特征曲线的温度传感器在同性介质中的质量检测方法 |
Non-Patent Citations (1)
Title |
---|
SATHISH RAJASEKARAN等: "Thermal Measurement of Cerebrospinal Fluid Flow Rate in Hydrocephalus Shunt", 《IEEE》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109979295A (zh) * | 2019-05-07 | 2019-07-05 | 莆田学院 | 一种脊髓型颈椎病的病变演示装置及演示方法 |
CN114366859A (zh) * | 2021-10-29 | 2022-04-19 | 首都医科大学宣武医院 | 一种脑脊液引流及脑温自动监测装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107822594B (zh) | 2020-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5738986B2 (ja) | 埋め込み型脳脊髄液(csf)シャントにおけるcsf流量を決定するための装置 | |
CN104812324B (zh) | 温度分布确定装置 | |
Kolios et al. | Experimental evaluation of two simple thermal models using transient temperature analysis | |
US20130035921A1 (en) | Method and apparatus for estimating temperature in a body | |
US8551011B2 (en) | Cerebrospinal fluid evaluation system having thermal flow and flow rate measurement pad using a plurality of control sensors | |
CN107822594B (zh) | 一种脑脊液分流流量检测系统及方法 | |
US7520862B2 (en) | Cerebral spinal fluid shunt evaluation system | |
US20150201882A1 (en) | Continuous real-time csf flow monitor and method | |
Valvano et al. | Analysis of the Weinbaum-Jiji model of blood flow in the Canine Kidney cortex for self-heated thermistors | |
Bronskill et al. | A quantitative measurement of peritoneal drainage in malignant ascites | |
Liu et al. | Sinusoidal heating method to noninvasively measure tissue perfusion | |
US20140081157A1 (en) | Apparatus and Method for determining a volume amount of a physiological volume | |
US9072866B2 (en) | Cerebrospinal fluid evaluation system having thermal flow and flow rate measurement pad using a plurality of control sensors | |
Delhomme et al. | Thermal diffusion probe and instrument system for tissue blood flow measurements: validation in phantoms and in vivo organs | |
CN115252961B (zh) | 一种临床点滴输液加温与流速监控系统 | |
Yue et al. | Determination of the 3D temperature distribution during ferromagnetic hyperthermia under the influence of blood flow | |
JP2013516199A5 (ja) | 生理的ボリューム量を決定するための装置と方法 | |
Schooneveldt et al. | Experimental validation of a thermophysical fluid model for use in a hyperthermia treatment planning system | |
Mudaliar et al. | A phantom tissue system for the calibration of perfusion measurements | |
WO2018101396A1 (ja) | 超音波画像化装置、超音波画像化方法、超音波画像化プログラム、及び超音波ct装置 | |
CN107713987B (zh) | 脑脊液分流流量检测方法 | |
Robinson et al. | Validation of methodologies for the estimation of blood perfusion using a minimally invasive probe | |
Tsang et al. | The assessment of hand volume using a modified volumetric technique | |
Stewart et al. | Investigation of a non-invasive venous blood flow measurement device: Using thermal mass measurement principles | |
EP4066261A1 (en) | Method for estimating heat transfer energy parameters in an encephalon |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201013 Termination date: 20210928 |