Methods and systems performing joint VVAz and AVAz inversion

Abstract

Methods for seismic exploration of an underground formation including at least one anisotropic layer perform a joint velocity-variation-with-azimuth, VVAz, and amplitude-variation-with-azimuth, AVAz, inversion using the azimuthal angle stacks to obtain a structural representation of the underground formation. The structural representation is used to generate scenarios for exploiting resources in at least one layer of the underground formation.

Claims

1. A method for seismic exploration of an underground formation including at least one anisotropic layer, the method comprising: receiving seismic data acquired by detecting seismic excitations traveling through the underground formation; grouping the seismic data in azimuthal angle stacks; performing a joint velocity-variation-with-azimuth, VVAz, and amplitude-variation-with-azimuth, AVAz, inversion using the azimuthal angle stacks to obtain a structural representation of the underground formation; and generating scenarios for exploiting resources in at least one layer of the underground formation based on the structural representation.

2. The method of claim 1, wherein the joint VVAz and AVAz inversion yields anisotropy properties consistent with both velocities and amplitudes.

3. The method of claim 1, wherein the performing of the joint VVAz and AVAz includes determining a group velocity and a phase velocity for a layer.

4. The method of claim 3, wherein the phase velocity is determined by rotating a stiffness matrix with the incidence angle and the azimuth angle.

5. The method of claim 3, wherein the group velocity, Vg, is obtained using the phase velocity according to: ${\left(\mathrm{Vg}\left(\right)\right)}^2={\left(v_P\left(\right)\right)}^2+{\left(\frac{v_P}{}\right)}^2$ where is a group angle, is a phase angle, and v.sub.p is the phase velocity.

6. The method of claim 1, wherein the performing of the joint VVAz and AVAz inversion uses different time axes for different azimuthal angle stacks.

7. The method of claim 6, wherein a difference in a two-way time-shift TWT.sub., for an incidence angle, , and an azimuth angle, , is given by ${\mathrm{TWT}}_{,}={\mathrm{TWT}}_0\left(\frac{{\mathrm{Vg}}_0}{\mathrm{Vg}\left(,\right)}-1\right)$ where Vg.sub.0 is a reference group velocity, Vg(,) is a direction-dependent group velocity and TWT.sub.0 is a reference layer thickness expressed as a two-way time.

8. A system for seismic exploration of an underground formation including at least one anisotropic layer, the system comprising: seismic data storing equipment configured store seismic data obtained by injecting seismic excitations into the underground formation, detecting the seismic excitations that travelled through the underground formation, and generating the seismic data based on the detected seismic excitations; a seismic data processing unit configured to group the seismic data in azimuthal angle stacks, to perform a joint velocity-variation-with-azimuth, VVAz, and amplitude-variation-with-azimuth, AVAz, inversion using the azimuthal angle stacks to obtain a structural representation of the underground formation, and to generate scenarios related to exploiting resources in at least one layer of the underground formation using the structural representation.

9. The system of claim 8, wherein the joint VVAz and AVAz inversion yields anisotropy's properties consistent with both velocities and amplitudes.

10. The system of claim 8, wherein the data processing unit includes determining a group velocity and a phase velocity for a layer in the joint VVAz and AVAz inversion.

11. The system of claim 10, wherein the data processing unit determines the phase velocity by rotating a stiffness matrix with the incidence angle and the azimuth angle.

12. The system of claim 11, wherein the data processing unit uses the phase velocity to obtain the group velocity, Vg, according to: ${\left(\mathrm{Vg}\left(\right)\right)}^2={\left(v_P\left(\right)\right)}^2+{\left(\frac{v_P}{}\right)}^2$ where is a group angle, is a phase angle, and v.sub.p is the phase velocity.

13. The system of claim 8, wherein the data processing unit uses different time axes for different azimuthal angle stacks when performing of the joint VVAz and AVAz inversion.

14. The system of claim 13, wherein a difference in a two-way time-shift TWT.sub., for an incidence angle, , and an azimuth angle, , is given by ${\mathrm{TWT}}_{,}={\mathrm{TWT}}_0\left(\frac{{\mathrm{Vg}}_0}{\mathrm{Vg}\left(,\right)}-1\right)$ where Vg.sub.0 is a reference group velocity, Vg(,) is a direction-dependent group velocity and TWT.sub.0 is a reference layer thickness expressed as two-way time.

15. A computer readable medium storing executable codes which, when executed by a processor, make the processor perform a method of seismic exploration of an underground formation including at least one anisotropic layer, the method comprising: receiving seismic data acquired by detecting seismic excitations traveling through the underground formation; grouping the seismic data in azimuthal angle stacks; performing a joint velocity-variation-with-azimuth, VVAz, and amplitude-variation-with-azimuth, AVAz, inversion using the azimuthal angle stacks to obtain a structural representation of the underground formation; and generating scenarios for exploiting resources in at least one layer of the underground formation based on the structural representation.

16. The computer readable medium of claim 15, wherein the joint VVAz and AVAz inversion yields anisotropy properties consistent with both velocities and amplitudes.

17. The computer readable medium of claim 15, wherein the performing of the joint VVAz and AVAz includes determining a group velocity and a phase velocity for a layer.

18. The computer readable medium of claim 17, wherein the group velocity, Vg, is obtained using the phase velocity according to: ${\left(\mathrm{Vg}\left(\right)\right)}^2={\left(v_P\left(\right)\right)}^2+{\left(\frac{v_P}{}\right)}^2$ where is a group angle, is a phase angle, and v.sub.p is the phase velocity.

19. The computer readable medium of claim 15, wherein the performing of the joint VVAz and AVAz inversion uses different time axes for different azimuthal angle stacks.

20. The computer readable medium of claim 19, wherein a difference in a two-way time-shift TWT.sub., for an incidence angle, , and an azimuth angle, , is given by ${\mathrm{TWT}}_{,}={\mathrm{TWT}}_0\left(\frac{{\mathrm{Vg}}_0}{\mathrm{Vg}\left(,\right)}-1\right)$ where Vg.sub.0 is a reference group velocity, Vg(,) is a direction-dependent group velocity and TWT.sub.0 is a reference layer thickness expressed as a two-way time.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate one or more embodiments and, together with the description, explain these embodiments. In the drawings:

(2) FIG. 1 illustrates seismic data of a COCA cube;

(3) FIG. 2 is a graph illustrating amplitude variation with azimuth for a given incidence angle and time sample;

(4) FIG. 3 is a flowchart of a conventional method sequentially performing VVAz and AVAz;

(5) FIG. 4 is a flowchart of a method performing joint VVAz and AVAz, according to an embodiment;

(6) FIG. 5 illustrates 4D time-shifts in a stratigraphic approach;

(7) FIG. 6 illustrates the use of the stratigraphic grid for multi-component data;

(8) FIG. 7 illustrates the stratigraphic approach for azimuthal anisotropy analysis;

(9) FIG. 8 illustrates variation with incidence angle only for the elements in FIG. 7;

(10) FIG. 9 illustrates azimuth angles for the elements in FIG. 7;

(11) FIG. 10 is a graph illustrating phase velocity dependence to azimuth angle for different incidence angles measured from the vertical;

(12) FIG. 11 is a graph illustrating group and phase angles variation with azimuth for different incidence angles measured from the vertical;

(13) FIG. 12 is a graph illustrating the difference between phase velocity and group velocity in an HTI medium as a function of incidence angle measured from the vertical for a given azimuth;

(14) FIG. 13 is a flowchart of a method according to an embodiment; and

(15) FIG. 14 is a schematic diagram of a seismic data processing apparatus according to an embodiment.

DETAILED DESCRIPTION

(16) The following description of the exemplary embodiments refers to the accompanying drawings. The same reference numbers in different drawings identify the same or similar elements. The following detailed description does not limit the invention. Instead, the scope of the invention is defined by the appended claims. The following embodiments are discussed, for simplicity, with regard to the terminology of seismic survey in general, the joint VVAz and AVAz inversion being usable for both land and marine seismic data.

(17) Reference throughout the specification to one embodiment or an embodiment means that a particular feature, structure or characteristic described in connection with an embodiment is included in at least one embodiment of the subject matter disclosed. Thus, the appearance of the phrases in one embodiment or in an embodiment in various places throughout the specification is not necessarily referring to the same embodiment. Further, the particular features, structures or characteristics may be combined in any suitable manner in one or more embodiments.

(18) FIG. 4 is a flowchart of a method performing joint VVAz and AVAz inversion, according to an embodiment. Seismic data 400 is grouped in unflattened azimuthal angle stacks at 410. As previously illustrated in FIG. 1, unflatten data has azimuthal variation of travel-times and amplitudes. Azimuthal normal-moveout correction (performed as part of the VVAz analysis) reveals that the travel-times anisotropy increases with the incidence angle. Both corrected and uncorrected traces reveal the amplitude anisotropy increases with incidence angle. The joint VVAz and AVAz inversion works with the unflatten traces and simultaneously translates/inverts the travel-times and amplitude changes in terms of anisotropy properties.

(19) The unflattened azimuthal angle stacks are then inverted using VVAz and AVAz jointly at 420, yielding anisotropy orientation and magnitude consistent with velocities and amplitudes 424.

(20) The joint VVAz and AVAz inversion uses a stratigraphic grid whose resolution may be as fine as a few ms, thereby greatly enhancing the resolution of the VVAz approach. In the stratigraphic grid can coexist multiple vertical axes (time or depth). The use of stratigraphic grid is explained first in the simpler 4D (i.e., surveys acquired in the same area with long time intervals there-between, the earliest survey being known as base and the later one(s) as monitor(s)), and PP-PS inversion contexts.

(21) In 4D, the multiple time axes allow handling production-induced time-shifts by defining a time axis per vintage without aligning the monitor(s) traces to the base traces. Neglecting compaction effects between vintages (i.e., at constant depth), 4D time-shifts (TWT.sub.4D) can be related to 4D changes in P-velocity (Vp.sub.4D) as illustrated in FIG. 5 and mathematically expressed in formula (1) below, where TWT.sub.B is the time-thickness of the layer for the base vintage, TWT.sub.M=TWT.sub.B+TWT.sub.4D is the time-thickness of the layer for the monitor vintage, Vp.sub.B is the P-velocity of the base and Vp.sub.M=Vp.sub.B+Vp.sub.4D is the P-velocity of the monitor vintage. The vertical arrows indicate that the PP seismic waves are traveling vertically (i.e., considering post-stack time-shifts only).

(22) $\begin{array}{cc}{\mathrm{TWT}}_{4D}={\mathrm{TWT}}_B\left(\frac{{\mathrm{Vp}}_B}{{\mathrm{Vp}}_M}-1\right)& \left(1\right)\end{array}$

(23) In the case of multi-component (PP-PS) data, the stratigraphic grid handles the difference in travel time between PP and PS waves by defining a time axis for each wave mode without squeezing the PS seismic traces to match PP times. Since the depth is the same in the PP and PS domains, similar to formula (1), formula (2) below links the layer thickness in PS times (TWT.sub.PS), in PP times (TWT.sub.PP), and the P- to S-velocity ratio Vp/Vs. Here again, post-stack travel times are considered. This correspondence is illustrated in FIG. 6.

(24) $\begin{array}{cc}{\mathrm{TWT}}_{\mathrm{PS}}=\frac{1}{2}{\mathrm{TWT}}_{\mathrm{PP}}\left(\frac{\mathrm{Vp}}{\mathrm{Vs}}+1\right)& \left(2\right)\end{array}$

(25) The same stratigraphic grid approach can be used for azimuthal anisotropy, but the travel-times considered are now pre-stack as non-vertical propagation that varies with incidence angle and azimuth as illustrated in FIG. 7. The arrows on top of the volume elements A, B and C (the middle element being reference) illustrate the different incidence angles and azimuth angles. The variation with incidence angle only for the elements in FIG. 7 is illustrated in FIG. 8, and the azimuth angles thereof are illustrated in FIG. 9.

(26) Unlike the previous cases (i.e., 4D and PP-PS) where the medium was assumed to be isotropic, an anisotropic medium is now considered. As a result, the phase and group velocities become different. AVAz (amplitude) is related to phase velocities, whereas VVAz (travel-time) is related to group velocities. In order to invert jointly velocities/travel-times and amplitudes, the phase velocities extracted from the amplitudes are converted into group velocities.

(27) FIG. 10 is a graph illustrating the phase velocity as a function of azimuth angle for different incidence angles measured from the vertical. The azimuthal variations increase with incidence angles. The previously cited 1998 Rger article provides a function linking the phase velocity to incidence angle and azimuth. Another form of the function may also be derived by rotating the stiffness matrix with rotation angles corresponding to the incidence angle and azimuth.

(28) Phase and group angles are also different. The phase angle is the angle of the wave-front normal with reference to the symmetry axis, whereas the group or ray angle is the angle along which energy propagates. The article entitled, Weak elastic anisotropy, by Thomsen, L. published in Geophysics, Vol. 51, No. 10 (October 1986), pp. 1954-1966, and the article entitled, Reflection moveout and parameter estimation for horizontal transverse isotropy, by Tsvankin, I. published in Geophysics, Vol. 62, No. 2 (March-April 1997), pp. 614-629, provide illustration of the difference between phase and group angles in the case of transverse isotropy media with vertical axis (VTI) and with horizontal axis respectively. Thomsen's and Tsvankin's articles are incorporated herein by reference in their entirety.

(29) The phase angle between the phase-velocity vector and the symmetry axis of the medium can be obtained from the group angle (angle between the reflected group-velocity vector and the symmetry axis) using the derivative of the phase velocity Vp with regards to the phase angle according to the following formula (3) derived from Tsvankin 1997:

(30) $\begin{array}{cc}\tan \left(-\right)=\frac{{\mathrm{dv}}_p/d}{v_p}.& \left(3\right)\end{array}$

(31) Formula (3) is a non-linear equation easily solved using a few iterations of Newton's method. FIG. 11 shows a comparison of group angles (solid lines) converted to phase angles (dashed lines) in HTI media as a function of azimuth angle for different incidence angles (0, 30, 60 and 90) as labeled.

(32) Then, the group velocity Vg may be obtained using the phase velocity Vp, the group and phase angles ( and respectively) using the following formula (4) derived from Thomsen 1986:

(33) $\begin{array}{cc}{\left(\mathrm{Vg}\left(\right)\right)}^2={\left(v_P\left(\right)\right)}^2+{\left(\frac{v_P}{}\right)}^2& \left(4\right)\end{array}$

(34) FIG. 12 is a graph illustrating the difference between phase velocity and group velocity in an HTI medium. Group velocities are used to model azimuthal travel-times (i.e., pre-stack time-shifts), and phase velocities are used to model the azimuthal amplitudes. The AVAz part of the joint VVAz-AVAz inversion can be solved using existing methods in the azimuthal angle stack domain or more efficient AVAz techniques such as Fourier coefficients and statistical moments.

(35) Similar to formula (1) and (2), the difference of travel-time as a function of incidence angle (i.e. group angle in the incidence plane using 1997 Tsvankin's notation) and azimuth (i.e. the source-receiver azimuth ) TWT.sub., is given by formula (5) below, where Vg.sub.0 is the reference group velocity, Vg(,) is the direction-dependent group velocity and TWT.sub.0 is the reference layer thickness in time (the middle volume element in FIG. 7). The stratigraphic grid therefore has a time axis per incidence angle and azimuth as illustrated in FIGS. 7-9.

(36) $\begin{array}{cc}{\mathrm{TWT}}_{,}={\mathrm{TWT}}_0\left(\frac{{\mathrm{Vg}}_0}{\mathrm{Vg}\left(,\right)}-1\right)& \left(5\right)\end{array}$

(37) Joint VVAz and AVAz inversion enhances the results of both VVAz and AVAz simultaneously because each method benefits from the other. Several time axes are defined in the stratigraphic grid, their azimuthal and incidence angle variations being controlled by the group velocity derived from the phase velocity. This method is formulated for HTI anisotropy but can easily be extended to orthorhombic anisotropy.

(38) FIG. 13 is a flowchart of method 1300 according to an embodiment. Method 1300 includes receiving seismic data acquired by detecting seismic excitations traveling through the underground formation at 1310, grouping the seismic data in azimuthal angle stacks at 1320, and performing a joint velocity-variation-with-azimuth, VVAz, and amplitude-variation-with-azimuth, AVAz, inversion using the azimuthal angle stacks to obtain a structural representation of the underground formation at 1330. Various aspects of steps 1310 and 1320 have been discussed above.

(39) The method further includes generating scenarios for exploiting resources in at least one layer of the underground formation based on the structural representation at 1340. The outcomes of inversions are structural representation of the explored underground formation, which means, in other words, they map properties such as velocity, locating interfaces between different materials, fractures, porosity, etc. These types of information allow evaluation of the presence of hydrocarbons and various minerals and developing exploitation plans (e.g., a path to be followed when drilling to reach the oil or sought-after minerals, to pump liquids when extracting gas by hydraulic fracturing, etc.) called scenarios. The methods presented in this section increase consistency between seismic measurements and, thus, reliability of the structural representations.

(40) The joint VVAz and AVAz inversion may be performed so as to yield anisotropy's properties consistent with both velocities and amplitudes. The joint VVAz and AVAz inversion may include determining a group velocity and a phase velocity for a layer. The phase velocity may be determined by rotating a stiffness matrix with the incidence angle and the azimuth angle or other methods. Group velocity may be determined using formula (3). Different time axes may be used for different azimuthal angle stacks. A difference in a two-way time-shift may be given by formula (5).

(41) The above-discussed methods may be implemented in a computing device 1400 connected to or including seismic data storing equipment as illustrated in FIG. 14. Hardware, firmware, software or a combination thereof may be used to perform the various steps and operations described herein.

(42) Exemplary computing device 1400 suitable for performing the activities described in the exemplary embodiments may include a server 1401. Server 1401 may include a central processor (CPU) 1402 coupled to a random-access memory (RAM) 1404 and to a read-only memory (ROM) 1406. ROM 1406 may also be other types of storage media to store programs, such as programmable ROM (PROM), erasable PROM (EPROM), etc. Processor 1402 may communicate with other internal and external components through input/output (I/O) circuitry 1408 and bussing 1410 to provide control signals and the like. Processor 1402 carries out a variety of functions as are known in the art, as dictated by software and/or firmware instructions.

(43) Server 1401 may also include one or more data storage devices, including hard drives 1412, CD-ROM drives 1414 and other hardware capable of reading and/or storing information, such as DVD, etc. In one embodiment, software for carrying out the above-discussed steps may be stored and distributed on a CD-ROM or DVD 1416, a USB storage device 1418 or other form of media capable of portably storing information such as seismic data. These storage media may be inserted into, and read by, devices such as CD-ROM drive 1414, disk drive 1412, etc. Server 1401 may be coupled to a display 1420, which may be any type of known display or presentation screen, such as LCD, plasma display, cathode ray tube (CRT), etc. A user input interface 1422 is provided, including one or more user interface mechanisms such as a mouse, keyboard, microphone, touchpad, touch screen, voice-recognition system, etc.

(44) Server 1401 may be coupled to other devices, such as sources, detectors, etc. The server may be part of a larger network configuration as in a global area network (GAN) such as the internet 1428, which allows ultimate connection to various computing devices.

(45) The disclosed exemplary embodiments provide methods for designing or adjusting a survey plan that achieves inline and cross-line seismic data irregularities suitable for compressive sensing reconstruction. It should be understood that this description is not intended to limit the invention. On the contrary, the exemplary embodiments are intended to cover alternatives, modifications and equivalents, which are included in the spirit and scope of the invention as defined by the appended claims. Further, in the detailed description of the exemplary embodiments, numerous specific details are set forth in order to provide a comprehensive understanding of the claimed invention. However, one skilled in the art would understand that various embodiments may be practiced without such specific details.

(46) Although the features and elements of the present exemplary embodiments are described in the embodiments in particular combinations, each feature or element can be used alone without the other features and elements of the embodiments or in various combinations with or without other features and elements disclosed herein.

(47) This written description uses examples of the subject matter disclosed to enable any person skilled in the art to practice the same, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the subject matter is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims.