Oilfield chemical-formation interaction and the effects on petrophysical properties: a review

Oil and gas recovery may cause formation damage during drilling, completion, and production phases. As a result of fundamental chemical, thermal, mechanical, and biological interactions, formation damage can occur due to impairment of permeability and porosity, causing undesirable operational and economic problem. The fluid-rock interactions resulting from oilfield chemicals injection during drilling, enhanced oil recovery (EOR) such as chemical flooding, or formation treatment could negatively impact on the formation properties such as geomechanical and geochemical, leading to alteration of the rock’s petrophysical properties. These chemical-rock interactions induce changes in both pore space geometry and rock strength. The resultant impact includes weakening of the formation bonding materials, formation damage, reduced production and consequently sand production simultaneously with reservoir fluids. It is therefore critical to evaluate these variables prior to designing any geo-sequestration, reservoir stimulation or EOR projects. Studies have shown that rock properties, especially permeability, porosity and strength, are altered or damaged during drilling, cementing, perforating, producing, stimulating, and injecting water or chemicals for EOR. Petroleum companies are likely to suffer significant financial losses due to this. This study provides a review on the influence of oilfield chemical-formation interactions on the formation rock properties both geophysical and mechanical, leading to formation damage and sand production. This study aims to provide researchers with a single document that gives insight and new perspectives on oilfield chemical-rock interactions through compilation of recent studies relating to the effect of chemical-rock interactions on rock's petrophysical properties, as well as geomechanical properties due to geochemical reactions that cause formation damage and eventually sand production. Having a solid understanding of fluid-rock interactions and how they impact petrophysical properties and cause formation damage is essential in predicting sand production and would help in minimizing economic losses, downtime and technicalities.


Introduction
Petroleum has since the mid-1950s become a major source of energy to the world (UKOG 2015). As the world's primary source of energy, it plays a major and influential role in global economy (Burclaff 2020). Petroleum, which is a naturally occurring hydrocarbon, can be found buried deep in sedimentary rocks (Priscilla 2021). To access the hydrocarbon (mainly petroleum), the sedimentary rock otherwise known as a reservoir rock has to be drilled (Alyafei 2019). However, in drilling the reservoir rock, stress is exerted on the formation, which may cause it to fail once stress exerted on the reservoir rock surpasses the reservoir rock strength (Oyeneyin et al. 2005). In addition to drilling and completion, damage to formations can result from reservoir rockfluid interactions during reservoir stimulation or enhanced oil recovery (EOR) operations. Hence, predicting petrophysical properties of reservoir formation rocks is important for petroleum exploration, development, and production (Yang et al. 2017a(Yang et al. , 2017b(Yang et al. , 2017c. It is therefore obvious that the mechanism of reservoir rock failure can be categorized as either chemically, mechanically, or operationally induced (Abass et al. 2002;Zhou and Sun 2016). Figure 1 shows a graphical process for fluid injection into an injection well for the purpose of enhancing oil production. Fluid injected into the reservoir interacts with the formation rock and impacts on the geochemical and geomehcanical properties of the rock resulting in the weakening of grain particles and consequently the production of sand while producing reservoir fluids.
This review focuses on formation damage caused by oilfield chemical-rock interaction and their impacts on formation properties and the potential for sanding. Fluidrock interactions resulting from carbon dioxide (CO 2 ) geosequestration, EOR fluids, and reservoir stimulation projects are impacted by the nature of formation rocks, the fluid type and properties, and the evolution of petrophysical properties of the formation under in situ conditions. The chemicals used for these oilfield operations include biocides (Wuyep et al. 2018), surfactants, scale inhibitors (Tomson et al. 2008;Oluyemi 2014aOluyemi , 2014b, scavengers, demulsifiers, and defoamers. The mobility and distribution of these chemicals in the reservoir depend on relative permeability and the porosity of the formation rock. Mineral dissolution and precipitation processes, as well as fine particle migration caused by fluidrock interactions, are known to cause rock porosity, permeability, and strength alterations. Furthermore, oil production from a reservoir is strongly influenced by the permeability near the wellbore; therefore, formation damage caused by fluid-rock interactions may reduce productivity significantly (Faegestad 2016). Reservoirs formation experiences numerous fluid-rock interactions; thus, it is crucial to identify the specific geochemical reactions and processes governing the evolution of porosity and permeability (Yang et al. 2017a(Yang et al. , 2017b(Yang et al. , 2017c. Investigations have shown that various chemicals would have different properties, which would have varying influences on the strength of the reservoir formation when they are injected (Wuyep et al. 2020). Oluyemi 2014a, 2014b) experimentally investigated the sanding potential of using scale inhibitors (PTEMP) on a clashach reservoir rock under dynamic condition. It was discovered that detached sand grains were found in the flowing stream. This suggests that fluid-rock interactions can potentially weaken formation bond resulting in sand production with the reservoir fluids. More recently, Wuyep et al. (2018) studied the physical and chemical impacts of oilfield chemical and rock interface on reservoir rock strength using biocide to control and treat microorganisms such as sulfate reducing bacteria (SRB). The results of the work showed that chemical dissolution/precipitation reaction occurred, and rock fine particles were transported resulting in the weakening of the formation strength. SRB is widely dispersed in the oilfield waters and their activities is the cause of reservoir souring. The activities of microorganisms in oilfield reservoirs and sulfide accumulation due to bacterial sulfate reduction can cause serious damage such as H 2 S production and formation plugging (Larsen 2002). Not only does H 2 S cause damage to the reservoir, it also presents life-threatening danger to oil industry workers when exposed it. Its corrosiveness also has a major negative effect on subsea infrastructures, thereby increasing maintenance cost and increasing environmental disaster in relation to oil leakages and spills (Hubert 2017). The injection of biocides for oilwell treatment as used by Wuyep et al. (2018) is the most common method for reservoir treatment (Gieg et al. 2011). But these chemicals are unsafe to both the environment and oilfield workers. However, another approach that is used in controlling the actions of sulfate reducing bacteria is the application of nitrates, which is able to display diverse influence on the formation (Davidova et al. 2001). Most research on nitrate usage has focused on oil recovery (Hitzman and Sperl 1994;Sandbeck and Hitzman 1995;Giangiacomo and Dennis 1997; Zhao et al. 2016) and reservoir souring control and treatment (Reinsel et al. 1996;Telang et al. 1998;Sunde et al. 2004;Greene et al. 2006;Grigoryan et al. 2008;Hubert 2010;Voordou et al. 2011;Tabari et al. 2011).
There are three kinds of chemical damage mechanisms: fluid-rock interactions, fluid-fluid interactions, and wettability alteration (Faegestad 2016). Therefore, another major concern regarding formation damage is wettability alteration. In addition to determining the capillary pressure, wettability also determines whether the pressure is a barrier to fluid flow or a driving force for production (Deglint 2018). Hydrocarbon recovery from rock matrix by imbibition has been greatly facilitated with the use of surfactants as wettability modifiers (Yang et al. 2021). Consequently, surfactants are frequently used in carbonate reservoirs to convert oilwetted to water-wetted wettability to enhance spontaneous imbibition. In this scenario, the pore geometry and wettability-related fluid−formation matrix interactions, as a result of mechanisms for wettability alteration, can induce changes in porosity, pore connectivity, permeability, and flow. Pal et al. (2018) published a detailed review on surfactant-based chemical EOR techniques applied in carbonate formations. Although chemical-EOR is a typical chemical-rock interaction, the focus of the review was on surfactant-driven oil recovery mechanisms, like altering wettability, reducing interfacial tension, understanding microemulsion phase behavior, analyzing surfactant adsorption and mitigation, and using foams. It is important to note that oilfield chemicals injected into the formation is absorbed and upon interaction with the rock, poses a weakening influence formation rock strength (Wuyep et al. 2019). Injecting oilfield chemicals into reservoir rocks is to improve oil and gas production, clean formations that are blocked, reengineer formation damage, improve oil flow assurance, and prevent scale formation and corrosion can be harmful to the rocks petrophysical and geomechanical properties. Investigation by Oluyemi (2014aOluyemi ( , 2014b revealed that the geomechanical and geophysical impact of oilfield chemicals on the rock properties could lead to sand production, which has negative cost implication in the petroleum industry. In the petroleum industry, sanding negatively impacts flow lines and separators as a result of rock and fines plugging perforated zones (Abass et al. 2002;Willson et al. 2002). This undesirable effect also increase cost of operations (Rahmati et al. 2012;Wang et al. 2016), credited to pipeline corrosion and other surface facilities, a decrease in production rate and expensive intervention measures (Tronvoll and Fjaer 1994;Rahmati et al. 2012 ). The corrosion of well integrity and well equipment's have an increasing and negative implications which includes wellbore failure and a rise in disposal costs and downtime (Penberthy and Jr Shaughnessy 1992;Eshiet 2012 ). The implication of producing sand is also unfavorable to the environment. For instance, the tear and wear of surface and downhole equipment, as a result of the extremely corrosive nature of fluids carrying sand, may be sufficiently adverse; it could lead to total collapse of surface and downhole facilities, resulting to severe health, safety, and environmental problems (Penberthy and Jr Shaughnessy 1992).
Geo-sequestration of CO 2 , injected fluid for EOR, reservoir acidizing, and reservoir stimulation projects demand an analytical assessment of rock properties, fluid properties, and the evolution of rock petrophysical properties as a result of fluid-formation interactions. For instance, the injected and in situ-produced solutes will interact chemically and physically with the storage formation rock during many underground carbon dioxide (CO 2 ) geologic carbon sequestration. The extent of the chemical-formation interactions depends on the concentration of the propensity minerals such as carbonates and clays in the pores of the formation rock (Saeedi et al. 2016). It has been found that the injection of CO 2 into a brine saturated reservoir resulted in carbonic acid formation, which in turn, react with the minerals constituting the rock's porous framework, resulting in mineral dissolution and/or precipitation (Saeedi et al. 2016;Zhao and Yin 2021). These chemical reactions change the petrophysical properties of the formation rock. Several studies have been conducted on the influence of fluid and rock interfaces on reservoir formation properties. Yang et al. (2017aYang et al. ( , 2017bYang et al. ( , 2017c investigated the transformation of reservoir formation properties of sandstone reservoir core samples, and they discovered that the alterations of the properties of the formation was due to the dissolution of feldspar, calcite, and other cement bonding minerals assembly present within the mineral content of the sandstone when reacting with brine and ethanoic acid. In a series of studies reported by Kan et al. (2005), Tomson et al. (2008), Bybee (2010), Egermann and Vizika (2011), Mohamed and Nasr-El-Din (2013), Oluyemi (2014aOluyemi ( , 2014b, Wilson (2016), and Wuyep et al. (2019), they found that the application of oilfield chemicals resulted in changes in the formation rocks geomechanical properties and geophysical properties, respectively. On the other hand, high-temperature operations, such as steam injection and in situ combustion to mobilize oil for production, can lead to the development of thermal damage mechanisms. Minerals can also dissolve and transform at elevated temperatures when they are catalyzed and transformed from nonreactive clays to reactive products, which can swell, merge, and reduce permeability (Faegestad 2016). Also, the introduction of bacteria and nutrients into the formation can induce biological damage. Thermal gradient and pH influences dissolution/precipitation of minerals which modify the rock's porosity and permeability (Fritz et al. 2010). The mineralogy, pore geometry, thermal maturation, and wettability of the formation are factors that are thought to control fluid flow in unconventional reservoirs (Wang et al. 2021). One of the established methods for studying fluid-rock interactions is core flooding, which is a powerful approach for elucidating fluid-formation interactions. Thereafter, imaging techniques such as scanning electron microscopy (SEM) and X-ray diffraction can be used to quantify fluid-rock interaction and characterize petrophysical properties due to dissolution/precipitation, fracturing, and fine particles migration (Deglint 2018). One of the major challenges is the scale-dependence of petrophysical and geomechanical properties of rocks. The imaging of the micro-and nano-scale should be the way forward in order to detect fluid-rock interactions at the micro-and nano-pore level, so that the distribution and flow of multi-phase fluids can be quantified at this scale. This is particularly important since the dominant pore and pore throat sizes in reservoir rocks are in the meso-and micropore range. However, the prediction of the possible evolution of petrophysics of rocks would involve geochemical modeling of fluid-rock interactions resulting in dissolution and precipitation of minerals, which would influence the mineralogy and petrography of the rocks (Fritz et al. 2010). Another approach would be a combination of experimental and analytical techniques to evaluate impacts of fluid-rock interactions on formation petrophysical properties.
Therefore, evaluating the thermos-chemo-mechanical processes resulting from fluid-formation interactions and their impact on the petrophysical properties of the rock during production is of critical importance (Faegestad 2016). This reactive transport through oil-formation causes rock microstructural properties evolution, such as porosity and permeability as a result of the flow and chemical reactions. This information can be used to formulate chemical and simulation studies, rock evolution models, and formation failure predictions. It is crucial to understand the evolution of formation petrophysical properties because the changes affect reservoir quality. An integrated model combining kinetic rate laws for minerals dissolution/precipitation, particle migration, and fluid flow is needed to predict the geochemical evolution of the fluid interactions and also sanding. It is therefore necessary to fully understand the performance of the reservoir formation and the tendency to produce sand or rock failure upon interaction with oilfield chemicals, in order to avoid or minimize losses resulting from negative impacts of oilfield chemicals on the petrophysical and geomechanical properties of the formation. Therefore, this paper reviews previous studies on oilfield chemicals and reservoir rocks interaction and their impact on reservoir petrophysical and geomechanical properties, analyzing and evaluating findings relevant to developing new perspectives for further investigation. In this review, we examine the effects of chemical injection on petrophysical and geomechanical properties of rock, geochemical reactions such as mineral dissolution and precipitation, and interactions between chemicals and rock, which result in formation damage and sanding. It is the aim of this review to provide insight into oilfield chemical-rock interactions, and to offer new perspectives on its impact on petrophysical and geomechanical properties of rocks by critical evaluating recent studies on the subject in order to guide future research.

Effects of oilfield chemical on petrophysical properties
It is well known that injected fluids interact with rock minerals and/or formation fluids and generate fine particles in addition to dissolution/precipitation which can cause damage to the formation, either from fluids in the reservoir or generated in situ as result of fluid-rock interactions. These processes induce properties changes to the rock permeability, porosity, and strength, which are considered to be very important factors for reservoir rock petrophysical properties (Ahmad et al. 2018), indicating that it can be adversely affected by geochemical and geomechanical phenomenon. Crawford (2011) developed a model which relates petrophysical and petrographic properties to rock strength. This suggests it is plausible to predict and estimate mechanical properties from geophysical properties of the rocks due to fluid-rock interactions. As a consequence, formation damage and onset of sanding can be predicted also. Using rock properties to calculate onset collapse of a pore can be helpful approximation to predicting onset collapse of the pore and changes in compressibility of the reservoir formation with the reduction in pressure of the fluid at uniaxial strain boundary settings. A review of selected approaches and models that have been developed for sanding prediction has been published (Rahmati et al. 2013). The majority of these models are built on a continuum assumption, but a few have recently been developed on the basis of discrete element models. The analytical models are commonly used for estimation of the onset of sanding, while numerical models are suitable for predicting sanding rate. Overall, a realistic sanding model should include all failure mechanisms (shear, tensile, compression), as well as fluid flow effects. Discrete element method (DEM) is another valuable tool to simulate sand production specially to understand the mechanism of sanding, but it cannot be used for large-scale problems because of large computational time required and difficulty in the calibration of model (Rahmati et al. 2013). In the investigation as described by Hill et al. (2018), data on the formations geomechanical and petrophysical properties like the rock mineralogy and grain size distribution were combined to explain the impact of fluid-rock interaction, the formation stability and the prospect of sand failure and its production. The pressure gradient developed in the pore channels due to fluid flow also facilitates the detachment of sand grains. Additionally, the flow of fluid leads to the transport and production of detached sand aggregates, as well as disaggregated sand grains. Therefore, sand management and modeling are always cost-effective if they are implemented early before a well is drilled. This complex phenomenon depends on various factors such as the stress distribution, the properties of the rock, and fluid-rock interactions in the reservoir (Rahmati et al. 2013). The fluid flow is mostly modeled using the Navier-Stokes equations solved by a finite volume method. However, it can be challenging to account for all these factors and mechanisms in the numerical models, resulting in limitations on their applicability.

Pore volume, porosity, and permeability
Petroleum reservoirs which are made of sedimentary rocks usually have porosity with ranges of 10-40% in sandstone formations and between 5 and 25% in carbonate formations. Sand production is a common occurrence in weakly consolidated sandstone reservoirs, which host a substantial portion of the world's oil and gas reserves (Rahmati et al. 2013). The strength-weakening effect of fluid-rock interactions may gradually lead to sandstone reservoir formation degradation. Additionally, rock porosity and permeability evolution influence chemical distribution and mobility in the reservoir. Porosity, as shown in Fig. 2, determines the rocks storing capacity which is described as the ratio of its void spaces, generally termed pore volume, to bulk volume and is described either as a percentage or a fraction (Coneybeare 1967;Keelan 1982).
The study of fluid-rock interaction requires an understanding of transport and geochemical kinetics. Presently, reactive transport is faced with the most challenging problem of how to treat cases where the petrophysical properties of the porous reservoir formation change due to geochemical reactions (Steefel and Mäher 2009). The most frequently identified evolution in petrophysical properties are porosity and permeability which are induced by chemical reactions such as dissolution or precipitation. The modeling approach involves coupled flow and chemical reaction. In reactive transport simulations, porosity is typically the first-order parameter predicted, as it is directly related to the total volume of minerals precipitated or dissolved (Steefel and Mäher 2009). Since the detailed pore geometry is required, it is more difficult to predict transport parameters, in particular permeability and diffusivity. It is possible to circumvent this problem by combining modeling and microscopic characterization of the porous reservoir formation as mineral dissolution or precipitation occurs. The governing advectiondispersion-reaction model is shown in the equation below (equation 1).
The aqueous (homogeneous) reactions, R r , and the rates of the N m solid-phase (mineral) reactions, R m , u is the fluid velocity, D is diffusivity, and stoichiometry given in the coefficients ν ir and ν im .
Currently, reactive transport in porous media is described by three types of models: (1) continuum models, (2) porescale models, and (3) multiple continuum or hybrid models involving a combination of scales (Steefel and Mäher 2009). Among the transport processes relevant for the study of chemical-rock interaction are advection, which occurs most commonly due to fluid flow through porous media, molecular diffusion, which is modified by electrochemical migration in the presence of charged species, and hydrodynamic dispersion.
The phenomena of oilfield chemical-rock interaction can actually cause alterations to the formation's porosity and permeability rock (Xu et al. 2018;Lamy-Chappuis et al. 2014). Therefore, this modifies the rock permeability and porosity characteristics. Also, a scientific model that describes the evolution mechanism of porosity and permeability triggered by several chemical and physical interactions among reservoir rocks and fluids has been investigated by Chang and Civan (1991). There is however not one fit all relationship among permeability and porosity applied to all porous medium. This can be credited to the added complication introduced by precipitation and the release of fine grains which is likely to plug pore throats, swelling of clay and improved compaction (Chang and Civan 1991). Thus, the link between permeability, porosity, and rock strength has been a focus of wide study, soon to be (1) Fig. 2 A graphic view of a porous media (Liou 2005) completely understood. Conversely, sedimentary rocks such as carbonate reservoirs are known from studies for their huge commercial hydrocarbon accumulation around the world with altered petrophysical properties, mainly due to precipitation and dissolution of minerals leading to increased porosity. The dissolution/precipitation of chemical specie that leads to changes in porosity depends on parameters like mineral composition, shape of grain, grain size distribution, and pore size (Wuyep et al. 2018). Additionally, Egermann et al. (2010) studied the petrophysical properties as a result of fluid (CO 2 injection) and rock interface on the capacity to hold and store carbon dioxide. The rock-CO 2 interaction was observed to alter the petrophysical properties such as permeability of the carbonate rock. This occurred as a result of the pore network framework imbalance and the degree of mineral dissolution in the rock. In light of this, many permeability and porosity experimental equations have been reported in literatures. Porosity is primarily associated with tortuosity and a large number of small pores, so dissolution of small sand particles or distinct calcite may have only a small impact on the connectivity of pores, as well as a decrease in tortuosity and/or dissolution at pore throats. As a result, the maximum strength of siltstones, sandstones, and a few limestones will reduce with increased porosity ( To study rock's petrophysical properties, the Kozeny-Carman (K-C) equation is usually fitted with experimental data to define the relation between porosity (ϕ) and permeability (κ) in a fluid-rock interface system (Walsh and Brace 1994; Lamy-Chappuis et al. 2014). This advocates powerlaw relationships among the appropriate pore geometry, such as tortuosity and surface area, hydraulic radius, porosity, and permeability. But Lamy-Chappuis and co-workers (Lamy-Chappuis et al. 2014) validated their experimental data with the use of equation (1), whereas equation 3 as established by Chang and Civan (1991) can be used to evaluate permeability alteration due to formation damage. Hence, equation 4 is the general form of equations (1) and (3).
In the equations, ϕ denotes porosity, k permeability (ϕ o and κ o represent initial porosity and permeability respectively), n ranges between 3 to 7, the efficiency factor f stated as the unplugged pore throat area fraction (determined stochastically) and the experimental matching parameters are represented as a, b, and c.
Furthermore, Yang et al. (2017aYang et al. ( , 2017bYang et al. ( , 2017c investigated sandstone rock and brine interface, examine rock properties during diagenesis, and the results revealed that minerals such as calcite dissolution leads to porosity increase. Using empirical, numerical simulations, petrographic observations, and the ability to restructure the diagenetic process, the authors concluded that reservoir formation petrophysical properties can be predicted. Additionally, Liu et al. (2019) examined the effect shale and water interface on the sensitivity of stress on shale rock rich in organic materials. Shale reservoirs were found to have natural fractures and low permeability based on the study results. However, when the fluid used is water, the chemistry between water and rock is rather different from that between oilfield chemicals and the rock interface, which is the focus of this study. On the other hand, the interaction between iron carbonate and corrosion inhibitors showed that imidazoline can bind to, adsorb to, and develop a film on non-conductive materials such as iron carbonate (FeCO3), resulting in a reduction in porosity and delaying the transportation of reactants to the corrosion surfaces has reported in the literature. This suggests that porosity is a function of the materials strength. In addition, Wuyep et al. (2019) found that a decrease in permeability and porosity suggests a decrease in pore space due to new minerals formation in the carbonate rock when treated with Betaine and Aminotri methylphosphonic acid (ATMP). It is therefore obvious that there is a strong connection between uniaxial strength and the porosity of formation materials, which is a function of the composition and internal structures (Li and Aubertin (2003). In other words, porosity is an important parameter when dealing with uniaxial strength. In a study by Bertaux and Lemanczyk (1987), the behavior of sandstone reservoir rock treated with potassium hydroxide and sodium hydroxide (KOH and NaOH) was examined. It was discovered that the interaction between the sandstone rock and alkali fluid at 90°C and 150°C, produced a new mineral crystalline zeolite and precipitated amorphous at 200°C. The amount of precipitated and dissolved minerals was found to be responsible for changes in permeability and clay reactivity with alkali chemical.
A typical example where the impact of oilfield chemical-rock interaction is experience is in chemical flooding techniques, such as alkali-surfactant-polymer (ASP) flooding for EOR (Pal et al. 2018). A study on the incremental oil recovery factor and degree of permeability damage in heterogeneous sandstone reservoirs subjected to strong-base and weak-base ASP flooding processes were also evaluated and compared (Zhong et al. 2020). It was found that the permeability-damage ratio can be reduced by approximately 15% with the weak-base ASP flooding compared to the strong-base flooding, and the reservoir flow assurance issues related to chemical loss can be addressed. This shows that the formulation of oilfield chemicals and their concentration. As reservoir rock minerals dissolve, the grains are carried away by fluid flow, which settle on the surface of the pore walls. The induced porosity change during fluid flow and rock-fluid interface is as a result of the precipitation of insoluble salts, formation swelling, clay particles deposition, mineral dissolution, and surface solid entrainment (Chang and Civan 1991). Shebl and Surdam (1995) empirically studied the possibility of oilwater-rock reacting. The authors utilized oxidized sandstone rock which originally contains anhydrite and intergranular clay cements with 6-15% range porosity and about 10-25% carbonate. After the experiment, it was found that there was a major possibility for a redox reaction to occur between an oxidized mineral and crude, and this reaction impacted permeability and porosity (increasing it from a range of 12-20%) of an elastic reservoir formation. Based on the results from literature, it can be concluded that chemical and rock interfaces cause mineral precipitation/dissolution, changing the porosity and permeability of the formation. But, changes in the permeability and porosity of the rock depends on the rock's mineralogy, particle size distribution, shape of the grain, and pore size (Lamy-Chappuis et al. 2014;Wuyep et al. 2019).

Effect of chemical-rock interaction on particle size distribution
When a reactive fluid flows through a geologic porous medium, such as a reservoir rock, it is possible for minerals embedded in the reservoir formation to dissolve, releasing fine particles. While the dissolved mineral species could be in ionic form, the precipitated counterparts and the sand grains comprise a range of particle size distribution. More also, the minerals (ions) may as well precipitate from other ions dissolved in the fluid or they may also be conveyed with the dissolved minerals to form precipitates when the concentration of the ions is greater than the solubility of the precipitates (Chang and Civan 1991). The solid particle size distribution in solution (see Fig. 3) may influence the ionic exchange between ions in the solution and rock surface and the possibility to precipitate. This is because several materials make up the reservoir rock, which differ in their porosity, micro-voids, and boundary structure, and thus are characterized by a fragile structure (Jongpradist et al. 2011). Therefore, the redistribution of grains and their arrangement in the reservoir formation is closely related to the characteristics of the rock (Kwan and Surip 2020). It can increase the deformation of rocks, affect their porosity and permeability, and lead to structural uncertainty (Alonso 2010).
The reactions that occur when oilfield chemicals are introduced into a geological formation are classified as follows: (1) mineral dissolution and precipitation of the rock, (2) precipitation and dissolution of metallic ions and ionic compounds in the chemical, (3) ion exchange between the rock surface and the chemical, and (4) crystal growth and nucleation of the precipitate (Chang and Civan 1991). Dvorkin and Gutrierrez proposed that particle size distribution (PSD) is an important feature of rocks, as it directly affects the rocks geomechanical properties and geochemical characteristics. It is possible therefore for chemical-reservoir rock interaction to cause distortions in mineral configuration, causing mineral dissolution and precipitation and, in turn, decreasing or increasing rocks petrophysical properties depending on particle size distributions (Lamy-Chappuis et al. 2014). But when interparticle pores dominate fluid flow channels, such  (Oluyemi 2014a(Oluyemi , 2014b as sandstone, they are moderately consistent in their spread. Therefore, it is possible for all pores of the bonding material to be completely sealed off by its range of particle sizes distribution (He and Stephens 2011). According to Oluyemi et al. (2006), changes in the PSD of a sand producing formation during production may change the uniaxial compressive strength (UCS) of the rock. This changes UCS can weaken the rocks resulting formation damage. Consequently, Israeli and Emmanuel (2018) with the use of MATLAB simulated the effect of PSD on formation weathering. It was observed that rocks with fine grains rapidly withered more but slows down when the grain size is increased. This reveals that particle size distribution strongly influences the rate of fluid-rock interactions. Even so, the study ignored physical processes such as fluid flow, diffusion, and the impact of oilfield chemicals and rock interface, making extrapolation difficult. On the other hand, Wuyep et al. (2019) reported the impact of oilfield chemical and rock interface on the rocks mechanical strength with a combined use of uniaxial compressive and analytical test. The results showed that there is a major difference between particle size of brine and its effluent from reservoir samples. The increase in PSD observed may be caused by the dissolution/precipitation reactions, which can weaken the grain-grain bonding and result in detachment. However, results from other reported researches showed that rock strength can be linked with the composition of its minerals (Li et al. 2014;Pan et al. 2016;Sun et al. 2017;Zhang et al. 2018), grain shape and size (Sun et al. 2017;Jin et al. 2017), and its particle size distribution (Hussain et al. 2006;Shimizu et al. 2010 ). According to Jin et al. (2017 ), the shapes of rock materials and grain size contribute to the strength of its cemented soil-rock mixtures. It is found that the cementcemented soil-rock mixtures could have localized shear-and strain-softening bands with a better strength and modulus than soil-rock mixtures without cement. In other words, the cement materials help to improved strength and modulus of formation soil-rock mixture. In general, the PSD is critical in evaluating the response of the reservoir rock formation to fluid interaction (Lade et al. 1996;Yu et al. 2016). By using DEM code to investigate the effect of PSD on the mechanical properties of the reservoir, it was found that the mechanical properties are influenced by the porosity of the formation and PSD (Shimizu et al. 2010 ). As the small grains of the rock occupies the space between the larger grains, they are restricted from dislodging by opposing grains. Furthermore, Wu et al. (2018 ) investigated the effects of grain size distribution on the structural and mechanical properties of cemented rockfills based on loading, the results obtained showed that cemented rockfill with rough particles displays a decreasing microstructure with an early failure. This is attributed to the effect of PSD composed within the rock. In essence, particle size distribution plays a very important role on rock deformation (Wang 1994;Zhang et al. 2018), as well as changes in permeability and porosity (Seely and Johns 2018).

Geomechanical properties
There is the possibility of sand production when there is not enough strength in the formation to withstand the destabilizing forces generated by reservoir fluid flow (Gharagheizi et al. 2017). Therefore, predicting sanding conditions in advance is necessary to provide technical support for sand control strategies. Therefore, a complete knowledge of rocks' mechanical properties such as elastic properties and the inelastic properties (e.g., fracture gradient and compressive strengths) are crucial for the successful design of oil recovery and a safe oilwell drilling and stimulation of petroleum reservoirs. In order to predict distortion and failure of reservoir rock, reservoir fluid flow through porous formations and fluid-rock interfaces must be modeled. The general behavior of the rock as a result of external load or stress majorly depends on the rocks geomechanical properties. It is also dependent on other parameters which include non-hydraulic stress, temperature, pressure, and the formation failure history. The compressive strengths, Young's modulus, and elastic modulus have been considered as important parameters for relating rocks' mechanical behavior under stress. Makani noted a strong connection between elastic modulus and compressive strength of the rock. Similarly, it has been reported that the rocks petrophysical properties, geochemical properties, and composition exact much greater influence on the mechanical behavior. Therefore, the mechanical and physicochemical properties of the rock controls its behavior in addition to the interaction between the reactive chemical and its flow. Therefore, rock texture characterization is needed to provide insight on microstructure and morphology. Engebretson et al. investigated the effect of polyethylene oxide (PEO), dodecyltrimethyl ammonium bromide (DTAB), and aliminium chloride (AICl 3 ) on sandstone reservoir rock strength with the aim of establishing fundamental knowledge that can be used to optimize chemically assisted fracking. This study was however limited in scope to static tests, may not adequately predict the dynamic processes occurring during reservoir stimulation or production. Jegarluei and coworkers (2010) reported the impact of cementation on the reservoir rocks mechanical properties where it was noted that the rocks mechanical properties are significant parameters for modeling wellbore stability and sand production.
In terms of understanding the in situ loads or stresses experienced in reservoir rocks, geomechanics can provide insight on the rock's stability and strength. By using the petrophysical and mechanical properties of the formation, we can predict the stability of the rock and its ability to withstand the changes in stress as a result of fluid flow. Oluyemi (2014aOluyemi ( , 2014b empirically investigated a clashach reservoir core to evaluate the interaction that takes place between a clastic reservoir formation and a chemical inhibitor. The results of this study showed that the interaction between a chemical inhibitor (phosphate scale inhibitor) and clastic reservoir rock could result in a physical failure of the rock and the release of sand during reservoir fluid production. Although numerical analysis was not fully taken into account in this study, it is recommended that further empirical studies be conducted in order to confirm the physicochemical models. Oilfield chemicals have been shown to have a large impact on the formation's geomechanical properties when injected into the formation (Wuyep et al. 2018 and. In 2018, Giani et al. investigated the effects of gas production-induced compaction and subsidence analyses, coupling fluid-flow and geomechanics (i.e., stress-strain phenomena) by combining three numerical models: (1) geological model of the rock structure and its geologic features; (2) geomechanical model that predicts the rocks stress-strain behavior; and (3) fluid flow model that simulates the evolution of pressure in space and time. As part of the study, an experimental scale that represents an offshore, weakly consolidated, shallow, multilayered gas formation was constructed, utilizing already prevailing data from offshore gas bearing formations. It was found that the evolution of the stress-strain on the rock remained within its elastic limits. The changes were small, the compaction-induced porosity reduction was negligible, and the maximum permeability variation reached 12%. The reason for the minimal impact on the rock geomechanical and petrophysical properties could be ascribed to the gas fluid and the weak compaction induced by the flow field. However, mineral precipitation and dissolution processes may alter permeability, flow pathways, and therefore deform a formation by fluid-rock interaction. The flow of fluids in porous rock formation occurs either through interconnected pores or through a network of natural pre-existing or induced rock fractures. Obviously, rock permeability, porosity, and geomechanical properties are altered by reactive fluid-rock interaction. Therefore, there need to be considerations for the coupling of flow, reaction, and rock deformation in modeling the processes. Furthermore, the chemicals could change the consistency of the stressed rock by mineralizing or eroding contact surfaces, resulting in chemically motivated deformation. This is due to changes in volume triggered by processes, such as oxidation, hydration, precipitation, and dissolution, carbonation of minerals in the rock that is restricted under stress. However, geochemical and mechanical disturbances to reservoir formation can occur independently, the processes are often in combination with fluid and mineral interactions (i.e., fluid-rock interaction). Though geochemical reactions may start and promulgate crack, change fracture opening because of mineral dissolution, pore closing due to mineral precipitation, lead to volume change and rock deformation. At the same time, it could affect the geomechanical processes. Hence, models that combine chemical-mechanical process in stressed rocks are critical and valuable for predicting and monitoring dynamic changes to the state of the rock and deformation evolution. Finally, the petrophysical, geomechanical properties, and morphology of the formation develop vigorously because of mineral precipitation and dissolution caused by pressure variations, temperature, fluid motion, and composition.
The compressive strength refers to the rock's ability to withstand stress intended to reduce its size. The compressive strength is a very significant geomechanical property of the formation mostly used to determine the stability of the structure when exposed to stress (Romana and Vásárhelyi 2007;Gholami and Fakhari 2017). The test determination approach involving core samples of rocks is destructive, costly, and time-consuming. To estimate UCS, various models have been proposed using mineralogical and petrographic characteristics, and rock's physical properties. In light of this, the strength of a rock is exponentially determined by its porosity, and this increases with a decrease in porosity, as demonstrated by Sarda et al. (1993). Thus, it is important not only for engineering designs, but also an important element of an effective design and implementation (Aladejare 2020). A good knowledge and precise evaluation would improve the proficiency of oilfield stimulations and other oilfield operations (Negara et al. 2017). Rashidi et al. (2008) in their effort to enhance bit real-time wear with the use of an intelligent drilling advisory system was able to demonstrate that uniaxial compression strength (UCS) test may be affected by particle shape, grain size, sorting, minerology, size and shape of test sample, and rate of loading. On the other hand, Sonmez et al. (2006) noted that compressive strength which is an input in rock design is essential in the prediction of the deformation modulus. Nevertheless, Palmstrom (1996), Aladejare and Wang (2019) utilized rock mass classifications system which combines rock mass index (RMI) and rock mass rating (RMR) system in evaluating the formation strength properties and rock mass deformation connected to response to load and its design. Tixier et al. (1973) proposed a simple correlation to estimate formation strength from its mechanical property, especially UCS and rock strength; where ρb and ∆t can be obtained through sonic and density logs.
Using a triaxial test, Yadav et al. (2016) investigated the changes in geomechanical properties of Berea shale and sandstone reservoir rock samples when they were interacting with oil and water based muds. It was shown that the oilbased-mud (OBM) preserves the shale strength better than the water-based-mud (WBM). Xu et al. (2018) empirically (5) b∕Δt 2 studied the influence of drilling fluid on tight sandstone hardness using WBM (water-based mud) and OBM (oil-based mud), and observed that after 2 h, the hardness of the rock (sandstone) decreased by 22.9% and 10.1% with WBM and OBMs respectively. It also shown that after 15 days, the sandstone hardness with WBM was observed to have decreased by 33.1% but remained unchanged with OBM. It was also noted that temperature has but only little influence on sandstone hardness with WBM but at temperatures of 50°C and above decreases the sandstone hardness with OBM. This means water exerts more impact on the geomechanical properties of the rock formation than oil. According to Wuyep et al. (2019), fluid-rock interaction resulting from mineral dissolution, precipitation, and geochemical reactions weakens reservoir rock grain fabrics of carbonate and sandstone. This decreased the UCS under atmospheric and ambient conditions of the experiment. Kitamura and Hirose (2017) investigated the impact of distillated water on the strength of several reservoir rocks (sandstone) with the use of indentation test in studying hardness of the different sandstone rock samples. The Young's modulus and compressive strength were found to increase with a decrease in porosity of the rock when the compressive strength was correlated with the ultrasonic wave velocity. Experimentally, Al-Osta (2018) evaluated the compressive strength of rocks using a correlation developed between the compressive strength and some important parameters gotten through indirect methods on the rock such as ultrasonic pulse velocity, mass density, and point load strength index. The obtained results showed that the connection between ultrasonic pulse velocity and UCS compressive strength displayed the best match, followed by compressive strength and the mass density, and lastly point load index and compressive strength. Likewise, Muqtadir et al. (2018) explored how the mechanical properties of a sandstone formation with low permeability is affected by fluid saturation. The authors observed that the sandstone formation when saturated in brine becomes weaker compared to when it is saturated in oil. The investigation also showed a 9% reduction in the UCS and 40%d tensile strength (TS) decrease when saturated with brine. But when the formation is saturated with oil, the reduction in the TS and UCS were discovered to be 25% and 10% respectively. This suggests that oilfield chemicals exert significant impact on the geomechanical properties of the reservoir rocks. Similarly, Smorodinov et al. (1970) explored correlations between some physical characteristics of rocks and rock strength, developing two relationships to estimate UCS based on porosity and density as a direct determinant of rock strength which are exponential functions of density and porosity. This indicates the possibility of estimating geomechanical properties of the rocks from their petrophysical properties. For instance, it has been generally shown that the rock's compressive strength is a function of mass density or porosity. Conversely, other correlations are based on the Poisson's ratio, the rock strength, sonic wave velocity, and bulk density (Fjaer et al. 1992;Moos et al. 1999).
It is therefore expected that the UCS decrease as the rock porosity increases; hence, at a certain porosity, the rock will become weakened and merely a loose aggregate. The confining pressure (stress) below which a reservoir formation deforms is acknowledged as its confined compressive strength (CCS), at which point, a breakdown in the rock gain structure begins, leading to a compaction of the rock and a loss in porosity. The CCS of a formation is a very significant parameter in petroleum engineering for selecting drill bits and drilling process optimization, predicting the rate of penetration (ROP) and performance analysis (Shi et al. 2015). It is also general understanding that the heterogeneity at rock grain scale, such as mineral stiffness and grain geometry, as well as particle and particle boundary micro-cracks (called microdefects) will reduce the UCS of the rock (Bahrani and Kaiser 2016). In considering rock strength, Rampersad et al. (1994) used a power function to develop a link between the CCS of the rock, strength of the rock, and confining pressure. Nonetheless, the model proposed neglected the impact of porosity on the strength of the rock, which is a critical petrophysical property. Shi et al. (2015) considered the effect of porosity and nonlinear properties with increasing confining pressure for selecting bits, predict rate of penetration (ROP), drilling optimization, and developed a compressive strength model. The authors' model when in combination with equation for mechanical specific energy was used to optimize drilling parameters. Figure 4 shows the relationship between confining pressure and rock strength for different types of rock. It can be observed that the model can recognize unproductive drilling circumstances and is appropriate for both overbalanced drilling (OBD) and underbalanced drilling (UBD). However, the study was carried out on a dry rock sample and therefore makes the result difficult for extrapolation in accounting for the rock behavior interacting with fluids.
The rock porosity is a function of peak stress at varying confining pressure as presented in Fig. 5. It is clear the rock porosity plays a significant role in both the elastic parameters of the rock and the strength of the rock. As a result, the load-bearing capability of the formation changes with change in the rock porosity, of which it has been reported that oilfield chemicals interaction with reservoir formation induces significant alterations to petrophysical properties of the rock, and by extension to the geomechanical properties, such as rock strength and UCS.

Chemical-rock interaction mechanisms
The main mechanisms when fluid and rock interact are mineral phase precipitation and dissolution (William et al. 2016). It is therefore important to understand the mechanisms occurring at the chemical-rock interfaces in order to determine what impact they have on the formation strength.
To develop the governing equations, momentum and mass conservation, continuity, reaction rate laws, and other constitutive relationships will be considered. Thus, the formationdamaged model for the rock during chemical-rock interaction should be the combination of fine particles migration, chemical reactions (dissolution and precipitation), and fluid flow models. Wuyep et al. (2019) evaluated biocides, scale, and corrosion inhibitors effects on carbonate and sandstone rocks. It was found that ionic exchange, mineral precipitation, and mineral dissolution contributed to the weakening of the rock grain fabrics, leading to rock failure. Guodong et al. (2017) examined the influence of CO 2 -water interface on permeability and porosity. It was observed that increasing CO 2 injection without proportionate pressure build-up influences the permeability-porosity relationship. It was however resolved that an increase in CO 2 injectivity was accountable for the increased rock mineral dissolution (dolomite). Bertaux and Lemanczy reported the interaction mechanism between potassium hydroxide (KOH), sodium hydroxide (NaOH), and the minerals of a sandstone reservoir formation. The results from these investigations confirmed that the primary mechanisms due to fluid-rock interactions are precipitation and dissolution. Horsrud et al. (1998) carried out a study to predict borehole instability by exposing shale sample rich in smectite to KCL. It was discovered that the shrinkage in the shale sample was as a result of cation exchange. The effect of the shrinkage was simulated, and it was observed that increasing the concentration of KCL resulted in a decreased compressive stress. This suggests that the chemical interaction exerted a significant effect on the rock's strength due to ion precipitation. In another investigation by El-Hajj et al. (2013), a carbonate reservoir sample was flooded with carbon dioxide at reservoir condition and assessed the CO 2 -carbonate formation interaction with an X-ray fluorescence (XRF) analyzer, scanning electron microscope (SEM), and an atomic force microscope (AFM). The authors noted that at high-pressure aging, CO 2 dissolved some of the rock minerals and also precipitated some of the dissolved carbonate. Equally, Shiraki and Dunn (2000) empirically studied the effect of brine and saturated CO 2 on dolomite-anhydrite-cemented sandstone rock. A week after, the authors discovered that the dissolution of dolomite and K-feldspar led to the formation of kaolinite in the pore throat and consequently reduced the permeability. Muller et al. (2009) hydrochemicaly experimented on Berea sandstone rock by using dry CO 2 to flush a saturated NaCl core, and discovered a reduction in permeability which was as a result of halite mineral precipitation. On the other hand, a calcite-dolomite cemented sandstone was reported by Ross et al. (1981), it was observed that the increase in permeability of the rock was due to carbonate mineral dissolution. This was credited to an increased pore space alteration due to the carbonate cement bond dissolution. In summary, dissolution and precipitation are the most common mechanism induced by chemical-rock interactions, which remarkably impact of the rock's petrophysical and geomechanical properties.

Geochemical reactions
Chemical-rock interface induces numerous chemical and physical processes as the oilfield chemicals flow through the porous rock system. The oilfield chemicals may dissolve  (Shi et al. 2015) some solid minerals embedded in the rock (Samuel and Thomas 2018). In addition to ion precipitation in the fluid, downstream transport of precipitates and dissolved minerals can also occur if the concentration exceeds the solubility of the precipitates (Chang and Civan 1991). In most cases, however, only precipitation and dissolution are considered as the geochemical reactions occurring. Investigation on the effects of mineral dissolution on reservoir rock petrophysical properties have been earlier defined by Colón et al. (2004) and Hu et al. (2010). The authors in their conclusion reported that permeability and porosity will increase upon these chemical interactions.

Chemical reaction (dissolution-precipitation)
Mineral precipitation and dissolution involve a concurrent process of solute transport in addition to adsorption and desorption of surface solutes, mineral surfaces and bulk solution, migration of solute on surface and dehydration, and hydration of ions (Appello and Postma 2005). Nonetheless, some of these processes may proceed faster than others. Consequently, clay particles residing in the rock pore spaces are thought to be responsible for the release of fines during fluid-formation interactions. But changes in pH and salinity could result in clay swelling and clay deflocculation. According to hypothesis, fines migration reduces permeability by causing formation damage (Saeedi et al. 2016). As a result of mineral precipitation and fine particles migrating within the host rock pores, formation permeability decreases (Zhao and Yin 2021). There are several types of mineral fines that contribute to permeability reduction, including clay, quartz, feldspar, and carbonate as they can fill or bridge pores. It is well known that rock permeability plays a significant role in fluid flow. On the other hand, the dissolution of rock minerals would result in an increase in formation porosity. The dissolution of formation may lead to a collapse of pore network structures. In addition to the changes in pore space geometry, the fluid-rock interactions also modify rock strength as stated earlier. Nonetheless, the changes in porosity and permeability of the underground formation rock-fluid interaction influences the injection and subsurface movement of the chemicals. Due to fluid-rock interaction, the petrophysical properties of rock may be altered, resulting in lower rates of production than expected. Therefore, the source of the reduction must be identified before any corrective measures can be taken.
The reactive fluid flow momentum through a porous material is accountable for the heterogeneous geochemical reactions like precipitation and dissolution that occur at the rock minerals and the pore-fluids interface (Lamy- Chappuis et al. 2014;Wuyep et al. 2020). In other words, fluid flow can be described as a reactive transport system. Thus, the model development for the chemical-rock interactions comprise precipitation, dissolution, and the transportation of minerals dissolved. As oilfield chemicals flow through a porous medium, it is adsorbed on the rock surface and due to the mineralogical components, it distorts the equilibrium system. As a result, mineral dissolution occurs resulting to a fluid configuration that is saturated with added minerals (Putnis 2015). Although, the rate of rock dissolution is a function of the chemical composition, temperature, ionic strength, chemical concentration, fluid motion, pH, and the hydrodynamic interface between the oilfield chemical and the rock phase that depends on the surface area of the mineral-fluid interface (Lasaga et al. 1994;Lasaga 1998). A rock's surface roughness is also a significant factor determining its reactivity, with rough surfaces dissolving more rapidly than smooth surfaces (Levenson and Emmanuel 2013).
The dissolution and precipitation reactions that occur as a result of chemical-rock interaction when fluid flows through a porous medium can also be modeled by the discretization of rate laws and chemical law of mass action. Liu et al. (1997) characterized the reactions in two groups: mineral and aqueous reactions. While the aqueous reaction only entails aqueous species such as acid/base dissociation, the mineral counterpart includes both mineral and aqueous species like precipitation and dissolution. The rate of mineral reactions is quick when compared aqueous reactions (Liu et al. 1997). Furthermore, the rate of reactions in reactive chemical-rock surfaces can also be influenced by the concentration of catalysts such as precipitates/dissolves at the surface, or the quantity of microorganisms acting on the rock surface. Saaltink et al. (1998) observed that the concentration of species in the medium contributes to the rate of reaction. Therefore, the rate of reaction is related to change in specie concentration. For instance, considering compound A to be converted to B (Scheme 1), the rate of reaction, r, may be modeled as expressed in Eq. 8. Equation 9 shows the equilibrium constant for the two reacting compounds A and B (scheme 2). The chemical reaction rate (r) describes the quantity of reactants that evolves to products per unit time.
where K eq represents the equilibrium constant, k the rate constant, n is reaction order. For simple cases, n is equals 0, 1, 2, for zero, first and second order reactions, respectively. Injecting reactive fluids in a porous medium will lead to the dissolution of several minerals that originally exists and can produce new species, dependent on temperature and chemical composition, the new chemical species may precipitate. The mineral dissolution and precipitation due to geochemical reaction with the reservoir formation is complicated and involves several mechanism rates law. Hence, a rate law would model the combined impact of fluid flow, chemical reaction such as overall rates of precipitation and dissolution of the minerals and factors like ionic strength, temperature, saturation conditions, and pH (Lasaga et al. 1994). Notably, significant changes in permeability and porosity are experienced in the rock due to mineral dissolution and precipitation, thereby modifying fluid flow path features. This establishes a link between the fluid flow and chemistry taking place in the porous medium. To simulate these coupled processes occurring concurrently, different numerical software tools have been developed. TOUGHRE-ACT is a numerical simulation tool developed and deployed to solve the reaction kinetics of precipitation and dissolution resulting from chemical--rock interaction kinetics in fluid flow through a porous medium, and also to study reactive chemical transport in a porous medium and geochemistry. The overall reaction rate R i is a combination of heterogeneous mineral reactions which is precipitation-dissolution reaction represented as (R i min ) and homogeneous aqueous reaction (R i aq ) as given in Eq. 10. The heterogeneous reaction (R i min ) represents the amount of all the fluid and rock reactions as shown in Eq. 11, whereas the rate of species production j from the aqueous reaction (R i aq ) phase as a result of homogeneous reaction is given in Eq. 12.
where r m signifies the rate of dissolution (R m less than 0) or precipitation (R m greater than 0) of minerals m per unit volume of formation, v im the mole number of i present in the mineral m, N m the amount of minerals that exists in the formation.
Mineral precipitation or dissolution from fluid-rock interface occurs simultaneously and through a series of developments like solutes transport between the rock surface and bulk fluid, surface detachment (desorption), and attachment (adsorption) of solutes, surface migration, dehydration, and hydration of ions. Steefel and Lasaga and Lasaga et al. (8) 1994) established a rate law for kinetically controlled mineral precipitation and dissolution reactions which is applied to reactive flow in a porous medium. In this study, the reaction rate is basically a function of all species concentration as expressed in Eq. 13. The overall rate of dissolution or precipitation of rock mineral in aqueous solution is given in Eq. 14.
where f (a i ) signifies the function of the activities of different ions in solution, A represents the reactive surface of the mineral (m 2 /m 3 total rock), the rate constant k (mol/m 2 s) and f (∆G) denotes the free energy of solution.
Based on computational viewpoint, the reaction rate law is governed by fluid flow with surface-controlled reactions like dissolution or precipitation as restrictive steps are much easier to determine. Appello and Postma (2005) also suggested a generic rate law, similar to that shown in equation (12) that takes into consideration changes in solute concentration due to mineral dissolution or precipitation reactions (Eq. 15). Notably, the rate is dependent on the ratio of the surface area A 0 to the chemical solution volume V and specific rate.
Given, R represents the general reaction rate (mol/L.s), specific rate k (mol/m 2 .s), volume of chemical fluid V (m 3 ), initial surface area of the rock A 0 (m 2 ), the original moles of solid m 0 and the moles of solid at a specific time m. From equation (14), the term (m/m 0 ) n defines change in the surface area, especially the reactive surface sites throughout mineral precipitation/dissolution. Nonetheless, for uniform size of grains precipitating or dissolving cubes and spheres, n = 2/3. For different particles size population, the finer grains will dissolve selectively. Thus, n is a function of the original PSD with a value ≤3.4 (Appello and Postma 2005). The function g(C) considers the effect of the chemical composition on rate, pH, the effects of catalysis and inhibition and the distance from equilibrium (Lasaga 1998). Similarly, Steefel and Lasaga presented another simplified expression that simulates the kinetics of precipitation and dissolution of reservoir rock minerals due to oilfield chemical and reservoir formation interface is shown in Eq. 16. Similarly, Rimstidt and Barnes and Saaltink and co-workers (1998) also reported a similar rate equation. However, since most rate constants are obtained at ambient temperature, the Arrhenius equation as presented in Eq. 17 can be used to study the effect of temperature on reaction rate constant k.
where −1 and + 1 are dissolution and precipitation respectively, k in mol/m 2 s denotes kinetic rate constant, σ specific reactive surface medium of rock mineral, ψ solution state with regards to the rock mineral (ψ = Q/K eq , K eq , and Q are equilibrium constant and ionic activity product of the formation, and θ and η are empirical parameters determine experimentally, activation energy E a and the rate of constant k25 given at 25°C, temperature T in degree Kelvin and R1 is gas constant. Note, when dissolution occurs, ψ < 1, > 1 for precipitation, the system will be at equilibrium if equal to 1 and 0 no precipitation/dissolution. In general, there are many studies on chemical-rock interaction focusing on the influence of chemical interactions on rock geochemicals, and geomechanics (Kan et al. 2005;Tomson et al. 2008;Bybee 2010;Egermann and Vizika 2011;Mohamed and Nasr El-Din 2013;Oluyemi 2014aOluyemi , 2014bWilson 2016;You and Lee 2020;You and Lee 2020). Tables 1 and 2 show a summary of previous empirical and numerical studies on chemical-rock interaction.

Direction for future research
Despite several works on sand production and oilfield interaction, fundamental research is still lacking. The fluid-rock interaction is known to cause several formation damages (14) r k = ±k − 1 due to fine migration, clay swelling, geochemical reactions between the invading chemical and the rock, resulting in rock failure (Chang and Civan 1991). Also, in changing the petrophysical properties of the reservoir formation, fluid-rock interaction can significantly alter the rock strength (Corapcioglu et al. 2014;Iriarte and Tutuncu 2018).
Although several studies have been conducted regarding oilfield chemical-rock interactions, there is still a serious void that needs to be investigated further. First, further exploration is needed to evaluate the host rock mechanical properties relative to its mineralogy and cementitious content such as clay and the composition of the interacting chemicals. The study must be extended to cover the correlation between rock composition, alterations in the formation properties, formation and fluid interface, and composition of the injected chemical. The reason is that the formation and fluid interaction is subjected to rock mineralogy and the composition of the injected fluid. Moreover, how fluid flow rate, pH, fluid composition, concentration, and temperature affect both petrophysical and geomechanical properties of the formation, and consequently rock damage, failure of the rock, and sand production, needs to be explored. Further investigation with residence time distribution tool to quantify the impact of flow rate on solute and geochemical reaction profiles during the fluid-rock interaction will help provide meaningful insight on pore-scale alterations. As a result, the real-time changes in permeability due to the injection of chemical can be elucidated.
Another area worth investigating is nanofluid formulation and its impact on rock mineral dissolution, resultant properties, and rock formation failure, including experimentation of formulated nanofluids (NFs) which will consist of dispersion of nanoparticles in the base fluid/chemical and their impact on formation upon interaction with the rock relative to standalone counterpart chemical. Without advancing characterization methods for petrophysical and geomechanical properties, it will be challenging to predict fluid-rock interactions and its impact on production and rock formation. Hence, the evolution of the formation microstructure, morphology, and grain-scale changes due to alteration in the formation physical properties as a result of chemical reaction occurring needs to be further examined using the SEM and X-ray microtomography methods. The pore network structure changes, causing changes to the petrophysical properties in the rock needs further investigation with nitrogen adsorption-desorption technique. Furthermore, a combination of modeling and microscopic characterization of the rock is required to predict the transport parameters such as porosity, diffusivity, and permeability. On the other hand, the ability to thoroughly simulate all aspects of fluid-rock interaction is expanding because of sophisticated numerical simulators capable of implementing several methods in modeling multi-scale porous in order to predict changes in the rock properties. The use of numerical simulator provides additional insight to analyze and interpret the fluid-rock interaction. In addition to insufficient database, this study approach faces challenges of limited information on the reaction mechanisms and processes happening at the surface of the mineral. Hence, the experimental measurements should be validated with numerical predictions. Consequently, the kinetic reactions occurring at the interface of the fluid rock need to be studied. Fluid-rock interactions such as wettability, porosity, and permeability changes in sandstone and carbonate reservoir rocks should be addressed experimentally. This should be validated with the combination of core flooding and pore-scale imaging simulation of the fluid-rock interaction. Zhong et al. (2019) investigated and reported the impacts of fluid flow and rock interactions on the mechanical and petrophysical properties using core samples rich in clay and another formation sample rich in calcareous material. Findings show that the extent of changes in the formation mechanical properties due to the fracking fluid-rock interaction is dependent on factors such as grain cementation and composition of the rock minerals. Although, there has been many studies on the influence of oilfield chemical and reservoir rock interface on the strength of the formation, permeability, and porosity of a rock, there is limited study on the correlation between rock composition, changes in rock geomechanical properties composition of rock minerals and chemical-rock interactions. Since different rocks will interact differently with injected fluids, it will be a good objective to evaluate the deformation and petrologic features of each type of rock and their likely interactions with various chemicals in the petroleum industry to stimulate the rock. Gomez and He (2012) proposed that the immersion test is the proper approach considered to investigate and observe the evolution of fluid and rock interface and the development of fracture.
The objective of injecting fluids into reservoir is to increase productivity. Therefore, the interface between fluid and reservoir rocks have been modeled to evaluate their contributions to formation damage (Chang and Civan 1991). The authors (Cheng and Civan, Chang and Civan 1991), in their study determined the model parameters through method of optimization by experimental data fitting. Most modeling efforts directed at formation damage addresses aspects such as permeability and porosity alteration and dissolution/precipitation phenomena. In 2019, Gardner and Sitar modeled the interface between rock formations and water by modeling distinct polyhedral blocks with the use of a discrete element method (DEM) and water was modeled with Lattice Boltzmann method (LBM). Though the model captured transient and fluid pressure, it was discovered that additional performance enhancements is required to model an accurate filed scale difficulty (Gardner and Sitar 2019). Hence, the future outlook will be to develop models that couple impacts of the prevailing processes during chemicalrock interface on geomechanical properties of the formation. This approach provides a more realistic and holistic impact of the processes occurring concurrently on the rock properties and the resultant failure rather than modeling single process, whereas the ability to apply single models is restricted in some cases than that of combined effects which cover wider areas. Notwithstanding, combined effects models will demand too much computational effort compared to single models. Consequently, the concentration of ionic species and pH on fine grains release, migration, and deposition also need to be investigated and integrated into the models. Likewise, the impact of precipitates from dissolution/precipitation reactions on rock porosity and network alteration is also another area needing further considerations. Most instabilities experienced in reservoir performance during production have been ascribed to fluid-rock interactions (Lal 1999;Muniz et al. 2005;Gomez and He 2012). It will be worthwhile therefore to develop criteria for formulating and selecting fluids that will enhance production while minimizing impacts on rock properties leading to formation damage, rock failure, and sand production. Furthermore, studies on mass transport induce changes in pore pressures due to fluid transport and ionic diffusion, and their impacts on rock formation are necessary to obtain a comprehensive stability analysis. Consequently, the microstructural and morphological changes developed within the rock can be evaluated to enhance the interpretation of impact of processes as a result of fluid-rock interaction. However, the intricate relationship between chemical reaction, fluid properties, and rock strength, the evolution of petrophysical properties and mechanical properties leading to rock failure is still poorly understood (Tenthorey et al. 2003). Hence, it will be insightful to further investigate the microstructural changes in response to stress-strain induced by the fluidrock interaction to understand hydro-mechanical changes during rock failure. On the other hand, it has been found that nanoparticle adsorption on rock surfaces in low salinity water (LSW) and flooding with nanoparticles have proven to minimize cementing material loss, mineral dissolution, ion exchange, and resistance to flow when related to only low salinity water (LSW), thereby minimizing rock formation failure (Abhishek et al. 2018). Hence, instead of using standalone chemical, formulation of nanofluids (NFs) which consist of dispersions of NPs in the base fluid/chemical can be studied as an injection fluid for enhancing production while minimizing the release of fine particles and formation damage. The NPs can alter surface and wettability of the rock (Abhishek et al. 2018); hence, the application of NFs can be considered as an area which further investigation will serve as a guide towards inhibiting the detrimental effects of formation damage associated fluid-rock interactions.

Conclusion
The interaction of oilfield chemicals and rock during enhanced oil recovery (EOR), production, reservoir stimulation, or formation treatment can affect the rock's petrophysical and geomechanical properties. These chemical injections include chemical flooding (e.g., alkali-surfactant-polymer, surfactant, surfactant-polymer), carbon dioxide flooding and sequestration, reservoir stimulation and formation treatment such as carbonate reservoir acidizing, and hydraulic fracture. It is crucial to know how specific geochemistry reactions and processes determine porosity and permeability evolution as reservoirs are subject to a wide range of chemical-rock interactions. Hence, the review of studies on the impact of oilfield chemicals on rock interface, mechanisms of the interface, and how their interactions affect reservoir rock petrophysical and geomechanical properties was reported in this paper. It was found that the oilfield chemicals exhibited strong interactions on the formation rock resulting in strong cement dissolution, which increases porosity and permeability of the reservoir formation. According to studies, chemically treated formations will release a variety of rock materials when they undergo reactions like adsorption, diffusion, precipitation, and dissolution. The extent of these reactions depends on the particle size distribution within the reservoir rock formation itself. Also, the chemical-rock interactions would vary considerably depending on temperature, the reservoir rock type, pressure, chemical composition and concentration. During these geochemical and petrophysical experiments, the problems of potential chemical-rock interactions were investigated and the consequences of such interactions on rock properties and production were reported. Resultantly, the rock petrophysical and geomechanical properties changes due to the chemical-rock interactions and fluid flow. Fine particles migration along with the fluid flow may occur due to geochemical reactions (dissolution and precipitation). This may either block the rock formation pores or create new pores resulting in changes in the mechanical strength of the reservoir formation. However, it is observed that fluid used in most of the studies are not typical oilfield chemicals and so would present unrelated chemistries since reactions would impact differently on different formation strength after chemical injection. Therefore, oilfield operators who wish to maximize oil production should study the impact of oilfield chemicals on weakening rock UCS in order to assess rock failure. In the sand production process, two stages are involved: firstly, the failure of rock matrix, followed by the detachment of sand grains and their transportation in fluid flow. Furthermore, the prediction of chemical-rock interactions and their impact on production and rock formation will be inconclusive without advances in characterization methods for petrophysical and geomechanical properties. Conclusively, it was discovered that a combination of petrographic observations with aid of advanced characterization techniques, laboratory experiments, and numerical simulations is capable of reconstructing the diagenetic process, as well as evaluating and predicting reservoir petrophysical properties due to chemical-rock interactions.
Funding The authors recognize the Tertiary Education Trust Fund (TETFund), Nigeria, for their sponsorship of the Research study. The Engineering school, Robert Gordon University, Aberdeen is also acknowledged for resources made available towards this review.

Conflict of interest
The authors declare that they have no competing interests.
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/.