用微动中的面波信息探测地壳浅部的速度结构

本文介绍了利用地脉动中的面波信息探测地壳浅部横波速度结构的观测方法、仪器设备、野外数据采集和用空间自相关方法测定相速度及相关的数据处理方法,并给出了在两个场地的探测实例. 通过与钻探结果对比表明,由该方法推测的地壳浅部1～2 km 深度范围的浅层横波速度结构,对于确定物性差异大的地层界面是有效和可靠的. 由于该方法使用的是地脉动噪声中的面波信息,不需要专门的震源, 因此具有经济快捷、应用范围广泛等优点,将可能成为一种新的具有很好应用前景的地球物理勘探方法.     

Correlation of cylindrical tank wall buckling with an earthquake motion recorded at a small distance from the tank

The bulging observed on the wall of a ground supported cylindrical storage tank, located at the Chiba Experimental Station of the University of Tokyo, is investigated. This damage is correlated with a ground motion recording obtained very near the tank base during a moderately strong earthquake that occurred on September 25, 1980 not far from the Chiba Experimental Station. From this correlation the importance of the horizontal excitation in amplitude and direction is clearly established. In addition, good agreement can be seen between the predicted tank wall earthquake performance, derived by an empirical approach proposed by the first author, and the shell deformation observed at the Chiba Experimental Station.     

地震波衰减规律及其恢复方法

针对地面地震分辨率低,不能有效识别薄层储层、薄层地层等问题,本文通过地震波衰减规律的分析,提出了一种恢复地震波高频衰减获得宽频带地震剖面的方法,从而大幅度地提高了地震分辨率.应用双井微地震测井资料对松辽盆地地震波动力学特征研究表明,地震波衰减规律是在近地表低速层和近地表低速层的近震源区地震波高频衰减巨大,而在高速层地震波高频衰减很小.药量大近震源区地震波高频衰减大,药量小近震源区地震波高频衰减小.那么,近震源区和近地表低速层对地震波的衰减是地面地震资料频带窄、分辨率低的主要原因.据此提出了如下确定性反褶积方法,用双井微地震测井资料求取近震源区、近地表低速层和虚反射等滤波因子,用其对地面地震资料作确定性反褶积处理,从而恢复近震源区、近地表低速层等几种因素的地震波衰减,将大药量激发地表接收的地面地震延拓成小药量激发高速层接收的宽频带地震.应用该方法对松辽盆地优势频带宽5~90 Hz,视主频50 Hz的地面地震资料处理后,地震剖面优势频带宽达5~360 Hz,视主频达180 Hz,使常规地震剖面分辨率提高2倍.具体的说松辽盆地中部含油组合的地面地震分辨能力由9~15 m提高到3~5 m.宽频带地震剖面与160 Hz的人工合成地震记录对比符合的很好,表明其处理结果是正确.该成果在油气勘探开发中的油气储层预测、构造学研究、沉积学研究等方面有重要的意义.     

Determination of the material damping ratio and the shear wave velocity with the Seismic Cone Penetration Test

The Seismic Cone Penetration Test (SCPT) is used to determine the variation with depth of the shear (and longitudinal) wave velocity at various sites in Belgium. The cross correlation technique together with a posteriori increase of the sampling frequency proves to be a reliable method to determine wave velocities by the SCPT. Additionally, the Spectral Ratio Slope (SRS) method is applied to Fourier transforms of the measured response to determine the variation of the material damping ratio with depth. This method is applied to synthetic signals computed for three different soil profiles with known characteristics and to the horizontal acceleration time histories that have been measured during SCPT on two sites. The influence of some signal processing techniques commonly applied in the SRS method is studied. A remarkable influence of the window length, applied on the predicted signals, on the computed material damping ratio is found, especially its scattering. It is therefore concluded that the use of a window function should be avoided wherever possible. On one of the two test sites, results of Resonant Column Tests and Free Torsion Pendulum Tests on undisturbed soil samples have been used to evaluate the values of the shear wave velocity and the material damping ratio as obtained in the in situ test.     

水库蓄放水对库底岩石介质弹性波速影响的数值模拟

为了解水库水体加、卸载及渗透过程中,库底岩石介质有效应力变化对弹性波速的影响,以及弹性波速变化与孔隙压力、有效应力变化之间的定量关系,以四川省紫坪铺水库为例,在前期已建立库区地质构造与水文地质结构模型,并推导出岩石介质变形与流体渗流耦合数学模型的基础上,进一步建立了有效应力对岩石弹性波速影响的数学模型,利用有限元方法定量计算了库底岩石介质的孔隙压力、附加有效应力和弹性波速的变化。结果表明:在水库水体加、卸载及渗透过程中,地下岩石介质弹性波速增量值的变化并不明显,最大变幅仅为±0.013km/s,变化范围主要集中在库底5km以上的局部区域;弹性波速增量值与附加有效应力变化同步,与水库水位变化形态相似;不同观测点弹性波速变化量的大小与埋藏位置、深度及附加有效应力的变化幅度有关,同一观测点,弹性纵波波速变化幅度大于弹性横波波速。文中数值模拟的结果与一些学者通过环境噪声方法获得的紫坪铺水库附近区域相对波速变化的反演结果在变化形态上大体一致,但变化幅度略有差异。     

Forward modeling of Rayleigh surface waves for analytical characterization of dominant dispersion trends

Forward modeling is of critical importance for inversion analysis of surface wave methods to obtain shear-wave velocity (V_S) profiles of soil sites. The dynamic stiffness matrix (DSM) method can provide forward modeling of Rayleigh surface waves to simulate complex wave propagation in layered soil sites. However, contamination from body waves and interference of multiple Rayleigh wave modes can reduce the accuracy of theoretical dispersion curves, especially at irregular soil sites with embedded low-velocity or high-velocity layers. An analytical method is developed herein to combine the techniques of the multichannel analysis of surface waves method with the DSM method to improve the accuracy of the theoretical dispersion analysis for soil sites. The proposed method implements multichannel analysis of the analytical displacement responses to capture dominant dispersion trends. Comparison of the results obtained with the new method against those from the transfer matrix method and the literature indicates that the new method can (1) effectively minimize the effects of contamination caused by body waves and interference from several Rayleigh wave modes, and (2) generate accurate dominant dispersion trends for soil sites with various stiffness profiles, especially for the high-frequency dispersion characteristics of the profiles with embedded low-velocity layers.     

S wave inelastic attenuation and site effect in the northern Tianshan area

In this paper,according to the relationship of observation spectrum with source spectrum,travel-path attenuationfactor and site effect,we use the digital seismic data of moderate-small earthquakes to invert the soft rock site ef-fect and S wave inelastic attenuation under the constraint of site effect,and determine the parameters of path andsite.The parameters obtained by this method exclude the topographic effect of seismic stations,and truly reflectthe seismic effect of free soft rock sites and the S wave inelastic attenuation in the studied region.The average softrock site effect is about 1.5 times in the frequency domain of 2～4 Hz and the S wave quality factor is obtained asQS=278 f0.346,which can be directly applied to the stochastic modeling of ground motions in the studied region.     

北天山地区S波非弹性衰减和场地效应研究

采用中小地震的数字观测资料,依据地震观测谱和震源谱、传播路径衰减因子及场地效应的关系式,在对场地效应约束的条件下,反演S波非弹性衰减和软基岩场地效应,确定研究区路径、 场地参数. 该方法确定的参数, 避免了台站地形效应的影响,真实地反映了研究区自由软基岩场地的地震效应和传播路径非弹性衰减的影响. 研究区的平均软基岩场地效应, 在2~4 Hz频率范围为1.5倍左右; S波品质因子为QS=278 f 0.346.这些参数可以直接用于研究区的地面运动随机模拟.      

用噪声成像方法研究云南地区强震前地震面波速度变化

本文利用2008年1月至2010年1月云南地震台网46个台站的连续波形资料, 采用背景噪声方法计算云南地区不同台站对之间的面波延迟时间, 通过层析成像方法反演得到不同时间段相同时间窗长的面波波速分布。 根据试验, 将三个月的时间窗长作为云南地区得到稳定信噪比的延迟时间的叠加天数。 结果表明, 不同时间段得到的波速结构具有明显的变化。 这些变化表明: ① 2008年汶川8.0级地震后, 程海断裂带波速增加, 能量积累加剧; ② 2008年攀枝花6.1级地震前攀枝花地区处于速度升高和速度降低的分界带上; ③ 2009年姚安6.0级地震前震中所在区域面波波速升高, 而震后出现降低。 与传统的动态成像相比, 背景噪声方法由于采用固定的台站对, 不同时间段的射线固定, 这样可以减小由于不同时段地震分布不均一和震源定位而引起的误差, 使动态反演结果更加可靠。     

Identification of soil degradation during earthquake excitations by Bayesian inference

A Bayesian inference approach is introduced to identify soil degradation behaviours at four downhole array sites. The approach of inference is based on a parametric time‐varying infinite impulse response filter model. The approach is shown to be adaptive to the changes of filter parameters and noise amplitudes. Four sites, including the Lotung (Taiwan), Chiba (Japan), Garner Valley (California), and Treasure Island (California) sites with downhole seismic arrays are analysed. Our results show two major types of soil degradation behaviour: the well‐known strain‐dependent softening, and reduction in stiffness that is not instantaneously recoverable. It is also found that both types of soil degradation are more pronounced in sandy soils than in clayey soils. The mechanism for the second type of soil degradation is not yet clear to the authors and suggested to be further studied. Copyright © 2003 John Wiley & Sons, Ltd.     

The use of MASW method in the assessment of soil liquefaction potential

The multi-channel analysis of surface wave (MASW) method is a non-invasive method recently developed to estimate shear wave velocity profile from surface wave energy. Unlike conventional SASW method, multi-station recording permits a single survey of a broad depth range and high levels of redundancy with a single field configuration. An efficient and unified wavefield transform technique is introduced for dispersion analysis and on site data quality control. The technique was demonstrated in the assessment of soil liquefaction potential at a site in Yuan Lin, Taiwan. The shear wave velocity and liquefaction potential assessments based on MASW method compares favorable to that based on SCPT shear wave measurements. Two-dimensional shear wave velocity profiles were estimated by occupying successive geophone spreads at several sites in central western Taiwan, at some of which sand boils or ground cracks occurred during 1999 Chi Chi earthquake. Liquefaction potential analysis based on MASW imaging was shown to be effective for estimating the extent of potential liquefaction hazard.     

基于正则化方法的高频地波雷达海浪方向谱反演

提出了一种从高频地波雷达海面回波谱中提取海浪方向谱的新方法:基于高频电波海洋探测基本原理将线性化后的非线性积分方程离散为矩阵方程组,通过引入正则化数学模型将不适定方程转化为正规方程,然后采用奇异值分解法求解.对于正则化方法中正则化参数,采用L曲线法来确定.为了验证本文算法的有效性,分别在单部雷达和双部雷达探测下进行了不同噪声水平下的数值模拟,结果表明了该正则化反演方法的有效性.文中还给出了实测数据分析初步结果.该算法的工程性应用还需进一步研究.对于实际的地波雷达海浪反演具有良好的应用前景.     

A simplified input output relation method using AR model for earthquake wave propagation analysis

This paper introduces a Simplified Input Output Relation Method (SIORM) using multiple-variable autoregression (AR) model which can be used to determine ground wave propagation properties. Using the AR model, a method is developed to establish the basis for the formulation of SIORM. The degree of accuracy of this method is evaluated. SIORM is applied to actual ground acceleration records taken at Minamisuna and Etchujima sites in Japan and its results are discussed in this paper. © 1997 John Wiley & Sons, Ltd.     

A strong motion database for the chiba seismometer array and its engineering analysis

A three-dimensional seismometer array was installed in the Chiba Experiment Station of the Institute of Industrial Science, University of Tokyo in 1982. The array system consists of 44 three-component accelerometers densely placed both on the ground surface and in boreholes. A complementary system for the measurement of ground and buried pipe strains was also installed at the same site. The array system has been successfully in operation, and more than 160 earthquakes have been recorded. Considering a wide use of these seismograms, the Chiba array database has recently been created comprising twenty-seven major events. This paper describes the Chiba array system and its strong motion database. Results of engineering analysis using the selected records are also presented.     

瑞利波法在工程抗震中的应用

工程抗震工作需要掌握各层土的剪切波速度资料,以便进行地震危险性分析和场地震害小区规划等项工作。采用瑞利波法进行岩土力学特性指标测定,通过实际性能应用表明该方法是简便易行的。     

Derivation of a curve number and kinematic-wave based flood frequency distribution

Abstract

The physically-based flood frequency models use readily available rainfall data and catchment characteristics to derive the flood frequency distribution. In the present study, a new physically-based flood frequency distribution has been developed. This model uses bivariate exponential distribution for rainfall intensity and duration, and the Soil Conservation Service-Curve Number (SCS-CN) method for deriving the probability density function (pdf) of effective rainfall. The effective rainfall-runoff model is based on kinematic-wave theory. The results of application of this derived model to three Indian basins indicate that the model is a useful alternative for estimating flood flow quantiles at ungauged sites.     

The effects of torsion and motion coupling in site response estimation

Soil amplification characteristics are investigated using data from the Chibaken‐Toho‐Oki earthquake and its aftershocks recorded at Chiba dense array in Japan. The frequency‐dependent amplification function of soil is calculated using uphole‐to‐downhole spectral ratio analysis, considering the horizontal components of shear wave. The identified spectral ratios consistently demonstrate the splitting of peaks in their resonance frequencies and low amplification values in comparison with a 1D model. The torsional behaviour and horizontal ground motion coupling are clarified as the reasons for these phenomena at the site. To prove the hypothesis, the torsional motion is directly evaluated using the data of the horizontal dense array in different depths at the site. The comparison between Fourier spectra of torsional motion and identified transfer functions reveals the peaks at the same frequencies. The wave equation including torsion and horizontal motion coupling is introduced and solved for the layered media by applying wave propagation theory. Using the developed model, the effects of torsional motion with horizontal motion coupling on soil transfer function are numerically examined. Splitting and low amplification at resonance frequencies are confirmed by the results of numerical analysis. Furthermore, the ground motion in two horizontal directions at the site is simulated using site geotechnical specification and optimizing the model parameters. The simulated and recorded motions demonstrate good agreement that is used to validate the hypothesis. In addition, the spectral density of torsional ground motions are compared with the calculated one and found to be well predicted by the model. Finally, the results are used to explain the overestimation of damping in back‐calculation of dynamic soil properties using vertical array data in small strain level. Copyright © 2003 John Wiley & Sons, Ltd.     

吴忠—灵武地区MS≥5.0地震前波速比变化特征分析

利用1980年1月-2008年12月宁夏地震台网记录的中小地震事件直达波到时数据,采用单台多震和达法,对吴忠-灵武地区该时段内M_S ≥ 5.0地震发生前波速比变化进行回溯性分析,并对该区未来短期内震情趋势进行研判,所得结果表明：①吴忠-灵武地区M_S ≥ 5.0地震前波速比存在正常→降低→上升(恢复)→发震的异常变化过程;②研究区台站出现波速比异常,其空间分布对发震位置具有一定指示意义,即发震位置可能位于波速比下降台站或存在异常台站所连线的多边形边缘或内部;③据近期扫描结果,海原地震台波速比异常较显著,需密切关注宁夏南部及邻区地球物理观测资料的变化。     

用浅层地震折射波法探测桥梁隧道所在地的基岩和构造

本文以广(州)—深(圳)高速公路(E)段的桥梁、隧道工程勘测为例,阐述了利用浅层地震折射波法在勘查断裂位置、基岩面产状,波速及岩性分布方面的应用。所得结果表明,该方法用于工程建设的选线、选址和概查是有效的。     

Rayleigh wave separation using high-resolution time-frequency polarization filter

Obtaining geological information from seismic data motivates researchers to innovate and improve seismic wave processing tools. Polarization-based methods have received much attention regarding their ability to discriminate between different phases of the seismic wave based on polarity. Combining the intuitive definition of polarity in the frequency domain (monochromatic waves) with the non-stationary properties provided by time-domain methods, time-frequency approaches are attracting widespread interest because they localize the information extracted from the seismic waves in the joint time and frequency domains. Due to the lack of high-resolution time-frequency maps, the time-frequency polarization approach was not able to resolve specific temporal polarity changes in the seismic signal. The main objective of this study was to devise a robust time-frequency-based polarization filtering method using high-resolution polarization attributes obtained directly from the sparse time-frequency map without using Eigen analysis or analytic signals. The method proposed here utilizes a computationally effective sparsity-based adaptive S-transform to obtain a high-resolution polarization map of an inherently non-stationary seismogram for the entire frequency content of the signal at different times. The superiority of the proposed method over the S-transform method was verified using synthetic and real data sets to calculate the polarization attributes in the time-frequency domain and separate the Rayleigh waves from the seismogram.