Influence of depth on induced geo-mechanical, chemical, and thermal poromechanical effects

Delivering efficient and cost-effective drilled and excavated holes require effective prediction of instability along the hole profile. Most drilled and excavated hole stability analyses in the literature are performed for a given zone without considering the influence of depth. This study focused on determining the influence of depth on induced geo-mechanical, chemical, and thermal stresses and strains in drilled or excavated holes. To this end, a new porochemothermoelastic model was developed based on extended poroelastic theory, and the developed model was employed in determining induced strains and stresses for an oil and gas well case study, using data from the literature. The study delineated the different significance levels of geo-thermal-, chemical-, and thermal-induced strains and stresses as depth increased. From the results obtained, it was clear that at shallow depths, chemically induced strains and stress were the most significant formation perturbations responsible for instability of drilled and excavated holes. On the other hand, at deeper depths, geo-mechanical-induced strains and stress were the most predominant. Comparatively, thermally induced strains and stresses were found to be the least significant formation perturbations responsible for instability of drilled and excavated holes. For this case study, the results indicated that chemical strains and stresses were more prominent at depths below 170 m, accounting for more than 50% of the total stresses and strains. At 170 m, both chemical and geo-mechanical stress and strain had equal contributions to the overall stress and strain. However, as depth increased, the percentage contribution of the geo-mechanical component increased and accounted for about 80% of the total strains and stresses at 1000 m, which increased to 98.48% at depths of 6000 m and beyond. The findings of this study will provide guide for future studies on the application of extended poroelasticity theory in solving instability problems of drilled and excavated holes.


Introduction
Apart from rare surface deposits, most precious minerals and resources on earth are usually found in subsurface deposits around the world. They are normally explored using a combination of gravimetric, magnetic, and acoustic techniques. For minerals found as fluids like oil and gas, wells are usually drilled from the surface to intercept target subsurface reservoirs. However, open-pit mining was the first method of exploitation adopted for the solid minerals (Ranjith et al. 2017), as some of them were found in deposits very close to the surface. As minerals that are found relatively close to the surface were quickly depleted, exploitation has gradually become more of underground mining. Currently, there are more active underground mines of different minerals scattered around the world than open-pit mines (Ranjith et al. 1 3 2017). Figure 1 shows current deepest depths for different subsurface resources.
Before exploitation using any method, the earth barrier obstructing easy recovery of the minerals needs to be removed. This entails reaching the target deposit or reservoir through drilled holes or excavations. This action is usually achieved either by creating boreholes in the earth or by excavating large tons of earth in open-pits or tunnels. While making boreholes using drill bits, making excavations manually or using tunneling machines, a hole of a desirable size is usually made. This involves cutting through layers of earth and rocks of varying strength and weakness. One of the commonest challenges of making drilled holes or excavations in the earth is the tendency for different sections to collapse or fail. This manifests in the form of landslides in openpit mines, collapse of underground mines (Niu et al. 2011), reduction or expansion in drilled hole sizes while drilling oil and gas wells, etc. (Weijermars et al. 2020). The failure of different hole sections during or after drilling or excavation is generally called drilled or excavated hole instability.
Drilled or excavated hole instability is a condition where the size of a drilled or excavated hole becomes markedly different from the intended size, thereby resulting in the hole losing its structural integrity (Li and Weijermars 2019;Mehrabian et al. 2020). It is caused when the mechanical stress induced by drilling or excavation exceeds formation rock strength. The factors affecting the stability of drilled or excavated holes can be grouped into controllable and uncontrollable factors. Uncontrollable factors result from naturally induced abnormal formation in situ stress patterns like naturally fractured formations, tectonically stressed formations, etc. Controllable factors are usually the side effects of drilling or excavation operations, including bottomhole pressures (mud density), rock-fluid interactions, hole inclination and trajectory, etc. Figure 2 shows caliper log representations of different forms of borehole instability, including breakout, washout, and keyseat compared to an in-gauge hole. Interestingly, about 90% of drilled or excavated hole instability problems occur in clay or shale because this type of formation makes up about three-fourths of all drilled or excavated formations (Chen et al. 2003).
Clayey or shale formation is described as low-permeability sedimentary rocks with small pore radii characterized by medium porosity and medium to high clay content. Shale is known as a problematic rock in engineering applications with delicate and complicated behavior. Compared to other sedimentary rock types like limestone and sandstone, it is more complicated because a good knowledge of site engineering geology and the rock mechanical properties of shale specimen is not usually enough for successful engineering design (Farrokhrouz and Asef 2013 Farrokhrouz and Asef (2013), the different forms of drilled or excavated hole instability associated with shale or clay can be grouped into three types. The first group comprises mechanical in situ stresses and rock strength. Also, there are chemical sensitivity and swelling resulting from interaction with an invading fluid, in addition to thermal-induced stresses. Rock failure mechanisms become complicated when thermally induced stresses due to invading fluid-formation temperature differences are considered. Carter and Booker (1982) carried out one of the earliest pore pressure and time-dependent stress analyses for excavations. They developed analytical solutions for circular tunnels in fluid-saturated media. Detournay and Cheng (1988) also developed non-hydrostatic analytical solutions for vertical wells. They showed that shear failure could also be initiated both on wellbore walls and within the formation. Cui et al. (1997) provided analytical solutions for isotropic porous media while considering wellbore inclination. The solutions were later extended to account for transitropic conditions (Abousleiman and Cui 1998) and time-dependent effects (Ekbote et al. 2004). Wilson and Aifantis (1982) developed hydrostatic dual-porosity analytical solutions for vertical wellbores in naturally fractured formations, while Li (2003) developed similar solutions but considered nonhydrostatic stress conditions. Abousleiman and Nguyen (2005) presented threedimensional analytical solutions for inclined wells in naturally fractured formation. Dual-porosity consideration predicted narrower mud-weight windows than single-porosity Nguyen and Abousleiman 2009). Other analytical solutions were equally developed for inclined wells while considering chemo-electrical and thermal effects (Nguyen and Abousleiman 2010;Abousleiman 2005, 2006).
The building block of most of these borehole instability models is Biot's poroelasticity theory of 1941. Although the theory was initially used in describing three-dimensional soil consolidation (Biot 1941), it was hugely successful that it had found widespread applications in diverse areas. These new areas of application are usually replete with porous materials having fluid-filled pores, including soil, rocks, bones, etc. (Cheng 2016). Refer to Fig. 2 description: Caliper log response due to borehole enlargement (Reinecker et al. 2003). Following the successful application of poroelasticity theory to borehole stability analysis, other extended forms of the theory quickly emerged. This was borne out of the fact that pure poroelastic analysis could not possibly explain borehole instability associated with chemically reactive formations and regions with unusually high subsurface temperatures (Cheng 2016). Numerous experimental pieces of evidence describing how chemical inhomogeneity and thermal imbalances downhole equally cause induced pore pressures have been presented (Kanfar et al. 2017). These various dimensions of addressing borehole instability led to poromechanics analyses under different names like porochemoelasticity (Ekbote and Abousleiman 2006), porothermoelasticity (Jaeger et al. 2010), porochemothermoelasticity (Ekbote and Abousleiman 2005), etc. The title of each suggesting the different types of porophysical phenomena is considered (Figs. 3, 4).
One of the commonest poromechanics phenomenon now considered in many extended forms of Biot's poromechanics theory is the chemical expansivity of argillaceous formations like shale on exposure to water. The chemical reactivity in these formations has been attributed to the presence of clay particles, which are responsible for volumetric changes that depends on water activity imbalances. When the formation water activity is higher than the water activity of invading fluid, positive volumetric changes (swelling) occur and vice versa (Luo et al. 2017). However, according to Cheng (2016), the mechanism of clay swelling can be divided into inner crystalline and osmotic swelling. Inner crystalline swelling occurs when water molecules invade the interlayer spaces between the silicate sheets of chemically sensitive formations in a process called absorption. This absorption process is known to result in swelling stresses of up to 100 N/mm 2 . But osmotic swelling takes place when electric double layers of silicate sheets repel themselves causing swelling, while water molecules trapped in this cloud of ions is said to be adsorbed. Osmotic swelling occurs over a larger distance compared to inner crystalline swelling, but its generated swelling stress is smaller than inner crystalline swelling stress-less than 2 N/mm 2 (Cheng 2016).
While it is true that one porophysical phenomenon might be significant in one area and insignificant in another, the fact remains that a complete analysis should include as many phenomena as possible. More so, a tripartite analysis involving geo-mechanical, chemical, and thermal variables becomes imperative given that the effects of these phenomena change as depth increases (Han et al. 2019). For instance, geological evidence abounds on the gradual transformation of compacted shale (found at shallower depths) to cemented shale (commonly found at deeper depths) as depth increases (Farrokhrouz and Asef 2013). This is exemplified by the conversion of smectite to illite with increased temperature and pressure at greater depth. Also, considering common geological temperature gradients, the effects of drilling mud temperature in oil and gas wells become more significant as greater depths are drilled (Jaeger et al. 2010).
However, most research on the analysis of borehole instability based on poroelasticity theory or one of its extended forms like Ekbote and Abousleiman (2005), Kanfar et al. (2017) and Li and Weijermars (2019), only consider effective hoop, radial, and tangential stresses at a given zone without considering the influence of depth. In addition, other poromechanicsbased treatments of drilled or excavated hole instability in the industry usually do not isolate induced stresses and strains from residual stresses and strains. This seemingly overlooked void is what this study seeks to fill using isotropic conditions and a simplified procedure. By focusing on induced stresses and strains in this study, due attention is given to the part of the total stresses and strains really responsible for instability. First, a novel model for borehole stability is developed, and the effect of depth on induced stress and strain is evaluated.

Model development
When carrying out poromechanical studies, a number of assumptions are frequently made. Following the classical Biot's poroelasticity theory, the following assumptions were made in this study: An isotropic, porous, fluid-filled continuous solid medium is assumed with a Cauchy stress distribution. Hence, average stresses on planes within and on its surface are treated as forces per unit area. Essentially, the stresses are also assumed to be borne by both the infiltrating fluid and the solid skeleton. Only small strains are considered, meaning that components of strain and the projections of displacement with its first derivatives in a given direction are linear. The solid medium follows the generalized Hooke's law and, therefore, its stress-strain relations are equally linear. Apart from mechanical homogeneity, the solid also exhibits chemical and thermal homogeneity, such that it deforms in a geometrically self-similar way. The fluid occupying the pore spaces does not chemically react with the solid skeleton. Upon imposition of the various porophysical stresses, only the final equilibrium states are important. As such, viscoelastic effects and time-dependency are not considered.
Assuming only small strains, the general Hooke's law for an isotropic elastic body can be written in principal coordinates as (Jaeger et al. 2010): where 1 , 2 , 3 are three principal stresses and 1 , 2 , 3 are the principal strains. Also, G is the shear modulus, and is the first Lamé parameter.
Recall that volumetric strain equals the sum of three normal strains, Then Eqs. (1)-(3), on substituting for v , can be written as: Consequently, after dropping i in favor of the more descriptive matrix terms, ij , for stress, Eqs. (5)- (7) is written more compactly in a matrix form as: When written out term by term in the general case when the coordinate system is not aligned with principal axes, Eq. (8) takes the following form: On inverting the above expressions of Hooke's law and making strains the subject, we have: where E = modulus of elasticity, v = Poisson's ratio and

Constitutive equations
To develop a linearized, non-hydrostatic theory of porochemothermoelasticity, we start by recasting Hooke's law in terms of strain using shear modulus G and Poisson ratio v: The above equations can be reduced in matrix notation: where the identity matrix, I , is given by:

Components of total strain
Stresses and strains are considered thermodynamic conjugates, such that a given strain is the result of its corresponding stress (Jaeger et al. 2010). Hence, for a tripartite treatment involving geo-mechanical, chemical, and thermal considerations with their respective stresses, there are corresponding geo-mechanical, chemical, and thermal strains. Therefore, the total strain, comprising effects of geo-mechanical, chemical, and thermal stresses, is taken as the summation of pore strain, hydration strain, and thermal strains. Pore strain refers to both positive and negative reductions in pore volume in a given direction due to positive and negative changes in the pore structure under geo-mechanical effects. The pore strain was used to represent geo-mechanical strain since for mechanically homogenous materials, pore strains are consequences of changes in geo-mechanical stresses (Jaeger et al. 2010). Also, chemical strains describe both positive and negative reductions in volume in a given direction due to invading fluid (water or drilling mud) and formation chemical interactions. Thermal strains are both positive and negative reductions in volume in a given direction as a result of the temperature difference between invading fluid and formation (Cheng 2016).
If pore strain is given as −C bp P p ∕3 (Jaeger et al. 2010), and chemical strain as − ( s − d )I (Luo et al. 2017), then thermal strain is given as (T − T o )I (Farrokhrouz and Asef 2013). Here, is a property of a given set of shale formation and the invading fluid. The invading fluid here could be drilling mud used in oil and gas drilling, rain water infiltrating exposed open-pit mine walls or subsurface aquifer water interacting with the walls of underground mines. The parameter also incorporates the membrane efficiency (I m ), measuring the suitability of the shale or clay as a semipermeable membrane (Farrokhrouz and Asef 2013). Combining the respective strains with the expression for strain, Or expressed differently in terms of effective strain, On successful inversion operation to express Eq. (22)  = thermal expansivity (K −1 ); = temperature difference between drilling fluid and formation (K); R = universal gas constant (J mol −1 K −1 ); V mw = molar volume of water (m 3 / mol); a wd , a ws = drilling mud and shale water activities (-).

Materials and methods
From the developed model, the terms on the left-hand side of Eqs. (25)-(27) are called effective stresses (Jaeger et al. 2010). The first term on the left-hand side of the equations is the stress in a given direction, and the remaining terms are induced pore pressure, chemical, and thermal stresses, which change with depth. Hence, the effective stress can be separated into a residual stress part and an induced stress part (Ranjith et al. 2017), where the residual stress represents the virgin in situ stress in a given direction before the initial stress equilibrium was upset through drilling or excavation, with the induced stress capturing the magnitude of geo-mechanical, chemical, and thermal imbalances introduced by the drilling or excavation operation. Similarly, for the effective strain equation, the residual and induced parts can equally be identified in Eq. (23). It can also be noticed from Eqs. (23) and (25)-(27) that the induced stresses appear only in the principal stress directions. This is because induced stresses are purely volumetric, and as result do not support shear stresses (Jaeger et al. 2010). From the signs of the induced stresses and strains, the equations show they are generally tensile (extensive). Hence, the induced stresses and strains were calculated from the induced parts of the effective stress and strain equations as given in Eqs. (23) and (24), while the percentage contributions of each porophysical phenomena considered were calculated from their respective contribution to the total stress and strain.
In addition, for the stress and strain sign convention, compressive stresses and strains were treated as positive, while tensile (extensive) stresses and strains were taken as negative in agreement with geomechanics convention. Table 1 shows the parameters used for the analysis in this study. The assumed values are consistent with field experience. Given available data, the case study analyzed is that of oil and gas drilling, but the results easily extend to excavations and tunnels in mines and construction works. Also, the case study involves a vertical well, but the same treatment can be extended to inclined wellbores.
Moreover, a complete description of how the total pore pressure varies depends on another set of phenomenological equations (Demirel 2007) needed for a mathematically well-posed borehole stability problem (Jaeger et al. 2010). These phenomenological equations and the accompanying Navier-type equations are beyond this study's scope. A complete porochemothermoelastic solution to borehole instability can also be found in other studies (Ekbote and Abousleiman 2005;Kanfar et al. 2017).

Results and discussion
Figures 5, 6, 7, and 8 show the geo-mechanical-, chemical-, and thermal-induced strain and stress variations with depth. In other to fully represent the observed trends, the depths were split into two: shallow and deep depths. For this study, shallow depths were taken as depths below 1000 m and deep depths were taken as depths above 1000 m. Hence, Figs. 5 and 6 are used to represent variations of geo-mechanical, chemical and thermal strains and stresses with depth, respectively, for shallow depths, while Figs. 7 and 8 are used to represent variations of geo-mechanical, chemical, and thermal strains and stresses with depth, respectively, for deep depths.  (2013) From Figs. 5 and 6, it can be seen that total induced strains and stresses as a summation of three induced strains and stresses components vary directly with depth. This clearly shows that induced strains and stresses become significant, with increase in depth. It can equally be seen that both geo-mechanical and chemical strain and stress constituted a major part of the total induced strain and stress in drilled holes and excavations at shallow depths compared to thermal strain and stress. In fact, at depths below 170 m, chemical strains and stresses accounted for the majority of the total strain and stress. Beyond this depth, geo-mechanical-induced strains and stresses clearly increased as depth was increasing, following the trend of the total strain and stress. The rate at which geo-mechanical strain and stress was increasing was fast enough that geo-mechanical strains and stresses quickly dwarfed the rates at which both chemical and thermal strains and stresses were changing. Making chemical and thermal strains and stresses look seeming constant over the depth range being considered.
Similar trends were also noticed for deep depths in as captured in Figs. 7 and 8. Here, geo-mechanical strain and stress followed the total strain and stress more closely compared to shallow depth, with chemical and thermal strain and stress accounting for even smaller proportions of the total strain and stress. This shows that at deeper depths, geomechanical strain and stress account for over 90% of induced stresses and strains in drilled holes and excavations. This observation agrees with typical geological gradients, which indicate that geo-mechanical variables generally increase with depth (Ranjith et al. 2017). Also, the disparity between the induced strain and stress responses of the different poromechanical components at shallow and deep depths can be traced primarily to increase in overburden stress accompanying increase in depth. A careful look at depth profiles of drilled holes and excavations will reveal different layers of formation systematically placed one above the other, with the weight of overlying sections carried by lower ones. This overlying weight becomes increasing higher at increasing depth as underlying formation carries increasing overlying weight, giving rise to geo-mechanical effects as represented by in situ stresses in different directions (Dosunmu et al. 2020). Since the underlying rocks cannot undergo full strain due to adjacent rocks, this overburden stress is also translated into both maximum and minimum horizontal geomechanical stresses.
In order to get a clearer picture of the variation of chemical and thermal strain and stress with depth over the whole depth profile, Figs. 9 and 10 are plotted. Figures 9 and 10 reveal that chemical strain and stress were actually increasing with depth, but the rate of increase was below the rate of geo-mechanical strain and stress increase. On the other hand, Figs. 9 and 10 equally reveals that thermal-induced strains and stresses, however, reduced as depth increased. This reduction in thermal-induced strains and stresses does not entirely signify reduction in magnitude per se. For this case study, it was due to the changing signs of drilling mud-formation temperature difference, and not necessary a decrease in magnitude. Despite there was a decrease in magnitude of the thermal strain and stress as the tensile strain and stress reduced to zero. But as soon as the strain and stress turned negative signifying the onset of compressive strain and

Chemical
Thermal Total stress, the strain and stress magnitudes then started increasing with depth. The drilling mud-formation temperature difference gradually changed sign from positive to negative as depth increased, thereby changing thermal strain and stress from tensile or extensive at lower depths to compressive at deeper depths. The drilling mud temperature will likely be higher than formation temperature at few meters depths, especially in cold or temperate regions. Another good reason for possible higher mud temperatures at shallow depths could equally be the pumping energy delivered by the pumping system on a drilling platform. Some of which would end up being utilized to overcome frictional resistance to flow, thereby increasing the mud temperature. Since thermal strains and stresses result from thermal imbalance between drilling mud and formation, higher mud temperature will lead to tensile strains and stresses on the formation at shallow depths. At deeper depths, the thermal strain and stress become compressive (negative strains and stresses), due to likely higher temperatures. Given that Figs. 5,6,7,8,9,and 10 also show that the rates of change in the magnitude of strain and stress for various components of induced strains and stresses were different with increase in depth, which is because each strain and stress component has a depth range for which its impact is significant (Bourgoyne et al. 1991). These varying rates of change in strain and stress informed the choice to investigate the percentage contribution of each induced strain and stress component to the total as depth increases. Figure 11 shows the relative contribution of the three induced strain and stress components.
From the statistics, it is evident that the total percentage contributions of all induced strains and stresses were equal to 100% at all depths considered. This total percentage induced strain and stress were used as checks for accuracy, as summation values different from 100% would have flagged

Depth (m)
Chemical Thermal calculation errors. However, chemical-induced strains and stresses accounted for significantly high percentage of total induced strains and stresses below 1000 m. This was because thermally and geo-mechanically induced strain and stress percentage contributions were comparatively low below this depth. At greater depths, the percentage contributions of chemical-induced strains and stresses to total strain and stress gradually decreased, while geo-mechanical strain and stress astronomically increased with depth. This phenomenon can be explained by the gradual transition from compacted or consolidated formation (including smectites, which are more chemically sensitive) to cemented formation (like illite, which is less chemically sensitive) as drilling depth increases. As more cemented formation is encountered at deeper depths, the percentage contributions of chemical strains and stresses to total strain and stress decrease due to a substantial reduction in porosity (Farrokhrouz and Asef 2013). This is because most wellbore stability problems are traceable to compacted chemically reactive smectite group of shales, which are usually found at shallower depths (Weijermars et al. 2020). With an increase in depth, both temperature and pressure increase, thereby increasing formation pore pressure simultaneously. Pore pressure serves as a good measure of geo-mechanical-induced strains and stresses (Jaeger et al. 2010). This increase in pore pressure due to increasing pressure and temperature following higher geo-mechanical stress at greater depths explains why the percentage contributions of geomechanical-induced strains and stresses to total induced strain and stress increased.
The percentage contributions of thermal-induced strains and stresses were the least among the three induced strain and stress components. The relatively low mud-formation temperature difference is perhaps not high enough to significantly contribute to total induced strain and stress compared to chemical and geo-mechanical components. Except where drilling mud is continuously cooled to achieve certain drilling objectives, the average mud temperature is not usually very far from average formation temperatures due to continuous conductive heat transfer along the wellbore wall. This possibly explains why thermal strains and stresses are sometimes ignored in some borehole stability analyses, especially in shallow holes and excavations. However, thermally induced strain and stress are known to be significant around radioactive waste canisters or during injection of cold water in geo-thermal wells for pressure maintenance (Jaeger et al. 2010).
For this case study, it was noticed that the chemical component of induced strain and stress accounted for more than 50% of total strains and stresses, at depths below 170 m. Also, at 170 m, the percentage contributions of chemicalinduced strains and stresses were the same with that of geo-mechanical strains and stresses. Particularly, the two poromechanics effects both had 46% contributions each to the total strains and stresses, with thermal-induced strains and stresses accounting for the remainder. At 1000 m, geomechanical strain and stress component were by far the main contributors to total strains and stresses, accounting for about 70% of the total strain and stress. Chemical strains and stresses accounted for less than 30% and thermal-induced strains and stresses accounted for the balance. However, the percentage contribution of the geo-mechanical component increased to 98.48% at 6000 m and beyond. At this depth, the percentage contributions of the poromechanics effects of chemical and thermal strains and stresses were 2.56% and 1.03%, respectively. The results support general drilling engineering rule of thumb that most borehole stability problems encountered at mid to shallow depths are likely attributable to chemical reactivity of shale. Consequently, this study showed that at shallow depths, chemically induced strains and stresses are more likely to contribute to borehole instability, while geo-mechanically induced strains and stresses are mainly responsible for borehole instability at greater depths. Also, thermal-induced strains and stresses are less likely to result in drilled or excavated hole instability problems.

Conclusion
From the discussion above, the following conclusion can be drawn: (a) A comprehensive drilled or excavated hole stability analysis should involve as many poromechanics phenomena as possible. Neglecting any of them in other to simplify analysis might be misleading, as even the smallest stress unaccounted for could still pose tremendous borehole stability problems. (b) Before any given poromechanics phenomenon is declared insignificant, a stress/strain percentage contribution analysis for the intended depth profile should be done. This will reveal the significant impact of each of them on borehole stability at each depth profile zone. (c) At shallow depths, chemically induced strains and stress are the most significant formation perturbations responsible for drilled or excavated hole instability. (d) At deep depths, geo-mechanical-induced strains and stress are the most significant formation perturbations accounting for drilled or excavated hole instability in drilling operations or excavations. (e) Thermally induced strains and stresses, although significant as depth increases, are the least significant formation perturbations responsible for drilled or excavated hole instability. (f) However, the remarkably similar results of strain and stress analyses support the fact that stress and strain are exact thermodynamic conjugates.

Appendix
See Tables 2, 3, 4, and 5.  Author's contribution All the authors contributed to the final version of this manuscript. ENE did model development, formal analysis, investigation, and wrote the first draft. The first and final drafts of the article were reviewed and edited by AVJ. Research conceptualization and funding was secured by JFO. At the same time, supervision was done by AVJ, DA, and JFO.
Funding This research benefitted from sponsorship of Petroleum Technology Development Fund (PTDF), Nigeria, with Grant Number PTDF-SP&D-AOGRG-V.III-117-011B. Grant is for experiment, data gathering, processing, and report writing.

Availability of data and materials
The data associated with this manuscript was obtained from journal articles published under creative commons license. And as such, the use of data thereof still falls under the applicable approval.

Conflict of interest
The authors hereby declare no inherent conflict of interest in this study.
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/.