روش سیگنال تحلیلی در برآورد پارامتر‌‌های منشأ مغناطیسی ساختار‌‌های دوبُعدی

نویسندگان

1 استادیار، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران، ایران

2 دانشجوی کارشناسی ارشد ژئومغناطیس، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران، ایران

3 دانشیار، گروه فیزیک زمین، مؤسسة ژئوفیزیک دانشگاه تهران و قطب علمی مهندسی نقشه‌برداری و مقابله با سوانح طبیعی، تهران، ایران

چکیده

روش سیگنال تحلیلی را برای تفسیر داده‌‌های بی‌هنجاری مغناطیسی و گرانی‌سنجی به کار می‌‌برند. مزیت استفاده از روش سیگنال تحلیلی برای تفسیر بی‌هنجاری‌‌های مغناطیسی، مستقل بودن شکل دامنه سیگنال تحلیلی ازپارامترهای جهت‌دار مانند جهت مغناطیدگی و شیب منشأ مغناطیسی است. دامنه سیگنال تحلیلی یک تابع متقارن زنگوله‌‌ای شکل است که بیشینه‌‌ آن دقیقاً بالای گوشه بی‌هنجاری قرار گرفته است و پهنای منحنی دامنه آن با عمق سطح رویه بی‌هنجاری نسبت مستقیم دارد. از این ویژگی سیگنال تحلیلی در برآورد کردن پارامترهای منشأ مغناطیسی، مانند عمق و ضخامت استفاده می‌‌شود. در این مقاله در محیط نرم‌افزار مت‌لَب MATLAB به منظور محاسبه سیگنال تحلیلی بی‌هنجاری‌‌های دوبُعدی، تعیین پارامترهای هندسی منشأ و مدل‌سازی مستقیم برنامه‌نویسی شده و برای بررسی کارایی برنامه‌‌های نوشته شده ، ابتدا آنها را روی داده‌‌های مصنوعی همراه با نوفه اِعمال می‌‌کنیم. سپس از این روش برای تفسیر داده‌‌های صحرایی منطقه سرخ دیزج زنجان استفاده کرده‌ایم. به منظور برآورد عمق و ضخامت منشأ مغناطیسی منطقه مورد بررسی، در نقشه مغناطیسی از 5 نیم‌رخ‌‌، عمود بر راستای امتداد بی‌هنجاری مغناطیسی استفاده کرده‌‌ایم. این نیم‌رخ‌‌‌‌ها در طول جغرافیایی m 305986 تا m306176 و عرض جغرافیایی m4051813 تا m4052013 بر مبنای واحد اندازه‌گیری UTM در نظر گرفته شده‌اند. که نتایج حاصل از این نیم‌رخ‌‌‌‌ها، برآورد عمق و ضخامت بی‌هنجاری است. در پایان نتایج به دست آمده از روش سیگنال تحلیلی را با نتایج حاصل از روش اویلر مقایسه کرده‌ایم.

کلیدواژه‌ها


عنوان مقاله [English]

Analytic signal method to estimate the magnetic source parameters of 2D structures

نویسندگان [English]

  • Behrouz Oskooi 1
  • Ali Karimi Kalaye 2
  • Vahid Ebrahimzadeh Ardestani 3
چکیده [English]

The analytic signal method is usually used to interpret the magnetic and gravity anomalies. The main advantage of using analytic signal method in interpretation of magnetic source parameters from magnetic anomalies is that the amplitude dip of the analytic signal is independent of directional parameters like magnetization and the dip of the source. Using the coefficient of the amplitude to estimate the depth can make the interpretation of the magnetic data easier. By this method the location depth and angle of the anomalies and even the dip thickness and the direction of its extension can be estimated. The amplitude of the signal is a symmetric bell-shaped function that the maximum of the curve is located exactly above the edge of the source and its width is proportional to the depth of the top of the magnetic source. This analytic signal characteristic is used to estimate the magnetic source parameters like depth location and the thickness.
We use these characteristics of the analytic signal function to estimate the depth, horizontal location and thickness of the magnetic sources. In this study, we use MATLAB in order to write codes for calculating the analytic signal of two dimensional anomalies, determine the geometrical parameters of the source and forward modeling of magnetic anomalies.
For evaluating the efficiency of the codes, first we applied them on the synthetic noise contaminated data. We then used analytic signal method to interpret the field magnetic data from Sorkh-Dizej, Zanjan. To estimate the depth and thickness of the magnetic source, we used 5 profiles perpendicular to the strike of the magnetic anomaly. These profiles are located in the longitude from 305986 m to 306176 m and latitude from 4051813 m to 4052013 m. Finally, we compared the results of analytic signal with the results of Euler method.
To consider the fact that the horizontal and vertical differentials are each other's Hilbert conversion, in each step of 2D analytic signal differentiation of the anomalies of 3D sources can be interpreted by extending 2D analytic signals.
The function of the analytic signal can be calculated easily in frequency region that its real and imaginary parts are horizontal and vertical differentials of the magnetic field in sequence. Another advantage of using analytic signal differential is the interpretation of magnetic and gravity data in common projects by Poisson equation.
The previous data in analytic signal method is not essential while in other methods the interpretation of the magnetic data needs some data such as the geology of the structure and other parameters. For example in Euler deconvolution method, that it is one of the depth estimation methods, knowing the structural index is vital and difficult.
In Sorkhdizaj in Zanjanin an iron mine with apatite content the results of analyzing data by analytic signal method show that the structure of the anomalies is too similar to a thin dike so that the method recognizes the thin dike with high accuracy and low uncertainty. Besides, it is shown that the depth of the top of the anomalies are changing from 10 to 55 m so that the results are too close to the ones obtained by Euler method calculated by Geosoft software. Also dikes with thickness less than 1 m to about 10 m have been estimated.

کلیدواژه‌ها [English]

  • Analytic signal
  • Depth
  • Horizontal location
  • Magnetization
  • Thickness