1
|
In Situ Chemically-Selective Monitoring of Multiphase Displacement Processes in a Carbonate Rock Using 3D Magnetic Resonance Imaging. Transp Porous Media 2018; 121:15-35. [PMID: 31983793 PMCID: PMC6954023 DOI: 10.1007/s11242-017-0945-6] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2017] [Accepted: 09/23/2017] [Indexed: 11/24/2022]
Abstract
Accurate monitoring of multiphase displacement processes is essential for the development, validation and benchmarking of numerical models used for reservoir simulation and for asset characterization. Here we demonstrate the first application of a chemically-selective 3D magnetic resonance imaging (MRI) technique which provides high-temporal resolution, quantitative, spatially resolved information of oil and water saturations during a dynamic imbibition core flood experiment in an Estaillades carbonate rock. Firstly, the relative saturations of dodecane (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$S_{\mathrm{o}})$$\end{document}So) and water (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$S_{\mathrm{w}})$$\end{document}Sw), as determined from the MRI measurements, have been benchmarked against those obtained from nuclear magnetic resonance (NMR) spectroscopy and volumetric analysis of the core flood effluent. Excellent agreement between both the NMR and MRI determinations of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$S_{\mathrm{o}}$$\end{document}So and \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$S_{\mathrm{w}}$$\end{document}Sw was obtained. These values were in agreement to 4 and 9% of the values determined by volumetric analysis, with absolute errors in the measurement of saturation determined by NMR and MRI being 0.04 or less over the range of relative saturations investigated. The chemically-selective 3D MRI method was subsequently applied to monitor the displacement of dodecane in the core plug sample by water under continuous flow conditions at an interstitial velocity of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$1.27\times 10^{-6}\,\hbox {m}\,\hbox {s}^{-1}$$\end{document}1.27×10-6ms-1 (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$0.4\,\hbox {ft}\,\hbox {day}^{-1})$$\end{document}0.4ftday-1). During the core flood, independent images of water and oil distributions within the rock core plug at a spatial resolution of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$0.31\,\hbox {mm}\times 0.39\,\hbox {mm} \times 0.39\,\hbox {mm}$$\end{document}0.31mm×0.39mm×0.39mm were acquired on a timescale of 16 min per image. Using this technique the spatial and temporal dynamics of the displacement process have been monitored. This MRI technique will provide insights to structure–transport relationships associated with multiphase displacement processes in complex porous materials, such as those encountered in petrophysics research.
Collapse
|
2
|
Xiao D, Balcom BJ. BLIPPED (BLIpped Pure Phase EncoDing) high resolution MRI with low amplitude gradients. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2017; 285:61-67. [PMID: 29112892 DOI: 10.1016/j.jmr.2017.10.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/18/2017] [Revised: 10/30/2017] [Accepted: 10/30/2017] [Indexed: 06/07/2023]
Abstract
MRI image resolution is proportional to the maximum k-space value, i.e. the temporal integral of the magnetic field gradient. High resolution imaging usually requires high gradient amplitudes and/or long spatial encoding times. Special gradient hardware is often required for high amplitudes and fast switching. We propose a high resolution imaging sequence that employs low amplitude gradients. This method was inspired by the previously proposed PEPI (π Echo Planar Imaging) sequence, which replaced EPI gradient reversals with multiple RF refocusing pulses. It has been shown that when the refocusing RF pulse is of high quality, i.e. sufficiently close to 180°, the magnetization phase introduced by the spatial encoding magnetic field gradient can be preserved and transferred to the following echo signal without phase rewinding. This phase encoding scheme requires blipped gradients that are identical for each echo, with low and constant amplitude, providing opportunities for high resolution imaging. We now extend the sequence to 3D pure phase encoding with low amplitude gradients. The method is compared with the Hybrid-SESPI (Spin Echo Single Point Imaging) technique to demonstrate the advantages in terms of low gradient duty cycle, compensation of concomitant magnetic field effects and minimal echo spacing, which lead to superior image quality and high resolution. The 3D imaging method was then applied with a parallel plate resonator RF probe, achieving a nominal spatial resolution of 17 μm in one dimension in the 3D image, requiring a maximum gradient amplitude of only 5.8 Gauss/cm.
Collapse
Affiliation(s)
- Dan Xiao
- Department of Physics, University of Windsor, Canada; MRI Research Center, Department of Physics, University of New Brunswick, Canada.
| | - Bruce J Balcom
- MRI Research Center, Department of Physics, University of New Brunswick, Canada.
| |
Collapse
|
3
|
Shakerian M, Marica F, Afrough A, Goora FG, Li M, Vashaee S, Balcom BJ. A high-pressure metallic core holder for magnetic resonance based on Hastelloy-C. THE REVIEW OF SCIENTIFIC INSTRUMENTS 2017; 88:123703. [PMID: 29289224 DOI: 10.1063/1.5013031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
A metallic core holder, fabricated from non-magnetic Hastelloy-C276, has been designed for Magnetic Resonance (MR) and Magnetic Resonance Imaging (MRI) of core plug samples at high pressures and temperatures. Core plug samples, 1.5″ in diameter and 2″ in length, can be tested in the core holder at elevated pressures and temperatures, up to 5000 psi and 80 °C. These are conditions commonly found in petroleum reservoirs. A radio frequency probe, which excites and detects magnetic resonance signals, was placed inside the metal vessel. Proximity to the sample improves the signal to noise ratio of the resulting measurements. The metallic core holder is positioned between the poles of a 0.2 T permanent magnet and subjected to rapidly switched magnetic field gradients as part of the imaging process. This switching induces eddy currents on the conductive core holder, which degrades the magnetic field gradient waveform in the sample space. The low electrical-conductivity of Hastelloy-C276 minimizes the duration and the magnitude of such eddy currents. A recently developed pre-equalization technique was employed to ensure that magnetic field gradient pulses, required for MRI, are near ideal in the sample space. A representative core flooding experiment was undertaken in conjunction with MR/MRI measurements.
Collapse
Affiliation(s)
- M Shakerian
- Department of Physics, UNB MRI Centre, University of New Brunswick, Fredericton, New Brunswick E3B 5A3, Canada
| | - F Marica
- Department of Physics, UNB MRI Centre, University of New Brunswick, Fredericton, New Brunswick E3B 5A3, Canada
| | - A Afrough
- Department of Physics, UNB MRI Centre, University of New Brunswick, Fredericton, New Brunswick E3B 5A3, Canada
| | - F G Goora
- Department of Physics, UNB MRI Centre, University of New Brunswick, Fredericton, New Brunswick E3B 5A3, Canada
| | - M Li
- Department of Physics, UNB MRI Centre, University of New Brunswick, Fredericton, New Brunswick E3B 5A3, Canada
| | - S Vashaee
- Department of Physics, UNB MRI Centre, University of New Brunswick, Fredericton, New Brunswick E3B 5A3, Canada
| | - B J Balcom
- Department of Physics, UNB MRI Centre, University of New Brunswick, Fredericton, New Brunswick E3B 5A3, Canada
| |
Collapse
|
4
|
Xiao D, Balcom BJ. T 2 selective π Echo-Planar Imaging for porous media MRI. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2017; 277:52-58. [PMID: 28232286 DOI: 10.1016/j.jmr.2017.02.004] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/16/2016] [Revised: 02/04/2017] [Accepted: 02/06/2017] [Indexed: 06/06/2023]
Abstract
The π Echo Planar Imaging (PEPI) method has recently been modified to permit proton density imaging of fluids in porous media with moderate T2 and short T2∗ signal components. In many applications, it is desirable to discriminate multiple T2 components within each image voxel. T2 selective imaging is explored in this paper through adiabatic inversion as a magnetization preparation with PEPI readout. When prior information of the sample relaxation times is known, responses of different species to broadband adiabatic inversion pulses can be predicted by Bloch equation simulation. Different relaxation components can be acquired by combining the images with and without inversion preparation pulses. T2 weighting can be easily introduced in the PEPI sequence by shifting the spatial encoding gradients based on its spin echo nature. T2 decay curves can be extracted for each image voxel from a series of T2 weighted images and spatially resolved T2 distributions can be generated. This method is reliable but slow. The two methods were implemented to image porous media samples with PEPI the common basis of spatial resolution. The results of both methods agree remarkably well.
Collapse
Affiliation(s)
- Dan Xiao
- MRI Research Center, Department of Physics, University of New Brunswick, 8 Bailey Drive, Fredericton NB E3B 5A3, Canada.
| | - Bruce J Balcom
- MRI Research Center, Department of Physics, University of New Brunswick, 8 Bailey Drive, Fredericton NB E3B 5A3, Canada.
| |
Collapse
|
5
|
Ramskill NP, Bush I, Sederman AJ, Mantle MD, Benning M, Anger BC, Appel M, Gladden LF. Fast imaging of laboratory core floods using 3D compressed sensing RARE MRI. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2016; 270:187-197. [PMID: 27500742 DOI: 10.1016/j.jmr.2016.07.017] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/16/2016] [Revised: 07/26/2016] [Accepted: 07/27/2016] [Indexed: 06/06/2023]
Abstract
Three-dimensional (3D) imaging of the fluid distributions within the rock is essential to enable the unambiguous interpretation of core flooding data. Magnetic resonance imaging (MRI) has been widely used to image fluid saturation in rock cores; however, conventional acquisition strategies are typically too slow to capture the dynamic nature of the displacement processes that are of interest. Using Compressed Sensing (CS), it is possible to reconstruct a near-perfect image from significantly fewer measurements than was previously thought necessary, and this can result in a significant reduction in the image acquisition times. In the present study, a method using the Rapid Acquisition with Relaxation Enhancement (RARE) pulse sequence with CS to provide 3D images of the fluid saturation in rock core samples during laboratory core floods is demonstrated. An objective method using image quality metrics for the determination of the most suitable regularisation functional to be used in the CS reconstructions is reported. It is shown that for the present application, Total Variation outperforms the Haar and Daubechies3 wavelet families in terms of the agreement of their respective CS reconstructions with a fully-sampled reference image. Using the CS-RARE approach, 3D images of the fluid saturation in the rock core have been acquired in 16min. The CS-RARE technique has been applied to image the residual water saturation in the rock during a water-water displacement core flood. With a flow rate corresponding to an interstitial velocity of vi=1.89±0.03ftday(-1), 0.1 pore volumes were injected over the course of each image acquisition, a four-fold reduction when compared to a fully-sampled RARE acquisition. Finally, the 3D CS-RARE technique has been used to image the drainage of dodecane into the water-saturated rock in which the dynamics of the coalescence of discrete clusters of the non-wetting phase are clearly observed. The enhancement in the temporal resolution that has been achieved using the CS-RARE approach enables dynamic transport processes pertinent to laboratory core floods to be investigated in 3D on a time-scale and with a spatial resolution that, until now, has not been possible.
Collapse
Affiliation(s)
- N P Ramskill
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Pembroke Street, Cambridge CB2 3RA, UK.
| | - I Bush
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Pembroke Street, Cambridge CB2 3RA, UK
| | - A J Sederman
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Pembroke Street, Cambridge CB2 3RA, UK
| | - M D Mantle
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Pembroke Street, Cambridge CB2 3RA, UK
| | - M Benning
- Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, UK
| | - B C Anger
- Shell Technology Centre, 3333 Highway 6 S, Houston, TX, USA
| | - M Appel
- Shell Technology Centre, 3333 Highway 6 S, Houston, TX, USA
| | - L F Gladden
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Pembroke Street, Cambridge CB2 3RA, UK
| |
Collapse
|
6
|
Li M, Xiao D, Romero-Zerón L, Marica F, MacMillan B, Balcom BJ. Mapping three-dimensional oil distribution with π-EPI MRI measurements at low magnetic field. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2016; 269:13-23. [PMID: 27208417 DOI: 10.1016/j.jmr.2016.05.008] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/13/2016] [Revised: 05/10/2016] [Accepted: 05/11/2016] [Indexed: 06/05/2023]
Abstract
Magnetic resonance imaging (MRI) is a robust tool to image oil saturation distribution in rock cores during oil displacement processes. However, a lengthy measurement time for 3D measurements at low magnetic field can hinder monitoring the displacement. 1D and 2D MRI measurements are instead often undertaken to monitor the oil displacement since they are faster. However, 1D and 2D images may not completely reflect the oil distribution in heterogeneous rock cores. In this work, a high-speed 3D MRI technique, π Echo Planar Imaging (π-EPI), was employed at 0.2T to monitor oil displacement. Centric scan interleaved sampling with view sharing in k-t space was employed to improve the temporal resolution of the π-EPI measurements. A D2O brine was employed to distinguish the hydrocarbon and water phases. A relatively homogenous glass bead pack and a heterogeneous Spynie core plug were employed to show different oil displacement behaviors. High quality 3D images were acquired with π-EPI MRI measurements. Fluid quantification with π-EPI compared favorably with FID, CPMG, 1D-DHK-SPRITE, 3D Fast Spin Echo (FSE) and 3D Conical SPRITE measurements. π-EPI greatly reduced the gradient duty cycle and improved sensitivity, compared to FSE and Conical SPRITE measurements, enabling dynamic monitoring of oil displacement processes. For core plug samples with sufficiently long lived T2, T2(∗), π-EPI is an ideal method for rapid 3D saturation imaging.
Collapse
Affiliation(s)
- Ming Li
- MRI Centre, Department of Physics, P.O. Box 4400, University of New Brunswick, Fredericton, NB E3B 5A3, Canada; Department of Chemical Engineering, P.O. Box 4400, University of New Brunswick, Fredericton, NB E3B 5A3, Canada
| | - Dan Xiao
- MRI Centre, Department of Physics, P.O. Box 4400, University of New Brunswick, Fredericton, NB E3B 5A3, Canada
| | - Laura Romero-Zerón
- Department of Chemical Engineering, P.O. Box 4400, University of New Brunswick, Fredericton, NB E3B 5A3, Canada
| | - Florea Marica
- MRI Centre, Department of Physics, P.O. Box 4400, University of New Brunswick, Fredericton, NB E3B 5A3, Canada
| | - Bryce MacMillan
- MRI Centre, Department of Physics, P.O. Box 4400, University of New Brunswick, Fredericton, NB E3B 5A3, Canada
| | - Bruce J Balcom
- MRI Centre, Department of Physics, P.O. Box 4400, University of New Brunswick, Fredericton, NB E3B 5A3, Canada.
| |
Collapse
|