CN115963541A - Brittleness index calculation method and system for high-porosity reservoir and electronic equipment - Google Patents
Brittleness index calculation method and system for high-porosity reservoir and electronic equipment Download PDFInfo
- Publication number
- CN115963541A CN115963541A CN202111193071.4A CN202111193071A CN115963541A CN 115963541 A CN115963541 A CN 115963541A CN 202111193071 A CN202111193071 A CN 202111193071A CN 115963541 A CN115963541 A CN 115963541A
- Authority
- CN
- China
- Prior art keywords
- modulus
- brittleness index
- young
- rock
- wave velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention provides a method, a system and electronic equipment for calculating a brittleness index of a high-porosity reservoir stratum, wherein the method comprises the steps of calculating the Young modulus, the Poisson ratio and the shear modulus of logging data according to the longitudinal wave velocity, the transverse wave velocity and the density in the logging data; calculating to obtain a Lame coefficient according to the Young modulus and the shear modulus; and establishing a rock brittleness index equation based on the Young modulus, the Poisson ratio and the Lame coefficient, and obtaining the rock brittleness index through the rock brittleness index equation. The technical scheme of the invention has the beneficial effects that a new brittleness index is provided for a brittle gas-containing area with higher porosity, the new brittleness index is only related to the longitudinal and transverse wave velocity and the density, the influence of factors such as rock porosity, fluid, organic matters and the like on the Young modulus and Poisson ratio prediction is avoided, and a better indication effect is realized on the favorable brittle gas-containing area.
Description
Technical Field
The invention belongs to the technical field of geophysical exploration, and particularly relates to a method and a system for calculating a brittleness index of a high-porosity reservoir stratum and electronic equipment.
Background
The brittleness of a rock is a property in which the rock breaks suddenly (undergoes little plastic deformation before breaking) when subjected to a certain limit, and is released entirely in the form of elastic energy upon breaking. The brittleness coefficient (or brittleness index) is generally used to describe the strength of rock brittleness. A large number of researches find that the rock brittleness coefficient (B) can be determined according to a rock mechanical elasticity parameter method and a rock mineral composition method, and the calculation methods comprise the following two methods:
(1) Rock mechanical elastic parameter method: how to predict brittleness characteristics by using the statistical relationship among the rock physical and mechanical parameters is a problem which is tried to be solved by many petroleum companies and scientific research institutions at present. At present, the brittleness of the rock is calculated by adopting elastic modulus and Poisson ratio, and the elastic modulus and the Poisson ratio can better reflect the destructive capacity of the rock under the action of stress and during the formation of micro cracks. After the rock is fractured, the poisson's ratio can reflect the change of stress, and the elastic modulus reflects the ability to maintain fracture propagation.
(2) Rock mineral composition method: the researchers also suggest that the brittleness evaluation method based on the mineral composition is adopted, for example, the ratio of the brittle minerals to the total mineral content is considered, and the method has certain advantages in practical application. The brittleness can be roughly judged according to the content of the brittle mineral components by measuring the mineral content of the rock. The amount of brittleness is generally expressed in terms of the percentage of quartz in the total minerals (quartz + clay + carbonate rock, etc).
Scientific evaluation of brittleness requires the establishment of mechanical mechanisms for brittle fracture and failure of rocks. The existing research shows that the brittleness of the rock is not related to the peak strength of the rock, but is related to the slope and the residual strength of a stress-strain curve, and the brittleness characterization needs to consider two stages before and after the peak. Therefore, the brittleness evaluation method based on the full stress-strain mechanical characteristics can well express the macro-micro characteristics (pre-peak characteristics and post-peak characteristics) of the brittleness damage, is the main method of the brittleness test, and is the most direct method. However, although the brittleness evaluation method based on the full stress-strain is more reasonable in terms of rock fracture mechanism, the brittleness evaluation method is difficult to measure through an instrument in practical application, and the practical application of the brittleness evaluation method is difficult to realize.
Disclosure of Invention
The invention provides a method and a system for calculating a brittleness index of a high-porosity reservoir and electronic equipment, and provides a new brittleness index for a brittle gas-containing area with higher porosity, wherein the new brittleness index is only related to the longitudinal and transverse wave velocity and the density, so that the method has a better indication effect on the brittle gas-containing area, and has obvious advantages in distinguishing the brittle gas-containing area.
In order to achieve the above object, the present invention provides a method for calculating a brittleness index of a high porosity reservoir, comprising:
s1, calculating to obtain Young modulus, poisson ratio and shear modulus of logging data according to longitudinal wave velocity, transverse wave velocity and density in the logging data;
s2, calculating according to the Young modulus and the shear modulus to obtain a Lame coefficient;
and S3, establishing a rock brittleness index equation based on the Young modulus, the Poisson ratio and the Lame coefficient, and obtaining the rock brittleness index through the rock brittleness index equation.
Preferably, the rock brittleness index equation is expressed as: BI = E/λ σ;
wherein BI is the rock brittleness index, E is the Young's modulus, lambda is the Lame coefficient, and sigma is the Poisson's ratio.
Preferably, the young 'S modulus, the poisson' S ratio and the shear modulus in the step S1 are calculated by the following formula (1), formula (2) and formula (3), respectively:
μ＝V 2 s *ρ (3)
wherein E is the same asModulus, σ is the Poisson's ratio, μ is the shear modulus, V p Is the velocity of the longitudinal wave, V s ρ is the density, which is the shear wave velocity.
In a preferred example, the said ramen coefficient in step S2 is calculated from the following formula (4) based on formula (1) and formula (3):
wherein λ is the Lame coefficient, E is the Young's modulus, μ is the shear modulus, V p Is said longitudinal wave velocity, V s ρ is the density, which is the shear wave velocity.
Preferably, the rock brittleness index equation is obtained based on equation (1), equation (2) and equation (4):
wherein BI is the rock brittleness index, E is the Young modulus, lambda is the Lame coefficient, sigma is the Poisson's ratio, mu is the shear modulus, and V is p Is the velocity of the longitudinal wave, V s P is the density, which is the shear wave velocity.
The invention also provides a system for realizing the method for calculating the brittleness index of the high-porosity reservoir stratum, which comprises the following steps:
the first calculation module is used for calculating and obtaining Young modulus, poisson ratio and shear modulus of the logging data according to the longitudinal wave velocity, the transverse wave velocity and the density in the logging data;
the second calculation module is used for calculating a Lame coefficient according to the Young modulus and the shear modulus;
and the third calculation module is used for obtaining the rock brittleness index through a rock brittleness index equation according to the Young modulus, the Poisson ratio and the Lame coefficient.
Preferably, the rock brittleness index equation is expressed as: BI = E/λ σ;
wherein BI is the rock brittleness index, E is the Young's modulus, lambda is the Lame coefficient, and sigma is the Poisson's ratio.
Preferably, the first calculation module includes:
a first calculation unit that calculates the Young's modulus using the following equation:
a second calculating unit that calculates the poisson's ratio using the following equation:
a third calculating unit, wherein the third calculating unit calculates the shear modulus by using the following formula: μ = V 2 s *ρ；
Wherein E is the Young's modulus, σ is the Poisson's ratio, μ is the shear modulus, V p Is said longitudinal wave velocity, V s ρ is the density, which is the shear wave velocity.
Preferably, the second calculating module calculates the ramee coefficient by using the following formula:
and the third calculation module calculates the rock brittleness index by using the following formula:
wherein BI is the rock brittleness index, E is the Young's modulus, lambda is the Lame coefficient, sigma is the Poisson's ratio, muIs said shear modulus, V p Is the velocity of the longitudinal wave, V s ρ is the density, which is the shear wave velocity.
The present invention also provides an electronic device, including:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor, wherein,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method for calculating a brittleness index of a high porosity reservoir described above.
The technical scheme of the invention has the beneficial effects that:
the invention provides a new brittleness index E/lambda sigma for a brittle gas-containing area with higher porosity, the new brittleness index is only related to the speed and the density of longitudinal and transverse waves, the influence of factors such as rock porosity, fluid, organic matters and the like on the Young modulus E and Poisson ratio prediction is avoided, a good indication effect is also realized on the brittle gas-containing area, and meanwhile, the new brittleness index has obvious advantages in the judgment of the brittle gas-containing area.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in detail exemplary embodiments thereof with reference to the attached drawings, wherein like reference numerals generally represent like parts in the exemplary embodiments of the present invention.
FIG. 1 is a flow chart of a method of calculating a brittleness index of a high porosity reservoir according to the present invention;
FIG. 2 is a schematic structural diagram of a brittleness index calculation system of a high porosity reservoir according to the present invention;
fig. 3 is a diagram illustrating the effect of the method for calculating the brittleness index of the high-porosity reservoir according to the present invention.
Description of the reference numerals:
1. a first calculation module; 2. a second calculation module; 3. a second calculation module; 11. a first calculation unit; 12. a second calculation unit; 13. and a third calculation unit.
Detailed Description
In order to quantitatively describe the brittleness degree of the reservoir, the brittleness index directly representing the brittleness degree of the reservoir can be further calculated through rock mechanical parameters obtained by a prestack inversion method. The conventional method can better predict the area with stronger brittleness, but in the actual development process of shale gas, the Young modulus is found to increase along with the increase of the content of stones in the shale, but the Young modulus is reduced along with the increase of the porosity. Meanwhile, the Young modulus is also reduced along with the increase of the organic matter and pore gas content in the reservoir. Therefore, the conventional brittleness calculation mode has no effective characterization on the areas where high-porosity, high organic matter content and high-gas-content high-quality shale layers are developed, and in order to take both the areas into consideration, the invention provides a new brittleness index for the brittle gas-containing areas with higher porosity.
Preferred embodiments of the present invention will be described in more detail below. While the following describes preferred embodiments of the present invention, it should be understood that the present invention may be embodied in various forms and should not be limited by 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 invention to those skilled in the art.
It should be noted that the embodiments and features of the embodiments may be combined with each other without conflict.
The invention is further described with reference to the following drawings and specific examples, which are not intended to be limiting.
Referring to fig. 1, the invention provides a method for calculating a brittleness index of a high porosity reservoir, which includes:
s1, according to the longitudinal wave velocity V in the logging data p Transverse wave velocity V s And the density rho, and calculating to obtain the Young modulus E, the Poisson ratio sigma and the shear modulus mu of the logging data;
s2, calculating according to the Young modulus E and the shear modulus mu to obtain a Lame coefficient lambda;
and S3, establishing a rock brittleness index equation based on the Young modulus E, the Poisson ratio sigma and the Lame coefficient lambda, and obtaining the rock brittleness index through the rock brittleness index equation.
In a preferred example, the rock brittleness index equation is expressed as: BI = E/λ σ;
wherein BI is a rock brittleness index, E is a Young modulus, lambda is a Lame coefficient, and sigma is a Poisson's ratio.
As a preferred example, the young' S modulus E, the poisson ratio σ, and the shear modulus μ in step S1 are calculated by the following formula (1), formula (2), and formula (3), respectively:
μ＝V 2 s *ρ (3)
wherein E is Young's modulus, σ is Poisson's ratio, μ is shear modulus, and V p Is the velocity of the longitudinal wave, V s The shear wave velocity is denoted by ρ as the density.
As a preferred example, the ramen coefficient λ in step S2 is calculated from the following formula (4) based on the formulas (1) and (3):
wherein λ is Lame coefficient, E is Young's modulus, μ is shear modulus, V p Is the velocity of longitudinal wave, V s The shear wave velocity is denoted by ρ as the density.
A preferred example, the rock brittleness index equation is obtained based on formula (1), formula (2), and formula (4):
wherein BI is a rock brittleness index, E is a Young modulus, lambda is a Lame coefficient, sigma is a Poisson's ratio, mu is a shear modulus, and V is p Is the velocity of longitudinal wave, V s The shear wave velocity is denoted by ρ as the density.
Specifically, the logging data includes measured actual data: velocity V of longitudinal wave p Transverse wave velocity V s And the density rho, the Young modulus E can be calculated by the formula (1), then the Poisson ratio sigma is calculated by the formula (2), the Poisson ratio sigma is an elastic parameter and can macroscopically reflect underground lithology, and is an essential important parameter for predicting lithology and oil gas, the shear modulus mu is the ratio of shear stress to shear strain of the rock in the elastic deformation proportion limit range under the action of the shear stress, and represents the capability of the rock for resisting the shear strain, and the larger the modulus is, the rock rigidity is strong, and the rock rigidity is strong together with the longitudinal wave velocity V p Transverse wave velocity V s The relationship between the density and the density rho is expressed by a formula (3), the Lame coefficient is expressed by a formula (4), the Young modulus E and the shear modulus mu are calculated, a rock brittleness index equation is obtained based on a formula (1), a formula (2) and a formula (4), and then a new brittleness index BI = E/lambda sigma is obtained.
New brittleness index BI is only related to the longitudinal wave velocity V p Transverse wave velocity V s The method is related to the density rho, so that the influence of factors such as rock porosity, fluid and organic matters on the prediction of the Young modulus E and the Poisson ratio sigma is avoided, a favorable brittle gas-containing area is also well indicated, in addition, the new brittleness index BI has an obvious advantage in distinguishing the brittle gas-containing area, and the target interval can be more clearly identified through the new brittleness index BI.
Referring to fig. 2, the present invention further provides a system for implementing the method for calculating the brittleness index of the high porosity reservoir, including:
the first calculation module 1 is used for calculating and obtaining Young modulus, poisson ratio and shear modulus of the logging data according to the longitudinal wave velocity, the transverse wave velocity and the density in the logging data;
the second calculation module 2 is used for calculating a Lame coefficient according to the Young modulus and the shear modulus;
and the third calculating module 3 is used for obtaining the rock brittleness index through a rock brittleness index equation according to the Young modulus, the Poisson ratio and the Lame coefficient.
In a preferred example, the rock brittleness index equation is expressed as: BI = E/λ σ.
A preferred example, the first computing module 1 comprises:
the first calculation unit 11, the first calculation unit 11 calculates the young's modulus E using the following formula:
the second calculating unit 12, the second calculating unit 12 calculates the poisson's ratio σ by using the following formula:
the third calculating unit 13, the third calculating unit 13 calculates the shear modulus μ by using the following formula: μ = V 2 s *ρ；
Wherein E is Young's modulus, σ is Poisson's ratio, μ is shear modulus, and V p Is the velocity of the longitudinal wave, V s Is the shear wave velocity and ρ is the density.
In a preferred example, the second calculation module 2 calculates the lame coefficient λ using the following formula:
wherein λ is Lame coefficient, E is Young's modulus, μ is shear modulus, V p Is the velocity of the longitudinal wave, V s The shear wave velocity is denoted by ρ as the density.
In a preferred example, the third calculation module 3 calculates the rock brittleness index BI by using the following formula:
wherein BI is a rock brittleness index, E is a Young modulus, lambda is a Lame coefficient, sigma is a Poisson's ratio, mu is a shear modulus, and V is p Is the velocity of longitudinal wave, V s Is the shear wave velocity and ρ is the density.
In particular, in order to verify the feasibility of the brittleness index calculation method for the high-porosity reservoir stratum, the longitudinal wave velocity V of the logging data is utilized p Transverse wave velocity V s And density p were subjected to simulation testing. Through the first calculation module 1, the method is used for calculating the velocity V of longitudinal waves in the well logging data p Transverse wave velocity V s And the density rho, and calculating to obtain the Young modulus E, the Poisson ratio sigma and the shear modulus mu of the logging data; the second calculation module 2 is used for calculating a Lame coefficient lambda according to the Young modulus E and the shear modulus mu; the third calculation module 3 is configured to obtain a rock brittleness index BI = E/λ σ through a rock brittleness index equation according to the young's modulus E, the poisson ratio σ, and the ramen coefficient λ.
As shown in FIG. 3, the target layer is in the black box, and the new brittleness index BI is only related to the longitudinal wave velocity V p Transverse wave velocity V s And the method is related to the density rho, so that the influence of factors such as rock porosity, fluid and organic matters on the Young modulus E and Poisson ratio sigma prediction is avoided, a good indication effect is also achieved on a brittle gas-containing area, in addition, the new brittleness index BI has an obvious advantage in distinguishing the brittle gas-containing area, and a target layer section can be identified more clearly through the new brittleness index BI.
The present invention also provides an electronic device comprising:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor, wherein,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method for calculating a brittleness index of a high porosity reservoir described above.
Specifically, at least one processor executes executable instructions in the memory to perform the following steps:
s1, according to the longitudinal wave velocity V in the logging data p Transverse wave velocity V s And the density rho, and calculating to obtain the Young modulus E, the Poisson ratio sigma and the shear modulus mu of the logging data;
s2, calculating according to the Young modulus E and the shear modulus mu to obtain a Lame coefficient lambda;
and S3, establishing a rock brittleness index equation based on the Young modulus E, the Poisson ratio sigma and the Lame coefficient lambda, and obtaining a rock brittleness index BI through the rock brittleness index equation.
The new brittleness index BI obtained based on the above steps is only related to the longitudinal wave velocity V p Transverse wave velocity V s The method is related to the density rho, so that the influence of factors such as rock porosity, fluid and organic matters on the prediction of the Young modulus E and the Poisson ratio sigma is avoided, a good indication effect is achieved on a favorable brittle gas-containing area, and in addition, the new brittleness index BI has an obvious advantage in the discrimination of the brittle gas-containing area.
Example 1
Referring to fig. 1, the present embodiment provides a method for calculating a brittleness index of a high porosity reservoir, including:
s1, according to longitudinal wave velocity V in logging data p Transverse wave velocity V s And the density rho, and calculating to obtain the Young modulus E, the Poisson ratio sigma and the shear modulus mu of the logging data;
s2, calculating according to the Young modulus E and the shear modulus mu to obtain a Lame coefficient lambda;
and S3, establishing a rock brittleness index equation based on the Young modulus E, the Poisson ratio sigma and the Lame coefficient lambda, and obtaining the rock brittleness index through the rock brittleness index equation.
In this embodiment, the rock brittleness index equation is expressed as: BI = E/λ σ;
wherein BI is a rock brittleness index, E is a Young modulus, lambda is a Lame coefficient, and sigma is a Poisson ratio.
In this embodiment, the young' S modulus E, the poisson ratio σ, and the shear modulus μ in step S1 are calculated by the following formulas (1), (2), and (3), respectively:
μ＝V 2 s *ρ (3)
wherein E is Young's modulus, σ is Poisson's ratio, μ is shear modulus, V p Is the velocity of longitudinal wave, V s Is the shear wave velocity and ρ is the density.
In this embodiment, the ramen coefficient λ in step S2 is calculated from the following formula (4) based on the formulas (1) and (3):
wherein λ is Lame coefficient, E is Young's modulus, μ is shear modulus, and V is p Is the velocity of longitudinal wave, V s The shear wave velocity is denoted by ρ as the density.
In this embodiment, the rock brittleness index equation is obtained based on the formula (1), the formula (2), and the formula (4):
wherein BI is a rock brittleness index, E is a Young modulus, lambda is a Lame coefficient, sigma is a Poisson's ratio, mu is a shear modulus, and V is p Is the velocity of the longitudinal wave, V s The shear wave velocity is denoted by ρ as the density.
Example 2
Referring to fig. 2, the present embodiment provides a system for implementing the method for calculating a brittleness index of a high-porosity reservoir, including:
a first calculating module 1, wherein the first calculating module 1 is used for calculating the well logging number according to the well logging numberAccording to the longitudinal wave velocity V p Transverse wave velocity V s And the density rho, and calculating to obtain the Young modulus E, the Poisson ratio sigma and the shear modulus mu of the logging data;
the second calculation module 2 is used for calculating a Lame coefficient lambda according to the Young modulus E and the shear modulus mu;
the third calculating module 3 is used for obtaining the rock brittleness index BI through a rock brittleness index equation according to the Young modulus E, the Poisson ratio sigma and the Lame coefficient lambda.
In this embodiment, the rock brittleness index equation is expressed as: BI = E/λ σ.
In this embodiment, the first calculating module 1 includes:
the first calculation unit 11, the first calculation unit 11 calculates the young's modulus E using the following formula:
the second calculating unit 12, the second calculating unit 12 calculates the poisson's ratio σ by using the following formula:
the third calculating unit 13, the third calculating unit 13 calculates the shear modulus μ by using the following formula: μ = V 2 s *ρ；
Wherein E is Young's modulus, σ is Poisson's ratio, μ is shear modulus, and V p Is the velocity of longitudinal wave, V s Is the shear wave velocity and ρ is the density.
In this embodiment, the second calculating module 2 calculates the ramee coefficient λ by using the following formula:
wherein λ is Lame coefficient, E is Young's modulus, μ is shear modulus, and V is p Is the velocity of longitudinal wave, V s Is the transverse wave velocity, ρ isDensity.
In this embodiment, the third calculating module 3 calculates the rock brittleness index BI by using the following formula:
wherein BI is a rock brittleness index, E is a Young modulus, lambda is a Lame coefficient, sigma is a Poisson's ratio, mu is a shear modulus, and V is p Is the velocity of longitudinal wave, V s The shear wave velocity is denoted by ρ as the density.
In conclusion, the invention provides a new brittleness index BI = E/lambda sigma for a brittle gas-containing zone with higher porosity. New brittleness index BI is only related to the longitudinal wave velocity V p Transverse wave velocity V s The method is related to the density rho, so that the influence of factors such as rock porosity, fluid and organic matters on the prediction of the Young modulus E and the Poisson ratio sigma is avoided, a favorable brittle gas-containing area is also well indicated, in addition, the new brittleness index BI has an obvious advantage in distinguishing the brittle gas-containing area, and the target interval can be more clearly identified through the new brittleness index BI.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.
Claims (10)
1. A brittleness index calculation method of a high porosity reservoir stratum is characterized by comprising the following steps:
s1, calculating to obtain Young modulus, poisson ratio and shear modulus of logging data according to longitudinal wave velocity, transverse wave velocity and density in the logging data;
s2, calculating according to the Young modulus and the shear modulus to obtain a Lame coefficient;
and S3, establishing a rock brittleness index equation based on the Young modulus, the Poisson ratio and the Lame coefficient, and obtaining the rock brittleness index through the rock brittleness index equation.
2. The brittleness index calculation method of claim 1, wherein the rock brittleness index equation is expressed as: BI = E/λ σ;
wherein BI is the rock brittleness index, E is the Young's modulus, lambda is the Lame coefficient, and sigma is the Poisson's ratio.
3. The brittleness index calculation method according to claim 1, wherein the young 'S modulus, the poisson' S ratio, and the shear modulus in step S1 are calculated by the following formula (1), formula (2), and formula (3), respectively:
μ＝V 2 s *ρ (3)
wherein E is the Young's modulus, σ is the Poisson's ratio, μ is the shear modulus, V p Is the velocity of the longitudinal wave, V s P is the density, which is the shear wave velocity.
4. The brittleness index calculation method according to claim 3, wherein the Lame coefficient in step S2 is calculated from the following formula (4) based on formula (1) and formula (3):
wherein λ is the Lame coefficient, E is the Young's modulus, μ is the shear modulus, V p Is said longitudinal wave velocity, V s ρ is the density, which is the shear wave velocity.
5. The brittleness index calculation method according to claim 4, wherein the rock brittleness index equation is obtained based on formula (1), formula (2), and formula (4):
wherein BI is the rock brittleness index, E is the Young modulus, lambda is the Lame coefficient, sigma is the Poisson's ratio, mu is the shear modulus, and V is p Is the velocity of the longitudinal wave, V s ρ is the density, which is the shear wave velocity.
6. A brittleness index calculation system of a high porosity reservoir realizes the brittleness index calculation method of the high porosity reservoir according to any one of claims 1 to 5, and is characterized by comprising the following steps:
the first calculation module is used for calculating and obtaining Young modulus, poisson ratio and shear modulus of the logging data according to the longitudinal wave velocity, the transverse wave velocity and the density in the logging data;
the second calculation module is used for calculating a Lame coefficient according to the Young modulus and the shear modulus;
and the third calculation module is used for obtaining the rock brittleness index through a rock brittleness index equation according to the Young modulus, the Poisson ratio and the Lame coefficient.
7. The brittleness index calculation system of claim 6, wherein the rock brittleness index equation is expressed as: BI = E/λ σ;
wherein BI is the rock brittleness index, E is the Young's modulus, lambda is the Lame coefficient, and sigma is the Poisson's ratio.
8. The friability index calculation system of claim 6, wherein the first calculation module comprises:
a first calculation unit that calculates the Young's modulus using the following equation:
a second calculating unit that calculates the poisson's ratio using the following equation:
a third calculation unit that calculates the shear modulus using the following equation:
μ＝V 2 s *ρ；
wherein E is the Young's modulus, σ is the Poisson's ratio, μ is the shear modulus, V p Is said longitudinal wave velocity, V s ρ is the density, which is the shear wave velocity.
9. The brittleness index calculation system of claim 8, wherein the second calculation module calculates the lamel coefficient using the formula:
and the third calculation module calculates the rock brittleness index by using the following formula:
wherein BI is the rock brittlenessNumber, E is the Young's modulus, λ is the Lame coefficient, σ is the Poisson's ratio, μ is the shear modulus, V p Is said longitudinal wave velocity, V s P is the density, which is the shear wave velocity.
10. An electronic device, characterized in that the electronic device comprises:
at least one processor; and (c) a second step of,
a memory communicatively coupled to the at least one processor, wherein,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of calculating a brittleness index of a high porosity reservoir according to any one of claims 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111193071.4A CN115963541A (en) | 2021-10-13 | 2021-10-13 | Brittleness index calculation method and system for high-porosity reservoir and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111193071.4A CN115963541A (en) | 2021-10-13 | 2021-10-13 | Brittleness index calculation method and system for high-porosity reservoir and electronic equipment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115963541A true CN115963541A (en) | 2023-04-14 |
Family
ID=87360392
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111193071.4A Pending CN115963541A (en) | 2021-10-13 | 2021-10-13 | Brittleness index calculation method and system for high-porosity reservoir and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115963541A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116165054A (en) * | 2023-04-21 | 2023-05-26 | 煤炭科学研究总院有限公司 | Rock mechanical parameter acquisition method and device and electronic equipment |
-
2021
- 2021-10-13 CN CN202111193071.4A patent/CN115963541A/en active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116165054A (en) * | 2023-04-21 | 2023-05-26 | 煤炭科学研究总院有限公司 | Rock mechanical parameter acquisition method and device and electronic equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Investigation on the effect of freeze-thaw on fracture mode classification in marble subjected to multi-level cyclic loads | |
| Diederichs | The 2003 Canadian Geotechnical Colloquium: Mechanistic interpretation and practical application of damage and spalling prediction criteria for deep tunnelling | |
| Zhou et al. | Study of the effect of loading rates on crack propagation velocity and rock fracture toughness using cracked tunnel specimens | |
| CN105221141B (en) | A kind of mud shale brittleness index Forecasting Methodology | |
| Shukla et al. | Mechanical behaviour of reservoir rock under brine saturation | |
| Sarkar et al. | An empirical correlation of index geomechanical parameters with the compressional wave velocity | |
| Wen et al. | Dynamic compression characteristics of layered rock mass of significant strength changes in adjacent layers | |
| Liao et al. | Numerical simulation and experimental study of parallel seismic test for piles | |
| Zhou et al. | Study of the fracture behavior of mode I and mixed mode I/II cracks in tunnel under impact loads | |
| Damjanac et al. | Mechanical degradation of emplacement drifts at Yucca Mountain—A modeling case study: Part II: Lithophysal rock | |
| Gao et al. | Estimation of effect of voids on frequency response of mountain tunnel lining based on microtremor method | |
| Chen et al. | The ultrasonic P-wave velocity-stress relationship of rocks and its application | |
| Zhou et al. | The influence of impacting orientations on the failure modes of cracked tunnel | |
| Chen et al. | Experimental study on uniaxial compression failure modes and acoustic emission characteristics of fissured sandstone under water saturation | |
| Dong et al. | Study of mode I crack dynamic propagation behaviour and rock dynamic fracture toughness by using SCT specimens | |
| Zhou et al. | The effect of a pre‐existing crack on a running crack in brittle material under dynamic loads | |
| Du et al. | Dynamic compression properties and deterioration of red-sandstone subject to cyclic wet-dry treatment | |
| Vinck et al. | Advanced in situ and laboratory characterisation of the ALPACA chalk research site | |
| Wang et al. | Study on the effects of joints orientation and strength on failure behavior in shale specimen under impact loads | |
| Du et al. | A damage constitutive model of red sandstone under coupling of wet-dry cycles and impact load | |
| CN115963541A (en) | Brittleness index calculation method and system for high-porosity reservoir and electronic equipment | |
| Ali et al. | Investigation of the acoustic emission and fractal characteristics of coal with varying water contents during uniaxial compression failure | |
| Lee et al. | Effect of artificial cementation on cone tip resistance and small strain shear modulus of sand | |
| Niu et al. | Effect of water content on dynamic fracture characteristic of rock under impacts | |
| He et al. | Rock burst criterion based on clay mineral content |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |