基于L-BFGS算法和同时激发震源的频率多尺度全波形反演

全波形反演可以利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,庞大的计算量和存储空间需求,限制了全波形反演的发展。在频率多尺度全波形反演中将L BFGS数值优化算法与同时激发震源技术相结合的方法来改善这一现状。首先,对Marmousi模型进行了速度反演：在计算过程中明显发现对计算机内存的占用减少,最终反演结果与实际Marmousi模型的拟合误差为0.095 9,较小;采用10个频带单炮震源正演384炮所需时间约为32 640 s,而采用同时激发震源（384炮）正演一次所需时间仅约为700 s。然后,基于高速楔形体模型进行了抗噪能力研究：原始含噪地震记录信噪比为11.147 3 dB;对反演得到的速度模型进行正演,其地震记录信噪比为22.251 8 dB。最后,基于逆冲断层模型进行了反演速度扰动能力研究,反演得到的最终模型很清晰,与具有速度扰动特性的实际模型非常接近,拟合误差仅为0.036 0。数值模拟试验结果表明：此方法反演精度高,内存开销较小,能够显著提高计算效率,并且具有良好的抗噪能力,能够反演出具有速度扰动特性的介质。     

利用地震属性解释煤层冲刷带

把地震属性技术和地震相分析技术应用于煤层冲刷带的解释工作中。首先从煤层反射波中提取了多种地震属性;然后分析对比地震属性与煤层冲刷带的响应特征;最后利用响应特征明显的地震属性对煤层冲刷带进行了综合解释,确定了煤层冲刷带的范围。与传统钻探解释方法相比,地震属性方法能够更准确地圈定煤层冲刷带的边界。     

陆相汕藏中砂泥岩互层时移地震遮掩效应及其对时移地震资料解释的影响

从差异波形特征的角度出发,分析时移地震差异波形特征与砂泥岩互层油藏中油层速度变化的关系.通过建立由5层地层组成的楔形地质模型,运用褶积公式及时移地震求差公式,得到泥岩速度为2 400 m/s和砂岩速度为2 700 m/s条件下,互层中砂岩速度变化±150 m/s的时移地震响应.研究结果表明砂泥岩互层的差异波形是一个多峰复合波,泥岩夹层的厚度是决定差异波形特征的一个主要因素,并且存在一个临界值,即四分之一主波长;如果泥岩夹层的厚度小于这一临界值,就会发生遮掩效应,即顶部油层的时移地震响应可以掩盖中、下部油层的时移地震响应.     

陆相油藏中砂泥岩互层时移地震遮掩效应及其对时移地震资料解释的影响

从差异波形特征的角度出发,分析时移地震差异波形特征与砂泥岩互层油藏中油层速度变化的关系。通过建立由5层地层组成的楔形地质模型,运用褶积公式及时移地震求差公式,得到泥岩速度为2 400 m/s和砂岩速度为2 700m/s条件下,互层中砂岩速度变化±150 m/s的时移地震响应。研究结果表明:砂泥岩互层的差异波形是一个多峰复合波,泥岩夹层的厚度是决定差异波形特征的一个主要因素,并且存在一个临界值,即四分之一主波长;如果泥岩夹层的厚度小于这一临界值,就会发生遮掩效应,即顶部油层的时移地震响应可以掩盖中、下部油层的时移地震响应。     

山西数字地震台网仿真前后震级偏差对比分析

选取山西数字地震台网2010年不同震级,不同地点的地震事件,按照中国地震局的编目要求对波形进行WA仿真处理,对比仿真前后的震级大小,找出其中存在的偏差,为地震资料的衔接提供参考.     

景泰地震台短周期地震仪记录的几种异常波形特征与景泰6.2级地震的关系

本文研究了景泰地震台记录到的几种异常地震波形的形态特征及随时间分布特征,采用测震学的方法和数理统计的方法,对各个阶段的异常波形出现的数量、类型和振幅进行了定量分析,发现在1990年10月20日景泰6.2级地震前,异常波形数量明显地呈平静一密集交替出现的现象;临震前各种异常波形的数量增多,振幅增大,异常持续时间增长,震后异常波形很快衰减并消失.与其他前兆对比,可以认为异常波形的出现可能与景泰地震的孕育过程有关,根据该类波形的变化特征有可能对地震作出预报.     

Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.     

利用背景噪声研究滇西北波速异常与M≥5.0地震关系

提出基于背景噪声波速测量的综合预测指标法,将测震资料更好地应用于震情跟踪和地震短临预报。利用滇西北5个台站2012-01~2020-11宽频带连续波形资料,基于背景噪声互相关及傅立叶变换等方法,提取10个台站对当天经验格林函数与参考经验格林函数的直达瑞利波走时偏移时间序列,设定±1.5倍标准差作为异常阈值,并以其间发生的6次M≥5.0地震为样本,采用R值评分法对每个台站对的映震能力进行效能检验,最后基于自适应加权综合预测方法提取适合于滇西北地区的地震短临异常识别指标（综合指标）。结果表明,利用该综合指标对滇西北2012年以来发生的6次M≥5.0地震进行90 d短临预报,异常指标共出现8次,其中准确预报地震5次,漏报1次,虚报4次,预报效能评分R为0.692,R₀为0.475。该综合指标的地震对应率为62.50%,概括率为83.33%。     

基于SEED波形文件的运行率统计

介绍了基于SEED波形文件运行率计算的技术和实现方法,介绍了SEED文件中标识时间的特殊掩码、时间片的具体格式和程序的实现,并给出了与其相关的部分源代码.