三维叠后地震数据可纵向显示的处理方法

在三维地震资料处理中,为使叠后相邻测线的地震数据在绘图纸上能够纵向同时显示,达到节约绘图纸、降低成本的目的,本文提出了一种利用地震资料处理模块生成可纵向显示的三维地震绘图数据的处理方法.该方法通过进行叠后三维数据体的分离、延长测线数据的记录长度、修改线号和数据体舍并等处理工作,产生可纵向显示的绘图数据体.经实际绘图显示,不同测线的地震剖面由原来的横向单张显示,转为了纵向的多张绘图.因此,大大节约了纸张,降低了地震资料处理的成本.     

断层识别技术及其在MB油气田的应用

结合MB凸起油气勘查项目,总结断层解释方法:既要依据地震时间剖面特征来分析判断,更要了解盆地范围内应力场、变形机理以及工区范围内的构造样式.针对微小断层识别困难这一问题,采取多种技术手段,包括:(1)叠后目标处理提高地震资料对小断层的分辨能力;(2)正演模型分析指导小断层的解释;(3)相干数据体水平切片和近反射道几何特征分析识别微小断层;(4)三维数据体任意线剖面检验微小断层,识别了常规地震剖面难以确定的小断层,达到了断层精细解释的目的.     

基于MATLAB的煤矿TEM数据体三维可视化技术

针对煤矿地质勘探中TEM面积性探测获取的三维数据体,基于Matlab语言对其三维可视化技术进行研究与实践.阐述煤矿TEM三维数据体格式及其可视化要求,提出‘立体图’、‘切片图’和‘等值面图’三种综合性可视化方式,论述使用Matlab语言调用Surfer库函数进行二次开发实现TEM测线数据的批量网格化,解决了煤矿TEM数据体的不规则性问题.同时论述三维可视化技术中关于表面图、三维曲线图、等值线、等值面等三维图形的关键绘制技术与应用函数,实现规则TEM数据体的基本图形表达.在以上工作基础上,作者开发出相应的三维可视化应用软件,通过实例展示其丰富、实用的显示功能.该软件为煤矿TEM资料的处理与解释工作提供了新手段,帮助提高对TEM解释结果的认知能力和资料解释水平.     

基于稀疏反演三维表面多次波压制方法

三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3DSRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.     

拉东投影法三维叠前深度偏移

对地下地质构造进行正确成像是地震勘探的最终目的,由于三维地震资料采集不可能都沿垂直构造走向的方向进行,为地震资料的三维处理带来了许多困难. 本文将三维叠后拉东投影偏移思想应用于三维叠前处理中,提出了三维叠前投影偏移算法. 利用拉东投影变换的原理,将整个三维叠前数据体投影到一系列各方向的径向线上,各方位角的构造都包含在其中某条或多条径向剖面上. 投影完成后,形成一系列的独立的二维叠前测线,可采用二维叠前深度偏移成像方法来实现各径向线的叠前偏移,当各径向剖面偏移完成后,在时间切片上进行反投影,从而最终形成三维叠前深度偏移结果. 实际应用表明,用本方法进行三维叠前深度偏移时,深度偏移剖面对横向分辨率有所提高,对陡地层和小断层的成像效果有所改善.     

三维多方向边缘探测技术识别盐丘的新方法（英文）

对3D地震数据进行盐丘构造的精确识别在地震数据分析应用中非常重要。本文提出了对三维偏移地震数据进行盐丘构造边界识别的新方法。本算法克服了目前边界识别技术中的不足,即二维数据只考虑x方向(测线方向)和y方向(垂直测线方向)的边界问题,三维数据只考虑x方向(测线方向)、y方向(垂直测线方向)和z方向(时间)的边界问题。算法结合了沿对角线方向和沿x、y、z方向计算得出的三维梯度,以达到对盐丘边界的精确探测。沿x、y、z和沿对角线方向的边界的联合检测保证了盐丘边界倾面表现为弱反射时本方法仍能获得较好的效果。不同与其它盐丘的边界及构造检测技术,本文提出的算法不依赖地震数据中的振幅变化。将本方法应用至公开的荷兰海洋F3号区块数据测试中,结果证明,本文方法对盐丘边界的探测精度高于单独使用基于梯度或基于结构的方法。更重要的是,算法计算效率高,可以满足实时实现及调度。     

基于数据规则化和稀疏反演的三维表面多次波压制方法

表面多次波是海洋地震勘探中的主要问题.目前,二维数据驱动的表面多次波压制技术(SRME)已经比较成熟,并且已经成为工业界压制海洋表面多次波的主流方法.但是由于二维SRME算法没有考虑横测线方向上多次波的贡献,导致在处理实际三维海洋资料时存在比较大的误差.将二维SRME算法扩展到三维空间后可以得到三维SRME算法,但是由于目前实际采集的三维海洋资料的观测系统存在拖缆漂移,而且横测线方向采样过于稀疏,直接应用三维SRME算法无法准确预测表面多次波.本文提出的通过数据规则化配合稀疏反演的三维表面多次波压制方法能够解决这种实际资料和三维SRME算法之间的矛盾.本文通过研究数据规则化与反规则化技术,使得数据分布满足三维SRME的要求;通过研究稀疏反演技术,有效解决了横测线方向采样稀疏对于多次波预测的影响,三维实际海洋资料的应用结果验证了方法的有效性和可行性.     

三维OBS探测实验中炮点和OBS位置校正新方法

炮点和海底地震仪(ocean bottom seismometer,OBS)位置校正是三维地震数据处理的基本环节,也是获取高精度三维速度结构的关键所在.本文基于南海洋陆转换带(Continental-Oceanic-Transition zone,COT)IODP367/368钻探区开展的三维OBS深地震探测数据,开展了炮点及OBS位置校正研究,新的校正方法主要体现在三个方面:(1)利用连续三个炮点的平均航向对中间炮点进行位置校正,更真实地反映气枪枪阵与船体之间的软连接状态;(2)根据"滑动窗口"思想将海水声学速度阈值划分成N等份,通过循环测试获得全局最优的OBS位置校正结果,改进了前人方法只能获得局部最优解的问题;(3)针对单条测线穿过的OBS,通过加入其临近测线的直达水波走时,构成视双测线OBS位置校正法,提高了校正精度.49台OBS位置校正结果表明,除3台单测线法校正的OBS在垂直测线方向存在较大不确定性外,其余误差范围均为35m左右.本文改进的OBS位置校正方法,不仅提高了单条测线穿过的OBS位置校正精度,保证后续三维地震结构研究的可靠性,而且为今后类似的OBS位置校正提供了经验和借鉴.     

基于GIS的当雄同震InSAR形变场分析

InSAR技术对同震形变的量测达到了厘米级的精度,但在数据后处理和结果分析方面仍然存在明显的不足。采用可对空间信息进行存储、管理及分析等功能的GIS,对InSAR数据获取的2008年10月6日16时30分西藏自治区当雄县MW6.3地震的同震形变场进行分析。结果表明:(1)GIS可对多源数据进行有效管理;(2)可以确定当雄地震的震中位置;(3)可获取沉降区的最大形变量;(4)可确定主要的形变区间;(5)可将形变结果进行三维展示。GIS可有效地弥补InSAR数据后处理、数据分析及成果展示方面的不足。     

源于地震沉积学的高精度层序地层学应用：以路易斯安那近海老虎滩地区50号中新统弗米利恩区块为例

介绍了应用钻井和三维地震资料来划分高精度(四级)层序和体系域的地震沉积学方法。其关键技术包括:(1)将地震数据进行90°相位调整来记录岩性,从而获得更完善的测井资料;(2)假设并解释沉积体系中连续的平面地貌特征。下面我们推荐一种新的解释方法,就是将高精度层序地层学的研究重点从解释垂直地震剖面转变为分析更加水平、分辨率更高的地震地貌信息。这方面的研究表明,有岩性限制的地震资料中,地层切片可提供通常为同沉积体系的连续地震图像。反过来,这种图像又可作为识别和划分高频体系域、层序边界以及层序地质年代范围的基础。在路易斯安那近海的中新统地层中,运用测井资料划分的四级层序或层序组,可通过地震资料来划分,地震数据体的主频为30Hz,相当于分辨率为30ft(10m)。所用的这一分辨率在井控外的地震覆盖区,对高精度层序地层格架进行精确再现非常有效。     

地震数据采集器非线性指标测试方法及其影响分析

为了分析总谐波失真度和互调失真度这两个非线性指标的测试方法、输入信号幅值等变化对其测试结果的影响,本文从总谐波失真信号、互调失真信号的产生原理出发,详细地给出了两种指标的测试方法、数据处理方法及相关影响分析,并通过实际应用中不同类型的输入信号进行了总谐波失真度和互调失真度试验。测试结果表明:信号模式和信号幅值对总谐波失真度的影响显著;输入信号幅值对互调失真度的影响较大;输入信号的频率变化对两指标的测试结果无影响。本研究进一步细化了地震数据采集器性能指标的标准化测试方法,该方法可广泛应用于地震数据采集器的非线性指标测试。      

基于弹性波全波形反演的主被动源多分量混采地震数据速度建模

在常规地震采集中,被动源地震波场往往被视为噪声而去除,这就造成了部分有用信息的丢失.在目标区进行主动源和被动源弹性波地震数据的多分量混合采集,并对两种数据进行联合应用,使其在照明和频带上优势互补,能显著提高成像和反演的质量.本文针对两种不同类型的主被动源混采地震数据,分别提出了相应的联合全波形反演方法.首先,针对主动源与瞬态被动源弹性波混采地震数据,为充分利用被动源对深部照明的优势,同时有效压制被动震源点附近的成像异常值,提出了基于动态随机组合的弹性波被动源照明补偿反演策略.然后,针对低频缺失主动源与背景噪声型被动源弹性波混采地震数据,为充分利用被动源波场携带的低频信息,并避免对被动源的定位和子波估计,提出了基于地震干涉与不依赖子波算法的弹性波主被动源串联反演策略.最后,分别将两种方法在Marmousi模型上进行反演测试.结果说明,综合利用主动源和被动源弹性波混采地震数据,不仅能增强深部弹性参数反演效果,还能更好地构建弹性参数模型的宏观结构,并有助于缓解常规弹性波全波形反演的跳周问题.     

基于SOM和PSO的非监督地震相分析技术

地震相分析技术是储层预测的一种重要方法,可以用来描述有利沉积相带的分布规律.传统的地震相聚类分析方法对大数据的处理运算速度较慢,且容易陷入局部极小值,造成聚类分析的结构不准确.本文提出基于自组织神经网络(SOM)和粒子群优化方法(PSO)相结合的地震相分析技术,利用自组织神经网络能够保持原始地震数据的拓扑结构特性的特点,将大量冗余样本压缩为小样本数据,再通过粒子群的全局寻优能力改善K均值聚类的效果.理论模型和实际应用表明该方法能既有效实现数据压缩,又能提供较为准确的全局解,在地震相预测中兼顾计算效率和计算精度.     

多震源地震数据偏移成像方法

多震源地震技术是一种高效的地震数据采集方法技术,得到的地震记录是来自多个震源的混合地震数据.本文在多震源波场传播理论和地震波场满足线性叠加原理的基础上,提出了两种多震源地震数据的偏移成像方法.第一种方法是首先对多震源地震数据进行分离,得到各个单震源的地震数据,然后再利用常规的偏移成像方法进行处理;第二种方法是多震源地震数据的直接偏移成像.把本文提出的多震源偏移成像方法应用于数值模拟的多震源地震数据,验证了本文方法的正确性和有效性,直接偏移成像方法较分离后再偏移方法具有更高的计算效率.     

基于噪声谱比法的福建台网观测台站的场地响应研究

应用噪声谱比法, 收集了福建台网观测台站记录的地脉动噪声数据, 计算分析了85个观测台站的场地响应． 根据场地响应曲线形态将福建台网观测台站分为4类, 并简要分析了影响台站场地响应结果的因素． 计算结果表明, 福建台网观测台站台基情况普遍较好, 在台站观测频带内场地响应曲线平坦, 大部分台站均没有明显的频率放大点． 本文采用H/V谱比法与Moya方法所获得的37个台站的场地响应结果基本一致, 仅有少数台站存在较大差异． 应用场地响应分析结果, 结合计算震级时所量取的周期, 对福建台网观测台站的单台震级校正值进行了估计, 并与统计方法获得的单台震级校正值进行了比对． 结果表明, 两种方法获得的台站震级校正值存在一定的相关性, 但由于考虑的影响因素不同, 具体的台站震级校正值尚存在差别．      

主动源与被动源地震数据插值及联合数据成像

本文提出了两种情况下主动源数据和被动源数据的插值方法,并研究了两种数据在偏移成像中的互补效果.基于互相关法被动源数据重构原理,本文提出了结合多域迭代去噪技术的重构方法.提出了两种时间域主动源和被动源数据的插值方法,分别是共炮点域能量匹配插值和共检波点域最小平方匹配插值.然后对获得的主动源和被动源联合地震数据进行叠前深度偏移成像.在被动源活跃度不是很高的地区进行被动源地震勘探时,少量的主动源地震数据可以有效控制和补充被动源数据的成像效果.在稀疏炮点的主动源勘探中,有效利用被动源的信息能够在成像中增加更多的细节信息,提高成像质量.     

预测多次波的逆散射级数方法与SRME方法及比较

逆散射级数方法和SRME方法是典型的两类基于波动方程的无需地下结构信息的自由表面多次波预测方法,本文详细讨论了这两类方法的基本思想、方法原理和实现思路,并对2维逆散射级数方法进行了降维简化,推导给出了1.5维逆散射级数衰减自由表面多次波的方法,减少了计算量,并降低了对3维规则观测系统的要求.利用经典SMAART模型数据测试比较了2维逆散射级数方法、1.5维逆散射级数方法及2维SRME方法预测自由表面多次波的效果及优缺点,比较了三者在方法实现、并行计算效率及对数据要求上的异同,并结合分析比较结果,给出了实际应用中方法选择的建议.     

Ground roll attenuation using a time-frequency dependent polarization filter based on the S transform

The ground roll and body wave usually show significant differences in arrival time, frequency content, and polarization characteristics, and conventional polarization filters that operate in either the time or frequency domain cannot consider all these elements. Therefore, we have developed a time-frequency dependent polarization filter based on the S transform to attenuate the ground roll in seismic records. Our approach adopts the complex coefficients of the S transform of the multi-component seismic data to estimate the local polarization attributes and utilizes the estimated attributes to construct the filter function. In this study, we select the S transform to design this polarization filter because its scalable window length can ensure the same number of cycles of a Fourier sinusoid, thereby rendering more precise estimation of local polarization attributes. The results of applying our approach in synthetic and real data examples demonstrate that the proposed polarization filter can effectively attenuate the ground roll and successfully preserve the body wave.     

基于互信息量的地震信号检测和初至提取方法

本文将现代信息科学中的重要基础理论——信息论引入到地震信号的分析和处理中,提出了一种利用互信息量检测地震信号并进行初至提取的方法,给出了该方法的基本原理和具体算法.计算结果表明,互信息量随滑动窗口移动的变化情况可以明确标示出待测地震信号的存在,通过提取互信息量的突变点或峰值点,能够有效拾取地震波的初至时刻.这种方法还可以特别针对非传统的地震信号,如连续随机信号、随机脉冲序列等进行处理.在这些信号条件下,往往信噪比很低,波形随传播距离发生较大变化,传统的信号检测方法效果不佳,而这种以随机信号检测为基础的互信息量探测方法具有很好的性能,很有发展前景.     

利用单台波形资料计算大柴旦6．3级地震余震震源机制

以2008年11月10日大柴旦Ms6．3级地震余震数字地震波观测资料为基础,运用矩张量反演方法,利用单台三分量波形资料计算余震的震源机制。收集整理大柴旦Ms6．3级地震公认的震源参数结果和青海区域地震台网数字地震波资料,以大柴旦Ms6．3级地震的震源参数作为约束,计算理论地震波来确定该区域的地壳速度模型。以此模型计算ML≥2．0级余震震源机制解。共收集160个余震资料,最后确定了85个余震震源机制解,其中最大震级为朋。4．7,最小震级为ML2.0。进行单台与多台联合矩张量反演结果对比,两种方法所确定的结果相差不大。实际计算表明,在观测资料信噪比较高的情况下,运用单台数字地震波观测资料可以较好地确定中小地震的震源机制。