GLORIA

GEOMAR Library Ocean Research Information Access

feed icon rss

Your email was sent successfully. Check your inbox.

An error occurred while sending the email. Please try again.

Proceed reservation?

Export
  • 1
    Electronic Resource
    Electronic Resource
    Oxford, UK : Blackwell Publishing Ltd
    Geophysical journal international 101 (1990), S. 0 
    ISSN: 1365-246X
    Source: Blackwell Publishing Journal Backfiles 1879-2005
    Topics: Geosciences
    Notes: The anisotropic linear viscoelastic rheological relation constitutes a suitable model for describing the variety of phenomena which occur in seismic wavefields. This rheology, known also as Boltzmann's superposition principle, expresses the stress as a time convolution of a fourth rank tensorial relaxation function with the strain tensor.The first problem is to establish the time dependence of the relaxation tensor in a general and consistent way. Two kernels based on the general standard linear solid are identified with the mean stress and with the deviatoric components of the stress tensor in a given coordinate system, respectively. Additional conditions are that in the elastic limit the relaxation matrix must give the elasticity matrix, and in the isotropic limit the relaxation matrix must approach the isotropic-viscoelastic matrix. The resulting rheological relation provides the framework for incorporating anelasticity in time-marching methods for computing synthetic seismograms. Through a plane wave analysis of the anisotropic-viscoelastic medium, the phase, group and energy velocities are calculated in function of the complex velocity, showing that those velocities are in general different from each other. For instance, the energy velocity which represents the wave surface, is different from the group velocity unlike in the anisotropic-elastic case. The group velocity loses its physical meaning at the cusps where singularities appear. Each frequency component of the wavefield has a different non-spherical wavefront. Moreover, the quality factors for the different propagating modes are not isotropic. Examples of these physical quantities are shown for transversely isotropic-viscoelastic clayshale and sandstone.As in the isotropic-viscoelastic case, Boltzmann's superposition principle is implemented in the equation of motion by defining memory variables which circumvent the convolutional relation between stress and strain. The numerical problem is solved by using a new time integration technique specially designed to deal with wave propagation in linear viscoelastic media. As a first application snapshots and synthetic seismograms are computed for 2-D transversely isotropic-viscoelastic clayshale and sandstone which show substantial differences in amplitude, waveform and arrival time with the results given by the isotropic and elastic rheologies.
    Type of Medium: Electronic Resource
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 2
    Electronic Resource
    Electronic Resource
    Springer
    Pure and applied geophysics 157 (2000), S. 811-827 
    ISSN: 1420-9136
    Keywords: Key Words: Wave velocities, porous medium, clay content, poroelasticity.
    Source: Springer Online Journal Archives 1860-2000
    Topics: Geosciences , Physics
    Notes: Abstract —The feasibility of modeling elastic properties of a fluid-saturated sand-clay mixture rock is analyzed by assuming that the rock is composed of macroscopic regions of sand and clay. The elastic properties of such a composite rock are computed using two alternative schemes.¶The first scheme, which we call the composite Gassmann (CG) scheme, uses Gassmann equations to compute elastic moduli of the saturated sand and clay from their respective dry moduli. The effective elastic moduli of the fluid-saturated composite rock are then computed by applying one of the mixing laws commonly used to estimate elastic properties of composite materials.¶In the second scheme which we call the Berryman-Milton scheme, the elastic moduli of the dry composite rock matrix are computed from the moduli of dry sand and clay matrices using the same composite mixing law used in the first scheme. Next, the saturated composite rock moduli are computed using the equations of Brown and Korringa, which, together with the expressions for the coefficients derived by Berryman and Milton, provide an extension of Gassmann equations to rocks with a heterogeneous solid matrix.¶For both schemes, the moduli of the dry homogeneous sand and clay matrices are assumed to obey the Krief’s velocity-porosity relationship. As a mixing law we use the self-consistent coherent potential approximation proposed by Berryman.¶The calculated dependence of compressional and shear velocities on porosity and clay content for a given set of parameters using the two schemes depends on the distribution of total porosity between the sand and clay regions. If the distribution of total porosity between sand and clay is relatively uniform, the predictions of the two schemes in the porosity range up to 0.3 are very similar to each other. For higher porosities and medium-to-large clay content the elastic moduli predicted by CG scheme are significantly higher than those predicted by the BM scheme.¶This difference is explained by the fact that the BM model predicts the fully relaxed moduli, wherein the fluid can move freely between sand and clay regions. In contrast, the CG scheme predicts the no-flow or unrelaxed moduli. Our analysis reveals that due to the extremely low permeability of clays, at seismic and higher frequencies the fluid has no time to move between sand and clay regions. Thus, the CG scheme is more appropriate for clay-rich rocks.
    Type of Medium: Electronic Resource
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 3
    Electronic Resource
    Electronic Resource
    New York, NY [u.a.] : Wiley-Blackwell
    International Journal for Numerical and Analytical Methods in Geomechanics 22 (1998), S. 263-275 
    ISSN: 0363-9061
    Keywords: non-ideal interface ; crack scattering ; domain decomposition ; collocation methods ; Engineering ; Civil and Mechanical Engineering
    Source: Wiley InterScience Backfile Collection 1832-2000
    Topics: Architecture, Civil Engineering, Surveying , Geosciences
    Notes: This work presents a numerical algorithm for solving crack scattering in a transversely isotropic medium whose symmetry axis is perpendicular to the crack surface. The crack is modelled as boundary discontinuities in the displacement u and the particle velocity v, of the stresses [κu+ζv], where the brackets denote discontinuities across the interface. The specific stiffness κ introduces frequency-dependence and phase changes in the interface response and the specific viscosity ζ is related to the energy loss.The numerical method is based on a domain decomposition technique that assignes a different mesh to each side of the interface, that includes the crack plane. As stated above, the effects of the crack on wave propagation are modelled through the boundary conditions, that require a special boundary treatment based on characteristic variables. The algorithm solves the particle velocity-stress wave equations and two additional first-order differential equations (two-dimensional case) in the displacement discontinuity. For each mesh, the spatial derivatives normal to the interface are solved by the Chebyshev method, and the spatial derivatives parallel to the interface are computed with the Fourier method. They allow a highly accurate implementation of the boundary conditions and computation of the spatial derivatives, and an optimal discretization of the model space. Moreover, the algorithm allows general material variability. © 1998 John Wiley & Sons, Ltd.
    Additional Material: 9 Ill.
    Type of Medium: Electronic Resource
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 4
    Publication Date: 2017-07-21
    Description: High-resolution seismic experiments, employing arrays of closely spaced, four-component ocean-bottom seismic recorders, were conducted at a site off western Svalbard and a site on the northern margin of the Storegga slide, off Norway to investigate how well seismic data can be used to determine the concentration of methane hydrate beneath the seabed. Data from P-waves and from S-waves generated by P–S conversion on reflection were inverted for P- and S-wave velocity (Vp and Vs), using 3D travel-time tomography, 2D ray-tracing inversion and 1D waveform inversion. At the NW Svalbard site, positive Vp anomalies above a sea-bottom-simulating reflector (BSR) indicate the presence of gas hydrate. A zone containing free gas up to 150-m thick, lying immediately beneath the BSR, is indicated by a large reduction in Vp without significant reduction in Vs. At the Storegga site, the lateral and vertical variation in Vp and Vs and the variation in amplitude and polarity of reflectors indicate a heterogeneous distribution of hydrate that is related to a stratigraphically mediated distribution of free gas beneath the BSR. Derivation of hydrate content from Vp and Vs was evaluated, using different models for how hydrate affects the seismic properties of the sediment host and different approaches for estimating the background-velocity of the sediment host. The error in the average Vp of an interval of 20-m thickness is about 2.5%, at 95% confidence, and yields a resolution of hydrate concentration of about 3%, if hydrate forms a connected framework, or about 7%, if it is both pore-filling and framework-forming. At NW Svalbard, in a zone about 90-m thick above the BSR, a Biot-theory-based method predicts hydrate concentrations of up to 11% of pore space, and an effective-medium-based method predicts concentrations of up to 6%, if hydrate forms a connected framework, or 12%, if hydrate is both pore-filling and framework-forming. At Storegga, hydrate concentrations of up to 10% or 20% were predicted, depending on the hydrate model, in a zone about 120-m thick above a BSR. With seismic techniques alone, we can only estimate with any confidence the average hydrate content of broad intervals containing more than one layer, not only because of the uncertainty in the layer-by-layer variation in lithology, but also because of the negative correlation in the errors of estimation of velocity between adjacent layers. In this investigation, an interval of about 20-m thickness (equivalent to between 2 and 5 layers in the model used for waveform inversion) was the smallest within which one could sensibly estimate the hydrate content. If lithological layering much thinner than 20-m thickness controls hydrate content, then hydrate concentrations within layers could significantly exceed or fall below the average values derived from seismic data.
    Type: Article , PeerReviewed
    Format: text
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 5
    Publication Date: 2014-01-09
    Description: We present a novel approach for the comprehensive, flexible and accurate simulation of poroelastic wave propagation in 3-D cylindrical coordinates. An important application of this method is the realistic modelling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, as of yet largely unresolved, problem in exploration geophysics. To this end, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the centre followed by the mudcake and/or casing, and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction and Fourier expansions in the vertical and azimuthal directions as well as a Runge–Kutta integration scheme for the time evolution. Trigonometric interpolation and a domain decomposition method based on the method of characteristics are used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces as well as to reduce the number of gridpoints in the innermost domain for computational efficiency. We apply this novel modelling approach to the particularly challenging scenario of near-surface borehole environments. To this end, we compare 3-D heterogeneous and corresponding rotationally invariant simulations, assess the sensitivity of Stoneley waves to formation permeability in the presence of a casing and evaluate the effects of an excavation damage zone behind a casing on sonic log recordings. Our results indicate that only first arrival times of fast modes are reasonably well described by rotationally invariant approximations of 3-D heterogenous media. We also find that Stoneley waves are indeed remarkably sensitive to the average permeability behind a perforated PVC casing, and that the presence of an excavation damage zone behind a casing tends to dominate the overall signature of recorded seismograms.
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 6
    Publication Date: 2016-02-08
    Description: Wave-induced local fluid flow is known as a key mechanism to explain the intrinsic wave dissipation in fluid-saturated rocks. Understanding the relationship between the acoustic properties of rocks and fluid patch distributions is important to interpret the observed seismic wave phenomena. A triple-layer patchy (TLP) model is proposed to describe the P -wave dissipation process in a double-porosity media saturated with two immiscible fluids. The double-porosity rock consists of a solid matrix with unique host porosity and inclusions which contain the second type of pores. Two immiscible fluids are considered in concentric spherical patches, where the inner pocket and the outer sphere are saturated with different fluids. The kinetic and dissipation energy functions of local fluid flow (LFF) in the inner pocket are formulated through oscillations in spherical coordinates. The wave propagation equations of the TLP model are based on Biot's theory and the corresponding Lagrangian equations. The P -wave dispersion and attenuation caused by the Biot friction mechanism and the local fluid flow (related to the pore structure and the fluid distribution) are obtained by a plane-wave analysis from the Christoffel equations. Numerical examples and laboratory measurements indicate that P -wave dispersion and attenuation are significantly influenced by the spatial distributions of both, the solid heterogeneity and the fluid saturation distribution. The TLP model is in reasonably good agreement with White's and Johnson's models. However, differences in phase velocity suggest that the heterogeneities associated with double-porosity and dual-fluid distribution should be taken into account when describing the P -wave dispersion and attenuation in partially saturated rocks.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 7
    Publication Date: 2016-05-11
    Description: The present study evaluates the capacity of the Boom Clay as a host rock for disposal purposes, more precisely its seismic characterization, which may assess its long-term performance to store radioactive wastes. Although the formation is relatively uniform and homogeneous, there are embedded thin layers of septaria (carbonates) that may affect the integrity of the Boom Clay. Therefore, it is essential to locate these geobodies. The seismic data to characterize the Boom Clay has been acquired at the Kruibeke test site. The inversion, which allowed us to obtain the anisotropy parameters and seismic velocities of the clay, is complemented with further information such as log and laboratory data. The attenuation properties have been estimated from equivalent formations (having similar composition and seismic velocities). The inversion yields quite consistent results although the symmetry of the medium is unusual but physically possible, since the anisotropy parameter is negative. According to a time-domain calculation of the energy velocity at four frequency bands up to 900 Hz, velocity increases with frequency, a behaviour described by the Zener model. Then, we use this model to describe anisotropy and anelasticity that are implemented into the equation of motion to compute synthetic seismograms in the space–time domain. The technique is based on memory variables and the Fourier pseudospectral method. We have computed reflection coefficients of the septaria thin layer. At normal incidence, the P -wave coefficient vanishes at specific thicknesses of the layer and there is no conversion to the S wave. For example, calculations at 600 Hz show that for thicknesses of 1 m the septarium can be detected more easily since the amplitudes are higher (nearly 0.8). Converted PS waves have a high amplitude at large offsets (between 30° and 80°) and can be useful to identify the target on this basis. Moreover, we have investigated the effect of septaria embedded in the Boom Clay with several simulations, by considering a lateral partial continuity of the calcareous thin inclusions. The simulations with layers of calcareous material show continuity of the reflections even when the percentage of carbonate within the layer is very small (5–15 per cent), while for low content of the calcareous material, isolated septaria boulders generate diffraction events. We have also simulated the stacked seismic section obtained from processing of the field data. The matching between the field and synthetic sections is acceptable.
    Keywords: Seismology
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 8
    Publication Date: 2014-02-13
    Description: We derive a time-domain differential equation for modelling seismic wave propagation in constant- Q viscoelastic media based on fractional spatial derivatives, specifically Laplacian differential operators of fractional order. The stress–strain relation is derived from the classical equation expressed in terms of fractional time derivatives. The new formulation has the advantage of not requiring additional field variables that increase the computer time and storage significantly. The spatial derivatives are calculated with a generalization of the Fourier pseudospectral method to the fractional-derivative case. The accuracy of the numerical solution is verified against an analytical solution in a homogeneous medium. An example shows that the proposed wave equation describes the constant- Q attenuation and velocity dispersion behaviour observed in Pierre Shale. Finally, we consider a plane-layer model and the Marmousi model to show how the new formulation applies to inhomogeneous media.
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 9
    Publication Date: 2014-04-16
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
  • 10
    Publication Date: 2014-04-16
    Print ISSN: 0956-540X
    Electronic ISSN: 1365-246X
    Topics: Geosciences
    Published by Oxford University Press on behalf of The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
    Location Call Number Limitation Availability
    BibTip Others were also interested in ...
Close ⊗
This website uses cookies and the analysis tool Matomo. More information can be found here...