A New Method for Predicting Acid Penetration Distance

1975 ◽  
Vol 15 (04) ◽  
pp. 277-286 ◽  
Author(s):  
L.D. Roberts ◽  
J.A. GUIN

Abstract A new method for calculating acid penetration distance in fractures bas been developed and tested experimentally. The method combines spending-time data from rotating-disk reaction pots with mass-transfer data obtained from laboratory fractures, thus allowing for both the effects of surface reaction kinetics and actual mixing patterns in the fracture. It is shown that the new method successfully predicts the acid spending obtained in laboratory fractures in both turbulent and laminar flow, using a reaction-rate constant obtained with a rotating-disk apparatus, This appears to be the first method that is easily applicable to small core samples and it allows properly for acid mixing in the fracture. Introduction Recently, there has been considerable interest in and research toward developing a more accurate method for calculating the acid penetration distance in a reservoir fracture. The acid penetration distance, defined as the distance the acid will travel before spending to some predetermined degree, is essential for estimating the production improvement obtainable by fracture acidizing. The first and probably most widely used method for calculating the penetration distance was based on the static-reaction test, in which a small core sample and a known quantity of acid were allowed to react for a given time in a small pot. By equating the spending time of acid in the pot to the residence time of acid flowing down the fracture (t = L/vi), a penetration distance was calculated. it has become penetration distance was calculated. it has become apparent that, because of the extremely fast surface reaction occurring in many acid-rock systems, the over-all acid spending rate to a large degree depends on the extent of fluid mixing at the rock surface. Since fluid-mixing patterns in small reaction pots may not be necessarily the same as those occuring in fractures, several experiments have been performed using actual laboratory model fractures. performed using actual laboratory model fractures. Recent investigations of this type have shown that, because of variable fluid properties, the mixing patterns in real fractures are very complex. To patterns in real fractures are very complex. To allow for this mixing and thereby to calculate more accurately the penetration distance in real fractures, design methods based on experiments in laboratory model fractures have been developed. Therefore, there appear to be two basic approaches to calculating the acid penetration distance, one using data obtained in small reaction pots, and the other using data gathered from model laboratory fractures. Both methods have some advantages. The former is quick, simple to operate, and applicable to the small core samples usually available for tests, while the latter method is more costly, more time consuming, requires special equipment, and is not applicable for use with small core samples. However, for reasons noted above, the latter method is probably more representative of mixing in actual reservoir fractures. In this paper we present a new method for calculating acid penetration distance that combines the advantages of both the above methods without incurring the disadvantages. The new method combines data from both reaction-pot experiments and laboratory-model fracture tests in a manner such that both the reaction rate of the actual rock (obtained conveniently from a small core sample) and the mixing occurring in an actual fracture are allowed for. Reaction-rate constants are obtained using a small batch reaction pot containing a rotating-disk core sample. These rate constants are then used with mass-transfer coefficients obtained from laboratory fractures to predict the acid penetration distance. penetration distance.The combined mass-transfer coefficient/ rate-constant method proposed here has several advantages over existing methods for predicting penetration distance. Since general correlations penetration distance. Since general correlations can be developed for mass-transfer coefficients (in fact, many applicable correlations already exist, most notably in the related field of heat transfer), A is not necessary, nor is it usually possible, to perform experiments in laboratory fractures for each perform experiments in laboratory fractures for each new field core sample obtained. SPEJ P. 277

SPE Journal ◽  
2014 ◽  
Vol 19 (06) ◽  
pp. 1192-1202 ◽  
Author(s):  
Ahmed I. Rabie ◽  
Daniel C. Shedd ◽  
Hisham A. Nasr-El-Din

Summary Lactic acid has been examined in various laboratories and applied in the oil field for acid fracturing and drilling-fluid-filter-cake removal, and as an iron-control agent during acid treatments. However, the reaction of lactic acid with calcite has not been addressed before. Determination of the reaction rate and the acid-diffusion properties is a critical step for successful treatments in matrix acidizing and acid fracturing. Therefore, the objective of this work is to conduct a detailed study on the reaction of lactic acid with calcite. Mass transfer and reaction kinetics are reported for the lactic acid/calcite system by use of the rotating-disk apparatus. Disk samples were cut from Indiana limestone or Silurian dolomite and were used in the reaction-rate experiments. The effect of lactic acid concentration (1, 5, and 10 wt%), temperature (80–250°F), disk rotational speed (100–1,800 rev/min), and different inorganic salts on the reaction rate was investigated. The diffusion coefficient of 5 wt% lactic acid was determined at low disk rotational speeds and reported at 80, 200, and 250°F. A model that accounts for the effect of the kinetics of the surface reactions and the transport of reactants and products was developed. The activation energy and the rate constant at 80, 150, and 250°F for the reaction of lactic acid with Indiana limestone were reported. Reaction experiments of lactic acid with dolomite at 150°F over disk rotational speeds of 100–1,800 rev/min, and at 1,500 rev/min over a temperature range of 80–250°F, were conducted and the results were compared with those obtained for the calcite reaction. At 80°F, the reaction of lactic acid with calcite was controlled by mass transfer at low disk rotational speeds (up to 500 rev/min) and was surface reaction limited at higher speeds. At higher temperatures (150, 200, and 250°F), both mass transfer and surface reaction influence the overall calcite dissolution. The kinetics of the surface reaction were influenced by both forward and backward reactions. At 80°F, the surface reaction contributes to 28% of the overall resistance. This dependence becomes much less (13 and 10%) at higher temperatures (150 and 250°F, respectively). The reaction of lactic acid with dolomite at 150°F was mainly controlled by mass transfer up to 1,000 rev/min and by the kinetics of the surface reaction after 1,000 rev/min. At 80 and 150°F, the rate of reaction of lactic acid with calcite was an order of magnitude higher than that with dolomite. At temperatures of 200 and 250°F, the rate of reaction of lactic acid with calcite is twice the rate of reaction with dolomite. The presence of Ca2+, Mg2+, and SO42− ions reduced the reaction rate, which is most likely because of the reduction in the concentration gradient of the products. The reduction in the concentration gradient will cause a reduction in the rate of diffusion of the generated calcium away from the surface, and hence a lower rate of dissolution.


SPE Journal ◽  
2011 ◽  
Vol 16 (04) ◽  
pp. 981-992 ◽  
Author(s):  
Ahmed I. Rabie ◽  
Ahmed M. Gomaa ◽  
Hisham A. Nasr-El-Din

Summary In-situ-gelled acids have been used extensively in matrix acidizing and acid fracturing for acid diversion and reducing the leakoff rate, respectively. A few studies investigated the rate of dissolution of calcite in polymer-based acids, yet none has addressed in detail the in-situ-gelled acids. Therefore, the aim of this work is to examine the mass transfer and the kinetics of the reaction of 5 wt% HCl in-situ-gelled acids with calcite and determine the effect of Fe crosslinker on the rate of calcite dissolution. The rate of reaction of 5 wt% HCl in-situ-gelled acid was measured using the rotating-disk apparatus. Rock samples of 1.5in. diameter and 1-in. length were used. The effect of temperature (100-250°F) and disk-rotational speed (100-1,800 rev/min) was investigated using Pink Desert limestone rock samples. Calcium concentration was measured in the collected samples and was used to determine the acid-reaction rate. Experimental results showed that the rate of calcite dissolution at 150°F was controlled mainly by the rate of mass transfer of the acid to the surface up to a disk rotational speed of 1,000 rev/min and by the rate of the surface reaction above this value. On the basis of the results obtained, the diffusion coefficient of 5 wt% HCl in in-situ-gelled acid at 150°F; the activation energy; and the reaction rate constant at 150, 200, and 250°F were determined for the first time. A power-law kinetic model was used to determine the kinetics parameters. The presence of Fe3+ crosslinker had a significant effect on the rate of dissolution in comparison with reactions with gelled acid (no crosslinker) at the same condition. The reaction rate decreased by a factor of 2.2 and by a factor of 1.4 when the reaction was conducted at 100 and 1,500 rev/min, respectively. A gel layer, formed on the surface, acted as a barrier between the acid and the rock, which reduced the rate of calcite dissolution.


1978 ◽  
Vol 33 (11) ◽  
pp. 1463-1470 ◽  
Author(s):  
P. Mishra ◽  
Prakash Chandra Singh

2015 ◽  
Vol 782 ◽  
pp. 260-299 ◽  
Author(s):  
Preyas N. Shah ◽  
Eric S. G. Shaqfeh

Surfaces that include heterogeneous mass transfer at the microscale are ubiquitous in nature and engineering. Many such media are modelled via an effective surface reaction rate or mass transfer coefficient employing the conventional ansatz of kinetically limited transport at the microscale. However, this assumption is not always valid, particularly when there is strong flow. We are interested in modelling reactive and/or porous surfaces that occur in systems where the effective Damköhler number at the microscale can be $O(1)$ and the local Péclet number may be large. In order to expand the range of the effective mass transfer surface coefficient, we study transport from a uniform bath of species in an unbounded shear flow over a flat surface. This surface has a heterogeneous distribution of first-order surface-reactive circular patches (or pores). To understand the physics at the length scale of the patch size, we first analyse the flux to a single reactive patch. We use both analytic and boundary element simulations for this purpose. The shear flow induces a 3-D concentration wake structure downstream of the patch. When two patches are aligned in the shear direction, the wakes interact to reduce the per patch flux compared with the single-patch case. Having determined the length scale of the interaction between two patches, we study the transport to a periodic and disordered distribution of patches again using analytic and boundary integral techniques. We obtain, up to non-dilute patch area fraction, an effective boundary condition for the transport to the patches that depends on the local mass transfer coefficient (or reaction rate) and shear rate. We demonstrate that this boundary condition replaces the details of the heterogeneous surfaces at a wall-normal effective slip distance also determined for non-dilute patch area fractions. The slip distance again depends on the shear rate, and weakly on the reaction rate, and scales with the patch size. These effective boundary conditions can be used directly in large-scale physics simulations as long as the local shear rate, reaction rate and patch area fraction are known.


SPE Journal ◽  
2016 ◽  
Vol 21 (01) ◽  
pp. 280-292 ◽  
Author(s):  
John Lyons ◽  
Hadi Nasrabadi ◽  
Hisham A. Nasr-El-Din

Summary Fracture acidizing is a well-stimulation technique used to improve the productivity of low-permeability reservoirs and to bypass deep formation damage. The reaction of injected acid with the rock matrix forms etched channels through which oil and gas can then flow upon production. The properties of these etched channels depend on the acid-injection rate, temperature, reaction chemistry, mass-transport properties, and formation mineralogy. As the acid enters the formation, it increases in temperature by heat exchange with the formation and the heat generated by acid reaction with the rock. Thus, the reaction rate, viscosity, and mass transfer of acid inside the fracture also increase. In this study, a new thermal-fracture-acidizing model is presented that uses the lattice Boltzmann method to simulate reactive transport. This method incorporates both accurate hydrodynamics and reaction kinetics at the solid/liquid interface. The temperature update is performed by use of a finite-difference technique. Furthermore, heterogeneity in rock properties (e.g., porosity, permeability, and reaction rate) is included. The result is a model that can accurately simulate realistic fracture geometries and rock properties at the pore scale and that can predict the geometry of the fracture after acidizing. Three thermal-fracture-acidizing simulations are presented here, involving injection of 15 and 28 wt% of hydrochloric acid into a calcite fracture. The results clearly show an increase in the overall fracture dissolution because of the addition of temperature effects (increasing the acid-reaction and mass-transfer rates). It has also been found that by introducing mineral heterogeneity, preferential dissolution leads to the creation of uneven etching across the fracture surfaces, indicating channel formation.


Sign in / Sign up

Export Citation Format

Share Document