CN101082676A - Earthquake post-stack forward method of undulating surface - Google Patents

Earthquake post-stack forward method of undulating surface Download PDF

Info

Publication number
CN101082676A
CN101082676A CNA200710049482XA CN200710049482A CN101082676A CN 101082676 A CN101082676 A CN 101082676A CN A200710049482X A CNA200710049482X A CN A200710049482XA CN 200710049482 A CN200710049482 A CN 200710049482A CN 101082676 A CN101082676 A CN 101082676A
Authority
CN
China
Prior art keywords
wave field
relief surface
field
continuation
wave
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.)
Granted
Application number
CNA200710049482XA
Other languages
Chinese (zh)
Other versions
CN100545679C (en
Inventor
熊晓军
贺振华
黄德济
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CNB200710049482XA priority Critical patent/CN100545679C/en
Publication of CN101082676A publication Critical patent/CN101082676A/en
Application granted granted Critical
Publication of CN100545679C publication Critical patent/CN100545679C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a post-stack forward modeling method of fluctuating ground, which comprises the following steps: proceeding rectangular packing for geological mould of the fluctuating ground; filling the fluctuating ground at random speed value or density value; adopting Snail elementary to calculate the reflective factor model; eliminating the reflective influence of fluctuating ground to transform to frerquency area as original wave field; utilizing Fuliye limit differential continuation operator of single-pass fluctuating equation to continue the wave field along depth axle from upwards to downwards; keeping the wave field between input datum face (fluctuating ground) and output datum surface at zero; transforming the wave field of frequency area after continuing to the time domain to image; processing the imaging result into earthquake section image through normal display software; possessing good earthquake wave field analoging effect for geological model of the fluctuating ground.

Description

The earthquake poststack forward modeling method of relief surface
Technical field
The present invention relates to the positive artistic skills art of the earthquake poststack category in the reflection wave seismic wave field numerical simulation process, is a kind of poststack forward modeling method of relief surface seismic wave field numerical simulation.
Background technology
The numerical simulation of seismic wave field is exactly under the situation of the structural model of the known underground medium of supposition and respective physical parameter, the propagation law of simulated earthquake ripple in underground various media, and calculate on ground or underground each observation station a kind of earthquake simulation method of the numerical value seismologic record that should observe.It on the one hand can be for earthquake data acquisition, processing, explanation provide theoretical foundation, the science of appraisal procedure, feasibility and advance; Can be used for checking the credibility of various interpretation results on the other hand, and the correctness of various inversion algorithms and reliability.From the angle of seismic wave propagation, the positive artistic skills art of the earthquake poststack of routine can be divided into ray-tracing scheme and wave equation method.Ray-tracing scheme can be handled relief surface neatly, but it is based on the kinematics of seismic event, the seismic event wave theory is reduced to ray theory, the main kinematics character of considering seismic wave propagation, lack the dynamic information of seismic event, result of calculation is similar to a certain extent, and using has the certain limitation condition, low to the complex geological tectonic imaging precision, the blind area appears sometimes.The wave equation method is directly found the solution wave equation by numerical method, the dynamic characteristic that has comprised seismic wave propagation, wave field abundant information, analog result are comparatively accurate, but the situation that it can only the processing horizontal face of land, the reference field of wave field extrapolation is necessary for a horizontal interface.
Summary of the invention
The present invention will provide a kind of earthquake poststack forward modeling method that is applicable to relief surface based on one-way wave equation, it not only can adapt to the surface conditions of fluctuating, and adopt one-way wave equation to calculate, have the imaging precision of very fast computing velocity and very high complex structure.
The positive artistic skills art of earthquake poststack of the relief surface based on one-way wave equation of the present invention, concrete steps comprise:
(1) adopt rectangle frame that the target geologic model is carried out overall package, the end face of rectangle frame is the reference field of output, and the relief surface top adopts arbitrary speed value or density value to carry out filling, and adopts regular grid that rate pattern and density model are carried out gridding.
(2) adopt the Si Naier theorem to calculate reflection coefficient, and utilize Ricker wavelet to carry out convolution and handle, obtain reflectivity model.
(3) eliminate the reflections affect of relief surface, and transform to frequency field as initial continuation wave field.
(4) utilize the Fourier finite difference continuation operator of one-way wave equation to carry out wave field extrapolation from the bottom up, and the wave field between maintenance input reference face (relief surface) and the output reference face is zero along degree of depth axle.
(5) wave field transformation with the frequency field after the continuation carries out imaging to time domain.
(6) software for display by routine is treated to seismic section image with imaging results.
The earthquake poststack forward modeling method of the relief surface based on one-way wave equation of the present invention, to keep a numerical value between input reference face and output reference face be that zero wave field layer solves the relief surface problem of wave equation earthquake poststack in just drilling, and drawn the advantage that ray-tracing scheme is handled relief surface flexibly.
The earthquake poststack forward modeling method of the relief surface based on one-way wave equation of the present invention, utilizing Fourier finite difference continuation operator to carry out wave field extrapolation in frequency field calculates, ability with stronger processing lateral velocity variation, have very high imaging precision for complex structure, improved ray-tracing scheme for the low shortcoming of complex structure imaging precision.
Specific implementation principle of the present invention is as follows:
The end face of selected rectangle frame is the reference field of wave field extrapolation, and it is a horizontal interface, and is positioned on the peak of relief surface; The right margin of rectangle frame, left margin and bottom boundaries overlap fully with the geologic model of input.Geologic model after the encapsulation is carried out the mesh generation of rule, so the surface configuration that rises and falls becomes the part of regular grid after by discretize, the grid on the relief surface adopts arbitrarily that velocity amplitude and density value carry out filling.Utilize the Si Naier theorem to calculate reflection coefficient,
R ( x , y , z i , t ) = Den ( x , y , z i + 1 , t ) × Vel ( x , y , z i + 1 , t ) - Den ( x , y , z i , t ) × Vel ( x , y , z i , t ) Den ( x , y , z i + 1 , t ) × Vel ( x , y , z i + 1 , t ) + Den ( x , y , z i , t ) × Vel ( x , y , z i , t )
And adopt Ricker wavelet (dominant frequency is selected according to the stratum characteristic in target work area, generally between 10Hz～90Hz) to carry out convolution and handle.
Eliminate the reflections affect of relief surface, and adopt Fourier transform with it transform to frequency field R (x, y, z, ω), as the wave field of initial continuation.Given frequencies omega is extracted initial wave field W (x, y, the z at this frequency place, ω), adopt Fourier finite difference continuation operator to carry out wave field extrapolation from the bottom up with up-going wave equation along degree of depth axle, the wave field after the continuation is W ' (x, y, 0, ω), this process can be expressed as
W′(x，y，0，ω)＝FFD z[W(x，y，z，ω)]
X wherein, y is respectively the coordinate of two horizontal directions in the rectangular coordinate system in space, and z is the coordinate of depth direction, and ω is a circular frequency, and FFDz is a Fourier finite difference continuation operator, (z t) is the density value of gridding to Den for x, y, (z t) is the velocity amplitude of gridding to Vel for x, y.The input wave field of next depth z i is a last depth z in the continuation process I+1The continuation result, one zero wave field layer is set between output reference face and input reference face simultaneously participates in normal wave field and propagate and the continuation process.
All frequencies are carried out cycle calculations, obtain on the output reference face frequency field wave field U (x, y, 0, ω).To U (x, y, 0, ω) carry out inverse-Fourier transform and transform to the time domain imaging, (t), the software for display by routine is treated to the stratigraphic section image with imaging results to simulation wave field I that must be to the end then for x, y.
Above-mentioned continuation operator obtains according to the one-way wave equation close approximation, adopts the stable Fourier finite difference operator of frequency field to realize (Ristow, 1994).This continuation operator can be approximately,
ω 2 v 2 + ∂ 2 ∂ x 2 = A 1 + A 2 + A 3
In the formula,
A 1 = ω 2 c 2 + ∂ 2 ∂ x 2 ; A 2 = ( ω v - ω c ) ; A 3 = ω v [ 1 - c v ] a v 2 ω 2 ∂ 2 ∂ x 2 1 + b v 2 ω 2 ∂ 2 ∂ x 2
Wherein v is the actual speed field, and c is the reference velocity field, and ω is an angular frequency, parameter a, and the value of b has influence on the inclination maximum of wave field imaging:
Figure A20071004948200062
The present invention adopts the seismic wave field poststack method for numerical simulation that is adapted to relief surface based on one-way wave equation, has following characteristics, mainly shows as:
(1) can handle the relief surface condition of earthquake poststack in just drilling neatly.
(2) the Fourier finite difference continuation operator with relief surface, one-way wave equation and frequency field organically combines.
(3) has accurate wave field imaging effect for the complex structure under the relief surface condition.
Description of drawings
Fig. 1 is the two-dimentional relief surface model of an acline, and the maximum fluctuating discrepancy in elevation on the face of land is 380 meters.
Fig. 2 carries out result behind the earthquake poststack forward simulation with present technique to the rate pattern among Fig. 1, and the reflection line-ups of acline is no longer straight as can be seen, but opposite with surface configuration; In addition, simulation only comprises primary reflection, signal to noise ratio (S/N ratio) height.
Fig. 3 is the two-dimentional relief surface model of a complex geological structure, the low-velocity anomal body that it is 2000 meter per seconds that the near surface place comprises two speed.
Fig. 4 carries out result behind the earthquake poststack forward simulation with present technique to the rate pattern among Fig. 3, and basically identical reflection interval of two of near surface low-velocity anomal bodies is not subjected to the influence of rugged topography as can be seen; In addition, because the influence of rugged topography and low-velocity anomal body, serious distortion has taken place in the reflection line-ups that is positioned at the acline at model middle part.
Embodiment
Earthquake poststack forward modeling method based on the relief surface of one-way wave equation is specially following steps:
(1) adopt rectangle frame that the target geologic model is carried out overall package, the end face of rectangle frame is the reference field of output, and the relief surface top adopts arbitrary speed value or density value to carry out filling, and adopts regular grid that rate pattern and density model are carried out gridding.If have only rate pattern, then adopt the experimental formula bulk density model of density and speed.
(2) adopt the Si Naier theorem to calculate reflection coefficient, and utilize Ricker wavelet to carry out convolution and handle, obtain reflectivity model.The dominant frequency of Ricker wavelet is selected (generally between 10Hz～90Hz) according to the stratum characteristic in work area.
(3) eliminate the reflections affect of relief surface, and transform to frequency field as initial continuation wave field.
(4) utilize the Fourier finite difference continuation operator of one-way wave equation to carry out wave field extrapolation from the bottom up, and the wave field between maintenance input reference face (relief surface) and the output reference face is zero along degree of depth axle.
Specifically: the wave field extrapolation of carrying out in the frequency circulation is calculated, and, with reference to primary wave with stop time of ripple, determine the frequency range f1 and the f2 of wave field extrapolation according to the stratum characteristic in work area, generally between 1Hz～5Hz, f2 is generally between 50Hz～100Hz for f1; In the single frequency circulation, carry out wave field extrapolation from the bottom up along degree of depth axle and calculate, determine Fourier finite difference continuation operator ω 2 v 2 + ∂ 2 ∂ x 2 = A 1 + A 2 + A 3 Wherein:
A 1 = ω 2 c 2 + ∂ 2 ∂ x 2 ; A 2 = ( ω v - ω c ) ; A 3 = ω v [ 1 - c v ] a v 2 ω 2 ∂ 2 ∂ x 2 1 + b v 2 ω 2 ∂ 2 ∂ x 2
Wherein v is the actual speed field, and c is the reference velocity field, and ω is an angular frequency, and selects parameter a according to the inclined degree on stratum, the value of b, and parameter a, the corresponding relation of the value of b and the inclination maximum of migration imaging is:
Figure A20071004948200075
Figure A20071004948200076
Figure A20071004948200077
Next depth z in the continuation process iThe input wave field be a last depth z I+1The continuation result, one zero wave field layer is set between output reference face and input reference face simultaneously participates in normal wave field and propagate and the continuation process.Wave field value between intact f1 of cycle calculations and the f2 is to the processing that zeroizes of other wave field value.
(5) wave field transformation with the frequency field after the continuation carries out imaging to time domain.
(6) software for display by routine is treated to seismic section image with imaging results.
Fig. 1 and Fig. 2 are the test examples of the two-dimentional relief surface model of an acline.
(1) rate pattern to the rectangle encapsulation carries out gridding, and the velocity setting of relief surface top is 500 meter per seconds (selected arbitrarily values), and density model is a constant 1; Determine that according to track pitch the spacing of grid directions X is 20 meters, determine that according to the degree of depth spacing of grid Z direction is 5 meters.
(2) adopt the Si Naier theorem to calculate reflection coefficient, and utilize dominant frequency to handle, obtain reflectivity model for the Ricker wavelet of 50Hz carries out convolution.
(3) eliminate the reflections affect of relief surface, and transform to frequency field as initial continuation wave field.
(4) utilize the Fourier finite difference continuation operator of one-way wave equation to carry out wave field extrapolation from the bottom up, and the wave field between maintenance input reference face (relief surface) and the output reference face is zero along degree of depth axle.Specifically,, determine the frequency range f1 and the f2 of cycle calculations according to the stratum characteristic in work area with reference to primary wave with stop time of ripple, general it be set to the effective band of geological data, f1=1 here, f2=100.Because the stratum is an acline, parameter a, b, the value inclination angle is the situation of 45 degree, is a=0.5, b=0.25.The initial wave field of each frequency between f1 and the f2 is carried out wave field extrapolation from the bottom up along degree of depth axle calculate, its continuation operator is the Fourier finite difference operator of the stable frequency field narrated previously.Next depth z in the continuation process iThe input wave field be a last depth z I+1The continuation result, one zero wave field layer is set between output reference face and input reference face simultaneously participates in normal wave field and propagate and the continuation process.Wave field value between intact f1 of cycle calculations and the f2 is to the processing that zeroizes of other wave field value.
(5) wave field transformation with the frequency field after the continuation carries out imaging to time domain.
(6) software for display by routine is treated to seismic section image with imaging results.Fig. 2 is the seismic section of earthquake poststack simulation.
Fig. 3 and Fig. 4 are the test examples of the two-dimentional relief surface model of an acline.
(1) rate pattern to the rectangle encapsulation carries out gridding, and the velocity setting of relief surface top is 300 meter per seconds (selected arbitrarily values), and density model is a constant 1; Determine that according to track pitch the spacing of grid directions X is 20 meters, determine that according to the degree of depth spacing of grid Z direction is 5 meters.
(2) adopt the Si Naier theorem to calculate reflection coefficient, and utilize dominant frequency to handle, obtain reflectivity model for the Ricker wavelet of 50Hz carries out convolution.
(3) eliminate the reflections affect of relief surface, and transform to frequency field as initial continuation wave field.
(4) utilize the Fourier finite difference continuation operator of one-way wave equation to carry out wave field extrapolation from the bottom up, and the wave field between maintenance input reference face (relief surface) and the output reference face is zero along degree of depth axle.Specifically,, determine the frequency range f1 and the f2 of cycle calculations according to the stratum characteristic in work area with reference to primary wave with stop time of ripple, general it be set to the effective band of geological data, f1=1 here, f2=100.Because the stratum is an acline, parameter a, b, the value inclination angle is the situation of 45 degree, is a=0.5, b=0.25.The initial wave field of each frequency between f1 and the f2 is carried out wave field extrapolation from the bottom up along degree of depth axle calculate, its continuation operator is the Fourier finite difference operator of the stable frequency field narrated previously.The input wave field of next depth z i is a last depth z in the continuation process I+1The continuation result, one zero wave field layer is set between output reference face and input reference face simultaneously participates in normal wave field and propagate and the continuation process.Wave field value between intact f1 of cycle calculations and the f2 is to the processing that zeroizes of other wave field value.
(5) wave field transformation with the frequency field after the continuation carries out imaging to time domain.
(6) software for display by routine is treated to seismic section image with imaging results.Fig. 4 is the seismic section of earthquake poststack simulation.

Claims (4)

1, the earthquake poststack forward modeling method of relief surface is characterized in that adopting following steps:
(1) adopt rectangle frame that the target geologic model is carried out overall package, the end face of rectangle frame is the reference field of output, and the relief surface top adopts arbitrary speed value or density value to carry out filling, and adopts regular grid that rate pattern and density model are carried out gridding.
(2) adopt the Si Naier theorem to calculate reflection coefficient, and utilize Ricker wavelet to carry out convolution and handle, obtain reflectivity model.
(3) eliminate the reflections affect of relief surface, and transform to frequency field as initial continuation wave field.
(4) utilize the Fourier finite difference continuation operator of one-way wave equation to carry out wave field extrapolation from the bottom up, and the wave field between maintenance input reference face (relief surface) and the output reference face (rectangular top) is zero along degree of depth axle.
(5) wave field transformation with the frequency field after the continuation carries out imaging to time domain.
(6) software for display by routine is treated to seismic section image with imaging results.
2, the earthquake poststack forward modeling method of relief surface according to claim 1, it is characterized in that: object module is carried out the rectangle encapsulation, and participate in normal wave field propagation and continuation between input reference face (relief surface) and output reference face (rectangular top) by one zero wave field layer is set.
3, the earthquake poststack forward modeling method of relief surface according to claim 1, it is characterized in that: adopt the Si Naier theorem to calculate reflection coefficient, and utilize Ricker wavelet to carry out convolution and handle, obtain reflectivity model, the dominant frequency of Ricker wavelet is selected according to the stratum characteristic in target work area, generally between 10Hz～90Hz.
4, the earthquake poststack forward modeling method of relief surface according to claim 1 and 2 is characterized in that: adopt the Fourier finite difference method operator of frequency field to carry out wave field extrapolation calculating from the bottom up along degree of depth axle.
CNB200710049482XA 2007-07-11 2007-07-11 The earthquake poststack forward modeling method of relief surface Expired - Fee Related CN100545679C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200710049482XA CN100545679C (en) 2007-07-11 2007-07-11 The earthquake poststack forward modeling method of relief surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200710049482XA CN100545679C (en) 2007-07-11 2007-07-11 The earthquake poststack forward modeling method of relief surface

Publications (2)

Publication Number Publication Date
CN101082676A true CN101082676A (en) 2007-12-05
CN100545679C CN100545679C (en) 2009-09-30

Family

ID=38912345

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200710049482XA Expired - Fee Related CN100545679C (en) 2007-07-11 2007-07-11 The earthquake poststack forward modeling method of relief surface

Country Status (1)

Country Link
CN (1) CN100545679C (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183790A (en) * 2011-02-12 2011-09-14 中国石油大学(华东) Elastic wave forward simulation technology based on space-time dual-variable grid
CN101738637B (en) * 2008-11-06 2012-01-04 北京北方林泰石油科技有限公司 Velocity change along with frequency information-based oil-gas detection method
CN102565856A (en) * 2010-12-29 2012-07-11 中国石油天然气集团公司 Near-surface noise suppression method on basis of wave equation forward modeling
CN102798887A (en) * 2011-05-27 2012-11-28 中国石油天然气集团公司 Filtration method for improving first arrival quality of controllable hypocenter record
CN102914799A (en) * 2012-10-12 2013-02-06 中国石油天然气股份有限公司 Forward modeling method and device for nonequivalent wave field
CN103091713A (en) * 2011-10-28 2013-05-08 中国石油化工股份有限公司 Automatic optimization method for undulating surface observation system
CN103592685A (en) * 2013-10-22 2014-02-19 中国石油天然气股份有限公司 Method and device for removing wave equation simulation direct wave in full waveform inversion
CN106772590A (en) * 2017-03-17 2017-05-31 中国地质科学院地球物理地球化学勘查研究所 A kind of free earth's surface finite-difference forward modeling system and method that acutely rises and falls
CN106950598A (en) * 2017-03-27 2017-07-14 中国科学院地质与地球物理研究所 A kind of migration velocity field method for evaluating reliability
CN107193043A (en) * 2017-05-15 2017-09-22 中国石油大学(华东) A kind of subsurface structure imaging method of relief surface
CN107479092A (en) * 2017-08-17 2017-12-15 电子科技大学 A kind of frequency domain high order ACOUSTIC WAVE EQUATION the Forward Modeling based on directional derivative
CN108828659A (en) * 2018-07-12 2018-11-16 中国石油天然气集团有限公司 Seismic wave field continuation method and device
CN109507722A (en) * 2017-09-15 2019-03-22 中国石油化工股份有限公司 Interbed multiple prediction technique and system based on model and dual wavefield continuation
CN109782354A (en) * 2019-02-27 2019-05-21 西安石油大学 Collaboration differential evolution algorithm based on direction guidance and its application in ray tracing
CN112147691A (en) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 Quick coding sequencing-free reference surface correction method and system

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102768367B (en) * 2012-07-04 2014-12-31 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Two-phase medium amplitude versus offset (AVO) forward modeling method based on triple constraints

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101738637B (en) * 2008-11-06 2012-01-04 北京北方林泰石油科技有限公司 Velocity change along with frequency information-based oil-gas detection method
CN102565856A (en) * 2010-12-29 2012-07-11 中国石油天然气集团公司 Near-surface noise suppression method on basis of wave equation forward modeling
CN102565856B (en) * 2010-12-29 2013-11-13 中国石油天然气集团公司 Near-surface noise suppression method on basis of wave equation forward modeling
CN102183790A (en) * 2011-02-12 2011-09-14 中国石油大学(华东) Elastic wave forward simulation technology based on space-time dual-variable grid
CN102798887A (en) * 2011-05-27 2012-11-28 中国石油天然气集团公司 Filtration method for improving first arrival quality of controllable hypocenter record
CN102798887B (en) * 2011-05-27 2016-03-30 中国石油天然气集团公司 A kind of filtering method improving vibroseis record first arrival quality
CN103091713A (en) * 2011-10-28 2013-05-08 中国石油化工股份有限公司 Automatic optimization method for undulating surface observation system
CN103091713B (en) * 2011-10-28 2015-11-18 中国石油化工股份有限公司 A kind of method to relief surface observing system Automatic Optimal
CN102914799A (en) * 2012-10-12 2013-02-06 中国石油天然气股份有限公司 Forward modeling method and device for nonequivalent wave field
CN102914799B (en) * 2012-10-12 2015-05-06 中国石油天然气股份有限公司 Forward modeling method and device for nonequivalent wave field
CN103592685A (en) * 2013-10-22 2014-02-19 中国石油天然气股份有限公司 Method and device for removing wave equation simulation direct wave in full waveform inversion
CN103592685B (en) * 2013-10-22 2016-04-06 中国石油天然气股份有限公司 Method and device for removing wave equation simulation direct wave in full waveform inversion
CN106772590A (en) * 2017-03-17 2017-05-31 中国地质科学院地球物理地球化学勘查研究所 A kind of free earth's surface finite-difference forward modeling system and method that acutely rises and falls
CN106950598A (en) * 2017-03-27 2017-07-14 中国科学院地质与地球物理研究所 A kind of migration velocity field method for evaluating reliability
CN106950598B (en) * 2017-03-27 2019-01-29 中国科学院地质与地球物理研究所 A kind of migration velocity field method for evaluating reliability
CN107193043A (en) * 2017-05-15 2017-09-22 中国石油大学(华东) A kind of subsurface structure imaging method of relief surface
CN107193043B (en) * 2017-05-15 2019-03-29 中国石油大学(华东) A kind of subsurface structure imaging method of relief surface
CN107479092A (en) * 2017-08-17 2017-12-15 电子科技大学 A kind of frequency domain high order ACOUSTIC WAVE EQUATION the Forward Modeling based on directional derivative
CN107479092B (en) * 2017-08-17 2019-02-12 电子科技大学 A kind of frequency domain high order ACOUSTIC WAVE EQUATION the Forward Modeling based on directional derivative
CN109507722B (en) * 2017-09-15 2020-11-13 中国石油化工股份有限公司 Model and dual-wavefield continuation-based interlayer multiple prediction method and system
CN109507722A (en) * 2017-09-15 2019-03-22 中国石油化工股份有限公司 Interbed multiple prediction technique and system based on model and dual wavefield continuation
CN108828659A (en) * 2018-07-12 2018-11-16 中国石油天然气集团有限公司 Seismic wave field continuation method and device
CN108828659B (en) * 2018-07-12 2020-02-14 中国石油天然气集团有限公司 Seismic wave field continuation method and device based on Fourier finite difference low-rank decomposition
CN109782354B (en) * 2019-02-27 2020-09-01 西安石油大学 Collaborative differential evolution algorithm based on direction guiding and application thereof in ray tracing
CN109782354A (en) * 2019-02-27 2019-05-21 西安石油大学 Collaboration differential evolution algorithm based on direction guidance and its application in ray tracing
CN112147691A (en) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 Quick coding sequencing-free reference surface correction method and system
CN112147691B (en) * 2019-06-28 2024-05-07 中国石油化工股份有限公司 Quick coding ordering-free reference plane correction method and system

Also Published As

Publication number Publication date
CN100545679C (en) 2009-09-30

Similar Documents

Publication Publication Date Title
CN100545679C (en) The earthquake poststack forward modeling method of relief surface
Boore et al. Comparison of two independent methods for the solution of wave‐scattering problems: Response of a sedimentary basin to vertically incident SH waves
CN102890290B (en) Pre-stack depth migration method under condition of undulating surface
CN102749643B (en) Method and device for acquiring frequency dispersion response of surface wave seismic record
CN103513277B (en) Seismic stratum fracture crack density inversion method and system
CN104570082B (en) Extraction method for full waveform inversion gradient operator based on green function characterization
CN108508482A (en) A kind of subterranean fracture seismic scattering response characteristic analogy method
CN107329171A (en) Depth domain reservoir stratum seismic inversion method and device
CN101545986A (en) Tridimensional integral prestack depth migration method based on maximum energy travel calculation
CN102116869A (en) High-precision prestack domain least square migration seismic imaging technology
CN102636809B (en) Method for generating spreading angle domain common image point gathers
CN104991268B (en) A kind of true amplitude migration imaging method
CN100349009C (en) Method for shifting depth before superposition in seismic data process of undulating the earth's surface
CN101021568A (en) Three-dimensional integral prestack depth migration method
CN103576200A (en) Low signal-to-noise ratio zone shallow wave impedance interface static correction method
CN105093319A (en) Ground micro-seismic static correction method based on three-dimensional seismic data
CN104570124A (en) Continuation imaging method suitable for cross-well seismic large-angle reflection conditions
CN104199088B (en) Incident angle gather extraction method and system
CN104166159A (en) Fracture conformation processing method of four-dimensional microseismic monitoring and system thereof
CN102778689B (en) Wide curved line seismic data underground reflection line building method
CN102565852B (en) Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir
CN106125139A (en) A kind of D seismic modeling method and system
CN104570090B (en) The extraction of full waveform inversion noise filter operator and the method filtered using its noise
CN102053269A (en) Analysis method of speed in seismic data
CN107942388A (en) A kind of triangle gridding reverse-time migration method in the case of mountain area earth's surface

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090930

Termination date: 20130711