scholarly journals Mineral Filling Pattern in Complex Fracture System of Carbonate Reservoirs: Implications from Geochemical Modeling of Water-Rock Interaction

Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-19
Author(s):  
Jianwei Feng ◽  
Xu Zhang ◽  
Peng Luo ◽  
Xizhe Li ◽  
He Du

Much research has been conducted on physical and numerical modeling that focus on stress state and structural controls on subsurface geofluid flow, yet very few attempts have been made to discover and quantify the mineral precipitation/dissolution kinetics in complex fracture system such as Tarim Basin of China. We conducted a geochemical simulation study using the outcrop fracture networks in Ordovician carbonate rocks in Tabei Outcrop Area of Tarim Basin. Structural analysis, filling analysis within the fracture networks and surrounding rocks were used to constrain the generation and geochemical evolution of the geofluids. Using an advanced reactive transport simulation platform TOUGHREACT, a pertinent thermodynamic system was applied to establish the geological model of the fracture-surrounding rock, where the corresponding calcium carbonate (CaCO3) solution was configured to replace the deep saturated hydrothermal fluids. Different types of mineral parameters were considered with material balance and phase equilibrium calculation to perform numerical simulation of multi-field, e.g., pressure field, temperature field, seepage field and chemical field under formation conditions. The simulation results were consistent with field observations. The major findings of this simulation study include: (1) Along with fluid injection, local dissolution occurred within the fractures and matrix, but with the gradual saturation of calcium ions and the increasing pH value, considerable calcite precipitation occurred. (2) The dissolution/precipitation in different fractures was mainly affected by their structure and physical properties, resulting in changes in fluid flow rate, temperature, pressure and ion concentration over time. (3) In the same group, the degree of mineral filling of small-aperture fractures, low-angle fractures and shallow fractures was significantly higher than other types of fractures. (4) The better the connectivity between reticular fractures and the higher the linear density of fractures, the lower the mineral filling degree. (5) Dissolution phenomenon strengthened within large-aperture conjugated fractures gradually along the flow direction. The proposed methodologies in this study can be applied to model effective fracture filling of other deep reservoirs.

Water ◽  
2019 ◽  
Vol 11 (12) ◽  
pp. 2502 ◽  
Author(s):  
Phuong Thanh Vu ◽  
Chuen-Fa Ni ◽  
Wei-Ci Li ◽  
I-Hsien Lee ◽  
Chi-Ping Lin

Fractures are major flow paths for solute transport in fractured rocks. Conducting numerical simulations of reactive transport in fractured rocks is a challenging task because of complex fracture connections and the associated nonuniform flows and chemical reactions. The study presents a computational workflow that can approximately simulate flow and reactive transport in complex fractured media. The workflow involves a series of computational processes. Specifically, the workflow employs a simple particle tracking (PT) algorithm to track flow paths in complex 3D discrete fracture networks (DFNs). The PHREEQC chemical reaction model is then used to simulate the reactive transport along particle traces. The study illustrates the developed workflow with three numerical examples, including a case with a simple fracture connection and two cases with a complex fracture network system. Results show that the integration processes in the workflow successfully model the tetrachloroethylene (PCE) and trichloroethylene (TCE) degradation and transport along particle traces in complex DFNs. The statistics of concentration along particle traces enables the estimations of uncertainty induced by the fracture structures in DFNs. The types of source contaminants can lead to slight variations of particle traces and influence the long term reactive transport. The concentration uncertainty can propagate from parent to daughter compounds and accumulate along with the transport processes.


2021 ◽  
Author(s):  
Cai Deng ◽  
Hua Shi ◽  
Jiali Zhang ◽  
Xinwei Liao ◽  
Zhiming Chen ◽  
...  

Abstract Fracture geometries and drainage radius are important parameters for developing a reasonable development plan of a single fractured well. In some unconventional gas reservoir, some scholars observed the phenomenon of single well controlled reserves increasing through the material balance curve, and put forward the idea of district supply. In addition, owing to fracture hits, the fracture geometries of fractured wells are sometimes more complex. Thus, those complex factors bring challenges for parameter estimations. In order to study the variation of the drainage radius and complex fracture geometries in the single model, a well testing based model for a finite-conductivity fractured vertical well in radial composite reservoirs with dynamic supply and fracture networks is established. Based on "successive steady state method", the point source function, pressure superposition principle and boundary element method are used to solve the reservoir model, and the methods of discrete fracture and pressure superposition are used to solve the fracture model. By introducing the rate normalized pseudo-pressure and material balance time, the variable fluid flux is equivalent to the constant fluid flux. Combined with the inversion idea of well test, the drainage radius value and fracture geometries are solved by fitting the log-log curves of pressure response, and case studies are performed. The results show that the drainage radius increases with the increase of production time and finally tends to a certain value, and it has a good exponential relationship with time. Also, the fracture geometries of the typical well are multiple-radial fracture networks. Through the study of dynamic drainage radius, the controlled reserves of single well in unconventional gas reservoir can be better determined, and it can also provide theoretical basis for fracture evaluation, productivity prediction and enhanced recovery study of the same type of unconventional gas reservoir.


2002 ◽  
Vol 5 (01) ◽  
pp. 49-59 ◽  
Author(s):  
J.L. Pletcher

Summary Experience with material-balance data sets from the field and from simulation has revealed some procedures that can be used to improve analysis of both oil and gas reservoirs:Failure to account for a weak waterdrive can result in significant material-balance errors.The assertion of previous authors that weak waterdrive exhibits a negative slope on the Cole (gas) and Campbell (oil) plots has been confirmed. A weak waterdrive is much more unambiguous on these plots than on commonly used plots, such as the p/z plot for gas.A modified version of the Cole plot is proposed to account for formation compressibility.The reservoir drive indices are a useful tool for determining the correctness of the material-balance solution because they must sum to unity. The drive indices should never be normalized to sum to unity because this obscures their usefulness and leads to a false sense of security.A modified version of the Roach plot (for gas) is proposed that improves interpretation in some waterdrive situations.Material balance has not been replaced by reservoir simulation; rather, it is complementary to simulation and can provide valuable insights to reservoir performance that cannot be obtained by simulation. Introduction Classical material balance is one of the fundamental tools of reservoir engineering. Many authors have addressed the difficult problem of solving the material balance in the presence of a waterdrive (Refs. 1 through 5 are just a few of the more significant ones). The emphasis in the literature has been on strong and moderate waterdrives. In this paper, examples of weak waterdrives are shown in which the effects on the material balance are significant. All aquifers studied here are of the "pot aquifer" type, which is time-independent. In gas reservoirs, the plot of p/z vs. cumulative gas production, Gp, is a widely accepted method for solving the gas material balance1 under depletion-drive conditions. Extrapolation of the plot to atmospheric pressure provides a reliable estimate of original gas in place (OGIP). If a waterdrive is present, the plot often appears to be linear, but the extrapolation will give an erroneously high value for OGIP. Many authors have addressed this problem (including those in Refs. 2 and 5 through 8), especially in cases of strong or moderate waterdrives. The p/z plot is actually more ambiguous in weak waterdrives than in strong or moderate ones. The Cole plot7,9 has proven to be a valuable diagnostic tool for distinguishing between depletion-drive gas reservoirs and those that are producing under a waterdrive. The analogous plot for oil reservoirs is the Campbell plot.10 The literature has emphasized strong and moderate waterdrives, the signature shapes of which are a positive slope and a hump-shaped curve, respectively, on these plots. Previous authors have recognized that weak waterdrives can produce negative slopes on these two diagnostic plots, but this author is not aware of any example plots in the literature. This paper shows examples, using simulation and actual field data, wherein a negative slope clearly reveals a weak waterdrive. These plots are much more diagnostic than the p/z plot. Once a weak waterdrive has been diagnosed, the appropriate steps can be taken in the material-balance equations to yield more accurate results. The Cole plot assumes that formation compressibility can be neglected, which is frequently the case with gas. However, in those reservoirs in which formation compressibility is significant, a modification to the Cole plot is presented that incorporates formation compressibility and gives more accurate results. The reservoir drive indices have been used to quantify the relative magnitude of the various energy sources active in a reservoir. It is shown here that the drive indices are also a useful diagnostic tool for determining the correctness of a material balance solution because they must sum to unity. If they do not sum to unity, a correct solution has not been obtained. In some commercial material-balance software, the drive indices are automatically normalized to sum to unity, which not only obscures their usefulness but also leads to the false impression of having achieved a correct solution. The Roach plot has been presented11 as a tool for solving the gas material balance when formation compressibility is unknown, with or without the presence of waterdrive. This paper shows that for waterdrives that fit the small pot aquifer model, incorporating cumulative water production into the x-axis plotting term improves the linearity of the Roach plot and gives more accurate values for OGIP. Finally, it is argued that even in those reservoirs for which a simulation study is performed, classical material-balance evaluation should be performed on a stand-alone basis. Simulation should not be viewed as a replacement for material balance because the latter can yield valuable insights that can be obscured during simulation. Performing a separate material balance study usually will improve overall reservoir understanding and enhance any subsequent simulation study. Material balance should be viewed as a complement to simulation, not as a competing approach. In this paper, formation compressibility, cf, is assumed to be constant and unchanging over the reservoir life under investigation. References are given for recommended methods to be used in those cases in which cf is variable.


Author(s):  
Hannes Hofmann ◽  
Tayfun Babadagli ◽  
Günter Zimmermann

The creation of large complex fracture networks by hydraulic fracturing is imperative for enhanced oil recovery from tight sand or shale reservoirs, tight gas extraction, and Hot-Dry-Rock (HDR) geothermal systems to improve the contact area to the rock matrix. Although conventional fracturing treatments may result in bi-wing fractures, there is evidence by microseismic mapping that fracture networks can develop in many unconventional reservoirs, especially when natural fracture systems are present and the differences between the principle stresses are low. However, not much insight is gained about fracture development as well as fluid and proppant transport in naturally fractured tight formations. In order to clarify the relationship between rock and treatment parameters, and resulting fracture properties, numerical simulations were performed using a commercial Discrete Fracture Network (DFN) simulator. A comprehensive sensitivity analysis is presented to identify typical fracture network patterns resulting from massive water fracturing treatments in different geological conditions. It is shown how the treatment parameters influence the fracture development and what type of fracture patterns may result from different treatment designs. The focus of this study is on complex fracture network development in different natural fracture systems. Additionally, the applicability of the DFN simulator for modeling shale gas stimulation and HDR stimulation is critically discussed. The approach stated above gives an insight into the relationships between rock properties (specifically matrix properties and characteristics of natural fracture systems) and the properties of developed fracture networks. Various simulated scenarios show typical conditions under which different complex fracture patterns can develop and prescribe efficient treatment designs to generate these fracture systems. Hydraulic stimulation is essential for the production of oil, gas, or heat from ultratight formations like shales and basement rocks (mainly granite). If natural fracture systems are present, the fracturing process becomes more complex to simulate. Our simulation results reveal valuable information about main parameters influencing fracture network properties, major factors leading to complex fracture network development, and differences between HDR and shale gas/oil shale stimulations.


Lithosphere ◽  
2021 ◽  
Vol 2021 (Special 1) ◽  
Author(s):  
Xin Cai ◽  
Wei Liu

Abstract Hydraulic fracturing experiments with low-viscosity fluids, such as supercritical CO2, demonstrate the formation of complex fracture networks spread throughout the rocks. To study the influence of viscosity of the fracturing fluids on hydraulic fracture propagation, a hydromechanical-coupled cohesive zone model is proposed for the simulation of mechanical response of rock grains boundary separation. This simulation methodology considers the synergistic effects of unsteady flow in fracture and rock grain deformation induced by hydraulic pressure. The simulation results indicate a tendency of complex fracture propagation with more branches as the viscosity of fracturing fluids decrease, which is in accord with experimental results. The low-viscosity fluid can flow into the microfractures with extremely small aperture and create more shear failed fracture. This study confirms the possibility of effective well stimulations by hydraulic fracturing with low-viscosity fluids.


Sign in / Sign up

Export Citation Format

Share Document