Search NASA⌕ Search

Engineering topics

Zhou, J

Publications and source records attributed to Zhou, J.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hydraulic fracturing experiments at 1500 m depth in a deep mine: Highlights from the kISMET project

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the Sanford Underground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the Precambrian Poorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the two outermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert (floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20 ft), depending on the unknown amount of leak-off.

Oldenburg, C↗
DOE OSTI2016

SU-C-202-03: A Tool for Automatic Calculation of Delivered Dose Variation for Off-Line Adaptive Therapy Using Cone Beam CT

Purpose: Monitoring the delivered dose is an important task for the adaptive radiotherapy (ART) and for determining time to re-plan. A software tool which enables automatic delivered dose calculation using cone-beam CT (CBCT) has been developed and tested. Methods: The tool consists of four components: a CBCT Colleting Module (CCM), a Plan Registration Moduel (PRM), a Dose Calculation Module (DCM), and an Evaluation and Action Module (EAM). The CCM is triggered periodically (e.g. every 1:00 AM) to search for newly acquired CBCTs of patients of interest and then export the DICOM files of the images and related registrations defined in ARIA followed by triggering the PRM. The PRM imports the DICOM images and registrations, links the CBCTs to the related treatment plan of the patient in the planning system (RayStation V4.5, RaySearch, Stockholm, Sweden). A pre-determined CT-to-density table is automatically generated for dose calculation. Current version of the DCM uses a rigid registration which regards the treatment isocenter of the CBCT to be the isocenter of the treatment plan. Then it starts the dose calculation automatically. The AEM evaluates the plan using pre-determined plan evaluation parameters: PTV dose-volume metrics and critical organ doses. The tool has been tested for 10 patients. Results: Automatic plans are generated and saved in the order of the treatment dates of the Adaptive Planning module of the RayStation planning system, without any manual intervention. Once the CTV dose deviates more than 3%, both email and page alerts are sent to the physician and the physicist of the patient so that one can look the case closely. Conclusion: The tool is capable to perform automatic dose tracking and to alert clinicians when an action is needed. It is clinically useful for off-line adaptive therapy to catch any gross error. Practical way of determining alarming level for OAR is under development.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-32: Do We Need KV Imaging During CBCT Based Patient Set-Up for Lung Radiation Therapy?

Purpose: To evaluate the role of 2D kilovoltage (kV) imaging to complement cone beam CT (CBCT) imaging in a shift threshold based image guided radiation therapy (IGRT) strategy for conventional lung radiotherapy. Methods: A retrospective study was conducted by analyzing IGRT couch shift trends for 15 patients that received lung radiation therapy to evaluate the benefit of performing orthogonal kV imaging prior to CBCT imaging. Herein, a shift threshold based IGRT protocol was applied, which would mandate additional CBCT verification if the applied patient shifts exceeded 3 mm to avoid intraobserver variability in CBCT registration and to confirm table shifts. For each patient, two IGRT strategies: kV + CBCT and CBCT alone, were compared and the recorded patient shifts were categorized into whether additional CBCT acquisition would have been mandated or not. The effectiveness of either strategy was gauged by the likelihood of needing additional CBCT imaging for accurate patient set-up. Results: The use of CBCT alone was 6 times more likely to require an additional CBCT than KV+CBCT, for a 3 mm shift threshold (88% vs 14%). The likelihood of additional CBCT verification generally increased with lower shift thresholds, and was significantly lower when kV+CBCT was used (7% with 5 mm shift threshold, 36% with 2 mm threshold), than with CBCT alone (61% with 5 mm shift threshold, 97% with 2 mm threshold). With CBCT alone, treatment time increased by 2.2 min and dose increased by 1.9 cGy per fraction on average due to additional CBCT with a 3mm shift threshold. Conclusion: The benefit of kV imaging to screen for gross misalignments led to more accurate CBCT based patient localization compared with using CBCT alone. The subsequently reduced need for additional CBCT verification will minimize treatment time and result in less overall patient imaging dose.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-192: A Quick and Effective Method to Validate Patient’s Daily Setup and Geometry Changes Prior to Proton Treatment Delivery Based On Water Equivalent Thickness Projection Imaging (WETPI) for Head Neck Cancer (HNC) Patient

Purpose: With the implementation of Cone-beam Computed-Tomography (CBCT) in proton treatment, we introduces a quick and effective tool to verify the patient’s daily setup and geometry changes based on the Water-Equivalent-Thickness Projection-Image(WETPI) from individual beam angle. Methods: A bilateral head neck cancer(HNC) patient previously treated via VMAT was used in this study. The patient received 35 daily CBCT during the whole treatment and there is no significant weight change. The CT numbers of daily CBCTs were corrected by mapping the CT numbers from simulation CT via Deformable Image Registration(DIR). IMPT plan was generated using 4-field IMPT robust optimization (3.5% range and 3mm setup uncertainties) with beam angle 60, 135, 300, 225 degree. WETPI within CTV through all beam directions were calculated. 3%/3mm gamma index(GI) were used to provide a quantitative comparison between initial sim-CT and mapped daily CBCT. To simulate an extreme case where human error is involved, a couch bar was manually inserted in front of beam angle 225 degree of one CBCT. WETPI was compared in this scenario. Results: The average of GI passing rate of this patient from different beam angles throughout the treatment course is 91.5 ± 8.6. In the cases with low passing rate, it was found that the difference between shoulder and neck angle as well as the head rest often causes major deviation. This indicates that the most challenge in treating HNC is the setup around neck area. In the extreme case where a couch bar is accidently inserted in the beam line, GI passing rate drops to 52 from 95. Conclusion: WETPI and quantitative gamma analysis give clinicians, therapists and physicists a quick feedback of the patient’s setup accuracy or geometry changes. The tool could effectively avoid some human errors. Furthermore, this tool could be used potentially as an initial signal to trigger plan adaptation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-191: 4D Dose Reconstruction of Intensity Modulated Proton Therapy (IMPT) Based On Breathing Probability Density Function (PDF) From 4D Cone Beam Projection Images: A Study for Lung Treatment

Purpose: With energy repainting in lung IMPT, the dose delivered is approximate to the convolution of dose in each phase with corresponding breathing PDF. This study is to compute breathing PDF weighted 4D dose in lung IMPT treatment and compare to its initial robust plan. Methods: Six lung patients were evaluated in this study. Amsterdam shroud image were generated from pre-treatment 4D cone-beam projections. Diaphragm motion curve was extract from the shroud image and the breathing PDF was generated. Each patient was planned to 60 Gy (12GyX5). In initial plans, ITV density on average CT was overridden with its maximum value for planning, using two IMPT beams with robust optimization (5mm uncertainty in patient position and 3.5% range uncertainty). The plan was applied to all 4D CT phases. The dose in each phase was deformed to a reference phase. 4D dose is reconstructed by summing all these doses based on corresponding weighting from the PDF. Plan parameters, including maximum dose (Dmax), ITV V100, homogeneity index (HI=D2/D98), R50 (50%IDL/ITV), and the lung-GTV’s V12.5 and V5 were compared between the reconstructed 4D dose to initial plans. Results: The Dmax is significantly less dose in the reconstructed 4D dose, 68.12±3.5Gy, vs. 70.1±4.3Gy in the initial plans (p=0.015). No significant difference is found for the ITV V100, HI, and R50, 92.2%±15.4% vs. 96.3%±2.5% (p=0.565), 1.033±0.016 vs. 1.038±0.017 (p=0.548), 19.2±12.1 vs. 18.1±11.6 (p=0.265), for the 4D dose and initial plans, respectively. The lung-GTV V12.5 and V5 are significantly high in the 4D dose, 13.9%±4.8% vs. 13.0%±4.6% (p=0.021) and 17.6%±5.4% vs. 16.9%±5.2% (p=0.011), respectively. Conclusion: 4D dose reconstruction based on phase PDF can be used to evaluate the dose received by the patient. A robust optimization based on the phase PDF may even further improve patient care.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-205: Effectiveness of Robust Treatment Planning to Account for Inter- Fractional Variation in Intensity Modulated Proton Therapy for Head Neck Cancer

Purpose: To evaluate the potential benefits of robust optimization in intensity modulated proton therapy(IMPT) treatment planning to account for inter-fractional variation for Head Neck Cancer(HNC). Methods: One patient with bilateral HNC previous treated at our institution was used in this study. Ten daily CBCTs were selected. The CT numbers of the CBCTs were corrected by mapping the CT numbers from simulation CT via Deformable Image Registration. The planning target volumes(PTVs) were defined by a 3mm expansion from clinical target volumes(CTVs). The prescription was 70Gy, 54Gy to CTV1, CTV2, and PTV1, PTV2 for robust optimized(RO) and conventionally optimized(CO) plans respectively. Both techniques were generated by RayStation with the same beam angles: two anterior oblique and two posterior oblique angles. The similar dose constraints were used to achieve 99% of CTV1 received 100% prescription dose while kept the hotspots less than 110% of the prescription. In order to evaluate the dosimetric result through the course of treatment, the contours were deformed from simulation CT to daily CBCTs, modified, and approved by a radiation oncologist. The initial plan on the simulation CT was re-replayed on the daily CBCTs followed the bony alignment. The target coverage was evaluated using the daily doses and the cumulative dose. Results: Eight of 10 daily deliveries with using RO plan achieved at least 95% prescription dose to CTV1 and CTV2, while still kept maximum hotspot less than 112% of prescription compared with only one of 10 for the CO plan to achieve the same standards. For the cumulative doses, the target coverage for both RO and CO plans was quite similar, which was due to the compensation of cold and hot spots. Conclusion: Robust optimization can be effectively applied to compensate for target dose deficit caused by inter-fractional target geometric variation in IMPT treatment planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-382: Volumetric Modulated Arc Therapy (VMAT) Beam Angle Optimization in Pulsed Partial Brain Irradiation (PPBI) for Newly Diagnosed Glioblastoma

Purpose: To optimize VMAT beam parameters in PPBI to minimize treatment time. We investigate the coverage and organs at risk (OR) avoidance capability of shorter arcs with shorter treatment times. Methods: We evaluated the treatment plans for eleven previously treated PPBI patients. Each patient received 46Gy (2Gy×23) to the initial target and an additional 14Gy (2Gy×7) as a sequential boost. Each daily 2-Gy fraction was delivered as ten 0.2-Gy pulses separated by 3-minute intervals using VMAT. Each pulse was delivered using the same arc and covered at least 95% of the PTV with at least 95% of the prescription dose. To optimize the VMAT beam angle, an initial 360° full-arc VMAT plan was implemented. Beam control points and their corresponding dose rates were exported. A curve of the product of control point and dose rate was plotted against treatment beam angle. The optimum angle range was determined from this relationship. We chose the minimum continuous angle range that covered 85% of the area under the curve. Planning parameters, including treatment time for each pulse (T-pulse), PTV coverage, maximum dose (Dmax), homogeneity index (HI=D5/D95), R50 (50%IDL/PTV), and Dmax to ORs, were compared. Results: Mean PTV volume was 364.1±181.5cc. Mean T-pulse of partial-arc beams was 34.3±10.6s, vs. 63.0±1.7s (p<0.001) for that of full-arc beams. No significant differences were found for PTV V95, Dmax and R50, 99.4%±1.2% vs. 99.7%±0.5% (p=0.066), 108.0%±1.2% vs. 107.5%±1.1% (p=0.107), 2.95±0.38 vs. 2.87±0.35 (p=0.165), for the plans with partial-arc and full-arc beams, respectively. However, plans using full-arc do provide better PTV V100 and HI, 96.0%±3.0% vs. 97.2%±2.0% (p=0.025) and 1.06±0.03 vs. 1.04±0.01 (p=0.009). No significant difference was found on Dmax to ORs. Conclusion: PPBI with optimized partial-arc plans are clinically comparable to full-arc plans, while treatment time be significantly reduced, average saving of 287s for a 10-pulse treatment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-202-02: Incorporating Regional Ventilation Function in Predicting Radiation Fibrosis After Concurrent Chemoradiotherapy for Lung Cancer

Purpose: To investigate the incorporation of pre-therapy regional ventilation function in predicting radiation fibrosis (RF) in stage III non-small-cell lung cancer (NSCLC) patients treated with concurrent thoracic chemoradiotherapy. Methods: 37 stage III NSCLC patients were retrospectively studied. Patients received one cycle of cisplatin-gemcitabine, followed by two to three cycles of cisplatin-etoposide concurrently with involved-field thoracic radiotherapy between 46 and 66 Gy (2 Gy per fraction). Pre-therapy regional ventilation images of the lung were derived from 4DCT via a density-change-based image registration algorithm with mass correction. RF was evaluated at 6-months post-treatment using radiographic scoring based on airway dilation and volume loss. Three types of ipsilateral lung metrics were studied: (1) conventional dose-volume metrics (V20, V30, V40, and mean-lung-dose (MLD)), (2) dose-function metrics (fV20, fV30, fV40, and functional mean-lung-dose (fMLD) generated by combining regional ventilation and dose), and (3) dose-subvolume metrics (sV20, sV30, sV40, and subvolume mean-lung-dose (sMLD) defined as the dose-volume metrics computed on the sub-volume of the lung with at least 60% of the quantified maximum ventilation status). Receiver operating characteristic (ROC) curve analysis and logistic regression analysis were used to evaluate the predictability of these metrics for RF. Results: In predicting airway dilation, the area under the ROC curve (AUC) values for (V20, MLD), (fV20, fMLD), and (sV20, and sMLD) were (0.76, 0.70), (0.80, 0.74) and (0.82, 0.80), respectively. The logistic regression p-values were (0.09, 0.18), (0.02, 0.05) and (0.004, 0.006), respectively. With regard to volume loss, the corresponding AUC values for these metrics were (0.66, 0.57), (0.67, 0.61) and (0.71, 0.69), and p-values were (0.95, 0.90), (0.43, 0.64) and (0.08, 0.12), respectively. Conclusion: The inclusion of regional ventilation function improved predictability of radiation fibrosis. Dose-subvolume metrics provided a promising method for incorporating functional information into the conventional dose-volume parameters for outcome assessment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-C-207-04: Reconstruction Artifact Reduction in X-Ray Cone Beam CT Using a Treatment Couch Model

Purpose: to mitigate artifacts induced by the presence of the RT treatment couch in on-board CBCT and improve image quality Methods: a model of a Varian IGRT couch is constructed using a CBCT scan of the couch in air. The model is used to generate a set of forward projections (FP) of the treatment couch at specified gantry angles. The model couch forward projections are then used to process CBCT scan projections which contain the couch in addition to the scan object (Catphan phantom), in order to remove the attenuation component of the couch at any given gantry angle. Prior to pre-processing with the model FP, the Catphan projection data is normalized to an air scan with bowtie filter. The filtered Catphan projections are used to reconstruct the CBCT with an in-house FDK algorithm. The artifact reduction in the processed CBCT scan is assessed visually, and the image quality improvement is measured with the CNR over a few selected ROIs of the Catphan modules. Results: Sufficient match between the forward projected data and the x-ray projections is achieved to allow filtering in attenuation space. Visual improvement of the couch induced artifacts is achieved, with a moderate expense of CNR. Conclusion: Couch model-based correction of CBCT projection data has a potential for qualitative improvement of clinical CBCT scans, without requiring position specific correction data. The technique could be used to produce models of other artifact inducing devices, such as immobilization boards, and reduce their impact on patient CBCT images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-122: The CBCT Dose Calculation Using a Patient Specific CBCT Number to Mass Density Conversion Curve Based On a Novel Image Registration and Organ Mapping Method in Head-And-Neck Radiation Therapy

Purpose: To develop a CBCT HU correction method using a patient specific HU to mass density conversion curve based on a novel image registration and organ mapping method for head-and-neck radiation therapy. Methods: There are three steps to generate a patient specific CBCT HU to mass density conversion curve. First, we developed a novel robust image registration method based on sparseness analysis to register the planning CT (PCT) and the CBCT. Second, a novel organ mapping method was developed to transfer the organs at risk (OAR) contours from the PCT to the CBCT and corresponding mean HU values of each OAR were measured in both the PCT and CBCT volumes. Third, a set of PCT and CBCT HU to mass density conversion curves were created based on the mean HU values of OARs and the corresponding mass density of the OAR in the PCT. Then, we compared our proposed conversion curve with the traditional Catphan phantom based CBCT HU to mass density calibration curve. Both curves were input into the treatment planning system (TPS) for dose calculation. Last, the PTV and OAR doses, DVH and dose distributions of CBCT plans are compared to the original treatment plan. Results: One head-and-neck cases which contained a pair of PCT and CBCT was used. The dose differences between the PCT and CBCT plans using the proposed method are −1.33% for the mean PTV, 0.06% for PTV D95%, and −0.56% for the left neck. The dose differences between plans of PCT and CBCT corrected using the CATPhan based method are −4.39% for mean PTV, 4.07% for PTV D95%, and −2.01% for the left neck. Conclusion: The proposed CBCT HU correction method achieves better agreement with the original treatment plan compared to the traditional CATPhan based calibration method.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-225: CEST Imaging in Head and Neck Cancer Patients

Purpose: Chemical Exchange Saturation Transfer (CEST) imaging is an MRI technique enables the detection and imaging of metabolically active compounds in vivo. It has been used to differentiate tumor types and metabolic characteristics. Unlike PET/CT,CEST imaging does not use isotopes so it can be used on patient repeatedly. This study is to report the preliminary results of CEST imaging in Head and Neck cancer (HNC) patients. Methods: A CEST imaging sequence and the post-processing software was developed on a 3T clinical MRI scanner. Ten patients with Human papilloma virus positive oropharyngeal cancer were imaged in their immobilized treatment position. A 5 mm slice CEST image was acquired (128×128, FOV=20∼24cm) to encompass the maximum dimension of tumor. Twenty-nine off-set frequencies (from −7.8ppm to +7.8 ppm) were acquired to obtain the Z-spectrum. Asymmetry analysis was used to extract the CEST contrasts. ROI at the tumor, node and surrounding tissues were measured. Results: CEST images were successfully acquired and Zspectrum asymmetry analysis demonstrated clear CEST contrasts in tumor as well as the surrounding tissues. 3∼5% CEST contrast in the range of 1 to 4 ppm was noted in tumor as well as grossly involved nodes. Injection of glucose produced a marked increase of CEST contrast in tumor region (∼10%). Motion and pulsation artifacts tend to smear the CEST contrast, making the interpretation of the image contrast difficult. Field nonuniformity, pulsation in blood vesicle and susceptibility artifacts caused by air cavities were also problematic for CEST imaging. Conclusion: We have demonstrated successful CEST acquisition and Z-spectrum reconstruction on HNC patients with a clinical scanner. MRI acquisition in immobilized treatment position is critical for image quality as well as the success of CEST image acquisition. CEST images provide novel contrast of metabolites in HNC and present great potential in the pre- and post-treatment assessment of patients undergoing radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-169: Characterization of Pacemaker/ICD Dose in SAVI HDR Brachytherapy

Purpose: It is important to estimate dose to pacemaker (PM)/Implantable Cardioverter Defibrillator (ICD) before undertaking Accelerated Partial Breast Treatment using High Dose Rate (HDR) brachytherapy. Kim et al. have reported HDR PM/ICD dose using a single-source balloon applicator. To the authors knowledge, there have so far not been any published PM/ICD dosimetry literature for the Strut Adjusted Volume Implant (SAVI, Cianna Medical, Aliso Viejo, CA). This study aims to fill this gap by generating a dose look up table (LUT) to predict maximum dose to the PM/ICD in SAVI HDR brachytherapy. Methods: CT scans for 3D dosimetric planning were acquired for four SAVI applicators (6−1-mini, 6−1, 8−1 and 10−1) expanded to their maximum diameter in air. The CT datasets were imported into the Elekta Oncentra TPS for planning and each applicator was digitized in a multiplanar reconstruction window. A dose of 340 cGy was prescribed to the surface of a 1 cm expansion of the SAVI applicator cavity. Cartesian coordinates of the digitized applicator were determined in the treatment leading to the generation of a dose distribution and corresponding distance-dose prediction look up table (LUT) for distances from 2 to 15 cm (6-mini) and 2 to 20 cm (10–1).The deviation between the LUT doses and the dose to the cardiac device in a clinical case was evaluated. Results: Distance-dose look up table were compared to clinical SAVI plan and the discrepancy between the max dose predicted by the LUT and the clinical plan was found to be in the range (−0.44%, 0.74%) of the prescription dose. Conclusion: The distance-dose look up tables for SAVI applicators can be used to estimate the maximum dose to the ICD/PM, with a potential usefulness for quick assessment of dose to the cardiac device prior to applicator placement.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-192: Commissioning of a Commercial 3D Dose Calculation Program

Purpose: To commission a commercial software package (CSP) that is used as secondary dose calculation check. The CSP uses an independent golden data beam model. However, some parameters can be modified to generate a customer specific model. Plan comparisons and point dose measurements were performed to test if and to what extent the beam model needed adjustment to optimize results. Methods: Beam parameter configurations were compared between the CSP and both TPS. Twelve phantom test plans ranging from simple to complex were generated in two treatment planning systems (TPS). Tests included small field, off axis, EDW, IMRT and VMAT plans. For each plan a point dose was measured to establish ground truth. Lastly, patient plans were compared for both TPS systems and the CSP. Results: Beam parameters agreed within 2%. The output factors for small fields were changed for the 15 MV beam by 2 and 1.5 % for the 1 cm and 2 cm field sizes, respectively. For the 6 MV beam output factors were adjusted by 3−0.8% for field sizes ranging from 1 to 5 cm. The MLC dynamic leaf gap was adjusted by 1.5 mm for 18 MV beam. Differences between the CSP and the TPS were noted in the built-up region. These differences affected the gamma pass rate in the surface region, however this effect is reduced with increasing number of beam angles and does not affect point dose calculations at depth. All IMRT and VMAT plans agreed with the CSP using a gamma pass rate of 95% (3%, 3mm). Conclusion: The CSP is used to verify point doses for all 3D plans generated in our clinic for the last 6 months. No point dose mismatches were encountered since the CSP was implemented. Next, the CSP will be adapted for secondary checks of all IMRT plans. KL had a beta tester agreement with Mobius Medical for an in-kind equipment and software loan.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-AB-201-04: Optimizing the Number of Catheter Implants and Their Tracks for Prostate HDR Brachytherapy

Purpose: In prostate HDR brachytherapy, interstitial implants are placed manually on the fly. The aim for this research is to develop a computer algorithm to find optimal and reliable implant trajectories using minimal number of implants. Methods: Our new algorithm mainly uses these key ideas: (1) positive charged static particles are uniformly placed on the surface of prostate and critical structures such as urethra, bladder, and rectum. (2) Positive charged kinetic particles are placed at a cross-section of the prostate with an initial velocity parallel to the principal implant direction. (3) The kinetic particles move through the prostate, interacting with each other, spreading out, while staying away from the prostate surface and critical structures. The initial velocity ensures that the trajectories observe the curvature constraints of typical implant procedures. (4) The finial trajectories of kinetic particles are smoothed using a third-degree polynomial regression, which become the implant trajectories. (5) The dwelling times and final dose distribution are calculated using least-distance programming. Results: (1) We experimented with previously treated cases. Our plan achieves all prescription goals while reducing the number of implants by 41%! Our plan also has less uniform target dose, which implies a higher dose is delivered to the prostate. (2) We expect future implant procedures will be performed under the guidance of such pre-calculated trajectories. To assess the applicability, we randomly perturb the tracks to mimic the manual implant errors. Our studies showed the impact of these perturbations are negligible, which is compensated by the least distance programming. Conclusions: We developed a new inverse planning system for prostate HDR therapy that can find optimal implant trajectories while minimizing the number of implants. For future work, we plan to integrate our new inverse planning system with an existing needle tracking system.

60 APPLIED LIFE SCIENCES↗