CN109324343A - A kind of analogy method and system of thin layer displacement multi-wave seismic wave field - Google Patents
A kind of analogy method and system of thin layer displacement multi-wave seismic wave field Download PDFInfo
- Publication number
- CN109324343A CN109324343A CN201710647906.6A CN201710647906A CN109324343A CN 109324343 A CN109324343 A CN 109324343A CN 201710647906 A CN201710647906 A CN 201710647906A CN 109324343 A CN109324343 A CN 109324343A
- Authority
- CN
- China
- Prior art keywords
- wave
- thin layer
- frequency domain
- seismic
- cos
- 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 32
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 21
- 230000005540 biological transmission Effects 0.000 claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 238000004590 computer program Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 22
- 238000006243 chemical reaction Methods 0.000 abstract description 5
- 230000009466 transformation Effects 0.000 abstract description 3
- 230000005428 wave function Effects 0.000 abstract description 2
- 239000010410 layer Substances 0.000 description 96
- 239000003245 coal Substances 0.000 description 15
- 238000005516 engineering process Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 238000005096 rolling process Methods 0.000 description 4
- 101100234408 Danio rerio kif7 gene Proteins 0.000 description 3
- 101100221620 Drosophila melanogaster cos gene Proteins 0.000 description 3
- 101100398237 Xenopus tropicalis kif11 gene Proteins 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000005266 casting Methods 0.000 description 3
- 101150118300 cos gene Proteins 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 241000219000 Populus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000000985 reflectance spectrum Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000003595 spectral 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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
Disclose the analogy method and system of a kind of thin layer displacement multi-wave seismic wave field.Method includes the following steps: 1) obtaining frequency domain thin layer based on thin layer displacement-reflection transmission matrix equation is displaced PP wave reflection coefficients RppWith PS wave reflection coefficients Rps；2) FFT transform is carried out to seismic wavelet w (t), obtained frequency domain seismic wavelet w (f)；3) it obtains frequency domain PP and involves PS wave reflection seismic wave field；4) frequency domain PP obtained to step 3) involves PS wave reflection seismic wave field and obtain time-domain PP against FFT transform and involve PS wave reflection seismic wave field.The present invention precisely simulation based on the thin-layered medium reflected P P of displacement wave function can involve PS wave, and reflected wave field includes that conversion involves interbed multiple information, can comprehensively, accurately reflect that thin layer acts on the transformation of seismic wavelet.
Description
Technical field
The present invention relates to seismic exploration technique fields, more particularly, to a kind of mould of thin layer displacement multi-wave seismic wave field
Quasi- method and system.
Background technique
With deepening continuously for the developments of resources such as China's oil gas field, coalfield, thin (mutual) reservoir, lithology, carbonate rock, very
Advising the complex reservoirs such as oil gas and coal seam reservoirs becomes the emphasis of exploration.In east China area, most of Mesozoic and Cenozoic time land
Phase oil-gas bearing basin accompanies a small amount of carbonate rock, shale and gypsum-salt bed mostly based on thin layer sand, mud stone deposition, these
Formation lithology and thickness cross directional variations are larger, and thickness is well below the vertical resolution that common seismic is explored, with thin layer
Form exist.China coal measure strata overwhelming majority coal seam and country rock physical difference are larger, and coal seam top, bed rock layer are good wave
Impedance interface can form stronger back wave detection.Since coal seam thickness is mostly between 2~10m, and often with thin interbed
Form exists, and the non-single interface of reflex response generates, but push up bottom reflection, interbed multiple, converted wave be formed by stacking it is compound
Wave.Therefore, qualitative or quantitative analysis is carried out to the seismic wave field of thin layer, thin interbed, is the basis of oil gas, coal seismic prospecting.
Wherein, the forward numerical simulation of seismic wave field is propagation characteristic, clothes of the full appreciation seismic wave in complex dielectrics
It is engaged in explanation observation data and gets a kind of architectonic effective means clear.Seismic forward simulation mainly include ray casting,
Finite difference (member) method and reflectivity method.Cost is relatively low for ray casting calculating, is capable of handling the geologic body of complicated form；Disadvantage
It is that precision is low and be difficult to handle the coupled problem in direction anisotropy medium between different type shear wave.Finite difference calculus application
Range can relatively accurately simulate the seismic wave field in any non-uniform dielectric than wide；The disadvantage is that higher cost is calculated, it is difficult
To directly apply to the inverting of seismic wave field, while needing using absorbing boundary condition.Reflectivity method is a kind of numerical transformation side
Method, it be realize layered half space medium in seismic full-field field simulation most effectual way, precision be higher than finite difference calculus and
Ray casting calculates cost therebetween；Shortcoming is to be only applicable to vertical heterogeneous media.For the mould of thin layer
The characteristics of quasi-, reflectivity method is able to satisfy required precision, and high with computational efficiency, may be directly applied to inverting.
Seimic wave propagation theory and wave-field simulation technology are constantly subjected to the attention of people in thin layer.Meissner and
Meixner (Meissner R, Meixner E.Deformation of seismic wavelets by thin layers
And layered boundaries.GeophysProsp, 1969,17 (1): 1-27) give the thin of specified type back wave
Layer reflection coefficient formula, Liu and Schmitt give thin layer top, bottom reflection superposition approximate acoustic reflection coefficient (Liu Y,
Schmitt D.R.Amplitude and AVO responses of a single thin bed.Geophysics, 2003,
68 (4): 1161-1168), Pan and Kristopher discuss the AVO feature of thin layer elastic wave on the basis of Liu et al. people
(Pan,W.,and A.Kristopher,2013,AVO/AVF analysis of thin-bed in elastic media:
83rd Annual International Meeting,SEG,Expanded Abstracts,373–377.).Boulez Hough this
Ji He gives the bit shift propogator matrix of layered medium, and Wang Enhua etc., Chen and Liu are used for thin-layered medium situation, and
(the conspicuous work of Paderewski, poplar instruction people translate to JI.M. Boulez to the seismic wave AVO feature and wave-field simulation for specifically discussing thin layer bit shift suddenly
Beijing wave in layered medium: Science Press, 1960；The layered medium reflection coefficient that Wang grace China is suitable for AVO analysis is public
Formula and model calculate geophysical prospecting for oil, 1991,4 (26): 425~429；Wang Enhua, He Zhenhua, Li Qing loyalty are based on thin
The reflection coefficient spectral theory and forward modeling Chengdu University of Technology journal of layer, 2001,28 (1): 70~74；Chen T S, Liu
Y.Multi-component AVO response of thin beds based on reflectance spectrum
Theory.Applied Geophysics, 2006,3 (1): 27-36).But the above thin layer seismic wave field communication theory and simulation skill
Art cannot achieve the simulation of thin layer displacement multi-wave seismic wave field.More waves are displaced therefore, it is necessary to develop one kind and can simulate thin layer
The method and system of seismic wave field.
The information for being disclosed in background of invention part is merely intended to deepen the reason to general background technique of the invention
Solution, and it is known to those skilled in the art existing to be not construed as recognizing or imply that the information is constituted in any form
Technology.
Summary of the invention
For the above problem in the prior art, it is an object of the invention to be directed to thin-layered medium, it is anti-to provide thin layer displacement
Transmission matrix equation is penetrated, and realizes the forward simulation of the displacement multi-wave seismic wave field in thin-layered medium on its basis.
According to an aspect of the invention, it is proposed that a kind of analogy method of thin layer displacement multi-wave seismic wave field, may include with
Lower step:
1) frequency domain thin layer is obtained based on thin layer displacement-reflection transmission matrix equation and is displaced PP wave reflection coefficients RppWith PS wave
Reflection Rps；
2) FFT transform is carried out to seismic wavelet w (t), obtained frequency domain seismic wavelet w (f)；
3) PP wave reflection coefficients R is displaced based on the frequency domain thin layer that step 1) obtainsppWith PS wave reflection coefficients RpsAnd step
The rapid frequency domain seismic wavelet w (f) 2) obtained obtains frequency domain PP and involves PS wave reflection seismic wave field；
4) frequency domain PP obtained to step 3) involves PS wave reflection seismic wave field and carries out obtaining the time against FFT transform
Domain PP involves PS wave reflection seismic wave field.
Preferably, the thin layer displacement-reflection transmission matrix equation are as follows:
In formula (1), RppPP wave reflection coefficient, R are displaced for frequency domain thin layerpsIt is anti-that PS wave is displaced for frequency domain thin layer
Penetrate coefficient, TppPP wave transmission coefficient, T are displaced for frequency domain thin layerpsIt is displaced PP wave transmission coefficient for frequency domain thin layer,It is defined as follows:
In formula (2)-(6),
a11=a44=2sin2β2cos P+cos2β2Cos Q,
a12=a34=-j (tan α2cos2β2sin P-sin2β2Sin Q),
a22=a33=cos2 β2cos P+2sin2β2Cos Q,
a31=a42=2j ρ2ωvS2sinβ2cos2β2(cos Q-cos P),
Wherein, α1、α2、α3The wave ray of longitudinal wave and the angle of normal, β respectively in the 1st, 2,3 layer1、β2、β3Respectively
The wave ray of shear wave and the angle of normal in the 1st, 2,3 layer, vPi,vSi,ρiIt (i=1,2,3) is respectively the 1st, 2,3 layer of longitudinal wave
Speed, shear wave velocity and density, the 1st layer of overlying strata for thin layer, the 2nd layer is thin layer, and the 3rd layer is stratum under thin layer,H is thickness of thin layer, and the π of ω=2 f, f are the frequency of incidence wave,
According to another aspect of the invention, it is proposed that a kind of simulator of thin layer displacement multi-wave seismic wave field.Described device
Including memory, processor and store the computer program that can be run on a memory and on a processor, which is characterized in that
The processor performs the steps of when executing described program
1) frequency domain thin layer is obtained based on thin layer displacement-reflection transmission matrix equation and is displaced PP wave reflection coefficients RppWith PS wave
Reflection Rps；
2) FFT transform is carried out to seismic wavelet w (t), obtained frequency domain seismic wavelet w (f)；
3) PP wave reflection coefficients R is displaced based on the frequency domain thin layer that step 1) obtainsppWith PS wave reflection coefficients RpsAnd step
The rapid frequency domain seismic wavelet w (f) 2) obtained obtains frequency domain PP and involves PS wave reflection seismic wave field；
4) frequency domain PP obtained to step 3) involves PS wave reflection seismic wave field and carries out obtaining the time against FFT transform
Domain PP involves PS wave reflection seismic wave field.
Preferably, the thin layer displacement-reflection transmission matrix equation are as follows:
In formula (1), RppPP wave reflection coefficient, R are displaced for frequency domain thin layerpsIt is anti-that PS wave is displaced for frequency domain thin layer
Penetrate coefficient, TppPP wave transmission coefficient, T are displaced for frequency domain thin layerpsIt is displaced PP wave transmission coefficient for frequency domain thin layer,It is defined as follows:
In formula (2)-(6),
a11=a44=2sin2β2cos P+cos2β2Cos Q,
a12=a34=-j (tan α2cos2β2sin P-sin2β2Sin Q),
a22=a33=cos2 β2cos P+2sin2β2Cos Q,
a31=a42=2j ρ2ωvs2sinβ2cos2β2(cos Q-cos P),
Wherein, α1、α2、α3The wave ray of longitudinal wave and the angle of normal, β respectively in the 1st, 2,3 layer1、β2, β 3 be respectively
The wave ray of shear wave and the angle of normal in the 1st, 2,3 layer, vPi, vSi, ρi(i=1,2,3) is respectively the 1st, 2,3 layer of longitudinal wave
Speed, shear wave velocity and density, the 1st layer of overlying strata for thin layer, the 2nd layer is thin layer, and the 3rd layer is stratum under thin layer,H is thickness of thin layer, and the π of ω=2 f, f are the frequency of incidence wave,
The analogy method and system of thin layer displacement multi-wave seismic wave field according to the present invention precisely simulation can be based on position
The thin-layered medium reflected P P for moving wave function involves PS wave, and reflected wave field includes that conversion involves interbed multiple information, Ke Yiquan
Face, accurately reflection thin layer act on the transformation of seismic wavelet.
Methods and apparatus of the present invention has other characteristics and advantages, these characteristics and advantages are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing and subsequent specific implementation being incorporated herein
It is stated in detail in example, these the drawings and specific embodiments are used together to explain specific principle of the invention.
Detailed description of the invention
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other
Purpose, feature and advantage will be apparent, wherein in exemplary embodiments of the present invention, identical reference label is usual
Represent same parts.
Fig. 1 is the process that the thin layer of exemplary implementation scheme according to the present invention is displaced the analogy method of multi-wave seismic wave field
Figure；
Reflection and transmission schematic diagram of the Fig. 2 for thin layer in the case of longitudinal wave incidence wave；
Fig. 3 a-3f is the PP wave synthetic seismic record in coal seam under different-thickness；
Fig. 4 a-4f is the PS wave synthetic seismic record in coal seam under different-thickness；
Fig. 5 a shows that the rolling land the single interface PP shake reflected wave field of the method for the present invention simulation, Fig. 5 b show conventional method simulation
The rolling land single interface PP shakes reflected wave field, and Fig. 5 c shows that the rolling land the single interface PS shake reflected wave field of the method for the present invention simulation, Fig. 5 d are aobvious
Show that the rolling land the single interface PS shake reflected wave field of conventional method simulation, Fig. 5 e show that the PP wave seismic wave field of two methods simulation is residual
Difference；Fig. 5 f shows the PS wave seismic wave field residual error of two methods simulation.
Specific embodiment
The present invention will be described in more detail below with reference to accompanying drawings.Although showing the preferred embodiment of the present invention in attached drawing,
However, it is to be appreciated that may be realized in various forms the present invention and should not be limited by the embodiments set forth herein.On the contrary, providing
These embodiments are of the invention more thorough and complete in order to make, and can will fully convey the scope of the invention to ability
The technical staff in domain.
The analogy method of the thin layer displacement multi-wave seismic wave field of exemplary implementation scheme according to the present invention described in detail below,
It is specifically included that
Step 1: frequency domain thin layer being obtained based on thin layer displacement-reflection transmission matrix equation and is displaced PP wave reflection coefficients RppWith
PS wave reflection coefficients Rps。
With reference to Fig. 1, PP wave reflection coefficients R is obtained in step 1ppWith PS wave reflection coefficients RpsSpecific steps are as follows: input
Thin layer elastic parameter, thickness and frequency, are substituted into thin layer displacement-reflection matrix equation, then can solve PP wave reflection system
Number RppWith PS wave reflection coefficients RpsValue.
Fig. 2 is reflection and the transmission schematic diagram of thin layer under longitudinal wave condition of incidence.For thin layer model as shown in Figure 2,
Shown in the accurate reflection and transmission matrix equation of thin layer such as formula (1) based on displacement function:
In formula (1), RppPP wave reflection coefficient, R are displaced for frequency domain thin layerpsIt is anti-that PS wave is displaced for frequency domain thin layer
Penetrate coefficient, TppPP wave transmission coefficient, T are displaced for frequency domain thin layerpsPP wave transmission coefficient is displaced for frequency domain thin layer.
It is defined as follows:
In formula (2)-(6),
a11=a44=2sin2β2cos P+cos2β2Cos Q,
a12=a34=-j (tan α2cos2β2sin P-sin2β2Sin Q),
a22=a33=cos2 β2cos P+2sin2β2Cos Q,
a31=a42=2j ρ2ωvS2sinβ2cos2β2(cos Q-cos P)
Wherein, as shown in Fig. 2, α1、α2、α3The wave ray of longitudinal wave and the angle of normal, β respectively in the 1st, 2,3 layer1、
β2、β3Respectively the wave ray of shear wave and the angle of normal, h are thickness of thin layer in the 1st, 2,3 layer.vPi, vSi, ρi(i=1,2,
It 3) is respectively the 1st, 2,3 layer of velocity of longitudinal wave, shear wave velocity and density, the 1st layer of overlying strata for thin layer, the 2nd layer is thin layer, the
3 layers are stratum under thin layer, and ω is circular frequency, and the π of ω=2 f, f are the frequency of incidence wave,P is longitudinal wave vertical wavenumber
With the product of thickness of thin layer, Q is the product of shear wave vertical wavenumber and thickness of thin layer, specifically,
When taking thickness of thin layer is 0, it is Zoeppritz equation that formula (1), which is degenerated, as shown in formula (7).
Step 2: FFT transform being carried out to seismic wavelet w (t), is obtained frequency domain seismic wavelet w (f).
Step 3: PP wave reflection coefficients R is displaced based on the frequency domain thin layer that step 1) obtainsppWith PS wave reflection coefficients RpsWith
And the frequency domain seismic wavelet w (f) that step 2) obtains, it obtains frequency domain PP and involves PS wave reflection seismic wave field.
In one example, by by frequency domain seismic wavelet w (f) and PP wave reflection coefficients RppWith PS wave reflection coefficient
RpsDot product is carried out, respectively to obtain frequency domain PP wave and PS wave reflection seismic wave field.
Step 4: frequency domain PP obtained to step 3) involves PS wave reflection seismic wave field and carries out inverse FFT transform, obtains
Time-domain PP involves PS wave reflection seismic wave field.
It, can be by PP wave reflection coefficients R since thin bed reflection coefficient is to rely on the function of frequencyppIt is denoted as Rpp
(f), by PS wave reflection coefficients RspIt is denoted as Rps(f)。
Time-domain PP wave synthetic seismic record are as follows:
Spp(t)=FFT-1[Rpp(f)·w(f)] (8)
Time-domain PS wave synthetic seismic record are as follows:
Sps(t)=FFT-1[Rps(f)·w(f)] (8)
By above step, then obtains time-domain PP and involve PS wave reflection seismic wave field.
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Bog-head coal layer model as shown in Table 1 is established, and utilizes the mould of thin layer according to the present invention displacement multi-wave seismic wave field
Quasi- method realizes thin layer multi-wave seismic wave-field simulation.Since thin response is that a top bottom reflection, interbed multiple, converted wave are folded
Complex wave made of adding, the aliasing effect of wave field makes can not to involve interlayer more for Direct Recognition conversion on the thin layer wave field of simulation
Subwave, therefore the thickness by modifying intermediate coal seam gradually becomes thin layer by thick-layer, analyzes analogy method of the present invention to layer transition
Involve reporting situations for multiple wave.
1 coal seam model parameter of table
Coal seam thickness is modified, the displacement earthquake PP in the case of simulation different-thickness involves PS wave field.Fig. 3 a, 3b, 3c,
3d, 3e, 3f are respectively PP wave synthetic seismic record of the thickness in 2 λ, λ, λ/2, λ/4, λ/6, λ/12, Fig. 4 a, 4b, 4c,
4d, 4e, 4f are respectively PS wave synthetic seismic record of the thickness in 2 λ, λ, λ/2, λ/4, λ/6, λ/12.As seen from the figure, exist
In the biggish situation of thickness, coal seam top, bottom reflection can be clearlyed distinguish, conversion, which involves interbed multiple, in waveform to know
Not, this is because it is then different larger with coal seam top, bottom reflection arrival time difference.And with the reduction of thickness, thin layer top, bottom are reflected,
Involve multiple wave including conversion, waveform moves closer to, and forms interference, is in aliasing.In the lesser situation of coal seam thickness, nothing
Method identifies thin layer top, bottom interface.
In order to further verify effectiveness of the invention, the thickness of thin layer in table 1 is set as 0 meter, formed Sandy Silt with
On the one hand PP wave and PS in the case of above-mentioned single interface special case are simulated in mud stone list impedance contrast interface using analogy method of the invention
Wave seismic wave field, as shown in Fig. 5 (a), 5 (c)；On the other hand skill is simulated using existing single impedance contrast interface echo field
Art --- the PP wave at the above-mentioned single interface of time-domain convolution simulated implementation and the Simulation of Seismic Wave of PS wave, such as Fig. 5 (b), 5 (d) institutes
Show.Analog result of the invention and conventional single interface are fitted seismic wave field residual error such as Fig. 5 (e), 5 (f).It can from analog result comparison
Both see it is completely the same, to demonstrate the system to the correctness and validity of single impedance contrast interface Simulation of Seismic Wave,
From correctness, the universality of side illustration the method for the present invention simulation wave field.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.The selection of term used herein, purport
In the principle, practical application or improvement to the technology in market for best explaining each embodiment, or make the art
Other those of ordinary skill can understand each embodiment disclosed herein.
Claims (6)
1. a kind of analogy method of thin layer displacement multi-wave seismic wave field, which comprises the following steps:
1) frequency domain thin layer is obtained based on thin layer displacement-reflection transmission matrix equation and is displaced PP wave reflection coefficients RppWith PS wave reflection
Coefficients Rps；
2) FFT transform is carried out to seismic wavelet w (t), obtained frequency domain seismic wavelet w (f)；
3) PP wave reflection coefficients R is displaced based on the frequency domain thin layer that step 1) obtainsppWith PS wave reflection coefficients RpsAnd step 2)
The frequency domain seismic wavelet w (f) of acquisition obtains frequency domain PP and involves PS wave reflection seismic wave field；
4) frequency domain PP obtained to step 3) involves PS wave reflection seismic wave field and carries out obtaining time-domain PP against FFT transform
Involve PS wave reflection seismic wave field.
2. the analogy method of thin layer displacement multi-wave seismic wave field according to claim 1, which is characterized in that the thin layer position
Move reflection and transmission matrix equation are as follows:
In formula (1), RppPP wave reflection coefficient, R are displaced for frequency domain thin layerpsPS wave reflection system is displaced for frequency domain thin layer
Number, TppPP wave transmission coefficient, T are displaced for frequency domain thin layerpsIt is displaced PP wave transmission coefficient for frequency domain thin layer,It is defined as follows:
In formula (2)-(6),
a11=a44=2sin2β2cos P+cos2β2cos Q
a12=a34=-j (tan α2cos 2β2sin P-sin 2β2Sin Q),
a22=a332 β of=cos2cos P+2sin2β2cos Q
a31=a42=2j ρ2ωvS2sinβ2cos 2β2(cos Q-cos P),
Wherein, α1、α2、α3The wave ray of longitudinal wave and the angle of normal, β respectively in the 1st, 2,3 layer1、β2、β3Respectively
1, the wave ray of shear wave and the angle of normal in 2,3 layers, vPi,vSi,ρiIt (i=1,2,3) is respectively the 1st, 2,3 layer of longitudinal wave speed
Degree, shear wave velocity and density, the 1st layer of overlying strata for thin layer, the 2nd layer is thin layer, and the 3rd layer is stratum under thin layer,H is thickness of thin layer, and the π of ω=2 f, f are the frequency of incidence wave,
3. the analogy method of thin layer displacement multi-wave seismic wave field according to claim 1, which is characterized in that in step 3)
In, by by frequency domain seismic wavelet w (f) and PP wave reflection coefficients RppWith PS wave reflection coefficients RpsDot product is carried out, respectively to obtain
Obtain frequency domain PP wave and PS wave reflection seismic wave field.
4. a kind of simulator of thin layer displacement multi-wave seismic wave field, which is characterized in that described device includes memory, processor
And store the computer program that can be run on a memory and on a processor, which is characterized in that the processor executes institute
It is performed the steps of when stating program
1) frequency domain thin layer is obtained based on thin layer displacement-reflection transmission matrix equation and is displaced PP wave reflection coefficients RppWith PS wave reflection
Coefficients Rps；
2) FFT transform is carried out to seismic wavelet w (t), obtained frequency domain seismic wavelet w (f)；
3) PP wave reflection coefficients R is displaced based on the frequency domain thin layer that step 1) obtainsppWith PS wave reflection coefficients RpsAnd step 2)
The frequency domain seismic wavelet w (f) of acquisition obtains frequency domain PP and involves PS wave reflection seismic wave field；
4) frequency domain PP obtained to step 3) involves PS wave reflection seismic wave field and carries out obtaining time-domain PP against FFT transform
Involve PS wave reflection seismic wave field.
5. the simulator of thin layer displacement multi-wave seismic wave field according to claim 4, which is characterized in that the thin layer position
Move reflection and transmission matrix equation are as follows:
In formula (1), RppPP wave reflection coefficient, R are displaced for frequency domain thin layerpsPS wave reflection system is displaced for frequency domain thin layer
Number, TppPP wave transmission coefficient, T are displaced for frequency domain thin layerpsIt is displaced PP wave transmission coefficient for frequency domain thin layer,It is defined as follows:
In formula (2)-(6),
a11=a44=2sin2β2cos P+cos 2β2cos Q
a12=a34=-j (tan α2cos 2β2sin P-sin2β2Sin Q),
a22=a332 β of=cos2cos P+2sin2β2cos Q
a31=a42=2j ρ2ωvS2sinβ2cos 2β2(cos Q-cos P),
Wherein, α1、α2、α3The wave ray of longitudinal wave and the angle of normal, β respectively in the 1st, 2,3 layer1、β2、β3Respectively
1, the wave ray of shear wave and the angle of normal in 2,3 layers, vPi, vSi, ρi(i=1,2,3) is respectively the 1st, 2,3 layer of longitudinal wave speed
Degree, shear wave velocity and density, the 1st layer of overlying strata for thin layer, the 2nd layer is thin layer, and the 3rd layer is stratum under thin layer,H is thickness of thin layer, and the π of ω=2 f, f are the frequency of incidence wave,
6. the simulator of thin layer displacement multi-wave seismic wave field according to claim 4, which is characterized in that in step 3)
In, by by frequency domain seismic wavelet w (f) and PP wave reflection coefficients RppWith PS wave reflection coefficients RpsDot product is carried out, respectively to obtain
Obtain frequency domain PP wave and PS wave reflection seismic wave field.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710647906.6A CN109324343A (en) | 2017-08-01 | 2017-08-01 | A kind of analogy method and system of thin layer displacement multi-wave seismic wave field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710647906.6A CN109324343A (en) | 2017-08-01 | 2017-08-01 | A kind of analogy method and system of thin layer displacement multi-wave seismic wave field |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109324343A true CN109324343A (en) | 2019-02-12 |
Family
ID=65246042
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710647906.6A Pending CN109324343A (en) | 2017-08-01 | 2017-08-01 | A kind of analogy method and system of thin layer displacement multi-wave seismic wave field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109324343A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110441817A (en) * | 2019-06-27 | 2019-11-12 | 西北大学 | The method changed in pore media based on 4D PP-PS wave inverting CCUS pressure and saturation degree |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104950332A (en) * | 2015-06-18 | 2015-09-30 | 河海大学 | Method for calculating plane wave reflection coefficients in elastic multi-layered medium |
| CN105629301A (en) * | 2016-03-29 | 2016-06-01 | 中国地质大学（北京） | Thin layer elastic wave reflection coefficient fast solving method |
| CN106772578A (en) * | 2016-12-07 | 2017-05-31 | 中国矿业大学(北京) | A kind of method and apparatus of synthetic seismogram |
-
2017
- 2017-08-01 CN CN201710647906.6A patent/CN109324343A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104950332A (en) * | 2015-06-18 | 2015-09-30 | 河海大学 | Method for calculating plane wave reflection coefficients in elastic multi-layered medium |
| CN105629301A (en) * | 2016-03-29 | 2016-06-01 | 中国地质大学（北京） | Thin layer elastic wave reflection coefficient fast solving method |
| CN106772578A (en) * | 2016-12-07 | 2017-05-31 | 中国矿业大学(北京) | A kind of method and apparatus of synthetic seismogram |
Non-Patent Citations (1)
| Title |
|---|
| 王书伦: "薄互层合成地震记录制作方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110441817A (en) * | 2019-06-27 | 2019-11-12 | 西北大学 | The method changed in pore media based on 4D PP-PS wave inverting CCUS pressure and saturation degree |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Amplitude and AVO responses of a single thin bed | |
| Maraschini et al. | A new misfit function for multimodal inversion of surface waves | |
| Boaga et al. | Mode misidentification in Rayleigh waves: Ellipticity as a cause and a cure | |
| CN101329405B (en) | Simple method of multi-parameter seismic inversion | |
| Wang et al. | Current developments on micro-seismic data processing | |
| CN102736107B (en) | Energy constraint heterogeneous reservoir thickness identification system | |
| Cheng et al. | Automated data selection in the Tau–p domain: Application to passive surface wave imaging | |
| CN104570072A (en) | Method for modeling reflection coefficient of spherical PP wave in viscoelastic medium | |
| CN105652316A (en) | Intelligent optimized earthquake multi-attribute fusion method based on crack model | |
| WO2008024150A2 (en) | Converted mode seismic survey design | |
| Masoni et al. | Alternative misfit functions for FWI applied to surface waves | |
| CN101359056A (en) | Method for generating longitudinal wave time domain high precision converted wave sectional plane | |
| Ning et al. | High-frequency surface-wave imaging from traffic-induced noise by selecting in-line sources | |
| Pan et al. | AVO/AVF analysis of thin beds in elastic media | |
| Jiang et al. | Pronounced seismic anisotropy in Kanto sedimentary basin: A case study of using dense arrays, ambient noise seismology, and multi‐modal surface‐wave imaging | |
| CN105092343A (en) | Method for eliminating thin layer tuning effect based on prestack gather | |
| CN109324343A (en) | A kind of analogy method and system of thin layer displacement multi-wave seismic wave field | |
| Ursenbach et al. | Efficient spherical-wave AVO modeling | |
| Tadeu et al. | Wave motion between two fluid-filled boreholes in an elastic medium | |
| Ávila-Carrera et al. | Simulating elastic wave propagation in boreholes: Fundamentals of seismic response and quantitative interpretation of well log data | |
| Li et al. | Elastic transmitted wave reverse time migration for imaging Earth’s interior discontinuities: a numerical study | |
| Hayashi | Near-surface engineering | |
| Farfour et al. | A review about recent seismic techniques in shale-gas exploration | |
| Shi et al. | Data processing workflow for reflections from cement-formation interface in ultrasonic pitch-catch measurement in cased hole | |
| Sinha et al. | Hydraulic fracture characterization using cased-hole sonic data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190212 |