Utilizing Machine Learning to Evaluate the Connection between Poisson's Ratio and the Petrophysical Properties of Reservoir Rocks

Al-Obaidi SH, Khalaf FH and Viktor S

Published on: 2025-04-19

Abstract

The Poisson's ratio is a crucial cornerstone, illuminating our understanding of geomechanically behaviour in wells during the dynamic drilling process and the inspiring recovery journey. This research rigorously employs machine learning methods to analyse the significant impact of geophysical parameters on the Poisson ratio in hydrocarbon reservoirs found in oil fields. The analysis utilized data from multiple oil and gas fields, highlighting the crucial relationships between the Poisson ratio, the natural radioactivity of rocks, and the velocity of the longitudinal wave. Understanding these dependencies is essential for optimizing extraction processes and improving resource management. The PIK-UIDK/PL unit was crucial in conducting triaxial tests on samples in reservoir conditions. The insights gained from these tests have led to essential dependencies that underscore their significance.

Numerous geophysical well surveys are utilized to establish robust equations that define the relationship between the Poisson ratio and geophysical parameters. This is achieved through effective linear regression and advanced machine-learning methods. Based on these dependencies, the Poisson ratio can be assessed more accurately for various rocks and fields. In addition to improving efficiency in hydrocarbon production and optimizing oil industry operations, the findings can be used to improve forecasts and modelling for oil field development processes.

Keywords

Poisson's ratio; Longitudinal wave; Machine learning; Petrophysical properties; PIK-UIDK/PL unit

Introduction

The static Poisson's ratio (ν) is one of the most significant and at the same time the least defined parameter in the calculation of the stress state. The study of the static Poisson's ratio is carried out everywhere during the implementation of almost any project on geomechanics [1-4]. Still, there is no proven method for distributing this parameter along a geological section based on borehole and seismic 3D studies. Therefore, despite the experiments carried out, in most cases, researchers use the dynamic Poisson's ratio, which practically does not change within one lithotype and is used, as a rule, by geophysicists as a characteristic determining the lithofacies composition of the strata [5-8]. At the same time, almost all experts point to the undeniable importance of the static Poisson's ratio in geomechanical calculations, particularly in hydraulic fracturing design [9-12].

At present, neither domestic nor foreign literature contains reliable universal relationships between Poisson's ratio and the parameters most frequently used in geomechanics, such as the elastic modulus, compressive and tensile strength limits, adhesion, internal friction angle, and longitudinal and transverse wave velocities [13-16]. There are relationships for individual deposits that cannot be extended even to neighbouring deposits due to the low tightness of the relationship between the obtained relationships. One reason for the lack of reliable dependencies is the complexities in determining this parameter, which is regulated by various domestic and international standards.

According to Method for Determining Deformation Characteristics under Uniaxial Compression, you can determine the deformation characteristics of rocks only when they are tested in a uniaxial mode [17-20]. ASTM standards describe uniaxial and triaxial tests [21,22], and ISRM recommendations have separate documents on determining strength under triaxial loading [23,24]. Many methods do not contain information on determining deformation characteristics. The issue of the difference in determining elastic parameters under uniaxial and triaxial loading is not considered in these methods. According to [25-28], tests are carried out to determine the modulus of elasticity, Poisson's ratio, deformation modulus, and transverse deformation coefficient. It is indicated that these characteristics should be defined in the range of required stresses from 5% to 50% of the ultimate strength under uniaxial compression. The ISRM and ASTM methodologies state that the modulus of elasticity can be determined using any method employed in engineering practice and outline several more commonly used techniques [29-32]. The first option is the tangential modulus of elasticity, determined at a fixed level from the ultimate strength (usually 50% of the ultimate strength in uniaxial compression). The second is the average modulus of elasticity on the linear portion of the stress-strain diagram. The third option is the secant modulus, determined at a fixed level from the tensile strength (also usually 50%). Poisson's ratio is calculated based on the ratio of the axial-to-radial deformation in the section where the elastic modulus is determined. ASTM standards state that elastic moduli can be calculated directly and by the best approximation of the section of the stress-strain diagram, which is under consideration by a straight line using the least squares method.

Thus, the spread of methods for determining the loading section of the sample, where the elastic modulus and Poisson's ratio are located, can lead to a significant difference in the parameters being set. And if this is not so significant for the elastic modulus due to its large boundaries - from several to hundreds of GPa- then for Poisson's ratio, a difference of 0.1 leads to the impossibility of obtaining reliable correlation dependencies.

Review of Previously Obtained Dependencies

It should be noted that in some cases it is possible to obtain relatively stable dependencies on individual geophysical parameters, primarily on the longitudinal wave velocity and the parameter characterizing the natural radioactivity of the section. Figure 1 shows the dependence of the static Poisson's ratio of sandstones and siltstones of oil fields in the Perm Territory on the intensity of gamma radiation [33,34]. As we can see, the distribution of the Poisson's ratio along the productive section of the object can be constructed based on the intensity of gamma radiation of terrigenous rocks.

Figure 1: Dependence of the static Poisson's ratio of sandstones and siltstones of oil fields in the Perm region on the intensity of gamma radiation.

Figure 2 shows a similar dependence of the Poisson ratio of productive formations of the Sarmatskoye oil field on the indicator characterizing the natural radioactivity of rocks. Samples for testing are represented by various lithotypes (argillites, dolomites, clayey limestone, etc.). A pattern was revealed when the test results were processed: the Poisson ratio for argillite in reservoir conditions reaches 0.35, and it does not exceed 0.25 for limestone. Moreover, the lower the clay content, the lower the Poisson ratio. Based on this, an attempt was made to determine the relationship between the static Poisson ratio and wells' gamma-ray logging (GR) data. As a result, for all tested samples in three wells, a relationship was built between the static Poisson's ratio and the GR data with a fairly high indicator of the tightness of the relationship. It should be noted that in earlier studies, a relationship was found between the Poisson's ratio and the volume content of clays for the Adamtash deposit (Central Asia), as well as with the intensity of gamma radiation in sandstones, siltstones, and argillites, which is determined by the content of clay minerals [35-38].

Figure 2: Dependence of the static Poisson ratio on the indicator characterizing the natural radioactivity of rocks.

Figure 3 shows the dependence of Poisson's ratio on the longitudinal wave velocity for samples of productive objects of the Salmanovskoye gas condensate field. The samples are represented by formations with average occurrence depths of 900 m (PK1-9), 1600 m (XM6-10) and 1600-2900 m (TP). As can be seen from the graph, the dependence is characterized by a relatively low connection (R2 = 0.42). Note that this dependence can be used for approximate calculations of lateral rock pressure.

Figure 3: Dependence of Poisson's ratio of samples of objects (formations) of the Salmanovskoye deposit on the velocity of the longitudinal wave. The point highlighted in red is excluded when constructing the dependence.

A similar dependence is typical for samples of Bazhenov deposits of the Sredne-Nazymskoye field, oriented perpendicular to the bedding. Most samples have longitudinal wave velocity Vp values in the range of 2500–4000 m/s. For this section, Figure 4 shows the dependence of Poisson's ratio on Vp. The points excluded from the construction of dependencies are marked in red in Figures 3 and 4. They may be outliers for the following reasons: significant heterogeneities in the samples, the influence of other factors unrelated to the longitudinal wave velocity, etc.

Figure 4: Dependence of Poisson's ratio of Bazhenov sediment samples of the Sredne-Nazymskoye field, oriented perpendicular to the bedding, on the velocity of the longitudinal wave. The points highlighted in red are excluded when constructing the dependence.

Thus, for individual deposits, the dependences of the Poisson ratio on the longitudinal wave travel speed and the GR log parameter are constructed. However, it is not possible to obtain sufficiently reliable dependences for a group of deposits using conventional statistical methods. Due to the limited dependencies of the Poisson ratio on the geological, geophysical, and geomechanical characteristics of the section, individual researchers are exploring the relationship between geophysical parameters and the Poisson ratio using machine learning techniques [39-42]. This paper presents the results of similar studies using machine learning to establish the dependence of the Poisson ratio on the geophysical characteristics of the studied sections of Western Urals deposits.

Methodology and Materials

Finding Dependencies Using Machine Learning Methods

A total of 156 tests were conducted on different sections of the Western Urals deposits; however, due to an uncertain signal when determining the longitudinal wave velocity, some samples were excluded from the final analysis. The final selection consisted of 145 samples. It contains samples of both terrigenous and carbonates rocks. Triaxial compression experiments were performed on the unique scientific installation PIK-UIDK/PL. The tests were conducted at effective pressure; pore pressure was not specified. The effective pressure corresponding to reservoir conditions was calculated by considering the dependencies for the Biot coefficient previously obtained [43-46]. Axial loading and unloading were performed with a constant rate of axial deformations, the same for all experiments, which was 10–5s–1. The rate of change of lateral pressure did not exceed 1 MPa/min.
The general program of triaxial testing of samples for elastic properties under reservoir conditions consists of the following stages;

  • Loading the sample to a specific all-around pressure that matches the reservoir conditions for the given sample.
  • Holding the sample until volumetric deformations stabilize

                           --------    (1)

  • Measurement of longitudinal and transverse wave velocities Vp, Vs.
  • Production of two sample loading/unloading cycles by creating an axial load without exceeding the elastic limit.

Obtaining the straightest possible line on the stress-strain graph on the load branch determines a sufficient load level. The elastic deformation limit is taken to be the bend in the curve in the axes "axial load" - "volume deformations". Volume deformations are calculated as the sum of axial and radial deformations. Continuous recording of axial stress, axial and transverse deformations, and all other parameters is carried out throughout the experiments. The results of all tests were reprocessed using a single method to exclude the influence of some processing features, such as the section size for determining elastic parameters. Differences may arise if different people process the test results at different times. Additional parameters, including the load level for determining elastic parameters, were also obtained for plotting the dependence of Poisson's ratio.

Results And Discussions

The elastic modulus and Poisson's ratio were calculated on the straight sections of the branches of the repeated loading of the sample on the diagram "axial stresses - deformations" in the region of elastic deformations. An example of the diagram of the deformation of the sample with the indication of the section for determining the elastic parameters is shown in Figure 5.

Figure 5: Example of a sample deformation diagram indicating the area for determining elastic parameters.

The dynamic elastic moduli and Poisson's ratio were determined from the velocities of longitudinal and transverse waves under reservoir conditions using established relationships [47-50]. The next step in data preparation for training involves selecting parameters and refining them into the required format. The parameters chosen for training were the well-logging data obtained throughout the wellbore, which can potentially be normalized. Before utilizing the well-logging data, it was essential to normalize them collectively. However, before normalization, the data needed to be adjusted for a more accurate comparison with the results from point tests. Figure 6 shows a heat map of pairwise correlation coefficients, which shows all the selected parameters, their cross-correlations, and their correlations with the Poisson ratio. It is clear that, for parameters with low correlation coefficients with the Poisson ratio, their cross-pairwise correlations are quite high in absolute terms, which is not always acceptable when applying machine learning methods.

Figure 6: Heat map of pair correlations: GR- gamma ray, v–propagation velocity of longitudinal wave, w – hydrogen content, error – error in determining Poisson’s ratio, load1 – load at which it was determined, poisson – Poisson’s ratio.

In regression analysis, multicollinearity is critical because it affects the model's ability to distinguish between unique values of the regression coefficients, which can distort their estimates. Small changes in the data used for estimation can cause large fluctuations in the estimated coefficients. This creates difficulties with inference: if two variables are highly correlated, a change in one variable can be offset by a change in the other, cancelling their effect. In gradient optimization methods, an unlimited increase in the model weights is also possible, which makes the search for optimal parameters of the optimization function impossible. In this case, it was decided to exclude the hydrogen content (w), as the parameter most correlated with acoustics (Vp), while having a lower correlation coefficient value with the target parameter. The analysis also used the GR, errors in determining the Poisson ratio when approximating a section of the curve with a straight line (error) and the load at which it was determined (load1) were taken into account. The search for a relationship was carried out using linear regression with the Huber loss function, which allows for excluding outliers. Various combinations of the parameters shown in Figure 6 were used as input parameters. The parameters were added to the model one by one, after which the model was retrained and the cross-validation error was estimated. Parameters were added until the error significantly decreased. The final cross-plot is shown in Figure 7. The resulting equation included the longitudinal wave travel speed and GR data in various variations. The error in determining the Poisson ratio (error), which was present in the regression model, is taken to be zero in the resulting equation. The correlation coefficient of the resulting regression dependence was 0.74, which is a good result in this case (for such an uncertain parameter). Figure 7 shows the cloud of obtained values, according to which confidence intervals were constructed.

Figure 7: Cross-plot of the determined from tests (measured Poisson's ratio) and calculated Poisson's ratio for terrigenous and carbonate productive objects of the Western Urals deposits: A – Tournaisian stage of the lower section of the Carboniferous system together with the Famennian stage of the upper section of the Devonian system, B – Tulskian and Bobrikovian horizons of the Visean stage of the lower section of the Carboniferous system, C – Bashkirian stage of the middle section of the Carboniferous system.

The lower confidence interval estimate should be used as the final result since it is the best approximation of the experiment's drained conditions. Such conditions imply a sufficiently low loading rate at which the sample's pore pressure does not increase. Under undrained conditions, the liquid does not have time to flow out of the sample due to low permeability or a high rate of load application.

The equation for the Poisson ratio, obtained from the lower confidence interval estimate, is

                                                                    v=2.34x10^(-5) V_p-6.14x10^(-2) Ln(GR+1)+3.29x10^(-10) V_p^2 GR+0.166  --------    (2)

Where GR is the natural radioactivity of rocks determined using GR log, μR/h, VP is the longitudinal wave velocity determined using well logging, m/s. Thus, using machine learning methods and processing heterogeneous test data using a single methodology, it is possible to obtain fairly stable correlations between the Poisson ratio and the geophysical characteristics of the productive section.

Conclusions

The Studies Conducted Lead Us to the Following Conclusions

Currently, there are no reliable dependencies of the Poisson ratio on the geomechanical and geophysical characteristics of productive objects. One reason is the lack of a strict, standardized method for determining this parameter. For specific productive types of deposits, correlations have been identified that connect the deposit properties with the velocity of longitudinal waves, the natural radioactivity of rocks, and the clay content. However, these correlations do not apply to other types of deposits. This paper proposes a machine learning-based method for establishing the Poisson ratio's dependence on the main geophysical parameters of the geological section. With this method, the correlation dependencies of this parameter on the longitudinal wave velocity and the GR data were obtained for terrigenous and carbonate-producing objects of the Western Urals deposits. The correlation coefficient of the obtained regression dependence was 0.74, which is a good result.

References

  1. Tao Z, Zhu C, He M, Karakus M. A physical modeling-based study on the control mechanisms of Negative Poisson's ratio anchor cable on the stratified toppling deformation of anti-inclined slopes. International J Rock Mechanics and Mining Sciences. 2021; 138: 104632.
  2. Obaidi SHA, Chang WJ. Techniques to Boost Oil Production in the Development of Multi-Reservoir Fields. J Earth and Environmental Waste Management. 2024; 2: 1-7.
  3. Kamensky IP, Al-Obaidi SH, Khalaf FH. Scale effect in laboratory determination of the properties of complex carbonate reservoirs. Int Res J Mod Eng Technol Sci. 2020; 2: 1-6.
  4. Al-Obaidi SH, Patkin AA, Guliaeva NI. Advance Use for the NMR relaxometry to investigate reservoir rocks. J Petrol Eng Technol. 2003; 2: 45-48.
  5. Fanny D, Gomord OF, Vandycke S, Schroeder C, Swennen R, Tshibangu JP. Relationships between geomechanical properties and lithotypes in NW European chalks. Geological Society, London, Special Publications. 2017.
  6. Al-Obaidi S. Investigation of rheological properties of heavy oil deposits. In Conference of the Arabian J of Geosciences. Cham: Springer International Publishing. 2019; 399-402.
  7. Miel H, Hameed AOS, Hussein KF. Modeling and monitoring the development of an oil field under conditions of mass hydraulic fracturing. Trends in Sciences. 2022; 19: 3436.
  8. Al-Obaidi SH, Kamensky IP, Hofmann M. Changes in the physical properties of hydrocarbon reservoir as a result of an increase in the effective pressure during the development of the field. engrXiv. 2021.
  9. Truong TV, Nguyen VT, Nguyen LK, Nguyen TV, Nguyen HT, Nguyen TT, et al. Prediction of Poisson's ratio for hydraulic fracturing operations in the Oligocene formations in the Bach Ho field. J of Mining and Earth Sciences. 2024; 65; 47-57.
  10. Al-Obaidi SH. Experimental study of the influence of fluid flow rate on the risk of rock destruction. J Pet Eng Emerg Technol. 2009; 9: 84-89.
  11. Chang WJ, Al-Obaidi S, Smirnov VI. Application of Multi-Stage Hydraulic Fracturing (MSFH) in the Development of Deep-Lying Deposits with Low Permeability of Hydrocarbon Fields. J Geol Geophys 2023; 12: 1154.
  12. Al-Obaidi S. Analysis of hydrodynamic methods for enhancing oil recovery. J of Petroleum Engineering and Technology. 2007; 6: 20-26.
  13. Gonzaga GG, Leite MH, Corthésy R. Determination of anisotropic deformability parameters from a single standard rock specimen, International J Rock Mechanics and Mining Sciences. 2008; 45: 1420-1438.
  14. Al-Obaidi SH, Khalaf FH. Acoustic logging methods in fractured and porous formations. J Geol Geophys. 2017; 6:2.
  15. Hofmann M, Al-Obaidi SH, Chang WJ. Evaluation Of Quantitative Criteria for Triassic Reservoirs in The South Mangyshlak Basin. Natural Science & Advanced Technology Education. 2023; 32: 7-24.
  16. Al-Obaidi SH, Guliaeva NI. Determination of flow and volumetric properties of core samples using laboratory NMR relaxometry. J Petrol Eng Technol. 2020; 1: 20-23.
  17. Xie QW, Liu XL, Zhang XP, Liu QS, Wang EZ. A review of test methods for uniaxial compressive strength of rocks: Theory, apparatus and data processing. J Rock Mechanics and Geotechnical Engineering. 2025; 17:1889-1905.
  18. Al-Obaidi SH, Khalaf FH, Alwan HH. Performance analysis of hydrocarbon wells based on the skin zone. Technium. 2012; 3: 50-56.
  19. Smirnov VI, Al-Obaidi S. Innovative methods of enhanced oil recovery. Oil Gas Res. 2008; 1: 1.
  20. Al-Obaidi SH, Chang W. Evaluation of the Quantitative Criteria of Triassic Carbonate Rocks Reservoirs. J Geol Geophys. 2023; 12: 1067.
  21. ASTM D3967-16. Standard test method for splitting tensile strength of intact rock core specimens. 2023; 4: 5.
  22. Al-Obaidi SH. Calculation Improvement of the Clay Content in the Hydrocarbon Formation Rocks. Oil Gas Res. 2017; 3: 130.
  23. Muralha J, Grasselli G, Tatone B, Blümel M, Chryssanthakis P, Yujing J. ISRM Suggested Method for Laboratory Determination of the Shear Strength of Rock Joints: Revised Version. Rock Mech Rock Eng. 2014; 47: 291-302.
  24. Khalaf FH, Al-Obaidi SH. The Effects of Hydro Confining Pressure on the Flow Properties of Sandstone and Carbonate Rocks. International J scientific & technology research. 2018; 7: 283-286.
  25. Cózar IR, Toro JJA, Maimí P, Otero F, González EV, Turon A. A novel methodology to measure the transverse Poisson’s ratio in the elastic and plastic regions for composite materials, Composites Part B: Engineering 2024; 272: 111098.
  26. Al-Obaidi S, Smirnov VI, Khalaf F. New technologies to improve the performance of high water cut wells equipped with ESP. Technium 2024; 3: 104-113.
  27. Chang WJ, Al-Obaidi SH, Patkin AA. The use of oil-soluble polymers to enhance oil recovery in hard to recover hydrocarbons reserves. International Research J Modernization in Engineering Technology and Science. 2021; 3: 982-987.
  28. AL-Obaidi SH, Chang WJ, Khalaf FH. Determination of the Upper Limit Up to Which the Linear Flow Law (Darcy’s Law) Can Be Applied. J Xidian University. 2012;15: 277-286.
  29. Wei-Qiang Xie, et al. A review of test methods for uniaxial compressive strength of rocks: Theory, apparatus and data processing, J Rock Mechanics and Geotechnical Engineering. 2024; 1674-7755.
  30. AL-Obaidi SH, Hofmann M, Smirnov VI, Khalaf FH, Alwan HH, et al. A Study of Compositions Relevant to Selective Water Isolation in Gas Wells. J Geol Geophys 2021; 10: 1000.
  31. Patkin A, Al-Obaidi SH. Influence of Temperature and Pressure of Incoming Oil- Containing Liquid from Field Wells on the Gas Separation Process. J Petroleum Engineering and Emerging Technology 2001; 3: 20-24.
  32. Al-Obaidi SH. High oil recovery using traditional water-flooding under compliance of the planned development mode. J Petroleum Engineering & Technology 2020; 6: 48-53.
  33. Kashnikov YUA, Ashihmin SG. Rock mechanics in petroleum industry. Moscow: Gornaya kniga. 2019; 496.
  34. Al-Obaidi SH, Khalaf FH. A New Approach for Enhancing Oil and Gas Recovery of the Hydrocarbon Fields with Low Permeability Reservoirs. Pet Petro Chem Eng J. 2023; 7: 000343.
  35. Hussain W, Luo M, Ali M, Khafaji HFA, Hussain I, Hussain M, et al. A Gamma-ray spectroscopy approach to evaluate clay mineral composition and depositional environment: A case study from the lower Goru Formation, Southern Indus Basin, Pakistan. J Applied Geophysics. 2024; 226: 105414.
  36. Al-Obaidi SH. Improve the efficiency of the study of complex reservoirs and hydrocarbon deposits – East Baghdad field. International J scientific & technology research. 2020; 5: 129-131.
  37. Galkin A, Al-Obaidi SH. Dependences of reservoir oil properties on surface oil. Jo Pet. Eng. Emerg. 2005; 5: 74-77.
  38. Al-Obaidi S, Hofmann M, Kamensky IV. A Study on How Permeability Heterogeneity Influences Initial Water Cut in Undersaturated Oil Formations. Natural Sciences and Advanced Technology Education 2024; 33: 347-361.
  39. Riyadi ZA, Olutoki JO, Hermana M, Latif AHA, Yogi IBS, Kadir SJA. Machine learning prediction of permeability distribution in the X field Malay Basin using elastic properties, Results in Engineering. 2024; 24: 103421.
  40. Al-Obaidi SH, Chang WJ, Hofmann M. Development of Oil Fields Using Science Artificial Intelligence and Machine Learning, Natural Sciences and Advanced Technology Education. 2023; 32: 3-4.
  41. Alexander G, Sudad AO, Viktor S. Numerical Simulation of Water Alternating Gas Injections (WAG) into Hydrocarbon Reservoirs: Factors Influencing Oil Recovery. J Geol Geophys. 2023; 12: 1147.
  42. Al-Obaidi SH, Khalaf FH. The Effect of anisotropy in formation permeability on the efficiency of cyclic water flooding. Int J Sci Technol Res. 2017; 6: 223-226.
  43. Nolte S, Fink R, Krooss BM, Littke R. Simultaneous determination of the effective stress coefficients for permeability and volumetric strain on a tight sandstone, J Natural Gas Science and Engineering. 2021; 95: 104186.
  44. Al-Obaidi SH, Galkin AP. The Assessment of Reservoir Pressure and Recoverable Hydrocarbon Reserves Using Production Well Research Data. J Geol Min. 2025; 2: 01-06.
  45. Hofmann M, Al-Obaidi S, Chang WJ. Evaluating The Efficacy of Water Alternating Gas Injection Technique in The Upper Jurassic Hydrocarbon Reservoirs, Considering Saturation Pressure. Natural Sciences and Advanced Technology Education 2024; 33: 3-4.
  46. Al-Obaidi SH, Wang JC, Hofmann M. Modelling the development of oil rim using water and gas injection. Nat Sci Adv Technol Edu 2022; 31: 279-292.
  47. Shijie S, Gao Y, Jia L. A Comparison of the Relationship between Dynamic and Static Rock Mechanical Parameters. Applied Sciences. 2024; 14: 4487.
  48. Al-Obaidi SH, Kamensky IP. Express study of rheological properties and group composition of oil and condensate using nuclear magnetic resonance–relaxometry. J Oil Gas Coal Technol. 2022; 4: 102.
  49. Chang WJ, Al-Obaidi SH, Patkin AA. The use of oil-soluble polymers to enhance oil recovery in hard to recover hydrocarbons reserves. Int Res J Mod Eng Technol Sci. 221; 3: 982-987.
  50. Al-Obaidi SH. Enhancing Hydrocarbon Field Recovery: Employing Multi-Stage Hydraulic Fracturing (MSHF) for the Development of Low-Permeability Deep-Lying Deposits. J Geol Earth Mar Sci. 2023; 1: 102.