CN104199088A - 一种提取入射角道集的方法及*** - Google Patents
一种提取入射角道集的方法及*** Download PDFInfo
- Publication number
- CN104199088A CN104199088A CN201410403449.2A CN201410403449A CN104199088A CN 104199088 A CN104199088 A CN 104199088A CN 201410403449 A CN201410403449 A CN 201410403449A CN 104199088 A CN104199088 A CN 104199088A
- Authority
- CN
- China
- Prior art keywords
- dtri
- wave field
- module
- shot point
- point wave
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例公开了一种提取入射角道集的方法及***,所述方法包括:读入子波数据、地震炮集数据和层速度模型数据,计算预设区域内的炮点波场和检波点波场;根据所述炮点波场和检波点波场,计算地震炮集的地震偏移成像结果,并基于所述地震偏移成像结果得到三维数据;基于所述炮点波场和预设频率值计算单频炮点波场;基于所述单频炮点波场计算入射方向矢量,基于所述三维数据计算法向矢量;根据所述入射方向矢量和所述法向矢量,计算所述预设区域中坐标位置对应的入射角,并将所述位置的偏移成像能量值映射至所述入射角形成的角度域,映射完成后得到入射角道集。本申请公开的提取入射角道集的方法及***可以提高提取的入射角道集的精度。
Description
技术领域
本申请涉及地震勘探技术领域,特别涉及一种提取入射角道集的方法及***。
背景技术
地震勘探是指利用地下介质弹性和密度的差异,通过观测和分析大地对人工激发地震波的响应,推断地下岩层的性质和形态的地球物理勘探方法。入射角道集是地震勘探过程中非常重要的地震数据,所述入射角道集可以包括:地层中的坐标位置、所述坐标位置对应的入射角,以及所述坐标位置的波场能量值。所述入射角道集对于地层中岩性分析、油藏分析一会波场传播建模都具有重要的意义。
现有的提取入射角道集的方法包括:
在波动方程偏移的成像阶段,以成像点为中心,沿水平方向分别以不同的半偏移距进行错位成像。根据错位成像的结果生成地下偏移距道集,然后使用拉冬变换将所述地下偏移道集转换为入射角道集。
在实现本申请过程中,发明人发现现有技术中至少存在如下问题:现有的提取入射角道集的方法直接根据错位成像的结果生成地下偏移道集,没有考虑地层倾角,该方法只适合应用于水平地层。但地震勘探过程中,地形一般比较复杂,通常地层都存在倾斜的情况,采用该方法在复杂地形中提取入射角道集时,由于没有考虑地层倾角,计时得到的入射角存在偏差,因此该方法产生的入射角道集计算精度不高。
发明内容
本申请实施例的目的是提供一种提取入射角道集的方法及***,以提高提取的入射角道集的精度。
为解决上述技术问题,本申请实施例提供一种提取入射角道集的方法及***是这样实现的:
一种提取入射角道集的方法,包括:
读入子波数据、地震炮集数据和层速度模型数据,根据子波数据和层速度模型数据计算预设区域内的炮点波场,根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场;
根据所述炮点波场和检波点波场,计算地震炮集的地震偏移成像结果,并基于所述地震偏移成像结果得到三维数据;
基于所述炮点波场和预设频率值计算单频炮点波场;
利用结构张量倾角扫描方法,基于所述单频炮点波场计算入射方向矢量,基于所述三维数据计算法向矢量;
根据所述入射方向矢量和所述法向矢量,计算所述预设区域中坐标位置对应的入射角,并将所述坐标位置的偏移成像能量值映射至所述入射角形成的角度域,映射完成后得到入射角道集。
优选方案中,所述根据子波数据和层速度模型数据计算预设区域内的炮点波场,包括:采用下述公式进行计算
上式中,P是需要计算的炮点波场;(x,y,z)为预设区域内的坐标,是变量;t为时间,是变量;S(t)是子波数据。
优选方案中,所述根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场,包括:采用下述公式进行计算
上式中,P是需要计算的检波点波场;(x,y,z)为预设区域内的坐标,是变量;t为时间,是变量;S(t)取地震炮集数据。
优选方案中,所述基于所述炮点波场和预设频率值计算单频炮点波场,包括:
将所述炮点波场乘以指数函数e-iωt得到乘积,其中ω=2πf,f是单频频率;
将计算的乘积的实部设置为单频炮点波场。
优选方案中,所述单频频率f的取值范围包括:15HZ~20HZ。
优选方案中,所述利用结构张量倾角扫描方法,基于所述单频炮点波场计算入射方向矢量,包括:
利用下述公式对所述单频炮点波场构造张量:
上式中,T表示构造张量;表示的转置;表示I的梯度;I表示单频炮点波场;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量;
对所述构造张量进行平滑,具体地,将所述构造张量T乘以平滑窗口算子S,如下述公式所示:
上式中,S表示平滑窗口算子;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量;
求解公式中矩阵ST的特征值,所得特征值包括:x方向特征值、y方向特征值和z方向特征值;
X方向特征值除以z方向特征值得到xdip1,y方向特征值除以z方向特征值得到ydip1;将所述xdip1和ydip1设置为入射方向矢量。
优选方案中,所述利用结构张量倾角扫描方法,基于所述三维数据计算法向矢量,包括:
利用下述公式对所述三维数据构造张量:
上式中,T表示构造张量;表示的转置;表示I的梯度;I表示三维数据;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量;
对所述构造张量进行平滑,具体地,将所述构造张量T乘以平滑窗口算子S,如下述公式所示:
上式中,S表示平滑窗口算子;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量;
求解公式中矩阵ST的特征值,所得特征值包括:x方向特征值、y方向特征值和z方向特征值;
X方向特征值除以z方向特征值得到xdip2,y方向特征值除以z方向特征值得到ydip2;将所述xdip2和ydip2设置为法向矢量。
优选方案中,所述入射角道集包括:坐标位置、所述坐标位置对应的入射角,以及所述坐标位置的能量值。
一种提取入射角道集的***,包括:波场计算模块、成像模块、单频炮点波场计算模块、矢量计算模块和入射角道集模块;其中,
所述波场计算模块,用于读入子波数据、地震炮集数据和层速度模型数据,根据子波数据和层速度模型数据计算预设区域内的炮点波场,根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场;
所述成像模块,用于根据所述炮点波场和检波点波场,计算地震炮集的地震偏移成像结果,并基于所述地震炮集的地震偏移成像结果得到三维数据;
所述单频炮点波场计算模块,用于基于所述炮点波场和预设频率值计算单频炮点波场;
所述矢量计算模块,用于利用结构张量倾角扫描方法,基于所述单频炮点波场计算入射方向矢量,基于所述三维数据计算法向矢量;
所述入射角道集模块,用于根据所述入射方向矢量和所述法向矢量,计算所述预设区域中坐标位置对应的入射角,并将所述坐标位置的偏移成像能量值映射至所述入射角形成的角度域,映射完成后提取得到入射角道集。
优选方案中,所述波场计算模块,包括:读入模块、炮点波场模块和检波点波场模块;其中,
所述读入模块,用于读入子波数据、地震炮集数据和层速度模型数据;
所述炮点波场模块,用于根据子波数据和层速度模型数据计算预设区域内的炮点波场;
所述检波点波场模块,用于根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场。
优选方案中,所述成像模块,包括:成像计算模块和三维数据模块;其中,
所述成像计算模块,用于根据所述波场计算模块中计算所得的炮点波场和检波点波场,计算地震炮集的地震偏移成像结果;
所述三维数据模块,用于根据成像计算模块得到的地震偏移成像结果得到三维数据。
优选方案中,所述矢量计算模块,包括:入射方向矢量模块和法向矢量模块;其中,
所述入射方向矢量模块,用于利用结构张量倾角扫描方法,基于单频炮点波场计算模块,所得的单频炮点波场计算入射方向矢量;
所述法向矢量模块,用于利用结构张量倾角扫描方法,基于所述成像模块,中得到的三维数据计算法向矢量。
优选方案中,所述入射角道集模块,包括:入射角计算模块、映射模块和道集提取模块;其中,
所述入射角计算模块,用于根据矢量计算模块计算得到的入射方向矢量和法向矢量计算所述预设区域中坐标位置对应的入射角;
所述映射模块,用于将所述坐标位置的偏移成像能量值映射至所述入射角计算模块计算的入射角所形成的角度域;
所述道集提取模块,用于根据映射模块的结果提取入射角道集。
由以上本申请实施例提供的技术方案可见,本申请实施例公开的提取入射角道集的方法及***,根据可以反映出波场的传播方向的单频炮点波场来计算入射方向矢量,根据三维数据计算法向矢量,再根据所述入射方向矢量和法向矢量来计算入射角。所述法向矢量可以表示地层的倾角。因此,采用本申请公开的方法计算出的入射角包含了复杂地形中的地层倾角,得到的入射角精度较高,从而可以提高提取的入射角道集的精度。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请一种提取入射角道集的方法实施例的流程图;
图2是本申请一种提取入射角道集的***实施例的模块图;
图3是本申请***实施例中波场计算模块的模块图;
图4是本申请***实施例中成像模块的模块图;
图5是本申请***实施例中矢量计算模块的模块图;
图6是本申请***实施例中入射角道集模块的模块图。
具体实施方式
本申请实施例提供一种提取入射角道集的方法及***。
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1是本申请一种提取入射角道集的方法实施例的流程图。如图1所示,所述提取入射角道集的方法可以包括:
S101:读入子波数据、地震炮集数据和层速度模型数据,根据子波数据和层速度模型数据计算预设区域内的炮点波场,根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场。
计算机可以读入子波数据、地震炮集数据和层速度模型数据。
所述子波数据可以用于模拟震源。所述子波数据中可以包括:位置参数(x,y)、时间参数t和频率F。所述频率F的取值范围可以是15HZ~25HZ。所述子波的个数可以与所述地震炮集数据中的炮数相同。
所述地震炮集数据可以是按照炮号排列的地震反射数据,所述地震反射数据可以是经过预处理的,例如是经过去噪等预处理的反射数据。所述地震炮集数据中的参数可以包括:位置参数(x,y)、时间参数t和频率F。
所述层速度模型用于描述地下岩层的纵波传播速度。所述层速度模型中的参数可以包括:位置参数(x,y)、深度参数z、和速度参数v。
根据子波数据和层速度模型数据,可以计算预设区域内的炮点波场。所述预设区域包括预设的位置和预设的深度。
所述计算预设区域内的炮点波场可以通过声波方程来实现,所述声波方程如公式(1)所示:
公式(1)中,P表示波场;(x,y,z)表示空间坐标,对应于预设区域内的位置参数(x,y)和深度参数z;t表示时间,v表示传播速度;S(t)表示波数据。
利用公式(1)计算炮点波场时,公式(1)中,P是需要计算的炮点波场;(x,y,z)为预设区域内的坐标,是变量;t为时间,是变量;S(t)取子波数据。
根据地震炮集数据和层速度模型数据还可以计算所述预设区域内的检波点波场。所述计时预设区域内的检波点波场也可以利用公式(1)计算得到。
利用公式(1)计算检波点波场是,公式(1)中,P是需要计算的检波点波场;(x,y,z)为预设区域内的坐标,是变量;t为时间,是变量;S(t)取地震炮集数据。
S102:根据所述炮点波场和检波点波场,计算地震炮集的地震偏移成像结果,并基于所述地震偏移成像结果得到三维数据。
根据计算得到的炮点波场和检波点波场,可以计算地震炮集的地震偏移成像结果。具体可以将所述地震炮集记录中每一炮在同一时刻的炮点波场和检波点波场进行互相关成像计算,得到地震偏移成像结果。
基于所述地震偏移成像结果可以得到三维数据,具体地,将炮集中每一炮计算所得的地震偏移成像结果进行叠加,所得的叠加结果即为三维数据。
S103:基于所述炮点波场和预设频率值计算单频炮点波场。
基于计算得到的炮点波场和预设的频率值,可以计算单频炮点波场。具体可以包括:
将计算所得的炮点波场乘以指数函数e-iωt得到乘积,其中ω=2πf,f是单频频率,所述f的取值范围可以包括:15HZ~20HZ。
将上述乘积的实部保留下来,可以将所述乘积的实部设置为单频炮点波场。所述单频炮点波场保存可以反映出波场的传播方向。
S104:利用结构张量倾角扫描方法,基于所述单频炮点波场计算入射方向矢量,基于所述三维数据计算法向矢量。
基于所述单频炮点波场可以计算入射方向矢量,所述入射方向矢量可以利用结构张量倾角扫描方法获得。
设所述单频炮点波场为I,则可以定义一构造张量T:
其中,表示的转置;表示I的梯度;所述具体如下式所示:
公式(3)中,表示x方向的变化量;表示y方向的变化量;表示z方向的变化量。
所述构造张量T可以用于表示区域的变化方向和沿变化方向的变化量大小,其中特征向量反映了局部区域变化的方向,特征值反映了变化的大小。
对所述单频炮点波场构造张量,可以如下式所示:
公式(4)中,T表示构造张量;表示的转置;表示I的梯度;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量。
可以对公式(4)的构造张量进行平滑,具体地,可以通过将所述构造张量T乘以平滑窗口算子S来实现,平滑后所得结果可以如公式(5)所示:
公式(5)中,S表示平滑窗口算子;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量。
求解公式(5)中矩阵的特征值,所得特征值可以包括:x方向特征值、y方向特征值和z方向特征值。
X方向特征值除以z方向特征值可以得到xdip1,y方向特征值除以z方向特征值可以得到ydip1。所述xdip1和ydip1可以用于表示入射方向矢量。
基于所述三维数据可以计算法向矢量,所述法向矢量也可以利用结构张量倾角扫描方法获得,计算方法与入射方向矢量的计算方法相同。计算得到的法向方向可以用xdip2和ydip2表示。所述法向矢量可以反映地质构造的地层倾角。
S105:根据所述入射方向矢量和所述法向矢量,计算所述预设区域中坐标位置对应的入射角,并将所述坐标位置的偏移成像能量值映射至所述入射角形成的角度域,映射完成后得到入射角道集。
根据所述计算得到的入射方向矢量和所述方向矢量,可以计算所述预设区域中每一坐标位置对应的入射角。具体地,(xdip1-xdip2)和(ydip1-ydip2)可以用于表示预设区域中某一位置的入射角。可以设置所述入射角形成的区域为角度域。可以将所述位置的能量值映射至所述角度域。例如所述预设区域中第一位置对应的入射角可以为第一入射角,可以将所述第一入射角形成的第一角度域中的能量值叠加上所述第一位置的偏移成像能量值。当所述预设区域中每一坐标位置的入射角都完成能量值的映射后,可以得到入射角道集。所述入射角道集可以包括:映射完成后,所述预设区域内每一坐标位置、所述坐标位置对应的入射角,以及所述坐标位置的能量值。
上述实施例公开的提取入射角道集的方法实施例,根据可以反映出波场的传播方向的单频炮点波场来计算入射方向矢量,根据三维数据计算法向矢量。再根据所述入射方向矢量和法向矢量来计算入射角。所述法向矢量可以表示地层的倾角。因此采用本申请公开的方法计算出的入射角包含了复杂地形中的地层倾角,得到的入射角精度较高,从而可以提高提取的入射角道集的精度。
图2是本申请一种提取入射角道集的***实施例的模块图。如图2所示,所述提取入射角道集的***可以包括:波场计算模块100、成像模块200、单频炮点波场计算模块300、矢量计算模块400、入射角道集模块500。其中,
所述波场计算模块100,可以用于读入子波数据、地震炮集数据和层速度模型数据,根据子波数据和层速度模型数据计算预设区域内的炮点波场,根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场。
所述成像模块200,可以用于根据所述炮点波场和检波点波场,计算地震炮集的地震偏移成像结果,并基于所述地震炮集的地震偏移成像结果得到三维数据。
所述单频炮点波场计算模块300,可以用于基于所述炮点波场和预设频率值计算单频炮点波场。
所述矢量计算模块400,可以用于利用结构张量倾角扫描方法,基于所述单频炮点波场计算入射方向矢量,基于所述三维数据计算法向矢量。
所述入射角道集模块500,可以用于根据所述入射方向矢量和所述法向矢量,计算所述预设区域中坐标位置对应的入射角,并将所述坐标位置的偏移成像能量值映射至所述入射角形成的角度域,映射完成后提取得到入射角道集。
图3是本申请***实施例中波场计算模块的模块图。如图3所示,所述波场计算模块100,可以包括:读入模块110、炮点波场模块120和检波点波场模块130。其中,
所述读入模块110,可以用于读入子波数据、地震炮集数据和层速度模型数据。
所述炮点波场模块120,可以用于根据子波数据和层速度模型数据计算预设区域内的炮点波场。
所述检波点波场模块130,可以用于根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场。
图4是本申请***实施例中成像模块的模块图。如图4所示,所述成像模块200,可以包括:成像计算模块210和三维数据模块220。其中,
所述成像计算模块210,可以用于根据所述波场计算模块100中计算所得的炮点波场和检波点波场,计算地震炮集的地震偏移成像结果。
所述三维数据模块220,可以用于根据成像计算模块210得到的地震偏移成像结果得到三维数据。
图5是本申请***实施例中矢量计算模块的模块图。如图5所示,所述矢量计算模块400,可以包括:入射方向矢量模块410和法向矢量模块420。其中,
所述入射方向矢量模块410,可以用于利用结构张量倾角扫描方法,基于单频炮点波场计算模块300所得的单频炮点波场计算入射方向矢量。
所述法向矢量模块420,可以用于利用结构张量倾角扫描方法,基于所述成像模块200中得到的三维数据计算法向矢量。
图6是本申请***实施例中入射角道集模块的模块图。如图6所示,所述入射角道集模块500,可以包括:入射角计算模块510、映射模块520和道集提取模块530。其中,
所述入射角计算模块510,可以用于根据矢量计算模块400计算得到的入射方向矢量和法向矢量计算所述预设区域中坐标位置对应的入射角。
所述映射模块520,可以用于将所述坐标位置的偏移成像能量值映射至所述入射角计算模块510计算的入射角所形成的角度域。
所述道集提取模块530,可以用于根据映射模块520的结果提取入射角道集。
上述公开的提取入射角道集的***实施例与本申请公开的提取入射角道集的方法实施例相对应,可以实现本申请中公开的方法实施例。
需要说明的是,本说明书中,对于***实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (13)
1.一种提取入射角道集的方法,其特征在于,包括:
读入子波数据、地震炮集数据和层速度模型数据,根据子波数据和层速度模型数据计算预设区域内的炮点波场,根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场;
根据所述炮点波场和检波点波场,计算地震炮集的地震偏移成像结果,并基于所述地震偏移成像结果得到三维数据;
基于所述炮点波场和预设频率值计算单频炮点波场;
利用结构张量倾角扫描方法,基于所述单频炮点波场计算入射方向矢量,基于所述三维数据计算法向矢量;
根据所述入射方向矢量和所述法向矢量,计算所述预设区域中坐标位置对应的入射角,并将所述坐标位置的偏移成像能量值映射至所述入射角形成的角度域,映射完成后得到入射角道集。
2.如权利要求1所述的一种提取入射角道集的方法,其特征在于,所述根据子波数据和层速度模型数据计算预设区域内的炮点波场,包括:采用下述公式进行计算
上式中,P是需要计算的炮点波场;(x,y,z)为预设区域内的坐标,是变量;t为时间,是变量;S(t)是子波数据。
3.如权利要求1所述的一种提取入射角道集的方法,其特征在于,所述根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场,包括:采用下述公式进行计算
上式中,P是需要计算的检波点波场;(x,y,z)为预设区域内的坐标,是变量;t为时间,是变量;S(t)取地震炮集数据。
4.如权利要求1所述的一种提取入射角道集的方法,其特征在于,所述基于所述炮点波场和预设频率值计算单频炮点波场,包括:
将所述炮点波场乘以指数函数e-iωt得到乘积,其中ω=2πf,f是单频频率;
将计算的乘积的实部设置为单频炮点波场。
5.如权利要求4所述的一种提取入射角道集的方法,其特征在于,所述单频频率f的取值范围包括:15HZ~20HZ。
6.如权利要求1所述的一种提取入射角道集的方法,其特征在于,所述利用结构张量倾角扫描方法,基于所述单频炮点波场计算入射方向矢量,包括:
利用下述公式对所述单频炮点波场构造张量:
上式中,T表示构造张量;表示的转置;表示I的梯度;I表示单频炮点波场;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量;
对所述构造张量进行平滑,具体地,将所述构造张量T乘以平滑窗口算子S,如下述公式所示:
上式中,S表示平滑窗口算子;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量;
求解公式中矩阵ST的特征值,所得特征值包括:x方向特征值、y方向特征值和z方向特征值;
X方向特征值除以z方向特征值得到xdip1,y方向特征值除以z方向特征值得到ydip1;将所述xdip1和ydip1设置为入射方向矢量。
7.如权利要求1所述的一种提取入射角道集的方法,其特征在于,所述利用结构张量倾角扫描方法,基于所述三维数据计算法向矢量,包括:
利用下述公式对所述三维数据构造张量:
上式中,T表示构造张量;表示的转置;表示I的梯度;I表示三维数据;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量;
对所述构造张量进行平滑,具体地,将所述构造张量T乘以平滑窗口算子S,如下述公式所示:
上式中,S表示平滑窗口算子;表示x方向的变化量;表示y方向的变化量;表示z方向的变化量;
求解公式中矩阵ST的特征值,所得特征值包括:x方向特征值、y方向特征值和z方向特征值;
X方向特征值除以z方向特征值得到xdip2,y方向特征值除以z方向特征值得到ydip2;将所述xdip2和ydip2设置为法向矢量。
8.如权利要求1-7中任意一项所述的一种提取入射角道集的方法,其特征在于,所述入射角道集包括:坐标位置、所述坐标位置对应的入射角,以及所述坐标位置的能量值。
9.一种提取入射角道集的***,其特征在于,包括:波场计算模块、成像模块、单频炮点波场计算模块、矢量计算模块和入射角道集模块;其中,
所述波场计算模块,用于读入子波数据、地震炮集数据和层速度模型数据,根据子波数据和层速度模型数据计算预设区域内的炮点波场,根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场;
所述成像模块,用于根据所述炮点波场和检波点波场,计算地震炮集的地震偏移成像结果,并基于所述地震炮集的地震偏移成像结果得到三维数据;
所述单频炮点波场计算模块,用于基于所述炮点波场和预设频率值计算单频炮点波场;
所述矢量计算模块,用于利用结构张量倾角扫描方法,基于所述单频炮点波场计算入射方向矢量,基于所述三维数据计算法向矢量;
所述入射角道集模块,用于根据所述入射方向矢量和所述法向矢量,计算所述预设区域中坐标位置对应的入射角,并将所述坐标位置的偏移成像能量值映射至所述入射角形成的角度域,映射完成后提取得到入射角道集。
10.如权利要求9所述的一种提取入射角道集的***,其特征在于,所述波场计算模块,包括:读入模块、炮点波场模块和检波点波场模块;其中,
所述读入模块,用于读入子波数据、地震炮集数据和层速度模型数据;
所述炮点波场模块,用于根据子波数据和层速度模型数据计算预设区域内的炮点波场;
所述检波点波场模块,用于根据地震炮集数据和层速度模型数据计算所述预设区域内的检波点波场。
11.如权利要求9所述的一种提取入射角道集的***,其特征在于,所述成像模块,包括:成像计算模块和三维数据模块;其中,
所述成像计算模块,用于根据所述波场计算模块中计算所得的炮点波场和检波点波场,计算地震炮集的地震偏移成像结果;
所述三维数据模块,用于根据成像计算模块得到的地震偏移成像结果得到三维数据。
12.如权利要求9所述的一种提取入射角道集的***,其特征在于,所述矢量计算模块,包括:入射方向矢量模块和法向矢量模块;其中,
所述入射方向矢量模块,用于利用结构张量倾角扫描方法,基于单频炮点波场计算模块,所得的单频炮点波场计算入射方向矢量;
所述法向矢量模块,用于利用结构张量倾角扫描方法,基于所述成像模块,中得到的三维数据计算法向矢量。
13.如权利要求9所述的一种提取入射角道集的***,其特征在于,所述入射角道集模块,包括:入射角计算模块、映射模块和道集提取模块;其中,
所述入射角计算模块,用于根据矢量计算模块计算得到的入射方向矢量和法向矢量计算所述预设区域中坐标位置对应的入射角;
所述映射模块,用于将所述坐标位置的偏移成像能量值映射至所述入射角计算模块计算的入射角所形成的角度域;
所述道集提取模块,用于根据映射模块的结果提取入射角道集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410403449.2A CN104199088B (zh) | 2014-08-15 | 2014-08-15 | 一种提取入射角道集的方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410403449.2A CN104199088B (zh) | 2014-08-15 | 2014-08-15 | 一种提取入射角道集的方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104199088A true CN104199088A (zh) | 2014-12-10 |
CN104199088B CN104199088B (zh) | 2017-02-08 |
Family
ID=52084400
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410403449.2A Active CN104199088B (zh) | 2014-08-15 | 2014-08-15 | 一种提取入射角道集的方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199088B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105372704A (zh) * | 2015-10-14 | 2016-03-02 | 中国石油天然气集团公司 | 一种获取地震波传播方向的方法及装置 |
CN107728196A (zh) * | 2016-08-10 | 2018-02-23 | 中国石油化工股份有限公司 | 获取角度域共成像道集的方法及*** |
CN107957234A (zh) * | 2016-10-17 | 2018-04-24 | 中国航空工业集团公司北京航空制造工程研究所 | 一种测量自由曲面任意点处法向矢量的方法及装置 |
CN110646840A (zh) * | 2018-06-27 | 2020-01-03 | 中国石油化工股份有限公司 | 角道集提取方法及*** |
CN112513925A (zh) * | 2018-08-01 | 2021-03-16 | 皇家飞利浦有限公司 | 为ct虚拟单色提供自动自适应能量设置的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
CN102841375A (zh) * | 2012-09-06 | 2012-12-26 | 中国石油大学(华东) | 一种复杂条件下基于角度域共成像点道集的层析速度反演方法 |
CN102901984A (zh) * | 2012-09-29 | 2013-01-30 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 真地表地震数据倾角道集构建方法 |
CN103472481A (zh) * | 2012-06-06 | 2013-12-25 | 中国石油化工股份有限公司 | 一种利用gpu进行逆时偏移提取角度道集的方法 |
-
2014
- 2014-08-15 CN CN201410403449.2A patent/CN104199088B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
CN103472481A (zh) * | 2012-06-06 | 2013-12-25 | 中国石油化工股份有限公司 | 一种利用gpu进行逆时偏移提取角度道集的方法 |
CN102841375A (zh) * | 2012-09-06 | 2012-12-26 | 中国石油大学(华东) | 一种复杂条件下基于角度域共成像点道集的层析速度反演方法 |
CN102901984A (zh) * | 2012-09-29 | 2013-01-30 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 真地表地震数据倾角道集构建方法 |
Non-Patent Citations (2)
Title |
---|
刘礼农 等: "波动方程叠前深度偏移直接产生角道集", 《地球物理学报》 * |
王保利 等: "逆时偏移中用Poynting矢量高效地提取角道集", 《地球物理学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105372704A (zh) * | 2015-10-14 | 2016-03-02 | 中国石油天然气集团公司 | 一种获取地震波传播方向的方法及装置 |
CN105372704B (zh) * | 2015-10-14 | 2018-06-01 | 中国石油天然气集团公司 | 一种获取地震波传播方向的方法及装置 |
CN107728196A (zh) * | 2016-08-10 | 2018-02-23 | 中国石油化工股份有限公司 | 获取角度域共成像道集的方法及*** |
CN107957234A (zh) * | 2016-10-17 | 2018-04-24 | 中国航空工业集团公司北京航空制造工程研究所 | 一种测量自由曲面任意点处法向矢量的方法及装置 |
CN107957234B (zh) * | 2016-10-17 | 2023-09-29 | 中国航空制造技术研究院 | 一种测量自由曲面任意点处法向矢量的方法及装置 |
CN110646840A (zh) * | 2018-06-27 | 2020-01-03 | 中国石油化工股份有限公司 | 角道集提取方法及*** |
CN110646840B (zh) * | 2018-06-27 | 2021-07-20 | 中国石油化工股份有限公司 | 角道集提取方法及*** |
CN112513925A (zh) * | 2018-08-01 | 2021-03-16 | 皇家飞利浦有限公司 | 为ct虚拟单色提供自动自适应能量设置的方法 |
CN112513925B (zh) * | 2018-08-01 | 2024-03-08 | 皇家飞利浦有限公司 | 图像处理***、成像装置和图像处理的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104199088B (zh) | 2017-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Maraschini et al. | A new misfit function for multimodal inversion of surface waves | |
Claerbout | Synthesis of a layered medium from its acoustic transmission response | |
Bormann et al. | Seismic wave propagation and earth models | |
US20140129479A1 (en) | Method to aid in the exploration, mine design, evaluation and/or extraction of metalliferous mineral and/or diamond deposits | |
Bednar | A brief history of seismic migration | |
EP3710867B1 (en) | Noise attenuation of multiple source seismic data | |
CN106094029A (zh) | 利用偏移距矢量片地震数据预测储层裂缝的方法 | |
CN103635831A (zh) | 在3d各向异性介质中的基于敏感度内核的偏移速度分析 | |
CN102053263B (zh) | 调查表层结构的方法 | |
CN107817526B (zh) | 叠前地震道集分段式振幅能量补偿方法及*** | |
CN104635261B (zh) | 山前带重磁电震联合分体建模方法 | |
CN102053261A (zh) | 一种地震数据处理方法 | |
CN104199088A (zh) | 一种提取入射角道集的方法及*** | |
CN102636809B (zh) | 一种传播角度域共成像点道集的生成方法 | |
CN104459777A (zh) | 基于流体体积模量avo反演的流体识别的方法及*** | |
Harding et al. | Velocity structure near IODP Hole U1309D, Atlantis Massif, from waveform inversion of streamer data and borehole measurements | |
CN103758511A (zh) | 一种井下逆时偏移成像识别隐蔽储层的方法及装置 | |
CN103926623A (zh) | 一种压制逆时偏移低频噪音的方法 | |
CN110850469A (zh) | 一种基于克希霍夫积分解的地震槽波深度偏移的成像方法 | |
CN113740906A (zh) | 一种水下垂直缆地震波干涉成像方法及装置 | |
Liu et al. | Full-waveform inversion method for tunnel seismic forward prospecting | |
Tonegawa et al. | Geographical distribution of shear wave anisotropy within marine sediments in the northwestern Pacific | |
CN105527648A (zh) | 用于各向异性参数反演的敏感度矩阵的计算方法及*** | |
CN103424775B (zh) | 基于地震波深度偏移理论的近地表多震源点位置定位方法 | |
Ensing et al. | A 3D crustal shear speed model of the Auckland volcanic field, New Zealand, from multi-component ambient noise tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |