Methods and devices using effective elastic parameter values for anisotropic media

<u>BACKGROUND</u>

<u>TECHNICAL FIELD</u>

0001Embodiments of the subject matter disclosed herein generally relate to methods and devices for seismic data modeling and inversion for anisotropic media and, more particularly, to using isotropic-type processing with effective elastic parameter values based on anisotropy and elastic parameter values.

<u>DISCUSSION OF THE BACKGROUND</u>

0002Seismic data is routinely acquired to obtain information about subsurface structure, stratigraphy, lithology and fluids contained in underground formations (sometimes called simply "earth"). Seismic data is a term used for recordings of an underground structure's response to seismic waves (also called "acoustic waves"). The seismic waves are reflected in underground formations at locations of changes in earth properties that impact the seismic wave propagation.

0003Modeling is used to gain understanding of seismic wave propagation and to help analyze seismic data. In modeling, seismic responses are estimated using models of earth properties and seismic wave propagation modeling algorithms. Models of earth properties are often specified in terms of physical parameters. For example, modeling methods are widely used to study changes in seismic reflection amplitudes with the changing angle of incidence of a plane wave reflecting from a flat interface as described in the article by<nplcit id="ncit0001" npl-type="s"><text> Castagna, J. P., and Backus, M. M., "Offset-dependent reflectivity - Theory and practice of AVO analysis," published in Investigations in Geophysics Series No. 8, 1993</text></nplcit>, Society of Exploration Geophysicists (incorporated herewith by reference in its entirety). In this article, two half-spaces above and below an interface are assumed to be homogeneous and isotropic so that each half-space can be described with just three earth parameters, e.g., p-wave velocity, s-wave velocity and density. In practice, alternative triplets of parameters may be used, e.g., p-wave impedance, s-wave impedance and density. These parameters are referred to as elastic parameters. In some cases, modeling methods start from other earth parameters, and the transforms to elastic parameters are included as part of the modeling.

0004In seismic modeling (often referred to as "forward modeling"), a source seismic wavelet may be used to model the earth response. In an opposite process to forward modeling called "inverse modeling" or "inversion," earth parameters are estimated based on the measured seismic responses, and an estimate of the source seismic wavelet. Inversion methods are generally based on forward models of seismic wave propagation, and some make use of some input elastic parameter data such as low frequency trend information or statistical distributions. Other inversion methods use some calibration of seismic amplitudes, performed in a pre-processing step or as part of an algorithm. Dependent on the seismic data acquisition geometries, estimates of earth properties obtained from the inversion methods are generally provided as a series of 2D sections or 3D volumes of elastic parameter values. Inversion is often followed by analysis and interpretation of the inversion results. Available borehole log measurements may be used during the analysis and interpretation.

0005Local isotropy is a core assumption of most conventional methods for forward modeling, wavelet estimation, inversion, analysis and interpretation of inversion results, as well as analysis and interpretation of seismic data. Local isotropy assumption considers the earth (i.e., the underground formation) to be a stack of layers, with each layer being described by a physical parameter or property with the same value when measured in different directions. Such methods are referred to hereinafter as "isotropic" methods.

0006In reality, the layers are generally anisotropic, with the acquired seismic data embedding the effects of anisotropy. Seismic data affected by the anisotropy is called "anisotropic seismic data." Earth parameters that describe anisotropy are referred to as anisotropy parameters. Anisotropy has to be taken into consideration in order to improve the accuracy of seismic modeling, wavelet estimation, inversion, and the analysis and interpretation of inversion results for anisotropic seismic data. Ruger's monography, "<nplcit id="ncit0002" npl-type="b"><text>Reflection Coefficients and Azimuthal AVO Analysis in Anisotropic Media," Geophysical Monograph Series No. 10, Society of Exploration Geophysicists</text></nplcit>, and <nplcit id="ncit0003" npl-type="b"><text>Thomsen's article, "Understanding Seismic Anisotropy in Exploration and Exploitation," published in Distinguished Instructor Series No. 5, Society of Exploration Geophysicists/European Association of Geoscientists and Engineers</text></nplcit> (incorporated herewith by reference in their entirety) describe some approaches for handling anisotropy.

0007Taking into consideration anisotropy makes seismic processing mathematically and numerically more complex than similar isotropic methods. Modeling methods (such as wavelet estimation and certain types of seismic data analysis and interpretation) must also be expanded to incorporate the anisotropy parameters, which makes them more challenging to use. Moreover, from the inversion perspective, explicit incorporation of anisotropy parameters is difficult. Performing inversion for the elastic parameters from amplitude-variation with offset (AVO) seismic data is recognized to be a difficult problem for most seismic data acquisition geometries. If the anisotropy parameters are also included in inversion, the number of unknown values increases, making it even harder to find a solution. Adding more parameters and coping with such difficulties also complicates the analysis and interpretation of inversion results.

0008<patcit id="pcit0001" dnum="US6901333B"><text>U.S. Patent No. 6,901,333</text></patcit> and <patcit id="pcit0002" dnum="WO2015014762A"><text>WO 2015/014762</text></patcit> (which are incorporated herein by reference) describe techniques for generating anisotropic elastic parameters (here called effective elastic parameters for anisotropic medium) which are used in isotropic modeling and isotropic inversion algorithms in particular anisotropic media: vertically transverse isotropic (VTI) media, which possess polar anisotropy characteristics, and horizontally transverse isotropic (HTI) media, which possess azimuthal anisotropy characteristics. However, underground formations may not be characterized by HTI or VTI anisotropy, being instead tilted transverse isotropic (TTI) and/or orthorhombic anisotropy (which are more general anisotropy types than VTI and HTI).

0009There is a need to develop methods and devices that incorporate and generate anisotropic (effective) elastic parameters in seismic modeling, wavelet estimation, inversion and the like, as well as analysis and interpretations of results for seismic data acquired over azimuthally anisotropic media such as TTI and orthorhombic.

<u>SUMMARY</u>

0010Various embodiments for modeling and seismic data inversion related to underground formation including anisotropic media use effective elastic parameter values based on anisotropy and elastic parameter values. Take into consideration anisotropy effects in this manner allows an efficient and accurate processing for orthorhombic and tilted transverse isotropic (TTI) media. Not taking into consideration the anisotropy effects may cause a significant error in locating hydrocarbon reservoirs.

0011According to an embodiment, there is a method for seismic exploration of an underground formation including an orthorhombic anisotropic medium. The method includes obtaining seismic data acquired for the underground formation, wherein the seismic data includes seismic amplitudes detected at plural values of an azimuth angle, ω-ϕ, of a source-receiver direction relative to a principal axis of the orthorhombic anisotropic medium. The method further includes calculating elastic parameter values and anisotropy parameter values related to the orthorhombic anisotropic medium based on the seismic data, and performing an isotropic-type processing technique using effective elastic parameter values calculated based on the elastic parameter values, the anisotropy parameter values, and the azimuth angle for the orthorhombic anisotropic medium. The isotropic-type processing technique leads to a quantitative model of the underground formation usable to locate hydrocarbon reservoirs.

0012According to another embodiment, there is a method for seismic exploration of an underground formation including a tilted transverse isotropic, TTI, medium having a symmetry axis tilted at a tilt angle <i>τ</i> relative to a horizontal plane thereof. The method includes obtaining seismic data acquired for the underground formation, wherein the seismic data includes seismic amplitudes detected at plural values of an azimuth angle, ω-ϕ, of a source-receiver direction relative to a projection of the symmetry axis in a horizontal plane. The method further includes calculating elastic parameter values and anisotropy parameter values related to the TTI medium based on the seismic data, and performing an isotropic-type processing technique using effective elastic parameter values calculated based on the elastic parameter values, the anisotropy parameter values, the values of the azimuth angle and of the tilt angle for the TTI medium. The isotropic-type processing technique leads to a quantitative model of the underground formation usable to locate hydrocarbon reservoirs.

0013According to yet another embodiment, there is a seismic data processing apparatus for exploration of an underground formation. The apparatus includes an interface configured to obtain seismic data acquired for the underground formation, wherein the seismic data includes seismic amplitudes detected at plural values of an azimuth angle, ω-ϕ, of a source-receiver direction relative to (1) a principal axis of the orthorhombic anisotropic medium if the underground formation includes an orthorhombic anisotropic medium, or (2) relative to a projection of a symmetry axis in a horizontal plane, if the underground formation includes a tilted transverse isotropic, TTI, medium having a symmetry axis tilted at a tilt angle <i>τ</i> relative to an isotropy plane thereof. The apparatus also includes a central processing unit connected to the interface and configured (A) to calculate elastic parameter values and anisotropy parameter values related to the orthorhombic anisotropic medium or to the TTI medium based on the seismic data, and (B) to perform an isotropic-type processing technique using effective elastic parameter values calculated based on the elastic parameter values, and the anisotropy parameter values and the azimuth angle. The isotropic-type processing technique leads to a quantitative model of the underground formation usable to locate hydrocarbon reservoirs, and the effective elastic parameter values depend also on the tilt angles for TTI medium.

<u>BRIEF DESCRIPTION OF THE DRAWINGS</u>

0014For a more complete understanding of the present inventive concept, reference is now made to the following descriptions taken in conjunction with the accompanying drawings, in which: <ul id="ul0001" list-style="none" compact="compact"><li><figref idref="f0001">Figure 1</figref> is an illustration of a seismic data acquisition system;</li><li><figref idref="f0002">Figure 2</figref> schematically illustrates an HTI model with the asymmetry being caused by a system of vertical parallel cracks;</li><li><figref idref="f0003">Figure 3</figref> schematically illustrates an orthorhombic model with the asymmetry being caused by two orthogonal systems of parallel cracks and thin horizontal layering;</li><li><figref idref="f0004">Figures 4A and 4B</figref> illustrate a VTI and a TTI model, respectively;</li><li><figref idref="f0005">Figure 5</figref> illustrates an orthorhombic medium with orthogonal symmetry planes, having a horizontal symmetry plane (x1, x2) and the other two planes (x2,x3) and (x1,x3) being vertical;</li><li><figref idref="f0006">Figure 6</figref> is a flowchart of a method for seismic exploration of an underground formation including an orthorhombic anisotropic medium according to an embodiment;</li><li><figref idref="f0007">Figure 7</figref> is a flowchart of a method for seismic exploration of an underground formation including an orthorhombic anisotropic medium according to an embodiment;</li><li><figref idref="f0008">Figures 8</figref> and <figref idref="f0009">9</figref> are graphs illustrating results of applying methods according to embodiments; and</li><li><figref idref="f0010">Figure 10</figref> is a block diagram of a seismic data processing apparatus according to an embodiment.</li></ul>

<u>DETAILED DESCRIPTION OF THE INVENTION</u>

0015The 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 using the terminology of seismic data processing for exploring underground structures, which maximizes the information associated with parameter variations and extracted from seismic data acquired with a given set of source-receiver pairs.

0016Reference 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.

0017Methods and devices that process seismic data taking into account the effects of anisotropy associated with tilted transverse isotropic (TTI) and orthorhombic media (i.e., subsurface layers being imaged which possess azimuthal anisotropy) according to various embodiments are describes in this section. In some embodiments, earth elastic and anisotropy parameters associated with TTI and orthorhombic media are converted into anisotropic parameters used in isotropic seismic data processing such as, but not limited to, seismic modeling, wavelet estimation, inversion, the analysis and interpretation of inversion results, and the analysis and interpretation of seismic data.

0018Prior to discussing such embodiments in detail, consider system 100 for land seismic data acquisition illustrated in <figref idref="f0001">Figure 1</figref>. This system is configured to generate, transmit and receive seismic waves thereby exploring underground formations. Such exploration may determine the absence or presence of hydrocarbon deposits (such as 10 in sediment layer 12) or at least the probability of the absence or presence of hydrocarbon deposits underground.

0019System 100 includes a seismic wave source 110, which may be a vibrator moved to different locations by the first vehicle/truck 112. The seismic source generates seismic waves. Some of the seismic waves' energy is detected by receivers 130 (e.g., geophones). The receivers convert the detected waves into electrical signals recorded by seismic data acquisition system 120 (which can be located, for example, on vehicle/truck 122). Source 110, receivers 130, and data acquisition system 120 are placed on ground surface 11 and are interconnected by one or more cables 132.

0020In <figref idref="f0001">Figure 1</figref>, there is a single seismic source with a single vibrator, but more sources and vibrators may be used. Vehicle 122 may communicate wirelessly with vehicle 112 via antenna 124 and 114, respectively. Antenna 134 may facilitate wireless communications of the receivers with vehicles 112 and 122.

0021Seismic waves generated by source 110 propagate on the ground surface 11 as surface wave 111 and underground as transmitted waves 116. The transmitted waves are partially reflected at interface 13 between layers 12 and 14 (which have different wave propagation properties) as reflected waves 118. The surface and the reflected waves may be detected by receivers 130 that convert the detected seismic energy into an electrical signal. The electrical signal represents both reflected waves and surface waves, which are undesirable and filtered out as much as is practically possible during seismic data processing. Data acquisition system 120 samples the electrical signals, yielding a series of sample amplitude versus time, which is the seismic data.

0022<figref idref="f0001">Figure 1</figref> illustrates a linear land data acquisition. However, receivers may be placed on a survey area enabling data acquisition for different azimuth angles (i.e., different angles of source-receiver direction relative to a reference direction) as illustrated, for example, in <patcit id="pcit0003" dnum="US9791578B"><text>U.S. Patent No. 9,791,578</text></patcit>.

0023In marine environment, seismic data acquired with streamer spreads towed yields at least a narrow range of azimuth angles. This limited range may be extended by using an assembly of vessels as described, for example, in <patcit id="pcit0004" dnum="US9482773B"><text>U.S. Patent 9,482,773</text></patcit>, of by placing a seismic source in a middle of a streamer spread.

0024<figref idref="f0002">Figure 2</figref> illustrates an HTI medium with vertical fractures that are substantially parallel with an isotropy plane and substantially perpendicular to symmetry axes in a symmetry plane. <figref idref="f0003">Figure 3</figref> illustrates an orthorhombic medium with two orthogonal systems of parallel fractures, thus having two vertical symmetry planes. <figref idref="f0002">Figures 2</figref> and <figref idref="f0003">3</figref> are adapted from <nplcit id="ncit0004" npl-type="s"><text>1997 article "Anisotropic parameters and P-wave velocity for orthorhombic media" by Tsvankin, I. published in Geophysics, vol. 62, no. 4, pp. 1292-1309</text></nplcit>.

0025<figref idref="f0004">Figure 4A</figref> illustrates HTI and <figref idref="f0004">Figure 4B</figref> illustrates TTI. Different from HTI, in TTI the symmetry axis is tilted with respect to a horizontal plane by a given angle τ.

0026<figref idref="f0005">Figure 5</figref> illustrates another example of an anisotropic medium. The illustrated orthorhombic medium's anisotropy is caused by two orthogonal systems of vertical parallel cracks and also by a horizontal fine layering. For this medium, three symmetry planes are orthogonal to one another, with one symmetry plane being horizontal and the other two being vertical. Both TTI and orthorhombic media (also known as media with azimuthal anisotropy) are frequently associated with cracks, fractures and stress that occur in underground formations where regional differential horizontal stress is the dominant stratification factor.

0027In seismic data, important information about earth elastic and anisotropy parameters is embedded in the change of seismic amplitudes as a function of the separation between sources and receivers (a feature known as Amplitude Versus Offset AVO). AVO may be converted, for example, to angles for Amplitude Versus Angle (AVA) analysis and interpretation. Rather than studying the data at the level of records, partial stacks of records are sometimes used to reduce the amount of data and improve robustness while retaining the structural information embedded in the amplitude changes with offset, angle or other parameters.

0028The link between the subsurface parameters and AVO is determined by seismic wave propagation modeling. For the most general case, this leads to a very complex relationship that can only be solved by numerical wave equation modeling. Most conventional modeling methods are based on a simplified wave propagation model. For example, the following assumption has been used to define a simplified seismic wave propagation model: the earth (i.e., underground formation) is stratified into isotropic parallel layers; seismic waves propagate as plane waves impinging on each interface at a constant angle; each interface acts as an independent reflector; transmission effects other than ray bending are neglected; and the calculated plane wave reflection coefficients assume that half-spaces above and below the reflector interface are isotropic. Even with this very simple model of the earth and of wave propagation, the resulting Zoeppritz equations that describe the AVO relationship are quite complex, making it necessary (for practical reasons) to approximate the Zoeppritz equations. The 3-term Aki-Richards, the 3-term Shuey (which is a rearrangement of the Aki-Richards equation) and 2-term Shuey approximations have been used to approximate the Zoeppritz equations. Conventional seismic AVO inversion methods are based on these forward modeling equations or the Zoeppritz equations.

0029The following embodiments use effective elastic parameters for anisotropic medium. The effective elastic parameters, which are based on the anisotropy parameters and the elastic parameters, are used in the isotropic-type techniques applied to the seismic data. These effective elastic parameters take into consideration azimuthal anisotropy to an acceptable level of accuracy.

0030It is well-known (as described, for example, in "<nplcit id="ncit0005" npl-type="b"><text>Reflection Coefficients and Azimuthal AVO Analysis in Anisotropic Media," by Ruger, A. published in Geophysical Monograph Series No 10, Society of Exploration Geophysicists, 2002</text></nplcit>) how to incorporate anisotropy in AVO modeling (i.e., generalizing the isotropic modeling to model plane wave reflection in the case of anisotropic media). Similar to the approximation of the Zoeppritz equations by the 3-term Aki-Richards or Shuey equations, a convenient approximation to the AVO relationship can be obtained for the AVO plane wave reflection coefficients on a flat interface bounded by anisotropic half-spaces. The 3-term approximation of the Zoeppritz <i>PP</i> reflectivity is given by: <maths id="math0001" num="(1)"><math display="block"><mrow><msub><mi>R</mi><mi>p</mi></msub><mfenced><mi>θ</mi></mfenced><mo>=</mo><msub><mi>R</mi><mn>0</mn></msub><mo>+</mo><msub><mi>R</mi><mn>2</mn></msub><msup><mi>sin</mi><mn>2</mn></msup><mi>θ</mi><mo>+</mo><msub><mi>R</mi><mn>4</mn></msub><msup><mi>sin</mi><mn>2</mn></msup><mi>θ</mi><msup><mi>tan</mi><mn>2</mn></msup><mi>θ</mi></mrow></math><img file="EP3772665A1_D0001.tif" /></maths> where <i>θ</i> is the incidence angle, while the intercept <i>R₀,</i> gradient <i>R₂</i> and curvature <i>R₄</i> coefficients are defined by the elastic parameters, Thomsen parameters and the orientation of the anisotropic media with respect to the source-receiver direction.

0031The effective elastic parameters V_P^eff, V_S^eff,and ρ^eff are defined such that the PP reflectivity coefficients have an identical form to the isotropic media case of Aki-Richards (see <nplcit id="ncit0006" npl-type="b"><text>Aki, K.T., and Richards, P.G., "Quantitative Seismology: Theory and Methods," Vol. 1, published in 1980 by W.H. Freeman and Co., page 153</text></nplcit>) when expressed in terms of the anisotropic elastic parameters: <maths id="math0002" num="(2)"><math display="block"><mrow><msub><mi>R</mi><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup></mrow><mrow><msubsup><mi>V</mi><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><mi mathvariant="normal">Δ</mi><msup><mi>ρ</mi><mi mathvariant="italic">eff</mi></msup></mrow><mrow><msup><mi>ρ</mi><mi mathvariant="italic">eff</mi></msup></mrow></mfrac></mfenced></mrow></math><img file="EP3772665A1_D0002.tif" /></maths><maths id="math0003" num="(3)"><math display="block"><mrow><msub><mi>R</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup></mrow><mrow><msubsup><mi>V</mi><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup></mrow></mfrac><mo>−</mo><mn>4</mn><mi>K</mi><mfrac><mrow><mi mathvariant="normal">Δ</mi><msup><mi>G</mi><mi mathvariant="italic">eff</mi></msup></mrow><mrow><msup><mi>G</mi><mi mathvariant="italic">eff</mi></msup></mrow></mfrac></mfenced></mrow></math><img file="EP3772665A1_D0003.tif" /></maths><maths id="math0004" num="(4)"><math display="block"><mrow><msub><mi>R</mi><mn>4</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mrow><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup></mrow><mrow><msubsup><mi>V</mi><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup></mrow></mfrac></mrow></math><img file="EP3772665A1_D0004.tif" /></maths> where <maths id="math0005"><math display="inline"><mrow><msup><mi>G</mi><mi mathvariant="italic">eff</mi></msup><mo>=</mo><msup><mi>ρ</mi><mi mathvariant="italic">eff</mi></msup><msup><mrow><mfenced><msubsup><mi>V</mi><mi>S</mi><mi mathvariant="italic">eff</mi></msubsup></mfenced></mrow><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0005.tif" /></maths> is the anisotropic shear modulus.

0032An expression for the <i>PP</i> reflectivity for orthorhombic medium is given by the following expressions: <maths id="math0006" num="(5)"><math display="block"><mrow><msub><mi>R</mi><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><msub><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><mfrac><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">ρ</mi></mrow><mrow><mover><mi>ρ</mi><mrow><mo>‾</mo></mrow></mover></mrow></mfrac></mfenced></mrow></math><img file="EP3772665A1_D0006.tif" /></maths><maths id="math0007" num="(6)"><math display="block"><mrow><msub><mi>R</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><msub><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>−</mo><mn>4</mn><mi>K</mi><mfrac><mrow><mi mathvariant="normal">Δ</mi><msup><mi>G</mi><mrow><mo>∥</mo></mrow></msup></mrow><mrow><msup><mi>G</mi><mrow><mo>∥</mo></mrow></msup></mrow></mfrac><mo>+</mo><mi mathvariant="normal">Δ</mi><msup><mi mathvariant="normal">Γ</mi><mrow><mo>*</mo></mrow></msup><mo>+</mo><mi mathvariant="normal">Δ</mi><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup></mfenced></mrow></math><img file="EP3772665A1_D0007.tif" /></maths><maths id="math0008" num="(7)"><math display="block"><mrow><msub><mi>R</mi><mn>4</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><mi mathvariant="normal">Δ</mi><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup></mrow></math><img file="EP3772665A1_D0008.tif" /></maths> where V<i>_P</i> is <i>P</i>-velocity along the vertical direction, <maths id="math0009"><math display="inline"><mrow><msup><mi>G</mi><mrow><mo>∥</mo></mrow></msup><mo>=</mo><mi>ρ</mi><msup><mrow><mfenced><msubsup><mi>V</mi><mi>S</mi><mrow><mo>∥</mo></mrow></msubsup></mfenced></mrow><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0009.tif" /></maths> is the vertical shear modulus that corresponds to the shear wave with polarization parallel to (x₂,x₃) plane, <maths id="math0010"><math display="inline"><mrow><msubsup><mi>V</mi><mi>S</mi><mrow><mo>∥</mo></mrow></msubsup></mrow></math><img file="EP3772665A1_D0010.tif" /></maths> is the reference shear velocity along the vertical direction, <i>ρ</i> is density, <maths id="math0011"><math display="inline"><mrow><mi>K</mi><mo>=</mo><mfenced><mfrac><mrow><msub><mrow><mover><mi>V</mi><mrow><mo>‾</mo></mrow></mover></mrow><mi>s</mi></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac></mfenced></mrow></math><img file="EP3772665A1_D0011.tif" /></maths> with <i><o ostyle="single">V</o>_P</i> and <i><o ostyle="single">V</o>_s</i> being the above and below averages across the interface for the reference velocities. Further, for the orthorhombic case <maths id="math0012" num="(8)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mo>=</mo><msup><mi mathvariant="italic">ϵ</mi><mfenced><mn>2</mn></mfenced></msup><msup><mi>cos</mi><mn>4</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced><mo>+</mo><mfenced><msup><mi>δ</mi><mfenced><mn>3</mn></mfenced></msup><mo>+</mo><mn>2</mn><msup><mi mathvariant="italic">ϵ</mi><mfenced><mn>2</mn></mfenced></msup></mfenced><msup><mi>cos</mi><mn>2</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced><msup><mi>sin</mi><mn>2</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced><mo>+</mo><msup><mi mathvariant="italic">ϵ</mi><mfenced><mn>1</mn></mfenced></msup><msup><mi>sin</mi><mn>4</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced></mrow></math><img file="EP3772665A1_D0012.tif" /></maths><maths id="math0013" num="(9)"><math display="block"><mrow><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup><mo>=</mo><msup><mi>δ</mi><mfenced><mn>2</mn></mfenced></msup><msup><mi>cos</mi><mn>2</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced><mo>+</mo><msup><mi>δ</mi><mfenced><mn>1</mn></mfenced></msup><msup><mi>sin</mi><mn>2</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced></mrow></math><img file="EP3772665A1_D0013.tif" /></maths><maths id="math0014" num="(10)"><math display="block"><mrow><mi mathvariant="normal">Γ</mi><mo>*</mo><mo>=</mo><mn>8</mn><msup><mi mathvariant="italic">Kγ</mi><mrow><mfenced><mn>3</mn></mfenced></mrow></msup></mrow></math><img file="EP3772665A1_D0014.tif" /></maths> where ϕ is the azimuth of the x₁ principal axis (relative to an arbitrary reference), ω is the azimuth of the source - receiver direction (relative to the same reference), while the Thomsen parameters are defined in following Table 1, where <i>c_ij</i> (with <i>i,j</i> = 1,2,3,4,5,6) denotes stiffness matrix components. <tables id="tabl0001" num="0001"><table frame="all"><title>Table 1</title><tgroup cols="4"><colspec colnum="1" colname="col1" colwidth="32mm" /><colspec colnum="2" colname="col2" colwidth="34mm" /><colspec colnum="3" colname="col3" colwidth="66mm" /><colspec colnum="4" colname="col4" colwidth="35mm" /><thead><row><entry>Thomsen Parameter</entry><entry><i>ε</i></entry><entry><i>δ</i></entry><entry><i>γ</i></entry></row></thead><tbody><row><entry valign="bottom">(<i>x</i>₂<i>,x</i>₃) plane</entry><entry valign="bottom"><maths id="math0015"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mn>1</mn></mfenced></mrow></msup><mo>=</mo><mfrac><mrow><msub><mi>c</mi><mn>22</mn></msub><mo>−</mo><msub><mi>c</mi><mn>33</mn></msub></mrow><mrow><mn>2</mn><msub><mi>c</mi><mn>33</mn></msub></mrow></mfrac></mrow></math><img file="EP3772665A1_D0015.tif" /></maths></entry><entry valign="bottom"><maths id="math0016"><math display="block"><mrow><msup><mi>δ</mi><mrow><mfenced><mn>1</mn></mfenced></mrow></msup><mo>=</mo><mfrac><mrow><msup><mrow><mfenced><msub><mi>c</mi><mn>23</mn></msub><mo>+</mo><msub><mi>c</mi><mn>44</mn></msub></mfenced></mrow><mn>2</mn></msup><mo>−</mo><msup><mrow><mfenced><msub><mi>c</mi><mn>33</mn></msub><mo>−</mo><msub><mi>c</mi><mn>44</mn></msub></mfenced></mrow><mn>2</mn></msup></mrow><mrow><mn>2</mn><msub><mi>c</mi><mn>33</mn></msub><mfenced><msub><mi>c</mi><mn>33</mn></msub><mo>−</mo><msub><mi>c</mi><mn>44</mn></msub></mfenced></mrow></mfrac></mrow></math><img file="EP3772665A1_D0016.tif" /></maths></entry><entry valign="bottom" /></row><row><entry valign="bottom">(<i>x</i>₁<i>,x</i>₃) plane</entry><entry valign="bottom"><maths id="math0017"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mn>2</mn></mfenced></mrow></msup><mo>=</mo><mfrac><mrow><msub><mi>c</mi><mn>11</mn></msub><mo>−</mo><msub><mi>c</mi><mn>33</mn></msub></mrow><mrow><mn>2</mn><msub><mi>c</mi><mn>33</mn></msub></mrow></mfrac></mrow></math><img file="EP3772665A1_D0017.tif" /></maths></entry><entry valign="bottom"><maths id="math0018"><math display="block"><mrow><msup><mi>δ</mi><mrow><mfenced><mn>2</mn></mfenced></mrow></msup><mo>=</mo><mfrac><mrow><msup><mrow><mfenced><msub><mi>c</mi><mn>13</mn></msub><mo>−</mo><msub><mi>c</mi><mn>55</mn></msub></mfenced></mrow><mn>2</mn></msup><mo>−</mo><msup><mrow><mfenced><msub><mi>c</mi><mn>33</mn></msub><mo>−</mo><msub><mi>c</mi><mn>55</mn></msub></mfenced></mrow><mn>2</mn></msup></mrow><mrow><mn>2</mn><msub><mi>c</mi><mn>33</mn></msub><mfenced><msub><mi>c</mi><mn>33</mn></msub><mo>−</mo><msub><mi>c</mi><mn>55</mn></msub></mfenced></mrow></mfrac></mrow></math><img file="EP3772665A1_D0018.tif" /></maths></entry><entry valign="bottom"><maths id="math0019"><math display="block"><mrow><msup><mi>γ</mi><mrow><mfenced><mn>2</mn></mfenced></mrow></msup><mo>=</mo><mfrac><mrow><msub><mi>c</mi><mn>66</mn></msub><mo>−</mo><msub><mi>c</mi><mn>44</mn></msub></mrow><mrow><mn>2</mn><msub><mi>c</mi><mn>44</mn></msub></mrow></mfrac></mrow></math><img file="EP3772665A1_D0019.tif" /></maths></entry></row><row><entry valign="bottom">(<i>x</i>₁,<i>x</i>₂) plane</entry><entry valign="bottom" /><entry valign="bottom"><maths id="math0020"><math display="block"><mrow><msup><mi>δ</mi><mrow><mfenced><mn>3</mn></mfenced></mrow></msup><mo>=</mo><mfrac><mrow><msup><mrow><mfenced><msub><mi>c</mi><mn>12</mn></msub><mo>+</mo><msub><mi>c</mi><mn>66</mn></msub></mfenced></mrow><mn>2</mn></msup><mo>−</mo><msup><mrow><mfenced><msub><mi>c</mi><mn>11</mn></msub><mo>−</mo><msub><mi>c</mi><mn>66</mn></msub></mfenced></mrow><mn>2</mn></msup></mrow><mrow><mn>2</mn><msub><mi>c</mi><mn>11</mn></msub><mfenced><msub><mi>c</mi><mn>11</mn></msub><mo>−</mo><msub><mi>c</mi><mn>66</mn></msub></mfenced></mrow></mfrac></mrow></math><img file="EP3772665A1_D0020.tif" /></maths></entry><entry valign="bottom"><maths id="math0021"><math display="block"><mrow><msup><mi>γ</mi><mrow><mfenced><mn>2</mn></mfenced></mrow></msup><mo>=</mo><mfrac><mrow><msub><mi>c</mi><mn>44</mn></msub><mo>−</mo><msub><mi>c</mi><mn>55</mn></msub></mrow><mrow><mn>2</mn><msub><mi>c</mi><mn>55</mn></msub></mrow></mfrac></mrow></math><img file="EP3772665A1_D0021.tif" /></maths></entry></row></tbody></tgroup></table></tables>

0033The effective elastic parameters in orthorhombic media are derived in the following manner. From equation (7) for the curvature term <i>R</i>₄<maths id="math0022" num="(11)"><math display="block"><mrow><mfrac><mrow><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup></mrow><mrow><mover><mrow><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><mi mathvariant="normal">Δ</mi><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup></mrow></math><img file="EP3772665A1_D0022.tif" /></maths> and from equation (5) for intercept <i>R</i>₀<maths id="math0023" num="(12)"><math display="block"><mrow><mfrac><mrow><mi mathvariant="normal">Δ</mi><msup><mi>ρ</mi><mi mathvariant="italic">eff</mi></msup></mrow><mrow><mover><mrow><msup><mi>ρ</mi><mi mathvariant="italic">eff</mi></msup></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>=</mo><mfrac><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">ρ</mi></mrow><mrow><mover><mi>ρ</mi><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>−</mo><mi mathvariant="normal">Δ</mi><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup></mrow></math><img file="EP3772665A1_D0023.tif" /></maths> which, when substituted in equation (6), gives the gradient term <i>R</i>₂<maths id="math0024" num="(13)"><math display="block"><mrow><mfrac><mrow><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mrow><mi>S</mi><mo>∥</mo></mrow><mi mathvariant="italic">eff</mi></msubsup></mrow><mrow><mover><mrow><msubsup><mi>V</mi><mrow><mi>S</mi><mo>∥</mo></mrow><mi mathvariant="italic">eff</mi></msubsup></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mrow><mi>S</mi><mo>∥</mo></mrow></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mrow><mi>S</mi><mo>∥</mo></mrow></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>4</mn><mi>K</mi><mo>+</mo><mn>1</mn></mrow><mrow><mn>8</mn><mi>K</mi></mrow></mfrac><mi mathvariant="normal">Δ</mi><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mo>−</mo><mfrac><mn>1</mn><mrow><mn>8</mn><mi>K</mi></mrow></mfrac><msup><mi mathvariant="normal">ΔΓ</mi><mrow><mo>*</mo></mrow></msup><mo>−</mo><mfrac><mn>1</mn><mrow><mn>8</mn><mi>K</mi></mrow></mfrac><mi mathvariant="normal">Δ</mi><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup><mn>.</mn></mrow></math><img file="EP3772665A1_D0024.tif" /></maths>

0034In the small contrast limit, integrating differential equations (11), (12) and (13) yield the following formulas for the anisotropic elastic parameters: <maths id="math0025" num="(14)"><math display="block"><mrow><msubsup><mrow><mi>ln</mi><mi>V</mi></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><msub><mrow><mi>ln</mi><mi>V</mi></mrow><mi>P</mi></msub><mo>+</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup></mrow></math><img file="EP3772665A1_D0025.tif" /></maths><maths id="math0026" num="(15)"><math display="block"><mrow><mi>ln</mi><msup><mi>ρ</mi><mi mathvariant="italic">eff</mi></msup><mo>=</mo><mi>ln</mi><mi>ρ</mi><mo>−</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup></mrow></math><img file="EP3772665A1_D0026.tif" /></maths><maths id="math0027" num="(16)"><math display="block"><mrow><msubsup><mrow><mi>ln</mi><mi>V</mi></mrow><mrow><mi>S</mi><mo>∥</mo></mrow><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><msub><mrow><mi>ln</mi><mi>V</mi></mrow><mrow><mi>S</mi><mo>∥</mo></mrow></msub><mo>+</mo><mfrac><mrow><mfenced><mn>4</mn><mi>K</mi><mo>+</mo><mn>1</mn></mfenced><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mo>−</mo><msup><mi mathvariant="normal">Γ</mi><mrow><mo>*</mo></mrow></msup><mo>−</mo><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup></mrow><mrow><mn>8</mn><mi>K</mi></mrow></mfrac><mn>.</mn></mrow></math><img file="EP3772665A1_D0027.tif" /></maths>

0035Equivalent to equations (14), (15) and (16), effective elastic parameters in an orthorhombic case may use P-impedance (<i>I</i>_P) and S-impedance (<i>I</i>_S) in addition to density (<i>ρ</i>), for seismic inversion: <maths id="math0028" num="(17)"><math display="block"><mrow><msubsup><mrow><mi>ln</mi><mi>I</mi></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><msub><mrow><mi>ln</mi><mi>I</mi></mrow><mi>P</mi></msub></mrow></math><img file="EP3772665A1_D0028.tif" /></maths><maths id="math0029" num="(18)"><math display="block"><mrow><msubsup><mrow><mi>ln</mi><mi>I</mi></mrow><mrow><mi>S</mi><mo>∥</mo></mrow><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><msub><mrow><mi>ln</mi><mi>I</mi></mrow><mrow><mi>S</mi><mo>∥</mo></mrow></msub><mo>+</mo><mfrac><mrow><mfenced><mn>1</mn><mo>−</mo><mn>4</mn><mi>K</mi></mfenced><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mo>−</mo><msup><mi mathvariant="normal">Γ</mi><mrow><mo>*</mo></mrow></msup><mo>−</mo><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup></mrow><mrow><mn>8</mn><mi>K</mi></mrow></mfrac></mrow></math><img file="EP3772665A1_D0029.tif" /></maths><maths id="math0030" num="(19)"><math display="block"><mrow><mi>ln</mi><msup><mi>ρ</mi><mi mathvariant="italic">eff</mi></msup><mo>=</mo><mi>ln</mi><mi>ρ</mi><mo>−</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mn>.</mn></mrow></math><img file="EP3772665A1_D0030.tif" /></maths>

0036Table 2 lists Thomsen parameters for the particular cases of HTI and VTI. <tables id="tabl0002" num="0002"><table frame="all"><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="34mm" /><colspec colnum="2" colname="col2" colwidth="20mm" /><colspec colnum="3" colname="col3" colwidth="24mm" /><thead><row><entry align="center" valign="top">Thomsen parameter</entry><entry align="center" valign="top">HTI</entry><entry align="center" valign="top">VTI</entry></row></thead><tbody><row rowsep="0"><entry align="center"><i>ε</i>⁽¹⁾</entry><entry align="center">0</entry><entry align="center"><i>ε</i> ≈ -<i>ε</i>^(<i>V</i>)</entry></row><row><entry align="center"><i>ε</i>⁽²⁾</entry><entry align="center"><i>ε</i>^(<i>V</i>)</entry><entry align="center"><i>ε</i> ≈ -<i>ε</i>^(<i>V</i>)</entry></row><row rowsep="0"><entry align="center"><i>δ</i>⁽¹⁾</entry><entry align="center">0</entry><entry align="center"><i>δ</i> ≈ <i>δ</i>^(<i>V</i>) - 2<i>ε</i>^(<i>V</i>)</entry></row><row rowsep="0"><entry align="center"><i>δ</i>⁽²⁾</entry><entry align="center"><i>δ</i>^(<i>V</i>)</entry><entry align="center"><i>δ</i> ≈ <i>δ</i>^(<i>V</i>) - 2<i>ε</i>^(<i>V</i>)</entry></row><row><entry align="center"><i>δ</i>⁽³⁾</entry><entry align="center"><i>δ</i>^(<i>V</i>) - 2<i>ε</i>^(<i>V</i>)</entry><entry align="center">0</entry></row><row rowsep="0"><entry align="center"><i>γ</i>⁽²⁾</entry><entry align="center">γ^(<i>V</i>)</entry><entry align="center"><i>γ</i> ≈ -<i>γ</i>^(<i>V</i>)</entry></row><row><entry align="center"><i>γ</i>⁽³⁾</entry><entry align="center">-<i>γ</i>^(<i>V</i>)</entry><entry align="center">0</entry></row></tbody></tgroup></table></tables>

0037The substitution of HTI Thomsen parameters from Table 2 into the general expressions (8), (9) and (10) yields the following equations valid in the HTI case: <maths id="math0031" num="(20)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mo>=</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><msup><mi>cos</mi><mn>4</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced><mo>+</mo><msup><mi>δ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><msup><mi>cos</mi><mn>2</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced><msup><mi>sin</mi><mn>2</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced></mrow></math><img file="EP3772665A1_D0031.tif" /></maths><maths id="math0032" num="(21)"><math display="block"><mrow><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup><mo>=</mo><msup><mi>δ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><msup><mi>cos</mi><mn>2</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced></mrow></math><img file="EP3772665A1_D0032.tif" /></maths><maths id="math0033" num="(22)"><math display="block"><mrow><msup><mi mathvariant="normal">Γ</mi><mrow><mo>*</mo></mrow></msup><mo>=</mo><mo>−</mo><mn>8</mn><msup><mi mathvariant="italic">Kγ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><msup><mi>cos</mi><mn>2</mn></msup><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced><mn>.</mn></mrow></math><img file="EP3772665A1_D0033.tif" /></maths>

0038Further substituting equations (20), (21) and (22) into the expressions (14), (15) and (16) gives for the effective elastic parameters in the HTI case: <maths id="math0034" num="(23)"><math display="block"><mrow><msubsup><mi>V</mi><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><msubsup><mi>δ</mi><mi>r</mi><mrow><msup><mi>cos</mi><mn>2</mn></msup><mfenced><mi>ω</mi><mo>−</mo><mi>ϕ</mi></mfenced></mrow></msubsup><msup><mfenced><mfrac><mrow><msub><mi mathvariant="normal">ϵ</mi><mi>r</mi></msub></mrow><mrow><msub><mi>δ</mi><mi>r</mi></msub></mrow></mfrac></mfenced><mrow><msup><mi>cos</mi><mn>4</mn></msup><mfenced><mi mathvariant="italic">ω</mi><mo>−</mo><mi>ϕ</mi></mfenced></mrow></msup><msub><mi>V</mi><mi>P</mi></msub></mrow></math><img file="EP3772665A1_D0034.tif" /></maths><maths id="math0035" num="(24)"><math display="block"><mrow><msup><mi>ρ</mi><mi mathvariant="italic">eff</mi></msup><mo>=</mo><msubsup><mi>δ</mi><mi>r</mi><mrow><mo>−</mo><msup><mi>cos</mi><mn>2</mn></msup><mfenced><mi mathvariant="italic">ω</mi><mo>−</mo><mi>ϕ</mi></mfenced></mrow></msubsup><msup><mfenced><mfrac><mrow><msub><mi>ϵ</mi><mi>r</mi></msub></mrow><mrow><msub><mi>δ</mi><mi>r</mi></msub></mrow></mfrac></mfenced><mrow><mo>−</mo><msup><mi>cos</mi><mn>4</mn></msup><mfenced><mi mathvariant="italic">ω</mi><mo>−</mo><mi>ϕ</mi></mfenced></mrow></msup><mi>ρ</mi></mrow></math><img file="EP3772665A1_D0035.tif" /></maths><maths id="math0036" num="(25)"><math display="block"><mrow><msubsup><mi>V</mi><mrow><mi>S</mi><mo>∥</mo></mrow><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><msup><mrow><mfenced><msub><mi>γ</mi><mi>r</mi></msub><msup><mrow><msub><mi>δ</mi><mi>r</mi></msub></mrow><mn>0.5</mn></msup></mfenced></mrow><mrow><msup><mi>cos</mi><mn>2</mn></msup><mfenced><mi mathvariant="italic">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced></mrow></msup><msup><mrow><mfenced><msub><mi mathvariant="italic">ϵ</mi><mi>r</mi></msub><mo>/</mo><msub><mi>δ</mi><mi>r</mi></msub></mfenced></mrow><mrow><mfenced><mfenced><mn mathvariant="normal">4</mn><mi mathvariant="normal">K</mi><mo>+</mo><mn mathvariant="normal">1</mn></mfenced><mo>/</mo><mn mathvariant="normal">8</mn><mi mathvariant="normal">K</mi></mfenced><msup><mi>cos</mi><mn>4</mn></msup><mfenced><mi mathvariant="italic">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced></mrow></msup><msub><mi>V</mi><mrow><mi>S</mi><mo>∥</mo></mrow></msub></mrow></math><img file="EP3772665A1_D0036.tif" /></maths> where <i>ε_r, δ_r</i> and <i>γ_r</i> are defined by <maths id="math0037" num="(26)"><math display="block"><mrow><msub><mrow><mi>ln</mi><mi mathvariant="italic">ϵ</mi></mrow><mi>r</mi></msub><mo>=</mo><msup><mi mathvariant="italic">ϵ</mi><mfenced><mi>V</mi></mfenced></msup></mrow></math><img file="EP3772665A1_D0037.tif" /></maths><maths id="math0038" num="(27)"><math display="block"><mrow><msub><mrow><mi>ln</mi><mi mathvariant="italic">δ</mi></mrow><mi>r</mi></msub><mo>=</mo><msup><mi mathvariant="italic">δ</mi><mfenced><mi>V</mi></mfenced></msup></mrow></math><img file="EP3772665A1_D0038.tif" /></maths><maths id="math0039" num="(28)"><math display="block"><mrow><msub><mrow><mi>ln</mi><mi mathvariant="italic">γ</mi></mrow><mi>r</mi></msub><mo>=</mo><msup><mi mathvariant="italic">γ</mi><mfenced><mi>V</mi></mfenced></msup><mn>.</mn></mrow></math><img file="EP3772665A1_D0039.tif" /></maths>

0039As in the orthorhombic case, the reflectivity expressions for TTI media (see <figref idref="f0004">Figure 4</figref>) are derived from the general anisotropic expression for PP wave reflectivity. <maths id="math0040" num="(29)"><math display="block"><mrow><msub><mi>R</mi><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><msub><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">α</mi></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mi>τ</mi></mfenced><mo>+</mo><mfrac><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">ρ</mi></mrow><mrow><mover><mi>ρ</mi><mrow><mo>‾</mo></mrow></mover></mrow></mfrac></mfenced></mrow></math><img file="EP3772665A1_D0040.tif" /></maths><maths id="math0041" num="(30)"><math display="block"><mrow><msub><mi>R</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><msub><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">α</mi></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mi>τ</mi></mfenced><mo>−</mo><mn>4</mn><mi>K</mi><mfrac><mrow><mi mathvariant="normal">Δ</mi><msup><mi>G</mi><mrow><mo>∥</mo></mrow></msup></mrow><mrow><mover><mrow><msup><mi>G</mi><mrow><mo>∥</mo></mrow></msup></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>−</mo><mn>4</mn><msubsup><mi mathvariant="italic">Kg</mi><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mi>τ</mi></mfenced><mo>+</mo><msup><mi mathvariant="normal">ΔΓ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mi>τ</mi><mi>φ</mi><mi>ω</mi></mfenced><mo>+</mo><mi mathvariant="normal">Δ</mi><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mi>τ</mi><mi>φ</mi><mi>ω</mi></mfenced></mfenced></mrow></math><img file="EP3772665A1_D0041.tif" /></maths><maths id="math0042" num="(31)"><math display="block"><mrow><msub><mi>R</mi><mn>4</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><msub><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">V</mi></mrow><mi>P</mi></msub></mrow><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><msubsup><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">α</mi></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mi>τ</mi></mfenced><mo>+</mo><mi mathvariant="normal">Δ</mi><msup><mi>ϵ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mi>τ</mi><mi>φ</mi><mi>ω</mi></mfenced></mfenced></mrow></math><img file="EP3772665A1_D0042.tif" /></maths> where <i>V_P</i> is the <i>P</i>-velocity as referenced within the isotropy symmetry plane, <i>G</i>^∥ is the shear modulus <maths id="math0043"><math display="inline"><mrow><msup><mi>G</mi><mrow><mo>∥</mo></mrow></msup><mo>=</mo><mi>ρ</mi><msup><mrow><mfenced><msubsup><mi>V</mi><mi>S</mi><mrow><mo>∥</mo></mrow></msubsup></mfenced></mrow><mn>2</mn></msup><mo>,</mo></mrow></math><img file="EP3772665A1_D0043.tif" /></maths><maths id="math0044"><math display="inline"><mrow><msubsup><mi>V</mi><mi>S</mi><mrow><mo>∥</mo></mrow></msubsup></mrow></math><img file="EP3772665A1_D0044.tif" /></maths> is the reference shear velocity that corresponds to the fast shear wave as defined within the isotropy symmetry plane, ρ is density, <i>K =</i> (<i><o ostyle="single">V_S</o></i>/<i><o ostyle="single">V_P</o></i>)² with <i><o ostyle="single">V_S</o></i> and <i><o ostyle="single">V_P</o></i> being the above and below averages across the interface for the reference velocities. Further, <maths id="math0045" num="(32)"><math display="block"><mrow><msubsup><mi>α</mi><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mi>τ</mi></mfenced><mo>=</mo><msup><mi>δ</mi><mfenced><mi>V</mi></mfenced></msup><msup><mfenced><mi>sin</mi><mi>τ</mi></mfenced><mn>2</mn></msup><mo>+</mo><mfenced><msup><mi>ϵ</mi><mfenced><mi>V</mi></mfenced></msup><mo>−</mo><msup><mi>δ</mi><mfenced><mi>V</mi></mfenced></msup></mfenced><msup><mfenced><mi>sin</mi><mi>τ</mi></mfenced><mn>4</mn></msup></mrow></math><img file="EP3772665A1_D0045.tif" /></maths><maths id="math0046" num="(33)"><math display="block"><mrow><msubsup><mi>g</mi><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mi>τ</mi></mfenced><mo>=</mo><mn>2</mn><msup><mi>γ</mi><mfenced><mi>V</mi></mfenced></msup><msup><mfenced><mi>sin</mi><mi>τ</mi></mfenced><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0046.tif" /></maths> where <i>τ</i> is the tilt angle between the symmetry axis and the horizontal plane (see <figref idref="f0004">Figure 4</figref>). The functions <i>ε*</i>(<i>τ, ϕ, ω</i>), <i>δ*</i>(<i>τ, ϕ, ω</i>) and Γ^* (<i>τ</i>, <i>ϕ</i>, <i>ω</i>) are similar to the ones in formula (8), (9) and (10) but also depend on the tilt angle <i>τ</i>. Note that here the <i>ϕ</i> is the azimuth angle of the projection of the symmetry axis in a horizontal plane. <maths id="math0047" num="(34)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mfenced><mn>1</mn></mfenced></msup><mfenced><mi>τ</mi></mfenced><mo>=</mo><mo>−</mo><mfenced><msup><mi>δ</mi><mfenced><mi>V</mi></mfenced></msup><msup><mfenced><mi>cos</mi><mi>τ</mi></mfenced><mn>2</mn></msup><mo>+</mo><msup><mi mathvariant="italic">ϵ</mi><mfenced><mi>V</mi></mfenced></msup><msup><mfenced><mi>sin</mi><mi>τ</mi></mfenced><mn>2</mn></msup></mfenced><msup><mfenced><mi>sin</mi><mi>τ</mi></mfenced><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0047.tif" /></maths><maths id="math0048" num="(35)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mn>2</mn></mfenced></mrow></msup><mfenced><mi>τ</mi></mfenced><mo>=</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><mi>cos</mi><mn>2</mn><mi>τ</mi></mrow></math><img file="EP3772665A1_D0048.tif" /></maths><maths id="math0049" num="(36)"><math display="block"><mrow><msup><mi>δ</mi><mfenced><mn>1</mn></mfenced></msup><mfenced><mi>τ</mi></mfenced><mo>=</mo><mo>−</mo><mfenced open="[" close="]"><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><mo>+</mo><mfenced><msup><mi>δ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><mo>−</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup></mfenced><mi>cos</mi><mn>2</mn><mi>τ</mi></mfenced><msup><mfenced><mi>sinτ</mi></mfenced><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0049.tif" /></maths><maths id="math0050" num="(37)"><math display="block"><mrow><msup><mi>δ</mi><mrow><mfenced><mn>2</mn></mfenced></mrow></msup><mfenced><mi>τ</mi></mfenced><mo>=</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><mi>cos</mi><mn>2</mn><mi>τ</mi><mo>+</mo><mfenced><msup><mi>δ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><mo>−</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup></mfenced><mi>cos</mi><mn>4</mn><mi>τ</mi></mrow></math><img file="EP3772665A1_D0050.tif" /></maths><maths id="math0051" num="(38)"><math display="block"><mrow><msup><mi>δ</mi><mfenced><mn>3</mn></mfenced></msup><mfenced><mi>τ</mi></mfenced><mo>=</mo><mo>−</mo><mfenced open="[" close="]"><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><mo>−</mo><mfenced><msup><mi>δ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup><mo>−</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mi>V</mi></mfenced></mrow></msup></mfenced><mi>cos</mi><mn>2</mn><mi>τ</mi></mfenced><msup><mfenced><mi>cos</mi><mi>τ</mi></mfenced><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0051.tif" /></maths><maths id="math0052" num="(39)"><math display="block"><mrow><msup><mi>γ</mi><mfenced><mn>3</mn></mfenced></msup><mfenced><mi>τ</mi></mfenced><mo>=</mo><mfenced><mo>−</mo><msup><mi>γ</mi><mfenced><mi>v</mi></mfenced></msup><mo>+</mo><mfrac><mn>1</mn><mi>K</mi></mfrac><mfenced><msup><mi>δ</mi><mfenced><mi>V</mi></mfenced></msup><mo>−</mo><msup><mi mathvariant="italic">ϵ</mi><mfenced><mi>V</mi></mfenced></msup></mfenced><msup><mfenced><mi>sin</mi><mi>τ</mi></mfenced><mn>2</mn></msup></mfenced><msup><mfenced><mi>cos</mi><mi>τ</mi></mfenced><mn>2</mn></msup><mn>.</mn></mrow></math><img file="EP3772665A1_D0052.tif" /></maths>

0040Formulas for the reflectivity coefficient were written by choosing the reference velocities within the isotropy plane. In turn, these formulas could be equivalently rewritten by choosing the reference <i>P</i>- and <i>S</i>- velocities along the anisotropic symmetry axis, by defining the so-called VTI reference system. In the VTI reference system, the reference p-wave <maths id="math0053"><math display="inline"><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>˜</mo></mrow></mover></mrow></math><img file="EP3772665A1_D0053.tif" /></maths> and s-wave <i>Ṽ_S</i> velocities are defined along the symmetry axis, and tilt angle <i>τ̃</i> of the symmetry axis is defined with respect to the vertical, i.e. <i>τ̃</i> = 90° - <i>τ</i>. <maths id="math0054" num="(40)"><math display="block"><mrow><msub><mi>R</mi><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><mi mathvariant="normal">Δ</mi><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>˜</mo></mrow></mover></mrow><mrow><mover><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>˜</mo></mrow></mover></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mover><mi>α</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mo>+</mo><mfrac><mrow><mi mathvariant="normal">Δ</mi><mi mathvariant="italic">ρ</mi></mrow><mi>ρ</mi></mfrac></mfenced></mrow></math><img file="EP3772665A1_D0054.tif" /></maths><maths id="math0055" num="(41)"><math display="block"><mrow><msub><mi>R</mi><mn>2</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><mi mathvariant="normal">Δ</mi><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>˜</mo></mrow></mover></mrow><mrow><mover><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>˜</mo></mrow></mover></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mover><mi>α</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>−</mo><mn>4</mn><mi>K</mi><mfrac><mrow><mi mathvariant="normal">Δ</mi><mover><mi>G</mi><mrow><mo>˜</mo></mrow></mover></mrow><mrow><mover><mrow><mover><mi>G</mi><mrow><mo>˜</mo></mrow></mover></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>−</mo><mn>4</mn><msubsup><mi mathvariant="italic">Kg</mi><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>+</mo><msup><mi mathvariant="normal">ΔΓ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mo>+</mo><mi mathvariant="normal">Δ</mi><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced></mfenced></mrow></math><img file="EP3772665A1_D0055.tif" /></maths><maths id="math0056" num="(42)"><math display="block"><mrow><msub><mi>R</mi><mn>4</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfenced><mfrac><mrow><mi mathvariant="normal">Δ</mi><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>˜</mo></mrow></mover></mrow><mrow><mover><mrow><mover><mrow><msub><mi>V</mi><mi>P</mi></msub></mrow><mrow><mo>˜</mo></mrow></mover></mrow><mrow><mo>‾</mo></mrow></mover></mrow></mfrac><mo>+</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mover><mi>α</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mi>τ</mi></mfenced><mo>+</mo><mi mathvariant="normal">Δ</mi><msup><mi>ϵ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced></mfenced></mrow></math><img file="EP3772665A1_D0056.tif" /></maths> where <maths id="math0057"><math display="inline"><mrow><mover><mi>G</mi><mrow><mo>˜</mo></mrow></mover><mo>=</mo><mi>ρ</mi><msubsup><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>S</mi><mn>2</mn></msubsup></mrow></math><img file="EP3772665A1_D0057.tif" /></maths> is the shear modulus, and <maths id="math0058" num="(43)"><math display="block"><mrow><msubsup><mrow><mover><mi>α</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>=</mo><mi>δ</mi><msup><mrow><mfenced><mi>sin</mi><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow><mn>2</mn></msup><mo>+</mo><mfenced><mi mathvariant="italic">ϵ</mi><mo>−</mo><mi>δ</mi></mfenced><msup><mrow><mfenced><mi>sin</mi><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow><mn>4</mn></msup></mrow></math><img file="EP3772665A1_D0058.tif" /></maths><maths id="math0059" num="(44)"><math display="block"><mrow><msubsup><mrow><mover><mi>g</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>=</mo><mn>2</mn><mi>γ</mi><msup><mrow><mfenced><mi>sin</mi><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow><mn>2</mn></msup><mn>.</mn></mrow></math><img file="EP3772665A1_D0059.tif" /></maths> Here, the tilt angle <i>τ̃</i> is the angle between the symmetry axis and the vertical direction. The functions <i>ε*</i>(<i>τ̃, ϕ, ω</i>), <i>δ*</i>(<i>τ̃, ϕ, ω</i>) and Γ* (<i>τ̃</i>, <i>ϕ, ω</i>) are given by formula similar with (34)-(39): <maths id="math0060" num="(45)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mn>1</mn></mfenced></mrow></msup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>=</mo><mfenced><mi mathvariant="italic">ϵ</mi><mo>+</mo><mfenced><mi mathvariant="italic">ϵ</mi><mo>−</mo><mi>δ</mi></mfenced><msup><mrow><mfenced><mi>sin</mi><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow><mn>2</mn></msup></mfenced><msup><mrow><mfenced><mi>cos</mi><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0060.tif" /></maths><maths id="math0061" num="(46)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mn>2</mn></msup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>=</mo><mi mathvariant="italic">ϵ</mi><mi mathvariant="normal">cos</mi><mn>2</mn><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mrow></math><img file="EP3772665A1_D0061.tif" /></maths><maths id="math0062" num="(47)"><math display="block"><mrow><msup><mi>δ</mi><mrow><mfenced><mn>1</mn></mfenced></mrow></msup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>=</mo><mfenced open="[" close="]"><mi mathvariant="italic">ϵ</mi><mo>+</mo><mfenced><mi>δ</mi><mo>−</mo><mi mathvariant="italic">ϵ</mi></mfenced><mi>cos</mi><mn>2</mn><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><msup><mrow><mfenced><mi>cos</mi><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0062.tif" /></maths><maths id="math0063" num="(48)"><math display="block"><mrow><msup><mi>δ</mi><mfenced><mn>2</mn></mfenced></msup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>=</mo><mi mathvariant="italic">ϵ</mi><mi>cos</mi><mn>2</mn><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mo>+</mo><mfenced><mi>δ</mi><mo>−</mo><mi mathvariant="italic">ϵ</mi></mfenced><mi>cos</mi><mn>4</mn><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mrow></math><img file="EP3772665A1_D0063.tif" /></maths><maths id="math0064" num="(49)"><math display="block"><mrow><msup><mi>δ</mi><mrow><mfenced><mn>3</mn></mfenced></mrow></msup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>=</mo><mfenced open="[" close="]"><mi mathvariant="italic">ϵ</mi><mo>−</mo><mfenced><mi>δ</mi><mo>−</mo><mi mathvariant="italic">ϵ</mi></mfenced><mi>cos</mi><mn>2</mn><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><msup><mrow><mfenced><mi>sin</mi><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow><mn>2</mn></msup></mrow></math><img file="EP3772665A1_D0064.tif" /></maths><maths id="math0065" num="(50)"><math display="block"><mrow><msup><mi>γ</mi><mrow><mfenced><mn>3</mn></mfenced></mrow></msup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced><mo>=</mo><mfenced><mi>γ</mi><mo>+</mo><mfrac><mn>1</mn><mi>K</mi></mfrac><mfenced><mi>δ</mi><mo>−</mo><mi mathvariant="italic">ϵ</mi></mfenced></mfenced><msup><mrow><mfenced><mi>sin</mi><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow><mn>2</mn></msup><mn>.</mn></mrow></math><img file="EP3772665A1_D0065.tif" /></maths>

0041Thus, the anisotropic effective elastic parameters for TTI media with VTI reference convention are: <maths id="math0066" num="(51)"><math display="block"><mrow><mi>ln</mi><msubsup><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><mi>ln</mi><msub><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>P</mi></msub><mo>+</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mo>+</mo><mi mathvariant="normal">Δ</mi><msubsup><mrow><mover><mi>α</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow></math><img file="EP3772665A1_D0066.tif" /></maths><maths id="math0067" num="(52)"><math display="block"><mrow><mi>ln</mi><msup><mrow><mover><mi>ρ</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi mathvariant="italic">eff</mi></msup><mo>=</mo><mi>ln</mi><mi>ρ</mi><mo>−</mo><msup><mi>ϵ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced></mrow></math><img file="EP3772665A1_D0067.tif" /></maths><maths id="math0068" num="(53)"><math display="block"><mrow><mi>ln</mi><msubsup><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>S</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><mi>ln</mi><msub><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>S</mi></msub><mo>+</mo><mn>1</mn><mo>/</mo><mn>8</mn><mi>K</mi><mfenced open="[" close="]"><mfenced><mn>4</mn><mi>K</mi><mo>+</mo><mn>1</mn></mfenced><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mo>−</mo><msup><mi mathvariant="normal">Γ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mo>−</mo><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mo>+</mo><mn>4</mn><mi>K</mi><msubsup><mrow><mover><mi>g</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mfenced><mn>.</mn></mrow></math><img file="EP3772665A1_D0068.tif" /></maths>

0042Effective elastic parameters in the TTI case may use P-impedance (<i>I</i>_P) and S-impedance (<i>I</i>_S) in addition to density (<i>ρ</i>), for seismic inversion, yielding the following formulas that are equivalent to equations (51), (52) and (53): <maths id="math0069" num="(54)"><math display="block"><mrow><mi>ln</mi><msubsup><mrow><mover><mi>I</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><msub><mrow><mi>ln</mi><mi>I</mi></mrow><mi>P</mi></msub><mo>+</mo><msubsup><mrow><mover><mi>α</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mrow></math><img file="EP3772665A1_D0069.tif" /></maths><maths id="math0070" num="(55)"><math display="block"><mrow><mi>ln</mi><msubsup><mrow><mover><mi>I</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>S</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><mi>ln</mi><msub><mi>I</mi><mi>S</mi></msub><mo>+</mo><mfrac><mn>1</mn><mrow><mn>8</mn><mi>K</mi></mrow></mfrac><mfenced open="[" close="]"><mfenced><mn>1</mn><mo>−</mo><mn>4</mn><mi>K</mi></mfenced><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mo>−</mo><msup><mi mathvariant="normal">Γ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mo>−</mo><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mo>+</mo><mn>4</mn><mi>K</mi><msubsup><mrow><mover><mi>g</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mfenced><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover></mfenced></mfenced></mrow></math><img file="EP3772665A1_D0070.tif" /></maths><maths id="math0071" num="(56)"><math display="block"><mrow><mi>ln</mi><msup><mrow><mover><mi>ρ</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi mathvariant="italic">eff</mi></msup><mo>=</mo><mi>ln</mi><mi>ρ</mi><mo>−</mo><msup><mi>ϵ</mi><mrow><mo>*</mo></mrow></msup><mfenced separators=",,"><mover><mi>τ</mi><mrow><mo>˜</mo></mrow></mover><mi>φ</mi><mi>ω</mi></mfenced><mn>.</mn></mrow></math><img file="EP3772665A1_D0071.tif" /></maths>

0043In the particular case of VTI media <i>τ̃</i> = 0 and therefore <maths id="math0072" num="(57)"><math display="block"><mrow><msubsup><mrow><mover><mi>α</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mo>=</mo><mn>0</mn></mrow></math><img file="EP3772665A1_D0072.tif" /></maths><maths id="math0073" num="(58)"><math display="block"><mrow><msubsup><mrow><mover><mi>g</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>r</mi><mrow><mo>*</mo></mrow></msubsup><mo>=</mo><mn>0</mn></mrow></math><img file="EP3772665A1_D0073.tif" /></maths><maths id="math0074" num="(59)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mfenced><mn>1</mn></mfenced></msup><mo>=</mo><mi mathvariant="italic">ϵ</mi></mrow></math><img file="EP3772665A1_D0074.tif" /></maths><maths id="math0075" num="(60)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mfenced><mn>2</mn></mfenced></msup><mo>=</mo><mi mathvariant="italic">ϵ</mi></mrow></math><img file="EP3772665A1_D0075.tif" /></maths><maths id="math0076" num="(61)"><math display="block"><mrow><msup><mi>δ</mi><mfenced><mn>1</mn></mfenced></msup><mo>=</mo><mi>δ</mi></mrow></math><img file="EP3772665A1_D0076.tif" /></maths><maths id="math0077" num="(62)"><math display="block"><mrow><msup><mi>δ</mi><mfenced><mn>2</mn></mfenced></msup><mo>=</mo><mi>δ</mi></mrow></math><img file="EP3772665A1_D0077.tif" /></maths><maths id="math0078" num="(63)"><math display="block"><mrow><msup><mi>δ</mi><mrow><mfenced><mn>3</mn></mfenced></mrow></msup><mo>=</mo><mn>0</mn></mrow></math><img file="EP3772665A1_D0078.tif" /></maths><maths id="math0079" num="(64)"><math display="block"><mrow><msup><mi>γ</mi><mrow><mfenced><mn>3</mn></mfenced></mrow></msup><mo>=</mo><mn>0</mn></mrow></math><img file="EP3772665A1_D0079.tif" /></maths> which leads to <maths id="math0080" num="(65)"><math display="block"><mrow><msup><mi mathvariant="italic">ϵ</mi><mrow><mo>*</mo></mrow></msup><mo>=</mo><mi mathvariant="italic">ϵ</mi></mrow></math><img file="EP3772665A1_D0080.tif" /></maths><maths id="math0081" num="(66)"><math display="block"><mrow><msup><mi>δ</mi><mrow><mo>*</mo></mrow></msup><mo>=</mo><mi>δ</mi></mrow></math><img file="EP3772665A1_D0081.tif" /></maths><maths id="math0082" num="(67)"><math display="block"><mrow><msup><mi mathvariant="normal">Γ</mi><mrow><mo>*</mo></mrow></msup><mo>=</mo><mn>0.</mn></mrow></math><img file="EP3772665A1_D0082.tif" /></maths>

0044The final formula for the effective elastic parameters for VTI media (as previously presented in <patcit id="pcit0005" dnum="US6901333B"><text>U.S. Patent No. 6,901,333</text></patcit>) are: <maths id="math0083" num="(68)"><math display="block"><mrow><mi>ln</mi><msubsup><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>P</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><mi>ln</mi><msub><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>P</mi></msub><mo>+</mo><mi mathvariant="italic">ϵ</mi></mrow></math><img file="EP3772665A1_D0083.tif" /></maths><maths id="math0084" num="(69)"><math display="block"><mrow><mi>ln</mi><msup><mrow><mover><mi>ρ</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi mathvariant="italic">eff</mi></msup><mo>=</mo><mi mathvariant="normal">ln</mi><mi mathvariant="italic">ρ</mi><mo>−</mo><mi mathvariant="italic">ϵ</mi></mrow></math><img file="EP3772665A1_D0084.tif" /></maths><maths id="math0085" num="(70)"><math display="block"><mrow><mi>ln</mi><msubsup><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>S</mi><mi mathvariant="italic">eff</mi></msubsup><mo>=</mo><mi>ln</mi><msub><mrow><mover><mi>V</mi><mrow><mo>˜</mo></mrow></mover></mrow><mi>S</mi></msub><mo>+</mo><mn>1</mn><mo>/</mo><mn>8</mn><mi>K</mi><mfenced open="[" close="]"><mfenced><mn>4</mn><mi>K</mi><mo>+</mo><mn>1</mn></mfenced><mi mathvariant="italic">ϵ</mi><mo>−</mo><mi>δ</mi></mfenced><mn>.</mn></mrow></math><img file="EP3772665A1_D0085.tif" /></maths>

0045Returning now to the effective elastic parameters in orthorhombic, any elastic parameter E has a corresponding effective elastic parameter E^eff expressed as: <maths id="math0086" num="(71)"><math display="block"><mrow><msup><mi>E</mi><mi mathvariant="italic">eff</mi></msup><mo>=</mo><msubsup><mi mathvariant="italic">ϵ</mi><mi>r</mi><mrow><msup><mfenced><mn>2</mn></mfenced><mi>X</mi></msup></mrow></msubsup><msubsup><mi>δ</mi><mi>r</mi><mrow><msup><mfenced><mn>3</mn></mfenced><mi>Y</mi></msup></mrow></msubsup><msubsup><mi>γ</mi><mi>r</mi><mrow><msup><mfenced><mn>3</mn></mfenced><mi>Z</mi></msup></mrow></msubsup><msubsup><mi mathvariant="italic">ϵ</mi><mi>r</mi><mrow><msup><mfenced><mn>1</mn></mfenced><mi>U</mi></msup></mrow></msubsup><msubsup><mi>δ</mi><mi>r</mi><mrow><msup><mfenced><mn>1</mn></mfenced><mi>V</mi></msup></mrow></msubsup><msubsup><mi>δ</mi><mi>r</mi><mrow><msup><mfenced><mn>2</mn></mfenced><mi>W</mi></msup></mrow></msubsup><mi>E</mi></mrow></math><img file="EP3772665A1_D0086.tif" /></maths> with In <maths id="math0087"><math display="inline"><mrow><msubsup><mi mathvariant="italic">ϵ</mi><mi>r</mi><mrow><mfenced><mn>2</mn></mfenced></mrow></msubsup><mo>=</mo><msup><mi mathvariant="italic">ϵ</mi><mrow><mfenced><mn>2</mn></mfenced></mrow></msup><mo>,</mo></mrow></math><img file="EP3772665A1_D0087.tif" /></maths><maths id="math0088"><math display="inline"><mrow><msubsup><mrow><mi>ln</mi><mi>δ</mi></mrow><mi>r</mi><mfenced><mn>3</mn></mfenced></msubsup><mo>=</mo><msup><mi>δ</mi><mfenced><mn>3</mn></mfenced></msup><mo>,</mo></mrow></math><img file="EP3772665A1_D0088.tif" /></maths><maths id="math0089"><math display="inline"><mrow><msubsup><mrow><mi>ln</mi><mi>γ</mi></mrow><mi>r</mi><mfenced><mn>3</mn></mfenced></msubsup><mo>=</mo><msup><mi>γ</mi><mfenced><mn>3</mn></mfenced></msup><mo>,</mo></mrow></math><img file="EP3772665A1_D0089.tif" /></maths><maths id="math0090"><math display="inline"><mrow><msubsup><mrow><mi>ln</mi><mi mathvariant="italic">ϵ</mi></mrow><mi>r</mi><mfenced><mn>1</mn></mfenced></msubsup><mo>=</mo><msup><mi mathvariant="italic">ϵ</mi><mfenced><mn>1</mn></mfenced></msup><mo>,</mo></mrow></math><img file="EP3772665A1_D0090.tif" /></maths><maths id="math0091"><math display="inline"><mrow><msubsup><mrow><mi>ln</mi><mi>δ</mi></mrow><mi>r</mi><mfenced><mn>1</mn></mfenced></msubsup><mo>=</mo><msup><mi>δ</mi><mfenced><mn>1</mn></mfenced></msup><mo>,</mo></mrow></math><img file="EP3772665A1_D0091.tif" /></maths><maths id="math0092"><math display="inline"><mrow><msubsup><mrow><mi>ln</mi><mi>δ</mi></mrow><mi>r</mi><mfenced><mn>2</mn></mfenced></msubsup><mo>=</mo><msup><mi>δ</mi><mfenced><mn>2</mn></mfenced></msup><mn>.</mn></mrow></math><img file="EP3772665A1_D0092.tif" /></maths>

0046For the effective elastic parameters in TTI medium, any elastic parameter E has a corresponding effective elastic parameter E^eff expressed as: <maths id="math0093" num="(72)"><math display="block"><mrow><msup><mi>E</mi><mi mathvariant="italic">eff</mi></msup><mo>=</mo><msubsup><mi mathvariant="italic">ϵ</mi><mi>r</mi><mi>X</mi></msubsup><msubsup><mi>δ</mi><mi>r</mi><mi>Y</mi></msubsup><msubsup><mi>γ</mi><mi>r</mi><mi>Z</mi></msubsup><mi>E</mi></mrow></math><img file="EP3772665A1_D0093.tif" /></maths> with <maths id="math0094"><math display="inline"><mrow><msubsup><mi mathvariant="italic">ϵ</mi><mi>r</mi><mi>X</mi></msubsup><mo>,</mo></mrow></math><img file="EP3772665A1_D0094.tif" /></maths><maths id="math0095"><math display="inline"><mrow><msubsup><mi>δ</mi><mi>r</mi><mi>Y</mi></msubsup><mo>,</mo></mrow></math><img file="EP3772665A1_D0095.tif" /></maths><maths id="math0096"><math display="inline"><mrow><msubsup><mi>γ</mi><mi>r</mi><mi>Z</mi></msubsup></mrow></math><img file="EP3772665A1_D0096.tif" /></maths> defined by equations (26), (27) and (28). Exponents X, Y, Z, U, V, W from equations (71) and (72) are independent from the elastic parameters and Thomsen parameters, but they may depend on <i>K.</i> It should be noted that these exponents are different for each of the elastic parameters and for each type of anisotropy.

0047Alternatively, the above effective elastic parameters for both orthorhombic and TTI media may be expressed as: <maths id="math0097" num="(73)"><math display="block"><mrow><mi>ln</mi><msup><mi>E</mi><mi mathvariant="italic">eff</mi></msup><mo>=</mo><msub><mi>b</mi><mn>0</mn></msub><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mi>cos</mi><mfenced><mn>2</mn><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced></mfenced><mo>+</mo><msub><mi>b</mi><mn>4</mn></msub><mi>cos</mi><mfenced><mn>4</mn><mfenced><mi mathvariant="normal">ω</mi><mo>−</mo><mi mathvariant="normal">φ</mi></mfenced></mfenced></mrow></math><img file="EP3772665A1_D0097.tif" /></maths> coefficients <i>b</i>₀<i>, b</i>₂ and <i>b</i>₄ depending on both elastic and Thomsen parameters, and additionally to the tilt angle <i>τ</i> in the case of TTI media. These coefficients are different for each elastic parameter and each type of anisotropy.

0048All of the above-derived formulas for effective elastic parameters are valid for arbitrary contrasts of azimuthal orientations for both orthorhombic and TTI anisotropy cases, and also valid for arbitrary contrast of tilt orientation for a TTI anisotropy case.

0049<figref idref="f0006">Figure 6</figref> is a flowchart illustrating a method 600 for seismic exploration of an underground formation including an orthorhombic anisotropic medium according to an embodiment. At 610, method 600 includes obtaining seismic data acquired for the underground formation, wherein the seismic data includes seismic amplitudes detected at plural values of an azimuth angle, ω-ϕ, of a source-receiver direction relative to a principal axis of the orthorhombic anisotropic medium. Method 600 further includes calculating elastic parameter values and anisotropy parameter values related to the orthorhombic anisotropic medium based on the seismic data at 620. The calculated values may include values for V_P, V_S∥, ρ, <i>ε</i>⁽¹⁾, <i>ε</i>⁽²⁾, <i>δ</i>⁽¹⁾, <i>δ</i>⁽²⁾, <i>δ</i>⁽³⁾, <i>γ</i>⁽³⁾, <i>ω</i>, <i>ϕ.</i> These values may be obtained via conventional seismic data processing, where offset and azimuth-dependent time shifts are indicative of changes in the propagation velocity. Additionally, at the position of the well borehole, one can perform additional seismic experiments to specifically recover TTI and/or orthorhombic parameters. Such experiments are known as walk-away and walk-around VSP (Vertical Seismic Profiling).

0050Method 600 further includes performing an isotropic-type processing technique using effective elastic parameter values at 630. The effective elastic parameter values are calculated based on the elastic parameter values, the anisotropy parameter values, and the azimuth angle (as in equations (14), (15) and (16)). The isotropic-type processing technique leads to a quantitative model of the underground formation usable to locate hydrocarbon reservoirs. That is, the output of the isotropic-type processing technique is a quantitative model of the orthorhombic anisotropic medium within the orthorhombic anisotropic medium.

0051The isotropic-type processing technique may be an isotropic seismic modeling method, an isotropic seismic analysis and interpretation method, an isotropic seismic wavelet estimation method, an isotropic seismic inversion method, or an isotropic method for analysis and interpretation of inversion results to produce processed effective elastic parameter data.

0052An elastic parameter E may have a corresponding effective elastic parameter E^eff expressed as in equation (71). Alternatively, a elastic parameter E may have a corresponding effective elastic parameter E^eff expressed as in equation (73). The elastic parameters values may represent p-wave velocity, s-wave velocity and density. Alternatively, the elastic parameters values may represent p-wave impedance, s-wave impedance and density.

0053<figref idref="f0007">Figure 7</figref> is a flowchart illustrating a method 700 for seismic exploration of an underground formation including a TTI medium having a symmetry axis tilted at a tilt angle <i>τ</i> relative to a horizontal plane thereof according to an embodiment. Method 700 includes obtaining seismic data acquired for the underground formation at 710. Here, the seismic data includes seismic amplitudes detected at plural values of an azimuth angle, ω-ϕ, of a source-receiver direction relative to a projection of the symmetry axis in the horizontal plane.

0054Method 700 further includes calculating elastic parameter values and anisotropy parameter values related to the TTI medium based on the seismic data at 720. Step 720 is similar to step 620 of method 600. The calculated values may include values for V_P, V_S∥,ρ, <i>ε</i>^(<i>V</i>)<i>, δ</i>^(<i>V</i>)<i>, ω, ϕ, τ</i>.

0055Method 700 further includes performing an isotropic-type processing technique using effective elastic parameter values at 730. The effective elastic parameter values are calculated based on the elastic parameter values, the anisotropy parameter values, the values of the azimuth angle and of the tilt angle (equations (51), (52), (53)). As in the case of method 600, the isotropic-type processing technique leads to quantitative model of the underground formation usable to locate hydrocarbon reservoirs.

0056The isotropic-type processing technique in methods 600 and 700 may be an isotropic seismic modeling method to synthesize anisotropic seismic data. Additionally, the synthesized anisotropic seismic data may be used in an isotropic seismic analysis and interpretation method for analysis and interpretation of anisotropic seismic data. Accordingly, the azimuthal effective elastic parameter data may be substituted for the isotropic elastic parameters in any of the above-mentioned methods.

0057Another embodiment is directed to a method for approximating anisotropic seismic modeling by applying isotropic seismic modeling. The method includes an initial step of inputting earth elastic parameter data and earth anisotropy parameter data for an area of interest. Next, the earth elastic parameter data is transformed to obtain effective elastic parameter data based on the earth anisotropy parameter data. Isotropic seismic modeling is then applied to the transformed effective elastic parameter data. The resulting modeled anisotropic seismic data is an approximation of seismic data obtained by a corresponding anisotropic seismic modeling. The processed anisotropic seismic data is then output. The method may further include a step of substituting the effective elastic parameter data for isotropic elastic parameter data to synthesize the anisotropic seismic data. The synthesized anisotropic seismic data may be used in an isotropic analysis and interpretation method for analysis and interpretation of the anisotropic seismic data. The area of interest may be imaged by acquisition of borehole data, wide azimuth (WAZ) data, three-dimension (3D) earth models, or four-dimensional (4D) earth models. This step of transforming may further include applying appropriate transform functions that convert the earth elastic parameter data and earth anisotropy parameter data to the effective elastic parameter data.

0058<figref idref="f0008">Figures 8</figref> and <figref idref="f0009">9</figref> illustrate modeling results for layers with relatively strong anisotropy contrasts. <figref idref="f0008">Figure 8</figref> illustrates the effect of using the 3-term approximation for the Zoeppritz equations described above for TTI media, wherein curve 810 (thin line) shows the anisotropic response and curve 820 (thick line) shows the isotropic response. In addition, curve 830 (dotted line) shows the result of isotropic modeling with the effective elastic parameters specified in the table under the graph, where V<i>_P</i>^eff, V<i>_S</i>^eff and ρ^eff are the anisotropic elastic primary wave velocity, shear wave velocity and density, respectively, calculated from the transform expressions (51), (52) and (53). The results show that isotropic modeling with the effective elastic parameters closely approximates anisotropic modeling.

0059<figref idref="f0009">Figure 9</figref> illustrates the similar results as the ones in <figref idref="f0008">Figure 8</figref> (i.e., curve 910 is the anisotropic response, curve 920 is the isotropic response and curve 930 is the result of isotropic modeling with the effective elastic parameters), but for the orthorhombic anisotropy case. The effective elastic parameter values calculated with equations (14), (15) and (16) are shown in the table under the graph. Again, the results illustrate that isotropic modeling with the effective elastic parameters closely approximates anisotropic modeling.

0060Thus, the graphs in <figref idref="f0008">Figures 8</figref> and <figref idref="f0009">9</figref> show that the approximation of anisotropic modeling using the effective elastic parameters defined for orthorhombic and TTI media in isotropic modeling is accurate.

0061The equations (71) for orthorhombic and (72) for TTI cases show that the transform functions have parameters that control the transforms. In some cases, appropriate values for these transform parameters can be obtained analytically. However, in another embodiment, an iterative procedure can also be readily followed to obtain appropriate values for the transform parameters. This implies further important benefits such as: exact modeling methods rather than 3-term approximations can be used as reference for the case of the simple single flat interface model, or more complex modeling methods can be applied; and other transform functions with other functional forms and with other transform parameters than equations (71) and (72) can be applied to evaluate if a better approximation can be obtained.

0062<figref idref="f0010">Figure 10</figref> is a block diagram of a seismic data processing apparatus 1000 according to an embodiment. Apparatus 1000 includes a CPU 1002, a data storage device 1004 (that may include a database 1006), and an interface 1008. Interface 1008 is configured to obtain seismic data acquired for the underground formation. The CPU is configured to perform methods similar to 600 and/or 700 described above. The CPU may include a software or hardware interface to a computer readable recording medium such as the data storage device 1004 or an external device. The computer readable recording medium stores executable codes, which, when executed by the CPU, make the CPU perform methods similar with 600 and 700 described above.

0063Embodiments described in this section may perform in the context of an arbitrary stack of layers with different anisotropy types (e.g., orthorhombic, TTI, VTI and/or HTI).

0064The disclosed embodiments provide methods and systems for seismic exploration of an underground formation including an orthorhombic anisotropic medium or a TTI medium. It should be understood that this description is not intended to limit the invention. On the contrary, the 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 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.

0065Although the features and elements of the present 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.

0066This 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.