Patent Application
20100302906
The present invention relates generally to geophysical exploration and in particular to a method of migration and inversion of seismic data using multiple reflections in such signals or data to obtain characteristics of a subsurface region of interest.
Reverse time migration (RTM) has been applied to imaging complex structures for oil and gas exploration and development. Compared to one-way prior art imaging methods, prior art RTM is based on solving the two-way wave equation and can propagate wavefields in all directions. RTM also preserves propagation amplitude accurately. These advantages over one-way imaging often result in significantly improved images of complex structures, especially when using wide azimuth data.
The current practice of RTM is still limited to using data with free-surface multiples removed. In this way, RTM is primarily used to focus multiple bounces (so-called prism waves) from the same hard interface such as salt flanks when compared to the one-way imaging methods. In the presence of free-surface related multiples, prior art RTM methods generate spurious events in output images due to imperfect data recording geometry (Mittet, 2002). Similarly, internal multiples can also lead to spurious events based on the same workflow.
The present invention provides methods to mitigate the current limitations in handling multiples and can utilize data more fully in a constructive way.
One embodiment of the present invention includes a method for wavefield-based data processing including the use of free-surface and internal multiples to obtain characteristics of a subsurface region of interest. The method includes obtaining an earth model (for example, the earth model may define velocity, density, and anisotropy) and a migration model (for example, the earth model may define macro-scale migration velocity and anisotropy) related to the subsurface region of interest. The method further includes determining a modeling geometry related to the subsurface region of interest for the earth model and for the migration model, and propagating forward at least one wavefield in the earth model from at least one excitation source obtained from the modeling geometry. The method also includes propagating forward at least one wavefield in the migration model from at least one excitation source obtained from the modeling geometry. The method also includes propagating backward at least one wavefield in the earth model utilizing at least one receiver location obtained from the modeling geometry. The method additionally includes determining at least one composite wavefield from the previous forward propagated source wavefield(s) (accessed in reverse time order through either storage or re-computation) and the backward propagated receiver wavefield(s) from the earth model. The method additionally includes applying imaging conditions to the forward propagated source wavefield (but accessed in reverse time order through either storage or re-computation) from the migration model and the composite wavefield from the earth model, wherein the imaging conditions utilize the multiples present in the composite wavefield to determine characteristics of the subsurface region of interest without generating corresponding spurious events of the multiples.
It is an object of the present invention to provide a method for utilizing multiples to determine characteristics of a subsurface region of interest wherein the multiples include at least one of free-surface multiples and/or internal multiples.
It is an object of the present invention to have embodiments utilizing multiples to obtain characteristics of a subsurface region which can be used for two-way propagation methods, waveform inversion, model building or property estimation.
It is an object of the present invention to have embodiments utilizing multiples to obtain characteristics of a subsurface region of interest in the frequency or wavelet domain.
It is an object of the present invention to utilize wavefields including derivative quantities, such as, but not limited to, residual wavefields.
Another embodiment of the present invention includes a migration or inversion method which includes establishing a data set, an estimated earth model, and a migration model corresponding to an exploration volume. The method also includes setting boundary or initial conditions of wavefield propagation, and propagating wavefields from a source governed by an appropriate wave equation using the earth model. The method further includes propagating wavefields from the source again, using the migration model, and back propagating the measured traces from receivers and concurrently back propagating the earth model-based source wavefields to construct composite wavefields. The method additionally includes applying imaging conditions such as, but not limited to cross correlation to the migration model-based source wavefields and earth model-based composite wavefields to obtain subsurface images or properties.
The present invention differs from prior art methods in that the input seismic used in the present invention doesn't require preprocessing to remove or suppress multiples. If the method of the prior art takes input data without multiples removal, spurious events will be present in final images. In contrast, the present invention can constructively use multiples in the data for imaging and inversion in that artificial transmission or reflection events from multiples are eliminated or largely reduced in the wave extrapolation process to avoid spurious images. As a result, the limited surface acquisition geometry is compensated by utilizing a good estimate of the earth properties to fully utilize two-way wave propagation for various applications.
Although the above-described embodiment, by way of example, requires a good estimate of the true earth model, this condition can be relaxed to various degrees in practice and can be substituted by other approximations to result in equivalent elimination/reduction of artificial events. In addition, any imperfect elimination of spurious events is also an indication of errors in the estimated earth model which can be leveraged to improve model building. Therefore, the present invention can also be used to improve model building and properties estimation.
It should also be appreciated by one skilled in the art that the present invention is intended to be used with a system which includes, in general, an electronic configuration including at least one processor, at least one memory device for storing program code or other data, a video monitor or other display device (i.e., a liquid crystal display) and at least one input device. The processor is preferably a microprocessor or microcontroller-based platform which is capable of displaying images and processing complex mathematical algorithms. The memory device can include random access memory (RAM) for storing event or other data generated or used during a particular process associated with the present invention. The memory device can also include read only memory (ROM) for storing the program code for the controls and processes of the present invention.
As an example, one embodiment of the present invention includes a system configured to perform wavefield-based seismic data processing including utilizing multiples to obtain characteristics of a subsurface region of interest. The system includes a data storage device having computer readable data including an earth model and a migration model related to the subsurface region of interest. The system also includes a processor, configured and arranged to execute machine executable instructions stored in a processor accessible memory for performing a method. The method includes determining a modeling geometry related to the subsurface region of interest for the earth model and for the migration model, and propagating forward at least one wavefield in the earth model from at least one excitation source obtained from the modeling geometry. The method also includes propagating forward at least one wavefield in the migration model from the at least one excitation source obtained from the modeling geometry, and propagating backward at least one wavefield in the earth model utilizing at least one receiver location obtained from the modeling geometry. The method further includes determining at least one composite wavefield from the forward and the backward propagated wavefields from the earth model, and applying imaging conditions to the forward propagated wavefield accessed in reverse time order from the migration model and the composite wavefield from the earth model, wherein the imaging conditions utilize the multiples present in the composite wavefield to determine characteristics of the subsurface region of interest without generating corresponding spurious events of the multiples.
It will also be appreciated that such a system described-above may also include a display device which displays the characteristics of the subsurface region of interest. These and other objects, features, and characteristics of the present invention, as well as the methods of operation and functions of the related elements of structure and the combination of parts and economies of manufacture, will become more apparent upon consideration of the following description and the appended claims with reference to the accompanying drawings, all of which form a part of this specification, wherein like reference numerals designate corresponding parts in the various Figures. It is to be expressly understood, however, that the drawings are for the purpose of illustration and description only and are not intended as a definition of the limits of the invention. As used in the specification and in the claims, the singular form of “a”, “an”, and “the” include plural references unless the context clearly dictates otherwise.
These and other objects, features and advantages of the present invention will become better understood with regard to the following description, pending claims and accompanying drawings where:
RTM is one kind of adjoint state problem. On the one hand, the source wavefield is propagated forward over time and accessed in reverse order through either state recording or re-computation. On the other hand, seismic data are back extrapolated and correlated with the source wavefield at the times when reflections occurred. However, prior art RTM requires that free-surface multiples be removed prior to migration otherwise multiples will be focused into spurious reflections in images.
The present invention provides methods to eliminate or significantly reduce spurious transmissions/reflections which can result in ghost images.
Using the methods in the present invention, free-surface multiple removal is no longer a data preprocessing requirement. Instead, free-surface and internal multiples can be used constructively towards imaging in addition to contributions from primaries. The inclusion of multiples in a constructive way can lead to improved imaging aperture, improved subsurface illumination, and improved solvability of inversion problems.
The above-described method is preferably implemented on either co-processor accelerated architectures, such as Field-Programmable-Gate-Arrays (FPGAs), Graphics-Processing-Units (GPUs), Cells, or general-purpose computers. The present invention provides apparatus and general-purpose computers and/or co-processors programmed with instructions to perform a method for the present invention, as well as computer-readable media encoding instructions to perform a method of the present invention.
An example of a system for performing the present invention is schematically illustrated in
It will be clear to one skilled in the art that the above embodiments may be altered in many ways without departing from the scope of the invention. For example, as is apparent to the skilled artisan, different initial conditions or boundary conditions or a different linear combination of the PDEs in the present invention can be used in modeling and migration as convenient.
While in the foregoing specification this invention has been described in relation to certain preferred embodiments thereof, and many details have been set forth for purpose of illustration, it will be apparent to those skilled in the art that the invention is susceptible to alteration and that certain other details described herein can vary considerably without departing from the basic principles of the invention.
