Produção Científica
Artigo em Revista
Total variation regularization for depth-to-basement estimate: Part 2 — Physicogeologic meaning and comparisons with previous inversion methods We applied the mathematical basis of the total variation (TV) regularization to analyze the physicogeologic meaning of the TV method and compared it with previous gravity inversion methods (weighted smoothness and entropic Regularization) to estimate discontinuous basements. In the second part, we analyze the physicogeologic meaning of the TV method and compare it with previous gravity inversion methods (weighted smoothness and entropic regularization) to estimate discontinuous basements. Presenting a mathematical review of these methods, we show that minimizing the TV stabilizing function favors discontinuous solutions because a smooth solution, to honor the data, must oscillate, and the presence of these oscillations increases the value of the TV stabilizing function. These three methods are applied to synthetic data produced by a simulated 2D graben bordered by step faults. TV regularization and weighted smoothness are also applied to the real anomaly of Steptoe Valley, Nevada, U.S.A. In all applications, the three methods perform similarly. TV regularization, however, has the advantage, compared with weighted smoothness, of requiring no a priori information about the maximum depth of the basin. As compared with entropic regularization, TV regularization is much simpler to use because it requires, in general, the tuning of just one regularization parameter. ©2011 Society of Exploration Geophysicists |
Artigo em Revista
Partitioned least-squares operator for large-scale geophysical inversion Geophysics 75, R121, 2010 Least-squares (LS) problems are encountered in many geophysical estimation and data analysis problems where a large number of observations (data) are combined to determine a model (some aspect of the earth structure). Examples of least squares in seismic exploration include several data processing algorithms, theoretically accurate LS migration, inversion for reservoir parameters, and background velocity estimation. A frequently encountered problem is that the volume of data in 3D is so large that the matrices required for the LS solution cannot be stored within the memory of a single computer. A new technique is described for parallel computation of the LS operator that is based on a partitioned-matrix algorithm. The classical LS method for solution of block-Toeplitz systems of normal equation (NE) to the general case of block-Hermitian and non-Toeplitz systems of NE. is generalized. Specifically, a solution of a block-Hermitian system of NE is shown that may be obtained recursively by linearly combining the solutions of lesser order that are related to the forward and backward subsystems of equations. This results in an efficient parallel algorithm in which each partitioned system can be evaluated independently. The application of the algorithm to the problem of 3D plane wave transformation is demonstrated. ©2010 Society of Exploration Geophysicists |
Apresentação
Ground-roll attenuation based on SVD filtering We present a singular value decomposition (SVD) filtering method for attenuation of the ground roll. Before the SVD computation, the normal move-out (NMO) correction is applied to the seismograms, with the purpose of flattening the reflections. SVD is performed on a small number of traces in a sliding window. The output trace is the central trace of the first few eigenimages. These contains mostly horizontally aligned signals, and other noise in the data will be suppressed. The new method preserves the character and frequency content of the horizontal reflections and attenuates all other type of events. We illustrate the method using land seismic data of the Tacutu basin, located in the north-east part of Brazil. The results show that the proposed method is effective and is able to reveal reflections masked by the ground-roll. ©2009 Society of Exploration Geophysicists |
Artigo em Revista
Simultaneous 3D depth-to-basement and density-contrast estimates using gravity data and depth control at few points We have developed a gravity-inversion method for simultaneously estimating the 3D basement relief of a sedimentary basin and the parameters defining a presumed parabolic decay of the density contrast with depth in a sedimentary pack, assuming prior knowledge about the basement depth at a few points. The sedimentary pack is approximated by a grid of 3D vertical prisms juxtaposed in both horizontal directions of a right-handed coordinate system. The prisms' thicknesses represent the depths to the basement and are the parameters to be estimated from the gravity data. To estimate the parameters defining the parabolic decay of the density contrast with depth and to produce stable depth-to-basement estimates, we imposed smoothness on the basement depths and proximity between estimated and known depths at boreholes. We applied our method to synthetic data from a simulated complex 3D basement relief with two sedimentary sections having distinct parabolic laws describing the density-contrast variation with depth. The results provide good estimates of the true parameters of the parabolic law of density-contrast decay with depth and of the basement relief. Inverting the gravity data from the onshore and part of the shallow offshore Almada Basin on Brazil's northeastern coast shows good correlation with known structural features. |
Dissertação de Mestrado
Atenuação de reflexão múltipla do fundo marinho utilizando deconvolução preditiva adaptativa A atenuação das reflexões múltiplas de curto, médio ou longo perÃodos em dados sÃsmicos marÃtimos represnta um grande e importante problema do procesamento de dados sÃsmicos enfrentado pela indústria do petróleo.Tais reflexões múltiplas estão presntes nos dados sÃmicos adquirdos sobre lâmina de água e prejudicam sobremaneira oprocesamento dos dados e qualidafefinal da imagem sÃmica. O método de deconvolução preditiva de Wiener-Levinson (WL) visa a predição e subtração desses eventos periódicos.Os filtros preditvos são obtidos com o método dos mÃnimos quadrados e implicitamente utilzam a premissa usual de que os dados sÃsmicos possuem valores de amplitude nula fora da jqnela de tempo considerada. Desta forma os filtros de WL perdem eficácia na predição do evento periódico penalizando o resultado final da atenuação da múltipla. Outro fator limitante do sucesso da deconvolução preditva de WL é a ausência de periodicidade das reflexõesmúltiplas ao longo do traço sÃmico. A perfeita periodicidade ocorre apenas em registros de incidência normal, quando os traço são de afastamento nulo entre fonte-receptor, o que faz a deconvolução preditiva de WL nãoser tão eficaz. Na presente dissertação foi investigada a aplicação do método de deconvolução preditiva utilzando o algoritmo de Morf modificado (predição unitária L = 1 para predição arbitráia L > 1), que não faz considerações quanto à s amplitudes das amostras fora da janela especificada. Este algoritmo é apropriado para trabalhar no interior de janelas do traço sÃmicos sem que o filtro perca sua eficácia de predizer e, consequentemente, atenuar a reflexão múltipla. Utilzando o algoritmo de Morf modificado sobre janelas deslizantes de tempo foi definido um procedimento que denominou-se de deconvolução adaptativa de Morf modicada (DAMM). Este método foi comparado com resultados obtidos com o método de deconvolução adaptativa de WL(DAWL). Resultados com dados sintéticos e reais demonstrarama eficácia e superioridade da DAMM na atenuação de reflexões múltiplas. O emprego do método sobre dados reais da Bacia do Jequitinhonha mostraram que o método DAMM atenuou as múltiplas de maneira mais efetiva que o método DAWL. Os grandes problemas e dificuldades que são encontrados na atenuação das reflexões múltiplas de curto perÃodo na parte de lâmica d'água rasa da plataforma continental e as múltiplas de médio perÃodo na região da quebra da talude foram resolvidos satisfatoriamente com o método DAMM. Os resultados com a DAMM são considerados bons, inclusive na comparação com resultados obtidos com o método SRME, bastante utilizaod na indústria do petróleo. Os resultados promissores obtidos sugerem o emprego da DAMM no processamento sÃsmico bem como a extensão do método DAMM para a forma multicanal, esperando-se assim tornar o método ainda mais robusto e eficaz. |
Apresentação
An adaptive local-slope SVD filtering approach to enhance events on seismic sections We present an adaptive singular value decomposition (SVD) filtering method for enhancement of the spacial coherence of the reflections and for the attenuation of the uncorrelated noise. The SVD filtering is performed on a small number of traces and a small number of samples collected around each data component. The method uses the local slope of the reflections to re-sample the data set surrounding each data component and the SVD filtering is locally applied to compute the filtered data. |
Artigo em Revista
Métodos estocásticos para modelagem de escoamento estacionário e transiente em meios porosos. Uma maneira de incorporar a incerteza das medidas de campo e a variabilidade espacial nas propriedades hidráulicas de aquÃferos e estabelecer distribuições de probabilidade para os parâmetros fÃsicos do meio. As propriedades estatÃsticas do potencial hidráulico são então determinadas pela solução numérica das equações diferenciais estocásticas que regem o regime de escoamento no aquÃfero. Neste trabalho descrevemos a utilização de dois métodos de elementos finitos estocásticos (o método de Monte Carlo e o método da Colocação) para estimar a média e a variância do potencial hidráulico para fluxo saturado em meio estatisticamente heterogêneo, supondo que o coeficiente de transmissividade hidráulica é descrito por um processo lognormal. Um dos fatores decisivos na precisão numérica dos métodos é o comprimento de correlação associado à transmissividade. Discutimos também algumas configurações de baterias de extração que foram recentemente propostas para uma exploração adequada do aquÃfero Recôncavo na bacia do rio Capivara (Bahia, Brasil), comparando dois modelos clássicos e introduzindo a aleatoriedade da transmissividade em um dos arranjos de poços que foram propostos. |
Artigo em Revista
Influence of Sea Water Resistivity on MCSEM Data. O Marine Controlled Source ElectroMagnetic (MCSEM) é um método geofÃsico para a detecção de camadas resistivas contidas abaixo do assoalhooceânico. Neste trabalho nós modelamos os dados do MCSEM em um ambiente unidimensional, incluindo variações na resistividade da água do mar, na forma de camadas de resistividades uniformes. Estas variações na resistividade da água podem surgir devido à influência de correntes marinhas, gradientes de temperatura ou mudança na salinidade da água. Nós estudamos o efeito destas variações nos dados do método MCSEM. Nossos resultados mostram que a interpretação pode ser fortemente influenciada, principalmente quando são analisados os dados normalizados. Vimos que as mudanças na resistividade da água têm um efeito sobre os dados similar à quele de mudanças na profundidade da lâmina d'água, sendo que ambos influenciam na atenuação da chamada "air-wave". |
Artigo em Revista
Estudo de formas implementacionais da equação da onda imagem para remigração na profundidade. Neste trabalho estudamos teórica e numericamente novas formas implementacionais da equação da onda imagem para remigração em profundidade. Trata-se de uma equação diferencial parcial de segunda ordem semelhante à equação da onda acústica. Determinamos, além da consistência e da estabilidade, as condições de dispersão e dissipação de seis esquemas de diferenças finitas para esta equação em diversas formas obtidas por mudanças de variáveis. Estas condições não podem ser satisfeitas simultaneamente, ie., um resultado preciso não é de fácil obtenção. Testes numéricos confirmam os resultados teóricos da estabilidade para três esquemas, mas falham nos outros. Os esquemas avançados na variável da velocidade e avançados ou atrasados em profundidade foram os que tiveram os resultados mais proveitosos, pois estes esquemas satisfazem as previsões teóricas quanto à sua estabilidade e podem ser utilizados para realizar a propagação da onda imagem, sendo o primeiro para velocidades crescentes e o segundo para velocidades decrescentes. |
Artigo em Revista
Wave equation depth migration using complex Padé approximation. Propomos um novo método de migração em profundidade baseado na solução da equação da onda com densidade constante no domÃnio da freqüência. Uma aproximação de Padé complexa é usada para aproximar o operador de evolução aplicado na extrapolação do campo de ondas. Esse método reduz as imprecisões e instabilidades devido à s ondas evanescentes e produz imagens com menos ruÃdos numéricos que aquelas obtidas usando-se a aproximação de Padé real para o operador exponencial, principalmente em meios com fortes variações de velocidades. Testes em dados de afastamento nulo do modelo de sal SEG/EAGE e nos dados de tiro comum 2-D Marmousi foram realizados. Os resultados obtidos mostram que o método de migração proposto consegue lidar com fortes variações laterais e também tem uma boa resposta para refletores com mergulhos Ãngremes. Os resultados foram comparados à queles resultados obtidos com os métodos split-step Fourier (SSF), phase shift plus interpolarion (PSPI) e Fourier diferenças-finitas (FFD). |
<< < 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 > >>