Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-20T14:53:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Ausência de palavras-chave on 2014-08-07T12:41:40Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-13T14:27:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T14:39:54Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Made available in DSpace on 2014-08-20T14:39:54Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) Previous issue date: 2003 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CTPETRO - Plano Nacional de Ciência e Tecnologia do Setor Petróleo e Gás Natural / FINEP - Financiadora de Estudos e Projetos / O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais. / The CRS stacking method simulates ZO seismic sections from multi-coverage data and does not dependente on a macro-velocity model. For 2-D medium the stacking traveltime depends on three parameters: the emergence angle of the normal ray (with respect to the measurement surface normal) and the wavefront curvatures of two hypothetical waves, called Normal-Incidence-Point (NIP) wave and Normal (N) wave. The CRS method consists of summing the amplitudes of the seismic traces in the multicoverage data along the surface defined by CRS stacking traveltime which that fits best the data set. The result of the CRS stack is assigned to points of a grid pre-defined in the ZO section. As the result obtain a simulated ZO section. This means that for each point of the ZO section must be estimated the three optimal parameters that yield the maximum coherence between the events of seismic reflection. In this Thesis I present formulae for the 2-D CRS method and for the NMO velocity that consider the topography of the measurement surface. The algorithm is based on the optimization strategy divided into three steps: 1) To search for the emergence angle and the curvature of the NIP wave, by applying a global optimization, 2) to search for the curvature of the N wave, by applying global optimization, and 3) to refine the initial parameters estimated in first two steps by applying local optimization. In the first two steps is used the Simulated Annealing (SA) algorithm and in the third step the Variable Metric (VM) algorithm. For the case of a measurement surface with smooth topography the curvature of this surface is included in the 2-D CRS stack formalism. This CRS algorithm implemented was applied to synthetic data set. The result is a simulated ZO section of high quality, with a high signal-to-noise ratio, and the estimative of the parameter triplet. It is performed a sensibility analysis for the new CRS stacking traveltime with respect to the curvature in several points of the curved measurement surface. This study showed that the CRS traveltime is more sensitive for fast midpoints of the central points and larger offsets. The expressions for the NMO velocities presented here is applied to estimate the interval velocities and the depth of the reflectors for 2-D model with a smooth topography. For the inversion of the velocities and the depth of the reflectors is considered the Dix-type inversion algorithm. The NMO velocity for a curved measurement surface deserves to best estimate the velocities and the depths of the reflectors than NMO velocities referred to planar surfaces. Also, I present an introduction to 3-D stack. In this case, the stacking traveltime depends on eight parameters. These parameters can be obtained by using some parameter-search strategies that I have showed in this Thesis. The combination of the strategy of the Traveltime Approximations and the strategy of Arbitrary Curvatures is used to apply 3-D CRS stack successful in synthetic and real data sets, respectively.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-09T12:43:26Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodosEletrico.pdf: 10195935 bytes, checksum: a6776849fa107eb79f60ab784dd5f935 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:13:45Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T16:47:31Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodosEletrico.pdf: 10195935 bytes, checksum: a6776849fa107eb79f60ab784dd5f935 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T17:00:47Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodosEletrico.pdf: 10195935 bytes, checksum: a6776849fa107eb79f60ab784dd5f935 (MD5) / Made available in DSpace on 2014-08-20T17:00:47Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodosEletrico.pdf: 10195935 bytes, checksum: a6776849fa107eb79f60ab784dd5f935 (MD5) Previous issue date: 1998 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / DCT - Departamento de Cooperação Científica, Técnica e Tecnológica / Este trabalho trata da aplicação dos métodos geofísicos eletroresistivo e eletromagnético no município de Altamira, região central do Estado do Pará, com a finalidade de fazer o mapeamento da espessura das camadas aqüíferas, bem como a ocorrência de intrusões de diabásio na área do município, permitindo uma melhor locação de pontos para perfuração com fins de captação à água subterrânea. Na área estudada, verificam-se afloramentos de rochas cristalinas do embasamento coberto por sedimentos com idade variando desde o Paleozóico ao Quaternário, com esperança de se obter água nos sedimentos do Devoniano Inferior, mais precisamente na formação Maecuru. Foram executadas 34 sondagens elétricas verticais com arranjo Schlumberger com AB/2 variando de 100 a 500 metros. 108 sondagens eletromagnéticas transientes também foram executadas com arranjo "loop coincidente", com 50 metros de lado. Por problemas eletrônicos adicionados aos ruídos de fenômenos supraparamagnéticos, os resultados das sondagens eletromagnéticas apresentaram-se muito distorcidos, não tendo sido possível extrair deles informação geológica de relevância. Por isso, foram criados modelos eletromagnéticos sintéticos representando a situação geológica da área. Os dados obtidos no campo em forma de curvas de resistividade aparente, foram plotados em papel bilogarítmo e interpretados com a ajuda de algoritmos direto e inverso os quais considera modelos de camadas horizontalmente estratificadas, homogêneas e isotrópicas. Os resultados da interpretação permitiram estabelecer modelos geoelétricos aproximados refletindo a situação geológica na subsuperfície. Os resultados da interpretação das sondagens elétricas, foram confrontados com os resultados da interpretação qualitativa do mapa de isogamas do campo magnético total, permitindo delimitar com maior precisão a zona de incidência de soleiras de diabásio, situada a norte da área. Dos quatro horizontes geoelétricos distinguidos, o terceiro é essencialmente constituído por arenitos, apresentando, por isso, melhores condições para armazenamento de água subterrânea. / This work due application of the electroresistivity and electromagnetic geophysical methods in the municipal district of Altamira, in the central area of the State of Pará, with the purpose of mapping the thickness of the aquifer layers and occurrence of diabase intrusions for finding a groundwater wells. In this area are many outcrops of basement crystalline rocks covered with sediments of ages from the Paleozioc to Quaternary. The objective is to obtain water from the sediments of Inferior Devonian, more exactly in the Maecuru Formation. There were made thirty four vertical electric soundings using the Shlumberger arrangement, the AB/2 distance varying from 100 to 500 meters. Also there were made 108 transient electromagnetic soundings with "coincident loop" arrangement, using an square loop of 50 meters. Due to electronic problems and supraparamagnetic noises, the results of this electromagnetic soundings became very distorted, impossible to get from them important geologic information. Then, electromagnetic synthetic models were created to represent the geology of the area. The data were obtained from the field were presented as apparent resistivity curves. They were plotted using a log-log scale and interpreted using the EGSLIB/SEV1D first and then with the EGSLIB/SEV1DIN program. The models used were horizontal, isotropic and homogeneous layers. The interpretation results allowed to determine approximate geoelectrical models which reflects the underground geological situation. The interpretation results from the electric surveys were confronted with the results of the qualitative interpretation from the isogams map of the total magnetic field, allowing to delimit with more precision the zone of incidence of diabase sills, located at the north of the area. From the four distinguished geoeletric horizons, the third is essentially constituted by sands, therefore it is best for groundwater storage.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T12:28:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:13:21Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-14T14:23:37Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T17:04:51Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Made available in DSpace on 2014-08-22T17:04:51Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) Previous issue date: 1991 / Utilizando-se dados do campo de Camorim (Bacia de Sergipe-Alagoas), foi testado e aplicado um conjunto de técnicas estatísticas multivariantes (análises de agrupamentos, de componentes principais e discriminante) no intuito de identificar as fácies, previamente definidas em poços testemunhados, por meio dos perfis, viabilizando o reconhecimento das mesmas nos demais poços não testemunhados da área. A segunda etapa do processo de determinação das fácies consistiu no emprego de métodos auxiliares (análises composicional e de seqüência de fácies), que, combinados com as técnicas multivariantes, propiciaram melhores resultados na calibração rocha-perfil. A determinação das fácies, uma vez estabelecida, possibilitou o refinamento do processo de avaliação de formações ao viabilizar o exame de cada fácies-reservatório, isoladamente. Assim, esse procedimento tornou possível a escolha, para cada litologia, dos parâmetros utilizados na interpretação dos perfis ao mesmo tempo em que permitiu a totalização em separado dos valores de espessura, porosidade e saturação dos fluidos, bem como a adoção de diferentes valores de corte (cut-offs) para cada grupo considerado. Outras aplicações incluíram a melhoria na estimativa da porosidade e da permeabilidade, a adaptação de algoritmos para o cálculo preliminar de porosidade, a confecção de mapas de fácies e a geração automática de seções estratigráficas. Finalmente, foram destacadas a perspectiva de integração desse estudo com sistemas estatísticos de descrição de reservatórios, outras técnicas de determinação de fácies em desenvolvimento e a retomada da utilização de métodos estatísticos multivariantes em dados de perfis, como ferramenta de exploração. / A series of multivariate statistical techniques (cluster, principal component and discriminant analysis) was tested and applied to well log data from the Camorim field (offshore Sergipe State, Brazil) in order to identify facies previously defined through core description. The second step in the process of facies determination was supported by auxiliary methods (compositional and sequence facies analysis), which produced better results in the calibration of rock-log data, when combined with the multivariate techniques. The facies determination, once established, permits the refinement of the process of formation evaluation, enabling the examination of each facies separately. This procedure made it possible to choose, for each lithology, the parameters used in log interpretation. In parallel, this process allowed the summation of thickness, porosity, fluid saturation and the adoption of different cut-offs for each group, separately. Other applications included: enhancement in the estimation of porosity and permeability, the adaptation of algorithms designed for fast porosity estimation, the mapping of variables useful in the characterization of the vertical variability of the reservoir rocks and the automatic generation of stratigraphic sections. Finally, the possibility of integration of the work from this study with statistical systems of reservoir description, other facies determination techniques currently being developed and the utilization of multivariate statistical methods in well log data, as an auxiliary exploratory tool, were illustrated.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T16:47:14Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoGeoeletricaArea.pdf: 3315484 bytes, checksum: 279133bbca3dc766b7e2c2c5bc9bb8b2 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:30:33Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T13:24:15Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoGeoeletricaArea.pdf: 3315484 bytes, checksum: 279133bbca3dc766b7e2c2c5bc9bb8b2 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T17:09:08Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoGeoeletricaArea.pdf: 3315484 bytes, checksum: 279133bbca3dc766b7e2c2c5bc9bb8b2 (MD5) / Made available in DSpace on 2014-08-22T17:09:08Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoGeoeletricaArea.pdf: 3315484 bytes, checksum: 279133bbca3dc766b7e2c2c5bc9bb8b2 (MD5) Previous issue date: 2005 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Foi desenvolvido um estudo da subsuperfície rasa para investigar a possibilidade de contaminação do subsolo por líquidos usados no tratamento de couros em área industrial. O estudo foi realizado na área da indústria Couro do Norte, Bairro de Maracacuera, no Distrito Industrial de Icoaraci (município de Belém-Pará). Durante o estudo foram feitas medidas para a determinação do fluxo de água do primeiro aqüífero e aplicada metodologia geofísica com a finalidade de detectar a possível contaminação do subsolo por poluentes liberados pelo curtume. As metodologias geofísicas utilizadas foram: eletrorresistividade (na forma de imageamento e través de sondagens verticais), potencial espontâneo, eletromagnético (Slingram) e radar de penetração no solo (GPR). A interpretação dos dados eletromagnéticos permitiu identificar a presença de zonas mais condutivas (provável pluma de contaminação) e zonas mais resistivas (blocos lateríticos). A distribuição dos potenciais medidos através do Método do Potencial Espontâneo indicou o sentido do fluxo subterrâneo local. Com as Sondagens Elétricas Verticais determinaram-se os valores de resistividade em subsuperfície para os diferentes horizontes correlacionando-os sempre com as informações obtidas no Caminhamento Elétrico. Através dos dados de GPR identificou-se a presença de blocos lateríticos, tubulações de metal e zonas com elevada atenuação (zona de sombra) que indicam a presença de contaminação. / A study covering the shallow subsurface was developed in order to investigate the possibility of underground contamination by liquids used in the treatment of animal skin. The study was carried out inside the area of the Couro do Norte industry, located in the Icoaraci Industrial District (Belém – Pará). During the study they were made measures for the determination of the flow of water of the first aquifer and applied geophysical methodology with the purpose of detecting the possible contamination of the underground for pollutant liberated by the tanning. The geophysical methodology applied was: resistivity (imaging and vertical electric soundings), spontaneous potential, electromagnetic (Slingram) and ground penetrating radar (GPR). Analyses of electromagnetic data lead to the identification of conductive zones (Contamination) and resistive zones (sandy material). Spontaneous potential gave information on the direction of local underground flow of water. Vertical electric soundings allowed to know the resistivity distribution in subsurface. The resistivity sounding models could be correlated to the electrical images obtained in electrical profiling. Through the GPR data was possible to identify the presence of lateritic blocks, metal piping and zones with high attenuation (shadow zones) that indicate the presence of contamination.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-26T13:04:20Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Migracao3DKirchhoff.pdf: 13210015 bytes, checksum: 55e96c2ad9e534d8dce3117b9d9f3623 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T14:52:53Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-08T12:41:18Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Migracao3DKirchhoff.pdf: 13210015 bytes, checksum: 55e96c2ad9e534d8dce3117b9d9f3623 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-08T14:44:47Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Migracao3DKirchhoff.pdf: 13210015 bytes, checksum: 55e96c2ad9e534d8dce3117b9d9f3623 (MD5) / Made available in DSpace on 2014-09-08T14:44:47Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Migracao3DKirchhoff.pdf: 13210015 bytes, checksum: 55e96c2ad9e534d8dce3117b9d9f3623 (MD5) Previous issue date: 2013 / CENPES - Centro de Pesquisas Leopoldo Américo Miguez de Mello / PETROBRAS - Petróleo Brasileiro S.A. / O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa na solução de problemas de modelagem e imageamento sísmicos. Nesta Tese, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB). As principais características que diferenciam a migração KGB, durante a realização do primeiro empilhamento, de outros métodos de migração que também utilizam a teoria dos Feixes Gaussianos, são o uso da primeira zona de Fresnel projetada para limitar a largura do feixe e a utilização, no empilhamento do feixe, de uma aproximação de segunda ordem do tempo de trânsito de reflexão. Como exemplos são apresentadas aplicações a dados sintéticos para modelos bidimensionais (2-D) e tridimensionais (3-D), correspondentes aos modelos Marmousi e domo de sal da SEG/EAGE, respectivamente. / 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.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-28T13:29:48Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoKirchhoffPreEmpilhamento.pdf: 7520606 bytes, checksum: 04380685465a2a6e5ef4c422ca54216a (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir assuntos on 2014-08-06T16:21:01Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-08T14:11:40Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoKirchhoffPreEmpilhamento.pdf: 7520606 bytes, checksum: 04380685465a2a6e5ef4c422ca54216a (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-09T12:08:59Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoKirchhoffPreEmpilhamento.pdf: 7520606 bytes, checksum: 04380685465a2a6e5ef4c422ca54216a (MD5) / Made available in DSpace on 2014-09-09T12:08:59Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_MigracaoKirchhoffPreEmpilhamento.pdf: 7520606 bytes, checksum: 04380685465a2a6e5ef4c422ca54216a (MD5) Previous issue date: 2007 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CTPETRO - Plano Nacional de Ciência e Tecnologia do Setor Petróleo e Gás Natural / FINEP - Financiadora de Estudos e Projetos / A teoria dos feixes gaussianos foi introduzida na literatura sísmica no início dos anos 80 por pesquisadores russos e tchecos, e foi originalmente utilizada no cálculo do campo de ondas eletromagnéticas, baseado na teoria escalar da difração. Na teoria dos feixes gaussianos, o campo de ondas sísmicas é obtido por uma integral, cujo o integrando é constituído de duas partes, a saber: (1) as amplitudes dos campos das ondas na vizinhança do ponto de observação e (2) a função fase de cada um desses campos de ondas, que neste caso é representada por um tempo de trânsito paraxial complexo. Como ferramenta de imageamento, mais precisamente como operador de migração, os primeiros trabalhos usando feixes gaussianos datam do final da década de 80 e início dos anos 90. A regularidade dos campos de ondas descritos pelos feixes gaussianos, além de sua alta precisão em regiões singulares do modelo de velocidades, tornaram o uso de feixes gaussianos como uma alternativa híbrida viável para a migração. Nesse trabalho, unimos a flexibilidade da migração tipo Kirchhoff em profundidade em verdadeira amplitude com a regularidade da descrição do campo de ondas, representado pela sobreposição de feixes gaussianos. Como forma de controlar de forma estável quantidades usadas na construção de feixes gaussianos, utilizamos informações advindas do volume de Fresnel, mais precisamente a zona de Fresnel ao redor do ponto de reflexão e a zona de Fresnel projetada, localizada ao redor do ponto de registro do sismograma e cuja a informação se encontra nas curvas de reflexão de dados sísmico. Nosso processo de migração pode ser chamado como uma migração Kirchhoff em verdadeira amplitude usando um operador de feixes gaussianos. / The Gaussian Beam (GB) concept was introduced in the seismic literature by Russian and Czech researchers in the begining of the 80’s. This theory, which by its turn was based on the scalar electromagnetic diffraction theory, is in fact a (zero order) complex paraxial ray theory, designed to satisfactorilly describe the seismic wavefield propagation beyond the standard zero order ray theory, up to then the only theory used to describe the high frequency seismic wavefield propagation in smoothed velocity models. As an imaging tool, the first works to deal with GB’s were published in the end of the 80’s and in the begining of the 90’s. The regularity in the description of the wavefield by GB’ s, as well as its high accuracy in some singular regions of the velocity model, transformed the use of GB’s into a viable hybrid alternative in the migration theory. In this work, we unite the flexibility in imaging of the true amplitude prestack Kirchhoff depth migration with the regularity in the description of the wavefield by a superposition of GB’s. As a way of controlling in a very stable way some quantities used in the construction of the beams, we have made use of some informations based on the Fresnel volume elements, more especifically speaking the Fresnel zone radius around the reflection point in depth and its counterpart, projected towards the acquisition surface. This information is centred around the recording point of the seismogram and is also present in the seismic data reflection traveltime curves. Our migration process can be named a true amplitude prestack Kirchhoff depth migration using GB’s as Green function, namely KGB-PSDM.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-05T13:45:51Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseComparativaTres.pdf: 6132554 bytes, checksum: 3f873423cc9d66ad912a70a0cf13c3b9 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:17:46Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-29T13:39:59Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseComparativaTres.pdf: 6132554 bytes, checksum: 3f873423cc9d66ad912a70a0cf13c3b9 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-09T14:38:20Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseComparativaTres.pdf: 6132554 bytes, checksum: 3f873423cc9d66ad912a70a0cf13c3b9 (MD5) / Made available in DSpace on 2014-09-09T14:38:20Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AnaliseComparativaTres.pdf: 6132554 bytes, checksum: 3f873423cc9d66ad912a70a0cf13c3b9 (MD5) Previous issue date: 1985 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / FINEP - Financiadora de Estudos e Projetos / Tradicionalmente, o método dos mínimos quadrados tem sido empregado na inversão não linear de dados de campo potencial. No caso em que as observações dos campos gravimétrico ou magnético contém apenas ruído Gaussiano. O método dos mínimos quadrados não apresenta problemas. Entretanto, quando as observações são perturbadas por ruído não Gaussiano, ou mesmo por ruído não aleatório, como é o caso de muitos ruídos geológicos, o método dos mínimos quadrados torna-se bastante ineficiente, e métodos alternativos devem ser empregados a fim de produzir interpretações realísticas. Neste trabalho, uma comparação é feita entre os métodos dos mínimos quadrados, dos mínimos absolutos e do ajuste-M, aplicados à inversão não linear de dados de campo potencial. A comparação é efetuada usando-se dados teóricos, onde diversas situações geológicas são simuladas. Os resultados mostram que na presença de ruído geológico, caracterizado por pequeno corpo raso acima do corpo principal, ou por corpo grande, adjacente ao corpo principal, o ajuste-M apresenta desempenho muito superior ao dos mínimos quadrados e dos mínimos absolutos. Na presença de ruído Gaussiano, entretanto, o ajuste-M tem um desempenho inferior aos outros dois métodos. Como o ruído Gaussiano é um ruído branco, parte dele pode ser removido por um filtro passa baixa adequado, sem muita perda do sinal, o que não ocorre com o ruído geológico que contém componentes importantes de baixo número de onda. Desse modo o ajuste-M se torna uma ferramenta importante na interpretação de áreas geologicamente complexas, onde é comum a contaminação das anomalias por ruído geológico. Os três métodos em estudo são aplicados a uma anomalia magnética real causada por uma intrusão de diabásio em forma de dique, em sedimentos arenosos da formação Piauí na Bacia do Parnaíba. Os três métodos apresentaram resultados semelhantes indicando que tanto o nível de ruído Gaussiano como geológico são baixos nesta anomalia. / Nonlinear inversion of potential field data has been traditionally accomplished by the least squares method. As far as the anomalous field is corrupted by Gaussian random noise. Least squares inversion has a good performance. However, when the data are contaminated by non Gaussian noise, which is the case of most geological noise, the least squares method presents an extremely poor performance. As a result, alternative methods must be employed in this case in order to produce realistic and meaningful interpretations. This paper presents a comparison among the least squares method, the minimum absolute error and M-fitting applied to non linear inversion of potential field data. The analysis is performed using theoretical data generated by synthetic models simulating several geological settings. The results show that in the presence of geological noise represented either by small shallow bodies above the main body, or by large interfering bodies adjacent to the main body, M-fitting presents a much batter performance as compared with the least squares or the minimun absolute error methods. In the presence of Gaussian random noise, however M-fitting has a poor performance. Since a Gaussian noise is a white noise, a law pass filter applied to the observed data would remove part of the Gaussian noise with a minimum loss of the low wavenumber signal. On the other hand, most geological noise have important low wavenumber spectral components so that this noise cannot be eliminated without a significant loss of signal. Therefore, the M-fitting method may become an important interpretation tool when applied to complex areas (where anomalies are usually contaminated by geological noise) provided the data have been previously filtered by a suitable law pass filter. All three methods analysed in the paper are applied to a real magnetic anomaly due to a dike of diabasic rock intruded in sandstones and shales from the Piauí Formation, in Parnaíba Basin, Brazil. All three methods yielded similar interpretations which are consistent with the available a priori geological information. The fact that all methods produced similar results indicates a low level of geological and Gaussian random noise in the data.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T15:42:02Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoAlgoritmosConversao.pdf: 3250632 bytes, checksum: a0e70466bb0b3ec4d4b51f6d03df6a9b (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Incluir palavras-chave on 2014-08-07T13:31:22Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T12:58:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoAlgoritmosConversao.pdf: 3250632 bytes, checksum: a0e70466bb0b3ec4d4b51f6d03df6a9b (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T14:49:27Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoAlgoritmosConversao.pdf: 3250632 bytes, checksum: a0e70466bb0b3ec4d4b51f6d03df6a9b (MD5) / Made available in DSpace on 2014-09-10T14:49:27Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoAlgoritmosConversao.pdf: 3250632 bytes, checksum: a0e70466bb0b3ec4d4b51f6d03df6a9b (MD5) Previous issue date: 2013 / Ainda hoje, a migração em tempo é o processo de imageamento substancialmente empregado na indústria do petróleo. Tal popularidade é devida ao seu alto grau de eficiência e robustez, além de sua habilidade em focalizar refletores nos mais variados ambientes geológicos. Entretanto, em áreas de alta complexidade geológica a migração em tempo falha de tal forma que a migração em profundidade e um campo de velocidade em profundidade são indispensáveis. Esse campo é geralmente obtido através de processos tomográficos partindo de um campo de velocidade inicial. A conversão de campos de velocidade de tempo para profundidade é uma forma rápida de se obter um campo inicial mais consistente geologicamente para tais processos. Alguns algoritmos de conversão tempo-profundidade recentemente desenvolvidos baseados no traçamento de raios-imagem são revistos e um algoritmo alternativo baseado na propagação da frente de onda-imagem é proposto. Os algoritmos são aplicados a dados sintéticos bidimensionais e avaliados de acordo com suas eficiência e acurácia, destacando suas vantagens, desvantagens e limitações na obtenção de campos de velocidade em profundidade. / Even today, time migration is the imaging process substantially employed in the oil industry. Such popularity is due to its high degree of efficiency and robustness, and its ability to focus reflectors in various geological environments. However, in areas of high geological complexity, time migration fails so that depth migration and a depth velocity field are indispensable. This field is usually obtained through tomographic iterative processes starting from an initial velocity field. The conversion of velocity fields from time to depth is a quick way of obtaining a geologically more consistent initial velocity field for such processes. Some algorithms for time-to-depth conversion based on newly developed image-ray tracing are reviewed and an alternative algorithm based on propagation of the image-wavefront is proposed. The algorithms are applied to two-dimensional synthetic data and evaluated according to their efficiency and accuracy, highlighting their advantages, disadvantages and limitations in obtaining depth velocity fields.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-12T13:56:10Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar as palavras-chave on 2014-08-07T13:38:08Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T13:51:07Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T14:49:50Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) / Made available in DSpace on 2014-09-10T14:49:50Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ComparacaoFiltrosVelocidade.pdf: 4324452 bytes, checksum: 90dbfc8b340951c647cbf95c43f4bab3 (MD5) Previous issue date: 2004 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t. / The geological motivation of this work is the imaging of sedimentary basin structures of the Amazon region, where the generation and accumulation of hydrocarbons is related to the existence of diabase sills. The seismic motivation is the fact that these intrusive rocks present a great impedance contrast with respect to the host rock, what gives rise to external and internal multiples, with primary-like amplitudes. The seismic signal of the multiples can predominate over the primary reflection signals from deeper interfaces, making difficult the processing, interpretation and imaging of seismic sections. In this work we study the attenuation of multiples in common shot (CS) sections by the comparison of two methods. The first one is the combination of the Wiener-Hopf-Levinson (WHLP) and the common reflection surface (CRS) stacking techniques, here called WHLPCRS, where the operator is exclusively designed in the space-time domain. The second method is a velocity filter (ω-k), applied after the CRS stacking, where the operator is exclusively designed in the frequency-wavenumber domain. The identification of the multiples is performed on the zero-offset (ZO) section simulated by the CRS stacking, using the periodicity between primaries and its multiples. The wavefront attributes, obtained through CRS stacking, are inserted on movable space-time windows, used to calculate the WHLP-CRS operator. The ω-k filter calculations are performed in the frequency-wavenumber domain, where the events of interest are selected for cutting or passage. The ω-k filter is classified as a cut-off filter, with amplitude alteration and preservation of phase, the limits of it are imposed by a space-time sampling. In practical aspects we conclude that for the case of multiples, separated events on the x-t domain are not necessarily separated on the ω-k domain, which raise difficulties in the designing of a ω-k operator with a similar performance when compared to the x-t operator.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-28T15:36:55Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Modelagem1d2.5dDados.pdf: 6036247 bytes, checksum: 09857d990695852fd78d1d492571ec9e (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:23:28Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-27T14:13:55Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Modelagem1d2.5dDados.pdf: 6036247 bytes, checksum: 09857d990695852fd78d1d492571ec9e (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-10T16:36:35Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Modelagem1d2.5dDados.pdf: 6036247 bytes, checksum: 09857d990695852fd78d1d492571ec9e (MD5) / Made available in DSpace on 2014-09-10T16:36:35Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Modelagem1d2.5dDados.pdf: 6036247 bytes, checksum: 09857d990695852fd78d1d492571ec9e (MD5) Previous issue date: 2014 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / INCT/GP - Instituto Nacional de Ciência e Tecnologia de Geofísica do Petróleo / Neste trabalho apresentamos a solução do campo eletromagnético gerado por um dipolo elétrico horizontal em meios transversalmente isotrópicos com eixo de simetria vertical (TIV) e com eixo de simetria inclinado (TII). Para modelos unidimensionais, o campo eletromagnético foi obtido por duas metodologias distintas: (1) solução semi-analítica das equações de Maxwell com auxílio de potenciais vetores no caso TIV e (2) em modelos com anisotropia transversal inclinada o campo eletromagnético foi separado em primário e secundário, e então, o campo secundário foi calculado pelo método de elementos finitos no domínio (kx, ky, z) da transformada de Fourier. Para estruturas bidimensionais, foi aplicada a mesma metodologia usado nos modelos TII unidimensionais, onde o campo secundário foi calculado pelo método de elementos finitos no domínio (x, ky, z), da transformada de Fourier, com a utilização de malhas não estruturadas para discretização dos modelos. Estas respostas foram usados para avaliar os efeitos da anisotropia elétrica nos dados CSEM marinho 1D e 2,5D. / In this work I present the solution to the electromagnetic field generated by a horizontal electric dipole in transversally isotropic media with vertical (TIV) as well as inclined (TII) symmetry axis. In one-dimensional models the electromagnetic field was obtained with two distinct methods: (1) For the TIV case, I have written a semi-analytical solution to the Maxwell’s equations, by using a vector potential formulation; (2) For the TII case, the field was represented as the composition of primary and secondary fields, where primary fields are those found in an underlying isotropic layered medium, and the secondary field is calculated numerically via the finite element method in the spatial Fourier transform domain (kx, ky, z). This last methodology was also used to calculate the fields in two-dimensional structures, including inclined anisotropy in any region of the models. In this 2,5D case, I have applied the finite element method in the (x, ky, z) do main. Here I have used unstructured meshes to discretise the media, and parallel programming to solve the linear systems of equations. The responses were used to study the effects of electrical anisotropy in marine CSEM data.