Programa de Pós-Graduação em Geofísica - CPGF/IG

URI Permanente desta comunidadehttps://repositorio.ufpa.br/handle/2011/2355

O Programa de Pós-Graduação em Geofísica da UFPA (CPGF) do Instituto de Geociências (IG) da Universidade Federal do Pará (UFPA). Foi o segundo no Brasil a formar recursos humanos em Geofísica em nível de pós-graduação stricto sensu. Criado em 1972, funcionou até 1992 junto com os Cursos de Pós-Graduação em Geoquímica e Geologia.

Navegar

Agora exibindo 1 - 20 de 42

Acesso aberto (Open Access)
Uma análise de métodos de supressão de múltiplas de superfície livre aplicados a um dado real
(Universidade Federal do Pará, 2011-10-14) OLIVEIRA, Andrei Gomes de; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
The seismic record in 2D linear events are common features in the seismograms. These events are overlap with interesting reflections (early reflections) complicating the stages of seismic processing and identification of possible subsurface geological horizons. It possible highlight multiple surface that distinguished from others by a low frequency, high amplitude and low velocity of propagation. The identification and removal of multiple reflflections in seismic data is a real step of great importance in seismic processing, because the attenuation can prevent many of these errors in the interpretation of geological images seismic. This study is the determination of a processing flow for multiple attenuation free surface. Methods of elimination of multiple surface (SRME, Radon filtering, predictive deconvolution domain τ - p and fk filltering) applied to a real data 2D were used in different combinations. These are analyzed and presented the results are the best combinations.
Acesso aberto (Open Access)
Análise de um esquema de inversão para estimativa de anisotropia local a partir de dados de VSP
(Universidade Federal do Pará, 2012-05-18) BARRETO, Adriano César Rodrigues; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
Inversion scheme for the determination of arbitrary in a vicinity of a receiver situated inside a borehole from the data obtained during a multi-offset multi-azimuthal VSP experiment is studied. The data consist of vertical components of the slowness vector and all components of the polarization vector of P waves direct and reected. The inversion is made from a first order approximation around an reference arbitrary isotropic medium, that does the relation between the data of slowness and polarization of P waves parameters and the weakly anisotropic medium (WA parameters). This inversion scheme applies to arbitrary anisotropy medium, provided that the anisotropy is weak and independent of downhole type (vertical, horizontal or inclined). The results show that reversing is sensitive to the number and orientation of the proles of sources that are distributed on the surface. An analysis of the design of an experiment and discussion about the choice of parameters of the reference isotropic medium are presented. Within the realistic assumptions about the level of noise is possible to estimate the phase velocity for a limited angular aperture around the well.
Acesso aberto (Open Access)
Aplicação de Well-Tie no imageamento de uma área da bacia Tumbes-progresso
(Universidade Federal do Pará, 2014-08-29) CHAGAS, Tiago Paulo Câncio das; SILVA, Carolina Barros da; http://lattes.cnpq.br/5306784916926352; GOMES, Ellen de Nazaré Souza; http://lattes.cnpq.br/1134403454849660
Geological areas with large lateral velocity variations, such as in areas with faults and salt domes are complex areas, in general, seismic imaging time is limited. In these cases, the seismic imaging depth can provide better results. However, a velocity model is needed in depth with precision, which is not an easy task. An alternative to estimate the depth velocity model is provided by Well-tie technique. With the calibration of the synthetic trace data obtained from wells, seismic section with the depth map of velocity can be obtained through the time-depth curve. In this study, we used the Well-Tie to improve the model of RMS velocity obtained from the processing of seismic data. The updated velocity model was converted to depth and used to generate an in depth migrated seismic image of a region of Tumbes Basin-Progress.
Acesso aberto (Open Access)
Aproximações não-hiperbólicas do tempo de trânsito utilizando aproximantes de Padé
(Universidade Federal do Pará, 2017-07-03) NEVES, Rodolfo André Cardoso; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
Traveltime approximation is a fundamental tool of the stack and migration steps in seismic data processing. To increase the accuracy of these approximations, we propose new traveltime approximations based on Padé approximants, to CMP and CRS gathers. Hyperbolic approximations such as normal moveout (NMO) and comom reflection surface (CRS) are taylor series approximations of second order of the reflection traveltime. Padé approximants appear as an alternative to Taylor series, because they converge quickly to the desired function, and they have a major radius of convergence improving approximations acuracy. They can be obtained through the proper Taylor serie of the approximated function. This new approximation is obtained from the [2/2] Padé approximation of the generalized moveout equation; and from [2/2] Padé approximation of the Taylor series expansions of fourth order of the CRS surface. The acuracy of Padé approximation is superior when compared with other convencional approximations: normal moveout, shifted hyperbola and Transversal isotropic medium with vertical symetry axis (VTI). CMP gather Padé approximations depend just only one more parameter than normal moveout approximation and they keep the acuracy for long offsets. CRS gather non hyperbolic approximations, non hyperbolic CRS, fourth order CRS and Padé CRS, have major acuracy than hyperbolic CRS, increasing the convergence of the approximation for offset and CMP domain. The quadratic approximation of fourth order CRS is superior than non hyperbolic CRS approximation, producing less error in least square CRS parameter inversion.
Acesso aberto (Open Access)
Avaliação do modelo hidrodinâmico SMC: uma aplicação na Costa Norte Brasileira
(Universidade Federal do Pará, 2016-05-24) MOTA, Fabíola Cardoso da; ROLLNIC, Marcelo; http://lattes.cnpq.br/6585442266149471
Coastal environments are highly dynamic because of complex processes that occur in the transition between the ocean and the continent. Climate changes induce changes in the coastal zone that cause decline in economic, environmental and subsistence values. The integrated coastal zone management is a continuous and dynamic process, in which decisions should be taken in order to sustainable development and protection of resources of coastal and marine areas. The main goal of integrated coastal zone management is to improve the life quality of human communities that are dependent of coastal resources, taking into account the maintenance of biological diversity and ecosystems productivity. The coastal modelling system called SMC Brazil was introduced in Brazil by the Ministry of Environment in order to standardize the modeling throughout the coastal zone of the country. This study aims to implement the model SMCBrasil in North Coast of Brazil (Princesa beach- Algodoal Island), evaluate it and analyze the hydrodynamic answers, in order to know more about coastal processes that govern this region. In order to understand the processes in the study area, the methodology consisted in collect in situ data to know about the intensity of coastal parameters. 6 cases were created in the model to represent the wave theories (Linear, Stokes and Composite) and sediment transport (Bailard and Soulsby) for two moments of tide: low tide and high tide. Two mesh grids were created in the study area, and one entry point that would fit better to it and the conditions limits of the model was chosen. The propagated waves by the composed theory represented the study area in question more than the others, and the simulated sediment transport simulated by Soulsby theory suffered less interference of interpolation than the Bailard, since it considers the bed transport and suspension. The simulated cases in this study were based on boundary conditions that more often occur on the study area, so the results can be used as a basin to understanding how the coastal processes occur in it in general way.
Acesso aberto (Open Access)
Cálculo da porosidade: identificação do argilomineral
(Universidade Federal do Pará, 2017-04-20) ALMEIDA, Thales Luiz Pinheiro de; ANDRADE, André José Neves; http://lattes.cnpq.br/8388930487104926
In the daily practice of Formation Evaluation it is common the adoption of simplifying hypotheses or premises about the physical properties of the reservoir rock constituent materials to allow the porosity calculation. The knowledge of the physical properties of the clay in the reservoir rock is fundamental for porosity calculation. In this work it is argued that the physical properties of the clay mineral in the reservoir rock are different from the mean physical properties of the nearby shale layers. Geologically, the clay is one of the constituent materials of shale rock and to admit that the physical properties of the shale are equal to the physical properties of the clay in the reservoir rock means to disregard all the other constituents and to assume a sedimentary continuity that due to numerous postpositional processes may not occur. In this work, we apply the angular competitive neural network to the Density- Neutron Plot to show that if a reservoir rock and a shale present in the basin have the same clay, they have the same angular pattern. This methodology is presented with synthetic data and evaluated with actual well logs and core analysis from borehole drilled in the Namorado’s field, in the Campos Basin, Brazil.
Acesso aberto (Open Access)
Cálculo de porosidade com a rede neural competitiva
(Universidade Federal do Pará, 2015-10-26) ROSELLÓN GUZMAN, Laura Yesenia; ANDRADE, André José Neves; http://lattes.cnpq.br/8388930487104926
Porosity is the petrophysical property that quantifies the fluid volume in the reservoir rock under for subsurface original condition. However, its calculation by the densityneutron method is extremely difficult in non cored borehole by the lack of the knowledge about the matrix physical properties (density and neutron porosity). This work presents a method for enabling the use of density-neutron Method in non cored boreholes, showing a realistic estimate of the matrix physical properties for each reservoir layer, using a angular competitive neural network. For each layer, network training is performed in the density-neutron plot built with the points of this layer and the information about the grain density (matrix density), obtained in the core analysis. This method is presented with synthetic data, which satisfy the petrophysical model and real data from two cored wells in the Namorado field, Campos basin.
Acesso aberto (Open Access)
Conversão tempo-profundidade de seções sísmicas empilhadas por raio imagem e raio normal
(Universidade Federal do Pará, 2014-06-24) REZENDE, Diogo Pena; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
In general, the oil industry makes the time-to-depth conversion of seismic data by the image ray tracing method. This method takes time into the depth, point to point, the amplitudes of the time migrated seismic section. For each point of the migrated time section, it is necessary to trace a ray perpendicular to the surface. After this, the amplitude of the migrated point of the section takes place in depth. The seismic migration method pre- or post-stack consists of placing seismic events in the correct positions in time or depth sections. Seismic depth sections provide an image near of the subsurface, in order to facilitate the identification of possible oil accumulating geological structures. The conversion of sections from the time to the depth domain is an intermediate step in the construction of seismic images in depth. This work developed and tested a method of converting time to depth the zero-offset seismic sections. In this case, the construction of sections in depth uses normal ray tracing method. The proposed method makes use of the (slowness versus time of intersection) transformation on the zero-offset section. Each point in the domain provides initial conditions for the normal ray tracing: a start position of the initial rays and initial angles formed with the normal to the surface, i.e., the slowness initial parameter. Unlike ray image method, several rays use the same travel time and the same initial position, defining an isochronous curve. The amplitude of each point in the zero-offset section takes place to depth from the distribution of values along each isochronous curve in depth. The image ray based time-to-depth conversion has good recovery of the depths of reflectors as well as low computational cost, since it is necessary only one ray to convert each point of the section in time. However, the reflector continuity may be damaged in case of sharpened curvature. In turn, by normal ray the time-to-depth conversion correctly approximate the depth of the reflectors, since the same point in time assigns several times in depth. However, it has a higher computational cost, because it is necessary many rays to convert one point in time.
Acesso aberto (Open Access)
Estimativas de umidade no solo usando GPR
(Universidade Federal do Pará, 2011) MESQUITA, Marcelo Jorge Luz; LUIZ, José Gouvêa; http://lattes.cnpq.br/8676530374310847
This work aims to analyze the effectiveness of GPR (Ground Penetrating Radar) as a measure of moisture in comparison to the established system TDR (Time Domain Reflectometry) and introduce changes in the calibration equations that allow estimation of the moisture from the wave velocity of the GPR. This work presents the methodologies for collecting and processing the GPR data, aimed at determining the moisture. One of the methods was applied in Cuiarana (Salinópolis - PA) with 400 MHz antenna and compared with data from TDR. The relative permittivity obtained by means of the velocities of the waves of GPR were used to infer the values of moisture from Topp and Roth's equations and thus to compare them to the values obtained with TDR. The data were statistically correlated and analyzed, enabling the establishment of a new calibration equation as a function of relative permittivity and other in terms of the speed of the GPR wave. The study successfully showed the feasibility and limitations of estimating water content with GPR.
Acesso aberto (Open Access)
Exemplo de bom condicionamento inconveniente causado numericamente na inversão gravimétrica para a estimação das densidades de uma camada
(Universidade Federal do Pará, 2014) SOARES, William Pareschi; MARTINS, Cristiano Mendel; http://lattes.cnpq.br/8303640454649778
A well-conditioned sensitivity matrix can be inconvenient for estimating densities of a layer? We found an example of this inconvenience numerically caused in the gravimetric inversion when the horizontal dimensions of the elementary sources that make up the interpretive model are very small. Amazingly in this case a gravity inversion to obtain the equivalent layer is not effective and does not fit the data. We found that in this case the well-conditioning of the sensitivity matrix occurs simultaneously with low singular values. This feature of the sensitivity matrix leads to severe loss of resolution and leads to biased estimates and very smooth. It happens that part of the resolution would be mathematically possible to be obtained is lost due to this phenomenon numerical computational degradation of the sensitivity matrix . We present a procedure for repayment of the resolution for mapping the density distribution of a layer, which enables new perspectives on gravimetric applications, including environmental studies. We skirted the numerical problem with a semi-heuristic approach which extends the horizontal dimensions of the elementary sources and subsequently corrects the estimates. We obtained with this new procedure in synthetic tests the distribution of the density contrast outlining lateral contacts between regions of different density contrasts, which would only be possible to recover elemental sources larger. We apply this methodology to the dataset of the Thomas Farm landfill site landfill.
Acesso aberto (Open Access)
Experimental verification of Hudson and Eshelby-Chen’s effective crack theory
(Universidade Federal do Pará, 2015-12-04) HENRIQUES, Jéssica Pena; FIGUEIREDO, José Jadsom Sampaio de; http://lattes.cnpq.br/1610827269025210
Physical modeling of cracked / fractured media in downscaled laboratory experiments has served as a great alternative for understanding the anisotropic media behavior. In this work, it was performed ultrasonic measurements on samples with low crack densities and different aspect ratios. The main goal was to investigate the experimental behavior of elastic parameters, such as: waves velocities and Thomsen parameters ε and 𝛾 and elastic stiffness coefficients for transverse isotropic media. Comparison of the results with the predictions made by the effective models of Hudson (1981) and Eshelby-Cheng (1993), it was also investigated in this work. Twelve samples were prepared with two types of cracks density, 5 and 8%. The cracks that have three different aspect ratios (0.133, 0.1778 and 0.2667) were formed by rubber inclusions in a homogeneous isotropic matrix resin. Moreover, an isotropic matrix sample was constructed by only epoxy resin. Among all samples, six (three for each density) have only one aspect ratio type (samples with single crack), while another six (three for each density) have three types of different aspect ratio (mixed samples). Among the predictions of the models, the Eshelby-Cheng (1993) shows a better fit in relation to the experimental results for samples with single crack (for the two densities of inclusions). However, none of the models predicts accurately with minimal tendency for the mixed samples.
Acesso aberto (Open Access)
High resolution GPR applied to the “Capela Pombo”
(Universidade Federal do Pará, 2016-02-29) ESPEJO ZEBALLOS, Manuel Alejandro Junior; SILVA, Marcos Welby Correa; http://lattes.cnpq.br/3213216758254128
In this work was applied electromagnetic method (GPR), to identify surface layers in a speci c urban environment. The measurements were performed to map the surface layers and artifacts that make up the study area "Capela de Nosso Senhor dos Passos - a Capela Pombo" in the Barrio of Campina-Bel em. This chapel has a previous work that rescues the importance y architectural features that it has, together with the evidence recollected by the author Dominic Savio de Castro Oliveira, who noted the historical importance that has not only to be the last private chapel in the city of Belém, this author also acknowledges her design as an Italian architect's work possible Guissepi Jos e Antonio Landi. The pro les obtained after corresponding processing procedure, show that this probably presences of structures, rest of the previous building structure, which previously was in this space. The study aimed to identify anomalies. The study aimed to identify structural abnormalities, as well as possible characteristic burials of the time when it was actively used. The results of the GPR measurements were encouraging, since the method presented response of approximately 1.80 meters deep, one can identify anomalies control for the reasons as presented Chapel, and present possible anomalies related to structural foreign bodies.
Acesso aberto (Open Access)
Identificação de fácies em perfis de poço com algoritmo inteligente
(Universidade Federal do Pará, 2014-02-17) SANTOS, Renata de Sena; ANDRADE, André José Neves; http://lattes.cnpq.br/8388930487104926
Facies identification in an uncored borehole is a classic problem in formation evaluation. In this study, this problem is treated as the extraction of geological information or facies descriptions from a cored borehole in terms of their physical properties registered in well logs and perform their encoding through the parameters L and K calculated from porosity logs, and shaliness calculated using the natural gamma ray log to construct the Vsh-L-K plot. For interpretation is presented an intelligent algorithm based on the competitive generalized angular neural network, built for angular pattern classification or data clustering in ndimensional space that have an approximately ellipsoidal envelope, which are the characteristics of clusters in the Vsh-L-K plot and make your visual interpretation extremely complex. The application of intelligent algorithm is able to identify and classify the layers present in uncored boreholes, in terms of the facies identified in the cored borehole or in terms of its main mineral, when it is absent in the cored borehole. This methodology is presented with synthetic data and well logs from cored boreholes in Namorado oil field, in the Campos Basin, located on the continental shelf of Rio de Janeiro, Brazil.
Acesso aberto (Open Access)
Identificação litológica com affinity propagation
(Universidade Federal do Pará, 2013-03-26) CALDAS, Nayara Safira da Silva; GUERRA, Carlos Eduardo; http://lattes.cnpq.br/7633019987920516; ANDRADE, André José Neves; http://lattes.cnpq.br/8388930487104926
This work presents a methodology for solving the problem of extracting geological information, such as identification of lithologies at depth, directly from wireline logs. For this propose, the Vsh-M-N plot is used as the formation evaluation technique for identifying the lithologies in the logged borehole, in terms of the physical properties of the main mineral in each lithology. However, the visual interpretation of this graphic is limited by the large spread of the points in the graph. To evaluate a form of computational interpretation of the Vsh-M-N, the affinity propagation algorithm is used in reason of its characteristic transport of information among similar objects, which enables the interpretation of Vsh-M-N plot considering the physical properties and continuity in the depth. This methodology is presented using synthetic data and well logs from one borehole in the Namorado field. Campos basin. Brazil.
Acesso aberto (Open Access)
Inversão de velocidades por otimização global usando a aproximação superfície de reflexão comum com afastamento finito
(Universidade Federal do Pará, 2016-08-25) MESQUITA, Marcelo Jorge Luz; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
The recent geophysical literature has shown the building of an accurate initial model is the more appropriate way to reduce the ill-posedness of the Full Waveform Inversion, providing the necessary convergence of the misfit function toward the global minimum. Optimized models are useful as initial guess for more sophisticated velocity inversion and migration methods. I developed an automatic P-wave velocity inversion methodology using pre-stack two-dimensional seismic data. The proposed inversion strategy is fully automatic, based on the semblance measurements and guided by the paraxial traveltime approximation, so-called Finite-Offset Common-Reflection-Surface. It is performed in two steps, at first using image rays and an a priori known initial velocity model we determine the reflector interfaces in depth from time migrated section. The generated depth macro-model is used as input at the second step, where the parametrization of the velocity model is made layer by layer. Each layer is separated from each other by smoothed interfaces. The inversion strategy is based on the scan of semblance measurements in each common-midpoint gather guided by the Finite-Offset Common-Reflection-Surface traveltime paraxial approximations. For beginning the inversion in the second step, the finite-offset common-midpoint central rays is built by ray tracing from the velocity macro-model obtained in the first step. By using the arithmetic mean of total semblance calculated from the whole common-midpoint gathers as objective function, layer after layer, a global optimization method called Very Fast Simulated Annealing algorithm is applied in order to obtain the convergence of the objective function toward the global maximum. By applying to synthetic and real data, I showed the robustness of the inversion algorithm for yielding an optimized P-wave velocity macro-model from pre-stack seismic data.
Acesso aberto (Open Access)
Investigação geofísica da contaminação por derivados de hidrocarbonetos nas redondezas do Posto Cristal (Belém-Pa)
(Universidade Federal do Pará, 2013) SILVA, Márcio Antonio dos Santos; SILVA, Lúcia Maria da Costa e; http://lattes.cnpq.br/8612431024609774
Acesso aberto (Open Access)
Investigações ambiental e forense com os métodos geofísicos radar de penetração do solo, polarização induzida e eletrorresistividade no Cemitério do Tapanã, Belém/Pará
(Universidade Federal do Pará, 2010-04-30) CATETE, Clístenes Pamplona; SILVA, Lúcia Maria da Costa e; http://lattes.cnpq.br/8612431024609774
The cemeteries are considered potential sources of underground water contamination. Since the use of geophysics in the mapping of the hydrostatic level (HL) has become a very important tool in the environmental monitoring, it makes possible to verify adequate or inadequate areas for burials in cemeteries. In the other hand, the localization of underground clandestine hollows and tunnels, can be of utility for the Forensic Science. In this work, geophysical methods were applied for mapping the HL in a cemetery and to find buried structures related to the Forensic Science. Ground Penetrating Radar (GPR) was used in the mapping of the HL in the whole Tapanã cemetery, while the resistivity and induced polarization (IP) methods were applied on a test site built in that cemetery for forensic, environmental and resgate related geophysical studies. The results with the GPR showed that the HL depth in the Tapanã cemetery varies from 7 m, at high topographic levels, to 1 m, at low topographic levels. The IP was sensible to the presence of a body buried in a shallow hollow as the clandestine hollows are and showed only fair results for a tunnel. The resistivity was effective for both the buried body and the tunnel, with better definition of the tunnel. The methods gave good results even within the weathering layer rich in clay that normally occurs in the Amazon region.
Acesso aberto (Open Access)
Modelagem do método MT 2D usando elementos finitos isoparamétricos
(Universidade Federal do Pará, 2017-04-07) SOUZA JUNIOR, Ivaldevingles Rodrigues de; SILVA, Marcos Welby Correa; http://lattes.cnpq.br/3213216758254128
The numerical modeling of the magnetotelluric (MT) method is of great importance for geophysics, since this method can be applied for several purposes, for example it can be used to study the crust of the earth as well as it can contribute to oil and gas exploration. For this reason, several researches have been carried out in recent decades to further develop MT. One of the major advantages of magnetotelluric is its relatively easy modeling, because in this technique the source is a plane wave. Existem varias ferramentas numericas que podem ser usadas para modelar o MT, There are several numerical tools that can be used to model the MT, among which stands out the nite element method (FEM). In this work the eciency of the isoparametric MEF for 2D modeling of the MT will be tested, whose main characteristic is the accomplishment of a transformation, with the purpose of mapping a distorted element to a coordinate system (natural coordinates ξ and η) where it becomes regular.
Acesso aberto (Open Access)
Modelagem numérica bidimensional de dados CSAMT com fonte dipolar elétrica usando elementos finitos
(Universidade Federal do Pará, 2016-08-20) PEREZ, Jarol David Garcia; SILVA, Valdelírio da Silva e; http://lattes.cnpq.br/9701131346395521; RÉGIS, Cícero Roberto Teixeira; http://lattes.cnpq.br/7340569532034401
This work presents the application of a 2,5D numerical modelling scheme to simulate data of the CSAMT (Controlled Source Audio MagnetoTellurics) electromagnetic method using the Finite Element method. We have applied the technique of separating the total electromagnetic field into two parts: the so called primary and secondary fields. The set of equations that govern the system has to be expressed in the spectral domain of the Fourier Transform. The solution is obtained as a composition of 2D solutions. In the first part of the dissertation, we present a validation of the answer generated for primary fields and a comparison of the 1D responses with those from the MT method. The second part of the text presents a validation of the 2,5D responses, as well as an analysis of the results, in terms of apparent resistivity and the measured electrical field related to an isolated body in a homogeneous half-space. We analyze the changes in the results resulting from changes in model parameters like the contrast of resistividade between the half-space and the anomalous body, the depth of the body and the frequency. The code developed is able to simulate with good precision the measurements of the CSAMT method in any survey configuration.
Acesso aberto (Open Access)
Modelagem sísmica acústica e elástica por diferenças finitas e imageamento do depósito de minério de ferro N4WS no estado do Pará
(Universidade Federal do Pará, 2017-07-03) HOLANDA, Rafael Mansano; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
The finite difference method consists in the numerical resolution of differential equations and one of its applications is a solution of the wave equations or elastodynamics, by replacing the partial derivatives in space and time by their approximations of finite differences. The present work consisted in the accomplishment of the seismic modeling using the method of the finite differences, to acoustic and elastic case. Then, we performed the seismic imaging, to acoustic case, in a model that simulates an iron ore deposit area of N4WS, located in the mineral province of Carajás, in the state of Pará. We found a heterogeneous, isotropic and bidimensional velocity model to acoustic and elastic cases. In order to generate the synthetic seismogram, we executed the program FDSKALAR, to acoustic case, and SOFI2D, to elastic case, using fourth-order approximation of the equation of acoustic and elastic wave in space and time. For a validation of the programs of finite differences and aid in the interpretation of the events, we compared the results with the travel times obtained by the ray theory. After the seismic modeling, we accomplished the processing of the data generate, using the open-source software package Seismic Unix, to obtain an image of the reflector contained in the mineral model. In this way, we searched a better understanding and control on the problems of seismic modeling and imaging, contributing to an interpretation of seismic data and understanding about a seismic wave propagation in environments where there are significant heterogeneities. In addition, we sought to demonstrate the application of seismic methods in the study and delimitation of mineral bodies.