FWI with areal and point sources
    1.
    发明授权

    公开(公告)号:US10338244B2

    公开(公告)日:2019-07-02

    申请号:US15477417

    申请日:2017-04-03

    IPC分类号: G06F11/30 G01V1/28 G01V1/30

    摘要: A method, including performing, with a computer, up/down separation of geophysical data, which produces an approximate up-going wavefield and an approximate down-going wavefield; creating an areal source based at least in part on the down-going wavefield; and performing, with a computer, a full wavefield inversion process with the areal source, and an objective function measuring a misfit between modeled up-going wavefields and recorded up-going wavefields, wherein the full wavefield inversion process generates a final subsurface physical property model.

    Joint full wavefield inversion of P-wave velocity and attenuation using an efficient first order optimization

    公开(公告)号:US10459096B2

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

    申请号:US15614960

    申请日:2017-06-06

    摘要: 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.

    Joint Full Wavefield Inversion of P-Wave Velocity and Attenuation Using an Efficient First Order Optimization

    公开(公告)号:US20180017690A1

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

    申请号:US15614960

    申请日:2017-06-06

    摘要: 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.

    Full-Wavefield Inversion Using Mirror Source-Receiver Geometry
    6.
    发明申请
    Full-Wavefield Inversion Using Mirror Source-Receiver Geometry 审中-公开
    使用镜像源接收几何的全波场反演

    公开(公告)号:US20160061974A1

    公开(公告)日:2016-03-03

    申请号:US14801409

    申请日:2015-07-16

    IPC分类号: G01V1/28 G06F17/10

    摘要: Method for performing a full wavefield inversion (FWI) without simulating free-surface multiple reflections. The free-surface multiples are removed from the field gathers of seismic data, which are then used to generate a subsurface velocity model by FWI. In the FWI, the field monopole sources and receivers are replaced with dipole (actual and mirror image) sources and receivers (21) when model-simulating (23) synthetic survey data. Also, direct arrivals at the mirror receiver locations are preferably simulated (25) with the dipole sources for each shot location and added (26) to the synthetic survey data (24) for that shot location, resulting in corrected synthetic survey data (27), which is used in the FWI to generate residuals. A model update may be computed by back-propagating the residuals by injecting them at both mirror and actual receiver locations.

    摘要翻译: 在不模拟自由表面多重反射的情况下执行全波场反演(FWI)的方法。 自由表面倍数从地震数据的场收集中去除,然后用FWI生成地下速度模型。 在FWI中,当模拟(23)综合测量数据时,场单极子源和接收器被偶极子(实际和镜像)源和接收器(21)所取代。 此外,优选地,在镜接收器位置处的直达到达用于每个镜头位置的偶极子源(25),并且将(26)添加到用于该镜头位置的综合测量数据(24),导致校正的综合测量数据(27) ,在FWI中用于产生残差。 模型更新可以通过在两个镜像和实际的接收器位置注入来反向传播残差来计算。

    Method for Temporal Dispersion Correction for Seismic Simulation, RTM and FWI
    7.
    发明申请
    Method for Temporal Dispersion Correction for Seismic Simulation, RTM and FWI 审中-公开
    用于地震模拟,RTM和FWI的时间色散校正方法

    公开(公告)号:US20150355356A1

    公开(公告)日:2015-12-10

    申请号:US14670085

    申请日:2015-03-26

    IPC分类号: G01V1/32 G01V1/48 G01V1/36

    摘要: Method for correcting seismic simulations, RTM, and FWI for temporal dispersion due to temporal finite difference methods in which time derivatives are approximated to a specified order of approximation. Computer-simulated seismic data (51) are transformed from time domain to frequency domain (52), and then resampled using a mapping relationship that maps, in the frequency domain, to a frequency at which the time derivative exhibits no temporal dispersion (53), or to a frequency at which the time derivative exhibits a specified different order of temporal dispersion. Alternatively, measured seismic data from a field survey (61) may have temporal dispersion of a given order introduced, by a similar technique, to match the order of approximation used to generate simulated data which are to be compared to the measured data.

    摘要翻译: 用于由时间有限差分法校正地震模拟,RTM和FWI用于时间色散的方法,其中时间导数近似为指定的近似次数。 计算机模拟地震数据(51)从时域(52)变换到频域(52),然后使用在频域中映射到时间导数不显示时间色散的频率的映射关系重新采样(53) 或时间导数呈现指定的不同时间色散顺序的频率。 或者,来自现场勘测(61)的测量的地震数据可以通过类似的技术具有引入的给定次序的时间色散,以与用于生成将与测量数据进行比较的模拟数据的近似的顺序相匹配。

    Poynting Vector Minimal Reflection Boundary Conditions
    10.
    发明申请
    Poynting Vector Minimal Reflection Boundary Conditions 审中-公开
    坡印矢量最小反射边界条件

    公开(公告)号:US20160223697A1

    公开(公告)日:2016-08-04

    申请号:US14887621

    申请日:2015-10-20

    IPC分类号: G01V1/28 G06F17/16

    摘要: A method for exploring for hydrocarbons, including: simulating a seismic waveform, using a computer, wherein computations are performed on a computational grid representing a subsurface region, said computational grid using perfectly matched layer (PML) boundary conditions that use an energy dissipation operator to minimize non-physical wave reflections at grid boundaries; wherein, in the simulation, the PML boundary conditions are defined to reduce computational instabilities at a boundary by steps including, representing direction of energy propagation by a Poynting vector, and dissipating energy, with the dissipation operator, in a direction of energy propagation instead of in a phase velocity direction; and using the simulated waveform in performing full waveform inversion or reverse time migration of seismic data, and using a physical property model from the inversion or a subsurface image from the migration to explore for hydrocarbons.

    摘要翻译: 一种用于探索碳氢化合物的方法,包括:使用计算机模拟地震波形,其中在表示地下区域的计算网格上执行计算,所述计算网格使用完全匹配层(PML)边界条件,其使用能量耗散算子 最小化网格边界处的非物理波反射; 其中,在模拟中,定义了PML边界条件以减少边界处的计算不稳定性,所述步骤包括表示坡坡廷矢量的能量传播方向,以及耗散算子在能量传播方向上耗散能量,而不是 在相速度方向上; 并使用模拟波形进行地震数据的全波形反演或反向时间偏移,并使用来自反演的物理性质模型或来自迁移的地下图像来探索碳氢化合物。