per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
1
21
47668
Research Article
ساختار هندسی جت آفریقا-آسیا در وردسپهر زبرین و پاسخ آن به گرمایش زمین در مدلهای CMIP5
The geometrical structure of the upper-tropospheric African–Asian jet and its response to global warming in the CMIP5 Models
شاهین عالمزاده
sh.alemzadeh@ut.ac.ir
1
فرهنگ احمدی گیوی
ahmadig@ut.ac.ir
2
علیرضا محب الحجه
amoheb@ut.ac.ir
3
دانیل یازجی
daniel.yazgi@ut.ac.ir
4
موسسه ژئوفیزیک دانشگاه تهران، تهران ایران
موسسه ژئوفیزیک دانشگاه تهران، تهران ایران
عضو هیئت علمی موسسه ژئوفیزیک دانشگاه تهران، تهران، ایران
موسسه ژئوفیزیک دانشگاه تهران، تهران، ایران
در این پژوهش، چگونگی شبیهسازی جت جنبحارهای آفریقا-آسیا در وردسپهر زبرین در مقایسه با دادههای بازتحلیل و همچنین شبیهسازی پاسخ این جت به تغییر اقلیم در خروجی مدلهای مرحله پنجم از «پروژه مقایسه متقابل مدلهای جفت شده» موسوم به CMIP5 ارزیابی و تحلیل میشود. دادههای بازتحلیل مورد استفاده در این کار برای بررسی خطای مدلها، ERA-Interim میباشد. این مطالعه بر مبنای شاخصهای هندسی جت شامل عرض جغرافیایی یا بهاختصار عرض جت، پهنا و تندی جت است، و منطقه جت به دو قطاع مجزای آسیا و آفریقا تقسیم میشود. مهمترین ویژگیهای چرخه فصلی جت در دوره گذشته که در میانگین چند مدلی CMIP5 شبیهسازی میشود، عبارت از شیب تغییرات مثبت (منفی) از فصل سرد (گرم) به گرم (سرد) در شاخص عرض جت و برعکس این چرخه در دو شاخص تندی و پهنای جت برای هر دو قطاع است. بر اساس نتایج بررسی، عرض جت در همه فصلها و تندی جت در همه فصلها بهجز بهار در قطاع آسیا بیشتر از آفریقا است. این ویژگیها بهخوبی با دادههای بازتحلیل همخوانی دارد و برای بیشتر فصلها و شاخصها، میانگین چند مدلی خطا استوار نمیباشد. خطاهای استوار برای عرض جت در میانگین چند مدلی، به تابستان و زمستان قطاع آفریقا (بهترتیب 86/1- و 90/0- درجه)، برای پهنای جت به تابستان قطاع آسیا (91/0 درجه) و برای تندی جت به زمستان و پاییز قطاع آفریقا (بهترتیب 92/2- و m/s 82/1) و تابستان آسیا (m/s 99/2-) مربوط میشود. پاسخ جت به گرمایش زمین، در همه شاخصها و در بیشتر فصلهای دو قطاع، در سناریوی RCP4.5 ضعیفتر از RCP8.5 است. در سناریوهای RCP4.5 و RCP8.5 میانگین جهانی واداشت تابشی کل ناشی از گسیل انسانی گازهای گلخانهای تا سال 2100 (پایان قرن بیستویکم) بهترتیب 5/4 و W/m2 5/8 فرض میشود. تغییرات برخی از شاخصها در بعضی از فصلها از دوره گذشته به آینده در میانگین همادی چند مدلی استوار است. تغییرات استوار مربوط به سناریوی RCP8.5 برای میانگین چند مدلی بدینصورت است که در فصل بهار، در هر دو قطاع آفریقا و آسیا، پهنای جت بهترتیب حدود 5/0 و 2/0 درجه و تندی جت در قطاع آسیا حدود m/s 1/1 افزایش مییابد. در تابستان، در قطاع آفریقا، تندی جت کاهش مییابد (حدود m/s 7/0) ولی در قطاع آسیا، جت علاوه بر تقویت (m/s 4/0)، حدود 8/0 درجه به سمت استوا جابهجا میشود. در فصل پاییز، در قطاع آفریقا، تندی جت کاهش (m/s 39/1-) یافته و عرض جت حدود 3/0 درجه به سمت قطب جابهجا میشود. دیگر نکته قابلتوجه آن است که برای فصل زمستان و هر دو قطاع، میانگین چند مدلی تغییر استواری را در هیچیک از شاخصهای جت، بهجز افزایش جزئی و نسبتاً استوار در پهنا برای قطاع آفریقا (حدود 2/0 درجه) نشان نمیدهد.
In this study, a set of models from the phase 5 of the Coupled Model Intercomparison Project (CMIP5) is used to examine the simulation of the upper-tropospheric subtropical African–Asian jet and its response to global warming. The ERA-Interim re-analysis dataset is used here to assess the model biases in representing the seasonal-mean jet features in the historical period (1980–2005). This study analyzes the geometrical parameters of the jet including “latitude”, “speed” and “width” in each season and for two separate sectors of the jet region: “North Africa” and “Southwest Asia”, which is briefly named “African” and “Asian” hereafter. The main features of the observed seasonal cycle of the jet in the re-analysis data is well captured in ensemble multimodel mean historical simulations: jet latitude increase (decrease) from cold (warm) to warm (cold) season and vice versa are correctly simulated for jet speed and width. In addition, in all seasons, the jet latitude and speed is greater in Asian sector than the African except for springtime jet speed. Despite the large inter-model spread in the historical jet simulations, the models do not show large systematic biases in most cases (seasons). However, systematic biases in each of the geometrical jet indices are found in some seasons: most models exhibit equatorward jet biases in summertime and wintertime of the African sector (about 1.8° and 0.9° of latitude respectively, in multimodel mean), positive biases in jet width in summertime Asia (0.9° in multimodel mean), negative biases in jet speed in summertime Asia and wintertime of the African sector (approximately 2.9 m/s) and positive jet speed biases in autumntime of the African sector (1.8 m/s). There is large spread across the models in the historical jet simulations and finding the sources of this spread and the model biases is a significant challenge that should be addressed in future works. In almost all seasons and for all of the geometrical jet indices, the multimodel mean jet response to climate change is stronger in RCP8.5 than RCP4.5 integrations. Robustness and the quantitative value of the multimodel mean jet response in each of the jet indices vary among different seasons and sectors. In winter months, we found no robust response in any of the geometrical jet indices in African or Asian sector except for a slight and relatively robust increase in jet width (0.2° of latitude in RCP8.5) in African sector. However, in other seasons, we found robust multimodel mean changes in jet indices between the historical period and the end of twenty first century (2076–2099) in the RCP8.5 scenario: In spring, models predict a robust increase in jet width of about 0.5° and 0.2° of latitude in African and Asian sectors, respectively, and also a robust increase in jet speed of 1.1 m/s for Asian sector. In summer, in the African sector, the jet speed is found to be decreased (0.7 m/s), whereas in the Asian sector, jet speed will increase (0.4 m/s), and it will move equatorward by 0.8° of latitude.
https://www.ijgeophysics.ir/article_47668_03a485f8445608ab40e929b608dd07b2.pdf
جت آفریقا-آسیا
گرمایش زمین
مدلهای CMIP5
خطای مدلها
شاخصهای هندسی
میانگین چند مدلی
African–Asian jet
Global warming
CMIP5 models
model biases
geometric jet indices
multimodel mean
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
22
32
47669
Research Article
تعیین لبههای بیهنجاریهای گرانی با استفاده از فیلتر ریختشناسی ریاضی ارتقا یافته
Edge detection of gravity anomalies by enhanced mathematical morphology filter
مجتبی بابایی
m_babaei@iauh.ac.ir
1
وحید ابراهیم زاده اردستانی
ebrahimz@ut.ac.ir
2
الهام امیریان
amirianelham@yahoo.com
3
گروه ژئوفیزیک، دانشگاه آزاد اسلامی واحد تویسرکان، تویسرکان، ایران
گروه فیزیک زمین، مؤسسه ژئوفیزیک دانشگاه تهران، تهران، ایران
گروه ژئوفیزیک، دانشگاه آزاد اسلامی واحد همدان، همدان، ایران
تشخیص لبههای چشمههای بیهنجاری، یک ابزار ضروری در تفسیر دادههای میدان پتانسیل است. روشهای زیادی برای تشخیص لبهها وجود دارد که بسیاری از آنها شامل فیلترهای بالا گذر بر اساس مشتقات دادههای میدان پتانسیل است. در این مقاله از یک روش تشخیص لبهی جدید به نام فیلتر ریختشناسی ریاضی پیشرفته (Enhanced Mathematical Morphology) یا EMM استفاده میشود. فیلتر EMM بهصورت نسبت فرسایش مشتق افقی کل به اتساع مشتق افقی کل و یا ترکیبی از این دو عملگر تعریف میشود. این فیلتر میتواند لبههای منابع کمعمق و عمیق حتی با چگالی کم را بهطور همزمان نمایش دهد. فیلتر EMM به محاسبه مشتقات عمودی نیازی ندارد که همین امر باعث میشود این روش از نظر محاسباتی پایدار باشد. در این مقاله ضمن ارائه رابطه جدیدی برای فیلتر EMM، ابتدا این فیلتر بهمنظور اعتبار سنجی روش، بر روی دادههای مصنوعی با و بدون نوفه اجرا و سپس بر روی دادههای واقعی گنبد نمکی قم آزمایش شد و مرز بیهنجاری معلوم گردید.
Enhancement on the edges of the causative source is a tool in the interpretation of potential field data. There are many methods for recognizing the edges, most of which involve high pass filters based on derivatives of potential field data. In this paper, a new edge detection method is introduce, called the enhanced mathematical morphology (EMM) filter for interpretation of field data. The EMM filter uses the ratio of the erosion of the total horizontal derivative to the dilation of the total horizontal derivative to recognize the edges of the sources, and can display the edges of the bodies simultaneously. Edge detection of the potential field data has been widely used as a significant tool for geophysical exploration technologies, which can delineate the horizontal locations of causative sources. Normally, various high-pass filters are used to recognize the edges of the potential field data (Evjen 1936; Fedi and Florio 2001; Verduzco et al. 2004; Cooper and Cowan 2006; Cooper and Cowan 2011; Ma and Li 2012; Ma 2013).
The EMM filter uses the ratio of the erosion of THD to the dilation of THD to recognize the edges of the source. Mathematical morphology was developed by Matheron and Serra in 1964 (SERRA, 1983); which is an image analysis and recognition tool. The structuring element (SE) is a basic operator in mathematical morphology, used to interact with an image and to draw conclusions about how a shape fits or misses the shapes in the image. SE consists of a matrix of 0s and 1s that can have any arbitrary shape and size. The basic operations of mathematical morphology are dilation and erosion. Dilation is defined as the maximum value in the window ascertained by the SE. Erosion is defined as the minimum value in the window ascertained by the SE. The EMM filter is expressed as (Lili et al., 2013):
where imerode (F,SE) and imdilate (F,SE) represent the erosion and dilation of the THD, respectively.
In this paper, a new relationship is presented for EMM filter that is tested on synthetic data with and without noise as well as the real potential field data in Qom salt dome. The EMM method successfully delineates the edges of the causative sources, which gives better resolution of the deeper source than other filters, and can display the edges of the bodies in a more centralized way. In this article, a new relationship is defined for the EMM filter as:
The EMM filter was used to recognize the edges of the sources. It can display the edges of the shallow and deep bodies simultaneously. The EMM filter does not require the computation of vertical derivatives, which makes this method computationally stable. The EMM filter is tested on synthetic, and real potential field data in Qom salt dome and the edge detection was done with reasonable results.
https://www.ijgeophysics.ir/article_47669_e2e07782b9f83528761ef259aec4abdf.pdf
تشخیص لبه
فیلتر ریختشناسی ریاضی ارتقا یافته
فرسایش
اتساع
گنبد نمکی قم
Edge detection
enhanced mathematical morphology
Filter
Erosion
dilation
Qom salt dom
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
33
44
47667
Research Article
کاربرد معادله اویلر دامنه سیگنال تحلیلیهای جهتی در تعیین مکان و نوع چشمههای آنومالی گرانی
Application of the Euler deconvolution of directional analytic signal amplitudes in determination of location and type of gravity anomaly sources
محمد برازش
mohammadbarazesh5300@gmail.com
1
موسسه ژئوفیزیک دانشگاه تهران، تهران، ایران
روش استاندارد اویلر روشی خودکار در تفسیر دادههای میدان پتانسیل است که در سالهای اخیر استفاده از آن گسترش زیادی پیدا کرده است. نتایج این روش وابسته به دقت پارامتر اندیس ساختاری فرض شده دارد و در حضور چشمههای تداخلی دقت نتایج کاهش مییابد. بهمنظور از بین بردن این مشکل روشهای متعددی بر اساس مشتقات این معادله طراحی شدهاند که یکی از آنها دامنه سیگنال تحلیلیهای جهتی است. دامنه سیگنال تحلیلیهای جهتی توابعی همگن هستند و در معادله اویلر صدق میکنند که بدینوسیله میتوان علاوه بر تخمین موقعیت چشمه بهصورت همزمان اندیس ساختاری را نیز تخمین زد. سیگنال تحلیلیهای جهتی از مؤلفههای تانسور گرادیان گرانی به دست میآیند و در این مقاله در حالت کلی اثبات شده است که هرکدام از توابع دامنه سیگنال تحلیلیهای جهتی در جهات x، y و z در معادله اویلر صدق میکنند. از ترکیبات آنها دو معادله جدید نتیجه شد که در تعیین موقعیت و اندیس ساختاری چشمه بسیار مؤثر است. کاربرد هرکدام از سیگنال تحلیلیهای جهتی بهصورت جداگانه و توأمان در تعیین موقعیت و اندیس ساختاری به مدلی در حضور چشمههای تداخلی و نوفه گاوسی نسبتاً بالا اعمال شد. دامنه سیگنال تحلیلی در جهت z و استفاده همزمان از هر سه دامنه سیگنال تحلیلیهای جهتی مشابه با دو معادله جدید معرفی شده در این مقاله نسبت به سایر معادلات نتایج بهتری ارائه دادند. سرانجام این روشها بر رویدادههای واقعی گرانی کانسار منگنزی صفو واقع در 25 کیلومتری شهرستان چالدران بکار رفت و برای این معدن مقداری برای عمق (حدوداً 6 متر) و اندیس ساختاری با مقداری منفی و نزدیک به صفر تخمین زده شد.
The components of Gravity Gradient Tensor (GGT) is used for second-order derivatives of the gravitational potential field in the directions x, y and z in a Cartesian coordinate system. The third column of the gravity gradient tensor is Hilbert transform pairs of the first and the second columns. Many methods have been designed to estimate the depth, the horizontal position and the type of the sources from gravity gradient tensor components. Often, these methods are used derivatives of potential field data or their compounds in directions x, y, and z. Standard Euler deconvolution method is an approach in the interpretation of potential field data. It is able to locate the sources and to estimate the regional parameters with the assumption of the structural index. This approach is an automated method that has seen rapid development in recent years. The result of this method closely related to the precision of the assumed structural index parameter, and the accuracy is reduced in the presence of interference sources. Euler deconvolution of the directional analytic signal amplitudes is one of many methods to eliminate this problem. It is shown that the components of the gravity vector satisfy Euler's equation. Thus, it is proved that the amplitudes of directional analytic signal are homogenous and by putting in Euler's equation can estimate the location and the structural index of the gravity anomalies. In addition, two new equations were obtained from the combination of directional analytic signal amplitudes that is very effective in locating and estimating the structural index of gravity sources.
This paper was examined the application of Euler's equation of the directional analytic signal amplitudes to determine the location and the structural index of gravity anomaly sources. First, it is proved that each of directional analytic signal amplitudes in directions x, y, and z satisfy Euler's equation. Second, using the combination of directional analytic signal amplitudes derived two new equations that is more successful in determining the location (horizontal positions and depth) and source type (structural index) directly over the edges of gravity anomaly sources. The maxima of analytic signal amplitude in the z- direction place directly on the edge of the anomaly sources, but the maxima of analytic signal amplitude in the x- and y- directions deviate from the edges. That is why the simultaneous use of two or three-directional analytic signal amplitude can provide more accurate solutions.
The method described above was tested on the synthetic model in the presence of relatively high level Gaussian noise and interference sources. Finally, the method was applied to the Safoo manganese ore and obtained horizontal position, depth (~6 m) and structural index. MATLAB software was used to apply the above-mentioned methods.
https://www.ijgeophysics.ir/article_47667_d683b36142828b7131dbe113dd8037d3.pdf
اویلر
واهمامیخیت
سیگنال تحلیلیهای جهتی
تانسور گرادیان گرانی
Euler deconvolution
directional analytic signal amplitudes
gravity gradient tensor
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
45
56
53424
Research Article
تعیین محدوده مخزن با استفاده از وارونسازی همزمان دادههای لرزهای
Determining the area of reservoir using simultaneous inversion of seismic data
حامد قنبر نژاد مغانلو
hqanbarnejad@ut.ac.ir
1
محمد علی ریاحی
mariahi@ut.ac.ir
2
مجید باقری
majidbagheri@ut.ac.ir
3
سید محسن سیدعلی
mohsen.seyedali@gmail.com
4
موسسه ژئوفیزیک دانشگاه تهران، تهران، ایران
مؤسسه ژئوفیزیک دانشگاه تهران، ایران
موسسه ژئوفیزیک دانشگاه تهران، تهران، ایران
شرکت نفت فلات قاره ایران، تهران، ایران
لرزهشناسی اکتشافی بهطور گستردهای مورد توجه مهندسین نفت بوده است، بهطوریکه بیشترین میزان هیدروکربنهای شناخته شده در سالهای اخیر با استفاده از روشهای لرزهای بوده است. تعیین محدوده مخزن با استفاده از روش لرزهشناسی بازتابی بهطور جدی از سال 1930 شروع شده و تاکنون مورد توجه بوده است. قدرت تفکیک بالای دادهها و هزینه پایین در مقایسه با حفاری، این روش غیرمستقیم را در اکتشاف هیدروکربنها بسیار ضروری نموده است. تاکنون روشهای مختلفی برای تعیین محدوده مخازن هیدروکربنی معرفی شده است که یکی از معروفترین آنها روشهای وارونسازی دادههای لرزهای است.
در وارونسازی پس از برانبارش، از دادههای لرزهای پس از برانبارش استفاده میشود. در این حالت تنها اطلاعاتی که میتواند از روی دادهها تخمین زده شود، مقاومت صوتی موج تراکمی میباشد. در وارونسازی همزمان پیش از برانبارش که از دادههای لرزهای پیش از برانبارش استفاده میکند، علاوه بر مقاومت صوتی موج تراکمی، اطلاعات موج برشی، چگالی و همچنین نسبت پواسون از روی دادهها قابل استخراج میباشد. بنابراین وارونسازی همزمان پیش از برانبارش اطلاعات بیشتری را نسبت به وارونسازی پس از برانبارش نتیجه میدهد. در این مقاله از روش وارونسازی همزمان دادههای لرزهای پیش از برانبارش برای تعیین محدوده مخزنی در یکی از میادین هیدروکربنی خلیجفارس استفاده شده است. با استفاده از این روش، کمیتهای مقاومت صوتی موج تراکمی، مقاومت صوتی موج برشی و چگالی تخمین زده شدند. در محدوده مخزنی کمیتهای مورد نظر با کاهش غیرعادی مواجه شدند که دلیلی بر وجود ناهنجاری هیدروکربنی در این محدوده است و بهاینترتیب محدوده مخزن مورد نظر شناسایی شد.
Poststack seismic inversion generally transmutes seismic amplitude to P-wave acoustic impedance, which lacks low-frequency component due to the stacking process. This component should be compensated using well logs as a priori constraint. If this low-frequency trend is known with adequate accuracy, poststack inversion could produce precise results. Nevertheless, in most cases, the mentioned information are far from the true model. In such cases, poststack inversion results could have high uncertainty. Because there is no mode conversion at normal incidence, postsatck inversion is completely acoustic, hence P-wave impedance is the only information which can be extracted from poststack inversion of P-wave data.
In simulations prestack inversion, in addition to the P-wave acoustic impedance, S-wave information, density, and Poisson’s ratio can also be derived from prestack data. Thus, prestack inversion can be used to get more information than poststack inversion. The two-step process of acoustic impedance and shear impedance by model-based inversion is replaced by one-step pre-stack simultaneous inversion.
In order to apply simultaneous inversion method to our prestack seismic data, the data should be transformed from offset domain to angle domain as the first step. A useful approach is to calculate offset as a function of incidence angle, using Snell’s Law to follow the ray path through the layers if velocity information is available. The next step is to build initial models of acoustic impedance, shear impedance, and density. We built these initial models using sonic log, Delta-Time Shear (DTSM) log and RHOB log which were available in the interest area.
There are two relationships that should hold for these wet rocks. The first relationship uses this fact that in wet clastics the ratio of the s-wave velocity over p-wave velocity should be constant within a rock layer. After reformulation of the mentioned trend, one can understand that the natural logarithm of shear impedance has a linear relationship with the natural logarithm of acoustic impedance. The second fact uses Gardner equation. After reformulation of the Gardner relation, it is understandable that the natural logarithm of density has a linear relationship with the natural logarithm of acoustic impedance, too.
We determined k , kc , m and mc which respectively are slope of the natural logarithm of shear impedance against natural logarithm of acoustic impedance, intercept of the natural logarithm of shear impedance against natural logarithm of acoustic impedance, slope of the natural logarithm of density against natural logarithm of acoustic impedance, intercept of the natural logarithm of shear impedance against natural logarithm of acoustic impedance.
Besides, we need a set of angle-dependent wavelets which are derived from angle stacks. Hence, we built three angle stacks; near-angle stack (0 to 11 degrees), middle-angle stack (11 to 20 degrees) and far-angle stack (20 to 29 degrees). Using these angle stacks, we built three statistical angle-dependent wavelets from three angle stacks. Finally, with log information, we built an initial model for acoustic impedance and tried to solve the inversion matrix using conjugate gradient method. Solving the equation, we can derive acoustic impedance, shear impedance, and density sections simultaneously from prestack data. Using simultaneous inversion, we identified hydrocarbon reservoir.
https://www.ijgeophysics.ir/article_53424_7c64171a9e726d4426bb26f5c1458080.pdf
وارونسازی همزمان پیش از برانبارش
وارونسازی پس از برانبارش
مقاومت صوتی موج تراکمی
محدوده مخزن
مقاومت صوتی موج برشی
چگالی
prestack simultaneous inversion
poststack inversion
P-wave acoustic impedance
reservoir area
S-wave acoustic impedance
Density
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
57
74
54700
Research Article
کاربست روش مککورمک فشرده مرتبه چهارم با پیمایش زمانی رونگه-کوتا چهار مرحلهای برای مسئله تنظیم راسبی غیرخطی ناپایا
Application of the fourth-order compact MacCormack scheme with a four-stage Runge–Kutta time marching for numerical solution of unsteady and non-linear Rossby adjustment problem
رضا جوان نژاد
rjavannezhad@gmail.com
1
سرمد قادر
sghader@ut.ac.ir
2
واحد علوم و تحقیقات، دانشگاه آزاد اسلامی، تهران، ایران
گروه فیزیک فضا، موسسه ژئوفیزیک دانشگاه تهران، تهران ایران
در پژوهش حاضر، به بررسی خطای بریدگی و آهنگ همگرایی روش مککورمک فشرده مرتبه چهارم با پیمایش زمانی رونگه-کوتا چهار مرحلهای پرداخته میشود. برای انجام این تحلیل از معادله فرارفت خطی یکبعدی استفاده شده است که دارای حل تحلیلی است. خطای بریدگی برای روشهای مککورمک مرتبه دوم، مککورمک فشرده مرتبه چهارم با پیمایش زمانی مرتبه دوم و مککورمک فشرده مرتبه چهارم با پیمایش زمانی رونگه-کوتا چهار مرحلهای استخراج و بیان شده است. برای بهدست آوردن خطای بریدگی از معادله فرارفت خطی یکبعدی با ضریب ثابت استفاده شده است. همچنین برای بررسی دقت آهنگ همگرایی برای روشهای متفاوتی از جمله لکسوندرف، لیپفراگ، بیموارمینگ و همچنین روشهای مککورمک کلاسیک با گسستهسازیهای مکانی مرتبهدوم و فشرده مرتبهچهارم و گسستهسازی زمانی رونگه-کوتا چهار مرحلهای بهدست آمده است. نتایج نشان میدهد آهنگ همگرایی بهدستآمده برای مسئله خطی متناسب با آهنگ همگرایی نظری میباشد. در ادامه نتایج حل عددی برای مسئله تنظیم راسبی غیرخطی ناپایا و برای حالتهای یکبعدی و دو بعدی با استفاده از روش مککورمک فشرده مرتبه چهارم با پیمایش زمانی رونگه-کوتا چهار مرحلهای مورد بررسی قرار گرفته است. مقایسه نتایج بهدستآمده از روش مذکور برای دو مسئله یک و دو بعدی با نتایج حاصل از کارهای سایر محققان گویای عملکرد مناسب این روش است بهویژه هنگامیکه میدان حل با ناپیوستگی همراه باشد.
By increasing the computing power of computers, the advantage of high-resolution numerical methods for numerical simulation of the governing equations of fluid flow is further emphasized. Recently, increasing the accuracy of numerical methods used for simulation of fluid dynamics problems, particularly the geophysical fluid dynamics problems (e.g., shallow water equations) has been the subject of many research works.
The compact finite difference schemes can provide a simple way to reach the main objectives in the development of numerical algorithms, i.e., having a low cost on the one hand and a highly accurate computational method on the other hand. These methods have also been used for numerical simulation of some geophysical fluid dynamics problems.
However, by splitting the derivative operator of a l compact centra method into one-sided forward and backward operators, a family of compact MacCormack-type schemes can be derived (Hixon and Turkel, 2000). While these classes of compact methods are as accurate as the original compact central methods used to derive the one-sided forward and backward operators, they need less computational work per grid point.
The present work is devoted to the assessment of the accuracy of different methods. The one-dimensional advection equation with the known analytical solution is employed as a prototype model. Also, the truncation error of the traditional second-order MacCormack scheme, the standard fourth-order compact Mac-Cormack scheme, and a fourth-order compact MacCormack scheme with a four-stage Runge–Kutta time marching method are studied. Furthermore, to be able to examine the accuracy, the Lax–Wendroff, the leap-frog and the Beam–Warming methods combined with the second-order and fourth-order compact finite difference methods for spatial differencing are also used. In addition, the convergence rates of different methods are studied. It can be seen that the convergence rates are in agreement with the theoretical order of convergence.
In this work, the traditional second-order MacCormack scheme (MC2), the standard fourth-order compact Mac-Cormack scheme (MC4) developed by Hixon and Turkel (2000) and a fourth-order compact MacCormack scheme with a four-stage Runge–Kutta time marching method (MCRK4) are used for numerical solution of the unsteady and non-linear Rossby adjustment problem (one- and two-dimensional cases). In the one-dimensional case, a single layer shallow water model is used to study the unsteady and nonlinear Rossby adjustment problem. The conservative form of the two-dimensional shallow water equations is used to study the unsteady and nonlinear Rossby adjustment problem in the two-dimensional case. For both cases, the time evolution of a fluid layer initially at rest with a discontinuity in height filed is considered for numerical simulations.
https://www.ijgeophysics.ir/article_54700_b6e316eb4967588d185c3c6aaa207527.pdf
روش مککورمک فشرده
دقت عددی
خطای بریدگی
تنظیم راسبی غیرخطی
compact MacCormack scheme
numerical accuracy
Truncation error
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
75
89
54701
Research Article
توفانهای خاک در منطقه غرب و جنوب غرب ایران و تأثیر آنها بر شارهای تابشی: مطالعه موردی
Dust storms in west and southwest Iran and their impact on radiation fluxes:
A case study
سعید فرهادی پور
saeed.farhadypoor@gmail.com
1
مجید آزادی
azadi68@hotmail.com
2
عباسعلی علی اکبری بیدختی
bidokhti@ut.ac.ir
3
امید علیزاده چوبری
omid.alizade@gmail.com
4
حبیب الله سیاری
farhadypour.s@asmerc.ir
5
پژوهشکده هواشناسی، تهران، ایران
پژوهشکده هواشناسی، تهران، ایران
موسسه ژئوفیزیک دانشگاه تهران، تهران، ایران
موسسه ژئوفیزیک دانشگاه تهران، تهران، ایران
دانشگاه عالی دفاع ملی، تهران، ایران
توفانهای خاک که امروزه در بسیاری از مناطق ایران (بهخصوص در غرب و جنوب غرب ایران) فراوانی وقوع بالایی دارد، علاوه بر تأثیر نامطلوب بر سلامتی بشر از طریق آلودگی هوا، تأثیر قابل ملاحظهای بر خواص نوری و توازن تابشی منطقه میگذارد. در مطالعه حاضر اثرات تابشی ناشی از توفان خاک در بازه زمانی 16 تا 21 ژوئن 2012 در منطقه غرب و جنوب غرب ایران با استفاده از مدل عددی WRF-Chem بررسی شده است. ابتدا عملکرد مدل با استفاده از دادههای اندازهگیری ایستگاهی (ایستگاههای اندازهگیری آلودگی وابسته به سازمان محیطزیست و ایستگاههای AERONET) و دادههای ماهوارهای MODIS، OMI و CALIPSO مورد ارزیابی قرار گرفت. نتایج ارزیابی مدل نشان از برآورد بیشتر غلظت PM10 در ایستگاه اهواز و در اغلب موارد برآورد کمتر مقادیر عمق نوری هواویزها (aerosol optical depth) در ایستگاههای AERONET دارد. با این حال، عملکرد مدل در شبیهسازی روند تغییرات و میزان گرد و خاک در طی توفان مذکور قابل قبول است، بهطوریکه توزیع افقی و قائم گرد و خاک شبیهسازیشده توسط مدل و مشاهده شده توسط ماهواره الگوهای تقریباً مشابهی را نشان میدهند. ذرات گرد و خاک در سقف جوّ و سطح زمین دارای اثرهای سرمایشی، اما در میانه جوّ دارای اثر گرمایشی هستند. میانگین پریشیدگی تابش طول موج کوتاه توسط گرد و خاک در منطقه غرب و جنوب غرب ایران در بازه زمانی 17 تا 20 ژوئن 2012 در سطح زمین، میانه جوّ و سقف جوّ بهترتیب 27/7-، 79/1 و W m-2 47/5 برآورد شد.
Dust aerosols make a considerable contribution to the climate system through their radiative effects due to their abundance in the atmosphere. Recent observations suggest that over the past decade, dust events have become more frequent in many parts of Iran, especially in the west and southwest. Through their radiative forcing, dust aerosols have significant effects on the regional radiation budget of the atmosphere, while their adverse effects on human health have also raised serious concerns. The primary aim of the present study is to examine the radiation effects associated with a severe dust storm that occurred in west and southwest Iran on 16 to 21 June 2012. To this end, the Weather Research and Forecasting with Chemistry (WRF-Chem) model was used. Two simulations were conducted: a model setup that did not include dust aerosols, and the one that included dust aerosols and their feedback to the atmosphere. A two-way interactive nested domain (nesting ratio:1:3) simulations were performed using 98 Í 90 and 151 Í 139 horizontal grid points, respectively. In the vertical, 27 σ-levels were used. The grid spacing for the two domains were 45 and 15 km, respectively. Simulations ran from 16 to 22 June 2012, and the first 24 hours was considered as the spin-up time. Meteorological initial conditions were obtained from the Global Forecast System (GFS) data at 0.5˚Í 0.5˚ resolution. The performance of the model was evaluated using the available observed data, including PM10 observations in Ahwaz located in southwest Iran, available AErosol RObotic NETwork (AERONET) data in nearby areas, and aerosol products of the Moderate Resolution Imaging Spectroradiometer (MODIS), the Ozone Monitoring Instrument (OMI) and the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) carried on board the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) spacecraft. Results indicate that PM10 concentration in Ahwaz is overestimated by the model, while simulated aerosol optical depth (AOD) is underestimated compared to the observed AERONET data. Relatively, good agreement is found between the model results and satellite products, and temporal evolution of the dust events is also well-simulated. Thus, generally, the performance of the model is acceptable for accurate simulation of the dust event. Our analysis indicated that radiative effects of dust particles cause cooling at the surface and top of the atmosphere, but warming in the middle of the troposphere. On average, perturbation of shortwave radiation by dust aerosols in the west, and southwest Iran is estimated to be -7.27, 1.79 and -5.47 W m-2 at the surface, in the middle and at the top of the atmosphere, respectively. Average perturbation of the longwave radiation by dust aerosols over the same region was estimated to be 2.2, -1.61 and 0.59 W m-2 at the surface, in the middle and at the top of the atmosphere, respectively. Thus, the net (shortwave + longwave) radiative effect of dust aerosols averaged in west and southwest Iran is found to be -5.07, 0.19 and -4.88 W m-2 at the surface, in the middle and at the top of the atmosphere, respectively.
https://www.ijgeophysics.ir/article_54701_e17fbb43dc8dfcf1a8a1c29ec0b83783.pdf
توفان گرد و خاک
مدل WRF-Chem
دادههای ماهواره
پریشیدگی تابش
Dust Storm
WRF-Chem
Satellite data
perturbation of radiation
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
90
104
49848
Research Article
تفسیر بیهنجاریهای میدان پتانسیل بهمنظور بررسی ساختارهای زمینشناسی و اکتشاف تلههای نفتی احتمالی، بررسی موردی: جنوب سمنان
Interpretation of potential field anomalies to investigate geological structures and oilfield exploration: A local study in the south Semnan
یاسین ناسوتی
yasin.nasuti@gmail.com
1
براتعلی فیض آبادی
afeizabady@yahoo.com
2
سید مجتبی تاجور
mojtaba.tajvar@yahoo.com
3
عزیز ناسوتی
aziz.nasuti@ngu.no
4
گروه فیزیک، دانشکده علوم پایه، دانشگاه حکیم سبزواری، سبزوار، ایران
گروه فیزیک، دانشکده علوم پایه، دانشگاه حکیم سبزواری، سبزوار، ایران
دانشکده علوم پایه، دانشگاه آزاد اسلامی واحد علوم و تحقیقات تهران، ایران
محقق زمینشناسی و ژئوفیزیک، گروه ژئوفیزیک، سازمان زمینشناسی نروژ، ترندهیم، نروژ
برای تعیین موقعیت منابع مولد بیهنجاریهای میدان پتانسیل نیاز به پردازش و تفسیر دقیق نقشههای بیهنجاری حاصل از برداشت دادههای میدان پتانسیل است. تفسیر بیهنجاریهای میدان پتانسیل با استفاده از روشهای مختلفی انجام میگیرد. در این مقاله جهت تفسیر و جداسازی بیهنجاریهای میدان پتانسیل از فیلترهای ژئوفیزیکی مانند روش روند سطحی، زاویه تیلت و ادامه فراسو استفاده گردید. در نهایت نقشههای هر یک از این روشها تهیه و نتایج با هم مقایسه شد. نتایج حاصل از این نقشهها نشان میدهد که هشت ساختمان تاقدیسی، پنج ساختمان ناودیسی، تعداد زیادی گسل و گنبد نمکی در ناحیه مورد بررسی وجود دارد. از بین تاقدیسهای تفسیر شده، تاقدیسB از عمق و ابعاد قابل توجهتری برخوردار است، بهطوریکه با ادامه فراسوی 3000 متر اثر این تاقدیس باقی مانده است. همچنین هیچگونه توده آذرینی بر روی نقشه مغناطیسی آن مشاهده نشده است. از آنجا که هیدروکربورها در سنگهای رسوبی یافت میشوند و در سنگهای آذرین وجود ندارند، بنابراین میتوان این تاقدیس را از نظر تجمع منابع هیدروکربوری، تله مناسبی تلقی کرد. این نقشهها یکی شدن گسل جنوب تاقدیس کوه سرخ با گسل جنوب تاقدیس سیاه کوه را تأیید میکنند. همچنین این نقشهها نشان میدهند که سیاه کوه ادامه کوه سرخ است که بهوسیله فعالیت گسلهای منطقه، فرسایش و نهشته شدن رسوبات جدید از هم جدا شدهاند.
South of Semnan province is an important region in Iran for the formation of oil traps, which has been received much attention for several decades. Many sulfur mines have also been discovered in this area. The sulfur is most probably derived from natural gases that are guided by regional faults. Numerous anticlines and salt domes are also present in this area, playing an important role in the formation of oil traps. As a consequence, this region seems to have a great potential to form hydrocarbon traps. Since the area is covered by Tertiary, Quaternary to current sediments, gravity, and magnetic studies are very useful to investigate and explore the geological structures in this region. Gravity and magnetic studies, in a trapezium grid, were performed to investigate the underground structures, sediment thicknesses, plutonic and volcanic igneous rocks and also hydrocarbon traps. Data acquisition was carried out on 86 profiles with a northwest-southeast trend. The distance between profiles and acquisition stations were selected to be 1000 meter. Due to the presence of swamps, mountainous areas, river, terrain, etc., some stations were removed from the survey plan. In addition, some transects measured in a direction deviating from the straight line. As a result, in some of the profiles, a number of missing stations can be observed. In some of the transects, profiling was not conducted in a regular 1000 m line spacing too. The average magnetic inclination and declination are 53.8 and 4.2 degrees, respectively, and the average total magnetic field is estimated 48181 nT as well. The gravity and magnetic data were collected simultaneously. In this study, to interpret and discern potential field anomalies, we applied geophysical filters such as surface trend, the tilt angle, and upward continuation. After producing a geophysical map for each method, the results were jointly interpreted. Joint interpretation demonstrated that eight anticlines, five synclines, several faults and salt domes could be recognized. Among these anticlines, anticline B has a considerable depth and dimension so that with 3000 m upward continuation, the gravity field can still be seen. Meanwhile, no igneous rocks were observed on the magnetic maps. Therefore, this anticline can be considered as an appropriate trap for the accumulation of hydrocarbon. These maps confirm the joining fault to the south of Kohe- Sorkh anticline with the fault to the south of Siah-Koh anticline, which are located in the south of Abulabad. The results show that the potential methods appear to be promising to characterize subsurface structures for the initial phase of hydrocarbon explorations. These maps also show that the Siah-Kuh is a continuation of Kuh-Sorkh, which are separated by regional fault activities, erosion, and new sedimentation. Consequently, it is suggested to use numerical modeling to define the shape, dimension and depth of anomalies, especially for the interpreted anticline B. Finally, a seismic survey can be performed over the potential anomalies that might have hydrocarbon accumulation.
https://www.ijgeophysics.ir/article_49848_f7cb0b4f2bb769d25873a61f7ea2e413.pdf
روند سطحی
زاویه تیلت
ادامه فراسو
تفسیر زمینشناسی
تلههای نفتی
ایران مرکزی
surface trend
tilt angle
upward continuation
geological interpretation
oil traps
Central Iran
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
105
118
54702
Research Article
بررسی پارامترهای لرزهای سیستمهای فولادی برون مرکز با تیر پیوند متوسط
A comparative study of seismic performance evaluation existing methods EBF systems with intermediate link beam
علی اکبر میمندی پاریزی
:milad.mohammadian@mail.sbu.ac.ir
1
میلاد محمدیان
milad.engeener@gmail.com
2
عباس مهدویان
ahmad.nik1395@gmail.com
3
مهندسی عمران گرایش زلزله، دانشگاه شهید بهشتی، تهران، ایران
مهندسی عمران گرایش زلزله، دانشگاه شهید بهشتی، تهران، ایران
دانشکده مهندسی عمران، آب و محیط زیست، دانشگاه شهید بهشتی، تهران، ایران
فلات ایران از نظر وقوع زلزله یکی از فعالترین مناطق دنیا بوده و هر از چند گاهی زمینلرزههای مخرب و مصیبتبار با آسیبهای جانبی و مالی وسیع در کشورمان به وقوع میپیوندد که پیشگیری از وقوع یا کاهش تلفات جانی و مالی ناشی از این زلزلهها در کشورمان ایران از اهمیت خاصی برخوردار است. خرابیهای وارده به سازه در زلزلههای گذشته که مطابق با آییننامهها طراحی گردیدند نقاط ضعف شیوههای موجود در طراحی لرزهای سازهها را بهخوبی مشخص ساخته است. لذا در حال حاضر بیشتر آییننامههای طراحی لرزهای بهسرعت به سمت طراحی بر اساس عملکرد حرکت میکنند که ساختمان با توجه به عملکردی که در زلزله از خود نشان میدهد طراحی گردد. با این اوصاف در سالهای اخیر بحث رفتار غیرخطی و روش عملکردی در ارزیابی لرزهای سازهها مطرح گردیده است. اگرچه دقیقترین شیوه در ارزیابی رفتار لرزهای سازهها استفاده از تحلیل دینامیکی غیرخطی تاریخچه زمانی است، اما دستورالعملهای معتبر (دستورالعملهایی همچون ATC و یا FEMA، آییننامههای بهسازی لرزهای ساختمان)، استفاده از روش تحلیل استاتیکی غیرخطی فزاینده بار افزون را به دلیل سادگی در بهکارگیری و تفسیر نتایج توصیه کردهاند. در پژوهش پیش رو ضوابط دستورالعملهای مذکور در خصوص ارزیابی لرزهای سیستمهای برونمحور (EBF) با الگوهای مختلف بار جانبی در روش بار افزون سنتی و بار افزون پیشرفته در مقایسه با تحلیل دینامیکی غیرخطی مورد ارزیابی قرارگرفته است؛ و نتایج بهصورت پیشنهادهایی ارائه شدهاند؛ اما در کل نتایج حاکی از آن است که اختلاف میان پاسخهای حاصل از روشهای مختلف بار افزون، در سازههای کوتاهمرتبه ناچیز است.
Natural disasters and their harmful impacts have always been one of the most challenging problems all over the world. As such phenomena are inevitable since the distant past mankind has been trying to predict them spatially and temporally and evaluate their loss. Earthquake, as a natural disaster, could not be predicted in time. However, its magnitude and location are predictable to some extent and so is the corresponding loss. Theoretical and computational advances in civil engineering lead to a precise understanding of structures’ behavior and earthquakes. Therefore, in the recent decades, nonlinear behavior and performance-based method have been introduced in the seismic evaluation of structures. Many studies have been carried out in this field by research centers and agencies like FEMA and ATC, resulted in useful guidelines. Eccentrically braced frames have high stiffness and suitable energy damping against the lateral forces like the earthquake. In this bracing system, the required stiffness and formability of the frame is provided by the link beam, and are dependent on the details and characteristics of the link beam. In recent years, Eccentrically Braced Frames (EBF) has been utilized as a resistant system against the earthquake lateral forces. The research has shown that the EBF have the ability to combine a high stiffness in the elastic range as well as an excellent ductility and energy dissipation in the inelastic range. Currently, seismic design provisions of most building codes are based on strength or force (base shear) considerations. These building codes are generally regarding the seismic effects as equivalent static forces with a height wise distribution, which is consistent with the first vibration mode shape. However, the design basis is being shifted from strength to deformation in modern performance-based design codes. Determining the shear story and overturning moment under earthquake excitation is an important problem in the seismic design of structures. There are several approaches in order to estimate an acceptable accurate response for the shear story and overturning moment of the structure in the nonlinear region. Both ATC and FEMA approaches are good ideas to evaluate the seismic performance, but more simplified approaches should be applied in seismic design codes. Most of the seismic design codes suggest a very simple relation for estimating the shear story in design base earthquake. In this study, some criterions of the mentioned guidelines are studied, which are about seismic evaluation of the eccentric braced frame (EBF) systems, then the suggestions are offered. In this research, a comparative study has been done to analyze the behavior of regular steel building structures of 4, 8, 12 and 16 stories, located in zones of high seismic hazard and soil type 2. Three-dimensional building systems composed of steel frame system with Intermediate Link Beam (EBF) have been selected for investigation. These 3D building structures have been considered with 4, 8, 12 and 16 stories. Then, the performance level of all regular structures is evaluated in one hazard level (with the return period of 475 years). In order to evaluate the performance level of the aforementioned structures, they were modeled three dimensionally using SAP V14.00 software for both nonlinear static and dynamic analysis. The criteria for predicting the target location guidelines ATC-40, FEMA-356 and modified methods in FEMA-440 were used. The loading pattern design for nonlinear static analysis of single-mode and modal pushover (MPA) was used. For nonlinear time history dynamic analysis out of nine coupled ground motion accelerations from the strong motion database of PEER, with a minimum of 20 km and maximum 45 km from the source and magnitude range of 6 to 7.5 were selected. The performed procedures in FEMA-356 and proposed plastic hinges in this guideline are utilized for performing the static nonlinear analysis. The soil type II was considered having the shear wave limit between 375 to 750 m/sec. The result and the accuracy of pushover analysis has been compared with the nonlinear time history analysis. This indicates that the results obtained by FEMA-440, are closer to the results of the nonlinear time history dynamic analysis. It is also concluded by the investigating of the shear story and overturning moment of the mentioned models that these parameters are dependent greatly on the length of the link beam and inadequacy of push-over analysis in demonstrating tall buildings performance are other results of this study.
https://www.ijgeophysics.ir/article_54702_c943db1b5e8e895f338f8e0aad91f717.pdf
تحلیل پوش آور
تحلیل دینامیکی غیرخطی
الگوی بار جانبی
تیر پیوند
دستورالعمل فما
lateral load patterns
link beam
Pushover analysis
FEMA
Nonlinear dynamic analysis
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
119
127
54708
Research Article
تلفیق دادههای مقاومت ویژه و مغناطیس جهت تعیین شرایط کانهزایی و گسترش کانسار آهن موئیل مشکینشهر
Combination of Magnetic and Resistivity methods for determination of the mineralization extension in Moeil Iron ore, Meshkin-Shahr
امیر امیرپور اصل میاندوآب
am_amirpour@sut.ac.ir
1
قهرمان سهرابی
q_sohrabi@yahoo.com
2
نوید شاد منامن
shadmanaman@sut.ac.ir
3
هاشم طباطبایی رئیسی
dr.tabatabaee@yahoo.com
4
دانشکده اکتشاف معدن، دانشگاه صنعتی سهند، تبریز، ایران
دانشکده علوم دانشگاه محقق اردبیلی، اردبیل، ایران
دانشکده معدن دانشگاه صنعتی سهند، تبریز، ایران
واحد دماوند، دانشگاه آزاد اسلامی، تهران، ایران
کانسار آهن گوتیتی مرتبط با سیستم چشمههای آبگرم منطقه موئیل در یک کیلومتری جنوب روستای موئیل، واقع در 17 کیلومتری جنوب خاوری مشکینشهر از توابع استان اردبیل، و در دامنه غربی کالدرای سبلان در شمال باختری ایران قرار دارد. در این تحقیق با استفاده از مطالعات ژئوفیزیکی انجام گرفته در منطقه مذکور شامل روش مغناطیسسنجی در مرحله نخست، به حجم 750 نقطه بهعنوان روش اکتشافی غیر مستقیم و روش مقاومت ویژه الکتریکی در مرحله بعد، بهعنوان روش مستقیم (همزمان با استفاده از سه آرایش گرادیان، شلومبرگر و ونر) به حجم هفت پروفیل به طولهای 150 تا 350 متر و نیز حفاری اکتشافی به روش گمانهزنی پودری به تعداد 16 گمانه به عمق متوسط 15 متر، نحوه گسترش قائم و افقی و شرایط کانهزایی آهن مورد بررسی قرار گرفته است. بر اساس اطلاعات روش مغناطیسسنجی، اولاً احتمال وجود کانهزایی مگنتیتی در زیر رخنمون گوتیتی بررسی و منتفی شد. همچنین با مقایسه اطلاعات مذکور با شواهد زمینشناسی موجود، مشخص شد نواحی با آنومالی مثبت بالا مرتبط با تودههای گرانیتی، آنومالی حد واسط مرتبط با کانهزایی آهن و نواحی با آنومالی منفی مرتبط با چشمههای مولد کانهزایی میباشند. وجود یک گسل اصلی با امتداد شمال غرب – جنوب شرقی و یک گسل متقاطع فرعی عمود بر گسل اصلی در نقشه مغناطیسی قابل تشخیص است. همچنین شواهد گسل اصلی و نیز یک لایه آبدار در نیمرخهای مقاومت ویژه الکتریکی آشکار است. لایه آهندار با مقاومت الکتریکی متوسط بین لایه آبدار و لایه آذرآواری با مقاومت بالا قرار گرفته است. نتایج بهدستآمده از اندازهگیریهای مغناطیسی و مقاومت ویژه الکتریکی تطابق خوبی با اطلاعات زمینشناسی و حفاریهای اکتشافی انجام گرفته در منطقه مورد مطالعه دارد.
Moeil Goethite iron mineralization zone related to Moeil hot streams system is located in 1 Km south of the Moeil village, at 17 km southeast of Meshkin Shahr town in Ardebil province. This zone is at the western slopes of Sabalan caldera in the northwest part of Iran. The studied area is located in Central Iran zone according to Geological Structures’ classification by Stoklin (1988) and in Alborz-Azerbayjan zone according to Nabavi (1355). General lithology of this area consists of Cenozoic and Quaternary volcanic rocks and pyroclastic material related to the Sabalan volcanism. The reason for the high thermal gradient of the region is because of hot intrusive bodies in depth. In this area, geothermal liquids move upward via fractures and fault systems when the atmospheric water is penetrated and contacted with deep intrusive bodies then hot water springs (e.g., Geinarja and Moeil) have been generated consequently. Chlorinated hot water leaches iron of Ferro-magnesium minerals, through moving beside the mafic and intermediate rocks, and deposits iron hydroxides (Goethite and Limonite) on the surface. Moeil iron ore in the south of Meshkin Shahr is considered as a prominent iron ore related to the hot water springs in Iran. In this ore, the average amount of iron oxide (Fe2O3) is 70 percent. In this paper, extension of Iron ore body, and mineralization situation was studied in depth and horizon by implementation of two geophysical approaches, including magnetic method as an indirect way for exploration of Hematite and as a method for checking possibility of magnetite mineralization in depth (about 750 data points) and electrical resistivity method as a direct measurement (seven data profiles all with lengths between 150 to 350 meters) by utilizing three different arrays at the same time (Gradient, Schlumberger and Wenner arrays) and exploration boreholes (16 boreholes with depth of 15 meters). As a result of the magnetic survey, firstly, the possibility of magnetite mineralization has been rejected. Secondly, high positive magnetic anomalies are related to granitic rocks, intermediate positive anomalies are related to Iron ore body, and low anomalies are related to springs that the minerals become precipitated. Besides, the existence of two conjugate faults with NW-SE and NE-SW directions is clear in magnetic anomaly map. In electrical resistivity pseudo sections, the presence of a nearly vertical fault is obvious. In addition, the existence of an aquifer with a very low resistivity at the bottom is detectable. Iron mineralization as a nearly horizontal layer is located between this aquifer in deep and high resistivity volcanic sediments at the top. Finally, as an outcome of this study, it can be mentioned that results of both magnetic and electrical resistivity activities, geology evidences and exploration drilling results confirm each other. In general, the below results have been achieved through Magnetic and resistivity approaches:
- Magnetite mineralization is not expected in this ore due to the relatively low intensity of recorded magnetic signals.
- From the magnetism viewpoint, iron layers do not have a sharp difference with base rocks (tuff).
- In some points, non-ferrous volcanic rocks (granitic) show high magnetism in comparison with mineralized zones. High intensity (more than 48550 nT) is related to these granitic bodies.
- The intensity between 48550 nT and 48800 nT is related to the iron mineralized zone.
- Results from geoelectric studies show that the specific resistivity between 150 and 700 ohm-meter is related to the iron layers, and the figures between 700 and 800 ohm-meter is associated with pyroclastic materials and volcanic rocks.
https://www.ijgeophysics.ir/article_54708_207d5d630c866dc9876e4456f0a48fda.pdf
ژئوفیزیک
مغناطیسسنجی
مقاومت ویژه
گوتیت
معدن موئیل
Geophysics
Magnetic method
resistivity method
Goethite
Moeil mine
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
128
143
54709
Research Article
کاهش نوفه تصاویر نجومی با استفاده از معادلات مشتقات جزئی
Nonlinear diffusion noise reduction in astronomical images
محمدعلی شریفی
sharifi@ut.ac.ir
1
سعید فرزانه
farzaneh@ut.ac.ir
2
مونا کوثری
mona.kosary@ut.ac.ir
3
دانشکده مهندسی نقشهبرداری و اطلاعات مکانی، پردیس دانشکدههای فنی، دانشگاه تهران، تهران، ایران
دانشکده مهندسی نقشهبرداری و اطلاعات مکانی، پردیس دانشکدههای فنی، دانشگاه تهران، تهران، ایران
دانشکده مهندسی نقشهبرداری و اطلاعات مکانی، پردیس دانشکدههای فنی، دانشگاه تهران، تهران، ایران
استفاده از سیستمهای بینایی مبنای نجومی بهعنوان روشی ارزان و مناسب بهمنظور تعیین مختصات نقاط میتواند بهعنوان روشی کمکی و همچنین جایگزین برای سیستمهای تعیین موقعیت جهانی در نظر گرفته شود. علاوه بر آن با استفاده از این سیستم میتوان حرکت ماهوارههای جاسوسی را نیز رصد نمود. با توجه به استفاده این سیستمها از تصاویر رقومی، کیفیت حاصله نقش مهمی در کیفیت خروجی نهایی خواهد داشت. نوفههای ایجاد شده در مرحله اخذ تصویر ازجمله اثراتی هستند که باعث مخدوش شدن کیفیت تصویر میشوند. در این مقاله الگوریتم جدیدی با الهام از روشهای مبتنی بر معادلات نفوذ گرما برای کاهش نوفه تصاویر نجومی ارائه گردیده است. از ویژگیهای مهم این روش جلوگیری از تغییر مکان عوارض در تصاویر یا بهعبارتدیگر حفظ ساختار تصویر است زیرا تشخیص بسیار دقیق مراکز ستارگان تأثیر بسزایی در انطباق آنها با کاتالوگ خواهد داشت. این روش بر این اصل استوار است که روشنایی لبههای موجود در تصاویر واقعی بهتدریج تغییر مییابد. بر این اساس برای حل عددی معادلات نفوذ گرما بهمنظور حفظ لبهها، لازم است همسایگیهایی با عمق بیش از یک پیکسل در نظر گرفته شود. با این کار موقعیتهای واقعی حفظ میشوند. این الگوریتم بر روی تصاویر مختلف آزمایش شده و عملکرد آن با الگوریتمهای موجود مقایسه شده است. نتایج عددی حاصل از مقایسه چهار معیار بین الگوریتم پیشنهادی با فیلتر میانگین و میانه نشان میدهد که روش ارائه شده در این مقاله، از عملکرد بهتری برخوردار است.
From ancient times, celestial bodies have been used by travelers and scientists for positioning and routing. By developing sciences, it was found that the celestial bodies form an accurate inertial system to use in navigation applications. In this system, each star is considered as a reference point in determining reference coordinate frame of the system. Due to the visibility of the satellite motion trace and the fundamental need to determine and modify satellites’ orbital parameters as well as to identify and locate espial satellites, determining the positional parameters of the satellite is also one of the modern and important applications of vision-based astronomical systems. In the modern vision-based astronomical systems, data collection is done using charge-coupled device (CCD) array. During the process of light collision to the surface of the CCD and then reading and measuring the number of photoelectrons as well as converting them to the digital numbers to store them as grey degree in each pixel, the smallest mistakes that result in lost or added electrons on each pixel can lead to distortion and noise in the image. The process of noise elimination should not only eliminate or reduce the noise but also avoid blurring the image and removing or relocating the edges of the image. To determine the primary orbit of the satellite using an optical method, the streak of the satellite must be extracted accurately because the misdiagnosis of the beginning and end points of the streak directly affects the accuracy of the determined orbit. Therefore, we need to find noise elimination methods that impose the lowest possible effects to the key complications of the astronomical images such as star and satellite streak. In this study, it is attempted to eliminate the noises using diffusion equation and solving it numerically. On the other hand, to identify the accurate position of the edges, the gradient is calculated by through convoluting the main image by the Gaussian filter. In this study, a numerical method is used to solve diffusion equation. Heat diffusion equation is an iteration-based method. It is obvious that the more the paces and iterations in the equation, the smoother the image. This factor must be chosen such that the image brightness does not exceed the main range. For this purpose, the noise must be eliminated from the image by choosing appropriate number of iterations. In this research, the structural similarity index (SSI) is used to select the optimum number of iterations. As a result, in this research, it is attempted to use noise elimination methods that impose the lowest changes to the existing satellite’s streak.
https://www.ijgeophysics.ir/article_54709_55c13f7991792283e8f8b8794b7e5857.pdf
تشخیص مراکز ستارگان
معادلات نفوذ گرما
کاهش نوفه
رد ماهواره
star centroiding
nonlinear diffusion
Noise Reduction
satellite streak
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
144
156
54710
Research Article
مدلسازی عددی دادههای رادار نفوذی زمین (GPR) با استفاده از روش اجزاء محدود
Numerical modeling of round-penetrating radar (GPR) usingfinite-element method
سجاد زارعی
1
امین رحیمی دلخانی
2
نوید امینی
navidamini@ut.ac.ir
3
گروه فیزیک زمین، مؤسسة ژئوفیزیک، دانشگاه تهران، تهران، ایران
دانشکده مهندسی نفت، دانشگاه صنعتی امیر کبیر، تهران، ایران
موسسه ژئوفیزیک، دانشگاه تهران، تهران، ایران
در این پژوهش از روش اجزاء محدود (FEM) بهمنظور مدلسازی امواج الکترومغناطیس استفاده شده است. با توجه به قابلیتهای این روش، ابتدا معادلات ماکسول در حیطه مکان گسستهسازی میشوند، سپس شرایط مرزی بهمنظور جذب امواج در کرانههای مدل اعمال میشود که از روش مرز جاذب مرتبه اول کلایتون و انگکوئیست استفاده شده است. در روش FEM عبارت مرز یک جمله جداگانه میباشد، به همین دلیل اعمال شرایط مرزی در این روش بسیار آسانتر از روش تفاضل محدود (FDM) است. پس از گسستهسازی مکانی با استفاده از روش FEM، گسستهسازی زمانی معادلات با استفاده از روش تفاضل محدود مرکزی صورت میگیرد. گسستهسازی زمانی معادلات، حجیمترین و زمانبرترین بخش محاسبات در مدلسازی هستند که نحوه گسستهسازی مکانی نقش بسزایی در این فرآیند ایفا میکند. با توجه به تُنُک بودن و متقارن بودن ماتریسهای تشکیل شده در روش FEM، درصورتیکه از الگوریتمهای بهینه بهمنظور محاسبات و ذخیرهسازی ماتریسها در این روش استفاده شود، زمان و هزینه محاسباتی بهطور قابل ملاحظهای کاهش خواهد یافت که در این تحقیق چند تکنیک بهمنظور کاهش حجم و زمان محاسبات در نرمافزار متلب ارائه شده است. فرمولها و روابط ارائه شده در این تحقیق به شکل ماتریسی هستند که بهراحتی در نرمافزار متلب قابل کد نویسی میباشند. بهمنظور بررسی روش FEM در مدلسازی دادههای GPR، الگوریتمهای توسعه داده شده بر روی مدلهای زمینشناسی فرضی آزمایش شده است که نتایج حاصل از مدلسازی دارای دقت قابل قبولی هستند.
Ground-penetrating radar (GPR) is a popular geophysical method for high-resolution imaging of the shallow subsurface structures. Numerical modeling of radar waves plays a significant role in interpretation, processing, and imaging of GPR data. A number of different approaches have been presented for the numerical modeling of GPR data. The most common approach for GPR modeling is the finite-difference method (FDM) because the FDM approach is conceptually simple and easy to program. The difficulties in applying boundary conditions at non-linear boundaries and the lack of sufficient accuracy in complex geometries are the most important drawbacks of FDM.
This paper presents a finite-element method, for simulation of ground-penetrating radar (GPR) in two dimensions in the time-domain. FEM is a powerful and versatile numerical technique for handling problems involving complex geometries and inhomogeneous media. The technique is based on a weak formulation of Maxwell’s equations. In the FEM method, the wavefield is discretized on the elements using Lagrange interpolation, and integration over an element is accomplished based upon the Gaussian-quad integration rule. The major difference between the various numerical methods is in the spatial discretization. In the elemental-based methods, the complex geometry of the problem is divided into several smaller and simpler elements, then the integrals are calculated for each element. These methods have no with any regular or irregular geometry. The responses of the model in the finite-element methods are approximated in nodal points, so nodal polynomials of Lagrange are used for interpolation of the model response. Besides, the systematic generality of the method makes it possible to construct general-purpose computer programs for solving a wide range of problems. In this paper, at first, Maxwell’s equations are discretized, then the boundary condition is applied to minimize artificial reflections from the edges of the computation domain. Although the governing equations and mechanisms are completely different between radar and seismic waves, most of GPR data processing approaches are derived from seismic data processing. Due to similarities in these two techniques, accordingly, we implement the first-order Clayton and Engquist absorbing boundary conditions (firstorder CE-ABC) introduced in the numerical finite-difference modeling of seismic wave propagation. This boundary condition is simple to apply. The presented formulations are in matrix notation that simplifies the implementation of the relations in computer programs, especially in MATLAB application. After spatial discretization with FEM, time discretization is done by Finite-Central Difference (FCD). The time discretization is the most massive and time-consuming step in modeling, which spatial discretization has an important role in this process. The stiffness, mass and damping matrices are sparse and symmetrical in FEM; so if we use the optimized numerical algorithms and storages strategies, computational costs and processing-time can be reduced significantly. To investigate the efficiency of FEM, the computer program has been written in MATLAB and has been tested on two models. The results show that the radar wave simulation via FEM is an accurate and effective approach in complex models.
https://www.ijgeophysics.ir/article_54710_c85154fac6e5c3cd4811aa7db6783ab8.pdf
رادار نفوذی زمین (GPR)
روش اجزاء محدود (FEM)
مدلسازی عددی
شرایط مرزی جاذب (ABC)
ground-penetrating radar(GPR)
finite-element method (FEM)
Numerical modeling
absorbing boundary conditions (ABC)
per
انجمن ملی ژئوفیزیک ایران
مجله ژئوفیزیک ایران
2008-0336
2783-168X
2017-11-22
11
3
157
173
51828
Research Article
بررسی هندسه مرزهای ناپیوستگی در مکران ساحلی با استفاده از مهاجرت و مدلسازی توابع گیرنده
Geometry of deep velocity discontinuities in the coastal Makran using
receiver functions
توحید نوزاد خلیل
t.nozad@gmail.com
1
سیدخلیل متقی
kmotaghi@iasbs.ac.ir
2
دانشکده علوم زمین، دانشگاه تحصیلات تکمیلی علوم پایه، زنجان، ایران
دانشکده علوم زمین، دانشگاه تحصیلات تکمیلی علوم پایه، زنجان، ایران
منطقه مکران بخشی از فرورانش فعال حاصل از همگرایی صفحات عربی و اوراسیا است. زاویه مخروطی کم و ضخامت زیاد رسوبات در منشور افزایشی از ویژگیهای متمایز کننده این منطقه هستند. در این مطالعه، با استفاده از دادههای جمعآوری شده در منطقه مکران ساحلی توسط ایستگاه باند پهن چابهار توابع گیرنده P و S و با استفاده از چهار ایستگاه کوتاه دوره، تابع گیرنده P محاسبه شد. مهاجرت به عمق توابع گیرنده P نشان میدهد که یک مرز (احتمالاً مرز موهو) در عمق 27 کیلومتری در زیر ایستگاه چابهار قرار دارد که شیبی اندک، در حدود 5/2 درجه، رو به شمال دارد. بررسی تغییرات دامنه فازهای تبدیلی بر روی مؤلفههای شعاعی و مماسی تابع گیرنده P در ایستگاه چابهار وجود یک لایه ناهمسانگرد تا عمق حدود 9 کیلومتر که احتمالاً مرتبط با رسوبات فلسی شکل در گوه برافزایشی است را تأیید میکند. مقدار عمقهای تخمینی 9 و 27 کیلومتر برای رسوبات و مرز موهو از مدلسازی همزمان تابع گیرنده و منحنی پاشش محاسبه شده است. محاسبه و برانبارش توابع گیرنده S برای ایستگاه CHBR یک پالس منفی با تأخیر زمانی 10 تا 13 ثانیه (متناسب با عمق 80 تا 100 کیلومتر) را در زیر این ایستگاه نشان میدهد که میتواند مرتبط با مرز لیتوسفر اقیانوسی با استنوسفر باشد.
We analyzed the teleseismic data gathered by a broad-band (CHBR) and four short-period (CDK, CNT, KHB, KSM) seismometers, located in western coastal Makran, north of Chabahar, Iran. The data were gathered by the roughly north-south direction quasi-linear profile and used to calculate P (for all stations) and S (only for CHBR) receiver functions utilizing iterative deconvolution technique of Ligorria and Ammon (1999). Because of backazimuth gaps in south and western directions, we used PKiKP and Pdiff phases to calculate receiver functions in a similar processing approach. Calculated P receiver functions are migrated to depth to clarify the geometry of velocity boundaries at the base of sediments and Moho. The result shows that there is a dipping interface lying at a depth of 27 km (beneath CHBR) to 31 km (beneath CDK), which imply a 2.5o dipping Moho boundary beneath the study region. To avoid the trade-off between velocity model and reported depth, we jointly modeled the stacked receiver function, and group velocity dispersion curve for CHBR and the output model was considered for any time to depth migration of receiver functions.
We analyzed the effects of P and S anisotropy on teleseismic converted waves to map the presence, the strike, and the depth of anisotropic structures. High-resolution PRFs are considered for such analysis. The following criteria are considered to select the high-quality receiver function (Schulte-Pelkum and Mahan, 2014): the signal-to-noise ratio of the three components of the seismograms is at least 1.5; the convolution of the PRF with the vertical component of the seismogram reproduces at least 60% of the horizontal component (defined as variance reduction by Ligorria and Ammon, 1999); the PRF shows a positive polarity direct P arrival; the receiver function amplitude does not exceed 1; any arrivals’ pulse length does not exceed 3.5 s. The latter two criteria are employed because very high amplitudes and long oscillatory pulses are typical characteristics of an unstable deconvolution (Schulte-Pelkum and Mahan, 2014). The calculated PRFs were then binned in 5° azimuthal groups with 5° overlap. In CHBR station, we recognized signs of the top (at 1 km depth) and bottom (at 9 km depth) of an anisotropic layer with almost north-south anisotropic symmetry axis. In addition, we recognized a flat interface beneath CHBR station at 27 km depth that is not in consistency with the result of migration to a depth of RFs showing a 2.5o dip Moho at the same place. For this reason, we utilize forward modelling to calculate synthetic PRFs to explain periodic amplitude variation of P to S converted phases with back-azimuths in each station that could be a signature for anisotropic velocity features. The forward modeling indicates that the horizontal interface makes a similar pattern on simulated PRfs as a low angle dipping interface with dip less than 10o.
Migration of S receiver functions reveals a deep velocity discontinuity at depth around 80 to 100 km that might be considered as a shallow lithosphere-asthenosphere boundary beneath the study region.
https://www.ijgeophysics.ir/article_51828_d4711e2773e37da92d15ba671cd8115d.pdf
مکران ساحلی
تابع گیرنده
ساختار عمیق
coastal Makran
receiver function
velocity discontinuity
Moho