Anisotropic dispersion and attenuation due to wave induced fluid flow: Quasi static finite element modeling in poroelastic solids
|
|
- Everett Dickerson
- 5 years ago
- Views:
Transcription
1 Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi: /2009jb006644, 2010 Anisotropic dispersion and attenuation due to wave induced fluid flow: Quasi static finite element modeling in poroelastic solids F. Wenzlau, 1,2 J. B. Altmann, 1 and T. M. Müller 1 Received 25 May 2009; revised 15 February 2010; accepted 26 February 2010; published 14 July [1] Heterogeneous porous media such as hydrocarbon reservoir rocks are effectively described as anisotropic viscoelastic solids. They show characteristic velocity dispersion and attenuation of seismic waves within a broad frequency band, and an explanation for this observation is the mechanism of wave induced pore fluid flow. Various theoretical models quantify dispersion and attenuation of normal incident compressional waves in finely layered porous media. Similar models of shear wave attenuation are not known, nor do general theories exist to predict wave induced fluid flow effects in media with a more complex distribution of medium heterogeneities. By using finite element simulations of poroelastic relaxation, the total frequency dependent complex stiffness tensor can be computed for a porous medium with arbitrary internal heterogeneity. From the stiffness tensor, velocity dispersion and frequency dependent attenuation are derived for compressional and shear waves as a function of the angle of incidence. We apply our approach to the case of layered media and to that of an ellipsoidal poroelastic inclusion. In the case of the ellipsoidal inclusion, compressional and shear wave modes show significant attenuation, and the characteristic frequency dependence of the effect is governed by the spatiotemporal scale of the pore fluid pressure relaxation. In our anisotropic examples, the angle dependence of the attenuation is stronger than that of the velocity dispersion. It becomes clear that the spatial attenuation patterns show specific characteristics of waveinduced fluid flow, implying that anisotropic attenuation measurements may contribute to the inversion of fluid transport properties in heterogeneous porous media. Citation: Wenzlau, F., J. B. Altmann, and T. M. Müller (2010), Anisotropic dispersion and attenuation due to wave induced fluid flow: Quasi static finite element modeling in poroelastic solids, J. Geophys. Res., 115,, doi: /2009jb Introduction [2] In generally anisotropic, dissipative media waves are attenuated differently in different directions. There has been a recent interest in the anisotropy of seismic attenuation as it may provide additional information about subsurface elastic properties. In particular, velocity and attenuation anisotropy have important implications for fracture characterization in exploration seismology [Liu et al., 1993, 2007]. Anisotropic dispersion in anisotropic media is also relevant for upscaling ultrasonic and sonic borehole measurements [Wong et al., 2008]. The observability of attenuation anisotropy is demonstrated by laboratory measurements of ultrasonic waves [Best et al., 2007; Chichinina et al., 2008; Zhu et al., 2007] as well as by VSP field measurements [Maultzsch et al., 2007]. 1 CSIRO Earth Science and Resource Engineering, Perth, Western Australia, Australia. 2 Now at Fraunhofer ITWM, Kaiserslautern, Germany. Copyright 2010 by the American Geophysical Union /10/2009JB [3] Wave attenuation in linear, isotropic viscoelastic media is studied in the classical works of Buchen [1971] and Borcherdt and Wennerberg [1985]. Variation of plane wave attenuation with the angle of incidence in layered, anelastic media has been analyzed by various authors [e.g., Winterstein, 1987; Krebes, 1983]. Recently, these formulations have been generalized to anisotropic viscoelasticity by Červený and Pšenčík [2005, 2008], Sharma [2008], or Borcherdt [2009], revealing the directionality of compressional and shear wave attenuation. Simplified formulas of attenuation anisotropy in the form of a so called Q matrix were introduced by Zhu and Tsvankin [2006] in analogy to the Thomsen velocity anisotropy parameters [Thomsen, 1986]. [4] Whereas anisotropic viscoelasticity is useful as a descriptive or phenomenological concept of wave propagation in real rocks [Carcione, 1992, 2001], it does not allow any inference about the underlying physical attenuation mechanism, and therefore cannot be directly used for rock property inversion. In order to interpret the observed velocity or attenuation anisotropy suitable rock physics models are required that link acoustic signatures to a 1of12
2 physical process. In other words, to get further insight from these laboratory observations it is important to analyze the attenuation anisotropy for a specific attenuation mechanism. [5] A very important seismic attenuation mechanism occurs in porous media that are heterogeneous on mesoscopic scales, i.e., on length scales larger than the pore size d but still smaller than the seismic wavelength l [Pride et al., 2004]. Denoting the scale of the heterogeneity by a, the relation of scales is expressed by d a : While grain sizes and pore diameters on the microscale d are typically several micrometers, and seismic wave length l is of the order of tens of meters, a is usually assumed to be on the centimeter scale. In a medium with mesoscopic heterogeneities, the propagating wave causes local fluid pressure differences across internal interfaces. This drives local fluid flow and governs the overall relaxation behavior with a characteristic frequency that depends explicitly on the scale of the heterogeneity and on rock permeability. Often, this local flow mechanism is referred to as wave induced fluid flow [e.g., Müller and Gurevich, 2005]. Since multiscale heterogeneities are present in porous rocks, seismic attenuation due to wave induced flow affects a large frequency range and plays a major role at seismic frequencies. This is particularly important in fractured porous rocks. [6] The assumption of scale separation as introduced in equation (1) allows to apply Biot equations of poroelasticity [Biot, 1962] with varying coefficients to describe the effects of mesoscopic flow fields on the seismic velocities and attenuation. Applications of Biot s theory include waveinduced flow in fractured rocks [Berryman, 2007; Gurevich et al., 2009] and rock with patchy saturation [Toms et al., 2007]. In the seismic frequency range, wave induced flow can be modeled using the quasi static or low frequency version of Biot s theory. This is possible since inertial effects of fluid flow in the pore space can be neglected and the pressure equilibration process is governed by fluid pressure diffusion [Pride, 2005]. Since the seismic wavelength is assumed to be much larger than the mesoscopic heterogeneities the effective material behavior of the poroelastic composite can be described using viscoelasticity [Masson and Pride, 2007; Rubino et al., 2009; Wenzlau and Müller, 2009]. In the most general case the effective material behavior requires the concept of anisotropic viscoelasticity which has not yet been analyzed for attenuation and dispersion resulting from wave induced flow. [7] For simple shapes of heterogeneities such as spheres and layers or particular random arrangements of heterogeneities analytical solutions for attenuation and dispersion of compressional waves are available [White et al., 1975; Norris, 1993; Müller and Gurevich, 2005]. However, for arbitrary internal medium geometry attenuation and dispersion are unknown. One possibility to estimate dispersion and attenuation is to perform numerical wave simulations using the spectral ratio method or frequency shift method [e.g., Helle et al., 2003; Carcione and Picotti, 2006; Rubino et al., 2007]. In these simulations, wave signals are recorded along the propagation path, then dispersion and attenuation are calculated by measuring traveltime and the log amplitude decay. A different approach was suggested by Masson ð1þ and Pride [2007] who simulate quasistatic relaxation based on a finite differences solver for the dynamic Biot equations. Using their approach, the model size is considerably reduced to only one representative elementary volume (REV). Also, finite element simulations are used to model attenuation and dispersion [Rubino et al., 2007, 2009]. [8] In this work, the ideas of Masson and Pride [2007] and Rubino et al. [2009] are further developed and generalized. We propose a new and efficient strategy for estimating effective elastic moduli from time domain, quasi static relaxation experiments using a finite element solver for poroelastic rheology. By applying the Fourier transform we obtain complex, frequency dependent elastic stiffness tensor components, which are then used to derive velocity dispersion and attenuation for both, compressional and shear waves. To validate the approach we compute compressional wave attenuation and dispersion for two simple composites and compare it with known, exact solutions. For the case of a 3 D model containing an ellipsoidal inclusion we extract 5 components of the viscoelastic stiffness tensor thus demonstrating the applicability of the approach in effectively anisotropic media. We compute the angle and frequency dependence of attenuation and dispersion for all wave modes. 2. Effective Viscoelastic Rheology for Porous Media [9] At seismic frequencies, the effects of wave induced fluid flow are completely described by the quasistatic Biot theory, i.e., by two coupled equations for the stress equilibrium within the porous matrix and the Darcy flow of pore fluid. Writing spatial derivatives in the i direction as i and using Einstein repeated indices convention they are given by 0 j ij ðx k ; tþ; ð2þ _w iðx k ; tþ i px ð k ; tþ; ð3þ where t ij and p are the components of the stress tensor and the fluid pressure, respectively. In the Darcy equation (3), _w i are the components of the filtration velocity, defined as the relative fluid volume flux per unit area. The hydraulic parameters are permeability and fluid viscosity h. Equations (2) and (3) are coupled through the poroelastic stress strain relations, since both the stress and the strain depend on the displacement fields of the matrix and the fluid, u i and w i [Biot, 1962] ij ¼ 2" ij þ k u k ij þ k w k ij ; p ¼ M@ k u k M@ k w k ; with the strain " ij =(u i,j + u j,i )/2 and d ij being the Kronecker symbol. In these poroelastic constitutive relations, l u and m are the undrained Lamé parameters, a is the Biot coefficient and M is the so called pore space modulus. We provide a definition of the poroelastic parameters in Appendix C. [10] The calculation of effective elastic properties of a porous medium involves an upscaling process in which the medium properties on the macroscale are calculated from the properties on the mesoscale. In order to emphasize this ð4þ ð5þ 2of12
3 Figure 1. Elementary volume of a heterogeneous porous material containing mesoscale heterogeneities. change of scales between the fields on the mesoscale and on the macroscale, we use capital letters for the upscaled quantities. In the case of stress and strain, both scales are related through a spatial average over one REV as Z T ij ðþ¼ t 1 V REV ij ðx k ; tþdv; ð6þ E ij ðþ¼ t 1 Z " ij ðx k ; tþdv: ð7þ V REV with the REV sample volume V = R dv. [11] On the macroscale, the stress strain relation of a linear viscoelastic solid is expressed by T ij ðþ¼ t _R ijkl ðþ*e t kl ðþ; t such that the stress field T ij is calculated from a convolution of the strain E kl and the so called tensor of relaxation rate _R ijkl [Carcione, 2001]. Adapting the notation of Aki and Richards [1980], this quantity is conveniently written as a sum of an instantaneous part and a viscoelastic part as _R ijkl ðþ¼ t : Cijkl u ðþþ_y t ijkl ðþ t ; ð9þ u where C ijkl is the unrelaxed stiffness tensor, d(t) the delta distribution and _Y ijkl the viscoelastic part of the relaxation rate tensor. Applying the Fourier transform to equation (8), one obtains the stress strain relation in the frequency domain as ~T ij ð! Þ ¼ ~C ijkl ð! Þ~E kl ð! Þ; ð8þ ð10þ with the circular frequency w, and C ~ ijkl denoting the effective stiffness tensor. ~C ijkl is a complex, frequency dependent quantity that is related to the relaxation rate tensor by Fourier transform (denoted by F) as : ~C ijkl ð! Þ ¼ F _R ijkl ðþ t ¼ Cijkl u 1 þf _Y ijkl ðþ t : ð11þ [12] Note that in equations (9) and (11), as well as later on in equation (17), the usual summation convention does not apply, but all operations are carried out elementwise. This is indicated by the ¼ : symbol. 3. Numerical Simulation of Poroelastic Relaxation [13] In order to obtain the effective stiffness tensor of a given porous sample, we perform numerical simulations in the time domain using the finite element solver Abaqus for poroelastic rheology. First, the term relaxation experiment and the transformation from the time domain to the frequency domain are explained. Then we give the necessary boundary conditions that allow us to assess every single component of the effective stiffness tensor for the case of vertical transverse isotropy. Finally, rock physics formulas are given to compute dispersion and attenuation of compressional and shear waves from the stiffness tensor. [14] Consider a heterogeneous, fluid saturated, porous rock sample with the height L that is confined laterally, sealed hydraulically and represents one REV of a porous 0 rock, as shown in Figure 1. Initially, a global mean strain E kl is imposed instantaneously on the sample, resulting into a heterogeneous fluid pressure distribution and an initial mean stress state Tij 0 ¼ T ij ðt ¼ 0Þ ð12þ within the sample. The following pressure diffusion process governs the relaxation of the whole sample, such that while maintaining the overall strain, the mean stress decreases until the sample reaches the relaxed state T 1 ij ¼ lim T ij ðþ: t t!1 ð13þ This so called relaxation experiment allows the effective viscoelastic behavior of the heterogeneous porous medium to be characterized. [15] Inserting a Heaviside step function as strain input signal E ij ðþ¼e t ij 0 Ht ðþ; ð14þ 3of12
4 Figure 2. Schematic view of a relaxation experiment. (a) A step change in strain E is imposed in the form of a Heaviside function H(t), such that (b) the stress response is identical to the relaxation curve R(t). (c) The relaxation rate is the sum of an instantaneous part C u d(t) and a viscoelastic part C u _Y(t). (d) From the relaxation rate, the complex stiffness is obtained by Fourier transform, as indicated by the symbol F. into equation (8) and applying the theorem stiffness tensor ~C ijkl is then obtained by time differentiation and Fourier transform as indicated in equation (11). The procedure is depicted schematically in Figure 2. The input strain signal is shown in Figure 2a, at positive times it is has the constant value E 0, which is used for normalization. The stress response, divided by the input strain yields the relaxation curve R, given in Figure 2, varying between the unrelaxed limit C u and the relaxed limit C r, calculated similar to equation (18), but with the relaxed stress state T ij. Time differentiation of the relaxation curve yields the relaxation rate (Figure 2c) which is transformed into the frequency domain to calculate the complex stiffness, shown in Figure 2d. As the relaxation function, the real part of the stiffness varies between the relaxed and the unrelaxed limits, which correspond to low and high frequency material behavior, respectively. [16] The stiffness tensor ~C ijkl of an anisotropic medium without symmetry has 21 independent components which cannot be computed from a single 3 D relaxation experiment. In the most general case, it is necessary to apply six basic deformation states and measure in each experiment all components of the stress tensor in order to fully describe the material behavior [Barbero, 2008]. Here, we define a basic deformation state such the strain tensor has only one single nonzero component. [17] The situation simplifies in the case of a medium with vertical transverse isotropy (VTI) which we are analyzing in the following. A viscoelastic VTI medium has only 5 independent stiffness tensor components and these are found by using three basic relaxation experiments, two uniaxial strain tests and one shear test, as depicted in Figure 3. [18] This means that in order to obtain the effectively VTI properties of a poroelastic medium, one has to solve equations (2) and (3) within the computational domain W under the following boundary conditions: 9 u 3 ¼ E 0 L on G T ; >= u i n i ¼ 0 on G n G T ; ð19þ w i n i ¼ 0 on G; >; E 33 _fðþ*h t ðþ¼f t ðþ; t ð15þ one finds that the stress response is directly proportional to the relaxation tensor R ijkl, such that E 11 ð20þ T ij ðþ¼r t ijkl ðþe t kl 0 : ð16þ 9 u 1 ¼ E 0 L on G L ; >= u i n i ¼ 0 on G n G L ; >; w i n i ¼ 0 on G; Integrating equation (9) and inserting the result into equation (16) yields an expression for Y ijkl as with Y ijkl ðþ¼ t : T ij t ðþekl 0 1 Cijkl u Ht ðþ; ð17þ C u ijkl ¼ T 0 ij E0 kl 1 : ð18þ This means that the viscoelastic relaxation tensor Y ijkl (t) is directly calculated from the applied mean deformation state E 0 ij and the transient mean stress response T ij (t). The complex u 1 ¼ E 0 L on 9 G T ; >= u 1 ¼ 0 on G B ; u 2 ¼ u 3 ¼ w i n i ¼ 0 on G; >; E 13 ð21þ for uniaxial vertical compression, for uniaxial horizontal compression and for simple shear, respectively. Here, L is the sample edge length, G is the sample boundary, the indices T, B and L refer to the top (x 3 = L), bottom (x 3 =0) and left side (x 1 = 0) of the sample. The vector n i denotes the unit vector pointing in outward normal direction with respect to the boundary. 4of12
5 two independent components of the stiffness tensor. For example, the compressional and the shear wave relaxation functions are computed as R 3333 ðþ¼ t T 33ðÞ t ; ð29þ E 33 Figure 3. Three deformation states used to obtain all five components of the effective stiffness tensor (assuming effective transversal isotropy). (a) Vertical compression. (b) Horizontal compression. (c) Simple shear deformation. [19] In order to save computation time, two basic deformation states, uniaxial compression and simple shear, can be combined to form one deformation state that allows all but one component of the VTI stiffness tensor to be obtained by using just one numerical experiment. The result is u 3 ¼ u 1 ¼ E 0 L on G T 9 ; u 3 ¼ u 1 ¼ 0 on G B ; >= u 2 ¼ 0 on G F [ G K ; u i j G L ¼ u i j G R; w i n i ¼ 0 on G; >; E 33 þ E 13 ð22þ where the indices R, F and K refer to the right, the front and the back of the computational domain W, respectively. The five components of a VTI medium are then found using the relations R 1111 ðþe t 0 ¼ T 11 ðþ; t ð23þ R 1313 ðþ¼ t T 33ðÞ T t 11 ðþ t : ð30þ 2E 33 [21] It is well known that in anisotropic, viscoelastic solids, dispersion relations of compressional and shear waves are found by solving an eigenvalue problem of the Christoffel matrix where the eigenvalues correspond to complex velocities [e.g., Červený and Pšenčík, 2005]. In general, these velocities depend on the wave polarization, frequency, angle of incidence and on the degree of wave inhomogeneity, i.e., on the angle between the wave propagation direction and the direction of amplitude decrease. For homogeneous body waves, the directions of wave propagation and attenuation coincide, and Carcione [1992] gives simplified expressions for the complex velocities in this case. We provide the corresponding formulas in Appendix B. From the complex velocities ~V, dispersion and attenuation are derived as h i 1; V m ð!; Þ ¼ Re ~V 1 ð31þ Q 1 m m Im ~V 2 ð!; Þ ¼ m ; ð32þ Re ~V 2 m where the attenuation is given in the form of the inverse quality factor Q 1 and m ={P, SV, SH} for compressional, vertically polarized shear and horizontally polarized shear waves, respectively. R 1122 ðþe t 0 ¼ T 22 ðþ; t R 1133 ðþe t 0 ¼ T 33 ðþ; t R 1313 ðþe t 0 ¼ T 13 ðþ; t ð24þ ð25þ ð26þ 4. Results for Isotropic and Transverse Isotropic Media [22] In order to exemplify the applicability and accuracy of the proposed method, results are shown for a porous rock (consolidated sandstone (see Table 1)) with fluid pressure diffusivity D given by 2R 1212 ðþe t 0 ¼ T 11 ðþ T t 22 ðþ: t ð27þ D ¼ N : ð33þ The vertical P wave relaxation R 3333 requires an additional uniaxial compression experiment, as given by equation (19). Its value is computed as R 3333 ðþe t 0 ¼ T 33 ðþ: t ð28þ Note that in equation (22), we make use of a periodic boundary condition by tying together the displacements of the left and right edge of the model. [20] In the case of an effectively isotropic material, one uniaxial relaxation experiment is sufficient to calculate the Here, and h are the hydraulic permeability and the fluid viscosity, respectively. For a definition of the poroelastic parameter N, see Appendix C. [23] The synthetic sample is saturated with two fluids, 50% water and 50% gas. A periodically layered fluid distribution with spatial period L is assumed. For this case, an analytical solution is provided by White et al. [1975]. Results for the velocity dispersion and attenuation of a P wave are shown in Figure 4. We find that our numerically obtained estimates are in excellent agreement with the predicted values. The dispersion curve shows a smooth crossover from 5of12
6 Table 1. Material Properties Used in the Numerical Relaxation Experiments Background Soft Inclusion Gas Layer Matrix Grain bulk modulus K g (GPa) Matrix bulk modulus K d (GPa) Matrix shear modulus G d (GPa) Porosity Permeability (Darcy) Grain density r s (t/m 3 ) Fluid Bulk modulus K f (GPa) Density r f (t/m 3 ) Viscosity h (mpas) velocities around 3020 m/s at frequencies below the characteristic frequency! c ¼ 9D L 2 ð34þ to 3190 m/s at high frequencies above the characteristic frequency. Using 50 elements and 100 time steps, the standard deviation from the exact solution is less than 0.1%. In addition to this, the attenuation behavior is predicted with the same accuracy. [24] Finally, as a consequence of Hill s theorem, a partially saturated rock with isotropic matrix properties is always effectively isotropic, regardless of the distribution of the fluid phases inside the sample [Hill, 1963]. In the case of the previous example, this is easily checked by calculating the shear relaxation function according to equation (30). We confirmed that for the layered model depicted in Figure 5a, the shear modulus is indeed time independent, implying dispersion free propagation of shear waves in partially saturated rocks (results not shown). [25] The situation changes if the matrix shear modulus is inhomogeneous, which in general results in an anisotropic effective behavior of the heterogeneous medium. As pointed out earlier, in this case, multiple relaxation experiments are necessary to determine all effective elastic moduli. [26] The special case of anisotropy is now considered where the geometry of the heterogeneity has a vertical symmetry axis, as e. g. in the case horizontal layering. Then, the corresponding effective medium possesses vertical transverse isotropy (VTI) and it is described by five independent elastic constants ~C 1111, ~C 1122, ~C 1133, ~C 3333 and ~C 1313 that are found as explained previously. [27] In order to exemplify the VTI case, two heterogeneous double porosity media will now be considered, both with an effectively anisotropic, viscoelastic material behavior. The first medium is a water saturated, 20% porosity consolidated sandstone with thin, soft embedded porous unconsolidated layers with 30% porosity. Material parameters for these two rock types are given in Table 1 and a picture of an REV is shown in Figure 5b. A second model is considered with the same material parameters, but with a soft 3 D ellipsoidal inclusion instead of the porous layer, see Figure 5c. The 3 D ellipsoid has an axis of rotational symmetry in vertical direction and an aspect ratio of 1:4. Both models are chosen such that the volume fraction of the inclusion is 2.83%. [28] For the layered case, exact theoretical solutions are available for the effective P wave modulus for normal incidence ~c 3333 (w) [Norris, 1993; Brajanovski et al., 2005]. In addition to that, the poroelastic Backus averaging gives the high and low frequency limits for all three wave modes (P, SV, SH). In the case of the ellipsoidal inclusion, a theoretical solution for the effective P wave modulus has been reported recently by Galvin and Gurevich [2006]. Their model includes fast to slow mode conversion but it is restricted to the case when the wave propagation direction is parallel to the vertical symmetry axis. In addition to that, their solution assumes a fluid filled inclusion and a weak shear modulus of the surrounding rock matrix, such that G d /K f < a/b. The advantage of numerical results if compared to the previously mentioned theoretical solutions is that it is possible to obtain the total complex, frequency dependent VTI stiffness tensor. Using numerical modeling, dispersion and attenuation of all three wave modes are therefore obtained for an arbitrary angle of incidence. [29] First of all, velocity and attenuation of all wave modes in the anisotropic models are a function of the angle of incidence. An overview of the degree of velocity anisotropy for the two double porosity media as compared to the partially saturated model is given in Figure 6. The curves show slowness surfaces of shear and compressional wave modes. For comparison, we also show the slowness surfaces for the effectively isotropic partially saturated model in Figure 6a. The strongest velocity anisotropy is observed for the 1 D layered double porosity model, with local minima of P wave slowness at zero and grazing angle of incidence. The SV wave slowness has a minimum at 45 Figure 4. (a) Velocity dispersion and (b) attenuation for the 1 D partially saturated rock sample shown in Figure 5a. Frequencies are normalized using the background diffusivity D and layer period L. 6of12
7 Figure 5. Three model geometries used in the relaxation experiments. (a) Water saturated consolidated sandstone model with a gas saturated layer of thickness h. It is used in the partial saturation benchmark test example. Double porosity models are (b) a water saturated consolidated sandstone containing a weak layer or (c) an ellipsoidal inclusion with rotational symmetry about the vertical axis. Total model dimension is L L L. and that of the SV wave at 0. The slowness anisotropy of the 3 D ellipsoidal model is qualitatively similar but less pronounced. The line thickness in the polar plots reflects the amount of dispersion, which is most clearly seen in the partially saturated model that is characterized by a strong frequency dependence of the P wave. [30] In addition to the phase velocity, the attenuation characteristics of the double porosity models are also anisotropic, as shown in the polar plots of Figure 7. While the P wave attenuation of the partially saturated model is independent of the propagation direction (Figure 7a), in the layered double porosity model (Figure 7b), only P and SV waves are attenuated, and the attenuation of the pure shear wave mode is zero. P wave attenuation has a maximum at 0 and it attains zero at approximately 55. SV wave attenuation is zero at = 0 and at = 90 and maximal at 45. [31] In the 3 D case, the P and SV waves behave qualitatively like in the 1 D case. However, there is neither for the P nor for the SV wave an angle of incidence where no attenuation occurs (Figure 7c). In addition to that, attenuation is also observed for the horizontally polarized SH wave, which is, however, always smaller than the attenuation of the P and SV waves. [32] A more detailed picture of P, SV and SH wave velocities is depicted in Figure 8 for the 1 D layered model and in Figure 9 for the 3 D ellipsoidal inclusion model. The amount of velocity dispersion is indicated by the grey shaded areas between the unrelaxed or high frequency limit and the relaxed or low frequency limit. The P wave velocities show local maxima at normal (0 ) and grazing (90 ) angle of incidence, with highest values at 90. A maximum in dispersion is found at 0. The maximum SV wave velocity is found at 45 in both models. [33] The frequency dependence of P and SH wave attenuation at normal and grazing incidence are shown in Figure 10. The first observation is that maximum attenuation throughout the whole frequency band is encountered for the P wave at vertical incidence ( = 0 ) for both models, but the values are higher for the layered case by a factor of 2.5 (note the different axes scaling in Figures 10a and 10b). Second, the depicted SH wave attenuation is zero for the layered case. This should be expected, since the dispersion for horizontally polarized shear waves is zero, see Figure 8c. Interestingly, SH wave attenuation is nonzero in the 3 D case of a model with the elliptical inclusion. Actually, the SH wave attenuation of vertically incident waves is comparable to that Figure 6. Slowness surfaces obtained from numerical simulations. The blue line corresponds to P waves, the red line corresponds to SV waves, and the black line corresponds to SH waves. Results are given for (a) the 1 D gas layer model, (b) the 1 D soft layer model, and (c) the 3 D model with a soft ellipsoidal inclusion. The vertical axis corresponds to normal incident waves while the horizontal axis corresponds to grazing incidence. The line thickness indicates the amount of velocity dispersion, which is strongest in the partially saturated case (Figure 6a). 7of12
8 Figure 7. The same as Figure 6 but for the maximum wave attenuation. (a) In the case of the partially saturated sample, the P wave (blue line) shows isotropic attenuation, and the shear waves propagate without attenuation. (b and c) In the double porosity models, the SV waves (red lines) are attenuated as well and attenuation shows anisotropic behavior. In the 3 D double porosity model (Figure 7c), attenuation occurs for the SH wave as well (black line). of vertically incident P waves, attaining a value of Q 1 SH = if compared to Q 1 P = Discussion [34] Patchy saturation and its influence on seismic attenuation has been analyzed extensively in the past, e.g., by White et al. [1975] and Norris [1992]. Toms et al. [2006] give a review of different models. Two important features of wave induced flow in porous media with partial saturation is the effective isotropy of the bulk medium and the absence of shear wave attenuation, which facilitates the theoretical analysis. Both features are also observed in the layered example shown in Figure 5a. The reason for the loss free shear wave propagation lies in the fact that shear waves do not change the pore volume and therefore do not excite the fluid pressure in the porous medium. [35] There is a qualitative difference between patchy saturated media and media with a heterogeneous rock matrix. In contrast to the patchy saturation case, the effective elastic behavior of the latter is in general anisotropic and an efficient coupling between shear and compressional waves occurs. One has to keep in mind that the fluid pressure is excited only if the strain field has locally dilatational components, as seen from the divergence terms in equation (5). In homogeneous porous media, the strain associated with shear wave motion is purely deviatoric, but this is no longer the case in heterogeneous double porosity media. In other words, in heterogeneous media, shear waves may perturb the fluid pressure distribution which entails pressure relaxation and therefore acoustic attenuation. The coupling can be zero as in the case of a horizontally polarized shear wave in layered media or very strong as in the case of a vertically polarized shear wave at 45 incidence angle. In fact, in the layered model, vertically polarized shear wave attenuation exceeds attenuation of compressional waves for angles of incidence between 25 and 80. [36] Our results are in agreement with the P wave attenuation measurement reported by Zhu et al. [2007, Figure 6], with the strongest attenuation occurring in the direction parallel to the symmetry axis. It remains to be investigated how the attenuation anisotropy manifests for damping mechanisms other than wave induced flow. Given that other attenuation mechanisms produce different directionality patterns, one might potentially infer from attenuation anisotropy measurements some information about the dominating loss mechanism. [37] Further insight into the behavior of heterogeneous porous solids is found by a more detailed comparison of the 1 D and the 3 D results. One fundamental observation is Figure 8. Phase velocity of compressional and shear waves as a function of angle of incidence obtained numerically for the 1 D layered double porosity model. The plots show (a) P wave velocities, as well as those of (b) vertically polarized and (c) horizontally polarized shear waves. The grey shaded areas between the limits of unrelaxed and the relaxed fluid pressure indicate the amount of dispersion. 8of12
9 Figure 9. The same as Figure 8 but for the 3 D ellipsoidal double porosity model. While the velocities are qualitatively similar to the 1 D case, the anisotropy and the dispersion are less pronounced in the 3 D case. that in the case of the layered medium, the frequencyattenuation curves (Figure 10a) are similar such that maximum attenuation occurs at the same frequency for horizontally and vertically propagating P waves. An explanation for this phenomenon can be made if one considers the possible paths of fluid pressure equilibration. In the 1 D case, equilibration can only occur vertically between internal layers within the model, independent of what deformation has been causing the fluid pressure disequilibrium in the unrelaxed state. Consequently, there exists only one characteristic time scale of fluid pressure relaxation, scaling with the layer period. [38] The situation is different in the case of a 3 D inclusion, which can be seen from the unrelaxed fluid pressure distributions corresponding to the different applied deformation states. In Figure 11, these distributions are given for uniaxial compression and pure shear deformation such as associated with P and SH waves under normal and grazing angles of incidence. The E 12 deformation state was obtained as the difference between the two horizontal compression states as p E 12 = p E 11 p E 22. [39] While in the case of uniaxial compression, fluid pressure gradients appear mainly across the interface between the inclusion and the surrounding matrix, in the case of shear deformation, double couple distributions of fluid pressure appear and fluid flow is induced mainly between different poles within the rock matrix. In addition to that, the spatial scale of fluid pressure distribution under horizontal shear is larger than in the case of vertical shear, entailing larger relaxation times for horizontally propagating SH waves and a shift of the characteristic frequency to lower frequencies. [40] In contrast to the scheme proposed by Masson and Pride [2007], the application of a fully implicit solver allows the free choice of the time step size, independent of numerical stability restrictions. Therefore, at the beginning of the experiment, a small time step is chosen and at larger simulation times, the time step is exponentially increased. The time step size is thus adapted optimally to the shape of the relaxation curve. Another advantage of the quasistatic finite element model if compared to a fully dynamic finite difference scheme, is that it excludes inertial effects and therefore, no resonance can occur in the model. As a consequence, a true step function can be applied without using a low frequency filter as in the work by Masson and Pride [2007], thus reducing the total simulation time. [41] Rubino et al. [2009] use a frequency domain finite element Biot solver to calculate effective viscoelastic properties of randomly heterogeneous porous rocks. They use time harmonic stress fields as boundary conditions. This setup allows the calculation of the bulk and shear moduli of a heterogeneous medium directly in the frequency domain, without calculating the relaxation tensor in the first place. Figure 10. Frequency dependent attenuation Q 1 derived from double porosity relaxation experiments, (a) for the layered model and (b) for the model with an ellipsoidal inclusion. P and SH waves are considered with normal (0 ) and grazing (90 ) angle of incidence. In contrast to the layered case, where shear wave attenuation is zero, it is nonzero for the 3 D geometry. 9of12
10 Figure 11. Unrelaxed fluid pressure distributions in the 3 D finite element model for different deformation states. Shown are profile cuts through the symmetry planes of the ellipsoidal inclusion, projected to the edges of modeling domain. Fluid pressure response due to uniaxial compression in (a) vertical and (b) horizontal directions are shown as well as pure shear within the (c) vertical and (d) horizontal planes. The fluid pressure distributions correspond to P 0, P 90, SH 0,andSH 90 waves. Rubino et al. [2009] demonstrate the accuracy of their approach by a comparison with the predictions of White s model. The new aspect of our simulations if compared to the previous works is the consideration of anisotropic behavior of double porosity media. This reveals for the first time the magnitude of SH and SV wave attenuation in addition to the compressional P wave mode. [42] The application of finite element modeling for the estimation of effective elastic properties is, however, still numerically expensive if 3 D geometries are considered. For the 3 D double porosity model with an ellipsoidal inclusion, an irregular mesh with tetrahedral elements was applied. The simulation of 60 time steps required to resolve the relaxation behavior takes 14 h on 10 nodes of the SGI Altix 350 located at the Geophysical Institute of Universität Karlsruhe. The model demands 8.9GB main memory. 6. Conclusions [43] We propose a methodology to efficiently compute macroscopic viscoelastic properties of a poroelastic solid containing arbitrary mesoscopic heterogeneities. Time domain finite element relaxation experiments yield all components of the relaxation rate tensor from which the complex, frequency dependent stiffness tensor is obtained by Fourier transform. Velocity dispersion and attenuation of compressional and shear wave modes, caused by induced fluid flow, are then calculated straightforwardly. A comparison of the simulation results for a 1 D layered model with the corresponding analytical solution reveals that velocity dispersion and attenuation are modeled with very high accuracy, i.e., with a standard deviation of less than 1%. [44] Examples of a layered double porosity medium and a 3 D medium with an ellipsoidal soft poroelastic inclusion have been analyzed. Both models show effectively vertical transverse isotropy (VTI). Three deformation states are applied to estimate the total VTI stiffness tensor. In the case of a 3 D ellipsoidal inclusion, we find that the amount of SV wave attenuation is comparable with that of P waves. In contrast to the 1 D case, the SH wave is dispersive and attenuating as well, although the attenuation is smaller than that of the P or SV waves. Each wave mode induces a characteristic spatial fluid pressure distribution governing the characteristic frequency of the relaxation process. A general feature of our attenuation results for both the 1 D and the 3 D model is that maximum SV wave attenuation occurs at 45 and the local maxima of P wave attenuation occurring for normal and grazing incidence (see Figures 7b and 7c). We conclude that these patterns are characteristic for wave induced flow in VTI media. [45] The approach presented here can be easily applied to other types of internal medium heterogeneity. However, if the medium is not effectively VTI, up to six linear independent deformation states are required to estimate the full effective stiffness tensor. Appendix A: Integrals [46] In this paper, we use the Fourier transform according to the definitions F fðþ t g ð! Þ ¼ Z 1 1 F 1 ff ð! ÞgðÞ¼ t 1 2 ðþexp t ð {!tþdt; ða1þ Z 1 1 F ð! Þexp ð{!tþd!: ða2þ Since here, f is a causal function, integration (A1) is performed for t > 0 only. The convolution integral is defined as fðþ*g t ðþ¼ t Z 1 1 f ðþgt ð Þd: ða3þ Appendix B: Dispersion and Attenuation of Plane Time Harmonic Waves in VTI Media [47] Due to the vertical symmetry, a layered medium is VTI and thus characterized by five independent elastic moduli In compact notation, the stiffness tensor ~C ijkl has the form 0 1 ~C 1111 ~C 1122 ~C ~C 1122 ~C 1111 ~C ~C 1133 ~C 1133 ~C ðb1þ ~C B ~C A ~C of 12
11 where ~C 1212 is computed as (~C 1111 ~C 1122 )/2. In the isotropic case, one has furthermore ~C 1111 ¼ ~C 3333 ¼ ~P; ðb2þ ~C 1313 ¼ ~C 1212 ¼ ~G; ðb3þ they are given as [Carcione, 1992] 2~V 2 P ¼ p ~C 1111 sin 2 þ ~C 3333 cos 2 þ ~C 1313 þ ffiffiffiffi X ; ðb10þ 2~V 2 SV ¼ p ~C 1111 sin 2 þ ~C 3333 cos 2 þ ~C 1313 ffiffiffiffi X ; ðb11þ ~C 1122 ¼ ~C 1133 ¼ ~P 2~G: ðb4þ In the isotropic case, ~P and ~G are the P wave modulus and the shear wave modulus, respectively. [48] If the stiffness tensor of an effectively elastic or viscoelastic medium is known, In the general case, solving the equation of motion using a plane waves ansatz provides velocity dispersion and attenuation of three wave modes, i.e., one quasi compressional and two quasi shear modes. Following Červený and Pšenčík [2008], the calculation of the signatures involves solving the three systems for the polarization U k ðg ik ik ÞU k ¼ 0; i ¼ 1; 2; 3; ðb5þ ~V 2 SH ¼ ~C 1212 sin 2 þ ~C 1313 cos 2 : ðb12þ Please note that the structure of these equations is exactly the same as in the purely elastic case, where all elastic moduli and the velocities are real valued [Mavko et al., 1998]. Since in the viscoelastic case, the components of the stiffness tensor are complex, so are the velocities derived from equations (B10) (B12). Appendix C: Poroelastic Parameters [50] The poroelastic quantities a, M, N and the Lamé parameters l u and m are computed from the fundamental rock properties listed in Table 1 as where G ik the generalized Christoffel matrix given by G ik ¼ 1 ~ C ijkl p j p l : ðb6þ ¼ 1 K d K g ; ðc1þ The solution of this eigenvalue problem are three slowness vectors p i and polarization vectors U i corresponding to the quasi longitudinal P wave and horizontally and vertically polarized quasi shear waves SH and SV. The quantity r is averaged bulk density. From the slownesses and the polarizations, dispersion and attenuation in terms of the inverse quality factor are calculated as M ¼ þ 1 ; ðc2þ K g K f K d þ 4 3 N ¼ M G d K d þ 4 3 G d þ 2 M ; ðc3þ V ð!; ; DÞ ¼ 1 jre p i j ; ðb7þ u ¼ K d 2 3 G d þ 2 M; ðc4þ Q 1 ð!; ; DÞ ¼ 2S jim p j ; S i Re p i ðb8þ ¼ G d : ðc5þ where S i denotes the Poynting vector that can be obtained for each wave mode from the corresponding slowness and polarization. In general, the seismic signatures are frequency dependent, change with respect of the angle of incidence. In addition to that, the inhomogeneity of the wave, as controlled by the parameter D [Červený and Pšenčík, 2008], influences the dispersion and attenuation characteristics. [49] The situation is simplified when we restrict the attention to homogeneous waves, i.e., to the case when the real part and the imaginary part of the complex slowness vector are pointing in the same direction. Then, explicit expressions for V and Q are obtained for an effective viscoelastic VTI medium by first calculating three frequencydependent, complex velocities ~V P, ~V SV and ~V SH as a function of angle of incidence. Abbreviating X ¼ ~C 1111 ~C 1313 sin 2 ~C 3333 ~C 1313 cos 2 2 2sin þ ~C 1133 þ ~C ; ðb9þ [51] Acknowledgments. We thank the Deutsche Forschungsgemeinschaft and the sponsors of the PHASE consortium for supporting the research presented in this paper. References Aki, K., and P. Richards (1980), Quantitative Seismology: Theory and Methods, W. H. Freeman, San Francisco, Calif. Barbero, E. J. (2008), Finite Element Analysis of Composite Materials, CRC Press, Boca Raton, Fla. Berryman, J. G. (2007), Seismic waves in rocks with fluids and fractures, Geophys. J. Int., 171, Best, A., J. Sothcott, and C. McCann (2007), A laboratory study of seismic velocity and attenuation anisotropy in near surface sedimentary rocks, Geophys. Prospect., 55(5), Biot, M. A. (1962), Mechanics of deformation and acoustic propagation in porous media, J. Appl. Phys., 33(4), Borcherdt, R. D. (2009), Viscoelastic Waves in Layered Media, Cambridge Univ. Press, Cambridge, U. K. Borcherdt, R. D., and L. Wennerberg (1985), General P, type I S, and type II S waves in anelastic solids: Inhomogeneous wave fields in low loss solids, Bull. Seismol. Soc. Am., 75(6), Brajanovski, M., B. Gurevich, and M. Schoenberg (2005), A model for P wave attenuation and dispersion in a porous medium permeated by aligned fractures, Geophys. J. Int., 163, of 12
12 Buchen, P. (1971), Plane waves in linear viscoelastic media, Geophys. J. Int., 23(5), Carcione, J. M. (1992), Anisotropic Q and velocity dispersion of finely layered media, Geophys. Prospect., 40(7), Carcione, J. M. (2001), Wave Fields in Real Media, Pergamon, New York. Carcione, J. M., and S. Picotti (2006), P wave seismic attenuation by slowwave diffusion: Effects of inhomogeneous rock properties, Geophysics, 71, O1 O8. Červený, V., and I. Pšenčík (2005), Plane waves in viscoelastic anisotropic media I. Theory, Geophys. J. Int., 161(1), Červený, V., and I. Pšenčík (2008), Quality factor Q in dissipative anisotropic media, Geophysics, 73, T63 T75. Chichinina, T. I., I. R. Obolentseva, and G. Ronquillo Jarillo (2008), Anisotropy of seismic attenuation in fractured media: Theory and ultrasonic experiment, Transp. Porous Media, 79, 1 14, doi: / s Galvin, R. J., and B. Gurevich (2006), Interaction of an elastic wave with a circular crack in a fluid saturated porous medium, Appl. Phys. Lett., 88, Gurevich, B., M. Brajanovski, R. J. Galvin, T. M. Müller, and J. Toms Stewart (2009), P wave dispersion and attenuation in fractured and porous reservoirs Poroelasticity approach, Geophys. Prospect., 57(2), Helle, H. B., N. H. Pham, and J. M. Carcione (2003), Velocity and attenuation in partially saturated rocks: poroelastic numerical experiments, Geophys. Prospect., 51, Hill, R. (1963), Elastic properties of reinforced solids: Some theoretical principles, J. Mech. Phys. Solids, 11(5), Krebes, E. (1983), The viscoelastic reflection/transmission problem: Two special cases, Bull. Seismol. Soc. Am., 73(6), Liu, E., S. Crampin, J. Queen, and W. Rizer (1993), Velocity and attenuation anisotropy caused by microcracks and macrofractures in azimuthal reverse VSPs, Can. J. Explor. Geophys., 29, Liu, E., M. Chapman, I. Varela, X. Li, J. Queen, and H. Lynn (2007), Velocity and attenuation anisotropy: Implication of seismic fracture characterizations, Leading Edge, 26(9), Masson, Y. J., and S. R. Pride (2007), Poroelastic finite difference modeling of seismic attenuation and dispersion due to mesoscopic scale heterogeneity, J. Geophys. Res., 112, B03204, doi: /2006jb Maultzsch, S., M. Chapman, E. Liu, and X. Y. Li (2007), Modelling and analysisof attenuation anisotropy in multi azimuth VSP data from the Clair field, Geophys. Prospect., 55, Mavko, G., T. Mukerji, and J. Dvorkin (1998), The Rock Physics Handbook: Tools for Seismic Analysis of Porous Media, Cambridge Univ. Press, Cambridge, U. K. Müller, T. M., and B. Gurevich (2005), Wave induced fluid flow in random porous media: Attenuation and dispersion of elastic waves, J. Acoust. Soc. Am., 117(5), Norris, A. N. (1992), On the correspondence between poroelasticity and thermoelasticity, J. Appl. Phys., 71(3), Norris, A. N. (1993), Low frequency dispersion and attenuation in partially saturated rocks, J. Acoust. Soc. Am., 94, Pride, S. R. (2005), Relationships between seismic and hydrological properties, Hydrogeophysics, 50, , doi: / _9. Pride, S. R., J. G. Berryman, and J. M. Harris (2004), Seismic attenuation due to wave induced flow, J. Geophys. Res., 109, B01201, doi: / 2003JB Rubino, J. G., J. E. Santos, S. Picotti, and J. M. Carcione (2007), Simulation of upscaling effects due to wave induced fluid flow in Biot media using the finite element method, J. Appl. Geophys., 62(3), Rubino, J. G., C. L. Ravazzoli, and J. E. Santos (2009), Equivalent viscoelastic solids for heterogeneous fluid saturated porous rocks, Geophysics, 74, N1 N13. Sharma, M. D. (2008), Propagation of inhomogeneous plane waves in anisotropic viscoelastic media, Acta Mech., 200(3 4), Thomsen, L. (1986), Weak elastic anisotropy, Geophysics, 51(10), Toms, J., T. M. Müller, R. Ciz, and B. Gurevich (2006), Comparative review of theoretical models for elastic wave attenuation and dispersion in partially saturated rocks, Soil Dyn. Earthquake Eng., 26(6 7), Toms, J., T. M. Müller, and B. Gurevich (2007), Seismic attenuation in porous rocks with random patchy saturation, Geophys. Prospect., 55, , doi: /j x. Wenzlau, F., and T. M. Müller (2009), Finite difference modeling of wave propagation and diffusion in poroelastic media, Geophysics, 74(4), T55 T66. White, J. E., N. G. Mikhaylova, and F. M. Lyakhovitskiy (1975), Lowfrequency seismic waves in fluid saturated layered rocks, Phys. Solid Earth, 11, Winterstein, D. F. (1987), Vector attenuation: Some implications for plane waves in anelastic layered media, Geophysics, 52, Wong, R. C. K., D. R. Schmitt, D. Collis, and R. Gautam (2008), Inherent transversely isotropic elastic parameters of over consolidated shale, J. Geophys. Eng., 5, Zhu, Y., and I. Tsvankin (2006), Plane wave propagation in attenuative transversely isotropic media, Geophysics, 71, T17 T30. Zhu, Y., I. Tsvankin, P. Dewangan, and K. van Wijk (2007), Physical modeling and analysis of P wave attenuation anisotropy in transversely isotropic media, Geophysics, 72(1), D1 D7. J. B. Altmann and T. M. Müller, CSIRO Earth Science and Resource Engineering, 26 Dick Perry Ave., Kensington, Perth, WA 6151, Australia. F. Wenzlau, Fraunhofer ITWM, Fraunhofer Platz 1, D Kaiserslautern, Germany. (fabian.wenzlau@itwm.fraunhofer.de) 12 of 12
6298 Stress induced azimuthally anisotropic reservoir - AVO modeling
6298 Stress induced azimuthally anisotropic reservoir - AVO modeling M. Brajanovski* (Curtin University of Technology), B. Gurevich (Curtin University of Technology), D. Nadri (CSIRO) & M. Urosevic (Curtin
More informationGeophysical Journal International
Geophysical Journal International Geophys J Int 203 95, 679 688 Advance Access publication 203 October 3 doi: 0093/gji/ggt354 Effect of fracture fill on seismic attenuation and dispersion in fractured
More informationWave Propagation in Fractured Poroelastic Media
Wave Propagation in Fractured Poroelastic Media WCCM, MS170: Advanced Computational Techniques in Geophysical Sciences, Barcelona, Spain, July 2014 Juan E. Santos Instituto del Gas y del Petróleo (IGPUBA),
More informationEffects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses
PROC. ITB Eng. Science Vol. 38 B, No. 2, 2006, 159-170 159 Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses Fatkhan Program Studi Teknik Geofisika FIKTM-ITB
More informationWave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media
SECTION I. SEISMIC EXPLORATION Volume 38 WAVE FIELDS IN REAL MEDIA: Wave Propagation in Anisotropic, Anelastic, Porous and Electromagnetic Media (SECOND EDITION, REVISED AND EXTENDED) by Jose M. CARCIONE
More informationOrder of Authors: Juan E Santos, PhD; Robiel Martinez Corredor; Jose M Carcione, Dr.
Sciences Elsevier Editorial System(tm) for International Journal of Rock Mechanics and Mining Manuscript Draft Manuscript Number: Title: Seismic velocity and Q anisotropy in fractured poroelastic media
More informationAcoustics in Porous Media.
Acoustics in Porous Media. JUAN E. SANTOS work in collaboration with J. M. Carcione, S. Picotti, P. M. Gauzellino and R. Martinez Corredor. Department of Mathematics, Purdue University, W. Lafayette, Indiana,
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationStress-induced transverse isotropy in rocks
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 318 Stress-induced transverse isotropy in rocks Lawrence M. Schwartz, 1 William F. Murphy, III, 1 and James G. Berryman 1 ABSTRACT The application
More informationCHARACTERIZATION OF SATURATED POROUS ROCKS WITH OBLIQUELY DIPPING FRACTURES. Jiao Xue and Robert H. Tatham
CHARACTERIZATION OF SATURATED POROUS ROCS WITH OBLIQUELY DIPPING FRACTURES Jiao Xue and Robert H. Tatham Department of Geological Sciences The University of Texas at Austin ABSTRACT Elastic properties,
More informationSeismic Velocity Dispersion and the Petrophysical Properties of Porous Media
Seismic Velocity Dispersion and the Petrophysical Properties of Porous Media L. Flora Sun* University of Toronto, Toronto, ON lsun@physics.utoronto.ca and B. Milkereit University of Toronto, Toronto, ON,
More informationInterface conditions for Biot s equations of poroelasticity Boris Gurevich The Geophysical Institute of Israel, P.O. Box 2286, Holon 58122, Israel
Interface conditions for Biot s equations of poroelasticity Boris Gurevich The Geophysical Institute of Israel, P.O. Box 2286, Holon 58122, Israel Michael Schoenberg Schlumberger-Doll Research, Old Quarry
More informationWaveform inversion for attenuation estimation in anisotropic media Tong Bai & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines
Waveform inversion for attenuation estimation in anisotropic media Tong Bai & Ilya Tsvankin Center for Wave Phenomena, Colorado School of Mines SUMMARY Robust estimation of attenuation coefficients remains
More informationLawrence Berkeley National Laboratory
Lawrence Berkeley National Laboratory Peer Reviewed Title: Fracture permeability and seismic wave scattering--poroelastic linear-slip interface model for heterogeneous fractures Author: Nakagawa, S. Publication
More informationAngular and Frequency-Dependent Wave Velocity and Attenuation in Fractured Porous Media
Pure Appl. Geophys. 170 (2013), 1673 1683 Ó 2013 Springer Basel DOI 10.1007/s00024-012-0636-8 Pure and Applied Geophysics Angular and Frequency-Dependent Wave Velocity and Attenuation in Fractured Porous
More informationRESEARCH PROPOSAL. Effects of scales and extracting methods on quantifying quality factor Q. Yi Shen
RESEARCH PROPOSAL Effects of scales and extracting methods on quantifying quality factor Q Yi Shen 2:30 P.M., Wednesday, November 28th, 2012 Shen 2 Ph.D. Proposal ABSTRACT The attenuation values obtained
More informationTHE ROCK PHYSICS HANDBOOK
THE ROCK PHYSICS HANDBOOK TOOLS FOR SEISMIC ANALYSIS IN POROUS MEDIA Gary Mavko Tapan Mukerji Jack Dvorkin Stanford University Stanford University Stanford University CAMBRIDGE UNIVERSITY PRESS CONTENTS
More informationRadiation pattern in homogeneous and transversely isotropic attenuating media
Radiation pattern in homogeneous and transversely isotropic attenuating media Satish Sinha*, Sergey Abaseyev** and Evgeni Chesnokov** *Rajiv Gandhi Institute of Petroleum Technology, Rae Bareli, UP 229010
More informationPEAT SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity
PEAT8002 - SEISMOLOGY Lecture 9: Anisotropy, attenuation and anelasticity Nick Rawlinson Research School of Earth Sciences Australian National University Anisotropy Introduction Most of the theoretical
More informationUnphysical negative values of the anelastic SH plane wave energybased transmission coefficient
Shahin Moradi and Edward S. Krebes Anelastic energy-based transmission coefficient Unphysical negative values of the anelastic SH plane wave energybased transmission coefficient ABSTRACT Computing reflection
More informationFar-field radiation from seismic sources in 2D attenuative anisotropic media
CWP-535 Far-field radiation from seismic sources in 2D attenuative anisotropic media Yaping Zhu and Ilya Tsvankin Center for Wave Phenomena, Department of Geophysics, Colorado School of Mines, Golden,
More informationNumerical Modeling for Different Types of Fractures
umerical Modeling for Different Types of Fractures Xiaoqin Cui* CREWES Department of Geoscience University of Calgary Canada xicui@ucalgary.ca and Laurence R. Lines Edward S. Krebes Department of Geoscience
More informationSeismic inversion for the parameters of two orthogonal fracture sets in a VTI background medium
GEOPHYSICS, VOL. 67, NO. 1 (JANUARY-FEBRUARY 2002); P. 292 299, 3 FIGS. 10.1190/1.1451801 Seismic inversion for the parameters of two orthogonal fracture sets in a VTI background medium Andrey Bakulin,
More informationLinearized AVO in viscoelastic media Shahpoor Moradi,Kristopher A. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada
Shahpoor Moradi,Kristopher. Innanen, University of Calgary, Department of Geoscience, Calgary, Canada SUMMRY Study of linearized reflectivity is very important for amplitude versus offset VO) analysis.
More informationComparative review of theoretical models for elastic wave attenuation and dispersion in partially saturated rocks
1 Comparative review of theoretical models for elastic wave attenuation and dispersion in partially saturated rocks J. Toms a*, T. Müller b, R. Ciz c and B. Gurevich a,c a Dept. Exp. Geophys., Curtin University
More informationDownloaded 08/30/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Modeling the effect of pores and cracks interactions on the effective elastic properties of fractured porous rocks Luanxiao Zhao*, De-hua Han, Qiuliang Yao and Fuyong Yan, University of Houston; Mosab
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationP314 Anisotropic Elastic Modelling for Organic Shales
P314 Anisotropic Elastic Modelling for Organic Shales X. Wu* (British Geological Survey), M. Chapman (British Geological Survey), X.Y. Li (British Geological Survey) & H. Dai (British Geological Survey)
More informationJournal of Geophysical Research: Solid Earth
RESEARCH ARTCLE Key Points: Anisotropy of seismic attenuation correlates to anisotropy of the heterogeneity Ratios of peak attenuation in different directions depend on correlation lengths Attenuation
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationPARTICLE MOTION OF PLANE WAVES IN VISCOELASTIC ANISOTROPIC MEDIA
Geologiya and Geophysics i Geofizika Vol. 47, No. 5, pp. 557-567, 26 UDC 55.834 PARTICLE MOTION OF PLANE WAVES IN VISCOELASTIC ANISOTROPIC MEDIA V. C v erveny and I. Ps v enc v ík* Department of Geophysics,
More informationQVOA analysis: P-wave attenuation anisotropy for fracture characterization
GEOPHYSCS, VOL. 7, NO. 3 MAY-JUNE 2006 ; P. C37 C48, FGS., 2 TABLES. 0.90/.29453 QVOA analysis: P-wave attenuation anisotropy for fracture characterization Tatiana Chichinina, Vladimir Sabinin, and Gerardo
More informationP235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties
P235 Modelling Anisotropy for Improved Velocities, Synthetics and Well Ties P.W. Wild* (Ikon Science Ltd), M. Kemper (Ikon Science Ltd), L. Lu (Ikon Science Ltd) & C.D. MacBeth (Heriot Watt University)
More informationCorrespondence between the low- and high-frequency limits for anisotropic parameters in a layered medium
GEOPHYSICS VOL. 74 NO. MARCH-APRIL 009 ; P. WA5 WA33 7 FIGS. 10.1190/1.3075143 Correspondence between the low- and high-frequency limits for anisotropic parameters in a layered medium Mercia Betania Costa
More informationMcMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007
McMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007 RANDOM POLYCRYSTALS OF GRAINS WITH CRACKS: MODEL OF ELASTIC BEHAVIOR FOR FRACTURED SYSTEMS James G. Berryman Earth Sciences Division
More informationReceiver. Johana Brokešová Charles University in Prague
Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate
More informationSENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE. Juan E. Santos
SENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE Juan E. Santos Instituto del Gas y del Petróleo, Facultad de Ingeniería UBA and Department
More informationEstimation of fractured weaknesses and attenuation factors from azimuthal seismic data
Estimation of fractured weaknesses and attenuation factors from azimuthal seismic data Huaizhen Chen and Kris Innanen CREWES project Department of Geoscience University of Calgary Summary Seismic wave
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 1.510 Introduction to Seismology Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 1.510 Introduction to
More informationAnisotropic poroelasticity and wave-induced fluid flow. Harmonic finite-element simulations.
Geophys. J. Int. (2011) 000 1?? Anisotropic poroelasticity and wave-induced fluid flow. Harmonic finite-element simulations. J. M. Carcione 1 J. E. Santos 2 and S. Picotti 1 1 Istituto Nazionale di Oceanografia
More informationPhysical modeling and analysis of P-wave attenuation anisotropy in transversely isotropic media
GEOPHYSCS, VOL. 72, NO. 1 JANUARY-FEBRUARY 27 ; P. D1 D7, 8 FGS. 1.119/1.2374797 Physical modeling and analysis of P-wave attenuation anisotropy in transversely isotropic media Yaping Zhu 1, lya Tsvankin
More informationConstitutive model and wave equations for linear, viscoelastic, anisotropic media
GEOPHYSICS, VOL. 60, NO. 2 (MARCH-APRIL 995); P. 537-548, 5 FIGS., TABLE. Constitutive model and wave equations for linear, viscoelastic, anisotropic media Jose M. Carcione* ABSTRACT Rocks are far from
More informationSeismic wave attenuation and dispersion resulting from wave-induced flow in porous rocks A review
GEOPHYSICS, VOL. 75, NO. 5 SEPTEMBER-OCTOBER 2010 ; P. 75A147 75A164, 9 FIGS. 10.1190/1.3463417 Seismic wave attenuation and dispersion resulting from wave-induced flow in porous rocks A review Tobias
More informationThe title of my presentation is Constraints on C13 for shales. -2-
1-1- 1. The title of my presentation is Constraints on C13 for shales. -2-2. Outline: C13 in shales and physical models. Different Poisson's ratios in VTI rock. Compliance tensor for VTI shale. 2-3- 3.It
More informationSummary. Simple model for kerogen maturity (Carcione, 2000)
Malleswar Yenugu* and De-hua Han, University of Houston, USA Summary The conversion of kerogen to oil/gas will build up overpressure. Overpressure is caused by conversion of solid kerogen to fluid hydrocarbons
More informationThe signature of attenuation and anisotropy on AVO and inversion sensitivities
The signature of attenuation and anisotropy on AVO and inversion sensitivities Shahpoor Moradi and Kristopher Innanen University of Calgary, Department of Geoscience, Calgary, Canada Summary We study the
More informationGEOPHYSICAL PROSPECTING: DYNAMIC RESERVOIR CHARACTERIZATION AND TIME-LAPSE MULTICOMPONENT SEISMOLOGY FOR RESERVOIR MONITORING UNESCO EOLSS
GEOPHYSICAL PROSPECTING: DYNAMIC RESERVOIR CHARACTERIZATION AND TIME-LAPSE MULTICOMPONENT SEISMOLOGY FOR RESERVOIR MONITORING Steven L. Roche CGGVeritas, Multicomponent Processing & Technology Group Thomas
More informationJuan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a
Juan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a Universidad de Buenos Aires, Fac. Ing., IGPUBA, ARGENTINA b Department of Mathematics, Purdue University, USA c Universidad
More informationSENSITIVITY ANALYSIS OF AMPLITUDE VARIATION WITH OFFSET (AVO) IN FRACTURED MEDIA
SENSITIVITY ANALYSIS OF AMPLITUDE VARIATION WITH OFFSET AVO) IN FRACTURED MEDIA Mary L. Krasovec, William L. Rodi, and M. Nafi Toksoz Earth Resources Laboratory Department of Earth, Atmospheric, and Planetary
More informationD scattering of obliquely incident Rayleigh waves by a saturated alluvial valley in a layered half-space
1842. 3-D scattering of obliquely incident Rayleigh waves by a saturated alluvial valley in a layered half-space Zhenning Ba 1, Jianwen Liang 2 Department of Civil Engineering, Tianjin University, Tianjin
More informationReflection of quasi-p and quasi-sv waves at the free and rigid boundaries of a fibre-reinforced medium
Sādhan ā Vol. 7 Part 6 December 00 pp. 63 630. Printed in India Reflection of quasi-p and quasi-sv waves at the free and rigid boundaries of a fibre-reinforced medium A CHATTOPADHYAYRLKVENKATESWARLU and
More informationReflection moveout and parameter estimation for horizontal transverse isotropy
GEOPHYSICS, VOL. 62, NO. 2 (MARCH-APRIL 1997); P. 614 629, 7 FIGS. Reflection moveout and parameter estimation for horizontal transverse isotropy Ilya Tsvankin ABSTRACT Transverse isotropy with a horizontal
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. 2016) 206, 1677 1694 Advance Access publication 2016 June 27 GJI Seismology doi: 10.1093/gji/ggw245 Dynamic transverse shear modulus for a heterogeneous
More informationFrequency-Dependent Amplification of Unsaturated Surface Soil Layer
Frequency-Dependent Amplification of Unsaturated Surface Soil Layer J. Yang, M.ASCE 1 Abstract: This paper presents a study of the amplification of SV waves obliquely incident on a surface soil layer overlying
More informationABSTRACT 1. INTRODUCTION
ABSTRACT In this paper the three-dimensional transient wave propagation is investigated due to a point force applied at the interface of a fluid and a poroelastic solid. Using the total response, it is
More informationLECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # & 2 #
LECTURE 5 - Wave Equation Hrvoje Tkalčić " 2 # "t = ( $ + 2µ ) & 2 # 2 % " 2 (& ' u r ) = µ "t 2 % & 2 (& ' u r ) *** N.B. The material presented in these lectures is from the principal textbooks, other
More informationSUMMARY INTRODUCTION EXPERIMENTAL PROCEDURE
Frequency dependent elastic properties and attenuation in heavy-oil sands: comparison between measured and modeled data Agnibha Das, and Michael Batzle, Colorado School of Mines SUMMARY We have measured
More informationFractured reservoirs: An analysis of coupled elastodynamic and permeability changes from pore-pressure variation
GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 2006 ; P. O33 O41, 12 FIGS. 10.1190/1.2231108 Fractured reservoirs: An analysis of coupled elastodynamic and permeability changes from pore-pressure variation
More informationBody-wave radiation patterns and AVO in transversely isotropic media
GEOPHYSICS, VOL. 60, NO. 5 (SEPTEMBER-OCTOBER 1995); P. 1409-1425, 11 FIGS. Body-wave radiation patterns and AVO in transversely isotropic media llya Tsvankin* ABSTRACT The angular dependence of reflection
More informationProceedings, International Snow Science Workshop, Banff, 2014
SNOW SLAB FAILURE DUE TO BIOT-TYPE ACOUSTIC WAVE PROPAGATION Rolf Sidler 2 Simon Fraser University, 8888 University Drive, Burnaby V5A1S6, Canada ABSTRACT: Even though seismic methods are among the most
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationHydrogeophysics - Seismics
Hydrogeophysics - Seismics Matthias Zillmer EOST-ULP p. 1 Table of contents SH polarized shear waves: Seismic source Case study: porosity of an aquifer Seismic velocities for porous media: The Frenkel-Biot-Gassmann
More informationLINK BETWEEN ATTENUATION AND VELOCITY DISPERSION
LINK BETWEEN ATTENUATION AND VELOCITY DISPERSION Jack Dvorkin Stanford University and Rock Solid Images April 25, 2005 SUMMARY In a viscoelastic sample, the causality principle links the attenuation of
More informationEffects of VTI Anisotropy in Shale-Gas Reservoir Characterization
P-014 Summary Effects of VTI Anisotropy in Shale-Gas Reservoir Characterization Niranjan Banik* and Mark Egan, WesternGeco Shale reservoirs are one of the hottest plays in the oil industry today. Our understanding
More informationModeling elastic wave velocities and attenuation in rocks saturated with heavy oil
GEOPHYSICS, OL. 73, NO. 4 JULY-AUGUST 2008 ; P. E115 E122, 8 FIGS. 10.1190/1.2940341 Modeling elastic wave velocities and attenuation in rocks saturated with heavy oil Boris Gurevich 1,2, Konstantin Osypov
More informationP-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry
GEOPHYSICS, VOL. 62, NO. 3 (MAY-JUNE 1997); P. 713 722, 7 FIGS., 1 TABLE. P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry Andreas Rüger ABSTRACT
More informationAn acoustic wave equation for orthorhombic anisotropy
Stanford Exploration Project, Report 98, August 1, 1998, pages 6?? An acoustic wave equation for orthorhombic anisotropy Tariq Alkhalifah 1 keywords: Anisotropy, finite difference, modeling ABSTRACT Using
More informationFocal mechanisms produced by shear faulting in weakly transversely isotropic crustal rocks
GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 2006 ; P. D145 D151, 6 FIGS., 2 TABLES. 10.1190/1.2236378 Focal mechanisms produced by shear faulting in weakly transversely isotropic crustal rocks Václav
More informationSEG/San Antonio 2007 Annual Meeting. γ Q in terms of complex-valued weaknesses that relates
Chichinina* G Ronquillo-Jarillo abinin nstituto Mexicano del etróleo; Obolentseva L Gik and B Bobrov nstitute of etroleum Geology and Geophysics of iberian Branch of Russian Academy of ciences ummary We
More informationSaturation Effects of Soils on Ground Motion at Free Surface Due to Incident SV Waves
Saturation Effects of Soils on Ground Motion at Free Surface Due to Incident SV Waves Jun Yang, M.ASCE 1 Abstract: A study is presented of saturation effects of subsoil on seismic motions at the free surface
More informationMicroseismic monitoring of borehole fluid injections: Data modeling and inversion for hydraulic properties of rocks
GEOPHYSICS, VOL. 68, NO. 2 (MARCH-APRIL 2003); P. 685 689, 5 FIGS. 10.1190/1.1567239 Short Note Microseismic monitoring of borehole fluid injections: Data modeling and inversion for hydraulic properties
More informationEvaluation Of The Fracture Sizes By Using The Frequency Dependence Of Anisotropy
Evaluation Of The Fracture Sizes By Using The Frequency Dependence Of Anisotropy Vahid Mahmoudian Ataabadi Arvandan Oil and Gas Company Saeed Sajjadian Hafizollah Kashani Birgani Abstract: Frequency-dependent
More informationScienceDirect. Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including a poroelastic approach
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 63 (2014 ) 4305 4312 GHGT-12 Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including
More informationPhase-shift modelling for HTI media
HTI-Modelling Phase-shift modelling for HTI media R. K. Sharma and R. J. Ferguson ABSTRACT Fractures play an important role in hydrocarbon production as they determine the pathways and volume of crustal
More informationReflection and Transmission coefficient for VTI media
Reflection and Transmission coefficient for VTI media R. K. Sharma and R. J. Ferguson ABSTRACT VTI, R & T coefficients Presently, we obtain the reflection (R and transmission (T coefficients of plane waves
More informationReflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress
Mathematica Aeterna, Vol. 4, 4, no. 8, 877-93 Reflection of SV- Waves from the Free Surface of a Magneto-Thermoelastic Isotropic Elastic Half-Space under Initial Stress Rajneesh Kakar Faculty of Engineering
More informationTheory - SBRC fundamentals
Chapter 2 Theory - SBRC fundamentals The approach proposed to estimate hydraulic properties of rocks using microseismicity is called Seismicity Based Reservoir Characterization (SBRC). It uses a spatio-temporal
More informationScale-up in poroelastic systems and applications to reservoirs
Scale-up in poroelastic systems and applications to reservoirs James G. Berryman 1 ABSTRACT A fundamental problem of heterogeneous systems is that the macroscale behavior is not necessarily well-described
More informationWenyong Pan and Lianjie Huang. Los Alamos National Laboratory, Geophysics Group, MS D452, Los Alamos, NM 87545, USA
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 019 SGP-TR-14 Adaptive Viscoelastic-Waveform Inversion Using the Local Wavelet
More informationRelationships between the velocities and the elastic constants of an anisotropic solid possessing orthorhombic symmetry
Relationships between the velocities and the elastic constants of an anisotropic solid possessing orthorhombic symmetry R. James Brown ABSTRACT This paper reviews the equations of body-wave propagation
More informationMoveout approximation for P waves in a homogeneous VTI medium
Moveout approximation for P waves in a homogeneous VTI medium Véronique Farra 1 and Ivan Pšenčík 2 1 Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, UMR 7154 CNRS,
More informationSUMMARY NUMERICAL SIMULATION OF BRINE-CO 2 FLOW AND WAVE PROPAGATION
Sensitivity analysis of the petrophysical properties variations on the seismic response of a CO 2 storage site Juan E. Santos, Universidad de Buenos Aires, Facultad Ingeniería, Instituto del Gas y del
More informationSeismic waves in finely layered VTI media: Poroelasticity, Thomsen parameters, and fluid effects on shear waves
Stanford Exploration Project, Report 5, May 22, 24, pages 444 Seismic waves in finely layered VTI media: Poroelasticity, Thomsen parameters, and fluid effects on shear waves James G. Berryman ABSTRACT
More informationTu P8 08 Modified Anisotropic Walton Model for Consolidated Siliciclastic Rocks: Case Study of Velocity Anisotropy Modelling in a Barents Sea Well
Tu P8 08 Modified Anisotropic Walton Model for Consolidated Siliciclastic Rocks: Case Study of Velocity Anisotropy Modelling in a Barents Sea Well Y. Zhou (Rock Solid Images), F. Ruiz (Repsol), M. Ellis*
More informationComposite memory variables for viscoelastic synthetic seismograms
Geophys. J. Int. (1995) 121,634-639 RESEARCH NOTE Composite memory variables for viscoelastic synthetic seismograms Tong Xu George A. McMechan Center for Lithospheric Studies, The University of Texas at
More informationUnderstanding hydraulic fracture variability through a penny shaped crack model for pre-rupture faults
Penny shaped crack model for pre-rupture faults Understanding hydraulic fracture variability through a penny shaped crack model for pre-rupture faults David Cho, Gary F. Margrave, Shawn Maxwell and Mark
More informationEffective attenuation anisotropy of thin-layered media
GEOPHYSICS, VOL. 72, O. 5 SEPTEMBER-OCTOBER 27; P. D93 D6, 3 FIGS., TABLE..9/.275485 Effective attenuation anisotropy of thin-layered media Yaping Zhu, Ilya Tsvankin 2, Ivan Vasconcelos 3 ABSTRACT One
More informationWorkflows for Sweet Spots Identification in Shale Plays Using Seismic Inversion and Well Logs
Workflows for Sweet Spots Identification in Shale Plays Using Seismic Inversion and Well Logs Yexin Liu*, SoftMirrors Ltd., Calgary, Alberta, Canada yexinliu@softmirrors.com Summary Worldwide interest
More informationAcoustic Anisotropy Measurements and Interpretation in Deviated Wells
Acoustic Anisotropy Measurements and Interpretation in Deviated Wells X. M. Tang, and D. Patterson, Houston Technology Center, Baker Atlas, Houston, Texas, USA ABSTRACT Many acoustic anisotropy measurements
More informationCrosswell tomography imaging of the permeability structure within a sandstone oil field.
Crosswell tomography imaging of the permeability structure within a sandstone oil field. Tokuo Yamamoto (1), and Junichi Sakakibara (2) (1) University of Miami and Yamamoto Engineering Corporation, (2)
More informationEstimating rock porosity and fluid saturation using only seismic velocities
Stanford Exploration Project, Report, October 5, 999, pages 4 57 Estimating rock porosity and fluid saturation using only seismic velocities James G. Berryman, Patricia A. Berge, and Brian P. Bonner keywords:
More informationSeismic wave propagation in Kelvin visco-elastic VTI media*
APPLIED GEOPHYSICS, Vol.7, No.4 (December 21), P. 368-375, 5 Figures. DOI: 1.17/s1177-1-256-y Seismic wave propagation in Kelvin visco-elastic VTI media* Lu Jun 1 and Wang Yun 2 Abstract: In this article,
More informationA THEORETICAL MODEL FOR SITE COEFFICIENTS IN BUILDING CODE PROVISIONS
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 3029 A THEORETICAL MODEL FOR SITE COEFFICIENTS IN BUILDING CODE PROVISIONS Roger D. Borcherdt 1 SUMMARY
More informationBorehole Acoustics and Logging Consortium. Annual Report
Borehole Acoustics and Logging Consortium Annual Report 1994 EARTH RESOURCES LABORATORY Department of Earth, Atmospheric, and Planetary Sciences Massachusetts Institute of Technology Cambridge, Mass. 02139
More informationExact Seismic Velocities for VTI and HTI Media and Extended Thomsen Formulas for Stronger Anisotropies. Abstract
Submitted to : Geophysics Exact Seismic Velocities for VTI and HTI Media and Extended Thomsen Formulas for Stronger Anisotropies James G. Berryman 1, 1 University of California, Lawrence Berkeley National
More informationWe Challenges in shale-reservoir characterization by means of AVA and AVAZ
We-06-15 Challenges in shale-reservoir characterization by means of AVA and AVAZ N.C. Banik* (WesternGeco), M. Egan (WesternGeco), A. Koesoemadinata (WesternGeco) & A. Padhi (WesternGeco) SUMMARY In most
More informationEstimation of interval anisotropic attenuation from reflection data
GEOPHYSCS, VOL. 74, NO. 6 NOVEMBER-DECEMBER 9 ; P. A69 A74, 7 FGS..9/.39733 Estimation of interval anisotropic attenuation from reflection data Jyoti Behura and lya Tsvankin ABSTRACT Knowledge of interval
More informationElastic wavefield separation for VTI media
CWP-598 Elastic wavefield separation for VTI media Jia Yan and Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT The separation of wave modes from isotropic elastic wavefields is typically
More informationSeismic Exploration of Hydrocarbons in Heterogeneous Reservoirs
Seismic Exploration of Hydrocarbons in Heterogeneous Reservoirs New Theories, Methods, and Applications Jing Ba Department of Earth and Atmospheric Sciences University of Houston, Houston, Texas, US Qizhen
More information