Abstract: A method, including: storing, in a computer storage device, geophysical seismic data that has been separated into a multiple-free component and a multiple contaminated component; performing, with a processor, a first full wavefield inversion process on the multiple-free component of the seismic data, wherein a first subsurface physical property model is generated; determining, with a processor, an extended target reflectivity, wherein the extended target reflectivity includes a reflectivity for each of a plurality of shots; separately performing, with a processor, a second full wavefield inversion process with the multiple contaminated component of the seismic data for each of the plurality of shots using the reflectivity corresponding to each of the plurality of shots, wherein a second subsurface physical property model is generated; and generating, with a processor, multiple-free final subsurface physical property model by combining the first subsurface physical property model and the second subsurface phys

Abstract: Method for performing simultaneous encoded-source inversion of geophysical data to estimate parameters of a physical property model (41), especially adapted for surveys without fixed-receiver acquisition geometry, such as marine seismic surveys with moving source and receivers. The encoding functions (32) used on the sources to generate one or more simultaneous encoded-source gathers of data (35), as well as to simulate the same (34), are orthogonal or pseudo-orthogonal with respect to cross-correlation. In addition, receivers are also encoded, with the receiver encoding being designed to make a given receiver less sensitive to sources to which it was not listening during the survey (38). The encoding functions may be temporal bandpass filters differing one from another by central frequency, phase, or both. Efficiency of the method may be further improved by grouping several sources into a super-source, grouping the corresponding gathers into a super-gather, and then applying the above encoding strategy.

Abstract: Method for reducing artifacts in a subsurface physical properties model (120) inferred by iterative inversion (140) of geophysical data (130), wherein the artifacts are associated with some approximation (110) made during the iterative inversion. In the method, some aspect of the approximation is changed (160) as the inversion is iterated such that the artifacts do not increase by coherent addition.

Abstract: The invention includes a method for reducing noise in migration of seismic data, particularly advantageous for imaging by simultaneous encoded source reverse-time migration (SS-RTM). One example embodiment includes the steps of obtaining a plurality of initial subsurface images; decomposing each of the initial subsurface images into components; identifying a set of components comprising one of (i) components having at least one substantially similar characteristic across the plurality of initial subsurface images, and (ii) components having substantially dissimilar characteristics across the plurality of initial subsurface images; and generating an enhanced subsurface image using the identified set of components. For SS-RTM, each of the initial subsurface images is generated by migrating several sources simultaneously using a unique random set of encoding functions. Another embodiment of the invention uses SS-RTM for velocity model building.

Abstract: A method, including: storing, in a computer storage device, geophysical seismic data that has been separated into a multiple-free component and a multiple contaminated component; performing, with a processor, a first full wavefield inversion process on the multiple-free component of the seismic data, wherein a first subsurface physical property model is generated; determining, with a processor, an extended target reflectivity, wherein the extended target reflectivity includes a reflectivity for each of a plurality of shots; separately performing, with a processor, a second full wavefield inversion process with the multiple contaminated component of the seismic data for each of the plurality of shots using the reflectivity corresponding to each of the plurality of shots, wherein a second subsurface physical property model is generated; and generating, with a processor, multiple-free final subsurface physical property model by combining the first subsurface physical property model and the second subsurface phys

Abstract: Method for reducing the time needed to perform geophysical inversion by using simultaneous encoded sources in the simulation steps of the inversion process. The geophysical survey data are prepared by encoding (3) a group of source gathers (1), using for each gather a different encoding signature selected from a set (2) of non-equivalent encoding signatures. Then, the encoded gathers are summed (4) by summing all traces corresponding to the same receiver from each gather, resulting in a simultaneous encoded gather. (Alternatively, the geophysical data are acquired from simultaneously encoded sources.) The simulation steps needed for inversion are then calculated using a particular assumed velocity (or other physical property) model (5) and simultaneously activated encoded sources using the same encoding scheme used on the measured data. The result is an updated physical properties model (6) that may be further updated (7) by additional iterations.

Abstract: Method for reducing artifacts in a subsurface physical properties model (120) inferred by iterative inversion (140) of geophysical data (130), wherein the artifacts are associated with some approximation (110) made during the iterative inversion. In the method, some aspect of the approximation is changed (160) as the inversion is iterated such that the artifacts do not increase by coherent addition.

Abstract: Method for reducing artifacts in a subsurface physical properties model (120) inferred by iterative inversion (140) of geophysical data (130), wherein the artifacts are associated with some approximation (110) made during the iterative inversion. In the method, some aspect of the approximation is changed (160) as the inversion is iterated such that the artifacts do not increase by coherent addition.

Abstract: A method for efficient inversion of measured geophysical data from a subsurface region to prospect for hydrocarbons. Gathers of measured data (40) are encoded (60) using a set of non-equivalent encoding functions (30). Then all data records in each encoded gather that correspond to a single receiver are summed (60), repeating for each receiver to generate a simultaneous encoded gather (80). The method employs iterative, local optimization of a cost function to invert the encoded gathers of simultaneous source data. An adjoint method is used to calculate the gradients of the cost function needed for the local optimization process (100). The inverted data yields a physical properties model (110) of the subsurface region that, after iterative updating, can indicate presence of accumulations of hydrocarbons.

Abstract: Method for simultaneous full-wavefield inversion of gathers of source (or receiver) encoded geophysical data to determine a physical properties model (118) for a subsurface region, especially suitable for surveys where fixed receiver geometry conditions were not satisfied in the data acquisition. Simultaneous source separation (104) is performed to lessen any effect of the measured geophysical data's not satisfying the fixed-receiver assumption. A data processing step (106) coming after the simultaneous source separation acts to conform model-simulated data (105) to the measured geophysical data (108) for source and receiver combinations that are missing in the measured geophysical data.

Abstract: Method for reducing artifacts in a subsurface physical properties model (120) inferred by iterative inversion (140) of geophysical data (130), wherein the artifacts are associated with some approximation (110) made during the iterative inversion. In the method, some aspect of the approximation is changed (160) as the inversion is iterated such that the artifacts do not increase by coherent addition.

Abstract: Method for reducing artifacts in a subsurface physical properties model (120) inferred by iterative inversion (140) of geophysical data (130), wherein the artifacts are associated with some approximation (110) made during the iterative inversion. In the method, some aspect of the approximation is changed (160) as the inversion is iterated such that the artifacts do not increase by coherent addition.

Abstract: Method for simultaneous full-wavefield inversion of gathers of source (or receiver) encoded (30) geophysical data (80) to determine a physical properties model (20) for a subsurface region, especially suitable for surveys where fixed-receiver geometry conditions were not satisfied in the data acquisition (40). The inversion involves optimization of a cross-correlation objective function (100).

Abstract: Method for converting seismic data to obtain a subsurface model of, for example, bulk modulus or density. The gradient of an objective function is computed (103) using the seismic data (101) and a background subsurface medium model (102). The source and receiver illuminations are computed in the background model (104). The seismic resolution volume is computed using the velocities of the background model (105). The gradient is converted into the difference subsurface model parameters (106) using the source and receiver illumination, seismic resolution volume, and the background subsurface model. These same factors may be used to compensate seismic data migrated by reverse time migration, which can then be related to a subsurface bulk modulus model. For iterative inversion, the difference subsurface model parameters (106) are used as preconditioned gradients (107).

Abstract: Method for performing simultaneous encoded-source inversion of geophysical data to estimate parameters of a physical property model (41), especially adapted for surveys without fixed-receiver acquisition geometry, such as marine seismic surveys with moving source and receivers. The encoding functions (32) used on the sources to generate one or more simultaneous encoded-source gathers of data (35), as well as to simulate the same (34), are orthogonal or pseudo-orthogonal with respect to cross-correlation. In addition, receivers are also encoded, with the receiver encoding being designed to make a given receiver less sensitive to sources to which it was not listening during the survey (38). The encoding functions may be temporal bandpass filters differing one from another by central frequency, phase, or both. Efficiency of the method may be further improved by grouping several sources into a super-source, grouping the corresponding gathers into a super-gather, and then applying the above encoding strategy.

Abstract: Method for reducing the time needed to perform geophysical inversion by using simultaneous encoded sources in the simulation steps of the inversion process. The geophysical survey data are prepared by encoding (3) a group of source gathers (1), using for each gather a different encoding signature selected from a set (2) of non-equivalent encoding signatures. Then, the encoded gathers are summed (4) by summing all traces corresponding to the same receiver from each gather, resulting in a simultaneous encoded gather. (Alternatively, the geophysical data are acquired from simultaneously encoded sources.) The simulation steps needed for inversion are then calculated using a particular assumed velocity (or other physical property) model (5) and simultaneously activated encoded sources using the same encoding scheme used on the measured data. The result is an updated physical properties model (6) that may be further updated (7) by additional iterations.

Abstract: Method for simultaneous full-wavefield inversion of gathers of source (or receiver) encoded geophysical data to determine a physical properties model for a subsurface region, especially suitable for surveys where fixed receiver geometry conditions were not satisfied in the data acquisition. First, a shallow time window of the data (202) where the fixed receiver condition is satisfied is inverted by simultaneous encoded (203) source inversion (205). Then, the deeper time window of the data (208) is inverted by sparse sequential source inversion (209), using the physical properties model from the shallow time window (206) as a starting model (207). Alternatively, the shallow time window model is used to simulate missing far offset data (211) producing a data set satisfying the stationary receiver assumption, after which this data set is source encoded (212) and inverted by simultaneous source inversion (214).

Abstract: Method for reducing the time needed to perform geophysical inversion by using simultaneous encoded sources in the simulation steps of the inversion process. The geophysical survey data are prepared by encoding (3) a group of source gathers (1), using for each gather a different encoding signature selected from a set (2) of non-equivalent encoding signatures. Then, the encoded gathers are summed (4) by summing all traces corresponding to the same receiver from each gather, resulting in a simultaneous encoded gather. (Alternatively, the geophysical data are acquired from simultaneously encoded sources.) The simulation steps needed for inversion are then calculated using a particular assumed velocity (or other physical property) model (5) and simultaneously activated encoded sources using the same encoding scheme used on the measured data. The result is an updated physical properties model (6) that may be further updated (7) by additional iterations.

Abstract: The invention includes a method for reducing noise in migration of seismic data, particularly advantageous for imaging by simultaneous encoded source reverse-time migration (SS-RTM). One example embodiment includes the steps of obtaining a plurality of initial subsurface images; decomposing each of the initial subsurface images into components; identifying a set of components comprising one of (i) components having at least one substantially similar characteristic across the plurality of initial subsurface images, and (ii) components having substantially dissimilar characteristics across the plurality of initial subsurface images; and generating an enhanced subsurface image using the identified set of components. For SS-RTM, each of the initial subsurface images is generated by migrating several sources simultaneously using a unique random set of encoding functions. Another embodiment of the invention uses SS-RTM for velocity model building.

Abstract: Provided is a method for processing seismic data. One exemplary embodiment includes the steps of obtaining a plurality of initial subsurface images; decomposing each of the initial subsurface images into components; identifying a set of components comprising one of (i) components having at least one substantially similar characteristic across the plurality of initial subsurface images, and (ii) components having substantially dissimilar characteristics across the plurality of initial subsurface images; and generating an enhanced subsurface image using the identified set of components. Each of the initial subsurface images is generated using a unique random set of encoding functions.