CN103149588B - 一种利用井震标定计算vti各向异性参数的方法及系统 - Google Patents

Abstract

本发明公开了一种利用井震标定计算VTI各向异性参数的方法及系统。其中，所述方法包括：沿地震剖面上存在各向异性的强反射界面拾取地震反射层位；根据地震反射层位，利用速度分析方法计算各层位的层速度，建立速度深度模型；根据速度深度模型计算得出各向异性地层厚度；根据速度深度模型，对地震数据进行叠前偏移处理，生成叠前偏移数据；利用各层位的层速度将叠前偏移数据和速度深度模型比例回时间域，生成时间域地震数据；在时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度；根据各向异性地层厚度以及各向同性地层厚度计算生成VTI各向异性参数。

Description

一种利用井震标定计算VTI各向异性参数的方法及系统

技术领域

本发明是关于地震勘探技术的方法，更具体地说，本发明是关于地震数据处理中一种通过井震标定手段来计算VTI介质各向异性参数δ的方法和系统。

背景技术

地层具有各向异性特性，地震资料处理时如果能进行准确的各向异性校正，不仅能改善地震资料的成像质量，而且处理结果能更准确地反映地层深度。VTI介质各向异性叠前深度偏移是现今地震资料处理的主流技术，所使用的偏移频散方程为： $\frac{{4\mathrm{\ω}}^2}{k_y^2+k_z^2}={\mathrm{\α}}_0^2[1+\frac{k_y^2({\mathrm{\δk}}_z^2+{\text{\ϵk}}_y^2)}{(k_z^2+k_y^2)}{]}^2,$ 各向异性校正的准确程度主要由参数δ和ε决定。由于地层是弱各向异性的，各向异性校正时δ的影响大于ε。计算参数ε时，要以先计算出的准确的δ为基础，因此，获得准确的δ是地震各向异性校正的关键。

现有的一种计算VTI各向异性参数的方式是根据测井分层数据近似计算δ的方法。这种方法存在两方面的局限：一是必须以已取得分层数据的地质界面作为各向异性分界面，实际上，分层数据所在的地质界面与各向异性分界面并不完全对应，导致各向异性速度模型与实际地质情况有较大差异；二是由于测井分层数据比较稀疏，这种方法只能计算出为数很少的几个大套地层的等效δ值，使精度受到限制。

发明内容

本发明所解决的技术问题克服了现有技术的缺点，提供了一种不是利用测井分层数据，而是利用井震标定计算VTI介质各向异性参数δ的方法。

本发明提供了一种利用井震标定计算VTI各向异性参数的方法，包括：沿地震剖面上存在各向异性的强反射界面拾取地震反射层位；根据所述地震反射层位，利用速度分析方法计算各层位的层速度，建立速度深度模型；根据所述速度深度模型计算得出各向异性地层厚度；根据所述速度深度模型，对地震数据进行叠前偏移处理，生成叠前偏移数据；利用所述各层位的层速度将所述叠前偏移数据和速度深度模型比例回时间域，生成时间域地震数据；在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度；根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数。

进一步的，所述利用速度分析方法计算各层位的层速度包括：利用Dix转换、相干反演来计算各层位的层速度。

进一步的，所述根据所述地震反射层位，利用速度分析方法计算各层位的层速度，建立速度深度模型还包括：利用剩余时差分析、反射旅行时层析成像技术，对所述速度深度模型进行优化迭代。

进一步的，所述获取的测井数据包括VSP走廊叠加数据或合成记录。

进一步的，所述制作合成记录包括收集包括声波测井曲线、密度测井曲线在内的测井曲线，进行井径校正、环境校正后制作合成地震记录。

进一步的，所述在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度，包括：将所述井旁道与所述测井数据进行波组关系、波形、相位的比对。

进一步的，所述根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数，包括根据关系式：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{{\mathrm{\ΔZ}}^I}{{\mathrm{\ΔZ}}^A}\right)}^2-1]$ 计算得出；

其中，ΔZ^A为所述各向异性地层厚度，ΔZ^I为所述各向同性地层厚度，δ为VTI各向异性参数。

本发明还提供了一种利用井震标定计算VTI各向异性参数的系统，包括：地震反射层位拾取单元，用于沿地震剖面上存在各向异性的强反射界面拾取地震反射层位；速度深度模型生成单元，用于根据所述地震反射层位，利用速度分析工具计算各层位的层速度，建立速度深度模型；各向异性地层厚度计算单元，用于根据所述速度深度模型计算得出各向异性地层厚度；叠前偏移数据生成单元，用于根据所述速度深度模型，对地震数据进行叠前深度偏移处理，生成叠前偏移数据；时间域地震数据生成单元，用于利用所述各层位的层速度将所述叠前偏移数据和速度深度模型比例回时间域，生成时间域地震数据；各向同性地层厚度计算单元，用于在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度；各向异性参数计算单元，用于根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数。

进一步的，所述速度深度模型生成单元中的速度分析工具包括Dix转换工具以及相干反演工具。

进一步的，所述速度深度模型生成单元还利用剩余时差分析、反射旅行时层析成像技术，对所述速度深度模型进行优化迭代。

进一步的，所述各向同性地层厚度计算单元获取的测井数据包括VSP走廊叠加数据或合成记录。

进一步的，所述合成记录包括收集包括声波测井曲线、密度测井曲线在内的测井曲线，进行井径校正、环境校正后制作合成地震记录。

进一步的，所述各向同性地层厚度计算单元，用于在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度，包括：将所述井旁道与所述测井数据进行波组关系、波形、相位的比对。

进一步的，所述各向异性参数计算单元，用于根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数，包括根据关系式：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{{\mathrm{\ΔZ}}^I}{{\mathrm{\ΔZ}}^A}\right)}^2-1]$ 计算得出；

其中，ΔZ^A为所述各向异性地层厚度，ΔZ^I为所述各向同性地层厚度，δ为VTI各向异性参数。

本发明与现有技术相比较，不局限于选择有分层数据的界面作为速度分界面，可根据实际需要，选择真正的各向异性分界面作为速度模型的分界面，使各向异性速度建模过程变得方便灵活，得到的速度模型与实际地质情况更符合；使用本方法计算的各向异性参数δ更均匀、更密集、更精细，可以更方便地利用叠前偏移手段，对δ进行迭代优化。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：

图1为本发明实施例利用井震标定计算VTI各向异性参数的方法流程图；

图2为本发明实施例利用井震标定计算VTI各向异性参数的系统示意图；

图3为本发明一利用井震标定计算VTI各向异性参数的具体实施例的流程图；

图4为本发明具体实施例的拾取的地震反射层位和速度深度模型示意图；

图5为本发明具体实施例的叠前深度偏移成像和比例回时间域与时间速度模型叠合示意图；

图6为本发明具体实施例的合成记录标定和反射层深度读取示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。

图1为本发明实施例利用井震标定计算VTI各向异性参数的方法流程图。如图所示，本发明实施例利用井震标定计算VTI各向异性参数的方法包括：

步骤S101，沿地震剖面上存在各向异性的强反射界面拾取地震反射层位；步骤S102，根据所述地震反射层位，利用速度分析方法计算各层位的层速度，建立速度深度模型；步骤S103，根据所述速度深度模型计算得出各向异性地层厚度；步骤S104，根据所述速度深度模型，对地震数据进行叠前偏移处理，生成叠前偏移数据；步骤S105，利用所述各层位的层速度将所述叠前偏移数据和速度深度模型比例回时间域，生成时间域地震数据；步骤S106，在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度；步骤S107，根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数。

在本实施例的步骤S101中，拾取的地震反射层位一般是重要的地质界面，可以作为理想的各向异性界面；拾取地震反射层位时尽量选择反射能量强、便于对比追踪的层位，同时，需要处理解释协同，以确保不遗漏重要的各向异性界面。在拾取的时间间隔上，每500ms-1000ms拾取一个层位比较适宜。

在本实施例中，步骤S102利用速度分析方法计算各层位的层速度是利用利用Dix转换、相干反演来计算。在建立速度深度模型时，还可利用剩余时差分析，反射旅行时层析成像技术，对速度深度模型进行优化迭代。

在本实施例中，根据叠前深度偏移原理，步骤S105中叠前偏移数据和速度深度模型应该有准确的对应的关系，如果速度深度模型中的速度界面和叠前偏移数据的偏移成像的反射同相轴之间不能很好的对应，则意味着速度深度模型存在误差，需要对速度深度模型进行优化迭代。

在本实施例中，步骤S106拾取测井数据包括通过输入VSP走廊叠加数据或制作合成记录得到。

其中，制作合成记录需要收集声波、密度等测井曲线，并进行井径校正、环境校正后得到合成记录。

步骤S106中将测井数据与井旁地震道数据进行比对包括：进行波组关系、波形、相位的比对，以确定步骤S101中拾取的地震反射层位的深度，也就是各岩层顶底界面埋藏深度。

在步骤S106中，井震标定是地震与测井联系的桥梁，是地震构造解释的基础，本领域技术人员可以据此计算出地震反射层位的准确深度，应用井震标定技术是本实施例的核心。

步骤S107中各向异性参数由所述各向异性地层厚度以及所述各向同性地层厚度依据关系式：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{{\mathrm{\ΔZ}}^I}{{\mathrm{\ΔZ}}^A}\right)}^2-1]$ 计算得出。

其中，ΔZ^A为所述各向异性地层厚度，ΔZ^I为所述各向同性地层厚度，δ为VTI各向异性参数。

图2为本发明实施例利用井震标定计算VTI各向异性参数的系统示意图。如图所示，本发明实施例利用井震标定计算VTI各向异性参数的系统包括：地震反射层位拾取单元201，用于沿地震剖面上存在各向异性的强反射界面拾取地震反射层位；速度深度模型生成单元202，用于根据所述地震反射层位，利用速度分析工具计算各层位的层速度，建立速度深度模型；各向异性地层厚度计算单元203，用于根据所述速度深度模型计算得出各向异性地层厚度；叠前偏移数据生成单元204，用于根据所述速度深度模型，对地震数据进行叠前深度偏移处理，生成叠前偏移数据；时间域地震数据生成单元205，用于利用所述各层位的层速度将所述叠前偏移数据和速度深度模型比例回时间域，生成时间域地震数据；各向同性地层厚度计算单元206，用于在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度；各向异性参数计算单元207，用于根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数。

在本实施例中，地震反射层位拾取单元201拾取的地震反射层位一般是重要的地质界面，可以作为理想的各向异性界面；拾取地震反射层位时尽量选择反射能量强、便于对比追踪的层位，同时，需要处理解释协同，以确保不遗漏重要的各向异性界面。在拾取的时间间隔上，每500ms-1000ms拾取一个层位比较适宜。

在本实施例中，速度深度模型生成单元202利用速度分析工具计算各层位的层速度包括利用Dix转换工具、相干反演工具。在建立速度深度模型时，还可利用剩余时差分析，反射旅行时层析成像技术，对速度深度模型进行优化迭代。

在本实施例中，根据叠前深度偏移原理，井旁地震道提取单元205中的叠前偏移数据和速度深度模型应该有准确的对应的关系，如果速度深度模型中的速度界面和叠前偏移数据的偏移成像的反射同相轴之间不能很好的对应，则意味着速度深度模型存在误差，需要对速度深度模型进行优化迭代。

在本实施例中，各向同性地层厚度计算单元206拾取的测井数据由VSP走廊叠加模块或合成记录模块得到。

其中，合成记录模块通过收集声波，密度等测井曲线，进行井径校正，环境校正后得到合成地震记。

在本实施例中，各向同性地层厚度计算单元206用于获取测井数据，并将测井数据与井旁地震道数据进行比对包括：进行波组关系、波形、相位的比对包括：进行波组关系、波形、相位的比对，以确定地震反射层位拾取单元201中拾取的地震反射层位的深度，也就是各岩层顶底界面埋藏深度。

各向异性参数计算单元207根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数，包括根据关系式：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{{\mathrm{\ΔZ}}^I}{{\mathrm{\ΔZ}}^A}\right)}^2-1]$ 计算得出。

其中，ΔZ^A为所述各向异性地层厚度，ΔZ^I为所述各向同性地层厚度，δ为VTI各向异性参数。

下面将结合图1及图2，以一具体实施例来说明本发明的实际应用。

图3至图6为本发明一具体实施例的示意图。其中，图3为本发明一利用井震标定计算VTI各向异性参数的具体实施例的流程图。

如图3所示，步骤S301，沿地震剖面上存在各向异性的强反射界面拾取X、Y、Z三个地震反射层位，作为各向异性分界面。图4所示为本实施例中拾取的地震反射层位和速度深度模型示意图，如图4中a部分所示，X、Y、Z为三个地震强反射层位，XY层位之间为火成岩地层，YZ之间为灰岩，两套岩层可能存在各向异性。

步骤S302，沿图4中a部分所拾取的层位，利用Dix转换、相干反演等速度分析工具，计算各层的层速度，建立速度深度模型(如图4中b部分所示)。其中图4中b部分的速度界面对应图4中a部分的各向异性界面，从速度深度模型中可以读出X、Y、Z三个层位在H7井处(图中黑点标记为H7)深度分别为4721米、4952米、5254米。由于计算层速度时，主要使用中远偏移距信息，因此可以底界面深度减去顶界面深度得到各向异性地层厚度ΔZ^A。在本实施例中，XY层位之间火成岩各向异性厚度为4952-4721=231米，YZ之间灰岩各向异性厚度为5254-4952=302米。

步骤S303，利用步骤S302中建立的速度深度模型，对地震数据进行叠前深度偏移处理，生成叠前偏移数据(如图5中a部分所示)，并利用所使用的层速度将叠前偏移数据和速度深度模型比例回时间域，生成时间域地震数据(如图5中b部分所示)。

步骤S304，从时间域地震数据中提取井旁地震道数据，将合成记录与井旁道井旁地震道数据对比（如图6中a部分所示），通过井震标定确定个个界面的深度。其中，合成记录的制作过程和获得的深度关系如图6中b部分所示，从图6中b部分的时间深度关系可以读出X、Y、Z三个反射层位的测量深度分别为5616米、5824米和6124米。由于测井解释中使用的深度为测量深度，而地震资料使用的是海拔深度，因此测井深度与地震深度之间整体相差一个井口海拔值。

步骤S305，对步骤S304中各界面标定的深度值，由底界面深度减去顶界面深度，获得地层的各向同性地层厚度ΔZ^I。因为VSP走廊叠加和合成记录都来自于测井资料，速度观测方向垂直于地层，不受地震远偏移距各向异性影像，因此计算的ΔZ^I是各向同性的地震厚度。在本实施例中，XY层位之间的火成岩厚度为5824-5616=208米，YZ层位之间的灰岩厚度为6124-5824=300米。

步骤S306，由步骤S302以及步骤S305计算VTI各向异性参数δ，根据关系式：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{{\mathrm{\ΔZ}}^I}{{\mathrm{\ΔZ}}^A}\right)}^2-1]$

计算可得出各地层的δ值。在本实施例中，XY之间的火成岩各向异性参数δ为：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{208}{231}\right)}^2-1]=-0.09107$

YZ之间的灰岩各向异性参数δ为：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{300}{302}\right)}^2-1]=-0.0066$

本实施例选择火成岩界面和灰岩界面作为速度分界面，使地质界面、速度界面和各向异性界面得到统一，使速度建模过程变得方便灵活，得到的速度模型与实际地质情况更符合，计算的各向异性参数δ更精细。

下面是本发明的另一具体实施例，对比现有技术和本发明的计算结果，可以了解本发明的应用效果。从浅到深有E₁、E₂、E₃、E₄、E₅、E₆、E₇共7套速度差异明显、具有VTI各向异性性质的地层，通过岩石物理测量方法，七套地层准确的δ值分别为0.1、0.1、0.05、-0.08、0.02、0.08、-0.05。

利用现有技术计算的数据：

已有E₂、E₅、E₆底界面和E₄、E₇中间含油层的分层数据，根据现有的分层数据法，由于E₄、E₇含油层的分层数据在地层中间而不在地层顶底界面上，不能用来计算各向异性参数δ值，能用于δ值计算的分层数据仅剩E₂、E₅、E₆底界面，只能计算出三套地层的等效δ值，分别是E₁E₂合并地层的等效δ值为0.091，E₃E₄E₅合并地层的等效δ值为0.012，E₆地层的δ值为0.078。

根据本发明，利用VSP走廊叠加资料计算的数据：

根据E₁到E₇的VSP走廊叠加资料，利用本发明的井震标定法计算出了E₁到E₇共7套地层的δ值，分别为0.110、0.101、0.049、-0.08、0.019、0.080、-0.049。

根据本发明，利用制作合成记录计算的数据：

根据E₁到E₇的声波和密度资料，制作合成记录，利用本发明的井震标定法计算出E₁到E₇共7套地层的δ值，分别为0.109、0.102、0.050、-0.081、0.021、0.078、-0.048。

利用现有技术计算的数据只得到了E₁E₂、E₃E₄E₅和E₆三套地层的等效δ值，与实际测量的δ值和真实地质情况均有很大差异，利用本发明计算的两组数据都得到了七套地层的δ值，与实测δ值比较，除E₁误差在10%外，E₂、E₃、E₄、E₅、E₆、E₇误差均小于5%，可见本发明计算的δ值更均匀、更密集、更精细。

本发明实施例不局限于选择有分层数据的界面作为速度分界面，还可根据实际需要，选择真正的各向异性分界面作为速度模型的分界面，使各向异性速度建模过程变得方便灵活，得到的速度模型与实际地质情况更符合；使用本方法计算的各向异性参数δ更均匀、更密集、更精细，可以更方便地利用叠前偏移手段，对δ进行迭代优化。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (14)

1.一种利用井震标定计算VTI各向异性参数的方法，其特征在于，包括：

沿地震剖面上存在各向异性的强反射界面拾取地震反射层位；

根据所述地震反射层位，利用速度分析方法计算各层位的层速度，建立速度深度模型；

根据所述速度深度模型计算得出各向异性地层厚度；

根据所述速度深度模型，对地震数据进行叠前偏移处理，生成叠前偏移数据；

利用所述各层位的层速度将所述叠前偏移数据和速度深度模型比例回时间域，生成时间域地震数据；

在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度；

根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数。

2.如权利要求1所述的方法，其特征在于，所述利用速度分析方法计算各层位的层速度包括：

利用Dix转换、相干反演来计算各层位的层速度。

3.如权利要求1所述的方法，其特征在于，所述根据所述地震反射层位，利用速度分析方法计算各层位的层速度，建立速度深度模型还包括：

利用剩余时差分析、反射旅行时层析成像技术，对所述速度深度模型进行优化迭代。

4.如权利要求1所述的方法，其特征在于，所述获取的测井数据包括VSP走廊叠加数据或合成记录。

5.如权利要求4所述的方法，其特征在于，所述合成记录包括收集包括声波测井曲线、密度测井曲线在内的测井曲线，进行井径校正、环境校正后制作的合成地震记录。

6.如权利要求1所述的方法，其特征在于，所述在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度，包括：将所述井旁道与所述测井数据进行波组关系、波形、相位的比对。

7.如权利要求1所述的方法，其特征在于，所述根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数，包括根据关系式：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{{\mathrm{\ΔZ}}^I}{{\mathrm{\ΔZ}}^A}\right)}^2-1]$ 计算得出；

其中，ΔZ^A为所述各向异性地层厚度；

ΔZ^I为所述各向同性地层厚度；

δ为VTI各向异性参数。

8.一种利用井震标定计算VTI各向异性参数的系统，其特征在于，包括：

地震反射层位拾取单元，用于沿地震剖面上存在各向异性的强反射界面拾取地震反射层位；

速度深度模型生成单元，用于根据所述地震反射层位，利用速度分析工具计算各层位的层速度，建立速度深度模型；

各向异性地层厚度计算单元，用于根据所述速度深度模型计算得出各向异性地层厚度；

叠前偏移数据生成单元，用于根据所述速度深度模型，对地震数据进行叠前深度偏移处理，生成叠前偏移数据；

时间域地震数据生成单元，用于利用所述各层位的层速度将所述叠前偏移数据和速度深度模型比例回时间域，生成时间域地震数据；

各向同性地层厚度计算单元，用于在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度；

各向异性参数计算单元，用于根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数。

9.如权利要求8所述的系统，其特征在于，所述速度深度模型生成单元中的速度分析工具包括Dix转换工具以及相干反演工具。

10.如权利要求8所述的系统，其特征在于，所述速度深度模型生成单元还利用剩余时差分析、反射旅行时层析成像技术，对所述速度深度模型进行优化迭代。

11.如权利要求8所述的系统，其特征在于，所述各向同性地层厚度计算单元获取的测井数据包括VSP走廊叠加数据或合成记录。

12.如权利要求11所述的系统，其特征在于，所述合成记录包括收集包括声波测井曲线、密度测井曲线在内的测井曲线，进行井径校正、环境校正后制作的合成地震记录。

13.如权利要求8所述的系统，其特征在于，所述各向同性地层厚度计算单元，用于在所述时间域地震数据中提取井旁道，与获取的测井数据进行比对，通过井震标定计算生成各向同性地层厚度，包括：将所述井旁道与所述测井数据进行波组关系、波形、相位的比对。

14.如权利要求8所述的系统，其特征在于，所述各向异性参数计算单元，用于根据所述各向异性地层厚度以及所述各向同性地层厚度计算生成VTI各向异性参数，包括根据关系式：

$\mathrm{\δ}=\frac{1}{2}[{\left(\frac{{\mathrm{\ΔZ}}^I}{{\mathrm{\ΔZ}}^A}\right)}^2-1]$ 计算得出；

其中，ΔZ^A为所述各向异性地层厚度；

ΔZ^I为所述各向同性地层厚度；

δ为VTI各向异性参数。