Research Article  Open Access
Jianye Mou, Lei Wang, Shicheng Zhang, Xinfang Ma, Boyang Li, "A Research on the Effect of Heterogeneities on Sandstone Matrix Acidizing Performance", Geofluids, vol. 2019, Article ID 6328909, 14 pages, 2019. https://doi.org/10.1155/2019/6328909
A Research on the Effect of Heterogeneities on Sandstone Matrix Acidizing Performance
Abstract
Matrix acidizing is one of the common methods to enhance production in sandstone reservoirs. Conventional acidizing designs generally neglected the effect of heterogeneities of mineral and flow field distributions both in areal and vertical directions and assumed that the acid front propagates with a pistonlike style. However, sandstone formations inevitably have smallscale heterogeneities of minerals and flow properties that may give rise to acid propagation in a manner much different from what is predicted based on homogeneous assumptions. In this paper, we conduct a research to numerically investigate how the heterogeneities affect acidizing performance under in situ conditions. Firstly, a heterogeneity model is built for mineral and porosity distributions by using the semivariogram model of geological statistics, based on which we generate spatially correlated porosity and mineral distributions. Next, a model of radial acid flooding is developed based on mass balance and the chemical reactions between the acids and minerals occurring during the acidizing process. The model is numerically solved to investigate the permeability response, acid distributions, precipitate distributions, and the effect of the heterogeneities on acidizing. The results show that the heterogeneities both in areal and vertical directions have a significant effect on acidizing. The flow field heterogeneities have a more serious impact than the mineral heterogeneities. In a plane, strong porosity heterogeneity can give rise to acid fingering and even channeling, which make the acid penetration distance longer than the homogeneous cases. The secondary precipitate has a significant effect when fastreacting mineral content is high. Vertically, severalfold permeability contrast creates the acid break through the highperm zone leaving the lowperm zone understimulated. Both flow field and mineral heterogeneities make it possible to create highpermeability channels during the acidizing process and to obtain a longer acid penetration distance.
1. Introduction
Drilling, completion, or some other well operations inevitably damage formations, which may seriously decrease the well productivity in sandstone reservoirs. Matrix acidizing is a common method to remove the formation damage and recover well productivity. In acidizing, the acid flows into the porous media, reacts with the minerals, and increase formation porosity as well as permeability. Many researchers conducted studies on acid flooding in experiments. One of the unignorable risks in acidizing in sandstone is the secondary precipitate. Crowe [1], Almond et al. [2], and Underdown et al. [3] investigated the conditions that produce precipitate. They found that silica gel is produced at a slow rate due to the slow reaction at room temperature. When the temperature is higher than 55°C, the fast reaction will generate a lot of precipitate and damage the formation seriously.
Due to the small size of the core sample and the complexity of parallel core flooding, fieldscale simulation and design necessitates numerical modeling. Sandstone acidizing models include the capillary model, the microscopic model, the kinetic model, the lumped parameter model, the distributed parameter model, the oneacid twomineral model, the twoacid threemineral model, and the generalized geochemistry model [4–10], most of which made the assumption of homogenous formation. The twoacid threemineral model developed by Bryant [4] takes into account the effect of precipitates and H_{2}SiF_{6} and can simulate well the physical phenomenon of acidizing. At present, the most sophisticated model [9] is the generalized geochemistry model. The model accounts for the reaction of multiple kinds of minerals and the acid and can analyze the effect of mineral content, reactant, product, and reaction kinetics. They found that dissolution and precipitation do not occur simultaneously. Li et al. [11] built a linear acid flooding model to simulate linear core flooding at an experimental scale with the consideration of the heterogeneities of minerals and flow fields. Li et al. [12] built a radial acidizing model and considered the effect of the temperature field on the reaction, but they did not consider the effect of the heterogeneities. Leong et al. [13, 14] did a modeling of the acidizing of HBF_{4} in linear cores with COMSOL.
In this paper, we modeled radial acid flooding based on a twoacid, threemineral model and considered the heterogeneities of minerals and flow fields in both planar direction and vertical directions. Based on the model, we did extensive numerical simulations to analyze the effect of heterogeneities on acidizing performance. Also considered are multilayer flooding and secondary precipitates.
2. Mathematical Model
In the twoacid threemineral model, the two acids mean HF and , and the three minerals mean fastreacting mineral (mineral 1), slowreacting mineral (mineral 2), and silica gel (mineral 3). The chemical reactions are as follows:
2.1. Governing Equations
Based on material balance, Darcy’s law, and acidrock reaction kinetics, we establish our governing equations as equations 2–7 with the following assumptions: (1) singlephase flow, (2) Darcy flow, (3) incompressible fluid and rock, and (4) neglecting gravity effect. A 2D model is developed as shown in Figure 1 with the wellbore at the center of the domain. (1)Cylindrical coordinate system is used for the radial flooding(2)Flow equations (3)Hydrofluoric acid concentration distribution equation (4)Fluosilicic acid concentration distribution equation (5)Mineral content equation (6)Porosity variation equation where is the porosity; is the hydrofluoric acid concentration; is fluosilicic acid concentration; is the velocity vector ( is the velocity in the radial direction, and is the velocity in the circumferential direction); is the mineral volume fraction; is the fastreacting mineral, slowreacting mineral, and silica gel; is the mole weight of acid, is the specific surface area of minerals, is the hydrofluoric acid or fluosilicic acid; is the reaction rate constant; is the reaction order; is the density of minerals; and is the gravimetric dissolving power.
2.2. Boundary and Initial Conditions
The initial conditions are as follows: where , , and are initial mineral 1, mineral 2, and mineral 3 distribution, respectively. is initial porosity distribution. is reservoir pressure.
The boundary conditions are as follows: where is pump rate. and are the injection concentrations of acids 1 and 2, respectively. is the pressure of the outer boundary.
3. Porosity and Mineral Distributions
3.1. Porosity Distribution
Change of porosity directly affects the permeability of the sandstone formation, thus affecting the distribution and reaction rate of acid fluid. Therefore, it is very important to generate the porosity distribution of a real formation for the accuracy of the whole model. In terms of geostatistics, porosity distribution shows characteristics of directionality [15–17]; that is, they are spatially correlated instead of completely random. A semivariogram model [8, 18], which is an important tool in geostatistics to describe spatially correlated distributions, is used to generate porosity distribution. We use the geostatistical software GSLIB [19], which incorporates a semivariogram model for spatial correlation, to generate spatially correlated numbers, with which we generate porosity distribution as where is the spatially correlated number output from GSLIB and represent the dimensionless correlation lengths of direction and direction, respectively. is the coefficient of variation of porosity. Bigger means strong inhomogeneity; otherwise, it means slight inhomogeneity. is the average porosity.
Figures 2 and 3 show the porosity distributions for two layers with . The figures indicate that porosity distributions vary a lot with different geostatistic parameters so as to affect flow fields significantly.
3.2. PorePermeability Correlation
Acidrock reaction increases porosity, resulting in a change of permeability, which in turn affects acid flow. We use the following model to describe the relationships [20] between permeability and local porosity. where is a constant which is dependent on the structure of the medium.
3.3. Mineral Distribution
Using the same method for porosity distribution, we generate the spatial association distribution function which is defined as follows: where and are the average volume fractions of the fastreacting mineral and slowreacting mineral, respectively. Figure 4 is an example of the distributions of the volume fraction () of fast and slowreacting minerals.
4. Numerical Solution
The equations are discretized with the finite volume method and solved sequentially. Firstly we generate initial porosity and mineral distributions. Then, we solve equations (2) and (3) to obtain the pressure and velocity fields. Next, we solve equations (4) and (5) to get the acid concentration distribution. Finally, we update the mineral content and porosity distributions based the acidrock reaction (equations (6) and (7)).
5. Result and Discussion
In this section, extensive numerical simulations are conducted to analyze acid distribution patterns as well as the effect of heterogeneities on acidizing performance with the parameters in Table 1. And acidrock reaction parameters are determined through references [21–23].

5.1. Comparison of Homogenous and Heterogeneous Cases
The parameters used in the simulation are listed in Table 2. Other parameters are kept the same for the two cases.

Figure 5 shows acid concentration and mineral content distribution in the radial direction. Since they are different in different directions for the heterogeneous cases, only 0 degree (relative to the black line in Figure 5, for example) is chosen for displaying. VF_M1, VF_M2, and VF_M3 stand for fastreacting mineral, slowreacting mineral, and silica gel precipitate content. Compared to the homogeneous case, the curves of the heterogeneous case fluctuate apparently, the acid front is longer, and precipitation is deeper. Since the fastreacting mineral has a much higher reaction rate than the slowreacting mineral, most the fastreacting mineral is spent, but a small amount the slowreacting mineral is dissolved in the live acid zone.
(a) Homogeneous case
(b) Heterogeneous case
Figures 6–10 show acid and mineral concentration distribution after 90minute acid injection. The homogeneous case gives a regular shape of distribution. Meanwhile, the heterogeneous case gives rise to an irregular distribution shape. The porosity distribution is strongly correlated in the horizontal direction. The horizontal direction is the leading flow direction, so acid penetrates deeper in this direction. The acid fingers through the highperm zone, leaving the acid front to propagate nonuniformly. Some lowperm zones are bypassed by the acid. For the homogeneous case, from the wellbore to the outer zone, there exist zones of high perm, low perm, and original formation perm. Near the wellbore zone, the permeability is enhanced due to rock dissolution by the acid. In the live acid zone, most of the fastreacting mineral is dissolved. The lowperm zone is cause by the secondary precipitate of spent acid. Silica gel precipitation due the reaction between H_{2}SiF_{6} and the fastreacting mineral is inevitable [24], which is unfavorable for sandstone acidizing. For the heterogeneous case, the reduced perm zone is not distinctly apparently. The precipitate does exist, but the drastic porosity variation from one place to another makes it hard to identify. The significant difference of the penetration distance between the homogeneous and heterogeneous cases shows just how important it is to consider heterogeneity in field treatments. The acid penetration distance is much longer than the one predicted by the homogenous model, especially when there exists major flow channels.
5.2. Effect of Vertical Heterogeneity
A well often consists of multilayers with much different properties. For this kind of well, an important aspect in acidizing is acid placement. Here, the difference of the properties of the layers is called vertical heterogeneity. Vertical heterogeneity considers the difference of permeability and the mineral mainly in this simulation.
5.3. Effect of Vertical Permeability Heterogeneity
The vertical permeability difference of layers can be characterized by the permeability ratio, which is defined as the ratio of the average permeability of layer A to layer B as follows: where is the permeability ratio and and represent the average permeabilities of layer A and layer B, respectively. We did simulations for , 5, and 10, respectively. Parameters are listed in Table 3. Except for the different porosity and permeability distribution values, other parameters are kept the same. Results after 90minute acid injection are shown as follows.

As shown in Figure 11, vertical permeability heterogeneity has a notable effect on the acid distribution of the vertical profile. The higher the permeability ratio, the more seriously acid intake differs. When , the acid entering into the lowperm zone is negligible compared to that entering into the highperm zone. Except for the permeability difference which causes the discrepancy of acid intake, the permeability increase due to rock dissolution also magnifies the difference. The highpermeability zone takes more acid and raises permeability faster, so it accepts more and more acid. This demonstrates the necessity of acid diversion in acidizing multilayer formations.
(a)
(b)
(c)
Figure 12 shows fastreacting mineral dissolution in two layers. In line with the acid distribution, the dissolution is greater in the higher perm layer, and the dissolutions in the lower perm layer decreases with the increase of the permeability ratio. Figure 13 shows the distribution of precipitation. In the higher perm layer, the precipitate spreads wider because of more spent acid generated by the acid entering into the layer. Precipitate will decrease the permeability; more precipitate in the higher perm layer will counteract the permeability increase caused by rock dissolution, but due to areal heterogeneity, the acid can channel through the lowered perm zone.
(a)
(b)
(c)
(a)
(b)
(c)
5.4. Effect of Vertical Mineral Content Heterogeneity
When analyzing the effect of mineral content in the two layers, other parameters are kept the same. The fastreacting mineral of 0.1 and 0.2 and the slowreacting mineral of 0.85 and 0.75 are set for layers A and B, respectively.
Figures 14–16 show the results after acidizing for 90 minutes. Since the two layers have identical initial porosity and permeability distribution, the dissolution geometry is similar, but the extent of layer A is notably larger than layer B. Layer B has a much higher fastreacting content than layer A. Due to the relative fast reaction rate of the fastreacting mineral, most of the fastreacting mineral is dissolved in the area covered by the live acid. The acid front propagates much slower in layer B. Figure 17 shows the pump rate allocation with time for the two layers. The initial rate is the same. As acidizing proceeds, the rate of layer A increases and the rate of layer B decreases, since the total rate is fixed. There are two reasons for this phenomenon. The first is that the acid front moves faster in layer A than layer B. The second is more serious precipitation in layer B due to the higher content of the fastreacting mineral. Even though the vertical heterogeneity of the mineral content is obvious, its effect is less remarkable than the effect of permeability.
5.5. Effect of Areal Heterogeneity
Areal heterogeneity can be expressed by the coefficient of variation of porosity (). The higher the , the more serious the heterogeneity is. means homogeneous case. Figures 18–20 show the results for . The dimensionless correlation lengths are 0.08 and 0.0001 in the and directions. means very weak heterogeneity. The porosity distribution is very close to the homogeneous case. means very strong heterogeneity, which generates some leading flow paths similar to the channels in the reservoir. The more serious the heterogeneity, the more irregular the geometry of acid coverage is. The case of makes the acid finger through some highperm paths, making the acid penetration distance much longer than the other two cases. Another areal heterogeneity is natural fractures. For example, there is a fracture connecting to the wellbore as shown in Figure 21. The acid mainly flows through the fracture to reach a deeper distance. These indicate that in strong heterogeneous reservoirs, the acid penetration distance may be much deeper than the one predicted by the conventional homogeneous model.
6. Conclusion
In this paper, we built a sandstone matrix acidizing model accounting for areal and vertical heterogeneities and the effect of secondary precipitates. Based on the model, extensive numerical simulation was conducted to investigate the effect of heterogeneity on acidizing performance. The study reached the following conclusions: (1)The heterogeneities both in areal and vertical directions have a significant effect on acidizing performance. The areal heterogeneity influences the acid dissolution pattern and live acid penetration distance. Heterogeneous cases give longer live acid penetration distance than the homogeneous cases. The vertical heterogeneities determine the vertical acid intake profile(2)Both flow field and mineral heterogeneities make it possible to create highpermeability channels in acidizing and to obtain longer acid penetration distance. The heterogeneity of a flow field has a more serious effect on the acid dissolution pattern than the mineral heterogeneities(3)Serious vertical heterogeneities make an unfavorable acid intake profile, with the highpermeability zone taking excessive acid and the lowpermeability zone taking insufficient acid. In multilayer formations, heterogeneity creates acid fingering instead of a pistonlike pattern, which indicates the necessity of an acid diversion(4)The secondary precipitate has a significant effect on acidizing performance when fastreacting mineral content is high
Nomenclature
:  Hydrofluoric acid concentration 
:  Fluosilicic acid concentration 
:  Injection concentration of acid 1 
:  Injection concentration of acid 2 
:  Coefficient of variation of porosity 
:  Reaction rate constant 
:  Spatially correlated number output from GSLIB 
:  Average permeability of Layer A 
:  Average permeability of Layer A 
:  Dimensionless correlation lengths of direction and direction, respectively 
:  Mole weight of acid 
, :  Average volume fractions of fastreacting mineral and slowreacting mineral, respectively 
:  Pump rate 
:  Reservoir pressure 
:  Pressure at outer boundary 
:  Specific surface area of minerals 
:  Velocity vector 
:  Velocity in the radial direction 
:  Velocity in the circumferential direction 
:  Mineral volume fraction 
, , :  Initial distribution of the volume fractions of mineral 1, mineral 2, and mineral 3, respectively 
:  Reaction order 
:  Permeability ratio 
:  Gravimetric dissolving power 
:  Constant dependent on the structure of the medium 
:  Porosity 
:  Initial porosity distribution 
:  Average porosity 
:  Density. 
Data Availability
The numerical simulation data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the support of the National Science and Technology Major Project of a Research on Controlling Factors of TargetOriented Acid Fracturing and Experimental Study on Diversion, SINOPEC Northwest Oilfield Branch (2016ZX05014005012).
References
 C. W. Crowe, “Precipitation of hydrated silica from spent hydrofluoric acid: how much of a problem is it? (includes associated papers 16441 and 16444),” Journal of Petroleum Technology, vol. 38, no. 11, pp. 1234–1240, 1986. View at: Publisher Site  Google Scholar
 S. W. Almond, J. L. Brady, and D. R. Underdown, “Return fluid analysis from the Sadlerochit Formation, Prudhoe Bay, Alaska: A Laboratory studypart II,” in International Arctic Technology Conference, SPE 18223, Anchorage, Alaska, May 1991. View at: Publisher Site  Google Scholar
 D. R. Underdown, J. J. Hickey, and S. K. Kalra, “Acidization of analcimecemented sandstone, Gulf of Mexico,” in Annual Technical Conference and Exhibition, SPE 20624, New Orleans, LA, USA, September 1990. View at: Publisher Site  Google Scholar
 S. L. Bryant, “An improved model of mud acid/sandstone chemistry,” in Annual Technical Conference and Exhibition, SPE 22855, Dallas, TX, USA, October 1991. View at: Publisher Site  Google Scholar
 R. Gdanski, “Fractional pore volume acidizing flow experiments,” in European Formation Damage Conference, SPE 30100, The Hague, Netherlands, May 1995. View at: Publisher Site  Google Scholar
 R. D. Gdanski, “Kinetics of the primary reaction of HF on aluminosilicateS,” SPE Production & Facilities, vol. 15, no. 4, pp. 279–287, 2000. View at: Publisher Site  Google Scholar
 Y. Hekim, H. S. Fogler, and C. C. Mccune, “Radial movement of permeability fronts and multiple reaction zones in porous media: Soc Pet Engr J, V22, N1, Feb 1982, P99107,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 19, no. 6, pp. 130–130, 1982. View at: Publisher Site  Google Scholar
 J. E. McElhiney, A. L. Schalge, R. S. McKnight, and J. R. Robinson, “Coreflood acidization: model, measurement, and methodology,” in Oilfield and Geothermal Chemistry Symposium, SPE 7891, Houston, TX, USA, January 1979. View at: Publisher Site  Google Scholar
 S. D. Sevougian, L. W. Lake, and R. S. Schechter, “A new geochemical simulator to design more effective sandstone acidizing treatments,” SPE Production & Facilities, vol. 10, no. 1, pp. 13–19, 1995. View at: Publisher Site  Google Scholar
 B. F. Wang, “Study of sandstone acidification design model,” Henan Petroleum, vol. 2, pp. 28–31, 1999. View at: Google Scholar
 C. Li, T. Xie, M. Pournik, D. Zhu, and A. D. Hill, “Finescale simulation of sandstone acidizing,” Journal of Energy Resources Technology, vol. 127, no. 3, pp. 225–232, 2005. View at: Publisher Site  Google Scholar
 S. Y. Li, Z. M. Li, and B. F. Li, “Acid/sandstone reaction modeling for sandstonematrix acidizing,” Journal of China University of Mining & Technology, vol. 41, no. 2, pp. 236–241, 2012. View at: Google Scholar
 V. H. Leong, H. B. Mahmud, M. C. Law, H. C. Y. Foo, and I. S. Tan, “A numerical modelling and simulation of corescale sandstone acidizing process: a study on the effect of temperature,” Journal of Petroleum Exploration and Production Technology, vol. 9, no. 1, pp. 483–516, 2018. View at: Publisher Site  Google Scholar
 V. H. Leong, H. B. Mahmud, M. C. Law, H. C. Y. Foo, and I. S. Tan, “A comparison and assessment of the modelling and simulation of the sandstone matrix acidizing process: a critical methodology study,” Journal of Natural Gas Science and Engineering, vol. 57, pp. 52–67, 2018. View at: Publisher Site  Google Scholar
 K. R. Kunze and C. M. Shaughnessy, “Acidizing sandstone formations with fluoboric acid,” SPE Journal, vol. 23, no. 1, pp. 65–72, 1983. View at: Publisher Site  Google Scholar
 M. U. Shafiq, H. K. B. Mahmud, and M. Arif, “Mineralogy and pore topology analysis during matrix acidizing of tight sandstone and dolomite formations using chelating agents,” Journal of Petroleum Science and Engineering, vol. 167, pp. 869–876, 2018. View at: Publisher Site  Google Scholar
 R. L. Thomas and C. W. Crowe, “Matrix treatment employs new acid system for stimulation and control of fines migration in sandstone formations,” Journal of Petroleum Technology, vol. 33, no. 8, pp. 1491–1500, 1981. View at: Publisher Site  Google Scholar
 H. H. Hardy and R. A. Beier, “Integration of large and smallscale data using Fourier transforms,” in ECMOR IV  4th European Conference on the Mathematics of Oil Recovery, Roros, Norway, 1994. View at: Publisher Site  Google Scholar
 C. V. Deutsch and A. G. Journel, “GSLIB: Geostatistical Software Library and User’s Guide,” Tech. Rep., Oxford University Press, New York, Oxford, 1998. View at: Google Scholar
 M. K. R. Panga, M. Ziauddin, and V. Balakotaiah, “Twoscale continuum model for simulation of wormholes in carbonate acidization,” AICHE Journal, vol. 51, no. 12, pp. 3231–3248, 2005. View at: Publisher Site  Google Scholar
 H. R. Gates and J. H. Lybarger, “SGMA performance, design. Conclusion. Proper reservoir choice, fluid design keys to SGMA success,” Oil & Gas Journal, vol. 76, no. 43, 1978. View at: Google Scholar
 N. Kume, R. V. Melsen, L. Erhahon, and A. F. Abiodun, “New HF acid system improves sandstone matrix acidizing success ratio by 400% over conventional mud acid system in Niger Delta Basin,” in Annual Technical Conference and Exhibition, SPE56527, pp. 185–199, Houston, TX, USA, October 1999. View at: Publisher Site  Google Scholar
 A. N. Martin and K. L. Smith, “New HF acid system produces significant benefits in Nigerian sandstones,” in Annual Technical Conference and Exhibition, SPE38595, San Antonio, TX, USA, October 1997. View at: Publisher Site  Google Scholar
 H. Xue, L. Zhao et al., “Acidification simulation for inhomogeneous pore type sandstone oil reservoir,” Oil Drilling & Production Technology, vol. 37, no. 4, pp. 100–104, 2015. View at: Google Scholar
Copyright
Copyright © 2019 Jianye Mou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.