Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-26T12:54:11Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MapeamentoRelevoEmbasamento.pdf: 11738916 bytes, checksum: bd08738ac99694bdec77ea1cc8ca6a2c (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-27T12:30:47Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MapeamentoRelevoEmbasamento.pdf: 11738916 bytes, checksum: bd08738ac99694bdec77ea1cc8ca6a2c (MD5) / Made available in DSpace on 2014-06-27T12:30:47Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MapeamentoRelevoEmbasamento.pdf: 11738916 bytes, checksum: bd08738ac99694bdec77ea1cc8ca6a2c (MD5) Previous issue date: 1998 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Apresentamos três novos métodos estáveis de inversão gravimétrica para estimar o relevo de uma interface arbitrária separando dois meios. Para a garantia da estabilidade da solução, introduzimos informações a priori sobre a interface a ser mapeada, através da minimização de um (ou mais) funcional estabilizante. Portanto, estes três métodos se diferenciam pelos tipos de informação físico-geológica incorporados. No primeiro método, denominado suavidade global, as profundidades da interface são estimadas em pontos discretos, presumindo-se o conhecimento a priori sobre o contraste de densidade entre os meios. Para a estabilização do problema inverso introduzimos dois vínculos: (a) proximidade entre as profundidades estimadas e verdadeiras da interface em alguns pontos fornecidas por furos de sondagem; e (b) proximidade entre as profundidades estimadas em pontos adjacentes. A combinação destes dois vínculos impõe uma suavidade uniforme a toda interface estimada, minimizando, simultaneamente em alguns pontos, os desajustes entre as profundidades conhecidas pelas sondagens e as estimadas nos mesmos pontos. O segundo método, denominado suavidade ponderada, estima as profundidades da interface em pontos discretos, admitindo o conhecimento a priori do contraste de densidade. Neste método, incorpora-se a informação geológica que a interface é suave, exceto em regiões de descontinuidades produzidas por falhas, ou seja, a interface é predominantemente suave porém localmente descontínua. Para a incorporação desta informação, desenvolvemos um processo iterativo em que três tipos de vínculos são impostos aos parâmetros: (a) ponderação da proximidade entre as profundidades estimadas em pontos adjacentes; (b) limites inferior e superior para as profundidades; e (c) proximidade entre todas as profundidades estimadas e um valor numérico conhecido. Inicializando com a solução estimada pelo método da suavidade global, este segundo método, iterativamente, acentua as feições geométricas presentes na solução inicial; ou seja, regiões suaves da interface tendem a tornar-se mais suaves e regiões abruptas tendem a tornar-se mais abruptas. Para tanto, este método atribui diferentes pesos ao vínculo de proximidade entre as profundidades adjacentes. Estes pesos são automaticamente atualizados de modo a acentuar as descontinuidades sutilmente detectadas pela solução da suavidade global. Os vínculos (b) e (c) são usados para compensar a perda da estabilidade, devida à introdução de pesos próximos a zero em alguns dos vínculos de proximidade entre parâmetros adjacentes, e incorporar a informação a priori que a região mais profunda da interface apresenta-se plana e horizontal. O vínculo (b) impõe, de modo estrito, que qualquer profundidade estimada é não negativa e menor que o valor de máxima profundidade da interface conhecido a priori; o vínculo (c) impõe que todas as profundidades estimadas são próximas a um valor que deliberadamente viola a profundidade máxima da interface. O compromisso entre os vínculos conflitantes (b) e (c) resulta na tendenciosidade da solução final em acentuar descontinuidades verticais e apresentar uma estimativa suave e achatada da região mais profunda. O terceiro método, denominado mínimo momento de inércia, estima os contrastes de densidade de uma região da subsuperfície discretizada em volumes elementares prismáticos. Este método incorpora a informação geológica que a interface a ser mapeada delimita uma fonte anômala que apresenta dimensões horizontais maiores que sua maior dimensão vertical, com bordas mergulhando verticalmente ou em direção ao centro de massa e que toda a massa (ou deficiência de massa) anômala está concentrada, de modo compacto, em torno de um nível de referência. Conceitualmente, estas informações são introduzidas pela minimização do momento de inércia das fontes em relação ao nível de referência conhecido a priori. Esta minimização é efetuada em um subespaço de parâmetros consistindo de fontes compactas e apresentando bordas mergulhando verticalmente ou em direção ao centro de massa. Efetivamente, estas informações são introduzidas através de um processo iterativo inicializando com uma solução cujo momento de inércia é próximo a zero, acrescentando, em cada iteração, uma contribuição com mínimo momento de inércia em relação ao nível de referência, de modo que a nova estimativa obedeça a limites mínimo e máximo do contraste de densidade, e minimize, simultaneamente, os desajustes entre os dados gravimétricos observados e ajustados. Adicionalmente, o processo iterativo tende a "congelar" as estimativas em um dos limites (mínimo ou máximo). O resultado final é uma fonte anômala compactada em torno do nível de referência cuja distribuição de constraste de densidade tende ao limite superior (em valor absoluto) estabelecido a priori. Estes três métodos foram aplicados a dados sintéticos e reais produzidos pelo relevo do embasamento de bacias sedimentares. A suavidade global produziu uma boa reconstrução do arcabouço de bacias que violam a condição de suavidade, tanto em dados sintéticos como em dados da Bacia do Recôncavo. Este método, apresenta a menor resolução quando comparado com os outros dois métodos. A suavidade ponderada produziu uma melhoria na resolução de relevos de embasamentos que apresentam falhamentos com grandes rejeitos e altos ângulos de mergulho, indicando uma grande potencialidade na interpretação do arcabouço de bacias extensionais, como mostramos em testes com dados sintéticos e dados do Steptoe Valley, Nevada, EUA, e da Bacia do Recôncavo. No método do mínimo momento de inércia, tomou-se como nível de referência o nível médio do terreno. As aplicações a dados sintéticos e às anomalias Bouguer do Graben de San Jacinto, California, EUA, e da Bacia do Recôncavo mostraram que, em comparação com os métodos da suavidade global e ponderada, este método estima com excelente resolução falhamentos com pequenos rejeitos sem impor a restrição da interface apresentar poucas descontinuidades locais, como no método da suavidade ponderada. / We present three new stable gravity inversion methods to estimate the relief of an interface separating two media. Solution stability is attained by introducing a priori information about the interface, through the minimization of one (or more) stabilizing functional. These methods are, therefore, characterized by the physical and geological information incorporated to the problem. The first method, named global smoothness, estimates the depths to the interface at discrete points by assuming that the density contrast between the media is known. To stabilize the inverse problem, we introduce two different constraints: (a) proximity between the true and estimated interface depths at a few isolated points, and (b) proximity between the estimated depths at adjacent points. The combination of these two constraints impose a uniform degree of smoothness all over the estimated interface, minimizing, simultaneously, the misfit between the known and estimated depths at a few boreholes, for example. The second method, named weighted smoothness, estimates the interface depths at discrete points, assuming that the density contrast is known a priori. In this method, it is incorporated the information that the interface is smooth almost everywhere, but at a few fault discontinuities. To incorporate this attribute to the estimated relief, we developed an iterative process where three kinds of constraints are imposed on parameters: (a) weighted smoothness between values of adjacent parameters, (b) lower and upper bounds on the estimated depths, and (c) proximity between the values of the parameters and a known numerical value. Starting with an initial solution produced by the global smoothness method, this method enhances initially estimated geometric features of the interface; that is, flat areas will tend to become flatter and steep areas will tend to become steeper. This is accomplished by weighting the constraints which require proximity between adjacent parameters. The weights are updated at each iteration so as to enhance the discontinuities detected in a subtle way by the global smoothness method. Constraints (b) and (c) are used both to compensate for the decrease in solution stability due to the introduction of small weights, and to reinforce flatness at the basin bottom. Constraint (b) imposes that any depth be nonnegative and smaller than an a priori known maximum depth value whereas constraint (c) imposes that all depths be closest to a value deliberately violating the maximum depth. The trade-off between these conflicting constraints is attained with a final relief presenting fiat bottom and steep borders. The third method, named minimum moment of inertia, estimates the density contrasts of a subsurface region discretized into elementary prismatic cells. It incorporates the geological information that the interface to be mapped encompasses an anomalous source which besides presenting horizontal extents much larger than its largest vertical extent, exhibits bordes dipping either vertically or toward the center of mass, and that most of the anomalous mass (or mass deficiency) is concentrated, in a compact way, about a reference level. Conceptually, these information are introduced through the minimization of the moment of inertia of the anomalous sources with respect to a reference level coinciding with the mean topographic surface. This minimization is performed in a subspace of parameters consisting of compact sources and presenting bordes which dip either vertically or toward the ce4ter of mass. Effectivelly, these informations are introduced by means of an iterative process starting with a tentative solution dose to the null solution, and adds, at each iteration, a contribution which has minimum moment of inertia with respect to the reference level, in such a way that the estimate of the next iteration does not violate the bounds on the density contrast and minimizes, at the same time, the misfit between the observed and the fitted data. Additionally, the iterative process "freezes" a density estimate if it becomes very dose to either bound. The final solution at the end of the iterative process is an estimated solution exhibiting a compact mass distribution concentrated about the reference level, whose density contrast distribution is dose to the upper (in absolute value) bound established a priori. All three methods were applied to synthetic and field gravity data, produced, respectively, by simulated and real sedimentary basins. The global smoothness method produced a good reconstruction of the basin structural framework even when the true basements were not globally smooth, as was the case of the Recôncavo Basin, Brazil. This method presents, however, the lowest resolution as compared with the other two methods. The weighted smoothness method improved the resolution of basements presenting disontinuities produced by gravity faults with large vertical offsets. It is, therefore, potentially useful in interpreting the structural framework of extensional basins as illustrated both with synthetic data and data from the Steptoe Valley, Nevada, USA and from Recôncavo Basin, Brazil. The minimum moment of inertia method was also applied to synthetic data and data from Recôncavo Basin and from San Jacinto Graben, California, USA. The results showed that, as compared with the other two methods, this method produces excellent estimates of a basement relief consisting of several adjacent discontinuities with small vertical offsets. This is a remarkable advantage over the weighted smoothness method which requires that the interface present few, local discontinuities with large vertical offsets.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-13T11:57:49Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ConstrucaoBancoDados.pdf: 10717820 bytes, checksum: 1039fea1b17d0457f0726cbb2cd35a4d (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-30T12:52:51Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ConstrucaoBancoDados.pdf: 10717820 bytes, checksum: 1039fea1b17d0457f0726cbb2cd35a4d (MD5) / Made available in DSpace on 2014-06-30T12:52:51Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ConstrucaoBancoDados.pdf: 10717820 bytes, checksum: 1039fea1b17d0457f0726cbb2cd35a4d (MD5) Previous issue date: 1994 / PETROBRAS - Petróleo Brasileiro S.A. / Neste trabalho desenvolveu-se uma solução para o clássico problema da determinação automática de litologias. Construiu-se um banco de dados utilizando-se a técnica de análise discriminante aplicada sobre trechos testemunhados de poços perfilados de um campo situado em lâmina d'água profunda localizado na bacia de Campos. A técnica de análise de agrupamento, utilizada por um grande número de pesquisadores, foi testada para o mesmo conjunto de dados porém não forneceu resultados satisfatórios em comparação com os obtidos com a análise discriminante. Os dados de perfis, após sofrerem correções ambientais, são colocados em profundidade com os dados de testemunhagem e determina-se uma função discriminante que classifica as respostas dos perfis de acordo com as litofácies descritas nos testemunhos. Dessa forma obtém-se um banco de dados, a partir do qual se faz a discriminação de litologias em poços, testemunhados ou não, da mesma área. O banco de dados assim construído permite sua atualização através da incorporação de trechos testemunhados de novos poços. Mostra-se a necessidade de agrupar as litologias em quatro litofácies principais para este campo de maneira a padronizar os resultados, diminuir as indeterminações e melhorar a apresentação final. Como produto final obtém-se uma curva discriminada de eletrofácies que reproduz, com algum grau de certeza, as litofácies que reproduzem as litologias esperadas em subsuperfície. Como exemplo de aplicação deste método, utilizou-se as curvas discriminadas para correlacionar três poços, utilizando para isto um algoritmo de correlação automática. / In this work we develop a solution for the classic problem of automatic lithology determination from wireline logs. A data bank generated using the discriminant analysis technique, applied over cored intervals of logged wells from a deep water field situated in the Campos basin. The clustering technique, although preferred by many workers in this field, was found to give inferior results compared to discriminant analysis. The log data, after going through environmental corrections, were placed on depth with the cored data, and a discriminant function was computed. This function classifies the log response according to the facies described from the cores. Thus we obtain a data bank, whitch allow us to perform lithologic discrimination on cored or uncored wells, in the same area. We show that the litology has to be clustered into four facies in order to get optimum results. The final product is a discrete (discriminated) curve which reproduces, reasonably well the subsurface litology. As an aplication of this method, we use the discrete (discriminated) curves to correlate three wells using an automatic correlation algorithm.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-16T16:17:08Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPocos.pdf: 13431901 bytes, checksum: 51e2f3d03ff48d53ea67eaaa7480b3b3 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-30T13:25:37Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPocos.pdf: 13431901 bytes, checksum: 51e2f3d03ff48d53ea67eaaa7480b3b3 (MD5) / Made available in DSpace on 2014-06-30T13:25:37Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPocos.pdf: 13431901 bytes, checksum: 51e2f3d03ff48d53ea67eaaa7480b3b3 (MD5) Previous issue date: 1994 / UFPA - Universidade Federal do Pará / PETROBRAS - Petróleo Brasileiro S.A. / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FINEP - Financiadora de Estudos e Projetos / As medidas de resistividade são de fundamental importância para o cálculo da saturação de óleo em reservatórios potencialmente produtores. A combinação das medidas de resistividade rasa e profunda permite a obtenção dos parâmetros Rt, Rxo e di. Mas, em reservatórios complexos existem dificuldades em se obter leituras confiáveis de Rt, devido à baixa resolução vertical das ferramentas de investigação profunda. Em reservatórios laminados, por exemplo, as leituras obtidas pela ferramenta de indução profunda (ILD) podem levar a uma interpretação errônea das mesmas, levando a acreditar que as medidas obtidas do perfil referem-se a uma única camada. Este problema pode ser em parte resolvido através de uma metodologia que melhore a resolução vertical dos perfis de investigação profunda, valendo-se do uso de informações obtidas de um perfil de alta resolução vertical, i.e; a curva de resistividade rasa. Uma abordagem neste sentido seria usar um perfil de alta resolução que apresente uma boa correlação com o perfil de investigação profunda. Esta correlação pode ser melhor avaliada se aplicarmos um filtro no perfil de alta resolução, de tal forma que o perfil resultante tenha teoricamente a mesma resolução vertical do perfil de baixa resolução. A obtenção deste filtro, porém, recai na premissa de que as funções respostas verticais para as ferramentas de alta e baixa resolução são disponíveis, o que não ocorre na prática. Este trabalho se propõe mostrar uma nova abordagem onde o filtro pode ser obtido a partir de um tratamento no domínio da freqüência. Tal tratamento visa igualar a energia espectral do perfil de alta resolução à energia espectral do perfil de baixa resolução tendo como base o Teorema de Parseval. Será mostrado que a resolução vertical depende fundamentalmente da energia espectral do perfil em questão. A seguir, uma regressão linear é aplicada sobre os perfis de alta resolução filtrado e de baixa resolução. Para cada ponto amostrado dos perfis, uma rotina de minimização é aplicada visando escolher o melhor intervalo de correlação entre os perfis. Finalmente, um fator de correção é aplicado sobre cada ponto do perfil de baixa resolução. Os resultados obtidos com os perfis de indução são promissores, demonstrando a eficácia da abordagem e mesmo quando aplicada em perfis com diferentes propriedades petrofísicas, a metodologia funcionou satisfatoriamente, sem danificar os perfis originais. / Resistivity measurements are of fundamental importance for the calculation of oil saturation in potentially producing reservoirs. The combined measurement of shallow and deep resistivities enables the determination of the parameters Rt, Rxo and di. But in complex reservoirs we have difficulty in obtaining a confident reading of Rt, due to the low vertical resolution of deep reading tools. In laminated reservoirs, for example, the deep induction reading ILD can be interpreted erroneously with the belief that the measurement refers to one bed. This may be true for extreme case of thick beds, but more often is not. This problem can be partly resolved by enhancement of the vertical resolution of the deep reading log through comparison with the high resolution (shallow resistivity) log. One approach is to use a high resolution log where there is good correlation with the deep reading log. This correlation can be better evaluated if we apply a filter to the high resolution log such that theoretically the resultant log has the same vertical resolution as the low resolution log. However, this assumes that the vertical response of the high and low resolution tools are available, and in practice this is often not the case. In this study we attempt to demonstrate an alternative approach where the filter can be obtained from consideration of the frequency domain. The technique compares the spectral energy of the high and low resolution logs. It is shown that the vertical resolution depends fundamentally on the spectral energy of the actual log based on Parseval Theorem. Next a linear regression is applied to the filtered high resolution and low resolution logs and for each depth a minimisation routine is applied to determine the best correlation interval between the logs. Finally a correction factor is applied to each point on the low resolution log. This factor is considered by the correlation coefficients over the interval minimised for each point. The results obtained with induction logs are promising and the metodology shoud be aplyed on diferent logs techniques.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-16T13:40:45Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPoco.pdf: 7959727 bytes, checksum: 95aa85fbc47f9a900fe029997cb02558 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-30T13:33:02Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPoco.pdf: 7959727 bytes, checksum: 95aa85fbc47f9a900fe029997cb02558 (MD5) / Made available in DSpace on 2014-06-30T13:33:02Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPoco.pdf: 7959727 bytes, checksum: 95aa85fbc47f9a900fe029997cb02558 (MD5) Previous issue date: 1992 / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FINEP - Financiadora de Estudos e Projetos / A maioria dos perfis de poço utilizados nas avaliações petrofísicas de reservatórios possuem uma resolução vertical na ordem de um metro. Isto cria um problema quando as espessuras típicas das camadas são inferiores a um metro, uma vez que não há correção das leituras. Os perfis de alta resolução vertical como da ferramenta de propagação eletromagnética (EPT, Schlumberger), o dipmeter (SHDT, Schlumberger) ou das ferramentas de varredura acústica ou elétrica possuem uma resolução vertical da ordem de centimetros, mas apresentam uma limitada aplicação para as avaliações petrofísicas. Nós apresentamos um método para a deconvolução de um perfil de baixa resolução vertical que utiliza informações de um perfil de alta resolução vertical para identificar uma nítida interface entre camadas que apresentam valores da propriedade petrofísica contrastante, mas localmente constante em ambos os lados. A partir desse intervalo de controle, nós determinamos a função resposta vertical da ferramenta sob as condições atuais do poço com base no teorema da convolução. Utilizamos várias interfaces de modo a obter valores mais representativos da resposta da ferramenta. O perfil de baixa resolução é então deconvoluido utilizando a transformada discreta de Fourier (FFT) sobre todo o intervalo de interesse. É importante destacar que a invasão do filtrado da lama e a presença do bolo de lama não produzem efeitos danosos sobre o método, que foi aplicado a perfis sintéticos e a dados de campo, onde a aplicação de filtros com um correto ajuste de profundidade, bem como a própria escolha do intervalo de controle, antes da deconvolução, são de extrema importância para o sucesso do método. / Most wireline logs used for the petrophysical evaluation of reservoirs have a vertical resolution in the order of one meter. This poses a problem when the typical layer thickness is less than one meter, since no correct reading will be obtained. High-resolution logs like the electromagnetic propagation tool (EPT, Schlumberger) the dipmeters (SHDT, Schlumberger) or acoustic and electrical imaging devices have a resolution well below one meter, sometimes as high as one centimeter, but their application to petrophysical reservoir evaluation is limited. We present an approach which uses information from high-resolution logs to deconvolve low-resolution logs. From the high-resolution log we first identify a sharp bed boundary, on both sides of which the petrophysical properties assume a locally constant but different value (an interface, for instance). From this control interval we determine the tool response function under actual borehole conditions of the low-resolution log using the basic convolution theorem. It can be obtained at various sharp bed boundaries in order to obtain a more representative tool response. The low-resolution log is then deconvolved using Fast Fourier Transforms (FFT) over the entire interval of interest. Neither the invasion of mud filtrate or the presence of the mud cake affect this method adversely. We apply it to a number of synthetic and field data sets. Noise filtering prior to deconvolution is found to be important, as is correct depth shifting and a proper choice of the control interval.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-16T13:48:22Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPocoRede.pdf: 1991292 bytes, checksum: 11d42098d8613d8bb1eb827e88db1f71 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-30T14:50:28Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPocoRede.pdf: 1991292 bytes, checksum: 11d42098d8613d8bb1eb827e88db1f71 (MD5) / Made available in DSpace on 2014-06-30T14:50:28Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeconvolucaoPerfisPocoRede.pdf: 1991292 bytes, checksum: 11d42098d8613d8bb1eb827e88db1f71 (MD5) Previous issue date: 2006 / ANP - Agência Nacional do Petróleo / Para a indústria do petróleo, a interpretação dos perfis de poço é a principal fonte de informação sobre a presença e quantificação de hidrocarbonetos em subsuperfície. Entretanto, em duas situações as novas tecnologias, tanto em termos do processo construtivo das ferramentas, quanto da transmissão dos dados não têm justificativa econômica, ensejando a utilização de um conjunto de perfis convencionais: reavaliações de campos maduros e avaliações de campos marginais. Os procedimentos de aquisição dos perfis convencionais podem alterar o valor da propriedade física bem como a localização dos limites verticais de uma camada rochosa. Este é um antigo problema na geofísica de poço – o paradoxo entre a resolução vertical e a profundidade de investigação de uma ferramenta de perfilagem. Hoje em dia, isto é contornado através da alta tecnologia na construção das novas ferramentas, entretanto, este problema ainda persiste no caso das ferramentas convencionais como, a ferramenta de raio gama natural (GR). Apresenta-se, neste trabalho, um novo método para atenuar as alterações induzidas no perfil pela ferramenta, através da integração do clássico modelo convolucional do perfil com as redes neurais recorrentes. Assume-se que um perfil de poço pode ser representado através da operação de convolução em profundidade entre a variação da propriedade física da rocha (perfil ideal) e uma função que representa a alteração produzida sobre a propriedade física, chamada como resposta vertical da ferramenta. Assim, desenvolve-se um processamento iterativo dos perfis, o qual atua na forma da operação de deconvolução, composto por três redes neurais recorrentes. A primeira visa estimar a resposta vertical da ferramenta; a segunda procura definir os limites verticais de cada camada rochosa e a última é construída para estimar o valor real da propriedade física. Este processamento é iniciado com uma estimativa externa tanto para o perfil ideal, quanto para a resposta vertical da ferramenta. Finalmente, mostram-se as melhorias na resolução vertical e na avaliação da propriedade física produzida por esta metodologia em perfis sintéticos e em perfis reais da formação Lagunillas, bacia do Lago Maracaibo, Venezuela. / For oil industry, the logs analysis is the main information source about the presence and quantification of hydrocarbon in subsurface. However, in two situations the new logging technologies are not economically viable and conventional logging tools must be used: The reevaluation of mature oil fields and evaluation of marginal oil fields. In conventional logs its data acquisition procedure may blur the value of physical property and the vertical limits of a rock layer. We are talking about an old problem in well logging – The paradox between vertical resolution and depth of investigation of a logging tool. Nowadays it is well handling by the high technology of new tools, but this problem persists in conventional old tools, e.g. natural gamma ray log (GR). Here, we present a method to smooth this kind of linear distortion in well logs by an integration of classical well log convolution model with recurrent neural networks. We assume that a well log can be well represented by an in depth convolution operation between the variation of rock physical property (ideal log) and a function that causes the distortion, called as vertical tool response. Thus, we develop an iterative data processing, which acts as a deconvolution operation, composed by three recurrent neural networks. The first one seeks to estimate the vertical tool response; the second one search for the vertical limits definition of each rock layer and the last one is constructed to estimate the actual physical property. To start this process we supply an appropriated first guess of ideal log and vertical tool response. Finally, we show the improvements in vertical resolution and in the physical property evaluation produced by this methodology in synthetic logs and actual well log data from Lagunillas formation, Maracaibo basin, Venezuela.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-30T12:55:35Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ModelagemSismicaInversao.pdf: 14582851 bytes, checksum: 072884b8cac95466365f0e2a5433aedc (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexação de assuntos on 2014-08-06T16:11:05Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-07T13:23:40Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ModelagemSismicaInversao.pdf: 14582851 bytes, checksum: 072884b8cac95466365f0e2a5433aedc (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-08T14:46:56Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ModelagemSismicaInversao.pdf: 14582851 bytes, checksum: 072884b8cac95466365f0e2a5433aedc (MD5) / Made available in DSpace on 2014-08-08T14:46:56Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ModelagemSismicaInversao.pdf: 14582851 bytes, checksum: 072884b8cac95466365f0e2a5433aedc (MD5) Previous issue date: 1993 / UFPA - Universidade Federal do Pará / PETROBRAS - Petróleo Brasileiro S.A. / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A necessidade da adoção de modelos elásticos anisotrópicos, no contexto da sísmica de exploração, vem crescendo com o advento de novas técnicas de aquisição de dados como VSP, walkway VSP, tomografia poço a poço e levantamentos sísmicos com grande afastamento. Meios anisotrópicos, no contexto da sísmica de exploração, são modelos efetivos para explicar a propagação de ondas através de meios que apresentam padrões de heterogeneidade em escala muito menor que o comprimento de onda das ondas sísmicas. Particularmente, estes modelos são muito úteis para explicar o dado sísmico mais robusto que são as medidas de tempo de trânsito. Neste trabalho, são investigados aspectos da propagação de ondas, traçado de raios e inversão de tempos de trânsito em meios anisotrópicos. É estudada a propagação de ondas SH em meios anisotrópicos estratificados na situação mais geral onde estas ondas podem ocorrer, ou seja, em meios monoclínicos com um plano vertical de simetria especular. É mostrado que o campo de ondas SH refletido a partir de um semi-espaço estratificado, não apresenta qualquer informação sobre a possível presença de anisotropia em subsuperfície. São apresentados métodos simples e eficientes para o traçado de raios em 3D através de meios anisotrópicos estratificados, baseados no princípio de Fermat. Estes métodos constituem o primeiro passo para o desenvolvimento de algoritmos de inversão de tempos de trânsito para meios anisotrópicos em 3D, a partir de dados de VSP e walkaway VSP. Esta abordagem é promissora para determinação de modelos de velocidade, que são necessários para migração de dados sísmicos 3D na presença de anisotropia. É efetuada a análise da inversão tomográfica não linear, para meios estratificados transversalmente isotrópicos com um eixo de simetria vertical(TIV). As limitações dos dados de tempo de trânsito de eventos qP para determinação das constantes elásticas, são estabelecidas e caracterizados os efeitos da falta de cobertura angular completa na inversão tomográfica. Um algoritmo de inversão foi desenvolvido e avaliado em dados sintéticos. A aplicação do algoritmo a dados reais demonstra a consistência de meios TIV. Esta abordagem é útil para casos onde há informação a priori sobre a estratificação quase plana das formações e onde os próprios dados do levantamento poço a poço apresentam um alto grau de simetria especular em relação a um plano vertical. Também pode ser útil em interpretações preliminares, onde a estimativa de um meio estratificado, serve como modelo de fundo para se efetuar análises mais detalhadas, por exemplo, como um modelo de velocidades anisotrópico para migração, ou como um modelo de calibração para análises de AVO. / The assumption of anisotropic elastic models, in the context of exploration seismology, has been increasing since the advent of new aquisition tecniques: VSP, walkway VSP, crosshole tomography and large offset seismic surveys. Effective anisotropic elastic models can account for patterns of inhomogeneities at a scale much lower than the wavelength of the propagating seismic energy. Particularly, effective media can account for the most robust seismic data, i.e., traveltime measurements. This work investigates some aspects of wave propagation, raytracing and traveltime inversion in anisotropic media. The propagation of SH waves in a layered anisotropic medium in the most general case where these waves can occur, i.e., propagation in the mirror symmetry plane of a monoclinic medium, is studied. It is shown that SH reflected field from a stratified half-space is 'blind' to the possible anisotropy below. Fast ray tracing procedures for stratified anisotropic media in 3D are presented. Such ray tracing codes are a first step towards the development of traveltime inversion algorithms for non-flat anisotropic layers in 3D. Traveltime inversion of VSP and walkway VSP is a promissing approach to 3D velocity model building, starting from a well location. These velocity models are necessary for the migration of seismic data in the presence of anisotropy. The analysis of the nonlinear tomographic inversion, for a vertically inhomogeneous transversally isotropic medium with vertical axis of symmetry (TIV), is presented. The limitations of qP traveltime data are pointed as well as the consequences of the lack of full angular ray coverage for tomographic inversion. An algorithm for tomographic inversion is presented and evaluated in synthetic data seis. Application to real data is presented. This approach is attractive for cases where the formations are known a priori to be relatively flat, and where crosswell data itself shows a high degree of left-right symmetry. It also may be of use for preliminary surveys, where the layered estimate can be used as a background model to carry out more detailed analysis, e.g., as an anisotropic velocity model for migration, or as a calibration model for AVO analysis.

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.

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.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-06T16:39:28Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoAnaliseSeries.pdf: 9563203 bytes, checksum: 19126d425ed2de3d0aeb62c84c46a64e (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:42:22Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-20T12:25:05Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoAnaliseSeries.pdf: 9563203 bytes, checksum: 19126d425ed2de3d0aeb62c84c46a64e (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T13:20:01Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoAnaliseSeries.pdf: 9563203 bytes, checksum: 19126d425ed2de3d0aeb62c84c46a64e (MD5) / Made available in DSpace on 2014-08-20T13:20:01Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoAnaliseSeries.pdf: 9563203 bytes, checksum: 19126d425ed2de3d0aeb62c84c46a64e (MD5) Previous issue date: 1993 / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FINEP - Financiadora de Estudos e Projetos / Esta tese descreve a aplicação de análise de séries temporais em perfis de poço. Através desta técnica é possível avaliar-se a repetição e a resolução vertical de perfis, e determinar-se o intervalo de amostragem e a velocidade de perfilagem ideais para diferentes perfis. A comparação entre três poços é também feita, baseada num mesmo tipo de perfil. Para tanto, na seqüência utilizada, procurou-se manter num mesmo domínio os dados cuja quantidade total de amostras (N) por perfil não ultrapassou 2048. Desses dados, foram inicialmente retirados o valor médio das amostras e o alinhamento polinomial algébrico eventualmente nelas embutido. Em seguida, foi efetuada a aplicação do ponderador cossenoidal, do filtro passa-alta, da janela Hanning, do cálculo da função coerência, do espectro de fase, da razão sinal-ruído e dos espectros de potência do sinal e do ruído, nesta ordem. Para a função coerência, fez-se necessário o cálculo dos níveis de confiança de 50%, 90% e 95%. O cálculo do primeiro nível teve por base a necessidade de se determinar a resolução vertical de alguns perfis, e dos demais, a fim de que fosse obtida uma informação referente à localização daqueles níveis para a coerência calculada. Em relação ao espectro de fase, seu cálculo surgiu da necessidade de se obter uma informação adicional a respeito dos perfis manipulados, ou seja, o conhecimento da ocorrência ou não de deslocamento relativo de profundidade entre a seção principal e a seção repetida. A razão sinal-ruído foi calculada no sentido de possibilitar a comparação, como elemento avaliador dos diversos tipos de perfis, com a coerência e o cálculo dos espectros de potência. Os espectros de potência do sinal e do ruído foram calculados para se ter mais um parâmetro de avaliação da seção repetida, já que em tese, os espectros de potência do sinal e do ruído da seção repetida devem ser iguais aos respectivos espectros da seção principal. Os dados utilizados na aplicação da metodologia proposta foram fornecidos pela PETROBRÁS e oriundos de quatro poços da Bacia Potiguar emersa. Por questões de sigilo empresarial, os poços foram identificados como poços A, B, C e D. A avaliação da repetição entre diferentes tipos de perfis indica que, para o poço A, o perfil micro-esférico (MSFL) tem repetição melhor do que o perfil de porosidade neutrônica (CNL), o qual tem, por sua vez, repetição melhor do que o perfil de raios gama normal (GR). Para os perfis do poço D, uma diminuição da velocidade de perfilagem de 550 m/h para 275 m/h é vantajosa apenas para o perfil de porosidade neutrônica. Já a velocidade de perfilagem de 920 m/h, utilizada, na obtenção dos perfis do poço C, é totalmente inadequada para os perfis de resistividade (MSFL, ILD e ILM). A diminuição do intervalo de amostragem de 0,20 m para 0,0508 m, nos perfis de raios gama e de porosidade neutrônica, e 0,0254 m para o perfil de densidade, apresenta bons resultados quando aplicada no poço D. O cálculo da resolução vertical indica, para o perfil de porosidade neutrônica, uma superioridade qualitativa em relação ao perfil de raios gama normal, ambos pertencentes ao poço A. Para o poço C, o perfil micro-esférico apresenta uma resolução vertical na mesma ordem de grandeza da resolução do perfil de raios gama do poço B, o que evidencia ainda mais a inconveniência da velocidade de perfilagem utilizada no poço C. Já para o poço D, o cálculo da resolução vertical indica uma superioridade qualitativa do perfil de densidade de alta resolução em relação ao perfil de raios gama de alta resolução. A comparação entre os poços A, B e D, levada a efeito através dos respectivos perfis de porosidade neutrônica normais, comprova que a presença de ruído aleatório, em geral, está diretamente ligada à porosidade da formação - uma maior porosidade indica uma presença maior de ruído e, por conseguinte, uma queda qualitativa no perfil obtido. A análise do espectro de fase de cada perfil indica um deslocamento em profundidade, existente entre as seções principal e repetida de todos os perfis do poço C. E isto pôde ser confirmado com a posterior superposição das seções. / This thesis describes the application of time seris analysis to wireline logs. Through this technique it is possible to evaluate both their repeatability and vertical resolution, and determine the optimum sampling interval and acquisition speed for different logs. A comparison between three wells is also made, based on the same type of log. The sequence used, is to obtain data, in the same domain, for which the number of samples (N) does not exceed 2048. For these data the mean sample value and the algebraic polynomial are determined. The following were then applied, in order: cosine taper, high pass filter, Hanning window, the calculation of the coherence function, the phase spectra, and the signal to noise ratio of the power spectra for both the signal and the noise. For the coherence function, it was necessary to calculate the level of confidence for 50 necessary to determine the vertical resolution of some logs. The others were calculated to provide information concerning the position of the levels of coherence calculated. In relation to the phase spectra, its calculation suggests it is necessary to obtain additional information in respect of the processed logs, or in other words knowledge of any relative depth shifts made between the principal and repeat sections. The signal to noise ratio was calculated to investigate the possibility of evaluating the different types of logs by making a comparison with the coherence and the calculated power spectra. The power spectra of the signal and the noise were calculated to provide one additional parameter to evaluate the repeat section. In theory the power spectra of the signal and the noise of the repeat section should be equal to their respective spectra for the principal sections. The data used in this work were provided by PETROBRAS and originated in four boles of the Potiguar Basin. These are referred to as Holes A, B, C, and D. The evaluation of the repeatability between different types of logs indicates that, for Hole A, the microspherical log (MSFL) has better repeatability than the porosity log (CNL), and that this repeats better than the gamma ray log (GR). For the logs run in Hole D, a decrease in the speed at which the logs are run, from 550 m/hr to 275 m/hr, is advantageous only for the neutron porosity log. The velocity of 920 m/hr used in Hole C is totally inadequate for the resistivity logs (MSFL, ILD, and ILM). A decrease in the sampling interval from 0.2m to 0.0508m, for the gamma ray and neutron porosity logs, and to 0.0254m for the density log, gives good results when applied to Hole D. The calculation of vertical resolution indicates, that the neutron porosity is of superior quality when compared with the natural gamma ray in Hole A. In hole C the microspherical log (MSFL) has a vertical resolution of the same order as the natural gamma ray log (GR) in Hole B. This confirms the inappropriate speed for the logs acquired in Hole C. For Hole D, the calculation of vertical resolution indicates a superior quality for the high resolution density log when compared with the high resolution gamma ray log. The comparison between Holes A, B and D, uses the respective neutron porosity logs, confirming that the presence of noise, in general, is directly linked with the porosity of the formation - a higher porosity indicates more noise and consequently, a lower quality log is obtained. Analysis of phase spectra of each log indicates a depth shift exists between the principal and repeat sections for all the logs of Hole C. This could be confirmed through superposition of the sections.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-18T12:37:29Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoRefletoresSismicos.pdf: 19890641 bytes, checksum: 386a6dda416849e7792da409aaa75ec8 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Incluir assuntos on 2014-08-07T16:02:06Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-21T13:19:11Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoRefletoresSismicos.pdf: 19890641 bytes, checksum: 386a6dda416849e7792da409aaa75ec8 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T14:02:41Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoRefletoresSismicos.pdf: 19890641 bytes, checksum: 386a6dda416849e7792da409aaa75ec8 (MD5) / Made available in DSpace on 2014-08-22T14:02:41Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeteccaoRefletoresSismicos.pdf: 19890641 bytes, checksum: 386a6dda416849e7792da409aaa75ec8 (MD5) Previous issue date: 1999 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / As redes neurais artificiais têm provado serem uma poderosa técnica na resolução de uma grande variedade de problemas de otimização. Nesta dissertação é desenvolvida uma nova rede neural, tipo recorrente, sem realimentação (self-feedback loops) e sem neurônios ocultos, para o processamento do sinal sísmico, para fornecer a posição temporal, a polaridade e as amplitudes estimadas dos refletores sísmicos, representadas pelos seus coeficientes de reflexão. A principal característica dessa nova rede neural consiste no tipo de função de ativação utilizada, a qual permite três possíveis estados para o neurônio. Busca-se estimar a posição dos refletores sísmicos e reproduzir as verdadeiras polaridades desses refletores. A idéia básica desse novo tipo de rede, aqui denominada rede neural discreta (RND), é relacionar uma função objeto, que descreve o problema geofísico, com a função de Liapunov, que descreve a dinâmica da rede neural. Deste modo, a dinâmica da rede leva a uma minimização local da sua função de Liapunov e consequentemente leva a uma minimização da função objeto. Assim, com uma codificação conveniente do sinal de saída da rede tem-se uma solução do problema geofísico. A avaliação operacional da arquitetura desta rede neural artificial é realizada em dados sintéticos gerados através do modelo convolucional simples e da teoria do raio. A razão é para explicar o comportamento da rede com dados contaminados por ruído, e diante de pulsos fonte de fases mínima, máxima e misturada. / The artificial neural networks have proven to be a powerful technique in solving a wide variety of optimization problems. In this work, we develop a new recurrent network, with no self-feedback loops, and no hidden neurons, for seismic signal processing, where this neural network gives the true polarity, reflectors location and magnitude estimations. The main characteristic of this neural network is the use of a type of activation function which permits three possible states of neurons, to estimate the position of the seismic reflectors in such way to reproduce its true polarities. The basic idea of this new neural network type, denominated here by discrete neural network (DNN), is to relate a cost function, that describes the geophysical problem, with the Liapunov function, that describes the neural network dynamics. In this way, the dynamics of the network leads to a local minimization of the Liapunov function, and will consequently lead to a minimization of the cost function. Thus, with a convenient output signal codification of the neural network a geophysical problem solution is obtained. The operational evaluation of this neural network architecture is performed with synthetic data obtained through the simple convolutional model and seismic ray theory, and its behavior explained with additive noise in the data with minimum, maximum and mixed phase time source pulses.