二维批量叠前反演方法和白云岩储层预测方法及系统技术方案

本发明专利技术提供了一种二维批量叠前反演方法和白云岩储层预测方法及系统。该二维批量叠前反演方法包括：对二维分偏移距叠加数据体进行融合处理得到分偏移距叠加体融合地震数据体；基于融合地震数据体及已钻井测井数据建立初始常速低频模型；对初始常速低频模型进行叠前反演求解；判断反演结果与已钻井测井曲线之间的误差是否在预设误差范围内；若不在，基于叠前反演结果重新修改层位；以修改后的层位为约束并利用叠前反演结果更新初始常速低频模型；对更新后的初始常速低频模型重新进行叠前反演求解更新叠前反演结果；直至反演结果与已钻井测井曲线之间的误差在预设误差范围内，输出此时的叠前反演结果。此时的叠前反演结果。此时的叠前反演结果。

【技术实现步骤摘要】
二维批量叠前反演方法和白云岩储层预测方法及系统


[0001]本专利技术涉及石油地球物理勘探领域，特别涉及一种二维批量叠前反演方法和白云岩储层预测方法及系统。

技术介绍

[0002]在油田的勘探、开发生产和科研过程中，有时会面临只有大量二维地震资料，而缺乏三维地震资料的状况，在这种状况下需要充分利用这些二维地震资料进行储层和油藏参数的叠后或者叠前反演和预测。在这些参数当中纵横波阻抗等又是极其重要的参数，有时甚至是其它参数反演和预测的基础，如何得到高质量的纵横波反演结果是技术人员必须重视的问题。
[0003]在二维工区内，在过井的二维地震测线上建模较为简单，但是由于技术人员不可能保证工区内所有测线上的每一个层位都是闭和连通的，那么对于非过井测线的建模就成为了一个难题。距钻井较远的二维测线需制作伪井进行反演，这使得可靠度大幅下降，最终影响反演结果插值的准确性。另外，二维地震单条线反演工作量大，耗时长。
[0004]国内大部分盆地碳酸盐岩储层埋藏较深，大部分地区地震资料信噪比较低，碳酸盐岩储层非均质性较强，相比于碎屑岩在储层预测方面有天然的难度，基于此急需开发能够利用二维测线进行碳酸盐岩叠前储层有效预测的方法。
[0005]总之，现有方法无法实现有效地利用二维叠前资料进行深层碳酸盐岩储层预测工作。

技术实现思路

[0006]针对现有技术的不足，本专利技术的目的在于提供一种基于多次迭代二维批量叠前反演方法，该叠前反演方法得到的叠前反演结果能够更好的适用于碳酸盐岩储层中白云岩有效储层预测，提高深层白云岩有效储层预测精度。
[0007]为了实现上述目的，本专利技术提供了一种二维批量叠前反演方法，其中，该方法包括：
[0008]获取工区二维分偏移距叠加数据体(即分偏移距叠加二维地震数据体)；
[0009]对二维分偏移距叠加数据体进行融合处理，得到分偏移距叠加体融合地震数据体；
[0010]利用分偏移距叠加体融合地震数据体进行迭代二维批量叠前反演得到反演结果，具体包括：
[0011]基于分偏移距叠加体融合地震数据体及已钻井测井数据建立初始常速低频模型；
[0012]对初始常速低频模型进行叠前反演求解得到叠前反演结果；
[0013]判断反演结果与已钻井测井曲线(即过井测井曲线)之间的误差是否在预设误差范围内；
[0014]若误差不在预设误差范围内，基于叠前反演结果重新修改层位；以修改后的层位
为约束并利用叠前反演结果更新初始常速低频模型；对更新后的初始常速低频模型重新进行叠前反演求解更新叠前反演结果；直至反演结果与已钻井测井曲线之间的误差在预设误差范围内，迭代结束，输出此时的叠前反演结果作为利用分偏移距叠加体融合地震数据体进行迭代二维批量叠前反演的最终结果；
[0015]其中，反演结果包括纵波阻抗叠前反演结果、横波阻抗叠前反演结果和密度叠前反演结果。
[0016]在上述二维批量叠前反演方法中，优选地，叠前反演求解基于全Knott-zoeppritz方程进行。
[0017]在一具体实施方式中，利用全Knott-zoeppritz方程计算横波速度、纵波速度和密度；其中，全Knott-zoeppritz方程如下所示：
[0018][0019]其中，r
pp
为纵波入射时的纵波反射系数；r
ps
为纵波入射时的横波反射系数；r
sp
为横波入射时的纵波反射系数；r
ss
为横波入射时的横波反射系数；V
p1
为反射界面的上方的纵波速度；V
s1
为反射界面的上方的横波速度；ρ1为反射界面的上方的介质体密度；V
p2
为反射界面的下方的纵波速度；V
s2
为反射界面的下方的横波速度；ρ2为反射界面的下方的介质体密度；θ1为反射界面的上方的纵波角度；θ2为反射界面的下方的纵波角度；为反射界面的上方的横波角度；为反射界面的下方的横波角度；满足如下公式：
[0020]利用计算得到的横波速度、纵波速度和密度计算横波阻抗、纵波阻抗；其中，横波阻抗、纵波阻抗的计算公式如下所示：
[0021]I
P
＝V
P
×
ρ；I
S
＝V
S
×
ρ；
[0022]其中，I
P
为纵波阻抗；V
P
为纵波速度；I
S
为横波阻抗；V
S
为横波速度；ρ为介质密度。
[0023]在上述二维批量叠前反演方法中，优选地，初始常速低频模型包括纵波阻抗低频参数、横波阻抗低频参数和密度低频参数。
[0024]在上述二维批量叠前反演方法中，优选地，二维分偏移距叠加数据体包括近二维分偏移距叠加数据体、中二维分偏移距叠加数据体和远二维分偏移距叠加数据体。
[0025]在上述二维批量叠前反演方法中，优选地，该方法进一步包括：在二维分偏移距叠加数据体进行融合处理前先对二维分偏移距叠加数据进行闭合差校正。
[0026]在上述二维批量叠前反演方法中，优选地，对二维分偏移距叠加数据体进行融合处理为分别对近二维分偏移距叠加数据体、中二维分偏移距叠加数据体和远二维分偏移距叠加数据体进行融合处理。
[0027]在上述二维批量叠前反演方法中，优选地，融合处理包括：以追加融合方法按坐标逐线加载二维分偏移距叠加数据体；其中，逐线加载为逐线追加不覆盖。在一具体实施方式中，融合处理包括：每一条二维分偏移距叠加数据按先后顺序逐个加载到工区中，形成分偏移距叠加体融合地震数据体。
[0028]本专利技术还提供了一种白云岩储层预测方法，其中，该方法包括：
[0029]采用上述二维批量叠前反演方法得到工区纵波阻抗叠前反演结果和横波阻抗叠前反演结果；
[0030]基于纵波阻抗叠前反演结果和横波阻抗叠前反演结果确定工区横波阻抗分布和泊松比分布；
[0031]确定工区白云岩横波阻抗门槛值、含气白云岩泊松比门槛值；
[0032]基于白云岩横波阻抗门槛值以及确定的工区横波阻抗分布，采用门槛值法确定工区白云岩厚度分布；基于含气白云岩泊松比门槛值以及确定的工区泊松比分布，采用门槛值法确定工区含气白云岩厚度分布。
[0033]在上述白云岩储层预测方法中，优选地，确定工区白云岩横波阻抗门槛值、含气白云岩泊松比门槛值基于工区岩石物理量分析结果进行。
[0034]本专利技术还提供了一种二维批量叠前反演系统，其中，该系统包括：
[0035]数据获取模块：用于获取工区二维分偏移距叠加数据体(即分偏移距叠加二维地震数据体)；
[0036]融合处理模块：用于对二维分偏移距叠加数据体进行融合处理，得到分偏移距叠加体融合地震数据体；
[0037]模型建立模块：用于基于分偏移距叠加体融合地震数据体及已钻井测井数据建立初始常速低频模型；
[0038]反本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种二维批量叠前反演方法，其中，该方法包括：获取工区二维分偏移距叠加数据体；对所述二维分偏移距叠加数据体进行融合处理，得到分偏移距叠加体融合地震数据体；基于分偏移距叠加体融合地震数据体及已钻井测井数据建立初始常速低频模型；对初始常速低频模型进行叠前反演求解得到叠前反演结果；判断反演结果与已钻井测井曲线之间的误差是否在预设误差范围内；若误差不在预设误差范围内，基于叠前反演结果重新修改层位；以修改后的层位为约束并利用叠前反演结果更新初始常速低频模型；对更新后的初始常速低频模型重新进行叠前反演求解更新叠前反演结果；直至反演结果与已钻井测井曲线之间的误差在预设误差范围内，迭代结束，输出此时的叠前反演结果作为利用分偏移距叠加体融合地震数据体进行迭代二维批量叠前反演的最终结果；其中，所述反演结果包括纵波阻抗叠前反演结果、横波阻抗叠前反演结果和密度叠前反演结果。2.根据权利要求1所述的叠前反演方法，其中，所述叠前反演求解基于全Knott-zoeppritz方程进行。3.根据权利要求1所述的叠前反演方法，其中，所述初始常速低频模型包括纵波阻抗低频参数、横波阻抗低频参数和密度低频参数。4.根据权利要求1所述的叠前反演方法，其中，所述二维分偏移距叠加数据体包括近二维分偏移距叠加数据体、中二维分偏移距叠加数据体和远二维分偏移距叠加数据体。5.根据权利要求4所述的叠前反演方法，其中，该方法进一步包括：在二维分偏移距叠加数据体进行融合处理前先对二维分偏移距叠加数据进行闭合差校正。6.根据权利要求4或5所述的叠前反演方法，其中，对所述二维分偏移距叠加数据体进行融合处理为分别对近二维分偏移距叠加数据体、中二维分偏移距叠加数据体和远二维分偏移距叠加数据体进行融合处理。7.根据权利要求1-6任一项所述的叠前反演方法，其中，所述融合处理包括：以追加融合方法按坐标逐线加载二维分偏移距叠加数据体；其中，所述逐线加载为逐线追加不覆盖。8.一种白云岩储层预测方法，其中，该方法包括：采用权利要求1-7任一项所述的二维批量叠前反演方法得到工区纵波阻抗叠前反演结果和横波阻抗叠前反演结果；基于纵波阻抗叠前反演结果和横波阻抗叠前反演结果确定工区横波阻抗分布和泊松比分布；确定工区白云岩横波阻抗门槛值、含气白云岩泊松比门槛值；基于白云岩横波阻抗门槛值以及确定的工区横波阻抗分布，采用门槛值法确定工区白云岩厚度分布；基于含气白云岩泊松比门槛值以及确定的工区泊松比分布，采用门槛值法确定工区含气白云岩厚度分布。9.根据权利要求8所述的白云岩储层预测方法，其中，所述确定工区白云岩横波阻抗门槛值、含气白云岩泊松比门槛值基于工区岩石物理量分析结果进行。10.一种二维批量叠前反演系统，其中，该系统包括：
数据获取模块：用于获取工区二维分偏移距叠加数据体；融合处理模块：用于对所述二维分偏移距叠加数据体进行融合处理，得到分偏移距叠加体融合地震数据体；模型建立模块：用于基于分偏移距叠加体...

【专利技术属性】
技术研发人员：窦玉坛，刘峰，王永刚，杜广宏，张盟勃，张杰，朱军，边东辉，周齐刚，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：