研究局部场地条件对地震动影响的响嘡遥测台阵

对原有的响嘡三维场地条件影响台阵进行了扩建改造。增设了一个深度为47m的测井,对钻孔土样进行了土动力学测试,弥补了缺乏详细土层特性资料的缺陷。新老测井观测点均架设了电话线路,实现了观测中心与台阵之间的远程监控和数据传输。已经获取的一批中小地震近场记录可以用来进行盲测研究。     

Shear-wave velocity estimation using a combination of ambient noise from small aperture array and small-scale active seismic measurements: a case study in the area of the natural gas fields of Northern Germany

In recent years, numerous induced seismic events have occurred in the proximity of the natural gas field in Northern Germany. To monitor the seismicity and to asses the seismic hazard potential, a local monitoring network was installed in the area. Focusing on the seismicity hazard assessment, a major challenge is the characterisation of potential site effects due to local soil characteristics. This is quantitatively performed by estimating the shear-wave velocity (V _s) variation with subsurface layer thickness. Such local effects can only be covered with a coarse spatial resolution due to the limited number of monitoring stations. Profiles were determined at three test sites (Langwedel, Walle and Bomlitz) by using a combined approach of small aperture 2D array ambient noise and small-scale active 1D measurements. The high-resolution frequency-wavenumber (HRFK), spatial autocorrelation (SPAC) and multichannel analysis of surface waves (MASW) methods were applied to the recorded ambient noise and active seismic data using various array sizes supplemented by the active measurements. This jointly allowed obtaining phase velocity dispersion curves covering a wide frequency range from 2 up to 32 Hz. The inversion of the obtained dispersion curves results in average S-wave velocity profiles down to depths of 70 m, identifying thin near-surface layers of a few meters as well as thicker layers of tens of meters in greater depth. A comparison with available borehole data shows a good correlation with the layering. Additionally, to asses the impact of a seismic event at the test sites, PGV estimations for various seismic events were performed. The final results of the test surveys demonstrate that the combined approach represents a suitable tool for near-surface characterisation, which can be used to improve the seismic hazard assessment in the area of the natural gas fields in Northern Germany.     

场地土层模型参数的地震动记录反演方法

考虑场地土层不同物理参数的综合反演,改进了水平与竖向谱比(Horizontal-to-Vertical Spectral Ratio,简称HVSR)混合全局优化反演方法,进一步通过土层的S、P波波速、厚度、泊松比、密度和S、P波品质因子等土层参数反演的敏感性分析,形成了可同时反演场地土层厚度及剪切波速的混合全局优化反演方法.以美国GVDA和日本IWTH27竖向强震动观测台阵为例,分别以理论HVSR及加速度观测记录获得的HVSR曲线为目标,反演获得了场地浅层速度结构,并与观测台阵场地钻孔揭示的土层模型进行比较,验证了发展的反演方法的合理性和适用性.本文研究表明,基于地震加速度记录的HVSR全局优化反演方法是获取场地土层速度结构的一种有效的途径.     

Designing and testing a network of passive seismic surveying and monitoring in Dehdasht (South Western Iran)

From August 2016 to July 2017, a passive seismic survey was conducted in South Western Iran as a part of a pilot project aimed to improve the imaging in geologically complex areas. Passive seismic methods have shown to be a useful tool to infer the physical properties of the underground geological structures where traditional hydrocarbon exploration methods are challenging. For this purpose, a dense passive seismic network consisting of 119 three-component borehole seismic stations was deployed over an area of 400 km² around the city of Dehdasht. This paper focuses on the details of the network design, which was devoted to high-resolution seismological applications, including local earthquake tomography and seismic attenuation imaging. In this regard, we describe the instrument types and the station installation procedures used to obtain high-quality data that were used to retrieve three-dimensional models of P- and S-wave velocity and P-wave attenuation in the area using tomographic inversion techniques. We also assess the network performance in terms of the seismic ambient noise levels recorded at each station site, and we revise the horizontal orientation of the sensors using surface waves from teleseismic earthquakes.     

Inverse analysis of weak and strong motion downhole array data from the Mw7.0 Sanriku-Minami earthquake

A hybrid optimization scheme, comprising a genetic algorithm in series with a local least-squares fit operator, is used for the inversion of weak and strong motion downhole array data obtained by the Kik-Net Strong Motion Network during the M_w7.0 Sanriku-Minami Earthquake. Inversion of low-amplitude waveforms is first employed for the estimation of low-strain dynamic soil properties at five stations. Successively, the frequency-dependent equivalent linear algorithm is used to predict the mainshock site response at these stations, by subjecting the best-fit elastic profiles to the downhole-recorded strong motion. Finally, inversion of the mainshock empirical site response is employed to extract the equivalent linear dynamic soil properties at the same locations. The inversion algorithm is shown to provide robust estimates of the linear and equivalent linear impedance profiles, while the attenuation structures are strongly affected by scattering effects in the near-surficial heterogeneous layers. The forward and inversely estimated equivalent linear shear wave velocity structures are found to be in very good agreement, illustrating that inversion of strong motion site response data may be used for the approximate assessment of nonlinear effects experienced by soil formations during strong motion events.     

Application of passive source surface-wave method in site engineering seismic survey

Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m underground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes, shorten survey period, and reduce engineering cost to some extent.     

An inversion of site response and Lg attenuation using Lg waveform

Based on spectral ratio method, a joint inversion method was used to obtain parameters of Lg wave attenuation and site response. The inversion method allows simple and direct (two-parameter) determination of Lg wave attenua- tion and site response from sparse spectral data, which are not affected by radiation pattern factor and different response of same instrument after geometrical spreading. The method was used successfully for estimating site re- sponse of stations of Zhejiang Seismic Network and measuring Lg wave attenuation. The study is based on 20 earth- quakes occurred in northeast of Taiwan with magnitude MS5.0~6.7 and 960 seismic wave records from 16 stations in Zhejiang area from 2002 to 2005. The parameters of site response and Lg attenuation were calculated with a fre- quency interval of 0.2 Hz in the range of 0.5 Hz to 10 Hz. Lg wave attenuation coefficient corresponding to U-D, E-W and N-S components are γ ( f )=0.00175 f 0.43485, γ ( f )=0.00145f 0.48467 and γ ( f )=0.0021f 0.41241, respectively. It is found that the site response is component-independent. It is also found that the site response of QIY station is significant above the frequency of 1.5 Hz, and that the site response of NIB station is low for most frequency     

Assessing uncertainty in refraction seismic traveltime inversion using a global inversion strategy

To analyse and invert refraction seismic travel time data, different approaches and techniques have been proposed. One common approach is to invert first‐break travel times employing local optimization approaches. However, these approaches result in a single velocity model, and it is difficult to assess the quality and to quantify uncertainties and non‐uniqueness of the found solution. To address these problems, we propose an inversion strategy relying on a global optimization approach known as particle swarm optimization. With this approach we generate an ensemble of acceptable velocity models, i.e., models explaining our data equally well. We test and evaluate our approach using synthetic seismic travel times and field data collected across a creeping hillslope in the Austrian Alps. Our synthetic study mimics a layered near‐surface environment, including a sharp velocity increase with depth and complex refractor topography. Analysing the generated ensemble of acceptable solutions using different statistical measures demonstrates that our inversion strategy is able to reconstruct the input velocity model, including reasonable, quantitative estimates of uncertainty. Our field data set is inverted, employing the same strategy, and we further compare our results with the velocity model obtained by a standard local optimization approach and the information from a nearby borehole. This comparison shows that both inversion strategies result in geologically reasonable models (in agreement with the borehole information). However, analysing the model variability of the ensemble generated using our global approach indicates that the result of the local optimization approach is part of this model ensemble. Our results show the benefit of employing a global inversion strategy to generate near‐surface velocity models from refraction seismic data sets, especially in cases where no detailed a priori information regarding subsurface structures and velocity variations is available.     

Comparison of spectral ratio techniques for estimation of site effects using microtremor data and earthquake motions recorded at the surface and in boreholes

The authors apply soil response estimation techniques employing accelerograms for fifteen earthquakes recorded at the Yokohama Strong Motion Network and its vertical array of nine sites, plus microtremor data recorded at all 150 sites. Assessment of the reliability of surface to reference site spectral ratios, and horizontal to vertical spectral ratios of S‐waves, coda, and microtremors, relied on cross‐validation with both surface to borehole spectral ratios, and theoretical soil response functions for vertically incident SH plane waves calculated from the vertical array recordings and logging data, respectively. The results indicate the following: (1) the reference site method provides more accurate soil response estimations; (2) the H/V of microtremors provides relatively better soil response estimations than its S‐wave and coda counterparts; and (3) the H/V of microtremors provides coincident soil response estimations at two‐thirds of the sites when considering no more than ±30% mismatch in the determination of the predominant periods. The matching of parameters obtained from the analysis of microtremors, combined with their relatively low economic cost and the possibility of recording without strict spatial or time restrictions turns microtremors into a particularly appealing approach. Copyright © 2003 John Wiley & Sons, Ltd.     

浅硬场地剪切波速测试标准差对地震动影响初探

采用Shake 2000程序,以Turkey Flat试验场地为模型,通过输入不同类型下多种强度的地震波,计算研究多工况下剪切波速测试标准差对地表加速度反应谱和峰值加速度的影响。结论为:(1)浅硬场地上剪切波速测试标准差对地震动的影响很大,影响程度与输入地震波的强度和频率以及场地剪切波速计算值有关;(2)如果将反应谱残差大于20%或加速度峰值差别大于20%定义为统计意义上的不可忽略,那么剪切波速测试标准差对计算结果的影响在大多数情况下均不可忽略;(3)当输入波的卓越周期与场地特征周期接近时,浅硬场地上剪切波速测试标准差引起的反应谱变化非常显著;(4)只有当输入波的卓越周期与场地特征周期相差较大且输入波强度偏小时,剪切波速测试标准差引起的反应谱变化才可略去;(5)当浅硬场地上剪切波速实测结果低于统计均值时,地震动计算结果的偏差一般明显大于剪切波速实测结果,高于统计均值时引起的偏差,且地震输入越强表现越明显。     

渤海海域黏性土剪切波速与抗剪强度统计关系的初步研究

剪切波速与地基土的抗剪强度、剪切模量和卓越周期等参数密切相关,是地震安全性评价中判定场地类别的一个主要指标和参数。鉴于海域工程中剪切波速往往难以直接由原位测得,而室内实验结果又常常与野外现场物探测试值存在较大差异,因此,如何通过其他途径有效获取满足工程需要的剪切波速参数在海域工程的地震安全性评价等方面具有迫切的实用需求。为此本文通过对渤海海域数十个石油平台项目中一系列饱和黏性土样品的剪切波速与抗剪强度实验数据的统计分析,尝试采用多种可能的函数来拟合确定二者之间的经验关系。结果表明:对于渤海海域黏性土剪切波速V_s与抗剪强度S_u之间的最佳统计经验关系为幂函数V_s=53.751S_u~(0.376)。此关系可为渤海海域工程中通过不排水抗剪强度估算剪切波速提供一种简便可行的实用性方法。     

Multiwell study of seismic attenuation at the CO2CRC Otway project geosequestration site: Comparison of amplitude decay,centroid frequency shift and 1D waveform inversion methods

At the CO2CRC Otway geosequestration site, the abundance of borehole seismic and logging data provides a unique opportunity to compare techniques of Q (measure of attenuation) estimation and validate their reliability. Specifically, we test conventional time-domain amplitude decay and spectral-domain centroid frequency shift methods versus the 1D waveform inversion constrained by well logs on a set of zero-offset vertical seismic profiles. The amplitude decay and centroid frequency shift methods of Q estimation assume that a seismic pulse propagates in a homogeneous medium and ignore the interference of the propagating wave with short-period multiples. The waveform inversion explicitly models multiple scattering and interference on a stack of thin layers using high-resolution data from sonic and density logs. This allows for stable Q estimation in small depth windows (in this study, 150 m), and separation of the frequency-dependent layer-induced scattering from intrinsic absorption. Besides, the inversion takes into account band-limited nature of seismic data, and thus, it is less dependent on the operating frequency bandwidth than on the other methods. However, all considered methods of Q estimation are unreliable in the intervals where subsurface significantly deviates from 1D geometry. At the Otway site, the attenuation estimates are distorted by sub-vertical faults close to the boreholes. Analysis of repeated vertical seismic profiles reveals that 15 kt injection of the CO₂-rich fluid into a thin saline aquifer at 1.5 km depth does not induce detectable absorption of P-waves at generated frequencies 5–150 Hz, most likely because the CO₂ plume in the monitoring well is thin, <15 m. At the Otway research site, strong attenuation Q ≈ 30–50 is observed only in shaly formations (Skull Creek Mudstone, Belfast Mudstone). Layer-induced scattering attenuation is negligible except for a few intervals, namely 500–650 m from the surface, and near the injection interval, at around 1400–1550 m, where Q_scat ≈ 50–65.     

唐山地区强震动记录的应用研究初探

本文给出了唐山地区强震动记录应用研究的两个实例,提出了建筑结构采用时程分析时选用强震动记录的原则和方法,通过对唐山地区强震动记录的分析处理,得到了其峰值加速度及加速度反应谱,确定了本地区进行弹性时程分析时选用的强震动记录;研究了局部场地条件对地震动影响的唐山响堂三维强震动观测台阵,以唐山响堂台阵2号测井(地下32m)的基岩强震动作为输入,通过2号测井的土层剖面,利用2个一维土层地震反应分析程序,分别计算得到地表的峰值加速度和加速度反应谱,并把计算结果与同次地震相应的地表强震动记录峰值加速度与加速度反应谱进行了对比分析。     

Multi‐variable relations for soil effects on seismic ground motion

Soil effects on peak ground acceleration, velocity and elastic response spectra (5% damping) are expressed by simple approximate relations in terms of five key parameters: (a) the fundamental vibration period of the non‐linear soil, T_S, (b) the period of a bedrock site of equal thickness, T_b, (c) the predominant excitation period, Te, (d) the peak seismic acceleration at outcropping bedrock, a, and (e) the number of significant excitation cycles, n. Furthermore, another relation is proposed for the estimation of T_S in terms of the soil thickness H, the average shear wave velocity of the soil V?_S,o and a. The aforementioned parameters were first identified through a simplified analytical simulation of the site excitation. The multivariable approximate relations were then formulated via a statistical analysis of relevant data from more than 700 one‐dimensional equivalent‐linear seismic ground response analyses, for actual seismic excitations and natural soil conditions. Use of these relations to back‐calculate the numerical results in the database gives an estimate of their error margin, which is found to be relatively small and unbiased. The proposed relations are also independently verified through a detailed comparison with strong motion recordings from seven well‐documented case studies: (a) two sites in the San Fernando valley during the Northridge earthquake, and (b) five different seismic events recorded at the SMART‐1 accelerometer array in Taiwan. It is deduced that the accuracy of the relations is comparable to that of the equivalent‐linear method. Hence, they can be readily used as a quick alternative for routine applications, as well as for spreadsheet computations (e.g. GIS‐aided seismic microzonation studies) where numerical methods are cumbersome to implement. Copyright © 2003 John Wiley & Sons, Ltd.     

InterPACIFIC project: Comparison of invasive and non-invasive methods for seismic site characterization. Part II: Inter-comparison between surface-wave and borehole methods

The InterPACIFIC project was aimed at assessing the reliability, resolution, and variability of geophysical methods in estimating the shear-wave velocity profile for seismic ground response analyses. Three different subsoil conditions, which can be broadly defined as soft-soil, stiff-soil, and hard-rock, were investigated. At each site, several participants performed and interpreted invasive measurements of shear wave velocity (Vs) and compression wave velocity (Vp) in the same boreholes. Additionally, participants in the project analysed a common surface-wave dataset using their preferred strategies for processing and inversion to obtain Vs profiles. The most significant difference between the invasive borehole methods and non-invasive surface wave methods is related to resolution of thin layers and abrupt contrasts, which is inherently better for invasive methods. However, similar variability is observed in the estimated invasive and non-invasive Vs profiles, underscoring the need to account for such uncertainty in site response studies. V_S,30 estimates are comparable between invasive and non-invasive methods, confirming that the higher resolution provided by invasive methods is quite irrelevant for computing this parameter.     

辽东湾近海海底土小应变动力特性试验研究

开展场地地震动反应分析是众多海洋工程活动顺利实施和长期安全运行的重要保障,其中查明海底土小应变动力特性是十分关键的。文章借助海上原位钻孔取样和室内共振柱仪对辽东湾近海海底土小应变动力特性开展试验研究,探讨前人总结提出的滨海海底土初始动剪切模量G_max预测经验公式在辽东湾海域的适用性,对比文章试验结果与周边海域已有海底土剪切波速测试结果,分析小应变范围内海底土动剪切模量的衰减特征和阻尼比发展演变趋势,并同时进行定量化模拟预测,评价海底土动剪切模量衰减曲线的归一化特征。研究结果表明：(1)共振柱试验得到的原状海底土剪切波速与周边海域已有海底土剪切波速测试结果吻合良好;(2)相比起陆地土,海洋土在小应变范围内归一化动剪切模量G/G_max衰减得更慢,阻尼比增长得也更慢;(3)通过引入临界剪应变,不同土类和埋深对应的海底土小应变动剪切模量衰减曲线可以进行归一化。研究成果对合理开展海洋工程场地地震动反应分析评价具有一定的指导意义。     

唐山响堂三维场地影响观测台阵

１９９４年７月中国地震局工程力学研究所在唐山余震区响堂镇建成了我国第一个三维场地影响观测台阵。该台阵目前有四个测点,分别布设在基岩地表、土层地表、地下１７ｍ和地下３２ｍ处。它安装了分辨率为１６位的井下数字观测系统,主机和从机同步运行,系统时钟采用Ｏｍｅｇａ导航信号自动校对（精度１ｍｓ）,地震数据采用固态方式（ＣＭＯＳ）存储,每个测点均布设一组 量力平衡式加速度计。该台阵自投入运行至１９９７年１２月     

土层平均剪切波速确定及其可靠性验证——基于钻井台阵强震动记录

震害资料显示,场地条件对地震动特性以及工程结构破坏程度影响显著。为减少因场地效应而造成的经济损失和社会影响,在进行场地地震反应分析时,需最大限度地减小因场地土层模型参数的不确定性引起的地震动评估偏差,为工程结构地震反应分析选取并生成适当的地震动输入。随着强震动观测技术的逐渐发展,大量可靠的钻井台阵记录为地震过程中场地观测点的动力反应提供了直接数据。以美国加州地区La Cienega钻井台阵强震动观测数据为基础,利用互相关函数,对不同强度地震作用下场地土层的平均剪切波速进行分析,并在此基础上,以Cyclic 1D为模拟平台,建立一维自由场地地震反应有限元分析模型。分析结果表明：通过钻井台阵地震动观测数据识别,得到场地平均剪切波速,能够反映该场地的动力特性,数值模拟计算结果和台阵地震动记录基本吻合,可为数值模型参数选取提供依据。     

Non-uniqueness in surface-wave inversion and consequences on seismic site response analyses

Surface-wave tests are based on the solution of an inverse problem for shear-wave velocity profile identification from the experimentally measured dispersion curve. The main criticisms for these testing methodologies are related to the inverse problem solution and arise from the possible equivalence of different shear-wave velocity profiles. In this paper, some implications of solution non-uniqueness for seismic response studies are investigated using both numerical simulations and experimental data. A Monte Carlo approach for the inversion problem has been used to obtain a set of equivalent shear-wave velocity models. This selection is based on a statistical test which takes into account both data uncertainty and model parameterization. This set of solutions (i.e., soil profiles) is then used to evaluate the seismic response with a conventional one-dimensional analysis. It is shown that equivalent profiles with respect to surface-wave testing are equivalent also with respect to site amplification, thus countering the criticism related to inversion uncertainty for the engineering use of surface-wave tests.     

Surface and downhole shear wave seismic methods for thick soil site investigations

Shear wave velocity–depth information is required for predicting the ground motion response to earthquakes in areas where significant soil cover exists over firm bedrock. Rather than estimating this critical parameter, it can be reliably measured using a suite of surface (non-invasive) and downhole (invasive) seismic methods. Shear wave velocities from surface measurements can be obtained using SH refraction techniques. Array lengths as large as 1000 m and depth of penetration to 250 m have been achieved in some areas. High resolution shear wave reflection techniques utilizing the common midpoint method can delineate the overburden-bedrock surface as well as reflecting boundaries within the overburden. Reflection data can also be used to obtain direct estimates of fundamental site periods from shear wave reflections without the requirement of measuring average shear wave velocity and total thickness of unconsolidated overburden above the bedrock surface. Accurate measurements of vertical shear wave velocities can be obtained using a seismic cone penetrometer in soft sediments, or with a well-locked geophone array in a borehole. Examples from thick soil sites in Canada demonstrate the type of shear wave velocity information that can be obtained with these geophysical techniques, and show how these data can be used to provide a first look at predicted ground motion response for thick soil sites.