Abstract: A method and apparatus for imaging seismic data includes obtaining an initial model of a subsurface formation, wherein the model includes a plurality of nodes that form at least part of a grid; an initial dip value for the nodes; and a set of origin coordinates for each of the nodes; performing bottom-up ray tracing for each node in the model, resulting in a set of arrival coordinates for each node; identifying a plurality of gathers from the seismic data; for each gather: calculating a set of midpoint coordinates; defining a midpoint vicinity surrounding the set of midpoint coordinates; identifying the nodes having arrival coordinates within the midpoint vicinity; and estimating a unique aperture for each of the gathers based on the respective origin coordinates; storing the estimated apertures in a table; and generating a subsurface volume or image with subsurface reflectors determined with apertures of the respective gathers.

Abstract: Method for locating fault lines or surfaces in 2-D or 3-D seismic data based on the fact that fault discontinuities in the space domain span a wide range in a local slowness (slope) domain, whereas other dipping events in the space domain data, such as noise, tend to be coherent, and hence will appear focused in the slowness dimension. Therefore, the method comprises decomposing the seismic data (102) by a transformation to the local slowness domain, preferably using Gaussian slowness period packets as the local slowness or slope decomposition technique, thereby avoiding problems with the data stationary assumption. In the local slowness domain, faults may be identified (104) using the principle mentioned above, i.e. that faults are represented as a truncation in the space domain data, hence they will appear broadband in the slowness dimension.

Abstract: A method and apparatus for imaging seismic data includes obtaining an initial model of a subsurface formation, wherein the model includes a plurality of nodes that form at least part of a grid; an initial dip value for the nodes; and a set of origin coordinates for each of the nodes; performing bottom-up ray tracing for each node in the model, resulting in a set of arrival coordinates for each node; identifying a plurality of gathers from the seismic data; for each gather: calculating a set of midpoint coordinates; defining a midpoint vicinity surrounding the set of midpoint coordinates; identifying the nodes having arrival coordinates within the midpoint vicinity; and estimating a unique aperture for each of the gathers based on the respective origin coordinates; storing the estimated apertures in a table; and generating a subsurface volume or image with subsurface reflectors determined with apertures of the respective gathers.

Abstract: Automated detection of microseismic events is performed by a system and method that provides automated calculation of one or more statistical parameters for statistics-based detection of microseismic events of interest. An event detector to perform signal subspace-based detection may automatically, e.g. iteratively, be configured based on the one or more automatically calculated and/or estimated statistical parameters, which may include a signal subspace dimension, an effective embedding dimension, a detection threshold, and/or a false alarm rate. Some of the statistical parameters may be calculated based on custom-generated synthetic seismic data.

Abstract: Method for locating fault lines or surfaces in 2-D or 3-D seismic data based on the fact that fault discontinuities in the space domain span a wide range in a local slowness (slope) domain, whereas other dipping events in the space domain data, such as noise, tend to be coherent, and hence will appear focused in the slowness dimension. Therefore, the method comprises decomposing the seismic data (102) by a transformation to the local slowness domain, preferably using Gaussian slowness period packets as the local slowness or slope decomposition technique, thereby avoiding problems with the data stationary assumption. In the local slowness domain, faults may be identified (104) using the principle mentioned above, i.e. that faults are represented as a truncation in the space domain data, hence they will appear broadband in the slowness dimension.

Abstract: Automated detection of microseismic events is performed by a system and method that provides automated calculation of one or more statistical parameters for statistics-based detection of microseismic events of interest. An event detector to perform signal subspace-based detection may automatically, e.g. iteratively, be configured based on the one or more automatically calculated and/or estimated statistical parameters, which may include a signal subspace dimension, an effective embedding dimension, a detection threshold, and/or a false alarm rate. Some of the statistical parameters may be calculated based on custom-generated synthetic seismic data.