Xiao-Bi Xie University of California at Santa Cruz Sanya , China July 24-28, 2011

1. Sensitivity kernels for finite-frequency signals: Applications in migration velocity updating and tomography Xiao-Bi Xie University of California at Santa Cruz Sanya, China July 24-28, 2011

2. Outline • A brief introduction • Data domain vs. depth domain • Sensitivity Kernel for Migration Velocity Analysis • The Inversion System • Velocity model partitioning and sensitivity kernel storage • Numerical Result • Conclusions

3. Outline • A brief introduction • Data domain vs. depth domain • Sensitivity Kernel for Migration Velocity Analysis • The Inversion System • Velocity model partitioning and sensitivity kernel storage • Numerical Result • Conclusions

7. df =arg(u/u0) u u0 arg(u/u0) = imag(U/u0) U Imag U / u0 u / u0 u0/ u0 Real

10. In applied seismology Huge data size. Efficiency is crucial. Suggested methods could be one-way propagator or Gaussian beam method. Complex background models. The velocity perturbations overlapped on the initial model are large (some times are more than 100%). Including not only transmitted observations, where the information is from the surface data, but also reflection type observation, where the information is collected in image domain.

11. Outline • A brief introduction • Data domain vs. depth domain • Sensitivity Kernel for Migration Velocity Analysis • The Inversion System • Velocity model partitioning and sensitivity kernel storage • Numerical Result • Conclusions

12. Complexity in data domain ▼ A synthetic shot record. The shot is located above relatively complicated structures. There are many complicated features in this synthetic section. Simplicity in depth domain

13. Outline • A brief introduction • Data domain vs. depth domain • Sensitivity Kernel for Migration Velocity Analysis • The Inversion System • Velocity model partitioning and sensitivity kernel storage • Numerical Result • Conclusions

14. Data Source Back project to modify the velocity Target Measuring incoherence in image The incoherence information are RMOs from different common image gathers. The methods that converting the RMO into velocity corrections. • Offset index CIG • Shot index CIG • Angle index CIG • Parameterized semblance • Ray-based tomography • Wave-equation based inversion Migration velocity updating

15. How the RMO sense the velocity perturbation: --- Direct measurements

16. The actual sensitivity map for a shot image (how the depth image senses the V-model error). To generate this map, we use an velocity error patch to scan the model. At each location, we conduct a migration and measure the RMO from the depth image. The RMOs are then presented in the model to show the sensitivity of the depth image to the velocity error. The sensitivity map is very complex. A positive error can generate either positive or negative RMOs; the sensitivity area is much broader than the ray based theory predicted. Our goal is to derive theoretical equations to express this sensitivity map and use it for velocity updating.

17. Direct measured sensitivity maps for shots at different locations in 2D SEG/EAGE salt model

18. Source side kernel Receiver side kernel Source Image point Imaginary source

19. The GB Green’s functions used to construct the sensitivity kernel for migration velocity analysis. (a) Down-going Green’s function , (b) up-going Green’s function , and (c) Green’s function .

20. Comparison of different kernels for a shot gather The sketch of a ray-based kernel The sensitivity kernel calculated using the finite-frequency theory The actual sensitivity map directly measured from migration

21. Outline • A brief introduction • Data domain vs. depth domain • Sensitivity Kernel for Migration Velocity Analysis • The Inversion System • Velocity model partitioning and sensitivity kernel storage • Numerical Result • Conclusions

22. ∫ dV = δv/v RMO Velocity model error Sensitivity kernel

23. Comparison between inversions using a finite-frequency sensitivity kernel and a ray kernel Sketch illustrating the relative residual moveout measurement from a pair of shots The differential sensitivity kernel for a pair of shots. Note the complexity and volumetric distribution of a finite-frequency kernel The ray based kernel for a pair of shots. Note the sensitivity distribution is unrealistic and the uneven ray distribution can cause singularities in inversions.

24. Differential RMO Differential kernel

25. Outline • A brief introduction • Data domain vs. depth domain • Migration Velocity Analysis • The Inversion System • Velocity model partitioning and sensitivity kernel storage • Numerical Result • Conclusions

26. A 5-layer velocity model used to demonstrate the migration velocity analysis.

27. How to partition the model? How to store huge amount of kernels?

29. Actual output and stored kernels Unknown perturbation at cell corners Parameter matrix

30. Stored parameter a1. The 4 groups of kernels are for 4 reflectors; the horizontal coordinate is for different image points and the vertical coordinate is for different sources. Model grid size 10m x 10m Cell size 500m x 500m 31shot x 31 imaging point x 4 reflectors, 32x10 cells spend 286Mb.

31. Outline • A brief introduction • Data domain vs. depth domain • Sensitivity Kernel for Migration Velocity Analysis • The Inversion System • Velocity model partitioning and sensitivity kernel storage • Numerical Result • Conclusions

32. migration velocity updating process (1) Conduct the migration using an initial model. (2) Calculate the RMOs from the shot-index CIGs. (3) Pick the reflector position from the initial depth image. (4) Use the initial model and reflector locations to calculate sensitivity kernels. (5) Input the RMOs and the sensitivity kernels to the inversion system to do the tomography. (6) Use the inverted errors to update the initial model and use it for the next iteration.

33. A 5-layer velocity model used to demonstrate the migration velocity analysis.

34. Comparison between the theoretically calculated kernels (left column) and actually measured sensitivity maps (right column). From top to bottom are for different reflectors.

35. Coverage of sensitivity kernels in the model. Panels (a) to (d) are kernel coverage for image points on the 4 reflectors. Panel (e) is the coverage from all kernels. Shown here is the summed positive parameter FK1.

36. Velocity models in updating process, with (a) initial model and (b) model after two iterations.

37. Depth image improved in the velocity updating process. (a) Image calculated using the initial model and (b) image calculated using the updated velocity model.

38. CIGs before and after the velocity updating, with (a) CIGs in the initial model and (b) CIGs in the updated velocity model.

39. Outline • A brief introduction • Data domain vs. depth domain • Sensitivity Kernel for Migration Velocity Analysis • The Inversion System • Velocity model partitioning and sensitivity kernel storage • Numerical Result • Conclusions

40. Summary • Based on the finite-frequency sensitivity theory, we present a migration velocity analysis method. The new approach is a wave-equation based method which naturally incorporates the wave phenomena and is best teamed with the wave-equation based migration for velocity analysis. • The finite-frequency sensitivity kernel is used to link the observed shot gather RMO with the errors in the migration velocity model. Angle domain decomposition is not required. • (3) We developed method to calculate the broadband sensitivity kernel in complex velocity models and for irregular reflectors.

41. Summary (continues) (4) A new velocity model partitioning approach is tested. This method partitions the model into small cells and uses interpolation function to represent the velocity model within cells. (5) To store the sensitivity kernels, we use interpolation functions as basis and expanded kernels to these basis. Thus we only need to store the expansion coefficients. The accuracy of the kernel is adaptive to the required accuracy of the velocity model. In this way, we significantly reduce the storage space of sensitivity kernels while without losing the required accuracy.

42. Summary (continues) (6) Using this approach, we demonstrate the velocity modelupdating. The updated velocity model improves the depth image by both flattened the common image gather and bring the image to the original location of reflectors.

43. EOF