Migração FD 3D em profundidade usando aproximação de Padé complexa
Conventional implementations of 3D finite-difference (FD) migration use splitting techniques to accelerate performace and save computational cost. However, such techniques are plagued with numerical anisotropy that jeopardizes correct positioning of dipping reflectors in the directions not used f...
| Autor principal: | COSTA, Carlos Alexandre Nascimento da |
|---|---|
| Grau: | Dissertação |
| Idioma: | por |
| Publicado em: |
Universidade Federal do Pará
2014
|
| Assuntos: | |
| Acesso em linha: |
http://repositorio.ufpa.br/jspui/handle/2011/5974 |
| id |
ir-2011-5974 |
|---|---|
| recordtype |
dspace |
| spelling |
ir-2011-59742018-01-03T12:54:30Z Migração FD 3D em profundidade usando aproximação de Padé complexa COSTA, Carlos Alexandre Nascimento da COSTA, Jessé Carvalho http://lattes.cnpq.br/7294174204296739 Equação de onda Diferenças finitas Algoritmos Sismologia CNPQ::CIENCIAS EXATAS E DA TERRA::GEOCIENCIAS::GEOFISICA::SISMOLOGIA Conventional implementations of 3D finite-difference (FD) migration use splitting techniques to accelerate performace and save computational cost. However, such techniques are plagued with numerical anisotropy that jeopardizes correct positioning of dipping reflectors in the directions not used for the splitted operators. We implement 3D downward continuation migration without splitting in the space coordinates using a complex Padé approximation and implicit finite differences. In this way, the numerical anisotropy is eliminated at the expense of a computationally more intensive solution of a large banded linear system. We compare the performance of the iterative stabilized biconjugate gradient (Bi-CGSTAB) and the multifrontal massively parallel direct solver (MUMPS). It turns out that the use of the complex Padé approximation provides an effective preconditioner for the Bi-CGSTAB, reducing the number of iterations relative to the real Padé expansion of the FD operator. As a consequence, the iterative Bi-CGSTABmethod ismore efficient than the directMUMPSmethodwhen solving for a single termin the Padé expansion. Forwide angle approximations direct methods are required. These algorithms are validated and the properties evaluated computing themigration impulse response in the SEG/EAGE salt model. CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico A implementação convencional do método de migração por diferenças finitas 3D, usa a técnica de splitting inline e crossline para melhorar a eficiência computacional deste algoritmo. Esta abordagem torna o algoritmo eficiente computacionalmente, porém cria anisotropia numérica. Esta anisotropia numérica por sua vez, pode levar a falsos posicionamentos de refletores inclinados, especialmente refletores com grandes ângulos de mergulho. Neste trabalho, como objetivo de evitar o surgimento da anisotropia numérica, implementamos o operador de extrapolação do campo de onda para baixo sem usar a técnica splitting inline e crossline no domínio frequência-espaço via método de diferenças finitas implícito, usando a aproximação de Padé complexa. Comparamos a performance do algoritmo iterativo Bi-gradiente conjugado estabilizado (Bi-CGSTAB) com o multifrontal massively parallel solver (MUMPS) para resolver o sistema linear oriundo do método de migração por diferenças finitas. Verifica-se que usando a expansão de Padé complexa ao invés da expansão de Padé real, o algoritmo iterativo Bi-CGSTAB fica mais eficientes computacionalmente, ou seja, a expansão de Padé complexa atua como um precondicionador para este algoritmo iterativo. Como consequência, o algoritmo iterativo Bi-CGSTAB é bem mais eficiente computacionalmente que o MUMPS para resolver o sistema linear quando usado apenas um termo da expansão de Padé complexa. Para aproximações de grandes ângulos, métodos diretos são necessários. Para validar e avaliar as propriedades desses algoritmos de migração, usamos o modelo de sal SEG/EAGE para calcular a sua resposta ao impulso. 2014-11-03T12:25:54Z 2014-11-03T12:25:54Z 2011-06-21 Dissertação COSTA, Carlos Alexandre Nascimento da. Migração FD 3D em profundidade usando aproximação de Padé complexa. 2011. 74 f. Dissertação (Mestrado) - Universidade Federal do Pará, Instituto de Geociências, Belém, 2011. Programa de Pós-Graduação em Geofísica. http://repositorio.ufpa.br/jspui/handle/2011/5974 por Acesso Aberto application/pdf Universidade Federal do Pará Brasil Instituto de Geociências UFPA Programa de Pós-Graduação em Geofísica |
| institution |
Repositório Institucional - Universidade Federal do Pará |
| collection |
RI-UFPA |
| language |
por |
| topic |
Equação de onda Diferenças finitas Algoritmos Sismologia CNPQ::CIENCIAS EXATAS E DA TERRA::GEOCIENCIAS::GEOFISICA::SISMOLOGIA |
| spellingShingle |
Equação de onda Diferenças finitas Algoritmos Sismologia CNPQ::CIENCIAS EXATAS E DA TERRA::GEOCIENCIAS::GEOFISICA::SISMOLOGIA COSTA, Carlos Alexandre Nascimento da Migração FD 3D em profundidade usando aproximação de Padé complexa |
| topic_facet |
Equação de onda Diferenças finitas Algoritmos Sismologia CNPQ::CIENCIAS EXATAS E DA TERRA::GEOCIENCIAS::GEOFISICA::SISMOLOGIA |
| description |
Conventional implementations of 3D finite-difference (FD) migration use splitting techniques
to accelerate performace and save computational cost. However, such techniques
are plagued with numerical anisotropy that jeopardizes correct positioning of dipping reflectors
in the directions not used for the splitted operators. We implement 3D downward
continuation migration without splitting in the space coordinates using a complex Padé
approximation and implicit finite differences. In this way, the numerical anisotropy is
eliminated at the expense of a computationally more intensive solution of a large banded
linear system. We compare the performance of the iterative stabilized biconjugate gradient
(Bi-CGSTAB) and the multifrontal massively parallel direct solver (MUMPS). It turns out
that the use of the complex Padé approximation provides an effective preconditioner for
the Bi-CGSTAB, reducing the number of iterations relative to the real Padé expansion of
the FD operator. As a consequence, the iterative Bi-CGSTABmethod ismore efficient than
the directMUMPSmethodwhen solving for a single termin the Padé expansion. Forwide
angle approximations direct methods are required. These algorithms are validated and
the properties evaluated computing themigration impulse response in the SEG/EAGE salt
model. |
| author_additional |
COSTA, Jessé Carvalho |
| author_additionalStr |
COSTA, Jessé Carvalho |
| format |
Dissertação |
| author |
COSTA, Carlos Alexandre Nascimento da |
| title |
Migração FD 3D em profundidade usando aproximação de Padé complexa |
| title_short |
Migração FD 3D em profundidade usando aproximação de Padé complexa |
| title_full |
Migração FD 3D em profundidade usando aproximação de Padé complexa |
| title_fullStr |
Migração FD 3D em profundidade usando aproximação de Padé complexa |
| title_full_unstemmed |
Migração FD 3D em profundidade usando aproximação de Padé complexa |
| title_sort |
migração fd 3d em profundidade usando aproximação de padé complexa |
| publisher |
Universidade Federal do Pará |
| publishDate |
2014 |
| url |
http://repositorio.ufpa.br/jspui/handle/2011/5974 |
| _version_ |
1832604161496055808 |
| score |
11.755432 |