Method and system for calculating brittleness index of high porosity reservoirs and electronic device

Abstract

The application provides a method and system for calculating a brittleness index of a high porosity reservoir and an electronic device, comprising calculating Young's modulus, Poisson's ratio and shear modulus of logging data according to P-wave velocity, S-wave velocity and density in the logging data; calculating Lame's coefficient according to the Young's modulus and the shear modulus; establishing a rock brittleness index equation based on the Young's modulus, the Poisson's ratio and the Lame's coefficient, and obtaining a rock brittleness index through the rock brittleness index equation. The technical scheme of the application has the beneficial effect of providing a new brittleness index for a brittle gas-bearing area with high porosity, which is only related to P-wave velocity, S-wave velocity and density, thereby avoiding the influence of rock porosity, fluid, organic matter and other factors on the prediction of Young's modulus and Poisson's ratio, and having a good indicating effect on a favorable brittle gas-bearing area.

Description

Technical Field

[0001] This invention belongs to the field of geophysical exploration technology, and more specifically, relates to a method, system and electronic equipment for calculating the brittleness index of high-porosity reservoirs. Background Technology

[0002] The brittleness of rocks refers to their property of suddenly fracturing when subjected to a certain stress limit (with minimal plastic deformation occurring before fracturing), releasing all energy as elastic energy upon fracture. The brittleness coefficient (or brittleness index) is generally used to describe the degree of rock brittleness. Extensive research has found that the brittleness coefficient (B) can be determined using the rock mechanics elastic parameter method and the rock mineral composition method. The calculation methods are as follows:

[0003] (1) Rock mechanics elastic parameter method: How to predict brittle characteristics using the statistical relationship between rock physical and mechanical parameters is a problem that many oil companies and research institutions are trying to solve. Currently, elastic modulus and Poisson's ratio are mostly used to calculate rock brittleness, as it is believed that elastic modulus and Poisson's ratio can better reflect the failure capacity of rocks under stress and microcrack formation. After cracks are formed in the rock, Poisson's ratio can reflect the change in stress, while elastic modulus reflects the ability to maintain crack propagation.

[0004] (2) Rock and mineral composition method: Some scholars have also suggested using a brittleness evaluation method based on mineral composition, such as considering the proportion of brittle minerals in the total mineral content. This type of method has certain advantages in practical applications. By measuring the rock and mineral content, the strength of brittleness can be roughly judged based on the content of brittle mineral components. Usually, the percentage of quartz in the total minerals (quartz + clay + carbonate rocks, etc.) is used to represent the degree of brittleness.

[0005] Scientific evaluation of brittleness requires understanding the mechanical mechanisms of brittle fracture and failure in rocks. Existing research indicates that the strength of rock brittleness is independent of peak strength but related to the slope of the stress-strain curve and residual strength. Characterizing brittleness necessitates considering both the pre-peak and post-peak stages. Therefore, brittleness evaluation methods based on full stress-strain mechanical characteristics can better represent the macroscopic and microscopic features (pre-peak and post-peak characteristics) of brittle failure, making them the primary and most direct method for brittleness testing. However, while full stress-strain-based brittleness evaluation methods are more reasonable from the perspective of rock fracture mechanisms, they are difficult to measure with instruments in practical applications, posing certain challenges to their practical use. Summary of the Invention

[0006] This invention provides a method, system, and electronic device for calculating the brittleness index of high-porosity reservoirs. A new brittleness index is proposed for brittle gas-bearing areas with high porosity. The new brittleness index is only related to the P-wave and S-wave velocities and density, and has a good indicative effect on favorable brittle gas-bearing areas. It has a significant advantage in identifying brittle gas-bearing areas.

[0007] To achieve the above objectives, the present invention provides a method for calculating the brittleness index of high-porosity reservoirs, comprising:

[0008] Step S1: Based on the P-wave velocity, S-wave velocity, and density in the well logging data, calculate the Young's modulus, Poisson's ratio, and shear modulus of the well logging data.

[0009] Step S2: Calculate the Lamé coefficient based on the Young's modulus and the shear modulus;

[0010] Step S3: Based on the Young's modulus, the Poisson's ratio, and the Lamé coefficient, establish a rock brittleness index equation, and obtain the rock brittleness index through the rock brittleness index equation.

[0011] Preferably, the rock brittleness index equation is expressed as: BI=E / λσ;

[0012] Wherein, BI is the rock brittleness index, E is the Young's modulus, λ is the Lamé coefficient, and σ is the Poisson's ratio.

[0013] Preferably, the Young's modulus, the Poisson's ratio, and the shear modulus in step S1 are calculated by the following formulas (1), (2), and (3), respectively:

[0014]

[0015]

[0016] μ = V 2 s *ρ (3)

[0017] Where E is the Young's modulus, σ is the Poisson's ratio, μ is the shear modulus, and V p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0018] In a preferred example, the Lamé coefficient in step S2 is calculated based on formulas (1) and (3) by the following formula (4):

[0019]

[0020] Where λ is the Lamé coefficient, E is the Young's modulus, μ is the shear modulus, and V is the shear modulus. p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0021] Preferably, the rock brittleness index equation is obtained based on formulas (1), (2), and (4):

[0022]

[0023] Wherein, BI is the rock brittleness index, E is the Young's modulus, λ is the Lamé coefficient, σ is the Poisson's ratio, μ is the shear modulus, and V... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0024] The present invention also provides a system for implementing the above-described method for calculating the brittleness index of high-porosity reservoirs, comprising:

[0025] The first calculation module is used to calculate the Young's modulus, Poisson's ratio and shear modulus of the logging data based on the P-wave velocity, S-wave velocity and density in the logging data.

[0026] The second calculation module is used to calculate the Lamé coefficient based on the Young's modulus and the shear modulus.

[0027] The third calculation module is used to obtain the rock brittleness index based on the Young's modulus, the Poisson's ratio, and the Lamé coefficient through the rock brittleness index equation.

[0028] Preferably, the rock brittleness index equation is expressed as: BI=E / λσ;

[0029] Wherein, BI is the rock brittleness index, E is the Young's modulus, λ is the Lamé coefficient, and σ is the Poisson's ratio.

[0030] Preferably, the first calculation module includes:

[0031] The first calculation unit calculates the Young's modulus using the following formula:

[0032]

[0033] The second calculation unit calculates the Poisson's ratio using the following formula:

[0034] The third calculation unit calculates the shear modulus using the following formula: μ = V2 s *ρ;

[0035] Where E is the Young's modulus, σ is the Poisson's ratio, μ is the shear modulus, and V p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0036] Preferably, the second calculation module calculates the Lamé coefficient using the following formula:

[0037]

[0038] Furthermore, the third calculation module calculates the rock brittleness index using the following formula:

[0039]

[0040] Wherein, BI is the rock brittleness index, E is the Young's modulus, λ is the Lamé coefficient, σ is the Poisson's ratio, μ is the shear modulus, and V... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0041] The present invention also provides an electronic device, the electronic device comprising:

[0042] At least one processor; and,

[0043] The memory communicatively connected to the at least one processor, wherein,

[0044] The memory stores instructions that can be executed by the at least one processor, which enables the at least one processor to perform the brittleness index calculation method for high-porosity reservoirs described above.

[0045] The beneficial effects of the technical solution of the present invention are as follows:

[0046] This invention proposes a new brittleness index E / λσ for brittle gas-bearing regions with high porosity. The new brittleness index is only related to P-wave and S-wave velocities and density. It not only avoids the influence of rock porosity, fluid, organic matter and other factors on the prediction of Young's modulus E and Poisson's ratio, but also has a good indicative effect on favorable brittle gas-bearing regions. At the same time, the new brittleness index has a significant advantage in identifying brittle gas-bearing regions. Attached Figure Description

[0047] The above and other objects, features and advantages of the present invention will become more apparent from the more detailed description of exemplary embodiments of the invention in conjunction with the accompanying drawings, wherein the same reference numerals generally represent the same components in the exemplary embodiments of the invention.

[0048] Figure 1 This is a flowchart of a method for calculating the brittleness index of a high-porosity reservoir according to the present invention;

[0049] Figure 2 This is a schematic diagram of the structure of a brittleness index calculation system for high-porosity reservoirs according to the present invention;

[0050] Figure 3 This is a diagram illustrating the implementation effect of the brittleness index calculation method for high-porosity reservoirs according to the present invention.

[0051] Explanation of reference numerals in the attached figures:

[0052] 1. First calculation module; 2. Second calculation module; 3. Second calculation module; 11. First calculation unit; 12. Second calculation unit; 13. Third calculation unit. Detailed Implementation

[0053] To quantitatively describe the brittleness of reservoirs, rock mechanics parameters obtained through pre-stack inversion methods can be further used to calculate a brittleness index that directly characterizes the brittleness of the reservoir. Conventional methods can predict areas with high brittleness relatively well. However, in actual shale gas development, it has been found that while Young's modulus increases with increasing stony content in shale, it decreases with increasing porosity. Simultaneously, increasing organic matter and gas content in the reservoir also leads to a decrease in Young's modulus. Therefore, conventional brittleness calculation methods lack effective characterization for areas with high porosity, high organic matter content, and high gas content in high-quality shale formations. To balance these two aspects, this invention proposes a new brittleness index for brittle gas-bearing areas with high porosity.

[0054] Preferred embodiments of the invention will now be described in more detail. While preferred embodiments of the invention are described below, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that the invention will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.

[0055] It should be noted that, unless otherwise specified, the embodiments and features described in the present invention can be combined with each other.

[0056] The present invention will be further described below with reference to the accompanying drawings and specific embodiments, but this is not intended to limit the scope of the invention.

[0057] Reference Figure 1 As shown, the present invention provides a method for calculating the brittleness index of high-porosity reservoirs, comprising:

[0058] Step S1, based on the P-wave velocity V in the well logging data p Shear wave velocity V s Based on the density ρ, the Young's modulus E, Poisson's ratio σ, and shear modulus μ of the well logging data are calculated.

[0059] Step S2: Calculate the Lamé coefficient λ based on Young's modulus E and shear modulus μ;

[0060] Step S3: Based on Young's modulus E, Poisson's ratio σ and Lamé coefficient λ, establish the rock brittleness index equation, and obtain the rock brittleness index through the rock brittleness index equation.

[0061] In a preferred example, the rock brittleness index equation is expressed as: BI=E / λσ;

[0062] Where BI is the rock brittleness index, E is Young's modulus, λ is Lamé coefficient, and σ is Poisson's ratio.

[0063] In a preferred example, the Young's modulus E, Poisson's ratio σ, and shear modulus μ in step S1 are calculated by the following formulas (1), (2), and (3), respectively:

[0064]

[0065]

[0066] μ = V 2 s *ρ (3)

[0067] Where E is Young's modulus, σ is Poisson's ratio, μ is shear modulus, and V is... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0068] In a preferred example, the Lamé coefficient λ in step S2 is calculated from the following formula (4) based on formulas (1) and (3):

[0069]

[0070] Where λ is the Lamé coefficient, E is Young's modulus, μ is the shear modulus, and V is the shear modulus. p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0071] A preferred example is to obtain the rock brittleness index equation based on formulas (1), (2), and (4):

[0072]

[0073] Where BI is the rock brittleness index, E is Young's modulus, λ is Lamé coefficient, σ is Poisson's ratio, μ is shear modulus, and V is... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0074] Specifically, the logging data includes the actual measured data: P-wave velocity V p transverse wave velocity V s Given the density ρ, the Young's modulus E can be calculated using formula (1), and then the Poisson's ratio σ can be obtained using formula (2). The Poisson's ratio σ is an elastic parameter that can macroscopically reflect underground lithology and is an essential parameter for predicting lithology and oil and gas. The shear modulus μ is the ratio of shear stress to shear strain of a rock under shear stress within the elastic deformation proportional limit. It characterizes the rock's ability to resist shear strain. The larger the modulus, the stronger the rock's rigidity. It is related to the longitudinal wave velocity V. p Shear wave velocity V s The relationship between density ρ and the coefficient of Lamé is expressed by formula (3), and the coefficient of Lamé is expressed by formula (4). It is calculated from Young's modulus E and shear modulus μ. Based on formula (1), formula (2) and formula (4), the rock brittleness index equation is obtained, and then the new brittleness index BI=E / λσ is obtained.

[0075] The new brittleness index BI is only related to the longitudinal wave velocity V. p Shear wave velocity V s It is related to density ρ, which avoids the influence of factors such as rock porosity, fluid, and organic matter on the prediction of Young's modulus E and Poisson's ratio σ. At the same time, it also has a good indicative effect on favorable brittle gas-bearing areas. In addition, the new brittleness index BI has a significant advantage in identifying brittle gas-bearing areas. Through the new brittleness index BI, the target layer can be identified more clearly.

[0076] Reference Figure 2 As shown, the present invention also provides a system for implementing the above-described method for calculating the brittleness index of high-porosity reservoirs, comprising:

[0077] The first calculation module 1 is used to calculate the Young's modulus, Poisson's ratio and shear modulus of the logging data based on the P-wave velocity, S-wave velocity and density in the logging data.

[0078] The second calculation module 2 is used to calculate the Lamé coefficient based on Young's modulus and shear modulus.

[0079] The third calculation module 3 is used to obtain the rock brittleness index based on Young's modulus, Poisson's ratio, and Lamé coefficient through the rock brittleness index equation.

[0080] In a preferred example, the rock brittleness index equation is expressed as: BI = E / λσ.

[0081] In a preferred example, the first computing module 1 includes:

[0082] The first calculation unit 11 calculates Young's modulus E using the following formula:

[0083]

[0084] The second calculation unit 12 calculates the Poisson's ratio σ using the following formula:

[0085] The third calculation unit 13 calculates the shear modulus μ using the following formula: μ=V 2 s *ρ;

[0086] Where E is Young's modulus, σ is Poisson's ratio, μ is shear modulus, and V is... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0087] In a preferred example, the second calculation module 2 calculates the Lamé coefficient λ using the following formula:

[0088]

[0089] Where λ is the Lamé coefficient, E is Young's modulus, μ is the shear modulus, and V is the shear modulus. p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0090] In a preferred example, the third calculation module 3 calculates the rock brittleness index BI using the following formula:

[0091]

[0092] Where BI is the rock brittleness index, E is Young's modulus, λ is Lamé coefficient, σ is Poisson's ratio, μ is shear modulus, and V is... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0093] Specifically, to verify the feasibility of this method for calculating the brittleness index of high-porosity reservoirs, the P-wave velocity V from well logging data was used. p Shear wave velocity V sSimulation tests were conducted using density ρ. The first calculation module 1 was used to calculate the P-wave velocity V from the well logging data. p Shear wave velocity V s The first calculation module 1 calculates Young's modulus E, Poisson's ratio σ, and shear modulus μ from the well logging data based on the density ρ. The second calculation module 2 calculates the Lamé coefficient λ based on Young's modulus E and shear modulus μ. The third calculation module 3 obtains the rock brittleness index BI = E / λσ based on Young's modulus E, Poisson's ratio σ, and Lamé coefficient λ through the rock brittleness index equation.

[0094] like Figure 3 As shown, the target layer is within the black box. The new brittleness index BI is only related to the P-wave velocity V. p Shear wave velocity V s It is related to density ρ, which avoids the influence of factors such as rock porosity, fluid, and organic matter on the prediction of Young's modulus E and Poisson's ratio σ. At the same time, it also has a good indicative effect on favorable brittle gas-bearing areas. In addition, the new brittleness index BI has a significant advantage in identifying brittle gas-bearing areas. Through the new brittleness index BI, the target layer can be identified more clearly.

[0095] The present invention also provides an electronic device, comprising:

[0096] At least one processor; and,

[0097] A memory that is communicatively connected to at least one processor, wherein,

[0098] The memory stores instructions that can be executed by at least one processor, which enables the at least one processor to perform the above-described method for calculating the brittleness index of high-porosity reservoirs.

[0099] Specifically, at least one processor executes executable instructions in memory to perform the following steps:

[0100] Step S1, based on the P-wave velocity V in the well logging data p Shear wave velocity V s Based on the density ρ, the Young's modulus E, Poisson's ratio σ, and shear modulus μ of the well logging data are calculated.

[0101] Step S2: Calculate the Lamé coefficient λ based on Young's modulus E and shear modulus μ;

[0102] Step S3: Based on Young's modulus E, Poisson's ratio σ and Lamé coefficient λ, establish the rock brittleness index equation, and obtain the rock brittleness index BI through the rock brittleness index equation.

[0103] The new brittleness index BI obtained based on the above steps is only related to the longitudinal wave velocity V. p Shear wave velocity V sIt is related to density ρ, which avoids the influence of factors such as rock porosity, fluid, and organic matter on the prediction of Young's modulus E and Poisson's ratio σ. At the same time, it also has a good indicative effect on favorable brittle gas-bearing areas. In addition, the new brittleness index BI has a significant advantage in identifying brittle gas-bearing areas.

[0104] Example 1

[0105] Reference Figure 1 As shown in the figure, this embodiment provides a method for calculating the brittleness index of high-porosity reservoirs, including:

[0106] Step S1, based on the P-wave velocity V in the well logging data p Shear wave velocity V s Based on the density ρ, the Young's modulus E, Poisson's ratio σ, and shear modulus μ of the well logging data are calculated.

[0107] Step S2: Calculate the Lamé coefficient λ based on Young's modulus E and shear modulus μ;

[0108] Step S3: Based on Young's modulus E, Poisson's ratio σ and Lamé coefficient λ, establish the rock brittleness index equation, and obtain the rock brittleness index through the rock brittleness index equation.

[0109] In this embodiment, the rock brittleness index equation is expressed as: BI=E / λσ;

[0110] Where BI is the rock brittleness index, E is Young's modulus, λ is Lamé coefficient, and σ is Poisson's ratio.

[0111] In this embodiment, the Young's modulus E, Poisson's ratio σ, and shear modulus μ in step S1 are calculated by the following formulas (1), (2), and (3):

[0112]

[0113]

[0114] μ = V 2 s *ρ (3)

[0115] Where E is Young's modulus, σ is Poisson's ratio, μ is shear modulus, and V is... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0116] In this embodiment, the Lamé coefficient λ in step S2 is calculated from the following formula (4) based on formulas (1) and (3):

[0117]

[0118] Where λ is the Lamé coefficient, E is Young's modulus, μ is the shear modulus, and V is the shear modulus. p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0119] In this embodiment, the rock brittleness index equation is obtained based on formulas (1), (2), and (4):

[0120]

[0121] Where BI is the rock brittleness index, E is Young's modulus, λ is Lamé coefficient, σ is Poisson's ratio, μ is shear modulus, and V is... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0122] Example 2

[0123] Reference Figure 2 As shown, this embodiment provides a system for implementing the above-described method for calculating the brittleness index of high-porosity reservoirs, comprising:

[0124] First calculation module 1, the first calculation module 1 is used to calculate the P-wave velocity V in the well logging data. p Shear wave velocity V s Based on the density ρ, the Young's modulus E, Poisson's ratio σ, and shear modulus μ of the well logging data are calculated;

[0125] The second calculation module 2 is used to calculate the Lamé coefficient λ based on Young's modulus E and shear modulus μ.

[0126] The third calculation module 3 is used to obtain the rock brittleness index BI based on Young's modulus E, Poisson's ratio σ and Lamé coefficient λ through the rock brittleness index equation.

[0127] In this embodiment, the rock brittleness index equation is expressed as: BI=E / λσ.

[0128] In this embodiment, the first calculation module 1 includes:

[0129] The first calculation unit 11 calculates Young's modulus E using the following formula:

[0130]

[0131] The second calculation unit 12 calculates the Poisson's ratio σ using the following formula:

[0132] The third calculation unit 13 calculates the shear modulus μ using the following formula: μ=V 2s *ρ;

[0133] Where E is Young's modulus, σ is Poisson's ratio, μ is shear modulus, and V is... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0134] In this embodiment, the second calculation module 2 calculates the Lamé coefficient λ using the following formula:

[0135]

[0136] Where λ is the Lamé coefficient, E is Young's modulus, μ is the shear modulus, and V is the shear modulus. p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0137] In this embodiment, the third calculation module 3 calculates the rock brittleness index BI using the following formula:

[0138]

[0139] Where BI is the rock brittleness index, E is Young's modulus, λ is Lamé coefficient, σ is Poisson's ratio, μ is shear modulus, and V is... p V is the longitudinal wave velocity. s Let ρ be the transverse wave velocity and ρ be the density.

[0140] In summary, this invention proposes a new brittleness index BI = E / λσ for brittle gas-bearing regions with high porosity. The new brittleness index BI is only related to the longitudinal wave velocity V. p Shear wave velocity V s It is related to density ρ, which avoids the influence of factors such as rock porosity, fluid, and organic matter on the prediction of Young's modulus E and Poisson's ratio σ. At the same time, it also has a good indicative effect on favorable brittle gas-bearing areas. In addition, the new brittleness index BI has a significant advantage in identifying brittle gas-bearing areas. Through the new brittleness index BI, the target layer can be identified more clearly.

[0141] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for calculating the brittleness index of high-porosity reservoirs, characterized in that, include: Step S1: Based on the P-wave velocity, S-wave velocity, and density in the well logging data, calculate the Young's modulus, Poisson's ratio, and shear modulus of the well logging data. Step S2: Calculate the Lamé coefficient based on the Young's modulus and the shear modulus; Step S3: Based on the Young's modulus, the Poisson's ratio, and the Lamé coefficient, establish a rock brittleness index equation, and obtain the rock brittleness index through the rock brittleness index equation; The rock brittleness index equation is expressed as: ; Wherein, BI is the brittleness index of the rock. The Young's modulus, Let Lamé coefficient be the coefficient. Let be the Poisson's ratio.

2. The method for calculating the brittleness index according to claim 1, characterized in that, The Young's modulus, Poisson's ratio, and shear modulus in step S1 are calculated by the following formulas (1), (2), and (3): （1） （2） （3） in, The Young's modulus, Let Poisson's ratio be the given value. The shear modulus, The longitudinal wave velocity, The transverse wave velocity, The density is stated.

3. The method for calculating the brittleness index according to claim 2, characterized in that, The Lamé coefficient in step S2 is calculated from the following formula (4) based on formulas (1) and (3): （4） in, Let Lamé coefficient be the coefficient. The Young's modulus, The shear modulus, The longitudinal wave velocity, The transverse wave velocity, The density is stated.

4. The method for calculating the brittleness index according to claim 3, characterized in that, The rock brittleness index equation is obtained based on formulas (1), (2), and (4): ； Wherein, BI is the brittleness index of the rock. The Young's modulus, Let Lamé coefficient be the coefficient. Let Poisson's ratio be the given value. The shear modulus, The longitudinal wave velocity, The transverse wave velocity, The density is stated.

5. A brittleness index calculation system for high-porosity reservoirs, implementing the brittleness index calculation method for high-porosity reservoirs as described in any one of claims 1-4, characterized in that, include: The first calculation module is used to calculate the Young's modulus, Poisson's ratio and shear modulus of the logging data based on the P-wave velocity, S-wave velocity and density in the logging data. The second calculation module is used to calculate the Lamé coefficient based on the Young's modulus and the shear modulus. The third calculation module is used to obtain the rock brittleness index based on the Young's modulus, the Poisson's ratio, and the Lamé coefficient through the rock brittleness index equation. The rock brittleness index equation is expressed as: ； Wherein, BI is the brittleness index of the rock. The Young's modulus, Let Lamé coefficient be the coefficient. Let be the Poisson's ratio.

6. The fragility index calculation system according to claim 5, characterized in that, The first computing module includes: The first calculation unit calculates the Young's modulus using the following formula: ； The second calculation unit calculates the Poisson's ratio using the following formula: ; The third calculation unit calculates the shear modulus using the following formula: ; in, The Young's modulus, Let Poisson's ratio be the given value. The shear modulus, The longitudinal wave velocity, The transverse wave velocity, The density is stated.

7. The brittleness index calculation system according to claim 6, characterized in that, The second calculation module calculates the Lamé coefficient using the following formula: ； Furthermore, the third calculation module calculates the rock brittleness index using the following formula: ； Wherein, BI is the brittleness index of the rock. The Young's modulus, Let Lamé coefficient be the coefficient. Let Poisson's ratio be the given value. The shear modulus, The longitudinal wave velocity, The transverse wave velocity, The density is stated.

8. An electronic device, characterized in that, The electronic device includes: At least one processor; and, The memory communicatively connected to the at least one processor, wherein, The memory stores instructions that can be executed by the at least one processor, which, when executed by the at least one processor, enables the at least one processor to perform the brittleness index calculation method for high-porosity reservoirs according to any one of claims 1-4.

Images