BR112017020991B1

BR112017020991B1 - Método e sistema de processamento de dados sísmicos, e, mídia de gravação legível por computador

Info

Publication number: BR112017020991B1
Application number: BR112017020991-8A
Authority: BR
Inventors: Matteo Giboli; Cyril AGUT; Reda Baina
Original assignee: Total Sa
Priority date: 2015-03-31
Filing date: 2016-03-31
Publication date: 2022-11-16
Also published as: BR112017020991A2; EP3259619A1; US20180106919A1; US10557956B2; EP3259619B1; WO2016156530A1; MX2017012567A

Abstract

MÉTODO E SISTEMA DE PROCESSAMENTO DE DADOS SÍSMICOS ATRAVÉS DE FORNECIMENTO DE COLETAS DE IMAGENS COMUNS DE ABERTURA DE SUPERFÍCIE. O método processa entrada (DS,G[t]) incluindo, para cada um de uma pluralidade de pulsos nas respectivas localizações de fonte, traços sísmicos registrados em uma pluralidade de localizações de receptores. Dados modulados de ponto médio comum (D'S,G[t]) são também computados ao multiplicar os dados sísmicos em cada traço sísmico através de um ponto médio horizontal. Um processo de migração de profundidade é aplicado (i) aos dados sísmicos para obter um primeiro conjunto de dados migrados (MS[x,y,z]), e (ii) para dados modulados de ponto médio para obter um segundo conjunto de dados migrados (M'S[x,y,z]). Para cada pulso, valores de abertura (âS[x,y,z]) são estimados e associados com posições de subsuperfície respectivas, através de um processo de divisão aplicado ao primeiro e ao segundo conjuntos de dados migrados. Um valor migrado (RX,Y[z,a]) para uma profundidade "z" e uma abertura "a" em uma coleta de imagens comuns de abertura de superfície (CIG) em uma posição horizontal x, y é um valor migrado do primeiro conjunto de dados (...).

Description

ANTECEDENTES DA INVENÇÃO

[0001] A presente invenção se refere a métodos de processamento de dados sísmicos. Esta é aplicável no campo de prospecção geofísica quando imagens da subsuperfície são necessárias.

[0002] Para obter imagens da subsuperfície, um método sísmico é comumente utilizado, o qual consiste em criar e enviar ondas sísmicas no solo utilizando fontes tais como veículos explosivos ou vibratórios em terra, ou armas pneumáticas em regiões offshore. As ondas sísmicas penetram o solo e são pulsadas, ou são refletidas para fora em descontinuidades geológicas na subsuperfície. Como um resultado, estas voltam à superfície onde são registradas através do uso de arranjos de três componentes de geofones (em terra) ou hidrofones (casos offshore) os quais são regularmente distribuídos para cobrir áreas de centenas de quilômetros quadrados.

[0003] Avanços recentes nas técnicas de aquisição sísmica tem aperfeiçoado consideravelmente a iluminação da subsuperfície. Esta capabilidade de melhor compensar uma iluminação desigual representa uma vantagem para produção de novas imagens da subsuperfície. Aquisições sísmicas são tipicamente projetadas para serem fortemente redundantes, com intuito de aperfeiçoar a razão sinal para ruído, entretanto, mais dados potencialmente significam mais diversidade na qualidade sobre diferentes dados particionados. Uma ampla aquisição de azimutes pode resultar em áreas de dados localizados que são "pobres" em termos de qualidade para certos azimutes, e irão penalizar os dados "bons" resultantes daquelas localizações em outros cubos de azimute quando processados no empilhamento final.

[0004] Adicionalmente, a inadequação da física utilizada em diversos estágios do fluxo de trabalho do processamento sísmico não permite um benefício completo de todas as novas medições disponíveis. A imprecisão dos modelos de velocidade produz erros nos períodos de trajeto computados e estes levam a inconsistências no posicionamento dos eventos sísmicos em diferentes imagens pré-empilhamento; se não estão adequadamente alinhados antes da compactação, podem interferir, destrutivamente, resultando em uma imagem migrada desfocada. Avanços nas técnicas de estimativa de velocidade certamente levarão a modelos de subsuperfície mais adequados ao longo do tempo, porém uma certa quantidade de imprecisões no modelo de velocidade é algo que sempre terá que ser levado em consideração.

[0005] Os últimos anos deflagraram uma grande quantidade de atividades em diversos temas de aperfeiçoamento de imagens através de soma otimizada de pré- empilhamentos de imagens migradas. A maioria delas é oriunda da observação de que o empilhamento convencional é proficiente somente quando todos os traços possuem amplitudes similares e número de série e quando os padrões de ruído são estatisticamente independentes do ruído de qualquer outro traço e do sinal (Mayne, W. H., 1962, "Técnicas de empilhamento de dados horizontais de ponto de reflexão comum", 27, 927-938, doi: 10.1190/1.1439118; Neelamani, R., T. A. Dickens, e M. Deffenbaugh, 2006; "Empilhamento e não ruído: Um novo método para empilhar conjuntos de dados sísmicos: 76o Encontro Internacional Anual, SEG, Resumos Expandidos, 2827-2831; Robinson, J. C., 1970, "Empilhamento otimizado estatisticamente de dados sísmicos: Geofísica", 35, 436-446, doi: 10.1190/1.1440105). A ideia básica é aumentar a razão sinal para ruído através somente do empilhamento destes volumes que contêm informação relevante e consistente. Liu et al. (Liu, G., S. Fomel, L. Jin, e X. Chen, 2009, "Empilhamento de dados sísmicos utilizando correlação local: Geofísica", 74, no.: 3, V 43 a 48) e Compton et al (Compton, S., e C., Stork, "Aperfeiçoamento de empilhamento não linear 3D: Correlação com base em empilhamento", Resumos expandidos de programa técnico de SEG, 2012, pp 1-5) proposta que considera representar uma medida da correlação de local entre cada traço de entrada e um traço de referência de empilhamento convencional. Kun et al. (Kun, J., X. Cheng, D., Sun, e D., Vigh, 2014, "Aperfeiçoamento de imagens de migração através de alinhamento otimizado de partições de imagens de vetor", Resumos expandidos de programa técnico de SEG, 2014, pp. 3699-3703. doi: 10.1190/segam 2014-1648.1) aplica uma técnica de compatibilização modelo entre as imagens de pré empilhamento e uma imagem de referência espelhada para encontrar um alinhamento otimizado e os coeficientes de ponderação. Estes procedimentos são semi-automáticos e a qualidade do empilhamento final é, de alguma maneira, dependente da qualidade do traço de referência. O método de empilhamento de correlação local pode se beneficiar de técnicas mais avançadas para obter o traço de referência (Sanchis, C., e A. Hanssen, 2011, "Método de empilhamento de correlação local aperfeiçoada: Geofísica, 76, no.: 3, V 33 a 45).

[0006] A figura 1 ilustra diagramaticamente uma pesquisa de dados sísmicos com uma fonte S de ondas sísmicas e pelo menos um receptor G. Isto também mostra um ponto B da subsuperfície o qual é assumido para contribuir para o sinal medido por um dos receptores G. As coordenadas horizontais do ponto B da subsuperfície são denotadas por x, y (ou somente uma coordenada espacial se a imagem em 2D é considerada no lugar da imagem 3D), enquanto a profundidade é denotada por z. A figura 1 também fornece uma representação simplificada (linhas pontilhadas 101 e 102) da propagação de ondas sísmicas a partir da fonte S para o ponto B e do ponto B para o receptor G. As ondas são refratadas em descontinuidades das camadas geológicas onde a impedância acústica muda, e refletidas ou difratadas em diferentes posições incluindo aquela do ponto B.

[0007] A letra "a" representa a abertura e "m" o ponto médio comum.

[0008] Os dados registrados em uma inspeção sísmica incluem, para cada pulso de uma fonte S e para cada receptor G, um traço sísmico o qual é uma série de tempo do sinal medido pelo receptor G. Os traços para um número de pulsos devem ser transformados para fornecer uma imagem da subsuperfície a qual será o resultado do empilhamento ou integração de uma grande quantidade de informações. Uma etapa importante da transformação é a migração na qual consiste no rearranjo dos dados com relação a um modelo no qual o empilhamento pode ser realizado de forma coerente. O modelo é normalmente um mapa da velocidade de propagação das ondas acústicas na subsuperfície. Não é conhecido, a princípio, e se configura como um principal desafio de todas as tecnologias de imagens sísmicas determinar um modelo que considere adequadamente os dados de campo após o empilhamento.

[0009] Nos métodos de migração de profundidade de pré-empilhamento (PSDM), dados migrados são computados para cada pulso utilizando o modelo de velocidade e disposto em um cubo de saída contendo valores migrados associados com posições na subsuperfície. Os cubos obtidos para diferentes pulsos são, então, analisados para verificar a consistência do modelo: estes podem ser utilizados tanto para obter uma imagem final ou para produzir Coletas de Imagens Comuns. O modelo pode ser corrigido e o processo é iterado até que uma imagem satisfatória seja obtida.

[0010] Coletas de Imagens Comuns (CIGs) são ferramentas populares para avaliar o campo de velocidade de migração e para aperfeiçoamento de imagens. Estas são feitas de dados extraídos dos cubos de saída, combinados de uma maneira conveniente para análise de modo a verificar o modelo de velocidade. Uma CIG é uma estrutura de dados bidimensional definida por uma dada posição horizontal x, y, com um primeiro eixo representando a profundidade "z" e o segundo eixo representando um parâmetro de domínio "A" referido para combinar os dados dos cubos de saída. Este contêm valores de refletividade obtidos dos cubos de saída resultantes da migração, formando uma imagem que pode ser analisada para verificar e/ou corrigir o modelo de velocidade. Nesta imagem, um valor de pixel em um ponto (z, A) representa um valor migrado derivado como uma contribuição da posição de subsuperfície x, y, z em um traço sísmico associado com o parâmetro de domínio A. Exemplos de parâmetros A de domínio comumente utilizados incluem a abertura a, nominalmente a distância entre o centro C da localização de fonte e o localização do receptor G e a projeção ortogonal do ponto B na superfície.

[0011] A computação de imagens comuns não é direta em todos os métodos de extrapolação de campos de onda.

[0012] A abertura de migração é um parâmetro crítico em migração de Kirchhoff para obter a melhor qualidade de imagem de uma dada folha de dados. A seleção da abertura ótima é o resultado de um ajuste. Reduzir a abertura de migração geralmente aumenta a razão sinal para ruído, porém em detrimento da formação de imagens de eventos rebaixados. Infelizmente, métodos de migração de equação de onda necessitam desviar desta metodologia.

[0013] Entretanto, tais ferramentas foram utilizadas, em sua maioria, em métodos de migração com base na estimativa de tempos de trajetos entre refletores e a superfície. Métodos de migração mais sofisticados foram desenvolvidos para construir imagens de PSDM através da resolução da equação de onda, de modo a obter amplitudes de refletor e posicionamento estrutural mais precisos. Por exemplo, a migração de tempo reverso (RTM) é uma solução de migração de duas formas, a qual descreve precisamente a propagação de onda em meio complexo. Isto é utilizado, de forma crescente, em explorações sísmicas em virtude dos avanços em potência computacional e programação.

[0014] As ferramentas de análise mencionadas acima não são utilizadas com os métodos PSDM de equação de onda, incluindo RTM, pois não é conhecida uma maneira de computar as CIGs indexadas de abertura nestes métodos (precisamente, as CIGs indexadas de abertura podem ser computadas somente, até então, através do método de Kirchhoff).

[0015] É desejável, portanto, obter CIGs indexadas de abertura com diferentes tipos de métodos de migração, em particular métodos de equação de onda incluindo RTM ou métodos de migração de uma forma, de modo a manter as vantagens dos métodos de campo de onda e, ao mesmo tempo, focar na limitação suposição assintótica de métodos com base em raios, enquanto combina os cubos migrados da mesma maneira que as coletas de deslocamento de superfície clássicas.

SUMÁRIO DA INVENÇÃO

[0016] Um método de processamento de dados sísmicos é proposto para fornecer coletas de imagens comuns no domínio de abertura de superfície. O método compreende: receber dados sísmicos, para cada um de uma pluralidade de pulsos em suas respectivas localizações de fonte, traços sísmicos registrados em uma pluralidade de localizações de receptor; aplicar um processo de migração de profundidade aos dados sísmicos para obter um primeiro conjunto de dados migrados incluindo, para cada pulso, primeiros valores migrados respectivamente associados com uma pluralidade de posições de subsuperfície; computar dados modulados de ponto médio ao multiplicar os dados sísmicos em cada traço sísmico através de uma função do centro da fonte e localizações de receptor para o dito traço sísmico; aplicar o processo de migração de profundidade aos dados de modulados de ponto médio para obter um segundo conjunto de dados migrados incluindo, para cada pulso, segundos valores migrados respectivamente associados com a pluralidade de posições de subsuperfície; para cada pulso, estimar valores de abertura respectivamente associados com pelo menos alguma das posições de subsuperfície, através de um processo de divisão aplicado ao primeiro e ao segundo conjuntos de dados migrados; e estimar uma coleta de imagens comuns indexadas de abertura em uma posição horizontal, compreendendo os respectivos valores migrados R(x,y,z,a) para pares de parâmetros, cada qual incluindo um parâmetro de profundidade e um parâmetro de abertura;

[0017] O processo de divisão pode ser realizado em um domínio de Radon.

[0018] O valor migrado para um par de parâmetros na coleta de imagens comuns na dita posição horizontal é um primeiro valor migrado do primeiro conjunto de dados migrados associado com uma posição de subsuperfície determinada pela dita posição horizontal e o parâmetro de profundidade do dito par de parâmetros para um pulso tal que o valor de abertura estimado associado com a dita posição de subsuperfície é o parâmetro de abertura do dito par de parâmetros.

[0019] As coletas indexadas de abertura de superfície podem ser obtidas sem se referir ao paradigma de Kirchhoff tradicional, o qual tem como base uma aproximação de frequência infinita da equação de onda e traçamento de raios. O método é, portanto, aplicável, a diversos tipos de processos de migração, incluindo processos de migração de profundidade pré-empilhamento de equação de onda (PSDM).

[0020] Se um método para determinar o deslocamento das CIGs pode ser conhecido, este método permite uma determinação direta de uma abertura indexada às CIGs sem computação adicional.

[0021] Em uma modalidade preferida, o processo de migração de profundidade é um processo de migração de tempo reverso (RTM).

[0022] O RTM fornece um melhor propagador do que o traçamento de raios. Portanto, este permite produzir CIGs em áreas complexas como de pré-sal, sobrecarregamentos complexos, etc. Além disso, uma vez que a combinação do resultado é a mesma, as coletas indexadas de abertura de Kirchhoff convencionais podem ser substituídas pelas CIGs obtidas através do método como entrada na maioria das ferramentas de pós-processamento padrão e, portanto, podem ser utilizadas para um empilhamento otimizado. Devido aos anos de know-how na migração do tipo Kirchhoff, este método permite uma interpretação mais fácil dos resultados migrados.

[0023] Para obter estabilidade aperfeiçoada, o processo de divisão utilizado para estimar os valores de abertura associados com uma posição de subsuperfície pode compreender uma função de custo definida por uma variável de abertura e valores locais do primeiro e do segundo valores migrados em uma vizinhança da posição de subsuperfície.

[0024] Outro aspecto da invenção se refere a um sistema para processamento de dados sísmicos, compreendendo recursos de computação configurados para executar um método de análise de dados sísmicos conforme definido acima.

[0025] Ainda outro aspecto da invenção se refere a um produto de programa de computador para um sistema de processamento de dados sísmicos, compreendendo instruções para executar um método conforme definido acima quando o dito produto de programa está sendo executado em uma unidade de processamento de computador do sistema.

[0026] Outras características e vantagens do método e aparelho revelados neste documento se tornarão aparentes a partir da descrição a seguir de modalidades não limitativas, com referência às figuras em anexo.

BREVE DESCRIÇÃO DAS FIGURAS

[0027] A figura 1 é um diagrama esquemático ilustrando a aquisição de dados sísmicos.

[0028] A figura 2 é um fluxograma de um método de processamento de dados sísmicos de acordo com uma modalidade da invenção.

[0029] A figura 3 é um diagrama ilustrando a derivação das CIGs indexadas de abertura de acordo com o método.

DESCRIÇÃO DAS MODALIDADES

[0030] Uma forma de obter CIGs indexadas de abertura para PSDM (ex.: métodos de migração de RTM) seria computar uma migração por pulso e por traço, exigindo um número de migrações igual ao número total de vezes de pulsos o número médio de receptores por pulso. Isto é claramente impraticável para o tempo, especialmente para casos de 3D. Ao invés disso, é proposto utilizar uma solução mais viável, conhecida como migração de atributo, também chamada de migração dupla.

[0031] No método da migração dupla conforme introduzido por N. Bleistein ("Na produção de imagens de refletores na terra", Geofísica, Vol. 52, No.: 7, Julho de 1987, pp. 931-942), duas migrações são computadas com os mesmos dados, a segunda envolvendo um operador de migração multiplicado pelo ângulo de reflexão especular. A divisão de duas imagens de migração, portanto, proporciona o ângulo especular ao longo dos refletores.

[0032] Um método similar pode ser utilizado com o ponto médio comum (conforme definido acima) no lugar do ângulo de reflexão especular conforme o atributo migrado. A migração pode ser realizada utilizando diversos métodos de PSDM incluindo RTM de gravação de pulso padrão.

[0033] A metodologia, portanto, compreende (ex: para casos 2D, 3D existem três imagens onde duas imagens são codificadas com CMPx e CMPy respectivamente): realizar uma primeira migração de RTM padrão (etapa 20 na figura 3); realizar uma segunda migração de RTM de atributo de ponto médio, onde os dados são multiplicados pelos valores (m) de pontos médios comuns antes da migração (etapas 30 e 40 na figura 3); realizar uma divisão dos dois dados migrados obtidos, em pelo menos um sentido quadrado, para obter o ponto médio comum (etapa 50 na figura 3); para cada ponto da imagem, subtrair a coordenada lateral da localização correspondente para obter o atributo migrado de abertura; adicionar a refletividade obtida do primeiro dado migrado para a subsuperfície correspondente e localização de coleta indexada de abertura dada pelo mapa de atributo obtido para reconstruir a CIGs de abertura migrada de atributo (etapa 60 na figura 3).

[0034] Nesta maneira, é possível obter CIGs indexados de abertura clássica utilizando uma migração de PSDM tal como RTM, a qual é o melhor método de campo de onda disponível atualmente para migração sísmica.

[0035] Na figura 2, os dados de entrada de traços sísmicos na etapa 10 as medições de campo são notadas como DS,G[t], onde S denota uma localização de fonte, G denota uma localização de receptor e t denota tempo. Cada traço é modulado na etapa 30 ao multiplicá-la pelo valor de ponto médio correspondente, nominalmente a distância de ponto médio horizontal entre o centro do segmento definido pela localização de fonte S e a localização de receptor G e a localização de receptor G ou a localização de fonte S. Os traços modulados de ponto médio são D'SG[t]=m.DSG[t].

[0036] Os dados sísmicos DSG[t] e os dados modulados de ponto médio D'SG[t] são migrados, respectivamente, nas etapas 20 e 40 para fornecer imagens/dados de PSDM MS[x,y,z] e M'S[x,y,z]. O primeiro conjunto de dados migrados produzidos na etapa 20 inclui, para cada pulso em uma localização de fonte S, um cubo de valores migrados MS[x,y,z] associado com as posições de subsuperfície x, y, z. Da mesma forma, o segundo conjunto de dados migrados obtidos na etapa 40 utilizando o mesmo processo de migração de profundidade inclui outro cubo de valores migrados M'S[x,y,z] para cada pulso.

[0037] Antes de acessar o valor de abertura, é mandatório aplicar um processo de subtração compreendendo a projeção da dita distância lateral entre a dita posição de imagem e a localização de ponto médio comum.

[0038] Com intuito de estimar o valor de abertura âs[x,y,x] para um valor migrado MS[x,y,z], ou seja, um valor para um pulso S e uma posição de subsuperfície x, y, z, um processo de divisão é realizado na etapa 50 para avaliar M'S[x,y,z]/MS[x,y,z].

[0039] Uma divisão bruta dos dois números pode gerar dúvidas sobre sua estabilidade, porém é possível. Ao invés disso, é melhor elencar a divisão como um conjunto dos locais de problemas de quadrado mínimo. O valor de abertura âS[x,y,z] é, então, encontrado ao minimizar uma função de custo JS,X,Y,Z(a) definido em uma vizinhança

centralizada na localização x,y,z. Uma expressão possível da função de custo JS,X,Y,Z(a) é:

onde MaS[x,y,z] = MS[x,y,z] + i.H(MS[x,y,z]) é o sinal analítico da refletividade, H denotando a transformada de Hilbert, e M'aS [x, y, z] = M'S[x,y,z] + i.H(M'S[x,y,z]). O tamanho da vizinhança (x,y,z) é variável e pode depender da aplicação. É selecionado de tal modo que o valor do atributo migrado a pode razoavelmente ser assumido como constante ao longo de (x,y,z) para um dado pulso.

[0040] Os valores de abertura âs[x,y,z] pode, portanto, ser contido como:

onde o sobrescrito (.)* denota o Hermitian. Note que a parte superior do lado direito da equação (2) é a correlação cruzada das duas imagens, a imagem migrada padrão e a imagem migrada de atributo. A parte inferior é a auto-correlação da imagem migrada padrão, ou o envelope. Desta forma, a estabilidade da divisão é aumentada.

[0041] Os valores de abertura âs[x,y,z] logicamente obtidos são utilizados para mapear os valores de refletividade MS[x,y,z] para aberturas correspondentes a, em que os valores de refletividade podem, portanto, ser arranjados como CIGs indexados de abertura de superfície nas posições horizontais x, y:

[0042] Este processo de computação das CIGs indexadas de abertura de superfície é ilustrado no diagrama da figura 3. O atributo de abertura é obtido através da divisão de mapa "migrado duplo". Os dois painéis superiores do diagrama descrevem duas imagens migradas de pulso comum mostrando um refletor horizontal. A primeira imagem é obtida por meio de Migração de Tempo Reverso de dados sísmicos puros, enquanto a segunda é o resultado de migração onde os dados têm sido multiplicados pela distância de ponto médio comum "m". Neste painel, amplitudes variam lateralmente ao longo do refletor e são proporcionais à posição de receptor e, portanto, à abertura. Em particular, iniciando pelo lado esquerdo, a amplitude é negativa, esta alcança zero exatamente abaixo do ponto C conforme definido na figura 1 e se torna positiva, posteriormente.

[0043] Para uma posição específica (x0,y0) na imagem migrada de pulso (nesta modalidade, a posição horizontal x0 pode ser 2D, com componentes x e y), existe um valor particular de refletividade R. Na mesma posição na imagem migrada de atributo, o valor da refletividade é R multiplicado pelo ponto médio comum m. A abertura é função da divisão das duas quantidades. Com estes quatro valores (x0,z0,a,R), pode ser construído o cubo de abertura migrada. O painel inferior do diagrama da figura 3 representa uma seção de ponto médio comum na localização x0. A contribuição da refletividade R é adicionada à posição cujas coordenadas são dadas pela dupla profundidade/abertura (z0,a).

[0044] O método proposto para obter coletas de superfície clássica para várias técnicas de PSDM incluindo Migração de Tempo Reverso mostrada como sendo aplicável de forma bem sucedida no contexto da propagação de banda limitada. Os métodos propostos não dependem de uma implementação particular do método de extrapolação de campo de onda, uma vez que podem ser realizados após migração de registro curto no domínio de Fourier, domínio de tempo, etc.

[0045] Uma vantagem deste método é que este permite o uso de propagadores melhores do que raios para propagar o campo de ondas, tornando possível adotar o arsenal completo de ferramentas padrão para processamento posterior desenvolvido de esquemas de migração assintótica. Comparados ao Deslocamento Subsuperficial e Deslocamento Subsuperficial de Ângulo de Dispersão, CIGs de Abertura para RTM são muito menos onerosos para computar, e estes preservar o movimento cinético.

[0046] As modalidades do método descrito neste documento podem ser implementadas em qualquer forma de computador ou computadores e os componentes podem ser implementados como aplicações dedicadas ou em arquiteturas de servidor de cliente, incluindo uma arquitetura com base em rede, e pode incluir programas funcionais, códigos e segmentos de código. Qualquer um dos computadores podem compreender um processador, uma memória para armazenar dados de programa e executá-los, um armazenamento permanente tal como um acionamento de disco, um pórtico de comunicações para manuseio de comunicações com dispositivos externos, e dispositivos de interface de usuário, incluindo um monitor, um teclado, mouse, etc.

[0047] Tipicamente, o método é realizado utilizando módulos de software os quais são armazenados como instruções de programa ou códigos legíveis por computador executáveis no processador em uma mídia legível por computador tal como uma memória somente de leitura (ROM), memória de acesso aleatório (RAM), CD-ROMs, fitas magnéticas, discos floppy, e dispositivos de armazenamento de dados ópticos. A mídia de armazenamento legível por computador pode também ser distribuída ao longo da rede acoplada aos sistemas de computadores de tal modo que o código legível por computador seja armazenado e executado de uma maneira distribuída. Esta mídia é legível pelo computador, armazenada na memória, e executada pelo processador.

[0048] Será verificado que as modalidades descritas acima são ilustrativas da invenção revelada neste documento e diversas modificações podem ser feitas sem se afastar do escopo tal como definido nas reivindicações em anexo.

Claims

1. Método de processamento de dados sísmicos, caracterizado por compreender: - receber dados sísmicos (DS,G[t]) incluindo, para cada um de uma pluralidade de pulsos em suas respectivas localizações de fonte, traços sísmicos registrados em uma pluralidade de localizações de receptor; - aplicar um processo de migração de profundidade aos dados sísmicos para obter um primeiro conjunto de dados migrados incluindo, para cada pulso, primeiros valores migrados (MS[x,y,z]) respectivamente associados com uma pluralidade de posições de subsuperfície; - computar dados modulados de ponto médio (D'S,G[t]) multiplicando os dados sísmicos em cada traço sísmico através de uma função de um centro da fonte e localizações de receptor para o dito traço sísmico; - aplicar o processo de migração de profundidade aos dados modulados de ponto médio para obter um segundo conjunto de dados migrados incluindo, para cada pulso, segundos valores migrados (M'S[x,y,z]) respectivamente associados com a pluralidade de posições de subsuperfície; - aplicar um processo de subtração compreendendo a projeção de uma distância lateral entre uma posição de imagem e a localização de ponto médio comum; - para cada pulso, estimar valores de abertura (âS[x,y,z]) respectivamente associados com pelo menos alguma das posições de subsuperfície, através de um processo de divisão realizado em um domínio Radon aplicado ao primeiro e ao segundo conjunto de dados migrados; e - estimar uma coleta de imagens comuns indexadas de abertura em uma posição horizontal, compreendendo respectivos valores migrados (Rx,y,[z,a]) para pares de parâmetros, cada qual incluindo um parâmetro de profundidade e um parâmetro de abertura; em que os valores migrados para um par de parâmetros na coleta de imagens comuns na dita posição horizontal é um primeiro valor migrado do primeiro conjunto de dados migrados associados com uma posição de subsuperfície determinada pela dita posição horizontal e o parâmetro de profundidade do dito par de parâmetros para um pulso tal que o valor de abertura estimada com a dita posição de subsuperfície é o parâmetro de abertura do dito par de parâmetros.

2. Método, de acordo com a reivindicação 1, caracterizado pelo fato de o processo de divisão utilizado para estimar os valores de abertura (âS[x,y,z]) associados com uma posição de subsuperfície compreender minimizar uma função de custo definida por uma variável de abertura e valores locais dos primeiro e segundo valores migrados (MS[x,y,z], M'S[x,y,z]) em uma vizinhança da dita posição de subsuperfície.

3. Método, de acordo com qualquer uma das reivindicações 1 ou 2, caracterizado pelo fato de o processo de migração de profundidade ser um processo de migração de tempo reverso (RTM).

4. Método, de acordo com qualquer uma das reivindicações 1 ou 2, caracterizado pelo fato de o processo de migração de profundidade ser um processo de migração de profundidade de pré-empilhamento de equação de onda.

5. Sistema para processamento de dados sísmicos, caracterizado pelo fato de compreender recursos de computação configurados para executar um método conforme definido em qualquer uma das reivindicações 1 a 4.

6. Mídia de gravação legível por computador caracterizada pelo fato de compreender instruções para executar um método conforme definido em qualquer uma das reivindicações 1 a 4, quando as ditas instruções são executadas em uma unidade de processamento de computador do sistema para processamento de dados sísmicos.