Method and device for measuring gas saturation of shale reservoir

By generating the target density curve and gas-containing indicator curve, combined with the density logging and total organic carbon content curve of the shale reservoir, the problem of inaccurate calculation of gas-containing saturation in the shale reservoir in the prior art is solved, and a more accurate gas-containing saturation reflection and reservoir evaluation are achieved.

CN119920350APending Publication Date: 2025-05-02YANGTZE UNIVERSITY
View PDF 0 Cites 0 Cited by

Patent Information

Application Number
CN202411992506.5
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2024-12-31
Publication Date
2025-05-02

AI Technical Summary

Technical Problem

When the existing calculation of the gas saturation of shale reservoirs based on non-electrical logging information, the calculation results cannot accurately reflect the gas saturation.

Method used

By obtaining the skeleton density curve, density logging curve and total organic carbon content curve of the shale reservoir, the target density curve and gas-containing indicator curve are generated based on the difference of these curves, and finally the gas-containing indicator curve is normalized to obtain the gas-containing saturation curve.

Benefits of technology

This method can accurately reflect the gas saturation of the shale reservoir, peel off the influence of the skeleton and organic matter, and provide more accurate reservoir evaluation and reserve calculation.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119920350A_ABST
    Figure CN119920350A_ABST
Patent Text Reader

Abstract

The invention relates to a shale reservoir gas saturation measuring method and device, and belongs to the technical field of shale gas exploration, the method comprises the following steps: obtaining a skeleton density curve, a density logging curve and a total organic carbon content curve of a shale reservoir; a target density curve is obtained according to the difference value of the skeleton density curve and the density logging curve; according to a difference value between the target density curve and the total organic carbon content curve, obtaining a gas-containing indication curve; and carrying out normalization processing on the gas-bearing indication curve to obtain a gas saturation curve. According to the gas saturation curve obtained through calculation, the influence of a framework and organic matter is stripped, and the gas saturation can be accurately reflected.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to the technical field of shale gas exploration, and in particular to a method and device for measuring gas saturation of a shale reservoir. Background Art

[0002] Shale gas is natural gas that exists in organic-rich shale reservoirs in free and adsorbed states. China has abundant shale gas resources and huge development potential. Shale gas saturation is a crucial parameter for evaluating the amount of free gas in shale, and has important practical significance for shale gas reservoir evaluation and reserve calculation.

[0003] In recent years, Chinese and foreign scholars have conducted a lot of research on the calculation of water saturation in low-resistance gas reservoirs in shale reservoirs. One of the commonly used methods is to calculate the gas saturation of shale reservoirs based on non-electrical logging information. For example, the gas saturation is calculated by the ratio of the shale gas organic matter background value to the measured or calculated total organic carbon content of the shale gas layer, the gas saturation is calculated by the intersection diagram of density logging and neutron logging, and the gas saturation is calculated by equivalent medium theory. However, the gas saturation calculated by these methods does not only reflect the level of gas saturation, or it cannot accurately reflect the gas saturation. Summary of the invention

[0004] In view of this, it is necessary to provide a shale reservoir gas saturation measurement method and device to solve the problem that the calculation results of the existing shale reservoir gas saturation based on non-electrical logging information cannot accurately reflect the gas saturation.

[0005] In order to solve the above problems, in a first aspect, the present invention provides a method for measuring gas saturation of a shale reservoir, comprising: Obtain skeleton density curve, density logging curve and total organic carbon content curve of shale reservoir; Obtaining a target density curve according to a difference between the skeleton density curve and the density logging curve; Obtaining a gas content indicating curve according to a difference between the target density curve and the total organic carbon content curve; The gas content indicator curve is normalized to obtain a gas saturation curve.

[0006] Optionally, the shale reservoir includes a non-source rock region; and obtaining a gas content indicator curve according to a difference between the target density curve and the total organic carbon content curve includes: Adjusting the target density curve and the total organic carbon content curve so that the target density curve and the total organic carbon content curve overlap in a curve segment of the non-source rock region; The adjusted target density curve and the total organic carbon content curve are normalized and then subtracted to obtain a gas content indicator curve.

[0007] Optionally, the adjusting the target density curve and the total organic carbon content curve includes: The ordinates of the target density curve and the total organic carbon content curve are adjusted.

[0008] Optionally, the normalizing the adjusted target density curve and the total organic carbon content curve and then subtracting the two curves to obtain a gas content indicator curve includes: The adjusted target density curve and each point in the total organic carbon content curve are normalized and subtracted by the following formula to obtain the gas content indication value corresponding to each point;

[0009]

[0010]

[0011] in, Indicates the gas content indication value. Indicates the target density corresponding to the point in the adjusted target density curve, Indicates the minimum target density value in the adjusted target density curve. Indicates the maximum value of target density in the adjusted target density curve. Indicates the total organic carbon content corresponding to the point in the adjusted total organic carbon content curve, Indicates the minimum value of total organic carbon content in the adjusted total organic carbon content curve. Indicates the maximum value of total organic carbon content in the adjusted total organic carbon content curve; A gas content indication curve is obtained according to the gas content indication values ​​corresponding to the points.

[0012] Optionally, the normalizing the gas content indicator curve to obtain a gas saturation curve includes: Each point in the gas content indicator curve is normalized by the following formula to obtain the gas saturation corresponding to each point;

[0013] in, Gas saturation, Indicates the gas content indication value. represents the left scale value of the gas content indicator curve, Indicates the right scale value of the gas content indicator curve; A gas content indication curve is obtained according to the gas content indication values ​​corresponding to the points.

[0014] Optionally, obtaining a skeleton density curve of a shale reservoir includes: Obtaining well logging values ​​of target elements at multiple depth points in a shale reservoir; Determining target skeleton densities at the plurality of depth points according to the well logging values ​​of the target elements at the plurality of depth points; A skeleton density curve is generated according to the target skeleton densities at the plurality of depth points.

[0015] Optionally, obtaining the logging values ​​of the target element at multiple depth points in the shale reservoir includes: Obtain core samples at multiple depths in shale reservoirs; The core skeleton density is determined according to the following formula:

[0016] in, represents the core skeleton density, represents the content of each mineral component in the core sample, represents the density of each mineral in the core sample; The logging values ​​of the multiple elements at the multiple depth points are obtained, and the logging value of the target element is determined from the element logging values ​​of the multiple elements according to the correlation between the logging values ​​of the multiple elements and the core skeleton density.

[0017] Optionally, the target elements include aluminum, iron, potassium and sulfur; and determining the target skeleton density at the multiple depth points according to the logging values ​​of the target elements at the multiple depth points includes: The target skeleton density at the multiple depth points is determined by the following formula:

[0018] in, represents the target skeleton density, Indicates the logging value of aluminum element, Indicates the logging value of iron element, represents the logging value of potassium, Indicates the logging value of sulfur element, is the correlation coefficient.

[0019] Optionally, obtaining a total organic carbon content curve of a shale reservoir includes: Obtaining total carbon content curve and inorganic carbon content curve of shale reservoir; A total organic carbon content curve is generated according to the total carbon content curve and the inorganic carbon content curve.

[0020] In a second aspect, the present invention further provides a shale reservoir gas saturation measurement device, comprising: An information acquisition module, used to obtain a skeleton density curve, a density logging curve and a total organic carbon content curve of a shale reservoir; A target density curve generating module, used for obtaining a target density curve according to the difference between the skeleton density curve and the density logging curve; A gas content indication curve generating module, used for obtaining a gas content indication curve according to the difference between the target density curve and the total organic carbon content curve; The gas saturation curve generating module is used to perform normalization processing on the gas content indicator curve to obtain a gas saturation curve.

[0021] The beneficial effects of the present invention are: After obtaining the skeleton density curve, density logging curve and total organic carbon content curve of the shale reservoir, the present invention first generates a target density curve according to the difference between the skeleton density curve and the density logging curve; then generates a gas content indicator curve according to the difference between the target density curve and the total organic carbon content curve; finally, the gas content indicator curve is normalized to obtain a gas saturation curve. The shale reservoir is composed of pores, organic matter and a skeleton, wherein the pores are filled with gas and water. The density logging curve reflects the comprehensive density of the entire shale reservoir. Based on the density logging curve of the shale reservoir, the present invention gradually removes the influence of the skeleton and organic matter of the shale reservoir from the density logging curve to obtain a curve that only retains the pore information, namely, the gas content indicator curve. Finally, the gas content indicator curve is normalized to obtain a gas saturation curve. The gas saturation curve calculated by the present invention removes the influence of the skeleton and organic matter and can accurately reflect the gas saturation. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] Figure 1 A schematic diagram of a shale reservoir volume model provided by the present invention; Figure 2 A schematic flow chart of an embodiment of a method for measuring gas saturation of a shale reservoir provided by the present invention; Figure 3 A schematic diagram of the relationship between skeleton density and element logging response provided by the present invention; Figure 4 A total carbon content curve graph and an inorganic carbon content curve graph provided by the present invention; Figure 5 A two-step stripping shale gas saturation measurement flow chart provided by the present invention; Figure 6A schematic structural diagram of an embodiment of a shale reservoir gas saturation measuring device provided by the present invention. DETAILED DESCRIPTION

[0023] The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0024] In the description of the embodiments of the present invention, unless otherwise specified, the meaning of "plurality" is two or more. The "first", "second", etc. involved in the embodiments of the present invention are used to distinguish similar objects, and are not used to describe a specific order or sequence, nor are they used to indicate or imply their relative importance or implicitly indicate the number of technical features indicated. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application can be implemented in an order other than those illustrated or described here, and the objects distinguished by "first", "second", etc. are generally of one type, and the number of objects is not limited. For example, the first object can be one or more.

[0025] Reference to "embodiments" herein means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present invention. The appearance of the phrase in various places in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0026] Reference Figure 1 , shows a schematic diagram of a shale reservoir volume model provided by the present invention. Shale reservoirs are composed of pores (filled with gas and water), organic matter and skeletons, and the density logging value reflects the comprehensive density response of the entire shale reservoir.

[0027] Reference Figure 2 , shows a schematic flow chart of an embodiment of a method for measuring gas saturation of a shale reservoir provided by the present invention, the method comprising: S201, obtaining a skeleton density curve, a density logging curve and a total organic carbon content curve of a shale reservoir.

[0028] The horizontal axis of the skeleton density curve is depth, and the vertical axis is skeleton density, which is used to indicate the change of skeleton density with depth in shale reservoirs. The horizontal axis of the density logging curve is depth, and the vertical axis is density logging value, which is used to indicate the change of density logging value with depth in shale reservoirs. The density logging curve can be obtained by density logging instrument testing. The horizontal axis of the total organic carbon content curve is depth, and the vertical axis is total organic carbon content, which is used to indicate the change of total organic carbon content with depth in shale reservoirs.

[0029] S202, obtaining a target density curve according to the difference between the skeleton density curve and the density logging curve.

[0030] Align the abscissas of the skeleton density curve and the density logging curve, then subtract the ordinates, and use the difference between the ordinates as the new ordinate to obtain the target density curve (the abscissa of the target density curve is the depth, and the ordinate is the difference calculated above).

[0031] Specifically, the vertical coordinates can be subtracted by the following formula:

[0032] Where: is the target density, g / cm 3 ; is the density logging value, g / cm 3 .

[0033] The target density curve is the density curve without the influence of the skeleton.

[0034] S203, obtaining a gas content indicator curve according to the difference between the target density curve and the total organic carbon content curve.

[0035] Align the target density curve with the abscissa of the total organic carbon content curve, and then use the difference between the ordinates of the target density curve and the total organic carbon content curve as the new ordinate to obtain the gas content indicator curve. The gas content indicator curve is a curve that further strips away the influence of organic matter on the basis of the target density curve and only retains the pore information.

[0036] S204, normalizing the gas content indicator curve to obtain a gas saturation curve.

[0037] The gas saturation curve can be obtained by normalizing the ordinate of the gas content indicator curve. The abscissa of the gas saturation curve is depth, and the ordinate is gas saturation.

[0038] The present invention is based on the density logging curve of the shale reservoir, and gradually removes the influence of the skeleton and organic matter of the shale reservoir from the density logging curve to obtain a curve that only retains pore information, that is, a gas content indicator curve, and finally normalizes the gas content indicator curve to obtain a gas saturation curve. The gas saturation curve calculated by the present invention removes the influence of the skeleton and organic matter and can accurately reflect the gas saturation.

[0039] In one embodiment, the shale reservoir includes a non-source rock area. In the shale reservoir, the non-source rock refers to the rock part that exists in the shale reservoir but does not have the ability to generate hydrocarbons. Shale reservoirs are usually dominated by organic-rich source rocks, but may also contain some non-source rock interlayers or lithological units. Then S203 may include: adjusting the target density curve and the total organic carbon content curve so that the target density curve and the total organic carbon content curve overlap in the non-source rock area; normalizing the adjusted target density curve and the total organic carbon content curve and then subtracting them to obtain a gas content indicator curve.

[0040] The target density curve and the total organic carbon content curve are aligned on the horizontal axis, and then the vertical scales of the two curves are adjusted. For example, the two curves can be moved along the vertical axis so that the curve segments of the two curves in the non-source rock area overlap. Then, the adjusted target density curve and the total organic carbon content curve are normalized and subtracted to obtain a gas content indicator curve.

[0041] In one embodiment, the step of normalizing the adjusted target density curve and the total organic carbon content curve and then subtracting them to obtain the gas content indication curve may include: normalizing the adjusted target density curve and each point in the total organic carbon content curve and then subtracting them according to the following formula to obtain the gas content indication value corresponding to each point;

[0042]

[0043]

[0044] In the formula, Indicates the gas content indication value. Indicates the target density corresponding to the point in the adjusted target density curve, Indicates the minimum target density value in the adjusted target density curve. Indicates the maximum value of target density in the adjusted target density curve. Indicates the total organic carbon content corresponding to the point in the adjusted total organic carbon content curve, Indicates the minimum value of total organic carbon content in the adjusted total organic carbon content curve. Indicates the maximum value of total organic carbon content in the adjusted total organic carbon content curve; Then, a gas content indication curve is obtained according to the gas content indication value corresponding to each point.

[0045] In one embodiment, step S204 may include: normalizing each point in the gas content indicator curve by the following formula to obtain the gas saturation corresponding to each point;

[0046] in, Gas saturation, Indicates the gas content indication value. Indicates the left scale value of the gas content indicator curve. Indicates the right scale value of the gas content indicator curve; A gas content indication curve is obtained according to the gas content indication value corresponding to each point.

[0047] Each gas content indicator curve has a scale range at the top of the graph, the smallest scale display value is defined as the left scale value, and the largest scale display is defined as the right scale value. It can also be understood as the minimum value of the gas content indication value. It can also be understood as the maximum value of the gas content indication value.

[0048] In one embodiment, the step of obtaining a skeleton density curve of a shale reservoir may include: obtaining logging values ​​of a target element at multiple depth points in the shale reservoir through a formation element capture logging instrument; determining target skeleton densities at multiple depth points based on the logging values ​​of the target element at multiple depth points; and generating a skeleton density curve based on the target skeleton densities at multiple depth points.

[0049] In one embodiment, the step of obtaining the logging values ​​of the target element at multiple depth points in the shale reservoir may include: obtaining core samples at multiple depth points in the shale reservoir through a coring well; then performing an X-ray diffraction experiment on the core samples to obtain the mineral composition and mineral content of each core sample, and determining the core skeleton density of the core sample according to the following formula:

[0050] In the formula, represents the core skeleton density, Indicates the content of each mineral component in the core sample, represents the density of each mineral in the core sample, i Represents quartz, plagioclase, calcite, dolomite, pyrite, and clay minerals (illite).

[0051] Afterwards, the formation element capture logging instrument is used to measure the logging values ​​of multiple elements at the same depth as the core sample, including aluminum, calcium, iron, potassium, magnesium, manganese, sodium, silicon, sulfur, etc. And according to the correlation between the logging values ​​of the multiple elements and the core skeleton density, the logging value of the target element is determined from the element logging values ​​of the multiple elements.

[0052] In one embodiment, the target elements include aluminum, iron, potassium and sulfur; the step of determining the target skeleton density at multiple depth points according to the logging values ​​of the target elements at multiple depth points may include: determining the target skeleton density at multiple depth points by the following formula:

[0053] in, represents the target skeleton density, Indicates the logging value of aluminum element, Indicates the logging value of iron element, represents the logging value of potassium, Indicates the logging value of sulfur element, is the correlation coefficient.

[0054] Reference Figure 3 , showing a schematic diagram of the relationship between skeleton density and element logging response provided by the present invention. Figure 3 From the last column of values, we can see that aluminum, iron, potassium and sulfur have a good correlation with the core skeleton density.

[0055] In one embodiment, the step of obtaining a total organic carbon content curve of a shale reservoir may include: obtaining a total carbon content curve and an inorganic carbon content curve of the shale reservoir; and subtracting the total carbon content curve from the inorganic carbon content curve to generate a total organic carbon content curve.

[0056] The formation element capture logging tool can simultaneously measure the inelastic scattering gamma spectrum and the captured gamma spectrum. The inelastic scattering gamma spectrum can be analyzed to obtain the total organic carbon content (TOC) of the formation. The total organic carbon content can be obtained by deducting the inorganic carbon content (TIC) from the total carbon content (TC). The calculation process is as follows:

[0057] The depth is preset at each interval, such as measuring TOC once every 0.125 m, and a total organic carbon content curve can be drawn based on each measurement value.

[0058] In addition to being directly measured and calculated by formation element capture logging, total organic carbon content can also be calculated based on conventional logging data using the ΔlogR method or curve overlap method. For wells that have been measured using formation element capture logging or conventional logging curves, the total organic carbon content can be obtained using the above method.

[0059] Reference Figure 4 , showing a total carbon content curve and an inorganic carbon content curve provided by the present invention. The TC curve is subtracted from the TIC curve to obtain the TOC curve.

[0060] Reference Figure 5 , showing a two-step stripping shale gas saturation measurement flow chart provided by the present invention.

[0061] (1) First, 50 core samples were drilled from the coring well (Well A) of the shale reservoir, and X-ray diffraction experiments were performed on the core samples to obtain the mineral composition and mineral content of each core sample, and the core skeleton density of each core sample was calculated based on the mineral composition and mineral content. Secondly, the logging values ​​of the formation elements at the same depth as the core samples were obtained. The formation elements included aluminum (Al), calcium (Ca), iron (Fe), potassium (K), magnesium (Mg), manganese (Mn), sodium (Na), silicon (Si), sulfur (S), etc. The correlation between the core skeleton density at each depth and the logging values ​​of the multiple elements was analyzed, and the logging values ​​of sensitive elements, such as aluminum, iron, potassium, and sulfur, were optimized. A statistical model of the target skeleton density and the logging values ​​of sensitive elements was established. The statistical model is as follows:

[0062] Based on the statistical model and the logging values ​​of aluminum, iron, potassium and sulfur measured at various depths, the target density at each depth can be calculated; according to the target density at each depth, the skeleton density curve of the shale reservoir can be obtained ( curve).

[0063] (2) Generate the TOC curve based on the TC curve and TIC curve.

[0064] (3) Taking Well A as an example, the gas-bearing indicator curve is obtained by stripping in two steps. First, Curve and The curve (density logging curve) is subtracted to remove the influence of the skeleton and obtain Curve. Secondly, TOC curve and The vertical scale of the curve is adjusted so that , -0.1, 0.2 respectively; , are -0.5 and 0.9 respectively, so that the TOC curve is The curves overlap in the non-source rock section. Then the adjusted TOC curve is compared with After the curves are normalized, they are subtracted to remove the influence of organic matter and obtain the gas indicator curve (HI curve). Finally, the scale of the ordinate of the HI curve is adjusted so that , -0.5 and 0.2 respectively, to make it consistent with the gas saturation of the core. Then the HI curve is normalized to obtain the gas saturation curve ( curve).

[0065] Reference Figure 6 , shows a schematic structural diagram of an embodiment of a shale reservoir gas saturation measurement device provided by the present invention, the device 60 comprises: Information acquisition module 601, used to obtain the skeleton density curve, density logging curve and total organic carbon content curve of the shale reservoir; The target density curve generating module 602 is used to obtain the target density curve according to the difference between the skeleton density curve and the density logging curve; A gas content indication curve generating module 603 is used to obtain a gas content indication curve according to the difference between the target density curve and the total organic carbon content curve; The gas saturation curve generating module 604 is used to normalize the gas content indicator curve to obtain a gas saturation curve.

[0066] It should be noted that the implementation principles or implementation processes of the above modules can refer to the above-mentioned embodiments of the shale reservoir gas saturation measurement method, and will not be described one by one here.

[0067] Those skilled in the art will appreciate that all or part of the processes of the above-mentioned embodiments can be implemented by instructing related hardware through a computer program, and the program can be stored in a computer-readable storage medium, wherein the computer-readable storage medium is a disk, an optical disk, a read-only storage memory, or a random access memory, etc.

[0068] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily conceived by any technician familiar with the technical field within the technical scope disclosed by the present invention should be covered within the protection scope of the present invention.

Claims

1. A method for measuring gas saturation in shale reservoirs, characterized in that: include: Obtain skeleton density curve, density logging curve and total organic carbon content curve of shale reservoir; Obtaining a target density curve according to a difference between the skeleton density curve and the density logging curve; Obtaining a gas content indicating curve according to a difference between the target density curve and the total organic carbon content curve; The gas content indicator curve is normalized to obtain a gas saturation curve.

2. The method for measuring gas saturation of shale reservoir according to claim 1, characterized in that: The shale reservoir includes a non-source rock region; the gas content indication curve is obtained according to the difference between the target density curve and the total organic carbon content curve, including: Adjusting the target density curve and the total organic carbon content curve so that the target density curve and the total organic carbon content curve overlap in a curve segment of the non-source rock region; The adjusted target density curve and the total organic carbon content curve are normalized and then subtracted to obtain a gas content indicator curve.

3. The method for measuring gas saturation of shale reservoir according to claim 2, characterized in that: The adjusting the target density curve and the total organic carbon content curve includes: The ordinates of the target density curve and the total organic carbon content curve are adjusted.

4. The method for measuring gas saturation of shale reservoir according to claim 2, characterized in that: The adjusted target density curve is normalized and then subtracted from the total organic carbon content curve to obtain a gas content indicator curve, including: The adjusted target density curve and each point in the total organic carbon content curve are normalized and subtracted by the following formula to obtain the gas content indication value corresponding to each point; in, Indicates the gas content indication value. Indicates the target density corresponding to the point in the adjusted target density curve, Indicates the minimum target density value in the adjusted target density curve. Indicates the maximum value of target density in the adjusted target density curve. Indicates the total organic carbon content corresponding to the point in the adjusted total organic carbon content curve, Indicates the minimum value of total organic carbon content in the adjusted total organic carbon content curve. Indicates the maximum value of total organic carbon content in the adjusted total organic carbon content curve; A gas content indication curve is obtained according to the gas content indication values ​​corresponding to the points.

5. The method for measuring gas saturation of shale reservoir according to claim 1, characterized in that: The normalizing the gas content indicator curve to obtain a gas saturation curve includes: Each point in the gas content indicator curve is normalized by the following formula to obtain the gas saturation corresponding to each point; in, Gas saturation, Indicates the gas content indication value. represents the left scale value of the gas content indicator curve, Indicates the right scale value of the gas content indicator curve; A gas content indication curve is obtained according to the gas content indication values ​​corresponding to the points.

6. The method for measuring gas saturation of shale reservoir according to claim 1, characterized in that: The step of obtaining a skeleton density curve of a shale reservoir comprises: Obtaining well logging values ​​of target elements at multiple depth points in a shale reservoir; Determining target skeleton densities at the plurality of depth points according to the well logging values ​​of the target elements at the plurality of depth points; A skeleton density curve is generated according to the target skeleton densities at the plurality of depth points.

7. The method for measuring gas saturation of shale reservoirs according to claim 6, characterized in that: The method of obtaining the well logging values ​​of the target element at multiple depth points in the shale reservoir comprises: Obtain core samples at multiple depths in shale reservoirs; The core skeleton density is determined according to the following formula: in, represents the core skeleton density, represents the content of each mineral component in the core sample, represents the density of each mineral in the core sample; The logging values ​​of the multiple elements at the multiple depth points are obtained, and the logging value of the target element is determined from the element logging values ​​of the multiple elements according to the correlation between the logging values ​​of the multiple elements and the core skeleton density.

8. The method for measuring gas saturation of shale reservoirs according to claim 6, characterized in that: The target elements include aluminum, iron, potassium and sulfur; and determining the target skeleton density at the multiple depth points according to the well logging values ​​of the target elements at the multiple depth points includes: The target skeleton density at the multiple depth points is determined by the following formula: in, represents the target skeleton density, Indicates the logging value of aluminum element, Indicates the logging value of iron element, represents the logging value of potassium, Indicates the logging value of sulfur element, is the correlation coefficient.

9. The method for measuring gas saturation of shale reservoirs according to claim 1, characterized in that: The method of obtaining a total organic carbon content curve of a shale reservoir comprises: Obtaining total carbon content curve and inorganic carbon content curve of shale reservoir; A total organic carbon content curve is generated according to the total carbon content curve and the inorganic carbon content curve.

10. A shale reservoir gas saturation measurement device, characterized in that: include: An information acquisition module, used to obtain a skeleton density curve, a density logging curve and a total organic carbon content curve of a shale reservoir; A target density curve generating module, used for obtaining a target density curve according to the difference between the skeleton density curve and the density logging curve; A gas content indication curve generating module, used for obtaining a gas content indication curve according to the difference between the target density curve and the total organic carbon content curve; The gas saturation curve generating module is used to perform normalization processing on the gas content indicator curve to obtain a gas saturation curve.