利用地震属性控制初始波阻抗的建模方法及系统的制作方法
【技术领域】
[0001]本发明涉及地震油气勘探技术领域，尤其涉及一种利用地震属性控制初始波阻抗的建模方法及系统。
【背景技术】
[0002]传统的建模方法横向上插值权系数一般采用到井距离的倒数(反距离方式，weight；= 1./dis ；)，这是一种线性插值方式，一般难以模拟沉积相变化。结果往往围绕井画圆圈，地质体的边界反映不准确。
【发明内容】
[0003]本发明要解决的技术问题在于:针对现有技术中在横向上采用传统的插值方法难以模拟沉积相变化的缺陷，提供一种横向可以反映沉积相的变化，且可以提高地震反演的精度的利用地震属性控制初始波阻抗的建模方法及系统。
[0004]本发明解决其技术问题所采用的技术方案是:
[0005]提供一种利用地震属性控制初始波阻抗的建模方法，该方法包括以下步骤:
[0006]第一步:获取测井波阻抗值、坐标、地震解释层位以及断层数据；
[0007]第二步:读取建模地震道的坐标、层位以及地震属性信息；
[0008]第三步:在地震解释层位的控制下采用顶平行型、底平行型或者过渡型三种模式计算当前建模地震道时间t。与第i 口井对应时间t 1；
[0009]第四步:计算建模地震道时间t。与第i 口井的权系数weight 对所有井按插值公式(I)完成计算，得到该道该时刻初始模型值mod (t0):
[0010]mod(t0) = Swell (t;).weight；(I)
[0011]式中，mod(tO)表示建模道t。时刻的阻抗值，well (t ；)表示t。时刻对应第i 口井^时刻的阻抗值，weight i表示第i 口井的横向插值权系数，该权系数具体通过公式(2)得到:
[0012]Weighti=Corri(2)
[0013]其中，Corr1为建模位置地震属性与第i 口井位置地震属性的相似系数；
[0014]第五步:按照第三、四步完成当前建模地震道所有点的插值，用于地震反演。
[0015]本发明所述的方法中，第三步中在地震解释层位的控制下计算插值对应时间的基础上，利用地震相关分析进行修正，用δ i代替t 1；
[0016]其中，滑动时间δ i的确定，是通过在时间t Jii近滑动，计算每个滑动点的地震相关系数，取相关系数最大时的滑动时间δ i作为修正插值时间，再用t f δ 4寸刻代替t 1样点。
[0017]本发明所述的方法中，第四步中利用地震属性的相似系数代替传统的“反距离”插值权系数。
[0018]本发明所述的方法中，第四步中，初始模型值mod(tO)利用式(3)计算，
[0019]mod (tO) = Swell (t；+ δ.Corri0(3)
[0020]其中α为反映各向异性的参数，当地质现象横向变化较大时，通过增大参数a提高模型反映横向的变化的能力。
[0021]本发明产生的有益效果是:本发明在横向上，采用井旁地震记录与插值位置地震记录波形或者振幅等属性的相似系数作为横向插值权系数，替代传统的反距离作为井插值权系数，这样结果在横向上能够反映沉积相的变化；在纵向上，除了与传统方法一样，按照三种方式模拟沉积规律外，还根据地震资料的相关分析计算最佳的插值对应时间点。因此本发明改进后，建立的模型在一定程度上横向可以反映沉积相的变化，纵向减少解释层位的影响和准确刻画地质体，提高了地震反演的精度。
【附图说明】
[0022]下面将结合附图及实施例对本发明作进一步说明，附图中:
[0023]图1是本发明实施例利用地震属性控制初始波阻抗的建模方法的流程图；
[0024]图2是本发明实施例纵向上的三种模式定义插值方式；
[0025]图3是本发明实施例初始模型建立插值示意图；
[0026]图4是反映沉积相的分布图；
[0027]图5a是反距离建模的反演结果；
[0028]图5b是利用本发明的建模方法得到的反演结果。
【具体实施方式】
[0029]为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0030]地震反演需要建立一个初始的波阻抗模型，而且初始模型的质量关系到地震反演的精度。
[0031]如图1所示，本发明实施例利用地震属性控制初始波阻抗的建模方法主要包括以下步骤:
[0032]S1、获取测井波阻抗值、坐标、地震解释层位以及断层数据；
[0033]S2、读取当前地震道的坐标、层位以及地震属性信息；
[0034]S3、纵向上，在地震解释层位的控制下，计算当前建模地震道时间t。与第i 口井的对应时间ti;横向上，计算当前地震道时间t。、第i 口井对应时间tdP权系数weight ;，对所有井按插值公式(I)完成计算，得到该道该时刻初始模型值mod (tO):
[0035]mod (tO) = Swell (t;).weight；(I)
[0036]式中，mod (tO)表示当前道t。时刻的阻抗值，well (t ；)表示t。时刻对应第i 口井h时刻的阻抗值，Weighti表示第i 口井的横向插值权系数，该横向权系数具体通过公式(2)得到:
[0037]Weighti=Corri(2)；
[0038]其中，Corr1为当前插值位置地震道与各井井旁地震的地震属性相似系数；
[0039]S4、按照步骤S3完成当前地震道所有点的插值；
[0040]初始波阻抗模型建立方法如图3所示，其中Wl (tl)、W2 (t2)即为well (tx)，M(t0)即为mod (tO);建模过程中如何确定时间对应点和插值权系数是影响初始模型质量的两个最重要因素。本发明是利用井处的波阻抗值在层位如hl、h2的控制下加权插值得到。传统的插值方法权系数采用当