1D Thermoembolization Model Using CT Imaging Data for Porcine Liver

Rohan Amare,Danielle Stolley,Steve Parrish,Megan Jacobsen,Rick Layman,Chimamanda Santos,Beatrice Riviere,Natalie Fowlkes,David Fuentes,Erik Cressman

2024-09-11

Abstract:Objective: Innovative therapies such as thermoembolization are expected to play an important role in improvising care for patients with diseases such as hepatocellular carcinoma. Thermoembolization is a minimally invasive strategy that combines thermal ablation and embolization in a single procedure. This approach exploits an exothermic chemical reaction that occurs when an acid chloride is delivered via an endovascular route. However, comprehension of the complexities of the biophysics of thermoembolization is challenging. Mathematical models can aid in understanding such complex processes and assisting clinicians in making informed decisions. In this study, we used a Hagen Poiseuille 1D blood flow model to predict the mass transport and possible embolization locations in a porcine hepatic artery. Method: The 1D flow model was used on in vivo embolization imaging data of three pigs. The hydrolysis time constant of acid chloride chemical reaction was optimized for each pig, and LOOCV method was used to test the model's predictive ability. Conclusion: This basic model provided a balanced accuracy rate of 66.8% for identifying the possible locations of embolization in the hepatic artery. Use of the model provides an initial understanding of the vascular transport phenomena that are predicted to occur as a result of thermoembolization.

Medical Physics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

The paper aims to predict possible embolization locations in porcine livers by constructing a 1D thermo - embolization model based on CT image data. Specifically, the researchers used the Hagen - Poiseuille 1D blood flow model to simulate and predict the mass transfer of acid chloride in the porcine hepatic artery and possible embolization locations during the thermo - embolization process.

### Background and motivation

1. **Hepatocellular carcinoma (HCC)**: Hepatocellular carcinoma is a global health problem with an annual incidence of more than 850,000 cases. Due to the increase in risk factors such as non - alcoholic fatty liver disease, its incidence is still rising.
2. **Treatment challenges**: Treatment methods for HCC vary depending on the tumor stage and the extent of underlying liver disease. Surgery (partial hepatectomy) may be curative under ideal conditions, but it is only applicable to about 5% of cases. Most HCC patients are diagnosed at an advanced stage because early symptoms are not obvious.
3. **Minimally invasive treatment methods**: For HCC patients who are not suitable for surgical resection, ablation and embolization are the two most commonly used minimally invasive treatment methods. These methods have been proven to have a survival advantage, but incomplete ablation is relatively common.
4. **Thermo - embolization**: Thermo - embolization is a new trans - arterial cancer treatment method. By delivering acid chloride dissolved in an inert oil solvent to the target tumor, an exothermic chemical reaction is induced. This method combines the advantages of embolization, thermal ablation and chemical ablation, and can simultaneously cause high temperature, ischemia and chemical denaturation of the target tissue.

### Research objectives

1. **Understand the complex biophysics of thermo - embolization**: Understand and predict vascular transport phenomena during the thermo - embolization process through a mathematical model.
2. **Optimize treatment parameters**: Optimize the infusion rate, dose and other treatment parameters of acid chloride through the model to improve the treatment effect.
3. **Predict embolization locations**: Use CT image data and 1D blood flow model to predict possible embolization locations and help clinicians make more informed decisions.

### Methods

1. **Animal experiments**: Conduct thermo - embolization experiments on three pigs and obtain CT image data.
2. **Image processing**: Use the Hessian - based vesselness filter to segment blood vessels and extract the blood vessel centerline.
3. **1D blood flow model**: Use the Hagen - Poiseuille equation to model 1D blood flow and simulate the transport of acid chloride in blood vessels.
4. **Uncertainty analysis**: Conduct uncertainty analysis by changing the hepatic artery blood flow and the minimum threshold to evaluate the predictive ability of the model.

### Results

1. **Model predictive ability**: The basic model achieved a balanced accuracy rate of 66.8% in identifying possible damage locations in the hepatic artery.
2. **Uncertainty analysis**: Through uncertainty analysis, it was found that there is a strong positive correlation between the hepatic artery blood flow (\( q_{\text{ref}} \)) and the optimal hydrolysis time constant (\( \gamma_t \)), while the correlation between the minimum threshold (\( \delta \)) and \( \gamma_t \) is weak and unstable.

### Discussion

1. **Model assumptions**: Although the model has been simplified, the overall balanced accuracy rate reaches 66.8%, showing good predictive ability.
2. **Future work**: Future research can further optimize the model and combine 3D tissue models for more complex coupled simulations to improve the prediction accuracy.

### Formulas

1. **Hagen - Poiseuille equation**:
\[
q_{ij} = k_{ij}(P_i - P_j)
\]
where,
\[
k_{ij} = \frac{\pi R_{ij}^4}{8 \mu L_{ij}}
\]

2. **Mass conservation equation**:
\[
\sum_{i \in N_j} k_{ij}(P_i - P_j) = 0
\]

3. **Inlet boundary condition**:
\[
P_{\text{root}} = P_{\text{MAP}}
\]

4. **Outlet boundary condition**:
\[
\sum_{j }