A shale gas reservoir young's modulus determination method and device based on element logging

Abstract

The application discloses a shale gas reservoir Young's modulus determination method and device based on element logging. The method comprises the following steps: based on the mineral content, rock porosity, TOC content, water saturation, rock skeleton density and rock density of the shale gas reservoir rock determined according to element logging data, rock skeleton P-wave and S-wave velocities and rock P-wave and S-wave velocities are determined in sequence by using a volume model method; and rock skeleton Young's modulus and rock Young's modulus are determined respectively according to the rock skeleton P-wave and S-wave velocities, the rock P-wave and S-wave velocities, the rock skeleton density and the rock density. The method can utilize element logging data to quickly and accurately determine the dynamic Young's modulus of the rock skeleton and the rock.

Description

Technical Field

[0001] This invention relates to the field of comprehensive application of elemental logging data and calculation of logging rock mechanical parameters in the field of oil and gas exploration and development, and particularly to a method and apparatus for determining the Young's modulus of shale gas reservoirs based on elemental logging. Background Technology

[0002] The Young's modulus of rocks used in the logging industry mainly comes from two sources: core laboratory analysis and logging Young's modulus. The Young's modulus obtained through laboratory physical experiments is called the static Young's modulus, while the Young's modulus calculated through logging P-wave (S-wave) transit time is called the dynamic Young's modulus. Generally, the static Young's modulus is 0.6 to 0.8 times that of the dynamic Young's modulus.

[0003] The currently disclosed invention patents related to elemental logging and Young's modulus of rocks are mainly as follows:

[0004] The first invention patent, titled "A Quantitative Evaluation Method for Young's Modulus of Oil-Bearing Shale," with patent application number 201510409190.7 and application date of July 14, 2015, is authored by Li Junqian, Lu Shuangfang, and others from China University of Petroleum (East China). This invention patent establishes a method for fitting static Young's modulus related to mineral composition and confining pressure by performing correlation analysis between the mineral composition, confining pressure, and corresponding laboratory-obtained static Young's modulus of core samples. However, this method obtains only the static Young's modulus and does not consider rock porosity or fluid properties. The second invention patent, titled "A Method for Testing Young's Modulus of Rock Using Drill Cuttings," with patent application number 201510167358.8 and application date of April 10, 2015, is authored by Chen Chaowei, Xiang Degui, and others from the Drilling Engineering Technology Research Institute of China National Petroleum Corporation. This invention patent uses scanning electron microscopy to determine the particle size of the rock cuttings, selects a probe with an appropriate indentation angle based on different particle sizes, determines the maximum indentation depth of the probe, and then vertically indents the rock cuttings using the probe to obtain the load-displacement curve of the rock cuttings. Finally, it calculates the Young's modulus of the rock using the static Young's modulus theory formula. This method uses micron-level indentation technology to conduct physical experiments and obtains the static Young's modulus. The third invention patent is titled "A Method for Calculating Young's Modulus of Rock Based on Low-Field Nuclear Magnetic Resonance Response," with patent application number 201910286882.5 and application date of April 11, 2019. The inventors are Ge Xinmin, Zhao Jier, and others from China University of Petroleum (East China). This invention patent mainly utilizes correlation analysis between nuclear magnetic resonance logging parameters and laboratory static Young's modulus, employing a multiple regression method to establish a Young's modulus calculation model based on nuclear magnetic resonance response. The advantage of this method is that it does not require consideration of the physical meaning and theoretical calculation formula of Young's modulus; the static Young's modulus is directly obtained through data fitting. The disadvantages of this method are that the key parameters and fitting coefficients in the static Young's modulus fitting formula may vary depending on the region and the nuclear magnetic resonance logging instrument used, resulting in poor universality. Furthermore, current unconventional oil and gas development primarily employs horizontal well hydraulic fracturing, which does not involve coring, making it impossible to obtain nuclear magnetic resonance logging data. The fourth invention patent is titled "A Method and Device for Correcting Young's Modulus of Tight Sandstone," with invention patent application number 201811186028.3 and application date of October 11, 2018. The inventors are Ma Hui, Wang Feilong, and others from Sinopec North China Petroleum Engineering Co., Ltd. This invention patent corrects the dynamic Young's modulus parameters provided by well logging by using the silicon content, iron content, clay content, and porosity of tight sandstone collected by elemental logging. This makes the corrected dynamic Young's modulus more accurate and has a higher correlation with the static Young's modulus.This method assumes that the dynamic Young's modulus provided by well logging is inaccurate and corrects it. However, currently only well logging can provide dynamic Young's modulus, so this invention cannot prove that its corrected dynamic Young's modulus is more accurate or reasonable. In addition, this method relies on the dynamic Young's modulus provided by well logging, and it is invalid without relevant well logging data.

[0005] The first three invention patents mentioned above determine the static Young's modulus. Although the last one determines the dynamic Young's modulus, it has poor versatility, limited accuracy, and cannot be verified.

[0006] In addition, the article "Real-time Calculation Method of Rock Mechanical Parameters in Complex Strata—Taking Bachu Uplift in Tarim Basin as an Example," published by Luo Faqiang, Wang Zhizhan, and others from Sinopec in June 2020 in *Science, Technology and Engineering*, proposed two methods for calculating Poisson's ratio and Young's modulus. One method is the forward modeling method, which calculates minerals based on element logging and then calculates P-wave and S-wave velocities based on the minerals. However, this method does not consider reservoir fluids and porosity, resulting in low accuracy in calculating Poisson's ratio and Young's modulus. The other method is the direct method, which establishes a fitting formula for Poisson's ratio and Young's modulus by examining the correlation between single elements obtained from element logging and the logging Poisson's ratio and Young's modulus, and directly calculates Poisson's ratio and Young's modulus. This method is basically consistent with the principle mentioned in the invention patent "A Method for Obtaining Rock Poisson's Ratio Based on XRF Element Logging While Drilling" filed by Liang Bo et al. Meanwhile, because shale gas horizontal wells lack core data, and due to cost-saving or engineering factors, some horizontal wells are not logged, making it impossible to obtain rock mechanical parameters such as Poisson's ratio and Young's modulus. Therefore, logging cannot accurately evaluate the rock mechanical parameters such as Poisson's ratio and Young's modulus of shale gas reservoirs.

[0007] In summary, existing methods for calculating dynamic Young's modulus suffer from poor versatility and limited calculation accuracy. Summary of the Invention

[0008] In view of the above problems, the present invention is proposed to provide a method and apparatus for determining the Young's modulus of shale gas reservoirs based on elemental logging, which overcomes or at least partially solves the above problems. It can quickly and accurately determine the rock skeleton and dynamic Young's modulus of the rock using elemental logging data, and has a wide range of applications.

[0009] In a first aspect, embodiments of the present invention provide a method for determining the Young's modulus of shale gas reservoirs based on elemental logging, comprising:

[0010] Based on the mineral content, porosity, TOC content, water saturation, rock skeleton density, and rock density of the shale gas reservoir rocks determined by elemental logging data, the volumetric model method was used to sequentially determine the P- and S-wave velocities of the shale gas reservoir rock skeleton and the P- and S-wave velocities of the rocks.

[0011] Based on the longitudinal and transverse wave velocities of the rock skeleton, the longitudinal and transverse wave velocities of the rock, the density of the rock skeleton, and the density of the rock, the Young's modulus of the rock skeleton and the Young's modulus of the rock are determined respectively.

[0012] Secondly, embodiments of the present invention provide a device for determining the Young's modulus of shale gas reservoirs based on elemental logging, comprising:

[0013] The P-wave and S-wave velocity determination module is used to determine the P-wave and S-wave velocities of the shale gas reservoir rock skeleton and the P-wave and S-wave velocities of the rock skeleton sequentially using the volume model method, based on the mineral content, rock porosity, TOC content, water saturation, rock skeleton density and rock density determined from the elemental logging data.

[0014] The Young's modulus determination module is used to determine the Young's modulus of the rock skeleton and the Young's modulus of the rock based on the longitudinal and transverse wave velocities of the rock skeleton, the longitudinal and transverse wave velocities of the rock skeleton, the density of the rock skeleton, and the density of the rock.

[0015] Thirdly, embodiments of the present invention provide a computer storage medium storing computer-executable instructions, which, when executed by a processor, implement the above-mentioned method for determining the Young's modulus of shale gas reservoirs based on element logging.

[0016] Fourthly, this disclosure provides a server, including: a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, it implements the above-described method for determining the Young's modulus of shale gas reservoirs based on element logging.

[0017] The beneficial effects of the above-described technical solutions provided in the embodiments of the present invention include at least the following:

[0018] (1) The method for determining the Young's modulus of shale gas reservoirs based on elemental logging provided in this embodiment of the invention, based on the mineral content, porosity, TOC content, water saturation, rock skeleton density, and rock density of the shale gas reservoir rock determined by elemental logging data, sequentially determines the P-wave and S-wave velocities of the rock skeleton and the rock using a volumetric model method; then, based on the P-wave and S-wave velocities of the rock skeleton and the rock, the Young's modulus of the rock skeleton and the rock are determined. Utilizing elemental logging data, considering the properties of pores and fluids, the dynamic Young's modulus of the rock skeleton and the rock is determined quickly and accurately, providing data support for calculating parameters such as pore pressure, minimum horizontal principal stress, maximum horizontal principal stress, and fracture pressure of the shale gas reservoir using elemental logging, laying a solid foundation for logging rock mechanics evaluation.

[0019] (2) The method for determining the Young's modulus of shale gas reservoirs based on elemental logging provided in this embodiment of the invention uses only elemental logging data as the basis for calculation, without the need for core data and related logging data, so it has a wide range of applications.

[0020] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description, claims, and drawings.

[0021] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0022] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings:

[0023] Figure 1 This is a flowchart of the method for determining the Young's modulus of shale gas reservoirs based on elemental logging in Embodiment 1 of the present invention;

[0024] Figure 2 for Figure 1 The detailed implementation flowchart of step S11 is shown below;

[0025] Figure 3 for Figure 1 The detailed implementation flowchart of step S12 is shown below;

[0026] Figure 4 This is a flowchart illustrating the specific implementation of the method for determining the Young's modulus of shale gas reservoirs based on elemental logging in Embodiment 2 of the present invention.

[0027] Figure 5 This is a comparison chart of Young's modulus calculated by element logging and Young's modulus calculated by well logging in Embodiment 3 of the present invention;

[0028] Figure 6 This is a schematic diagram of the device for determining the Young's modulus of shale gas reservoirs based on elemental logging in an embodiment of the present invention. Detailed Implementation

[0029] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0030] It should be understood that the terminology used in this invention is merely for describing particular embodiments and is not intended to limit the invention. Furthermore, with respect to numerical ranges in this invention, it should be understood that each intermediate value between the upper and lower limits of the range is also specifically disclosed. Every smaller range between any stated value or intermediate value within a stated range, and any other stated value or intermediate value within said range, is also included in this invention. The upper and lower limits of these smaller ranges may be independently included or excluded from the range.

[0031] Unless otherwise stated, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains. While only preferred methods and materials have been described herein, any methods and materials similar or equivalent to those described herein may be used in the implementation or testing of this invention. All references to this specification are incorporated by way of citation to disclose and describe methods and / or materials associated with those references. In the event of any conflict with any incorporated reference, the content of this specification shall prevail.

[0032] This invention provides a method and apparatus for determining the Young's modulus of shale gas reservoirs based on elemental logging, which can quickly and accurately determine the rock skeleton and the dynamic Young's modulus of the rock using elemental logging data.

[0033] Example 1

[0034] Embodiment 1 of this invention provides a method for determining the Young's modulus of shale gas reservoirs based on elemental logging, the process of which is as follows: Figure 1 As shown, it includes the following steps:

[0035] Step S11: Based on the mineral content, rock porosity, TOC content, water saturation, rock skeleton density, and rock density determined from the elemental logging data, the P- and S-wave velocities of the shale gas reservoir rock skeleton and the P- and S-wave velocities of the rock are determined sequentially using the volumetric model method.

[0036] Shale gas reservoirs consist of a rock framework, total organic carbon (TOC), shale gas, and formation water. Assuming the shale gas reservoir is a homogeneous body, see [link to relevant documentation]. Figure 2 As shown, the rock skeleton and P- and S-wave velocities of the shale gas reservoir are determined through the following steps:

[0037] Step S111: Obtain the mineral content, rock porosity, TOC content, water saturation, rock skeleton density, and rock density of the shale gas reservoir rocks determined based on elemental logging data.

[0038] Specifically, the above parameters are obtained in the following way:

[0039] (1) Calculate the mineral content in shale gas reservoir rocks according to the method of elemental logging inversion of sedimentary rock minerals, including the content of at least one mineral among quartz, potassium feldspar, sodium feldspar, analcime, kaolinite, montmorillonite, illite, mica, calcite, dolomite, gypsum and pyrite.

[0040] Furthermore, the content here refers to the percentage of minerals by mass relative to the rock skeleton.

[0041] Mineral content can be a decimal, that is, a number less than 1, and the sum of the contents of all minerals equals 1 (the unit of mineral content can also be %, and the sum of the contents of all minerals equals 100). However, in actual calculations, it is almost impossible to calculate all the mineral components in the rock. Therefore, minerals with extremely low content and less important content do not need to be calculated. This makes the calculated sum of the contents of all minerals equal to 1 an approximation.

[0042] (2) Calculate rock porosity using the method of calculating shale gas reservoir porosity based on element logging.

[0043] (3) Calculate the total organic carbon (TOC) content, water saturation, rock skeleton density and rock density based on the method of calculating shale gas reservoir density by element logging.

[0044] The methods for calculating various parameters based on element logging described above can all be implemented using existing technologies, so they will not be described in detail here.

[0045] Step S112: Determine the P- and S-wave velocities of the rock skeleton using the volumetric model method based on the mineral content, rock skeleton density, and predetermined mineral P- and S-wave velocities.

[0046] Based on the mineral content, rock skeleton density, and predetermined mineral P-wave velocities, the following formula (1) is used to determine the rock skeleton P-wave velocity V. ps-el :

[0047]

[0048] Based on the mineral content, rock skeleton density, and predetermined mineral transverse wave velocities, the rock skeleton transverse wave velocity V is determined using the following formula (2). ss-el :

[0049]

[0050] In formulas (1) and (2), ρ s-el For the density of the rock skeleton, ω i-el ρ i-el V ps-i-el and V ss-i-elThese represent the content, density, P-wave velocity, and S-wave velocity of the i-th mineral, respectively, where i is the index of the mineral type (i = 1, 2, ..., n), and n is the number of mineral types, with a mineral content less than 1. Furthermore, the units for the rock skeleton density and mineral density are not specifically limited, as long as they are consistent; the units for the P-wave and S-wave velocities of the minerals are also not specifically limited, as long as the units for the P-wave and S-wave velocities of various minerals are consistent, and the units of the calculated rock skeleton P-wave and S-wave velocities are consistent with the units of the mineral P-wave and S-wave velocities; optionally, the unit for mineral content can also be %, which can be converted to a decimal by dividing by 100 during the calculation.

[0051] Step S113: Based on the P and S wave velocities of the rock skeleton, rock porosity, rock density, TOC content and water saturation, as well as the pre-determined P and S wave velocities and TOC density of shale gas, TOC and formation water, the volumetric model method is used to determine the P and S wave velocities of the shale gas reservoir rocks.

[0052] Specifically, the longitudinal wave velocity V of the shale gas reservoir rock is determined using the following formula (3). pb-el The shear wave velocity V of shale gas reservoir rocks is determined using the following formula (4). sb-el :

[0053]

[0054]

[0055] In formulas (3) and (4), V ps-el and V ss-el These represent the longitudinal and transverse wave velocities of the rock skeleton, V. p-fw and V s-fw These are the longitudinal wave velocity and transverse wave velocity of formation water, respectively, V pg-el and V sg-el These are the P-wave velocity and S-wave velocity of shale gas, respectively, V pb-TOC-ef and V sb-TOC-ef These are the TOC longitudinal wave velocity and transverse wave velocity, respectively, ρ b-el and ρ TOC These are rock density and TOC density, respectively. el φ el and S w-el These are TOC content, rock porosity, and water saturation, respectively; TOC el φ el and S w-el All values ​​are less than 1. This is optional, and the unit can also be %. In the calculation process, simply divide by 100 to convert it to a decimal.

[0056] Furthermore, the units for rock density and TOC density are not specifically limited, as long as they are consistent; V p-fw V s-fwV pg-el V sg-el V pb-TOC-ef and V sb-TOC-ef The units are not specifically limited, as long as they are consistent. The units of the calculated rock longitudinal and transverse wave velocities are consistent with those of the former.

[0057] Step S12: Determine the Young's modulus of the rock skeleton and the Young's modulus of the rock skeleton based on the longitudinal and transverse wave velocities of the rock skeleton, the density of the rock skeleton, and the density of the rock.

[0058] See Figure 3 As shown, the specific steps include the following:

[0059] Step S121: Determine the time difference between P and S waves of the rock skeleton based on the P and S waves velocity of the rock skeleton.

[0060] The reciprocal of the speed can be used to determine the corresponding time difference.

[0061] Step S122: Determine the Young's modulus of the rock skeleton and the Young's modulus of the rock based on the rock skeleton density and P-wave transit time, and the rock density and P-wave transit time.

[0062] Specifically, the Young's modulus of the rock skeleton and the Young's modulus of the rock are determined by the following formulas (5) and (6):

[0063]

[0064]

[0065] In formulas (5) and (6), E s-el and E b-el These are the Young's modulus of the rock framework and the Young's modulus of the rock, ρ. s-el and ρ b-el These are the rock skeleton density and rock density, respectively, DT ps-el DT ss-el DT pb-el and DT sb-el These are the longitudinal wave time difference and transverse wave time difference of the rock skeleton, and the longitudinal wave time difference and transverse wave time difference of the rock itself.

[0066] Furthermore, regarding the units of the parameters in the above formulas (5) and (6), the unit of Young's modulus can be GPa, and the unit of density can be g / cm³. 3 The unit of the P-wave and S-wave time difference is μs / m. At this time, since the above formulas (5) and (6) need to be transformed into:

[0067]

[0068]

[0069] Optionally, other units can be used for the above-mentioned parameters. The specific form of the formula needs to be further multiplied or divided by a certain constant based on formulas (5) and (6) to achieve matching of the units of various parameters.

[0070] The method for determining the Young's modulus of shale gas reservoirs based on elemental logging provided in Embodiment 1 of this invention utilizes elemental logging data, considering the properties of pores and fluids, to quickly and accurately determine the rock skeleton and the dynamic Young's modulus of the rock. This provides data support for calculating parameters such as pore pressure, minimum horizontal principal stress, maximum horizontal principal stress, and fracture pressure of shale gas reservoirs using elemental logging, laying a solid foundation for logging-based rock mechanics evaluation. This method uses only elemental logging data as the basis for calculation, without requiring core data or related logging data, thus having a wide range of applications.

[0071] Example 2

[0072] Embodiment 2 of the present invention provides a specific implementation flow of a method for determining the Young's modulus of shale gas reservoirs based on elemental logging, as follows: Figure 4 As shown, it includes the following steps:

[0073] Step S41: Obtain the mineral content, rock porosity, TOC content, water saturation, rock skeleton density, and rock density of the shale gas reservoir rocks determined based on elemental logging data.

[0074] The mineral content includes the content of quartz, potassium feldspar, sodium feldspar, analcime, kaolinite, montmorillonite, illite, mica, calcite, dolomite, gypsum, and pyrite, specifically the mass percentage of the corresponding mineral in the rock framework, in percentages (%).

[0075] Step S42: Determine the P and S wave velocities of the rock skeleton using the volumetric model method based on the mineral content, density, rock skeleton density, and predetermined P and S wave velocities of the minerals.

[0076] V ps-el =(ω quartz-el ×ρ s-el / ρ quartz-el ×V ps-quartz-el +ω pf-el ×ρ s-el / ρ pf-el ×V ps-pf-el +ω albite-el ×ρ s-el / ρ albite-el ×V ps-albite-el +ω analcite-el ×ρ s-el / ρ analcite-el ×V ps-analcite-el +ω kaolinite-el ×ρ s-el / ρkaolinite-el ×V ps-kaolinite-el +oh montmorillonite-el ×r s-el / r montmorillonite-el ×V ps-montmorillonite-el +oh illite-el ×r s-el / r illite-el ×V ps-illite-el +oh mica-el ×r s-el / r mica-el ×V ps-mica-el +oh calcite-el ×r s-el / r calcite-el ×V ps-calcite-el +oh dolomite-el ×r s-el / r dolomite-el ×V ps-dolomite-el +oh plaster-el ×r s-el / r plaster-el ×V ps-plaster-el +oh pyrite-el ×r s-el / r pyrite-el ×V ps-pyrite-el ) / 100 (9)

[0077] V ss-el =(ω quartz-el ×r s-el / r quartz-el ×V ss-quartz-el +oh pf-el ×r s-el / r pf-el ×V ss-pf-el +oh albite-el ×r s-el / r albite-el ×V ss-albite-el +oh analcite-el ×r s-el / r analcite-el ×V ss-analcite-el +oh kaolinite-el ×r s-el / r kaolinite-el ×V ss-kaolinite-el +oh montmorillonite-el ×r s-el / r montmorillonite-el ×V ss-montmorillonite-el +oh illite-el ×r s-el / r illite-el ×V ss-illite-el +oh mica-el ×r s-el / r mica-el ×V ss-mica-el +oh calcite-el ×rs-el / ρ calcite-el ×V ss-calcite-el +ω dolomite-el ×ρ s-el / ρ dolomite-el ×V ss-dolomite-el +ω plaster-el ×ρ s-el / ρ plaster-el ×V ss-plaster-el +ω pyrite-el ×ρ s-el / ρ pyrite-el ×V ss-pyrite-el ) / 100 (10)

[0078] In formulas (9) and (10), V ps-el V represents the longitudinal wave velocity of the rock skeleton. ss-el For the transverse wave velocity of the rock skeleton, ρ s-el ρ is the density of the rock skeleton. quartz-el ρ is the density of quartz. pf-el The density of potassium feldspar, ρ albite-el The density of albite is ρ. analcite-el ρ is the density of anabolic zeolite. kaolinite-el ρ is the density of kaolinite. montmorillonite-el ρ is the density of montmorillonite. illite-el The density of illite, ρ mica-el ρ is the density of mica. calcite-el ρ is the density of calcite. dolomite-el ρ is the density of dolomite. plaster-el ρ is the density of gypsum. pyrite-el The density of pyrite is 2.5, the density of potassium feldspar is 2.57, the density of sodium feldspar is 2.62, the density of analcime is 2.6, the density of kaolinite is 2.63, the density of montmorillonite is 2.6, the density of illite is 2.8, the density of mica is 3.05, the density of calcite is 2.71, the density of dolomite is 2.87, the density of gypsum is 2.4, and the density of pyrite is 5.2; ω quartz-el ω represents the percentage of quartz mass. pf-el ω represents the mass percentage of potassium feldspar. albite-el ω represents the mass percentage of albite. analcite-el ω represents the mass percentage of anabolic zeolite. kaolinite-el ω represents the percentage of kaolinite by mass. montmorillonite-el ω represents the percentage of montmorillonite by mass. illite-el ω represents the percentage of illite by mass. mica-el ω represents the percentage of mica mass. calcite-el Let be the mass percentage of calcite, and be ω. dolomite-el Dolomite mass percentage, ω plaster-el ω represents the percentage of gypsum by mass. pyrite-elV represents the percentage of pyrite by mass. ps-quartz-el V is the longitudinal wave velocity of the quartz skeleton. ps-pf-el V represents the longitudinal wave velocity of the potassium feldspar framework. ps-albite-el V represents the longitudinal wave velocity of the albite framework. ps-analcite-el V represents the longitudinal wave velocity of the anabolic zeolite framework. ps-kaolinite-el V represents the longitudinal wave velocity of the kaolinite skeleton. ps-montmorillonite-el V represents the longitudinal wave velocity of the montmorillonite skeleton. ps-illite-el V represents the longitudinal wave velocity of the illite framework. ps-mica-el V is the longitudinal wave velocity of the mica skeleton. ps-calcite-el Let V be the longitudinal wave velocity of the calcite framework. ps-dolomite-el Longitudinal wave velocity of dolomite skeleton, V ps-plaster-el V is the longitudinal wave velocity of the plaster skeleton. ps-pyrite-el V represents the longitudinal wave velocity of the pyrite skeleton. ss-quartz-el V is the transverse wave velocity of the quartz skeleton. ss-pf-el V represents the transverse wave velocity of the potassium feldspar framework. ss-albite-el V represents the transverse wave velocity of the albite framework. ss-analcite-el V represents the transverse wave velocity of the anabolic zeolite framework. ss-kaolinite-el V represents the transverse wave velocity of the kaolinite skeleton. ss-montmorillonite-el V represents the transverse wave velocity of the montmorillonite skeleton. ss-illite-el V is the shear wave velocity of the illite skeleton. ss-mica-el V is the shear wave velocity of the mica skeleton. ss-calcite-el Let V be the transverse wave velocity of the calcite framework. ss-dolomite-el Transverse wave velocity of dolomite skeleton, V ss-plaster-el V is the transverse wave velocity of the plaster skeleton. ss-pyrite-el The transverse wave velocity is the pyrite skeleton velocity; the longitudinal (transverse) wave velocity is in m / s, the mineral content is in %, and the density is in g / cm³. 3 .

[0079] Step S43: Based on the P and S wave velocities of the rock skeleton, rock porosity, rock density, TOC content and water saturation, as well as the pre-determined P and S wave velocities and TOC density of shale gas, TOC and formation water, the P and S wave velocities of the shale gas reservoir rock are determined using the volumetric model method.

[0080] V pb-el =V ps-el ×(1-ρ b-el ×TOC el / ρ TOC / 100-φ el / 100)+S w-el ×φ el / 10000×V p-fw +(1-S w-el / 100)×φ el / 100×Vpg-el +ρ b-el ×TOC el / ρ TOC / 100×V pb-TOC-ef (11)

[0081] V sb-el =V ss-el ×(1-ρ b-el ×TOC el / ρ TOC / 100-φ el / 100)+S w-el ×φ el / 10000×V s-fw +

[0082] (1-S w-el / 100)×φ el / 100×V sg-el +ρ b-el ×TOC el / ρ TOC / 100×V sb-TOC-ef (12)

[0083] Longitudinal (transverse) wave velocity is expressed in m / s; porosity, TOC, and water saturation are expressed in %, and density is expressed in g / cm³. 3 , where V pb-el V is the longitudinal wave velocity of the rock. sb-el TOC represents the transverse wave velocity of the rock. el The TOC content of shale gas reservoirs calculated from elemental logging, ρ TOC For TOC density, φ el V is the porosity of shale gas reservoir rocks calculated using elemental logging. p-fw S represents the longitudinal wave velocity of formation water. w-el V is the water saturation of the shale gas reservoir calculated from elemental logging. pg-el V represents the longitudinal wave velocity of shale gas. pb-TOC-ef V represents the TOC P-wave velocity of the shale gas reservoir. s-fw V is the shear wave velocity of the formation water. sg-el V represents the shear wave velocity of shale gas. sb-TOC-ef The TOC shear wave velocity is the velocity of the shale gas reservoir.

[0084] Step S44: Calculate the P-wave (S-wave) transit time of the shale gas reservoir rock skeleton and the P-wave (S-wave) transit time of the rock skeleton based on the P-wave (S-wave) velocity and the P-wave (S-wave) velocity of the rock.

[0085] DT ps-el =1000000 / V ps-el (13)

[0086] DT ss-el =1000000 / V ss-el (14)

[0087] DT pb-el =1000000 / V pb-el (15)

[0088] DT sb-el =1000000 / V sb-el (16)

[0089] The unit for P-wave (S-wave) velocity is m / s, and the unit for P-wave (S-wave) time difference is μs / m, where DT ps-el P-wave transit time of shale gas reservoir rock skeleton calculated for elemental logging, DT ss-el Shear wave time difference of shale gas reservoir rock skeleton calculated for elemental logging, DT pb-el P-wave transit time of shale gas reservoir rocks calculated for elemental logging, DT sb-el Shear wave time difference of shale gas reservoir rock calculated for elemental logging.

[0090] Step S45: Based on the rock skeleton density, rock density, shale gas reservoir rock skeleton P-wave (S-wave) transit time, and rock P-wave (S-wave) transit time calculated from elemental logging, use the Young's modulus acoustic calculation formula to calculate the Young's modulus of the shale gas reservoir rock skeleton and the Young's modulus of the rock.

[0091] E s-el =ρ s-el / (DT ss-el ×DT ss-el )×(3×DT ss-el ×DT ss-el -4×DT ps-el ×DT ps-el )×

[0092] 1000000 / (DT ss-el ×DT ss-el -DT ps-el ×DT ps-el ) (17)

[0094] E b-el =ρ b-el / (DT sb-el ×DT sb-el )×(3×DT sb-el ×DT sb-el -4×DT pb-el ×DT pb-el )×

[0095] 1000000 / (DT sb-el ×DT sb-el-DT pb-el ×DT pb-el ) (18)

[0097] Young's modulus is measured in GPa, and density is measured in g / cm³. 3 E s-el E represents the Young's modulus of the shale gas reservoir rock framework. b-el ρ is the Young's modulus of shale gas reservoir rocks. s-el ρ is the density of the rock skeleton. b-el This represents the density of the rock.

[0098] Example 3

[0099] This invention provides a specific application of a method for determining the Young's modulus of shale gas reservoirs based on elemental logging. Taking a horizontal shale gas well B in the xx area of ​​the xx basin as an example, the mineral contents of quartz, calcite, dolomite, feldspar, illite, kaolinite, montmorillonite, pyrite, etc., obtained from elemental logging inversion, along with the TOC content, water saturation, and porosity calculated from elemental logging, are substituted into the aforementioned formula for calculating the Young's modulus of rock from elemental logging. Within the well section from 4000m to 5840m, 773 elemental samples were analyzed. The minimum calculated Young's modulus was 37.04 GPa, the maximum was 73.79 GPa, and the average was 44.71 GPa. The minimum logging Young's modulus corresponding to the elemental analysis well depth was 36.32 GPa, the maximum was 69.21 GPa, and the average was 46.92 GPa. Comparing the Young's modulus calculated from elemental logging with that from well logging, the maximum positive error was 39.67%, with an average of 10.33%, while the maximum negative error was -38.61%, with an average of -9.60%. Overall, the errors remained largely within ±10% (see...). Figure 5 (See Table 1). It can be seen that the Young's modulus calculated through the embodiments of the present invention and the error of the Young's modulus in well logging meet the production requirements, and the overall trend is consistent.

[0100] Table 1B shows the Young's modulus calculated from well element logging and the statistics of Young's modulus from well logging.

[0101]

[0102]

[0103] Note: Table 1 is only an illustration and does not fully show the details of the 773 sample points.

[0104] Example 4

[0105] Embodiment 4 of this invention provides a method for calculating the Young's modulus of shale gas reservoirs based on elemental logging. Taking a horizontal shale gas well in the xx area of ​​the xx basin as an example, the method includes the following steps:

[0106] (1) Calculate the content of minerals such as quartz, potassium feldspar, sodium feldspar, analcime, kaolinite, montmorillonite, illite, mica, calcite, dolomite, gypsum, and pyrite based on the method of inverting sedimentary rock minerals by element logging;

[0107] Within the well section with a depth of 4000–5840 m, elemental logging measurements and analysis were performed on the samples. Taking the well depth of 4168 m as an example, based on the inversion calculation of the content of each element, the mass fractions of the above minerals are as follows: quartz 55.27%, calcite 22.16%, dolomite 7.45%, clay 9.59%, pyrite 3.73%, and albite 1.80%. The clay is further divided into illite 0.00%, montmorillonite 0.69%, kaolinite 0.00%, and mica 8.90%.

[0108] (2) The porosity of the shale gas reservoir was calculated using the method of calculating the porosity of the shale gas reservoir based on element logging. The porosity of the well at a depth of 4168m was 5.57%.

[0109] (3) The total organic carbon (TOC) content, water saturation, framework density, and rock density of the shale gas reservoir were calculated using the method of calculating shale gas reservoir density by element logging. The TOC content of the well at a depth of 4168m was 4.78%, the water saturation was 15.62%, the framework density was 2.67 g / cm3, and the rock density was 2.51 g / cm3.

[0110] (4) Based on the mineral content, framework density, rock density and P (S) wave velocity corresponding to single minerals such as quartz, potassium feldspar, sodium feldspar, zeolite, kaolinite, montmorillonite, illite, mica, calcite, dolomite, gypsum, and pyrite calculated by element logging, the P (S) wave velocity of the shale gas reservoir framework was calculated using a volumetric model. The P (S) wave velocity of the framework at a well depth of 4168m was 5882.35m / s and the S (S) wave velocity of the framework was 3412.39m / s.

[0111] (5) Based on the longitudinal (transverse) wave velocity, porosity, TOC content, water saturation, and the corresponding densities of TOC, shale gas, and formation water, the longitudinal (transverse) wave velocity of the shale gas reservoir rock was determined. At a well depth of 4168m, the longitudinal wave velocity of the rock was 4511.82m / s and the transverse wave velocity was 2769.09m / s.

[0112] (6) Based on the longitudinal (transverse) wave velocity of the shale gas reservoir skeleton and the longitudinal (transverse) wave velocity of the rock, the longitudinal (transverse) wave transit time of the shale gas reservoir skeleton and the longitudinal (transverse) wave transit time of the rock were calculated. At a well depth of 4168m, the longitudinal wave transit time of the skeleton was 170.00μs / m, the transverse wave transit time of the skeleton was 293.05μs / m, the longitudinal wave transit time of the rock was 221.64μs / m, and the transverse wave transit time of the rock was 361.13μs / m.

[0113] (7) Based on the calculation formulas for the longitudinal (horizontal) transit time of the shale gas reservoir skeleton, the longitudinal (horizontal) transit time of the rock, and the Young's modulus sonic parameters calculated from elemental logging (Well Logging Principles and Comprehensive Interpretation / Hong Youmi, ed. — Dongying: China University of Petroleum Press, 1998 (reprinted in February 2007)), the Young's modulus of the shale gas reservoir skeleton and the Young's modulus of the rock are calculated. According to the calculation formulas, the Young's modulus E of the skeleton at a well depth of 4168m is... s-el =2.67 / (293.05×293.05)×(3×293.05×293.05-4×170.00×170.00)×1000000 / (293.05×293.05-170.00×170.00)=77.47GPa, Young's modulus of rock E s-el =2.67 / (361.13×361.13)×(3×361.13×361.13-4×221.64×221.64)×1000000 / (361.13×361.13-221.64×221.64)=46.04GPa.

[0114] In this embodiment, at the same depth and location of the shale gas well, the logging Young's modulus is 43.11 GPa, and the error between the elemental logging-calculated Young's modulus and the logging Young's modulus is 6.80%.

[0115] Example 5

[0116] Embodiment 5 of the present invention provides a method for calculating the Young's modulus of shale gas reservoirs based on elemental logging. Taking a horizontal shale gas well in the xx area of ​​the xx basin as an example, the method includes the following steps:

[0117] (1) Calculate the content of minerals such as quartz, potassium feldspar, sodium feldspar, analcime, kaolinite, montmorillonite, illite, mica, calcite, dolomite, gypsum, and pyrite based on the method of inverting sedimentary rock minerals by element logging;

[0118] Within the well section with a depth of 4000–5840 m, elemental logging measurements and analysis were performed on the samples. Taking the well depth of 4388 m as an example, based on the inversion calculation of the content of each element, the mass fractions of the above minerals are as follows: quartz 55.31%, calcite 19.76%, dolomite 6.99%, clay 12.13%, pyrite 4.09%, and albite 1.72%. Among them, the clay is divided into illite 1.80%, montmorillonite 0.75%, kaolinite 0.00%, and mica 9.58%.

[0119] (2) The porosity of the shale gas reservoir was calculated using the method of calculating the porosity of the shale gas reservoir based on element logging. The porosity of the well at a depth of 4388m was 5.60%.

[0120] (3) The total organic carbon (TOC) content, water saturation, framework density, and rock density of the shale gas reservoir were calculated using the method of calculating shale gas reservoir density by element logging. The TOC content of the well at a depth of 4388m was 5.30%, the water saturation was 28.62%, the framework density was 2.68 g / cm3, and the rock density was 2.50 g / cm3.

[0121] (4) Based on the mineral content, framework density, rock density and P (S) wave velocity corresponding to single minerals such as quartz, potassium feldspar, sodium feldspar, zeolite, kaolinite, montmorillonite, illite, mica, calcite, dolomite, gypsum, and pyrite calculated by element logging, the P (S) wave velocity of the shale gas reservoir framework was calculated using a volumetric model. The P (S) wave velocity of the framework at a well depth of 4388m was 5765.35m / s and the S (S) wave velocity of the framework was 3301.2m / s.

[0122] (5) Based on the longitudinal (transverse) wave velocity, porosity, TOC, water saturation and TOC content of the shale gas reservoir skeleton, the corresponding densities of shale gas and formation water, and the longitudinal (transverse) wave velocity, the longitudinal (transverse) wave velocity of the shale gas reservoir rock was determined. At a well depth of 4388m, the longitudinal wave velocity of the rock was 4494.18m / s and the transverse wave velocity of the rock was 2775.08m / s.

[0123] (6) Based on the longitudinal (transverse) wave velocity of the shale gas reservoir skeleton and the longitudinal (transverse) wave velocity of the rock, the longitudinal (transverse) wave transit time of the shale gas reservoir skeleton and the longitudinal (transverse) wave transit time of the rock were calculated. At a well depth of 4388m, the longitudinal wave transit time of the skeleton was 173.45μs / m, the transverse wave transit time of the skeleton was 302.92μs / m, the longitudinal wave transit time of the rock was 222.51μs / m, and the transverse wave transit time of the rock was 360.35μs / m.

[0124] (7) Based on the sonic calculation formulas for the longitudinal (horizontal) transit time of the shale gas reservoir framework, the longitudinal (horizontal) transit time of the rock, and the Young's modulus of the shale gas reservoir calculated from elemental logging (Well Logging Principles and Comprehensive Interpretation / Hong Youmi, ed. — Dongying: China University of Petroleum Press, 1998 (reprinted in February 2007)), the Young's modulus of the shale gas reservoir framework and the Young's modulus of the rock are calculated. According to the calculation formulas, the Young's modulus E of the framework at a well depth of 4388m is... s-el =2.68 / (302.92×302.92)×(3×302.92×302.92-4×173.45×173.45)×1000000 / (302.92×302.92-173.45×173.45)=73.27GPa, Young's modulus of rock E s-el =2.50 / (360.35×360.35)×(3×360.35×360.35-4×222.51×222.51)×1000000 / (360.35×360.35-222.51×222.51)=45.88GPa.

[0125] In this embodiment, at the same depth and location of the shale gas well, the logging Young's modulus is 44.58 GPa, and the error between the elemental logging-calculated Young's modulus and the logging Young's modulus is 2.92%.

[0126] Based on the inventive concept of this invention, embodiments of this invention also provide a device for determining the Young's modulus of shale gas reservoirs based on elemental logging. The structure of this device is as follows: Figure 6 As shown, it includes:

[0127] The P-wave and S-wave velocity determination module 61 is used to determine the P-wave and S-wave velocities of the shale gas reservoir rock skeleton and the P-wave and S-wave velocities of the rock skeleton in sequence using the volume model method based on the mineral content, rock porosity, TOC content, water saturation, rock skeleton density and rock density in the shale gas reservoir rock determined by elemental logging data.

[0128] The Young's modulus determination module 62 is used to determine the Young's modulus of the rock skeleton and the Young's modulus of the rock skeleton based on the longitudinal and transverse wave velocities of the rock skeleton, the longitudinal and transverse wave velocities of the rock skeleton, the density of the rock skeleton, and the density of the rock.

[0129] In some embodiments, the P-wave and S-wave velocity determination module 61 is specifically used for:

[0130] Obtain the mineral content, porosity, TOC content, water saturation, rock skeleton density, and rock density of the shale gas reservoir rocks based on elemental logging data; determine the P- and S-wave velocities of the rock skeleton using a volumetric model method based on the mineral content, rock skeleton density, and pre-determined P- and S-wave velocities of the minerals; determine the P- and S-wave velocities of the shale gas reservoir rocks using a volumetric model method based on the P- and S-wave velocities of the rock skeleton, rock porosity, rock density, TOC content, water saturation, and pre-determined P- and S-wave velocities and TOC density of shale gas, TOC, and formation water.

[0131] In some embodiments, the P-wave and S-wave velocity determination module 61 is specifically used for:

[0132] Based on the mineral content, rock skeleton density, and predetermined mineral P-wave velocities, the rock skeleton P-wave velocity V is determined using the following formula (1). ps-el :

[0133]

[0134] Based on the mineral content, rock skeleton density, and predetermined mineral transverse wave velocities, the rock skeleton transverse wave velocity V is determined using the following formula (2). ss-el :

[0135]

[0136] In formulas (1) and (2), ρ s-el For the density of the rock skeleton, ω i-el ρ i-el V ps-i-el and V ss-i-el These are the content, density, P-wave velocity, and S-wave velocity of the i-th mineral, respectively, where i is the index of the mineral type, i = 1, 2, ..., n, and n is the number of mineral types, with a mineral content less than 1.

[0137] In some embodiments, the P-wave and S-wave velocity determination module 61 is specifically used for:

[0138] Based on the P and S wave velocities of the rock skeleton, rock porosity, rock density, TOC content, and water saturation, and the pre-determined P and S wave velocities and TOC density of shale gas, TOC, and formation water, the P wave velocity V of the shale gas reservoir rock is determined using the following formula (3). pb-el The shear wave velocity V of shale gas reservoir rocks is determined using the following formula (4). sb-el :

[0139]

[0140]

[0141] In formulas (3) and (4), V ps-el and V ss-el These represent the longitudinal and transverse wave velocities of the rock skeleton, V. p-fw and V s-fw These are the longitudinal wave velocity and transverse wave velocity of formation water, respectively, V pg-el and V sg-el These are the P-wave velocity and S-wave velocity of shale gas, respectively, V pb-TOC-ef and V sb-TOC-ef These are the TOC longitudinal wave velocity and transverse wave velocity, respectively, ρ b-el and ρ TOC These are rock density and TOC density, respectively. el φ el and S w-el These are TOC content, rock porosity, and water saturation, respectively. el φ el and S w-el All are less than 1.

[0142] In some embodiments, the Young's modulus determination module 62 is specifically used for:

[0143] The P and S wave transit times of the rock skeleton are determined based on the P and S wave velocities of the rock skeleton, and the P and S wave transit times of the rock are determined based on the P and S wave velocities of the rock skeleton. The Young's modulus of the rock skeleton and the Young's modulus of the rock are determined based on the rock skeleton density and the P and S wave transit times, and the rock density and the P and S wave transit times, respectively.

[0144] In some embodiments, the Young's modulus determination module 62 is specifically used for:

[0145] Based on the rock skeleton density and P-wave transit time, and the rock density and P-wave transit time, the Young's modulus of the rock skeleton and the Young's modulus of the rock are determined by the following formulas (5) and (6), respectively:

[0146]

[0147]

[0148] In formulas (5) and (6), E s-el and E b-el These are the Young's modulus of the rock framework and the Young's modulus of the rock, ρ. s-el and ρ b-el These are the rock skeleton density and rock density, respectively, DT ps-el DT ss-el DT pb-el and DT sb-el These are the longitudinal wave time difference and transverse wave time difference of the rock skeleton, and the longitudinal wave time difference and transverse wave time difference of the rock itself.

[0149] Regarding the apparatus in the above embodiments, the specific manner in which each module performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated upon here.

[0150] Based on the inventive concept of the present invention, the embodiments of the present invention also provide a computer storage medium storing computer-executable instructions, which, when executed by a processor, implement the above-mentioned method for determining the Young's modulus of shale gas reservoirs based on element logging.

[0151] Based on the inventive concept of this invention, this embodiment of the invention also provides a server, including: a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, it implements the above-mentioned method for determining the Young's modulus of shale gas reservoirs based on element logging.

[0152] Unless otherwise specifically stated, terms such as processing, calculation, operation, determination, display, etc., may refer to the actions and / or processes of one or more processing or computing systems or similar devices that represent the manipulation and conversion of data representing physical (e.g., electronic) quantities within the registers or memory of the processing system into other data similarly representing physical quantities within the memory, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals can be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description can be represented by voltage, current, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0153] It should be understood that the specific order or hierarchy of steps in the disclosed process is an example of an exemplary method. Based on design preferences, it should be understood that the specific order or hierarchy of steps in the process may be rearranged without departing from the scope of this disclosure. The appended method claims provide elements of various steps in an exemplary order and are not intended to limit the scope to the specific order or hierarchy described.

[0154] In the detailed description above, various features are combined together in a single embodiment to simplify this disclosure. This approach to disclosure should not be construed as reflecting an intention that embodiments of the claimed subject matter require more features than are explicitly stated in each claim. Rather, as reflected in the appended claims, the invention is presented with fewer features than all of the features in a single disclosed embodiment. Therefore, the appended claims are hereby explicitly incorporated into the detailed description, with each claim representing a separate preferred embodiment of the invention.

[0155] Those skilled in the art will also understand that the various illustrative logic blocks, modules, circuits, and algorithm steps described in conjunction with the embodiments herein can be implemented as electronic hardware, computer software, or a combination thereof. To clearly illustrate the interchangeability between hardware and software, the various illustrative components, blocks, modules, circuits, and steps described above are generally described in terms of their functionality. Whether such functionality is implemented as hardware or software depends on the specific application and the design constraints imposed on the overall system. Those skilled in the art can implement the described functionality in alternative ways for each specific application; however, such implementation decisions should not be construed as departing from the scope of this disclosure.

[0156] The steps of the methods or algorithms described in conjunction with the embodiments herein can be directly embodied in hardware, software modules executed by a processor, or a combination thereof. The software modules can reside in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disks, removable disks, CD-ROMs, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor, enabling the processor to read information from and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and storage medium can reside in an ASIC. The ASIC can reside in a user terminal. Alternatively, the processor and storage medium can exist as discrete components in the user terminal.

[0157] For software implementation, the techniques described in this application can be implemented using modules (e.g., procedures, functions, etc.) that perform the functions described in this application. This software code can be stored in memory units and executed by a processor. The memory units can be implemented within the processor or outside the processor; in the latter case, they are communicatively coupled to the processor via various means, as is well known in the art.

[0158] The foregoing description includes examples of one or more embodiments. It is certainly impossible to describe all possible combinations of components or methods in order to describe the above embodiments, but those skilled in the art will recognize that further combinations and arrangements of the various embodiments are possible. Therefore, the embodiments described herein are intended to cover all such changes, modifications, and variations that fall within the scope of the appended claims. Furthermore, the term "comprising" as used in the specification or claims is interpreted in a manner similar to the term "including," as interpreted when used as a conjunction in the claims. Additionally, the use of any term "or" in the specification of the claims is intended to mean "non-exclusive or."

Claims

1. A method for determining the Young's modulus of shale gas reservoirs based on elemental logging, characterized in that, include: Obtain mineral content, rock porosity, TOC content, water saturation, rock skeleton density, and rock density in shale gas reservoir rocks based on elemental logging data; Based on the mineral content, rock skeleton density, and predetermined mineral P-wave velocity, the rock skeleton P-wave velocity V is determined using the following formula (1). ps-el : （1）； Based on the mineral content, rock skeleton density, and predetermined mineral transverse wave velocity, the rock skeleton transverse wave velocity V is determined using the following formula (2). ss-el : （2）； In formulas (1) and (2), Density of the rock skeleton , , and These are the content, density, P-wave velocity, and S-wave velocity of the i-th mineral, respectively, where i is the index of the mineral type, i=1,2……n, n is the number of mineral types, and the mineral content is less than 1. Based on the P and S wave velocities of the rock skeleton, rock porosity, rock density, TOC content, and water saturation, and the pre-determined P and S wave velocities and TOC density of shale gas, TOC, and formation water, the P wave velocity of the shale gas reservoir rock is determined using the following formula (3). The shear wave velocity of shale gas reservoir rocks is determined using the following formula (4). : （3）； （4）； In formulas (3) and (4), and These are the longitudinal wave velocity and transverse wave velocity of the rock skeleton, respectively. and These are the longitudinal wave velocity and transverse wave velocity of formation water, respectively. and These are the P-wave velocity and S-wave velocity of shale gas, respectively. and These are the TOC longitudinal wave velocity and transverse wave velocity, respectively. and These are rock density and TOC density, respectively. , and These are TOC content, rock porosity, and water saturation, respectively. , and All less than 1; Based on the longitudinal and transverse wave velocities of the rock skeleton, the longitudinal and transverse wave velocities of the rock, the density of the rock skeleton, and the density of the rock, the Young's modulus of the rock skeleton and the Young's modulus of the rock are determined respectively.

2. The method as described in claim 1, characterized in that, The mineral content refers to the percentage of minerals by mass relative to the rock skeleton.

3. The method as described in claim 1, characterized in that, The determination of the rock skeleton modulus and the rock Young's modulus based on the rock skeleton P-wave velocity, rock P-wave velocity, rock skeleton density, and rock density, respectively, specifically includes: Determine the P- and S-wave transit times of the rock skeleton based on the P- and S-wave velocities of the rock skeleton. The Young's modulus of the rock skeleton is determined based on the rock skeleton density and the transit time of the rock skeleton P and S waves; the Young's modulus of the rock is determined based on the rock density and the transit time of the rock P and S waves.

4. The method as described in claim 3, characterized in that, The determination of the Young's modulus of the rock skeleton based on the rock skeleton density and the P-wave and S-wave transit times of the rock skeleton specifically includes: Based on the rock skeleton density and the longitudinal and transverse wave transit times of the rock skeleton, the Young's modulus of the rock skeleton is determined by the following formula (5): （5）； The Young's modulus of the rock is determined by the following formula (6) based on the rock density and the transit time of the longitudinal and transverse waves: （6）； In formulas (5) and (6), and These are the Young's modulus of the rock skeleton and the Young's modulus of the rock, respectively. and These are the rock skeleton density and the rock density, respectively. , , and These are the longitudinal wave time difference and transverse wave time difference of the rock skeleton, and the longitudinal wave time difference and transverse wave time difference of the rock itself.

5. The method according to any one of claims 1 to 4, characterized in that, The minerals include at least one of quartz, potassium feldspar, sodium feldspar, analcime, kaolinite, montmorillonite, illite, mica, calcite, dolomite, gypsum, and pyrite.

6. A device for determining the Young's modulus of shale gas reservoirs based on elemental logging, characterized in that, The apparatus is used to perform the method for determining the Young's modulus of shale gas reservoirs based on elemental logging as described in claim 1. The apparatus comprises: The P-wave and S-wave velocity determination module is used to determine the P-wave and S-wave velocities of the shale gas reservoir rock skeleton and the P-wave and S-wave velocities of the rock skeleton sequentially using the volume model method, based on the mineral content, rock porosity, TOC content, water saturation, rock skeleton density and rock density determined from the elemental logging data. The Young's modulus determination module is used to determine the Young's modulus of the rock skeleton and the Young's modulus of the rock based on the longitudinal and transverse wave velocities of the rock skeleton, the longitudinal and transverse wave velocities of the rock skeleton, the density of the rock skeleton, and the density of the rock.

7. A computer storage medium, characterized in that, The computer storage medium stores computer-executable instructions, which, when executed by a processor, implement the method for determining the Young's modulus of shale gas reservoirs based on element logging as described in any one of claims 1 to 5.

8. A server, characterized in that, include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the program, implements the method for determining the Young's modulus of shale gas reservoirs based on element logging as described in any one of claims 1 to 5.

Images