海城地震震源区的地震尾波衰减特征

本文使用简化的能流模型(Frankel，A．et al.,1987)研究了海城地震震源区的地震尾波衰减及介质Q值的时空分布特征。结果表明,使用该方法求得的尾波Q值基本上是介质非弹性吸收的度量；海城7．3级地震前震的Q值较高，约为500，震后显著降低至380，但5．9级最大余震震源区的Q值仍与大震前相同。这可能说明大震发生时该处介质并未破裂。最大余震发生后，整个震源区的平均Q值降至344。但地震前后介质Q值发生明显变化的范围并不大，直径约为20km。     

粤西地区的地震尾波Q值和尾波衰减特征

根据Aki的尾波理论,利用短周期仪器记录的地震资料,用两种方法估测了粤西地区的尾波Q值,得到了如下结果;Q_0值为260—370,频率相关系数η值为0.0—0.2。阳江地区的Q_0值较低,其它地区大致相同.Q_0值的分布与地质构造和地震活动有较好的相关性. 研究了粤西地区的尾波衰减特征,得到尾波衰减率K为1.88—3.18,K值与Q值的分布彼此有较好的对应性。以上结果可为地震烈度衰减和地震危险性分析等方面的研究提供参考。     

大柴旦地震序列尾波衰减Q值变化特征研究

利用大柴旦数字地震台记录的2008年6.3级和2009年6.4级的数字地震资料,系统地分析了大柴旦地震序列尾波Q值随时间的变化,着重分析了尾波Q值在强震前的短期变化特征.结果表明,Q值在大震前后有较明显的变化,基本趋势是震前Q值增高,震后降低.     

含衰减地层微地震震源机制反演及其反演分辨率

微地震震源机制的反演对于非常规油气开发具有至关重要的作用.微地震信号主频高、能量小,容易受地层吸收衰减作用的影响使其波形发生畸变,本文提出了一种考虑地层吸收衰减作用的微震源机制反演方法,并利用费雷谢偏导矩阵的SVD分解（特征值分解）方法,分析研究了地层的吸收衰减因子的变化对于微地震震源机制反演分辨率的影响,根据理论计算给出了不同地震数据对各种震源机制反演的适用条件.理论计算证明,采用直达P波和S波数据联合反演震源的T值,单独利用直达P波反演震源的k值,可以有效降低地层吸收衰减作用对反演结果的影响.     

地震定量的国际现状

地球物理观测是全球性的,国际地震学和地球内部物理学联合会(IASPEI)下设的地震定量分委员会(主席:D.Denham)对世界各国和地震机构进行了地震定量情况的全面调查。我们想把该委员会的调查结果刊登出来,使大家对世界上测定震级的情况有一全面地了解,并促进我国的地震观测尽快与世界各国协调一致。     

陕西省地震动衰减及台站场地响应研究

利用陕西省测震台网的25个台站和2008—2010年所记录到的19个地震进行了陕西省Q值及台站场地响应值的计算,得出非弹性衰减Q值随频率f的关系为Q（f）=826.57f 0.3538。这说明陕西省处于高Q值区域,地壳相对较稳定,地震活动性相对较弱。计算得到的陕西省数字测震台网25个台站的场地响应值,其中,只有铜川台（...     

地震模拟城动台波形再现定量判别法

本文给出了地震模拟振动台数字迭代控制精度的定量判据以及实际应用的效果。     

用于深度域地震剖面衰减与频散补偿的反Q滤波方法

根据地震波在非弹性介质中的传播规律，本文提出了一种在深度域地震剖面进行反Q滤波的新方法. 在深度域的Q滤波算子符合地震波在衰减介质中的传播规律，不仅考虑了介质吸收对地震波振幅的影响，而且还保证了所造成的波形畸变满足因果规律，即地震体波具有某种频散性质. 利用该方法处理了2条总长度约60km的某浅海地区叠后深度偏移地震剖面. 合成记录剖面的正演和反演计算结果表明，在对地震剖面进行深度域反Q滤波之后 ，不仅达到了压缩中、深部目的层子波宽度提高分辨率和信噪比的目的，而且保持了原始地震剖面的整体形态和能量分布.     

云南地区尾波衰减及其与地震前兆的关系

本文用38个遍布云南地区的短周期地震台站的地震仪器记录的160多个区域地震的尾波观测资料,详细研究了云南地区尾波衰减的某些特征,给出了全省地壳尾波Q分布图。在考虑Q值与频率相点的情况下,发现频率相关值η在地震活动频繁的地区比较高,而在相对平静的地区比较低。从尾波Q值的空间分布来看,滇西北地区最低,在70—180之间,东部地区最高,约300—400左右,云南地区总的趋势是自西向东Q值逐渐增高。根据这些研究结果,讨论了尾波衰减的某些特征与地震活动前兆关系的理论基础,把尾波Q值作为地震前兆的可能性研究是具有实际意义的。     

利用互相关函数求地震波衰减

与传统利用单事件波形求介质衰减的方法不同,我们提出一种以一系列事件为源,以源与远处记录波形的互相关函数为基础计算介质品质因子的方法.依据野外现场实验的数据,我们分别以单次和多次激发为源,讨论了事件波形和互相关函数波形的振幅变化及相应的介质品质因子,结果发现,与事件波形一样,互相关函数波形也可以用于计算介质的衰减系数.计算结果表明：（1）实验场地浅层介质在原场条件下的品质因子Q值在40左右;（2）互相关函数的最大幅度与事件波形的振幅随距离具有相同的衰减规律.     

沙漠地区地震勘探随机噪声建模及其在噪声压制中的应用

随机噪声是影响地震勘探有效信号的主要因素,其存在大大降低了地震记录的信噪比.在噪声压制方法不断被改进的同时,对随机噪声特性进行研究,了解噪声的产生机制是对其进行压制的先决条件,目前对噪声的研究主要是特性研究以寻找规律性,对其进行定性定量的分析还比较少.本文根据塔里木沙漠地区实际采集环境,考虑到噪声的连续性给计算带来的不便,假设各类噪声源以点源的形式分布在检波器周围,依据相应理论确定各类噪声源的源函数,其激发的噪声经由波动方程传播,将随机噪声作为各类噪声源共同作用的综合波场,建立随机噪声的理论模型.通过分析不同种噪声对地震记录的影响,选取合适的滤波方法对其进行压制,实验结果表明,通过建立沙漠地区随机噪声的理论模型,为选择有效的滤波方法,提高地震记录信噪比起到理论指导作用.     

海南岛介质非弹性衰减及场地响应反演

利用海南数字地震台网18个地震台274条波形资料,运用遗传算法反演海南岛非弹性衰减系数及台站场地响应,计算得到海南岛非弹性衰减Q值随频率f的关系为Q（f）=343.7f^0.42,说明海南岛适中Q值区域,地壳相对稳定。18个地表台站的场地响应场与频率有关,其中：7个台站的场地响应接近1,平稳变化;2个台站在高频部分衰减;7个台站的场地响应小于1,变化平稳; 2个台站低频部分场地响应接近1,高频部分放大。     

地震尺度下碳酸盐岩储层的岩石物理建模方法(英文)

碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性,使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型,以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上,研究采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体,根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量,并根据不同孔隙连通性进行流体替换,计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量,采用Hashin-Shtrikman-Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析,明确不同孔隙对岩石弹性参数的影响特征,模拟分析结果与实际资料认识一致。     

基于门控循环单元主导的深度网络半监督地震衰减补偿方法

地震波在地下传播时存在能量衰减,呈现非平稳性,导致地震资料的分辨率降低,不利于储层刻画.为了解决这个问题,需要对地震记录进行衰减补偿来消除非平稳性.考虑到GRU(Gated Recurrent Unit)即门控循环单元在处理时间序列时具有长时和短时记忆的优势,本文构建了以GRU为主导的网络结构(简称GRU网络), 提出了基于GRU网络的地震衰减补偿方法.实际地震数据处理中,测井资料往往非常有限,导致标签数据较少.为解决小样本问题,本文借鉴自编码器的思想,除井旁道和测井合成地震记录的标签外,将其余地震记录作为无标签数据引入训练,半监督地学习非平稳数据到平稳数据的非线性映射,实现了地震记录的衰减补偿.最后利用含噪合成地震记录和实际地震资料测试了本文提出的方法,证明了其有效性.

     

A compound method for random noise attenuation

Random noise attenuation, preserving the events and weak features by improving signal‐to‐noise ratio and resolution of seismic data are the most important issues in geophysics. To achieve this objective, we proposed a novel seismic random noise attenuation method by building a compound algorithm. The proposed method combines sparsity prior regularization based on shearlet transform and anisotropic variational regularization. The anisotropic variational regularization which is based on the linear combination of weighted anisotropic total variation and anisotropic second‐order total variation attenuates noises while preserving the events of seismic data and it effectively avoids the fine‐scale artefacts due to shearlets from the restored seismic data. The proposed method is formulated as a convex optimization problem and the split Bregman iteration is applied to solve the optimization problem. To verify the effectiveness of the proposed method, we test it on several synthetic seismic datasets and real datasets. Compared with three methods (the linear combination of weighted anisotropic total variation and anisotropic second‐order total variation, shearlets and shearlet‐based weighted anisotropic total variation), the numerical experiments indicate that the proposed method attenuates random noises while alleviating artefact and preserving events and features of seismic data. The obtained result also confirms that the proposed method improves the signal‐to‐noise ratio.     

中国东北及周边地区地壳横波衰减的成像研究

利用146个台站记录的中国地震年报及黑龙江、吉林、辽宁和内蒙古四省地震台网观测报告的6371个近震46652条横波振幅与周期资料,采用二维层析成像方法反演东北地区地壳介质品质因子Q0空间分布图像.在绝大部分研究区都可达到2°×2°的分辨率.结果表明,东北及周边地区地壳介质品质因子Q0平均值为323,其Q0值分布及其所揭示的衰减变化特征与研究区的地表构造具有明显的相关性,如华北平原、松辽盆地和海拉尔平原等沉积层较厚的区域通常呈现为低Q0值;而燕山、太行、兴蒙褶皱带等沉积层薄或基岩出露的区域则呈现为相对高的Q0值.整个东北地区呈现高低衰减带相间分布.     

模拟地震波传播的三维逐元并行谱元法

高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre (GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性.

     

基于多辛结构谱元法的保结构地震波场模拟

近年来构造高精度、高效且具有长时程跟踪能力的保结构算法已逐渐成为地震波模拟算法发展的重要方向之一. 本文基于谱元法(SEM)进行空间域离散结合新推导的三阶辛算法(NTSTO)进行时间域离散,构造了一种具有时-空保结构特性的新算法. 本文给出的多组数值试验对比结果表明,本算法无论在内存消耗、稳定性及计算耗时,还是长时程跟踪能力方面都有上佳的表现; 另外,本文给出的起伏地表多层介质模型的数值算例验证了该算法处理复杂几何形状和复杂介质时的有效性. 该多辛结构谱元法的发展将为长时程地震波传播的计算及模拟提供更为广泛而有效的选择.     

井间地震物理模拟研究

井间地震资料具有极高的分辨率,但是,其波场十分复杂,特别是有效一次反射波场能量较弱,且往往被管波等强相干干扰淹没,波场的识别与分离难度较大.地震物理模拟技术是认识复杂地震波场的有效手段.我们通过单一地质体模型、复杂地质体模型和真实井间模拟模型等多种逼真地质模型的物理模拟,揭示了井间地震观测下直达波、透射波、反射波、折射波,多次波和导波等各类波型的特征.物理模拟试验表明,在炮检域平面上分析直达波的信噪比与能量分布,可以使处理人员对旅行时的误差有更直观的认识.所以识别初至与拾取旅行时间时,不仅要在共震源点道集(CSG)上分析拾取,而且还要在共炮检距道集(COG)上分析拾取,这样就提高了直达波识别的可靠性与拾取时间的精度.