一种井震联合的沉积旋回匹配分析方法
【技术领域】
[0001] 本发明属于地球物理勘探地震数据解释及综合研究领域，具体涉及一种井震联合 的沉积旋回匹配分析方法。
【背景技术】
[0002] 旋回是地层学和沉积学中表示具有重视性或周期性沉积记录的术语，包含描述性 和解释性的双重含义。地层记录的旋回性是沉积作用、沉积环境或驱动机制的不同，地层旋 回通常区分为自旋回和它旋回两类。自旋回主要由环境旋回所致，具有局域性和区域性，如 河道、三角洲朵叶体的侧向迁移形成的地层旋回。它旋回通常是由环境旋回因素以外的驱 动机制所致，具有大区域、跨地区，甚至全球可对比性，如海平面变换旋回和米兰科维奇旋 回等。在沉积旋回划分对比过程中，常常受到很多问题的制约。应用地震资料划分沉积旋 回，由于受分辨率的限制。只能实现大尺度旋回的追踪，如典型的不整合面或大的海（湖） 泛面，不能实现精细的等时对比，而测井数据是载有地层沉积旋回且分辨率更高的信号，但 其横向对比难度大。因此，需要通过井震联合的对比方法将两者结合起来，取长补短，才能 更好的进行沉积旋回的解释。井震联合对比面临的主要难题是不同尺度的地震旋回和测井 旋回该如何匹配？本发明提出了一种井震联合的沉积旋回匹配方法，能够较好解决这一难 题，有助于开展井震联合的沉积旋回对比解释研究工作。
【发明内容】
[0003]本发明的目的是针对现有技术的不足，提供一种井震联合的沉积旋回匹配分析方 法，以解决由于地震与测井数据分辨率不同造成的对应沉积旋回对比困难的问题。
[0004]本发明采用的技术方案是：
[0005] 首先通过数据的预处理得到相同采样率的地震和测井数据，在此基础上分别对两 者采用经验模态分解的方法进行处理，得到不同频谱范围的分解曲线，对分解曲线进行相 关性分析，并通过地震相位旋转调整优化相关性，最后得到一个较为满意的匹配结果。
[0006] 依据上述技术方案的包括以下步骤：
[0007] A :加载测井数据与对应井旁道地震数据，作为沉积旋回匹配分析的初始数据；
[0008] B:选择一个时深关系将测井数据转换到时间域，按照井旁道地震数据采样率将测 井数据重采样；
[0009] C:通过经验模态分解方法分别对测井和地震数据进行处理得到不同频率范围的 分析曲线；
[0010]D:对各个分解曲线进行两两相关性分析，如果得到的相关性较差，旋转地震相位 后再重复步骤C，直到得到一个比较满意的匹配结果；
[0011] E:利用得到的分解曲线，进行地层沉积旋回对比。
[0012] 上述步骤C包括下列步骤：
[0013] C1 :找出测井地震信号的局部最大最小值，将这些最大、最小值用三次样条插值函 数拟合形成信号的上下包络；
[0014] C2 :计算上下包络的平均值，将测井地震信号与平均值相减，得到一个新的信号；
[0015] C3 :重复上述步骤Cl、C2直到新的信号不存在负的局部极大值和正的局部极小 值。
[0016] 上述步骤D的相关性分析方法包括如下步骤：
[0017]D1 :分别计算井分解曲线、地震分解曲线的数学期望；
[0018] D2 :分别计算井分解曲线、地震分解曲线的方差；
[0019] D3 :计算井分解曲线与地震分解曲线的协方差；
[0020] D4 :计算井分解曲线与地震分解曲线的相关系数。
[0021] 上述步骤D的地震相位旋转方法包括如下步骤：
[0022] DD1 :对数据进行离散傅里叶变换，得到数据的实部和虚部；
[0023] DD2 :将实部数据的前半部分乘以cos(Θ)，虚部数据的前半部分乘以Sin(Θ)，其 中Θ为旋转角度；
[0024] DD3 :实部和虚部后半部分赋零值；
[0025] DD4:将实部和虚部进行离散傅里叶反变换，其中实部数据为相位旋转后的数据。
[0026] 本发明利用经验模态分解将不同频率的地震数据和测井数据进行处理，得到能够 反映不同级别沉积旋回的分解曲线，通过给定的时深关系对井震各个分解曲线进行相关性 分析，选取相关性最高的井震分解曲线可用于联合开展沉积旋回的对比划分。该方法能够 解决不同尺度资料的地层联合对比问题，可用于井震联合地层对比等勘探综合研究工作。
【附图说明】
[0027] 图1井震联合沉积旋回匹配分析流程
[0028] 图2井震数据分解曲线对比图
[0029] 图3地震相位旋转前井第4分量与地震第2分量对比图
[0030] 图4地震相位90°旋转后井第4分量与地震第2分量对比图
[0031] 图5地震道第2分量曲线地层旋回解释示意图
【具体实施方式】
[0032] 本发明的一种井震联合的沉积旋回匹配分析方法较佳实施方式为：
[0033] a.加载测井数据与对应井旁道地震数据，作为沉积旋回匹配分析的初始数据；选 择一个时深关系将测井数据转换到时间域，按照井旁道地震数据采样率将测井数据重采 样。具体步骤如下：首先将地震起始时间为开始点，查找离该点最近的对应深度点，将该深 度点对应的曲线值作为井数据的起始点，将起始时间按照地震采样率依次查找，并记录对 应深度点的曲线值，直到完成井曲线的重采样。
[0034]b.通过经验模态分解方法分别对测井和地震数据进行处理得到不同频率范围的 分析曲线。主要输入参数是样点个数、数据数组、采样率和终止条件值（SD)。具体步骤如 下：找出数据的局部最大最小值，如果极大极小值个数小于4则终止，将这些最大、最小值 用三次样条插值函数拟合形成信号的上下包络；求取上下包络的平均值并用原始数据减去 该平均值，利用积分法计算终止条件值（SD)，如果不满足终止条件则将减去包络平均值的 数据作为原始数据继续上一步操作，直到满足终止条件，得到最终分解曲线。
[0035]c.对各个分解曲线进行两两相关性分析，如果得到的相关性较差，旋转地震相位 后再重复步骤C。直到得到一个比较满意的