Pore-Structural Characteristics of Tight Fractured-Vuggy Carbonates and Its Effects on the P- and S-Wave Velocity: A Micro-CT Study on Full-Diameter Cores

Pore structure has been widely observed to affect the seismic wave velocity of rocks. Although taking lab measurements on 1.0-inch core plugs is popular, it is not representative of the fractured-vuggy carbonates because many fractures and vugs are on a scale up to several hundred microns (and greater) and are spatially heterogeneous. To overcome this shortage, we carried out the lab measurements on full-diameter cores (about 6.5–7.5 cm in diameter). The micro-CT (micro computed tomography) scanning technique is used to characterize the pore space of the carbonates and image processing methods are applied to filter the noise and enhance the responses of the fractures so that the constructed pore spaces are reliable. The wave velocities of P- and S-waves are determined then and the effects of the pore structure on the velocity are analyzed. The results show that the proposed image processing method is effective in constructing and quantitatively characterizing the pore space of the full-diameter fractured-vuggy carbonates. The porosity of all the collected tight carbonate samples is less than 4%. Fractures and vugs are well-developed and the spatial distributions of them are heterogeneous causing, even the samples having similar porosity, the pore structure characteristics of the samples being significantly different. The pores and vugs mainly contribute to the porosity of the samples and the fractures contribute to the change in the wave velocities more than pores and vugs.

A Model To Simulate Matrix-Acid Stimulation for Wells in Dolomite Reservoirs with Vugs and Natural Fractures

Summary Designing matrix-acid stimulation treatments in vuggy and naturally fractured carbonate reservoirs is a challenging problem in the petroleum industry. It is often difficult to physically model this process, and current mathematical models do not consider vugs or fractures. There is a significant gap in the literature for models that design and evaluate matrix-acid stimulation in vuggy and naturally fractured carbonate reservoirs. The objective of this work is to develop a new model to simulate matrix acidizing under field conditions in vuggy and naturally fractured carbonates. To obtain accurate and reliable simulation parameters, acidizing coreflood experiments were modeled using a reactive-flow simulator. A 3D radial field-scale model was used to study the flow of acid in the presence of vugs (pore spaces that are significantly larger than grains) and natural fractures (breaks in the reservoir that were formed naturally by tectonic events). The vugs’ size and distribution effects on acid propagation were studied under field conditions. The fracture length, conductivity, and orientation, and the number of fractures in the formation, were studied by the radial model. The results of the numerical simulation were used to construct Gaussian-process (GP)-based surrogate models for predicting acid propagation in vuggy and naturally fractured carbonates. Finally, the acid propagation in vuggy/naturally fractured carbonates was evaluated, as well.The simulation results of vuggy carbonates show that the presence of vugs in carbonates results in faster and deeper acid propagation in the formation when compared with homogeneous reservoirs at injection velocities lower than 8×10–4 m/s. Results also revealed that the size and density of the vugs have a significant impact on acid consumption and the overall performance of the acid treatment. The output of the fracture model illustrates that under field conditions, fracture orientations do not affect the acid-propagation velocity. The acid does not touch all of the fractures around the well. The GP model predictions have an accuracy of approximately 90% for both vuggy and naturally fractured cases. The vuggy/naturally fractured model simulations reveal that fractures are the main reason behind the fast acid propagation in these highly heterogeneous reservoirs.

Modeling Flow in Porous Media With Double Porosity/Permeability: Mathematical Model, Properties, and Analytical Solutions

Geomaterials such as vuggy carbonates are known to exhibit multiple spatial scales. A common manifestation of spatial scales is the presence of (at least) two different scales of pores with different hydromechanical properties. Moreover, these pore-networks are connected through fissures and conduits. Although some models are available in the literature to describe flows in such media, they lack a strong theoretical basis. This paper aims to fill this gap in knowledge by providing the theoretical foundation for the flow of incompressible single-phase fluids in rigid porous media that exhibit double porosity/permeability. We first obtain a mathematical model by combining the theory of interacting continua and the maximization of rate of dissipation (MRD) hypothesis, and thereby provide a firm thermodynamic underpinning. The governing equations of the model are a system of elliptic partial differential equations (PDEs) under a steady-state response and a system of parabolic PDEs under a transient response. We then present, along with mathematical proofs, several important mathematical properties that the solutions to the model satisfy. We also present several canonical problems with analytical solutions which are used to gain insights into the velocity and pressure profiles, and the mass transfer across the two pore-networks. In particular, we highlight how the solutions under the double porosity/permeability differ from the corresponding ones under Darcy equations.