انجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Study of WRF model predictions for precipitable water and its relation with estimated precipitation by Tehran weather Radarبررسی پیشبینیهای مدل WRFبرای آب بارششو و ارتباط آن بابرآوردبارش به کمک دادههای رادار تهران33549FAامیر محمدیهامحمد حسین معماریان0000-0002-1899-8302مجید آزادیمحمد ریحانی پروریJournal Article20150419Precipitable water (PW) is an important meteorological quantity that cloud physics scientists have paid special attention to it. In fact, PW can be an estimation of the total column water vapor in the atmosphere which is the result of convergence of water vapor in the lower atmosphere. Awareness of the relationship between this quantity and rainfall amount and intensity is important in predicting the atmospheric conditions. In most part of the world, PW is measured with weather radars and satellites.
   The first purpose of this research is to find a relationship between precipitation and cloud precipitable water by utilizing Tehran weather radar data. The second purpose is checking the accuracy and skill of WRF model in forecasting and simulating the PW value and its changes. For fullfilling this purpose, particular rain gauges with the ability to record rainfall in short intervals (e.g. 1 hour intervals) are needed. Thus, three automatic rain gauges of Meteorological Organization (Pakdasht, Kooshk-e-Nosrat and Qom) were chosen. The Rainbow software was used to illustrate the Pakdasht, Kooshk-e-Nosrat and Qom Radar data. Desired products of radar data for this research were Surface Precipitation Intensity (SRI) and Vertical Integrated Liquid (VIL) (Actually VIL product is cloud precipitable water). For running the WRF model, three domains were considered including a 36 km horizontal range of parent domain and two nest domains with ranges of 12 and 4 km. It is worth noting that the relation between cloud precipitable water (CPW) and precipitation of radar measurements was investigated in the time period of 11/1/2010-11/4/2010.
   The following results were obtained by comparing the quantity of the cloud precipitable water (CPW) and radar surface rainfall intensity (SRI) in three considered stations. It was found that before the beginning of the rainfall, the amount of CPW of cloud always increases, and then with passing the amount of CPW through the 0.1 mm at that station, the radar shows SRI product on its display screen; this means that the rainfall is recorded by radar. After starting of the precipitation, the general trends of rainfall amounts would follow the trends in precipitable water values. However, the radar recorded a less SRI product at Kooshk-e-Nosrat and Qom stations compared with Pakdasht station for similar amounts of CPW. Meanwhile, the WRF model predictions for PW were compared with observations of this quantity by the weather Radar. Finally, it should be noted that the outcomes of model predictions for PW and Radar observations for CPW in the region under the cover of Tehran Weather Radar showed a very high accuracy (significant correlation at 5% significance level) in temporal-spatial forecasting and also in model simulation outputs for the changes of PW quantity. The results also revealed that whenever predicted PW by model having values more than 20 mm, radar observations displayed values more than 0.1 mm for the CPW.آب بارششو (PW) ، یکی از کمیتهای مهم هواشناختی است که همواره مورد توجه ویژه دانشمندان فیزیک ابر بوده است. در حقیقت آب بارششو برآوردی از ستون بخار آب در کل جوّ و حاصل همگرایی بخار آب در سطوح پایین جوّ است. آگاهی از ارتباط بین این کمیت با مقدار و شدت بارش، در پیشبینی وضعیت جوّ، از جمله امکان پیشبینی کوتاهمدت احتمال وقوع سیل دارای اهمیت زیادی است.
هدف از این پژوهش بررسی ارتباط بین دو کمیت بارش و آب بارششوی ابر به کمک دادههای رادار هواشناسی تهران و سپس بررسی دقت و مهارت مدل در پیشبینی و شبیهسازی میزان و تغییرات کمیت آب بارششو است. بدینمنظور ابتدا در سه ایستگاه بارانسنجی پاکدشت، کوشکنصرت و قم به بررسی ارتباط بین محتوای آب مایع و بارش سنجش شده با رادار در دوره زمانی 1 تا 4 نوامبر 2010 پرداخته شد.
نتایج نشان میدهد که اگر مقدار آب بارششوی ابر اندازهگیری شده با رادار از 1/0 میلیمتر بیشتر شود، بارش در آن منطقه آغاز میشود. سپس پیشبینیهای جوّی مدل عددی WRFاز کمیت آب بارششو با دیدبانیهای متناظر از رادار هواشناسی مقایسه شد. نتایج بررسی مهارت قابلقبول مدل عددی WRFرا در پیشبینی مقدار و تغییرات زمانی و مکانی این کمیت نشان میدهد که برای مقادیر بیشتر از 1/0 میلیمتر در آب بارششوی ابر سنجش شده با رادار، مدل نیز مقادیر بیشتر از 20 میلیمتر را در کل آب بارششوی جوّ نشان میدهد.https://www.ijgeophysics.ir/article_33549_9dc69c10a12392dfef7598b5bc086fd2.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Resolution improvement of velocity spectra using weighted semblance methodبهبود توان تفکیک طیف سرعت با استفاده از روش شباهت وزنی33550FAمرتضی رحمانیفاطمه خجستهامین روشندل کاهوشاÙرÙد0000-0002-2214-2558Journal Article20150419Velocity analysis is one of the most important steps of seismic reflection data processing which estimates the appropriate velocity for stacking and normal moveout correction. Taner and Koehler (1969) introduced the conventional semblance method for velocity analysis which is a normalized coherency measure. Semblance-based velocity analysis generates a velocity spectra for a common midpoint gather. Various parameters such as offset, signal to noise ratio, window length and etc. influence the resolution of semblance velocity spectra.
   For velocity analysis of a single reflector, if we consider the true velocity, then the hyperbolic reflection event aligns to the horizontal event in an analysis window and computed value of semblance is equal to one. If the operator selects an incorrect velocity, the hyperbolic event does not align to the horizontal event and the semblance value is reduced. If the selected velocity is close to the true value, then the computed semblance has a value close to one. Therefore, the resolution of the velocity spectra is reduced.
   Various methods have been introduced to improve the resolution of the semblance velocity spectra. Biondi and Kostov (1989) used the eigenstructure method to increase the resolution of velocity spectra. Sachi (1998) obtained the high resolution velocity spectra using the bootstrap algorithm. Roshandel et al (1387) improved the resolution of velocity spectra by multiplying the ratio of the two first eigen values of analysis window and the semblance value.
   In this study, we used the weighted semblance value introduced by Luo and Hale (2012) to obtain the high resolution velocity spectra. In this method, the weight functions are added to the semblance equation. The weight functions make the semblance value more sensitive to velocity and increase the resolution of the velocity spectra. Therefore, the weight functions must be calculated in a way that it is sensitive to the temporal change in normal moveout correction for every velocity at every offset.
   To illustrate the efficiency of the weighting function on the resolution of semblance spectra, weighted semblance and conventional semblance implemented on both synthetic and real CMP gathers. In the first synthetic example, we generated a CMP gather for a 37-layered earth model with linearly increased velocity. We computed its velocity spectra by both conventional and weighted semblance methods. Then, we added random noise to the synthetic CMP gather and re-calculated the velocity spectra. In these two synthetic examples, we have observed that the weighted semblance had a better resolution than conventional form. In the next step, we added multiple to CMP gather with near velocity to the primary event. Comparing the velocity spectra obtained from two mentioned methods showed that the conventional semblance cannot completely distinguish between the primary and multiple velocities. But they were easily separated from each other in velocity spectra of weighted semblance. Finally, both methods were applied to real data relating to the South West of Iran. The obtained result showed that extension of peak in the velocity spectra in the weighted semblance decreased. Thus, we can say the weighted semblance method has a better resolution than that of conventional methods.از مهمترین مراحل پردازش دادههای لرزهای بازتابی، تحلیل سرعت است که سرعت مناسب برای تصحیح برونراند نرمال و برانبارش را برای بازتابندهها، برآورد میکند. یکی از روشهای مرسوم تحلیل سرعت، محاسبه شباهت است. عوامل گوناگونی از جمله گسترش دوراُفت، نسبت سیگنال به نوفه، طول پنجره و مانند آن روی توان تفکیک طیف سرعت و پهنشدگی آن تاثیرگذار است. افزایش تفکیک در محاسبه شباهت سبب میشود که دقت برآورد سرعت تصحیح برونراند نرمال و برانبارش بهبود یابد و همچنین میتوان رویدادهای بازتابی را که از نظر سرعت و زمان دوطرفه به یکدیگر نزدیک هستند از یکدیگر تفکیک کرد. کیفیت طیف سرعت به حساسیت محاسبه شباهت به تغییرات سرعت، بستگی دارد. با محاسبه عبارات وزنی که به تغییرات سرعت حساستر هستند، میتوان کیفیت طیف سرعت شباهت را افزایش داد. پیادهسازی روش وزنی رابطه مستقیمی با بسط تیلور شباهت دارد. میتوان میزان شباهت را با استفاده از کمینه کردن تابع وزنی آن افزایش داد. در مقایسه، شباهت وزنی بهتر میتواند قلههای شباهت را برای تشخیص رویدادهای تداخلی متمایز کند. کارایی روش شباهت وزنی با اِعمال بر دادههای لرزهای مصنوعی و واقعی بررسی شده و نتایج آن با روش شباهت متداول مورد مقایسه قرار گرفته است. نتایج بهدست آمده بیانگر آن است که روش شباهت وزنی نسبت به روش شباهت متداول دارای توان تفکیک بیشتری است.https://www.ijgeophysics.ir/article_33550_77fd0f5b0a62adb13fab5379cc4c216b.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Random noise attenuation by rational-dilation wavelet transformتضعیف نوفه تصادفی با تبدیل موجک گسسته ضریب اتساع گویا33551FAمحمد ایرانیمهرمحمدعلی ریاحیJournal Article20150419The purpose of the geophysical activity is to recognize the underlying phenomena with a precision as high as possible. The record of noise in seismic data is inevitable. Different noises affect seismic data among which random noise is one of the most important ones, resulting from random oscillation of particles during the sampling. Random noise is present at all times and all frequencies; it decreases the seismic data quality. To achieve a trustworthy interpretation, geophysicists try to prepare data with the least random and coherent noise. There are many methods for analyzing seismic signals. In recent years, discrete wavelet transform has been introduced as a suitable choice in seismic signal processing. Bayram and Selesnick (2009) developed a family of wavelets transforms, which have a wide range of the Q-factor(in bandpass filters,Q-factor is defined as wavelet centre frequency/bandwidth). This research has tested the rational dilation wavelet transform (RDWT) with enhanced time-frequency discrimination. In this research the rational dilation wavelet transform was used for random noise attenuation from synthetic and real seismic data. The capability to apply differentQ-factors is the main advantage of the RDWT compared to the dyadic DWT. This transform overcame the limitation of conventional wavelet transform with a constant quality factor.
   This method is robust in the task of de-noising with fewer signal distortion effects because of its ability to choose an appropriate Q-factor and the degree of over-completeness.These WTs provide a rich range of redundancy andQ-factors. Moreover, the RDWT used in this study is based on rational (non-dyadic) dilation and attains over-completeness by increasing sampling in both time and frequency (Bayram and Selesnick, 2009).The RMS error is a measurement of the differences between values of a trace after an RDWT inverse transform without any de-noising and the input trace shown in Table 1 of the present article.
   The capability to apply differentQ-factors let the user choose the appropriate Q-factor and desired frequency and time resolution. The user can choose appropriate Q-factor setting some parameters (p, q, s, j),wherein âq/p'is the dilation factor and âsâis the sampling factor of the high-pass filter and the number of processing levels is determined by â jâ. These parameters must be chosen such that the time frequency representation of each level of decomposition will be enhanced. In an earlier work, this method was used with the purpose of ground roll attenuation (Iranimehr et al., 2013).
   Using the appropriate quality factor according to data type, the used wavelet will match well with the seismic signal and the signal obtains a large amount of wavelet coefficients but the undesirable random noise would attainsmall amount of wavelet coefficients.
   In the next step, noise was separated from signal using soft thresholding or other authorized thresholding methods. The wavelet transform was applied to synthetic data with different noise levels and a real marine data with high frequency content and the results of the random noise attenuation were compared with the dyadic wavelet transform. Table 2 of the present article compares the results of random noise attenuation by dyadic discrete wavelet transform and rational dilation wavelet transform by different type threshold.هدف از اجرای عملیات ژئوفیزیکی تشخیص پدیدههای زیر زمین با بیشترین دقت ممکن است. ثبتنوفهدربرداشتدادههای لرزهای اجتنابناپذیراست. نوفههای گوناگونیدادههای لرزهراتحتتاثیرقرارمیدهند،ازمهمترینآنهانوفهتصادفی است کهازحرکتنوسانی بینظمذراتدرهنگامدادهبرداریحاصل میشود. نوفهتصادفیدرهمهزمانها ودرهمهبسامدها حضورداردوکیفیتدادههای لرزهای راکاهشمیدهد. بهمنظور رسیدن به تفسیر قابلاعتماد، ژئوفیزیکدانان سعی بر رسیدن به دادهای با کمترین نوفه تصادفی و همدوس دارند. درسالهای اخیرتبدیلموجکگسستهدرحکمگزینهمناسبیدرامرپردازشسیگنالمعرفیشدهاست. دراینتحقیقازتبدیلجدید موجکگسسته ضریب اتساع گویا، باقابلیتبهبودتفکیکپذیریزمانی- بسامدیبرای تضعیف نوفهتصادفیاستفادهشده است.تبدیل موجک ضریب اتساع گویا دارای قابلیت بهدست آوردن محدودهای از ضریب کیفیتها است.در این تبدیل کاربر میتواند با تغییر چند متغیر به ضریب کیفیت مناسب برای پردازش داده دست یابد. این تبدیل با انتخاب ضریب کیفیت مناسب، با توجه به نوع داده، بهتر میتواند سیگنال را تجزیه کند و نسبت به تبدیل موجک گسسته معمولی با ضریب کیفیت ثابت، نتایج بهتری بهدست آورد.درمرحلهبعدبااستفادهازآستانهگیرینرمیا سایر آستانهگیریهای مجاز،امکان جداسازی نوفهازسیگنالمیسر میشود. تبدیل موجک جدید روی داده مصنوعی با سطح نوفههای تصادفی متفاوت و همچنین داده واقعی از نوع داده دریایی با محتوای بسامدی زیاد اِعمال شدو نتایج حاصل با نتایج تضعیف نوفه تصادفی با تبدیل موجک گسسته دوتایی مورد مقایسه قرار گرفت.https://www.ijgeophysics.ir/article_33551_cefcf59495dd6dafedd7809906e7807e.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Synoptic analysis of the onset of the earliest widespread winter precipitation in Iran(except the Caspian Sea coastal region)تحلیل همدیدی آغاز بارشهای فراگیر زودرس دوره سرد سال در ایرانauthor:33552FAارغوان رفیعاییبهلول علیجانیمحمدرضا یزدانیJournal Article20150419Different local and temporal distribution of rainfall across the Earth has been amongst the determining parameters in many civilizationsâ fate through history. Rainfall fluctuations and changes in climate phenomena are the monumental causes of these anomalies. Scientists have done a large number of researches to find the causes of the onset of widespread precipitation and its effects on agricultural crops, forecasting the onset of widespread rainfall, the teleconnection on widespread rain and the presentation of proper cultivation calendar for winter cultivation. Cold season brings about a great amount of widespread rainfall in many places in Iran. However, the basic consumption of water is during the warm season not only in civil but also in agricultural sectors. Therefore, being aware of the onset of widespread winter precipitation (OWWP) and its year to year fluctuations is a pressing issue. To examine and analyze the pattern of atmospheric circulation for the earliest OWWP, a 33-year period of rainfall data from 50 synoptic stations and one rain gauge station within four areas of Iran were analyzed. To find the early onset of the widespread rainfall, the average daily precipitations of all stations through each area were calculated. Furthermore, OWWR has two particular features per year. Firstly, the day having the higher standard deviation than the average rainfall of all the areasâ stations would be the onset day. Secondly, there should be at least two consecutive days having precipitation in more than 50% of all stations in the given area. As a result, there was a specific onset of rainfall for each area during 33 years. Moreover, if the OWWR in one year had started before the rainfall onset of that year, it would have been the year with the early onset. Consequently, the geopotential and vorticity maps of 500 and 1000 hPa showed the presence of a significant mridional component of atmospheric circulation in the earliest years with abundance of blocking highs and cut-off lows especially in the East Mediterranean and Iran. In addition, the anomalous positive relative vorticity in the earliest years lasted until the end of winter but the negative anomaly of relative vorticity showed the lowest amount compared with the other years. Finally, the result of Tukeyâs test (HSD) showed a meaningful correlation (at 5% confidence level) between the OWWP and the total amount of rainfall in those above mentioned earliest years. In conclusion, 1974 in the northwest part, 1976 in the west, 1991 in the south and southwest and 1982 in the center and east parts of Iran presented the earliest widespread winter rainfall during those 33 years. In those mentioned years, the polar vortex had the highest positive anomaly and the high pressure of Saudi Arabia had its most optimal conditionon the Arabian Sea which caused a wet year for Iran. Also in the mentioned years, heavy rain and flood have been seen in the semi arid areas of the East Mediterranean, Saudi Arabia and Sienna desert.توزیع ناهمگون زمانی و مکانی نزولات جوّی در سطح کره زمین از عاملهای تعیینکننده سرنوشت بسیاری از تمدنها در طول تاریخ بوده است. در ایران، فصل سرد سال در برگیرندة عمدة بارشهای فراگیر در اکثر نقاط است، درحالیکه عمدة مصارف آب شیرین حاصل از بارندگی در بخش خانگی و یا در بخش کشاورزی در فصل گرم سال است. ازاینرو، آگاهی از زمان آغاز بارشهای فراگیر دورة سرد سال و نوسان آن در سالهای متفاوت، بهلحاظ در برداشتن بیشترین بارشها اهمیت ویژهای دارد. در این تحقیق با هدف بررسی و تحلیل همدیدی الگوها و بی هنجاریهای گردش جوّ در سالهای همراه با آغاز زودرس بارشهای فراگیر در 4 منطقه بارشی در ایران، 33 سال آمار بارش روزانه دوره سرد سال (سپتامبر - مارس) در 50 ایستگاه همدیدی و 1 ایستگاه بارانسنجی، مورد بررسی قرار گرفت. نقشههای فشار و تاوایی تراز 500 و 1000 هکتوپاسکال برای اولین بارش فراگیر در سالهای زودرس، الگوی گردش جوّ درهمریخته و دوشاخه و نصفالنهاری را نشان داد. درنهایت آزمون آماری توکی همبستگی معنیداری را در سطح 5 درصد بین آغاز زودرس بارشهای دوره سرد سال در هر ناحیه و کل مقادیر بارش دریافتی آن دوره نشان داد.https://www.ijgeophysics.ir/article_33552_2f775cb9cdd67bcd731e35cceeba3f9d.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Atmospheric stability analysis and its correlation with the concentration of air pollutants: A case study of a critical air pollution episode in Tehranتحلیل پایداری جوّی و همبستگی آن با غلظت آلایندههای هوابررسی موردی بحران آلودگی هوای تهران33553FAخسرو اشرفیداÙشگا٠تÙراÙعلی احمدی ارکمیتÙراÙJournal Article20150419In the present survey, it has been tried to demonstrate a significant correlation between air pollutant concentrations and meteorological parameters, by the study of the atmospheric conditions in an air pollution episode. Carbon monoxide and nitrogen dioxide were selected as air pollutants. Temperature, wind speed and direction and the dimensionless Richardson number were considered as the meteorological parameters in this study. To this aim, a period of eleven days from 30 November to 10 December 2012 was considered. In this period, a severe temperature inversion and consequently a long-term stable situation was seen in the Tehran weather. The governing stability classes were studied using the Turner algorithm and vertical temperature gradient scheme. For the mentioned period, the sun altitude, wind speed at the height of 10 meters, total cloud cover and ceiling were extracted from the U.S. Naval Observatory web page (http://aa.usno.navy.mil/data/docs/AltAz.php) to be used in the Turner algorithm. The stability analysis showed that for this time interval, stable conditions prevailed in 89% of cases. The sounding data was collected from the Mehrabad weather station. By the sounding data, the daily maximum mixing depth was calculated and the vertical trend of the temperature was plotted in each day during this period. The morning vertical temperature gradients showed a strong temperature inversion in these eleven days. Also, the maximum mixing depth decreased to the minimum value of 600 meters on the 6th day of the time interval. The maximum mixing height and visibility graphs also confirmed that the atmospheric mixing decreased and consequently more pollutants were trapped in the middle of the period. Once confirming critical conditions by comparing the stability classes, maximum mixing depth and the visibility, the correlation between atmospheric parameters and pollutants concentration was obtained by multiple regression method. Among the meteorological parameters that were considered, temperature, wind speed, wind direction and the Richardson number had a greater correlation with the pollutant concentration. The bulk Richardson number was calculated from the sounding data and used in the regression. The CO and NO2 concentrations were gathered from the pollutant concentration measuring station which is located on the region-10 municipality building, 2500 meters far from the Mehrabad weather station. The coefficients of correlation between the normalized CO concentrations and temperature, West-East component of wind, South-North component of wind and the Richardson number were 0.707, 0.078, 0.028 and -0.019, respectively. And for NO2 concentration, the correlation coefficients were 0.353, 0.016, 0.015 and -0.019, respectively. It could be observed that the coefficients of temperature and surface wind for the CO concentration were almost twice the corresponding coefficients for the NO2 concentration. The reason is that the NO2 concentration has been influenced by chemical and photochemical reactions and the mixing depth, while the CO concentration is only affected mainly by the rate of transportation and the mixing depth. In fact, the CO concentration had a stronger functionality to temperature and wind speed than the NO2 concentration. In the most cases, the Richardson number had a positive value and by increasing its value (in the early morning hours or night hours), the pollutant concentration were reduced. Therefore, as expected, its coefficient were negative in both regressions.در تحقیق حاضر سعی شده است تا با بررسی شرایط جوّی در بحرانهای آلودگی هوا، رابطه معنیداری بین غلظت آلایندههای کربن مونوکسید و ازت دیاکسید و پارامترهای هواشناسی مانند دما، سرعت و جهت باد و عدد بیبُعد ریچاردسون حاصل شود. برای این مهم، یک دوره یازده روزه، از 30 نوامبر تا 10 دسامبر 2012 در نظر گرفته شد که در این دوره یک وارونگی دمای شدید و به تبع آن شرایط پایداری طولانیمدت برای هوای شهر تهران مشاهده شده بود. ردههای پایداری حاکم با استفاده از روش ترنر (Turner)و روش گرادیان دمایی مورد بررسی قرار گرفت. پس از تایید بحرانی بودن شرایط حاکم از روی مقایسه ردههای پایداری، بیشینه عمق آمیختگی و میدان دید، ارتباط بین پارامترهای جوّی و غلظت آلایندهها با بهکارگیری روش وایازش (رگرسیون) چندمتغیره بهدست آمد. نتایج نشان میدهد که در دوره مورد بررسی، شرایط پایدار در 89 درصد مواقع غالب بوده و کمترین مقدار بیشینه عمق آمیختگی در روز ششم و به میزان 600 متر رسیده است. همچنین ضرایب همبستگی غلظت بهنجار شده آلاینده کربن مونوکسید به پارامترهای دما، سرعت باد در راستای غرب به شرق، سرعت باد در راستای جنوب به شمال و عدد ریچاردسون بهترتیب برابر 707/0 و 078/0 و 028/0 و 019/0- است و این همبستگیها برای آلاینده ازت دیاکسید بهترتیب برابر 353/0 و 016/0 و 015/0 و 019/0- بهدست آمد.https://www.ijgeophysics.ir/article_33553_17a6bfc3aa330fcc0b32da96c283caaa.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Investigation of tectonic lineaments by using magnetic data in Poshtbadam area, Yazd provinceبررسی خطوارههای زمینساختی در ناحیه پشتبادام با استفاده از روش مغناطیسسنجی33554FAمسلم جهانتیغغلامرضا نوروزیفرشاد ژولیدهسرJournal Article20150419Magnetic method is an important passive surface geophysical method which uses the Earthâs natural magnetic fields to investigate the susceptibility structure of the subsurface material. The geomagnetic field measurements can be used to determine the structure of the earth, because the rocks often contain magnetic minerals. The interpretation of data of this nature, in some cases, will determine geological characteristics that would help contribute to the success of mining or oil exploration. In this study, for the review of more details and achieving an idea about accurate location, dip, depth and mineralization spread, geomagnetic surveys were applied and then data processing method were used. Preparing a suitable map which contains the informative data without interference from noise is the main purpose of the data processing. To reach this aim, the acquired data should be corrected for different non-geological effects.
These are: 1. Diurnal correction to remove time varying parts of the magnetic field. 2. IGRF correction to remove the field of earth core and upper mantle response
Many interpreting methods, that estimate the depth, location and the shape of a potential source, are based on using the gradients of potential fields. Derivatives are high pass filters. They intrinsically amplify any noise and shallow anomalies present in the data. Therefore, using high order derivatives would be less common. Tilt angle has many interesting properties. For example, due to the nature of arctan trigonometric function, all tilt amplitude are restricted to values between -90° and +90° (Salem et al. 2007). Another property of tilt angle is that the value of tilt angle above the edges of the contacts is 0°, and it is equal to 45° when h = zc and -45° when h = zc. This suggests that contours of magnetic tilt angle can identify the locations of contact-like structures
The studied area is located in Yazd province and it is part of micro continent Iran central. The intrusive rock in this area is genais and the sedimentary rock is schist and the host rock of studied area is limestone. Igneous units of the trust are located on limestone. Also outcrop of hematite  has been observed in this area. Designing the field operation grid were carried out based on geology and primary studies of region and therefore 2400 points were along 60 profiles. After diurnal correction, total magnetic intensity has been corrected by using IGRF 2000 so that components of Earth's magnetic fields have been eliminated. Also reduction to pole operation was done. In this paper, first derivative, second derivative, analytical signal tilt angle and Euler deconvolution methods were used to study fractured structures. The edges of these faults were detected by using analytical signal, first derivatives, second derivatives and the tilt angle method
Then, we use the standard Euler deconvolution method. Reid et al (1990) used this method as Grade analysis. Simple geometric models for magnetic field source (Balkely, 1995) and determination of structural index with earlier data are disadvantages of the Euler method. However, Thompson (1982) and Reid et al (1990) determined the optimum structural index for different structures; nevertheless the Euler method estimated quite well the location and depth of magnetic sources (Aboud et al., 2005)
The results of the first vertical derivative showed three major faults in this area strikes of which were east-west in north, northeastern-southwestern in middle and northwestern-southeastern in south of the area and the other interpretation methods have an acceptable degree of reliability. The faults specified by using these methods that are offsprings of mineralization have good compatibility with geological information. Also, the depths of the faults in the area were estimated by the Euler method. The estimated depths of the top of major faults in the area were similar at 5 m and 15 m.بهمنظور پیجویی خطوارههای زمینساختی در منطقه پشتبادام و نقش آن در کانیزایی این محدوده دادههای مغناطیسی در این گستره مورد پردازش و تفسیر قرار گرفته است. بهاینمنظور پس از مشاهدات صحرایی، شبکه برداشت10×20 در 60 نیمرُخ و در 2400 ایستگاه طراحی شد. در همین راستا با استفاده از فیلترهای مشتق قائم مرتبه اول، مشتق قائم مرتبه دوم، سیگنال تحلیلی و زاویه کجی امتداد گسلها بهطور دقیق مشخص و همچنین با استفاده از روش اویلر عمق گسلهای محدوده مورد بررسی برآورد شد. در این محدوده سه گسل با راستاهای، گسل1، شمال غرب-جنوب شرق، گسل 2، شرقی-غربی و گسل 3، شمال شرق-جنوب غرب مشاهده میشود. در واقع به نوعی میتوان عامل اصلی کانیسازی را زون خُرد شده ایجاد شده از گسلهای نامبرده معرفی کرد. پس از پردازشهای لازم روی دادهها، روش اویلر برای برآورد عمق، مورد استفاده قرار گرفت. اندیس ساختاری و اندازه پنجره با بررسیهای زمینشناسی 9/0 و 10×10 تعیین شد که درنهایت عمق سطح فوقانی گسلها 5 تا 15 متر برآورد شد که با نتایج حاصل از تفسیر کیفی تطابق خوبی دارد.https://www.ijgeophysics.ir/article_33554_f3972f5473ee7513f2b2fd52e2ee7893.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Site characterization of the source microtremors using the H/V methodاثر چشمه میکروترمورها در تعیین مشخصات ساختگاهی با استفاده از روش H/V33555FAرضا حیدریمحمدرضا قائمقامیانظاهر حسین شمالی0000-0001-6254-7560Journal Article20150419In this study, the local site amplification characteristics are investigated at 10 accelerometer stations organized by Tehran Disaster Mitigation and Management Organization (TDMMO) in Tehran city. The microtremor data were selected from 40 days contious recording at each station. The site amplification functions are calculated using horizontal to vertial spectral ratio (H/V) in the frequecny domain. The results revealed large sensitivy to the microteremor sources. When compared the results with previous studies using earthquake data, the computed amplification functions showed the peaks not related to the site or different amplification values at predominant peak frequencies at the sites. To find the reasons, the amplification function at the sites was calculated using microteremor data recorded at weekdays and holidays. we calculated Meanwhile, the results of calculated amplification functions for holidays showed a good agreement with those of previous studies that  meansthe industrial noise largly affect the results, especially in the south of Tehran.Furthermore, It was found that the industrial noise mostly affected the frequency range of 1 to 1.5 Hz. This also led to increase the amplification coefficient in an order of 1.5 to 2.2. This emphesized on the fact that the microtremor data recorded in the mega-cities like Tehran should be carefully analyzed  to prevent misleading results. Furthermore, it is suggested to use microtremor data recored in the holidays to avoid such heavy industrial noise contamination.در این مقاله تابع بزرگنمایی ساختگاه در 10 ایستگاه شتابنگاری متعلق به شبکه شتابنگاری سازمان پیشگیری و مدیریت بحران شهر تهران در گسترده شهر تهران تعیین شده است. بدین منظور از ثبت نوفههای پیوسته در 40 روز متوالی و تحلیل آنها به روش نسبت طیفی مولفههای افقی به قائم (H/V)، استفاده شده است. نتایج بهدست آمده نشاندهنده حساسیت تابعهای بزرگنمایی ساختگاه به میزان و منشأ نوفه صنعتی بوده است. این تابعها در برخی از محدودههای بسامدی بهشدت تحتتأثیر نوفههای صنعتی حاصل از کارخانهها بوده و هر دو پارامتر بسامد غالب و ضریب بزرگنمایی ساختگاه را در آن محدودههای بسامدی متأثر کرده است. بررسی تابعهای بزرگنمایی نشان داد که در شهر تهران و بهویژه در مناطق جنوبی آن کارخانههای صنعتی بسیاری وجود دارد که دامنههای نسبتاً قوی در محدوده بسامدی 1 تا 5/1 هرتز را بهوجود میآورند و باعث افزایش کاذب ضرایب بزرگنمایی ساختگاه به میزان 5/1 تا 2 برابر در این محدوده بسامدی شدهاند. این در حالی است که بسامد غالب و ضرایب بزرگنمایی بهدست آمده در شهر تهران در غیاب این نوفههای صنعتی مقادیری متفاوت و منطبق با نتایج تحقیقات قبلی دارند. بنابراین، در روش نسبت طیفی مولفههای افقی به قائم در تعیین اثرهای ساختگاه برای شهرهای صنعتی (با نوفههای صنعتی شدید) مشابه تهران بایستی نقش نوفههای صنعتی را بهمنزلة عاملی مهم و تأثیرگذار مشخص نموده و در نتایج تحلیل مدنظر قرار داد. این بدان دلیل است که دامنه قلههای موجود در بسامدهای مربوط به نوفههای صنعتی گاهی ممکن است درحکم بسامد غالب ساختگاه تعبیر و منجر به ارائه نتایج گمراه کننده در تفسیر اثرهای ساختگاهی در شهر تهران گردد.https://www.ijgeophysics.ir/article_33555_df73c84150382846a2ab31caec444c30.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Determination of ambient seismic noise directionality in the Central-Alborz/Iran using cross-correlation functionsتعیین جهتیافتگی نوفه لرزهای محیطی در منطقه البرز مرکزی/ایران با استفاده از تابعهای همبستگی33556FAتقی شیرزادظاهرحسین شمالی0000-0001-6254-7560Journal Article20150419Recent progress in seismology has demonstrated that empirical Greenâs functions (EGFs) of inter-station distances can be extracted using cross correlation of ambient seismic noise recorded in the similar time at two stations (Weaver and Lobkis, 2002; Shapiro and Compillo, 2004; Wapenaar, 2004). Consequently, this method provides a great set of data even in low seismicity regions to apply in the tomographic studies. Thus, the resulted tomographic images using the ambient seismic noise method (hereafter ANT) can show interior earth structures with a higher resolution compared to classical tomography methods (Shapiro et al, 2005; Lin et al., 2007; Shirzad et al, 2013).
Diffused signals are the main assumption in the ANT method (Snider, 2004). Ambient seismic noise sources generate a coherent and transient noise wavefield with random amplitude and phase in a medium (Van-Tighelen, 2003; Gorin et al., 2006). Reconstruction of the propagating path information using the amplitude of the recorded noise wavefield is impossible, but coherent information provided by propagating path can be extracted using cross correlation of long time ambient seismic noise recorded (Weaver and Lobkis, 2004; Gorin et al., 2006). This coherent information is called elastic response of medium or empirical Greenâs functions (Shapiro and Compillo, 2004; Roux et al., 2005; Sabra et al., 2005).
   Generally, the ambient seismic noise recorded for each station is composed of surface waves (Rayleigh and Love) with random amplitude and phase (Aki and Richards, 1980). Cross correlation function of these data will be symmetric if the ambient seismic noise wavefields generated by random sources are distributed uniformly (Snider, 2004). Earth structures can be studied using travel-time of extracted EGFs such as Rayleigh wave fundamental mode (Shapiro et al., 2005). Some studies (e.g. Stehly et al., 2006; Pedersen et al., 2007) indicate that the inhomogeneous distribution of the signal energy in various azimuths, which results in directionality of ambient seismic noise, produces deviation in tomography results and causes incorrect interpretations. Consequently, optimization of extracted tomographic maps based on the ANT method needs comprehensive knowledge of spatial and seasonal distribution of the noise wavefield in study areas (Stehly et al., 2006; Pedersen et al, 2007).
   Gutenberg (1936) suggested that the sources of primary and secondary oceanic microseisms observed throughout the Europe are located in the northeastern Atlantic Ocean. Primary and secondary microseisms dominate the noise wavefield in certain frequency ranges. The interaction between the swell and the sea bottom generates the primary microseisms which are dominated by periods of 12â25 s. Also, interfering water wavefield components travelling in opposite directions generate the secondary microseisms which are dominated by periods of 5-10 s (Gutenberg, 1936).
   In this study, we analyzed three-component recordings of continuous data from 30 stations in the Central Alborz region depicted in Figure 1. The Alborz Mountain range in the southern margin of the Caspian Sea is a part of the AlpineâHimalayan orogenic belt. The Alborz Mountain range resulted from a stress state derived from the horizontal compressive forces of the Central Iran Plateau has been induced by the collision of the Arabian plateau and the Asian continent (Berberian and King, 1981; Zanchi et al., 2006). The dataset used in this study consisted of 10 digital accelerometers with CMG-5TD sensors operated by the Tehran Disaster Mitigation and Management Organization (TDMMO), 18 digital narrow-band seismometers with SS1 seismometer sensors (corner frequency â¥1 Hz) operated by the Iranian Seismological Center (IRSC) at the University of Tehran, and two digital broadband instruments with a CMG-3T sensor operated by the International Institute of Earthquake Engineering and Seismology (IIEES). For the TDMMO acceleration network, the IRSC and the IIEES seismic networks continuous data from 2010 were analyzed.
   In the case of azimuthal distribution of the ambient noise, normalized amplitude of the cross-correlations versus azimuth (rose-diagram) constrained the direction to the sources of the ambient seismic noise, based on all available station-pairs. The average fractions (the number of Love/Rayleigh path with a SNR>10 in a given 20° azimuthal bin was normalized to the total number of Love/Rayleigh paths in that given bin) of the Love and Rayleigh yearly empirical Greenâs functions with a SNR>10 were in the orders of 0.78 and 0.73, respectively, at the period band of 1â10 s. Our final results indicated that the average fractions per cent of Love and Rayleigh paths with SNR>10 were above 68% and 64% on a yearly scale, and never decreased to 45% and 50% on a monthly scale at the period band of 1-10 s, respectively.امروزه با بهکارگیری نوفه لرزهای محیطی، مجموعه بزرگی از دادهها در راستای تعیین تصاویر توموگرافی با کیفیت بالا در مقیاسهای محلی و منطقهای، حتی در مناطقی با لرزهخیزی کم، فراهم شده است. پیش از استفاده از تابعهای گرین تجربی بینایستگاهی حاصل از همبستگی نوفههای لرزهای لازم است که توزیع فضایی و زمانی نوفه لرزهای بررسی شود. در این تحقیق جهتیافتگی نوفه لرزهای در البرز مرکزی/ایران با استفاده از دادههای پیوسته سه مولفهای دائم ثبت شده در30 ایستگاه لرزهنگاری با حسگرهای متفاوت عرضه میشود. همچنین برای تعیین تابعهای گرین تجربی بینایستگاهی فقط از همبستهساختن نوفههای لرزهای محیطی ثبت شده در ایستگاههایی با حسگر مشابه استفاده شده است. نتایج نشان میدهد که بیش از 68% و 64% انرژی سیگنالهای امواج بهترتیب لاو و ریلی سال 2010 در همه سمت (آزیموت)ها در بازه تناوبی 1 تا 10 ثانیه توزیع میشود. بررسی توزیع ماهیانه سیگنالها نشان میدهد که حداقل 45% و 50% انرژی سیگنالهای بهترتیب لاو و ریلی در همه ماههای سال 2010 در همه سمتها در بازه تناوبی 1 تا 10 ثانیه توزیع شدهاند.https://www.ijgeophysics.ir/article_33556_7b979cd61673c4453ca8fdeacb38f089.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Porosity estimation of pore volume in one of the Persian Gulfâs hydrocarbon fields use multi-attribute analysesبرآورد تخلخل مخزن از راه تحلیل چندنشانگری در یکی از میادین نفتی خلیج فارس33557FAسالار زحمتکشعبدالرحیم جواهریانJournal Article20150419Integration of 3D seismic data with petrophysical measurements gives a better vision to reservoir characterization. The integration of well-logs and seismic data has been a consistent aim of geoscientists which become increasingly important and successful. In recent years, because of the shift from exploration to development of existing fields with a large number of wells penetrating them, improving reservoir study has been the most important pre-drilling activity. One type of integration is forward modelling of synthetic seismic data from the logs. A second type of integration is inverse modelling of the logs from the seismic data. It is called seismic inversionSeismic inversion, in geophysics, is the process of transforming seismic reflection data into a quantitative rock-property description of the reservoir. Another method is to estimate the log properties by seismic attributes. In this study, linear multi-attribute transform, and non-linear multi-attribute transform were used for predicting porosity in one of the Iranian hydrocarbon fields. The analysis data consisted of the target log (in this study, the porosity log) from wells tied with 3D seismic volume. From the 3D seismic volume, a series of sample-based attributes was calculated. The objective was to derive a multi-attribute transform, which was a linear or nonlinear transform between a subset of attributes and target log values. The selected subset was determined by a process of forward stepwise regression, which derived increasingly larger subsets of attributes. In the linear mode, the transform consisted of a series of weights derivedby least-squares minimization. These weights are coefficients of the selected attributes in a linear multi-attribute transform. In the nonlinear mode, a neural network was trained using the selected attributes as the input.Two methods of neural network used in this study include probabilistic neural network and multi-layer feed-forward network. The basic idea behind the general regression probabilistic neural network is to use a set of one or more measured values, called independent variables, to predict the value of a single dependent variable. The multi-layer feed-forward network method consists of a set of neurons, arranged into two or more layers. There is always an input layer and an output layer, each containing at least one neuron. Between them, there are one or more âhiddenâ layers. The neurons are connected in the following fashion: inputs to neurons in each layer come from outputs of the previous layer, and outputs from these neurons are passed to neurons in the next layer, and each connection represents a weight.
   To estimate the reliability of the derived multi-attribute transform, cross-validation was used. In this process, each well was systematically removed from the training set, and the transform was rederived from the remaining wells. Then, the prediction error was calculated for the hidden well. The validation error, which is the average error for all hidden wells, was used as a measure of the likely prediction error when the transform was applied to the seismic volume. There was a continuous improvement in predictive power as it was progressed from a single-attribute regression to a linear multi-attribute prediction. This improvement was evident not only in the training data but, more importantly, in the validation data. In addition, the neural network did not show a significant increase in resolution over that from the linear regression. As a conclusion, the best result of porosity estimation in this field was provided by the linear multi-attribute transform.همراهسازی دادههای سهبُعدی لرزهای با اندازهگیریهای پتروفیزیکی حاصل از چاهها راه مناسبی برای دستیابی به توصیف مناسب از مخزن است. روش مرسوم برای همراهسازی دادههای لرزهای و دادههای چاه، وارونسازی دادههای لرزهای است. روش دیگر، برآورد خواص مخزنی از قبیل تخلخل با استفاده از نشانگرهای لرزهای است. در این تحقیق، تبدیل تکنشانگری، چندنشانگری خطی و تبدیل چندنشانگری غیرخطی برای برآورد تخلخل در یکی از میادین نفتی ایران در خلیج فارس مورد استفاده قرار گرفته است. این روش بر پایه تعریف رابطه ریاضی بین چند نشانگر لرزهای و خاصیت مخزنی موردنظر، مانند تخلخل است. نشانگرهایی که با خاصیت مخزنی موردنظر (تخلخل) رابطه خطی قابلقبولی دارند، در کنار یکدیگر مورد استفاده قرار میگیرند و بهبود قابلتوجهی در انطباق بین تخلخل پیشبینی شده از این روش، با تخلخل واقعی حاصل میشود. در انتها مدل سهبُعدی مخزن تهیه شده است و فضای متخلخل مخزن برای مکعبهای تخلخل که با استفاده از تبدیل تکنشانگری و چندنشانگری خطی و غیرخطی بهدست آمدهاند، محاسبه و مقایسه شده است. تبدیل چندنشانگری خطی با بیشترین میزان همبستگی در برآورد لاگ، صحیحترین فضای متخلخل را پیشبینی میکند. همبستگی اعتبارسنجی در برآورد تخلخل، تاثیر بسزایی در میزان فضای متخلخل مخزن دارد.https://www.ijgeophysics.ir/article_33557_af25cfccb0934c9f4386176c531295f2.pdfانجمن ملی ژئوفیزیک ایرانمجله ژئوفیزیک ایران2008-03368320140923Improving Ground Penetrating Radar (GPR) forward modeling approach using the numerical finite difference methodبهبود مدلسازی پیشرو دادههای رادار نفوذی به زمین (GPR)به روش عددی تفاضل متناهی33558FAرضا احمدینادر فتحیانپورغلامحسین نوروزیJournal Article20150419The behaviour of Ground Penetrating Radar (GPR) electromagnetic field can be simulated using Maxwellâs equations and associated boundary conditions. So far a number of numerical methods for modeling GPR data have been proposed including the popular Time Domain Finite Difference (TDFD) technique. The popularity of TDFD is mainly due to being relatively simple to implement, its high flexibility and capability to simulate complex subsurface geology. Also, the TDFD approach is not only conceptually accurate for complex geological models but also enables us to design realistic antenna and to study physical electromagnetic phenomena such as dispersion in electrical properties. Despite having these advantages, the finite difference method has pitfalls such as becoming very time consuming in simulating the most common media especially with high dielectric permittivity causing the forward modeling process to become very time consuming even by modern high-speed computers.
Synthetic GPR responses are useful for predicting expected GPR data over known geometries such as horizontal cylinders and prisms. The GPR data and subsurface lithological and hydrogeological properties are related through a numerical forward modeling engine.Â
Therefore, the GPR forward modeling engine can transform the subsurface electrical properties into expected GPR responses which in turn can be used for optimizing data acquisition procedures on pre-defined subsurface targets. Since any efficient inversion routine requires a fast forward modeling engine, this study aimed the development of a fast forward modeling algorithm capable of being implemented in any inversion routine.To have efficient numerical forward modeling algorithm, we have adapted a leap-frog, staggered-grid approach introduced by Yee, which incorporates offsetting the electric and magnetic field components in both space and time in such a way that the finite difference approximations of the governing partial derivatives in each equation are centered on the same spatiotemporal location.
Obviously 2-D modeling is limited and cannot fully account for antenna behaviour and out-of-plane variations in material properties; however, many important features common to most GPR responses can be identified via employing a computationally cost effective 2D algorithm.
In the current study, the patterns of GPR responses that are well known to be hyperbola in shape are used as leading models in order to reduce the execution time. A GPR system collects the reflected pulses coming from different depths in the form of traces which when gathered along with a profile, they make a GPR section called radargram. In general, the simulated GPR traces of common reflected objects are time shifted like the Normal MoveOut (NMO) traces encountered in seismic reflection responses. This property suggests the application of Fourier transform to the GPR traces and the use of time shifting property of such transformation to interpolate traces between the adjusted traces in frequency domain. Therefore, the lateral resolution of GPR traces computed along with any profile is enhanced using a linear interpolation in the Fourier domain resulting in an increased speed of the forward modeling algorithm. Selecting the minimum lateral trace to trace interval with the appropriate sampling frequency of the signal, prevent any aliasing to occur. It is shown that such methodology can significantly decrease the computing time by more than 12.5 times.در تحقیق حاضر از بین روشهای عددی گوناگون دردسترس، از روش تفاضل متناهی دوبُعدی، به علت سادگی درک مفاهیم، انعطافپذیری، قابلیت شبیهسازی و مدلسازی محیطهای پیچیده و قابلقبول بودن پاسخهای آن در موارد بهکار رفته، برای شبیهسازی رفتار حاکم بر انتشار امواج الکترومغناطیسی (EM) در روش رادار نفوذی به زمین (GPR)، از راه حل معادلههای ماکسول و شرایط مرزی مناسب، استفاده شده است. روش تفاضل متناهی در حوزه زمان، علیرغم فایدههای بسیار، دارای نقاط ضعفی نیز میباشد ازجمله آنکه شبیهسازی مدلهای گوناگون بهویژه در محیطهای با گذردهی دیالکتریک زیاد، حتی با رایانههای نوین پُرسرعت نیز بسیار زمانبر است. لذا در تحقیق حاضر، این نقص نیز با اِعمال رویکردی جدید، بهبود یافت.
هدف از مدلسازی پیشرو پاسخ GPRیک ساختار زیرسطحی شناخته شده، بررسی تأثیر ویژگیهای الکتریکی و درنتیجه سنگشناختی یا آبشناختی آن ساختار روی پاسخهای GPRبرداشت شده از آن است بنابراین مدلسازیهای عددی پاسخ GPRمدلهای مصنوعی با هندسه مشخص، میتوانند در تعیین ارتباط بین ویژگیهای زیرسطحی و دادههای GPR، مفید واقع شوند. در این تحقیق برای مدلسازی دادههای GPRاز ثبت تپ (پالس)های بازتاب یافته به سطح (که فرستنده و گیرنده هر دو، روی سطح زمین قرار دارند) در مُد TMاستفاده شده است؛ بهطوریکه آنتنها عمود بر صفحه برداشت (x-z) توجیه شدهاند. به علت مدلسازی دوبُعدی، منبع و گیرنده بهصورت اجزای خطی هستند که عمود بر صفحه برداشت تا بینهایت گسترش مییابند؛ با این تفاسیر مدلسازی دوبُعدی صورت گرفته، قادر به مدلسازی بسیاری از ویژگیهای بازتابهای سامانههای GPRاست. نظر به اینکه فرآیند مدلسازی وارون دادههای GPR، نیازمند یک هسته مدلسازی پیشرو پُرسرعت است و ازآنجاییکه پاسخ GPRاغلبِ هدفهای مدفون در زیر زمین در نگاشتهای راداری حاصله که متشکل از تعداد زیادی تریس هستند، بهصورت یک رخداد شبههذلولی است، لذا در این تحقیق با استفاده از جابهجایی زمانی تریسها در الگوهای پاسخ GPRو تبدیل فوریه آنها در حوزه بسامدی و سپس درونیابی یک_بعدی در فاصله بین هر دو تریس متوالی، قدرت تفکیک جانبی دادهها را در امتداد پیمایش نیمرُخ، افزایش میدهد و درنتیجه پاسخهای با تفکیکپذیری بیشتر، حاصل میشود. همچنین با بکارگیری چنین روشی زمان محاسباتی شبیهسازیها، بهطور قابلتوجهی تا بیش از 5/12 برابر نسبت به حالتبهبود نیافته، کاهش مییابد. بهمنظور راستیآزمایی الگوریتم بهبودیافته، پاسخ GPRشیءهای مصنوعی نظیر استوانه افقی منفرد، منشور مربعی دوبُعدی و نیز استوانههای افقی جفت، مدلسازی شد.https://www.ijgeophysics.ir/article_33558_22a54a6f5b20403cc778dc5129b21e49.pdf