CN110764144A - Shale reservoir shear wave velocity calculation method and system - Google Patents

Shale reservoir shear wave velocity calculation method and system Download PDF

Info

Publication number
CN110764144A
CN110764144A CN201810828120.9A CN201810828120A CN110764144A CN 110764144 A CN110764144 A CN 110764144A CN 201810828120 A CN201810828120 A CN 201810828120A CN 110764144 A CN110764144 A CN 110764144A
Authority
CN
China
Prior art keywords
shale reservoir
skeleton model
rock skeleton
stiffness tensor
calculating
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
Application number
CN201810828120.9A
Other languages
Chinese (zh)
Inventor
郑四连
刘百红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
China Petrochemical Corp
Original Assignee
Sinopec Geophysical Research Institute
China Petrochemical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sinopec Geophysical Research Institute, China Petrochemical Corp filed Critical Sinopec Geophysical Research Institute
Priority to CN201810828120.9A priority Critical patent/CN110764144A/en
Publication of CN110764144A publication Critical patent/CN110764144A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

A shale reservoir shear wave velocity calculation method and system are disclosed. The method can comprise the following steps: establishing a rock skeleton model, and calculating the elastic modulus and the density of the rock skeleton model; calculating the elastic stiffness tensor of the single-layer isotropy of the rock skeleton model according to the elastic modulus and the density of the rock skeleton model; calculating the transverse isotropic elastic stiffness tensor of the rock skeleton model according to the single-layer isotropic elastic stiffness tensor of the rock skeleton model; adding fractures into the rock skeleton model to obtain a shale reservoir model, and calculating the elastic stiffness tensor of the shale reservoir model; and calculating the shear wave velocity of the shale reservoir model according to the elastic stiffness tensor and the logging density of the shale reservoir model. According to the method, the shale reservoir transverse wave velocity is quickly, simply and conveniently calculated through element logging and rock physics, and necessary parameters are provided for shale gas reservoir prediction, brittleness evaluation and fracture prediction.

Description

Shale reservoir shear wave velocity calculation method and system
Technical Field
The invention relates to the field of oil and gas seismic exploration and development, in particular to a shale reservoir shear wave velocity calculation method and a shale reservoir shear wave velocity calculation system.
Background
The seismic prestack inversion and prestack attribute analysis are important means for reservoir prediction and fluid prediction at present, and the longitudinal wave velocity, the transverse wave velocity and the density are necessary parameters for the seismic prestack inversion and prestack attribute analysis. On the other hand, along with the continuous rising of energy demand of countries in the world, unconventional energy sources such as shale gas, dense gas, coal bed gas and the like are used as the supplement of conventional energy sources, and the attention of people is gradually drawn. The shale gas reservoir belongs to a typical unconventional natural gas reservoir with low permeability and low porosity, and occupies a great proportion in oil and gas resources in China. However, the development cost is high, the development difficulty is high, and the special reservoir characteristics determine that a strengthening means, namely a reservoir fracturing modification technology, is required to be adopted for developing the reservoir, the oil-gas flow seepage condition is improved, so that the aim of effective exploitation is fulfilled, and the transverse wave velocity is a key parameter for researching the mechanical properties of the reservoir. In addition, compared to conventional sandstone reservoirs, shales have a more complex mineral composition and pore structure and tend to develop natural microfractures. The micro-fractures in the shale are both reservoir spaces for oil and gas and channels for oil and gas migration. The high and low output of the shale gas is directly related to the development degree of natural microcracks in the shale, and the existence of the microcracks improves the effectiveness of a hydraulic fracturing effect to a certain extent, so that the seepage capability of the shale is greatly improved, and a necessary migration channel is provided for the shale gas to enter a well hole from a bedrock pore. And the shear wave velocity is an important parameter for evaluating the development degree of reservoir fractures.
Typically, both velocity and density of longitudinal waves can be obtained from well log data. Because the shear wave logging technology is developed late and the cost of shear wave logging is high, shear wave velocity logging information does not exist in a plurality of work areas. The conventional shear wave velocity calculation method is based on empirical formulas, and the empirical formulas are obtained based on sandstone reservoirs. Therefore, it is necessary to develop a method and a system for calculating the shear wave velocity of the shale reservoir.
The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides a shale reservoir shear wave velocity calculation method and a shale reservoir shear wave velocity calculation system, which can quickly, simply and conveniently calculate the shale reservoir shear wave velocity through element logging and petrophysics and provide necessary parameters for shale gas reservoir prediction, brittleness evaluation and fracture prediction.
According to one aspect of the invention, a shale reservoir shear wave velocity calculation method is provided. The method may include: establishing a rock skeleton model, and calculating the elastic modulus and the density of the rock skeleton model; calculating the elastic stiffness tensor of the single-layer isotropy of the rock skeleton model according to the elastic modulus and the density of the rock skeleton model; calculating the transverse isotropic elastic stiffness tensor of the rock skeleton model according to the single-layer isotropic elastic stiffness tensor of the rock skeleton model; adding fractures into the rock skeleton model to obtain a shale reservoir model, and calculating the elastic stiffness tensor of the shale reservoir model; and calculating the shear wave velocity of the shale reservoir model according to the elastic stiffness tensor and the logging density of the shale reservoir model.
Preferably, the building of the rock skeleton model, and the calculating of the elastic modulus of the rock skeleton model comprises: establishing the rock skeleton model according to mineral components and volume content of rock obtained by element logging; and calculating the Voigt average modulus and the Reuss average modulus of the rock skeleton model, and further calculating the elastic modulus of the rock skeleton model.
Preferably, the Voigt average modulus is:
Figure BDA0001742951260000021
wherein M isVIs the Voigt mean modulus, fiAnd MiThe volume content and the elastic modulus of the ith mineral component, i ═ 1, 2, …, N, respectively.
Preferably, the reus mean modulus is:
Figure BDA0001742951260000031
wherein M isRIs the reus mean modulus.
Preferably, the elastic modulus of the rock skeleton model is:
wherein M is the elastic modulus of the rock skeleton model.
Preferably, the single layer isotropic elastic stiffness tensor is:
Figure BDA0001742951260000033
wherein, ICijThe elastic stiffness tensor is a single-layer isotropic elastic stiffness tensor, wherein a, b, c, d, f and m are parameters of the single-layer isotropic elastic stiffness tensor, a is c and lambda +2 mu, b is f and lambda, and d is m and mu.
Preferably, the transversely isotropic elastic stiffness tensor is:
Figure BDA0001742951260000034
wherein, VCijIs the transversely isotropic elastic stiffness tensor, A, B, C, D, F, M is a parameter of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000035
M＝<m>wherein the operator<x>Represents a pair of brackets<>The mineral composition x of the rock in (a) is weighted averaged by volume content.
Preferably, the elastic stiffness tensor of the shale reservoir model is:
Figure BDA0001742951260000041
wherein, TCijElastic stiffness tensor, Delta, for shale reservoir modelNIn order to provide rigidity for crack initiation,
Figure BDA0001742951260000042
ΔHin order to provide the horizontal stiffness of the crack,
Figure BDA0001742951260000043
g is the parametric ratio of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000044
and e represents the crack density.
Preferably, the shear wave velocity of the shale reservoir model is:
Figure BDA0001742951260000045
wherein v issThe shear wave velocity is denoted by ρ as the density.
According to another aspect of the invention, a shale reservoir shear wave velocity calculation system is provided, which is characterized by comprising: a memory storing computer-executable instructions; a processor executing computer executable instructions in the memory to perform the steps of: establishing a rock skeleton model, and calculating the elastic modulus and the density of the rock skeleton model; calculating the elastic stiffness tensor of the single-layer isotropy of the rock skeleton model according to the elastic modulus and the density of the rock skeleton model; calculating the transverse isotropic elastic stiffness tensor of the rock skeleton model according to the single-layer isotropic elastic stiffness tensor of the rock skeleton model; adding fractures into the rock skeleton model to obtain a shale reservoir model, and calculating the elastic stiffness tensor of the shale reservoir model; and calculating the shear wave velocity of the shale reservoir model according to the elastic stiffness tensor and the logging density of the shale reservoir model.
Preferably, the building of the rock skeleton model, and the calculating of the elastic modulus of the rock skeleton model comprises: establishing the rock skeleton model according to mineral components and volume content of rock obtained by element logging; and calculating the Voigt average modulus and the Reuss average modulus of the rock skeleton model, and further calculating the elastic modulus of the rock skeleton model.
Preferably, the Voigt average modulus is:
Figure BDA0001742951260000051
wherein M isVIs the Voigt mean modulus, fiAnd MiThe volume content and the elastic modulus of the ith mineral component, i ═ 1, 2, …, N, respectively.
Preferably, the reus mean modulus is:
Figure BDA0001742951260000052
wherein M isRIs the reus mean modulus.
Preferably, the elastic modulus of the rock skeleton model is:
Figure BDA0001742951260000053
wherein M is the elastic modulus of the rock skeleton model.
Preferably, the single layer isotropic elastic stiffness tensor is:
Figure BDA0001742951260000054
wherein, ICijThe elastic stiffness tensor is a single-layer isotropic elastic stiffness tensor, wherein a, b, c, d, f and m are parameters of the single-layer isotropic elastic stiffness tensor, a is c and lambda +2 mu, b is f and lambda, and d is m and mu.
Preferably, the transversely isotropic elastic stiffness tensor is:
Figure BDA0001742951260000061
wherein, VCijIs the transversely isotropic elastic stiffness tensor, A, B, C, D, F, M is a parameter of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000062
Figure BDA0001742951260000063
M＝<m>wherein the operator<x>Represents a pair of brackets<>The mineral composition x of the rock in (a) is weighted averaged by volume content.
Preferably, the elastic stiffness tensor of the shale reservoir model is:
Figure BDA0001742951260000064
wherein, TCijElastic stiffness tensor, Delta, for shale reservoir modelNIn order to provide rigidity for crack initiation,
Figure BDA0001742951260000065
ΔHin order to provide the horizontal stiffness of the crack,
Figure BDA0001742951260000066
g is the parametric ratio of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000067
and e represents the crack density.
Preferably, the shear wave velocity of the shale reservoir model is:
Figure BDA0001742951260000068
wherein v issThe shear wave velocity is denoted by ρ as the density.
The present invention has other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts.
Fig. 1 shows a flow chart of the steps of a shale reservoir shear wave velocity calculation method according to the present invention.
Detailed Description
The invention will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Fig. 1 shows a flow chart of the steps of a shale reservoir shear wave velocity calculation method according to the present invention.
In this embodiment, the shale reservoir shear wave velocity calculation method according to the present invention may include: step 101, establishing a rock skeleton model, and calculating the elastic modulus and the density of the rock skeleton model; 102, calculating the elastic stiffness tensor of the single-layer isotropy of the rock skeleton model according to the elastic modulus and the density of the rock skeleton model; 103, calculating the transverse isotropic elastic stiffness tensor of the rock framework model according to the single-layer isotropic elastic stiffness tensor of the rock framework model; 104, adding fractures into the rock skeleton model to obtain a shale reservoir model, and calculating the elastic stiffness tensor of the shale reservoir model; and 105, calculating the shear wave velocity of the shale reservoir model according to the elastic stiffness tensor and the logging density of the shale reservoir model.
In one example, building a rock skeleton model, calculating an elastic modulus of the rock skeleton model comprises: establishing a rock skeleton model according to mineral components and volume content of rocks obtained by element logging; and calculating the Voigt average modulus and the Reuss average modulus of the rock skeleton model, and further calculating the elastic modulus of the rock skeleton model.
In one example, the Voigt average modulus is:
Figure BDA0001742951260000081
wherein M isVIs the Voigt mean modulus, fiAnd MiThe volume content and the elastic modulus of the ith mineral component, i ═ 1, 2, …, N, respectively.
In one example, the reus mean modulus is:
Figure BDA0001742951260000082
wherein M isRIs the reus mean modulus.
In one example, the modulus of elasticity of the rock skeleton model is:
Figure BDA0001742951260000083
wherein M is the elastic modulus of the rock skeleton model.
In one example, the elastic stiffness tensor for the single layer isotropy is:
Figure BDA0001742951260000084
wherein, ICijThe elastic stiffness tensor is a single-layer isotropic elastic stiffness tensor, wherein a, b, c, d, f and m are parameters of the single-layer isotropic elastic stiffness tensor, a is c and lambda +2 mu, b is f and lambda, and d is m and mu.
In one example, the transversely isotropic elastic stiffness tensor is:
Figure BDA0001742951260000091
wherein, VCijIs the transversely isotropic elastic stiffness tensor, A, B, C, D, F, M is a parameter of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000092
Figure BDA0001742951260000093
M＝<m>wherein the operator<x>Represents a pair of brackets<>The mineral composition x of the rock in (a) is weighted averaged by volume content.
In one example, the elastic stiffness tensor of the shale reservoir model is:
Figure BDA0001742951260000094
wherein, TCijElastic stiffness tensor, Delta, for shale reservoir modelNIn order to provide rigidity for crack initiation,
Figure BDA0001742951260000095
ΔHin order to provide the horizontal stiffness of the crack,
Figure BDA0001742951260000096
g is the parametric ratio of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000097
and e represents the crack density.
In one example, the shear wave velocity of the shale reservoir model is:
wherein v issThe shear wave velocity is denoted by ρ as the density.
In particular, formation element logging is a neutron-gamma spectroscopy logging technology which obtains element content by measuring gamma rays emitted by the action of neutrons and formation element nuclei so as to determine mineral content. Fast neutrons emitted by the high-energy neutron source are slowed down and absorbed by nuclear reactions such as fast neutron inelastic scattering, neutron activation, and thermal neutron capture. The gamma rays emitted by different elements are characteristic during the nuclear reaction with neutrons. Therefore, the element content in the rock can be obtained through the analysis of secondary gamma energy spectrum such as inelastic scattering spectrum, thermal neutron capture spectrum and the like. According to the mineral combination model established by different lithological profiles, the rock mineral composition and content can be calculated according to the element content. And conventional logging information such as density, neutron, photoelectric index, natural gamma and the like is combined to obtain lithologic components, porosity, water saturation and kerogen content of the shale reservoir.
Common minerals that constitute shale reservoirs are feldspar, quartz, calcite, dolomite, kaolinite, illite, montmorillonite, pyrite, mica, anhydrite, and the like. There is also kerogen, and fluid in the reservoir, and structurally shale also has a pronounced thin lamellar structure and near vertical fissures. The elastic parameters of these minerals, such as modulus, density and velocity, can be measured indoors and are therefore considered known. The volume content of the corresponding mineral can be obtained from the formation element log.
Firstly, the minerals and kerogen are taken as matrix, the organic porosity and the intergranular porosity are taken as matrix porosity, a rock skeleton model is established, and M of the rock skeleton model is calculated according to the volume content and the elastic modulus of the mineralsV、MRRespectively in formula (1) and formula (2), and further calculating the elastic modulus and density of the rock skeleton model by using Voigt-reus-Hill average calculation, wherein the density can be obtained by conventional well logging, or can be obtained by using formula (1) and formula (2), namely after the volume content of each mineral in the rock is obtained by element well logging, the elastic modulus M of each mineral in formula (1) and formula (2) is measurediAnd replacing the density of the mineral to calculate the density of the rock skeleton model. Due to the elastic modulus of pure mineralsAnd density are generally invariant to the environment, so known laboratory measurements can also be used.
For the matrix pores, assuming that all shale gas is contained in the matrix pores, the gas phase can be contained in the rock skeleton by means of Kuster-Toksoz formula or self-compatible approximation (SCA) or Differential Equivalence (DEM) and the like, and the elastic modulus and the density of a single layer in the shale with a thin physical structure are obtained. This step thus results in an isotropic monolayer having an isotropic elastic stiffness tensor of formula (4). Each monolayer having an elastic stiffness tensor IC of the form described aboveijHowever, the elastic modulus and thus the elemental value may be different.
Then, regarding the single layer as a transverse isotropic (VTI) medium, the whole transverse isotropic medium is obtained by combining a plurality of single layers, and the elastic stiffness tensor of the VTI medium can be calculated to be a formula (5) by utilizing the elastic modulus and the density obtained by the previous step, wherein the operator < x > represents that the weighted average of the mineral component x of the rock in the bracket < > is carried out according to the volume content, namely the Backus average.
And adding approximately vertical fractures into the shale reservoir (background medium) with the thin physical structure to form a shale reservoir model. Still assuming that all the fractures of the shale gas reservoir contain shale gas, after parameters such as fracture density, fracture radius, fracture aspect ratio and the like are given, the elastic stiffness tensor of the shale gas reservoir can be obtained by using a Hudson fracture model or a Schoenberg fracture model as a formula (6).
And (3) calculating the shear wave velocity of the shale reservoir model as formula (7) according to the elastic stiffness tensor and the logging density of the shale reservoir model.
According to the method, the shale reservoir transverse wave velocity is quickly, simply and conveniently calculated through element logging and rock physics, and necessary parameters are provided for shale gas reservoir prediction, brittleness evaluation and fracture prediction.
Application example
To facilitate understanding of the solution of the embodiments of the present invention and the effects thereof, a specific application example is given below. It will be understood by those skilled in the art that this example is merely for the purpose of facilitating an understanding of the present invention and that any specific details thereof are not intended to limit the invention in any way.
The shale reservoir shear wave velocity calculation method comprises the following steps:
establishing a rock skeleton model according to mineral components and volume content of rocks obtained by element logging; calculating the Voigt average modulus and the Reuss average modulus of the rock skeleton model respectively through a formula (1) and a formula (2), wherein the Voigt average modulus is as follows:
wherein M isVIs the Voigt mean modulus, fiAnd MiThe volume content and the elastic modulus of the ith mineral component are respectively, i is 1, 2, …, N, and the average modulus of Reuss is:
Figure BDA0001742951260000122
wherein M isRCalculating the elastic modulus of the rock skeleton model according to the formula (3) as the average modulus of Reuss:
Figure BDA0001742951260000123
wherein M is the elastic modulus of the rock skeleton model, and the density is obtained by conventional well logging.
And (3) calculating the elastic stiffness tensor of the single-layer isotropy of the rock skeleton model as a formula (4) according to the elastic modulus and the density of the rock skeleton model:
Figure BDA0001742951260000124
wherein, ICijThe elastic stiffness tensor is a single-layer isotropic elastic stiffness tensor, wherein a, b, c, d, f and m are parameters of the single-layer isotropic elastic stiffness tensor, a is c and lambda +2 mu, b is f and lambda, and d is m and mu.
Step 103, calculating the transverse isotropic elastic stiffness tensor of the rock skeleton model as a formula (5) according to the single-layer isotropic elastic stiffness tensor of the rock skeleton model:
Figure BDA0001742951260000131
wherein, VCijIs the transversely isotropic elastic stiffness tensor, A, B, C, D, F, M is a parameter of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000132
Figure BDA0001742951260000133
M＝<m>wherein the operator<x>Represents a pair of brackets<>The mineral composition x of the rock in (a) is weighted averaged by volume content.
Adding fractures into a rock skeleton model to obtain a shale reservoir model, and calculating the elastic stiffness tensor of the shale reservoir model as a formula (6):
wherein, TCijElastic stiffness tensor, Delta, for shale reservoir modelNIn order to provide rigidity for crack initiation,
Figure BDA0001742951260000135
ΔHin order to provide the horizontal stiffness of the crack,g is the parametric ratio of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000137
and e represents the crack density.
According to the elastic stiffness tensor and the logging density of the shale reservoir model, calculating the shear wave velocity of the shale reservoir model as a formula (7):
Figure BDA0001742951260000138
wherein v issThe shear wave velocity is denoted by ρ as the density.
In conclusion, the shale gas reservoir transverse wave velocity is quickly, simply and conveniently calculated through element logging and rock physics, and necessary parameters are provided for shale gas reservoir prediction, brittleness evaluation and fracture prediction.
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
The shale reservoir shear wave velocity calculation system is characterized by comprising the following components: a memory storing computer-executable instructions; a processor executing computer executable instructions in the memory to perform the steps of: establishing a rock skeleton model, and calculating the elastic modulus and the density of the rock skeleton model; calculating the elastic stiffness tensor of the single-layer isotropy of the rock skeleton model according to the elastic modulus and the density of the rock skeleton model; calculating the transverse isotropic elastic stiffness tensor of the rock skeleton model according to the single-layer isotropic elastic stiffness tensor of the rock skeleton model; adding fractures into the rock skeleton model to obtain a shale reservoir model, and calculating the elastic stiffness tensor of the shale reservoir model; and calculating the shear wave velocity of the shale reservoir model according to the elastic stiffness tensor and the logging density of the shale reservoir model.
In one example, building a rock skeleton model, calculating an elastic modulus of the rock skeleton model comprises: establishing a rock skeleton model according to mineral components and volume content of rocks obtained by element logging; and calculating the Voigt average modulus and the Reuss average modulus of the rock skeleton model, and further calculating the elastic modulus of the rock skeleton model.
In one example, the Voigt average modulus is:
Figure BDA0001742951260000141
wherein M isVIs the Voigt mean modulus, fiAnd MiThe volume content and the elastic modulus of the ith mineral component, i ═ 1, 2, …, N, respectively.
In one example, the reus mean modulus is:
Figure BDA0001742951260000151
wherein M isRIs the reus mean modulus.
In one example, the modulus of elasticity of the rock skeleton model is:
wherein M is the elastic modulus of the rock skeleton model.
In one example, the elastic stiffness tensor for the single layer isotropy is:
Figure BDA0001742951260000153
wherein, ICijThe elastic stiffness tensor is a single-layer isotropic elastic stiffness tensor, wherein a, b, c, d, f and m are parameters of the single-layer isotropic elastic stiffness tensor, a is c and lambda +2 mu, b is f and lambda, and d is m and mu.
In one example, the transversely isotropic elastic stiffness tensor is:
wherein, VCijIs the transversely isotropic elastic stiffness tensor, A, B, C, D, F, M is a parameter of the transversely isotropic elastic stiffness tensor,
Figure BDA0001742951260000155
Figure BDA0001742951260000156
M＝<m>wherein the operator<x>Represents a pair of brackets<>The mineral composition x of the rock in (a) is weighted averaged by volume content.
In one example, the elastic stiffness tensor of the shale reservoir model is:
Figure BDA0001742951260000161
wherein, TCijElastic stiffness tensor, Delta, for shale reservoir modelNIn order to provide rigidity for crack initiation,
Figure BDA0001742951260000162
ΔHin order to provide the horizontal stiffness of the crack,
Figure BDA0001742951260000163
g is the parametric ratio of the transversely isotropic elastic stiffness tensor,and e represents the crack density.
In one example, the shear wave velocity of the shale reservoir model is:
Figure BDA0001742951260000165
wherein v issThe shear wave velocity is denoted by ρ as the density.
The system can quickly and simply calculate the shale reservoir transverse wave velocity through element well logging and rock physics, and provides necessary parameters for shale gas reservoir prediction, brittleness evaluation and fracture prediction.
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 shale reservoir shear wave velocity calculation method is characterized by comprising the following steps:
establishing a rock skeleton model, and calculating the elastic modulus and the density of the rock skeleton model;
calculating the elastic stiffness tensor of the single-layer isotropy of the rock skeleton model according to the elastic modulus and the density of the rock skeleton model;
calculating the transverse isotropic elastic stiffness tensor of the rock skeleton model according to the single-layer isotropic elastic stiffness tensor of the rock skeleton model;
adding fractures into the rock skeleton model to obtain a shale reservoir model, and calculating the elastic stiffness tensor of the shale reservoir model;
and calculating the shear wave velocity of the shale reservoir model according to the elastic stiffness tensor and the logging density of the shale reservoir model.
2. The shale reservoir shear wave velocity calculation method of claim 1, wherein building a rock skeleton model, calculating an elastic modulus of the rock skeleton model comprises:
establishing the rock skeleton model according to mineral components and volume content of rock obtained by element logging;
and calculating the Voigt average modulus and the Reuss average modulus of the rock skeleton model, and further calculating the elastic modulus of the rock skeleton model.
3. The shale reservoir shear wave velocity calculation method of claim 2, wherein the Voigt average modulus is:
Figure FDA0001742951250000011
wherein M isVIs the Voigt mean modulus, fiAnd MiThe volume content and the elasticity of the ith mineral componentAnd (3) the modulus of elasticity, i ═ 1, 2, …, N.
4. The shale reservoir shear wave velocity calculation method of claim 3, wherein the reus average modulus is:
Figure FDA0001742951250000021
wherein M isRIs the reus mean modulus.
5. The shale reservoir shear wave velocity calculation method of claim 4, wherein the modulus of elasticity of the rock skeleton model is:
Figure FDA0001742951250000022
wherein M is the elastic modulus of the rock skeleton model.
6. The shale reservoir shear wave velocity calculation method of claim 1, wherein the elastic stiffness tensor of the single layer isotropy is:
Figure FDA0001742951250000023
wherein, ICijThe elastic stiffness tensor is a single-layer isotropic elastic stiffness tensor, wherein a, b, c, d, f and m are parameters of the single-layer isotropic elastic stiffness tensor, a is c and lambda +2 mu, b is f and lambda, and d is m and mu.
7. The shale reservoir shear wave velocity calculation method of claim 6, wherein the transversely isotropic elastic stiffness tensor is:
Figure FDA0001742951250000031
wherein, VCijBeing transversely isotropicElastic stiffness tensor A, B, C, D, F, M is a parameter of the transverse isotropic elastic stiffness tensor, where A is<a-f2·c-1>+<c-1>-1·<f·c-1>2，B＝<b-f2·c-1>+<c-1>-1·<f·c-1>2，C＝<c-1>-1，F＝<c-1>-1·<f·c-1>，D＝<d-1>-1，M＝<m>Wherein the operator<x>Represents a pair of brackets<>The mineral composition x of the rock in (a) is weighted averaged by volume content.
8. The shale reservoir shear wave velocity calculation method of claim 7, wherein the elastic stiffness tensor of the shale reservoir model is:
Figure FDA0001742951250000032
wherein, TCijElastic stiffness tensor, Delta, for shale reservoir modelNIn order to provide rigidity for crack initiation,
Figure FDA0001742951250000033
ΔHin order to provide the horizontal stiffness of the crack,
Figure FDA0001742951250000034
g is the parametric ratio of the transversely isotropic elastic stiffness tensor,
Figure FDA0001742951250000035
and e represents the crack density.
9. The shale reservoir shear wave velocity calculation method of claim 7, wherein the shear wave velocity of the shale reservoir model is:
Figure FDA0001742951250000041
wherein v issThe shear wave velocity is denoted by ρ as the density.
10. A shale reservoir shear wave velocity calculation system, the system comprising:
a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
establishing a rock skeleton model, and calculating the elastic modulus and the density of the rock skeleton model;
calculating the elastic stiffness tensor of the single-layer isotropy of the rock skeleton model according to the elastic modulus and the density of the rock skeleton model;
calculating the transverse isotropic elastic stiffness tensor of the rock skeleton model according to the single-layer isotropic elastic stiffness tensor of the rock skeleton model;
adding fractures into the rock skeleton model to obtain a shale reservoir model, and calculating the elastic stiffness tensor of the shale reservoir model;
and calculating the shear wave velocity of the shale reservoir model according to the elastic stiffness tensor and the logging density of the shale reservoir model.
CN201810828120.9A 2018-07-25 2018-07-25 Shale reservoir shear wave velocity calculation method and system Pending CN110764144A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810828120.9A CN110764144A (en) 2018-07-25 2018-07-25 Shale reservoir shear wave velocity calculation method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810828120.9A CN110764144A (en) 2018-07-25 2018-07-25 Shale reservoir shear wave velocity calculation method and system

Publications (1)

Publication Number Publication Date
CN110764144A true CN110764144A (en) 2020-02-07

Family

ID=69327246

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810828120.9A Pending CN110764144A (en) 2018-07-25 2018-07-25 Shale reservoir shear wave velocity calculation method and system

Country Status (1)

Country Link
CN (1) CN110764144A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113705868A (en) * 2021-08-16 2021-11-26 长江大学 Shale gas state quantitative prediction method, system, equipment and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105095631A (en) * 2014-05-21 2015-11-25 中国石油化工股份有限公司 Shale anisotropic rock physical modeling method
CN105203739A (en) * 2015-10-29 2015-12-30 中国石油大学(北京) Rock physical model-based shale stratum anisotropy parameter prediction method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105095631A (en) * 2014-05-21 2015-11-25 中国石油化工股份有限公司 Shale anisotropic rock physical modeling method
CN105203739A (en) * 2015-10-29 2015-12-30 中国石油大学(北京) Rock physical model-based shale stratum anisotropy parameter prediction method

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
M . SCHOENBERG ET AL: "ELASTIC WAVE PROPAGATION IN MEDIA WITH PARALLEL FRACTURES AND ALIGNED CRACKS", 《GEOPHYSICAL PROSPECTING》 *
XINPENG PAN ET AL.: "Estimation of effective geostress parameters driven by anisotropic stress and rock physics models with orthorhombic symmetry", 《JOURNAL OF GEOPHYSICS AND ENGINEERING》 *
ZHIQI GUO ET AL.: "Seismic rock physics characterization of anisotropic shale—a Longmaxi Shale case study", 《JOURNAL OF GEOPHYSICS AND ENGINEERING》 *
刘百红等: "HTI介质的方位AVO正演研究", 《石油物探》 *
张广智等: "基于页岩岩石物理等效模型的地应力预测方法研究", 《地球物理学报》 *
杨展等: "《地理学大辞典》", 31 August 1992, 安徽人民出版社 *
汪云甲等: "《基于多传感器及多元监测数据的瓦斯预警理论与方法研究》", 31 May 2017, 徐州：中国矿业大学出版社 *
董宁等: "《页岩气地震勘探技术》", 31 December 2016, 上海：华东理工大学出版社 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113705868A (en) * 2021-08-16 2021-11-26 长江大学 Shale gas state quantitative prediction method, system, equipment and storage medium
CN113705868B (en) * 2021-08-16 2024-02-09 长江大学 Shale gas state quantitative prediction method, system, equipment and storage medium

Similar Documents

Publication Publication Date Title
Wyllie et al. Elastic wave velocities in heterogeneous and porous media
CA2580570C (en) Integrated anisotropic rock physics model
Kuster et al. Velocity and attenuation of seismic waves in two-phase media: Part I. Theoretical formulations
US8498853B2 (en) Petrophysical method for predicting plastic mechanical properties in rock formations
Jeppson et al. Geophysical properties within the San Andreas Fault Zone at the San Andreas Fault Observatory at Depth and their relationships to rock properties and fault zone structure
Ali et al. 3D static modeling and petrographic aspects of the Albian/Cenomanian Reservoir, Komombo Basin, Upper Egypt
Xie et al. Construction of a 3D geomechanical model for development of a shale gas reservoir in the Sichuan Basin
CN104007482B (en) A kind of mud shale petrophysical model method based on anisotropy effective field
CN104636588A (en) Method and device for calculating content of organic carbon in source rock
Ali et al. 3D structural, facies and petrophysical modeling of C member of six hills formation, Komombo Basin, Upper Egypt
CN110851937B (en) Shale reservoir brittleness index calculation method and system
Amjad et al. Sembar formation as an unconventional prospect: new insights in evaluating shale gas potential combined with deep learning
Al-Qayim et al. Reservoir characterization of an intra-orogenic carbonates platform: Pila Spi Formation, Taq Taq oil field, Kurdistan, Iraq
CN110764144A (en) Shale reservoir shear wave velocity calculation method and system
Porlles et al. Integrated porosity methods for estimation of gas-in-place in the Muerto Formation of Northwestern Peru
Franquet et al. Integrated acoustic, mineralogy, and geomechanics characterization of the Huron shale, Southern West Virginia, USA
CN116794746A (en) Geophysical quantitative prediction method for sea shale lithofacies
McCormack et al. Discrete measurements of the least horizontal principal stress from core data: An application of viscoelastic stress relaxation
Kadhim et al. Correlation between cementation factor and carbonate reservoir rock properties
Qi et al. Rock physics modeling and well-log practice for unconventional reservoirs
Perez Application of LMR inversion and clustering analysis in the Barnett Shale
Reyer Outcrop analogue studies of rocks from the Northwest German Basin for geothermal exploration and exploitation: fault zone structure, heterogeneous rock properties, and application to reservoir conditions
Lobo et al. Methodology for Petrophysical and Geomechanical Analysis of Shale Plays. Case Study: La Luna and Capacho Formations, Maracaibo Basin.
Narasimhan et al. Constructing High Resolution, Inch Scale Continuous Logs via a Multi Domain Approach to Improve Hydraulic Fracturing by Capturing Thin Beds in Bone Spring, Delaware Basin, Reeves County, TX
Ma et al. Petroleum Science

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200207