Sensitivity kernal-based migration velocity analysis in 3D anisotropic media

• US 8,830,788 B2
• Filed: 02/24/2011
• Issued: 09/09/2014
• Est. Priority Date: 02/24/2011
• Status: Active Grant

• Current Assignee: Landmark Graphics Corporation (Halliburton Company)
• Original Assignee: Landmark Graphics Corporation (Halliburton Company)
• Inventors: Xia, Fan, Ren, Yiqing, Jin, Shengwen
• Primary Examiner(s): Ratcliffe, Luke
• Assistant Examiner(s): Baghdasaryan, Hovhannes
• Application Number:
  US13/033,768
• Publication Number:
  US 20120218861A1
• Time in Patent Office: 1,293 Days
• Field of Search: None
• US Class Current: 367/50
• CPC Class Codes: G01V 1/282 Application of seismic models, synthetic seismograms
• CPC Class Codes: G01V 1/301 for determining seismic cross-sections or geostructures
• CPC Class Codes: G01V 1/303 for determining velocity profiles or travel times
• CPC Class Codes: G01V 11/00 Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
• CPC Class Codes: G01V 2210/47 Slowness, e.g. tau-pi
• CPC Class Codes: G01V 2210/51 Migration
• CPC Class Codes: G01V 2210/6222 Velocity; travel time
• CPC Class Codes: G01V 2210/626 with anisotropy
• CPC Class Codes: G01V 2210/64 Geostructures, e.g. in 3D data cubes

Abstract

Seismic imaging systems and methods that employ sensitivity kernel-based migration velocity analysis in 3D anisotropic media may demonstrate increased stability and accuracy. Survey data analysts employing the disclosed systems and methods are expected to provide better images of the subsurface and be better able to identify reservoirs and deposits for commercial exploitation. Certain embodiments migrate seismic survey data with an anisotropic velocity model to obtain common angle image gathers. These gathers are processed to obtain depth residuals along one or more horizons. Angle-domain sensitivity kernels are used to convert the depth residuals into velocity errors, which are then used to refine the velocity model. A user is then able to view a representation of the subsurface structure determined in part from the refined velocity model. Multiple iterations may be needed for the velocity model to converge. The velocity model may be a layered to have constant velocity values between formation boundaries.

Claims

1. A seismic imaging method that comprises:

migrating seismic survey data with an anisotropic velocity model to obtain common angle image gathers;

processing the common angle image gathers to obtain depth residuals along one or more horizons;

determining subsurface opening angle sensitivity kernels for the horizons independent of shot location and based on the anisotropic velocity model and a Green'"'"'s function at an image point;

inverting the depth residuals using the determined sensitivity kernels to obtain velocity errors;

applying updates to the velocity model based on the obtained velocity errors; and

displaying a representation of the survey region to a user based at least in part on the updated velocity model.

2. The method of claim 1, further comprising re-migrating the seismic survey data using the updated velocity model.

3. The method of claim 2, further comprising repeating said processing, determining, inverting, applying, and re-migrating operations until a velocity model convergence is achieved.

4. The method of claim 1, wherein said processing comprises:

stacking the common image gathers to identify said one or more horizons; and

determining local dip angles based on the one or more horizons.

5. The method of claim 4, wherein said processing further comprises:

measuring similarity between image gathers at different assumed velocity errors to obtain the depth residuals Δ

z.

6. The method of claim 5, wherein said measuring similarity comprises determining semblance of the gathers, and wherein the depth residuals correspond to the assumed velocity error that yields the highest semblance.

7. The method of claim 1, wherein said displaying includes determining a layered velocity model having constant velocity values between horizons.

8. A seismic imaging system that comprises:

at least one storage device that stores traces from a seismic survey of a given survey region;

a memory that stores seismic imaging software; and

at least one processor coupled to the memory to execute the seismic imaging software, wherein the software causes the at least one processor to;

migrate seismic survey data with an anisotropic velocity model to obtain image gathers;

process the image gathers to obtain depth residuals along one or more horizons;

determine subsurface opening angle sensitivity kernels for the horizons independent of shot location and based on the anisotropic velocity model and a Green'"'"'s function at an image point;

invert the depth residuals using the determined sensitivity kernels to obtain velocity errors;

apply updates to the velocity model based on the obtained velocity errors; and

display a representation of the survey region to a user based at least in part on the updated velocity model.

9. The system of claim 8, wherein the software further causes the at least one processor to re-migrate the seismic survey data using the updated velocity model.

10. The system of claim 9, wherein the software further causes the at least one processor to repeat said processing, determining, inverting, applying, and re-migrating operations until a velocity model convergence is achieved.

11. The system of claim 8, wherein as part of said processing operation, the software causes the at least one processor to:

stack the image gathers to identify said one or more horizons; and

determine local dip angles based on the one or more horizons.

12. The system of claim 11, wherein as part of said processing operation, the software further causes the at least one processor to:

measure similarity between image gathers at different assumed velocity errors to obtain the depth residuals Δ

z.

13. The system of claim 12, wherein as part of said measuring, the software causes the processor to determine semblance of the image gathers, and wherein said depth residuals correspond to the assumed velocity error that yields the highest semblance.

14. The system of claim 8, wherein the software further causes the processor to determine a layered velocity model having constant velocity values between horizons.

15. The system of claim 8, wherein the velocity model is 3D.

16. A seismic imaging method that comprises:

migrating seismic data with an anisotropic velocity model to obtain image gathers;

processing the image gathers to obtain depth residuals along one or more horizons;

determining subsurface opening angle sensitivity kernels for the horizons independent of shot location and based on the anisotropic velocity model and a Green'"'"'s function at an image point;

employing the determined sensitivity kernels to convert the depth residuals into velocity errors;

applying updates to the velocity model based on the velocity errors; and

displaying a representation of the survey region to a user based at least in part on the updated velocity model.

17. The method of claim 16, wherein the velocity model is 3D.

18. The method of claim 16, wherein said processing comprises:

stacking the common image gathers to identify said one or more horizons; and

determining local dip angles based on the one or more horizons.

19. The method of claim 18, wherein said processing further comprises:

measuring similarity between image gathers at different assumed velocity errors to obtain the depth residuals Δ

z.

20. The method of claim 16, further comprising determining a layered velocity model having constant velocity values between horizons.