Abstract: A multi-stage FWI workflow uses multiple-contaminated FWI models to predict surface-related multiples. A method embodying the present technological advancement, can include: using data with free surface multiples as input into FWI; generating a subsurface model by performing FWI with the free-surface boundary condition imposed on top of the subsurface model; using inverted model from FWI to predict multiples; removing predicted multiples from the measured data; using the multiple-free data as input into FWI with absorbing boundary conditions imposed on top of the subsurface model; and preparing a multiple free data set for use in conventional seismic data processing.

Abstract: A method, including: obtaining a seismic dataset that is separated into subsets according to predetermined subsurface reflection angle ranges; performing, with a computer, an acoustic full wavefield inversion process on each of the subsets, respectively, to invert for density and generate respective density models; generating acoustic impedances for each of the subsets, as a function of reflection angle, using the respective density models; and transforming, using a computer, the acoustic impedances for each of the subsets into reflectivity sections, wherein the transforming includes normalizing the reflectivity sections by their respective bandwidth.

Abstract: A deterministic method for selecting a set of encoding weights for simultaneous encoded-source inversion of seismic data that will cause the iterative inversion to converge faster than randomly chosen weights. The encoded individual source gathers are summed (83), forming a composite gather, and simulated in a single simulation operation. The invention creates multiple realizations of the simulation (84), each with its own encoding vector (82) whose components are the weights for the shots in the composite gather. The encoding vectors of the invention are required to be orthogonal (82), which condition cannot be satisfied by random weights, and in various embodiments of the invention are related to eigenvectors of a Laplacian matrix, sine or cosine functions, or Chebyshev nodes as given by the roots of Chebyshev polynomials. For non-fixed receiver geometry, an encoded mask (61) may be used to approximately account for non-listening receivers.

Abstract: A method, including: obtaining a velocity model generated by an acoustic full wavefield inversion process; generating, with a computer, a variable Q model by applying pseudo-Q migration on processed seismic data of a subsurface region, wherein the velocity model is used as a guided constraint in the pseudo-Q migration; and generating, with a computer, a final subsurface velocity model that recovers amplitude attenuation caused by gas anomalies in the subsurface region by performing a visco-acoustic full wavefield inversion process, wherein the variable Q model is fixed in the visco-acoustic full wavefield inversion process.

Abstract: The present disclosure provides a system and method for estimating fracture density within a subsurface formation from S-wave seismic data. In one embodiment, the S-wave seismic data is separated into fast (“S1”) and slow (“S2”) data. A computer is used to compute local similarity of the S1 and S2 data and to compute a cumulative time-difference by which the S2 data lags the S1 data from the local similarity. Based on the computed cumulative time-difference, the fracture density of a subsurface formation is estimated.

Abstract: Method for using the full wavefield (primaries, internal multiples and free-surface multiples) in inversion of marine seismic data, including both pressure and vertical velocity data (21), to infer a subsurface model of acoustic velocity or other physical property. The marine seismic data are separated (22) into up-going (23) and down-going (24) wavefields, and both wavefields are inverted in a joint manner, in which the final model is impacted by both wavefields. This may be achieved by inverting both wavefields simultaneously (25), or one after the other, i.e. in a cascaded approach (35?37, or 45?47), for the subsurface properties (26, 38, 48).

Abstract: A method, including: obtaining a seismic dataset that is separated into subsets according to predetermined subsurface reflection angle ranges; performing, with a computer, an acoustic full wavefield inversion process on each of the subsets, respectively, to invert for density and generate respective density models; generating acoustic impedances for each of the subsets, as a function of reflection angle, using the respective density models; and transforming, using a computer, the acoustic impedances for each of the subsets into reflectivity sections, wherein the transforming includes normalizing the reflectivity sections by their respective bandwidth.

Abstract: A method, including: obtaining a velocity model generated by an acoustic full wavefield inversion process; generating, with a computer, a variable Q model by applying pseudo-Q migration on processed seismic data of a subsurface region, wherein the velocity model is used as a guided constraint in the pseudo-Q migration; and generating, with a computer, a final subsurface velocity model that recovers amplitude attenuation caused by gas anomalies in the subsurface region by performing a visco-acoustic full wavefield inversion process, wherein the variable Q model is fixed in the visco-acoustic full wavefield inversion process.

Abstract: A multi-stage FWI workflow uses multiple-contaminated FWI models to predict surface-related multiples. A method embodying the present technological advancement, can include: using data with free surface multiples as input into FWI; generating a subsurface model by performing FWI with the free-surface boundary condition imposed on top of the subsurface model; using inverted model from FWI to predict multiples; removing predicted multiples from the measured data; using the multiple-free data as input into FWI with absorbing boundary conditions imposed on top of the subsurface model; and preparing a multiple free data set for use in conventional seismic data processing.

Abstract: The present disclosure provides a system and method for inferring one or more physical property parameters of a subsurface media by inverting converted wave data acquired during a seismic survey. Composite seismic traces are generated at a plurality of survey azimuths (step 609). These composite traces are composed such that their amplitudes are free of effects of subsurface anisotropy. At least one of the generated composite seismic traces is then inverted by isotropic inversion to determine a property parameter of the subsurface media (step 623).

Abstract: 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.

Abstract: A deterministic method for selecting a set of encoding weights for simultaneous encoded-source inversion of seismic data that will cause the iterative inversion to converge faster than randomly chosen weights. The encoded individual source gathers are summed (83), forming a composite gather, and simulated in a single simulation operation. The invention creates multiple realizations of the simulation (84), each with its own encoding vector (82) whose components are the weights for the shots in the composite gather. The encoding vectors of the invention are required to be orthogonal (82), which condition cannot be satisfied by random weights, and in various embodiments of the invention are related to eigenvectors of a Laplacian matrix, sine or cosine functions, or Chebyshev nodes as given by the roots of Chebyshev polynomials. For non-fixed receiver geometry, an encoded mask (61) may be used to approximately account for non-listening receivers.

Abstract: Method for using the full wavefield (primaries, internal multiples and free-surface multiples) in inversion of marine seismic data, including both pressure and vertical velocity data (21), to infer a subsurface model of acoustic velocity or other physical property. The marine seismic data are separated (22) into up-going (23) and down-going (24) wavefields, and both wavefields are inverted in a joint manner, in which the final model is impacted by both wavefields. This may be achieved by inverting both wavefields simultaneously (25), or one after the other, i.e. in a cascaded approach (35?37, or 45?47), for the subsurface properties (26, 38, 48).

Abstract: The present disclosure provides a system and method for inferring one or more physical property parameters of a sub-surface media by inverting converted wave data acquired during a seismic survey. Composite seismic traces are generated at a plurality of survey azimuths (step 609). These composite traces are composed such that their amplitudes are free of effects of subsurface anisotropy. At least one of the generated composite seismic traces is then inverted by isotropic inversion to determine a property parameter of the subsurface media (step 623).

Abstract: The present disclosure provides a system and method for estimating fracture density within a subsurface formation from S-wave seismic data. In one embodiment, the S-wave seismic data is separated into fast (“S1”) and slow (“S2”) data. A computer is used to compute local similarity of the S1 and S2 data and to compute a cumulative time-difference by which the S2 data lags the S1 data from the local similarity. Based on the computed cumulative time-difference, the fracture density of a subsurface formation is estimated.

Abstract: Method for determining fracture orientation and fracture intensity in multiple fractured layers in the subsurface in a layer-stripping manner. Multi-component, multi-azimuth seismic data are required (31), from which the horizontal, primarily converted wave, components are selected, and these data are further reduced by selecting only the data for which the survey azimuths are either parallel or perpendicular to the general fracture strike (33). If the general fracture trend is unknown, such selective data may be determined by an azimuth-offset scanning process. Layer stripping is performed on azimuth/offset stacks (42) to produce fracture parameter maps (43). All offsets are stacked in those azimuths that produce consistent fracture parameter maps (44), then layer stripping is performed (45) on the stacks to produce final fracture orientation and S-wave time difference maps (46). These maps can be used to produce true amplitude fast and slow S-waves (56).

Abstract: Method for aligning converted wave seismic reflection data (PS data) with conventional PP seismic reflection data so that both data types may be used to more accurately image the subsurface for hydrocarbon exploration or field development. Amplitude vs. angle (AVA) or amplitude vs. offset (AVO) attributes of PP and PS seismic data are identified and defined, which attributes are theoretically expected to be in phase and optimize seismic resolution in the data. In one embodiment of the invention, such attributes are calculated (310), then the same horizons are identified in a series of PP attributes and in a series of PS attributes, then the second series is aligned with the first at the horizon locations (316, 320), then a time transfer function is generated and applied to the PS mode data (322), and the aligned joint-mode data are inverted (326) using, for example, AVA attributes.

Abstract: Method for aligning converted wave seismic reflection data (PS data) with conventional PP seismic reflection data so that both data types may be used to more accurately image the subsurface for hydrocarbon exploration or field development. Amplitude vs. angle (AVA) or amplitude vs. offset (AVO) attributes of PP and PS seismic data are identified and defined, which attributes are theoretically expected to be in phase and optimize seismic resolution in the data. In one embodiment of the invention, such attributes are calculated (310), then the same horizons are identified in a series of PP attributes and in a series of PS attributes, then the second series is aligned with the first at the horizon locations (316, 320), then a time transfer function is generated and applied to the PS mode data (322), and the aligned joint-mode data are inverted (326) using, for example, AVA attributes.