CN108957538A - A kind of virtual focus two dimension wavefront construction seimic travel time calculation method - Google Patents
A kind of virtual focus two dimension wavefront construction seimic travel time calculation method Download PDFInfo
- Publication number
- CN108957538A CN108957538A CN201810641651.7A CN201810641651A CN108957538A CN 108957538 A CN108957538 A CN 108957538A CN 201810641651 A CN201810641651 A CN 201810641651A CN 108957538 A CN108957538 A CN 108957538A
- Authority
- CN
- China
- Prior art keywords
- wavefront
- ray
- quadrangle
- travel time
- point
- 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
- 238000010276 construction Methods 0.000 title claims abstract description 27
- 238000004364 calculation method Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000005070 sampling Methods 0.000 claims description 4
- 230000002596 correlated effect Effects 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/125—Virtual source
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/57—Trace interpolation or extrapolation, e.g. for virtual receiver; Anti-aliasing for missing receivers
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- 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/67—Wave propagation modeling
- G01V2210/671—Raytracing
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of virtual focus wavefront construction seimic travel time calculation methods, comprising the following steps: reads in relevant parameter, rate pattern；Ray is tracked using runge kutta method from focus along different directions, and is inserted into new ray according to correlated judgment condition during ray tracing；Multiple wavefront quadrangles are divided into space is calculated by the discrete point on before adjacent ray and adjacent wave；Search out wavefront quadrangle intranet lattice point；When wavefront quadrangle intranet lattice point be calculated by virtual focus method walking；Complete calculating when all wavefront quadrangle intranet lattice points are walked.When the present invention is walked by using virtual focus method calculating wavefront quadrangle intranet lattice point, computational accuracy when mesh point is walked is improved, a kind of high-precision wavefront construction seimic travel time calculation method is realized.
Description
Technical field
The present invention relates to seimic travel time calculating field, especially a kind of two-dimentional wavefront construction seimic travel time calculating side
Method.
Background technique
" Jilin University's journal " (geoscience version) the 2nd phase in 2008 discloses Han Fuxing etc. " based on two-dimentional cubic convolution
The ray-Tracing of Wavefront Construction of interpolation algorithm " describes a kind of improved wavefront construction seimic travel time calculation method, by two
Dimension cubic convolution interpolation method has been applied to when mesh point is walked in calculating, calculates essence so as to improve the seimic travel time obtained
The computational efficiency of degree and algorithm.And by uniform dielectric to the wavefront construction side based on two-dimentional cubic convolution interpolation algorithm
Method has carried out error analysis, and analysis result has obtained relatively good effect.
" calculating physics " the 2nd phase in 2008 discloses Han Fuxing etc., and " different interpolation algorithms are in ray-Tracing of Wavefront Construction
Application and comparison ", comparative analysis adjacent domains interpolation method, bilinear interpolation hair, fragment interpolation method and two-dimentional three secondary volumes
Application effect of the product interpolation method in ray-Tracing of Wavefront Construction, the model calculation reflect two-dimentional cubic convolution interpolation side
Method is relative to other available more accurate ray paths of three kinds of methods.
" Advances in Geophysics " the 5th phase in 2009 discloses " the mesh point relative positioning in wavefront construction method such as Han Fuxing
And attribute calculates research ", it describes and how vector multiplication cross method to be used to judge rectangular mesh node and irregular wavefront quadrangle
Relative positional relationship, and according to the positional relationship of mesh point and irregular wavefront quadrangle, give corresponding interpolation method
Calculate the attribute information of mesh point.And the method for proposition is verified by homogeneous model and calculated examples, is achieved
Good calculated result.
By example above as can be seen that existing two dimension wavefront construction seimic travel time calculation method to a certain extent
It is able to ascend computational accuracy, but there is no the propagation law for considering seismic wave, the computational accuracies of promotion for interpolation method therein
Also limited.
Summary of the invention
Technical problem to be solved by the invention is to provide a kind of virtual focus wavefront construction seimic travel time calculating sides
Method, it is contemplated that the propagation law of seismic wave in the medium uses during when grid node is walked in calculating wavefront quadrangle
A kind of virtual hypocenter computing process improves wavefront construction seismic wave to replace original two-dimentional cubic convolution interpolation method
The computational accuracy and stability of calculation method when walking.
In order to solve the above technical problems, the technical solution adopted by the present invention is that:
A kind of high-precision wavefront construction seimic travel time calculation method, comprising the following steps:
Step 1: reading in relevant parameter file, rate pattern, wherein the Parameter File includes the grid of rate pattern
Points, grid spacing, hypocentral location, tracking step-length, x-ray angle range, ray sampling interval；
Step 2: tracking ray, and in tracing process it is inserted into new ray and guarantees with this coverage rate of ray.Tracking is penetrated
The essence of line is to solve kinematics ray tracing equation group using runge kutta method, is shown below:
Wherein, xiIndicate location components, piIndicate that slowness component, τ indicate that seimic travel time, v indicate seimic wave velocity.
Step 3: the spatial positional information by calculating discrete point on the ray obtained is divided into multiple waves for space is calculated
Preceding quadrangle；
Step 4: judging the positional relationship of mesh point Yu wavefront quadrangle by vector multiplications method, find wavefront quadrangle
The mesh point for being included；
Step 5: virtual hypocentral location corresponding to each point is gone out by wavefront quadrangle vertex associated information calculation, and
The seimic travel time of the included mesh point of wavefront quadrangle is calculated based on these virtual hypocentral locations；
Step 6: completing calculating when walking of all grid nodes, and export final calculated result when walking.
Compared with prior art, it the beneficial effects of the present invention are: considering the propagation law of seismic wave in the medium, adopts
When being walked with virtual focus method calculating wavefront quadrangle intranet lattice point, computational accuracy when mesh point is walked is improved, and then improve
The computational accuracy of wavefront construction seimic travel time method.
Detailed description of the invention
Fig. 1 is the virtual focus wavefront construction seimic travel time calculation method flow chart of the present invention.
Divide schematic diagram in the zoning Fig. 2.
Fig. 3 is that virtual focus seimic travel time calculates schematic diagram, and A, B, C, D are four vertex of wavefront quadrangle, OA、OB、
OC、ODThe respectively corresponding virtual hypocentral location A, B, C, D, R are the mesh point that wavefront quadrangle ABCD includes.
Fig. 4 is relative error when conventional wavefront construction method is walked in uniform dielectric.
Fig. 5 is virtual focus wavefront construction seimic travel time calculation method relative error in uniform dielectric.
Specific embodiment
The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.
Fig. 1 is virtual focus wavefront construction seimic travel time calculation method flow chart, shows the method for the present invention in figure
Implementation process, specific as follows:
1) relevant parameter file, rate pattern are read in, wherein Grid dimension of the Parameter File comprising rate pattern,
Grid spacing, hypocentral location, tracking step-length, x-ray angle range, ray sampling interval；
2) ray is tracked, and in tracing process is inserted into new ray and guarantees with this coverage rate of ray.Ray emission angle
Spend range are as follows: -80 ° to+80 °, the sampling interval is 3 ° to 6 °, and ray tracing step-length is 2ms to 6ms.Tracking ray essence be
Kinematics ray tracing equation group is solved using runge kutta method, is shown below:
Wherein, xiIndicate location components, piIndicate that slowness component, τ indicate that seimic travel time, v indicate seimic wave velocity；
3) by calculating the spatial positional information of discrete point on the ray obtained, on before adjacent ray adjacent wave
Four points are wavefront quadrangle vertex, will calculate space and are divided into multiple wavefront quadrangles (as shown in Figure 2)；
4) mesh point for being likely located at its coverage area is substantially filtered out by wavefront quadrangle vertex position first, with
The positional relationship for judging mesh point Yu wavefront quadrangle by vector multiplications method afterwards determines the net that wavefront quadrangle is included
Lattice point.
5) these are extrapolated by information such as the directions of rays on wavefront quadrangle vertex, seimic travel time, seimic wave velocities
The position of the virtual focus of vertex correspondence, as shown in Figure 3: assuming that four vertex of a certain wavefront quadrangle are respectively A, B, C, D, R
For the mesh point of its inside.Go out the corresponding virtual focus O of A, B, C, D by associated information calculationA、OB、OC、ODBehind position, R point
Seimic travel time expression formula are as follows:
Wherein | OBR|、|OBR|、|OCR|、|ODR | respectively indicate OA、OB、OC、ODTo the distance of R point, VRIndicate what R was pointed out
Seimic wave velocity.
6) calculating when walking of all grid nodes is completed, and exports final calculated result when walking.
Analysis verifying is carried out below by computational accuracy of the homogeneous model to the method for the present invention.
Fig. 4, Fig. 5 are respectively that conventional wavefront construction method and virtual focus wavefront construction method are exhausted in uniform dielectric model
To error, model transverse grid points are 761, and longitudinal grid points are 777, and lateral longitudinal grid spacing is 10m, speed
For 1000m/s, focus is located at transverse direction 3800m.As can be seen from the figure virtual focus wavefront construction seimic travel time calculating side
Method has large increase relative to conventional wavefront construction method computational accuracy.
When the present invention is walked by using virtual focus method calculating wavefront quadrangle intranet lattice point, four side of wavefront is improved
Computational accuracy when shape intranet lattice point is walked realizes a kind of high-precision wavefront construction seimic travel time calculation method.
Claims (2)
1. a kind of virtual focus two dimension wavefront construction seimic travel time calculation method, which comprises the following steps:
Step 1: reading in relevant parameter file, rate pattern, wherein Grid dimension of the Parameter File comprising rate pattern,
Grid spacing, hypocentral location, tracking step-length, x-ray angle range, ray sampling interval；
Step 2: tracking ray, and in tracing process it is inserted into new ray and guarantees with this coverage rate of ray.Track ray
Essence is to solve kinematics ray tracing equation group using runge kutta method, is shown below:
Wherein, xiIndicate location components, piIndicate that slowness component, τ indicate that seimic travel time, v indicate seimic wave velocity.
Step 3: the spatial positional information by calculating discrete point on the ray obtained is divided into multiple wavefront four for space is calculated
Side shape；
Step 4: the positional relationship of mesh point Yu wavefront quadrangle is judged by vector multiplications method, finding wavefront quadrangle is included
Mesh point；
Step 5: virtual hypocentral location corresponding to each point being gone out by wavefront quadrangle vertex associated information calculation, and is based on this
A little virtual hypocentral locations calculate the seimic travel time of the included mesh point of wavefront quadrangle；
Step 6: completing calculating when walking of all grid nodes, and export final calculated result when walking.
2. a kind of virtual focus two dimension wavefront construction seimic travel time calculation method as described in claim 1, which is characterized in that
In steps of 5, the seimic travel time of the included mesh point of wavefront quadrangle is calculated by virtual focus method.Assuming that a certain wavefront
Four vertex of quadrangle are respectively A, B, C, D, and R is its internal mesh point.It is corresponding to go out A, B, C, D by associated information calculation
Virtual focus OA、OB、OC、ODBehind position, the seimic travel time expression formula of R point are as follows:
Wherein | OBR|、|OBR|、|OCR|、|ODR | respectively indicate OA、OB、OC、ODTo the distance of R point, VR indicates the earthquake that R is pointed out
Wave velocity.
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810641651.7A CN108957538A (en) | 2018-06-21 | 2018-06-21 | A kind of virtual focus two dimension wavefront construction seimic travel time calculation method |
| PCT/CN2018/094459 WO2019242045A1 (en) | 2018-06-21 | 2018-07-04 | Method for calculating virtual source two-dimensional wavefront construction seismic wave travel time |
| LU100878A LU100878B1 (en) | 2018-06-21 | 2018-07-04 | Method for calculating the 2-D seismic propagation time on the basis of a wavefront construction with virtual sources |
| NL2021354A NL2021354B1 (en) | 2018-06-21 | 2018-07-19 | 2-D Seismic Travel Time Calculation Method Based on Virtual Source Wavefront Construction |
| BE2018/5525A BE1025828B1 (en) | 2018-06-21 | 2018-07-19 | Method for calculating 2-D seismic propagation time based on a virtual source wavefront construction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810641651.7A CN108957538A (en) | 2018-06-21 | 2018-06-21 | A kind of virtual focus two dimension wavefront construction seimic travel time calculation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108957538A true CN108957538A (en) | 2018-12-07 |
Family
ID=63141839
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810641651.7A Pending CN108957538A (en) | 2018-06-21 | 2018-06-21 | A kind of virtual focus two dimension wavefront construction seimic travel time calculation method |
Country Status (5)
| Country | Link |
|---|---|
| CN (1) | CN108957538A (en) |
| BE (1) | BE1025828B1 (en) |
| LU (1) | LU100878B1 (en) |
| NL (1) | NL2021354B1 (en) |
| WO (1) | WO2019242045A1 (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110568496A (en) * | 2019-09-26 | 2019-12-13 | 核工业北京地质研究院 | ray tracing method under complex medium condition |
| CN110660135A (en) * | 2019-09-20 | 2020-01-07 | 西南石油大学 | Method for realizing wavefront construction by utilizing triangular gridding rays |
| CN111257939A (en) * | 2020-03-26 | 2020-06-09 | 中国石油大学(北京) | Time-lapse seismic virtual source bidirectional wave field reconstruction method and system |
| CN112379413A (en) * | 2020-10-28 | 2021-02-19 | 中国石油天然气集团有限公司 | Irregular seismic source characterization method and device based on energy spectrum equivalence |
| CN114924312A (en) * | 2022-05-10 | 2022-08-19 | 吉林大学 | Gaussian beam offset method and device based on wave-front initialization ray tracing technology |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108957538A (en) * | 2018-06-21 | 2018-12-07 | 成都启泰智联信息科技有限公司 | A kind of virtual focus two dimension wavefront construction seimic travel time calculation method |
| CN115201901A (en) * | 2022-06-30 | 2022-10-18 | 中铁第四勘察设计院集团有限公司 | Method, device and equipment for determining tunnel wave front travel time and readable storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003023450A1 (en) * | 2001-09-07 | 2003-03-20 | Shell Internationale Research Maatschappij B.V. | Seismic imaging a subsurface formation by means of virtual sources |
| CN108072897A (en) * | 2018-01-23 | 2018-05-25 | 西南交通大学 | It is a kind of to mix computational methods when two-dimensionally seismic wave is walked |
| NL2021354A (en) * | 2018-06-21 | 2018-08-14 | Chengdu Qitai Zhilian Information Tech Co Ltd | 2-D Seismic Travel Time Calculation Method Based on Virtual Source Wavefront Construction |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6035256A (en) * | 1997-08-22 | 2000-03-07 | Western Atlas International, Inc. | Method for extrapolating traveltimes across shadow zones |
| US20170363759A1 (en) * | 2016-06-17 | 2017-12-21 | Cgg Services Sa | System and method for seismic interferometry optimized data acquisition |
| CN105974470B (en) * | 2016-07-04 | 2017-06-16 | 中国石油大学(华东) | A kind of multi-component seismic data least square reverse-time migration imaging method and system |
-
2018
- 2018-06-21 CN CN201810641651.7A patent/CN108957538A/en active Pending
- 2018-07-04 LU LU100878A patent/LU100878B1/en active IP Right Grant
- 2018-07-04 WO PCT/CN2018/094459 patent/WO2019242045A1/en active Application Filing
- 2018-07-19 NL NL2021354A patent/NL2021354B1/en not_active IP Right Cessation
- 2018-07-19 BE BE2018/5525A patent/BE1025828B1/en not_active IP Right Cessation
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003023450A1 (en) * | 2001-09-07 | 2003-03-20 | Shell Internationale Research Maatschappij B.V. | Seismic imaging a subsurface formation by means of virtual sources |
| CN108072897A (en) * | 2018-01-23 | 2018-05-25 | 西南交通大学 | It is a kind of to mix computational methods when two-dimensionally seismic wave is walked |
| NL2021354A (en) * | 2018-06-21 | 2018-08-14 | Chengdu Qitai Zhilian Information Tech Co Ltd | 2-D Seismic Travel Time Calculation Method Based on Virtual Source Wavefront Construction |
Non-Patent Citations (2)
| Title |
|---|
| 何洋: "基于波前构建的射线走时和振幅计算", 《中国优秀硕士学位论文全文数据库（硕士）》 * |
| 王东鹤: "地震波射线追踪方法研究综述", 《地球物理学进展》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110660135A (en) * | 2019-09-20 | 2020-01-07 | 西南石油大学 | Method for realizing wavefront construction by utilizing triangular gridding rays |
| CN110568496A (en) * | 2019-09-26 | 2019-12-13 | 核工业北京地质研究院 | ray tracing method under complex medium condition |
| CN111257939A (en) * | 2020-03-26 | 2020-06-09 | 中国石油大学(北京) | Time-lapse seismic virtual source bidirectional wave field reconstruction method and system |
| CN112379413A (en) * | 2020-10-28 | 2021-02-19 | 中国石油天然气集团有限公司 | Irregular seismic source characterization method and device based on energy spectrum equivalence |
| CN114924312A (en) * | 2022-05-10 | 2022-08-19 | 吉林大学 | Gaussian beam offset method and device based on wave-front initialization ray tracing technology |
| CN114924312B (en) * | 2022-05-10 | 2024-03-12 | 吉林大学 | Gaussian beam migration method and device based on wavefront initialization ray tracing technology |
Also Published As
| Publication number | Publication date |
|---|---|
| NL2021354A (en) | 2018-08-14 |
| WO2019242045A1 (en) | 2019-12-26 |
| NL2021354B1 (en) | 2019-04-26 |
| WO2019242045A9 (en) | 2021-02-18 |
| BE1025828A1 (en) | 2019-07-18 |
| BE1025828B1 (en) | 2019-07-25 |
| LU100878B1 (en) | 2019-12-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108957538A (en) | A kind of virtual focus two dimension wavefront construction seimic travel time calculation method | |
| Li et al. | Fourier neural operator with learned deformations for pdes on general geometries | |
| CN108180918B (en) | Point cloud geodesic path forward tracking generation method and device | |
| CN106918838A (en) | Gaussian beam offset imaging method and device under the conditions of relief surface | |
| CN104570106A (en) | Near-surface tomographic velocity analysis method | |
| CN106501852A (en) | A kind of multiple dimensioned full waveform inversion method of three-dimensional acoustic wave equation arbitrarily-shaped domain and device | |
| CN108267781B (en) | Ray tracing algorithm for solving fast travel function equation of non-uniform medium of any curved surface | |
| CN106886980A (en) | A kind of enhanced method of point cloud density based on three-dimensional laser radar target identification | |
| CN108303736B (en) | Ray tracing forward method for shortest path of anisotropic TI medium | |
| CN109799530A (en) | Rayleigh waves dispersion curve inversion method for seismic surface wave exploration | |
| CN104360396B (en) | A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well | |
| CN102692644A (en) | Depth domain common-image gather generation method | |
| CN102338887B (en) | Irregular-size space-variant grid tomography imaging statics correction method | |
| CN107886569A (en) | It is a kind of that controllable surface parameterization method and system are estimated based on discrete Lie derivatives | |
| WO2013152221A1 (en) | Converting a first acquired data subset to a second acquired data subset | |
| CN110045333A (en) | A kind of sound source three-dimensional positioning method based on Kalman filtering | |
| Zhou et al. | An iterative factored topography-dependent eikonal solver for anisotropic media | |
| CN106842314B (en) | The determination method of formation thickness | |
| CN108802819B (en) | A kind of trapezoidal grid finite difference Simulation of Seismic Wave method of uniform depth sampling | |
| CN113221395A (en) | Construction method and application of seismic travel time parameterized grid model based on layered medium | |
| Yu et al. | A minimum traveltime ray tracing global algorithm on a triangular net for propagating plane waves | |
| CN109387872B (en) | Surface multiple prediction method | |
| Amodei et al. | Computation of uniform wave forms using complex rays | |
| CN107561583A (en) | Local angle computational methods and imaging method for Gaussian beam pre-stack depth migration | |
| CN108280863A (en) | TIN model equivalence line search methods based on bilateral structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181207 |