RF ablation thermal simulation model: Parameter sensitivity analysis
Abstract
OBJECTIVE:
The aim of the research is to obtain the relative influences of some critical electro-thermal parameters on the ablation temperature and lesion volume during temperature-controlled radiofrequency ablation (RFA) of liver tumor by parameter sensitivity analysis.
METHODS:
The finite element method (FEM) has been used to establish the simulation model of RFA temperature field, and the sensitivity of the tissue parameters has been analyzed. The effects of six parameters have been taken into account, including the thermal specific capacity (Cp), the thermal conductivity (k), the electrical conductivity (Sigma), the density (rho), the dielectric constant (Epsilon) and the resistance (R). The simulation processes based on different parameter values have been accomplished with Comsol Multiphysics software, and the sensitivity parameters have been obtained utilizing the variance contribution rate (SS%) or the main effects.
RESULTS:
It was found that the ablation temperature and lesion volume increased with increasing the values of Rand Sigma, but was a reverse situation for Cp and rho. Besides, the influence of k on ablation volume was relatively small and Epsilon had a negligible effect on ablation temperature.
CONCLUSIONS:
It is concluded that these parameter sensitivity results can provide scientific and reliable reference for the specificity analysis of the RF ablation models.
1.Introduction
RFA technology is to apply the AC electric field to the tumor tissue, thus achieving tumor coagulation necrosis by the resistance heating. This technique has been widely used in the treatment of liver tumors due to its minimal invasiveness. However, the ablation quality still depends on the experiences of clinicians. Some studies [1, 2, 3] have shown that RFA outcomes mainly rely on the lesion volume size, shape of the ablation zone and its internal temperature distribution. Thence, the precise characterization of the RFA temperature field is particularly important.
Computer modeling and simulation technology can not only obtain the temperature changes in the ablation zone, but also help doctors to make a reasonable surgical plan before surgery. Zhang [4] established the temperature field models based on the Pennes and Hyperbolic heat transfer equations and discussed the effects of different heat transfer models on the temperature distributions. Nie [5] studied the temperature field of RFA in the spinal tumor conformal therapy. Jin et al. [6] found a RFA model of thyroid treatment based on MRI. Huang et al. [7] attempted to adjust the electrode needle voltage based on the proportional integral feedback control method, to build temperature-controlled RFA simulation model, but did not consider the parameter changes over time. In these literatures, however, the electro-thermal parameters were regarded as fixed values, and the specificity of the biological tissue parameters was not considered. Actually, these parameters will vary with temperature. Rossmann and Haemmerich [8] presented a review of temperature dependence of tissue electro-thermal properties. The study also suggested that all these parameters were necessary for determining absorption of electromagnetic energy (e.g. microwaves or RF current), and estimating the temperature profile. Trujillo and Berjano [9] also investigated the mathematical functions used to describe temperature dependence of electrical and thermal conductivities in RFA models. Thus, using a single reference value in RFA simulation model is not sufficiently accurate. In order to improve the accuracy of the simulation model of liver tumor temperature field, it is necessary to research the influence of changes of tissue parameters on the lesion volume in RFA.
In view of the above, a sensitivity analysis method is required to precisely characterize how the temperature dependence of parameters will influence the simulation temperature and lesion volume. Sensitivity analysis is a quantitative description of the input parameters on the output response [10], and its quantitative index is the variance contribution rate of error. The greater the SS% is, the higher the effects of the parameter on the model responses will be [11]. The purpose of this paper is to achieve the sensitivity of each characteristic parameter by sensitivity analysis method. In the feedback application, the less sensitive parameters can be removed, so as to effectively reduce the complexity of the model, simplify data operations, and further improve the accuracy of the simulation model. This paper mainly analyzes the characteristic parameters of the temperature distribution model by variance analysis method, and obtains the sensitivity parameters which have significant influences on the ablation temperature and lesion volume.
In clinical RF thermal ablation surgery, the temperature-controlled ablation instrument is commonly used. The working principle is to adjust the output power in real time according to the set parameters (center temperature, temperature-raising rate, etc.), acquire a constant central temperature by power compensation, and finally achieve the required thermal coagulation effects. Temperature-controlled RFA model was established based on the ablation instrument RFA-I (Blade Co., Ltd., Beijing, China) and the liver tissue phantom. The temperature field modeling was performed using COMSOL Multiphysics software version 5.1 (COMSOL Inc., Palo Alto, CA, USA). The sensitivity of each characteristic parameter was studied by factor analysis (DOE) in Minitab (State College, PA, USA).
2.Materials and methods
2.1Geometric modeling of RFA
The RFA simulation model was established and solved by the finite element method (FEM) in Comsol software. FEM was used to solve the diffusion equation numerically [12, 13], which offered advantages in speed and flexibility.
Due to the symmetry of the electrode, only half of the model needed to be considered to save computation time. The symmetrical model was constructed into a semi-cylinder of homogenized tissue (50 mm * 25 mm * 70 mm), and a clinical multi-polar electrode was embedded in the tissue as shown in Fig. 1. The shapes and dimensions of the electrode were designed based on the real multi-polar electrode used in the experiments (Fig. 2). Six sub-electrodes were distributed on one side of the multi-polar electrode uniformly. Therefore, the angle between two sub-electrodes was about 55∘. Each sub-electrode had a different degree of deflection, and the deflections of the adjacent two sub-electrodes were 104∘ and 57∘.
To obtain the solution of FEM model, the boundary conditions and initial conditions of the simulation model were set. For the current field, the boundary of the multi-polar RF electrode was set as the voltage source, and the boundary of the liver tissue was set as the ground; for the bio-heat transfer physical field, the initial temperature of the model was 28∘C, the external thermal boundary condition was 28∘C. At the same time, we assumed that the tissue perfusion rate ω had no effect on the simulated temperature in the liver phantom.
2.2Mathematical modeling of RFA
The construction of RFA mathematical model was based on electromagnetic governing equation and bio-heat transfer equation. The heat source of each point was created by electromagnetic energy, and the temperature distribution was formed by bio-heat transfer mechanism.
In the process of RFA simulation, the current field is regarded as a heat source, and the control equation is Laplace equation:
where J represents the current density (A/m2), σ represents the conductivity (S/m), E represents the electric field strength (V/m), V represents voltage source (V) and can be determined by the following formula:
where R(Ω) is the equivalent resistance of the tissue, P(W) is gained according to the ablator schedule, in this study as a function of time:
The bio-heat transfer equation employed is the classical Pennes bio-heat transfer equation, which is described as follows:
where T and Tb are the temperature of tissue and blood (∘C); t is the time (s); ρ and ρb are the density of tissue and blood, respectively (kg/m3); c and cb are the specific heat capacity of tissue and blood, respectively (J/(Kg⋅∘C)); k is the thermal conductivity of tissue (J/(s⋅m⋅∘C)); ω is the volume blood perfusion rate (Kg/(m
2.3Sensitivity analysis method
Thermal and electric parameters were selected as the influence variables in this study. Response surface methodology (RSM) is a sensitivity analysis technique, which can present the correlations between several input variables and response variables. RSM include 3k full factorial experimental design, Central Composite Design (CCD) and Box-Behnken Design (BBD), where CCD is an experimental design based on a two level full factorial and fractional experimental design. By adding a design point to the 2-level test (equivalent to increase one level), CCD can utilize to investigate the non-linear relationships between the input variables and response variable.
Parameter | Unit | Symbol | -1 | 0 | 1 |
---|---|---|---|---|---|
Resistance, R | Ω | A | 16 | 16.8 | 17.6 |
Specific heat capacity, Cp | (J/(Kg⋅∘C)) | B | 3400 | 3570 | 3740 |
Thermal conductivity, k | (W/(m⋅k)) | C | 0.48 | 0.504 | 0.528 |
Electrical conductivity, Sigma | (S/m) | D | 0.2 | 0.21 | 0.22 |
Density, rho | (kg/m3) | E | 1000 | 1050 | 1100 |
Dielectric constant, Epsilon | 1 | F | 40 | 42 | 44 |
Run order | Code value | |||||
---|---|---|---|---|---|---|
A | B | C | D | E | F | |
1 | -1 | 1 | 1 | -1 | -1 | -1 |
2 | 1 | 1 | -1 | -1 | 1 | 1 |
3 | -1 | -1 | -1 | 1 | 1 | -1 |
4 | -1 | -1 | -1 | 1 | -1 | 1 |
5 | 1 | -1 | 1 | -1 | 1 | 1 |
6 | -1 | -1 | 1 | -1 | -1 | 1 |
7 | 1 | -1 | -1 | 1 | 1 | 1 |
8 | 1 | 1 | 1 | -1 | 1 | -1 |
9 | -1 | 1 | 1 | -1 | 1 | 1 |
10 | -1 | 1 | -1 | -1 | 1 | -1 |
11 | -1 | -1 | 1 | 1 | 1 | 1 |
12 | 1 | 1 | 1 | 1 | 1 | 1 |
13 | -1 | 1 | 1 | 1 | 1 | -1 |
14 | 0 | 0 | 0 | 0 | 0 | 0 |
15 | -1 | 1 | -1 | -1 | -1 | 1 |
16 | 0 | 0 | 0 | 0 | 0 | 0 |
17 | 0 | 0 | 0 | 0 | 0 | 0 |
18 | -1 | 1 | -1 | 1 | 1 | 1 |
19 | 0 | 0 | 0 | 0 | 0 | 0 |
20 | -1 | -1 | 1 | 1 | -1 | -1 |
21 | 1 | -1 | 1 | 1 | -1 | 1 |
22 | -1 | 1 | -1 | 1 | -1 | -1 |
23 | -1 | -1 | -1 | -1 | 1 | 1 |
24 | 1 | -1 | -1 | -1 | -1 | 1 |
25 | 1 | -1 | -1 | -1 | 1 | -1 |
26 | 0 | 0 | 0 | 0 | 0 | 0 |
27 | 0 | 0 | 0 | 0 | 0 | 0 |
28 | 0 | 0 | 0 | 0 | 0 | 0 |
29 | -1 | 1 | 1 | 1 | -1 | 1 |
30 | -1 | -1 | -1 | -1 | -1 | -1 |
31 | 1 | 1 | -1 | 1 | 1 | -1 |
32 | 1 | 1 | 1 | -1 | -1 | 1 |
33 | 1 | 1 | 1 | 1 | -1 | -1 |
34 | 1 | 1 | -1 | -1 | -1 | -1 |
35 | 1 | -1 | 1 | -1 | -1 | -1 |
36 | 1 | 1 | -1 | 1 | -1 | 1 |
37 | 1 | -1 | 1 | 1 | 1 | -1 |
38 | 0 | 0 | 0 | 0 | 0 | 0 |
39 | -1 | -1 | 1 | -1 | 1 | -1 |
40 | 1 | -1 | -1 | 1 | -1 | -1 |
41 | 0 | 0 | 0 | 0 | 0 | 1 |
42 | 0 | 0 | 0 | 0 | 0 | 0 |
43 | 1 | 0 | 0 | 0 | 0 | 0 |
44 | -1 | 0 | 0 | 0 | 0 | 0 |
45 | 0 | 0 | 0 | -1 | 0 | 0 |
46 | 0 | 1 | 0 | 0 | 0 | 0 |
47 | 0 | 0 | 0 | 0 | 0 | 0 |
48 | 0 | -1 | 0 | 0 | 0 | 0 |
49 | 0 | 0 | 1 | 0 | 0 | 0 |
50 | 0 | 0 | 0 | 0 | 1 | 0 |
51 | 0 | 0 | -1 | 0 | 0 | 0 |
52 | 0 | 0 | 0 | 0 | 0 | -1 |
53 | 0 | 0 | 0 | 0 | -1 | 0 |
54 | 0 | 0 | 0 | 1 | 0 | 0 |
Taking into account the non-linearity of the parametric effects during RFA, in the present study, CCD experimental design was used to analyze the influences of thermal and electric parameters on the ablation volume. Six parameters were selected at three levels to analyze parameters sensitivity and the statistical model table (with 54 runs) was listed in Table 2. Parameters and their levels were summarized in Table 1, where (-1), (0) and (+1) represented low, middle and high levels respectively. The ranges of these parameters were collected from references [8, 14]. As shown in Table 2, statistical analysis of the 54 runs was performed according to the defined conditions by Minitab software.
3.Results
Some representative points including the far field point (0 mm, 15 mm, -19 mm) and the near field point (0 mm, -5 mm, -19 mm) as well as the 54∘C isothermal surface volume of the ablation zone were analyzed, so as to study the effects of the characteristic parameters on the temperature distributions in RFA model. 54∘C isothermal surface volume here represents the voulme enclosed by the 54∘C isotherm around the electrode and is used interchangeably with the term “ablation volume” or “lesion volume”. Figure 3 showed the positions of representative points in simulation model.
According to the 54 runs in Table 2, the temperature field simulation was implemented in Comsol software, and the simulation time was 360 s. Figures 4 and 5 showed the main effects of the parameters R, Cp, Sigma, rho and Epsilon on the far field temperature and near field temperature, respectively at t= 50 s and t= 360 s. As could be seen in Figs 4 and 5, the influences of Cp and rho on the far field temperature were getting smaller and smaller over time, and these parameters were negatively correlated with the simulation temperature. The effects of R and Sigma on the far field temperature were always obvious, so R and Sigma were significant influencing factors. With regard to Epsilon, there was no effect on ablation temperature. Besides, the SS% trendlines of the parameter k in the near field and far field were different.
Source | DOF | Adj_SS | Adj_MS | F-value | P-value | SS% |
t= 50 s | ||||||
Model | 27 | 0.138602 | 0.005133 | 119.47 | 0.0 | 99.23% |
Linear | 6 | 0.136185 | 0.022698 | 528.25 | 0.0 | 97.50% |
A | 1 | 0.032656 | 0.032656 | 760.01 | 0.0 | 23.38% |
B | 1 | 0.033395 | 0.033395 | 777.22 | 0.0 | 23.91% |
C | 1 | 0.001724 | 0.001724 | 40.12 | 0.0 | 1.23% |
D | 1 | 0.034916 | 0.034916 | 812.62 | 0.0 | 25.00% |
E | 1 | 0.033395 | 0.033395 | 777.22 | 0.0 | 23.91% |
F | 1 | 0.000099 | 0.000099 | 2.29 | 0.142 | 0.07% |
Square | 6 | 0.00026 | 0.000043 | 1.01 | 0.443 | 0.19% |
Interaction | 15 | 0.001595 | 0.000106 | 2.47 | 0.022 | 1.14% |
Error | 25 | 0.001074 | 0.000043 | – | – | – |
Lack-of-fit | 16 | 0.001074 | 0.000067 | – | – | – |
Pure error | 9 | 0.0 | 0.0 | – | – | – |
Total | 52 | 0.139676 | – | – | – | – |
t= 360 s | ||||||
Model | 27 | 11.1469 | 0.41285 | 1420.17 | 0.0 | 99.94% |
Linear | 6 | 11.0018 | 1.83364 | 6307.62 | 0.0 | 98.63% |
A | 1 | 4.728 | 4.72796 | 16263.93 | 0.0 | 42.39% |
B | 1 | 0.539 | 0.53898 | 1854.05 | 0.0 | 4.83% |
C | 1 | 0.1752 | 0.17517 | 602.58 | 0.0 | 1.57% |
D | 1 | 5.0201 | 5.02006 | 17268.73 | 0.0 | 45.01% |
E | 1 | 0.539 | 0.53898 | 1854.05 | 0.0 | 4.83% |
F | 1 | 0.0007 | 0.00069 | 2.38 | 0.136 | 0.01% |
Square | 6 | 0.0 | 0.0 | 0.01 | 1 | 0.0% |
Interaction | 15 | 0.0272 | 0.00181 | 6.23 | 0 | 0.24% |
Error | 25 | 0.0073 | 0.00029 | – | – | – |
Lack-of-fit | 16 | 0.0073 | 0.00045 | – | – | – |
Pure error | 9 | 0.0 | 0.0 | – | – | – |
Total | 52 | 11.1541 | – | – | – | – |
In order to acquire the quantitative analysis of the influences of the parameters on the ablation temperature, the variance analysis was carried out. As shown in Table 3, the F-value is the variance of the mean of the test data. The larger the value is, the more significant the parameter is. The P-value is used to determine the result of the hypothesis test [10]. If the P-value < 0.05, the parameter is significant. DOF represents the degree of freedom; Adj_MS represents the mean of the variance of the parameters; SS% represents the square sum of the variance deviations from the mean, and its value quantitatively reflects the significance of the parameters on the ablation temperature, as shown in Eq. (7):
where i can take A, B, C, D, E, and F; 𝐴𝑑𝑗_SSi represents the variance of the ith parameter.
Source | DOF | Adj_SS | Adj_MS | F-value | P-value | SS% |
---|---|---|---|---|---|---|
t= 50 s | ||||||
Model | 27 | 2.67669 | 0.099137 | 118.29 | 0.0 | 99.22% |
Linear | 6 | 2.62753 | 0.437921 | 522.51 | 0.0 | 97.40% |
A | 1 | 0.7665 | 0.766497 | 914.56 | 0.0 | 28.41% |
B | 1 | 0.51714 | 0.517135 | 617.03 | 0.0 | 19.17% |
C | 1 | 0.0107 | 0.010699 | 12.77 | 0.0 | 0.40% |
D | 1 | 0.81414 | 0.814139 | 971.4 | 0.0 | 30.18% |
E | 1 | 0.51714 | 0.517135 | 617.03 | 0.0 | 19.17% |
F | 1 | 0.00192 | 0.001921 | 2.29 | 0.143 | 0.07% |
Square | 6 | 0.00374 | 0.000624 | 0.74 | 0.619 | 0.14% |
Interaction | 15 | 0.03105 | 0.00207 | 2.47 | 0.022 | 1.15% |
Error | 25 | 0.02095 | 0.000838 | – | – | – |
Lack-of-fit | 16 | 0.02095 | 0.00131 | – | – | – |
Pure error | 9 | 0.0 | 0.0 | – | – | – |
Total | 52 | 2.69765 | – | – | – | – |
t= 360 s | ||||||
Model | 27 | 125.645 | 4.6535 | 383.11 | 0.0 | 99.76% |
Linear | 6 | 123.651 | 20.6085 | 1696.65 | 0.0 | 98.18% |
A | 1 | 58.098 | 58.0984 | 4783.12 | 0.0 | 46.13% |
B | 1 | 1.587 | 1.5868 | 130.64 | 0.0 | 1.26% |
C | 1 | 0.842 | 0.8415 | 69.28 | 0.0 | 0.67% |
D | 1 | 61.51 | 61.5095 | 5063.95 | 0.0 | 48.84% |
E | 1 | 1.587 | 1.5868 | 130.64 | 0.0 | 1.26% |
F | 1 | 0.028 | 0.028 | 2.31 | 0.141 | 0.02% |
Square | 6 | 0.006 | 0.001 | 0.08 | 0.997 | 0.0% |
Interaction | 15 | 0.597 | 0.0398 | 3.27 | 0.004 | 0.47% |
Error | 25 | 0.304 | 0.0121 | – | – | – |
Lack-of-fit | 16 | 0.304 | 0.019 | – | – | – |
Pure error | 9 | 0.0 | 0.0 | – | – | – |
Total | 52 | 125.949 | – | – | – | – |
The SS% values of single parameter in the far field and near field were shown in Tables 3 and 4. The Interaction and Square results between mentioned parameters were also listed. The sensitivity of single parameter, interaction and square item at t= 50 s was: Sigma>Cp (rho) >R>k> Interaction > Square >Epsilon and at t= 360 s was: Sigma> R>Cp (rho) >k> Interaction >Epsilon> Square; where Epsilon and its Interaction and Square item had little influence on RFA temperature.
This paper also showed the SS% trendlines to more clearly illustrate the changes of the parameter sensitivity during the RFA period. As could be seen in Figs 6 and 7, in the far field, the SS% of Cp and rho ranged from 33.4% to 4.8%; the SS% of Sigma and R ranged from 14.4% to 45%; in the near field, the SS% of Cp and rho ranged from 23% to 3.2%; the SS% of Sigma and R ranged from 22.8% to 48.8%. Therefore, the ablation temperature was sensitive to Cp and rho at early ablation stage; Sigma and R had important effects on ablation temperature at the middle and later period.
However, the sensitivity analysis merely based on individual points had some particularity. This study further discussed the sensitivity of the parameters on thermal coagulation zone, which was the main output of RFA. 50∘C isothermal surface [9, 15, 16] and 60∘C isothermal surface [17, 18] were usually used as the ablation zone boundary during the RFA. In another word, 50∘C isothermal surface volume or 60∘C isothermal surface volume could be regarded as the thermal damage volume. In order to highlight the effects of the parameters on the lesion volume, the zone with the temperature more than 54∘C was considered as the ablation volume.
Source | DOF | Adj_SS | Adj_MS | F-value | P-value | SS% |
---|---|---|---|---|---|---|
t= 100 s | ||||||
Model | 27 | 398686 | 14766 | 4427.41 | 0.0 | 99.98% |
Linear | 6 | 379695 | 63282 | 18974.31 | 0.0 | 95.22% |
A | 1 | 140829 | 140829 | 42225.47 | 0.0 | 35.32% |
B | 1 | 37234 | 37234 | 11163.93 | 0.0 | 9.34% |
C | 1 | 13748 | 13748 | 4122.26 | 0.0 | 3.45% |
D | 1 | 150648 | 150648 | 45169.57 | 0.0 | 37.78% |
E | 1 | 37234 | 37234 | 11163.93 | 0.0 | 9.34% |
F | 1 | 2 | 2 | 0.7 | 0.409 | 0.00% |
Square | 6 | 2358 | 393 | 117.83 | 0.0 | 0.59% |
Interaction | 15 | 13620 | 908 | 272.24 | 0.0 | 3.42% |
Error | 25 | 83 | 3 | – | – | – |
Lack-of-fit | 16 | 83 | 5 | – | – | – |
Pure error | 9 | 0.0 | 0.0 | – | – | – |
Total | 52 | 398769 | – | – | – | – |
t= 360 s | ||||||
Model | 27 | 3496230 | 129490 | 5531.04 | 0.0 | 99.98% |
Linear | 6 | 3447637 | 574606 | 24543.72 | 0.0 | 98.59% |
A | 1 | 1594073 | 1594073 | 68089.24 | 0.0 | 45.59% |
B | 1 | 57771 | 57771 | 2467.62 | 0.0 | 1.65% |
C | 1 | 40849 | 40849 | 1744.82 | 0.0 | 1.17% |
D | 1 | 1697170 | 1697170 | 72492.9 | 0.0 | 48.53% |
E | 1 | 57771 | 57771 | 2467.62 | 0.0 | 1.65% |
F | 1 | 3 | 3 | 0.14 | 0.71 | 0.00% |
Square | 6 | 385 | 64 | 2.74 | 0.035 | 0.0% |
Interaction | 15 | 2604 | 174 | 7.42 | 0.0 | 0.07% |
Error | 25 | 585 | 23 | – | – | – |
Lack-of-fit | 16 | 585 | 37 | – | – | – |
Pure error | 9 | 0.0 | 0.0 | – | – | – |
Total | 52 | 3496816 | – | – | – | – |
Figure 8 illustrated the main effects of the parameters on 54∘C isothermal surface volume during the ablation period. The results showed that the influences of R and Sigma on the 54∘C isothermal surface volume were significant. Epsilon had little effect on the isothermal surface volume. With the passing of time, the sensitivity of Cp and rho to the isothermal surface volume became inconspicuous. These conclusions were consistent with the results of single-point sensitivity analysis. Table 5 also showed the variance analysis results of 54∘C isothermal surface volume.
As shown in Fig. 9, the SS% trendlines of 54∘C isothermal surface volume were plotted. The SS% trendlines of parameters Cp and rho were consistent. Figure 8 indicated that the lesion volume was significantly sensitive to Sigma and R during RFA period. Their contribution rate was about 35% at early ablation stage, and was above 46% at the middle and later period of RFA.
4.Discussion
Based on the RFA experiment and FEM model of vitro phantom, the parameters that have a significant influence on ablation zone size of the liver tissue were analyzed. There are two primary contributions: (1) this study has established an accurate RFA model based on the actual power profile and equivalent resistance; (2) because the non-linearity variation of the liver tissue parameters during RFA process, six parameters with three levels each were selected to analyze the sensitivity to ablation temperature and lesion volume. Parameter sensitivity analysis results, presented in this study, could provide a priori information for clinical practitioners and simplify the finite element modeling of liver tumor RFA.
Although many methods can carry out parameter sensitivity analysis, the efficiency of different methods varies widely. In order to estimate the effects of tissue parameters to ablation zone, full factorial experimental design can be used to provide the optimum results by testing all the possible combinations of the parameters. Johansson et al. [19] used a 24-factorial design to investigate the influences of tissue parameters on lesion zone of the brain. Nevertheless, the authors did not consider the non-linearity effects of the parameters during RFA. Furthermore, the more experimental combinations will be performed with increasing input parameters. For instance, with regard to 5 parameters each having three levels, 3
At present, some scholars have conducted sensitivity analysis on the parameters during RFA process. However, there are still some differences between the present sensitivity analysis results. Hall et al. [21] indicated that blood perfusion, electrical conductivity and the cell death model were the most important factors for estimating the ablation zone size in liver RFA models. Similarly, our results obviously demonstrated that electric conductivity was one of the important parameters. Monsalvo et al. [22] analyzed how the variations of the different parameters influenced the temperature in liver tissue through using the complex Taylor series expansion (CTSE) finite element method. The authors have found that the thermal and electrical conductivity of the healthy liver tissue have extremely important influence on the final temperature of RFA simulation models. According to what we has mentioned above, it is clear that electric conductivity is one of the most important parameters, thus confirming the results of our study. Moreover, the insensitivity parameters of the lesion size cannot be focused to simplify the RFA model of liver tumor. However, due to RFA model differences, analysis parameter differences and experimental environment differences (e.g. in vitro experiments or in vivo experiments) during sensitivity analysis, there are a few differences between our results and the previous studies.
Although some important discoveries have been revealed in this paper, there are also limitations. Firstly, the proposed RFA model for predicting the temperature did not consider the cooling effects of large blood vessels. Actually, during RFA of liver tumor, the cooling influence of the large blood vessel (i.e., the blood vessel with the diameter larger than 3 mm) has a remarkable impact on lesion volume [23]. This issue will be taken into account in the future work. Secondly, CCD experimental design needed complex experimental combinations in contrast with the Taguchi’s orthogonal array design. Thirdly, 54∘C isotherm contours was used to assess the extent of the ablation zone in the RFA computer models. However, the death of tissue during RFA depends on the local temperature as well as the RF ablation time [24]. Therefore, Arrhenius [25] models ought to be used in further study. Finally, due to the variability between in vitro and vivoexperiment, the subsequent step is to establish the RF finite element simulation model with blood vessels for sensitivity analysis.
5.Conclusions
In this paper, based on the temperature-controlled RFA simulation model, the influences of the electro-thermal parameters on the temperature distributions and ablation volume were analyzed. According to the main effects and the SS% values, the sensitive parameters were obtained. It can be concluded that the size of the RFA zone was mainly dependent on Sigma and R, and these parameters were positively correlated with ablation zone. The sensitivity of Cp and rho on the ablation zone was reduced over time, and these parameters were negatively correlated with the ablation zone; Epsilon had little effect on the ablation zone during the RFA period. There was no good agreement for the parameter k on the temperature changes. In the near field, the heat was absorbed by the heat exchange at the beginning of RFA, so the temperature will increase with the increase of the k. At the later stage of ablation, heat is released by heat exchange, so k was negatively correlated with temperature. In the far field, heat was mainly absorbed by heat conduction and, as a result, was positively correlated with k. The sensitivity of the parameter k on some points and 54∘C isothermal surface volume was slightly different in the early stage of ablation, and the trends of change in the later ablation session were basically the same.
In conclusion, the most important parameters in predicting the RFA volume are electric conductivity and bioelectric impedance. These results can provide a scientific basis for the feedback regulation of the sensitivity parameters. These extended applications would be the topic of our future study.
Acknowledgments
This work was supported by Beijing Natural Science Foundation [Grant No. 7174279]; National Natural Science Foundation of China [Grant No. 71661167001]; China P-ostdoctoral Science Foundation [Grant No. 2017M620566]; and Postdoctoral Research Fund of Chaoyang District, Beijing [172858].
Conflict of interest
None to report.
References
- 1.
ClasenS ,RemppH ,SchmidtD , et al. Multipolar radiofrequency ablation using internally cooled electrodes in ex vivo bovine liver: Correlation between volume of coagulation and amount of applied energy. European Journal of Radiology2012; 81(1): 11113. doi: 10.1016/j.ejrad.2010.10.031[Google Scholar] - 2.
BerjanoEJ ,HorneroF . Thermal-electrical modeling for epicardial atrial radiofrequencyablation. IEEE Transactions on Biomedical Engineering2004; 51(8): 1348-1357. [PubMed][Google Scholar] - 3.
MatschekJ ,HimmelA ,HaeselerFV , et al. Mathematical modelling and sensitivity analysis of multipolar radiofrequency ablation in the spine. IFAC Papers Online2015; 48(20): 243-248. doi: 10.1016/j.mbs.2016.06.008[Google Scholar] - 4.
ZhangM . Study on surgical planning system and key technology of tumor thermal ablation. Beijing University of Technology, 2016. - 5.
NieXH . Study on temperature field of radiofrequency ablation of spinal tumor. Beijing University of Technology, 2016. - 6.
JinC ,HeZ ,LiuJ . MRI-based finite element simulationon radiofrequency ablation of thyroid cancer. Computer Methods & Programs in Biomedicine2014; 113(2): 529-538. [PubMed][Google Scholar] - 7.
HuangW ,LuoHY ,PanJH , et al. Construction model of temperature-controlled radiofrequency ablation based on CT images of liver. Journal of Interventional Imaging and Therapeutics2014; 11(8): 532-536. [Google Scholar] - 8.
RossmannC ,HaemmerichD . Review of temperature dependence of thermal properties, dielectric properties, and perfusion of biological tissues at hyperthermic and ablation temperatures. Critical Reviews™ in Biomedical Engineering2014; 42(6). doi: 10.1615/CritRevBiomedEng.2015012486[Google Scholar] - 9.
TrujilloM ,BerjanoE . Review of the mathematical functions used to model the temperature dependence of electrical and thermal conductivities of biological tissue in radiofrequency ablation. International Journal of Hyperthermia2013; 29(6): 590-597. doi: 10.3109/02656736.2013.807438[PubMed][Google Scholar] - 10.
AkbarzadehM ,RashidiaS ,BovandM ,EllahiR . A sensitivity analysis on thermal and pumping power for the flow of nanofluid inside a wavy channel. Journal of Molecular Liquids2016; 220: 1-13. [Google Scholar] - 11.
CaiY ,XingY ,HuD . Study on sensitivity analysis. Journal of Beijing Normal University (Natural Science Edition)2008; 44(1): 1. [Google Scholar] - 12.
SchweigerM ,ArridgeSR ,HiraokaM , et al. The finite element method for the propagation of light in scattering media: Boundary and source conditions. Medical Physics1995; 22(11): 1779-1792. doi: 10.1118/1.597634[PubMed][Google Scholar] - 13.
LiJX ,ZhangWJ ,ZhangSX . Finite element method and application status. Science and Technology Innovation Guide2012; 31: 120-121. [Google Scholar] - 14.
CavagnaroM ,PintoR ,LoprestoV .Numerical models of microwave thermal ablation procedures Microwave Conference, IEEE2014; 480-483. doi: 10.1109/eumc.2014.6986475 - 15.
HaemmerichD .Mathematical modeling of impedance controlled radiofrequency tumor ablation and ex-vivo validation International Conference of the IEEE Engineering in Medicine and Biology, IEEE2010; 1605-1608. doi: 10.1109/iembs.2010.5626659 - 16.
PengT ,O’NeillD ,PayneS .Mathematical study of the effects of different intrahepatic cooling on thermal ablation zones International Conference of the IEEE Engineering in Medicine and Biology Society, IEEE2011; 6866-6869. doi: 10.1109/iembs.2011.6091693 - 17.
SunDH ,QianF ,SunG , et al. Clinical application of radiofrequency ablation. Jilin Medical Journal2012; 33(13): 2823-2825. [Google Scholar] - 18.
HaaseS ,SüssP ,SchwientekJ , et al. Radiofrequency ablation planning: An application of semi-infinite modelling techniques. European Journal of Operational Research2012; 218(3): 856-864. doi: 10.1016/j.ejor.2011.12.014[Google Scholar] - 19.
JohanssonJD ,ErikssonO ,WrenJ , et al. Radio-frequency lesioning in brain tissue with coagulation-dependent thermal conductivity: Modelling, simulation and analysis of parameter influence and interaction. Medical and Biological Engineering and Computing2006; 44(9): 757-766. doi: 10.1007/s11517-006-0098-1[PubMed][Google Scholar] - 20.
JamilM ,NgEYK . Quantification of the effect of electrical and thermal parameters on radiofrequency ablation for concentric tumour model of different sizes. Journal of Thermal Biology2015; 51: 23-32. doi: 10.1016/j.jtherbio.2015.03.002[PubMed][Google Scholar] - 21.
HallSK ,OoiEH ,PayneSJ . Cell death, perfusion and electrical parameters are critical in models of hepatic radiofrequency ablation. International Journal of Hyperthermia2015; 31(5): 538-550. doi: 10.3109/02656736.2015.1032370[PubMed][Google Scholar] - 22.
MonsalvoJF ,GarcíaMJ ,MillwaterH , et al. Sensitivity analysis for radiofrequency induced thermal therapies using the complex finite element method. Finite Elements in Analysis and Design2017; 135: 11-21. doi: 10.1016/j.finel.2017.07.001[Google Scholar] - 23.
PochFGM ,RiederC ,BallhausenH , et al. The vascular cooling effect in hepatic multipolar radiofrequency ablation leads to incomplete ablation ex vivo. International Journal of Hyperthermia2016; 32(7): 749-756. [PubMed][Google Scholar] - 24.
ZhangB ,MoserMAJ ,ZhangEM , et al. A review of radiofrequency ablation: Large target tissue necrosis and mathematical modelling. Physica Medica2016; 32(8): 961-971. doi: 10.1016/j.ejmp.2016.07.092[PubMed][Google Scholar] - 25.
PearceJA . Comparative analysis of mathematical models of cell death and thermal damage processes. International Journal of Hyperthermia2013; 29(4): 262-280. doi: 10.3109/02656736.2013.786140[PubMed][Google Scholar]