本发明专利技术提出了一种基于递推方式的实时快速磁纳米温度测量方法，步骤为：将磁纳米粒子样品放置在待测对象处，对磁纳米粒子样品所在的区域同时施加交流和直流磁场，利用空心线圈实时获取磁纳米粒子样品在施加磁场激励下的磁化信息，通过低噪声前置放大器进行信号放大，采用数据采集卡对信号进行离散采集，利用谐波提取算法检测出各次谐波信号的幅值信息；根据郎之万函数和二次谐波与一次谐波幅值之比导函数的特性建立递推公式，构建实时温度测量递推模型，获取到温度信息。本发明专利技术可以实时快速的获取到被测非接触式物体温度信息，从根本上解决了磁纳米温度测量实时性较低的问题，同时又避开了高次谐波测量困难的难题，实现快速实时温度测量。
【技术实现步骤摘要】
本专利技术涉及纳米测试技术和非侵入式温度测量的
，特别是指一种基于递推方式的实时快速磁纳米温度测量方法，适用于肿瘤癌症等细胞内部温度信息测量。
技术介绍
肿瘤是严重威胁人类身体健康甚至危及生命的疾病之一，据世界卫生组织统计，肿瘤癌症已经成为人类死亡的疾病之首。世界科研与技术工作者为此提出治疗手段，如手术切除、放射与化学疗法。然而手术切除只能针对早期患者，放射与化学疗法给病人带来极大痛苦的同时治愈率较低。世界各国都对热疗进行了大量的临床试验，“肿瘤热疗法”被称为治愈恶性肿瘤最有希望的手段之一，同时也被誉称为“绿色疗法”。肿瘤细胞的耐热性较低，一般在温度为42.5℃左右时，保持一定时间肿瘤组织细胞即可死亡，而正常细胞则可以较长时间耐热此温度。但是，目前一些问题还有待解决，如肿瘤癌症细胞温度的实时精密测量。临床常用于肿瘤热疗的测温方法如光钎温度传感器、热电阻和热电偶等，存在以下局限：首先易受传感器数量的限制，无法精确掌握整个肿瘤组织的温度场；其次温度测量探头本身会对肿瘤组织加热产生影响，使加热不均匀；最后这些温度传感器需要植入患者体内，给患者带来极大的痛苦体验，使治疗风险加大。此外，药物治疗肿瘤细胞时，需要保证肿瘤组织区域在合适的药物浓度下并维持一定时间，才能达到最佳的治疗效果，而温度则是控制靶向磁纳米药物热敏载体释放药物的主要因素。药物热敏载体在温度为42℃左右时，释放药物速率能达到峰值，可以迅速使肿瘤组织的药物浓度上升，此时肿瘤药物的毒性最强，当温度变化幅度在2℃时，释放药物的速率将会下降超出峰值的50％左右。当然，如果温度一直保持较高，会使得肿瘤组织药物浓度过高，肿瘤细胞来不及吸收，多余的药物将会随着血液流向其他正常组织器官，造成对身体的损害，所以需要在一定时间内将温度降下来。所以，如果温度的掌控精度不够高，将直接影响药物的释放，导致肿瘤组织药物浓度不足或者长时间过高，不能有效杀死肿瘤细胞或损伤其他组织器官而影响肿瘤的治疗效果，增加医治风险；而药物浓度不足很容易诱使它们产生抗体，导致产生抗药性，增加医治难度；加热温度过高，又容易发生热损伤情况。因此，体内组织细胞温度实时测量技术是目前肿瘤热疗方法亟待解决的重要研究课题。磁纳米材料的出现给上述问题提供一种解决方法。美国学者J.B.Weaver于2009年通过实验的方式验证了磁纳米粒子的温度敏感性，即磁纳米粒子在不同温度下的磁化响应不同。研究发现在单频交变磁场激励下，磁纳米粒子的交流磁化强度信息中的三次谐波幅值和五次谐波幅值的比值与温度具有较好的线性度，遗憾的是其缺少相关的理论依据支撑。华中科技大学刘文中教授领导的科研小组于2011年提出了磁纳米测温模型，从理论上给磁纳米粒子温度测量方法提供依据，研究发现，在直流磁场激励下的磁纳米粒子磁化强度倒数具有极强的温度敏感性，并提出了根据磁化率倒数进行温度测量的理论模型，并成功测量到磁纳米温度的信息。只是该测量方法测量时间较长，无法满足医疗应用要求。钟景博士通过对磁纳米粒子在单频交流磁场激励下的磁化强度温度敏感性进行研究，提出了单频交变磁场激励下的温度测量模型，该方法在一定程度上提高了测温实时性，然而该方法要求的磁纳米粒子磁化响应的高次谐波信息测量难度较大。
技术实现思路
本专利技术为了解决了磁纳米温度测量实时性较低的技术问题，提出一种基于递推方式的实时快速磁纳米温度测量方法，能够对磁纳米实现快速实时的温度测量，从而满足肿瘤癌症热疗中所要求的快速实时温度测量。为了解决上述技术问题，本专利技术的技术方案是：一种基于递推方式的实时快速磁纳米温度测量方法，其步骤在于：步骤一：将磁纳米样品放置于待测对象处；步骤二：对磁纳米样品所在区域施加直流激励和交流激励磁场：H＝Hdc+H1sin(ω1t)，其中，Hdc是直流磁场强度，H1是频率为ω1的交流磁场强度；步骤三：采用空心线圈作为磁探测传感器采集直流激励磁场和交流激励磁场共同激励下磁纳米样品的磁化强度信号；步骤四：提取k时刻温度Tk、k+1时刻温度Tk+1、k+2时刻温度Tk+2下的磁纳米样品以频率ω1为基频下的磁化强度信号的一次谐波的幅值M1(Tk)、M1(Tk+1)、M1(Tk+2)，提取k时刻温度Tk、k+1时刻温度Tk+1、k+2时刻温度Tk+2下的磁纳米样品以频率ω1为基频下的磁化强度信号的二次谐波的幅值M2(Tk)、M2(Tk+1)、M2(Tk+2)；其中，k＝1,2,3，…；步骤五：根据郎之万函数建立二次谐波幅值与一次谐波幅值之比和温度之间的函数关系构建其导函数：其中，T为待测对象温度，Ms为磁纳米粒子有效磁矩，kB为波尔兹曼常数；对进行泰勒级数展开，整理得到温度递推公式：其中，Tk为k时刻的温度；步骤六：根据步骤四中获取的k、k+1、k+2时刻的一次谐波和二次谐波幅值之比以及温度值Tk和Tk+1作为步骤五所得到递推公式的初值，带入步骤五的递推公式获取到k+2时刻的磁纳米温度值Tk+2；步骤七：k的值加1，循环步骤四-五获取到k+2时刻的磁纳米温度值Tk+2。所述直流激励和交流激励磁场利用通电的亥姆霍兹线圈或者螺线管产生，交流激励磁场强度H1的范围在5Gs以下，直流激励磁场强度Hdc大于10Gs。所述空心线圈为空心式差分线圈，空心线圈的采集的磁化强度信号送入低噪声前置放大器对其进行信号放大和滤波预处理，然后通过数据采集卡进行离散采集。所述步骤四中利用数字相敏检波算法或快速傅立叶变换算法提取磁化强度信号的各次谐波幅值。所述郎之万函数泰勒展开项数m一般取值范围为2-8。本专利技术的有益效果：将磁纳米样品(固体/液体)放置在待测对象处，然后施加直流磁场和交流磁场激励共同激励，磁纳米粒子的磁化强度信息包含丰富的各次谐波信息，通过数字相敏检波算法提取出所需各次谐波的幅值，采用谐波幅值比(如二次谐波与一次谐波幅值比，三次谐波与一次谐波幅值比等)构建磁纳米温度递推公式，根据前面两个时刻的温度和谐波幅值之比值通过温度递推公式测量当前时刻的温度值。本专利技术仅仅利用磁纳米粒子前两个时刻的温度值和谐波比值进行当前温度的递推，大大提高了磁纳米温度测量的时间分辨率，有望解决肿瘤癌症热疗中实时快速温度测量问题，而且还适用于其他非侵入式快速温度测量场合。本专利技术利用递推方式进行温度信息测量，不需要进行复杂的温度反演求解过程，大大提高了温度测量的实时性，单个温度点测量时间仅仅取决于数据采集卡对信号离散采集速度，而以往单个磁纳米温度测量时间为1秒。本专利技术可以实时快速的获取到被测物体温度信息，尤其是可以实现非接触式物体内部温度信息的测量；从根本上解决了磁纳米温度测量实时性较低的问题，同时又避开了高次谐波测量困难的难题，实现快速实时温度测量。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术的具体流程图。图2为310-320K温度范围内(步进为1K)二次与一次谐波幅值之比随着温度变化的仿真图。图3为310-311K温度范围内(步进为0.1K)二次与一次谐波幅值之比随着温度变化的本文档来自技高网...
【技术保护点】
一种基于递推方式的实时快速磁纳米温度测量方法，其特征在于，其步骤在于：步骤一：将磁纳米样品放置于待测对象处；步骤二：对磁纳米样品所在区域施加直流激励和交流激励磁场：H＝Hdc+H1sin(ω1t)，其中，Hdc是直流磁场强度，H1是频率为ω1的交流磁场强度；步骤三：采用空心线圈作为磁探测传感器采集直流激励磁场和交流激励磁场共同激励下磁纳米样品的磁化强度信号；步骤四：提取k时刻温度Tk、k+1时刻温度Tk+1、k+2时刻温度Tk+2下的磁纳米样品以频率ω1为基频下的磁化强度信号的一次谐波的幅值M1(Tk)、M1(Tk+1)、M1(Tk+2)，提取k时刻温度Tk、k+1时刻温度Tk+1、k+2时刻温度Tk+2下的磁纳米样品以频率ω1为基频下的磁化强度信号的二次谐波的幅值M2(Tk)、M2(Tk+1)、M2(Tk+2)；其中，k＝1,2,3，…；步骤五：根据郎之万函数建立二次谐波幅值与一次谐波幅值之比和温度之间的函数关系构建其导函数：M2M1(T)=F(T)=H12Hdc30(MskBT)3-(H14Hdc189+2H12Hdc3189)(MskBT)5+...H13MskBT-(H1360+H1Hdc215)(MskBT)3+(H15756+2H1Hdc4189+H13Hdc263)(MskBT)5+...2(M2M1(T))dT=f(T)=d(H12Hdc30(MskBT)3-(H14Hdc189+2H12Hdc3189)(MskBT)5+...H13MskBT-(H1360+H1Hdc215)(MskBT)3+(H15756+2H1Hdc4189+H13Hdc263)(MskBT)5+...)/dT;]]>其中，T为待测对象温度，Ms为磁纳米粒子有效磁矩，kB为波尔兹曼常数；对进行泰勒级数展开，整理得到温度递推公式：Tk+2=M2M1(Tk+2)-M2M1(Tk+1)M2M1(Tk+1)-M2M1(Tk)(Tk+1-Tk)+Tk,k=1,2,3...;]]>其中，Tk为k时刻的温度；步骤六：根据步骤四中获取的k、k+1、k+2时刻的一次谐波和二次谐波幅值之比以及温度值Tk和Tk+1作为步骤五所得到递推公式的初值，带入步骤五的递推公式获取到k+2时刻的磁纳米温度值Tk+2；步骤七：k的值加1，循环步骤四‑五获取到k+2时刻的磁纳米温度值Tk+2。...
【技术特征摘要】
1.一种基于递推方式的实时快速磁纳米温度测量方法，其特征在于，其步骤在于：步骤一：将磁纳米样品放置于待测对象处；步骤二：对磁纳米样品所在区域施加直流激励和交流激励磁场：H＝Hdc+H1sin(ω1t)，其中，Hdc是直流磁场强度，H1是频率为ω1的交流磁场强度；步骤三：采用空心线圈作为磁探测传感器采集直流激励磁场和交流激励磁场共同激励下磁纳米样品的磁化强度信号；步骤四：提取k时刻温度Tk、k+1时刻温度Tk+1、k+2时刻温度Tk+2下的磁纳米样品以频率ω1为基频下的磁化强度信号的一次谐波的幅值M1(Tk)、M1(Tk+1)、M1(Tk+2)，提取k时刻温度Tk、k+1时刻温度Tk+1、k+2时刻温度Tk+2下的磁纳米样品以频率ω1为基频下的磁化强度信号的二次谐波的幅值M2(Tk)、M2(Tk+1)、M2(Tk+2)；其中，k＝1,2,3，…；步骤五：根据郎之万函数建立二次谐波幅值与一次谐波幅值之比和温度之间的函数关系构建其导函数：M2M1(T)=F(T)=H12Hdc30(MskBT)3-(H14Hdc189+2H12Hdc3189)(MskBT)5+...H13MskBT-(H1360+H1Hdc215)(MskBT)3+(H15756+2H1Hdc4189+H13Hdc263)(MskBT)5+...2(M2M1(T))dT=f(T)=d(H12Hdc30(MskBT)3-(H14Hdc189+2H12Hdc3189)(MskBT)5+...H13MskBT-(H1360+H1Hdc215)(MskBT)3+(H15756+2...
【专利技术属性】
技术研发人员：杜中州，苏日建，史江山，刘文中，甘勇，皮仕强，叶娜，郭功兵，张涛，
申请(专利权)人：郑州轻工业学院，
类型：发明
国别省市：河南;41
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。