Focal beam distortion and treatment planning in abdominal focused ultrasound surgery

Transducer module apodization to reduce bone heating during focused ultrasound uterine fibroid ablation with phased arrays: A numerical study

BACKGROUND

During magnetic resonance-guided focused ultrasound (MRgFUS) surgery for uterine fibroids, ablation of fibrous tissues in proximity to the hips and spine is challenging due to heating within the bone that can cause patients to experience pain and potentially damage nerves. This far-field bone heating limits the volume of fibroid tissue that is treatable via MRgFUS.

PURPOSE

To investigate transducer module apodization for improving the ratio of focal-to-bone heating (Δ T ratio $\Delta T_{\mathrm{ratio}}$ ) when targeting fibroid tissue close to the hips and spine, to enable MRgFUS treatments closer to the bone.

METHODS

Acoustic and thermal simulations were performed using 3D magnetic resonance imaging (MRI)-derived anatomies of ten patients who underwent MRgFUS ablation for uterine fibroids using a low-frequency (0.5 MHz $0.5 \ \text{MHz}$ ) 6144-element flat fully-populated modular phased array system (Arrayus Technologies Inc., Burlington, Canada) at our institution as part of a larger clinical trial (NCT03323905). Transducer modules (64 elements $64 \ \text{elements}$ per module) whose beams intersected with no-pass zones delineated within the field were identified, their output power levels were reduced by varying blocking percentage levels, and the resulting temperature field distributions were evaluated across multiple sonications near the hip and spine bones in each patient. Acoustic and thermal simulations took approximately20 min $20 \ \text{min}$ (7 min $7 \ \text{min}$ ) and1 min $1 \ \text{min}$ (30 s $30 \ \text{s}$ ) to run for a single near-spine (near-hip) target, respectively.

RESULTS

For all simulated sonications, transducer module blocking improvedΔ T ratio $\Delta T_{\mathrm{ratio}}$ compared to the no blocking case. In just over half of sonications, full module blocking maximizedΔ T ratio $\Delta T_{\mathrm{ratio}}$ (increase of 82% ± $\pm$ 38% in 50% of hip targets and 49% ± $\pm$ 30% in 62% of spine targets vs. no blocking; mean ± SD), at the cost of more diffuse focusing (focal heating volumes increased by 13% ± 13% for hip targets and 39% ± 27% for spine targets) and thus requiring elevated total (hip: 6% ± 17%, spine: 37% ± 17%) and peak module-wise (hip: 65% ± 36%, spine: 101% ± 56%) acoustic power levels to achieve equivalent focal heating as the no blocking control case. In the remaining sonications, partial module blocking provided further improvements in bothΔ T ratio $\Delta T_{\mathrm{ratio}}$ (increased by 29% ± 25% in the hip and 15% ± 12% in the spine) and focal heating volume (decrease of 20% ± 10% in the hip and 34% ± 17% in the spine) relative to the full blocking case. The optimal blocking percentage value was dependent on the specific patient geometry and target location of interest. Although not all individual target locations saw the benefit, element-wise phase aberration corrections improved the averageΔ T ratio $\Delta T_{\mathrm{ratio}}$ compared to the no correction case (increase of 52% ± 47% in the hip, 35% ± 24% in the spine) and impacted the optimal blocking percentage value. Transducer module blocking enabled ablative treatments to be carried out closer to both hip and spine without overheating or damaging the bone (no blocking:42 ± 1 mm $42\pm 1 \ \text{mm}$ /17 ± 2 mm $17 \pm 2 \ \text{mm}$ , full blocking:38 ± 1 mm $38\pm 1 \ \text{mm}$ /8 ± 1 mm $8\pm 1 \ \text{mm}$ , optimal partial blocking:36 ± 1 mm $36\pm 1 \ \text{mm}$ /7 ± 1 mm $7\pm 1 \ \text{mm}$ for hip/spine).

CONCLUSION

The proposed transducer apodization scheme shows promise for improving MRgFUS treatments of uterine fibroids, and may ultimately increase the effective treatment envelope of MRgFUS surgery in the body by enabling tissue ablation closer to bony structures.

Treatment Planning and Aberration Correction Algorithm for HIFU Ablation of Renal Tumors

High-intensity focused ultrasound (HIFU) applications for thermal or mechanical ablation of renal tumors often encounter challenges due to significant beam aberration and refraction caused by oblique beam incidence, inhomogeneous tissue layers, and presence of gas and bones within the beam. These losses can be significantly mitigated through sonication geometry planning, patient positioning, and aberration correction using multielement phased arrays. Here, a sonication planning algorithm is introduced, which uses the simulations to select the optimal transducer position and evaluate the effect of aberrations and acoustic field quality at the target region after aberration correction. Optimization of transducer positioning is implemented using a graphical user interface (GUI) to visualize a segmented 3-D computed tomography (CT)-based acoustic model of the body and to select sonication geometry through a combination of manual and automated approaches. An HIFU array (1.5 MHz, 256 elements) and three renal cell carcinoma (RCC) cases with different tumor locations and patient body habitus were considered. After array positioning, the correction of aberrations was performed using a combination of backpropagation from the focus with an ordinary least squares (OLS) optimization of phases at the array elements. The forward propagation was simulated using a combination of the Rayleigh integral and k-space pseudospectral method (k-Wave toolbox). After correction, simulated HIFU fields showed tight focusing and up to threefold higher maximum pressure within the target region. The addition of OLS optimization to the aberration correction method yielded up to 30% higher maximum pressure compared to the conventional backpropagation and up to 250% higher maximum pressure compared to the ray-tracing method, particularly in strongly distorted cases.

Aberration correction in abdominal histotripsy

In transabdominal histotripsy, ultrasound pulses are focused on the body to noninvasively destroy soft tissues via cavitation. However, the ability to focus is limited by phase aberration, or decorrelation of the ultrasound pulses due to spatial variation in the speed of sound throughout heterogeneous tissue. Phase aberration shifts, broadens, and weakens the focus, thereby reducing the safety and efficacy of histotripsy therapy. This paper reviews and discusses aberration effects in histotripsy and in related therapeutic ultrasound techniques (e.g., high intensity focused ultrasound), with an emphasis on aberration by soft tissues. Methods for aberration correction are reviewed and can be classified into two groups: model-based methods, which use segmented images of the tissue as input to an acoustic propagation model to predict and compensate phase differences, and signal-based methods, which use a receive-capable therapy array to detect phase differences by sensing acoustic signals backpropagating from the focus. The relative advantages and disadvantages of both groups of methods are discussed. Importantly, model-based methods can correct focal shift, while signal-based methods can restore substantial focal pressure, suggesting that both methods should be combined in a 2-step approach. Aberration correction will be critical to improving histotripsy treatments and expanding the histotripsy treatment envelope to enable non-invasive, non-thermal histotripsy therapy for more patients.

Management of unresectable and recurrent intra-abdominal desmoid tumors treated with ultrasound-guided high-intensity focused ultrasound: A retrospective single-center study

To assess the efficacy and safety of ultrasound (US)-guided high-intensity focused ultrasound (HIFU) ablation for treatment of unresectable and recurrent intra-abdominal desmoid tumors. From June 2014 to March 2020, 15 patients with consecutive unresectable and recurrent diseases that pathologically proven to be intra-abdominal desmoid tumors had undergone the treatment of US-guided HIFU ablation. All patients underwent contrast-enhanced magnetic resonance imaging before and after HIFU treatment. Nonperfused volume ratio was used to evaluate the effect of HIFU therapy. Intraprocedural and postprocedural adverse effects and complications are recorded to assess the safety of the therapy. Outcome of HIFU ablation has been investigated through serial contrast-enhanced imaging examinations during follow up. Out of 15 patients 14 of them have successfully completed the whole therapy, 1 patient is ineffective and gives up further treatment. The mean nonperfused volume ratio is 71.1% (95% confidence interval, 3% to 88.2%). During a mean follow up of 29 months (range from 8 to 61 months), the mean tumor volume was reduced by 59% (95% confidence interval, +49% to -100%). No tumor spreads along the treated area in all patients except one. Complications have occurred in 5 patients (33.3%), including bowel rupture (1 case), intra-abdominal abscess (1 case), slight injury to the femoral nerve (1 case), and bone injury (2 cases), the bowel rupture patient underwent surgery; the others have been cured during the follow up. US-guided HIFU ablation is an effective treatment modality for patients suffered from unresectable and recurrent intra-abdominal desmoid tumors.

Effects of phase aberration on transabdominal focusing for a large aperture, lowf-number histotripsy transducer

Objective. Soft tissue phase aberration may be particularly severe for histotripsy due to large aperture and lowf-number transducer geometries. This study investigated how phase aberration from human abdominal tissue affects focusing of a large, strongly curved histotripsy transducer.Approach.A computational model (k-Wave) was experimentally validated withex vivoporcine abdominal tissue and used to simulate focusing a histotripsy transducer (radius: 14.2 cm,f-number: 0.62, central frequencyfc: 750 kHz) through the human abdomen. Abdominal computed tomography images from 10 human subjects were segmented to create three-dimensional acoustic property maps. Simulations were performed focusing at 3 target locations in the liver of each subject with ideal phase correction, without phase correction, and after separately matching the sound speed of water and fat to non-fat soft tissue.Main results.Experimental validation in porcine abdominal tissue showed that simulated and measured arrival time differences agreed well (average error, ∼0.10 acoustic cycles atfc). In simulations with human tissue, aberration created arrival time differences of 0.65μs (∼0.5 cycles) at the target and shifted the focus from the target by 6.8 mm (6.4 mm pre-focally along depth direction), on average. Ideal phase correction increased maximum pressure amplitude by 95%, on average. Matching the sound speed of water and fat to non-fat soft tissue decreased the average pre-focal shift by 3.6 and 0.5 mm and increased pressure amplitude by 2% and 69%, respectively.Significance.Soft tissue phase aberration of large aperture, lowf-number histotripsy transducers is substantial despite low therapeutic frequencies. Phase correction could potentially recover substantial pressure amplitude for transabdominal histotripsy. Additionally, different heterogeneity sources distinctly affect focusing quality. The water path strongly affects the focal shift, while irregular tissue boundaries (e.g. fat) dominate pressure loss.

Comparison of computer simulations and clinical treatment results of magnetic resonance-guided focused ultrasound surgery (MRgFUS) of uterine fibroids

PURPOSE

Magnetic resonance-guided focused ultrasound surgery (MRgFUS) can be used to noninvasively treat symptomatic uterine fibroids by heating with focused ultrasound sonications while monitoring the temperature with magnetic resonance (MR) thermometry. While prior studies have compared focused ultrasound simulations to clinical results, studies involving uterine fibroids remain scarce. In our study, we perform such a comparison to assess the suitability of simulations for treatment planning.

METHODS

Sonications (N = 67) were simulated retrospectively using acoustic and thermal models based on the Rayleigh integral and Pennes bioheat equation followed by MR-thermometry simulation in seven patients who underwent MRgFUS treatment for uterine fibroids. The spatial accuracy of simulated focus location was assessed by evaluating displacements of the centers of mass of the thermal dose distributions between simulated and treatment MR thermometry slices. Temperature-time curves and sizes of 240 equivalent minutes at 43°C (240EM43 ) volumes between treatment and simulation were compared.

RESULTS

The simulated focus location showed errors of 2.7 ± 4.1, -0.7 ± 2.0, and 1.3 ± 1.2 mm (mean ± SD) in the anterior-posterior, foot-head, and right-left directions for a fibroid absorption coefficient of 4.9 Np m-1  MHz-1 and perfusion parameter of 1.89 kg m-3 s-1 . Linear regression of 240EM43 volumes of 67 sonications of patient treatments and simulations utilizing these parameters yielded a slope of 1.04 and a correlation coefficient of 0.54. The temperature rise ratio of simulation to treatment near the end of sonication was 0.47 ± 0.22, 1.28 ± 0.60, and 1.49 ± 0.71 for 66 sonications simulated utilizing fibroid absorption coefficient of 1.2, 4.9, and 8.6 Np m-1  MHz-1 , respectively, and the aforementioned perfusion value. The impact of perfusion on peak temperature rise is minimal between 1.89 and 10 kg m-3 s-1 , but became more substantial when utilizing a value of 100 kg m-3 s-1 .

CONCLUSIONS

The results of this study suggest that perfusion, while in some cases having a substantial impact on thermal dose volumes, has less impact than ultrasound absorption for predicting peak temperature elevation at least when using perfusion parameter values up to 10 kg m-3 s-1 for this particular array geometry, frequencies, and tissue target which is good for clinicians to be aware of. The results suggest that simulations show promise in treatment planning, particularly in terms of spatial accuracy. However, in order to use simulations to predict temperature rise due to a sonication, knowledge of the patient-specific tissue parameters, in particular the absorption coefficient is important. Currently, spatially varying patient-specific tissue parameter values are not available during treatment, so simulations can only be used for planning purposes to estimate sonication performance on average.

Receiver array design for sonothrombolysis treatment monitoring in deep vein thrombosis

High intensity focused ultrasound (HIFU) can disintegrate blood clots through the generation and stimulation of bubble clouds within thrombi. This work examined the design of a device to image bubble clouds for monitoring cavitation-based HIFU treatments of deep vein thrombosis (DVT). Acoustic propagation simulations were carried out on multi-layered models of the human thigh using two patient data sets from the Visible Human Project. The design considerations included the number of receivers (32, 64, 128, 256, and 512), their spatial positioning, and the effective angular array aperture (100° and 180° about geometric focus). Imaging array performance was evaluated for source frequencies of 250, 750, and 1500 kHz. Receiver sizes were fixed relative to the wavelength (pistons, diameter  =  λ/2) and noise was added at levels that scaled with receiver area. With a 100° angular aperture the long axis size of the  -3 dB main lobe was ~1.2λ-i.e. on the order of the vessel diameter at 250 kHz (~7 mm). Increasing the array aperture to span 180° about the geometric focus reduced the long axis by a factor of ~2. The smaller main lobe sizes achieved by imaging at higher frequencies came at the cost of increased levels of sensitivity to phase aberrations induced during acoustic propagation through the intervening soft tissue layers. With noise added to receiver signals, images could be reconstructed with peak sidelobe ratios  <  -3 dB using single-cycle integration times for source frequencies of 250 and 750 kHz (NRx  ⩾  128). At 1500 kHz, longer integration times and/or higher element counts were required to achieve similar peak sidelobe ratios. Our results suggest that a modest number of receivers(i.e. NRx  =  128) arranged on a semi-cylindrical shell may be sufficient to enable passive acoustic imaging with single-cycle integration times (i.e. volumetric rates up to 0.75 MHz) for monitoring cavitation-based HIFU treatments of DVT.

Patient Specific Simulation of HIFU Kidney Tumour Ablation

High Intensity Focussed Ultrasound (HIFU) is emerging as a non-invasive treatment for localised renal tumours. However, challenges remain in the delivery of the treatment to tumours at depth, with clinical results showing a variation in the ablation efficacy. One clinical trial conducted at the Churchill hospital, Oxford, to investigate the applicability of HIFU for renal tumour ablation found that in 4/10 patients less than 5% of the tumour volume was ablated successfully. The current study looks at the role tissue geometry has on the resulting focal pressure and focal heating. CT scans from 4 patients within the trial were selected, who experienced 70%, < 5%, < 5% and 95% ablation of the target tumour. The CT scans were segmented into bone, fat, kidney, and generic tissue. Full three-dimensional ultrasound simulations were carried out using k-Wave (an open source Matlab toolbox) and for three patients a tight focus was achieved in the kidney but peak pressures varied by 20%. While in the fourth patient there was significant fragmentation of the -6 dB focal volume due to the intervening ribcage. Thermal simulations were used to compare the temperature rise induced across the different patient models. For the three patients with a tight focus, the predicted 47°C iso-volume of the patient with 70% ablation was 2-3 times larger than the two patients with < 5% ablation. For the patient in which the ribcage resulted in focal fragmentation the thermal simulation predicted just a 1°C temperature rise.

Equivalent-Source Acoustic Holography for Projecting Measured Ultrasound Fields Through Complex Media

Holographic projections of experimental ultrasound measurements generally use the angular spectrum method or Rayleigh integral, where the measured data are imposed as a Dirichlet boundary condition. In contrast, full-wave models, which can account for more complex wave behavior, often use interior mass or velocity sources to introduce acoustic energy into the simulation. Here, a method to generate an equivalent interior source that reproduces the measurement data is proposed based on gradient-based optimization. The equivalent-source can then be used with full-wave models (for example, the open-source k-Wave toolbox) to compute holographic projections through complex media including nonlinearity and heterogeneous material properties. Numerical and experimental results using both time-domain and continuous-wave sources are used to demonstrate the accuracy of the approach.

Estimation of the distribution of low-intensity ultrasound mechanical index as a parameter affecting the proliferation of spermatogonia stem cells in vitro

Considering the use of physical and mechanical stimulation, such as low-intensity ultrasound for proliferation and differentiation of stem cells, it is essential to understand the physical and acoustical mechanisms of acoustic waves in vitro. Mechanical index is used for quantifying acoustic cavitation and the relationship between acoustic pressure and the frequency. In this study, modeling of the mechanical index was applied to provide treatment protocol and to understand the effective physical processes on reproducibility of stem cells. Due to low intensity of ultrasound, Rayleigh integral model has been used for acoustic pressure computation. The acoustic pressure and mechanical index equations are modeled and solved to estimate optimal mechanical index for 28, 40, 150kHz and 1MHz frequencies. This model are solved in different intensities and distances from transducer in cylindrical coordinates. Based on the results of the mechanical index, regions with threshold mechanical index of 0.7 were identified for extracting of radiation arrangement to cell medium. Acoustic pressure distribution along the axial and radial was extracted. In order to validate the results of the modeling, the acoustic pressure in the water and near field depth was measured by a piston hydrophone. Results of modeling and experiments show that the model is consistent well to experimental results with 0.91 and 0.90 correlation of coefficient (p<0.05) for 1MHz and 40kHz. Low-intensity ultrasound with 0.40 mechanical index is more effective on enhancing the proliferation rate of the spermatogonia stem cells during the seven days of culture. In contrast, higher mechanical index has a harmful effect on the spermatogonial stem cells. Thus, considering cavitation threshold of different materials is necessary to find effective mechanical index ranges on proliferation for the used frequencies. This acoustic propagation model and ultrasound mechanical index assessments can be used with acceptable accuracy, for the extraction special arrangement of acoustic exposure used in biological conditions in vitro. This model provides proper treatment planning in vitro and in vivo by estimating the cavitation phenomenon.

Frequency considerations for deep ablation with high-intensity focused ultrasound: A simulation study

PURPOSE

The objective of this study was to explore frequency considerations for large-volume, deep thermal ablations with focused ultrasound. Though focal patterns, focal steering rate, and the size of focal clusters have all been explored in this context, frequency studies have generally explored shallower depths and hyperthermia applications. This study examines both treatment efficiency and near-field heating rate as functions of frequency and depth.

METHODS

Flat, 150 mm transducer arrays were simulated to operate at frequencies of 250, 500, 750, 1000, 1250, and 1500 kHz. Each array had λ2 interelement spacing yielding arrays of 2000-70 000 piston-shaped elements arranged in concentric rings. Depths of 50, 100, and 150 mm were explored, with attenuation (α) values of 2.5-10 (Np/m)/MHz. Ultrasound propagation was simulated with the Rayleigh-Sommerfeld integral over a volume of homogeneous simulated tissue. Absorbed power density was determined from the acoustic pressure which, in turn, was modeled with the Pennes bioheat transfer equation. Using this knowledge of temperature over time, thermal dose function of Sapareto and Dewey was used to model the resulting bioeffect of each simulated sonication. Initially, single foci at each depth, frequency, and α were examined with either fixed peak temperatures or fixed powers. Based on the size of the resulting, single foci lesions, larger compound sonications were designed with foci packed together in multiple layers and rings. For each depth, focal patterns were chosen to produce a similar total ablated volume for each frequency. These compound sonications were performed with a fixed peak temperature at each focus. The resulting energy efficiency (volume ablated per acoustic energy applied), near-field heating rate (temperature increase in the anterior third of the simulation space per unit volume ablated), and near- and far-field margins were assessed.

RESULTS

Lesions of comparable volume were created with different frequencies at different depths. The results reflect the interconnected nature of frequency as it effects focal size (decreasing with frequency), peak pressure (generally increasing with frequency), and attenuation (also increasing with frequency). The ablation efficiency was the highest for α = 5 (Np/m)/MHz at a frequency of 750 kHz at each depth. For α = 10 (Np/m)/MHz, efficiency was the highest at 750 kHz for a depth of 50 mm, and 500 kHz at depths of 100 and 150 mm. At all sonication depths, near-field heating was minimized with lower frequencies of 250 and 500 kHz.

CONCLUSIONS

Large-volume ablations are most efficient at frequencies of 500-750 kHz at depths of 100-150 mm. When one considers that near-field heat accumulation tends to be the rate limiting factor in large-volume ablations like uterine fibroid surgery, the results show that frequencies as low as 500 kHz are favored for their ability to reduce heating in the near-field.

Cavitation-enhanced back projection for acoustic rib detection and attenuation mapping

High-intensity focused ultrasound allows for minimally invasive, highly localized cancer therapies that can complement surgical procedures or chemotherapy. For high-intensity focused ultrasound interventions in the upper abdomen, the thoracic cage obstructs and aberrates the ultrasonic beam, causing undesired heating of healthy tissue. When a phased array therapeutic transducer is used, such complications can be minimized by applying an apodization law based on analysis of beam path obstructions. In this work, a rib detection method based on cavitation-enhanced ultrasonic reflections is introduced and validated on a porcine tissue sample containing ribs. Apodization laws obtained for different transducer positions were approximately 90% similar to those obtained using image analysis. Additionally, the proposed method provides information on attenuation between transducer elements and the focus. This principle was confirmed experimentally on a polymer phantom. The proposed methods could, in principle, be implemented in real time for determination of the optimal shot position in intercostal high-intensity focused ultrasound therapy.

MR-guided high-focused ultrasound for renal sympathetic denervation-a feasibility study in pigs

Background

Renal sympathetic denervation has recently gained clinical relevance for the treatment of therapy-resistant hypertension. Denervation is currently mainly performed using catheter-based transarterial radiofrequency ablation of periarterial sympathetic nerve fibers. Since this approach has numerous limitations, we conducted a study to evaluate the feasibility, safety, and efficacy of magnetic resonance-guided high-focused ultrasound (MRgHiFUS) for renal sympathetic denervation in pigs as an alternative to catheter-based ablation.

Methods

Renal periarterial MRgHiFUS was performed under general anesthesia in ten pigs. Blood pressure measurements and magnetic resonance imaging (MRI) of the kidneys, renal arteries, and surrounding structures were obtained immediately before and after the interventions and after 4 weeks. Histological examinations of periarterial tissues and determination of renal norepinephrine (NE) concentration were performed to assess treatment efficacy.

Results and discussion

In each pig, 9.8 ± 2.6 sonications with a mean energy deposition of 2,670 ± 486 J were performed. The procedure was well tolerated by all pigs. No major complications occurred. MRgHiFUS induced periarterial edema in three pigs, but only one pig showed corresponding histological changes. The NE level of the treated kidney was lower in five pigs (-8% to -38%) compared to the untreated side. Overall, there was no significant difference between the NE values of both kidneys in any of the treated pigs. Postinterventional MRI indicated absorption of ultrasound energy at the transverse process and fascia.

Conclusion

MRgHiFUS had some thermal periarterial effects but failed to induce renal denervation. Insufficient energy deposition is most likely attributable to a small acoustic window with beam path impediment in the porcine model. Since HiFUS treatment in humans is expected to be easier to perform due to better access to renal sympathetic nerves, further studies of this method are desirable to investigate the potential of MRgHiFUS as an alternative for patients not suitable for catheter-based renal sympathicolysis.

Design and experimental evaluation of a 256-channel dual-frequency ultrasound phased-array system for transcranial blood-brain barrier opening and brain drug delivery

Focused ultrasound (FUS) in the presence of microbubbles can bring about transcranial and local opening of the blood-brain barrier (BBB) for potential noninvasive delivery of drugs to the brain. A phased-array ultrasound system is essential for FUS-BBB opening to enable electronic steering and correction of the focal beam which is distorted by cranial bone. Here, we demonstrate our prototype design of a 256-channel ultrasound phased-array system for large-region transcranial BBB opening in the brains of large animals. One of the unique features of this system is the capability of generating concurrent dual-frequency ultrasound signals from the driving system for potential enhancement of BBB opening. A wide range of signal frequencies can be generated (frequency = 0.2-1.2 MHz) with controllable driving burst patterns. Precise output power can be controlled for individual channels via 8-bit duty-cycle control of transistor-transistor logic signals and the 8-bit microcontroller-controlled buck converter power supply output voltage. The prototype system was found to be in compliance with the electromagnetic compatibility standard. Moreover, large animal experiments confirmed the phase switching effectiveness of this system, and induction of either a precise spot or large region of BBB opening through fast focal-beam switching. We also demonstrated the capability of dual-frequency exposure to potentially enhance the BBB-opening effect. This study contributes to the design of ultrasound phased arrays for future clinical applications, and provides a new direction toward optimizing FUS brain drug delivery.

Stochastic ray tracing for simulation of high intensity focal ultrasound therapy

An algorithm is presented for rapid simulation of high-intensity focused ultrasound (HIFU) fields. Essentially, the method combines ray tracing with Monte Carlo integration to evaluate the Rayleigh-Sommerfeld integral. A large number of computational particles, phonons, are distributed among the elements of a phase-array transducer. The phonons are emitted into random directions and are propagated along trajectories computed with the ray tracing method. As the simulation progresses, an improving stochastic estimate of the acoustic field is obtained. The method can adapt to complicated geometries, and it is well suited to parallelization. The method is verified against reference simulations and pressure measurements from an ex vivo porcine thoracic tissue sample. Results are presented for acceleration with graphics processing units (GPUs). The method is expected to serve in applications, where flexibility and rapid computation time are crucial, in particular clinical HIFU treatment planning.

MR-guided focused ultrasound surgery, present and future

MR-guided focused ultrasound surgery (MRgFUS) is a quickly developing technology with potential applications across a spectrum of indications traditionally within the domain of radiation oncology. Especially for applications where focal treatment is the preferred technique (for example, radiosurgery), MRgFUS has the potential to be a disruptive technology that could shift traditional patterns of care. While currently cleared in the United States for the noninvasive treatment of uterine fibroids and bone metastases, a wide range of clinical trials are currently underway, and the number of publications describing advances in MRgFUS is increasing. However, for MRgFUS to make the transition from a research curiosity to a clinical standard of care, a variety of challenges, technical, financial, clinical, and practical, must be overcome. This installment of the Vision 20∕20 series examines the current status of MRgFUS, focusing on the hurdles the technology faces before it can cross over from a research technique to a standard fixture in the clinic. It then reviews current and near-term technical developments which may overcome these hurdles and allow MRgFUS to break through into clinical practice.

Computational study on the propagation of strongly focused nonlinear ultrasound in tissue with rib-like structures

The presence of a rib cage is a significant hindrance to the potential applications of focused ultrasound as a noninvasive extracorporeal surgery modality for various internal organs. Here the influence of ribs on the propagation of strongly focused high-intensity nonlinear ultrasound beam inside the body is studied. Based on the spheroidal beam equation, a three-dimensional numerical algorithm is developed to solve the nonlinear acoustic field generated by a focused ultrasonic transducer with a large aperture angle. Idealized ribs, of rectangular cross sectional, with high absorption and impedance, and various dimensions, are used to simulate human anatomical configurations. The changes in the spatial distribution of acoustic intensity and the reduction of the acoustic pressure amplitude and heat deposition rate due to the presence of "ribs" are investigated. It is somewhat surprising that in some cases, the axial peak positions shift less than 2 mm and more than 80% of the sound energy can propagate through the space of the rib cage in the strongly focused sound field. This study also includes quantitative analyses of the effects of different rib configurations and transducers of various f-numbers. The results can be used as reference information for further study and clinical applications.

Pre-clinical study of in vivo magnetic resonance-guided bubble-enhanced heating in pig liver

Bubble-enhanced heating (BEH) can be exploited to increase heating efficiency in treatment of liver tumors with non-invasive high-intensity focused ultrasound (HIFU). The objectives of this study were: (i) to demonstrate the feasibility of increasing the heating efficiency of sonication exploiting BEH in pig liver in vivo using a clinical platform; (ii) to determine the acoustic threshold for such effects with real-time, motion-compensated magnetic resonance-guided thermometry; and (iii) to compare the heating patterns and thermal lesion characteristics resulting from continuous sonication and sonication including a burst pulse. The threshold acoustic power for generation of BEH in pig liver in vivo was determined using sonication of 0.5-s duration ("burst pulse") under real-time magnetic resonance thermometry. In a second step, experimental sonication composed of a burst pulse followed by continuous sonication (14.5 s) was compared with conventional sonication (15 s) of identical energy (1.8 kJ). Modification of the heating pattern at the targeted region located at a liver depth between 20 and 25 mm required 600-800 acoustic watts. The experimental group exhibited near-spherical heating with 40% mean enhancement of the maximal temperature rise as compared with the conventional sonication group, a mean shift of 7 ± 3.3 mm toward the transducer and reduction of the post-focal temperature increase. Magnetic resonance thermometry can be exploited to control acoustic BEH in vivo in the liver. By use of experimental sonication, more efficient heating can be achieved while protecting tissues located beyond the focal point.

Simulation techniques in hyperthermia treatment planning

Abstract Clinical trials have shown that hyperthermia (HT), i.e. an increase of tissue temperature to 39-44 °C, significantly enhance radiotherapy and chemotherapy effectiveness [1]. Driven by the developments in computational techniques and computing power, personalised hyperthermia treatment planning (HTP) has matured and has become a powerful tool for optimising treatment quality. Electromagnetic, ultrasound, and thermal simulations using realistic clinical set-ups are now being performed to achieve patient-specific treatment optimisation. In addition, extensive studies aimed to properly implement novel HT tools and techniques, and to assess the quality of HT, are becoming more common. In this paper, we review the simulation tools and techniques developed for clinical hyperthermia, and evaluate their current status on the path from 'model' to 'clinic'. In addition, we illustrate the major techniques employed for validation and optimisation. HTP has become an essential tool for improvement, control, and assessment of HT treatment quality. As such, it plays a pivotal role in the quest to establish HT as an efficacious addition to multi-modality treatment of cancer.

Large improvement of the electrical impedance of imaging and high-intensity focused ultrasound (HIFU) phased arrays using multilayer piezoelectric ceramics coupled in lateral mode

With a change in phased-array configuration from one dimension to two, the electrical impedance of the array elements is substantially increased because of their decreased width (w)-to-thickness (t) ratio. The most common way to compensate for this impedance increase is to employ electrical matching circuits at a high cost of fabrication complexity and effort. In this paper, we introduce a multilayer lateral-mode coupling method for phased-array construction. The direct comparison showed that the electrical impedance of a single-layer transducer driven in thickness mode is 1/(n²(1/(w/t))²) times that of an n-layer lateral mode transducer. A large reduction of the electrical impedance showed the impact and benefit of the lateral-mode coupling method. A one-dimensional linear 32-element 770-kHz imaging array and a 42-element 1.45-MHz high-intensity focused ultrasound (HIFU) phased array were fabricated. The averaged electrical impedances of each element were measured to be 58 Ω at the maximum phase angle of -1.2° for the imaging array and 105 Ω at 0° for the HIFU array. The imaging array had a center frequency of 770 kHz with an averaged -6-dB bandwidth of approximately 52%. For the HIFU array, the averaged maximum surface acoustic intensity was measured to be 32.8 W/cm² before failure.

Acoustic accessibility investigation for ultrasound mediated treatment of glycogen storage disease type Ia patients

Glycogen storage disease type Ia (GSDIa) is caused by an inherited defect in the glucose-6-phosphatase gene. The recent advent of targeted ultrasound-mediated delivery (USMD) of plasmid DNA (pDNA) to the liver in conjunction with microbubbles may provide an alternative treatment option. This study focuses on determining the acoustically accessible liver volume in GSDIa patients using transducer models of various geometries with an image-based geometry-driven approach. Results show that transducers with longer focal lengths and smaller apertures (up to an f/number of 2) are able to access larger liver volumes in GSDIa patients while still being capable of delivering the required ultrasound dose in situ (2.5 MPa peak negative pressure at the focus). With sufficiently large acoustic windows and the ability to use glucose to easily assess efficacy, GSD appears to be a good model for testing USMD as proof of principle as a potential therapy for liver applications in general.

Focus shift and phase correction in soft tissues during focused ultrasound surgery

During the treatment of soft tissue tumors using focused ultrasound surgery (FUS), the focus can shift away from the desired point due to tissue inhomogeneity. In this paper, a numerical method to calculate the focus shift in multiple-layered tissues and a faster phase-correction method to restore the focus were developed. Data from the simulations showed that the focus shifted about 2 mm along the transducer axis in multiple-layered soft tissues. After phase correction, the focus was restored at the desired point. The ex vivo experiments were conducted to verify the simulations, and the results agreed well with those of the simulations. An empirical formula was obtained to estimate the focus shift in a two-layered water-tissue model and was verified by numerical calculations. Moreover, the focus shift in multiple-layered tissues can be summed by the shifts in the component of each layer of tissues. The factors affecting the focus shift were studied. The focus shift varied linearly with the tissue acoustic speed and tissue thickness, whereas it slightly changed with transducer F number (radius of curvature/diameter). Overall, the findings of this study can help in the development of a better treatment plan for FUS in soft tissues.

A method for MRI guidance of intercostal high intensity focused ultrasound ablation in the liver

PURPOSE

High intensity focused ultrasound (HIFU) is a promising method for the noninvasive treatment of liver tumors. However, the presence of ribs in the HIFU beam path remains problematic since it may lead to adverse effects (skin burns) by absorption and reflection of the incident beam at or near the bone surface. This article presents a method based on magnetic resonance (MR) imaging for identification of the ribs in the HIFU beam, and for selection of the transducer elements to deactivate.

METHODS

The ribs are visualized on anatomical images acquired prior to heating and manually segmented. The resulting regions of interest surrounding the ribs are projected onto the transducer surface by ray tracing from the focal point. The transducer elements in the "shadow" of the ribs are then deactivated. The method was validated ex vivo and in vivo in pig liver during breathing under multislice real-time MR thermometry, using the proton resonance frequency shift method.

RESULTS

Ex vivo and in vivo temperature data showed that the temperature increase near the ribs was substantial when HIFU sonications were performed with all elements active, whereas the temperature was reduced with deactivation of the transducer elements located in front of the ribs. The temperature at the focal point was similar with and without deactivation of the transducer elements, indicative of no loss of heat efficiency with the proposed technique.

CONCLUSIONS

This method is simple, rapid, and reliable, and enables intercostal HIFU ablation while sparing ribs and their surrounding tissues.

Focal beam distortion and treatment planning for transrib focused ultrasound thermal therapy: a feasibility study using a two-dimensional ultrasound phased array

PURPOSE

The purpose of this study is to numerically investigate the feasibility of employing a spherical-section ultrasound phased array for transrib thermal ablation of liver tumors.

METHODS

Based on CT images, the authors performed a 3D reconstruction of the ribs and the surrounding soft tissues. A 3D pseudospectral time-domain (PSTD) solver was used to assess wave propagation and the distribution of pressure, with the aim of determining the specific absorption rate (SAR) and the resulting thermal doses and dynamics. Phase aberrations caused by the interposed ribs were corrected to assess the efficacy of the device in improving the SAR gain between the ribs and the target positions.

RESULTS

Experimental results supported the usefulness of the PSTD solver for predicting the pressure distribution due to the interfering obstacle. In addition, the method allowed the correction of phase aberrations caused by the ribs, and a significant improvement (176%) in the SAR gain between the ribs and the target points was observed at specific frequencies.

CONCLUSIONS

The method allowed successful tissue targeting without causing overheating of the ribs. One main advantage of this approach is the accurate estimation of phase aberration caused by heterogeneously porous ribs and other interposed tissues. This strategy might prove useful to assess the effectiveness and safety of focused ultrasound thermal ablation prior to transrib treatment.

MRI-guided focused ultrasound treatments

Focused ultrasound (FUS) allows noninvasive focal delivery of energy deep into soft tissues. The focused energy can be used to modify and eliminate tissue for therapeutic purposes while the energy delivery is targeted and monitored using magnetic resonance imaging (MRI). MRI compatible methods to deliver these exposures have undergone rapid development over the past 10 years such that clinical treatments are now routinely performed. This paper will review the current technical and clinical status of MRI-guided focused ultrasound therapy and discuss future research and development opportunities.

Computational aspects in high intensity ultrasonic surgery planning

Therapeutic ultrasound treatment planning is discussed and computational aspects regarding it are reviewed. Nonlinear ultrasound simulations were solved with a combined frequency domain Rayleigh and KZK model. Ultrasonic simulations were combined with thermal simulations and were used to compute heating of muscle tissue in vivo for four different focused ultrasound transducers. The simulations were compared with measurements and good agreement was found for large F-number transducers. However, at F# 1.9 the simulated rate of temperature rise was approximately a factor of 2 higher than the measured ones. The power levels used with the F# 1 transducer were too low to show any nonlinearity. The simulations were used to investigate the importance of nonlinarities generated in the coupling water, and also the importance of including skin in the simulations. Ignoring either of these in the model would lead to larger errors. Most notably, the nonlinearities generated in the water can enhance the focal temperature by more than 100%. The simulations also demonstrated that pulsed high power sonications may provide an opportunity to significantly (up to a factor of 3) reduce the treatment time. In conclusion, nonlinear propagation can play an important role in shaping the energy distribution during a focused ultrasound treatment and it should not be ignored in planning. However, the current simulation methods are accurate only with relatively large F-numbers and better models need to be developed for sharply focused transducers.

Current status and future potential of MRI-guided focused ultrasound surgery

The combination of the imaging abilities of magnetic resonance imaging (MRI) with the ability to delivery energy to targets deep in the body noninvasively with focused ultrasound presents a disruptive technology with the potential to significantly affect healthcare. MRI offers precise targeting, visualization, and quantification of temperature changes and the ability to immediately evaluate the treatment. By exploiting different mechanisms, focused ultrasound offers a range of therapies, ranging from thermal ablation to targeted drug delivery. This article reviews recent preclinical and tests clinical of this technology.

Feasibility of transrib focused ultrasound thermal ablation for liver tumors using a spherically curved 2D array: a numerical study

The use of focused ultrasound thermal ablation to treat hepatocarcinoma and other liver tumors produces promising clinical results. However, one of the major drawbacks is the high absorption of ultrasonic energy by the rib, making partial rib removal necessary in many cases. This study numerically investigated the feasibility of using a spherical ultrasound phased array for transrib liver-tumor thermal ablation. An independently array-element activitation scheme, which switches off the transducer elements obstructed by the ribs based on feedback anatomical medical imaging, was proposed to reduce the rib-overheating problem. The numerical results showed that the proposed treatment planning strategy can effectively reduce the specific energy absorbed by the rib while maintaining the energy at the target position, which both reduces the rib-overheating problem and increases the possibility of treating a target lesion under an intact rib. The analysis also demonstrated that the target position and the ultrasound frequency play key roles in the treatment. Patients with diverse characteristics were also tested to show the generality of the proposed strategy. The proposed treatment planning strategy also provides useful information for evaluating the treatment effectiveness prior to clinically performing transrib ultrasound liver-tumor thermal ablation.

A fast and conformal heating scheme for producing large thermal lesions using a 2D ultrasound phased array

The treatment conformability and the total treatment time of large tumors are both important issues in ultrasound thermal therapy. Previous heating strategies all show their restrictions in achieving these two issues to satisfactory levels simultaneously. This work theoretically presents a new heating strategy which is capable of both increasing the treatment conformability and shortening the treatment time, when using a 2D ultrasound phased array transducer. To perform this, a set of the multiple-foci patterns (considered the basic heating units) were temporally switched to steer the beam at different focal planes with the lesion length being well-controlled. Then, to conformally cover an irregular target volume, the 2D phased array was laterally shifted by a positioning system to deposit a suitable heating unit to cover a subvolume part. Results demonstrated that the totally treatment time can be largely reduced. The heating rate can be increased up to 0.96 cm3/min compared to the previously reported 0.26 cm3/min. Also, the proposed scheme showed that the tumor regions can be completely treated with the normal tissue damage at satisfactory level. The feasibility of the proposed strategy for irregular tumor treatment was also demonstrated. This study offers useful information in large tumor treatment in ultrasound thermal therapy.

A novel strategy to increase heating efficiency in a split-focus ultrasound phased array

Focus splitting using sector-based phased arrays increases the necrosed volume in a single sonication and reduces the total treatment time in the treatment of large tumors. However, split-focus sonication results in a lower energy density and worse focal-beam distortion, which limits its usefulness in practical treatments. Here, we propose a new heating strategy involving consecutive strongly focused and split-focus sonications to improve the heating efficiency. Theoretical predictions including linear and thermal-dose-dependent attenuation change were employed to investigate potential factors of this strategy, and ex vivo tissue experiments were conducted to confirm its effectiveness. Results showed that the thermal lesions produced by the proposed strategy could be increased when comparing with the previous reported strategies. The proposed heating strategy also induces a thermal lesion more rapidly, and exhibits higher robustness to various blood perfusion conditions, higher robustness to various power/heating time combinations, and superiority to generate deep-seated lesions through tissues with complex interfaces. Possible mechanisms include the optimization of the thermal conduction created by the strongly focused sonication and the temperature buildup gained from thermally induced tissue attenuation change based on the theoretical analysis. This may represent a useful technique for increasing the applicability of split-focus and multiple-focus sonication techniques, and solve the obstacles encountered when attempting to use these methods to shorten the total clinical treatment time.

Image-based control of the magnetic resonance imaging-guided focused ultrasound thermotherapy

Magnetic resonance imaging (MRI)-guided focused ultrasound surgery (FUS) is a full noninvasive approach for localized thermal ablation of deep tissues, coupling the following: (1) a versatile, nonionizing physical agent for therapy and (2) a state-of-the art diagnosis and on-line monitoring tool. A commercially available, Food and Drug Administration-approved device using the MRI-guided FUS exists since 2004 for the ablation of benign tumors (uterine fibroids); however, the ultimate goal of the technological, methodological, and medical research in this field is to provide a clinical-routine tool for fighting localized cancer. When addressing cancer applications, the accurate spatial control of the delivered thermal dose is mandatory. Contiguous destruction of the target volume must be achieved in a minimum time, whereas sparing as much as possible the neighboring healthy tissues and especially when some adjacent regions are critical. This paper reviews some significant developments reported in the literature related to the image-based control of the FUS therapy for kidney, breast, prostate, and brain, including the own experience of the authors on the active feedback control of the temperature during FUS ablation. In addition, preliminary results of an original study of MRI-guided FUS ablation of VX2 carcinoma in kidney, under active temperature control, are described here.

A real-time measure of cavitation induced tissue disruption by ultrasound imaging backscatter reduction

A feedback method for obtaining real-time information on the mechanical disruption of tissue through ultrasound cavitation is presented. This method is based on a substantial reduction in ultrasound imaging backscatter from the target volume as the tissue structure is broken down. Ex-vivo samples of porcine liver were exposed to successive high-intensity ultrasound pulses at a low duty cycle to induce mechanical disruption of tissue parenchyma through cavitation (referred to as histotripsy). At the conclusion of treatment, B-scan imaging backscatter was observed to have decreased by 22.4 +/- 2.3 dB in the target location. Treated samples of tissue were found to contain disrupted tissue corresponding to the imaged hypoechoic volume with no remaining discernable structure and a sharp boundary. The observed, substantial backscatter reduction may be an effective feedback mechanism for assessing treatment efficacy in ultrasound surgery using pulsed ultrasound to create cavitation.

Magnetic resonance-guided focused ultrasound surgery for uterine fibroids: relationship between the therapeutic effects and signal intensity of preexisting T2-weighted magnetic resonance images

OBJECTIVE

This study was undertaken to clarify the relationship between the signal intensity of T2-weighted magnetic resonance images and the therapeutic effect of magnetic resonance-guided focused ultrasound surgery (MRgFUS) on uterine fibroids.

STUDY DESIGN

Ninety-five fibroids in 63 patients were classified into 3 types based on the signal intensity of T2-weighted magnetic resonance images as follows: type 1, low intensity; type 2, intermediate intensity; type 3, high intensity. The treated area ratio of MRgFUS and the volume reduction ratio 6 months after treatment were used as the indices of therapeutic effect.

RESULTS

The treated area ratio of type 3 fibroids was the lowest among the 3 types (P < .01). The volume reduction ratio correlated with the treated area ratio (r = 0.64; P < .01).

CONCLUSION

The efficacy of MRgFUS correlates with the signal intensity of T2-weighted magnetic resonance images. Type 1 and type 2 fibroids are suitable candidates for MRgFUS, whereas type 3 fibroids are not.

Uterine leiomyomas: MR imaging-based thermometry and thermal dosimetry during focused ultrasound thermal ablation

PURPOSE

To retrospectively evaluate magnetic resonance (MR) imaging-based thermometry and thermal dosimetry during focused ultrasound treatments of uterine leiomyomas (ie, fibroids).

MATERIALS AND METHODS

All patients gave written informed consent for the focused ultrasound treatments and the current HIPAA-compliant retrospective study, both of which were institutional review board approved. Thermometry performed during the treatments of 64 fibroids in 50 women (mean age, 46.6 years +/- 4.5 [standard deviation]) was used to create thermal dose maps. The areas that reached dose values of 240 and 18 equivalent minutes at 43 degrees C were compared with the nonperfused regions measured on contrast material-enhanced MR images by using the Bland-Altman method. Volume changes in treated fibroids after 6 months were compared with volume changes in nontreated fibroids and with MR-based thermal dose estimates.

RESULTS

While the thermal dose estimates were shown to have a clear relationship with resulting nonperfused regions, the nonperfused areas were, on average, larger than the dose estimates (means of 1.9 +/- 0.7 and 1.2 +/- 0.4 times as large for areas that reached 240- and 18-minute threshold dose values, respectively). Good correlation was observed for smaller treatment volumes at the lower dose threshold (mean ratio, 1.0 +/- 0.3), but for larger treatment volumes, the nonperfused region extended to locations within the fibroid that clearly were not heated. Variations in peak temperature increase were as large as a factor of two, both between patients and within individual treatments. On average, the fibroid volume reduction at 6 months increased as the ablated volume estimated by using the thermal dose increased.

CONCLUSION

Study results showed good correlation between thermal dose estimates and resulting nonperfused areas for smaller ablated volumes. For larger treatment volumes, nonperfused areas could extend within the fibroid to unheated areas.

Magnetic Resonance–Guided Focused Ultrasound Surgery of Breast Cancer: Reliability and Effectiveness

BACKGROUND

Magnetic resonance-guided focused ultrasound surgery (MRgFUS) is a noninvasive technique that has been shown to coagulate benign and malignant tumors. The purpose of this study was to evaluate MRgFUS safety and effectiveness for the ablation of breast carcinomas.

STUDY DESIGN

Thirty women with biopsy-proved breast cancer underwent MRgFUS treatment. Gadolinium-enhanced MR images were used for treatment planning and posttreatment radiologic assessment of treated tissue, and temperature-sensitive MR images provided real-time treatment monitoring. After MRgFUS, all 30 women underwent wide excision or mastectomy. The extent of thermal ablation was assessed with tumor histology.

RESULTS

Treatment was well tolerated, with a minimum of adverse effects, especially when performed under local anesthesia. On pathologic examination, mean (+/-SD) necrosis of the targeted breast tumors was 96.9 +/- 4% (median 100%, range 78% to 100%) of tumor volume. Fifteen (53.5%) of 28 evaluable patients had 100% necrosis of the ablated tumor; only 3 patients (10.7%) had less than 95% necrosis. In 28 (93.3%) patients, 100% of the malignancy was within the treatment field, and 98% and 95% of tumor lay within the treatment field in 2 remaining patients. Retrospective analysis in two patients with residual tumor showed treatment was not delivered to the full recommended area, reaffirming the need for precise localization and the value of contrast-enhanced images for treatment planning.

CONCLUSIONS

MRgFUS has great potential to become a viable noninvasive replacement for lumpectomy. Additional studies focusing on posttreatment image-based evaluation are needed.

Annular phased-array high-intensity focused ultrasound device for image-guided therapy of uterine fibroids

An ultrasound (US), image-guided high-intensity focused ultrasound (HIFU) device was developed for noninvasive ablation of uterine fibroids. The HIFU device was an annular phased array, with a focal depth range of 30-60 mm, a natural focus of 50 mm, and a resonant frequency of 3 MHz. The in-house control software was developed to operate the HIFU electronics drive system for inducing tissue coagulation at different distances from the array. A novel imaging algorithm was developed to minimize the HIFU-induced noise in the US images. The device was able to produce lesions in bovine serum albumin-embedded polyacrylamide gels and excised pig liver. The lesions could be seen on the US images as hyperechoic regions. Depths ranging from 30 to 60 mm were sonicated at acoustic intensities of 4100 and 6100 W/cm2 for 15 s each, with the latter producing average lesion volumes at least 63% larger than the former. Tissue sonication patterns that began distal to the transducer produced longer lesions than those that began proximally. The variation in lesion dimensions indicates the possible development of HIFU protocols that increase HIFU throughput and shorten tumor treatment times.

Heating efficiency improvement by using a spherically-concaved sectored array in focused ultrasound thermal therapy

Focus splitting by using sector-sectioned phased arrays is one of effective methods to increase the necrosed volume in single sonication and to reduce the total treatment time in large tumor treatment. However, the split focus contains less concentrated energy and severer focal beam distortion, which limits its usefulness in practical treatments. In this study, we proposed a new heating strategy by combining sonications of strongly-focused and split-focused patterns to increase the heating efficiency. Theoretical predictions and ex-vivo tissue experiments showed that thermal lesions can be enlarged in dimensions after applying the proposed strategy. This may provide a useful way to solve current obstacles in low heating efficiency of split-focus sonications that attempted to shorten the total treatment time in current clinical application.