本发明专利技术公开了一种用于进行全方位纵波速度分析的方法和装置。该方法包括：选择[δ(V),α]数据对，并得到相应的方位动校正速度进一步得到相应的短排列中各道的旅行时tS(x)；得到短排列中各道在与该旅行时tS(x)相关联的时窗内的总能量信息；选择使总能量最大的[δ(V),α]数据对；从而可确定特定自激自收时间样点下的方位动校正速度应用本公开可以对短排列的各道进行分析，以得到任意方位任意偏移距的方位动校正速度，提高对HTI介质进行速度分析的精度和效率，以便于更为直观和准确地了解地址构造和油气储藏情况等信息。应用本公开还可对全排列各道进行动校正。
【技术实现步骤摘要】
本公开涉及地球物理领域，更具体地，涉及一种用于进行全方位纵波速度分析的方法和一种用于进行全方位纵波速度分析的装置。
技术介绍
裂缝型储层是油气资源勘探的主要对象之一。据不完全统计，我国目前和在未来的几十年中，主要油气储量将来自裂缝型储层，约占总储量的60％以上。随着油气勘探深度的加大，这类储层在我国各油田均有分布。裂缝(例如垂直裂缝)会引起强的方位各向异性，地震波传播振幅、传播速度和衰减等属性都随着传播方位的变化而变化。通过研究介质的各向异性，可了解地质构造以及油气储藏情况等。可通过在激发点(也可称为炮点)激发地震波，在观测点(也可称为检波点)接收地震波并对其进行分析来研究介质的各向异性等。当观测点和激发点位于同一位置(也称为自激自收)时，观测点观测得的地震波信息能直观地反映界面的情况。实际中，观测点和激发点通常并不位于同一位置，因此可对观测值进行分析和处理，以消除炮检距(即激发点和观测点之间的距离)对反射波(例如反射波的速度和旅行时等)的影响，以便于进一步了解界面情况。纵波全方位地震数据的速度分析和动校正中主要需要考虑HTI的方位各向异性影响。目前对HTI方位各向异性介质进行速度分析和动校正的方法通常是将原始地震数据分解成多个方位数据体，并且在对每个方位数据体进行单独速度分析后，对多个方位的方位动校正速度场进行椭圆拟合，以建立方位各向异性速度场并进行全方位的旅行时动校正。这种方法不但计算量大，而且分方位后每个方位数据体的覆盖次数低，导致信号的信噪比低、速度谱分辨率低，使得容易受速度谱解释误差影响，不便于准确了解地址构造信息和油气储藏情况等。专利技术人意识到，提高HTI方位各向异性介质的全方位纵波速度分析的精度和效率是非常有必要的。
技术实现思路
本公开介绍了一种不需分方位即可对纵波进行全方位分析的方法和装置，以提高针对纵波的分析精度和效率。根据本公开的一方面，提出了一种用于进行全方位纵波速度分析的方法，该方法包括：从预定范围内选择至少两个[δ(V),α]数据对，并基于下式得到对应于每个[δ(V),α]数据对的方位动校正速度其中δ(V)为各向异性参数，α为裂缝方位角，为观测方位，为各向同性动校正速度；对应每个[δ(V),α]数据对，基于下式得到短排列的各道中每一道的旅行时tS(x)：其中tp0为自激自收时间样点，x为偏移距；对应每个[δ(V),α]数据对，得到短排列的各道在与旅行时tS(x)相关联的时窗内的总能量信息E1；从所述至少两个[δ(V),α]数据对中选择使总能量信息E1最大的[δd(V),αd]；基于下式确定自激自收时间样点为tp0时的方位动校正速度根据本公开的另一方面，提出了一种用于进行全方位纵波速度分析的装置，该装置包括：用于从预定范围内选择至少两个[δ(V),α]数据对、并基于下式得到对应于每个[δ(V),α]数据对的方位动校正速度的单元：其中δ(V)为各向异性参数，α为裂缝方位角，为观测方位，为各向同性动校正速度；对应每个[δ(V),α]数据对，用于基于下式得到短排列的各道中每一道的旅行时tS(x)的单元：其中tp0为自激自收时间样点，x为偏移距；对应每个[δ(V),α]数据对，用于得到短排列的各道在与旅行时tS(x)相关联的时窗内的总能量信息E1的单元；用于从所述至少两个[δ(V),α]数据对中选择使总能量信息E1最大的[δd(V),αd]的单元；用于基于下式确定自激自收时间样点为tp0时的方位动校正速度的单元：本公开的各方面可对短排列的各道进行分析以得到任意方位任意偏移距的方位动校正速度，提高对HTI介质进行速度分析的精度和效率，以便于更为直观和准确地了解地址构造和油气储藏情况等信息。进一步地，本公开还可得到任意方位和任意偏移距的旅行时动校正量，以提高针对HTI介质的动校正精度和效率。附图说明通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。图1示出了根据本公开的一个实施方式的用于进行全方位纵波速度分析的方法的流程图。图2示出了根据本公开的一个实施方式的用于进行全方位纵波速度分析和动校正的方法的流程图。图3示出了从共反射点道集中抽取的共偏移距道集。图4示出了根据本公开的一个实施例得到的[δ(V),α]谱的示意图。图5示出了根据本公开的一个实施例得到的椭圆动校正速度的示意图。图6示出了根据本公开的一个实施例对共反射点道集进行动校正后得到的共偏移距道集。具体实施方式下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。以下介绍了本公开所基于的原理。应当理解的是，本公开中涉及的术语“长排列”和“短排列”是本领域技术人员根据需要和现场情况具体设定的。针对单层介质，VTI介质的纵波长排列旅行时tV(x)可表示为(AlkhalifahandTsvankin，1995)：tV2(x)=tp02+x2Vnmo2-2ηx4Vnmo2[Vnmo2tp02+(1+2η)x2],---(1)]]>其中，tp0是自激自收双程旅行时(也可称为自激自收时间样点)，x是偏移距，η是纵波各向异性参数，Vnmo是纵波动校正速度。根据等效VTI介质理论，HTI介质对称轴面的反射波旅行时方程为：tH2(x)=tp02+x2Vnmo2(0)-2η(V)x4Vnmo2(0)[Vnmo2(0)tp02+(1+2η(V))x2],---(2)]]>其中η(V)是HTI介质的等效纵波的各向异性参数，Vnmo(0)是HTI对称轴面(垂直裂缝面)的纵波动校正速度。在HTI介质的各向同性面内(平行裂缝面)，地震波传播速度为各向同性，其动校正速度为Vnmo(90)。在HTI介质中，方位动校正速度呈椭圆特征。当裂缝的对称轴方位角(本公开中可简称为裂缝方位角)为α时，对应于任意观测方位的动校正速度为：设纵波垂直速度为Vp0，则在各向同性面内有在对称轴面内有δ(V)为各项异性参数。因此式(3)可简化为：在弱各向异性条件下：根据现有技术可得到：针对HTI介质，任意观测方位的旅行时方程(包括长排列和短排列)为：其中，任意观测方位的短排列旅行时方程可简化表示为：上述式(1)-(7)是本领域已知的公式。专利技术人经过深入研究后得到其中表示各向同性动校正速度，详细推导步骤可参见下列描述。当输入全方位短排列时，可用各向同性动校正速度对全方位数据进行速度分析时，即使得用计算得到的旅行时与实际全方位旅行时误差最小，即：式(8)中i可表示偏移距序号，j可表示观测方位序号，xj,i表示偏移距序号i、观测方位序号j所对应的观测点的偏移距，表示观测方位序号j所对应的观测方位。可得到(例如采用最小二乘原理)：式(9)可写成积分形式，以得到积分公式：把式(4)代入式(10)可得的表达式：V~nmo=Vp01+2δ(V本文档来自技高网...
【技术保护点】
一种用于进行全方位纵波速度分析的方法，该方法包括：从预定范围内选择至少两个[δ(V),α]数据对，并基于下式得到对应于每个[δ(V),α]数据对的方位动校正速度其中δ(V)为各向异性参数，α为裂缝方位角，为观测方位，为各向同性动校正速度；对应每个[δ(V),α]数据对，基于下式得到短排列的各道中每一道的旅行时tS(x)：其中tp0为自激自收时间样点，x为偏移距；对应每个[δ(V),α]数据对，得到短排列的各道在与旅行时tS(x)相关联的时窗内的总能量信息E1；从所述至少两个[δ(V),α]数据对中选择使总能量信息E1最大的[δd(V),αd]；基于下式确定自激自收时间样点为tp0时的方位动校正速度
【技术特征摘要】
1.一种用于进行全方位纵波速度分析的方法，该方法包括：从预定范围内选择至少两个[δ(V),α]数据对，并基于下式得到对应于每个[δ(V),α]数据对的方位动校正速度其中δ(V)为各向异性参数，α为裂缝方位角，为观测方位，为各向同性动校正速度；对应每个[δ(V),α]数据对，基于下式得到短排列的各道中每一道的旅行时tS(x)：其中tp0为自激自收时间样点，x为偏移距；对应每个[δ(V),α]数据对，得到短排列的各道在与旅行时tS(x)相关联的时窗内的总能量信息E1；从所述至少两个[δ(V),α]数据对中选择使总能量信息E1最大的[δd(V),αd]；基于下式确定自激自收时间样点为tp0时的方位动校正速度2.根据权利要求1所述的方法，该方法还包括通过下面的步骤得到基于得到：其中i表示偏移距序号，j表示观测方位序号，xj,i表示偏
\t移距序号i、观测方位序号j所对应的观测点的偏移距，表示观测方位序号j所对应的观测方位；进一步得到积分公式：将代入所述积分公式，得到的表达式：V~nmo=Vp01+2δ(V)1+23δ(V);]]>将所述的表达式代回得到：3.根据权利要求1所述的方法，该方法还包括：从预定范围内选择至少两个η(V)，其中η(V)为各向异性参数；对应每个η(V)，基于下式得到长排列的各道中每一道的旅行时tL(x)：对应每个η(V)，得到长排列的各道在与旅行时tL(x)相关联的时窗内的总能量信息E2；从所述至少两个η(V)中选择使总能量信息E2最大的ηd(V)；基于下式确定自激自收时间样点为tp0时全排列的各道中每一道的动校正量Δt，所述全排列包括短排列和长排列：4.根据权利要求1所述的方法，其中，各向同性动校正速度是基于短排列的各道确定的。5.根据权利要求1所述的方法，其中，与旅行时ts(x)相关联的时窗是以旅行时tS(x)为中心的时窗。6.根据权利要求3所述的方法，其中，与旅行时tL(x...
【专利技术属性】
技术研发人员：陈天胜，刘春园，刘韬，刘炯，张雅勤，魏修成，季玉新，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。