基于小波分析的断层检测

断层在地震数据中显示出奇异性,经过小波变换可以得到断层位置处地震数据的奇异性属性.对包含断层的地震数据进行小波分析处理,能够得到断层的垂直和水平位置.在实际地震资料中将地震信号表示成不同尺度和不同位置的基本单元,然后对变换系数进行极值提取,检测出不同尺度下的地震信号突变特征,从而进行断层检测.对实际地震资料进行地震信号奇异性检测时,首先将地震剖面划分成层,然后在每一层内将尺度参数进行离散化,计算地震记录的小波变换系数,对于某一个尺度求取每一道小波变换系数的最大值,将每一道地震记录小波变换系数的最大值根据原地震道的位置进行排列,得到奇异性曲线.对于某一尺度,断层所在的水平位置对应着奇异性曲线的最值位置,最后绘出整个剖面的极值点检测结果.     

中国东部活动断裂的现代构造运动

本文分析了中国东部第四纪以来活动断裂的展布和基本特征。根据断层位移测量把活动断裂的现代运动分为四级:第一级为≥10毫米/年,第二级为<10—1毫米/年,第三级为<1—0.1毫米/年,第四级<0.1毫米/年。并根据活动断裂的展布和现代活动程度,把该区分为五个相对活动区和三个相对稳定区。指出了断裂活动的强弱与地震活动有一定的相关性,近年来强震则多发生在第一、二级断裂带上,相反,活动速率较低的断裂,近期地震活动水平不高,也未发生过强震。本文最后讨论了应力作用方式     

多尺度奇异性检测及其在地震相分析中的应用

本文从地震道的奇异属性出发,利用连续小波变换求取地震道的小波变换模极大值曲线,并沿此曲线提取地震道的小波变换系数,称之为小波变换模极大值连线振幅.此属性不仅可以代表信号本身,而且可以最大程度的区别于相邻道,并且具有地震道多尺度的特征,即兼具时频域的特征.由此结合自组织神经网络,我们提出了一种新的地震相分析方法.通过模型合成地震记录实验分析,证明此方法是可行的,且对地质层位的解释误差具有一定的容许度.最后,将此方法应用于了实际资料,取得了良好的效果.     

地震资料预测薄层厚度方法研究与应用

薄层厚度预测一直是公认的难题之一,其难度就在于如何准确地识别和提取薄层的地震属性.常规方法是利用时间域或频率域地震属性与薄层厚度的线性关系计算.但是理论与实际资料表明,不同的薄层和地层组合对地震波的动力学的信息影响很大,各种参数与薄层厚度成非线性关系,使用单一的信息不可能准确预测薄层厚度.本文利用三种线性预测原理（模型）,经数学变换为属性参数,采用非线性BP网络预测薄层厚度,取得了令人满意的效果.     

利用叠前地震数据预测裂缝储层的应用研究

叠前地震资料储层预测技术是在Zoeppritz方程基础上发展起来的,通过处理地震数据随着不同入射角地震反射属性,得到地震属性随着入射角变化而改变,研究分析得到反映岩性变化的纵波速度、横波速度、泊松比和截距与梯度剖面,预测裂缝储层的发育及分布.同时根据不同方位角地震资料属性,计算得到不同方位角目的层的属性差异,使用各向异性椭圆公式作拟合,求出背景趋势A和各向异性因子B,利用最大振幅包络方位和对应θ,求出裂缝发育优势方向,及B/A各向异性因子,实现对裂缝储层预测.     

Contourlet变换在局部角度属性提取中的应用

为了直观展示地质的构造走向、暗示细小裂缝信息,便于地震资料解释,应用Contourlet变换和无抽样Contourlet变换提取地震剖面的局部角度属性,提出了基于Contourlet变换的局部角度属性提取方法.首先用Contourlet变换或无抽样Contourlet变换将地震剖面分解成为不同尺度下的方向子带,然后将不...     

用平面波延拓方程进行地震数据的叠前速度反演

本文讨论地震勘探数据的叠前速度反演方法及其在海洋地震勘探数据上的反演试验.反演主要的计算步骤是:1.采用Fourier-Hankel变换把球面波分解为平面谐波;2.用平面谐波的延拓方程将上行波与下行波同时向下延拓,并计算每一层底部的反射系数和下一层的波阻抗;3.用最小二乘法从波阻抗中确定该层的声波速度.重复第2步与第3步,直到某一预定深度时结束.通过反演试验,对地震振幅比例的改变,子波变形,以及第1层速度和密度的误差对反演方法的稳定性及其精度的影响进行了分析.还通过实际海洋地震勘探数据的反演试验,对这一方法在地震勘探中的应用前景作了论述.     

用逻辑树方法估计地震年发生率的不确定性

逻辑树方法可用于定量评定震级-频度关系,推断低震级地震年发生率时产生的不确定性．这些不确定性因素包括:震级-频度关系对自相似性的偏离;不同的拟合方法;确定年平均发生率的方法;拟合时的震级步长; 起算震级及震级误差等．本文以鲜水河震源区为例,研究强地震序列外推地震危险性评定所需要的4级以上地震年发生率时可能产生的不确定性．敏感性分析表明,同一种拟合方法中最敏感的因子是确定年发生率的方法,其次是震级步长和起算震级,震级误差影响相对较小． 逻辑树评定震级-频度关系不确定性的步骤为: ① 根据研究区已有资料和认识建立表述震级-频度关系不确定性的逻辑树;② 计算所有末端分支的震级-频度关系;③ 分析各不确定性因子的敏感性,修正逻辑树的结构和调整初始权重;④ 再计算末端分支的震级-频度关系,计算每一震级档的余累积分布函数; ⑤ 给定分位数,求得与此分位水准相应的地震年频数． 本文以80％和20％的分位水准,得到鲜水河震源区4级地震（不含余震）的年频数为0.643 0,而与分位数50％相应的年频数为0.631 8,两者十分接近．      

地震勘探信号时频分析方法对比与应用分析

针对地震勘探信号,对非平稳信号时频分析几种适效方法:短时Fourier变换、小波变换、S变换、Wigner分布、平滑伪Wigner分布、锥形核时频分布、AOK(adaptive optimum kernel,自适应最优核函数)分布等进行对比与应用研究.在阐明各种方法基本原理的基础上,进行数值分析与应用研究.首先对非平稳地震勘探模拟信号进行试算及时频属性提取,结合各类方法的信号表示理论,在时频局部化的精度和交叉项抑制等方面对计算结果进行对比分析;进一步应用于实际二维地震数据,提取瞬时频率和瞬时带宽等时频属性,进行比较研究.研究表明:对于地震勘探信号,就本文涉及的几种时频分析方法而言,AOK分布是时频局部化精度最高、交叉项抑制最好、时频匹配最优的方法,值得在地震勘探信号分析和地震属性提取、频谱分解等应用中深入研究和应用.     

地震勘探中加速度信号转化为速度信号的分析及应用

陆上地震勘探中多采用速度型检波器接收,随着单点高密度地震勘探的开展,加速度型检波器也被日趋广泛的应用.这两种检波器接收的地震信号具有不同的响应特征,在地震勘探中经常会遇到这样的情况:两个相邻区块分别应用速度型检波器和加速度型检波器接收,需要将这两个区块的地震资料进行连片处理,两类检波器接收资料如何处理拼接?本文在理论上分析了加速度检波器测量的加速度信号和速度检波器测量的速度信号之间的数学关系;并通过数值数据模拟分析、实际的地震数据资料分析证明了加速度信号通过积分计算以及去除积分误差后转化为速度信号的可行性和正确性,并进行了相应的对比和验证.在不同类型检波器接收地震资料的拼接处理过程中应用这种转换方法进行资料处理,得到一致性很好的地震剖面.     

利用地震P波方位属性评价瓦斯富集区

以往研究表明煤层具有弱各向异性,为正交各向异性介质,因此利用P波的方位各向异性属性来进行煤层裂隙裂缝检测、预测煤层瓦斯富集带是行之有效的方法.本文对山西阳泉新景佛洼区的三维地震资料进行处理,得到4个方位的偏移数据体,并提取与煤层裂隙裂缝相关的多种地震属性.通过计算得到裂隙裂缝的相对密度,进而对煤层瓦斯的富集带进行预测.     

一种新的地震属性:基于灰—消除趋势波动法的分形标度指数（英文）

地震属性技术广泛应用于油气勘探与开发,但因地震波在地球介质中传播的复杂性、地震数据采集观测系统的局限性以及噪音干扰等原因,地震属性用于地震资料解释仍然具有很大的不确定性,而地震属性的抗噪性直接影响着地震解释的可信度。灰色理论用于处理时间序列信号具有弱化随机性、增加规律性的能力。消除趋势波动法能有效地消除数据中的各种不确定性未知趋势本文将灰色理论和消除趋势波动法相结合,提出了基于灰色消除趋势波动法计算分形标度指数地震属性的新方法。论文以Weierstrass函数产生的非线性时间序列以及实际地震资料增加随机噪音为例,讨论了灰色-消除趋势波动法计算分形标度指数的抗噪能力,研究表明新方法计算的分形标度指数具有良好的抗噪能力。将该方法应用于四川东北部地区三维地震叠后偏移数据,与传统的消除趋势波动法提取的分形标度指数地震属性进行对比,结果表明新方法计算的分形标度指数地震属性与钻井等资料结合分析,能够更好地区分沉积相带的差异,揭示生物礁的发育带。     

基于逆时偏移的井间地震角道集提取方法

井间地震角道集是地震资料偏移成像和地震资料解释的重要资料,在井间地震成像速度建模中有重要用途.针对井间地震观测系统的特殊性,研究了井间地震角道集提取方法.考虑到井间地震记录中包含上行反射波和下行反射波,从逆时偏移方程出发,推导了井间地震观测系统下的局部偏移距成像点道集提取公式,进而给出了角度计算公式和角道集提取流程.基于模型数据测试在速度正确时得到了拉平的共成像点道集,测试了速度误差对井间共成像点道集的影响,实际资料试处理也得到了较好的效果,表明本方法是合理可行的,为后续的速度分析等地震数据处理与解释工作奠定了基础.     

基于地震属性的煤层厚度预测模型及其应用

地震属性技术在岩性和构造解释等方面得到了越来越广泛的应用,特别是在煤、油气资源勘探中具有重要的作用.基于淮南矿区谢桥1区13 1煤层地震勘探资料,提取了28种地震属性数据;通过地震属性的分析,优选出平均峰值振幅、振幅的峰态、最大绝对振幅、瞬时频率斜率等4种地震属性作为煤层厚度预测模型基本参数,结合已知钻孔资料,利用多元多项式回归以及BP人工神经网络方法,求出了各属性与煤厚之间的多元多项式回归模型及人工神经网络模型,并对模型进行了误差分析和应用结果对比分析,反映出人工神经网络模型在煤厚预测中具有好的应用效果.     

紫坪铺水库诱发地震与构造地震波谱时频特征差异性研究

利用紫坪铺水库2004～2007年记录的实际地震数据,通过比较多种滤波方法,发展出一种既能有效滤除地震数据中本底噪音、又对有效信号损伤较小的比值滤波方法。对去噪后的水库诱发地震和构造地震数据,利用小波变换方法,分别进行时频分析获得两类地震的时频谱,并提取可以反映地震发生前后能量分布和聚集情况的时频属性,最后根据各个属性的响应效果,组合得到新的属性。将时频谱与时频属性相结合,总结归纳紫坪铺水库诱发地震与构造地震波谱时频特征的差异性,为地震监测分析提供基础依据。     

AVO截距和梯度属性的多波联合提取方法

AVO截距和梯度及其交会图是岩性变化识别和油气异常检测的有效手段.但从纵波地震资料中提取的梯度值精确度不高,这是由于梯度属性中含有横波速度和密度成份.为此,本文提出联合使用纵波和转换横波来提取截距和梯度属性方法.首先,本文基于Aki-Richards建立的转换横波反射系数近似式,将上下层纵波速度、横波速度、密度参数进行等量代换推导出以AVO截距和梯度为参数的转换横波反射系数近似式.然后,基于多波多分量地震资料实现了纵波和转换横波联合提取AVO截距和梯度属性参数.层状模型测试表明联合使用纵波和转换横波比仅用纵波提取的AVO梯度参数精度高,并且对背景噪声有一定的压制而获得更稳定的梯度属性结果.多分量实际资料提取的截距属性可指示目标层的流体信息,而通过截距和梯度差属性显示目标层的流体呈局部富集,更符合测试工区的流体分布状况.     

横波速度预测方法

准确的横波测井速度是叠前地震反演和叠前地震属性分析的必要参数,然而实际生产中往往缺乏横波速度信息.采用经验公式往往精度有限,本文采用Biot-Gassmann低频速度模型,以及Pride公式建立起基质弹性模量与骨架弹性模量关系,证明了纵波速度大小随固结系数的增大而减小,因而可以通过迭代方式计算出合适的固结系数,进而得到横波速度大小.通过两个实例说明该方法能得到很高的预测精度.     

基于相对走时计算的地面微地震监测SET定位算法

基于绝对走时计算的地面微地震监测地震发射层析成像(Seismic emission tomography,SET)定位算法的定位精度受地层速度模型的准确性影响很大,速度模型中的误差会使定位结果出现偏差.本文提出了一种基于相对走时计算的地面微地震监测SET定位算法,以克服绝对走时计算对速度模型精度的高度依赖.首先,通过分析走时误差推导了相对走时计算公式.然后,利用该公式对射孔初至进行校正,得到射孔点附近其他待监测网格点的走时.通过数值模拟和实际资料测试,结果表明在基于SET定位算法的前提下,本文提出的相对走时计算方法能得到比绝对走时计算方法更准确的定位结果.     

地震几何属性的快速算法实现

地震几何属性被广泛用于提取地震数据中的几何结构特征,从而辅助解释相关的地质沉积和构造过程.本文提出基于多种递归滤波和构造导向滤波的地震几何属性快速算法,能显著提高地震相干和体曲率的计算效率和分辨率.递归滤波的计算效率远高于传统的加权求和计算,且其计算成本与平均窗口的大小无关;同时,利用高维滤波的可分离特性可以将其分解为多个一维滤波来进一步提升计算效率并有利于多线程并行运算.此外,使用构造导向滤波实现相干计算,可以有效地消除倾斜构造对结果的影响而无需传统方法中耗时的插值运算.三维实际地震数据的应用结果表明,本文提出的快速算法能将传统的相干和曲率属性计算速度提高10～30倍且对断层的刻画更加完整和连续.     

瞬变电磁频-时转换混合优化算法研究

瞬变电磁正演的常规算法是把频率域响应转化到时间域,通常有三种方法:Guptasma算法、余弦变换折线逼近算法和Gaver-Stehfest逆拉普拉斯变换算法.在计算精度上,余弦变换折线逼近算法难以求出符合精度的感应电动势,误差集中在晚期,但是当频率范围取值合理时,可以得到精度较高的磁感应强度.Guptasma算法和Gaver-Stehfest逆拉普拉斯变换算法在晚期通常有较大误差,而早期的感应电动势和磁感应强度精度都较高.在计算速度上,Guptasma算法最快,折线逼近算法由于需要计算大量频点导致速度很慢.因此,为了能够既快速又有效的完成频率域到时间域的转换,为反演奠定基础,本文以电性源为例,设计出一种新的混合算法,该算法在早期采用Guptasma算法,在误差较大的晚期采用余弦变换多项式近似算法,首先求出磁感应强度,然后采用数值差分的方法求得感应电动势.通过对Guptasma算法、余弦变换折线逼近算法以及混合优化算法的最终结果的误差和运算时间进行比较,说明混合优化算法充分发挥了Guptasma算法和余弦变换折线逼近算法的各自的优势,计算速度较快且精度较高,能够取得较好的效果.