基于快速流式算法的局部余弦相似属性

相似性是一种常用的衡量不同图像之间差异程度的属性,广泛应用于地震数据处理环节.由于地震数据本质上是非平稳的,局部相似性比全局相似性更适用于刻画地震数据的时-空变化特征.现有的局部相似属性可以通过正则化最小二乘问题进行计算,但是其计算过程需要大量的计算时间和数据存储空间,难以适应当前的海量数据处理任务.本文提出了一种基于...     

巴颜喀拉块体东缘形变及九寨沟地震孕震环境数值分析

2017年8月8日四川发生九寨沟M7.0地震,是继2008年汶川M8.0地震后发生在巴颜喀拉块体东部的又一强震.现今GPS速度观测数据显示,2008年汶川地震前后的1999-2007年和2011-2016年两个时间段内巴颜喀拉块体东部地表速度场存在明显的差异.本文以实际GPS速度观测资料为约束,构建三维有限元地球动力学模型,分别计算分析了两个时段内震源区及周边现今地壳形变、弹性应变能和应力积累特征,进一步探讨汶川地震的发生对九寨沟地区变形及应力的影响.数值模拟结果显示,汶川地震之后(2011-2016年)巴颜喀拉块体东部的地壳形变、弹性应变能积累及应力积累速率均明显大于震前,增加量值达1.5-3倍;九寨沟地震发震断裂上库仑应力增长率在1999-2007年约为0.7 kPa·a~(-1),2011-2016年间增至1.2 kPa·a~(-1).上述结果表明,现今巴颜喀拉块体东部地壳应力积累过程有利于左旋走滑型九寨沟地震的发生,汶川地震的发生调整了区域应力状态,加速了九寨沟地震的孕育过程.     

转换波波动方程叠前深度偏移

克希霍夫叠前时间偏移被广泛运用于转换波地震成像领域,该方法具有计算效率高、速度模型易于获取的优点。但是,纵横波速度差使纵波和转换波的成像在时间域不匹配,且克希霍夫偏移也无法较好地处理多次到达和保幅等问题。这里提出一种基于波动方程的转换波叠前深度偏移方法,以改善克希霍夫时间偏移的不足。     

近断层脉冲型地震动重要参数的识别方法

在拟合或合成近断层脉冲型地震动的研究中,速度等效脉冲中的参数需要依据实际地震动确定,因此,如何从地震动中识别这些参数就成为一个需要解决的问题.详细分析了已有脉冲周期和脉冲速度峰值的识别方法,分析结果表明这些方法均具有局限性,尤其是对于不规则地震动.提出了基于平滑后的地震动识别脉冲周期和脉冲速度峰值的方法,并验证了这种识别方法的有效性.结果表明: 提出的脉冲周期和脉冲速度峰值识别方法可以克服以往方法的局限性,可以实现脉冲周期和脉冲速度峰值识别的程序化.     

响应谱地震动预测方程中场地效应参数Vs30与场地周期的比较

最近的许多地震动预测方程（GMPE）都使用上30m土层的平均剪切波速度Vs30表示场地效应。然而,在一些研究发现V鹦。是表征场地效应合理参数的同时,另一些研究则给出了相反的证据。本研究使用日本的大地震动数据集对这两个场地效应参数的预测能力进行系统的比较。采用该方法的基础是使用场地周期（Ts,4倍剪切波从基岩到地表的走时）或Vs30对经验模拟场地效应的标准偏差和放大比振幅进行比较。模拟的场地效应除了包括地震动预测方程的场地效应项,还特别包括KiK-net台网的地表与井下记录之间的场地放大比。就KiK—net台网数据而言,对Ts〉0．6s的土层场地,Ts被确定是比Vs。。好的预测参数;而对T。〈0．6s的场地,这两个参数对放大比产生了相似的变异性。在所有场地类别的多数周期上,由地震动预测方程获得的场地效应,Vs30和Ts在统计上相等;而在一些谱周期上,Vs30产生了比Ts小的变异性。KiK-net台网地表一井下记录与地震动预测方程结果的矛盾可能由地震动预测方程中包括震源变异性、路径变异性和场地变异性的大变异性所致。相比之下,KiK-net台网地表一井下数据对的变异性小多了。虽然VS30。和Ts对地震动预测方程的数据产生了统计上相似的标准偏差,但瓦仍然获得了比Vs30更好的中值放大比。     

浅层结构不确定性对广角地震走时正演模拟结果的影响——以西沙地块OBS2011-1测线为例

海洋地球物理勘探中常使用多道地震数据建立浅部地层模型,结合海底地震仪广角地震数据进行射线追踪和走时反演求取地层速度特征.受采集环境和成本限制,部分海底地震仪测线缺少同步采集的多道地震数据,建立的初始模型其浅层结构有较大不确定性,可能会对正演模拟结果造成严重影响.本研究针对海底地震仪广角地震走时模拟中浅部结构的不确定性因素,使用RayInvr软件构建了理论模型和走时数据体,分析讨论了初始模型沉积层厚度与速度变化对正演模拟结果的影响.之后进行实测数据处理,在MCS2019-3测线数据的约束下建立初始模型,对OBS2011-1测线震相重新进行了走时拟合,将得到的模型结果与先前结果进行对比,得到以下结论:沉积层厚度及速度不确定性对走时模拟结果的浅部结构影响较大,深部结构影响较小;使用错误的速度与基底深度组合来拟合沉积层反射震相PsP,会使台站下方的地壳和地幔折射震相Pg和Pn,以及莫霍面反射震相PmP都出现走时小起伏,走时提前对应沉积层速度与厚度偏大,走时延后对应沉积层速度与厚度偏小;PsP震相缺失时,不同的速度-深度组合都可以很好地拟合震相,但会增加结果的不确定性.     

邻近算法在微地震震源位置和一维速度模型同时反演中的应用

微地震事件的空间分布可以用来监测水力压裂过程中裂缝的发育情况.因此,震源定位是微震监测中重要的环节.震源定位依赖准确的速度模型,而震源位置和速度模型的耦合易导致线性迭代的同时反演方法陷入局部极小值.邻近算法作为一种非线性全局优化算法,能够最大程度地避免陷入局部最优解.本文将邻近算法应用于单井监测的微震定位和一维速度模型...     

利用背景噪声层析成像方法反演新疆天山及邻区S波速度结构

利用新疆地震台网52个固定台站和天山地区新布设的11个流动台为期1年的观测数据,采用背景噪声层析成像方法获得了天山及邻区(41°～48°N,79°～91°E)10～50s的瑞利面波相速度分布图像,使用基于贝叶斯的马尔科夫链蒙特卡洛(MCMC)方法反演得到研究区地壳上地幔S波速度结构.研究结果表明,准噶尔盆地沉积盖层南深...     

地震波在孔隙介质中低频衰减现象的粘弹性特征分析及近似（英文）

由于介观尺度的孔隙流体流动,弹性波传播过孔隙岩层时在地震频段表现出较强的频散和衰减。Johnson理论给出了在任意孔隙形状的条件下,部分气水饱和孔隙介质的理论相速度和品质因子的解析解。本文在Johnson模型的基础上,通过对Q值曲线的低频和高频近似,推导了Q值曲线的近似公式,以及基于孔隙介质基本地球物理参数和孔隙斑块几何形态参数T和比表面积S/V的最大衰减Qmin近似公式。通过与理论值的对比,对Qmin近似公式存在的线性误差进行改正,进一步提高了精度。复杂的斑块形态对最大衰减Qmin和过渡频率ftr的都产生一定影响,且对ftr影响更大。因为数值模拟直接求解介观尺度的Biot孔隙介质方程需要极大的计算量,我们使用Zener模型建立了等效粘弹模型,有效地模拟了地震频带内的衰减和频散现象。     

汶川地震前后地震波速比和视速度变化特征的研究

2008年5月12日在中国四川省境内龙门山断裂带发生了8级巨大地震．在其主震前后随着震源区附近应力状态的改变,地壳介质的物性是否也存在相应的异常变化过程？该问题成为汶川地震研究的热点科学问题之一．本研究搜集和整理了中国地震局四川数字地震台网2001年1月1日-2010年5月31日产出的震相观测报告,和中国地震局地球物理研究所流动数字地震台网2008年5月12日——9月30日产出的震相观测报告,采用多台和达法和多地震联合测定法,根据入选的直达波Pg和Sg的震相数据,重点研究了龙门山断裂及其附近地区的P波和S波视速度的时间变化特征,并与波速比的时间变化特征进行了对比分析．为保证计算结果的可靠性和稳定性,对所收集到的数据进行了较为严格的筛选和限定．研究发现,龙门山断裂及其附近地区,在汶川地震前P波和S波视速度存在4年左右的明显降低过程,到震前约半年出现快速恢复．在震前P波和S波的视速度明显降低的过程中,波速比的低值异常过程仅为2年多,且幅度小持续时间较短．在波速比和视速度异常发展的进程中,进入发震中短期异常快速恢复是两者的共同特征．