青藏高原地区瑞利波群速度和地壳构造

本文用单台多重滤波方法测定了经过青藏高原地区瑞利波群速度频散曲线。所得基阶瑞利波的观测周期为5.0-56.0秒,速度标准偏差为0.08-0.15公里/秒,一阶瑞利波的观测周期范围为10-16秒,速度标准偏差为0.05-0.13公里/秒。利用广义线性反演方法对频散曲线进行反演,可得出一个由五层构成的地震横波速度地壳模型。在27-40公里之间存在低速层,其横波速度为3.29公里/秒,比上一层低0.21公里/秒。     

Travel-time residuals and three-dimensional velocity structure of Italy

Jeffreys-Bullen P and PKP travel-time residuals observed at more than 50 seismic stations distributed along Italy and surrounding areas in the time interval 1962–1979, indicate the complex velocity pattern of this region. Strong lateral velocity inhomogeneities and low velocity zones are required to explain the observed pattern of residuals. In particular, late arrivals of about 1 sec are observed in the Apenninic mountain range, requiring both greater crustal thickness and low velocity layers, coherent with seismic refraction data and surface wave dispersion measurements. The seismic stations located in the Western and Eastern Alps indicate the presence of high velocities. In the Western Alps the strong azimuthal variation of residuals and the high values of early arrivals have a close relationship to the Ivrea body, an intrusive crustal complex characterized by a velocity as high as 7–7.2 km/sec.A travel-time inversion performed with theAki et al. (1977) block model, confirms the peculiar characteristics and the sharp variations in the lithosphere of the whole Italian region, with values of velocity perturbations between many adjacent blocks, ranging in size from 50 to 100 km, and independent from the earth parametrization chosen, reaching values up to 10% in the lithospheric part and 5% in the asthenosphere. 3-D inversion requires also high velocity along the Tyrrhenian coastal margin, equivalent to an uprise of major crustal and lithospheric discontinuities along this part of the Italian peninsula. Moreover low velocity material must be present in the northern part of the Adriatic foreland, in the lithosphere-asthenosphere system, closely related to the stress and seismicity pattern, and the lateral bending of the lithosphere in the same region.     

DECOMPOSITION AND INVERSION OF SEISMIC DATA–AN INSTANTANEOUS SLOWNESS APPROACH*

Interpretation techniques are presented that aim at the estimation of seismic velocities. The application of localized slant stacks, weighted by coherency, produces a decomposition of multichannel seismic data into single trace instantaneous slowness p(x, t) components. Colour displays support the interpretation of seismic data relevant to the near surface velocity structure. Since p(x, t) is directly related to stacking velocities and the depth of reflection, or bottoming points, in the subsurface, this data transformation provides a powerful tool for the inversion of reflection and refraction data.     

THE TERRAPAK *

The Terrapak is a new seismic energy source of the surface energy source type. Essentially, the Terrapak consists of a compressed air piston which transmits a shock to the ground by means of a striker plate attached to the piston rod. The Terrapak is truck-mounted together with its ancillary equipment. Shocks can be generated at 8 second intervals, and although the weight of the striker is relatively low (approximately 1400 Ibs) the impact velocity exceeds 66 ft/sec. The energy available is about 85,000 ft. lb., which is equivalent to that produced by conventional weight-dropping methods. The power spectrum of the seismic wave produced by the Terrapak is spread out towards the higher frequencies because its impact velocity is greater than that of a dropping weight. There is a noticeable difference between noise records made with the Terrapak and noise records made with a conventional weight-dropping technique. The same noise frequencies are found in both cases, but low frequency noise is stronger with weight dropping than with the Terrapak, whereas the intermediate frequency reflections obtained using the Terrapak are more suitable for correlation. The expansion of the compressed air generates a sound wave, the effects of which must be eliminated by the use of multiple geophone arrays. A spectrum analysis of the records has confirmed that the waves produced by the Terrapak are rich in intermediate frequency components. This facilitates the picking of reflections in normal areas where the apparent frequency of the reflections is about 30 to 35 c/s.     

考虑震源破裂方向的地震动衰减模型

本文基于Haskell的3源辐射谱模型,加上介质非弹性吸收项,几何扩散项和自由面放大效应,推演出峰值加速度和均方根加速度表达式。据此,建议一种考虑震源破裂方向的地震动衰减模型。用Morgan Hill和Imperial Valley两地震的峰值加速度资料,单震级统计分析结果表明,文中建议的地震动衰减模型是合理的。与Joyner和Boore所用的衰减模型相比较,剩余标准差可减少0.1——0.2。本文所考虑的震源破裂的地震动衰减模型,可分别用于不同断层类型的地震,如单侧和双侧走滑型地震及倾滑型地震,分别建立地震动衰减关系,进行地震危险性分析。文章最后给出了美国西部考虑震源破裂方向的地震动峰值加速度和均方根加速度的衰减关系。      

SEISMIC REFRACTION MEASUREMENTS IN THE BALTIC SEA*

Seismic refraction measurements were made in two areas of the Baltic Sea in June 1967. The refraction data were obtained in the course of the transmission measurements program of Operation MILOC BALTIC 67. Three profile pairs were obtained, two in the area south of Öland Island, Sweden, the third to the east, north of the peninsula of Hel, Poland. The water depths vary from 60 to 90 meters between the areas. The receiving positions for the two profiles south of Oland Island are only 46 km apart but the structures differ markedly. The northern section, only 500 m thick, shows a rather thin sedimentary cover above 350 m of 3.7 km/sec material that in turn overlies 5.6 - 5.9 km/sec velocity material. The southern section, almost 2 km thick, has an equivalent amount of low velocity material, and a layer about 1 km thick having velocity 4.8 km/sec that overlies 6.0 km/sec velocity material. The eastern profile shows high velocity material, 5.6 km/sec, at 2.5 km depth. Correlation of the layers determined by seismic refraction with nearby geology suggests that the structural change south of Öland Island may represent the boundary of the Sarmatian Shield in this region.     

VERTICAL HETEROGENEITY AND MOVEOUT IN THE ONE-DIMENSIONAL MEDIUM*

Since the important contributions of Dürbaum and Dix, 30 years ago, velocity profile estimation procedures on horizontally layered and vertically heterogeneous media from seismic probing data have been based largely on hyperbolic moveout models and RMS and stacking velocity concepts. Re-examination of the fundamentals reveals that quantitative velocity heterogeneity and canonical valocity profiles have been implicit factors for moveout modelling and for profile inversion in the use of the Dix procedure. Heterogeneity h is the ratio (and v_RMS the geometric or harmonic mean) of the path-average and time-average velocities for a raypath or, in a more restricted sense, for the normal ray belonging to a velocity profile. The canonical profile for a given velocity profile or profile segment is a moveout-equivalent monotonically increasing ramp-like profile. The ramp or constant gradient in depth is the simplest velocity profile approximator which can explicitly accommodate velocity heterogeneity. A ramp model structure is detailed which facilitates moveout simulation and model parameter estimation, and the parametric effects are explored. The horizontal offset range is quantified for which this model can give good moveout approximations.     

地震波速与地壳岩石生热率关系的研究评述

本文概述了由地震波速（Ｖｐ）计算地壳深部岩石生热率（Ａ）这一方法产生的历史背景，介绍了提出Ｖｐ－Ａ关系的初始理论，Ｖｐ－Ａ关系的各种计算式、适用区域及温压校正方法，综述了近十年来围绕这一关系而展开的争论，报道了研究新进展，提出了地壳深部岩石生热率的研究方法。     

祁连山—河西走廊地壳速度结构及速度与电性的联合解释

联合利用天然地震和人工地震测深资料对祁连山—河西走廊地区进行地壳三维速度结构反演,获得了该区地壳三维速度结构图像．给出了按5km间隔的中、上地壳速度切面和地壳速度结构与莫霍界面等深线图．综合大地热流实测结果和深部岩石物性等,建立了下地壳部分熔融的韧性模型．给出了计算地震波速度和电阻率的方法与公式,进行了部分熔融状态下深部电性与地震波速度的联合数值模拟和联合解释,其结果与实测的基本一致．分析了各构造单元之间下地壳电阻率、速度与热流的关系,对下地壳部分熔融程度进行了估计．     

Handlotungen und Echolotungen in der Seevermessung geringer Wassertiefen in strömenden Gewässern — eine kritische Betrachtung

Zusammenfassung Vergleichslotungen zwischen Echolot und Handlot führen zu einer mathematischen Ableitung der Handlotleinen-Kurve in strömenden Gewässern. Durch Einführung von zwei Vereinfachungen — Vernachlässigung des Lotleinen-Gewichts und Annahme der Konstanz der Strömungsgeschwindigkeit nach der Tiefe — wird theoretisch die Kettenlinie, praktisch eine Parabel ermittelt. Als bestimmend für den Parameter der Kettenlinie ergeben sich das Lotgewicht und der Staudruck auf die Lotleine. Aus einem praktischen Beispiel wird der Parameter entnommen, und es werden für diese Verhältnisse die für die Tiefenwerte bis 26 m zu erwartenden Fehlmessungen des Handlotes sowie die Größe der Ablagen der Lotbirne vom Lotpunkt berechnet. Für andere Werte der Strömungsgeschwindigkeit werden die Veränderungen auf die Ausbuchtung der Lotleine ebenfalls rechnerisch bestimmt und in Kurven dargestellt. Aus diesen Kurven ergeben sich zwangsläufig die größten Tiefen, die man bei verschiedenen Strömungsgeschwindigkeiten mit dem schweren Handlot vermessen darf; sie liegen bei v=1 m/sec bei 8 m, bei v=0,7 m/sec bei 12 m und bei v=0,5 m/sec bei 20 m Wassertiefe.Praktische Beispiele der Seevermessung aus der jüngsten Zeit zeigen deutlich die Notwendigkeit, zukünftig nur noch mit dem Echolot zu vermessen.

---

A comment on the applicability of hand leads and echo sounders to surveying operations in flowing shallow waters

Summary From comparisons between hand-lead soundings and echo soundings in streaming waters the divergency of the lead line from the vertical is mathematically derived. In introducing two simplifications, i. e. in neglecting the weight of the lead line and in assuming the velocity of current to be constant at all depths, the catenary curve which practically speaking is a parabola, is theoretically ascertained. The weight of the lead and the dynamic pressuro of the water (the pressure as exerted on a jet and measured by it) are found to be decisive for the parameter of the catenary-curve. The parameter is derived from a practical example and, according to the conditions of same, the errors to be expected in sounding with the hand lead at depths up to 26 m (= 79 ft) as well as the divergency of the sounding lead from the vertical are calculated. The various degrees of the lead line's divergency caused by a greater or lesser intensity of the current are likewise calculated and represented by curves. These curves necessarily indicate the maximum depths that can be sounded with sufficient accuracy with the heavy hand lead (5 kilos) at different velocities of current; the approximate values are as follows: 8 m for 1,0 m/sec, 12 m for 0,7 m/sec, 20 m for 0,5 m/sec (or: 25 ft for 2,0 knots, 37 ft for 1,5 knots, 61 ft for 1,0 knot). As a result of practical experience in recent survey operations the echo sounder should be regarded as the only reliable instrument to be used in future surveying.

---

Note critique sur l'applicabilité de la sonde au plomb-poisson et celle du sondeur acoustique aux levés hydrographiques effectués en faibles profondeurs dans les eaux courantes

Résumé Les comparaisons des sondages au plomb-poisson avec ceux effectués par écho permettent d'en dériver mathématiquement la courbe formée par la ligne de sonde dans des eaux courantes. En introduisant deux simplifications, c'est à dire en négligeant le poids de la ligne et en supposant que la vitesse de courant soit constante dans toutes les profondeurs on arrive à déterminer théoriquement la chaînette qui est en effet une parabole. Le paramètre est fonction du poids du plomb-poisson et celle de la poussée sensiblement horizontale sur le fil due au frottement du liquide. On se sert d'un exemple pratique pour en dériver le paramètre, pour calculer, dans les conditions de cet exemple, les sondages faux de la sonde au plomb-poisson dont il faut tenir compte dans les profondeurs jusqu'à 26 m, et pour déterminer la déviation entre le plomb au bout de la ligne et la verticale. Les diverses variations de la courbe, que la ligne subit et qui sont causées par d'autres valeurs de la vitesse de courant sont aussi calculées et représentées à l'aide des courbes. Ces courbes indiquent forcément les profondeurs permettant encore une précision suffisante du sondage au plombpoisson lourd (5 kg) aux différentes variations de la vitesse de courant. Ces maxima sont à peu près comme suit:v=1,0 m/sec à 8 m de profondeur,v=0,7 m/sec à 12 m de profondeur,v=0,5 m/sec à 20 m de profondeur. Des exemples pratiques offerts par des récentes levés hydrographiques prouvent que le sondeur acoustique doit être considéré comme le seul instrument assurant le précis levé hydrographique.

     

Near surface velocity and QS structure of the Quaternary sediment in Bohai basin,China

Heavily populated by Beijing and Tianjin cities, Bohai basin is a seismically active Cenozoic basin suffering from huge lost by devastating earthquakes, such as Tangshan earthquake. The attenuation (Q_P and Q_S) of the surficial Quaternary sediment has not been studied at natural seismic frequency (1?10 Hz), which is crucial to earthquake hazards study. Borehole seismic records of micro earthquake provide us a good way to study the velocity and attenuation of the surficial structure (0?500 m). We found that there are two pulses well separated with simple waveforms on borehole seismic records from the 2006 M_W4.9 Wen'an earthquake sequence. Then we performed waveform modeling with generalized ray theory (GRT) to confirm that the two pulses are direct wave and surface reflected wave, and found that the average v_P and v_S of the top 300 m in this region are about 1.8 km/s and 0.42 km/s, leading to high v_P/v_S ratio of 4.3. We also modeled surface reflected wave with propagating matrix method to constrain QS and the near surface velocity structure. Our modeling indicates that Q_S is at least 30, or probably up to 100, much larger than the typically assumed extremely low Q (~10), but consistent with Q_S modeling in Mississippi embayment. Also, the velocity gradient just beneath the free surface (0?50 m) is very large and velocity increases gradually at larger depth. Our modeling demonstrates the importance of borehole seismic records in resolving shallow velocity and attenuation structure, and hence may help in earthquake hazard simulation.     

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.     

Soil-structure interaction effects on the seismic response of tall chimneys

The effects of soil-structure interaction on the seismic response of tall (>100m) steel and reinforced concrete chimneys are described. Detailed models of a 130m high steel chimney and a 150m high reinforced concrete chimney are used as structural models. The foundations are represented as rigid blocks resting on a uniform viscoelastic soil model. Perfect bonding between the foundation and the soil is assumed. Parametric studies of the interaction effects on the magnitude and distribution of bending moments and shear forces include four soil rigidities and two seismic excitations characterized by very different frequency contents. The results obtained indicate strong interaction effects for intermediate and soft soils (V_s500 m/sec). The extent of the interaction effects are highly dependent on the characteristics of the seismic excitation.     

Site Response from Ambient Vibrations in the Towns of Lod and Ramle (Israel) and Earthquake Hazard Assessment

Throughout their long history, the towns of Lod and Ramle have been severely affected by strong earthquakes. The last destructive earthquake occurred on July 11, 1927 and caused the destruction of large parts of these cities, reaching a seismic intensity of VIII–IX on the MSK scale. Such a high intensity from a relatively distant earthquake (about 70 km) of magnitude 6.2 is likely to be the result of local site effects of the sedimentary layers that may have significantly enhanced earthquake ground motions. This study is focused on estimating the seismic hazard to Lod and Ramle by implementing a three-step process (1) detailed mapping of the characteristics of the H/V spectral ratios from ambient noise, (2) incorporating geological information and well data to construct subsurface models for different sites within the investigated area and (3) estimating the seismic hazard in terms of uniform hazard site-specific accelerations. The horizontal-to-vertical spectral ratios of ambient noise were used to approximate the fundamental resonance frequencies of the subsurface and their associated amplitudes. About 360 sites in Lod and Ramle were instrumented for varying periods. The soil sites exhibits H/V peak amplitudes ranging from 4 to 6 in the frequency range 0.5–2.5 Hz. These data were used to constrain 1-D subsurface models that were developed using geological data and borehole information. H/V spectral ratio observations were checked against theoretical subsurface transfer functions at locations where borehole information is available farther constraint the range of possible V_s velocities of the different layers and thus, by means of trial an error it was possible to conclude a systematic spatial distribution of the V_s velocity and thickness in the substrata that are also consistent with the spatial distribution of the fundamental resonance frequencies of the soft sediments obtained by means of the H/V spectral ratios, and other geological and geophysical information available at different locations in the study area. The evaluated subsurface models are introduced using the SEEH procedure of Shapira and van~Eck [(1993) Natural Hazards 8, 201–205] to assess Uniform Hazard Site-Specific Acceleration Spectra for different zones within the towns of Lod and Ramle. These evaluations are very important for realistic assessment of the vulnerabilities of all types of existing and newly designed structures and for urban and land use planning.     

我国西北地区地壳中的高速夹层

在我国西北地区的柴达木盆地东部和甘肃地区,在距离炮点40互100公里处,能够接收到不少能量较强的地壳深界面反射波。另外还发现一种与一般反射波性质不同的波,其视速度特大,视速度随距离的变化不大,而且有较明显的终点;其吋距曲线与一般深界面反射波的时距曲线相交。根据它的特征可以判断地壳中存在具有速度梯度的高速夹层.求得的夹层参数为: 甘肃地区柴达木盆地东部覆盖层厚度 18.8公里 30.5公里覆盖层平均速度 5.5公里/秒 5.3公里/秒夹层厚度 6.0公里 3.2公里夹层速度 7.5-8.5公里/秒 7.5-8.0公里/秒夹层的上下界面均为强反射面,可以产生多次反射波。分別利用相邻两个反射波可以求得各层参数,并能避免射线折射的影响。甘肃地区和柴达木盆地东部的地壳厚度分別为51和52公里。地壳中有高速夹层的存在,可以更好地说明P~*速度分散的原因,而且也能够解释Lg波的传播机制。     

Relative source time functions of seismic events at the Rudna copper mine, Poland: estimation of inversion uncertainties

The relative source time function (RSTF) inversion uncertainty assessment was performed for two small, mining-induced seismic events (M _W =2.9 and 3.0) that occurred at Rudna copper mine in Poland. The seismograms of selected events were recorded by the seismic net work composed of over 60, short-period, vertical seismometers, recording ground velocity, located in the distance ranging from 400 m up to 8 km from their hypocenters. The RSTFs were calculated for each seismic station independently, using the empirical Green’s function technique. The pseudospectral approximation of the sought RSTF by a finite sum of Gaussian kernel functions was used and the inverse problem was solved with the adaptive simulated annealing algorithm. Both methods improved the stability of the deconvolution procedure and physical correctness of the final solution in comparison to the classical deconvolution methods. To estimate the inversion uncertainties, classical Markov-chain Monte-Carlo techniques were used. The uncertainty analysis allows for improved selection of a priori data to the following inversion for kinematic rupture process.     

Crustal structure and temporal velocity change in southern California

Summary About twenty blasts are used to determine crustal structure and to monitor temporal seismic velocity changes in southern California. The shot time is determined up to 10 msec by using a disposable pick-up placed directly on the explosive. About 17 permanent stations and 20 temporary stations are used for the recordings. With a fast paper speed (typically 1 cm/sec) and the WWVB radio signals superposed on the seismic trace, absolute timing accuracy of up to 10 msec is achieved. A representative structure thus determined consists of a 4 km thick 5.5 km/sec layer underlain successively by 23.4 km thick 6.3 km/sec layer, 5.0 km thick 6.8 km/sec layer and 7.8 km/sec half space. The details of the lower crust are somewhat uncertain. This structure can explain the travel time data, corrected for the station and source elevations and for the station delays, to ±0.15 sec. Small but systematic temporal velocity changes up to 3% have been found for some of the profiles. If the effect of the migration of the shot point is small enough, these changes are larger than experimental errors and represent real temporal change in the material property between the shot point and the stations.Contribution No. 2530, Division of Geological Sciences, California Institute of Technology, Pasadena, California.     

DERIVATION OF LAYER PARAMETERS OF AN ELASTIC MEDIUM FROM REFLECTION COEFFICIENT MATRICES*

In this paper we develop a recursive algorithm to obtain the layer parameters of an elastic medium (density, P-wave velocity, S-wave velocity) from reflection coefficient matrices in terms of energy flux ratios for a non-normal incidence case. We define a layer impedance matrix, analogous to the impedance of an acoustic medium. Next we derive a matrix relationship between the layer impedance matrix of the n+ 1st layer and the reflection coefficient and parameter matrices of the nth layer. This relationship leads to recursively computing the parameters of the subsurface. We show that the elastic case—unlike the acoustic case—allows one to recover the layer parameters from the impedance matrix for non-normal incidence. The results of this work play a key role in the solution of the inverse problem with non-normal-incidence plane-wave seismic data when using a downward continuation technique.     

REFLECTION-WINDOW MAPPING OF SHALLOW BEDROCK*

A field study was undertaken to evaluate the effectiveness of the high frequency seismic reflection technique for mapping of shallow and irregular bedrock. Bedrock reflections were obtained using a hammer source with both in-line and common offset field layouts. The recording equipment included 12-channel enhancement seismographs, 28 Hz vertical geo-phones and a microcomputer. The latter increased the overall versatility of the seismic system. Field sites for this study are typical of the geological settings of the tin mining areas of Malaysia. The topographical ‘lows’ of the irregular bedrock control the localization of tin ore. The subsurface geology consists of a thin low velocity layer (± 300 m/s) overlying the compact overburden (± 1700 m/s) which in turn lies on bedrock. This paper discusses various criteria for designing an optimum window for obtaining usable reflections between the first arrival and the leading edge of the ground roll cone. Detailed mapping of the overburden and the bedrock interface by the reflection method can be useful in delineating areas for exploratory drilling and for optimum planning of mining operations.     

P‐wave velocity distribution in basalt flows of the Enni Formation in the Faroe Islands from refraction seismic analysis

The main objective of this work is to establish the applicability of shallow surface‐seismic traveltime tomography in basalt‐covered areas. A densely sampled ～1300‐m long surface seismic profile, acquired as part of the SeiFaBa project in 2003 ( Japsen et al. 2006 ) at Glyvursnes in the Faroe Islands, served as the basis to evaluate the performance of the tomographic method in basalt‐covered areas. The profile is centred at a ～700‐m deep well. V_P, V_S and density logs, a zero‐offset VSP, downhole‐geophone recordings and geological mapping in the area provided good means of control. The inversion was performed with facilities of the Wide Angle Reflection/Refraction Profiling program package ( Ditmar et al. 1999 ). We tested many inversion sequences while varying the inversion parameters. Modelled traveltimes were verified by full‐waveform modelling. Typically an inversion sequence consists in several iterations that proceed until a satisfactory solution is reached. However, in the present case with high velocity contrasts in the subsurface we obtained the best result with two iterations: first obtaining a smooth starting model with small traveltime residuals by inverting with a high smoothing constraint and then inverting with the lowest possible smoothing constraint to allow the inversion to have the full benefit of the traveltime residuals. The tomogram gives usable velocity information for the near‐surface geology in the area but fails to reproduce the expected velocity distribution of the layered basalt flows. Based on the analysis of the tomogram and geological mapping in the area, a model was defined that correctly models first arrivals from both surface seismic data and downhole‐geophone data.