next up previous print clean
Next: Imaging of GPR data Up: Diffraction-focusing WEMVA example Previous: Diffraction-focusing WEMVA example

Delineation of rough salt bodies

BPAITdat shows the zero-offset data we use for velocity analysis to delineate the top of the rough salt body. The section contains a large number of diffractors, whose focusing allows us to constrain the overburden velocity model.

 
BPAITdat
BPAITdat
Figure 38
Zero-offset synthetic data used for focusing migration velocity analysis.
view burn build edit restore

BPAITimg1(a) depicts the starting velocity model, and BPAITimg1(b) depicts the initial image obtained by zero-offset migration. The starting velocity is a typical Gulf of Mexico v(z) function hanging from the water bottom. Uncollapsed diffractions are visible at the top of the salt, indicating that the velocity in the overburden is not accurate. Such defocusing also hampers our ability to pick accurately the top of the salt and, therefore, degrades imaging at depth.

 
BPAITimg1
BPAITimg1
Figure 39
Zero-offset migrated image for the synthetic data in BPAITdat: velocity model (a), and migrated image (b). Migration using the initial v(z) velocity model.


[*] view burn build edit restore

As for the preceding synthetic example, I run residual migration on the background image (BPAITimg1). BPAITsrm shows this image after residual migration with various velocity ratios (storm). From top to bottom, the ratios are: 1.04, 1.00, 0.96, 0.92, 0.88. At $\rho=1.00$ we recover the initial image. Different parts of the image come into focus at different values of the velocity ratio.

 
BPAITsrm
BPAITsrm
Figure 40
Residual migration applied to the image migrated with the initial velocity model, BPAITimg1. From top to bottom, the images correspond to the ratios $\rho=1.04, 1.00, 0.96, 0.92, 0.88$.
[*] view burn build edit restore

BPAITpck(a) shows the picked velocity ratios at various locations in the image. The white background corresponds to picked $\Delta \rho=0$, and the gray shades correspond to $\Delta \rho$ between and 0.08. BPAITpck(b) shows a map of the weights ($\Wop$) associated to each picked value. The white background corresponds to $\Wop=0$, indicating low confidence in the picked values, and the dark regions correspond to $\Wop=1$, indicating high confidence in the picked values. In this example, we disregard regions where we did not pick any diffractions. All other regions receive an arbitrary ratio value ($\rho=1.0$), but also a low weight such that they do not contribute to the inversion. Exceptions include the water bottom, for which we assign a high weight of the picked ratio $\rho=1.0$,and a few other reflectors for which we did not have any diffraction focusing information.

 
BPAITpck
BPAITpck
Figure 41
Residual migration picks (a) and the associated confidence weights (b).
[*] view burn build edit restore

BPAITmva(a) shows the slowness perturbation obtained after 20 iterations of zero-offset inversion from the image perturbation in BPAITmva(b). The image perturbation is non-zero only in the regions where we had diffractions we could pick, as indicated by BPAITpck. The smooth slowness perturbation is further constrained by the regularization operator we use, which is a simple Laplacian penalizing the rough portions of the model.

 
BPAITmva
BPAITmva
Figure 42
Slowness perturbation (a), derived from an image perturbation (b) derived from the background image in BPAITimg1 and the velocity ratio picks in BPAITpck.
[*] view burn build edit restore

BPAITimg2(a) shows the updated slowness model and BPAITimg2(b) shows the zero-offset migrated image corresponding to the updated model. Most of the diffractions at the top of the salt have been collapsed, and the rough top of the salt can be easily picked. The diffractions corresponding to the salt bodies at x=2000-4000 ft, z=3500 ft are not fully collapsed, indicating that another nonlinear iteration involving residual migration and picking might be necessary.

 
BPAITimg2
BPAITimg2
Figure 43
Zero-offset migrated image for the synthetic data in BPAITdat: velocity model (a), and migrated image (b). Migration using the updated velocity.
[*] view burn build edit restore

Finally, BPAITpre shows prestack migrated images using the initial velocity model (a) and the one updated using zero-offset focusing (b). The top panels depict stacks, and the bottom panels depict angle-domain common image gathers (adcig). The ADCIGs show substantial bending after migration with the initial velocity, but they are mostly flat after migration with the updated velocity, although none of the moveout information has been used for velocity update. BPAITcig shows two ADCIGs at x=-2350 ft from the images obtained with the initial velocity model (a) and the updated velocity model (b). The ADCIG in panel (a) corresponds to a notch in the top of the salt and is complicated to use for velocity analysis. However, after migration with the updated velocity model, panel (b), the ADCIG is much simpler, and the small residual moveouts can be picked for velocity updates.

 
BPAITpre
BPAITpre
Figure 44
Prestack migrated images using the initial velocity model (a) and the updated velocity model (b). The top panels depict image stacks and the bottom panels depict angle-domain common image gathers.


view burn build edit restore

 
BPAITcig
Figure 45
Angle-domain common image gather obtained after migration with the initial velocity model (a) and the updated velocity model (b).

BPAITcig
view burn build edit restore

A comparison of BPAITimg2(b) and BPAITpre(b) shows a potential limitation of this technique in the presence of prismatic waves (8). Both images are obtained with the same velocity, the first one with zero-offset data and the second one with prestack data. The imaging artifacts visible at the bottom of the deep canyons at the top of the salt in BPAITimg2(b) are created by prismatic waves that are not properly imaged from zero-offset data. Prismatic waves are better (though not perfectly) handled by full prestack migration, and thus the artifacts are not visible in the prestack-migrated image shown in BPAITpre(b). Since these artifacts resemble uncollapsed diffractions, they may mislead the analysis of the residual migrated images and be interpreted as symptoms of velocity inaccuracies.


next up previous print clean
Next: Imaging of GPR data Up: Diffraction-focusing WEMVA example Previous: Diffraction-focusing WEMVA example
Stanford Exploration Project
11/4/2004