CN109363826B - 一种基于植入线圈的肿瘤热疗装置及其参数优化方法 - Google Patents

一种基于植入线圈的肿瘤热疗装置及其参数优化方法 Download PDF

Info

Publication number
CN109363826B
CN109363826B CN201811123473.5A CN201811123473A CN109363826B CN 109363826 B CN109363826 B CN 109363826B CN 201811123473 A CN201811123473 A CN 201811123473A CN 109363826 B CN109363826 B CN 109363826B
Authority
CN
China
Prior art keywords
coil
vivo
radius
formula
implanted
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.)
Active
Application number
CN201811123473.5A
Other languages
English (en)
Other versions
CN109363826A (zh
Inventor
程瑜华
陈国雄
王高峰
李文钧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wenzhou Huidian Technology Co ltd
Original Assignee
Hangzhou Dianzi University Wenzhou Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University Wenzhou Research Institute Co Ltd filed Critical Hangzhou Dianzi University Wenzhou Research Institute Co Ltd
Priority to CN201811123473.5A priority Critical patent/CN109363826B/zh
Publication of CN109363826A publication Critical patent/CN109363826A/zh
Application granted granted Critical
Publication of CN109363826B publication Critical patent/CN109363826B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F7/00—Heating or cooling appliances for medical or therapeutic treatment of the human body
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N2/00—Magnetotherapy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
    • A61F7/00—Heating or cooling appliances for medical or therapeutic treatment of the human body
    • A61F2007/009—Heating or cooling appliances for medical or therapeutic treatment of the human body with a varying magnetic field acting upon the human body, e.g. an implant therein

Abstract

本发明公开了一种基于植入线圈的肿瘤热疗装置及其参数优化方法。传统的肿瘤磁介导热疗法中需要向人体注射产热的磁性颗粒；且随着颗粒的扩散，热疗的效果降低。本发明一种基于植入线圈肿瘤磁介导热疗装置，包括体内线圈和体外线圈。所述的体外线圈为开环的单圈圆形线圈。所述的体内线圈为闭环的单圈圆形线圈。体外线圈、体内线圈的半径分别为r1和r2，线半径分别为rw1和rw2。
Figure DDA0001811785330000011
其中，d为无线传输距离。
Figure DDA0001811785330000012
本发明能够对肿瘤达到可持续性热疗的效果。

Description

一种基于植入线圈的肿瘤热疗装置及其参数优化方法
技术领域
本发明属于生物医学电子技术领域，具体涉及一种基于植入线圈肿瘤磁介导热疗的加热效率优化方法。
背景技术
传统的肿瘤磁介导热疗法是一种通过直接注射、静脉注射或介入等方式使产热材料定向聚集在肿瘤部位，并在交变磁场的作用下发生磁致产热效应将肿瘤组织加热到41℃以上治疗肿瘤的方法。磁性纳米或微米颗粒是当前广泛使用的一种产热介质。尽管通过注射磁性颗粒的方式是一种较小侵入性的方式，但是热疗的效果依赖于磁性颗粒的剂量，而且随着颗粒的扩散，热疗的效果也会随着降低，并且磁性颗粒扩散到人体其他部位后容易造成二次伤害。
发明内容
本发明的目的在于提供一种基于植入线圈的肿瘤热疗装置及其参数优化方法。
本发明给出了一种基于植入线圈的肿瘤热疗装置，包括体内线圈和体外线圈。所述的体外线圈为开环的单圈圆形线圈。所述的体内线圈为闭环的单圈圆形线圈。体外线圈、体内线圈的半径分别为r1和r2，线半径分别为rw1和rw2。r2与rw2满足等式
Figure BDA0001811785310000011
其中，f为装置工作频率，根据应用背景在1MHz至20MHz之间进行选取；μ0为真空中磁导率；σ2为体内线圈的电导率；
Figure BDA0001811785310000021
其中，d为无线传输距离(工作状态下，发射线圈与接收线圈之间的距离)。线半径rw1的取值范围为
Figure BDA0001811785310000022
进一步地，本发明的一种基于植入线圈肿瘤磁介导热疗装置还包括信号发生器。信号发生器的信号输出端、地线端与体外线圈的两个接线端分别相连。所述的体内线圈采用生物金属材料。所述的生物金属材料为钛合金。
进一步地，所述体内线圈的半径r2等于待植入肿瘤块直径的1/5。
进一步地，所述体内线圈的半径r2等于5mm。无线传输距离d等于30mm。
该基于植入线圈的肿瘤热疗装置的参数优化方法具体如下：
步骤一、求取体外线圈半径
Figure BDA0001811785310000023
其中，r2为体内线圈的半径；d为无线传输距离。
步骤二、计算rw1＝s·r1；其中
Figure BDA0001811785310000024
步骤三、联立Q2＝1与
Figure BDA0001811785310000025
其中，μ0为真空中磁导率；f为装置工作频率；σ2为体内线圈的电导率。得到体内线圈的最优线半径rw2。
进一步地，步骤一执行前，根据待植入肿瘤的直径确定体内线圈的半径r2，使得体内线圈能够完全没入待植入肿瘤内，使得r2等于待植入肿瘤直径的1/5。根据待植入肿瘤的位置，确定无线传输距离d，使得d等于待植入肿瘤几何中心到人体外表面的最小距离加上5mm所得值。
进一步地，步骤三中，f取f1、f2、…、fn；f1、f2、…、fn均小于20MHz；分别进行体内线圈的线半径rw2的计算；得到n个体内线圈的线半径候选值rw2(1)、rw2(2)、…、rw2(n)。之后执行以下步骤。
3-1.i＝1,2,…,n依次执行步骤3-2至3-6。
3-2.计算体外线圈的电感值L1如式(1)所示、体内线圈的电感值L2如式(2)所示；
L1＝μ0r1(ln(8r1/rw1)-2) 式(1)
L2＝μ0r2(ln(8r2/rw2(i))-2) 式(2)
3-3.计算体外线圈的寄生电阻值R1如式(3)所示、体内线圈的寄生电阻值R2如式(4)所示；
R1＝r1/σ1δ1rw1 式(3)
R2＝r2/σ2δ2rw2(i) 式(4)
式(3)和式(4)中，δ1为体外线圈的趋肤深度，其表达式为
Figure BDA0001811785310000031
δ2为体内线圈的趋肤深度，其表达式为
Figure BDA0001811785310000032
其中σ1为体外线圈的电导率；σ2为体内线圈的电导率。
3-4建立体外线圈与体内线圈之间的耦合系数为k与体外线圈的半径r1、线半径rw1以及体内线圈的半径r2、线半径rw2的关系式如式(5)所示。
Figure BDA0001811785310000033
式(5)中，d为无线传输距离；c1的表达式为c1＝ln(8r1/rw1)-2；c2的表达式为c2＝ln(8r2/rw2(i))-2。
3-5、建立体外线圈的品质因数Q1的表达式如式(6)所示、体内线圈的品质因数Q2的表达式如式(7)所示
Q1＝ωL1/Rx1 式(6)
Q2＝ωL2/R2 式(7)
式(6)和式(7)中，ω＝2πf。Rx1为体外线圈的总和阻抗。(Q1的大小不影响体外线圈、体内线圈形状参数的计算，仅影响传输效率计算。)
3-6.建立耦合系数为k、体外线圈的品质因数Q1、体内线圈的品质因数Q2与体内线圈、体外线圈之间电能传输效率ηi之间的关系式如式(9)所示；
Figure BDA0001811785310000041
3-7.取η1、η2、…、ηn中的最大值对应的那个体内线圈的线半径候选值作为最终的体内线圈的线半径rw2。
进一步地，n＝4，且f1＝1MHz，f2＝2.5MHz，f3＝5MHz，f4＝10MHz。
进一步地，Rx1＝R1。
进一步地，Rx1＝R1+；其中ΔZ1的表达式如下：
Figure BDA0001811785310000042
其中，j为复数符号；J1(k)是以k为变量的1阶第一类贝塞尔函数；
Figure BDA0001811785310000043
是以
Figure BDA0001811785310000044
为变量的1阶第一类贝塞尔函数；k′2的表达式为k′2＝ωμ0σtr1 2；σt是人体组织的电导率；装置工作频率为1MHz时，σt＝0.50268；装置工作频率为2.5MHz时，σt＝0.55928；装置工作频率为5MHz时，σt＝0.59008；装置工作频率为10MHz时，σt＝0.61683。
本发明具有的有益效果是：
1、本发明通过向肿瘤中植入线圈的方式，能够对肿瘤达到可持续性热疗的效果。
2、本发明通过用线圈代替磁性颗粒的方式避免了磁性颗粒扩散对人体造成伤害的问题。
3、本发明通过对体外线圈的半径、线半径以及体内线圈的线半径进行优化，具有较高的无线传输效率，能够有效减少目标区域外人体组织的热损耗。
附图说明
图1为本发明的工作示意图；
图2为本发明的等效电路图；
图3为本发明的一个实例中在装置工作频率分别为1MHz、2.5MHz、5MHz、10MHz计算出的传输效率与体内线圈线半径的函数关系图；
图4为本发明的一个实例中以1MHz作为装置工作频率时得到的热分布仿真图；
图5为本发明的一个实例中以2.5MHz作为装置工作频率时得到的热分布仿真图；
图6为本发明的一个实例中以5MHz作为装置工作频率时得到的热分布仿真图；
图7为本发明的一个实例中以10MHz作为装置工作频率时得到的热分布仿真图。
具体实施方式
以下结合附图对本发明作进一步说明。
如图1所示，一种基于植入线圈肿瘤磁介导热疗装置，包括信号发生器、体内线圈2和体外线圈1。体外线圈1为开环的单圈圆形线圈(即有一处断开)。体外线圈1断开的两个接线端与信号发生器的信号输出端、地线端分别相连。体内线圈2为闭环的单圈圆形线圈(即体内线圈首尾相连)。体外线圈1、体内线圈2的半径分别为r1和r2，所使用的导线的线半径分别为rw1和rw2。体内线圈采用生物金属材料(本实施例中采用钛合金)。r2等于待植入肿瘤直径的1/5。r2与rw2满足等式
Figure BDA0001811785310000061
其中，f为装置工作频率，取值为1MHz、2.5MHz、5MHz或10MHz；μ0为真空中磁导率，取值为4π×10-7H/m；σ2为体内线圈的电导率；
Figure BDA0001811785310000062
其中，d为无线传输距离，取值等于待植入肿瘤3几何中心到人体外表面的最小距离。
Figure BDA0001811785310000063
体内线圈2埋设在待植入肿瘤3内。
该基于植入线圈的肿瘤热疗装置的参数优化方法具体如下：
步骤一、根据待植入肿瘤的直径确定体内线圈的半径r2，使得体内线圈能够完全没入待植入肿瘤内，本实施例中，使得r2等于待植入肿瘤直径的1/5。在待植入肿瘤未知的情况下，r2取5mm。根据待植入肿瘤的位置，确定无线传输距离d，使得d等于待植入肿瘤几何中心到人体外表面的最小距离加上5mm所得值。在待植入肿瘤几何中心到人体外表面的最小距离未知的情况下，d取30mm。
建立体外线圈的电感值L1与体外线圈的半径r1、线半径rw1的关系式如式(1)所示，体内线圈的电感值L2与体内线圈的半径r2、线半径rw2的关系式如式(2)所示；
L1＝μ0r1(ln(8r1/rw1)-2) 式(1)
L2＝μ0r2(ln(8r2/rw2)-2) 式(2)
式(1)和式(2)中，μ0为真空中磁导率，取值为4π×10-7H/m。
步骤二、建立体外线圈的寄生电阻值为R1与体外线圈的半径r1、线半径rw1的关系式如式(3)所示；体内线圈的寄生电阻值为R2与体外线圈的半径r2、线半径rw2的关系式如式(4)所示；
Figure BDA0001811785310000071
Figure BDA0001811785310000072
式(3)和式(4)中，m1的表达式为
Figure BDA0001811785310000073
μr1为体外线圈的磁导率；σ1为体外线圈的电导率；RDC1为体外线圈直流电阻(即通直流电时的电阻值)，其表达式为
Figure BDA0001811785310000074
m2的表达式为
Figure BDA0001811785310000075
μr2为体内线圈的磁导率；σ2为体内线圈的电导率；RDC2为体内线圈直流电阻(即通直流电时的电阻值)，其表达式为
Figure BDA0001811785310000076
ber(m)是以m为变量的0阶第一类开尔文函数的实部；ber′(m)是以m为变量的0阶第一类开尔文函数实部的导函数。bei(m)是以m为变量的0阶第一类开尔文函数的虚部；bei′(m)是以m为变量的0阶第一类开尔文函数虚部的导函数。
式(3)、式(4)中不容易直观发现体外线圈的寄生电阻值为R1与rw1之间的关系。由于单圈线圈的寄生电阻主要是由趋肤效应引起的，因此可以将寄生电阻R1的表达式简化为：R1≈2πr1/σ1δ1(2πrw1)＝r1/σ1δ1rw1；寄生电阻R2的表达式简化为：R2≈2πr2/σ2δ2(2πrw2)＝r2/σ2δ2rw2；δ1为体外线圈的趋肤深度，其表达式为
Figure BDA0001811785310000077
δ2为体内线圈的趋肤深度，其表达式为
Figure BDA0001811785310000078
其中，f为装置工作频率，取值为1MHz、2.5MHz、5MHz或10MHz；
步骤三、建立体外线圈与体内线圈之间的耦合系数为k与体外线圈的半径r1、线半径rw1以及体内线圈的半径r2、线半径rw2的关系式如式(5)所示。
Figure BDA0001811785310000079
式(5)中，M为互感系数，其表达式如式(6)所示；d为无线传输距离；c1的表达式为c1＝ln(8r1/rw1)-2；c2的表达式为c2＝ln(8r2/rw2)-2。
Figure BDA0001811785310000081
步骤四、建立体外线圈的品质因数Q1的表达式如式(7)所示、体内线圈的品质因数Q2的表达式如式(8)所示
Q1＝ωL1/R1 式(7)
Q2＝ωL2/R2 式(8)
式(7)和式(8)中，ω为装置工作角频率，取值为ω＝2πf。
由于体外线圈通交流电后会在人体组织内会产生涡流。体外线圈在人体组织中引起的涡流损耗可以折算到体外线圈的阻抗变化。体外线圈在有人体组织和无人体组织两种情况下阻抗的变化量ΔZ1可以计算为
Figure BDA0001811785310000082
其中，j为复数符号；J1(k)是以k为变量的1阶第一类贝塞尔函数；
Figure BDA0001811785310000083
是以
Figure BDA0001811785310000084
为变量的1阶第一类贝塞尔函数；k′2的表达式为k′2＝ωμ0σtr1 2；σt是人体组织的电导率(用一种典型的人体组织的电导率代替实际情况中多层人体组织情况，装置工作频率为1MHz时，σt＝0.50268；装置工作频率为2.5MHz时，σt＝0.55928；装置工作频率为5MHz时，σt＝0.59008；装置工作频率为10MHz时，σt＝0.61683)。ΔZ1的实部会增加寄生电阻值大小，ΔZ1的虚部会减小电感值大小。利用此阻抗的变化值，得到Q1的表达式：Q1＝ωL1/(R1+ΔZ1)；进而得到更精准的效率值。
步骤五、建立耦合系数为k、体外线圈的品质因数Q1、体内线圈的品质因数Q2与体内线圈、体外线圈之间电能传输效率η之间的关系式如式(9)所示；
Figure BDA0001811785310000091
步骤六、最大化耦合系数k，得到最优体外线圈半径r1。从式(9)可以看出，电能传输效率η是耦合系数k的单调递增函数，因此，增大耦合系数k，能够增大电能传输效率η。
求取耦合系数k对体外线圈半径r1的导函数，并对所得导函数求零点，得到耦合系数k取最大值时对应体外线圈半径
Figure BDA0001811785310000092
由于体内线圈半径r2和传输距离d均根据肿瘤的情况决定，为已知值，故能够计算出确定的体外线圈半径r1。
步骤七、最大化体外线圈的品质因数Q1，得到最优体外线圈线半径rw1。从式(9)可以看出，电能传输效率η是体外线圈的品质因数Q1的单调递增函数，因此，增大体外线圈的品质因数Q1，能够增大电能传输效率η。
根据R1＝r1/σ1δrw1可以得到简化的
Figure BDA0001811785310000093
式中可以看出，体外线圈的品质因数Q1是体外线圈线半径rw1的单调递增函数；因此，体外线圈线半径rw1越大，则电能传输效率η越大。但是，过大的rw1将导致过小的电感值，使得前级驱动电路难以驱动。而且，为了保证安全性，rw1增大将增大发送线圈的中心位置与体内线圈中心位置的间距(原因在于：rw1越大，则体外线圈的中心位置越远离人体表面)，从而导致耦合系数k的降低，加热效率将会降低。因此，限定rw1的范围为
Figure BDA0001811785310000094
步骤八、将装置工作频率f的取值定为1MHz。
步骤九、确定体内线圈的线半径rw2。
求取
Figure BDA0001811785310000101
对Q2的偏导函数；并对所得导函数求零点，得到耦合系数η取最大值时对应的Q2＝1。联立Q2＝1与
Figure BDA0001811785310000102
得到体内线圈的线半径rw2。
步骤九、将装置工作频率f依次更换为2.5MHz、5MHz、10MHz；并分别执行步骤九，从而得到四个不同的体内线圈的线半径rw2。将所得的四个体内线圈的线半径rw2分别代入
Figure BDA0001811785310000103
得到四个传输效率η1、η2、η3、η4；取η1、η2、η3、η4中的最大值对应的那个装置工作频率作为最终的装置工作频率f。
本发明在一个实例的步骤八中分别为1MHz、2.5MHz、5MHz、10MHz作为装置工作频率，得到的半径rw2与工作频率的关系图如图3所示。实例中，体内线圈2的半径定为5mm；传输距离d定为30mm。从而计算得到体外线圈1的半径为30.55mm，线半径rw1为5mm；体内线圈2的线半径rw2为随装置工作频率不同而不同，具体取值如下表所示：
装置工作频率(MHz) 体内线圈2的线半径r<sub>w2</sub>(mm)
1 0.1765
2.5 0.1039
5 0.0700
10 0.0474
图3中，四个黑色圆点分别为计算得到的1MHz、2.5MHz、5MHz、10MHz对应的传输效率最高点；四个黑色方点分别为仿真得到的1MHz、2.5MHz、5MHz、10MHz对应的传输效率最高点。可以看出，尽管计算出的传输效率与仿真出的传输效率有一定的误差，但是计算得到的传输效率最大值所对应的横坐标(体内线圈线半径rw2)与测试得到的传输效率最大值所对应的横坐标(体内线圈线半径rw2)非常接近。由此可以看出本发明计算的传输效率最大值所对应的体内线圈线半径rw2与实际的传输效率最大值所对应的体内线圈线半径rw2基本相等。故本发明计算得到的体内线圈线半径rw2有助于提高本发明提出的一种基于植入线圈肿瘤磁介导热疗装置的发热效率。
针对所述实例进行0.2小时的加热模拟，以1MHz作为装置工作频率得到的热分布仿真图如图4所示；以2.5MHz作为装置工作频率得到的热分布仿真图如图5所示；以5MHz作为装置工作频率得到的热分布仿真图如图6所示；以10MHz作为装置工作频率得到的热分布仿真图如图7所示；从中可以看出，四种装置工作频率均能够达到无线传输加热的效果，且以5MHz作为装置工作频率时，加热效果最佳。

Claims (8)

1.一种基于植入线圈的肿瘤热疗装置，包括体内线圈和体外线圈；其特征在于：所述的体外线圈为开环的单圈圆形线圈；所述的体内线圈为闭环的单圈圆形线圈；体外线圈、体内线圈的半径分别为r1和r2，线半径分别为rw1和rw2；r2与rw2满足等式
Figure FDA0002613823620000011
其中，f为装置工作频率，f的取值范围为1MHz≤f≤20MHz；μ0为真空中磁导率；σ2为体内线圈的电导率；
Figure FDA0002613823620000012
其中，d为无线传输距离；
Figure FDA0002613823620000013
所述体内线圈的半径r2等于待植入肿瘤直径的1/5；无线传输距离d等于待植入肿瘤几何中心到人体外表面的最小距离。
2.根据权利要求1所述的一种基于植入线圈的肿瘤热疗装置，其特征在于：还包括信号发生器；信号发生器的信号输出端、地线端与体外线圈的两个接线端分别相连；所述的体内线圈采用生物金属材料；所述的生物金属材料为钛合金。
3.根据权利要求1所述的一种基于植入线圈的肿瘤热疗装置，其特征在于：所述体内线圈的半径r2等于5mm；无线传输距离d等于30mm。
4.如权利要求1所述的一种基于植入线圈的肿瘤热疗装置的参数优化方法，其特征在于：步骤一、根据待植入肿瘤的直径确定体内线圈的半径r2，使得体内线圈能够完全没入待植入肿瘤内，使得r2等于待植入肿瘤直径的1/5根据待植入肿瘤的位置，确定无线传输距离d，使得d等于待植入肿瘤几何中心到人体外表面的最小距离加上5mm所得值；求取体外线圈半径
Figure FDA0002613823620000014
步骤二、计算rw1＝s·r1；其中
Figure FDA0002613823620000015
步骤三、联立Q2＝1与
Figure FDA0002613823620000021
其中，μ0为真空中磁导率；f为装置工作频率，取值为1MHz、2.5MHz、5MHz或10MHz；σ2为体内线圈的电导率；得到体内线圈的线半径rw2。
5.根据权利要求4所述的一种基于植入线圈的肿瘤热疗装置的参数优化方法，其特征在于：步骤三中，f取f1、f2、…、fn；f1、f2、…、fn均小于20MHz；分别进行体内线圈的线半径rw2的计算；得到n个体内线圈的线半径候选值rw2(1)、rw2(2)、…、rw2(n)；之后执行以下步骤；
3-1.i＝1,2,…,n依次执行步骤3-2至3-6；
3-2.计算体外线圈的电感值L1如式(1)所示、体内线圈的电感值L2如式(2)所示；
L1＝μ0r1(ln(8r1/rw1)-2) 式(1)
L2＝μ0r2(ln(8r2/rw2(i))-2) 式(2)
3-3.计算体外线圈的寄生电阻值R1如式(3)所示、体内线圈的寄生电阻值R2如式(4)所示；
R1＝r1/σ1δ1rw1 式(3)
R2＝r2/σ2δ2rw2(i) 式(4)
式(3)和式(4)中，δ1为体外线圈的趋肤深度，其表达式为
Figure FDA0002613823620000022
δ2为体内线圈的趋肤深度，其表达式为
Figure FDA0002613823620000023
其中，σ1为体外线圈的电导率；σ2为体内线圈的电导率；
3-4建立体外线圈与体内线圈之间的耦合系数为k与体外线圈的半径r1、线半径rw1以及体内线圈的半径r2、线半径rw2的关系式如式(5)所示；
Figure FDA0002613823620000031
式(5)中，d为无线传输距离；c1的表达式为c1＝ln(8r1/rw1)-2；c2的表达式为c2＝ln(8r2/rw2(i))-2；
3-5、建立体外线圈的品质因数Q1的表达式如式(6)所示、体内线圈的品质因数Q2的表达式如式(7)所示
Q1＝ωL1/Rx1 式(6)
Q2＝ωL2/R2 式(7)
式(6)和式(7)中，ω＝2πf；Rx1为体外线圈的总和阻抗；
3-6.建立耦合系数为k、体外线圈的品质因数Q1、体内线圈的品质因数Q2与体内线圈、体外线圈之间电能传输效率ηi之间的关系式如式(9)所示；
Figure FDA0002613823620000032
3-7.取η1、η2、…、ηn中的最大值对应的那个体内线圈的线半径候选值作为最终的体内线圈的线半径rw2。
6.根据权利要求5所述的一种基于植入线圈的肿瘤热疗装置的参数优化方法，其特征在于：n＝4，且f1＝1MHz，f2＝2.5MHz，f3＝5MHz，f4＝10MHz。
7.根据权利要求5所述的一种基于植入线圈的肿瘤热疗装置的参数优化方法，其特征在于：Rx1＝R1。
8.根据权利要求6所述的一种基于植入线圈的肿瘤热疗装置的参数优化方法，其特征在于：Rx1＝R1+ΔZ1；其中ΔZ1的表达式如下：
Figure FDA0002613823620000033
其中，j为复数符号；J1(k)是以k为变量的1阶第一类贝塞尔函数；
Figure FDA0002613823620000041
是以
Figure FDA0002613823620000042
为变量的1阶第一类贝塞尔函数；k′2的表达式为k′2＝ωμ0σtr1 2；σt是人体组织的电导率；装置工作频率为1MHZ时，σt＝0.50268；装置工作频率为2.5MHZ时，σt＝0.55928；装置工作频率为5MHZ时，σt＝0.59008；装置工作频率为100MHZ时，σt＝0.61683。
CN201811123473.5A 2018-09-26 2018-09-26 一种基于植入线圈的肿瘤热疗装置及其参数优化方法 Active CN109363826B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811123473.5A CN109363826B (zh) 2018-09-26 2018-09-26 一种基于植入线圈的肿瘤热疗装置及其参数优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811123473.5A CN109363826B (zh) 2018-09-26 2018-09-26 一种基于植入线圈的肿瘤热疗装置及其参数优化方法

Publications (2)

Publication Number Publication Date
CN109363826A CN109363826A (zh) 2019-02-22
CN109363826B true CN109363826B (zh) 2020-10-16

Family

ID=65401959

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811123473.5A Active CN109363826B (zh) 2018-09-26 2018-09-26 一种基于植入线圈的肿瘤热疗装置及其参数优化方法

Country Status (1)

Country Link
CN (1) CN109363826B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110354392A (zh) * 2019-08-02 2019-10-22 北京大学第三医院（北京大学第三临床医学院） 一种用于热疗的线圈机构、发热装置及热疗设备
CN111398724B (zh) * 2020-04-20 2021-12-21 杭州电子科技大学温州研究院有限公司 一种无线供电用植入式集成线圈的检测方法
CN111759460B (zh) * 2020-07-07 2021-05-11 杭州电子科技大学温州研究院有限公司 基于微型植入式自谐振线圈的热疗系统的参数优化方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4740747A (en) * 1983-06-15 1988-04-26 Nippon Steel Corporation Method of and apparatus for measuring transformation degree
CN101682216A (zh) * 2007-03-27 2010-03-24 麻省理工学院 无线能量传输
CN102333567A (zh) * 2009-02-27 2012-01-25 皇家飞利浦电子股份有限公司 利用磁性纳米颗粒处置对象的治疗设备
CN104010693A (zh) * 2011-08-26 2014-08-27 海岩群落生物系统有限公司 用于产生用于治疗体腔或类似体腔的器官内的癌症的能量场的装置
CN107546866A (zh) * 2017-08-04 2018-01-05 河南师范大学 正向并联线圈电磁谐振能量传输系统设计方法
CN207184168U (zh) * 2017-08-23 2018-04-03 池州睿成微电子有限公司 一种用于传感器多节点的能量供应系统
CN108110908A (zh) * 2018-01-23 2018-06-01 中国矿业大学(北京) 非对称线圈磁耦合谐振无线电能传输系统及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6850804B2 (en) * 2002-01-18 2005-02-01 Calfacior Corporation System method and apparatus for localized heating of tissue

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4740747A (en) * 1983-06-15 1988-04-26 Nippon Steel Corporation Method of and apparatus for measuring transformation degree
CN101682216A (zh) * 2007-03-27 2010-03-24 麻省理工学院 无线能量传输
CN102333567A (zh) * 2009-02-27 2012-01-25 皇家飞利浦电子股份有限公司 利用磁性纳米颗粒处置对象的治疗设备
CN104010693A (zh) * 2011-08-26 2014-08-27 海岩群落生物系统有限公司 用于产生用于治疗体腔或类似体腔的器官内的癌症的能量场的装置
CN107546866A (zh) * 2017-08-04 2018-01-05 河南师范大学 正向并联线圈电磁谐振能量传输系统设计方法
CN207184168U (zh) * 2017-08-23 2018-04-03 池州睿成微电子有限公司 一种用于传感器多节点的能量供应系统
CN108110908A (zh) * 2018-01-23 2018-06-01 中国矿业大学(北京) 非对称线圈磁耦合谐振无线电能传输系统及方法

Also Published As

Publication number Publication date
CN109363826A (zh) 2019-02-22

Similar Documents

Publication Publication Date Title
CN109363826B (zh) 一种基于植入线圈的肿瘤热疗装置及其参数优化方法
Tang et al. Intermediate range wireless power transfer with segmented coil transmitters for implantable heart pumps
Xiao et al. An LCC-C compensated wireless charging system for implantable cardiac pacemakers: Theory, experiment, and safety evaluation
Yang et al. Inductor modeling in wireless links for implantable electronics
Xiao et al. Wireless charging system considering eddy current in cardiac pacemaker shell: Theoretical modeling, experiments, and safety simulations
Li et al. An automatic impedance matching method based on the feedforward-backpropagation neural network for a WPT system
Chow et al. An investigation into the use of orthogonal winding in loosely coupled link for improving power transfer efficiency under coil misalignment
CN108110908A (zh) 非对称线圈磁耦合谐振无线电能传输系统及方法
Stauffer et al. Practical induction heating coil designs for clinical hyperthermia with ferromagnetic implants
Duan et al. Rectangular coils optimization for wireless power transmission
Danilov et al. An algorithm for the computer aided design of coil couple for a misalignment tolerant biomedical inductive powering unit
Campi et al. Wireless powering of next-generation left ventricular assist devices (LVADs) without percutaneous cable driveline
CN104274902B (zh) 一种mri相容的植入式电极及其制造方法
Ko et al. A novel hybrid resonator for wireless power delivery in bio-implantable devices
Bao et al. Optimization of the coupling coefficient of the inductive link for wireless power transfer to biomedical implants
Lu et al. S‐PS resonant topology of WPT system for implantable spinal cord stimulator
CN212038613U (zh) 磁感应热疗控温装置
CN105326638A (zh) 一种低功耗可穿戴式电磁超声药灸装置
CN111759460B (zh) 基于微型植入式自谐振线圈的热疗系统的参数优化方法
Meng et al. Design considerations for ultrasonic power transmission to millimeter-sized implantable microelectronics devices
CN210992621U (zh) 一种用于热疗的线圈机构、发热装置及热疗设备
Campi et al. Feasibility study of a wireless power transfer system applied to a left ventricular assist device
Kavitha et al. A handy wireless charging prototype for bio-medical implantable devices
CN110890795B (zh) 一种基于非接触供电的肿瘤热消融装置
Jena et al. Efficient wireless power transfer system for biomedical applications

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201014

Address after: 325024 Floor 3 and 4 of Block B of Nanyang Avenue, Yaoxi Street, Longwan District, Wenzhou City, Zhejiang Province

Patentee after: HANGZHOU DIANZI UNIVERSITY WENZHOU RESEARCH INSTITUTE Co.,Ltd.

Patentee after: HANGZHOU DIANZI University

Address before: 325024 Floor 3 and 4 of Block B of Nanyang Avenue, Yaoxi Street, Longwan District, Wenzhou City, Zhejiang Province

Patentee before: HANGZHOU DIANZI UNIVERSITY WENZHOU RESEARCH INSTITUTE Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201202

Address after: Room 325 B, building 405, Nanyang Avenue, Wenzhou, Zhejiang Province

Patentee after: Wenzhou Huidian Technology Co.,Ltd.

Address before: 325024 Floor 3 and 4 of Block B of Nanyang Avenue, Yaoxi Street, Longwan District, Wenzhou City, Zhejiang Province

Patentee before: HANGZHOU DIANZI UNIVERSITY WENZHOU RESEARCH INSTITUTE Co.,Ltd.

Patentee before: HANGZHOU DIANZI University

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20190222

Assignee: Wenzhou Yueshi Technology Co.,Ltd.

Assignor: Wenzhou Huidian Technology Co.,Ltd.

Contract record no.: X2021330000315

Denomination of invention: A tumor hyperthermia device based on implanted coil and its parameter optimization method

Granted publication date: 20201016

License type: Common License

Record date: 20210924