基于奇异谱分析的静态相对重力观测重力固体潮提取

针对相对重力测量数据中存在的重力固体潮信号,本文提出将奇异谱分析方法（Singular Spectrum Analysis,SSA）应用到相对重力测量数据处理中,在不需要测站坐标等先验信息的条件下从相对重力数据中提取重力固体潮,提供了一种获取重力固体潮的新思路.采用模拟的相对重力数据进行实验,利用SSA方法和小波变换方法分别从模拟信号中提取重力固体潮并进行结果对比,SSA获取的重力固体潮与理论值残差RMS为0.3 μGal,小波方法获取的残差RMS为1.6 μGal.利用CG-5相对重力仪实测数据进行实验,提出一种利用SSA外推时间序列来削弱边界效应的新思路,实验结果显示采用这种方法后重力固体潮值与理论值残差序列的RMS和STD均有所减小.通过实验发现削弱边界效应后SSA提取的重力固体潮与采用Tamura潮波表计算的重力固体潮理论值残差RMS值为2.2 μGal.利用SSA提取的零点漂移值与最小二乘拟合得到的结果基本一致,十天内的差值小于0.4 μGal/d.     

历史大地震断层滑动模型的建立及其对同震数值计算的影响——以1920年宁夏海原MS8.5大地震为例

大地震的发生会引起区域位移场和应力场发生变化,进而改变区域内及临近断层的应力状态和地震活动性.目前,研究学者可据已有的断层滑动模型来计算分析大地震同震应力变化,同时采用库仑应力触发理论来进一步分析震后余震分布和断层危险性.然而,历史上曾经发生过不少大地震,例如,1920年的海原M_S8.5大地震,是全球范围内少见的特大地震之一.局限于无确切的地震台站地震波等资料,前人在研究历史地震的影响时往往给出一些简单的断层滑动模型,将断层面上错动量视为均匀分布.为更准确地了解历史地震对后续地震的影响,基于前人研究和一般地震滑动形态分布规律及地震反射剖面等资料,以海原M_S8.5大地震为例,探讨了如何建立海原大地震断层滑动模型,并分别搭建了简单断层滑动模型和复杂断层滑动模型的全球同震横向不均匀并行椭球型地球模型.通过对海原M_S8.5地震的同震位移场和应力场的计算,发现采用复杂断层滑动模型比简单断层滑动模型地表位错分布更切合实际.同时,进一步计算和分析了此次大地震对青藏高原东北缘近100年历史地震和周围断层的应力触发作用,得出断层滑动模型对同震计算结果的影响集中在发震断层附近而对远场影响较小.     

2001年昆仑山口西MS8.1地震对周围断层的应力影响数值分析

大地震的发生往往会引起周围区域形变场和应力场变化,且对临近断层上的应力状态也有影响.2001年11月4日,昆仑山口西发生了半个世纪以来中国最大的M_S8.1级地震.本文基于已有的滑动模型,建立了三维含地形高程的横向不均匀性椭球型地球有限元模型,采用等效体力方法,分析了此次M_S8.1地震产生的全球同震位移和应力场变化.与解析方法相比,该模型考虑了地形、Moho面起伏和地球介质横向不均匀性;与一般的有限元数值模拟相比,该模型考虑了地球曲率和椭率,合理地规避了有限块体模型假定边界位移为零所引入的误差.计算得出同震位移与GPS观测数据可以很好地吻合.据库仑破裂应力准则和震源参数,计算得出昆仑山口西M_S8.1地震的发生造成了汶川、芦山、改则和当雄地震的发震断层上库仑应力增加,对这些地震的发生起促进作用;而造成玉树和德令哈地震发震断层上的库仑应力变化为负值,在一定程度上抑制了这些断层的地震活动性.此外,计算结果显示地球地形高程、介质非均匀性和椭率对昆仑山口西M_S8.1地震同震变化计算有一定的影响,其中地形和椭率造成的同震位移场相对误差约10%.     

浙江省梅汛期洪涝灾情分析和预测

通过所构设的雨涝强度、雨涝影响面积等统计量来估计灾情，有较高的正确率，能估计出极端情况，为灾害评估工作提供了新的思路。     

气汞的第一个同震效应记录

利用杭州超距科技有限公司最近研发的ATG-6138型测汞仪在云南省弥勒井气汞（Hg）观测中记录到尼泊尔8.1级地震的同震效应。这是我国气汞观测的第一个同震效应的纪录。这一事实表明了进一步提升地震地下流体化学量观测仪器的观测精度与采样率,有可能记录到固体潮与同震效应等更多的地壳动力作用信息,从而进一步提升其地震前兆监测能力。     

基于二次场的可控源电磁法三维有限元-无限元数值模拟

本文提出了一种新的混合有限元-无限元三维可控源电磁法(CSEM)问题快速高精度正演模拟算法.首先从电场双旋度方程出发,推导了水平电偶极子源的二次场边值问题,采用无限元代替截断边界条件和有限元离散内部计算区域的新策略,达到减小计算区域的目的,基于并行直接求解技术,实现多源CSEM问题的快速精确求解.其次,通过层状解析模型测试,一方面验证了新算法的正确性,另一方面通过与其他三种已知CSEM问题求解策略进行对比,表明了本文提出的基于二次场有限元-无限元算法具有离散区域小、求解速度快和计算精度高等优点.最后,通过3D模型计算,清晰直观地模拟了场源阴影效应,为野外数据的处理与解释提供指导.     

东北冷涡对海河流域初夏降水异常的影响

东北冷涡是东亚中高纬度地区的重要天气系统,其频繁活动可能导致显著的气候效应,不仅对东北局地,而且对海河流域地区初夏降水产生影响.本文利用海河流域34个地面气象台站逐日降水资料以及NCEP/NCAR再分析资料,对东北冷涡与海河流域初夏降水的关系进行研究,结果表明:初夏东北冷涡活动与海河流域降水存在显著的相关关系,对其空间分布亦有重要影响.当东北冷涡活动偏多(少)时,海河流域北部地区(尤其是海河下游)初夏降水可能偏多(少),对应时间系数为正(负)时海河流域初夏降水EOF第二模态北多(少)南少(多)的空间分布.在东北冷涡活动偏多年,东北冷涡引导北方"干冷"空气南下推进至海河流域北部,伴随该地区源自西风带水汽输送和水汽辐合的增强,形成了"上干下湿"的不稳定层结,在上升运动的触发下,导致海河流域北部初夏降水偏多;而东北冷涡活动偏少年情况正好相反.东北冷涡活动为海河流域尤其是海河下游地区初夏降水预测提供了具有参考意义的结果,其活动指数与海河下游初夏降水距平的符号一致率在近30年高达83%.     

场地效应研究的传递台站线性反演法——以芦山地震为例

采用线性反演法研究了地震动的场地效应,并针对线性反演法的局限性,提出了"传递台站"的方法来扩展其应用.通过传递台站,地震事件不再局限于必须被参考台记录到,台站也不再局限于须与参考台站同时获取到地震记录,改进后可大大增加数据样本,使得结果更趋稳定.同时,对于原来难以参与分析的台站,也可得到其场地反应.研究中,还采用了基于遗传算法的非参考台方法,进一步分析了场地反应和验证所提方法的合理性.将所提方法应用于2013年芦山7.0级地震周边区域,分析了这一地区的强震台站的场地反应和地震波衰减因子Q_s值.研究结果显示,当采用一个和二个传递台站进行分析时,所能利用的地震事件和所能求解场地反应的台站数目逐步增加,随着数据样本的增加,分析结果趋于稳定.研究区域的Q_s值拟合结果为Q_s(f)=334.16f~(0.73).遗传算法结果表明,线性反演法中所选取的参考台站有其自身的场地反应.当利用遗传算法获得的参考场的场地反应去校正线性反演法的结果时,对于多数场地,两者的结果符合得很好,说明了遗传算法的结果可以合理地认为是场地较为真实的反应.     

福建南一水库气枪震源实验激发效果最优分析与研究

气枪主动源探测实验需要获得信噪比较优的波形数据,因此实验前很有必要进行不同工况的激发实验。通过综合分析台站记录特征,寻找最优激发效果的工况条件,为后续定点激发奠定重要基础。本文在福建南一水库大容量气枪震源实验时通过单次激发分析、叠加分析、信噪比系数等方法综合研判工况激发效果,确定出最优工况,再利用台网噪声评估对得到的最优工况做进一步验证。结果表明:南一水库最优工况为:水深25m、枪阵沉放深度12 m、枪阵尺寸7 m×7 m。该最优工况激发时段未处在低噪声时刻,说明该工况的激发效果确为优越,结果较为可靠。