CN1186647C - 处理和探测地震信号的方法及装置 - Google Patents
处理和探测地震信号的方法及装置 Download PDFInfo
- Publication number
- CN1186647C CN1186647C CNB971903492A CN97190349A CN1186647C CN 1186647 C CN1186647 C CN 1186647C CN B971903492 A CNB971903492 A CN B971903492A CN 97190349 A CN97190349 A CN 97190349A CN 1186647 C CN1186647 C CN 1186647C
- Authority
- CN
- China
- Prior art keywords
- mentioned
- eigenvalue
- covariance matrix
- vector
- cube
- 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
- 238000000034 method Methods 0.000 title claims abstract description 85
- 238000012545 processing Methods 0.000 title claims description 9
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 7
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 53
- 238000003325 tomography Methods 0.000 claims description 18
- 238000005259 measurement Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 239000004215 Carbon black (E152) Substances 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 3
- 230000008021 deposition Effects 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 claims description 2
- 238000004141 dimensional analysis Methods 0.000 claims 4
- 238000010586 diagram Methods 0.000 description 11
- 230000014509 gene expression Effects 0.000 description 9
- 230000008901 benefit Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 239000000523 sample Substances 0.000 description 4
- 230000035508 accumulation Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000003973 irrigation Methods 0.000 description 2
- 230000002262 irrigation Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000003129 oil well Substances 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- SGTNSNPWRIOYBX-UHFFFAOYSA-N 2-(3,4-dimethoxyphenyl)-5-{[2-(3,4-dimethoxyphenyl)ethyl](methyl)amino}-2-(propan-2-yl)pentanenitrile Chemical compound C1=C(OC)C(OC)=CC=C1CCN(C)CCCC(C#N)(C(C)C)C1=CC=C(OC)C(OC)=C1 SGTNSNPWRIOYBX-UHFFFAOYSA-N 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 208000035346 Margins of Excision Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/40—Transforming data representation
- G01V2210/48—Other transforms
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)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
一种用于探测碳氢化合物的方法及装置,包括步骤:获取分布在地球预定的三维体积上的地震信号轨迹分布组;将三维体积分割成多个具有位于此的至少两个地震轨迹部分的分析单元;计算每个单元内的地震轨迹向量;从这些向量形成对于每个单元的协方差矩阵;计算每个单元的协方差矩阵的主本征值和本征值之和并且由每个单元的协方差矩阵的主本征值与本征值之和的比率计算地震分布;和形成单元选择的组群的地震分布图。
Description
本发明总体上涉及一种地震信号探测方面的内容,尤其涉及一种通过在三维空间中识辨地层的结构和特点来对石油和气体进行探测和开发的方法和装置。
在地震探测中,沿着由陆地检波器组合或水上射光横向地震检波器组成的按线获取地震数据。地震检波器和水中地震检波器作为接收传递到地面并从地表下岩石的介面反射到表面的能量的传感器。能量通常由通过以一预定的间隔和频率摇动大地来把脉冲传递到表面的Vibroseis传运器提供到岸上。在水上的情况,通常采用气枪源。返回表面的能量的细微改变通常反映储库的地层、结构和流量。
在进行三维(3-D)地震探测中,原理是类似的,但是按线和阵列排的越密,就可提供越详细的地下覆盖面。对于这种高密度的覆盖,在最后的数据结果分析之前需要记录、储存和处理大量的数字数据。处理需要广泛的计算机资源和复杂的软件以加强从地下接收到的信号并抑制掩遮信号的伴随噪声。
一旦数据被处理,地球物理工作者以能够有效地表示地下结构显示的3-D数据立方(见图1)的形式汇编并分析地震信息。利用计算机工作站,解译员还可以将地层切片以不同的地震层研究库区问题。利用地震或信息源的数据还可以在任何方向进行竖直切片或分区。可把时间图转换成深度,以在特定的水平上提供结构的分析。
传统上为了用图象表示地震回波而去获取并处理地震数据。但是,由于在截面上显示的地层特点信息量有限,地层上的变化通常很难用传统的地震显示来探测。虽然这些观点提供了一个观察这些结构的绝大部分的机会,但要在没有记录到断层回波的3-D空间内识辨断层面是很困难的。
相关和相似(一种多波道相关测量)是两个地震轨迹相似性或相异性测量。当两个地震轨迹的相关性增加,它们就越相似。在一个由0到1的尺度上指定一种相关测量时,“0”表示最大的缺少相似性,而值“1”表示全部或完全的相似性(即两个理想轨迹)。多于两个轨迹相关性可类似定义。
计算相关性的一种方法是由Bahorich和Farmer(Amoco公司指定)在以序号08/353,934和1994年12月12日为申请日的美国专利申请中公开。一种用于计算相似的方法是由Marfurt等(Amoco公司指定)在以申请号60/005,032和1995年10月6日为申请日的美国专利申请中公开。Marfurt等的发明包括一个在后选倾角和方位角的强力搜索。
两种方法的优点已得到证明,但它们同样具有一些局限性。改进结果和计算速度总是期望的。
根据本发明,公开一种更为健全的多轨迹本征分解方法,它具有比原有已知方法较高的分辨率。在本发明的一个实施例中,公开了一种用于探测气体和石油的方法。方法包括下列步骤:搜集供给预定的地球三维体积的地震信号轨迹数据集;在运行的窗口决定的从至少两个地震轨迹形成的至少两个数据矢量的向量;通过增加向量形成协方差矩阵;计算作为协方差矩阵的至少主本征函数的地震特征;和形成从至少一部分预定的地球三维体积计算的地震特征的图象。
在本发明的另一个实施例中,本发明方法在用于引导计算机操作到计算地震特征的计算机可读介质(如磁盘、磁带、CD-ROM,等等)上编码。在本发明的另一个实施例中,图象从上面描述方法准备并被用于气体和石油沉积的定位。
本技术尤其适合解释3-D地震体积内的断层面和探测细微的3-D地层特征。这是因为被断层线切割的地震轨迹在断层两边一般较之地震轨迹具有不同的地震特征。沿时间切面测量地震轨迹相似性(即相关或3-D连续性)揭示沿断层线的低相关的特征。这样的相关值可揭示对于传统地震截面不显著的临界地层下表面细节。同样沿一系列时间切面计算相关性,这些断层特征确定了断层面或表面。
从下面的本发明的详细描述、实施例的描述,从权利要求书和从附图所表明,本发明其它大量的优点和特征将变得更为显著。
图1表示从处理3-D地震数据中得到的信息图示;
图2是根据现有技术中处理的3-D地震数据的水平时间图示;
图3是两个相邻的地震轨迹图示;
图4到图8是根据本发明的一对地震轨迹的相关性曲线;
图9是表示运行窗口分析立方的图示;
图10A、10B和10C为地震轨迹群的图示;
图11A、11B和11C为两维分析窗口的图示;
图12是基本方法流程图;和
图13A、13B和13C是根据Bahorich等、Marfurt等发明和本发明相同水平时间切面图示。
当本发明可以有许多不同形式的实施例时,几个特定的本发明实施例在图中表示和将详细描述。可是应该明白,本发明的公开考虑了本发明原理的例证和不期望把本发明限制为所描述的特定实施例或算法。
在描述本发明方法之前,将描述操作的基本原理。
在给定的欲检验相关性的N个样本的时间窗口或深度窗口考虑两个轨迹t1和t2。在图3中表示轨迹和相关分析窗口的代表性曲线。第一轨迹t1包括时间序列(t11,t12,…,t1N)和第二轨迹t2包括时间序列(t21,t22,…,t2N)。在这两个时间序列中,第一个下标表示轨迹数(即,轨迹1或轨迹2),而第二个下标表示样本数。
通过在同样的二维笛卡尔坐标系中绘制对应于另一个的轨迹,可更好地理解本发明的文中得到的相关性的含义。从两个轨迹绘制等同的时间样本(即,点对(t11,t21),(t12,t22),∧,(t1N,t2N)),产生一个交叉的两个时间序列。令x轴表示第一轨迹t1和y轴表示第二轨迹t1,结果绘于图4中。这是通过这些点表示两个轨迹相关的点形成的图形。
这两个相关轨迹的一般形状是通过椭圆形表示的点集。这个椭圆形是广义的因为它不代表每一单个点但是描述了所有点的“全部”特征。这个椭圆的长轴和短轴将指向由对点的几何学决定的方向。两个轴的长度同样由几何学决定。表示这些点的典型代表和对应的椭圆由图5表示。
椭圆长轴和短轴的方向和大小通过长矢量沿长轴指向和短矢量沿短轴指向的两个比例化的矢量表示。这两个矢量大小对应于数据协方差矩阵两个本征值和规一化的矢量对应于本征矢量。由各自本征值比例化的本征矢量表示长轴和短轴的大小和方向。“基本分量”对应于与主要本征值相联系的本征矢量。
附图(图6到图8)期望给出对以前讨论中隐含的力学的本质理解。在这些图中,使用简单的振幅和相位变化建立了轨迹,以及观察到这些变化对相联的本征值和本征矢量的作用。
图6表明两个相同的轨迹退化为一条45度的线(即,具有短轴的椭圆退化为零)。“零长度”表示第二本征值为“零”和表示对应于主本征值的本征矢量与长轴准直。当值为1时,相关性最大。
具有相等的振幅和45度相位差的两个轨迹的情况在图7中表示。表明了相位转换如何把长度变为短轴和因此增加了第二本征值的大小。同样表示了通过各自本征值两个本征矢量被比例化。由于这些轨迹中的差异,相关性值降到小于1。
最后,在图8中容许相位和振幅变化。两个轨迹具有45度相位差和2∶1振幅比。得出的椭圆有一个非零的短轴(第二本征值非零),反映了相位差。另外,由于振幅差异椭圆和本征矢量被旋转。再则,相位和振幅产生相关性的下降。
前述的要点是渐进给出相关可作为本征值λ1和λ2,和本征向量v1和v2的函数。相关的函数表示为:
相关=f(λ1,λ2,V1V2) (1)
两个轨迹的过程可容易地扩展为任意数量的轨迹。从实际的观点出发,随着轨迹数的增加计算负担加重以及可得到的的计算资源受到限制。对于一个3-D地震数据集,分析可在空间和时间(或空间和深度)内的移动或运行分析窗口重复,结果为在移动窗口中央测量相关。输出结果是包含由原有的数据体积定义相关值的3-D数据集。
本方法的一个优点是不同方面数据在本征值和本征向量中被分配,这将从下面的讨论变得更加清晰。信息如振幅和位相可用健全和严密的方式被分析和处理。例如分辨率可通过控制本征值和本征向量而得到改进。其好处可在计算相关切面中看到。根据本发明的相关值的另一个重要方面是它们对于原有数据不连续表现敏感和揭示细微的地质特征,如断层和沟渠。
回到本发明方法,第一步骤是得到或取得一组三维的地震数据。这样的数据是以分布在地球三维体积上的地震信号轨迹形式。处理诸如3-D地震数据的获取和降低为数字形式的方法是在本领域技术人员熟知的。这样的数据可通过在土地测量或海洋调查专业的地球物理方面的商人获得。这样的数据还通过商人出售或得到商人许可以及通常通过磁带传递或存储以输入地震工作站的存储中。
下一步骤是把3-D数据集分解为许多的单元或分析立方20(见图9)。这些立方20起着分解或把地震数据存储进组群或单元中以便进一步处理。结果是,一个分析立方总括了整个地震数据集或3D数据立方30。每一分析立方20包括一个一般的平板矩形时间层22的堆积。
为了简化之目的,3-D数据立方30、分析立方20和时间层22以直角的平行六面体或立方形式表示(一般为一个“单元”)。本领域中的这些技术在重复操作3-D数据集时倾向于矩形几何的简化。其它几何学和单元形状是可能的并在某些情况下是可以得到保证的。
回到图10,每一时间层22或堆积具有通过的地震轨迹t1(只给出了一个以避免图形混乱)的一部分。以3×3等间距排列的9个轨迹,在图9中给出。也可使用以星状的5个轨迹(见图10A)或三个轨迹(见图10B)。最好是对称排列。每个时间层中9个轨迹通常比3个轨迹要好。至少要用2个轨迹。
为帮助图解本发明概念,读者可参考图11A和图11B。尤其是,在9维(如9个轨迹情况)中图解通常很困难。在图11A和图11B中,两维的运行窗口24(或分析窗口)只给出包含的两个轨迹t1和t2。每一轨迹(见图11C)包括一个N个样本的时间序列
t1=(t11,t12,…,t1N)
t2=(t21,t22,…,t2N)
图11C的分析窗口将被分成多层垂直堆积的矩形时间层22。
应当明白,在处理选择窗口尺度和距离之前,要在分辨率和稳定性之间权衡。换言之,在参数估计中小的分析窗口或立方容许高的空间或时间频率,但给出小的统计稳定性或对于这些估计较少的自由度。另一方面,非常大的窗口具有低分辨率和趋于模糊数据以致于重要的地质特征被丢失。
回到图11C,每个时间层22内的数据点定义地震数据矢量(此处为1×2矩阵,N=2)。如此,图11C的两维时间层(或图9的三维时间层)形成或定义矢量窗口。因此,9轨迹(N=9)将得出具有9个元素的数据矢量。
在每个矢量窗口或时间层22中的数据矢量的外乘积产生了一个N×N矩阵。因此9元素数据矢量给出了9×9矩阵。相加这些矩阵(一个矩阵对每个矢量窗口)产生一个对于整个分析窗口24(或分析立方20)的N×N协方差矩阵。因此如果分析立方包括9个时间层,则9个N×N矩阵加在一起形成了一个N×N协方差矩阵。
根据本发明,通过计算协方差矩阵本征值可得到一个非常有用的和非显而易见的轨迹相关性方法。尤其是,计算了协方差矩阵最大或主要本征值和协方差矩阵本征值之和。这两种数字的比率表示了与本征值之和相关的主本征值的大小。同样还表明在分析立方中轨迹的变率。一个有用的地震分布可用数学式子来表示:
其中λ1是协方差矩阵本征值,λ1是主本征值。因此,ψ是为地震数据矢量的点元素正确程度的指标(见图4和5)。
协方差矩阵的每个和所有本征值不需要计算表示。本领域的技术人员知道已有的计算主本征值的方法(例如乘方法、Rayleigh系数、(两个中较快者)、等等)。而且,还知道协方差矩阵对角线元素之和等于协方差矩阵本征值之和。为方便起见,方程(2)的比率(即,相关性)可被指定为分析立方20或分析窗口的中点。应该理解,使用主本征值的一个优点是其趋于更直接表示分析窗口内的轨迹变率。主本征矢量将不是变率(即,相关)的度量。
在进行相关性的度量时,分析立方20或分析窗口22有效地扫描(即,横向和垂直向)图9所示的整个3-D体积或图11A的所有轨迹。最好是相邻的分析立方20或矢量窗口22相互叠置(见图11A的窗口24、24′和24″)。叠置改进了空间分辨率。
扫描整个3D数据体积30的分析窗口20和方程(2)相关性程度指定的结果是在每个数据矢量的位置沿每个轨迹分配的相关值的排列。如此,3D数据体积被转换成一种3D“相关立方”。
包含在相关立方中的相关度量的数据是通过以地震分布图的形式显示相关数据而得到更方便的解释。这样的图常常采用沿通过相关立方表面的相关值显示的形式。两个例子分别是通过普通水平时间切面的平板,和通过由地震解码器选择的地震地层线的曲面。另一个例子是地质沉积时间的线性表示以便从相同的地质年代捕获特征。
当以灰色阴影(如,白色表示高相关性和黑色表示低相关性)或另一彩色标尺来显示,相关值容易被显示而被解释。(见美国专利USPatent4,979,699“用于彩色绘制地球物理数据的方法”)。
例如,Landmark和GeoQuest解释工作站通过载人作为地震体积的相关立方可用于显示和解释断层和地质学特征。这样的工作站通常被本领域的技术人员使用。未处理的3D地震数据通过磁带或磁盘方法可被方便地载人工作站,这些磁带或磁盘被指示计算机执行上述程序的指令所编码。可视软件(如,Landmark′s SeisCube软件)可通过相关立方快速切面以帮助理解复杂的断层关系。包括以地震分布图印刷输出的相关显示,当用于选择何地震线可解释时,可减少解释循环时间,使得解码器围绕弱数据区工作。另外,对于传统地震显示不明显的细微的地质学特征和复杂的断层,可很快得到证实和解释。
图13A、13B和13C提供了通过另外的处理和根据本发明的相同地震信息显示和处理的比较。差异是显著的。
当然,本发明方法通常很方便地通过编写计算机程序来执行和完成刚才描述的步骤。这样的处理被先前确认的工作站日常执行。
一个基本的处理方法流程图示于图12。在本发明的一个实施例中,用FORTRAN77编写的软件程序执行刚才描述的过程。3-D地震数据被写入存贮中。根据得到的地震数据的大小和内容,执行起始步骤32,程序参数初始化,数据区间建立,进行基本检查,和设定窗口大小。读出错误或读出用户选择路径。下一过程34开始。尤其是,子程序36被调入扫描具有一个分析立方的数据体积。在每一分析立方中,子程序38计算协方差矩阵,在另一个子程序40计算主本征值,本征值之和及综合的相关值。最后,结合步骤42和计算值的结果以相关立方的形式被存储在步骤44。因此,工作站操作员可得到相关立方用于在CRT46显示选择点(例如,通过立方的时间切面),以产生印刷输出或地震分布图48,开始进一步分析或传递给存储或磁带50以便在其它地方进一步处理。
由于要考虑到被3-D测量覆盖的区域地层学和地质学的已知部分,本领域的技术人员要小心使用本发明的方法。
相关图在几个3-D测量中被运行。在合理的数据质量当中,大约90%的断层可被确认。在相关图上被识辨的断层在地震部分非常细微,但是清晰地展现在相关图中,因为方法较为健全和图象展现了断层图形。因为相关图可在未解码的时间切面运行,本发明提供了一种结构框架的最大加速绘图和揭示断层关系细节的方法,否则只有通过长时间的断层探测才被解码。
特例
相关图沿探测到的地层产生并清晰地识辨挤人海上的盐/页岩。
在其它地点利用本发明的方法清楚地标识火山泥浆和气体。
几个时间相关切面显著地表示了地层特征,如废弃的河道、泥石流、尖沙洲和海底峡谷。在地震剖面上这些特征有时明显,但在有些情况下即使近距观察也不能识辨。
本发明,类似于Bahorich等和Marfurt等的发明,提供了一种在3D体积中揭示断层平面方法,而在3D体积中没有记录断层反射。断层通常对于石油的积累非常重要。一个断层通过切断地形结构或地层结构可形成一个隔离层,使得石油再次被断层包围。另一方面,如果断层平面包含没有完全结成的石块,它将形成流体的沟渠。这将容许碳氢化合物沿断层平面漂移进入地形并被包围或通过沿断层平面漂移出去。因此,例如断层线可预告储库中的流动图象并在喷口和产油井之间联系。
地震的不连续性还可提供与油井之间储库预测的必然联系和建立储库连续性以及沿油田的流动图象。相关技术可被用于地表下结构和沉积学地形如断层、盐的挤入、不相称性、沟渠系统、岩溶和通常与捕获和储存碳氢化合物相联系的碳酸盐礁表面的寻找、识辨和绘图。因此,本技术有助于碳氢化合物的寻找、分离和生产。另外,它用于识辨浅和深的钻探的危险性(例如,气体离地面很近的地方或不稳定的地方)。另一个例子是应用本发明来寻找已知的储库或地下洞穴的泄漏路径。
3-D地震的相关绘图是用于结构和地层学绘图的功能极其强大和非常有效的工具。新方法对于小波特征的横向变化尤其敏感,并且对于小波中横向变化的普通原因(例如,断层位移或地层变化)尤其敏感。因此,本主题发明围绕一种处理方法,及其以计算机指令形式被记录的装置、该方法的产物(如,绘图)和用于气体和石油探测的产物的手段。
由前述,对于本领域的技术人员,可发现大量的变化、替换和改进是显而易见的。因此,本描述可只作为解释以及用于告诉本领域的技术人员完成本发明的方法。例如,地震轨迹被描绘成等间距。不等间距的轨迹(见图10C)通过插值可方便转化为等间距。如另一个例子,对于滤除形成数据矢量的轨迹以分隔远离部分是有用的。中等滤除可用于每一数据矢量元素的分类。切除边缘可达到进一步的平滑。
另外,其它算法可用于测量地震数据邻近区域的相似性或产生“相关立方”。相关值或方程(2)的地震分布用作在地质方程内信号不连续性以及沿断层截面和风化的不一致性的更为健全的估算或计算。本方法还建议协方差矩阵本征值的其它结合(如算术平均、平方根平均、平均、中值、平方之和的平方根、平方积的平方根、最小、最大、和、积,等等)。另外,本发明过程可与其它分布(如,AVO斜面)相结合并且应用于多向量地震数据。
本发明某种特征还可独立用于本发明的其它特征。例如,根据本发明识辨的地质特征可与速度图叠置以提供交叉检查速度的方法。
此外,相关切面图是功能非常强大的绘图工具,当与振幅的勘测绘图和倾角图结合使用时,在用易于获得的3-D地震数据对墨西哥湾或类似盆地绘制的效果有希望成为该技术里程碑。可以相信,通过以地图的角度绘制并通过较少的传统逐线择取将会加速结构和地层的详细绘制。以地图的角度解释“勘测”数据在解释的质量和数量上提供了显著的改进。另外,本发明的处理方法将随之加快速度。当这种勘测成为可能时,这种速度有助于迅速作出投标的选择。
最后,还应理解本发明的原理也可同样地用于其它领域(如声波定位仪,在这种情况中传感器是音响装置并且信号源是敌方的潜艇;地震和核武器引爆探测系统,在这种情况中,传感器是测震设备而信号源是地震或爆破震中;天文干涉仪,其中的传感器应是无线天文望远镜而信号源应是遥远的星系或类星体;相控阵雷达,其中的传感器应是阵列天线),其中信号(如雷达,声纳,无线频率能等)被处理形成图象或在由这些图象代表的结构中的位置变化。
因此,应该理解可以在不脱离如所附的权利要求限定的本发明实质和范围内作各种修改、变化。当然包含在权利要求范围内的所有的这些修改也被所附的权利要求所覆盖。
Claims (46)
1.一种用于探测碳氢化合物的方法,包括如下步骤:
a)得到分布在地球预定的三维体积的地震轨迹组;
b)将上述三维体积分割成多个的垂直堆积和一般间隔的水平时间层并将上述时间层排列进入为大量的横向和垂直向扩展的单元,具有至少两个地震轨迹部分的每一个上述时间层位于定义的数据矢量处;
c)在上述单元每一上述时间层内计算上述数据矢量的向量;
d)结合上述向量得到对于每一上述单元的协方差矩阵;
e)在每一上述单元内计算上述地震轨迹的相关性度量,其特征在于,上述相关性度量是至少为上述协方差矩阵的最大本征值的函数;和
f)从多个的上述地震轨迹的上述相关性度量形成地震分布图。
2.根据权利要求1的方法,其特征在于,在步骤(f)中上述图通过显示相对于经过预定的地震地层的表面的上述相关性度量来形成。
3.根据权利要求1的方法,其特征在于,在步骤(f)中上述图通过显示相对于经过预定的时间线的表面的上述相关性度量来形成。
4.根据权利要求1的方法,其特征在于,在步骤(b)中上述单元包括具有至少位于立方内的5个地震轨迹的分析立方;和在步骤(c)中每一向量是以5×5矩阵的形式。
5.根据权利要求1的方法,其特征在于,在步骤(b)中上述单元包括具有至少位于立方内的9个地震轨迹部分的分析立方;和上述数据向量有9个元素。
6.根据权利要求5的方法,其特征在于,在步骤(b)中上述9个地震轨迹以3×3格点排列。
7.根据权利要求1的方法,其特征在于,在步骤(b)中上述单元厚度至少小于100毫秒。
8.根据权利要求1的方法,其特征在于,在步骤(c)中以时间为主畴。
9.根据权利要求1的方法,其特征在于,在步骤(e)中执行包括:计算上述协方差矩阵的上述最大本征值,计算上述协方差矩阵的上述本征值之和,和计算上述协方差矩阵的上述最大本征值与上述本征值之和的比率。
10.根据权利要求9的方法,其特征在于,上述协方差矩阵的上述本征值之和是通过计算上述协方差矩阵的对角元素之和形成。
11..根据权利要求1的方法,其特征在于,在执行步骤(b)中在每一个单元中的上述两个地震轨迹的其中一个位于一个相邻的单元,因此上述单元相互叠置。
12.一种定位洞穴地形、断层和外形的方法,包括如下步骤:
a)获取覆盖一个预定的地球体积的3-D地震数据,上述数据包括具有时间、位置和振幅特征的地震轨迹;
b)将上述体积至少一部分分割成至少一个相对小的、相邻的、叠置的、三维的分析立方,其特征在于,每一上述分析立方包括至少三个横向分离的地震轨迹,和将每一个分析立方分割成多个的样本间隔,其中每一样本间隔定义多个1×3数据向量;
c)对作为上述数据矢量的向量形成的协方差矩阵主本征值函数的上述立方计算地震分布;和
d)将上述地震分布存储为用于显示的形式。
13.根据权利要求12的方法,其特征在于,在步骤(c)中上述地震分布是上述主本征值与上述立方及上述协方差矩阵至少两个上述本征值之和的比值的函数。
14.根据权利要求13的方法,其特征在于,在步骤(c)中上述地震分布是上述主本征值与上述立方及上述协方差矩阵所有对角元素之和的比值的函数。
15.根据权利要求14的方法,其特征在于,上述地震分布被指定给分析立方的中央。
16.根据权利要求15的方法,其特征在于,步骤(b)在多个的时间切面的执行进一步包括步骤:
e)显示经过上述分析立方及上述中央的连续时间切面的上述地震分布,识辨相对于空间和时间的不变特征。
17.一种地震探测方法,在地震探测中3-D地震数据包括反射的地震能量被作为时间的函数记录以产生一系列地震轨迹,和采用计算机处理这样的地震轨迹,该方法包括:
被计算机可读的介质和用于上述计算机执行过程的指令包括步骤:
(a)在预定的体积上获取3-D地震数据,上述数据包括具有时间、位置和振幅特征的地震信号矢量;和
(b)确定上述体积及上述3-D地震数据的邻近区域的相似性,通过:
(1)将上述数据至少一部分分割成相对小的、相邻的、叠置的、三维的分析立方,其中每一上述分析立方包括至少两个数据矢量;和
(2)对以上述立方上述矢量的向量之和形成的协方差矩阵基本本征值为函数的每一立方计算地震分布。
18.根据权利要求17的方法,其特征在于,上述介质载有使计算机执行步骤(2)来计算上述基本本征值与上述协方差矩阵本征值之和的比值的指令。
19.根据权利要求17的方法,其特征在于,上述介质载有使计算机执行步骤(2)来计算上述基本本征值与上述协方差矩阵对角元素之和的比值的指令。
20.根据权利要求19的处理方法,其特征在于,上述介质载有使计算机执行步骤(1)来形成具有位于此的至少5个地震轨迹的通常矩形排列的分析立方的指令;并且其中上述协方差矩阵为至少5×5矩阵和至少由3个向量矩阵形成。
21.根据权利要求20的处理方法,其特征在于,上述介质载有使计算机指定地震分布给其分析立方的中央的指令。
22.在地震探测中其中反射的地震能量被作为时间的函数记录以产生一系列地震轨迹,该方法包括步骤:
(a)确定由至少两个地震轨迹形成的两个向量数据矢量;
(b)通过增加步骤(a)的向量形成一个协方差矩阵;
(c)计算步骤(b)的以协方差矩阵的至少基本本征值为函数的地震分布;
(d)沿至少一个时间窗口的至少一部分重复步骤(a)到(b);
(e)在上述时间窗口形成上述地震分布的一个图。
23.根据权利要求22的方法,其特征在于,在步骤(c)通过计算上述基本本征值与至少上述协方差矩阵的本征值之和的一部分的比率来执行。
24.根据权利要求22的方法,其特征在于,在步骤(c)通过计算上述基本本征值与至少上述协方差矩阵的对角元素之和的一部分的比率来执行。
25.根据权利要求22的方法,其特征在于,步骤(d)通过从先前的步骤(a)操作使用至少一个地震轨迹和至少位于相邻于上述至少一个地震轨迹的两个新地震轨迹。
26.根据权利要求22的方法,其特征在于,步骤(a)包括步骤:
(1)在预定的地球体积上获取3-D地震数据,上述3-D地震数据具有时间、位置和振幅特征;和
(2)将上述体积一部分分割成包含相对小的、叠置的、包括上述至少两个地震轨迹的三维分析立方的排列的至少一个时间窗口。
27.一种地震探测方法,包括如下步骤:
a)读取包括在地球体积上地震信号轨迹分布的3-D地震数据组;
b)从上述体积选择至少一个地层切面并形成以横向扩展行和列排列的单元,每一个上述单元包括至少3个一般扩展的地震轨迹;
c)计算每一上述单元:
(1)数据矢量向量由多个的在上述单元中央每一边的时间间隔定义;
(2)由步骤(1)上述向量的协方差矩阵;
(3)上述协方差矩阵至少最大的本征值;和
d)沿上述至少一个地层切面检查上述单元上述本征值。
28.根据权利要求27的方法,其特征在于,步骤(3)是通过显示经过至少一个水平时间切面的上述单元的上述最大本征值图象而被执行。
29.根据权利要求28的方法,其特征在于,图象是上述单元上述协方差矩阵之和及上述单元上述最大本征值的函数。
30.一种地震探测方法,在地震探测中反射的地震能量被作为时间的函数记录以产生一系列地震轨迹,其方法包括步骤:
(a)将3-D地震数据安排在相对小的、叠置的、包括多个地震轨迹的三维分析立方内;
(b)确定由上述分析单元定义的数据矢量的向量;
(c)对每一立方通过增加步骤(b)的上述向量形成一个协方差矩阵;
(d)计算以协方差矩阵的所有本征值之和与每一协方差矩阵的基本本征值比率为函数的地震分布;和
(e)安排上述地震分布作为一个图来显示。
31.一种装置,包括:
用于一个过程的载有指令的计算机可读方法包括步骤:
(1)把覆盖地球预定体积的3-D地震数据读入存储;
(2)把上述3-D地震数据数字分类给相对小的三维单元排列,其中每一上述单元包含至少三个地震轨迹;
(3)从由多个的上述至少三个轨迹形成的协方差矩阵的本征值计算每一上述单元相关值。
(4)存储用于显示由上述相关值表示的与洞穴地形二维图形式相同的上述单元上述相关值。
32.根据权利要求31的装置,其特征在于,在步骤(3)中上述相关值至少是上述协方差矩阵的最大上述本征值的函数。
33.根据权利要求32的装置,其特征在于,上述相关值是上述本征值之和及上述最大本征值的函数。
34.根据权利要求32的装置,其特征在于,上述计算机可读方法是从包括磁盘、磁带、光盘和CD-ROM等组中选择的。
35.一种定位洞穴地形、断层和外形的方法,包括如下步骤:
a)获取覆盖一个预定的地球体积的地震数据;
b)将上述体积分割成相对小的三维单元排列,其中每一上述单元具有至少两个位于此的地震矢量;
c)从上述数据矢量的向量计算一个协方差矩阵;
d)绘制上述协方差矩阵的本征值图象。
36.根据权利要求35的方法,其特征在于,步骤(c)是通过使用增加多个向量形成协方差矩阵而被执行。
37.根据权利要求35的方法,其特征在于,步骤(d)是通过绘制本征值之和与最大本征值的比率而被执行。
38.一种预测碳氢化合物沉积的方法,包括如下步骤:
a)得到用于地球预定的三维体积的3-D地震数据相关值的地震分布图,上述图通过使用计算机得到并且一种用于引导上述计算机的计算机的程序执行下面步骤:
(1)读取上述数据并将上述体积分类到相对小的三维单元的排列中,其中每一上述单元具有至少两个位于此的地震数据矢量;和
(2)对以上述数据矢量的向量形成的协方差矩阵的本征值为函数的上述地震轨迹计算每一上述单元相关值。
b)使用上述图识辨地表下结构和通常与捕获和存储碳氢化合物相联系的沉积特征。
39.根据权利要求38的方法,进一步包括使用上述图以识辨钻探危险性。
40.根据权利要求39的方法,进一步包括在步骤(b)中识辨位置的钻探步骤。
41.根据权利要求38的方法,其特征在于,上述程序引导上述计算机执行步骤(a)(2),通过:
(i)计算每一协方差矩阵的最大本征值和上述协方差矩阵本征值之和;和
(ii)计算上述最大本征值与上述之和的比率。
42.根据权利要求41的方法,其特征在于,在步骤(i)中上述程序通过计算上述协方差矩阵对角元素之和引导上述计算机计算上述本征值之和。
43.一种被计算机工作站使用并用于石油和气体探测的装置,包括被计算机可读的介质和包含相关立方的图象,上述相关立方包括3D地震数据相关测量,每个上述测量为由两个地震数据矢量的至少两个向量相加形成的协方差矩阵本征值的函数。
44.根据权利要求43的装置,其特征在于,上述数据矢量具有空间和时间坐标特征;和其中相关的上述测量被指定给空间和时间坐标。
45.根据权利要求44的装置,其特征在于,每个上述测量为至少各自协方差矩阵主本征值的函数。
46.根据权利要求45的装置,其特征在于,每个上述测量为至少本征值之和的函数。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US63178896A | 1996-04-12 | 1996-04-12 | |
US08/631,788 | 1996-04-12 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1189218A CN1189218A (zh) | 1998-07-29 |
CN1186647C true CN1186647C (zh) | 2005-01-26 |
Family
ID=24532737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB971903492A Expired - Fee Related CN1186647C (zh) | 1996-04-12 | 1997-01-02 | 处理和探测地震信号的方法及装置 |
Country Status (9)
Country | Link |
---|---|
US (1) | US5892732A (zh) |
EP (1) | EP0832442B1 (zh) |
CN (1) | CN1186647C (zh) |
AR (1) | AR006341A1 (zh) |
AU (1) | AU710968B2 (zh) |
CA (1) | CA2220274C (zh) |
NO (1) | NO319052B1 (zh) |
RU (1) | RU2169931C2 (zh) |
WO (1) | WO1997039367A1 (zh) |
Families Citing this family (107)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6765570B1 (en) * | 1998-07-21 | 2004-07-20 | Magic Earth, Inc. | System and method for analyzing and imaging three-dimensional volume data sets using a three-dimensional sampling probe |
FR2786294B1 (fr) * | 1998-11-19 | 2000-12-22 | Elf Exploration Prod | Methode de detection de structures chaotiques dans un milieu donne |
DE19904347C2 (de) * | 1999-02-03 | 2002-08-14 | Henning Trappe | Verfahren zur seismischen Datenverarbeitung |
US6332037B1 (en) * | 1999-04-09 | 2001-12-18 | Board Of Regents Of The University Of Nebraska | Invariant, Eigenvalue based, non-degenerate data structure characterization, storage and retrieval indexing method |
US6625311B1 (en) | 1999-04-09 | 2003-09-23 | The Board Of Regents Of The University Of Nebraska | Methodology for data structure characterization, indexing, storage and retrieval |
US6611609B1 (en) * | 1999-04-09 | 2003-08-26 | The Board Of Regents Of The University Of Nebraska | Method of tracking changes in a multi-dimensional data structure |
DK1058131T3 (da) | 1999-06-03 | 2004-09-20 | Jason Geosystems B V | Fremgangsmåde til bestemmelse af rumlige ændringer i underjordisk struktur, stratigrafi, litologi og fluidum-indhold og reduktion af seismisk stöj |
US6594585B1 (en) | 1999-06-17 | 2003-07-15 | Bp Corporation North America, Inc. | Method of frequency domain seismic attribute generation |
DE19933717C1 (de) * | 1999-07-19 | 2001-01-11 | Henning Trappe | Verfahren zur seismischen Datenverarbeitung |
DE19943325C2 (de) * | 1999-09-10 | 2001-12-13 | Trappe Henning | Verfahren zur Bearbeitung seismischer Meßdaten mit einem neuronalen Netzwerk |
US6418381B1 (en) | 1999-12-06 | 2002-07-09 | Landmark Graphics Corp. | Transformation method applied to seismic data traces |
US6490526B2 (en) * | 2000-03-20 | 2002-12-03 | Exxonmobil Upstream Research Company | Method for characterization of multi-scale geometric attributes |
US6490528B2 (en) | 2000-04-17 | 2002-12-03 | Exxonmobil Upstream Research Company | Method for imaging discontinuites in seismic data |
WO2002003099A2 (en) | 2000-06-30 | 2002-01-10 | Exxonmobil Upstream Research Company | Method for imaging discontinuities in seismic data using dip-steering |
US6571177B1 (en) | 2000-09-18 | 2003-05-27 | Conoco Inc. | Color displays of multiple slices of 3-D seismic data |
US7006085B1 (en) * | 2000-10-30 | 2006-02-28 | Magic Earth, Inc. | System and method for analyzing and imaging three-dimensional volume data sets |
US6487502B1 (en) | 2000-12-01 | 2002-11-26 | Rdsp I, L.P. | System for estimating the locations of shaley subsurface formations |
US6597994B2 (en) * | 2000-12-22 | 2003-07-22 | Conoco Inc. | Seismic processing system and method to determine the edges of seismic data events |
US6850845B2 (en) * | 2001-07-20 | 2005-02-01 | Tracy Joseph Stark | System for multi-dimensional data analysis |
US6597992B2 (en) * | 2001-11-01 | 2003-07-22 | Soil And Topography Information, Llc | Soil and topography surveying |
US6766252B2 (en) | 2002-01-24 | 2004-07-20 | Halliburton Energy Services, Inc. | High resolution dispersion estimation in acoustic well logging |
AU2003223297B2 (en) * | 2002-03-27 | 2006-03-09 | Westerngeco Seismic Holdings Limited | Parametric fk techniques for seismic applications |
US20060122780A1 (en) * | 2002-11-09 | 2006-06-08 | Geoenergy, Inc | Method and apparatus for seismic feature extraction |
GB2403803B (en) * | 2003-07-05 | 2005-09-21 | Westerngeco Ltd | Processing seismic data |
US7298376B2 (en) * | 2003-07-28 | 2007-11-20 | Landmark Graphics Corporation | System and method for real-time co-rendering of multiple attributes |
CA2485761C (en) * | 2003-10-24 | 2015-11-24 | Bernd Milkereit | Resonance scattering seismic method |
US7796468B2 (en) * | 2004-02-26 | 2010-09-14 | Saudi Arabian Oil Company | Prediction of shallow drilling hazards using seismic refraction data |
US8234923B2 (en) | 2004-09-20 | 2012-08-07 | Innervision Medical Technologies Inc. | Systems and methods for ultrasound imaging |
US7554883B2 (en) * | 2004-10-11 | 2009-06-30 | Landmark Graphics Corporation | Fault filter for seismic discontinuity data |
GB2422012B (en) * | 2005-01-11 | 2008-09-10 | Statoil Asa | Method of seismic signal processing |
GB2429278B (en) * | 2005-08-15 | 2010-08-11 | Statoil Asa | Seismic exploration |
US7914451B2 (en) * | 2005-09-15 | 2011-03-29 | Innervision Medical Technologies Inc. | Determining attributes using ultrasound |
US7333392B2 (en) | 2005-09-19 | 2008-02-19 | Saudi Arabian Oil Company | Method for estimating and reconstructing seismic reflection signals |
WO2007092054A2 (en) | 2006-02-06 | 2007-08-16 | Specht Donald F | Method and apparatus to visualize the coronary arteries using ultrasound |
CN103185896B (zh) | 2006-09-01 | 2016-08-10 | 哈利伯顿兰德马克绘图公司 | 用于对波形体进行成像的系统和方法 |
WO2008051639A2 (en) | 2006-10-25 | 2008-05-02 | Maui Imaging, Inc. | Method and apparatus to produce ultrasonic images using multiple apertures |
GB2443843B (en) * | 2006-11-14 | 2011-05-25 | Statoil Asa | Seafloor-following streamer |
US8638328B2 (en) * | 2007-01-05 | 2014-01-28 | Landmark Graphics Corporation | Systems and methods for visualizing multiple volumetric data sets in real time |
US20080165185A1 (en) * | 2007-01-05 | 2008-07-10 | Landmark Graphics Corporation, A Halliburton Company | Systems and methods for selectively imaging objects in a display of multiple three-dimensional data-objects |
WO2008111037A2 (en) * | 2007-03-12 | 2008-09-18 | Geomage 2003 Ltd | A method for identifying and analyzing faults/fractures using reflected and diffracted waves |
EP2150841A1 (en) * | 2007-05-17 | 2010-02-10 | Spectraseis AG | Seismic attributes for reservoir localization |
US7724608B2 (en) * | 2007-07-20 | 2010-05-25 | Wayne Simon | Passive reflective imaging for visualizing subsurface structures in earth and water |
US9171391B2 (en) | 2007-07-27 | 2015-10-27 | Landmark Graphics Corporation | Systems and methods for imaging a volume-of-interest |
WO2009029135A1 (en) * | 2007-08-24 | 2009-03-05 | Exxonmobil Upstream Research Company | Method for predicting well reliability by computer simulation |
US8548782B2 (en) | 2007-08-24 | 2013-10-01 | Exxonmobil Upstream Research Company | Method for modeling deformation in subsurface strata |
US8768672B2 (en) * | 2007-08-24 | 2014-07-01 | ExxonMobil. Upstream Research Company | Method for predicting time-lapse seismic timeshifts by computer simulation |
CA2690991C (en) * | 2007-08-24 | 2013-12-24 | Exxonmobil Upstream Research Company | Method for multi-scale geomechanical model analysis by computer simulation |
US7630865B2 (en) * | 2007-09-11 | 2009-12-08 | Geomage (2003) Ltd | Complex analysis of kinematics for non-hyperbolic moveout corrections |
US9788813B2 (en) | 2010-10-13 | 2017-10-17 | Maui Imaging, Inc. | Multiple aperture probe internal apparatus and cable assemblies |
US9282945B2 (en) | 2009-04-14 | 2016-03-15 | Maui Imaging, Inc. | Calibration of ultrasound probes |
GB0722469D0 (en) | 2007-11-16 | 2007-12-27 | Statoil Asa | Forming a geological model |
US7702463B2 (en) | 2007-12-12 | 2010-04-20 | Landmark Graphics Corporation, A Halliburton Company | Systems and methods for enhancing a seismic data image |
GB0724847D0 (en) | 2007-12-20 | 2008-01-30 | Statoilhydro | Method of and apparatus for exploring a region below a surface of the earth |
GB0803701D0 (en) * | 2008-02-28 | 2008-04-09 | Statoilhydro Asa | Improved interferometric methods and apparatus for seismic exploration |
US8209126B2 (en) * | 2008-04-01 | 2012-06-26 | Geo{umlaut over (m)}age (2003) Ltd. | Wavefront-defined Radon transform |
EA026356B1 (ru) * | 2008-05-22 | 2017-03-31 | Эксонмобил Апстрим Рисерч Компани | Способ (варианты) и машиночитаемый носитель для сейсмического исследования глубинной зоны |
US8213261B2 (en) * | 2008-05-22 | 2012-07-03 | Exxonmobil Upstream Research Company | Method for geophysical and geological interpretation of seismic volumes in the domains of depth, time, and age |
CA2723381C (en) * | 2008-06-06 | 2017-02-07 | Landmark Graphics Corporation, A Halliburton Company | Systems and methods for imaging a three-dimensional volume of geometrically irregular grid data representing a grid volume |
US8805654B2 (en) * | 2008-06-09 | 2014-08-12 | Landmark Graphics Corporation | Distribution of properties in a 3D volumetric model using a maximum continuity field |
US8295124B2 (en) * | 2008-08-15 | 2012-10-23 | Bp Corporation North America Inc. | Method for separating independent simultaneous sources |
EA024624B1 (ru) * | 2008-11-14 | 2016-10-31 | Эксонмобил Апстрим Рисерч Компани | Способ (варианты) для обнаружения аномалий в наборах геофизических данных с использованием выполняемого в окне статистического анализа и способ добычи углеводородов из подземной области |
WO2010120907A2 (en) | 2009-04-14 | 2010-10-21 | Maui Imaging, Inc. | Multiple aperture ultrasound array alignment fixture |
US8463551B2 (en) * | 2009-11-17 | 2013-06-11 | Schlumberger Technology Corporation | Consistent dip estimation for seismic imaging |
EP2536339B1 (en) | 2010-02-18 | 2024-05-15 | Maui Imaging, Inc. | Point source transmission and speed-of-sound correction using multi-aperture ultrasound imaging |
GB2479200A (en) | 2010-04-01 | 2011-10-05 | Statoil Asa | Interpolating pressure and/or vertical particle velocity data from multi-component marine seismic data including horizontal derivatives |
US8380435B2 (en) | 2010-05-06 | 2013-02-19 | Exxonmobil Upstream Research Company | Windowed statistical analysis for anomaly detection in geophysical datasets |
US8757270B2 (en) | 2010-05-28 | 2014-06-24 | Statoil Petroleum As | Subsea hydrocarbon production system |
US8747325B2 (en) | 2010-07-16 | 2014-06-10 | Fundacao De Amparo A Pesquisa Do Estado De Sao Paulo (Fapesp) | Non-invasive method for diagnosing the severity of heart failure by extracting and analyzing acetone concentrations in captured exhaled breath |
EP2601574A1 (en) * | 2010-08-05 | 2013-06-12 | ExxonMobil Upstream Research Company | Obtaining data from an earth model using functional decriptors |
JP6092109B2 (ja) | 2010-10-13 | 2017-03-08 | マウイ イマギング,インコーポレーテッド | 凹面超音波トランスデューサ及び3dアレイ |
RU2458363C1 (ru) * | 2011-03-16 | 2012-08-10 | Сергей Алексеевич Бахарев | Способ прямого поиска углеводородов |
US8886503B2 (en) | 2011-04-19 | 2014-11-11 | Schlumberger Technology Corporation | Discontinuity detection |
EP2724181B1 (en) * | 2011-06-24 | 2020-08-19 | ION Geophysical Corporation | Method and apparatus for seismic noise reduction |
EP2785253B1 (en) | 2011-12-01 | 2023-11-15 | Maui Imaging, Inc. | Motion detection using ping-based and multiple aperture doppler ultrasound |
JP2015503404A (ja) | 2011-12-29 | 2015-02-02 | マウイ イマギング,インコーポレーテッド | 任意経路のmモード超音波イメージング |
JP6438769B2 (ja) | 2012-02-21 | 2018-12-19 | マウイ イマギング,インコーポレーテッド | 多数開口超音波を用いた物質の硬度の決定 |
US9542507B2 (en) | 2012-02-23 | 2017-01-10 | Schlumberger Technology Corporation | Feature detection in seismic volumes |
KR102103137B1 (ko) | 2012-03-26 | 2020-04-22 | 마우이 이미징, 인코포레이티드 | 가중 인자들을 적용함으로써 초음파 이미지 품질을 향상시키는 시스템들 및 방법들 |
US20130261981A1 (en) * | 2012-04-03 | 2013-10-03 | Westerngeco L.L.C. | Covariance estimation using sparse wavelet representation |
US9261615B2 (en) * | 2012-06-15 | 2016-02-16 | Exxonmobil Upstream Research Company | Seismic anomaly detection using double-windowed statistical analysis |
EP2883079B1 (en) | 2012-08-10 | 2017-09-27 | Maui Imaging, Inc. | Calibration of multiple aperture ultrasound probes |
EP2887879B1 (en) | 2012-09-06 | 2021-05-26 | Maui Imaging, Inc. | Method of ultrasound imaging |
US9303499B2 (en) | 2012-10-18 | 2016-04-05 | Elwha Llc | Systems and methods for enhancing recovery of hydrocarbon deposits |
US9134444B2 (en) * | 2012-11-27 | 2015-09-15 | Chevron U.S.A., Inc. | System and method for deducing cavern properties |
US9213117B2 (en) | 2013-01-02 | 2015-12-15 | Schlumberger Technology Corporation | Dip seismic attributes |
US9829591B1 (en) | 2013-01-07 | 2017-11-28 | IHS Global, Inc. | Determining seismic stratigraphic features using a symmetry attribute |
US9510806B2 (en) | 2013-03-13 | 2016-12-06 | Maui Imaging, Inc. | Alignment of ultrasound transducer arrays and multiple aperture probe assembly |
CN103217719B (zh) * | 2013-04-11 | 2015-07-15 | 中国矿业大学 | 基于单一炮检对观测系统的煤巷超前探测断失翼煤层的方法 |
US9883848B2 (en) | 2013-09-13 | 2018-02-06 | Maui Imaging, Inc. | Ultrasound imaging using apparent point-source transmit transducer |
US10401493B2 (en) | 2014-08-18 | 2019-09-03 | Maui Imaging, Inc. | Network-based ultrasound imaging system |
CN104360386B (zh) * | 2014-12-04 | 2016-09-21 | 中国海洋石油总公司 | 一种针对花岗岩地层划分对比的面元法 |
CN104459801B (zh) * | 2014-12-10 | 2017-04-12 | 中国石油天然气集团公司 | 用于识别断层的相干增强处理方法 |
CN105093298B (zh) * | 2015-07-10 | 2017-06-13 | 北京派特森科技股份有限公司 | 一种微地震数据四维搜索逆时叠加的快速计算方法 |
RU2603828C1 (ru) * | 2015-08-21 | 2016-11-27 | Федеральное государственное бюджетное учреждение "Всероссийский научно-исследовательский геологический нефтяной институт" (ФГБУ "ВНИГНИ") | Способ региональной сейсморазведки слабо изученных осадочных бассейнов для выявления и локализации нефтегазовых зон и объектов |
CN113729764A (zh) | 2016-01-27 | 2021-12-03 | 毛伊图像公司 | 具有稀疏阵列探测器的超声成像 |
CN106094051B (zh) * | 2016-06-08 | 2017-12-15 | 中国海洋石油总公司 | 一种花岗岩地层划分对比方法 |
JP6722521B2 (ja) * | 2016-06-23 | 2020-07-15 | 古野電気株式会社 | 水中探知システム |
US10782433B2 (en) * | 2016-07-29 | 2020-09-22 | Baker Hughes Holdings Llc | Method for an automatic detection of acoustic reflectors and their parameters from borehole acoustic array data |
US10634805B2 (en) * | 2017-05-26 | 2020-04-28 | Chevron U.S.A. Inc. | System and method for predicting fault seal from seismic data |
US11041976B2 (en) | 2017-05-30 | 2021-06-22 | Exxonmobil Upstream Research Company | Method and system for creating and using a subsurface model in hydrocarbon operations |
RU2664503C1 (ru) * | 2017-12-20 | 2018-08-17 | Дмитрий Юрьевич Степанов | Способ формирования куба или разреза площадок, способ автоматического прослеживания горизонтов/годографов и способ автоматического выявления зон тектонических деформаций и зон трещиноватости |
US11609355B2 (en) | 2018-10-02 | 2023-03-21 | Chevron U.S.A. Inc. | System and method for generating an earth model |
US11282269B2 (en) * | 2020-05-13 | 2022-03-22 | Chevron U.S.A. Inc. | Image-comparison based analysis of subsurface representations |
CN113740908B (zh) * | 2020-05-29 | 2024-05-07 | 中国石油化工股份有限公司 | 地震切片的二维变差分析方法、电子设备及介质 |
CN114063154B (zh) * | 2020-07-29 | 2024-05-07 | 中国石油化工股份有限公司 | 走滑断层位移量计算方法、装置、电子设备及介质 |
CN112433248B (zh) * | 2020-11-10 | 2022-09-30 | 成都理工大学 | 一种碳酸盐岩沉积环境下隐蔽储层的探测方法 |
US11906679B2 (en) * | 2021-09-13 | 2024-02-20 | Saudi Arabian Oil Company | Systems and methods for detecting seismic discontinuities by coherence estimation |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3981443A (en) * | 1975-09-10 | 1976-09-21 | Northrop Corporation | Class of transform digital processors for compression of multidimensional data |
US4467461A (en) * | 1981-01-05 | 1984-08-21 | Conoco Inc. | Interactive color analysis of geophysical data |
US4633400A (en) * | 1984-12-21 | 1986-12-30 | Conoco Inc. | Method for waveform feature extraction from seismic signals |
USH374H (en) * | 1987-02-09 | 1987-11-03 | The United States Of America As Represented By The Secretary Of The Army | Optimum multiple target detection and resolution |
US4894795A (en) * | 1987-04-28 | 1990-01-16 | The United States Of America As Represented By The Secretary Of The Navy | High-resolution technique for time-frequency signal analysis using modified wigner-ville analysis |
US4937747A (en) * | 1988-02-16 | 1990-06-26 | Amoco Corporation | Iterative disjoint cluster and discriminant function processing of formation log responses and other data |
US4858199A (en) * | 1988-09-06 | 1989-08-15 | Mobile Oil Corporation | Method and apparatus for cancelling nonstationary sinusoidal noise from seismic data |
US4910716A (en) * | 1989-01-31 | 1990-03-20 | Amoco Corporation | Suppression of coherent noise in seismic data |
US5047991A (en) * | 1989-04-28 | 1991-09-10 | Schlumberger Technology Corporation | Lithology identification using sonic data |
US4951266A (en) * | 1989-04-28 | 1990-08-21 | Schlumberger Technology Corporation | Method of filtering sonic well logging data |
US5031155A (en) * | 1989-04-28 | 1991-07-09 | Schlumberger Technology Corporation | Compression and reconstruction of sonic data |
GB9011836D0 (en) * | 1990-05-25 | 1990-07-18 | Mason Iain M | Seismic surveying |
US5056066A (en) * | 1990-06-25 | 1991-10-08 | Landmark Graphics Corporation | Method for attribute tracking in seismic data |
US5245587A (en) * | 1990-12-14 | 1993-09-14 | Hutson William H | Multi-dimensional signal processing and display |
US5175710A (en) * | 1990-12-14 | 1992-12-29 | Hutson William H | Multi-dimensional data processing and display |
US5253192A (en) * | 1991-11-14 | 1993-10-12 | The Board Of Governors For Higher Education, State Of Rhode Island And Providence Plantations | Signal processing apparatus and method for iteratively determining Arithmetic Fourier Transform |
US5563949A (en) * | 1994-12-12 | 1996-10-08 | Amoco Corporation | Method of seismic signal processing and exploration |
-
1997
- 1997-01-02 RU RU97119642/28A patent/RU2169931C2/ru not_active IP Right Cessation
- 1997-01-02 CN CNB971903492A patent/CN1186647C/zh not_active Expired - Fee Related
- 1997-01-02 WO PCT/US1997/000249 patent/WO1997039367A1/en active IP Right Grant
- 1997-01-02 CA CA002220274A patent/CA2220274C/en not_active Expired - Fee Related
- 1997-01-02 AU AU18247/97A patent/AU710968B2/en not_active Ceased
- 1997-01-02 EP EP97903759A patent/EP0832442B1/en not_active Expired - Lifetime
- 1997-03-11 US US08/815,496 patent/US5892732A/en not_active Expired - Lifetime
- 1997-03-21 AR ARP970101134A patent/AR006341A1/es unknown
- 1997-12-09 NO NO19975784A patent/NO319052B1/no not_active IP Right Cessation
Also Published As
Publication number | Publication date |
---|---|
US5892732A (en) | 1999-04-06 |
CA2220274C (en) | 2005-06-28 |
MX9710022A (es) | 1998-07-31 |
AU1824797A (en) | 1997-11-07 |
NO975784L (no) | 1998-02-12 |
RU2169931C2 (ru) | 2001-06-27 |
NO319052B1 (no) | 2005-06-13 |
EP0832442B1 (en) | 2005-03-16 |
WO1997039367A1 (en) | 1997-10-23 |
AU710968B2 (en) | 1999-09-30 |
NO975784D0 (no) | 1997-12-09 |
AR006341A1 (es) | 1999-08-25 |
CA2220274A1 (en) | 1997-10-23 |
CN1189218A (zh) | 1998-07-29 |
EP0832442A1 (en) | 1998-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1186647C (zh) | 处理和探测地震信号的方法及装置 | |
AU612552B2 (en) | Model-based depth processing of seismic data | |
CN1275048C (zh) | 用于浅流探测的方法 | |
Mallick et al. | Reflection/transmission coefficients and azimuthal anisotropy in marine seismic studies | |
CN102395902B (zh) | 使用快速面向目标照明计算的地震成像系统及方法 | |
US6957146B1 (en) | System for utilizing seismic data to estimate subsurface lithology | |
Bormann et al. | Seismic wave propagation and earth models | |
US8233351B2 (en) | Analyzing borehole seismic data to locate subsurface diffractors | |
US8289809B2 (en) | Common reflection azimuth migration | |
CN101551463B (zh) | 三维观测系统噪声压制估算方法 | |
US6487502B1 (en) | System for estimating the locations of shaley subsurface formations | |
CN1237256A (zh) | 用于测定油层渗透率、孔隙率和伪相对渗透率的分布的方法 | |
CN88101749A (zh) | 多分量地震数据深度成象法 | |
US20080137480A1 (en) | Method of Building a Subsurface Velocity Model | |
EP2548052B1 (en) | System and method of 3d salt flank vsp imaging with transmitted waves | |
AU2008216238B2 (en) | Performing 3-D scatterer imaging from 2-D seismic data | |
CN1886678A (zh) | 向量3分量3维克希霍夫预堆栈迁移 | |
MX2014010606A (es) | Extraccion de datos de cizallamiento sv a partir de datos marinos de la onda p. | |
Pavlis | Imaging the earth with passive seismic arrays | |
Carrion et al. | Tomographic processing of Antartica's data | |
Talagapu | 2D and 3D land seismic data acquisition and seismic data processing | |
Talwani et al. | Some recent developments in the acquisition and processing of seismic data | |
Ganssle | Seismic explorations | |
Vesnaver | Industrial Tomography, a Snapshot | |
MXPA97010022A (en) | Method and apparatus for processing and exploring sismi signals |
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 | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1050492 Country of ref document: HK |
|
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20050126 Termination date: 20110102 |