CN106959469A - Method and device for simulating and analyzing velocity and attenuation of seismic wave - Google Patents
Method and device for simulating and analyzing velocity and attenuation of seismic wave Download PDFInfo
- Publication number
- CN106959469A CN106959469A CN201710242793.1A CN201710242793A CN106959469A CN 106959469 A CN106959469 A CN 106959469A CN 201710242793 A CN201710242793 A CN 201710242793A CN 106959469 A CN106959469 A CN 106959469A
- Authority
- CN
- China
- Prior art keywords
- medium
- seismic wave
- wave
- saturation
- frequency
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000011435 rock Substances 0.000 claims abstract description 101
- 239000011159 matrix material Substances 0.000 claims abstract description 65
- 239000007787 solid Substances 0.000 claims description 47
- 230000000644 propagated effect Effects 0.000 claims description 27
- 239000012530 fluid Substances 0.000 claims description 22
- 238000003780 insertion Methods 0.000 claims description 17
- 230000037431 insertion Effects 0.000 claims description 17
- 238000004088 simulation Methods 0.000 claims description 14
- 239000011148 porous material Substances 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 abstract description 4
- 229920006395 saturated elastomer Polymers 0.000 abstract 4
- 230000001902 propagating effect Effects 0.000 abstract 1
- 230000008859 change Effects 0.000 description 15
- 238000003860 storage Methods 0.000 description 14
- 238000004590 computer program Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 8
- 230000007246 mechanism Effects 0.000 description 8
- 238000012545 processing Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the application provides a method and a device for simulating and analyzing the velocity and attenuation of seismic waves, wherein the method comprises the following steps: establishing a saturated pore rock model and a saturated fracture rock model of a medium in a work area; embedding the saturated fracture rock model into the saturated pore rock model under a set condition to obtain an elastic coefficient matrix of the medium; the elasticity coefficient matrix comprises elasticity coefficients of the medium in all directions; determining the complex velocity of seismic waves propagating in the medium according to the elastic coefficient matrix of the medium; and determining the inverse quality factor of the medium and the propagation speed of the seismic waves in the medium according to the complex speed. The analysis result obtained by the embodiment of the application can more objectively reflect the influence of the underground medium on the seismic wave velocity and the energy attenuation.
Description
Technical field
The application is related to the speed of seismic wave and energy attenuation analysis technical field in geophysical exploration, more particularly, to
The speed and decay analog analysing method and device of a kind of seismic wave.
Background technology
In seismic prospecting, the analysis that seimic wave velocity decays in actual seismic data is most important, and is situated between based on equivalent
The velocity attenuation analytical technology of matter modeling is important method therein.When seismic wave is propagated in Fracture-Pore rock, its frequency
The solid relative motion of the stream being primarily due to when seismic wave passes through produced by reservoir particle (ripple causes stream) for dissipating and decaying, and one
The characteristic that the crack that group is closely aligned can cause reservoir to have equivalent transverse isotropy (TI) medium again.Therefore, split
The solid relative motion of stream and Fracture orientation arrangement in seam porosity rock cause the velocity attenuation and frequency and incidence angle of seismic wave
Degree is closely related to meet actual Fracture-Pore petrophysical model, it is necessary to set up, and based on the model, to the speed and frequency of medium
It is scattered to be analyzed.
It is considered as earthquake that ripple under meso-scale (observing result is much larger than pore-size, much smaller than earthquake wavelength), which causes stream,
Seismic-wave dispersion and the main cause of decay, can be divided into two classes at present under frequency band:The first kind is due to institute in crack and hole
Caused by the elasticity of fluid attribute difference of filling, Equations of The Second Kind is due to the difference of crack and hole in terms of solid elastic attribute
Cause.Current most of petrophysical models are all that both of these case is separated into discussion, but in actual reservoir, two types
Attenuating mechanism is simultaneous, therefore, is difficult to objectively respond underground medium to seimic wave velocity based on current currently existing scheme
With the influence of energy attenuation.
The content of the invention
The purpose of the embodiment of the present application is the speed and decay analog analysing method and device for providing a kind of seismic wave, with
Underground medium is objectively responded to seimic wave velocity and the influence of energy attenuation.
To reach above-mentioned purpose, the embodiment of the present application provides the speed and decay analog analysing method of a kind of seismic wave,
Comprise the following steps:
Set up the saturation porous rocks model and saturation crack rock model of medium in work area;
In the case where imposing a condition, by the saturation crack rock model insertion into the saturation porous rocks model, to obtain
Take the elastic coefficient matrix of the medium;The elastic coefficient matrix includes coefficient of elasticity of the medium in all directions;
The complex velocity that seismic wave is propagated in the medium is determined according to the elastic coefficient matrix of the medium;
The inverse quality factor and the seismic wave for determining the medium according to the complex velocity are propagated in the medium
Speed.
The speed and decay analog analysing method of the seismic wave of the embodiment of the present application, it is described to set up medium in target work area
Saturation porous rocks model and saturation crack rock model, including:
Determine the elastic modelling quantity and density of solid dielectric in work area;The solid dielectric includes the first solid dielectric and second
Solid dielectric;
The porosity of first solid dielectric is determined, after interstitial space is introduced into first solid dielectric, to
Filling first fluid in the interstitial space, forms saturation porous rocks model;The first fluid has the first bulk modulus
With the first density;
The crevice volume percentage of second solid dielectric is determined, particular volume is being introduced into second solid dielectric
After the fracture spaces of product module amount, second fluid is filled into the fracture spaces, saturation crack rock model is formed;Described second
Fluid has the second bulk modulus and the second density.
The speed and decay analog analysing method of the seismic wave of the embodiment of the present application, it is described in the case where imposing a condition, will be described
Saturation crack rock model insertion is into the saturation porous rocks model, to obtain the elastic coefficient matrix of the medium, bag
Include:
When the frequency of seismic wave in the medium is the low frequency less than lower frequency limit, by the saturation crack rock mould
Type is embedded in the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under low frequency;
When the frequency of seismic wave in the medium is the high frequency higher than upper limiting frequency, by the saturation crack rock mould
Type is embedded in the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under high frequency;
When the frequency of seismic wave in the medium is the intermediate frequency between the lower frequency limit and the upper limiting frequency,
On fracture surface direction, by the saturation crack rock model insertion into the saturation porous rocks model, to obtain
The medium is taken under intermediate frequency in the compressional wave elastic modelling quantity on fracture surface direction;
According to formulaObtain the 3rd coefficient of elasticity of the medium under intermediate frequency
Matrix Cij(ω);
Wherein,For first coefficient of elasticity matrix,For first coefficient of elasticity matrix, R (ω) is the frequency of the medium
Rate rely on function andC33(ω) is the medium under intermediate frequency on fracture surface direction
Compressional wave elastic modelling quantity,It is the medium under low frequency in the compressional wave elastic modelling quantity on fracture surface direction,For height
The medium is in the compressional wave elastic modelling quantity on fracture surface direction under frequency.
The speed and decay analog analysing method of the seismic wave of the embodiment of the present application, the elasticity system according to the medium
Matrix number determines the complex velocity that seismic wave is propagated in the medium, including:
3rd elastic coefficient matrix of the medium under intermediate frequency is substituted into below equation and obtains seismic wave in the medium
The complex velocity of propagation:
Wherein, D
=[(C11(ω)-C44(ω))sin2θ-(C33(ω)-C44(ω))cos2θ]2+4(C13(ω)+C44(ω))2sin2θcos2θ, ω
For the frequency of seismic wave in the medium, θ is the incidence angle of seismic wave in the medium, C11(ω) is the medium flat
Compressional wave elastic modelling quantity, C on row fracture surface direction44(ω) be the medium in the shear wave elastic modelling quantity on fracture surface direction,
C13(ω) be the medium with fracture surface into the elastic modelling quantity on 45 angular direction, vp(ω, θ) be seismic wave in the medium
Compressional wave complex velocity, vsv(ω, θ) is the shear wave complex velocity of seismic wave in the medium.
The speed and decay analog analysing method of the seismic wave of the embodiment of the present application, it is described that institute is determined according to the complex velocity
The speed that the inverse quality factor and the seismic wave of matter are propagated in the medium is given an account of, including:
According to formulaIt is determined that described
The speed that the inverse quality factor and the seismic wave of medium are propagated in the medium;
Wherein, Vp' (ω, θ) be the velocity of longitudinal wave of seismic wave in the medium, Vs'v(ω, θ) is seismic wave described
Shear wave complex velocity in medium, Qp -1(ω, θ) is the inverse quality factor of compressional wave in the medium, Qsv -1(ω, θ) exists for shear wave
Inverse quality factor in the medium, real represents real part, and image represents imaginary part.
On the other hand, the embodiment of the present application additionally provides the speed and decay simulation analytical equipment of a kind of seismic wave, including:
Petrophysical model sets up module, saturation porous rocks model and saturation crack rock mould for setting up medium in work area
Type;
Coefficient of elasticity acquisition module, in the case where imposing a condition, the saturation crack rock model insertion to be satisfied to described
In porous rocks model, to obtain the elastic coefficient matrix of the medium;The elastic coefficient matrix exists including the medium
The coefficient of elasticity of all directions;
Complex velocity acquisition module, for determining that seismic wave is passed in the medium according to the elastic coefficient matrix of the medium
The complex velocity broadcast;
Speed and inverse quality factor acquisition module, for according to the complex velocity determine the inverse quality factor of the medium with
And the speed that the seismic wave is propagated in the medium.
The speed and decay simulation analytical equipment of the seismic wave of the embodiment of the present application, it is described to set up medium in target work area
Saturation porous rocks model and saturation crack rock model, including:
Determine the elastic modelling quantity and density of solid dielectric in work area;The solid dielectric includes the first solid dielectric and second
Solid dielectric;
The porosity of first solid dielectric is determined, after interstitial space is introduced into first solid dielectric, to
Filling first fluid in the interstitial space, forms saturation porous rocks model;The first fluid has the first bulk modulus
With the first density;
The crevice volume percentage of second solid dielectric is determined, particular volume is being introduced into second solid dielectric
After the fracture spaces of product module amount, second fluid is filled into the fracture spaces, saturation crack rock model is formed;Described second
Fluid has the second bulk modulus and the second density.
The speed and decay simulation analytical equipment of the seismic wave of the embodiment of the present application, it is described in the case where imposing a condition, will be described
Saturation crack rock model insertion is into the saturation porous rocks model, to obtain the elastic coefficient matrix of the medium, bag
Include:
When the frequency of seismic wave in the medium is the low frequency less than lower frequency limit, by the saturation crack rock mould
Type is embedded in the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under low frequency;
When the frequency of seismic wave in the medium is the high frequency higher than upper limiting frequency, by the saturation crack rock mould
Type is embedded in the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under high frequency;
When the frequency of seismic wave in the medium is the intermediate frequency between the lower frequency limit and the upper limiting frequency,
On fracture surface direction, by the saturation crack rock model insertion into the saturation porous rocks model, to obtain
The medium is taken under intermediate frequency in the compressional wave elastic modelling quantity on fracture surface direction;
According to formulaObtain the 3rd coefficient of elasticity of the medium under intermediate frequency
Matrix Cij(ω);
Wherein,For first coefficient of elasticity matrix,For first coefficient of elasticity matrix, R (ω) is the medium
Frequency rely on function andC33(ω) is the medium under intermediate frequency on fracture surface direction
Compressional wave elastic modelling quantity,It is the medium under low frequency in the compressional wave elastic modelling quantity on fracture surface direction,For
The medium is in the compressional wave elastic modelling quantity on fracture surface direction under high frequency.
The speed and decay simulation analytical equipment of the seismic wave of the embodiment of the present application, the elasticity system according to the medium
Matrix number determines the complex velocity that seismic wave is propagated in the medium, including:
3rd elastic coefficient matrix of the medium under intermediate frequency is substituted into below equation and obtains seismic wave in the medium
The complex velocity of propagation:
Wherein, D
=[(C11(ω)-C44(ω))sin2θ-(C33(ω)-C44(ω))cos2θ]2+4(C13(ω)+C44(ω))2sin2θcos2θ, ω
For the frequency of seismic wave in the medium, θ is the incidence angle of seismic wave in the medium, C11(ω) is the medium flat
Compressional wave elastic modelling quantity, C on row fracture surface direction44(ω) be the medium in the shear wave elastic modelling quantity on fracture surface direction,
C13(ω) be the medium with fracture surface into the elastic modelling quantity on 45 angular direction, vp(ω, θ) be seismic wave in the medium
Compressional wave complex velocity, vsv(ω, θ) is the shear wave complex velocity of seismic wave in the medium.
The speed and decay simulation analytical equipment of the seismic wave of the embodiment of the present application, it is described that institute is determined according to the complex velocity
The speed that the inverse quality factor and the seismic wave of matter are propagated in the medium is given an account of, including:
According to formulaIt is determined that described
The speed that the inverse quality factor and the seismic wave of medium are propagated in the medium;
Wherein, Vp' (ω, θ) be the velocity of longitudinal wave of seismic wave in the medium, Vs'v(ω, θ) is seismic wave described
Shear wave complex velocity in medium, Qp -1(ω, θ) is the inverse quality factor of compressional wave in the medium, Qsv -1(ω, θ) exists for shear wave
Inverse quality factor in the medium, real represents real part, and image represents imaginary part.
The technical scheme provided from above the embodiment of the present application, the embodiment of the present application initially sets up medium in work area
Saturation porous rocks model and saturation crack rock model;Secondly in the case where imposing a condition, by saturation crack rock model insertion extremely
In saturation porous rocks model, to obtain the elastic coefficient matrix of medium;Elastic coefficient matrix includes medium in all directions
Coefficient of elasticity;Then the complex velocity that seismic wave is propagated in media as well is determined according to the elastic coefficient matrix of medium;So as to according to multiple
Speed determines the speed that the inverse quality factor and seismic wave of medium are propagated in the medium, wherein being reflected against quality factor
The attenuation degree of the energy of seismic wave in media as well, therefore, the embodiment of the present application can consider two kinds of attenuating mechanisms simultaneously over the ground
Seismic velocity and the influence of decay.Therefore, the analysis result that the embodiment of the present application is obtained can more objectively reflect that underground is situated between
Confrontation seimic wave velocity and the influence of energy attenuation.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, are not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.In the accompanying drawings:
Fig. 1 is the speed of the seismic wave of the embodiment of the application one and the method flow diagram for analog analysing method of decaying;
Fig. 2 a-2c are used by vertical fracture face direction attenuation of P-wave result (2b) in the embodiment of the application one and prior art
The result (2a) of model calculating, the comparison diagram of (2c);
Fig. 3 a-3b be respectively in the embodiment of the application one velocity of longitudinal wave when incidence angle is 30 ° with frequency change curve
With the change curve (3b) of incidence angle when (3a) and dominant frequency are 30Hz;
Fig. 4 a-4b be respectively in the embodiment of the application one attenuation of P-wave when incidence angle is 30 ° with frequency change curve
With the change curve (4b) of incidence angle when (4a) and dominant frequency are 30Hz;
Fig. 5 is change curve of the velocity of longitudinal wave in the embodiment of the application one with frequency and incidence angle;
Fig. 6 is change curve of the attenuation of P-wave in the embodiment of the application one with frequency and incidence angle;
Fig. 7 a-7b be respectively in the embodiment of the application one shear wave velocity when incidence angle is 30 ° with frequency change curve
With the change curve (7b) of incidence angle when (7a) and dominant frequency are 30Hz;
Fig. 8 a-8b be respectively in the embodiment of the application one shear wave decay when incidence angle is 30 ° with frequency change curve
With the change curve (8b) of incidence angle when (8a) and dominant frequency are 30Hz;
Fig. 9 is change curve of the shear wave velocity in the embodiment of the application one with frequency and incidence angle;
Figure 10 is change curve of the shear wave decay in the embodiment of the application one with frequency and incidence angle;
Figure 11 is the speed of the seismic wave of the embodiment of the application one and the structured flowchart of decay simulation analytical equipment.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.Based on the embodiment in the application, this area is common
The every other embodiment that technical staff is obtained under the premise of creative work is not made, should all belong to the application protection
Scope.
With reference to shown in Fig. 1, the speed and decay analog analysing method of the seismic wave of the embodiment of the present application comprise the following steps:
Step S101, the saturation porous rocks model for setting up medium in work area and saturation crack rock model.
In one embodiment of the application, the saturation porous rocks model and saturation for setting up medium in target work area
Crack rock model may include steps of:
First, the elastic modelling quantity and density p of solid dielectric in work area are determinedg;The solid dielectric is situated between including the first solid
Matter and the second solid dielectric.Wherein, the elastic modelling quantity of the solid dielectric can include bulk modulus KgWith modulus of shearing μg。
Then, it is determined that the porosity φ of first solid dielectricm, being introduced into first solid dielectric, hole is empty
Between after, fill first fluid into the interstitial space, form saturation porous rocks model;The first fluid has the first body
Product module amount KfbWith the first density pfb。
Then, it is determined that the crevice volume percentage υ of second solid dielectricc, draw into second solid dielectric
Enter designated volume modulus LcFracture spaces after, fill second fluid into the fracture spaces, form saturation crack rock mould
Type;The second fluid has the second bulk modulus KfcWith the second density pfc。
, can also be by introducing the solid modular ratio F=K of crack stream in another embodiment of the applicationfc/Lc, fracture parameter
Simplified, and form saturation crack rock model on this basis.
It should be noted that in the embodiment of the present application, saturation porous rocks model and saturation crack rock model have no
Sequence requirement.
Step S102, in the case where imposing a condition, by the saturation crack rock model insertion to the saturation porous rocks mould
In type, to obtain the elastic coefficient matrix of the medium;The elastic coefficient matrix includes bullet of the medium in all directions
Property coefficient.
It is described in the case where imposing a condition in one embodiment of the application, by the saturation crack rock model insertion extremely
In the saturation porous rocks model, to obtain the elastic coefficient matrix of the medium, it can include:
When the frequency of seismic wave in the medium is the low frequency less than lower frequency limit, by the saturation crack rock mould
Type is embedded in the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under low frequencyIn one embodiment of the application, when the frequency of seismic wave in the medium is less than lower limit frequently
During rate, the pressure differential between the crack in hole and the saturation crack rock model in the saturation porous rocks model becomes
In 0, now, it can obtain under this limiting condition, coefficient of elasticity of the medium in all directions.
When the frequency of seismic wave in the medium is the high frequency higher than upper limiting frequency, by the saturation crack rock mould
Type is embedded in the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under high frequencyIn one embodiment of the application, when the frequency of seismic wave in the medium is higher than the upper limit frequently
During rate, the pressure differential between the crack in hole and the saturation crack rock model in the saturation porous rocks model becomes
In infinity, now, it can obtain under this limiting condition, coefficient of elasticity of the medium in all directions.
When the frequency of seismic wave in the medium is the intermediate frequency between the lower frequency limit and the upper limiting frequency,
On fracture surface direction, by the saturation crack rock model insertion into the saturation porous rocks model, to obtain
The medium is taken under intermediate frequency in the compressional wave elastic modelling quantity on fracture surface direction;
According to formulaObtain the 3rd coefficient of elasticity of the medium under intermediate frequency
Matrix Cij(ω);
Wherein,For first coefficient of elasticity matrix,For first coefficient of elasticity matrix, R (ω) is the medium
Frequency rely on function andC33(ω) is the medium under intermediate frequency on fracture surface direction
Compressional wave elastic modelling quantity,It is the medium under low frequency in the compressional wave elastic modelling quantity on fracture surface direction,For
The medium is in the compressional wave elastic modelling quantity on fracture surface direction under high frequency.
, can be theoretical by the saturation crack rock model based on poroelasticity Backus in one embodiment of the application
It is embedded in the saturation porous rocks model.Wherein, poroelasticity Backus theories are on basis theoretical elastic Backus
On grow up, the demand to meet pore media modeling.
Step S103, the complex-velocity that seismic wave is propagated in the medium is determined according to the elastic coefficient matrix of the medium
Degree.
In one embodiment of the application, the elastic coefficient matrix according to the medium determines seismic wave described
The complex velocity of Propagation, can include:
3rd elastic coefficient matrix of the medium under intermediate frequency is substituted into below equation and obtains seismic wave in the medium
The complex velocity of propagation:
Wherein, D
=[(C11(ω)-C44(ω))sin2θ-(C33(ω)-C44(ω))cos2θ]2+4(C13(ω)+C44(ω))2sin2θcos2θ, ω
For the frequency of seismic wave in the medium, θ is the incidence angle of seismic wave in the medium, C11(ω) is the medium flat
Compressional wave elastic modelling quantity, C on row fracture surface direction44(ω) be the medium in the shear wave elastic modelling quantity on fracture surface direction,
C13(ω) be the medium with fracture surface into the elastic modelling quantity on 45 angular direction, vp(ω, θ) be seismic wave in the medium
Compressional wave complex velocity, vsv(ω, θ) is the shear wave complex velocity of seismic wave in the medium.
Step S104, according to the complex velocity determine that the inverse quality factor and the seismic wave of the medium are being given an account of
The speed propagated in matter.
In one embodiment of the application, it is described according to the complex velocity determine the medium inverse quality factor and
The speed that the seismic wave is propagated in the medium, can include:
According to formulaIt is determined that described
The speed that the inverse quality factor and the seismic wave of medium are propagated in the medium;
Wherein, Vp' (ω, θ) be the velocity of longitudinal wave of seismic wave in the medium, Vs'v(ω, θ) is seismic wave described
Shear wave complex velocity in medium, Qp -1(ω, θ) is the inverse quality factor of compressional wave in the medium, Qsv -1(ω, θ) exists for shear wave
Inverse quality factor in the medium, real represents real part, and image represents imaginary part.Wherein, inverse quality factor reflects seismic wave
Energy attenuation degree in media as well.
Because the embodiment of the present application has obtained the complete elastic coefficient matrix Cij of the medium under intermediate frequency, therefore, this Shen
Please example except result in velocity of longitudinal wave and decay with frequency change, moreover it is possible to calculate velocity of longitudinal wave and decay with incidence angle
Influence.That is, the embodiment of the present application can obtain P-wave And S the speed of the media and energy attenuation with frequencies omega and
The change of incidence angle θ.For example in one exemplary embodiment of the application, such as Fig. 3 a-3b, Fig. 4 a-4b, Fig. 5-Fig. 6, Fig. 7 a-
7b, Fig. 8 a-8b, and shown in Fig. 9-Figure 10.
As shown in Fig. 2 a- Fig. 2 c, in one exemplary embodiment of the application, Fig. 2 b are hanging down based on the embodiment of the present application
Straight fracture surface direction attenuation of P-wave result;Fig. 2 a are the attenuation of P-wave results using existing White models, and the model is only accounted for
First kind attenuating mechanism in medium produced by fluid properties difference, with the increase of frequency, its attenuation curve, which is presented, first to be increased
The process reduced again;Fig. 2 c are the attenuation of P-wave results using existing Brajanovski models, and the model is only accounted in medium
Equations of The Second Kind attenuating mechanism produced by crack and hole solid properties difference, with the increase of frequency, its attenuation curve is presented first
Rapid increase, the again gentle trend for increasing, being finally gradually reduced.And as can be seen that with the increase of frequency, being based on from Fig. 2 b
The attenuation of P-wave that the model that the present invention is set up is obtained is while have the curve feature of two class attenuating mechanisms.When solid model amount is flowed in crack
When smaller than F, the curve feature of attenuation of the first kind mechanism is presented;With F gradually increase, attenuation curve from the first type by
Gradually converted to second of type;When F is larger, the curve feature of attenuation of the second kind mechanism is presented.
As can be seen that the embodiment of the present application can consider two kinds of attenuating mechanisms to earthquake simultaneously from above-mentioned comparative analysis
Wave velocity and the influence of decay, and P- and S-wave velocity can simultaneously be analyzed with the change of frequency and incidence angle, it is actual to answer
It is wider with scope.
Although procedures described above flow includes the multiple operations occurred with particular order, it should however be appreciated that understand,
These processes can include more or less operations, and these operations can sequentially be performed or performed parallel (such as using parallel
Processor or multi-thread environment).
Shown in Figure 11, the speed and decay simulation analytical equipment of the seismic wave of the embodiment of the present application can include:
Petrophysical model sets up module 111, and the saturation porous rocks model and saturation that can be used for setting up medium in work area split
Stitch petrophysical model;
Coefficient of elasticity acquisition module 112, can be used in the case where imposing a condition, by the saturation crack rock model insertion extremely
In the saturation porous rocks model, to obtain the elastic coefficient matrix of the medium;The elastic coefficient matrix includes described
Coefficient of elasticity of the medium in all directions;
Complex velocity acquisition module 113, can be used for determining seismic wave described according to the elastic coefficient matrix of the medium
The complex velocity of Propagation;
Speed and inverse quality factor acquisition module 114, can be used for the inverse product that the medium is determined according to the complex velocity
The speed that prime factor and the seismic wave are propagated in the medium.
The speed of the speed of the seismic wave of the embodiment of the present application and the seismic wave of decay simulation analytical equipment and above-described embodiment
Degree and decay analog analysing method correspondence, therefore, the speed and decay simulation analytical equipment for being related to the seismic wave of the application are thin
Section, refers to the speed and decay analog analysing method of the seismic wave of above-described embodiment, will not be repeated here.
For convenience of description, it is divided into various units during description apparatus above with function to describe respectively.Certainly, this is being implemented
The function of each unit can be realized in same or multiple softwares and/or hardware during application.
It should be understood by those skilled in the art that, embodiments of the invention can be provided as method, system or computer program
Product.Therefore, the present invention can be using the reality in terms of complete hardware embodiment, complete software embodiment or combination software and hardware
Apply the form of example.Moreover, the present invention can be used in one or more computers for wherein including computer usable program code
The computer program production that usable storage medium is implemented on (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.)
The form of product.
The present invention is the flow with reference to method according to embodiments of the present invention, equipment (system) and computer program product
Figure and/or block diagram are described.It should be understood that every one stream in flow chart and/or block diagram can be realized by computer program instructions
Journey and/or the flow in square frame and flow chart and/or block diagram and/or the combination of square frame.These computer programs can be provided
The processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices is instructed to produce
A raw machine so that produced by the instruction of computer or the computing device of other programmable data processing devices for real
The device for the function of being specified in present one flow of flow chart or one square frame of multiple flows and/or block diagram or multiple square frames.
These computer program instructions, which may be alternatively stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that the instruction being stored in the computer-readable memory, which is produced, to be included referring to
Make the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one square frame of block diagram or
The function of being specified in multiple square frames.
These computer program instructions can be also loaded into computer or other programmable data processing devices so that in meter
Series of operation steps is performed on calculation machine or other programmable devices to produce computer implemented processing, thus in computer or
The instruction performed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in individual square frame or multiple square frames.
In a typical configuration, computing device includes one or more processors (CPU), input/output interface, net
Network interface and internal memory.
Internal memory potentially includes the volatile memory in computer-readable medium, random access memory (RAM) and/or
The forms such as Nonvolatile memory, such as read-only storage (ROM) or flash memory (flash RAM).Internal memory is computer-readable medium
Example.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology realizes information Store.Information can be computer-readable instruction, data structure, the module of program or other data.
The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM (SRAM), moved
State random access memory (DRAM), other kinds of random access memory (RAM), read-only storage (ROM), electric erasable
Programmable read only memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc read-only storage (CD-ROM),
Digital versatile disc (DVD) or other optical storages, magnetic cassette tape, the storage of tape magnetic rigid disk or other magnetic storage apparatus
Or any other non-transmission medium, the information that can be accessed by a computing device available for storage.Define, calculate according to herein
Machine computer-readable recording medium does not include temporary computer readable media (transitory media), such as data-signal and carrier wave of modulation.
It should also be noted that, term " comprising ", "comprising" or its any other variant are intended to nonexcludability
Comprising so that process, method, commodity or equipment including a series of key elements are not only including those key elements, but also wrap
Include other key elements being not expressly set out, or also include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that wanted including described
Also there is other identical element in process, method, commodity or the equipment of element.
It will be understood by those skilled in the art that embodiments herein can be provided as method, system or computer program product.
Therefore, the application can be using the embodiment in terms of complete hardware embodiment, complete software embodiment or combination software and hardware
Form.Deposited moreover, the application can use to can use in one or more computers for wherein including computer usable program code
The shape for the computer program product that storage media is implemented on (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.)
Formula.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes execution particular task or the routine for realizing particular abstract data type, program, object, group
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Remote processing devices connected by communication network perform task.In a distributed computing environment, program module can be with
Positioned at including in the local and remote computer-readable storage medium including storage device.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.It is real especially for system
Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
Embodiments herein is the foregoing is only, the application is not limited to.For those skilled in the art
For, the application can have various modifications and variations.It is all any modifications made within spirit herein and principle, equivalent
Replace, improve etc., it should be included within the scope of claims hereof.
Claims (10)
1. the speed and decay analog analysing method of a kind of seismic wave, it is characterised in that comprise the following steps:
Set up the saturation porous rocks model and saturation crack rock model of medium in work area;
In the case where imposing a condition, by the saturation crack rock model insertion into the saturation porous rocks model, to obtain
Give an account of the elastic coefficient matrix of matter;The elastic coefficient matrix includes coefficient of elasticity of the medium in all directions;
The complex velocity that seismic wave is propagated in the medium is determined according to the elastic coefficient matrix of the medium;
The speed that the inverse quality factor and the seismic wave of the medium are propagated in the medium is determined according to the complex velocity
Degree.
2. the speed and decay analog analysing method of seismic wave according to claim 1, it is characterised in that described to set up mesh
The saturation porous rocks model and saturation crack rock model of medium in work area are marked, including:
Determine the elastic modelling quantity and density of solid dielectric in work area;The solid dielectric includes the first solid dielectric and the second solid
Medium;
The porosity of first solid dielectric is determined, after interstitial space is introduced into first solid dielectric, to described
Filling first fluid in interstitial space, forms saturation porous rocks model;The first fluid has the first bulk modulus and the
One density;
The crevice volume percentage of second solid dielectric is determined, particular volume product module is being introduced into second solid dielectric
After the fracture spaces of amount, second fluid is filled into the fracture spaces, saturation crack rock model is formed;The second fluid
With the second bulk modulus and the second density.
3. the speed and decay analog analysing method of seismic wave according to claim 1, it is characterised in that described in setting
Under the conditions of, by the saturation crack rock model insertion into the saturation porous rocks model, to obtain the bullet of the medium
Property coefficient matrix, including:
It is when the frequency of seismic wave in the medium is the low frequency less than lower frequency limit, the saturation crack rock model is embedding
Enter into the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under low frequency;
It is when the frequency of seismic wave in the medium is the high frequency higher than upper limiting frequency, the saturation crack rock model is embedding
Enter into the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under high frequency;
When the frequency of seismic wave in the medium is the intermediate frequency between the lower frequency limit and the upper limiting frequency, hanging down
Directly on fracture surface direction, by the saturation crack rock model insertion into the saturation porous rocks model, with acquisition
The medium is in the compressional wave elastic modelling quantity on fracture surface direction under frequency;
According to formulaObtain the 3rd elastic coefficient matrix of the medium under intermediate frequency
Cij(ω);
Wherein,For first coefficient of elasticity matrix,For first coefficient of elasticity matrix, R (ω) for the medium frequency according to
Rely function andC33(ω) is the medium under intermediate frequency in the compressional wave on fracture surface direction
Elastic modelling quantity,It is the medium under low frequency in the compressional wave elastic modelling quantity on fracture surface direction,For under high frequency
The medium is in the compressional wave elastic modelling quantity on fracture surface direction.
4. the speed and decay analog analysing method of seismic wave according to claim 3, it is characterised in that described according to institute
The elastic coefficient matrix for giving an account of matter determines the complex velocity that seismic wave is propagated in the medium, including:
3rd elastic coefficient matrix of the medium under intermediate frequency is substituted into below equation acquisition seismic wave to propagate in the medium
Complex velocity:
Wherein, D=[(C11(ω)-C44(ω))sin2θ-(C33(ω)-C44(ω))cos2θ]2+4(C13(ω)+C44(ω))2sin2θcos2θ, ω are the frequency of seismic wave in the medium, and θ is the incidence angle of seismic wave in the medium, C11(ω)
For the medium on the direction of parallel fracture face compressional wave elastic modelling quantity, C44(ω) is the medium on fracture surface direction
Shear wave elastic modelling quantity, C13(ω) be the medium with fracture surface into the elastic modelling quantity on 45 angular direction, vp(ω, θ) is earthquake
The compressional wave complex velocity of ripple in the medium, vsv(ω, θ) is the shear wave complex velocity of seismic wave in the medium.
5. the speed and decay analog analysing method of seismic wave according to claim 4, it is characterised in that described according to institute
State complex velocity and determine the speed that the inverse quality factor and the seismic wave of the medium are propagated in the medium, including:
According to formulaDetermine the medium
The speed that inverse quality factor and the seismic wave are propagated in the medium;
Wherein, Vp' (ω, θ) be the velocity of longitudinal wave of seismic wave in the medium, Vs'v(ω, θ) is seismic wave in the medium
In shear wave complex velocity, Qp -1(ω, θ) is the inverse quality factor of compressional wave in the medium, Qsv -1(ω, θ) is shear wave described
Inverse quality factor in medium, real represents real part, and image represents imaginary part.
6. the speed and decay simulation analytical equipment of a kind of seismic wave, it is characterised in that including:
Petrophysical model sets up module, saturation porous rocks model and saturation crack rock model for setting up medium in work area;
Coefficient of elasticity acquisition module, in the case where imposing a condition, by the saturation crack rock model insertion to the saturation hole
In gap petrophysical model, to obtain the elastic coefficient matrix of the medium;The elastic coefficient matrix includes the medium at each
The coefficient of elasticity in direction;
Complex velocity acquisition module, for determining what seismic wave was propagated in the medium according to the elastic coefficient matrix of the medium
Complex velocity;
Speed and inverse quality factor acquisition module, inverse quality factor and institute for determining the medium according to the complex velocity
State the speed that seismic wave is propagated in the medium.
7. the speed and decay simulation analytical equipment of seismic wave according to claim 6, it is characterised in that described to set up mesh
The saturation porous rocks model and saturation crack rock model of medium in work area are marked, including:
Determine the elastic modelling quantity and density of solid dielectric in work area;The solid dielectric includes the first solid dielectric and the second solid
Medium;
The porosity of first solid dielectric is determined, after interstitial space is introduced into first solid dielectric, to described
Filling first fluid in interstitial space, forms saturation porous rocks model;The first fluid has the first bulk modulus and the
One density;
The crevice volume percentage of second solid dielectric is determined, particular volume product module is being introduced into second solid dielectric
After the fracture spaces of amount, second fluid is filled into the fracture spaces, saturation crack rock model is formed;The second fluid
With the second bulk modulus and the second density.
8. the speed and decay simulation analytical equipment of seismic wave according to claim 6, it is characterised in that described in setting
Under the conditions of, by the saturation crack rock model insertion into the saturation porous rocks model, to obtain the bullet of the medium
Property coefficient matrix, including:
It is when the frequency of seismic wave in the medium is the low frequency less than lower frequency limit, the saturation crack rock model is embedding
Enter into the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under low frequency;
It is when the frequency of seismic wave in the medium is the high frequency higher than upper limiting frequency, the saturation crack rock model is embedding
Enter into the saturation porous rocks model, to obtain the first coefficient of elasticity matrix of the medium under high frequency;
When the frequency of seismic wave in the medium is the intermediate frequency between the lower frequency limit and the upper limiting frequency, hanging down
Directly on fracture surface direction, by the saturation crack rock model insertion into the saturation porous rocks model, with acquisition
The medium is in the compressional wave elastic modelling quantity on fracture surface direction under frequency;
According to formulaObtain the 3rd elastic coefficient matrix of the medium under intermediate frequency
Cij(ω);
Wherein,For first coefficient of elasticity matrix,For first coefficient of elasticity matrix, R (ω) for the medium frequency according to
Rely function andC33(ω) is the medium under intermediate frequency in the compressional wave on fracture surface direction
Elastic modelling quantity,It is the medium under low frequency in the compressional wave elastic modelling quantity on fracture surface direction,For under high frequency
The medium is in the compressional wave elastic modelling quantity on fracture surface direction.
9. the speed and decay simulation analytical equipment of seismic wave according to claim 8, it is characterised in that described according to institute
The elastic coefficient matrix for giving an account of matter determines the complex velocity that seismic wave is propagated in the medium, including:
3rd elastic coefficient matrix of the medium under intermediate frequency is substituted into below equation acquisition seismic wave to propagate in the medium
Complex velocity:
Wherein, D=[(C11(ω)-C44(ω))sin2θ-(C33(ω)-C44(ω))cos2θ]2+4(C13(ω)+C44(ω))2sin2θcos2θ, ω are the frequency of seismic wave in the medium, and θ is the incidence angle of seismic wave in the medium, C11(ω)
For the medium on the direction of parallel fracture face compressional wave elastic modelling quantity, C44(ω) is the medium on fracture surface direction
Shear wave elastic modelling quantity, C13(ω) be the medium with fracture surface into the elastic modelling quantity on 45 angular direction, vp(ω, θ) is earthquake
The compressional wave complex velocity of ripple in the medium, vsv(ω, θ) is the shear wave complex velocity of seismic wave in the medium.
10. the speed and decay simulation analytical equipment of seismic wave according to claim 9, it is characterised in that the basis
The complex velocity determines the speed that the inverse quality factor and the seismic wave of the medium are propagated in the medium, including:
According to below equationIt is determined that described
The speed that the inverse quality factor and the seismic wave of medium are propagated in the medium;
Wherein, Vp' (ω, θ) be the velocity of longitudinal wave of seismic wave in the medium, Vs'v(ω, θ) is seismic wave in the medium
In shear wave complex velocity, Qp -1(ω, θ) is the inverse quality factor of compressional wave in the medium, Qsv -1(ω, θ) is shear wave described
Inverse quality factor in medium, real represents real part, and image represents imaginary part.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710242793.1A CN106959469A (en) | 2017-04-14 | 2017-04-14 | Method and device for simulating and analyzing velocity and attenuation of seismic wave |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710242793.1A CN106959469A (en) | 2017-04-14 | 2017-04-14 | Method and device for simulating and analyzing velocity and attenuation of seismic wave |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106959469A true CN106959469A (en) | 2017-07-18 |
Family
ID=59484385
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710242793.1A Pending CN106959469A (en) | 2017-04-14 | 2017-04-14 | Method and device for simulating and analyzing velocity and attenuation of seismic wave |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106959469A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109116420A (en) * | 2018-10-16 | 2019-01-01 | 河海大学 | A kind of pore media velocity of longitudinal wave containing crack and decaying prediction technique |
CN109143340A (en) * | 2018-08-20 | 2019-01-04 | 中国海洋石油集团有限公司 | A kind of visco-elastic medium Simulating Seismic Wave method and system based on normal Q model |
CN109521463A (en) * | 2017-09-20 | 2019-03-26 | 中国石油化工股份有限公司 | Determine the method and system of the best earthquake excitation li thology of igneous rock near surface |
CN111751886A (en) * | 2020-08-12 | 2020-10-09 | 重庆地质矿产研究院 | Shale gas reservoir fracture modeling method based on microseism monitoring data |
CN112363213A (en) * | 2020-12-03 | 2021-02-12 | 中国石油大学(华东) | Frequency dispersion and attenuation calculation method for layered rock containing randomly distributed cracks |
CN112799127A (en) * | 2020-12-30 | 2021-05-14 | 成都理工大学 | Seismic wave frequency dispersion and attenuation numerical calculation method considering non-uniform difference of seepage of fractured pore rock |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104570084A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Cross-scale seismic rock physical attenuation model and method for predicating attenuation and dispersion |
US9158017B2 (en) * | 2011-03-22 | 2015-10-13 | Seoul National University R&Db Foundation | Seismic imaging apparatus utilizing macro-velocity model and method for the same |
CN106054242A (en) * | 2016-05-04 | 2016-10-26 | 中国地质大学(北京) | Three dimensional anisotropy attenuation medium wave field simulation method |
CN106125135A (en) * | 2016-06-12 | 2016-11-16 | 成都理工大学 | Gas-bearing sandstone reservoir seismic response method for numerical simulation based on petrophysical model |
CN106154331A (en) * | 2016-06-29 | 2016-11-23 | 中国石油化工股份有限公司 | Orthogonal medium Simulation of Seismic Wave frequency dispersion drawing method |
-
2017
- 2017-04-14 CN CN201710242793.1A patent/CN106959469A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9158017B2 (en) * | 2011-03-22 | 2015-10-13 | Seoul National University R&Db Foundation | Seismic imaging apparatus utilizing macro-velocity model and method for the same |
CN104570084A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Cross-scale seismic rock physical attenuation model and method for predicating attenuation and dispersion |
CN106054242A (en) * | 2016-05-04 | 2016-10-26 | 中国地质大学(北京) | Three dimensional anisotropy attenuation medium wave field simulation method |
CN106125135A (en) * | 2016-06-12 | 2016-11-16 | 成都理工大学 | Gas-bearing sandstone reservoir seismic response method for numerical simulation based on petrophysical model |
CN106154331A (en) * | 2016-06-29 | 2016-11-23 | 中国石油化工股份有限公司 | Orthogonal medium Simulation of Seismic Wave frequency dispersion drawing method |
Non-Patent Citations (2)
Title |
---|
LIYUN KONG ET AL.: "Effect of fracture fill on frequency-dependent anisotropy of fractured porous rocks", 《GEOPHYSICAL PROSPECTING》 * |
孔丽云 等: "裂缝诱导TTI双孔隙介质波场传播特征", 《物理学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109521463A (en) * | 2017-09-20 | 2019-03-26 | 中国石油化工股份有限公司 | Determine the method and system of the best earthquake excitation li thology of igneous rock near surface |
CN109143340A (en) * | 2018-08-20 | 2019-01-04 | 中国海洋石油集团有限公司 | A kind of visco-elastic medium Simulating Seismic Wave method and system based on normal Q model |
CN109143340B (en) * | 2018-08-20 | 2020-03-10 | 中国海洋石油集团有限公司 | Viscoelastic medium seismic wave simulation method and system based on constant Q model |
CN109116420A (en) * | 2018-10-16 | 2019-01-01 | 河海大学 | A kind of pore media velocity of longitudinal wave containing crack and decaying prediction technique |
CN111751886A (en) * | 2020-08-12 | 2020-10-09 | 重庆地质矿产研究院 | Shale gas reservoir fracture modeling method based on microseism monitoring data |
CN112363213A (en) * | 2020-12-03 | 2021-02-12 | 中国石油大学(华东) | Frequency dispersion and attenuation calculation method for layered rock containing randomly distributed cracks |
CN112799127A (en) * | 2020-12-30 | 2021-05-14 | 成都理工大学 | Seismic wave frequency dispersion and attenuation numerical calculation method considering non-uniform difference of seepage of fractured pore rock |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106959469A (en) | Method and device for simulating and analyzing velocity and attenuation of seismic wave | |
US20170038489A1 (en) | Fracture-Size-Correlated Aperture Mapping for Localized Porosity and Permeability Determination | |
AU2008205064B8 (en) | Systems and methods for selectively imaging objects in a display of multiple three-dimensional data-objects | |
US10495768B2 (en) | Method of operating a data-processing system for the simulation of the acoustic wave propagation in the transversely isotropic media comprising an hydrocarbon reservoir | |
EA025919B1 (en) | Systems and methods for modeling 3d geological structures | |
US20110054857A1 (en) | Gridless geological modeling | |
CN103119471A (en) | System and method for generating images of subsurface structures | |
CN107272058A (en) | Imaging method, imaging device and computer-readable storage medium | |
CN108983285A (en) | A kind of a variety of source wavefield analogy methods and device based on moment tensor | |
CN108415073B (en) | Angle domain back scattering offset imaging method and device | |
CN106896406A (en) | Method and device for predicting high-quality reservoir stratum based on impedance domain attributes | |
US10718188B2 (en) | Method for evaluation of fluid transport properties in heterogenous geological formation | |
CN108897036A (en) | A kind of seismic data processing technique and device | |
CN109597964A (en) | Theoretical flow through oil reservoir parameter value determines method, system, equipment, readable medium | |
CN111323814B (en) | Method for quantitatively determining water saturation of sandstone reservoir based on rock physical template | |
WO2013033651A1 (en) | Full elastic wave equation for 3d data processing on gpgpu | |
CN109738944B (en) | Wide-angle reflection-based seismic acquisition parameter determination method and device | |
CN109521470B (en) | Method for analyzing influence of geological structure on seismic inversion crack density | |
US10454713B2 (en) | Domain decomposition using a multi-dimensional spacepartitioning tree | |
CN112764105B (en) | HTI medium quasi-longitudinal wave forward simulation method and device, storage medium and processor | |
CN107390271B (en) | The vibration amplitude compensation method and device of seismic data | |
Yang et al. | Unstructured mesh based elastic wave modelling on GPU: a double-mesh grid method | |
CN109581484A (en) | P wave component multiple wave general ambient light degree index analysis method and device | |
Chou et al. | A crystal defect theory approach to faulting in geological structures | |
AU2013406187B2 (en) | Geocellular modeling |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170718 |
|
WD01 | Invention patent application deemed withdrawn after publication |