不同支座布置下组合基础隔震结构能量反应分析

通过对5种不同支座布置方案的组合基础隔震结构施加不同场地条件下的双向水平地震波,研究了场地类别和地震波加速度峰值及支座布置方式对结构总输入能量的影响,并进行了耗能分析。结果表明,随着地震波峰值的增大,结构总输入能量增大;相同地震峰值条件下,场地土越软,总的输入能量越大;组合隔震结构应用于软土场地时隔震层滞回耗能比更大,隔震效果更好;组合隔震结构两种支座中选择分散布置方案时,隔震层耗能比更大,隔震效果优于支座集中布置方案。     

2017年6月16日湖北秭归Ms4.3地震成因初探

北京时间2017年6月16日19时48分, 湖北省秭归县发生M_S4.3地震(震中位置31.06°N, 110.48°E)。 此次地震导致当地160多间建筑物出现不同程度损坏, 100多人生活受到影响, 震区中危岩体出现一定程度破坏。 分析此次地震的震源参数将有助于了解其震源特征和发震机制。 本文首先获取了震中距在3°以内的40个宽频带固定台站的三分量数据, 然后利用FK(Frequency wavenumber)方法计算了在不同深度下Crust2.0模型和改进的1-D模型的格林函数, 分别使用Cut And Paste(CAP)方法反演得到本次地震的震源机制解和深度; 并采用基阶Rayleigh波振幅谱进一步约束了质心深度。 结果表明此次地震的最佳震源机制解的节面Ⅰ: 68°/59°/163°; 节面Ⅱ: 166°/75°/32°, 最佳矩震级为M_W4.3, 最佳深度约为5 km。 对比2013年巴东地震、 2014年秭归地震以及2011年湖北阳新—江西瑞昌地震的震源参数和余震信息, 此次地震有可能是水库诱发地震, 但要判断成因还需进一步研究。     

Some remarks on Q-compensated sparse deconvolution without knowing the quality factor Q

The subsurface media are not perfectly elastic, thus anelastic absorption, attenuation and dispersion (aka Q filtering) effects occur during wave propagation, diminishing seismic resolution. Compensating for anelastic effects is imperative for resolution enhancement. Q values are required for most of conventional Q-compensation methods, and the source wavelet is additionally required for some of them. Based on the previous work of non-stationary sparse reflectivity inversion, we evaluate a series of methods for Q-compensation with/without knowing Q and with/without knowing wavelet. We demonstrate that if Q-compensation takes the wavelet into account, it generates better results for the severely attenuated components, benefiting from the sparsity promotion. We then evaluate a two-phase Q-compensation method in the frequency domain to eliminate Q requirement. In phase 1, the observed seismogram is disintegrated into the least number of Q-filtered wavelets chosen from a dictionary by optimizing a basis pursuit denoising problem, where the dictionary is composed of the known wavelet with different propagation times, each filtered with a range of possible values. The elements of the dictionary are weighted by the infinity norm of the corresponding column and further preconditioned to provide wavelets of different values and different propagation times equal probability to entry into the solution space. In phase 2, we derive analytic solutions for estimates of reflectivity and Q and solve an over-determined equation to obtain the final reflectivity series and Q values, where both the amplitude and phase information are utilized to estimate the Q values. The evaluated inversion-based Q estimation method handles the wave-interference effects better than conventional spectral-ratio-based methods. For Q-compensation, we investigate why sparsity promoting does matter. Numerical and field data experiments indicate the feasibility of the evaluated method of Q-compensation without knowing Q but with wavelet given.     

南海海底地震仪异常数据的分析和处理

海底地震仪(Ocean Bottom Seismometer, OBS)数据处理至关重要,是获取深部地壳结构的基础与前提.2006年实施OBS2006-2测线时,有2台OBS(OBS03,OBS06)数据出现异常,无法使用.由于海上航次花费巨大,采集到的数据弥足珍贵.本文采用数据格式检查、邻近台站对比分析、重采样等方法,成功地对这2台OBS数据进行了解编处理,得到了这两个台站的综合地震记录剖面;利用上述方法对2011年实施的OBS973-3测线中的异常台站OBS03进行了分析处理,同样得到了OBS03台站的综合地震剖面;通过查看两次海上实验班报发现,OBS2006-2测线之OBS06与OBS973-3测线之OBS03内部Sedis编号相同,为同一台记录仪器,再一次验证上述处理方法正确可行;然后对OBS2006-2测线2个台站进行震相识别与走时拾取后,利用前人纵波速度模型开展了射线追踪与走时模拟.此次对异常OBS数据的重新处理工作,不仅为OBS探测提供了宝贵的数据处理经验,而且将提高OBS2006-2测线地壳结构的可靠性和约束性,具有重要的研究意义.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

青藏高原及周边2001年以来三次特大地震引起的形变场分布特征

基于地壳黏弹模型在GPS观测资料和地震位错数据为约束条件下,应用有限元数值模拟方法对2001年昆仑山8.1级、2008年汶川8.0级和2015年尼泊尔8.1级特大地震引起的地壳形变分布特征进行了模拟计算,获得了地震位移场和形变场.这3次大地震分别发生在青藏高原的北部、东部边界和南部边界,尽管震级大小基本相当,但发震区域和发震断层性质都各不相同,其各自产生的地壳形变场分布特征存在明显的差异,主要表现为形变场区域大小、幅值大小等的差异,以及在不同地壳深度也存在明显的差异,这些差异主要取决于地震断层性质和地下介质结构的不同.通过模拟计算,可以进一步了解大地震产生的应力加卸载区分布特征,对预测未来地震发生区域范围提供重要参考依据.     

甘孜-玉树断裂带晚第四纪走滑速率与滑动分解作用

滑动速率是研究断裂运动学特征、地震活动性和区域应变分配的重要参数和依据。前人关于甘孜-玉树断裂带滑动速率的研究结果存在较大差异,因此,其晚第四纪滑动速率有待进一步调查研究。本文基于卫星影像解译和野外实地考察,对甘孜-玉树断裂带西段（玉树断裂）上典型断错地貌点进行测量分析,得到玉树断裂晚第四纪走滑速率为6.6±0.1-7.4±1.2mm/a。通过与前人对甘孜-玉树断裂带东段（甘孜断裂）滑动速率的研究结果进行对比,发现甘孜-玉树断裂带东、西段滑动速率不一致,其原因是甘孜断裂的左旋滑移在向西传递的过程中,一部分应变被分配到了巴塘盆地南缘断裂上。巴塘盆地南缘断裂的存在很好地解释了玉树断裂的走滑速率比甘孜断裂偏低的原因。但是,从区域变形来看,巴塘盆地南缘断裂分配的滑动速率恰好说明了甘孜-玉树断裂带东、西段及鲜水河断裂带的水平构造变形是协调一致的。     

长岛震群震源谱参数反演和震源参数特征

采用高频截止(High-Cut)震源模型,以均方根误差最小原则稳健地求解震源谱参数,并由此推算震源尺度和静态应力降。实际应用显示,该模型的理论谱对观测谱有很好的拟合,可明显改善拐角频率识别准确度。计算了长岛震群内71次M_L≥2.5地震事件的震源参数,结果表明：(1)拐角频率处于2～10 Hz范围,与震级大小存在一定的相关性,截止频率范围处于10～30 Hz之间,与地震大小的相关性不明显;(2)地震矩M₀分布在10¹²～10¹⁴ N·m,与震级M_L存在正相关关系：logM₀=0.977M_L+10.186;M_L与矩震级M_W之间的关系为：M_W=0.651M_L+0.766;(3)根据相对应力降时域演化发现,自2017年3月3日M_L4.5地震之后应力快速释放,应力降水平在均值附近波动,而且多数M_L≥3.5地震发生于...