发明名称 一种无损测量生物活体组织热参数的方法
摘要 本发明提供了一种无损测量生物活体组织热参数的方法,用于通过三点表面测温确定柱状生物活体组织的热参数。该方法包括:建立柱坐标下柱状生物活体组织的二维Pennes传热模型;采用逐行法求解在阶跃热流或方波热流的干扰作用下,柱状生物活体组织内的温度变化;通过灵敏度计算,分析热参数的灵敏度系数之间的相关性以及热参数对温度变化的影响,优化三个测量点的选择;开发生物活体组织表面温度的测量系统,实时采集温度的瞬态响应数据;应用改进的高斯参数估计方法,同时获得柱状生物活体组织的导热系数、热扩散系数、血液灌注率以及体积热容量等热参数。其优点在于:实现了方便、步骤简单的无损测量。
申请公布号 CN1555759A 申请公布日期 2004.12.22
申请号 CN200410000021.X 申请日期 2004.01.02
申请人 北京科技大学 发明人 张欣欣;乐恺;于帆;姜泽毅;冯研卉;林林
分类号 A61B5/00;A61B5/01 主分类号 A61B5/00
代理机构 北京科大华谊专利代理事务所 代理人 刘月娥
主权项 1.一种无损测量生物组织热参数的方法,用于通过三点表面测温确定柱状生物活体组织的热参数,其特征在于:对长度为0.10m~0.75m、半径为0.01m~0.20m的柱状生物活体组织,建立了特定形式热干扰作用下的二维Pennes柱坐标传热模型;具体步骤如下:a、建立特定形式热干扰作用下,柱坐标下柱状生物活体组织的二维Pennes传热模型其控制方程式:<math> <mrow> <msub> <mrow> <mo>(</mo> <mi>&rho;c</mi> <mo>)</mo> </mrow> <mi>t</mi> </msub> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>T</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mi>k</mi> <mo>[</mo> <mfrac> <mrow> <msup> <mo>&PartialD;</mo> <mn>2</mn> </msup> <mi>T</mi> </mrow> <msup> <mrow> <mo>&PartialD;</mo> <mi>x</mi> </mrow> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mi>r</mi> </mfrac> <mfrac> <mo>&PartialD;</mo> <mrow> <mo>&PartialD;</mo> <mi>r</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>r</mi> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>T</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>r</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>]</mo> <mo>+</mo> <msub> <mi>w</mi> <mi>b</mi> </msub> <msub> <mi>c</mi> <mi>b</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>a</mi> </msub> <mo>-</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>q</mi> <mi>m</mi> </msub> </mrow> </math> (式1)其边界条件如式2所示:<math> <mrow> <mi>x</mi> <mo>=</mo> <mn>0,0</mn> <mo>&lt;</mo> <mi>r</mi> <mo>&lt;</mo> <mi>R</mi> <mo>,</mo> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>T</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>x</mi> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>;</mo> <mi>x</mi> <mo>=</mo> <mi>L</mi> <mo>,</mo> <mn>0</mn> <mo>&lt;</mo> <mi>r</mi> <mo>&lt;</mo> <mi>R</mi> <mo>,</mo> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>T</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>x</mi> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>;</mo> <mi>r</mi> <mo>=</mo> <mn>0,0</mn> <mo>&le;</mo> <mi>x</mi> <mo>&le;</mo> <mi>L</mi> <mo>,</mo> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>T</mi> </mrow> <mrow> <mo>&PartialD;</mo> <mi>r</mi> </mrow> </mfrac> <mo>=</mo> <mn>0</mn> <mo>;</mo> </mrow> </math> (式2)其中,(ρc)t为生物组织的体积热容量,k为导热系数,wb为血液灌注率,cb为血液比热容,Ta为动脉血液温度,qm为体积代谢热,qo为加热带热流密度,hA为边界的综合换热系数,T∞为环境温度。并且采用无其它热干扰的二维稳态计算结果作为初始温度场;b、采用数值方法求解柱状生物活体组织内的温度变化;c、计算参数的灵敏度系数,进行灵敏度分析;d、实时采集温度的瞬态响应数据;e、应用参数估计方法,确定柱状生物活体组织的导热系数、热扩散系数、血液灌注率以及体积热容量等热参数。
地址 100083北京市海淀区学院路30号
您可能感兴趣的专利