BR112012028653B1

BR112012028653B1 - método para análise sísmica de sistema de hidrocarbonetos

Info

Publication number: BR112012028653B1
Application number: BR112012028653-6A
Authority: BR
Inventors: Matthias G. Imhof; Pavel Dimitrov; Kelly Wrobel; Krishnan Kumaran; Martin J. Terrell; Stefan Hussenoeder
Original assignee: Exxonmobil Upstream Research Company
Priority date: 2010-05-28
Filing date: 2011-04-22
Publication date: 2020-11-10
Also published as: AU2011258764A1; WO2011149609A1; MY162927A; US20130064040A1; CN102918423B; RU2012157782A; CA2799504C; CN102918423A; AU2011258764B2; EP2577356B1; CA2799504A1; US9194968B2; EP2577356A4; BR112012028653A2; EP2577356A1; RU2573166C2

Abstract

MÉTODO PARA ANÁLISE SÍSMICA DE SISTEMA DE HIDROCARBONETOS Método para analisar dados sísmicos representando uma região de subsuperfície quanto à presença de um sistema de hidrocarbonetos ou um play (um modelo conceptual de um estilo de acumulação de hidrocarbonetos com frequência usado para desenvolver prospecção em uma bacia) particular. Atributos sísmicos são computados, os atributos sendo selecionados em relação aos elementos clássicos de um sistema de hidrocarbonetos, isto é, reservatório, selagem, trape, fonte, maturação, e migração. Preferivelmente, os atributos são computados ao longo de tecidos estruturas (1) da região de subsuperfície, e são atenuados através de pelo menos dezenas ou centenas de voxels de dados. Os atributos geológicos resultantes (2) são usados para analisar os dados quanto aos elementos do sistema de hidrocarbonetos e/ou reconhecimento de plays específicos, e para classificar e anotar regiões divididas (3) do volume de dados com base no tamanho, qualidade, e confiança na predição de prospectividade (5). Um catálogo (80) de configurações de trape de hidrocarbonetos pode ser criado e usado para identificar a presença potencial de trapes de hidrocarbonetos e/ou ajudar no score (4) e/ou classificação das regiões divididas como prospectos de hidrocarboneto.

Description

[0001] Este pedido reivindica o benefício do Pedido de Patente Provisória U.S. 61/349.534, depositado em 28 de Maio de 2010, intitulado SYSTEM FOR SEISMIC HYDROCARBON SYSTEM ANALYSIS, cuja totalidade é incorporada aqui por referência.

CAMPO DA INVENÇÃO

[0002] Esta invenção refere-se genericamente ao campo de prospecção geofísica e, mais particularmente, à interpretação de dados sísmicos. Especificamente, a descrição descreve um método para detectar e classificar oportunidades potenciais de hidrocarbonetos utilizando dados sísmicos.

FUNDAMENTOS DA INVENÇÃO

[0003] Um sistema de hidrocarboneto ativo é definido pela presença de uma formação de reservatório porosa que provê espaço de armazenagem para hidrocarbonetos, uma selagem que evita que hidrocarbonetos escapem do reservatório, uma boa geometria de aprisionamento, e uma formação de fonte que contém uma alta porcentagem de material biogênico. Sob a influência de elevada temperatura e aumentada pressão, o material biogênico é maturado (ou cozido) para formar hidrocarbonetos incluindo gás, óleo bruto, asfaltos e alcatrão. Acionado por flutuação e pressões diferenciais, os hidrocarbonetos migram e uma fração destes hidrocarbonetos acumula-se em trapes formados por arranjos geométricos acidentais de formações de reservatório (isto é, geometrias de aprisionamento) e selagens. Os trapes têm um volume finito, entretanto, e podem derramar ou vazar algum dos hidrocarbonetos acumulados, cuja parte pode então acumular-se em outros tanques.

[0004] Imagens sísmicas da subsuperfície permitem que intérpretes identifiquem alguns trapes potenciais baseados em geometrias práticas e sugestivas. Às vezes, os dados sísmicos podem prover uma indicação direta quanto à presença de hidrocarbonetos. As práticas de interpretação padrão, entretanto, são de intensa mão-de-obra e com frequência focaliza-se em áreas em que o intérprete coleta alguma indicação de prospectividade. Muitas oportunidades, portanto, permanecem não detectadas em razão das indicações serem muito sutis ou ocultas, por exemplo, por ruído sísmico. Mesmo se sinais de prospectividade forem observados, eles podem não ser examinados quando na presença de mais óbvias oportunidades, ou quando o intérprete está limitado pelas restrições de tempo. Assim, algumas acumulações de hidrocarbonetos são descobertas tarde ou permanecem não descobertas.

[0005] Experiências publicadas para resolver problemas similares incluem as seguintes:

[0006] “Method for Seismic Interpretation Using Seismic Texture Attributes” (Publicação do Pedido de Patente PCT WO 2010/053618), de Imhof, descreve urn método para a computação de atributos de textura que pode ser usado para score e segmentação de dados sísmicos baseados em sua aparência local. A textura pode ser usada para definir fácies sísmicas.

[0007] “Windowed Statistical Analysis for Anomaly Detection in Geophysical Datasets” (Publicação do Pedido de Patente PCT WO 2010/056424), de Kumaran et al., descreve um método para examinar dados sísmicos em regiões que são estatisticamente anômalas no contexto dos dados e, assim, servem para realçar estatisticamente áreas não usuais ou salientes.

[0008] “Method For Geophysical and Geological Interpretation of Seismic Volumes In Depth, Time, and Age” (Publicação do Pedido de Patente U.S. No. 2010/0149917), de Imhof et al., descreve um método para a transformação de dados sísmicos a partir de domínios geofísicos de profundidade ou tempo de deslocamento de duas direções para um domínio de idade geológica em que todas as reflexões sísmicas são aproximadamente horizontais e comparáveis a seu estado no tempo geológico de sua formação. Os atributos sísmicos formados a partir deste domínio de idade podem melhorar a definição dos elementos do sistema de hidrocarboneto.

[0009] A Patente Européia EP1110103 B1 (“Method Of Seismic Signal Processing”), Meldahl et al., descreve um método para encontrar áreas em dados sísmicos que são similares às áreas especificadas pelo intérprete. Além disso, ela descreve o uso deste método para encontrar indicações sísmicas de hidrocarbonetos de percolação.

[0010] A Patente U.S. No. 6.226.596 B1 (“Method for analyzing and classifying three dimensional seismic information”), de Gao, descreve um método para gerar atributos de textura sísmica que podem ser usados para segmentação, score, ou definição de fácies sísmicas.

[0011] A Patente U.S. No. 6.438.493 B1 (“Method for seismic facies interpretation using textural analysis and neural networks”) de West e May, descreve um método para gerar atributos de textura sísmica que são usados em um score supervisionado para designar atributos de fácies sísmicas.

[0012] A Patente U.S. No. 6.516.274 B2 (“Method for imaging discontinuities in seismic data using dip-steering”) de Cheng et al., descreve um método para a detecção de bordas ou descontinuidades em dados sísmicos que com frequência indicam rompimento e deslocamento de falhas na subsuperfície da camada regular.

[0013] A Patente U.S. No. 6.725.174 B2 (“Edge-preserving enhancement of seismic images by nonlinear anisotropic diffusion”) de Bouts et al., descreve urn método de processamento de dados sísmicos que acentua os dados sísmicos por remoção de ruído desconexo, enquanto preservando descontinuidades pronunciadas correspondendo à falhas ou transições estratigráficas rápidas, tais como canais entrincheirados.

[0014] A Patente U.S. No. 6.850.845 B2 (“System for multi-dimensional data analysis”), de Stark, descreve um método de processamento de dados sísmicos que permite a computação de um atributo de tempo geológico usado para atributos derivatórios de aplanamento e associados.

[0015] A Patente U.S. No. 6.850.864 B2 (“Method for analyzing dip in seismic data volumes”), de Gillard et al., descreve um método de processamento de dados sísmicos para a estimativa de inclinações de reflexão que definem a estrutura local.

[0016] A Patente U.S. No. 7.203.342 B2 (“Image feature extraction”), de Pedersen, descreve um método de processamento de dados sísmicos para detectar e melhorar falhas ou horizontes em dados relativamente ruidosos.

[0017] A Patente U.S. No. 7.424.367 B2 (“Method for predicting lithology and porosity from seismic reflection data”), de Saltzer et al., descreve um método para predizer litologia e porosidade de rochas de subsuperficie a partir de dados sísmicos e, assim, permite a diferenciação entre rochas reservatório e formações de rochas de selagem.

[0018] A Patente U.S. No. 7.454.292 B2 (“Inverse-vector method for smoothing dips and azimuths”), de Wang et al., descreve um método para a computação robusta de inclinações de horizonte e orientações que definem a estrutura local.

[0019] A Publicação do Pedido de Patente PCT WO 2009/01 1735 (“Geologic Features From Curvelet Based Seismic Attributes”), de Neelamani e Converse, descreve um método para a computação de indicadores de hidrocarboneto ou atributos de textura que podem ser usados para a identificação de aspectos de subsuperficie.

[0020] A Publicação do Pedido de Patente PCT WO 2009/082545 Al (“Detection Of Features In Seismic Images”), de Kumaran e Wang, descreve um método para a detecção de falhas, canais, e aspectos similares nos dados sísmicos.

[0021] A Publicação do Pedido de Patente PCT WO 2009/137150 Al (“Method For Geophysical And Stratigraphic Interpretation Using Waveform Anomalies”), de Imhof, descreve um método de processamento de dados sísmicos para mapear terminações estratigráficas e pinch outs (afinamento da camada de reservatório).

[0022] A Publicação do Pedido de Patente PCT WO 2009/137228 A2 (“Transport Property Data Calculated From Derivative Seismic Rock Property Data For Transport Modeling”), de Oppert et al., descreve um método para avaliar propriedades, tais como fluxo de calor ou permeabilidade de fluido, que afetam elementos do sistema de hidrocarboneto.

[0023] A Publicação do Pedido de Patente PCT W02009/142872 A (“Seismic Horizon Skeletonization”), de Imhof et al., descreve um método automático para extrair urn grande número de horizontes a partir de um conjunto de dados sísmicos. Além disso, ela descreve um fluxo de trabalho de reconhecimento de amplo padrão, que divide um conjunto de dados, analisa as regiões, e classifica-as de acordo com seu potencial de conter hidrocarbonetos.

[0024] “A New Class of Large-scale Attributes for Seismic Stratigraphy”, de Gesbert et al, 71st EAGE Conference & Exhibition, (2009) descreve um conjunto de atributos estratigráficos computados a partir de dados sísmicos bidimensionais que acentuam desconformidades regionais e tendências de afinamento regional e quantifica tendências de fácies sísmicas regionais.

[0025] “Applications of plane-wave destruction filters” de Fomel, Geophysics 67, 1946-1960, (2002), descreve um método para avaliar a inclinação e orientação de horizonte sísmico que definem a estrutura local.

[0026] “Estimating Seismic Heterogeneity with the Structure Tensor” de Imhof, 67a EAGE Conference & Exhibition, (2005), descreve um método para avaliar a inclinação e orientação de horizonte sísmico que definem a estrutura local e atributos de textura sísmica que caracteriza a heterogeneidade local.

[0027] “Flattening without picking” by Lomask et al, Geophysics 71, P13-P20 (2006), descreve um método de processar dados sísmicos para nivelar aproximadamente os dados que perm item a caracterização de alguns elementos do sistema de hidrocarboneto.

[0028] “Hydrocarbon leakage interpreted on seismic data” de Loseth et al, Marine and Petroleum Geology 26, 1304-1319, (2009), descreve métodos de acionamento- intérprete para detecção de filtração de hidrocarbonetos através da subsuperfície.

[0029] “Hydrocarbon Traps, K.T. Biddle and C.C. Wielchowsky, The Petroleum System - From Source to Trap, AAPG Memoir 60, págs. 219-235, (1994), apresenta uma coleção de tipos de trape de hidrocarbonetos.

[0030] “Imaging Vector Fields Using Line Integral Convolution” by Cabral and Leedom, Proceedings of ACM SigGraph 93, 263-270, (1993), descreve um método de visualizar campos vetoriais de linhas de fluxo.

[0031] “Lithofacies Prediction in Deep Water Water Reservoirs” de Oppert et al, Society of Exploration Geophysicists, Expanded Abstracts, 1708-171 1, (2006), descreve um método para avaliar a litologia da subsuperfície empregando dados sísmicos e com-fio.

[0032] Ata da conferência de 2008 sobre “Seismic Rock-Property Inversion and Lithofacies Prediction at Erha Field, Nigeria” de Xu et al, Nigerian Association of Petroleum Explorationists (NAPE), descreve um método para avaliar a litologia da subsuperfície empregando dados sísmicos e com-fio.

[0033] Randen and Sonneland (“Atlas of 3D Seismic Attributes”, in Mathematical Methods and Modeling in Hydrocarbon Exploration and Production, Iske e Randen (editores), Springer, págs. 23-46 (2005)), apresenta um resumo de atributos sísmicos tridimensionais que caracterizam aspectos de textura sísmica ou sismo- estratigráficos.

[0034] O que é necessário é um sistema automatizado que explore um conjunto de dados inteiro para os elementos de um sistema de hidrocarboneto e emita uma lista de prospectos para o intérprete examinar. Preferivelmente, esta lista de alvos potenciais é classificada por volume esperado, presença e qualidade de elementos do sistema de hidrocarboneto, e confiança em sua detecção e identificação. Preferivelmente, a lista de prospectos é também anotada. A presente invenção satisfaz pelo menos estas necessidades.

SUMÁRIO DA INVENÇÃO

[0035] A invenção, em um de seus aspectos, é um método que computa atributos sísmicos de multiescalas e de estrutura tipicamente orientada que se referem aos elementos clássicos de um sistema de hidrocarbonetos, isto é, reservatório, selagem, trape, fonte, maturação, e migração. Os atributos são espacialmente correlacionados e comparados com um catálogo de configurações de trape de hidrocarbonetos para determinar a presença potencial de trapes de hidrocarbonetos e para avaliar a confiança de sua existência.

[0036] Em uma forma de realização, a invenção é um método implementado por computador para analisar um volume composto de voxeis de dados sísmicos representando uma região de subsuperfície quanto à presença de um sistema de hidrocarboneto ou um play particular, compreendendo: - dividir o volume de dados sísmicos para formar uma pluralidade de segmentos; e - classificar a pluralidade de segmentos quanto à presença de um sistema de hidrocarbonetos ou play particular, com base pelo menos parcialmente nos scores de prospectividade para os voxeis de dados sísmicos em cada segmento; em que o score de prospectividade é baseado na computação de pelo menos dois atributos que se referem a diferentes elementos de um sistema de hidrocarbonetos ou play particular.

[0037] Como com qualquer método de processamento de dados geofísico, a invenção em aplicações práticas é altamente automatizada, isto é, é realizada com a ajuda de um computador programado de acordo com as descrições aqui.

BREVE DESCRIÇÃO DOS DESENHOS

[0038] A presente invenção e suas vantagens serão melhor compreendidas por referência à seguinte descrição detalhada e aos desenhos anexos, em que:

[0039] A Fig. 1 ilustra elementos do Sistema de Hidrocarbonetos para um trape anticlinal;

[0040] A Fig. 2 é um fluxograma mostrando etapas básicas em uma forma de realização da invenção;

[0041] A Fig. 3 é um diagrama esquemático ilustrando uma forma de realização do presente método inventivo;

[0042] A Fig. 4 ilustra como o estêncil de filtro de um filtro de convolução integral de linha segue as direções de tangente;

[0043] A Fig. 5 ilustra como o filtro de convolução integral de linha pode ser feito mais robusto;

[0044] A Fig. 6 ilustra que o estêncil de filtro de convolução integral de linha pode ter uma espessura estendida.

[0045] A Fig. 7 mostra um exemplo de conversão de um atributo local em um regional (à esquerda: dados sísmicos, metade-esquerda: atributo de convergência local, metade-direita: campo vetorial tangente, e à direita: convergência regional);

[0046] A Fig. 8 mostra um exemplo dos atributos de convergência tridimensional com um componente de magnitude de convergência (claro: convergência forte, escuro: nenhuma convergência) e um componente de orientação de convergência;

[0047] A Fig. 9 define relações de terminação estratigráfica;

[0048] A fig. 10 representa um trape estrutural delimitado por dobra (anticlinal);

[0049] A pig. 11 representa um trape estrutural relacionado com falha;

[0050] A fig. 12 mostra um trape estrutural relacionado com sal;

[0051] A pig- 15 ilustra um trape estratigráfico relacionado com um pinchout do reservatório;

[0052] A Fig- 14 mostra um trape estratigráfico relacionado com uma incorform idade;

[0053] A Fig- 15 representa um trape estratigráfico formado por relevo erosivo soterrado;

[0054] A Fig- 16 representa um trape estratigráfico formado por diferenças diagenéticas;

[0055] A Fig- 17 apresenta um modelo de sequência deposicional esquemática, isto é, diagrama de lingote;

[0056] A Fig- 18 ilustra como os diferentes elementos podem ser integrados para avaliar a chance de sucesso para uma acumulação de hidrocarbonetos;

[0057] A Fig. 19 mostra uma aplicação esquemática de classificação de quatro alvos potenciais; e

[0058] A Fig. 20 ilustra uma aplicação exemplar do presente método inventivo.

[0059] A invenção será descrita em combinação com as formas de realização exemplares. Entretanto, na medida em que a seguinte descrição detalhada é específica a uma forma de realização particular ou a um uso particular da invenção, esta destina-se a ser somente ilustrativa e não é para ser interpretada como limitante do escopo da invenção. Ao contrário, destina-se a abranger todas as alternativas, modificações e equivalentes que possam ser incluídos dentro do escopo da invenção, como definido pelas reivindicações anexas.

DESCRIÇÃO DETALHADA DAS FORMAS DE REALIZAÇÃO EXEMPLARES

[0060] A invenção, em um de seus aspectos, compreende computação de atributos sísmicos que se relacionam com elementos do sistema de hidrocarbonetos, análise dos dados para os elementos do sistema de hidrocarbonetos e/ou reconhecimento de plays específicos ou estilos conceituais, e classificação e anotação destas regiões com base no tamanho, qualidade, e confiança dos prospectos. Algumas definições são dadas a seguir.

[0061] Embora o termo possa ser usado mais ampla ou estritamente em outro lugar, um sistema de petróleo ou hidrocarbonetos é geralmente empregado aqui para significar um sistema natural que abrange uma cápsula de rocha matriz ativa e todo óleo e gás relacionados. Isso inclui todos os elementos e processos geológicos que são essenciais para um acúmulo de hidrocarboneto existir, como ilustrado na Fig. 1. Os hidrocarbonetos encontrados na realidade incluem altas concentrações de gás térmico e/ou biogênico encontrados em reservatórios convencionais ou em hidratas gasosos, reservatórios apertados, xisto fraturado, ou hulha; e condensados, óleos brutos, óleos pesados, asfaltas e alcatrões. O termo “sistema” descreve os elementos e processos independentes que formam a unidade funcional que cria acumulações de hidrocarbonetos. Os elementos essenciais incluem uma rocha matriz de petróleo (source), rocha reservatório (reservoir), rocha de selagem (seal), e rocha de cobertura (overburden). Os processos são a formação do trape e a maturação (geração), migração, e acumulação de hidrocarbonetos. Está implícita nesses processos uma sequência ou distribuição de eventos.

[0062] Uma definição alternativa do sistema de hidrocarbonetos pode incluir somente a rocha matriz, os processos de maturação e migração, e sua distribuição; neste caso, reservatório, selagem e trape podem ser definidos para formar um play. Para fins de explicação do presente método inventivo, o termo sistema de hidrocarboneto é definido abranger fonte, reservatório, selagem, trape, maturação, migração e distribuição. Além disso, o termo play é geralmente usado aqui para indicar uma combinação específica e arranjo de reservatório, selagem e geometria de aprisionamento.

[0063] A rocha matriz é uma rocha rica em matéria orgânica que, se aquecida suficientemente, gerará óleo e/ou gás durante o tempo. As rochas matrizes comuns incluem xistos ou calcários. As rochas de origem marinha tendem a ser propensas a óleo, enquanto as rochas matrizes terrestres (tais como hulha) tendem a ser propensas a gás. A preservação de matéria orgânica sem degradação é crítica para criar uma rocha de boa fonte, e necessária para um completo sistema de petróleo.

[0064] O reservatório é um corpo de rocha de subsuperficie tendo suficiente porosidade e permeabilidade para receber, armazenar, e transmitir fluidos. As rochas sedimentares são a maioria das rochas reservatório comuns, em razão de terem mais porosidade do que a maior parte das rochas ígneas ou metamórficas e formarem condições de temperatura inferior em que hidrocarbonetos podem ser preservados. Um reservatório é um componente crítico de um sistema de petróleo completo.

[0065] Rocha de selagem é uma rocha relativamente impermeável, geralmente xisto, anidrito, ou sal, que forma uma barreira ou tampa acima e parcialmente em tomo da rocha reservatório, de modo que fluidos não possam migrar além do reservatório. Uma rocha de selagem é um componente crítico de um sistema de petróleo completo.

[0066] Rocha de cobertura é a rocha sobre o topo da fonte e reservatório. No contexto do sistema de petróleo, sua função principal é formar um cobertor grosso por cima da fonte, onde aumenta a temperatura e pressão até o grau necessário para converter matéria orgânica em hidrocarbonetos.

[0067] Trape é uma configuração de rochas adequadas para conter hidrocarbonetos e é vedada por uma formação relativamente impermeável através da qual hidrocarbonetos não migrarão. Os trapes são descritos como trapes estruturais (nos estratos deformados, tais como dobras e falhas) ou trapes estratigráficos (em áreas em que tipos de rochas mudam, tais como desconformidades, pinch-outs e recifes) ou suas combinações. Para trapes estruturais, a deformação deve ocorrer antes da migração de hidrocarbonetos, ou os hidrocarbonetos não se acumularão. Um trape é um componente essencial de um sistema de petróleo.

[0068] Geração ou maturação é a formação de hidrocarbonetos de uma rocha matriz, visto que betume se forma de querogênio e acumula-se como óleo ou gás. A geração depende de três fatores principais: presença de matéria orgânica suficiente para produzir hidrocarbonetos, temperatura adequada, e suficiente tempo para trazer a rocha matriz à maturidade. A pressão e a presença de bactérias e catalisadores também afetam a geração. As pressão e temperatura insuficientes, causadas, por exemplo, por um enterramento superficial com uma cobertura fina, renderão uma fonte imatura e a geração será escassa ou incompleta. As pressão e temperatura excessivas, causadas, por exemplo, por enterramento profundo sob uma cobertura grossa, causarão degradação do óleo gerado em gás e, subsequentemente, em dióxido de carbono e água. A geração é uma fase crítica no desenvolvimento de um sistema de petróleo.

[0069] Migração é o movimento de hidrocarbonetos de sua fonte para dentro das rochas reservatório. O movimento de hidrocarbonetos recentemente gerados para fora de sua rocha matriz é a migração primária, também chamada de expulsão. O outro movimento de hidrocarbonetos para dentro da rocha reservatório em um trape de hidrocarbonetos ou outra área de acumulação é a migração secundária. A migração tipicamente ocorre de uma área estruturalmente baixa para uma área mais elevada, em razão da flutuação relativa de hidrocarbonetos em comparação com a rocha circundante. A migração pode ser local ou pode ocorrer ao longo de distâncias de centenas de quilômetros em grandes bacias sedimentares e é crítica para a formação de um sistema de petróleo viável.

[0070] Acumulação refere-se tanto a uma ocorrência de hidrocarbonetos aprisionados, isto é, um play ou um campo de óleo ou gás, como à fase de desenvolvimento de um sistema de petróleo, durante a qual hidrocarbonetos migram para dentro e permanecem aprisionados nas rochas reservatório.

[0071] A distribuição refere-se à ordem relativa em que elementos são formados ou modificados, ou a ordem em que os processos ocorrem. Um trape pode acumular hidrocarbonetos migrando somente se for formado antes da migração. Um trape pode não ser preenchido se a migração ainda não tiver alcançado sua localidade. Um trape pode perder sua carga, pelo menos parcialmente, se a rocha de selagem for rompida após acumulação.

[0072] Um play é um modelo conceptual para um estilo de acumulação de hidrocarbonetos, com frequência usado para desenvolver prospectos em uma bacia, região ou tendência, ou usado para continuar a explorar uma tendência identificada. Um play (ou um grupo de plays interrelacionados) geralmente ocorre em um único sistema de hidrocarbonetos e pode ser compreendido por um grupo de prospectos similares.

[0073] Um prospecto é uma área em que hidrocarbonetos foram preditos existirem em quantidade econômica. Um prospecto é com frequência uma anomalia, tal como uma estrutura geológica ou uma anomalia de amplitude sísmica que é recomendada como um local para perfuração de um poço para verificação de quantidades econômicas de hidrocarbonetos. A justificativa para a perfuração de um prospecto é feita unindo-se evidências para um sistema de hidrocarboneto ativo, ou demonstrando-se probabilidades razoáveis de encontrar rocha reservatório de boa qualidade, um trape de suficiente tamanho, rocha de selagem adequada, e condições apropriadas para geração e migração de hidrocarbonetos para preencher o reservatório. Para fins do método inventivo, o prospecto é usado amplamente para indicar uma área que é recomendada para outra análise detalhada.

[0074] Como citado acima, o que é necessário é um método que automaticamente analise dados sísmicos quanto à presença de elementos do sistema de hidrocarbonetos, regiões de sinalizadores onde elementos de play são justapostos em configurações favoráveis ou consistentes com um play conhecido ou especificado, e classifique estes prospectos com referência a seu potencial de acumulação de hidrocarbonetos. Um tal sistema focaliza a análise e interpretação em áreas mais prospectivas. Além disso, o sistema pode reconhecer o tipo de play e prover um score de confiança para elementos individuais. Para cada prospecto, elementos não identificados ou elementos com baixa confiança são vitais e demandam atenção especial durante subsequentes análises para não expor ao perigo o prospecto. O sistema pode ser empregado várias vezes durante o ciclo de vida de uma região ou ativo. Primeiro, o sistema pode ser usado em dados regionais, tipicamente bidimensionais, para identificar áreas de prospecção, por exemplo, para preparar uma oferta de um bloco ou para localizar uma pesquisa de exploração tridimensional. Na fase de exploração, o sistema pode ser usado para examinar cuidadosamente, para focalizar e guiar o intérprete para subconjuntos do volume de dados sísmicos menores e mais controláveis. No estágio de produção, o sistema pode ser usado para localizar menores prospectos próximos de infreestrutura anteriormente localizada. Finalmente, o sistema pode ser usado pare garantir que nenhum prospecto foi inspecionado, antes de abandonar ou negociar o ativo.

[0075] A Fig. 2 é um fluxograma mostrando etapas básicas, em uma forma de realização do método inventivo, que incluem duas etapas determinadas (caixas tendo linhas contínuas) e seis opcionais (caixas de linhas tracejadas) que dependem de entradas, suposições feitas, e o modo em que o sistema é empregado. Na etapa 2, atributos são gerados que se referem a elementos do sistema de hidrocarbonetos ou elementos de play específicos. Na etapa 5, os dados são analisados e locais são classificados com respeito aos atributos, assim identificando prospectos potenciais que são então armazenados pare outra análise ou visualização.

[0076] A Etapa 1 opcional define um tecido da subsuperficie, isto é, as geometrias dos estratos de subsuperficie e sua deformação por falhas e dobres, que formam a base para a computação de alguns atributos ou podem ser necessárias pare outras computações, pare agregar informações de estratos similares. Muitos atributos, que são úteis pare o sistema inventivo, são formados integrando-se ou calculando-se a média dos atributos mais tradicionais ao longo do tecido. Na Etapa opcional 3, os dados são separados em pelo menos dois segmentos ou divisões pare análise e definição de prospecto. A etapa opcional 4 é a formação de scores normalizados para elementos do sistema de hidrocarbonetos empregando um ou mais dos atributos. A etapa opcional 6 é a análise dos prospectos formados e classificados; por exemplo, identificação do menor elemento(s) expresso que deve ser ainda examinado. A etapa 7 é a definição ou seleção de pelo menos um tipo de conceito específico ou play que é pare ser procurado em vez de uma busca genérica de elementos vizinhos do sistema de hidrocarbonetos. Finalmente, a etapa 8 é a definição, criação e controle de um catálago de tipos de play e a configuração de seus elementos.

[0077] Uma aplicação esquemática do método inventivo é apresentada na Fig. 3, em que atributos sísmicos 32 são computados a partir de dados sísmicos 31. Os atributos são então combinados 33 em regiões realçadas (por exemplo, 36) que podem constituir plays estratigráficos, neste caso o play de monte-encoberto 34 e o play pinch-out estratigráfico 35.

[0078] Definição de Tecido

[0079] Um atributo é uma propriedade mensurável de dados sísmicos, tal como amplitude, inclinação, frequência, fase e polaridade. Os atributos podem ser medidos em um instante do tempo ou através de uma janela de tempo, e podem ser medidos em um traço único ou conjunto de traços, ou em uma superfície interpretada de dados sísmicos, uma janela de dados, ou mesmo em múltiplos volumes sísmicos simultaneamente. O presente método inventivo emprega muitos atributos clássicos, que são bem conhecidos àqueles experientes na arte de analisar atributo ou interpretação sísmica. Para o método inventivo, entretanto, alguns destes atributos são modificados, por exemplo, integrando-se ao longo estratos ou estrutura, para enfatizar variações regionais através de certos locais. Além disso, novos atributos geológicos (2) são descritos que diretamente se referem a elementos do sistema de hidrocarbonetos ou elementos de play, e a configurações particulares destes elementos. Uma vez que atributos geológicos podem ser abrangidos em modos alternativos, e diferentes atributos geológicos podem referir-se ao mesmo elemento ou configuração de elementos, scores (4) podem ser computados que sintetizem as diferentes formas de realização e atributos entre si para uma medida que avalie a chance de encontrar um elemento de play particular ou do sistema de hidrocarbonetos, ou um arranjo geométrico particular de tais elementos. Os scores são então combinados (6) em probabilidades para as configurações que podem ser usadas para visualização, análise, ou taxação e classificação.

[0080] Com respeito à Fig. 2, menos do que todas as etapas ilustradas podem ser requeridas para implementar uma forma de realização particular da invenção. As etapas individuais podem ser combinadas ou separadas em múltiplos componentes. Além disso, metodologias adicionais e/ou alternativas podem empregar etapas adicionais não mostradas aqui. Embora o fluxograma ilustre várias ações ocorrendo consecutivamente, deve-se também observar que algumas ações poderiam ocorrer em série, substancialmente em paralelo, e/ou em pontos de tempo consideravelmente diferentes. As etapas podem também ser repetidas. Um exemplo destas variações é a etapa 1, a definição de tecido, estrutura ou estrutura de camada.

[0081] O tecido pode ser necessário para alguns atributos empregados na invenção e, assim, quando se empregando somente tais atributos, a definição de tecido pode ser ignorada. A maioria dos atributos usados para o método inventivo, entretanto, requererá um tecido para sua formação. O tecido poderia ser computado em movimento para todo o atributo quando necessário. Preferivelmente, entretanto, o tecido é gerado somente uma vez, em seguida armazenado e usado repetidamente para gerar, modificar, ou integrar atributos. Na Etapa 1, a formação de tecido é mostrada como uma etapa opcional, porque, dependendo dos atributos empregados nas diferentes formas de realização do método inventivo, o tecido deveria ser computado uma vez, repetidamente, ou de modo algum.

[0082] Uma vantagem da geração de tecido uma vez é a consistência entre diferentes atributos, porque todos são baseados no mesmo tecido. Uma desvantagem de gerá-lo somente uma vez é que o tecido resultante precisa ser representado e armazenado de alguma maneira que possa ser ineficiente ou mesmo inadequado para algumas aplicações de geração, modificação, ou integração de atributo. Altemativamente, cada aplicação pode gerar seu próprio tecido em uma representação adequada e eficaz que pode causar deficiências para o sistema total, porque o mesmo tecido com as mesmas representações pode ser computado múltiplas vezes. Já outra opção é computar o tecido com alguns métodos e representações escolhidos, que são armazenados e assim reutilizáveis. Com este método alternativo, as desvantagens são exigências de armazenagem aumentadas e repetidos acessos aos meios de armazenagem. Na prática, o método preferido pode ser computar e armazenar o tecido ou seus componentes para alguns dos métodos e representações mais comuns, e computar menos tecidos quando necessário.

[0083] Os métodos para computar o tecido podem ser baseados em estimativas de inclinação, orientação ou colisão; tensores de estrutura, correlações de forma de onda, ou esqueletizações. Um modo particular de estimar inclinação e orientação de inclinação é pelo uso de gradientes em dados sísmicos nas duas direções horizontais e na direção vertical (por exemplo, a Patente U.S. No. 6.850.864 B2, “Method for analyzing dip in seismic data volumes” de Gillard et al). Outro estimador de inclinação é baseado em filtros de desconstrução de onda plana (por exemplo, “Applications of plane-wave destruction filters” de Fomel, Geophysics 67(6), 1946- 1960, (2002)). Os métodos baseados no tensor de estrutura permitem a computação de vetores normais e tangentes em reflexões sísmicas (por exemplo, “Estimating Seismic Heterogeneity with the Structure Tensor” de Imhof, 67th EAGE Conference & Exhibition, (2005); ou a Patente U.S. No. 7.454.,292 B2, “Inverse-vector method for smoothing dips and azimuths”, de Wang et al). As inclinações podem ser resolvidas em deslocamentos de tempo para as colhetas automáticas de superfície (“Flattening without picking”, de Lomasket al, Geophysics 71(4), P13-P20 (2006)).

[0084] Um método alternativo de computar inclinação e azimute de reflexão é baseado no gradiente de isocontornos. O gradiente é um vetor local que caracteriza a direção de incremento mais íngreme de uma função. Considerando dados sísmicos, amplitude, por exemplo, como uma função f em espaço tridimensional (x,y,t), o gradiente é definido como Vf. Uma vez que há três componentes no vetor de gradiente resultante Vf, somente três amostras vizinhas são necessárias para avaliar Vf para primeira ordem. Empregando-se amostras mais circundantes, entretanto, permite-se a estimativa do gradiente de um sistema superdeterminado, por exemplo, no sentido dos quadrados mínimos, ou empregando-se uma aproximação de ordem mais elevada para a aproximação de diferença-finita. Em qualquer dos dois casos, o gradiente resultante será mais uniforme. A seleção de pontos vizinhos não é limitada a qualquer sequência ou padrão particular. De preferência, a vizinhança pode ser definida como pontos ou voxeis dispostos em uma esfera, um cilindro, uma caixa, ou qualquer outro tipo de formato circundando o ponto de interesse. Além disso, este formato pode ser alinhado com o tecido ou com uma sua estimativa em todo o ponto. Os atributos derivativos podem então ser derivados deste gradiente, incluindo magnitude e direção, isto é, inclinação e azimute em 3D. O gradiente, como ordinariamente definido, sempre aponta na direção de amplitude aumentando. Em dados sísmicos com ciclagens de amplitude de vales negativos para picos positivos e vice-versa, o resultado seria que as camadas paralelas exibiriam valores de inclinação opostos. Uma solução é normalizar o gradiente com respeito a um meio-espaço virando em sua direção, por exemplo, para inverter o sinal do gradiente, de modo que todos os pontos gradientes sejam orientados na direção positiva y:

[0085] Este cálculo de gradiente é uma formulação mais robusta do que os métodos de cálculo de gradiente existentes. Também tem a capacidade de multiescalas, uma vez que o tamanho de janela (número de pontos vizinhos), dentro do qual o gradiente é computado, pode ser variado.

[0086] A correlação lateral de formas de onda é outro método de gerar um tecido. Em uma primeira etapa, eventos, tais como picos, vales, e/ou cruzamentos zero, são determinados. Em uma segunda etapa, formas de onda em tomo destes eventos são comparadas com os eventos de traços vizinhos para encontrar e conectar os mais similares. Observa-se a assimetria que ocorre pelo fato de cada evento conectar-se somente com um outro evento, isto é, o mais similar a ele. Múltiplos eventos, entretanto, podem conectar-se a um e ao mesmo evento. O resultado é uma estrutura de gráfico direcionada em que os vértices correspondem a eventos, tais como picos, vales, e/ou cruzamentos zero; e bordas conectam eventos individuais ao mais similar dos traços vizinhos. Partir de um evento em um traço particular é com frequência possível para seguir bordas direcionadas entre eventos e traços e depois retornar para um diferente evento no traço original, que indica inconsistências estratigráficas (ou topológicas). Para distinguir este método de outros, os termos esqueleto bruto ou esqueletização serão usados, se necessários.

[0087] Um método preferido para computar um tecido é por esqueletização topológica, que automaticamente cria superfícies baseadas em reflexão de uma maneira topologicamente consistente, em que superfícies individuais não se sobrepõem e conjuntos de múltiplas superfícies são consistentes com princípios de sobreposição estratigráficos (Publicação do Pedido de Patente PCT WO 2009/142872 Al, “Seismic Horizon Skeletonization”, de Imhof et al.). Para distinguir este método de outros, os termos esqueleto topológico ou esqueletização topológica são usados, se necessários. Neste método, as superfícies topologicamente consistentes são rotuladas monotonicamente em um modelo de topo-descendente, que permite a designação de uma pseudoidade para todas as amostras de dados sísmicos e a conversão dos dados sísmicos da profundidade geofísica tradicional ou domínio de tempo-de-duas-direções para um domínio de idade geológica (por exemplo, Patente U.S. No. 6.850.845 B2, “System for multidimensional data analysis”, de Stark; Publicação do Pedido de Patente PCT WO 2009/142872 Al, “Seismic Horizon Skeletonization”, de Imhof et al.; e continua na Publicação do Pedido de Patente U.S. No. 2010/0149917 “Method for Geophysical and Geological Interpretation of Seismic Volumes in Depth, Time, and Age”, de Imhof et al.). Um volume de mapeamento de idade resultante, e/ou volume de mapeamento de profundidade, permite a transferência de dados entre o domínio de profundidade (ou tempo-de-duas-direções) geofísica e o domínio de idade geológica construído e vice- versa. Cada fatia horizontal através de urn volume de mapeamento de profundidade, corresponde a um mapa de estrutura de profundidade (ou tempo) para um horizonte particular.

[0088] Já outro método de esqueletização é baseado em adelgaçamento morfológico; o resultado será designado esqueleto morfológico. Dados sísmicos são bloqueados para imagens binárias, por exemplo, reduzindo-se os dados para apenas suas polaridades, por exemplo, ±1. As faixas de valor +1 são reduzidas a linhas com um valor de +1 com igual conectividade, por exemplo, por aplicação de adelgaçamento morfológico. As faixas de valor -1 são reduzidas a linhas com um valor de -1. Todas as outras amostras são determinadas em zero. Um resultado similar é obtido do atributo de polaridade evidente que é formado pela polaridade de fase instantânea computada na extrema amplitude local.

[0089] Outro aspecto de computação do tecido refere-se a falhas que podem ser detectadas como descontinuidades em dados sísmicos (por exemplo, Patente U.S. No. 6.516.274 B2, “Method for imaging discontinuities in seismic data using dip- steering”, de Cheng et al.; ou Publicação do Pedido de Patente PCT WO 2009/082545 Al (“Detection Of Features In Seismic Images”, de Kumaran e Wang). Segmentos de falha detectados podem ser ainda completamente limpos ou refinados (por exemplo, (Patente U.S. No. 7.203.342 B2, “Image feature extraction”, de Pedersen).

[0090] Atributos Geológicos

[0091] 1.a Geração de Atributos Regionais de Estrutura de Longa Distância Suavemente Orientados

[0092] Um atributo geológico é um atributo sísmico que realça ou quantifica algum aspecto do sistema de hidrocarbonetos ou um play. Na, realidade, os atributos geológicos são com frequência regionais. Eles podem, por exemplo, estimar uma propriedade de dados sísmico através de dezenas ou centenas de voxeis. Através de tais distâncias, os atributos são preferivelmente computados ao longo do tecido em uma estrutura-orientada de maneira a evitar mistura e, assim, contaminação de outros cortes de estratos através da janela de análise.

[0093] Dado o tecido como definido pelo volume de dados sísmicos original, qualquer atributo pode ser convertido em um atributo regional, integrando-se ou calculando-se a média dele ao longo do tecido. Um tal método é baseado em filtragem de difusão anisotrópica não linear (por exemplo, Patente U.S. No. 6.725.174 B2, “Edge-preserving enhancement of seismic images by nonlinear anisotropic diffusion”, de Bouts et al.), que aplaina os dados sísmicos predominantemente ao longo da direção das reflexões. Tradicionalmente, o tecido ou as direções das reflexões são computados empregando-se tensores de estrutura sísmica para o dado conjunto de dados sísmicos que é para ser aplainado, porém o tecido pode ser definido com qualquer método. Alguns são descritos aqui.

[0094] Na presente invenção, a difusão anisotrópica não linear pode ser usada para computar o tecido para o volume de amplitude sísmica primário, porém, em seguida, o tecido é usado para aplainar um conjunto de dados secundários ao longo da direção das reflexões do volume de amplitude sísmica primária. Em alguns casos, este conjunto de dados secundários não pode mesmo ser estendido em camadas ou faixas, o que evita a computação de seu próprio tensor de estrutura. Em casos em que o conjunto de dados secundário é estendido em faixas, ele pode, contudo, ser vantajoso para aplainar ao longo do tecido de um volume sísmico de referência, por exemplo, para aumentar a consistência entre conjuntos de dados diferentes após aplainamento de estrutura-orientada.

[0095] Um método alternativo de estrutura orientada, aplainando através de longas distâncias da presente invenção, é baseado em convolução integral de linha (LIC). A convolução integral de linha é um algoritmo de síntese de textura bem conhecido usado em processamento de imagem ou visualização de dados (“Imaging Vector Fields Using Line Integral Convolution” by Cabral and Leedom, Proceedings of ACM SigGraph 93, 263-270, (1993)). Um filtro de baixa-passagem é usado para envolver uma textura de interferência de entrada, por exemplo, uma imagem de interferência aleatória, ao longo de formas aerodinâmicas simetricamente bidirecionais de pixel-centralizado para explorar a correlação espacial na direção de fluxo. Com frequência, as formas aerodinâmicas não são representadas por curvas, porém por um campo vetorial cujos vetores são tangentes às formas aerodinâmicas, e as formas aerodinâmicas são computadas por integração.

[0096] A substituição de alguns dados de atributo sísmico para a textura de interferência, e usando o tecido para definir as formas aerodinâmicas, permite usar convolução integral de linha para filtrar o atributo dado ao longo do tecido, que é uma nova abordagem para filtragem de estrutura-orientada. As definições de tecido preferidas são baseadas na inclinação, por exemplo, quando computadas de tensores de estrutura. Mais eficiência é ganha pela discretização do tecido ou, mais especificamente, pelas inclinações ou tangentes de reflexão, a 0, 45, 90, 135, 180, 225, 270, ou 315°; isto é, puramente lateral, puramente vertical, ou puramente diagonal, o que substitui a integração numérica e interpolação associada com movimento ao longo de uma lista associada. Uma vez que os dados sísmicos são tipicamente mais ou menos horizontalmente estendidos em faixas, mais eficiência é ganha forçando-se todas as tangentes para irem da esquerda para a direita, o que reduz as inclinações para 0, +45, e -45°. A Figura 4 ilustra a aplicação de um filtro de convolução de linha em que as cruzes cinzas indicam o local das amostras e as setas indicam as direções de inclinação discretizadas. Os filtros são aplicados aos locais indicados pelos pontos pretos encorpados. Os desenhos mostram como os estênceis de filtro centralizados nos pontos seguem os vetores e, assim, ilustram o tecido colhendo os valores da amostra, marcados com cruzes pretas, ao longo destes trajetos.

[0097] Para tornar o filtro de convolução integral de linha mais robusto na presença de falhas e outras complexidades, pode ser feita a preservação de borda aplicando-se o filtro em múltiplas configurações de estêncil, que não são necessariamente centralizadas no local de amostra atual (ilustrado na Fig. 5). Fora de todas as configurações de estêncil testadas, a que exibe a menor variação entre suas amostras é preferível e finalmente aplicada em um local particular. Dependendo dos coeficientes de filtro e da natureza do atributo que está sendo tratado com um tal filtro de convolução integral de linha de preservação de borda, critérios que não a variação mínima podem ser aplicados. Outros possíveis critérios incluem variação máxima, saída mínima ou máxima, e/ou evitação de locais de amostra especialmente sinalizados.

[0098] Os coeficientes de filtro não precisam ser uniformes, o que permite aplicação de diferenciadores (filtro de alta-passagem), integradores (filtro de baixa- passagem), ou suas combinações. Os filtros podem ser lineares ou não lineares, por exemplo, filtros medianos ou filtros de modo. A Figura 6 ilustra que o estêncil de filtro pode ter um formato arbitrário que é deformado para conformar-se ao trajeto de vetor através do local de amostra atual. A Figura 7 mostra uma aplicação exemplar de tais ferramentas de filtragem, convertendo um atributo local em um atributo regional.

[0099] Já outro método usado para realizar filtração de estrutura orientada através de longas distâncias envolve a exploração de superfícies obtida por esqueletização sísmica (Publicação da Patente PCT No. WO 2009/142872 Al, “Seismic Horizon Skeletonization”, de Imhof et al.), em que as superficies de esqueleto proveem o tecido ao longo do qual os atributos são filtrados. Em vez de aplainamento ao longo de diferentes superfícies, que podem requerer interpolação, o aplainamento pode ser realizado ao longo de isosuperfícies do volume de mapeamento de idade construído das superfícies de esqueleto. Alternativamente, o atributo pode ser transformado pelo domínio tradicional de profundidade ou tempo de deslocamento de duas direções geofísico para o domínio de idade geológica. No domínio de idade geológica, uma fatia horizontal corresponde a um horizonte e, assim, aplainamento ou filtração de fatias horizontais é, por definição, estrutura- orientada (Pedido de Patente U.S. No. 12/623.034 “Method For Geophysical and Geological Interpretation of Seismic Volumes In Depth, Time, and Age”, de Imhof et al.). Após filtragem no domínio de idade, os resultados podem ser transformados de volta para o domínio de profundidade ou tempo de deslocamento de duas direções.

[00100] Independente do método utilizado para filtragem de estrutura-orientada, atributos geológicos podem ser computados como atributos de multiescalas em razão do tamanho de filtro ser especificado pelo intérprete ou sistema. Em vez de computar um atributo geológico em somente uma escala, um intérprete ou sistema pode escolher múltiplos tamanhos de filtro, permitindo computação do mesmo atributo geológico em diferentes escalas.

[00101] 1.b Atributos Hessianos

[00102] Um novo atributo sísmico da presente invenção é o Hessiano ou tensor derivado (dyadic) espacial de segunda ordem H = VVf, que é computado em um modelo similar que o do gradiente modificado. Em todo local de interesse, o Hessiano pode ser representado nas três dimensões por uma matriz 3x3 simétrica, isto é, por seis componentes independentes. Para calcular o Hessiano em qualquer local de voxel, pelo menos seis voxeis vizinhos são necessários, uma vez que ele tem seis componentes independentes. O cálculo, entretanto, pode envolver mais voxeis para obter aplainamento e mais estimativas estáveis para o Hessiaπo: calculando-se seus componentes de um sistema superdeterminado, por exemplo, no sentido dos quadrados mínimos, ou empregando-se uma aproximação de ordem mais elevada para a aproximação de diferença-finita das derivadas espaciais. A seleção de pontos vizinhos não é limitada a qualquer sequência ou padrão particular. De preferência, a vizinhança pode ser definida como pontos ou voxeis dispostos em uma esfera, um cilindro, uma caixa, ou qualquer outro tipo de formato circundando o ponto de interesse. Além disso, este formato pode ser alinhado com o tecido ou com uma sua estimativa em todo o ponto.

[00103] De um ponto de vista geométrico, o volume de dados sísmicos consiste de aspectos pontuais (por exemplo, ruído), aspectos curvilíneos (por exemplo, terminações ou canais), aspectos de superfície (por exemplo, falhas), ou aspectos de volume (por exemplo, correias de canal ou corpos salinos) embutidos em estrutura de segundo plano encarnada. Em práticas recentes, derivadas de amplitude de segunda ordem direcionais ou sinais de correlação de traço são geralmente calculados para detectar heterogeneidades que tipicamente ocorrem em aspectos geológicos, tais como falhas, limites de erosão, e pinch-outs. Faltando derivadas-cruzadas, atributos baseados em derivadas de segunda ordem tradicionais representam variações locais dos dados incompletamente. Com seus seis componentes, o Hessiano caracteriza a convexidade ou concavidade local dos dados sísmicos. As derivadas de segunda-ordem tradicionais são encontradas na diagonal do Hessiano, enquanto as derivadas-cruzadas são encontradas fora da diagonal. Os atributos podem assim ser formados destes seis componentes e suas combinações. Os componentes e, portanto, os atributos, são formados de derivadas de segunda ordem de um volume sísmico e mostram mudança espacial em gradiente. Eles refletem mudanças espaciais em ângulos de inclinação e azimute do volume sísmico, o que permite realçar regiões anômalas e grupos de aspectos interrelacionados. Realçando-se diferentes tipos de “descontinuidades” e suas interrelações espaciais, volumes sísmicos podem ser melhor caracterizados e visualizados.

[00104] Outro uso do Hessiano é para servir como uma etapa intermediária para calcular curvaturas espaciais, como ilustrado na sexta seção. Além disso, eigenvalores da matriz Hessiana e projeções do Hessiano em certas superfícies também formam medidas de quanto curvados os isocontornos são.

[00105] A Curvatura mede o grau em que um objeto, uma superfície, por exemplo, desvia de ser plana. Entretanto, definições específicas dependem do contexto. Por exemplo, a curvatura pode ser definida para um plano ou uma curva de espaço, em cujo caso ela é uma escalar. Pode também ser definida para uma superfície parametrizada embutida no espaço tridimensional, em cujo caso ela é um tensor. Com o objetivo de descobrir atributos relacionados à estratigrafia, dois novos tipos de curvatura são definidos aqui. O primeiro é curvatura de isocontorno em uma imagem sísmica bidimensional. Uma seção cruzada bidimensional arbitrária através de um volume sísmico é tratada como um mapa de uma função f(xy). Um isocontorno é a curva desta seção cruzada cujas coclassificadas satisfazem f(x,y)= constante. Sua curvatura é definida como

em que todas as quantidades são definidas por componentes do vetor gradiente de isocontorno ou tensor Hessiano. Este atributo de curvatura pode ser particularmente útil em realçar aspectos geométricos de alta curvatura, tais como terminações, cantos, e junções.

[00106] Outra definição de curvatura de isocontorno é baseada em uma generalização tridimensional. Os dados sísmicos são tratados como uma função de três variáveis independentes f(x,y). Todos os voxeis que satisfazem f(x,y) = constante constituem uma isosuperfície, isto é, uma superfície curvada embutida no espaço tridimensional. Para um espaço curvado, a curvatura é uma quantidade tensora que permite extração de múltiplas quantidades escalares relacionadas à curvatura, tal como a curvatura principal, curvatura média, e curvatura Gaussiana. Para uma tal superfície curvada, pode-se definir uma superfície normal por normalização do gradiente de isocontorno n = Vf /|Vf|. Tangentes de superfície podem ser definidas usando-se a condição t n = 0. Observa-se que há um número finito de tangentes que satisfazem esta condição em qualquer ponto dado sobre a superfície. E há igualmente muitos planos ortogonais à superfície definidos por uma tangente t e a normal n. A interseção de tal plano com a superfície é uma curva, cuja curvatura é chamada de a curvatura normal. O máximo e mínimo de todas as curvaturas normais são designados como as curvaturas principais K1 e K2, em que a curvatura Gaussiana é computada como o produto K=K1 K2, enquanto sua média algébrica obtém a curvatura média M=(K1 + K2)/2.

[00107] Todas estas quantidades são definidas para cada ponto do espaço, e podem ser expressas compactamente como uma matriz S conhecida como o Operador de Formato. Para construção de uma forma de realização conveniente de S, observa-se que o espaço de todas as tangentes do ponto dado pode ser expresso com os eigenvetores de uma matriz T obtida subtraindo-se o produto externo do normal com a própria matriz de identidade, T = (I - nm). T eficazmente projeta qualquer vetor tridimensional sobre o plano de tangente. Assim, pode-se medir a mudança do gradiente restrito para aquele plano multiplicando-se a matriz Hessiana por T ou, simbolicamente, H-H T, em que ■ indica o interior, ou matriz, produto. Finalmente, S é obtido normalizando-se o Hessiano: S=H’/|Vf|2. As curvaturas principais são os eigenvalores de S, dividindo o traço de S por dois campos, a curvatura média, M=traço(S)/2, e a determinante det(S)=K definem a curvatura Gaussiana.

[00108] Todas estas medidas de curvatura podem ser tratadas como atributos sísmicos e podem ser usadas para detectar terminações de horizonte e junções que definem aspectos de importância geológica, tais como falhas, truncamentos, pinchouts, e outras relações geométricas.

[00109] Os atributos de isocontornos são na forma de vetores e tensores de segunda ordem. As projeções destes atributos em direções e superfícies particulares podem caracterizar o volume sísmico melhor do que qualquer seu componente individual. A projeção do gradiente sobre qualquer direção é simplesmente o ponto- produt entre o vetor gradiente e um vetor unidirecional unitário. Uma tal projeção define a derivada direcional de uma função de isocontorno ao longo de uma direção particular. Quanto maior a projeção (ou derivada direcional) for, mais íngreme a mudança de amplitude ao longo daquela direção. Similarmente, a projeção Hessiana sobre um vetor direcional unitário d é a derivada de segunda ordem de uma função de isocontorno ao longo daquela direção e pode ser computada como fdd = d H d, em que ■ indica o interior, ou ponto, produto. Para o Hessiano, projeções de particular interesse estão ao longo da direção gradiente, direção de inclinação, ou direção ortogonal ao plano atravessado tanto por gradiente como por inclinação.

[00110] O gradiente do volume de inclinação de isocontorno e suas projeções também contêm informações úteis da estrutura de subsuperficie. O gradiente de inclinação quantifica variações de direção de inclinação no volume original. A projeção de gradiente de inclinação ao longo da direção de inclinação é uma medida da taxa de mudança de inclinação. Isto ajuda a detectar aspectos curvados nos refletores, tais como cristas, colinas, cavidades e vales. A projeção de gradiente de inclinação em um vetor que é normal à direção de inclinação define uma medida da taxa de convergência ou divergência do refletor.

[00111] Embora estes atributos de isocontorno (gradiente, inclinação, orientação, Hessiano, curvatura, e projeções) possam ser usados como atributos locais, os resultados podem ser afetados por artefatos em áreas contaminadas com ruído ou áreas com complexidade elevada. O aplainamento reduz tais artefatos e pode ser realizado, por exemplo, por convolução com um amplo filtro de baixa passagem ou aplicação de um filtro mediano provido de grande janela. Preferivelmente, entretanto, o aplainamento é realizado de uma maneira de estrutura-orientada através de distâncias maiores ao longo do tecido para evitar mistura através dos estratos. Por aplicação de uma longa-distância, filtro de estrutura-orientada, atributos de isocontorno locais podem ser convertidos em atributos geológicos.

[00112] 1c. Convergência

[00113] Outro novo atributo geológico da presente invenção é a convergência, que detecta onde reflexões sísmicas convergem para formar uma tendência de adelgaçamento regional de pacotes de reflexão, realçar pinchouts, cunhas, ou onlaps e downlaps estratigráficos. O atributo de convergência pode ser usado para enfatizar aspectos estratigráficos a um intérprete ou como entrada para um sistema de reconhecimento padrão assistido ou automatizado, tal como o sistema de análise de sistema de hidrocarbonetos sísmico da presente invenção. Dois métodos de computar convergência são: pelo uso de Hessiano ou pelo uso de gradiente de inclinação em direções normais ao tecido ou inclinação. Dependendo da vizinhança usada para computar o Hessiano ou gradiente, o aplainamento de longa distância de estrutura-orientada ao longo do tecido pode ser empregado usando-se convolução integral de linha, um filtro de difusão, um filtro guiado por esqueleto, ou filtragem no domínio de idade. Cada um destes utiliza os dados sísmicos originais para guiar o filtro ao longo dos refletores.

[00114] Outro método para computar o atributo de convergência é primeiro medir a espessura do refletor local, em seguida, estimar uma derivada de espessura lateral local e, finalmente, realizar uma operação de aplainamento lateral ao longo das reflexões para obter tendências de convergência regional. Se as tendências forem estimadas em seções transversais bidimensionais únicas, então a convergência pode ocorrer à esquerda ou à direita ou, alternativamente, ser chamada de convergência e divergência ao longo de uma direção de referência, por exemplo, esquerda para direita. Se as tendências forem estimadas de fatias ortogonais extraídas de um conjunto de dados 3D, então convergência pode ser associada com uma orientação e magnitude. Além disso, as tendências podem também ser estimadas em um sentido volumétrico, empregando-se janelas de análise tridimensional em vez de múltiplas análises bidimensionais ortogonais.

[00115] Os métodos para estimar espessura de refletor local incluem: duração de circuito, a distância entre extremos justapostos ou cruzamentos zero; frequência instantânea, ou a diferença de tempo entre duas superfícies justapostas do esqueleto sísmico. Um atributo contínuo, tal como frequência instantânea, pode precisar ser bloqueado. Qualquer tal atributo de espessura pode ser diferenciado lateralmente para avaliar a taxa de mudanças de espessura local. A taxa de mudanças de espessura local pode também ser estimada pela divergência de um campo vetorial, por exemplo, normais de reflexão, ou pela divergência de um campo sensor, por exemplo, tangentes de reflexão.

[00116] Para obterem-se tendências de convergência regional, derivadas de espessura local são aplainadas com um operador de filtro de longa-escala de estrutura-orientada ao longo do tecido, empregando-se convolução integral de linha, um filtro de difusão, um filtro guiado por esqueleto, filtragem do domínio de idade, ou qualquer outro método que use os dados sísmicos originais para guiar o filtro ao longo dos refletores.

[00117] Para eficiência computacional, os rótulos ou indicadores de horizonte podem ser codificados com as derivadas.

[00118] A Figura 7 mostra os atributos de convergência baseados na duração de circuito, derivada lateral, e empregando a convolução integral de linha para aumentar as tendências de convergência de longa distância. Realizando-se este procedimento tanto em linhas como em linhas cruzadas, um atributo de convergência tridimensional com magnitude e orientação pode ser construído (Fig. 8).

[00119] 1.d Confluência

[00120] Já outro novo atributo geológico da presente invenção é chamado de confluência. As reflexões sísmicas são tratadas como uma rede de caminhos transportando tráfego, e alguns destes caminhos são congestionados em razão deles estarem unidos com muitos outros. A confluência mede esta congestão, estimando como muitos trajetos passam através de qualquer voxel dado. Em termos de estratigrafia sísmica, ramificar e unir reflexões produz terminações nas formas de downlaps, onlaps, toplaps, e truncações, como representado na Fig. 9. Todas estas terminações são potencialmente associadas com desconformidades; downlaps podem também ser associados com superfícies de downlaps, isto é, superfícies de inundação ou transgressivas.

[00121] Uma família de métodos usados para computar confluência é baseada na densidade das linhas de fluxo especificadas pelas tangentes de reflexão. As linhas de fluxo são computadas por integração do campo de tangente ou simplesmente seguindo-se uma cadeia de direções de tangente discretizadas. Partindo-se, por exemplo, com o traço mais à esquerda, uma nova linha de fluxo pode ser iniciada e propagada para a direita. As novas linhas de fluxo são iniciadas em incrementos de profundidade regularmente afastados ou tempo de deslocamento de duas-direções, ou em eventos específicos, tais como picos (máximos), vales (mínimos), e/ou cruzamentos zero. As amostras em uma linha de fluxo e, preferivelmente, em uma zona pequena em torno da linha de fluxo, são marcadas para pertencer àquela linha de fluxo particular. Uma vez todas as linhas de fluxo do primeiro traço sejam iniciadas e marcadas, o algoritmo prossegue para o traço vizinho à direita e começa a iniciar e marcar linhas de fluxo adicionais, em incrementos regulares em áreas até agora não marcadas ou em eventos até agora não marcados. As linhas de fluxo podem ser traçadas e marcadas à esquerda ou à direita; preferivelmente, as linhas de fluxo são computadas e marcadas tanto à esquerda como à direita do traço atual. Uma vez que todas as áreas ou eventos não marcados foram usados para iniciar uma linha de fluxo e, assim, serem marcados, o processo se repete com o próximo traço.

[00122] Quando toda a amostra ou evento tiver sido marcado pelo menos uma vez, a confluência é determinada computando-se quantas vezes cada amostra foi marcada. Preferivelmente, este cálculo é realizado simultaneamente com a marcação de amostra/evento, simplesmente incrementando-se contadores associados nos locais de amostra. Embora estas etapas possam ser realizadas da esquerda para a direita, preferivelmente elas são repetidas da direita para a esquerda também e, em seguida, adicionadas ao resultado anterior, garantindo simetria. A marcação das amostras em uma região em tomo de uma linha de fluxo pode ser realizada com um indicador constante ou com valores que variam dependendo de sua distância da linha de fluxo. Os pesos indicadores variáveis, por exemplo, poderiam parecer-se com um triângulo ou um Gaussiano. Desprezando os problemas de quando e onde iniciar novas linhas de fluxo, poder-se-ia simplesmente marcar e contar as linhas de fluxo sem marcar quaisquer regiões e, então, em seguida, realizar uma operação de aplainamento, por exemplo, com um vagonete, um triângulo ou um filtro Gaussiano.

[00123] Em vez de traçar linhas de fluxo de todos os sentidos, da esquerda para direita (e vice-versa), as linhas de fluxo podem ser computadas para somente uma distância finita do ponto de iniciação. Esta computação pode ser realizada por convolução integral de linha ou, preferivelmente, uma sua variação que simplesmente incremente os contadores das amostras tocadas pelo estêncil.

[00124] Generalizações destas formas de realização baseadas em linha de fluxo de confluência para dados tridimensionais podem ser obtidas realizando-se computações independentemente ao longo das direções em linha e em linha cruzada, ou ao longo de um conjunto de direções arbitrárias. Alternativamente, linhas de fluxo (curvas) podem ser generalizadas a lâminas de fluxo (superfícies), ou integrais de linha podem ser generalizadas a integrais de superfície.

[00125] Outras formas de realização de confluência podem ser obtidas usando-se esqueletos sísmicos. O esqueleto bruto é representado por um gráfico de eventos direcionado (picos, vales, e cruzamentos zero: vértices) e uma correlação de forma de onda com eventos em traços vizinhos (bordas). De um vértice particular (evento), o gráfico pode ser atravessado para a esquerda e/ou para a direita, seguindo-se as bordas (correlações). Iniciando-se as travessias de gráfico de diferentes vértices, permite-se contar quantas vezes qualquer vértice fornecido obteve passando através. Todo vértice é usado para iniciar uma travessia ou, preferivelmente, somente vértices que não foram passados através anteriormente são usados para iniciar uma nova travessia. No caso preferido, primeiro, todo o vértice iniciará uma travessia, porém progressivamente, vértices terão passado através pelo menos uma vez e o início de novas travessias tornar-se-á aumentadamente raro. Esta forma de realização particular de travessia-de-gráfico de confluência pode também ser baseada em esqueletos morfológicos.

[00126] Duas outras formas de realização de confluência são baseadas no esqueleto topológico. Primeiro, uma pode contar quantas vezes diferentes superfícies sobrepõem ou subpõem qualquer superfície dada. Se uma superfície abaixo de outra superfície particular terminar, então uma nova superfície deve existir 3M57 abaixo da superfície original (sem-terminação). O score de quantas diferentes superfícies estão acima e abaixo de uma dada superfície, portanto, aproxima-se do número de terminações junto a ela. Por sua natureza, grandes superfícies tendem a ter altos scores. Normalizar os scores por tamanho de superfície, entretanto, tende a designar altos scores relativos para as superfícies mais diminutas, em razão de um pequeno número de scores ser dividido por uma pequena área ou um pequeno número de amostras formar a superfície. Assim, os scores precisam ser normalizados. Um método de normalização envolve usar a área de superfície aumentada de alguma capacidade, por exemplo, uma metade. Uma normalização preferida, entretanto, é baseada na análise da significância de encontrar um certo número de terminações em uma superfície de um tamanho particular. Para uma superfície pequena, por exemplo, há uma chance de encontrar um grande número de terminações em relação ao tamanho de superfície e, assim, este score não pode ser estatisticamente significativo. Para uma grande superfície, contudo, mesmo um score de pequena terminação pode ser acima da média de uma maneira estatisticamente significativa. A probabilidade de terminação de uma superfície pode ser estimada calculando-se a média dos scores de terminação individuais relativas, alguma sua versão ponderada, ou do número total de superfícies divididas pela área total destas superfícies. Uma vez que uma probabilidade de terminação (ou densidade) esperada é calculada, um teste binomial pode então ser usado para computar a significância estatística de um desvio do número de terminações esperado para uma superfície de um dado tamanho. Para grandes superfícies, a distribuição binomial é bem aproximada por convenientes distribuições contínuas, e esta pode ser usadas como a base para testes alternativos que são muito mais rápidos para computar, isto é, teste do x ao quadrado de Pearson e o teste-G. Entretanto, para pequenas amostras, estas aproximações se esgotam e não há alternativa para o teste binomial.

[00127] Outra forma de realização de confluência derivada de esqueletização topológica é baseada no volume de mapeamento de idade associado com as superfícies. Áreas de elevada confluência têm uma elevada densidade de linhas de fluxo ou lâminas de fluxo. Um isocontorno ou uma isosuperfície do volume de mapeamento de idade é similar a uma linha ou uma lâmina de fluxo derivada de tangentes de reflexão. Assim, uma alta densidade de isosuperfícies indica uma área de alta confluência. Uma alternativa preferida para a construção de isosuperfícies e estimativa de sua densidade é a computação de um gradiente ou derivada vertical do volume de mapeamento de idade. As áreas com elevados gradientes ou elevadas derivadas têm alta confluência.

[00128] Quaisquer formas de realização de confluência podem ainda ser homogeneizadas por aplicação de um aplainador de longa distância, por exemplo, uma integral de convolução de linha.

[00129] 1.e Pontos de Fechamento e Derramamento

[00130] Outro conjunto de atributos geológicos, descrito na Publicação do Pedido de Patente U.S. No. 2010/0149917, “Method for Geophyisical and Geological Interpretation of Seismic Volumes in Depth, Time, and Age”, de Imhof et al, é o envoltório e pontos de derramamento. Em práticas de interpretação tradicionais, um horizonte é analisado para determinar a existência de contornos fechados envolvendo um alto sítio topográfico, formando um fechamento, que poderia aprisionar hidrocarbonetos. Nota-se que o fechamento poderia ser implícito, por exemplo, onde as linhas de contorno terminam em uma falha ou contra um domo salino e são implicitamente admitidos seguir estes limites. Para cada contorno fechado, envoltórios especificam a área contida dentro. Assim, qualquer local em uma superfície pode ser analisado para determinar se ou não se situa em um contorno fechado rodeando um elevado sítio e, se assim, a área cercada pode ser computada. Para cada elevado sítio, seu fechamento máximo especifica a extensão máxima de trape de hidrocarboneto potencial.

[00131] Os fechamentos e envoltórios podem ser determinados por qualquer horizonte único ou por qualquer conjunto de horizontes, por exemplo, superfícies mapeadas por meios tradicionais ou por esqueletização automática. Estimando-se o envoltório para cada superfície em um conjunto de dados do volume de idade, ou preferivelmente, do volume de profundidade, um volume de envoltório pode ser computado que designe a cada amostra um valor de zero, se não partir de um contorno de profundidade fechado que circula um elevado sítio. De outro modo, a área, da área circundada, é designada para o local de amostra.

[00132] Usando o volume de idade, uma superfície pode ser construída por seleção de uma idade, ou um local, ou alguns outros critérios. Esta superfície pode ser analisada para elevados sítios, contornos que circulam os elevados sítios, e suas áreas. Preferivelmente, entretanto, o volume de profundidade é empregado. Cada fatia horizontal em um volume de profundidade representa uma superfície em alguma idade. De fato, um volume de profundidade não é mais do que uma pilha de superfícies classificadas por idade. Para cada fatia (ou superfície de alguma idade), seus valores correspondem à profundidade e, assim, cada fatia constitui um mapa de profundidade. Cada fatia pode ser analisada para elevados sítios, contornos fechados, e áreas envolvidas que permitem a computação de um volume envoltório no domínio de idade. Se desejado, o volume de envoltório inteiro ou uma sua parte, uma fatia, por exemplo, pode ser facilmente transformado em domínio de profundidade.

[00133] Uma extensão de volumes de fechamento ou envoltório são pontos de derramamento e volumes de pontos de derramamento. Os pontos de derramamento são um local próximo ao contorno de fechamento máximo onde os contornos são rompidos e, assim, o trape potencial vaza (Fig. 1). A determinação de pontos de derramamento permite a criação de um volume de pontos de derramamento e exame de como diferentes trapes potenciais derramam e se alimentam entre si. Os locais e número de pontos de derramamento podem ser usados para arriscar um prospecto ou para guiar uma investigação quanto às regiões onde hidrocarbonetos vazaram de e acumularam-se em. Embora pontos de derramamento possam ser identificados a partir de volumes de idade por extração de superfícies de isoidades e seu exame, os pontos de derramamento e volumes de pontos de derramamento são preferivelmente gerados de volumes de profundidade que correspondem a pilhas de mapas de profundidade de diferentes idades (isto é, pilhas de mapas de profundidade para diferentes horizontes). Se necessário, os pontos de derramamento e volumes de pontos de derramamento, determinados pelos volumes de profundidade, podem ser facilmente transformados em domínio de profundidade, empregando-se o volume de profundidade como uma tabela de consulta ou por interpolação.

[00134] 1.f Resíduos de Fase

[00135] Outros atributos geológicos podem ser formados de anomalias de fase instantânea, cujos pontos de fora se localizam onde a ondinha sísmica se divide (Publicação do Pedido de Patente PCT WO 2009/137150 Al, “Method For Geophysical And Stratigraphic Interpretation Using Waveform Anomalies”, de Imhof). O atributo de anomalia de fase instantânea é um atributo diferente indicando se ou não uma divisão de onda pequena ocorre em um local de amostra. A convolução tradicional com um filtro de pulsos longos permite a computação de uma densidade local de divisões de ondinha. O aplainamento de estrutura-orientada de longa distância, por exemplo, com uma convolução integral de linha, permite a computação de densidade de divisões de ondinha, que é consistente com o tecido sísmico subjacente. Preferivelmente, os circuitos e cordões de divisão de ondinha são primeiro classificados com base em seu formato e orientação para suprimir divisões de ondinha causadas por ruído ou artefatos.

[00136] Os resíduos de fase ocorrem em locais onde formas de onda se dividem, o que implica em que um evento de reflexão adicional comece ou que uma reflexão incorpore-se em outra. Em outras palavras, uma terminação ocorreu. As terminações são com frequência classificadas como onlap, downlap, toplap e truncação (Fig. 9), dependendo das relações entre as reflexões. Assim, os resíduos de fase podem ser classificados como onlaps, downlaps, toplaps, e truncações. Como uma alternativa para resíduos de fase, as terminações podem também ser detectadas e classificadas de esqueletos sísmicos. A menos que a informação de inclinação seja usada, downlaps e onlaps não podem ser distinguidos e podem precisar ser agrupados entre si. Sem informação de inclinação, toplap e truncação não podem ser distinguidos e precisam ser agrupados entre si. As terminações e seu score podem ser afixados em superfícies ou linhas de fluxo de confluência. Preferivelmente, as terminações e seus scores são convertidos em atributos de densidades ou geológicos por aplanamento. Os atributos geológicos particularmente úteis baseados em terminações são densidades de terminação, tais como a densidade onlap/downlap, a densidade truncação/toplap, ou a diferença entre a densidade onlap/downlap e a densidade de truncação/toplap que indica qual conduta domina.

[00137] 1.g Textura

[00138] Outros atributos geológicos são baseados em textura sísmica. A Patente U.S. No. 6.438.493 B I, “Method for seismic facies interpretation using textural analysis and neural networks”, de West e May, descreve um método, para a identificação de fácies sísmicas, baseado em atributos têxteis computados com matriz de coocorrência de nível-cinza (GLCM). Os atributos de textura locais derivados de GLCM incluem, mas não são restritos a: textura de homogeneidade, inércia (também conhecida como o momento ou contraste de diferença de elemento), entropia, e energia (também conhecida como uniformidade). Empregando-se aplanamento de longa distância ao longo do tecido, estes atributos podem ser convertidos em atributos geológicos regionais relacionados à textura de fácies sísmicas. Os atributos de textura baseados em GLCM podem ainda ser generalizados em atributos de textura tridimensionais (por exemplo, Patente U.S. No. 6.226.596 B1), “Method for analyzing and classifying three dimensional seismic information”, de Gao).

[00139] Outro atributo de textura (Publicação do Pedido de Patente PCT WO 2009/011735 “Geologic Features From Curvelet Based Seismic Attributes”, de Neelamani ande Converse) identifica aspectos estratigráficos de dados sísmicos, coletando-se uma transformação de dados de uma curva pequena. Desta representação de curva pequena, atributos de dados geofísicos selecionados e suas interdependências são extraídos e usados para identificar aspectos geológicos. Empregando-se aplainamento de longa distância ao longo do tecido, estes atributos podem ser convertidos em atributos geológicos regionais relacionados à textura e fácies sísmicas.

[00140] Já outros atributos de textura são baseados em medidas Fourier multidimensionais providas de janela (Publicação do Pedido de Patente PCT WO 2010/053618, “Method for Seismic Interpretation Using Seismic Texture Attributes”, de Imhof). Uma tal medida particular é a regularidade, um atributo de textura que mede como dados sísmicos providos de faixa (ou regulares) parecem ser. As áreas livres, reflexões atravessantes, exibem um alto grau de regularidade, enquanto áreas de ruído com reflexões desorganizadas exibem baixa regularidade. Com janelas de pequena análise, a regularidade é uma medida de descontinuidade. Com grandes janelas de análise, a regularidade atua como uma medida de desordem. Assim, a regularidade pode também ser designada caos de descontinuidade de multitraços. Usando-se um aplanamento de longa distância ao longo do tecido, estes atributos Fourier multidimensionais providos de janelas podem ser convertidos em atributos geológicos regionais relacionados à textura e fácies sísmicas. Por exemplo, reflexões regionais atravessantes, comumente denominadas rastros de ferrovia, com frequência correspondem às superfícies transgressivas (ou superfícies de inundação) que frequentemente formam vedações. Assim, o aplanamento de longa distância de regularidade (ou seu inverso, caos) cria um atributo geológico relacionado à vedações de hidrocarbonetos e, em razão dos xistos, formam tanto vedações como atuam como rochas matriz, rochas matriz de hidrocarboneto.

[00141] Randen e Sonneland (“Atlas of 3D Seismic Attributes”, in Mathematical Methods and Modeling in Hydrocarbon Exploration and Production, Iske e Randen (editores), Springer, págs. 23-46 (2005)) apresentam um resumo dos atributos sísmicos tridimensionais adicionais que caracterizam a textura sísmica ou aspectos sismo-estratigráficos.

[00142] 1.h Anomalias WPCA

[00143] Outra família de atributos relevantes para esta invenção é descrita na Publicação do Pedido de Patente PCT WO 2010/056424 “Windowed Statistical Analysis for Anomaly Detection in Geophysical Datasets”, de Kumaran et al. Estes atributos realçam locais em um ou múltiplos conjuntos de dados onde os dados sísmicos são estatisticamente anômalos em comparação com outros locais. Além disso, estes atributos geram uma base de padrões vetoriais variando a partir da maioria comum para a maioria anômala. Projetar os dados em uma combinação linear destes padrões enfatiza alguns padrões, enquanto suprimindo outros. Por exemplo, o enfaixamento sísmico ubíquo pode ser suprimido. Empregando-se aplainamento de longa distância ao longo do tecido, o filtrado sísmico pode ser convertido em atributos geológicos regionais. Outra aplicação é a decomposição local dos dados nestes padrões e a determinação do padrão dominante, cuja score de efeitos das amostras sísmicas baseia-se nos padrões da amostra circundante. Empregando-se um filtro mediano ou de modo de longa distância ao longo do tecido, a classe local pode ser convertida a uma regional.

[00144] 1.i Outros Atributos

[00145] Outro tipo de atributo é baseado na inversão de dados sísmicos multideslocados e o comportamento observado de amplitude versus deslocamento que permite a predição de porosidade e teor de argila, e permite o score de litofácies, (por exemplo, Patente U.S. No. 7.424.367 B2, Ata da conferência de 2008 sobre “Method for predicting lithology and porosity from seismic reflection data”, de Saltzer et al., “Seismic Rock-Property Inversion and Lithofacies Prediction at Erha Field, Nigeria”, Xu et al, Nigerian Association of Petroleum Explorationists (NAPE); “Lithofacies Prediction in Deep Water Water Reservoirs”, Oppert et al, Society of Exploration Geophysicists, Expanded Abstracts, 1708-1711, (2006)). Empregando-se aplanamento de longa distância ao longo do tecido, estes atributos podem ser convertidos em atributos geológicos regionais relacionados à litofácies, teor de argila, e porosidade.

[00146] A Patente Européia No. EP 1110103 Bl “Method Of Seismic Signal Processing”, tode Meldahl et al., descreve um método para gerar atributos que permitem a detecção de chaminés de gás potenciais. Uma chaminé de gás é um distúrbio vertical de resposta sísmica de filtração de gás que degrada os dados sísmicos devido a limitações de aquisição e processamento. O atributo de cubo de chaminé resultante realça distúrbios verticais de sinais sísmicos que são com frequência associados com chaminés de gás. Ele revela informações da história dos hidrocarbonetos e fluxo de fluido. Em outras palavras, o cubo de chaminé pode revelar onde hidrocarbonetos originaram-se, como eles migraram para dentro de um prospecto e como eles deslocam-se destes prospectos. Como tal, um cubo de chaminé pode ser visto como uma nova ferramenta indicadora de hidrocarbonetos indiretos.

[00147] Além disso, Loseth et al. (“Hydrocarbon leakage interpreted on seismic data, “Marine and Petroleum Geology 26(7), 1304-1319, (2009)) provê uma revisão de interpretação de escapamento de hidrocarbonetos sísmico.

[00148] Finalmente, os atributos sísmicos podem ser relacionados com as propriedades da subsuperfície que podem ser usadas para simular fenômenos de transporte, tais como fluxo de calor e temperatura, que afetam a maturação de hidrocarbonetos, ou permeabilidade e fluxo de fluido, que afetam a migração de hidrocarbonetos (Publicação do Pedido de Patente PCT WO 2009/137228 A2, “Transport Property Data Calculated From Derivative Seismic Rock Property Data For Transport Modeling”, de Oppert et al).

[00149] 1 .J Atributos de Combinação I Atributos Play

[00150] Em algumas formas de realização do presente método inventivo, os dados sísmicos são examinados para um ou múltiplos plays específicos, em vez de elementos de um sistema de hidrocarbonetos mais genérico. Exemplos poderiam incluir: plays anticlinais (Fig. 10), plays de falha-normal (Fig. 11), plays de flanco de domo-salino (Fig. 12), plays de canal de cordão de sapata, ou outros, tais como ilustrados nas Figs. 13 a 16. Nestes casos, atributos, detectores ou fluxos de trabalho são necessários para distinguir aspectos específicos, tais como falhas, sal, ou canais.

[00151] A Publicação do Pedido de Patente PCT WO 2009/082545 “Detection of Features in Seismic Images”, de Kumaran e Wang, descreve tal método para a detecção de canais ou falhas em dados sísmicos. Para detecção de falhas, bordas são identificadas em imagens sísmicas aplainadas e as intensidades da borda são integradas em múltiplas direções, por exemplo, empregando transformação de Radon para detectar a presença e a orientação de linhas de falha. Para canais, as bordas são detectadas e convertidas em curvas suaves para identificar bordas de canal. Conjuntos de curvas suaves paralelas são então examinados para encontrar pares de curvas que correspondem às bordas de canal à esquerda e à direita, assim definindo o canal.

[00152] A Patente U.S. No. 7,203,342 B2, “Image Feature Extraction,”, de Pedersen, descreve outro tal método designado para extrair falhas de dados de atributo sísmico, embora possa também ser usado para outros problemas de extração de linha ou superfície também. O método é baseado em rastreio de antena, onde antenas ou agentes numéricos vagueiam através de uma descontinuidade ou volume de detecção de borda, lentamente conectando peças próximo de bordas alinhadas.

[00153] Um último exemplo é a detecção de corpos salinos empregando-se atributos de regularidade ou caos descritos na Publicação do Pedido de Patente PCT WO 2010/053618 “Method for Seismic Interpretation Using Seismic Texture Attributes”, de Imhof; ou Randen e Sonneland, “Atlas of 3D Seismic Attributes” (in Mathematical Methods and Modeling in Hydrocarbon Exploration and Production, Iske e Randen (editores), Springer, págs. 23-46 (2005)).

[00154] Estes são apenas três métodos exemplares para identificar aspectos específicos. Muitos outros foram descritos e são bem conhecidos a praticantes da arte.

[00155] Partição

[00156] A fim de detectar a presença simultânea de múltiplos elementos proximais do sistema de hidrocarbonetos ou múltiplos elementos de play, o volume de dados pode opcionalmente ser dividido em pelo menos um segmento para análise, talvez em combinação com um segmento de segundo plano que não seja analisado. Modos típicos de partição são a análise de: voxeis individuais, pequenos blocos, blocos de tecido-alinhados, camadas, ou grupos de voxeis contíguos. A partição não precisa ser mutuamente exclusiva. Divisões individuais podem sobrepor-se.

[00157] A partição mais simples é a análise de voxel-por-voxel, porém os resultados podem ser feitos por remendos, porque raramente todos os elementos são reconhecidos no mesmo voxel e é improvável encontrar uma região estendida de voxeis contíguos com todos os elementos presentes. Tais técnicas, tais como aplainar atributos ou elementos detectados estendidos em voxeis circundantes, podem ser empregadas na etapa 5 do presente método inventivo (Fig. 2), avaliação, para criar prospectos contíguos maiores. A partição em voxeis é um método de divisão preferido, em razão dos dados serem com frequência prontamente representados como voxeis. Por esta razão, a etapa de partição é considerada opcional. Além disso, partições (e assim voxeis) podem ser agregadas em maiores regiões durante a etapa de avaliação. Assim, para fins de ensinamento do método inventivo, se nenhuma divisão tiver sido explicitamente realizada, então cada voxel é considerado formar sua própria divisão.

[00158] Outro esquema de partição envolve interromper o volume de dados em blocos ou tijolos Cartesianos regulares, por exemplo, amostras de tamanho 20x20x20, em que a expectativa é que algum ou todos os elementos estejam presentes dentro de um bloco prospective.

[00159] Tijolos Cartesianos cortarão através dos estratos e tecido. Uma alternativa para regular tijolos ou blocos Cartesianos é alinhar os blocos com o tecido. Neste esquema, haverá diferenças de tamanho e formato entre as divisões em razão de conform arem-se ao tecido. Um modo particular de gerar tal divisão é usar o próprio tecido para definir uma estrutura de camada grossa. Especificamente, por exemplo, uma seleção de partes de superfície criada por esqueletização sísmica poderia ser estendida verticalmente para criar corpos ou segmentos.

[00160] Outro método preferido para divisão é baseado em um ou múltiplos atributos sísmicos como gerados, por exemplo, na etapa 2. As partições são criadas por limiares dos atributos, seguidos por análise de componente conectado ou similar para gerar regiões contíguas embutidas em um segundo plano. Este processo é também chamado de detecção semente (multivolumes). Um atributo preferido para controlar a divisão é a saliência, um atributo que realça locais em um ou múltiplos conjuntos de dados em que os dados sísmicos são estatisticamente anômalos em comparação com outros locais. Os atributos de saliência são descritos na Publicação do Pedido de Patente PCT WO 2010/056424 “Windowed Statistical Analysis for Anomaly Detection in Geophysical Datasets”, de Kumaran et al. A saliência é um nome genérico para um atributo que realça anomalias estatísticas em dados. Kumaran et al. descreve “inverse covariance”, “(WPCA) residue”, ou “(WPCA) anomaly,” terminologia para formas de realização específicas de saliência. Atributos de saliência adicionais podem ser computados empregando-se métodos de Kumanan. A saliência pode ser tomada limiar para realizar a análise somente de regiões anômalas contíguas.

[00161] Um método alternativo de usar atributos para divisão de controle é primeiro realizar um método de divisão diferente, por exemplo, blocos Cartesianos ou corpos de tecido alinhados e, em seguida, preservar somente divisões onde um atributo computado para cada divisão exceder um valor absoluto ou um valor relativo, ou satisfizer alguma condição prescrita, por exemplo, onde a saliência máxima é no quintil de topo ou onde a regularidade média situa-se dentro da faixa entre 0,3 e 0.6.

[00162] Um último exemplo de um método de divisão envolve subdividir o volume de dados em pequenas regiões que são classificadas por um(ns) atributo(s), por exemplo, saliência. Para muitas formas de realização da presente invenção, a ordem em que alguém avalia a prospectividade das diferentes divisões ou voxeis individuais não importa. Este último exemplo demonstra o uso de um atributo, tal como saliência (ou tamanho, outro atributo, ou uma combinação de outros atributos) para determinar a ordem em que as divisões são analisadas. Usar uma tal priorização pode não ser necessário para analisar todas as regiões divididas. Esta forma de realização do presente método inventivo começa a avaliação do sistema de hidrocarbonetos (etapa 5) com as regiões mais salientes e continua com progressivamente menos regiões salientes até um número prescrito de regiões ter sido analisado, um número prescrito de prospectos ter sido constatado, ou um tempo prescrito permitido para análise ter sido excedido.

[00163] Scores de Elemento

[00164] O score pode incluir normalização, condicionamento, combinação, ou 42157 escalação. Para alguns elementos do sistema play ou de hidrocarbonetos, esta etapa admite pelo menos um atributo por elemento disponível para formar um score para cada partição (ou no caso mais simples, para cada voxel) que expressa probabilidade ou expectativa de um elemento particular ser contido em uma divisão particular. O score é suposto uma etapa opcional, porque ele pode ser realizado como uma etapa explícita separada 4, como mostrado no fluxograma da Fig. 2. Em algumas formas de realização da invenção, entretanto, ele pode ser realizade em combinação com a formação de atributo (etapa 2) ou avaliação de prospecto (etapa 5).

[00165] Na forma mais simples do método inventivo, um atributo é diretamente usado como um score para indicar se ou não um elemento particular está presente em qualquer dado local ou não. Diferentes atributos, entretanto, podem ter variações para seus valores e são com frequência vantajosos para normalizar seus valores, por exemplo, entre 0 e 1 ou 0 a 255, para mais fáceis comparações entre atributos ou para facilitar a interpretação probabilística. A Normalização é apenas uma transformação linear dos valores de atributo.

[00166] Uma alternativa para normalização é a calibração ou a aplicação de uma transformação não linear de valores de atributo que é acionada por dados, por exemplo, equalização de histograma ou transformação de histograma, ou acionada por um modelo que expressa como um atributo relaciona-se com a probabilidade para um dado elemento existir. Um tal modelo pode basear-se na teoria, medições, conhecimento anterior ou contextuai, experiência, ou intuição.

[00167] Os atributos sísmicos são observações, medições, ou computações realizadas em dados sísmicos. Eles podem relacionar-se com elementos do sistema de hidrocarbonetos ou elementos de play, porém, com frequência eles não medem realmente estes elementos. Além disso, eles podem não ser únicos. Diferentes atributos podem referir-se ao mesmo elemento. O mesmo atributo pode relacionar- se com múltiplos elementos, enquanto também afetados por ruído, aquisição, processamento de dados, e com o algoritmo e parâmetros usados para gerar o atributo. O mesmo tipo de atributo pode ser computado com diferentes algoritmos. A combinação de atributos permite a formação de um score que indica a probabilidade ou expectativa para a presença de um elemento em um dado local.

[00168] Se para um dado elemento, nenhum atributo direto tiver sido computado que defina um score, então o score pode precisar ser definido indiretamente, empregando-se um atributo substituto ou score, um heurístico, um conceito, ou uma expectativa anterior (tal como um valor constante de 0,1 que simplesmente indica uma chance de 10 % para o elemento). O reservatório, por exemplo, pode ser diretamente indicado por uma estimativa de relação areia-xisto e/ou uma estimativa de porosidade. Se nenhum tal atributo indicador existir, então um substituto conceituai, derivado, por exemplo, de conceitos estratigráficos de sequência, pode precisar ser usado. Um tal atributo substituto poderia ser uma combinação de um fã de piso de bacia lowstand que sobrepõe uma desconformidade (um limite de sequência) como mostrado na Fig. 17. Fãns de piso de bacia lowstand com frequência exibem boa porosidade e permeabilidade e, assim, têm o potencial para formar o elemento de reservatório. Elevando-se o nível do mar relativo cobrir-se-á esse fã com os dedos da cunha lowstand, provavelmente para consistir de matéria mais fina, menos permeável, que é coberta por uma superfície transgressiva. Mais elevação do nível do mar relativo enterrará o lowstand abaixo dos tratos de sistemas transgressivos e hightstand. Em locais distantes, o fã pode parecer ser intercalado entre uma “desconformidade” e uma superfície “transgressiva”. Assim, um corpo de rocha entre uma desconformidade e uma superfície transgressiva, localizado no lado profundo de dita desconformidade, poderia ser um fã lowstand e seria, portanto, designado um elevado score de reservatório.

[00169] Outro aspecto de designação de score é a escalação dos segmentos ou partições. Muitos atributos são definidos em cada local de voxel. Outros podem somente estar disponíveis em superfícies ou em alguns locais selecionados. Uma partição particular pode corresponder a um voxel ou a conjuntos contíguos de voxeis. A fim de designar um score para cada partição (ou no caso cada voxel), atributos podem precisar ser interpolados ou reduzidos para obter-se um score único por partição. Para uma partição que contém múltiplas amostras, um score representativo pode ser encontrado computando-se a média, número médio, ou modo; a aplicação de um procedimento de votação; ou seleção do valor de amostra mínimo, máximo ou um aleatório.

[00170] É vantojoso aumentar os scores com uma confiança. O score representa uma probabilidade de encontrar um elemento particular em um dado local. Um score, entretanto, não indica quanto verossímil esta probabilidade é. Um local onde um atributo sísmico prediz elevada porosidade pode receber o mesmo score que um local com um fã lowstand, como sugerido pela proximidade de uma desconformidade e uma superfície transgressiva. Porém, no primeiro local, há uma medição direta relativa de porosidade e, assim, a presença de reservatório. No segundo local, a existência de reservatório é inteirada por um modelo conceituai. A medição direta inspira mais elevada confiança do que uma predição baseada em um modelo conceituai. Esta diferença é expressa no valor de confiança associado com um score.

[00171] Avaliação de Prospecto

[00172] Usar os scores para os elementos do sistema play ou de hidrocarbonetos permite a avaliação dos prospectos de partição (ou voxel) (etapa 5) e a identificação de suas deficiências, isto é, elementos que são fracamente expressos, deficientes ou não resolvidos.

[00173] Muitos dos atributos geológicos descritos representam tendências e, assim, têm uma resolução relativamente fraca. Em algumas formas de realização do sistema inventivo, as necessidades para o arranjo espacial de elementos, por exemplo, rocha de selagem acima do reservatório, são relaxadas e substituídas por um teste quanto à presença de todos os elementos dentro de uma partição singular. Este teste quanto à presença de todos os elementos é chamado de um teste lógico provido de portas, em razão de uma divisão ter que ser aprovada no teste por todos os elementos a fim de ser julgada prospective. A abordagem lógica provida de portas dos elementos combinando, portanto, é uma abordagem de “tudo ou nada”. Satisfazer um elemento com frequência significa exceder um limiar de score que tipicamente é especificado pelo intérprete ou codificado no sistema. Elementos insatisfatórios podem ser indicados para outro exame porque constituem as ligações fracas de uma partição. A lógica provida de portas é uma forma de realização preferida de avaliação de prospecto.

[00174] Um teste de prospectividade mais gradual é um voto que simplesmente conta quantos elementos estão presentes em um dado voxel ou partição. Quanto mais elementos presentes, mais uma divisão é julgada prospective. Um elemento está presente ou ausente baseado em uma comparação de seu score em relação a um limiar especificado pelo intérprete ou codificado no sistema. Elementos ausentes ou fracamente expressos podem ser sinalizados. A votação é outra forma de realização preferida da avaliação de prospecto. Uma variação da votação é votação ponderada em que alguns elementos recebem um peso mais elevado ou mais votos do que outros. A votação ponderada poderia ser usada, por exemplo, para enfatizar elementos que são detectáveis com mais elevada confiança.

[00175] As combinações de lógica provida de portas com votação são possíveis, por exemplo, requerendo que um número de elementos especificado esteja presente em vez de requerer a presença de todos os elementos. Outro exemplo é a presença necessária de alguns elementos especificados e a presença desejada de outros.

[00176] A avaliação de prospecto pode incluir aspectos da etapa 4, score. Um exemplo é a combinação de diferentes atributos, por exemplo, atributos computados com diferentes algoritmos ou parametrizações, que se referem ao(s) mesmo(s) elemento(s). Em vez de primeiro combinar os atributos em scores para elementos individuais e, em seguida, avaliar o prospecto baseado nos scores, a prospectividade de partição pode ser avaliada pelos diferentes atributos, por exemplo, votando-se ou contando-se quantos atributos limiares o usuário ou sistema-especificado excedeu.

[00177] Limiares rígidos podem ser evitados pelo uso de uma lógica multiestimada, por exemplo, lógica difusa, que é derivada da teoria de conjunto nebuloso para lidar com o raciocínio, que é aproximado em vez de preciso. A lógica provida de portas é clara e binária, com valores de associação 0 e 1 representando abaixo do limiar (elemento ausente) e acima do limiar (elemento presente). Os scores de lógica difusa, em vez disso, têm valores de associação que variam de 0 a 1 e representam o grau de verdade de uma afirmação. Tanto o grau de verdade quanto o de probabilidades podem variar entre 0 e 1 e, portanto, podem parecer similares. Entretanto, eles são conceitualmente distintos. A verdade representa a associação em conjuntos vagamente definidos, não probabilidades de algum evento ou condição como na teoria de probabilidade. Obtendo, por exemplo, uma camada estratigráfica que contém 70% de sedimento e 30 % de xisto. Podemos considerar dois conceitos: reservatório e selagem. O significado de cada um deles pode ser representado por um certo conjunto nebuloso. Um poderia definir a formação como sendo 0,7 reservatório e 0,3 selagem. Observa-se que o conceito de reservatório seria subjetivo e, assim, dependeria do observador ou designer. Outro designer poderia igual e corretamente construir um conjunto de funções associadas em que a formação seria considerada reservatório se a parte de sedimento excedesse 50%. A lógica difusa usa graus de verdade como um modelo matemático do fenômeno de incerteza, enquanto a probabilidade é um modelo matemático de aleatoriedade. Um cenário probabilístico definiria primeiro uma variável escalar para a fração de sedimento e, depois, definiria distribuições condicionais descrevendo a probabilidade de alguém chamar a formação de um dado reservatório de uma específica fração de sedimento. Observa-se que o condicionamento pode ser obtido tendo-se um observador específico aleatoriamente selecionando o rótulo para a camada, uma distribuição através de observadores deterministas, ou ambos. Consequentemente, a probabilidade não tem nada em comum com a difusibilidade; estes são simplesmente diferentes conceitos que superficialmente parecem similares em razão de usarem o mesmo intervalo de números reais entre 0 e 1. Além disso, a confusão surge porque propriedades de variáveis aleatórias são análogas à propriedades de estados lógicos binários, e teoremas, tais como o de De Morgan, que se refere a operadores lógicos “e” e “ou” em termos entre si via negação, têm dupla aplicabilidade.

[00178] Uma alternativa para lógica difusa é a lógica Bayesiana, que é baseada na teoria de probabilidade Bayesiana, que possibilita raciocinar com declarações incertas. Para avaliar a probabilidade de uma hipótese ou configuração, um intérprete humano ou máquina especifica alguma probabilidade anterior, que é então atualizada na luz de novos dados relevantes. A interpretação Bayesiana provê um conjunto padrão de procedimentos e fórmulas para realizar este cálculo. Um método para integrar as necessidades de uma configuração particular e para determinar sua probabilidade é usar uma Bayesian Belief Network (BBN). Uma rede Bayesiana é um modelo baseado em gráfico probabilístico que representa um conjunto de variáveis aleatórias e suas dependências condicionais, via um gráfico acíclico direcionado (DAG). Este gráfico indica a estrutura de independência condicional entre variáveis aleatórias representando os elementos diferentes e seus arranjos espaciais. Por exemplo, uma rede Bayesiana poderia representar as relações probabilísticas entre elementos de play e tipos de play. Dado certos elementos de play, a rede poderia ser usada para computar a probabilidade da presença de vários tipos de play. Formalmente, redes Bayesianas são gráficos acíclicos direcionados cujos nodos representam variáveis aleatórias no sentido Bayesiano: podem ser quantidades observáveis, variáveis latentes, parâmetros desconhecidos ou hipóteses. As bordas representam dependências condicionais; nodos que não estão conectados representam variáveis que são condicionalmente independentes entre si. A Figura 18 apresenta um gráfico de rede exemplar para o sistema de hidrocarbonetos inteiro. Os scores para o reservatório e selagem são combinados por análise de sua sobreposição, que indica se o reservatório está coberto pela rocha de selagem. Dada a geometria local, o sistema reservatório/selagem é designado um score de trape. Preferivelmente, o trape é enchido com hidrocarbonetos e, assim, o score de trape é então combinado com scores para a presença de uma fonte e indicações de pelo menos o potencial para trajetos de migração que se direcionam da fonte para o trape. De forma ideal, há indicação direta para hidrocarbonetos (DHI) nos dados sísmicos também. Portanto, os scores para trape, fonte, migração, e DHI são todos combinados para produzir um score de acumulação. Todas os scores de acumulação podem agora ser ponderados por tamanho e confiança para arriscar e classificar alvos de hidrocarbonetos potenciais.

[00179] Mesmo quando empregando atributos geológicos que representem tendências, as partições podem ser muito pequenas para conter todos os elementos, por exemplo, quando se dividindo em voxeis individuais. Duas abordagens para testar a coexistência ou colocação dentro de pequenas partições ou entre voxeis são mudanças espaciais relativas e extensão das regiões de influência. A Mudança de scores converte espacialmente um teste para uma dada relação espacial em um teste para colocação. Por exemplo, o score de selagem pode ser mudado descendentemente em relação ao score do reservatório para examinar a existência de rocha de selagem através do reservatório. Sem a aplicação de tal mudança, precisa-se testar quanto à selagem sobre o reservatório. Com aplicação de tal mudança, pode-se examinar a presença colocada de selagem e reservatório, por exemplo, com uma lógica provida de portas. Em vez de usar uma única mudança, pode ser vantajoso avaliar a coexistência de duas mudanças para uma faixa de modificações, dando preferência a menores mudanças; por exemplo, designar uma confiança mais elevada para menores mudanças para indicar a mais estreita proximidade dos elementos requeridos ou reduzir a prospectividade como uma função de distância. As mudanças laterais podem ser usadas para detectar elementos de play, tais como falhas, próximo ao reservatório. Em vez de aplicar mudanças verticais ou laterais, as mudanças podem ser aplicadas ao longo do tecido ou perpendicular ao tecido. Além disso, podem também ser aplicadas em direções arbitrárias.

[00180] A mudança é ineficiente se muitas diferentes variações e orientações ou direções precisarem ser examinadas. Uma alternativa preferida para a mudança é estender a região de influência para scores brutos ou limiares, por exemplo, por convolução ou dilatação morfológica. No primeiro caso, scores são prejudicados ou estendidos ao longo de direções especificadas, por exemplo, por convolução com um filtro direcional afunilado que decai com a faixa aumentando, codificando a diminuição de confiança. Em uma forma de realização particular da extensão de região de influência particular por convolução, o resultado da convolução é adicionado ao score original para formar um score atualizado que é usado para avaliação de prospecto. A direção de extensão é especificada por intérprete ou determinada pelo tecido.

[00181] Scores brutos ou limiares podem ser estendidos por dilatação morfológica direcional que arrasta elevados scores ao longo de uma direção especificada em áreas com scores inferiores. A direção pode ser especificada pelo intérprete, ser codificada no sistema, ou ser derivada do tecido. Operações morfológicas podem ser aplicadas em scores brutos ou limiares e partições arbitrárias.

[00182] A Figura 19 apresenta uma aplicação esquemática do método inventivo para quatro partições ou áreas alvo potenciais. Admitindo-se que nem o score da fonte nem o score da migração podem ser determinados pelos dados fornecidos que não direcionem indicações de hidrocarbonetos sendo observadas, e que os elementos remanescentes sejam independentes, a prospectividade para cada alvo poderia ser computada multiplicando-se confiança com tamanho e a média dos scores de reservatório, rochas de selagem, e trape. A tabela 1 apresenta um exemplo deste processo de score.

[00183] Tabela 1: classificação exemplar para os alvos da Fig. 19, de tamanho, confiança, e os scores médios para o reservatório, rocha de selagem e trape.

[00184] Em vez de representar os scores para cada elemento com valores únicos, como empregado no exemplo da Tabela 1, os scores podem ser representados por distribuições que capturam os scores medidos e suas incertezas e/ou confianças. Os scores elevadamente certos associados com alta confiança têm uma distribuição em pico com menos scores certos tendo uma ampla distribuição. Os scores que são necessários, mas não medidos por algum atributo, terão uma distribuição uniforme. Em alguns casos, as distribuições podem ser propagadas e combinadas através do sistema (Fig. 18). Os scores poderiam também ser integrados em um modelo de Monte Cario em que scores individuais são aleatoriamente retirados das correspondentes distribuições e propagados através do sistema. Repetir estas etapas muitas vezes permite a formação de subsequentes distribuições.

[00185] Análise e Visualização

[00186] A análise e visualização (etapa 6) é uma etapa opcional que pode ser combinada com a etapa 5, avaliação de prospecto. Os métodos de análise incluem: combinação de partições vizinhas (repartição), classificação de prospectos, validação dos prospectos primários, análise de prospectos secundários, e análise de ligação mais fraca.

[00187] As partições vizinhas podem ser combinadas para formar partições contíguas maiores. Um primeiro exemplo é a combinação de partições de voxel- único em corpos contíguos maiores, por exemplo, por análise de componente conectado ou detecção semente de multivolumes.

[00188] A invenção pode usar milhões de voxeis ou dezenas a milhões de partições que contêm pelo menos alguns elementos do sistema de hidrocarbonetos ou pelo menos alguns elementos de play. Preferivelmente, os prospectos em que todos os elementos estão presentes são classificados a fim de facilitar a validação do prospecto. Preferivelmente, partições que apresentam pelo menos alguns elementos são classificadas a fim de confirmar ou refutar suas faltas. A classificação pode ser realizada, por exemplo, pelo tamanho ou confiança.

[00189] As grandes partições que contém todos os elementos são de interesse primário para outra avaliação, em razão delas representarem prospectos, isto é, áreas em que o sistema inventivo prediz a existência de hidrocarbonetos, e podem finalmente ser recomendadas para perfuração. A justificativa para perfurar um prospecto será feita por outras análises tradicionais que são especialmente focalizadas nos elementos de menor confiança. Também de interesse são as partições que carecem do menor número de elementos, por exemplo, em razão dos baixos valores de confiança. Tais partições podem ser analisadas com métodos tradicionais para ver se podem ser reclassificadas como prospectos ou não. No caso, outra investigação é preferivelmente focalizada na ligação mais fraca de um sistema de hidrocarbonetos potencial ou prospecto associado com uma partição.

[00190] Prospectividade, confiança, scores, elementos sinalizados, e atributos, são tipicamente armazenados na memória ou no disco para mais análise e visualização.

[00191] Seleção de Play

[00192] Em algumas formas de realização do método inventivo, os dados são examinados não para sistemas de hidrocarbonetos genéricos, consistindo de reservatório, trape, rocha de selagem etc., porém para plays específicos (por exemplo, um play de flanco-salino) que contém elementos de play específicos. Para um play particular a ser considerado, o intérprete precisa especificar uma configuração com os elementos de definição, suas relações espaciais e apropriados atributos geológicos que se relacionam com estes elementos específicos.

[00193] As Figuras 10-16 representam alguns plays comuns. Em todas estas figuras, as rochas porosas (indicadas em 101 na Fig. 10) constituem reservatórios que potencialmente contêm hidrocarbonetos (indicados por pontos). Para evitar que os hidrocarbonetos escapem verticalmente, uma rocha de selagem é necessária (indicada por traços). Um trape ou mecanismo de trape é necessário para evitar que os hidrocarbonetos escapem lateralmente. Os plays estruturais são criados por deformação de estratos geológicos, que incluem rocha de selagem e formações de reservatório, em geometrias (ou estruturas) que permitem a acumulação de hidrocarbonetos. Tais geometrias resultantes, que envolvem reservatórios e rochas de selagem, são dominadas por dobras (Fig. 10), falhas (Fig. 11) ou diápiros salinos (Fig. 12 - onde cruzes marcam o diápiro salino). Em plays estratigráficos, a geometria de trape é formada por variações nas rochas (ou estratigrafia) que se relacionam com sua deposição. Um exemplo de plays estratigráficos é representado πa Fig. 13, onde a formação de reservatório, areia por exemplo, afina na formação de rocha de selagem, xisto por exemplo. Outros plays estratigráficos relacionam-se com a erosão dos reservatórios e formação de desconformidades que são cobertas por estratos de rocha de selagem. A Figura 14 representa um play estratigráfico clássico, onde os estratos acima da desconformidade (indicada pela linha serpeante cheia) proveem a rocha de selagem para formações de reservatório mergulhantes que foram aplainadas pela erosão. A Figura 15 representa outro play estratigráfico, em que o relevo erosional nas formações de reservatório porosas é enterrado embaixo das formações de rocha de selagem de trapes. A Figura 16 representa uma configuração em que o reservatório é formado por um calcário poroso, que é coberto por uma formação de xisto de rocha de selagem. Variações laterais em processos diagenéticos preservaram a porosidade do reservatório, porém obstruíram o espaço poroso na direção de mergulho para cima, assim evitando que os hidrocarbonetos escapem lateralmente. Os plays mostrados nas Figs. 10-16 não são de modo algum exaustivos, porém simplesmente representam exemplos. Muitos outros casos foram descritos e são conhecidos de práticos da arte (por exemplo, Hydrocarbon Traps, K.T. Biddle e C.C. Wielchowsky. The Petroleum System - From Source to Trap, AAPG Memoir 60, páginas 219 - 235 (1994)).

[00194] Para cada play, seus elementos de definição, suas relações espaciais e apropriados atributos ou scores geológicos precisam ser especificados. Para o play anticlinal, preso em dobra (Fig. 10), a exigência mínima é geométrica por natureza, isto é, a presença de estratos com estrutura anticlinal parecendo-se com um copo de cabeça para baixo. Preferivelmente, existe uma formação de reservatório com suficiente porosidade e permeabilidade para armazenar e transmitir fluidos que é coberta diretamente por uma formação de rocha de selagem, que é impermeável e forma uma barreira para evitar que estes fluidos vazem. De forma ideal, há uma indicação geofísica direta de hidrocarbonetos ou pelo menos alguma indicação de rochas fonte e trajetos de migração de fluido para dentro do trape.

[00195] Para o play preso em falha (Fig. 11), as exigências mínimas são também geométricas por natureza, isto é, a existência de estratos mergulhando para cima em uma falha e a formação de um compartimento de reservatório tridimensionalmente fechado. Se possível, uma formação pode ser identificada constituir um reservatório, enquanto outra formação em estreita proximidade acima dela é identificável como uma rocha de selagem. Idealmente, existe uma indicação direta ou indireta de hidrocarbonetos, uma fonte próxima potencial ou trajetos de migração de fluido.

[00196] Um play de flanco salino (Fig. 12) é definido por estratos mergulhando para cima dentro de um domo salino. Um play afinado estratigráfico (Figf. 13) requer uma formação de mergulho ascendente que termina em ponta. Um trape de desconformidade estratigráfica (Fig. 14) requer pelo menos formações de mergulho que são aplainadas por uma desconformidade. O play estratigráfico de relevo erosional enterrado (Fig. 15) requer pelo menos uma desconformidade com um formato em arco tridimensional. Para todas estas configurações, critérios adicionais incluem: a presença de estratos formando o reservatório coberto por estratos formando a rocha de selagem e a presença de uma fonte potencial, sinais de migração de fluido, ou mesmo indicação direta de hidrocarbonetos aprisionados.

[00197] Esquemáticos tais como aqueles representados nas Figs. 10-16 são úteis para desenvolver, definir e transmitir configurações. Os plays prospectivos podem ser encontrados e analisados em uma maneira ad-hoc com um fluxo de trabalho (ou script) específico feito sob medida ou interativamente por um intérprete. Algumas formas de realização do sistema inventivo podem ser baseadas em um ou múltiplos fluxos de trabalho que são na maior parte independentes.

[00198] Preferivelmente, os plays e atributos específicos usados para caracterizar seus elementos são armazenados em um catálogo de configuração para reutilização em situações similares.

[00199] Definição de um Catálogo de Configuração

[00200] E m vez de definir plays e sua configuração para cada aplicação da invenção novamente, pode ser vantajoso criar um catálogo ou biblioteca de plays e configurações para uso repetido. O intérprete do sistema pode então selecionar um ou múltiplos tipos de play com suas configurações associadas do catálogo. A confecção de um tal catálogo incentiva a reutilização de configurações, promove a reprodutibilidade e facilita a validação. Além disso, o catálogo permite formas de realização automatizadas do sistema inventivo, que compara os dados ou partições em relação a um conjunto de configurações potencialmente grandes tiradas do catálogo. Poderia mesmo ser desejável terem-se diferentes formas de realização do mesmo play do catálogo, para permitir variações na qualidade dos dados, ou expressão do play. Neste caso, o usuário da invenção poderia procurar um play de flanco-salino pegando todas as formas de realização de flanco-salino do catálogo, realizar a análise utilizando cada um e usar o melhor resultado de qualquer dada partição. O usuário poderia mesmo percorrer o inteiro catálogo através das partições e então designar para cada partição um tipo de play de acordo com as probabilidades.

[00201] Para análise automatizada, reconhecimento, ou comparação com um grande conjunto, ou mesmo o inteiro catálogo, as configurações são preferivelmente representadas em uma maneira formal, em vez de ad-hoc, para cada configuração. Um exemplo é uma representação gráfica para as configurações, em que vértices ou nodos correspondem aos elementos requeridos e expressam os atributos preferidos, enquanto as bordas gráficas marcam suas relações espaciais. Outra representação das entradas do catálogo é na forma de uma base de dados relacional.

[00202] As vantagens da utilização de uma representação mais formal é a reusabilidade dos componentes subjacentes, rápida modificação das configurações existentes, rápida adição de novas configurações e consistência entre os resultados. Além disso, uma representação mais formal incentiva a definição rigorosa de concepções de arranjo espacial, tais como acima ou próximo a. Por fim, o sistema pode potencialmente ser aumentado proporcionalmente em uma maneira recursiva.

[00203] Um método menos formal, porém com frequência preferido de representar configurações no catálogo, é pela designação de pesos aos atributos geológicos. Cada configuração de play é caracterizada pelos pesos atribuídos aos atributos ou scores geológicos associados. Os atributos ou scores que são importantes para os elementos particulares recebem um grande peso, enquanto os não-importantes recebem um baixo peso, zero por exemplo. As relações espaciais requeridas entre os elementos são capturadas por regiões de influência ou mudanças designadas para os atributos ou scores.

[00204] Exemplo

[00205] O exemplo é baseado em um cubo de dados sísmicos com um tamanho de 1426 por 1057 por 131 amostras (também chamadas voxeis). Três atributos geológicos são computados: regularidade, convergência e anexo. A regularidade é computada com um tamanho de janela relativamente grande de 61 por 41 por 41 amostras. A elevada regularidade demarcará os estratos planares, enquanto baixa regularidade demarcará estratos complexos com estrutura complicada ou estratigrafia complicada com erosão, entrincheiramento e enchimento. Dada a grande janela usada para a computação de regularidade, valores intermediários de regularidade correspondem a transições entre estratos planos e complexos, isto é, a zonas em que estratos planos e complexos são justapostos. Preferivelmente, portanto, associação de estratos planos com rochas de selagem e estratos complexos com reservatórios conduz à identificação de áreas com regularidade intermediária, como aquelas com justaposição de reservatório de rocha de selagem potencial.

[00206] O segundo atributo geológico usado neste exemplo é convergência. Nenhum afinamento em uma área indica que os estratos são localmente planos. Elevadas magnitudes de convergência ou afinamento indicam que os estratos estão localmente mudando sua espessura e que as superfícies sísmicas estão convergindo ou divergindo. A convergência é orientada e, assim, uma quantidade de vetor, porém o presente exemplo descarta orientação ou sinal e somente utiliza a magnitude de convergência. Áreas com convergência marcada contêm potencialmente trapes estratigráficos, por exemplo, por estratos que se afinam (por exemplo, Fig. 13) ou são chanfrados por uma desconformidade (por exemplo, Fig. 14). Assim, o atributo de convergência diferencia as áreas potenciais de coexistência de trape-reservatório.

[00207] O terceiro atributo geológico usado neste exemplo é uma inclusão. Cada voxel de volume sísmico é usado como um ponto semente para a formação de uma superfície sísmica que segue a estrutura sísmica. A superfície é examinada para determinar se o ponto semente é localizado em um contorno fechado circundando um local elevado da superfície, ou em um contorno fechado circundando um ponto baixo local ou um contorno aberto intersectando o limite de volume. No primeiro caso, o ponto tem fechamento, significando uma partícula fluida flutuante movendo- se para cima a partir deste local potencialmente ser aprisionada em um tape. No último caso, a partícula flutuante pode ascender ao longo da superfície e deixar o volume sem ser aprisionado. Em consequência, o atributo de inclusão delineia áreas de trape potenciais.

[00208] Os três atributos usados não são independentes entre si. Valores intermediários de regularidade referem-se a rocha de selagem e reservatório; elevada convergência refere-se a reservatório e trape; e inclusão somente refere-se a trape. Em vez de desmisturá-los por formação de scores separados para rocha de selagem, reservatório e trape, os atributos são limiarizados para formar scores binários quanto à presença potencial de reservatório-rocha de selagem, trape- reservatório e trape. Preferivelmente, múltiplos atributos seriam decompostos em separados scores para os diferentes elementos do sistema de hidrocarbonetos e/ou configuração de trape, por exemplo, por análise de componentes de princípio ou agrupamento. Porém dada a aplicação de costume em conjunto com a lógica de tudo-ou-nada provida de porta, os resultados desta forma de realização preferível da invenção serão muito similares em um esforço computacional muito mais elevado.

[00209] As três condições definidas serão satisfeitas em diferentes locais. Para suprimir pequenas áreas e com frequência isoladas, a rotulação do componente conectado é realizada através da intersecção dos três scores para encontrar regiões maiores, contíguas onde amostras de todas as três condições são satisfeitas simultaneamente. A Figura 20 apresenta uma fatia extraída dos três dados dimensionais movendo-se através do sistema. Para fins ilustrativos, o sistema é mostrado em um modo em cascata ou serial, onde o score de cada atributo é usado para sequencialmente destruir os voxeis não tendo um aspecto particular. Alternativamente, o sistema poderia ser implementado em um modo paralelo por intersecção simultânea de todos os scores imediatamente. 200 representam uma fatia de dados de 131 fatias. A primeira etapa é computação do tecido ou estrutura que é realizada por esqueletização do horizonte sísmico, seguido por construção de um volume de mapeamento de profundidade, 201, onde listas zebra são usadas para ilustrar o tecido detectado. O atributo de regularidade, 202, é computado e áreas com baixa e elevada regularidades são suprimidas (cinza claro), enquanto áreas com regularidade intermediária são enfatizadas (cinza escuro). Estes scores de regularidade agora-binários são usados para suprimir voxeis sísmicos com baixa ou elevada regularidade 203. O tecido sísmico, 201, é usado para computar o atributo de convergência 204 que é convertido em um score por limiarização. Suprimindo-se os voxeis sísmicos com scores de baixa convergência, destroi-se mais áreas dos dados sísmicos 205. O tecido sísmico 201 é ainda usado para computar o atributo de inclusão 206, que é então diretamente usado como um score. Suprimindo-se os voxeis sísmicos sem inclusão, reduz-se ainda mais áreas dos dados sísmicos 207. Muitos voxeis que obtêm todos três critérios são isolados ou conectados a muito poucos similares. Assim, a etapa final é considerar o tamanho das áreas conectadas restantes e suprimir todas exceto as maiores 208. Neste exemplo, somente quatro regiões permanecem, uma das quais é um reservatório de hidrocarbonetos conhecido.

[00210] O pedido de patente precedente é dirigido a formas de realização particulares da presente invenção para a finalidade de ilustrá-la. Será evidente, entretanto, para uma pessoa hábil na arte, que muitas modificações e variações nas formas de realização descritas aqui são possíveis. Todas tais modificações e variações são destinadas a situarem-se dentro do escopo da presente invenção, como definido nas reivindicações anexas. As pessoas hábeis na arte prontamente reconhecerão que em aplicações práticas da invenção, pelo menos algumas das etapas do presente método inventivo são realizadas em ou com o auxílio de um computador, isto é, a invenção é implementada por computador.

Claims

1. Método implementado por computador para analisar um volume composto de voxeis de dados (31) representando uma região de subsuperfície quanto à presença de um sistema de hidrocarbonetos ou um play (um modelo conceptual de um estilo de acumulação de hidrocarbonetos) particular, compreendendo: dividir (3) o volume de dados sísmicos para formar uma pluralidade de segmentos; e classificar a pluralidade de segmentos quanto à presença de um sistema de hidrocarbonetos ou o play particular com base pelo menos parcialmente nos scores (4) de prospectividade para os voxeis de dados sísmicos em cada segmento; em que o score (4) de prospectividade é baseado na computação de pelo menos dois atributos (2, 32) que se referem a diferentes elementos de um sistema de hidrocarbonetos ou do play particular; caracterizado pelo fato de que: o aplainamento ser aplicado aos pelo menos dois atributos (2, 32), desse modo gerando atributos geológicos (2), e o aplainamento ser aplicado através de distâncias de pelo menos dez voxeis de dados sísmicos; e um dos pelo menos dois atributos (2, 32) relacionar-se com a existência de um trape de hidrocarbonetos e a existência de um trape de hidrocarbonetos ser predita correlacionando-se espacialmente os atributos geológicos e comparando-os com um catálogo de configurações de trape de hidrocarbonetos (8).

2. Método de acordo com a reivindicação 1, caracterizado pelo fato de a partição (3) compreender combinar voxeis individuais contíguos com base em um atributo de saliência ou um score de prospectividade computado para cada voxel.

3. Método de acordo com a reivindicação 1, caracterizado pelo fato de os elementos serem selecionados de um grupo consistindo de reservatório, selagem, trape, fonte, carga, cobertura, maturação, migração, acumulação, e distribuição.

4. Método de acordo com a reivindicação 1, caracterizado pelo fato de pelo menos dois elementos serem selecionados para a computação do score (4) de prospectividade para cada voxel, e para cada elemento pelo menos um atributo sísmico (2, 32) ser selecionado e, em seguida, um score de elemento ser calculado para cada elemento selecionado com base no pelo menos um atributo sísmico (2, 32) selecionado, e o score de prospectividade ser formado combinando-se os scores do elemento.

5. Método de acordo com a reivindicação 1, caracterizado pelo fato de ainda compreender determinar dos dados sísmicos uma estrutura de camada, chamada um tecido (1), representativa da região de subsuperficie, e em que o aplainamento para gerar os atributos geológicos (2, 32) compreende integrar ou calcular por média os atributos selecionados ao longo do tecido (1) para evitar mistura de outros estratos.

6. Método de acordo com a reivindicação 5, caracterizado pelo fato de o aplainamento ao longo do tecido (1) ser realizado por filtragem de difusão anisotrópica, não linear, ou por convolução integral de linha.

7. Método de acordo com a reivindicação 6, caracterizado pelo fato de o tecido (1) ser determinado de esqueletização sísmica.

8. Método de acordo com a reivindicação 1, caracterizado pelo fato de os pelo menos dois atributos (2, 32) compreenderem pelo menos um de um grupo consistindo de atributos Hessianos, atributos de isocontorno, convergência (204), e confluência.

9. Método de acordo com a reivindicação 1, caracterizado pelo fato de os pelo menos dois atributos (2, 32) compreenderem pelo menos um de um grupo consistindo de envoltório (206) e ponto de derramamento; regularidade (202); resíduos de fase; textura; anomalias de análise de componente principal provido de janelas; atributos play; terminações; densidades de terminação; mapas de difusão; e suas combinações.

10. Método de acordo com a reivindicação 1, caracterizado pelo fato de os segmentos não se sobreporem.

11. Método de acordo com a reivindicação 1, caracterizado pelo fato de os segmentos serem classificados baseados em um único score para cada segmento gerado combinando-se os scores de prospectividade para cada voxel do segmento.

12. Método de acordo com a reivindicação 1, caracterizado pelo fato de cada voxel do volume de dados sísmico ser um segmento separado.

13. Método de acordo com a reivindicação 1, caracterizado pelo fato de os valores de atributo computados para diferentes atributos (2, 32) serem normalizados em uma escala comum.

14. Método de acordo com a reivindicação 1, caracterizado pelo fato de ainda compreender estimar um valor de confiança para pelo menos um score (4) de prospectividade.

15. Método de acordo com a reivindicação 2, caracterizado pelo fato de o atributo de saliência ser computado empregando-se análise de componente principiado provido de janelas ou mapeamento de difusão aplicado a pelo menos um atributo sísmico (2, 32).

16. Produto de programa de computador, compreendendo um computador utilizável médio tendo um código de programa legível por computador embutido nele, o dito código de programa legível por computador adaptado para ser executado para implementar um método para analisar dados sísmicos (31) representando uma região de subsuperfície quanto à presença de um sistema de hidrocarbonetos ou um play particular, onde compreende: dividir (3) o volume de dados sísmicos para formar uma pluralidade de segmentos; e classificar a pluralidade de segmentos quanto à presença de um sistema de hidrocarbonetos ou o play particular com base pelo menos em parte nos scores (4) de prospectividade para os voxeis de dados sísmicos em cada segmento; em que o score (4) de prospectividade é baseado na computação de pelo menos dois atributos (2, 32) que se referem a diferentes elementos de um sistema de hidrocarbonetos ou do play particular; caracterizado pelo fato de que o aplainamento ser aplicado aos pelo menos dois atributos (2, 32), desse modo gerando atributos geológicos (2), e o aplainamento ser aplicado através de distâncias de pelo menos dez voxeis de dados sísmicos; e um dos pelo menos dois atributos (2, 32) relacionar-se com a existência de um trape de hidrocarbonetos e a existência de um trape de hidrocarbonetos ser predita correlacionando-se espacialmente os atributos geológicos e comparando-os com um catálogo de configurações de trape de hidrocarbonetos (8).

17. Método para produzir hidrocarbonetos, caracterizado pelo fato de compreender: conduzir uma inspeção sísmica de uma região de subsuperficie; obter uma análise de prospectividade dos dados sísmicos (31) da inspeção em que a análise foi conduzida por um método como definido na reivindicação 1, que é incorporada aqui por referência; perfurar um poço na região de subsuperficie com base pelo menos em parte na análise de prospectividade, e produzir hidrocarbonetos a partir do poço.