CN111505712A - 地震层速度计算方法及装置 - Google Patents

地震层速度计算方法及装置 Download PDF

Info

Publication number
CN111505712A
CN111505712A CN201910096375.5A CN201910096375A CN111505712A CN 111505712 A CN111505712 A CN 111505712A CN 201910096375 A CN201910096375 A CN 201910096375A CN 111505712 A CN111505712 A CN 111505712A
Authority
CN
China
Prior art keywords
velocity
layer
stacking
spectrum
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910096375.5A
Other languages
English (en)
Inventor
杜磊
郭平
刘宝鸿
柳世光
冉波
宫振超
李理
殷敬红
韩敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201910096375.5A priority Critical patent/CN111505712A/zh
Publication of CN111505712A publication Critical patent/CN111505712A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362—Effecting static or dynamic corrections; Stacking
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/52—Move-out correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种地震层速度计算方法及装置，该方法包括：根据地震处理工区的叠加速度谱，建立层速度模型；对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。本发明可以计算地震层速度，准确性高。

Description

地震层速度计算方法及装置
技术领域
本发明涉及地震勘探领域，尤其涉及一种地震层速度计算方法及装置。
背景技术
目前，地震勘探中常采用DIX公式(迪克斯公式)计算层速度，DIX公式是一种用均方根速度来求解层速度的公式，因为层速度可以反映和求得某一地层很多信息，该方法计算简单、快捷，但是，采用DIX公式来计算层速度，是基于地下地层为水平层状均匀介质的，该假设过于理想化，无法适应地下地层速度反转的真实情况，因为地下地层确实存在高速层在上，低速层在下的情况(如火山岩覆盖区速度比下伏地层速度高)，也存在高速层在低速层的侧面，速度横向变化剧烈(如盐体刺穿区速度比围岩地层速度高)的情况，因此DIX公式计算的层速度，准确性不高。
发明内容
本发明实施例提出一种地震层速度计算方法，用以计算地震层速度，准确性高，该方法包括：
根据地震处理工区的叠加速度谱，建立层速度模型；
对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；
若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。
本发明实施例提出一种地震层速度计算装置，用以计算地震层速度，准确性高，该装置包括：
层速度模型建立模块，用于根据地震处理工区的叠加速度谱，建立层速度模型；
叠加速度谱获得模块，用于对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；
层速度模型校正模块，用于若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。
本发明实施例还提出了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述地震层速度计算方法。
本发明实施例还提出了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述地震层速度计算方法的计算机程序。
在本发明实施例中，根据地震处理工区的叠加速度谱，建立层速度模型；对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。在本发明实施例中，若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型，因此，获得的层速度模型的准确性高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中地震层速度计算方法的流程图；
图2为本发明实施例中层速度模型的示意图；
图3为采用DIX公式计算得到的层速度的目标测线速度剖面图；
图4为本发明实施例中采用本发明的方法计算得到的层速度的目标测线速度剖面图；
图5为采用DIX公式计算得到的层速度剖面图；
图6为本发明实施例中采用本发明的方法计算得到的层速度剖面图；
图7为本发明实施例的地震层速度计算装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1为本发明实施例中地震层速度计算方法的流程图，如图1所示，该方法包括：
步骤101，根据地震处理工区的叠加速度谱，建立层速度模型；
步骤102，对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；
步骤103，若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。
在本发明实施例中，若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型，因此，获得的层速度模型的准确性高。
具体实施时，首先需要获取地震处理工区的叠加速度谱，要求该叠加速度谱分布均匀且变化规律，无明显异常值，每一个叠加速度谱点对应一个叠前CDP道集(共反射点)道集，叠加速度谱可以用数据集D＝{(t0,i，Vsi，αi)，i＝1,2,…,N}表示，其中，t0,i为叠加速度谱中第i层的自激自收时刻，Vsi为叠加速度谱中第i层的叠加速度，αi为t0,i对应的反射波同相轴的时间梯度，N为地震层的数目。
在一实施例中，对层速度模型进行正演模拟，获得正演模拟的叠加速度谱，可以包括：
根据层速度模型，建立反射波时距曲线；
根据反射波时距曲线，获得叠加速度谱。
具体实施时，可以根据层速度模型，采用射线追踪法求出反射波时距曲线，时距曲线指地震波走时与距离的关系曲线，即地震波到达各检波点的时间同检波点到***点的距离之间的关系曲线，射线追踪法是一种给定发射点和接收点位置及介质的波速，求解两点间射线路径及电磁波传播时间的方法，主要分为正向算法和反向算法两种。正向算法就是将发射出的电磁波等分为均匀的射线束，对每条射线进行追踪，射线经过数次的反射绕射后被接收点处接受，未到达或未在接收点范围内的射线舍弃。反向算法就是根据几何光学原理，从接收点开始反过来搜寻每一条能从接收点到达发射点的路径。
根据反射波时距曲线，获得叠加速度谱的具体过程可以包括：
对反射波时距曲线进行动校正，获得叠加速度谱，动校正用的速度合适，则叠加后有效波的能量最强，这时的速度称为最佳叠加速度。反之，动校正的速度用的不合适，叠加后有效波能量就要削弱。叠加记录的振幅随叠加速度而变化,这就是叠加速度谱。
根据反射波时距曲线，还可以采用最小平方法获得叠加速度谱，具体过程可以包括：
将反射波时距曲线与如下双曲线进行拟合：
Figure BDA0001964632960000041
由最小平方准则，得到如下目标函数：
Figure BDA0001964632960000042
令
Figure BDA0001964632960000043
达到极小，则：
Figure BDA0001964632960000044
Figure BDA0001964632960000045
根据以上两个公式可得：
Figure BDA0001964632960000046
Figure BDA0001964632960000051
对于第i个地震层的反射界面，可根据如下公式计算对应的叠加速度：
Figure BDA0001964632960000052
其中，
Figure BDA0001964632960000053
为叠加速度谱中第i层的叠加速度；
Figure BDA0001964632960000054
为层速度模型中第i层的层速度。
将层速度模型中的αi用层速度模型中的倾角替换，即可得到叠加速度谱的新的数据集D＝{(t0,i，Vsi，αi)，i＝1,2,…,N}。
在一实施例中，层速度模型可以包括多个地震层的层速度、倾角和共中心点处深度，图2为本发明实施例中层速度模型的示意图，如图2所示，在一个CDP道集的局部范围内，可以将地下介质抽象为倾斜层状介质，该模型由参数集M＝{hi，Vi，βi}，i＝1,2,…,N}唯一确定，其中hi为共中心点处界面的深度，Vi为第i层的层速度，βi为倾角，N为地震层的数目。
在一实施例中，根据地震处理工区的叠加速度谱，建立层速度模型，可以包括：
根据地震处理工区的叠加速度谱，采用倾角校正的DIX公式，计算层速度模型中每一层的层速度：
Figure BDA0001964632960000055
θi-1＝arcsin(0.5αi-1Vsi-1)
θi＝arcsin(0.5αiVsi)
其中，Vi为层速度模型中第i层的层速度；
t0,i和t0,i-1分别为叠加速度谱中第i层和第i-1层的自激自收时刻；
Vsi和Vsi-1分别为叠加速度谱中第i层的叠加速度和第i-1层的叠加速度；
αi为t0,i对应的反射波同相轴的时间梯度，αi-1为t0,i-1对应的反射波同相轴的时间梯度。
在现有技术中，常采用如下的DIX公式计算层速度模型中每一层的层速度：
Figure BDA0001964632960000056
其中，Vi为层速度模型中第i层的层速度；
t0,i和t0,i-1分别为叠加速度谱中第i层和第i-1层的自激自收时刻；
Vsi和Vsi-1分别为叠加速度谱中第i层的叠加速度和第i-1层的叠加速度。
该DIX公式需满足三个假设条件：一是水平层状介质，二是射线垂直入射，三是速度横向无变化，该公式没有考虑炮检距和倾角，计算精度不高，而采用倾角校正的DIX公式计算时，地震层无需是是水平层状介质，计算精度高于DIX公式。
在一实施例中，根据地震处理工区的叠加速度谱，建立层速度模型，包括：
根据地震处理工区的叠加速度谱，采用如下公式，计算层速度模型中每一层的倾角：
Figure BDA0001964632960000061
Figure BDA0001964632960000062
Figure BDA0001964632960000063
其中，βi为层速度模型中第i层的倾角；
V1、Vj-1和Vj分别为层速度模型中第1层、第j-1层和第j层的层速度；
αi为叠加速度谱中第i层的自激自收时刻t0,i对应的反射波同相轴的时间梯度。
在一实施例中，根据地震处理工区的叠加速度谱，建立层速度模型，包括：
根据地震处理工区的叠加速度谱，采用如下公式，计算层速度模型中每一层的共中心点处深度：
Figure BDA0001964632960000064
Figure BDA0001964632960000065
其中，hi为层速度模型中第i层的共中心点处深度；
βK和βj分别为层速度模型中第K层和第j层的倾角；
Vj为层速度模型中第j层的层速度；
t0,K为叠加速度谱中第K层的自激自收时刻。
下面给出一具体实施例，说明本发明提出的地震层速度计算方法的具体应用，首先获取地震处理工区的叠加速度谱，该叠加速度谱分布均匀且变化规律，无明显异常值，叠加速度谱用数据集D＝{(t0,i，Vsi，αi)，i＝1,2,…,N}表示，其中，t0,i为叠加速度谱中第i层的自激自收时刻，Vsi为叠加速度谱中第i层的叠加速度，αi为t0,i对应的反射波同相轴的时间梯度，N为地震层的数目。
根据地震处理工区的叠加速度谱，采用倾角校正的DIX公式，计算层速度模型中每一层的层速度：
Figure BDA0001964632960000071
θi-1＝arcsin(0.5αi-1Vsi-1)
θi＝arcsin(0.5αiVsi)
其中，Vi为层速度模型中第i层的层速度；
t0,i和t0,i-1分别为叠加速度谱中第i层和第i-1层的自激自收时刻；
Vsi和Vsi-1分别为叠加速度谱中第i层的叠加速度和第i-1层的叠加速度；
αi为t0,i对应的反射波同相轴的时间梯度，αi-1为t0,i-1对应的反射波同相轴的时间梯度。
根据地震处理工区的叠加速度谱，采用如下公式，计算层速度模型中每一层的倾角：
Figure BDA0001964632960000072
Figure BDA0001964632960000073
Figure BDA0001964632960000074
其中，βi为层速度模型中第i层的倾角；
V1、Vj-1和Vj分别为层速度模型中第1层、第j-1层和第j层的层速度；
αi为叠加速度谱中第i层的自激自收时刻t0,i对应的反射波同相轴的时间梯度。
根据地震处理工区的叠加速度谱，采用如下公式，计算层速度模型中每一层的共中心点处深度：
Figure BDA0001964632960000075
Figure BDA0001964632960000076
其中，hi为层速度模型中第i层的共中心点处深度；
βK和βj分别为层速度模型中第K层和第j层的倾角；
Vj为层速度模型中第j层的层速度；
t0,K为叠加速度谱中第K层的自激自收时刻。
然后根据层速度模型，采用射线追踪法求出反射波时距曲线，对反射波时距曲线进行动校正，获得叠加速度谱。
将叠加速度谱中的αi用层速度模型中的倾角替换，即可得到叠加速度谱的新的数据集D1＝{(t0,i，Vsi *，αi *)，i＝1,2,…,N}。
若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，即
Figure BDA0001964632960000081
其中ε为预设范围，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。
图3为采用DIX公式计算得到的层速度的目标测线速度剖面图，图4为本发明实施例中采用本发明的方法计算得到的层速度的目标测线速度剖面图，参考图3和图4可知，采用本发明的方法计算得到的层速度渐变趋势更加合理，更接近地下真实构造形态，有益于地震处理的叠前偏移速度建模。
图5为采用DIX公式计算得到的层速度剖面图，图6为本发明实施例中采用本发明的方法计算得到的层速度剖面图，参考图5和图6的椭圆形标出的部分可知，DIX公式计算得到的层速度在速度剖面上渐变不柔和，说明层速度突变现象较严重，不利用偏移成像，采用本发明的方法计算得到的层速度过渡自然，有利用后续偏移成像。
在本发明实施例中，根据地震处理工区的叠加速度谱，建立层速度模型；对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型，获得的层速度模型的准确性高。
另外，采用本发明的方法计算不受水平层状介质的限制，得到的层速度渐变趋势更加合理，更接近地下真实构造形态，有益于地震处理的叠前偏移速度建模；且层速度过渡自然，有利用后续偏移成像。
基于同一发明构思，本发明实施例还提供了一种地震层速度计算装置，如下面的实施所述。由于这些解决问题的原理与地震层速度计算方法相似，因此装置的实施可以参见方法的实施，重复之处不在赘述。
图7为本发明实施例的地震层速度计算装置的结构示意图，如图7所示，该装置包括：
层速度模型建立模块701，用于根据地震处理工区的叠加速度谱，建立层速度模型；
叠加速度谱获得模块702，用于对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；
层速度模型校正模块703，用于若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。
在一实施例中，叠加速度谱获得模块702可以具体用于：
根据层速度模型，建立反射波时距曲线；
根据反射波时距曲线，获得叠加速度谱。
综上所述，在本发明实施例中，根据地震处理工区的叠加速度谱，建立层速度模型；对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型，获得的层速度模型的准确性高。
另外，采用本发明的方法计算不受水平层状介质的限制，得到的层速度渐变趋势更加合理，更接近地下真实构造形态，有益于地震处理的叠前偏移速度建模；且层速度过渡自然，有利用后续偏移成像。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种地震层速度计算方法，其特征在于，包括：
根据地震处理工区的叠加速度谱，建立层速度模型；
对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；
若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。
2.如权利要求1所述的地震层速度计算方法，其特征在于，对层速度模型进行正演模拟，获得正演模拟的叠加速度谱，包括：
根据层速度模型，建立反射波时距曲线；
根据反射波时距曲线，获得叠加速度谱。
3.如权利要求1所述的地震层速度计算方法，其特征在于，层速度模型包括多个地震层的层速度、倾角和共中心点处深度。
4.如权利要求3所述的地震层速度计算方法，其特征在于，根据地震处理工区的叠加速度谱，建立层速度模型，包括：
根据地震处理工区的叠加速度谱，采用倾角校正的DIX公式，计算层速度模型中每一层的层速度：
Figure FDA0001964632950000011
θi-1＝arcsin(0.5αi-1Vsi-1)
θi＝arcsin(0.5αiVsi)
其中，Vi为层速度模型中第i层的层速度；
t0,i和t0,i-1分别为叠加速度谱中第i层和第i-1层的自激自收时刻；
Vsi和Vsi-1分别为叠加速度谱中第i层的叠加速度和第i-1层的叠加速度；
αi为t0,i对应的反射波同相轴的时间梯度，αi-1为t0,i-1对应的反射波同相轴的时间梯度。
5.如权利要求3所述的地震层速度计算方法，其特征在于，根据地震处理工区的叠加速度谱，建立层速度模型，包括：
根据地震处理工区的叠加速度谱，采用如下公式，计算层速度模型中每一层的倾角：
Figure FDA0001964632950000021
Figure FDA0001964632950000022
Figure FDA0001964632950000023
其中，βi为层速度模型中第i层的倾角；
V1、Vj-1和Vj分别为层速度模型中第1层、第j-1层和第j层的层速度；
αi为叠加速度谱中第i层的自激自收时刻t0,i对应的反射波同相轴的时间梯度。
6.如权利要求3所述的地震层速度计算方法，其特征在于，根据地震处理工区的叠加速度谱，建立层速度模型，包括：
根据地震处理工区的叠加速度谱，采用如下公式，计算层速度模型中每一层的共中心点处深度：
Figure FDA0001964632950000024
Figure FDA0001964632950000025
其中，hi为层速度模型中第i层的共中心点处深度；
βK和βj分别为层速度模型中第K层和第j层的倾角；
Vj为层速度模型中第j层的层速度；
t0,K为叠加速度谱中第K层的自激自收时刻。
7.一种地震层速度计算装置，其特征在于，包括：
层速度模型建立模块，用于根据地震处理工区的叠加速度谱，建立层速度模型；
叠加速度谱获得模块，用于对层速度模型进行正演模拟，获得正演模拟的叠加速度谱；
层速度模型校正模块，用于若叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值不在预设范围内，则将地震处理工区的叠加速度谱更新为正演模拟的叠加速度谱，重新执行以上步骤，直至叠加速度谱中的叠加速度与正演模拟的叠加速度谱中的叠加速度的差值在预设范围内，输出层速度模型。
8.如权利要求7所述的地震层速度计算装置，其特征在于，叠加速度谱获得模块进一步用于：
根据层速度模型，建立反射波时距曲线；
根据反射波时距曲线，获得叠加速度谱。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至6任一所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至6任一所述方法的计算机程序。
CN201910096375.5A 2019-01-31 2019-01-31 地震层速度计算方法及装置 Pending CN111505712A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910096375.5A CN111505712A (zh) 2019-01-31 2019-01-31 地震层速度计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910096375.5A CN111505712A (zh) 2019-01-31 2019-01-31 地震层速度计算方法及装置

Publications (1)

Publication Number Publication Date
CN111505712A true CN111505712A (zh) 2020-08-07

Family

ID=71877367

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910096375.5A Pending CN111505712A (zh) 2019-01-31 2019-01-31 地震层速度计算方法及装置

Country Status (1)

Country Link
CN (1) CN111505712A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114114392A (zh) * 2020-09-01 2022-03-01 中国石油天然气股份有限公司 层速度模型建立方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4984220A (en) * 1989-03-06 1991-01-08 Amoco Corporation Geophysical exploration using velocity spectra regional coherency peaks
CN104122582A (zh) * 2014-06-24 2014-10-29 夏正元 利用叠加速度求取高精度地震波速度的方法
CN105093281A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 一种反演框架下的地震多波建模方法
CN106405651A (zh) * 2016-11-14 2017-02-15 中国石油化工股份有限公司 一种基于测井匹配的全波形反演初始模型构建方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4984220A (en) * 1989-03-06 1991-01-08 Amoco Corporation Geophysical exploration using velocity spectra regional coherency peaks
CN105093281A (zh) * 2014-05-16 2015-11-25 中国石油化工股份有限公司 一种反演框架下的地震多波建模方法
CN104122582A (zh) * 2014-06-24 2014-10-29 夏正元 利用叠加速度求取高精度地震波速度的方法
CN106405651A (zh) * 2016-11-14 2017-02-15 中国石油化工股份有限公司 一种基于测井匹配的全波形反演初始模型构建方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
PING GUO ET AL.: "Wide-azimuth Acquisition and Azimuthal Anisotropy Processing for Complex Structural Imaging in Liaohe Oilfield", 《PROCEEDINGS OF THE 11TH SEGJ INTERNATIONAL SYMPOSIUM, 2013》 *
WEIJIAN MAO ET AL.: "Automated Interval Velocity Inversion", 《SEG 2000 EXPANDED ABSTRACTS》 *
杨金华等: "复杂构造区 t0 时间漂移和叠加速度特征分析及其校正", 《石油地球物理勘探》 *
王宏语: "《地震资料地质解释与应用》", 30 April 2016, 北京：地质出版社 *
甘永忠: "广角层速度分析", 《石油物探》 *
甘永忠等: "由叠加速度谱反演层速度", 《1992年中国地球物理学会第八届学术年会》 *
钱绍瑚: "《地震勘探》", 30 April 1989, 中国地质大学出版社 *
陈传仁等: "一种精确层速度反演方法", 《物探化探计算技术》 *
黄延祜等: "水平层状介质中多次叠加的某些理论问题", 《地球物理学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114114392A (zh) * 2020-09-01 2022-03-01 中国石油天然气股份有限公司 层速度模型建立方法及装置
CN114114392B (zh) * 2020-09-01 2023-11-28 中国石油天然气股份有限公司 层速度模型建立方法及装置

Similar Documents

Publication Publication Date Title
EP3259619B1 (en) Method and system of processing seismic data by providing surface aperture common image gathers
EP2715405B1 (en) Method of processing seismic data by providing surface offset common image gathers
CN105093292A (zh) 一种地震成像的数据处理方法和装置
CN105277978A (zh) 一种确定近地表速度模型的方法及装置
CN105093319B (zh) 基于三维地震数据的地面微地震静校正方法
CN103149586A (zh) 一种倾斜层状粘弹性介质中波场正演模拟方法
CN112698390B (zh) 叠前地震反演方法及装置
CN101614826A (zh) 三维地震数据处理中实现面元均化的方法和装置
CN105629299A (zh) 角度域叠前深度偏移的走时、角度表获取方法及成像方法
CN106199704B (zh) 一种三维三分量海底电缆地震资料速度建模方法
CN108303736B (zh) 各向异性ti介质最短路径射线追踪正演方法
CN106338766B (zh) 基于分步傅里叶算法的叠前时间偏移方法
CN104199088B (zh) 一种提取入射角道集的方法及***
CN111505712A (zh) 地震层速度计算方法及装置
Tessmer Reverse-time migration in TTI media
WO2014149561A1 (en) Migration velocity analysis method for vertical seismic profile data
CN107561583A (zh) 用于高斯束叠前深度偏移的局部角度计算方法及成像方法
CN113866827A (zh) 一种解释性速度建模地震成像方法、***、介质和设备
CN104502972B (zh) 三分量地震波积分偏移方法和装置
CN104101901A (zh) 一种转换波弯曲射线保幅各向异性叠前时间偏移方法
US20170336523A1 (en) Seismic signal processing method, apparatus and system
CN103424775B (zh) 基于地震波深度偏移理论的近地表多震源点位置定位方法
CN113126152B (zh) 深度域速度模型构建方法及装置
O'Brien Mounting an offense against poor-quality shear data
CN113093279B (zh) 转换波静校正方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20200807

RJ01 Rejection of invention patent application after publication