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). 相似文献
As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques
have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K2[6](Al3Mg)[4](Si7Al)O20(OH)4. Values of the atomic exchange interaction parameters Jn used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied
and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral
sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst
we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets
alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the
T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range
order was observed.
Received: 8 August 2002 / Accepted: 14 February 2003
Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed
on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility. 相似文献
The synthetic amphibole Na0.95(Li0.95Mg1.05)Mg5Si8O22(OH)2 was studied in situ at high-T, using IR OH-stretching spectroscopy and synchrotron X-ray powder diffraction. At room-T the sample has P21/m symmetry, as shown by the FTIR spectrum. It shows in the OH region two well-defined and intense absorptions at 3,748 and
3,712 cm−1, respectively, and two minor bands at 3,667 and 3,687 cm−1. The main bands are assigned to the two independent O–H groups in the primitive structure. The two minor bands evidencing
the presence of small amount of vacant A-site (A□0.05). With increasing T, these bands shift continuously and merge into a unique absorption at high temperature. A change as a function of increasing
T is revealed by the evolution of the refined unit-cell parameters, whose trend shows a transition to C2/m at about 320–330°C. The spontaneous scalar strain, fitted with a tricritical 2–6 Landau potential, gives a Tc of 325(10)°C (β parameter = 0.27). Comparison with the second-order P21/m ⇔ C2/m phase transition at 255°C for synthetic amphibole ANa0.8B(Na0.8Mg1.2)CMg5Si8O22(OH)2 indicates that the substitution of Na with Li at the B-sites strongly affects the thermodynamic character and the Tc of the phase transition. The comparison of LNMSH amphiboles with cummingtonitic ones shows that the high-T thermodynamic behaviour is affected by A-site occupancy. 相似文献