Lessons for machine learning from the analysis of porosity-permeability transforms for carbonate reservoirs
Not provided.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not provided.
Not provided.
Geothermal formations often contain extensive fracture networks. These fracture networks contribute to the significant loss of drilling fluids during geothermal drilling. Multiple loss circulation materials (LCM) such as fiber, granules, and pills have been proposed to tackle this problem but with only limited success. Recent advances in materials science have led to the development of thermoset shape memory polymers (SMP) to address the lost circulation problem. In this paper, we evaluate a thermoset SMP performance in sealing near wellbore fractures of different sizes in geothermal wells. The SMP performance was assessed using granite disks and cylindrical granite cores having fracture sizes of 1000 μm and 3000 μm. A static filtration test was performed using cedar fiber, CaCO 3 , and SMP. Results showed cedar fiber performed better than the CaCO 3 ., reducing fluid loss by 89% and improving sealing pressure by 200 psi. A novel dynamic testing unit that allows for high-temperature testing under flowing conditions was used in this study. The analysis showed that 3% by weight SMP and fiber blends could bridge and plug the 1000 μm fracture. For a larger fracture of 3000 μm width, there was a need to increase the weight concentration of the SMP to 6% to plug the fracture opening effectively. We showed the influence of key parameters such as the type of LCM, concentration, and particle size distribution in optimizing the performance of drilling fluid loss treatment.
Development and preservation of carbonate porosity and permeability are critical to characterizing reservoirs. However, secondary porosity, such as vugs and fractures, are difficult to identify and require collection of expensive wireline tools or core sampling. Wireline logs and cores have traditionally been used to identify the presence of secondary porosity but fail to quantify the contribution to total reservoir porosity. Additionally, advanced wireline logs and core are not readily available for most wells. Dual energy CT scans were collected on whole core from the A-1 Carbonate and Brown Niagaran formations drawn from six wells in northern Michigan. 3D analysis techniques were applied to identify and isolate secondary porosity features. A series of machine learning and data analytics techniques were applied to the dataset to predict secondary porosity features on basic wireline logs. The predictive model successfully predicted secondary porosity with high confidence.
We propose a novel approach for rate optimization during a waterflood under geologic uncertainty in reservoir properties such as permeability and porosity. The traditional approach typically involves several runs of the forward simulator. This may not scale well when the optimization is to be performed at the full field-level and over multiple geologic realizations. A machine-learning (ML) based approach which is quick and scalable for rate optimization over multiple geologic realizations is proposed instead. The training data for the model is generated by running the forward simulator with randomly assigned well rates using multiple geologic realizations. A reduced order representation of the permeability heterogeneity in each of the realizations is derived using a grid connectivity transformation (GCT). This step involves finding basis functions corresponding to the different modal frequencies of the grid connectivity represented by the grid Laplacian. The projection of the heterogeneous property field along these basis functions gives the basis coefficients that form the reduced order representation. Subsequently, for each training datapoint, streamlines are traced and the minimum time of flight (TOF) representing the tracer breakthrough time at each producer is recorded. The basis coefficients and well rates are fed to a machine learning model as input and the minimum TOF at the producers forms the output of the model. This trained model can then be used along with an optimizer for computing the optimal injection rates to maximize the injection sweep efficiency. This corresponds to minimizing the variance in the minimum TOF within each well group. Different architectures of neural network are tested using 5-fold cross validation to decide the best ML model to compute the streamline time of flight. The trained model is used to perform well rate optimization over multiple realizations of geology by using a risk tolerance penalty. The optimal well rates thus obtained are compared with two cases: a) equal well rates assigned to all injectors and producers and b) well rates obtained by optimizing over a single realization without considering the uncertainty in geology. The optimal well rates are seen to offer better oil recovery and sweep efficiency than both cases.
Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.
This article studies the rheology and annular flow of a smart lost circulation material (LCM) at 160 °C, using a high-temperature flow loop. The smart LCM is a shape memory polymer (SMP) that activates at high temperatures, and its particle size increases to seal fractures in geothermal formations. SMP was mixed with the base fluid in two different concentrations, 1.0 and 3.0 wt%, to study rheology, wellbore hydraulics, activation process, and settling behavior under different testing conditions. The results of this study showed that the SMP could be activated at high temperatures. An increase of 80–100% in the particle size was observed at 160 °C. The mud samples showed a high shear-thinning behavior at the two concentrations with a Power-law flow index (n) ranging between 0.025 and 0.101. No additional frictional pressure losses were observed when SMP was added to the base fluid. SMP particles showed an excellent suspension at 1.0 wt% while, at 3.0 wt%, a bed was formed at a low flow rate and without pipe rotation. Increasing the drill pipe rotational speed or flow rate effectively removed the bed and homogeneously dispersed the SMP particles, ensuring a better sealing efficiency. SMP particle dispersion in inclined wells was better than in horizontal wells. Moreover, the findings of this study help optimize the lost circulation treatment by considering a wide range of operating parameters that can further be extended to different systems and geometries.
Gel treatment is an effective way to attack excessive water production in many mature oilfields around the world. Selective penetration is desired for successful gel treatments. That is, gel materials should easily penetrate the target zones (i.e., channeling features such as superpermeable channels) without entering/damaging the nontarget zones (i.e., reservoir matrices or oil zones). This study revealed that presence of threshold penetration pressure (ΔP th ) was responsible for selective penetration behavior of tested microgels. The concept of ΔP th was utilized to figure out favorable working conditions for effective gel treatments. Microgel dispersions were injected into superpermeable (super-k) sandpacks (mimicking super-k channels in reservoirs, 60–221 darcies), heterogeneous models with super-k channels (79–230 darcies), and sandstone cores (mimicking reservoir matrices, 50–5000 md). The results demonstrated that a minimum differential driving pressure (i.e., threshold penetration pressure, ΔP th ) was required to push microgel particles to penetrate channels or matrices. The critical penetration behavior was closely related to the particle/pore size ratio. Low ΔP th at smaller particle/pore ratios was beneficial to allow easy penetration of gel materials into the channeling zones. On the contrary, high ΔPth at larger particle/pore ratios was desirable to prevent gel materials from massively invading and damaging the matrices. Instead, the gel particles accumulated at the inlet surface, and a gel cake was gradually formed. The cake further prevented the invasion of the gels. The cake could be removed by chemical breakers to resume the injectivity/productivity of the matrices. Correlations were developed to describe the relationship between ΔP th and particle/pore ratio. A distinct transition was identified at the particle/pore ratio of about 3. This research could help identify the favorable conditions to achieve successful gel treatments. In an effective conformance treatment, the particle/pore ratio in the channel should be sufficiently low to allow easy penetration of gel materials into the channel (e.g., particle/pore ratio<2 in this study). Meanwhile, the particle/pore ratio in the matrix should be large enough to support a high ΔP th and thus prevent massive gel invasion into the matrix. This study advances the current pore scale studies (a single particle passing through a single channel) to Darcy-scale characterization.
Production decline from wells producing from shale rocks is typically quite rapid. One possible factor responsible for such decline is loss in fracture conductivity. This is true in the case of the Tuscaloosa Marine shale reservoir. In this work, the stress-dependent fracture conductivity of TMS core samples is investigated using five propped cores under varying stress conditions. In addition, the effects of static rock mechanical properties and rock mineralogy on conductivity decline are investigated. The mineralogy of each core sample was determined by X-ray Diffraction analysis. Fracture conductivity was measured using a Hassler-type pressure core holder with a pressure rating of 15,000 psi (103.42 MPa). Estimated rock mechanical properties showed that shale anisotropy could affect fracture conductivity. The time decline in conductivity revealed two possible decline trends, a) decline from the reduction in induced fracture width, and b) decline from the healing of developed micro-cracks. An exponential decline in fracture conductivity with increasing confinement pressure was observed with an average decline rate constant of 3.15 x 10 -10 mPa -1 . While rock mineralogy did not have any clear relationship with conductivity decline, the effects of rock mechanical properties were only important at 10.34 MPa confinement pressure and above. The observation showed that a high Young's modulus correlate to high conductivity, whereas a high Poisson's ratio correlate to low conductivity. This study provides useful observations for optimizing fracture design to enhance well productivity in the TMS.
Unique characteristics of clay (e.g., high surface area, deficiency in positive charges) make any clay abundant shales a good source of adsorbed gas. However, the presence of clay in porous media can also affect shale rock's producibility. The clay-water interaction can cause significant swelling and fines migration in the formation and impair the overall hydrocarbon recovery from shale. Geomaterial micromodels, developed by functionalizing traditional glass or PDMS surfaces with geomaterials (e.g., calcite, Quartz, clay), can represent the physicochemical properties of the natural porous media and have been used during the last few years to understand and evaluate the solid-fluid physicochemical interactions. Here, this study focuses on developing a clay-coated geomaterial surface and investigates the effect of base fluid's salinity on clay adsorption to the glass surface. Glass capillary tubes and straight channel borosilicate glass micromodels are coated with Illite clay to represent the pore-scale clay chemistry of Caney Shale, a Mississippian unconventional play in Southern Oklahoma, USA. 10 wt% of Illite clay slurries made with brines of four different salinities are used to coat the glass-capillary tubes to understand the effect of salinity in clay adsorption on the glass surface and overall coating quality. To achieve a stable coating and evaluate the impact of heat treatment on coating's stability, straight channel glass flow cells coated with Illite clay are heat-treated at low (25 °C) and high temperature (125 °C), and the results are compared with an untreated coated surface. Flooding tests were carried out with brines of different salinities on the coated surface to evaluate the stability of the coating. The experimental results indicated a strong relationship between the brine's salinity and the adsorption of clay particles on the glass surface. An increase in brine concentration resulted in improved adsorption of clay particles on the glass surface. Experiments involving heat treating the glass surface following the coating demonstrated significant improvement in the stability of the coating.
Optimum solutions for geologic modeling and reservoir simulation in industries such as oil and gas recovery and carbon capture and storage require accurate characterization of reservoir properties, which are often heterogeneous. In this study, high-quality micro-computed tomography (CT) images (1.475-μm/pixel resolution) of a sandstone core acquired from the Bell Creek oil field, USA, were used to provide nondestructive analysis of pore- and core-scale heterogeneity across measurement scales of 94–566 μm. In addition to characterizing the as-received sample, the core sample was flooded with brine to evaluate the capacity of the core sample to receive injected fluids. The micro-CT images were systematically segmented into pore spaces and grains via machine learning (ML) steps including image preprocessing, label creation using a traditional ML method based on limited manual image annotation, and finally U-Net segmentation. The segmented image stacks were reconstructed into digital cubes of various scales of voxel lengths. The 3D porosity values were calculated for all the digital cubes, and the fractal dimensions of the cubes were estimated using a box-counting method. The results showed that smaller cubes had greater heterogeneity and that the porosity values could be accurately estimated by fractal dimension and voxel lengths using ML models. For the core sample with brine flooding, the ratio of pores filled by brine to the total pore space was related to the porosity and could also be accurately estimated by porosity, fractal dimension, and voxel lengths using ML models. In conclusion, the results of this study demonstrate that the concept of fractal dimension can be a useful vector to perform upscaling analysis of sandstone rock heterogeneity from the pore to core scale and that fractal dimensions can be used to estimate porosity values and pore space-filling capacity across those scales.