Two‐dimensional terrain correction in magnetotelluric surveys

Geophysics ◽  
1988 ◽  
Vol 53 (6) ◽  
pp. 854-862 ◽  
Author(s):  
Michel Chouteau ◽  
Karl Bouchard

Field distortions caused by topography hamper the interpretation of magnetotelluric (MT) data. Topographic features that can be simulated by two‐dimensional models seriously affect the H-polarization results. A technique to reduce those effects in MT data uses a finite‐element program to compute correction coefficients. After correction, the resulting data can be interpreted as if they were obtained over a flat surface and depended only on the subsurface structure. The technique is applied to four examples representative of MT survey targets in high‐relief terrain. Results indicate that terrain correction removes the misleading topographic anomalies and improves the quality of subsurface interpretation in regions where the surface relief is two‐dimensional. The correction technique yields some geometrical distortion of the original subsurface structure, but the distortion is usually of small importance. In practice, telluric dipoles of realistic length do not smooth out topographic effects having wavelengths longer than the telluric dipole. A practical procedure derived from the proposed technique allows terrain correction when the relief is approximately two‐dimensional.

2004 ◽  
Vol 120 ◽  
pp. 727-735
Author(s):  
F. Frerichs ◽  
Th. Lübben ◽  
U. Fritsching ◽  
H. Lohner ◽  
A. Rocha ◽  
...  

The prediction of mechanical behaviour of specimen during heat treatment by means of numerical simulation requires numerous modules e.g. for heat transfer and mechanical behaviour. The quality of predictions depend on the quality of the applied models within the modules. In this paper the strain hardening model used in the mechanical module will be investigated. For simulation of mechanical behaviour during gas quenching it is first of all necessary to calculate the interaction between gas and specimen. Using simulated flow field and temperature distribution within the gas, the heat transfer coefficient is calculated from computational fluid dynamics. The cooling and further the mechanical behaviour e.g. residual stresses and distortion of the specimen are simulated by a commercial Finite Element program. To investigate strain hardening it is helpful to choose in a first step a material that will not show phase transformations due to heat treatment. Therefore simulation of mechanical behaviour of austenitic cylinders (SAE30300) is investigated. The required thermo-physical properties such as thermal conductivity, density, and specific heat are taken from literature. With the exception of Poisson’s ratio the mechanical properties are measured and calculated by own investigations. For description of the temperature dependent stress strain curves the Ramberg-Osgood model is used. The simulated results are compared with experimental data in order to decide which model better describes the mechanical response, whether the kinematic or isotropic strain hardening.


Author(s):  
H. Al Nageim ◽  
F. Mohammad ◽  
Lewis Lesley

The finite element method is used to determine the response of a new lightweight rail track system (LR55) to various loading environments. To calculate the vertical displacement of the rail track system and the vertical pressure in the sub-base and sub-grade layers due to wheel loads and loads exerted by transversely passing vehicles across the track, the global stiffness matrix of the structure is determined. This is done by using one and two dimensional finite element programs. In the two dimensional finite element analysis, linear isoparametric elements with 4-node quadrilateral and 3-node triangular shapes in the discretised mesh of the whole structure are used, also non-homogeneous materials with isotropic of linear elastic behaviour are assumed for all the components forming the track system and surrounding media. The results of the one dimensional finite element program are compared with those predicted from analytical approach in order to validate the finite element program developed. From the various examples presented the LR55 proves to withstand the main-line railway loading and the vertical pressure distribution in the sub-base and sub-grade of the pavement are within the allowable limit.


Author(s):  
Jyoti Agarwal ◽  
Ankush Chaudhary

For dynamic analysis, it is required to provide viscous boundaries in PLAXIS to reduce the boundary effects and to prevent the reflection of waves from boundaries. So, a study has been carried out to compare the various methods of providing silent boundaries and to see the effectiveness of viscous boundaries used in PLAXIS. In this work, three methods of providing silent boundaries, which are viscous boundaries, local damping, and extended boundary, are analyzed using a 2D finite element program in FORTRAN by considering the simple problem of a two-dimensional vertical bar. Parameters, such as, normal stress at the bottom, vertical displacement at the top, potential energy, kinetic energy, strain energy, and total energy of bar are determined with and without using the above three methods of providing silent boundaries. Results are compared using graphs. It was observed that standard viscous boundaries are not much effective for static analysis but most effective for dynamic analysis.


1988 ◽  
Vol 110 (2) ◽  
pp. 68-76
Author(s):  
E. Evgin ◽  
S. Lord ◽  
I. Konuk

The results of the two and three-dimensional, nonlinear analyses of a caisson retained island in the Beaufort Sea are presented. A hyperbolic stress-strain relation for soils has been implemented in a general purpose finite element program, ADINA, to carry out the analysis. Based on the comparisons of the results, conclusions were drawn about whether a two-dimensional analysis can be safely used to replace a three-dimensional analysis.


Geophysics ◽  
1980 ◽  
Vol 45 (1) ◽  
pp. 75-93 ◽  
Author(s):  
Richard C. Fox ◽  
Gerald W. Hohmann ◽  
Terry J. Killpack ◽  
Luiz Rijo

We have made a systematic study of dipole‐dipole apparent resistivity anomalies due to topography and of the effect of irregular terrain on induced‐polarization (IP) anomalies, using a two‐dimensional (2-D), finite‐element computer program. A valley produces a central apparent resistivity low in the resistivity pseudosection, flanked by zones of higher apparent resistivity. A ridge produces just the opposite anomaly pattern—a central high flanked by lows. A slope generates an apparent resistivity low at its base and a high at its top. Topographic effects are important for slope angles of 10 degrees or more and for slope lengths of one dipole‐length or greater. The IP response of a homogeneous earth is not affected by topography. However, irregular terrain does affect the observed IP response of a polarizable body due to variations in the distance between the electrodes and the body. These terrain‐induced anomalies can lead to erroneous interpretations unless topography is included in numerical modeling. A field case demonstrates the importance of including topography, where it is significant, in interpretation models. A technique for correcting apparent resistivity for topographic effects uses the finite‐element program to compute correction factors.


Geophysics ◽  
1986 ◽  
Vol 51 (11) ◽  
pp. 2131-2144 ◽  
Author(s):  
Philip E. Wannamaker ◽  
John A. Stodt ◽  
Luis Rijo

We simulate the magnetotelluric response to two‐dimensional earth topography using finite elements. Linear interpolation of the secondary field parallel to strike over triangular elements allows accurate modeling of inclined resistivity boundaries, including topographic surfaces. To avoid discontinuities in field derivatives or resistivity, care must be taken that the nodal values of the field parallal to strike used to obtain the auxiliary secondary fields are kept within uniform earth media. The nodal locations may be shifted, but the derivatives still are evaluated at the field points of interest. Correct values may be returned at gentle breaks in slope as well as along straight surfaces. The finite‐element program is verified by comparison with the analytic transverse magnetic response of a hemicylindrical depression and with Rayleigh scattering and transmission surface results for transverse electric and transverse magnetic polarization. Agreement with the other methods generally is excellent, with the exception of some results of the transmission surface technique (especially the transverse magnetic mode). A result presented shows that placing the H-field sensors horizontally reduces topographic anomalies compared to locating sensors parallel to the slope. Moreover, if earth resistivity increases with elevation, the apparent resistivity is relatively nonanomalous near the base of the topography.


Sign in / Sign up

Export Citation Format

Share Document