@article { author = {Ghari, Hosseinali and Ghorbani, Ahmad and Ansari, Abdol hamid}, title = {Inversion of spectral induced polarization data for Cole–Cole parameters using simulated annealing algorithm}, journal = {Iranian Journal of Geophysics}, volume = {7}, number = {3}, pages = {78-92}, year = {2013}, publisher = {Iranian Geophysical Society}, issn = {2008-0336}, eissn = {2783-168X}, doi = {}, abstract = {Spectral induced polarization (SIP) method is an exploration geophysical method. The induced polarization (IP) method has been used extensively in mine prospecting and increasingly in environmental investigations because IP measurements are very sensitive to the low frequency capacitive properties of rocks and soils. Many different models have been proposed for the description of the dispersive behaviour of the IP. However, the most widely used model is the Cole-Cole model. This model describes the resistivity dispersion observed in field data from areas with metallic mineral content. It is also used to estimate various subsurface properties of nonmetallic soil and rocks in IP frequency domain investigations (SIP). A multiple Cole-Cole model is typically a more general and proper model than a single Cole-Cole model for describing IP data with various dispersion ranges caused either by multiple-length scales in sediments or by coupling effects in the IP measurements. The Cole-Cole model parameters are widely used to interpret both time- and frequency-domain induced polarization data. Among many studies in which the Cole-Cole parameters are estimated from SIP measurements on soils and rocks, a majority of them use least squares (deterministic) methods. In this work, we have developed a Bayesian method with simulated annealing sampling algorithm to invert single and double Cole-Cole parameters from SIP data. We have also reproduced the Gibbs sampling algorithm developed by Chen et al in 2008. The Bayesian approach consists of propagating the information provided by measurements throughout the model and combining this information with a priori knowledge of the data. A major limitation towards a more widespread implementation of Bayesian approaches is that obtaining the posterior distribution often requires the integration of high-dimensional functions. This can be computationally very difficult, but several approaches have been proposed which fall short of direct integration functions. Here, we focused on Markov Chain Monte Carlo (MCMC) methods, which attempt to simulate direct draws from some complex distribution of interest. The simulated annealing and Gibbs sampling are particular MCMC methods widely applicable to a broad class of Bayesian problems and have sparked a major increase in the application of Bayesian analysis. This interest is likely to continue in the future. In the Gibbs sampling method, the acceptance probability is 1; that is, Gibbs sampler candidates are always accepted. On the other hand, the acceptance probability is between 0 and 1 for the simulated annealing method.The MCMC-based inversion method provides extensive global information on unknown parameters such as the marginal probability distribution functions from which we can obtain better estimates and tighter uncertainty bounds of the parameters compared to deterministic methods. We compared performance of the the simulated annealing method with the Gibbs sampling algorithm method through an inversion of synthetic and real data. Both methods provide a global approach to invert SIP data for the Cole-Cole parameters in which the obtained estimates are independent of the initial values. In addition, this approach has the ability to avoid getting trapped in a local optimum. For synthetic data with 10% random noise and 50 second run-time, the results of the simulated annealing are more similar to the results of the Gibbs sampling method. The latter needs 6 minutes of time. Also, in the case of the real data obtained in laboratory, the simulated annealing method has provided a more suitable fit in comparison with the Gibbs sampling method. As a result, for a problem with more parameters, the required time is increased in the Gibbs sampling method, while the simulated annealing method reduces the required time.  }, keywords = {Spectral induced polarization,Cole-Cole model,Inversion,Bayesian,Simulated Annealing,Sampling}, title_fa = {معکوس‌سازی پارامترهای مدل کول-کول با استفاده از الگوریتم شبیه‌سازی بازپخت برای داده‌های قطبش القایی طیفی (SIP)}, abstract_fa = {قطبش القایی طیفی (SIP) شاخه‌ای از روش‌‌های ژئوفیزیکی است که به‌طور گسترده در پی‌جویی‌های معدنی و زیست‌محیطی مورد استفاده قرار می‌گیرد. برای تفسیر و وارون‌سازی داده‌های قطبش القایی طیفی با مدل کول-کول، چهار پارامتر r0، m، t و c بازیابی می‌شود. تحقیقات زیادی در زمینه بازیابی پارامترهای کول-کول از داده‌های قطبش القایی طیفی صورت گرفته است که اکثراً براساس روش‌های کمترین مربعات خطا استوار است. در پژوهش حاضر از توسعه استنباط بیزی (Bayesian) برای برآورد پارامترهای دو کول-کول استفاده شده است. در این استنباط دو روشِ مرسوم برای نمونه‌برداری از تابع توزیع پسین وجود دارد: 1.روش شبیه‌سازی بازپخت‌ (Simulated Annealing (SA)) 2. روش نمونه‌برداری گیبس(Gibss Sampling (GS)) . در مقاله حاضر، ابتدا نمونه‌برداری با استفاده از روش شبیه‌سازی بازپخت صورت گرفته است‌. این الگوریتم در محیط نرم‌افزار مَت‌لَب برنامه‌نویسی شده و با رایانه‌ای به بسامد 2.53GHz و حافظه‌ زنده 4Gb اجرا شده است. سپس الگوریتم نمونه‌برداری گیبس، که چن و همکاران آن را گسترش داده‌اند، بازسازی شده است. همچنین، نتایج وارون‌سازی حاصل از دو الگوریتم برای داده‌های قطبش القایی طیفی مصنوعی و واقعی به‌دست آمده در آزمایشگاه، با هم مقایسه شده‌ است.     نتایج نشان می‌دهد که هر دو روش یک رهیافت کلی برای وارون‌سازی پارامترهای مدل کول-کول از داده‌های قطبش القایی طیفی تامین می‌کنند، در یافتن کمینه‌‌ واقعی موفق بوده‌اند، درگیر کمینه‌‌‌های محلی نمی‌شوند و برآورد‌های به‌دست آمده از آنها مستقل از مقادیر اولیه‌ پارامترها است. برای داده‌های مصنوعی با نوفه‌‌ تصادفی 10% و زمان 50 ثانیه، نتایج شبیه‌سازی بازپخت نسبت به نمونه‌برداری گیبس به داده‌های واقعی نزدیک‌تر است، درحالی‌که روش نمونه‌برداری گیبس برای رسیدن به چنین تقریبی به 6 دقیقه زمان نیاز دارد. برای داده‌های آزمایشگاهی نیز روش شبیه‌سازی بازپخت نسبت به روش نمونه‌برداری گیبس در زمان کمتر، برازش مناسب‌تری بر داده‌ها به‌دست می‌دهد. در نتیجه برای مسئله‌ای با پارامترهای بیشتر، زمان صرف شده در روش نمونه‌برداری گیبس به مراتب افزایش می‌یابد، درصورتی‌که روش شبیه‌سازی بازپخت، این زمان را به حداقل می‌رساند.    }, keywords_fa = {قطبش القایی طیفی,مدل کول-کول,معکوس‌سازی,روش بیزی,شبیه‌سازی بازپخت,نمونه‌برداری}, url = {https://www.ijgeophysics.ir/article_40359.html}, eprint = {https://www.ijgeophysics.ir/article_40359_bf981a9e631d019fd093230bb971bcf8.pdf} }