| Summary: | We attempt to address two issues in seismic data processing: 1) quantifying the various forms of error that enter into the seismic data processing work-flow and relating them to uncertainty on imaged structures; and, 2) the data fusion problem, i.e. combining different sources of information, each related to seismic velocity. To begin addressing these issues a synthetic model was generated consisting of 4 tilted layers (3 interfaces), each with a different isotropic P-wave velocity. A synthetic well log was extracted from this model to be incorporated later. Synthetic shot gathers were also created. Following the standard seismic processing work-flow, stacking velocities were estimated. Uncertainty on these velocities was incorporated by under- and over-picking the velocities by ±10% and examining the effects on the final image. The stacking velocity information was then converted to interval velocity and fused with the well velocity information. Along with the under-, over-, and best picked velocities, realizations of the velocity field were created via geostatistical methods according to an assumed correlation structure. By further applying time migration and time to depth conversion, equiprobable realizations of the subsurface structure were generated along with upper and lower bounds on their locations. The realizations honor all the existing data sets and give a visual representation of the uncertainty on the spatial location of imaged structures.
|
|---|