Teses em Geofísica (Doutorado) - CPGF/IG

URI Permanente para esta coleçãohttps://repositorio.ufpa.br/handle/2011/2357

O Doutorado Acadêmico pertente a o Programa de Pós-Graduação em Geofísica (CPGF) do Instituto de Geociências (IG) da Universidade Federal do Pará (UFPA).

Navegar

Agora exibindo 1 - 12 de 12

Acesso aberto (Open Access)
Espalhamento elástico em meios anisotrópicos estratificados
(Universidade Federal do Pará, 1994-06-09) PROTÁZIO, João dos Santos; SCHOENBERG, Michael
AVO analysis is an importante tool for extracting lithological information from seismic data using the contrast in acoustic impedance at the lithological boundaries. The isotropic assumption behind this analysis does not hold in many cases. The advent of large offset surveys and multi-component data has revealed the presence of subsurface anisotropy. To interpret such data, the AVO analysis must include anisotropy. This work presents an AVO theory and numerical results for an anisotropic layered medium. This thesis contains three contributions. First, a new approach to the study of reflection-transmission at a plane interface between anisotropic media with a horizontal mirror plane is presented. The Zoeppritz equations are generalized to include anisotropy by the introduction of impedance matrices which greatly simplifies the previous formalisms. Second, the study of the P-wave reflection at an interface between an isotropic and a transversally isotropic medium is described and it is show that the reflected P-wave does not have information about the underlying anisotropic subspace for pre-critical incidence. Finally, the behavior of post-critical reflected and transmitted pulses through a stack of anisotropic layers is discussed. The post-critical pulses are show to carry valuable information on the anisotropy of the structure through which the waves propagate.
Acesso aberto (Open Access)
Estimativa de parâmetros elásticos em meios anisotrópicos
(Universidade Federal do Pará, 2003-06-20) GOMES, Ellen de Nazaré Souza; PROTÁZIO, João dos Santos; http://lattes.cnpq.br/4210442535067685
Amplitude, polarization and the slowness vector measurements carry information about the medium where wave propagation occurs. This thesis investigates these data aiming at the recovery of elastic properties in anisotropic media. Reflection coefficients can be estimated from amplitude data and depend nonlinearly on elastic and density contrasts across an interface. When the impedance contrast is weak, the linear approximations for the qP reflectivity are more convenient for inversion of density and elastic parameters using analysis of amplitude versus the angle of incidence (AVO) and amplitude versus the direction of the incidence plane (AVD). Partitioning the linear system defined by Zoepprittz equations allows one to write the solution of these equations in terms of impedance and polarization matrices. Using this solution, linear approximations for the qP reflectivity are derived for weak impedance contrasts and arbitrary symmetry classes of anisotropy. The linear approximations are evaluated for different acquisition geometries and choice of the reference medium. The approximations for the reflection coefficients of the reflected qP and the converted waves are in good agreement with the exact solution for incidence angles up to 30° for media that satisfy the weak impedance assumption. If a single oriented set of fractures is represented by a transversely isotropic effective medium, the linear approximations for qP reflectivity can be used to estimate the fractures orientation. Under these assumptions this problem is reframed as the estimation of the symmetry axis orientation from qP reflectivity data. This work shows the requirement of multiple components and multiple azimuthal data and quantifies the minimum amount of data for stable estimation. Also it is shown that the reflection coefficients of converted waves qS and qT only are sensitive to fractures dip. The inversion of polarization and slowness from multiazimutal VSP data are investigated for the estimation of local anisotropy. We use measurements of the vertical component of the slowness vector and the qP polarization data of direct and reflected waves. The inversion algorithm is validated in synthetic data sets for different choices of the wave front normal, reference medium and acquisition geometries. This analysis shows that only a subset of elastic parameters is recovered. An important application of this approach is its potential to determine the class of anisotropy. The application of this methodology to the Java Sea data set shows that isotropy and transversely isotropic models are inadequate to fit the data.
Acesso aberto (Open Access)
Interpolação de dados de campo potencial através da camada equivalente
(Universidade Federal do Pará, 1992-09-15) MENDONÇA, Carlos Alberto; SILVA, João Batista Corrêa da; http://lattes.cnpq.br/1870725463184491
The equivalent layer technique is an useful tool to incorporate (in the process of interpolation of potential field data) the constraint that the anomaly is a harmonic function. However, this technique can be applied only in surveys with small number of data points because it demands the solution of a least-squares problem involving a linear system whose order is the number of data. In order to make feasible the application of the equivalent layer technique to surveys with large data sets we developed the concept of equivalent data and the EGTG method. Basically, the equivalent data principle consists in selecting a subset of the data such that the least-squares fitting obtained using only this selected subset will also fit all the remaining data within a threshold value. The selected data will be called equivalent data and the remaining data, redundant data. This is equivalent to splitting the original linear systems in two sub-systems. The first one related with the equivalent data and, the second one, with the redundant data in such way that, the least-squares solution obtained by the first one, will reproduce all the redundant data. This procedure enables fitting all the measured data using only the equivalent data (and not the entire data set) reducing, in this way, the amount of operations and the demand of computer memory. The EGTG method optimizes the evaluation of dot products in solving least-squares problems. First, the dot product is identified as being a discrete integration of a known analytic integral. Then, the evaluation of the discrete integral is approximated by the evaluation of the analytic integral. This method should be applied when the evaluation of analytic integral needs less computational efforts than the discrete integration. To determine the equivalent data we developed two algorithms namely DOE and DOEg. The first one identifies the equivalent data of the whole linear systems while the second algorithm identifies the equivalent data in sub-systems of the entire linear systems. Each DOEg's iteration consists of one application of the DOE algorithm in a given subsystem. The algorithm DOE yields an interpolating surface that fits all data points allowing a global interpolation. On the other hand, the algorithm DOEg optimizes the local interpolation because it employs only the equivalent data while the other current algorithms for local interpolation employ all data. The interpolation methods using the equivalent layer technique was comparatively tested with the minimum curvature method by using synthetic data produced by prismatic source model. The interpolated values were compared with the true values evaluated from the source model. In all tests, the equivalent layer method had a better performance than the minimum curvature method. Particularly, in the case of bad sampled anomaly, the minimum curvature method does not recover the anomalies at the points where the anomaly presents high curvature. For data acquired at different levels, the minimum curvature method presented the worse performance while the equivalent layer produced very good results. By applying the DOE algorithm, it was possible to fit, using an equivalent layer model, 3137 gravity free-air data and 4941 total field anomaly data from the marine Equant-2 Project and the aeromagnetic Carauari-Norte Project, respectively. The DOEg algorithm was also applied in the same data sets optimizing the local interpolation. It is important to stress that none of these applications would have been possible without the concept of equivalent data. The ratio between CPU times (executing the programs with the same memory allocation) required by the minimum curvature method and the equivalent layer method in global interpolation was 1:31. This ratio was 1:1 in local interpolation.
Acesso aberto (Open Access)
Inversão da forma de onda orientada ao alvo
(Universidade Federal do Pará, 2016-09-16) COSTA, Carlos Alexandre Nascimento da; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
We propose a new target-oriented waveform inversion to estimate the physical parameters from a specific target in the subsurface from observed data from deviated-VSP acquisition or surface seismic data. Furthermore, we investigate a strategy to estimate the impulse responses from a local target in the subsurface from deviated-VSP acquisition or surface seismic data as an iterative sparse inversion approach, where the main feature of this strategy is that all multiple scattering in the data is used to enhance the illumination at target level. In these approaches we fit the upgoing wavefields observed at a specific level near the local target with the upgoing wavefields estimated at same depth level through convolution-type representation for the Green’s function. The main feature of the target-oriented waveform inversion is that we just need to know the up- and downgoing wavefields at the depth level above the target area to estimate the physical parameters for the area of interest. We show through numerical tests that the iterative sparse inversion approach does not require dense sources sampling to estimate the impulse responses from a target below a complex overburden, because of all the extra illumination via multiples. The physical parameters above the target area is not necessary to know if we use the data from deviated-VSP geometry of acquisition, but for surface seismic data we need to know a smooth physical parameter above the target area to estimate the up- and downgoing wavefields at depth level nearby the local target. For surface seismic data we used Joint Migration Inversion to estimate the up- and downgoing wavefields at depth level near the target area.
Acesso aberto (Open Access)
Migração 3-D Kirchhoff-Gaussian-Beam (KGB) pré-empilhamento no domínio da profundidade
(Universidade Federal do Pará, 2013-06-24) PEREIRA, Glauco Lira; CRUZ, João Carlos Ribeiro; http://lattes.cnpq.br/8498743497664023
The Gaussian Beam (GB) is an asymptotic solution of the elastodynamic equation in the paraxial vicinity of a central ray, which approaches better the wave field than the standard zero-order ray theory. The GB regularity in the description of the wave field, as well as its high accuracy in some singular regions of the propagation medium, provide a strong alternative to solve seismic modeling and imaging problems. In this thesis, i presenty a new procedure for pre-stack depth migration with true-amplitude, combining the flexibility and robustness of Kirchhoff migration type using superposition of Gaussian beams to represent the wave field. The proposed migration algorithm comprises in two stacking process: the first is the beam stack is applied to subsets of seismic data multiplied by a weight function defined such that stack operator has the same formulation of the integral of the Gaussian beams superposition; the second is a weighted diffraction stack by means of the Kirchhoff type integral having as input the GB stacked data. For these reasons it is called Kirchhoff-Gaussian-Beam (KGB) migration. The main characteristics that distinguish the KGB migration, during the first stage stacking, with other migration methods that also use the theory of Gaussian beams, is the use of the first Fresnel zone projected to limit the width of the subset of seismic traces (beam) using a second-order approximation of the reflection travel time. Examples are shown for applications on two-dimensional (2-D) and three-dimensional (3-D) synthetic seismic data, respectively, to the models Marmousi and SEG/EAGE salt dome data sets.
Acesso aberto (Open Access)
Modelagem 2,5D dos campos usados no Método Eletromagnético a Multi-Frequência - EMMF
(Universidade Federal do Pará, 2012-12-19) SILVA, Valdelírio da Silva e; DIAS, Carlos Alberto; http://lattes.cnpq.br/9204009150155131; RÉGIS, Cícero Roberto Teixeira; http://lattes.cnpq.br/7340569532034401
This thesis shows the 2.5D calculation of synthetic data from the Multifrequency Electromagnetic Method (EMMF, from the portuguese name). The work is presented in two main parts: the first one presents the details of the numerical methods used in the calculation of the fields generated by the Horizontal Loop of Current, placed on the surface of bi-dimensional earth models; the second part uses the achieved results to simulate EMMF method data, which are the real and imaginary parts of the radial magnetic field from the loop. In the second part, we observe the computed fields from several different models, including variations in their physical properties as well as in the geometry of their features. We aim at estimating the sensibility of these fields to the structures that can be found in a sedimentary basin. With this modeling, we are able to study the characteristics of the data and to find how the two components of the radial field, real and the imaginary, contribute with separate and complimentary information about the structures in the models. The results show that the radial magnetic field presents very good lateral resolution, for a measure taken from a fixed source. The observed capability of these data to resolve target structures will be extremely important for the future work with the inversion of such data, as well as for the building of apparent resistivity sections.
Acesso aberto (Open Access)
Modelagem e imageamento 2.5D no domínio do tempo através de diferenças finitas
(Universidade Federal do Pará, 2010) SILVA NETO, Francisco de Assis; SCHLEICHER, Maria Amélia Novais; http://lattes.cnpq.br/4767998352165705; COSTA, Jessé Carvalho; http://lattes.cnpq.br/7294174204296739
This thesis discuss modeling and imaging of seismic wavefields in 2.5D using finite-differences to solve numerically the wave equation. Modeling in 2.5D is extended to anisotropic elastic media with an arbitrary class of symmetry. The sources of the wavefield are generalized to simulate of explosive, dipole and double-couple distributions. The acquisition geometry is not required to coincide with a symmetry plane. Reverse time migration in 2.5D is implemented in conjunction with a new imaging condition based on the asymptotic analysis of the classical correlation imaging condition. The new imaging condition is designed to improve the amplitudes in reverse time migration (RTM) images, and to reduce back-scattering artifacts. Numerical experiments indicate that 2.5D RTM improves the resolution of the migrated images when compared to its 2D counterpart, and that the proposed imaging condition was effective improving the amplitudes and reducing back-scattering artifacts.
Acesso aberto (Open Access)
Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais
(Universidade Federal do Pará, 1993-02-05) GONZÁLEZ CARRASQUILLA, Antonio Abel; RIJO, Luiz; http://lattes.cnpq.br/3148365912720676
South America presents several unique geomagnetic features, one of which is the Equatorial Electrojet, a current system which extends east - west in Northern Brazil, for almost 3500 km. Considering the fact that the influence of the Equatorial Electrojet can be detected at great distances from its central axis, it is important to understand its effects on magnetotelluric exploration in Brazil. In magnetotelluric prospecting, the influence of the equatorial electrojet has been modelled for both one and two - dimensional geological media, employing both analytical and numerical techniques solutions such as finite element and finite difference. Three-dimensional geological media have been modeled by "thin – layers "using a "thin sheet" algorithm. Lines of current, gaussian electrojets and undulatory electrojets have been used as induction sources to simulate the equatorial electrojet in these algorithms. In this thesis the principal objective is to model the effects of the equatorial electrojet on three - dimensional structures commonly found in geophysical prospecting. To accomplish this, we have computed numerical solutions of the integral equation for three - dimensional media using the inductive sources mentioned before. As previous works have indicated, our results also show that the influence of the equatorial electrojet is prominant only for frequencies lower than 10-1 Hz. This influence decreases with distance but is detectable up to as far as 3000 km from the center of the electrojet. Thus, the presence of peaks in the apparent resistivity profiles over a homogeneous half - space indicates that the influence of the electrojet is more noticeable for this kind of medium. These peaks display different amplitudes for each type of simulated electrojet, and the peak locations also change from one electrojet to the other. However, when we use more geologically realistic one - dimensional models such as a stratified media, the effect of the electrojet source diminishes considerably and the results do not vary greatly for the different kinds of electrojet employed in the model. This effect is caused by the electromagnetic energy dissipation due to the presence of stratified conductive layers within the media. Within the 3000 km region, the three - dimensional electromagnetic response can be larger or smaller than the plane wave response, depending on location body, frequency, kind of the electrojet and geology. When the apparent resistivity is larger than the plane wave response, there is a spread between the one and the three-dimensional sounding curves caused by the electrojet, as well as a widening of the profile anomaly caused by the three - dimensional inhomogeneities relative to the one due to a plane wave. When the apparent resistivity is less the two kinds of sources yield anomalies approximately equal as well as a shortening of the profile anomaly. On the other hand, the phase usually shows an inverted way related to the apparent resistivity. This means that when one phase goes up the apparent resistivity goes down, and vice - versa. Similarily at high frequencies, the one and three - dimensional phases are offset, while at low frequencies they are the same, except for the undulating electrojet with undulation factor α = -2 and -3. Our results also show that the geometric characteristics of three-dimensional structures, such as their orientation relative to the direction of the electrojet and the dimensions of their principal direction, change the response due to the electrojet source as compared to a plane wave source. For example, when the three - dimensional structure is rotated 90 degrees, relative to the direction of the electrojet, there is a change in polarization (of the electric and magnetic fields) but there are no changes in the values of apparent resistivity at the center of symmetry of the structure. When the measurements are taken away from the center of symmetry changes in the apparent resistivity are shown when compared with the unrotated structure. This is due to the persistente of the galvanic effects at the center of the structure and the presence of inductive effects at the periphery of the three - dimensional body. When we elongate the principal direction of the three - dimensional structure, the magnetotelluric soundings start to approach those of two-dimensional structures. This is more noticiable in the XY polarization. Even so, the responses of the models tested are still quite different from the responses of the quasi two - dimensional structures. Nevertheless, the effects of the electrojet in structures with elongated principal direction are very similar to those present in smaller structures, considering the differences between the soundings of both kinds of structures. On the other hand, the apparent resistivities of this kind of elongated structures show a great increase at the edges of these structures, both for the plane wave or electrojet sources. This effect is caused by the chanelling of the current along the main direction of the structure. The modelling of gelogical structures in the Marajó Basin, confirms that the effects of the electrojet can be detected even in small horst and graben structures located at great distances from electrojet center. Electrojet effects are noticeable in both one and three - dimensional media for two frequency bands, one near 10-1 Hz and another band, between 10-3 Hz to DC. This possibly reflects the influence of a resistive geological basement and a resistive lower crust, respectively. The results of the analysis, using the electrojet as an inductive source show that at low frequencies the computed responses from the magnetotelluric soundings can be strongly distorted, both by galvanic effects in the three - dimensional structure as well as the presence of the electrojet. Therefore, if the equatorial electrojet effects are not accounted for, a misleading interpretation magnetotelluric data will result. Similarily, a three - dimensional model should be used to interpret the data, instead of one-dimensional Tikhonov-Cagniard model.
Acesso aberto (Open Access)
Modelagem sísmica e inversão na presença de anisotropia
(Universidade Federal do Pará, 1993-11-04) COSTA, Jessé Carvalho; SCHOENBERG, Michael
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.
Acesso aberto (Open Access)
Modelo fractal para resistividade complexa de rochas: interpretação petrofísica e aplicação à exploração geoelétrica
(Universidade Federal do Pará, 1995-12-21) ROCHA, Brígida Ramati Pereira da; HABASHY, Tarek Mohamed
Rocks containing disseminated metallics or clay particles in natural environment where electrolytic solutions fill the pore spaces, show a certain type of polarization at low frequencies known as induced electrical polarization. In this thesis, a new model to describe the electrical polarization on rocks was developed, not only for low frequencies, but spanning the entire electromagnetic spectrum used in geolectric prospection. This new model encompasses most of the other commonly used models as special cases, and overcomes some of the known limitations. The proposed circuit analog includes a non-linear impedance r(iwtf)-1 which simulates the effects of the rough surface of the interfaces between the blocking grains (metallic or clay particles) and the electrolyte. This generalized Warburg impedance is in series with the resistance of the blocking grains and both are shunted by the double layer capacitance. This combination is in series with the resistance of the electrolyte in the blocked pore passages. The unblocked pore paths are represented by a. resistance which corresponds to the normal DC resistivity of the rock. The parallel combination of this resistance with the "bulk" sample capacitance is finally connected in parallel to the rest of the above-mentioned circuit. The parameters of this model include the DC resistivity (p0), the chargeability (m), three relaxation times (T, Tf and T0), a grain resistivity factor (δr) and the frequency exponent (η). The fractal relaxation time (Tf) and the frequency exponent (η) are related to the fractal geometry of the rough pore interfaces between the conductive grains (metallic or clay minerals which are blocking the pore paths) and the electrolyte. The relaxation time T is a result of the low-frequency relaxation of the electrical double layers formed between the electrolyte and the crystals, whereas T0 is a macroscopic relaxation time of the "bulk" sample. The grain resistivity factor (δr) relates the resistivity of the conductive grains with the DC resistivity value of the rock. The DC resistivity of the rock and δr are related to the porosity, the electrolyte conductivity and the volumetric ratios between the matrix and the conductive grains. The model was tested over a wide range of frequencies against experimental data obtained for amplitude and phase of resistivity or conductivity as well as for the complex dielectric constant. The data used in this work were obtained from digitizing published experimental data, obtained by several authors from sedimentary, metamorphic and igneous rocks. The results show that the parameters of this model are related to textural and mineralogical aspects of the rocks. This model was introduced firstly as the intrinsic electric property of a homogeneous and polarizable half-space, and it was demonstrated, in this thesis, that the response observed at the surface is equivalent to the intrinsic property of the polarizable medium, been the electromagnetic coupling irrelevant to frequencies lower than 104 Hz. Next, the polarizable medium was embedded as an intermediate layer between two non-polarizable layers with the same De resistivity. The response obtained shows that the frequency exponent of the fractal medium could be determined even when the polarizable medium is at a considerable depth in relation to the dipole-dipole length. This justifies the use of simple models developed to explain the response of laboratory samples to fit field data, and that is being used without a right justification. These results shows the importance of the proposed model to the geoelectric prospection.
Acesso aberto (Open Access)
Região do espaço que mais influencia em medidas eletromagnéticas no domínio da frequência: caso de uma linha de corrente sobre um semi-espaço condutor
(Universidade Federal do Pará, 1994-07-28) BRITO, Licurgo Peixoto de; DIAS, Carlos Alberto; http://lattes.cnpq.br/9204009150155131
One of the major interpretation problems in geophysics is to determine the region in the subsurface which generates the main part of the signal. In this thesis, the position and size of this region, hereinafter called the main zone, have been found by modelling an electromagnetic system in which the source is an infinite line of electric current, extended over a conductive half-space. The earth has been modelled as a conductive half-space with an inhomogeneity in it as being an infinite layer or a prism of infinite length in the direction of the source line. The signal in the receiver of an electromagnetic system over a conductive homogeneous half-space is different from the one taken over the half space including an inhomogeneity. This difference is a function of the position of the inhomogeneity in relation to the transmiter-receiver system, besides other parameters. Therefore, with the other parameters fixed, there will be a specific position where this difference will maximize. Since this position depends on conductivity contrast, inhomogeneity dimensions and on source frequency, instead of a single position one will have a region where the inhomogeneity will give the maximum contribution to the measured signal. This region is called the main zone. Once the main zone is identified, the targets in the subsurface can be more precisely located. Usually they are conductive parts of the earth with some specific interest. One can facilitate the exploration and reduce production costs if these conductors are well identified during prospecting. A detectability function (∆) has been defined to measure the contribution to the signal due to the inhomogeneity. The ∆ function has been computed using amplitude and phase of the magnetic field components: Hx and Hz which are, respectively, the tangential and the normal to the earth's surface. The size and position of the main zone has been identified using the extremais of the ∆ function, which change with conductivity contrast, and the inhomogeneities' size and depth. Electromagnetic fields for one-dimensional models were calculated using a hybrid form, numerically solving the integrals that were obtained analytically. Two-dimensional models were computed numerically, by the finite elements technique. The maximum values of ∆ function, computed with amplitude of Hx, have been chosen to locate the main zone. This shows more stable results than other amplitude and phase components, both for one and two-dimensional models, when physical properties and geometric dimensions are changed. For the one-dimensional model, where the inhomogeneity is an infinitely extended horizontal layer, the depth of its central plane was found to be po = 0.17 δo, where po is the depth of this central plane and δo is the skin depth for the plane wave (in an homogeneous half-space having a conductivity σ1 equal to that of the backgound, and the frequency w corresponding to the maximum value of ∆ calculatede for the amplitude of Hx). For two-dimensional inhomogeneities, the co-ordinates of the main zone central axis was found to be do = 0,77 r0 (where do is the horizontal distance from this axis to the source) and po = 0,36 δo (where po is the depth of this central axis), with r0 being the source-receiver separation and δo the skin depth in the same conditions as in the one-dimensional case. If the values of r0 and δo are known, it is possible to determine (do, po). Associating each value of ∆ function (calculated using the amplitude of Hx) with the values of d = 0,77 r and p = 0,36 δ for each r and w used to generate ∆, a method to locate the main zone is sugested. The isovalue curves of ∆ are plotted to construct sections of ∆. These sections indicate the conductors position and provide some helpful insight into their geometric forms when the values of ∆ get dose to the maximum.
Acesso aberto (Open Access)
Simulação de perfis nucleares de poço em formações complexas
(Universidade Federal do Pará, 1993-04-26) SILVA, Jadir da Conceição da; EVANS, Hilton Bernard
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.

Navegar

Navegando Teses em Geofísica (Doutorado) - CPGF/IG por Agência de fomento "PETROBRAS - Petróleo Brasileiro S.A."