基于稀疏反演理论的自动叠加速度反演方法

叠加速度分析技术是常规地震资料处理中的重要环节,也是经典的时间域速度建模方法.叠加速度分析技术主要包括速度谱计算和拾取两个步骤.至今为止,多数研究工作通过提高速度谱的分辨率以及抗噪声能力,获得高质量的速度谱从而有利于拾取.本文的目标是将叠加速度分析技术转为一个全自动化的处理流程.从参数估计的角度出发,将叠加速度估计转化为稀疏反演框架下的模型参数估计问题,并通过稀疏反演算法自动反演叠加速度,进而提高叠加速度建模的效率.为实现这一目标,首先给出了正问题的定义,即层状介质中CMP道集的预测模型,利用叠加速度、垂向双程走时(t_0)以及反射子波以及CMP道集时距关系(如双曲时距关系)可以预测CMP道集.接着,速度分析反问题可以描述为已知观测的CMP道集,估计模型参数(叠加速度及t_0时间等).利用模型参数的稀疏性作为约束条件并用L_0范数作为模型稀疏性的度量准则,叠加速度分析可以转化为L_0范数约束下的稀疏反演问题.本文提出了一种基于预测校正思想的匹配追踪算法求解上述反问题,实现了自动叠加速度建模并为后续的高精度速度反演方法提供较好的初始模型.理论和实际资料的测试结果证明了本文方法的有效性.     

ANALYSIS ON THE CHARACTERISTICS OF DURATION AND PERIOD OF GROUND MOTION OF THE LUSHAN EARTHQUAKE BASED ON THE STATION RECORDS

Based on 60 records from the 20 stations within 100km to the epicenter of Lushan earthquake, the predominant period, period of peak response spectrum, duration of ground motion, and source duration are investigated. By the study, we conclude that within 100km to the epicenter, the scope of predominant period is 0.013~0. 275s in EW, NS and UD direction; the scope of period of peak response spectrum for 5% damping ratio is 0.03~0.65s; the scope of 90% ground motion durations is 5. 1~35. 9s; the scope of averaging source duration is 6.41 (EW), 6.05 (NS) and 5.47s (UD). Furthermore, based on the ground motion duration calculated by 20 stations dada, the predictive curve and equation of ground motion duration is regressed and compared with the recent equation of ground motion duration by Bommer (2009). We find that the ground motion durations of most stations are larger than predictively mean value by Bommer (2009), which means that the source effect of Lushan earthquake is relatively larger. Lastly, by the contour figure of ground motion duration and source duration, we conclude that the directivity character of duration is relatively apparent in NE direction. The relatively larger source duration and ground motion duration in NE direction indicate more energy release in the main shock of Lushan earthquake, which perhaps causes the relatively less aftershocks in this direction. The duration has no hanging wall effect, which perhaps results from the blind-reverse fault structure of Lushan earthquake.     

全张量磁梯度数据解释的均衡边界识别及深度成像技术

全张量磁梯度数据具有高精度、高分辨率、多参量的优点,能更加清晰地刻画地质体的分布特征,综合利用磁张量梯度数据准确地获得地质体水平位置和深度信息是解释的主要目的.磁张量数据的方向解析信号具有减小倾斜磁化干扰的优点,常被用来圈定磁源体的水平位置,但解析信号强度随着地质体埋深的增加急剧衰减,难以有效识别较深的地质体.张量数据均衡边界识别技术,利用不同方向解析信号的比值函数,能有效地均衡不同深度地质体的响应,同时显示不同深度地质体的边界,提高了对较深地质体的分辨率.磁张量数据深度成像技术根据实测张量数据与假定模型张量数据的相关系数来给定地质体的深度,综合利用多参量数据联合反演提高了反演结果的准确性,且无需进行复杂的反演运算,是大数据量张量数据解释的有效方法.理论模型试验证明：磁张量数据均衡边界识别技术可清晰和准确地识别地质体的水平范围,受倾斜磁化干扰小;磁张量数据深度成像技术可准确地获得地质体的深度信息,具有较强的抗噪性.将上述方法应用于铁矿区实测航磁张量梯度数据解释,获得了铁矿体水平分布与埋深,深度结果与张量欧拉反褶积法计算结果一致.     

红山基准地震台ELF电磁观测资料异常特征分析

对河北省红山基准地震台2014年3月7日至2015年12月31日ELF观测的天然电磁场数据进行处理,研究电磁场强度与视电阻率变化规律,分析2014年8月3日鲁甸6.5级地震、2015年的4月19日文安3.0级地震、4月25日尼泊尔8.1级地震、11月15日滦县2.9级地震发生前后电磁短临异常现象,结果表明,1 Hz、39 Hz、178 Hz电磁场自功率谱有明显的台阶、突跳和集丛现象,视电阻率异常变化幅度达到1—2个数量级。     

一种深厚覆盖层场地地震动参数的确定方法

基于欧美规范确定了坐落在深厚覆盖层上KH抽水蓄能电站上、下库场地基本运行和最大设计地震动峰值加速度、反应谱和时程等动参数。首先依据场地区域地震烈度区划图、特征周期区划图和依据场地地质地震条件选取的5条种子实测地震动确定场地基岩输入加速度时程、峰值加速度和设计反应谱,进而基于各土层地质参数和一维弹性波传播模拟程序确定覆盖层表面的平均峰值加速度、平均反应谱和5条地震动时程,对所得到的平均反应谱和峰值加速度进行光滑处理后确定可用于各建筑物结构抗震设计的地震动参数,包括覆盖层表面水平向动力响应加速度时程、峰值加速度和设计反应谱。该方法可较好地保留输入地震动的真实动力特性,如持时、相位和频率等,为我国规范中建议的确定场地地震动参数的方法提供有益的补充。     

岷县漳县6.6级地震及余震强震动记录特征分析

对2013年7月22日甘肃岷县—漳县地震获取的273条主余震加速度记录进行格式转换、基线校正和滤波等常规处理,分析该地震主震(M_S6.6)和余震(M_S5.6)两次地震记录的幅值、持时以及反应谱特征,发现M_S6.6主震记录的PGA范围在0.728～177.5 gal间,M_S5.6余震记录的PGA范围在0.732～69.3 gal间;将观测数据与霍俊荣和第五代《中国地震动参数区划图》地震动衰减关系进行对比,发现霍俊荣衰减关系更吻合于此次地震的主余震加速度衰减;绘制主余震5%～95%重要持时分布图,并针对主震62MXT反应谱和本地设计谱以及近些年国内主要强震震中反应谱开展比较分析;最后研究土层台、基岩台、相同台站各个震级反应谱的特征。     

基于FLAC3D路桥不稳定边坡治理工程的动力响应研究

抗滑桩广泛应用于边坡支护工程,但路桥不稳定边坡抗滑桩对其进行加固时,其在地震作用下的动力响应问题仍未定论。以道真至瓮安高速公路河槽大桥不稳定边坡治理工程为依托,以SIMQKE_GR软件拟合地震设计反应谱,利用FLAC3D有限差分软件建立数值模型,对比地震作用下抗滑桩等支挡结构设立与否对路桥不稳定边坡的动力响应变化规律的影响。结果表明:（1）加速度放大系数无论是沿坡内还是抗滑桩,其值沿竖向均存在一定的震荡变化,但在抗滑桩支护情况下,其震荡幅度明显变小,加速度放大系数变化更平稳。（2）速度放大系数随着高度的增加逐渐减小,未加抗滑桩时速度放大系数在顶层附近出现大幅度波动,添加了抗滑桩的情况下其震荡幅度明显趋缓,变化更稳定。（3）坡体在施加抗滑桩后,治理效果明显,残余位移最大值由无抗滑桩的16.9 cm变为有抗滑桩的1 cm,但最大值均位于边坡表层碎石土和粉质黏土附近。     

俯冲带板间地震水平加速度谱阻尼修正系数研究

基于日本K-net和Kik-net台站的俯冲带地区板间地震的数据,分SCⅠ~SCⅣ四种场地分别建立了水平向加速度反应谱阻尼修正系数模型。利用最小二乘法进行回归拟合,拟合公式采用对数三次项提高拟合精度。通过统计检验和残差分析,验证模型的拟合效果。研究表明:低于5%与高于5%阻尼比的阻尼修正系数变化趋势差异明显,分别在周期点0.3~0.5s区间达到极大值和极小值;不同场地类别的阻尼修正系数统计上显著不同,建立不同场地类别的模型是必要的;模型残差的标准差分布在0.012~0.712之间,建立的阻尼修正系数模型较好的拟合了阻尼修正系数随谱周期和阻尼比的变化趋势。     

2018年2月台湾花莲震群热红外遥感亮温异常分析

应用中国静止气象卫星FY-2G的TBB产品为基础数据,针对2018年2月上旬台湾花莲附近海域震群前出现在震中及附近区域的显著热红外异常,运用小波变换和相对功率谱法进行异常提取和资料分析。结果表明：花莲震群热红外异常范围主要分布在震中SW区域,震前约8个月出现功率谱增强现象,随时间推移,异常范围延纵谷断裂带逐渐扩大,持续近2个月,最大异常幅度超过14。     

地动噪声RMS值计算方法的适用性及误差概率分析

对比两种计算RMS值的方法,结合误差概率统计给出结果的误差概率分布。在使用有限脉冲响应（FIR）带通滤波器进行滤波并计算RMS值的过程中,滤波器窗口函数和阶数是影响计算结果误差的主要因素,其中阶数的作用更大。通过计算功率谱密度（PSD）反算RMS值的误差主要受pwelch函数的窗口函数类型、窗口长度及重叠率等参数影响,其中窗口长度作用更大。从同等误差水平的概率分布看,在合理设置参数的前提下,使用PSD反算RMS值的方法更优。