专利名称：用于地震烃体系分析的方法
技术领域：
本发明一般涉及地球物理勘测领域，更具体而言，涉及地震数据的诠释。具体地，本公开描述利用地震数据来检测和分级潜在的烃机会(opportunities)的方法。
背景技术：
活性烃体系通过提供存储空间给烃类的多孔储油岩层、防止烃类逃逸储层的封闭层(seal)、良好的圈闭几何形态(trapping geometry)和含有高百分比生源物质(biogenic material)的源层(source formation)的存在而被限定。在高温高压的影响下，生源物质被成熟(或被蒸热(cooked))，以形成烃类，包括气、原油、浙青和焦油。由浮力和压力差驱动，烃类运移并且一部分这些烃类聚集在圈闭和封闭层内，所述圈闭由储油岩层的不规则几何排列(即，圈闭几何形态)形成。圈闭具有有限的体积，然而，其可以溢出或泄露一些聚集的烃类，其中一部分然后可以在其它圈闭聚集。地下地震图像允许解释者(interpreter)根据实践和提示性的几何形态鉴别一些潜在的圈闭。有时，地震数据可以提供烃类存在的直接指示。然而，标准解释实践是劳动密集型的，并常常集中于解释者收集某些远景性指示的区域。因此，许多机会仍未，例如通过地震噪音被检测到，因为指示太微妙或隐蔽。即使观察到远景性的线索，在多个明显机会存在的情况下或在解释者被时间约束限制时，它们可能不被检测。因此，一些烃聚集物后来才被发现或仍未被发现。公布的解决类似问题的尝试包括如下“Method for Seismic Interpretation Using Seismic TextureAttributes” (PCT专利申请公开W02010/053618)，Imhof,公开了用于计算纹理属性的方法，所述纹理属性可用于根据其局部外观分类和分段地震数据。纹理可用于限定地震相。“Windowed Statistical Analysis for Anomaly Detection in GeophysicalDatasets” (PCT专利申请公开W02010/056424)，Kumaran等，公开了一种方法，用于针对数据背景下统计学上异常的区域检验地震数据，从而用于突出统计学上不一般的或显著的区域。“Method For Geophysical and Geological Interpretation of SeismicVolumesIn Depth, Time and Age”(美国专利申请公开号 2010/0149917)，Imhof 等，公开了一种方法，用于将来自地球物理深度域或双程走时域的地震数据转换成地质年代域，在所述地质年代域中，所有地震反射大致是水平的并可以与它们在其岩层的地质年代的状态相t匕。从该年代域形成的地震属性可以改进烃体系要素的限定。欧洲专利EP1110103B1 ( “Method Of Seismic Signal Processing”)，Meldahl等，公开了一种方法，用于在地震数据中寻找类似于由解释者指定的区域的区域。此外，它公开了该方法用于寻找渗滤烃类的地震指示中的应用。美国专利号6，226，596B1 ( “Method for analyzing and classifying threedimensional seismic information” )，Gao,公开了一种方法，用于产生地震纹理属性，其可以用于分段、分类或限定地震相。美国专利号6，438，493B1 (“Method for seismic facies interpretation usingtextural analysis and neural networks”)，West 和 May，公开了一种方法，用于产生地震纹理属性，其可以用于监督分类，以指定地震相属性。美国专利号6，516，274B2 ( “Method for imaging discontinuities in seismicdatausing dip-steering”)，Cheng等，公开了一种方法，用于检测地震数据中的边缘(edge)或间断，其常常指示断裂和移动正常地下分层的断层。
美国专利号6，725，174B2 ( “Edge-preserving enhancement of seismic imagesby nonlinear anisotropic diffusion”)，Bouts 等，公开了一种地震数据处理方法，其通过去除固有噪音同时保持与断层或快速地层过渡相应的明显间断诸如嵌入通道(entrenched channel)，来增强地震数据。美国专利号6，850，845B2 ( “System for multi-dimensional data analysis”)，Stark，公开了一种地震数据处理方法，其允许计算用于修整(flattening)的地质年代属性和相关的衍生(derivatory)属性。美国专利号6，850，864B2( “Method for analyzing dip in seismic datavolumes”) ,Gillard等，公开了一种地震数据处理方法，用于评估限定局部结构的反射偏角(reflection dip)。美国专利号7，203，342B2 ( “Image feature extraction，，)，Pedersen，公开了一种地震数据处理方法，用于在相对有噪音的数据中检测和改进断层或层位(horizon)。美国专利号7，424，367B2 ( “Method for predicting lithology and porosityfrom seismic reflection data”)，Saltzer等，公开了一种方法，用于从地震数据中预测地下岩石的岩性学和孔隙度，从而允许区分储集岩和封闭地层。美国专利号7，454，292B2 ( “Inverse-vector method for smoothing dips andazimuths”)，Wang等，公开了一种方法，用于有力地计算限定局部结构的层位偏角和定向。PCT 专利申请公开 W02009/011735 ( “Geologic Features From Curvelet BasedSeismic Attributes”)，Neelamani和Converse，公开了一种方法，用于计算径指不物或纹理属性，其可用于鉴定地下特征。PCT 专利申请公开 W02009/082545A1 ( “Detection Of Features In SeismicImages”)，Kumaran和Wang，公开了一种方法，用于在地震数据中检测断层、通道和类似特征。PCT 专利申请公开 W02009/137150A1( “Method For Geophysical AndStratigraphic Interpretation Using Waveform Anomalies”)，Imhof，公开了一种地震数据处理方法，用于绘制地层终端(termination)和尖灭的图。PCT 专利申请公开 W02009/137228A2 ( “Transport Property Data CalculatedFrom Derivative Seismic Rock Property Data For Transport Modeling”)，Oppert 等，公开了一种方法，用于评估影响烃体系要素的特性诸如热流或流体渗透率。PCT 专利申请公开 W02009/142872A( “Seismic Horizon Skeletonization，，)，Imhof等，公开了一种自动化方法，用于从地震数据集中提取大量层位。此外，它还公开了广泛模式(broad pattern)识别工作流程，其划分数据集、分析区域并根据它们含有烃类的潜在性进行分级。“A New Class of Large-scale Attributes for SeismicStratigraphy”, Gesbert 等，71stEAGE Conference&Exhibition, (2009)公开了一组由二维地震数据计算的地层属性，其突出区域不整合和区域变薄趋势并量化区域地震相趋势。“Applications of plane-wave destruction filters”，Fomel, Geophysics67, I946-1960, (2002)公开了一种方法,用以评估限定局部结构的地震层位偏角和方向。 “Estimating Seismic Heterogeneity with the Structure Tensor”，Imhof, 67thEAGE Conference&Exhibition, (2005)公开了一种方法，用以评估限定局部结构的地震层位偏角和方向和以局部异质性为特征的地震纹理属性。“Flattening without picking”，Lomask 等，Geophysics71, P13-P20 (2006)公开了一种处理地震数据的方法，用以大致修整允许表征烃体系的一些要素的数据。“Hydrocarbon leakage interpreted on seismic data，，，I.oseth等，Marine andPetroleum Geology26, 1304-1319，(2009)公开了解释者驱动的方法,用于检测通过地下渗
漏的烃类。“Hydrocarbon Traps, K. T. Biddle and C. C. ffielchowsky, The PetroleumSystem-From Sourceto Trap, MPG Memoir60, 219-235页，(1994)提出一系列经圈闭类型。“Imaging Vector Fields Using Line Integral Convolution，，，Cabral andLeedom, Proceedings of ACM SigGraph93, 263-270, (1993)公开了一种使流线(flowline)的向量场可见的方法。^Lithofacies Prediction in Deep Water Water Reservoirs”， Oppert 等，Society of Exploration Geophysicists, Expanded Abstracts, 1708-1711, (2006)公开了一种方法，以利用地震和有线数据(wireline data)评估地下岩性学。“Seismic Rock-Property Inversion and Lithofacies Prediction at ErhaField, Nigeria，，, Xu 等,Nigerian Association of Petroleum Explorationists (NAPE) 2008conference proceedings,公开了一种方法，以利用地震和有线数据评估地下岩性学。Randen 和 Sonneland( “Atlas of3D Seismic Attributes”，in MathematicalMethods and ModelinginHydrocarbon Exploration and Production, Iske andRanden (编著)，Springer,第23-46页(2005))提出表征地震纹理或地震地层特征的三维地震属性的综述。所需要的是自动化的系统，其扫描烃体系要素的整个数据集，并输出远景构造列表给解释者，以进行检验。优选地，该潜在目标的列表通过烃体系要素的期望体积、存在和性质以及其检测和鉴别的置信度(confidence)进行分级。优选地，远景构造列表也被评注。本发明满足至少这些要求
发明内容
本发明的一方面是一种方法，其计算多尺度(multi-scale)和通常为结构导向的(structure-oriented)地震属性,所述地震属性涉及典型烃体系的要素，即储层、封闭层、圈闭、源、熟化和运移。属性在空间上相关，并与烃圈闭构造的目录(catalogue)进行比较，以确定烃圈闭的潜在存在并评估其存在的置信度。在一个实施方式中，本发明是一种计算机执行的方法，用于针对烃体系或特定远景带(particular play)的存在分析由代表地下区域的地震数据的体素组成的体(volume),所述方法包括划分地震数据体，以形成多个段；和至少部分根据每一段中地震数据体素的远景性分数(prospectivity score),针对烃体系或特定远景带的存在对多个段分级；其中所述远景性分数是基于对至少两个属性的计算，所述至少两个属性涉及烃体系或特定远景带的不同要素。 与任何地球物理学数据处理方法一样，本发明在实际应用中是高度自动化的，SP，是在根据本文公开内容而程序化的计算机的辅助下执行的。
附图简介本发明及其优势将通过参考以下详细描述和附图而被更好地理解，所述附图中图I图解背斜圈闭的烃体系要素；图2是流程图，显示在本发明一个实施方式中的基本步骤；图3是示意图，图解本发明方法的一个实施方式；图4图解线积分卷积过滤器的过滤模板(filter stencil)如何遵循切线方向；图5图解如何可以使得线积分卷积过滤器更有力；图6图解线积分卷积过滤模板可以具有扩充的厚度；图7显示将局部属性转换成区域属性的实例(左地震数据；左中局部收敛(convergence)属性；右中切线向量场；和右区域收敛；图8显示具有收敛量度(magnitude)分量(component)(浅的强收敛；深的没有收敛)和收敛方向分量的三维收敛属性的实例；图9限定地层终端关系；
图10描述褶皱界定的(背斜)构造圈闭；图11描述断层相关的构造圈闭；图12显示盐相关的构造圈闭；图13图解与储层尖灭相关的地层圈闭；图14显示与不整合(unconformity)有关的地层圈闭；图15描述由埋藏的侵蚀地形(relief)形成的地层圈闭；图16描述由成岩差异形成的地层圈闭；图17呈现示意性沉积层序模型，即，段塞(slug)图；图18图解不同的要素如何可以被积分，以评估烃聚集的成功机会；图19显示分级四个潜在目标的示意性应用；和图20图解本发明方法的实例应用。
本发明将结合实例的实施方式被描述。然而，在一定程度上，以下详细说明是针对本发明的具体实施方式
或具体用途的，这意图仅是说明性的而不应该被解释为对本发明范围的限制。相反，意图涵盖由所附权利要求书限定的、可包含在本发明范围内的所有替代方案、修饰和等同物。实例的实施方式的详细描述本发明的一个方面包括计算涉及烃体系要素的地震属性、分析烃体系要素的数据和/或识别具体远景带或概念(conceptual)类型以及根据远景构造(prospect)的大小、性质和置信度对这些区域进行分级和评注。接下来给出一些定义。尽管术语可能被更广义或狭义地用于其它地方，但石油或烃体系在本文中通常用于表示天然体系，其包括活性源岩矿体(pod)和所有相关的油和气。它包括对于烃聚集来说必需要存在的所有地质要素(地质单元,geologicelement)和过程，如图I所示。在自·然界中发现的烃类包括高浓度的热和/或生物成因气，其发现于常规储层或气水合物中、紧密的储层、断裂的页岩或煤；以及冷凝物、原油、重油、浙青和焦油。术语“体系”描述相互依赖的要素和过程，其形成产生烃聚集的功能单元。基本要素包括石油源岩(源)、储集岩(储层)、密封岩(封闭层)和上覆岩层岩(上覆岩层)。该过程是圈闭和熟化(生油(generation))的形成、运移和烃类的聚集。这些过程暗含事件的顺序或时间(timing)。烃体系的可选定义可以仅包括源岩、熟化和运移过程及其时间，在这种情况下，储层、封闭层和圈闭可以被限定为形成远景带。为了说明本发明方法的目的，术语烃体系被限定为涵盖源、储层、封闭层、圈闭、熟化、运移和时间。此外，术语远景带在本文中通常被用于表示储层、封闭层和圈闭几何形态的特定组合和排列。源岩是富含有机物质的岩石，其如果被充分加热，将随着时间产生油和/或气。普通源岩包括页岩或石灰岩。海成岩趋于易生油，而陆地源岩(如煤)趋于易生气。保存有机物质而不降解对于产生好的源岩是关键的，并且对于完整的石油系统是必要的。储层是地下的岩石主体，其具有足够的孔隙度和渗透性，以接收、存储和输送流体。沉积岩是最常见的储集岩，因为它们比大部分火成岩和变质岩具有更大的孔隙度，并且在可以保存烃类的温度条件下形成。储层是完整石油系统的重要组成。封闭层是相对不可渗透的岩石，通常为页岩、无水石膏或盐，其在储集岩上并部分围绕储集岩形成屏障或覆盖层，以便流体不能运移到储层外面。封闭层是完整石油系统的重要组成。上覆岩层是在源和储层上面的岩石。在石油系统的背景下，其主要功能是在源上形成厚的覆层，其中它将温度和压力提高到将有机物质转化成烃类必需的程度。圈闭是适合包含烃类并被相对不可渗透的地层封闭的岩石构造，烃类不能通过相对不可渗透的地层运移。圈闭被描述为构造圈闭(在变形的岩层诸如褶皱和断层中)或地层圈闭(在岩石类型变化的区域，如不整合、尖灭和矿脉(reef))或其组合。对于构造圈闭，变形必需发生在烃运移之前，或者烃类将不会聚集。圈闭是石油系统的重要组成。生油或熟化是从源岩形成烃类，此时浙青从干酪根形成，并聚集为油或气。生油取决于三个主要因素足以产生烃类的丰富有机物质的存在、足够的温度，以及足以使源岩成熟的温度和时间。压力以及细菌和催化剂的存在也影响生油。例如通过薄的上覆岩层进行浅埋藏引起的不充足的压力和温度致使源不成熟，而且生油将缺乏或不完整。例如通过在厚的上覆岩层下进行的深埋藏引起的过高的压力和温度将导致产生的油降解成气，随后成为二氧化碳和水。生油是石油系统开发中的关键阶段。运移是烃类从其源到储集岩中的运动。新产生的烃类在其源岩外的运动是初级运移，也称为驱逐(expulsion)。在烃圈闭或其它聚集区域中烃类向储集岩的进一步运移是次级运移。运移通常从构造上低的区域到较高的区域发生，这是因为与周围岩石相比烃类的相对浮力。运移可以是局部的，或者可以在大的沉积盆地中沿着数百公里的距离发生，其对于形成可行的石油系统是重要的。聚集既指圈闭的烃类，S卩，远景带或油田或气田的发生，又指开发石油系统的阶段，在该阶段中，烃类运移到储集岩中并保持圈闭在储集岩中。时间指其中要素被形成或改变的相对顺序，或其中过程发生的顺序。只要在运移之前形成，圈闭就可以聚集运移的烃类。如果运移尚未到达其位置，圈闭可以不被填充。如果封闭层在聚集之后破坏，则圈闭可能至少部分地失去其填充物(charge)，。
远景带是烃聚集方式的概念模型，其常常用于在盆地、区域或趋势(trend)中发展远景构造或用于继续开采鉴定的趋势。远景带(或一组相关的远景带)通常出现在单一的烃体系中，并可以由一组类似的远景构造组成。远景构造是其中已经预测存在经济量的烃类的区域。远景构造常常是一种异常，诸如地质构造或地震振幅异常，其被推荐为进行钻井的位置，以探知经济量的烃类。通过收集活性烃体系的证据或证明遇见良好性质的储集岩、足够大小的圈闭、足够密封的岩石和用于产生和运移烃类以填充储层的适当条件的合理可能性，作出钻凿远景构造的论证。对于本发明方法的目的，远景构造被广泛地用于表示被推荐用于进一步详细分析的区域。如上所述，所需要的是一种方法，其针对烃体系要素的存在自动分析地震数据、标记远景带要素以有利的构造并置或与已知的或指定的远景带一致的区域，并且针对其烃聚集的潜在性对这些远景构造进行分级。这种体系关注于对更多远景区的分析和解释。此外，该系统还可以识别远景带类型并提供单个要素的置信分数。对于每一远景构造，未被鉴定的要素或具有低置信度的要素是重要的，其在随后的分析中需要特别关注，以化解远景构造的风险。该系统在区域或资源的周期中可以被应用很多次。首先，该系统可用在区域性、通常为二维的数据上，以鉴定远景区；例如，针对区块(block)进行投标(bid)或确定三维勘探调查的地点。在勘探阶段，该系统可用于范围界定，以使解释者关注并引向更小、更容易处理的地震数据体的子集。在生产阶段，该系统可用于将更小的远景构造定位于已经定位的基础结构的有效范围内。最后，该系统可用于确保在放弃或出售资源之前，没有远景构造被疏忽。图2是流程图，显示本发明方法的一个实施方式中的基本步骤，包括两组步骤(具有实线的方块)和六个可选步骤(具有虚线的方块)，这取决于输入、进行的假设和该系统被使用的模式。在步骤2中，产生属性，其涉及烃体系要素或具体远景带要素。在步骤5中，分析数据并针对属性对位置进行分级，从而鉴定潜在的远景构造，其然后被储存用于进一步分析或可视化。可选步骤I通过断层和褶皱限定地下组构(fabric)，即，地下岩层的几何形态及其变形，其形成用于计算一些属性的基础或可能被其它计算方法需要以从类似岩层中集合信息。许多属性对于本发明系统是有用的，其通过沿着组构积分或平均化更多的传统属性而形成。在可选步骤3中，数据被至少分成两个段或部分，用于进行分析和远景构造限定。可选步骤4是利用一个或多个属性形成烃体系要素的归一化分数。可选步骤6是分析形成的和分级的远景构造；例如，鉴别应该被进一步检验的、最少表达的要素(一个或多个)。可选步骤7是限定或选择至少一个将被搜索的具体概念或远景带类型，而不是烃体系邻近要素的通用搜索。最后，可选步骤8是限定、产生和管理远景带类型及其要素构造的目录。图3呈现本发明方法的示意性应用，其中地震属性32由地震数据31计算。然后，合并属性33，以突出可能构成地层远景带的区域(例如36)，在这种情况下，地层远景带为埋藏山远景带34和地层尖灭远景带35。组构限定属性是地震数据的可测量的性质，诸如振幅、偏角、频率、相和极性。属性可以在瞬间被及时测量或在时间窗口(window)内被测量,并且可以在单一迹(trace)或一组迹上或在由地震数据、数据窗口解释的面上，或者甚至在多个地震体上同时进行测量。本发明方法 应用许多典型的属性，其对于属性分析或地震解释领域的有经验的人来说是悉知的。然而，对于本发明方法，这些属性中的一些被修改，例如通过沿着岩层或结构进行积分，以突出相对于局部的区域变化。此外，公开了新的地质属性(2)，其直接涉及烃体系要素或远景带要素以及这些要素的具体构造。由于地质属性可以以交替方式进行体现，并且不同的地质属性可能涉及相同的要素或要素构造，所以可以计算分数(4)，其将不同实施方式和属性一起综合成一次测量，评估遇见特定远景带或烃体系要素或这些要素的具体几何排列的机会。然后，将分数合并(6)成构形的似然性，其可用于可视化、分析或评价(rating)和分级。对于图2，可能需要比示出的所有步骤少的步骤来实施本发明的具体实施方式
。单独的步骤可以合并或分成多个分量。此外，另外的和/或可选的方法可应用本文没有显示的另外的步骤。虽然流程图显示顺次发生的各种操作(action)，但也应该理解，一些操作可以连续发生、基本上并列和/或在相当不同的时间点发生。也可以重复步骤。这些变化的一个实例是步骤1，对组构、结构或分层构造的限定。本发明中应用的一些属性可能不需要组构，因此，在仅利用这些属性时，可以忽略组构限定。然而，用于本发明方法的大部分属性将需要用于其形成的组构。对于每一需要的属性，可以快速(on the fly)计算组构。然而，优选地，组构仅被产生一次，然后被存储并重复使用，以产生、修改或积分属性。步骤1，组构的形成，其被显示为任选的步骤，因为根据在本发明方法的不同实施方式中使用的属性，组构可以被计算一次、被重复计算或根本不计算。产生一次组构的优势是不同属性之间的一致性，因为所有均基于相同的组构。产生它仅一次的缺点是所得组构需要以某些方式被呈现和储存，这会是效率低的，或者甚至不适于属性产生、修改或积分的一些应用。可选地，每一应用均可以以合适和有效的表不产生其自身的组构，这会引起整体系统的低效，因为具有相同表示的相同组构可能被计算多次。再一个选择是用少数选择方法和表示计算组构，所述选择方法和表示被储存，因而是可重复使用的。对于这种可选方法，缺点是提高的储存要求和重复访问存储介质。在实践中，优选的方法可以是计算和存储一些更常见的方法和表示的组构或其分量，在需要时计算较不常见的那些。计算组构的方法可以基于偏角、方向或走向(strike)评估、构造张量(tensor)、波形相关性或骨架化(skeletonization)。估计偏角和偏角方向的一个具体方式是通过在两个水平方向和垂直方向的地震数据上使用梯度(gradient)(例如，美国专利号 6,850,864B2, “Method for analyzing dip in seismic data volumes，，，Gillard等)。另一偏角估计器是基于平面波解构过滤器(deconstruction filter)(例如，“Applications of plane-wave destruction filters^ by Fomel, Geophysics67(6), 1946-1960，(2002))。基于构造张量的方法允许计算地震反射的向量法线(normal)和切线(例如，“Estimating Seismic Heterogeneity with the Structure Tensor”，Imhof, 67thEAGE Conference&Exhibition, (2005);或美国专利号 7，454, 292B2, ^ Inverse-vectormethod for smoothing dips and azimuths”, Wang 等)。偏角可以被分解成(resolved)时移，用于面的自动采集(例如，“Flattening without picking”，Lomask等，Geophysics71(4)，P13-P20(2006))。计算反射偏角和方位的可选方法是基于等高线(isocontour)的梯度。梯度是局
部向量，其表征函数的最速增量方向。处理地震数据，例如振幅，其在三维空间(x，y，t)中·作为函数f，梯度被定义为因为在所得梯度向量中有三个分量，所以仅需要三个邻近
样品来评估一阶(first order) Y/:然而,利用更多的周围样品允许评估超定系统的梯度，例如，
在最小平方的意义上，或应用高阶逼近到有限差分逼近。在任意情况下，所得梯度将更平
滑。邻近点的选择不限于任何具体的层序或模型。相反地，邻域可被限定为点或体素，其以
球形、圆柱形、框或围绕感兴趣的点的任何其它类型的形状被排列。此外，该形状可以在每
一点与组构或其估计量排列。然后，可以从该梯度产生衍生的属性，包括量度和方向，即，3D
的偏角和方位。如通常所限定的，梯度通常指向振幅增加的方向。在振幅从负值谷循环到
正值峰并返回的地震数据中，结果可能是平行层显示相反的偏角值。一个方案是通过翻转
(flipping)其方向针对一个半空间归一化梯度，例如转换梯度的符号，以便所有梯度点定
向在正的y方向
a* /:+ sgn(/r)f · = Sgn(A)- mctm^=====^ = arctan^==^ = arctani^ e
机 +Λ+ΛΛ该梯度计算是比现有的梯度计算方法更有力的公式。它还具有多尺度能力，因为其中梯度被计算的窗口尺寸(邻近点的数目)可以变化。波形的横向相关是产生组构的另一方法。在第一步骤中，事件诸如峰、谷和/或零交叉被测定。在第二步骤中，将围绕这些事件的波形与邻迹的事件相比较，以找到并关联最相似的那些。注意，在每一事件中发生的不对称性仅与一个其它事件相关联，即，最相似的那个。然而，多个事件可以与一个和相同的事件关联。结果是有向图结构(directed graphstructure),其中顶点(vertice)与事件诸如峰、谷和/或零交叉对应；和边缘(edge)使单个事件与最相似的事件在邻迹中关联。以具体迹中的一个事件开始，常常可能沿着事件和迹之间定向的边缘，并稍后返回到原始迹中的不同事件，这表明地层(或拓扑)不一致性。为了区分该方法与其它方法，在必要时将使用术语原始骨架(skeleton)或原始骨架化。用于计算组构的优选方法是通过拓扑学骨架化来进行的，其以拓扑学上一致的方式自动产生基于反射的面，其中，单个面本身并不重叠，并且多个面的集合与地层的重叠原理一致(PCT专利申请公开W02009/142872A1, “Seismic HorizonSkeletonization”，Imhof等，)。为了区分该方法与其它方法,在必要时将使用术语拓扑骨架或拓扑骨架化。在该方法中，拓扑学上一致的面以自上而下的方式被单调地标记，这允许将假定的年代指定给地震数据的每一样本以及将地震数据从传统的地球物理学深度或双程时间域转换为地质年代域(例如，美国专利号6，850，845B2，“System formulti-dimensional dataanalysis”，Stark ;PCT专利申请公开W02009/142872A1,“SeismicHorizonSkeletonization” , Imhof等；和继续在美国专利申请公开号2010/0149917中的“Method for Geophysical and Geological Interpretation of Seismic VolumesinDepth, Time, and Age”，Imhof 等)中。所得年代映射体(mapping volume)和 / 或深度映射体允许数据在地球物理学深度(或双程时间)域和构建的地质年代域之间转换，反之亦然。对于具体层位，经过深度映射体的每一水平切面对应于深度(或时间)结构图。再一个骨架化方法是基于地貌变薄；结果将被称为地貌骨架。地震数据被区块化成二值图像，例如，通过将数据减少到仅仅它们的极性，例如，±1。具有+1值的带(Band)以相等的连通性被减小到具有+1值的线(line)，例如通过应用地貌变薄。具有-I值的带被较小到具有-I值的线。所有其它样品被设为零。从明显的极性属性中获得类似的结果，所述极性属性通过在局部振幅极值处计算的瞬间相的极性形成。
计算组构的另一方面涉及断层，其被检测为地震数据的间断性(例如，美国专利号 6，516,274B2, “Method for imaging discontinuities in seismic data usingdip-steering”，Cheng 等；或 PCT 专利申请公开 W02009/082545A1 ( “Detection OfFeatures In Seismic Images”，Kumaran和Wang)。检测的断层段可以进一步被清除或细化(例如，美国专利号 7，203，342B2, “Image feature extraction，，，Pedersen)。地质属性I. a通过长距离结构导向的修勻(Oriented Smoothing)产生区域属性地质属性是地震属性，其突出或量化烃体系或远景带的一些方面。地质属性在自然界中常常是区域性的。例如，它们可以相对于数十百或数百的体素评估地震数据性质。相对于这种距离，优选以结构导向的方式沿着组构计算属性，以防止混合，从而防止通过分析窗口切割的其它岩层的污染。在给出由原始地震数据体所限定的组构的情况下，任何属性均可通过沿着组构积分或平均化而被转换成区域属性。一个这样的方法是基于非线性、各向异性扩散过滤(例如，美国专利号 6，725，174B2, “Edge-preserving enhancement of seismic images bynonlinear anisotropic diffusion”, Bouts等)，其主要沿着反射的方向使地震数据修勻。传统地，利用给定的待被修匀的地震数据集的地震构造张量计算组构或反射方向，但组构可以通过任何方法进行限定。本文公开了一些方法。在本发明中，非线性、各向异性扩散可用于计算初级(primary)地震振幅体的组构，但是，随后该组构被用于沿着初级地震振幅体中反射的方向修勻次级(secondary)数据集。在一些情况下，该次级数据集可能甚至不被分层或片段化(banded)，这防止计算其本身的正确构造张量。在次级数据集片段化的情况下，然而，沿着参考地震体的组构使其修匀可能是有利的，例如在结构导向的修匀之后增强数据集之间的一致性。本发明中在长距离上结构导向的修匀的可选方法是基于线积分卷积(LIC)。线积分卷积是悉知的纹理合成算法，其用于图像处理或数据可视化中("Imaging VectorFields Using Line Integral Convolution", Cabral 和 Leedom, Proceedings of ACMSigGraph93, 263-270，(1993))。低通过滤器(low-pass filter)被用于沿着以像素为对称中心的双向流线来卷积输入的噪音纹理，例如随机噪音图像，以在流动方向上发掘空间相关性。通常，流线不是由曲线表示而是由向量场表示，所述向量场的向量是流线的切线并且流线通过积分进行计算。用一些地震属性数据取代噪音纹理并利用组构限定流线允许利用线积分卷积来沿着组构过滤给定的属性，这是进行结构导向的过滤的一种新方法。优选的组构限定是基于偏角，例如，由构造张量计算得到的。通过离散组构，或者更具体而言，使偏角或反射切线离散成0、45、90、1 35、180、225、270或315°，获得进一步的效率；即，完全横向的、完全垂直的或完全对角线的，其用沿着链表的运动取代数值积分和关联的插值(associatedinterpolation)。由于地震数据通常或多或少被水平片段化，所以通过使所有切线从左到右延伸来获得进一步的效率，这将偏角减小到O、+45和-45°。图4图解线卷积过滤器的应用，其中灰色十字表示样品的位置，而箭头表示离散的偏角方向。过滤器在由黑色粗点表示的位置处被应用。附图显示以点为中心的过滤模板如何遵循向量，因而显示获取样品值的组构，其沿着这些途径用黑色十字进行标记。在断层和其它复杂性的存在下为了使得线积分卷积过滤器更有力，可以通过在多个模板构造中应用过滤器而进行边缘保存，所述模板构造不一定以当前的样本位置为中心(如图5所示)。在所有测试的模板构造中，在其样本中显示最少方差的那个模板构造优选最终应用在特定位置。取决于过滤系数和通过这种边缘保存的线卷积积分过滤器处理的属性的性质，除了最小方差以外的标准可以是合适的。其它可能的标准包括最大方差、最小或最大输出值和/或避免特别标记的样本位置。过滤系数并不需要是一致的，这允许应用微分器(differentiator)(高通过滤器)、积分器(integrator)(低通过滤器)或其组合。过滤器可以是线性的或非线性的，例如中值过滤器或模式过滤器。图6图解过滤模板可以具有任意的形状，所述任意的形状被变形，以顺应经过当前样本位置的向量路径。图7显示这种过滤工具的实例应用，其将局部属性转换成区域属性。用于在长距离上执行结构导向的过滤的再一个方法包括利用(exploitation)通过地震骨架化获得的面(PCT专利申请号W02009/142872A1, “Seismic HorizonSkeletonization”，Imhof等),其中，骨架面提供组构，沿着组构，属性被过滤。不同于沿着离散面进行修匀一这可能需要插值，可以沿着由骨架面构建的年代映射体的等值面(isosurface)进行修勻。可选地，属性可以从传统的深度或地球物理双程走时域转换到地质年代域。在地质年代域中，水平切面对应于一个层位，因而，通过定向的限定结构对水平切面进行修勻或过滤(美国专利申请号12/623，034 “Method For Geophysical andGeological Interpretation of Seismic Volumes In Depth, Time, and Age”，Imhof等)。在年代域中进行过滤之后，可将结果转换回到深度或双程走时域。与用于结构导向的过滤的方法无关，地质属性可以作为多尺度属性被计算，因为过滤器大小由解释者或该系统指定。不同于仅以一个尺度计算地质属性，解释者或系统可以选择多种过滤器尺寸，这允许以不同尺度计算相同的地质属性。Lb黑塞属性本发明的新地震属性是黑塞或二阶空间导数(双积)张量//=▽▽；其以类似于修改的梯度的方式被计算。在每一个感兴趣的位置处，三维的黑塞可以通过对称的3X3矩阵来表示，即，通过六个独立的分量来表示。为了计算任何体素位置处的黑塞，需要至少六个相邻体素，因为它具有六个独立分量。然而，计算可以包括更多个体素，以获得对黑塞更平滑和更稳定的评估通过从超定系统计算其分量，例如从最小平方的意义上，或应用空间导数的高阶逼近到有限差分逼近。邻近点的选择并不限于任何具体的层序或模式。相反，邻域可被限定为以球形、圆柱形、框或围绕感兴趣的点的任何其它类型的形状排列的点或体素。此外，该形状可以在每一个点被排列成组构或其估计量。从几何观点来看，地震数据体由点特征(例如，噪音)、曲线特征(例如，终端或通道)、面特征(例如，断层)或体特征(例如，通道带或盐体)组成，其嵌入在分层的背景结构(background structure)中。在当前的实践中,振幅或迹相关性信号的方向二阶导数通常被计算，以检测异质性，所述异质性通常发生于地质特征诸如断层、侵蚀边界和尖灭。缺乏交叉导数，基于传统二阶导数的属性不完全地代表数据中的局部变化。通过其六个分量，黑塞表征地震数据的局部凸度或凹度。在黑塞的对角线上发现传统的二阶导数，而在非对角线上发现交叉导数。因而，属性可由这六个分量及其组合形成。分量，以及因而属性由地震体的二阶导数形成，并显示梯度的空间变化。它们反映地震体中偏角和方位角的空间变化，这允许突出异常区域和相互关联的特征组。通过突出不同类型的“间断性”及其空间相 互关系，地震体可以被更好地表征和可视化。黑塞的另一用途是用作计算空间曲率的中间步骤，如在下一部分所阐明的。另外，黑塞矩阵的本征值以及黑塞在一定面上的投影也形成对等高线有多弯曲的度量。曲率测量对象，例如面偏离平面的程度。然而，具体的限定取决于上下文。例如，可以针对平面曲线或空间曲线限定曲率，在这种情况下，曲率是标量(scalar)。也可以针对嵌入到三维空间中的参数化的面限定曲率，在这种情况下，曲率是张量。以发现与地层学有关的属性为目标，在本文中限定两种新类型的曲率。第一种是二维地震图像中的等高线曲率。经过地震体的任意二维截面被处理为函数映射(map of a function) f (x, y)。等高线是该截面中的曲线，其坐标满足f U，y)=常数。其曲率被限定为
.......fjl +/;,/； m
ιο 131 κ =-:-----—
(Λ2+Λ2)3+2其中，所有量均通过等高线梯度向量或黑塞张量的分量限定。该曲率属性在突出高曲率的几何特征诸如终端、角(corner)和接合点(junction)中特别有用。.等高线曲率的另一限定是基于三维的广义化(generalization)。地震数据被处理为三个独立变量的函数f (X，y, t)。所有满足f (x, y, t)=常数的体素构成等值面，即，嵌入到三维空间的曲面。对于弯曲空间，曲率是张量量,其允许对与曲率诸如主曲率、中曲率和高斯曲率有关的多个标量量开方。对于这种曲面，可以通过归一化等高线梯度限定曲面法线n = V//|V/|。可以利用条件t ·η=0限定表面切线。注意，存在无穷数目的切线，其在面上任何给定的点处满足该条件。并且，存在恰好一样多的、与由切线t和法线η限定的面正交的平面。这种平面与面的相交是曲线，其曲率被称为法曲率。所有法曲率中的最大值和最小值被称为主曲率K1和K2，通过K1和K2将高斯曲率计算为积K = K1^K2,同时其代数平均数获得中曲率M= (K^K2)/2ο所有这些量均针对空间中的每一点被限定，并且它们可以被简洁地表达为矩阵S，称为形状算子。为了建立S的一个常规实施方式，注意在给定点的所有切线的空间(space)可以被表示为矩阵T的本征向量，该矩阵T通过从单位矩阵减去法线与其本身的外积而获得，T=(I-nn)。T有效地使任何三维向量在切线平面上投影。因而，通过以T乘以黑塞矩阵，或象征性地H' =H ·Τ，可以测量限制在该平面的梯度的变化，其中， 表示内部的或矩阵、积。最后，通过归一化黑塞获得S :S=H7|砂I2.主曲率是S的本征值，S的迹除以2获得中曲率，M=迹(S)/2,并且行列式(determinant) det (S) =K限定高斯曲率。所有这些曲率测量均可被处理为地震属性，并可用于检测层位终端和接合点，其限定具有地质意义的特征，诸如断层、削蚀作用、尖灭和其它几何关系。等高线属性的形式为向量和二阶张量。这些属性在特定方向和面上的投影可以比任何其单个分量更好地表征地震体。梯度在任何方向上的投影仅仅是梯度向量和单位方向向量之间的点积。这种投影沿着特定方向限定等高线函数的方向导数。该投影(或方向导数)越大，振幅沿该方向的变化越大。类似地，黑塞在单位方向向量d上的投影是等高线函数沿该方向的二阶导数，并可以被计算为fdd=d · H · d，其中·表示内部的或点、积。对于黑塞，特别感兴趣的投影沿着梯度方向、偏角方向或与由梯度和偏角均跨越的平面正交的方向。等高线偏角体(dip volume)的梯度及其投影也包含有用的地下构造信息。偏角梯度量化原始体中倾斜(dipping)方向的变化。偏角梯度沿偏角方向的投影是对偏角变化速率的度量。这帮助检测反射物(reflector)诸如脊、丘、洞和谷上的弯曲特征。偏角梯度在向量——其是偏角方向的法线——上的投影限定对反射物收敛或发散(divergence)速
率的度量。虽然这些等高线属性(梯度、偏角、方位、黑塞、曲率和投影)可用作局部属性，但·结果可被噪音污染的区域或具有高度复杂性的区域中的假象(artifact)影响。修匀减少这种假象，并且例如可以通过用宽的低通过滤器进行卷积或应用有大窗口的中值过滤器执行。然而，优选地，沿着组构在较长距离上以结构导向的方式进行修匀，以防止跨层混合。通过应用长距离、结构导向的过滤器，局部等高线属性可被转换成地质属性。I. c 收敛本发明的另一新的地质属性是收敛，其检测地震反射收敛以形成反射组的区域变薄的趋势的地方，从而突出尖灭、楔形(wedge)或地层的上超和下超。收敛属性可用于向解释者强调地层特征或作为输入数据进入到辅助或自动模式识别系统，诸如本发明的地震烃体系分析系统。计算收敛的两种方法是通过使用黑塞或通过使用组构或偏角法线方向上的偏角的梯度。根据用以计算黑塞或梯度的邻域，沿着组构的结构导向的长距离修匀可以利用线积分卷积、扩散过滤器、骨架引导的过滤器或在年代域中的过滤被应用。这些当中的每一个均利用原始地震数据，以沿着反射物引导过滤器。计算收敛属性的另一方法是首先测量局部反射物厚度，然后评估局部横向的厚度导数，最后沿着反射进行横向修匀操作，以获得区域收敛趋势。如果在单一的二维截面上评估趋势，则收敛可以向左或向右发生，或可选地，被称为沿着一个参考方向的收敛和发散，例如左到右。如果通过从3-D数据集提取的正交切面(slice)评估趋势，则收敛可以与方位和量度有关。此外，也可以在体积意义上利用三维分析窗口而不是多重正交二维分析窗口来评估趋势。
评估局部反射物厚度的方法包括回路持续时间(loop duration),并置的极值或零交叉之间的距离；瞬时频率，或地震骨架的两个并置的面之间的时间差异。连续属性诸如瞬时频率可能需要被区块化。任何这种厚度属性均可被横向区别，以评估局部厚度变化的速率。也可以通过向量场的发散，例如反射法线，或通过张量场的发散，例如反射切线来评估局部厚度变化的速率。为了获得区域收敛趋势，利用线积分卷积、扩散过滤器、骨架引导的过滤器、在年代域中过滤或使用原始地震数据以沿着反射物引导过滤的任何其它方法，用结构导向的长尺度过滤器算子(long-scale filter operator)沿着组构对局部厚度导数进行修勻。对于计算效率，标记或层位指示物可用导数进行编码。图7显示收敛属性，其基于回路持续时间、横向导数并利用线积分卷积，以增强长距离收敛趋势。通过在纵测线和横测线上同时执行该程序，可以构建具有量度和方位的三维收敛属性(图8)。 I. d 汇流(confluence)本发明的再一个新的地质属性被称为汇流。地震反射被处理为携带信息量(traffic)的途径网络，并且这些途径中的一些被拥塞，因为它们与许多其它途径连接。汇流通过评估有多少路径经过任何给定的体素来测量该拥塞。就地震地层学而言，分支和合并反射产生终端，其形式为下超、上超、顶超和削截(truncation)，如图9所描述。所有这些终端潜在地与不整合有关；下超也可能与下超面有关，即，溢流(flooding)或海侵面。用于计算汇流的一组方法是基于由反射切线指定的流线(flow line)的密度。流线或者通过切向场的积分进行计算或者仅通过遵循离散的切线方向链进行计算。例如，以最左边的迹开始，新的流线可以朝右被起动(initiated)和扩展(propagated)。新的流线或者在有规律地间隔开的深度或双程走时的增量处，或者在特定的事件诸如峰(最大)、谷(最小)和/或零交叉处被起动。流线上的样本，优选在围绕流线的小地带(zone)中的样本被标记成属于该特定的流线。在第一迹中的所有流线均被起动和标记后，算法进行到右边的邻迹并开始起动和标记另外的流线——或者在尚未标记的区域中的有规律的增量处或者在尚未标记的事件中的有规律的增量处。流线可以被追踪和标记到左边或右边；优选地，流线被计算和标记到当前迹的左边和右边。在所有未标记的区域或事件已经被用于起动流线并因而被标记后，该过程在下一个迹进行重复。当每一个样本或事件已经被标记至少一次时，通过计算每一样本被标记多少次来测定汇流。优选地，通过仅仅增加与样本位置有关的计数器，使该计算与样本/事件标记同时进行。虽然这些步骤可以从左到右进行，但优选地，它们也从右到左被重复，然后被加入到之前的结果中，确保对称性。在围绕流线的区域中的样本的标记可以用常数指示物或用根据其与流线的距离而变化的值进行。变量指示物权重，例如可以类似于三角形或高斯。忽略什么时候以及在哪里起动新的流线的问题，可以简单地标记和计数流线而不用标记任何区域，之后，进行修匀操作，例如通过矩形波、三角形或高斯曲率过滤器。不总是从左到右追踪流线(反之亦然)，可以从起动点仅计算流线有限的距离。该计算可以通过线积分卷积进行，或者优选地，通过其变化——仅仅增加模板接触的样本的计数器——来进行。对三维数据的这些基于流线的汇流的实施方式的广义化可以通过沿着纵测线和横测线方向或沿着一组任意的方向独立地进行计算来获得。可选地，流线(曲线)可以被广义化成流片(flow sheet)(面)，或线积分可以被广义化成面积分。汇流的其它实施方式可以利用地震骨架获得。原始骨架由事件(峰、谷和零交叉顶点)的有向图和与邻迹(边缘)中的事件的波形相关来表示。通过沿着边缘(相关性)，图可以从特定的顶点(事件)遍历到左边和/或到右边。从不同的顶点起动图形遍历允许计数任何给定的顶点多久能被通过。每一个顶点均被用于起动遍历，或者优选地，仅尚未被通过的顶点被用于起动新的遍历。在优选的情况下，首先，每个顶点将起动遍历，但逐渐地，顶点将被通过至少一次，并且新遍历的起动将变得越来越少。汇流的这种特定的图形遍历实施方式也可以基于地貌骨架。汇流的两个其它实施方式是基于拓扑学骨架。首先，可以计数有多少个不同的面在任何给定的面上面或下面。如果在另一个特定面下面的一个面终止，则新的面必需存在于原始(未终止的)面的下面。对有多少个不同的面在给定面的上面和下面的计数因而接近于针对它的终端的数目。由于其性质，大的面趋于具有高的计数。然而，通过面尺寸归一化计数趋于指定高的相关计数给最小的面，因为少数目的计数除以小的面积或形成面的小 数目的样本。因此，计数需要被归一化。一种归一化方法包括利用被提高到特定幂例如二分之一的表面积。然而，优选的归一化是基于对在具有特定尺寸的面上遇见一定数目的终端的意义的分析。例如，对于小的面，存在遇见相对于面尺寸大量的终端的机会，因此，该计数可能在统计学上无意义。可是，对于大的面，即使小的终端计数也可能以统计学上有意义的方式大于平均值。面的终端概率可以通过平均化单个相对终端计数、其一些加权形式或通过面的总数目除以这些面的总面积进行评估。一旦预期的终端概率(或密度)被计算，则可应用二项检验来计算与给定尺寸的面的预期终端数目偏差的统计学意义。对于大的面，二项分布通过方便的连续分布被充分逼近，并且这些可被用作可选检验的基础，所述可选检验计算快得多，即，皮尔森卡方检验(Pearson’schi-square test)和G检验。然而,对于少的样本，这些逼近崩溃(break down)，并且没有其它可选的二项检验。源自拓扑学骨架化的另一汇流实施方式是基于与面有关的年代映射体。高汇流的区域具有高密度的流线或流片。年代映射体中的等高线或等值面类似于源自反射切线的流线或片。因此，高密度的等值面表示高汇流的区域。对等值面的构建和其密度的评估的优选方式是计算年代映射体的梯度或垂直导数。具有高梯度或高导数的区域具有高汇流。汇流的任何实施方式均可以进一步通过应用长距离修勻器(smoother),例如线卷积积分来修匀。I. e闭合(closure)和溢出点另一组地质属性——公开于美国专利申请公开号2010/0149917，“Method forGeophysical and Geological Interpretation of Seismic Volumes in Depth, Time, and
Age”，Imhof等中-是包体(enclosure)和溢出点。在传统的解释实践中,层位被分析，
以确定闭合的轮廓(contour)的存在，所述闭合的轮廓环绕地形高点，形成可以圈闭烃类的闭合。注意，闭合可以是暗示性的，例如其中轮廓线在断层处或靠近盐丘处终止，并且被暗示性地假定遵循这些界限。对于每一个闭合的轮廓，包体规定其中包含的区域。因此，在面上的任何位置均可以被分析，以确定它是否位于环绕高点的闭合的轮廓上，如果是如此，则被环绕的区域可以被计算。对于每一个高点，其最大包体规定最大程度的潜在烃圈闭。
可以针对任意单个层位或针对任意一组层位，例如通过传统手段或通过自动骨架化映射的面，测定闭合和包体。通过从年代体(age volume),或者优选地从深度体针对数据集中的每一个面评估包体，可以计算包围体，如果样本不是环绕高点的、闭合的深度轮廓的部分，则所述包围体分配给每一个这样的样本零值。否则，环绕的区域的面积被分配给样本位置。利用年代体，可以通过选择年代、或位置或一些其它标准构建面。该面可以针对高点，环绕高点的轮廓及其面积进行分析。然而，优选使用深度体。深度体中的每一水平切面表示处于一定年代的面。事实上，深度体只不过是一堆以年代排序的面。对于每一切面(或一定年代的面)，其值对应于深度，因而，每一切面构成深度映射。每一切面均可针对高点、闭合的轮廓和闭合的区域进行分析，这允许计算年代域中的包围体。如果期望，整个包围体或其部分，例如切面可以容易地变换成深度域。闭合或包围体的延伸是溢出点和溢出点体。溢出点是在最大闭合轮廓附近的位·置，其中轮廓被破坏，因而潜在的圈闭泄露(图I)。对所有溢出点的测定允许产生溢出点体和检查不同的潜在圈闭如何溢出和彼此流入。溢出点的位置和数目可被用于赌注远景构造或用于引导对其中烃类从中泄露并聚集到其中的区域的研究。尽管通过提取等年代(isoage)面和对其检查可以从年代体鉴定溢出点，但溢出点和溢出点体优选从深度体产生，所述深度体对应于不同年代的深度映射堆(即，不同层位的深度映射堆)。如果需要，从深度体测定的溢出点和溢出点体可通过使用深度体作为查找表或通过插值被容易地变换成深度域。l.f 相留数(Residue)可以由瞬时相异常形成另一地质属性，其指出地震子波分裂(split)的位置(PCT 专利申请公开 W02009/137150A1, “Method For Geophysical AndStratigraphicInterpretation Using Waveform Anomalies”, Imhof)。瞬时相异常属性是离散属性，其指示子波分裂是否发生在样本位置处。用矩形波过滤器的传统卷积允许计算子波分裂的局部密度。例如用线积分卷积的长距离结构导向的修匀允许计算与潜在的地震组构一致的子波分裂的密度。优选地，子波分裂环(loop)和线(string)首先基于其形状和方位被分类，以抑制由噪音或假象引起的子波分裂。相留数发生在波形分裂的位置，其暗示另外的反射事件出现，或者暗示一个反射合并到另一个反射。换言之，终端发生了。终端常常被分类为上超、下超、顶超和削截(图9)，这取决于反射之间的关系。因此，相留数可以被分类为上超、下超、顶超和削截。对于相留数可选的是，还可以从地震骨架检测并分类终端。除非使用偏角信息，否则不能区分下超和上超并且可能需要分在一起。在没有偏角信息的情况下，不能区分顶超和削截并且需要分在一起。终端及其分类可以被布置在面或汇流流线上。优选地，终端及其分类通过修匀被转换成密度或地质属性。基于终端的特别有用的地质属性是终端密度诸如上超/下超密度、削截/顶超密度，或上超/下超密度和削截/顶超密度之间的差异，其指示哪一种行为占优势。l.g 纹理其它地质属性是基于地震纹理。美国专利号6，438，493B1，“Method for seismicfacies interpretation using textural analysis and neural networks”，West 和May，公开了一种方法，用于根据纹理属性鉴定地震相，所述纹理属性通过灰度共生矩阵(gray-level co-occurrence matrix) (GLCM)被计算。源自GLCM的局部纹理属性包括,但不限于纹理一致性、惯性(也称为要素差异矩(moment)或对比(contrast))、熵和能(也称为一致性)。沿着组构利用长距离修匀，可将这些属性转换成区域地质属性，使纹理与地震相相关。基于GLCM的纹理属性可进一步被广义化成三维纹理属性(例如，美国专利号6,226,596B1, “Method for analyzing and classifying three dimensional seismicinformation，，, Gao)。另一纹理属性(PCT专利申请公开 W02009/011735 “Geologic Features FromCurvelet Based Seismic Attributes”，Neelamani 和 Converse)通过米用数据的小曲线变换从地震数据鉴定地层特征。由该小曲线表示，选择的地球物理数据属性及其相互依赖性被提取并用于鉴定地质特征。沿着组构利用长距离修匀，可将这些属性转换成区域地质属性，其涉及纹理和地震相。再一纹理属性是基于窗口化的多维傅里叶度量(multidimensional Fouriermeasures) (PCT 专利申请公开 W02010/053618, “Method for SeismicInterpretation Using Seismic Texture Attributes”，Imhof)。一个具体的这种度量是规律性，测量带状的(banded)(或有规律的)地震数据看起来如何的纹理属性。具有纯净、通透反射的区域显示高度的规律性，而具有紊乱反射的噪音区域显示低规律性。对于小的分析窗口，规律性是间断性度量。对于大的分析窗口，规律性充当混沌度量。因此，规律性也可被称为多迹间断性混沌。沿着组构利用长距离修匀，可将这些窗口化的多维傅里叶属性转换成区域地质属性，其涉及纹理和地震相。例如，通透区域反射一通常被称为铁路轨迹(railroadtrack)—常常与海侵面(或溢流面)相应，其频繁地形成封闭层。因此，规律性的长距离修匀(或与其相反，混沌)产生与烃封闭层有关的地质属性，并且因为页岩均形成封闭层并充当源岩、烃源岩。Randen 和 Sonneland( “Atlas of 3D Seismic Attributes”，在MathematicalMethods and Modeling in Hydrocarbon Exploration and Production 中Iske andRanden(editors)，Springer, 23-46页(2005))提出对表征地震纹理或地震-地层特征的其它三维地震属性的综述。l.h WPCA 异常与本发明有关的另一组属性公开在PCT专利申请公开W02010/056424 “WindowedStatistical Analysis for Anomaly Detection in Geophysical Datasets，，, Kumaran 等中。这些属性突出一个或多个数据集中与其它位置相比地震数据在统计学上异常的位置。此外，这些属性产生模型的向量基础，其范围从最常见的到最异常的。将数据投影到这些模型的线性组合上强调一些模型同时抑制其它模型。例如，普遍存在的地震条带可以被抑制。沿着组构利用长距离修匀，过滤后的地震可以被转换成区域地质属性。另一应用是将数据局部分解成这些模型并确定占据支配地位的那个模型，其影响根据周围样本的模型对地震样本的分类。沿着组构利用长距离中值过滤或模式过滤，局部类别可以被转换成区域类别。Li其它属性另一种属性是基于多偏移(multi-offset)地震数据的反演和观察到的振幅对比偏移的行为，允许预测孔隙度和粘土含量和允许对岩相进行分类(例如，美国专利号 7,424,367B2 “Method for predicting lithology and porosity from seismicreflection data，，, Saltzer等；“Seismic Rock-Property Inversion and岩相Predictionat Erha Field, Nigeria，，，Xu 等，Nigerian Association of Petroleum Explorationists (NAPE) 2008conference proceedings ;“Lithofacies Prediction in Deep Water WaterReservoirs，，，Oppert 等，Society ofExploration Geophysicists, ExpandedAbstracts, I708-1711，(2006))。沿着组构利用长距离修匀，可将这些属性转换成区域地质属性，其涉及岩相、粘土含量和孔隙度。欧洲专利号EP1110103Bl“MethodOf Seismic Signal Processing”, Meldahl 等,公开了一种产生属性的方法，其允许检测潜在的气烟道(chimney)。气烟道是来自气渗漏的地震响应的垂直干扰，所述气渗漏使地震数据恶化，这是由于在获得和处理中的局限性。所得烟道立方体(cube)属性突出地震信号的垂直干扰，这常常与气烟道有关。它揭示烃历史和流体流动的信息。换言之，烟道立方体可以揭示烃类起源自哪里，它们如何被运移到远景构造中以及它们如何从该远景构造中溢出。因此，烟道立方体可以被视为新的间接烃指示工具。 此外，I.oseih等，(“Hydrocarbon leakage interpreted on seismic data”, Marineand Petroleum Geology26 (7)，1304-1319，(2009))提供地震烃泄露解释的综述。最后，地震属性可以与地下特性有关，可用于模拟运输现象诸如热流和影响烃熟化的温度，或影响烃运移的渗透性和流体流动(PCT专利申请公开W02009/137228A2, “Transport Property Data Calculated From Derivative SeismicRock Property Data For Transport Modeling，，，Oppert 等)。I. j组合属性/远景带属性在本发明方法的一些实施方式中，地震数据针对一个或多个具体远景带而不是更一般的烃体系的要素被检验。实例可以包括背斜远景带(图10)、法线-断层远景带(图11)、盐丘侧翼远景带(图12)、鞋带状通道远景带，或其它诸如图13到16中所示出的那些。在这些情况下，需要属性、检测器或工作流程来区分具体特征，诸如断层、盐或通道。PCT 专利申请公开 W02009/082545 “Detection of Features in SeismicImages”，Kumaran和Wang，公开了一种这样的方法，其用于检测地震数据中的通道或断层。对于断层检测，在修匀的地震图像中鉴别边缘，并且边缘强度在多个方向被积分，例如利用拉东变换，以检测断层线的存在和方位。对于通道，边缘被检测并转换成修匀的曲线，以鉴别通道边缘。然后，检验平行的修匀曲线组，以发现对应于左和右通道边缘的成对曲线，从而限定通道。美国专利号7，203，342B2, “Image Feature Extraction，，，Pedersen,公开了另一种这样的方法，其被设计来从地震属性数据中提取断层，尽管它也可以被用于其它线(line)或面提取问题。该方法是基于蚂蚁追踪(anttracking),其中数字蚁群(numericalants)或媒介(agent)在间断性或边缘检测体上漫游，缓慢地连接上排列的边缘的附近块(piece)。最后的实例是检测盐体，其利用在PCT专利申请公开W02010/053618 “Methodfor Seismic Interpretation Using Seismic Texture Attributes，，，Imhof ; 或Randen 和 Sonneland, “Atlas of3D Seismic Attributes” (在 Mathematical Methodsand ModelinginHydrocarbon Exploration and Production 中，Iske 和 Randen(编著)，Springer, 23-46页(2005))中公开的规律性或混沌属性。这些属性突出地震数据中看起来混杂有不规则纹理的区域，其常常与盐体对应。这些仅是鉴别特定特征的三个实例方法。已经公开了许多其它方法，它们对于本领域的从业者来说是悉知的。划分为了检测多个近源烃体系要素或多个远景带要素的同时存在，数据体可任选地被划分成至少一个用于分析的段，可以结合未被分析的背景段。典型的划分模式是分析单个体素、小的区块、组构-排列的区块、层或邻近的体素组。划分不需要是相互排除的。单个划分可以重叠。最简单的划分是体素-对-体素分析，但结果可能是不调和的，因为很少是在相同体素中识别的所有要素，并且不可能遇见所有要素都存在的邻近体素的延伸区域。诸如修 匀属性或扩充检测的要素到周围的体素中的技术可用于本发明方法的步骤5中(图2)，评价，以产生较大的邻近的远景构造。划分成体素是优选的划分方法，因为数据常常已经被表示为体素。由于这个理由，划分步骤被认为是任选的。此外，在评价步骤过程中，划分(以及因而体素)可以被聚集成较大的区域。因此，为了教导本发明方法的目的，如果没有明确进行划分，那么每一个体素均被认为形成其自身的划分。另一划分方案包括将数据体打破成规则的卡迪尔区块或程序块(brick)，例如大小为20X20X20的样本，其中预期是一些或所有要素存在于预期的程序块内。卡迪尔程序块将通过岩层和组构切割。规则卡迪尔程序块或区块的替代方案是将区块排列成组构。在该方案中，在划分之间将存在尺寸和形状的差异，因为它们与组构一致。产生这种划分的具体方式是通过使用组构本身来限定总的层结构。具体地，例如，对由地震骨架化产生的面块的选择可以被垂直扩充，以产生体或段。用于划分的另一优选的方法是基于一个或多个产生的地震属性，例如，在步骤2中。分区通过属性的阈值化接下来进行连通域分析(connected component analysis)而产生或类似于产生嵌入背景中的邻近区域。该方法也被称为(多体)种子检测。控制划分的优选属性是凸起(saliency)，突出一个或多个数据集中相比于其它位置地震数据统计学上异常的位置的属性。凸起属性公开在PCT专利申请公开W02010/056424 “WindowedStatistical Analysis for Anomaly Detection in Geophysical Datasets，，，Kumaran等中，凸起是突出数据统计学异常的属性的总称。Kumaran等公开了“逆协方差(inversecovariance) ”、“ (WPCA)留数”或“ (WPCA)异常”，这些是凸起的具体实施方式
的术语。可以利用Kumaran的方法计算另外的凸起属性。凸起可以被阈值化以仅在邻近的异常区域上进行分析。利用属性控制划分的可选方法是首先进行不同的划分方法，例如卡迪尔区块或组构排列的体，然后仅保存这样的分区其中针对每一分区计算的属性超过绝对值或相对值或满足一些规定的条件，例如，其中最大凸起处于上五分位数(top quintile)处或其中平均规律性在O. 3和O. 6之间的范围内。划分方法的最后实例包括将数据体再分成小的区域，其通过属性(一个或多个)，例如凸起进行分级。对于本发明的许多实施方式，评估不同分区或单个体素的远景性的顺序并不重要。该最后的实例显示属性诸如凸起(或尺寸、另一属性或其它属性的组合)用于确定分析分区的顺序的用途。利用这种优先性，可能不必要分析所有划分的区域。本发明方法的该实施方式开始用最凸起的区域评估烃体系(步骤5)，并继续用逐渐较少凸起的区域进行评估，直到已经分析了规定数目的区域、已经发现规定数目的远景构造或已经超过允许分析的规定时间。要素评分评分可以包括归一化、调节、组合或定标(scaling)。对于一些远景带或烃体系要素，该步骤对于每一个可利用要素采用至少一个属性，以形成每一分区的评分(或者在最简单的情况下，每一体素的评分)，其表示特定要素包含在特定分区中的似然性或预期。评分被认为是可选步骤，因为它可以作为图2所示流程图中显示的明确的、独立步骤4进行。然而，在本发明的一些实施方式中，它可以属性形成(步骤2)或远景构造评估(步骤5)组合进行。在本发明方法的最简单的形式中，属性被直接用作评分，以指示特定要素是否存在于任何给定的位置。然而，不同属性的值可以具有不同的范围，并且，将它们的值归一化，·例如在O和I或O到255之间常常有利于属性之间更容易的比较或有利于促进概率解释。归一化仅仅是属性值的线性转换。归一化的替代方案是校准或应用属性值的非线性转换，其或者由数据，例如直方图均衡化或直方图变换驱动，或者由模型驱动，所述模型表达属性如何与给定要素存在的似然性相关。这种模型可以基于理论、测量结果、现有知识或背景知识、经验或直觉。地震属性是在地震数据上进行的观察、测量或计算。它们可以涉及烃体系要素或远景带要素，但是，它们常常并不实际测量这些要素。此外，它们可能不是唯一的。不同属性可以涉及相同的要素。相同属性可以涉及多个要素，同时也受噪音、采集、数据处理以及用于产生属性的算法和参数影响。相同类型的属性可以用不同的算法计算。属性的组合允许形成一个评分，其指示要素在给定位置存在的似然性或预期。如果对于给定要素，没有计算限定评分的直接属性，则需要利用替代属性或评分、启发、概念或现有预期(如常数值O. 1，其仅仅表示要素的10%机会)间接限定评分。例如，储层可以通过砂-页岩比例评估和/或孔隙度评估间接表示。如果没有这种指示属性存在，那么可能需要应用例如源自层序地层概念的概念替代。这种替代可以是覆盖不整合(层序边界)的低水位期盆地底部扇形地的组合，如图17所示。低水位期盆地底部扇形地常常显示良好的孔隙度和渗透性，并因而具有形成储层要素的潜力。提高相对海平面将通过低水位期楔形的下端(toe)覆盖该扇形地，所述低水位期楔可能由被海侵面封住的、更细的渗透性较小的物质组成。进一步提高相对海平面将使低水位期埋藏在海侵和高水位期系统区块下面。在远端位置，扇形地可以看起来被夹在“不整合”和“海侵”面之间。因此，不整合和海侵面之间的岩石体一位于所述不整合的深侧上一可以是低水位期扇形地，并因而可以被分配高的储层分数。评分分配的另一方面是段或分区的定标。在每一体素位置限定许多属性。其它的仅在面上或少数选择位置处可得。具体的分区可以对应于一个体素或邻近的体素组。为了分配分数给每一个分区(或者在最简单的情况下，每一个体素)，属性可能需要被插入或减少，以获得每一分区的单一分数。对于包含多个样本的分区，可以通过计算平均数、中值或模式；应用表决程序；或选择最小、最大或随机样本值而找到代表分数。用置信度增大分数是有优势的。分数代表在给定位置遇到特定要素的似然性。然而，分数并不表示该似然性有多可信。地震属性预测高孔隙度的位置可以接受与具有低水位期扇形地的位置相同的分数，如通过不整合和海侵面的接近性所暗示的。但是在第一位置，存在对孔隙度相对直接的测量，因而对储层存在也存在相对直接的测量。在第二位置，由概念模型假定储层的存在。直接测量产生比基于概念模型的预测更高的置信度。该差异表达在与分数有关的置信值中。远景构造评价利用远景带或烃体系要素的分数允许评价分区(或体素)远景构造(步骤5)并识别其缺点，即，弱表达、缺乏或未分辩的要素。许多公开的地质属性表示趋势，因而具有相对低的分辨率。在本发明系统中的一些实施方式中，对于要素空间排列，例如封闭层在储层上面，的要求是放松的，并用针对所 有要素在单一分区内的存在的检验进行替代。针对所有要素的存在的该检验被称为门控逻辑检验(gatedlogic test)，因为分区必需通过针对所有要素的检验，以被视为是远景的。因此，组合要素的门控逻辑方法是一种“不全则无的(all or nothing)”方法。满足要素常常意味着超过分数阈值，该分数阈值通常由解释者指定或被编码在系统中。不满足的要素可以被标记用于进一步检验，因为它们构成分区的弱连接。门控逻辑是远景构造评价的优选实施方式。更渐进的远景性检验是表决(vote)，其仅仅计数有多少要素存在于给定的体素或分区中。存在的要素越多，分区越被视为是远景的。根据其评分与由解释者指定的或在系统中编码的阈值的比较，要素或者存在或者不存在。不存在或弱表达的要素可以被标记。表决是远景构造评价的另一优选的实施方式。表决的变化是加权表决，其中一些要素被给予更高的权重或比其它更多的表决。可以应用加权表决，例如强调可用更高置信度检测的要素。门控逻辑与表决组合是可能的，例如通过要求指定数目的要素存在而不是要求所有要素存在。另一实例是一些指定要素的必需存在和其它要素的期望存在。远景构造评价可以包括步骤4的方面，评分。实例是涉及相同要素(一个或多个)的不同属性的组合，例如用不同算法或参数化法计算的属性。不同于将属性首先组合成单个要素的分数然后根据分数评价远景构造，例如通过表决或计数有多少属性超过用户或系统规定的阈值从不同属性可以直接评价分区远景性。通过应用多值逻辑，例如源自模糊集合论的模糊逻辑可以避免硬阈值，以处理近似而不是精确的推理。门控逻辑是明确的和二元的，其中隶属值O和I表示低于阈值(要素不存在)和高于阈值(要素存在)。相反地，模糊逻辑分数具有范围从O到I的隶属值，并且代表陈述的真值度。真值度和概率均可以在O和I的范围内，因而可以看起来是类似的。然而，它们在概念上不同。真值代表模糊限定的集合中的隶属关系，而不是概率理论中某一事件或条件的似然性。例如地层层，其包括70%粉砂和30%页岩。我们可以考虑两个概念储层和封闭层。它们每一个的意思均可以通过一定的模糊集合表示。可以限定地层为O. 7储层和O. 3封闭层。注意，储层的概念可以是主观性的，因而会取决于观察者或设计者。另一设计者可以等同地、很好地建立集合隶属函数，其中如果粉砂部分超过50%，地层可以被认为是储层。模糊逻辑利用真值度作为模糊现象的数学模型，而概率是随机性的数学模型。概率设置会首先限定粉砂部分的标量变量，其次限定条件分布，描述某人在特定粉砂部分的情况下将称地层为储层的概率。注意，条件化可以通过使特定的观测器随机为层选择标记、在确定性观测器中分布或两者来实现。因此，概率与模糊度没有任何共同点，这些仅仅是表面上看起来类似的不同概念，因为它们使用O和I之间相同的实数区间。尽管如此，仍存在混乱，因为随机变量的特性类似于二值逻辑态和定理诸如德摩根定律的特性，其通过否定就彼此而言使逻辑操作子“和”和“或”关联，具有双重应用性。模糊逻辑的替代方案是 贝叶斯逻辑，其是基于贝叶斯概率理论，该理论能够用不确定的陈述推理。为了评估假设或构形的概率，人或机器解释者规定一些先验概率，其然后按照新的相关数据被更新。贝叶斯解释提供标准程序和公式集合来进行该计算。用于汇集对特定构形的要求和测定其似然性的一个方法是通过使用贝叶斯信念网络(BBN)。贝叶斯网络是一种基于概率图的模型，其通过无回路有向图(directed acyclic graph) (DAG)表示一组随机变量及其条件依赖性。该图表示随机变量之间的条件独立结构，代表不同要素及其空间排列。例如，贝叶斯网络可以代表远景带要素和远景带类型之间的概率关系。在给出某些远景带要素的情况下，网络可被用于计算各种远景带类型存在的似然性。在形式上，贝叶斯网络是无回路有向图，其节点表示贝叶斯意义上的随机变量它们可以是可观察的量、隐性变量、未知参数或假定。边缘代表条件依赖性，未被连接的节点代表条件上彼此独立的变量。图18代表整个烃体系的实例网络图。储层和封闭层的分数可以通过分析其并置而被组合，所述并置指示储层是否被封闭层封住。在给出局部几何形状的情况下，储层/封闭层系统被分配圈闭分数。优选地，圈闭填充有烃类，因而，圈闭分数然后与源存在的分数和从源引导到圈闭的运移路径的至少可能的指示的分数组合。理想地，在地震数据中也有对烃类的直接指示(DHI)。因此，圈闭、源、运移和DHI的分数均被组合，以产生聚集分数。所有聚集分数现在可以通过大小和置信度被加权，以赌注和分级潜在的烃目标。甚至在应用代表趋势的地质属性时，分区可能太小，以至于不能包含所有要素，例如当划分成单个体素时。检验小的分区内或体素之间的共存或并列的两种方法是相对空间移位和影响区域的扩充。在空间上移位分数将对给定空间关系的检验转换成对并列的检验。例如，封闭层分数可以相对于储层分数被向下移位，以检查封闭层在储层上的存在。不应用这种移位，则需要对封闭层在储层上的检验。通过应用这种移位，可以检查封闭层和储层的并列存在，例如用门控逻辑。与利用单一移位不同，针对一定范围的移位，优选较小的移位，评价两个分数的共存可能是有优势的；例如，分配较高的置信度给较小的移位，以指示更接近需要的要素或作为距离的函数减小远景性。横向位移可用于检测远景带要素诸如邻近储层的断层。不同于应用垂直或横向位移，可沿着组构或垂直于组构应用移位。此外，移位还可以以任意方向被应用。如果需要检查许多不同的范围和方位或方向，移位可能是低效率的。移位的优选替代方案是针对原分数或阈值化分数扩充影响区域，例如通过卷积或地貌扩张。在前者的情况下，分数被抹掉或沿着指定的方向扩充，例如通过锥形定向滤器进行卷积，所述锥形定向滤器随着范围增加而衰变，编码置信度的降低。在通过卷积的影响区域扩充的一个具体实施方式
中，卷积的结果被加入到原始分数中，以形成更新的分数，用于远景构造评价。扩充方向是从组构由解释者规定的或决定的。
原始分数或阈值化分数可以通过定向地貌扩张被扩充，所述定向地貌扩张拖动高分数沿着指定的方向进入到具有较低分数的区域。方向可以由解释者规定，被编码在该系统中或源自组构。地貌操作可应用于原始分数或阈值化分数以及任意的划分。图19显示本发明方法对四个分区或潜在目标区域的示意性应用。假定源分数和运移分数均不能由给定的数据确定，也没有观察到直接的烃指示并且剩余要素是独立的，则每一目标的远景性可以通过使置信度与储层、封闭层和圈闭分数的大小和平均值相乘来计算。表I示出该评分方法的实例。
权利要求
1.计算机执行的方法，用于针对烃体系或特定远景带的存在分析由代表地下区域的地震数据的体素组成的体，所述方法包括 划分所述地震数据体，以形成多个段；和 至少部分根据每一段中地震数据体素的远景性分数，针对烃体系或所述特定远景带的存在对所述多个段分级； 其中所述远景性分数是基于对至少两个属性的计算，所述至少两个属性涉及烃体系或所述特定远景带的不同要素。
2.权利要求I所述的方法，其中所述划分包括根据针对每一体素计算的凸起属性或远景性分数组合邻近的单个体素。
3.权利要求I所述的方法，其中所述要素选自储层、封闭层、圈闭、源、填充物、上覆岩层、熟化、运移、聚集和时间。
4.权利要求I所述的方法，其中所述要素选自储层、封闭层、圈闭和填充物。
5.权利要求I所述的方法，其中，为了计算每一体素的所述远景性分数，选择至少两个要素，并且对于每一要素，选择至少一个地震属性，然后，根据选择的至少一个地震属性计算每一选择的要素的要素分数，并且通过组合所述要素分数形成所述远景性分数。
6.权利要求I所述的方法，其中修匀被应用于所述至少两个属性，从而产生地质属性，并且所述修匀被应用于至少十个地震数据体素的距离上。
7.权利要求6所述的方法，进一步包括从代表所述地下区域的所述地震数据确定层结构，所述层结构被称为组构，并且其中产生所述地质属性的所述修匀包括沿着所述组构积分或平均化所述选择的属性，以防止与其它岩层混合。
8.权利要求7所述的方法，其中沿着所述组构的所述修匀通过非线性各向异性扩散过滤或通过线积分卷积进行。
9.权利要求7所述的方法，其中所述组构由地震骨架化确定。
10.权利要求I所述的方法，其中所述至少两个属性包括黑塞属性；等高线属性；收敛；以及汇流中的至少一个。
11.权利要求I所述的方法，其中所述至少两个属性包括包体和溢出点；规律性；相留数；纹理；窗口化的主分量分析异常；远景带属性；终端；终端密度；扩散映射；以及其组合中的至少一个。
12.权利要求I所述的方法，其中所述段彼此不重叠。
13.权利要求I所述的方法，其中根据通过组合所述段中每一体素的所述远景性分数而产生的每一段的单个分数，对所述段进行分级。
14.权利要求I所述的方法，其中所述地震数据体的每一体素是单独的段。
15.权利要求I所述的方法，其中针对不同的属性计算的属性值被归一化成常用尺度。
16.权利要求I所述的方法，其中测定所述远景性分数利用阈值、模糊逻辑和贝叶斯逻辑中的一个。
17.权利要求I所述的方法，进一步包括评估至少一个远景性分数的置信度值。
18.权利要求6所述的方法，其中所述至少两个属性中的一个涉及烃圈闭的存在，并且所述烃圈闭的存在通过空间关联所述地质属性并与烃圈闭构形的目录进行比较进行预测。
19.权利要求2所述的方法，其中所述凸起属性利用窗口化的主分量分析或应用于至少一个地震属性的扩散映射进行计算。
20.计算机程序产品，包括计算机可用介质，所述计算机可用介质上嵌入计算机可读程序代码，所述计算机可读程序代码适于被执行，以实施针对烃体系或特定远景带的存在分析代表地下区域的地震数据的方法，包括 划分地震数据体，以形成多个段；和 至少部分根据每一段中地震数据的体素的远景性分数，针对烃体系或所述特定远景带的存在对所述多个段进行分级； 其中所述远景性分数是基于对涉及烃体系或所述特定远景带的不同要素的至少两个属性的计算。
21.生产烃类的方法，包括 进行地下区域的地震调查； 从调查获得对地震数据的远景性分析，其中所述分析通过权利要求I所述的方法进行，该方法通过引用被并入此处； 至少部分根据远景性分析向所述地下区域中钻井，并从所述井生产烃类。
全文摘要
针对烃体系或特定远景带的存在分析代表地下区域的地震数据的方法。计算地震属性，所述属性被选择以与烃体系的典型要素，即储层、封闭层、圈闭、源、熟化和运移相关。优选地，沿着地下区域的结构组构(1)计算所述属性，并且在至少上十或上百的数据体素上修匀所述属性。所得地质属性(2)用于分析烃体系要素的数据和/或识别具体远景带，并且用于在远景性预测(5)中根据大小、性质和置信度分级和评注划分的数据体区域(3)。可以生成烃圈闭构形的目录(8)，并用于鉴定烃圈闭的潜在存在和/或帮助评分(4)并分级划分的区域为烃远景构造。
文档编号G01V1/40GK102918423SQ201180026539
公开日2013年2月6日 申请日期2011年4月22日 优先权日2010年5月28日
发明者M·G·英霍夫, P·迪米特罗夫, K·罗贝尔, K·库马兰, M·J·特雷尔, S·胡斯瑙德 申请人:埃克森美孚上游研究公司