SEISMIC AND ELECTRICAL PROPERTIES OF UNCONSOLIDATED PERMAFROST1

A model has been developed to relate the velocities of acoustic waves V_p and V_s in unconsolidated permafrost to the porosity and extent of freezing of the interstitial water. The permafrost is idealized as an assemblage of spherical quartz grains embedded in a matrix composed of spherical inclusions of water in ice. The wave-scattering theory of Kuster and Toksoz is used to determine the effective elastic moduli, and hence the acoustic velocities. The model predicts V_p and V_s to be decreasing functions of both the porosity and the water-to-ice ratio. The theory has been applied to laboratory measurements of V_p and V_s in 31 permafrost samples from the North American Arctic. Although no direct measurements were made of the extent of freezing in these samples, the data are consistent with the predictions of the model. Electrical resistivity measurements on the permafrost samples have demonstrated their essentially resistive behaviour. The ratio of resistivity of permafrost in its frozen state to that in its unfrozen state has been related to the extent of freezing in the samples. Electromagnetic and seismic reflection surveys can be used together in areas of permafrost: firstly an EM survey to determine the extent of freezing and then the acoustic velocity model to predict the velocities in the permafrost. The necessary transit time corrections can thus be made on seismic reflection records to compensate for the presence of permafrost.     

Determination of shale stiffness tensor from standard logs

A technique allowing inversion of the shale stiffness tensor from standard logging data: sonic velocities, density, porosity and clay content is developed. The inversion is based on the effective medium theory. The testing of the technique on laboratory measurements of the elastic wave velocities in shale samples shows that the inversion makes it possible to predict the elastic wave velocities V_P, V_S1 and V_S2 in any direction within an error of a few per cent. The technique has been applied for the stiffness tensor inversion along a well penetrating a shale formation of the Mississippian age altered by thin layers of limestone. It is demonstrated that the symmetry of a stiffness tensor inverted at the sonic frequency (2 kHz) is slightly orthorhombic and taking into account the experimental errors, can be related to the vertical transverse isotropy symmetry. For the productive interval of the shale formation, the Thomsen parameters ?, γ, and δ average, respectively, 0.32, 0.25 and 0.21, which indicate anelliptic behaviour of the velocities in this shale. The coefficients of anisotropy of this shale interval are around 24% and 20% for the compressional and shear waves, respectively. The values of the inverted velocities in the bedding plane for this interval are in good agreement with the laboratory measurements. The technique also allows inversion of the water saturation of the formation (Sw) and the inverted values are in agreement with the Sw values available for this formation. A Backus‐like upscaling of the inverted stiffness tensors is carried out for the lower and upper bounds of the frequency band used in the crosswell tomography (100 Hz and 500 Hz). These results can serve as an initial velocity model for the microearthquake location during hydrofracking of the shale formation.     

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.     

Comparing shear-wave velocity profiles inverted from multichannel surface wave with borehole measurements

Recent field tests illustrate the accuracy and consistency of calculating near-surface shear (S)-wave velocities using multichannel analysis of surface waves (MASW). S-wave velocity profiles (S-wave velocity vs. depth) derived from MASW compared favorably to direct borehole measurements at sites in Kansas, British Columbia, and Wyoming. Effects of changing the total number of recording channels, sampling interval, source offset, and receiver spacing on the inverted S-wave velocity were studied at a test site in Lawrence, Kansas. On the average, the difference between MASW calculated V_s and borehole measured V_s in eight wells along the Fraser River in Vancouver, Canada was less than 15%. One of the eight wells was a blind test well with the calculated overall difference between MASW and borehole measurements less than 9%. No systematic differences were observed in derived V_s values from any of the eight test sites. Surface wave analysis performed on surface data from Wyoming provided S-wave velocities in near-surface materials. Velocity profiles from MASW were confirmed by measurements based on suspension log analysis.     

In situ measurements of sound velocities and densities across the orthopyroxene → high-pressure clinopyroxene transition in MgSiO3 at high pressure

Using acoustic measurement interfaced with a large volume multi-anvil apparatus in conjunction with in situ X-radiation techniques, we are able to measure the density and elastic wave velocities (V_P and V_S) for both ortho- and high-pressure clino-MgSiO₃ polymorphs in the same experimental run. The elastic bulk and shear moduli of the unquenchable high-pressure clinoenstatite phase were measured within its stability field for the first time. The measured density contrast associated with the phase transition OEN → HP-CEN is 2.6-2.9% in the pressure of 7-9 GPa, and the corresponding velocity jumps are 3-4% for P waves and 5-6% for S waves. The elastic moduli of the HP-CEN phase are K_S=156.7(8) GPa, G = 98.5(4) GPa and their pressure derivatives are K_S′=5.5(3) and G′ = 1.5(1) at a pressure of 6.5 GPa, room temperature. In addition, we observed anomalous elastic behavior in orthoenstatite at pressure above 9 GPa at room temperature. Both elastic wave velocities exhibited softening between 9 and 13-14 GPa, which we suggest is associated with a transition to a metastable phase intermediate between OEN and HP-CEN.     

Experimental verification of the fracture density and shear‐wave splitting relationship using synthetic silica cemented sandstones with a controlled fracture geometry

We present laboratory ultrasonic measurements of shear‐wave splitting from two synthetic silica cemented sandstones. The manufacturing process, which enabled silica cementation of quartz sand grains, was found to produce realistic sandstones of average porosity 29.7 ± 0.5% and average permeability 29.4 ± 11.3 mD. One sample was made with a regular distribution of aligned, penny‐shaped voids to simulate meso‐scale fractures in reservoir rocks, while the other was left blank. Ultrasonic shear waves were measured with a propagation direction of 90° to the coincident bedding plane and fracture normal. In the water saturated blank sample, shear‐wave splitting, the percentage velocity difference between the fast and slow shear waves, of <0.5% was measured due to the bedding planes (or layering) introduced during sample preparation. In the fractured sample, shear‐wave splitting (corrected for layering anisotropy) of 2.72 ± 0.58% for water, 2.80 ± 0.58% for air and 3.21 ± 0.58% for glycerin saturation at a net pressure of 40 MPa was measured. Analysis of X‐ray CT scan images was used to determine a fracture density of 0.0298 ± 0.077 in the fractured sample. This supports theoretical predictions that shear‐wave splitting (SWS) can be used as a good estimate for fracture density in porous rocks (i.e., SWS = 100ε_f, where ε_f is fracture density) regardless of pore fluid type, for wave propagation at 90° to the fracture normal.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

ACOUSTIC MODELING OF THE SEAFLOOR*

Data from routine seismic surveys contain considerable information about the geo-acoustic properties of the seafloor. Waves are reflected at a wide range of angles of incidence from near-vertical reflections (higher multiples) to supercritical reflections (primary and lower multiples). The reflection coefficient is approximately constant for small angles of incidence (< 10°) but varies greatly for larger angles of incidence. Near-vertical reflections are used to determine the seafloor density. The P-velocity in the seafloor is determined in advance from the critical distance using the amplitude variation of the primary as well as the multiples. The V_p/V_S ratio is determined by modeling the amplitude variation with the angle of incidence. The primary reflection from the seafloor and the first three multiples are included in the modeling. Seismic data obtained with both conventional and superlong airgun arrays have been modeled. Data collected from the Barents Sea show that even if the P-velocity is the same at different sites, the V_p/V_s ratio, density and Poisson's ratio vary significantly. The most extreme example shows that for a P-velocity of 2.80 km/s the V_p/V_s ratio varies between 1.9 and 6.0. The corresponding densities vary from 2.36 g/cm³ to 1.80 g/cm³ and the Poisson's ratio varies from 0.31 to 0.49. The acoustic modeling offers a method of assessing the mean geotechnical or mechanical properties of larger volumes of marine sediments in terms of incompressibility, shear modulus and Poisson's ratio.     

砂-粉混合料剪切波速的弯曲元及共振柱试验对比试验研究

为探讨测试方法、试验条件以及级配特征对砂-粉混合料剪切波速的影响，对具有不同细粒含量FC，相对密度D_r以及初始有效围压$\sigma_{3 \mathrm{c}}^{\prime}$的砂-粉混合料进行弯曲元和共振柱试验。结果表明:当D_r =35%或50%时，剪切波速V_s随FC的增大先减小后增大；当D_r =60%时，V_s随FC的增大而减小；弯曲元试验测得的V_s明显大于共振柱试验测得的V_s，随着FC的增大，弯曲元试验与共振试验得到的V_s差值逐渐减小，而当FC＞20%时，两种试验得到的V_s基本相同。在考虑V_s弥散性之后，不同FC的混合料弯曲元与共振柱试验得到的V_s结果具有较好的一致性。基于Hardin模型建立的砂-粉混合料V_s预测方法具有较好的预测效果。     

THE USE OF SEISMIC SHEAR WAVES AND COMPRESSIONAL WAVES FOR LITHOLOGICAL PROBLEMS OF SHALLOW SEDIMENTS*

From a great variety of in situ shear wave experiments, i.e., reflection, refraction and borehole surveys in the shallow sediments of the north German plains, several specific properties have been derived. Shear waves (S) differ from compressional waves (P) in that:

• 1 they are not affected by the degree of water saturation. Thus, they provide a better correlation between the velocity V_s and (solid) lithology;
• 2 they generally have lower frequencies, but shorter wavelength and, hence, a better resolution of thin layers;
• 3 they have lower absorption Q_s^?1 and hence a better penetration in partially saturated and gas-containing sediments than P-waves.

Correlations have been established between V_s and the confining pressure and between reduced V_s values and several lithological parameters like the grain size of sandy material. More lithological and hydrological information is obtained by using S- and P-wave surveys along the same profile. The best information on a sedimentological structure is obtained by the simultaneous observation of V_s, V_p, Q_s and Q_p.     

Effective velocities in fractured media: a numerical study using the rotated staggered finite‐difference grid

The modelling of elastic waves in fractured media with an explicit finite‐difference scheme causes instability problems on a staggered grid when the medium possesses high‐contrast discontinuities (strong heterogeneities). For the present study we apply the rotated staggered grid. Using this modified grid it is possible to simulate the propagation of elastic waves in a 2D or 3D medium containing cracks, pores or free surfaces without hard‐coded boundary conditions. Therefore it allows an efficient and precise numerical study of effective velocities in fractured structures. We model the propagation of plane waves through a set of different, randomly cracked media. In these numerical experiments we vary the wavelength of the plane waves, the crack porosity and the crack density. The synthetic results are compared with several static theories that predict the effective P‐ and S‐wave velocities in fractured materials in the long wavelength limit. For randomly distributed and randomly orientated, rectilinear, non‐intersecting, thin, dry cracks, the numerical simulations of velocities of P‐, SV‐ and SH‐waves are in excellent agreement with the results of the modified (or differential) self‐consistent theory. On the other hand for intersecting cracks, the critical crack‐density (porosity) concept must be taken into account. To describe the wave velocities in media with intersecting cracks, we propose introducing the critical crack‐density concept into the modified self‐consistent theory. Numerical simulations show that this new formulation predicts effective elastic properties accurately for such a case.     

Modeling of mineralogical composition,density and elastic wave velocities in anhydrous magmatic rocks

We use the technique of direct minimization of the Gibbs free energy of the 8-component (K₂O-Na₂O-Fe₂O₃-FeO-CaO-MgO-Al₂O₃-SiO₂) multiphase system in order to determine the equilibrium mineral assemblages of rocks of different bulk chemical compositions equilibrated at various P-T conditions. The calculated modal compositions of rocks and experimental data on elastic moduli of single crystals are then used to calculate densities and isotropic elastic wave velocities of rocks together with their pressure and temperature derivatives. Sufficient accuracy of the calculations is confirmed by comparison with experimental data on the gabbro-eclogite transformation and precise ultrasonic measurements of elastic wave velocities in a number of magmatic and metamorphic rocks.We present calculated phase diagrams with isolines of density, elastic wave velocities, and their pressure and temperature derivatives for several anhydrous magmatic rocks, from granite to lherzolite. Density and elastic properties of rocks are controlled by their chemical compositions, especially the SiO₂ content, and by P-T of equilibration, and they increase with pressure due to mineral reactions changing mineral assemblages from plagioclase-bearing and garnet-free to garnetbearing and plagioclase-free. TheV _p -density correlation is high, and shows two clear trends: one for iron-poor ultramafic rocks and another for all the other rocks considered. Mineral reactions, which occur at high pressures, changeV _p and density of anhydrous magmatic rocks following the well-known Birch (or a similar) law.Felsic, intermediate and mafic rocks can be well distinguished in theV _p -V _p /V _s - diagram, although their values ofV _p can be close to one another. TheV _p -V _p /V _s -density diagrams together with calculated phase diagrams can serve as efficient instruments for petrologic interpretation of seismic velocities.     

基于典型地质特征分析的大同盆地VS30预测研究

土体剪切波速是进行土层地震反应分析的动力学参数，对场地地震动参数确定具有重要意义。基于地质地貌分析，将大同盆地划分为5类典型地质单元。对盆地1429个钻孔剪切波速资料进行分析，探讨V_S30与V_S20的相关性，研究土体埋深、岩性、地质单元、标贯击数及密实度等地质特征对V_S的影响，并基于地质单元、剪切波速比、密实度系数及第四系上部覆盖层厚度相关性分析给出土体V_S30预测模型。研究结果表明，基于典型地质特征的V_S30预测模型拟合优度R2>0.90，预测精度很高，对于离散性较大、直接拟合估算较差及无剪切波速场地来说，以区分地质单元及土体类型的方式进行V_S30分解预测是良好的研究思路。首次在区分地质单元及土体类型的前提下提出剪切波速比及密实度系数，并将其与第四系上部覆盖层厚度综合应用于V_S30预测研究。研究结果可为大同盆地城市防震减灾规划、震害预测、区域性地震安全评价提供重要技术支撑。     

Lithology and hydrocarbon mapping from multicomponent seismic data

Elastic rock properties can be estimated from prestack seismic data using amplitude variation with offset analysis. P‐wave, S‐wave and density ‘reflectivities’, or contrasts, can be inverted from angle‐band stacks. The ‘reflectivities’ are then inverted to absolute acoustic impedance, shear impedance and density. These rock properties can be used to map reservoir parameters through all stages of field development and production. When P‐wave contrast is small, or gas clouds obscure reservoir zones, multicomponent ocean‐bottom recording of converted‐waves (P to S or Ps) data provides reliable mapping of reservoir boundaries. Angle‐band stacks of multicomponent P‐wave (Pz) and Ps data can also be inverted jointly. In this paper Aki‐Richards equations are used without simplifications to invert angle‐band stacks to ‘reflectivities’. This enables the use of reflection seismic data beyond 30° of incident angles compared to the conventional amplitude variation with offset analysis. It, in turn, provides better shear impedance and density estimates. An important input to amplitude variation with offset analysis is the V_s/V_p ratio. Conventional methods use a constant or a time‐varying V_s/V_p model. Here, a time‐ and space‐varying model is used during the computation of the ‘reflectivities’. The V_s/V_p model is generated using well log data and picked horizons. For multicomponent data applications, the latter model can also be generated from processing V_s/V_p models and available well data. Reservoir rock properties such as λρ, μρ, Poisson's ratio and bulk modulus can be computed from acoustic impedance, shear impedance and density for pore fill and lithology identification. λ and μ are the Lamé constants and ρ is density. These estimations can also be used for a more efficient log property mapping. V_p/V_s ratio or Poisson's ratio, λρ and weighted stacks, such as the one computed from λρ and λ/μ, are good gas/oil and oil/water contact indicators, i.e., pore fill indicators, while μρ mainly indicates lithology. μρ is also affected by pressure changes. Results from a multicomponent data set are used to illustrate mapping of gas, oil and water saturation and lithology in a Tertiary sand/shale setting. Whilst initial log crossplot analysis suggested that pore fill discrimination may be possible, the inversion was not successful in revealing fluid effects. However, rock properties computed from acoustic impedance, shear impedance and density estimates provided good lithology indicators; pore fill identification was less successful. Neural network analysis using computed rock properties provided good indication of sand/shale distribution away from the existing wells and complemented the results depicted from individual rock property inversions.     

The geophysical significance of oceanic plagiogranite

Rocks ranging in composition from trondhjemite to diorite (plagiogranites) have been recovered from ocean ridges and are common constituents of ophiolites. Velocities and densities of diorite and trondhjemite from the Mid-Atlantic Ridge are shown to differ significantly from similar properties of metadolerite and gabbro. Compressional (V_p) and shear (V_s) velocities of plagiogranites are relatively low (V_p = 4.78–5.91km/s at1kbar,V_s = 2.81–3.37km/s at1kbar), as are densities (2.57–2.64 g/cm³) and Poisson's ratios (0.24–0.27). These data lend strong support to the probable existence of a low-velocity/density zone within layer 3 of the oceanic crust. Based on observations in ophiolites, it is postulated that this zone can be up to 1 km in thickness and is laterally discontinuous.     

IMPROVEMENTS IN THE OBSERVATION OF SHEAR WAVES*

The observation of shear waves is finding more interest in seismic exploration, especially for the determination of additional lithologic parameters, e.g. V^p/V^s. We demonstrate the observation of converted waves in routine seismic work by means of horizontal geophones. Field technique and data processing permit acquisition and interpretation of converted waves in connection with routine seismic measurements. A special interpretation of the recorded converted waves results in the shear wave velocity and the V^p/V^s relation.     

Stress wave patterns and the size of the non-elastic zones of an explosive source

Summary The travel-time curves of characteristic moments of stress wave patterns were investigated. The arrival time of the onset of stress wavet ₀, of the maximum amplitudet ₁ and of the moment terminating the pressure part (the first half-wave) of this wavet ₂ were taken to be representative. These travel-time curves were used to determine the velocities of propagation of these moments (V ₀,V ₁ andV ₂) as functions of the distance from the source. According to their variations it is possible to appoint the size of the cavity created by the explosion and to determine the distance at which the elastic source surface of stress waves is to be found. The radius of the cavity is given by 1.5 times the distance at whichV ₁ separates fromV ₀. The elastic source surface of stress waves is defined by the distance from whichV ₀ is constant and at whichV ₁ andV ₂ have minimum values. These two distances determine the points at which the character of the stress wave pattern changes: from shock wave in gases, through shock wave in a solid medium to a seismic wave.Paper presented at the General Assembly of the European Seismological Commission, Luxembourg, Sept. 21–29, 1970.     

Seismic Characterization of Neapolitan Soils

— Detailed shear-wave velocity profiles versus depth have been obtained in typical lithostratigraphies of Napoli. FTAN and hedgehog methods have been applied to Rayleigh surface waves recorded in refraction seismic surveys. The comparison with literature measurements shows good agreement with nearby down- and cross-hole tests. The pumiceous and lapilli content, and the different welding and alteration degree of the Neapolitan pyroclastic soils cause a strong scattering of the shear wave velocities (V_S) from bore-hole measurements, even for the same formation. Surface measurements, based on FTAN-hedgehog methods, determine average V_S along travel paths of about 100 m, give results that are comparable with down- and cross-hole velocity profiles, and have the additional advantage of being less scattered, and thus more representative of average properties than bore-hole measurements. The results of surface measurements should be preferred in the computation of realistic seismograms and are particularly suitable in urban areas, as they are not destructive and need just one receiver.     

P- and S-wave velocities of consolidated sediments from a seafloor seismic survey in the North Celtic Sea Basin, offshore Ireland

A geophysical survey was conducted over a hydrocarbon prospect in the North Celtic Sea Basin using a small array of ocean‐bottom seismographs (OBSs). The purpose of this study was to determine the ratio of compressional (P)‐ to shear (S)‐wave velocity of consolidated sedimentary rocks in order to constrain possible subsurface variations in pore‐fluid content. The ratio of V_P and V_S is known to be particularly sensitive to lithology, porosity and pore‐fluid content, making it a useful parameter for evaluating hydrocarbon prospects. OBSs offer a relatively cheap and time‐effective means of acquiring multi‐component data compared with ocean‐bottom cables. In this contribution, we demonstrate the ability of an OBS survey comprising three pairs of two OBSs spaced at 1.6 km to recover lateral variations in the V_P/V_S ratio. A key requirement of this type of study is that S waves will be generated by mode conversions in the subsurface, since they cannot be generated in nor travel through fluids. In this survey, the contrast in physical properties of the hard seabed of the North Celtic Sea Basin provided a means of generating converted S waves. Two‐dimensional ray‐tracing and forward modelling was used to create both V_P and V_S models along a profile crossing the Blackrock prospect in the North Celtic Sea Basin. These models comprise four layers and extend to a maximum depth of 1.1 km. The observed northward decrease in the V_P/V_S ratio at depths of 500–1000 m below the seafloor in the study area is interpreted to represent lateral variation in the amount of gas present in the pore space of Upper Cretaceous chalks and shales overlying the prospective reservoir.