基于萤火虫和蝙蝠群智能算法的瑞雷波频散曲线反演

反演瑞雷波频散曲线能有效获取地层横波速度和厚度.但由于其高度的非线性、多参数、多极值等特点,传统的全局搜索方法易出现收敛速度慢、早熟收敛及搜索精度低的问题.鉴于此,本文提出并测试了基于萤火虫优化算法（FA）和带惯性权重的蝙蝠优化算法（WBA）的新的瑞雷波频散曲线反演策略.在瑞雷波频散曲线反演中,FA全局搜索能力强,但后期搜索精度低,而WBA局部搜索能力强,搜索精度高,但易出现早熟收敛.故本文将二者结合,提出了一种新的优化策略,称其为WFBA,即在反演前期使用FA,后期使用WBA,很好地解决了FA后期搜索精度低及WBA早熟收敛的问题.本文首先反演了三个典型理论模型的无噪声、含噪声的数据,验证了WFBA对瑞雷波数据反演的有效性与稳定性.然后将WFBA与WBA、FA单独反演以及不含惯性权重的FBA和粒子群优化算法（PSO）反演的结果进行了对比,说明了WFBA相对于WBA、FA、FBA和PSO具有更稳定、收敛速度更快、求解精度更高等优点.最后,反演了来自美国怀俄明地区的实测资料,检验了WFBA对瑞雷波数据反演的实用性.理论模型试算和实测资料分析表明,WFBA很适用于瑞雷波频散曲线的定量解释,具有很高的实用性价值.     

基于背景噪声研究青藏高原东北缘瑞利波相速度和方位各向异性

本研究收集了甘肃、青海、宁夏等118个宽频带数字地震台站的连续波形资料,利用噪声互相关,经过计算和筛选,在5~38 s范围内,共得到5773条瑞利波相速度频散曲线.然后采用1°×1°的网格划分,反演获得青藏高原东北缘相速度和方位各向异性分布.结果表明：短周期8~12 s内,鄂尔多斯从低速异常变为高速异常;该周期范围内各向异性结果与区域断裂走向有很好的一致性.18~25 s周期内,祁连地块、松潘-甘孜地块、羌塘地块低速异常范围逐渐变大,随周期增加地壳低速异常与人工探测结果相符;鄂尔多斯表现为速度随周期增加逐渐变大,说明其中下地壳速度相对偏高,不存在低速异常;该周期范围内的各向异性特征表现为,祁连地块和松潘甘孜地块大致呈NW-SE方向,而青藏高原内部快波方向显示了顺时针旋转的形态.在30~35 s范围内面波速度主要受莫霍面深度和莫霍面附近介质速度的影响,与地壳厚度分布有非常好的吻合.综合不同方法获得的各向异性研究结果,支持印度-欧亚板块的碰撞使青藏高原东北缘地壳发生缩短和逐渐隆升的观点,认为整个岩石圈的垂直缩短变形是青藏高原东北缘的主要形成机制.     

地震槽波的数学-物理模拟初探

针对地震槽波在低速层的传播特性,开展了煤层内地震槽波勘探的数值模拟和物理模拟研究的初探工作.在数值模拟研究方面,采用交错网格有限差分法对煤层中的地震槽波进行三分量全波场模拟.基于波场快照和人工合成地震记录研究了不同模型中的波场特征和各种波型的传播规律.在物理模拟方面,通过选用不同配比的环氧树脂和硅橡胶类材料构建地震槽波物理模型,利用透射法和反射法观测系统获得了清晰的地震槽波记录以研究槽波的地震学特征.研究表明,在煤层内槽波的地震波场中,Love型槽波的能量小于Rayleigh型槽波的SV分量,大于Rayleigh型槽波的SH分量.相对于Love型槽波和Rayleigh型槽波的SH分量,Rayleigh型槽波的SV分量在围岩中的泄露能量较强.在煤层界面附近的围岩中,地震波仍以槽波形式传播,随着距离的增加能量逐渐衰减.随着煤层变薄,煤层槽波主频向高频方向移动,频散现象增强,传播速度增大.     

青藏高原东北缘及周边地区基于程函方程的面波层析成像

利用"中国地震科学台阵探测"在南北地震带北段布设的670套宽频带地震台站记录到的面波资料,使用新近发展的程函方程面波层析成像方法,获得了青藏高原东北缘及周边地区12~60 s周期范围比以往成像结果具有更高分辨率的瑞利面波相速度分布图像.青藏高原东北缘的祁连褶皱系西段、秦岭褶皱系西段和松潘一甘孜褶皱系,在16~60s周期范围内均显示出明显的低速异常分布,表明该地区的地壳力学强度较低,在强烈的构造应力作用下易发生形变.与西段不同,祁连褶皱系东段和秦岭褶皱系中段的相速度分布特征揭示,其中下地壳的速度明显高于高原内部区域.鄂尔多斯块体整体上表现为稳定块体具有的高速特征,但其西部边缘在中上地壳的速度比块体中部地区偏低,且存在一定的不均匀性.鄂尔多斯块体西北缘的临河断陷盆地和西缘的银川断陷盆地,在较短的周期范围内(12~20 s)表现为局部低速特征,但与银川断陷盆地不同,临河断陷盆地的低速特征可一直延续至60 s周期以上,表明该盆地下方地壳及上地幔速度明显偏低,可能与深部热作用有关.阿拉善块体与其北部地区的速度差异主要表现在中上地壳,这一现象值得今后进一步探讨.基于程函方程面波层析成像方法给出了青藏高原东北缘及周边地区高分辨率的成像结果,揭示了以往面波层析成像难以获得的深部细节特征,为该地区的深部构造研究提供了新的信息.     

Spatial behaviour of Rayleigh waves in layered half‐spaces under active surface sources

In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms.     

Experimental study of deformable connection consisting of buckling‐restrained brace and rubber bearings to connect floor system to lateral force resisting system

This paper presents experimental and numerical studies of a full‐scale deformable connection used to connect the floor system of the flexible gravity load resisting system to the stiff lateral force resisting system (LFRS) of an earthquake‐resistant building. The purpose of the deformable connection is to limit the earthquake‐induced horizontal inertia force transferred from the floor system to the LFRS and, thereby, to reduce the horizontal floor accelerations and the forces in the LFRS. The deformable connection that was studied consists of a buckling‐restrained brace (BRB) and steel‐reinforced laminated low‐damping rubber bearings (RB). The test results show that the force–deformation responses of the connection are stable, and the dynamic force responses are larger than the quasi‐static force responses. The BRB+RB force–deformation response depends mainly on the BRB response. A detailed discussion of the BRB experimental force–deformation response is presented. The experimental results show that the maximum plastic deformation range controls the isotropic hardening of the BRB. The hardened BRB force–deformation responses are used to calculate the overstrength adjustment factors. Details and limitations of a validated, accurate model for the connection force–deformation response are presented. Numerical simulation results for a 12‐story reinforced concrete wall building with deformable connections show the effects of including the RB in the deformable connection and the effect of modeling the BRB isotropic hardening on the building seismic response. Copyright © 2016 John Wiley & Sons, Ltd.     

RESEARCH ON PASSIVE SERVO SEISMIC ROTATIONAL ACCELERATION SENSOR WITH LARGE DAMPING

Many strong motion records show that under the strong seismic vibration of, the torsional disfigurement of building structures is a common and serious damage. At present, there are no special sensors for measuring seismic rotation in the world. Most of the experts obtain rotational components through observing deformation, theoretical analysis and calculation. The theory of elastic wave and source dynamics also prove the conclusion that the surface of the earth will rotate when an earthquake occurs. Based on a large number of investigations and experiments, a rotational acceleration sensor was developed for the observation of the rotational component of strong ground motions. This acceleration sensor is a double-pendulum passive servo large-damped seismic rotational acceleration sensor with the moving coil transducer. When an earthquake occurs, the seismic rotational acceleration acts on the bottom plate at the same time. The magnetic circuit system and the middle shaft fixedly connected to the bottom plate follow the bottom plate synchronous vibration, and the moving part composed of the mass ring, the swing frame and the moving ring produces relative corners to the central axis. The two working coils mounted on the two pendulums produce the same relative motion with respect to the magnetic gaps of the two magnetic circuits. Both working coils at this time generate an induced electromotive force by cutting magnetic lines of force in the respective magnetic gaps. The generated electromotive forces are respectively input to respective passive servo large damper dynamic ring transducer circuits and angular acceleration adjusting circuits, and the signals are simultaneously input to the synthesizing circuit after conditioning. Finally, the composite circuit outputs a voltage signal proportional to the seismic rotational acceleration to form a seismic rotational acceleration sensor. The paper presents the basic principles of the rotational acceleration sensor, including its mechanical structure diagram, circuit schematic diagram and mathematical models. The differential equation of motion and its circuit equation are derived to obtain the expressions of the main technical specifications, such as the damping ratio and sensitivity. The calculation shows that when the damping ratio is much larger than 1, the output voltage of the passive servo large damping dynamic coil transducer circuit is proportional to the ground rotation acceleration, and the frequency characteristic of bandpass is wider when the damping ratio is larger. Based on the calibration test, the dynamic range is greater than or equal to 100dB and the linearity error is less than 0.05%. The amplitude-frequency characteristics, the phase-frequency characteristics and their corresponding curves of the passive servo rotational acceleration sensor are acquired through the calculations. Based on the accurate measurement of the micro-vibration of the precision rotating vibration equipment, the desired result is obtained. The measured data are presented in the paper, which verify the correctness of the calculation result. The passive servo large damping rotational acceleration sensor has simple circuit design, convenient operation and high resolution, and can be widely applied to seismic acceleration measurement of earthquake or structure.     

利用随机噪声研究华南地区瑞利波相速度特征

基于2009~2010年国家、区域测震台网的387个宽频带台站数据,通过互相关方法提取到可靠的瑞利波经验格林函数,利用相位匹配滤波时频分析技术测量瑞利波相速度频散曲线,最后采用噪声层析成像方法获得了华南地区不同周期的瑞利波相速度分布图。结果显示,华南地区速度结构横向变化幅度较小,反映了华南地区作为一个整体较为稳定,与华南地区自晚中生代以来未发生过强烈构造活动的特征基本一致;虽然华南地区整个岩石圈速度结构较为均匀,但扬子块体西部、四川盆地与扬子块体东部、华夏块体间存在明显的速度差异,体现在周期为8~10s时华夏块体相速度大于扬子块体西部、川滇块体以及四川盆地,由于沉积层较厚四川盆地速度最低;周期为10~30s时华夏块体面波相速度大于扬子块体西部和四川盆地,而川滇块体速度最低;周期为35s时扬子块体、华夏块体、四川盆地的速度基本一致,且高于川滇块体,这与华南地区地壳厚度明显小于川滇块体的特征相符。     

Evaluation of wind effects on a supertall building based on full‐scale measurements

This paper describes the results obtained from the full‐scale measurements of wind effects on a 70‐storey building in Hong Kong. The building which has a height of approximately 370 m is the second tallest structure in Hong Kong. The field data such as wind speed, wind direction and wind‐induced acceleration responses have been measured since 1995 including the close passage of two typhoons; typhoon Sally and typhoon Kent. Detailed analysis of the field data is conducted. The full‐scale measurements are compared with the wind tunnel results obtained in the Boundary Layer Wind Tunnel Laboratory at Western Ontario University. The amplitude‐dependent characteristics of damping and natural frequency that were obtained by using the random decrement technique are investigated. Copyright © 2000 John Wiley & Sons, Ltd.