起伏地表条件下各向异性地震波最短路径射线追踪

在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度.     

Magnetic fabric of Plio-Pleistocene sediments from the Crotone fore-arc basin: Insights on the recent tectonic evolution of the Calabrian Arc (Italy)

Low-field anisotropy of magnetic susceptibility (AMS) analyses were performed on 532 samples collected in 36 (mostly lower Pliocene to lower Pleistocene) marine clay sites from the Crotone basin, a fore-arc basin located on top of the external Calabrian accretionary wedge. The Crotone basin formed since mid-late Miocene under a predominant extensional tectonic regime, but it was influenced thereafter by complex interactions with NW–SE left-lateral strike-faults bounding the basin, which also yielded post-1.2 Ma ∼30° counterclockwise block rotations. The basin is filled by continental to marine sediments yielding one of the thickest and best-exposed Neogene succession available worldwide. The deep-marine facies – represented by blue-grey marly clays gave the best results, as they both preserved a clear magnetic fabric, and provided accurate chronology based on previously published magnetostratigraphy and calcareous plankton (i.e. foraminifers and nannofossils) biostratigraphy. Magnetic susceptibility range and rock magnetic analyses both indicate that AMS reflects paramagnetic clay matrix crystal arrangement. The fabric is predominantly oblate to triaxial, the anisotropy degree low (<1.06), and the magnetic foliation mostly subparallel to bedding. Magnetic lineation is defined in 30 out of 36 sites (where the e12 angle is <35°). By also considering local structural analysis data, we find that magnetic fabric was generally acquired during the first tectonic phases occurring after sediment deposition, thus validating its use as temporally dependent strain proxy. Although most of the magnetic lineations trend NW–SE and are orthogonal to normal faults (as observed elsewhere in Calabria), few NE–SW compressive lineations show that the Neogene extensional regime of the Crotone basin was punctuated by compressive episodes. Finally, compressive lineations (prolate magnetic fabric) documented along the strike-slip fault bounding the basin to the south support the significance of Pleistocene strike-slip tectonics. Thus the Crotone basin shows a markedly different tectonics with respect to other internal and western basins of Calabria, as it yields a magnetic fabric still dominated by extensional tectonics but also revealing arc-normal shortening episodes and recent strike-slip fault activity. The tectonics documented in the Crotone basin is compatible with a continuous upper crustal structural reorganization occurring during the SE-migration of the Calabria terrane above the Ionian subduction system.     

Anisotropic cosmological models with bulk viscosity for variable <Emphasis Type="Italic">G</Emphasis> and Λ

In this paper we have considered axially symmetric Bianchi-I, Kantowski Sachs and Bianchi-III space-time models with bulk viscosity, where the gravitational constant G and the cosmological term Λ vary with time. In Einstein equations this variation in G and Λ are taken in such a way as to preserve the energy momentum tensor. Solutions are obtained with the cosmological term varying inversely with square of time.     

The anisotropy of cosmic ray arrival directions around 10 eV

Anisotropy in the arrival directions of cosmic rays with energies above 10¹⁷ eV is studied using data from the Akeno 20 km² array and the Akeno Giant Air Shower Array (AGASA), using a total of about 114 000 showers observed over 11 years. In the first harmonic analysis, we have found a strong anisotropy of 4% around 10¹⁸ eV, corresponding to a chance probability of 0.2% after taking the number of independent trials into account. with two-dimensional analysis in right ascension and declination, this anisotropy is interpreted as an excess of showers near the directions of the Galactic Center and the Cygnus region.     

Seismic structure of the upper mantle beneath the western Philippine Sea

We investigate the seismic structure of the western Philippine Sea using two sets of seismological observations: ScS reverberations, which provide the layering framework for a regional upper mantle model, and observations of frequency-dependent phase delays for direct S waves, surface-reflected phases (sS, SS, sSS), and surface waves (R₁, G₁), which constrain the velocity and anisotropy structure within the layers. The combined data set, comprising 17 discontinuity amplitudes and layer travel times from the ScS-reverberation stack and more than 1000 frequency-dependent phase delays, was inverted for a path-averaged, radially anisotropic model. Mineralogical estimates of the bulk sound velocity and density are incorporated as complementary constraints. The final model, PHB3, is characterized by a 11.5-km thick crust, an anisotropic lid bounded by a sharp negative G discontinuity at 89 km, an anisotropic low-velocity layer extending to 166 km, a subjacent high-gradient region, and transition-zone discontinuities at depths of 408 km, 520 km and 664 km. The lid is slower than in a comparable model for the Tonga–Hawaii corridor (PA5), but also significantly thicker, requiring a compositional variation between the two regions. We explore the hypothesis that the thickness of the oceanic lid is controlled by the melting depth at the spreading centers during crust formation, and that the thicker crust and lid in the Philippine Sea results from deeper melting owing to a higher potential temperature and perhaps a higher water content in the upper mantle.     

Travel time computations using a compact eikonal equation for vertical transverse isotropic media

Eikonal solvers often have stability problems if the velocity model is mildly heterogeneous. We derive a stable and compact form of the eikonal equation for P‐wave propagation in vertical transverse isotropic media. The obtained formulation is more compact than other formulations and therefore computationally attractive. We implemented ray shooting for this new equation through a Hamiltonian formalism. Ray tracing based on this new equation is tested on both simple as well as more realistic mildly heterogeneous velocity models. We show through examples that the new equation gives travel times that coincide with the travel time picks from wave equation modelling for anisotropic wave propagation.     

Study on the Anisotropy of the Underground Medium in Hohhot

In this paper, the data of earthquake events of magnitude M_S6.0 and above produced in Hohhot Seismic Station from 2008 to 2015 and the data of M_L ≥ 1.0 seismic events from 2015 to 2016 in Horinger Seismic Station and the surrounding mobile stations in southern Hohhot are selected. Using Splitlab and SAM software, the spin-correlation method, the least-energy method and the cross-correlation coefficient method are used to analyze the teleseismic and near-seismic phases (SKS, S). The results of this study are in good agreement with the results previously obtained by other researchers. The study of teleseismic SKS splitting reflects the characteristics of the anisotropy of the upper mantle beneath Hohhot, that is, the anisotropy of the upper mantle shows NW, which reflects "fossil" Anisotropy, mainly in the continental structure of stable units and preserves the history of mantle deformation information. The crustal anisotropy reflected by the near-earthquake S-wave splitting study is similar to that of the active fault zone, trending NE as a whole and is consistent with the tectonic stress field background of the northeastern margin of Ordos block.     

Lithology-dependent seismic anisotropic amplitude variation with offset correction in transversely isotropic media

Analysis of amplitude variation with offset is an essential step for reservoir characterization. For an accurate reservoir characterization, the amplitude obtained with an isotropic assumption of the reservoir must be corrected for the anisotropic effects. The objective is seismic anisotropic amplitude correction in an effective medium, and, to this end, values and signs of anisotropic parameter differences (Δδ and Δε) across the reflection interfaces are needed. These parameters can be identified by seismic and well log data. A new technique for anisotropic amplitude correction was developed to modify amplitude changes in seismic data in transversely isotropic media with a vertical axis of symmetry. The results show that characteristics of pre-stack seismic data, that is, amplitude variation with offset gradient, can be potentially related to the sign of anisotropic parameter differences (Δδ and Δε) between two layers of the reflection boundary. The proposed methodology is designed to attain a proper fit between modelled and observed amplitude variation with offset responses, after anisotropic correction, for all possible lithofacies at the reservoir boundary. We first estimate anisotropic parameters, that is, δ and ε, away from the wells through Backus averaging of elastic properties resulted from the first pass of isotropic pre-stack seismic inversion, on input data with no amplitude correction. Next, we estimate the anisotropic parameter differences at reflection interfaces (values and signs of Δδ and Δε). We then generate seismic angle gather data after anisotropic amplitude correction using Rüger's equation for the P-P reflection coefficient. The second pass of isotropic pre-stack seismic inversion is then performed on the amplitude-corrected data, and elastic properties are estimated. Final outcome demonstrates how introduced methodology helps to reduce the uncertainty of elastic property prediction. Pre-stack seismic inversion on amplitude-corrected seismic data results in more accurate elastic property prediction than what can be obtained from non-corrected data. Moreover, a new anisotropy attribute (ν) is presented for improvement of lithology identification.