CN104597488A - Optimum design method of finite difference template of non-equiangular long-grid wave equation - Google Patents
Optimum design method of finite difference template of non-equiangular long-grid wave equation Download PDFInfo
- Publication number
- CN104597488A CN104597488A CN201510029000.9A CN201510029000A CN104597488A CN 104597488 A CN104597488 A CN 104597488A CN 201510029000 A CN201510029000 A CN 201510029000A CN 104597488 A CN104597488 A CN 104597488A
- Authority
- CN
- China
- Prior art keywords
- finite difference
- wave
- optimization
- difference coefficient
- beta
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention discloses an optimum design method of a finite difference template of a non-equiangular long-grid wave equation. The method comprises the following steps: according to a time sampling interval and a spatial sampling interval, mesh subdivision is executed to the simulating region of an actual geologic model; according to the maximum appointed permissible error and the wave number range, the corresponding operator lengths of different acoustic velocities are calculated; based on the finite difference method of the least square optimization time space domain and the corresponding operator length, the optimized finite difference coefficient of the grid is acquired; the acquired optimized finite difference coefficient is brought into the wave equation with a difference scheme to perform forward modeling of the wave equation. According to the method, the finite difference method of the time space domain only applied to square and cube grids is expanded into the rectangular or rectangular parallelepiped grids, thus a special simulation precision requirement and a requirement for saving a calculated amount are met in actual production.
Description
Technical field
The present invention relates to seismic forward modeling numerical simulation technology field, particularly the non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method of one.
Background technology
Seismic forward modeling numerical simulation technology is when complex underground structure group (comprising isotropic medium, anisotropic medium, the heterogeneous anisotropic medium of Biot, random hole medium etc.) is known, numerical computation method is utilized to make ripple at this Propagation, through the repeatedly transmission of subsurface geological structure, reflection, scattering, by the process that the wave detector of earth's surface or underground receives.Accurate wave equation numerical solution is utilized to carry out the seismic response of simulate formation complex geological structure, for many aspects such as Study of Seismic wave traveling mechanism, the special treatment method of seismic data and the explanation of bad ground provide the mathematical physics foundation of more science.In recent years, Wave Equation Numerical method is widely used in reverse-time migration and full waveform inversion.
Wave equation forward modeling has multiple method, comparatively common are: finite difference method, pseudo-spectrometry, finite element method, boundary element method, spectral element method etc.Wherein finite difference method because its calculated amount is little, counting yield is high, more complicated rate pattern can be adapted to and be widely used.Method of finite difference can be divided into according to different standards: Explicit finite difference and implicit expression finite difference; Regular grid finite difference, staggering mesh finite-difference and rotationally staggered grid finite difference.In method of finite difference, difference coefficient can be tried to achieve by Taylor series expansion or optimization method, respectively corresponding finite difference based on Taylor series expansion and the finite difference based on optimization.In conventional finite method of difference, difference coefficient is obtained by the dispersion relation of minimization spatial domain.In recent years, occurred a kind of time-space domain method of finite difference, the method asks for difference coefficient by the dispersion relation of minimization time domain and spatial domain, has higher simulation precision and better stability.
Current time-space domain finite difference method requires that the spatial sampling interval in all directions is equal, and namely needing model facetization is square or square grid.And in actual production, in order to meet specific accuracy requirement or in order to save calculated amount, we usually need to be rectangle or rectangular parallelepiped grid (normally the mesh spacing of depth direction is different from horizontal direction) by model facetization.Be applicable to the time-space domain finite difference method of square and square grid, specific accuracy requirement can not be met and maybe can not save calculated amount.
Summary of the invention
Embodiments provide a kind of non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method, the time-space domain finite difference method being only applicable to square and square grid is expanded in rectangle or rectangular parallelepiped grid, meet the demand of specific simulation precision requirement and saving calculated amount in actual production, the method comprises:
According to time sampling interval and spatial sampling interval, mesh generation is carried out to the simulated domain of actual geologic model;
According to the given limits of error and wave-number range, corresponding operator length is obtained to different acoustic velocity;
Based on the time-space domain method of finite difference of Least-squares minimization and the operator length of described correspondence, obtain the optimization finite difference coefficient of grid;
The optimization finite difference coefficient obtained is brought into the wave equation of difference scheme, carry out Wave equation forward modeling.
In one embodiment, describedly according to time sampling interval and spatial sampling interval, difference gridding subdivision is carried out to the simulated domain of actual geologic model, comprising: the simulated domain of actual geologic model is split into rectangular parallelepiped grid or rectangular node;
Based on the time-space domain method of finite difference of Least-squares minimization and the operator length of described correspondence, obtain the optimization finite difference coefficient of grid, comprising: obtain the optimization finite difference coefficient of rectangular parallelepiped grid or the optimization finite difference coefficient of rectangular node.
In one embodiment, the optimization finite difference coefficient of described acquisition rectangular parallelepiped grid, according to following formulae discovery:
Wherein,
Wherein, b is wave number, and M is operator length, a
mfor the finite difference coefficient after optimization, θ is the angle of plane wave propagation direction and surface level, θ ∈ [0, π]; φ is the position angle of plane wave propagation, φ ∈ [0,2 π];
v is acoustic velocity, and τ is time sampling interval, and h is x, y direction sampling interval; C=Δ z/h, c are parametric variable, and Δ z is z direction sampling interval; β=kh, k are parametric variable, and β is wave-number range, β ∈ [0, b]; M is parametric variable, and m is integer, m ∈ [1, M]; N is parametric variable, and n is integer, n ∈ [1, M].
In one embodiment, the maximum error that the optimization finite difference coefficient of described rectangular parallelepiped grid is corresponding meets following constraint condition:
ξ
1max<η;
Wherein, ξ
1maxfor the maximum error that the optimization finite difference coefficient of rectangular parallelepiped grid is corresponding; η is that maximum permission mistake 5 is poor.
In one embodiment, the maximum error ξ that the optimization finite difference coefficient of described rectangular parallelepiped grid is corresponding
1maxbe calculated as follows:
Wherein,
In one embodiment, the described wave equation optimization finite difference coefficient obtained being brought into difference scheme, carry out Wave equation forward modeling, comprise: the optimization finite difference coefficient of the rectangular parallelepiped grid obtained is substituted in the three-dimensional acoustic wave wave equation of difference scheme, carries out three-dimensional acoustic wave Wave equation forward modeling; The three-dimensional acoustic wave wave equation of described difference scheme is:
Wherein, P is acoustic pressure.The optimization finite difference coefficient of described acquisition rectangular node, according to following formulae discovery:
Wherein,
Wherein, b is wave number, and M is operator length, a
mfor the finite difference coefficient after optimization, θ is the angle of plane wave propagation direction and surface level, θ ∈ [0, π];
v is acoustic velocity, and τ is time sampling interval, and h is x direction sampling interval; C=Δ z/h, c are parametric variable, and Δ z is z direction sampling interval; β=kh, k are parametric variable, and β is wave-number range, β ∈ [0, b]; M is parametric variable, and m is integer, m ∈ [1, M]; N is parametric variable, and n is integer, n ∈ [1, M].
In one embodiment, the maximum error that the optimization finite difference coefficient of described rectangular node is corresponding meets following constraint condition:
ξ
2max<η;
Wherein, ξ
2maxfor the maximum error that the optimization finite difference coefficient of rectangular node is corresponding; η is the limits of error.
In one embodiment, the maximum error ξ that the optimization finite difference coefficient of described rectangular node is corresponding
2maxbe calculated as follows:
Wherein,
In one embodiment, described the optimization finite difference coefficient obtained is brought in the wave equation of difference scheme, carry out Wave equation forward modeling, also comprise: the optimization finite difference coefficient of the rectangular node obtained is substituted in the two-dimentional Acoustic Wave-equation of difference scheme, carries out two-dimentional Acoustic Wave-equation forward simulation; The two-dimentional Acoustic Wave-equation of described difference scheme is:
Wherein, P is acoustic pressure.
In embodiments of the present invention, propose a kind of non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method, the time-space domain finite difference method being only applicable to square and square grid extend in rectangle or rectangular parallelepiped grid by the method, meets the demand of specific simulation precision requirement and saving calculated amount in actual production.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:
Fig. 1 is the non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method flow diagram of one that the embodiment of the present invention provides;
Fig. 2 is that the graph of errors obtained by traditional finite difference method that provides of the embodiment of the present invention is with the Changing Pattern figure of the direction of propagation;
Fig. 3 is the Changing Pattern figure of the graph of errors obtained by non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method that provides of the embodiment of the present invention with the direction of propagation;
Fig. 4 is that the graph of errors obtained by traditional finite difference method that provides of the embodiment of the present invention is with the Changing Pattern figure of operator length;
Fig. 5 is the Changing Pattern figure of the graph of errors obtained by non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method that provides of the embodiment of the present invention with operator length;
Fig. 6 is the wave field snapshot comparison diagram of the even sound wave medium obtained respectively by traditional finite difference method and non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method that provides of the embodiment of the present invention in the 1.0s moment;
Fig. 7 is the comparison of wave shape figure of dotted line position in Fig. 6;
Fig. 8 is the wave field snapshot comparison diagram of the Marmousi model obtained respectively by traditional finite difference method and non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method that provides of the embodiment of the present invention in the 1.0s moment;
Embodiment
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with embodiment and accompanying drawing, the present invention is described in further details.At this, exemplary embodiment of the present invention and illustrating for explaining the present invention, but not as a limitation of the invention.
Inventor finds, current time-space domain finite difference method requires that the spatial sampling interval in all directions is equal, namely needing model facetization is square or square grid, but this kind of method can not meet specific accuracy requirement maybe can not save calculated amount.If be rectangle or rectangular parallelepiped grid by model facetization, then can meet specific accuracy requirement or save calculated amount.Based on this, the present invention proposes a kind of non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method.
Fig. 1 is the non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method flow diagram of one that the embodiment of the present invention provides, and as shown in Figure 1, the method comprises:
Step 101: mesh generation is carried out to the simulated domain of actual geologic model according to time sampling interval and spatial sampling interval;
Step 102: according to the given limits of error and wave-number range, obtains corresponding operator length to different acoustic velocity;
Step 103: based on the time-space domain method of finite difference of Least-squares minimization and the operator length of described correspondence, obtains the optimization finite difference coefficient of grid;
Step 104: the wave equation optimization finite difference coefficient obtained being brought into difference scheme, carries out Wave equation forward modeling.
During concrete enforcement, need to be rectangular node or rectangular parallelepiped grid (normally the mesh spacing of depth direction is different from horizontal direction) by the simulated domain of actual geologic model (or zoning) subdivision.
Three-dimensional acoustic wave wave equation is as follows:
In formula, P represents acoustic pressure, and V represents acoustic velocity.
Be different from the situation in x, y direction for z direction mesh spacing, the optimization finite difference scheme of concrete rectangular parallelepiped grid is as follows:
Wherein, τ is time sampling interval; H is x, y direction sampling interval; C=Δ z/h, c are parametric variable, and Δ z is z direction sampling interval; a
mfor the finite difference coefficient after optimization, M represents operator length; M is parametric variable, and m is integer, m ∈ [1, M].
In order to contracted notation, we will
be abbreviated as
according to Plane wave theory, Wo Menling:
Wherein,
Wherein, i, l, j, n are parametric variable, and n is integer, n ∈ [1, M]; K is parametric variable; The angle that θ ∈ [0, π] is plane wave propagation direction and surface level, the position angle that φ ∈ [0,2 π] is plane wave propagation.Equation (3), (4) are substituted into equation (2), and suitable abbreviation, can obtain:
Wherein,
again according to constraint condition;
Can obtain:
Wherein, β=kh, β are wave-number range, β ∈ [0, b].Optimization finite-difference algorithm will be obtained a set of difference coefficient exactly and make equation (7) two ends error minimum.According to Least Square Theory, this problem can be converted into and solve following system of equations:
Wherein,
By solving system of linear equations (8), we can obtain the optimization difference coefficient of rectangular parallelepiped grid.Again difference coefficient is updated in equation (2), just can carries out solving of three-dimensional acoustic wave wave equation.
In order to ensure the validity of the inventive method, need the constraint condition added in the optimization procedure above below, the maximum error that namely the optimization finite difference coefficient of rectangular parallelepiped grid is corresponding meets following constraint condition:
ξ
1max<η (10)
We pass through maximum error corresponding to following formula calculation optimization difference coefficient:
Wherein,
Wherein, η is the limits of error.
ξ
1maxonly relevant with b and M.When M fixes, ξ
1maxbe decided by b, ξ
1maxincrease along with b and become large.Therefore, if initial b does not meet equation (10), we should by reducing b until ξ gradually
1max< η obtains best b.In addition, when b fixes, ξ
1maxonly relevant with M, ξ
1maxincrease along with M and diminish.So if initial M does not meet equation (10), we should by increasing M until ξ gradually
1max< η obtains best M.
For two-dimensional rectangle grid, be different from the situation in x direction for z direction mesh spacing, concrete optimization finite difference scheme is as follows:
In formula, P represents acoustic pressure, and V represents acoustic velocity; Wherein, τ is time sampling interval; H is x direction sampling interval; C=Δ z/h, c are parametric variable, and Δ z is z direction sampling interval; a
mfor the finite difference coefficient after optimization, M represents operator length; M is parametric variable, and m is integer, m ∈ [1, M].
We can obtain the difference coefficient of the optimization of rectangular node by solving following linear equation:
Wherein,
By solving system of linear equations (15), we can obtain the optimization difference coefficient of rectangular node.Again difference coefficient is updated in equation (13), just can carries out solving of two-dimentional Acoustic Wave-equation.
Same, in order to ensure the validity of the inventive method, need the constraint condition added in the optimization procedure above below, the maximum error that namely the optimization finite difference coefficient of rectangular node is corresponding meets following constraint condition:
ξ
2max<η (16)
We pass through maximum error corresponding to following formula calculation optimization difference coefficient:
Wherein,
Wherein, η is the limits of error.
For an even sound wave medium, advantage of the present invention is described below.Velocity of longitudinal wave is 1500m/s, and time sampling interval is 1ms, x direction sampling interval be 10m, z direction sampling interval is 5m, operator length M ∈ [2,10].
Regulation b=2.74 and M=8, the inventive method (i.e. non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method) and traditional finite difference method is adopted to obtain the Changing Pattern of graph of errors with the direction of propagation, respectively as shown in Figures 2 and 3.Regulation η=10
-2.5≈ 0.003, adopts the inventive method and traditional finite difference method to obtain the Changing Pattern of graph of errors with operator length, respectively as shown in Figures 4 and 5.From Fig. 2 to Fig. 5, compared with the time-space domain finite difference method (traditional finite difference method) based on Taylor, the inventive method has wider effective band and less numerical solidification.Regulation η=10
-2.5≈ 0.003, adopts traditional finite difference method and the inventive method to obtain the wave field snapshot comparison diagram of even sound wave medium, as shown in Figure 6.Wherein, the left side figure in Fig. 6 is the wave field snapshot adopting traditional finite difference method to obtain even sound wave medium, and the right figure is the wave field snapshot of the even sound wave medium adopting the inventive method to obtain.Fig. 7 is the comparison of wave shape figure of dotted line position in Fig. 6.From Fig. 6 and Fig. 7, when operator length is identical, the inventive method simulation precision is higher.
In order to further illustrate validity of the present invention, we just drill with Marmousi model.Time sampling interval is 1ms, x direction sampling interval be 12m, z direction sampling interval is 9m.The operator length M=9 of traditional method of finite difference.Time-space domain finite difference adopts the algorithm becoming operator length, selects the operator of different length, M ∈ [2,9] according to different speed.Fig. 8 is the wave field snapshot comparison diagram of the Marmousi model obtained respectively by traditional finite difference method and the inventive method that provides of the embodiment of the present invention in the 1.0s moment, upper graph be the Marmousi model that obtained by traditional finite difference method at the wave field snapshot in 1.0s moment, lower edge graph is the Marmousi model that obtained respectively by the inventive method wave field snapshot in the 1.0s moment; As shown in Figure 8, the inventive method is applicable to the forward simulation of complex dielectrics, and its simulate effect is obviously better than traditional method of finite difference.
In sum, the inventive method has the following advantages: 1. the numerical solidification reducing medium-high frequency section.2. more realistic based on the dispersion relation of time-space domain, simulation precision is higher.3. be applicable to rectangle and rectangular parallelepiped grid, the demand of specific accuracy requirement and saving calculated amount in actual production can be met.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. a non-equilateral fourdrinier wire lattice wave equation finite difference optimum formwork design method, is characterized in that, comprising:
According to time sampling interval and spatial sampling interval, mesh generation is carried out to the simulated domain of actual geologic model;
According to the given limits of error and wave-number range, corresponding operator length is obtained to different acoustic velocity;
Based on the time-space domain method of finite difference of Least-squares minimization and the operator length of described correspondence, obtain the optimization finite difference coefficient of grid;
The optimization finite difference coefficient obtained is brought into the wave equation of difference scheme, carry out Wave equation forward modeling.
2. the method for claim 1, it is characterized in that, describedly according to time sampling interval and spatial sampling interval, difference gridding subdivision is carried out to the simulated domain of actual geologic model, comprising: the simulated domain of actual geologic model is split into rectangular parallelepiped grid or rectangular node;
Based on the time-space domain method of finite difference of Least-squares minimization and the operator length of described correspondence, obtain the optimization finite difference coefficient of grid, comprising: obtain the optimization finite difference coefficient of rectangular parallelepiped grid or the optimization finite difference coefficient of rectangular node.
3. method as claimed in claim 2, is characterized in that, the optimization finite difference coefficient of described acquisition rectangular parallelepiped grid, according to following formulae discovery:
Wherein,
Wherein, b is wave number, and M is operator length, a
mfor the finite difference coefficient after optimization, θ is the angle of plane wave propagation direction and surface level, θ ∈ [0, π]; φ is the position angle of plane wave propagation, φ ∈ [0,2 π];
v is acoustic velocity, and τ is time sampling interval, and h is x, y direction sampling interval; C=Δ z/h, c are parametric variable, and Δ z is z direction sampling interval; β=kh, k are parametric variable, and β is wave-number range, β ∈ [0, b]; M is parametric variable, and m is integer, m ∈ [1, M]; N is parametric variable, and n is integer, n ∈ [1, M].
4. method as claimed in claim 3, it is characterized in that, the maximum error that the optimization finite difference coefficient of described rectangular parallelepiped grid is corresponding meets following constraint condition:
ξ
1max<η;
Wherein, ξ
1maxfor the maximum error that the optimization finite difference coefficient of rectangular parallelepiped grid is corresponding; η is the limits of error.
5. method as claimed in claim 4, is characterized in that, the maximum error ξ that the optimization finite difference coefficient of described rectangular parallelepiped grid is corresponding
1maxbe calculated as follows:
Wherein,
6. method as claimed in claim 5, it is characterized in that, the described wave equation optimization finite difference coefficient obtained being brought into difference scheme, carry out Wave equation forward modeling, comprise: the optimization finite difference coefficient of the rectangular parallelepiped grid obtained is substituted in the three-dimensional acoustic wave wave equation of difference scheme, carries out three-dimensional acoustic wave Wave equation forward modeling; The three-dimensional acoustic wave wave equation of described difference scheme is:
Wherein, P is acoustic pressure.
7. method as claimed in claim 2, is characterized in that, the optimization finite difference coefficient of described acquisition rectangular node, according to following formulae discovery:
Wherein,
Wherein, b is wave number, and M is operator length, a
mfor the finite difference coefficient after optimization, θ is the angle of plane wave propagation direction and surface level, θ ∈ [0, π];
v is acoustic velocity, and τ is time sampling interval, and h is x direction sampling interval; C=Δ z/h, c are parametric variable, and Δ z is z direction sampling interval; β=kh, k are parametric variable, and β is wave-number range, β ∈ [0, b]; M is parametric variable, and m is integer, m ∈ [1, M]; N is parametric variable, and n is integer, n ∈ [1, M].
8. method as claimed in claim 7, it is characterized in that, the maximum error that the optimization finite difference coefficient of described rectangular node is corresponding meets following constraint condition:
ξ
2max<η;
Wherein, ξ
2maxfor the maximum error that the optimization finite difference coefficient of rectangular node is corresponding; η is the limits of error.
9. method as claimed in claim 8, is characterized in that, the maximum error ξ that the optimization finite difference coefficient of described rectangular node is corresponding
2maxbe calculated as follows:
Wherein,
10. method as claimed in claim 9, it is characterized in that, described the optimization finite difference coefficient obtained is brought in the wave equation of difference scheme, carry out Wave equation forward modeling, also comprise: the optimization finite difference coefficient of the rectangular node obtained is substituted in the two-dimentional Acoustic Wave-equation of difference scheme, carries out two-dimentional Acoustic Wave-equation forward simulation; The two-dimentional Acoustic Wave-equation of described difference scheme is:
Wherein, P is acoustic pressure.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510029000.9A CN104597488B (en) | 2015-01-21 | 2015-01-21 | Optimum design method of finite difference template of non-equiangular long-grid wave equation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510029000.9A CN104597488B (en) | 2015-01-21 | 2015-01-21 | Optimum design method of finite difference template of non-equiangular long-grid wave equation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104597488A true CN104597488A (en) | 2015-05-06 |
CN104597488B CN104597488B (en) | 2017-05-24 |
Family
ID=53123393
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510029000.9A Active CN104597488B (en) | 2015-01-21 | 2015-01-21 | Optimum design method of finite difference template of non-equiangular long-grid wave equation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104597488B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106646597A (en) * | 2016-12-14 | 2017-05-10 | 中国石油大学(北京) | Forward modeling method and device based on spring network model |
CN107179549A (en) * | 2017-07-11 | 2017-09-19 | 中海石油(中国)有限公司 | A kind of acoustic wave equation in time domain Explicit finite difference seismic response analogy method |
CN107976710A (en) * | 2017-11-17 | 2018-05-01 | 河海大学 | A kind of implicit time-space domain finite difference numerical simulation method of linear optimization based on ACOUSTIC WAVE EQUATION |
CN108279437A (en) * | 2018-01-17 | 2018-07-13 | 中国石油大学(华东) | Variable density ACOUSTIC WAVE EQUATION time higher order accuracy staggering mesh finite-difference method |
CN109490955A (en) * | 2018-11-14 | 2019-03-19 | 深圳市勘察研究院有限公司 | A kind of the Acoustic Wave-equation the Forward Modeling and device of rule-based grid |
CN109490947A (en) * | 2018-10-16 | 2019-03-19 | 中国科学院地质与地球物理研究所 | A kind of high-temperature medium seimic wave propagation analogy method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102062875A (en) * | 2010-11-30 | 2011-05-18 | 中国石油集团川庆钻探工程有限公司 | Forward modeling method for fluctuating surface elastic wave equation |
CN102183790A (en) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | Elastic wave forward simulation technology based on space-time dual-variable grid |
US20120243371A1 (en) * | 2011-03-23 | 2012-09-27 | Chevron U.S.A. Inc. | System and method for seismic data modeling and migration |
CN103630933A (en) * | 2013-12-09 | 2014-03-12 | 中国石油天然气集团公司 | Nonlinear optimization based time-space domain staggered grid finite difference method and device |
-
2015
- 2015-01-21 CN CN201510029000.9A patent/CN104597488B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102062875A (en) * | 2010-11-30 | 2011-05-18 | 中国石油集团川庆钻探工程有限公司 | Forward modeling method for fluctuating surface elastic wave equation |
CN102183790A (en) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | Elastic wave forward simulation technology based on space-time dual-variable grid |
US20120243371A1 (en) * | 2011-03-23 | 2012-09-27 | Chevron U.S.A. Inc. | System and method for seismic data modeling and migration |
CN103630933A (en) * | 2013-12-09 | 2014-03-12 | 中国石油天然气集团公司 | Nonlinear optimization based time-space domain staggered grid finite difference method and device |
Non-Patent Citations (2)
Title |
---|
李胜军 等: "声波方程有限差分数值模拟的变网格步长算法", 《工程地球物理学报》 * |
金其虎: "弹性波有限差分数值模拟及井间地震逆时偏移成像研究", 《中国优秀硕士学位论文全文数据库·基础科学辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106646597A (en) * | 2016-12-14 | 2017-05-10 | 中国石油大学(北京) | Forward modeling method and device based on spring network model |
CN107179549A (en) * | 2017-07-11 | 2017-09-19 | 中海石油(中国)有限公司 | A kind of acoustic wave equation in time domain Explicit finite difference seismic response analogy method |
CN107976710A (en) * | 2017-11-17 | 2018-05-01 | 河海大学 | A kind of implicit time-space domain finite difference numerical simulation method of linear optimization based on ACOUSTIC WAVE EQUATION |
CN107976710B (en) * | 2017-11-17 | 2019-05-28 | 河海大学 | A kind of implicit time-space domain finite difference numerical simulation method of linear optimization based on ACOUSTIC WAVE EQUATION |
CN108279437A (en) * | 2018-01-17 | 2018-07-13 | 中国石油大学(华东) | Variable density ACOUSTIC WAVE EQUATION time higher order accuracy staggering mesh finite-difference method |
CN109490947A (en) * | 2018-10-16 | 2019-03-19 | 中国科学院地质与地球物理研究所 | A kind of high-temperature medium seimic wave propagation analogy method |
CN109490947B (en) * | 2018-10-16 | 2020-01-10 | 中国科学院地质与地球物理研究所 | High-temperature medium seismic wave propagation simulation method |
CN109490955A (en) * | 2018-11-14 | 2019-03-19 | 深圳市勘察研究院有限公司 | A kind of the Acoustic Wave-equation the Forward Modeling and device of rule-based grid |
Also Published As
Publication number | Publication date |
---|---|
CN104597488B (en) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104597488A (en) | Optimum design method of finite difference template of non-equiangular long-grid wave equation | |
Zhang et al. | Parameter prediction of hydraulic fracture for tight reservoir based on micro-seismic and history matching | |
CN102222365B (en) | Method for reconstructing curved surface of complex space | |
CN107203002A (en) | The preparation method of the picture of inversion speed model and its method for building up and underground structure | |
CN105182444A (en) | High resolution sequence stratigraphic framework constraint geostatistical inversion method | |
CN103901478A (en) | Method for determining deposition characteristics and distribution of reservoirs by combining logging and seismic information | |
CN102636809B (en) | Method for generating spreading angle domain common image point gathers | |
CN103033846A (en) | Seismic inversion system and seismic inversion method controlled by geologic facies | |
CN103913774A (en) | Reservoir stratum geological mechanics parameter retrieval method based on micro seismic event | |
CN106814391A (en) | Ground micro-seismic state event location method based on Fresnel zone tomographic inversion | |
CN105549077A (en) | Micro-earthquake epicenter positioning method calculated based on multilevel multi-scale grid similarity coefficient | |
CN103838936A (en) | High-precision tectonic stress field simulation method applicable to turbidite sand low-permeability reservoirs | |
Yao et al. | Tuning fractures with dynamic data | |
CN104316958A (en) | Coherent processing method for identifying different scales of formation fractures | |
CN105676280A (en) | Two-phase medium geological data obtaining method and device based on rotationally staggered grids | |
CN104330822A (en) | Method and device for determining residual oil-gas distribution by adopting coupled four-dimensional seismic inversion | |
CN104297792B (en) | The phased inversion method of water channel reservoir is stacked on a kind of fan | |
CN105573963A (en) | Reconstruction method for horizontal nonuniform structure of ionized layer | |
CN108279437A (en) | Variable density ACOUSTIC WAVE EQUATION time higher order accuracy staggering mesh finite-difference method | |
Ou et al. | 3D visualization of hydraulic fractures using micro-seismic monitoring: Methodology and application | |
CN105607119A (en) | Near-surface model construction method and static correction value calculation method | |
CN102768367B (en) | Two-phase medium amplitude versus offset (AVO) forward modeling method based on triple constraints | |
CN106443793A (en) | Space-time bivariant forward modeling method | |
CN103217715B (en) | Multiple dimensioned regular grid Static Correction of Tomographic Inversion method | |
CN102841374B (en) | Pseudo three-dimensional fast microseism forward modeling method based on scanning surface forward modeling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |