WO2019071504A1 - 一种基于两点射线追踪的地震走时层析反演方法 - Google Patents

一种基于两点射线追踪的地震走时层析反演方法 Download PDF

Info

Publication number
WO2019071504A1
WO2019071504A1 PCT/CN2017/105817 CN2017105817W WO2019071504A1 WO 2019071504 A1 WO2019071504 A1 WO 2019071504A1 CN 2017105817 W CN2017105817 W CN 2017105817W WO 2019071504 A1 WO2019071504 A1 WO 2019071504A1
Authority
WO
WIPO (PCT)
Prior art keywords
travel time
wave
velocity
layer
seismic
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.)
Ceased
Application number
PCT/CN2017/105817
Other languages
English (en)
French (fr)
Inventor
方鑫定
陈晓非
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southern University of Science and Technology
Original Assignee
Southern University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southern University of Science and Technology filed Critical Southern University of Science and Technology
Priority to CN201780001180.7A priority Critical patent/CN108064348B/zh
Priority to PCT/CN2017/105817 priority patent/WO2019071504A1/zh
Priority to US15/950,350 priority patent/US20190113641A1/en
Publication of WO2019071504A1 publication Critical patent/WO2019071504A1/zh
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

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
    • 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. 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
    • 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 exploration, in particular to a seismic travel time tomography inversion method based on two-point ray tracing.
  • Seismic tomography is a method of inverting the velocity structure of a region using seismic observation data.
  • the principle of seismic tomography is similar to medical CT technology.
  • the inversion calculation of the travel time or waveform of the observed elastic wave in the rock and soil medium is carried out.
  • the ray tracing method based on ray theory is a kind of forward algorithm.
  • the ray tracing method is based on high frequency approximation, and only calculates the ray path from the source point to the receiving point, and the calculation efficiency is high.
  • the inversion calculation generally uses an optimization algorithm, such as the steepest descent method, the conjugate gradient method, the Newton iterative method, the random search, etc., to solve the optimal solution satisfying the given minimum error criterion.
  • seismic tomography can describe the velocity structure of the study region using a one-dimensional velocity model, that is, the velocity is only a function of depth.
  • the formation velocity variation in a small area can generally be approximated by a one-dimensional model.
  • the engineering geophysical will generally adopt a seismic layer based on the one-dimensional model. Analyze the imaging method. The existing one-dimensional travel time tomographic velocity inversion assumes that the velocity in each layer is uniform. Through multiple iterations, a layered velocity model is obtained. Generally, it is necessary to divide more layers to describe the actual formation more accurately. The characteristics of speed changes.
  • the object of the present invention is to provide a seismic walk-time tomography inversion method based on two-point ray tracing, which aims to solve the problem of the computational complexity of seismic tomography in the prior art and the existing two-point traveltime ray tracing Newton iterative method.
  • a method for seismic traveltime tomography inversion based on two-point ray tracing comprising the following steps:
  • the ray parameter p is represented by the variable q
  • the ray path is uniquely determined by the ray parameter p
  • the theoretical direct wave travel time and the reflected wave travel time are calculated;
  • the one-dimensional continuous layered model with continuously variable velocity in the layer is adjusted by an optimization algorithm until the calculated theoretical direct wave travel time and reflected wave travel time and actual direct travel time data and reflected waves The difference in travel time data satisfies the given error criteria and the model output is performed.
  • the two-point ray tracing-based seismic traveltime tomography inversion method has the following beneficial effects compared with the prior art: by establishing a one-dimensional continuous layered model in which the formation velocity can be continuously changed, The inverse evolution quantity based on this model can be greatly reduced, and the actual formation velocity structure can be described more accurately, and the inversion calculation efficiency can be significantly improved.
  • the variable q is used indirectly to solve the ray parameter p.
  • the iterative solution process is stable, the convergence is fast, and the problem that the iteration does not converge at the large incident angle is effectively avoided.
  • FIG. 1 is a flow chart of a method for seismic travel time tomography inversion based on two-point ray tracing according to an embodiment of the present invention
  • 2 is a parameter definition of a continuous layered model of a seismic travel time tomography inversion method based on two-point ray tracing according to an embodiment of the present invention
  • FIG. 3 is a schematic diagram of a stratigraphic model of a seismic traveltime tomography inversion method based on two-point ray tracing according to an embodiment of the present invention
  • FIG. 4 is a schematic diagram of data acquisition of seismic exploration surface data based on a two-point ray tracing seismic travel time tomography inversion method according to an embodiment of the present invention
  • FIG. 5 is a schematic diagram of a surface-excited vertical seismic section of a seismic travel time tomography inversion method based on two-point ray tracing according to an embodiment of the present invention
  • FIG. 6 is a schematic diagram of a simulated vertical seismic section in a well based on a two-point ray tracing-based seismic travel time tomography inversion method according to an embodiment of the present invention
  • FIG. 7 is a schematic diagram of cross-well imaging of a seismic travel time tomography inversion method based on two-point ray tracing according to an embodiment of the present invention.
  • an embodiment of the present invention provides a method for seismic traveltime tomography inversion based on two-point ray tracing, the method comprising the following steps:
  • Step S1 Collect actual travel time data. Perform seismic data acquisition in the study area to obtain direct wave travel time data and reflected wave travel time data;
  • Step S2 Model parameterization. Modeling the study area and establishing an initial one-dimensional continuous layered model with continuous changes in the velocity within the layer;
  • the velocity in the layer may be a longitudinal wave or a transverse wave velocity, which is increased or decreased with depth;
  • Step S3 Calculate the ray parameters.
  • the ray parameter p is represented by the variable q
  • Step S4 Calculate the theoretical travel time. Calculating the direct wave travel time and the reflected wave travel time according to the ray parameter p;
  • Step S5 Compare the theoretical travel time with the actual travel time. Comparing the calculated theoretical direct wave travel time and reflected wave travel time with the actual direct travel time data and reflected wave travel time data obtained by the seismic data acquisition, and judging the theoretical direct travel time and reflected wave travel time and the actual data obtained by the seismic data acquisition Straight Whether the difference between the arrival time data of the wave and the reflected wave travel time meets the given error criterion, if yes, proceed to step S7, otherwise proceed to step S6;
  • Step S6 Optimize the model.
  • the optimization algorithm is used to adjust the one-dimensional continuous layered model whose velocity in the layer can be continuously changed until the calculated theoretical direct wave travel time and reflected wave travel time and the actual direct wave travel time data and the reflected wave travel time data satisfy the given The error standard is up to now.
  • the optimization algorithm used may be a steepest descent method, a conjugate gradient method, a Newton iterative method, a random search, or the like;
  • Step S7 Model output. After the calculated theoretical time-to-time and actual data arrival time meets the given error criteria, the model output is performed.
  • the two-point ray tracing-based seismic travel time tomography inversion method provides a one-dimensional continuous layered model in which the intra-layer velocity can be continuously changed, and the intra-layer velocity is expressed as a function of depth, and the model allows intra-layer
  • the continuous change of speed greatly reduces the number of layers, which greatly reduces the amount of anti-evolution, can more accurately describe the actual formation velocity structure, and significantly improve the efficiency of inversion calculation; and by expressing the ray parameter p with the variable q, indirectly Using the variable q to solve the ray parameter p, the iterative solution process is stable and the convergence is fast, effectively avoiding the problem that the iteration does not converge under the condition of large incident angle.
  • the intra-layer velocity Vk is a function of the depth z
  • the k-th layer velocity function is expressed as:
  • V k a k z+b k ,
  • a k and b k are the model parameters to be inverted, respectively the gradient and intercept of the k-th velocity function, and when the formation longitudinal velocity model is inverted, the intra-layer velocity V k
  • the intra-layer velocity V k For the longitudinal wave velocity; when the formation shear wave velocity model is inverted, the intra-layer velocity V k is the shear wave velocity; when the formation is converted to the wave velocity model, the intra-layer velocity V k is determined by the intra-layer converted wave property as the longitudinal wave velocity.
  • the shear wave velocity where the converted wave is a longitudinal wave to a transverse wave converted wave or a transverse wave to a longitudinal wave converted wave, and the converted wave property in the layer is a longitudinal wave property or a transverse wave property.
  • the source distance X is expressed as a function of the ray parameter p, and the source distance X is expressed as:
  • ⁇ k , ⁇ k , h k , ⁇ k and ⁇ k are intermediate parameters
  • the ray parameter p is represented by the variable q as: Where V M is the maximum velocity of the simulated ray path through the formation.
  • the source distance X is expressed as a function of the variable q, and the source distance X is expressed as:
  • the coefficients ⁇ 1 , ⁇ 2 , ⁇ 1 and ⁇ 2 are obtained by the Taylor expansion of the rational function formula and the function formula of the source X.
  • the initial value estimation formula of q is obtained:
  • the initial value is obtained by using the initial value estimation formula of the q, and the iterative calculation is performed, and the exact solution of q is obtained through iteration.
  • the initial value obtained by this method can generally reach more than 95%, and only 1 to 3 iterations can be obtained to obtain the exact solution of q.
  • the calculation formula of the direct wave travel time and the reflected wave travel time is:
  • the steps of seismic data acquisition in the study area are specifically as follows:
  • the source and the receiving detector are arranged on the surface;
  • the seismic signal used for seismic data acquisition is a longitudinal wave, or a transverse wave, or a converted wave of a longitudinal wave to a transverse wave, or a converted wave of a transverse wave to a longitudinal wave.

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

一种基于两点射线追踪的地震走时层析反演方法，包括步骤：地震数据采集（S1）；建立层内速度可连续变化的初始一维连续层状模型（S2）；将射线参数p用变量q表示，震源距X表示为变量q的函数X＝f(q)，利用牛顿迭代法求解X＝f(q) （S3），根据射线参数p计算直达波走时和反射波走时（S4）；比较理论到时和实际数据到时（S5），用最优化算法调整模型速度参数，直到理论到时和实际到时的差别满足给定的误差标准为止（S6）。通过建立地层速度可连续变化的一维连续层状模型，大大减少了划分层数，使实际地层速度结构描述地更加准确，提高了反演计算效率；并且通过利用变量q求解射线参数p，保证了在大入射角的情况下也能快速稳定收敛。

Description

一种基于两点射线追踪的地震走时层析反演方法 技术领域
本发明涉及地震勘探领域，尤其涉及一种基于两点射线追踪的地震走时层析反演方法。
背景技术
地震层析成像是指利用地震观测数据反演研究区域速度结构的一种方法。地震层析成像的原理类似于医学CT技术，是根据弹性波理论及其在地层介质中的传播规律，对观测得到的弹性波在岩土体介质中的走时或波形进行反演计算，重建被测范围内岩体弹性波速度分布规律的图像，从而达到确定地层结构或圈定地质异常体的一种物探反演解释方法。基于射线理论的射线追踪法是正演算法的一种，射线追踪法基于高频近似，只计算源点到接收点的射线路径，计算效率高。反演计算一般是采用最优化算法，例如最速下降法、共轭梯度法、牛顿迭代法、随机搜索等，求解满足给定最小误差准则的最优解。
当地层速度只在深度方向上变化，水平方向上无变化时，地震层析成像可用一维速度模型描述研究区域的速度结构，即速度只是深度的函数。在近地表浅层的层析成像问题中，小区域范围内的地层速度变化一般可用一维模型近似，在地震数据采集密度不高的情况下，工程物探一般会采用基于一维模型的地震层析成像方法。现有的一维走时层析速度反演假设每层层内速度都是均匀的，通过多次迭代来得到一个层状的速度模型，一般需要划分比较多的层数才能比较准确地描述实际地层速度变化的特征。实际地层的速度变化多为非均匀且呈现出连续变化的特征，因此利用现有技术得到的结果存在固有近似性。并且层数越多需要反演的模型参数就越多，因而地震层析成像的计算量越大。另 外，在大入射角的情况下(一般震源与接收点的距离比地震信号反射深度大很多的时候会出现此情况)，利用现有技术中的两点走时射线追踪牛顿迭代方法存在收敛慢或不收敛的问题。
技术问题
本发明的目的在于提供一种基于两点射线追踪的地震走时层析反演方法，旨在解决现有技术中地震层析成像的计算量繁重和现有的两点走时射线追踪牛顿迭代方法存在的收敛慢或不收敛的问题。
技术解决方案
为解决上述技术问题，本发明是这样实现的，一种基于两点射线追踪的地震走时层析反演方法，包括下述步骤:
在研究区域进行地震数据采集，获取直达波走时数据和反射波走时数据；
对研究区域进行模型参数化，建立层内速度可连续变化的初始一维连续层状模型；
根据所述层内速度可连续变化的一维连续层状模型，将射线参数p用变量q表示，震源距X表示为变量q的函数X＝f(q)，利用牛顿迭代法求解X＝f(q)，进而得到射线参数p，射线路径由射线参数p唯一确定，得到射线参数p后，计算得到理论的直达波走时和反射波走时；
比较计算得到的理论直达波走时和反射波走时与所述地震数据采集获得的实际的直达波走时数据和反射波走时数据，判断理论直达波走时和反射波走时与所述地震数据采集获得的实际的直达波走时数据和反射波走时数据的差别是否满足给定的误差标准，是则进行模型输出，否则进行下一步；
用最优化算法调整所述层内速度可连续变化的一维连续层状模型，直到计算得到的理论的直达波走时和反射波走时和实际的直达波走时数据和反射波 走时数据的差别满足给定的误差标准为止，并进行模型输出。
有益效果
与现有技术相比，本发明提供的基于两点射线追踪的地震走时层析反演方法与现有技术相比存在以下有益效果：通过建立地层速度可连续变化的一维连续层状模型，使基于此模型的反演变量数可以大幅度减少，可以更加准确地描述实际地层速度结构，显著提高反演计算效率；通过将射线参数p用变量q表示，间接利用变量q求解射线参数p，使得迭代求解过程稳定，收敛快速，有效避免在大入射角的情况下迭代不收敛的问题。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1是本发明实施例提供的基于两点射线追踪的地震走时层析反演方法流程图；
图2是本发明实施例提供的基于两点射线追踪的地震走时层析反演方法连续层状模型的参数定义；
图3是本发明实施例提供的基于两点射线追踪的地震走时层析反演方法地层模型示意图；
图4是本发明实施例提供的基于两点射线追踪的地震走时层析反演方法地震勘探地表数据采集示意图；
图5是本发明实施例提供的基于两点射线追踪的地震走时层析反演方法地表激发垂直地震剖面示意图；
图6是本发明实施例提供的基于两点射线追踪的地震走时层析反演方法井中激发垂直地震剖面示意图；
图7是本发明实施例提供的基于两点射线追踪的地震走时层析反演方法井间成像示意图。
本发明的实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
请参阅图1至图3，本发明实施例提供了一种基于两点射线追踪的地震走时层析反演方法，该方法包括下述步骤：
步骤S1：采集实际走时数据。在研究区域进行地震数据采集，获取直达波走时数据和反射波走时数据；
步骤S2：模型参数化。对研究区域进行模型参数化，建立层内速度可连续变化的初始一维连续层状模型；
具体地，该层内速度可以是纵波或横波速度，随深度递增或递减；
步骤S3：计算射线参数。根据层内速度可连续变化的一维连续层状模型，将射线参数p用变量q表示，震源距X表示为变量q的函数X＝f(q)，利用牛顿迭代法求解X＝f(q)，进而得到射线参数p，射线路径由射线参数p唯一确定；
步骤S4：计算理论走时。根据射线参数p计算直达波走时和反射波走时；
步骤S5：比较理论走时与实际走时。比较计算得到的理论直达波走时和反射波走时与所述地震数据采集获得的实际的直达波走时数据和反射波走时数据，判断理论直达波走时和反射波走时与所述地震数据采集获得的实际的直 达波走时数据和反射波走时数据的差别是否满足给定的误差标准，是则进行步骤S7，否则进行步骤S6；
步骤S6：优化模型。用最优化算法调整层内速度可连续变化的一维连续层状模型，直到计算得到的理论的直达波走时和反射波走时和实际的直达波走时数据和反射波走时数据的差别满足给定的误差标准为止。
具体地，使用的最优化算法可以为最速下降法、共轭梯度法、牛顿迭代法、随机搜索等方法；
步骤S7：模型输出。计算得到的理论到时和实际数据到时的差别满足给定的误差标准后，进行模型输出。
本发明实施例提供的基于两点射线追踪的地震走时层析反演方法建立了层内速度可连续变化的一维连续层状模型，将层内速度表示为深度的函数，此模型允许层内速度连续变化，大大减少了划分层数，从而大幅度减少了反演变量数，可以更加准确地描述实际地层速度结构，显著提高反演计算效率；并且通过将射线参数p用变量q表示，间接利用变量q求解射线参数p，使得迭代求解过程稳定，收敛快速，有效避免了在大入射角的情况下迭代不收敛的问题。
进一步地，请参考图2，在层内速度可连续变化的一维连续层状模型中，层内速度Vk为深度z的函数，第k层速度函数表示为：
Vk＝akz+bk，
其中，下标k表示第k层，ak和bk是需要反演的模型参数，分别为第k层速度函数的梯度和截距，当反演地层纵波速度模型时，层内速度Vk为纵波速度；当反演地层横波速度模型时，则层内速度Vk为横波速度；当反演地层 转换波速度模型时，则层内速度Vk由层内转换波属性决定其为纵波速度或横波速度，这里的转换波为纵波到横波转换波或横波到纵波转换波，该层内转换波属性为纵波属性或横波属性。
当ak＝0时，表示在第k层为层内速度不变的均匀层。
进一步地，根据斯涅尔定律，震源距X表示为射线参数p的函数，震源距X表示为：
Figure PCTCN2017105817-appb-000001
其中，各项定义为：
Figure PCTCN2017105817-appb-000002
Figure PCTCN2017105817-appb-000003
Figure PCTCN2017105817-appb-000004
Figure PCTCN2017105817-appb-000005
Figure PCTCN2017105817-appb-000006
其中，εk，ωk，hk，μk和δk为中间参数，下标s表示震源所在的层，可取值范围为1至n，n为反射波反射发生所在层的标号，zs为震源深度，z(k)表示第k层的深度，下标k表示第k层，k＝0项是关于震源位置的一个修正项。
进一步地，为了避免在大入射角情况下迭代不收敛的问题，将射线参数p用变量q表示为：
Figure PCTCN2017105817-appb-000007
其中，VM为模拟射线路径经过地层的最大速度。
进一步地，将震源距X表示为变量q的函数，震源距X表示为：
Figure PCTCN2017105817-appb-000008
其中，各项定义为：
Figure PCTCN2017105817-appb-000009
Figure PCTCN2017105817-appb-000010
Figure PCTCN2017105817-appb-000011
Figure PCTCN2017105817-appb-000012
给定震源距X，利用牛顿迭代法求解X＝f(q)，得到参数q的值，将参数q代回到与射线参数p的关系式中，即可得到射线参数p的值。
进一步地，牛顿迭代法求解方程X＝f(q)时，q的初值通过以下方法获得：将震源距X和变量q用有理函数公式逼近，有理函数公式为：
Figure PCTCN2017105817-appb-000013
其中，系数α1，α2，β1和β2通过有理函数公式和震源X的函数公式的泰勒展开式获得，利用上述有理函数公式，获得q的初值估算公式：
Figure PCTCN2017105817-appb-000014
其中，各项系数α1，α2，β1和β2的表达式分别为：
Figure PCTCN2017105817-appb-000015
Figure PCTCN2017105817-appb-000016
Figure PCTCN2017105817-appb-000017
Figure PCTCN2017105817-appb-000018
Figure PCTCN2017105817-appb-000019
Figure PCTCN2017105817-appb-000020
Figure PCTCN2017105817-appb-000021
Figure PCTCN2017105817-appb-000022
Figure PCTCN2017105817-appb-000023
Figure PCTCN2017105817-appb-000024
利用所述q的初值估算公式获得初值进行迭代计算，经过迭代得到q的精确解。利用此方法获得的初值精度一般能达到95％以上，只需要1至3次迭代就可以得到q的精确解。
进一步地，得到射线参数p的值后，直达波走时和反射波走时的计算公式为：
Figure PCTCN2017105817-appb-000025
进一步地，请参考图4至图7，在研究区域进行地震数据采集的步骤具体为：
如图4所示的地表地震勘探，震源和接收检波器都布置在地表；
或者为：如图5所示的垂直地震剖面勘探，震源在地表，接收器在井中；
或者为：如图6所示的垂直地震剖面勘探，震源在井中，接收器在地表；
或者为：如图7所示的井间成像勘探，震源和接收检波器在两口不同的井中。
进一步地，地震数据采集采用的地震信号是纵波、或横波、或纵波到横波的转换波，或横波到纵波的转换波。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (10)

  1. 一种基于两点射线追踪的地震走时层析反演方法，其特征在于，包括下述步骤:
    在研究区域进行地震数据采集，获取直达波走时数据和反射波走时数据；
    对研究区域进行模型参数化，建立层内速度可连续变化的初始一维连续层状模型；
    根据所述层内速度可连续变化的一维连续层状模型，将射线参数p用变量q表示，震源距X表示为变量q的函数X＝f(q)，利用牛顿迭代法求解X＝f(q)，进而得到射线参数p，射线路径由射线参数p唯一确定，得到射线参数p后，计算得到理论的直达波走时和反射波走时；
    比较计算得到的理论直达波走时和反射波走时与所述地震数据采集获得的实际的直达波走时数据和反射波走时数据，判断理论直达波走时和反射波走时与所述地震数据采集获得的实际的直达波走时数据和反射波走时数据的差别是否满足给定的误差标准，是则进行模型输出，否则进行下一步；
    用最优化算法调整所述层内速度可连续变化的一维连续层状模型，直到计算得到的理论的直达波走时和反射波走时和实际的直达波走时数据和反射波走时数据的差别满足给定的误差标准为止，并进行模型输出。
  2. 根据权利要求1所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，在所述层内速度可连续变化的一维连续层状模型中，层内速度Vk为深度z的函数，第k层速度函数表示为：
    Vk＝akz+bk，
    其中，下标k表示第k层，ak和bk是需要反演的模型参数，分别为所述第k层速度函数的梯度和截距，当反演地层纵波速度模型时，所述层内速度Vk为纵波速度；当反演地层横波速度模型时，则所述层内速度Vk为横波速度；当反演地层转换波速度模型时，则所述层内速度Vk由层内转换波属性决定其为纵波速度或横波速度，所述转换波为纵波到横波转换波或横波到纵波转换波，所述层内转换波属性为纵波属性或横波属性。
  3. 根据权利要求2所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，当ak＝0时，所述层内速度可连续变化的一维连续层状模型在第k层为层内速度不变的均匀层。
  4. 根据权利要求1所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，根据斯涅尔定律，震源距X表示为射线参数p的函数，震源距X表示为：
    Figure PCTCN2017105817-appb-100001
    其中，各项定义为：
    Figure PCTCN2017105817-appb-100002
    Figure PCTCN2017105817-appb-100003
    Figure PCTCN2017105817-appb-100004
    Figure PCTCN2017105817-appb-100005
    Figure PCTCN2017105817-appb-100006
    其中，εk，ωk，hk，μk和δk为中间参数，下标s表示震源所在的层，可取值范围为1至n，n为反射波反射发生所在层的标号，zs为震源深度，z(k)表示第k层的深度，下标k表示第k层，k＝0项是关于震源位置的一个修正项。
  5. 根据权利要求4所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，所述射线参数p用变量q表示为：
    Figure PCTCN2017105817-appb-100007
    其中，VM为模拟射线路径经过地层的最大速度。
  6. 根据权利要求5所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，将所述震源距X表示为所述变量q的函数，震源距X表示为：
    Figure PCTCN2017105817-appb-100008
    其中，各项定义为：
    Figure PCTCN2017105817-appb-100009
    Figure PCTCN2017105817-appb-100010
    Figure PCTCN2017105817-appb-100011
    Figure PCTCN2017105817-appb-100012
    给定震源距X，利用牛顿迭代法求解X＝f(q)，得到所述参数q的值，将所述参数q代回到与射线参数p的关系式
    Figure PCTCN2017105817-appb-100013
    中，即可得到射线参数p的值。
  7. 根据权利要求6所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，牛顿迭代法求解方程X＝f(q)时，q的初值通过以下方法获得：将所述震源距X和所述变量q用有理函数公式逼近，所述有理函数公式为：
    Figure PCTCN2017105817-appb-100014
    其中，系数α1，α2，β1和β2通过对比所述有理函数公式和所述震源距X的函数公式的泰勒展开式获得，利用所述有理函数公式，获得q的初值估算公式：
    Figure PCTCN2017105817-appb-100015
    其中，各项系数α1，α2，β1和β2的表达式分别为：
    Figure PCTCN2017105817-appb-100016
    Figure PCTCN2017105817-appb-100017
    Figure PCTCN2017105817-appb-100018
    Figure PCTCN2017105817-appb-100019
    Figure PCTCN2017105817-appb-100020
    Figure PCTCN2017105817-appb-100021
    Figure PCTCN2017105817-appb-100022
    Figure PCTCN2017105817-appb-100023
    Figure PCTCN2017105817-appb-100024
    Figure PCTCN2017105817-appb-100025
    利用所述q的初值估算公式获得初值进行迭代计算，经过迭代得到q的精确解。
  8. 根据权利要求1至7所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，所述直达波走时和反射波走时的计算公式为：
    Figure PCTCN2017105817-appb-100026
  9. 根据权利要求1所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，所述在研究区域进行地震数据采集的步骤具体为：
    地表地震勘探，震源和接收检波器都布置在地表；
    或者为：垂直地震剖面，震源在地表，接收器在井中；
    或者为：垂直地震剖面，震源在井中，接收器在地表；
    或者为：井间成像，震源和接收检波器在两口不同的井中。
  10. 根据权利要求1所述的一种基于两点射线追踪的地震走时层析反演方法，其特征在于，所述地震数据采集采用的地震信号是纵波、或横波、或纵波到横波的转换波，或横波到纵波的转换波。
PCT/CN2017/105817 2017-10-12 2017-10-12 一种基于两点射线追踪的地震走时层析反演方法 Ceased WO2019071504A1 (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201780001180.7A CN108064348B (zh) 2017-10-12 2017-10-12 一种基于两点射线追踪的地震走时层析反演方法
PCT/CN2017/105817 WO2019071504A1 (zh) 2017-10-12 2017-10-12 一种基于两点射线追踪的地震走时层析反演方法
US15/950,350 US20190113641A1 (en) 2017-10-12 2018-04-11 Seismic travel time tomographic inversion method based on two point ray tracing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2017/105817 WO2019071504A1 (zh) 2017-10-12 2017-10-12 一种基于两点射线追踪的地震走时层析反演方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US15/950,350 Continuation-In-Part US20190113641A1 (en) 2017-10-12 2018-04-11 Seismic travel time tomographic inversion method based on two point ray tracing

Publications (1)

Publication Number Publication Date
WO2019071504A1 true WO2019071504A1 (zh) 2019-04-18

Family

ID=62141990

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2017/105817 Ceased WO2019071504A1 (zh) 2017-10-12 2017-10-12 一种基于两点射线追踪的地震走时层析反演方法

Country Status (3)

Country Link
US (1) US20190113641A1 (zh)
CN (1) CN108064348B (zh)
WO (1) WO2019071504A1 (zh)

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110618460A (zh) * 2019-07-22 2019-12-27 中国石油化工股份有限公司 一种结合层位信息的微测井方位加权插值建模方法
CN111580157A (zh) * 2020-06-08 2020-08-25 石川泰克(北京)能源有限公司 一种叠前深度偏移的近似真地表速度模型建立的方法
CN111650638A (zh) * 2020-05-21 2020-09-11 长江大学 一种地震波旅行时间计算方法
CN112083486A (zh) * 2019-06-14 2020-12-15 中国石油天然气集团有限公司 一种低速层速度获取方法及装置
CN112180441A (zh) * 2019-07-03 2021-01-05 中国石油天然气集团有限公司 转换波初始速度建模方法及装置
CN112305595A (zh) * 2019-07-24 2021-02-02 中国石油化工股份有限公司 基于折射波分析地质体结构的方法及存储介质
CN112485825A (zh) * 2019-09-11 2021-03-12 中国石油化工股份有限公司 一种基于初至波走时层析的微测井解释方法
CN112526610A (zh) * 2019-09-17 2021-03-19 中国石油化工股份有限公司 一种约束表层建模的三维地震采集激发井深设计方法
CN112596103A (zh) * 2020-11-24 2021-04-02 中国地质科学院地球物理地球化学勘查研究所 射线追踪方法、装置和电子设备
CN113589375A (zh) * 2020-04-30 2021-11-02 中国石油化工股份有限公司 基于倾斜层约束旅行时计算的vsp层速度反演方法
CN113777654A (zh) * 2021-08-06 2021-12-10 同济大学 一种基于伴随状态法初至波走时层析的海水速度建模方法
CN113805232A (zh) * 2020-06-17 2021-12-17 中国石油化工股份有限公司 浅层地表的品质因子的估计方法、系统及存储介质
CN113885076A (zh) * 2021-09-30 2022-01-04 吉林大学 一种微地震地面监测速度模型校正方法
CN114428318A (zh) * 2020-10-15 2022-05-03 中国石油化工股份有限公司 一种初至层析近地表速度高效建模方法及系统
CN114814949A (zh) * 2021-01-21 2022-07-29 中国石油化工股份有限公司 一种浅层逆vsp初至层析及地层预测方法
CN114839675A (zh) * 2021-01-31 2022-08-02 中国石油化工股份有限公司 一种建立三维速度模型的方法
CN115308801A (zh) * 2022-08-29 2022-11-08 南方海洋科学与工程广东省实验室(广州) 一种利用直达波走时和地形数据定位海底地震仪的方法及处理终端
CN115358087A (zh) * 2022-09-06 2022-11-18 中国地质科学院地球物理地球化学勘查研究所 基于多尺度、多方向逐段迭代和反向追踪的联合射线追踪方法
CN115755183A (zh) * 2022-11-15 2023-03-07 哈尔滨工程大学 一种基于尾波干涉的震源位置追踪方法
CN116009067A (zh) * 2021-10-21 2023-04-25 中国石油化工股份有限公司 自适应阻尼正则化的近地表速度建模方法及装置
CN116009072A (zh) * 2021-10-21 2023-04-25 中国石油化工股份有限公司 基于射线密度均衡的小尺度异常体层析反演方法和装置
CN116340710A (zh) * 2023-05-30 2023-06-27 中国科学院精密测量科学与技术创新研究院 基于分层快速三维射线追踪的中性大气斜延迟计算方法
CN116879950A (zh) * 2023-07-12 2023-10-13 成都理工大学 基于直达波和sPL初动极性与振幅比的震源机制反演方法
CN117784223A (zh) * 2022-09-21 2024-03-29 新奥(天津)能源技术有限公司 射线追踪方法、装置、设备及计算机可读存储介质
CN118244355A (zh) * 2024-05-30 2024-06-25 山东省科学院海洋仪器仪表研究所 一种基于重建观测地震数据的反射波走时反演方法
CN118330732A (zh) * 2024-04-30 2024-07-12 中国地震局地球物理研究所 一种基于三维tti介质模型的地震定位方法
CN118363062A (zh) * 2024-04-26 2024-07-19 中海石油(中国)有限公司海南分公司 一种高精度射线参数域纵横波联合叠前反演方法
CN119716987A (zh) * 2024-12-24 2025-03-28 西安交通大学 一种走时层析成像方法及系统
CN119781035A (zh) * 2023-10-08 2025-04-08 中国石油化工股份有限公司 基于节点数据的三角网格层析速度反演优化方法
CN119783424A (zh) * 2023-10-08 2025-04-08 中国石油化工股份有限公司 基于海水中波传播射线理论模拟的水体体波现象识别方法
CN119960018A (zh) * 2023-11-08 2025-05-09 中国石油天然气集团有限公司 一种近地表网格速度模型的建立方法、装置和电子设备
CN120028831A (zh) * 2023-11-21 2025-05-23 中国石油天然气集团有限公司 一种射线追踪区域分析方法和装置

Families Citing this family (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR112018070565A2 (pt) 2016-04-07 2019-02-12 Bp Exploration Operating Company Limited detecção de eventos de fundo de poço usando características de domínio da frequência acústicas
AU2017246520B2 (en) 2016-04-07 2022-04-07 Bp Exploration Operating Company Limited Detecting downhole events using acoustic frequency domain features
EP3608503B1 (en) 2017-03-31 2022-05-04 BP Exploration Operating Company Limited Well and overburden monitoring using distributed acoustic sensors
BR112020003742A2 (pt) 2017-08-23 2020-09-01 Bp Exploration Operating Company Limited detecção de localizações de ingresso de areia em fundo de poço
CA3078842C (en) 2017-10-11 2024-01-09 Bp Exploration Operating Company Limited Detecting events using acoustic frequency domain features
CN111164462B (zh) * 2018-08-06 2022-05-06 南方科技大学 一种人工源面波勘探方法、面波勘探装置及终端设备
US20200174149A1 (en) 2018-11-29 2020-06-04 Bp Exploration Operating Company Limited Event Detection Using DAS Features with Machine Learning
GB201820331D0 (en) 2018-12-13 2019-01-30 Bp Exploration Operating Co Ltd Distributed acoustic sensing autocalibration
CN110187382B (zh) * 2019-03-05 2020-10-13 中国石油大学(华东) 一种回折波和反射波波动方程旅行时反演方法
WO2021052602A1 (en) 2019-09-20 2021-03-25 Lytt Limited Systems and methods for sand ingress prediction for subterranean wellbores
CN110568496B (zh) * 2019-09-26 2021-02-09 核工业北京地质研究院 一种复杂介质条件下射线追踪方法
EP4045766A1 (en) 2019-10-17 2022-08-24 Lytt Limited Fluid inflow characterization using hybrid das/dts measurements
CA3154435C (en) 2019-10-17 2023-03-28 Lytt Limited Inflow detection using dts features
CN110879412A (zh) * 2019-10-31 2020-03-13 南方科技大学 地下横波速度反演方法、装置、计算设备及存储介质
WO2021093974A1 (en) 2019-11-15 2021-05-20 Lytt Limited Systems and methods for draw down improvements across wellbores
CN113156495B (zh) * 2020-01-07 2024-06-25 中国石油天然气集团有限公司 网格层析反演反射点确定方法及装置
CN111273344B (zh) * 2020-03-02 2022-01-25 广州海洋地质调查局 一种基于续至折射波的层析反演方法及处理终端
CN113534250A (zh) * 2020-04-18 2021-10-22 中国石油化工股份有限公司 一种基于快速匹配追踪的多尺度地震反演方法
WO2021249643A1 (en) 2020-06-11 2021-12-16 Lytt Limited Systems and methods for subterranean fluid flow characterization
CA3182376A1 (en) 2020-06-18 2021-12-23 Cagri CERRAHOGLU Event model training using in situ data
WO2021254633A1 (en) 2020-06-18 2021-12-23 Lytt Limited Event model training using in situ data
CN113970789B (zh) * 2020-07-24 2024-04-09 中国石油化工股份有限公司 全波形反演方法、装置、存储介质及电子设备
CN112068185B (zh) * 2020-08-24 2022-09-13 东南大学 一种融合球谐函数和近似Chapman函数的电离层层析方法
CN112255671A (zh) * 2020-08-28 2021-01-22 长江大学 一种两点间地震波正演方法及装置
CN113761462B (zh) * 2021-09-10 2022-05-31 山东大学 基于弦截法的初始入射角迭代计算改进方法
CN113791447B (zh) * 2021-10-12 2023-06-20 同济大学 一种反射结构导引的反射波层析反演方法
CN115993632B (zh) * 2021-10-20 2025-12-30 中国石油化工股份有限公司 一种优化关键采集参数的方法及系统
CN114047549B (zh) * 2021-10-28 2025-01-24 中国石油化工股份有限公司 一种用于多工区的融合层析静校正处理方法
CN114966831B (zh) * 2022-01-28 2025-01-21 东北石油大学 一种基于速度-衰减解耦的粘声全波形反演方法
CN114879249B (zh) * 2022-04-13 2023-04-28 中国海洋大学 基于四面体单元走时扰动插值的地震波前走时计算方法
CN115032694B (zh) * 2022-04-14 2025-01-24 海南浙江大学研究院 一种vsp初至旅行时层析成像方法及系统
CN117518261B (zh) * 2022-07-29 2026-03-10 中国石油天然气集团有限公司 地震数据处理方法和装置
CN116009085B (zh) * 2023-02-02 2024-03-12 哈尔滨工业大学 一种基于全波形反演的软地层横波速度测量方法及装置
CN119620170A (zh) * 2023-09-12 2025-03-14 中国石油天然气股份有限公司 初至波走时层析速度反演方法、系统、设备和介质
CN119805566B (zh) * 2023-10-09 2025-11-11 中国石油天然气集团有限公司 地震数据处理方法、装置、设备及存储介质
CN118011483B (zh) * 2024-02-04 2024-12-13 湖南工商大学 一种考虑地形效应的面波层析成像方法及系统
CN117724166A (zh) * 2024-02-07 2024-03-19 中国石油大学(华东) 基于大炮初至的近地表三维速度建模方法
CN118095666B (zh) * 2024-04-29 2024-06-21 山东科岳科技有限公司 一种台网监测能力评估方法
CN119535561B (zh) * 2024-11-26 2025-09-23 兖矿能源集团股份有限公司 一种基于主被动双源联合震动波ct反演方法及冲击地压预警方法
CN119903583A (zh) * 2024-12-31 2025-04-29 江汉大学 基坑爆破影响下围护桩基保护减振孔设计方法
CN119781006B (zh) * 2025-01-24 2025-10-31 广东省水利电力勘测设计研究院有限公司 基于跨孔地震波ct检查技术的地下管线探测系统及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102841376A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种基于起伏地表的层析速度反演方法
CN105425286A (zh) * 2015-10-30 2016-03-23 中国石油天然气集团公司 地震走时获取方法及基于其的井间地震走时层析成像方法
CN105549081A (zh) * 2016-01-29 2016-05-04 中国石油大学(华东) 各向异性介质共炮域高斯束偏移成像方法
CN105589100A (zh) * 2014-10-21 2016-05-18 中国石油化工股份有限公司 一种微地震震源位置和速度模型同时反演方法
CN106353793A (zh) * 2015-07-17 2017-01-25 中国石油化工股份有限公司 一种基于走时增量双线性插值射线追踪的井间地震层析反演方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6278948B1 (en) * 1999-04-02 2001-08-21 Conoco Inc. Method for gravity and magnetic data inversion using vector and tensor data
US7663972B2 (en) * 2008-02-22 2010-02-16 Pgs Geophysical As Method for three dimensional seismic travel time tomography in transversely isotropic media
US8553497B2 (en) * 2008-08-11 2013-10-08 Exxonmobil Upstream Research Company Removal of surface-wave noise in seismic data
US8861309B2 (en) * 2011-01-31 2014-10-14 Chevron U.S.A. Inc. Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context
WO2014111431A2 (en) * 2013-01-15 2014-07-24 Cgg Services Sa System and method for ray based tomography guided by waveform inversion

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102841376A (zh) * 2012-09-06 2012-12-26 中国石油大学(华东) 一种基于起伏地表的层析速度反演方法
CN105589100A (zh) * 2014-10-21 2016-05-18 中国石油化工股份有限公司 一种微地震震源位置和速度模型同时反演方法
CN106353793A (zh) * 2015-07-17 2017-01-25 中国石油化工股份有限公司 一种基于走时增量双线性插值射线追踪的井间地震层析反演方法
CN105425286A (zh) * 2015-10-30 2016-03-23 中国石油天然气集团公司 地震走时获取方法及基于其的井间地震走时层析成像方法
CN105549081A (zh) * 2016-01-29 2016-05-04 中国石油大学(华东) 各向异性介质共炮域高斯束偏移成像方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
TIAN, YUE ET AL.: "A Rapid and Accurate Two-Point Ray Tracing Method in Horizontally Layered Velocity Model", ACTA SEISMOLOGICA SINICA, vol. 2, no. 27, 31 March 2015 (2015-03-31), pages 147 - 154, ISSN: 0253-3782 *

Cited By (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112083486A (zh) * 2019-06-14 2020-12-15 中国石油天然气集团有限公司 一种低速层速度获取方法及装置
CN112180441A (zh) * 2019-07-03 2021-01-05 中国石油天然气集团有限公司 转换波初始速度建模方法及装置
CN112180441B (zh) * 2019-07-03 2024-03-26 中国石油天然气集团有限公司 转换波初始速度建模方法及装置
CN110618460A (zh) * 2019-07-22 2019-12-27 中国石油化工股份有限公司 一种结合层位信息的微测井方位加权插值建模方法
CN112305595A (zh) * 2019-07-24 2021-02-02 中国石油化工股份有限公司 基于折射波分析地质体结构的方法及存储介质
CN112305595B (zh) * 2019-07-24 2024-05-17 中国石油化工股份有限公司 基于折射波分析地质体结构的方法及存储介质
CN112485825B (zh) * 2019-09-11 2024-04-09 中国石油化工股份有限公司 一种基于初至波走时层析的微测井解释方法
CN112485825A (zh) * 2019-09-11 2021-03-12 中国石油化工股份有限公司 一种基于初至波走时层析的微测井解释方法
CN112526610A (zh) * 2019-09-17 2021-03-19 中国石油化工股份有限公司 一种约束表层建模的三维地震采集激发井深设计方法
CN113589375A (zh) * 2020-04-30 2021-11-02 中国石油化工股份有限公司 基于倾斜层约束旅行时计算的vsp层速度反演方法
CN111650638A (zh) * 2020-05-21 2020-09-11 长江大学 一种地震波旅行时间计算方法
CN111580157A (zh) * 2020-06-08 2020-08-25 石川泰克(北京)能源有限公司 一种叠前深度偏移的近似真地表速度模型建立的方法
CN113805232A (zh) * 2020-06-17 2021-12-17 中国石油化工股份有限公司 浅层地表的品质因子的估计方法、系统及存储介质
CN113805232B (zh) * 2020-06-17 2024-04-09 中国石油化工股份有限公司 浅层地表的品质因子的估计方法、系统及存储介质
CN114428318A (zh) * 2020-10-15 2022-05-03 中国石油化工股份有限公司 一种初至层析近地表速度高效建模方法及系统
CN112596103A (zh) * 2020-11-24 2021-04-02 中国地质科学院地球物理地球化学勘查研究所 射线追踪方法、装置和电子设备
CN114814949A (zh) * 2021-01-21 2022-07-29 中国石油化工股份有限公司 一种浅层逆vsp初至层析及地层预测方法
CN114814949B (zh) * 2021-01-21 2023-09-01 中国石油化工股份有限公司 一种浅层逆vsp初至层析及地层预测方法
CN114839675A (zh) * 2021-01-31 2022-08-02 中国石油化工股份有限公司 一种建立三维速度模型的方法
CN114839675B (zh) * 2021-01-31 2023-09-05 中国石油化工股份有限公司 一种建立三维速度模型的方法
CN113777654A (zh) * 2021-08-06 2021-12-10 同济大学 一种基于伴随状态法初至波走时层析的海水速度建模方法
CN113885076A (zh) * 2021-09-30 2022-01-04 吉林大学 一种微地震地面监测速度模型校正方法
CN116009072A (zh) * 2021-10-21 2023-04-25 中国石油化工股份有限公司 基于射线密度均衡的小尺度异常体层析反演方法和装置
CN116009067A (zh) * 2021-10-21 2023-04-25 中国石油化工股份有限公司 自适应阻尼正则化的近地表速度建模方法及装置
CN115308801A (zh) * 2022-08-29 2022-11-08 南方海洋科学与工程广东省实验室(广州) 一种利用直达波走时和地形数据定位海底地震仪的方法及处理终端
CN115358087A (zh) * 2022-09-06 2022-11-18 中国地质科学院地球物理地球化学勘查研究所 基于多尺度、多方向逐段迭代和反向追踪的联合射线追踪方法
CN117784223A (zh) * 2022-09-21 2024-03-29 新奥(天津)能源技术有限公司 射线追踪方法、装置、设备及计算机可读存储介质
CN115755183A (zh) * 2022-11-15 2023-03-07 哈尔滨工程大学 一种基于尾波干涉的震源位置追踪方法
CN116340710B (zh) * 2023-05-30 2023-09-12 中国科学院精密测量科学与技术创新研究院 基于分层快速三维射线追踪的中性大气斜延迟计算方法
CN116340710A (zh) * 2023-05-30 2023-06-27 中国科学院精密测量科学与技术创新研究院 基于分层快速三维射线追踪的中性大气斜延迟计算方法
CN116879950B (zh) * 2023-07-12 2024-03-08 成都理工大学 基于直达波和sPL初动极性与振幅比的震源机制反演方法
CN116879950A (zh) * 2023-07-12 2023-10-13 成都理工大学 基于直达波和sPL初动极性与振幅比的震源机制反演方法
CN119781035A (zh) * 2023-10-08 2025-04-08 中国石油化工股份有限公司 基于节点数据的三角网格层析速度反演优化方法
CN119783424A (zh) * 2023-10-08 2025-04-08 中国石油化工股份有限公司 基于海水中波传播射线理论模拟的水体体波现象识别方法
CN119781035B (zh) * 2023-10-08 2025-12-09 中国石油化工股份有限公司 基于节点数据的三角网格层析速度反演优化方法
CN119960018A (zh) * 2023-11-08 2025-05-09 中国石油天然气集团有限公司 一种近地表网格速度模型的建立方法、装置和电子设备
CN120028831A (zh) * 2023-11-21 2025-05-23 中国石油天然气集团有限公司 一种射线追踪区域分析方法和装置
CN118363062A (zh) * 2024-04-26 2024-07-19 中海石油(中国)有限公司海南分公司 一种高精度射线参数域纵横波联合叠前反演方法
CN118330732A (zh) * 2024-04-30 2024-07-12 中国地震局地球物理研究所 一种基于三维tti介质模型的地震定位方法
CN118244355A (zh) * 2024-05-30 2024-06-25 山东省科学院海洋仪器仪表研究所 一种基于重建观测地震数据的反射波走时反演方法
CN119716987A (zh) * 2024-12-24 2025-03-28 西安交通大学 一种走时层析成像方法及系统

Also Published As

Publication number Publication date
CN108064348A (zh) 2018-05-22
CN108064348B (zh) 2020-05-05
US20190113641A1 (en) 2019-04-18

Similar Documents

Publication Publication Date Title
WO2019071504A1 (zh) 一种基于两点射线追踪的地震走时层析反演方法
CN104133245B (zh) 一种地震资料的静校正方法及系统
KR102021276B1 (ko) 진폭 보존을 갖는 fwi 모델 도메인 각도 스택들
CN105277978B (zh) 一种确定近地表速度模型的方法及装置
CN108919354B (zh) 近地表q偏移方法及装置
CN108363101B (zh) 一种斜井井间地震高斯束叠前深度偏移成像方法
CN110007340B (zh) 基于角度域直接包络反演的盐丘速度密度估计方法
CN111123359B (zh) 随钻测井与地层格架约束的井周地震成像探测方法及装置
CN109444956B (zh) 三维起伏观测面地震斜率层析成像方法
CN107462924A (zh) 一种不依赖于测井资料的绝对波阻抗反演方法
CN109633749A (zh) 基于散射积分法的非线性菲涅尔体地震走时层析成像方法
CN102901984B (zh) 真地表地震数据倾角道集构建方法
CN115903042A (zh) 一种基于构造整形正则化的波形反演方法及其设备
CN113466933B (zh) 基于深度加权的地震斜率层析成像方法
CN109655890A (zh) 一种深度域浅中深层联合层析反演速度建模方法及系统
CN104199088B (zh) 一种提取入射角道集的方法及系统
CN115032694B (zh) 一种vsp初至旅行时层析成像方法及系统
CN117991377B (zh) 一种基于多源信息融合的初至波走时层析成像方法及系统
CN115685329A (zh) 一种初始速度模型的构建方法、装置、电子设备和介质
CN119939886A (zh) 一种层析速度建模方法、装置及设备
Zeng et al. Subsalt high steep structure imaging technique: A case study of Keshen area in Kuqa depression, Tarim Basin, NW China
CN109085644A (zh) 基于双射线走时的真地表成像方法
CN108254787A (zh) 波的走时获得方法及装置、成像方法及装置
Balkaya et al. A comparison of two travel-time tomography schemes for crosshole radar data: Eikonal-equation-based inversion versus ray-based inversion
Jingcheng et al. An Integrated Refraction and Tomography Static Correction Method and its Application in Sichuan Changning, China

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: 17928766

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: 17928766

Country of ref document: EP

Kind code of ref document: A1