跳至主要内容

A MATLAB-Based Numerical and GUI Implementation of Cross-Gradients Joint Inversion of Gravity and Magnetic Data

Read  full  paper  at:
http://www.scirp.org/journal/PaperInformation.aspx?PaperID=54101#.VO17AyzQrzE

Author(s)  
 
ABSTRACT
The cross-gradients joint inversion technique has been applied to multiple geophysical data with a significant improvement on compatibility, but its numerical implementation for practical use is rarely discussed in the literature. We present a MATLAB-based three-dimensional cross-gradients joint inversion program with application to gravity and magnetic data. The input and output information was examined with care to create a rational, independent design of a graphical user interface (GUI) and computing kernel. For 3D visualization and data file operations, UBC-GIF tools are invoked using a series of I/O functions. Some key issues regarding the iterative joint inversion algorithm are also discussed: for instance, the forward difference of cross gradients, and matrix pseudo inverse computation. A synthetic example is employed to illustrate the whole process. Joint and separate inversions can be performed flexibly by switching the inversion mode. The resulting density model and susceptibility model demonstrate the correctness of the proposed program.
 
Cite this paper
Zhou, J. , Zhang, X. and Xiu, C. (2015) A MATLAB-Based Numerical and GUI Implementation of Cross-Gradients Joint Inversion of Gravity and Magnetic Data. Journal of Software Engineering and Applications, 8, 93-101. doi: 10.4236/jsea.2015.82010.
 
References
[1]Zeyen, H. and Pous, J. (1993) 3-D Joint Inversion of Magnetic and Gravimetric Data with a Priori Information. Geophysical Journal International, 112, 244-256.
http://dx.doi.org/10.1111/j.1365-246X.1993.tb01452.x
 
[2]Gallardo, L.A., Perez-Flores, M.A. and Gomez-Trevino, E. (2005) Refinement of Three-Dimensional Multilayer Models of Basins and Crustal Environments by Inversion of Gravity and Magnetic Data. Tectonophysics, 397, 37-54.http://dx.doi.org/10.1016/j.tecto.2004.10.010
 
[3]Bosch, M., Meza, R. and Jimenez, R. (2006) Joint Gravity and Magnetic Inversion in 3D Using Monte Carlo Methods. Geophysics, 71, G163-G156. http://dx.doi.org/10.1190/1.2209952
 
[4]Pilkington, M. (2006) Joint Inversion of Gravity and Magnetic Data for Two-Layer Models. Geophysics, 71, L35-L42.http://dx.doi.org/10.1190/1.2194514
 
[5]Shamsipour, P., Marcotte, D. and Chouteau, M. (2012) 3D Stochastic Joint Inversion of Gravity and Magnetic Data. Journal of Applied Geophysics, 79, 27-37.
http://dx.doi.org/10.1016/j.jappgeo.2011.12.012
 
[6]Zhdanov, M.S., Gribenko, A. and Wilson, G. (2012) Generalized Joint Inversion of Multimodal Geo-physical Data Using Gramian Constraints. Geophysical Research Letters, 39, L09301.
http://dx.doi.org/10.1029/2012GL051233
 
[7]Williams, N.C. (2008) Geologically-Constrained UBC-GIF Gravity and Magnetic Inversions with Ex-amples from the Agnew-Wiluna Greenstone Belt, Western Australia. Ph. D. Dissertation. The University of British Columbia, Vancouver.
 
[8]Frogoso, E. and Gallardo, L.A. (2009) Cross-Gradients Joint 3D Inversion with Applications to Gravity and Magnetic Data. Geophysics, 74, L31-L42. http://dx.doi.org/10.1190/1.3119263
 
[9]Gallardo, L.A. and Meju, M.A. (2004) Joint Two-Dimensional DC Resistivity and Seismic Travel Time Inversion with Cross-Gradients Constraints. Journal of Geophysical Research: Solid Earth (1978-2012), 109, B03311. http://dx.doi.org/10.1029/2003JB002716
 
[10]Oldenburg, D.W. and Jones, F.H.M. (2007) Inversion for Applied Geophysics: Learning Resources about Geophysical Inversion. The University of British Columbia, Vancouver.
http://www.eos.ubc.ca/ubcgif/iag/index.htm
 
[11]Özgü Arisoy, M. and Dikmen, ü. (2011) Potensoft: MATLAB-Based Software for Potential Field Data Processing, Modeling and Mapping. Computers & Geosciences, 37, 935-942.
http://dx.doi.org/10.1016/j.cageo.2011.02.008
 
[12]Lelièvre, P.G. and Oldenburg, D.W. (2009) A Comprehensive Study of Including Structural Orientation Information in Geophysical Inversions. Geophysical Journal International, 178, 623-637.
http://dx.doi.org/10.1111/j.1365-246X.2009.04188.x
 
[13]Hahn, B.H. and Valentine, D.T. (2013) Essential Matlab for Engineers and Scientists. 5th Edition, Academic Press, Boston, 129-160. http://dx.doi.org/10.1016/B978-0-12-394398-9.00006-X
 
[14]Eshagh, M. and Abdollahzadeh, M. (2012) Software for Generating Gravity Gradients Using a Geopotential Model Based on an Irregular Semivectorization Algorithm. Computers & Geosciences, 39, 152-160. http://dx.doi.org/10.1016/j.cageo.2011.06.003
 
[15]Pidlisecky, A., Haber, E. and Knight R. (2007) RESINVM3D: A 3D Resistivity Inversion Package. Geo-physics, 72, H1-H10. http://dx.doi.org/10.1190/1.2402499
 
[16]Hansen, P.C. (1994) Regularization Tools: A Matlab Package for Analysis and Solution of Discrete Ill-Posed Problems. Numerical Algorithm, 6, 1-35. http://dx.doi.org/10.1007/BF02149761        eww150225lx
标签：

评论

发表评论

此博客中的热门博文

A Comparison of Methods Used to Determine the Oleic/Linoleic Acid Ratio in Cultivated Peanut (Arachis hypogaea L.)

Cultivated peanut ( Arachis hypogaea L.) is an important oil and food crop. It is also a cheap source of protein, a good source of essential vitamins and minerals, and a component of many food products. The fatty acid composition of peanuts has become increasingly important with the realization that oleic acid content significantly affects the development of rancidity. And oil content of peanuts significantly affects flavor and shelf-life. Early generation screening of breeding lines for high oleic acid content greatly increases the efficiency of developing new peanut varieties. The objective of this study was to compare the accuracy of methods used to classify individual peanut seed as high oleic or not high oleic. Three hundred and seventy-four (374) seeds, spanning twenty-three (23) genotypes varying in oil composition (i.e. high oleic (H) or normal/not high oleic (NH) inclusive of all four peanut market-types (runner, Spanish, Valencia and Virginia), were individually tested ...
发表评论
阅读全文

Location Optimization of a Coal Power Plant to Balance Costs against Plant’s Emission Exposure

Fuel and its delivery cost comprise the biggest expense in coal power plant operations. Delivery of electricity from generation to consumers requires investment in power lines and transmission grids. Placing a coal power plant or multiple power plants near dense population centers can lower transmission costs. If a coalmine is nearby, transportation costs can also be reduced. However, emissions from coal plants play a key role in worsening health crises in many countries. And coal upon combustion produces CO 2 , SO 2 , NO x , CO, Metallic and Particle Matter (PM10 & PM2.5). The presence of these chemical compounds in the atmosphere in close vicinity to humans, livestock, and agriculture carries detrimental health consequences. The goal of the research was to develop a methodology to minimize the public’s exposure to harmful emissions from coal power plants while maintaining minimal operational costs related to electric distribution losses and coal logistics. The objective was...
发表评论
阅读全文

Evaluation of the Safety and Efficacy of Continuous Use of a Home-Use High-Frequency Facial Treatment Appliance

At present, many home-use beauty devices are available in the market. In particular, many products developed for facial treatment use light, e.g., a flash lamp or a light-emitting diode (LED). In this study, the safety of 4 weeks’ continuous use of NEWA TM , a high-frequency facial treatment appliance, every alternate day at home was verified, and its efficacy was evaluated in Japanese individuals with healthy skin aged 30 years or older who complained of sagging of the facial skin.  Transepidermal water loss (TEWL), melanin levels, erythema levels, sebum secretion levels, skin color changes and wrinkle improvement in the facial skin were measured before the appliance began to be used (study baseline), at 2 and 4 weeks after it had begun to be used, and at 2 weeks after completion of the 4-week treatment period (6 weeks from the study baseline). In addition, data obtained by subjective evaluation by the subjects themselves on a visual analog scale (VAS) were also analyzed. Fur...
发表评论
阅读全文