[发明专利]一种应用于陆相盆地储层反演的曲线频率融合方法

摘要

本发明公开一种应用于陆相盆地储层反演的曲线频率融合方法，包括如下步骤：（a）基础曲线预处理；（b）地震合成记录，建立时深关系曲线；（c）曲线频谱分析及滤波处理；（d）多频率曲线融合；（e）基于频率融合曲线的地震反演。本发明提出了一种基于频率融合的曲线融合方法，在陆相盆地储层储层反演中根据测井曲线的不同频率尺度反映不同的地质与沉积特征的特点，利用频率融合的方法建立拟波阻抗测井曲线。将曲线频率融合成果应用于地震反演，解决了地震反演精细储层刻画的难题，有利于提高地震反演储层预测的精度和效率。

主权项

1.一种应用于陆相盆地储层反演的曲线频率融合方法，包括如下步骤：(a)基础曲线预处理；(b)地震合成记录，建立时深关系曲线；(c)曲线频谱分析及滤波处理；(d)多频率曲线融合；(e)基于频率融合曲线的地震反演；其中，所述步骤(a)包括如下内容：(a1)对声波时差曲线开展标准化处理，选取全区稳定发育的泥岩段作为标准层，计算标准井标准层段的测井曲线平均值同时计算每口井目的层段测井曲线平均值求出校正差，从测井曲线中减去校正量，即完成校正，校正公式为：其中，AC_标为标准化后测井曲线值；AC为标准化前测井曲线值；为目的层段测井曲线平均值；为标准井标准层段的测井曲线平均值；(a2)对自然电位曲线开展基线校正，选取全区稳定发育、测井响应具有一致性的稳定的泥岩段作为标准层，读取标准层自然电位曲线平均值作为自然电位曲线基线漂移的基值，选取目的层内各个泥岩层深度及其对自然电位幅值，按照自然电位曲线采样率进行垂向上线性插值，生成校正基线曲线SP基，对自然电位曲线进行基线校正，利用公式：其中，为第k个采样点处基线校正后的自然电位幅度值；Sp^k为第k个采样点处基线校正前的自然电位幅度值；为第k个采样点处校正基准曲线幅度值；SP_标准为标准泥岩段自然电位平均幅值；(a3)将步骤(a2)校正结果SP_NEW曲线归一化到波阻抗曲线量纲范围内，根据读取波阻抗曲线开展频率直方图统计，得出波阻抗曲线的最大幅值最小幅值由于自然电位幅值极性与波阻抗曲线幅值极性相反，对步骤(a2)校正结果SP_NEW曲线开展反极性求取SP_‑NEW＝SP_NEW，并对SP_‑NEW曲线进行频率直方图统计，统计SP_‑NEW曲线最大值和最小值将SP_NEW曲线量纲到波阻抗曲线的范围之内，均一化公式如下：其中，为第k个采样点处均一化后的自然电位幅度值；为第k个采样点处均一化前的自然电位幅度值；和为波阻抗曲线的最大幅值和最小幅值；和为井待校正自然电位曲线的最大幅值和最小幅值；(a4)引用步骤(a3)所读取的波阻抗量曲线最大幅值和最小幅值对步骤(a2)校正结果SP_NEW曲线进行频率直方图统计，得出电阻率曲线最大值和最小值将电阻率曲线量纲到波阻抗曲线的范围之内，均一化公式如下：其中，为第k个采样点处均一化后的电阻率幅度值；为第k个采样点处均一化前的电阻率幅度值；和为波阻抗曲线的最大幅值和最小幅值；和为井待校正电阻率曲线的最大幅值和最小幅值；通过步骤(a3)和步骤(a4)将完成自然电位曲线和电阻率曲线的归一化处理，将自然电位曲线和电阻率曲线归一到波阻抗曲线量纲范围，以此达到各井曲线的一致性和可对比性；所述步骤(b)包括如下内容：(b1)制作地震合成记录，建立时深关系曲线，利用步骤(a1)标准化的声波时差曲线和密度曲线，根据公式：其中，R^k为第k个采样点处换算反射系数值；ρ^k为第k个采样点处密度曲线幅度值；v^k为第k个采样点处速度曲线幅值，其大小为换算出反射系数后，并选用合适频率的雷克子波δ，将反射系数曲线与子波褶积得到合成记录：其中，Xi为第i个采样点合成记录值；δk为第k个采样点处子波值；Ri为第i个采样点处反射系数；M1、M2分别为子波的起始样点和末尾样点序号；以合成记录为基础，建立地震数据的时间域与测井深度域之间的对应关系；(b2)应用步骤(b1)所建立地震数据的时间域与测井深度域之间的对应关系，将步骤(a3)中所得到SPlog曲线转换为地震反射时间与自然电位幅值相对应的SPtime曲线；(b3)应用步骤(b1)所建立地震数据的时间域与测井深度域之间的对应关系，将步骤(a3)中所得到RLLDlog曲线转换为地震反射时间与电阻率幅值相对应的RLLDtime曲线；所述步骤(c)包括如下内容：(c1)对步骤(b1)中的得到SPtime曲线进行快速傅里叶变换(FFT)，得到SPtime曲线频率谱H(f)SP，利用低通滤波器公式：H(f)SPlow＝C(f)×H(f)SP其中，C(f)为低通滤波函数；f为频率；fL为低通滤波器的截止频率；H(f)SP为自然电位曲线频谱函数；H(f)SPlow为自然电位曲线低频分量频谱函数；通过低通滤波处理得到自然电位曲线低频分量频谱函数H(f)SPlow；(c2)对步骤(b2)中的得到RLLDlog曲线进行快速傅里叶变换(FFT)得到RLLDtime曲线频谱H(f)RLLD，根据高通滤波器公式：H(f)RLLDhigh＝C(f)×H(f)RLLD其中，C(f)为高通滤波函数；f为频率；fH为高通滤波器的截止频率；H(f)RLLD为自然电位曲线频谱函数；H(f)RLLDhigh为电阻率曲线高频分量频谱函数；通过高通滤波处理得到电阻率曲线高频分量频谱函数H(f)RLLDhigh；截止频率fL和fH的选取为地震反演所应用到的地震数据的有效频带的最高值；所述步骤(d)包括如下内容：通过融合技术将步骤(c1)得到的自然电位曲线低频分量频谱函数H(f)SPlow和步骤(c2)得到的电阻率曲线高频分量频谱函数H(f)RLLDhigh进行融合生成频谱曲线H(f)Ip；针对频谱曲线H(f)Ip开展反傅里叶变换，即得到基于频率融合拟波阻抗曲线Ip拟合；所述步骤(e)包括：将步骤(d)所得频率融合拟波阻抗曲线Ip拟合应用于地震反演运算过程。