Abstract: A method is described for predicting and preventing wellbore interactions at wells that are near the injection well. The method includes receiving fiber optics data; performing object detection by detecting object-like events in the fiber optic data; and sending instructions to a hydraulic fracturing system based on the object detection. The method is executed by a computer system.

Abstract: A method is described for monitoring a stimulated reservoir volume (SRV) including receiving simulation parameters, performing 3D fully coupled quasi-static poro-elastic finite difference modeling using the simulation parameters, wherein the 3D fully coupled quasi-static poro-elastic finite difference modeling is based on a rescaling of solid rock and fluid flow density parameters and generates simulated temporal quasi-static stresses, and pore pressure. In addition, simulated stresses may be used for performing calculation of the 3D rotation of the simulated stresses to principal directions; performing calculation of the temporal 3D Mohr-Coulomb (MC) failure criteria from the calculated principal stresses and the simulated pore pressure for all or selected time steps; and displaying the computed temporal MC failure criteria results on a graphical display. The method may also be used in time-lapse monitoring of the reservoir for microseismic depletion delineation.

Abstract: A method is described for generating a subsurface model using stochastic full waveform inversion by receiving a seismic dataset representative of a subsurface volume of interest; performing stochastic full waveform inversion of the seismic dataset to generate a long wavelength subsurface model; and performing full waveform inversion of the seismic dataset using the long wavelength subsurface model as a starting model to generate an improved subsurface model. The method may further include performing seismic imaging of the seismic dataset using the improved subsurface model to generate a seismic image and identifying geologic features based on the seismic image. The method may be executed by a computer system.

Abstract: A method is described for predicting and preventing wellbore interactions at wells that are near the injection well. The method includes receiving fiber optics data; performing object detection by detecting object-like events in the fiber optic data; and sending instructions to a hydraulic fracturing system based on the object detection. The method is executed by a computer system.

Abstract: A method is described for property estimation including receiving a seismic dataset representative of a subsurface volume of interest and a well log from a well location within the subsurface volume of interest; identifying seismic traces in the seismic dataset that correspond to the well location to obtain a subset of seismic traces; windowing the subset of seismic traces and the well log to generate windowed seismic traces and a windowed well log; multiplying the windowed seismic traces and the windowed well log by a random matrix to generate a plurality of training datasets; and training a neural network using the plurality of training datasets. The method may be executed by a computer system.

Abstract: A method is described for property estimation including receiving a seismic dataset representative of a subsurface volume of interest and a well log from a well location within the subsurface volume of interest; identifying seismic traces in the seismic dataset that correspond to the well location to obtain a subset of seismic traces; windowing the subset of seismic traces and the well log to generate windowed seismic traces and a windowed well log; multiplying the windowed seismic traces and the windowed well log by a random matrix to generate a plurality of training datasets; and training a neural network using the plurality of training datasets. The method may be executed by a computer system.

Abstract: Geophysical data processing is remotely controlled and monitored over a wide-area network such as the Internet. A customer using a client computer builds geophysical data processing flows (concatenations of geophysical data processing modules or filters) and enters parameter values required for flow execution. The flow descriptions and associated parameter values are then transferred from the client to a geophysical data processing server, for example a parallel supercomputer. The flows (jobs) are executed on the server, typically over periods ranging from hours to weeks. Intermediate or partial results are made available to the customer for visualization before the processing of a flow is complete. The customer can then modify the flow before its complete execution. Data-entry windows are automatically generated for geophysical processing modules by parsing the source code of the modules.

Abstract: A method of performing migration velocity analysis may include: obtaining seismic data and an initial velocity model; determining reflection points; deriving a wavepath backprojection operator based on the initial velocity model and the reflection points by constructing wavepaths from each reflection point of the reflection points; and performing a traveltime inversion using the wavepath backprojection operator. The initial velocity model may be updated based on the traveltime inversion. Determining reflection points may be automated by calculating reflection points based on results from a depth migration algorithm performed on the initial velocity model. Selection of residual moveout values may be automated by selecting based on a dip field for each prestack gather obtained from a depth migration algorithm performed on the initial velocity model. Residual traveltimes may be calculated using the selected residual moveout values. The residual traveltimes may be used in the traveltime inversion.

Abstract: In one embodiment, a computer-implemented seismic viewing/editing collaboration method includes performing real-time collaborative cursor tracking, copaging, picking, and image manipulation in a distributed-display-processing, peer-to-peer architecture. A parameterized, minimal set of information required to update a display is transferred directly between different clients. A group state containing events generated by different clients is enforced to be synchronized on the different clients.

Abstract: In one embodiment, a computer-implemented common azimuth migration seismic data processing method comprises: providing a common-azimuth input data set for a geophysical data processing volume of interest; providing a velocity model for the volume; applying an offset antialiasing operator to the input data set; and performing a recursive downward-continuation of the common-azimuth input data set to a plurality of successive common-azimuth surfaces to generate an image of the volume of interest. In one embodiment, the present invention further provides for selecting a depth dependence of an offset range employed in the downward continuation; selecting a frequency-dependence of a depth step size employed in the downward continuation; selecting a frequency dependence of a cutoff depth employed in the downward continuation; and adding reciprocal traces to the data around zero offset, for reducing imaging artifacts introduced by data edge effects.

Abstract: A seismic velocity analysis method includes tying velocity parameter values such as residual velocity values to geological horizons (reflectors) within a seismic exploration volume. Common image gathers (CIGs) such a common reflection point (CRP) gathers or angle-domain common image gathers (ACIGs) are generated for a set of CIG grid points. Computed best-fit residual velocity values are then snapped to a neighboring horizon or vertically interpolated to the horizon, to generate residual velocity values along the horizon. The residual velocity values for points along the horizon are then selectively employed in updating the velocity model for the volume of interest.

Abstract: Input seismic data are re-gridded to an arbitrary output grid by output-based azimuth moveout. An input seismic data set corresponding to an input grid is used to generate an equivalent output seismic data set corresponding to an output grid different from the input grid. Preferably, the output grid is divided into blocks, and each output grid block is assigned to one of a plurality of independent parallel processors. For each output trace corresponding to an output location, the contributions of plural input traces to the output trace are computed according to an azimuth moveout operator. The contributions are then summed into the output trace.

Abstract: Geophysical data processing is remotely controlled and monitored over a wide-area network such as the Internet. A customer using a client computer builds geophysical data processing flows (concatenations of geophysical data processing modules or filters) and enters parameter values required for flow execution. The flow descriptions and associated parameter values are then transferred from the client to a geophysical data processing server, for example a parallel supercomputer. The flows (jobs) are executed on the server, typically over periods ranging from hours to weeks. Intermediate or partial results are made available to the customer for visualization before the processing of a flow is complete. The customer can then modify the flow before its complete execution. Data-entry windows are automatically generated for geophysical processing modules by parsing the source code of the modules.

Abstract: Migration velocity analysis is performed using Angle-Domain Common Image Gathers (ACIGs). When the correct velocity model is employed for migration, all ACIG events corresponding to a subsurface location are aligned along a horizontal line. Residual moveout can be performed on each ACIG with a suite of trial residual velocity values, according to an angle-domain residual moveout equation. A best-fit residual velocity value that leads to horizontally-aligned events upon moveout can be selected by generating a distribution of semblance (amplitude summed over a given depth) over residual velocity. Best-fit residual velocity values corresponding to selected subsurface points can be employed to update the initial velocity model using a vertical update, normal ray update, or tomographic update method.

Abstract: Geophysical data processing is remotely controlled and monitored over a wide-area network such as the Internet. A customer using a client computer builds geophysical data processing flows (concatenations of geophysical data processing modules or filters) and enters parameter values required for flow execution. The flow descriptions and associated parameter values are then transferred from the client to a geophysical data processing server, for example a parallel supercomputer. The flows (jobs) are executed on the server, typically over periods ranging from hours to weeks. Intermediate or partial results are made available to the customer for visualization before the processing of a flow is complete. The customer can then modify the flow before its complete execution. Data-entry windows are automatically generated for geophysical processing modules by parsing the source code of the modules.

Abstract: Migration velocity analysis is performed using Angle-Domain Common Image Gathers (ACIGs). When the correct velocity model is employed for migration, all ACIG events corresponding to a subsurface location are aligned along a horizontal line. Residual moveout can be performed on each ACIG with a suite of trial residual velocity values, according to an angle-domain residual moveout equation. A best-fit residual velocity value that leads to horizontally-aligned events upon moveout can be selected by generating a distribution of semblance (amplitude summed over a given depth) over residual velocity. Best-fit residual velocity values corresponding to selected subsurface points can be employed to update the initial velocity model using a vertical update, normal ray update, or tomographic update method.

Abstract: A seismic velocity analysis method includes tying velocity parameter values such as residual velocity values to geological horizons (reflectors) within a seismic exploration volume. Common image gathers (CIGs) such a common reflection point (CRP) gathers or angle-domain common image gathers (ACIGs) are generated for a set of CIG grid points. Computed best-fit residual velocity values are then snapped to a neighboring horizon or vertically interpolated to the horizon, to generate residual velocity values along the horizon. The residual velocity values for points along the horizon are then selectively employed in updating the velocity model for the volume of interest.

Abstract: A computationally efficient semi-recursive Kirchhoff migration algorithm is capable of obtaining accurate images of complex structures in complex geological environments which commonly give rise to multivalued arrivals. The method alternates wave-equation datuming (downward continuation) and Kirchhoff migration. Breaking up the complicated velocity structure into small depth regions allows traveltimes to be calculated in regions where the computation is well behaved and where the computation corresponds to energetic arrivals. Because traveltimes are computed for small depth domains, the effects of multivalued arrivals are implicitly accounted for in the migration and the adverse effects of caustics and headwaves do not develop. The invention can use any simple first-arrival traveltime algorithm, and thus benefits from the computational efficiency, robustness, and simplicity of such algorithms.