Next Article in Journal
Poly-arginine-18 (R18) Confers Neuroprotection through Glutamate Receptor Modulation, Intracellular Calcium Reduction, and Preservation of Mitochondrial Function
Next Article in Special Issue
Analysis of the Stochastic Quarter-Five Spot Problem Using Polynomial Chaos
Previous Article in Journal
Hybrid Nanocellulose-Copper (II) Oxide as Engine Oil Additives for Tribological Behavior Improvement
Previous Article in Special Issue
Exposure Time Impact on the Geomechanical Characteristics of Sandstone Formation during Horizontal Drilling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Approach to Improve Acid Diversion in Carbonate Rocks Using Thermochemical Fluids: Experimental and Numerical Study

by
Mustafa Ba Alawi
1,
Amjed Hassan
1,
Murtada Saleh Aljawad
1,*,
Muhammad Shahzad Kamal
1,
Mohamed Mahmoud
1 and
Ayman Al-Nakhli
2
1
College of Petroleum Engineering & Geosciences, King Fahd University of Petroleum & Minerals, Dhahran 31261, Saudi Arabia
2
Advanced Research Center (EXPEC ARC), Dhahran 31311, Saudi Arabia
*
Author to whom correspondence should be addressed.
Molecules 2020, 25(13), 2976; https://doi.org/10.3390/molecules25132976
Submission received: 1 June 2020 / Revised: 24 June 2020 / Accepted: 26 June 2020 / Published: 28 June 2020
(This article belongs to the Special Issue Fluids in Porous Media)

Abstract

:
The distribution of acid over all layers of interest is a critical measure of matrix acidizing efficiency. Chemical and mechanical techniques have been widely adapted for enhancing acid diversion. However, it was demonstrated that these often impact the formation with damage after the acid job is completed. This study introduces, for the first time, a novel solution to improve acid diversion using thermochemical fluids. This method involves generating nitrogen gas at the downhole condition, where the generated gas will contribute in diverting the injected acids into low-permeability formations. In this work, both lab-scale numerical and field-scale analytical models were developed to evaluate the performance of the proposed technique. In addition, experimental measurements were carried out in order to demonstrate the application of thermochemical in improving the acid diversion. The results showed that a thermochemical approach has an effective performance in diverting the injected acids into low-permeability rocks. After treatment, continuous wormholes were generated in the high-permeability rocks as well as in low-permeability rocks. The lab-scale model was able to replicate the wormholing impact observed in the lab. In addition, alternating injection of thermochemical and acid fluids reduced the acid volume 3.6 times compared to the single stage of thermochemical injection. Finally, sensitivity analysis indicates that the formation porosity and permeability have major impacts on the acidizing treatment, while the formations pressures have minor effect on the diversion performance.

1. Introduction

Matrix acidizing in carbonates aims to both remove damage near the wellbore area and further improve its permeability by creating conductive flow pathways (wormholes) around the wellbore [1,2,3]. Acid treating fluids designs have diverse compositions and formulations directed towards the targeted formation nature [4,5]. In addition, the acid fluids must be injected below the fracturing pressure to avoid the reservoir to be fractured and prevent the stimulating fluids from being lost consequently [6]. However, the conventional acidizing treatments are less effective with high permeability contrast between zones with large thickness [7,8,9]. With high degree of heterogeneity, it is required that the injected acid treatment is effectively distributed over both the less and more damaged layers [10,11]. Hence, the industry has adapted widely both mechanical and chemical diverters depending on the effectiveness of each approach on the different lithology [12,13].
Mechanical diverters are mainly utilized in sandstone reservoirs and occasionally used for carbonate formations [5,9]. The common mechanical diverters are ball sealers, coiled tubing (CT), and inflow control devices and valves [9,14,15]. Among all the mechanical methods, CT has become a very advantageous tool for enhancing stimulation fluid placement [9]. Acid can be injected through the tubing and terminated conveniently if the treatment is not performing as expected. Additionally, since the size of the tubing is small, the stimulating fluids can be quickly displaced. However, due to the small sizes of CT compared to production tubing and drill pipe, the injection rates can be limited and restricted due the high frictional losses. In addition, corrosion is always with acidizing applications using CT and it is considered disastrous when it is encountered [9]. The mechanical methods are known for effectiveness; however, high costs and risks are always associated whenever they are introduced into the wellbore. Additionally, mechanical diverters such as ball sealers and solid diverting agents are considered as temporary solutions due to their direct effects on downhole temperature that leads them to eventually dissolve.
In addition, there are many designed chemical diverters that are being utilized in the industry [5,7,16]. The injected materials are designed to be present high viscosity contrast to clog the formed wormholes and divert the acid to untreated zones. The widely used chemical diverters are emulsified acids, foams, and in situ gelled acid-based polymers or surfactants [17,18,19,20,21]. Emulsified acids were used to improve the acidizing treatment in high temperature reservoirs by reducing the fluid loss and increasing the acid penetration [22,23]. This has been served by diesel emulsified acid (DEA) which consists of acid, diesel, and an emulsifying agent. The injected emulsion will encroach into the high-permeability layers and divert the acids into the low-permeability layers, utilizing the viscosity contrast in different layers [17].
Moreover, acids can be foamed by introducing gas or by slugging foams into them [16]. Foam enables the placement of a large volume of gas in the formation that is later trapped when acid is injected. Consequently, zones with high gas saturation will clog the pore throats and will require high pressure to be mobilized. Thus, these zones will have low liquid relative permeability that results in acid diverting from them [7]. Recently, nitrogen gas has been adapted for this technique as the discontinuous phase, more so than carbon dioxide, for both economic and environmental reasons. Foam quality is attributed to the volumetric fraction of its gas content, and is an important factor to optimize to attain successful acid diverting performance [21,24]. However, foams require surfactants to keep them stable. In addition, they require the formation to be injected with surfactants as well to ensure its stability in situ [24]. In oil reservoir cases, oil can break the foam stability, so in many cases, it is advised to pump a mutual solvent ahead of the foam [21]. Overall, the main limitation of using foam as an acid diverter is that the foam quality will be significantly reduced when it gets to the wellbore. Therefore, in situ foam generation could be a good solution in order to preserve the advantages of foam treatment and avoid the limitations of foam stability.
Furthermore, viscoelastic surfactant (VES) can be used to divert the injected acid by increasing the fluid viscosity upon the reaction of acid with the formation matrix [20]. The high-viscosity system will block the pathways made inside the rock by the acid and divert it into the uninvaded area to effectively act on the low-permeability zone [25]. Viscoelastic surfactants exhibit efficient and clean wormhole creation inside the treated formation and do not require any flowback processes; solvents or post-flush can be used to retract the viscosity to the gelled acid. However, the effect of this system in causing damage inside the created pathways is possible. In addition, the solvents and VES viscosity contrast can lead to fingering and cause poor formation clean-up [19].
Overall, mechanical approaches have shown good performance in acid diversion, however, they may require high operational cost. In addition, deep or high temperature reservoirs may restrict the application of mechanical diverters due to corrosion issues. On the other hand, chemical diverters provide effective acid placement solutions corresponding to the formation conditions. However, they still have either high potentiality of causing severe damage to the formation or issues in stability; both require further pre- and post-processing. Injection of gas prior to the acid can help significantly in achieving effective acid diversion. However, the major concern is that gas injection may be associated with high operational costs due to gas transportation and treatment. Therefore, in situ gas generation at the reservoir condition can preserve the advantages of using gas injection for acid diversion applications. In addition, this will minimize the operational cost since the gas will be generated only in the wellbore. Similarly, in situ foam generation can result in better foam quality compared to foam injection from the surface.
This study presents an effective method to place acid uniformly in the damaged formations using in situ generated gas by thermochemical reactions. This method provides the advantage of not requiring an external source of gas for diversion. In addition, the induced pressure from the thermochemical reaction will force the generated gas to flow into the high-permeability zones, thereby improving the performance of acid diversion treatment. In this work, experimental and modeling studies were carried out in order to evaluate the efficiency of the proposed technique. The experiments work constitutes of parallel coreflooding tests using rock samples that have similar geometry and petrophysical properties but different permeability values to simulate heterogeneous formations. A new analytical model was developed to predict the performance of acid diversion using thermochemical fluids. Finally, sensitivity analysis was conducted to assess the impact of formation parameters on the treatment efficiency. The influences of formation porosity, layer permeability, and pressure difference between the two layers on the acid diversion performance were studied.

2. Results and Discussion

2.1. Impact of Acid Diverter on Acidizing Performance

2.1.1. Experimental Outcomes

Acidizing experiment was conducted using two rock samples with permeability contrast to study the acid distribution when the acid is injected solely without any diverter. Two carbonate rocks were treated with 15 wt% HCl acid, and the used samples had initial permeability values of 90 and 27 mD, respectively. Parallel coreflood system was used, and the two rock samples were treated using the same experimental conditions, with a similar injection pressure and acid flow rate applied to the two rock samples. However, when acid was injected, the majority of the acid invaded into the high-permeability core (90 mD) and only small amounts flowed into the lower permeability core (27 mD). Therefore, a continuous wormhole was generated in the high-permeability rock, while limited wormhole generation was induced in the low-permeability rock. Figure 1 shows the CT-scan images of the tested rocks after acidizing treatment. It should be noted that this acidizing treatment was conducted without using any acid diverter. Therefore, the propagation of wormholes in the high-permeability core achieved breakthrough successfully while it did not in the low-permeability core. Furthermore, acidizing treatment was performed, and thermochemical fluids were injected to witness its diversion performance. Two rock samples with original permeability of 94 and 26 mD were used in this experiment. Thermochemical fluids were injected into the parallel core system prior to acid treatment. In this case, the injected acid was diverted to the lower permeability core. Therefore, the high- and low-permeability rocks were successfully stimulated. Figure 2 shows CT-scan images for the treated rocks after the treatment using HCl acid and thermochemical fluids. Continuous wormholes were observed in the high- and low-permeability rocks, indicating the effective performance of thermochemical fluids in diverting the injected acids into the low-permeability rocks.

2.1.2. Model Match

The abovementioned experimental conditions were used as input in the model to replicate the impact of permeability contrast on acid treatment. The cross-sectional area of the core was modeled to be square (1.33 in edge length), which has a similar area to the 1.5 in diameter core. The model assumes that the two core samples were placed on top of each other with a thin low-permeability layer in the middle acting as a flow barrier. This is to account for the actual separation between the cores in the experimental setup. The porosity and permeability were randomly populated in the domain using Gaussian distribution. Both cores were subjected to the same 15 wt% HCl acid which was injected at 1 cc/min. Figure 3 shows the permeability distribution where the high-permeability core sample was placed on the top. The mean permeability of the top domain was 90 mD, while it was 27 mD for the bottom domain; this was to replicate the behavior in Figure 1. The porosity distribution of both samples is shown in Figure 4 where the mean porosity is 0.2. Porosity lab measurements of the high- and low-permeability rock samples were similar; hence, the same distribution was assigned during simulations. Some statistical parameters, such as the correlation length (Ix, Iy), were assumed to be similar for both porosity and permeability. The standard deviation determines how far the distribution is from the mean. The standard deviation for permeability distribution in the domain, σ k , was 3, while for porosity, σ φ was 0.04, as illustrated in the figure captions.
Figure 5 shows the final concentration profile after wormhole breakthrough assuming HCl acid was injected with no diversion. It is noted that wormhole propagated fully in the high-permeability core (the top domain of Figure 5) while only limited propagation was observed in the low-permeability core. Various simulations of different permeability contrast were conducted, and all showed that acid preferentially invaded the high-permeability core. This replicated the experimental outcomes obtained in Figure 1 which also built confidence in the fine-scale modeling approach.
Once the thermochemical fluids react, they generate a tremendous amount of nitrogen that invades the high-permeability zone, diverting acid to the low-permeability one. Modeling the diversion impacts of thermochemical fluids requires a two-phase flow model which is out of the current research scope. Nevertheless, the simulations in Figure 6 show that the wormhole propagated equally only when the permeability of the two cores to the HCl acid was within the same range. This indicates that the equal acid stimulation in Figure 2 was generated due to the thermochemical fluids’ ability to equate the relative permeability of the two cores.

2.2. Alternating Injection of Thermochemical and Acid Fluids

The thermochemical approach showed good performance in diverting the injected acids and then enhancing the efficiency of acidizing treatment. However, the treatment efficiency can be further improved by alternately injecting thermochemical and acid fluids. The alternating injection can lead to a reduced amount of injected fluids and it can also result in deeper acid penetration compared to the single-stage injection. Therefore, in this work, the technique of alternating injection of thermochemical and HCl fluids into the rock samples of different permeability (11 and 102 mD) was examined. A total of six injection cycles were applied: three thermochemical fluid and three HCl injection cycles. Figure 7 and Figure 8 show the pressure profiles for the low- and high-permeability rocks, respectively. The duration of each cycle was determined based on the stability of the pore pressure along the treated rock. During the thermochemical cycles, the pressures at the cores’ inlet and outlet are almost flat for all samples. However, during the HCl injection, the pressures fluctuate due to the reaction between the injected fluid and the carbonate matrix. At the final stage of acid injection, the pressure fluctuation is very high, which can indicate significant solid migration. Injecting HCl into carbonate rock can dissolve part of the carbonate matrix and induce solid migration [26]. Moreover, Figure 9 compares the pressure-drop profiles for the low- and high-permeability rocks during the alternating injection of thermochemical and HCl fluids. The pressure-drop profiles indicate that wormholes were generated in both rocks, since the final pressure-drop is almost zero. The wormhole was created in the high-permeability rock during the second cycle, while it was generated during the fourth injection cycles in the low-permeability rock. Ultimately, continuous wormholes, from the core inlet to the outlet, were achieved in both samples after the thermochemical/HCl treatment. Figure 10 and Figure 11 show images for the cores’ inlet and outlet after the alternating injection of thermochemical and HCl fluids.

2.3. Modeling of Acid Diversion Using Thermochemical Fluids

2.3.1. Single-Stage Injection

The developed model was used to assess the acidizing treatment in a heterogenous reservoir with two layers of different permeability values (20 and 100 mD). The injected acid will be diverted into the low-permeability layer utilizing the in situ gas generation technique. The description of wellbore and reservoir formations during this process are provided in Section 3.1 and Section 3.2. The properties of the reservoir, well, and fluid are listed in Table 1. The top layer is characterized by considerably high permeability compared to the bottom layer. Table 2 summarizes the properties of top formation (layers 1) and bottom formation (layer 2). The permeability of layer 1 is five times higher than the permeability of layer 2. Therefore, during the thermochemical injection, most of the generated gas will invade into layer 1. In this work, the performance of thermochemical fluids in diverting the injected acids was quantified by determining the flow rates into the top (Q1) and bottom (Q2) layers. Then, the ratio Q1/Q2 was used as an indicator of diversion efficiency. If the Q1/Q2 ratio equals one, this indicates that the injected fluids are equally distributed into the two layers. Meanwhile, if the flow rate ratio is less than one, this reveals that the injected acid was successfully diverted into the low-permeability layer [8], and the acid flow rate into the low-permeability layer (Q2) is higher than into the high-permeability layer (Q1).
Figure 12 shows the profiles of Q1/Q2 ratio as a function of the injected acid volume, with thermochemical reaction and N2 injection used as acid diverter. The flow rate ratio (Q1/Q2) is decreasing, with acid injection indicating that the injected acid is invading into the low-permeability layer and this will result in reducing the skin factor in the low-permeability layer. Figure 13 shows the formation skin for the top and bottom layers against the acid volume. This figure indicates that thermochemical reaction has a better diversion performance than nitrogen injection since the thermochemical reaction showed lower skin values compared to the N2 injection. Thermochemical approach can reduce the acid volume by around 23% compared to the N2 injection, and the required acid volume was reduced from around 22 to less than 18 ft3/ft. This can be attributed to the injection of thermochemical fluids providing more pressure difference and viscous contrast compared to the N2 injection. Ultimately, the thermochemical approach showed effective acid diversion as indicated by the reduction in the Q1/Q2 ratio to less than one, revealing that more acid is invaded into the low-permeability layer.

2.3.2. Alternating Thermochemical and Acid Injection

Alternating thermochemical and acid fluids injection can be used to improve the efficiency of acid diversion and minimize the amount of injected acids. The developed model was used to predict the performance of alternating thermochemical and acid injection. Figure 14 shows the reduction in skin factors in high- and low-permeability layers during the alternating thermochemical and acid injection. The alternating injection approach showed better performance compared to the single stage of thermochemical injection, with the treatment volume reduced from 20 to around 5.6 ft3/ft. Alternating the thermochemical fluids will advance the treatment by the advantages of increasing the pressure and temperature in every cycle and then reduce the acid volume by around 3.6 times compared to single-stage injection.

2.4. Sensitivity Analysis

Sensitivity analysis was conducted to assess the impact of formation parameters on the treatment efficiency. The influences of formation porosity, layer permeability, and pressure difference between the two layers on acid diversion performance were studied. Table 3 lists the model parameters of the different cases used in the sensitivity analysis work. Figure 15 shows the impact on the pressure difference on the acid diversion performance. The pressure difference between the two layers (low- and high-permeability layers) varied from 600 to 3100 psi, and the diversion performance was evaluated based on the Q1/Q2 ratio. Recalling that, a lower Q1/Q2 value indicates that the injected acids are invading into the low-permeability layers and better acid diversion is then obtained. The pressure difference showed a minor effect on the diversion performance. Similar profiles of Q1/Q2 were observed using different pressure values in the treated formations. This is caused by the major difference of viscosity between the reservoir liquid and injected gas. Thus, it has substantially overcome the effect of pressure difference between the reservoir and bottomhole ones. Although pressure difference is a major contributor to the injected fluids to encroach further to the wellbore adjacent layers, the model honors the viscosity contrast and adequately illustrates almost equal flowrates to both layers.
Additionally, Figure 16 examines the response of diversion performance with different porosity values. The highest porosity values exhibited the best diversion performance compared to the other two cases. The bigger porosity value with gas invading the top layer pores assisted in providing an additional drive for the acid to be diverted to the lower layer. By contrast, the smaller porosity value allowed for small volumes that acid can saturate and, consequently, the diversion effect was slightly less. Moreover, the effect of layer permeability on the acid diversion performance is shown in Figure 17, where permeability ratios of 5 and 50 were used. As the results convey, the model manifested permeability as the major driving parameter for acid flow between the two layers. As the permeability of the top layer increases, it accumulates most of the injected acid correspondingly. This has also resulted in delaying the treatment of the bottom layer massively in comparison to the case with less permeability value of the top layer. Finally, while the bottom layer permeability was decreasing down to the tight levels, and the diversion effect responded with weaker performance. Additionally, the small permeability draws slight volumes of acid resulting in the top layer drawing most of this. As an ultimate response, the top layer gets a significant lead to be treated ahead of the bottom layer.

3. Methods and Materials

3.1. Description of the Proposed Approach

Thermochemical fluids are used to improve acid diversion, utilizing the in situ gas generation mechanism. Figure 18 shows a schematic of the acid diversion using thermochemical fluids, where two layers of different permeability are shown (K1 and K2). Initially, the two layers are saturated with oil. Then, thermochemical fluids are injected to react downhole and generate nitrogen gas. Due to the flow resistance, most of the generated gas will flow into the high-permeability layer while a small amount will flow into the low-permeability zone. The chemical reaction can be triggered by the reservoir temperature and acid can also be injected to reduce the solution pH and accelerate the thermochemical reaction [27,28]. Thereafter, acid is injected to stimulate the reservoir layers and reduce the formation skin. During the acidizing treatment, the generated nitrogen will block the high-permeability zones and divert the injected acids into the low-permeability zones. Consequently, both the high-permeability and low-permeability zones will be stimulated.
The thermochemical fluids reaction can be represented by Equation (1) [27,28]:
N H 4 C l + N a N O 2 N a C l + 2 H 2 O + N 2   ( g ) + Δ H ( h e a t )
Nitrogen gas (N2) will be generated in a form of chemically induced pulse that ought to significantly increase the pressure of the invaded zones, and the generated pressure can reach up to 5000 psi [29]. The induced pressure in the wellbore will force the generated gas to flow into the reservoir formations. Due to the flow resistance, most of the generated gas will flow into the high-permeability zones while small volume of N2 gas will invade into the low-permeability zones. Hence, the generated gas will promote the acid to be diverted toward the low-permeability layers for which treatment is desired. Additionally, the thermochemical reaction will increase the temperature substantially, and could reach up to 600 °F [29], which will reduce the viscosity of the formation fluids. This will also promote the viscous contrast diversion to distribute the acid efficiently over the formation zones. However, the thermochemical treatment should be precisely designed in order to optimize the amounts of generated gas and to ultimately attain successful treatment. The volumes and concentrations of thermochemical fluids should be accurately selected to ensure that the generated pressure and temperature will be within the specifications of downhole completions.

3.2. Model Development

3.2.1. Numerical Fine-Scale Model

A finite volume based average continuum, also known as the two-scale model, was developed in this research to observe the impact of permeability contrast on acid distribution [30]. It is assumed that at any point in the domain, both fluids and solids exist in a proportion defined by the porosity. The elements in the domain are large enough to contain both solids and liquids but small enough to have homogenous properties within the element. A correlated random Gaussian distribution of porosity and permeability was implemented to ensure a heterogonous system within the full domain.
First, the continuity equation is implemented which conserve mass in the liquid phase, expressed as
· ( ρ f u ) = ( ϕ ρ f ) t
where ρ f is the fluid density, u is the velocity vector, ϕ is the porosity, and t is time. The superficial velocity defined by Darcy’s law is written as
u = k μ f · p
where k is the permeability tensor,     μ f is the fluid viscosity, and p is the pressure. The fluid superficial velocity is obtained from solving the system of equations above. Then, the convective, diffusive, and reactive transport equation of acid in porous media is solved:
( ϕ ρ f C A ) t + . ( ρ f u C A ) = . ( ϕ ρ f D e C A ) + a v ( k r k c ) k r + k c C A
where C A is the HCl acid concentration represented in mass fraction, D e is the effective acid diffusion, k r is the reaction rate constant, k c is the mass transfer coefficient, and a v is the carbonate mineral-specific area. The first term in the equation above represents acid accumulation, the second term represents acid convection, the third term represents acid diffusion, and the last term represents acid/rock reaction. The acid/rock reaction causes the porosity of the rock to increase according to the following:
( ϕ ) t = a v ( k r k c ) k r + k c C A β 100 ( ρ f / ρ r )
where β 100 is the dissolving power of the full strength HCl acid, and ρ r is the rock density.

Constitutive Equations

The algorithm used in this model is sequential, which starts by solving Equations (2) and (3) for the velocity profile, Equation (4) for the acid concentration, and Equation (5) for acid/rock dissolution. To solve the system above, some properties should be updated at each time step, such as k c , D e , a v , and k . The effective acid diffusion in porous media can be related the molecular diffusion,   D m , through the following:
D e = ϕ D m
The mass transfer coefficient can be obtained through relating Schmidt number,   S c , Reynold number of flow in pores, R e p , and Sherwood number,   S h . The following are the definitions of these numbers:
S h = S h + b R e p 1 2   S c 1 3
R e p = 2 ρ f r p | u | μ f
Sc = μ f ρ f D m
where r p is the pore radius, S h is the asymptotic Sherwood number, and b is a constant. The mass transfer coefficient can be obtained by solving the equation below once Sherwood number is determined from Equation (7).
S h = 2 k c r p D m
The equations below are used to update the pore radius, permeability, and specific surface area based on the change in porosity [31]:
r p = r p o ( ϕ ϕ o ) ε
k = k o [ ( ϕ ϕ o ) 3 ( 1 ϕ o 1 ϕ ) 2 ] γ
a v = a v o [ ( ϕ ϕ o ) ( 1 ϕ o 1 ϕ ) ] η
where ε , γ , and η are empirical parameters that could be used to tune the model with experiment outcomes.

Initial and Boundary Conditions

To solve Equations (2) and (4), initial and boundary conditions should be implemented. The following conditions were implemented for Equation (2):
p ( x , y , z ,   t = 0 ) = p i I n i t i a l l y
( p x ) ( x = 0 , y , z , t ) = q B μ f k x A c a t   t h e   i n l e t p ( x = L , y , z ,   t ) = p o u t a t   t h e   o u t l e t n . p = 0 a t   t h e   o t h e r   b o u n d a r i e s
where p i is the initial pressure, p o u t is the outlet pressure, q is the injection rate, B is the formation volume factor,   A c is the grid cross-sectional area, L is the core length, and n is the normal vector to the surface. The conditions for Equation (4) are as follows:
C A ( x , y , z ,   t = 0 ) = 0 I n i t i a l l y
C A ( x = 0 , y , z ,   t ) = C i n i a t   t h e   i n l e t n . C A = 0 a t   t h e   o t h e r   b o u n d a r i e s
where C i n i is the initial acid concentration.

3.2.2. Skin Revolution Model

An analytical model was developed to evaluate the performance of using thermochemical fluids to improve the acidizing treatment in heterogenous reservoirs. MATLAB program was developed to estimate the skin factor for the permeable and tight zones around the wellbore. The skin factor for each zone will be determined as a function of the injected acids, assuming incompressible and immiscible fluids. In addition, due to the substantial mobility ratio of the generated gas with comparison to the injected acid, the flow of acid is assumed to be acting with piston-like displacement. Moreover, it is assumed that nitrogen gas is encroaching the permeable layer in significant amounts, while negligible amounts of the nitrogen gas will invade into the low-permeability layer. Finally, the layer skin is used as an indicator of the acidizing performance, with the layer skin being reduced with higher volumes of acid invading into the treated layer. The equations used to develop the acid diversion model are described below.
The pressure difference due to formation skin can defined as follows [32]:
Δ p s = Δ p r e a l Δ p i d e a l
where Δ p s is the pressure drop due to skin, Δ p r e a l is the real or actual pressure drop, and Δ p i d e a l is the ideal pressure drop where no damage was induced. The pressure drop can be determined assuming steady state flow as follows:
Δ P = q μ 2 π k h ln ( r e r w )  
where q is the flow rate, μ is the fluid viscosity, k is the formation permeability, h is the formation thickness, re is the drainage radius, and rw is the wellbore radius. Figure 19 shows the different regions around the wellbore during acidizing treatment.
Considering the regions around the wellbore that are saturated with spacer and acids fluids, the pressure drop for real and ideal cases can be expressed as:
Δ P r e a l = q μ s 2 π k h ln ( r s r a ) + q μ a 2 π k h ln ( r a r w )
Δ P i d e a l = q μ 2 π k h ln ( r s r w )
where rs is the distance to which spacer penetrated, ra is the acid-saturated zone,   μ s is the viscosity of spacer fluid, and μ a is the acid viscosity. Substituting Equations (16) and (17) in Equation (18), the pressure drop due to formation skin can be given by
Δ P s = Δ P r e a l Δ P i d e a l = q μ s 2 π k h ln ( r s r a ) + q μ a 2 π k h ln ( r a r w ) q μ 2 π k h ln ( r s r w )
Since the treated zones have similar thickness and flow rate, the pressure drop can be expressed as
Δ P s = q μ 2 π k h [ μ s μ ln ( r s r a ) + μ a μ ln ( r a r w ) ln ( r s r w ) ]
Then, the skin factor ( S v i s ) can be defined as
S v i s = [ μ s μ ln ( r s r a ) + μ a μ ln ( r a r w ) ln ( r s r w ) ]        
and the pressure drop due to formation skin will be given as
Δ P s = q μ 2 π k h S v i s .
Moreover, in this work, the diversion effect was quantified as a viscous skin factor Svis for gas with viscosity µg that is being displaced by an acid with viscosity µacid as the following [32]:
S v i s = ( μ a c i d μ g 1 ) ln r a c i d r w .
The radius of acid penetration racid towards the formation targeted zone can be calculated as
V a c i d = π ϕ h   ( r a c i d 2 r w 2 )
where Vacid is the volume of injected acid,   ϕ is the formation porosity, and h is the formation thickness. Then, the radius of acid penetration (racid) can be determined by rearranging Equation (25).
r a c i d = r w 2 + V ¯ j π ϕ
where V ¯ j is volume per unit thickness (V/h) of acid invaded into a particular zone (j). Equation (26) can be substituted in Equation (24) to get
S v i s = ( μ a c i d μ g 1 ) [ 0.5 ln ( r w 2 + V ¯ j π ϕ ) ln r w ]
Equation (27) indicates that the viscous skin factor ( S v i s ) is changing with the volume of the injected acid ( V ¯ j ). Assuming steady state flow, the acid flow rate into a particular layer (qi) can be determined as follows [8]:
q i h = 4.92 × 10 6 k ( P w P e ) μ ( ln ( r e r w ) + S + S v i s ) .
Then, for a certain zone, the cumulative volume of injected acid can be estimated utilizing the acid flow rate (qi) as follows:
V 1 = V 0 + ( q i h ) 0 ( t 1 t 0 )
Acid-injected volumes, in turn, will be increasing until the removal of all damage, at which point the formation skin (S) can be deduced depending on the acid volume (Vc) [8]:
S = S 0 C   V ¯ c
where S0 is the original skin, C is the rate of skin decrease, and V ¯ c is the cumulative injected acid at a particular time.
Ultimately, the developed model can estimate the formation skin in each layer utilizing the radius of the acid and spacer zones around the wellbore as presented in Equation (22). In addition, the developed model can estimate the formation skin as a function of the injected acid and the original formation skin (before treatment) as presented in Equation (30). Therefore, MATLAB program was developed to estimate the skin factor using Equations (22) and (30).

3.3. Experimental Measurements

Coreflooding experiments were conducted in order to determine the performance of the developed approach of using thermochemical fluids in diverting the injected acids. A parallel coreflood setup was used, and carbonate rock samples of different permeability values were used. The acidizing experiments were conducted under high pressure and high temperature conditions. An overburden pressure of 1000 psi and backpressure of 400 psi were applied. In addition, HCl of 15wt% concentration was used, the injection rate was maintained at 1 cm3/min, and the system temperature was 80 °C. Moreover, the rock samples were prepared by measuring the rock permeability and porosity. Then, the treatment was implemented using alternating injection of HCl acid and thermochemical fluids. The pressure drop across the high- and low-permeability rock samples were monitored using pressure transducers at the core’s inlet and outlet. The properties of the used rock samples are listed in Table 4, and Figure 20 shows the experimental setup used in this work.

4. Conclusions

This study presents a novel solution to improve acid diversion utilizing in situ nitrogen generation. This method involves injecting thermochemical fluids during the acidizing treatment. In this work, numerical and analytical models were developed, and experimental measurements were carried out in order to evaluate the performance of the proposed technique. Based on this study, the following conclusions were drawn:
  • Injecting thermochemical fluids before acidizing treatment showed an effective performance in diverting the injected acids into the low-permeability rocks and then enhancing the efficiency of acidizing treatment.
  • After the treatment, CT-scan images showed that continuous wormholes were generated in the high-permeability rocks as well as in low-permeability rocks.
  • The two-scale continuum model was able to reproduce the experimental outcomes of wormhole propagation.
  • Alternating injection of thermochemical and acid fluids showed better acidizing performance compared to the single stage of thermochemical injection.
  • Alternating the thermochemical fluids will advance the treatment through the advantages of increasing the pressure and temperature in every cycle and then reduce the acid volume by around 3.6 times compared to single-stage injection.
  • The thermochemical approach showed better diversion performance compared to nitrogen injection, and lower values of formation skin were obtained using the thermochemical fluids.
  • Injecting thermochemical fluids can provide greater pressure difference and viscous contrast compared to the N2 injection. Hence, a thermochemical approach can reduce the required acid volume by around 23% compared to N2 injection.
  • Sensitivity analysis indicates that the formation porosity and permeability have great impact on the acidizing treatment, while the pressure difference between the different layers has a minor effect on diversion performance.

Author Contributions

Conceptualization, M.M and M.S.A.; methodology, M.B.A., A.H., M.S.A., M.S.K., M.M. and A.A.-N.; validation, A.H., M.S.A., and M.M.; formal analysis, M.B.A., A.H., M.S.A., and M.M.; investigation, M.S.A., M.M. and A.A.-N.; data curation, M.S.A., M.M. and A.A.-N.; writing—original draft preparation, M.B.A., A.H., and M.S.A.; writing—review and editing, M.S.A., and M.M.; visualization, M.S.A., M.S.K., M.M., and A.A.-N.; supervision, M.S.A., and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The college of Petroleum and Geoscience at King Fahd University of Petroleum & Minerals is acknowledged for the support and permission to publish this work.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

A c Cross-sectional area, L2, ft2 [m2]
a v Minerals specific surface are, L2/L3, ft2/ft3 [m2/m3]
B Formation volume factor, dimensionless, ft3/scf [m3/m3]
β 100 Acid gravimetric dissolving power
C A Acid concentration, wt%, mol/L3, lbmol/ft3 [mol/m3]
C i n i Initial acid concentration, wt%, mol/L3, lbmol/ft3 [mol/m3]
CRate of skin decrease for layer j, ft/ft3 [m/m3]
D e Effective acid diffusion, L2/t [cm2/s]
I Correlation length, dimensionless, dimensionless
k Permeability tensor, L2, mD [m2]
k r Reaction rate constant, kg-mol mineral/m2/s/(kg∙mol HCl/m3 acid)n
k c Mass transfer coefficient
hFormation or layer thickness, L, ft [m]
p Pressure, M/Lt2, psia [Pa]
pwWellbore pressure, M/Lt2, psia [Pa]
pePressure at drainage drainage, M/Lt2, psia [Pa]
qiInjection rate, L3/t, bbl/day [m3/day]
racidAcid penetration radius, L, ft [m]
reDrainage radius, L, ft [m]
RepFlow in porous media Reynold number
rwWellbore radius, L, ft [m]
rpPore radius, L, microns [µm]
ScSchmidt number
ShSherwood number
SvisViscous skin factor
u Velocity vector, L/t, ft/min [m/s]
Vacid, aAcid volume, L3, gal [m3]
V ¯ c Acid volume to reduce skin to zero, L3/L, ft3/ft [m3/m]
VtTotal acid volume, L3/L, gal/ft [m3/m]
ϕPorosity
ε Empirical parameter for pore radius, dimensionless
η Empirical parameter for surface area, dimensionless
γ Empirical parameter for permeability, dimensionless
σ Standard deviation
μ f Fluid viscosity, M/Lt, cp [Pa-s]
ρ f Fluid density, M/L3, lbm/ft3 [kg/m3]
µacid, aAcid viscosity, M/Lt, cp [Pa-s]

References

  1. Alghamdi, A.; Nasr-El-Din, M.A.; Hill, A.D.; Nasr-El-Din, H.A. Diversion and propagation of viscoelastic surfactant-based acid in carbonate cores. In Proceedings of the SPE International Symposium on Oilfield Chemistry, The Woodlands, TX, USA, 20–22 April 2009. [Google Scholar]
  2. Abdelgawad, K.Z.; Mahmoud, M.; Elkatatny, S.; Abdulraheem, A.; Patil, S. Reaction Kinetics and Coreflooding Study of High-Temperature Carbonate Reservoir Stimulation Using GLDA in Seawater. Energies 2019, 12, 3407. [Google Scholar] [CrossRef] [Green Version]
  3. Yuan, T.; Wei, C.; Zhang, C.S.; Qin, G. A Numerical Simulator for Modeling the Coupling Processes of Subsurface Fluid Flow and Reactive Transport Processes in Fractured Carbonate Rocks. Water 2019, 11, 1957. [Google Scholar] [CrossRef] [Green Version]
  4. Hansen, J.H.; Nederveen, N. Controlled Acid Jet (CAJ) technique for Effective Single Operation Stimulation of 14,000+ ft long reservoir sections. In Proceedings of the European Petroleum Conference, Society of Petroleum Engineers, Aberdeen, Scotland, UK, 29–31 October 2002. [Google Scholar]
  5. Nasr-El-Din, H.A.; Hill, A.D.; Chang, F.F.; Sultan, A.S. Chemical diversion techniques used for carbonate matrix acidizing: An overview and case histories. In Proceedings of the International Symposium on Oilfield Chemistry, Houston, TX, USA, 28 February–2 March 2007. [Google Scholar]
  6. Paccaloni, G.; Tambini, M. Advances in matrix stimulation technology. J. Petro. Tech. 1993, 45, 256–263. [Google Scholar] [CrossRef]
  7. Kalfayan, L.; Anthony, N. The Art and Practice of Acid Placement and Diversion: History, Present State, and Future. In Proceedings of the SPE Annual Technical Conference and Exhibition, New Orleans, LA, USA, 4–7 October 2009. [Google Scholar]
  8. Fadele, O.; Zhu, D.; Hill, A.D. Matrix acidizing in gas wells. In Proceedings of the SPE/CERI Gas Technology Symposium, Calgary, AB, Canada, 3–5 April 2000. [Google Scholar]
  9. Bybee, K. Acid Placement and Diversion. J. Petro. Tech. 2010, 62, 53–54. [Google Scholar] [CrossRef]
  10. Wang, D.; Sun, J.; Li, Y.; Peng, H. An Efficient Hybrid Model for Nonlinear Two-Phase Flow in Fractured Low-Permeability Reservoir. Energies 2019, 12, 2850. [Google Scholar] [CrossRef] [Green Version]
  11. Taylor, K.C.; Nasr-El-Din, H.A. Laboratory evaluation of in-situ gelled acids for carbonate reservoirs. SPE J. 2003, 8, 426–434. [Google Scholar] [CrossRef]
  12. Huang, T.; Hill, A.D.; Schechter, R.S. Reaction rate and fluid loss: The keys to wormhole initiation and propagation in carbonate acidizing. SPE J. 2000, 5, 287–292. [Google Scholar] [CrossRef]
  13. Hung, K.M.; Hill, A.D.; Sepehrnoori, K. A mechanistic model of wormhole growth in carbonate matrix acidizing and acid fracturing. J. Petro. Tech. 1989, 41, 59–66. [Google Scholar] [CrossRef]
  14. Erbstoesser, S.R. Improved ball sealer diversion. J. Petro. Tech. 1980, 32, 1–903. [Google Scholar] [CrossRef]
  15. Hajjar, H.; Al-Azizi, B. Use of Sliding Sleeve Valves in Long Laterals for Control of Production, Injection and Stimulation. In Proceedings of the SPE/IADC Middle East Drilling Technology Conference and Exhibition, Abu Dhabi, UAE, 26–28 January 2016. [Google Scholar]
  16. Kovscek, A.R.; Radke, C.J. Fundamentals of foam transport in porous media. Adv. Chem. 1994, 242, 115–163. [Google Scholar]
  17. Buijse, M.A.; Van Domelen, M.S. Novel application of emulsified acids to matrix stimulation of heterogeneous formations. In Proceedings of the SPE Formation Damage Control Conference, Lafayette, LA, USA, 26–28 January 1998. [Google Scholar]
  18. Eoff, L.S.; Dalrymple, E.D.; Reddy, B.R. Development of associative polymer technology for acid diversion in sandstone and carbonate lithology. SPE Prod. Facil. 2005, 20, 250–256. [Google Scholar] [CrossRef]
  19. Lungwitz, B.R.; Fredd, C.N.; Brady, M.E.; Miller, M.J.; Ali, S.A.; Hughes, K.N. Diversion and cleanup studies of viscoelastic surfactant-based self-diverting acid. SPE Prod. Opera. 2007, 22, 121–127. [Google Scholar] [CrossRef]
  20. Nasr-El-Din, H.A.; Al-Ghamdi, A.H.; Al-Qahtani, A.A.; Samuel, M.M. Impact of acid additives on the rheological properties of a viscoelastic surfactant and their influence on field application. SPE J. 2008, 13, 35–47. [Google Scholar] [CrossRef]
  21. Farajzadeh, R.; Andrianov, A.; Bruining, J.; Zitha, P.L.J. New Insights into Application of Foam for Acid Diversion. In Proceedings of the 8th European Formation Damage Conference, Scheveningen, The Netherlands, 27–29 May 2009. [Google Scholar]
  22. Nasr-El-Din, H.A.; Al-Anazi, H.A.; Mohamed, S.K. Stimulation of water-disposal wells using acid-in-diesel emulsions: Case histories. SPE Prod. Facil. 2000, 15, 176–182. [Google Scholar] [CrossRef]
  23. Al-Mutairi, S.H.; Nasr-El-Din, H.A.; Hill, A.D.; Al-Aamri, A. Effect of droplet size on the reaction kinetics of emulsified acid with calcite. SPE J. 2009, 14, 606–616. [Google Scholar] [CrossRef]
  24. Nguyen, Q.P.; Zitha, P.L.; Currie, P.K.; Rossen, W.R. CT study of liquid diversion with foam. SPE Prod. Opera. 2009, 24, 12–21. [Google Scholar] [CrossRef]
  25. Chang, F.; Qu, Q.; Frenier, W. A novel self-diverting-acid developed for matrix stimulation of carbonate reservoirs. In Proceedings of the SPE International Symposium on Oilfield Chemistry, Houston, TX, USA, 13–16 February 2001. [Google Scholar]
  26. Civan, F. Reservoir Formation Damage; Gulf Professional Publishing: Houston, TX, USA, 2015. [Google Scholar]
  27. Hassan, A.M.; Mahmoud, M.A.; Al-Majed, A.A.; Al-Nakhli, A.R.; Bataweel, M.A.; Elkatatny, S. Mitigation of Condensate Banking Using Thermochemical Treatment: Experimental and Analytical Study. Energies 2019, 12, 800. [Google Scholar] [CrossRef] [Green Version]
  28. Hassan, A.M.; Mahmoud, M.A.; Al-Majed, A.A.; Elkatatny, S.; Al-Nakhli, A.R.; Bataweel, M.A. Novel technique to eliminate gas condensation in gas condensate reservoirs using thermochemical fluids. Energy Fuels 2018, 32, 12843–12850. [Google Scholar] [CrossRef]
  29. Al-Nakhli, A.R.; Sukkar, L.A.; Arukhe, J.; Mulhem, A.; Mohannad, A.; Ayub, M.; Arifin, M. In-situ steam generation a new technology application for heavy oil production. In Proceedings of the SPE Heavy Oil Conference and Exhibition, Kuwait City, Kuwait, 6–8 December 2016. [Google Scholar]
  30. Panga, M.K.; Ziauddin, M.; Balakotaiah, V. Two-scale continuum model for simulation of wormholes in carbonate acidization. AIChE J. 2005, 51, 3231–3248. [Google Scholar] [CrossRef]
  31. Schwalbert, M.P.; Zhu, D.; Hill, A.D. Extension of an Empirical Wormhole Model for Carbonate Matrix Acidizing Through Two-Scale Continuum 3D Simulations. In Proceedings of the SPE Europec Featured at 79th EAGE Conference and Exhibition, Paris, France, 12–15 June 2017. [Google Scholar]
  32. Hill, A.D.; Rossen, W.R. Fluid placement and diversion in matrix acidizing. In Proceedings of the University of Tulsa Centennial Petroleum Engineering Symposium, Tulsa, OK, USA, 29–31 August 1994. [Google Scholar]
Sample Availability: Not Available.
Figure 1. CT-scan images for the rock samples after acidizing treatment without using thermochemical fluids (no acid diverter).
Figure 1. CT-scan images for the rock samples after acidizing treatment without using thermochemical fluids (no acid diverter).
Molecules 25 02976 g001
Figure 2. CT-scan images for the rock samples after diverting the injected acid using a thermochemical approach.
Figure 2. CT-scan images for the rock samples after diverting the injected acid using a thermochemical approach.
Molecules 25 02976 g002
Figure 3. Spatially correlated permeability ( I x = 3.0 ,   I y = 0.8 ,   σ k = 3.0 ) along the core dimensions before treatment.
Figure 3. Spatially correlated permeability ( I x = 3.0 ,   I y = 0.8 ,   σ k = 3.0 ) along the core dimensions before treatment.
Molecules 25 02976 g003
Figure 4. Spatially correlated porosity ( I x = 3.0 ,   I y = 0.8 , σ φ = 0.04 ) along the core dimensions before treatment.
Figure 4. Spatially correlated porosity ( I x = 3.0 ,   I y = 0.8 , σ φ = 0.04 ) along the core dimensions before treatment.
Molecules 25 02976 g004
Figure 5. Concentration profile at the breakthrough for the HCl acid treatment with no diversion, showing a wormhole growing mainly in the high-permeability core sample.
Figure 5. Concentration profile at the breakthrough for the HCl acid treatment with no diversion, showing a wormhole growing mainly in the high-permeability core sample.
Molecules 25 02976 g005
Figure 6. Concentration profile at the breakthrough for the thermochemical treatment case, showing a wormhole growing almost equally.
Figure 6. Concentration profile at the breakthrough for the thermochemical treatment case, showing a wormhole growing almost equally.
Molecules 25 02976 g006
Figure 7. Profiles of inlet (Pin) and outlet (Pout) pressures and the pressure-drop across the low-permeability rock during alternating injection of thermochemical and HCl fluids. The pressure drop is read from the y-axis in the right-hand side.
Figure 7. Profiles of inlet (Pin) and outlet (Pout) pressures and the pressure-drop across the low-permeability rock during alternating injection of thermochemical and HCl fluids. The pressure drop is read from the y-axis in the right-hand side.
Molecules 25 02976 g007
Figure 8. Profiles of inlet (Pin) and outlet (Pout) pressures and the pressure-drop across the high-permeability sample during the thermochemical/HCl treatment.
Figure 8. Profiles of inlet (Pin) and outlet (Pout) pressures and the pressure-drop across the high-permeability sample during the thermochemical/HCl treatment.
Molecules 25 02976 g008
Figure 9. Profiles of pressure-drop across the low-permeability (10 mD) and high-permeability (100 mD) rocks during the alternating injection of thermochemical and HCl fluids.
Figure 9. Profiles of pressure-drop across the low-permeability (10 mD) and high-permeability (100 mD) rocks during the alternating injection of thermochemical and HCl fluids.
Molecules 25 02976 g009
Figure 10. Images for the core inlet (left) and outlet (right), for the low-permeability sample after alternating treatment.
Figure 10. Images for the core inlet (left) and outlet (right), for the low-permeability sample after alternating treatment.
Molecules 25 02976 g010
Figure 11. Images for the core inlet (left) and outlet (right) for the high- permeability sample after thermochemical/HCl treatment.
Figure 11. Images for the core inlet (left) and outlet (right) for the high- permeability sample after thermochemical/HCl treatment.
Molecules 25 02976 g011
Figure 12. The injection rate ratio between top (Q1) and bottom (Q2) layer using thermochemical reaction and N2 injection.
Figure 12. The injection rate ratio between top (Q1) and bottom (Q2) layer using thermochemical reaction and N2 injection.
Molecules 25 02976 g012
Figure 13. The skin behavior during acidizing and diversion process using thermochemical reaction and N2 injection.
Figure 13. The skin behavior during acidizing and diversion process using thermochemical reaction and N2 injection.
Molecules 25 02976 g013
Figure 14. Skin reduction during alternating thermochemical and acid injection into high- (zone 1) and low-permeability layers (zone 2).
Figure 14. Skin reduction during alternating thermochemical and acid injection into high- (zone 1) and low-permeability layers (zone 2).
Molecules 25 02976 g014
Figure 15. Effect of pressure difference (between layers 1 and 2) on the acid diversion performance.
Figure 15. Effect of pressure difference (between layers 1 and 2) on the acid diversion performance.
Molecules 25 02976 g015
Figure 16. Effect of formation porosity on acid diversion performance.
Figure 16. Effect of formation porosity on acid diversion performance.
Molecules 25 02976 g016
Figure 17. Effect of layer permeability on the acid diversion performance.
Figure 17. Effect of layer permeability on the acid diversion performance.
Molecules 25 02976 g017
Figure 18. Acid diversion using thermochemical approach.
Figure 18. Acid diversion using thermochemical approach.
Molecules 25 02976 g018
Figure 19. A schematic of the different regions around the wellbore during acidizing treatment.
Figure 19. A schematic of the different regions around the wellbore during acidizing treatment.
Molecules 25 02976 g019
Figure 20. A schematic of the parallel coreflood setup.
Figure 20. A schematic of the parallel coreflood setup.
Molecules 25 02976 g020
Table 1. Reservoir, well, and fluid properties for liquid formation.
Table 1. Reservoir, well, and fluid properties for liquid formation.
ParameterValue
Porosity, ϕ (%)20
Viscosity of reservoir fluid, μf (cP)1
Viscosity of injected acid, μa (cP)1.5
Viscosity of thermochemical fluids, (cP)1.2
Wellbore pressure, pw (psi)3600
Pressure around wellbore, pe (psi)1600
Drainage radius, re (ft)1676
Wellbore radius, rw (ft)0.25
Table 2. Input data for liquid formation.
Table 2. Input data for liquid formation.
ParameterLayer 1Layer 2
Permeability, k (mD)10020
Thickness, h (ft)2020
Damaged radius, rd (in)1212
Skin, S012.512.5
Table 3. Model parameters for the sensitivity analysis scenarios.
Table 3. Model parameters for the sensitivity analysis scenarios.
ParameterBase CaseCase 1Case 2Case 3Case 4Case 5
Pressure Difference, DP (psi)20006003100200020002000
Porosity, ϕ (fraction)0.20.20.20.30.10.2
Permeability ratio, k1/k25555550
Table 4. Properties of the rock samples used in this study.
Table 4. Properties of the rock samples used in this study.
Sample IndexSample TypeDiameter (in)Length (in)Porosity (%)Permeability (mD)
Core1High-permeability rocks1.5122094
Core21.5122090
Core31.5419102
Core4Low-permeability
rocks
1.5122026
Core51.5122027
Core61.541811

Share and Cite

MDPI and ACS Style

Ba Alawi, M.; Hassan, A.; Aljawad, M.S.; Kamal, M.S.; Mahmoud, M.; Al-Nakhli, A. A Novel Approach to Improve Acid Diversion in Carbonate Rocks Using Thermochemical Fluids: Experimental and Numerical Study. Molecules 2020, 25, 2976. https://doi.org/10.3390/molecules25132976

AMA Style

Ba Alawi M, Hassan A, Aljawad MS, Kamal MS, Mahmoud M, Al-Nakhli A. A Novel Approach to Improve Acid Diversion in Carbonate Rocks Using Thermochemical Fluids: Experimental and Numerical Study. Molecules. 2020; 25(13):2976. https://doi.org/10.3390/molecules25132976

Chicago/Turabian Style

Ba Alawi, Mustafa, Amjed Hassan, Murtada Saleh Aljawad, Muhammad Shahzad Kamal, Mohamed Mahmoud, and Ayman Al-Nakhli. 2020. "A Novel Approach to Improve Acid Diversion in Carbonate Rocks Using Thermochemical Fluids: Experimental and Numerical Study" Molecules 25, no. 13: 2976. https://doi.org/10.3390/molecules25132976

APA Style

Ba Alawi, M., Hassan, A., Aljawad, M. S., Kamal, M. S., Mahmoud, M., & Al-Nakhli, A. (2020). A Novel Approach to Improve Acid Diversion in Carbonate Rocks Using Thermochemical Fluids: Experimental and Numerical Study. Molecules, 25(13), 2976. https://doi.org/10.3390/molecules25132976

Article Metrics

Back to TopTop