LU100878B1 - Method for calculating the 2-D seismic propagation time on the basis of a wavefront construction with virtual sources - Google Patents

Method for calculating the 2-D seismic propagation time on the basis of a wavefront construction with virtual sources Download PDF

Info

Publication number
LU100878B1
LU100878B1 LU100878A LU100878A LU100878B1 LU 100878 B1 LU100878 B1 LU 100878B1 LU 100878 A LU100878 A LU 100878A LU 100878 A LU100878 A LU 100878A LU 100878 B1 LU100878 B1 LU 100878B1
Authority
LU
Luxembourg
Prior art keywords
wavefront
propagation time
quadrilaterals
calculating
seismic
Prior art date
Application number
LU100878A
Other languages
French (fr)
Inventor
Jiabin Li
Original Assignee
Chengdu Qitai Zhilian Information Tech Co Ltd
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 Qitai Zhilian Information Tech Co Ltd filed Critical Chengdu Qitai Zhilian Information Tech Co Ltd
Application granted granted Critical
Publication of LU100878B1 publication Critical patent/LU100878B1/en

Links

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. analysis, for interpretation, for correction
    • 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
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G01V20/00—
    • 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

Abstract

L'invention concerne un procédé de calcul de temps de propagation sismique 2-D basé sur une construction de front d'onde à sources virtuelles. Le procédé de calcul du temps de propagation sismique comprend les étapes suivantes : lire des paramètres pertinents et un modèle de vitesse ; tracer des rayons dans différentes directions provenant d'une source sismique à l'aide d'un procédé de Runge-Kutta, et insérer de nouveaux rayons pendant le traçage de rayons en fonction des conditions d'évaluation pertinentes dans le processus de traçage des rayons ; diviser un espace de calcul en une pluralité de quadrilatères de front d'onde via des points discrets sur les rayons adjacents et sur des surface de front d'onde adjacentes ; trouver des points de grille dans les quadrilatères de front d'onde ; calculer le temps de propagation du point de grille dans chacun des quadrilatères de front d'onde par un procédé à sources virtuelles ; et terminer le calcul du temps de propagation des points de la grille dans tous les quadrilatères de front d'onde. Grâce à l'adoption du procédé à sources virtuelles, lors du calcul du temps de propagation des points de grille dans les quadrilatères de front d'onde de l'invention, la précision de calcul du temps de propagation des points de grille est améliorée et un procédé de calcul sismique basé sur la construction de front d'onde est réalisé.The invention relates to a 2-D seismic propagation time calculation method based on a wavefront construction with virtual sources. The method for calculating the seismic propagation time comprises the following steps: reading relevant parameters and a speed model; trace rays in different directions from a seismic source using a Runge-Kutta process, and insert new rays during ray tracing according to the relevant evaluation conditions in the ray tracing process ; dividing a computing space into a plurality of wavefront quadrilaterals via discrete points on the adjacent rays and on adjacent wavefront surfaces; find grid points in the wavefront quadrilaterals; calculating the propagation time of the grid point in each of the wavefront quadrilaterals by a method with virtual sources; and finish calculating the propagation time of the grid points in all the wavefront quadrilaterals. Thanks to the adoption of the method with virtual sources, during the calculation of the propagation time of the grid points in the wavefront quadrilaterals of the invention, the accuracy of calculation of the propagation time of the grid points is improved and a seismic calculation method based on the wavefront construction is carried out.

LU100878A 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 LU100878B1 (en)

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
LU100878B1 true LU100878B1 (en) 2019-12-30

Family

ID=63141839

Family Applications (1)

Application Number Title Priority Date Filing Date
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

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)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
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
CN110660135A (en) * 2019-09-20 2020-01-07 西南石油大学 Method for realizing wavefront construction by utilizing triangular gridding rays
CN110568496B (en) * 2019-09-26 2021-02-09 核工业北京地质研究院 Ray tracing method under complex medium condition
CN111257939B (en) * 2020-03-26 2021-06-01 中国石油大学(北京) Time-lapse seismic virtual source bidirectional wave field reconstruction method and system
CN114924312B (en) * 2022-05-10 2024-03-12 吉林大学 Gaussian beam migration method and device based on wavefront initialization ray tracing technology

Family Cites Families (6)

* Cited by examiner, † Cited by third party
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
EP1423735B1 (en) * 2001-09-07 2010-08-18 Shell Internationale Research Maatschappij B.V. Seismic imaging a subsurface formation by means of virtual sources
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
CN108072897A (en) * 2018-01-23 2018-05-25 西南交通大学 It is a kind of to mix computational methods when two-dimensionally seismic wave is walked
CN108957538A (en) * 2018-06-21 2018-12-07 成都启泰智联信息科技有限公司 A kind of virtual focus two dimension wavefront construction seimic travel time calculation method

Also Published As

Publication number Publication date
BE1025828A1 (en) 2019-07-18
CN108957538A (en) 2018-12-07
WO2019242045A9 (en) 2021-02-18
NL2021354B1 (en) 2019-04-26
NL2021354A (en) 2018-08-14
BE1025828B1 (en) 2019-07-25
WO2019242045A1 (en) 2019-12-26

Similar Documents

Publication Publication Date Title
LU100878B1 (en) Method for calculating the 2-D seismic propagation time on the basis of a wavefront construction with virtual sources
Kimmel et al. Skeletonization via distance maps and level sets
KR101305305B1 (en) System and method for measuring surface image velocity using correlation analysis of spatio-temporal images
CN110717231B (en) Sub-basin confluence simulation method based on slope channel river channel three-level structure
CN105761310B (en) A kind of sunykatuib analysis and image display method of sky visible range numerical map
CN104181593B (en) Three-dimensional ray-free tracing contraflexure wave tomographic imaging method and device
CN113313826B (en) Scale-adaptive visual domain analysis method
CN103997636B (en) A kind of method calculating hiding relation between hologram three-dimensional display object
ATE421737T1 (en) METHOD FOR GENERATING A 2D IMAGE OF A 3D OBJECT
CN110674593A (en) Method for automatically predicting over-short excavation in tunnel excavation process
CN116774292A (en) Seismic wave travel time determining method, system, electronic equipment and storage medium
CN108680968B (en) Evaluation method and device for seismic exploration data acquisition observation system in complex structural area
ES2845280T3 (en) Post-stacking Kirchhoff depth demigration method for transverse inclined isotropic (TTI) and heterogeneous media based on ray tracing on migrated data
CN110263428B (en) Riverbed evolution analysis method based on flow weighted average streamline length index
CN106780704A (en) Based on the direct lighting effect proximity rendering method of three-dimensional scenic that observability is reused
CN106844889B (en) Quadrilateral structure grid division method based on non-uniform block principle
CN109359360B (en) Structural stress processing method based on local features
CN112002012A (en) Visibility analysis method for urban area
Bento et al. Image-based techniques for the advanced characterization of scour around bridge piers in laboratory
Kim et al. A leak detection and 3D source localization method on a plant piping system by using multiple cameras
KR101237128B1 (en) Development of tomographic ptv
CN110378917B (en) Tooth segmentation method based on peak value clustering
CN107229067B (en) Stratified model anisotropy ray-tracing procedure and system
CN109614585B (en) Novel ionosphere region reconstruction method
CN111899297A (en) Line structure light stripe center extraction method

Legal Events

Date Code Title Description
FG Patent granted

Effective date: 20191230