3D resistivity mapping of airborne EM data

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 1896-1905 ◽  
Author(s):  
Zhiyi Zhang

A 3D resistivity mapping technique has been developed to provide fast estimates of resistivity distributions in airborne electromagnetic surveys. This proposed 3D mapping method consists of an approximate 3D linear inverse operator and a generalized subspace solver. The 3D inverse operator can be generated using any forward approximation that is linear in resistivity. The generalized subspace method is an alternative to the conjugate gradient method, and it reduces the original large linear system of equations to a much smaller but nonlinear one that is solved iteratively. The major benefit of using generalized subspace methods is that subspace vectors can be built based upon physical principles such as skin and investigation depths. Since the 3D mapping is a linear inverse problem, no iteration, and thus no forward modeling nor sensitivity updating, is needed. The 3D resistivity‐mapping technique can be used directly to estimate 3D resistivity distribution or to provide a model update during an intermediate iteration in a nonlinear 3D inversion. Synthetic and field data examples indicate that the 3D mapping can provide quantitative information about the resistivity and spatial distributions of the 3D targets.

2019 ◽  
Vol 92 (1102) ◽  
pp. 20190109 ◽  
Author(s):  
Fengdan Wang ◽  
Haiping Zhang ◽  
Chanyuan Wu ◽  
Qian Wang ◽  
Bo Hou ◽  
...  

Objective: Dermatomyositis (DM) and polymyositis (PM) make up the largest group of potentially treatable myopathies and require early diagnosis. This study investigates whether the edema of thigh muscles in DM/PM can be quantitatively assessed by a novel accelerated T2 mapping technique—GRAPPATINI. Methods: Three conventional MR sequences and GRAPPATINI accelerated T2 mapping of bilateral thighs from 20 patients (7 DM and 13 PM) and 10 healthy volunteers were prospectively carried out on a 3 T MR scanner. Afterwards, T2 values of 477 thigh muscles from the patients and the healthy controls were manually measured. In addition, the correlations between T2 values and serum muscle enzymes in patients were also analyzed. Results: The new GRAPPATINI technique made quantitative T2 mapping of bilateral thighs feasible with a scanning time of only 2 min 18 s. Moreover, GRAPPATINI-generated T2 values of muscles from patients were markedly higher than those from healthy subjects (p < 0.001). GRAPPATINI accelerated T2 mapping appeared a more sensitive technique in that some DM/PM muscles appearing normal per conventional MRI had increased T2 relaxation time. Furthermore, GRAPPATINI-generated T2 values of DM/PM thigh muscles positively correlated with serum enzyme levels (p < 0.001), which reflected the severity of myopathy. Conclusion: GRAPPATINI can significantly shorten acquisition time of T2 mapping and may potentially be applied clinically in DM and PM. Advances in knowledge: GRAPPATINI acceleration makes T2 mapping feasible in clinical practice in providing quantitative information regarding thigh muscle inflammation in DM and PM.


Geophysics ◽  
2007 ◽  
Vol 72 (4) ◽  
pp. F189-F195 ◽  
Author(s):  
Changchun Yin ◽  
Greg Hodges

The traditional algorithms for airborne electromagnetic (EM) inversion, e.g., the Marquardt-Levenberg method, generally run only a downhill search. Consequently, the model solutions are strongly dependent on the starting model and are easily trapped in local minima. Simulated annealing (SA) starts from the Boltzmann distribution and runs both downhill and uphill searches, rendering the searching process to easily jump out of local minima and converge to a global minimum. In the SA process, the calculation of Jacobian derivatives can be avoided because no preferred searching direction is required as in the case of the traditional algorithms. We apply SA technology for airborne EM inversion by comparing the inversion with a thermodynamic process, and we discuss specifically the SA procedure with respect to model configuration, random walk for model updates, objective function, and annealing schedule. We demonstrate the SA flexibility for starting models by allowing the model parameters to vary in a large range (far away from the true model). Further, we choose a temperature-dependent random walk for model updates and an exponential cooling schedule for the SA searching process. The initial temperature for the SA cooling scheme is chosen differently for different model parameters according to their resolvabilities. We examine the effectiveness of the algorithm for airborne EM by inverting both theoretical and survey data and by comparing the results with those from the traditional algorithms.


Geophysics ◽  
1981 ◽  
Vol 46 (9) ◽  
pp. 1278-1290 ◽  
Author(s):  
L. E. Reed

In June 1974, a diamond drill operated for Selco Mining Corp. intersected zinc‐copper sulfides in Brouillan Township in northwestern Québec. To date, two bodies have been outlined. These bodies were discovered during a ground follow‐up of a Mark VI Input® electromagnetic (EM) survey. The Input survey covered an area selected on the basis of regional geology and local outcrops of acid volcanic rocks. Conductors were identified that appeared to be associated with potentially favorable geology. They were selected for ground follow‐up. One was the discovery zone. The airborne responses over the zone were less encouraging than those often observed over highly conductive massive sulfides. The low apparent conductivity‐thickness (5 mhos) was suggestive of conductive overburden. However, the character of the profiles suggested a bedrock source. Ground geophysical confirmation identified a drill target. Subsequent to the discovery, more intensive geophysical surveys, both ground and airborne, were carried out. The best EM response suggested a confined source within a much larger mineralized halo. Weaker ground EM response from the halo correlated with the early channel response of the Input system. An airborne EM survey conducted in 1958 over the same area identified both conductive zones. However, they were not followed up. Only with later advances in exploration philosophy, geologic appreciation, and instrumentation were the conductive zones recognized as viable exploration targets.


Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 492-500 ◽  
Author(s):  
James E. Reid ◽  
James C. Macnae

When a confined conductive target embedded in a conductive host is energized by an electromagnetic (EM) source, current flow in the target comes from both direct induction of vortex currents and current channeling. At the resistive limit, a modified magnetometric resistivity integral equation method can be used to rapidly model the current channeling component of the response of a thin-plate target energized by an airborne EM transmitter. For towed-bird transmitter–receiver geometries, the airborne EM anomalies of near-surface, weakly conductive features of large strike extent may be almost entirely attributable to current channeling. However, many targets in contact with a conductive host respond both inductively and galvanically to an airborne EM system. In such cases, the total resistive-limit response of the target is complicated and is not the superposition of the purely inductive and purely galvanic resistive-limit profiles. Numerical model experiments demonstrate that while current channeling increases the width of the resistive-limit airborne EM anomaly of a wide horizontal plate target, it does not necessarily increase the peak anomaly amplitude.


2020 ◽  
Vol 9 (6) ◽  
pp. 258-267
Author(s):  
Xiang Yao ◽  
Kaihua Zhou ◽  
Bin Lv ◽  
Lei Wang ◽  
Jun Xie ◽  
...  

Aims Tibial plateau fractures (TPFs) are complex injuries around the knee caused by high- or low-energy trauma. In the present study, we aimed to define the distribution and frequency of TPF lines using a 3D mapping technique and analyze the rationalization of divisions employed by frequently used classifications. Methods In total, 759 adult patients with 766 affected knees were retrospectively reviewed. The TPF fragments on CT were multiplanar reconstructed, and virtually reduced to match a 3D model of the proximal tibia. 3D heat mapping was subsequently created by graphically superimposing all fracture lines onto a tibia template. Results The cohort included 405 (53.4%) cases with left knee injuries, 347 (45.7%) cases with right knee injuries, and seven (0.9%) cases with bilateral injuries. On mapping, the hot zones of the fracture lines were mainly concentrated around the anterior cruciate ligament insertion, posterior cruciate ligament insertion, and the inner part of the lateral condyle that extended to the junctional zone between Gerdy’s tubercle and the tibial tubercle. Moreover, the cold zones were scattered in the posteromedial fragment, superior tibiofibular syndesmosis, Gerdy’s tubercle, and tibial tubercle. TPFs with different Orthopaedic Trauma Association/AO Foundation (OTA/AO) subtypes showed peculiar characteristics. Conclusion TPFs occurred more frequently in the lateral and intermedial column than in the medial column. Fracture lines of tibial plateau occur frequently in the transition zone with marked changes in cortical thickness. According to 3D mapping, the four-column and nine-segment classification had a high degree of matching as compared to the frequently used classifications. Cite this article: Bone Joint Res 2020;9(6):258–267.


IEEE Access ◽  
2019 ◽  
Vol 7 ◽  
pp. 73593-73601 ◽  
Author(s):  
Binghua Guo ◽  
Hongyue Dai ◽  
Zhonghua Li ◽  
Wei Huang

2010 ◽  
Vol 3 (1) ◽  
pp. 13-41 ◽  
Author(s):  
T. J. Reerink ◽  
M. A. Kliphuis ◽  
R. S. W. van de Wal

Abstract. Here, we present a mapping method OBLIMAP, which projects and interpolates fields like surface temperature, surface mass balance, and surface height between a geographical based coordinate system of a General Circulation Model (GCM) and a rectangular based Ice Model (IM). We derive an oblique stereographic projection and its inverse, which holds for any area at the Earth's surface, and which can be combined with two different interpolation methods. The first one is suited to interpolate the projected fields of a coarse GCM grid on a fine meshed IM grid. The second one is appropriate for the opposite case. Both grids are allowed to be arbitrary and irregularly spaced. Therefore the OBLIMAP technique is suitable for any GCM-IM combination. After a first scan of the GCM grid coordinates and the specification of the IM grid, fast mapping of various fields is possible. To and fro (GCM-IM-GCM) mapping tests with the Climate Community System Model (CCSM) at T42 resolution (~313 km) and the Regional Atmospheric Climate Model (RACMO) at ~11 km and ~55 km, show average temperature differences of less than 0.1 K with small standard deviations. OBLIMAP, available at GMD, is an accurate, robust and well-documented mapping method for coupling an IM with a GCM or to map state of the art initial and forcing fields available at geographical coordinates to any local IM grid with an optimal centered oblique projection. Currently, the oblique stereographic and the oblique Lambert azimuthal equal-area projections for both the sphere and the ellipsoid are implemented in OBLIMAP.


Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 484-491 ◽  
Author(s):  
Gregory A. Newman ◽  
David L. Alumbaugh

A 3‐D finite‐difference solution is implemented for simulating induction log responses in the quasi‐static limit that include the wellbore and bedding that exhibits transverse anisotropy. The finite‐difference code uses a staggered grid to approximate a vector equation for the electric field. The resulting linear system of equations is solved to a predetermined error level using iterative Krylov subspace methods. To accelerate the solution at low induction numbers (LINs), a new preconditioner is developed. This new preconditioner splits the electric field into curl‐free and divergence‐free projections, which allows for the construction of an approximate inverse operator. Test examples show up to an order of magnitude increase in speed compared to a simple Jacobi preconditioner. Comparisons with analytical and mode matching solutions demonstrate the accuracy of the algorithm.


Pomorstvo ◽  
2021 ◽  
Vol 35 (2) ◽  
pp. 353-364
Author(s):  
A. Mertcan Yasa ◽  
Abdi Kükner

In this paper, a method to predict slamming pressures and pressure distribution at the time of water entry for 2D sections is presented. The mathematical model is based on the Schwarz-Christoffel conformal mapping method. This conformal mapping technique has been used to calculate slamming loads during water entry. The pile-up of water during motion is also considered and an alternative pile-up coefficient is assumed against Wagner’s generalized method. A simplified and accurate method is presented, which does not include non-linear terms and jet flow in the calculated pressure distribution on monotonically increasing 2D sections like wedge shapes. Comparison with real ship sections has been done to show accuracy of the results. Finally, a simple yet powerful method is obtained to aid the initial design stage of ships.


Sign in / Sign up

Export Citation Format

Share Document