CN102008351A - Acquisition method for radio frequency ablation instrument temperature field distribution - Google Patents

Acquisition method for radio frequency ablation instrument temperature field distribution Download PDF

Info

Publication number
CN102008351A
CN102008351A CN 201010520109 CN201010520109A CN102008351A CN 102008351 A CN102008351 A CN 102008351A CN 201010520109 CN201010520109 CN 201010520109 CN 201010520109 A CN201010520109 A CN 201010520109A CN 102008351 A CN102008351 A CN 102008351A
Authority
CN
China
Prior art keywords
temperature
damage
damaging
central
central point
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
Application number
CN 201010520109
Other languages
Chinese (zh)
Other versions
CN102008351B (en
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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201010520109XA priority Critical patent/CN102008351B/en
Publication of CN102008351A publication Critical patent/CN102008351A/en
Application granted granted Critical
Publication of CN102008351B publication Critical patent/CN102008351B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Surgical Instruments (AREA)

Abstract

The invention discloses an acquisition method for radio frequency ablation instrument temperature field distribution, belonging to the technical field of biological tissue heat conduction. In the method, the classical biological tissue heat conduction equation (Pennes equation) is optimized to a more concrete and higher-usability temperature field distribution model. By adopting the temperature field distribution model, the transient temperature in any place of the biological tissue during the ablation process can be conveniently, quickly and accurately obtained; and the transient temperature distribution has significant meaning on calculating the size of the valid lesion during the biological tissue ablation process.

Description

The acquisition methods that a kind of radio-frequency ablation instrument temperature field distributes
Technical field:
The present invention relates to a kind of acquisition methods that distributes at the temperature field of radio-frequency ablation instrument of optimization, especially utilize fitting formula to draw the method that the temperature field distributes in the radio-frequency ablation biological tissue process.
Background technology:
Radio-frequency ablation is to utilize the hot condenser of radio frequency that tissue is damaged, its principle is to produce high frequency electric by radio-frequency signal generator, electric current arrives by lead and damages electrode, and link to each other with radio-frequency signal generator by the reference electrode that invests body surface, form the radio-frequency current loop, produce big electric current at the tip of radio-frequency electrode, thereby being heated up, tissue destroys neurocyte and fiber, simultaneously, the electrode of organizing and make conversely that is heated up heats up, the temperature of electrode tip exposed end can detect by radio-frequency signal generator, with this temperature as near tissue.Damage the size and thickness, the length of implanting the damage electrode of kitchen range, the organizational parameter of damaging the electrode effect is relevant.
The radio system that the present invention relates to is to originate from the neural radiofrequency generator of Leksell that Sweden's medical courses in general reach company, is applied to the holonomic system of nerve stimulation, stereotaxis damage, pain therapy and bipolar coagulation.Leksell can carry out continuous monitoring to various parameters in the therapeutic process, shows that Current Temperatures changes and rate of change.Undergo surgery when damaging experiment in the treatment parkinson disease, working electrode is inserted target spot, neutral electrode is affixed on body surface, configures the radio frequency temperature and radio frequency was damaged after the persistent period.Neural radiofrequency generator continues to send altofrequency radio electric current, radio-frequency current produces the electric field of an altofrequency alternate between neutral electrode and working electrode, because working electrode end electric lines of force high concentration, density are big, the alternate electric field is enough to make hydrone along with the electric field of alternate produces the higher-order of oscillation and frictional heat, thereby this portion of tissue temperature is raise; Allowing cell produce the hot difference of coagulating, lose biological activity or the variation of generation physicochemical properties at different temperature values changes.
Can be during radio-frequency ablation in the inner temperature field that forms of the biological tissue of treatment target area with gradient, be characterized in, the temperature that the near more position of distance damage electrode obtains is high more, electrode position temperature far away more is low more apart from damaging, and distance damage electrode unlimited distance temperature is biological tissue's normal body temperature.For example, when damaging cerebral tissue, then the body temperature of distance damage electrode unlimited distance temperature behaviour is 37 ℃.What the temperature field distributed obtains, and is the successful key of radio-frequency ablation operation.
At present, main utilization is classical biological heat diffusivity equation---Pennes equation in the distribution of the temperature field of radio-frequency ablation instrument is obtained, because the Pennes equation has certain limitation in actual applications, for example amount of calculation is bigger, in embedded system, be difficult to real-time implementation etc., so can't obtain the problem of the inner transient temperature field distribution of biological tissue quickly and easily in real time.Therefore, with the temperature field distributed model that the Pennes equation is optimized for more specifically, availability is higher, by the convenient, fast and accurate transient temperature that draws place, the arbitrary locus of biological tissue in the damage process of this temperature field distributed model, this transient temperature distributes and has important directive significance for calculating radio-frequency ablation instrument effective size of damaging kitchen range in biological tissue's damage process.
Summary of the invention
Technical problem to be solved of the present invention is at the problem that can't obtain the inner transient temperature field distribution of biological tissue in the present radio-frequency ablation instrument practical work process quickly and easily in real time, proposes the acquisition methods that a kind of radio-frequency ablation instrument treatment temperature field distributes.
The present invention adopts following technical scheme for solving the problems of the technologies described above:
The acquisition methods that a kind of radio-frequency ablation instrument treatment temperature field distributes may further comprise the steps:
Step 1 on organism equation of heat conduction basis, utilizes limited element analysis technique emulation to draw in the radio-frequency ablation process, in given damage central temperature Be respectively 60 ℃, 65 ℃, 70 ℃, 75 ℃, 80 ℃, 85 ℃, 90 ℃ and damage persistent period 7 * 4 groups of biological tissue temperature fields that are respectively under 30s, 60s, 90s, the 120s distribute;
Step 2, extraction step one described 7 * 4 groups of Temperature Distribution fields are apart from damaging the central point distance
Figure 902294DEST_PATH_IMAGE004
Be respectively the temperature data at 0.5mm, 1mm, 1.5mm, 2mm, 2.5mm, 3mm, 3.5mm, 4mm, 4.5mm, 5mm place, add up different given damage central temperatures Following Temperature Distribution field and damage central point distance
Figure 382002DEST_PATH_IMAGE004
The temperature at place
Figure DEST_PATH_IMAGE005
, set up following relational expression:
Figure 643219DEST_PATH_IMAGE006
(1)
In the formula Be the time that continues to damage;
Step 3 adopts match software, respectively each is damaged central temperature
Figure 430096DEST_PATH_IMAGE001
Lasting damage persistent period when constant
Figure 246742DEST_PATH_IMAGE003
, the Temperature Distribution field with damage the central point distance
Figure 51930DEST_PATH_IMAGE004
With
Figure DEST_PATH_IMAGE007
Substitution relational expression (1) draws following relational expression:
Figure 637632DEST_PATH_IMAGE008
(2)
Wherein, coefficient
Figure DEST_PATH_IMAGE009
,
Figure 192110DEST_PATH_IMAGE010
, Be respectively about damaging central temperature Function; , , With the damage central temperature
Figure 728003DEST_PATH_IMAGE001
Functional relationship draw by above-mentioned match software, relational expression is respectively as follows:
Figure 869134DEST_PATH_IMAGE012
(5)
Figure DEST_PATH_IMAGE013
(6)
Figure 491745DEST_PATH_IMAGE014
(7)
Step 4 is updated to relational expression (2) in the relational expression (1), obtains following relational expression:
Figure DEST_PATH_IMAGE015
(3)
With formula (5), (6), (7) substitution formula (3), obtain following relational expression respectively:
Figure 60130DEST_PATH_IMAGE016
(8)
Obtaining the damage central temperature by formula (8) is
Figure 338665DEST_PATH_IMAGE001
The time, arbitrary time
Figure 221170DEST_PATH_IMAGE003
Following distance is damaged the central point distance
Figure 635971DEST_PATH_IMAGE004
The temperature at place;
Step 5 according to formula (3), is damaged central temperature at difference , according to damaging the persistent period
Figure 78770DEST_PATH_IMAGE003
With the Temperature Distribution field apart from damaging the central point distance
Figure 499388DEST_PATH_IMAGE004
Span, try to achieve the temperature in pairing temperature field, obtain the distribution in temperature field.
Further, the damage persistent period of the acquisition methods of radio-frequency ablation instrument treatment temperature field distribution of the present invention
Figure 51852DEST_PATH_IMAGE003
With the Temperature Distribution field apart from damaging the central point distance
Figure 962039DEST_PATH_IMAGE004
Span as follows:
A,, the temperature of damaging central point damages central temperature when reaching to set
Figure 887270DEST_PATH_IMAGE001
The time, the selected damage persistent period Should be less than or equal to effective damage persistent period
Figure DEST_PATH_IMAGE017
, and Temperature Distribution field and the distance of damaging central point
Figure 563288DEST_PATH_IMAGE004
Should be less than or equal in the Temperature Distribution field and the maximum distance of damaging central point , wherein effectively damage the persistent period
Figure 322482DEST_PATH_IMAGE017
, in the Temperature Distribution field with the maximum distance of damaging electrode
Figure 147219DEST_PATH_IMAGE018
Solution formula respectively as follows:
Figure DEST_PATH_IMAGE019
?; (12)
; (11)
B,, the temperature of damaging central point do not reach given damage central temperature when being in ascent stage The time, the selected damage persistent period
Figure 808510DEST_PATH_IMAGE003
Should be less than or equal to described in the A step and effectively damage the persistent period
Figure 436938DEST_PATH_IMAGE017
, and Temperature Distribution field and the distance of damaging central point
Figure 473027DEST_PATH_IMAGE004
Should be less than or equal in the Temperature Distribution field and the coverage of damaging central point
Figure DEST_PATH_IMAGE021
,
Figure 161497DEST_PATH_IMAGE021
With
Figure 610933DEST_PATH_IMAGE003
Relational expression be:
Figure 715155DEST_PATH_IMAGE022
?; (9)
Further, the organism equation of heat conduction is the Pennes equation described in the step 1 of the acquisition methods of radio-frequency ablation instrument treatment temperature field distribution of the present invention.
Further, match software is 1stOpt match software described in the step 3 of the acquisition methods of radio-frequency ablation instrument treatment temperature field distribution of the present invention.
Further, limited element analysis technique adopts finite element analysis software Femlab described in the step 1 of the acquisition methods of radio-frequency ablation instrument treatment temperature field distribution of the present invention.
The present invention adopts technique scheme to have following beneficial effect:
1, real-time, accuracy is high.The temperature field distribution acquiring method at radio-frequency ablation instrument of a kind of optimization that proposes among the present invention can obtain radio-frequency ablation instrument in real time and damage the transient temperature that object is pointed out arbitrarily , obtaining in real time for calculating biological tissue of this temperature damages that stove is long-pending to have an important directive significance.
2, workable.The method that the present invention proposes is compared the easier realization in the radio-frequency ablation instrument work process of classical Pennes equation, can be easily uses in system such as embedded.
Description of drawings:
Fig. 1 is phantom of the present invention and damages the central temperature ascending curve.
Eletrode tip was in the damage situation of different time points when Fig. 2 was 70 ℃.
Fig. 3 is the temperature curve that extracts in the emulation.
Fig. 4 is the time dependent curve of temperature.
Fig. 5 is greater than 45 ℃ of temperature curves.
Fig. 6 is
Figure 995964DEST_PATH_IMAGE001
Temperature curve in the time of=65 ℃.
Specific embodiments:
Be described in further detail below in conjunction with the enforcement of accompanying drawing technical scheme:
Step 1, each parameter of rat cerebral tissue related among the present invention is as shown in table 1.
Figure 793019DEST_PATH_IMAGE024
Table 1
Finite element analysis software Femlab is selected in the emulation of radio-frequency ablation process for use, draws biological tissue's temperature field distribution by setting up geometric model, definition physical parameter, division finite element grid and finding the solution.The present invention utilizes Femlab software that real rat cerebral tissue is simplified and is assumed to be uniform, an infinitely-great spheroid, and the radio-frequency ablation eletrode tip is positioned at ball centre, as accompanying drawing 1(1) shown in, horizontal vertical coordinate unit is m among the figure, color bar unit is Kelvin.Bulb tissue's diameter is 20mm among the figure, and ball centre is positioned at (0,0) to be located, probe tip be positioned at (0.5mm ,-1mm) locate, probe diameter is 1mm, exposed length is 2mm.The ball centre rectangle part is the electrode exposed part, its temperature curve is as accompanying drawing 1(2) shown in, electrode is under the temperature rise of 5 ℃/s, from room temperature undertissue temperature (37 ℃), linearity rises to radio-frequency ablation central temperature (70 ℃), and it is constant to remain on radio-frequency ablation central temperature point then.Classical Pennes equation (4) is adopted in emulation,
Figure DEST_PATH_IMAGE025
(4)
In the formula, T, ρ, cWith kBe respectively temperature, density, specific heat and the pyroconductivity of biological tissue; W b Be the blood perfusion rate of biological tissue, C b Specific heat for blood; T a Be the arterial blood temperature; q m Metabolic heat production rate for tissue; q r Be outside heat supply heat flow.Under the hypothesis of relative idealization, this paper has ignored tissue metabolism's heat production rate in research process q m With outside heat supply heat flow q r , and to choose rat cerebral tissue's temperature be 37 ℃ of normal body temperatures.
Choose 60 ℃, 65 ℃, 70 ℃, 75 ℃, 80 ℃, 85 ℃, 90 ℃ seven kinds different damage central temperatures respectively and carry out 7 * 4 groups of emulation with 30s, 60s, four kinds of different damage persistent period of 90s, 120s, four kinds of different design sketchs of damaging the persistent period when now being 70 ℃ with the damage central temperature are that example describes, and see accompanying drawing 2.Probe is when heating 30s down for 70 ℃ among Fig. 2 (a), the surrounding tissue temperature field presents a less ellipsoid shape, along with the prolongation of time, enlarging gradually by Fig. 2 (b), (c), (d) three width of cloth figure semielliptical shape that the surrounding tissue temperature field forms of popping one's head in as can be seen.It can be said that brightly under the certain situation of temperature, the volume of damaging kitchen range increases along with the prolongation of time, and this conforms to practical situation.
Step 2 is damaged electrode surrounding tissue temperature in the statistic procedure one, shown in the straight line of Fig. 1 ball centre position, extract and respectively organize among the simulated effect figure (5cm, 0) point to the temperature data between (5cm, 0) point.
The temperature curve of transverse diameter direction when 7 curves are 60 ℃, 65 ℃, 70 ℃, 75 ℃, 80 ℃, 85 ℃, 90 ℃ respectively for damaging central temperature in the accompanying drawing 3.Add up different given damage central temperatures
Figure 464171DEST_PATH_IMAGE001
Down, distance damage electrode position is in the temperature field
Figure 699981DEST_PATH_IMAGE004
The temperature at place
Figure 710662DEST_PATH_IMAGE023
, set up following relational expression:
Figure 740935DEST_PATH_IMAGE008
(1)
In the formula
Figure 905243DEST_PATH_IMAGE003
For damaging the persistent period;
Irreversible change all takes place when cerebral tissue is damaged during the operation of clinical radio-frequency ablation, in addition according to the situation of clinical practice, during radio-frequency ablation temperature all greater than 45 ℃, so mainly to data that the irreversible change group takes place ( 〉=45 ℃) analyze.Now choose and damage central temperature when being 60 ℃, one group of data of temperature gradient distribution describe under each time point, and data are as shown in table 2.
Figure DEST_PATH_IMAGE027
In the table 2, damaging central temperature is 60 ℃, has extracted the damage persistent period
Figure 544352DEST_PATH_IMAGE003
Be respectively 30s, 60s, 90s, 120s, distance is damaged the central point distance
Figure 745526DEST_PATH_IMAGE004
Be respectively the temperature data at 0.5mm, 1mm, 1.5mm, 2mm, 2.5mm, 3mm, 3.5mm, 4mm, 4.5mm, 5mm place.
Match software 1stOpt is in nonlinear regression, and it is very extensive that application such as the field of grade is found the solution in curve fitting, non-linear complex engineering model parameter estimation.Utilize this software, can draw under this damage central temperature according to this group temperature data, each point place temperature with
Figure 391271DEST_PATH_IMAGE003
With
Figure 640987DEST_PATH_IMAGE004
Relation such as formula (2).Wherein, coefficient ,
Figure 998336DEST_PATH_IMAGE010
,
Figure 131377DEST_PATH_IMAGE011
Be to damage central temperature
Figure 184784DEST_PATH_IMAGE001
Function.
Figure 86881DEST_PATH_IMAGE008
(2)
Step 3, emulation have selected for use 7 kinds to damage central temperature , be respectively 60 ℃, 65 ℃, 70 ℃, 75 ℃, 80 ℃, 85 ℃, 90 ℃, so can draw seven in the step 2 and have different parameters
Figure 922298DEST_PATH_IMAGE007
With
Figure 841713DEST_PATH_IMAGE004
With the damage persistent period
Figure 535999DEST_PATH_IMAGE003
Relation, the mathematical model of relational expression is formula (2).Coefficient in the relational expression ,
Figure 295194DEST_PATH_IMAGE010
,
Figure 752720DEST_PATH_IMAGE011
Analysis and statistics see Table 3.
Figure 301513DEST_PATH_IMAGE028
Step 4 utilizes match software 1stOpt to draw parameter by step 3
Figure 451872DEST_PATH_IMAGE009
,
Figure 781222DEST_PATH_IMAGE010
,
Figure 714543DEST_PATH_IMAGE011
With the damage central temperature
Figure 445738DEST_PATH_IMAGE001
Relation be respectively:
Figure 704682DEST_PATH_IMAGE012
(5)
Figure 255749DEST_PATH_IMAGE013
(6)
Figure 320657DEST_PATH_IMAGE014
(7)
Formula (2) is updated in the formula (1), obtains formula (3):
Figure 844042DEST_PATH_IMAGE015
(3)
With formula (5), (6), (7) substitution formula (3), obtain:
Figure 64764DEST_PATH_IMAGE016
(8)
Like this, can obtain the damage central temperature is The time, arbitrary time
Figure 909410DEST_PATH_IMAGE003
Following distance is damaged central point
Figure 349618DEST_PATH_IMAGE004
The temperature of place's point.The present invention has carried out the significance checking in fit procedure, significance level is a critical probability value.It is illustrated in " assumed statistical inspection ", infers with sample data and makes the probability size that refusal " is supposed " mistake when overall.Significance level is more little, and it is more little to make the wrong probability that refusal " supposes ".For formula (8), significance level
Figure DEST_PATH_IMAGE029
=0, therefore above-mentioned equation is set up.
The transmission in vivo of step 5, temperature is dynamic, and under given damage central temperature, the temperature of biological tissue can raise along with the prolongation of damaging the persistent period.Be example now with accompanying drawing 4, given damage central temperature among the figure
Figure 12681DEST_PATH_IMAGE001
Be 70 ℃, the situation of change of the biological curve of biological tissue when four curves are respectively and damage duration and be 30s, 60s, 90s, 120s, transverse axis is the distance (the damage central point is at 0 place, and unit is mm) at biological tissue and center, the longitudinal axis is a temperature.
At among the figure with damage the biological tissue of central point distance for a place, when the damage persistent period was respectively 30s, 60s, 90s, 120s, the tissue temperature at this place was respectively
Figure 476023DEST_PATH_IMAGE030
, ,
Figure 945050DEST_PATH_IMAGE032
, , and
Figure 352264DEST_PATH_IMAGE032
Figure 234769DEST_PATH_IMAGE033
, analyze as can be known, when the time infinitely increased, temperature can reach balance in this biological tissue's body, and promptly the arbitrfary point temperature tends towards stability in the temperature field, no longer rises.Since in the invention main by research temperature field inner tissue temperature greater than 45 ℃ zone to judge the radio-frequency ablation situation, the data of other temperature provinces are meaningless to the present invention, therefore, mainly in this zone
Figure 649570DEST_PATH_IMAGE004
With
Figure 326539DEST_PATH_IMAGE003
Analyze and illustrate, see accompanying drawing 4.
Formula in the step 4 (8) should divide following two kinds of situation discussion:
A kind of is temperature when not reaching balance, With
Figure 575304DEST_PATH_IMAGE034
Range specification.Suppose this moment For damaging the persistent period be
Figure 958061DEST_PATH_IMAGE034
The time, the tissue temperature that makes that the damage central temperature can be transmitted reaches 45 ℃ farthest.For example, when damaging 30s among Fig. 5, be organized in
Figure DEST_PATH_IMAGE035
The place can reach 45 ℃ of critical temperatures, when damaging 120s, is organized in
Figure 945608DEST_PATH_IMAGE036
The place reaches 45 ℃ of critical temperatures.The damage central temperature of this moment
Figure 169916DEST_PATH_IMAGE001
,
Figure 559309DEST_PATH_IMAGE021
With
Figure 905977DEST_PATH_IMAGE003
Between have certain relation,
Figure 52924DEST_PATH_IMAGE021
Along with
Figure 149100DEST_PATH_IMAGE003
Prolongation and increase.If the temperature field reaches stable state and needs the time
Figure DEST_PATH_IMAGE037
, then ought select
Figure 455316DEST_PATH_IMAGE003
Figure 910568DEST_PATH_IMAGE037
The time,
Figure 872708DEST_PATH_IMAGE004
Must be in corresponding scope value, differ bigger result otherwise can produce with practical situation.
Through utilizing match software 1stOpt to analyze and statistics, can draw in the dynamic temperature field
Figure 501135DEST_PATH_IMAGE001
, With
Figure 225695DEST_PATH_IMAGE003
Relational expression be:
Figure 347235DEST_PATH_IMAGE022
(9)
Promptly be when damaging central temperature ,
Figure 997845DEST_PATH_IMAGE003
Figure 794899DEST_PATH_IMAGE037
The time,
Figure 466052DEST_PATH_IMAGE037
Value can obtain with reference to the derivation of formula (12), choose
Figure 373965DEST_PATH_IMAGE004
Should less than
Figure 446963DEST_PATH_IMAGE021
Another kind is a temperature when reaching balance,
Figure 477236DEST_PATH_IMAGE004
With
Figure 573368DEST_PATH_IMAGE003
Range specification.No matter this moment
Figure 347289DEST_PATH_IMAGE003
How to increase,
Figure 478056DEST_PATH_IMAGE021
Value can infinitely not increase because of its increase yet.Promptly Can infinitely increase, and
Figure 59396DEST_PATH_IMAGE021
Can only be in limited scope value.From simulation result, for damaging central temperature
Figure 574691DEST_PATH_IMAGE001
=70 ℃, when damaging the persistent period during greater than 120s, the temperature field has reached balance substantially.Therefore, can think ought
Figure 622282DEST_PATH_IMAGE001
=70 ℃,
Figure 932040DEST_PATH_IMAGE003
Data during 〉=120s are nonsignificant data, promptly damage central temperature at one, and the effective time of a maximum should be arranged Corresponding with it.
Finding the solution
Figure 852909DEST_PATH_IMAGE001
With Relation before, need know and damage central temperature
Figure 315561DEST_PATH_IMAGE001
Down, temperature can obtain according to simulation result along with the variation relation of time
Figure 873581DEST_PATH_IMAGE026
With
Figure 527416DEST_PATH_IMAGE003
Relational expression be:
Figure 549599DEST_PATH_IMAGE038
(10)
Damage central temperature to seven respectively at this
Figure 201160DEST_PATH_IMAGE001
Under the formula differentiate, draw each temperature variations.It is 65 ℃ now to damage central temperature
Figure DEST_PATH_IMAGE039
The pass is that example describes, at distance center distance in the field
Figure 43214DEST_PATH_IMAGE040
Place's point temperature data statistics over time is as follows.
Figure DEST_PATH_IMAGE041
Utilize the form match of match software 1stOpt according to formula (10) Curve draws =48,
Figure 262209DEST_PATH_IMAGE010
=-42.Formula after the differentiate is:
Figure 529242DEST_PATH_IMAGE042
(11)
Accompanying drawing 6 is that the damage central temperature is 65 ℃, damages temperature variation curve in the process of persistent period from 40s to 150s.In actual applications, consider the operability of instrument, can think when variations in temperature within the 10s during smaller or equal to 1 ℃ the temperature field in a basic balance, therefore need find the solution and work as
Figure DEST_PATH_IMAGE043
The time
Figure 790459DEST_PATH_IMAGE003
Value.For formula (11), obtain Promptly (perhaps be called effective time when damaging the persistent period
Figure 515019DEST_PATH_IMAGE017
) when being 35.4s, the temperature field reaches balance substantially.Each damages the data statistics of central temperature and effective time:
Can draw according to this rule
Figure 393982DEST_PATH_IMAGE001
With
Figure 396573DEST_PATH_IMAGE017
Relational expression:
Figure 982275DEST_PATH_IMAGE019
(12)
Accordingly
Figure 146540DEST_PATH_IMAGE018
Span with
Figure 450482DEST_PATH_IMAGE001
Data statistics:
The pass is:
Figure 431394DEST_PATH_IMAGE020
(13)
Promptly be when damaging central temperature
Figure 94456DEST_PATH_IMAGE001
The time, the selected time
Figure 557799DEST_PATH_IMAGE003
Should smaller or equal to
Figure 970369DEST_PATH_IMAGE017
, and with the distance of damaging central point
Figure 202767DEST_PATH_IMAGE004
Should smaller or equal to
Figure 36731DEST_PATH_IMAGE018

Claims (5)

1. a radio-frequency ablation instrument is treated the acquisition methods that the temperature field distributes, and it is characterized in that: may further comprise the steps:
Step 1 on organism equation of heat conduction basis, utilizes limited element analysis technique emulation to draw in the radio-frequency ablation process, in given damage central temperature Be respectively 60 ℃, 65 ℃, 70 ℃, 75 ℃, 80 ℃, 85 ℃, 90 ℃ and damage persistent period
Figure DEST_PATH_IMAGE002
7 * 4 groups of biological tissue temperature fields that are respectively under 30s, 60s, 90s, the 120s distribute;
Step 2, extraction step one described 7 * 4 groups of Temperature Distribution fields are apart from damaging the central point distance
Figure 872978DEST_PATH_IMAGE003
Be respectively the temperature data at 0.5mm, 1mm, 1.5mm, 2mm, 2.5mm, 3mm, 3.5mm, 4mm, 4.5mm, 5mm place, add up different given damage central temperatures
Figure 278814DEST_PATH_IMAGE001
Following Temperature Distribution field and damage central point distance
Figure 674023DEST_PATH_IMAGE003
The temperature at place , set up following relational expression:
(1)
In the formula
Figure 186224DEST_PATH_IMAGE002
For damaging the lasting time;
Step 3 adopts match software, respectively each is damaged central temperature
Figure 79356DEST_PATH_IMAGE001
Damage persistent period when constant
Figure 402890DEST_PATH_IMAGE002
, the Temperature Distribution field with damage the central point distance
Figure 71768DEST_PATH_IMAGE003
With
Figure DEST_PATH_IMAGE006
Substitution relational expression (1) draws following relational expression:
Figure 815865DEST_PATH_IMAGE007
(2)
Wherein, coefficient ,
Figure 297489DEST_PATH_IMAGE009
,
Figure DEST_PATH_IMAGE010
Be respectively about damaging central temperature
Figure 785233DEST_PATH_IMAGE001
Function;
Figure 308618DEST_PATH_IMAGE008
,
Figure 364561DEST_PATH_IMAGE009
,
Figure 730821DEST_PATH_IMAGE010
With the damage central temperature
Figure 632043DEST_PATH_IMAGE001
Functional relationship draw by above-mentioned match software, relational expression is respectively as follows:
(5)
Figure DEST_PATH_IMAGE012
(6)
Figure 361413DEST_PATH_IMAGE013
(7)
Step 4 is updated to relational expression (2) in the relational expression (1), obtains following relational expression:
Figure DEST_PATH_IMAGE014
(3)
With formula (5), (6), (7) substitution formula (3), obtain following relational expression respectively:
Figure 841067DEST_PATH_IMAGE015
(8)
Obtaining the damage central temperature by formula (8) is
Figure 44515DEST_PATH_IMAGE001
The time, arbitrary time
Figure 168591DEST_PATH_IMAGE002
Following distance is damaged the central point distance The temperature at place;
Step 5 according to formula (3), is damaged central temperature at difference
Figure 15510DEST_PATH_IMAGE001
, according to damaging the persistent period With the Temperature Distribution field apart from damaging the central point distance
Figure 1232DEST_PATH_IMAGE003
Span, try to achieve the temperature in pairing temperature field, obtain the distribution in temperature field.
2. the acquisition methods that radio-frequency ablation instrument treatment according to claim 1 temperature field distributes is characterized in that: the described damage persistent period With the Temperature Distribution field apart from damaging the central point distance
Figure 132447DEST_PATH_IMAGE003
Span as follows:
A,, the temperature of damaging central point damages central temperature when reaching to set
Figure 179163DEST_PATH_IMAGE001
The time, the selected damage persistent period
Figure 651732DEST_PATH_IMAGE002
Should be less than or equal to effective damage persistent period
Figure DEST_PATH_IMAGE016
, and Temperature Distribution field and the distance of damaging central point
Figure 515914DEST_PATH_IMAGE003
Should be less than or equal in the Temperature Distribution field and the maximum distance of damaging central point
Figure 565779DEST_PATH_IMAGE017
, wherein effectively damage the persistent period
Figure 416185DEST_PATH_IMAGE016
, in the Temperature Distribution field with the maximum distance of damaging electrode
Figure 602316DEST_PATH_IMAGE017
Solution formula respectively as follows:
Figure DEST_PATH_IMAGE018
?; (12)
Figure 512766DEST_PATH_IMAGE019
; (11)
B,, the temperature of damaging central point do not reach given damage central temperature when being in ascent stage
Figure 49926DEST_PATH_IMAGE001
The time, the selected damage persistent period
Figure 704024DEST_PATH_IMAGE002
Should be less than or equal to described in the A step and effectively damage the persistent period
Figure 885606DEST_PATH_IMAGE016
, and Temperature Distribution field and the distance of damaging central point
Figure 199913DEST_PATH_IMAGE003
Should be less than or equal in the Temperature Distribution field and the coverage of damaging central point , With
Figure 42677DEST_PATH_IMAGE002
Relational expression be:
?; (9)。
3. the acquisition methods that radio-frequency ablation instrument treatment according to claim 1 temperature field distributes, it is characterized in that: the equation of heat conduction of organism described in the step 1 is the Pennes equation.
4. the acquisition methods that radio-frequency ablation instrument treatment according to claim 1 temperature field distributes, it is characterized in that: the software of match described in the step 3 is 1stOpt match software.
5. the acquisition methods that radio-frequency ablation instrument treatment according to claim 1 temperature field distributes, it is characterized in that: limited element analysis technique described in the step 1 adopts finite element analysis software Femlab.
CN201010520109XA 2010-10-27 2010-10-27 Acquisition method for radio frequency ablation instrument temperature field distribution Expired - Fee Related CN102008351B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010520109XA CN102008351B (en) 2010-10-27 2010-10-27 Acquisition method for radio frequency ablation instrument temperature field distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010520109XA CN102008351B (en) 2010-10-27 2010-10-27 Acquisition method for radio frequency ablation instrument temperature field distribution

Publications (2)

Publication Number Publication Date
CN102008351A true CN102008351A (en) 2011-04-13
CN102008351B CN102008351B (en) 2012-05-30

Family

ID=43838832

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010520109XA Expired - Fee Related CN102008351B (en) 2010-10-27 2010-10-27 Acquisition method for radio frequency ablation instrument temperature field distribution

Country Status (1)

Country Link
CN (1) CN102008351B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107049476A (en) * 2017-04-20 2017-08-18 北京工业大学 A kind of voltage source characterizing method that feedback is surveyed based on single needle
CN107578820A (en) * 2017-08-21 2018-01-12 北京工业大学 A kind of biological properties parameter characterization method based on single needle actual measurement

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6706040B2 (en) * 2001-11-23 2004-03-16 Medlennium Technologies, Inc. Invasive therapeutic probe
CN1771882A (en) * 2005-09-02 2006-05-17 武汉市昊博科技有限公司 Method and device for obtaining internal heat source information from the surface temperature distribution of living body

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6706040B2 (en) * 2001-11-23 2004-03-16 Medlennium Technologies, Inc. Invasive therapeutic probe
CN1771882A (en) * 2005-09-02 2006-05-17 武汉市昊博科技有限公司 Method and device for obtaining internal heat source information from the surface temperature distribution of living body

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《工程热物理学报》 20040731 赵镇南等 肿瘤射频热疗SAR场与瞬时温度分布模拟 667-669 1-5 第25卷, 第4期 2 *
《江苏医药》 20100531 陶振玉等 体内外射频热凝毁损灶体积的相关性分析 1058-1061+988 1-5 第36卷, 第9期 2 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107049476A (en) * 2017-04-20 2017-08-18 北京工业大学 A kind of voltage source characterizing method that feedback is surveyed based on single needle
CN107049476B (en) * 2017-04-20 2019-07-26 北京工业大学 A kind of voltage source characterizing method based on single needle actual measurement feedback
CN107578820A (en) * 2017-08-21 2018-01-12 北京工业大学 A kind of biological properties parameter characterization method based on single needle actual measurement

Also Published As

Publication number Publication date
CN102008351B (en) 2012-05-30

Similar Documents

Publication Publication Date Title
CN108830017B (en) Radio frequency heating temperature field prediction system based on individual impedance
CN110151309B (en) Preoperative planning method and equipment for multi-modal ablation therapy
Dodde et al. Monopolar electrosurgical thermal management for minimizing tissue damage
CN107526928B (en) Microwave thermal ablation temperature field simulation method based on characteristic parameter feedback
KR20140022883A (en) Devices and methods for shaping therapy in fluid enhanced ablation
CA2733619A1 (en) Monitoring tissue temperature while using an irrigated catheter
TW200835462A (en) Thermal imaging feedback for optimizing radio frequency ablation therapy
JP2020508152A (en) System and method for ablation status monitoring and custom ablation molding
Shao et al. Studying the thermal performance of a bipolar radiofrequency ablation with an improved electrode matrix system: In vitro experiments and modelling
Watanabe et al. Temperature dependence of thermal conductivity of liver based on various experiments and a numerical simulation for RF ablation
CN102008351A (en) Acquisition method for radio frequency ablation instrument temperature field distribution
Ranjbartehrani et al. Characterization of miniature probes for cryosurgery, thermal ablation, and irreversible electroporation on small animals
CN107157577B (en) A kind of large area discharge ablating device
CN107578820A (en) A kind of biological properties parameter characterization method based on single needle actual measurement
Ulucakli Simulation of radiofrequency ablation and thermal damage to tissue
Garcia et al. An experimental investigation of temperature changes during electroporation
Watanabe et al. Validation of accuracy of liver model with temperature-dependent thermal conductivity by comparing the simulation and in vitro RF ablation experiment
Isobe et al. Real-time temperature control system based on the finite element method for liver radiofrequency ablation: Effect of the time interval on control
CN207041577U (en) A kind of tumour radio-frequency ablation electrode
Watanabe et al. Modeling the temperature dependence of thermophysical properties: Study on the effect of temperature dependence for RFA
TW201721144A (en) Measuring device for thermal physical parameter and needle
CN205814414U (en) Piercing assembly and the electric fire needle minimally-invasive treatment system with this piercing assembly
Yamazaki et al. Modeling the internal pressure dependence of thermal conductivity and in vitro temperature measurement for lung RFA
CN215306633U (en) Integrative radio frequency electrode of temperature measurement radio frequency
Lu et al. Method for estimating the temperature distribution associated with the vessel cooling effect in radio frequency ablation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120530

CF01 Termination of patent right due to non-payment of annual fee