基于扩展弹性阻抗的地质储层检测方法
【技术领域】
[0001] 本发明总体来说涉及油气田勘探技术领域，属于地震资料解释范畴，更具体地讲， 涉及一种基于扩展弹性阻抗的地质储层检测方法。
【背景技术】
[0002] 地球物理勘探利用常规的纵波叠后偏移地震资料处理中损失的信息来解决构造 油气藏问题不大，但是要解决岩性复杂的油气藏就显得力不从心。随着叠前偏移技术的发 展，叠前道集中蕴藏的不同方位角、入射角地震信息可用来解决裂缝、岩性和流体识别技术 问题。
[0003] 弹性阻抗是声阻抗的概念延伸和推广，它建立在非零偏移距的基础上，是纵波速 度、横波速度、密度以及入射角的函数。弹性波阻抗反演属于叠前反演技术，它包含了丰富 的岩性以及流体信息，但是由于AVO效应的影响，常规弹性阻抗值随入射角发生剧烈变化 (即，弹性阻抗值随入射角的增大而急剧减小），可能会掩盖远、近角弹性阻抗值分析的某 些信息（例如，包含流体或岩性变化差异的信息）。
【发明内容】
[0004] 本发明的目的是提出一种基于扩展弹性阻抗的地质储层检测方法，以解决弹性阻 抗值随入射角的增大而急剧减小，无法有效识别地下地质储层的岩性和流体性质的技术问 题。
[0005] 本发明的一方面提供一种基于扩展弹性阻抗的地质储层检测方法，所述方法包 括：(a)获取目的层段的测井数据，在所述目的层段中选取多个采样点，并从所述测井数据 中提取与所述多个采样点中的每个采样点对应的纵波速度、横波速度和密度；(b)选取多 个实际入射角，并根据提取的纵波速度、横波速度和密度，来计算所述多个实际入射角中的 每个实际入射角在所述每个采样点处的反射系数；（c)选取多个理论入射角，并根据提取 的纵波速度、横波速度和密度，来计算所述多个理论入射角中的每个理论入射角在所述每 个采样点处的反射系数；（d)根据所述多个理论入射角的反射系数和所述多个实际入射角 的反射系数，分别确定与所述多个实际入射角中的每个实际入射角对应的等效理论入射 角；（e)利用基于所述等效理论入射角的扩展弹性阻抗来进行岩性和流体检测，以识别地 下地质储层的状态。
[0006] 可选地，在步骤（b)中，选取多个实际入射角的步骤可包括：获取实际入射角的初 始值，并将获取的实际入射角的初始值作为最大实际入射角，在零度至最大实际入射角的 范围之内以第一增量步长为间隔，选取所述多个实际入射角。
[0007] 可选地，在步骤（C)中，选取所述多个理论入射角的步骤可包括：在预定角度范围 内以第二增量步长为间隔，选取所述多个理论入射角。
[0008] 可选地，第二增量步长的数值可小于第一增量步长的数值。
[0009] 可选地，在步骤（C)中，根据提取的纵波速度、横波速度和密度，来计算所述多个 理论入射角中的每个理论入射角在所述每个采样点处的反射系数的步骤可包括：（Cl)根 据提取的纵波速度、横波速度和密度，计算所述多个理论入射角中的所述每个理论入射角 在所述每个采样点处的扩展弹性阻抗；(C2)根据所述每个理论入射角在所述每个采样点 处的扩展弹性阻抗，来计算所述每个理论入射角在所述每个采样点处的反射系数。
[0010] 可选地，在步骤（Cl)中，可利用下面的公式的来计算所述每个理论入射角在所述 每个采样点处的扩展弹性阻抗，
[0012] 其中，
[0013] Pj= cos X j+sin X j
[0014] qj= _8Ksin x j
[0015] rj= cos x j-4Ksin x s
[0016] 其中，EEIk ( x ,)为所述多个理论入射角中的第j个理论入射角x ,在所述多个采 样点中的第k个采样点处的扩展弹性阻抗，I < j <M，M为选取的理论入射角的总个数， I < k < N，N为选取的采样点的总个数，a k为与第k个采样点对应的纵波速度，β k为与第 k个采样点对应的横波速度，Pk为与第k个采样点对应的密度，α。为所述目的层段中的纵 波速度的初始值，为所述目的层段中的横波速度的初始值，P ^为所述目的层段中的密 度的初始值，K为纵横波速度比的平均值，K= [(0k+1/ak+1)2+(0k/a k)2]/2，ak+1为与第 k+1个采样点对应的纵波速度，β k+1为与第k+1个采样点对应的横波速度。
[0017] 可选地，在步骤（c2)中，可利用下面的公式的来计算所述每个理论入射角在所述 每个采样点处的反射系数，
[0019] 其中，Rk ( X ,)为所述多个理论入射角中的第j个理论入射角X ,在所述多个采样 点中的第k个采样点处的反射系数，I < j <M，M为选取的理论入射角的总个数，I < N， N为选取的采样点的总个数，EEIk+1( X ,)为所述多个理论入射角中的第j个理论入射角X , 在第k+Ι个采样点处的扩展弹性阻抗，EEIk ( X ,)为所述多个理论入射角中的第j个理论入 射角 '在第k个采样点处的扩展弹性阻抗。
[0020] 可选地，在步骤（b)中，根据提取的纵波速度、横波速度和密度，来计算所述多个 实际入射角中的每个实际入射角在所述每个采样点处的反射系数的步骤可包括：(bl)根 据提取的纵波速度、横波速度和密度，计算所述多个实际入射角中的每个实际入射角在所 述每个采样点处的弹性阻抗；(b2)根据所述每个实际入射角在所述每个采样点处的弹性 阻抗，来计算所述每个实际入射角在所述每个采样点处的反射系数。
[0021] 可选地，在步骤（d)中，根据所述多个理论入射角的反射系数和所述多个实际入 射角的反射系数，分别确定与所述多个实际入射角中的每个实际入射角对应的等效理论入 射角的步骤可包括：（dl)针对所述多个实际入射角中的任一实际入射角，分别计算所述任 一实际入射角的反射系数与所述多个理论入射角中的每个理论入射角的反射系数的相关 系数；(d2)确定所述相关系数中的最大值，并将所述相关系数中的最大值所对应的理论入 射角作为所述任一实际入射角对应的等效理论入射角。
[0022] 可选地，在步骤（dl)中，可利用下面的公式的来计算所述任一实际入射角的反射 系数与所述多个理论入射角中的每个理论入射角的反射系数的相关系数，
[0024] 其中，Corr[R( Θ丄R( X )]为所述多个实际入射角中的第i个实际入射角Θ ^勺 反射系数R(S1)与所述多个理论入射角中的第j个理论入射角X,的反射系数R(x ,)的 相关系数，I < i < η，η为选取的实际入射角的总个数，I < j < M，M为选取的理论入射角 的总个数，I < k < N，N为选取的采样点的总个数。
[0025] 采用本发明示例性实施例的上述基于扩展弹性阻抗的地质储层检测方法，可准确 地预测地下地质储层的岩性和流体性质，以真实反映地下地质储层的真实状态，有效降低 地质勘探风险，提尚钻探成功率。
【附图说明】
[0026] 通过下面结合附图进行的详细描述，本发明的上述和其它目的、特点和优点将会 变得更加清楚，其中：
[0027] 图1示出根据本发明示例性实施例的基于扩展弹性阻抗的地质储层检测方法的 流程图；
[0028] 图2示出根据本发明示例性实施例的计算每个实际入射角在每个采样点处的反 射系数的步骤的流程图；
[0029] 图3示出根据本发明示例性实施例的计算每个理论入射角在每个采样点处的反 射系数的步骤的流程图；
[0030] 图4示出根据本发明示例性实施例的确定与每个实际入射角对应的等效理论入 射角的步骤的流程图；
[0031] 图5示出根据本发明示例性实施例的弹性阻抗与扩展弹性阻抗分别随入射角变 化的关系图；
[0032] 图6示出根据本发明示例性实施例的实际入射角与理论入射角的等效换算图。
【具体实施方式】
[0033] 下面，将参照附图详细描述本发明的示例性实施例。
[0034] 提供参照附图的以下描述以帮助对由权利要求及其等同物限定的本发明的示例 性实施例的全面理解。包括各种特定细节以帮助理解，但这些细节仅被视为是示例性的。因 此，本领域的普通技术人员将认识到在不脱离本发明的范围和精神的情况下，可对描述于 此的示例性实施例进行各种改变和修改。此外，为了清楚和简洁，省略对公知的功能和结构 的描述。
[0035] 图1示出根据本发明示例性实施例的基于扩展弹性阻抗的地质储层检测方法的 流程图。
[0036] 在步骤Sl