Accurate phase-shift velocimetry in rock

Optimal control flow encoding for time-efficient magnetic resonance velocimetry

Phase contrast velocimetry relies on bipolar gradients to establish a direct and linear relationship between the phase of the magnetic resonance signal, and the corresponding fluid motion. Despite its utility, several limitations and drawbacks have been reported, the most important being the extended echo time due to the encoding after the excitation. In this study, we elucidate a new approach based on optimal control theory that circumvents some of these disadvantages. An excitation pulse, termed FAUCET (flow analysis under controlled encoding transients), is designed to encode velocity into phase already during the radiofrequency excitation. As a result of concurrent excitation and flow encoding, and hence elimination of post-excitation flow encoding, FAUCET achieves a shorter echo time than the conventional method. This achievement is a matter of significance not only because it decreases the loss of signal due to spin-spin relaxation and B₀ inhomogeneity, but also because a shorter echo time is always preferred in order to reduce the dimensionless dephasing parameter and the required residence time of the flowing sample in the detection coil. The method is able to establish a non-linear bijective relationship between phase and velocity, which can be employed to enhance the resolution over a specific range of velocities, for example along flow boundaries. A computational comparison between the phase contrast and optimal control methods reveals that the latter's encoding is more robust against remnant higher-order-moment terms of the Taylor expansion for faster voxels, such as acceleration, jerk, and snap.

Integrating Pore-Scale Flow MRI and X-ray μCT for Validation of Numerical Flow Simulations in Porous Sedimentary Rocks

Single-phase fluid flow velocity maps in Ketton and Estaillades carbonate rock core plugs are computed at a pore scale, using the lattice Boltzmann method (LBM) simulations performed directly on three-dimensional (3D) X-ray micro-computed tomography (µCT) images (≤ 7 µm spatial resolution) of the core plugs. The simulations are then benchmarked on a voxel-by-voxel and pore-by-pore basis to quantitative, 3D spatially resolved magnetic resonance imaging (MRI) flow velocity maps, acquired at 35 µm isotropic spatial resolution for flow of water through the same rock samples. Co-registration of the 3D experimental and simulated velocity maps and coarse-graining of the simulation to the same resolution as the experimental data allowed the data to be directly compared. First, the results are demonstrated for Ketton limestone rock, for which good qualitative and quantitative agreement was found between the simulated and experimental velocity maps. The flow-carrying microstructural features in Ketton rock are mostly larger than the spatial resolution of the µCT images, so that the segmented images are an adequate representation of the pore space. Second, the flow data are presented for Estaillades limestone, which presents a more heterogeneous case with microstructural features below the spatial resolution of the µCT images. Still, many of the complex flow patterns were qualitatively reproduced by the LBM simulation in this rock, although in some pores, noticeable differences between the LBM and MRI velocity maps were observed. It was shown that 80% of the flow (fractional summed z-velocities within pores) in the Estaillades rock sample is carried by just 10% of the number of macropores, which is an indication of the high structural heterogeneity of the rock; in the more homogeneous Ketton rock, 50% of the flow is carried by 10% of the macropores. By analysing the 3D MRI velocity map, it was found that approximately one-third of the total flow rate through the Estaillades rock is carried by microporosity—a porosity that is not captured at the spatial resolution of the µCT image.

Impact of turbulence-induced asymmetric propagators on the accuracy of phase-contrast velocimetry

Phase-contrast magnetic resonance velocimetry (PC-MRI) has been widely used to investigate flow properties in numerous systems. In a horizontal cylindrical pipe (3 mm diameter), we investigated the accuracy of PC-MRI as the flow transitioned from laminar to turbulent flow (Reynolds number 352-2708). We focus primarily on velocimetry errors introduced by skewed intra-voxel displacement distributions, a consequence of PC-MRI theory assuming symmetric distributions. We demonstrated how rapid fluctuations in the velocity field, can produce broad asymmetric intravoxel displacement distributions near the wall. Depending on the shape of the distribution, this resulted in PC-MRI measurements under-estimating (positive skewness) or over-estimating (negative skewness) the true mean intravoxel velocity, which could have particular importance to clinical wall shear stress measurements. The magnitude of these velocity errors was shown to increase with the variance and decrease with the kurtosis of the intravoxel displacement distribution. These experimental results confirm our previous theoretical analysis, which gives a relationship for PC-MRI velocimetry errors, as a function of the higher moments of the intravoxel displacement distribution (skewness, variance, and kurtosis) and the experimental parameters q and Δ. This suggests that PC-MRI errors in such unsteady/turbulent flow conditions can potentially be reduced by employing lower q values or shorter observation times Δ.

Limits to flow detection in phase contrast MRI

Pulsed gradient spin echo (PGSE) complex signal behavior becomes dominated by attenuation rather than oscillation when displacements due to flow are similar or less than diffusive displacements. In this "slow-flow" regime, the optimal displacement encoding parameter q for phase contrast velocimetry depends on the diffusive length scale q s l o w = 1 / l D = 1 / 2 D Δ rather than the velocity encoding parameter v _enc = π/(qΔ). The minimum detectable mean velocity using the difference between the phase at +q _slow and -q _slow is 〈 v m i n 〉 = 1 / SNR D / Δ . These theories are then validated and applied to MRI by performing PGSE echo planar imaging experiments on water flowing through a column with a bulk region and a beadpack region at controlled flow rates. Velocities as slow as 6 μm/s are detected with velocimetry. Theories, MRI experimental protocols, and validation on a controlled phantom help to bridge the gap between porous media NMR and pre-clinical phase contrast and diffusion MRI.

Accuracy of phase-contrast velocimetry in systems with skewed intravoxel velocity distributions

Phase contrast velocimetry (PCV) has been widely used to investigate flow properties in numerous systems. Several authors have reported errors in velocity measurements and have speculated on the sources, which have ranged from eddy current effects to acceleration artefacts. An often overlooked assumption in the theory of PCV, which may not be met in complex or unsteady flows, is that the intravoxel displacement distributions (propagators) are symmetric. Here, the effect of the higher moments of the displacement distribution (variance, skewness and kurtosis) on the accuracy of PCV is investigated experimentally and theoretically. Phase and propagator measurements are performed on tailored intravoxel distributions, achieved using a simple phantom combined with a single large voxel. Asymmetric distributions (Skewness ≠ 0) are shown to generate important phase measurement errors that lead to significant velocimetry errors. Simulations of the phase of the spin vector sum, based on experimentally measured propagators, are shown to quantitatively reproduce the relationship between measured phase and experimental parameters. These allow relating the observed velocimetry errors to a discrepancy between the average phase of intravoxel spins considered in PCV theory and the vector phase actually measured by a PFG experiment. A theoretical expression is derived for PCV velocimetry errors as a function of the moments of the displacement distribution. Positively skewed distributions result in an underestimation of the true mean velocity, while negatively skewed distributions result in an overestimation. The magnitude of these errors is shown to increase with the variance and decrease with the kurtosis of the intravoxel displacement distribution.

Adapted MR velocimetry of slow liquid flow in porous media

MR velocimetry of liquid flow in opaque porous filters may play an important role in better understanding the mechanisms of deep bed filtration. With this knowledge, the efficiency of separating the suspended solid particles from the vertically flowing liquid can be improved, and thus a wide range of industrial applications such as wastewater treatment and desalination can be optimized. However, MR velocimetry is challenging for such studies due to the low velocities, the severe B₀ inhomogeneity in porous structures, and the demand for high spatial resolution and an appropriate total measurement time during which the particle deposition will change velocities only marginally. In this work, a modified RARE-based MR velocimetry method is proposed to address these issues for velocity mapping on a deep bed filtration cell. A dedicated RF coil with a high filling factor is constructed considering the limited space available for the vertical cell in a horizontal MR magnet. Several means are applied to optimize the phase contrast RARE MRI pulse sequence for accurately measuring the phase contrast in a long echo train, even in the case of a low B₁ homogeneity. Two means are of particular importance. One uses data acquired with zero flow to correct the phase contrast offsets from gradient imperfections, and the other combines the phase contrast from signals of both odd and even echoes. Results obtained on a 7T preclinical MR scanner indicate that the low velocities in the heterogeneous system can be correctly quantified with high spatial resolution and an adequate total measurement time, enabling future studies on flow during the filtration process.