Simultaneous Source Encoding and Source Separation As A Practical Solution For Full Wavefield Inversion

• US 20120073825A1
• Filed: 09/09/2011
• Published: 03/29/2012
• Est. Priority Date: 09/27/2010
• Status: Active Grant

• Current Assignee: Carey Marcinkovich, Jerome R. Krebs, Partha S. Routh, Ramesh Neelamani, Spyridon Lazaratos, Sunwoong Lee
• Original Assignee: Carey Marcinkovich, Jerome R. Krebs, Partha S. Routh, Ramesh Neelamani, Spyridon Lazaratos, Sunwoong Lee
• Inventors: Routh, Partha S., Lee, Sunwoong, Neelamani, Ramesh, Krebs, Jerome R., Lazaratos, Spyridon, Marcinkovich, Carey
• Granted Patent: US 8,775,143 B2
• Time in Patent Office: Days
• US Class Current: 166/369
• CPC Class Codes: G01V 1/28 Processing seismic data, e.g. for interpretation or for event detection G01V1/48 takes precedence
• CPC Class Codes: G01V 2210/67 Wave propagation modeling

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.

Claims

1. A computer-implemented method for simultaneous inversion of measured geophysical data from multiple encoded sources to determine a physical properties model for a subsurface region, the measured geophysical data resulting from a survey for which the fixed-receiver assumption of simultaneous encoded-source inversion may not be valid, said method comprising using a computer to perform simultaneous-source separation to lessen any effect of the measured geophysical data'"'"'s not satisfying the fixed-receiver assumption, wherein a data processing step coming after a simultaneous source separation acts to conform model-simulated data to the measured geophysical data for source and receiver combinations that are missing in the measured geophysical data.

2. The method of claim 1, wherein the method comprises:

constructing an initial physical properties model and using it to simulate synthetic data corresponding to the measured geophysical data, wherein sources and/or receivers in the simulation are encoded and simulated simultaneously, thereby generating simulated encoded data;

separating the simulated encoded data according to source, resulting in simulated sequential-source data;

processing the simulated sequential-source data to make them conform to the measured geophysical data for source and receiver combinations that are missing in the measured geophysical data;

computing differences, called data residuals, between the processed simulated sequential source data and the measured geophysical data for source and receiver combinations that are present in the measured geophysical data;

encoding the data residuals using the same or different encoding as was used in generating the simulated encoded data;

using the encoded data residuals to compute an update to the initial physical properties model; and

combining the update with the initial model to form an updated physical properties model.

3. The method of claim 2, wherein the encoded data residuals are used to compute a gradient of an objective function, said gradient being with respect to parameters of the physical properties model, which is then used in a selected optimization scheme to update the physical properties model.

4. The method of claim 2, further comprising iterating the method steps at least one time, using the updated physical properties model from one iteration as the initial model for a next iteration.

5. The method of claim 4, wherein the sources are encoded using a selected set of encoding functions, and wherein a different set of encoding functions is selected for at least one of the iterations.

6. The method of claim 2, further comprising at any time before the computing differences step:

encoding the measured geophysical data, using the same or different source encoding as was used in generating the simulated encoded data;

separating the encoded measured data according to source using a source separation algorithm that was also used in separating the simulated encoded data according to source; and

using the separated measured data in computing the data residuals.

7. The method of claim 2, further comprising after the separating step and before the computing differences step, removing source-receiver locations from the simulated sequential-source data to correspond to source-receiver locations with missing data in the measured geophysical data.

8. The method of claim 7, wherein the removing source-receiver locations from the simulated sequential-source data is accomplished by one or more data processing steps from a group consisting of offset muting i.e. removing of near and far offsets;

muting based on type of arrivals including transmission versus reflection component of data, and removing the data'"'"'s elastic component where forward simulation is based on an acoustic assumption.

9. The method of claim 2, wherein each gather of the simulated encoded data is encoded by temporally convolving all traces from the gather with an encoding signature selected for the gather.

10. The method of claim 1, wherein the method comprises:

constructing an initial physical properties model and using it to simulate synthetic data corresponding to the measured geophysical data, wherein sources and/or receivers in the simulation are encoded using selected encoding functions and simulated simultaneously, thereby generating simulated encoded data;

separating the simulated encoded data according to source and/or receiver using a source or receiver separation algorithm, resulting in simulated sequential source and/or receiver data;

processing the simulated sequential source and/or receiver data to mute those data, referred to as the missing data, corresponding to source and receiver combinations that are missing in the measured geophysical data;

encoding the missing data, using encoding functions that are the same or different than said selected encoding functions;

encoding the measured geophysical data, using encoding functions that are the same as used in the encoding the missing data, and subtracting them, according to source and receiver position, from the simulated encoded data, thereby generating encoded data residuals, which will be incorrect due to the missing source and receiver combinations in the measured geophysical data;

subtracting the encoded missing data from the encoded data residuals to generate corrected encoded data residuals;

using the corrected encoded data residuals to compute an update to the initial physical properties model; and

combining the update with the initial model to form an updated physical properties model.

11. The method of claim 10, further comprising iterating the method steps at least one time, using the updated physical properties model from one iteration as the initial model for a next iteration.

12. The method of claim 11, wherein the encoding functions or signatures are changed for at least one of the iterations.

13. The method of claim 10, wherein fewer than all source and/or receiver locations in the measured geophysical data are used in the generating simulated encoded data in order to improve accuracy of the source or receiver separation algorithm by making the separation well-posed, but not reducing sampling of the measured geophysical data to an extent that causes aliasing.

14. The method of claim 10, wherein the missing data occur at near and far ends of the measured geophysical data'"'"'s offset range.

15. The method of claim 10, wherein each gather of the simulated encoded data is encoded by temporally convolving all traces from the gather with an encoding signature selected for the gather.

16. The method of claim 10, wherein the encoding the measured geophysical data is accomplished by obtaining gathers of data from a geophysical survey in which data are acquired from a plurality of simultaneously operating, uniquely encoded source devices.

17. The method of claim 1, wherein the method comprises:

constructing an initial physical properties model and using it to simulate synthetic data corresponding to the measured geophysical data, wherein sources and/or receivers in the simulation are encoded, using selected encoding functions, and simulated simultaneously, thereby generating simulated encoded data;

using said selected encoding functions to encode the measured geophysical data and subtracting them, according to source and receiver position, from the simulated encoded data, thereby generating encoded data residuals, which will include incorrect values due to the missing source and receiver combinations in the measured geophysical data;

separating the encoded data residuals according to source and/or receiver using a source or receiver separation algorithm, resulting in sequential residual data;

processing the sequential residual data to estimate simulated data, referred to as missing data, corresponding to source and receiver combinations that are missing in the measured geophysical data;

encoding the missing data, using said selected encoding functions;

subtracting the encoded missing data from the encoded data residuals to generate corrected encoded data residuals;

using the corrected encoded data residuals to compute an update to the initial physical properties model; and

combining the update with the initial model to form an updated physical properties model.

18. The method of claim 1, wherein the geophysical data are seismic data, of which a full wavefield is inverted.

19. The method of claim 1, wherein a plurality of encoded data sets are used to perform the simultaneous-source separation.

23. A method for producing hydrocarbons from a subsurface region, comprising:

performing a geophysical survey of the subsurface region resulting in measured geophysical data;

processing the measured geophysical data on a computer by a method of claim 1 to generate a physical properties model of the subsurface region;

assessing hydrocarbon potential of the subsurface region using the physical properties model; and

drilling a well into the subsurface region based at least in part on the assessment of hydrocarbon potential, and producing hydrocarbons from the well.

20. A computer program product, comprising a non-transitory computer usable medium having a computer readable program code embodied therein, said computer readable program code adapted to be executed to implement a method for full-wavefield inversion of measured geophysical data to determine a physical properties model for a subsurface region, said method comprising:

constructing or inputting an initial physical properties model and using it to simulate synthetic data corresponding to the measured geophysical data, wherein sources and/or receivers in the simulation are encoded and simulated simultaneously, thereby generating simulated encoded data;

separating the simulated encoded data according to source, resulting in simulated sequential-source data;

computing differences, called data residuals, between the simulated sequential source data and the corresponding measured geophysical data;

encoding the data residuals according to source using the same encoding as was used in generating the simulated encoded data;

using the encoded data residuals to compute an update to the initial physical properties model; and

combining the update with the initial model to form an updated physical properties model.

21. A computer program product, comprising a non-transitory computer usable medium having a computer readable program code embodied therein, said computer readable program code adapted to be executed to implement a method for full-wavefield inversion of measured geophysical data to determine a physical properties model for a subsurface region, said method comprising:

Constructing or inputting an initial physical properties model and using it to simulate synthetic data corresponding to the measured geophysical data, wherein sources and/or receivers in the simulation are encoded using selected encoding functions and simulated simultaneously, thereby generating simulated encoded data;

separating the simulated encoded data according to source and/or receiver using a source or receiver separation algorithm, resulting in simulated sequential source and/or receiver data;

processing the simulated sequential source and/or receiver data to mute those data, referred to as the missing data, corresponding to source and receiver combinations that are missing in the measured geophysical data;

encoding the missing data, using encoding functions that are the same or different than said selected encoding functions;

encoding the measured geophysical data, using encoding functions that are the same as used in the encoding the missing data, and subtracting them, according to source and receiver position, from the simulated encoded data, thereby generating encoded data residuals;

subtracting the encoded missing data from the encoded data residuals to generate corrected encoded data residuals;

using the corrected encoded data residuals to compute an update to the initial physical properties model; and

combining the update with the initial model to form an updated physical properties model.

22. A computer program product, comprising a non-transitory computer usable medium having a computer readable program code embodied therein, said computer readable program code adapted to be executed to implement a method for full-wavefield inversion of measured geophysical data to determine a physical properties model for a subsurface region, said method comprising:

Constructing or inputting an initial physical properties model and using it to simulate synthetic data corresponding to the measured geophysical data, wherein sources and/or receivers in the simulation are encoded, using selected encoding functions, and simulated simultaneously, thereby generating simulated encoded data;

using said selected encoding functions to encode the measured geophysical data and subtracting them, according to source and receiver position, from the simulated encoded data, thereby generating encoded data residuals;

separating the encoded data residuals according to source and/or receiver using a source or receiver separation algorithm, resulting in sequential residual data;

processing the sequential residual data to estimate simulated data, referred to as missing data, corresponding to source and receiver combinations that are missing in the measured geophysical data;

encoding the missing data, using said selected encoding functions;

subtracting the encoded missing data from the encoded data residuals to generate corrected encoded data residuals;

using the corrected encoded data residuals to compute an update to the initial physical properties model; and

combining the update with the initial model to form an updated physical properties model.