WO2019242045A1 - Procédé de calcul de temps de déplacement d'onde sismique de construction de front d'onde bidimensionnel de source virtuelle - Google Patents

Procédé de calcul de temps de déplacement d'onde sismique de construction de front d'onde bidimensionnel de source virtuelle Download PDF

Info

Publication number
WO2019242045A1
WO2019242045A1 PCT/CN2018/094459 CN2018094459W WO2019242045A1 WO 2019242045 A1 WO2019242045 A1 WO 2019242045A1 CN 2018094459 W CN2018094459 W CN 2018094459W WO 2019242045 A1 WO2019242045 A1 WO 2019242045A1
Authority
WO
WIPO (PCT)
Prior art keywords
wavefront
travel time
seismic wave
ray
calculating
Prior art date
Application number
PCT/CN2018/094459
Other languages
English (en)
Chinese (zh)
Other versions
WO2019242045A9 (fr
Inventor
孟繁昌
孙辉
孟恩泓
黎佳宾
Original Assignee
成都启泰智联信息科技有限公司
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 成都启泰智联信息科技有限公司 filed Critical 成都启泰智联信息科技有限公司
Publication of WO2019242045A1 publication Critical patent/WO2019242045A1/fr
Publication of WO2019242045A9 publication Critical patent/WO2019242045A9/fr

Links

Images

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

Definitions

  • the invention relates to the field of seismic wave travel time calculation, in particular to a method for calculating seismic wave travel time of constructing a two-dimensional wavefront.
  • the existing two-dimensional wavefront construction method for seismic travel time calculation can improve the calculation accuracy to a certain extent, but the interpolation method does not consider the propagation law of the seismic wave, and the improved calculation accuracy is also limited.
  • the technical problem to be solved by the present invention is to provide a method for calculating the travel time of a seismic wavefront by constructing a virtual source wave, taking into account the propagation law of the seismic wave in the medium, and using a virtual source to calculate the travel time of grid nodes within the wavefront quadrilateral
  • the calculation method replaces the original two-dimensional cubic convolution interpolation method, and improves the calculation accuracy and stability of the travel time calculation method for seismic wavefront construction.
  • the technical solutions adopted by the present invention are:
  • a method for calculating travel time of a seismic wavefront with high accuracy including the following steps:
  • Step 1 Read in the relevant parameter file and velocity model, where the parameter file contains the grid points, grid spacing, source location, tracking step, ray angle range, and ray sampling interval of the speed model;
  • Step 2 Trace the ray and insert a new ray during the tracing process to ensure the coverage of the ray.
  • the essence of tracing rays is to use the Runge-Kutta method to solve the kinematic ray tracing equations, as shown in the following formula:
  • x i represents a position component
  • p i represents a slowness component
  • ⁇ represents a seismic wave travel time
  • v represents a seismic wave velocity
  • Step 3 Divide the computing space into multiple wavefront quads by calculating the spatial position information of the discrete points on the ray;
  • Step 4 Use the vector product method to determine the positional relationship between the grid points and the wavefront quadrilateral, and find the grid points contained in the wavefront quadrilateral;
  • Step 5 Calculate the virtual source position corresponding to each point by using the information about the apex of the wavefront quadrilateral, and calculate the seismic wave travel time of the grid points included in the wavefront quadrant based on these virtual source positions;
  • Step 6 Complete the travel time calculation of all grid nodes, and output the final travel time calculation result.
  • the present invention has the beneficial effects that, considering the propagation law of seismic waves in the medium, the virtual source method is used to calculate the travel time of the grid points inside the wavefront quadrangle, which improves the calculation accuracy of the travel times of the grid points, thereby further improving The calculation accuracy of the method for constructing seismic wave travel time of wavefront is presented.
  • FIG. 1 is a flowchart of a method for calculating travel time of constructing a seismic wavefront of a virtual source wavefront according to the present invention.
  • FIG. 2 is a schematic diagram of calculation region segmentation.
  • FIG 3 is a virtual computing a schematic view of seismic waves source, A, B, C, D quadrangular four vertices wavefront, O A, O B, O C, O D are A, B, C, D corresponding virtual source location , R is a grid point contained in the wavefront quadrilateral ABCD.
  • Figure 4 shows the relative travel time error of the conventional wavefront construction method in a homogeneous medium.
  • Figure 5 shows the relative error of the travel time calculation method for the construction of a virtual source wavefront in a homogeneous medium.
  • Fig. 1 is a flowchart of a method for calculating the travel time of constructing seismic wavefronts from a virtual source wavefront.
  • the figure shows the implementation process of the method of the present invention, which is as follows:
  • the essence of tracing rays is to use the Runge-Kutta method to solve the kinematic ray tracing equations, as shown in the following formula:
  • x i represents a position component
  • p i represents a slowness component
  • ⁇ represents a seismic wave travel time
  • v represents a seismic wave velocity
  • the position of the virtual source corresponding to these vertices is derived from the ray direction, seismic travel time, and seismic wave velocity of the vertices of the wavefront quadrilateral, as shown in Figure 3: Assume that the four vertices of a wavefront quadrilateral are A, B, C, D, R are the grid points inside. After calculating the positions of the virtual sources O A , O B , O C , and O D corresponding to A, B, C, and D through relevant information, the travel time expression of the seismic wave at point R is:
  • the calculation accuracy of the method of the present invention is analyzed and verified by a uniform model below.
  • Figures 4 and 5 show the absolute error of the conventional wavefront construction method and the virtual source wavefront construction method in a homogeneous medium model.
  • the model has a horizontal grid number of 761, a vertical grid number of 777, and a horizontal and vertical grid spacing. 10m, velocity is 1000m / s, the source is located at 3800m in the transverse direction. It can be seen from the figure that the calculation method of seismic wave travel time for virtual source wavefront construction has greatly improved the calculation accuracy compared with the conventional wavefront construction method.
  • the invention uses the virtual source method to calculate the travel time of the internal grid point of the wavefront quadrangle, improves the calculation accuracy of the travel time of the internal grid point of the wavefront quadrangle, and realizes a high-precision method for calculating the travel time of the seismic wavefront.

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

L'invention concerne un procédé de calcul d'un temps de déplacement d'onde sismique de construction de front d'onde de source virtuelle, consistant à: lire un fichier de paramètre associé et d'un modèle de vitesse; suivre un rayon, et insérer un nouveau rayon pendant le processus de suivi pour assurer le taux de couverture des rayons; diviser un espace de calcul en une pluralité de quadrilatères de front d'onde au moyen des informations de position spatiale de points discrets sur les rayons obtenus par calcul; déterminer la relation d'emplacement entre des points de grille et les quadrilatères de front d'onde au moyen d'un procédé de produit de vecteur, et trouver les points de grille compris dans les quadrilatères de front d'onde; calculer les emplacements de source virtuelle correspondant aux points au moyen d'informations associées à des sommets de quadrilatère de front d'onde, et calculer les temps de déplacement d'onde sismique des points de grille compris dans les quadrilatères de front d'onde sur la base de ces emplacements de source virtuelle; et achever le calcul de temps de déplacement de tous les nœuds de grille, et émettre un résultat de calcul de temps de déplacement final. Le procédé améliore la précision de calcul du temps de déplacement d'un point de grille.
PCT/CN2018/094459 2018-06-21 2018-07-04 Procédé de calcul de temps de déplacement d'onde sismique de construction de front d'onde bidimensionnel de source virtuelle WO2019242045A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201810641651.7 2018-06-21
CN201810641651.7A CN108957538A (zh) 2018-06-21 2018-06-21 一种虚拟震源二维波前构建地震波走时计算方法

Publications (2)

Publication Number Publication Date
WO2019242045A1 true WO2019242045A1 (fr) 2019-12-26
WO2019242045A9 WO2019242045A9 (fr) 2021-02-18

Family

ID=63141839

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/094459 WO2019242045A1 (fr) 2018-06-21 2018-07-04 Procédé de calcul de temps de déplacement d'onde sismique de construction de front d'onde bidimensionnel de source virtuelle

Country Status (5)

Country Link
CN (1) CN108957538A (fr)
BE (1) BE1025828B1 (fr)
LU (1) LU100878B1 (fr)
NL (1) NL2021354B1 (fr)
WO (1) WO2019242045A1 (fr)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108957538A (zh) * 2018-06-21 2018-12-07 成都启泰智联信息科技有限公司 一种虚拟震源二维波前构建地震波走时计算方法
CN110660135A (zh) * 2019-09-20 2020-01-07 西南石油大学 一种利用三角网格化射线实现波前构建的方法
CN110568496B (zh) * 2019-09-26 2021-02-09 核工业北京地质研究院 一种复杂介质条件下射线追踪方法
CN111257939B (zh) * 2020-03-26 2021-06-01 中国石油大学(北京) 一种时移地震虚拟震源双向波场重构方法和系统
CN114924312B (zh) * 2022-05-10 2024-03-12 吉林大学 一种基于波前初始化射线追踪技术的高斯波束偏移方法及装置

Citations (3)

* 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
CN105974470A (zh) * 2016-07-04 2016-09-28 中国石油大学(华东) 一种多分量地震资料最小二乘逆时偏移成像方法及系统
US20170363759A1 (en) * 2016-06-17 2017-12-21 Cgg Services Sa System and method for seismic interferometry optimized data acquisition

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1423735B1 (fr) * 2001-09-07 2010-08-18 Shell Internationale Research Maatschappij B.V. Representation en image d'une formation sous la surface au moyen de sources virtuelles
CN108072897A (zh) * 2018-01-23 2018-05-25 西南交通大学 一种混合二维地震波走时计算方法
CN108957538A (zh) * 2018-06-21 2018-12-07 成都启泰智联信息科技有限公司 一种虚拟震源二维波前构建地震波走时计算方法

Patent Citations (3)

* 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
US20170363759A1 (en) * 2016-06-17 2017-12-21 Cgg Services Sa System and method for seismic interferometry optimized data acquisition
CN105974470A (zh) * 2016-07-04 2016-09-28 中国石油大学(华东) 一种多分量地震资料最小二乘逆时偏移成像方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HAN, FUXING ET AL.: "Research Status of the Wavefront Construction Method", PROGRESS IN GEOPHYSICS, vol. 26, no. 3, 30 June 2011 (2011-06-30), pages 1045 - 1051 *
N. ETTRICH ET AL.: "Wave Front Construction in Smooth Media for Prestack Depth Migration", PURE AND APPLIED GEOPHYSICS, vol. 148, 30 September 1996 (1996-09-30), pages 3 - 4 *
WANG, DONGHE ET AL.: "Review of the Seismic Ray Tracing Method", PROGRESS IN GEOPHYSICS, vol. 31, no. 1, 31 January 2016 (2016-01-31), pages 344 - 353 *

Also Published As

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

Similar Documents

Publication Publication Date Title
WO2019242045A1 (fr) Procédé de calcul de temps de déplacement d'onde sismique de construction de front d'onde bidimensionnel de source virtuelle
CN107038308B (zh) 一种基于线性内插的规则格网地形建模方法
CN109116413B (zh) 成像域立体层析速度反演方法
CN116774292B (zh) 一种地震波走时确定方法、系统、电子设备及存储介质
CN108267781B (zh) 任意曲面非均匀介质快速行进程函方程求解射线追踪算法
CN108303736B (zh) 各向异性ti介质最短路径射线追踪正演方法
CN112099082B (zh) 一种共面元共方位角道集的地震回折波走时反演方法
CN116359981A (zh) 一种地震波走时快速确定方法
CN108845355A (zh) 地震偏移成像方法及装置
JP6019230B2 (ja) 地震水平線の決定を向上させるための方法
CN111948708B (zh) 一种浸入边界起伏地表地震波场正演模拟方法
CN109581494B (zh) 叠前偏移方法及装置
Chambers et al. A practical implementation of wave front construction for 3-D isotropic media
CN109387872B (zh) 表面多次波预测方法
Li et al. A simple and stable centeredness measure for 3D curve skeleton extraction
CN114924312A (zh) 一种基于波前初始化射线追踪技术的高斯波束偏移方法及装置
Jensen et al. Ray methods
CN110390702B (zh) 基于地质走向分析的变面元平面插值方法
Dietz B-spline approximation with energy constraints
Wang et al. Analytic double product integrals for all-frequency relighting
WO2002065159A1 (fr) Procede d'imagerie sismique dans lequel on utilise le calcul du temps de parcours direct
CN112596103A (zh) 射线追踪方法、装置和电子设备
CN109636839A (zh) 局部泊松曲面约束的实物表面采样点集法向估计方法
CN114429047B (zh) 一种基于三角网格的二次圆方程旅行时插值方法
CN114428292B (zh) 近地表速度模型的构建方法和存储介质

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18923542

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18923542

Country of ref document: EP

Kind code of ref document: A1