测井信号驱动的地震记录有色补偿处理方法技术

测井信号驱动的地震记录有色补偿处理方法，应用于油气勘探地震资料处理。特征：人工激发并记录地震信号；对地震记录进行噪声压制、反褶积、速度分析、偏移等常规地震信号处理；采集声波测井信号和密度测井信号，计算反射系数序列的振幅谱；利用反射系数序列的振幅谱确定反射系数的颜色参数；对地震记录进行有色补偿；将有色补偿之后的地震记录绘制为能够反映地下结构特征的地震剖面图像。效果是：补偿白噪反射系数假设对地震记录反褶积的影响，提高地震记录的分辨率，增强了地震记录反映地下结构的能力。

【技术实现步骤摘要】

本专利技术涉及油气地球物理勘探
，特别涉及地震资料处理，是一种。
技术介绍
地震勘探是一种利用人工地震技术探测地下构造的勘探方法。它按照一定的方式在地表附近激发地震波，产生称之为地震子波的地震信号，地震子波由震源开始向地下传播，遇到地质界面之后，在界面处发生透射和反射，透射的地震子波继续向下传播，而反射子波在界面位置向上传播，其强度和极性依赖于界面反射系数的大小和符号。来自不同深度界面的地震子波以不同的时间到达地表，通过布设在地表的一种称之为检波器的接收装置，接收来自不同深度地质界面的反射波，接收到的数字信号称为地震记录。地震子波具有一定的延续时间，当地层较薄时，来自薄层底面和顶面的反射子波相互干涉，影响了地震信号对薄层的分辨能力。现有的地震勘探理论认为，当薄层的厚度小于地震信号四分之一波长时，在地震记录不能有效分辨。反褶积技术是目前提高地震信号分辨率最为主要的方法，它通过提高地震信号的有效频率，压缩子波的延续时间，减小地震波长，弱化薄层厚度与地震波长的矛盾，增强地震信号反映薄层结构的能力。反褶积方法属于数学上的欠定问题，为此，现有的反褶积方法需要两个基本的数学假设，一是假设地震子波的相位为最小相位，二是假设反射系数序列为白噪，。Walden等(1988)在对不同地区反射系数序列进行统计分析后发现，反射系数序列的频谱并非白谱，整体趋势表现为低频弱、高频强的蓝谱。Rosa等(1991)进一步发现实际反射系数的频谱特征可以模拟为指数为O. 5-1. 5的指数函数。Okaya (1995)基于反射系数的频谱特征对薄储层的分布规律进行了预测和描述。李国发(2009)就反射系数白噪假设对反褶积的影响进行了深入分析，指出了白噪假设对高分辨率地震资料处理的潜在风险。为克服白噪假设对地震信号分辨率的影响，需要对现有反褶积方法处理后的地震记录进行有色补偿，目前最为常用的有色补偿算法为，y (z) =c (Z) χ (Z)其中，z是Z变换符号，无量纲，X(Z)是地震记录的Z变换，无量纲，y (Z)是有色补偿之后地震记录的Z变换，无量纲，C(Z)是有色补偿算子的Z变换，无量纲，其具体表达式为，本文档来自技高网...

【技术保护点】
一种测井信号驱动的地震记录有色补偿处理方法，其特征在于：实施步骤如下：步骤A、人工激发地震波，利用检波器探测来自地下的反射地震信号，形成地震记录并记录在磁带上。步骤B、从磁带上读取地震记录，对地震记录进行噪声压制、反褶积、速度分析、偏移等常规地震资料处理，得到常规处理之后的地震记录并输出到步骤E所在的处理单元。由于常规处理方法隐含有反射系数为白噪的数学假设，与实际地层反射系数的颜色特征不一致，需要利用测井信号对实际反射系数的颜色特征进行模拟分析，根据测井分析的结果对常规处理的地震记录进行有色补偿。步骤C、利用测井信号分析得到实际地层反射系数的振幅谱，具体步骤为：步骤C1、采集并记录声波测井信号和密度测井信号，将其数字化后存储在磁盘上。步骤C2、从磁盘上读取声波测井信号和密度测井信号，计算地层波阻抗并输入到步骤C3所在的处理单元，采用的算法为：a(i)=ρ(i)p(i)Δz其中，a(i)是波阻抗，单位是（克/立方厘米）×（米/秒），p(i)为声波测井信号，单位是秒，ρ(i)为密度测井信号，单位是克/立方厘米，Δz为深度采样间隔，单位是米，i是深度采样序号。步骤C3、接收步骤C2输出的波阻抗信号，计算实际地层反射系数并输出到步骤C4所在的处理单元，采用的公式为：r(i)=a(i)-a(i-1)a(i)+a(i-1)式中，r(i)为反射系数，无量纲。步骤C4、接收步骤C3输出的反射系数，计算反射系数的振幅谱并输出到步骤D所在的处理单元。计算反射系数振幅谱之前，需要将反射系数由深度域转换到时间域，由于速度随深度是变化的，深度域等间隔采样的反射系数变换到时间域之后，其采样间隔是不相等的，需要进行等间隔重新采样，由此引发的假频问题会严重影响反射系数振幅谱的估算精度。为此，本专利技术采用下式在不进行时深转换和重采样的情况下直接计算反射系数的振幅谱，避免了采样假频对振幅谱的影响。R(f)=||Σi=1mr(i)e-j2πft(i)||式中，R(f)为反射系数的振幅谱，无量纲，j是单位虚数，||||表示复数取模运算，m是反射系数采样个数，f是频率，单位是赫兹，t(i)是利用公式由声波测井信号p(i)得到的不同深度的反射时间，单位是秒。步骤D、接收步骤C4输出的实际反射系数的振幅谱，对其进行统计模拟，得到描述实际反射系数颜色特征的两个无量纲参数α和β。具体过程为：设定α起始值α0=0.0，最大值αm=1.0，扫描增量Δα=0.1，β起始值β0=0.0，最大值βm=1.0，扫描增量Δβ=0.1。计算所有由αi=α0+iΔα和βk=β0+kΔβ所确定的有色补偿算子的振幅谱Cik(f)，采用的算法为：Cik(f)=||1-αie-j2πf1-βke-j2πf||其中，Cik(f)是有色补偿算子的振幅谱，无量纲，f是频率，单位是赫兹，j是单位虚数，||||表示复数取模运算。计算有色补偿算子振幅谱Cij(f)与反射系数振幅谱R(f)的误差eij，采用的算法为：eij=∫[Cij(f)?R(f)]2df比较所有误差eij的大小，最小误差所对应的颜色参数α和β即为反射系数的实际颜色参数。由反射系数振幅谱估算颜色参数α和β在数学上属于典型的非线性最优化问题，该最优化问题很容易掉进局部极值问题的陷进。本专利技术所述的计算方法有效地避免了局部极值问题，能够得到可靠的实际反射系数颜色参数。步骤E、接收步骤B输出的地震记录，利用步骤D得到的反射系数颜色参数α和β，对地震记录进行有色补偿，得到有色补偿之后的地震记录。采用的算法为：y(Δt)=x(Δt),i=1y(iΔt)=x(iΔt)-αx[(i-1)Δt]+βy[(i-1)Δt],1<i≤N其中，x是补偿之前的地震记录，y是补偿之后的地震记录，若采用压电检波器接收地震信号，地震记录的单位是帕，若采用速度检波器接收地震信号，则地震记录的单位是米/秒，Δt是时间采样间隔，单位是毫秒，N是时间采样个数。步骤F、将有色补偿之后的地震记录绘制为能够反映地下结构特征的地震剖面图像。FDA00002406388100022.jpg...

【技术特征摘要】
1.一种测井信号驱动的地震记录有色补偿处理方法，其特征在于实施步骤如下步骤A、人工激发地震波，利用检波器探测来自地下的反射地震信号，形成地震记录并记录在磁带上。步骤B、从磁带上读取地震记录，对地震记录进行噪声压制、反褶积、速度分析、偏移等常规地震资料处理，得到常规处理之后的地震记录并输出到步骤E所在的处理单元。由于常规处理方法隐含有反射系数为白噪的数学假设，与实际地层反射系数的颜色特征不一致，需要利用测井信号对实际反射系数的颜色特征进行模拟分析，根据测井分析的结果对常规处理的地震记录进行有色补偿。步骤C、利用测井信号分析得到实际地层反射系数的振幅谱，具体步骤为 步骤Cl、采集并记录声波测井信号和密度测井信号，将其数字化后存储在磁盘上。步骤C2、从磁盘上读取声波测井信号和密度测井信号，计算地层波阻抗并输入到步骤C3所在的处理单元，采用的算法为 a =Az P{i) 其中，a⑴是波阻抗，单位是(克/立方厘米)X (米/秒)，p⑴为声波测井信号，单位是秒，P (i)为密度测井信号，单位是克/立方厘米，Az为深度采样间隔，单位是米，i是深度米样序号。步骤C3、接收步骤C2输出的波阻抗信号，计算实际地层反射系数并输出到步骤C4所在的处理单元，采用的公式为 r(/) = —H- 、)a(i) + a(j-l) 式中，r(i)...

【专利技术属性】
技术研发人员：李国发，周新源，王招明，彭更新，杨海军，梁向豪，李皓，段文胜，满益志，赵锐锐，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：