Summary of the invention:
The present invention overcomes in the prior art combination that two above microwave intervention treatment antennas produce and solidifies the technological deficiency that the real-time three-dimensional of thermal field distributes and can't learn, a kind of thermal field combined method of many active antennas is provided, and the real-time three-dimensional that can obtain various combination thermal fields distributes.
The present invention is achieved in that 1, analogue measurement and calculate the real-time three-dimensional heat of solidification field distribution of single active antenna under different radiant power effects by the composite equation formula; 2, the real-time three-dimensional thermal field with single antenna is distributed as the basis, two above antennas are carried out permutation and combination by following six kinds of conditions, set up computer model, utilize finite element analysis software to obtain three-dimensional heat of solidification field distribution data base under the microwave different condition, six kinds of conditions are: A), antenna type condition B) the antenna amount condition C) antenna space locality condition D) power condition .E), the time conditions F of different capacity effect), the blood perfusion rate condition of cancerous issue; 3, will treat required temperature boundary condition carries out adding behind the three-dimensional reconstruction among the three-dimensional heat of solidification field distribution thermal field data base under the described different condition, and in computer, form effective microwave thermal field data storehouse of satisfying the said temperature boundary condition under six kinds of conditions, what this data base comprised the real-time distribution data of solidifying thermal field and each data correspondence solidifies thermal field graphics spectrum.
Above-mentioned measurement can be an analogue measurement, promptly sets up the corresponding microwave radiation experimental provision of organizing phantom, measures the variations in temperature of different capacity under the microwave field, time and position, and with experimental data substitution composite equation formula; Described composite equation formula is meant the distribution function equation of biological heat transfer equation and microwave radiation energy SAR.Described biological heat transfer equation is:
ρc dT/dt=SAR+KΔT-M
bC
b(T-T
b)
Wherein: ρ is the (kg/m of tissue density
3), c is a specific heat, and T is a temperature, and t is the time, and K is the heat conductivity (W/m ℃) of tissue, M
b, C
bBe respectively the filling rate (kg/m of blood
3S) and specific heat (J/kg ℃), T
bFor this regional arterial blood temperature (℃).
The distribution function equation of described microwave radiation energy SAR is:
SAR=ρce
-αr[x
0z
3+x
1z
2+x
2z+x
3][cos(x
5θ)+x
4]
Wherein, r, θ, z are coordinate variable, α, X
0, X
1, X
2, X
3, X
4, X
5Be coefficient.
Learn that through preliminary experiment the SAR distribution function is the plane distribution function:
SAR=ρce
-αr[X
0z
3+X
1z
2+X
2z+X
3]
It is that the thermal field emittance SAR/ ρ c that solidifies of 60W microwave radiation partly is distributed as at the forward direction of microwave antenna central point apart from antenna pointed end 15mm that above-mentioned real-time three-dimensional thermal field distributes for power:
2.2667×2.9e
-0.3164r(-0.0002z
3+0.0042z
2-0.0344z+0.4864)
The microwave antenna central point is distributed as to part apart from the back of other end 12mm:
2.2667×2.9e
-0.3164r(-0.00003z
3+0.0013z
2-0.0306z+0.4864)
It is that the thermal field emittance SAR/ ρ c that solidifies of 50W microwave radiation partly is distributed as at the forward direction of microwave antenna central point apart from antenna pointed end 15mm that above-mentioned real-time three-dimensional thermal field distributes for power:
1.4903×3.02e
-0.3002r(8E-5z
3-0.003z
2+0.0091z+0.4152)
The microwave antenna central point is distributed as to part apart from the back of other end 12mm:
1.4903×2.98e
-0.3002r(4E-5z
3-0.0012z
2-0.0066z+0.4168)
It is that the thermal field emittance SAR/ ρ c that solidifies of 40W microwave radiation partly is distributed as at the forward direction of microwave antenna central point apart from antenna pointed end 15mm that above-mentioned real-time three-dimensional thermal field distributes for power:
1.1264×5.3e
-0.2829r(0.00003z
3-0.0016z
2+0.0042z+0.3651)
The microwave antenna central point is distributed as to part apart from the back of other end 12mm:
1.1264×5.2e
-0.2829r(0.00003z
3-0.0009z
2-0.0117z+0.3692)
The temperature range of said temperature boundary condition is 54-60 ℃.
Good effect of the present invention:
It is A that thermal field combined method of the present invention has been set up with described six kinds of conditions in computer), antenna type condition B) power condition .E antenna amount condition C) antenna space locality condition D)), the time conditions F of different capacity effect), the blood perfusion rate condition of cancerous issue is a parameter, and effective microwave thermal field data storehouse of satisfied temperature boundary condition, what this data base comprised the real-time three-dimensional data of microwave coagulation thermal field and each corresponding data correspondence solidifies thermal field graphics spectrum, effective microwave thermal field data storehouse of the three-dimensional real-time distribution of this combination thermal field provides multiple alternative therapeutic scheme for clinical micro-wave therapeutic tumor, and can be in clinical treatment select therapeutic regimen according to patient's different situations, make the normal tissue injury that produces tumor little and therapeutic effect is good.
Method of the present invention is applicable to that various tumor microwaves import the thermal field combination of treatment, and different tumors only need carry out analogue measurement according to its phantom, draw the thermal field data of single antenna, and in the combination of carrying out various conditions, the real-time three-dimensional thermal field of obtaining after the combination distributes.
The specific embodiment:
Be distributed as example in the hope of the three-dimensional thermal field in real time of many combination of antennas microwave intervention treatment of hepatocarcinoma.
At first measure and calculate the real-time three-dimensional heat of solidification field distribution of single radiation electrode under different radiant power effects by the composite equation formula;
Learn one pole implanted microwave antenna in the r direction by the Hepar Sus domestica that exsomatizes is carried out one pole microwave radiation preliminary experiment, the θ direction is symmetrical distribution, thereby SAR can be found the solution and be used as the two-dimensional problems processing.The analogue measurement device of design drawing 1 in view of the above, with thick 1.5cm, 15 * 15cm is square, flat poly (methyl methacrylate) plate center hollows out, form the framework 20 of all hem width 1.5cm, this framework 20 is just put into the D=15.1cm that is used for adorning phantom, high 20cm is as in the lucite cylinder barrel 21 of support cover, adopt lucite to be because this kind material is noiseless to the microwave thermal field, and it is transparent, visual good, two opposite side at organic glass frame 20 are symmetrical arranged the guid needle through hole 22 that a pair of diameter is 1.7mm, just put into the microwave guide needle, each makes a call to the backguy through hole 24 of 5 D<1mm up and down at two opposite side of its right side decentre 5mm, distance between borehole 5mm, pass last pull-down line through hole 24 with 1 trumpeter's art stitching thread, form 5 backguys 25 with microwave antenna parallel interval 5mm along framework 20; Described microwave guide needle guiding microwave antenna 26 is inserted in the framework 1 from guid needle through hole 22; Certain position apart from microwave antenna 26 in framework 20 is provided with measurement thermocouple 27, the diameter that thermocouple 27 adopts U.S. Ω company to produce is 0.003 inch an iron-constantan thermocouple silk, all thermocouples 27 are linked to each other with described data acquisition unit, data acquisition unit links to each other with computer, by preliminary experiment as can be known, because thermal field is along the r direction, the θ direction is symmetrical distribution, thereby SAR can be found the solution and be used as the two-dimensional problems processing, thermocouple 27 carries out the plane by Fig. 2 layouts, and Fig. 3 is the spatial value of Fig. 2 each point correspondence.And be separately fixed in the backguy 25 among Fig. 1 by this cloth point diagram, each 16 points of cloth are after cloth is better; together put into cylinder barrel 21 together with microwave antenna 26; the slow base of liver MICROWAVE PHANTOM MUSCLE TISSUE that has prepared is gone into, left standstill 12-24 hour, until phantom temperature and indoor temperature balance.Adopt Jiang-Han Area guarantor of Sichuan University to wait the MICROWAVE PHANTOM MUSCLE TISSUE prescription of development in this liver MICROWAVE PHANTOM MUSCLE TISSUE, the complex permittivity of this phantom is in the 300-2500MHz scope, with error≤5% of actual human body liver organization.Check cold junction temperature before measuring, and the temperature by the preceding measurement point of computer acquisition microwave action, with in the temperature fluctuation 5 minutes<0.01 ℃ as the stable state index.(60W300s, 50W10min 40W30min) divide above-mentioned 16 points to carry out 16-18 time and measure each energy group, and 3 energy groups are carried out 53 times altogether and measured each site repetition 5-9 time, the meansigma methods that each point is got in calculating at last.Energy gives 60W300s (every 100s is 10s at interval) respectively; 50W10min; 40W30min; described on data acquisition unit is with every 500ms (60W300s) or 1s (40W30min; 50W10min) 16 passages of the continuous scanning of speed; the computer real-time record; shutting down the back observes continuously; with temperature drop to (initial temperature 21-24 ℃) below 43 ℃ and in 1 minute variations in temperature<1 ℃ be the index of shutdown; the 60W effect is organized general 15-20 minute and is reached stable, reaches stable state in 50W10 minute general 30-35 minute, and 40W30min organizes general 35-45 minute and reaches relative stable state; deposit resultant data in computer, use to do to analyze.
In the r direction, the θ direction is symmetrical distribution by preliminary experiment notice of invitation utmost point implanted microwave antenna, and the Z direction is for being the center distribution or asymmetrical distribution apart from antenna pointed end 15mm place promptly with the microwave antenna central point.The distribution function equation simplification of SAR is: SAR=ρ Ce
-α Z[X
0r
3+ X
1r
2+ X
2R+X
3], so only need measure 7 some place variation of temperature, separate the Simultaneous Equations of forming by 7 equations, just can obtain every coefficient of distribution function.Obtain the distribution function of SAR/ ρ c, see the following form.
60W |
Forward direction |
2.2667×2.9e
-0.3164r(-0.0002z
3+0.0042z
2-0.0344z+0.4864)
|
The back to |
2.2667×2.9e
-0.3164r(-0.00003z
3+0.0013z
2-0.0306z+0.4864)
|
50W |
Forward direction |
1.4903×3.02e
-0.3002r(8E-5z
3-0.003z
2+0.0091z+0.4152)
|
The back to |
1.4903×2.98e
-0.3002r(4E-5z
3-0.0012z
2-0.0066z+0.4168)
|
40W |
Forward direction |
1.1264×5.3e
-0.2829r(0.00003z
3-0.0016z
2+0.0042z+0.3651)
|
The back to |
1.1264×5.2e
-0.2829r(0.00003z
3-0.0009z
2-0.0117z+0.3692)
|
When known density and specific heat, just can obtain the distribution function of whole microwave thermal field SAR.
Utilize U.S. ANSYS (V5.5 version) finite element analysis software to finish whole computational process and (comprise pre-treatment-----spatial spreading based on FInite Element, the derivation of finite element equation and post processing-----further calculation of parameter, data output, the drafting of form and curve chart etc.), calculate the dynamic 3 D thermal field of singly leading implanted microwave solidification treatment liver tumor.
Secondly, after trying to achieve the radiation thermal field of singly leading, can make up, try to achieve the distributed data base of combination thermal field by following six conditions and value thereof.
(1) power condition:
Adopt 40W, 50W, three kinds of microwave powers of 60W.
(2) be used for the antenna amount condition of hepatocarcinoma solidification treatment
According to the size of different tumors, applicable number of electrodes is: 1,2,3,4, be respectively singly and lead, twoly lead, three lead, four lead.
(3) hepatocarcinoma local organization blood perfusion rate (Mb) condition is divided into 10 grades:
Mb is respectively: 1.5,2.5,3.0,3.5,4.0,4.5,5.0,6.0,7.0,8.0 (units: kg/m
3S).
(4) used time conditions:
Corresponding to different capacity, adopt the different time.
40W——600、900、1200、1800s
50W——300、400、500、600、700、800、900s
60W——300、400、500、600s
(5) antenna space locality condition:
Singly lead situation slightly.
Under two leading (timesharing is two to be led and the two simultaneously leading) situation, two antenna distances adopt 1.0cm, 1.2cm, 1.5cm, 1.6cm, 2.0cm.(annotate: two the leading of timesharing all adopted 20s: the timesharing ratio of 20s, i.e. and antenna effect 20s, another root antenna remakes uses 20s).Fig. 4 has shown that two antennas 31 and 32 parallel arrangements act on the sketch map on the tumor 33.
Three lead under the situation, can be in two kinds of situation, and promptly wherein two electrodes are that timesharing double electrode or three electrodes act on simultaneously.When two electrodes wherein are that the two of same microwave source lead the timesharing electrode, promptly two electrodes can act on different time respectively, and being increases thermocoagulation efficient, and another single electrode adopts continuous action.For obtaining better therapeutic effect, adopt the classification of big freezing range: two timesharing electrodes of leading are with spacing 1.5cm, 2.0cm, make another single electrode be positioned among above two electrode connecting lines on the vertical line spacing 1.5cm or 2.0cm.
Four lead under the situation, can divide three kinds of situations, and promptly wherein having only two electrodes is that timesharing double electrode or four electrodes are respectively two timesharing double electrodes or four electrodes are effect simultaneously.When solidification mode for two when leading the effect of timesharing electrode and forming in twos, can be with four electrodes by square arrangement, the adjacent electrode spacing is 1.0cm, 1.2cm, 1.5cm, 1.6cm, 2.0cm.And the mode that adopts adjacent electrode to act on simultaneously.
(6) antenna type condition
Refer to that antenna produces difformity and solidifies the dissimilar antennas of thermal field correspondence, the antenna or the sphere of solidifying thermal field such as the generation spindle-type are solidified the thermal field antenna.
Above 6 kinds of situations are carried out permutation and combination, set up computer model, and utilize finite element analysis software to obtain microwave coagulation thermal field under the different condition.
To treat required boundary condition at last carries out adding among the thermal field data base of the three-dimensional heat of solidification field distribution under the described different condition behind the three-dimensional reconstruction, the preferred temperature boundary condition of rule of thumb learning treatment hepatocarcinoma is 54 ℃ and/or 60 ℃, the temperature boundary condition determines that the back just can form effective microwave thermal field data storehouse of satisfying the said temperature boundary condition under six kinds of conditions respectively in computer, this data base comprises the real-time distribution data of solidifying thermal field and the corresponding thermal field graphics spectrum of solidifying, the real-time distribution data of this thermal field are the temperature profile data of each point of space on different time points, just form many 54 and 60 isothermal surfaces of spending two temperature spots at one time.Like this, only need before operation, to have measured the size of tumor, obtain the space coordinates and the three-dimensional collection of illustrative plates of form at tumor edge, record its blood perfusion rate again, solidify the real-time distribution data of thermal field and corresponding known isothermal surface data compare with above-mentioned, just obtain several and the close thermal field of space coordinates its edge, and the three-dimensional collection of illustrative plates of the three-dimensional collection of illustrative plates of the form of tumor and several approximate thermal fields covered contrast, determine preferred thermal field at last, the quantity term of other conditions of this thermal field such as power condition, antenna and type and locus condition are just directly as can be known.
Further, effective microwave thermal field data storehouse that combinations thereof is tried to achieve can be divided into following plurality of data storehouse: 1) thermal field data base: the distributed in three dimensions data of following 54 degrees centigrade and the 60 degrees centigrade isothermal levels of different parameters condition, comprise corresponding thermal field number, r, θ, z coordinate, temperature, tumor blood flow filling rate, corresponding scheme number.Include the corresponding described thermal field graphics spectrum 2 of solidifying in this thermal field data base) the therapeutic scheme data base: the therapeutic scheme storehouse of thermal field distributed in three dimensions data correspondence comprises the locus of scheme number, power, time, antenna number, antenna.3) property data base: according to case statistics in the past tumor is done the form classification by certain sorting technique, size fractionation, the blood flow classification is convenient to sort out relatively contrast, comprises tumor form classification, big or small rank and blood perfusion rate rank.
For the tumor of other types, only need when analogue measurement, to select for use the phantom of corresponding tissue to carry out analogue measurement, other programs of setting up effective microwave thermal field are identical with hepatocarcinoma with step.