2017年8月8日四川九寨沟7.0级地震强震记录及特征分析

2017年8月8日四川九寨沟县发生7.0级地震,中国数字强震动台网布设在四川、甘肃、陕西、宁夏的66个强震台获得主震加速度记录。本文首先对198条三分向强震记录进行常规处理,计算出近场强震记录的加速度峰值随震中距的分布情况;再根据2个典型台站的加速度时程记录,通过计算其加速度反应谱并与设计反应谱比较,分析本次地震的基本特征;然后将实际观测数据与意大利新一代地震动衰减公式对比,分析峰值加速度（PGA）及谱加速度的衰减关系;最后结合已有的工程场地钻探资料,采用H/V谱比法对2个不同类别的典型台站进行地场地效应分析,发现该方法能很好地反映实际台站场地的反应特征。     

爆炸波土层动力响应对比分析

以人工爆炸波为震源,通过现场测试获取基岩及土层的地震动参数,并采用等效线性化分析方法计算相应的地震动参数用于与实测结果进行对比分析。峰值加速度对比结果表明,等效线性化分析方法对于Ⅱ类场地的适应性较好,计算结果与实测结果非常接近,而Ⅲ类场地的计算结果与实测值之间存在较大的误差。加速度反应谱的对比结果表明,无论是计算结果还是实测记录,加速度反应谱的峰值均比基岩输入的要大,且土层反应计算的结果小于实际记录;加速度反应谱的宽度与场地类别关系密切,Ⅲ类场地明显比Ⅱ类场地要大,两类场地的计算结果也均小于实测值。     

地震动峰值加速度衰减关系对中国西部地区的适用性分析

地震动参数衰减关系描述了地震动参数随震级和距离等因素的变化情况,选取合适的地震动衰减关系是地震危险性分析中确定工程场地地震动参数的关键环节。本文收集、整理和分析了中国西部地区2012年3月9日—2017年9月16日间发生的42次地震事件中获取的强震动记录资料,并选择了3个国际上基于不同地区强震动记录数据建立的地震动加速度衰减关系(ASB11、SCEMY97和LLCS11),开展了加速度衰减关系计算结果与中国西部强震动记录值的对比研究。比较了衰减关系峰值加速度的预测中值与实际记录值之间的差异,并分析了其残差随震级和距离的变化。得到了以下结论:①相比于SCEMY97和LLCS11,ASB11更适合于中国西部地区;②研究中国西部地区的衰减关系时,有必要考虑高频地震动的震级饱和现象。同时,对于如何修正已有衰减关系以使其更适合于中国西部地区给出了相应的建议。     

Comparing the performances of four stochastic optimisation methods using analytic objective functions, 1D elastic full‐waveform inversion,and residual static computation

We compare the performances of four stochastic optimisation methods using four analytic objective functions and two highly non‐linear geophysical optimisation problems: one‐dimensional elastic full‐waveform inversion and residual static computation. The four methods we consider, namely, adaptive simulated annealing, genetic algorithm, neighbourhood algorithm, and particle swarm optimisation, are frequently employed for solving geophysical inverse problems. Because geophysical optimisations typically involve many unknown model parameters, we are particularly interested in comparing the performances of these stochastic methods as the number of unknown parameters increases. The four analytic functions we choose simulate common types of objective functions encountered in solving geophysical optimisations: a convex function, two multi‐minima functions that differ in the distribution of minima, and a nearly flat function. Similar to the analytic tests, the two seismic optimisation problems we analyse are characterised by very different objective functions. The first problem is a one‐dimensional elastic full‐waveform inversion, which is strongly ill‐conditioned and exhibits a nearly flat objective function, with a valley of minima extended along the density direction. The second problem is the residual static computation, which is characterised by a multi‐minima objective function produced by the so‐called cycle‐skipping phenomenon. According to the tests on the analytic functions and on the seismic data, genetic algorithm generally displays the best scaling with the number of parameters. It encounters problems only in the case of irregular distribution of minima, that is, when the global minimum is at the border of the search space and a number of important local minima are distant from the global minimum. The adaptive simulated annealing method is often the best‐performing method for low‐dimensional model spaces, but its performance worsens as the number of unknowns increases. The particle swarm optimisation is effective in finding the global minimum in the case of low‐dimensional model spaces with few local minima or in the case of a narrow flat valley. Finally, the neighbourhood algorithm method is competitive with the other methods only for low‐dimensional model spaces; its performance sensibly worsens in the case of multi‐minima objective functions.     

Seismic analysis of a tall metal wind turbine support tower with realistic geometric imperfections

The global growth in wind energy suggests that wind farms will increasingly be deployed in seismically active regions, with large arrays of similarly designed structures potentially at risk of simultaneous failure under a major earthquake. Wind turbine support towers are often constructed as thin‐walled metal shell structures, well known for their imperfection sensitivity, and are susceptible to sudden buckling failure under compressive axial loading. This study presents a comprehensive analysis of the seismic response of a 1.5‐MW wind turbine steel support tower modelled as a near‐cylindrical shell structure with realistic axisymmetric weld depression imperfections. A selection of 20 representative earthquake ground motion records, 10 ‘near‐fault’ and 10 ‘far‐field’, was applied and the aggregate seismic response explored using lateral drifts and total plastic energy dissipation during the earthquake as structural demand parameters. The tower was found to exhibit high stiffness, although global collapse may occur soon after the elastic limit is exceeded through the development of a highly unstable plastic hinge under seismic excitations. Realistic imperfections were found to have a significant effect on the intensities of ground accelerations at which damage initiates and on the failure location, but only a small effect on the vibration properties and the response prior to damage. Including vertical accelerations similarly had a limited effect on the elastic response, but potentially shifts the location of the plastic hinge to a more slender and, therefore, weaker part of the tower. The aggregate response was found to be significantly more damaging under near‐fault earthquakes with pulse‐like effects and large vertical accelerations than far‐field earthquakes without these aspects. Copyright © 2016 John Wiley & Sons, Ltd.     

Median floor acceleration spectra of linear structures with uncertain properties

Response parameters used to estimate nonstructural damage differ depending on whether deformation‐sensitive or acceleration‐sensitive components are considered. In the latter case, seismic demand is usually represented through floor spectra, that is response spectra in terms of pseudo‐acceleration, which are calculated at the floor levels of the structure where the nonstructural components are attached to. Objective of this paper is to present a new spectrum‐to‐spectrum method for calculating floor acceleration spectra, which is able to explicitly account for epistemic uncertainties in the modal properties of the supporting structure. By using this method, effects on the spectra of possible variations from nominal values of the periods of vibration of the structure can be estimated. The method derives from the extension of closed‐form equations recently proposed by the authors to predict uniform hazard floor acceleration spectra. These equations are built to rigorously account for the input ground motion uncertainty, that is the record‐to‐record variability of the nonstructural response. In order to evaluate the proposed method, comparisons with exact spectra obtained from a standard probabilistic seismic demand analysis, as well as spectra calculated using the Eurocode 8 equation, are finally shown. Copyright © 2017 John Wiley & Sons, Ltd.     

谱强度SI和地震动峰值加速度PGA参数在高铁列车地震预警中的应用

为了降低没有破坏性的高频小震对高铁地震报警的干扰,引入谱强度SI参数作为报警参数,分析其在高铁地震报警中的适用性。研究发现低频振动对SI的影响比高频振动大,对于高铁线路附近的高频小震,SI参数能有效的排除,但对于破坏性不大的远震大震可能会引发报警,从而产生误报现象,影响高铁列车运行效率。为了降低SI在破坏性小的远震大震中的误报率,本文引入PGA与SI作为报警参数。最后给出了不同时速下基于PGA与SI联合报警的阈值供参考。     

近钻头井斜动态测量重力加速度信号提取方法研究

利用加速度计进行近钻头井斜动态测量时,钻具的高速旋转,井下强振动、强冲击环境给重力加速度测量带来极大干扰,如何从干扰噪声中有效提取重力加速度信号对于提高井斜角和工具面角的测量精度至关重要.根据重力加速度径向和切向分量为周期性信号,轴向分量为近似直流信号,离心加速度为缓慢变化信号,振动和冲击加速度为随机信号的特征,提出一种基于互相关检测的重力加速度信号提取方法,选择径向或切向磁力计信号作为参考信号.利用相关检测方法分别对仿真数据、实验室旋转测试数据和振动测试数据进行处理,并计算了井斜角和工具面角,结果表明相关检测方法可以有效提取重力加速度信号,基于MEMS加速度计的井斜角测量精度优于0.5°,工具面角变化连续平滑,满足井斜动态测量需求.通过进一步开展实际钻井环境测试,该方法有望应用于近钻头地质导向系统或旋转导向钻井系统等需要近钻头井斜动态测量的场合.     

山东省场地类别分布及地震动峰值加速度区划调整

首先,本文搜集了山东省5220个建设工程场地的资料,引入"分布区"的概念,统计分析了山东省场地类别的区域分布并进行了分区。然后,根据《中国地震动参数区划图（GB18306-2015）》中的山东省地震动峰值加速度区划以及Ⅲ类场地调整系数,对山东省Ⅲ类场地分布区的地震动峰值加速度区划进行了调整,得到了考虑场地类别分区的山东省地震动峰值加速度区划图,增进了对山东省地震灾害风险空间分布的认识。结果表明,山东省场地类别分布与地貌分区有很强的相关性,Ⅰ-Ⅱ类场地分布区包括鲁中南山地丘陵区（西南部的山间平原地带除外）以及鲁东丘陵区,Ⅲ类场地分布区包括鲁西北-鲁西南平原区和鲁中南山地丘陵区西南部山间平原地带。Ⅰ-Ⅱ类场地分布区约占全省陆地面积的59.5%,Ⅰ-Ⅱ类场地分布区与Ⅲ类场地分布区的面积比例约为1.47。地震动峰值加速度区划图经过调整后,Ⅶ度及以上设防区域占全省陆地面积的比例由79%提高到90%,10个地级市的峰值加速度有提高。     

Displacement–acceleration control via stiffness–damping collaboration

An approach is presented to stiffness–damping simultaneous optimization for displacement–acceleration simultaneous control. To make a shear building model stiffer, the sum of mean-square interstorey drifts to stationary random excitations is minimized or the mean-square top-floor absolute acceleration is maximized subject to the constraints on total storey stiffness capacity and total damper capacity. Optimality conditions are derived and a two-step optimization method using the optimality conditions is devised. In the first step, the optimal design is found for a specified set of total storey stiffness capacity and total damper capacity. In the second step, a series of optimal designs is found with respect to a varied set of total storey stiffness capacity and total damper capacity. While increase of total stiffness capacity and increase of total damper capacity are both effective in reduction of deformation, only increase of total damper capacity is effective in reduction of acceleration. Acceleration control is carried out in the second step via increase of total damper capacity. It is shown through numerical examples that the proposed method is efficient and reliable. Copyright © 1999 John Wiley & Sons, Ltd.