Resumos
O algoritmo desenvolvido neste trabalho é baseado no método das diferenças finitas aplicado à equação acústica da onda, assumindo que a Terra se comporta como um meio acústico, que permite implementar numericamente uma modelagem sísmica empregando malhas regulares em modelos que representam meios geológicos bidimensionais. As derivadas segundas presentes na equação da onda foram obtidas por expansões da série de Taylor de quarta ordem, para o caso espacial, e de segunda, para o temporal. Os modelos foram representados por dois tipos distintos de parametrizações: por blocos e polinomial trigonométrica. São apresentadas simulações do fenômeno de propagação de ondas em vários modelos, representados pelas duas parametrizações mencionadas, possibilitando, assim, a geração de sismogramas sintéticos para comparação. Os sismogramas obtidos, com a parametrização polinomial, mostram alguns comportamentos não desejados como: produção de reflexões artificiais, atenuação de reflexões, supressão de difrações, e pequenas alterações nos tempos de trânsito e amplitudes calculados. Por outro lado, como vantagens desta parametrização proposta, temos: economia de espaço de memória (pois um modelo complicado de velocidade pode ser representado por uma pequena quantidade de coeficientes, isto é: o modelo não é armazenado em um arquivo, mas é comprimido, ou contido, numa fórmula matemática ou numa representação analítica); produção de modelos de velocidades suavizadas (úteis na migração sísmica); sintetização do modelo num pequeno conjunto de coeficientes que se tornam parâmetros do modelo a serem estimados em possíveis procedimentos de inversão (independentemente do grau de complicação da geometria do modelo e da variação da velocidade sísmica); melhor representação de modelos realísticos por polinômios trigonométricos em comparação com a parametrização por blocos; e redução das ambigüidades normalmente presente nas técnicas de inversão, se uma quantidade suficiente de dados for disponível. Cabe acrescentar que: o efeito de Gibbs, presente em representações polinomiais, é contornado, no método das diferenças finitas, através da escolha conveniente dos nós da malha ou pelo aumento do número de coeficientes do polinômio.
Campo de Velocidade; Parametrização; Equação da Onda; Método das Diferenças Finitas; Sismogramas
The algorithm developed, in this work, is based on the finite difference method applied to the acoustic wave equation, assuming that the Earth has an acoustic behavior, that allow us to implement, in a numerical way, a seismic modeling employing regular grids in models representing two-dimensional geological media. Second derivatives with respect to space and time of wave equation are obtained by Taylor series expansion of fourth and second orders, respectively. The models are represented by two different kinds of parameterizations: blocks and trigonometric polynomials. Simulations of the wave propagation phenomenon are accomplished in several models represented by the two mentioned parameterizations, making possible to generate synthetic seismograms to be compared. Seismograms, obtained when the polynomial parameterization is used, show some undesired behaviors such as: production of artificial reflections, weakening of reflection events, suppression of diffractions, and a little alteration in calculated traveltimes and amplitudes. By the other hand, some advantages of the proposed polynomial parameterization are: economy of computer memory space (because a complicated velocity model can be represented by a few quantity of coefficients, it means: the model is not more stored in a file, but it is compressed, or contained, in a mathematical formula or an analytical representation); production of a smooth velocity model (useful to generate a time field to be used in seismic migration, for example, in reverse time migration); the polynomial coefficients are just the parameters of the model to be estimated by an inversion procedure (independently of the degree of geometrical complication of the model and how is varying the seismic velocity on it); realistic models are better represented by trigonometric polynomial than block parameterization; and the ambiguity commonly present in inversion techniques results, is reduced if sufficient data is available. Gibbs effects, present in polynomial representations, is avoided by the finite differences method by choosing conveniently the knots of the mesh or increasing the number of polynomial coefficients.
Velocity Field; Parameterization; Wave Equation; Finite-Difference Method; Seismograms
Modelagem acústica bidimensional usando diferentes parametrizações de campos de velocidades
Roberto Hugo M. dos Santos; Wilson M. Figueiró
CPGG-UFBA, Centro de Pesquisa em Geofísica e Geologia da Universidade Federal da Bahia, LAGEP, Instituto de Geociências, Rua Barão de Jeremoabo, sn, 40170-290 Federação, Salvador, BA, Brasil. Fones: (55)(71) 3203-8522 e 3203-8520; Fax: (55)(71) 3203-8501 - E-mails: rhms@cpgg.ufba.br e figueiro@cpgg.ufba.br
RESUMO
O algoritmo desenvolvido neste trabalho é baseado no método das diferenças finitas aplicado à equação acústica da onda, assumindo que a Terra se comporta como um meio acústico, que permite implementar numericamente uma modelagem sísmica empregando malhas regulares em modelos que representam meios geológicos bidimensionais. As derivadas segundas presentes na equação da onda foram obtidas por expansões da série de Taylor de quarta ordem, para o caso espacial, e de segunda, para o temporal. Os modelos foram representados por dois tipos distintos de parametrizações: por blocos e polinomial trigonométrica. São apresentadas simulações do fenômeno de propagação de ondas em vários modelos, representados pelas duas parametrizações mencionadas, possibilitando, assim, a geração de sismogramas sintéticos para comparação. Os sismogramas obtidos, com a parametrização polinomial, mostram alguns comportamentos não desejados como: produção de reflexões artificiais, atenuação de reflexões, supressão de difrações, e pequenas alterações nos tempos de trânsito e amplitudes calculados. Por outro lado, como vantagens desta parametrização proposta, temos: economia de espaço de memória (pois um modelo complicado de velocidade pode ser representado por uma pequena quantidade de coeficientes, isto é: o modelo não é armazenado em um arquivo, mas é comprimido, ou contido, numa fórmula matemática ou numa representação analítica); produção de modelos de velocidades suavizadas (úteis na migração sísmica); sintetização do modelo num pequeno conjunto de coeficientes que se tornam parâmetros do modelo a serem estimados em possíveis procedimentos de inversão (independentemente do grau de complicação da geometria do modelo e da variação da velocidade sísmica); melhor representação de modelos realísticos por polinômios trigonométricos em comparação com a parametrização por blocos; e redução das ambigüidades normalmente presente nas técnicas de inversão, se uma quantidade suficiente de dados for disponível. Cabe acrescentar que: o efeito de Gibbs, presente em representações polinomiais, é contornado, no método das diferenças finitas, através da escolha conveniente dos nós da malha ou pelo aumento do número de coeficientes do polinômio.
Palavras-chave: Campo de Velocidade, Parametrização, Equação da Onda, Método das Diferenças Finitas, Sismogramas.
ABSTRACT
The algorithm developed, in this work, is based on the finite difference method applied to the acoustic wave equation, assuming that the Earth has an acoustic behavior, that allow us to implement, in a numerical way, a seismic modeling employing regular grids in models representing two-dimensional geological media. Second derivatives with respect to space and time of wave equation are obtained by Taylor series expansion of fourth and second orders, respectively. The models are represented by two different kinds of parameterizations: blocks and trigonometric polynomials. Simulations of the wave propagation phenomenon are accomplished in several models represented by the two mentioned parameterizations, making possible to generate synthetic seismograms to be compared. Seismograms, obtained when the polynomial parameterization is used, show some undesired behaviors such as: production of artificial reflections, weakening of reflection events, suppression of diffractions, and a little alteration in calculated traveltimes and amplitudes. By the other hand, some advantages of the proposed polynomial parameterization are: economy of computer memory space (because a complicated velocity model can be represented by a few quantity of coefficients, it means: the model is not more stored in a file, but it is compressed, or contained, in a mathematical formula or an analytical representation); production of a smooth velocity model (useful to generate a time field to be used in seismic migration, for example, in reverse time migration); the polynomial coefficients are just the parameters of the model to be estimated by an inversion procedure (independently of the degree of geometrical complication of the model and how is varying the seismic velocity on it); realistic models are better represented by trigonometric polynomial than block parameterization; and the ambiguity commonly present in inversion techniques results, is reduced if sufficient data is available. Gibbs effects, present in polynomial representations, is avoided by the finite differences method by choosing conveniently the knots of the mesh or increasing the number of polynomial coefficients.
Keywords: Velocity Field, Parameterization, Wave Equation, Finite-Difference Method, Seismograms.
INTRODUÇÃO
Métodos de soluções numéricas para a equação da onda foram estudadas por Mitchell (1969) e Boor (1972). Alford et al. (1974) comparou os resultados dos trabalhos usando diferenças finitas (DF) com resultados analíticos em modelos com quinas. Estes métodos atingem um alto grau de aproximação, desde que se tenha uma grade suficientemente fina.
Vista como solução do problema direto na metodologia sísmica, a modelagem sísmica numérica pode ser empregada em: formulação da inversão sísmica não-linear, que requer uma teoria para a simulação de ondas sísmicas; processos de migração; auxílio à interpretação sísmica de regiões com geologia complexa; e geração de dados para testes em algoritmos de processamento.
Tradicionalmente, utiliza-se simplificações da equação da onda, como aproximações assintóticas e linearizações. Por exemplo: a teoria do raio ( et al., 1977; Keller, 1962) tem como base a equação eiconal, a qual é uma aproximação de alta freqüência (ótica geométrica), sendo inadequada quando o meio é altamente heterogêneo, apresentando problemas em modelos contendo interfaces com curvaturas abruptas (Trorey, 1970; Hilterman, 1970), entretanto com base no teorema de Green, ela permite a geração de resultados mais abrangentes, sendo útil para a modelagem de uma reflexão específica (evento). Como outro exemplo, temos os métodos espectrais, que apresentam limitações para modelos com geometria complexa e apresentam fortes gradientes de propriedades físicas, como acontece em contatos de litologias distintas.
O uso do método das diferenças finitas, MDF, (Alterman & Karal, 1968; Ottaviane, 1971) permite o desenvolvimento de técnicas com base na equação da onda sem qualquer aproximação ou limitação na geometria dos modelos. Desta forma, o método leva em conta não apenas as ondas diretas e refletidas (primárias e múltiplas), como também, ondas: superficiais, refratadas, convertidas, difratadas, e ondas observadas nas zonas de sombra da teoria do raio. O MDF permite, ainda, a obtenção de fotografias ("snapshots") do campo de ondas em qualquer instante da propagação.
O MDF usa uma discretização do meio e, consequentemente, discretiza o campo de onda (amplitudes) que atravessa tal meio, usando uma malha de pontos. Esse método aproxima as derivadas espaciais e temporais por operadores de DF obtidos por uma expressão de Taylor truncada, usando pontos vizinhos em diferenças ponderadas. Seus resultados, tão precisos quanto desejado, apresentam acurácia muito superior em comparação aos outros métodos de modelagem sísmica.
EQUAÇÃO ACÚSTICA DA ONDA
Para os objetivos deste trabalho, é suficiente considerar meios acústicos na realização do estudo da influência que diferentes tipos de parametrizações do meio exercem na modelagem. Características distintivas básicas devem ser investigadas, preliminarmente em modelos representativos de meios menos complicados, além disso, espera-se que tais características se tornem cada vez mais acentuadas com o aumento do grau de complicação do meio e, consequentemente, do modelo que o representa.
Uma equação que relaciona a pressão com sua velocidade de propagação em meios acústicos é a chamada equação da onda de pressão e é dada por:
onde p = p(x, t) representa a pressão, x é o vetor posição, t é o tempo, c é a velocidade de propagação da onda de pressão em (x, t), r é a densidade em (x, t), Ñ2 o laplaciano relativo às dimensões espaciais, o operador gradiente, e ln o logaritmo neperiano. Se r é constante, então a Equação 1 ganha a forma:
Esta equação é chamada de equação da onda escalar acústica. Temos que a relação pressão-deslocamento é dada por:
onde x = (x, z), k é o módulo de compressão volumétrica, e u(x, t) = (u1(x, t),u2(x, t)) é o vetor deslocamento do ponto material situado em x no instante t. Substituindo-se (3) em (2), chega-se à seguinte equação:
sendo assim,
e
para f1 e f2 tais que:
Portanto, vamos considerar a equação escalar da onda, no caso bidimensional acústico, como:
onde u = u2(x, z, t) representa o deslocamento vertical do ponto material do meio relativamente à sua posição de equilíbrio x = (x, z) no instante t, c = c(x, z) a velocidade de propagação da onda em (x, z), e f(x, z, t) a função fonte dada por:
f(x, z, t) = 0 se (x, z) ¹ (xs, 0) e f(xs, 0, t) = f(t)
com xs sendo a coordenada x da posição da fonte.
MÉTODO DAS DIFERENÇAS FINITAS
Para substituir uma equação diferencial por uma equação que envolva somente diferenças e quocientes finitos (que não envolvem valores infinitamente grandes ou pequenos), que é o cerne do método, é necessário obter expressões que permitam calcular numericamente as derivadas. Isto pode ser feito via série de Taylor, que nos permite escrever, para a função u = u(x), as seguintes expansões:
e
onde u(n)(x) representa a n-ésima derivada de u relativamente a x.
Rearranjando os termos das Equações (9) e (10), obtemos então as seguintes séries para u(1)(x):
e
Truncando-se as Equações (11) e (12) a partir dos termos nos quais Dx tem grau maior ou igual a 1, obtemos as seguintes aproximações para a derivada primeira:
e
com ui+1 = u(x + Dx), ui = u(x), ui-1 = u(x - Dx) e h = Dx.
Podemos obter uma aproximação com menor erro de truncamento subtraindo-se a Equação (10) da Equação (9) e explicitando u(1)(x):
Truncando a série acima a partir do termo com h2, obtemos a seguinte aproximação para u(1)(x):
Podemos obter a derivada segunda subtraindo-se a Equação (12) da Equação (11) e explicitando u(2)(x), assim:
Usando as Equações (9) e (10) obtemos a seguinte expressão para a derivada segunda:
onde ui+2 = u(x + 2Dx), ui-2 = u(x - 2Dx) e (uxx)i é a derivada parcial segunda da função u com respeito a x, calculada no nó i da malha.
Usando a Equação (17), obtemos a seguinte aproximação para as derivadas segundas temporais:
onde (utt)k é a derivada segunda da função u relativamente à variável t, calculada em t = kDt, onde Dt é o intervalo temporal.
Os operadores usados em DF executam uma ação específica sobre a função nos pontos da malha.
Em geral, quanto mais precisas as aproximações para as derivadas, mais pontos vizinhos são requeridos e, portanto, mais dispendioso torna-se o cálculo. No entanto, não só a precisão tem de ser levada em consideração, devemos, também, considerar critérios de estabilidade que são importantes para garantir convergência, que é assegurada por determinada largura da malha espacial, assim como pelo intervalo temporal a ser usado.
As expressões (18) e (19) são, respectivamente, chamadas de aproximações de quarta ordem espacial e segunda ordem temporal da derivada segunda. A ordem da aproximação depende do número de pontos nas quais a função é amostrada. Tais pontos estão regularmente e simetricamente espaçados em relação àquele no qual se deseja conhecer a derivada segunda da função considerada.
PROCEDIMENTO RECURSIVO DO MDF
Utilizando-se as aproximações para as derivadas segundas representadas pelas Equações (18) e (19), obtemos a seguinte expressão recursiva necessária à implementação computacional do MDF aplicado à equação da onda:
onde os índices i, j, e k representam, respectivamente, as variáveis x, z e t; e g = -(c · Dt)2 f.
CONDIÇÃO PARA ESTABILIDADE E NÃO DISPERSÃO NUMÉRICA
No MDF, o cálculo das derivadas envolvidas nas equações diferenciais são realizados segundo aproximações com maior ou menor precisão, o que gera erro no resultado numérico. No caso particular da equação da onda, este erro se apresenta sobre a forma de dispersão numérica. Este problema geralmente se manifesta gerando oscilações intermitentes na forma temporal do pulso sísmico. Outro problema causado é a presença de ruído em um determinado ponto antes da chegada do sinal.
Para que haja estabilidade, ou seja, para que o erro na determinação das amplitudes das ondas permaneça pequeno e não cresça com o tempo, faz-se necessário calcular o intervalo de amostragem temporal (Dt) máximo que seja mais conveniente.
Segundo Cunha & Mufti (1997), a modelagem por DF emprega freqüentemente malha uniforme (Dx = Dz), que, visando evitar dispersão ou instabilidade numéricas, deve atender às relações que seguem:
e
onde Dx e Dz são os intervalos espaciais; Dt é o intervalo temporal utilizado; cmin, cmax são, respectivamente, as velocidades mínima e máxima presentes no modelo; fmax é a freqüência máxima da fonte; e as constantes a e b são parâmetros que dependem da ordem de aproximação usada para cálculo das derivadas espaciais e temporais. A qualidade do resultado melhora com o aumento de a.
FUNÇÃO FONTE
Um conceito muito discutido quando se fala em resolução vertical é o poder de resolução, Pa, do pulso sísmico. Fisicamente é uma medida do quanto um pulso de banda limitada tende a se aproximar de uma função impulso perfeita. Para estudos de resolução é mais conveniente que seja usado umpulso de fase zero. Este pode ser dado diretamente no domínio do tempo ou pode ser definido a partir de seu espectro de freqüências. Porém, é conveniente que seja uma função bem comportada, com espectro limitado ou que decaia rapidamente a partir de um determinado ponto, de modo que uma freqüência dominante, assim como a de corte, fiquem bem caracterizadas. É bom enfatizar ser necessário estrito controle sobre o conteúdo de freqüência da fonte quando se trabalha com DF. Desta forma, se pode saber qual a largura da malha que mantém a dispersão sob controle. Neste trabalho foi utilizada como função fonte o pulso que possui a forma da segunda derivada da função gaussiana, definida pela fórmula:
onde fp é a freqüência de pico. A Figura 1 ilustra a função f(t), que foi usada como fonte.
Nos experimentos realizados neste trabalho a fonte foi colocada na posição 550m na superfície de observação com uma freqüência máxima de 80Hz.
REFLEXÃO NAS BORDAS
Qualquer aproximação por DF, bem como qualquer outro método de solução numérica da equação da onda, enfrenta problemas com as bordas do modelo. Uma solução, no caso de fonte pontual, é a aproximação explícita, que consiste em estender os limites do modelo, de modo que estas reflexões indesejadas não tenham tempo de estar presente no sismograma. Porém, este método é muito caro computacionalmente. Outra maneira de atacar este problema é tentar anular as reflexões impondo-se condições de contorno apropriadas nas referidas bordas. O esquema aqui utilizado, Cerjan et al. (1985), aplica, a cada iteração, uma diminuição gradativa das amplitudes em uma faixa nas proximidades de cada borda. Pontos a uma distância i da borda tem sua amplitude P0 reduzida para P através de:
onde F é o fator de atenuação; N é o número de pontos da malha entre a fonte e a borda do modelo; e i é a distância, em pontos da malha, até a borda do modelo. Cerjan et al. (1985) sugere N = 20 e F = 0,015.
DEFINIÇÃO DE MODELOS SÍSMICOS 2-D
Os modelos considerados têm dimensão horizontal de 1100m e profundidade de 600m sendo os modelos sísmicos com graus de complexibilidade crescente na seguinte seqüência: MI (de duas camadas homogêneas separadas por uma interface horizontal, Figura 2), modelo MII (de seis camadas homogêneas, Figura 3), MIII (apresenta variações laterais e verticais, Figura 4), MIV (constituído por quatro camadas e dotado de uma intrusão, Figura 5), MV (de três camadas com uma falha normal, Figura 6), e MVI (de três camadas e uma bacia, Figura 7). Sismogramas são gerados para cada modelo, que são representados por duas diferentes parametrizações: por blocos (PB) e pela polinomial trigonométrica (PPT).
PARAMETRIZAÇÕES DE CAMPOS DE VELOCIDADES SÍSMICAS
Divisão do meio em Células
A utilização de células (ou blocos) na representação do modelo é o tipo de parametrização mais utilizado na modelagem sísmica, que é, simplesmente, a divisão do modelo em blocos (ou células), onde cada um possui uma velocidade constante. As Figuras 8, e 9 exibem, respectivamente, os modelos MIV, e MVI; parametrizados por blocos. Nesta parametrização o número de blocos totalizam 19 × 26 = 494, onde cada um tem, aproximadamente, 58m de comprimento e 23m de largura. Cada bloco é homogêneo, isotrópico, e com velocidade igual àquela da região do modelo por ele coberta.
Polinômio Trigonométrico
A parametrização polinomial trigonométrica consiste em representar o modelo geofísico através de uma única função polinomial trigonométrica. Tal representação é obtida através da particularização do teorema que fornece séries baseadas em funções ortogonais para funções de duas variáveis (Kreider et al., 1972). Isto é, o campo de velocidade sísmica V(x, z) é representado por:
onde
onde R é a região espacial na qual o modelo é considerado, e l e d representam, respectivamente, o comprimento e a profundidade máximas do modelo. Tal função representa o campo de velocidades do modelo no qual a onda se propaga.
Para todos os modelos sísmicos parametrizados pelo polinômio trigonométrico foram utilizados 80 coeficientes. Para modelos mais simples tal número pode parecer um exagero, entretanto, à medida que vamos tendo modelos mais complicados, ou seja, com um número cada vez maior de parâmetros, o número de coeficientes utilizados vai ajustar-se à necessidade de representação desses modelos. O número 80 foi mantido em todos os modelos visando uma interpretação uniformizada dos resultados.
Com o objetivo de facilitar os cálculos das integrais, foram feitas aproximações com pequenos quadrados nas proximidades das interfaces inclinadas tanto para o modelo MV como para o modelo MVI, o que nos dá a certeza de que, independentemente da complexibilidade dos modelos, este tipo de parametrização pode ser utilizada.
As Figuras 10 e 11 mostram, respectivamente, os modelos MIV, e MVI; parametrizados por polinômios trigonométricos. Neles observamos, no interior das regiões homogêneas, oscilações características das funções trigonométricas. Além disso, nas interfaces podemos observar efeitos devido ao fenômeno de Gibbs.
Dentre as vantagens da parametrização polinomial temos: economia de espaço de memória; os polinômios podem ser, facilmente e infinitamente, derivados e integrados; todos os aspectos do modelo ficam resumidos nos coeficientes do polinômio; possíveis trabalhos de inversão ficam sintetizados na estimativa dos coeficientes do polinômio; e se o campo de velocidades for contínuo (ou possuir descontinuidades de salto finito), então existe uma seqüência de funções polinomiais trigonométricos que para ele converge.
MODELAGEM ACÚSTICA 2-D
Foram realizadas modelagens acústicas baseadas em operadores escalares gerando doze sismogramas para os seis modelos propostos utilizando-se as duas parametrizações consideradas para cada um deles. Foram implementados programas Fortran para simulação do campo de onda acústico. Tal simulação foi realizada adotando-se: o MDF com aproximação de segunda ordem no tempo e quarta ordem no espaço, malhas regulares, condições de contorno não reflexivas (versão "one-way" da equação da onda) conjugadas com bordas de absorção como propostas por Cerjan et al. (1985), e termo fonte proveniente da segunda derivada da função gaussiana, como definido em Cunha & Mufti (1997). Além disso, para garantir estabilidade numérica e evitar dispersão de energia, os intervalos de amostragem foram definidos por 10 amostras por comprimento de onda (temporal e espacial).
Para os modelos considerados na seqüência de MI a MVI, as Figuras 12, 13, 14, 15, 16, e 17 mostram em (a) os sismogramas sintéticos obtidos por modelagem acústica aplicando o MDF à equação da onda, usando blocos (PB), enquanto em (b), mostram aqueles produzidos quando os modelos são representados pela função polinomial trigonométrica.
ANÁLISE DOS RESULTADOS
Na Figura 12 os eventos se manifestam com maior clareza em (a) do que em (b). A primeira reflexão na interface de MI revela um pequeno deslocamento vertical quando comparamos (a) com (b). Pode-se observar também que as oscilações do polinômio trigonométrico cria fracas reflexões.
Na Figura 13, a chegada da onda direta apresenta-se com maior clareza em (b), entretanto as primeiras reflexões em cada interface manifestam-se de modo mais nítido em (a), além disso, no sismograma (a), obtido com a PB, são visíveis as múltiplas referentes às duas primeiras interfaces, enquanto que na PPT, sismograma (b), estas múltiplas não são facilmente visíveis.
Na Figura 14, em (a), além das ondas diretas, primeiras reflexões nas interfaces e a múltipla da primeira interface, também podem ser observadas as difrações nas quinas do modelo MIII. Entretanto, em (b), estas difrações são atenuadas devido à suavização das quinas, do mesmo modelo sísmico, causada pela PPT.
Na Figura 15, todos os eventos em (a) são facilmente identificados enquanto em (b) o refletor acima da base e do topo da intrusão não são identificados.
Observa-se em (a) da Figura 16, que as reflexões nas interfaces são prolongadas pelas difrações nas quinas do modelo MV e pelas reflexões no plano de falha, ocorrendo dessa forma o cruzamento de eventos. Já em (b), onde as quinas difratoras são suavizadas, não ocorre este prolongamento e sim uma descontinuidade causada pelo plano de falha.
Na Figura 17 , os resultados, de certa forma, são semelhantes àqueles mostrados na Figura 16, pois de forma análoga o sismograma (a) apresenta eventos contínuos causados pelas difrações nas quinas do modelo MVI e reflexões nas paredes da bacia, enquanto que em (b) as difrações são atenuadas e as reflexões do fundo da bacia só são registradas nos geofones centrais.
É importante citar que nas comparações entre os sismogramas dos diferentes modelos, existe uma diferença no posicionamento dos refletores causada pelas diferentes velocidades nas camadas em cada um dos modelos, pois, diferentes parametrizações aplicadas ao mesmo modelo sísmico geram modelos parametrizados diferentes. Em outras palavras, diferentes parametrizações de um mesmo modelo podem gerar novos modelos que, no nível do detalhe, podem ser distintos. Emtodos os sismogramas gerados através da PPT, não foi solucionado, satisfatoriamente, o problema das reflexões nas bordas do modelo, apesar de ser utilizada a mesma técnica tanto no caso da PB como no da PPT, sendo que para a primeira obtivemos um bom resultado.
Como os modelos implementados neste trabalho possuem uma forte variação de velocidade, que dificulta o processo de modelagem pelo MDF, os modelos de velocidades poderiam ser suavizados, utilizando-se, por exemplo, média móvel, de modo a obtermos melhores resultados na modelagem. Entretanto, embora a parametrização possa implicar em suavização, como ocorre com a PPT, o contrário não ocorre. Isto é, a suavização não é um modo de representação de modelos, portanto está fora do objetivo deste trabalho, que consiste em avaliar os efeitos que diferentes representações (parametrizações), de um mesmo modelo, produzem na modelagem.
A PPT exibe ondulações na representação do modelo, isto produz falsos eventos no sismograma. Por outro lado, o uso da PB gera sismogramas com eventos bem definidos. Isto pode nos fazer pensar que a PB é superior a PPT na representação dos modelos, mas não podemos perder de vista o fato dos modelos usados neste trabalho serem esquemáticos (com extensas regiões homogêneas). Para estes modelos, como este trabalho mostra, a PB parece ser mais adequada. Entretanto, para modelos mais realísticos não podemos dizer o mesmo, pois a PB pode não ser capaz de atender ao elevado grau de complexidade do modelo e, além disso, neste caso, não existem garantias de que o sismograma exibirá eventos menos espúrios. Deve-se acrescentar ainda, que as oscilações da PPT podem ser reduzidas melhorando-se a acurácia dos cálculos para a determinação dos coeficientes ou aumentando-se o número destes, principalmente se o modelo for contínuo ou com fracas descontinuidades.
O uso da PPT na representação de vagarosidades pode ser menos oscilante do que na representação de velocidades, desde que estas nunca sejam menores do que 1,0 km/s, pois, neste caso, as oscilações serão insuportavelmente altas. Além disso, para velocidades altas, a vagarosidade diminui drasticamente, e ela encontra-se elevada ao quadrado no denominador da fórmula recursiva (20) usada no MDF. Isto pode causar grandes instabilidades numéricas no referido método.
CONCLUSÕES
Em termos gerais, o algoritmo implementado por diferenças finitas (DF) produziu bons resultados tanto para a parametrização por blocos (PB) quanto para a parametrização polinomial trigonométrica (PPT).
Tanto na PB como na PPT, os sismogramas obtidos apresentaram similaridades gerais nos resultados. Reflexões artificiais presentes nos sismogramas gerados com a PPT foram capazes de esconder muitos eventos que se manifestaram de forma clara quando a PB foi usada. Tais reflexões artificiais se devem ao fato da PPT apresentar oscilações.
A PPT proporciona uma maior economia de espaço de memória do computador, devido ao fato de não necessitar ter, como entrada no programa de DF, um arquivo de velocidades, pois estas, podem ser calculadas durante a execução do programa, bastando entrar com a função polinomial trigonométrica que é definida por um número relativamente pequeno decoeficientes.
Pelo menos todo o modelo heterogêneo e isotrópico pode ser representado, com o grau de exatidão desejado, pela PPT.
A PPT simplifica a representação de modelos mais complicados, de maneira que, independentemente de sua complexibilidade, estes modelos podem ser divididos em pequenos blocos, onde se realiza as integrações necessárias para obtenção dos coeficientes, como foi feito para os modelos MV e MVI.
A PPT para modelos com extensas regiões de velocidade homogênea, tais como os estudados neste trabalho, revelou-se menos adequado do que a PB. Entretanto, acredita-se que esta situação possa ser revertida ao considerar-se modelos mais complicados. Um estudo deste gênero poderá ser considerado em desenvolvimentos futuros.
AGRADECIMENTOS
Agradecemos à Agência Nacional de Petróleo (ANP) pelo apoio financeiro em forma de bolsa de estudos.
Recebido em 19 dezembro, 2005 / Aceito em 14 março, 2006
Received on December 19, 2005 / Accepted on March 14, 2006
NOTAS SOBRE OS AUTORES
Roberto Hugo Melo dos Santos graduou-se em Matemática pela UNIT-SE (Universidade Tiradentes, Sergipe) em 1998. Obteve seu mestrado em Geofísica de Exploração de Petróleo na Universidade Federal da Bahia (UFBA) em 2002. Atualmente é aluno do curso de doutorado em Geofísica da UFBA. Suas principais áreas de interesse são: Modelagem, Inversão e Processamento de Dados Sísmicos.
Wilson Mouzer Figueiró bacharelou-se em Matemática na Universidade Federal do Rio de Janeiro (UFRJ) em 1983. Ingressou como docente no Departamento de Geologia e Geofísica Aplicada (DGGA) do Instituto de Geociências (IGEO) da Universidade Federal da Bahia (UFBA) em 1992, onde, desde então, leciona Geofísica Matemática. Doutorou-se em Geofísica na UFBA em 1994. Realizou estudos de Pós-doutorado no "Département Terre-Atmosphère-Océan de L'Ecole Normale Supérieure de Paris" (DTAO-ENS), França, no período 1997-1999. É pesquisador do Centro de Pesquisa em Geofísica e Geologia (CPGG) da UFBA. Suas áreas de interesses são: Modelagem, Inversão, Traçamento de Raios Sísmicos, e Matemática Aplicada às Geociências.
- ALFORD R, KELLY K & BOOR D. 1974. Accuracy of finite-difference modeling of the acoustic wave equation. Geophysics, 6: 834-842.
- ALTERMAN Z & KARAL FC. 1968. Propagation of elastic wave in layered media by finite difference methods. Bull. Seis. Soc., 58: 367-398.
- BOOR DM. 1972. Finite-Difference Methods for Seismic Wave Propagation in Heterogeneous Materials. Academic Press, New York.
- CERJAN J, KOSLOFF D, KOSLOFF R & RESHEF M. 1985. A nonreflecting boundary condition for discrete acoustic and elastic wave equations. Geophysics, 50: 705-708.
- CERVENÝ V, MOLOTOV I & PSENCÍK I. 1977. Ray Method in Seismology. Univ. Karlova, Praha.
- CUNHA PEM & MUFTI I. 1997. Um Método Eficiente para Migraçăo Pré-Empilhamento 3-D em Profundidade Reversa no Tempo pelo Método das Diferenças Finitas. Diss. de Mestrado, Universidade Federal da Bahia, Salvador, Brasil.
- HILTERMAN FJ. 1970. Three-dimensional seismic modeling. Geophysics, 35: 1020-1037.
- KELLER JB. 1962. Geometrical theory of diffraction. J. Opt. Soc. Am., 52: 116-130.
- KREIDER D, KULLER RC, OSTBERG DR & PERKINS FW. 1972. Introduçăo ŕ Análise Linear 2. Editora Universal de Brasília, Rio de Janeiro, 436-442.
- MITCHELL AR. 1969. Computational Methods in Partial Differential Equations. John Wiley and Sons, London.
- OTTAVIANE M. 1971. Elastic wave propagation in two evenly welded quarter spaces. Bull. Seis. Soc. Am., 61: 1119-1152.
- TROREY AW. 1970. A simple theory for seismic diffractions. Geophysics, 35: 762-784.
Datas de Publicação
-
Publicação nesta coleção
15 Maio 2007 -
Data do Fascículo
Mar 2006
Histórico
-
Recebido
19 Dez 2005 -
Aceito
14 Mar 2006