[发明专利]测井信号驱动的地震记录有色补偿处理方法

摘要

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

主权项

1.一种测井信号驱动的地震记录有色补偿处理方法，其特征在于：实施步骤如下：步骤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，最大值α_m=1.0，扫描增量Δα=0.1，β起始值β₀=0.0，最大值β_m=1.0，扫描增量Δβ=0.1。计算所有由α_i=α₀+iΔα和β_k=β₀+kΔβ所确定的有色补偿算子的振幅谱C_ik(f)，采用的算法为：Cik(f)=||1-αie-j2πf1-βke-j2πf||]]>其中，C_ik(f)是有色补偿算子的振幅谱，无量纲，f是频率，单位是赫兹，j是单位虚数，||||表示复数取模运算。计算有色补偿算子振幅谱C_ij(f)与反射系数振幅谱R(f)的误差e_ij，采用的算法为：e_ij=∫[C_ij(f)-R(f)]²df比较所有误差e_ij的大小，最小误差所对应的颜色参数α和β即为反射系数的实际颜色参数。由反射系数振幅谱估算颜色参数α和β在数学上属于典型的非线性最优化问题，该最优化问题很容易掉进局部极值问题的陷进。本发明所述的计算方法有效地避免了局部极值问题，能够得到可靠的实际反射系数颜色参数。步骤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、将有色补偿之后的地震记录绘制为能够反映地下结构特征的地震剖面图像。