CN102282481B - 基于地震能见度分析的数据采集和叠前偏移 - Google Patents
基于地震能见度分析的数据采集和叠前偏移 Download PDFInfo
- Publication number
- CN102282481B CN102282481B CN200980154752.0A CN200980154752A CN102282481B CN 102282481 B CN102282481 B CN 102282481B CN 200980154752 A CN200980154752 A CN 200980154752A CN 102282481 B CN102282481 B CN 102282481B
- Authority
- CN
- China
- Prior art keywords
- visibility
- offset
- seismic
- object event
- receiver
- 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.)
- Expired - Fee Related
Links
- 238000013508 migration Methods 0.000 title claims abstract description 17
- 230000005012 migration Effects 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 48
- 238000003384 imaging method Methods 0.000 claims abstract description 19
- 238000004088 simulation Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 2
- 230000002708 enhancing effect Effects 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 230000006870 function Effects 0.000 description 13
- 238000012545 processing Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 238000013459 approach Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000005452 bending Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 241000193935 Araneus diadematus Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- 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/67—Wave propagation modeling
-
- 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/67—Wave propagation modeling
- G01V2210/673—Finite-element; Finite-difference
-
- 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/67—Wave propagation modeling
- G01V2210/675—Wave equation; Green's functions
-
- 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/67—Wave propagation modeling
- G01V2210/679—Reverse-time modeling or coalescence modelling, i.e. starting from receivers
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
利用所选地下结构的地震能见度分析确定能够提供目标事件的高能见度的表面位置。然后这些位置可用作采集额外的地震勘测数据的基础和/或为更复杂的偏移方法来选择用于再偏移的现有轨迹的基础。对于任一用法,期望使用新偏移数据提供目标事件的增强图像。在一些实施例中,能见度的确定包括:针对多个模拟炮点的每一个,使用基于波动方程的传播函数在地震模型中找到目标事件的反射波场;以及针对多个接收器位置的每一个,计算每个反射波场的贡献信号。能见度的确定过程还进一步包括将每一个贡献信号转化为源-接收器能见度值。因为数据采集和/或再偏移被限制在所选区域中,所以极大的减少了目标事件成像的工作。
Description
技术领域
本发明涉及地震能见度分析,更具体地,涉及基于地震能见度分析的数据采集和叠前偏移(prestack migration)。
背景技术
科学家和工程师通常使用地震勘测来进行勘探、地球物理学研究以及工程规划。地震勘测可以提供有关地下结构的信息,包括地层边界、岩石类型以及是否存在流体储层。上述信息对搜寻水、地热储层、以及例如油气和矿石的矿物沉积层提供了极大的帮助。尤其是石油公司经常投资于昂贵的地震勘测来选择用于勘探油井的地点。
传统的地震勘测使用了人工地震能量源,例如炮弹、***或者能够产生地震波的震动源。上述这些源在点火时,会制造地震“事件”,即如地震波一样从源处向下传入地球的地震能量脉冲。位于不同地层之间的断层和边界产生不同的声阻抗,这引起地震波的部分反射。地震传感器阵列可以检测并记录这些反射来用于后续分析。而后对记录的信号应用复杂的处理技术来提取地下(subsurface)结构的图像或其他表示。
遗憾的是,地震分析人员经常发现地下特征成像质量差或无法充分的辨识。在这种情况下,唯一的解决方案要么是追寻更加复杂的处理技术,要么是在先前勘测的地区采集额外的数据。这些方案的每一个在时间以及金钱上开销过大。
发明内容
因此,本发明于此公开了用于对所选择的地下结构进行地震能见度分析的***与方法。所述***与方法确定出能够最好的展示地下结构的细节的地震源位置和接收器位置。然后这些位置可以作为采集额外地震数据的基础和/或对所选择的现有数据的子集进行更复杂数据处理的基础。因为采集和/或处理的数据区域大大减少了,所以相关的开销也最小化了。
一些示例的方法实施例包括了一种地震勘测方法,其包括:确定目标事件的能见度,能见度根据地震源位置和接收器位置而变化;以及在所选区域中采集地震数据,所述所选区域应至少部分包括能见度值高于阈值的位置。然后基于新采集的地震数据对目标事件成像。示例的方法实施例还包括一种地震偏移方法,其包括:在现有地震勘测中的轨迹的源位置和接收器位置来确定目标事件的能见度;以及将能见度值高于阈值的轨迹进行再偏移以对目标事件成像。在上述两种情况下,能见度的确定可以包括针对多个模拟炮点的每一个,使用基于波动方程的传播函数在地震模型中寻找目标事件的反射波场;以及针对多个接收器位置的每一个,计算来自从每个反射波场的贡献信号。能见度的确定还可进一步包括将贡献信号转化为源-接收器能见度值。
附图说明
根据附图以及详细描述将对多个公开的实施例有更好的了解:
图1为示例的地震勘测环境;
图2为示例的地震源和接收器的几何分布;
图3为示例的地震轨迹;
图4为对示例地震模型的偏移深度图像;
图5a显示了从源到目标的地震能量传播;
图5b为接收器能见度的示例图;
图6a-6c为在不同假设下源能见度的示例图;
图7为示例的源-接收器能见度函数;
图8为示例的对地震数据的偏移子集的偏移深度图像;
图9为示例的能见度分析方法的流程图;以及
图10为示例的能见度分析***。
由于本发明适用于多种修改及替代形式,因此其具体的实施例通过附图举例并进行了详细说明。然而,应理解附图及详细描述不意指将本发明限制于所公开的特殊形式,而是相反,本发明覆盖落入到权利要求范围中的所有修改、等同或替代形式。
具体实施方式
本发明公开了多种能见度分析方法和***,所述方法和***确定能够最优的测量一个或多个所选地下特征的属性的地震勘测源和/或接收器位置。分析人员从而可以对这些区域进行采集与处理工作来改善这些所选特征的成像细节。借助示例用法内容中的描述,可以最好地理解公开的***和方法。
相应的,图1为示例的地震勘测环境,其中,地震接收器102的阵列间隔布置在地球表面104来探测地震波。接收器102经无线或电缆连接至数据采集单元106,所述数据采集单元106接收、处理并存储由接收器收集的地震信号数据。地震能量源108(例如震源车)在多个位置被触发以产生地震能量波,地震能量波通过地球110传播并且从声阻抗不连续处反射至接收器102。上述声阻抗不连续处可能是由断层、地层层床之间的边界、或地层流层间的边界造成的。这些不连续处在由地震信号数据生成的地下结构表示中将显示为亮点。
图1进一步表示了一示例地下模型,其可以用作本公开的一个例子。在该模型中,地球具有四个相互关联的平面地层,并且在第三层和第四层之间具有陡峭弯曲的边界。在每一层中的声速由上至下分别为2000m/s,3500m/s,2800m/s和4000m/s。作为举例应注意到的是,图1并没有按比例绘制。接收器阵列的水平宽度方向应限制为几千米或更少,而地下模型的水平平尺寸为16千米,垂直尺寸为5千米。
图2为源位置202和接收器位置204的示例几何分布,可看做是俯视图。可行的地震勘测的几何分布具有无穷的变化形式,并且可基于区块-区块、滚动的几何分布、移动阵列、及其各种组合来聚集。在这里,主要信息是响应于地震源的每次点火所获得的接收器信号的数量相当大,并且当考虑源的点火次数时,所获得的轨迹很容易到达几百万条。
图3为示例的接收到的地震信号轨迹S0-S2。这些轨迹可以表示位移、速度、加速度、压力或为时间的函数的地震能量强度的其他测量值。通常在每次点火之后,以每秒大约500个样本的速率持续大约30秒钟,将每个接收器接收的信号采样和数字化为8至32比特的分辨率。在一些情况下,接收器感测到多个分量的数据,这为每一条轨迹进一步增加了样本数据量。轨迹数据在存储前可进行滤波或压缩。存储的地震勘测数据被传送或通信至数据处理装置。
数据处理装置处的计算机网络利用已知技术来处理数据以估计声速的容积分布。参见,例如Jon F.Claerbout,Fundamentals of GeophysicalData Processing,p.246-56,其以引用方式结合在此。可选地,速度分布可以从其他源独立的得到,例如,测井记录(well log)。根据估计的速度分布,数据处理装置将地震轨迹“偏移”,即将它们从时间的函数转化为深度的函数。
现有多种偏移技术,包括那些基于单程波动方程(单程WEM)的偏移技术,以及基于全程波动方程的逆时偏移(RTM)。单程WEM是更普遍、广泛使用的技术,因为在很多情况下是有效的并且具有相对较小的计算复杂度。然而,在具有复杂结构的区域,尤其是那些生成了强倒转波(例如棱波(prism waves))以及多个反射(可能产生双波(duplex waves))的区域,单程WEM无法简单地使复杂结构成像。这一情况在海域地震勘测中将会加重,因为窄方位角接收器的使用限制了交叉线偏移的数量。
RTM可以解决上述成像问题(参见,例如,E.Baysal,D.D.Kosloff,和J.W.C.Sherwood,“Reverse time migration”,Geophysics,48,1514[1983];G.A.McMechan,“Migration by extrapolation of time-dependentboundary values”,Geophysical Prospecting,31,413-420[1983];和N.D.Whitmore,“Iterative depth imaging by backward time propagation”,SEG Expanded Abstracts,2,382-385[1983].)。近年来,RTM在复杂媒介以及亚盐构造的叠前深度成像处理上更受青睐。然而,与单程WEM方法相比,RTM计算开销更大并且需要数据处理装置具有大的内存与大硬盘容量的计算机。而且,当将波场的高频分量由于有限差分模式的数值弥散(numerical dispersion)而偏移时,RTM就变得更加具有挑战性了。
图4为图1中示例地震模型的偏移深度图像402。利用有限差分前向模型模拟了100个炮点。炮点位置范围从6km到14km,并且具有80米的炮点间隔。每个炮点的接收器孔径覆盖了整个模型。接收器间隔为10米。图4的深度图像是利用全孔径来偏移所有炮点得到的。总之,弯曲事件(curved event)成像较好,但陡沉事件(steep dip event)404的振幅较弱。为了使该目标事件更好地成像,没有必要获取大量的新勘测数据或使用RTM再处理所有的现有数据。相反,我们仅关注那些对目标事件成像贡献较大的地震轨迹。下面描述的地震能见度分析方法和***将对提供这些贡献的轨迹进行定量确定。
图5a是地震能量传播的射线追踪的示意图。射线502显示了从源到目标事件一端并返回表面的能量传播,而射线504显示了从源到目标事件另一端并返回表面的能量传播。实际上,波动方程偏移用来向下传播地震能量并返回表面。更具体的,优选地使用基于波动方程的传播函数而不是基于高频渐近射线的方法来进行能见度分析。波动方程方法更精确并提供了更可靠的结果。
能见度分析有两个阶段。首先,模拟炮点的波场向下传播并且用软件测量来自目标事件的反射波场。在第二个阶段中,反射波场向回传播并且用软件测量目标对每个接收器记录的信号的贡献。通过将测量的贡献信号csr(t)的平方进行积分而获得目标事件的源-接收器能见度V(s,r)(近似于将图3轨迹中每个样本值平方以及求和):
其中,r是接收器位置,s是源位置,并且csr(t)是测量的贡献信号,其为时间的函数,时间位于t=0的炮点点火时刻和t=T的记录间隔结束之间。模拟炮点和接收器位置可以在模型区域均匀分布,或者在整个勘测环境(例如海上汽船几何分布)中由用户定制。
接收器能见度VR(r)定义为针对给定源位置s=S的源-接收器能见度V(s,r):
VR(r)=V(S,r)
图5b为针对图5a中所示的源点火位置的目标事件404的接收器能见度506。从图中可以看出对于该源位置,接收器能见度被极大的限定在了目标事件上方清晰可辨的邻域内。
此处词语源能见度VS(s)定义为在所有接收器位置{R}的源-接收器能见度V(s,r)的和:
图6为在三种不同假设下,目标事件404的源能见度。图6a为具有全接收器孔径(例如,所有接收器对所有源点火进行响应)的源能见度。图6b为假设±2km孔径(例如,只有源位置的2km内的接收器对源点火进行响应)的源能见度。图6c为零偏移源能见度(例如,接收器与源位于同一位置)。(应注意可通过“***反射器”假设来简化零偏移能见度计算,其中目标事件看作是点***的分布集,每个点***的强度与目标事件反射率相等)。在每种情况中,能见度随着源位置的不同而变化,并且在位于目标事件上方以及右侧的邻域中的炮点对目标事件的能见度贡献最大。这些源能见度函数之间的对比使得易于量化接收器孔径的作用。如果只考虑到零偏移情况,将损失大部分的能见度。
图7为目标事件404的源-接收器能见度函数V(s,r)的图示。水平轴表示接收器位置,而垂直轴表示源位置。在区域702得到的源-接收器位置的轨迹对事件404的能见度贡献最大。沿着线704的能见度函数是针对图5a中所示的源位置的接收器能见度函数VR(r)。图6c中显示的零偏移源能见度函数是沿着45°线708建立的。图6c中显示的源能见度是沿着整个图水平积分得到的。图6b中显示的源能见度函数是在线706之间经水平积分得到。
最大接收器孔径(在地震勘测中对应于电缆长度)可以通过调整线706的间距来捕获能见度函数下非零区域的大部分。然后通过选择源位置来捕获能见度函数下非零区域的大部分。使用这种策略来选择用于叠前深度偏移的轨迹(必要的话采集数据)能极大地减少为提高目标事件成像质量所需的工作量。图8显示了仅使用所选数据的深度偏移图像,其说明了能见度控制的图像的振幅行为更加平衡。
图9为示例的能见度分析方法,该方法可通过计算机***自动实现或在储层分析人员等用户的引导下实现。由模块902开始,***获得用来表示深度偏移图像的数据。在很多情况下,这些数据将由***本身基于以前的地震勘测自动生成。在模块904中,***确定一个或更多个目标事件,例如,未充分成像的特征。在一些实现中,***通过向用户显示深度偏移图像来确定目标事件,并且请求用户输入显示为充分定义或未充分定义的区域。
在模块906中,***选择了一种比生成原始数据偏移更复杂的偏移方法。例如,原始偏移可采用单程WEM,但是***也可以使用全波动RTM。当可以采用多重增强时,用户能够选择期望的偏移方法。
在模块908中,***使用所选的偏移方法来确定源-接收器能见度V(s,r),以在原始深度偏移数据确定的试验结构中模拟炮点。正如前面所提到的,通过计算来自多个源位置的每个的目标事件的反射波场,然后测量从这些反射波场到在多个接收器位置的每个测得的信号的贡献,来确定能见度。
在模块910中,***确定了那些现有的轨迹,这些现有的轨迹的源-接收器位置具有高于给定阈值的目标事件能见度。可以基于能见度值的峰值来预定所述阈值,或者选择所述阈值以捕获在多维能见度表面下方区域的预定比例(例如90%)。在模块912中,***将所选的偏移方法应用到所确定的高能见度轨迹。因为期望所确定的轨迹能够表示可用数据的小的子集,所以使用较复杂偏移方法更加可行。
在模块914中,***确定目标事件是否已充分成像,如果可以,方法跳转到模块922。在一些实现中,***通过向用户显示深度-偏移数据并请求用户反馈来进行该确定。如果目标事件仍没有充分成像,那么就需要采集额外的数据。所以,在模块916中,***至少部分基于源-接收器能见度的计算来确定勘测区域和其他勘测参数。在一些实施方式中,期望的源和接收器位置的区域能够通过画矩形(用于陆地勘测)或平行四边形(用于接收器位置随源位置变化的海域勘测)来确定,上述矩形或平行四边形基本涵盖了高能见度区域的大部分。
在模块918中,***从新勘测中得到了轨迹数据,并且在模块920中所选偏移方法被用于生成包含目标事件的区域的新的深度偏移数据图像。在模块922中,合成并显示一组合图像。该组合图像包括从原始偏移数据确定的整体结构,但也包括新偏移数据中的目标事件图像。然后储层工程师可以评估针对分析所充分确定的感兴趣结构的可能产量。
图10为示例的计算机50形式的能见度分析***,其具有能够与用户交互的一个或多个输入设备54以及一个或多个输出设备56。软件(如示例的便携式信息存储盘52)配置计算机处理器以接收用户指令并且相应的从网络或内部存储中获取数据,处理数据,并且生成图像供用户观察或分析。当实施所公开的方法时,软件通常将数据处理的负荷分布到通过网络互连的多个计算机。
上述描述是基于用于勘探目的的2D地震模型。在实践中,期望是3D成像,并且每个源和接收器位置由至少两个空间坐标确定。结果,源-接收器能见度图(见图7)应含有至少四个空间维度。然而,基本原理是相同的。
综上所述,本发明公开了多个地震能见度分析方法和***。这些***和方法针对复杂的媒介,在表面定量的确定了目标事件。对于给定的源-接收器几何分布,能见度强度指出了该几何分布下目标事件是否可见。这一知识可以应用于采集勘测设计和叠前深度偏移。能见度试验提供了下列突破:
·对于给定的目标事件,高能见度区域外的勘测数据是不必要的。
·对于给定的目标事件,对高能见度区域外的轨迹进行再偏移是不必要的。
·如果给定的目标事件没有高能见度区域,利用所选偏移方法不能重构。
当充分理解上述公开内容后对于本领域技术人员而言多种变形以及修改是显而易见的。下述权利要求目的在于包括上述所有的变形和修改。
Claims (17)
1.一种地震勘测方法包括:
确定目标事件的能见度,所述目标事件的能见度根据地震源位置和接收器位置而变化;以及
在所选择的区域中采集地震数据,该所选择的区域应至少部分包含能见度值高于阈值的位置,
其中所述确定目标事件的能见度的步骤包括:
针对多个模拟炮点的每一个,在地震模型中找到来自目标事件的反射波场;
针对多个接收器位置的每一个,计算每个反射波场的贡献信号;以及
将每个贡献信号转化为源-接收器能见度值,以及
其中所述阈值被选择为捕获在多维能见度函数下方的区域的预定比例。
2.如权利要求1所述的方法,进一步包括:
基于所采集的地震数据将目标事件成像。
3.如权利要求1所述的方法,其中至少部分的基于地震勘测的偏移数据来确定所述地震模型。
4.如权利要求1所述的方法,其中所述目标事件是由用户选择的未充分成像的地下特征。
5.如权利要求1所述的方法,其中所述能见度值为所述贡献信号平方的积分。
6.一种地震偏移方法包括:
在现有地震勘测的轨迹中确定源位置和接收器位置处的目标事件的能见度;以及
将能见度值高于阈值的轨迹进行再偏移以对目标事件成像,
其中确定目标事件的能见度的方法包括:
针对多个模拟炮点的每一个,在地震模型中找到所述目标事件的反射波场;
针对多个接收器位置的每一个,计算来自每个反射波场的贡献信号;以及
将每个贡献信号转化为源-接收器能见度值,以及
其中所述阈值被选择为捕捉在能见度函数下方的区域的预定比例。
7.如权利要求6所述的方法,其中至少部分的基于现有地震勘测的偏移数据来确定所述地震模型,其中通过第一个偏移方法来获得所述偏移数据。
8.如权利要求7所述的方法,其中所述转化包括确定贡献信号平方的积分。
9.如权利要求7所述的方法,其中使用不同于第一次偏移方法的第二次偏移方法确定所述能见度。
10.如权利要求7所述的方法,其中使用不同于第一次偏移方法的第二次偏移方法来进行所述再偏移。
11.如权利要求6所述的方法,其中所述目标事件是由用户选择的未充分成像的地下特征。
12.一种能见度分析***,包括:
存储能见度分析软件的存储器;以及
与所述存储器连接的至少一个处理器,以执行所述能见度分析软件,其中软件将所述至少一个处理器配置为:
针对多个模拟炮点的每一个,确定地震模型中至少一个目标事件的反射波场;
针对多个接收器位置的每一个,计算每个反射波场的贡献信号;以及
将每个贡献信号转化为源-接收器能见度值,
其中所述阈值被选择为捕捉在多维能见度函数下方的区域的预定比例。
13.如权利要求12所述的***,其中所述软件进一步配置所述处理器以图形化所述能见度,所述能见度根据源位置和接收器位置中的至少一个而变化。
14.如权利要求12所述的***,其中至少部分的基于地震勘测的偏移数据来确定所述地震模型,其中所述偏移数据通过第一次偏移方法获得。
15.如权利要求14所述的***,其中使用不同于第一次偏移方法的第二次偏移方法确定所述能见度。
16.如权利要求12所述的***,其中所述软件进一步配置所述处理器以从现有地震勘测中确定轨迹,用于使用增强偏移方法进行再偏移,其中所述轨迹至少部分的基于它们的能见度值来确定。
17.如权利要求12所述的***,其中所述软件进一步配置所述处理器,以针对额外地震数据采集来确定源位置和接收器位置,其中所述位置至少部分的基于源-接收器能见度值来确定。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2009/031400 WO2010082938A1 (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102282481A CN102282481A (zh) | 2011-12-14 |
CN102282481B true CN102282481B (zh) | 2014-11-12 |
Family
ID=42340028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200980154752.0A Expired - Fee Related CN102282481B (zh) | 2009-01-19 | 2009-01-19 | 基于地震能见度分析的数据采集和叠前偏移 |
Country Status (8)
Country | Link |
---|---|
US (1) | US9329288B2 (zh) |
EP (1) | EP2376947B1 (zh) |
CN (1) | CN102282481B (zh) |
AU (1) | AU2009337134B2 (zh) |
BR (1) | BRPI0924056A2 (zh) |
CA (1) | CA2747146C (zh) |
MX (1) | MX2011007084A (zh) |
WO (1) | WO2010082938A1 (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2443436B (en) * | 2006-11-04 | 2010-07-21 | Westerngeco Seismic Holdings | Method of determing properties of the earth using turning wave data |
US8116168B1 (en) | 2008-06-18 | 2012-02-14 | Halliburton Energy Services, Inc. | Hybrid one-way and full-way wave equation migration |
CN102282481B (zh) | 2009-01-19 | 2014-11-12 | 兰德马克图形公司 | 基于地震能见度分析的数据采集和叠前偏移 |
WO2010120301A1 (en) | 2009-04-16 | 2010-10-21 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
US8830788B2 (en) | 2011-02-24 | 2014-09-09 | Landmark Graphics Corporation | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media |
US9164184B2 (en) | 2011-05-27 | 2015-10-20 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
WO2012166733A1 (en) | 2011-05-27 | 2012-12-06 | Conocophillips Company | Two-way wave equation targeted data selection for improved imaging of prospects among complex geologic structures |
EP2715404A4 (en) * | 2011-05-27 | 2015-01-28 | Conocophillips Co | ABOUT A BIDIRECTIONAL WAVE EQUALIZATION DATA SELECTION FOR THE SEISMIC COLLECTION OF COMPLEX GEOLOGICAL STRUCTURES |
WO2013048585A1 (en) * | 2011-09-28 | 2013-04-04 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for improved imaging of complex geologic structures |
US10061046B2 (en) * | 2012-02-06 | 2018-08-28 | Ion Geophysical Corporation | Integrated passive and active seismic surveying using multiple arrays |
GB2510873A (en) | 2013-02-15 | 2014-08-20 | Total Sa | Method of modelling a subsurface volume |
GB2510872A (en) * | 2013-02-15 | 2014-08-20 | Total Sa | Method of modelling a subsurface volume |
US20150226868A1 (en) * | 2014-02-10 | 2015-08-13 | Microseismic, Inc. | Method for microseismic event moment magnitude estimation |
CN104570055A (zh) * | 2015-02-04 | 2015-04-29 | 中国海洋石油总公司 | 一种***震源控制*** |
EP3371749A1 (en) | 2015-11-06 | 2018-09-12 | Google LLC | Regularizing machine learning models |
FR3059700A3 (fr) * | 2016-12-02 | 2018-06-08 | Al Khatib Habib Agissant Au Nom Et Pour Le Compte De La Soc Spotlight En Cours De Formation | Procede d’amelioration des systemes actifs ultralegers de detection sismique |
US11137509B2 (en) * | 2018-02-02 | 2021-10-05 | Fairfield Industries, Inc. | Seismic imaging with a temporal decomposition imaging condition |
CN112782754A (zh) * | 2019-11-06 | 2021-05-11 | 中国石油天然气集团有限公司 | 多源同步采集数据模拟方法及装置 |
US11474267B2 (en) | 2020-06-11 | 2022-10-18 | China Petroleum & Chemical Corporation | Computer-implemented method and system employing compress-sensing model for migrating seismic-over-land cross-spreads |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5784334A (en) * | 1996-03-13 | 1998-07-21 | Atlantic Richfield Company | Method and system for detecting hydrocarbon reservoirs using amplitude versus offset analysis of seismic signals |
US6311131B1 (en) * | 1996-05-08 | 2001-10-30 | Lloyd Peardon | Method of and apparatus for determining the quality of seismic data |
US20020033832A1 (en) * | 2000-09-18 | 2002-03-21 | Rafail Glatman | Method for computer modeling of visual images and wave propagation |
US20050207278A1 (en) * | 2002-11-12 | 2005-09-22 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
US20050270537A1 (en) * | 2004-06-08 | 2005-12-08 | Mian Zahid F | Image-based visibility measurement |
Family Cites Families (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5014230A (en) | 1986-06-05 | 1991-05-07 | Xox Corporation | Solid-modeling system using topology directed subdivision for determination of surface intersections |
US4890242A (en) | 1986-06-05 | 1989-12-26 | Xox Corporation | Solid-modeling system using topology directed subdivision for determination of surface intersections |
US4953142A (en) * | 1989-01-06 | 1990-08-28 | Marathon Oil Company | Model-based depth processing of seismic data |
US5138584A (en) | 1989-09-06 | 1992-08-11 | Chevron Research & Technology Company | Migration of seismic turning waves |
US5274605A (en) | 1992-06-26 | 1993-12-28 | Chevron Research And Technology Company | Depth migration method using Gaussian beams |
US5490120A (en) | 1993-04-30 | 1996-02-06 | Union Oil Company Of California | Method of imaging overturned waves |
US5530679A (en) | 1993-05-10 | 1996-06-25 | Western Atlas International, Inc. | Method for migrating seismic data |
US5544126A (en) | 1993-11-04 | 1996-08-06 | Berryhill; John R. | System and method of seismic shot-record migration |
US6002642A (en) | 1994-10-19 | 1999-12-14 | Exxon Production Research Company | Seismic migration using offset checkshot data |
US6446007B1 (en) | 1998-11-05 | 2002-09-03 | Exxonmobil Upstream Research Company | Method for controlled-amplitude prestack time migration of seismic data |
US6021094A (en) | 1998-12-03 | 2000-02-01 | Sandia Corporation | Method of migrating seismic records |
US6687618B2 (en) | 2000-08-07 | 2004-02-03 | 3D Geo Development, Inc. | Typing picks to horizons in migration velocity analysis |
WO2002023222A1 (en) | 2000-09-15 | 2002-03-21 | Nutec Sciences, Inc. | Illumination weighted imaging condition for migrated seismic data |
US6611761B2 (en) | 2000-12-19 | 2003-08-26 | Schlumberger Technology Corporation | Sonic well logging for radial profiling |
US6864890B2 (en) | 2002-08-27 | 2005-03-08 | Comoco Phillips Company | Method of building and updating an anisotropic velocity model for depth imaging of seismic data |
US6778909B1 (en) | 2002-10-05 | 2004-08-17 | 3Dgeo Development, Inc. | Seismic data processing systems and methods |
US6819628B2 (en) | 2003-04-07 | 2004-11-16 | Paradigm Geophysical (Luxembourg) S.A.R.L. | Wave migration by a krylov space expansion of the square root exponent operator, for use in seismic imaging |
FR2858063B1 (fr) * | 2003-07-21 | 2005-09-16 | Geophysique Cie Gle | Procede d'estimation du taux de couverture d'illumination dans le domaine migre |
US7065004B2 (en) | 2004-04-22 | 2006-06-20 | Pgs Americas, Inc. | Horizon-based residual depth migration velocity analysis |
US7376517B2 (en) | 2005-05-13 | 2008-05-20 | Chevron U.S.A. Inc. | Method for estimation of interval seismic quality factor |
US20060270537A1 (en) * | 2005-05-27 | 2006-11-30 | Karol Fisher | Exercising system with trellis |
US20070162249A1 (en) | 2006-01-06 | 2007-07-12 | Min Lou | Traveltime calculation in three dimensional transversely isotropic (3D TTI) media by the fast marching method |
US7196969B1 (en) | 2006-02-09 | 2007-03-27 | Pgs Geophysical As | Three-dimensional two-way acoustic wave equation pre-stack imaging systems and methods |
US7952960B2 (en) | 2006-10-03 | 2011-05-31 | Bp Corporation North America Inc. | Seismic imaging with natural Green's functions derived from VSP data |
US8120991B2 (en) | 2006-11-03 | 2012-02-21 | Paradigm Geophysical (Luxembourg) S.A.R.L. | System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems |
US7400553B1 (en) * | 2006-11-30 | 2008-07-15 | Shengwen Jin | One-return wave equation migration |
US20080137480A1 (en) | 2006-12-07 | 2008-06-12 | Woodside Energy Limited | Method of Building a Subsurface Velocity Model |
US8194498B2 (en) | 2008-01-18 | 2012-06-05 | Westerngeco L.L.C. | Using a wave propagator for transversely isotropic media |
US7663972B2 (en) | 2008-02-22 | 2010-02-16 | Pgs Geophysical As | Method for three dimensional seismic travel time tomography in transversely isotropic media |
US8116168B1 (en) | 2008-06-18 | 2012-02-14 | Halliburton Energy Services, Inc. | Hybrid one-way and full-way wave equation migration |
US8296069B2 (en) | 2008-10-06 | 2012-10-23 | Bp Corporation North America Inc. | Pseudo-analytical method for the solution of wave equations |
CN102282481B (zh) | 2009-01-19 | 2014-11-12 | 兰德马克图形公司 | 基于地震能见度分析的数据采集和叠前偏移 |
WO2010120301A1 (en) | 2009-04-16 | 2010-10-21 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
US8332156B2 (en) | 2009-07-10 | 2012-12-11 | Chevron U.S.A. Inc. | Method for propagating pseudo acoustic quasi-P waves in anisotropic media |
US8406081B2 (en) | 2009-09-25 | 2013-03-26 | Landmark Graphics Corporation | Seismic imaging systems and methods employing tomographic migration-velocity analysis using common angle image gathers |
AU2009354764B2 (en) | 2009-11-02 | 2014-03-20 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a 3D reverse time migration with tilted transverse isotropy |
US8830788B2 (en) | 2011-02-24 | 2014-09-09 | Landmark Graphics Corporation | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media |
-
2009
- 2009-01-19 CN CN200980154752.0A patent/CN102282481B/zh not_active Expired - Fee Related
- 2009-01-19 EP EP09838520.6A patent/EP2376947B1/en active Active
- 2009-01-19 AU AU2009337134A patent/AU2009337134B2/en not_active Ceased
- 2009-01-19 MX MX2011007084A patent/MX2011007084A/es active IP Right Grant
- 2009-01-19 US US13/142,478 patent/US9329288B2/en active Active
- 2009-01-19 BR BRPI0924056-0A patent/BRPI0924056A2/pt not_active IP Right Cessation
- 2009-01-19 CA CA2747146A patent/CA2747146C/en active Active
- 2009-01-19 WO PCT/US2009/031400 patent/WO2010082938A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5784334A (en) * | 1996-03-13 | 1998-07-21 | Atlantic Richfield Company | Method and system for detecting hydrocarbon reservoirs using amplitude versus offset analysis of seismic signals |
US6311131B1 (en) * | 1996-05-08 | 2001-10-30 | Lloyd Peardon | Method of and apparatus for determining the quality of seismic data |
US20020033832A1 (en) * | 2000-09-18 | 2002-03-21 | Rafail Glatman | Method for computer modeling of visual images and wave propagation |
US20050207278A1 (en) * | 2002-11-12 | 2005-09-22 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
US20050270537A1 (en) * | 2004-06-08 | 2005-12-08 | Mian Zahid F | Image-based visibility measurement |
Also Published As
Publication number | Publication date |
---|---|
EP2376947A1 (en) | 2011-10-19 |
EP2376947A4 (en) | 2015-09-02 |
US20110273959A1 (en) | 2011-11-10 |
BRPI0924056A2 (pt) | 2015-07-14 |
AU2009337134A1 (en) | 2011-07-07 |
CA2747146A1 (en) | 2010-07-22 |
US9329288B2 (en) | 2016-05-03 |
EP2376947B1 (en) | 2021-08-25 |
WO2010082938A1 (en) | 2010-07-22 |
CA2747146C (en) | 2016-10-04 |
AU2009337134B2 (en) | 2013-11-21 |
CN102282481A (zh) | 2011-12-14 |
MX2011007084A (es) | 2011-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102282481B (zh) | 基于地震能见度分析的数据采集和叠前偏移 | |
Foti et al. | Guidelines for the good practice of surface wave analysis: a product of the InterPACIFIC project | |
CN102395902B (zh) | 使用快速面向目标照明计算的地震成像***及方法 | |
US9766358B2 (en) | System and method for local attribute matching in seismic processing | |
US9874646B2 (en) | Seismic data processing | |
US6763305B2 (en) | Subsurface illumination, a hybrid wave equation-ray-tracing method | |
US7952960B2 (en) | Seismic imaging with natural Green's functions derived from VSP data | |
US8082106B2 (en) | 3D surface related multiple elimination for wide azimuth seismic data | |
US8659974B2 (en) | System and method of 3D salt flank VSP imaging with transmitted waves | |
EA026344B1 (ru) | Система и способ получения и обработки сейсмических данных о полях упругих волн | |
CN113805237B (zh) | 偏移陆地交叉排列地震的使用压缩感测模型的方法和*** | |
CN111352151B (zh) | 一种井中虚源法 | |
US9658354B2 (en) | Seismic imaging systems and methods employing correlation-based stacking | |
US20140165694A1 (en) | Methods and systems for quality control of seismic illumination maps | |
Bosch et al. | Full-waveform inversion of intensity-focused seismic data | |
US12050294B1 (en) | Method and apparatus for building envelope guided low frequency model for seismic exploration of subsurface formations | |
McDonald et al. | Intelligent Monitoring Systems and Advanced Well Integrity and Mitigation | |
Ganssle | Seismic explorations | |
Kumar | Seismic data processing: velocity analysis & migration |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141112 Termination date: 20170119 |