公开(公告)号：US10459096B2

公开(公告)日：2019-10-29

申请号：US15614960

申请日：2017-06-06

Applicant: Sirui Tan ,  Yaxun Tang ,  Anatoly I. Baumstein ,  Gboyega Ayeni ,  Tetyana Vdovina ,  Thomas A. Dickens

Inventor： Sirui Tan ,  Yaxun Tang ,  Anatoly I. Baumstein ,  Gboyega Ayeni ,  Tetyana Vdovina ,  Thomas A. Dickens

IPC: G01V1/28 ,  G01V1/30 ,  G01V1/36 ,  G06F17/11

Abstract: A method for iteratively inverting seismic data to jointly infer a model for at least P-wave velocity and attenuation parameters of the subsurface, the method including: jointly inverting the P-wave velocity and attenuation parameters with an iterative visco-acoustic full wavefield inversion process, wherein the iterative visco-acoustic full wavefield inversion process includes computing a gradient of an objective function, the objective function measuring a misfit between all or part of the seismic data and corresponding model-simulated seismic data; for each of the P-wave velocity and attenuation parameters, computing a search direction in model space from the gradient; determining line search step sizes α and β for the search directions for the P-wave velocity and attenuation parameters, respectively, wherein a ratio of the step sizes is a function of the P-wave velocity parameter; and using the step sizes α and β and the search directions for each of the P-wave velocity and attenuation parameters, computing a new search direction in model space, then performing a line search along the new search direction to arrive at a new step size, and using the new step size and the new search direction to generate an updated model for a current iteration of the iterative visco-acoustic full wavefield inversion process.

公开(公告)号：US11340366B2

公开(公告)日：2022-05-24

申请号：US16571950

申请日：2019-09-16

Applicant: Anatoly I. Baumstein ,  Gboyega Ayeni ,  Carey M. Marcinkovich ,  Jaewoo Park ,  Sirui Tan ,  Peter Traynin

Inventor： Anatoly I. Baumstein ,  Gboyega Ayeni ,  Carey M. Marcinkovich ,  Jaewoo Park ,  Sirui Tan ,  Peter Traynin

IPC: G01V1/28 ,  G01V1/30

Abstract: Iterative methods for inversion of seismic data to update a physical property model are disclosed. Such methods may comprise iteratively updating the model until a first predetermined resolution is achieved, using full wavefield inversion of the seismic data up to a first frequency threshold and assuming the seismic data is free of attenuation effects; extracting geobodies from the updated model; obtaining a Q model using the geobodies; and updating the physical property model using an inversion process, wherein the Q model is incorporated into the inversion process. These steps may be repeated until a second predetermined resolution of the physical property model is achieved, wherein the first frequency threshold is progressively increased in each repetition. The Q model may be updated with seismic data at all available frequencies to obtain a full-band Q model; and the physical property model may be updated using full-band migration and the full-band Q model.

公开(公告)号：US20160047924A1

公开(公告)日：2016-02-18

申请号：US14800147

申请日：2015-07-15

Applicant: Christine Krohn ,  Partha S. Routh ,  John E. Anderson ,  Gboyega Ayeni ,  Kenneth E. Green

Inventor： Christine Krohn ,  Partha S. Routh ,  John E. Anderson ,  Gboyega Ayeni ,  Kenneth E. Green

IPC: G01V1/28 ,  G06F17/10

CPC classification number: G01V1/282 ,  G01V1/303 ,  G01V2210/614 ,  G01V2210/67

Abstract: A method, including: obtaining an initial model of a subsurface property; simulating synthetic data from the initial model; obtaining recorded borehole seismic data, wherein the recorded borehole seismic data was obtained with a seismic source or receiver located in a well; and inverting, with a computer, the recorded borehole seismic data by full wavefield inversion, wherein the full wavefield inversion includes comparing the synthetic data to the recorded borehole seismic data, and computing a cost function, obtaining a gradient function from the cost function, wherein the gradient function is related to a change in the objective function with an incremental change in model parameters, using the initial model to compute an illumination function or a resolution function for seismic sources and receivers, and obtaining a conditioned gradient function by conditioning the gradient function with the illumination function or the resolution function.

Abstract translation: 一种方法，包括：获得地下属性的初始模型; 从初始模型模拟合成数据; 获得记录的井眼地震数据，其中记录的井眼地震数据是用位于井中的地震源或接收器获得的; 并且利用计算机反转通过全波场反演的记录的井眼地震数据，其中全波场反演包括将合成数据与记录的井眼地震数据进行比较，并计算成本函数，从成本函数获得梯度函数，其中 梯度函数与目标函数的变化与模型参数的增量变化相关，使用初始模型计算地震源和接收机的照明功能或分辨率函数，并通过调整梯度函数获得调节梯度函数 具有照明功能或分辨率功能。

公开(公告)号：US09977142B2

公开(公告)日：2018-05-22

申请号：US14670039

申请日：2015-03-26

Applicant: Yaxun Tang ,  Gboyega Ayeni

Inventor： Yaxun Tang ,  Gboyega Ayeni

IPC: G01V1/30 ,  G01V1/32 ,  G01V1/28

CPC classification number: G01V1/32 ,  G01V1/282 ,  G01V1/303 ,  G01V1/306 ,  G01V2210/42 ,  G01V2210/63

Abstract: Method for simultaneously inverting full-wavefield seismic data (51) for multiple classes of physical property parameters (e.g., velocity and anisotropy) by computing the gradient (53), i.e. search direction (54), of an objective function for each class of parameters, then applying (preferably exhaustive) first-pass independent line searches to each parameter class to obtain the corresponding step size (55) along the search direction for each parameter class; then without yet updating the model, using the step sizes to define a relative scaling between gradients of all parameter classes. Next, each scaled search direction is recombined to form a new search direction (56), and a new second-pass line search is performed along the new search direction (57), and all parameters are simultaneously updated with the resulting step size (58). Alternatively to the preceding alternating two-pass embodiment, the model may be updated after each first-pass line search, and no second-pass line search is performed.

公开(公告)号：US20180017690A1

公开(公告)日：2018-01-18

申请号：US15614960

申请日：2017-06-06

Applicant: Sirui Tan ,  Yaxun Tang ,  Anatoly I. Baumstein ,  Gboyega Ayeni ,  Tetyana Vdovina ,  Thomas A. Dickens

Inventor： Sirui Tan ,  Yaxun Tang ,  Anatoly I. Baumstein ,  Gboyega Ayeni ,  Tetyana Vdovina ,  Thomas A. Dickens

IPC: G01V1/28 ,  G01V1/30 ,  G01V1/36 ,  G06F17/11

CPC classification number: G01V1/282 ,  G01V1/303 ,  G01V1/368 ,  G06F17/11

Abstract: A method for iteratively inverting seismic data to jointly infer a model for at least P-wave velocity and attenuation parameters of the subsurface, the method including: jointly inverting the P-wave velocity and attenuation parameters with an iterative visco-acoustic full wavefield inversion process, wherein the iterative visco-acoustic full wavefield inversion process includes computing a gradient of an objective function, the objective function measuring a misfit between all or part of the seismic data and corresponding model-simulated seismic data; for each of the P-wave velocity and attenuation parameters, computing a search direction in model space from the gradient; determining line search step sizes α and β for the search directions for the P-wave velocity and attenuation parameters, respectively, wherein a ratio of the step sizes is a function of the P-wave velocity parameter; and using the step sizes α and β and the search directions for each of the P-wave velocity and attenuation parameters, computing a new search direction in model space, then performing a line search along the new search direction to arrive at a new step size, and using the new step size and the new search direction to generate an updated model for a current iteration of the iterative visco-acoustic full wavefield inversion process.

公开(公告)号：US20170010373A1

公开(公告)日：2017-01-12

申请号：US14738471

申请日：2015-06-12

Applicant: Gboyega Ayeni ,  Partha S. Routh ,  John E. Anderson

Inventor： Gboyega Ayeni ,  Partha S. Routh ,  John E. Anderson

IPC: G01V1/28

CPC classification number: G01V1/282 ,  G01V2210/624 ,  G01V2210/626 ,  G01V2210/67

Abstract: A method, including: obtaining initial estimates of a plurality of subsurface parameters; obtaining a recorded wavefield decomposed into a plurality of discrete components; performing, with a computer, a cascaded inversion where the initial estimates of the subsurface parameters are individually updated, wherein each of the subsurface parameters are updated using a different discrete component of the recorded wavefield of the plurality of discrete components; and generating, with the computer, updated subsurface models from the cascaded inversion for each of the subsurface parameters.

Abstract translation: 一种方法，包括：获得多个地下参数的初始估计; 获得分解成多个分立分量的记录波场; 使用计算机执行级联反演，其中地下参数的初始估计被单独更新，其中使用所述多个离散分量中记录的波场的不同离散分量更新每个所述地下参数; 并且利用计算机从每个地下参数的级联反演生成更新的地下模型。

公开(公告)号：US20150323689A1

公开(公告)日：2015-11-12

申请号：US14670039

申请日：2015-03-26

Applicant: Yaxun Tang ,  Gboyega Ayeni

Inventor： Yaxun Tang ,  Gboyega Ayeni

IPC: G01V1/32 ,  G01V1/30 ,  G01V1/28

CPC classification number: G01V1/32 ,  G01V1/282 ,  G01V1/303 ,  G01V1/306 ,  G01V2210/42 ,  G01V2210/63

Abstract: Method for simultaneously inverting full-wavefield seismic data (51) for multiple classes of physical property parameters (e.g., velocity and anisotropy) by computing the gradient (53), i.e. search direction (54), of an objective function for each class of parameters, then applying (preferably exhaustive) first-pass independent line searches to each parameter class to obtain the corresponding step size (55) along the search direction for each parameter class; then without yet updating the model, using the step sizes to define a relative scaling between gradients of all parameter classes. Next, each scaled search direction is recombined to form a new search direction (56), and a new second-pass line search is performed along the new search direction (57), and all parameters are simultaneously updated with the resulting step size (58). Alternatively to the preceding alternating two-pass embodiment, the model may be updated after each first-pass line search, and no second-pass line search is performed.

Abstract translation: 通过计算每个参数类别的目标函数的梯度（53），即搜索方向（54）来同时反转全波场地震数据（51），用于多类物理参数（如速度和各向异性） ，然后对每个参数类进行（优选穷尽的）首过独立行搜索，以获得每个参数类沿着搜索方向的相应步长（55）; 然后没有更新模型，使用步长来定义所有参数类的梯度之间的相对缩放。 接下来，每个缩放的搜索方向被重组以形成新的搜索方向（56），并且沿着新的搜索方向（57）执行新的二次线搜索，并且所有参数同时被更新，所得到的步长（58 ）。 作为上述交替双路实施例的替代方案，可以在每次首次通过行搜索之后更新模型，并且不执行二次线搜索。