ESTRATÉGIA HÍBRIDA DE OTIMIZAÇÃO SEM DERIVADAS PARA A INVERSÃO COMPLETA DA FORMA DA ONDA
DOI:
https://doi.org/10.15628/holos.2017.5445Palavras-chave:
Inversão sísmica, Inversão completa da forma da onda, Otimização sem derivadasResumo
Na exploração sísmica investigamos as características de subsuperfície usando técnicas de inversão completa da forma da onda (Full Waveform Inversion - FWI) a qual foi abordada como um problema de otimização não linear. A técnica FWI tradicionalmente usa métodos matemáticos baseados em derivadas e, portanto, falha quando a função objetiva é não diferenciável. Ademais, isto acarreta um alto custo computacional e uma precisão limitada a mínimos locais. Portanto, neste trabalho foi adotada uma metodologia sem derivadas, Derivative-Free Optimization (DFO) para encontrar o mínimo global. Neste tipo de abordagem utilizamos as técnicas Salto Aleatório (RJT), Busca Aleatória Controlada (CRS) e Simplex Adaptativo de Nelder-Mead (ANMS). Desenvolvemos um algoritmo FWI-DFO que resolve numericamente a equação da onda acústica 2D pelo Método das Diferenças Finitas (FDM) e se utiliza de um método híbrido RJT-CRS-ANMS como técnica de otimização para a inversão sísmica. A estratégia é balancear automaticamente as buscas globais e locais iterativamente pelo CRS e ANMS, respectivamente. Aplicamos a metodologia a cinco modelos reais de subsuperfície. Os resultados mostraram uma ótima concordância com os modelos reais. O tempo computacional apresentou valores razoáveis e a função objetivo mostrou ser bastante sensível a pequenas alterações nos parâmetros do modelo para os casos aqui analisados. Em síntese, a metodologia FWI-DFO empregada se mostrou bastante promissora na inversão sísmica.Downloads
Referências
Alford, R., Kelly, K., Boore, D. (1974). Accuracy of finite-difference modeling of the acoustic wave equation. Society of Explorations Geophysicists, 39 (6), 834-842. doi:10.1190/1.1440470.
Almeida, V. T. X. D. (2013). Tomografia sísmica de superfície: modelagem direta calculando as trajetórias de ondas refratadas. (Dissertação de mestrado, Universidade Federal Fluminense). Recuperado de http://www.isis.uff.br/wp-content/uploads/2014/01/dissertacao_victor_almeida.pdf
Bakr, M. (2013). Nonlinear optimization in electrical engineering with applications in matlab. London: The Institution of Engineering and Technology.
Conn, A.R., Scheinberg, K., Vicente, L.N. (2009). Introduction to derivative-free optimization. Society for Industrial and Applied Mathematics. doi: 10.1137/1.9780898718768.
Cunha, M. C. C. (2009). Métodos numéricos. Campinas: UNICAMP.
Diniz-Ehrhardt, M.A. (2010). Métodos sem derivadas para minimização irrestrita. Notas em Matemática Aplicada, 49 (3), 87. Recuperado de https://www.sbmac.org.br/eventos/cnmac/xxxiii_cnmac/pdf/2049.pdf
Diniz-Ehrhardt, M.A., Martinez, J.M., Pedroso, L.G. (2011). Derivative-free methods for nonlinear programming with general lower-level constraints. Computational & Applied Mathematics, 30 (1), 19-52. Recuperado de http://www.scielo.br/scielo.php?pid=S1807- 03022011000100003&script=sci_arttext&tlng=es
Dourado, W. B., Imada, R. N., Pagamisse, A. (2013). Comparação entre a máscara de nitidez cúbica e o laplaciano para realce de imagens digitais. Anais do Congresso de Matemática Aplicada e Computacional Sudeste. Bauru/SP, Setembro. Recuperado de http://www.sbmac.org.br/cmacs/cmac-se/2013/trabalhos/
Fichtner, A. (2011). Full seismic waveform modelling and inversion. [Springer-Verlag].
doi: 10.1007/978-3-642-15807-0.
Gao, F. Han, L. (2012). Implementing the Nelder-Mead simplex algorithm with adaptive parameters. Computational Optimization and Applications. 51 (1), 259-277. doi: 10.1007/s10589-010-9329-3.
Lagarias, J. C., Reeds, J. A., Wright, M. H., Wright, P. E. (1998). Convergence properties of the Nelder-Mead simplex method in low dimensions. SIAM Journal on Optimization, 9(1), 112-147. doi: 10.1137/s1052623496303470.
Lima, M.V.A.G. (2006). Reflexão sísmica rasa: estimativa de parâmetros elásticos pela inversão dos coeficientes de reflexão da onda P acima do ângulo crítico. (Dissertação de mestrado -Universidade de São Paulo). Recuperado de ftp://www.grec.iag.usp.br/private/adm/katia/teses/agg/m_marcus_vag_lima.pdf
Manzanares-Filho, N., Sousa, B.S., Camacho, R.G.R., Albuquerque, R.B.F., Jorge, A.B. (2007). Airfoil shape optimization using a controlled random search algorithm. Anais do Congresso de THE INTERNATIONAL CONGRESS OF MECHANICAL ENGINEERING, Brasília, DF. Recuperado de http://www.abcm.org.br/app/webroot/anais/cobem/2007/pdf/COBEM2007-2028.pdf
Martins, M.A. (2015). Estratégias para inversão do campo de ondas completo associado a fontes sísmicas reais e virtuais. (Tese de doutorado - Universidade Federal do Rio de Janeiro). Recuperado de http://www.coc.ufrj.br/index.php?option=com_docman&view=document&alias=2696-martins-ma-tm-14-pdf&category_slug=2015-1&Itemid=428
Menke, W. (2012). Geophysical data analysis: discrete inverse theory. New York: Academic Press.
Mrinal, K. S. (2006). Seismic inversion. Society of Petroleum Engineers. Richardson, TX: Society of Petroleum Engineers.
Nelder, J. A., Mead, R. (1965). A Simplex Method for Function Minimization. The Computer Journal, 7 (4), 308-313. doi: 10.1093/comjnl/7.4.308.
Odile Gauthier, Jean Virieux, and Albert Tarantola (1986). Two?dimensional nonlinear inversion of seismic waveforms: Numerical results. GEOPHYSICS, 51(7), 1387-1403. doi: 10.1190/1.1442188
Pedroso, L. G. (2005) Sobre o desempenho de métodos de busca direta para minimização irrestrita. (Dissertação de mestrado – UNICAMP). Recuperado de http://repositorio.unicamp.br/handle/REPOSIP/306047
Price, W. L. (1977). A Controlled Random Search Procedure for Global Optimisation. The Computer Journal, 20 (4), 367-370. Recuperado de https://academic.oup.com/comjnl/article/20/4/367/393971/A-controlled-random-search-procedure-for-global
Reynolds, A. C. (1978). Boundary conditions for the numerical solution of wave propagation problems. Society of Explorations Geophysicists, 43 (6), 1099-1110. doi:10.1190/1.1440881.
Ricker, N. H. (1940). The form and nature of seismic waves and the structure of seismograms. Geophysics, 5 (4), 348-366. doi: 10.1190/1.1441816.
Rocha Jr., D. C. (2013). Inversão linearizada com regularização para obtenção de velocidades intervalares. (Monografia de Graduação - Universidade Federal da Bahia). Recuperado de http://www.cpgg.ufba.br/gr-geof/geo213/trabalhos-graducao/Daniel-Rocha.pdf
Tarantola, A. (2006). Inverse Problem Theory and Methods for Model Parameter Estimation. Paris: Society for Industrial and Applied Mathematics. doi: 10.1137/1.9780898717921.
Tran, K. T., Hiltunen, D. R. (2012). One-Dimensional Inversion of Full Waveforms using a Genetic Algorithm. Journal of Environmental & Engineering Geophysics, 17 (4), 197-213. doi: 10.2113/JEEG17.4.197.
Van Leeuwen, T. & Herrmann, F. J. (2013), Fast waveform inversion without source-encoding. Geophysical Prospecting, 61 (1), 10–19. doi:10.1111/j.1365-2478.2012.01096.x