专利名称：地震勘探烃类检测方法
地震勘探烃类检测方法，是一种利用数字地震记录，通过适当的数字处理来直接找油找气的方法。
关于这种方法，大庆石油管理局地球物理勘探公司，在1987年11月1日申请一个名称为“地震勘探烃类检测方法”的专利，申请号为87107616。
专利申请号为87107616的“地震勘探烃类检测方法”的基本内容和公式为将野外数字磁带记录，在地震处理机上，做基本处理，使地震记录道变成做过动静校正的道、迭加道、迭偏道。在同一条测线上，对上述三种道任选一种或分选后靠近炮点道，按如下公式处理在第i个采样点，设
其中，n＝2k，k＝0，1，2，……。
在所选道的时窗内，对每个采样点都计算出λ(t)。对所有的λ(t)≥1的值分成等级，每一级用一个机器符号代表，形成符号道，把符号道输出出来就是烃类检测图。但该方法不是唯一的方法，根据该方法的理论依据，还存在另外的处理方法。
本发明仍根据专利申请号为87107616的“地震勘探烃类检测方法”的理论依据，提出另一“地震勘探烃类检测方法”。
根据方法的理论依据，可画出如附
图1、2、3、4的示意图。设地下有一储油气水平单层，它的反射模型如图1所示。(射线〔1〕，〔2〕，〔3〕都是垂直层面并重合在一起的，为说明方便，而把它们分开画出。)入射波方向〔1〕，波形函数b(t);层顶反射波方向〔2〕，波形函数b1(t);层底反射波方向〔3〕，波形函数b2(t)。b(t)、b1(t)、b2(t)的频谱为B(f);B1(f);B2(f)。
图2是水平单层顶、底的反射波形。〔5〕是b2(t)的形状，〔7〕是b1(t)的形状，〔6〕是采样点ti的位置。
设此层反射记录道主要是b1(t)、b2(t)组合的复合振动。反射记录为它们的代数和。则d(t)＝b1(t)+b2(t)。
假设，在做地震处理时，在带通滤波B(f1，f2，f3，f4)中，选f1＝5，f2＝7，f3＝12，f4＝15时，地震道被这样低频带通滤波器滤波后d(t)变为d′低(t)，b1(t)被滤波后变为b′1(t)，b2(t)被滤波后变为b′2(t)。图3中〔8〕，〔10〕分别是b′2(t)，b′1(t)的波形。〔9〕是采样点ti的位置。在带通滤波B(f1，f2，f3，f4)中，选f1＝5，f2＝7，f3＝55，f4＝60，地震道被这样宽频带滤波器滤波后d(t)变为d′全(t)。
假设，在地震资料处理时，在带通滤波B(f1，f2，f3，f4)中，选f1＝12，f2＝15，f3＝55，f4＝60时，地震道被这样躲开10HZ附近的低频的宽频带滤波器滤波后d(t)变为d′高(t)。同时，b1(t)被滤波后变为b″1(t)，b2(t)被滤波后变为b″2(t)。
图4中，〔11〕，〔13〕分别为b″2(t)，b″1(t)的波形。〔12〕是采样点ti的位置。
新方法的理论分析。B1(f)的10HZ附近的低频能量不会明显地增加;而B2(f)的10HZ附近的低频能量，相对于B1(f)的将增加很多。这样，d(t)被躲开10HZ附近低频的带通滤波器滤波后b2(t)明显地被压缩变成b″2(t);而b1(t)将保留相当强的振幅，变成b″1(t)。在ti附近将有｜d′高(t)｜＝｜b″1(t)+b″2(t)｜＞｜d(t)｜我们的讨论是针对地震有效波的，d(t)就应是d′全(t)，即｜d′高(t)｜＞｜d′全(t)｜，在ti附近。
设λ′(t)＝｜d′高(t)｜/｜d′全(t)｜设在ti附近计算λ′(t)时，分母小，分子变大，λ′(t)更大。λ′(t)能灵敏地反映地震波穿过油气层时d′全(t)的特点。
另外，从图3，图4中看出d′低(t)与d′高(t)比较，在ti附近d′低(t)，d′高(t)值的符号将相反。
设λ″(t)＝d′高(t)/d′低(t)取λ″(t)的符号。在储有油气的地层中，λ″(t)取负数的机会很多;否则，λ″(t)取负数的机会少。
根据以上分析，该地震烃类检测方法可这样描述。
本方法的内容是将野外数字磁带记录，在地震机上做基本处理，使地震记录道变成做过动静校正的道，迭加道，迭偏道。在同条测线上，对上述三种道任选一种或分选后靠近炮点道。
对上述的道，在每个采样点上，可用两种公式处理。
(1).在第i个采样点，设
……(D)其中，K＝0，1，2，……，n。n是小于道的采样点个数的正整数。
d′高(ti-j)是选的所有道被躲开10HZ附近低频的宽带滤波器滤波后的第i-j个采样点的值。
例如，对迭加道，在带通滤波为B(f1，f2，f3，f4)中，选f1＝12，f2＝15，f3＝55，f4＝60，对所选的每道做滤波。
d′全(ti-j)是选的所有道，在带通滤波器B(f1，f2，f3，f4)中，选如f1＝5，f2＝7，f3＝55，f4＝60，地震道被这样宽频带滤波器滤波后的第i-j个采样点的值。
在每个采样点计算λ′(t)值。把λ′(t)值用不等间隔，从数值1开始，按数值大小分成数十个等级，用机器符号代表每一级。把这些符号道输出出来，就形成烃类检测图。
(2).方法的另外一个公式，是在第i个采样上做
的符号比较，不妨设
……(E)其中，K＝0，1，2，……，n。n是小于采样点个数的正整数。
取λ″(ti)的符号。若λ″(ti)＜0，用地震机上任一符号表示;λ″(ti)＞0，用空白表示。对每个采样点都这样做，形成符号道，把符号道输出出来，也形成烃类检测图。
d′低(ti-j)是选的所有道被低频带通滤波器滤波后的第i-j个采样点的值例如，对迭加道，在带通滤波B(f1，f2，f3，f4)中，选f1＝6，f2＝8，f3＝12，f4＝14，对每个道做滤波。
(1)和(2)所述的公式，实践表明，当n＝0时，效果最好。此时，在第i个采样点，(D)式变成λ′(ti)＝｜d′高(ti)｜/｜d′全(ti)｜
然后，把λ(t
)按(1)中所述处理。
当n＝0时，在第i个采样点，(E)式变成λ″(ti)＝d′高(ti)/d′低(ti)然后，把在每个采样点计算出的λ″(t
)按(2)所述方法处理，得烃类检测图。
本发明的优点本发明的方法，在任何地震机上均可实现，公式简单易做，除野外数字地震记录外，不需要其他附加条件，把相对迭加剖面上储油气带成面积的显示出来。
看检测图6，7。从右到左是CDP点数增加方向。CDP点数在图的上方显示出来。时间在图的右边显示。
用此方法已经做了三条过井测线，都有较好的符合。
图5是齐家176.0测线，CDP点数500到639的水平迭加剖面。A处是金6井的位置。CDP点数为556。
图6是176.0测线，CDP点数从492到608的烃类检测图。竖直线过CDP556。本图用新方法(1)做的。两箭头所夹的带为油气显示。
图7是176.0测线，CDP点从450到623烃类检测图。竖直穿过556。本图是用专利申请号为87107616所提出的方法做的。两箭头所夹的带为油气显示。
从图5看，金6井的产油层是箭头所指的强层上面一层。在井处to时间是1.235秒支1.28秒之间。产油井段(加上油水合层)共76米，试油日产38吨，正常生产估计为6吨。
补正
权利要求
1.一种利用数字地震记录资料，直接找油气的地震勘探烃类检测方法，先将野外数字磁带记录，在地震机上做基本处理，使地震记录道变成做过动静校正后的道，迭加道，迭偏道。其特征在于在同条测线上，任选上述三种道的一种或分选后靠近炮点道，按如下公式处理在第i个采样点，设
其中，K=0，1，2，……，n。n是小于每道采样点个数的正整数。(d′高)是选的所有道被躲开10HZ附近低频的宽带通滤波器滤波后的第i-j个采样点的值；例如，对迭加道，在带通滤波B(f，f2，f3，f4)中，选f1=12，f2=15，f3=55，f4=60，对每个道做滤波。(d′全)是选的所有道被宽频带通滤波器滤波后的第i-j个采样点的值例如，对迭加值，在带通滤波B(f1，f2，f3，f4)中，选f1=6，f2=8，f3=55，f4=60，对每个道做滤波。把λ′(t)值用不等间隔，从数1开始，按数值大小分成数十个等级。用机器内符号代表每一级，把这些符号道都绘图显示出来，就形成烃类检测图。
2.据权利要求1所述的地震勘探烃类检测方法，其特征在于在同条测线上，任选权利要求1所述的四种道的一种，按如下办法处理在每个采样点上，做
的符号比较。例如，在第i个采样点上，设
其中，k＝0，1，2，……，n。n是小于采样点个数的正整数。d′低(ti-j)是选的所有道被低频带通滤波器滤波后的第i-j个采样点的值例如，对迭加道，在带通滤波B(f1，f2，f3，f4)中，选f1＝6，f2＝8，f3＝12，f4＝14，对每个道做滤波。看λ″(ti)的符号。如果λ″(ti)＜0，用机器内任何符号代替λ″(ti)的负符号。如果λ″(ti)＞0，用空白表示，形成符号道，把这样符号道输出出来，也得烃类检测图。
3.椐权利要求1所述的地震勘探烃类检测方法，其特征在于当n＝0时效果最好，此时按如下公式处理在第i个采样点，公式为λ′(ti)＝｜d′高(ti)｜/｜d′全(ti)｜
4.据权利要求2的地震勘探烃类检测方法，其特征在于当n＝0时效果最好，此时按如下方式处理在第i个采样点，做d′高(ti)与d′低(ti)的符号比较。设λ″(ti)＝d′高(ti)/d′低(ti)
全文摘要
地震勘探烃类检测方法，是一种用数字地震记条来直接找油气的方法。按一种公式，在每个采样点上计算出一种油气存在的地震特性参数，用符号表示不同参数值。把符号道输出出来，形成烃类检测图；按另外一种公式计算，若采样点附近地层界面的反射系数有极性反转，用一个机器符号表示，否则为空白。把这样符号道输出出来，也形成烃类检测图。这二种图都带有统计作用，把相对迭加剖面上储油气带成面积的显示出来。
文档编号G01V1/28GK1049411SQ8910557
公开日1991年2月20日 申请日期1989年8月10日 优先权日1989年8月10日
发明者黄克有 申请人:大庆石油管理局地球物理勘探公司