Ray-based seismic modeling methods can be applied at various stages of the exploration and production process. The standard ray method has several advantages, e.g., computational efficiency and the possibility of simulating propagation of elementary waves. As a high-frequency approximation, the method also has a number of limitations, particularly with respect to a lack of amplitude reliability in the presence of rapid changes of the model functions representing elastic parameters and interfaces. Given the objective of improving the applicability of the standard ray method, we present a strategy that does not require specific extension to finite frequencies. Instead, we define each ray-based process as an element of a system that, as a composite process, is able to obtain better results than the ray-based process applied alone. Other elements of the composite process can be finitedifference modeling or numerical solutions for surface and volume integrals, which are basic constituents of Kirchhoff modeling and imaging. We also include among the process elements recently developed techniques for simulating the migration amplitude on a target reflector and in a local volume, e.g., a reservoir zone. The model is decomposed according to its complexity into volume elements, surface elements, or a combination. The composite process consists of a specified interaction between process elements and model elements, which fits well with the philosophy of modern software design. Model elements that will be exposed to ray-tracing algorithms may need appropriate preparation, e.g., smoothing and resampling. We demonstrate specifically, in a tutorial example, that simulating the seismic response from a reflector by ray-based composite processes can yield better results than applying standard ray tracing alone.

Free-surface-related multiples in marine seismic data are commonly attenuated using adaptive subtraction of the predicted multiple energy. An alternative method, based on deconvolution of the upgoing wavefield by the downgoing wavefield, was previously applied to ocean-bottom data. We apply the deconvolution method to towed-streamer data acquired in an over/under configuration. We also use direct arrival deconvolution that results in source wavelet designature only, as a benchmark to verify the full multiple deconvolution result. Detailed synthetic data analysis, including sensitivity tests, explains each data processing step and its effects on the final result. We then apply this verified preprocessing sequence to field data from the Kristin area of the North Sea, with a focus on the direct arrival prediction using the near-field hydrophone method. Prestack evaluation of the results shows that the method applied to the field data provides designature, source-side deghosting, and attenuation of multiples. We show comparable stacked results from our method and from 2D iterative surface-related multiple elimination. The workflow has the benefit that it does not require an adaptive subtraction step or iterative application. However, an accurate direct arrival prediction is essential for the successful application of the method. This prediction is obtained using near-field hydrophone measurements that can be recorded with some commercial acquisition systems.

No abstract

No abstract

No abstract