CN112649854A - Seismic wave frequency dispersion and attenuation prediction method and system based on dual-scale model - Google Patents
Seismic wave frequency dispersion and attenuation prediction method and system based on dual-scale model Download PDFInfo
- Publication number
- CN112649854A CN112649854A CN201910963929.7A CN201910963929A CN112649854A CN 112649854 A CN112649854 A CN 112649854A CN 201910963929 A CN201910963929 A CN 201910963929A CN 112649854 A CN112649854 A CN 112649854A
- Authority
- CN
- China
- Prior art keywords
- rock
- modulus
- frequency
- seismic wave
- 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
Links
- 239000006185 dispersion Substances 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 22
- 239000011435 rock Substances 0.000 claims abstract description 128
- 239000012530 fluid Substances 0.000 claims abstract description 80
- 239000011148 porous material Substances 0.000 claims abstract description 66
- 230000001419 dependent effect Effects 0.000 claims abstract description 20
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 10
- 239000011707 mineral Substances 0.000 claims abstract description 10
- 229920006395 saturated elastomer Polymers 0.000 claims description 16
- 239000002245 particle Substances 0.000 claims description 12
- 239000007787 solid Substances 0.000 claims description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 230000035699 permeability Effects 0.000 claims description 7
- 239000003570 air Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 238000009736 wetting Methods 0.000 claims description 6
- 230000009977 dual effect Effects 0.000 claims 8
- 238000002834 transmittance Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
A seismic wave frequency dispersion and attenuation prediction method and system based on a dual-scale model are disclosed. The method can comprise the following steps: calculating the bulk modulus and the shear modulus of the rock frequency-changing wet framework according to the bulk modulus and the shear modulus of rock minerals and framework pore fluid; calculating the frequency-dependent bulk modulus of the pore fluid; calculating the effective volume modulus and the effective shear modulus of the rock according to the frequency-variable volume modulus of the pore fluid and the volume modulus and the shear modulus of the rock frequency-variable wet skeleton; and calculating the seismic wave velocity and the attenuation factor according to the effective volume modulus and the effective shear modulus of the rock. The method utilizes the characteristics of high-frequency attenuation and low-frequency enhancement of oil gas, simultaneously reflects the attenuation information of the seismic frequency band and the ultrasonic frequency band, improves the accuracy of prediction, and provides more comprehensive explanation for the reasons of velocity dispersion and attenuation of elastic waves in a full frequency band range.
Description
Technical Field
The invention relates to the field of petroleum exploration and development, in particular to a seismic wave frequency dispersion and attenuation prediction method and system based on a dual-scale model.
Background
To study seismic wave dispersion and attenuation, understanding the attenuation mechanism of seismic waves is the first and several relatively mature attenuation mechanisms have been developed.
Jet Flow (also called local Flow) is usually generated by pore fluid when seismic waves propagate in a subsurface reservoir, which is considered by the scholars to be the main factor causing seismic wave velocity dispersion, and Mavko-Jizba (1974) proposes a jet Flow (Squirt Flow) model. Seismic wave propagation can cause pore pressure gradients in the rock pores, with different pores in the rock having different compressibilities. When seismic waves propagate in a saturated fluid porous medium, pore fluid in the soft pores flows to the hard pores; when seismic waves are transmitted in a saturated part of gas porous medium, the existence of pore fluid pressure can enable fluid with poor compression capacity to flow to gas with strong compression capacity, and the two flows are jet flows; since the model is only suitable for high-frequency conditions, Gurevich (2009) optimizes the model for the limitation, and establishes a high-frequency limit jet flow model which can be suitable for saturated gas rock pore space, fluid and dry rock. The model only reflects the phenomena of seismic wave velocity dispersion and energy attenuation in an ultrasonic frequency band caused by rock skeleton saturated fluid.
White (1975a) and the like establish a periodic spherical plaque saturation model for the first time. The model assumes that the fluid in each unit cube is distributed in the shape of water-in-air spot blocks, and the stress on the fluid is uniformly distributed. White (1975b) and the like establish a periodic layered plaque saturation model on the basis of a spherical plaque saturation model, the White model is only applicable to the condition of fluid pressure discontinuity between saturated two different fluid stratum interfaces, and Dutta (1979) solves the White spherical plaque saturation model again on the basis of the Biot theory so as to be applicable to the condition of continuous fluid pressure at different fluid saturated stratum interfaces; johnson D L (2001) further optimizes the model so that the model is suitable for various patch shapes, and Muller T M (2009) studies the relationship between physical parameters and seismic wave velocity and attenuation factors. The model only reflects the phenomena of velocity dispersion and capability attenuation of seismic waves in a seismic frequency band caused by mutual flow of plaque saturated fluids in rocks.
In actual measurement, the velocity dispersion and energy attenuation phenomena exist in the whole frequency band range, the two models can only show the dispersion and attenuation phenomena in a single frequency band, and the dispersion and attenuation phenomena in the whole seismic frequency band range cannot be reflected well at the same time, so that the simulation result and the actual result have great difference. Therefore, it is necessary to develop a seismic wave dispersion and attenuation prediction method and system based on a dual-scale model.
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 seismic wave frequency dispersion and attenuation prediction method and system based on a dual-scale model, which utilize the characteristics of high-frequency attenuation and low-frequency enhancement of oil gas, simultaneously embody the attenuation information of a seismic frequency band and an ultrasonic frequency band, improve the prediction accuracy and give a more comprehensive explanation to the reasons of the velocity frequency dispersion and the attenuation of elastic waves in a full-frequency band range.
According to one aspect of the invention, a seismic wave frequency dispersion and attenuation prediction method based on a dual-scale model is provided. The method may include: calculating the bulk modulus and the shear modulus of the rock frequency-changing wet framework according to the bulk modulus and the shear modulus of rock minerals and framework pore fluid; calculating the frequency-dependent bulk modulus of the pore fluid; calculating the effective volume modulus and the effective shear modulus of the rock according to the frequency-variable volume modulus of the pore fluid and the volume modulus and the shear modulus of the frequency-variable wet skeleton of the rock; and calculating the seismic wave velocity and the attenuation factor according to the effective volume modulus and the effective shear modulus of the rock.
Preferably, the bulk modulus of the rock frequency-change wet skeleton is calculated by formula (1):
calculating the shear modulus of the rock frequency-wetting framework by formula (2):
wherein, Kmf(omega) is the bulk modulus, mu, of the frequently wetted skeleton of the rockmf(omega) is the shear modulus of the frequently wetted framework of the rock, KhIs the bulk modulus, K, of the dry rock skeleton under high pressurem、μmRespectively the bulk modulus and shear modulus of the dry rock skeleton, KfIs the bulk modulus, K, of the pore fluidsIs the bulk modulus of the solid particles,
is the soft porosity of the skeleton.
Preferably, the frequency-dependent bulk modulus of the pore fluid is calculated by equation (3):
wherein the content of the first and second substances,
is the frequency-dependent bulk modulus of the pore fluid, k is the wavenumber of the pressure dissipation waves in the soft pores,
a is the soft pore radius, η is the fluid viscosity, and κ is the permeability of the saturated fluid rock.
Preferably, the rock effective bulk modulus is calculated by equation (4):
K(ω)＝Kmf(ω)+2α2(ω)D(ω) (4)
calculating the rock effective shear modulus by equation (5):
μ(ω)＝μmf(ω) (5)
wherein K (omega) is the effective bulk modulus of the rock, mu (omega) is the effective shear modulus of the rock, alpha (omega) and D (omega) are both calculation parameters,
preferably, the seismic wave velocity is calculated by equation (6):
calculating the attenuation factor by equation (7):
wherein Vp (omega) is the seismic wave velocity, 1/Q is the attenuation factor, rhos、ρw、ρg、ρfDensity, p, of solid particles, water, air, mixed fluid, respectivelyf＝ρwSw+ρgSg，
For rock hard porosity, Imag (Vp (ω)) is the imaginary part of the seismic wave velocity and Real (Vp (ω)) is the Real part of the seismic wave velocity.
According to another aspect of the invention, a seismic wave dispersion and attenuation prediction system based on a dual-scale model 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: calculating the bulk modulus and the shear modulus of the rock frequency-changing wet framework according to the bulk modulus and the shear modulus of rock minerals and framework pore fluid; calculating the frequency-dependent bulk modulus of the pore fluid; calculating the effective volume modulus and the effective shear modulus of the rock according to the frequency-variable volume modulus of the pore fluid and the volume modulus and the shear modulus of the frequency-variable wet skeleton of the rock; and calculating the seismic wave velocity and the attenuation factor according to the effective volume modulus and the effective shear modulus of the rock.
Preferably, the bulk modulus of the rock frequency-change wet skeleton is calculated by formula (1):
calculating the shear modulus of the rock frequency-wetting framework by formula (2):
wherein, Kmf(omega) is the bulk modulus, mu, of the frequently wetted skeleton of the rockmf(omega) is the shear modulus of the frequently wetted framework of the rock, KhIs the bulk modulus, K, of the dry rock skeleton under high pressurem、μmRespectively the bulk modulus and shear modulus of the dry rock skeleton, KfIs the bulk modulus, K, of the pore fluidsIs the bulk modulus of the solid particles,
is the soft porosity of the skeleton.
Preferably, the frequency-dependent bulk modulus of the pore fluid is calculated by equation (3):
wherein the content of the first and second substances,
is the frequency-dependent bulk modulus of the pore fluid, k is the wavenumber of the pressure dissipation waves in the soft pores,
a is the soft pore radius, η is the fluid viscosity, and κ is the permeability of the saturated fluid rock.
Preferably, the rock effective bulk modulus is calculated by equation (4):
K(ω)＝Kmf(ω)+2α2(ω)D(ω) (4)
calculating the rock effective shear modulus by equation (5):
μ(ω)＝μmf(ω) (5)
wherein K (omega) is the effective bulk modulus of the rock, mu (omega) is the effective shear modulus of the rock, alpha (omega) and D (omega) are both calculation parameters,
preferably, the seismic wave velocity is calculated by equation (6):
calculating the attenuation factor by equation (7):
wherein Vp (omega) is the seismic wave velocity, 1/Q is the attenuation factor, rhos、ρw、ρg、ρfAre respectively provided withDensity of solid particles, water, air and mixed fluid, rhof＝ρwSw+ρgSg，
For rock hard porosity, Imag (Vp (ω)) is the imaginary part of the seismic wave velocity and Real (Vp (ω)) is the Real part of the seismic wave velocity.
The method and apparatus of the present invention have 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 is a flow chart illustrating the steps of a seismic wave dispersion and attenuation prediction method based on a dual-scale model according to the invention.
Fig. 2a and 2b show a schematic diagram of the trend of the attenuation factor with frequency and the trend of the velocity with frequency, respectively, calculated from the jet model.
Fig. 3a and 3b show a diagram of the trend of the attenuation factor with frequency and the trend of the velocity with frequency, respectively, calculated from the plaque saturation model.
Fig. 4a and 4b show a diagram of the trend of the attenuation factor with frequency and the trend of the velocity with frequency, respectively, calculated by the dual-scale model 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 is a flow chart illustrating the steps of a seismic wave dispersion and attenuation prediction method based on a dual-scale model according to the invention.
In this embodiment, the seismic wave dispersion and attenuation prediction method based on the dual-scale model according to the present invention may include: step 101, calculating the bulk modulus and the shear modulus of a rock frequency-changing wet framework according to the bulk modulus and the shear modulus of rock minerals and framework pore fluid; 102, calculating the frequency-dependent bulk modulus of the pore fluid; 103, calculating the effective volume modulus and the effective shear modulus of the rock according to the frequency-variable volume modulus of the pore fluid and the volume modulus and the shear modulus of the rock frequency-variable wet skeleton; and step 104, calculating the seismic wave velocity and the attenuation factor according to the effective volume modulus and the effective shear modulus of the rock.
In one example, the bulk modulus of a rock frequency-changing wet skeleton is calculated by equation (1):
calculating the shear modulus of the rock frequency-wetting framework by the formula (2):
wherein, Kmf(omega) is the bulk modulus, mu, of the frequently wetted skeleton of the rockmf(omega) is the shear modulus of the frequently wetted framework of the rock, KhIs the bulk modulus, K, of the dry rock skeleton under high pressurem、μmRespectively the bulk modulus and shear modulus of the dry rock skeleton, KfIs the bulk modulus, K, of the pore fluidsIs the bulk modulus of the solid particles,
is the soft porosity of the skeleton.
In one example, the frequency-dependent bulk modulus of the pore fluid is calculated by equation (3):
wherein the content of the first and second substances,
is the frequency-dependent bulk modulus of the pore fluid, k is the wavenumber of the pressure dissipation waves in the soft pores,
a is the soft pore radius, η is the fluid viscosity, and κ is the permeability of the saturated fluid rock.
In one example, the effective bulk modulus of the rock is calculated by equation (4):
K(ω)＝Kmf(ω)+2α2(ω)D(ω) (4)
calculating the effective shear modulus of the rock by formula (5):
μ(ω)＝μmf(ω) (5)
wherein K (omega) is the effective bulk modulus of the rock, mu (omega) is the effective shear modulus of the rock, alpha (omega) and D (omega) are both calculation parameters,
in one example, the seismic wave velocity is calculated by equation (6):
the attenuation factor is calculated by equation (7):
wherein Vp (omega) is the seismic wave velocity, 1/Q is the attenuation factor, rhos、ρw、ρg、ρfDensity, p, of solid particles, water, air, mixed fluid, respectivelyf＝ρwSw+ρgSg，
For rock hard porosity, Imag (Vp (ω)) is the imaginary part of the seismic wave velocity and Real (Vp (ω)) is the Real part of the seismic wave velocity.
Specifically, in a porous medium saturated with fluid, two different flow modes of fluid relative flow and jet flow exist, which correspond to different scales and frequency band ranges respectively. In the range of seismic frequency bands, the main reason for causing the velocity dispersion and energy attenuation of seismic waves is that plaques for saturating two different fluids (water and air) are distributed between a rock matrix and pores, and the different fluids can generate relative flow, and the mechanism is a plaque saturation model and belongs to a mesoscopic scale; in the ultrasonic frequency range, the main reason for seismic wave velocity dispersion and energy attenuation is that soft pore fluid existing in the rock skeleton generates jet flow, namely a jet flow model, and belongs to the micro scale.
In fact, when the seismic wave propagates in the rock, both attenuation mechanisms exist, and both models only consider the attenuation mechanism in a single frequency band, namely, the velocity dispersion phenomenon and the energy attenuation phenomenon of the seismic wave cannot be comprehensively analyzed by only using a single Squirt model or a plaque saturation model. Therefore, the velocity dispersion and energy attenuation phenomena of the seismic waves in the full frequency band can be more comprehensively analyzed by establishing a double-scale model.
The seismic wave frequency dispersion and attenuation prediction method based on the dual-scale model can comprise the following steps:
according to the bulk modulus and shear modulus of rock minerals and framework pore fluid, a frequency-varying Squirt model is adopted, the bulk modulus of the rock frequency-varying wet framework is calculated through a formula (1), and the shear modulus of the rock frequency-varying wet framework is calculated through a formula (2), wherein the rock matrix minerals are quartz, and the framework pore fluid is water.
Calculating the frequency-dependent bulk modulus of the pore fluid by formula (3); substituting the frequency-variable volume modulus of pore fluid, the volume modulus and the shear modulus of the rock frequency-variable wet skeleton into a full-band White layered plaque model, calculating the effective volume modulus of the rock through a formula (4), and calculating the effective shear modulus of the rock through a formula (5).
And (3) calculating the seismic wave velocity through a formula (6) according to the effective volume modulus and the effective shear modulus of the rock, and further calculating the attenuation factor through a formula (7).
The change of the seismic wave velocity along with the frequency can be observed according to the image, namely the seismic wave frequency dispersion phenomenon is obtained; the attenuation factor changes with the frequency and has an extreme point, namely the frequency with the maximum attenuation degree
The method utilizes the characteristics of high-frequency attenuation and low-frequency enhancement of oil gas, simultaneously reflects the attenuation information of the seismic frequency band and the ultrasonic frequency band, improves the accuracy of prediction, and provides more comprehensive explanation for the reasons of velocity dispersion and attenuation of elastic waves in a full frequency band range.
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 invention is illustrated by comparing experimental measurement results and model calculation results of partially saturated (water, air) sanders:
the rock sample is water-gas partially saturated Fengdan white dew pure sandstone, stress strain of the rock sample under different frequencies (1-2000Hz) is measured through DAS facilities, the stress strain is converted into seismic wave velocity and attenuation factors under different frequencies according to a rock physical formula, and actually measured data are shown as black dots in figures 2a-b, 3a-b and 4 a-b.
The physical property parameters of the mineral components of the rock sample, and the parameters such as porosity and permeability of the rock sample are specifically shown in table 1.
TABLE 1
Fig. 2a and 2b show a schematic diagram of the trend of the attenuation factor with frequency and the trend of the velocity with frequency, respectively, calculated from the jet model. The model only reflects the attenuation information of the ultrasonic frequency band, the numerical value of the attenuation factor is obviously larger than that of the measured data, and the seismic wave speed is consistent with the measured data at the low-frequency limit and the high-frequency band.
Fig. 3a and 3b show a diagram of the trend of the attenuation factor with frequency and the trend of the velocity with frequency, respectively, calculated from the plaque saturation model. The model only reflects attenuation information of a seismic frequency band, the data goodness of fit is high, the seismic wave speed is matched with measured data in a low frequency band, but the high frequency band is far smaller than the measured data.
Fig. 4a and 4b show a diagram of the trend of the attenuation factor with frequency and the trend of the velocity with frequency, respectively, calculated by the dual-scale model according to the present invention. The model can simultaneously reflect the attenuation information of the seismic frequency band and the ultrasonic frequency band, and the seismic wave speed can be well matched with the actually measured data in the ranges of the seismic frequency band and the ultrasonic frequency band.
In conclusion, the method utilizes the characteristics of high-frequency attenuation and low-frequency enhancement of oil gas, simultaneously reflects the attenuation information of the seismic frequency band and the ultrasonic frequency band, improves the accuracy of prediction, and provides a more comprehensive explanation for the reasons of velocity dispersion and attenuation of elastic waves in a full frequency band range.
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.
According to an embodiment of the invention, a seismic wave frequency dispersion and attenuation prediction system based on a dual-scale model 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: calculating the bulk modulus and the shear modulus of the rock frequency-changing wet framework according to the bulk modulus and the shear modulus of rock minerals and framework pore fluid; calculating the frequency-dependent bulk modulus of the pore fluid; calculating the effective volume modulus and the effective shear modulus of the rock according to the frequency-variable volume modulus of the pore fluid and the volume modulus and the shear modulus of the rock frequency-variable wet skeleton; and calculating the seismic wave velocity and the attenuation factor according to the effective volume modulus and the effective shear modulus of the rock.
In one example, the bulk modulus of a rock frequency-changing wet skeleton is calculated by equation (1):
calculating the shear modulus of the rock frequency-wetting framework by the formula (2):
wherein, Kmf(omega) is the bulk modulus, mu, of the frequently wetted skeleton of the rockmf(omega) is the shear modulus of the frequently wetted framework of the rock, KhIs the bulk modulus, K, of the dry rock skeleton under high pressurem、μmRespectively the bulk modulus and shear modulus of the dry rock skeleton, KfIs the bulk modulus, K, of the pore fluidsIs the bulk modulus of the solid particles,
is the soft porosity of the skeleton.
In one example, the frequency-dependent bulk modulus of the pore fluid is calculated by equation (3):
wherein the content of the first and second substances,
is the frequency-dependent bulk modulus of the pore fluid, k is the wavenumber of the pressure dissipation waves in the soft pores,
a is the soft pore radius, η is the fluid viscosity, and κ is the permeability of the saturated fluid rock.
In one example, the effective bulk modulus of the rock is calculated by equation (4):
K(ω)＝Kmf(ω)+2α2(ω)D(ω) (4)
calculating the effective shear modulus of the rock by formula (5):
μ(ω)＝μmf(ω) (5)
wherein K (omega) is the effective bulk modulus of the rock, mu (omega) is the effective shear modulus of the rock, alpha (omega) and D (omega) are both calculation parameters,
in one example, the seismic wave velocity is calculated by equation (6):
the attenuation factor is calculated by equation (7):
wherein Vp (omega) is the seismic wave velocity, 1/Q is the attenuation factor, rhos、ρw、ρg、ρfDensity, p, of solid particles, water, air, mixed fluid, respectivelyf＝ρwSw+ρgSg，
For rock hard porosity, Imag (Vp (ω)) is the imaginary part of the seismic wave velocity and Real (Vp (ω)) is the Real part of the seismic wave velocity.
The system utilizes the characteristics of high-frequency attenuation and low-frequency enhancement of oil gas, simultaneously reflects the attenuation information of seismic frequency bands and ultrasonic frequency bands, improves the accuracy of prediction, and provides more comprehensive explanation for the reasons of velocity dispersion and attenuation of elastic waves in a full frequency band range.
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.
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 seismic wave frequency dispersion and attenuation prediction method based on a dual-scale model is characterized by comprising the following steps:
calculating the bulk modulus and the shear modulus of the rock frequency-changing wet framework according to the bulk modulus and the shear modulus of rock minerals and framework pore fluid;
calculating the frequency-dependent bulk modulus of the pore fluid;
calculating the effective volume modulus and the effective shear modulus of the rock according to the frequency-variable volume modulus of the pore fluid and the volume modulus and the shear modulus of the frequency-variable wet skeleton of the rock;
and calculating the seismic wave velocity and the attenuation factor according to the effective volume modulus and the effective shear modulus of the rock.
2. The dual scale model-based seismic wave dispersion and attenuation prediction method of claim 1, wherein the bulk modulus of the rock frequency-wet framework is calculated by formula (1):
calculating the shear modulus of the rock frequency-wetting framework by formula (2):
wherein, Kmf(omega) is the bulk modulus, mu, of the frequently wetted skeleton of the rockmf(omega) is the shear modulus of the frequently wetted framework of the rock, KhIs the bulk modulus, K, of the dry rock skeleton under high pressurem、μmRespectively the bulk modulus and shear modulus of the dry rock skeleton, KfIs the bulk modulus, K, of the pore fluidsIs the bulk modulus of the solid particles,
is the soft porosity of the skeleton.
3. The dual scale model-based seismic wave dispersion and attenuation prediction method of claim 1, wherein the frequency-dependent bulk modulus of the pore fluid is calculated by equation (3):
wherein the content of the first and second substances,
is the frequency-dependent bulk modulus of the pore fluid, k is the wavenumber of the pressure dissipation waves in the soft pores,
a is the soft pore radius, η is the fluid viscosity, and κ is the permeability of the saturated fluid rockAnd (4) the transmittance.
4. The dual scale model-based seismic wave dispersion and attenuation prediction method of claim 1, wherein the effective rock bulk modulus is calculated by equation (4):
K(ω)＝Kmf(ω)+2α2(ω)D(ω) (4)
calculating the rock effective shear modulus by equation (5):
μ(ω)＝μmf(ω) (5)
wherein K (omega) is the effective bulk modulus of the rock, mu (omega) is the effective shear modulus of the rock, alpha (omega) and D (omega) are both calculation parameters,
5. the dual scale model-based seismic wave dispersion and attenuation prediction method of claim 1, wherein the seismic wave velocity is calculated by equation (6):
calculating the attenuation factor by equation (7):
wherein Vp (omega) is the seismic wave velocity, 1/Q is the attenuation factor, rhos、ρw、ρg、ρfDensity, p, of solid particles, water, air, mixed fluid, respectivelyf＝ρwSw+ρgSg，
For rock hard porosity, Imag (Vp (ω)) is the imaginary part of the seismic velocityReal (Vp (ω)) is the Real part of the seismic wave velocity.
6. A seismic wave frequency dispersion and attenuation prediction system based on a dual-scale model is characterized by comprising:
a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
calculating the bulk modulus and the shear modulus of the rock frequency-changing wet framework according to the bulk modulus and the shear modulus of rock minerals and framework pore fluid;
calculating the frequency-dependent bulk modulus of the pore fluid;
calculating the effective volume modulus and the effective shear modulus of the rock according to the frequency-variable volume modulus of the pore fluid and the volume modulus and the shear modulus of the frequency-variable wet skeleton of the rock;
and calculating the seismic wave velocity and the attenuation factor according to the effective volume modulus and the effective shear modulus of the rock.
7. The dual scale model-based seismic wave dispersion and attenuation prediction system of claim 6, wherein the bulk modulus of the rock frequency-wet framework is calculated by equation (1):
calculating the shear modulus of the rock frequency-wetting framework by formula (2):
wherein, Kmf(omega) is the bulk modulus, mu, of the frequently wetted skeleton of the rockmf(omega) is the shear modulus of the frequently wetted framework of the rock, KhIs the bulk modulus, K, of the dry rock skeleton under high pressurem、μmRespectively the bulk modulus and shear modulus of the dry rock skeleton, KfIs the bulk modulus, K, of the pore fluidsIs the bulk modulus of the solid particles,
is the soft porosity of the skeleton.
8. The dual scale model-based seismic wave dispersion and attenuation prediction system of claim 6, wherein the frequency-dependent bulk modulus of the pore fluid is calculated by equation (3):
wherein the content of the first and second substances,
is the frequency-dependent bulk modulus of the pore fluid, k is the wavenumber of the pressure dissipation waves in the soft pores,
a is the soft pore radius, η is the fluid viscosity, and κ is the permeability of the saturated fluid rock.
9. The dual scale model-based seismic wave dispersion and attenuation prediction system of claim 6, wherein the effective rock bulk modulus is calculated by equation (4):
K(ω)＝Kmf(ω)+2α2(ω)D(ω) (4)
calculating the rock effective shear modulus by equation (5):
μ(ω)＝μmf(ω) (5)
wherein K (omega) is the effective bulk modulus of the rock, mu (omega) is the effective shear modulus of the rock, alpha (omega) and D (omega) are both calculation parameters,
10. the dual scale model-based seismic wave dispersion and attenuation prediction system of claim 6, wherein the seismic wave velocity is calculated by equation (6):
calculating the attenuation factor by equation (7):
wherein Vp (omega) is the seismic wave velocity, 1/Q is the attenuation factor, rhos、ρw、ρg、ρfDensity, p, of solid particles, water, air, mixed fluid, respectivelyf＝ρwSw+ρgSg，
For rock hard porosity, Imag (Vp (ω)) is the imaginary part of the seismic wave velocity and Real (Vp (ω)) is the Real part of the seismic wave velocity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910963929.7A CN112649854A (en) | 2019-10-11 | 2019-10-11 | Seismic wave frequency dispersion and attenuation prediction method and system based on dual-scale model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910963929.7A CN112649854A (en) | 2019-10-11 | 2019-10-11 | Seismic wave frequency dispersion and attenuation prediction method and system based on dual-scale model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112649854A true CN112649854A (en) | 2021-04-13 |
Family
ID=75343410
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910963929.7A Pending CN112649854A (en) | 2019-10-11 | 2019-10-11 | Seismic wave frequency dispersion and attenuation prediction method and system based on dual-scale model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112649854A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114236609A (en) * | 2021-12-17 | 2022-03-25 | 河海大学 | Prediction method for longitudinal wave velocity and attenuation of partially saturated hole fracture medium |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110058310A (en) * | 2019-04-16 | 2019-07-26 | 中国石油大学(华东) | The seimic wave velocity frequency dispersion and decaying prediction technique and device of earthquake rock |
-
2019
- 2019-10-11 CN CN201910963929.7A patent/CN112649854A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110058310A (en) * | 2019-04-16 | 2019-07-26 | 中国石油大学(华东) | The seimic wave velocity frequency dispersion and decaying prediction technique and device of earthquake rock |
Non-Patent Citations (8)
| Title |
|---|
| J. GERM´AN RUBINO ET AL.: "Research note: Seismic attenuation due to wave-induced fluid flow at microscopic and mesoscopic scales", 《GEOPHYSICAL PROSPECTING》 * |
| J. GERM´AN RUBINO ET AL.: "Research note: Seismic attenuation due to wave-induced fluid flow at microscopic and mesoscopic scales", 《GEOPHYSICAL PROSPECTING》, 31 December 2013 (2013-12-31), pages 882 * |
| YAN-XIAO HE ET AL.: "Modeling the Effect of Microscopic and Mesoscopic Heterogeneities on Frequency-Dependent Attenuation and Seismic Signatures", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
| YAN-XIAO HE ET AL.: "Modeling the Effect of Microscopic and Mesoscopic Heterogeneities on Frequency-Dependent Attenuation and Seismic Signatures", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》, vol. 15, no. 8, 31 August 2018 (2018-08-31), pages 1174 - 1175, XP011687741, DOI: 10.1109/LGRS.2018.2829837 * |
| 张广智等: "微观与介观波致流下的速度频散与衰减", 《石油地球物理勘探》 * |
| 张广智等: "微观与介观波致流下的速度频散与衰减", 《石油地球物理勘探》, vol. 52, no. 4, 31 August 2017 (2017-08-31), pages 745 - 746 * |
| 李龙: "碎屑岩致密储层岩石物理模量定量表征", 中国优秀博硕士学位论文全文数据库（硕士）, pages 24 - 25 * |
| 王琪: "考虑孔隙流体影响的岩石物理模型研究", 中国优秀博硕士学位论文全文数据库（硕士）, pages 58 - 64 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114236609A (en) * | 2021-12-17 | 2022-03-25 | 河海大学 | Prediction method for longitudinal wave velocity and attenuation of partially saturated hole fracture medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11098565B2 (en) | Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves | |
| Carcione et al. | Computational poroelasticity—A review | |
| CN108957542B (en) | Method for establishing seismic wave attenuation rock physical drawing board | |
| Yew et al. | The determination of biot's parameters for sandstones: Part 1: Static tests Experimental procedure for determining Biot's parameters for four rock specimens is presented and theoretical values of wave speeds propagating through them are compared with experimental values | |
| Sidler | A porosity-based Biot model for acoustic waves in snow | |
| CN109298443B (en) | Method for evaluating elastic wave of double-fluid jet flow model and computer readable storage medium | |
| Matuszyk et al. | Solution of coupled poroelastic/acoustic/elastic wave propagation problems using automatic hp-adaptivity | |
| CN114236609B (en) | Prediction method for longitudinal wave velocity and attenuation of partially saturated hole fracture medium | |
| Santos et al. | Numerical simulation in applied geophysics | |
| Ba et al. | Velocity field of wave-induced local fluid flow in double-porosity media | |
| CN112649854A (en) | Seismic wave frequency dispersion and attenuation prediction method and system based on dual-scale model | |
| Yu et al. | Rock physics modeling of heterogeneous carbonate reservoirs: porosity estimation and hydrocarbon detection | |
| Mesgouez et al. | Transient solution for multilayered poroviscoelastic media obtained by an exact stiffness matrix formulation | |
| Wuttke et al. | Influence of poroelasticity on the 3D seismic response of complex geological media | |
| Boadu | Predicting oil saturation from velocities using petrophysical models and artificial neural networks | |
| Gao et al. | Seismic wave equations in tight oil/gas sandstone media | |
| Santos et al. | A parametric analysis of waves propagating in a porous solid saturated by a three-phase fluid | |
| CN110967746B (en) | Fluid saturation seismic inversion method and system | |
| Barbosa et al. | Hybrid modeling of 3D fractured media based on 2D numerical simulations: Aligned fractures case | |
| Zhang et al. | Frequency-dependent energy attenuation and velocity dispersion in periodic layered media | |
| CN112649861A (en) | Fractured reservoir rock physical modeling method and system | |
| Bezemer-Krijnen et al. | Predicting sound absorption of stacked spheres: Combining an analytical and numerical approach | |
| CN111597763B (en) | Cross multi-scale flow pore medium full-band elastic wave dispersion attenuation analysis method | |
| Markov et al. | Elastic properties of rocks containing oriented systems of ellipsoidal inclusions | |
| Shi et al. | Seismic wave dispersion and attenuation resulting from multiscale wave-induced fluid flow in partially saturated porous media |
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 |