广东及邻区地震活动熵值的异常特征

应用文献」１「的方法，对广东及部分邻近地区进行了熵值计算和分析处理，结果表明，时间熵，频度熵和强度熵在ＭＬ≥４．７级地震之前，大都存在较明显的低值变化，即存在减熵现象，地震大多数以生在减熵的中，后期，少数发生在减熵初期。三种熵值变化基本同步，震级越大，熵值波动越大，其中尤以强度熵变化幅度最大，频度熵次之，时间熵最弱。     

Magnetic anisotropy produced by magma flow: theoretical model and experimental data from Ferrar dolerite sills (Antarctica)

Volcanic rocks forming sills, dykes or lava flows may display a magnetic anisotropy derived from the viscous flow during their emplacement. We model a sill as a steady-state flow of a Bingham fluid, driven by a pressure gradient in a horizontal conduit. The magma velocity as a function of depth is calculated from the motion and constitutive equations. Vorticity and strain rate are determined for a reference system moving with the fluid. The angular velocity and the orientation of an ellipsoidal magnetic grain immersed in the fluid are calculated as functions of time or strain. Magnetic susceptibility is then calculated for a large number of grains with a uniform distribution of initial orientations. It is shown that the magnetic lineation oscillates in the vertical plane through the magma flow direction, and that the magnetic foliation plane changes periodically from horizontal to vertical. The results are compared with the magnetic fabric of Ferrar dolerite sills (Victoria Land, East Antarctica) derived from low-field susceptibility measurements.     

Shear-wave anisotropy: spatial and temporal variations in time delays at Parkfield, Central California

Shear-wave splitting is analysed on data recorded by the High Resolution Seismic Network (HRSN) at Parkfield on the San Andreas fault, Central California, during the three-year period 1988-1990. Shear-wave polarizations either side of the fault are generally aligned in directions consistent with the regional horizontal maximum compressive stress, at some 70° to the fault strike, whereas at station MM in the immediate fault zone, shear-wave polarizations are aligned approximately parallel to the fault. Normalized time delays at this station are found to be about twice as large as those in the rock mass either side. This suggests that fluid-filled cracks and fractures within the fault zone are elastically or seismically different from those in the surrounding rocks, and that the alignment of fault-parallel shear-wave polarizations are associated with some fault-specific phenomenon.
Temporal variations in time delays between the two split shear-waves before and after a M_L = 4 earthquake can be identified at two stations with sufficient data: MM within the fault zone and VC outside the immediate fault zone. Time delays between faster and slower split shear waves increase before the M_L = 4 earthquake and decrease near the time of the event. The temporal variations are statistically significant at 68 per cent confidence levels. Earthquake doublets and multiplets also show similar temporal variations, consistent with those predicted by anisotropic poroelasticity theory for stress modifications to the microcrack geometry pervading the rock mass. This study is broadly consistent with the behaviour observed before three other earthquakes, suggesting that the build-up of stress before earthquakes may be monitored and interpreted by the analysis of shear-wave splitting.     

P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation

P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.     

Seismic anisotropy within the uppermost mantle of southern Germany

This paper presents an updated interpretation of seismic anisotropy within the uppermost mantle of southern Germany. The dense network of reversed and crossing refraction profiles in this area made it possible to observe almost 900 traveltimes of the P_n phase that could be effectively used in a time-term analysis to determine horizontal velocity distribution immediately below the Moho. For 12 crossing profiles, amplitude ratios of the P_n phase compared to the dominant crustal phase were utilized to resolve azimuthally dependent velocity gradients with depth. A P -wave anisotropy of 3–4 per cent in a horizontal plane immediately below the Moho at a depth of 30 km, increasing to 11 per cent at a depth of 40 km, was determined. For the axis of the highest velocity of about 8.03 km s⁻¹ at a depth of 30 km a direction of N31°F was obtained. The azimuthal dependence of the observed P_n amplitude is explained by an azimuth-dependent sub-Moho velocity gradient decreasing from 0.06 s⁻¹ in the fast direction to 0 s⁻¹ in the slow direction of horizontal P -wave velocity. From the seismic results in this study a petrological model suggesting a change of modal composition and percentage of oriented olivine with depth was derived.     

Three-dimensional Rayleigh hysteresis of oriented core samples from the German Continental Deep Drilling Program: susceptibility tensor, Rayleigh tensor, three-dimensional Rayleigh law

Rayleigh hysteresis, as defined by the well-known Rayleigh relations, has been observed not only when magnetization of pyrrhotite-bearing KTB-samples is measured in parallel to a weak dc magnetic field, but also in experiments where field and measuring directions have been adjusted strictly perpendicularly to each other. Nine-tupels of independent Rayleigh hysteresis loops could thus be compiled. Their characteristic coefficients X ^ijk of initial susceptibility together with the Rayleigh loss coefficients α^jk have been proved to determine completely the samples' weak-field magnetic anisotropy. Interpreting the coefficient matrices ( X ^ijk) and (α^jk) as the tensor of initial susceptibility and the Rayleigh tensor, respectively, generalization of the isotropic Rayleigh relations in terms of corresponding tensor relationships has been suggested for the anisotropic case. Application to the KTB samples showed 3-D Rayleigh hysteresis measurements to be an excellent tool for rock magnetic analysis in terms of ore content and crystalline texture. In particular, a magnetocrystalline double texture of the basal planes of pyrrhotite precipitates and their [1120] directions of easy magnetization have been clearly detected. Surprisingly, the welt-known theorem α= const. X ²_I, formulated by Néel (1942) for the isotropic case, has been found to hold true even in tensor generalization (α_jk) = const ( X ²_jk). To reach sufficient resolution for the measurements performed, a sensitive vibrating coil magnetometer (VCM) has been developed.     

Conditional curvilinear stochastic simulation using pixel-based algorithms

In geology, structures displaying differing local directions of continuity are widespread, a typical example being a flusial depositional system. Conventional pixel-based geostatistical algorithms, may fail to reproduce such curvilinear structures. Conversely, object-based algorithms can reproduce curvilinear shapes but are difficult to condition to dense local data. Local depositional directions as obtained from dipmeter data. 3D seismic data, and geological interpretation represent critical information. An improved pixel-based geostatistical algorithm is proposed to account for such directional information. Case studies demonstrate the potential and limitations of the algorithm.