Effect of acid treatment on the geomechanical properties of rocks: an experimental investigation in Ahdeb oil field

Acidizing is one of the most used stimulation techniques in the petroleum industry. Several reports have been issued on the difficulties encountered during the stimulation operation of the Ahdeb oil field, particularly in the development of the Mishrif reservoir, including the following: (1) high injection pressures make it difficult to inject acid into the reservoir formation, and (2) only a few acid jobs have been effective in Ahdeb oil wells, while the bulk of the others has been unsuccessful. The significant failure rate of oil well stimulation in this deposit necessitates more investigations. Thus, we carried out this experimental study to systematically investigate the influence of acid treatment on the geomechanical properties of Mi4 formation of the Mishrif reservoir. The acid core-flood experiments were performed on seven core samples from the oil reservoir in central Iraq. The porosity, permeability, acoustic velocities, rock strength, and dynamic elastic parameters were computed before and after the acidizing treatment. To determine the optimal acid injection rate, different injection flow rates were used in the core-flooding experiments. The propagation of an acid-induced wormhole and its effect on the rock properties were analyzed and compared to that of intact rocks. Computed tomography (CT) scan and a 3D reconstruction technique were also conducted to establish the size and geometry of the generated wormhole. To analyze the influence of mineralogical variation and heterogeneity and confirm the consistency of the outcomes, acidizing experiments on different rock samples were conducted. The results demonstrate that for all the rock samples studied, the mechanical properties exhibit rock weakening post-acid treatment. The Young’s modulus reduced by 26% to 37%, while the Poisson’s ratio, the coefficient of lateral earth pressure at rest, and the material index increased by 13% to 20%, 23% to 32%, and 28% to 125%, respectively. The CT scan visually confirmed that the acid treatment effectively creates a pathway for fluid flow through the core.


Introduction
Acid treatment is an effective and efficient well-stimulation technique, and it has been applied extensively to stimulate carbonate reservoirs (Gomaa et al. 2018;Guo et al. 2014;Ituen et al. 2017;Zhu et al. 2015). In acid stimulation, fracture acidizing and matrix acidizing are two distinct subcategories. On one hand, matrix acidizing involves removing formation damage that impairs the permeability around the wellbore caused by drilling fluid invasion/fines migration; on the other hand, acid fracturing enhances the connectivity between the well and reservoir by creating pathways that penetrate deep into the reservoir referred to as wormholes (Al-Arji et al. 2021;Shafiq et al. 2018). In this regard, wormholing has been experimentally investigated and summarized by many authors (Al-Arji et al. 2021;Dong 2018;Xue et al. 2018). Al-Arji et al. (2021) and Xue et al. (2018) studied wormholing mechanisms where the rock minerals were dissolved by the contacting acid to form cylindrical pores during the wormholing process. Typically, acid fracturing yields better results in formations with low permeability and low confining stress or hard rocks. Conversely, acidizing is more effective in soft formations with high permeability and high confining stress . Various acid treatment techniques, such as acid jetting and CO 2 energized acid treatment, have been presented to improve acidizing effectiveness (Beckham et al. 2015;Ndonhong et al. 2016). Several researchers investigated the mechanism of acidrock reactions, acidizing fluid efficiency, acid flowback mechanism, acid leakoff, and the acidizing models Ghommem et al. 2015;Li et al. 2015;Lungwitz et al. 2007;Yoo et al. 2018;Zhang et al. 2020a, b, c;Zhu et al. 2015). On the other hand, others investigated the effect of acid treatment on the mechanical rock structure (Zhang et al. 2020b) and the influence of the rock mineralogy on acidizing efficiency (Martyushev et al. 2022). Various acidizing fluids such as self-diverting acid (Bazin et al. 1999;Lungwitz et al. 2007), visco-elastic surfactant (VES) acid, gelled acid, self-generated acids, and recently, emulsified acid (Martyushev and Vinogradov 2021) and chelating agents (Tariq et al. 2021) have been developed and investigated both in the laboratory and in field-scale (Gou et al. 2021;Hassan and Al-Hashim 2017;Isah et al. 2021a, b;Kiani et al. 2021;Lai et al. 2021;Li and Shi 2021;Martyushev and Vinogradov 2021;Melendez et al. 2007;Taylor and Nasr-El-Din 2001).
In the acid treatment of naturally fractured carbonate formation, VES acid forms more complex fractures compared to self-generated acid and gelled acid. However, gelled acid is capable of decreasing the rock's breakdown pressure significantly (up to 57% less than that achievable by water fracturing), thus increasing fracture propagation and enhancing efficiency (Gou et al. 2021). Moreover, gelled acid creates larger wormholes compared to crosslinked acid; consequently, it weakens the rock's mechanical properties more than crosslinked acid . Stimulation success depends on the length and width of these wormholes (Al-Arji et al. 2021); thus, a successful acid treatment operation requires that the wormhole propagates deep into the formation.
In radial acid treatment experiments of hollow chalk samples, Walle and Papamichos (2015) demonstrated that acidizing rock samples causes a reduction in the rock's mechanical strength, and this was confirmed by comparing the mechanical properties of the acid-treated rocks and the intact ones. Mustafa et al. (2022) studied the impact of acid wormholes on the mechanical properties of carbonates (chalk, limestone, and dolomite), and their findings indicate that acidizing reduces the hardness and elastic modulus of the rocks. The authors noted that dolomite was the least impacted by the acid treatment, while chalk samples were affected the most. Barri et al. (2016) investigated the effect of acidizing using chelating agents (ethylenediaminetetraacetic acid [EDTA] and diethylenetriaminepentaacetic acid [DTPA]) on the mechanical properties of carbonate rocks. The outcome showed that the elastic properties of weak carbonates (such as Austin chalk) were most affected, while hard rocks such as Indiana limestones were not significantly affected. Zhou et al. (2021) experimentally researched fracture surface strength before and after acid treatment. The authors argued that several reported mechanical deteriorations of carbonate rocks after acid etching could not be applied to evaluate fracture conductivity since such investigations provide information on the mechanical properties of the rock mass rather than the surface of the fracture. They emphasized that fracture surface strength measurement data before and after acid etching is necessary for fracture optimization and conductivity evaluation of the acidizing job. The work of Li and Shi (2021) also showed that acid fracturing can alter rock strength.
Certainly, the acid dissolution of rock minerals modifies the rock structure, the mineralogy, as well as the mechanical properties of the artificial fracture surface (Liu and Mostaghimi 2017). Therefore, acidizing often leads to modifications in the rock's mechanical properties around the wormhole; consequently, it improves or impairs reservoir quality. A notable body of literature reported rock weakening due to acidizing as reservoir impairment; however, rock loosening may imply the generation of flow paths by the acid dissolution of the rock and the consequent wormhole propagation into the reservoir. Several reports have been issued on the difficulties encountered during the stimulation operation of the Ahdeb oil field, particularly in the development of the Mishrif reservoir, including the following: (1) high injection pressures make it difficult to inject acid into the reservoir formation; and (2) only a few acid jobs have been effective in Ahdeb oil wells, while the bulk of the others has been unsuccessful. This called for more investigations. Thus, we carried out this experimental study to systematically investigate the influence of acid treatment on the geomechanical properties of Mi4 formation of the Mishrif reservoir in Iraq. The propagation of an acid-induced wormhole and its effect on the rock strength were analyzed and compared to that of intact rocks. CT scan was also conducted to establish the size and shape of the fractures generated. The findings of this investigation provide a better understanding of the acid treatment of carbonates; thus, comprehensive planning and implementation of acidizing jobs in the Mishrif reservoir can be accomplished.

Geological and petrophysical properties of Mishrif reservoir
The thin section photomicrographs of the cored sections of Mishrif reservoir's well ADM12 used for this study are shown in Fig. 1. A total of seven (7) rock samples were used in this investigation: samples 1-4 were extracted from section A (composed mainly of echinoderm, green algae, and micrite limestone), while samples 5-7 were cored from section B (micrite, orbitolina sand grains, and limestone). Both sets (sections A and B) were composed mainly of calcite mineral and 1-1.5% clay content, and the filling material was micrite. Table 1 shows the mineral composition of the experimental core samples obtained using X-ray diffraction (XRD) analysis. There are 34 zones interpreted in this well, and an aggregate of 135.8 m is accumulated with the overall thickness of the 19 oil and oil-poor zones. The entire thickness is 57.2 m, with seven transition zones. There are 8 water zones with an overall thickness of 150.8 m. A comprehensive description of the geological and fluid properties of this reservoir is provided in the appendices.
The wireline logging data of Mishrif reservoir's well ADM12 for petrophysical evaluation is provided in Fig. 2. It is necessary to distinguish between oil-rich and oil-poor zones. The oil-rich zones exhibit higher resistivity, the oil saturation of these zones is higher than 47 percent, and its effective porosity is greater than 15 percent. The shale volume of most oil zones is typically less than 10 percent. Overall, the Mishrif Fig. 1 Photomicrographs for the two cored sections of the Mishrif reservoir's well ADM12. a corresponds to samples 1-4 from section A: The texture is characterized by biological burrows, the internal ring is composed of echinoderm, and the core is composed of gastropods and green algae. The overall texture consists of echinoderm, green algae, micrite, large globigerina, and beehive worms. b cor-responds to samples 5-7 from section B: The limestone is predominated by orbitolina debris by strong solution, also named orbitolina sand grains. Visible dividual orbitolina and echinoderm with a large number of sand grains with a strong solution and dense gathering can be observed

Samples preparation
After coring, seven core samples were cut into the required dimensions (2.526 cm in diameter and a length of 4.571 cm each). To prepare the samples for strength tests and coreflooding experiments, the ends of the cores were ground to a smooth and parallel surface using end-face grinding. Solvents (toluene and alcohol) were used to clean the core plugs and prepare them for the measurements. Toluene was used to clean

Basic properties of gelled acid
Acid diversion is critical for stimulating vertical wells with extended target zones or horizontal wells in carbonates (Bazin et al. 1999). Increased viscosity of the injected acid and a delay in the acid interaction with the formation are two benefits of in situ gelled acids that improve treatment efficiency. Additionally, the gel should break down quickly as the acid is depleted, allowing for better clean-up when the acid treatment is completed (Taylor and Nasr-El-Din 2001). Table 3 shows the properties of the gelled acid. It is noteworthy that the viscosity of reacted acid is only 3 mPa s (around that of water); therefore, it can be concluded that there was no adverse effect such as formation damage by the acidizing fluid to the fracture and formation after the acidizing operation.

Core-flooding experiment
Core-flooding tests were conducted based on acid injection to establish wormholes in the plugs. Figure 3 shows a schematic of the core-flooding system. The core holder is connected downstream of the ISCO pump and the fluids (acid and brine) accumulator, and it is provided with pressure gauges to measure confining pressure and pressure drop across the core plug. The data acquisition system was used to monitor and collect data on the core-flooding process. The optimal injection rate was determined after seven acidizing experiments using varying flow rates. The flow rate ranged from 0.667 to 6.67 cc/ min. The core samples were vacuumed with a vacuum pump to ensure full saturation and then saturated with brine (see Table 4 for brine specifications) for 2 h; afterward, the samples were left in the desiccator for 48 h. An ENERPAC pump was used to provide a confining pressure of 2000 psi to the core plugs for all core-flooding experiments. To maintain a constant flow with minimal pressure variations, brine was injected first through the ISCO pump. The inlet and outlet pressures at various injection rates were recorded to determine the liquid permeability of the core samples using Darcy's law. Acid was then injected until the wormhole was created as indicated by a sharp pressure drop. Pressure and injection volumes were monitored and controlled using the core-flooding system's monitoring and control unit.

Porosity and permeability measurements
Pre-and post-acidizing core plug porosity and permeability measurements were performed. The gravimetric or saturation technique for calculating porosity was employed. Pore volume is calculated (Eq. 2), and porosity was then computed using Eq. 3: W dry is the weight of the dry plugs, W sat is the weight of the saturated cores, and W brine is the weight of the brine saturating the rock's pores. Permeability was assessed using flowing fluids with known viscosity through a core sample with known dimensions and then measuring the flow rate and pressure drop. In this work, we utilized brine flow through the plug and employed the Darcy law (Eq. 4) to compute liquid permeability: where Q is the instantaneous flow rate (m 3 s −1 ), K is the permeability (µm 2 ), is the dynamic viscosity of brine (Pa s), ΔP is the pressure drop across the rock plug (Pa), and A (1) 2.5 wt % Corrosion inhibitor 1.5% ADT-1 and 6% ADZ-1 (m 2 ) and l (m) are the cross-sectional area and length of the core, respectively.

Determination of rock geomechanical properties
Surface rock hardness, acoustic velocities, and dynamic elastic parameters were determined before and after the acidizing treatment. To analyze the elastic and geotechnical features of the carbonate rocks in this study, longitudinal (V p ) and transverse (V s ) wave velocities were measured for all samples before and after acidizing. The measurements were performed using the Sonic viewer (Model-5217A). A core plug connects the two transducers, one of which acts as a transmitter and the other as a receiver. Ultrasonic transducers are installed on both parallel sides to measure the core's transit time (T), which is determined by delivering a series of ultrasonic impulses through the core sample. Sample length (L) is divided by the transit time (T) of the propagated waves to obtain the velocity of the compressional and shear waves (V p and V s ). All core samples were analyzed for significant geotechnical characteristics based on the values of the velocities (V p and V s ). The combination of these findings was utilized to evaluate the subsurface rock's properties. Analysis of elastic and petrophysical characteristics before and after acid injection was carried out using the mechanical properties of the rock. Young's modulus (E) is one of the most important geomechanical properties of rocks. It is a measure of the rock's stiffness, that is, the material's resistance against being compressed (or extended) due to applied stress. The Young's modulus of the rock samples is computed using Eq. 5 before and after acid treatment using the  compressional and shear waves velocities (V p and V s ) (Fjaer et al. 2008). Poisson's ratio (PR) indicates the ratio of the lateral expansion (or contraction) deformation to longterm extension deformation. For the majority of mineral deposits, the Poisson's ratio is in the range of 0 to 0.5, with more rigid, competent, and incompressible rocks having a lower Poisson's ratio and vice versa (Abd El-Rahman et al. 1992;Domenico 1984), and it is computed using Eq. 6.
The elasticity of materials is impacted by several factors including the content of the material, the degree of consolidation, joints, fractures, and the presence of fluids in porous structures, all of which influence the material index (Shirani et al. 2010). The lateral earth pressure at rest coefficient (K o ) shows the degree of the strength and consolidation of rocks. It is in the range of 0 to 1 for most materials. As with Poisson's ratio, this coefficient decreases with densely consolidated sediments and rises with loose, unconsolidated sediments (Al-Awsi and Khorshid 2021). The material index (IM) and lateral earth pressure at rest coefficient (K o ) were computed using Eqs. 7 and 8, respectively.

Computed tomographic (CT) scanning
CT scanning was used to examine the size and geometry of the wormholes in each rock sample. Helical acquisition 140 kV and 500 mAs resolution CT images were taken to see the micro-scale alterations generated by the chemical interaction between acid solution and the core materials. The mechanical characteristics of various rock types are affected by the size of the wormholes created. The overall workflow for the experimental approach is shown in Fig. 4.

Results and discussion
The porosity, permeability, acoustic velocities, rock strength, and dynamic elastic parameters were computed before and after the acidizing treatment to establish the reasons behind the difficulties encountered during the acid stimulation operation of the Ahdeb oil field, particularly in the development of the Mishrif reservoir. Thus, we systematically analyzed the influence of acid treatment on the petrophysical and geomechanical properties of Mi4

Fig. 4
Schematic diagram of the workflow for the experimental procedure formation of the Mishrif reservoir. Acid efficiency curves were constructed using core-flooding tests on core samples (see Fig. 5). Ultimately, the goal is to discover the injection rate that results in the least quantity of acid being used throughout the process (Safari et al. 2014). Typically, the best injection rate results in the most optimal wormhole, reducing the amount of acid pore volume injection. In this regard, the reaction rate and convective mass transfer are the factors that regulate the wormholing process. Due to face dissolution, substantial amounts of acid are consumed before the wormhole can break through at low injection rates. Consequently, conical-shaped wormholes are formed. However, when the injection rate is substantially greater than the reaction rate, branching wormholes are formed. This was corroborated in the work of Mustafa et al. (2022). Thus, the development of a wormhole is maximized at a certain Damkohler number (i.e., the ratio of reaction rate to mass transfer). In this study, the rock samples yielded the lowest breakthrough injected pore volume of 3.7 PV at 2.16 cc/min; thus, this was considered the optimum injection rate. However, at 0.67 cc/min and 6.67 cc/min, the pore volume at breakthrough (PVBT) measured around 24 and 43.5 PV, respectively. The effects of PVBT on the elastic characteristics of the rock samples are detailed and analyzed in the subsequent section. The acidizing treatment resulted in a considerable increase in porosity for each of the rock samples. A positive relationship between PVBT and the relative increase in porosity for all rock plugs was observed as can be seen in Fig. 6. The greater the PVBT was, the more likely the relative porosity increased. This finding can be attributed to increased rock dissolution/etching as more PV is injected.

Effect of acid treatment on geomechanical properties
The geomechanical properties including the acoustic parameters (V s and V p ) of the rock samples were determined after acidizing treatment. The following subsections provide a comparative analysis of the pre-and post-acid-treated rocks.

Impact of acid treatment on acoustic properties
Various factors affect the seismic velocities of rocks, among which are density, lithology, mineralogy, porosity, grain size, stress levels, joints and fractures, pore fluid, anisotropy, and temperature. These parameters (i.e., the seismic velocities) are used to compute the rock geomechanical properties such as the elasticity modulus or Young's modulus (E), Poisson's ratio (PR), and rock material index (IM) of rocks (Dobrin 1976).
The compressional (P) and shear (S) waves of all plugs were monitored before and after the acid experiment to confirm the alteration of the ultrasonic waves due to acid treatment. A comparison between the ultrasonic waveforms of acid-treated samples of rocks and those of intact saturated samples of rock was performed. Ultrasonic waveforms for both intact and treated cores are shown in Fig. 7, illustrating how wormholes and fluid flow alter velocity, amplitude, and frequency for plug sample No.1. The acid-treated sample's P-wave arrival time is longer compared to that of the intact rock, and the S-wave amplitude is more significant compared to the intact rock's relatively low amplitudes. P-wave signals from untreated and treated rock have somewhat differing frequencies.

Impact of porosity and wormhole on the elastic characteristics of rock
An inverse correlation between effective porosity and compressional velocity (V p ) was established in this work, as However, rocks with a high porosity before acid treatment but a persistent network of wormholes after acid treatment have a propagation velocity that is lower compared to that of rocks predominated by pores as compared to micropores. In general, the higher the porosity is, the lower the propagation velocity is. Figure 9 illustrates the relationship between porosity and shear wave. It is consistent with the trend in compressional velocity before acid treatment; however, the relationship between porosity and shear wave (V s ) is less obvious after acid treatment. This is due to the rock structural modification by acid dissolution. As the porosity increases, V s demonstrates a lower decline as compared to the slope before acid flooding. Based on the established positive relationship between PVBT and the relative increase in porosity of the rock plugs, it can be concluded that an increase in PVBT will likely result in an increase in porosity and a consequent decrease in the values of the acoustic properties of the rock.

Rock mechanical properties pre-and post-acid treatment
This section compares the findings of the dynamic elastic properties, that is, Young's modulus, Poisson's ratio, lateral earth pressure at rest coefficient, and material index obtained before and after acidizing. Formations can be damaged by drilling and production activities, and factors such as the migration of fines and mud invasions may alter the permeability around the wellbore (Shafiq et al. 2018). Completion activities may cause a change in pressure and flow rate in a wellbore, and this is known as a mechanical skin factor (Furui et al. 2008;Yildiz 2006). Damage around the wellbore formation might lead to further pressure drops in the surrounding region, reducing reservoir output. As a result, the acidizing process is used to minimize the damage and restore the reservoir's productivity by enabling the hydrocarbons to migrate easily (Ghommem et al. 2015). Highpermeability flow pathways are created, and the damage can be removed to increase the well's performance. However, during acidizing operations, acid dissolution of rock minerals can modify the rock's structure, the mineralogy, and the mechanical properties of the artificially fractured rocks (Liu and Mostaghimi 2017). Therefore, the effects of acid treatment on Young's modulus (E) of the seven rock samples studied are shown in Fig. 10. The Young's modulus (E) values before acid flooding are in the range (1.37 E + 10 -2.77 E + 10), while after acidizing, the E values are between (8.62 E + 09 -2.04 E + 10). This implies that the pre-acid treatment values of E showed a general reduction compared to the values after acidizing. As a measure of rock stiffness, materials with high resistance to being compressed (or extended) due to applied stress display larger values of E, while softer, less stiff rocks have lower E values. Thus, the results are indicative of the rocks' weakening post-acid flooding. Since this outcome is for the rock mass, that is, the whole core sample and not only the fractured surface, it implies that a wormhole is created, which causes an increase in void space in the rock and consequently lower seismic velocities; thus, Young's modulus reduces.
The Poisson's ratio (PR) values of the pre-and post-acidtreated rock samples are compared, and the results are shown in Fig. 11. Specifically, the findings revealed that the PR In geotechnical engineering, the coefficient of lateral earth pressure at rest (K o ) is crucial since it influences the strength and consolidation of rocks, as well as the design of engineered structures such as optimal wellbore direction to prevent wellbore instability issues. The knowledge of K o can be used to optimize stimulation jobs to avoid sanding after acidizing. From the K o values obtained, a similar pattern is noted: it decreases for densely consolidated sediments and rises in loosely consolidated strata. Rock samples pre-and post-acid treatment were examined for changes in the values of the coefficient of lateral earth pressure at rest, as illustrated in Fig. 12. A range of 0.43 to 0.59 was found for the K o values of the rock samples before acid flooding, corresponding to rocks that were moderately consolidated before flooding. The K o values for the acid-treated plug samples ranged from 0.57 to 0.73. Fractures or wormhole channels are indicated by rocks with higher values, while rocks with lower values exhibit moderately confined pores. Material index (IM) values before and after acid treatment are shown in Fig. 13. In this regard, the parameter's values for the rock samples pre-acid treatment range between -0.2 to -0.49, indicating the hardness of the rocks. On the other hand, the values of this parameter range from -0.45 to -0.63 for the acid-treated rock samples, indicating rock weakening as a result of rock mineral dissolution during the acid reaction. A summary of the mechanical parameters of sample 1 and sample 5 before and after acid treatment is provided in Table 5.

Computed tomography (CT)
In addition to the petrophysical and acoustic measurements, CT scanning was employed to visualize the results of the acid core flooding. Nondestructive 3D imaging was conducted using computed tomography. When X-rays travel through a substance, they are attenuated, and the corresponding material image is formed (Tembely et al. 2021). CT scan has been employed to monitor fracture openings (Martyushev and Yurikov 2021). In this study, acid core-flooding tests were carried out at various injection rates on different rock samples, causing the wormholes to differ in structure from one rock sample to the other and even within the same rock type. The size and nature of the wormhole may be responsible for the change in mechanical and physical characteristics of the rock (Artyushkova et al. 2021;Martyushev et al. 2019;Mustafa et al. 2022;Song et al. 2021;Zhang et al. 2020a). Micro-CT images of two of the acidized rock samples were reconstructed using the 3D Slicer software. A medical CT scanner was used to scan the chosen plug samples before and after acid treatment. Contrast and noise reduction were prioritized while scanning. Figures 14, 15 and 16 show the CT scan of plug sample No. 1 after acid treatment with a flow rate injection of 6.67 cc/min (Figs. 14b and  16a), where the wormhole channel was created as a result of convection, thus forming a wormhole, and plug sample number No. 5 after acid treatment with flow rate injection of 0.667 cc/min (Figs. 15b and 16b), where the acid reaction is the dominant mechanism, thus creating a conicalshaped wormhole (Fredd and Fogler 1999).

Recommendation
In acid flooding, rock weakening does not necessarily mean deterioration/sand production/instability; instead, it depends on other factors. For instance, it could be due to acid-induced fracture/wormhole, which corresponds to porosity enhancement that resulted in lower acoustic readings and thus manifested as reduced rock strength. Accordingly, the following recommendations are drawn: • We suggest that independent investigations of the bulk rock and fractured surface strength be conducted to be able to make a clear distinction between these occurrences. • Furthermore, the experimental results presented in this study for the MI4 formation of the Mishrif reservoir demonstrated that despite the observed rock's mechanical properties alterations, which show wormhole propagation, effective wormholes/flow paths were created for hydrocarbon flow from the reservoir into the wellbore.
The results presented implied a successful acid treatment of the Ahdeb oilfield despite the development problems of the Ahdeb oilfield reported (i.e., the high injection pressures that makes acid injection difficult in this reservoir formation and the numerous acidizing job failures) during the stimulation operation of this reservoir formation. • Formation permeability greatly governs the choice of reservoir acid stimulation technique to be employed, that is, either acid fracture or matrix acidizing. Generally, acid fracturing yield better results in low permeable hard rocks. On the other hand, matrix acidizing is more effective in soft, high permeability formations . Based on the core and well logs analysis for all of Mishrif's units, permeability values vary from 0.1 mD to 57 mD (see Fig. 2). This shows a broad range of permeability variations; thus, this petrophysical property should be prioritized when planning stimulation operations for the Ahdeb oilfield, particularly in the development of the Mishrif reservoir. • Therefore, the high injection pressures reportedly encountered are likely due to injections into zones with very low permeability, such as 0.1 mD zones, thus resulting in the injection difficulties. Moreover, the acidizing jobs failures could also be as a result of wrong choice of acid treatment technique, and it is thus recommended that acid fracturing should be applied for low permeability and hard rock zones, while matrix acidizing should be considered for zones with high permeability in the Mishrif reservoir. Therefore, stimulation design optimization is recommended to reduce these failures as well as injectivity issues.

Summary and conclusions
Significant challenges are encountered during acid stimulation operations in the Mishrif reservoir of the Ahdeb oil field, including high injection pressures, and several acid treatment failures are recorded. The significant failure rate of oil well stimulation in this deposit necessitates more research. Therefore, we presented experimental investigations of the effect of acidizing on the geomechanical properties of the Mi4 formation of the Mishrif reservoir. Several acid-flooding experiments were performed on different rock samples to study the influence of mineralogy and confirm the consistency of the outcomes. The porosity, permeability, acoustic velocities, rock strength, and dynamic elastic parameters were computed before and after the acidizing treatment. The propagation of an acid-induced wormhole and its effect on the rock strength was analyzed and compared to that of intact rocks. Thus, the following conclusions are drawn: • We found that the acid treatment resulted in a considerable increase in porosity for all the rock samples. • A positive relationship between PVBT and the relative increase in porosity for all rock plugs was established. Thus, we concluded that an increase in PVBT will likely results in an increase in porosity and a consequent decrease in the values of the acoustic properties of the rock. • The acid efficiency curve yielded the lowest pore volume injected at breakthrough of 3.7 PV at 2.16 cc/min; thus, the optimum injection rate that results in the most optimal wormhole and the least quantity of acid being used for this reservoir is 2.16 cc/min. • Furthermore, the results also demonstrated that the mechanical properties exhibit rock weakening post-acid treatment. The mechanical properties (i.e., elastic modulus, Poisson's ratio, material index, and the coefficient of lateral earth pressure at rest) showed that the rock weakened after acidizing, and this finding is in line with the work of other researchers Mustafa et al. 2022;Zhang et al. 2020a). • Young's modulus for the seven rock samples investigated exhibited a reduction of 26-37%. Conversely, the Poisson's ratio, the coefficient of lateral earth pressure at rest, and the material index increases by 13-20%, 23-32%, and 28-125%, respectively. • The CT scan visually confirmed that the acid treatment effectively creates a pathway for oil flow from the reservoir to the wellbore. • It further showed that at low acid flowrate (0.67 cc/ min), face dissolution occurred, and a substantial amount of acid was consumed before the wormhole could breakthrough at this rate, thus resulting in shorter, conical-shaped wormholes. At high acid injection rate (6.67 cc/min), a longer and wider wormhole was generated.

Appendix
Comprehensive geological properties of the Mishrif reservoir are provided in Fig. 17. The geology and mineralogy of Mishrif reservoir's well ADM12, Cenomanian Mi 4 formation is shown. Photomicrographs for the two cored sections of the Mishrif reservoir's well ADM12. (a) For samples 1-4: a texture with burrow texture, internal ring is composed of echinoderm. Overall texture is echinoderm, green algae micrite texture, large globigerina, and beehive worm; (b) for sample 5-7: Limestone composed mainly of orbitolina debris, also termed orbitolina sand grain, visible dividual orbitolina and echinoderm with large amount of sand grain with and dense gathering. Table 6. shows the petrophysical and fluids properties of the Mishrif reservoir's Mi 4 formation. The porosity and the effective permeability are 22.2% and 10.32 × 10 −3 m 2 , respectively. The residual oil and connate water saturations are 31.23% and 35.09%, respectively.
Funding There is no funding for this work.

Conflict of interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.  Effective permeability of oil,10 −3 m 2 10.32