Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-11-14T14:02:46Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemUnidimensionalDados.pdf: 3780841 bytes, checksum: 1ab8bca8ded68ad9b937cc088d05a174 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-17T15:54:06Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemUnidimensionalDados.pdf: 3780841 bytes, checksum: 1ab8bca8ded68ad9b937cc088d05a174 (MD5) / Made available in DSpace on 2014-11-17T15:54:06Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelagemUnidimensionalDados.pdf: 3780841 bytes, checksum: 1ab8bca8ded68ad9b937cc088d05a174 (MD5) Previous issue date: 2007 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O Marine Controlled Source ElectroMagnetic (MCSEM) é um novo método geofísico que tem sido empregado para a detecção de camadas resistivas no sub-solo oceânico e tem se mostrado uma ferramenta promissora para a indústria do petróleo. Neste trabalho, é modelado levantamentos do método MCSEM incluindo variações nos valores de resistividade da água do mar. Tais variações podem resultar da influência de correntes marinhas, gradientes de temperatura ou de qualquer fator que possa afetar a salinidade da água. Observamos que os dados do método são afetados de maneira significativa por tais variações. Tanto a variação na resistividade quanto na espessura da camada de água pode afetar os resultados do MCSEM de maneira semelhante. É modelado também o posicionamento aleatório dos receptores, mostrando um decaimento significativo da amplitude do sinal quando os receptores se posicionam a uma certa distância da linha de medidas. / The Marine Controlled Source ElectroMagnetic (MCSEM) is a new geophysical method that is very promising in the detection of resistive layers in the sediments under the oceanic floor, a useful information for the oil industry. In this work we model MCSEM data including variations in the resistivity of the ocean water. Such variations can arise from the influence of marine currents, temperature gradients or any other source of influence on the water salinity. We observe that the data from this method are significantly affected by those variations.We also simulate the random positioning of the receivers when deployed to the sea botton, our results show a strong decline in the amplitude of the signal when the receivers are located at a certain distance from the ideal array line.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-11-18T12:24:17Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoEletromagneticoAnalogico.pdf: 8231910 bytes, checksum: e06ef7849ac9f47689cabfb2f3939fe0 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-11-18T16:58:47Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoEletromagneticoAnalogico.pdf: 8231910 bytes, checksum: e06ef7849ac9f47689cabfb2f3939fe0 (MD5) / Made available in DSpace on 2014-11-18T16:58:47Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoEletromagneticoAnalogico.pdf: 8231910 bytes, checksum: e06ef7849ac9f47689cabfb2f3939fe0 (MD5) Previous issue date: 1989 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Nas últimas três décadas os métodos eletromagnéticos vem se desenvolvendo satisfatoriamente em função da aplicabilidade na prospecção de corpos de sulfetos maciços. Em regiões de climas tropicais, normalmente o manto de intemperismo apresenta-se condutivo. E este, na maioria das vezes, não é levado em conta na prospecção eletromagnética, causando portanto erros consideráveis de interpretações. Neste trabalho consideramos o manto de intemperismo em contato e sem-contato com o corpo condutor. Com o objetivo de estudar os efeitos dos mantos sobre anomalias EM de corpos tabulares inclinados, foram feitos vários experimentos utilizando modelamento analógico em escala reduzida, admitindo-se diferentes parâmetros de resposta para o corpo e o manto. Para simular o corpo foram utilizadas placas de aço inoxidável com as dimensões suficientemente grande em relação ao espaçamento entre as bobinas, de tal modo que simulassem um semi-plano. Para simular o manto foi utilizado uma solução de sais, sendo que para o caso de manto-condutivo, o corpo foi colocado em contato galvânico com a solução. Para o manto-indutivo foi considerado sem-contato galvânico, de tal forma que o corpo e o manto fossem acoplados apenas indutivamente. Onde 1) o corpo foi colocado totalmente sem contato com o manto 2) o corpo foi revestido por uma película resistiva e colocado em contato com o manto. Com a presença de manto-indutivo, observamos que a amplitude dos perfis é levemente atenuada. Além disso, observamos na quadratura a reversão e o aparecimento de um pico-extra nas inclinações do corpo θ≤60º e nos valores de número de indução αc≥78.16 e αm≥0.5 respectivamente do corpo e do manto. E a rotação de fase se dava no sentido horário, sendo mais intensa para altos valores de número de indução do corpo. No manto-indutivo o corpo parece estar a uma profundidade maior que a verdadeira, e ser mais condutivo do que realmente é. Com manto-condutivo, observamos que as amplitudes dos perfis são ligeiramente acrescidos assim como, a rotação de fase se dava no sentido anti-horário sendo mais intensa para pequenos valores de número de indução do corpo. Os demais efeitos tais como reversão na quadratura e presença de pico-extra ocorrem de modo análogo ao ocorrido no caso de manto-indutivo. No manto-condutivo, o corpo parece estar a uma profundidade inferior à verdadeira e ser menos condutivo. As anomalias EM são ligeiramente modificadas em função da rotação de fase e atenuação de amplitudes que ocorrem nos campos primário e secundário quando atravessam o manto, e também em consequência da interação indutiva de corrente induzida entre corpo e manto. Além disso ocorre a redistribuição de corrente no manto devido à presença de corpo dentro do manto. No manto-condutivo as correntes são canalizadas dentro do corpo que está em contato galvânico com o manto, enquanto no manto-indutivo ocorre um desvio de corrente, devido à película resistiva que envolve o corpo. / In the last three decades, electromagnetic methods are continuing to develop due to their proven usefulness in the search for massive sulfides. In tropical regions, the overburden is usually conducting. However the effect of this conducting layer is seldon taken into account while interpretating the EM data. this can cause considerable interpretation erros. In this work, we considered the overburden in-contact and not in-contact with the underlying conductor. Reduced scale model experiments were carried out to study the effects of the overburden on the EM anomalies of inclined tabular like bodies. With this objective, response parameters of the conductor and the overburden were varied. The target-conductor is represented in the model experiments by stainless steel sheets. These sheets are, larg compared the coil separation in order to simulated the behavior of a half-plane. The overburden is simulated by a conducting solution. To represented the case of a conductive-overburden the model is placed in galvanic contact with the solution. In case of the inductive-overburden the model is coupled to the overburden only inductively and did not have a galvanic contact. This situation is attained by keeping the model 1) totally out of the solution, and 2) partially in the solution but the model is covered by a resistive film to avoid the galvanic contact with the solution. In the presence of inductive-overburden, the anomalies are slightly attenuated. In the quadrature component the profile is inverted and an extra-peak appear in the case of low dipping models at higher induction numbers of the overburden and the model. Also the anomaly suffers a clockwise rotation which is intense for higher induction number of the conductor. Therefore, in the presence of an inductive-overburden the conductor appears to be more conducting and at higher depth than it really is. In the presence of a conductive-overburden, the anomaly amplitude are slightly enhanced and suffer an anticlockwise phase rotation, which is more intense at the lower induction number of the conductor. However, the other effects in the quadrature component, like the inversion and an appearance of the extra-peak, are similar to that of the inductive-overburden. Due to these effects of the conductive-overburden, the conductor appears to be at a shallower depth and less conducting than it actually is. These modifications in anomalies are caused by: a) The primary and secondary EM fields suffering attenuation and phase rotation when passing through a conducting overburden, b) an inductive interaction between the induced currents in the conductor and the overburden, and also, c) a redistribution of the currents in the overburden when the conductor is placed in it. In case of a conductive-overburden, currents are channeled in the conductor which is in galvanic contact with the less conducting overburden. On the other hand when a conductor covered with a resistive film is placed in the conducting solution, the currents are displaced and they crowd in the solution next to the boundary of the resistive film.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-12-03T14:22:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ModelamentoPermissividadeDieletrica.pdf: 9420207 bytes, checksum: 28f0ebee4414f2a902ff9e96c0b0b87e (MD5) / Rejected by Albirene Aires (albireneufpa@gmail.com), reason: Favor incluir as palavras chave on 2014-12-05T11:45:19Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-12-10T00:14:49Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPermissividadeDieletrica.pdf: 9420207 bytes, checksum: 28f0ebee4414f2a902ff9e96c0b0b87e (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-12-10T17:20:36Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPermissividadeDieletrica.pdf: 9420207 bytes, checksum: 28f0ebee4414f2a902ff9e96c0b0b87e (MD5) / Made available in DSpace on 2014-12-10T17:20:36Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPermissividadeDieletrica.pdf: 9420207 bytes, checksum: 28f0ebee4414f2a902ff9e96c0b0b87e (MD5) Previous issue date: 1990 / A ferramenta de propagação eletromagnética (EPT) fornece o tempo de propagação (T_pl) e a atenuação (A) de uma onda eletromagnética que se propaga num meio com perdas. Estas respostas da EPT são funções da permissividade dielétrica do meio. Existem vários modelos e fórmulas de misturas sobre a permissividade dielétrica de rochas reservatório que podem ser utilizados na interpretação da ferramenta de alta frequência. No entanto, as fórmulas de mistura não consideram a distribuição e a geometria do espaço poroso, e estes parâmetros são essenciais para que sejam obtidas respostas dielétricas mais próximas de uma rocha real. Foi selecionado um modelo baseado nos parâmetros descritos acima e este foi aplicado à dados dielétricos disponíveis na literatura. Foi obtida uma boa concordância entre as curvas teóricas e os dados experimentais, comprovando assim que a distribuição e a geometria dos poros têm que ser levadas em conta no desenvolvimento de um modelo realista. Foram conseguidas também funções de distribuição de razão de aspecto de poros, através das quais geramos várias curvas relacionando as respostas da EPT com diversas saturações de óleo/gás. Estas curvas foram aplicadas na análise de perfis. Como o modelo selecionado ajusta-se bem aos dados dielétricos disponíveis na literatura, torna-se atraente aplicá-lo à dados experimentais obtidos em rochas de campos brasileiros produtores de hidrocarbonetos para interpretação da EPT corrida em poços destes campos petrolíferos. / The electromagnetic propagation tool (EPT) provides the propagation time (T_pl) and the attenuation (A) of an electromagnetic wave propagating in a lossy medium. These EPT responses depend on the dielectric permittivity of the medium. There are several models and mixing equations concerning the dielectric permittivity of reservoir rocks that can be used in the interpretation of the high frequency tool. However, the mixing equations do not take into account the distribution and the geometry of the pore space, and these parameters are essential to obtaining dielectric responses approximating a true rock. A model based on the parameters described above was selected and this was applied to dielectric data available in the literature. A good agreement was reached between the theoretical curves and experimental data, confirming that the distribution and geometry of the pore space must be considered in the development of a realistic model. Aspect ratio distribution functions of the pores were also obtained, which were used for generating several curves relating the EPT responses to various oil/gas saturations. These curves were applied to the log analysis. The selected model fit the dielectric data available in the literature reasonably well, thus, making it suitable for application to experimental data of rock from Brazilian producing fields for the interpretation of the EPT in these fields.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-12-03T14:21:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ModelamentoPerfisInducao.pdf: 7044960 bytes, checksum: d5bab32aace9336f20b8ad536792da82 (MD5) / Rejected by Albirene Aires (albireneufpa@gmail.com), reason: Favor incluir as palavras chave on 2014-12-05T11:44:55Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-12-10T00:15:50Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPerfisInducao.pdf: 7044960 bytes, checksum: d5bab32aace9336f20b8ad536792da82 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-12-10T17:23:04Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPerfisInducao.pdf: 7044960 bytes, checksum: d5bab32aace9336f20b8ad536792da82 (MD5) / Made available in DSpace on 2014-12-10T17:23:04Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPerfisInducao.pdf: 7044960 bytes, checksum: d5bab32aace9336f20b8ad536792da82 (MD5) Previous issue date: 1990 / Neste trabalho foi implementado um algoritmo que permite o modelamento das respostas das sondas de indução (2C40-6FF40-ILD) em meios estratificados. O algoritmo é baseado na solução exata do problema de Green para a obtenção do campo eletromagnético em meios estratificados com a fonte no interior dos estratos. O código computacional desenvolvido neste trabalho para meios horizontalmente estratificados tem o objetivo de auxiliar o desenvolvimento de códigos mais avançados, uma vez que é difícil o acesso para a pesquisa dos códigos citados na literatura. Os códigos disponíveis são comerciais, mas estes têm o objetivo único de produção, além de serem, bastante dispendiosos. O código também foi utilizado para a obtenção da resposta da sonda ILD em reservatórios do tipo laminado, onde a sonda tem baixa resolução. E verificamos, a exemplo de outros autores, que a condutividade obtida pela sonda pode ser expressa pela soma das condutividades dos constituintes litológicos do reservatório ponderada pelos seus percentuais em volume. Esta relação foi definida para reservatórios ainda não analisados na literatura, tipo arenito moderadamente silicificado/arenito/folhelho, onde a lâmina de interesse tem resistividade superior a resistividade de uma das lâminas, e inferior a resistividade da outra. Foi analisado também reservatórios tipo arenito moderadamente silicificado/arenito, onde a lâmina de interesse tem a menor resistividade. Com esta relação pode-se obter para estes reservatórios laminados estimativas mais precisas de Rt e consequentemente determinar as saturações água/óleo e produtibilidade, também de forma mais precisa. / This thesis applies an algorithm that allows modeling induction sondes (2C40-6FF40-ILD) responses in multi layered configuration. The algorithm is based on the exact solution of the Green's problem to obtain the electromagnetic field in layered configuration with the source located in the layers. The code, developed for horizontal layers, can be used to develop advanced codes, because there are not codes available for research. The codes available are commercial, but these are very expensive and primarily useful for production processing. The code was applied to obtain ILD sonde response in laminated reservoirs, where the sonde resolution is reduced. Using this code, the conductivity from the sonde response, can be expressed by the sum of the conductivities of the lithologic constituints of the reservoirs, weighed by their volume contribution. This relation was defined for reservoir type not yet analyzed. The models are for reservoirs of silicified sandstone/sandstone/shale where the resistivity of the layer of interest is superior to the resistivity of one layer and inferior to the resistivity of another layer, and for reservoirs of silidfied sandstone/sandstone where the layer of interest has the lower resistivity. These models allow the establishment of better estimates of resistivity (Rt) in layered reservoirs, and consequently, to the hydrocarbon saturation and producibility.

Made available in DSpace on 2011-03-23T21:19:32Z (GMT). No. of bitstreams: 0 / Item created via OAI harvest from source: http://www.bdtd.ufpa.br/tde_oai/oai2.php on 2011-03-23T21:19:32Z (GMT). Item's OAI Record identifier: oai:bdtd.ufpa.br:233 / ANP - Agência Nacional do Petróleo / The Common-Reflection-Surface (CRS) stack is a new seismic processing method for simulating zero-offset (ZO) and common-offset (CO) sections. It is based on a second-order hyperbolic paraxial approximation of reflection traveltimes in the vicinity of a central ray. For ZO section simulation the central ray is a normal ray, while for CO section simulation the central ray is a finite-offset ray. In addition to the ZO section, the CRS stack method also provides estimates of wavefield kinematic attributes useful for solving interval velocity inversion, geometrical spreading calculation, Fresnel zone estimate, and also diffraction events simulation. In this work, Its proposed a new strategy to do a pre-stack depth migration by using the CRS derived wavefield kinematic attributes, so-called CRS based pre-stack depth migration (CRS-PSDM) method. The CRS-PSDM method uses the CRS results (ZO section and kinematic attributes) to construct an optimized stack traveltime surface along which the amplitudes of the multi-coverage seismic data are to be summed and the result is put in a point of the migration target zone in depth. In the same sense as in Kirchhoff type pre-stack depth migration (K-PSDM), the CRSPSDM method needs a migration velocity model. Unlike the K-PSDM method, the CRS-PSDM needs only to calculate the zero-offset traveltimes, i.e, along only ray conecting the considered point in depth to a given coincident position of source-receptor at surface. The final result is a zero-offset time-to-depth converted seismic image of reflectors from pre-stack seismic data. / O empilhamento por superfície de reflexão comum (ou empilhamento SRC), conhecido como empilhamento CRS, do inglês Commom reflection surface, constitui-se em um novo método para o processamento sísmico na simulação de seções afastamento nulo (AN) e afastamento comum (AC). Este método é baseado em uma aproximação paraxial hiperbólica de segunda ordem dos tempos de trânsito de reflexão na vizinhança de um raio central. Para a simulação de seção AN, o raio central é um raio normal, enquanto que para a simulação de uma seção AC o raio central é um raio de afastamento finito. Em adição à seção AN, o método de empilhamento SRC também fornece estimativas dos atributos cinemáticos do campo de onda, sendo aplicados, por exemplo, na determinação (por um processo de inversão) da velocidade intervalar, no cálculo do espalhamento geométrico, na estimativa da zona de Fresnel, e também na simulação de eventos de tempos de difrações, este último tendo uma grande importância para a migração pré-empilhamento. Neste trabalho é proposta uma nova estratégia para fazer uma migração em profundidade pré-empilhamento, que usa os atributos cinemáticos do campo de onda derivados do empilhamento SRC, conhecido por método CRS-PSDM, do inglês CRS based pre-stack depth migration. O método CRS-PSDM usa os resultados obtidos do método SRC, isto é, as seções dos atributos cinemáticos do campo de onda, para construir uma superfície de tempos de trânsito de empilhamento, ao longo da qual as amplitudes do dado sísmico de múltipla cobertura são somadas, sendo o resultado da soma atribuído a um dado ponto em profundidade, na zona alvo de migração que é definida por uma malha regular. Similarmente ao método convencional de migração tipo Kirchhoff (K-PSDM), o método CRS-PSDM precisa de um modelo de velocidade de migração. Contrário ao método K-PSDM, o método CRS-PSDM necessita apenas computar os tempos de trânsito afastamento nulo, ao seja, ao longo de um único raio ligando o ponto considerado em profundidade a uma dada posição de fonte e receptor coincidentes na superfície. O resultado final deste procedimento é uma imagem sísmica em profundidade dos refletores a partir do dado de múltipla cobertura.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-20T12:35:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AvaliacaoMetodosEletromagneticos.pdf: 5388098 bytes, checksum: e8b24df5c1300748048fd62b50531b18 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:10:21Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-22T12:41:38Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AvaliacaoMetodosEletromagneticos.pdf: 5388098 bytes, checksum: e8b24df5c1300748048fd62b50531b18 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T13:16:46Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AvaliacaoMetodosEletromagneticos.pdf: 5388098 bytes, checksum: e8b24df5c1300748048fd62b50531b18 (MD5) / Made available in DSpace on 2014-08-22T13:16:46Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AvaliacaoMetodosEletromagneticos.pdf: 5388098 bytes, checksum: e8b24df5c1300748048fd62b50531b18 (MD5) Previous issue date: 2006 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / Em investigações geofísicas rasas que empregam os métodos eletromagnéticos indutivos mais avançados, alvos com baixo número de indução (Low Induction Number – LIN) produzem anomalias eletromagnéticas muito baixas e de difícil interpretação. Para suprir esta deficiência, neste trabalho são estudados a aplicabilidade de campos eletromagnéticos polarizados e focalizados – POLFOCEM como fonte primária de indução. Os campos E.M. focalizados e polarizados, vertical e horizontalmente, são obtidos pelas combinações vetoriais de pares de dipolos transmissores e, ocorrem na região central entre eles. A focalização é observada nesta região na profundidade de 0,25 do espaçamento entre esses transmissores – L. Portanto, máximos acoplamentos podem ser obtidos através da seleção da polarização de acordo com a geometria do alvo, ocorrendo um aumento na densidade de fluxo magnético sobre ele e, máximas anomalias produzidas. É utilizada uma metodologia numérica para o cômputo dessas anomalias por meio da técnica dos elementos finitos para solução do problema 2,5-D. Em todos os experimentos numéricos são realizadas comparações qualitativas e quantitativas entre as respostas obtidas pelos sistemas POLFOCEM e convencional, o qual emprega um único dipolo como transmissor (dipolo-dipolo). As anomalias produzidas pelo sistema POLFOCEM, em que os dipolos transmissores são acionados simultaneamente, correspondem à soma das anomalias produzidas por cada um desses dipolos independentes, caracterizando, desta forma, a linearidade dos campos eletromagnéticos. Os experimentos numéricos são realizados para alvos prismáticos bidimensionais com três diferentes inclinações, inseridos num semi-espaço resistivo, e para as freqüências das fontes na faixa das ondas de rádio. As anomalias assimétricas no sistema convencional, que se tornam simétricas no sistema POLFOCEM, apresentam valores menores em amplitude. Contudo, aquelas anomalias tanto assimétricas quanto simétricas que se tornam anti-simétricas apresentam valores maiores. Em decorrência dessas diminuições e aumentos nas amplitudes ocorrem rotações nos diagramas de Argand, no sentido horário e anti-horário para alvos com baixos valores de condutividade, respectivamente. Em experimentos de identificação de presença de dois alvos próximos, o sistema convencional é capaz de identificá-los primeiramente, prevalecendo o seu uso. / In the near surface investigations employing advanced electromagnetic induction methods, the low induction number – LIN targets produce very low anomalies which are difficult to interpret. To overcome this difficulty, in this thesis I have studied the applicability of the Polarized Focused Electromagnetic – POLFOCEM fields as a primary inductive source. The horizontal end vertical E.M. field polarizations are obtained by the vectorial combinations of the fields due to the coil pairs in the central region between them; and the focalization occurs in this region at a depth 0.25 of the spacing between the transmitting coils – L. Therefore, maximum coupling can be obtained by selecting the field polarizations depending upon the geometry of the target and focusing the magnetic flux on it by varying L; consequently producing higher anomalies. This study has been realized by computing anomalies employing the Finite Element Technique for 2.5D problem. In all these numerical experiments a comparative qualitative and quantitative study has been done between the POLFOCEM anomalies and the conventional system anomalies where only one transmitter is used (dipole-dipole system). The POLFOCEM system anomalies, wherein both transmitters T1 and T2 are energized simultaneously, correspond to the sum of the anomalies produced by the T1-R and R- T2 dipoles independently; obeying in this way the linearity the electromagnetic fields. These numerical experiments have been carried out for a bidimensional prismatic model in a resistive half-space at three different inclinations and using radio frequencies. The asymmetrical anomalies obtained in a conventional coil system become symmetrical in POLFOCEM system, and have their amplitudes reduced in relation to that of the conventional system. On the other hand conventional anomalies, symmetrical or asymmetrical, become anti-symmetrical in POLFOCEM, and have their amplitude increased. As a consequence, the Argand diagrams of LIN targets for conventional coil system suffer respectively anticlockwise or clockwise rotation for POLFOCEM system. The experiments, to verify the quality in identifying the two nearby targets, shows that the conventional system is superior in relation to the POLFOCEM system.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-20T12:07:56Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Ausência de palavras-chave on 2014-08-07T12:42:09Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T15:30:42Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:32:11Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Made available in DSpace on 2014-09-18T11:32:11Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) Previous issue date: 2003 / O objetivo central deste trabalho é o estudo e a aplicação do método Kalman-Bucy no processo de deconvolução ao impulso e de deconvolução com predição, onde é considerado que os dados observados são classificados como não-estacionários. Os dados utilizados neste trabalho são sintéticos e, com isto, esta Tese tem características de um exercício numérico e investigativo. O operador de deconvolução ao impulso é obtido a partir da teoria de CRUMP (1974) fazendo uso das soluções das equações Wiener-Hopf apresentadas por KALMAN-BUCY (1961) nas formas contínuas e discretas considerando o processo como não estacionário. O operador de predição (KBCP) está baseado nas teorias de CRUMP (1974) e MENDEL ET AL (1979). Sua estrutura assemelha-se ao filtro Wiener-Hopf onde os coeficientes do operador (WHLP) são obtidos através da autocorrelação, e no caso (KBCP) são obtidos a partir da função bi(k). o problema é definido em duas etapas: a primeira consta da geração do sinal, e a segunda da sua avaliação. A deconvolução realizada aqui é classificada como estatística, e é um modelo fortemente baseado nas propriedades do sinal registrado e de sua representação. Os métodos foram aplicados apenas em dados sintéticos de seção fonte-comum obtida a partir dos modelos com interfaces contínuas e camadas homogêneas. / The main objective of this work is the study and the application of the Kalman-Bucy method in the processo f deconvolution to the impulse and deconvolution with prediction, considering the observed data as no stationary. The data used in this work are synthetic and, with this, this Thesis has characteristics of a numerical and search. The operator of deconvolution to the impulse is obtained from the Crump theory (1974), doing use of the solution of equation of Wiener-Holp presented by Kalman-Bucy in the continuoun and discrete forms considering the stacionary process. The prediction operator (KBCP) is based the Crump (1974) and Mendel et al (1979) theorics. Its structure resembles the Wiener-Hopf filter, where the coefficients of the operator are obtained through the autocorrelation, in the case (KBCP) are obtained from the function bi(k). The problem is defined in two steps: the first consists of the generation of the signal, and second of its evaluation. The deconvolution performed is classified as statistics, and is a model based in the properties of the registered signal and its representation. The method were applied only in synthetic data with common-shot section obtained from models with continuous interfaces and homogeneous layers.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-29T16:01:44Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ProcessamentoInterpretativoDados.pdf: 4112570 bytes, checksum: 6957010292e9aff837e775e8d78ca0ed (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-27T14:19:13Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ProcessamentoInterpretativoDados.pdf: 4112570 bytes, checksum: 6957010292e9aff837e775e8d78ca0ed (MD5) / Made available in DSpace on 2014-06-27T14:19:13Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ProcessamentoInterpretativoDados.pdf: 4112570 bytes, checksum: 6957010292e9aff837e775e8d78ca0ed (MD5) Previous issue date: 2007 / ANP - Agência Nacional do Petróleo / Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área. / We present two potential field data interpretation methods applied to hydrocarbon prospecting. The first one uses aeromagnetic data to estimate the horizontal projection of the limit between the continental and oceanic crusts. This method is based on the existence of geological sources, which are magnetized and belongs exclusively to the continental crust. In this way, the estimates of the sources’ ends are used as estimates of the continental crust limits. The total-field anomaly measured above the continental shelf, continental slope and part of the continental rise is amplified using the downward continuation operator implemented in two different ways: the equivalent layer principle, and the Dirichlet boundary condition. Most of the computational load in computing the downward continued anomaly comes from the solution of a large-scale system of linear equations. This computational effort has been reduced not only by processing the whole area by moving windows of smaller dimensions, but also by the use of the conjugate gradient in the solution of the system of equations. Because the downward continuation operator is unstable, it was stabilized through the first-order Tikhonov stabilizing functional. Tests with noisecorrupted synthetic data have shown the efficiency of both implementations to enhance the termination of magnetic sources belonging to the continental crust, allowing, in this way, the estimation of the limit between the continental and oceanic crusts. Both implementations have been applied to two different areas offshore the Brazilian coast: Foz do Amazonas and Bacia do Jequitinhonha. The second method simultaneously delineates the basement topography and the geometry of salt structures occurring within the sedimentary rocks using gravity data. The interpretation models consist of a set of vertical, 2D juxtaposed prisms for the sedimentary pack, and of 2D horizontal prisms with polygonal cross sections for the salt structures. The solution has been stabilized by incorporating the geometric characteristics of the basement relief and of the salt structures, which are compatible with the a priori knowledge about the geological setting. To this end, we imposed inequality constraints to the parameters of the interpretation model, and used the stabilizing functionals known as global smoothness, weighted smoothness, and mass concentration along selected directions. We applied the method to synthetic data produced by simulated intracratonic and marginal basins, presenting density contrast with the basement varying with depth and containing salt structures. The results have shown that the method is potentially useful in simultaneously delineating the faulted basement relief and the salt structures geometries. We applied the method to real data along two gravity profiles across the Campos and Jequitinhonha Basins, and obtained interpretations in accordance with the known geology of the area.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-02T15:50:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-30T13:51:25Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) / Made available in DSpace on 2014-06-30T13:51:25Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) Previous issue date: 1993 / PETROBRAS - Petróleo Brasileiro S.A. / A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P1, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm3 e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas. / The identification and description of lithological parameters of a formation are essential in the evaluation of complex formations. Based on this, the combination of the nuclear tool response in uncased wells has been used systematically. The resultant logs can be considered as the interaction between two distinct phases: • The radiation transport phase from a source to one or more detectors through the formation. • The detection phase that consists of the collection of radiation, its transformation into current pulses, and the spectral distribuition of these pulses. As the presence of the detector does not strongly affect the radiation transport result, each phase can be simulated independent of the other, which allows us to introduce a new type of model in which the transport phase and the detection phase are uncoupled. In this work, the final response is simulated combining transport numerical solutions with a library of the detector responses to different incident energies and for each specific source - detector array. The radiation is calculated by the Finite Elements Method (FEM), as a 2½-D scalar flux derived from the numerical solution of the multigroup diffusion approximation of the Boltzmann transport equation in phase space. This is known as P1 approximation, where the variable direction is expanded in terms of the Legendre orthogonal polymonials, leading to the dimensionality reduction of the problem in such a way as to let it be more consistent with the FEM, where the flux depends only on the spatial variable and the physical properties of the formation. The NaI(Tl) response function is obtained separately by the Monte Carlo method (MC) where the life of a particle within the scintillator crystal is reconstructed simulating interaction by interaction the position, direction and energy of the different particles using a random number technique with associated appropriate probabilities laws. Each type of interaction (e.g., Rayleigh, Photo-electric, Compton and Pair production) is determined similarly and the simulation is concluded when the detector response functions are convolved with the scalar flux. The final response is the pulse-height spectrum of the simulated system. From this spectrum, a set of channels called detection windows are then selected. The count rates in each window show different dependencies on density and lithology. This fact allows one to use a combination of these windows in determining the density and photoelectric absorption factor of the formation. According to the method developed in this work, the logs in both thin and thick layers can be simulated. The performance of the method has been tested in complex formations, mainly where the presence of clay minerals, feldspars and micas have produced considerable effects sufficient to perturb the final response of the sonde. The results show that it is possible to identify physical and lithological parameters in formations having densities between 1.8 and 4.0 g/cm3 and photoelectric absorption factors in the interval of 1.5 to 5.0 barns/e-. The concentrations of Potassium, Uranium and Thorium can be obtained through the introduction of a new system of calibration which corrects the effects due to high variances and negative correlations observed on the mass concentration of Uranium and Potassium. In the simulation of the CNL response using the Tittle polynomial regression algorithm, it is verified that due to the limited vertical resolution of this sonde, the porosity value is poorly measured for most layers of thickness less than the source - far detector spacing, thus it has application only in thick layers. A new method was developed to solve this problem; the contribution of the relative area of each layer within the maximum information zone is determined. Thus, this neutron porosity makes possible an in-depth evaluation of expected CNL porosity-lithology response, convolving that area factor with the local formation porosity index, considering only thick layers. The presence of perturbating minerals is solved by considering the formation as formed by a predominant base matrix mineral, totally saturated by fresh water; the rest of the components are then considered as a perturbation of this base case. These results enable the calculation of synthetic well logs that can be used in inversion schemes in order to get a more detailed quantitative evaluation of complex formations.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T12:53:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-07-03T14:41:41Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) / Made available in DSpace on 2014-07-03T14:41:41Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) Previous issue date: 2013 / UFPA - Universidade Federal do Pará / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / INCT - Instituto Nacional de Ciência e Tecnologia / PETROBRAS - Petróleo Brasileiro S.A. / A utilização dos métodos de reflexão sísmica na exploração e desenvolvimento de reservatórios de hidrocarbonetos ocorre devido à sua vasta e densa amostragem, tanto em área quanto em profundidade, aliada ao refinamento de técnicas para o tratamento dos dados de reflexão sísmica, a partir destes dados, são geradas seções sísmicas, que após a aplicação de tratamento adequado, são utilizadas na interpretação dos estratos e/ou estruturas geológicas da subsuperfície. Neste trabalho é feita uma análise Geofísica Geológica de duas linhas sísmicas reais 2D marinhas da porção de quebra de talude da Bacia do Jequitinhonha. Para tanto, foi realizado um conjunto de processamento sísmico com objetivo de atenuar as reflexões múltiplas comuns em dados marinhos, além disso, foram estimados os modelos de velocidade em profundidade, utilizados para determinação das seções sísmicas migradas em profundidade. Nestas foram identificadas as superfícies refletoras. Através da análise dessas superfícies foram feitas as marcações de sismofácies, com base nos conceitos iniciais da sismoestratigrafia, com a finalidade de avaliar a qualidade do produto derivado do processamento sísmico, empregado neste estudo, para uma interpretação sismoestratigráfica, a qual está fundamentada na análise dos padrões de terminações dos refletores e padrão interno das sismofácies. / The seismic reflection methods use at the exploration and development hydrocarbon reservoirs is due to its vast and dense sampling in both area and depth, combined with techniques refinement to processing seismic reflection data, from these data, are generated seismic sections, which after applying appropriate treatment, are used in the interpretation of the strata and/or the subsurface geologic structures. At this paper have made one Geophysics Geological analysis at two real marine seismic lines 2D of break portion embankment Jequitinhonha Basin. Therefore, we performed a set of seismic processing in order to mitigate the multiple reflections common marine data, moreover, were estimated velocity models in depth, used to determine the depth migrated seismic sections. At these were identified reflective surfaces. Through analysis at these surfaces markings of sismofacies. have made based on the initial concepts seismic stratigraphy in order to evaluate the product quality derived from the seismic processing used at this study for seismic stratigraphic interpretation, which is based on the analysis termination patterns of reflectors and internal standard of sismofacies.