No Arabic abstract
This article provides an overview of the current state of digital rock technology, with emphasis on industrial applications. We show how imaging and image analysis can be applied for rock typing and modeling of end-point saturations. Different methods to obtain a digital model of the pore space from pore scale images are presented, and the strengths and weaknesses of the different methods are discussed. We also show how imaging bridges the different subjects of geology, petrophysics and reservoir simulations, by being a common denominator for results in all these subjects. Network modeling is compared to direct simulations on grid models, and their respective strengths are discussed. Finally we present an example of digital rock technology applied to a sandstone oil reservoir. Results from digital rock modeling are compared to results from traditional laboratory experiments. We highlight the mutual benefits from conducting both traditional experiments and digital rock modeling.
Permeability is the key parameter for quantifying fluid flow in porous rocks. Knowledge of the spatial distribution of the connected pore space allows, in principle, to predict the permeability of a rock sample. However, limitations in feature resolution and approximations at microscopic scales have so far precluded systematic upscaling of permeability predictions. Here, we report fluid flow simulations in capillary network representations designed to overcome such limitations. Performed with an unprecedented level of accuracy in geometric approximation at microscale, the pore scale flow simulations predict experimental permeabilities measured at lab scale in the same rock sample without the need for calibration or correction. By applying the method to a broader class of representative geological samples, with permeability values covering two orders of magnitude, we obtain scaling relationships that reveal how mesoscale permeability emerges from microscopic capillary diameter and fluid velocity distributions.
Phase-field modeling -- a continuous approach to discontinuities -- is gaining popularity for simulating rock fractures due to its ability to handle complex, discontinuous geometry without an explicit surface tracking algorithm. None of the existing phase-field models, however, incorporates the impact of surface roughness on the mechanical response of fractures -- such as elastic deformability and shear-induced dilation -- despite the importance of this behavior for subsurface systems. To fill this gap, here we introduce the first framework for phase-field modeling of rough rock fractures. The framework transforms a displacement-jump-based discrete constitutive model for discontinuities into a strain-based continuous model, and then casts it into a phase-field formulation for frictional interfaces. We illustrate the framework by constructing a particular phase-field form employing a rock joint model originally formulated for discrete modeling. The results obtained by the new formulation show excellent agreement with those of a well-established discrete method for a variety of problems ranging from shearing of a single discontinuity to compression of fractured rocks. Consequently, our phase-field framework provides an unprecedented bridge between a discrete constitutive model for rough discontinuities -- common in rock mechanics -- and the continuous finite element method -- standard in computational mechanics -- without any algorithm to explicitly represent discontinuity geometry.
SNOLAB is one of the deepest underground laboratories in the world with an overburden of 2092 m. The SNO+ detector is designed to achieve several fundamental physics goals as a low-background experiment, particularly measuring the Earths geoneutrino flux. Here we evaluate the effect of the 2 km overburden on the predicted crustal geoneutrino signal at SNO+. A refined 3D model of the 50 x 50 km upper crust surrounding the detector and a full calculation of survival probability are used to model the U and Th geoneutrino signal. Comparing this signal with that obtained by placing SNO+ at sea level, we highlight a $1.4^{+1.8}_{-0.9}$ TNU signal difference, corresponding to the ~5% of the total crustal contribution. Finally, the impact of the additional crust extending from sea level up to ~300 m was estimated.
A multi-scale scheme for the invasion percolation of rock fracture networks with heterogeneous fracture aperture fields is proposed. Inside fractures, fluid transport is calculated on the finest scale and found to be localized in channels as a consequence of the aperture field. The channel network is characterized and reduced to a vectorized artificial channel network (ACN). Different realizations of ACNs are used to systematically calculate efficient apertures for fluid transport inside differently sized fractures as well as fracture intersection and entry properties. Typical situations in fracture networks are parameterized by fracture inclination, flow path length along the fracture and intersection lengths in the entrance and outlet zones of fractures. Using these scaling relations obtained from the finer scales, we simulate the invasion process of immiscible fluids into saturated discrete fracture networks, which were studied in previous works.
There have been multiple attempts to demonstrate that quantum annealing and, in particular, quantum annealing on quantum annealing machines, has the potential to outperform current classical optimization algorithms implemented on CMOS technologies. The benchmarking of these devices has been controversial. Initially, random spin-glass problems were used, however, these were quickly shown to be not well suited to detect any quantum speedup. Subsequently, benchmarking shifted to carefully crafted synthetic problems designed to highlight the quantum nature of the hardware while (often) ensuring that classical optimization techniques do not perform well on them. Even worse, to date a true sign of improved scaling with the number of problem variables remains elusive when compared to classical optimization techniques. Here, we analyze the readiness of quantum annealing machines for real-world application problems. These are typically not random and have an underlying structure that is hard to capture in synthetic benchmarks, thus posing unexpected challenges for optimization techniques, both classical and quantum alike. We present a comprehensive computational scaling analysis of fault diagnosis in digital circuits, considering architectures beyond D-wave quantum annealers. We find that the instances generated from real data in multiplier circuits are harder than other representative random spin-glass benchmarks with a comparable number of variables. Although our results show that transverse-field quantum annealing is outperformed by state-of-the-art classical optimization algorithms, these benchmark instances are hard and small in the size of the input, therefore representing the first industrial application ideally suited for testing near-term quantum annealers and other quantum algorithmic strategies for optimization problems.