基于Hartley变换的归一化总梯度法

归一化总梯度法对下半空间的重力总梯度模进行归一化计算,从而利用极值来获得地质体的分布,但现有计算方法存在计算精度较低的问题,并且计算结果不太稳定。本文采用Hartley变换,并引入新的余弦滤波因子,实现重力归一化总梯度,来压制高频干扰;此外,针对计算过程中下延深度越大干扰越严重的弊端,将余弦滤波因子的幂次设计为随深度变化的函数,从而获得更加稳定的计算结果。模型试验和在玲珑金矿区采空区的实际数据结果表明,所提出的重力归一化总梯度法具有计算速度快、精度高、稳定性强等优点。     

利用重力归一化总梯度及相位法研究断裂构造

用波数域的位场转换和快速傅立叶变换理论研究重力归一化总梯度和归一化相位;在波数域向下延拓和导数计算的滤波算子中,分别引入圆滑滤波因子,抑制了向下延拓计算中对噪声干扰的放大作用,增加了计算的稳定性;利用改进的方法,计算并阐述了4种断裂构造GH场等值线和相位曲线的变化特征;别列兹金重力归一化总梯度法中极大值理论在研究断裂构造中具有局限性,对于无限延伸断裂构造,应根据GH场等值线的走向、形状的特点,结合归一化相位曲线的转折点确定断裂空间位置.     

基于Facet模型梯度算子一致性的地震数据不连续性识别方法

存在断层、角度不整合面等不连续结构的地质体的自动识别在地震构造解释中具有重要的意义,这些地质特征的地震响应为同相轴不连续。常规的地震数据不连续性识别方法应用范围有限,参数设置依赖于人为经验,易导致识别结果辨识度差。本文将一致性作为新的不连续性识别方法引入到地震数据处理中。首先利用定位精度高、易于扩展的Facet模型梯度算子计算一致性,其次对一致性数据作阈值化处理,最后利用数学形态学中的腐蚀、膨胀及细化算法作进一步处理,实现了对地震数据不连续性信息的自动识别。经过理论和实际资料测试,并与C3相干算法和方差算法对比分析,证实了本文所提方法在地震数据不连续性识别方面具有更高的稳定性和辨识度,可以作为地层不连续性识别的有力工具。     

基于支持向量机法识别砂岩中流体类型

砂岩储层孔隙中的流体识别一直是石油勘探开发过程中重要的环节,传统方法主要依赖于测井数据,但是在测井数据缺失的条件下较难得到准确的流体识别结果.本文提出一种只依靠地震数据的砂岩中流体识别的新方法,并选择地球物理方法可测或可求的地球物理参数σ、ρλ和ρμ作为流体识别因子,然后进行模型实验.首先,设置典型流体状态,用Gassmann方程进行流体替换,将得到的流体识别因子作为支持向量机的训练集数据,并定义支持向量机的分类标签;之后,设置随机流体状态,利用Gassmann方程计算流体因子,将得到的结果作为支持向量机的测试集数据.将训练集、测试集数据集输入支持向量机,进行分类,得出测试集数据的分类结果.模型实验分类结果表明,支持向量机法可以判别砂岩孔隙中流体的主要属性.     

基于OC-Seislet变换的海洋涌浪噪声衰减方法

在海上地震资料采集过程中,涌浪噪声是一种常见的噪声干扰,这类噪声一般表现为低频、强能量、长周期,使用一般的方法很难在去噪的同时达到信号保真的效果。笔者针对该噪声的特征,选取了基于波动方程炮检距连续(offset continuation,OC)算子的OC-Seislet变换方法进行消噪处理;该方法应用OC算子来表征复杂波场,对含噪声数据进行压缩,在变换域通过软阈值处理实现信噪分离,再将结果反变换到数据域,从而达到去除涌浪噪声的目的。通过对模型数据和实际数据的处理,验证了OC-Seislet变换方法能够在去除原始数据中涌浪噪声的同时,最大程度地保护复杂构造下的地震波信息。     

基于四阶累积量的同相轴自动拾取方法

针对强随机噪声干扰情况下地震勘探资料中同相轴无法正确拾取的问题,提出了基于四阶累积量的同相轴自动拾取方法。该方法利用四阶累积量进行时间延迟估计,避免了相关高斯噪声的影响;同时在参考道的选择上,以一定的间隔设置多个参考道,减少了同相轴拾取过程中的累积误差。理论模型计算表明:该方法相对于互相关方法,对相关高斯噪声有更强的抑制能力,信噪比达到-5 dB。在实际资料应用中,拾取了与同相轴有关的参数：同相轴位置,该参数为地震资料的属性分析和精确定量解释提供了依据。     

连续缺失地震数据的高阶流式预测滤波插值方法

由于观测系统实施以及经济因素的限制,采集到的勘探地震数据在空间方向上总是不规则分布的,并且往往会出现数据的大范围连续缺失情况.许多后续地震数据处理方法(例如：多次波压制和波动方程偏移等)都需要空间上规则分布的数据.插值技术是一种解决地震数据缺失问题的有效手段,但是传统的数据插值方法在进行连续缺失数据重建时往往会出现失效的情况,尤其在处理非平稳地震同相轴时精度不高,并且大多数的方法需要迭代计算,在处理高维大规模数据时效率较低.针对连续缺失地震数据的快速插值问题,本文提出了一种非迭代的时空域高阶流式预测滤波插值方法,通过使用高阶限制条件来提高连续缺失数据的滤波器估计精度,提高局部约束条件的稳定性,改善低阶流式计算由于滤波器系数无法连续更新所造成的插值失效情况.同时,空间非因果滤波器和蛇形插值路径的设计方案可以有效减小大范围连续缺失数据和数据边界对于预测滤波器的计算误差,本方法能够有效处理包括近炮检距缺失情况在内的连续缺失数据插值重建.通过与工业标准傅里叶凸集投影(POCS)方法进行比较,理论模型和实际数据处理结果表明,本文提出的高阶流式预测滤波插值方法对高维连续缺失地震数据有较好的重建效果...     

大地电磁数据揭示中国东北地区牡丹江缝合带俯冲极性及活化特征

佳木斯地块和松嫩地块位于中亚造山带最东端,两者的汇聚受控于二叠纪至侏罗纪期间牡丹江洋的俯冲及闭合作用.牡丹江缝合带后期又受到古太平洋板块和太平洋板块俯冲作用的直接影响,因此是研究古缝合带俯冲极性和后期改造作用的理想地区.本文对跨越牡丹江缝合带的160km长大地电磁剖面数据进行三维反演,建立了缝合带及邻区的电阻率模型.结果揭示了牡丹江洋的俯冲极性和俯冲痕迹,并提供了(古)太平洋板块俯冲作用导致牡丹江缝合带活化的地电学证据.缝合带表现为复杂的高导结构,位于松嫩地块-佳木斯地块碰撞带下方的地壳尺度西倾高导体,代表牡丹江洋闭合的古俯冲带,表征了牡丹江洋西向俯冲的极性.同时,由地表地质定义的牡丹江断裂不能完全代表牡丹江缝合带的深部结构,对于该缝合带的定义应扩展至横向延伸约70km的整个高导区域.俯冲前缘弧下的固体导电矿物由深部被抬升至上地壳,其下与地幔连通的“烟囱”状高导结构代表了幔源物质的上侵通道,表明牡丹江缝合带受古太平洋板块和太平洋板块的俯冲作用影响发生活化,导致通道内原本冷却、结晶的物质发生重熔.因此,(古)太平洋板块的俯冲作用破坏了中亚造山带东段古碰撞带的岩石圈结构,其下方大尺度地壳高...     

高频电磁波在不同目标体上反射图像特征研究

针对高频电磁波技术在北京勘察设计行业应用的空白,本次研究利用探地雷达技术对北京地区典型工程目标体进行了研究.针对不同的地下目标体对高频电磁波所反映出的不同反射特性,总结出一些图像特征规律.通过建立对各类雷达探测图像进行地质解释的理论依据和基础,本次研究对层状介质、地下管道、CFG桩、碎石夯扩桩、水泥土桩和洞穴等目标体的雷达图像特征进行了较系统的研究和总结,为今后日常的检测工作奠定了很好的基础.     

松原地区地震发震机制与迁移特点研究

本研究重新拾取震相到时,对2003—2021年松原地区发生的地震进行了绝对定位,采用双差层析成像法在松原地区获得0～20 km深度的高分辨率三维V_P和V_S模型,对该区地震的发震机制、地震迁移与震源区速度结构之间的关系等进行了深入讨论.结果显示：2003—2021年松原地区地震总体沿着松原—肇东断裂呈北北东向展布.震源区(123°E—125°E,44°N—46°N,0～15 km)的泊松比变化范围为0.24～0.30,平均值为0.242,低于0.3,推断震源区地壳内部几乎不存在部分熔融.绝大多数地震发生在低P、低S波速度异常处,推测震源区下方热物质上涌所携带的流体注入到地壳后触发了松原地区地震.研究区下方5～10 km深度现存在一条连接松原市、前郭县和长岭地区,北北东走向的S波低速带,推断松原—肇东断裂为深部流体提供了通道,控制了地震呈北北东向串珠状展布,且地震随时间具有回跳的迁移特点.根据断裂弱化模型推测松原地区发生更大地震的可能性较低,但在长岭地区,小地震可能会持续发生.