Spelling suggestions: "subject:"inversão biofísica"" "subject:"cnversão biofísica""
1 |
Inversão de dados eletromagnéticos com o regularizador Variação Total e o uso da matriz de sensibilidade aproximadaLUZ, Edelson da Cruz 20 December 2012 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-22T12:11:30Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoDadosEletromagneticos.pdf: 3113358 bytes, checksum: 5be6334b23281e76cdc2db18d5c271cd (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar assuntos on 2014-08-06T16:08:31Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-28T16:30:54Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoDadosEletromagneticos.pdf: 3113358 bytes, checksum: 5be6334b23281e76cdc2db18d5c271cd (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T16:56:49Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoDadosEletromagneticos.pdf: 3113358 bytes, checksum: 5be6334b23281e76cdc2db18d5c271cd (MD5) / Made available in DSpace on 2014-09-10T16:56:49Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoDadosEletromagneticos.pdf: 3113358 bytes, checksum: 5be6334b23281e76cdc2db18d5c271cd (MD5)
Previous issue date: 2012 / ANP - Agência Nacional do Petróleo / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho apresentamos um estudo da aplicação do regularizador “Variação Total” (VT) na inversão de dados geofísicos eletromagnéticos. O regularizador VT reforça a proximidade entre os parâmetros adjacentes, mas, quando a influência de uma descontinuidade é sentida nos dados, este permite mudanças abruptas sobre os parâmetros. Isso faz com que o método seja uma alternativa válida, quando os dados observados usados na inversão provém de um ambiente geológico com uma distribuição suave de condutividade, mas que pode apresentar descontinuidades em lugares como as interfaces entre as camadas geoelétricas, como na margem de uma zona de óleo ou de um corpo de sal, que podem ser zonas muito resistivas no interior de sedimentos condutivos. Quando, devido a baixa resolução nos dados, o método não tem informações o suficiente para identificar a interface, o regularizador variação total reforça a proximidade entre os parâmetros adjacentes fazendo um transição suave entre as condutividades camadas, da mesma forma que é apresentado pela suavidade global. O método de Variação Total permite que modelos menos suaves sejam alcançados porque na norma L1 a medida de desajuste entre os pares de parâmetros adjacentes, dará o mesmo valor se a variação dos parâmetros é suave ou se a variação é abrupta, o que não é o caso se o mesmo desajuste é medido na norma L2, pois em uma distribuição suave a medida do desajuste é menor, sendo assim favorecida pela minimização desta norma. O uso deste regularizador permite uma melhor estimativa do tamanho de um corpo, seja ele resistivo ou condutivo. O trabalho está apresentado na forma de três artigos, cada um descrevendo uma etapa no desenvolvimento do problema da inversão, seguindo uma sequência de complexidade crescente no problema direto. O primeiro artigo neste trabalho é intitulado “Inversão de dados do CSEM marinho 1D de meio estratificado anisotrópico com o regularizador Variação Total”. Este descreve o passo inicial no desenvolvimento do problema: a inversão de dados do CSEM marinho de modelos estratificados 1D com anisotropia na condutividade das camadas. Este problema se presta bem para este desenvolvimento, porque tem solução computacional muito mais rápida do que o 2D, e nele já estão presentes as características principais dos dados do método CSEM marinho, como a largura muito grande da faixa de amplitudes medidas em um levantamento, e a baixa resolução, inerente às baixas frequências empregadas. A anisotropia acrescenta uma dificuldade a mais no problema, por aumentar o nível de ambiguidade nos dados e demandar ainda mais informação do que no caso puramente isotrópico. Os resultados mostram que a aplicação dos vínculos de igualdade do método VT permite a melhor identificação de uma camada alvo resistiva do que a simples aplicação dos vínculos tradicionais de suavidade. Até onde podemos aferir, esta solução se mostra superior a qualquer outra já publicada para este problema. Além de ter sido muito importante para o desenvolvimento de códigos em paralelo. O segundo artigo apresentado aqui, “Inversão de dados Magnetotelúricos com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada”, trata da inversão de dados do método Magnetotelúrico em ambientes 2D. Este problema demanda um esforço computacional muito maior do que o primeiro. Nele, estudamos a aplicação do método dos estados adjuntos para gerar uma boa aproximação para as derivadas necessárias para a construção da matriz de sensibilidade usada na inversão. A construção da matriz de sensibilidade é a etapa que demanda mais tempo no processo de inversão, e o uso do método de estados adjuntos foi capaz de reduzir muito este tempo, gerando derivadas com um bom nível de aproximação. Esta etapa da pesquisa foi fundamental pelo problema direto ser matematicamente e computacionalmente muito mais simples do que o do CSEM marinho 2D. Novamente em comparação com a aplicação do regularizador de suavidade global, o regularizador de Variação Total permitiu, neste problema, uma melhor delimitação das bordas de heterogeneidades bidimensionais. A terceira parte deste trabalho, apresentada no artigo “Inversão de dados do CSEM marinho 2.5D com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada”, apresenta a apliação do método de Variação Total ao problema da inversão de dados CSEM marinho 2.5D. Usamos o método dos estados adjuntos para gerar uma boa aproximação para as derivadas necessárias para a construção da matriz de sensibilidade usada na inversão, acelerando assim o processo de inversão. Para deixar o processo de inversão ainda mais rápido, lançamos mão da programação em paralelo com o uso de topologia. A comparação entre a aplicação do regularizador de suavidade global, e o regularizador de Variação Total permitiu, assim como nos casos anteriores, uma melhor delimitação das bordas de heterogeneidades bidimensionais.
|
2 |
Análise de velocidade por otimização do semblance na reflexão sísmicaVIEIRA, Wildney Wallacy da Silva January 2010 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-06T12:37:09Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseVelocidadeOtimizacao.pdf: 2382926 bytes, checksum: 413bbca447f7bac1271afc2aae31c08d (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T12:35:46Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-07T12:54:14Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseVelocidadeOtimizacao.pdf: 2382926 bytes, checksum: 413bbca447f7bac1271afc2aae31c08d (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-18T14:09:19Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseVelocidadeOtimizacao.pdf: 2382926 bytes, checksum: 413bbca447f7bac1271afc2aae31c08d (MD5) / Made available in DSpace on 2014-08-18T14:43:01Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseVelocidadeOtimizacao.pdf: 2382926 bytes, checksum: 413bbca447f7bac1271afc2aae31c08d (MD5)
Previous issue date: 2011 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / FINEP - Financiadora de Estudos e Projetos / PETROBRAS - Petróleo Brasileiro S.A. / Este trabalho teve como objetivo geral desenvolver uma metodologia sistemática para a inversão de dados de reflexão sísmica em arranjo ponto-médio-comum (PMC), partindo do caso 1D de variação vertical de velocidade e espessura que permite a obtenção de modelos de velocidades intervalares, vint,n, as espessuras intervalares, zn, e as velocidades média-quadrática, vRMS,n, em seções PMC individualizadas. Uma consequência disso é a transformação direta destes valores do tempo para profundidade. Como contribuição a análise de velocidade, foram desenvolvidos dois métodos para atacar o problema baseado na estimativa de velocidade intervalar. O primeiro método foi baseado na marcação manual em seções PMC, e inversão por ajuste de curvas no sentido dos quadrados-mínimos. O segundo método foi baseado na otimização da função semblance para se obter uma marcação automática. A metodologia combinou dois tipos de otimização: um Método Global (Método Price ou Simplex), e um Método Local (Gradiente de Segunda Ordem ou Conjugado), submetidos a informação à priori e vínculos. A marcação de eventos na seção tempo-distância faz parte dos processos de inversão, e os pontos marcados constituem os dados de entrada juntamente com as informações à priori do modelo a ser ajustado. A marcação deve, por princípio, evitar eventos que representem múltiplas, difrações e interseções, e numa seção pode ser feita mais de 50 marcações de eventos, enquanto que num mapa semblance não se consegue marcar mais de 10 eventos de reflexão. A aplicação deste trabalho é voltada a dados sísmicos de bacias sedimentares em ambientes marinhos para se obter uma distribuição de velocidades para a subsuperfície, onde o modelo plano-horizontal é aplicado em seções PMC individualizadas, e cuja solução pode ser usada como um modelo inicial em processos posteriores. Os dados reais da Bacia Marinha usados neste trabalho foram levantados pela PETROBRAS em 1985, e a linha sísmica selecionada foi a de número L5519 da Bacia do Camamu, e o PMC apresentado é a de número 237. A linha é composta de 1098 pontos de tiro, com arranjo unilateraldireito. O intervalo de amostragem é 4 ms. O espaçamento entre os geofones é 13,34 m com o primeiro geofone localizado a 300 m da fonte. O espaçamento entre as fontes é de 26,68 m. Como conclusão geral, o método de estimativa de velocidade intervalar apresentada neste trabalho fica como suporte alternativo ao processo de análise de velocidades, onde se faz necessário um controle sobre a sequência de inversão dos PMCs ao longo da linha sísmica para que a solução possa ser usada como modelo inicial ao imageamento, e posterior inversão tomográfica. Como etapas futuras, podemos propor trabalhos voltados direto e especificamente a análise de velocidade sísmica estendendo o caso 2D de otimização do semblance ao caso 3D, estender o presente estudo para o caso baseado na teoria do raio imagem com a finalidade de produzir um mapa continuo de velocidades para toda a seção sísmica de forma automática. / This work had as a general aim to develop a systematic methodology for the inversion of seismic reflection data organized in common-midpoint gathers (CMP), starting from 1D vertical variation of velocity and thickness that allows to obtain interval velocity, vint,n, in time, the correspondent interval thickness, zn, and the correspondent mean-square velocity, vRMS,n in individualized CMP gathers. A direct consequence of this work the transformation of these values from time to depth. Two methods were developed to attack the problem defined as velocity analysis based on the estimation of interval velocity. The first method was based on manual picking of reflection events on CMP gathers, and inversion by curve fitting in the least-square sense. The second method was based on the otimization of the semblance function to obtain an automatic picking. The methodology combined two types of optimization: a Global Method (Price or Simplex), and Local Method (second order gradient or cojugate), subject to a priori information and constraints. The picking of events in time-distance section is of fundamental importance in the process of inversion, and the picked points are the input data along with a priori information of the model to be adjusted. The picking must, in principle, avoid events that represent multiples, diffractions and intersections, and in a section over 50 pickings can be made, while in a semblance map not more than 10 events could usually be picked by eye. The application of this work is focused on seismic data of marine sedimentary basins to obtain a distribution of velocities for the subsurface, where a plane-horizontal model is applied for individual CMP sections, and that the solution can be used as an initial model in subsequent processes. The real data used in this study were collected by Petrobras in 1985, and the selected seismic line was of number L5519 of the Camamu Basin, and the CMP presented is of number 237. The line consists of 1098 shot points with right-lateral arrangement. The sampling interval is 4 ms. The spacing between the geophones is 13.34 m with the first geophone located at 300 m from the source. The spacing between the sources is 26.68 m. As a general conclusion, the method for estimating interval velocity in this work stands as an alternative support to velocity analysis, where it is necessary a control over the sequential inversion of CMP gathers along the seismic line such that the solution can be used as an initial model for imaging, and further tomographic inversion. As future work, we can be propose studies directely and specifically related to seismic velocity analysis by extending the 2D semblance optimization method to 3D, and extending the present studies to the method based on the image ray, aiming at producing a continuous velocity map for the entire section in an automatic way.
|
3 |
Inversão de momentos de fonte em métodos potenciaisMEDEIROS, Walter Eugênio de 16 August 1993 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-28T12:03:20Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoMomentosFonte.pdf: 26856422 bytes, checksum: 4c2d4fd4c3ccfc0072a0e1986811b2a7 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:14:12Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-19T12:49:34Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoMomentosFonte.pdf: 26856422 bytes, checksum: 4c2d4fd4c3ccfc0072a0e1986811b2a7 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T11:56:40Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoMomentosFonte.pdf: 26856422 bytes, checksum: 4c2d4fd4c3ccfc0072a0e1986811b2a7 (MD5) / Made available in DSpace on 2014-08-20T11:56:41Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoMomentosFonte.pdf: 26856422 bytes, checksum: 4c2d4fd4c3ccfc0072a0e1986811b2a7 (MD5)
Previous issue date: 1993 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / UFPA - Universidade Federal do Pará / A inversão de momentos de fonte gravimétrica tridimensional é analisada em duas situações. Na primeira se admite conhecer apenas a anomalia. Na segunda se admite conhecer, além da anomalia, informação a priori sobre o corpo anômalo. Sem usar informação a priori, mostramos que é possível determinar univocamente todo momento, ou combinação linear de momentos, cujo núcleo polinomial seja função apenas das coordenadas Cartesianas que definem o plano de medida e que tenha Laplaciano nulo. Além disso, mostramos que nenhum momento cujo núcleo polinomial tenha Laplaciano não nulo pode ser determinado. Por outro lado, informação a priori é implicitamente introduzida se o método de inversão de momentos se baseia na aproximação da anomalia pela série truncada obtida de sua expansão em multipolos. Dado um centro de expansão qualquer, o truncamento da série impõe uma condição de regularização sobre as superfícies equipotenciais do corpo anômalo, que permite estimar univocamente os momentos e combinações lineares de momentos que são os coeficientes das funções-bases da expansão em multipolos. Assim, uma distribuição de massa equivalente à real é postulada, sendo o critério de equivalência especificado pela condição de ajuste entre os campos observado e calculado com a série truncada em momentos de uma ordem máxima pré-estabelecida. Os momentos da distribuição equivalente de massa foram identificados como a solução estacionária de um sistema de equações diferenciais lineares de 1a. ordem, para a qual se asseguram unicidade e estabilidade assintótica. Para a série retendo momentos até 2a. ordem, é implicitamente admitido que o corpo anômalo seja convexo e tenha volume finito, que ele esteja suficientemente distante do plano de medida e que a sua distribuição espacial de massa apresente três planos ortogonais de simetria. O método de inversão de momentos baseado na série truncada (IMT) é adaptado para o caso magnético. Para este caso, mostramos que, para assegurar unicidade e estabilidade assintótica, é suficiente pressupor, além da condição de regularização, a condição de que a magnetização total tenha direção e sentido constantes, embora desconhecidos. O método IMT baseado na série de 2a. ordem (IMT2) é aplicado a anomalias gravimétricas e magnéticas tridimensionais sintéticas. Mostramos que se a fonte satisfaz as condições exigidas, boas estimativas da sua massa ou vetor momento de dipolo anômalo total, da posição de seu centro de massa ou de momento de dipolo e das direções de seus três eixos principais são obtidas de maneira estável. O método IMT2 pode falhar parcialmente quando a fonte está próxima do plano de medida ou quando a anomalia tem efeitos localizados e fortes de um corpo pequeno e raso e se tenta estimar os parâmetros de um corpo grande e profundo. Definimos por falha parcial a situação em que algumas das estimativas obtidas podem não ser boas aproximações dos valores verdadeiros. Nas duas situações acima descritas, a profundidade do centro da fonte (maior) e as direções de seus eixos principais podem ser erroneamente estimadas, embora que a massa ou vetor momento de dipolo anômalo total e a projeção do centro desta fonte no plano de medida ainda sejam bem estimados. Se a direção de magnetização total não for constante, o método IMT2 pode fornecer estimativas erradas das direções dos eixos principais (mesmo se a fonte estiver distante do plano de medida), embora que os demais parâmetros sejam bem estimados. O método IMT2 pode falhar completamente se a fonte não tiver volume finito. Definimos por falha completa a situação em que qualquer estimativa obtida pode não ser boa aproximação do valor verdadeiro. O método IMT2 é aplicado a dados reais gravimétricos e magnéticos. No caso gravimétrico, utilizamos uma anomalia situada no estado da Bahia, que se supõe ser causada por um batólito de granito. Com base nos resultados, sugerimos que as massas graníticas geradoras desta anomalia tenham sido estiradas na direção NNW e adelgaçadas na direção vertical durante o evento compressivo que causou a orogênese do Sistema de Dobramentos do Espinhaço. Além disso, estimamos que a profundidade do centro de massa da fonte geradora é cerca de 20 km. No caso magnético, utilizamos a anomalia de um monte submarino situado no Golfo da Guiné. Com base nos resultados, estimamos que o paleopolo magnético do monte submarino tem latitude 50°48'S e longitude 74°54'E e sugerimos que não exista contraste de magnetização expressivo abaixo da base do monte submarino. / The inversion of three-dimensional gravity source moments is analyzed in two situations. In the first one only the anomalous field is assumed to be known. In the second situation a priori information about the anomalous body is assumed to be known besides the field data. Without using a priori information, we show that it is possible to determine uniquely any moment, or linear combination of moments, whose polynomial kernel: (a) is not a function of the Cartesian coordinate which is orthogonal to the measuring plane and (b) has null Laplacian. Besides, we show that it is impossible to determine any moment whose polynomial kernel has non-null Laplacian. On the other hand, we show that a priori information is implicitly introduced if the source moment inversion method is based on the approximation of the anomalous field by the truncated series obtained from its multipole expansion. Given any center of expansion, the series truncation impores a regularization condition on the equipotential surfaces of the anomalous body that allows the moments and linear combination of moments (which are the coefficients of the multipole expansion basis function) to be uniquely estimated. So, a mass distribution equivalent to the real mass distribution is postulated, being the equivalence criterion specified by the fitting conditions between the observed anomaly and the anomaly calculated with the truncated multipole expansion series. The highest order for the retained terms in the truncated series is specified by the previously defined maximum order for the moments. The moments of the equivalent mass distribution were identified as the stationary solution of a system of first order linear differential equations, for which uniqueness and assymptotic stability are assured. For the series having moments up to 2nd order, it is implicitly assumed that the anomalous body: (1) has finite volume, (2) that it is sufficiently far from the measuring plane and (3) that its spatial naass distribution is convex and presents three orthogonal planes of symmetry. The source moment inversion method based on the approximation of the anomalous field by a truncated series (MIT) is adapted to the magnetic case. In this case, we show that in order to guarantee uniqueness and assymptotic stability it is sufficient to assume, besides the regularization condition, that the total magnetization has constant but unknown direction. The MIT method based on the 2nd order series (MIT2) is applied to three-dimensional synthetic gravity and magnetic anomalies. If the source satisfies all imposed conditions, we show that it is possible to obtain in a stable way good estimates of the total anomalous mass or dipole moment vector, of the position of center of mass or dipole moment and of the directions of all three principal axes. A partia' failure of MIT2 method may occur either if the source is dose to the measuring plane or if the anomaly presents a localized but strong effect due to a shallow and small body and an attempt is made to estimate the moments of a large and deep body. By partial failure we mean the situation when some of the estimates may be poor aproximations of the true values. In these two cases we show that the estimates of the depth and the directions of the principal axes of the (main) source may be poor but the estimates of the total anomalous mass or dipole moment vector and the projection on the measuring plane of the center of mass or dipole moment of the source are good. If the total magnetization direction is not constant, MIT2 method may produce poor estimates of the directions of the principal axes (even if the source is far from the measuring plane) but good estimates are obtained for the other parameters. A complete failure of MIT2 method may occur if the source does not have finite volume. By complete failure we mean the situation when any obtained estimate may be a poor aproximation of the true value. MIT2 method is applied to real gravity and magnetic data. In the gravimetric case we used an anomaly located in Bahia state, Brazil, which is assumed to be produced by the presence of a large granitic body. Based on the inversion results, we propose that the grafite was deformed into an oblate ellipsoid during the compressive event that generated the Middle Proterozoic Espinhaço orogeny. The center of mass estimated for this body is about 20 km. In the magnetic case, we used an anomaly produced by a seamount located in the Gulf of Guinea. Based on the inversion results, we estimate a magnetic palaeopole for the seamount at 50°48'S and 74°54'E and we suggest that no important magnetization contrast exists below the bottom of the seamount.
|
4 |
Inversão vinculada de dados de EMAP e MT-1D anisotrópicoRÉGIS, Cícero Roberto Teixeira 26 November 1999 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-22T14:44:50Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoVinculadaDados.pdf: 34199418 bytes, checksum: fccd8f42ebfa3799899d9531b88a0cb8 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:13:16Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-29T17:16:12Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoVinculadaDados.pdf: 34199418 bytes, checksum: fccd8f42ebfa3799899d9531b88a0cb8 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T12:55:01Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoVinculadaDados.pdf: 34199418 bytes, checksum: fccd8f42ebfa3799899d9531b88a0cb8 (MD5) / Made available in DSpace on 2014-09-10T12:55:01Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoVinculadaDados.pdf: 34199418 bytes, checksum: fccd8f42ebfa3799899d9531b88a0cb8 (MD5)
Previous issue date: 2000 / Uma técnica para a inversão de dados magnetotelúricos é apresentada neste trabalho. Dois tipos de dados são tratados aqui, dados gerados por modelos unidimensionais com anisotropia na condutividade das camadas e dados bi-dimensionais de levantamentos do método EMAP (ElectroMagnetic Array Profiling). Em ambos os casos fazemos a inversão usando vínculos aproximados de igualdade para estabilizar as soluções. Mostramos as vantagens e as limitações do uso destes vínculos nos processos de inversão. Mesmo vinculada a inversão ainda pode se tornar instável. Para inverter os dados 2-D do EMAP, apresentamos um processo que consiste de três partes: 1 – A construção de um modelo interpretativo e da aproximação inicial para a inversão a partir dos dados de seções de resistividade aparente filtradas pelo processo de filtragem do EMAP; 2 – a inclusão de uma camada de corpos pequenos aflorantes, chamada de camada destatic shift, aos modelos interpretativos para resolver as fontes de distorções estáticas que contaminam os dados; 3 – o uso dos vínculos aproximados de igualdade absoluta para estabilizar as soluções. Os dois primeiros passos nos permitem extrair o máximo de informação possível dos dados, enquanto que o uso dos vínculos de igualdade nos permite incluir informação a priori que possua significado físico e geológico. Com estes passos, obtemos uma solução estável e significativa. Estudaremos o método em dados sintéticos de modelos bi-dimensionais e em dados reais de uma linha de EMAP feita na Bacia do Paraná. / We present a technique to invert magnetotelluric data. Two kinds of MT data are treated: data generated by anisotropic one-dimensional models and two-dimensional data from EMAP (ElectroMagnetic Array Profiling) surveys. In both cases we perform the inversions using approximate equality constraints to stabilize the solutions. We show the advantages and the limitations of those constraints in the inversion process. Even when the inversion is constrained, it can still be unstable. To invert the two-dimensional EMAP data we present a process that consists of three steps: 1 – The interpretation model and the initial guess for each case are derived from the filtered apparent resistivity sections which are the result of the filtering process of the EMAP method; 2 – One layer formed by small outcropping bodies (termed static shift layer) is added to the interpretation model in order to resolve the sources of static distortions that contaminate the data; 3 – The inversion is constrained by approximate equality constraints. The first two steps enable us to extract the largest possible amount of information from the data, while the use of the approximate equality constraints allow us to incorporate a priori information which has true physical or geological meaning. With these steps we obtain solutions that are stable and meaningful. We study the method in two-dimensional models as well as in real data from an EMAP line surveyed at the Paraná Basin.
|
5 |
Inversão de dados de ondas de superfície RayleighSANTA ROSA, Antonio Nuno de Castro 26 February 1996 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-22T13:48:05Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoDadosOndas.pdf: 140994516 bytes, checksum: 401ef5e70ef9a780409ce40a7c36d018 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:22:02Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-01T16:35:48Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoDadosOndas.pdf: 140994516 bytes, checksum: 401ef5e70ef9a780409ce40a7c36d018 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T12:55:18Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoDadosOndas.pdf: 140994516 bytes, checksum: 401ef5e70ef9a780409ce40a7c36d018 (MD5) / Made available in DSpace on 2014-09-10T12:55:18Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Tese_InversaoDadosOndas.pdf: 140994516 bytes, checksum: 401ef5e70ef9a780409ce40a7c36d018 (MD5)
Previous issue date: 1996 / Neste trabalho compilamos informações sobre um grande número de medidas de velocidade de grupo para ondas Rayleigh do modo fundamental, com período até 100 segundos. Tais dados consistiram de informações retiradas da literatura geofísica e cobriram toda a Terra. Parte dos dados foi organizada em trabalhos anteriores e uma segunda parte foi apresentada aqui de forma inédita. Para a América do Sul, selecionamos os principais conjuntos de dados de tais ondas e elaboramos diversos perfis onde a distribuição de velocidade de ondas cisalhantes foi obtida a partir da inversão das curvas de dispersão de velocidade de grupo. Tais perfis serviram para termos uma ideia inicial da estrutura interna da Terra em nosso continente. Com o conjunto global de dados de velocidade de grupo foi possível obtermos os mapas de distribuição lateral de valores de velocidade para cada período referencial entre 20 e 100 segundos. Tais mapas foram produzidos da mesma forma que os mapas de velocidade de fase de ROSA (1986), onde a amostragem for para realizada para blocas medindo 10x10 graus, englobando toda a Terra, em projeção mercator. O valor de velocidade de grupo em cada bloco, para cada período, foi obtido a partir da inversão estocástica dos dados de anomalia de velocidade em relação aos modelos regionalizados de JORDAN (1981) com os valores de velocidade de grupo de ROSA et al. (1992). Os mapas de velocidade de grupo obtidos aqui foram então empregados, na América do Sul, com os valores de velocidade de fase dos mapas obtidos por ROSA (1986). Assim, foi possível determinarmos, em profundidade, os mapas de variação de velocidade de onda cisalhante e os mapas de distribuição de valores de densidade. Com isto, pudemos construir o primeiro mapa de profundidade do Moho (todo do Manto Superior) da América do Sul. / We have compiled a large set of group velocity measurements for fundamental mode Rayleigh waves with period smaller than 100 seconds. These data are composed of worldwide information gathered from the geophysical literature. Part of the data was compiled in previous works, and a second part of these data was measured in this thesis. For the South American plate, we have selected the main sets of measurements for such waves and gathered them in several profiles, for which the depth distribution of the shear wave velocity was determined from the inversion of the group velocity dispersion curves. These depth profiles were useful to have an idea of the internal structure of the Earth underneath the South American plate. Using the global group velocity data set, it was possible to determine the maps showing the lateral variation of group velocity for reference period values ranging from 20 to 100 seconds. These maps were produced in the same way that ROSA (1986) used to generate the corresponding phase velocity maps, using blocks measuring 10x10 degrees around the earth and considering the mercator projection. The group velocity value at each block was obtained, for each reference period value, from the stochastic inversion of the travel time delays from the JORDAN (1981) each model, considering the average group velocity values determined by ROSA et al. (1992) for this model. The group velocity maps obtained here were then jointly used with the phase velocity values determined by ROSA (1986) for the determination of the earth’s internal distribution of shear wave velocity as well as the density distribution, from the depth inversion of both phase and group velocity dispersion curves. The results enabled us to construct the first map showing the possible depth of the Moho discontinuity (upper mantle depth) for South America.
|
6 |
Análise de um esquema linear para estimativa de anisotropia local a partir de dados de onda P em experimentos de VSP multiazimutal / Analysis of a linear scheme for estimation of local anisotropy from P-wave data in multiazimuth VSP surveysMACAMBIRA, Raiza de Nazaré Assunção 08 August 2013 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T14:17:45Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseEsquemaLinear.pdf: 958488 bytes, checksum: d78a32e1815ce14cc9e903e1c610c694 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Incluir palavras-chave on 2014-08-07T14:17:45Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-11T12:42:20Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseEsquemaLinear.pdf: 958488 bytes, checksum: d78a32e1815ce14cc9e903e1c610c694 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-15T16:01:43Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseEsquemaLinear.pdf: 958488 bytes, checksum: d78a32e1815ce14cc9e903e1c610c694 (MD5) / Made available in DSpace on 2014-09-15T16:01:43Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseEsquemaLinear.pdf: 958488 bytes, checksum: d78a32e1815ce14cc9e903e1c610c694 (MD5)
Previous issue date: 2013 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Neste trabalho é apresentada uma análise do esquema de inversão linear para a
estimativa de anisotropia na vizinhança de um receptor situado em um poço a partir de
da componente vertical do vetor de vagarosidade e do vetor de polarização de ondas P
medidops em experimentos de VSP walkaway multiazimutal. Independente do meio
acima do geofone (homogêneo ou heterogêneo) e da forma do poço (pode ser direcional
ou curvado, vertical e inclinado), a inversão é feita a partir de uma aproximação de
primeira ordem em torno de um meio isotrópico de referência. O esquma da inversão é
analisado considerando fatores como: o nível de ruído nos dados, o tipo de onda P, o
grau de anisotropia do meio, a escolha dos parâmetros no meio isotrópico de referência
e grau de heterogeneidade do meio. Os resultados são apresentados. / This study presents an analysis of the linear inversion scheme for estimating anisotropy
in the neighborhood of a receiver placed in a well using the vertical components of the
slowness and polarization vectors of P-waves measured in multi-azimuth walkaway
vertical seismic profile (VSP) surveys. Independently of the medium above the
geophone (homogeneous or heterogeneous) and the shape of the well (directional or
curved, vertical or sloped), an inversion is performed from a first-order approximation
around a reference isotropic medium. The analysis of the inversion scheme considers
factors such as the noise level of the data, the type of P-wave, the degree of the
anisotropy of the medium, the choice of parameters in the reference isotropic medium
and the degree of heterogeneity of the medium.
|
7 |
Caracterização de reservatórios fraturados através de dados de ondas qP em levantamentos VSP WalkawaySILVA, Saulo da Costa e January 2008 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-12T11:26:19Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_CaracterizacaoReservatoriosFraturados.pdf: 1892695 bytes, checksum: c9b2c5a8a38811f366576b899c2b9603 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:23:47Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-04T13:34:51Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_CaracterizacaoReservatoriosFraturados.pdf: 1892695 bytes, checksum: c9b2c5a8a38811f366576b899c2b9603 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-15T16:04:02Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_CaracterizacaoReservatoriosFraturados.pdf: 1892695 bytes, checksum: c9b2c5a8a38811f366576b899c2b9603 (MD5) / Made available in DSpace on 2014-09-15T16:04:02Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_CaracterizacaoReservatoriosFraturados.pdf: 1892695 bytes, checksum: c9b2c5a8a38811f366576b899c2b9603 (MD5)
Previous issue date: 2008 / ANP - Agência Nacional do Petróleo / Esta dissertação apresenta um método para estimar a orientação de fratura de um meio partindo do pressuposto que este se comporta efetivamente como um meio transversalmente isotrópico (TI). Desta forma, a orientação de fratura é obtida a partir da estimativa do eixo de simetria de um meio TI. Esta estimativa é obtida a partir dos vetores de vagarosidade e polarização de ondas qP medidos em experimentos VSP Walkaway. O processo de inversão é baseado em equações linearizadas
dos vetores de vagarosidade e polarização de ondas qP e os parâmetros fracamente anisotrópicos do meio. São apresentados testes numéricos, em que a sensibilidade a fatores como grau de anisotropia, geometria do levantamento, tipo de onda utilizado e nível de ruído são analisados. São também apresentados os resultados para um conjunto de dados reais. / This dissertation presents a method to estimate the orientation of the axis of symmetry of a medium assuming it behaves effectively as a transversely isotropic medium (TI). The fracture orientation is
then achieved from the estimation of the axis of symmetry of a TI medium. This estimative is done with the slowness and polarization vectors from qP waves, measured in VSP Walkaway experiments. The inversion process is based on linearized equations of the slowness and polarization vectors from qP waves and the weakly anisotropic parameters of the medium. Numerical tests are presented, on which the sensibility to factors like the strength of anisotropy, survey geometry, type
of wave utilized and noise level are analyzed. Test results for a set of real data are also shown.
|
8 |
Análise das aproximações R<sub>PP</sub> e R<sub>SP</sub> para meios isotrópicosSANTOS, Darcicléa Ferreira 16 August 2002 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T12:13:44Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:07:51Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-05T13:52:00Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-22T11:21:42Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5) / Made available in DSpace on 2014-09-22T11:21:42Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_AnaliseAproximacoesRpp.pdf: 4103377 bytes, checksum: 5e03cda58fdd1cf4609469596ddec5ad (MD5)
Previous issue date: 2002 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / ANP - Agência Nacional do Petróleo / FINEP - Financiadora de Estudos e Projetos / Neste trabalho são apresentadas aproximações lineares e quadráticas das equações de Zoeppritz para a obtenção dos coeficientes de reflexão e transmissão de eventos P-P e P-S em função dos ângulos de incidência e da média angular, bem como a análise de inversão linear de AVO, considerando os eventos de reflexão P-P e P-S dissociados e combinados. O uso das
chamadas aproximações pseudo-quadráticas foi aplicado para obtenção de aproximações quadráticas apenas para eventos-PP, em torno dos contrastes médios das velocidades de ondas compressionais e cisalhantes e da razão Vs/Vp. Os resultados das aproximações desenvolvidas neste trabalho mostram que as aproximações quadráticas são mais precisas que as lineares, nas duas versões angulares. As comparações entre as aproximações em termos do ângulo de
incidência e da média angular mostram que as aproximações quadráticas são equivalentes dentro
do limite angular de [0º, 30º]. Por outro lado, as aproximações lineares em função do ângulo de
incidência mostram-se mais precisas que as aproximações lineares em função da média angular.
Na inversão linear, fez-se análises de sensibilidade e de ambigüidade e observou-se que, nos caso
de eventos de reflexão P-P e P-S dissociados, apenas um parâmetro pode ser estimado e que a
combinação destes eventos consegue estabilizar a inversão, permitindo a estimativa de dois dentre os parâmetros físicos dos meios (contrastes de impedância, de velocidade de onda P e de módulo de cisalhamento). / This work presents linear and quadratic approximations of the Zoeppritz equations for the
derivation of reflection and transmission coefficients from P-P and P-S events as a function of incidence angles and angular average, as well as the linear inversion analysis, AVO, in respect to the disassociated and combined P-P and P-S reflection events. The use of the so-called pseudoquadratic approximations was applied for the derivation of quadratic approximations only for PPevents,
around the average contrasts of compressional and shear waves velocities and Vs/Vp ratio. The results of approximations derived in this work show that quadratic approximations were more precise than the linear ones in the two angular versions. The comparisons between
these approximations in terms incidence angle and angular average show that the quadratic
approximations are equivalent within the angular limit of [0º to 30º]. In the other hand, the linear
approximations as a function of incidence angle are more precise than the linear approximations
as a function of the angular average. In the linear inversion, sensitivity and ambiguity analyses
were carried out and one could see that in the case of disassociated reflection P-P and P-S events,
just a parameter can be estimated and the combination of these events can stabilize the inversion
permitting the estimation of two physical parameters for the media (impedance, P-wave velocity
and shear bulk module contrasts).
|
9 |
Inversão da forma de onda orientada ao alvoCOSTA, Carlos Alexandre Nascimento da 16 September 2016 (has links)
Submitted by Cássio da Cruz Nogueira (cassionogueirakk@gmail.com) on 2017-08-10T15:47:06Z
No. of bitstreams: 2
license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5)
Tese_InversaoFormaOnda.pdf: 6710248 bytes, checksum: 11d58ee38a890375b5e4f94c3fd3e8d9 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2017-08-23T13:46:27Z (GMT) No. of bitstreams: 2
license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5)
Tese_InversaoFormaOnda.pdf: 6710248 bytes, checksum: 11d58ee38a890375b5e4f94c3fd3e8d9 (MD5) / Made available in DSpace on 2017-08-23T13:46:27Z (GMT). No. of bitstreams: 2
license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5)
Tese_InversaoFormaOnda.pdf: 6710248 bytes, checksum: 11d58ee38a890375b5e4f94c3fd3e8d9 (MD5)
Previous issue date: 2016-09-16 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FAPESP - Fundação de Amparo à Pesquisa do Estado de São Paulo / INCT/GP - Instituto Nacional de Ciência e Tecnologia de Geofísica do Petróleo / PETROBRAS - Petróleo Brasileiro S.A. / Propomos uma nova metodologia de inversão da forma de onda orientada ao alvo para estimar os parâmetros físicos de uma área alvo em subsuperfície para dados sísmicos adquiridos com aquisição VSP-desviado ou com aquisição com fontes e receptores localizados na superfície. Além disso, investigamos a importância de eventos de múltiplos espalhamentos no conjunto de dados usados como inputs para estimar as respostas ao impulso da área alvo em subsuperfície iterativamente através de um esquema de inversão esparso para as mesma geometria de aquisição anteriormente citadas. Essas metodologias são baseadas no ajuste entre os campos de onda ascendente observado e modelado em um específico nível em profundidade próximo da área alvo, onde o campo de onda ascendente modelado é estimado através da representação tipo-convolução para função de Green. A principal característica de nossa metodologia de inversão da forma de onda orientada ao alvo é usar como inputs os campos de onda descendente e ascendente para estimar os parâmetros físicos locais sem necessitar que estes campos de onda sejam locais, dessa maneira evitamos adotar qualquer metodologia de redatumação para estimar as respostas ao impulso locais. Mostramos através de exemplos numéricos que podemos relaxar a necessidade de uma densa amostragem de fontes e receptores, característicos de esquemas de redatumação baseados em interferometria, para estimar as respostas ao impulso da área alvo através do esquema de inversão esparso usando como inputs os campos de onda com eventos de múltiplos espalhamentos. Estas metodologia são atraentes para dados sísmicos de aquisição VSP-desviado, pois para este tipo de dado não é necessário conhecer a região localizada acima da área alvo para estimar os campos de onda ascendente e descendente usados como inputs em ambas as metodologias. No entanto, para dados sísmicos adquiridos com fontes e receptores localizados na superfície, em princípio é necessário conhecer um modelo de velocidade a priori para estimar os campos de onda ascendente e descendente próximo da área alvo. Para ambas as metodologias investigadas, usamos a inversão da forma de onda baseada em migração para estimar os campos de onda ascendente e descendente próximo da área alvo a partir de dados sísmicos adquiridos com fontes e receptores na superfície. / We propose a new target-oriented waveform inversion to estimate the physical parameters from a specific target in the subsurface from observed data from deviated-VSP acquisition or surface seismic data. Furthermore, we investigate a strategy to estimate the impulse responses from a local target in the subsurface from deviated-VSP acquisition or surface seismic data as an iterative sparse inversion approach, where the main feature of this strategy is that all multiple scattering in the data is used to enhance the illumination at target level. In these approaches we fit the upgoing wavefields observed at a specific level near the local target with the upgoing wavefields estimated at same depth level through convolution-type representation for the Green’s function. The main feature of the target-oriented waveform inversion is that we just need to know the up- and downgoing wavefields at the depth level above the target area to estimate the physical parameters for the area of interest. We show through numerical tests that the iterative sparse inversion approach does not require dense sources sampling to estimate the impulse responses from a target below a complex overburden, because of all the extra illumination via multiples. The physical parameters above the target area is not necessary to know if we use the data from deviated-VSP geometry of acquisition, but for surface seismic data we need to know a smooth physical parameter above the target area to estimate the up- and downgoing wavefields at depth level nearby the local target. For surface seismic data we used Joint Migration Inversion to estimate the up- and downgoing wavefields at depth level near the target area.
|
10 |
Inversão quadrática de dados de AVO usando algoritmo de MarquardtMENDES, Ana Cristina Marques 12 November 2004 (has links)
Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-24T12:47:44Z
No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_InversaoQuadraticaDados.pdf: 975991 bytes, checksum: 5c694f4239ad08ee42219b354ed909b9 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-30T13:41:17Z (GMT) No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_InversaoQuadraticaDados.pdf: 975991 bytes, checksum: 5c694f4239ad08ee42219b354ed909b9 (MD5) / Made available in DSpace on 2014-09-30T13:41:17Z (GMT). No. of bitstreams: 2
license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5)
Dissertacao_InversaoQuadraticaDados.pdf: 975991 bytes, checksum: 5c694f4239ad08ee42219b354ed909b9 (MD5)
Previous issue date: 2004 / Neste trabalho são apresentadas expressões exatas e aproximações quadráticas das conhecidas equações de Knott-Zöeppritz, que calculam as amplitudes dos coeficientes Rpp e Rsp em termos dos contrastes médios relativos, bem como é procedida a inversão dos parâmetros de dados de AVO a partir destas aproximações quadráticas. Nesta inversão é utilizado o algoritmo de Levenberg-Marquardt, e são considerados apenas os eventos refletidos Rpp e convertidos Rsp, não associados. Nos estudos dos parâmetros físicos dos meios contrastes de impedância (δz), módulo de cisalhamento (δμ) e velocidade da onda p (δα), verificou-se quais desses parâmetros podem ser invertidos. Os resultados obtidos mostram que o contraste de impedância (δz) é muito bem resolvido estando ele relacionado com o contraste de velocidade da onda p (δα) ou com o contraste de cisalhamento (δμ), no caso de eventos refletidos considerando modelos de alto,
moderado e baixo contrastes. Por outro lado ao se fixar o contraste de impedância e relacionar os outros dois parâmetros em consideração, os resultados mostram que esses são mal resolvidos, ou seja, a região de ambiguidade torna-se muito grande e os parâmetros tornam-se ambíguos e instáveis. No caso do evento convertido e na combinação do evento refletido com o convertido, para os modelos de baixo e moderado contrastes, (δz) é muito bem resolvido, caso que não acontece para modelo de alto contraste. Diante desses resultados verifica-se que no procedimento de inversão quadrática de dados de AVO, fixado (δμ), a recuperação dos dois parâmetros variados é muito boa, no caso do evento refletido, e razoavelmente boa no caso do evento convertido, por esse motivo optou-se pela fixação do módulo de cisalhamento. / In this work, the exact and quadratic approach expression of Knott-Zöeprtiz are presented, as well as its coefficients, amplitudes, Rpp and Rsp, respectively, all of them as a function of the relative contrast, as well as the inversion of quadratic parameters of AVO data. In this inversion of the quadratic parameters, we have used the non-linear Levenberg-Marquadt algorithm, considering only converted and reflected events, either in the non-combined form. In the study of the physical parameters impedance contrast, shear modulus and P wave velocity, we have werified which of those parameters can be effectively inverted. The results have shown that the impedance contrast (δz) is very well solved when it is related with the P wave velocity contrast (δα) or with the shear contrast (δμ), considering reflected events in the case of high, moderate and low model contrasts. On the other hand, when the impedance constrast is fixed and the two other parameters in consideration are related, the results show that these parameters are bad solved, e.g., the ambiguity region becomes large and the parameters become ambiguous and unstable. In the case of a converted event and in the combination of converted and reflected events, considering low and moderate contrast models (δz) is well solved, which is not the case when we consider a high contrast model.
|
Page generated in 0.0444 seconds