响应谱地震动预测方程中场地效应参数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年多,且幅度小持续时间较短．在波速比和视速度异常发展的进程中,进入发震中短期异常快速恢复是两者的共同特征．      

利用初至P震相测定石嘴山ML4.4地震序列的震源深度

基于初至P震相较易辨识、读取误差小、有效震相丰富等优点,提出联合Pg与Pn震相到时数据测定震源深度的方法,并将其应用于石嘴山M_L4.4地震序列(11次)的震源深度计算,结果显示该地震序列的震源较浅,这也是此次有感地震震感强烈的原因之一．然后利用初至P震相到时方法分别计算了双层和4层地壳速度模型下的震源深度,并与双差定位所得深度进行对比,结果显示: 利用初至P震相到时方法较双差定位方法得到的震源深度整体上一致性较好,但前者得到的震源深度较之后者略深;对于双层与4层地壳速度模型,二者的震源深度计算结果具有较好的一致性．因此,基于可靠区域地壳速度模型下的初至P震相到时方法可以应用于震源深度的计算,并能够取得较好的结果．另外,石嘴山M_L4.4地震序列的主震深度为7—8 km,最大余震深度为6 km,根据银川盆地的速度成像结果可知,主震和最大余震发生于银川盆地基底下部;其余地震的震级偏小,多集中于3—5 km深度,主要发生于银川盆地基底顶部或基底覆盖层内．      

云南漾濞M6.4地震震区三维速度结构

利用2015年1月至2021年5月28日期间我国云南省漾濞县及周边地区固定台站和漾濞地震后布设的流动台站所记录到的近震资料,使用双差层析成像方法获得了该地震震区的高分辨率地壳三维速度结构和震源位置。重定位结果显示,漾濞M6.4地震序列主要沿NW?SE向展布,与维西—乔后—巍山断裂走向一致,地震主要集中在4—10 km的深度范围,呈约80°高倾角分布。结合定位结果与三维速度结构显示:漾濞M6.4地震序列的空间分布与速度结构变化具有相关性,主震位于P波、S波高低速异常交界处,这种介质物性变化的交界地带可能有利于中强地震的孕育和发生,余震主要分布在低P波速度、高S波速度和低波速比的脆性区域;沿漾濞地震序列的分布走向,主震两侧呈现完全不同的速度结构,其西北部具有明显的高P波速度、低S波速度特征,该地区高密度、强韧性的地层可能是阻挡漾濞地震的NW向破裂而呈单向破裂特征的原因。      

中国南北带地壳上地幔三维面波速度结构和各向异性

本文用长周期763地震仪面波群速度资料反演了中国南北带及邻区的三维速度结构.其中采集238条瑞利波和358条勒夫波混合频散曲线,使用均等显示滤波方法,并以4°×4°为一格将我国境内分为147格.用随机逆反演方法得到了研究区16格的纯路径频散.面波速度结构及演结果表明:1.莫霍界面深度一般在40-50km之间,最深达65km.总趋势是从东到西加深,且在南北带西侧南北两端向中部明显加深,东侧变化小.2.地幔顶部普遍出现很厚的低速层,上界面一般埋深60-80km.上地幔顶盖厚度一般为20-60km,速度为4.30-4.50km/s.3.研究区普遍存在各向异性,而且勒夫波和瑞利波速度的差值(V_SV-V_SH)的绝对值随深度有增大的特点,在南北带南部和西北部V_SV-V_SH各向异性现象更为明显.     

用有限区域风速场准确求解流函数和速度势场的方法

流函数和速度势是气象业务和研究中常用于表述风速的一组变量。用有限区域风速场, 使用有限差分方法求解得到的流函数和速度势场重建初始风速场, 由于受区域边界的限制往往有明显的偏差。虽然有许多求解方法的研究, 但是, 至今仍尚未见到一种真正准确的求解计算方案。首先, 介绍用Arakawa A网格和D网格分布的有限区域风速场求解流函数和速度势场的一般有限差分计算方法, 探讨用它们的解重建风速场产生误差的原因。然后, 针对这些原因, 对给定的有限区域, 通过线性外推初始风速场, 扩展求解计算区域, 使用协调、一致的有限差分格式方案, 准确计算求解区域的边界有旋风速、散度风速和速度势的定解边界条件, 以及恰当选择流函数、速度势、涡度和散度等变量的分布网格, 设计了用上述两种网格分布的风速场准确求解流函数、速度势场的方案, 并对其正确性加以证明, 它们可以推广应用于其他Arakawa网格。用实际资料试验同样显示, 方案避免了重建风速场误差的出现, 与初始风速场相比, 全场风速最大偏差精度达到10-12m/s或以上, 在计算机精度造成的计算误差影响范围内。本文的研究很好解决了长期以来用有限区域风速场、 使用有限差分方法无法准确求解流函数和速度势场的问题。