انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Random noise attenuation in seismic data by non-local mean filteringتضعیف نوفههای تصادفی دادههای لرزهای با استفاده از فیلتر میانگین غیرمحلی11039519FAفرید کریمیدانشکده معدن، نفت و ژئوفیزیک، دانشگاه صنعتی شاهرود، ایرانامین روشندل کاهودانشکده معدن، نفت و ژئوفیزیک، دانشگاه صنعتی شاهرود، ایران0000-0002-2214-2558Journal Article20161116The reliability of seismic mapping is strongly dependent upon the quality of the records. Seismic records are usually affected by various types of noise such as ground rolls, multiples, random noise, reflection and reflected refraction from near surface structures, and so on. Random noise resulted from random oscillation during acquiring data is one of the most important and harmful noises that exist in seismic data over entire time and frequency. Random noise attenuation is an important step in seismic data processing affecting the data interpretation.
Many efforts have been made to remove this type of noise from the seismic data. The predictive filter is an ordinary method commonly used for random noise attenuation from seismic data. This filter can be used in various domains, such as the f-x domain (Haris and White, 1997) and the discrete cosine domain (Lu and Liu, 2007). Jones and Levy (1987) removed events which were not coherent trace-to-trace events by means of the Karhunen-Loeve transform. Karsli et al. (1996) applied complex trace analysis to seismic data to suppress random noise and improve the temporal resolution of the data. Terickett (2008) attenuated the random noise from the seismic data by Cadzow filtering of a constant frequency slice in an f-xy domain. Bekara and Baan (2008) attacked the random noise problem using the empirical mode decomposition technique.
In this study, we used non-local means filtering (Buades et al., 2005) developed for image denoising for random noise suppression in seismic data. In this method, a seismic record can be considered as an image. The non-local means method is based on the assumption that the image content is likely to repeat itself within some neighborhood. The neighborhood of a pixel is generally chosen to be a square with a dimension size of three to nine centered upon the pixel of interest. However, the size and shape of the neighborhood can vary. In this method, for each pixel with a neighborhood , the pixels with the neighborhoods similar to that of the interest pixel are found. The denoised value of a pixel is determined by a weighted average of all the pixels in the image. The weight of each pixel can be calculated by the similarity between the two pixels and computed using the Gaussian weighted Euclidean distance between the neighborhood around the pixel and the neighborhood around the pixel .To investigate the efficiency of the proposed method, we tested the non-local means algorithm on both synthetic and real seismic data. We also compared the obtained results by that of the traditional mean and median algorithm for seismic data denoising. To investigate further, we applied the three denoising methods to synthetic seismic data with different amounts of signal-to-noise ratios. After analyzing all of the results, the non-local means algorithm proved to be a better algorithm for seismic data denoising and had the best performance among the three methods.
دادههای لرزهای خام همواره شامل انرژیهای ناخواسته هستند. از این انرژیهای ناخواسته و نامطلوب در لرزهشناسی با عنوان نوفه یاد میشود. تضعیف نوفه تصادفی یکی از مراحل مهم در پردازش دادههای لرزهای است که تاثیر بسزایی در تعبیر و تفسیر دادهها دارد. تابهحال روشهای گوناگونی برای تضعیف نوفه تصادفی عرضه شده است که هرکدام مزایا و معایب مربوط به خود را دارد. در این مقاله از فیلتر میانگین غیرمحلی که بهمنظور تضعیف نوفه از تصاویر در فن پردازش تصویر توسعه یافته است، برای تضعیف نوفههای تصادفی در دادههای لرزهای استفاده میشود. این الگوریتم مقدار تضعیف نوفه هر پیکسل داخل تصویر را با استفاده از نمونهها یا پیکسلهای مشابه آن، بدون در نظر گرفتن موقعیت آنها نسبت به یکدیگر تعیین میکند. در واقع این الگوریتم، یک فرایند غیرمحلی است. نتایج بهدست آمده در مقایسه با دو روش متداول فیلتر میانگین و فیلتر میانه که بهصورت محلی عمل میکنند، حاکی از برتری روش فیلتر میانگین غیرمحلی در تضعیف نوفههای تصادفی در دادههای لرزهای است.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Lithofacies estimation by multi-resolution graph-based clustering of petrophysical well logs: a case study from one of the Persian Gulf’s gas fields, Iranتعیین رخسارههای سنگی با استفاده از خوشهسازی براساس نمودار با تفکیکپذیری چندتایی (MRGC) نمودارهای پتروفیزیکی چاه: بررسی یکی از میدانهای خلیج فارس113039520FAمنصور آقچه لوشرکت ملی مناطق نفت خیز جنوب، اهواز، ایرانحمیدرضا همتی آهوییدانشگاه بین المللی امام خمینی (ره)، قزوین، ایرانمجید نبی بیدهندیموسسه ژئوفیزیک دانشگاه تهران، ایرانعلی اکبر رحیمی بهارپژوهشگاه صنعت نفت، تهران، ایرانJournal Article20161116Located in the Persian Gulf, the gas field studied in this research is one of the largest gas fields in the world. Its gas-in-place is estimated to be about 14.2 trillion cubic meters while amount of its condensate-in-place might be around 18 billion barrels. This gas field has also an oil layer containing about 6 billion barrels of oil-in-place. In this study, Kangan and Dalan formations of this field were considered. Kangan formation has three main facies: clean carbonate facies, basic clay and shale facies, and evaporate carbonate facies. Dalan formation contains four facies: shore restricted carbonates, shore organic carbonates, carbonates of the open sea, in-shore carbonate-clastic. These two formations have gas & condensate fluids. Since this field is a heterogeneous carbonate system, lithofacies characterization is the best solution for overcoming the problem of heterogeneity in determining the petrophysical properties of the reservoir rock, reservoir modeling and identifying producing zones. However, coring as the most robust method of lithofacies identification is very expensive, time consuming and limited to a few number of wells. Therefore, this study is focused on determining the lithofacies of the study formations from available well logs.
For this purpose, multi-resolution graph-based clustering (MRGC) technique which is a dot-pattern recognition method based on non-parametric K-nearest neighbor and graph data representation was applied to sonic, density, neutron porosity and gamma ray logs to define electrofacies similar to core-derived facies determined as eight distinct facies.The cluster of the MRGC method is defined from a model with a specific character associated with the group of lithofacies. Then, Kernel representative index was used to calculate the optimal number of clusters. Small facies groups were formed based on utilizing the neighboring index to determine a K-nearest neighbor attraction for each point. At last, final clusters were constructed by combining the small clusters which lead to identifying eight facies of this gas field from well logs of high accuracy.
When electrofacies of one of the wells is built basd on its lithofacies, its cross-plots will be plotted and the certainty of electrofacies with respect to lithofacies will be checked. If the model is acceptable, it is applied to the data from two other wells and their electrofacies will be obtained. For testing facies, the cross-plots of these two wells were also drawn and painted based on facies. If there are similar petrophysical properties for each facies, the model created in the wells without cores is confirmed. MRGC is a fast method that allows the geologist or petrophysist to analyze and test different combinations of data in a short amount of time. It is also not limited by the dimensions of the data and number of the clusters. The method used in this study has obviated the need for extensive coring in this field which caused saving large amounts of money and time; and it can help to optimize the determination of new well locations and optimum pay zones.
میدان مورد بررسی در این پژوهش، یکی از بزرگترین میدانهای گازی جهان است که در خلیج فارس قرار دارد. ازآنجاکه سنگ مخزن این میدان کربنات نامتجانس است، تعیین رخسارههای سنگی آن میتواند راهحلی کارآمد برای طبقهبندی رخسارههای مخزن در جهت تعیین خواص پتروفیزیکی سنگ مخزن، مدلسازی مخزن و تفکیک کردن مناطقی که پتانسیل بالاتری برای تولید نفت دارند از مناطق کمپتانسیل باشد .دقیقترین روش تعیین رخساره استفاده از مغزه است اما مغزهگیری بسیار گران و وقتگیر است و بهصورت محدود در برخی چاهها قابل دسترسی است. بنابراین در این تحقیق بر تعیین رخساره سنگی با استفاده از دادههای نمودارهای چاه پرداخته شده است. در این راستا روش خوشهسازی براساس نمودار با تفکیکپذیری چندتایی (MRGC) که روشی براساس روش ناپارامتری K- امین همسایه نزدیک و نمایش نموداری داده است، بر روی نمودارهای صوتی، تخلخل نوترون، چگالی و پرتو گاما اِعمال شده تا رخساره الکتریکی معادل رخساره سنگی بهدست آمده از مغزه، بازسازی شود. خوشههای روش MRGC از مدلی که دارای خواص مشخص وابسته به گروهی از رخسارههای سنگی است، تعریف میشود. گروههای رخسارههای کوچک با استفاده از شاخص همسایگی برای تعیین جاذبه همسایة نزدیک kاُم برای هر نقطه، شکل میگیرند. درنهایت خوشههای نهایی با تجمع خوشههای کوچک شکل میگیرند که منجر به شناسایی 8 رخساره این میدان گازی از نمودارهای چاه با دقت زیاد میشود. روش مورداستفاده در این تحقیق نیاز به مغزهگیری گسترده در این میدان را رفع میکند و منجر به صرفهجویی بسیار در هزینه و زمان میشود.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Climatological study of Asian and European blockings by means of two indices in the period of 1950 -2010بررسی اقلیمشناختی بندالهای آسیا و اروپا با دو شاخص در دوره 1950- 2010315139521FAسارا فهیمیدانشگاه آزاد اسلامی، واحد علوم و تحقیقات تهران، ایرانفرهنگ احمدی گیویموسسه ژئوفیزیک دانشگاه تهران، ایران0000-0002-9487-4862مجید مزرعه فراهانیموسسه ژئوفیزیک دانشگاه تهران، ایرانJournal Article20161116Blocking is one of the most important weather systems in mid-latitudes that have significant effects on atmospheric air flow and regional weather conditions. In this study, climatology of Asian and European blockings was investigated for the period of 1950–2010. The data used in this study were 500-hPa geopotential heights from the National Centers for Environmental Prediction/ National Center for Atmospheric Research (NCEP/NCAR) reanalysis dataset with 2.5°×2.5° resolution. The blocking events were detected using two indices; the first one, a one-dimensional index, was the new version of the Tibaldi and Molteni (1990) index presented by Barriopedro et al. (2005) and the second one, a two-dimensional index, was presented by Diao et al. (2006). Then, applying these indices, the characteristics of detected blockings such as frequency, duration, intensity, position of center and extension were obtained and compared with each other.
In the period of this study, 1054 and 729 blocking events were identified using the first and second indices, respectively. The spatial and temporal distribution analysis showed that blocking was most prevalent in the longitude band between 30 °W and 47.5 °E and most of the maximum intensity occurrence was in the east of Atlantic sector between 30 °W and 5 °E.
It was found that blocking events were more frequent in the Atlantic-European sector than in the other sectors in winter, spring and autumn, but more persistent in the Atlantic-European sector than in the Asian sector for all seasons. Additionally, the maximum amplitude of intensity was 3.0 for the first index and 150 gpm for the second index. A comparison of the blocking events in the Asian sector with the Atlantic-European sector indicated that the intensity of winter events in the European sector was higher than that of the other seasons. This may be due to synoptic-scale eddies and planetary waves which are more active and stronger in winter compared to the other seasons.
The trend analysis displayed that the number of blocking events in the period of this study was decreasing slightly in the entire region, including Iran. The number of events was significantly more in the west and northwest of Iran than in the central and eastern parts.
Also, the results of two indices showed that the latitude band of blocking occurrence and its maximum frequency was similar in the different seasons for the first index, whereas this band changed for the second index. The maximum frequency of blocking event was seen in latitude 55 °N for winter, spring and autumn and in latitude 62.5 °N for summer. Thus, it indicated that blockings occurred in a higher-latitude band in summer than the other seasons. In addition, it was found that there was a good consistency between the longitude band of blocking event, including the Asian and European sectors, and its frequency in summer and autumn using the two indices. In the other seasons, the longitude band of the Asian and European blocking events, especially in the eastern part of the Asian sector, was different for the two indices. The other noticeable point was that the activity centers of blocking events in winter, spring, and autumn obtained by the second index were located more westward of the Atlantic-European sector and more eastward of the Asian sector than those obtained by the first index.
The main conclusion was that since the Diao et al. (2006) index is a two-dimensional one, it is able to identify better and more precisely the blocking events and their characteristics as compared to the Barriopedro et al. (2005) index.
در این مقاله، بررسی اقلیمشناختی پدیده بندال طی یک دوره 61 ساله (1950-2010) در بخشی از نیمکره شمالی، شامل آسیا و اروپا و قسمت شرقی اقیانوس اطلس با بهکارگیری دادههای شبکهای NCEP-NCAR صورت گرفته است. با استفاده از دو روش آشکارسازی بندال یکبُعدی و دوبُعدی، بهترتیب تعداد 1054 و 729 رویداد بندال در دوره زمانی و ناحیه موردنظر، شناسایی شد. براساس نتایج این بررسی، تعداد رویدادهای بندال این دوره دارای یک روند نسبتاً کاهشی برای کل منطقه مورد بررسی و منطقه ایران است. همچنین توسط این دو روش، یعنی شاخصهای باریو پدرو و همکاران (2005) و دیااو و همکاران (2006)، ویژگیهای رویدادهای بندال از قبیل بسامد، طول عمر، محل وقوع، شدت و مراکز آنها بررسی و نتایج ویژگیهای اقلیمی، با یکدیگر مقایسه شدهاند. طبق بررسیهای صورت گرفته، بیشترین وقوع بندالها در منطقه مورد بررسی از نظر توزیع مکانی و در همة فصلها، در محدوده طول جغرافیایی E° 5/47 تا °W30 و از نظر وقوع بیشینه شدتها روی شرق اطلس واقع در E ° 5 تا W° 30 قرار دارد. از طرفی، در فصلهای زمستان و بهار بسامد وقوع بندال و دوره استمرار آن مقادیر بیشتری دارند، درحالی که در فصل تابستان این مقادیر کمینه است. نتایج بسامد وقوع در عرضهای جغرافیایی و فصلهای گوناگون نشان میدهد که برای شاخص اول، عرض جغرافیایی ثابت همراه با بیشینههای بدون تغییر در فصلهای گوناگون است، اما برای شاخص دوم، این بازه متغیر است و محل بیشینه وقوع بندالها در عرض N° 55 در فصلهای بهار، زمستان و پاییز، و در فصل تابستان با جابهجایی به سمت عرضهای بالاتر، در عرضN ° 5/62 قرار دارد. بهعلاوه، در مقایسه نتایج دو شاخص دیده میشود که بندالهای شناسایی شده در بازه طول جغرافیایی با شاخص دوم، دارای درصد بسامد وقوع کمتری نسبت به شاخص اول در فصلهای زمستان و بهار هستند، حالآنکه الگوهای مربوط به دو شاخص در فصلهای تابستان و پاییز تقریبا مشابه هستند.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Depth detection of magnetic bodies by using the analytic signal derivativeبرآورد عمق تودههای مغناطیسی با استفاده از مشتقات سیگنال تحلیلی5236339522FAمسلم فاتحیدانشکده مهندسی معدن دانشگاه صنعتی اصفهان، ایرانغلامحسین نوروزیدانشکده مهندسی معدن، پردیس دانشکدههای فنی، دانشگاه تهران، ایرانفاطمه حاجیئیشرکت مجریان توسعه معادن آسیا، ایرانJournal Article20161116The analytic signal and its derivatives have been used in magnetic data interpretation such as estimation of the edge, depth and slope of the magnetic bodies for several years. Nabighian (1972, 1974, 1984) extensively described the principles of the analytic signal method for the location of the 2D sources. The analytic signal shape can be used to determine the depth of the magnetic sources. Atchuta et al (1981) used the anomaly width at half the amplitude to derive the depths. W. R. Roest and his collaborators used the 3D analytic signal amplitude to estimate the magnetic source depth (Roestetal.,1992) and to identify remanent magnetization (Roest and Pilkington, 1993). However, their applications and many conclusions are based on a 2D vertical- magnetic-contact model assumption. Further implementation of this technique was made by Hsu et al. (1996) who developed an enhanced analytic signal applied to a higher order vertical derivative of potential - field anomalies and thus providing a better visualization of the outlines of magnetic bodies. This enhanced analytic signal was adapted as an automatic interpretation tool (Debeglia and Corpel, 1997; Hsu et al. ,1998). Hsu et al, 1998 used the enhanced analytic signal to estimate the depth of step models and the thin dike. They also proposed an algorithm for the type of these models afterwards. Bastani and Pedersen (2001) have developed an algorithm for the automatic estimation of the source parameters, including the dip and susceptibility contrast, from the analytic signal in the case of magnetic profiles. Salem and Ravat (2003) proposed a combined method (AN-EUL) based on the Euler equation and the analytic signal. Li (2006) discussed a 3D analytic signal and proposed that, generally, in three dimensions, the analytic signal is dependent on everything that the total magnetic intensity (TMI) itself may depend on i.e. the direction of the inducing field, the direction of the remanent magnetization, the dipping angle of the source body, and the depths to the top and bottom of the source body.
In this study, the equations of analytic signal and its derivatives were developed based on their effective parameters for contact, thin dikes and horizontal cylinder models. Afterwards, using the maximum value of the ratio of analytic signal and its derivatives, an equation for depth and structural index estimation was obtained which was exactly similar to the equations in the AN–EUL method. Then, assuming a priori knowledge about the shape of magnetic body, other equations were obtained for the depth detection of contact, thin dike and horizontal cylinder. To evaluate the precision of these methods, we have used a dike model with constant width located at deferent depths. In this modeling study, we proved that the estimated depth from the AN-EUL method was dependent on the depth-to-width ratio. Using the AN-EUL method, it is not possible to detect the depth of bodies with their depth-to-width ratio less than two. To solve this problem, we used proposed equations for different models. This method was used for the depth detection of the Aliabad iron deposit, and its results were compared to the results from drilling exploration data. The correlation coefficient between the exact depth and the estimated depth was equal to 85 percent.
سالها است که از سیگنال تحلیلی و مشتقات آن در تفسیر دادههای مغناطیسی اعم از برآورد مرز، عمق و شیب تودهها استفاده میشود. در این نوشته رابطه سیگنال تحلیلی و مشتقات آن برحسب پارامترهای موثر بر آن برای تودههای گوناگون از قبیل کنتاکت، دایک نازک و استوانه افقی، توسعه داده شده است. با استفاده از نسبت بین بیشینه مقدار سیگنال تحلیلی ساده و مشتقات مرتبه اول و دوم آن، روابطی برای برآورد عمق تودهها و شاخص ساختاری آنها بهدست آمد که دقیقاً مشابه با روابط برآورد عمق و شاخص ساختاری روش تلفیق سیگنال تحلیلی و اویلر (AN-EUL) هستند. همچنین با دانش قبلی نسبت به مدل توده، روابط دیگری برای برآورد عمق هرکدام از مدلهای کنتاکت، دایک نازک و استوانه افقی بهدست آمد. برای بررسی دقت این روابط، از یک مدل دایکیشکل در عمقهای متفاوت استفاده شد. وابستگی روش تلفیق اویلر و سیگنال تحلیلی به نسبت عمق توده به عرض آن، نتیجه این بررسیها بود. بنابراین وقتی نسبت عمق توده به عرض آن کمتر از 2 باشد نمیتوان از این روش برای برآورد عمق استفاده کرد. برای حل این مشکل از روابط برآورد عمق عرضه شده برای مدلهای متفاوت استفاده شده است. از این روش برای برآورد عمق دادههای مغناطیسی کانسار آهن علیآباد واقع در استان زنجان استفاده شد و نتایج آن با اطلاعات حفاریهای اکتشافی مورد مقایسه قرار گرفت. ضریب همبستگی بین عمق برآوردی و عمق واقعی برابر با 85 درصد شد.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Statistical-dynamical analysis of the mutual effects of NAO and MJOتحلیل دینامیکی–آماری اثر متقابل نوسان اطلس شمالی (NAO) و نوسان مدن-جولیان (MJO)648039526FAشاهین عالمزادهمؤسسۀ ژئوفیزیک دانشگاه تهران، ایرانفرهنگ احمدی گیویمؤسسۀ ژئوفیزیک دانشگاه تهران، ایران0000-0002-9487-4862علیرضا محبالحجهمؤسسۀ ژئوفیزیک دانشگاه تهران، ایرانمحمدعلی نصراصفهانیدانشگاه شهرکرد، شهرکرد، ایران0000-0001-6249-3355Journal Article20161116Vast areas of the Northern Hemisphere are influenced by the North Atlantic Oscillation (NAO) and the Madden-Julian Oscillation (MJO). Revealing the interaction between them can help us predict the phase and amplitude of each one based on the conditions of the other.Statistical-dynamical study on critical events of MJO and NAO during the Northern Hemisphere winter in a 37-year period (1974-2011), reveals a considerable connection between NAO and MJO. In a statistical method, the computed mean value of NAO index for different phases of MJO associated with 5 to 25-day lags, showed a causal relation between the two phenomena. Phase 7 (3) of MJO was observed 20 to 25 days after the positive (negative) phase of NAO. On the other hand, 5 to 15 days after phase 3 (7) of MJO, positive (negative) phase of NAO occured.
Lagged composite maps of 200-hPa zonal wind anomaly and the outgoing long-wave radiation (OLR) anomaly were computed for different phases of critical NAO in order to analyze the dynamical impact of NAO on MJO. Within 5 to 25 days after the positive NAO, maximum 250-hPa westerly anomalies in the Africa-Asia jet region propagate eastward from Africa to central Indian Ocean. On the other hand, OLR anomaly patterns propagate eastward in the same period and same region after the positive NAO that might be induced by the propagation features of 250-hPa zonal wind anomalies. On the 25-day lag relative to the critical positive NAO, both 250-hPa westerly anomalies over the central and eastern Indian Ocean and the development of positive OLR anomalies in the same region were consistent with the conditions of phase 7 of MJO. However, on the 5 to 25-day lags relative to the critical negative NAO, 250-hPa zonal wind anomalies showed the attenuation of the westerly Africa-Asia jet leading to the favorable conditions for evolution of the easterly anomalies over the Indian Ocean. On the other hand, the OLR anomaly patterns displayed either little or no westward propagation. On the 25-day lag relative to the critical negative NAO, the development of both 250-hPa easterly anomalies and negative OLR anomalies over the central Indian Ocean corresponded to the conditions of phase 3 of MJO.
Lagged composites of 500-hPa geopotential height anomaly for different phases of MJO were also computed in order to analyze the dynamical impact of MJO on NAO. On the 5 to 15-day lags relative to phase 3 of MJO, a Rossby wave train was observed over the Pacific Ocean extending to the North America and Atlantic region. The waves were oriented in the South West – North East direction. On the 10-day and 15-day lags, the waves started to gradually break anticyclonically over the North America and North Atlantic regions. On the 15-day lag, a dipolar pattern formed in the Atlantic region associated with a negative anomaly center of 500-hPa geopotential height in the northern part of the North Atlantic and a positive anomaly center in the central part of the North Atlantic, leading to the formation of the positive NAO phase. On lag days after phase 7 of critical MJO, a similar wave train was also extended from the Pacific and North America to the Atlantic region; however, the orientation of the waves was South East-North to West at this time. On the 5 and 10-day lags, a cyclonic wave breaking was observed over the North America and North Atlantic regions. Ultimately, on the 10-day lag a dipole structure similar to the negative NAO phase was formed over the North Atlantic region.
نوسان اطلس شمالی (NAO) و نوسان مدن-جولین (MJO)، بخشهای وسیعی از نیمکرۀ شمالی زمین را تحت تأثیر قرار میدهند. با بررسی برهمکنش متقابل آنها و با داشتن اطلاعات کافی از هریک، میتوان فاز و شدت دیگری را پیشبینی کرد. بررسی آماری و دینامیکی رخدادهای NAO و MJO بحرانی زمستان نیمکرۀ شمالی در دورۀ 37 ساله (1974-2011)، رابطۀ قابلتوجهی بین این دو پدیده را نشان میدهد. در روش آماری، محاسبۀ شاخص میانگین NAO به تفکیک فازهای MJO و با تأخیرهای 5 تا 25 روزه، حاکی از رابطۀ علت و معلولی بین این دو پدیده است. برای بررسی دینامیکی اثر NAO بر MJO، نقشههای ترکیبی تأخیری بیهنجاری باد مداری تراز 200 هکتوپاسکال و بیهنجاری تابش موج بلند خروجی (OLR) به تفکیک فازهای فرین NAO و برای تحلیل اثر MJO بر NAO، بیهنجاری تأخیری ارتفاع ژئوپتانسیلی تراز 500 هکتوپاسکال به تفکیک فازهای MJO بهدست آمده است. 20 تا 25 روز پس از فاز مثبت (منفی) بحرانی NAO، فاز 7 (3) MJO دیده میشود. این اثر NAO بر MJO با ایجاد تغییر در بیهنجاری باد مداری تراز 200 هکتوپاسکال و نیز بیهنجاری OLR در عرضهای حارهای از اقیانوس اطلس تا اقیانوس هند به تفکیک فازهای بحرانی NAO همراه است. تحول زمانی الگوی بیهنجاری باد مداری تراز 200 هکتوپاسکال و جهت انتشار این بیهنجاریها در روزهای پس از فاز مثبت (منفی) NAO با تحول زمانی و جهت انتشار بیهنجاری OLR همخوانی بسیار دارد. بنابراین به نظر میرسد که اثر NAO بر همرفت MJO از راه القای تغییرات باد مداری تراز بالای جوّ تا اقیانوس هند منتقل میشود. ازسویدیگر، 5 تا 15 روز پس از فاز 3 (7) بحرانی MJO، NAO به فاز مثبت (منفی) میرود. این اثر MJO بر NAO، با جهتگیریهای متفاوت قطار امواج راسبی و در نتیجه شکست موج متفاوت در امریکای شمالی و اطلس شمالی به تفکیک فازهای 3 و 7 از MJO بحرانی قابل مشاهده است.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Crustal structure of the Western Alborz by joint inversion of receiver function and surface wave dispersion curveساختارسرعتی پوسته در البرز غربی به روش برگردان همزمان تابعهای انتقال گیرنده و منحنی پاشندگی امواج سطحی819439549FAسمیرا موسویانپژوهشگاه بین المللی زلزلهشناسی و مهندسی زلزله، تهران، ایرانمحمد تاتارپژوهشگاه بین المللی زلزلهشناسی و مهندسی زلزله، تهران، ایران0000-0003-2904-6195Journal Article20161116Crustal velocity structure of a seismically active region such as the Alborz located in North of Iran has a great influence on the precise location of earthquakes, attributing the seismicity to the active faults, and improving the reliability of the seismic hazard assessment. In spite of several researches carried out on the Central Alborz crustal structure, very little is known about the structure and thickness of the crust beneath the Western Alborz which is significant due to occurrence of the 1990 Manjil-Tarom earthquake with Ms = 7.3. In this study, we intend to examine the thickness and structure of the crust beneath three stations of Zanjan, Roudbar and Ghazvin located in Western Alborz by joint inversion of the receiver functions and Rayleigh wave group velocity dispersion measurements. A combined inversion of Rayleigh wave group velocities and body wave receiver functions increases the uniqueness of the solution over separate inversions and also facilitates explicit parameterization of layer thickness in the model space. The time-domain iterative deconvolution procedure, which has higher stability with noisy data compared to frequency-domain methods, is employed to deconvolve the vertical component of the teleseismic P waveforms from the corresponding horizontal components and obtain radial and transverse receiver functions for two broadband stations of ZNJK and RUD, and one short-period station of GZV. The waveforms were corrected from the instrument response before proceeding with the receiver function deconvolution. High frequencies were filtered using a Gaussian filter, at 2.5, 1.6, and 1.0, giving an effective high-frequency limit of about 1.2, 0.8 and 0.5 Hz, respectively. As the structure may vary with azimuth and with epicentral distance, all the observations were grouped by azimuth (< 10°) and distance (∆ < 10°). To increase the signal-to-noise ratio of the deconvolved traces, the individual receiver functions were aligned according to the P-wave arrival and point-to-point stacked waveforms. The stacked receiver function was then allocated an average slowness and the back-azimuth of every event was included in the stack. Rayleigh wave group velocity dispersion came from tomographic images in a time period between 10 s and 70 s produced by a study of regional fundamental modes of Rayleigh waves propagating across Iran and surrounding regions. Fundamental-mode Rayleigh wave group velocities to each stations were taken from the corresponding tomographic cell containing the station.
The results show that the crust beneath the Roudbar station has a thickness of 36 ±3 km. A shallow low velocity sedimentary layer, of about 4 km thickness, and a velocity discontinuity at a depth of ~12 km is observed in the crustal model of this station. Beneath the Zanjan Station, the Moho depth varies from 38 to 42 km. The same sedimentary layer as beneath the Roudbar station and an interface at about 15 km depth are observed. Toward east, beneath the Ghazvin station, the thickness of the crust increases up to 52 km which is close to that proposed for the crust of the Central Alborz. Our seismological results show that the Western Alborz has a moderate crustal root but of insufficient thickness to compensate the elevation of the range.
آگاهی از ساختار سرعتی و عمق موهو در یک منطقه فعال از نظر لرزهخیزی چون زون البرز واقع در شمال ایران، عامل مهم و مؤثری در مکانیابی دقیق زمینلرزههای آن منطقه، نسبت دادن لرزهخیزی به گسلهای فعال و درنهایت بهبود نتایج تحلیل خطر زمینلرزه است. اگرچه پژوهشهای معتبر متعددی در راستای برآورد عمق موهو در البرز مرکزی صورت گرفته است، اما درالبرز غربی که به علت رخداد زمینلرزه سال 1369 رودبار-منجیل با بزرگی Ms = 7.3 اهمیت فراوانی دارد، هنوز برآورد دقیقی از ساختار و ضخامت پوسته در دست نیست. در تحقیق حاضر تلاش بر این است تا با استفاده از دادههای مربوط به زمینلرزههای دور ثبت شده در سه ایستگاه قزوین، زنجان و رودبار واقع در زون البرز غربی، ضخامت و در صورت امکان ساختار سرعتی پوسته در این زون به روش برگردان همزمان تابعهای انتقال گیرنده و اطلاعات پاشندگی سرعت گروه امواج ریلی تعیین شود. نتایج پژوهش حاضر نشان میدهد که ضخامت پوسته در زیر ایستگاه رودبار 36 (3±) کیلومتر است. یک لایه رسوبی کمسرعت با ضخامت حدود 3 کیلومتر و یک سطح ناپیوستگی درون پوستهای در عمق حدود 13 کیلومتر قابل تشخیص است. در زیر ایستگاه زنجان، عمق موهو بین 38 تا 42 کیلومتر متغیر است. لایهای کمسرعت مشابه با ایستگاه رودبار به ضخامت 3 کیلومتر و یک ناپیوستگی درون پوستهای در عمق 15 کیلومتری در زیر این ایستگاه مشاهده میشود. در حرکت به سمت شرق، در زیر ایستگاه قزوین، ضخامت پوسته تا 52 (3±) کیلومتر افزایش مییابد که در واقع به نتایج بهدست آمده در البرز مرکزی نزدیک میشود. نتایج فوق کماکان بر نبود هرگونه ریشه در زیر زون البرز غربی دلالت دارد و نشان میدهد که این منطقه در تعادل همایستایی نیست.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Seismic imaging in geologically complex thrust belts by kinematic wavefield attributesتصویرسازی لرزهای در کمربندهای رورانده با ساختار پیچیده با نشانگرهای جنبشی میدان موج9511639550FAمهرداد سلیمانیدانشگاه صنعتی شاهرود، دانشکده مهندسی معدن، نفت و ژئوفیزیکاحسان ادیبیشرکت مدیریت اکتشاف، تهرانهاشم شاهسونیدانشکده مهندسی معدن، دانشگاه سنندج، کردستانمحمد رضا سکوتیشرکت مدیریت اکتشاف، تهرانJournal Article20161116Structural interpretation in geologically complex structures is a controversial task in the field of seismic interpretation. In recent years, Seismic imaging technology continues to make remarkable progress in imaging of complex structures such as sub-salt, fault-and-thrust areas and mud volcano bearing areas. The latest advances in the areas of migration and velocity estimation are the reverse time migration, pre-stack reverse time-depth migration, common reflection surface stack, partial common reflection surface and common diffraction surface stack methods. The continuous advances in computing facilities make such data-driven approaches feasible which have increasingly gained in relevance in recent years. In these methods, the subsurface structures are imaged without deriving a complete model for the elastic properties of the layers. In this case, just an implicit knowledge of the elastic properties directly derived from data would be sufficient. Common reflection surface stack method is one of the data-based seismic imaging methods which simulate the zero offset section. It has the great advantage that it is independent from the velocity model. By paraxial ray theory, the second order travel-time equation of the common reflection surface would be derived which gives a surface-shape operator that works as the stacking surface. The basic idea of this method is to take the kinematic reflection response of a segment of the reflector with defined curvature and orientation by two hypothetical experiments providing the wavefronts of the so-called eigenwaves. One eigenwave is obtained by placing a point source at the reflector that produces the upgoing normal incidence point wave. An exploding reflector experiment yields the second upgoing eigenwave called the normal wave. The curvature and the raypath of the wavefronts of these waves are known as the parameters of the stacking operator. To simulate the zero offset section, three parameters or kinematic wavefield attributes are needed one of which is related to the emergence angle of the central ray or to the orientation of the reflector segment, while two others are related to the curvature of the reflector. Therefore, the common reflection surface stack method, unlike the common mid-point stack method, is not restricted to a subset of multicoverage data, and it also works on the full data volume. Seismic data from a complex structure in the Southwest of Iran was processed by the common reflection surface stack method to overcome some of the ambiguities of seismic imaging in such regions. In the first step, some preprocessing was done on the data to be prepared for seismic imaging. Then, the data was processed with two imaging methods. The results obtained by the conventional method showed lots of ambiguities in the final seismic section. To solve these problems, the common reflection surface stack method was applied to the data. However, the result was not promising and some sort of optimization on the method was performed. This optimization gave the best kinematic wavefield attributes in complex structures. However, due to large computation time and a large amount of sections that this method gives, the optimized common reflection surface stack was used from the beginning. The result of the latter method proved that using this method in geologically complex regions could be trusted and helpful in obtaining high quality images.
در این تحقیق، امکان تهیه مقطع لرزهای با کیفیت زیاد با استفاده از روش نوین برانبارش سطح بازتاب مشترک مورد بررسی قرار میگیرد. بدینلحاظ، بهمنظور برطرف کردن برخی از مشکلات تصویرسازی لرزهای، دادههای بهدست آمده در یکی از مناطق جنوب غرب ایران که دارای ساختارهای زمینشناسی بسیار پیچیده است، مورد پردازش قرار میگیرد. در ابتدا دادهها تحت پیشپردازشهای مرسوم قرار گرفته و برای استفاده از تکنیکهای تصویرسازی زمانی آماده میشود. در ادامه دادهها به دو روش متفاوت تحت پردازش قرار خواهند گرفت. نتایج بهدست آمده از روش پردازش متداول نشان داد که این روشها در تصویرسازی لرزهای از ساختارهای پیچیده دارای ضعفهای متعددی هستند. بهمنظور برطرف کردن این مشکلات، از روش برانبارش سطح بازتاب مشترک استفاده شد. بااینحال بسیاری از روشهای نوین پردازش نیز مانند روش برانبارش سطح بازتاب مشترک علیرغم مزایای بسیار، در برخی مناطق با زمینشناسی پیچیده نیز در تصویرسازی لرزهای با مشکل روبهرو خواهند شد. به همین دلیل روند جستوجوی نشانگرهای جنبشی میدان موج و فرایند بهینهسازی آن در روش سطح بازتاب مشترک به گونهای تغییر کرد که بتوان در مناطق بسیار پیچیده نیز تصاویر لرزهای با کیفیت قابلقبول بهدست آورد. مقطعهای لرزهای و مقطعهای نشانگرهای جنبشی میدان موج بهدست آمده در این روش در مقایسه با مقطع لرزهای بهدست آمده در روش پردازش متداول نشان داد که میتوان در تصویرسازی لرزهای در ساختارهای پیچیده، به روشهای نوین اعتماد کرد.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Detection of subsurface installations and analysis of GPR pulse characteristicsآشکارسازی تأسیسات زیرسطحی و تحلیل مشخصههای تپ GPR11713339551FAمهدی محمدیویژهسازمان زمینشناسی و اکتشافات معدنی کشور، تهران، ایرانابوالقاسم کامکار روحانیدانشکده معدن، نفت و ژئوفیزیک، دانشگاه صنعتی شاهرود، ایرانJournal Article20161116Detection of subsurface installations and knowing their characteristics and accurate burial places are some challenges that the engineers involved always face with. Urban development, and thus, the necessity for the development of new and convenient installation networks and also the need for maintenance of old installation networks indicate the significance of accurate placement of the subsurface installations. In this regard, non-destructive methods have been considered to be used by many engineers in this field in order to save time and costs and also to reduce the disturbances due to unnecessary drillings in urban areas. In general, geophysical methods are among the most important non-destructive methods. One of shallow subsurface geophysical methods which has obtained some degree of success in this regard is the ground penetrating radar (GPR) method. In this research work, the GPR method is applied in several case studies having different subsurface conditions, and thus, the capability of the GPR method is investigated for accurate placement of subsurface installations in these case studies. Furthermore, we use a different procedure in data interpretation and analysis to analyze the characteristics of amplitude, phase and frequencies of GPR pulses in order to derive maximum information from GPR sections.
Three-dimensional (3D) data survey in one of case studies carried out in this research work had made it possible to process, interpret, and visualize subsurface structures in 3D form. Besides, resistivity surveys in the study area had provided us with the possibility of investigation and correlation of the GPR data results with the resistivity findings. Considering that a survey line had been carried out on the subsurface media containing high electrical conductivities, the unfavorable influence of the conductivities of survey media on the GPR measurements made along the survey line was indicated. Since the resistivities (or conductivities) of the subsurface media play an important role in the penetration of GPR waves, we can estimate the resistivities of the subsurface media by analyzing the frequency spectrum of the GPR pulses propagated in the media. Furthermore, the systematic noise affecting the GPR measurements is recognized by this frequency–spectrum analysis as this noise has frequencies different from those of signals. This study also considers GPR and resistivity investigations along a survey line crossing subsurface metallic and non-metallic installations, and thus, detailed investigations made on these subsurface targets. In this area, an estimation of the diameter of the non-metallic pipe installed for transferring of water in the subsurface has been made with a good accuracy from the acquired GPR time section along the survey line. Moreover, the analysis of amplitude and phase characteristics of the GPR pulses from some parts of the survey line has made it possible to distinguish or separate GPR responses of metallic installations from non-metallic ones. Finally, by considering GPR measurements along a survey line passing a subsurface non-metallic pipe, we investigated the effect of the fluid inside the pipe.
آشکارسازی تأسیسات زیرسطحی و آگاهی از مشخصات و محل دفن دقیق آنها یکی از چالشهایی است که مهندسان تأسیسات همواره با آن مواجهاند. گسترش شهرها و بهواسطه آن نیاز به توسعه شبکههای تأسیساتی و رفاهی جدید و همچنین نیاز به تعمیر و نگهداری از شبکههای قدیمیتر، تعیین محل دقیق آنها را بسیار پراهمیت ساخته است. در این راستا روشهای غیر مخرب بهمنظور صرفهجویی در وقت و هزینه و همچنین کاهش اختلالهای ناشی از حفاریهای بیمورد در مناطق شهری، مورد توجه بسیاری از مهندسان این رشته قرار گرفته است. یکی از روشهای ژئوفیزیک نزدیک سطح که در سالهای اخیر در این زمینه موفقیتهایی حاصل کرده، روش GPR است. در این تحقیق با در نظر گرفتن چند مورد متفاوت تحقیقاتی در شرایط متفاوت زیرسطحی قابلیت روش GPR در تعیین محل دقیق اینگونه ساختارها، مورد بررسی قرار گرفته است. بهعلاوه در این تحقیق با رویکردی متفاوت در تفسیر دادهها، با تحلیل مشخصههای دامنه، فاز و بسامد تپهای GPR سعی در دریافت حداکثر اطلاعات از مقطعهای GPR شده است.
برداشت سهبُعدی دادهها در یکی از سایتهای مورد بررسی، پردازش، تفسیر و نمایش سهبُعدی ساختارهای زیرسطحی را ممکن ساخته است. از طرفی برداشتهای مقاومت ویژه در این منطقه، بررسی و تطبیق نتایج حاصل از اندازهگیریهای GPR را فراهم کرده است. با در نظر گرفتن یک نیمرخ روی محیطهای با رسانندگی الکتریکی زیاد، تأثیر نامطلوب رسانندگی محیط برداشت در اندازهگیریهای GPR نشان داده شده است. ازآنجاکه مقاومت ویژه محیطهای زیرسطحی نقش مهمی در نفوذ امواج ایفا میکند، با تحلیل طیف بسامد تپهای GPR در این محیط، برآوردی از مقاومت ویژه محیط زیرسطحی بهدست آمده است. همچنین با تحلیل این طیف، نوفههای سامانمند که در اندازهگیریهای GPR تأثیرگذار و دارای بسامد متفاوتی از سیگنالها هستند، قابل شناساییاند. در ادامه با در نظر گرفتن یک نیمرخ روی تأسیسات فلزی و غیرفلزی دانشگاه صنعتی شاهرود، بررسیهای جزئیتر روی اینگونه هدفها صورت گرفته است. در این سایت برآورد قطر لوله غیرفلزی محتوی آب با استفاده از اطلاعات موجود در مقطع زمانی GPR، با دقت خوبی ممکن شده است. همچنین با تحلیل مشخصههای دامنه و فاز تپهای دریافتی در قسمتهایی از این نیمرخ، پاسخهای فلزی و غیرفلزی مورد بررسی، تفکیکپذیر است. در خاتمه این مقاله با در نظر گرفتن یک نیمرخ روی لولههای غیرفلزی، اثر شاره محتوی آنها در اندازهگیریهای GPR بررسی شده است.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116An investigation of factors affecting the inversion of a complex earthquake source using synthetic seismogramsبررسی عوامل موثر بر وارونسازی یک چشمه پیچیده زمینلرزهای با استفاده از لرزهنگاشتهای مصنوعی13414839552FAامیره سیالیموسسه ژئوفیزیک دانشگاه تهران، ایرانمحمدرضا حاتمیموسسه ژئوفیزیک دانشگاه تهران، ایرانظاهرحسین شمالیموسسه ژئوفیزیک دانشگاه تهران، ایران0000-0001-6254-7560Journal Article20161116In many earthquake source studies, the seismic source is assumed to have a double-couple (DC) source mechanism, matching shear motion on a planar fault. Observations of increasing quality and coverage, however, resolve departures from the DC model for many earthquakes and find some earthquakes, especially in volcanic and geothermal areas, that have strongly non-DC mechanisms. Burdick and Mellman (1976) used a time-domain inversion method to determine some of the complexity of the source time function. Several attempts have been made to explain the complexity of body waves from large earthquakes by using a multiple event model (Kikuchi and kanamori, 1982).
Deviation of earthquakes from the double-couple (DC) mechanism is an important, but delicate tool to study their source processes. It is believed that the earthquake focal mechanisms not only provide information regarding the stress field, but also give information about the rupture phenomena e.g. crack opening. Hence, the deviation from the double-couple mechanism is a matter of investigation. In fact, the study of the non-DC components makes opportunity to provide information about the earth processes. This information may be used to facilitate the operation of the geothermal energy and help predict volcanic activity. On the other side, the complex time history of energy release is a common attribute of large earthquake failures, as is the presence of non-uniform surface displacement along the outcrop of surface-breaking faults.
An interesting topic is the search for the connections between the non-DC events, multiple-double-couple events, segmentation of faults and their fractal properties. Critical papers emphasize difficulties in obtaining reliable non-DC components, e.g. due to noise, poor station coverage or incomplete structural models. Understanding non-DC earthquakes is important both for studying the process of faulting in detail and for identifying nonshear-faulting processes that apparently occur in some earthquakes.
To assess the non-double-couple component, a new method is suggested, i.e., a hierarchic grid search of the centroid position and time, during which the double-couple percentage (DC%) convergence is studied (Sokos and Zahradnik, 2008). This article uses the iterative deconvolution of multiple point sources, based on Kikuchi and Kanamori (1991), often used to study complexity of earthquakes .The method was modified for regional distances by Zahradnik <em>et al. </em>(2005). The modification involves the full Green’s functions (Bouchon 1981, 2003). Possibly complex events are represented by multiple point source models, which may represent their isolated asperities (Zahradnik and Sokos, 2008). ISOLA calculates the moment tensors (MT) by the least-square fitting of the complete waveforms in the time domain (Zahradnik and Sokos, 2005).
Due to the significance of an earthquake rupture process and investigating whether an event of high non-DC percentage can consist of several subevents with slight spatial and temporal intervals, this survey involves more accurate investigation using synthetic seismograms and their analysis. In this study, synthetic seismograms are produced using CPS software package and the seismograms of the main earthquake and its subevent were summed together and the resulting seismogram provided a basis for further studies. In fact, by modeling a source of high non-DC percentage (resulted from summing two 100% DC subevents), the ability of ISOLA to distinguish two subevents with different focal mechanisms is analyzed. Frequency band, crustal velocity model, moment ratio of earthquakes and the added noise value are determined as factors affecting successful retrieval of the two subevents.
<strong> </strong>
انحراف زمینلرزهها از سازوکار دو زوج نیرو بسیار مهم و یک ابزار حساس برای بررسی فرایندهای چشمه آنها است. این باور وجود دارد که سازوکار کانونی زلزلهها نه تنها اطلاعاتی در مورد میدان تنش فراهم میسازد بلکه در مورد پدیده گسیخت مانند بازشدگی ترکها نیز اطلاعاتی بهدست میدهد، به همین علت انحراف از سازوکار دو زوج نیرو، مورد بررسی قرار میگیرد. از طرف دیگر تغییرات زمانی پیچیده آزادسازی انرژی، مانند جابهجاییهای سطحی غیریکنواخت در امتداد برونزد گسلهای سطحی، یک خصوصیت گسیخت در زمینلرزههای بزرگ است.
به علت اهمیت نحوه گسیخت در یک زمینلرزه و بررسی اینکه یک رخداد با درصد غیر دو زوج نیروی زیاد میتواند شامل چند زیررخداد با فواصل مکانی و زمانی کم باشد، شناخت این موضوع نیازمند بررسی دقیقتر با استفاده از لرزهنگاشتهای مصنوعی و تحلیل آن است. در این تحقیق لرزهنگاشتهای مصنوعی با استفاده از مجموعه نرمافزاری CPS تولید و لرزهنگاشت حاصل از زمینلرزه اصلی و زیررخداد آن با همدیگر جمع و لرزهنگاشت حاصل مبنای کار واقع شده است. در واقع با مدلسازی یک چشمه با درصد غیر دو زوج نیروی زیاد (حاصل جمع دو زیررخداد 100 درصد دو زوج نیرو)، توانایی روش ISOLA برای تشخیص دو زیررخداد با سازوکارهای متفاوت، تحلیل میشود. نوار بسامدی، مدل سرعتی پوسته، نسبت گشتاور زمینلرزهها و میزان نوفه افزایشی بهمنزلة عوامل تاثیرگذار بر بازیابی موفیقت آمیز دو زیررخداد تعیین شد.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Comparison of high-resolution gridded monthly satellite and ground-based precipitation data over Iranمقایسه دادههای بارش ماهانه ماهوارهای و زمینی در شبکهای با تفکیک زیاد روی ایران14916039553FAپریسیما کتیرایی بروجردیدانشکده علوم و فنون دریایی، دانشگاه آزاد اسلامی واحد تهران شمال، ایرانJournal Article20161116In some countries such as Iran with large deserts and sparse rain gauge network satellite-based precipitation estimates have great potential for a wide range of applications. However, satellite-based precipitation estimates are not operational for decision making applications because of a lack of information regarding the associated uncertainties and reliability of these products. Obviously these data sets like others have error. To reduce the error, precipitation data users must have enough information about error characteristics over different parts of the world. In this analysis, satellite-based precipitation estimate data derived from the “Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks” (PERSIANN) is compared with ground-based data over Iran. The PERSIANN gridded precipitation data with 0.25°×0.25° latitude/longitude spatial and 3-hourly temporal resolution are used for evaluation of the satellite product. In this study, satellite-based precipitation data were accumulated to daily and then monthly totals for comparison with accumulated monthly gauge precipitation for the period 2003-2007 over Iran (25°-40° N, 44°-63°E). This rain gauge network included more than 2000 rain gauge stations over Iran. The arithmetic monthly mean of gauge precipitations calculated for every pixel. A comparison of the mean annual precipitation maps between these two networks over Iran shows that for the study period the spatial variation pattern of annual precipitation is reasonably accurate in PERSIANN, but it underestimates the precipitation amount over most parts of Iran. In the next step, the correlation coefficient and the scatter plot of the monthly precipitation for PERSIANN and gauge data were calculated for grid cells which included at least one gauge. Then in an attempt to better evaluate precipitation and reduce the effect of gauge uncertainties, the study was limited to pixels, each of which contained at least five rain gauges. The results for both of them (pixels which include one rain gauge and pixels which include five rain gauges) show that this satellite product underestimates the monthly precipitation. The correlation coefficient between satellite and gauge monthly precipitation for all pixels which include at least one gauge is equal to 0.3035 (99% significant). The correlation coefficient for pixels which include at least five gauges is 0.2598 (99% significant). For a comparison between these two data sets in different topography and climates, the time series of satellite and ground monthly precipitation are plotted for five cells which included at least five rain gauges. A pixel located in mountainous area of Zagros which most of annual precipitation falls in winter and summer is almost a dry season. Another pixel is in the most humid region in the coast of Caspian Sea which precipitation falls in almost all months of the year. One of the selected pixel located in the desert area of the eastern part of country. The two others located at North West and North East of country which have wet cold climate. The results show that PERSIANN underestimates the precipitation in Zagros, and it also has a very poor performance over the Caspian coast. On the other hand, this satellite product overestimates the precipitation over dry eastern area. بارش برآورد شده از دادههای ماهوارهای در زمینههای گوناگون از جمله بررسیهای هیدرولوژیکی کاربرد دارد. استفاده کنندگان این دادهها باید شناخت کافی از خصوصیات خطای این محصولات در مناطق گوناگون جهان داشته باشند. در این بررسی دادههای برآورد ماهوارهای بارش که به روش «برآورد بارش اطلاعات از راه دور با استفاده از نظریة شبکه عصبی» (PERSIANN) (Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks) بهدست آمده، با دادههای زمینی روی ایران مقایسه شده است. دادههای بارش ماهانه PERSIANN با تفکیک مکانی °25/0 طول و عرض جغرافیایی روی شبکهای با محدوده جغرافیایی بین °25 تا °40 شمالی و °44 تا °63 شرقی روی ایران طی سالهای 2003 تا 2007 با دادههای بارانسنجی با همین تفکیک مکانی شبکهبندی و مقایسه شد. این شبکه بارانسنجی متشکل از بیش از دو هزار ایستگاه اندازهگیری در سراسر ایران است. مقایسه نقشه میانگین بارش سالانه این دو شبکه روی ایران نشان میدهد که PERSIANN الگوی تغییرات مکانی بارش سالانه را بهخوبی میشناسد ولی مقدار آن را را در اغلب نقاط ایران کمبرآورد و در بعضی مناطق خشک کویری بیشبرآورد میکند. ضریب همبستگی و پراکنش بارش ماهانه PERSIANN و بارانسنجی برای همه سلولهای شبکه که حداقل شامل یک بارانسنج بودند و سپس برای سلولهایی که حداقل پنج بارانسنج داشتند بهدست آمد. هر دو مورد نشاندهنده کمبرآورد شدن بارش ماهانه این محصول ماهوارهای است. برای مقایسه این دو داده شبکهای در مناطق جغرافیایی و اقلیمهای متفاوت، نمودار سری زمانی بارش ماهانه ماهوارهای و زمینی پنج سلول که حداقل پنج بارانسنج داشتند رسم شد. نتایج نشان داد که PERSIANN بارش را در منطقه زاگرس کم و در کرانه دریای خزر بهشدت کم و همچنین سلول واقع در منطقه خشک شرق را زیاد برآورد میکند.
<strong> </strong>
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Improved particle swarm optimization for strong ground motion combinationارتقاء روش دسته ذرات بهینهیاب برای ترکیب مناسب نگاشتهای جنبش نیرومند زمین16116939554FAمحسن شهروزیدانشکده فنی و مهندسی دانشگاه خوارزمی، تهرانJournal Article20161116A step-by-step numerical solution to dynamic or nonlinear systems depends on the input time history record of the strong ground motion. Well-known seismic design codes recommend using a set of scaled records instead of only one. That is because a single input motion is not reliable enough to cover the seismic characteristics of all possible earthquakes for engineering design purposes. As a current state of practice, the scaling coefficients are taken similar and the result does not necessarily lead to a close compatibility of the resulted mean spectrum with the standard target. In this regard, the coefficients as floating-point numbers should be searched in the form of an optimization problem. Such a problem deals with a continuous type of search space with an infinite number of points. Thus complex optimizer engines are required to deal with it and seek a true global optimum set of scaling coefficients in order to suitably combine the selected ground motion records for further seismic design purposes. Particle Swarm Optimization, PSO, is one of the multi-agent meta-heuristic methods successfully applied to a variety of engineering design problems.
The present work formulates the problem of seeking optimal scaling factors using a utilized version of swarm intelligence as this method is best suited for continuous search spaces. In addition to inertial, cognitive and social terms of design point directions, an extra term is also utilized to improve the algorithm performance. In the first case, it is a new design vector with its elements selected one by one from other vectors in the current population of the swarm particles. The second modification is ordinating the velocity vector toward a randomly re-initiated particle position. In addition to the standard form of the particle swarm optimization, the aforementioned PSO versions are thus utilized, employed and further compared in the present work.
A variety of real-world ground motion records are picked from an available earthquake. The target spectrum is selected according to the Iranian Standard No.2800 for the soil type III and five percent damping considering the case of practical building structures. For each earthquake, the response spectra are generated for its longitude and transverse excitations and then combined using the Square Root of Sum of Squares (SRSS) method. The resulting SRSS acceleration response spectra are then averaged with the scaling factors with respect to the design target assigned by the developed optimization algorithms.
For the sake of true comparison, the same control parameters are chosen for three proposed versions of PSO, taken 2 for all the social, cognitive and extra congregation or combinatorial term except for the inertial coefficient which is taken 0.4. An error function is also defined to evaluate the compatibility of such weighted average spectrum with the target. The optimizer then samples various combinations of the scaling coefficients for the SRSS spectra and evaluated the error function for every such set. It seeks the optimal set among the continuous search space using the utilized artificial swarm intelligence for a pre-assigned number of particles and algorithm iterations. Consequently, the achieved results in the treated example show superiority of the optimized set over those in the code practice and considerable convergence/effectiveness improvement in the proposed methods over the standard particle swarm optimization.
ماهیت احتمالاتی زلزله مانع از آن میشود که بتوان به نگاشت یک زمینلرزه خاص برای تحلیل اعتماد و اکتفا کرد، لذا آییننامههای طرح لرزهای مقیاس کردن چندین نگاشت را برای کاهش حساسیت نتایج و امکان تصمیمگیری بهتر توصیه نمودهاند تا با مبنای تطبیق طیفی از وجاهت آییننامه طراحی برخوردار باشند. از سوی دیگر برای تحلیلهای کمّی دقیق در برآورد آسیبپذیری یا طراحی لرزهای تاریخچه زمانی جنبش نیرومند زمین موردنیاز است. پژوهش حاضر روش ذرات بهینهیاب را که طی سالهای اخیر در مسائل مهندسی توسعه یافته برای ترکیب بهینه شتابنگاشتهای زلزله فرمولبندی میکند و سپس دو شیوه ارتقا برای آن عرضه میکند که برپایه ترکیب جایگشتی متغیرهای حافظه موجود و مقداردهی احتمالاتی آنها استوار است. طی بحث نظری و بررسی مقایسهای صورت گرفته بین نتایج مزایای روش ابداعی حاضر و شیوه استاندارد در افزایش تطبیق طیفی شتابنگاشتها با طیف هدف آشکار میشود.
انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03367420161116Estimation of wind erosion emission from iron ore piles using CFD and Von Karman methodsبرآورد میزان بادبردگی از پشتههای سنگآهن با دو روش دینامک شارههای محاسباتی و ونکارمن17018039555FAخسرو اشرفیدانشکده محیط زیست دانشگاه تهران، ایرانمصطفی کلهردانشکده محیط زیست دانشگاه تهران، ایرانمجید شفیعپوردانشکده محیط زیست دانشگاه تهران، ایرانایوب ترکیانانستیتو آب و انرژی، دانشگاه صنعتی شریف، تهران، ایرانJournal Article20161116Open aggregate storage piles are used more and more in industrial sites. In industrial areas, emitted particulate matters from piles of row material can affect the quality of life of workers and employees and also the quality of the environment. Study of dispersion patterns and concentration of particulate matters over a landscape is important for the strategy of monitoring and controlling particulate matter. Within an industrial facility, dust emission may be generated by wind erosion of open aggregate storage piles and therefore, it pollutes the environment and wastes the row materials. Emission of particulate matters from surface of a pile depends on many parameters such as characteristics of wind (e.g. wind speed and wind direction), specifications of particles (e.g. particle diameter, density, shape, etc.) and erosion properties of surface. Therefore, for emission calculation of particulate matters from a pile and also for simulation of dispersion of emitted particles, it is necessary to simplify the physics of this phenomenon. Simplifications have been carried out based on governing equations and also applying the empirical relations obtained by field studies. Based on these theoretical and empirical investigations, a few methodologies are available for atmospheric wind erosion calculations from storage piles of row materials. The U.S. Environmental Protection Agency (EPA) method is one of the most famous approaches to this kind of calculation. It focuses on estimating the wind erosion and it cannot be used for dispersion pattern prediction. On the other hand, some models and methods have been developed to calculate the dispersion of pollutant in near and far distances from sources. One can combine the calculation of particulate matters emission with dispersion models in order to determine the particulate matters concentration at the environment.
In the present work, two methods including the U.S. EPA wind erosion estimating approach combined with Von-Karman's scheme for dust settlement and computational fluid dynamic (CFD) method using Fluent 6.3.2 Software are applied to predict particulate matters dispersion patterns from an iron ore pile. The Von-Karman's method is based on the length and time of the particulate matter settling. In the present work, the concentration of particulate maters in different distances from the source has been calculated using these parameters. In the CFD technique, the geometry of a pile is generated in Gambit Software using a structured mesh tool. The number of the generated mesh on the pile is 104,214. In this study, the flow condition is assumed to be incompressible, turbulent and steady state. Turbulence modeling is carried out based on two types of modeling namely and theories. Atmospheric wind profile is assumed to be in neutral conditions and defined by a user-defined function (UDF) tool from Fluent Software. The results from the two methods are compared with concentration of particulate matters measured based on 10 points. The maximum concentration position predicted by the CFD approach is more precise than that predicted by Von-Karman's method. Good quantitative and qualitative agreements are observed between the CFD predicted deposition and the measurement results. The determination coefficient for CFD and Von-Karman methods are 0.71 and 0.35, respectively. Also, Von-Karman method underestimates the concentration of particulate matters in all 10 measurement points.بررسی میزان و الگوی پخش ذرات در توپوگرافیهای گوناگون بهمنظور تعیین راهبردهای اندازهگیری و کنترل آلودگی ناشی از ذرات، دارای اهمیت بسیار زیادی است. در بسیاری از صنایع بهواسطه وجود دپوی مواد خام در نواحی گوناگون، شاهد پدیده بادبردگی هستیم که باعث آلودگی و اتلاف منابع مواد خام میشود. در بررسی پدیده بادبردگی، روشهای اندکی وجود دارد که یکی از بهترین این روشها، روش عرضه شده آژانس حفاظت از محیط زیست امریکا است. در این مرجع فقط به بررسی پتانسیل بادبردگی پشتههای گوناگون پرداخته شده و در مورد نحوه پخش و توزیع غلظت بحث نشده است. در این مقاله میزان بادبردگی از پشتههای سنگآهن صنایع فولادسازی با استفاده از روشهای عددی و تجربی مورد بررسی قرار گرفته است. از نرمافزار فلوئنت و روش تجربی ون کارمن در ترکیب با روش آژانس حفاظت از محیط زیست امریکا برای برآورد میزان بادبردگی در طول دوره یک ماهه از یک پشته هرمیشکل استفاده شده است. نتایج حاصل از هر دو روش با اندازهگیریهای صورت گرفته در محل در 10 نقطه اطراف پشته مقایسه شده است. بهمنظور برآورد عددی میزان بادبردگی، هندسه پشته و خصوصیات فیزیکی سنگآهن موجود در پشته شامل اندازه و توزیع دانهبندی ذرات در مدل دینامیکی شبیهسازی شده و نیمرخ باد در شرایط خنثی لایه مرزی جوّ بر آن اّعمال شده است. در حالت تجربی از روابط آژانس حفاظت از محیط زیست امریکا برای برآورد پتانسیل بادبردگی از پشته استفاده، و سپس نحوه پخش و توزیع آن با روش ونکارمن محاسبه شده است. نتایج هر دو روش در نقاط نمونهبرداری استخراج و با نتایج اندازهگیری مقایسه شده است. نتایج نشان میدهد که الگوی پخش کلی بهدست آمده از اندازهگیری و روشهای پیشگفته بهلحاظ کیفی مطابقت خوبی با یکدیگر دارند. روش دینامیک شارههای محاسباتی از نظر کمّی و کیفی تطابق بسیار خوبی با مقادیر اندازهگیری، بهخصوص در پایین دست جریان داشته است. ضریب تعیّن در حالت کلی برابر 71/0 برای مدل دینامیک شارههای محاسباتی و 35/0 برای روش ونکارمن بهدست آمده است. روش ونکارمن در همة نقاط، غلظت آلاینده PM<sub>10</sub> را کمتر از مقادیر اندازهگیری شده بهدست میدهد.