专利名称：地震勘探烃类检测方法
地震勘探烃类检测方法，是一种利用数字地震记录，通过适当的数字处理来直接找油找气的方法。
关于这种方法，大庆石油管理局地球物理勘探公司，已经获得二项专利，专利号分别为87107616和89105573.8。这二项专利基本内容和公式为将数字磁带记录，在地震处理机上，做基本处理，使地震记录道变成做过动静校正的道、迭加道、选偏道。在同条测线上，对上述三种道任选一种或分选后靠近炮点道，按专利号87107616公式处理。
在第i个采样点，设 其中，k＝0，1，2，……，n。n是小于道的采样点个数的正整数。
按专利号为89105573.8的公式处理在第i个采样点，设 其中，k＝0，1，2，……，n。n是小于道的采样点个数的正整数。
在所选道的时窗内，对每个采样点都计算出λ′(t)。对所有的λ′(t)≥1的值分成等级，每一级用一个机器符号代表，形成符号道，把符号道输出出来就是烃类检测图。
按专利号为89105573.8的另一个公式处理在第i个采样做 与 的符号比较，不妨设 其中，K＝0，1，2，……，n。n是小于道的采样点个数的正整数。
取λ″(ti)的符号。若λ″(ti)＜0，用地震机上任一符号表示；λ″(ti)≥0用空白表示。对每个采样点都这样做，形成符号道，把符号道输出出来，也形成烃类检测图。
本发明仍根据专利号为87107616，89105573.8的“地震勘探烃类检测方法”的理论依据，提出另一“地震勘探烃类检测方法”的申请。
根据方法的理论依据，可画出如附
图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，f3＝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＝12，f2＝15，f3＝55，f4＝60，地震道被这样躲开10HZ附近的低频的宽频滤波器滤波后d(t)变为d′高(t)。在带通滤波B(f1，f2，f3，f4)中，选如f1＝5，f2＝7，f3＝55，f4＝60，地震道被这样宽频带滤波器滤波后d(t)变为d′全(t)。
新方法的理论分析。从示意图3和图4可看出在ti附近，d′低(t)与d′高(t)值的符号相反。在ti附近，设λ-(t)＝|d′低(t)-d′高(t)|＝|d′低(t)|+|d′高(t)|则有λ-(t)＞|d′低(t)|，λ-(t)＞|d′高(t)|，λ-(t)＞|d′全(t)|。
在ti附近，设λ+(t)＝|d′低(t)+d′高(t)|则有 λ+(t)小于|d′低(t)|，|d′高(t)|中的至少一个值。如果，λ+(t)＞|d′高(t)|又λ+(t)＞|d′低(t)|，则在ti附近的层肯定不是油气层。
同样，在ti附近，设λ-(t)＝|a1d′低(t)-a2d′高(t)|λ+(t)＝|a3d′低(t)+a4d′(t)|a1，a2，a3，a4是实数。先取ai(1≤i≤4)为某实数，然后对其他三个系数根据ai进行调整确定，使λ-(t)与λ-(t)，λ+(t)与λ+(t)具有同样性质。在ti附近，用λ-(t)与λ+(t)进行比较，能把油气层做为异常带找出来。
在公式中加系数的办法，能使方法更具有普遍性，充分反映方法的思想。另外，对少数油气层显示不太好的，可调整系数a1，a2，a3，a4，使油气层显示的更明显，方便于应用。
d′高(ti-j)是选的所有道被躲开10HZ附近低频的宽频带滤波器滤波后的第i-j个采样点的值。
d′全(ti-j)是选的所有道，被包含地震记录全部有效频率的全频滤波器滤波后的第i-j个采样点的值。
d′低(ti-j)是选的所有道被低频带通滤波器滤波后的第i-j个采样点的值。
d′(ti-j)是选的所有道，被另一种带通滤波器滤波后的第i-j个采样点的值。这种带通滤波器分二种情况，使d′(ti-j)也分两种情况当λ+(t)式中a3＝0时，d′(ti-j)是用所选道的全部有效频率的带通滤波器对所选道滤波后的第i-j个采样点的值，等于d′全(ti-j)；当a3≠0时，d′(ti-j)是用不包含油气层影响所选道的低频带而包含所选道其余全部有效波频率的带通滤波器对所选道滤波后的第i-j个采样点的值，等于d′高(ti-j)。
上面所谈到的几种滤波器B(f1，f2，f3，f4)中的频率参数的选择，只需根据要求能把主要频带包括进去，就能计算出较好的烃类检测图。因此，具体频率参数的选择范围较大。它主要与储油层的岩性和野外采集接收频率有关。参照已知井，在处理机上做实验能较好地选取频带参数。
根据以上分析，该方法可这样描述。本方法的内容是将野外磁带记录，在地震处理机中，做基本处理，使地震记录道变成做过动静校正的道，迭加道，迭偏道。在同一条线上，对上述三种道任选一种或分选后靠近炮点道，做如下处理，处理公式为在第i个采样点，设λ′-(ti-j)＝a1d′低(ti-j)-a2d′高(ti-j)λ′+(ti-j)＝a3d′低(ti-j)+a4d′(ti-j)λ-1(ti)=Σj=0±K|λ--′(ti-j)|---λ-2(ti)=|Σj=0±Kλ--′(ti-j)|]]>λ-3(ti)=Σj=0±K|λ-+′(ti-j)|---λ-4(ti)=|Σj=0±Kλ-+′(ti-j)|]]>λ1(ti)=Σj=0±K|λ--′(ti-j)|/|λ-+′(ti-j)|]]>λ2(ti)＝λ1(ti)/λ3(ti) λ3(ti)＝λ1(ti)/λ4(ti)λ4(ti)＝λ2(ti)/λ3(ti) λ3(ti)＝λ2(ti)/λ4(ti)其中，K＝0，1，2，……，n。n是小于道的采样点个数的正整数。
(1)若λj(ti)≥b成立(1≤j≤5，b＞0，b是常数)，则在ti采样点保留原振幅值(用其他非零值也可)。否则，把ti采样点的值用零代替。同样，也可反之替代ti采样点的值。b的选择，根据系数a1，a2，a3，a4和已知井的烃类显示符合程度而定。
在选定的时窗内，对每道采样点都这样做。然后，用波形或波形加面积办法显示出来，形成烃类检测图。也可对每道在所选的时窗内，对每个采样点都计算λi(t)(1≤i≤5)。在某一采样点λi(t)＜b，此采样点用零值代替；对λi(t)≥b的λi(t)值分成等级，每个等级用一个机器符号替代采样点的值。在所选时窗内每道都这样做，形成符号道，把符号道输出出来，也形成烃类检测图。
(2)对(1)中所述，当n＝0时效果最好，此时λj(ti)＝|λ′-(ti)|/|λ′+(ti)| (1≤j≤5)其他处理和显示办法与(1)中所述相同。
(3)对公式λj(ti)(1≤j≤5)可取倒数，取倒数后的公式易从(1)和(2)中看出判断油气层的条件和显示办法，这里不再论述。
如何预测油气带。在烃类检测图上，除了把油气层显示出来外，还把储水层按有效厚度显示出来。用频率扫描的办法，因油气与水的频率适应性不同，把水区分开。在这方面大庆油田第九采油厂地质所有很成功的经验。
如何看检测图。图的上方所示是X方向CDP点数。纵向是t0时间。根据已知井，已知储油气带，烃类检测所示储油气层，符合的最佳，图的显示最强，符号显示最密集的为最佳。
权利要求
1.一种利用地震记录，通过适当的数字处理来直接找油找气的方法。将野外磁带记录，在地震处理机上，做基本处理，使地震记录道变成做过动静校正的道、迭加道、选偏道。其特征在于在同条测线上，对上述三种道任选一种或分选后靠近炮点道，做如下处理，处理公式为在第i个采样点，设λ′-(ti-j)＝a1d′低(ti-j)-a2d′高(ti-j)λ′+(ti-j)＝a3d′低(ti-j)+a4d′(ti-j)λ-1(ti)=Σj=0±K|λ--′(ti-j)|---λ-2(ti)=|Σj=0±Kλ--′(ti-j)|]]>λ-3(ti)=Σj=0±K|λ-+′(ti-j)|---λ-4λ(ti)=|Σj=0±Kλ-+′(ti-j)|]]>λ1(ti)=Σj=0±K|λ--′(ti-j)|/|λ-+′(ti-j)|]]>λ2(ti)＝λ1(ti)/λ3(ti) λ3(ti)＝λ1(ti)/λ4(ti)λ4(ti)＝λ2(ti)/λ3(ti) λ3(ti)＝λ2(ti)/λ4(ti)其中，K＝0，1，2，……，n.n是小于道的采样点个数的正整数。a1，a2，a3，a4是实数。先取ai(1≤i≤4)为某实数，然后对其他三个系数根据ai进行调整确定，使|λ′-(ti-j)|与|d′低(ti-j)-d′高(ti-j)|，|λ′+(ti-j)|与|d′低(ti-j)+d′高(ti-j)|具有同样性质。d′高(ti-j)是选的所有道被躲开10HZ附近低频的宽频带滤波器滤波后的第i-j个采样点的值。d′全(ti-j)是选的所有道，被包含地震记录全部有效频率的全频滤波器滤波后的第i-j个采样点的值。d′低(ti-j)是选的所有道被低频带通滤波器滤波后的第i-j个采样点的值。d′(ti-j)是选的所有道，被另一种带通滤波器滤波后的第i-j个采样点的值。这种带通滤波器分二种情况，使d′(ti-j)也分两种情况当λ4(t)式中a3＝0时，d′(ti-j)是用所选道的全部有效频率的带通滤波器对所选道滤波后的第i-j个采样点的值，等于d′全(ti-j)；当a3≠0时，d′(ti-j)是用不包含油气层影响所选道的低频带而包含所选道其余全部有效波频率的带通滤波器对所选道滤波后的第i-j个采样点的值，等于d′高(ti-j)。上面所谈到的几种滤波器B(f1，f2，f3，f4)中的频率参数的选择，只需根据要求能把主要频带包括进去，就能计算出较好的烃类检测图。因此，具体频率参数的选择范围较大。它主要与储油层的岩性和野外采集接收频率有关。参照已知井，在处理机上做实验能较好地选取频带参数。若λj(ti)≥b成立(1≤j≤5，b＞0，b是常数)，则在ti采样点保留原振幅值(用其他非零值也可)。否则，把ti采样点的值用零代替。同样，也可反之替代ti采样点的值。b的选择，根据系数a1，a2，a3，a4和已知井的烃类显示符合程度而定。在选定的时窗内，对每道采样点都这样做。然后，用波形或波形加面积办法显示出来，形成烃类检测图。也可对每道在所选的时窗内，对每个采样点都计算λi(t)(1≤i≤5)。在某一采样点λi(t)＜b，此采样点用零值代替；对λi(t)≥b的λi(t)值分成等级，每个等级用一个机器符号替代采样点的值。在所选时窗内每道都这样做，形成符号道，把符号道输出出来，也形成烃类检测图。
2.据权利要求1所述的地震勘探烃类检测方法，其特征在于当n＝0时效果最好，此时λj(ti)＝|λ′-(ti)|/|λ′+(ti)| (1≤j≤5)其他处理和显示办法与1.中所述相同。
3.对公式λj(ti)(1≤j≤5)可取倒数，取倒数后的公式判断油气层的条件和显示办法，易从1.和2.中看出。如何预测油气带。在烃类检测图上，除了把油气层显示出来外，还把储水层按有效厚度显示出来。用频率扫描的办法，因油气与水的频率适应性不同，把水区分开。在这方面大庆油田第九采油厂地质所有很成功的经验。如何看检测图。图的上方所示是X方向CDP点数。纵向是t0时间。根据已知井，已知储油气带，烃类检测所示储油气层，符合的最佳，图的显示最强，符号显示最密集的为最佳。
全文摘要
地震勘探烃类检测方法，是一种用野外地震记录来直接找油气的方法。利用不同带通滤波器，对做过基本处理的各种道滤波。在相同采样点上对不同滤波后的值加权系数做某种计算和比较。在选定的时窗内对所有道都这样做，把地下油气带的位置找出并显示出来，形成烃类检测图。此图能把相对迭加剖面上储油气带整个形态显示出来。
文档编号G01V1/00GK1136668SQ9610265
公开日1996年11月27日 申请日期1996年1月26日 优先权日1996年1月26日
发明者黄克有 申请人:大庆石油管理局地球物理勘探公司