基于分区多步快速行进法下3D起伏层状介质中多震相射线追踪

快速行进法(FMM)是一种求解程函方程数值解计算网格节点走时,然后向后处理进行射线追踪的方法.为了求取任意起伏界面下高精度多震相的地震走时与相应的射线路径,本文采用任意起伏地表条件下的的三维不等距上行差分公式结合分区多步计算技术实现了三维复杂层状起伏介质中多震相(透射、反射、转换波)地震走时的计算,利用上行有限差分公式逐次进行射线路径的追踪,并且通过与较为成熟的不规则最短路径法(ISPM)对比,验证了本算法的计算精度和有效性.数值模拟实例和对比结果表明该算法具有较高的计算精度,数值计算稳健,能灵活处理含任意三维起伏界面模型中多震相地震走时及相应射线路径的追踪问题.     

基于HAFMM的无射线追踪跨孔雷达走时层析成像

本文使用最小二乘线性迭代反演方法对跨孔雷达直达波初至时数据进行反演,每次迭代过程中,用有限差分法求解走时程函方程,并用高精度快速推进方法(HAFMM)进行波前扩展,通过追踪波前避免了进行射线追踪.为了验证该方案,我们对三组合成数据进行了测试,分析了单位矩阵算子、一阶差分算子和拉普拉斯算子等三种不同模型参数加权算子对模型的约束和平滑效果;讨论了FMM和HAFMM对反演精度的影响;测试了LSQR,GMRES和BICGSTAB等三种矩阵反演算法的反演效果.此外,我们还对一组野外实测数据进行了反演,对比了基于本方案以及基于平直射线追踪和弯曲射线追踪的走时层析成像反演效果.对比分析结果表明,使用拉普拉斯算子和HAFMM进行反演能较好地进行目标体重建,而三种矩阵反演方法对反演效果的影响差别不大;并且通过对波前等时线图的分析可以定性地判断异常体的性质和位置;而在对实测数据目标体的重建上,本方案能达到甚至优于弯曲射线算法的重建效果.     

基于MSFM的复杂近地表模型走时计算

地震走时层析成像方法是解决复杂近地表模型速度建模问题的重要技术.该方法是一种迭代反演方法,在反演过程中需要反复计算地震射线走时.故而,高效高精度且能适应复杂模型的走时计算方法是地震走时层析成像实用化的关键技术之一.本文引入医学成像领域研究的MSFM(Multi-stencils Fast Marching Methods)用于地震层析反演中的走时计算.该方法在标准FMM(Fast Marching Methods)基础上利用坐标旋转生成新的FMM计算模板,使计算网格点对角方向邻点参与计算,改善了标准FMM存在对角方向误差大的缺陷.本文分析对比了MSFM和标准FMM的计算精度和计算效率;针对地震层析成像技术解决的起伏地表模型建模问题,研究了起伏地表模型地震走时计算的MSFM实现方法;采用炮点邻近区域局部细分网格技术只需增加很少的计算量即可大幅提高计算精度.理论分析和模型试算表明MSFM算法明显改善了FMM的计算精度,同时保持了FMM算法的高效性.文章通过对崎岖地表模型的正演和层析反演试算,验证了基于MSFM的地震走时计算方法对复杂模型有很强的适应能力.研究表明该方法作为地震走时层析反演中高效高精度的正演算法,有很好的应用价值.     

慢度平方三角网格立体层析反演方法

三角剖分模型是一种稀疏的模型描述方式.在慢度平方三角剖分模型中,界面可以得到自然的描述,射线追踪在一个三角块内存在解析解.基于慢度平方三角剖分模型下的射线扰动理论,实现了立体层析数据空间对慢度平方三角剖分模型要求的所有模型分量的偏导数求取.相比常规的B样条或者大网格模型描述方式,慢度平方三角剖分模型下的立体层析FRECHET导数矩阵规模被大幅压缩,求解精度得到良好保证的同时计算成本得到大幅降低.理论数据算例证实了上述观点.     

基于二维三次卷积插值算法的辛几何射线追踪

射线追踪是地震波走时层析成像的基础,射线空间位置的准确性及射线走时的精度决定了层析成像的可靠性.本文根据哈密尔顿系统可以有效提高程函方程解稳定性的特性,采用辛几何算法（SAM-Symplectic Algorithm Method）及二维三次卷积插值技术进行地震波射线追踪.由于采用了SAM算法,保证了地震波波前精度,提高了射线空间位置的准确性.数值模拟结果表明SAM既能保证哈密尔顿系统的稳定性又具有运算速度快的特点,提高了射线追踪的计算精度.     

基于波前面三角形网格剖分的波前重建法三维射线追踪

在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用.     

最短路径射线追踪方法及其改进

综述了用网络最短路径算法求解地震射线追踪问题的原理、方法技术以及存在问题和改进措施。特别介绍了作者在最短路径算法基础上，提出的动态网络最短路径地震射线追踪方法。该方法先采集从炮点到整个模型所有节点上的初至旅行时，其中，在一个单元内，对相邻每对已计算出最小旅行时的节点进行线性插值，并利用Fermat原理计算未知节点的最小旅行时；然后，利用同样的方法，从接收点开始，反向追踪炮点到接收点的射线路径、该方法能适于各种复杂的非均匀介质，极大地提高了射线追踪的精度。     

TI介质射线追踪及其高斯束成像应用

射线追踪是一种高频近似前提下快速有效的波场近似计算方法,传统的基于弹性参数的各向异性介质运动学和动力学射线追踪方程,求解过程中需要处理Christoffel方程的特征值问题,因而计算效率比较低.为了解决这一问题,本文通过引入相速度和群速度,对运动学和动力学追踪方程进行修改和简化,有效地提高了各向异性介质射线追踪算法的计算效率.另外,我们将该算法应用到各向异性偏移中,实现了共炮域TI介质高斯束叠前深度偏移方法.VTI介质Hess模型和TTI介质洼陷模型的试算结果说明了该方法的正确性和有效性.     

基于图形结构的三维射线追踪方法

在地震层析成像研究中,为了克服最小走时射线路径追踪方法存在的问题,对该方法计算过程中的关键步骤进行了改进．在节点走时的计算中引入Ｂｒｅｓｅｎｈａｍ画线算法;在最小走时节点查寻中,结合使用快速排序算法与插入排序算法,替代以往方法中多采用的堆排序算法;所采用的节点设置方式,可以引入速度界面,还可以实现反射波射线追踪．模型计算证明,改进的最小走时射线路径方法具有精度高,速度快的特点,所提出的三维空间反射波射线追踪算法简便易行。     

论地震波初至走时的计算方法

本文主要论述了近来出现的计算地震波初至走时的一些新方法。这些方法可以分为两大类,一类是利用基于解决最短路径问题的方法来求解地震学中任意两点间的射线路径及走时;另一类是直接求解程函方程。这些方法不仅适用于任意几何位置的源点-检波点排列,而且适用于介质速度的任意变化,其计算精度和效率也比通常使用的射线追踪方法如试射法和弯曲法要高。     

兰州地震研究所发展的三段历史回顾

对兰州地震研究所的三段历史作回顾,分别为国家地震局、中国地震局兰州地震研究所的体制沿革变迁,中国科学院地球物理研究所;中国科学院兰州地球物理研究所观象台早期创建;兰州地震研究所硕士点早期阶段。     

山西大同数字遥测地震台网记录的震相特征

从地球及地壳的分层结构出发，使用较成熟的理论，从地震波传播的基本原理入手，结合多年的震相观察分析工作实践，选取了山西大同数字遥测地震台网近年来的一些典型震相进行了分析，得出本台网记录到震中距小于100km的地方震震相有Pg，P11，Sg，S11震中距在100km～500km的近震震相有Pn，Pg，P11，Sn，Sg,S11对比分析了在本台网内及周边地区发生的爆破与天然地震在震相上的差异性，并对其特征作了较详细的描述。     

淋巴结结核的影像学表现特点分析

目的:分析淋巴结结核的CT和MRI表现特点,并与病理学对照,探讨两种检查方法的应用价值。方法:回顾性分析2010年1月至2014年12月间34例经手术、穿刺病理或临床证实为淋巴结结核患者的CT或MRI影像学资料,分别从生长方式（融合/非融合）、边缘（光滑/欠光滑）、境界（清晰/欠清晰）、密度及信号、增强后强化方式（均匀/不均匀/环形规则/环形不规则）这5个方面对所受累淋巴结进行分析。结果:34例患者中共累及223枚淋巴结,其中颈部淋巴结67.7%,纵隔内淋巴结27.7%,肺门淋巴1.8%,肺内淋巴结0.5%,腹腔淋巴结0.5%,腹膜后淋巴结1.8%。所有受累淋巴结中,8.5%呈融合性生长,83.9%形态规则,79.4%边缘光滑。CT上发现的171枚病变淋巴结以等密度较多（53.8%）、低密度次之（43.3%）,MRI上发现的66枚病变淋巴结中98.5%为等T1长T2信号。所有CT及MRI观察到的病变淋巴结中,除57%的病灶密度或信号均匀外,仅伴钙化的为12.6%,仅伴部分坏死的26.9%,另外有3.6%同时伴钙化及坏死;增强扫描的204个病变淋巴结中,1%钙化显著、未见明显强化,38.7%均匀强化,伴坏死者内见不强化的低密度区,11.3%呈环形规则强化（均位于颈部）,8.3%环形不规则强化（纵隔内82.4%）,5.9%可见病变淋巴结内有分隔样强化（均位于颈部）。结论:CT与MRI对于淋巴结结核具有重要诊断价值,两者均能清晰显示病变部位、形态以及边缘情况,MRI对于病变边缘及周围情况显示较CT更明确,尤其是增强后,对于包膜及邻近组织受累情况的显示也更加敏感。      

Correction of twilight measurements of influence of night sky

Studia Geophysica et Geodaetica -     

山西夏县中心地震台温泉水化学特征

采集夏县中心地震台温泉水及周边水点样品进行水化学组分特征分析。利用矩形图、Na-K-Mg三角图等方法,初步分析温泉水的水质类型、水-岩平衡状态、热储温度以及循环深度等,并结合氢氧同位素组成特征,初步分析温泉水补给来源。研究表明,夏县中心地震台温泉水化类型属于Na-Cl·SO₄型,水-岩反应属于部分成熟水,热储温度为148.8℃,循环深度为3.91 km,补给源主要为大气降水,温泉为断裂型温泉。以上结果可为该台流体异常分析提供基础研究资料,为夏县地区水化研究奠定基础。     

山东数字地震台网子台背景噪声位移值的分布特征

利用山东数字地震台网的实时记录波形资料,采用在频率域积分方法,计算出各台站背景噪声的位移值。然后进行统计分析,得出了山东地区背景噪声位移值的空间分布和时间变化特征：台站背景噪声幅值与不同地区地质环境、岩性密切相关,位于第四系黄土覆盖层较厚的鲁西聊城、菏泽地区,鲁西北德州地区和鲁北东营、滨州地区各子台的背景噪声位移值较大,鲁中、鲁南山区和胶东半岛地区的基岩台站背景噪声位移值较小;同时,白天噪声值大于夜间噪声值,但深井观测台站的背景噪声值相对稳定。此项研究为绘制台网监测能力图提供基础数据,并为计算各台站的场地响应及震级修正值等测震学研究提供基础资料,促进台网的改建和发展。