一种基于等地质时间的致密薄储层预测方法
【技术领域】
[0001]本发明涉及石油地球物理勘探领域，尤指一种基于等地质时间的致密薄储层预测方法。
【背景技术】
[0002]随着地震勘探程度的不断深入，勘探目标逐步转向薄互层低孔、低渗致密储层，并成为油气勘探的重要领域，因此，薄互层致密储层预测技术成为国内外各油田的重要研究课题。一般把厚度小于1/4波长的储层称为薄储层，目前薄储层地震预测方法有两类，一类是采用叠后地震属性分析、叠后地震资料反演、可视化等技术手段；另一类是采用叠前属性、叠前反演等技术手段。
[0003]由于致密薄储层与围岩波阻抗差异小，有效储层预测难度大，而且气层厚度相对较薄(<10m)，物性相对较差(渗透率〈2 X 10 3 μ m2)，储层横向变化大，地震预测难度大，储层预测方法技术手段已经从叠后转向叠前，甚至多波技术。同时由于解释较薄的目的储层层位非常困难，利用上述两种方法进行预测时，层位基本上采用的是漂移层位。这种情况下，有些地方“穿时”会比较严重，那么利用“穿时”的层位约束的储层预测结果，会存在比较大的误差。
【发明内容】
[0004]为解决致密薄储层的预测问题，克服约束层位“穿时”存在的缺陷，本发明基于叠前多波资料(PP波+PS波地震资料)，利用相对等地质时间层位的约束，进行叠前PP+PS同时反演，以实现对致密薄储层进行有效储层预测。
[0005]本发明提出的基于等地质时间的致密薄储层预测方法，包括:步骤1，对包含PP波及PS波地震数据进行保幅处理，获取叠前道集的PP波分偏移距数据和PS波分偏移距数据；步骤2，联合全波列测井资料，分别对步骤1获取的PP波分偏移距数据和PS波分偏移距数据的大套标志层层位进行解释，获得PP波的大套标志层层位解释数据和PS波的大套标志层层位解释数据；步骤3，对步骤2获得的PP波的大套标志层层位解释数据和PS波的大套标志层层位解释数据进行反射等时性判断，获得判断结果，并进行层位反复修改，获得PP波和PS波的大套标志层层位等地质时间反射层位；步骤4，采用步骤3获得的等地质时间反射层位内插获得等地质时间的PP波目的层薄层反射层位和PS波目的层薄层反射层位；步骤5，利用步骤4获得的等地质时间的PP波目的层薄层反射层位约束，基于步骤1中PP波分偏移距数据进行叠前反演得到纵、横波速度比数据体，根据所述纵、横波速度比数据体建立初始走时数据体；步骤6，根据步骤2获得的PP波的大套标志层层位解释数据和PS波的大套标志层层位解释数据分别对所述PP波分偏移距数据和PS波分偏移距数据进行数据对齐，并将误差分配到所述初始走时数据体上，获得最终走时数据体，利用最终走时数据体进一步将PP波分偏移距数据和PS波分偏移距数据进行数据对齐；步骤7，在步骤4获得的等地质时间的PP波目的层薄层反射层位和PS波目的层薄层反射层位的约束下，对所述PP波分偏移距数据和PS波分偏移距数据进行叠前联合反演，并配合岩石物理解释量板进行预测，获得致密薄储层的预测结果。
[0006]进一步的，在步骤1中，至少获取3个以上的叠前道集的PP波分偏移距数据和PS波分偏移距数据。
[0007]进一步的，在步骤3中，对步骤2获得的PP波的大套标志层层位解释数据和PS波的大套标志层层位解释数据进行反射等时性判断，获得判断结果，并进行层位反复修改，获得PP波和PS波的大套标志层层位等地质时间反射层位，包括:步骤31，对步骤2获得的PP波的大套标志层层位解释数据和PS波的大套标志层层位解释数据分别进行小波分频，对小波分频后的数据进行地震倾角估算，获得高、低频的倾角；步骤32，当高、低频的倾角差值高于或等于阈值时，对层位进行修改，并重新计算高、低频的倾角，当高、低频的倾角差值低于阈值时，判定地震反射轴是等时的，并在等时性分析的控制下，获得PP波和PS波的大套标志层层位等地质时间反射层位。
[0008]本发明的基于等地质时间的致密薄储层预测方法可以在地震PP+PS分偏移距资料的基础上，引入地震等时性分析技术，解决目的层薄层的层位“穿时”问题，并且利用等地质时间的层位的约束，进行叠前PP+PS同时反演，可以进一步获得精度和分辨率较高的弹性参数体，结合岩石物理解释量板进行解释，从而可以提高地震薄储层预测的能力和钻井预测成功率，特别适用于致密气田薄层储层的预测。
【附图说明】
[0009]此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中:
[0010]图1为本发明一实施例的基于等地质时间的致密薄储层预测方法的流程图。
[0011]图2为本发明一具体实施例的PP、PS分偏移距数据的大套标志层层位解释结果示意图。
[0012]图3为本发明一具体实施例的PP波分偏移距剖面与倾角差体叠合示意图。
[0013]图4为本发明一具体实施例的在等时性分析的基础上内插得到的目的层等地质时间层位示意图。
[0014]图5A及图5B分别为本发明一具体实施例的叠前PP+PS联合反演获得的弹性参数结果及父会结果不意图。
【具体实施方式】
[0015]以下配合图示及本发明的较佳实施例，进一步阐述本发明为达成预定发明目的所采取的技术手段。
[0016]图1为本发明一实施例的基于等地质时间的致密薄储层预测方法的流程图。如图1所示，该方法包括:
[0017]1、一种基于等地质时间的致密薄储层预测方法，其特征在于，该方法包括:
[0018]步骤1，对包含PP波及PS波地震数据进行保幅处理，获取叠前道集的PP波分偏移距数据和PS波分偏移距数据。其中，至少获取n(n ^ 3)个叠前道集的PP波(纵波)分偏移距数据和η (η ^ 3)个叠前道集的PS波(转换波)分偏移距数据。
[0019]步骤2，联合全波列测井资料，分别对步骤1获取的PP波分偏移距数据和PS波分偏移距数据的大套标志层层位进行解释，获得PP波的大套标志层层位解释数据和PS波的大套标志层层位解释数据。
[0020]结合图2所示，图2为一具体实施例的PP、PS分偏移距数据的大套标志层层位解释结果示意图。利用全波列测井资料，分别标定PP和PS分偏移距数据，并分别解释地震反射轴相对连续较好的易于追踪的大套标志层层位分别为q5_PP、tc2_PP和q5_PS、tc2_PS。[0021 ] 步骤3，对步骤2获得的PP波的大套标志层层位解释数据和PS波的大套标志层层位解释数据进行反射等时性判断，获得判断结果，并进行层位反复修改，获得PP波和PS波的大套标志层层位等地质时间反射层位。
[0022]具体的，步骤3包括:
[0023]步骤31，对步骤2获得的PP波的大套标志层层位解释数据和PS波的大套标志层层位解释数据分别进行小波分频，对小波分频后的数据进行地震倾角估算，获得高、低频的倾角；
[0024]步骤32，当高、低频的倾角差值高于或等于阈值时，对层位进行修改，并重新计算高、低频的倾角，当高、低频的倾角差值低于阈值时，判定地震反射轴是等时的。结合图3所示，图3为一具体实施例的PP波分偏移距剖面与倾角差体叠合示意图。从图中可以看出q5_PP、tc2_PP倾角差较小，可以认为q5_PP