High-resolution seismic data (onshore and offshore), geophysical borehole data as well as detailed lithofacies from airlift boreholes were acquired in northern Netherlands on and around the island of Ameland. Marine and land seismic data combined with information from land boreholes have been explored with the objective of providing a sedimentary model. Qualitative seismic facies analysis of the valley fill commonly shows a thin unit with high amplitude reflectors at the base. Thick units of variable seismic facies (transparent to high amplitude) occur higher up in the sequence. Onlap is common at mid–upper levels within the sandy valley fill (with clay in mm layering), and a transparent seismic facies, corresponding to firm clays, is common at the top. Almost all lithological unit boundaries recognised within core parameters correspond with seismic unconformities within error margins. Subunits contain multiple cyclical trends in gamma ray and grain size. Cyclical trends show lower order fluctuations in gamma radiation on a scale of less than 1 m. Gamma-ray pattern variability between units, e.g. in general coarsening-up or fining-up units, suggests migration of subaqueous outwash fans or ice margin fluctuations. Seismic results could support a headward excavation and backfilling process suggested by Praeg [Morphology, stratigraphy and genesis of buried Elsterian tunnel valleys in the southern North Sea basin [PhD thesis]: University of Edinburgh, 207 pp.; Journal of Applied Geophysics, (this volume)] as being responsible for the formation of buried valleys. On a lithological scale, a more complicated, detailed and cyclical pattern arises. Catastrophic processes are considered unlikely as being responsible for the infill sequence because of the observed small-scale facies variability and because of the presence of diamicton layers. Diamicton layers at the base of basal unconformities as well as higher in sequence could suggest subglacial deformation by grounded ice before and during the valley-fill process. 相似文献
We prove that rays in linearly elastic anisotropic nonuniform media obey Fermat's principle of stationary traveltime. First, we formulate the concept of rays, which emerges from the Hamilton equations. Then, we show that these rays are solutions of the variational problem stated by Fermat's principle. This proof is valid for all rays except the ones associated with infection points on the phase-slowness surface. 相似文献
Transverse isotropy (TI) with a vertical symmetry axis (VTI) often provides an appropriate earth model for prestack imaging of steep-dip reflection seismic data. Exact P-wave and SV-wave phase velocities in VTI media are described by complicated equations requiring four independent parameters. Estimating appropriate multiparameter earth models can be difficult and time-consuming, so it is often useful to replace the exact VTI equations with simpler approximations requiring fewer parameters. The accuracy limits of different previously published VTI approximations are not always clear, nor is it always obvious how these different approximations relate to each other. Here I present a systematic framework for deriving a variety of useful VTI approximations. I develop first a sequence of well-defined approximations to the exact P-wave and SV-wave phase velocities. In doing so, I show how the useful but physically questionable heuristic of setting shear velocities identically to zero can be replaced with a more precise and quantifiable approximation. The key here to deriving accurate approximations is to replace the stiffness a13 with an appropriate factorization in terms of velocity parameters. Two different specific parameter choices lead to the P-wave approximations of Alkhalifah (Geophysics 63 (1998) 623) and Schoenberg and de Hoop (Geophysics 65 (2000) 919), but there are actually an infinite number of reasonable parametrizations possible. Further approximations then lead to a variety of other useful phase velocity expressions, including those of Thomsen (Geophysics 51 (1986) 1954), Dellinger et al. (Journal of Seismic Exploration 2 (1993) 23), Harlan (Stanford Exploration Project Report 89 (1995) 145), and Stopin (Stopin, A., 2001. Comparison of v(θ) equations in TI medium. 9th International Workshop on Seismic Anisotropy). Each P-wave phase velocity approximation derived this way can be paired naturally with a corresponding SV-wave approximation. Each P-wave or SV-wave phase velocity approximation can then be converted into an equivalent dispersion relation in terms of horizontal and vertical slownesses. A simple heuristic substitution also allows each phase velocity approximation to be converted into an explicit group velocity approximation. From these, in turn, travel time or moveout approximations can also be derived. The group velocity and travel time approximations derived this way include ones previously used by Byun et al. (Geophysics 54 (1989) 1564), Dellinger et al. (Journal of Seismic Exploration 2 (1993) 23), Tsvankin and Thomsen (Geophysics 59 (1994) 1290), Harlan (89 (1995) 145), and Zhang and Uren (Zhang, F. and Uren, N., 2001. Approximate explicit ray velocity functions and travel times for P-waves in TI media. 71st Annual International Meeting, Society of Exploration Geophysicists, Expanded Abstracts, 106–109). 相似文献
Seismic reflection methods provide continuous access both to stratigraphy (vertical) and to subsurface morphology (horizontal), for which the scales of interest may differ by orders of magnitude. Seismic surveys of Quaternary successions have generally sought to optimise vertical resolution, through the use of higher source frequency content. Here, I show that low-frequency bandwidth is not necessarily a limiting factor for the seismic resolution of glacigenic morpho-sedimentary features. Observations are presented from a buried network of large mid-Pleistocene (Elsterian) tunnel-valleys in the southern North Sea Basin, across a 100×130 km study area with water depths less than 30 m. Low-frequency 2D and 3D seismic multi-channel data, acquired for deeper hydrocarbon exploration, are compared with previously available high-frequency single- and multi-channel profiles (5–15 km grid spacing). The low-frequency data contribute to a new understanding of the basal morphology and fill stratigraphy of the tunnel-valleys, in part due to higher data densities (≥1 km grid spacing), but also to improved imaging of reflectors at depth. The tunnel-valleys are seen to be overdeepened troughs, shallow (≤0.5 km) relative to their widths (≤6 km). The basal unconformity defines a series of arborescent elements, convergent to the south; erosional overlap by younger elements to the north has resulted in anastomosing patterns in places. The fill is dominated by axially downlapping clinoforms, descending to the north, onlapped and overlain by subhorizontal reflectors. Well data show that sand-dominated glaciofluvial sediments are overlain by glaciolacustrine to marine muds. Better definition of the clinoforms on low- versus high-frequency multi-channel data is suggested to reflect the coarse spatial scale of the backset glaciofluvial strata. The results support a simple interpretation of time-transgressive tunnel-valley formation by coeval glaciofluvial erosion and backfill beneath the outer tens of kilometres of the northward receding Elsterian ice sheet margin. Comparable submarginal interpretations have been proposed for drainage features (tunnel-valleys and eskers) of the last deglaciation of both northern Europe and North America using integrated geomorphologic and stratigraphic methods. Commercial 2D and 3D seismic data are widely available from exploration areas such as the North Sea and are argued to constitute an underexploited resource for Quaternary research. 相似文献
A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys. 相似文献
The devastating earthquake on 26 January 2001 at Bhuj, India, resulted in large-scale death and destruction of properties of several million US dollars. The moment magnitude of the earthquake was 7.7 and its maximum focal intensity exceeded X in MM scale. The rate of aftershocks of this earthquake, recorded at Gauribidanur seismic array station (GBA), shows a monotonic decay with time superposed with oscillations. For the Indian continent the Lg phase is a prominent arrival at regional distances. The estimate of Lg amplitude is obtained by optimally fitting the Lg wave train to a exponential decay curve. The logarithm of these amplitudes and logarithm of root mean square (rms) value of actual amplitudes of the Lg are calibrated with USGS mb to create a local mbLg magnitude scale. The energy released from these aftershocks is calculated from the rms value of Lg phase. The plot of cumulative energy release with time follows the power law of the form tp, superposed with oscillations. The exponent of the power law, p, is estimated both by a time-window scanning method and by an interpolation method. The value of p is 0.434 for time-window scanning method and 0.432 for the interpolation method. The predominant periods found in the oscillatory part of the cumulative energy, obtained by differencing the observed from the power law fit, are 10.6, 7.9, 5.4, 4.6 and 3.5 h for time-window scanning method. The corresponding periods for interpolation method are 13.4, 11.5, 7.4, 4.2, 3.5, 2.6 and 2.4 h. 相似文献
The design of a drainage system for a roofing slate quarry was implemented by the enhancement of discharge peak estimation, and the uncertainty inevitably associated with the engineering model was reduced.
The development of a topographical, geological, and vegetation cover database developed from a Geographical Information System (GIS) allowed for the definition of the drainage network for a hydraulic system, along with the calculation of the runoff coefficient. This is applied to the digital model of accumulated flow (DMF) as a weight correction coefficient, using a matrix-based model at 5×5 m resolution. The new digital model of corrected accumulated flow (DMCF) is the result of combining the thematic maps with the map of slope <3%, which was previously created from the slope model. It is demonstrated that this new model allows to apply the “Rational Method” on cartographic units defined by the GIS.
The DMCF is compared with other traditional applications of the Rational Method based on the calculation of the discharge peak considering: (1) the drainage basin as a single watershed or (2) defining an average runoff coefficient in each sub-watershed. Both approaches have bigger discharge peaks than those obtained by the DMCF since the slope, lithology, and vegetation cover have average values, and the runoff coefficient is poorly defined, increasing the uncertainty in the discharge peak. 相似文献
Shallow seismic measurements in harzburgite from the Oman ophiolite performed in a zone where the maximum horizontal anisotropy is expected (vertical foliation and horizontal lineation) point to a dominant dependence of seismic properties on fracturing.
Optical microscopy studies show that microcracks are guided by the serpentine (lizardite) penetrative network oriented subparallel to the harzburgite foliation and subperpendicular to the mineral lineation, and that serpentine (lizardite) vein filling has a maximum concentration of (001) planes parallel to the veins walls. The calculated elastic properties of the oriented alteration veins filled with serpentine in an anisotropic matrix formed by oriented crystals of olivine and orthopyroxene are compared with seismic velocities measured on hand specimens.
Laboratory ultrasonic data indicate that open microcracks are closed at 100 MPa pressure, e.g. (J. Geophys. Res. 65, (1960) 1083) and (Proc. ODP Sci. Results Leg 118, (1990) 227). Above this pressure, laboratory measurements and modeling show that P-compressional and S-shear wave velocities are mainly controlled by the mineral preferred orientation. Veins sealed with serpentine are effective in slightly lowering P and S velocities and increasing anisotropy. The penetrative lizardite network does not affect directly the geometry of seismic anisotropy, but contributes indirectly in the fact that this network controls the microcrack orientations.
Comparison between seismic measurements of peridotite and gabbro in the same conditions suggest that P- and S-waves anisotropies are a possible discriminating factor between the two lithologies in the suboceanic lithosphere. 相似文献