Harmonization of PET image reconstruction parameters in simultaneous PET/MRI

Objective Simultaneous PET/MRIs vary in their quantitative PET performance due to inherent differences in the physical systems and differences in the image reconstruction implementation. This variability in quantitative accuracy confounds the ability to meaningfully combine and compare data across scanners. In this work, we define image reconstruction parameters that lead to comparable contrast recovery curves across simultaneous PET/MRI systems. Method The NEMA NU-2 image quality phantom was imaged on one GE Signa and on one Siemens mMR PET/MRI scanner. The phantom was imaged at 9.7:1 contrast with standard spheres (diameter 10, 13, 17, 22, 28, 37 mm) and with custom spheres (diameter: 8.5, 11.5, 15, 25, 32.5, 44 mm) using a standardized methodology. Analysis was performed on a 30 min listmode data acquisition and on 6 realizations of 5 min from the listmode data. Images were reconstructed with the manufacturer provided iterative image reconstruction algorithms with and without point spread function (PSF) modeling. For both scanners, a post-reconstruction Gaussian filter of 3–7 mm in steps of 1 mm was applied. Attenuation correction was provided from a scaled computed tomography (CT) image of the phantom registered to the MR-based attenuation images and verified to align on the non-attenuation corrected PET images. For each of these image reconstruction parameter sets, contrast recovery coefficients (CRCs) were determined for the SUVmean, SUVmax and SUVpeak for each sphere. A hybrid metric combining the root-mean-squared discrepancy (RMSD) and the absolute CRC values was used to simultaneously optimize for best match in CRC between the two scanners while simultaneously weighting toward higher resolution reconstructions. The image reconstruction parameter set was identified as the best candidate reconstruction for each vendor for harmonized PET image reconstruction. Results The range of clinically relevant image reconstruction parameters demonstrated widely different quantitative performance across cameras. The best match of CRC curves was obtained at the lowest RMSD values with: for CRCmean, 2 iterations-7 mm filter on the GE Signa and 4 iterations-6 mm filter on the Siemens mMR, for CRCmax, 4 iterations-6 mm filter on the GE Signa, 4 iterations-5 mm filter on the Siemens mMR and for CRCpeak, 4 iterations-7 mm filter with PSF on the GE Signa and 4 iterations-7 mm filter on the Siemens mMR. Over all reconstructions, the RMSD between CRCs was 1.8%, 3.6% and 2.9% for CRC mean, max and peak, respectively. The solution of 2 iterations-3 mm on the GE Signa and 4 iterations-3 mm on Siemens mMR, both with PSF, led to simultaneous harmonization and with high CRC and low RMSD for CRC mean, max and peak with RMSD values of 2.8%, 5.8% and 3.2%, respectively. Conclusions For two commercially available PET/MRI scanners, user-selectable parameters that control iterative updates, image smoothing and PSF modeling provide a range of contrast recovery curves that allow harmonization in harmonization strategies of optimal match in CRC or high CRC values. This work demonstrates that nearly identical CRC curves can be obtained on different commercially available scanners by selecting appropriate image reconstruction parameters.

Modified Method of Characteristics for Generating EOR Oil Recovery Curves

The paper describes a fast and approximate 1D simulation algorithm for calculating the percent recovery that can be obtained from an oil reservoir if gas injection is carried out at a pressure lower than the minimum miscibility pressure. The algorithm is based on the Method of Characteristics. While a conventional 1D reservoir simulation of a gas injection scenario may take minutes or even hours, the proposed algorithm allows a full evaluation of the recovery to be completed within seconds. To make the method numerically robust, a number of approximations were needed. The result is an extremely fast algorithm that not only provides a good estimate of the recovery obtained by gas injection, but also gives a good visualization of how the gas displaces the oil.

Exposure Stability of Quantitative Fluorine-18 Measurements On The Biograph Vision Quadra PET/CT System

Background: Our aim was to determine sets of reconstruction parameters for the Biograph Vision Quadra (Siemens Healthineers) PET/CT system that result in quantitative images compliant with the European Association of Nuclear Medicine Research Ltd. (EARL) criteria. Using the Biograph Vision 600 (Siemens Healthineers) PET/CT technology but extending the axial field of view to 106 cm, gives the Vision Quadra currently an around fivefold higher sensitivity over the Vision 600 with otherwise comparable spatial resolution. Therefore, we also investigated how the number of incident positron decays - i.e. exposure - affects EARL compliance. This will allow estimating a minimal acquisition time or a minimal applied dose in clinical scans while retaining data comparability. Methods: We measured activity recovery curves on a NEMA IEC body phantom filled with an aqueous 18 F solution and a sphere to background ratio of 10 to 1 according to the latest EARL guidelines. Reconstructing 3570 images with varying OSEM PSF iterations, post-reconstruction Gaussian filter full width at half maximum (FWHM), and varying exposure from 0.2 MDecays/ml (= 10 sec frame duration) to 59.2 MDecays/ml (= 1 h frame duration), allowed us to determine sets of parameters to achieve compliance with the current EARL 1 and EARL 2 standards. Recovery coefficients (RCs) were calculated for the metrics RC max , RC mean , and RC peak , and the respective recovery curves were analysed for monotonicity. Results: Using 6 iterations, 5 subsets and 7.8 mm Gauss filtering resulted in best EARL 1 compliance and recovery curve monotonicity in all analysed frames. Most robust EARL2 compliance and monotonicity was achieved with 4 iterations, 5 subsets, and 4.6 mm Gauss FWHM in frames with durations between 10 min and 30 sec. RC peak only impeded EARL2 compliance in the 10 sec frame. Conclusions: While EARL1 compliance proved to be robust over all exposure ranges, EARL2 compliance required exposures between 0.6 MDecays/ml to 11.5 MDecays/ml. The Biograph Vision Quadra’s high sensitivity makes frames as short as 10 sec feasible for comparable quantitative images. Lowering EARL2 RC max limits closer to unity would possibly even permit 10 sec EARL2 compliant frames.

Evaluating Burn Recovery Outcomes in Children with Neurodevelopmental Symptoms

Neurodevelopmental symptoms (NS) including attention and behavioral problems, developmental delays, intellectual disabilities and learning problems are prevalent in children with burn injuries. The presence of NS may predispose children to poorer burn injury recovery outcomes compared to children without these symptoms (non-NS). The Multi-Center Benchmarking Study (MCBS) monitored recovery outcomes in children with burn injuries in real time utilizing the Burn Outcomes Questionnaire (BOQ). The objective of this study was to retrospectively assess the long term burn recovery outcomes in NS patients versus non-NS patients from the MCBS population. This study assessed parent-reported BOQ outcomes in a sample of 563 patients aged 5 to 18 years who were admitted for burn injuries to a pediatric burn center. A subsample of patients had reported NS (n=181). Analyses compared BOQ outcomes within the NS subsample and the non-NS subsample (n=382) across three longitudinal points post-discharge. The prevalence rate of NS was 32.1% in the full sample. Findings revealed statistically significant improvement in the recovery curves in all five BOQ subscales for the non-NS group and all subscales except for Compliance for the NS group across all longitudinal points. When compared to non-NS patients, NS patients had significantly poorer burn recovery outcomes on the Satisfaction and Compliance subscales. Although it is important to educate all clinicians, parents, and children on burn prevention efforts, targeted education is necessary for children with NS since they may be at greater risk for burn injury as well as worse recovery outcomes.

Evaluation of Iodine-123 and Iodine-131 SPECT activity quantification: a Monte Carlo study

Purpose The quantitative accuracy of Nuclear Medicine images, acquired for both planar and SPECT studies, is influenced by the isotope-collimator combination as well as image corrections incorporated in the iterative reconstruction process. These factors can be investigated and optimised using Monte Carlo simulations. This study aimed to evaluate SPECT quantification accuracy for 123I with both the low-energy high resolution (LEHR) and medium-energy (ME) collimators and 131I with the high-energy (HE) collimator. Methods Simulated SPECT projection images were reconstructed using the OS-EM iterative algorithm, which was optimised for the number of updates, with appropriate corrections for scatter, attenuation and collimator detector response (CDR), including septal scatter and penetration compensation. An appropriate calibration factor (CF) was determined from four different source geometries (activity-filled: water-filled cylindrical phantom, sphere in water-filled (cold) cylindrical phantom, sphere in air and point-like source), investigated with different volume of interest (VOI) diameters. Recovery curves were constructed from recovery coefficients to correct for partial volume effects (PVEs). The quantitative method was evaluated for spheres in voxel-based digital cylindrical and patient phantoms. Results The optimal number of OS-EM updates was 60 for all isotope-collimator combinations. The CFpoint with a VOI diameter equal to the physical size plus a 3.0-cm margin was selected, for all isotope-collimator geometries. The spheres’ quantification errors in the voxel-based digital cylindrical and patient phantoms were less than 3.2% and 5.4%, respectively, for all isotope-collimator combinations. Conclusion The study showed that quantification errors of less than 6.0% could be attained, for all isotope-collimator combinations, if corrections for; scatter, attenuation, CDR (including septal scatter and penetration) and PVEs are performed. 123I LEHR and 123I ME quantification accuracies compared well when appropriate corrections for septal scatter and penetration were applied. This can be useful in departments that perform 123I studies and may not have access to ME collimators.

Interstitial and cerebrospinal fluid exchanging process revealed by phase alternate labeling with null recovery MRI

Purpose: To develop Phase Alternate LAbeling with Null recovery (PALAN) MRI methods for the quantification of interstitial to cerebrospinal fluid flow (ICF) and cerebrospinal to interstitial fluid flow (CIF) in the brain. Method: In both T1-PALAN and apparent diffusion coefficient (ADC)-PALAN MRI methods, the cerebrospinal fluid (CSF) signal was nulled, while the residual interstitial fluid (ISF) was labeled by alternating the phase of pulses. ICF was extracted from the difference between the recovery curves of CSF with and without labeling. Similarly, CIF was measured by the T2-PALAN MRI method by labeling CSF, which took advance of the significant T2 difference between CSF and parenchyma. Results: Both T1-PALAN and ADC-PALAN observed a rapid occurrence of ICF at 67±56 ms and 13±2 ms interstitial fluid transit times, respectively. ICF signal peaked at 1.5 s for both methods. ICF was 1153±270 ml/100ml/min with T1-PALAN in the third and lateral ventricles, which was higher than 891±60 ml/100ml/min obtained by ADC-PALAN. The results of the T2-PALAN suggested the ISF exchanging from ependymal layer to the parenchyma was extremely slow. Conclusion: The PALAN methods are suitable tools to study ISF and CSF flow kinetics in the brain.

Oil recovery enhancement by reservoir hydraulic compression technique employment

Industrial experiment works (IEW) were carried out to study the mechanism of filtration and reservoir properties changes (FRP) in the process of wells swabbing. Based on the hydrodynamic studies, the results of the works are analyzed. A method for oil production enhancing by reservoirs hydraulic compression has been worked out. In the process of well swabbing the barograms were recorded, pressure recovery curves were taken with the determination of hydraulic conductivity and piezoconductivity values, potential productivity coefficients, well flow rate, reservoir pressure before and after exposure. The interpretation of hydrodynamic studies was carried out by the deterministic analysis with subsequent modeling of the situation. The reservoir, opened by the perforation interval, is of complex structure, as a result of which the liquid was absorbed by the interlayer located above the area with newly formed microcracks. Keywords: hard-to-recover reserves; swabbing; carbonate reservoirs; filtration reservoir properties; pressure recovery curve.

Evaluation of Iodine-123 and Iodine-131 SPECT Activity Quantification: A Monte Carlo Study

Purpose: The quantitative accuracy of Nuclear Medicine images, acquired for both planar and SPECT studies, is influenced by the isotope-collimator combination as well as image corrections incorporated in the iterative reconstruction process. These factors can be investigated and optimised using Monte Carlo simulations. This study aimed to evaluate SPECT quantification accuracy for 123I with both the low energy high resolution (LEHR) and medium energy (ME) collimators, and 131I with the high energy (HE) collimator. Methods: Simulated SPECT projection images were reconstructed using the OS-EM iterative algorithm, which was optimised for the number of updates, with appropriate corrections for scatter, attenuation, and collimator detector response (CDR), including septal scatter and penetration compensation. An appropriate conversion factor (CF) was determined from four different source geometries (activity-filled: water-filled cylindrical phantom, sphere in water-filled (cold) cylindrical phantom, sphere in air and point-like source), investigated with different VOI diameters. Recovery curves were constructed from recovery coefficients to correct for partial volume effects (PVEs). The quantitative method was evaluated for spheres in voxel-based digital cylindrical and patient phantoms. Results: The optimal number of OS-EM updates was 60 for all isotope-collimator combinations. The CFpoint with a VOI diameter equal to the physical size plus a 3.0 cm margin was selected, for all isotope-collimator geometries. The spheres’ quantification errors in the voxel-based digital cylindrical and patient phantoms were less than 3.2% and 5.4%, respectively, for all isotope-collimator combinations. Conclusion: The study showed that quantification errors of less than 6.0% could be attained, for all isotope-collimator combinations, if corrections for; scatter, attenuation, CDR (including septal scatter and penetration) and PVEs, are performed. 123I LEHR and 123I ME quantification accuracies compared well when appropriate corrections for septal scatter and penetration were applied. This can be useful in departments that perform 123I studies and may not have access to ME collimators.

Effect of viscosity on oil recovery by spontaneous imbibition from partially water-covered matrix blocks

Spontaneous imbibition is an important mechanism of oil recovery from fractured reservoirs and unconventional reservoirs. Oil is produced by combining co- and counter-current imbibition when the matrix blocks was partially covered by water. In this paper, we focused on the effect of viscosity ratios on oil production by spontaneous imbibition and established the numerical model for one-dimensional linear imbibition with TEO-OW boundary conditions, which was validated by the experimental data. The effect of viscosity ratio on co- and counter-current imbibition is investigated and scaling result of the imbibition recovery curve for wide range of viscosity ratio using the conventional scaling equation was tested, which indicates that the close correlation was achieved only when oil-water viscosity ratios are higher. Then, a modified scaling equation was developed based on the piston-like assumption for one-dimensional co-current imbibition and close correlation of imbibition recovery curves was achieved when viscosity ratios are lower. Finally, correlation of imbibition recovery curves was improved for wide range of viscosity ratios by combining conventional and modified scaling equation. Results show that since the shape of imbibition recovery curves is not similar for different viscosity ratios, it is difficult to obtain the perfect correlation using the constant viscosity term.

APPLICABILITY INDICATORS AND IDENTIFICATION TECHNIQUE FOR A NONLINEAR MAXWELL-TYPE ELASTO-VISCOPLASTIC MODEL USING REPEATED CREEP RECOVERY TESTS

A physically non-linear Maxwell-type constitutive relation for non-aging elasto-viscoplastic materials is considered. General properties of repeated creep recovery curves generated by the relation under four-step uni-axial loadings consisting of two rectangular pulses of stress and two rest periods at zero stress are studied analytically assuming two material functions of the relation are arbitrary. The analysis reveals several characteristic features of the theoretic creep recovery curves and repeated creep recovery curves that can be employed as the relation applicability (or non- applicability) indicators which are convenient for check using test data of a material. Two effective general identification techniques are developed. The first one is based on a set of creep recovery tests at various stress levels and implies measurement of two strain magnitudes in each test. The second one is based on a set of repeated creep and recovery tests and implies measurement of four strain magnitudes in each test. More complex loading program enables to reduce the number of tests and samples twice and to obtain more information on a material behavior and to check a number of additional applicability indicators. The explicit expressions are derived in each case to determine the material functions values at arbitrarily chosen points in stress domain. The identification techniques enable separate and direct evaluation of the material functions values via test data without error accumulation. A number of the identification technique versions are considered and their advantages and shortcomings are discussed. In the case of materials exhibiting creep rate power dependence on stress, a specific rapid procedure is developed for the model identification in the class of power material functions and additional applicability indicators are found. In this case only one repeated creep recovery test and four measured magnitudes of strain are sufficient to determine four material parameters through the explicit expressions obtained.