BR112013009050B1

BR112013009050B1 - sistema e método para modelação e migração de dados sísmicos

Info

Publication number: BR112013009050B1
Application number: BR112013009050-2A
Authority: BR
Inventors: Guojian Shan; Linbin Zhang; Yue Wang
Original assignee: Chevron U.S.A. Inc.
Priority date: 2011-03-23
Filing date: 2012-02-21
Publication date: 2020-12-08
Also published as: EA201391376A1; US20120243371A1; CN103229075A; CA2817622A1; CN103229075B; WO2012128879A2; CA2817622C; AU2012231651A1; EP2689273B1; EP2689273A2; EP2689273A4; US8614930B2; AU2012231651B2; BR112013009050A2; WO2012128879A3

Abstract

SISTEMA E MÉTODO PARA MODELAÇÃO E MIGRAÇÃO DE DADOS SÍSMICOS. Um sistema e método implementado por computador para estimar a dispersão temporal em propagação de ondas sísmicos de diferenças finitas de baixa ordem são divulgados. Uma modalidade do método inclui a transformação de um conjunto de dados sísmicos do domínio de tempo para domínio de frequência para obter um conjunto de dados sísmicos no domínio de frequência, se aplicando um filtro de variação de tempo de domínio de frequência ao conjunto de dados sísmicos no domínio de frequência para obter um conjunto de dados sísmicos no domínio da frequência filtrada, e transformando o conjunto de dados sísmicos de frequência de domínio filtrado do domínio de freqüência para o domínio de tempo para obter um conjunto de dados sísmicos filtrado do domínio de tempo. O filtro de variação de tempo de domínio de frequência é baseado na velocidade de fase efetiva inerente a uma solução de diferenças finitas para a equação de onda.

Description

CAMPO TÉCNICO

[0001] A presente invenção se refere genericamente a métodos e sistemas para a modelação e migração de dados sísmicos, utilizando operadores de modelação de diferenças finitas, e, em particular métodos e sistemas para estimar a dispersão temporal causada pela utilização de operadores de modelação de diferenças finitas de baixa ordem para modelação e migração de dados sísmicos.

ANTECEDENTES DA INVENÇÃO

[0002] A migração exata de dados sísmicos 3D permite uma interpretação adequada de reservatórios de hidrocarbonetos de subsuperfície. A migração sísmica é essencialmente o processo de reversão de propagação de ondas sísmicas; assim sendo, o enfoque tem sido posto na modelação de propagação de ondas sísmicas a mais exata possível. Um método de modelação de propagação da onda sísmica é a modelação de diferenças finitas.

[0003] Na modelação de diferenças finitas, a solução para a equação de onda é aproximadamente a utilização de um método de diferenças finitas (FD). Este método pode produzir múltiplas soluções aproximadas que são referidas de acordo com a sua ordem, por exemplo, FD de segunda ordem e FD de quarta ordem, o que indica a precisão com que representam a verdadeira solução para a equação de onda. Quando uma solução FD de segunda ordem é utilizada para modelar a propagação de ondas sísmicas, os dados sísmicos resultantes têm uma dispersão temporal significativa devido à baixa precisão da solução FD. Soluções FD de ordem superior tais como as de quarta ordem para a produção de dados sísmicos sintéticos com menos dispersão temporal. A precisão melhorada de soluções de elevada ordem tem um grande custo computacional. Por exemplo, modelação sísmica FD de quarta ordem requer o dobro de operações de computação em relação à FD de segunda ordem.

[0004] Outro método que pode ser utilizado para melhorar a precisão de modelação é um operador pseudo-analítico, tal como um pseudo-Laplaciano. Este método é semelhante à modelação FD de segunda ordem excetuando que este modifica os derivados espaciais e temporais, de modo a que eles tenham sinais opostos, e, ajustando os coeficientes, os erros nos derivados irão se contrapor mutuamente, reduzindo assim a imprecisão no resultado. Este método é mais preciso do que a modelação FD de segunda ordem e não tão computacionalmente caro como a modelação FD de quarta ordem. No entanto, métodos pseudo- analíticos são ainda mais caros do que a modelação FD de segunda ordem.

[0005] O custo computacional associado à modelação FD de ordem elevada ou métodos pseudo-analíticos se torna ainda mais significativo quando consideradas em termos de migração reversa no tempo (RTM). Em migração reversa no tempo, um campo de onda de fonte é propagada para a frente na subsuperfície, muitas vezes usando modelação pseudo-analítica ou FD, enquanto um conjunto de dados sísmicos gravados é propagado para trás na subsuperfície. Os dois campos de ondas são combinados em locais de subsuperfície por meio de uma condição de imagem, freqüentemente atraso zero na correlação cruzada, para criar uma imagem. O conjunto de dados sísmicos gravados é o resultado real da energia sísmica que tenha passado através da subsuperfície e, como tal, não tem a dispersão temporal que surge a partir da modelação pseudo-analítica ou FD aproximada. Para que os dados sísmicos sintéticos se propaguem para a frente para coincidir com os dados sísmicos gravados que se propagam para trás, a dispersão temporal deve ser considerada. Portanto, os métodos atuais de RTM usam modelação pseudo- analítica ou FD de ordem superior.

[0006] A prática atual para modelação sísmica FD e RTM usa métodos pseudo- Laplaciano ou métodos FD de ordem superior. Estes métodos são mais precisos e computacionalmente mais caros do que a modelação FD de segunda ordem convencional.

SUMÁRIO DA INVENÇÃO

[0007] De acordo com uma implementação da presente invenção, um método para estimar a dispersão temporal inclui a transformação de um conjunto de dados sísmicos do domínio de tempo para o domínio de freqüência para obter um conjunto de dados sísmicos no domínio de freqüência, se aplicando um filtro de variação de tempo de domínio de freqüência para o conjunto de dados sísmicos no domínio de freqüência para obter um conjunto de dados sísmicos no domínio da freqüência filtrada, e transformando o conjunto de dados sísmicos de freqüência de domínio filtrado do domínio de freqüência para o domínio de tempo para obter um conjunto de dados sísmicos filtrados do domínio de tempo. O filtro de variação de tempo de domínio de freqüência é baseado na velocidade de fase efetiva inerente a uma solução de diferenças finitas para a equação de onda. o filtro de variação de tempo no domínio da freqüência pode ser aplicado a um conjunto de dados sísmicos sintéticos que foram gerados por modelação de diferenças finitas de baixa ordem.

[0008] Noutra modalidade, o filtro de variação de tempo do domínio de freqüência pode ser usado para gravar dados sísmicos anteriores à migração reversa no tempo. A migração reversa no tempo pode então usar modelação de diferenças finitas de baixa ordem para propagar para a frente o campo de ondas de fonte.

[0009] A presente invenção também pode ser praticada como um sistema para estimar a dispersão temporal em métodos de processamento sísmico usando modelação de diferenças finitas de baixa ordem. O sistema pode incluir um dispositivo para fornecer informação representativa da área de subsuperfície de interesse e de um processador de computador, em comunicação com o dispositivo e configurado de modo a receber os dados de para a execução de um código executável por computador de resposta aos dados. O código executável por computador pode incluir um módulo de transformação de domínio capaz de transformar a partir de um domínio do tempo para o domínio da freqüência e/ou a partir do domínio da freqüência para o domínio do tempo, e um módulo de aplicação de filtro de domínio de freqüência para aplicação de um filtro de variação de tempo de domínio de freqüência. O sistema pode também incluir uma interface de usuário. Numa modalidade, o sistema pode incluir um módulo de migração reversa no tempo.

[0010] A presente invenção também pode ser praticada como um artigo de fabrico incluindo um meio legível por computador tendo um código legível por computador, sendo o código legível por computador configurado para implementar um método para estimar a dispersão temporal em propagação de onda de diferenças finitas de baixa ordem. O método pode incluir a transformação de um conjunto de dados sísmicos de domínio de tempo para o domínio de freqüência para se obter um conjunto de dados sísmicos no domínio da freqüência, se aplicando um filtro de variação de tempo de domínio de freqüência para o conjunto de dados sísmicos no domínio da freqüência para obter um conjunto de dados sísmicos no domínio da freqüência filtrada, e transformando o conjunto de dados sísmicos no domínio da freqüência filtrada do domínio da freqüência para o domínio do tempo para obter um conjunto de dados sísmicos filtrados do domínio do tempo. O filtro de variação de tempo de domínio de freqüência é baseado na velocidade de fase eficaz inerente a uma solução de diferenças finitas para a equação de onda.

[0011] A secção de sumário acima é fornecida para apresentar uma seleção de conceitos de uma forma simplificada que estão descritos abaixo na secção de descrição detalhada. O sumário não pretende identificar as características essenciais ou chave da matéria reivindicada, nem se destina a ser utilizado para limitar o escopo de matéria reivindicada. Além disso, a matéria reivindicada não está limitada a aplicações que resolvem qualquer ou todas as desvantagens observadas, em qualquer parte da presente divulgação.

BREVE DESCRIÇÃO DOS DESENHOS

[0012] Estas e outras características da presente invenção serão melhor compreendidas no que se refere à descrição seguinte, reivindicações pendentes e desenhos anexos, em que:

[0013] A Figura 1 é um fluxograma que ilustra um método para a realização de modelação de dados sísmicos, de acordo com uma modalidade da invenção;

[0014] A Figura 2 mostra uma pequena onda que foi propagada utilizando um método FD de segunda ordem e o resultado da aplicação do filtro de variação do tempo no domínio da freqüência da presente invenção à pequena onda;

[0015] A Figura 3 é um exemplo de modelação de dados sísmicos, utilizando um método FD de segunda ordem convencional comparado com a modelação de dados sísmicos utilizando a presente invenção;

[0016] A Figura 4 é um fluxograma que ilustra um método para a realização de migração reversa no tempo (RTM) de acordo com uma modalidade da invenção;

[0017] A Figura 5 é um exemplo de migração reversa no tempo usando um método FD de segunda ordem convencional em relação ao método de migração reversa no tempo utilizando a presente invenção; e

[0018] A Figura 6 ilustra esquematicamente um sistema para a realização de um método de acordo com uma modalidade da invenção, DESCRIÇÃO DETALHADA DO INVENTO

[0019] A presente invenção pode ser descrita e implementada no contexto geral de um sistema e métodos de computador para serem executados por um computador. Tais instruções executáveis por computador podem incluir programas, rotinas, objetos, componentes, estruturas de dados e tecnologias de software de computador que podem ser usadas para executar tarefas específicas e para processar tipos de dados de processos abstratos. Implementações de software da invenção atual podem ser codificadas em diferentes línguas para a aplicação numa variedade de plataformas e ambientes de computação. Se faz observar que o escopo e princípios subjacentes da presente invenção não se encontram limitados a qualquer tecnologia especial de software de computador.

[0020] Além disso, os peritos na técnica apreciarão que a invenção atual pode ser posta em prática utilizando qualquer uma ou uma combinação de configurações de hardware e de software, incluindo mas não se limitando a um sistema que tem um sistema de processadores de computador de processadores simples e/ou múltiplos, dispositivos portáteis, electrônica de consumidor programável, minicomputadores, computadores “mainframe”, e similares. A invenção também pode ser posta em prática em ambientes de computação distribuída em que as tarefas são realizadas por servidores ou outros dispositivos de processamento que estão ligados através de uma ou mais redes de comunicações de dados. Num ambiente de computação distribuída, os módulos do programa podem estar localizados em meios de armazenamento de computadores locais e remotos, incluindo dispositivos de armazenamento de memória.

[0021] Além disso, um artigo de fabrico para uso com um processador de computador, tal como um CD, um disco pré-gravado ou outros dispositivos equivalentes, pode incluir um meio de armazenamento de computador e meios de programa gravados no mesmo para dirigir o processador de computador de modo a facilitar a aplicação e a prática da presente tecnologia. Tais dispositivos e artigos de fabrico também caem dentro do espírito e do escopo da presente invenção.

[0022] Referindo-nos agora aos desenhos, serão descritas modalidades da presente invenção. A invenção pode ser implementada de várias maneiras, incluindo, por exemplo, um sistema (incluindo um sistema de processamento de computador), um método (incluindo um método implementado por computador), um aparelho, um meio legível por computador, um produto de programa de computador, uma interface gráfica de usuário, um portal Web, ou uma estrutura de dados fixa de forma tangível numa memória legível por computador. Várias modalidades da presente invenção são discutidas abaixo. Os desenhos em anexo ilustram apenas modalidades típicas da presente invenção e, portanto, não devem ser consideradas como limitando o seu escopo e amplitude.

[0023] A presente invenção se refere à consideração da dispersão temporal na modelação de dados sísmicos de diferenças finitas, e a título de exemplo e não como limitação, pode ser utilizada para reduzir a dispersão temporal em dados sísmicos sintéticos gerados por modelação de diferenças finitas de baixa ordem ou introduzir dispersão temporal em dados sísmicos gravados antes de reverter a migração do tempo que usa modelação de diferenças finitas de baixa ordem como propagador para a frente. A dispersão temporal é caracterizada por um filtro de variação de tempo de domínio de freqüência.

[0024] O inventor determinou que é possível gerar um filtro de variação de tempo de domínio de freqüência baseado na inerente velocidade de fase eficaz numa solução de diferenças finitas para a equação de onda. As soluções de diferenças finitas de baixa ordem, tais como a solução de segunda ordem, serão imprecisas e irão introduzir dispersão temporal nos dados sísmicos sintéticos gerados por modelação de diferenças finitas (FD) de baixa ordem. Além disso, o inventor determinou que, ao aplicar o filtro de variação de tempo de domínio de freqüência para dados sísmicos sintéticos criados por modelação FD de baixa ordem, a dispersão temporal pode ser reduzida. Além disso, é possível utilizar o filtro de variação do tempo de domínio de freqüência para adicionar dispersão temporal para gravar dados sísmicos de modo que ele pode ser usado como entrada para um método de migração reversa no tempo utilizando modelação FD de baixa ordem como seu propagador para a frente.

[0025] A este respeito, um exemplo de um método 100 de acordo com a presente invenção é ilustrado no fluxograma da Figura 1. No passo 10, um conjunto de dados sísmicos sintético é transformado do domínio do tempo para o domínio da freqüência. Isto pode ser feito, por exemplo, por uma transformada rápida de Fourier. Numa modalidade, o conjunto de dados sísmicos sintético foi gerado por um operador de modelação de diferenças finitas (FD) de baixa ordem, tais como, e não se limitam a um operador de modelação de segunda ordem FD. Um operador de modelação FD de baixa ordem é aquele que se baseia numa solução de diferenças finitas para a equação de onda e é conhecido por não ser exacto, mas é utilizado porque é menos dispendioso computacionalmente do que um operador de modelação FD de alta ordem.

[0026] No passo 12, o conjunto de dados de sísmica sintética que está agora no domínio da freqüência tem um filtro de variação de tempo de domínio de freqüência variável aplicado a ele. Este filtro de variação de tempo de domínio de freqüência é projetado com base na velocidade de fase eficaz que é inerente a uma solução FD de baixa ordem para a equação de onda. Como um exemplo de como um filtro pode ser construído, para a equação de onda acústica 3D

[0027] Equação 1

[0028] em que P é o campo de ondas, t é o tempo de percurso e V é a velocidade, a solução de diferenças finitas de segunda ordem pode ser escrita como:

[0029] Equação 2

[0030] Considerando a transformada de Fourier da Equa. 2 Resulta em

[0031] Equação 3

[0032] onde K é o número de onda. A partir deste resultado, é possível determinar a velocidade de fase eficaz para a solução FD de segunda ordem. Neste caso, a velocidade de fase efetiva Va da Equa. 3 pode ser expressa como:

[0033] Equação 4

[0034] A partir da Equa. 4 é evidente que Va é freqüentemente dependente e diferente da velocidade real V. A diferença entre V e Va produz um atraso de tempo, ou dispersão temporal, num campo de onda que é propagado de acordo com esta solução FD. A dispersão temporal pode ser descrita matematicamente:

[0035] Equação 5

[0036] em que T representa o tempo de percurso calculado pelo modelo de diferenças finitas e True de percurso é o tempo de percurso que usa a velocidade real. Uma vez que a dispersão temporal é definida, um filtro de variação de tempo de domínio de freqüência F(w,T) pode ser concebido para compensar isso:

[0037] Equação 6

[0038] Aplicando este filtro para dados sísmicos sintéticos gerados por modelação FD de segunda ordem FD, como é feito no passo 12, reduz a dispersão temporal.

[0039] O efeito do filtro na Equa. 6 pode ser visto na Figura 2. Neste exemplo, uma única pequena onda é mostrada como curva 20 e foi propagada por modelação FD de segunda ordem. O resultado mostra a distorção de fase antecipada como um pequeno pico antes da chegada principal e os picos grandes negativos e desiguais positivos que são característicos de dispersão temporal causada pela solução FD imprecisa. Aplicando o filtro da Equa. 6 para a curva 20 que produz a curva 22, a qual não apresenta sinais de dispersão temporal.

[0040] Referindo-nos novamente à Figura 1, o filtro da Equa. 6 é aplicado aos dados sísmicos de domínio de freqüência, os dados são transformados de volta para o domínio do tempo, no passo 14. Um exemplo do método 100 pode ser visto na Figura 3. O painel 30 mostra um campo de onda que foi propagado por um operador de modelação FD de segunda ordem. A dispersão temporal é indicada pela seta 31. Painel 32 mostra dados após realização do método 100. A dispersão temporal é muito reduzida.

[0041] Uma outra modalidade da presente invenção está ilustrada no fluxograma da Figura 4. Nesta modalidade, um filtro de variação de tempo no domínio da freqüência é utilizada para preparar dados sísmicos gravados para a migração através de um algoritmo de migração inversa no tempo (RTM). O conjunto de dados sísmicos gravados representa energia sísmica que foi gerada por uma fonte sísmica, propagada através da subsuperfície, refração e reflexão dependendo das diferenças na velocidade sísmica e densidade das camadas de rocha de subsuperfície e foi gravada por receptores sísmicos. Como tal, os dados sísmicos gravados representam a verdadeira solução para a equação de onda e por isso não tem a dispersão temporal, que é causada por imprecisões na modelação FD. Na RTM, o conjunto de dados sísmicos gravados ou o campo de ondas é propagado para trás na subsuperfície, enquanto um campo de ondas de fonte é propagado para a frente usando um método FD ou pseudo-analítico. Os dois campos de ondas são submetidos a uma condição de imagem para criar a imagem migrada. Uma condição típica de imagem é O atraso zero na correlação cruzada:

[0042] Equação 7

[0043] onde I(x,y,z) é a imagem migrada, S(x,y,z,t) é o campo de ondas de fonte que está propagada para a frente por um método FD ou pseudo-analítico, e R(x,y,z,t) é o campo de ondas gravado que é propagado para trás. Devido ao atraso zero na correlação cruzada, qualquer dispersão temporal, que ocorre em apenas num dos campo de ondas de fonte propagadas para a frente irá provocar artefatos na imagem migrada e profundidade mal posicionada. Normalmente RTM usa a modelação pseudo-analítica ou FD de ordem elevada computacionalmente mais cara para reduzir a dispersão temporal. No entanto, se a dispersão temporal na fonte e campo de ondas gravados são o mesmo, a condição de imagem de atraso zero na correlação cruzada não irá causar artefatos. Esta modalidade da presente invenção, cria dispersão temporal nos dados sísmicos gravados, de modo que possam ser utilizados como entrada para o RTM usando computacionalmente um operador de modelação FD de baixa ordem mais barato para o campo de ondas de fonte, enquanto reduz artefatos.

[0044] No método 400, o passo 40 transforma um conjunto de dados sísmicos registados no domínio da freqüência. Isto pode ser feito, por exemplo, por uma transformada rápida de Fourier. Os dados de freqüência de domínio possuem um filtro de variação de tempo no domínio da freqüência no passo 42. Este filtro é gerado com base na velocidade de fase eficaz a partir da solução FD, como era o derivado de filtro para o método 100, mas, nesta modalidade, o filtro irá introduzir dispersão temporal, ao invés vez de reduzir. Para uma solução de FD de segunda ordem, o filtro de variação de tempo de domínio de freqüência que irá ser aplicado aos dados sísmicos registados de domínio de freqüência é

[0045] Equação 8

[0046] Se note que este filtro é quase idêntico ao filtro do método 100 mostrado na Equa. 6; apenas o sinal no expoente é diferente.

[0047] No passo 44, o conjunto de dados gravados filtrados é transformado de volta para o domínio do tempo. No passo 46, o conjunto de dados registados filtrado é usado como entrada para RTM que usa um operador de modelação FD da mesma ordem que a solução FD usada para derivar o filtro de variação de tempo de domínio de freqüência. Deve se notar que, se o método RTM está preparado para ter dados sísmicos gravados do domínio da freqüência como entrada, o passo 44 pode ser ignorado.

[0048] Um exemplo do método 400 pode ser visto na Figura 5. O painel 50 apresenta um resultado RTM que usou um segundo operador de modelação FD de segunda ordem para propagar campo de ondas de fonte e usar o conjunto de dados sísmicos gravados sem a aplicação de um filtro de variação de tempo de domínio de freqüência variável. O painel 52 mostra um resultado RTM produzido pelo método 400. A linha 51 é uma referência horizontal, que ajuda a demonstrar a diferença no posicionamento vertical do refletor sísmico. Este refletor está localizado na profundidade correta no painel 52; no painel 50, o uso do operador de modelação FD de baixa ordem e os dados sísmicos registados não filtrados resultaram em artefatos que fazem o refletor parecer ser mais raso do que deve ser e aparentemente alterou a fase da onda pequena. Os efeitos também podem ser vistos no pacote de refletores sísmicos profundos nas imagens. A diferença na profundidade entre os dois painéis é da ordem de dezenas de metros, o que é uma diferença significativa quando se tenta assentar um poço.

[0049] Um sistema 600 para realizar a presente invenção é ilustrado esquematicamente na Figura 6. O sistema inclui uma fonte de dados 60, que, para o método 400, pode conter um conjunto de dados sísmicos gravados. Para o método 100, a fonte de dados pode conter um conjunto de dados sísmicos sintéticos gerados por um operador de modelação FD de baixa ordem. A fonte de dados está em comunicação com o processador do computador 62. O processador 62 está configurado para receber os dados de para executar os módulos compilados a partir do código legível por computador. Estes módulos podem incluir o módulo de transformação de domínio 65, que podem ser capazes de transformar os dados do domínio do tempo para o domínio da freqüência e do domínio da freqüência para o domínio do tempo. A transformação pode ser realizada, por exemplo, por uma transformada rápida de Fourier inversa e para a frente. Os módulos podem também incluir o módulo de aplicação de filtro 66, que aplica um filtro de variação de tempo de domínio de freqüência. Para o método 400, os módulos podem também incluir o módulo de migração reversa no tempo 67. O processador 62 está em comunicação com a interface de usuário 69. A interface de usuário 69 pode ser usada tanto para exibir dados e produtos de dados processados e para permitir ao usuário selecionar entre as opções para implementar aspectos do método. Os produtos de dados processados do processador 62 podem ser armazenados em fonte de dados 60.

[0050] Enquanto que na descrição precedente a presente invenção tenha sido descrita em relação a certas modalidades preferidas da mesma, e muitos pormenores foram apresentados para fins de ilustração, será evidente para os peritos na técnica que o invento é suscetível a modificações e certos outros detalhes aqui descritos podem variar consideravelmente sem se afastarem dos princípios de base da invenção. Além disso, deve se notar que as características estruturais ou passos do método mostrados ou descritos em qualquer modalidade aqui pode ser utilizado em outras modalidades também.

Claims

1. Método implementado por computador para estimar a dispersão temporal numa propagação de onda sísmica de diferenças finitas de baixa ordem, caracterizado pelo fato de que compreende: a. transformação de um conjunto de dados sísmicos do domínio de tempo para o domínio de freqüência para obter um conjunto de dados sísmicos do domínio da freqüência; b. aplicação de um filtro de variação do tempo de domínio da freqüência para o conjunto de dados sísmicos do domínio de freqüência para obter um conjunto de dados sísmicos de freqüência de domínio filtrados c. transformação do conjunto de dados sísmicos de domínio da freqüência filtrados do domínio de freqüência para o domínio do tempo para obter um conjunto de dados sísmicos filtrados do domínio do tempo.

2. Método de acordo com a reivindicação 1, caracterizado pelo fato de que o filtro de variação de tempo de domínio de freqüência variável é baseado numa velocidade de fase efetiva.

3. Método de acordo com a reivindicação 1, caracterizado pelo fato de que o conjunto de dados sísmicos é um conjunto de dados sísmicos sintéticos gerados por um operador de modelação de diferenças finitas de baixa ordem.

4. Método de acordo com a reivindicação 3, caracterizado pelo fato de que um operador de modelação de diferenças finitas de baixa ordem é um operador de modelação de diferenças finitas de segunda ordem.

5. Método de acordo com a reivindicação 1, caracterizado pelo fato de que o conjunto de dados sísmicos é um conjunto de dados sísmicos gravados.

6. Método de acordo com a reivindicação 5, caracterizado pelo fato de que compreende ainda o uso do conjunto de dados sísmicos filtrados do domínio do tempo como entrada para a migração reversa no tempo.

7. Método de acordo com a reivindicação 6, caracterizado pelo fato de que o tempo de migração reversa utiliza os dados sísmicos do domínio de freqüência filtrados como entrada.

8. Método de acordo com a reivindicação 6, caracterizado pelo fato de que o tempo de migração reversa utiliza um operador de modelação de diferenças finitas de baixa ordem.

9. Método de acordo com a reivindicação 8, caracterizado pelo fato de que o operador de modelação de diferenças finitas de baixa ordem é um operador de modelação de diferenças finitas de segunda ordem.

10. Sistema para estimar a dispersão temporal em métodos de processamento sísmicos utilizando modelação de diferenças finitas de baixa ordem, caracterizado pelo fato de que compreende: a. um dispositivo para fornecimento de informação representativa da área de subsuperfície de interesse; e b. um processador de computador, em comunicação com o dispositivo e configurado para receber os dados e para executar um código executável por computador que responde aos dados, o código executável por computador compreende: i. um módulo de transformação de domínio para a transformação de domínio de tempo para o domínio de freqüência e/ou a partir do domínio da freqüência para o domínio do tempo, e ii. um módulo de aplicação de filtro de domínio de freqüência para aplicar um filtro de variação de tempo de domínio de freqüência.

11. Sistema de acordo com a reivindicação 10, caracterizado pelo fato de que compreende ainda uma interface de usuário.

12. Sistema de acordo com a reivindicação 10, caracterizado pelo fato de que a informação representativa da área de subsuperfície de interesse compreende um conjunto de dados sísmicos sintéticos gerados por um operador de modelação de diferenças finitas de baixa ordem.

13. Sistema de acordo com a reivindicação 10, caracterizado pelo fato de que a informação representativa da área de subsuperfície de interesse compreende um conjunto de dados sísmicos gravados.

14. Sistema de acordo com a reivindicação 10, caratacrizado pelo fato de que compreende adicionalmente um módulo de migração reversa no tempo.

15. Sistema de acordo com a reivindicação 14, caracterizado pelo fato de que o módulo de migração inversa utiliza um operador de modelação de diferenças finitas de baixa ordem.