On shear-wave splitting in the Los Angeles basin

Three-component seismograms at the three USC stations, PVP, GFP and DHB, have been examined. Most earthquakes, with magnitudes ranging from 1.4 to 5.0, within a period from 1985 to 1988, show evidence of shear-wave splitting. The preferred polarization of the first split-shear wave arrivals at PVP is nearly in N-S which is consistent with both regional maximum horizontal compressive stress direction and local subsurface fault strike, showing that shear-wave splitting is caused by liquid-filled cracks or fractures associated with the N-S faulting. The polarizations of first shear wave arrivals at GFP are roughly divided into two almost perpendicular directions, ENE-WSW and NNW-SSE, which are parallel or perpendicular to the strike of the geology or topography near the station. Because GFP is near the foothills of Santa Monica Mountains, the shear-wave arrivals may be disturbed by topographic irregularities and by subsurface dipping interfaces. Two examples at DHB clearly display shear-wave splitting. Their polarizations of shear wave are in the direction of N-S, which agree with the fragmentary surface and fracturing direction there. From some relatively reliable delay times, the crack densities at three stations are given, that is, 0.025 at PVP, 0.020 at GFP and 0.045 at DGB. No systematic change of shear-wave polarization is discovered in this study.     

Controlled sources for shear-wave surveys in mines

The ability to analyse shear-wave anisotropy in a mine environment is greatly aided by using multiple source orientations of a reproducible, impulsive shear-wave source. The analysis of what is probably the first controlled source shear-wave experiment in a mine environment demonstrates clearly that shear-wave polarizations and time delays between split shear-wave arrivals are reliably measured because of the use of multiple source orientations rather than a single shear-wave source. Reliability is further aided by modelling the shear-wave source radiation pattern, which allows for the unequivocal discrimination between seismic raypaths where shear-wave splitting did and did not occur. The analysis also demonstrates the great importance of high reproducibility of the seismic source for the use of shear waves in time-lapse surveys to monitor changes in a rockmass.     

SECONDARY SHEAR WAVES FROM SOURCE BOREHOLES1

Previous studies of radiation from point sources in fluid-filled boreholes have most often been based on far-field, stationary phase analysis. In these papers, the explicit contribution of the borehole itself acting as a waveguide has not been properly considered, with a few exceptions. In general, these studies accurately describe S-wave radiation in high-velocity rocks such as granites and limestones and P-wave radiation in most rocks, and experiments have confirmed this. However, tube waves directly influence the external wavefield and in fact create a shear-wave ‘wake’ outside the borehole due to constructive interference of tube-wave emission if a velocity condition is met. This constructive interference or wake is generated when the tube-wave velocity is greater than the shear-wave velocity. When this happens, a tube-wave complex pole invalidates the mathematical assumptions for stationary phase analysis and the stationary phase predictions do not agree with experimentally derived radiation patterns. Shales at shallow depths and other soft sediments characteristically have tube-wave velocities greater than shear-wave velocities. Because the tube-wave is of relatively high amplitude compared to body waves generated directly by the source, these secondary shear waves can be the highest amplitude arrivals on receiver arrays. The shape and properties of these secondary shear waves are calculated and shown to have identical properties to Mach waves of aerodynamics and seismology. For instance, these waves are geometrically conical and the aperture of the cone and the moveout velocity can be calculated. This paper also demonstrates the important effect that casing has on the Mach waves and provides predictions about when these waves are likely to be observed. Finally, evidence of Mach waves in data sets is examined and it is shown how these waves have been confused with receiver borehole tube waves. It is possible, though rare, that the tube-wave velocity of the borehole is greater than the compressional-wave velocity of the surrounding medium. In this case secondary compressional or compressional Mach waves would be generated although this problem is not addressed here.     

Key elements of regional seismic velocity models for long period ground motion simulations

Regional 3-D seismic velocity models used for broadband strong motion simulations must include compressional-wave velocity (Vp), shear-wave velocity (Vs), intrinsic attenuation (Qp, Qs), and density. Vs and Qs are the most important of these parameters because the strongest ground motions are generated chiefly by shear- and surface-wave arrivals. Because Vp data are more common than Vs data, many researchers first develop a Vp model and convert it to a Vs model. I describe recent empirical relations between Vs, Vp, Qs, Qp, and density that allow velocity models to be rapidly and accurately calculated. The U. S. government reserves the right to retain a non-exclusive, royalty free license in and to any copyright.     

Approximate traveltime inversion in downhole microseismic monitoring

In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media.     

Crustal seismic anisotropy in southeastern Capital area, China

Shear-wave splitting parameters of 24 stations in southeastern Capital area of North China (38.5°N~39.85°N, 115.5°E ~118.5°E) are obtained with systematic analysis method of shear-wave splitting (SAM) based on the data recorded by Capital Area Seismograph Network (CASN) from 2002 to 2005. The results show that the average polarization of fast shear-wave in southeastern Capital area is consistent with regional maximum horizontal prin- cipal compressive stress in the area, and is also consistent with maximum horizontal principal compressive strain from GPS in North China. The average shear-wave splitting in southeastern Capital area (in basin) is different from that in northwestern Capital area where uplifts and basin exist, which means that tectonics can be related to shear-wave splitting results. Research also shows that the distribution of faults around stations can obviously affect the shear-wave splitting results, and complicated distribution of faults can result in much more scatter of shear-wave splitting. Moreover, in the north and south of the studied area (southeastern Capital area), the polariza- tions of fast shear-wave are not very consistent, which may be related to differences in tectonic and stress for the two areas.     

Application research of seismic method in assessment of active fault

ＡｐｐｌｉｃａｔｉｏｎｒｅｓｅａｒｃｈｏｆｓｅｉｓｍｉｃｍｅｔｈｏｄｉｎａｓｅｓｍｅｎｔｏｆａｃｔｉｖｅｆａｕｌｔＭＩＮＧＣＡＩＸＵ（徐明才）ＪＩＮＧＨＵＡＧＡＯ（高景华）ＭＩＮＧＴＡＯＣＨＡＩ（柴铭涛）ＧＵＡＮＧＫＥＷＡＮＧ（王广科）ＪＩ...     

Fractured reservoir delineation using multicomponent seismic data

The characteristic seismic response to an aligned-fracture system is shear-wave splitting, where the polarizations, time-delays and amplitudes of the split shear waves are related to the orientation and intensity of the fracture system. This offers the possibility of delineating fractured reservoirs and optimizing the development of the reservoirs using shear-wave data. However, such applications require carefully controlled amplitude processing to recover properly and preserve the reflections from the target zone. Here, an approach to this problem is suggested and is illustrated with field data. The proposed amplitude processing sequence contains a combination of conventional and specific shear-wave processing procedures. Assuming a four-component recording (two orthogonal horizontal sources recorded by two orthogonal horizontal receivers), the split shear waves can be simulated by an effective eigensystem, and a linear-transform technique (LTT) can be used to separate the recorded vector wavefield into two principal scalar wavefields representing the fast and slow split shear waves. Conventional scalar processing methods, designed for processing P-waves, including noise reduction and stacking procedures may be adapted to process the separated scalar wavefields. An overburden operator is then derived from and applied to the post-stacked scalar wavefields. A four-component seismic survey with three horizontal wells drilled nearby was selected to illustrate the processing sequence. The field data show that vector wavefield decomposition and overburden correction are essential for recovering the reflection amplitude information in the target zone. The variations in oil production in the three horizontal wells can be correlated with the variations in shear-wave time-delays and amplitudes, and with the variations in the azimuth angle between the horizontal well and the shear-wave polarization. Dim spots in amplitude variations can be correlated with local fracture swarms encountered by the horizontal wells. This reveals the potential of shear waves for fractured reservoir delineation.     

剪切波速和含水率双指标评价黄土震陷

为了研究地震荷载作用下黄土场地的震陷预测问题,以兰州黄土为研究对象来分析剪切波速和含水率与震陷系数的关系,并分别得出拟合函数关系。通过分析剪切波速和含水率在震陷过程中的内在联系,建立剪切波速和含水率双指标预测黄土震陷系数的方法,并依据该方法预测永登地震区黄土场地震陷量。研究结果表明:在一定的剪切波速区间内,震陷系数随着含水率的增加而增大;在一定的含水率区间内,震陷系数随着剪切波速的增加而减小。同时提出剪切波速可以反映黄土不因外部条件而变化的潜在震陷能力为绝对震陷能力的概念,相应的震陷系数称为绝对震陷系数,含水率影响潜在震陷能力的释放。通过震害实例计算验证剪切波速和含水率双指标预测方法的可行性。     

Using refracted shear waves for velocity estimation

The most difficult part of multicomponent processing is the estimation of the shear-wave velocity map for migration. We used refracted shear waves and a simple iterative method called wavefield continuation (WFC) to evaluate the shallow shear-wave velocity profile on a real data example. The WFC was developed in 1981 by Clayton and McMechan to determine compressional-wave velocity profiles from refracted compressional waves. The application to refracted shear waves is straightforward. The real data example shows that shear structure can be easily determined independently of the compressional structure.     

卢龙地区S波偏振与上地壳裂隙各向异性

由三分量数字地震仪组成的小孔径流动台网记录了1982年10月19日河北卢龙M_s=6.1级地震的部分余震.用质点运动图的方法对横波的偏振进行了分析。研究结果表明,在横波窗内的各观测点都存在横波的分裂现象.不同离源角和方位角快波偏振的水平投影都具有近NE40°方向的优势取向,与根据卢龙地震两组断层错动在各向同性介质中所辐射的横波的偏振方向不一致.这可以由传播介质中应力所导致裂隙的定向排列来解释.这一观测结果提供了卢龙地区脆性上地壳大范围膨胀各向异性（EDA）的证据,并表明这一地区直立平行排列裂隙取向和水平主压应力的方向为NE40°.     

卢龙地区S波偏振与上地壳裂隙各向异性

由三分量数字地震仪组成的小孔径流动台网记录了1982年10月19日河北卢龙Ms=6.1级地震的部分余震.用质点运动图的方法对横波的偏振进行了分析。研究结果表明,在横波窗内的各观测点都存在横波的分裂现象.不同离源角和方位角快波偏振的水平投影都具有近NE40°方向的优势取向,与根据卢龙地震两组断层错动在各向同性介质中所辐射的横波的偏振方向不一致.这可以由传播介质中应力所导致裂隙的定向排列来解释.这一观测结果提供了卢龙地区脆性上地壳大范围膨胀各向异性（EDA）的证据,并表明这一地区直立平行排列裂隙取向和水平主压应力的方向为NE40°.     

基于多重互相关函数计算地探信号到时差

论文介绍了一种提高地震等地探信号到时差精度的新方法-多重互相关函数法. 从理论上推导出了多重互相关函数的公式,且通过实验验证了用此函数不仅能提高信号的信噪比,还能提取出相似信号所共有的频率信号,进而精确地识别出它们的到时差.探讨了此原理在场地、地震剪切波速信号到时差提取中的应用,分析结果显示:相比互相关函数,多重互相关函数能更准确地求取两相似信号的初至时刻差.     

剪切波分裂分析系统SAM(2007)——区域地震台网资料应用实例

通过对云南、辽宁、福建遥测地震台网的波形资料分析实例,介绍了"剪切波分裂系统分析方法"软件SAM(2007)的用法。通过对区域台网资料的处理分析,表明SAM(2007)是研究剪切波分裂的一个比较实用的软件,能够有效的处理快、慢剪切波识别的问题,研究地壳介质的地震各向异性问题,能普遍用于分析国内区域地震台网资料。     

基于约束条件横波速度反演和流体替代

依据辽东湾地区一组在地层条件和不同流体相态下（气饱和与水饱和）实验室测定的岩石纵波速度、横波速度及密度等数据,讨论Gassmann方程在横波速度反演和流体替代方面的应用效果,指出岩石骨架剪切模量和体积模量估算精度是影响纵波、横波速度预测精度的关键因素之一,提出基于约束条件横波速度反演和流体替代方法,速度预测精度的改善十分明显.     

Interpreting non-orthogonal split shear waves for seismic anisotropy in multicomponent VSPS

There are two main sources of non-orthogonality in multicomponent shear-wave seismics: inherent non-orthogonal split shear waves arising from substantial ray deviation in off-symmetry planes due to strong anisotropy or complex overburden, and apparent non-orthogonal split shear waves in the horizontal plane due to variation of the angle of incidence even if the two shear waves along the raypath are orthogonal. Many techniques for processing shear-wave splitting in VSP data ignore these kinds of non-orthogonality of the split shear waves. Assuming inherent non-orthogonality in zero-offset VSPs, and apparent non-orthogonality in offset VSPs, we derive equations for the four-component data matrix. These can be solved by extending the linear-transform technique (LTT) to determine the shear-wave polarizations in zero-offset and offset VSPs. Both full-wave synthetic and field data are used to evaluate the technique and to examine the effects of non-orthogonal polarized split shear waves. If orthogonality is incorrectly assumed, errors in polarization measurements increase with the degree of non-orthogonality, which introduces a consistent decreasing trend in the polarization measurements. However, the effect of non-orthogonality on the estimation of geophone orientation and time delays of the two split shear waves is small and negligible in most realistic cases. Furthermore, for most cases of weak anisotropy (less than 5% shear-wave anisotropy) apparent non-orthogonality is more significant than inherent non-orthogonality. Nevertheless, for strong anisotropy (more than 10% shear-wave anisotropy) with complicated structure (tilted or inclined symmetry axis), inherent non-orthogonality may no longer be negligible. Applications to both synthetic and real data show that the extended linear-transform techniques permit accurate recovery of polarization measurements in the presence of both significant inherent and apparent non-orthogonality where orthogonal techniques often fail.     

Source Directivity, Signal Decorrelation, Spectral Modulation and Analysis of Spatio-temporal Patterns of Multiple Explosions

—?Spectral modulation is used as a criterion for discriminating multiple mining explosions from earthquakes and single blasts. Features in the spectra of regional arrivals from quarry blasts may be explained by the spatio-temporal configuration of shot patterns and are explored via a model based on spatial waveform decorrelation and propagation delay (directivity) effects. The phenomenon of decreasing modulation with decreasing average phase velocities of the seismic wave arrivals can be attributed to a trend of decreasing similarities of waveforms with increasing distance between individual explosions making up the shot pattern for the respective arrivals. Rescaling of the modulation patterns with respect to frequency can be explained by the Doppler effect.¶The postulated relationship between shot pattern layouts and sizes, together with their firing sequences, on one hand, and the similarities and differences in modulation characteristics among the various regional arrivals, on the other hand, can be used to derive information about the mining practices in a region that is otherwise not accessible. Moreover, such analyses could be utilized for detecting unusual activities in the CTBT monitoring system.     

ANISOTROPY DETECTED IN SHALLOW CLAYS USING SHEAR-WAVE SPLITTING IN A VSP SURVEY1

Shear waves can provide valuable information about seismic anisotropy. On entering an anisotropic medium, a shear wave generally splits (shear-wave splitting) into a fast and a slow quasi-shear wave with polarizations fixed by the elastic properties of the medium and direction of travel. If the medium contains planar discontinuities with common normals, the fast shear wave will be suitably propagated if its polarization lies in the plane of the discontinuities. Measuring this polarization, using a VSP geometry with oriented three-component geophones in the borehole, offers the possibility of monitoring the orientation and density of the discontinuities as a function of depth. Such a shear-wave VSP was carried out in an uncased 0.3 m diameter borehole drilled to a depth of 120 m in the north of The Netherlands. The upper 80 m of the sequence, consisting of a glacial till and sands and clays of Pleistocene age, was studied. The clays in this sequence have been subjected to glacial deformation and as a result are overconsolidated and locally fissured. In our shallow VSP experiment, shear-wave splitting and therefore anisotropy was identified at various geophone depths for one source offset. Hodograms showed a consistent polarization of the fast shear-wave component over a large depth interval. Under the assumption that the anisotropy was caused by planar discontinuities with common normals, this polarization direction gives the strike of the fissures in this interval. The polarization direction of the fast S-wave did not correspond exactly with the strike which was obtained from geological information on the fissures. The geological information was from undisturbed oriented 70 mm core samples taken at 3 m intervals in the borehole. The discrepancy, however, could be explained in terms of dipping fissures, and such a dip was confirmed by the geological and geotechnical information. The orientation of fissures is an important factor in the directional deformation and strength characteristics of clays as far as geotechnical behaviour is concerned. This study thus illustrates a practical application of shear-wave splitting observed in shallow shear-wave VSP for geology and geotechnical engineering.     

EXAMINATION OF A SPECTRAL METHOD FOR MEASURING THE EFFECTS OF ANISOTROPY1

The most diagnostic effect of anisotropy on shear waves is shear-wave splitting. This phenomenon creates a distinctive signature in the 3D particle motion. Methods to extract the effects of anisotropy from shear-wave data attempt to measure details of this motion. Many techniques have been published recently which process the shear waves in the time or frequency domain. Here we examine the way in which information on the interference effects between the split shear waves is contained within the frequency domain, and suggest some criteria which may be used in future processing algorithms. The time-delay between the split shear waves, and the polarization direction of the leading shear wave can be converted into easily measured features from analysis of the Fourier spectrum of the shear-wave signal on each component of motion. These features arise in the spectral interference patterns which are formed by the interaction between the two closely-spaced and similar waveforms. The interference patterns are interpreted for synthetic and observed seismogram data.     

测井方法测定含煤地层的岩石弹性参数

本文论述了利用声波全波列数字测井资料提取横波的两种方法:子波逼近法和相似相关法。根据测井获得的纵波、横波和密度资料,计算了岩石的弹性模量,所得结果有足够的精度,为煤矿生产和建设提供了有用的岩石物理力学参数。文中还讨论了用静态法和动态法所测得的弹性模量存在的某些差别。