饱和土一维简谐响应解析解的求解和应用：II 应用

将相关论文[1]中得到的饱和土一维简谐响应解析解，应用到饱和土中两类压缩波的独立作用、饱和土中波的传播速度、u-p方程的适用范围等研究中。在两类压缩波的独立作用研究中，克服了以前文献中的数学不严密性。研究了饱和土各参数对三种体波波速的影响。通过比较u-w和u-p方程的解答，定量地说明了u-p方程的适用范围。     

马来西亚沙捞越邦达、什兰江控矿角砾岩筒构造对比研究及其找矿意义

马来西亚沙捞越邦达、什兰江控矿角砾岩筒构造均为典型的浅成—半浅成角砾岩筒构造，角砾岩及控矿类型具有明显的垂向分带特点。邦达控矿角砾岩筒构造发育角砾岩筒上段的塌陷角砾岩，角砾岩筒剥蚀深度较浅，控制典型的高硫化型浅成低温金铜矿床，且其深部具有较大的找矿潜力。什兰江控矿角砾岩筒构造发育角砾岩筒下段的爆发角砾岩和中段的流化角砾岩，分别控制斑岩型金矿化和低硫化型浅成低温热液型金矿化，角砾岩筒剥蚀深度较大，对矿床的保存和找矿极为不利。     

Modeling of Rayleigh waves dispersion in the Sannio region (Southern Italy) from seismic active refraction data

Short period surface waves, recorded during a seismic refractionsurvey in the Sannio region (Southern Italy), have been modeled to infera shallow velocity model for the area. Based on the decrease of resolutionwith depth, due to the bias on group velocity estimates arising frominterference of the Rayleigh waves with higher modes, we carried out aprocedure of fitting, with synthetic seismograms, of selected filtered traceswith a gaussian filter, having a width at half height equal to 1 Hz and acentral frequency lying in the range [1,4] Hz. We estimated the likelihoodbetween synthetic and observed seismograms by measuring their semblance.In this way we were able to infer a more refined local velocity modelcharacterized by a high Vp and Vs vertical gradient in the sedimentarycover. Two ad hoc resolution studies, based on group velocity andamplitude data respectively, indicate that the local velocity model is a goodvelocity model also for the entire studied area. The increase in the numberof available data when using amplitude information allows us to make amore selective choice in the model parameter space (Vp and Vs of eachlayer) and to solve for the Vp/Vs ratio. The inferred Vp velocity in thehalf-space is equal to 2.8 km/s. This value is in excellent agreement withthat inferred by other authors (3 km/s) by modeling P-wave travel timevs. distance. The best-fit model furnish low Vp/Vs for the sedimentarycover so indicating a high degree of the sediment's compaction in thestudied area. The inferred shallow high-velocity gradient indicates thatthe shallow sedimentary layer in the area could trap and focus the energytraveling into it.     

Practical VTI approximations: a systematic anatomy

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 a₁₃ 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 evidence for deep low-velocity anomalies in the transition zone beneath West Antarctica

We present a three-dimensional (3D) SV-wave velocity model of the upper mantle beneath the Antarctic plate constrained by fundamental and higher mode Rayleigh waves recorded at regional distances. The good agreement between our results and previous surface wave studies in the uppermost 200 km of the mantle confirms that despite strong differences in data processing, modern surface wave tomographic techniques allow to produce consistent velocity models, even at regional scale. At greater depths the higher mode information present in our data set allows us to improve the resolution compared to previous regional surface wave studies in Antarctica that were all restricted to the analysis of the fundamental mode. This paper is therefore mostly devoted to the discussion of the deeper part of the model. Our seismic model displays broad domains of anomalously low seismic velocities in the asthenosphere. Moreover, we show that some of these broad, low-velocity regions can be more deeply rooted. The most remarkable new features of our model are vertical low-velocity structures extending from the asthenosphere down to the transition zone beneath the volcanic region of Marie Byrd Land, West Antarctica and a portion of the Pacific-Antarctic Ridge close to the Balleny Islands hotspot. A deep low-velocity anomaly may also exist beneath the Ross Sea hotspot. These vertical structures cannot be explained by vertical smearing of shallow seismic anomalies and synthetic tests show that they are compatible with a structure narrower than 200 km which would have been horizontally smoothed by the tomographic inversion. These deep low-velocity anomalies may favor the existence of several distinct mantle plumes, instead of a large single one, as the origin of volcanism in and around West Antarctica. These hypothetical deep plumes could feed large regions of low seismic velocities in the asthenosphere.     

The submarine eruption of the Bombarda volcano,Milos Island,Cyclades, Greece

The Milos volcanic field includes a well-exposed volcaniclastic succession which records a long history of submarine explosive volcanism. The Bombarda volcano, a rhyolitic monogenetic center, erupted ∼1.7 Ma at a depth <200 m below sea level. The aphyric products are represented by a volcaniclastic apron (up to 50 m thick) and a lava dome. The apron is composed of pale gray juvenile fragments and accessory lithic clasts ranging from ash to blocks. The juvenile clasts are highly vesicular to non-vesicular; the vesicles are dominantly tube vesicles. The volcaniclastic apron is made up of three fades: massive to normally graded pumice-lithic breccia, stratified pumice-lithic breccia, and laminated ash with pumice blocks. We interpret the apron beds to be the result of water-supported, volcaniclastic mass-How emplacement, derived directly from the collapse of a small-volume, subaqueous eruption column and from syn-eruptive, down-slope resedimentation of volcaniclastic debris. During this eruptive phase, the activity could have involved a complex combination of phreatomagmatic explosions and minor submarine effusion. The lava dome, emplaced later in the source area, is made up of flow-banded lava and separated from the apron by an obsidian carapace a few meters thick. The near-vertical orientation of the carapace suggests that the dome was intruded within the apron. Remobilization of pyroclastic debris could have been triggered by seismic activity and the lava dome emplacement. Published online: 30 January 2003 Editorial responsibility: J. McPhie     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

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.