AN ULTRASONIC PROFILING INVESTIGATION ON SOME FRESH AND WEATHERED GRANITES OF HYDERABAD,INDIA*

The ultrasonic profiling method of measuring the compressional and shear wave velocities in cylindrical rock samples is extended to measurements in some weathered and fresh granite blocks collected from the Hyderabad (India) region. This possibility of the method provides a means of investigating the elastic properties of the less compact rocks, of which the near-surface formations are particularly important. In this article the important parts of the ultrasonic profiling instrument developed are described and the relevant aspects of the seismic wave fields and identification of the individual waves in the wavetrain responses to longitudinal excitation are considered. Compressional, shear and surface (Rayleigh) wave velocities in some fresh and weathered granites are detailed. The compressional velocities range from 4.8 km/s to 5.5 km/s in fresh granites and lie between 1.1 km/s and 2.5 km/s in weathered granites. Young's modulus and Poisson's ratios computed from the measured velocities are also presented. An empirical relation of the form log E= 4.27 + 2.11 log V_p between Young's modulus E and compressional velocities V_p in the fresh granites studied is deduced. The versatility of the approach is thus demonstrated.     

Application of the virtual refraction to near‐surface characterization at the Boise Hydrogeophysical Research Site‡

Seismic interferometry is a relatively new technique to estimate the Green's function between receivers. Spurious energy, not part of the true Green's function, is produced because assumptions are commonly violated when applying seismic interferometry to field data. Instead of attempting to suppress all spurious energy, we show how spurious energy associated with refractions contains information about the subsurface in field data collected at the Boise Hydrogeophysical Research Site. By forming a virtual shot record we suppress uncorrelated noise and produce a virtual refraction that intercepts zero offset at zero time. These two features make the virtual refraction easy to pick, providing an estimate of refractor velocity. To obtain the physical parameters of the layer above the refractor we analyse the cross‐correlation of wavefields recorded at two receivers for all sources. A stationary‐phase point associated with the correlation between the reflected wave and refracted wave from the interface identifies the critical offset. By combining information from the virtual shot record, the correlation gather and the real shot record we determine the seismic velocities of the unsaturated and saturated sands, as well as the variable relative depth to the water‐table. Finally, we discuss how this method can be extended to more complex geologic models.     

Point source radiation in inhomogeneous anisotropic structures

The ray formulae for the radiation from point sources in unbounded inhomogeneous isotropic as well as anisotropic media consist of two factors. The first one depends fully on the type and orientation of the source and on the parameters of the medium at the source. We call this factor the directivity function. The second factor depends on the parameters of the medium surrounding the source and this factor is the well-known geometrical spreading. The displacement vector and the radiation pattern defined as a modulus of the amplitude of the displacement vector measured on a unit sphere around the source are both proportional to the ratio of the directivity function and the geometrical spreading.For several reasons it is desirable to separate the two mentioned factors. For example, there are methods in exploration seismics, which separate the effects of the geometrical spreading from the observed wave field (so-called true amplitude concept) and thus require the proposed separation. The separation also has an important impact on computer time savings in modeling seismic wave fields generated by point sources by the ray method. For a given position in a given model, it is sufficient to calculate the geometrical spreading only once. A multitude of various types of point sources with a different orientation can then be calculated at negligible additional cost.In numerical examples we show the effects of anisotropy on the geometrical spreading, the directivity and the radiation pattern. Ray synthetic seismograms due to a point source positioned in an anisotropic medium are also presented and compared with seismograms for an isotropic medium.     

流体分布对松散介质中P波速度和衰减的影响

为了深入研究流体对岩石中弹性波速度和衰减的影响,必须考虑到流体的分布和粘性。引入气体包裹体模型来研究粘性流体的分布对松散介质中P波速度和衰减的影响,用气泡平均半径来描述流体分布的不均匀性,计算了不同气泡半径和频率下P波速度和衰减随饱和度变化的曲线,并与有效流体模型作了比较,由于流体喷流的存在会使Gassmann方程在高频下不适用,用干燥和饱和流体的P波、S波速度修正了理论曲线。测量了玻璃微珠中不同水饱和度下高频P波的速度和衰减,并尝试用峰值频率来计算衰减。此方法求出的Q和频谱比法求出的Q在干燥或饱和水时基本相同,随饱和度的变化规律也基本一致,但衰减峰的大小有差异。根据实测值得似事经修正的波速和衰减理论曲线从而估算出气泡平均半径,认为P波速度和衰减不仅与饱和度有关而且也与介质内部气体－液体压力平衡有关。     

Joint 3D processing of active and passive seismic data

Passive seismic provides additional illumination sources in producing reservoirs, improving the Earth's imaging obtained by standard 3D seismic surveys. The joint tomographic inversion of surface and borehole data, both active and passive, even allows the delineation of thin reservoirs that cannot be resolved by reflection tomography. As an application example, we present a feasibility study for a real case of CO₂ geological storage, showing that this operation may benefit both environment and reservoir monitoring. The origin time of micro‐earthquakes due to production operations is critical for merging active and passive data. We show here that the Wadati's method is not accurate for borehole data in a layered earth model, when the ratio between P and S velocities is not constant, as occurs in most hydrocarbon reservoirs. This drawback can be solved by deploying a few receivers at the surface close to the well.     

ANISOTROPIC Q AND VELOCITY DISPERSION OF FINELY LAYERED MEDIA1

When a seismic signal propagates through a finely layered medium, there is anisotropy if the wavelengths are long enough compared to the layer thicknesses. It is well known that in this situation, the medium is equivalent to a transversely isotropic material. In addition to anisotropy, the layers may show intrinsic anelastic behaviour. Under these circumstances, the layered medium exhibits Q anisotropy and anisotropic velocity dispersion. The present work investigates the anelastic effect in the long-wavelength approximation. Backus's theory and the standard linear solid rheology are used as models to obtain the directional properties of anelasticity corresponding to the quasi-compressional mode qP, the quasi-shear mode qSV, and the pure shear mode SH, respectively. The medium is described by a complex and frequency-dependent stiffness matrix. The complex and phase velocities for homogeneous viscoelastic waves are calculated from the Christoffel equation, while the wave-fronts (energy velocities) and quality factor surfaces are obtained from energy considerations by invoking Poynting's theorem. We consider two-constituent stationary layered media, and study the wave characteristics for different material compositions and proportions. Analyses on sequences of sandstone-limestone and shale-limestone with different degrees of anisotropy indicate that the quality factors of the shear modes are more anisotropic than the corresponding phase velocities, cusps of the qSV mode are more pronounced for low frequencies and midrange proportions, and in general, attenuation is higher in the direction perpendicular to layering or close to it, provided that the material with lower velocity is the more dissipative. A numerical simulation experiment verifies the attenuation properties of finely layered media through comparison of elastic and anelastic snapshots.     

MODELLING OF SOFT SEDIMENTS AND LIQUID-SOLID INTERFACES: MODIFIED WAVENUMBER SUMMATION METHOD AND APPLICATION1

Alekseev and Mikhailenko have developed a wavenumber-summation method which combines a finite integral transformation with a finite-difference calculation and involves no approximations other than numerical ones. However, numerical anisotropy causes velocity errors for shear waves which are unacceptable if Poisson's ratios are larger than 0.4 and unless the number of grid points per wavelength is chosen considerably higher than the value generally regarded as sufficient in finite-difference computations. To overcome this limitation in the applicability of the otherwise very powerful modelling scheme, the method is applied to the elastodynamic equations for the velocity vector. Thus, instead of solving a second-order hyperbolic system as in the case of the wave equation, solutions to a first-order hyperbolic system are computed. The finite-difference iteration is performed in a staggered grid. In addition to mastering the numerical difficulties in cases where the Poisson's ratio is unusually high, this approach results in a code which can be used for the modelling of liquid layers. With the new scheme, water reverberations are investigated in terms of normal modes. It is found that for realistic sea-bottom velocities the critical and supercritical cases exist only for P-waves. It means that compressional waves are trapped within the water layer but energy leaks into the substratum through converted shear waves. These leaky compressional normal modes attain properties similar to those of shear normal modes or Pseudo-Love waves. Due to their origin from conversion of dispersed multi-modal compressional waves the shear waves generated at the sea-bottom form a long complex wavetrain. They were found to mask the reflections from the target horizon in an offset-VSP field section.     

The effect of clay distribution on the elastic properties of sandstones

The shape and location of clay within sandstones have a large impact on the P‐wave and S‐wave velocities of the rock. They also have a large effect on reservoir properties and the interpretation of those properties from seismic data and well logs. Numerical models of different distributions of clay – structural, laminar and dispersed clay – can lead to an understanding of these effects. Clay which is located between quartz grains, structural clay, will reduce the P‐wave and S‐wave velocities of the rock. If the clay particles become aligned or form layers, the velocities perpendicular to the alignment will be reduced further. S‐wave velocities decrease more rapidly than P‐wave velocities with increasing clay content, and therefore Poisson's ratios will increase as the velocities decrease. These effects are more pronounced for compacted sandstones. Small amounts of clay that are located in the pore space will have little effect on the P‐wave velocity due to the competing influence of the density effect and pore‐fluid stiffening. The S‐wave velocity will decrease due to the density effect and thus the Poisson's ratio will increase. When there is sufficient clay to bridge the gaps between the quartz grains, P‐wave and S‐wave velocities rise rapidly and the Poisson's ratios decrease. These effects are more pronounced for under‐compacted sandstones. These general results are only slightly modified when the intrinsic anisotropy of the clay material is taken into account. Numerical models indicate that there is a strong, nearly linear relationship between P‐wave and S‐wave velocity which is almost independent of clay distribution. S‐wave velocities can be predicted reasonably accurately from P‐wave velocities based on empirical relationships. However, this does not provide any connection between the elastic and petrophysical properties of the rocks. Numerical modelling offers this connection but requires the inclusion of clay distribution and anisotropy to provide a model that is consistent with both the elastic and petrophysical properties. If clay distribution is ignored, predicting porosities from P‐wave or S‐wave data, for example, can result in large errors. Estimation of the clay distribution from P‐wave and S‐wave velocities requires good estimates of the porosity and clay volume and verification from petrographic analyses of core or cuttings. For a real data example, numerical models of the elastic properties suggest the predominance of dispersed clay in a fluvial sand from matching P‐wave and S‐wave velocity well log data using log‐based estimates of the clay volume and porosity. This is consistent with an interpretation of other log data.     

Convection in a rotating,laterally heated annulus pattern velocities and amplitude oscillations

Abstract

The velocities of the wave patterns relative to the rotating annulus have been measured with either increasing or decreasing positive radial temperature gradients and different rotation rates, with the fluid in thermal equilibrium and in contact with a rigid lid. The pattern velocities are dependent on initial conditions except in the unique areas of the stability diagram, where the velocities observed with either increasing or decreasing ΔT, overlap. The pattern velocities change discontinuously with each wave number transition, with a particularly large discontinuity at the transition from two to one wave. The frequency of the amplitude oscillations of the waves has been measured also. It has been found that the period of the oscillation of the three wave pattern is inversely proportional to the period of the pattern velocity, which means that in this case the ratio of the frequency of amplitude oscillation and the frequency of the pattern revolution is incommensurate.     

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.     

论组合法滤波的数学物理概念

本文以傅里叶变换为基础,用类似数学滤波中的数学分析方法,比较系统地研究了组合法滤波的数学物理概念及方法的内在规律。可以把方向特性函数P（ω）看作是一种周期“波动”函数,而组合的加权分布函数h（x）则是其离散“频谱”。因此可以引用符合于频谱分析的物理概念来研究反演问题。可以看出,组合滤波的加权分布函数及方向特性分别相当于数学滤波中的滤波因子及频率响应;而数字滤波的移门法等又可以引伸到构成组合法带通滤波中去。二者在数学上有其共性,可以互相引用,在物理上又有各自的特殊性。研究结果给组合法以比较全面的认识,扩大了组合法应用的功能。     

论组合法滤波的数学物理概念

本文以傅里叶变换为基础,用类似数学滤波中的数学分析方法,比较系统地研究了组合法滤波的数学物理概念及方法的内在规律。可以把方向特性函数P（ω）看作是一种周期“波动”函数,而组合的加权分布函数h（x）则是其离散“频谱”。因此可以引用符合于频谱分析的物理概念来研究反演问题。可以看出,组合滤波的加权分布函数及方向特性分别相当于数学滤波中的滤波因子及频率响应;而数字滤波的移门法等又可以引伸到构成组合法带通滤波中去。二者在数学上有其共性,可以互相引用,在物理上又有各自的特殊性。研究结果给组合法以比较全面的认识,扩大了组合法应用的功能。     

水饱和裂纹对地壳岩样中地震波速及各向异性的影响

选择４种地壳岩石样品,经干燥或水饱和处理后在不同围压条件下测量了在其中传播的纵、横波的速度及其各向异性．在大气压条件下低孔隙度(＜１％岩样中,水饱和样品中的纵波速度明显地比干燥样品中的高,但横波速度的差别不大．因为在低孔隙度岩样中纵波速度对孔隙流体的反应比横波速度敏感,可以用泊松比的变化来反映随着围压的增加晶粒间流体对弹性波传播特性的影响．根据实验数据,按Ｏ’Ｃｏｎｎｅｌｌ模型分别计算了干燥和水饱和岩样中的裂纹密度,与通过实测体应变曲线得到的裂纹孔隙度十分吻合．利用横波的速度和偏振特性可以推断岩样中定向排列微裂纹的空间取向情况．研究表明,同时测量在岩样中传播的纵、横波的速度,通过Ｖｐ／Ｖｓ比值可以给出有关颗粒边界流体的证据,也可以估计岩样中的裂纹密度．     

勒夫波群速度频散与太平洋地壳及上地幔三维构造

应用现代适配滤波频时分析技术对SRO/GDSN长周期面波记录进行处理,获得穿过太平洋盆地的117条波路径的勒夫波频散数据。我们用改进的分格频散反演方法,从混合路径频散数据中提取出对应于构成整个太平洋盆地的10°×10°分格的纯路径频散数据。所获得的18s到215s范围内的28个周期的太平洋勒夫波群速度的横向分布,揭示了速度随海底年龄的增大而逐渐增加的总趋势,以及其它变化细节。例如,对于同一等年龄线区而言,北太平洋及南太平洋的群速度要比中太平洋的低0.1-0.3km/s。在分格模型反演的基础上得出的深至300km的三维剪切波速度结构提供了太平洋盆地地壳及上地幔横向不均匀性的详尽的描述。     

Indirect Boundary Element Method applied to fluid–solid interfaces

In this paper scattering of elastic waves in fluid–solid interfaces is investigated. We use the Indirect Boundary Element Method to study this wave propagation phenomenon in 2D models. Three models are analyzed: a first one with an interface between two half-spaces, one fluid on the top part and the other solid in the bottom; a second model including a fluid half-space above a layered solid; and finally, a third model with a fluid layer bounded by two solid half-spaces. The source, represented by Hankel's function of the second kind, is always applied in the fluid. This indirect formulation can give to the analyst a deep physical insight on the generated diffracted waves because it is closer to the physical reality and can be regarded as a realization of Huygens' principle. In any event, mathematically it is fully equivalent to the classical Somigliana's representation theorem. In order to gauge accuracy we test our method by comparing with an analytical solution known as Discrete Wave Number. A near interface pulse generates scattered waves that can be registered by receivers located in the fluid and it is possible to infer wave velocities of solids. Results are presented in both time and frequency domain, where several aspects related to the different wave types that emerge from this kind of problems are pointed out.     

Effects of pore aspect ratios on velocity prediction from well‐log data

We develop a semi‐empirical model which combines the theoretical model of Xu and White and the empirical formula of Han, Nur and Morgan in sand–clay environments. This new model may be used for petrophysical interpretation of P‐ and S‐wave velocities. In particular, we are able to obtain an independent estimation of aspect ratios based on log data and seismic velocity, and also the relationship between velocities and other reservoir parameters (e.g. porosity and clay content), thus providing a prediction of shear‐wave velocity. To achieve this, we first use Kuster and Toksöz's theory to derive bulk and shear moduli in a sand–clay mixture. Secondly, Xu and White's model is combined with an artificial neural network to invert the depth‐dependent variation of pore aspect ratios. Finally these aspect ratio results are linked to the empirical formula of Han, Nur and Morgan, using a multiple regression algorithm for petrophysical interpretation. Tests on field data from a North Sea reservoir show that this semi‐empirical model provides simple but satisfactory results for the prediction of shear‐wave velocities and the estimation of reservoir parameters.     

土层的性质与剪切波速的关系

本文用多元逐步回归分析法和非线性模型参数确定法,对西安、徐州、鞍山等地的土层埋深及物理水理指标与剪切波速的关系进行了综合的分析和对比,结果显示深度和液性指数对剪切波速有明显的影响。作者认为,用非线性模型参数确定法有利于提高V_(?)经验公式的精度。     

STATIC CORRECTIONS FOR SHEAR WAVE SECTIONS*

The computation of static corrections requires information about subsurface velocities. This information can be obtained by different methods: surface wave analysis, short refraction lines, downhole times, uphole times and first arrivals from seismograms. For pure shear waves generated by SH sources the analysis of first arrivals from seismograms combined, if necessary, with short refraction lines has proved to be most accurate and economic. A comparison of first-arrival plots from P- and S-wave surveys of the same line measured in areas of unconsolidated sediments in northern Germany illustrates the characteristic differences between the two velocity models. P-waves show a marked velocity increase at the water table from about 600 to 1800 m/s. S-wave velocities of the same strata increase gradually from about 100 to 400 m/s. As a consequence, S-wave models are vertically and laterally more complex and, in general, show no significant velocity increase at a defined boundary as P-wave models do. Therefore, other suitable correction levels with specific velocities must be chosen. A comparison of “tgd-corrections” (correction time between geophone position and datum level) for P- and S-waves in areas of unconsolidated sediments shows that their ratio is different from the P-/S-velocity ratio for the respective correction level because of the greater depth of the S-wave refractor. Therefore, P- and S-waves are influenced by different near-surface anomalies, and time corrections calculated for both wave types are largely independent.     

ANALYSIS OF THE COMMON MIDPOINT GATHER BY DECOMPOSITION INTO ELEMENTARY WAVEFRONTS1

The imaginary wavefront of a common midpoint gather of seismic data is analysed in terms of elementary wavefronts. The wavefronts emitted by the sources are split into a sequence of elementary circular wavefronts. Each such wavefront is independently transmitted through the subsurface and received by a group of receivers associated with the sources. The theory is first demonstrated in the case of a one-layer medium and is then extended to the case of a multilayer medium. We show that pairs of imaginary stationary points F and F' occur in the first medium through which all the rays, corresponding to the group of sources and receivers associated with the elementary wavefront pass. The location of the points is determined by 2n equations with 2n unknowns; (2n - 1) of these equations correspond to n layers of the medium, and an additional equation is given by the geometry of the CMP gather. Applying Fermat's principle between the points F and F', we obtain an algorithm for CMP ray tracing. The wavefront sampling problem is considered but not theoretically solved. Results obtained in the general case (a multilayer medium with lateral velocity variation) are presented and compared with those obtained by other modelling techniques.     

哀牢山变质带元江──墨江剖面岩石的纵波波速特征及其地质意义

依据研究区的地热梯度（25℃／km）,在高温高压（最高温度为1050℃,最高压力为1.2GPa）条件下系统测量了横穿红河-哀牢山断裂带的元江-墨江地质剖面上的哀牢山岩群各类变质岩（千枚岩、片岩、浅粒岩、变粒岩、大理岩和片麻岩）的纵波速度．实验结果表明,不同岩类的纵波速度随温度压力变化的趋势不同．在相当于衷牢山岩群变质岩峰期变质温度和压力条件下（P=0．4-0.8GPa,T=35-700℃）,测得大部分岩石的纵波速度为5．50－5．80km／s,这一纵波速度值与区域地球物理测深揭示的中地壳低速层的纵波速度相当因此,结合该区变质岩、地壳内热状态及地球物理测深研究成果可初步认为：组成哀牢山岩群的浅粒岩、变粒岩、酸性片麻岩以及部分千枚岩、片岩为该地区中地壳低速层的主要岩石类型．