Embodiments of dynamic gain adjustments in seismic surveys are described. One method of acquiring a seismic survey includes determining an arrival time at a seismic receiver of a downgoing seismic wavefield associated with a seismic source based at least in part on an estimated position of the seismic source, an estimated position of the seismic receiver, or combinations thereof. The method also includes adjusting a gain of the seismic receiver based at least in part on the determined arrival time of the downgoing seismic wavefield in order to, for example, help prevent overdriving or clipping of the seismic receiver when the downgoing seismic wavefield arrives at or passes by the seismic receiver.

The orientation of the symmetry axis of an underground formation including an HTI layer is determined by comparing azimuthal Fourier coefficient of inversion results in distinct source-receiver azimuth ranges with values expected from the HTI assumption. A branch-stacking technique or prior knowledge may be used to select one of the anisotropy axis orientation values.

One embodiment includes receiving distributed acoustic sensing (DAS) data for responses associated with seismic excitations in an area of interest. The area of interest includes a sea surface, the water column, a seafloor, and a subseafloor. The seismic excitations are generated by at least one seismic source in the area of interest. The responses are detected by at least one fiber optic sensing apparatus configured for DAS that is in the water column, on the seafloor, in a wellbore drilled through the seafloor and into the subseafloor, or any combination thereof. The embodiment includes determining a function of speed of sound in water using the DAS data, and correcting a digital seismic image associated with the area of interest using the function of speed of sound in water to generate a corrected digital seismic image.

The present invention relates to a method of processing seismic data. The method may include calculating a number of calculated structure tensors for each of a number of seismic data lines, the seismic data lines being spatially distributed about an area of the surface of the Earth. The method also may include interpolating the calculated structure tensors to find interpolated structure tensors in a region of the area between the lines of the seismic data lines, and calculating calculated seismic data from the interpolated structure tensors.

Systems, media, and methods for processing seismic data are disclosed. For example, in one embodiment, the method may include receiving a plurality of partial image partitions of a migrated seismic image, and stacking the plurality of partial image partitions such that a first image is generated. The method may also include aligning the plurality of partial image partitions based at least partially on the first image. Aligning may include adjusting at least one of the plurality of partial image partitions and generating a displacement field. The method may also include, based at least in part on the displacement field, stacking the plurality of aligned partial image partitions to generate a second image. The method may further include based at least in part on the second image, realigning the plurality of aligned partial image partitions.

In some embodiments, a time migration diffraction imaging method includes computing a pseudo-depth characterizing a subsurface seismic event by scaling a vertical traveltime using a scaling velocity. A specularity value for a subsurface seismic event is determined according to the pseudo-depth, and a contribution weight for a corresponding seismic trace amplitude is determined according to the specularity value. The specularity value may be determined according to an angle between a traveltime gradient and a normal to a local reflector surface. A diffraction image is generated according to a weighted sum of seismic trace amplitudes. The weighted sum attenuates the contribution of specular events relative to diffraction events.

A method for generating seismic attribute gathers, the method including: computing, with a computer, seismic images with a field dataset; generating, with a computer, synthetic data corresponding to the seismic images; computing, with a computer, an attribute volume by applying an expectation method to the synthetic data; mapping, with a computer, the attribute volume to the seismic images; and generating, with a computer, seismic attribute gathers by stacking the seismic images mapped to the attribute volume.

Computing device, computer instructions and method for estimating a broadband wavelet associated with a given seismic data set. The method includes receiving broadband seismic data; constructing and populating a misfit function; calculating the broadband wavelet based on the misfit function and the broadband seismic data; and estimating physical reservoir properties of a surveyed subsurface based on the broadband wavelet. The broadband wavelet is constrained, through the misfit function, by (1) an amplitude only long wavelet, and (2) an amplitude and phase short wavelet. The amplitude and phase short wavelet is shorter in time than the amplitude only long wavelet.

A method is described that includes obtaining a multicomponent seismic data set for a subterranean region of interest and determining, using a computer processor, a PP stacked time-domain seismic image and a PS stacked time-domain seismic image from the multicomponent seismic data set. The method further includes transforming a recording-time axis of at least one of the PP stacked time-domain seismic image and the PS stacked time-domain seismic image to produce a pair of coarsely-registered PP and PS seismic images and filtering at least one of the pair to produce a pair of spectrally-matched PP and PS seismic images. Further, the method includes dynamically warping at least one of the pair of spectrally-matched PP and PS seismic images to produce a pair of fully-registered PP and PS seismic images.

Method for acquiring, at reduced acquisition cost, seismic data using simultaneous, field-encoded sources in the field (702), and then constructing pseudo source-records (703) that better meet the requirements for using additional simultaneous computer-encoded sourcing for computer simulations or forward modeling (706) as part of (707) iterative FWI (Full Wavefield Inversion) or RTM (Reverse Time Migration), with additional reduction in computational costs. By better meeting the requirements of simultaneous sourcing for FWI or RTM (701), artifacts and crosstalk are reduced in the output. The method can be used for marine streamer acquisition and other non-fixed spread geometries to acquire both positive and negative offsets and to mitigate the “missing data” problem for simultaneous-source FWI. It can also be used for land data to overcome issues with moving spreads and long continuous records.