CN104570084B - Across yardstick earthquake rock physicses attenuation model and the method for prediction decay and frequency dispersion - Google Patents
Across yardstick earthquake rock physicses attenuation model and the method for prediction decay and frequency dispersion Download PDFInfo
- Publication number
- CN104570084B CN104570084B CN201310520608.2A CN201310520608A CN104570084B CN 104570084 B CN104570084 B CN 104570084B CN 201310520608 A CN201310520608 A CN 201310520608A CN 104570084 B CN104570084 B CN 104570084B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- msup
- mfrac
- model
- 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.)
- Active
Links
Abstract
The invention discloses across yardstick earthquake rock physicses attenuation model and the method for prediction decay and frequency dispersion, to simulate the model of underground oil gas two-phase media characteristic, its construction method is the model：Under macro-scale Biot two-phase media elastic wave prorogation theory frameworks, it is introduced into the influence of " injection stream " in micro-scale crack in " BISQ " elastic wave prorogation theory, on this basis, the cycle stratification White model introduced under meso-scale, water-bearing layer is set to be stacked with gas-bearing bed interaction cycle, every layer of laterally unlimited extension, and be isotropic medium, every layer of thickness will be much smaller than earthquake wavelength and be more than particle size；Every layer of governing equation is " BISQ " Time Migration of Elastic Wave Equation, builds the earthquake rock physicses attenuation model across yardstick." injection stream " three kinds of yardstick attenuating mechanisms of " Biot " stream of the model including macro-scale, " the local stream " of meso-scale and microcosmic yardstick.The invention also discloses the method for predicting the model attenuation of P-wave and frequency dispersion, the pad value of the pad value predicted and experiment measurement coincide very much.
Description
Technical field
The present invention relates to exploration physical earth field, and in particular to across yardstick earthquake rock physicses attenuation model and prediction
Decay and the method for frequency dispersion.
Background technology
Underground oil and gas oil-containing gas medium shows obvious two-phase media feature in elasticity, studies earthquake in two-phase media
The propagation law of ripple is significant for the precision for improving oil-gas exploration.It can draw during propagation of the seismic wave in two-phase media
Flow of fluid (WIFF) is played, and then produces the attenuation of seismic wave and frequency dispersion (M ü ller, 2010).When seismic wave passes through, due to rock
Stone skeleton or pore-fluid skewness, barometric gradient is produced, so as to cause flow of fluid to cause decay and frequency dispersion.With fluid
Earthquake rock physicses theory about seismic wave attenuation by absorption can be divided into three classes from yardstick：Macro-scale attenuating mechanism, it is situated between
See yardstick and (be less than earthquake wavelength and be much larger than particle size) attenuating mechanism and micro-scale attenuating mechanism.Due to geological data master
For frequency at tens hertz, frequency bandwidth is limited, and the sight attenuating mechanism that is situated between is most important.The representative of macro-scale attenuation theory is Biot
(1956) elastic wave prorogation theory that the seismic wave derived is propagated in two-phase media, it is the classical theory frame for studying two-phase media
Frame, but the decay of the theoretical prediction is smaller compared with the decay actually measured, the permeability and the coefficient of viscosity tried to achieve are to decay
It is not inconsistent with the affecting laws of frequency dispersion with experiment.Hole of different shapes can produce different deformation, cause when seismic wave passes through
The flowing of fluid in hole, so as to produce decay, this mechanism is referred to as injection stream by Mavko and Nur (1975), micro-scale
Mainly as caused by " injection stream " effect in crack, its representational theory is that Dvorkin (1993) is set up for decay
" BISQ " elastic wave prorogation theory, the BISQ that Yang Dinghui (2000) have studied pore space anisotropy is theoretical, has expanded BISQ theories
Application, what firewood step on etc. (2000,2001,2003) by Biot is theoretical and injection stream theory based on, derived oil-containing water two
The Seismic Wave Propagation Equations of phase fluid pore media, and carry out forward simulation.Nie Jianxin (2004) is based on EFFECTIVE MEDIUM THEORY, carries
BISQ models in fractional saturation poroelasticity medium are gone out, and the propagation law of ripple have been analyzed.Tang Xiaoming (2011) considers
The interaction of hole and crack, it is theoretical to Biot and BISQ theories are promoted, and crack is analyzed to decay and frequency dispersion
The influence of feature.But its unpredictable seismic wave is in the decay of seismic band.The attenuating mechanism of meso-scale mainly has：
The cycle stratification patch saturation and spherical patchy saturation that White etc. (1975) is established, subsequent Dutta and Ode (1979) are used
More rigorous pore media mechanics is modified to the spherical patchy saturation that White is proposed.Vogelaar (2007) and Liu
Bright (2009,2010) solve the attenuation of P-wave and frequency dispersion for analyzing two kinds of White patchy saturations with the method for direct decoupling.
Dual porosity model that Pride and Berryman (2003) and Ba Jing (2011) are derived, contained using pseudo- spectrometry numerical simulation
Biot dissipates sees the seismic wave field flowed in porous media with middle, analyzes the attenuation of seismic wave caused by middle sight flow of fluid.Although
These models can portray decay of the seismic wave in seismic band, but not account for micro-scale crack to seismic wave in earthquake
The influence of frequency band.
The content of the invention
In view of the shortcomings of the prior art, it is an object of the invention to：Across yardstick earthquake rock physicses attenuation model, institute are provided
Established model more approaches the underground medium of complexity, can portray influence of the crack to seismic wave and seismic wave can be predicted in earthquake
The decay of frequency band.It is a further object of the present invention to provide utilize decay and frequency dispersion of the above-mentioned model prediction seismic wave in seismic band
Method.
The technical solution adopted by the present invention is：Across yardstick earthquake rock physicses attenuation model, the model contain for simulation underground
The model of oil gas two-phase media characteristic, it is characterised in that：
Building the method for the model is：
Under macro-scale Biot two-phase media elastic wave prorogation theory frameworks, it is introduced into microcosmic in " BISQ " elastic wave prorogation theory
The influence of " injection stream " in yardstick crack, on this basis, the cycle stratification White model introduced under meso-scale, make aqueous
Layer is stacked with gas-bearing bed interaction cycle, every layer of laterally unlimited extension, and is isotropic medium, and every layer of thickness will be much smaller than
Earthquake wavelength and it is more than particle size；Every layer of governing equation is " BISQ " Time Migration of Elastic Wave Equation, builds the earthquake across yardstick
Rock physicses attenuation model.
Further, the influence for being introduced into " injection stream " in micro-scale crack in " BISQ " elastic wave prorogation theory refers to：By spy
Spray penetration this parameter introducing " BISQ " Time Migration of Elastic Wave Equation is levied, portrays micro-scale crack " injection stream " effect size.
Further, comprising the concrete steps that for the model is built：
Step 1：By elastic constitutive relation and nonlinear New-tonian law, Biot models are established；
Step 2：The injection stream conservation of mass, the quality of pipeline stream of micro-crack are taken into full account on the basis of Biot models
Law of conservation, it is established that BISQ models；
Step 3:The meso-scale boundary condition of feature unit is taken into full account on the basis of BISQ models, is built across chi
Spend earthquake petrophysical model and its every layer of " BISQ " Time Migration of Elastic Wave Equation.
Further, in step 3, the feature unit refers to,
From the model each levels interface, take upper and lower two layers height it is half and as cylinder
Height, cylinder diameter take infinitely small, and such cylinder is referred to as feature unit；
The feature unit not only includes the macroscopic view " Biot streams " of longitudinal direction, the microcosmic " injection due to Effect of Fissure also be present
Stream ", and the different caused sights " local stream " that are situated between of Gas-Water Contant Fluid pressure；Numerous feature unit composition gas-bearing bed and water-bearing layer
It must be stacked in the cycle, every layer of laterally unlimited extension；
The meso-scale refers to the difference yardstick that can react the different caused sights " local stream " that are situated between of Gas-Water Contant Fluid pressure.
Further, described every layer " BISQ " Time Migration of Elastic Wave Equation be：
ρ11utt+ρ12Utt+b(ut-Ut)=P1uxx+Q1Uxx
ρ12utt+ρ22Utt-b(ut-Ut)=Q1uxx+R1Uxx(1)；
In formula, ρ11+ρ12=(1- φ) ρ1, ρ12+ρ22=φρf, Q1=(α-φ) F1, R1=φ F1,
Wherein, u, U are respectively solid and displacement of fluid, and subscript t, x is represented to its seeking time and directional derivative；Subscript s, f,
D represents solid, fluid and rock matrix respectively；φ is porosity, and ρ is density,For the structure factor of pore media, for ball
Shape solid particle r=0.5, η, κ are respectively the coefficient of viscosity and permeability of fluid；W is circular frequency, wBFor Biot characteristic frequencies；K
For bulk modulus, α is Biot coefficients, α=1-Kd/K1, J1、J0For single order and zero Bessel function, RBIt is characterized jet flow
Length, μ be skeleton modulus of shearing, ρ11、ρ12、ρ22For the quality coefficient of two-phase media.
Further, the model parameter of across the yardstick earthquake rock physicses attenuation model is explained by existing geologic information
Arrive.
The method for predicting the attenuation of P-wave and frequency dispersion of above-mentioned model：
This method includes：With the solid position of the upper and lower interface of feature unit in the method solving model of Dutta bilingual coupling
Move, and then try to achieve the plane wave modulus of feature unit, finally obtain decay and the frequency of across yardstick earthquake rock physicses attenuation model
Dissipate.
Further, in the solving model upper and lower interface of feature unit solid displacement, and then try to achieve feature unit
The specific method of plane wave modulus is：
Each feature unit in the model, when the external world applies a simple harmonic quantity powerWhen, feature unit can produce single shaft should
BecomeWherein (i)2=-1, it is unit imaginary number, it can thus be concluded that the equivalent plane ripple modulus of the model：
Uniaxial strain in formula is：
Wherein, ua、ubIt is characterized the solid displacement of the upper and lower interface of unit；
The frequency dispersion speed of feature unit can obtain：
Wherein, ρeFor equivalent density,
Subscript a, b represents the earthquake petrophysical parameter of levels respectively；
The quality factor of feature unit：
Wherein, real () represents the real part of access value, and imag () represents the imaginary part of access value；
Using Dutta to the method for the bilingual coupling of two-phase media Time Migration of Elastic Wave Equation of spherical White models to shown in formula 1
Time Migration of Elastic Wave Equation solved, Time Migration of Elastic Wave Equation is written as form first：
In formula, W=φ (U-u), ρb=[(1- φ) ρ1+φρf], H=P1+2Q1+R1, m=ρ22/φ
Subscript t represents time-derivative, ρ22For the quality coefficient of two-phase media,Q1=(α-φ)
F1, R1=φ F1；
The form for making the quasi-static solution of the equation shown in formula 6 is：
U=u (x) eiwtW=W (x) eiwt(7)；
Then solid and fluid structurecoupling displacement are decoupled, is written as form：
W=We+WdU=ue+ud(8)；
Wherein：ud=σdWd, ue=σeWe, the expression of subscript c, d is fast, Slow P-wave, brings formula (6) institute into together with formula (11)
Two Time Migration of Elastic Wave Equation as follows can be obtained in the equation shown：
Wherein, the displacement coefficient σ of Concerning With Fast-slow Wavesc、σdMeet following quadratic equation with one unknown：
Speed wave number kc、kdFor：
In formula, σ11=ρbσc+ρf, σ12=ρbσd+ρf, σ21=ρfσc+ m, σ11=ρfσd+ m, c1=H σc+(Qi+Ri)/φ, c2
=H σd+(Qi+Ri)/φ, c3=(Q1+R1)σe/φ+R1/φ2, c4=(Q1+R1)σd/φ+R1/φ2；
Fluid structurecoupling displacement can be obtained by being solved to formula 8 and solid displacement is：
W=Blcos (kcx)+B2sin(kcx)+B3cos((kdx)+B4sin(kdx)
U=σcB1cos(kcx)+σcB2sin(kcx)+σdB3cos(kdx)+σdB4sin(kdx) (12)；
Wherein, B1, B2, B3, B4 are undetermined coefficients；
According to the total stress of feature unit upper and lower interface and stream maintenance phase shift etc., at bed boundary, total stress is continuous, hole
Pressure is continuous, solid displacement is continuous and fluid flow is equal, it is assumed that the dielectric thickness comprising two kinds of fluids is da、db, establish as follows
Boundary condition：
1)τa=-PeeiwtX=-da
2)ua=UaX=-da
3)τb=-PeeiwtX=db
4)ub=UbX=db
5)τa=τbX=0
6)
7)ua=ubX=0
8)Wa=WbX=0 (13)
Wherein, subscript a, b represents saturation different fluid medium, pfFor Fluid pressure；
By eight absorbing boundary equations more than solving, eight undetermined coefficients of levels solid displacement can be tried to achieve, and then can obtain
Solid displacement, the attenuation of P-wave and frequency dispersion of across yardstick earthquake petrophysical model are calculated finally by formula (2)-formula (5).
In order that institute's established model more approaches the underground medium of complexity, influence of the crack to seismic wave can be portrayed and can be with
Decay of the seismic wave in seismic band is predicted, model of the present invention provides across yardstick earthquake rock physicses attenuation model,
Both the influence in micro-scale crack had been introduced, and can portrays decay and frequency dispersion of the seismic wave in seismic band.This across yardstick earthquake
The beneficial effect of rock physicses attenuation model be considered macro-scale " Biot streams ", meso-scale " local stream " with it is micro-
Seeing the attenuating mechanism of yardstick " injection stream " three yardsticks influences, due to the saturation fluid of underground medium when predicting seimic wave propagation
Decay and frequency dispersion of the seismic wave caused by flowing in seismic band.Earthquake petrophysical model not only allows for harder circular hole gap, also
Consider soft hole-crack so that model is more nearly underground Complicated Geologic Condition.The pad value predicted and laboratory measurement
Pad value more coincide.
Brief description of the drawings
Fig. 1 is the schematic diagram of across yardstick earthquake rock physicses attenuation model.
Fig. 2 is the technical scheme flow chart for building across yardstick earthquake rock physicses attenuation model.
Fig. 3 a are using White1975 formula, Biot equations as cycle stratification White model cootrols equation and across frequency
Attenuation of P-wave figure with earthquake petrophysical model.
Fig. 3 b are White1975 formula, Biot equations as cycle stratification White model cootrols equation and across frequency band
The compressional wave frequency dispersion figure of earthquake petrophysical model.
Fig. 4 a are the attenuation of P-wave figures of this attenuation model under the conditions of different permeabilities.
Fig. 4 b are the compressional wave frequency dispersion figures of this attenuation model under the conditions of different permeabilities.
The attenuation of P-wave figure of this attenuation model under the conditions of Fig. 5 a difference gas saturation.
The compressional wave frequency dispersion figure of this attenuation model under the conditions of Fig. 5 b difference gas saturation.
The attenuation of P-wave figure of this attenuation model under the conditions of Fig. 5 c difference gas saturation.
The compressional wave frequency dispersion figure of this attenuation model under the conditions of Fig. 5 d difference gas saturation.
Fig. 6 a are the attenuation of P-wave figures of this attenuation model under the conditions of different FRACTURE CHARACTERISTICS spray penetrations.
Fig. 6 b are the compressional wave frequency dispersion figures of this attenuation model under the conditions of different FRACTURE CHARACTERISTICS spray penetrations.
Embodiment
The present invention will be further described with reference to the accompanying drawings and examples.
Embodiment 1.Across yardstick earthquake rock physicses attenuation model, it is characterised in that：The model is simulation underground oil gas
The model of two-phase media characteristic, macro-scale, meso-scale, the attenuating mechanism of micro-scale are considered, in given thing
Under the conditions of property and fluid parameter, decay and frequency dispersion of the seismic wave in seismic band can be predicted.
Building the method for the model is：
Under macro-scale Biot two-phase media elastic wave prorogation theory frameworks, it is introduced into microcosmic in " BISQ " elastic wave prorogation theory
The influence of " injection stream " in yardstick crack, on this basis, the cycle stratification White model introduced under meso-scale, make aqueous
Layer is stacked with gas-bearing bed interaction cycle, every layer of laterally unlimited extension, and is isotropic medium, and every layer of thickness will be much smaller than
Earthquake wavelength and it is more than particle size；Every layer of governing equation is " BISQ " Time Migration of Elastic Wave Equation, builds the earthquake across yardstick
Rock physicses attenuation model.
The specific steps of the model are built as shown in Fig. 2 technical scheme flow chart, including：
Step 1：By elastic constitutive relation and nonlinear New-tonian law, Biot models are established；Biot models consider macroscopic view
The effect of " Biot streams ".
Step 2：The injection stream conservation of mass, the quality of pipeline stream of micro-crack are taken into full account on the basis of Biot models
Law of conservation, it is established that BISQ models.It is introduced into the influence of " injection stream " in micro-scale crack in " BISQ " elastic wave prorogation theory
Refer to：By feature spray penetration, this parameter introduces " BISQ " Time Migration of Elastic Wave Equation, portrays micro-scale crack " injection stream " work
Use size.BISQ models consider " injection stream " effect.
Step 3:The meso-scale boundary condition of feature unit is taken into full account on the basis of BISQ models, is built across chi
Spend the math equation of earthquake petrophysical model.Meso-scale just refers to the yardstick between both macro and micro；It is generally acknowledged that
Its yardstick is between nanometer and millimeter.Meso-scale of the present invention refers to react caused by Gas-Water Contant Fluid pressure difference
It is situated between and sees the difference yardstick of " local stream ".
The feature unit refers to, from the model each levels interface, takes the half of upper and lower two layers of height
And the height as cylinder, cylinder diameter take infinitely small, such cylinder is referred to as feature unit；The feature unit is not
The macroscopic view " Biot streams " of longitudinal direction is only included, is also existed due to microcosmic " injection stream " of Effect of Fissure, and Gas-Water Contant Fluid pressure
It is situated between caused by different and sees " local stream ".Fig. 1 illustrates the schematic diagram of across yardstick earthquake rock physicses attenuation model.Gassiness in figure
Layer and water-bearing layer cycle are stacked, and laterally unlimited extension, is made up of numerous cylinder feature unit.
Described every layer " BISQ " Time Migration of Elastic Wave Equation be：
ρ11μtt+ρ12Utt+b(ut-Ut)=P1uxx+Q1Uxx
ρ12utt+ρ22Utt-b(ut-Ut)=Q1uxx+R1Uxx(1)；
In formula, ρ11+ρ12=(1- φ) ρ1, ρ12+ρ22=φ ρf, Q1=(α-φ) F1, R1=φ F1,
Wherein, u, U are respectively solid and displacement of fluid, and subscript t, x is represented to its seeking time and directional derivative；Each coefficient table
Up in formula,
Subscript s, f, d represent solid, fluid and rock matrix respectively；φ is porosity, and ρ is density,For pore media
Structure factor, for spherical solid particles r=0.5, η, κ are respectively the coefficient of viscosity and permeability of fluid；W is circular frequency, wB
For Biot characteristic frequencies；K is bulk modulus, and α is Biot coefficients, α=1-Kd/K1, J1、J0For single order and zero Bessel function,
RBIt is characterized jet flow length, μ is the modulus of shearing of skeleton, ρ11、ρ12、ρ22For the quality coefficient of two-phase media.
The model parameter of across the yardstick earthquake rock physicses attenuation model can be explained to obtain by existing geologic information.
This across yardstick earthquake rock physicses attenuation model.It includes " Biot " stream, " local of meso-scale of macro-scale
" injection stream " three kinds of yardstick attenuating mechanisms of stream " and microcosmic yardstick, physical parameter (permeability, feature injection can be predicted
Length), fluid parameter (gas saturation) is decayed under seismic band by seismic wave and the influence of frequency dispersion.The gas saturation
SgThe ratio of feature unit gross thickness is accounted for for the thickness of gas-bearing bed.Specifically formula is：Sg=da/(da+db), wherein, daIt is characterized
The distance of unit top interface distance interface, dbWith the being characterized unit distance (referring to Fig. 1) of cross-sectional distance interface.
Embodiment 2.The method for predicting the attenuation of P-wave and frequency dispersion of above-mentioned across yardstick earthquake rock physicses attenuation model, bag
Include：With the solid displacement of the upper and lower interface of feature unit in the method solving model of Dutta bilingual coupling, and then try to achieve feature list
The plane wave modulus of member, finally obtain decay and the frequency dispersion of across yardstick earthquake rock physicses attenuation model.
The solid displacement of the upper and lower interface of feature unit in the solving model, and then try to achieve the plane wave mould of feature unit
The specific method of amount is：
Each feature unit in the model, when the external world applies a simple harmonic quantity powerWhen (simulation longitudinal wave propagation), feature list
Member can produce uniaxial strain, wherein (i)2=-1, it is unit imaginary number, it can thus be concluded that the equivalent plane ripple modulus of the model：
Uniaxial strain in formula is：
Wherein, ua、ubIt is characterized the solid displacement of the upper and lower interface of unit.
The frequency dispersion speed of feature unit can obtain：
Wherein, ρeFor equivalent density,
Subscript a, b represents the earthquake petrophysical parameter of levels respectively.
The quality factor of feature unit：
Wherein, real () represents the real part of access value, and imag () represents the imaginary part of access value.
Using Dutta to the method for the bilingual coupling of two-phase media Time Migration of Elastic Wave Equation of spherical White models to formula (1) institute
Show that Time Migration of Elastic Wave Equation is solved, Time Migration of Elastic Wave Equation is written as form first：
In formula, W=φ (U-u), ρb=[(1- φ) ρ1+φρf], H=P1+2Q1+R1, m=ρ22/φ
Subscript t represents time-derivative, ρ22For the quality coefficient of two-phase media, P1、Q1、R1Implication and formula (1) in one
Cause.
The form for making the quasi-static solution of equation shown in formula (6) is：
U=u (x) eiwtW=W (x) eiwt (7)。
Then solid and fluid structurecoupling displacement are decoupled, is written as form：
W=We+WdU=ue+ud(8)；
Wherein：ud=σdWd, ue=σeWe, the expression of subscript c, d is fast, Slow P-wave, is brought into together with formula (11) in equation (6)
Two Time Migration of Elastic Wave Equation as follows can be obtained：
Wherein, the displacement coefficient σ of Concerning With Fast-slow Wavesc、σdMeet following quadratic equation with one unknown：
Speed wave number ke、kdFor：
In formula, σ11=ρbσc+ρf, σ12=ρbσd+ρf, σ21=ρfσc+ m, σ22=ρfσd+ m, c1=H σc+(Qi+Ri)/φ, c2
=H σd+(Qi+Ri)/φ, c3=(Qt+Rt)σe/φ+R1/φ2, c4=(Qt+Rt)σd/φ+R1/φ2。
Fluid structurecoupling displacement can be obtained and solid displacement is by being solved to formula (8)：
W=Blcos (kcx)+B2sin(kcx)+B3cos((kdx)+B4sin(kdx)
U=σcBlcos(kcx)+σcB2sin(kcx)+σdB3cos(kdx)+σdB4sin(kdx) 12)；
Wherein, B1, B2, B3, B4 are undetermined coefficients.
According to the total stress of feature unit upper and lower interface and stream maintenance phase shift etc., at bed boundary, total stress is continuous, hole
Pressure is continuous, solid displacement is continuous and fluid flow is equal, it is assumed that the dielectric thickness comprising two kinds of fluids is da、db, establish as follows
Boundary condition "：
1)τa=-PeeiwtX=-da
2)ua=UaX=-da
3)τb=-PeeiwtX=db
4)ub=UbX=db
5)τa=τbX=0
6)
7)ua=ubX=0
8)Wa=WbX=0 (13)
Wherein, subscript a, b represents saturation different fluid medium, pfFor Fluid pressure；
By eight absorbing boundary equations more than solving, eight undetermined coefficients of levels solid displacement can be tried to achieve, and then can obtain
Solid displacement, the attenuation of P-wave and frequency dispersion of across yardstick earthquake petrophysical model are calculated finally by formula 2-5.
With it, physical parameter (such as permeability, feature spray penetration), fluid parameter (such as gassiness can also be predicted
Saturation degree) seismic wave is decayed under seismic band and the influence of frequency dispersion.
Across yardstick the earthquake rock physicses attenuation model and its prediction effect that the present invention designs, built compared to using conventional method
The pad value of mould, the pad value predicted and experiment measurement more coincide.
Fig. 3 a are using White1975 formula, Biot equations as cycle stratification White model cootrols equation and across frequency
Attenuation of P-wave figure with earthquake petrophysical model.Fig. 3 b are White1975 formula, Biot equations as cycle stratification
White model cootrols equation and the compressional wave frequency dispersion figure for shaking petrophysical model across band.It can be seen that add the " injection in crack
After stream " influences, seismic wave can occur to change accordingly in the decay of seismic band and frequency dispersion.
Fig. 4 a are the attenuation of P-wave figures of this attenuation model under the conditions of different permeabilities.Fig. 4 b are sheets under the conditions of different permeabilities
The compressional wave frequency dispersion figure of attenuation model.It can be seen that with the rise of permeability, decay dominant frequency is to high-frequency mobile, and damping peak is constant.
Fig. 5 a- Fig. 5 d are the attenuation of P-wave and frequency dispersion of this attenuation model under the conditions of different gas saturation.It can be seen that with containing
The increase of gas saturation, the decay of compressional wave and frequency dispersion first increase and then reduced.In gas saturation 0.1, this attenuation model
Attenuation of P-wave is maximum, consistent with Gautam (2003) experimental result.
Fig. 6 a are the attenuation of P-wave figures of this attenuation model under the conditions of different FRACTURE CHARACTERISTICS spray penetrations.Fig. 6 b are different cracks
The compressional wave frequency dispersion figure of this attenuation model under the conditions of feature spray penetration.Curve 1 represents to be used as the attenuation model by the use of " Biot " equation
The vertical frequency decay (a) and frequency dispersion (b) that every layer of governing equation calculates, curve 2,3,4,5,6,7,8,9 represent that injection characteristics are grown respectively
Spend for 1,0.4,0.01,0.05,0.02,0.01,0.001,0.0001m.It can be seen that when feature spray penetration gradually increases, this
The attenuation of P-wave and frequency dispersion of attenuation model gradually increase, and the effect of " injection stream " in crack reduces, when feature spray penetration reduces
When to a certain extent, the attenuation of P-wave and frequency dispersion of this attenuation model are constant, the influence of " injection stream " in crack to this attenuation model
Keep constant.
Above embodiment is to illustrate the invention and not to limit the present invention.
Claims (8)
1. a kind of method for building across yardstick earthquake rock physicses attenuation model, the model is simulation underground oil gas two-phase media
The model of characteristic, it is characterised in that：
Building the method for the model is：
Under macro-scale Biot two-phase media elastic wave prorogation theory frameworks, micro-scale in " BISQ " elastic wave prorogation theory is introduced into
The influence of " injection stream " in crack, on this basis, introduce meso-scale under cycle stratification White model, make water-bearing layer with
Gas-bearing bed interaction cycle is stacked, every layer of laterally unlimited extension, and is isotropic medium, and every layer of thickness will be much smaller than earthquake
Wavelength and it is more than particle size；Every layer of governing equation is " BISQ " Time Migration of Elastic Wave Equation, builds the earthquake rock across yardstick
Physical Attenuation model.
A kind of 2. method for building across yardstick earthquake rock physicses attenuation model as claimed in claim 1, it is characterised in that：
The influence for being introduced into " injection stream " in micro-scale crack in " BISQ " elastic wave prorogation theory refers to：By feature spray penetration this
One parameter introduces " BISQ " Time Migration of Elastic Wave Equation, portrays micro-scale crack " injection stream " effect size.
A kind of 3. method for building across yardstick earthquake rock physicses attenuation model as claimed in claim 1, it is characterised in that：Structure
Build comprising the concrete steps that for the model：
Step 1：By elastic constitutive relation and nonlinear New-tonian law, Biot models are established；
Step 2：The injection stream conservation of mass, the conservation of mass of pipeline stream of micro-crack are taken into full account on the basis of Biot models
Law, it is established that BISQ models；
Step 3:The meso-scale boundary condition of feature unit is taken into full account on the basis of BISQ models, is built across yardstick
Shake petrophysical model and its every layer of " BISQ " Time Migration of Elastic Wave Equation.
A kind of 4. method for building across yardstick earthquake rock physicses attenuation model as claimed in claim 3, it is characterised in that：Step
In rapid 3, the feature unit refers to,
From the model each levels interface, half and as cylinder the height of upper and lower two layers of height is taken, circle
Column diameter takes infinitely small, and such cylinder is referred to as feature unit；
The feature unit not only includes the macroscopic view " Biot streams " of longitudinal direction, also exists due to microcosmic " injection stream " of Effect of Fissure,
And the different caused sights " local stream " that are situated between of Gas-Water Contant Fluid pressure；Numerous feature unit composition gas-bearing bed and water-bearing layer obtain all
Phase is stacked, every layer of laterally unlimited extension；
The meso-scale refers to the difference yardstick that can react the different caused sights " local stream " that are situated between of Gas-Water Contant Fluid pressure.
5. across the yardstick earthquake rock physicses attenuation model of a kind of structure as described in any one claim of Claim 1-3
Method, it is characterised in that：Described every layer " BISQ " Time Migration of Elastic Wave Equation be：
ρ11utt+ρ12Utt+b(ut-Ut)=Pluxx+QlUxx
ρ12utt+ρ22Utt-b(ut-Ut)=Q1uxx+RlUxx(1)；
In formula,ρ12+ρ22=φ ρf, Q1=(α-φ) F+, R1=φ F1,
Wherein, u, U are respectively solid and displacement of fluid, and subscript t, x is represented to its seeking time and directional derivative；S, f, d points of subscript
Solid, fluid and rock matrix are not represented；φ is porosity, and ρ is density,For the structure factor of pore media, for spherical solid
Body particle r=0.5, η, κ are respectively the coefficient of viscosity and permeability of fluid；W is circular frequency, wBFor Biot characteristic frequencies；K is body
Product module amount, α are Biot coefficients,J1、J0For single order and zero Bessel function, RBIt is characterized jet flow length
Degree, μ be skeleton modulus of shearing, ρ11、ρ12、ρ22For the quality coefficient of two-phase media.
A kind of 6. method for building across yardstick earthquake rock physicses attenuation model as claimed in claim 5, it is characterised in that：Institute
The model parameter for stating across yardstick earthquake rock physicses attenuation model is explained to obtain by existing geologic information.
7. predict a kind of across yardstick earthquake rock physicses attenuation model of structure described in any one claim of claim 1 to 4
Attenuation of P-wave and frequency dispersion method, it is characterised in that：
This method includes：With the solid displacement of the upper and lower interface of feature unit in the method solving model of Dutta bilingual coupling, enter
And the plane wave modulus of feature unit is tried to achieve, finally obtain decay and the frequency dispersion of across yardstick earthquake rock physicses attenuation model.
8. the method for prediction attenuation of P-wave and frequency dispersion as claimed in claim 7, it is characterised in that：It is special in the solving model
The solid displacement of the upper and lower interface of unit is levied, and then the specific method for trying to achieve the plane wave modulus of feature unit is：
Each feature unit in the model, when the external world applies a simple harmonic quantity power PeeIwtWhen, feature unit can produce uniaxial strain θ
eIwt, wherein (i)2=-1, it is unit imaginary number, it can thus be concluded that the equivalent plane ripple modulus of the model：
<mrow>
<msup>
<mi>p</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>p</mi>
<mi>e</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>I</mi>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msup>
</mrow>
<mrow>
<msup>
<mi>&theta;e</mi>
<mrow>
<mi>I</mi>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Uniaxial strain in formula is：
<mrow>
<msup>
<mi>&theta;e</mi>
<mrow>
<mi>I</mi>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>u</mi>
<mi>b</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>I</mi>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>a</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>I</mi>
<mi>w</mi>
<mi>t</mi>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>d</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>b</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, ua、ubIt is characterized the solid displacement of the upper and lower interface of unit；
The frequency dispersion speed of feature unit can obtain：
<mrow>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>/</mo>
<msqrt>
<mrow>
<msup>
<mi>P</mi>
<mo>*</mo>
</msup>
<mo>/</mo>
<msub>
<mi>&rho;</mi>
<mi>e</mi>
</msub>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
Wherein, ρeFor equivalent density,
<mrow>
<msub>
<mi>&rho;</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<msub>
<mi>d</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&phi;</mi>
<mi>a</mi>
</msub>
<msubsup>
<mi>&rho;</mi>
<mi>f</mi>
<mi>a</mi>
</msubsup>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mi>a</mi>
</msub>
</mrow>
<mo>)</mo>
<msubsup>
<mi>&rho;</mi>
<mi>s</mi>
<mi>a</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>b</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&phi;</mi>
<mi>b</mi>
</msub>
<msubsup>
<mi>&rho;</mi>
<mi>f</mi>
<mi>b</mi>
</msubsup>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mi>b</mi>
</msub>
</mrow>
<mo>)</mo>
<msubsup>
<mi>&rho;</mi>
<mi>s</mi>
<mi>b</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msub>
<mi>d</mi>
<mi>b</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Subscript a, b represents the earthquake petrophysical parameter of levels respectively；
The quality factor of feature unit：
<mrow>
<mi>Q</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>a</mi>
<mi>l</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>P</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>P</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, real () represents the real part of access value, and imag () represents the imaginary part of access value；
Using Dutta to the method for the bilingual coupling of two-phase media Time Migration of Elastic Wave Equation of spherical White models to the bullet shown in formula 1
Property wave equation is solved, and Time Migration of Elastic Wave Equation is written as into form first：
<mrow>
<msub>
<mi>&rho;</mi>
<mi>b</mi>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>t</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<msub>
<mi>W</mi>
<mrow>
<mi>t</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>Hu</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>&phi;</mi>
</mfrac>
<msub>
<mi>W</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>t</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>mW</mi>
<mrow>
<mi>t</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>&phi;</mi>
</mfrac>
<msub>
<mi>u</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<msup>
<mi>&phi;</mi>
<mn>2</mn>
</msup>
</mfrac>
<msub>
<mi>W</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<mi>b</mi>
<msup>
<mi>&phi;</mi>
<mn>2</mn>
</msup>
</mfrac>
<msub>
<mi>W</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, W=φ (U-u),H=P1+2Q1+R1, m=p22/φ
Subscript t represents time-derivative, ρ22For the quality coefficient of two-phase media,Q1=(α-φ) F1, R1
=φ F1；
The form for making the quasi-static solution of the equation shown in formula 6 is：
U=u (x) eiwtW=W (x) eIwt(7)；
Then solid and fluid structurecoupling displacement are decoupled, is written as form：
W=We+WdU=ue+ud(8)；
Wherein：ud=σdWd, ue=σeWe, the expression of subscript c, d is fast, Slow P-wave, is brought into together with formula (11) shown in formula (6)
Two Time Migration of Elastic Wave Equation as follows can be obtained in equation：
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>k</mi>
<mi>c</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<msub>
<mi>W</mi>
<mi>x</mi>
</msub>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<msubsup>
<mi>k</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
<msub>
<mi>W</mi>
<mi>d</mi>
</msub>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
<mo>;</mo>
</mrow>
Wherein, the displacement coefficient σ of Concerning With Fast-slow Wavesc、σdMeet following quadratic equation with one unknown：
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mrow>
<mi>&phi;</mi>
</mfrac>
<msub>
<mi>&rho;</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<msub>
<mi>H&rho;</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>b</mi>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mrow>
<msup>
<mi>&phi;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>-</mo>
<mi>m</mi>
<mi>H</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mi>i</mi>
<mi>b</mi>
<mi>H</mi>
</mrow>
<mrow>
<msup>
<mi>&phi;</mi>
<mn>2</mn>
</msup>
<mi>w</mi>
</mrow>
</mfrac>
<mo>)</mo>
<mi>&sigma;</mi>
<mo>+</mo>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mrow>
<msup>
<mi>&phi;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>&phi;</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>i</mi>
<mi>b</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>w&phi;</mi>
<mn>1</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
<mo>;</mo>
</mrow>
Speed wave number kc、kdFor：
<mrow>
<mfrac>
<msubsup>
<mi>k</mi>
<mi>e</mi>
<mn>3</mn>
</msubsup>
<msup>
<mi>w</mi>
<mn>3</mn>
</msup>
</mfrac>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>11</mn>
</msub>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>31</mn>
</msub>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>ibc</mi>
<mn>3</mn>
</msub>
<mo>/</mo>
<msup>
<mi>&phi;</mi>
<mn>3</mn>
</msup>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<msubsup>
<mi>k</mi>
<mi>d</mi>
<mn>3</mn>
</msubsup>
<msup>
<mi>w</mi>
<mn>3</mn>
</msup>
</mfrac>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>13</mn>
</msub>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>33</mn>
</msub>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>ibc</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msup>
<mi>&phi;</mi>
<mn>3</mn>
</msup>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
In formula, σ11=ρbσc+ρf, σ13=ρbσd+ρf, σ21=ρfσc+ m, σ33=ρfσd+ m, c1=H σc+(Qi+Ri)/φ, c2=H σd
+(Qi+Ri)/φ,
Fluid structurecoupling displacement can be obtained by being solved to formula 8 and solid displacement is：
W=B1cos (kcx)+B2sin(kcx)+B3cos(kdx)+B4sin(kdx)
U=σcB1cos(kcx)+σcB2sin(kcx)+σdB3cos(kdx)+σdB4sin(kdx) (12)；
Wherein, B1, B2, B3, B4 are undetermined coefficients；
According to the total stress of feature unit upper and lower interface and stream maintenance phase shift etc., at bed boundary, total stress is continuous, pore pressure
Continuously, solid displacement is continuous and fluid flow is equal, it is assumed that the dielectric thickness of two kinds of fluids of saturation is da、db, establish such as lower boundary
Condition：
1)τa=-PeeiwtX=-da
2)ua=UaX=-da
3)τb=-PeeiwtX=db
4)ub=UbX=db
5)τa=τbX=0
6)X=0
7)ua=ubX=0
8)Wa=WbX=0 (13)
Wherein, subscript a, b represents saturation different fluid medium, pfFor Fluid pressure；
By eight absorbing boundary equations more than solving, eight undetermined coefficients of levels solid displacement can be tried to achieve, and then solid can be obtained
Displacement, the attenuation of P-wave and frequency dispersion of across yardstick earthquake petrophysical model are calculated finally by formula (2)-formula (5).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310520608.2A CN104570084B (en) | 2013-10-29 | 2013-10-29 | Across yardstick earthquake rock physicses attenuation model and the method for prediction decay and frequency dispersion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310520608.2A CN104570084B (en) | 2013-10-29 | 2013-10-29 | Across yardstick earthquake rock physicses attenuation model and the method for prediction decay and frequency dispersion |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104570084A CN104570084A (en) | 2015-04-29 |
| CN104570084B true CN104570084B (en) | 2018-01-05 |
Family
ID=53086598
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310520608.2A Active CN104570084B (en) | 2013-10-29 | 2013-10-29 | Across yardstick earthquake rock physicses attenuation model and the method for prediction decay and frequency dispersion |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104570084B (en) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106959469A (en) * | 2017-04-14 | 2017-07-18 | 中国石油天然气股份有限公司 | The speed and decay analog analysing method and device of seismic wave |
| CN109298443B (en) * | 2017-07-25 | 2020-12-11 | 中国石油化工股份有限公司 | Method for evaluating elastic wave of double-fluid jet flow model and computer readable storage medium |
| CN107389794B (en) * | 2017-08-14 | 2020-01-10 | 中国石油大学(北京) | Method and system for measuring rock attenuation coefficient |
| CN110967744A (en) * | 2018-09-28 | 2020-04-07 | 中国石油化工股份有限公司 | Oil-gas identification method based on biphase medium theory and storage medium |
| CN110058310B (en) * | 2019-04-16 | 2021-07-09 | 中国石油大学(华东) | Seismic wave velocity frequency dispersion and attenuation prediction method and device for seismic rock |
| CN110261898A (en) * | 2019-04-26 | 2019-09-20 | 中国石油化工股份有限公司 | Well logging and seismic velocity matching process based on the analysis of earthquake petrophysics experiment |
| CN110456412B (en) * | 2019-07-09 | 2021-06-08 | 河海大学 | Method for identifying fluid saturation of carbonate reservoir based on post-stack seismic data |
| CN112230278A (en) * | 2019-07-15 | 2021-01-15 | 中国石油天然气集团有限公司 | Seepage field characteristic parameter determination method and device |
| CN111562613B (en) * | 2020-04-15 | 2023-03-31 | 东华理工大学 | Method for analyzing seismic wave frequency-dependent reflection coefficient of thin reservoir or mutual reservoir model |
| CN111812709B (en) * | 2020-07-08 | 2022-03-04 | 中国石油大学(北京) | Method, device and equipment for establishing multi-scale wave induced flow model |
| CN113093276B (en) * | 2021-03-19 | 2022-03-22 | 中国石油大学(北京) | Method, device, equipment and system for predicting seismic wave velocity dispersion and attenuation |
| CN114114403B (en) * | 2021-12-22 | 2023-06-27 | 东北石油大学 | Anisotropic attenuation medium simulation method based on fractional order Laplace operator |
| CN114333248B (en) * | 2022-01-11 | 2022-09-23 | 西南交通大学 | High-speed rail earthquake early warning method, device, equipment and readable storage medium |
| CN115774286B (en) * | 2022-12-06 | 2023-09-05 | 上海电子信息职业技术学院 | Longitudinal wave speed prediction method for pore medium containing viscoelastic fluid and application of longitudinal wave speed prediction method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102141545A (en) * | 2010-11-27 | 2011-08-03 | 江西理工大学 | Method for testing rock mass mechanics parameters based on explosion seismic wave space-time attenuation law |
| CN102478670A (en) * | 2010-11-24 | 2012-05-30 | 中国石油天然气集团公司 | Method for forecasting reservoir fluid property through earthquake attenuation attribute |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7675818B2 (en) * | 2005-07-28 | 2010-03-09 | ExxonMobil Upstream Resarch Co. | Method for tomographic inversion by matrix transformation |
-
2013
- 2013-10-29 CN CN201310520608.2A patent/CN104570084B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102478670A (en) * | 2010-11-24 | 2012-05-30 | 中国石油天然气集团公司 | Method for forecasting reservoir fluid property through earthquake attenuation attribute |
| CN102141545A (en) * | 2010-11-27 | 2011-08-03 | 江西理工大学 | Method for testing rock mass mechanics parameters based on explosion seismic wave space-time attenuation law |
Non-Patent Citations (1)
| Title |
|---|
| 含孔隙、裂隙介质弹性波动的统一理论--Biot理论的推广;唐晓明;《中国科学：地球科学》;20111231;第41卷(第6期);第784-795页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104570084A (en) | 2015-04-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104570084B (en) | Across yardstick earthquake rock physicses attenuation model and the method for prediction decay and frequency dispersion | |
| Bear | Dynamics of fluids in porous media | |
| Zimmerman | Compressibility of sandstones | |
| CN106837315B (en) | Method for representing coupling effect of fractured carbonate rock matrix and fractures | |
| Zhao et al. | Simulation of a multistage fractured horizontal well in a tight oil reservoir using an embedded discrete fracture model | |
| Ma et al. | Torsional dynamic response of a pipe pile in homogeneous unsaturated soils | |
| Nian et al. | Submarine slope failure due to overpressure fluid associated with gas hydrate dissociation | |
| Chen et al. | Interaction between internal waves and a permeable seabed | |
| CN104298797A (en) | Method for determining remaining oil of fracture and hole type oil reservoir under encirclement of high flow guide channel | |
| Huang et al. | DEM-LBM simulation of stress-dependent absolute and relative permeabilities in porous media | |
| Pan et al. | Influence of mineral compositions of rocks on mechanical properties | |
| Mei et al. | Fractal analysis of shape factor for matrix-fracture transfer function in fractured reservoirs | |
| Liu et al. | Numerical simulation on hydromechanical coupling in porous media adopting three-dimensional pore-scale model | |
| Peng et al. | Study of the borehole hydraulic fracturing and the principle of gas seepage in the coal seam | |
| Du | Anisotropic Rock Poroelasticity Evolution in Ultra‐low Permeability Sandstones under Pore Pressure, Confining Pressure, and Temperature: Experiments with Biot's Coefficient | |
| CN108875115B (en) | Method for determining rock strength | |
| Wu et al. | Effect of stress and material barriers on hydraulic fracture height containment in layered formations | |
| Gupta et al. | Effect of irregularity on the propagation of torsional surface waves in an initially stressed anisotropic poro-elastic layer | |
| Wang et al. | Hydrodynamic Analysis of CO2 Migration in Heterogeneous Rocks: Conventional and Micro‐Bubble CO2 Flooding Experiments and Pore‐Scale Numerical Simulations | |
| Qiu et al. | Numerical analysis of 1‐D compression wave propagation in saturated poroelastic media | |
| Si et al. | Numerical simulation on mesoscale mechanism of seepage in coal fractures by fluid-sloid coupling method | |
| Qian et al. | Numerical simulation of direct shear tests on gas hydrate-bearing sediments using the discrete element method | |
| CN113250659A (en) | Gas injection experiment simulation method and device for improving displacement efficiency | |
| CN112649861A (en) | Fractured reservoir rock physical modeling method and system | |
| McAllister et al. | Nonlinear site response analyses and high frequency dilation pulses |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |