Characterizing pore-scale structure-flow correlations in sedimentary rocks using magnetic resonance imaging.

Quantitative, three-dimensional (3D) spatially resolved magnetic resonance flow imaging (flow MRI) methods are presented to characterize structure-flow correlations in a 4-mm-diameter plug of Ketton limestone rock using undersampled k- and q-space data acquisition methods combined with compressed sensing (CS) data reconstruction techniques. The acquired MRI data are coregistered with an X-ray microcomputed tomography (µCT) image of the same rock sample, allowing direct correlation of the structural features of the rock with local fluid transport characteristics. First, 3D velocity maps acquired at 35 µm isotropic spatial resolution showed that the flow was highly heterogeneous, with ∼10% of the pores carrying more than 50% of the flow. Structure-flow correlations were found between the local flow velocities through pores and the size and topology (coordination number) associated with these pores. These data show consistent trends with analogous data acquired for flow through a packing of 4-mm-diameter spheres, which may be due to the microstructure of Ketton rock being a consolidation of approximately spherical grains. Using two-dimensional and 3D visualization of coregistered µCT images and velocity maps, complex pore-scale flow patterns were identified. Second, 3D spatially resolved propagators were acquired at 94 µm isotropic spatial resolution. Flow dispersion within the rock was examined by analyzing each of the 331 776 local propagators as a function of observation time. Again, the heterogeneity of flow within the rock was shown. Quantification of the mean and standard deviation of each of the local propagators showed enhanced mixing occurring within the pore space at longer observation times. These spatially resolved measurements also enable investigation of the length scale of a representative elementary volume. It is shown that for a 4-mm-diameter plug this length scale is not reached.

Identification of sampling patterns for high-resolution compressed sensing MRI of porous materials: 'learning' from X-ray microcomputed tomography data.

There exists a strong motivation to increase the spatial resolution of magnetic resonance imaging (MRI) acquisitions so that MRI can be used as a microscopy technique in the study of porous materials. This work introduces a method for identifying novel data sampling patterns to achieve undersampling schemes for compressed sensing MRI (CS-MRI) acquisitions, enabling 3D spatial resolutions of 17.6 µm to be achieved. A data-driven learning approach is used to derive k-space undersampling schemes for 3D MRI acquisitions from 3D X-ray microcomputed tomography (µCT) datasets acquired at a higher spatial resolution than can be acquired using MRI. The performance of the new sampling approach was compared to other, well-established sampling strategies using simulated MRI data obtained from high-resolution µCT images of rock core plugs. These simulations were performed for a range of different k-space sampling fractions (0.125-0.375) using images of Ketton limestone. The method was then extended to consideration of imaging Estaillades limestone and Fontainebleau sandstone. The results show that the new sampling approach performs as well as or better than conventional variable density sampling and without need for time-consuming parameter optimisation. Further, a bespoke sampling pattern is produced for each rock type. The novel undersampling strategy was employed to acquire 3D magnetic resonance images of a Ketton limestone rock at spatial resolutions of 35 and 17.6 µm. The ability of the k-space sampling scheme produced using the new approach in enabling reconstruction of the pore space characteristics of the rock was then demonstrated by benchmarking against the pore space statistics obtained from high-resolution µCT data. The MRI data acquired at 17.6 µm resolution gave excellent agreement with the pore size distribution obtained from the X-ray microcomputed tomography dataset, while the pore coordination number distribution obtained from the MRI data was slightly skewed to lower coordination numbers. This approach provides a method of producing a k-space undersampling pattern for MRI acquisition at a spatial resolution for which a fully sampled acquisition at that spatial resolution would be impractically long. The approach can be easily extended to other CS-MRI techniques, such as spatially resolved flow and relaxation time mapping. LAY DESCRIPTION: Magnetic resonance imaging (MRI) is widely used to study the microstructure of, and fluid transport phenomena in porous media relevant for engineering applications. A major application is the study of water and hydrocarbon transport in porous sedimentary rocks, which typically have pore sizes smaller than 100 µm. The spatial resolution of routine MRI acquisitions, however, is limited to several hundred µm due to the relatively low sensitivity of the magnetic resonance method. Therefore, there exists a strong motivation to increase the spatial resolution of MRI by one to two orders of magnitude to be able to study these rocks at a pore scale. This work reports the initial step towards achieving this. Three-dimensional images of rock pore structure are acquired at both 35 and 17.6 µm spatial resolution. In ongoing work, these methods are now being incorporated into magnetic resonance velocity imaging methods, thereby enabling imaging of both pore structure and hydrodynamics at these much higher spatial resolutions than were hitherto possible. Although X-ray microcomputed tomography (µCT) produces high spatial resolution images, it is far more limited in being able to spatially map transport processes (i.e. flow) in porous media. This work reports a strategy for accelerating the image acquisition time such that sufficient signal-to-noise ratio (SNR) is achieved to increase the spatial resolution, that is, the voxel size within which there is sufficient SNR within the resulting image. To achieve this, a technique known as compressed sensing is used which exploits undersampling of the acquired data relative to the standard fully sampled image. In MRI, data are acquired in so-called k-space and Fourier transformed to yield the real space image. The challenge, when undersampling, is to optimise the specific points in k-space that are acquired because these will influence the quality of the resulting image. This work reports a straightforward, robust strategy for identifying the optimal sets of k-space points to acquire. The method introduced uses simulated MRI images calculated from high-resolution µCT images of the rocks of interest, from which optimised MRI sampling patterns are obtained. The method does not require any optimisation of parameters for its implementation, which is a significant advantage compared to other strategies. Moreover, we show that the pore space characteristics of the acquired MRI images are in excellent agreement with the same characteristics obtained from a high-resolution µCT image.

Accelerating the estimation of 3D spatially resolved T2 distributions.

Obtaining quantitative, 3D spatially-resolved T2 distributions (T2 maps) from magnetic resonance data is of importance in both medical and porous media applications. Due to the long acquisition time, there is considerable interest in accelerating the experiments by applying undersampling schemes during the acquisition and developing reconstruction techniques for obtaining the 3D T2 maps from the undersampled data. A multi-echo spin echo pulse sequence is used in this work to acquire the undersampled data according to two different sampling patterns: a conventional coherent sampling pattern where the same set of lines in k-space is sampled for all equally-spaced echoes in the echo train, and a proposed incoherent sampling pattern where an independent set of k-space lines is sampled for each echo. The conventional reconstruction technique of total variation regularization is compared to the more recent techniques of nuclear norm regularization and Nuclear Total Generalized Variation (NTGV) regularization. It is shown that best reconstructions are obtained when the data acquired using an incoherent sampling scheme are processed using NTGV regularization. Using an incoherent sampling pattern and NTGV regularization as the reconstruction technique, quantitative results are obtained at sampling percentages as low as 3.1% of k-space, corresponding to a 32-fold decrease in the acquisition time, compared to a fully sampled dataset.

Relationship between neighborhood socioeconomic status and venous thromboembolism: results from a population-based study.

Essentials Literature on socioeconomic status (SES) and incidence of venous thromboembolism (VTE) is scarce. We assessed neighborhood SES with VTE risk in a population of over 1.4 million inhabitants. Higher neighborhood SES was associated with lower incidence of VTE. These findings are helpful to inform policy and resource allocation in health systems. SUMMARY: Background The association between socioeconomic status and arterial cardiovascular disease is well established. However, despite its high burden of disability-adjusted life years, little research has been carried out to determine whether socioeconomic status is associated with venous thromboembolism. Objective To determine if neighborhood socioeconomic status is associated with venous thromboembolism in a population-based study from the Netherlands. Methods We identified all patients aged 15 years and older with a first event of venous thromboembolism from inhabitants who lived in the urban districts of The Hague, Leiden and Utrecht in the Netherlands in 2008-2012. Neighborhood socioeconomic status was based on the status score, which combines educational level, income and unemployment on a four-digit postal code level. Incidence rate ratios of venous thromboembolism were calculated for different levels of neighborhood socioeconomic status, with adjustments for age and sex. Results A total of 7373 patients with a first venous thromboembolism (median age 61 years; 50% deep vein thrombosis) were identified among more than 1.4 million inhabitants. Higher neighborhood SES was associated with lower incidence of VTE. In the two highest status score groups (i.e. the 95-99th and > 99th percentile), the adjusted incidence rate ratios were 0.91 (95% confidence interval [CI], 0.84-1.00) and 0.80 (95% CI, 0.69-0.93), respectively, compared with the reference status score group (i.e. 30-70th percentile). Conclusions High neighborhood socioeconomic status is associated with a lower risk of first venous thromboembolism.

Management of stab wounds to the pregnant uterus: a case report and a review of the literature.

Only 18 cases of stab wounds to the pregnant uterus have been reported to date. This paper presents an additional case--a patient at 26 weeks' gestation with multiple voluntary stab wounds to the abdomen and lower chest resulting in fetal death. The literature is reviewed and management of abdominal stab wounds in pregnancy is discussed.